Source author record

David Vanderbilt

David Vanderbilt appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.mtrl-sci cond-mat.mes-hall cond-mat.str-el cond-mat.other physics.comp-ph hep-th mtrl-th physics.app-ph

Catalog footprint

What is connected

80works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Direct visualization of surface spin-flip transition in MnBi4Te7

We report direct visualization of spin-flip transition of the surface layer in antiferromagnet MnBi4Te7, a natural superlattice of alternating MnBi2Te4 and Bi2Te3 layers, using cryogenic magnetic force microscopy (MFM). The observation of magnetic contrast across domain walls and step edges confirms that the antiferromagnetic order persists to the surface layers. The magnetic field dependence of the MFM images reveals that the surface magnetic layer undergoes a first-order spin-flip transition at a magnetic field that is lower than the bulk transition, in excellent agreement with a revised Mills' model. Our analysis indicates an enhancement of the order parameter in the surface magnetic layer, implying robust ferromagnetism in the single-layer limit. The direct visualization of surface spin-flip transition not only opens up exploration of surface metamagnetic transitions in layered antiferromagnets, but also provides experimental support for realizing quantized transport in ultra-thin films of MnBi4Te7 and other natural superlattice topological magnets.

preprint2022arXiv

Electronic correlation in nearly free electron metals with beyond-DFT methods

For more than three decades, nearly free electron elemental metals have been a topic of debate because the computed bandwidths are significantly wider in the local density approximation to density-functional theory (DFT) than indicated by angle-resolved photoemission (ARPES) experiments. Here, we systematically investigate this using first-principles calculations for alkali and alkaline-earth metals using DFT and various beyond-DFT methods such as meta-GGA, G$_0$W$_0$, hybrid functionals (YS-PBE0, B3LYP), and LDA+eDMFT. We find that the static non-local exchange, as partly included in the hybrid functionals, significantly increase the bandwidths even compared to LDA, while the G$_0$W$_0$ bands are only slightly narrower than in LDA. The agreement with the ARPES is best when the local approximation to the self-energy is used in the LDA+eDMFT method. We infer that even moderately correlated systems with partially occupied s-orbitals, which were assumed to approximate the uniform electron gas, are very well described in terms of short-range dynamical correlations that are only local to an atom.

preprint2022arXiv

Engineering magnetic topological insulators in Eu$_5M_2X_6$ Zintls

Magnetic topological insulator provide a prominent material platform for quantum anomalous Hall physics and axion electrodynamics. However, the lack of material realizations with cleanly gapped surfaces hinders technological utilization of these exotic quantum phenomena. Here, using the Zintl concept and the properties of non-symmorphic space groups, we computationally engineer magnetic topological insulators. Specifically, we explore Eu$_5M_2X_6$ ($M$=metal, $X$=pnictide) Zintl compounds and find that Eu$_5$Ga$_2$Sb$_6$, Eu$_5$Tl$_2$Sb$_6$ and Eu$_5$In$_2$Bi$_6$ form stable structures with non-trivial $\mathbb{Z}_2$ indices. We also show that epitaxial and uniaxial strain can be used to control the $\mathbb{Z}_2$ index and the bulk energy gap. Finally, we discuss experimental progress towards the synthesis of the proposed candidates and provide insights that can be used in the search for robust magnetic topological insulators in Zintl compounds.

preprint2022arXiv

Vibrational fingerprints of ferroelectric hafnia

Hafnia (HfO2) is a promising material for emerging chip applications due to its high-k dielectric behaviour, suitability for negative capacitance heterostructures, scalable ferroelectricity, and silicon compatibility. The lattice dynamics along with phononic properties such as thermal conductivity, contraction, and heat capacity are under-explored, primarily due to the absence of high quality single crystals. Herein, we report the vibrational properties of a series of HfO2 crystals stabilized with yttrium (chemical formula HfO2:xY, where x = 20, 12, 11, 8, and 0%) and compare our findings with a symmetry analysis and lattice dynamics calculations. We untangle the effects of Y by testing our calculations against the measured Raman and infrared spectra of the cubic, antipolar orthorhombic, and monoclinic phases and then proceed to reveal the signature modes of polar orthorhombic hafnia. This work provides a spectroscopic fingerprint for several different phases of HfO2 and paves the way for an analysis of mode contributions to high-k dielectric and ferroelectric properties for chip technologies.

preprint2022arXiv

Vibrational properties of CuInP2S6 across the ferroelectric transition

In order to explore the properties of a two-sublattice ferroelectric, we measured the infrared and Raman scattering response of CuInP2S6 across the ferroelectric and glassy transitions and compared our findings to a symmetry analysis, calculations of phase stability, and lattice dynamics. In addition to uncovering displacive character and a large hysteresis region surrounding the ferroelectric transition temperature T_C, we identify the vibrational modes that stabilize the polar phase and confirm the presence of two ferroelectric variants with opposite polarizations. Below TC, a poorly understood relaxational or glassy transition at Tg is characterized by local structure changes in the form of subtle peak shifting and activation of low frequency out-of-plane Cu- and In-containing modes. The latter are due to changes in the Cu/In coordination environments and associated order-disorder processes. Moreover, Tg takes place in two steps with another large hysteresis region and significant underlying scattering. Combined with imaging of the room temperature phase separation, this effort lays the groundwork for studying CuInP2S6 under external stimuli and in the ultra-thin limit.

preprint2021arXiv

Humble planar defects in SiGe nanopillars

We report a new \{001\} planar defect found in SiGe nanopillars. The defect structure, determined by atomic resolution electron microscopy, matches the Humble defect model proposed for diamond. We also investigated several possible variants of the Humble structure using first principles calculations and found that the one lowest in energy was also in best agreement with the STEM images. The pillar composition has been analyzed with electron energy loss spectroscopy, which hints at how the defect is formed. Our results show that the structure and formation process of defects in nanostructured group IV semiconductors can be different from their bulk counterparts.

preprint2021arXiv

Lattice dynamics and structural transition of the hyperhoneycomb iridate $β$-Li$_2$IrO$_3$ investigated by high-pressure Raman scattering

We report a polarized Raman scattering study of the lattice dynamics of $β$-Li$_2$IrO$_3$ under hydrostatic pressures up to 7.62 GPa. At ambient pressure, $β$-Li$_2$IrO$_3$ exhibits the hyperhoneycomb crystal structure and a magnetically ordered state of spin-orbit entangled Jeff = 1/2 moments that is strongly influenced by bond-directional (Kitaev) exchange interactions. At a critical pressure of ~ 4.1 GPa, the phonon spectrum changes abruptly consistent with the reported structural transition into a monoclinic, dimerized phase. A comparison to the phonon spectra obtained from density functional calculations shows reasonable overall agreement. The calculations also indicate that the high-pressure phase is a nonmagnetic insulator driven by the formation of Ir-Ir dimer bonds. Our results thus indicate a strong sensitivity of the electronic properties of $β$-Li$_2$IrO$_3$ to the pressure-induced structural transition.

preprint2021arXiv

Orbital-selective Mott phase and non-Fermi liquid in FePS$_3$

The layered metal phosphorous trisulfide FePS$_3$ is reported to be a Mott insulator at ambient conditions and to undergo structural and insulator-metal phase transitions under pressure. However, the character of the resulting metallic states has not been understood clearly so far. Here, we theoretically study the phase transitions of FePS$_3$ using first-principles methods based on density functional theory and embedded dynamical mean field theory. We find that the Mott transition in FePS$_3$ can be orbital-selective, with $t_{2g}$ states undergoing a correlation-induced insulator-to-metal transition while $e_g$ states remain gapped. We show that this orbital-selective Mott phase, which occurs only when non-hydrostatic pressure is used, is a bad metal (or non-Fermi liquid) with large fluctuating moments due to Hund's coupling. Further application of pressure increases the crystal-field splitting and converts the system to a conventional Fermi liquid with low-spin configurations dominant. Our results show that FePS$_3$ is a novel example of a system that realizes an orbital-selective Mott phase, allowing tuning between correlated and uncorrelated metallic properties in an accessible pressure range ($\leq$ 18 GPa).

preprint2021arXiv

Phonon spectrum of Pr$_2$Zr$_2$O$_7$ and Pr$_2$Ir$_2$O$_7$ as an evidence of coupling of the lattice with electronic and magnetic degrees of freedom

Magnetic materials with pyrochlore crystal structure form exotic magnetic states due to the high lattice frustration. In this work we follow the effects of coupling of the lattice and electronic and magnetic degrees of freedom in two Praseodymium-based pyrochlores Pr$_2$Zr$_2$O$_7$ and Pr$_2$Ir$_2$O$_7$. In both materials the presence of magnetic interactions does not lead to magnetically ordered low temperature states, however their electronic properties are different. A comparison of Raman phonon spectra of Pr$_2$Zr$_2$O$_7$ and Pr$_2$Ir$_2$O$_7$ allows us to identify magneto-elastic coupling in Pr$_2$Zr$_2$O$_7$ that elucidates its magnetic properties at intermediate temperatures, and allows us to characterize phonon-electron coupling in the semimetallic Pr$_2$Ir$_2$O$_7$. We also show that the effects of random disorder on the Raman phonon spectra is small.

preprint2021arXiv

Quadrupole moments, edge polarizations, and corner charges in the Wannier representation

The modern theory of polarization allows for the determination of the macroscopic end charge of a truncated one-dimensional insulator, modulo the charge quantum $e$, from a knowledge of bulk properties alone. A more subtle problem is the determination of the corner charge of a two-dimensional insulator, modulo $e$, from a knowledge of bulk and edge properties alone. While previous works have tended to focus on the quantization of corner charge in the presence of symmetries, here we focus on the case that the only bulk symmetry is inversion, so that the corner charge can take arbitrary values. We develop a Wannier-based formalism that allows the corner charge to be predicted, modulo $e$, only from calculations on ribbon geometries of two different orientations. We elucidate the dependence of the interior quadrupole and edge dipole contributions upon the gauge used to construct the Wannier functions, finding that while these are individually gauge-dependent, their sum is gauge-independent. From this we conclude that the edge polarization is not by itself a physical observable, and that any Wannier-based method for computing the corner charge requires the use of a common gauge throughout the calculation. We satisfy this constraint using two Wannier construction procedures, one based on projection and another based on a gauge-consistent nested Wannier construction. We validate our theory by demonstrating the correct prediction of corner charge for several tight-binding models. We comment on the relations between our approach and previous ones that have appeared in the literature.

preprint2020arXiv

Axion coupling in the hybrid Wannier representation

Many magnetic point-group symmetries induce a topological classification on crystalline insulators, dividing them into those that have a nonzero quantized Chern-Simons magnetoelectric coupling ("axion-odd" or "topological"), and those that do not ("axion-even" or "trivial"). For time-reversal or inversion symmetries, the resulting topological state is usually denoted as a "strong topological insulator" or an "axion insulator" respectively, but many other symmetries can also protect this "axion $Z_2$" index. Topological states are often insightfully characterized by choosing one crystallographic direction of interest, and inspecting the hybrid Wannier (or equivalently, the non-Abelian Wilson-loop) band structure, considered as a function of the two-dimensional Brillouin zone in the orthogonal directions. Here, we systematically classify the axion-quantizing symmetries, and explore the implications of such symmetries on the Wannier band structure. Conversely, we clarify the conditions under which the axion$Z_2$ index can be deduced from the Wannier band structure. In particular, we identify cases in which a counting of Dirac touchings between Wannier bands, or a calculation of the Chern number of certain Wannier bands, or both, allows for a direct determination of the axion $Z_2$ index. We also discuss when such symmetries impose a "flow" on the Wannier bands, such that they are always glued to higher and lower bands by degeneracies somewhere in the projected Brillouin zone, and the related question of when the corresponding surfaces can remain gapped, thus exhibiting a half-quantized surface anomalous Hall conductivity. Our formal arguments are confirmed and illustrated in the context of tight-binding models for several paradigmatic axion-odd symmetries including time reversal, inversion, simple mirror, and glide mirror symmetries.

preprint2020arXiv

Berry Flux Diagonalization: Application to Electric Polarization

The switching polarization of a ferroelectric is a characterization of the current that flows due to changes in polarization when the system is switched between two states. Computation of this change in polarization in crystal systems has been enabled by the modern theory of polarization, where it is expressed in terms of a change in Berry phase as the material switches. It is straightforward to compute this change of phase, but only modulo $2π$, requiring a branch choice from among a lattice of values separated by $2π$. The measured switching polarization depends on the actual path along which the material switches, which in general involves nucleation and growth of domains and is therefore quite complex. In this work, we present a physically motivated approach for predicting the experimentally measured switching polarization that involves separating the change in phase between two states into as many gauge-invariant smaller phase changes as possible. As long as the magnitudes of these smaller phase changes remain smaller than $π$, their sum forms a phase change which corresponds to the change one would find along any path involving minimal evolution of the atomic and electronic structure. We show that for typical ferroelectrics, including those that would have otherwise required a densely sampled path, this technique allows the switching polarization to be computed without any need for intermediate sampling between oppositely polarized states.

preprint2020arXiv

Engineering Weyl phases and nonlinear Hall effects in T$_d$-MoTe$_2$

MoTe$_2$ has recently attracted much attention due to the observation of pressure-induced superconductivity, exotic topological phase transitions, and nonlinear quantum effects. However, there has been debate on the intriguing structural phase transitions among various observed phases of MoTe$_2$, and their connection to the underlying topological electronic properties. In this work, by means of density-functional theory (DFT+U) calculations, we investigate the structural phase transition between the polar T$_d$ and nonpolar 1T$'$ phases of MoTe$_2$ in reference to a hypothetical high-symmetry T$_0$ phase that exhibits higher-order topological features. In the T$_d$ phase we obtain a total of 12 Weyl points, which can be created/annihilated, dynamically manipulated, and switched by tuning a polar phonon mode. We also report the existence of a tunable nonlinear Hall effect in T$_d$-MoTe$_2$, and propose the use of this effect as a probe for the detection of polarity orientation in polar (semi)metals. By studying the role of dimensionality, we identify a configuration in which a nonlinear surface response current emerges. The potential technological applications of the tunable Weyl phase and the nonlinear Hall effect are discussed.

preprint2020arXiv

Gapless hinge states from adiabatic pumping of axion coupling

We demonstrate that chiral hinge modes naturally emerge in insulating crystals undergoing a slow cyclic evolution that changes the Chern-Simons axion angle $θ$ by $2π$. This happens when the surface (not just the bulk) returns to its initial state at the end of the cycle, in which case it must pass through a metallic state to dispose of the excess quantum of surface anomalous Hall conductivity pumped from the bulk. If two adjacent surfaces become metallic at different points along the cycle, there is an interval in which they are in topologically distinct insulating states, with chiral modes propagating along the connecting hinge. We illustrate these ideas for a tight-binding model consisting of coupled layers of the Haldane model with alternating parameters. The surface topology is determined in a slab geometry using two different markers, surface anomalous Hall conductivity and surface-localized charge pumping (flow of surface-localized Wannier bands), and we find that both correctly predict the appearance of gapless hinge modes in a rod geometry. When viewing the axion pump as a four-dimensional (4D) crystal with one synthetic dimension, the hinge modes trace Fermi arcs in the Brillouin zone of the 2D hinge connecting a pair of 3D surfaces of the 4D crystal.

preprint2020arXiv

Magnetic phase transitions and spin density distribution in the molecular multiferroic GaV$_4$S$_8$ system

We have carried out neutron diffraction and small angle neutron scattering measurements on a high quality single crystal of the cubic lacunar spinel multiferroic, GaV$_4$S$_8$, as a function of magnetic field and temperature to determine the magnetic properties for the single electron that is located on the tetrahedrally coordinated V$_4$ molecular unit. Our results are in good agreement with the structural transition at 44 K from cubic to rhombohedral symmetry where the system becomes a robust ferroelectric, while long range magnetic order develops below 13 K in the form of an incommensurate cycloidal magnetic structure, which can transform into a Néel-type skyrmion phase in a modest applied magnetic field. Below 5.9(3) K, the crystal enters a ferromagnetic phase, and we find the magnetic order parameter indicates a long range ordered ground state with an ordered moment of 0.23(1) $μ_\mathrm{B}$ per V ion. Both polarized and unpolarized neutron data in the ferroelectric-paramagnetic phase have been measured to determine the magnetic form factor. The data are consistent with a model of the single spin being uniformly distributed across the V$_4$ molecular unit, rather than residing on the single apical V ion, in substantial agreement with the results of first-principles theory. In the magnetically ordered state, polarized neutron measurements are important since both the cycloidal and ferromagnetic order parameters are clearly coupled to the ferroelectricity, causing the structural peaks to be temperature and field dependent. For the ferromagnetic ground state, the spins are locked along the $[1,1,1]$ direction by a surprisingly large anisotropy.

preprint2020arXiv

Molecular Mott state in the deficient spinel GaV$_4$S$_8$

In this study, we investigated theoretically the Mott-insulating phase of a deficient spinel chalcogenide GaV$_4$S$_8$, which is known to form a tetrahedral V$_4$S$_4$ cluster unit that results in molecular orbitals (MOs) with a narrow bandwidth in the noninteracting limit. We used a cluster extension of charge self-consistent embedded dynamical mean-field theory to study the impact of strong intra-cluster correlations on the spectral properties as well as the structural degrees of freedom of the system. We found that the strong tetrahedral clustering renders the atomic Mott picture ineffective, and that the resulting MO picture is essential to describe the Mott phase. It was also found that, while the spectral properties can be qualitatively described by the truncation of the Hilbert space down to the lowest-energy MO, a proper description of the structural degrees of freedom requires the inclusion of multi-MO correlations that span a larger energy window. Specifically, we found that the lowest-energy MO description overemphasizes the clustering tendency, while the inclusion of the Hund's coupling between the lower- and higher-energy MOs corrects this tendency, bringing the theoretically predicted crystal structure into good agreement with the experiment.

preprint2020arXiv

Robust A-type order and spin-flop transition on the surface of the antiferromagnetic topological insulator MnBi$_2$Te$_4$

Here we present microscopic evidence of the persistence of uniaxial A-type antiferromagnetic order to the surface layers of MnBi$_2$Te$_4$ single crystals using magnetic force microscopy. Our results reveal termination-dependent magnetic contrast across both surface step edges and domain walls, which can be screened by thin layers of soft magnetism. The robust surface A-type order is further corroborated by the observation of termination-dependent surface spin-flop transitions, which have been theoretically proposed decades ago. Our results not only provide key ingredients for understanding the electronic properties of the antiferromagnetic topological insulator MnBi$_2$Te$_4$, but also open a new paradigm for exploring intrinsic surface metamagnetic transitions in natural antiferromagnets.

preprint2019arXiv

Systematic beyond-DFT study of binary transition metal oxides

Various methods going beyond density-functional theory (DFT), such as DFT+U, hybrid functionals, meta-GGAs, GW, and DFT-embedded dynamical mean field theory (eDMFT), have been developed to describe the electronic structure of correlated materials, but it is unclear how accurate these methods can be expected to be when applied to a given strongly correlated solid. It is thus of pressing interest to compare their accuracy as they apply to different categories of materials. Here, we introduce a novel paradigm in which a chosen set of beyond-DFT methods is systematically and uniformly tested on a chosen class of materials. For a first application, we choose the target materials to be the binary transition-metal oxides FeO, CoO, MnO, and NiO in their antiferromagnetic phase and present a head-to-head comparison of spectral properties as computed using the various methods. We also compare with available experimental angle-resolved photoemission spectroscopy (ARPES), inverse-photoemission spectroscopy, and with optical absorption. We find both B3LYP and eDMFT can reproduce the experiments quite well, with eDMFT doing best in particular when comparing with the ARPES data.

preprint2019arXiv

Wannier90 as a community code: new features and applications

Wannier90 is an open-source computer program for calculating maximally-localised Wannier functions (MLWFs) from a set of Bloch states. It is interfaced to many widely used electronic-structure codes thanks to its independence from the basis sets representing these Bloch states. In the past few years the development of Wannier90 has transitioned to a community-driven model; this has resulted in a number of new developments that have been recently released in Wannier90 v3.0. In this article we describe these new functionalities, that include the implementation of new features for wannierisation and disentanglement (symmetry-adapted Wannier functions, selectively-localised Wannier functions, selected columns of the density matrix) and the ability to calculate new properties (shift currents and Berry-curvature dipole, and a new interface to many-body perturbation theory); performance improvements, including parallelisation of the core code; enhancements in functionality (support for spinor-valued Wannier functions, more accurate methods to interpolate quantities in the Brillouin zone); improved usability (improved plotting routines, integration with high-throughput automation frameworks), as well as the implementation of modern software engineering practices (unit testing, continuous integration, and automatic source-code documentation). These new features, capabilities, and code development model aim to further sustain and expand the community uptake and range of applicability, that nowadays spans complex and accurate dielectric, electronic, magnetic, optical, topological and transport properties of materials.

preprint2016arXiv

Ferroelectricity in corundum derivatives

The search for new ferroelectric (FE) materials holds promise for broadening our understanding of FE mechanisms and extending the range of application of FE materials. Here we investigate a class of ABO3 and A2BB'O6 materials that can be derived from the X2O3 corundum structure by mixing two or three ordered cations on the X site. Most such corundum derivatives have a polar structure, but it is unclear whether the polarization is reversible, which is a requirement for a FE material. In this paper, we propose a method to study the FE reversal path of materials in the corundum derivative family. We first categorize the corundum derivatives into four classes and show that only two of these allow for the possibility of FE reversal. We then calculate the energy profile and energy barrier of the FE reversal path using first-principles density functional methods with a structural constraint. Furthermore, we identify several empirical measures that can provide a rule of thumb for estimating the energy barriers. Finally, the conditions under which the magnetic ordering is compatible with ferroelectricity are determined. These results lead us to predict several potentially new FE materials.

preprint2016arXiv

Flux states and topological phases from spontaneous time-reversal symmetry breaking in CrSi(Ge)Te$_3$-based systems

We study adatom-covered single layers of CrSiTe$_3$ and CrGeTe$_3$ using first-principles calculations based on hybrid functionals. We find that the insulating ground state of a monolayer of La (Lu) deposited on single-layer CrSiTe$_3$ (CrGeTe$_3$) carries spontaneously generated current loops around the Cr sites, These "flux states" induce antiferromagnetically ordered orbital moments on the Cr sites and are also associated with nontrivial topological properties. The calculated Chern numbers for these systems are predicted to be $\pm 1$ even in the absence of spin-orbit coupling, with sizable gaps on the order of 100 meV. The flux states and the associated topological phases result from spontaneous time-reversal symmetry breaking due to the presence of nonlocal Coulomb interactions.

preprint2016arXiv

Intertwined Rashba, Dirac and Weyl Fermions in Hexagonal Hyperferroelectrics

By means of density functional theory based calculations, we study the role of spin-orbit coupling in the new family of ABC hyperferroelectrics [Phys. Rev. Lett. 112, 127601 (2014)]. We unveil an extremely rich physics strongly linked to ferroelectric properties, ranging from the electric control of bulk Rashba effect to the existence of a three dimensional topological insulator phase, with concomitant topological surface states even in the ultrathin film limit. Moreover, we predict that the topological transition, as induced by alloying, is followed by a Weyl semi-metal phase of finite concentration extension, which is robust against disorder, putting forward hyperferroelectrics as promising candidates for spin-orbitronic applications.

preprint2016arXiv

Temperature effects in the band structure of topological insulators

We study the effects of temperature on the band structure of the Bi$_2$Se$_3$ family of topological insulators using first-principles methods. Increasing temperature drives these materials towards the normal state, with similar contributions from thermal expansion and from electron-phonon coupling. The band gap changes with temperature reach $0.3$ eV at $600$ K, of similar size to the changes caused by electron correlation. Our results suggest that temperature-induced topological phase transitions should be observable near critical points of other external parameters.

preprint2016arXiv

Topological nodal-line semimetals in alkaline-earth stannides, germanides and silicides

Based on first-principles calculations and an effective Hamiltonian analysis, we systematically investigate the electronic and topological properties of alkaline-earth compounds $AX_2$ ($A$=Ca, Sr, Ba; $X$=Si, Ge, Sn). Taking BaSn$_2$ as an example, we find that when spin-orbit coupling is ignored, these materials are three-dimensional topological nodal-line semimetals characterized by a snake-like nodal loop in three-dimensional momentum space. Drumhead-like surface states emerge either inside or outside the loop circle on the (001) surface depending on surface termination, while complicated double-drumhead-like surface states appear on the (010) surface. When spin-orbit coupling is included, the nodal line is gapped and the system becomes a topological insulator with Z$_2$ topological invariants (1;001). Since spin-orbit coupling effects are weak in light elements, the nodal-line semimetal phase is expected to be achievable in some alkaline-earth germanides and silicides.

preprint2016arXiv

Tunable inverse topological heterostructure utilizing $(Bi_{1-x}In_{x})_{2}Se_{3}$

In typical topological insulator (TI) systems the TI is bordered by a non-TI insulator, and the surrounding conventional insulators, including vacuum, are not generally treated as part of the TI system. Here, we implement the first material system where the roles are reversed, and the TSS form around the non-TI (instead of the TI) layers. This is realized by growing a layer of the tunable non-TI $(Bi_{1-x}In_{x})_{2}Se_{3}$ in between two layers of the TI $Bi_2Se_3$ using the atomically-precise molecular beam epitaxy technique. On this tunable inverse topological platform, we systematically vary the thickness and the composition of the $(Bi_{1-x}In_{x})_{2}Se_{3}$ layer and show that this tunes the coupling between the TI layers from strongly-coupled metallic to weakly-coupled, and finally to a fully-decoupled insulating regime. This system can be used to probe the fundamental nature of coupling in TI materials and provides a tunable insulating layer for TI devices.

preprint2015arXiv

Chiral degeneracies and Fermi-surface Chern numbers in bcc Fe

The degeneracies in the spinor bandstructure of bcc Fe are studied from first principles. We find numerous isolated band touchings carrying chiral charges of magnitude one (Weyl points) or two (double-Weyl nodes), as well as nonchiral degeneracy loops (nodal rings). Some degeneracies are located on symmetry lines or planes in the Brillouin zone and others at generic low-symmetry points, realizing all possible scenarios consistent with the magnetic point group. We clarify the general theory relating the chiral band touchings to the Chern numbers of the Fermi sheets enclosing them, and use this approach to determine the Chern numbers on the Fermi surface of bcc Fe. Although most Fermi sheets enclose Weyl nodes, in almost all cases the net enclosed charge vanishes for symmetry reasons, resulting in a vanishing Chern number. The exceptions are two inversion-symmetric electron pockets along the symmetry line $Δ$ parallel to the magnetization. Each of them surrounds a single Weyl point, leading to Chern numbers of $\pm 1$. These small topological pockets contribute a sizable amount to the nonquantized part of the intrinsic anomalous Hall conductivity, proportional to their reciprocal-space separation. Variation of the Fermi level (or other system parameters) may lead to a touching event between Fermi sheets, accompanied by a transfer of Chern number between them.

preprint2015arXiv

Emergence of a Chern-insulating state from a semi-Dirac dispersion

A Chern insulator (quantum anomalous Hall insulator) phase is demonstrated to exist in a typical semi-Dirac system, the TiO2/VO2 heterostructure. By combining first-principles calculations with Wannier-based tight-binding model, we calculate the Berry curvature distribution, finding a Chern number of -2 for the valence bands, and demonstrate the existence of gapless chiral edge states, ensuring quantization of the Hall conductivity to 2e^2/h. A new semi-Dirac model, where each semi-Dirac cone is formed by merging three conventional Dirac points, is proposed to reveal how the nontrivial topology with finite Chern number is compatible with a semi-Dirac electronic spectrum.

preprint2015arXiv

First-principles theory of flexoelectricity

In this Chapter we provide an overview of the current first-principles perspective on flexoelectric effects in crystalline solids. We base our theoretical formalism on the long-wave expansion of the electrical response of a crystal to an acoustic phonon perturbation. In particular, we recover the known expression for the piezoelectric tensor from the response at first order in wavevector ${\bf q}$, and then obtain the flexoelectric tensor by extending the formalism to second order in $\bf q$. We put special emphasis on the issue of surface effects, which we first analyze heuristically, and then treat more carefully by presenting a general theory of the microscopic response to an arbitrary inhomogeneous strain. We demonstrate our approach by presenting a full calculation of the flexoelectric response of a SrTiO$_3$ film, where we point out an unusually strong dependence of the bending-induced open-circuit voltage on the choice of surface termination. Finally, we briefly discuss some remaining open issues concerning the methodology and some promising areas for future research.

preprint2015arXiv

Gauge-discontinuity contributions to Chern-Simons orbital magnetoelectric coupling

We propose a new method for calculating the Chern-Simons orbital magnetoelectric coupling, conventionally parametrized in terms of a phase angle $θ$. According to previous theories, $θ$ can be expressed as a 3D Brillouin-zone integral of the Chern-Simons 3-form defined in terms of the occupied Bloch functions. Such an expression is valid only if a smooth and periodic gauge has been chosen in the entire Brillouin zone, and even then, convergence with respect to the $\mathbf{k}$-space mesh density can be difficult to obtain. In order to solve this problem, we propose to relax the periodicity condition in one direction (say, the $k_z$ direction) so that a gauge discontinuity is introduced on a 2D $\mathbf{k}$ plane normal to $k_z$. The total $θ$ response then has contributions from both the integral of the Chern-Simons 3-form over the 3D bulk BZ and the gauge discontinuity expressed as a 2D integral over the $\mathbf{k}$ plane. Sometimes the boundary plane may be further divided into subregions by 1D "vortex loops" which make a third kind of contribution to the total $θ$, expressed as a combination of Berry phases around the vortex loops. The total $θ$ thus consists of three terms which can be expressed as integrals over 3D, 2D and 1D manifolds. When time-reversal symmetry is present and the gauge in the bulk BZ is chosen to respect this symmetry, both the 3D and 2D integrals vanish; the entire contribution then comes from the vortex-loop integral, which is either 0 or $π$ corresponding to the $\mathbb{Z}_2$ classification of 3D time-reversal invariant insulators. We demonstrate our method by applying it to the Fu-Kane-Mele model with an applied staggered Zeeman field.

preprint2015arXiv

Magnetic charges and magnetoelectricity in hexagonal rare-earth manganites and ferrites

Magnetoelectric (ME) materials are of fundamental interest and show broad potential for technological applications. Commonly the dominant contribution to the ME response is the lattice-mediated one, which is proportional to both the Born electric charge $Z^{\rm e}$ and its analogue, the dynamical magnetic charge $Z^{\rm m}$. Our previous study has shown that exchange striction acting on noncollinear spins induces much larger magnetic charges than those that depend on spin-orbit coupling. The hexagonal manganites $R$MnO$_3$ and ferrites $R$FeO$_3$ ($R$ = Sc, Y, In, Ho-Lu) exhibit strong couplings between electric, magnetic and structural degrees of freedom, with the transition-metal ions in the basal plane antiferromagnetically coupled through super-exchange so as to form a 120$^\circ$ noncollinear spin arrangement. Here we present a theoretical study of the magnetic charges, and of the spin-lattice and spin-electronic ME constants, in these hexagonal manganites and ferrites, clarifying the conditions under which exchange striction leads to an enhanced $Z^{\rm m}$ values and anomalously large in-plane spin-lattice ME effects.

preprint2015arXiv

Surface polarization and edge charges

The term "surface polarization" is introduced to describe the in-plane polarization existing at the surface of an insulating crystal when the in-plane surface inversion symmetry is broken. Here, the surface polarization is formulated in terms of a Berry phase, with the hybrid Wannier representation providing a natural basis for study of this effect. Tight binding models are used to demonstrate how the surface polarization reveals itself via the accumulation of charges at the corners/edges for a two dimensional rectangular lattice and for GaAs.

preprint2014arXiv

Adiabatic Pumping of Chern-Simons Axion Coupling

We study the adiabatic pumping of the Chern-Simons axion (CSA) coupling along a parametric loop characterized by a non-zero second Chern number $C^{(2)}$ from the viewpoint of the hybrid Wannier representation, in which the Wannier charge centers (WCCs) are visualized as sheets defined over a projected 2D Brillouin zone. We derive a new formula for the CSA coupling, expressing it as an integral involving Berry curvatures and potentials defined on the WCC sheets. We show that a loop characterized by a non-zero $C^{(2)}$ requires a series of sheet-touching events at which $2π$ quanta of Berry curvature are passed from sheet to sheet, in such a way that $e^2/h$ units of CSA coupling are pumped by a lattice vector by the end of the cycle. We illustrate these behaviors via explicit calculations on a model tight-binding Hamiltonian and discuss their implications.

preprint2014arXiv

Chern insulator at a magnetic rocksalt interface

Considerable efforts have recently been devoted to the experimental realization of a two-dimensional Chern insulator, i.e., a system displaying a quantum anomalous Hall effect. However, existing approaches such as those based on magnetic doping of topological-insulator thin films have resulted in small band gaps, restricting the effect to low temperatures. We use first-principles calculations to demonstrate that an interface between thin films of the topologically trivial ferromagnetic insulators EuO and GdN can result in a band inversion and a non-zero Chern number. Both materials are stoichiometric and the interface is non-polar and lattice-matched, which should allow this interface to be achievable experimentally. We show that the band structure can be tuned by layer thickness or epitaxial strain, and can result in Chern insulators with gaps of over 0.1 eV.

preprint2014arXiv

Dynamical magnetic charges and linear magnetoelectricity

Magnetoelectric (ME) materials are of fundamental interest and have been investigated for their broad potential for technological applications. The search for, and eventually the theoretical design of, materials with large ME couplings present challenging issues. First-principles methods have only recently been developed to calculate the full ME response tensor $α$ including both electronic and ionic (i.e., lattice-mediated) contributions. The latter is proportional to both the Born dynamical electric charge $Z^{\rm e}$ and its analogue, the dynamical magnetic charge $Z^{\rm m}$. Here we present a theoretical study of the magnetic charge $Z^{\rm m}$ and the mechanisms that could enhance it. Using first-principles density-functional methods, we calculate the atomic $Z^{\rm m}$ tensors in $\rm{Cr_2O_3}$, a prototypical magnetoelectric, and in KITPite, a fictitious material that has previously been reported to show a strong ME response arising from exchange striction effects. Our results confirm that in $\rm{Cr_2O_3}$, the $Z^{\rm m}$ values and resulting ME responses arise only from spin-orbit coupling (SOC) and are therefore rather weak. In KITPite, by contrast, the exchange striction acting on the non-collinear spin structure induces much $Z^{\rm m}$ values that persist even when SOC is completely absent.

preprint2014arXiv

Quantum anomalous Hall phase in (001) double-perovskite monolayers via intersite spin-orbit coupling

Using tight-binding models and first-principles calculations, we demonstrate the possibility to achieve a quantum anomalous Hall (QAH) phase on a two-dimensional square lattice, which can be realized in monolayers of double perovskites. We show that effective intersite spin-orbit coupling between eg orbitals can be induced perturbatively, giving rise to a QAH state. Moreover, the effective spin-orbit coupling can be enhanced by octahedral rotations. Based on first-principles calculations, we propose that this type of QAH state could be realized in La2MnIrO6 monolayers, with the size of the gap as large as 26 meV in the ideal case. We observe that the electronic structure is sensitive to structural distortions, and that an enhanced Hubbard U tends to stabilize the nontrivial gap.

preprint2014arXiv

Spin-orbit spillage as a measure of band inversion in insulators

We propose a straightforward and effective approach for quantifying the band inversion induced by spin-orbit coupling in band insulators. In this approach we define a quantity as a function of wavevector in the Brillouin zone measuring the spillage between the occupied states of a system with and without SOC. Plots of the spillage throughout the BZ provide a ready indication of the number and location of band inversions driven by SOC. To illustrate the method, we apply this approach to the 2D Kane-Mele model, a 2D Bi bilayer with applied Zeeman field, and to first-principles calculations of some 3D materials including both trivial and \Z2\ topological insulators. We argue that the distribution of spillage in the BZ is closely related to the topological indices in these systems. Our approach provides a fresh perspective for understanding topological character in band theory, and should be helpful in searching for new materials with non-trivial band topology.

preprint2014arXiv

Weyl semimetals from noncentrosymmetric topological insulators

We study the problem of phase transitions from 3D topological to normal insulators without inversion symmetry. In contrast with the conclusions of some previous work, we show that a Weyl semimetal always exists as an intermediate phase regardless of any constriant from lattice symmetries, although the interval of the critical region is sensitive to the choice of path in the parameter space and can be very narrow. We demonstrate this behavior by carrying out first-principles calculations on the noncentrosymmetric topological insulators LaBiTe$_3$ and LuBiTe$_3$ and the trivial insulator BiTeI. We find that a robust Weyl-semimetal phase exists in the solid solutions LaBi$_{1-x}$Sb$_x$Te$_3$ and LuBi$_{1-x}$Sb$_x$Te$_3$ for $x\!\approx\!38.5-41.9$\% and $x\!\approx\!40.5-45.1$\% respectively. A low-energy effective model is also constructed to describe the critical behavior in these two materials. In BiTeI, a Weyl semimetal also appears with applied pressure, but only within a very small pressure range, which may explain why it has not been experimentally observed.

preprint2013arXiv

Canonical magnetic insulators with isotropic magnetoelectric coupling

We have performed an exhaustive representation-theory-based search for the simplest structures allowing isotropic magnetoelectric coupling. We find 30 such structures, all sharing a common pattern of atomic displacements in the direction of atomic magnetic moments. We focus on one of these 30 canonical structures and find that it is generically realized in a class of fractionally substituted pyrochlore compounds with an all-in-all-out magnetic order. Furthermore, we find that these substituted pyrochlore compounds have a substantial Chern-Simons orbital magnetoelectric component (θ=0.1--0.2). While this component is also formally present in strong Z_2 topological insulators (θ=π), its effects are observable there only if time-reversal symmetry is broken at the surface.

preprint2013arXiv

Comment on "Weyl fermions and the anomalous Hall effect in metallic ferromagnets"

We point out that, contrary to an assertion by Chen, Bergman and Burkov [Phys. Rev. B 88, 125110 (2013)], the non-quantized part of the intrinsic anomalous Hall conductivity can indeed be expressed as a Fermi-surface property even when Weyl points are present in the bandstructure.

preprint2013arXiv

Effective J=1/2 insulating state in Ruddlesden-Popper iridates: An LDA+DMFT study

Using ab-initio methods for correlated electrons in solids, we investigate the metal-insulator transition across the Ruddlesden-Popper (RP) series of iridates and explore the robustness of the J_eff=1/2 state against band effects due to itineracy, tetragonal distortion, octahedral rotation and Coulomb interaction. We predict the effects of epitaxial strain on the optical conductivity, magnetic moments, and J_eff=1/2 ground-state wave functions in the RP series. To describe the solution of the many-body problem in an intuitive picture, we introduce a concept of energy-dependent atomic states, which strongly resemble the atomic J_eff=1/2 states but with coefficients that are energy/time-dependent. We demonstrate that the deviation from the ideal J_eff=1/2 state is negligible at short time scales for both single- and double-layer iridates, while it becomes quite significant for Sr3Ir2O7 at long times and low energy. Interestingly, Sr2IrO4 is positioned very close to the SU(2) limit, with only 3% deviation from the ideal J_eff=1/2 situation.

preprint2013arXiv

First-principles theory and calculation of flexoelectricity

We develop a general and unified first-principles theory of piezoelectric and flexoelectric tensor, formulated in such a way that the tensor elements can be computed directly in the context of density-functional calculations, including electronic and lattice contributions. We introduce a practical supercell-based methods for calculating the flexoelectric coefficients from first principles, and demonstrate them by computing the coefficients for a variety of cubic insulating materials, including C, Si, MgO, NaCl, CsCl, BaZrO3, BaTiO3, PbTiO3 and SrTiO3.

preprint2013arXiv

Hyperferroelectrics: proper ferroelectrics with persistent polarization

All known proper ferroelectrics are unable to polarize normal to a surface or interface if the resulting depolarization field is unscreened, but there is no fundamental principle that enforces this behavior. In this work, we introduce hyperferroelectrics, a new class of proper ferroelectrics which polarize even when the depolarization field is unscreened, this condition being equivalent to instability of a longitudinal optic mode in addition to the transverse-optic-mode instability characteristic of proper ferroelectrics. We use first principles calculations to show that several recently discovered hexagonal ferroelectric semiconductors have this property, and we examine its consequences both in the bulk and in a superlattice geometry.

preprint2013arXiv

Topological phase transitions in (Bi$_{1-x}$In$_{x})_2$Se$_3$ and (Bi$_{1-x}$Sb$_{x})_2$Se$_3$

We study the phase transition from a topological to a normal insulator with concentration $x$ in (Bi$_{1-x}$In$_{x})_2$Se$_3$ and (Bi$_{1-x}$Sb$_{x})_2$Se$_3$ in the Bi$_2$Se$_3$ crystal structure. We carry out first-principles calculations on small supercells, using this information to build Wannierized effective Hamiltonians for a more realistic treatment of disorder. Despite the fact that the spin-orbit coupling (SOC) strength is similar in In and Sb, we find that the critical concentration $x_{\rm c}$ is much smaller in (Bi$_{1-x}$In$_{x})_2$Se$_3$ than in (Bi$_{1-x}$Sb$_{x})_2$Se$_3$. For example, the direct supercell calculations suggest that $x_{\rm c}$ is below 12.5% and above 87.5$%$ for the two alloys respectively. More accurate results are obtained from realistic disordered calculations, where the topological properties of the disordered systems are understood from a statistical point of view. Based on these calculations, $x_c$ is around 17% for (Bi$_{1-x}$In$_{x})_2$Se$_3$, but as high as 78%-83% for (Bi$_{1-x}$Sb$_{x})_2$Se$_3$. In (Bi$_{1-x}$Sb$_{x})_2$Se$_3$, we find that the phase transition is dominated by the decrease of SOC, with a crossover or "critical plateau" observed from around 78$%$ to 83$%$. On the other hand, for (Bi$_{1-x}$In$_{x})_2$Se$_3$, the In 5$s$ orbitals suppress the topological band inversion at low impurity concentration, therefore accelerating the phase transition. In (Bi$_{1-x}$In$_{x})_2$Se$_3$ we also find a tendency of In atoms to segregate.

preprint2013arXiv

Wannier Center Sheets in Topological Insulators

We argue that various kinds of topological insulators (TIs) can be insightfully characterized by an inspection of the charge centers of the hybrid Wannier functions, defined as the orbitals obtained by carrying out a Wannier transform on the Bloch functions in one dimension while leaving them Bloch-like in the other two. From this procedure, one can obtain the Wannier charge centers (WCCs) and plot them in the two-dimensional projected Brillouin zone. We show that these WCC sheets contain the same kind of topological information as is carried in the surface energy bands, with the crucial advantage that the topological properties of the bulk can be deduced from bulk calculations alone. The distinct topological behaviors of these WCC sheets in trivial, Chern, weak, strong, and crystalline TIs are first illustrated by calculating them for simple tight-binding models. We then present the results of first-principles calculations of the WCC sheets in the trivial insulator Sb$_2$Se$_3$, the weak TI KHgSb, and the strong TI Bi$_2$Se$_3$, confirming the ability of this approach to distinguish between different topological behaviors in an advantageous way.

preprint2012arXiv

Bloch-type Domain Walls in Rhombohedral BaTiO3

Ferroelectric domain walls (FDWs) are usually considered to be of Ising type, but there have been suggestions in recent years that Bloch-type FDWs are also possible in some cases, e.g., in the rhombohedral phase of BaTiO3. The mechanically compatible and electrically neutral FDWs in rhombohedral BaTiO3 are of 71-degree, 109-degree, and 180-degree type. We have investigated these FDWs based both on first-principles calculations and on a Ginzburg-Landau-Devonshire (GLD) model [P. Marton, I. Rychetsky, and J. Hlinka, Phys. Rev. B 81, 144125 (2010)]. The results from both approaches confirm the Ising nature of the 71-degree FDW and the Bloch nature of the 180-degree FDW, and predict both Ising-type and Bloch-type FDWs are possible for the 109-degree case. Considering the relatively small rhombohedral strain in BaTiO3, the competition between the energies of Bloch and Ising FDWs can be discussed in terms of a picture in which a Bloch wall is regarded as being composed of a pair of smaller-angle Ising ones. A reduction by 40% in the parameters describing the gradient term in the GLD model brings it into better agreement with the first-principles results for detailed properties such as the energies and widths of the FDWs.

preprint2012arXiv

Chern insulators from heavy atoms on magnetic substrates

We propose searching for Chern insulators by depositing atomic layers of elements with large spin-orbit coupling (e.g., Bi) on the surface of a magnetic insulator. We argue that such systems will typically have isolated surface bands with non-zero Chern numbers. If these overlap in energy, a metallic surface with large anomalous Hall conductivity (AHC) will result; if not, a Chern-insulator state will typically occur. Thus, our search strategy reduces to looking for examples having the Fermi level in a global gap extending across the entire Brillouin zone. We verify this search strategy and identify several candidate systems by using first-principles calculations to compute the Chern number and AHC of a large number of such systems on MnTe, MnSe, and EuS surfaces. Our search reveals several promising Chern insulators with gaps of up to 140\,meV.

preprint2012arXiv

Electrically driven octahedral rotations in SrTiO3 and PbTiO3

We investigate the oxygen octahedral rotations that occur in two perovskites, SrTiO3 and PbTiO3, as a function of applied three-dimensional electric displacement field, allowing us to map out the phase diagram of rotations in both the paraelectric and ferroelectric regions of the polar response. First-principles calculations at fixed electric displacement field are used to extract parameters of a Landau-Devonshire model that is analyzed to identify the phase boundaries between different rotational states. The calculations reveal a rich phase diagram of rotations versus applied field in both SrTiO3 and PbTiO3, although the details are quite different in the two cases.

preprint2012arXiv

Full magnetoelectric response of Cr2O3 from first principles

The linear magnetoelectric response of Cr2O3 at zero temperature is calculated from first principles by tracking the change in magnetization under a macroscopic electric field. Both the spin and the orbital contributions to the induced magnetization are computed, and in each case the response is decomposed into lattice and electronic parts. We find that the transverse response is dominated by the spin-lattice and spin-electronic contributions, whose calculated values are consistent with static and optical magnetoelectric measurements. In the case of the longitudinal response, orbital contributions dominate over spin contributions, but the net calculated longitudinal response remains much smaller than the experimentally measured one at low temperatures. We also discuss the absolute sign of the magnetoelectric coupling in the two time-reversed magnetic domains of Cr2O3.

preprint2012arXiv

Half-Heusler semiconductors as piezoelectrics

One of the central challenges in materials science is the design of functional and multifunctional materials, in which large responses are produced by applied fields and stresses. A rapidly developing paradigm for the rational design of such materials is based on the first-principles study of a large materials family, the perovskite oxides being the prototypical case. Specifically, first-principles calculations of structure and properties are used to explore the microscopic origins of the functional properties of interest and to search a large space of equilibrium and metastable phases to identify promising candidate systems. In this paper, we use a first-principles rational-design approach to demonstrate semiconducting half-Heusler compounds as a previously-unrecognized class of piezoelectric materials, and to provide guidance for the experimental realization and further investigation of high-performance materials suitable for practical applications.

preprint2012arXiv

Hexagonal $ABC$ as semiconducting ferroelectrics

We use a first-principles rational-design approach to identify a previously-unrecognized class of ferroelectric materials in the $P63mc$ LiGaGe structure type. We calculate structural parameters, polarization and ferroelectric well depths both for reported and as-yet hypothetical representatives of this class. Our results provide guidance for the experimental realization and further investigation of high-performance materials suitable for practical applications.

preprint2012arXiv

Maximally localized Wannier functions: Theory and applications

The electronic ground state of a periodic system is usually described in terms of extended Bloch orbitals, but an alternative representation in terms of localized "Wannier functions" was introduced by Gregory Wannier in 1937. The connection between the Bloch and Wannier representations is realized by families of transformations in a continuous space of unitary matrices, carrying a large degree of arbitrariness. Since 1997, methods have been developed that allow one to iteratively transform the extended Bloch orbitals of a first-principles calculation into a unique set of maximally localized Wannier functions, accomplishing the solid-state equivalent of constructing localized molecular orbitals, or "Boys orbitals" as previously known from the chemistry literature. These developments are reviewed here, and a survey of the applications of these methods is presented. This latter includes a description of their use in analyzing the nature of chemical bonding, or as a local probe of phenomena related to electric polarization and orbital magnetization. Wannier interpolation schemes are also reviewed, by which quantities computed on a coarse reciprocal-space mesh can be used to interpolate onto much finer meshes at low cost, and applications in which Wannier functions are used as efficient basis functions are discussed. Finally the construction and use of Wannier functions outside the context of electronic-structure theory is presented, for cases that include phonon excitations, photonic crystals, and cold-atom optical lattices.

preprint2012arXiv

Orthorhombic $ABC$ semiconductors as antiferroelectrics

We use a first-principles rational-design approach to identify a previously-unrecognized class of antiferroelectric materials in the $Pnma$ MgSrSi structure type. The MgSrSi structure type can be described in terms of antipolar distortions of the nonpolar $P6_{3}/mmc$ ZrBeSi structure type, and we find many members of this structure type are close in energy to the related polar $P6_{3}mc$ LiGaGe structure type, which includes many members we predict to be ferroelectric. We highlight known $ABC$ combinations in which this energy difference is comparable to the antiferroelectric-ferroelectric switching barrier of PbZrO$_{3}$. We calculate structural parameters and relative energies for all three structure types, both for reported and as-yet hypothetical representatives of this class. Our results provide guidance for the experimental realization and further investigation of high-performance materials suitable for practical applications.

preprint2012arXiv

Smooth gauge for topological insulators

We develop a technique for constructing Bloch-like functions for 2D Z_2-insulators (i.e., quantum spin-Hall insulators) that are smooth functions of k on the entire Brillouin-zone torus. As the initial step, the occupied subspace of the insulator is decomposed into a direct sum of two "Chern bands," i.e., topologically nontrivial subspaces with opposite Chern numbers. This decomposition remains robust independent of underlying symmetries or specific model features. Starting with the Chern bands obtained in this way, we construct a topologically nontrivial unitary transformation that rotates the occupied subspace into a direct sum of topologically trivial subspaces, thus facilitating a Wannier construction. The procedure is validated and illustrated by applying it to the Kane-Mele model.

preprint2012arXiv

Spin-phonon coupling effects in transition-metal perovskites:a DFT+$U$ and hybrid-functional study

Spin-phonon coupling effects, as reflected in phonon frequency shifts between ferromagnetic (FM) and G-type antiferromagnetic (AFM) configurations in cubic CaMnO$_3$, SrMnO$_3$, BaMnO$_3$, LaCrO$_3$, LaFeO$_3$ and La$_2$(CrFe)O$_6$, are investigated using density-functional methods. The calculations are carried out both with a hybrid-functional (HSE) approach and with a DFT+$U$ approach using a $U$ that has been fitted to HSE calculations. The phonon frequency shifts obtained in going from the FM to the AFM spin configuration agree well with those computed directly from the more accurate HSE approach, but are obtained with much less computational effort. We find that in the $A$MnO$_3$ materials class with $A$=Ca, Sr, and Ba, this frequency shift decreases as the A cation radius increases for the $Γ$ phonons, while it increases for R-point phonons. In La$M$O$_3$ with $M$=Cr, Fe, and Cr/Fe, the phonon frequencies at $Γ$ decrease as the spin order changes from AFM to FM for LaCrO$_3$ and LaFeO$_3$, but they increase for the double perovskite La$_2$(CrFe)O$_6$. We discuss these results and the prospects for bulk and superlattice forms of these materials to be useful as multiferroics.

preprint2012arXiv

Structure and energetics of a ferroelectric organic crystal of phenazine and chloranilic acid

We report first-principles calculations for a ferroelectric organic crystal of phenazine and chloranilic acid molecules. Weak intermolecular interactions are properly treated by using a second version of van der Waals density functional known as vdW-DF2 [K. Lee et al., Phys. Rev. B 82, 081101 (2010)]. Lattice constants, total energies, spontaneous electric polarizations, phonon modes and frequencies, and the energy barrier of proton transfer are calculated and compared with PBE and experiments whenever possible. We show that the donation of one proton from a chloranilic acid molecule to a neighboring phenazine molecule is energetically favorable. This proton transfer is the key structural change that breaks the centrosymmetry and leads to the ferroelectric structure. However, there is no unstable phonon associated with the proton transfer, and an energy barrier of 8 meV is found between the paraelectric and ferroelectric states.

preprint2011arXiv

Computing topological invariants without inversion symmetry

We consider the problem of calculating the weak and strong topological indices in noncentrosymmetric time-reversal (T) invariant insulators. In 2D we use a gauge corresponding to hybrid Wannier functions that are maximally localized in one dimension. Although this gauge is not smoothly defined on the two-torus, it respects the T symmetry of the system and allows for a definition of the Z_2 invariant in terms of time-reversal polarization. In 3D we apply the 2D approach to T-invariant planes. We illustrate the method with first-principles calculations on GeTe and on HgTe under [001] and [111] strain. Our approach differs from ones used previously for noncentrosymmetric materials and should be easier to implement in ab initio code packages.

preprint2011arXiv

First-principles theory of frozen-ion flexoelectricity

We demonstrate that the frozen-ion contribution to the flexoelectric coefficient is given solely in terms of the sum of third moments of the charge density distortions induced by atomic displacements, even for ferroelectric or piezoelectric materials. We introduce several practical supercell-based methods for calculating these coefficients from first principles, and demonstrate them by computing the coefficients for C, Si, MgO, NaCl, SrTiO$_3$, BaTiO$_3$, and PbTiO$_3$. Three important subtleties associated with pseudopotentials, the treatment of surfaces, and the calculation of transverse components are also discussed.

preprint2011arXiv

Interplay of epitaxial strain and rotations in PbTiO$_3$/PbZrO$_3$ superlattices from first principles

We present first-principles calculations of the structural phase behavior of the [1:1] \pzta\ superlattice and the \ptoa\ and \pzoa\ parent compounds as a function of in-plane epitaxial strain. A symmetry analysis is used to identify the phases and clarify how they arise from an interplay between different kinds of structural distortions, including out-of-plane and in-plane polar modes, rotation of oxygen octahedra around out-of-plane or in-plane axes, and an anti-polar mode. Symmetry-allowed intermode couplings are identified and used to elucidate the nature of the observed phase transitions. For the minimum-period [1:1] \pzta\ superlattice, we identify a sequence of three transitions that occur as the in-plane lattice constant is increased. All four of the phases involve substantial oxygen octahedral rotations, and an antipolar distortion is important in the high-tensile-strain phase. Inclusion of these distortions is found to be crucial for an accurate determination of the phase boundaries.

preprint2011arXiv

Mapping the energy surface of PbTiO3 in multidimensional electric-displacement space

In recent years, methods have been developed that allow first-principles electronic-structure calculations to be carried out under conditions of fixed electric field. For some purposes, however, it is more convenient to work at fixed electric displacement field. Initial implementations of the fixed-displacement-field approach have been limited to constraining the field along one spatial dimension only. Here, we generalize this approach to treat the full three-dimensional displacement field as a constraint, and compute the internal-energy landscape as a function of this multidimensional displacement-field vector. Using PbTiO3 as a prototypical system, we identify stable or metastable tetragonal, orthorhombic and rhombohedral structures as the displacement field evolves along [001], [110] and [111] directions, respectively. The energy minimum along [001] is found to be deeper than that along [110] or [111], as expected for a system having a tetragonal ground state. The barriers connecting these minima are found to be quite small, consistent with the current understanding that the large piezoelectric effects in PbTiO3 arise from the easy rotation of the polarization vector.

preprint2011arXiv

Orbital magnetoelectric coupling at finite electric field

We extend the band theory of linear orbital magnetoelectric coupling to treat crystals under finite electric fields. Previous work established that the orbital magnetoelectric response of a generic insulator at zero field comprises three contributions that were denoted as local circulation, itinerant circulation, and Chern-Simons. We find that the expression for each of them is modified by the presence of a dc electric field. Remarkably, the sum of the three correction terms vanishes, so that the total coupling is still given by the same formula as at zero field. This conclusion is confirmed by numerical tests on a tight-binding model, for which we calculate the field-induced change in the linear magnetoelectric coefficient.

preprint2011arXiv

Polar distortions in hydrogen bonded organic ferroelectrics

Although ferroelectric compounds containing hydrogen bonds were among the first to be discovered, organic ferroelectrics are relatively rare. The discovery of high polarization at room temperature in croconic acid [Nature \textbf{463}, 789 (2010)] has led to a renewed interest in organic ferroelectrics. We present an ab-initio study of two ferroelectric organic molecular crystals, 1-cyclobutene-1,2-dicarboxylic acid (CBDC) and 2-phenylmalondialdehyde (PhMDA). By using a distortion-mode analysis we shed light on the microscopic mechanisms contributing to the polarization, which we find to be as large as 14.3 and 7.0\,$μ$C/cm$^{2}$ for CBDC and PhMDA respectively. These results suggest that it may be fruitful to search among known but poorly characterized organic compounds for organic ferroelectrics with enhanced polar properties suitable for device applications.

preprint2011arXiv

Wannier-based calculation of the orbital magnetization in crystals

We present a first-principles scheme that allows the orbital magnetization of a magnetic crystal to be evaluated accurately and efficiently even in the presence of complex Fermi surfaces. Starting from an initial electronic-structure calculation with a coarse ab initio k-point mesh, maximally localized Wannier functions are constructed and used to interpolate the necessary k-space quantities on a fine mesh, in parallel to a previously-developed formalism for the anomalous Hall conductivity [X.Wang, J. Yates, I. Souza, and D. Vanderbilt, Phys. Rev. B 74, 195118 (2006)]. We formulate our new approach in a manifestly gauge-invariant manner, expressing the orbital magnetization in terms of traces over matrices in Wannier space. Since only a few (e.g., of the order of 20) Wannier functions are typically needed to describe the occupied and partially occupied bands, these Wannier matrices are small, which makes the interpolation itself very efficient. The method has been used to calculate the orbital magnetization of bcc Fe, hcp Co, and fcc Ni. Unlike an approximate calculation based on integrating orbital currents inside atomic spheres, our results nicely reproduce the experimentally measured ordering of the orbital magnetization in these three materials.

preprint2010arXiv

Chern-Simons orbital magnetoelectric coupling in generic insulators

We present a Wannier-based method to calculate the Chern-Simons orbital magnetoelectric coupling in the framework of first-principles density-functional theory. In view of recent developments in connection with strong Z2 topological insulators, we anticipate that the Chern-Simons contribution to the magnetoelectric coupling could, in special cases, be as large or larger than the total magnetoelectric coupling in known magnetoelectrics like Cr2O3. The results of our calculations for the ordinary magnetoelectrics Cr2O3, BiFeO3 and GdAlO3 confirm that the Chern-Simons contribution is quite small in these cases. On the other hand, we show that if the spatial inversion and time-reversal symmetries of the Z2 topological insulator Bi2Se3 are broken by hand, large induced changes appear in the Chern-Simons magnetoelectric coupling.

preprint2010arXiv

Interface enhancement of ferroelectricity in CaTiO$_3$/BaTiO$_3$ superlattices

We carry out first-principles calculations for CaTiO$_3$/BaTiO$_3$ superlattices with epitaxial strain corresponding to growth on a SrTiO$_3$ substrate, and consider octahedral rotations as well as ferroelectric distortions. The calculations are done as a function of electric displacement field, and both a macroscopic and a local electrostatic analysis are carried out. We find that strong octahedral rotations occur for TiO$_6$ octahedra sandwiched between CaO layers on both sides, but are strongly suppressed if either neighboring layer is a BaO layer. Due to the resulting enhancement of the ferroelectric instability in the BaO-neighboring octahedra, we find that overall the ferroelectric instability of the superlattice is enhanced by the interface. Thus, short-period superlattices in this system have a larger ferroelectric polarization than longer-period ones of the same average composition, contrary to the expected trend.

preprint2010arXiv

Orbital magnetoelectric coupling in band insulators

Magnetoelectric responses are a fundamental characteristic of materials that break time-reversal and inversion symmetries (notably multiferroics) and, remarkably, of "topological insulators" in which those symmetries are unbroken. Previous work has shown how to compute spin and lattice contributions to the magnetoelectric tensor. Here we solve the problem of orbital contributions by computing the frozen-lattice electronic polarization induced by a magnetic field. One part of this response (the "Chern-Simons term") can appear even in time-reversal-symmetric materials and has been previously shown to be quantized in topological insulators. In general materials there are additional orbital contributions to all parts of the magnetoelectric tensor; these vanish in topological insulators by symmetry and also vanish in several simplified models without time-reversal and inversion those magnetoelectric couplings were studied before. We give two derivations of the response formula, one based on a uniform magnetic field and one based on extrapolation of a long-wavelength magnetic field, and discuss some of the consequences of this formula.

preprint2010arXiv

Si-compatible candidates for high-K dielectrics with the Pbnm perovskite structure

We analyze both experimentally (where possible) and theoretically from first-principles the dielectric tensor components and crystal structure of five classes of Pbnm perovskites. All of these materials are believed to be stable on silicon and are therefore promising candidates for high-K dielectrics. We also analyze the structure of these materials with various simple models, decompose the lattice contribution to the dielectric tensor into force constant matrix eigenmode contributions, explore a peculiar correlation between structural and dielectric anisotropies in these compounds and give phonon frequencies and infrared activities of those modes that are infrared-active. We find that CaZrO_3, SrZrO_3, LaHoO_3, and LaYO_3 are among the most promising candidates for high-K dielectrics among the compounds we considered.

preprint2010arXiv

Theory of orbital magnetoelectric response

We extend the recently-developed theory of bulk orbital magnetization to finite electric fields, and use it to calculate the orbital magnetoelectric response of periodic insulators. Working in the independent-particle framework, we find that the finite-field orbital magnetization can be written as a sum of three gauge-invariant contributions, one of which has no counterpart at zero field. The extra contribution is collinear with and explicitly dependent on the electric field. The expression for the orbital magnetization is suitable for first-principles implementations, allowing to calculate the magnetoelectric response coefficients by numerical differentiation. Alternatively, perturbation-theory techniques may be used, and for that purpose we derive an expression directly for the linear magnetoelectric tensor by taking the first field-derivative analytically. Two types of terms are obtained. One, the `Chern-Simons' term, depends only on the unperturbed occupied orbitals and is purely isotropic. The other, `Kubo' term, involves the first-order change in the orbitals and gives isotropic as well as anisotropic contributions to the response. In ordinary magnetoelectric insulators both terms are generally present, while in strong Z2 topological insulators only the former type is allowed, and is quantized. In order to validate the theory we have calculated under periodic boundary conditions the linear magnetoelectric coupling for a 3-D tight-binding model of an ordinary magnetoelectric insulator, using both the finite-field and perturbation-theory expressions. The results are in excellent agreement with calculations on bounded samples.

preprint2010arXiv

Theory of prospective tetrahedral perovskite ferroelectrics

Using first-principles methods, we predict the energy landscape and ferroelectric states of double perovskites of the form AA$'$BB$'$O$_6$ in which the atoms on both the A and B sites are arranged in rock-salt order. While we are not aware of compounds that occur naturally in this structure, we argue that they might be realizable by directed synthesis. The high-symmetry structure formed by this arrangement belongs to the tetrahedral $F\bar{4}3m$ space group. If a ferroelectric instability occurs, the energy landscape will tend to have minima with the polarization along tetrahedral directions, leading to a rhombohedral phase, or along Cartesian directions, leading to an orthorhombic phase. We find that the latter scenario applies to CaBaTiZrO$_6$ and KCaZrNbO$_6$, which are weakly ferroelectric, and the former one applies to PbSnTiZrO$_6$, which is strongly ferroelectric. The results are modeled with a fourth- or fifth-order Landau-Devonshire expansion, providing good agreement with the first-principles calculations. Computations of zone-center soft modes are also carried out in order to characterize the polar and octahedral-rotation instabilities in more detail. Prospects for synthesis of ferroelectric materials belonging to this class are discussed.

preprint2010arXiv

Wannier representation of Z_2 topological insulators

We consider the problem of constructing Wannier functions for Z_2 topological insulators in two dimensions. It is well known that there is a topological obstruction to the construction of Wannier functions for Chern insulators, but it has been unclear whether this is also true for the Z_2 case. We consider the Kane-Mele tight-binding model, which exhibits both normal (Z_2-even) and topological (Z_2-odd) phases as a function of the model parameters. In the Z_2-even phase, the usual projection-based scheme can be used to build the Wannier representation. In the Z_2-odd phase, we do find a topological obstruction, but only if one insists on choosing a gauge that respects the time-reversal symmetry, corresponding to Wannier functions that come in time-reversal pairs. If instead we are willing to violate this gauge condition, a Wannier representation becomes possible. We present an explicit construction of Wannier functions for the Z_2-odd phase of the Kane-Mele model via a modified projection scheme followed by maximal localization, and confirm that these Wannier functions correctly represent the electric polarization and other electronic properties of the insulator.

preprint2009arXiv

Berry-phase theory of polar discontinuities at oxide-oxide interfaces

In the framework of the modern theory of polarization, we rigorously establish the microscopic nature of the electric displacement field D. In particular, we show that the longitudinal component of D is preserved at a coherent and insulating interface. To motivate and elucidate our derivation, we use the example of LAO/STO interfaces and superlattices, where the validity of the above conservation law is not immediately obvious. Our results generalize the "locality principle" of constrained-D density functional theory to the first-principles modeling of charge-mismatched systems.

preprint2009arXiv

First-principles modeling of ferroelectric capacitors via constrained-D calculations

First-principles modeling of ferroelectric capacitors presents several technical challenges, due to the coexistence of metallic electrodes, long-range electrostatic forces and short-range interface chemistry. Here we show how these aspects can be efficiently and accurately rationalized by using a finite-field density-functional theory formalism in which the fundamental electrical variable is the displacement field D. By performing calculations on model Pt/BaTiO3/Pt and Au/BaZrO3/Au capacitors we demonstrate how the interface-specific and bulk-specific properties can be identified and rigorously separated. Then, we show how the electrical properties of capacitors of arbitrary thickness and geometry (symmetric or asymmetric) can be readily reconstructed by using such information. Finally, we show how useful observables such as polarization and dielectric, piezoelectric and electrostrictive coefficients are easily evaluated as a byproduct of the above procedure. We apply this methodology to elucidate the relationship between chemical bonding, Schottky barriers and ferroelectric polarization at simple-metal/oxide interfaces. We find that BO2-electrode interfaces behave analogously to a layer of linear dielectric put in series with a bulk-like perovskite film, while a significant non-linear effect occurs at AO-electrode interfaces.

preprint2009arXiv

First-principles modeling of multferroic RMn$_2$O$_5$

We investigate the phase diagrams of RMn$_2$O$_5$ via a first-principles effective-Hamiltonian method. We are able to reproduce the most important features of the complicated magnetic and ferroelectric phase transitions. The calculated polarization as a function of temperature agrees very well with experiments. The dielectric-constant step at the commensurate-to -incommensurate magnetic phase transition is well reproduced. The microscopic mechanisms for the phase transitions are discussed.

preprint2009arXiv

First-principles study of high-field piezoelectricity in tetragonal PbTiO3

We calculate from first principles the nonlinear piezoelectric response of ferroelectric PbTiO3 for the case of a polarization-enhancing electric field applied along the tetragonal axis. We focus mainly on the case of fixed in-plane lattice constants, corresponding to epitaxially constrained thin films. We find that the dependence of the c/a ratio on electric field is almost linear in the range up to 500 MV/m, with little saturation. This result contrasts with expectations from Landau-Devonshire approaches based on experimental results obtained at lower fields, but is in qualitative agreement with a recent experiment in which higher fields were attained using pulsed-field methods. We also study cases in which the in-plane epitaxy constraint is removed, or an artificial negative pressure is applied, or both. These calculations demonstrate that PbTiO3 can show a strikingly non-linear piezoelectric response under modified elastic boundary conditions.

preprint2008arXiv

Electric displacement as the fundamental variable in electronic-structure calculations

Finite-field calculations in periodic insulators are technically and conceptually challenging, due to fundamental problems in defining polarization in extended solids. While significant progress has been made recently with the establishment of techniques to fix the electric field E or the macroscopic polarization P in first-principles calculations, both methods lack the ease of use and conceptual clarity of standard zero-field calculations. Here we develop a new formalism in which the electric displacement D, rather than E or P, is the fundamental electrical variable. Fixing D has the intuitive interpretation of imposing open-circuit electrical boundary conditions, which is particularly useful in studying ferroelectric systems. Furthermore, the analogy to open-circuit capacitors suggests an appealing reformulation in terms of free charges and potentials, which dramatically simplifies the treatment of stresses and strains. Using PbTiO3 as an example, we show that our technique allows full control over the electrical variables within the density functional formalism.

preprint2008arXiv

Enhancement of ferroelectricity at metal/oxide interfaces

By performing first-principles calculations on four capacitor structures based on BaTiO3 and PbTiO3, we determine the intrinsic interfacial effects that are responsible for the destabilization of the polar state in thin-film ferroelectric devices. We show that, contrary to the established interpretation, the interfacial capacitance depends crucially on the local chemical environment of the interface through the force constants of the metal-oxide bonds, and is not necessarily related to the bulk screening properties of the electrode material. In particular, in the case of interfaces between AO-terminated perovskites and simple metals, we demonstrate a novel mechanism of "interfacial ferroelectricity" that produces an overall enhancement of the ferroelectric instability of the film, rather than a suppression as is usually assumed. The resulting "negative dead layer" suggests a route to thin-film ferroelectric devices that are free of deleterious size effects.

preprint2008arXiv

Ideal barriers to polarization reversal and domain-wall motion in strained ferroelectric thin films

The ideal intrinsic barriers to domain switching in c-phase PbTiO_3 (PTO), PbZrO_3 (PZO), and PbZr_{1-x}Ti_xO_3 (PZT) are investigated via first-principles computational methods. The effects of epitaxial strain on the atomic structure, ferroelectric response, barrier to coherent domain reversal, domain-wall energy, and barrier to domain-wall translation are studied. It is found that PTO has a larger polarization, but smaller energy barrier to domain reversal, than PZO. Consequentially the idealized coercive field is over two times smaller in PTO than PZO. The Ti--O bond length is more sensitive to strain than the other bonds in the crystals. This results in the polarization and domain-wall energy in PTO having greater sensitivity to strain than in PZO. Two ordered phases of PZT are considered, the rock-salt structure and a (100) PTO/PZO superlattice. In these simple structures we find that the ferroelectric properties do not obey Vergard's law, but instead can be approximated as an average over individual 5-atom unit cells.

preprint2008arXiv

Theoretical investigation of polarization-compensated II-IV/I-V perovskite superlattices

Recent work suggested that head-to-head and tail-to-tail domain walls could be induced to form in ferroelectric superlattices by introducing compensating "delta doping" layers via chemical substitution in specified atomic planes [Phys. Rev. B 73, 020103(R), 2006]. Here we investigate a variation of this approach in which superlattices are formed of alternately stacked groups of II-IV and I-V perovskite layers, and the "polar discontinuity" at the II-IV/I-V interface effectively provides the delta-doping layer. Using first-principles calculations on SrTiO3/KNbO3 as a model system, we show that this strategy allows for the growth of a superlattice with stable polarized regions and large polarization discontinuities at the internal interfaces. We also generalize a Wannier-based definition of layer polarizations in perovskite superlattices [Phys. Rev. Lett. 97, 107602 (2006)] to the case in which some (e.g., KO or NbO2) layers are non-neutral, and apply this method to quantify the local variations in polarization in the proposed SrTiO3/KNbO3 superlattice system.

preprint2004arXiv

Amorphous ZrO2 from Ab-initio molecular dynamics: Structural, electronic and dielectric properties

Realistic models of amorphous ZrO2 are generated in a ``melt-and-quench'' fashion using ab-initio molecular dynamics in a plane-wave pseudopotential formulation of density-functional theory. The structural properties of the resulting amorphous models are analyzed, with special attention to coordination statistics. The vibrational and dielectric properties of one of these models are then investigated from first principles using linear-response methods. The electronic dielectric constant and Born effective charges are found to be very similar to those of the crystalline phases. Encouragingly, the predicted total static dielectric constant is about 22, comparable to that of the monoclinic phase. This work is motivated by the search for improved gate dielectric materials for sub-0.1 micron CMOS technology, and may also have implications for HfO2 and for silicates of ZrO2 and HfO2.

preprint1999arXiv

Theory of PbTiO3, BaTiO3, and SrTiO3 Surfaces

First-principles total-energy calculations are carried out for (001) surfaces of the cubic perovskite ATiO3 compounds PbTiO3, BaTiO3, and SrTiO3. Both AO-terminated and TiO2-terminated surfaces are considered, and fully-relaxed atomic configurations are determined. In general, BaTiO3 and SrTiO3 are found to have a rather similar behavior, while PbTiO3 is different in many respects because of the partially covalent character of the Pb-O bonds. PbTiO3 and BaTiO3 are ferroelectrics, and the influence of the surface upon the ferroelectric distortions is studied for the case of a tetragonal ferroelectric distortion parallel to the surface. The surface relaxation energies are found to be substantial, i.e., many times larger than the bulk ferroelectric well depth. Nevertheless, the influence of the surface upon the ferroelectric order parameter is modest, and is qualitatively as well as quantitatively different for the two materials. Surface energies and electronic properties are also computed. It is found that for BaTiO3 and SrTiO3 surfaces, both AO-terminated and TiO2-terminated surfaces can be thermodynamically stable, whereas for PbTiO3 only the PbO surface termination is stable.

preprint1995arXiv

First-principles theory of ferroelectric phase transitions for perovskites: The case of BaTiO3

We carry out a completely first-principles study of the ferroelectric phase transitions in BaTiO$_3$. Our approach takes advantage of two features of these transitions: the structural changes are small, and only low-energy distortions are important. Based on these observations, we make systematically improvable approximations which enable the parameterization of the complicated energy surface. The parameters are determined from first-principles total-energy calculations using ultra-soft pseudopotentials and a preconditioned conjugate-gradient scheme. The resulting effective Hamiltonian is then solved by Monte Carlo simulation. The calculated phase sequence, transition temperatures, latent heats, and spontaneous polarizations are all in good agreement with experiment. We find the transitions to be intermediate between order-disorder and displacive character. We find all three phase transitions to be of first order. The roles of different interactions are discussed.

David Vanderbilt

What is connected

Connect this record

See the researcher in context

Building this map preview

80 published item(s)

Direct visualization of surface spin-flip transition in MnBi4Te7

Electronic correlation in nearly free electron metals with beyond-DFT methods

Engineering magnetic topological insulators in Eu$_5M_2X_6$ Zintls

Vibrational fingerprints of ferroelectric hafnia

Vibrational properties of CuInP2S6 across the ferroelectric transition

Humble planar defects in SiGe nanopillars

Lattice dynamics and structural transition of the hyperhoneycomb iridate $β$-Li$_2$IrO$_3$ investigated by high-pressure Raman scattering

Orbital-selective Mott phase and non-Fermi liquid in FePS$_3$

Phonon spectrum of Pr$_2$Zr$_2$O$_7$ and Pr$_2$Ir$_2$O$_7$ as an evidence of coupling of the lattice with electronic and magnetic degrees of freedom

Quadrupole moments, edge polarizations, and corner charges in the Wannier representation

Axion coupling in the hybrid Wannier representation

Berry Flux Diagonalization: Application to Electric Polarization

Engineering Weyl phases and nonlinear Hall effects in T$_d$-MoTe$_2$

Gapless hinge states from adiabatic pumping of axion coupling

Magnetic phase transitions and spin density distribution in the molecular multiferroic GaV$_4$S$_8$ system

Molecular Mott state in the deficient spinel GaV$_4$S$_8$

Robust A-type order and spin-flop transition on the surface of the antiferromagnetic topological insulator MnBi$_2$Te$_4$

Systematic beyond-DFT study of binary transition metal oxides

Wannier90 as a community code: new features and applications

Ferroelectricity in corundum derivatives

Flux states and topological phases from spontaneous time-reversal symmetry breaking in CrSi(Ge)Te$_3$-based systems

Intertwined Rashba, Dirac and Weyl Fermions in Hexagonal Hyperferroelectrics

Temperature effects in the band structure of topological insulators

Topological nodal-line semimetals in alkaline-earth stannides, germanides and silicides

Tunable inverse topological heterostructure utilizing $(Bi_{1-x}In_{x})_{2}Se_{3}$

Chiral degeneracies and Fermi-surface Chern numbers in bcc Fe

Emergence of a Chern-insulating state from a semi-Dirac dispersion

First-principles theory of flexoelectricity

Gauge-discontinuity contributions to Chern-Simons orbital magnetoelectric coupling

Magnetic charges and magnetoelectricity in hexagonal rare-earth manganites and ferrites

Surface polarization and edge charges

Adiabatic Pumping of Chern-Simons Axion Coupling

Chern insulator at a magnetic rocksalt interface

Dynamical magnetic charges and linear magnetoelectricity

Quantum anomalous Hall phase in (001) double-perovskite monolayers via intersite spin-orbit coupling

Spin-orbit spillage as a measure of band inversion in insulators

Weyl semimetals from noncentrosymmetric topological insulators

Canonical magnetic insulators with isotropic magnetoelectric coupling

Comment on "Weyl fermions and the anomalous Hall effect in metallic ferromagnets"

Effective J=1/2 insulating state in Ruddlesden-Popper iridates: An LDA+DMFT study

First-principles theory and calculation of flexoelectricity

Hyperferroelectrics: proper ferroelectrics with persistent polarization

Topological phase transitions in (Bi$_{1-x}$In$_{x})_2$Se$_3$ and (Bi$_{1-x}$Sb$_{x})_2$Se$_3$

Wannier Center Sheets in Topological Insulators

Bloch-type Domain Walls in Rhombohedral BaTiO3

Chern insulators from heavy atoms on magnetic substrates

Electrically driven octahedral rotations in SrTiO3 and PbTiO3

Full magnetoelectric response of Cr2O3 from first principles

Half-Heusler semiconductors as piezoelectrics

Hexagonal $ABC$ as semiconducting ferroelectrics

Maximally localized Wannier functions: Theory and applications

Orthorhombic $ABC$ semiconductors as antiferroelectrics

Smooth gauge for topological insulators

Spin-phonon coupling effects in transition-metal perovskites:a DFT+$U$ and hybrid-functional study

Structure and energetics of a ferroelectric organic crystal of phenazine and chloranilic acid

Computing topological invariants without inversion symmetry

First-principles theory of frozen-ion flexoelectricity

Interplay of epitaxial strain and rotations in PbTiO$_3$/PbZrO$_3$ superlattices from first principles

Mapping the energy surface of PbTiO3 in multidimensional electric-displacement space

Orbital magnetoelectric coupling at finite electric field

Polar distortions in hydrogen bonded organic ferroelectrics

Wannier-based calculation of the orbital magnetization in crystals

Chern-Simons orbital magnetoelectric coupling in generic insulators

Interface enhancement of ferroelectricity in CaTiO$_3$/BaTiO$_3$ superlattices

Orbital magnetoelectric coupling in band insulators

Si-compatible candidates for high-K dielectrics with the Pbnm perovskite structure

Theory of orbital magnetoelectric response

Theory of prospective tetrahedral perovskite ferroelectrics

Wannier representation of Z_2 topological insulators

Berry-phase theory of polar discontinuities at oxide-oxide interfaces

First-principles modeling of ferroelectric capacitors via constrained-D calculations

First-principles modeling of multferroic RMn$_2$O$_5$

First-principles study of high-field piezoelectricity in tetragonal PbTiO3

Electric displacement as the fundamental variable in electronic-structure calculations