Showing 1–50 of 243 results for author: Vanderbilt, D

Inducing topological flat bands in bilayer graphene with electric and magnetic superlattices

Authors: Daniel Seleznev, Jennifer Cano, David Vanderbilt

Abstract: It was recently argued that Bernal stacked bilayer graphene (BLG) exposed to a 2D superlattice (SL) potential exhibits a variety of intriguing behaviors [Ghorashi et al., Phys. Rev. Lett. 130, 196201 (2023)]. Chief among them is the appearance of flat Chern bands that are favorable to the appearance of fractional Chern insulator states. Here, we explore the application of spatially periodic out-of… ▽ More It was recently argued that Bernal stacked bilayer graphene (BLG) exposed to a 2D superlattice (SL) potential exhibits a variety of intriguing behaviors [Ghorashi et al., Phys. Rev. Lett. 130, 196201 (2023)]. Chief among them is the appearance of flat Chern bands that are favorable to the appearance of fractional Chern insulator states. Here, we explore the application of spatially periodic out-of-plane orbital magnetic fields to the model of Ghorashi et al. to find additional means of inducing flat Chern bands. We focus on fields that vary on length scales much larger than the atomic spacing in BLG, generating what we refer to as magnetic SLs. The magnetic SLs we investigate either introduce no net magnetic flux to the SL unit cell, or a single quantum of flux. We find that magnetic SLs acting on their own can induce topological flat bands, but richer behavior, such as the appearance of flat and generic bands with high Chern numbers, can be observed when the magnetic SLs act in conjunction with commensurate electric SLs. Finally, we propose a method of generating unit-flux-quantum magnetic SLs along with concomitant electric SLs. The magnetic SL is generated by periodic arrays of flux vortices originating from type II superconductors, while the electric SL arises due to a magnetic SL-induced charge density on the surface of a magnetoelectric material. Tuning the vortex lattice and the magnetoelectric coupling permits control of both SLs, and we study their effects on the band structure of BLG. △ Less

Submitted 18 July, 2024; originally announced July 2024.

Magnetochiral Tunneling in Paramagnetic Co$_{1/3}$NbS$_2$

Authors: Seongjoon Lim, Sobhit Singh, Fei-Ting Huang, Shangke Pan, Kefeng Wang, Jaewook Kim, Jinwoong Kim, David Vanderbilt, Sang-Wook Cheong

Abstract: Electric currents have the intriguing ability to induce magnetization in nonmagnetic crystals with sufficiently low crystallographic symmetry. Some associated phenomena include the non-linear anomalous Hall effect in polar crystals and the nonreciprocal directional dichroism in chiral crystals when magnetic fields are applied. In this work, we demonstrate that the same underlying physics is also m… ▽ More Electric currents have the intriguing ability to induce magnetization in nonmagnetic crystals with sufficiently low crystallographic symmetry. Some associated phenomena include the non-linear anomalous Hall effect in polar crystals and the nonreciprocal directional dichroism in chiral crystals when magnetic fields are applied. In this work, we demonstrate that the same underlying physics is also manifested in the electronic tunneling process between the surface of a nonmagnetic chiral material and a magnetized scanning probe. In the paramagnetic but chiral metallic compound Co$_{1/3}$NbS$_2$, the magnetization induced by the tunneling current is shown to become detectable by its coupling to the magnetization of the tip itself. This results in a contrast across different chiral domains, achieving atomic-scale spatial resolution of structural chirality. To support the proposed mechanism, we used first-principles theory to compute the chirality-dependent current-induced magnetization and Berry curvature in the bulk of the material. Our demonstration of this magnetochiral tunneling effect opens up a new avenue for investigating atomic-scale variations in the local crystallographic symmetry and electronic structure across the structural domain boundaries of low-symmetry nonmagnetic crystals. △ Less

Submitted 3 February, 2024; originally announced February 2024.

Comments: 4 figures

Vacancy-induced tunable Kondo effect in twisted bilayer graphene

Authors: Yueqing Chang, Jinjing Yi, Ang-Kun Wu, Fabian B. Kugler, Eva Andrei, David Vanderbilt, Gabriel Kotliar, J. H. Pixley

Abstract: In single sheets of graphene, vacancy-induced states have been shown to host an effective spin-1/2 hole that can be Kondo-screened at low temperatures. Here, we show how these vacancy-induced impurity states survive in twisted bilayer graphene (TBG), which thus provides a tunable system to probe the critical destruction of the Kondo effect in pseudogap hosts. Ab-initio calculations and atomic-scal… ▽ More In single sheets of graphene, vacancy-induced states have been shown to host an effective spin-1/2 hole that can be Kondo-screened at low temperatures. Here, we show how these vacancy-induced impurity states survive in twisted bilayer graphene (TBG), which thus provides a tunable system to probe the critical destruction of the Kondo effect in pseudogap hosts. Ab-initio calculations and atomic-scale modeling are used to determine the nature of the vacancy states in the vicinity of the magic angle in TBG, demonstrating that the vacancy can be treated as a quantum impurity. Utilizing this insight, we construct an Anderson impurity model with a TBG host that we solve using the numerical renormalization group combined with the kernel polynomial method. We determine the phase diagram of the model and show how there is a strict dichotomy between vacancies in the AA / BB versus AB / BA tunneling regions. In AB / BA vacancies, we find that the Kondo temperature at the magic angle develops a broad distribution with a tail to vanishing temperatures due to multifractal wavefunctions at the magic angle. We argue that the scanning tunneling microscopy response in the vicinity of the vacancy can act as a non-trivial probe of both the critical single-particle states and the underlying many-body ground state in magic-angle TBG. △ Less

Submitted 14 December, 2023; originally announced December 2023.

Comments: 7 pages, 4 figures; including Supplementary Material (11 pages, 12 figures)

Non-collinear 2k antiferromagnetism in the Zintl semiconductor Eu$_5$In$_2$Sb$_6$

Authors: Vincent C. Morano, Jonathan Gaudet, Nicodemos Varnava, Tanya Berry, Thomas Halloran, Chris J. Lygouras, Xiaoping Wang, Christina M. Hoffman, Guangyong Xu, Jeffrey W. Lynn, Tyrel M. McQueen, David Vanderbilt, Collin L. Broholm

Abstract: Eu$_5$In$_2$Sb$_6$ is an orthorhombic non-symmorphic small band gap semiconductor with three distinct Eu$^{2+}$ sites and two low-temperature magnetic phase transitions. The material displays one of the greatest (negative) magnetoresistances of known stoichiometric antiferromagnets and belongs to a family of Zintl materials that may host an axion insulator. Using single crystal neutron diffraction… ▽ More Eu$_5$In$_2$Sb$_6$ is an orthorhombic non-symmorphic small band gap semiconductor with three distinct Eu$^{2+}$ sites and two low-temperature magnetic phase transitions. The material displays one of the greatest (negative) magnetoresistances of known stoichiometric antiferromagnets and belongs to a family of Zintl materials that may host an axion insulator. Using single crystal neutron diffraction, we show that the $T_{\mathrm{N1}}=14\mathrm{~K}$ second-order phase transition is associated with long-range antiferromagnetic order within the chemical unit cell $\left( k_1 = (000) \right)$. Upon cooling below $T_{\mathrm{N1}}$, the relative sublattice magnetizations of this structure vary until a second-order phase transition at $T_{\mathrm{N2}}=7\mathrm{~K}$ that doubles the unit cell along the $\hat{c}$ axis $\left( k_2 = \left(00\frac{1}{2}\right) \right)$. We show the anisotropic susceptibility and our magnetic neutron diffraction data are consistent with magnetic structures described by the $Γ_3$ irreducible representation with the staggered magnetization of the $k_1$ and $k_2$ components polarized along the $\hat{b}$ and $\hat{a}$ axis, respectively. As the $k_2$ component develops, the amplitude of the $k_1$ component is reduced, which indicates a 2k non-collinear magnetic structure. Density functional theory is used to calculate the energies of these magnetic structures and to show the $k_1$ phase is a metal so $T_{\mathrm{N1}}$ is a rare example of a unit-cell-preserving second-order phase transition from a paramagnetic semiconductor to an antiferromagnetic metal. DFT indicates the transition at $T_{\mathrm{N2}}$ to a doubled unit cell reduces the carrier density of the metal, which is consistent with resistivity data. △ Less

Submitted 1 November, 2023; originally announced November 2023.

Comments: 13 pages, 9 figures

Adiabatic dynamics of coupled spins and phonons in magnetic insulators

Authors: Shang Ren, John Bonini, Massimiliano Stengel, Cyrus E. Dreyer, David Vanderbilt

Abstract: In conventional \textit{ab initio} methodologies, phonons are calculated by solving equations of motion involving static interatomic force constants and atomic masses. The Born-Oppenheimer approximation, where all electronic degrees of freedom are assumed to adiabatically follow the nuclear dynamics, is also adopted. This approach does not fully account for the effects of broken time-reversal symm… ▽ More In conventional \textit{ab initio} methodologies, phonons are calculated by solving equations of motion involving static interatomic force constants and atomic masses. The Born-Oppenheimer approximation, where all electronic degrees of freedom are assumed to adiabatically follow the nuclear dynamics, is also adopted. This approach does not fully account for the effects of broken time-reversal symmetry in systems with magnetic order. Recent attempts to rectify this involve the inclusion of the velocity dependence of the interatomic forces in the equations of motion, which accounts for time-reversal symmetry breaking, and can result in chiral phonon modes with non-zero angular momentum even at the zone center. However, since the energy ranges of phonons and magnons typically overlap, the spins cannot be treated as adiabatically following the lattice degrees of freedom. Instead, phonon and spins must be treated on a similar footing. Focusing on zone-center modes, we propose a method involving Hessian matrices and Berry curvature tensors in terms of both phonon and spin degrees of freedom, and describe a first-principles methodology for calculating these. We then solve Lagrange's equations of motion to determine the energies and characters of the mixed excitations, allowing us to quantify, for example, the energy splittings between chiral pairs of phonons in some cases, and the degree of magnetically induced mixing between infrared and Raman modes in others. The approach is general, and can be applied to determine the adiabatic dynamics of any mixed set of slow variables. △ Less

Submitted 4 January, 2024; v1 submitted 11 July, 2023; originally announced July 2023.

Comments: Added a discussion section on the role of SOC and fixed some typos. This article has been accepted by Physical Review X

Journal ref: Phys. Rev. X 14, 011041 (2024)

Cyclic Ferroelectric Switching and Quantized Charge Transport in CuInP$_2$S$_6$

Authors: Daniel Seleznev, Sobhit Singh, John Bonini, Karin M. Rabe, David Vanderbilt

Abstract: The van der Waals layered ferroelectric CuInP$_2$S$_6$ has been found to exhibit a variety of intriguing properties arising from the fact that the Cu ions are unusually mobile in this system. While the polarization switching mechanism is usually understood to arise from Cu ion motion within the monolayers, a second switching path involving Cu motion across the van der Waals gaps has been suggested… ▽ More The van der Waals layered ferroelectric CuInP$_2$S$_6$ has been found to exhibit a variety of intriguing properties arising from the fact that the Cu ions are unusually mobile in this system. While the polarization switching mechanism is usually understood to arise from Cu ion motion within the monolayers, a second switching path involving Cu motion across the van der Waals gaps has been suggested. In this work, we perform zero-temperature first-principles calculations on such switching paths, focusing on two types that preserve the periodicity of the primitive unit cell: ``cooperative" paths preserving the system's glide mirror symmetry, and ``sequential" paths in which the two Cu ions in the unit cell move independently of each other. We find that CuInP$_2$S$_6$ features a rich and varied energy landscape, and that sequential paths are clearly favored energetically both for cross-gap and through-layer paths. Importantly, these segments can be assembled to comprise a globally insulating cycle with the out-of-plane polarization evolving by a quantum as the Cu ions shift to neighboring layers. In this sense, we argue that CuInP$_2$S$_6$ embodies the physics of a quantized adiabatic charge pump. △ Less

Submitted 1 May, 2023; originally announced May 2023.

Spin-phonon interactions and magnetoelectric coupling in Co$_4$$B_2$O$_9$ ($B$ = Nb, Ta)

Authors: K. Park, J. Kim, S. Choi, S. Fan, C. Kim, D. G. Oh, N. Lee, S. -W. Cheong, V. Kiryukhin, Y. J. Choi, D. Vanderbilt, J. H. Lee, J. L. Musfeldt

Abstract: In order to explore the consequences of spin-orbit coupling on spin-phonon interactions in a set of chemically-similar mixed metal oxides, we measured the infrared vibrational properties of Co$_4B_2$O$_9$ ($B$ = Nb, Ta) as a function of temperature and compared our findings with lattice dynamics calculations and several different models of spin-phonon coupling. Frequency vs. temperature trends for… ▽ More In order to explore the consequences of spin-orbit coupling on spin-phonon interactions in a set of chemically-similar mixed metal oxides, we measured the infrared vibrational properties of Co$_4B_2$O$_9$ ($B$ = Nb, Ta) as a function of temperature and compared our findings with lattice dynamics calculations and several different models of spin-phonon coupling. Frequency vs. temperature trends for the Co$^{2+}$ shearing mode near 150 cm$^{-1}$ reveal significant shifts across the magnetic ordering temperature that are especially large in relative terms. Bringing these results together and accounting for noncollinearity, we obtain spin-phonon coupling constants of -3.4 and -4.3 cm$^{-1}$ for Co$_4$Nb$_2$O$_9$ and the Ta analog, respectively. Analysis reveals that these coupling constants derive from interlayer (rather than intralayer) exchange interactions and that the interlayer interactions contain competing antiferromagnetic and ferromagnetic contributions. At the same time, beyond-Heisenberg terms are minimized due to fortuitous symmetry considerations, different than most other 4$d$- and 5$d$-containing oxides. Comparison with other contemporary oxides shows that spin-phonon coupling in this family of materials is among the strongest ever reported, suggesting an origin for magnetoelectric coupling. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Bismuth antiphase domain wall: A three-dimensional manifestation of the Su-Schrieffer-Heeger model

Authors: Jinwoong Kim, Cheng-Yi Huang, Hsin Lin, David Vanderbilt, Nicholas Kioussis

Abstract: The Su, Schrieffer and Heeger (SSH) model, describing the soliton excitations in polyacetylene due to the formation of antiphase domain walls (DW) from the alternating bond pattern, has served as a paradigmatic example of one-dimensional (1D) chiral topological insulators. While the SSH model has been realized in photonic and plasmonic systems, there have been limited analogues in three-dimensiona… ▽ More The Su, Schrieffer and Heeger (SSH) model, describing the soliton excitations in polyacetylene due to the formation of antiphase domain walls (DW) from the alternating bond pattern, has served as a paradigmatic example of one-dimensional (1D) chiral topological insulators. While the SSH model has been realized in photonic and plasmonic systems, there have been limited analogues in three-dimensional (3D) electronic systems, especially regarding the formation of antiphase DWs. Here, we propose that pristine bulk Bi, in which the dimerization of $(111)$ atomic layers renders alternating covalent and van der Waals bonding within and between successive $(111)$ bilayers, respectively, serves as a 3D analogue of the SSH model. First, we confirm that the two dimerized Bi structures belong to different Zak phases of 0 and $π$ by considering the parity eigenvalues and Wannier charge centers, while the previously reported bulk topological phases of Bi remain invariant under the dimerization reversal. Next, we demonstrate the existence of topologically non-trivial $(111)$ and trivial $(11\bar{2})$ DWs in which the number of in-gap DW states (ignoring spin) is odd and even respectively, and show how this controls the interlinking of the Zak phases of the two adjacent domains. Finally, we derive general criteria specifying when a DW of arbitrary orientation exhibits a $π$ Zak phase based on the flip of parity eigenvalues. An experimental realization of dimerization in Bi and the formation of DWs may be achieved via intense femtosecond laser excitations that can alter the interatomic forces and bond lengths. △ Less

Submitted 10 January, 2023; originally announced January 2023.

Journal ref: Phys. Rev. B 107, 045135 (2023)

Towards a theory of surface orbital magnetization

Authors: Daniel Seleznev, David Vanderbilt

Abstract: The theory of bulk orbital magnetization has been formulated both in reciprocal space based on Berry curvature and related quantities, and in real space in terms of the spatial average of a quantum mechanical local marker. Here we consider a three-dimensional antiferromagnetic material having a vanishing bulk but a nonzero surface orbital magnetization. We ask whether the surface-normal component… ▽ More The theory of bulk orbital magnetization has been formulated both in reciprocal space based on Berry curvature and related quantities, and in real space in terms of the spatial average of a quantum mechanical local marker. Here we consider a three-dimensional antiferromagnetic material having a vanishing bulk but a nonzero surface orbital magnetization. We ask whether the surface-normal component of the surface magnetization is well defined, and if so, how to compute it. As the physical observable corresponding to this quantity, we identify the macroscopic current running along a hinge shared by two facets. However, the hinge current only constrains the difference of the surface magnetizations on the adjoined facets, leaving a potential ambiguity. By performing a symmetry analysis, we find that only crystals exhibiting a pseudoscalar symmetry admit well-defined magnetizations at their surfaces at the classical level. We then explore the possibility of computing surface magnetization via a coarse-graining procedure applied to a quantum local marker. We show that multiple expressions for the local marker exist, and apply constraints to filter out potentially meaningful candidates. Using several tight-binding models as our theoretical test bed and several potential markers, we compute surface magnetizations for slab geometries and compare their predictions with explicit calculations of the macroscopic hinge currents of rod geometries. We find that only a particular form of the marker consistently predicts the correct hinge currents. △ Less

Submitted 15 February, 2023; v1 submitted 17 October, 2022; originally announced October 2022.

Comments: 26 pages, 12 figures

Magnetic phase crossover in strongly correlated EuMn2P2

Authors: Tanya Berry, Nicodemos Varnava, Dominic Ryan, Veronica Stewart, Riho Rästa, Ivo Heinmaa, Nitesh Kumar, Walter Schnelle, Rishi Bhandia, Christopher Pasco, N. P. Armitage, Raivo Stern, Claudia Felser, David Vanderbilt, Tyrel M. McQueen

Abstract: Strong electron correlations underlie a plethora of electronic and magnetic components and devices and are often used to identify and probe novel ground states in quantum materials. Herein we report a magnetic phase crossover in EuMn2P2, an insulator which shows Eu antiferromagnetism at TN=17K, but no phase transition attributed to Mn magnetism. The absence of a Mn magnetic phase transition contra… ▽ More Strong electron correlations underlie a plethora of electronic and magnetic components and devices and are often used to identify and probe novel ground states in quantum materials. Herein we report a magnetic phase crossover in EuMn2P2, an insulator which shows Eu antiferromagnetism at TN=17K, but no phase transition attributed to Mn magnetism. The absence of a Mn magnetic phase transition contrasts with the formation of long-range Mn order at T=130K in isoelectronic EuMn2Sb2 and EuMn2As2. Temperature-dependent specific heat and 31P NMR measurements provide evidence for the development of Mn magnetic correlations from T=250-100 K. Density functional theory calculations demonstrate an unusual sensitivity of the band structure to the details of the imposed Mn and Eu magnetic order, with antiferromagnetic Mn order required to recapitulate an insulating state. Our results imply a picture in which long range Mn magnetic order is suppressed by chemical pressure, but that magnetic correlations persist, narrowing bands and producing an insulating state. △ Less

Submitted 4 September, 2022; originally announced September 2022.

Frequency splitting of chiral phonons from broken time reversal symmetry in CrI$_3$

Authors: John Bonini, Shang Ren, David Vanderbilt, Massimiliano Stengel, Cyrus E. Dreyer, Sinisa Coh

Abstract: Conventional approaches for lattice dynamics based on static interatomic forces do not fully account for the effects of time-reversal-symmetry breaking in magnetic systems. Recent approaches to rectify this involve incorporating the first-order change in forces with atomic velocities under the assumption of adiabatic separation of electronic and nuclear degrees of freedom. In this work, we develop… ▽ More Conventional approaches for lattice dynamics based on static interatomic forces do not fully account for the effects of time-reversal-symmetry breaking in magnetic systems. Recent approaches to rectify this involve incorporating the first-order change in forces with atomic velocities under the assumption of adiabatic separation of electronic and nuclear degrees of freedom. In this work, we develop a first-principles method to calculate this velocity-force coupling in extended solids, and show via the example of ferromagnetic CrI$_3$ that, due to the slow dynamics of the spins in the system, the assumption of adiabatic separation can result in large errors for splittings of zone-center chiral modes. We demonstrate that an accurate description of the lattice dynamics requires treating magnons and phonons on the same footing. △ Less

Submitted 30 August, 2022; originally announced August 2022.

Integer quantum Hall effect and enhanced g-factor in quantum confined Cd3As2 films

Authors: Run Xiao, Junyi Zhang, Juan Chamorro, Jinwoong Kim, Tyrel M. McQueen, David Vanderbilt, Morteza Kayyalha, Yi Li, Nitin Samarth

Abstract: We investigate the integer quantum Hall effect in Cd3As2 thin films under conditions of strong to moderate quantum confinement (thicknesses of 10 nm, 12 nm, 15 nm). In all the films, we observe the integer quantum Hall effect in the spin-polarized lowest Landau level (filling factor ν = 1) and at spin-degenerate higher index Landau levels with even filling factors (ν = 2,4,6). With increasing quan… ▽ More We investigate the integer quantum Hall effect in Cd3As2 thin films under conditions of strong to moderate quantum confinement (thicknesses of 10 nm, 12 nm, 15 nm). In all the films, we observe the integer quantum Hall effect in the spin-polarized lowest Landau level (filling factor ν = 1) and at spin-degenerate higher index Landau levels with even filling factors (ν = 2,4,6). With increasing quantum confinement, we also observe a lifting of the Landau level spin degeneracy at ν = 3, manifest as the emergence of an anomaly in the longitudinal and Hall resistivity. Tight-binding calculations show that the enhanced g-factor likely arises from a combination of quantum confinement and corrections from nearby subbands. We also comment on the magnetic field induced transition from an insulator to a quantum Hall liquid when the chemical potential is near the charge neutrality point. △ Less

Submitted 18 July, 2022; originally announced July 2022.

Electronic structure of Humble defects in Ge and Ge$_{0.8}$Si$_{0.2}$

Authors: Shang Ren, Hongbin Yang, Sobhit Singh, Philip E. Batson, Eric L. Garfunkel, David Vanderbilt

Abstract: The group-IV diamond-structure elements are known to host a variety of planar defects, including {001} planar defects in C and {001}, {111} and {113} planar defects in Si and Ge. Among the {001} planar defects, the Humble defect, known for some time to occur in Ge, has recently also been observed in Si/Ge alloys, but the details of its electronic structure remain poorly understood. Here we perform… ▽ More The group-IV diamond-structure elements are known to host a variety of planar defects, including {001} planar defects in C and {001}, {111} and {113} planar defects in Si and Ge. Among the {001} planar defects, the Humble defect, known for some time to occur in Ge, has recently also been observed in Si/Ge alloys, but the details of its electronic structure remain poorly understood. Here we perform first-principles density functional calculations to study Humble defects in both Ge and Ge$_{0.8}$Si$_{0.2}$. We also measure the Si L$_{2,3}$-edge electron energy loss spectra both at the defect and in a bulk-like region far from the defect, and compare with theoretical calculations on corresponding Si sites in our first-principles calculations. We find that inclusion of core-hole effects in the theory is essential for reproducing the observed L$_{2,3}$ edge spectra, and that once they are included, the results provide a set of fingerprints for different types of local atomic bonding environments in Ge$_{0.8}$Si$_{0.2}$. Our first-principles calculations reveal that the Humble defects have a tendency to enlarge the electronic band gap, which may have potential uses in band engineering. The use of hybrid functionals for an improved description of the band gap in these systems is also discussed. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: 11 pages, 10 figures; Comments are welcome

Direct visualization of surface spin-flip transition in MnBi4Te7

Authors: Wenbo Ge, Jinwoong Kim, Ying-Ting Chan, David Vanderbilt, Jiaqiang Yan, Weida Wu

Abstract: 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 t… ▽ More 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. △ Less

Submitted 9 June, 2022; originally announced June 2022.

Comments: 6 pages, 4 figures

Dipolar Magnetic Interactions and A-type Antiferromagnetic Order in the Zintl Phase Insulator EuZn2P2

Authors: Tanya Berry, Veronica J. Stewart, Benjamin W. Y. Redemann, Chris Lygouras, Nicodemos Varnava, David Vanderbilt, Tyrel M. McQueen

Abstract: Zintl phases, containing strongly covalently bonded frameworks with separate ionically bonded ions, have emerged as a critical materials family in which to couple magnetism and strong spin-orbit coupling to drive diverse topological phases of matter. Here we report the single-crystal synthesis, magnetic, thermodynamic, transport, and theoretical properties of the Zintl compound EuZn2P2 that crysta… ▽ More Zintl phases, containing strongly covalently bonded frameworks with separate ionically bonded ions, have emerged as a critical materials family in which to couple magnetism and strong spin-orbit coupling to drive diverse topological phases of matter. Here we report the single-crystal synthesis, magnetic, thermodynamic, transport, and theoretical properties of the Zintl compound EuZn2P2 that crystallizes in the anti-La2O3 P-3m1 structure, containing triangular layers of Eu2+ ions. In-plane resistivity measurements reveal insulating behavior with an estimated bandgap of Eg=0.11eV. Comparing Eu magnetic ordering temperatures across trigonal EuM2X2 (M=divalent metal, X=pnictide) shows that EuZn2P2 exhibits the highest ordering temperature, with variations in TN correlating with changes in expected dipolar interaction strengths within and between layers and independent of the magnitude of electrical conductivity. These results provide experimental validation of the cytochemical intuition that the cation Eu2+ layers and the anionic (M2X2)2- framework can be treated as electronically distinct subunits, enabling further predictive materials design. △ Less

Submitted 23 March, 2022; originally announced March 2022.

Journal ref: Phys. Rev. B 106, 054420, 2022

Engineering magnetic topological insulators in Eu$_5M_2X_6$ Zintls

Authors: Nicodemos Varnava, Tanya Berry, Tyrel M. McQueen, David Vanderbilt

Abstract: 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 topologica… ▽ More 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. △ Less

Submitted 15 March, 2022; v1 submitted 11 March, 2022; originally announced March 2022.

Nonreciprocal directional dichroism at telecom wavelengths

Authors: K. Park, M. O. Yokosuk, M. Goryca, J. J. Yang, S. A. Crooker, S. -W. Cheong, K. Haule, D. Vanderbilt, H. -S. Kim, J. L. Musfeldt

Abstract: Magnetoelectrics with ultra-low symmetry and spin-orbit coupling are well known to display a number of remarkable properties including nonreciprocal directional dichroism. As a polar and chiral magnet, Ni$_3$TeO$_6$ is predicted to host this effect in three fundamentally different configurations, although only two have been experimentally verified. Inspired by the opportunity to unravel the struct… ▽ More Magnetoelectrics with ultra-low symmetry and spin-orbit coupling are well known to display a number of remarkable properties including nonreciprocal directional dichroism. As a polar and chiral magnet, Ni$_3$TeO$_6$ is predicted to host this effect in three fundamentally different configurations, although only two have been experimentally verified. Inspired by the opportunity to unravel the structure-property relations of such a unique light-matter interaction, we combined magneto-optical spectroscopy and first-principles calculations to reveal nonreciprocity in the toroidal geometry and compared our findings with the chiral configurations. We find that formation of Ni toroidal moments is responsible for the largest effects near 1.1 eV - a tendency that is captured by our microscopic model and computational implementation. At the same time, we demonstrate deterministic control of nonreciprocal directional dichroism in Ni$_3$TeO$_6$ across the entire telecom wavelength range. This discovery will accelerate the development of photonics applications that take advantage of unusual symmetry characteristics. △ Less

Submitted 4 March, 2022; originally announced March 2022.

Comments: accepted in npj. Quantum Materials

Band-Mott mixing hybridizes the gap in Fe$_2$Mo$_3$O$_8$

Authors: K. Park, G. L. Pascut, G. Khanal, M. O. Yokosuk, Xianghan Xu, Bin Gao, M. J. Gutmann, A. P. Litvinchuk, S. -W. Cheong, D. Vanderbilt, K. Haule, J. L. Musfeldt

Abstract: We combined optical spectroscopy and first principles electronic structure calculations to reveal the charge gap in the polar magnet Fe$_2$Mo$_3$O$_8$. Iron occupation on the octahedral site draws the gap strongly downward compared to the Zn parent compound, and subsequent occupation of the tetrahedral site creates a narrow resonance near the Fermi energy that draws the gap downward even further.… ▽ More We combined optical spectroscopy and first principles electronic structure calculations to reveal the charge gap in the polar magnet Fe$_2$Mo$_3$O$_8$. Iron occupation on the octahedral site draws the gap strongly downward compared to the Zn parent compound, and subsequent occupation of the tetrahedral site creates a narrow resonance near the Fermi energy that draws the gap downward even further. This resonance is a many-body effect that emanates from a flat valence band in a Mott-like state due to screening of the local moment - similar to expectations for a Zhang-Rice singlet, except that here, it appears in a semi-conductor. We discuss the unusual hybridization in terms of orbital occupation and character as well as the structure-property relationships that can be unveiled in various metal-substituted systems (Ni, Mn, Co, Zn). △ Less

Submitted 4 March, 2022; originally announced March 2022.

Journal ref: Phys. Rev. B. 104 (2021) 195143

Vibrational properties of CuInP2S6 across the ferroelectric transition

Authors: Sabine N. Neal, Sobhit Singh, Xiaochen Fang, Choongjae Won, Fei-ting Huang, Sang-Wook Cheong, Karin M. Rabe, David Vanderbilt, Janice L. Musfeldt

Abstract: 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 tra… ▽ More 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. △ Less

Submitted 31 January, 2022; originally announced February 2022.

Vibrational fingerprints of ferroelectric hafnia

Authors: Shiyu Fan, Sobhit Singh, Xianghan Xu, Kiman Park, Yubo Qi, S. W. Cheong, David Vanderbilt, Karin M. Rabe, J. L. Musfeldt

Abstract: 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 sing… ▽ More 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. △ Less

Submitted 29 January, 2022; originally announced January 2022.

Anomalous Residual Surface Conductivity in a Superconductor with Strong Spin-Orbit Coupling

Authors: J. R. Chamorro, P. Chauhan, C. Sun, N. Varnava, M. J. Winiarski, N. Ng, H. K. Vivanco, L. A. Pressley, C. M. Pasco, D. Vanderbilt, Yi Li, N. P. Armitage, T. M. McQueen

Abstract: Conventional BCS superconductors are expected to exhibit a conductivity with vanishing dissipation with decreasing temperature. While bulk physical properties measurements indicate PdPb$_{2}$ is a conventional superconductor with a $T_c$ of 3.0 K, measurements of surface impedance through the microwave cavity perturbation technique indicate a large, non-vanishing dissipative component below $T_c$… ▽ More Conventional BCS superconductors are expected to exhibit a conductivity with vanishing dissipation with decreasing temperature. While bulk physical properties measurements indicate PdPb$_{2}$ is a conventional superconductor with a $T_c$ of 3.0 K, measurements of surface impedance through the microwave cavity perturbation technique indicate a large, non-vanishing dissipative component below $T_c$ that is at odds with conventional superconductivity. We demonstrate PdPb$_2$ to be a possible topological superconductor with a fully gapped bulk and a dissipative Majorana fluid surface. △ Less

Submitted 9 May, 2022; v1 submitted 20 December, 2021; originally announced December 2021.

Humble planar defects in SiGe nanopillars

Authors: Hongbin Yang, Shang Ren, Sobhit Singh, Emily M. Turner, Kevin S. Jones, Philip E. Batson, David Vanderbilt, Eric Garfunkel

Abstract: 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… ▽ More 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. △ Less

Submitted 7 November, 2021; originally announced November 2021.

Polarization selectivity of aloof-beam electron energy-loss spectroscopy in one-dimensional ZnO nanorods

Authors: Yao-Wen Yeh, Sobhit Singh, David Vanderbilt, Philip E. Batson

Abstract: Orientation dependent electronic properties of wurtzite zinc oxide nanorods are characterized by aloof beam electron energy-loss spectroscopy (EELS) carried out in a scanning transmission electron microscope (STEM). The two key crystal orientation differentiating transitions specific to the in-plane (13.0 eV) and out-of-plane (11.2 eV) directions with respect to the wurtzite structure are examined… ▽ More Orientation dependent electronic properties of wurtzite zinc oxide nanorods are characterized by aloof beam electron energy-loss spectroscopy (EELS) carried out in a scanning transmission electron microscope (STEM). The two key crystal orientation differentiating transitions specific to the in-plane (13.0 eV) and out-of-plane (11.2 eV) directions with respect to the wurtzite structure are examined by first principles density-functional theory calculations. We note some degree of orientation dependence at the onset of direct band gap transition near 3.4 eV. We demonstrate that good polarization selectivity can be achieved by placing the electron probe at different locations around the specimen with increasing impact parameter while keeping the beam-specimen orientation fixed. The observed results are qualitatively elucidated in terms of the perpendicular electric fields generated by the fast electron (60 kV) used in the microscope. The fact that good polarization selectivity can be achieved by aloof beam EELS without the requirement of sample reorientation is an attractive aspect from the characterization method point of view in the STEM-EELS community. △ Less

Submitted 5 October, 2021; originally announced October 2021.

Comments: 9 pages and 9 figures

High-temperature phonon-mediated superconductivity in monolayer Mg2B4C2

Authors: Sobhit Singh, ? Aldo H. Romero, José D. Mella, Vitalie Eremeev, Enrique Muñoz Anastassia N. Alexandrova, Karin M. Rabe, David Vanderbilt, Francisco Muñoz

Abstract: A new two-dimensional material { Mg2B4C2, belonging to the family of the conventional superconductor MgB2, is theoretically predicted to exhibit superconductivity with critical temperature Tc estimated in the 47-48 K range (predicted using the McMillian-Allen-Dynes formula) without any tuning of external parameters such as doping, strain, or substrate-induced effects. The origin of such a high int… ▽ More A new two-dimensional material { Mg2B4C2, belonging to the family of the conventional superconductor MgB2, is theoretically predicted to exhibit superconductivity with critical temperature Tc estimated in the 47-48 K range (predicted using the McMillian-Allen-Dynes formula) without any tuning of external parameters such as doping, strain, or substrate-induced effects. The origin of such a high intrinsic Tc is ascribed to the presence of strong electron-phonon coupling and topological Dirac states (which are absent in MgB2) yielding a large density of states at the Fermi level. This material also features a nontrivial electronic band topology exhibiting Dirac points, practically gapless Dirac nodal lines, and topological nontrivial edge states. Consequently, it is a potential candidate for realization of topological superconductivity in 2D. This system is obtained after replacing the chemically active boron layers in MgB2 by chemically inactive boron-carbon layers. Hence, the surfaces of this material are inert. Our calculations confirm the stability of 2D Mg2B4C2. We also find that the key features of this material remain essentially unchanged when its thickness is increased by modestly increasing the number of inner MgB2 layers. △ Less

Submitted 3 October, 2021; originally announced October 2021.

Comments: 11 pages, 3 figures

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

Authors: Yuanyuan Xu, Huiyuan Man, Nan Tang, Takumi Ohtsuki, Santu Baidya, Satoru Nakatsuji, David Vanderbilt, Natalia Drichko

Abstract: 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… ▽ More 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. △ Less

Submitted 24 August, 2021; v1 submitted 3 August, 2021; originally announced August 2021.

Comments: 6 pages, 5 figures

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

Authors: Minsung Kim, Heung-Sik Kim, Kristjan Haule, David Vanderbilt

Abstract: 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 functi… ▽ More 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). △ Less

Submitted 11 June, 2021; originally announced June 2021.

Comments: 5 pages, 5 figures

Journal ref: Physical Review B 105, L041108 (2022)

Importance of dynamic lattice effects for crystal field excitations in quantum spin ice candidate Pr$_2$Zr$_2$O$_7$

Authors: Yuanyuan Xu, Huiyuan Man, Nan Tang, Santu Baidya, Hongbin Zhang, Satoru Nakatsuji, David Vanderbilt, Natalia Drichko

Abstract: Pr$_2$Zr$_2$O$_7$ is a pyrochlore quantum spin-ice candidate. Using Raman scattering spectroscopy we probe crystal electric field excitations of Pr$^{3+}$, and demonstrate the importance of their interactions with the lattice. We identify a vibronic interaction with a phonon that leads to a splitting of a doublet crystal field excitation at around 55~meV. We also probe a splitting of the non-Krame… ▽ More Pr$_2$Zr$_2$O$_7$ is a pyrochlore quantum spin-ice candidate. Using Raman scattering spectroscopy we probe crystal electric field excitations of Pr$^{3+}$, and demonstrate the importance of their interactions with the lattice. We identify a vibronic interaction with a phonon that leads to a splitting of a doublet crystal field excitation at around 55~meV. We also probe a splitting of the non-Kramers ground state doublet of Pr$^{3+}$ by observing a double line of the excitations to the first excited singlet state $E^0_g \rightarrow A_{1g}$. We show that the splitting has a strong temperature dependence, with the doublet structure most prominent between 50~K and 100~K, and the weight of one of the components strongly decreases on cooling. We suggest a static or dynamic deviation of Pr$^{3+}$ from the position in the ideal crystal structure can be the origin of the effect, with the deviation strongly decreasing at low temperatures. △ Less

Submitted 26 July, 2021; v1 submitted 15 February, 2021; originally announced February 2021.

Comments: 6 pages, 3 figures

Journal ref: Phys. Rev. B 104, 075125 (2021)

Mirror Chern numbers in the hybrid Wannier representation

Authors: Tomáš Rauch, Thomas Olsen, David Vanderbilt, Ivo Souza

Abstract: The topology of electronic states in band insulators with mirror symmetry can be classified in two different ways. One is in terms of the mirror Chern number, an integer that counts the number of protected Dirac cones in the Brillouin zone of high-symmetry surfaces. The other is via a $\mathbb{Z}_2$ index that distinguishes between systems that have a nonzero quantized orbital magnetoelectric coup… ▽ More The topology of electronic states in band insulators with mirror symmetry can be classified in two different ways. One is in terms of the mirror Chern number, an integer that counts the number of protected Dirac cones in the Brillouin zone of high-symmetry surfaces. The other is via a $\mathbb{Z}_2$ index that distinguishes between systems that have a nonzero quantized orbital magnetoelectric coupling ("axion-odd"), and those that do not ("axion-even"); this classification can also be induced by other symmetries in the magnetic point group, including time reversal and inversion. A systematic characterization of the axion $\mathbb{Z}_2$ topology has previously been obtained by representing the valence states in terms of hybrid Wannier functions localized along one chosen crystallographic direction, and inspecting the associated Wannier band structure. Here we focus on mirror symmetry, and extend that characterization to the mirror Chern number. We choose the direction orthogonal to the mirror plane as the Wannierization direction, and show that the mirror Chern number can be determined from the winding numbers of the touching points between Wannier bands on mirror-invariant planes, and from the Chern numbers of flat bands pinned to those planes. In this representation, the relation between the mirror Chern number and the axion $\mathbb{Z}_2$ index is readily established. The formalism is illustrated by means of $\textit{ab initio}$ calculations for SnTe in the monolayer and bulk forms, complemented by tight-binding calculations for a toy model. △ Less

Submitted 7 May, 2021; v1 submitted 26 January, 2021; originally announced January 2021.

Journal ref: Phys. Rev. B 103, 195103 (2021)

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

Authors: Sungkyun Choi, Heung-Sik Kim, Hun-Ho Kim, Aleksandra Krajewska, Gideok Kim, Matteo Minola, Tomohiro Takayama, Hidenori Takagi, Kristjan Haule, David Vanderbilt, Bernhard Keimer

Abstract: 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… ▽ More 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. △ Less

Submitted 19 January, 2021; originally announced January 2021.

Comments: 16 pages, 7 figures, 5 tables, Supplemental Material included

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

Authors: Subhasish Mandal, Kristjan Haule, Karin M. Rabe, David Vanderbilt

Abstract: 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 meta… ▽ More 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. △ Less

Submitted 29 July, 2022; v1 submitted 8 January, 2021; originally announced January 2021.

Lattice dynamics and magnetic exchange interactions in GeCo2O4, a spinel with S = 1/2 pyrochlore lattice

Authors: Prativa Pramanik, Sobhit Singh, Mouli Roy Chowdhury, Sayandeep Ghosh, Vasant Sathe, Karin M. Rabe, David Vanderbilt, Mohindar S. Seehra, Subhash Thota

Abstract: GeCo$_2$O$_4$ is a unique system in the family of cobalt spinels ACo$_2$O$_4$ (A= Sn, Ti, Ru, Mn, Al, Zn, Fe, etc.) in which magnetic Co ions stabilize on the pyrochlore lattice exhibiting a large degree of orbital frustration. Due to the complexity of the low-temperature antiferromagnetic (AFM) ordering and long-range magnetic exchange interactions, the lattice dynamics and magnetic structure of… ▽ More GeCo$_2$O$_4$ is a unique system in the family of cobalt spinels ACo$_2$O$_4$ (A= Sn, Ti, Ru, Mn, Al, Zn, Fe, etc.) in which magnetic Co ions stabilize on the pyrochlore lattice exhibiting a large degree of orbital frustration. Due to the complexity of the low-temperature antiferromagnetic (AFM) ordering and long-range magnetic exchange interactions, the lattice dynamics and magnetic structure of GeCo$_2$O$_4$ spinel has remained puzzling. To address this issue, here we present theoretical and experimental investigations of the highly frustrated magnetic structure, and the infrared (IR) and Raman-active phonon modes in the spinel GeCo$_2$O$_4$, which exhibits an AFM ordering below the Néel temperature $T_N$ ~21 K, followed by a cubic ($Fd{\bar 3}m$) to tetragonal ($I4_{1}/amd$) structural phase transition at $T_S$ ~16 K. Our density-functional theory (DFT+U) calculations reveal that one needs to consider magnetic-exchange interactions up to the third nearest neighbors to get an accurate description of the low-temperature AFM order in GeCo$_2$O$_4$. At room temperature three distinct IR-active modes ($T_{1u}$) are observed at frequencies 680, 413, and 325 cm$^{-1}$ along with four Raman-active modes $A_{1g}$, $T_{2g}(1)$, $T_{2g}(2)$, and $E_{g}$ at frequencies 760, 647, 550, and 308 cm$^{-1}$, respectively, which match reasonably well with our DFT+U calculated values. All the IR-active and Raman-active phonon modes exhibit signatures of moderate spin-phonon coupling. The temperature dependence of various parameters, such as the shift, width, and intensity, of the Raman-active modes, is also discussed. Noticeable changes around $T_N$ and $T_S$ are observed in the Raman line parameters of the $E_{g}$ and $T_{2g}$ modes, which are associated with the modulation of the Co-O bonds in CoO$_6$ octahedra during the excitations of these modes. △ Less

Submitted 29 July, 2021; v1 submitted 24 December, 2020; originally announced December 2020.

Comments: 9 figures and 11 pages

Journal ref: Phys. Rev. B 104, 014433 (2021)

Incommensurate magnetism mediated by Weyl fermions in NdAlSi

Authors: Jonathan Gaudet, Hung-Yu Yang, Santu Baidya, Baozhu Lu, Guangyong Xu, Yang Zhao, Jose A. Rodriguez, Christina M. Hoffmann, David E. Graf, Darius H. Torchinsky, Predrag Nikolić, David Vanderbilt, Fazel Tafti, Collin L. Broholm

Abstract: Emergent relativistic quasiparticles in Weyl semimetals are the source of exotic electronic properties such as surface Fermi arcs, the anomalous Hall effect, and negative magnetoresistance, all observed in real materials. Whereas these phenomena highlight the effect of Weyl fermions on the electronic transport properties, less is known about what collective phenomena they may support. Here, we rep… ▽ More Emergent relativistic quasiparticles in Weyl semimetals are the source of exotic electronic properties such as surface Fermi arcs, the anomalous Hall effect, and negative magnetoresistance, all observed in real materials. Whereas these phenomena highlight the effect of Weyl fermions on the electronic transport properties, less is known about what collective phenomena they may support. Here, we report a new Weyl semimetal, NdAlSi that offers an example. Using neutron diffraction, we report a long-wavelength magnetic order in NdAlSi whose periodicity is linked to the nesting vector between two topologically non-trivial Fermi pockets, which we characterize using density functional theory and quantum oscillation measurements. Our work provides a rare example of Weyl fermions driving collective magnetism. △ Less

Submitted 23 December, 2020; originally announced December 2020.

Comments: 6 main figures. Supplemental information also attached

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

Authors: Shang Ren, Ivo Souza, David Vanderbilt

Abstract: 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 tende… ▽ More 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. △ Less

Submitted 28 January, 2021; v1 submitted 26 October, 2020; originally announced October 2020.

Comments: 21 pages, 10 figures, typo fixed, references updated, published version

Journal ref: Phys. Rev. B 103, 035147 (2021)

First principles theory of Dirac semimetal Cd$_3$As$_2$ under Zeeman magnetic field

Authors: Santu Baidya, David Vanderbilt

Abstract: Time-reversal broken Weyl semimetals have attracted much attention recently, but certain aspects of their behavior, including the evolution of their Fermi surface topology and anomalous Hall conductivity with Fermi-level position, have remained underexplored. A promising route to obtain such materials may be to start with a nonmagnetic Dirac semimetal and break time-reversal symmetry via magnetic… ▽ More Time-reversal broken Weyl semimetals have attracted much attention recently, but certain aspects of their behavior, including the evolution of their Fermi surface topology and anomalous Hall conductivity with Fermi-level position, have remained underexplored. A promising route to obtain such materials may be to start with a nonmagnetic Dirac semimetal and break time-reversal symmetry via magnetic doping or magnetic proximity. Here we explore this scenario in the case of the Dirac semimetal Cd$_{3}$As$_{2}$, based on first-principles density-functional calculations and subsequent low-energy modeling of Cd$_{3}$As$_{2}$ in the presence of a Zeeman field applied along the symmetry axis. We clarify how each four$-$fold degenerate Dirac node splits into four Weyl nodes, two with chirality $\pm 1$ and two higher-order nodes with chirality $\pm 2$. Using a minimal \kdotp model Hamiltonian whose parameters are fit to the first-principles calculations, we detail the evolution of the Fermi surfaces and their Chern numbers as the Fermi energy is scanned across the region of the Weyl nodes at fixed Zeeman field. We also compute the intrinsic anomalous Hall conductivity as a function of Fermi-level position, finding a characteristic inverted-dome structure. Cd$_{3}$As$_{2}$ is especially well suited to such a study because of its high mobility, but the qualitative behavior revealed here should be applicable to other Dirac semimetals as well. △ Less

Submitted 24 August, 2020; originally announced August 2020.

Journal ref: Phys. Rev. B 102, 165115 (2020)

Controlling a quantum point junction on the surface of an antiferromagnetic topological insulator

Authors: Nicodemos Varnava, Justin H. Wilson, J. H. Pixley, David Vanderbilt

Abstract: The abstract notion of topology has led to profound insights into real materials. Notably, the surface and edges of topological materials can host physics, such as unidirectional charge or spin transport, that is unavailable in isolated one- and two-dimensional systems. However, to fully control the mixing and interference of edge-state wave functions, one needs robust and tunable junctions. We pr… ▽ More The abstract notion of topology has led to profound insights into real materials. Notably, the surface and edges of topological materials can host physics, such as unidirectional charge or spin transport, that is unavailable in isolated one- and two-dimensional systems. However, to fully control the mixing and interference of edge-state wave functions, one needs robust and tunable junctions. We propose to achieve this control using an antiferromagnetic topological insulator that supports two distinct types of gapless unidirectional channels on its surface, one from antiferromagnetic domain walls and the other from single-height steps. The distinct geometric nature of these edge modes allows them to intersect robustly to form quantum point junctions, and their presence at the surface makes them subject to control by magnetic and electrostatic tips like those used in scanning probe microscopes. Prospects for realizing such junctions are encouraged by recent material candidate proposals, potentially leading to exciting applications in quantum computing and sensing. △ Less

Submitted 7 August, 2020; originally announced August 2020.

The Joint Automated Repository for Various Integrated Simulations (JARVIS) for data-driven materials design

Authors: Kamal Choudhary, Kevin F. Garrity, Andrew C. E. Reid, Brian DeCost, Adam J. Biacchi, Angela R. Hight Walker, Zachary Trautt, Jason Hattrick-Simpers, A. Gilad Kusne, Andrea Centrone, Albert Davydov, Jie Jiang, Ruth Pachter, Gowoon Cheon, Evan Reed, Ankit Agrawal, Xiaofeng Qian, Vinit Sharma, Houlong Zhuang, Sergei V. Kalinin, Bobby G. Sumpter, Ghanshyam Pilania, Pinar Acar, Subhasish Mandal, Kristjan Haule , et al. (3 additional authors not shown)

Abstract: The Joint Automated Repository for Various Integrated Simulations (JARVIS) is an integrated infrastructure to accelerate materials discovery and design using density functional theory (DFT), classical force-fields (FF), and machine learning (ML) techniques. JARVIS is motivated by the Materials Genome Initiative (MGI) principles of developing open-access databases and tools to reduce the cost and d… ▽ More The Joint Automated Repository for Various Integrated Simulations (JARVIS) is an integrated infrastructure to accelerate materials discovery and design using density functional theory (DFT), classical force-fields (FF), and machine learning (ML) techniques. JARVIS is motivated by the Materials Genome Initiative (MGI) principles of developing open-access databases and tools to reduce the cost and development time of materials discovery, optimization, and deployment. The major features of JARVIS are: JARVIS-DFT, JARVIS-FF, JARVIS-ML, and JARVIS-Tools. To date, JARVIS consists of 40,000 materials and 1 million calculated properties in JARVIS-DFT, 1,500 materials and 110 force-fields in JARVIS-FF, and 25 ML models for material-property predictions in JARVIS-ML, all of which are continuously expanding. JARVIS-Tools provides scripts and workflows for running and analyzing various simulations. We compare our computational data to experiments or high-fidelity computational methods wherever applicable to evaluate error/uncertainty in predictions. In addition to the existing workflows, the infrastructure can support a wide variety of other technologically important applications as part of the data-driven materials design paradigm. The JARVIS datasets and tools are publicly available at the website: https://jarvis.nist.gov . △ Less

Submitted 11 July, 2021; v1 submitted 3 July, 2020; originally announced July 2020.

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

Authors: Rebecca L. Dally, William D. Ratcliff II, Lunyong Zhang, Heung-Sik Kim, Markus Bleuel, J. W. Kim, Kristjan Haule, David Vanderbilt, Sang-Wook Cheong, Jeffrey W. Lynn

Abstract: 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… ▽ More 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. △ Less

Submitted 1 July, 2020; originally announced July 2020.

Gapless hinge states from adiabatic pumping of axion coupling

Authors: Thomas Olsen, Tomáš Rauch, David Vanderbilt, Ivo Souza

Abstract: 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 pum… ▽ More 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. △ Less

Submitted 16 June, 2020; originally announced June 2020.

Comments: 11 pages, 11 figures

Journal ref: Phys. Rev. B 102, 035166 (2020)

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

Authors: Paul M. Sass, Jinwoong Kim, David Vanderbilt, Jiaqiang Yan, Weida Wu

Abstract: 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 corro… ▽ More 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. △ Less

Submitted 13 June, 2020; originally announced June 2020.

Comments: 11 pages, 4 figures

Journal ref: Physical Review Letters 125, 037201 (2020)

Berry Flux Diagonalization: Application to Electric Polarization

Authors: John Bonini, David Vanderbilt, Karin M. Rabe

Abstract: 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 c… ▽ More 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. △ Less

Submitted 16 July, 2020; v1 submitted 7 February, 2020; originally announced February 2020.

Journal ref: Phys. Rev. B 102, 045141 (2020)

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

Authors: Sobhit Singh, Jinwoong Kim, Karin M. Rabe, David Vanderbilt

Abstract: 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 d… ▽ More 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. △ Less

Submitted 21 July, 2020; v1 submitted 22 January, 2020; originally announced January 2020.

Journal ref: Phys. Rev. Lett. 125, 046402 (2020)

Axion coupling in the hybrid Wannier representation

Authors: Nicodemos Varnava, Ivo Souza, David Vanderbilt

Abstract: 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 insul… ▽ More 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. △ Less

Submitted 27 January, 2020; v1 submitted 26 December, 2019; originally announced December 2019.

Journal ref: Phys. Rev. B 101, 155130 (2020)

Influence of magnetic ordering on the spectral properties of binary transition metal oxides

Authors: Subhasish Mandal, Kristjan Haule, Karin M. Rabe, David Vanderbilt

Abstract: Using the ab initio embedded DMFT (eDMFT) approach, we study the effect of long-range magnetic ordering on the spectral properties in the binary transition metal oxides, and find that the most significant changes appear in the momentum resolved spectral functions, which sharpen into quite well-defined bands in the antiferromagnetic (AFM) phase. The strongest change across the transition is found a… ▽ More Using the ab initio embedded DMFT (eDMFT) approach, we study the effect of long-range magnetic ordering on the spectral properties in the binary transition metal oxides, and find that the most significant changes appear in the momentum resolved spectral functions, which sharpen into quite well-defined bands in the antiferromagnetic (AFM) phase. The strongest change across the transition is found at the topmost valence band edge (VBE), which is commonly associated with the Zhang-Rice bound state. This VBE state strengthens in the AFM phase, but only for the minority spin component, which is subject to stronger fluctuations. A similar hybridized VBE state also appears in the DFT single-particle description of the AFM phase, but gets much stronger and acquires a well-defined energy in the eDMFT description. △ Less

Submitted 19 September, 2019; originally announced September 2019.

Journal ref: Phys. Rev. B 100, 245109 (2019)

Polar and phase domain walls with conducting interfacial states in a Weyl semimetal MoTe2

Authors: Fei-Ting Huang, Seong Joon Lim, Sobhit Singh, Jinwoong Kim, Lunyong Zhang, Jae-Wook Kim, Ming-Wen Chu, Karin M. Rabe, David Vanderbilt, Sang-Wook Cheong

Abstract: Much of the dramatic growth in research on topological materials has focused on topologically protected surface states. While the domain walls of topological materials such as Weyl semimetals with broken inversion or time-reversal symmetry can provide a hunting ground for exploring topological interfacial states, such investigations have received little attention to date. Here, utilizing in-situ c… ▽ More Much of the dramatic growth in research on topological materials has focused on topologically protected surface states. While the domain walls of topological materials such as Weyl semimetals with broken inversion or time-reversal symmetry can provide a hunting ground for exploring topological interfacial states, such investigations have received little attention to date. Here, utilizing in-situ cryogenic transmission electron microscopy combined with first-principles calculations, we discover intriguing domain-wall structures in MoTe2, both between polar variants of the low-temperature(T) Weyl phase, and between this and the high-T high-order topological phase. We demonstrate how polar domain walls can be manipulated with electron beams and show that phase domain walls tend to form superlattice-like structures along the c axis. Scanning tunneling microscopy indicates a possible signature of a conducting hinge state at phase domain walls. Our results open avenues for investigating topological interfacial states and unveiling multifunctional aspects of domain walls in topological materials. △ Less

Submitted 8 August, 2019; originally announced August 2019.

Comments: 6 figures

Systematic beyond-DFT study of binary transition metal oxides

Authors: Subhasish Mandal, Kristjan Haule, Karin M. Rabe, David Vanderbilt

Abstract: 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… ▽ More 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. △ Less

Submitted 24 July, 2019; originally announced July 2019.

Journal ref: npj Comput. Mater. 5, 115 (2019)

Wannier90 as a community code: new features and applications

Authors: Giovanni Pizzi, Valerio Vitale, Ryotaro Arita, Stefan Blügel, Frank Freimuth, Guillaume Géranton, Marco Gibertini, Dominik Gresch, Charles Johnson, Takashi Koretsune, Julen Ibañez-Azpiroz, Hyungjun Lee, Jae-Mo Lihm, Daniel Marchand, Antimo Marrazzo, Yuriy Mokrousov, Jamal I. Mustafa, Yoshiro Nohara, Yusuke Nomura, Lorenzo Paulatto, Samuel Poncé, Thomas Ponweiser, Junfeng Qiao, Florian Thöle, Stepan S. Tsirkin , et al. (6 additional authors not shown)

Abstract: 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 n… ▽ More 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. △ Less

Submitted 23 July, 2019; originally announced July 2019.

Comments: 26 pages, 8 figures

Journal ref: J. Phys. Cond. Matt. 32, 165902 (2020)

Negative piezoelectric response of van der Waals layered bismuth tellurohalides

Authors: Jinwoong Kim, Karin M. Rabe, David Vanderbilt

Abstract: The polarization and piezoelectric response of the BiTe$X$ ($X$=Cl, Br, and I) layered tellurohalides are computed from first principles. The results confirm a mixed ionic-covalent character of the bonding, and demonstrate that the internal structure within each triple layer is only weakly affected by the external stress, while the changes in the charge distribution with stress produce a substanti… ▽ More The polarization and piezoelectric response of the BiTe$X$ ($X$=Cl, Br, and I) layered tellurohalides are computed from first principles. The results confirm a mixed ionic-covalent character of the bonding, and demonstrate that the internal structure within each triple layer is only weakly affected by the external stress, while the changes in the charge distribution with stress produce a substantial negative piezoelectric response. This suggests a new mechanism for negative piezoelectric response that should remain robust even in ultra-thin film form in this class of materials. △ Less

Submitted 17 June, 2019; originally announced June 2019.

Journal ref: Phys. Rev. B 100, 104115 (2019)

Ferromagnetic Anomalous Hall Effect in Cr-doped Bi$_2$Se$_3$ Thin Films via Surface-State Engineering

Authors: Jisoo Moon, Jinwoong Kim, Nikesh Koirala, Maryam Salehi, David Vanderbilt, Seongshik Oh

Abstract: The anomalous Hall effect (AHE) is a non-linear Hall effect appearing in magnetic conductors, boosted by internal magnetism beyond what is expected from the ordinary Hall effect. With the recent discovery of the quantized version of the AHE, the quantum anomalous Hall effect (QAHE), in Cr- or V-doped topological insulator (TI) (Sb,Bi)$_2$Te$_3$ thin films, the AHE in magnetic TIs has been attracti… ▽ More The anomalous Hall effect (AHE) is a non-linear Hall effect appearing in magnetic conductors, boosted by internal magnetism beyond what is expected from the ordinary Hall effect. With the recent discovery of the quantized version of the AHE, the quantum anomalous Hall effect (QAHE), in Cr- or V-doped topological insulator (TI) (Sb,Bi)$_2$Te$_3$ thin films, the AHE in magnetic TIs has been attracting significant interest. However, one of the puzzles in this system has been that while Cr- or V-doped (Sb,Bi)$_2$Te$_3$ and V-doped Bi$_2$Se$_3$ exhibit AHE, Cr-doped Bi$_2$Se$_3$ has failed to exhibit even ferromagnetic AHE, the expected predecessor to the QAHE, though it is the first material predicted to exhibit the QAHE. Here, we have successfully implemented ferromagnetic AHE in Cr-doped Bi$_2$Se$_3$ thin films by utilizing a surface state engineering scheme. Surprisingly, the observed ferromagnetic AHE in the Cr-doped Bi$_2$Se$_3$ thin films exhibited only positive slope regardless of the carrier type. We show that this sign problem can be explained by the intrinsic Berry curvature of the system as calculated from a tight-binding model combined with a first-principles method. △ Less

Submitted 3 May, 2019; originally announced May 2019.

Comments: 24 pages, 5 figures

Metric-wave approach to flexoelectricity within density-functional perturbation theory

Authors: Andrea Schiaffino, Cyrus E. Dreyer, David Vanderbilt, Massimiliano Stengel

Abstract: Within the framework of density functional perturbation theory (DFPT), we implement and test a novel "metric wave" response-function approach. It consists in the reformulation of an acoustic phonon perturbation in the curvilinear frame that is comoving with the atoms. This means that all the perturbation effects are encoded in the first-order variation of the real-space metric, while the atomic po… ▽ More Within the framework of density functional perturbation theory (DFPT), we implement and test a novel "metric wave" response-function approach. It consists in the reformulation of an acoustic phonon perturbation in the curvilinear frame that is comoving with the atoms. This means that all the perturbation effects are encoded in the first-order variation of the real-space metric, while the atomic positions remain fixed. This approach can be regarded as the generalization of the uniform strain perturbation of Hamann et al. [D. R. Hamann, X. Wu, K. M. Rabe, and D. Vanderbilt, Phys. Rev. B 71, 035117 (2005)] to the case of inhomogeneous deformations, and greatly facilitates the calculation of advanced electromechanical couplings such as the flexoelectric tensor. We demonstrate the accuracy of our approach with extensive tests on model systems and on bulk crystals of Si and SrTiO$_3$. △ Less

Submitted 30 November, 2018; originally announced November 2018.

Comments: 15 pages, 5 figures, 4 tables, 1 appendix

Journal ref: Phys. Rev. B 99, 085107 (2019)

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

Authors: Heung-Sik Kim, Kristjan Haule, David Vanderbilt

Abstract: 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 intr… ▽ More 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. △ Less

Submitted 10 August, 2020; v1 submitted 22 October, 2018; originally announced October 2018.

Comments: 5 pages, 5 figures

Journal ref: Phys. Rev. B 102, 081105 (2020)