Search NASA⌕ Search

SEARCH · Search NASA

Results for “Wannier functions”

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.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Nonempirical Prediction of the Length-Dependent Ionization Potential in Molecular Chains

The ionization potential of molecular chains is well-known to be a tunable nanoscale property that exhibits clear quantum confinement effects. State-of-the-art methods can accurately predict the ionization potential in the small molecule limit and in the solid-state limit, but for intermediate, nanosized systems prediction of the evolution of the electronic structure between the two limits is more difficult. Recently, optimal tuning of range-separated hybrid functionals has emerged as a highly accurate method for predicting ionization potentials. This was first achieved for molecules using the ionization potential theorem (IPT) and more recently extended to solid-state systems, based on an ansatz that generalizes the IPT to the removal of charge from a localized Wannier function. Here, we study one-dimensional molecular chains of increasing size, from the monomer limit to the infinite polymer limit using this approach. By comparing our results with other localization-based methods and where available with experiment, we demonstrate that Wannier-localization-based optimal tuning is highly accurate in predicting ionization potentials for any chain length, including the nanoscale regime.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Modeling Liquid Water by Climbing up Jacob’s Ladder in Density Functional Theory Facilitated by Using Deep Neural Network Potentials

Within the framework of Kohn–Sham density functional theory (DFT), the ability to provide good predictions of water properties by employing a strongly constrained and appropriately normed (SCAN) functional has been extensively demonstrated in recent years. Here, we further advance the modeling of water by building a more accurate model on the fourth rung of Jacob’s ladder with the hybrid functional, SCAN0. In particular, we carry out both classical and Feynman path-integral molecular dynamics calculations of water with the SCAN0 functional and the isobaric–isothermal ensemble. To generate the equilibrated structure of water, a deep neural network potential is trained from the atomic potential energy surface based on ab initio data obtained from SCAN0 DFT calculations. For the electronic properties of water, a separate deep neural network potential is trained by using the Deep Wannier method based on the maximally localized Wannier functions of the equilibrated trajectory at the SCAN0 level. The structural, dynamic, and electric properties of water were analyzed. The hydrogen-bond structures, density, infrared spectra, diffusion coefficients, and dielectric constants of water, in the electronic ground state, are computed by using a large simulation box and long simulation time. For the properties involving electronic excitations, we apply the GW approximation within many-body perturbation theory to calculate the quasiparticle density of states and bandgap of water. Compared to the SCAN functional, mixing exact exchange mitigates the self-interaction error in the meta-generalized-gradient approximation and further softens liquid water toward the experimental direction. For most of the water properties, the SCAN0 functional shows a systematic improvement over the SCAN functional. However, some important discrepancies remain. The H-bond network predicted by the SCAN0 functional is still slightly overstructured compared to the experimental results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

pyGWBSE: a high throughput workflow package for GW-BSE calculations

Abstract We develop an open-source python workflow package, py GWBSE to perform automated first-principles calculations within the GW-BSE (Bethe-Salpeter) framework. GW-BSE is a many body perturbation theory based approach to explore the quasiparticle (QP) and excitonic properties of materials. GW approximation accurately predicts bandgaps of materials by overcoming the bandgap underestimation issue of the more widely used density functional theory (DFT). BSE formalism produces absorption spectra directly comparable with experimental observations. py GWBSE package achieves complete automation of the entire multi-step GW-BSE computation, including the convergence tests of several parameters that are crucial for the accuracy of these calculations. py GWBSE is integrated with Wannier90 , to generate QP bandstructures, interpolated using the maximally-localized wannier functions. py GWBSE also enables the automated creation of databases of metadata and data, including QP and excitonic properties, which can be extremely useful for future material discovery studies in the field of ultra-wide bandgap semiconductors, electronics, photovoltaics, and photocatalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance

By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car–Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson’s equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. Furthermore, this is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500–1000 atoms (e.g., (H₂O)₂₅₆) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

First-principles calculation of gate-tunable ferromagnetism in magic-angle twisted bilayer graphene under pressure

Magic-angle twisted bilayer graphene (MATBG) is notable as a highly tunable platform for investigating strongly correlated phenomena such as unconventional superconductivity and quantum spin liquids, due to easy control of doping level through gating and sensitive dependence of the magic angle on hydrostatic pressure. Experimental observations of correlated insulating states, unconventional superconductivity and ferromagnetism in MATBG indicate that this system exhibits rich exotic phases. In this work, using density functional theory calculations in conjunction with the effective screening medium method, we find the MATBG under pressure at a twisting angle of 2.88 °and simulate how its electronic states evolve when doping level and electric field perpendicular to plane are tuned by gating. Our calculations show that, at doping levels between two electrons and four holes per moiré unit cell, a ferromagnetic (FM) solution with spin density localized at AA stacking sites is lower in energy than the nonmagnetic solution. The magnetic moment of this FM state decreases with both electron and hole doping and vanishes at four electrons/holes doped per moiré unit cell. Hybridization between the flat bands at the Fermi level and the surrounding dispersive bands can take place at finite doping. On increasing the out-of-plane electric field at zero doping, a transition from the FM state to the nonmagnetic one is seen. An investigation of impurity effects shows that both absorption of O 2 molecules and occurrence of Stone–Wales impurities suppress the FM state, and the mechanisms are understood from our calculations. We also analyze the interlayer bonding character due to flat bands via Wannier functions. Finally, we report trivial band topology of the flat bands in the FM state at a certain doping level.

Physics↗
DOE OSTI2022

Strongly anisotropic electronic and magnetic structures in oxide dichlorides RuOCl 2 and OsOCl 2

The van der Waals oxide dichlorides MOX 2 (M=V, Ta, Nb, Ru, and Os; X=halogen element), with different electronic densities, are attracting considerable attention. Ferroelectricity, spin-singlet formation, and orbital-selective Peierls phases were reported in this family with d 1 or d 2 electronic configurations, all believed to be caused by the strongly anisotropic electronic orbital degree of freedom. Here, using density functional theory and density matrix renormalization group methods, we investigate the electronic and magnetic properties of RuOCl 2 and OsOCl 2 with d 4 electronic configurations. Different from a previous study using VOI 2 with d 1 configuration, these systems with 4d4 or 5d4 do not exhibit a ferroelectric instability along the a axis. Due to the fully occupied dxy orbital in RuOCl 2 and OsOCl 2 , the Peierls instability distortion disappears along the b axis, leading to an undistorted Immm phase (No. 71). Furthermore, we observe strongly anisotropic electronic and magnetic structures along the a axis. For this reason, the materials of our focus can be regarded as “effective one-dimensional” systems even when they apparently have a dominant two-dimensional lattice geometry. The large crystal-field splitting energy (between dxz/yz and dxy orbitals) and large hopping between nearest-neighbor Ru and Os atoms suppresses the J = 0 singlet state in MOCl 2 (M=Ru or Os) with electronic density n = 4, resulting in a spin-1 system. Moreover, we find staggered antiferromagnetic order with π wave vector along the M-O chain direction (a axis) while the magnetic coupling along the b axis is weak. Based on Wannier functions from first-principles calculations, we calculated the relevant hopping amplitudes and crystal-field splitting energies of the t 2g orbitals for the Os atoms to construct a multiorbital Hubbard model for the M-O chains. Staggered AFM with ↑|-|↓-↑|-|↓ spin structure dominates in our density matrix renormalization group calculations, in agreement with density functional theory calculations. Our results for RuOCl 2 and OsOCl 2 provide guidance to experimentalists and theorists working on this interesting family of oxide dichlorides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Molecular Simulation of Covalent Bond Dynamics in Liquid Silicon

Many atomic liquids can form transient covalent bonds reminiscent of those in the corresponding solid states. These directional interactions dictate many important properties of the liquid state, necessitating a quantitative, atomic-scale understanding of bonding in these complex systems. A prototypical example is liquid silicon, wherein transient covalent bonds give rise to local tetrahedral order and consequent nontrivial effects on liquid-state thermodynamics and dynamics. To further understand covalent bonding in liquid silicon, and similar liquids, we present an ab initio-simulation-based approach for quantifying the structure and dynamics of covalent bonds in condensed phases. Through the examination of structural correlations among silicon nuclei and maximally localized Wannier function centers, we develop a geometric criterion for covalent bonds in liquid Si. We use this to monitor the dynamics of transient covalent bonding in the liquid state and estimate a covalent bond lifetime. We compare covalent bond dynamics to other processes in liquid Si and similar liquids and suggest experiments to measure the covalent bond lifetime.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Linear Combination of Atomic Dipoles to Calculate the Bond and Molecular Dipole Moments of Molecules and Molecular Liquids

Here we report a linear combination of atomic dipole (LCAD) method for calculating the bond dipole moments of molecules. We show that the LCAD method reproduces the known molecular dipole moments of small to large molecules with a small error with respect to experimental and benchmark ab initio calculations, and molecular dipole distributions of bulk water that agree with maximally localized Wannier functions. The bond dipole moments derived from LCAD are also chemically interpretable in terms of the trend in bond ionicity in going from neutral to charged molecules. Moreover, the LCAD method accurately captures the influence of electric fields, supported by the correct trend in the change of the dipole moment under a uniform external electric field. The better grounding of bond dipole calculations indicates that it should also serve as a useful approach to bond dipole-field models used in catalysis or to reconstruct the small dipole of a H-terminated graphene flake.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Topological flat bands in frustrated kagome lattice CoSn

Electronic flat bands in momentum space, arising from strong localization of electrons in real space, are an ideal stage to realize strongly-correlated phenomena. Theoretically, the flat bands can naturally arise in certain geometrically frustrated lattices, often with nontrivial topology if combined with spin-orbit coupling. Here, we report the observation of topological flat bands in frustrated kagome metal CoSn, using angle-resolved photoemission spectroscopy and band structure calculations. Throughout the entire Brillouin zone, the bandwidth of the flat band is suppressed by an order of magnitude compared to the Dirac bands originating from the same orbitals. The frustration-driven nature of the flat band is directly confirmed by the chiral d-orbital texture of the corresponding real-space Wannier functions. Spin-orbit coupling opens a large gap of 80 meV at the quadratic touching point between the Dirac and flat bands, endowing a nonzero Z2 invariant to the flat band. These findings demonstrate that kagome-derived flat bands are a promising platform for novel emergent phases of matter at the confluence of strong correlation and topology.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Origin of charge density wave in topological semimetals SrAl 4 and EuAl 4

Topological semimetals in BaAl 4 -type structure show many interesting behaviors, such as charge density wave (CDW) in SrAl 4 and EuAl 4 , but not the isostructural and isovalent BaAl 4 , SrGa 4 , and BaGa 4 . Here using Wannier functions based on density functional theory, we calculate the susceptibility functions with millions of k-points to reach the small q-vector and study the origin and driving force behind the CDW. Our comparative study reveals that the origin of the CDW in SrAl 4 and EuAl 4 is the strong electron-phonon coupling interaction for the transverse acoustic mode at small q-vector along the Γ-Z direction besides the maximum of the real part of the susceptibility function from the nested Fermi surfaces of the Dirac-like bands, which explains well the absence of CDW in the other closely related compounds in a good agreement with experiment. We also connect the different CDW behaviors in the Al compounds to the macroscopic elastic properties.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Microscopic scattering approach to in-gap states: Cr adatoms on superconducting 𝛽−Bi 2 ⁢Pd

We develop a microscopic scattering formalism to describe Yu-Shiba-Rusinov (YSR) states due to a single Cr adatom on the Bi-terminated surface of 𝛽−Bi 2 ⁢Pd, by combining ab initio Wannier functions with a real-space Green's function approach in the Bogoliubov–de Gennes formalism. Our framework reproduces key scanning tunneling spectroscopy features, including a single particle-hole asymmetric YSR peak and isotropic 𝑑⁢𝐼/𝑑⁢𝑉 maps around the impurity. Decomposing the YSR states reveals contributions from four nearly degenerate 𝐶 4⁢𝑣 representations, with energy broadening masking their individual signatures. Spin-orbit coupling induces partial spin polarization, while the spatial asymmetry between particle and hole components arises from Cr 𝑑−Bi 𝑝 hybridization. These results highlight the importance of realistic band structures and microscopic modeling for interpreting STM data and provide a foundation for studying impurity chains hosting topological excitations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Topology invisible to eigenvalues in obstructed atomic insulators

Here we consider the extent to which symmetry eigenvalues reveal the topological character of bands. Specifically, we compare distinct atomic limit phases (band representations) that share the same irreducible representations (irreps) at all points in the Brillouin zone and, therefore, appear equivalent in a classification based on eigenvalues. We derive examples where such “irrep-equivalent” phases can be distinguished by a quantized Berry phase or generalization thereof. These examples constitute a generalization of the Su-Schrieffer-Heeger chain: neither phase is topological, in the sense that localized Wannier functions exist, yet there is a topological obstruction between them. We refer to two phases as “Berry obstructed atomic limits” if they have the same irreps, but differ by Berry phases. This is a distinct notion from eigenvalue obstructed atomic limits, which differ in their symmetry irreps at some point in the Brillouin zone. We compute exhaustive lists of elementary band representations that are irrep equivalent, in all space groups, with and without time-reversal symmetry and spin-orbit coupling, and use group theory to derive a set of necessary conditions for irrep equivalence. Finally, we conjecture, and in some cases prove, that irrep-equivalent elementary band representations that are not equivalent can be distinguished by a topological invariant.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effective one-band models for the one-dimensional cuprate Ba 2-x Sr x CuO 3+δ

In this work, we consider a multiband Hubbard model H m for Cu and O orbitals in Ba 2-x Sr x CuO 3+δ similar to the three-band model for two-dimensional cuprates. The hopping parameters are obtained from maximally localized Wannier functions derived from ab initio calculations. Using the cell perturbation method, we derive both a generalized t–J model H tJ and a one-band Hubbard model H H to describe the low-energy physics of the system. H tJ has the advantage of having a smaller relevant Hilbert space, facilitating numerical calculations, while additional terms should be included in H H to accurately describe the multiband physics of H m . Using H tJ and the density matrix renormalization group method, we calculate the wave-vector-resolved photoemission and discuss the relevant features in comparison with recent experiments. In agreement with previous calculations, we find that the addition of an attractive nearest-neighbor interaction of the order of the nearest-neighbor hopping shifts the weight from the 3k F to the holon-folding branch. Kinetic effects also contribute to this process.

1-dimensional systems↗
DOE OSTI2023

Coupled topological flat and wide bands: Quasiparticle formation and destruction

Flat bands amplify correlation effects and are of extensive current interest. They provide a platform to explore both topology in correlated settings and correlation physics enriched by topology. Recent experiments in correlated kagome metals have found evidence for strange-metal behavior. A major theoretical challenge is to study the effect of local Coulomb repulsion when the band topology obstructs a real-space description. In a variant to the kagome lattice, we identify an orbital-selective Mott transition in any system of coupled topological flat and wide bands. This was made possible by the construction of exponentially localized and Kramers-doublet Wannier functions, which, in turn, leads to an effective Kondo-lattice description. Our findings show how quasiparticles are formed in such coupled topological flat-wide band systems and, equally important, how they are destroyed. Our work provides a conceptual framework for the understanding of the existing and emerging strange-metal properties in kagome metals and beyond.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Microscopic Scattering Approach to In-Gap States

We develop a microscopic scattering formalism to describe Yu-Shiba-Rusinov (YSR) states due to a single Cr adatom on the Bi-terminated surface of beta Bi2Pd, by combining ab initio Wannier functions with a real-space Green's function approach in the Bogoliubov-de Gennes formalism[1]. Our framework reproduces key scanning tunneling spectroscopy features, including a single particle-hole asymmetric YSR peak and isotropic dIdV maps around the impurity. Decomposing the YSR states reveals contributions from four nearly degenerate C4v representations, with energy broadening masking their individual signatures. Spin-orbit coupling induces partial spin polarization, while the spatial asymmetry between particle and hole components arises from Cr d-Bi p hybridization. These results highlight the importance of realistic band structures and microscopic modeling for interpreting STM data for magnetic in-gap states on superconductors. Further advances examining layered 2D material surfaces, such as NbSe2, will be described[2]. For this system the superconducting properties are obtained from a full anisotropic Eliashberg calculation of the superconducting order parameter along with the charge density wave gap. Additional features associated with proposals to measure the dynamics of these individual YSR states will be presented. [1] arXiv:2507.08740 [2] arXiv:2507.11856

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRSDate not supplied

Downfolding Complex Materials Problems Onto Model Hamiltonians for Quantum Computers

Simulating the properties of quantum materials is expected to be one of the exciting applications for quantum computers and where we hope to see advantages over classical hardware. The complexity of ab initio Hamiltonians describing the physics of application-relevant materials places them beyond the realm of possibility for solution on near-term hardware with a limited number of qubits. Various Hamiltonian approximations, including Hamiltonian downfolding, offers a possibility towards simulating complex materials on near-term hardware. This is accomplished by approximating the relevant physics of a given material through their representation by simpler model Hamiltonians, such as the Hubbard Hamiltonian or extensions of it. Here we employ a well-defined first-principles methodology for deriving downfolded multi-band extended Hubbard Hamiltonians of materials, capturing strong electronic correlation and electron-phonon coupling, based on the formalism of Wannier functions and the calculation of the screened Coulomb interaction. We demonstrate for a variety of systems that quantum simulation of these downfolded Hamiltonians reproduces key properties, thus establishing downfolding as a promising route to achieve near-term simulation of application-relevant systems on quantum hardware.

Antonios Markos Alvertis↗
NASA NTRSDate not supplied

Downfolding complex materials problems onto model Hamiltonians for quantum computers

Simulating the properties of quantum materials is expected to be one of the exciting applications for quantum computers and where we hope to see advantages over classical hardware. The complexity of ab initio Hamiltonians describing the physics of application-relevant materials places them beyond the realm of possibility for solution on near-term hardware with a limited number of qubits. Various Hamiltonian approximations, including Hamiltonian downfolding, offers a possibility towards simulating complex materials on near-term hardware. This is accomplished by approximating the relevant physics of a given material through their representation by simpler model Hamiltonians, such as the Hubbard Hamiltonian or extensions of it. Here we employ a well-defined first-principles methodology for deriving downfolded multi-band extended Hubbard Hamiltonians of materials, capturing strong electronic correlation and electron-phonon coupling, based on the formalism of Wannier functions and the calculation of the screened Coulomb interaction. We demonstrate for a variety of systems that quantum simulation of these downfolded Hamiltonians reproduces key properties, thus establishing downfolding as a promising route to achieve near-term simulation of application-relevant systems on quantum hardware.

Antonios M. Alvertis↗