Search NASA⌕ Search

SEARCH · Search NASA

Results for “Electronic Structure Calculations”

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 91 records · Page 5

Hybrid Quantum-Classical Neural Network for Calculating Ground State Energies of Molecules

We present a hybrid quantum-classical neural network that can be trained to perform electronic structure calculation and generate potential energy curves of simple molecules. The method is based on the combination of parameterized quantum circuits and measurements. With unsupervised training, the neural network can generate electronic potential energy curves based on training at certain bond lengths. To demonstrate the power of the proposed new method, we present the results of using the quantum-classical hybrid neural network to calculate ground state potential energy curves of simple molecules such as H 2 , LiH, and BeH 2 . The results are very accurate and the approach could potentially be used to generate complex molecular potential energy surfaces.

74 ATOMIC AND MOLECULAR PHYSICS↗

Competition between Hydrogen and Chalcogen Bonding in Homodimers of Chalcogen Hydrides (H 2 X) 2 , X = O, S, Se, Te

The structural and chemical bonding motifs manifested by the competition between hydrogen and chalcogen bonding in the homodimers of chalcogen hydrides (H 2 X) 2 , where X = O, S, Se, Te have been characterized using high-level electronic structure calculations and analysis of the electron density based on Quasi-atomic orbital (QUAO) and the Symmetry-adapted perturbation theory (SAPT) methods. The QUAO analysis clearly identifies a three-center interaction responsible for either hydrogen or chalcogen bonds: in the former, the σ-bond connecting the donor and hydrogen atom participating in the hydrogen bond interacts with the lone pair on the nucleophile acceptor via the hydrogen atom, while in the latter this same σ-bond interacts with the nucleophile lone pair via the donor chalcogen. The number of minimum energy structures increase dramatically from one for (H 2 O) 2 , three for (H 2 S) 2 , four for (H 2 Se) 2 , and finally six for (H 2 Te) 2 . Furthermore, the emergence of the chalcogen-bonded arrangements appears for (H 2 S) 2 with their subsequent energetic stabilization over the hydrogen-bonded minima manifesting in (H 2 Se) 2 and (H 2 Te) 2 . In particular, one of the (H 2 S) 2 , two of the (H 2 Se) 2 , and three of the (H 2 Te) 2 dimers are chalcogen bonded. Induction plays a small but important role in stabilizing hydrogen over chalcogen-bonded structures, while dispersion is more important for chalcogen bonds.

Chemical structure↗

Electronic structure and magnetism of pristine, defected, and strained $\mathrm{Ti_2N}$ $\mathrm{MX}$ene

Here, from first principles electronic structure calculations, we unravel the evolution of structural, electronic, and magnetic properties of pristine, defected, and strained titanium nitride MXene with different functional groups (-F, -O, -H, and -OH). The formation and cohesive energies reveal their chemical stability. The dynamical stability of Ti 2 N mono-layer is also confirmed by phonon calculations. The MAX phase and defect free functionalized MXenes are metallic except for oxygen terminated (Ti 2 NO 2 ) one which is 100% spin polarized half-metallic ferromagnet. The spin–orbit coupling significantly influences the bare MXene (Ti 2 N) to exhibit Dirac topology and band inversion near the high symmetry directions. The strain effect sways the Fermi level thereby shifting it towards lower energy state under compression and towards higher energy state under tensile strain in Ti 2 NH 2 . The Ti 2 NO 2 exhibits exotic electronic structure not only in pristine but also in strained and defected structures. Its half-metallic nature changes to semi-metallic under 1% compression and it is completely destroyed under 2% compression. In single vacancy defect, its band structure remarkably transforms from half-metallic to semi-conducting with large band gap in 12.5% Ti, weakly semi-conducting in 5.5% Ti, and semi-metallic in 12.5% O. The 25% N defect changes its half-metallic characteristic to metallic. Further, the 12.5% Co substitution preserves its half-metallic character, whereas Mn substitution allows it to convert half-metallic characteristic into weak semi-metallic characteristic preserving ferromagnetism. However, Cr substitution converts half-metallic ferromagnetic state to half-metallic anti-ferromagnetic state. The understanding made here on collective structural stability, and electronic band structure, and magnetic phenomena in novel 2D Ti 2 N derived MXenes open up their possibility in designing them for synthesis.

36 MATERIALS SCIENCE↗

Interaction of Th with H 0/–/+ : Combined Experimental and Theoretical Thermodynamic Properties

In this work, high-level electronic structure calculations of the lowlying energy electronic states for ThH, ThH – , and ThH + are reported and compared to experimental measurements. The inclusion of spin–orbit coupling is critical to predict the ground-state ordering as inclusion of spin–orbit switches the coupled-cluster CCSD(T) ordering of the two lowest energy states for ThH and ThH + . At the multireference spin–orbit SO-CASPT2 level, the ground states of ThH, ThH – , and ThH + are predicted to be the 2 Δ 3/2 , 3 Φ 2 , and 3 Δ 1 states, respectively. The adiabatic electron affinity is calculated to be 0.820 eV, and the vertical detachment energy is calculated to be 0.832 eV in comparison to an experimental value of 0.87 ± 0.02 eV. The observed ThH – photoelectron spectrum has many transitions, which approximately correlate with excitations of Th + and/or Th. The adiabatic ionization energy of ThH including spin–orbit corrections is calculated to be 6.181 eV. The natural bond orbital results are consistent with a significant contribution of the Th + H – ionic configuration to the bonding in ThH. The bond dissociation energies for ThH, ThH – , and ThH + using the Feller–Peterson–Dixon approach were calculated to be similar for all three molecules and lie between 259 and 280 kJ/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental and Computational Description of the Interaction of H and H – with U

In this work, the results of ab initio correlated molecular orbital theory electronic structure calculations for low-lying electronic states are presented for UH and UH – and compared to photoelectron spectroscopy measurements. The calculations were performed at the CCSD(T)/CBS and multireference CASPT2 including spin–orbit effects by the state interacting approach levels. The ground states of UH and UH – are predicted to be 4 Ι 9/2 and 5 Λ 6 , respectively. The spectroscopic parameters T e , r e , ω e , ω e x e , and Be were obtained, and potential energy curves were calculated for the low energy Ω states of UH. The calculated adiabatic electron affinity is 0.468 eV in excellent agreement with an experimental value of 0.462 ± 0.013 eV. The lowest vertical detachment energy was predicted to be 0.506 eV for the ground state, and the adiabatic ionization energy (IE) is predicted to be 6.116 eV. The bond dissociation energy (BDE) and heat of formation values of UH were obtained using the IE calculated at the Feller–Peterson–Dixon level. For UH, UH – , and UH + , the BDEs were predicted to be 225.5, 197.9, and 235.5 kJ/mol, respectively. The BDE for UH is predicted to be ~20% lower in energy than that for ThH. The analysis of the natural bond orbitals shows a significant U + H – ionic component in the bond of UH.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atomistic learning in the electronically grand-canonical ensemble

Abstract A strategy is presented for the machine-learning emulation of electronic structure calculations carried out in the electronically grand-canonical ensemble. The approach relies upon a dual-learning scheme, where both the system charge and the system energy are predicted for each image. The scheme is shown to be capable of emulating basic electrochemical reactions at a range of potentials, and coupling it with a bootstrap-ensemble approach gives reasonable estimates of the prediction uncertainty. The method is also demonstrated to accelerate saddle-point searches, and to extrapolate to systems with one to five water layers. We anticipate that this method will allow for larger length- and time-scale simulations necessary for electrochemical simulations.

36 MATERIALS SCIENCE↗

ThAu 2 – , ThAu 2 O – , and ThAuOH – anions: Photoelectron spectroscopic and theoretical characterization

The thorium–gold negative ions ThAu 2 – , ThAu 2 O – , and ThAuOH – have been observed and experimentally characterized by anion photoelectron spectroscopy. These experiments are accompanied by extensive ab initio electronic structure calculations using a relativistic composite methodology based primarily on coupled cluster singles and doubles with perturbative triples calculations. The theoretical electron affinities (EAs) at 0 K agree with the experimental adiabatic EAs to within 0.02 eV for all species. Two separate isomers were located in the calculations for ThAuOH – , and detachment from both of these appears to be present in the photoelectron spectrum. Excited electronic states of the neutral molecules are reported at the equation of motion-coupled cluster singles and doubles level of theory. Atomization energies and heats of formation are also calculated for each neutral species and have expected uncertainties of 3 and 4 kcal/mol, respectively. The σ bonds between Th and Au are determined by natural bond orbital analysis to consist of predominately sd hybrids on Th bonding with the Au 6s orbital. In order to investigate the correspondence between the bonding in Th–Au and Th–F molecules, a limited number of calculations were also carried out on most of the F-analogs of this study. These results demonstrate that Au does behave like F in these cases, although the Th–F σ bonds are much more ionic compared to Th–Au. This results in an EA for ThF 2 that is 10 kcal/mol smaller than that of ThAu 2 . The EA values for the Th(IV) species, i.e., ThX 2 O and ThXOH, only differed, however, by 3–4 kcal/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electronic structure modulation of MoS2 by substitutional Se incorporation and interfacial MoO3 hybridization: Implications of Fermi engineering for electrocatalytic hydrogen evolution and oxygen evolution

The design of earth-abundant electrocatalysts that can facilitate water splitting at low overpotentials, provide high current densities, and enable prolonged operational lifetimes is central to the production of sustainable fuels. The distinctive atomistic and electronic structure characteristics of the edges of MoS2 imbue high reactivity toward the hydrogen evolution reaction. MoS2 is nevertheless characterized by significantly high overpotentials as compared to platinum. Here, we demonstrate that modulation of the electronic structure of MoS2 through interfacial hybridization with MoO3 and alloying of selenium on the anion sublattice allows for systematic lowering of the conduction band edge and raising of the valence band edge, respectively. The former promotes enhanced electrocatalytic activity toward hydrogen evolution, whereas the latter promotes enhanced activity toward the oxygen evolution reaction. Such alloyed heterostructures prepared by sol-gel reactions and hydrothermal selenization expose a high density of edge sites. The alloyed heterostructures exhibit low overpotential, high current density, high turnover frequency, and prolonged operational lifetime. The mechanistic origins of catalytic activity have been established based on electronic structure calculations and x-ray absorption and emission spectroscopy probes of electronic structure, which suggest that interfacial hybridization at the MoO3 interface yields low-lying conduction band states that facilitate hydrogen adsorption. In contrast, shallow Se 4p-derived states give rise to a raised effective valence band maximum, which facilitates adsorption of oxygen intermediates and engenders a low overpotential for the oxygen evolution reaction. The findings illustrate the use of electronic structure modulation through interfacial hybridization and alloying to systematically improve electrocatalytic activity.

Parija, Abhishek↗

From Atom to Engine: Understanding Fundamental Effects of Structure on Combustion Using Tandem Experiment and Computation

The ability to tie structural features to a fuel candidate’s ignition properties provides a path for rational design of advantaged fuels that facilitate higher-efficiency, lower-emitting combustion in engines. The time evolution of a fuel’s radical population directly correlates with global combustion parameters, and because the generated radicals depend on the initial structure and reaction pathways of the fuel molecule, it follows that the molecular structure of a fuel has a direct impact on this cascade, and therefore salient ignition-tied fuel properties. This work aims to understand the impact of structural features in the context of global combustion properties for direct application to emerging engine technologies. Specifically, we highlight the ability to identify potential fuel blendstocks for specific engine strategies by examining and elucidating the radical cascade tied to ignition. A representative subset of alcohols possessing a variety of structural features (e.g. branching, chain length) was selected for study using a robust combination of experimental and computational methods. The approach begins with identifying short-lived radicals and likely abstraction sites upon pyrolysis in a microreactor coupled with photoionization mass spectrometry for direct detection of product species. In tandem, electronic structure calculations were performed to calculate the relevant potential energy surfaces for fuel pyrolysis as well as the subsequent O2 addition-isomerization reactions, which govern ignition behavior. Furthermore, the theoretical energetic and ro-vibrational data are used to calculate rate constants for direct ignition modeling to explore how structural effects impact NTC behavior. The microreactor data combined with electronic structure theory calculations provide a well-informed picture of the reactions relevant to ignition for each fuel in this study, as well as the most likely radicals generated under engine conditions. This thorough fundamental picture of fuel decomposition and ignition chemistry is used to correlate the observed ignition behavior to molecular structural effects.

advanced biofuels↗

Methods to Calculate Electronic Excited-State Dynamics for Molecules on Large Metal Clusters with Many States: Ensuring Fast Overlap Calculations and a Robust Choice of Phase

Here, we present an efficient set of methods for propagating excited-state dynamics involving a large number of configuration interaction singles (CIS) or Tamm-Dancoff approximation (TDA) single-reference excited states. Specifically, (i) following Head-Gordon et al., we implement an exact evaluation of the overlap of singly-excited CIS/TDA electronic states at different nuclear geometries using a biorthogonal basis and (ii) we employ a unified protocol for choosing the correct phase for each adiabat at each geometry. For many-electron systems, the combination of these techniques significantly reduces the computational cost of integrating the electronic Schrodinger equation and imposes minimal overhead on top of the underlying electronic structure calculation. As a demonstration, we calculate the electronic excited-state dynamics for a hydrogen molecule scattering off a silver metal cluster, focusing on high-lying excited states, where many electrons can be excited collectively and crossings are plentiful. Interestingly, we find that the high-lying, plasmon-like collective excitation spectrum changes with nuclear dynamics, highlighting the need to simulate non-adiabatic nuclear dynamics and plasmonic excitations simultaneously. In the future, the combination of methods presented here should help theorists build a mechanistic understanding of plasmon-assisted charge transfer and excitation energy relaxation processes near a nanoparticle or metal surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Validating first-principles molecular dynamics calculations of oxide/water interfaces with x-ray reflectivity data

Metal oxide/water interfaces play a crucial role in many electrochemical and photocatalytic processes, such as photoelectrochemical water splitting, the creation of fuel from sunlight, and electrochemical CO 2 reduction. First-principles electronic structure calculations can reveal unique insights into these processes, such as the role of the alignment of the oxide electronic energy levels with those of liquid water. An essential prerequisite for the success of such calculations is the ability to predict accurate structural models of these interfaces, which in turn requires careful experimental validation. Here we report a general, quantitative validation protocol for first-principles molecular dynamics simulations of oxide/aqueous interfaces. The approach makes direct comparisons of interfacial x-ray reflectivity (XR) signals from experimental measurements and those obtained from ab initio simulations with semilocal and van der Waals functionals. The protocol is demonstrated here for the case of the Al 2 O 3 (001)/water interface, one of the simplest oxide/water interfaces. We discuss the technical requirements needed for validation, including the choice of the density functional, the simulation cell size, and the optimal choice of the thermodynamic ensemble. Our results establish a general paradigm for the validation of structural models and interactions at solid/water interfaces derived from first-principles simulations. Furthermore, while there is qualitative agreement between the simulated structures and the experimental best-fit structure, direct comparisons of simulated and measured XR intensities show quantitative discrepancies that derive from both bulk regions (i.e., alumina and water) as well as the interfacial region, highlighting the need for accurate density functionals to properly describe interfacial interactions. Our results show that XR data are sensitive not only to the atomic structure (i.e., the atom locations) but also to the electron-density distributions in both the substrate and at the interface.

36 MATERIALS SCIENCE↗

Prediction of impact sensitivity, heat of formation and heat of explosion using atomic connectivity

In these proceedings we revisit a large collection of explosives and explosive descriptors with the goal of predicting impact sensitivity using only local atomic environments that can be deciphered from molecular SMILES strings as descriptors without utilizing empirically measured values or computationally expensive electronic structure calculations. From the original database of nearly 500 descriptors, removing empirically measured and electronic structure values decreased the number of descriptors to 135, which we reduced to 18 the most important descriptors using Random Forests. The condensed model predicted impact sensitivity with essentially the same accuracy as the existing, more complex model (R 2 = 0.788 and RSME = 0.312), while remaining applicable to all types of explosives (Peroxides, azides, C-Nitros, Nitroamines, Nitrate Esters, etc.). In addition to impact sensitivity, we proposed similar models to accurately predict values heat of formation (ΔH f ) and heat of explosion (Q), with R 2 = 0.966 and 0.916, respectively. In conclusion, the work in these proceedings allows for prediction of explosive performance and sensitivity with only chemical structure information and an estimate of density.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Mixed Precision Fermi-Operator Expansion on Tensor Cores from a Machine Learning Perspective

Here we present a second-order recursive Fermi-operator expansion scheme using mixed precision floating point operations to perform electronic structure calculations using tensor core units. A performance of over 100 teraFLOPs is achieved for half-precision floating point operations on Nvidia’s A100 tensor core units. The second-order recursive Fermi-operator scheme is formulated in terms of a generalized, differentiable deep neural network structure, which solves the quantum mechanical electronic structure problem. We demonstrate how this network can be accelerated by optimizing the weight and bias values to substantially reduce the number of layers required for convergence. We also show how this machine learning approach can be used to optimize the coefficients of the recursive Fermi-operator expansion to accurately represent the fractional occupation numbers of the electronic states at finite temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Scalable generalized screening for high-order terms in the many-body expansion: Algorithm, open-source implementation, and demonstration

The many-body expansion lies at the heart of numerous fragment-based methods that are intended to sidestep the nonlinear scaling of ab initio quantum chemistry, making electronic structure calculations feasible in large systems. In principle, inclusion of higher-order n-body terms ought to improve the accuracy in a controllable way, but unfavorable combinatorics often defeats this in practice and applications with n ≥ 4 are rare. Here, we outline an algorithm to overcome this combinatorial bottleneck, based on a bottom-up approach to energy-based screening. This is implemented within a new open-source software application (“Fragme∩t”), which is integrated with a lightweight semi-empirical method that is used to cull subsystems, attenuating the combinatorial growth of higher-order terms in the graph that is used to manage the calculations. This facilitates applications of unprecedented size, and we report four-body calculations in (H2O)64 clusters that afford relative energies within 0.1 kcal/mol/monomer of the supersystem result using less than 10% of the unique subsystems. We also report n-body calculations in (H2O)20 clusters up to n = 8, at which point the expansion terminates naturally due to screening. These are the largest n-body calculations reported to date using ab initio electronic structure theory, and they confirm that high-order n-body terms are mostly artifacts of basis-set superposition error.

Chemistry↗

Single crystal structure, vibrational spectroscopy, gas sorption and antimicrobial properties of a new inorganic acidic diphosphates material (NH 4 ) 2 Mg(H 2 P 2 O 7 ) 2 • 2 H 2 O

We report on the successful synthesis of diammonium magnesium dihydrogendiphosphate (V) dihydrate compound (NH 4 ) 2 Mg(H 2 P 2 O 7 ) 2 •2H 2 O using a wet chemical route. Single crystal X-ray diffraction analysis and micro Raman spectroscopy are employed to characterize the compound. We demonstrate, using a multidisciplinary approach, that this compound is ideal for carbon dioxide (CO 2 ) capture in addition to other anthropogenic gasses. We show here -from both an experimental as well as from a density functional theory (DFT) calculations routes- the potential for adopting this compound into domestic air-conditioning units (ACUs). From these experiments, the resistance to bacterial growth is also investigated, which is critical for the adoption of this compound in ACUs. Our compound exhibits a higher methane (CH 4 ) sorptivity as compared to CO 2 at 25 °C and 45 °C under pressures up to 50 bars. Furthermore, DFT electronic structure calculations are used to compute the main structural and electronic properties of the compound, taking into consideration the characteristics of the identified pores as a function of the progressive CO 2 vs. CH 4 loadings. Finally, the antibacterial assay reveals a strong antibacterial activity against the tested Gram-positive and Gram-negative bacteria, with a large zone of inhibition against the tested E. Coli, S. Aureus and K. Pneumonia

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Using a multistate mapping approach to surface hopping to predict the ultrafast electron diffraction signal of gas-phase cyclobutanone

Using the recently developed multistate mapping approach to surface hopping (multistate MASH) method combined with SA(3)-CASSCF(12,12)/aug-cc-pVDZ electronic structure calculations, the gas-phase isotropic ultrafast electron diffraction (UED) of cyclobutanone is predicted and analyzed. After excitation into the n-3s Rydberg state (S2), cyclobutanone can relax through two S2/S1 conical intersections, one characterized by compression of the CO bond and the other by dissociation of the α–CC bond. Subsequent transfer into the ground state (S0) is then achieved via two additional S1/S0 conical intersections that lead to three reaction pathways: α ring-opening, ethene/ketene production, and CO liberation. The isotropic gas-phase UED signal is predicted from the multistate MASH simulations, allowing for a direct comparison to the experimental data. This work, which is a contribution to the cyclobutanone prediction challenge, facilitates the identification of the main photoproducts in the UED signal and thereby emphasizes the importance of dynamics simulations for the interpretation of ultrafast experiments.

Chemistry↗

Can Third-Body Stabilization of Bimolecular Collision Complexes in Cold Molecular Clouds Happen?

For more than half a century, networks of radiative association, dissociative recombination, and bimolecular reactions have been postulated to drive the low-temperature chemistry of cold molecular clouds. Third-body stabilizations of collision complexes have been assumed to be 'irrelevant' due to short lifetime of such complexes. Here, we conduct crossed molecular beam studies of ground state atomic silicon with diacetylene in combination with electronic structure calculations and microcanonical kinetics models operating under cold molecular cloud conditions. Our combined experimental, electronic structure, and microcanonical kinetics modelling investigations provide compelling evidence that three-body collisions of molecular hydrogen with long-lived reaction intermediates accessed through intersystem crossing are prevalent deep inside molecular clouds. This concept might be exportable to reactions involving polycyclic aromatic hydrocarbons thus affording a versatile machinery to complex organics via third-body stabilizations of bimolecular collision complexes deep inside cold molecular clouds.

astrochemistry↗