Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Electronic structure 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
NASA NTRS1994

Theoretical determination of chemical rate constants using novel time-dependent methods

The work completed within the grant period 10/1/91 through 12/31/93 falls primarily in the area of reaction dynamics using both quantum and classical mechanical methodologies. Essentially four projects have been completed and have been or are in preparation of being published. The majority of time was spent in the determination of reaction rate coefficients in the area of hydrocarbon fuel combustion reactions which are relevant to NASA's High Speed Research Program (HSRP). These reaction coefficients are important in the design of novel jet engines with low NOx emissions, which through a series of catalytic reactions contribute to the deterioration of the earth's ozone layer. A second area of research studied concerned the control of chemical reactivity using ultrashort (femtosecond) laser pulses. Recent advances in pulsed-laser technologies have opened up a vast new field to be investigated both experimentally and theoretically. The photodissociation of molecules adsorbed on surfaces using novel time-independent quantum mechanical methods was a third project. And finally, using state-of-the-art, high level ab initio electronic structure methods in conjunction with accurate quantum dynamical methods, the rovibrational energy levels of a triatomic molecule with two nonhydrogen atoms (HCN) were calculated to unprecedented levels of agreement between theory and experiment.

Dateo, Christopher E.↗
DOE OSTI2024

The Role of Spin–Orbit Coupling in the Linear Absorption Spectrum and Intersystem Crossing Rate Coefficients of Ruthenium Polypyridyl Dyes

The successful use of molecular dyes for solar energy conversion requires efficient charge injection, which in turn requires the formation of states with sufficiently long lifetimes (e.g., triplets). The molecular structure elements that confer this property can be found empirically, however computational predictions using ab initio electronic structure methods are invaluable to identify structure-property relations for dye sensitizers. The primary challenge for simulations to elucidate the electronic and nuclear origins of these properties is a spin-orbit interaction which drives transitions between electronic states. Here, in this work, we present a computational analysis of the spin-orbit corrected linear absorption cross sections and intersystem crossing rate coefficients for a derivative set of phosphonated tris(2,2'-bipyridine)ruthenium(2+) dye molecules. After sampling the ground state vibrational distributions, the predicted linear absorption cross sections indicate that the mixture between singlet and triplet states plays a crucial role in defining the line shape of the metal-to-ligand charge transfer bands in these derivatives. Additionally, an analysis of the intersystem crossing rate coefficients suggests that transitions from the singlet into the triplet manifolds are ultrafast with rate coefficients on the order of 10 13 s -1 for each dye molecule.

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

An NV − center in magnesium oxide as a spin qubit for hybrid quantum technologies

Recent predictions suggest that oxides, such as MgO and CaO, could serve as hosts of spin defects with long coherence times and thus be promising materials for quantum applications. However, in most cases, specific defects have not yet been identified. Here, by using a high-throughput first-principles framework and advanced electronic structure methods, we identify a negatively charged complex between a nitrogen interstitial and a magnesium vacancy in MgO with favorable electronic and optical properties for hybrid quantum technologies. We show that this NV − center has stable triplet ground and excited states, with singlet shelving states enabling optical initialization and spin-dependent readout. We predict several properties, including absorption, emission, and zero-phonon line energies, as well as zero-field splitting tensor, and hyperfine interaction parameters, which can aid in the experimental identification of this defect. Our calculations show that due to a strong pseudo-Jahn Teller effect and low-frequency phonon modes, the NV − center in MgO is subject to a substantial vibronic coupling. We discuss design strategies to reduce such coupling and increase the Debye-Waller factor, including the effect of strain and the localization of the defect states. We propose that the favorable properties of the NV − defect, along with the technological maturity of MgO, could enable hybrid classical-quantum applications, such as spintronic quantum sensors and single qubit gates.

36 MATERIALS SCIENCE↗
NASA NTRS2016

Computing Highly Accurate Spectroscopic Line Lists for Characterization of Planetary Atmospheres: CO2 and SO2 Line Lists Needed for Modeling Venus

Over the last decade, it has become apparent that the most effective approach for determining highly accurate rotational and rovibrational line lists for molecules of interest in planetary atmospheres and other astrophysical environments is through a combination of highresolution laboratory experiments coupled with state-of-the art ab initio quantum chemistry methods. The approach involves computing the most accurate potential energy surface (PES) possible using state-of-the art electronic structure methods, followed by computing rotational and rovibrational energy levels using an exact variational method to solve the nuclear Schrödinger equation. Then, reliable experimental data from high-resolution experiments is used to refine the ab initio PES in order to improve the accuracy of the computed energy levels and transition energies. From the refinement step, we have been able to achieve an accuracy of approximately 0.015 cm-1 for rovibrational transition energies, and even better for purely rotational transitions. This combined "experiment / theory" approach allows for determination of essentially a complete line list, with hundreds of millions of transitions, and having the transition energies and intensities be highly accurate. Our group has successfully applied this approach to determine highly accurate line lists for NH3, CO2 and isotopologues, and SO2 and isotopologues. Here I will report our latest results for CO2 and SO2 including all isotopologues. Comparisons to the available data in HITRAN2012 and other available databases will be shown, though we note that our line lists for SO2 are significantly more complete than any other databases. Since it is important to span a large temperature range in order to model the spectral signature of Venus as well as exoplanets, we will demonstrate how the spectra change on going from low temperatures (100 K) to higher temperatures (500 K to 1500 K).

potential energy surface (PES)5405 PLANETARY SCIEN↗
DOE OSTI2025

Actinide 5f Occupations: The Case of PuO 2

In actinide chemistry, the formal number of open-shell 5f electrons, n open , is a well-defined quantity with an integer value. The effective 5f occupation, n f , additionally takes donation and back-donation into account, generally has a non-integer value, and has varying numerical definitions. The present study explores the important distinction between n f and n open in actinide chemistry with the example of PuO2, by using electronic structure methods with a relativistic Hamiltonian in combination with experimental Pu M 5 -edge high-energy-resolution X-ray absorption and emission spectroscopic data. The total donation to the metal in PuO 2 is between 3.1 and 2.4 electrons, depending on the type of calculation, most of which is to the Pu 6d and 5f shells. The donation into 5f is sensitive to the approximations in the electronic structure model but likely amounts to 1.6/0.8 electrons when the diffuse regions of the 5f shell are included/excluded. Valence band resonant inelastic X-ray scattering experiments demonstrate that Pu 5f electron density is present in the valence band; thus, there is a clear experimental signature of covalent bonding in PuO 2 . Pu M 5 -edge and M 3 -edge high-energy-resolution X-ray absorption near-edge structures for Pu 3+ and Pu 4+ in an aqueous solution are compared to PuO 2 , showing that Pu in PuO 2 has an nf closer to Pu 4+ (aq).

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

Gas-Phase Reaction Pathways and Rate Coefficients for the Dichlorosilane-Hydrogen and Trichlorosilane-Hydrogen Systems

As part of NASA Ames Research Center's Integrated Process Team on Device/Process Modeling and Nanotechnology our goal is to create/contribute to a gas-phase chemical database for use in modeling microelectronics devices. In particular, we use ab initio methods to determine chemical reaction pathways and to evaluate reaction rate coefficients. Our initial studies concern reactions involved in the dichlorosilane-hydrogen (SiCl2H2--H2) and trichlorosilane-hydrogen (SiCl2H-H2) systems. Reactant, saddle point (transition state), and product geometries and their vibrational harmonic frequencies are determined using the complete-active-space self-consistent-field (CASSCF) electronic structure method with the correlation consistent polarized valence double-zeta basis set (cc-pVDZ). Reaction pathways are constructed by following the imaginary frequency mode of the saddle point to both the reactant and product. Accurate energetics are determined using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations (CCSD(T)) extrapolated to the complete basis set limit. Using the data from the electronic structure calculations, reaction rate coefficients are obtained using conventional and variational transition state and RRKM theories.

Dateo, Christopher E.↗
DOE OSTI2025

Effective bands and band-like electron transport in amorphous solids

The localization of electrons caused by atomic disorder is a well-known phenomenon. However, under which circumstances electrons remain delocalized and retain band-like characteristics even when the crystal structure is completely absent, as found in certain amorphous solids, is less well understood. Here, in this study, to probe this phenomenon, we develop a fully first-principles description of the electronic structure and charge transport in amorphous materials, which combines a representation of the amorphous state as a composite (ensemble) of local environments and the state-of-the-art many-body electronic structure methods. Using amorphous In 2 O 3 as an example, we demonstrate the accuracy of our approach in reproducing the band-like nature of the conduction electrons as well as their disorder-limited mobility. Our approach reveals the physical origins responsible for the electron delocalization and survival of the band dispersions despite the absence of long-range order.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A Bottom-Up Approach to Rational Design of Crystalline Materials: Investigation of Vibronic Coherences Underlying Exciton Dynamics in Semiconductors

In this project we uncovered structure-function relationships of donor-acceptor co-crystals used to develop next-generation optoelectronic devices. Unraveling the photodynamics of molecular crystalline materials poses many challenges for spectroscopy due to broad, overlapping features representing numerous underlying dynamical processes. This leads researchers to make many assumptions about the dynamics of a system in choosing an appropriate kinetic fitting model. Computationally, electronic structure methods are either prohibitively expensive or underdeveloped for computing the excited state structure of molecular materials, especially states that exhibit charge transfer. Researchers must therefore perform calculations of excited electronic states using truncated models of molecular materials. Here we present a joint experimental-theoretical approach to bridging the gap between the photodynamics of a molecular material and its constituent molecules. We focus our efforts on quantifying the timescales and mechanisms of photoexcitation in donor-acceptor co-crystals and donor-acceptor dimers where the lowest-lying excited state is characterized by charge transfer from the donor to the acceptor. We employ ultrafast UV pump, UV-Vis probe transient absorption spectroscopy to unravel the time-resolved spectroscopic signatures of the photodynamics in both the crystalline material and donor-acceptor dimers in solution. We perform electronic structure and excited state dynamics calculations of the dimers to inform kinetic fitting models and assign the spectral features. The photodynamics of the crystal vs. dimer systems have many similarities, enabling unprecedented insights into the formation and evolution of charge transfer excitons in the crystalline systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Symmetry breaking forms split-off flat bands in quantum oxides controlling metal versus insulator phases

The crystal structure used as input to electronic structure calculation of conventionally bonded solids consists primarily of the standard crystallographic degrees of freedom. Quantum solids sometimes have additional microscopic degrees of freedom (m-DOF) nested within the crystallographic structure. These might include local motifs including positional (Peierls dimers, deformed octahedra, Jahn-Teller distortions), magnetic (local moment configurations) and dipolar (local ferroelectric configurations). Such motifs can be observed experimentally via local probes that avoid averaging, and theoretically as distinct total energy lowering features relative to more simplified “average crystallographic structures”. Here we examine the ability of electronic structure methods independent of strong correlation physics to explain the broad phenomenology of metal-insulator selectivity in quantum oxides by energy-lowering symmetry breaking. We do this by avoiding the restriction of considering only (i) electron-electron interactions in a fixed unresponsive lattice—as done in Mott strong correlation explanations—allowing, however, (ii) local positional and magnetic motifs that coexist within a crystallographic landscape. We find in a broad range of quantum oxides with different symmetry breaking modes the formation of split-off flat bands that can also control metallic versus insulating characteristics. The split-off band effect is common to magnetic and non-magnetic materials, including both binary and ternary oxides. One finds that when such local symmetry breaking motifs are considered in mean-field-like (e.g. density functional theory) electronic structure calculations of quantum oxides, they explain many of the observed trends in metal vs insulator phases discussed otherwise in prior literature via strong correlation in symmetry unbroken view.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

How Accurate Are Simulations and Experiments for the Lattice Energies of Molecular Crystals?

Molecular crystals play a central role in a wide range of scientific fields, including pharmaceuticals and organic semiconductor devices. However, they are challenging systems to model accurately with computational approaches because of a delicate interplay of intermolecular interactions such as hydrogen bonding and Van der Waals dispersion forces. Here, by exploiting recent algorithmic developments, we report the first set of diffusion Monte Carlo lattice energies for all 23 molecular crystals in the popular and widely used X23 dataset. Comparisons with previous state-of-the-art lattice energy predictions (on a subset of the dataset) and a careful analysis of experimental sublimation enthalpies reveals that high-accuracy computational methods are now at least as reliable as (computationally derived) experiments for the lattice energies of molecular crystals. Overall, this work demonstrates the feasibility of high-level explicitly correlated electronic structure methods for broad benchmarking studies in complex condensed phase systems, and signposts a route towards closer agreement between experiment and simulation. Published by the American Physical Society 2024

Physics↗
DOE OSTI2026

Quantum effects on the dynamics and properties of soft materials

The quantum effects of nuclear and electronic motion play an important role in the structure, dynamics, and function of soft materials, yet they are difficult to capture with conventional classical simulations or static electronic–structure methods. In this work several complementary approaches for treating quantum effects in polymeric and soft–matter systems are demonstrated, with a focus being on the hydrogen-bonded networks, ion and charge transport, and photoactive chromophores. The proton transfer, tunneling, and isotope effects are captured within the reduced-dimensionality models by implementing grid-based nuclear quantum dynamics in terms of the discrete variable and Fourier bases. The nuclear quantum dynamics is extended to larger systems by employing the quantum trajectories and quantum–thermal bath schemes combined with on-the-fly electronic structure, enabling the description of high-dimensional polymeric environments at feasible cost. The dynamics in the electronic degrees of freedom, simulating the optical response in large chromophores such as chlorophylls, is performed using the real-time time-dependent density functional theory implemented in the real-space multigrid (RMG) code. These approaches are demonstrated on case studies of the proton and hydroxide transport in hydrated polymer membranes, charge transfer in conjugated polymers, and the optical spectra of chlorophyll chromophores relevant to polymerized chlorophyll materials and chlorophyll–polymer hybrids. The reviewed methods and applications highlight practical routes of including quantum effects in simulations of soft functional materials.

Garashchuk, Sophya [Univ. of South Carolina, Colum↗
DOE OSTI2026

Low-Lying Excited States of Linear All- Trans Polyenes: Insights from Analytic Gradient and Nonadiabatic Coupling Calculations Based on Multireference Configuration Interaction

Polyenes serve as a rigorous test for theoretical models and electronic structure methods, playing a key role in advancing computational and theoretical chemistry. Here, we present a high-level theoretical investigation of linear, all-trans polyenes using energy gradients and nonadiabatic coupling vectors based on an MR-CISD wave function to describe electronic transitions involving the ground state (1 1 A g – ) and three low-lying excited states (2 1 A g – , 1 1 B u + , and 2 1 B u – ) of hexatriene, octatetraene, and decapentaene. This approach enables accurate evaluation of both adiabatic and vertical excitation and emission energies, yielding results in excellent agreement with experiment, as well as locating minima on the crossing seam between adiabatic states. Our results show that vertical excitation energies to the 1 1 B u + state are blue-shifted by 0.2–0.3 eV relative to the experimental absorption maximum, whereas the vertical emission energy from the 2 1 A g – state is red-shifted by ∼0.2 eV relative to the experimental emission maximum. Upon relaxation from the Franck–Condon geometry, the 2 1 A g – state stabilizes by around 1 eV, compared to 0.2–0.3 eV for the 1 1 B u + state. An analysis of the S 1 /S 0 crossing seam in hexatriene shows that its minimum involves asymmetric backbone deformations and provides an efficient channel for ultrafast internal conversion to the ground state, consistent with the absence of detectable fluorescence in this molecule. These results demonstrate the power of analytic gradients and nonadiabatic coupling vectors based on an MR-CISD wave function for accurately characterizing the electronic structure and photophysics of polyenes.

Excited states↗
DOE OSTI2025

Weighted active space protocol for multireference machine-learned potentials

Multireference methods such as multiconfiguration pair-density functional theory accurately capture electronic correlation in systems with strong multiconfigurational character, but their cost precludes direct use in molecular dynamics. Combining these methods with machine-learned interatomic potentials (MLPs) can extend their reach. However, the sensitivity of multireference calculations to the choice of the active space complicates the consistent evaluation of energies and gradients across structurally diverse nuclear configurations. To overcome this limitation, we introduce the weighted active space protocol (WASP), a systematic approach to assign a consistent active space for a given system across uncorrelated configurations. By integrating WASP with MLPs and enhanced sampling techniques, we propose a data-efficient active learning cycle that enables the training of an MLP on multireference data. We demonstrated the approach on the TiC + -catalyzed C–H activation of methane, a reaction that poses challenges for Kohn–Sham density functional theory due to its significant multireference character. This framework enables accurate and efficient modeling of catalytic dynamics, establishing a paradigm for simulating complex reactive processes beyond the limits of conventional electronic-structure methods.

enhanced sampling↗
DOE OSTI2025

Accelerating Embedding Potential Optimization by Reconstructing the Pseudo-Valence Electron Density

Density functional embedding theory (DFET) enables use of electronic structure methods with higher accuracy than density functional theory in a local region, with applications thus far ranging from (photo/electro)catalysis to reactions in solution. DFET partitions a large collection of atoms into smaller groups that interact via a shared embedding (interaction) potential V emb , determined via functional optimization. The optimized effective potential (OEP) process used to optimize V emb is time-consuming and becomes a computational bottleneck due to sharp, oscillating features of V emb near nuclei. Here, similar to pseudopotential theory, by reconstructing electron densities used in the OEP process from smoother pseudo-valence-only (PVO) electron densities as proxies for total densities of the full system and subsystems, we can retain accuracy in the embedded electronic structure calculations while potentially reducing the overhead of V emb construction, within the projector augmented-wave (PAW) formalism. We explore three different chemical reactions as exemplars to test PVO–DFET, namely, H 2 dissociative adsorption on a Cu(111) surface, H 2 O adsorption on a Pt(111) surface, and aqueous [Ca 2+ –SO 4 2– ] ion-pair formation. The PVO approximation works well for all three systems with minimal loss of accuracy (∼10–70 meV error relative to the original exact-derivative (ED) approach) while accelerating V emb generation for the Cu and Pt systems respectively by 20× and 5×. Given proper numerical convergence parameters, the spatial distributions of differences between PVO- and ED-based V emb outside the core regions are small, explaining the exceptional agreement between the two approaches. Finally, we anticipate that this more efficient PVO–DFET approximation will be useful whenever computation of V emb is much more expensive than subsequent embedded high-level electron correlation calculations.

approximation↗
DOE OSTI2025

Higher-order Van Hove singularities in kagome topological bands

Motivated by the growing interest in band structures featuring higher-order Van Hove singularities (HOVHS), we investigate a spinless fermion kagome system characterized by nearest-neighbor (NN) and next-nearest-neighbor (NNN) hopping amplitudes. While NN hopping preserves time-reversal symmetry, NNN hopping, akin to chiral hopping on the Haldane lattice, breaks time-reversal symmetry and leads to the formation of topological bands with Chern numbers ranging from 𝐶 = ±1 to ±4. We perform analytical and numerical analysis of the energy bands near the high-symmetry points Γ, ±𝐊, and 𝐌 𝑖 (𝑖 = 1, 2, and 3), which uncover a rich and complex landscape of HOVHS, controlled by the magnitude and phase of the NNN hopping. We observe power-law divergences in the density of states (DOS), 𝜌⁡(𝜀)∼|𝜀| −𝜈 , with exponents 𝜈 = 1/2, 1/3, 1/4, which can significantly affect the anomalous Hall response at low temperatures when the Fermi level crosses the HOVHS. Additionally, the NNN hopping induces the formation of higher Chern number bands 𝐶 = ±2, ±4 in the middle of the spectrum obeying a sublattice interference whereupon electronic states are maximally localized in each of the sublattices when the momentum approaches the three high-symmetry points 𝐌 𝑖 (𝑖 = 1, 2, and 3) on the Brillouin zone boundary. Finally, this classification of HOVHS in kagome systems provides a platform to explore unconventional electronic orders induced by electronic correlations.

Chern insulators↗
DOE OSTI2024

Light-induced electron dynamics of molecules in cavities: comparison of model Hamiltonians

The rapidly expanding field of polaritonic chemistry requires accurate theoretical simulations to understand new phenomena at the atomic scale. Computing the optoelectronic properties of molecules using established electronic structure methods is a careful balance of accuracy and computational expense, and expanding these methods to quantum electrodynamics to describe coupled cavity-molecule systems is an active topic of development. Key to these methods are the Hamiltonian operators representing the photon cavity modes. The recently introduced quantum electrodynamics time-dependent configuration interaction (QED-TDCI) method allows for the combination of electron dynamics simulations with quantum electrodynamics, enabling the simulation of time-dependent optoelectronic properties of cavity-molecule systems. Using this method, a comparison of two many-state QED Hamiltonians—the Pauli-Fierz and quantum Rabi model Hamiltonians—is presented, with a particular focus on time-dependent properties in applied electric fields.

Peyton, Benjamin G. (ORCID:0000000163190751)↗
DOE OSTI2025

The electronic structure, crystal fields, and magnetic anisotropy in RECo 5 magnets

The current progress in describing rare-earth-based magnets' electronic structure and magnetic properties is discussed. We use several currently popular electronic structure methods to show the typical values of critical parameters that define the physics of RECo 5 (RE = rare earth atom) materials. The magnetic moments and magnetic anisotropy of 4f atoms are obtained using several approaches, including anisotropic 4f-charge density-constrained DFT and DFT+HI methods. We also suggest the introduction of "penalty" functional for obtaining correct variational total energy in the traditional Hund's rule-constrained DFT-based techniques. The applicability and future extensions are discussed. The proposed combination of methods is potentially suitable for high-throughput computational searches of new rare-earth-containing magnetic materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Computing Reaction Kinetics with MC-PDFT–OPESf: Combining Multireference Electronic Structure Theory and Enhanced Sampling

Accurate rate constants are crucial for understanding and optimizing catalytic reactions mediated by enzymes, metalloproteins, and heterogeneous catalysts. These systems frequently present a dual computational challenge. Multiconfigurational reaction sites require multireference techniques for the accurate treatment of the electronic structure, and high activation barriers prevent efficient sampling of unbiased reactive transitions. In this work, we combine multiconfiguration pair-density functional theory (MC-PDFT) as an accurate and efficient multireference electronic structure method with on-the-fly probability-enhanced sampling flooding (OPESf) as an enhanced sampling method capable of accelerating reactive transitions. We demonstrate the approach on the Diels–Alder [4+2] cycloaddition between cis-butadiene and ethene as a reaction characterized by a large activation barrier and multireference character. MC-PDFT–OPESf provides reaction rates in agreement with experiments at a fraction of the computational cost required by conventional unbiased ab initio calculations. Here, we propose MC-PDFT–OPESf as an efficient approach for computing kinetics in strongly correlated molecular systems.

Chemical calculations↗