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 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 OSTI2014

Spin-forbidden and spin-allowed cyclopropenone (c-H{sub 2}C{sub 3}O) formation in interstellar medium

Three proposed mechanisms of cyclopropenone (c-H{sub 2}C{sub 3}O) formation from neutral species are studied using high-level electronic structure methods in combination with nonadiabatic transition state and collision theories to deduce the likelihood of each reaction mechanism under interstellar conditions. The spin-forbidden reaction involving the singlet electronic state of cyclopenylidene (c-C{sub 3}H{sub 2}) and the triplet state of atomic oxygen is studied using nonadiabatic transition state theory to predict the rate constant for c-H{sub 2}C{sub 3}O formation. The spin-allowed reactions of c-C{sub 3}H{sub 2} with molecular oxygen and acetylene with carbon monoxide were also investigated. The reaction involving the ground electronic states of acetylene and carbon monoxide has a very large reaction barrier and is unlikely to contribute to c-H{sub 2}C{sub 3}O formation in interstellar medium. The spin-forbidden reaction of c-C{sub 3}H{sub 2} with atomic oxygen, despite the high probability of nonadiabatic transition between the triplet and singlet states, was found to have a very small rate constant due to the presence of a small (3.8 kcal mol{sup –1}) reaction barrier. In contrast, the spin-allowed reaction between c-C{sub 3}H{sub 2} and molecular oxygen is found to be barrierless, and therefore can be an important path to the formation of c-H{sub 2}C{sub 3}O molecule in interstellar environment.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2010

Free Energies for Degradation Reactions of 1,2,3-Trichloropropane from ab initio Electronic Structure Theory

Electronic structure methods were used to calculate the gas-phase and aqueous phase reaction energies for reductive dechlorination (i.e. hydrogenolysis), reductive Beta-elimination, dehydrochlorination, and nucleophilic substitution by OH- of 1,2,3-trichloropropane. The thermochemical properties Delta Hof(298.15K), So(298.15K,1 bar), and Delta GS(298.15K, 1 bar) were calculated by using ab initio electronic structure calculations, isodesmic reactions schemes, gas-phase entropy estimates, and continuum solvation models for 1,2,3-trichloropropane and several likely metabolites. On the basis of these thermochemical estimates, together with a Fe(II)/Fe(III) chemical equilibrium model for natural reducing environments, all of the reactions studied were predicted to be very favorable in the standard state and under a wide range of pH conditions. The most favorable reaction was reductive Beta-elimination (Delta Gorxn ≈ -32 kcal/mol), followed closely by reductive dechlorination (Delta Gorxn ≈ -27 kcal/mol), dehydrochlorination (Delta Gorxn ≈ -27kcal/mol), and nucleophilic substitution by OH- (Delta Gorxn ≈ -25 kcal/mol). For both reduction reactions studied, it was found that the first electron-transfer step, yielding the intermediate CH2-CHCl-CH2Cl , and CH2Cl-CH-CH2Cl species, was not favorable in the standard state (Delta Gorxn ≈ +15 kcal/mol) and was predicted to occur only at relatively high pH values. This result suggests that reduction by natural attenuation is unlikely.

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

QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion quantum Monte Carlo

We review recent advances in the capabilities of the open source ab initio Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for greater efficiency and reproducibility. The auxiliary field QMC (AFQMC) implementation has been greatly expanded to include k-point symmetries, tensor-hypercontraction, and accelerated graphical processing unit (GPU) support. These scaling and memory reductions greatly increase the number of orbitals that can practically be included in AFQMC calculations, increasing the accuracy. Advances in real space methods include techniques for accurate computation of bandgaps and for systematically improving the nodal surface of ground state wavefunctions. Results of these calculations can be used to validate application of more approximate electronic structure methods, including GW and density functional based techniques. To provide an improved foundation for these calculations, we utilize a new set of correlation-consistent effective core potentials (pseudopotentials) that are more accurate than previous sets; these can also be applied in quantum-chemical and other many-body applications, not only QMC. These advances increase the efficiency, accuracy, and range of properties that can be studied in both molecules and materials with QMC and QMCPACK.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗