Dataset of tensorial optical and transport properties of materials from the Wannier function method
Not provided.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Not provided.
This paper presents a fully ab initio many-body photoemission framework that includes coherent three-body electron-photon-phonon scattering to predict the transverse momentum distributions and the mean transverse energies (MTEs) of bulk photoelectrons from single-crystal photocathodes. The need to develop such a theory stems from the lack of studies that provide complete understanding of the underlying fundamental processes governing the transverse momentum distribution of photoelectrons emitted from single crystals. For example, initial predictions based on density-functional theory calculations of effective electron masses suggested that the (111) surface of PbTe would produce very small MTEs (≤ 15 meV), whereas our experiments yielded MTEs ten to twenty times larger than these predictions and also exhibited a lower photoemission threshold than predicted. The ab initio framework presented in this paper correctly reproduces both the magnitude of the MTEs from our measurements in PbTe(111) and the observed photoemission below the predicted threshold. Our results show that both photoexcitations into states that propagate in the bulk of the material and coherent many-body electron-photon-phonon scattering processes, which initial predictions ignored, play surprisingly important roles in photoemission from PbTe(111). Lastly, from the lessons learned, we recommend a procedure for rapid computational screening of potential single-crystal photocathodes for applications in next-generation ultrafast electron diffraction and X-ray free-electron lasers, which will enable significant advances in condensed matter research.
We derive a “Wannier-Hubbard” model consisting of an array of overlapping atomic orbitals interacting via a local Coulomb interaction. Transforming to an orthogonal Wannier basis set, the resulting Hamiltonian displays long-range hopping and interactions, with new terms such as correlated hopping and ferromagnetic direct exchange, among others. We numerically study the one-dimensional version of the model at half-filling using the density matrix renormalization group method, unveiling a rich phase diagram as a function of the interaction U and the overlap s with metallic and ferromagnetic phases separated by a ferrimagnetic region. Our results indicate a path toward understanding emergent phases under pressure and beyond standard model Hamiltonians.
Transition metals in body-centered cubic (bcc) structures under compression can display several novel physical properties because of their complex electronic structures and electron-phonon interactions. Here, we used inelastic x-ray scattering experiments in a diamond-anvil cell up to ∼45 GPa and density-functional theory calculations up to 210 GPa to investigate the phonon dispersions, and electronic and elastic properties of single-crystal molybdenum (Mo). Our results show a pressure-induced Kohn anomaly at 𝑞∼0.5 along the [ξ00] direction in the longitudinal acoustic mode at ∼45 GPa; this anomaly is triggered by the pressure-enhanced Fermi-surface nesting effect. Theoretical calculations show that electron redistributions in the 𝑠-to-𝑑 orbitals of bcc-Mo contribute to the shear modulus anomaly at ∼50 GPa. In contrast, the Young's modulus anomaly in bcc-Mo at ∼210 GPa results from a Lifshitz-type electronic topological transition. In conclusion, our results shed light on the complex electronic behaviors that are associated with macroscopic elastic properties in typical bcc 𝑑-block transition metals under compression.
Here, the nature of the superconducting order parameter in Sr 2 RuO 4 has generated intense interest in recent years. Since the superconducting gap is very small, high-resolution methods such as scanning tunneling spectroscopy might be the best chance to directly resolve the gap symmetry. Recently, a Bogoliubov quasiparticle interference imaging (BQPI) experiment has suggested that the d x 2 –y 2 gap symmetry is appropriate for Sr 2 RuO 4 . In this work, we use a material-specific theoretical approach based onWannier functions of the surface of Sr 2 RuO 4 to calculate the continuum density of states as detected in scanning tunneling microscopy experiments. We examine several different proposed gap order parameters and calculate the expected BQPI pattern for each case. Comparing to the available experimental data, our results suggest that a s' + id xy gap order parameter is the most probable state, but the measured BQPI patterns still display features unaccounted for by the theory for any of the states currently under discussion.
We report a first-principles density-functional study of electron-phonon interactions in and thermoelectric transport properties of the full Heusler compounds Sr 2 BiAu and Sr 2 SbAu. Our results show that ultrahigh intrinsic bulk thermoelectric performance across a wide range of temperatures is physically possible and point to the presence of multiply degenerate and highly dispersive carrier pockets as the key factor for achieving this. Sr 2 BiAu, which features ten energy-aligned low-effective-mass pockets (six along Γ-X and four at L), is predicted to deliver n-type zT=0.4-4.9 at T=100-700 K. Comparison with the previously investigated compound Ba 2 BiAu shows that the additional L pockets in Sr 2 BiAu significantly increase its low-temperature power factor to a maximum value of 12 mW m -1 K -2 near T=300 K. However, at high temperatures the power factor of Sr 2 BiAu drops below that of Ba 2 BiAu because the L states are heavier and subject to strong scattering by phonon deformation, as opposed to the lighter Γ-X states, which are limited by polar-optical scattering. Sr2SbAu is predicted to deliver a lower n-type zT=3.4 at T=750 K due to appreciable misalignment between the L and Γ-X carrier pockets, generally heavier scattering, and a slightly higher lattice thermal conductivity. Soft acoustic modes, which are responsible for the low lattice thermal conductivity, also increase the vibrational entropy and high-temperature stability of these Heusler compounds, suggesting that their experimental synthesis may be feasible. The dominant intrinsic defects are found to be Au vacancies, which drive the Fermi level towards the conduction band and work in favor of n-doping.
Since the discovery of pressure-induced superconductivity in the two-leg ladder system BaFe 2 X 3 (X=S, Se), with the 3d iron electronic density n=6, quasi-one-dimensional iron-based ladders have attracted considerable attention. In this work, we use density-functional theory to predict that the novel n=6 iron ladder BaFe 2 Se 3 could be stable with a similar crystal structure as BaFe 2 Se 3 . Furthermore, our results also indicate that BaFe2Te3 will display a complex 2×2 block-type magnetic order. Due to the magnetic striction effects of this block order, BaFe2Te3 should be a magnetic noncollinear ferrielectric system with a net polarization 0.31μC/cm2. In summary, the similar electronic density and magnetic ground state of Te- and Se-based ladders indicates both should display similar properties. In particular, the physical and structural similarity with BaFe 2 Se 3 suggests that BaFe 2 Se 3 could become superconducting under high pressure.
Localized bases play an important role in understanding electronic structure. In periodic insulators, a natural choice of localized basis is given by the Wannier functions which depend on a choice of unitary transform known as a gauge transformation. Over the past few decades, there have been many works that have focused on optimizing the choice of the gauge so that the corresponding Wannier functions are maximally localized or reflect some symmetry of the underlying system. In this work, we consider fully nonperiodic materials where the usual Wannier functions are not well defined and gauge optimization is impractical. To tackle the problem of calculating exponentially localized generalized Wannier functions in both periodic and nonperiodic systems, we discuss the ‘iterated projected position (IPP)” algorithm. Here, the IPP algorithm is based on matrix diagonalization and therefore unlike optimization-based approaches, it does not require initialization and cannot get stuck at a local minimum. Furthermore, the IPP algorithm is guaranteed by a rigorous analysis to produce exponentially localized functions under certain mild assumptions. We numerically demonstrate that the IPP algorithm can be used to calculate exponentially localized bases for the Haldane model, the Kane-Mele model (in both Z 2 invariant even and Z 2 invariant odd phases), and the p x + ip y model on a quasicrystal lattice.
In this work, we study the different ways of introducing light-matter interaction in first-principle tight-binding (TB) models. The standard way of describing optical properties is the velocity gauge, defined by linear coupling to the vector potential. In finite systems a transformation to represent the electromagnetic radiation by the electric field instead is possible, albeit subtleties arise in periodic systems. The resulting dipole gauge is a multi-orbital generalization of the Peierls substitution. In this work we investigate the accuracy for both pathways, with particular emphasis on gauge invariance, for TB models constructed from maximally localized Wannier functions. This approach accurately captures the light-matter interaction close to the Fermi level. Focusing on paradigmatic two-dimensional materials, we construct first-principle models and calculate the response to electromagnetic fields in linear response and for strong excitations. Benchmarks against fully converged first-principle calculations allow for ascertaining the accuracy of the TB models. We find that the dipole gauge provides a more accurate description than the velocity gauge in all cases. The main deficiency of the velocity gauge is an imperfect cancellation of paramagnetic and diamagnetic current. Formulating a corresponding sum rule however provides a way to explicitly enforce this cancellation. This procedure corrects the TB models in the velocity gauge, yielding excellent agreement with dipole gauge and thus restoring gauge invariance.
Co 3 Sn 2 S 2 is a magnetic Weyl semimetal, in which ferromagnetic ordering at 177 K is predicted to stabilize Weyl points. We perform temperature and spatial dependent angle--resolved photoemission spectroscopy measurements through the Curie temperature (T c ), which show large band shifts and renormalization concomitant with the onset of magnetism. We argue that Co 3 Sn 2 S 2 evolves from a Mott ferromagnet below T c to a correlated metallic state above T c . To understand the magnetism, we derive a tight-binding model of Co-3d x 2 -y 2 orbitals on the kagome lattice. At the filling obtained by first-principles calculations, this model reproduces the ferromagnetic ground state, and results in the reduction of Coulomb interactions due to cluster effects. Using a disordered local moment simulation, we show how this reduced Hubbard U leads to a collapse of the bands across the magnetic transition, resulting in a correlated state which carries associated characteristic photoemission signatures that are distinct from those of a simple lifting of exchange splitting. Finally, the behavior of topology across T c is discussed in the context of this description of the magnetism.
Here we report the results of our experimental studies on the magnetic, transport, and thermoelectric properties of the ferromagnetic metal CoMnSb. A sizable anomalous Hall conductivity $σ_{yx}$ and transverse thermoelectric conductivity $α_{yx}$ are found experimentally and comparable in size to the values estimated from first-principles calculation. Our experiment further reveals that CoMnSb exhibits -T log T critical behavior in $α_{yx}$ (T), deviating from Fermi liquid behavior $α_{yx}$ ~ T over a decade of temperature between 10 and 400 K, similar to ferromagnetic Weyl and nodal-line semimetals. Our theoretical calculation for CoMnSb also predicts the - T log T behavior when the Fermi energy locates near the Weyl nodes in momentum space.
Here we study the conditions for Bloch bands to be spanned by symmetric and strictly compact Wannier states that have zero overlap with all lattice sites beyond a certain range. Similar to the characterization of topological insulators in terms of an algebraic (rather than exponential) localization of Wannier states, we find that there may be impediments to the compact localization even of topologically “trivial” obstructed atomic insulators. These insulators admit exponentially localized Wannier states centered at unoccupied orbitals of the crystalline lattice. First, we establish a sufficient condition for an insulator to have a compact representative. Second, for $\mathscr{C}_2$ rotational symmetry, we prove that the complement of fragile topological bands cannot be compact, even if it is an atomic insulator. Third, for $\mathscr{C}_4$ symmetry, our findings imply that there exist fragile bands with zero correlation length. Fourth, for a $\mathscr{C}_3$-symmetric atomic insulator, we explicitly derive that there are no compact Wannier states overlapping with less than 18 lattice sites. We conjecture that this obstruction generalizes to all finite Wannier sizes. Our results can be regarded as the stepping stone to a generalized theory of Wannier states beyond dipole or quadrupole polarization.
Recent experimental papers on hole-doped overdoped cuprates have argued that a series of observations showing unexpected behavior in the superconducting state imply the breakdown of the quasiparticle-based Landau–BCS paradigm in that doping range. In contrast, some of the present authors have argued that a phenomenological “dirty d-wave” theoretical analysis explains essentially all aspects of thermodynamic and transport properties in the superconducting state, provided the unusual effects of weak, out-of-plane dopant impurities are properly accounted for. Here we attempt to place this theory on a more quantitative basis by performing ab initio calculations of dopant impurity potentials for LSCO and Tl-2201. These potentials are more complex than the pointlike impurity models considered previously, and require calculation of forward scattering corrections to transport properties. Including realistic, ARPES-derived band structures, Fermi liquid renormalizations, and vertex corrections, we show that the theory can explain semiquantitatively the unusual superfluid density measurements of the two most studied overdoped materials.
Ab initio and continuum model studies predicted that the Γ valley transition metal dichalcogenide (TMD) homobilayers could simulate the conventional multiorbital Hubbard model on the moiré honeycomb lattice. Here, in this work, we perform the Wannierization starting from the continuum model and show that a more general moiré Kanamori-Hubbard model emerges, beyond the extensively studied standard multiorbital Hubbard model, which can be used to investigate the many-body physics in the Γ valley TMD homobilayers. Using the unrestricted Hartree-Fock and Lanczos techniques, we study these half-filled multiorbital moiré bands. By constructing the phase diagrams we predict the presence of an antiferromagnetic state and in addition we found unexpected and dominant states, such as a S = 1 ferromagnetic insulator and a charge density wave state. Our theoretical predictions made using this model can be tested in future experiments on the Γ valley TMD homobilayers.
Here, using resonant inelastic x-ray scattering, we measure the paramagnon dispersion and damping of undoped, antiferromagnetic Ca 2 CuO 2 Cl 2 as well as doped, superconducting Na x Ca 2–x CuO 2 Cl 2 . Our estimation of the spin-exchange parameter and width of the paramagnon peak at the zone boundary X=(0.5, 0) confirms that no simple relation can be drawn between these parameters and the critical temperature T c . Consistently with other cuprate compounds, we show that upon doping there is a slight softening at (0.25,0) but not at the zone boundary X. In combination with these measurements we perform calculations of the dynamical spin structure factor of the one-band Hubbard model using cluster dynamical mean-field theory. The calculations are in excellent agreement with the experiment in the undoped case, both in terms of energy position and width. While the increase in width is also captured upon doping, the dynamical spin structure factor shows a sizable softening at X, which provides insightful information on the length-scale of the spin fluctuations in doped cuprates.
We explore the electronic and structural properties of Lu-H-N structures with specific N,H ordering due to their remarkably narrow hydrogen-based bands at the Fermi level. Although exhibits a classical instability persisting up to 17 GPa, it is anharmonically stable near ambient pressure when accounting for quantum nuclear effects. The presence of flat bands near is understood to arise from destructive quantum interference between and surrounding orbitals, with certain types of defects leaving the flat bands unaffected. The results suggest that if this system is superconducting, there is an optimal pressure near ambient where the superconducting is maximized by anharmonically stabilized low-frequency and nonadiabatically coupled high-frequency hydrogen modes. Despite the metastability of this structure, its dynamical stability when modeled beyond a classical harmonic approach suggests the narrow bands near do not cause structural distortions in certain Lu-H-N stoichiometries.
The possibility of selecting magnetic space groups by orienting the magnetization direction or tuning magnetic orders offers a vast playground for engineering symmetry-protected topological phases in magnetic materials. In this work, we study how selective tuning of symmetry and magnetism can influence and control the resulting topology in a two-dimensional magnetic system, and we illustrate such a procedure in the ferromagnetic monolayer MnPSe 3 . Density functional theory calculations reveal a symmetry-protected accidental semimetallic (SM) phase for out-of-plane magnetization, which becomes an insulator when the magnetization is tilted in-plane, reaching band-gap values close to 100 meV. We identify an order-2 composite antiunitary symmetry and threefold rotational symmetry that induce the band crossing, and we classify the possible topological phases using symmetry analysis, which we support with tight-binding and k · p models. Breaking of inversion symmetry opens a gap in the SM phase, giving rise to a Chern insulator. We demonstrate this explicitly in the isostructural Janus compound Mn 2 P 2 S 3 Se 3 , which naturally exhibits Rashba spin-orbit coupling that breaks inversion symmetry. Our results map out the phase space of topological properties of ferromagnetic transition-metal phosphorus trichalcogenides, and they demonstrate the potential of the magnetization-dependent metal-to-insulator transition as a spin switch in integrated two-dimensional electronics.
Motivated by the recently proposed alternating single-layer trilayer stacking structure for the nickelate La 3 Ni 2 O 7 , we comprehensively study this system using ab initio and random-phase approximation techniques. Here, our analysis unveils similarities between this novel La 3 Ni 2 O 7 structure and other Ruddlesden-Popper nickelate superconductors, such as a similar charge-transfer gap value and orbital-selective behavior of the eg orbitals. Pressure primarily increases the bandwidths of the Ni e g bands, suggesting an enhancement of the itinerant properties of those e g states. By changing the cell volume ratio V/V 0 from 0.9 to 1.10, we found that the bilayer structure in La 3 Ni 2 O 7 always has lower energy than the single-layer trilayer stacking La 3 Ni 2 O 7 . In addition, we observe a “self-doping” effect (compared to the average 1.5 electrons per eg orbital per site of the entire structure) from the trilayer to the single-layer sublattices and this effect will be enhanced by overall electron doping. Moreover, we find a leading d x 2 -y 2 -wave pairing state that is restricted to the single layer. Because the effective coupling between the single layers is very weak, due to the nonsuperconducting trilayer in-between, this suggests that the superconducting transition temperature T c in this structure should be much lower than in the bilayer structure.