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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2019

Partition Density Functionals: Theory and Applications. Final Technical Report

Partition Density Functional Theory is a method to calculate the energy and density of a molecule via self-consistent calculations on isolated fragments. At the start of this project, Partition Density Functional Theory (PDFT) had only been applied to one-dimensional model systems of non-interacting electrons. Its most appealing feature (the reason why PDFT was born in the first place) was that it allowed one to construct an internally-consistent Density Functional Theory (DFT) valid for non-integer electron numbers to be used in the context of Chemical Reactivity Theory. Many formal properties of the quantities entering the theory, such as partition potentials, were unknown when the project started. The connections with other fragment-based electronic structure methods, and especially with embedding theories, were unclear, and the potential of PDFT for fixing problems of approximate exchange-correlation functionals was unanticipated. Supported by grant DE-SC0005031, we made progress in developing, understanding, implementing, and applying PDFT. We discovered its potential for applications in situations where approximate DFT fails, such as when stretching bonds or when calculating the energy of weakly-interacting systems, both cases of importance for energy research applications. We extended PDFT in several directions (spin-densities, current densities, external electric and magnetic fields), and carried out many calculations on diatomic molecules and small clusters to test ideas for approximate non-additive functionals. As a result, we were able to propose physically-motivated fragment-density approximations for both covalent and weak bonds. Because many applications will increasingly benefit from efficient and accurate electronic-structure calculations via density-embedding techniques, the research results enabled by this grant should prove useful for addressing complex problems in a wide variety of fields, from chemistry to molecular biology and materials engineering.

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

Low energy electron induced cytosine base release in 2′-deoxycytidine-3′-monophosphate via glycosidic bond cleavage: A time-dependent wavepacket study

Low energy electron (LEE) induced cytosine base release in a selected pyrimidine nucleotide, viz., 2′-deoxycytidine-3′-monophosphate is investigated using ab initio electronic structure methods and time dependent quantum mechanical calculations. It has been noted that the cytosine base scission is comparatively difficult process than the 3′ C–O bond cleavage from the lowest π{sup *} shape resonance in energy region <1 eV. This is mainly due to the high activation energy barrier associated with the electron transfer from the π{sup *} orbital of the base to the σ{sup *} orbital of the glycosidic N–C bond. In addition, the metastable state formed after impinging LEE (0–1 eV) has very short lifetime (10 fs) which may decay in either of the two competing auto-detachment or dissociation process simultaneously. On the other hand, the selected N–C mode may cleave to form the cytosine base anion at higher energy regions (>2 eV) via tunneling of the glycosidic bond. Resonance states generated within this energy regime will exist for a duration of ∼35–55 fs. Comparison of salient features of the two dissociation events, i.e., 3′ C–O single strand break and glycosidic N–C bond cleavage in 3′-dCMPH molecule are also provided.

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

Modeling Molecules under Pressure with Gaussian Potentials

The computational modeling of molecules under high pressure is a growing research area that augments experimental high-pressure chemistry. Here, a new electronic structure method for modeling atoms and molecules under pressure, Gaussians On Surface Tesserae Simulate HYdrostatic Pressure (GOSTSHYP) approach, is introduced. In this method, a set of Gaussian potentials is distributed evenly on the van der Waals surface of the investigated chemical system, leading to a compression of the electron density and the atomic scaffold. Since no parameters other than pressure need to be specified, GOSTSHYP allows straightforward geometry optimizations and ab initio molecular dynamics simulations of chemical systems under pressure for nonexpert users. In this work, calculated energies, bond lengths, and dipole moments under pressure fall within the range of established computational methods for high-pressure chemistry. A Diels-Alder reaction and the cyclotrimerization of acetylene showcase the ability of GOSTSHYP to model pressure-induced chemical reactions. The connection to mechanochemistry is pointed out.

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

The charge-asymmetric nonlocally determined local-electric (CANDLE) solvation model

Many important applications of electronic structure methods involve molecules or solid surfaces in a solvent medium. Since explicit treatment of the solvent in such methods is usually not practical, calculations often employ continuum solvation models to approximate the effect of the solvent. Previous solvation models either involve a parametrization based on atomic radii, which limits the class of applicable solutes, or based on solute electron density, which is more general but less accurate, especially for charged systems. We develop an accurate and general solvation model that includes a cavity that is a nonlocal functional of both solute electron density and potential, local dielectric response on this nonlocally determined cavity, and nonlocal approximations to the cavity-formation and dispersion energies. The dependence of the cavity on the solute potential enables an explicit treatment of the solvent charge asymmetry. With four parameters per solvent, this “CANDLE” model simultaneously reproduces solvation energies of large datasets of neutral molecules, cations, and anions with a mean absolute error of 1.8 kcal/mol in water and 3.0 kcal/mol in acetonitrile.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Ordered structures based on self-organized Au and CdSe nanoparticles

Methods for obtaining cadmium selenide and gold nanoparticles have been developed. The sizes of the nanoparticles are determined and the morphology, structure, and chemical composition of these nanoparticles and their ensembles are studied by a complex of structural methods: electron diffraction, X-ray diffraction, energy-dispersive X-ray analysis, high-resolution transmission electron microscopy, and small-angle X-ray scattering. Gold nanoparticles are mainly spherical and have an average size of 10 nm. They are single-phase and have an fcc crystal structure. Samples of synthesized CdSe nanoparticles contain monodisperse spherical particles 12 nm in size with a wurtzite structure. The deposition of nanoparticles on a carbon substrate is accompanied by their self-organization into a closely packed two-dimensional structure with a pronounced texture in which all nanoparticles are oriented in the [001] direction perpendicularly to the carbon substrate plane.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
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 OSTI2020

Computational Methods in Heterogeneous Catalysis

The unprecedented ability of computations to probe atomic-level details of catalytic systems holds immense promise for the fundamentals-based bottom-up design of novel heterogeneous catalysts, which are at the heart of the chemical and energy sectors of industry. Here, we critically analyze recent advances in computational heterogeneous catalysis. First, we will survey the progress in electronic structure methods and atomistic catalyst models employed, which have enabled the catalysis community to build increasingly intricate, realistic, and accurate models of the active sites of supported transition-metal catalysts. We then review developments in microkinetic modeling, specifically mean-field microkinetic models and kinetic Monte Carlo simulations, which bridge the gap between nanoscale computational insights and macroscale experimental kinetics data with increasing fidelity. Here, we finally review the advancements in theoretical methods for accelerating catalyst design and discovery. Throughout the review, we provide ample examples of applications, discuss remaining challenges, and provide our outlook for the near future.

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

Combining Dissociative Ionization Pump–Probe Spectroscopy and ab initio Calculations to Interpret Dynamics and Control Through Conical Intersections

Nonadiabatic processes play an important role in molecular dynamics, and understanding these processes better can help interpret and guide control over molecules. We are using high level electronic structure calculations in combination with intense, shaped, ultrafast laser pulses to study excited state dynamics in the nucleic acid bases, cytosine and uracil. These molecules have very short excited state lifetimes as they relax radiationless through conical intersections after absorption of UV radiation. The presence of more than one relaxation pathway provides the possibility to control which pathway can be involved in the dynamics. In our approach the molecules were excited using ultrafast laser pulses in the deep UV and then probed with strong field near infrared pulses which ionize and dissociate the molecules. Key to this approach is the fact that different fragments exhibit different dynamics and we can correlate these fragments, and their associated dynamics, to the various pathways involved in the neutral dynamics. Multiconfigurational electronic structure methods were used to calculate potential energy surfaces of the neutral and ionic states involved in the dynamics. Calculating mechanisms for fragmentation in the ion enables us to relate specific fragments to different neutral pathways, and use them as signatures to follow the dynamics. Possibilities for control are also discussed.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

How Well Can the M06 Suite of Functionals Describe the Electron Densities of Ne, Ne 6+ , and Ne 8+ ?

The development of better approximations to the exact exchange-correlation functional is essential to the accuracy of density functionals. A recent study indicated that functionals with few parameters offer more accurate electron densities than recently developed many-parameter functionals for light closed-shell atomic systems. In this study, we calculated electron densities, their gradients, and Laplacians of Ne, Ne 6+ , and Ne 8+ using 19 electronic structure methods, and we compared them to the CCSD reference results. Two basis sets, namely, aug-cc-pωCV5Z and aug-cc-pV5Z, are utilized in the calculations. We found that the choice of basis set has a significant impact on the errors and rankings of some of the selected methods. The errors of electron densities, their gradients, and Laplacians calculated with the aug-cc-pV5Z basis set are substantially reduced, especially for Minnesota density functionals, as compared to the results using the aug-cc-pωCV5Z basis set (a larger basis set utilized in earlier work (Medvedev et al. Science2017, 355, 49–52)). The rankings of the M06 suite of functionals among the 19 methods are greatly improved with the aug-cc-pV5Z basis set. Furthermore, the performances of the HSE06, BMK, MN12-L, and MN12-SX functionals are also improved with the aug-cc-pV5Z basis set. The M06 suite of functionals is capable of providing accurate electron densities, gradients, and Laplacians using the aug-cc-pV5Z basis set, and thus it is suitable for a wide range of applications in chemistry and physics.

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

Computational understanding of Li-ion batteries

Over the last two decades, computational methods have made tremendous advances, and today many key properties of lithium-ion batteries can be accurately predicted by first principles calculations. For this reason, computations have become a cornerstone of battery-related research by providing insight into fundamental processes that are not otherwise accessible, such as ionic diffusion mechanisms and electronic structure effects, as well as a quantitative comparison with experimental results. The aim of this review is to provide an overview of state-of-the-art ab initio approaches for the modelling of battery materials. Here, we consider techniques for the computation of equilibrium cell voltages, 0-Kelvin and finite-temperature voltage profiles, ionic mobility and thermal and electrolyte stability. The strengths and weaknesses of different electronic structure methods, such as DFT+U and hybrid functionals, are discussed in the context of voltage and phase diagram predictions, and we review the merits of lattice models for the evaluation of finite-temperature thermodynamics and kinetics. With such a complete set of methods at hand, first principles calculations of ordered, crystalline solids, i.e., of most electrode materials and solid electrolytes, have become reliable and quantitative. However, the description of molecular materials and disordered or amorphous phases remains an important challenge. We highlight recent exciting progress in this area, especially regarding the modelling of organic electrolytes and solid–electrolyte interfaces.

25 ENERGY STORAGE↗
DOE OSTI2015

Short-range stabilizing potential for computing energies and lifetimes of temporary anions with extrapolation methods

The energy of a temporary anion can be computed by adding a stabilizing potential to the molecular Hamiltonian, increasing the stabilization until the temporary state is turned into a bound state, and then further increasing the stabilization until enough bound state energies have been collected so that these can be extrapolated back to vanishing stabilization. The lifetime can be obtained from the same data, but only if the extrapolation is done through analytic continuation of the momentum as a function of the square root of a shifted stabilizing parameter. This method is known as analytic continuation of the coupling constant, and it requires—at least in principle—that the bound-state input data are computed with a short-range stabilizing potential. In the context of molecules and ab initio packages, long-range Coulomb stabilizing potentials are, however, far more convenient and have been used in the past with some success, although the error introduced by the long-rang nature of the stabilizing potential remains unknown. Here, we introduce a soft-Voronoi box potential that can serve as a short-range stabilizing potential. The difference between a Coulomb and the new stabilization is analyzed in detail for a one-dimensional model system as well as for the {sup 2}Π{sub u} resonance of CO{sub 2}{sup −}, and in both cases, the extrapolation results are compared to independently computed resonance parameters, from complex scaling for the model, and from complex absorbing potential calculations for CO{sub 2}{sup −}. It is important to emphasize that for both the model and for CO{sub 2}{sup −}, all three sets of results have, respectively, been obtained with the same electronic structure method and basis set so that the theoretical description of the continuum can be directly compared. The new soft-Voronoi-box-based extrapolation is then used to study the influence of the size of diffuse and the valence basis sets on the computed resonance parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗