Search NASA⌕ Search

SEARCH · Search NASA

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Association of Cl with C 2 H 2 by unified variable-reaction-coordinate and reaction-path variational transition-state theory

Significance A key issue in chemical kinetics is advancing the theoretical framework to handle reactions beyond the domain of textbook transition-state theory by including––for example––anharmonicity, barrierless transition states, transition states in series, and the effect of conformational flexibility on equilibrium constants. Other key issues are validating affordable electronic structure methods for direct dynamics and the direct calculation of high-pressure limiting rate constants, which are often obtainable experimentally only by extrapolation. The present article addresses all these issues for a prototype radical–molecule association reaction and thereby demonstrates how to combine improved theoretical methods to provide rate constants in cases where experimental data are uncertain or missing.

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

Refined energetic ordering for sulphate–water ( n = 3–6) clusters using high-level electronic structure calculations

This work reports refinements of the energetic ordering of the known low-energy structures of sulfate-water clusters SO 4 2- (H2O)n (n = 3-6) using high-level electronic structure methods. Coupled cluster singles and doubles with perturbative triples (CCSD(T)) is used in combination with an estimate of basis set effects up to the complete basis set limit using second order Møller-Plesset theory. Harmonic zero point energy (ZPE), included at the B3LYP/6-311++G** level, was found to have a significant effect on the energetic ordering. Limitations of the ZPE calculations, both due to electronic structure errors, and use of the harmonic approximation, probably constitute the largest remaining errors. Due to the often small energy differences between cluster isomers, and the significant role of ZPE, deuteration can alter the relative energies of low-lying structures, and, when it is applied in conjunction with calculated harmonic ZPE’s, even alters the global minimum for n = 4.

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

The role of Rydberg and continuum levels in computing high harmonic generation spectra of the hydrogen atom using time-dependent configuration interaction

We study the role of Rydberg bound-states and continuum levels in the field-induced electronic dynamics associated with the High-Harmonic Generation (HHG) spectroscopy of the hydrogen atom. Time-dependent configuration-interaction (TD-CI) is used with very large atomic orbital (AO) expansions (up to L= 4 with sextuple augmentation and off-center functions) to describe the bound Rydberg levels, and some continuum levels. To address the lack of ionization losses in TD-CI with finite AO basis sets, we employed a heuristic lifetime for energy levels above the ionization potential. The heuristic lifetime model is compared against the conventional atomic orbital treatment (infinite lifetimes), and a third approximation which is TD-CI using only the bound levels (continuum lifetimes go to zero). The results suggest that spectra calculated using conventional TD-CI do not converge with increasing AO basis set size, while the zero lifetime and heuristic lifetime models converge to qualitatively similar spectra, with implications for how best to apply bound state electronic structure methods to simulate HHG. The origin of HHG spectral features including the cutoff and extent of interference between peaks is uncovered by separating field-induced coupling between different types of levels (ground state, bound Rydberg levels, and continuum) in the simulated electronic dynamics. Thus the origin of deviations between the predictions of the semi-classical three step model and the full simulation can be associated with particular physical contributions, which helps to explain both the successes and the limitations of the three step model.

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

Kinetics of the Hydrogen Atom Abstraction Reactions from 1-Butanol by Hydroxyl Radical: Theory Matches Experiment and More

In the present work, we study the H atom abstraction reactions by hydroxyl radical at all five sites of 1-butanol. Multistructural variational transition state theory (MS-VTST) was employed to estimate the five thermal rate constants. MS-VTST utilizes a multifaceted dividing surface that accounts for the multiple conformational structures of the transition state, and we also include all the structures of the reactant molecule. The vibrational frequencies and minimum energy paths (MEPs) were computed using the M08-HX/MG3S electronic structure method. The required potential energy surfaces were obtained implicitly by direct dynamics employing interpolated variational transition state theory with mapping (IVTST-M) using a variational reaction path algorithm. The M08-HX/MG3S electronic model chemistry was then used to calculate multistructural torsional anharmonicity factors to complete the MS-VTST rate constant calculations. The results indicate that torsional anharmonicity is very important at higher temperatures, and neglecting it would lead to errors of 26 and 32 at 1000 and 1500 K, respectively. Our results for the sums of the site-specific rate constants agree very well with the experimental values of Hanson and co-workers at 896–1269 K and with the experimental results of Campbell et al. at 292 K, but slightly less well with the experiments of Wallington et al., Nelson et al., and Yujing and Mellouki at 253–372 K; nevertheless, the calculated rates are within a factor of 1.61 of all experimental values at all temperatures. Finally, this gives us confidence in the site-specific values, which are currently inaccessible to experiment.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2023

Elucidating the Role of Hydrogen Bonding in the Optical Spectroscopy of the Solvated Green Fluorescent Protein Chromophore: Using Machine Learning to Establish the Importance of High-Level Electronic Structure

Hydrogen bonding interactions with chromophores in chemical and biological environments play a key role in determining their electronic absorption and relaxation processes, which are manifested in their linear and multidimensional optical spectra. For chromophores in the condensed phase, the large number of atoms needed to simulate the environment has traditionally prohibited the use of high-level excited-state electronic structure methods. By leveraging transfer learning, we show how to construct machine-learned models to accurately predict the high-level excitation energies of a chromophore in solution from only 400 high-level calculations. Here, we show that when the electronic excitations of the green fluorescent protein chromophore in water are treated using EOM-CCSD embedded in a DFT description of the solvent the optical spectrum is correctly captured and that this improvement arises from correctly treating the coupling of the electronic transition to electric fields, which leads to a larger response upon hydrogen bonding between the chromophore and water.

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

Group-III quantum defects in diamond are stable spin-1 color centers

Color centers in diamond have emerged as leading solid-state “artificial atoms” for a range of quantum technologies, from quantum sensing to quantum networks. Concerted research activities are now underway to identify new color centers that combine stable spin and optical properties of the nitrogen vacancy (NV – ) with the spectral stability of the silicon vacancy (SiV – ) centers in diamond, with recent research identifying other group-IV color centers with superior properties. In this paper, we investigate a class of diamond quantum emitters from first principles, the group-III color centers, which we show to be thermodynamically stable in a spin-1, electric-field-insensitive structure. Further, from ab initio electronic structure methods, we characterize the product Jahn-Teller (pJT) effect present in the excited-state manifold of these group-III color centers, where we capture symmetry-breaking distortions associated with strong electron-phonon coupling. These predictions can guide experimental identification of group-III vacancy centers and their use in applications in quantum information science and technology.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Photodissociation of OCS: Deviations between theory and experiment, and the importance of higher order correlation effects

The photodissociation of carbonyl sulfide (OCS) was investigated theoretically in a series of studies by Schmidt and co-workers. Initial studies [J. A. Schmidt, M. S. Johnson, G. C. McBane, and R. Schinke, J. Chem. Phys. 136, 131101 (2012); J. A. Schmidt, M. S. Johnson, G. C. McBane, and R. Schinke, J. Chem. Phys. 137, 054313 (2012)] found photodissociation in the first UV-band to occur mainly by excitation of the 2{sup 1}A{sup ′} (A) excited state. However, in a later study [G. C. McBane, J. A. Schmidt, M. S. Johnson, and R. Schinke, J. Chem. Phys. 138, 094314 (2013)] it was found that a significant fraction of photodissociation must occur by excitation of 1{sup 1}A{sup ″} (B) excited state to explain the product angular distribution. The branching between excitation of the A and B excited states is determined by the magnitude of the transition dipole moment vectors in the Franck-Condon region. This study examines the sensitivity of these quantities to changes in the employed electronic structure methodology. This study benchmarks the methodology employed in previous studies against highly correlated electronic structure methods (CC3 and MRAQCC) and provide evidence in support of the picture of the OCS photodissociation process presented in [G. C. McBane, J. A. Schmidt, M. S. Johnson, and R. Schinke, J. Chem. Phys. 138, 094314 (2013)] showing that excitation of A and B electronic states both contribute significantly to the first UV absorption band of OCS. In addition, this study presents evidence in support of the assertion that the A state potential energy surface employed in previous studies underestimates the energy at highly bent geometries (γ ∼ 70°) leading to overestimated rotational energy in the product CO.

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

Anti-site disorder and improved functionality of Mn₂NiX (X = Al, Ga, In, Sn) inverse Heusler alloys

Recent first-principles calculations have predicted Mn₂NiX (X = Al, Ga, In, Sn) alloys to be magnetic shape memory alloys. Moreover, experiments on Mn₂NiGa and Mn₂NiSn suggest that the alloys deviate from the perfect inverse Heusler arrangement and that there is chemical disorder at the sublattices with tetrahedral symmetry. In this work, we investigate the effects of such chemical disorder on phase stabilities and magnetic properties using first-principles electronic structure methods. We find that except Mn₂NiAl, all other alloys show signatures of martensitic transformations in presence of anti-site disorder at the sublattices with tetrahedral symmetry. This improves the possibilities of realizing martensitic transformations at relatively low fields and the possibilities of obtaining significantly large inverse magneto-caloric effects, in comparison to perfect inverse Heusler arrangement of atoms. We analyze the origin of such improvements in functional properties by investigating electronic structures and magnetic exchange interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2010

X-ray resonant magnetic scattering and x-ray magnetic circular dichroism branching ratios, L[subscript 3] / L[subscript 2], for heavy rare earths

In this study we have used first principles electronic structure methods to investigate the detailed contributions to the L{sub 3}/L{sub 2} branching ratio in the heavy rare earth elements. The calculations use the full potential, relativistic, linear augmented plane wave method with the LSDA+U approach for consideration of the local 4f electronic orbitals. With no spin orbit coupling (SOC) in the conducting bands, and with the same radial function for the 2p{sub 3/2} and 2p{sub 1/2} core states, the branching ratio (BR) is exactly 1:-1 for the x-ray magnetic circular dichroism spectra of the ferromagnetic heavy rare earth metals. However, with full SOC the BR ranges from 1.5 to 6.0 in going from Gd to Er. The energy and spin dependence of the 5d radial functions are important. The results point to problems with modified atomic models which have been proposed to explain the BR. Recent x-ray resonant magentic scattering experiments on (Gd,Tb,Dy,Ho,Er,Tm)Ni{sub 2}Ge{sub 2} are discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Corresponding Orbitals Derived from Periodic Bloch States for Electron Transfer Calculations of Transition Metal Oxides

An approach for modeling electron transfer in solids and at surfaces of iron - (oxy- hydr)oxides and other redox active solids has been developed for electronic structure methods (i.e., plane-wave DFT) capable of performing ab initio molecular dynamics eciently with periodic cells and large system sizes while at the same time being accurate enough to be used in the estimation of the electron-transfer coupling matrix element, VAB, and the electron transfer transmission factor, el. This method is an extension of the valence bond theory ET method for molecules and clusters, imple- mented by Dupuis and co-workers and used extensively by Rosso and co-workers, in which scaled corresponding orbitals derived from the Bloch states are used to calculate the o-diagonal matrix elements HAB and SAB. A key development of the present work is the formulation of algorithms to improve of the accuracy of integration of the exact exchange integral in periodic boundary conditions. This method is demonstrated on model systems for electron small polaron transfer in iron-(oxyhydr)oxides, including bare Fe2+{Fe3+ ions, and in [Fe3+(OH2)2(OH)2)]n+ n chains representing the common edge-sharing Fe octahedral motif in these materials.

Bylaska, Eric J.↗
DOE OSTI2021

Direct coherent switching with decay of mixing for intersystem crossing dynamics of thioformaldehyde: The effect of decoherence

We evaluate the effect of electronic decoherence on intersystem crossing in the photodynamics of thioformaldehyde. First, we show that the state-averaged complete-active-space self-consistent field electronic structure calculations with a properly chosen active space of 12 active electrons in 10 active orbitals can predict the potential energy surfaces and the singlet–triplet spin–orbit couplings quite well for CH 2 S, and we use this method for direct dynamics by coherent switching with decay of mixing (CSDM). We obtain similar dynamical results with CSDM or by adding energy-based decoherence to trajectory surface hopping, with the population of triplet states tending to a small steady-state value over 500 fs. Without decoherence, the state populations calculated by the conventional trajectory surface hopping method or the semiclassical Ehrenfest method gradually increase. This difference shows that decoherence changes the nature of the results not just quantitatively but qualitatively.

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

Guest-Host Interactions in Clathrate Hydrates: Benchmark MP2 and CCSD(T)/CBS Binding Energies of CH4, CO2 and H2S in (H2O)20 Cages

We present benchmark binding energies of naturally occurring gas molecules CH4, CO2, and H2S in the small cage, namely the pentagonal dodecahedron (512) (H2O)20, which is one of the constituent cages of the 3 major lattices (structures I, II and H) of clathrate hydrates. These weak interactions require higher levels of electron correlation and converge slowly with increasing basis set to the Complete Basis Set (CBS) limit, necessitating the use of large basis sets up to the augcc- pV5Z and subsequent correction for Basis Set Superposition Error (BSSE). For the host hollow (H2O)20 cages we have identified a most stable isomer with binding energy of -200.8 ± 2.1 kcal/mol at the CCSD(T)/CBS limit (-199.2 ± 0.5 kcal/mol at the MP2/CBS limit). Additionally, we report converged second order Moller-Plesset (MP2) CBS binding energies for the encapsulation of guests in the (H2O)20 cage of -4.3 ± 0.1 for CH4@(H2O)20, -6.6 ± 0.1 for CO2@(H2O)20 and -8.5 ± 0.1 kcal/mol for H2S@(H2O)20, respectively. For CH4@(H2O)20, exhibiting the weakest encapsulation affinity among the three, we report CCSD(T)/aug-cc-pVTZ binding energies and, based on them, a CCSD(T)/CBS estimate of -4.75 ± 0.1 kcal/mol. To the best of our knowledge, the CCSD(T)/aug-cc-pVTZ calculation for CH4@(H2O)20 is the largest one reported to date (168 valence electrons, 1978 basis functions and the correlation of 84 doubly occupied and 1873 virtual orbitals) and required a scalable implementation of the (T) module on 6144 nodes (350208 cores) of the “Cori” supercomputer at the National Energy Research Supercomputing Center (NERSC) for a total execution time of 195 minutes (for the (T) part). These efficient scalable implementations of highly correlated methods offer the capability to obtain long-lasting benchmarks of intermolecular interactions in complex systems. They also provide a path towards parametrizing classical potentials needed to study the dynamical and transport properties in these complex systems as well as assess the accuracy of lower scaling electronic structure methods such as Density Functional Theory (DFT) and MP2 including its spin-biased variants.

Heindel, Joseph↗
DOE OSTI2018

Development of formalisms based on locally coupled open subsystems for calculations in molecular electronic structure and dynamics

Quantum embedding theories model a collection of interacting molecules as a set of subsystems, where each can be treated with a particular electronic structure method (wave function or density functional theory, for example); these theories can lead to computationally efficient and accurate algorithms. Motivated by challenges in the field, we previously described a formalism (that models two subsystems), which we call “locally coupled open subsystems” (LCOS), for the computation of ground-state energies, fractional electron-occupation numbers of the subsystems, and the size-consistent limit of subsystem dissociation. In this work we present the full (nonrelativistic) LCOS theory and the following extensions of our previous work: (i) the framework to study systems composed of multiple subsystems and a procedure to spin-adapt the auxiliary wave function that describes the partitioned system, so that its spin state matches that of the real system of interest; (ii) potential functionals and ideas to employ machine learning for the computation of ground-state densities and energies; (iii) formulation of two LCOS ground-state approaches where the fragments are assigned Kohn-Sham wave functions; and (iv) a time-dependent (TD) extension of these two ground-state formalisms in which the state of the subsystems evolves according to a unitary propagation; from this evolution we can extract TD electron populations of the fragments, for instance. We also discuss potential applications of the TD LCOS theory to linear photoabsorption and Raman spectroscopy. As a result, the developments presented in this work can lead to ground-state and TD electronic structure calculations where the computational scaling can be controlled, depending on the level of theory and the accuracy desired to model each one of the subsystems and their coupling.

74 ATOMIC AND MOLECULAR PHYSICS↗