Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles calculation”

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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Fermi surface reconstruction due to orthorhombic distortion in the Dirac semimetal YbMnSb 2

A Dirac semimetal with magnetic atoms as constituents delivers an interesting platform to investigate the interplay of Fermi surface (FS) topology, electron correlation, and magnetism. One such family of semimetals is YbMn⁢ 𝑃⁢𝑛 2 (𝑃⁢𝑛 = Sb,Bi), which is being actively studied due to the intertwined spin and charge degrees of freedom. Here, in this study, we investigate the relationship between the magnetic/crystal structures and FS topology of YbMnSb 2 using single crystal x-ray diffraction, neutron scattering, magnetic susceptibility, magnetotransport measurements, and complimentary density functional theory (DFT) calculations. Contrary to previous reports, the x-ray and neutron diffraction reveal that YbMnSb 2 crystallizes in an orthorhombic 𝑃⁢𝑛⁢𝑚⁢𝑎 structure with a notable antiphase displacement of the magnetic Mn ions that increases in magnitude upon cooling. As a result of this orthorhombicity, first-principles DFT calculation reveals a reduced Brillouin zone and more anisotropic FS of YbMnSb 2 compared to YbMnBi 2 . Moreover, the hole-type carrier density drops by two orders of magnitude as YbMnSb 2 orders antiferromagnetically, indicating band folding in a magnetic ordered state. In addition, the Landau level fan diagram yields a nontrivial nature of the Shubnikov–de Haas quantum oscillation frequency arising from the Dirac-like Fermi pocket. These results imply that YbMnSb 2 is an ideal platform to explore the interplay of subtle lattice distortion, magnetic order, and topological transport arising from relativistic quasiparticles.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Relativistic self-consistent G W : Exact two-component formalism with one-electron approximation for solids

We present a formulation of relativistic self-consistent GW for solids based on the exact two-component formalism with one-electron approximation (X2C1e) and nonrelativistic Coulomb interactions. Our theory allows us to study scalar relativistic effects, spin-orbit coupling, and the interplay of relativistic effects with electron correlation without adjustable parameters. Our all-electron implementation is fully ab initio and does not require a pseudopotential constructed from atomic calculations. We examine the effect of the X2C1e approximation by comparison to the established four-component formalism and reach an excellent agreement. The simplicity of X2C1e enables the construction of higher order theories, such as embedding theories, on top of perturbative calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Machine-learned quantum molecular dynamics calculations of warm dense equation of state and ionic transport coefficients of deuterated water

White dwarf models require accurate equations of state and ionic transport coefficients in the warm dense matter regime, where kinetic theory models and tabulated equations of state are often inaccurate. In this work, spectral-partitioned density functional theory and machine-learned interatomic potentials are combined to perform large-scale, first-principles quantum molecular dynamics simulations of deuterated water (D 2 O) near the principal Hugoniot. This approach retains Kohn-Sham accuracy while achieving orders-of-magnitude speedup, yielding converged equation of state and transport properties over a broad pressure and temperature range. The results reveal the thermodynamic conditions under which ionic transport models for interdiffusivity and shear viscosity converge and identify those in closest agreement with density functional theory benchmarks at temperatures in the warm dense matter regime. The present framework extends first-principles transport calculations to higher temperatures than previously achieved, and provides an efficient, scalable, and general approach for studying transport properties in complex multicomponent mixtures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Quantum computation of stopping power for inertial fusion target design

Stopping power is the rate at which a material absorbs the kinetic energy of a charged particle passing through it—one of many properties needed over a wide range of thermodynamic conditions in modeling inertial fusion implosions. First-principles stopping calculations are classically challenging because they involve the dynamics of large electronic systems far from equilibrium, with accuracies that are particularly difficult to constrain and assess in the warm-dense conditions preceding ignition. Here, we describe a protocol for using a fault-tolerant quantum computer to calculate stopping power from a first-quantized representation of the electrons and projectile. Our approach builds upon the electronic structure block encodings of Su et al. [ PRX Quant. 2 , 040332 (2021)], adapting and optimizing those algorithms to estimate observables of interest from the non-Born–Oppenheimer dynamics of multiple particle species at finite temperature. We also work out the constant factors associated with an implementation of a high-order Trotter approach to simulating a grid representation of these systems. Ultimately, we report logical qubit requirements and leading-order Toffoli costs for computing the stopping power of various projectile/target combinations relevant to interpreting and designing inertial fusion experiments. We estimate that scientifically interesting and classically intractable stopping power calculations can be quantum simulated with roughly the same number of logical qubits and about one hundred times more Toffoli gates than is required for state-of-the-art quantum simulations of industrially relevant molecules such as FeMoco or P450.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Influence of crystalline order and defects on the absolute work functions and electron affinities of TiO2- and SrO-terminated n-SrTiO3(001)

We establish how the detailed structure and composition of the terminal layer of n-SrTiO3 (001) determine key surface electronic properties. We have measured and calculated from first principles the absolute work functions and electron affinities of bulk SrNb0.01Ti0.99O3(001) with both TiO2 and nominally SrO terminal planes. The match between theory and experiment is excellent at 4.6(1) eV for the TiO2 termination if an ideal, bulk-truncated surface structure is assumed. In contrast, the ideal SrO-termination leads to a calculated work function that is ~1.3 eV lower than the experimental value of 3.6(1) eV. We show that this discrepancy is associated with defects on the SrO surface that serve as electron scavengers. These defects deplete the concentration of itinerant electrons in the subsurface region and increase the negative charge density on the surface, thus increasing the work function. Several different surface defect configurations were modeled; the ones that yield the best agreement with experiment involve Sr vacancies in the terminal layer along with OH–/H– pairs or O22– occupying anion sites adjacent to the Sr vacancies.

Chambers, Scott A.↗
DOE OSTI2021

Prediction of unconventional magnetism in doped FeSb 2

Significance For many decades, it has been commonly believed that all electronic states of a collinear antiferromagnet (AF) are spin-degenerate, unless the underlying crystal structure lacks centrosymmetry and has spin–orbit coupling. This has been essentially definitional for antiferromagnetism and is widely used experimentally to distinguish ferromagnets from AFs. Recently, it was demonstrated that a new class of magnets, possessing antiferromagnetic order and without net magnetization but showing a typical ferromagnetic response in many aspects, is possible. We predict that F e S b 2 , which is well known but poorly understood magnetically, is an incipient unconventional magnet of this type and can be pushed to become one by Co or Cr doping. Moreover, the calculated magnetic anisotropy is favorable for exhibiting various anomalous properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

An energy-modified quantum defect method for the analysis of Rydberg spectra: Application to 2-butyne

The high resolution Rydberg absorption spectrum of 2-butyne C 4 H 6 recorded previously at the SOLEIL synchrotron facility has been interpreted using multichannel quantum defect theory (MQDT). The calculations are based on the continuum scattering calculations of Xu et al., J. Chem. Phys. 136, 154303 (2012) and of Jacovella et al., J. Phys. Chem. A 119, 12339 (2015) pertaining to the dipole-allowed excited state symmetries in absorption from the ground state. In contrast to the traditional approach of calculating low-lying electronic states first and then attempting to extend the calculations to ever higher energy, here the analysis proceeds through the extension of these detailed calculations of the electronic continuum scattering down into the discrete region of the spectrum. The continuum reaction matrices and dipole transition moments are adapted to the discrete Rydberg region via the use of an energy-modified formulation of MQDT theory and associated energy dependences of the quantum defects. The analysis reproduces more than 40 Rydberg states from n ≈ 10 down to the 3 d and 4 s levels with an rms error of better than 20 cm −1 . Further, these belong to five Rydberg series with three different molecular symmetries. While the approach predicts many additional series, most of these are calculated and observed to carry only little oscillator strength. The analysis shows that the Rydberg spectrum is dominated by the excitation of an e″ symmetry electron of fδ and gπ type, in line with what previous studies of the above-threshold shape resonance of 2-butyne have shown. The present study is intended to serve as an example showing how first principles continuum calculations may be useful for the interpretation of highly bound discrete states in a range that poses problems for the standard ab initio techniques. The quantitative treatment of the dipole absorption cross sections is deferred to a future paper.

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

Full-potential KKR within the removed-sphere method: A practical and accurate solution to the Poisson equation

Here an efficient and accurate generalization of the removed-sphere method (RSM) to solve the Poisson equation for total charge density in a solid with space-filling convex Voronoi polyhedra (VPs) and any symmetry is presented. The generalized RSM avoids the use of multipoles and VP shape functions for cellular integrals, which have associated ill-convergent large, double-internal L sums in spherical-harmonic expansions, so that fast convergence in single-L sums is reached. Our RSM adopts full Ewald formulation to work for all configurations or when symmetry breaking occurs, such as for atomic displacements or elastic constant calculations. The structure-dependent coefficients A L that define RSM can be calculated once for a fixed structure and speed up the whole self-consistent-field procedure. The accuracy and rapid convergence properties are confirmed using two analytic models, including the Coulomb potential and energy. We then implement the full-potential RSM using the Green's function Korringa-Kohn-Rostoker (KKR) method for real applications and compare the results with other first-principle methods and experimental data, showing that they are equally as accurate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Starting-point-independent quantum Monte Carlo calculations of iron oxide

Quantum Monte Carlo (QMC) methods are useful for studies of strongly correlated materials because they are many body in nature and use the physical Hamiltonian. Typical calculations assume as a starting point a wave function constructed from single-particle orbitals obtained from one-body methods, e.g., density functional theory. However, mean-field-derived wave functions can sometimes lead to systematic QMC biases if the mean-field result poorly describes the true ground state. In this study, we examine the accuracy and flexibility of QMC trial wave functions using variational and fixed-node diffusion QMC estimates of the total spin density and lattice distortion of antiferromagnetic iron oxide (FeO) in the ground state B1 crystal structure. We found that for relatively simple wave functions the predicted lattice distortion was controlled by the choice of single-particle orbitals used to construct the wave function, rather than by subsequent wave function optimization techniques within QMC. By optimizing the orbitals with QMC, we then demonstrate starting-point independence of the trial wave function with respect to the method by which the orbitals were constructed by demonstrating convergence of the energy, spin density, and predicted lattice distortion for two qualitatively different sets of orbitals. The results suggest that orbital optimization is a promising method for accurate many-body calculations of strongly correlated condensed phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Dynamic mode decomposition for extrapolating nonequilibrium Green's-function dynamics

The Hartree-Fock generalized Kadanoff-Baym ansatz (HF-GKBA) offers an approximate numerical procedure for propagating the two-time nonequilibrium Green's function (NEGF). Here, using the GW self-energy, we compare the HF-GKBA to exact results for a variety of systems with long- and short-range interactions, different two-body interaction strengths, and various nonequilibrium preparations. We find excellent agreement between the HF-GKBA and exact time evolution in models when more realistic long-range exponentially decaying interactions are considered. Furthermore, this agreement persists for long times and for intermediate to strong interaction strengths. In large systems, HF-GKBA becomes prohibitively expensive for long-time evolutions. For this reason, we look at the use of dynamical mode decomposition (DMD) to reconstruct long-time NEGF trajectories from a sample of the initial trajectory. Using no more than 16% of the total time evolution, we reconstruct the total trajectory with high fidelity. Our results show the potential for DMD to be used in conjunction with HF-GKBA to calculate long-time trajectories in large-scale systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Metallization and spin fluctuations in Cu-doped lead apatite

An electronic structure and magnetic properties analysis of the recently proposed Cu-doped lead apatite is performed. Here, we show that electronic structures of differently Cu-substituted structures are characterized by localized molecular Cu-O bands at or near the Fermi level. The Cu substitutions can happen at both Pb1 and Pb2 sites, leading to metallic and semiconducting systems differently. The electronic systems in these bands are highly unstable magnetically and form clusters of rigidly ferromagnetically coupled magnetic moments on Cu and neighboring oxygen atoms with a total moment of about 1µ B . The ground state of uniformly Cu-doped lead apatite appears to be magnetic and semiconducting. The nonuniform distribution of two Cu atoms at the nearest Pb2 sites leads to an antiferromagnetic semiconducting state with formation energy close to uniformly distributed Cu configurations. The inclusion of quantum spin fluctuations confirms the stability of magnetic Cu-O clusters. Our calculations revealed the absence of the long-range magnetic order between uniformly distributed Cu-O clusters, creating the spin glass type of system.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Solving the electronic structure problem for over 100000 atoms in real space

Using a real-space high-order finite-difference approach, we investigate the electronic structure of large spherical silicon nanoclusters. Within Kohn-Sham density functional theory and using pseudopotentials, we report the self-consistent field convergence of a system with over 100000 atoms: a Si 107,641 ⁢H 9,084 nanocluster with a diameter of 16 nm. Our approach uses Chebyshev-filtered subspace iteration to speed up the convergence of the eigenspace, and blockwise Hilbert space-filling curves to speed up sparse matrix-vector multiplications, all of which are implemented in the parsec code. For the largest system, we utilized 2048 nodes (114 688 cores) on the Frontera machine in the Texas Advanced Computing Center. Our quantitative analysis of the electronic structure shows how it gradually approaches its bulk counterpart as a function of nanocluster size. The band gap is enlarged due to quantum confinement in nanoclusters, but decreases as the system size increases, as expected. In conclusion, our work serves as a proof of concept for the capacity of the real-space approach in efficiently parallelizing very large calculations using high-performance computer platforms, which can straightforwardly be replicated in other systems with more than 10 5 atoms.

0-dimensional systems↗
DOE OSTI2024

Charge self-consistent density functional theory plus ghost rotationally invariant slave-boson theory for correlated materials

We present a charge self-consistent density functional theory combined with the ghost rotationally invariant slave-boson (DFT+gRISB) formalism for studying correlated materials. Here, this method is applied to SrVO 3 and NiO, representing prototypical correlated metals and charge-transfer insulators. For SrVO 3 , we demonstrate that DFT+gRISB yields an accurate equilibrium volume and effective mass close to experimentally observed values. Regarding NiO, DFT+gRISB enables the simultaneous description of charge-transfer and Mott-Hubbard bands, significantly enhancing the accuracy of the original DFT+RISB approach. Furthermore, the calculated equilibrium volume and spectral function reasonably agree with experimental observations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Intrinsic thermal conductivity of ZrC from low to ultrahigh temperatures: A critical revisit

Current phonon transport theory based on ground-state calculations has been successful in predicting thermal conductivity at room and medium temperatures but may misrepresent behavior at high temperatures. Here, in this work, we predict the thermal conductivity (κ) of ZrC including electronic and phonon contributions from 300 to 3500 K, by including high-order phonon scattering; lattice expansion; temperature-dependent (TD) second-, third-, and fourth-order force constants (2FC, 3FC, and 4FC); and interband phonon conduction by using first principles. For the phonon transport, we find that four-phonon scattering (4ph) significantly reduces the phonon thermal conductivity (κ ph ), by as much as ~75% at 3500 K. After including 4ph scattering and all other factors, κ ph shows a ~T –1.5 rather than ~T –1 dependence. TD 2FC decreases three-phonon scattering rates but increases 4ph rates by decreasing and increasing the scattering phase spaces, respectively. For 4ph phase space, the TD 2FC flattens phonon bands, and allows more redistribution-4ph processes (1 + 2 → 3 + 4) to happen. The combination effect of TD 2FC and TD 4FC reduces 4ph rates of acoustic modes but increases those of optical modes. The TD 3FC and 4FC decrease the phonon scattering cross section and increase the κ ph significantly (by 52% at 3500 K). The contribution from interband (Wigner) phonon conduction is small, even at ultrahigh temperatures. For electronic thermal transport, we find that it is sensitive to and can be changed by 20% by the TD lattice constants. The Lorenz number varies from 1.6 to 3.3 × 10 –8 WΩK –2 at different temperatures. The theoretical prediction in the literature overpredicts κph (e.g., ~28%) and underpredicts the κ el (e.g., ~38%), resulting in an overall underprediction of κ (~26% at 1500 K). The impacts of grain size and defects are found to be strong, leading to the lower observed thermal conductivity in experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Asymmetric equilibrium core structures of pyramidal-II < c + a > dislocations in ten hexagonal-close-packed metals

The structures of pyramidal-II < c + a > dislocations, one of the most important defects in structural hexagonal-close-packed (HCP) metals, have not been fully characterized for many of the HCP metals in use today. Here, we employ ab initio informed phase-field dislocation dynamics to determine the minimum energy structure of pyramidal {1¯1¯22} < 11¯23 > dislocations in ten HCP metals, including Be, Co, Mg, Re, Ti, Zn, Cd, Hf, Y, and Zr. As input for the simulations, we calculate, using first-principles density functional theory, the {1¯1¯22} generalized stacking fault energy (GSFE) curves for all ten metals. From these calculations, it is found that magnetism in Co is necessary for achieving a local minimum in the GSFE curve. We observe in simulations that edge and screw character dislocations split into two partials separated by a low-energy intrinsic stacking fault. The splitting distance is shown to scale inversely with the local minimum energy normalized by the product of its shear modulus and Burgers vector. Interestingly, some HCP metals exhibit an asymmetric structure, with either unequal partial Burgers vectors or widths, in contrast to the symmetric configuration expected from linear elastic dislocation theory. We explain these structures by properties of the local maxima in their GSFE curves. Metals with larger degrees of elastic anisotropy result in dislocations with larger splitting distances than would be expected under the commonly used assumption of elastic isotropy. Furthermore, these findings on the sizes and asymmetry in the structures of pyramidal-II < c + a > dislocations are fundamental to understanding how these dislocations glide and interact or react with other defects when these metals are mechanically strained.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Unified Approach to Polarons and Phonon-Induced Band Structure Renormalization

Ab initio calculations of the phonon-induced band structure renormalization are currently based on the perturbative Allen-Heine theory and its many-body generalizations. These approaches are unsuitable to describe materials where electrons form localized polarons. Here, we develop a self-consistent, many-body Green’s function theory of band structure renormalization that incorporates localization and self-trapping. Here, we show that the present approach reduces to the Allen-Heine theory in the weak-coupling limit, and to total energy calculations of self-trapped polarons in the strong-coupling limit. To demonstrate this methodology, we reproduce the path-integral results of Feynman and diagrammatic Monte Carlo calculations for the Fröhlich model at all couplings, and we calculate the zero point renormalization of the band gap of an ionic insulator including polaronic effects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Changes in polarization dictate necessary approximations for modeling electronic deexcitation intensity: Application to x-ray emission

Accurate simulation of electronic excitations and deexcitations are critical for complementing complex spectroscopic experiments and can provide validation to theoretical approaches. Here, using a generalized framework, we contrast the accuracy and validity of orbital-constrained and linear-response approaches that build upon Kohn-Sham density functional theory (DFT) to simulate emission spectra of electronic origin and propose an efficient approximation, named many-body x-ray emission spectroscopy or MBXES, for simulating such processes. We show analytically as well as with computed examples that for electronic (de)excitation leading to an appreciable change in polarization (i.e., density rearrangement), the adiabatic approximation in a response-based formalism will be inadequate for the calculation of oscillator strength. Thus, such a change (e.g., in the net electrostatic dipole moment of a finite system) can be used as a metric for evaluating the applicability of the adiabatic response-based approach and can be particularly valuable in x-ray emission spectroscopy. On the other hand, MBXES, the flexible method introduced in this paper, can compute oscillator strengths accurately at a much lower computational expense on the basis of two DFT-based self-consistent field calculations. Using illustrative examples of emission spectra, the efficacy of the MBXES method is demonstrated by comparison with its parent theory, orbital-optimized DFT, and with experiments.

74 ATOMIC AND MOLECULAR PHYSICS↗