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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Toward improved property prediction of 2D materials using many-body quantum Monte Carlo methods

The field of 2D materials has grown dramatically in the past two decades. 2D materials can be utilized for a variety of next-generation optoelectronic, spintronic, clean energy, and quantum computing applications. These 2D structures, which are often exfoliated from layered van der Waals materials, possess highly inhomogeneous electron densities and can possess short- and long-range electron correlations. The complexities of 2D materials make them challenging to study with standard mean-field electronic structure methods such as density functional theory (DFT), which relies on approximations for the unknown exchange-correlation functional. To overcome the limitations of DFT, highly accurate many-body electronic structure approaches such as diffusion Monte Carlo (DMC) can be utilized. In the past decade, DMC has been used to calculate accurate magnetic, electronic, excitonic, and topological properties in addition to accurately capturing interlayer interactions and cohesion and adsorption energetics of 2D materials. Here, this approach has been applied to 2D systems of wide interest, including graphene, phosphorene, MoS 2 , CrI 3 , VSe 2 , GaSe, GeSe, borophene, and several others. In this review article, we highlight some successful recent applications of DMC to 2D systems for improved property predictions beyond standard DFT.

2D materials↗
DOE OSTI2024

Resonant inelastic x-ray scattering in the topological semimetal FeSi

The energy spectrum of topological semimetals contains protected degeneracies in reciprocal space that correspond to Weyl, Dirac, or multifold fermionic states. To exploit the unconventional properties of these states, one has to access the electronic structure of the three-dimensional bulk. In this paper, we present a joint theory-experiment study of the electronic structure of a candidate topological semimetal with resonant inelastic x-ray scattering (RIXS). We resolve the bulk electronic states of FeSi using momentum-dependent RIXS at the Fe L 3 edge. We observe a broad excitation continuum devoid of sharp features, consistent with particle-hole scattering in an underlying electronic band structure. Using density functional theory (DFT), we calculate the electronic structure of FeSi and derive a band theory formulation of RIXS in the fast collision approximation to model the scattering process with zero adjustable parameters. While band theory predicts an excitation continuum with broad spectral features similar to the observed ones, discrepancies between theory and experiment suggest the presence of low-energy processes that DFT alone does not account for. Further, this study of RIXS in a topological semimetal shows that RIXS is a useful tool for revealing unanticipated behavior of bulk electronic states in this class of materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

High-throughput ab initio reaction mechanism exploration in the cloud with automated multi-reference validation

Quantum chemical calculations on atomistic systems have evolved into a standard approach to studying molecular matter. These calculations often involve a significant amount of manual input and expertise, although most of this effort could be automated, which would alleviate the need for expertise in software and hardware accessibility. Here, we present the AutoRXN workflow, an automated workflow for exploratory high-throughput electronic structure calculations of molecular systems, in which (i) density functional theory methods are exploited to deliver minimum and transition-state structures and corresponding energies and properties, (ii) coupled cluster calculations are then launched for optimized structures to provide more accurate energy and property estimates, and (iii) multi-reference diagnostics are evaluated to back check the coupled cluster results and subject them to automated multi-configurational calculations for potential multi-configurational cases. All calculations are carried out in a cloud environment and support massive computational campaigns. Key features of all components of the AutoRXN workflow are autonomy, stability, and minimum operator interference. We highlight the AutoRXN workflow with the example of an autonomous reaction mechanism exploration of the mode of action of a homogeneous catalyst for the asymmetric reduction of ketones.

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

Electron avalanche structure determined by random walk theory

A self-consistent avalanche solution which accounts for collective long range Coulomb interactions as well as short range elastic and inelastic collisions between electrons and background atoms is made possible by a random walk technique. Results show that the electric field patterns in the early formation stages of avalanches in helium are close to those obtained from theory based on constant transport coefficients. Regions of maximum and minimum induced electrostatic potential phi are located on the axis of symmetry and within the volume covered by the electron swarm. As formation time continues, however, the region of minimum phi moves to slightly higher radii and the electric field between the extrema becomes somewhat erratic. In the intermediate formation periods the avalanche growth is slightly retarded by the high concentration of ions in the tail which oppose the external electric field. Eventually the formation of ions and electrons in the localized regions of high field strength more than offset this effect causing a very abrupt increase in avalanche growth.

Englert, G. W.↗
DOE OSTI2023

High-order finite element method for atomic structure calculations

We introduce featom, an open source code that implements a high-order finite element solver for the radial Schrödinger, Dirac, and Kohn-Sham equations. The formulation accommodates various mesh types, such as uniform or exponential, and the convergence can be systematically controlled by increasing the number and/or polynomial order of the finite element basis functions. The Dirac equation is solved using a squared Hamiltonian approach to eliminate spurious states. Here, to address the slow convergence of the $κ=±1$ states due to divergent derivatives at the origin, we incorporate known asymptotic forms into the solutions. We achieve a high level of accuracy (10 -8 Hartree) for total energies and eigenvalues of heavy atoms such as uranium in both Schrödinger and Dirac Kohn-Sham solutions. We provide detailed convergence studies and computational parameters required to attain commonly required accuracies. Finally, we compare our results with known analytic results as well as the results of other methods. In particular, we calculate benchmark results for atomic numbers (Z) from 1 to 92, verifying current benchmarks. We demonstrate significant speedup compared to the state-of-the-art shooting solver dftatom. An efficient, modular Fortran 2008 implementation, is provided under an open source, permissive license, including examples and tests, wherein particular emphasis is placed on the independence (no global variables), reusability, and generality of the individual routines.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Thermodynamics of Plutonium Monocarbide from Anharmonic and Relativistic Theory

Thermodynamics of plutonium monocarbide is studied from first-principles theory that includes relativistic electronic structure and anharmonic lattice vibrations. Density-functional theory (DFT) is expanded to include orbital-orbital coupling in addition to the relativistic spin-orbit interaction for the electronic structure and it is combined with anharmonic, temperature dependent, lattice dynamics derived from the self-consistent ab initio lattice dynamics (SCAILD) method. The obtained thermodynamics are compared to results from simpler quasi-harmonic theory and experimental data. Formation enthalpy, specific heat, and Gibbs energy calculated from the anharmonic model are validated by direct comparison with a calculation of phase diagram (CALPHAD) assessment of PuC and sub-stochiometric PuC0.896. Overall, the theory reproduces CALPHAD results and measured data for PuC rather well, but the comparison is hampered by the sub-stoichiometric nature of plutonium monocarbide. It was found that a bare theoretical approach that ignores spin-orbit and orbital-orbital coupling (orbital polarization) of the plutonium 5f electrons promotes too soft phonons and Gibbs energies that are incompatible with that of the CALPHAD assessment of the experimental data. The investigation of PuC suggests that the electronic structure is accurately described by plutonium 5f electrons as “band like” and delocalized, but correlate through spin polarization, orbital polarization, and spin-orbit coupling, in analogy to previous findings for plutonium metal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Pinning of energy transitions of defects, complexes, and surface states in AlGaN alloys

In this work, we determine the dependence of the defect transition energies, electronic bands, and surface charge neutrality levels in AlGaN. With Vacuum level as reference, we show that energy transitions of localized defects and the surface Fermi level are independent of the alloy composition as electronic bands diverge with the increase in the bandgap as a function of alloy composition. The invariance of localized states on the alloy composition creates a convenient internal reference energy with respect to which other energy states may be measured. We demonstrate a higher generality to the universality rule with the independence of deep transition states of otherwise shallow donor type defects [(+1/+3) transition for VN] and defect complexes (C N +Si III ) in addition to the earlier predicted independent nature of mid-gap states when they are either the antibonding state between cationic impurities and host anion or acceptors at anion sites.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Spatial locality of electronic correlations in LiFeAs

In this work, we address the question of the degree of spatial nonlocality of the self-energy in the iron-based superconductors, a subject which is receiving considerable attention. Using LiFeAs as a prototypical example, we extract the self-energy from angular-resolved photoemission spectroscopy data. We use two distinct electronic structure references: density functional theory in the local density approximation and linearized quasiparticle self-consistent GW (LQSGW). We find that with the LQSGW reference, spatially local dynamical correlations provide a consistent description of the experimental data, and account for some surprising aspects of the data such as the substantial out-of-plane dispersion of the electron Fermi surface having dominant xz/yz character. Hence, correlations effects can be separated into static nonlocal contributions well described by LQSGW and dynamical local contributions. Hall effect and resistivity data are shown to be consistent with this description.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Applications of Nickelate perovskites for neuromorphic computing from electronic structure and Machine Learning

While the limit of Moore's law is presently being reached with current microelectronic technologies, we need to develop new paradigms that overcome this limitation. In that respect, neuromorphic computing is a concept that emulates the neural behavior and response of the human brain, and it has been recognized as a promising alternative approach. In this research project, we will perform multi-fidelity scale bridging to explore the potential use of materials with metal to insulator transition for neuromorphic applications. In particular, rare earth nickelates are promising for such purposes, as the transition in these materials is quite sensitive to a broad set of different external stimuli. Our multi-fidelity approach will bridge the high-fidelity electronic structure calculations with classical potentials. We will bridge dynamical mean field theory with a classical atomistic representation via a deep learning force field. The neural network is trained with energies, charges, and forces obtained by accurate electronic structure theories based on Dynamical Mean Field Theory. The configurational space is generated from known crystal phases, ab initio molecular dynamics with exchange-correlation functionals corrected with the Hubbard model, disordered phases with different concentrations of oxygen vacancies, and nonsymmetrical positions and induced strain by grain interfaces or contact with a substrate. Strategies to train the model with a reduced number of training examples are obtained from active learning methods, and new structures for improving the learning process are generated by using machine learning autoencoders. This classical potential will be validated through a diversity of electronic structure methods and represents an important step to combine the flexibility and accuracy of first-principles with the speed of classical potentials. The generated multi-fidelity surrogate model will be used to understand the role of strain, oxygen vacancies, proton doping, the variation of the crystal phase, substrate effects, vibrational effects as the octahedral rotation, grain boundaries and defect effects on the response of a Metal to Insulator Transition (MIT) in correlated materials. Long time and large-scale simulations will help understand the role of different stimuli to control the hysteresis of the MIT, as it has been experimentally suggested. Selected configurations will be analyzed with higher-level theories to provide an accurate electronic description and to study how the orbitals and charges are rearranged under different conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

ELECTRONIC STRUCTURE METHODS AND PROTOCOLS WITH APPLICATION TO DYNAMICS, KINETICS AND THERMOCHEMISTRY

Hydrocarbon combustion involves the reaction dynamics of a tremendous number of species beginning with many-component fuel mixtures and proceeding via a complex system of intermediates to form primary and secondary products. Combustion conditions corresponding to new advanced engines and/or alternative fuels rely increasingly on autoignition and low-temperature-combustion chemistry. In these regimes various transient radical species such as HO2, ROO·, ·QOOH, HCO, NO2, HOCO, and Criegee intermediates play important roles in determining the detailed as well as more general dynamics. A clear understanding and accurate representation of these processes is needed for effective modeling. Given the difficulties associated with making reliable experimental measurements of these systems, computation can play an important role in developing these energy technologies. Accurate calculations have their own challenges since even within the simplest dynamical approximations such as transition state theory, the rates depend exponentially on critical barrier heights and these may be sensitive to the level of quantum chemistry. Moreover, it is well-known that in many cases it is necessary to go beyond statistical theories and consider the dynamics. Quantum tunneling, resonances, radiative transitions, and non-adiabatic effects governed by spin-orbit or derivative coupling can be determining factors in those dynamics. Building upon progress made during a period of prior support through the DOE Early Career Program, this project combines developments in the areas of potential energy surface (PES) fitting and multistate multireference quantum chemistry to allow spectroscopically and dynamically/kinetically accurate investigations of key molecular systems (such as those mentioned above), many of which are radicals with strong multireference character and have the possibility of multiple electronic states contributing to the observed dynamics. An ongoing area of investigation is to develop general strategies for robustly convergent electronic structure theory for global multichannel reactive surfaces including diabatization of energy and other relevant surfaces such as dipole transition. Combining advances in ab initio methods with automated interpolative PES fitting allows the construction of high-quality PESs (incorporating thousands of high-level data) to be done rapidly through parallel processing on high-performance computing (HPC) clusters. In addition, new methods and approaches to electronic structure theory will be developed and tested through applications. This project will explore limitations in traditional multireference calculations (e.g., MRCI) such as those imposed by internal contraction, lack of high-order correlation treatment and poor scaling. Methods such as DMRG-based extended active-space CASSCF and various Quantum Monte Carlo (QMC) methods will be applied (including VMC/DMC and FCIQMC). Insight into the relative significance of different orbital spaces and the robustness of application of these approaches on leadership class computing architectures will be gained. Synergy with other components of this research program such as automated PES fitting and multireference quantum chemistry will be used to address challenges encountered by the standard approaches to computational thermochemistry (those being single-reference quantum chemistry and perturbative treatments of the anharmonic vibrational energy, which break down for some cases of electronic structure or floppy strongly coupled vibrational modes).

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

The Li–F–H ternary system at high pressures

Evolutionary crystal structure prediction searches have been employed to explore the ternary Li-F-H system at 300 GPa. Metastable phases were uncovered within the static lattice approximation, with LiF 3 H 2 , LiF 2 H, Li 3 F 4 H, LiF 4 H 4 , Li 2 F 3 H and LiF 3 H lying within 50 meV/atom of the 0 K convex hull. All of these phases contain H n F¯ n+1 (n = 1; 2) anions, and Li + cations. Other structural motifs such as LiF slabs, $H$ $^{+}_{3}$ molecules and F δ- ions are present in some of the low enthalpy Li-F-H structures. The bonding within the H n F¯ n+1 molecules, which may be bent or linear, symmetric or asymmetric, is analyzed. The five phases closest to the hull are insulators, while LiF 3 H is metallic and predicted to have a vanishingly small superconducting critical temperature. Li 3 F 4 H is predicted to be stable at zero pressure. Furthermore, this study lays the foundation for future investigations of the role of temperature and anharmonicity on the stability and properties of compounds and alloys in the Li-F-H ternary system.

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

Observation of spatially resolved Rashba states on the surface of CH 3 NH 3 PbBr 3 single crystals

Hybrid organic-inorganic perovskites (HOIPs) are prime candidates for studying Rashba effects due to the heavy metal and halogen atoms in their crystal structure coupled with predicted inversion symmetry breaking. Nevertheless, observation of the Rashba effect in cubic CH 3 NH 3 PbBr 3 single crystals that possess bulk inversion symmetry is the subject of extensive debate due to the lack of conclusive experiments and theoretical explanations. In this work, we provide experimental evidence that Rashba state in cubic CH 3 NH 3 PbBr 3 single crystals at room temperature occurs exclusively on the crystal surface and depends on specific surface termination that results in local symmetry breaking. We demonstrate this using a suite of spatially resolved and depth-sensitive techniques, including circular photogalvanic effect, inverse spin Hall effect, and multiphoton microscopy, that are supported by first principle calculations. Our work suggests using surface Rashba states in these materials for spintronic applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Hierarchical Multi-agent Large Language Model Reasoning for Autonomous Heterogeneous Catalyst Discovery

Artificial intelligence is reshaping scientific exploration, but most methods automate procedural tasks without engaging in scientific reasoning, limiting autonomy in discovery. We demonstrate that hierarchical agentic large language model reasoning can efficiently drive simulation and scientific exploration. Across two chemical applications, CO adsorption on Cu surface transition metal adatoms and on M–N–C catalysts, reasoning-guided exploration reduces required atomistic simulations by up to 90% relative to heuristic or random selection. Comparisons across single-agent, multi-agent, and stochastic baselines show that hierarchical strategies yield more coherent and information-efficient search trajectories. Reasoning traces reveal chemically grounded decisions that cannot be explained by semantic bias or stochastic sampling. We realize these agentic reasoning strategies in Materials Agents for Simulation and Theory in Electronic-structure Reasoning (MASTER), a multimodal system that translates natural language into density functional theory workflows. Altogether, multi-agent collaboration accelerates heterogeneous catalyst discovery and marks a step toward more autonomous, reasoning-guided scientific exploration.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

Dirac's equation and its implications for density functional theory based calculations of materials containing heavy elements

Electronic structure calculations based on density functional theory (DFT) give quantitatively accurate predictions of properties of most materials containing light elements. For heavy materials, and in particular for f -electron systems, DFT based methods can fail both qualitatively and quantitatively for two distinct reasons: their failure to describe confinement effects arising from localized f -electron behavior and their incomplete or approximate treatment of relativity. In addition, different methods for incorporating relativistic effects, which give identical results in most light materials, can give different predictions in heavy elements. In order to develop a quantitative capability for calculating the properties of these materials, it is essential to separate the predictions of the underlying equations from the uncertainty introduced in approximations used in computation. Working toward that goal, here we have developed a code, called dirac-fp , which is based directly on solving the Dirac-Kohn-Sham equations and uses the full-potential linear muffin-tin orbital (FP-LMTO) approach to electronic structure. In order to assess the performance of dirac-fp , we perform calculations on three different face-centered cubic materials using different approximate treatments of relativity: the scalar relativistic (SR) approach commonly used in most solid-state DFT codes, the scalar relativistic plus spin-orbit coupling corrections (SR+SO) approach which includes spin-orbit coupling self-consistently using the SR states inside the muffin tins, and the Dirac-Kohn-Sham (Dirac) approach implemented in dirac-fp . Performing calculations on thorium, in which relativistic effects should be strong, aluminum, in which relativistic effects should be negligible, and gold, in which relativistic effects play an intermediate role, we find that the Dirac approach is able to provide theoretically consistent results in the electronic structure and ground-state properties across all three materials.

3-dimensional systems↗
DOE OSTI2024

Electronic Band Structure of a Superconducting Nickelate Probed by the Seebeck Coefficient in the Disordered Limit

Superconducting nickelates are a new family of strongly correlated electron materials with a phase diagram closely resembling that of superconducting cuprates. While analogy with the cuprates is natural, very little is known about the metallic state of the nickelates, making these comparisons difficult. We probe the electronic dispersion of thin-film superconducting five-layer ( n = 5 ) and metallic three-layer ( n = 3 ) nickelates by measuring the Seebeck coefficient S . We find a temperature-independent and negative S / T for both n = 5 and n = 3 nickelates. These results are in stark contrast to the strongly temperature-dependent S / T measured at similar electron filling in the cuprate La 1.36 Nd 0.4 Sr 0.24 CuO 4 . The electronic structure calculated from density-functional theory can reproduce the temperature dependence, sign, and amplitude of S / T in the nickelates using Boltzmann transport theory. This demonstrates that the electronic structure obtained from first-principles calculations provides a reliable description of the fermiology of superconducting nickelates and suggests that, despite indications of strong electronic correlations, there are well-defined quasiparticles in the metallic state. Finally, we explain the differences in the Seebeck coefficient between nickelates and cuprates as originating in strong dissimilarities in impurity concentrations. Our study demonstrates that the high elastic scattering limit of the Seebeck coefficient reflects only the underlying band structure of a metal, analogous to the high magnetic field limit of the Hall coefficient. This opens a new avenue for Seebeck measurements to probe the electronic band structures of relatively disordered quantum materials. Published by the American Physical Society 2024

Grissonnanche, G. (ORCID:0000000251153125)↗
DOE OSTI2025

Doping Effects on Multivalence States, Electronic Structure, and Optical Band Gap in LaCrO 3 under Varied Atmospheres: An Integrated Experimental and Density Functional Theory Study

Doping effects on the valence state, electronic structure, and optical band and the effects on electrical conductivity were studied on the doped lanthanum chromite (LaCrO 3 ) system. The specific compositions studied were La 1–x Ca x CrO 3 (LCCx), La 1–x Sr x CrO 3 (LSCx), and La 0.8 Sr 0.2 Cr 1–x Mn x O 3 (LSCMx) (0.1 ≤ x ≤ 0.4). The powders were synthesized using a modified Pechini sol–gel method, and the ceramic samples were densified using a reactive sintering method resulting in densities >96% theoretical. X-ray photoelectron spectroscopy (XPS) was completed to characterize the defect states and cationic valence compensation as a result of divalent (Ca 2+ or Sr 2+ ) and trivalent (Mn 3+ ) substitutions. XPS was completed for samples tested in oxidizing and reducing atmospheres (up to 1500 °C), which provided insights into the oxidation state transitions induced by the Ca 2+ and Sr 2+ dopants. The work notably demonstrated, for the first time, the oxidation/reduction transitions of Cr 4+ to Cr 3+ in Sr 2+ /Mn 3+ co-doped samples under reducing atmospheres. Reflectance UV–vis spectrophotometry optical band gap measurements were also completed for the same materials; a decrease in the optical band gap (2.81–3.12 eV) was shown with increased substitution, suggesting electronic structure modifications in the LaCrO 3 perovskite. Density functional theory calculations validated experimental trends, predicting a diminishing band gap with a rising dopant concentration. The transition in Cr oxidation states was attributed to the presence of divalent/trivalent cations. These findings contribute some insights into methods to tune the LaCrO 3 electrical properties for various low- and high-temperature applications.

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

Correlating Crystallography, Magnetism, and Electronic Structure Across Anhysteretic First-Order Phase Transition in Pr 2 In

Here, temperature-dependent powder X-ray diffraction and magnetization measurements of Pr 2 In conclusively prove that the unusual anhysteretic first-order paramagnetic-ferromagnetic phase transition in the compound is related to concurrent changes in both the magnetic and crystallographic lattices. At the same time, the hexagonal Ni 2 In-type structure is stable at least between 6 and 298 K, including at T C = ~57 K. From the density functional theory calculations, the electronic structure of the compound is extraordinarily sensitive to minor changes in lattice parameters that occur across the phase transition, revealing the origin of strong magnetoelastic coupling. In the vicinity of T C , the maximum entropy change, ΔS Max = -16 J Kg -1 K -1 induced by a moderate magnetic field change of 20 kOe (ΔS Max = -20 J Kg -1 K -1 for 50 kOe magnetic field change) is comparable to other known potentially functional materials demonstrating large cryogenic magnetocaloric effect.

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

High-Temperature Thermodynamics Modeling of Graphite

We present high-temperature thermodynamic properties for graphite from first-principles anharmonic theory. The ab initio electronic structure is obtained from density-functional theory coupled to a lattice dynamics method that is used to model anharmonic lattice vibrations. This combined approach produces free energies and specific heats for graphite that compare well with available experiments and results from models that empirically represent experimental data, such as CALPHAD. We show that anharmonic theory for the phonons is essential for accurate thermodynamic quantities above about 1000 K.

36 MATERIALS SCIENCE↗