Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wannier function methods”

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 OSTI2025

Thermal disorder and phonon softening in the ferroelectric phase transition of lead titanate

We report a molecular dynamics study of ab initio quality of the ferroelectric phase transition in crystalline PbTi⁢O3. We model anharmonicity accurately in terms of potential energy and polarization surfaces trained on density functional theory data with modern machine learning techniques. Our simulations demonstrate that the transition has a strong order-disorder character, in agreement with diffraction experiments, and provide fresh insight into the approach to equilibrium across the phase transition. We find that the emergence and disappearance of the macroscopic polarization is driven by dipolar switching at the nanometer scale. We also computed the infrared optical absorption spectra in both the ferroelectric and the paraelectric phases, finding good agreement with the experimental Raman frequencies. Often, the almost ideal displacive character of the soft mode detected by Raman scattering in the paraelectric phase has been contrasted with the order-disorder character of the transition suggested by diffraction experiments. We settle this issue by showing that the soft mode coexists with a strong Debye relaxation associated with thermal disordering of the dipoles. The Debye relaxation feature is centered at zero frequency and appears near the transition temperature in both the ferroelectric and the paraelectric phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Lone Pair Rotational Dynamics in Solids

Traditional classifications of crystalline phases focus on nuclear degrees of freedom. Through examination of both electronic and nuclear structure, we introduce the concept of an electronic plastic crystal. Such a material is classified by crystalline nuclear structure, while localized electronic degrees of freedom — here lone pairs — exhibit orientational motion at finite temperatures. This orientational motion is an emergent phenomenon arising from the coupling between electronic structure and polarization fluctuations generated by collective motions, such as phonons. Using ab initio molecular dynamics simulations, we predict the existence of electronic plastic crystal motion in halogen crystals and halide perovskites, and suggest that such motion may be found in a broad range of solids with lone pair electrons. Finally, such fluctuations in the charge density should be observable, in principle via synchrotron scattering.

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

Ab Initio Electron-Phonon Interactions in Correlated Electron Systems

Electron-phonon ( e - ph ) interactions are pervasive in condensed matter, governing phenomena such as transport, superconductivity, charge-density waves, polarons, and metal-insulator transitions. First-principles approaches enable accurate calculations of e - ph interactions in a wide range of solids. However, they remain an open challenge in correlated electron systems (CES), where density functional theory often fails to describe the ground state. Therefore reliable e - ph calculations remain out of reach for many transition metal oxides, high-temperature superconductors, Mott insulators, planetary materials, and multiferroics. Here we show first-principles calculations of e - ph interactions in CES, using the framework of Hubbard-corrected density functional theory ( DFT + U ) and its linear response extension ( DFPT + U ), which can describe the electronic structure and lattice dynamics of many CES. We showcase the accuracy of this approach for a prototypical Mott system, CoO, carrying out a detailed investigation of its e - ph interactions and electron spectral functions. While standard DFPT gives unphysically divergent and short-ranged e - ph interactions, DFPT + U is shown to remove the divergences and properly account for the long-range Fröhlich interaction, allowing us to model polaron effects in a Mott insulator. Our work establishes a broadly applicable and affordable approach for quantitative studies of e - ph interactions in CES, a novel theoretical tool to interpret experiments in this broad class of materials.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Prediction of s ± -Wave Superconductivity Enhanced by Electronic Doping in Trilayer Nickelates La 4 Ni 3 O 10 under Pressure

Motivated by the recently reported signatures of superconductivity in trilayer La 4 ⁢Ni 3 ⁢O 10 under pressure, here we comprehensively study this system using ab initio and random-phase approximation techniques. Without electronic interactions, the Ni d 3z 2 –r 2 orbitals show a bonding-antibonding and nonbonding splitting behavior via the O p z orbitals inducing a “trimer” lattice in La 4⁢ Ni 3 ⁢O 10 , analogous to the dimers of La 3 ⁢Ni 2 ⁢O 7 . The Fermi surface consists of three electron sheets with mixed e g orbitals, and a hole and an electron pocket made up of the d 3⁢z 2 –r 2 orbital, suggesting a Ni two-orbital minimum model. In addition, we find that superconducting pairing is induced in the s ± -wave channel due to partial nesting between the M = (π,π) centered pockets and portions of the Fermi surface centered at the Γ = (0,0) point. With changing electronic density n, the s ± instability remains leading and its pairing strength shows a domelike behavior with a maximum around n = 4.2 ( ~6.7% electron doping). The superconducting instability disappears at the same electronic density as that in the new 1313 stacking La 3 ⁢Ni 2 ⁢O 7 , correlated with the vanishing of the hole pocket that arises from the trilayer sublattice, suggesting that the high-T c superconductivity of La 3⁢ Ni 2 ⁢O 7 does not originate from a trilayer and monolayer structure. Furthermore, we confirm the experimentally proposed spin state in La 4 ⁢Ni 3 ⁢O 10 with an in-plane (π, π) order and antiferromagnetic coupling between the top and bottom Ni layers, and spin zero in the middle layer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

First-principles derivation of magnetic interactions in the triangular quantum spin liquid candidates KYb C h 2 ( C h = S , Se , Te ) and A YbSe 2 ( A = Na , Rb )

The AYb Ch 2 (A=alkali metal, Ch=chalcogen) delafossites are a family of crystals which have a triangular lattice of effective S=1/2 moments on their Yb atoms. They hold great promise for the realization of the triangular quantum spin liquid state because of their defect-minimized growth and the ability to interchange chemical constituencies among the family. Here we use ab initio computations to evaluate the exchange couplings of four realized (and one theoretical) rhombohedral delafossite structures and examine the influence of chemical substitution on promoting the development of a quantum spin liquid state. We find good agreement with experiment regarding the antiferromagnetic nearest-neighbor exchange J 1 , but our calculations underestimate J 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Superfluid Weight Bounds from Symmetry and Quantum Geometry in Flat Bands

Flat-band superconductivity has theoretically demonstrated the importance of band topology to correlated phases. In two dimensions, the superfluid weight, which determines the critical temperature through the Berezinksii-Kosterlitz-Thouless criteria, is bounded by the Fubini-Study metric at zero temperature. We show this bound is nonzero within flat bands whose Wannier centers are obstructed from the atoms—even when they have identically zero Berry curvature. Next, we derive general lower bounds for the superfluid weight in terms of momentum space irreps in all 2D space groups, extending the reach of topological quantum chemistry to superconducting states. We find that the bounds can be naturally expressed using the formalism of real space invariants (RSIs) that highlight the separation between electronic and atomic degrees of freedom. Finally, using exact Monte Carlo simulations on a model with perfectly flat bands and strictly local obstructed Wannier functions, we find that an attractive Hubbard interaction results in superconductivity as predicted by the RSI bound beyond mean field. Hence, obstructed bands are distinguished from trivial bands in the presence of interactions by the nonzero lower bound imposed on their superfluid weight.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Wannier–Koopmans method calculations for transition metal oxide band gaps

The widely used density functional theory (DFT) has a major drawback of underestimating the band gaps of materials. Wannier–Koopmans method (WKM) was recently developed for band gap calculations with accuracy on a par with more complicated methods. WKM has been tested for main group covalent semiconductors, alkali halides, 2D materials, and organic crystals. Here we apply the WKM to another interesting type of material system: the transition metal (TM) oxides. TM oxides can be classified as either with d 0 or d 10 closed shell occupancy or partially occupied open shell configuration, and the latter is known to be strongly correlated Mott insulators. We found that, while WKM provides adequate band gaps for the d0 and d10 TM oxides, it fails to provide correct band gaps for the group with partially occupied d states. This issue is also found in other mean-field approaches like the GW calculations. We believe that the problem comes from a strong interaction between the occupied and unoccupied d-state Wannier functions in a partially occupied d-state system. We also found that, for pseudopotential calculations including deep core levels, it is necessary to remove the electron densities of these deep core levels in the Hartree and exchange–correlation energy functional when calculating the WKM correction parameters for the d-state Wannier functions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Real-Space Quantification of Exciton Localization in Acene Crystals Using Wannier Function Decomposition

Here, we introduce the Wannier function decomposition of excitons (WFDX) method to quantify exciton localization in solids within the ab initio Bethe-Salpeter equation framework. By decomposing each Bloch exciton wave function into products of single-particle electron and hole maximally localized Wannier functions, this real-space approach provides well-defined orbital- and spatial-resolved measures of both Frenkel and charge-transfer excitons at low computational cost. We apply WFDX to excitons in acene crystals, quantifying how the number of rings, the exciton spin state, and the center-of-mass momentum affect spatial localization. Additionally, we show how this real-space representation reflects structural nonsymmorphic symmetries that are hidden in standard reciprocal-space descriptions. We demonstrate how the WFDX framework can be used to efficiently interpolate exciton expansion coefficients in reciprocal-space and outline how it may facilitate evaluation of observables involving position operators, highlighting its potential as a general tool for both analyzing and computing excitonic properties in solids.

Tao, Zui [University of California, Berkeley, CA (↗
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 OSTI2022

Competing magnetic and nonmagnetic states in monolayer VSe 2 with charge density wave

The field of two-dimensional ferromagnets has been reinvigorated by the discovery of VSe 2 monolayer grown on van der Waals substrates, which is reported to be ferromagnetic with a Curie point higher than 330 K. However, the ferromagnetic and nonmagnetic states of pristine monolayer VSe 2 are highly debated. Here, employing density functional theory, Wannier function calculations, and the band unfolding method, we explore the electronic structure of monolayer VSe 2 with a $\sqrt{3} × \sqrt{7}$ charge density wave (CDW). Certain qualitative aspects of the calculated unfolded band dispersion and unfolded Fermi surface of monolayer VSe 2 with the $\sqrt{3} × \sqrt{7}$ CDW in the nonmagnetic state agree well with previous angle-resolved photoemission spectroscopy results, albeit with uncertainty about whether these experiments probed single or multiple domains. Specifically, we find that an isolated CDW domain naturally induces a strong breaking of the threefold symmetry of the electronic structure. In addition we find that, relative to the undistorted structure, the CDW structure shows a strong competition between nonmagnetic and various magnetic states, with an energy difference less than 5 meV/f.u. For the CDW structure in the antiferromagnetic state, the band dispersions and Fermi surface are similar to those in the nonmagnetic state, while the unfolded bands of the ferromagnetic CDW state display a sizable exchange splitting. These results indicate the possibility of various antiferromagnetic fluctuations in VSe 2 to coexist and compete with ferromagnetic order and the experimentally reported CDW order. In this work, our calculations build insights for exploring the interplay between magnetism and CDW behaviors more generally in transition metal dichalcogenides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 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

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 OSTI2021

Prediction of intrinsic topological superconductivity in Mn-doped GeTe monolayer from first-principles

The recent discovery of topological superconductors (TSCs) has sparked enormous interest. The realization of TSC requires a delicate tuning of multiple microscopic parameters, which remains a great challenge. Here, we develop a first-principles approach to quantify realistic conditions of TSC by solving self-consistently Bogoliubov-de Gennes equation based on a Wannier function construction of band structure, in presence of Rashba spin-orbit coupling, Zeeman splitting and electron-phonon coupling. We further demonstrate the power of this method by predicting the Mn-doped GeTe (Ge 1- x Mn x Te) monolayer—a well-known dilute magnetic semiconductor showing superconductivity under hole doping—to be a Class D TSC with Chern number of -1 and chiral Majorana edge modes. By constructing a first-principles phase diagram in the parameter space of temperature and Mn concentration, we propose the TSC phase can be induced at a lower-limit transition temperature of ~40 mK and the Mn concentration of x ~0.015%. Our approach can be generally applied to TSCs with a phonon-mediated pairing, providing useful guidance for future experiments.

Chemistry↗
DOE OSTI2022

Optimally tuned starting point for single-shot GW calculations of solids

The dependence of ab initio many-body perturbation theory within the $\textit{GW}$ approximation on the eigensystem used in calculating quasiparticle corrections limits this method's predictive power. Here, we investigate the accuracy of the recently developed Wannier-localized optimally tuned screened range-separated hybrid (WOT-SRSH) functional as a generalized Kohn-Sham starting point for single-shot $\textit{GW}$ ($G_0W_0$) calculations for a range of semiconductors and insulators. Comparison to calculations based on well-established functionals, namely, PBE, PBE0, and HSE, as well as to self-consistent GW schemes and to experiment, shows that band gaps computed via $G_0W_0$@WOT-SRSH have a level of precision and accuracy that is comparable to that of more advanced methods such as quasiparticle self-consistent $\textit{GW}$ and eigenvalue self-consistent $\textit{GW}$. In this work, we also find that $G_0W_0$@WOT-SRSH improves the description of states deeper in the valence band manifold. Finally, we show that $G_0W_0$@WOT-SRSH significantly reduces the sensitivity of computed band gaps to ambiguities in the underlying WOT-SRSH tuning procedure.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Impact of hole polaron formation on excitonic transitions in MgO from first principles

Here, we present a first-principles investigation of the excitonic properties of magnesia (MgO), an ionic insulator known to host hole polarons. We combine a density functional theory-based approach for structural relaxation in the presence of the hole and many-body perturbation theory to describe the excitonic properties. We determine that the hole polaron introduces new in-gap occupied states 0.6–0.8 eV above the valence band maximum that lead to two low-energy peaks in the optical spectrum. The predicted redshift of the lowest-energy transition due to polaron formation of 0.8 eV agrees well with the experimental Stokes shift of 0.8–0.9 eV. Analysis of the exciton wave function indicates that the electron-hole pair consists of a localized hole and delocalized electron, but that the wave function retains its Wannier-Mott character even in the presence of the hole polaron. Our study demonstrates that combining these previously established methods allows for a relatively computationally inexpensive approach to studying the exciton polaron in materials where only one charge carrier forms a polaron.

electronic structure↗
DOE OSTI2023

Improved quasiparticle self-consistent electronic band structure and excitons in β – LiGaO 2

The band structure of β–LiGaO 2 is calculated using the quasiparticle self-consistent QSGWˆ method where the screened Coulomb interaction Wˆ is evaluated including electron-hole interaction ladder diagrams and G is the one-electron Green's function. Improved convergence compared to previous calculations leads to a significantly larger band gap of about 7.0 eV. However, exciton binding energies are found to be large and lead to an exciton gap of about 6.0 eV if also a zero-point-motion correction of about –0.4 eV is included. Furthermore, these results are in excellent agreement with recent experimental results on the onset of absorption. Besides the excitons observed thus far, the calculations indicate the existence of a Rydberg-like series of exciton excited states, which is however modified from the classical Wannier exciton model by the anisotropies of the material and the more complex mixing of Bloch states in the excitons resulting from the Bethe-Salpeter equation. The exciton fine structure and the exciton wave functions are visualized and analyzed in various ways.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effect of off-diagonal elements in the Wannier Hamiltonian on DFT + DMFT for low-symmetry materials: Study of Li 2 MnO 3

Here, we study the effect of the off-diagonal elements of the Wannier Hamiltonian on the electronic structure of the low-symmetry material Li 2 MnO 3 ( C2/m ), using dynamical mean field theory calculations with a continuous-time quantum Monte Carlo impurity solver. The presence of significant off-diagonal elements leads to a pronounced suppression of the energy gap. The off-diagonal elements are largest when the Wannier projection is used based on the global coordinate, and they remain substantial even with the projection using the local coordinate close to the direction of Mn-O bonds. We show that the energy gap is enhanced by the diagonalization of the Mn d block in the full p-d Hamiltonian with the application of a unitary rotation matrix. Additionally, the inclusion of small double counting energy is crucial for achieving the experimental gap by reducing p-d hybridization. Furthermore, we establish the efficiency of a low-energy (d-only basis) model for studying the electronic structure of Li 2 MnO 3 , as the Wannier basis represents a hybridized state of Mn d and O p orbitals. These findings suggest an appropriate approach for investigating low-symmetry materials using the density functional theory plus dynamical mean field theory (DFT + DMFT) method. We also find that the antiferromagnetic ground state $\Gamma$ 2u is stable with U ≤ 2 eV within density functional theory+$U$ calculations, which is much smaller than the widely used U = 5 eV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗