Showing 1–50 of 186 results for author: Pickett, W E

Symmetry Enforced Fermi Surface Degeneracies Observed in Time-Reversal Symmetry-Breaking Superconductor LaNiGa$_2$

Authors: Matthew Staab, Robert Prater, Sudheer Sreedhar, Journey Byland, Eliana Mann, Davis Zackaria, Yunshu Shi, Henry J. Bowman, Andrew L. Stephens, Myung-Chul Jung, Antia S. Botana, Warren E. Pickett, Valentin Taufour, Inna Vishik

Abstract: LaNiGa$_2$ is superconductor that breaks time-reversal symmetry in the superconducting state without any known nearby magnetism. Recently, single crystals of LaNiGa$_2$ have been synthesized, revealing a nonsymmorphic Cmcm space group. Here, we report measurements of the electronic structure of LaNiGa$_2$ throughout the three-dimensional Brillouin zone (BZ) using angle-resolved photoemission spect… ▽ More LaNiGa$_2$ is superconductor that breaks time-reversal symmetry in the superconducting state without any known nearby magnetism. Recently, single crystals of LaNiGa$_2$ have been synthesized, revealing a nonsymmorphic Cmcm space group. Here, we report measurements of the electronic structure of LaNiGa$_2$ throughout the three-dimensional Brillouin zone (BZ) using angle-resolved photoemission spectroscopy (ARPES). Our findings show broad consistency with density functional theory (DFT) calculations and provide evidence for degeneracies in the electronic structure that are predicted from the space group. The calculations also predict four Fermi surfaces which cross the purported nodal plane and should therefore form two degenerate pairs. We report evidence for those predicted symmetry enforced degeneracies as well as accidental near degeneracies throughout the BZ. These degeneracies and near-degeneracies may play a role in the pairing mechanism of LaNiGa$_2$. Our results provide insight into the interplay between structure, Fermiology, and superconductivity in unconventional superconductors with nonsymmorphic space group. △ Less

Submitted 18 December, 2023; originally announced December 2023.

Comments: Main: 8 pages, 6 figures. SI: 6 pages, 8 figures

Comment on "Nontrivial Quantum Geometry and the Strength of Electron-Phonon Coupling", arXiv:2305.02340, J. Yu, C. J. Ciccarino, R. Bianco, I. Errea, P. Narang, B. A. Bernevig

Authors: Warren E. Pickett

Abstract: The titular manuscript invites a description of the background of research on MgB$_2$,whose unique electron-phonon coupling strength and record superconducting T$_c$ for such a material at ambient pressure was demonstrated by the combined efforts of several groups two decades ago in terms of conventional but highly developed Migdal-Eliashberg theory. This Comment provides some of the basic feature… ▽ More The titular manuscript invites a description of the background of research on MgB$_2$,whose unique electron-phonon coupling strength and record superconducting T$_c$ for such a material at ambient pressure was demonstrated by the combined efforts of several groups two decades ago in terms of conventional but highly developed Migdal-Eliashberg theory. This Comment provides some of the basic features of the theoretical understanding of MgB$_2$ and draws contrasts with the model of Yu {\it et al.} △ Less

Submitted 4 June, 2023; originally announced June 2023.

Comments: 5 pages of text and references

Room Temperature Superconductivity: the Roles of Theory and Materials Design

Authors: Warren E. Pickett

Abstract: For half a century after the discovery of superconductivity, materials exploration for better superconductors proceeded without knowledge of the underlying mechanism. The 1957 BCS theory cleared that up: the superconducting state occurs due to pairing of electrons over the Fermi surface. Over the following half century higher critical temperature T$_c$ was achieved only serendipitously as new mate… ▽ More For half a century after the discovery of superconductivity, materials exploration for better superconductors proceeded without knowledge of the underlying mechanism. The 1957 BCS theory cleared that up: the superconducting state occurs due to pairing of electrons over the Fermi surface. Over the following half century higher critical temperature T$_c$ was achieved only serendipitously as new materials were synthesized. Meanwhile the formal theory of phonon-coupled superconductivity at the material-dependent level became highly developed: given a known compound, its value of T$_c$, the superconducting gap function, and several other properties of the superconducting state became available independent of further experimental input. More recently, density functional theory based computational materials design has progressed to a predictive level -- new materials can be predicted on the basis of various numerical algorithms. Taken together, these capabilities enable theoretical prediction of new superconductors. Here the process that resulted in three new highest temperature superconductors, predicted numerically, confirmed experimentally -- SH$_3$, LaH$_{10}$, and YH$_9$ -- is recounted. These hydrides have T$_c$ in the 200-280K range at megabar pressures, and here the development will be chronicled. Current activities and challenges are discussed, together with Regularities in compressed hydrides that can guide further exploration. △ Less

Submitted 5 January, 2023; v1 submitted 12 April, 2022; originally announced April 2022.

Comments: 29 pages, 12 figures, 104 references

A $d^8$ anti-Hund's Singlet Insulator in an Infinite-layer Nickelate

Authors: Hyo-Sun Jin, W. E. Pickett, K. -W. Lee

Abstract: The status of nickelate superconductors in relation to cuprate high temperature superconductors is one of the concepts being discussed in high temperature superconductivity in correlated transition metal oxides. New additions to the class of infinite layer nickelates can provide essential input relating to connections or distinctions. A recently synthesized compound \bnoas, which contains isolated… ▽ More The status of nickelate superconductors in relation to cuprate high temperature superconductors is one of the concepts being discussed in high temperature superconductivity in correlated transition metal oxides. New additions to the class of infinite layer nickelates can provide essential input relating to connections or distinctions. A recently synthesized compound \bnoas, which contains isolated `infinite layer' NiO$_2$ planes, may lead to new insights. Our investigations have discovered that, at density functional theory mean field level, the ground state consists of an unusual $e_g$ singlet on the Ni$^{2+}$ ion arising from large but separate Mott insulating gaps in both $e_g$ orbitals, but with different, anti-Hund's, spin directions of their moments. This textured singlet incorporates at the least new physics, and potentially a new platform for nickelate superconductivity, which might be of an unconventional form for transition metal oxides due to the unconventional undoped state. We include in this paper a comparison of electronic structure parameters of Ba$_2$NiO$_2$(AgSe)$_2$ with a better characterized infinite layer nickelate LaNiO$_2$. We provide more analysis of the $d^8$ anti-Hund's singlet that emerges in this compound, and consider a minimally correlated wavefunction for this singlet in an itinerant background, and begin discussion of excitations -- real or virtual -- that may figure into new electronic phases. △ Less

Submitted 12 April, 2022; v1 submitted 28 December, 2021; originally announced December 2021.

Comments: 11 pages, 8 figures

Journal ref: J. Phys. Materials 5, 024008 (2022)

The dawn of the nickel age of superconductivity

Authors: Warren E. Pickett

Abstract: An invited Commentary in: Nature Reviews Physics 3, 7 (2021). An invited Commentary in: Nature Reviews Physics 3, 7 (2021). △ Less

Submitted 15 December, 2021; originally announced December 2021.

Comments: 3 pages, 1 figure. An invited Commentary in Nature Reviews Physics

Journal ref: Nature Reviews Physics 3, 7 (2021)

Low Valence Nickelates: Launching the Nickel Age of Superconductivity

Authors: Antia S. Botana, Kwan-Woo Lee, Michael R. Norman, Victor Pardo, Warren E. Pickett

Abstract: The discovery of superconductivity in thin films ($\sim$10 nm) of infinite-layer hole-doped NdNiO$_2$ has invigorated the field of high-temperature superconductivity research, reviving the debate over contrasting views that nickelates that are isostructural with cuprates are either (1) sisters of the high-temperature superconductors, or (2) that differences between nickel and copper at equal band… ▽ More The discovery of superconductivity in thin films ($\sim$10 nm) of infinite-layer hole-doped NdNiO$_2$ has invigorated the field of high-temperature superconductivity research, reviving the debate over contrasting views that nickelates that are isostructural with cuprates are either (1) sisters of the high-temperature superconductors, or (2) that differences between nickel and copper at equal band filling should be the focus of attention. Each viewpoint has its merits, and each has its limitations, suggesting that such a simple picture must be superseded by a more holistic comparison of the two classes. Several recent studies have begun this generalization, raising a number of questions without suggesting any consensus. In this paper, we organize the findings of the electronic structures of $n$-layered NiO$_2$ materials ($n$= 1 to $\infty$) to outline (ir)regularities and to make comparisons with cuprates, with the hope that important directions of future research will emerge. △ Less

Submitted 11 February, 2022; v1 submitted 1 November, 2021; originally announced November 2021.

Journal ref: Front. Phys. 9, 813532 (2022)

Dirac lines and loop at the Fermi level in the Time-Reversal Symmetry Breaking Superconductor LaNiGa$_2$

Authors: Jackson R. Badger, Yundi Quan, Matthew C. Staab, Shuntaro Sumita, Antonio Rossi, Kasey P. Devlin, Kelly Neubauer, Daniel S. Shulman, James C. Fettinger, Peter Klavins, Susan M. Kauzlarich, Dai Aoki, Inna M. Vishik, Warren E. Pickett, Valentin Taufour

Abstract: Unconventional superconductors have Cooper pairs with lower symmetries than in conventional superconductors. In most unconventional superconductors, the additional symmetry breaking occurs in relation to typical ingredients such as strongly correlated Fermi liquid phases, magnetic fluctuations, or strong spin-orbit coupling in noncentrosymmetric structures. In this article, we show that the time-r… ▽ More Unconventional superconductors have Cooper pairs with lower symmetries than in conventional superconductors. In most unconventional superconductors, the additional symmetry breaking occurs in relation to typical ingredients such as strongly correlated Fermi liquid phases, magnetic fluctuations, or strong spin-orbit coupling in noncentrosymmetric structures. In this article, we show that the time-reversal symmetry breaking in the superconductor LaNiGa$_2$ is enabled by its previously unknown topological electronic band structure. Our single crystal diffraction experiments indicate a nonsymmorphic crystal structure, in contrast to the previously reported symmorphic structure. The nonsymmorphic symmetries transform the $k_z=π/c$ plane of the Brillouin zone boundary into a node-surface. Band-structure calculations reveal that distinct Fermi surfaces become degenerate on the node-surface and form Dirac lines and a Dirac loop at the Fermi level. Two symmetry related Dirac points remain degenerate under spin-orbit coupling. ARPES measurements confirm the calculations and provide evidence for the Fermi surface degeneracies on the node-surface. These unique topological features enable an unconventional superconducting gap in which time-reversal symmetry can be broken in the absence of other typical ingredients. LaNiGa$_2$ is therefore a topological crystalline superconductor that breaks time-reversal symmetry without any overlapping magnetic ordering or fluctuations. Our findings will enable future discoveries of additional topological superconductors. △ Less

Submitted 14 September, 2021; originally announced September 2021.

Journal ref: Commun Phys 5, 22 (2022) [featured article]

MoB2 under Pressure: Superconducting Mo Enhanced by Boron

Authors: Yundi Quan, Kwan-woo Lee, Warren E. Pickett

Abstract: The discovery of the first high critical temperature (Tc) transition metal diboride superconductor, MgB2 structure alpha-MoB2 under pressure with Tc up to 32 K at 100 GPa, provides new input into some unexplained aspects of electron-phonon coupling in intermetallic compounds. We establish that MoB2 is a phonon-mediated superconductor but has little in common with MgB2 (Tc=40 K at zero pressure). M… ▽ More The discovery of the first high critical temperature (Tc) transition metal diboride superconductor, MgB2 structure alpha-MoB2 under pressure with Tc up to 32 K at 100 GPa, provides new input into some unexplained aspects of electron-phonon coupling in intermetallic compounds. We establish that MoB2 is a phonon-mediated superconductor but has little in common with MgB2 (Tc=40 K at zero pressure). MoB2 is a strongly metallic, three dimensional, multi-Fermi surface material, becoming of additional interest because it displays a frequency separation of Mo and B vibrations that mirrors that of metal superhydrides with Tc approaching room temperature. This separation, which is unusual in intermetallic compounds, allows an analysis separately for Mo and B providing, amongst the other parameters essential for understanding phonon coupling, the matrix elements for scattering by the individual atoms. Strong coupling (lambda(Mo)=1.48) is provided by Mo (total lambda=1.67). A factor of 15 weaker coupling to each B atom is compensated by that coupling being to mean high frequency modes around 85-90 meV (maximum of 140 meV), versus 18-20 meV for Mo. As a result, B enhances Tc by 43% over the Mo-only value, to 33 K, corresponding to the experimental value. These results provide a guideline for designing higher Tc materials from a cooperation of strong coupling from heavy atoms with weakly coupled light atoms. The new high Tc paradigm discovered here highlights the need for studying and engineering larger ionic scattering matrix elements. △ Less

Submitted 3 September, 2021; originally announced September 2021.

Comments: 13 pages, 9 figures, 56 references

Field-induced Bose-Einstein condensation and supersolid in the two-dimensional Kondo necklace

Authors: Wei-Lin Tu, Eun-Gook Moon, Kwan-Woo Lee, Warren E. Pickett, Hyun-Yong Lee

Abstract: The application of an external magnetic field of sufficient strength to a spin system composed of a localized singlet can overcome the energy gap and trigger bosonic condensation and so provide an alternative method to realize exotic phases of matter in real materials. Previous research has indicated that a spin Hamiltonian with on-site Kondo coupling may be the effective many-body Hamiltonian for… ▽ More The application of an external magnetic field of sufficient strength to a spin system composed of a localized singlet can overcome the energy gap and trigger bosonic condensation and so provide an alternative method to realize exotic phases of matter in real materials. Previous research has indicated that a spin Hamiltonian with on-site Kondo coupling may be the effective many-body Hamiltonian for $\text{Ba}_2\text{NiO}_2\text{(AgSe)}_2$ (BNOAS) and here we study such a Hamiltonian using a tensor network ansatz in two dimensions. Our results unveil a phase diagram which indicates the underlying phases of BNOAS. We propose, in response to the possible doping-induced superconductivity of BNOAS, a fermionic model for further investigation. We hope that our discovery can bring up further interest in both theoretical and experimental researches for related nickelate compounds. △ Less

Submitted 26 May, 2022; v1 submitted 25 July, 2021; originally announced July 2021.

Comments: 11 pages, 6 figures, 2 tables

Journal ref: Communications Physics 5, 130 (2022)

Two-band Conduction and Nesting Instabilities in Superconducting Ba$_2$CuO$_{3+δ}$: a First Principles Study

Authors: Hyo-Sun Jin, W. E. Pickett, K. -W. Lee

Abstract: First principles investigations of the high temperature superconducting system Ba$_2$CuO$_{3+δ}$, recently discovered at $δ\approx0.2$ at $T_c=70$ K, are applied to demonstrate the effects of oxygen ordering on the electronic and magnetic properties. The observed `highly over-doped' superconducting phase displays stretched Cu-planar oxygen O$_{\rm P}$ distances and anomalously shortened Cu-apical… ▽ More First principles investigations of the high temperature superconducting system Ba$_2$CuO$_{3+δ}$, recently discovered at $δ\approx0.2$ at $T_c=70$ K, are applied to demonstrate the effects of oxygen ordering on the electronic and magnetic properties. The observed `highly over-doped' superconducting phase displays stretched Cu-planar oxygen O$_{\rm P}$ distances and anomalously shortened Cu-apical O$_{\rm A}$ separations compared with other cuprates. The stoichiometric system $δ=0$, with its strongly one-dimensional (1D) Cu-O$_{\rm P}$ chain structure, when nonmagnetic shows 1D Fermi surfaces that lead, within density functional theory, to antiferromagnetic Cu-O$_{\rm P}$ chains (a spin-Peierls instability). Accounting for 1D fluctuations and small interchain coupling according to the theory of Schulz indicates this system, like Sr$_2$CuO$_3$, is near the 1D Luttinger-liquid quantum critical phase. The unusual Cu-O bond lengths per se have limited effects on other properties for $δ$=0. We find that a `doubled bilayer' structure of alternating Cu-O$_{\rm P}$ chains and wide rung Cu$_3$O$_4$ ladders is the energetically preferred one of three possibilities where the additional oxygen ions bridge Cu-O$_{\rm P}$ chains in the superconducting phase $δ=1/4$. Nominal formal valences of the three Cu sites are discussed. The six-fold (octahedral) site is the most highly oxidized, accepting somewhat more holes in the $d_{z^2}$ orbital than in the $d_{x^2-y^2}$ orbital. The implication is that two-band physics is involved in the pairing mechanism and the superconducting carriers. The Fermi surfaces of this metallic bilayer structure show both 1D and 2D strong (incipient) nesting instabilities, possibly accounting for the lack of clean single-phase samples based on this structure and suggesting importance for the pairing mechanism. △ Less

Submitted 30 August, 2021; v1 submitted 15 April, 2021; originally announced April 2021.

Comments: 13 pages

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

UTe$_2$: a nearly insulating half-filled $j=\frac{5}{2}$ $5f^3$ heavy fermion metal

Authors: Alexander B. Shick, Shin-ichi Fujimori, Warren E. Pickett

Abstract: Correlated band theory implemented as a combination of density functional theory with exact diagonalization [DFT+U(ED)] of the Anderson impurity term with Coulomb repulsion $U$ in the open 14-orbital $5f$ shell is applied to UTe$_2$. The small gap for $U$=0, evidence of the half-filled $j=\frac{5}{2}$ subshell of $5f^3$ uranium, is converted for $U$=3 eV to a flat band semimetal with small heavy-c… ▽ More Correlated band theory implemented as a combination of density functional theory with exact diagonalization [DFT+U(ED)] of the Anderson impurity term with Coulomb repulsion $U$ in the open 14-orbital $5f$ shell is applied to UTe$_2$. The small gap for $U$=0, evidence of the half-filled $j=\frac{5}{2}$ subshell of $5f^3$ uranium, is converted for $U$=3 eV to a flat band semimetal with small heavy-carrier Fermi surfaces that will make properties sensitive to pressure, magnetic field, and off-stoichiometry, as observed experimentally. The predicted Kondo temperature around 100 K matches the experimental values from resistivity. The electric field gradients for the two Te sites are calculated by DFT+U(ED) to differ by a factor of seven, indicating a strong site distinction, while the anisotropy factor $η=0.18$ is similar for all three sites. The calculated uranium moment $<M^2>^{1/2}$ of 3.5$μ_B$ is roughly consistent with the published experimental Curie-Weiss values of 2.8$μ_B$ and 3.3$μ_B$ (which are field-direction dependent), and the calculated separate spin and orbital moments are remarkably similar to Hund's rule values for an $f^3$ ion. The $U$=3 eV spectral density is compared with angle-integrated and angle-resolved photoemission spectra, with agreement that there is strong $5f$ character at, and for several hundred meV below, the Fermi energy. Our results support the picture that the underlying ground state of UTe$_2$ is that of a half-filled $j=\frac{5}{2}$ subshell with two half-filled $m_j=\pm\frac{1}{2}$ orbitals forming a narrow gap by hybridization, then driven to a conducting state by configuration mixing (spin-charge fluctuations). UTe$_2$ displays similarities to UPt$_3$ with its $5f$ dominated Fermi surfaces rather than a strongly localized Kondo lattice system. △ Less

Submitted 21 March, 2021; originally announced March 2021.

Comments: 14 pages, 6 figures

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

Proposed ordering of textured spin singlets in a bulk infinite layer nickelate

Authors: Hyo-Sun Jin, W. E. Pickett, K. -W. Lee

Abstract: The infinite-layer structure nickelate Ba$_2$NiO$_2$(AgSe)$_2$ (BNOAS) with $d^8$ Ni ions and a peculiar susceptibility $χ(T)$ is studied with correlated density functional methods. The overriding feature of the calculations is violation of Hund's rule coupled with complete but unconventional spin-orbital polarization, leading to an unexpected low spin $^1B_1$, "off-diagonal singlet" (ODS) texture… ▽ More The infinite-layer structure nickelate Ba$_2$NiO$_2$(AgSe)$_2$ (BNOAS) with $d^8$ Ni ions and a peculiar susceptibility $χ(T)$ is studied with correlated density functional methods. The overriding feature of the calculations is violation of Hund's rule coupled with complete but unconventional spin-orbital polarization, leading to an unexpected low spin $^1B_1$, "off-diagonal singlet" (ODS) textured by an internal orbital structure of compensating $d_{x^2-y^2}^{\uparrow}$ and $d_{z^2}^{\downarrow}$ spins. This unconventional configuration has lower energy than conventional high-spin or low-spin alternatives. An electronic transition is obtained at a critical Ni-O separation $d_c^{Ni-O}=$2.03 Å, above which Ni becomes magnetic in square planar NiO$_2$ compounds. We propose scenarios for the signature of magnetic reconstruction in $χ(T)$ at $T_{m}$=130 K without any Curie-Weiss background (no moment) that invoke ordering of Ni $d^8$ moieties that are largely this generalized Kondo singlet. The underlying physics of this system is modeled by a Kondo sieve model (2D Kondo necklace) of a "Kondo" $d_{z^2}$ spin on each site, coupled to a $d_{x^2-y^2}$ spin that is itself strongly coupled to neighboring like-spins within the layer. The observed magnetic order places BNOAS below the quantum critical point of the Kondo sieve model, providing a realization of the previously unreported long-range ordered near-singlet weak antiferromagnetic phase. We propose electron doping experiments that would drive the system toward a $d^{9-δ}$ configuration and possible superconductivity with similarity to the recently reported Ba$_2$CuO$_{3.2}$ that superconducts at 73 K. △ Less

Submitted 12 August, 2020; v1 submitted 6 July, 2020; originally announced July 2020.

Comments: 12 pages (minor corrections)

Journal ref: Phys. Rev. Research 2, 033197 (2020)

Fluctuation-frustrated flat band instabilities in NdNiO2

Authors: Mi-Young Choi, W. E. Pickett, K. -W. Lee

Abstract: The discovery that Nd$_{1-x}$Sr$_x$NiO$_2$, with the CaCuO$_2$ infinite-layer structure, superconducts up to 15 K around the hole-doping level $x$=0.2 raises the crucial question of its fundamental electronic and magnetic processes. The unexplained basic feature that we address is that, for $x$=0 and as opposed to strongly antiferromagnetic (AFM) CaCuO$_2$, NdNiO$_2$ with the same structure and fo… ▽ More The discovery that Nd$_{1-x}$Sr$_x$NiO$_2$, with the CaCuO$_2$ infinite-layer structure, superconducts up to 15 K around the hole-doping level $x$=0.2 raises the crucial question of its fundamental electronic and magnetic processes. The unexplained basic feature that we address is that, for $x$=0 and as opposed to strongly antiferromagnetic (AFM) CaCuO$_2$, NdNiO$_2$ with the same structure and formal $d^9$ configuration does not undergo AFM order. We study this issue not in the conventional manner, as energetically unfavored or as frustrated magnetic order, but as an instability of the AFM phase itself. We are able to obtain the static AFM ordered state, but find that a flat-band, one-dimensional-like van Hove singularity (vHs) is pinned to the Fermi level. This situation is unusual in a non-half-filled, effectively two-band system. The vHs makes the AFM phase unstable to spin-density disproportionation, breathing and half-breathing lattice distortions, and (innate or parasitic) charge-density disproportionation. These flat-band instabilities, distant relatives of single band cuprate models, thereby inhibit but do not eliminate incipient AFM tendencies at low temperature. The primary feature is that a pair of active bands ($d_{x^2-y^2}$, $d_{z^2})$ eliminate half-filled physics and, due to instabilities, preclude the AFM phase seen in CaCuO$_2$. This strongly AFM correlated, conducting spin-liquid phase with strong participation of the Ni $d_{z^2}$ orbital, forms the platform for superconductivity in NdNiO$_2$. △ Less

Submitted 21 September, 2020; v1 submitted 6 May, 2020; originally announced May 2020.

Comments: 8 pages

Journal ref: Phys. Rev. Research 2, 033445 (2020)

Accidental Degeneracy in k-space, Geometrical Phase, and the Perturbation of $π$ by Spin-orbit Interactions

Authors: Philip B. Allen, Warren E. Pickett

Abstract: Since closed lines of {\it accidental} electronic degeneracies were demonstrated to be possible, even frequent, by Herring in 1937, no further developments arose for eight decades. The earliest report of such a nodal loop in a real material -- aluminum -- is recounted and elaborated on. Nodal loop semimetals have become a focus of recent activity, with emphasis on other issues. Band degeneracies a… ▽ More Since closed lines of {\it accidental} electronic degeneracies were demonstrated to be possible, even frequent, by Herring in 1937, no further developments arose for eight decades. The earliest report of such a nodal loop in a real material -- aluminum -- is recounted and elaborated on. Nodal loop semimetals have become a focus of recent activity, with emphasis on other issues. Band degeneracies are, after all, the origin of topological phases in crystalline materials. Spin-orbit interaction lifts accidental band degeneracies, with the resulting spectrum being provided here. The geometric phase $γ(C)=\pmπ$ for circuits $C$ surrounding a line of such degeneracy cannot survive completely unchanged. The change depends on how the spin is fixed during adiabatic evolution. For spin fixed along the internal spin-orbit field, $γ(C)$ decreases to zero as the circuit collapses around the line of lifted degeneracy. For spin fixed along a perpendicular axis, the conical intersection persists and $γ(C)=\pmπ$ is unchanged. △ Less

Submitted 3 April, 2020; originally announced April 2020.

Comments: 6 pages, 2 figures. arXiv admin note: substantial text overlap with arXiv:0709.1457

Journal ref: Physica C 549, 102 (2018)

Topological and Thermoelectric Properties of Double Antiperovskite Pnictides

Authors: Wen Fong Goh, Warren E. Pickett

Abstract: Doubling the perovskite cell (double perovskite) has been found to open new possibilities for engineering functional materials, magnetic materials in particular. This route should be applicable to the antiperovskite (aPV) class. In the pnictide based double aPV (2aPV) class introduced here magnetism is very rare, and we address them as new topological materials, possibly with thermoelectric intere… ▽ More Doubling the perovskite cell (double perovskite) has been found to open new possibilities for engineering functional materials, magnetic materials in particular. This route should be applicable to the antiperovskite (aPV) class. In the pnictide based double aPV (2aPV) class introduced here magnetism is very rare, and we address them as new topological materials, possibly with thermoelectric interest. We have found that the 2aPV supercell provides a systematically larger band gap that can serve to inhibit bulk conductivity, and also large spin-orbit coupling (SOC) for band inversion. We present examples from a broad study of double antiperovskites focusing on the X$_6$AA$'$B$_2$ configuration, where X is the alkaline earth element and A and B are the group 5A pnictogens. We find that an "extended s" state at the valence band minimum, described alternatively as a cation valence state or a modulated interstitial planewave state, plays a crucial role in both topological and thermoelectric properties. Several of these compounds may house topological phases, while transport calculations indicate they may also find themselves useful in thermoelectric applications. △ Less

Submitted 1 April, 2020; originally announced April 2020.

Role of $4f$ states in infinite-layer NdNiO$_2$

Authors: Mi-Young Choi, K. -W. Lee, W. E. Pickett

Abstract: Atomic $4f$ states have been found to be essential players in the physical behavior of lanthanide compounds, at the Fermi level $E_F$ as in the proposed topological Kondo insulator SmB$_6$, or further away as in the magnetic superconductor system ${\cal R}$Ni$_2$B$_2$C (${\cal R}$=rare earth ion) and in Y$_{1-x}$Pr$_x$Ba$_2$Cu$_3$O$_7$, where the $4f$ shell of Pr has a devastating effect on superc… ▽ More Atomic $4f$ states have been found to be essential players in the physical behavior of lanthanide compounds, at the Fermi level $E_F$ as in the proposed topological Kondo insulator SmB$_6$, or further away as in the magnetic superconductor system ${\cal R}$Ni$_2$B$_2$C (${\cal R}$=rare earth ion) and in Y$_{1-x}$Pr$_x$Ba$_2$Cu$_3$O$_7$, where the $4f$ shell of Pr has a devastating effect on superconductivity. In hole-doped ${\cal R}$NiO$_2$, the ${\cal R}$=Nd member is found to be superconducting while ${\cal R}$=La is not, in spite of the calculated electronic structures being nearly identical. We report first principles results that indicate that the Nd $4f$ moment affects states at $E_F$ in infinite-layer NdNiO$_2$, an effect that will not occur for LaNiO$_2$. Treating 20% hole-doping in the virtual crystal approach indicates that 0.15 holes empty the $Γ$-centered Nd-derived electron pocket while leaving the other electron pocket unchanged; hence Ni only absorbs 0.05 holes; the La counterpart would behave similarly. However, coupling of $4f$ states to the electron pockets at $E_F$ arises through the Nd intra-atomic $4f-5d$ exchange coupling $K\approx 0.5$ eV and is ferromagnetic (FM), i.e. anti-Kondo, in sign. This interaction causes spin-disorder broadening of the electron pockets and should be included in models of the normal and superconducting states of Nd$_{0.8}$Sr$_{0.2}$NiO$_2$ The Ni moments differ by 0.2$μ_B$ for FM and antiferromagnetic alignment (the latter are larger), reflecting some itineracy and indicating that Heisenberg coupling of the moments may not provide a quantitative modeling of Ni-Ni exchange coupling. △ Less

Submitted 12 January, 2020; v1 submitted 7 November, 2019; originally announced November 2019.

Comments: 13 pages with Supplemental Material

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

Spin-Orbit Coupling Induced Degeneracy in the Anisotropic Unconventional Superconductor UTe$_2$

Authors: Alexander B. Shick, Warren E. Pickett

Abstract: The orthorhombic uranium dichalcogenide UTe$_2$ displays superconductivity below 1.7 K, with the anomalous feature of retaining 50$\%$ of normal state (ungapped) carriers, according to heat capacity data from two groups. Incoherent transport that crosses over from above 50 K toward a low temperature, Kondo lattice Fermi liquid regime indicates strong magnetic fluctuations and the need to include c… ▽ More The orthorhombic uranium dichalcogenide UTe$_2$ displays superconductivity below 1.7 K, with the anomalous feature of retaining 50$\%$ of normal state (ungapped) carriers, according to heat capacity data from two groups. Incoherent transport that crosses over from above 50 K toward a low temperature, Kondo lattice Fermi liquid regime indicates strong magnetic fluctuations and the need to include correlation effects in theoretical modeling. We report density functional theory plus Hubbard U (DFT+U) results for UTe$_2$ to provide a platform for modeling its unusual behavior, focusing on ferromagnetic (FM, time reversal breaking) long range correlations along the ${\hat a}$ axis as established by magnetization measurements and confirmed by our calculations. States near the Fermi level are dominated by the $j=\frac{5}{2}$ configuration, with the $j_z=\pm\frac{1}{2}$ sectors being effectively degenerate and half-filled. Unlike the small-gap insulating nonmagnetic electronic spectrum, the FM Fermi surfaces are large (strongly metallic) and display low dimensional features, reminiscent of the FM superconductor UGe$_2$. △ Less

Submitted 5 August, 2019; originally announced August 2019.

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

Strong Particle-Hole Symmetry Breaking in a 200 kelvin Superconductor

Authors: Soham S. Ghosh, Yundi Quan, Warren E. Pickett

Abstract: The superconducting state of metals has long provided a classic example of particle-hole symmetry (PHS) at low energy. Fermionic self-energy results based on first principles theory for the electron-phonon coupling in H$_3$S presented here illustrate strong PHS-breaking dynamics arising from the underlying sharp structure in the fermionic density of states. Thus H$_3$S is not only the superconduct… ▽ More The superconducting state of metals has long provided a classic example of particle-hole symmetry (PHS) at low energy. Fermionic self-energy results based on first principles theory for the electron-phonon coupling in H$_3$S presented here illustrate strong PHS-breaking dynamics arising from the underlying sharp structure in the fermionic density of states. Thus H$_3$S is not only the superconductor with the highest critical temperature $T_c$ (through 2018), but its low energy, low temperature properties deviate strongly from textbook behavior. The minor momentum and band dependence of the fermionic self-energy allows evaluation of the momentum-resolved and zone-averaged spectral densities and interacting thermal distribution function, all of which clearly illustrate strong particle-hole asymmetry. △ Less

Submitted 6 June, 2019; originally announced June 2019.

Comments: 7 pages, 5 figures

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

How Compressed Hydrides Produce Room Temperature Superconductivity

Authors: Yundi Quan, Soham S. Ghosh, Warren E. Pickett

Abstract: The 2014-2015 prediction, discovery, and confirmation of record high temperature superconductivity above 200K in H$_3$S, followed by the 2018 extension to superconductivity in the 250-280K range in lanthanum hydride, marks a new era in the longstanding quest for room temperature superconductivity: quest achieved, at the cost of supplying 1.5-2 megabars of pressure. Predictions of numerous high tem… ▽ More The 2014-2015 prediction, discovery, and confirmation of record high temperature superconductivity above 200K in H$_3$S, followed by the 2018 extension to superconductivity in the 250-280K range in lanthanum hydride, marks a new era in the longstanding quest for room temperature superconductivity: quest achieved, at the cost of supplying 1.5-2 megabars of pressure. Predictions of numerous high temperature superconducting metal hydrides $XH_n$ ($X$=metal) have appeared, but are providing limited understanding of what drives the high transition temperature T$_c$, or what limits T$_c$. We apply an opportunistic atomic decomposition of the coupling function to show, first, that the $X$ atom provides coupling strength as commonly calculated, but is it irrelevant for superconductivity; in fact, it is important for analysis that its contribution is neglected. Five $X$H$_n$ compounds, predicted to have T$_c$ in the 150-300K range, are analyzed consistently for their relevant properties, revealing some aspects that confront conventional wisdom. A phonon frequency -- critical temperature ($ω_2$-T$_c$) phase diagram is obtained that reveals a common phase instability limiting T$_c$ at the {\it low pressure} range of each compound. The hydrogen scattering strength is identified and found to differ strongly over the hydrides. A quantity directly proportional to T$_c$ in these hydrides is identified. △ Less

Submitted 6 June, 2019; originally announced June 2019.

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

Superconducting Symmetry Phases and Dominant bands in (Ca-) Intercalated AA- Bilayer Graphene

Authors: Rouhollah Gholami, Rostam Moradian, Sina Moradian, Warren E. Pickett

Abstract: Built on a realistic multiband tight-binding model, mirror symmetry is used to map a calcium-intercalated bilayer graphene Hamiltonian into two independent single layer graphene-like Hamiltonians with renormalized hopping. The quasiparticles exhibit two types of chirality. Here a quasi-particle consists of two electrons from opposing layers where possess an additional quantum number called "cone i… ▽ More Built on a realistic multiband tight-binding model, mirror symmetry is used to map a calcium-intercalated bilayer graphene Hamiltonian into two independent single layer graphene-like Hamiltonians with renormalized hopping. The quasiparticles exhibit two types of chirality. Here a quasi-particle consists of two electrons from opposing layers where possess an additional quantum number called "cone index" which can be regarded as the eigenvalue of mirror symmetry operations. To obtain tight-binding parameters, both effective monolayer Schrodinger equations are solved analytically and fitted to first-principles band structure results. Two quasi-particles (four electrons) can team up to build a Cooper pair with even or odd chirality. Treatment of the pairing Hamiltonian leads to two decoupled gap equations. The pairing of quasi-particles with different cone indexes is forbidden. The decoupled gap equations are solved analytically to obtain all the possible superconducting phases. Two nearly "flat bands" crossing the Fermi energy, each related to the graphene-like structures, are responsible for two distinct superconductivity gaps that emerge. Depending on how much these bands are affected by the intercalant and which is closer to the Fermi energy, distorted s-wave or d-wave superconductivity may become dominant. Numerical calculations reveal that d-wave superconductivity is dominant in both sectors. For these two dominant phases, within the range of 0-6 K which superconductivity has been observed, numerically the transition from single-gap to dual-gap superconductivity is possible. Adopting the two-gap viewpoint of superconductivity in C$_6$CaC$_6$, the dominant $d$-wave states should have the same critical temperature. Around $T_c=2K$ these two relations intersect, otherwise, superconductivity has been realized just in one of these two sectors and disappears in the other one. △ Less

Submitted 31 January, 2020; v1 submitted 4 April, 2019; originally announced April 2019.

Noncentrosymmetric compensated half-metal hosting pure spin Weyl nodes, triple nodal points, nodal loops, and nexus fermions

Authors: Hyo-Sun Jin, Young-Joon Song, W. E. Pickett, K. -W. Lee

Abstract: Materials containing multiple topological characteristics become more exotic when combined with noncentrosymmetric crystal structures and unusual magnetic phases such as the compensated half-metal state, which is gapped in one spin direction and conducting in the other. First principles calculations reveal these multiple topological features in the compensated half-metal Cr$_2$CoAl having neither… ▽ More Materials containing multiple topological characteristics become more exotic when combined with noncentrosymmetric crystal structures and unusual magnetic phases such as the compensated half-metal state, which is gapped in one spin direction and conducting in the other. First principles calculations reveal these multiple topological features in the compensated half-metal Cr$_2$CoAl having neither time-reversal nor inversion symmetries. In the absence of (minor) spin-orbit coupling (SOC), there are (1) a total of twelve pairs of magnetic Weyl points, (2) three distinct sets of triple nodal points near the Fermi level that are (3) interconnected with six symmetry related nodal lines. This combination gives rise to fully spin polarized nexus fermions, in a system with broken time-reversal symmetry but negligible macroscopic magnetic field. The observed high Curie temperature of 750 K and calculated SOC hybridization mixing of several meV should make these nexus fermions readily measurable. Unlike topological features discussed for other Heuslers which emphasize their strong ferromagnetism, this compensated half-metal is impervious to typical magnetic fields, thus providing a complementary set of experimental phenomena. Making use of the soft calculated magnetic state, large magnetic fields can be used to rotate the direction of magnetism, during which certain topological features will evolve. Our results suggest that these features may be common in inverse-Heusler systems, particularly the isostructural and isovalent Ga and In analogs. △ Less

Submitted 1 February, 2019; v1 submitted 13 December, 2018; originally announced December 2018.

Comments: 10 pages, 8 figures (Typos are corrected.)

Journal ref: Phys. Rev. Materials 3, 021201 (2019)

Pressure-induced Frustration of Magnetic Coupling in Elemental Europium

Authors: Shu-Ting Pi, Sergey Y. Savrasov, Warren E. Pickett

Abstract: Applying linear response and the magnetic force theorem in correlated density functional theory, the inter-sublattice exchange constants of antiferromagnetic Eu are calculated and found to vanish near the pressure of P$_c$=82 GPa, just where magnetic order is observed experimentally to be lost. The Eu $4f^7$ moment remains unchanged at high pressure, again in agreement with spectroscopic measureme… ▽ More Applying linear response and the magnetic force theorem in correlated density functional theory, the inter-sublattice exchange constants of antiferromagnetic Eu are calculated and found to vanish near the pressure of P$_c$=82 GPa, just where magnetic order is observed experimentally to be lost. The Eu $4f^7$ moment remains unchanged at high pressure, again in agreement with spectroscopic measurements, leaving the picture of perfect frustration of interatomic Ruderman-Kittel-Kasuya-Yoshida couplings in a broad metallic background, leaving a state of electrons strongly exchange coupled to arbitrarily oriented, possibly quasistatic local moments. This strongly frustrated state gives way to superconductivity at T$_c$=1.7K, observed experimentally. These phenomena, and free energy considerations related to correlations, suggest an unusual phase of matter that is discussed within the scenarios of the Doniach Kondo lattice phase diagram, the metallic spin glass class, and itinerant spin liquid or spin gas systems. △ Less

Submitted 14 August, 2018; v1 submitted 5 June, 2018; originally announced June 2018.

Comments: 6 pages, 4 figures

Journal ref: Phys. Rev. Lett. 122, 057201 (2019)

Superconducting Phases in Lithium Decorated Graphene LiC 6

Authors: Rouhollah Gholami, Rostam Moradian, Sina Moradian, Warren E. Pickett

Abstract: A study of possible superconducting phases of graphene has been constructed in detail. A realistic tight binding model, fit to ab initio calculations, accounts for the Li-decoration of graphene with broken lattice symmetry, and includes $s$ and $d$ symmetry Bloch character that influences the gap symmetries that can arise. The resulting seven hybridized Li-C orbitals that support nine possible bon… ▽ More A study of possible superconducting phases of graphene has been constructed in detail. A realistic tight binding model, fit to ab initio calculations, accounts for the Li-decoration of graphene with broken lattice symmetry, and includes $s$ and $d$ symmetry Bloch character that influences the gap symmetries that can arise. The resulting seven hybridized Li-C orbitals that support nine possible bond pairing amplitudes. The gap equation is solved for all possible gap symmetries. One band is weakly dispersive near the Fermi energy along $Γ\rightarrow M$ where its Bloch wave function has linear combination of $d_{x^2-y^2}$ and $d_{xy}$ character, and is responsible for $ d_{x^{2}-y^{2}}$ and $d_{xy}$ pairing with lowest pairing energy in our model. These symmetries almost preserve properties from a two band model of pristine graphene. Another part of this band, along $K\rightarrow Γ$, is nearly degenerate with upper $s$ band that favors extended $s$ wave pairing which is not found in two band model. Upon electron doping to a critical chemical potential $μ_1=0.22 eV$ the pairing potential decreases, then increases until a second critical value $μ_2$=1.3 eV at which a phase transition to a new phase which is not appear in two band model. This phase in the pristine graphene converts to usual extended s-wave pairing. △ Less

Submitted 2 September, 2018; v1 submitted 10 April, 2018; originally announced April 2018.

Journal ref: Scientific Reports 8, 13795 (2018)

Coexistence of Triple Nodal Points, Nodal Links, and Unusual Flat Bands in intermetallic ${\cal A}$Pd$_3$ (${\cal A}$=Pb, Sn)

Authors: Kyo-Hoon Ahn, W. E. Pickett, K. -W. Lee

Abstract: We investigate the electronic structure and several properties, and topological character, of the cubic time-reversal invariant intermetallic compounds PbPd$_3$ and SnPd$_3$ using density functional theory based methods. These compounds have a dispersionless band along the $Γ-X$ line, forming the top of the Pd $4d$ bands and lying within a few meV of the Fermi level $E_F$. Effects of the flat band… ▽ More We investigate the electronic structure and several properties, and topological character, of the cubic time-reversal invariant intermetallic compounds PbPd$_3$ and SnPd$_3$ using density functional theory based methods. These compounds have a dispersionless band along the $Γ-X$ line, forming the top of the Pd $4d$ bands and lying within a few meV of the Fermi level $E_F$. Effects of the flat band on transport and optical properties have been inspected by varying the doping concentration treated with the virtual crystal approximation for substitution on the Pb site. In the absence of spin-orbit coupling (SOC), we find triple nodal points and three-dimensional nodal loops, which are known to lead to surface bands and drumhead states, respectively, which we discuss for PbPd$_3$. SOC removes degeneracy in most of the zone, providing a topological index $Z_2$=1 on the $k_z=0$ plane that indicates a topological character on that plane. The isovalent and isostructural compound SnPd$_3$ shows only minor differences in its electronic structures, so it is expected to display similar electronic, transport, and topological properties. △ Less

Submitted 15 July, 2018; v1 submitted 21 March, 2018; originally announced March 2018.

Comments: 13 pages

Journal ref: Phys. Rev. B 98, 035130 (2018)

Design of Chern Insulating Phases in Honeycomb Lattices

Authors: W. E. Pickett, K. -W. Lee, R. Pentcheva

Abstract: The search for robust examples of the magnetic version of topological insulators, referred to as quantum anomalous Hall insulators or simply Chern insulators, so far lacks success. Our groups have explored two distinct possibilities based on multiorbital 3d oxide honeycomb lattices. Each has a Chern insulating phase near the ground state, but materials parameters were not appropriate to produce a… ▽ More The search for robust examples of the magnetic version of topological insulators, referred to as quantum anomalous Hall insulators or simply Chern insulators, so far lacks success. Our groups have explored two distinct possibilities based on multiorbital 3d oxide honeycomb lattices. Each has a Chern insulating phase near the ground state, but materials parameters were not appropriate to produce a viable Chern insulator. Further exploration of one of these classes, by substituting open shell 3d with 4d and 5d counterparts, has led to realistic prediction of Chern insulating ground states. Here we recount the design process, discussing the many energy scales that are active in participating (or resisting) the desired Chern insulator phase. △ Less

Submitted 20 February, 2018; originally announced February 2018.

Comments: 4 pages

Journal ref: Phyica C 549, 99 (2018)

Atomic-layer-resolved composition and electronic structure of the cuprate Bi$_2$Sr$_2$CaCu$_2$O$_{8+δ}$ from soft x-ray standing-wave photoemission

Authors: Cheng-Tai Kuo, Shih-Chieh Lin, Giuseppina Conti, Shu-Ting Pi, Luca Moreschini, Aaron Bostwick, Julia Meyer-Ilse, Eric Gullikson, Jeffrey B. Kortright, Slavomir Nemšák, Julien E. Rault, Patrick Le Fèvre, François Bertran, Andrés F. Santander-Syro, Ivan A. Vartanyants, Warren E. Pickett, Romuald Saint-Martin, Amina Taleb, Charles S. Fadley

Abstract: A major remaining challenge in the superconducting cuprates is the unambiguous differentiation of the composition and electronic structure of the CuO$_2$ layers and those of the intermediate layers. The large c axis for these materials permits employing soft x-ray (930.3 eV) standing wave (SW) excitation in photoemission that yields atomic layer-by-atomic layer depth resolution of these properties… ▽ More A major remaining challenge in the superconducting cuprates is the unambiguous differentiation of the composition and electronic structure of the CuO$_2$ layers and those of the intermediate layers. The large c axis for these materials permits employing soft x-ray (930.3 eV) standing wave (SW) excitation in photoemission that yields atomic layer-by-atomic layer depth resolution of these properties. Applying SW photoemission to Bi$_2$Sr$_2$CaCu$_2$O$_{8+δ}$ yields the depth distribution of atomic composition and the layer-resolved densities of states. We detect significant Ca presence in the SrO layers and oxygen bonding to three different cations. The layer-resolved valence electronic structure is found to be strongly influenced by the supermodulation structure--as determined by comparison to DFT calculations, by Ca-Sr intermixing, and by the Cu 3d-3d Coulomb interaction, further clarifying the complex interactions in this prototypical cuprate. Measurements of this type for other quasi-two-dimensional materials with large-c represent a promising future direction. △ Less

Submitted 24 September, 2018; v1 submitted 16 January, 2018; originally announced January 2018.

Comments: 22 pages, 5 figures, plus Supplemental Material (15 pages)

Report number: arXiv:1801.05142

Journal ref: Phys. Rev. B 98, 155133 (2018)

Room Temperature Superconductivity Revolution: Foreshadowed by Victorians, Enabled by Millenials

Authors: Warren E. Pickett

Abstract: Room temperature superconductivity has been the most prominent, highly ambitious, but still imaginable, acme of materials physics for half a century. The struggle toward this revolution was foreshadowed by a Victorian novelist and championed,unsuccessfully, by dogged physicists in the 1960s to 1980s who had a workable theory but uncompliant materials. Discovery of superconductivity of H$_3$S at 20… ▽ More Room temperature superconductivity has been the most prominent, highly ambitious, but still imaginable, acme of materials physics for half a century. The struggle toward this revolution was foreshadowed by a Victorian novelist and championed,unsuccessfully, by dogged physicists in the 1960s to 1980s who had a workable theory but uncompliant materials. Discovery of superconductivity of H$_3$S at 200 K in the 160-200 GPa pressure range has renewed anticipation of yet higher values of the critical temperature T$_c$. With the several reports of metalization of hydrogen, and theoretical extensions enabled by modern algorithms and unprecedented computational hardware and spurred forward by the Materials Genome Initiative, it is possible that the room temperature precipice has thereby already been breached in a silent revolution. This concise note draws analogies of this development with an earlier revolution. △ Less

Submitted 30 December, 2017; originally announced January 2018.

Comments: 4 pages, 3 figures

Survey of the class of isovalent antiperovskite alkaline earth-pnictide compounds

Authors: Wen Fong Goh, Warren E. Pickett

Abstract: The few reported members of the antiperovskite structure class $Ae_3Pn_APn_B$ of alkaline earth ($Ae$ = Ca,Sr,Ba) pnictides ($Pn$ = N,P,As,Sb,Bi) compounds are all based on the B-site anion $Pn_B$=N. All fit can be categorized as narrow gap semiconductors, making them of interest for several reasons. Because chemical reasoning suggests that more members of this class may be stable, we provide here… ▽ More The few reported members of the antiperovskite structure class $Ae_3Pn_APn_B$ of alkaline earth ($Ae$ = Ca,Sr,Ba) pnictides ($Pn$ = N,P,As,Sb,Bi) compounds are all based on the B-site anion $Pn_B$=N. All fit can be categorized as narrow gap semiconductors, making them of interest for several reasons. Because chemical reasoning suggests that more members of this class may be stable, we provide here a density functional theory (DFT) based survey of this entire class of $3\times5\times5$ compounds. We determine first the relative energetic stability of the distribution of pairs of $Pn$ ions in the A and B sites of the structure, finding that the $B$ site always favors the small pnictogen anion. The trends of the calculated energy gaps with $Ae$ cation and $Pn$ anions are determined, and we study effects of spin-orbit coupling as well as two types of gap corrections to the conventional DFT electronic spectrum. Because there have been suggestions that this class harbors topological insulating phases, we have given this possibility attention and found that energy gap corrections indicate the cubic structures will provide at most a few topological insulators. Structural instability is addressed by calculating phonon dispersion curves for a few compounds, with one outcome being that distorted structures should be investigated further for thermoelectric and topological character. Examples of the interplay between spin-orbit coupling and strain on the topological nature are provided. A case study of $\mathrm{Ca_3BiP}$ including the effect of strain illustrates how a topological semimetal can be transformed into topological insulator and Dirac semimetal. △ Less

Submitted 29 October, 2017; originally announced October 2017.

Journal ref: Phys. Rev. B 97, 035202 (2018)

Local nematic susceptibility in stressed BaFe$_2$As$_2$ from NMR electric field gradient measurements

Authors: T. Kissikov, R. Sarkar, M. Lawson, B. T. Bush, E. I. Timmons, M. A. Tanatar, R. Prozorov, S. L. Bud'ko, P. C. Canfield, R. M. Fernandes, W. F. Goh, W. E. Pickett, N. J. Curro

Abstract: The electric field gradient (EFG) tensor at the $^{75}$As site couples to the orbital occupations of the As p-orbitals and is a sensitive probe of local nematicity in BaFe$_2$As$_2$. We use nuclear magnetic resonance to measure the nuclear quadrupolar splittings and find that the EFG asymmetry responds linearly to the presence of a strain field in the paramagnetic phase. We extract the nematic sus… ▽ More The electric field gradient (EFG) tensor at the $^{75}$As site couples to the orbital occupations of the As p-orbitals and is a sensitive probe of local nematicity in BaFe$_2$As$_2$. We use nuclear magnetic resonance to measure the nuclear quadrupolar splittings and find that the EFG asymmetry responds linearly to the presence of a strain field in the paramagnetic phase. We extract the nematic susceptibility from the slope of this linear response as a function of temperature and find that it diverges near the structural transition in agreement with other measures of the bulk nematic susceptibility. Our work establishes an alternative method to extract the nematic susceptibility which, in contrast to transport methods, can be extended inside the superconducting state. △ Less

Submitted 26 October, 2017; originally announced October 2017.

Comments: 6 pages, 4 figures

Journal ref: Phys. Rev. B 96, 241108 (2017)

Perovskite ThTaN3: a Large Thermopower Topological Crystalline Insulator

Authors: Myung-Chul Jung, K. -W. Lee, W. E. Pickett

Abstract: ThTaN$_3$, a rare cubic perovskite nitride semiconductor, has been studied using {\it ab initio} methods. Spin-orbit coupling (SOC) results in band inversion and a band gap of 150 meV at the zone center. In spite of the trivial $Z_2$ indices, two pairs of spin-polarized surface bands cross the gap near the zone center, indicating that this system is a topological crystalline insulator with the mir… ▽ More ThTaN$_3$, a rare cubic perovskite nitride semiconductor, has been studied using {\it ab initio} methods. Spin-orbit coupling (SOC) results in band inversion and a band gap of 150 meV at the zone center. In spite of the trivial $Z_2$ indices, two pairs of spin-polarized surface bands cross the gap near the zone center, indicating that this system is a topological crystalline insulator with the mirror Chern number of $|{\cal C}_m|=2$ protected by the mirror and $C_4$ rotational symmetries. Additionally, SOC doubles the Seebeck coefficient, leading to a maximum of $\sim$400 $μ$V/K at 150 K for carrier-doping levels of several $10^{17}$/cm$^3$. ThTaN$_3$ combines excellent bulk thermopower with parallel conduction through topological surface states that provide a platform for large engineering devices with ever larger figures of merit. △ Less

Submitted 15 March, 2018; v1 submitted 4 September, 2017; originally announced September 2017.

Comments: 8 pages, 7 figures

Journal ref: Phys. Rev. B 97, 121104 (2018)

All 3d electron-hole bilayers in CrN/MgO(111) multilayers

Authors: Antia S. Botana, Victor Pardo, Warren E. Pickett

Abstract: CrN/MgO(111) multilayers modeled via \textit{ab initio} calculations give rise to nanoscale, scalable, spatially separated two-dimensional electron and hole gases (2DEG+2DHG), each confined to its own CrN interface. Due to the Cr $3d^3$ configuration, both electron and hole gases are based on correlated transition metal layers involving bands of $3d$ character. Transport calculations predict each… ▽ More CrN/MgO(111) multilayers modeled via \textit{ab initio} calculations give rise to nanoscale, scalable, spatially separated two-dimensional electron and hole gases (2DEG+2DHG), each confined to its own CrN interface. Due to the Cr $3d^3$ configuration, both electron and hole gases are based on correlated transition metal layers involving bands of $3d$ character. Transport calculations predict each subsystem will have a large thermopower, on the order of 250 $μ$$V /K$ at room temperature. These heterostructures combine a large thermoelectric efficiency with scalable nanoscale conducting sheets; for example, operating at a temperature difference of 50K, 40 bilayers could produce a 1 V voltage with a film thickness of 100 nm. △ Less

Submitted 12 May, 2017; originally announced May 2017.

Comments: 8 pages, 6 figures

Journal ref: Phys. Rev. Applied 7, 024002 (2017)

A single nodal loop of accidental degeneracies in minimal symmetry: triclinic CaAs$_3$

Authors: Y. Quan, Z. P. Yin, W. E. Pickett

Abstract: The existence of closed loops of degeneracies in crystals has been intimately connected to associated crystal symmetries, raising the question: what is the minimum symmetry required for topological character, and can one find an example? Triclinic CaAs$_3$, in space group $P{\bar 1}$ with only a center of inversion, has been found to display, without need for tuning, a nodal loop of accidental deg… ▽ More The existence of closed loops of degeneracies in crystals has been intimately connected to associated crystal symmetries, raising the question: what is the minimum symmetry required for topological character, and can one find an example? Triclinic CaAs$_3$, in space group $P{\bar 1}$ with only a center of inversion, has been found to display, without need for tuning, a nodal loop of accidental degeneracies with topological character, centered on one face of the Brillouin zone that is otherwise fully gapped. The small loop is very flat in energy, yet is cut four times by the Fermi energy, a condition that results in an intricate repeated touching of inversion related pairs of Fermi surfaces at Weyl points. Spin-orbit coupling lifts the fourfold degeneracy along the loop, leaving trivial Kramers pairs. With its single nodal loop that emerges without protection from any point group symmetry, CaAs$_3$ represents the primal "hydrogen atom" of nodal loop systems. △ Less

Submitted 13 March, 2017; originally announced March 2017.

Magnetic Order-Disorder Transitions on a 1/3 - Depleted Square Lattice

Authors: H. -M. Guo, T. Mendes-Santos, W. E. Pickett, R. T. Scalettar

Abstract: Quantum Monte Carlo simulations are used to study the magnetic and transport properties of the Hubbard Model, and its strong coupling Heisenberg limit, on a one-third depleted square lattice. This is the geometry occupied, after charge ordering, by the spin-$\frac{1}{2}$ Ni$^{1+}$ atoms in a single layer of the nickelate materials La$_4$Ni$_3$O$_8$ and (predicted) La$_3$Ni$_2$O$_6$. Our model is a… ▽ More Quantum Monte Carlo simulations are used to study the magnetic and transport properties of the Hubbard Model, and its strong coupling Heisenberg limit, on a one-third depleted square lattice. This is the geometry occupied, after charge ordering, by the spin-$\frac{1}{2}$ Ni$^{1+}$ atoms in a single layer of the nickelate materials La$_4$Ni$_3$O$_8$ and (predicted) La$_3$Ni$_2$O$_6$. Our model is also a description of strained graphene, where a honeycomb lattice has bond strengths which are inequivalent. For the Heisenberg case, we determine the location of the quantum critical point (QCP) where there is an onset of long range antiferromagnetic order (LRAFO), and the magnitude of the order parameter, and then compare with results of spin wave theory. An ordered phase also exists when electrons are itinerant. In this case, the growth in the antiferromagnetic structure factor coincides with the transition from band insulator to metal in the absence of interactions. △ Less

Submitted 20 October, 2016; originally announced October 2016.

Journal ref: Phys. Rev. B 95, 045131 (2017)

Tuning ferromagnetic BaFe$_2$(PO$_4$)$_2$ through a high Chern number topological phase

Authors: Young-Joon Song, Kyo-Hoon Ahn, W. E. Pickett, K. -W. Lee

Abstract: There is strong interest in discovering or designing wide gap Chern insulators. Here we follow a Chern insulator to trivial Mott insulator transition versus interaction strength U in a honeycomb lattice Fe-based transition metal oxide, discovering that a spin-orbit coupling energy scale $ξ$=40 meV can produce and maintain a topologically entangled Chern insulating state against large band structur… ▽ More There is strong interest in discovering or designing wide gap Chern insulators. Here we follow a Chern insulator to trivial Mott insulator transition versus interaction strength U in a honeycomb lattice Fe-based transition metal oxide, discovering that a spin-orbit coupling energy scale $ξ$=40 meV can produce and maintain a topologically entangled Chern insulating state against large band structure changes arising from an interaction strength U up to 60 times as large. Within the Chern phase the minimum gap switches from the zone corner K to the zone center $Γ$ while maintaining the topological structure. At a critical strength $U_c$, the continuous evolution of the electronic structure encounters a gap closing then reopening, upon which the system reverts to a trivial Mott insulating phase. This Chern insulator phase of honeycomb lattice Fe$^{2+}$ BaFe$_2$(PO$_4$)$_2$ corresponds to a large Chern number C=-3 that will provide enhanced anomalous Hall conductivity due to the associated three edge states threading through the bulk gap of 80 meV. △ Less

Submitted 22 September, 2016; v1 submitted 31 August, 2016; originally announced September 2016.

Comments: 8 pages, 7 figures

Journal ref: Phys. Rev. B 94, 125134 (2016)

Role of Hydrogen in the Electronic Properties of CaFeAsH-based Superconductors

Authors: Y. N. Huang, D. Y. Liu, L. J. Zou, W. E. Pickett

Abstract: The electronic and magnetic properties of the new hydride superconductor CaFeAsH, which superconducts up to 47 K when electron-doped with La, and the isovalent alloy system CaFeAsH$_{1-x}$F$_x$ are investigated using density functional based methods. The $\vec Q = (π,π,0)$ peak of the nesting function $ξ(\vec q)$ is found to be extremely strong and sharp, and additional structure in $ξ(\vec q)$ as… ▽ More The electronic and magnetic properties of the new hydride superconductor CaFeAsH, which superconducts up to 47 K when electron-doped with La, and the isovalent alloy system CaFeAsH$_{1-x}$F$_x$ are investigated using density functional based methods. The $\vec Q = (π,π,0)$ peak of the nesting function $ξ(\vec q)$ is found to be extremely strong and sharp, and additional structure in $ξ(\vec q)$ associated with the near-circular Fermi surfaces (FSs) that may impact low energy excitations is quantified. The unusual band introduced by H, which shows strong dispersion perpendicular to the FeAs layers, is shown to be connected to a peculiar van Hove singularity just below the Fermi level. This band provides a three dimensional electron ellipsoid Fermi surface not present in other Fe-based superconducting materials nor in CaFeAsF. Electron doping by 25\% La or Co has minor affect on this ellipsoid Fermi surface, but suppresses FS nesting strongly, consistent with the viewpoint that eliminating strong nesting and the associated magnetic order allows high T$_c$ superconductivity to emerge. Various aspects of the isovalent alloy system CaFeAsH$_{1-x}$F$_x$ and means of electron doping are discussed in terms of influence of incipient bands. △ Less

Submitted 8 June, 2016; originally announced June 2016.

Comments: 9 pages, 8 figures

Journal ref: Physical Review B 93, 195148 (2016)

Charge ordering in Ni$^{1+}$/Ni$^{2+}$ nickelates: La$_4$Ni$_3$O$_8$ and La$_3$Ni$_2$O$_6$

Authors: A. S. Botana, V. Pardo, W. E. Pickett, M. R. Norman

Abstract: Ab initio calculations have been performed to unravel the origin of the recently found superlattice peaks in the trilayer nickelate La$_4$Ni$_3$O$_8$. These peaks arise from static charge ordering of Ni$^{2+}$/ Ni$^{1+}$ stripes oriented at 45$^{\circ}$ to the Ni-O bonds. An insulating state originates from a combination of structural distortions and magnetic order, with the gap being formed solel… ▽ More Ab initio calculations have been performed to unravel the origin of the recently found superlattice peaks in the trilayer nickelate La$_4$Ni$_3$O$_8$. These peaks arise from static charge ordering of Ni$^{2+}$/ Ni$^{1+}$ stripes oriented at 45$^{\circ}$ to the Ni-O bonds. An insulating state originates from a combination of structural distortions and magnetic order, with the gap being formed solely within the d$_{x^2-y^2}$ manifold of states. When doped, electrons or holes would go into these states, in a similar fashion to what occurs in the cuprates. Analogous calculations suggest that checkerboard charge order should occur in the bilayer nickelate La$_3$Ni$_2$O$_6$. These results reveal a close connection between La$_4$Ni$_3$O$_8$ and La$_3$Ni$_2$O$_6$ with La$_{2-x}$Sr$_x$NiO$_4$ for x=1/3 and x=1/2, respectively. △ Less

Submitted 30 August, 2016; v1 submitted 21 April, 2016; originally announced April 2016.

Journal ref: Phys. Rev. B 94, 081105(R) (2016)

Design of Mott and topological phases on buckled 3d-oxide honeycomb lattices

Authors: David Doennig, Santu Baidya, Warren E. Pickett, Rossitza Pentcheva

Abstract: Perovskite bilayers with (111)-orientation combine a honeycomb lattice as a key feature with the strongly correlated, multiorbital nature of electrons in transition metal oxides. In a systematic DFT+$U$ study of (111)-oriented (La$X$O$_3$)$_2$/(LaAlO$_3$)$_4$ superlattices, we establish trends in the evolution of ground states versus band filling in (111)-oriented (La$X$O$_3$)$_2$/(LaAlO$_3$)$_4$… ▽ More Perovskite bilayers with (111)-orientation combine a honeycomb lattice as a key feature with the strongly correlated, multiorbital nature of electrons in transition metal oxides. In a systematic DFT+$U$ study of (111)-oriented (La$X$O$_3$)$_2$/(LaAlO$_3$)$_4$ superlattices, we establish trends in the evolution of ground states versus band filling in (111)-oriented (La$X$O$_3$)$_2$/(LaAlO$_3$)$_4$ superlattices, with $X$ spanning the entire $3d$ transition metal series. The competition between local quasi-cubic and global triangular symmetry triggers unanticipated broken symmetry phases, with mechanisms ranging from Jahn-Teller distortions, to charge-, spin-, and orbital-ordering. LaMnO$_3$, where spin-orbit coupling opens a sizable gap in the Dirac-point Fermi surface, emerges as a topological Chern insulator. △ Less

Submitted 30 October, 2015; originally announced October 2015.

Comments: 6 pages, 3 figures

Journal ref: Phys. Rev. B 93, 165145 (2016)

Interplay between Spin-Orbit Coupling and Strong Correlation Effects: Comparison of Three Osmate Double Perovskites: Ba$_{2}A$OsO$_{6}$ ($A=$Na, Ca, Y)

Authors: Shruba Gangopadhyay, Warren E. Pickett

Abstract: High formal valence Os-based double perovskites are a focus of current interest because they display strong interplay of large spin orbit coupling and strong electronic correlation. Here we present the electronic and magnetic characteristics of a sequence of three cubic Os based double perovskites Ba$_{2}$$A$OsO$_{6}$, ($A=$Na, Ca, Y), with formal valences (charge states) of Os$^{+7}~(d^{1})$, Os… ▽ More High formal valence Os-based double perovskites are a focus of current interest because they display strong interplay of large spin orbit coupling and strong electronic correlation. Here we present the electronic and magnetic characteristics of a sequence of three cubic Os based double perovskites Ba$_{2}$$A$OsO$_{6}$, ($A=$Na, Ca, Y), with formal valences (charge states) of Os$^{+7}~(d^{1})$, Os$^{+6}(d^{2})$, Os$^{+5}(d^{3})$. For these first principles based calculations we apply an ``exact exchange for correlated electrons'' functional, with exact exchange applied in a hybrid fashion solely to the Os($5d$) states. While Ba$_{2}$NaOsO$_{6}$ is a ferromagnetic Dirac-Mott insulator studied previously, the other two show antiferromagnetic ordering while all retain the cubic (undistorted) structure. For comparison purposes we have investigated only the ferromagnetic ordered phase. A metal-insulator transition is predicted to occur upon rotating the direction of magnetization in all three materials, reflecting the central role of spin orbit coupling in these small gap osmates. Surprises arising from comparing formal charge states with the radial charge densities are discussed. Chemical shielding factors and orbital susceptibilities are provided for comparison with future nuclear magnetic resonance data. △ Less

Submitted 6 October, 2015; originally announced October 2015.

Nitride Multilayers as a Platform for Parallel Two-Dimensional Electron-Hole Gases: MgO/ScN(111)

Authors: Antia S. Botana, Victor Pardo, Warren E. Pickett

Abstract: At interfaces between insulating oxides LaAlO$_3$ and SrTiO$_3$, a two dimensional electron gas (2DEG) has been observed and well studied, while the predicted hole gas (2DHG) has not been realized due to the strong tendency of holes in oxygen $2p$ orbitals to localize. Here we propose, via ab initio calculations, an unexplored class of materials for the realization of parallel two dimensional (2D)… ▽ More At interfaces between insulating oxides LaAlO$_3$ and SrTiO$_3$, a two dimensional electron gas (2DEG) has been observed and well studied, while the predicted hole gas (2DHG) has not been realized due to the strong tendency of holes in oxygen $2p$ orbitals to localize. Here we propose, via ab initio calculations, an unexplored class of materials for the realization of parallel two dimensional (2D), two carrier (electron+hole) gases: nitride-oxide heterostructures, with (111)-oriented ScN and MgO as the specific example. Beyond a critical thickness of five ScN layers, this interface hosts spatially separated conducting Sc-$3d$ electrons and N-$2p$ holes, each confined to $\sim$two atomic layers -- the transition metal nitride provides both gases. A guiding concept is that the N$^{3-}$ anion should promote robust two carrier 2D hole conduction compared to that of O$^{2-}$; metal mononitrides are mostly metallic and even superconducting while most metal monoxides are insulating. A second positive factor is that the density of transition metal ions, hence of a resulting 2DG, is about three times larger for a rocksalt (111) interface than for a perovskite(001) interface. Our results, including calculation of Hall coefficient and thermopower for each conducting layer separately, provide guidance for new exploration, both experimental and theoretical, on nitride-based conducting gases that should promote study of long sought exotic states viz. new excitonic phases and distinct, nanoscale parallel superconducting nanolayers. △ Less

Submitted 25 January, 2016; v1 submitted 28 September, 2015; originally announced September 2015.

Comments: 6 pages, 4 figures

Journal ref: Phys. Rev. B 93, 085125 (2016)

Disturbing the dimers: electron- and hole-doping in the intermetallic insulator FeGa$_3$

Authors: Antia S. Botana, Yundi Quan, Warren E. Pickett

Abstract: Insulating FeGa$_3$ poses peculiar puzzles beyond the occurrence of an electronic gap in an intermetallic compound. This Fe-based material has a very distinctive structural characteristic with the Fe atoms occurring in dimers. The insulating gap can be described comparably well in either the weakly correlated limit or the strongly correlated limit within density functional theory viewpoints, where… ▽ More Insulating FeGa$_3$ poses peculiar puzzles beyond the occurrence of an electronic gap in an intermetallic compound. This Fe-based material has a very distinctive structural characteristic with the Fe atoms occurring in dimers. The insulating gap can be described comparably well in either the weakly correlated limit or the strongly correlated limit within density functional theory viewpoints, where the latter corresponds to singlet formation on the Fe$_2$ dimers. Though most of the calculated occupied Wannier functions are an admixture of Fe $3d$ and Ga $4s$ or $4p$ states, there is a single bonding-type Wannier function per spin centered on each Fe$_2$ dimer. Density functional theory methods have been applied to follow the evolution of the magnetic properties and electronic spectrum with doping, where unusual behavior is observed experimentally. Both electron and hole doping are considered, by Ge and Zn on the Ga site, and by Co and Mn on the Fe site, the latter introducing direct disturbance of the Fe$_2$ dimer. Results from weakly and strongly correlated pictures are compared. Regardless of the method, magnetism including itinerant phases appears readily with doping. The correlated picture suggests that in the low doping limit Mn (for Fe) produces an in-gap hole state, while Co (for Fe) introduces a localized electronic gap state. △ Less

Submitted 30 November, 2015; v1 submitted 28 September, 2015; originally announced September 2015.

Comments: 11 pages, 9 figures

Journal ref: Phys. Rev. B 92, 155134 (2015)

van Hove singularities and spectral smearing in high temperature superconducting H3S

Authors: Yundi Quan, Warren E. Pickett

Abstract: The superconducting phase of hydrogen sulfide at Tc=200 K observed by Drozdov and collaborators at pressures around 200 GPa is simple bcc Im-3m H3S, predicted beforehand by Duan {\it et al.}, has experimental confirmation. The various "extremes" that are involved -- pressure, implying extreme reduction of volume, extremely high H phonon energy scale around 1500K, extremely high temperature for a s… ▽ More The superconducting phase of hydrogen sulfide at Tc=200 K observed by Drozdov and collaborators at pressures around 200 GPa is simple bcc Im-3m H3S, predicted beforehand by Duan {\it et al.}, has experimental confirmation. The various "extremes" that are involved -- pressure, implying extreme reduction of volume, extremely high H phonon energy scale around 1500K, extremely high temperature for a superconductor -- necessitates a close look at new issues raised by these characteristics in relation to high Tc. We use first principles methods to analyze the H3S electronic structure, particularly the van Hove singularities (vHs) and the effect of sulfur. Focusing on the two closely spaced vHs near the Fermi level that give rise to the impressively sharp peak in the density of states, the implications of strong coupling Migdal-Eliashberg theory are assessed. The electron spectral density smearing due to virtual phonon emission and absorption needs to be included explicitly to obtain accurate theoretical predictions and current understanding. Means for increasing Tc in H3S-like materials are addressed. △ Less

Submitted 3 September, 2015; v1 submitted 18 August, 2015; originally announced August 2015.

Comments: 13 pages, 10 figures

Journal ref: Phys. Rev. B 93, 104526 (2016)

Spin-orbit interaction driven collective electron-hole excitations in a noncentrosymmetric nodal loop Weyl semimetal

Authors: Kyo-Hoon Ahn, K. -W. Lee, W. E. Pickett

Abstract: NbP is one member of a new class of nodal loop semimetals characterized by the cooperative effects of spin-orbit coupling (SOC) and a lack of inversion center. Here transport and spectroscopic properties of NbP are evaluated using density functional theory methods. SOC together with the lack of inversion symmetry splits degeneracies, giving rise to "Russian doll nested" Fermi surfaces containing 4… ▽ More NbP is one member of a new class of nodal loop semimetals characterized by the cooperative effects of spin-orbit coupling (SOC) and a lack of inversion center. Here transport and spectroscopic properties of NbP are evaluated using density functional theory methods. SOC together with the lack of inversion symmetry splits degeneracies, giving rise to "Russian doll nested" Fermi surfaces containing 4*10$^{-4}$ electron (hole) carriers/f.u. Due to the modest SOC strength in Nb, the Fermi surfaces map out the Weyl nodal loops. Calculated structure around T$^*$~100 K in transport properties reproduces well the observed transport behavior only when SOC is included, attesting to the precision of the (delicate) calculations and the stoichiometry of the samples. Low energy collective electron-hole excitations (plasmons) in the 20-60 meV range result from the nodal loop splitting. △ Less

Submitted 22 September, 2016; v1 submitted 13 July, 2015; originally announced July 2015.

Comments: 6 pages (corrected Fig. 6)

Journal ref: Phys. Rev. B 92, 115149 (2015)

Large orbital moment and spin-orbit enabled Mott transition in the Ising Fe honeycomb lattice BaFe2(PO4)2

Authors: Young-Joon Song, K. -W. Lee, W. E. Pickett

Abstract: BaFe2(PO4)2 is an unusual Ising insulating ferromagnet based on the Fe$^{2+}$ spin $S$ = 2 ion, the susceptibility of which suggests a large orbital component to the Fe local moment. We apply density functional theory based methods to obtain a microscopic picture of the competing interactions and the critical role of spin-orbit coupling (SOC) in this honeycomb lattice system. The low-temperature f… ▽ More BaFe2(PO4)2 is an unusual Ising insulating ferromagnet based on the Fe$^{2+}$ spin $S$ = 2 ion, the susceptibility of which suggests a large orbital component to the Fe local moment. We apply density functional theory based methods to obtain a microscopic picture of the competing interactions and the critical role of spin-orbit coupling (SOC) in this honeycomb lattice system. The low-temperature ferromagnetic phase displays a half-semimetallic Dirac point pinning the Fermi level and preventing gap opening before consideration of SOC, presenting a case in which correlation effects modeled by a repulsive Hubbard $U$ fail to open a gap. Simultaneous inclusion of both correlation and SOC drives a large orbital moment in excess of 0.7 $μ_B$ (essentially $L$ = 1) for spin aligned along the $\hat{c}$ axis, with a gap comparable with the inferred experimental value. The large orbital moment accounts for the large Ising anisotropy, in spite of the small magnitude of the SOC strength on the 3$d$ (Fe) ion. Ultimately, the Mott-Hubbard gap is enabled by degeneracy lifting by SOC and the large Fe moments, rather than by standard Hubbard interactions alone. We suggest that competing orbital occupations are responsible for the structural transitions involved in the observed re-entrant rhombohedral-triclinic-rhombohedral sequence. △ Less

Submitted 6 September, 2015; v1 submitted 3 June, 2015; originally announced June 2015.

Comments: 6+ pages, 5 figures

Journal ref: Phys. Rev. B 92, 125109 (2015)

Organometallic-like localization of 4d-derived spins in an inorganic conducting niobium suboxide

Authors: Kwan-Woo Lee, Warren E. Pickett

Abstract: Based on the refined crystal structure comprised of columns of 3*4 planar blocks of NbO6 octahedra and first principles electronic structure methods, we find that orthorhombic (o)-Nb$_{12}$O$_{29}$ introduces a new class of transition metal oxide. The electronic system consists of a large Nb dimer-based localized orbital comparable in size to those in organometallic compounds, yet is tightly bound… ▽ More Based on the refined crystal structure comprised of columns of 3*4 planar blocks of NbO6 octahedra and first principles electronic structure methods, we find that orthorhombic (o)-Nb$_{12}$O$_{29}$ introduces a new class of transition metal oxide. The electronic system consists of a large Nb dimer-based localized orbital comparable in size to those in organometallic compounds, yet is tightly bound and weakly interacting with itinerant electronic bands. These local moments, a rare occurrence for Nb, form one-dimensional spin chains that criss-cross perpendicularly oriented conducting nanowires. The local moment bandwidth is comparable to what is seen in rare earth compounds with extremely localized orbitals. The microscopic origin is traced to the local structure of the NbO6 octahedra and associated orbital+spin ordering. The resulting 1D$_s$*1D$_c$ anisotropic two dimensional Heisenberg-Kondo lattice model (s=spin, c=charge) provides a strongly anisotropic spin-fermion lattice system for further study. △ Less

Submitted 13 May, 2015; originally announced May 2015.

Comments: 6+ pages, 6 embedded figures

Journal ref: Phys. Rev. B 91, 195152 (2015)

Cubic H$_3$S around 200 GPa: an atomic hydrogen superconductor stabilized by sulfur

Authors: D. A. Papaconstantopoulos, B. M. Klein, M. J. Mehl, W. E. Pickett

Abstract: The multiple scattering-based theory of Gaspari and Gyorffy for the electron-ion matrix element in close packed metals is applied to Im-3m H3S, which has been predicted by Duan {\it et al.} and Bernstein {\it et al} to be the stable phase at this stoichiometry around 190 GPa, thus is the leading candidate to be the phase observed to superconduct at 190K by Drozdov, Eremets, and Troyan. The nearly… ▽ More The multiple scattering-based theory of Gaspari and Gyorffy for the electron-ion matrix element in close packed metals is applied to Im-3m H3S, which has been predicted by Duan {\it et al.} and Bernstein {\it et al} to be the stable phase at this stoichiometry around 190 GPa, thus is the leading candidate to be the phase observed to superconduct at 190K by Drozdov, Eremets, and Troyan. The nearly perfect separation of vibrational modes into those of S and of H character provides a simplification that enables identification of contributions of the two atoms separately. The picture that arises is basically that of superconducting atomic H stabilized by strong covalent mixing with S 3p and "3d" character. The reported isotope shift is much larger than the theoretical one,suggesting there is large anharmonicity in the H vibrations. Given the relative unimportance of sulfur, hydrides of lighter atoms at similarly high pressures may also lead to high temperature superconductivity. △ Less

Submitted 27 May, 2015; v1 submitted 16 January, 2015; originally announced January 2015.

Comments: 5 pages, 3 figures. The published version

Journal ref: Phys. Rev. B 91, 184511 (2015)

Spin-orbit coupling, strong correlation, and insulator-metal transitions: the J$_{\rm eff}$ =3\2 ferromagnetic Mott insulator Ba$_{2}$NaOsO$_{6}$

Authors: Shruba Gangopadhyay, Warren E. Pickett

Abstract: The double perovskite Ba$_{2}$NaOsO$_{6}$ (BNOO), an exotic example of a very high oxidation state (heptavalent) osmium $d^1$ compound and also uncommon by being a ferromagnetic Mott insulator without Jahn-Teller (JT) distortion, is modeled using the density functional theory (DFT) hybrid functional based exact exchange for correlated electrons (oeeHyb) method and including spin-orbit coupling (SO… ▽ More The double perovskite Ba$_{2}$NaOsO$_{6}$ (BNOO), an exotic example of a very high oxidation state (heptavalent) osmium $d^1$ compound and also uncommon by being a ferromagnetic Mott insulator without Jahn-Teller (JT) distortion, is modeled using the density functional theory (DFT) hybrid functional based exact exchange for correlated electrons (oeeHyb) method and including spin-orbit coupling (SOC). The experimentally observed narrow gap ferromagnetic insulating ground state is obtained, with easy axis along [110] in accord with experiment, providing support that this approach provides a realistic method for studying this system. The predicted spin density for [110] spin orientation is nearly cubic (unlike for other directions), providing an explanation for the absence of JT distortion. An orbital moment of -0.4$μ_B$ strongly compensates the +0.5$μ_B$ spin moment on Os, leaving a strongly compensated moment more in line with experiment. Remarkably, the net moment lies primarily on the oxygen ions. An insulator-metal transition by rotating the magnetization direction with an external field under moderate pressure is predicted as one consequence of strong SOC, and metallization under moderate pressure is predicted. Comparison is made with the isostructural, isovalent insulator Ba$_2$LiOsO$_6$ which however orders antiferromagnetically. △ Less

Submitted 12 October, 2014; originally announced October 2014.

Comments: 10 pages, 6 figues

Analysis of charge states in the mixed valent ionic insulator AgO

Authors: Yundi Quan, Warren E. Pickett

Abstract: The doubly ionized $d^9$ copper ion provides, originally in La$_2$CuO$_4$ and later in many more compounds, the platform for high temperature superconductivity when it is forced toward higher levels of oxidation. The nearest chemical equivalent is Ag$^{2+}$, which is almost entirely avoided in nature. AgO is an illustrative example, being an unusual nonmagnetic insulating compound with an open… ▽ More The doubly ionized $d^9$ copper ion provides, originally in La$_2$CuO$_4$ and later in many more compounds, the platform for high temperature superconductivity when it is forced toward higher levels of oxidation. The nearest chemical equivalent is Ag$^{2+}$, which is almost entirely avoided in nature. AgO is an illustrative example, being an unusual nonmagnetic insulating compound with an open $4d$ shell on one site. This compound has been interpreted in terms of one Ag$^{3+}$ ion at the fourfold site and one Ag$^{+}$ ion that is twofold coordinated. We analyze more aspects of this compound, finding that indeed the Ag$^{3+}$ ion supports only four occupied $4d$-based Wannier functions per spin, while Ag$^+$ supports five, yet their physical charges are nearly equal. The oxygen $2p$ Wannier functions display two distinct types of behavior, one type of which includes conspicuous Ag $4d$ tails. Calculation of the Born effective charge tensor shows that the mean effective charges of the Ag ions differ by about a factor of two, roughly consistent with the assigned formal charges. We analyze the $4d$ charge density and discuss it in terms of recent insights into charge states of insulating (and usually magnetic) transition metal oxides. What might be expected in electron- and hole-doped AgO is discussed briefly. △ Less

Submitted 22 September, 2014; originally announced September 2014.

Unquenched $e_g^1$ orbital moment in the Mott insulating antiferromagnet KOsO4

Authors: Young-Joon Song, Kyo-Hoon Ahn, K. -W. Lee, W. E. Pickett

Abstract: Applying the correlated electronic structure method based on density functional theory plus the Hubbard $U$ interaction, we have investigated the tetragonal scheelite structure Mott insulator KOsO$_4$, whose $e_g^1$ configuration should be affected only slightly by spin-orbit couping (SOC). The method reproduces the observed antiferromagnetic Mott insulating state, populating the Os $d_{z^2}$ majo… ▽ More Applying the correlated electronic structure method based on density functional theory plus the Hubbard $U$ interaction, we have investigated the tetragonal scheelite structure Mott insulator KOsO$_4$, whose $e_g^1$ configuration should be affected only slightly by spin-orbit couping (SOC). The method reproduces the observed antiferromagnetic Mott insulating state, populating the Os $d_{z^2}$ majority orbital. The quarter-filled $e_g$ manifold is characterized by a symmetry breaking due to the tetragonal structure, and the Os ion shows a crystal field splitting $Δ_{cf}$ = 1.7 eV from the $t_{2g}$ complex, which is relatively small considering the high formal oxidation state Os$^{7+}$. The small magnetocrystalline anisotropy before including correlation (i.e., in the metallic state) is increased by more than an order of magnitude in the Mott-insulating state, a result of a strong interplay between large SOC and a strong correlation. In contrast to conventional wisdom that the $e_g$ complex will not support orbital magnetism, we find that for the easy axis [100] direction the substantial Os orbital moment $M_L\approx-0.2 μ_B$ compensates half of the Os spin moment $M_S$ = 0.4$μ_B$. The origin of the orbital moment is analyzed and understood in terms of additional spin-orbital lowering of symmetry, and beyond that due to structural distortion, for magnetization along [100]. Further interpretation is assisted by analysis of the spin density and the Wannier function with SOC included. △ Less

Submitted 9 December, 2014; v1 submitted 18 August, 2014; originally announced August 2014.

Comments: 7 pages

Journal ref: Phys. Rev. B 90, 245117 (2014)

Dielectric Response of Electron-doped Ionic Superconductor Li$_x$ZrNCl

Authors: Antia S. Botana, Warren E. Pickett

Abstract: When electron doped, the layered transition metal nitrides ${\cal T}$NCl (${\cal T}$ = group IVB transition metal ion) become impressive superconductors with critical temperature T$_c$= 15-26K. Here we take the most studied member, ZrNCl, as a representative and calculate the dielectric response $ε(ω)$ versus frequency and concentration of doped electronic carriers. The static dielectric constant… ▽ More When electron doped, the layered transition metal nitrides ${\cal T}$NCl (${\cal T}$ = group IVB transition metal ion) become impressive superconductors with critical temperature T$_c$= 15-26K. Here we take the most studied member, ZrNCl, as a representative and calculate the dielectric response $ε(ω)$ versus frequency and concentration of doped electronic carriers. The static dielectric constant $ε_{\infty}$=5 is reproduced extremely well. We establish that the differences between rigid band modeling and virtual crystal treatment are small, and compare also with actual lithium doping using supercells. We obtain the variations upon changing the doping level of the reflectivity and energy loss function as well, many of which are found not to be correlated with the observed (non)variation of T$_c(x)$. The main plasmon peaks appear where the electron gas model suggests, in the range 1.2-2.0 eV for $x$ varying from 0.16 to 0.50. △ Less

Submitted 14 August, 2014; originally announced August 2014.

Strong interactions, narrow bands, and dominant spin-orbit coupling in Mott insulating quadruple perovskite CaCo3V4O12

Authors: H. B. Rhee, W. E. Pickett

Abstract: We investigate the electronic and magnetic structures and the character and direction of spin and orbital moments of the recently synthesized quadruple perovskite cobaltovanadate CaCo3V4O12 using a selection of methods from density functional theory. Implementing the generalized gradient approximation and the Hubbard U correction (GGA+U), ferromagnetic spin alignment leads to half-metallicity rath… ▽ More We investigate the electronic and magnetic structures and the character and direction of spin and orbital moments of the recently synthesized quadruple perovskite cobaltovanadate CaCo3V4O12 using a selection of methods from density functional theory. Implementing the generalized gradient approximation and the Hubbard U correction (GGA+U), ferromagnetic spin alignment leads to half-metallicity rather than the observed narrow gap insulating behavior. Including spin-orbit coupling (SOC) leaves a half-semimetallic spectrum which is essentially Mott insulating. SOC is crucial for the Mott insulating character of the V d^1 ion, breaking the d_{m=\pm 1} degeneracy and also giving a substantial orbital moment. Evidence is obtained of the large orbital moments on Co that have been inferred from the measured susceptibility. Switching to the orbital polarization (OP) functional, GGA+OP+SOC also displays clear tendencies toward very large orbital moments but in its own distinctive manner. In both approaches, application of SOC, which requires specification of the direction of the spin, introduces large differences in the orbital moments of the three Co ions in the primitive cell. We study a fictitious but simpler cousin compound Ca3CoV4O12 (Ca replacing two of the Co atoms) to probe in more transparent fashion the interplay of spin and orbital degrees of freedom with the local environment of the planar CoO4 units. The observation that the underlying mechanisms seems to be local to a CoO4 plaquette, and that there is very strong coupling of the size of the orbital moment to the spin direction, strongly suggest non-collinear spins, not only on Co but on the V sublattice as well. △ Less

Submitted 10 August, 2014; originally announced August 2014.

Comments: 9 pages, 6 figures, 3 tables

Journal ref: Phys. Rev. B 90, 205119 (2014)