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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Theory of the reaction dynamics of small molecules on metal surfaces

The objective of this project has been to develop realistic theoretical models for gas-surface interactions, with a focus on processes important in heterogeneous catalysis. The dissociative chemisorption of a molecule on a metal is a key step in many catalyzed reactions, and is often the rate-limiting step. We have explored the dissociative chemisorption of H 2 , H 2 O and CH 4 on a variety of metal surfaces. Most recently, our extensive studies of methane dissociation on Ni and Pt surfaces have fully elucidated its dependence on translational energy, vibrational state and surface temperature, providing the first accurate comparisons with experimental data. We have explored Eley-Rideal and hot atom reactions of H atoms with H- and C-covered metal surfaces. H atom interactions with graphite have also been explored, including both sticking and Eley-Rideal recombination processes. Again, our methods made it possible to explain several experiments studying these reactions. The sticking of atoms on metal surfaces has also been studied. To help elucidate the experiments that study these processes, we examine how the reaction dynamics depend upon the nature of the molecule-metal interaction, as well as experimental variables such as substrate temperature, beam energy, angle of impact, and the internal states of the molecules. Electronic structure methods based on Density Functional Theory are used to compute each molecule-metal potential energy surface. Both time-dependent quantum scattering techniques and quasi-classical methods are used to examine the reaction or scattering dynamics. Much of our effort has been directed towards developing improved quantum methods that can accurately describe reactions, as well as include the effects of substrate temperature (lattice vibration).

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2010

Combinatorial–Computational–Chemoinformatics (C³) Approach to Finding and Analyzing Low-Energy Tautomers

Finding the most stable tautomer or a set of low-energy tautomers of molecules is critical in many aspects of molecular modelling or virtual screening experiments. Enumeration of low-energy tautomers of neutral molecules in the gas-phase or typical solvents can be performed by applying available organic chemistry knowledge. This kind of enumeration is implemented in a number of software packages and it is relatively reliable. However, in esoteric cases such as charged molecules in uncommon, non-aqueous solvents there is simply not enough available knowledge to make reliable predictions of low energy tautomers. Over the last few years we have been developing an approach to address the latter problem and we successfully applied it to discover the most stable anionic tautomers of nucleic acid bases that might be involved in the process of DNA damage by low-energy electrons and in charge transfer through DNA. The approach involves three steps: (1) combinatorial generation of a library of tautomers, (2) energy-based screening of the library using electronic structure methods, and (3) analysis of the information generated in step (2). In steps 1–3 we employ combinatorial, computational and chemoinformatics techniques, respectively. Therefore, this hybrid approach is named “Combinatorial*Computational*Chemoinformatics”, or just abbreviated as C³ (or C-cube) approach. This article summarizes our developments and most interesting methodological aspects of the C³ approach. It can serve as an example how to identify the most stable tautomers of molecular systems for which common chemical knowledge had not been sufficient to make definite predictions.

59 BASIC BIOLOGICAL SCIENCES↗
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 OSTI2020

TeraChem: A graphical processing unit-accelerated electronic structure package for large-scale ab initio molecular dynamics

TeraChem was born in 2008 with the goal of providing fast on-the-fly electronic structure calculations to facilitate ab initio molecular dynamics studies of large biochemical systems such as photoswitchable proteins and multichromophoric antenna complexes. Originally developed for videogaming applications, graphics processing units (GPUs) offered a low-cost parallel computer architecture that became more accessible for general-purpose GPU computing with the release of CUDA in 2007. The evaluation of the electron repulsion integrals (ERIs) is a major bottleneck in electronic structure codes and provides an attractive target for acceleration on GPUs. Thus, highly efficient routines for evaluation of and contractions between the ERIs and density matrices were implemented in TeraChem. Here, electronic structure methods were developed and implemented to leverage these integral contraction routines, resulting in the first quantum chemistry package designed from the ground up for GPUs. This GPU acceleration makes TeraChem capable of performing large-scale ground and excited state calculations in the gas and condensed phase. Today, TeraChem's speed forms the basis for a suite of quantum chemistry applications, including optimization and dynamics of proteins, automated and interactive chemical discovery tools, and large-scale nonadiabatic dynamics simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2014

Simplified continuum solvent model with a smooth cavity based on volumetric data

We present a continuum solvent model (CSM) with a smooth cavity for the application in grid-based electronic structure methods. The cavity is identified with the inherently smooth distribution function of a binary mixture at infinite dilution. We obtain a cavity model based on atomic van der Waals radii and one free parameter controlling the overall size. This single parameter is sufficient to adequately reproduce experimental partial molar volumes. The CSM based on this cavity is of similar accuracy in the prediction of aqueous solvation Gibbs energies of small neutral molecules and ions as other CSMs with a smooth cavity. We apply the model to systems in non-aqueous solution, i.e., spiropyran/merocyanin energetics, a proton transfer reaction in dimethyl sulfoxide, and the electrostatic screening of charged gold clusters in an ionic liquid.

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

Potential energy surfaces for the HBr{sup +} + CO{sub 2} → Br + HOCO{sup +} reaction in the HBr{sup +} {sup 2}Π{sub 3/2} and {sup 2}Π{sub 1/2} spin-orbit states

Quantum mechanical (QM) + molecular mechanics (MM) models are developed to represent potential energy surfaces (PESs) for the HBr{sup +} + CO{sub 2} → Br + HOCO{sup +} reaction with HBr{sup +} in the {sup 2}Π{sub 3/2} and {sup 2}Π{sub 1/2} spin-orbit states. The QM component is the spin-free PES and spin-orbit coupling for each state is represented by a MM-like analytic potential fit to spin-orbit electronic structure calculations. Coupled-cluster single double and perturbative triple excitation (CCSD(T)) calculations are performed to obtain “benchmark” reaction energies without spin-orbit coupling. With zero-point energies removed, the “experimental” reaction energy is 44 ± 5 meV for HBr{sup +}({sup 2}Π{sub 3/2}) + CO{sub 2} → Br({sup 2}P{sub 3/2}) + HOCO{sup +}, while the CCSD(T) value with spin-orbit effects included is 87 meV. Electronic structure calculations were performed to determine properties of the BrHOCO{sup +} reaction intermediate and [HBr⋯OCO]{sup +} van der Waals intermediate. The results of different electronic structure methods were compared with those obtained with CCSD(T), and UMP2/cc-pVTZ/PP was found to be a practical and accurate QM method to use in QM/MM direct dynamics simulations. The spin-orbit coupling calculations show that the spin-free QM PES gives a quite good representation of the shape of the PES originated by {sup 2}Π{sub 3/2}HBr{sup +}. This is also the case for the reactant region of the PES for {sup 2}Π{sub 1/2} HBr{sup +}, but spin-orbit coupling effects are important for the exit-channel region of this PES. A MM model was developed to represent these effects, which were combined with the spin-free QM PES.

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 OSTI2018

Design of high-strength refractory complex solid-solution alloys

Nickel-based superalloys and near-equiatomic high-entropy alloys containing molybdenum are known for higher temperature strength and corrosion resistance. Yet, complex solid-solution alloys offer a huge design space to tune for optimal properties at slightly reduced entropy. For refractory Mo-W-Ta-Ti-Zr, we showcase KKR electronic structure methods via the coherent-potential approximation to identify alloys over five-dimensional design space with improved mechanical properties and necessary global (formation enthalpy) and local (short-range order) stability. Deformation is modeled with classical molecular dynamic simulations, validated from our first-principle data. We predict complex solid-solution alloys of improved stability with greatly enhanced modulus of elasticity (3× at 300 K) over near-equiatomic cases, as validated experimentally, and with higher moduli above 500 K over commercial alloys (2.3× at 2000 K). We also show that optimal complex solid-solution alloys are not described well by classical potentials due to critical electronic effects.

36 MATERIALS SCIENCE↗
DOE OSTI2011

Photoisomerization mechanism of 4-methylpyridine explored by electronic structure calculations and nonadiabatic dynamics simulations

In the present paper, different electronic structure methods have been used to determine stationary and intersection structures on the ground (S{sub 0}) and {sup 1}{pi}{pi}* (S{sub 2}) states of 4-methylpyridine, which is followed by adiabatic and nonadiabatic dynamics simulations to explore the mechanistic photoisomerization of 4-methylpyridine. Photoisomerization starts from the S{sub 2}({sup 1}{pi}{pi}*) state and overcomes a small barrier, leading to formation of the prefulvene isomer in the S{sub 0} state via a S{sub 2}/S{sub 0} conical intersection. The ultrafast S{sub 2}{yields} S{sub 0} nonradiative decay and low quantum yield for the photoisomerization reaction were well reproduced by the combined electronic structure calculation and dynamics simulation. The prefulvene isomer was assigned as a long-lived intermediate and suggested to isomerize to 4-methylpyridine directly in the previous study, which is not supported by the present calculation. The nonadiabatic dynamics simulation and electronic structure calculation reveal that the prefulvene isomer is a short-lived intermediate and isomerizes to benzvalene form very easily. The benzvalene form was predicted as the stable isomer in the present study and is probably the long-lived intermediate observed experimentally. A consecutive light and thermal isomerization cycle via Dewar isomer was determined and this cycle mechanism is different from that reported in the previous study. It should be pointed out that formation of Dewar isomer from the S{sub 2}({sup 1}{pi}{pi}*) state is not in competition with the isomerization to the prefulvene form. The Dewar structure observed experimentally may originate from other excited states.

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

Energies of the X- and L-valleys in In{sub 0.53}Ga{sub 0.47}As from electronic structure calculations

Several theoretical electronic structure methods are applied to study the relative energies of the minima of the X- and L-conduction-band satellite valleys of In{sub x}Ga{sub 1−x}As with x = 0.53. This III-V semiconductor is a contender as a replacement for silicon in high-performance n-type metal-oxide-semiconductor transistors. The energy of the low-lying valleys relative to the conduction-band edge governs the population of channel carriers as the transistor is brought into inversion, hence determining current drive and switching properties at gate voltages above threshold. The calculations indicate that the position of the L- and X-valley minima are ∼1 eV and ∼1.2 eV, respectively, higher in energy with respect to the conduction-band minimum at the Γ-point.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

From one to three, exploring the rungs of Jacob’s ladder in magnetic alloys

Magnetic systems represent an important challenge for electronic structure methods, in particular Density Functional Theory (DFT), which uses a single determinant wavefunction. To assess the predictions obtained by DFT in this type of materials, we benchmark different exchange correlation functionals with respect to each other, and with respect to available experimental data, on two families of binary iron alloys which are metallic and magnetic. We climb three rungs in Jacob’s ladder of DFT (i) the local density approximation, (ii) the industry standard approximation due to Perdew, Burke and Ernzerhof, and the revised version for solids, PBEsol (iii) and finally a very accurate meta-GGA functional SCAN, which corresponds to the third rung. More than 350 structures in ferromagnetic and antiferromagnetic configurations were considered. We compare the Convex Hull, the calculated magnetic moment, crystal structure, formation energy and electronic gap if present. We conclude that none of the functionals work in all conditions: whereas PBE and PBEsol can give a fair description of the crystal structure and the energetics, SCAN strongly overestimates the formation energy – giving values which are at least twice as large as PBE (and experiment). Magnetic moments are better predicted by PBE as well. Finally, our results show that magnetic and strongly correlated materials are a tough litmus test for DFT, and that they represent the next frontier in the development of a truly universal exchange correlation functionals.

Romero, Aldo H.↗
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↗