Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “electronic structure theory”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

MISPR : an open-source package for high-throughput multiscale molecular simulations

Computational tools provide a unique opportunity to study and design optimal materials by enhancing our ability to comprehend the connections between their atomistic structure and functional properties. However, designing materials with tailored functionalities is complicated due to the necessity to integrate various computational-chemistry software (not necessarily compatible with one another), the heterogeneous nature of the generated data, and the need to explore vast chemical and parameter spaces. The latter is especially important to avoid bias in scattered data points-based models and derive statistical trends only accessible by systematic datasets. Here, we introduce a robust high-throughput multi-scale computational infrastructure coined MISPR (Materials Informatics for Structure–Property Relationships) that seamlessly integrates classical molecular dynamics (MD) simulations with density functional theory (DFT). By enabling high-performance data analytics and coupling between different methods and scales, MISPR addresses critical challenges arising from the needs of automated workflow management and data provenance recording. The major features of MISPR include automated DFT and MD simulations, error handling, derivation of molecular and ensemble properties, and creation of output databases that organize results from individual calculations to enable reproducibility and transparency. In this work, we describe fully automated DFT workflows implemented in MISPR to compute various properties such as nuclear magnetic resonance chemical shift, binding energy, bond dissociation energy, and redox potential with support for multiple methods such as electron transfer and proton-coupled electron transfer reactions. The infrastructure also enables the characterization of large-scale ensemble properties by providing MD workflows that calculate a wide range of structural and dynamical properties in liquid solutions. MISPR employs the methodologies of materials informatics to facilitate understanding and prediction of phenomenological structure–property relationships, which are crucial to designing novel optimal materials for numerous scientific applications and engineering technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Out-of-Plane Magnetic Anisotropy Engineered via Band Distortion in Two-Dimensional Materials

We unveil the connection between magnetic anisotropy and band structures in Γ-specialized In 2 Se 3 /Fe 3 GeTe 2 heterostructure, through first-principles calculations. A significant band distortion occurs near the Γ point with a formed valley. As the band-distortion-caused valley moves relative to the quasivalence band maximum, the magnetic anisotropy switches between in-plane and out-of-plane magnetic anisotropy. Such a rule applies to either the strained or polarization-switched In 2 Se 3 /Fe 3 GeTe 2 heterostructure. Furthermore, these findings demonstrate the feasibility of predicting the magnetic anisotropy by band structure modification in two-dimensional magnetic systems.

2-dimensional systems↗
DOE OSTI2023

Impact of the Electronic Properties of Chalcogenide Ligands in their Complexation with Uranyl Nitrate Complexes

The idea of covalency being directly correlated with selectivity in f-element separations has gained attention over the years. However, “covalency” as a concept is still under debate as it is not an experimental nor quantum mechanical observable. Regardless, covalency along with ionicity are still the main concepts used to explain the nature of the chemical bond of coordination complexes. Another concept that has gained the attention in the separations community is the hard-soft acid-base (HSAB) theory, which favors soft-donor ligands having the potential to increase selectivity for actinides (An) over the lanthanides (Ln) due to their ability to engage more effectively f-electrons in covalent interactions. However, it has been shown that while using softer donor ligands in fact increased orbital mixing and covalency with the actinides, it compromises the strength of the interaction owing to a decrease in the metal-ligand electrostatic interactions.1 Herein, we study the electronic structure of different types of chalcogenide ligands (LE, E = O, S, Se), their interaction with the uranyl ion, and their complexation free energies. Surprisingly, contrary to the expected order of covalency, i.e. An–LSe > An–LS > An–LO, our results suggest that the degree of covalency and associated complexation free energies depend on the nature of the coordinating moiety rather than solely on the directly coordinated atom. 1) Sadhu, Biswajit, and Michael Dolg. Enhancing actinide (III) over lanthanide (III) selectivity through hard-by-soft donor substitution: exploitation and implication of near-degeneracy-driven covalency. Inorganic Chemistry, 2019, 58, 9738 – 9748.

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

A nearsighted force-training approach to systematically generate training data for the machine learning of large atomic structures

A challenge of atomistic machine-learning (ML) methods is ensuring that the training data are suitable for the system being simulated, which is particularly challenging for systems with large numbers of atoms. Most atomistic ML approaches rely on the nearsightedness principle (“all chemistry is local”), using information about the position of an atom’s neighbors to predict a per-atom energy. Here, in this work, we develop a framework that exploits the nearsighted nature of ML models to systematically produce an appropriate training set for large structures. We use a per-atom uncertainty estimate to identify the most uncertain atoms and extract chunks centered around these atoms. It is crucial that these small chunks are both large enough to satisfy the ML’s nearsighted principle (that is, filling the cutoff radius) and are large enough to be converged with respect to the electronic structure calculation. We present data indicating when the electronic structure calculations are converged with respect to the structure size, which fundamentally limits the accuracy of any nearsighted ML calculator. These new atomic chunks are calculated in electronic structures, and crucially, only a single force—that of the central atom—is added to the growing training set, preventing the noisy and irrelevant information from the piece’s boundary from interfering with ML training. The resulting ML potentials are robust, despite requiring single-point calculations on only small reference structures and never seeing large training structures. We demonstrated our approach via structure optimization of a 260-atom structure and extended the approach to clusters with up to 1415 atoms.

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

Assessing Nontrivial Topology in Weyl Semimetals by Dichroic Photoemission

The electronic structure of Weyl semimetals features Berry flux monopoles in the bulk and Fermi arcs at the surface. While angle-resolved photoelectron spectroscopy (ARPES) is successfully used to map the bulk and surface bands, it remains a challenge to explicitly resolve and pinpoint these topological features. Here we combine state-of-the-art photoemission theory and experiments over a wide range of excitation energies for the Weyl semimetals TaAs and TaP. Further, our results show that simple surface-band-counting schemes, proposed previously to identify nonzero Chern numbers, are ambiguous due to pronounced momentum-dependent spectral weight variations and the pronounced surface-bulk hybridization. Instead, our findings indicate that dichroic ARPES provides an improved approach to identify Fermi arcs but requires an accurate description of the photoelectron final state.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Interpretable machine learning-guided design of Fe-based soft magnetic alloys

Here, we present a machine learning (ML) guided approach to predict saturation magnetization (𝑀 S ) and coercivity (𝐻 C ) in Fe-rich soft magnetic alloys, particularly Fe-Si-B systems. ML models trained on experimental data reveal that increasing Si and B content reduces 𝑀 S from 1.81 T (DFT ≈ 2.04 T) to ≈1.54 T (DFT ≈ 1.56T) in Fe-Si-B, which is attributed to decreased magnetic density and structural modifications. Experimental validation of ML predicted magnetic saturation on Fe-1Si-1B (2.09 T), Fe-5Si-5B (2.01 T), and Fe-10Si-10B (1.54 T) alloy compositions further supports our findings. These trends are consistent with density functional theory predictions, which link increased electronic disorder and band broadening to lower 𝑀 S values. Experimental validation on selected alloys confirms the predictive accuracy of the ML model, with good agreement across compositions. Beyond predictive accuracy, detailed uncertainty quantification and model interpretability including through feature importance and partial dependence analysis reveal that 𝑀 S is governed by a nonlinear interplay between Fe content and early transition metal ratios, while 𝐻 C is more sensitive to processing conditions such as ribbon thickness and thermal treatment windows. The ML framework was further applied to Fe-Si-B/Cr/Cu/Zr/Nb alloys in a pseudoquaternary compositional space, which shows comparable magnetic properties to NANOMET (Fe 84.8 ⁢Si 0.5 ⁢B 9.4 ⁢Cu 0.8⁢ P 3.5 ⁢C 1 ), FINEMET (Fe 73.5 ⁢Si 13.5 ⁢B 9 Cu 1 ⁢Nb 3 ), NANOPERM (Fe 88 ⁢Zr 7⁢ B 4 ⁢Cu 1 ), and HITPERM (Fe 44 ⁢Co 44 ⁢Zr 7⁢ B 4 ⁢Cu 1 . Our findings demonstrate the potential of the ML framework for accelerated search of high-performance soft magnetic materials.

density functional theory↗
DOE OSTI2025

Relationship between atomic and electronic structure in Ln-bearing oxides

Metastable states of matter are of great interest as they offer the promise of novel functionality. They are often a natural consequence of exposure to nonequilibrium environments. In oxides, metastability can take the form of new polymorphs, chemical disorder, and even amorphization. While significant attention has been given to the impact those changes have on the atomic properties of the material, the corresponding changes in the electronic structure have received less attention. Here, using density functional theory, we consider how the electronic structure varies with potential metastable structures in two classes of lanthanide-bearing oxides—pyrochlores and interlanthanide sesquioxides. We find that the changes depend strongly on both the crystal structure and crystal chemistry of the compound with, for example, disordering and amorphization either increasing or decreasing the bandgap depending on the chemistry. For the 𝐴 2 ⁢𝐵 2 ⁢O 7 pyrochlores, we find different dependencies of the bandgap on the 𝐴 = 𝐿⁢𝑛 cations as the 𝐵 cation is changed, which we relate to the nature of the density of states at the conduction band minimum for different 𝐵 chemistries. Our calculations are validated by electron energy loss spectroscopy measurements for two pyrochlore compounds in which amorphization does reduce the bandgap, consistent with our calculations on these two compounds. In conclusion, our results highlight the relationship between atomic and electronic structure and how radiation can be used to modify and potentially control the electronic properties of oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nernst power factor and figure of merit in the compensated semimetal ScSb

Recently, topological semimetals have emerged as strong candidates for solid-state thermomagnetic refrigerators due to their enhanced Nernst effect. This enhancement arises from the combined contributions of the Berry-curvature-induced anomalous Nernst coefficient associated with topological bands and the normal Nernst effect resulting from synergistic electron-hole compensation. Generally, these two effects are intertwined in topological semimetals, making it challenging to evaluate them independently. Here, we report the observation of a high Nernst effect in the electron-hole compensated semimetal ScSb with topologically trivial electronic band structures. Remarkably, we find a high maximum Nernst power factor of 𝑃⁢𝐹 𝑁 ∼ 35 × 10 −4 W m −1 K −2 in ScSb. The Nernst thermopower (𝑆 𝑥⁢𝑦 ) exhibits a peak of ∼ 47 µ⁢V/K at 12 K and 14 T, yielding a Nernst figure of merit (𝑧 𝑁 ) of ∼ 28 × 10 −4 K −1 . Notably, despite its trivial electronic band structure, both the 𝑃⁢𝐹 𝑁 and 𝑧 𝑁 values of ScSb are comparable to those observed in topological semimetals with Dirac band dispersions. In conclusion, the origin of the large Nernst signal in ScSb is explained by compensated electron and hole carriers, through Hall resistivity measurements, angle-resolved photoemission spectroscopy, and density functional theory calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Effects of strain and defects on the topological properties of HfTe 5

Manipulating topological states is a topic of vigorous research. However, the impact of such topological phase transitions on the quasiparticle dynamics remains elusive. Here, in this work, we present the effects of a transition from a strong to weak topological insulator in HfTe 5 as a function of Te vacancy concentration. We observed a significant transition from distinct sharp surface states and Dirac crossing to a Fermi-liquid-like quasiparticle state in which these surface-localized features are heavily suppressed. Additionally, by inducing the same effect through applied uniaxial stress, we demonstrate that changes in the lattice constants play the foremost role in determining the electronic structure, self-energy, and topological states of HfTe 5 . Our results demonstrate the possibility of using both defect chemistry and strain as control parameters for topological phase transitions and associated many-body physics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics

Extended Lagrangian Born–Oppenheimer molecular dynamics [ Phys. Rev. Lett. 2008, 100, 123004] is presented for Hartree–Fock theory, where the extended electronic degrees of freedom are represented by a density matrix, including fractional occupation numbers at elevated electronic temperatures. In contrast to regular direct Born–Oppenheimer molecular dynamics simulations, no iterative self-consistent field optimization is required prior to the force evaluations. To sample regions of the potential energy landscape where the gap is small or vanishing, which leads to particular convergence problems in regular direct Born–Oppenheimer molecular dynamics simulations, an adaptive integration scheme for the extended electronic degrees of freedom is presented. The integration scheme is based on a tunable, low-rank approximation of a fourth-order kernel, K, that determines the metric tensor, T ≡ K T K, used in the extended harmonic oscillator of the Lagrangian that generates the dynamics of the electronic degrees of freedom. Here, the formulation and algorithms provide a general guide to implement extended Lagrangian Born–Oppenheimer molecular dynamics for quantum chemistry, density functional theory, and semiempirical methods using a density matrix formalism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2023

Influence of Cation Species on Thermal Expansion of Y2Si2O7–Gd2Si2O7 Solid Solutions

Mixtures of Y 2 Si 2 O 7 and Gd 2 Si 2 O 7 were synthesized by solid-state reaction at 1600°C and characterized via in situ x-ray diffraction (XRD) to determine their coefficients of thermal expansion (CTE). All solid solutions within the system exhibited the orthorhombic δ-RE 2 Si 2 O 7 (Pna2 1 ) structure. Thermal expansion measurements of Y 2 Si 2 O 7 and Gd 2 Si 2 O 7 correlated well with reported values in literature, and all synthesized solid solutions exhibited CTEs between Y 2 Si 2 O 7 and Gd 2 Si 2 O 7 . Generally, there was a slight decrease in CTE exhibited by the materials with increasing Gd 2 Si 2 O 7 content, with Gd 2 Si 2 O 7 having the lowest CTEs and Y 2 Si 2 O 7 the highest CTEs. The decrease in CTE was attributed to stronger bonds of Gd-O over Y-O, as determined by calculated crystal orbital Hamilton populations using density functional theory. However, such differences were very small and crystal structure was the dominating factor in CTE trends.

rare earth silicates↗
DOE OSTI2026

Disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y alloys

This dataset contains input and output files from density functional theory (DFT) simulations used to study the disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y (0 ≤ x + y ≤ 1) alloys and their binary end members MnTe, MnSb, and MnBi. The alloys adopt the hexagonal NiAs-type (nickeline) structure and span ternary (MnTexSb1-x, MnTexBi1-x, MnBixSb1-x), and quaternary compositions across the full MnTe–MnSb–MnBi composition triangle. For each alloy composition, the dataset provides DFT calculations in three magnetic configurations: A-type antiferromagnetic (AFM), C-type AFM, and ferromagnetic (FM). Every magnetic configuration folder contains the fully relaxed crystal structure (CONTCAR), VASP input parameters (INCAR), and the main VASP output file (OUTCAR), from which total electronic energies, Mn magnetic moments, lattice parameters, and percent volume changes between magnetic states are extracted. These data are used to construct compositional phase diagrams, evaluate thermodynamic stability (formability), and map magnetoelastic responses across the alloy space. For A-type AFM and FM configurations, additional data are provided as follows: (i) FORCE_CONSTANTS and thermal_properties.yaml files at the top level of A-type_AFM/ and FM/ folders — present only for compositions marked with an asterisk (*) in Table I of the main text. These are derived from Phonopy finite-displacement calculations on full disordered 128-atom supercells and provide vibrational free energy, entropy (Svib)contribution from explicit disorder calculations. (Table I of the associated main manuscript) (ii) A VCA/ subfolder within A-type_AFM/ and FM/, containing FORCE_CONSTANTS and thermal_properties.yaml from Virtual Crystal Approximation phonon calculations (without spin-orbit coupling). VCA data are available for all compositions and are used to estimate vibrational contributions to the Gibbs free energy across the full composition space. (iii) A SOC/ subfolder containing CONTCAR, INCAR, and OUTCAR from spin-orbit coupling calculations, providing relativistic corrections to electronic energies and lattice parameters (Tables S2–S3 of the SM, and Table I of the main manuscript). (iv) A SOC/VCA/ subfolder containing FORCE_CONSTANTS and thermal_properties.yaml from VCA phonon calculations performed within the SOC framework, combining relativistic and vibrational thermodynamic corrections. The computed properties are used to map the AFM–FM magnetic crossover near MnTe0.75Sb0.25, demonstrate disorder- and spin-induced phonon broadening, identify a semiconductor-to-metal crossover, and quantify the pronounced magnetoelastic volume response near the magnetic phase boundary.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Substrate-Dependent Band Structures in Trilayer Graphene / h - BN Heterostructures

The tight-binding model has been spectacularly successful in elucidating the electronic and optical properties of a vast number of materials. Within the tight-binding model, the hopping parameters that determine much of the band structure are often taken as constants. Here, using ABA-stacked trilayer graphene as the model system, we show that, contrary to conventional wisdom, the hopping parameters and therefore band structures are not constants, but are systematically variable depending on their relative alignment angle between h-BN. Moreover, the addition or removal of the h-BN substrate results in an inversion of the K and K' valley in trilayer graphene’s lowest Landau level. Our work illustrates the oft-ignored and rather surprising impact of the substrates on band structures of 2D materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

One-dimensionality signature in optical conductivity of heavy-fermion CeIr 3 B 2

In low dimensions, the combined effects of interactions and quantum fluctuations can lead to dramatically new physics distinct from that existing in higher dimensions. Here, we investigate the electronic and optical properties of CeIr 3 ⁢B 2 , a quasione-dimensional (1D) Kondo lattice system, using ab initio calculations. The Ce atoms in the hexagonal crystal structure form 1D chains along the c axis, with extremely short Ce-Ce distances. The quasi-1D nature of the crystal structure is well reflected in its electronic structure. Extremely flat bands emerge within the ab plane of the Brillouin zone, yielding sharp optical transitions in the corresponding optical conductivity. Further, our calculations indicate that these prominent peaks in the optical conductivity provide a clear signature of quasi-1D heavy fermion systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Hourglass-type bulk Ni 3 d band and Ce 4 f Kondo resonance states in the potential topological Kondo semimetal CeNiSn via angle-resolved photoemission spectroscopy

The electronic structure of CeNiSn, which is a potential topological Kondo insulator and a Dirac nodal-loop semimetal, has been investigated by employing temperature (T) dependent angle-resolved photoemission spectroscopy (ARPES). The Fermi surfaces (FSs) and the band structures of CeNiSn for three orthogonal crystallographic planes are measured, in which both the very dispersive bands and the flat bands are observed, having mainly the Ni 3d character and the Ce 4f character, respectively. The measured FSs and ARPES bands agree reasonably well with the density functional theory (DFT) calculations. The Fermi-edge (E F ) photon energy (hν) map along k b (=k (010) ) shows that the metallic E F -crossing states on the (010) surface have the three-dimensional character, suggesting that the observed E F -crossing metallic states do not correspond to the topological surface states of the two-dimensional character. On the other hand, albeit weak, the features of the hourglass-type bulk band crossings are observed along SXS, with the energies and the slopes being similar to those predicted by the DFT calculations, supporting the Dirac semimetallic nature of CeNiSn. In T-dependent ARPES, the Ce 4f Kondo resonance states are clearly revealed at low T, which become much suppressed above ~80 K. This feature is consistent with the Kondo temperature of CeNiSn, estimated from its ρ(T) data. This work demonstrates the importance of the coherent Kondo states in determining the topological properties of CeNiSn.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Performance of Heterostructural TaC/AlGaN Schottky Diodes Based on First Principles Electronic Structure Properties

Advances in ultra-wide bandgap materials, such as high Al-content AlxGa1-xN (AlGaN), are essential for next generation power electronics, but the requirement for lattice matched substrates is currently a significant obstacle. Recently, conductive TaC has emerged as a promising virtual substrate for AlGaN heteroepitaxy, with wurtzite (0001) AlxGa1-xN lattice-matched to rocksalt (111) TaC at x ~ 0.5. Thus, understanding and controlling the electronic properties of the TaC/AlGaN interface is key for developing technological applications based on TaC/AlGaN devices. Using density functional theory and electronic structure calculations, we here investigate TaC/Al0.5Ga0.5N interfaces, where we include explicit alloy models in the slab calculations. We predict the Schottky barrier height and the electric field discontinuity resulting from interface charges. Considering all possible combinations of (Ta, C) substrate termination, (Al/Ga, N) nucleation, and (Al/Ga, N) polarity, we construct a chemical potential phase diagram to identify the stable interfaces that can be accessed through variation of the synthesis conditions. The predicted interface electronic properties are implemented in device performance simulations to demonstrate a practical design for a strain-free, high-efficiency TaC/AlGaN Schottky diode with a low barrier height and without interface charges, underscoring the potential of TaC as a substrate for ultra-wide bandgap devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Electronic structure and two-band superconductivity in unconventional high- T c cuprates Ba 2 CuO 3+δ

We report that the recently discovered cuprate superconductor Ba 2 CuO 3+δ exhibits a high T c ≃73 K at δ ≃ 0.2. The polycrystal grown under high pressure has a structure similar to La 2 CuO 4 but with dramatically different lattice parameters due to the CuO 6 octahedron compression. The crystal field in the compressed Ba 2 CuO 4 leads to an inverted Cu 3d e g complex with the d x 2 -y 2 orbital sitting below the d 3z 2 - r 2 and an electronic structure highly unusual compared to the conventional cuprates. We construct a two-orbital Hubbard model for the Cu d 9 state at hole doping x = 2δ and study the orbital-dependent strong correlation and superconductivity. For the undoped case at x = 0 , we found that strong correlation drives an orbital-polarized Mott-insulating state with the spin-1/2 moment of the localized d 3z 2 - r 2 orbital. In contrast to the single-band cuprates where superconductivity is suppressed in the overdoped regime, hole doping the two-orbital Mott insulator leads to orbital-dependent correlations and the robust spin and orbital exchange interactions produce a high-T c antiphase d-wave superconductor even in the heavily doped regime at x = 0.4 . We conjecture that Ba 2 CuO 3+δ realizes mixtures of such heavily hole-doped superconducting Ba 2 CuO 4 and disordered Ba 2 CuO 3 chains in a single-layer or predominately separated bilayer structure. Our findings suggest that unconventional cuprates with liberated orbitals as doped two-band Mott insulators can be a direction for realizing high-T c superconductivity with enhanced transition temperature T c .

36 MATERIALS SCIENCE↗
DOE OSTI2023

Shadow energy functionals and potentials in Born–Oppenheimer molecular dynamics

In Born–Oppenheimer molecular dynamics (BOMD) simulations based on the density functional theory (DFT), the potential energy and the interatomic forces are calculated from an electronic ground state density that is determined by an iterative self-consistent field optimization procedure, which, in practice, never is fully converged. The calculated energies and forces are, therefore, only approximate, which may lead to an unphysical energy drift and instabilities. Here, we discuss an alternative shadow BOMD approach that is based on backward error analysis. Instead of calculating approximate solutions for an underlying exact regular Born–Oppenheimer potential, we do the opposite. Instead, we calculate the exact electron density, energies, and forces, but for an underlying approximate shadow Born–Oppenheimer potential energy surface. In this way, the calculated forces are conservative with respect to the approximate shadow potential and generate accurate molecular trajectories with long-term energy stabilities. We show how such shadow Born–Oppenheimer potentials can be constructed at different levels of accuracy as a function of the integration time step, δt, from the constrained minimization of a sequence of systematically improvable, but approximate, shadow energy density functionals. For each energy functional, there is a corresponding ground state Born–Oppenheimer potential. These pairs of shadow energy functionals and potentials are higher-level generalizations of the original “zeroth-level” shadow energy functionals and potentials used in extended Lagrangian BOMD. The proposed shadow energy functionals and potentials are useful only within this extended dynamical framework, where also the electronic degrees of freedom are propagated as dynamical field variables together with the atomic positions and velocities. The theory is quite general and can be applied to MD simulations using approximate DFT, Hartree–Fock, or semi-empirical methods, as well as to coarse-grained flexible charge models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗