Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “electronic structure method”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
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 OSTI2022

Electronic Structure Theory and Novel Materials

This grant supported research on electronic structure and materials theory, with focus on three main issues: (i) novel techniques to deal with correlation in the electronic ground-state, (ii) topological materials, (iii) the phase diagram of lattice spin models. Regarding (i), we applied to the homogeneous electron liquid an approach that we previously developed in the context of molecular systems. In this scheme the electronic occupation probabilities and the natural spin orbitals are used to construct an approximate two-body density matrix for the electronic ground-state. Regarding (ii) we used standard electronic structure methods based on density functional theory to model topological materials and interpret experimental observations. Finally, regarding (iii) we further developed a numerical approach to compute the renormalized couplings within real space renormalization group theory in the context of lattice spin models. The main findings were the following. (i) We found that with our approximate two-body density matrix, which works well for small molecules, is not sufficiently accurate for condensed phase systems. Missing a systematic way of improving on the adopted approximations, we decided not to pursue this approach. (ii) We performed two studies. In one, we investigated the influence of Te defects on the topological properties of a WTe2 monolayer, finding that while Te vacancies, even in modest concentration, destroy the topological character, Te adatoms do not, consistent with a recent experiment. In another study, we predicted Weyl semimetal character and strong anomalous Hall effect in the Heusler compensated ferrimagnet Ti2MnAl. (iii) We developed a new Monte Carlo method to do real space renormalization group calculations for lattice spin models. We subsequently extended the scheme to deal with lattice spin models in presence of quenched disorder, finding that the approach can distinguish systems with finite and strong disorder. In the finite disorder case, the method allows one to find with good approximation the critical coupling distribution and the critical exponents.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Multi-site reaction dynamics through multi-fragment density matrix embedding

The practical description of disordered chemical reactions, where the reactions involve multiple species at multiple sites, is presently a challenge using correlated electronic structure methods due to their high computational cost and steep scaling. Here, we describe the gradient theory of multi-fragment density matrix embedding theory, which potentially provides a minimal computational framework to model such processes at the correlated electron level. We present the derivation and implementation of the gradient theory, its validation on model systems and chemical reactions using density matrix embedding, and its application to a molecular dynamics simulation of proton transport in a small water cluster, a simple example of multi-site reaction dynamics.

Chemistry↗
DOE OSTI2024

Diffusion quantum Monte Carlo approach to the polaritonic ground state

Making and using polaritonic states (i.e., hybrid electron-photon states) for chemical applications has recently become one of the most prominent and active fields that connects the communities of chemistry and quantum optics. Modeling of such polaritonic phenomena using ab initio approaches calls for new methodologies, leading to the reinvention of many commonly used electronic structure methods, such as Hartree-Fock, density functional, and coupled cluster theories. Here in this work, we explore the formally exact diffusion quantum Monte Carlo approach to obtain numerical solutions to the polaritonic ground state during the dissociation of the H 2 molecular system. Here, we examine various electron-nuclear-photon properties throughout the dissociation, such as changes to the minimum of the cavity Born-Oppenheimer surface, the localization of the electronic wave function, and the average mode occupation. Finally, we directly compare our results to that obtained with state-of-the-art, yet approximate, polaritonic coupled cluster approaches.

74 ATOMIC AND MOLECULAR PHYSICS↗
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 OSTI2022

Improving the Accuracy of Atomistic Simulations of the Electrochemical Interface

Atomistic simulation of the electrochemical double layer is an ambitious undertaking, requiring quantum mechanical description of electrons, phase space sampling of liquid electrolytes, and equilibration of electrolytes over nanosecond time scales. All models of electrochemistry make different trade-offs in the approximation of electrons and atomic configurations, from the extremes of classical molecular dynamics of a complete interface with point-charge atoms to correlated electronic structure methods of a single electrode configuration with no dynamics or electrolyte. Here, we review the spectrum of simulation techniques suitable for electrochemistry, focusing on the key approximations and accuracy considerations for each technique. Additionally, we discuss promising approaches, such as enhanced sampling techniques for atomic configurations and computationally efficient beyond density functional theory (DFT) electronic methods, that will push electrochemical simulations beyond the present frontier.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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

Predicting the X-ray Absorption Spectrum of Ozone with Single Configuration State Functions

X-ray absorption spectra (XAS) of biradicaloid species are often thought to represent a challenge to theoretical methods. This has led to the testing of recently developed multireference techniques on the XAS of ozone, but reproduction of the experimental spectral profile has proven difficult. We utilize a minimal model consisting of a single configuration state function (CSF) per excited state to model core-level excitations of ozone, with the orbitals of each CSF optimized using the restricted open-shell Kohn–Sham (ROKS) method. This protocol leads to semiquantitative agreement with experimental XAS. In fact, we find that low-lying core-hole excited states in biradicaloids can be approximated with individual CSFs, despite the presence of multireference character in the ground state. Here, we also report that the 1s → π* and 1s → σ* transitions have quite distinct widths for O 3 . This reveals the importance of sampling over a representative range of geometries from the vibrational ground state for properly assessing the accuracy of electronic structure methods against experiments instead of the popular procedure of uniformly broadening stick spectra at the equilibrium geometry.

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

The symmetric quasi-classical model using on-the-fly time-dependent density functional theory within the Tamm–Dancoff approximation

The primary computational challenge when simulating nonadiabatic ab initio molecular dynamics is the unfavourable compute costs of electronic structure calculations with molecular size. Simple electronic structure theories, like time-dependent density functional theory within the Tamm–Dancoff approximation (TDDFT/TDA), alleviate this cost for moderately sized molecular systems simulated on realistic time scales. Although TDDFT/TDA does have some limitations in accuracy, an appealing feature is that, in addition to including electron correlation through the use of a density functional, the cost of calculating analytic nuclear gradients and nonadiabatic coupling vectors is often computationally feasible even for moderately sized basis sets. Here in this work, some of the benefits and limitations of TDDFT/TDA are discussed and analysed with regard to its applicability as a ‘back-end’ electronic structure method for the symmetric quasi-classical Meyer–Miller model (SQC/MM). In order to investigate the benefits and limitations of TDDFT/TDA, SQC/MM is employed to predict and analyse a prototypical example of excited-state hydrogen transfer in gas-phase malonaldehyde. Then, the ring-opening dynamics of selenophene are simulated, which highlight some of the deficiencies of TDDFT/TDA. Additionally, some new algorithms are proposed that speed up the calculation of analytic nuclear gradients and nonadiabatic coupling vectors for a set of excited electronic states.

molecular dynamics↗
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↗
DOE OSTI2022

Nonadiabatic Molecular Dynamics by Multiconfiguration Pair-Density Functional Theory

Herein, we present the first implementation of multiconfiguration pair-density functional theory (MC-PDFT) ab initio molecular dynamics. MC-PDFT is a multireference electronic structure method that in many cases has a similar accuracy (or even better accuracy) the complete active space second-order perturbation theory (CASPT2) at a significantly lower computational cost. In this study, we introduced MC-PDFT analytical gradients into the SHARC molecular dynamics program for ab initio, nonadiabatic molecular dynamics simulations. We verify our implementation by examining the intersystem crossing dynamics of thioformaldehyde, and we observe excellent agreement with recent CASPT2 and experimental findings. Moreover, with MC-PDFT, we could perform dynamics simulations with the 12 electron in 10 orbitals active space that was computationally too expensive for direct dynamics with CASPT2.

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

Cartesian equivariant representations for learning and understanding molecular orbitals

Qualitative and quantitative orbital properties such as bonding/antibonding character, localization, and orbital energies are critical to how chemists understand reactivity, catalysis, and excited-state behavior. Despite this, representations of orbitals in deep learning models have been very underdeveloped relative to representations of molecular geometries and Hamiltonians. Here, we apply state-of-the-art equivariant deep learning architectures to the task of assigning global labels to orbitals, namely energies characterizations, given the molecular coefficients from Hartree–Fock or density functional theory. The architecture we have developed, the Cartesian Equivariant Orbital Network (CEONET), shows how molecular orbital coefficients are readily featurized as equivariant node features common to all graph-based machine-learned potentials. We find that CEONET performs well at predicting difficult quantitative labels such as the orbital energy and orbital entropy. Furthermore, we find that the CEONET representation provides an intuitive latent space for differentiating orbital character for the qualitative assignment of e.g. bonding or antibonding character. In addition to providing a useful representation for further integrating deep learning with electronic structure theory, we expect CEONET to be useful for automatizing and interpreting the results of advanced electronic structure methods such as complete active space self-consistent field theory. In particular, the ability of CEONET to infer multireference character via the orbital entropy paves the way toward the machine-learned selection of active spaces.

chemical reactions↗
DOE OSTI2025

Advances and Challenges of SCAN and r 2 SCAN Density Functionals in Transition‐Metal Compounds

ABSTRACT Transition‐metal compounds (TMCs) with open‐shell d ‐electrons are characterized by a complex interplay of lattice, charge, orbital, and spin degrees of freedom, giving rise to various fascinating applications. Often exhibiting exotic properties, these compounds are commonly classified as correlated systems due to strong inter‐electronic interactions called Hubbard U . This inherent complexity presents significant challenges to Kohn‐Sham density functional theory (KS‐DFT), the most widely used electronic structure method in condensed matter physics and materials science. While KS‐DFT is, in principle, exact for the ground‐state total energy, its exchange‐correlation energy must be approximated in practice. The mean‐field nature of KS implementations, combined with the limitations of current exchange‐correlation density functional approximations, has led to the perception that DFT is inadequate for correlated systems, particularly TMCs. Consequently, a common workaround involves augmenting DFT with an on‐site Hubbard‐like U correction. In recent years, the strongly constrained and appropriately normed (SCAN) density functional, along with its refined variant r 2 SCAN, has achieved remarkable progress in accurately describing the structural, energetic, electronic, magnetic, and vibrational properties of TMCs, challenging the traditional perception of DFT's limitations. This review explores the design principles of SCAN and r 2 SCAN, highlights their key advancements in studying TMCs, explains the mechanisms driving these improvements, and addresses the remaining challenges in this evolving field.

Zhang, Yubo [Minjiang Collaborative Center for The↗
DOE OSTI2024

Machine learning for accuracy in density functional approximations

Machine learning techniques have found their way into computational chemistry as indispensable tools to accelerate atomistic simulations and materials design. In addition, machine learning approaches hold the potential to boost the predictive power of computationally efficient electronic structure methods, such as density functional theory, to chemical accuracy and to correct for fundamental errors in density functional approaches. In this paper, recent progress in applying machine learning to improve the accuracy of density functional and related approximations is reviewed. Promises and challenges in devising machine learning models transferable between different chemistries and materials classes are discussed with the help of examples applying promising models to systems far outside their training sets.

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

Insulating band gaps both below and above the Néel temperature in d -electron LaTi O 3 , LaV O 3 , SrMn O 3 , and LaMn O 3 perovskites as a symmetry-breaking event

Metal d-electron oxides having an odd number of electrons per cell should exhibit band degeneracy at the Fermi energy, making them, in band theory, formally metallic. In many cases, however, these are “false metals”, as evidenced by the observation that many ABO 3 oxide perovskites with a magnetic 3d B-atom are observed to be insulators both below and above the Néel temperature. These inconsistencies between experimental observation and expectation have historically been resolved by invoking degeneracy-breaking physics, based largely on -pure electron effects, such as strong interelectronic correlation for d-electron compounds (the Mott mechanism). Such explanations generally consider microscopic lattice or magnetic degrees of freedom (m-DOF) as largely passive spectators, not a cause of the formally metallic being an insulator. Yet, it has long been known that ABO 3 perovskites can manifest an arrangement of m-DOFs in the form of octahedral tilting, bond dimerization, Jahn-Teller distortions, and ordering of local magnetic moments. It appears reasonable that such structural and magnetic local degrees of freedom need to be allowed to compete with purely electronic strong correlation. To answer this question, we explored a range of d-electron oxide perovskites exemplified by the archetypes LaTiO 3 , LaVO 3 , SrMnO 3 , and LaMnO 3 with 1, 2, 3, or 4 d -electrons, respectively. Using a mean-field-like electronic structure method (here, density functional theory), we find that a combination of magnetic symmetry breaking with structural distortions can account for insulating band gaps in this series while at the same time correctly predicting for the control case an intrinsic paramagnetic metal in SrVO 3 as symmetry-breaking is insufficiently strong to remove the degeneracy. Finally, this indicates that calculating quantitatively local magnetic and positional symmetry-breaking motifs in unit cells that avoid averaging at the outset over the low symmetry motifs can provide consistent trends in a Mott transition without Mott U.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Exploring the scaling limitations of the variational quantum eigensolver with the bond dissociation of hydride diatomic molecules

Abstract Materials simulations involving strongly correlated electrons pose fundamental challenges to state‐of‐the‐art electronic structure methods but are hypothesized to be the ideal use case for quantum computing algorithms. To date, no quantum computer has simulated a molecule of a size and complexity relevant to real‐world applications, despite the fact that the variational quantum eigensolver (VQE) algorithm can predict chemically accurate total energies. Nevertheless, because of the many applications of moderately sized, strongly correlated systems, such as molecular catalysts, the successful use of the VQE stands as an important waypoint in the advancement toward useful chemical modeling on near‐term quantum processors. In this paper, we take a significant step in this direction. We lay out the steps, write, and run parallel code for an (emulated) quantum computer to compute the bond dissociation curves of the TiH, LiH, NaH, and KH diatomic hydride molecules using the VQE. TiH was chosen as a relatively simple chemical system that incorporates d orbitals and strong electron correlation. Because current VQE implementations on existing quantum hardware are limited by qubit error rates, the number of qubits available, and the allowable gate depth, recent studies using it have focused on chemical systems involving s and p block elements. Through VQE + UCCSD calculations of TiH, we evaluate the near‐term feasibility of modeling a molecule with d‐orbitals on real quantum hardware. We demonstrate that the inclusion of d‐orbitals and the use of the UCCSD ansatz, which are both necessary to capture the correct TiH physics, dramatically increase the cost of this problem. We estimate the approximate error rates necessary to model TiH on current quantum computing hardware using VQE + UCCSD and show them to likely be prohibitive until significant improvements in hardware and error correction algorithms are available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗