Search NASA⌕ Search

SEARCH · Search NASA

Results for “Chemical 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Direct time-resolved detection and quantification of key reactive intermediates in diethyl ether oxidation at T = 450–600 K

High-pressure multiplexed photoionization mass spectrometry (MPIMS) with tunable vacuum ultraviolet (VUV) ionization radiation from the Lawrence Berkeley Labs Advanced Light Source is used to investigate the oxidation of diethyl ether (DEE). Kinetics and photoionization (PI) spectra are simultaneously measured for the species formed. Several stable products from DEE oxidation are identified and quantified using reference PI cross-sections. In addition, we directly detect and quantify three key chemical intermediates: peroxy (ROO˙), hydroperoxyalkyl peroxy (˙OOQOOH), and ketohydroperoxide (HOOP[double bond, length as m-dash]O, KHP). These intermediates undergo dissociative ionization (DI) into smaller fragments, making their identification by mass spectrometry challenging. Furthermore, with the aid of quantum chemical calculations, we identify the DI channels of these key chemical species and quantify their time-resolved concentrations from the overall carbon atom balance at T = 450 K and P = 7500 torr. This allows the determination of the absolute PI cross-sections of ROO˙, ˙OOQOOH, and KHP into each DI channel directly from experiment. The PI cross-sections in turn enable the quantification of ROO˙, ˙OOQOOH, and KHP from DEE oxidation over a range of experimental conditions that reveal the effects of pressure, O 2 concentration, and temperature on the competition among radical decomposition and second O 2 addition pathways.

42 ENGINEERING↗
DOE OSTI2024

State Preparation in Quantum Algorithms for Fragment-Based Quantum Chemistry

State preparation for quantum algorithms is crucial for achieving high accuracy in quantum chemistry and competing with classical algorithms. The localized active space–unitary coupled cluster (LAS–UCC) algorithm iteratively loads a fragment-based multireference wave function onto a quantum computer. Here, in this study, we compare two state preparation methods, quantum phase estimation (QPE) and direct initialization (DI), for each fragment. We test the two state preparation methods on three systems, ranging from a model system, a set of interacting hydrogen molecules, to more realistic chemical problems, like the C–C double bond breaking in transbutadiene and the spin ladder in a bimetallic system. We analyze the impact of QPE parameters, such as the number of ancilla qubits and Trotter steps, on the prepared state. We find a trade-off between the methods, where DI requires fewer resources for smaller fragments, while QPE is more efficient for larger fragments. Our resource estimates highlight the benefits of system fragmentation in state preparation for subsequent quantum chemical calculations. These findings have broad applications for preparing multireference quantum chemical wave functions on quantum circuits that can be used for realistic chemical applications.

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

Thioester synthesis by a designed nickel enzyme models prebiotic energy conversion

The formation of carbon–carbon bonds from prebiotic precursors such as carbon dioxide represents the foundation of all primordial life processes. In extant organisms, this reaction is carried out by the carbon monoxide dehydrogenase (CODH)/acetyl coenzyme A synthase (ACS) enzyme, which performs the cornerstone reaction in the ancient Wood–Ljungdahl metabolic pathway to synthesize the key biological metabolite, acetyl-CoA. Despite its significance, a fundamental understanding of this transformation is lacking, hampering efforts to harness analogous chemistry. To address these knowledge gaps, we have designed an artificial metalloenzyme within the azurin protein scaffold as a structural, functional, and mechanistic model of ACS. We demonstrate the intermediacy of the Ni I species and requirement for ordered substrate binding in the bioorganometallic carbon–carbon bond-forming reaction from the one-carbon ACS substrates. The electronic and geometric structures of the nickel-acetyl intermediate have been characterized using time-resolved optical, electron paramagnetic resonance, and X-ray absorption spectroscopy in conjunction with quantum chemical calculations. Moreover, we demonstrate that the nickel-acetyl species is chemically competent for selective acyl transfer upon thiol addition to biosynthesize an activated thioester. Drawing an analogy to the native enzyme, a mechanism for thioester generation by this ACS model has been proposed. The fundamental insight into the enzymatic process provided by this rudimentary ACS model has implications for the evolution of primitive ACS-like proteins. Ultimately, these findings offer strategies for development of highly active catalysts for sustainable generation of liquid fuels from one-carbon substrates, with potential for broad applications across diverse fields ranging from energy storage to environmental remediation.

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

A rotationally invariant approach based on Gutzwiller wave function for correlated electron systems

Here, we introduce a rotationally invariant approach combined with the Gutzwiller conjugate gradient minimization method to study correlated electron systems. In the approach, the Gutzwiller projector is parametrized based on the number of electrons occupying the onsite orbitals instead of the onsite configurations. The approach efficiently groups the onsite orbitals according to their symmetry and greatly reduces the computational complexity, which yields a speedup of $20 \sim 50 \times $ in the minimal basis energy calculation of dimers. The computationally efficient approach promotes more accurate calculations beyond the minimal basis that is inapplicable in the original approach. A large-basis energy calculation of F 2 demonstrates favorable agreements with standard quantum-chemical calculations Bytautas et al (2007 J. Chem. Phys. 127 164317).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Quantum Calculations of VX Ammonolysis and Hydrolysis Pathways via Hydrated Lithium Nitride

Recently, lithium nitride (Li 3 N) has been proposed as a chemical warfare agent (CWA) neutralization reagent for its ability to produce nucleophilic ammonia molecules and hydroxide ions in aqueous solution. Quantum chemical calculations can provide insight into the Li 3 N neutralization process that has been studied experimentally. Here, we calculate reaction-free energies associated with the Li3N-based neutralization of the CWA VX using quantum chemical density functional theory and ab initio methods. We find that alkaline hydrolysis is more favorable to either ammonolysis or neutral hydrolysis for initial P-S and P-O bond cleavages. Reaction-free energies of subsequent reactions are calculated to determine the full reaction pathway. Notably, products predicted from favorable reactions have been identified in previous experiments.

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

Long-range exciton diffusion in molecular non-fullerene acceptors

Abstract The short exciton diffusion length associated with most classical organic semiconductors used in organic photovoltaics (5-20 nm) imposes severe limits on the maximum size of the donor and acceptor domains within the photoactive layer of the cell. Identifying materials that are able to transport excitons over longer distances can help advancing our understanding and lead to solar cells with higher efficiency. Here, we measure the exciton diffusion length in a wide range of nonfullerene acceptor molecules using two different experimental techniques based on photocurrent and ultrafast spectroscopy measurements. The acceptors exhibit balanced ambipolar charge transport and surprisingly long exciton diffusion lengths in the range of 20 to 47 nm. With the aid of quantum-chemical calculations, we are able to rationalize the exciton dynamics and draw basic chemical design rules, particularly on the importance of the end-group substituent on the crystal packing of nonfullerene acceptors.

14 SOLAR ENERGY↗
DOE OSTI2023

Uncovering novel liquid organic hydrogen carriers: a systematic exploration of chemical compound space using cheminformatics and quantum chemical methods

We present a comprehensive, in silico-based discovery approach to identifying novel liquid organic hydrogen carrier (LOHC) candidates using cheminformatics methods and quantum chemical calculations. We screened over 160 billion molecules from ZINC15 and GDB-17 chemical databases for structural similarity to known LOHCs and employed a data-driven selection criterion connecting molecular features with dehydrogenation enthalpy. This scoring criterion effectively predicts dehydrogenation enthalpies from SMILES strings, streamlining the LOHC screening process. After rigorous screening and down-selection, we compiled a database of 3000 dehydrogenation reactions for the most promising LOHC candidates, setting the stage for future selection based on kinetics and catalysis. This work demonstrates the significant impact of integrating quantum chemistry and cheminformatics in materials discovery, accelerating the selection process while reducing experimental efforts and time. By proposing new molecules as prospective LOHC candidates, our study provides a valuable resource for researchers and engineers in the development of advanced LOHC systems and showcases a successful approach for high-throughput discovery, contributing to more efficient and sustainable energy storage solutions.

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

Ground and excited state properties of ThBe and AcBe

In this work, the ground and excited states of ThBe and AcBe were investigated by performing high-level multireference and single-reference coupled-cluster quantum chemical calculations with large correlation consistent basis sets. Full potential energy curves (PECs), chemical bonding patterns, energetics, spectroscopic parameters (T e , r e , ω e , and ω e x e ), and spin–orbit effects of 13 and 8 electronic states of ThBe and AcBe, respectively, are reported. The ground electronic states of ThBe and AcBe are single-reference 1 3 Σ − (1σ 2 2σ 2 1π 2 ) and 1 2 Π (1σ 2 2σ 2 1π 1 ), respectively, and originate from their corresponding ground state fragments. The chemical bonding of ThBe (1 3 Σ − ) and AcBe (1 2 Π) are π-dative in character and are formed by d-electron transfers from Th/Ac to the empty 2p x and 2p y of the Be atom. The electron populations of the f-orbitals of both ThBe (1 3 Σ − ) and AcBe (1 2 Π) are minor which exhibit their “transition-metal-like” nature. The estimated bond energies of the spin–orbit ground states of ThBe (1 3 Σ − 0+ ) and AcBe (1 2 Π 1/2 ) are 12.79 and 11.02 kcal mol −1 , respectively. Finally, the bond energy of ThBe was used to estimate its heat of formation ΔH 0 f (298 K) of 869.61 ± 6 kJ mol −1 .

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Determining Partial Atomic Charges for Liquid Water: Assessing Electronic Structure and Charge Models

Partial atomic charges provide an intuitive and efficient way to describe the charge distribution and the resulting intermolecular electrostatic interactions in liquid water. Many charge models exist and it is unclear which model provides the best assignment of partial atomic charges in response to the local molecular environment. In this work, we systematically scrutinize various electronic structure methods and charge models (Mulliken, natural population analysis, CHelpG, RESP, Hirshfeld, Iterative Hirshfeld, and Bader) by evaluating their performance in predicting the dipole moments of isolated water, water clusters, and liquid water as well as charge transfer in the water dimer and liquid water. Although none of the seven charge models is capable of fully capturing the dipole moment increase from isolated water (1.85 D) to liquid water (about 2.9 D), the Iterative Hirshfeld method performs best for liquid water, reproducing its experimental average molecular dipole moment, yielding a reasonable amount of intermolecular charge transfer, and showing modest sensitivity to the local water environment. The performance of the charge model is dependent on the choice of the density functional and the quantum treatment of the environment. The computed molecular dipole moment of water generally increases with the percentage of the exact Hartree–Fock exchange in the functional, whereas the amount of charge transfer between molecules decreases. For liquid water, including two full solvation shells of surrounding water molecules (within about 5.5 Å of the central water) in the quantum chemical calculation converges the charges of the central water molecule. Furthermore, our final pragmatic quantum chemical charge-assigning protocol for liquid water is the Iterative Hirshfeld method with M06-HF/aug-cc-pVDZ and a quantum region cutoff radius of 5.5 Å.

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

Potassium iodide cluster based superhalogens and superalkalis: Theoretical calculations and experimental confirmation

We studied a series of potassium-iodide clusters with formulas of (KI)nI- and (KI)nK+ (n = 1-3) by quantum chemical calculations. The calculated vertical detachment energies (VDEs) of (KI)nI- clusters are all higher than those of halogen anions, and thus can be classified as superhalogen anions; while the calculated vertical electron affinities (VEAs) of (KI)nK+ clusters are all lower than those of alkali metal cations, and thus they can be recognized as superalkali cations. Our calculated VEA for (KI)K+, i.e., K2I+, agrees well with the previously measured value by thermal ionization mass spectrometry. We confirmed the theoretical predictions by measuring the negative ion photoelectron (NIPE) spectra of (KI)nI- anions, from which the measured VDEs are found to be in good agreement with the theoretically calculated values. In addition, we also explored the possibility of supersalts formed by these superhalogen anions and superalkali cations, for example, K3I3 can be regarded as a supersalt formed by K2I+ and KI2– subunits.

potassium-iodide clusters, superhalogen, superalka↗
DOE OSTI2021

Crossover from hydrogen to chemical bonding

Hydrogen bonds (H-bonds) can be interpreted as a classical electrostatic interaction or as a covalent chemical bond if the interaction is strong enough. As a result, short strong H-bonds exist at an intersection between qualitatively different bonding descriptions, with few experimental methods to understand this dichotomy. The [F-H-F] − ion represents a bare short H-bond, whose distinctive vibrational potential in water is revealed with femtosecond two-dimensional infrared spectroscopy. It shows the superharmonic behavior of the proton motion, which is strongly coupled to the donor-acceptor stretching and disappears on H-bond bending. In combination with high-level quantum-chemical calculations, we demonstrate a distinct crossover in spectroscopic properties from conventional to short strong H-bonds, which identify where hydrogen bonding ends and chemical bonding begins.

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

Ab Initio Direct Dynamics

The reactivity and dynamics of molecular systems can be explored computationally by classical trajectory calculations. The traditional approach involves fitting a functional form of a potential energy surface (PES) to the energies from a large number of electronic structure calculations and then integrating numerous trajectories on this fitted PES to model the molecular dynamics. The ever-decreasing cost of computing and continuing advances in computational chemistry software have made it possible to use electronic structure calculations directly in molecular dynamics simulations without first having to construct a fitted PES. In this “on-the-fly” approach, every time the energy and its derivatives are needed for the integration of the equations of motion, they are obtained directly from quantum chemical calculations. This approach started to become practical in the mid-1990s as a result of increased availability of inexpensive computer resources and improved computational chemistry software. The application of direct dynamics calculations has grown rapidly over the last 25 years and would require a lengthy review article. The present Account is limited to some of our contributions to methods development and various applications. To improve the efficiency of direct dynamics calculations, we developed a Hessian-based predictor-corrector algorithm for integrating classical trajectories. Hessian updating made this even more efficient. Furthermore, this approach was also used to improve algorithms for following the steepest descent reaction paths. For larger molecular systems, we developed an extended Lagrangian approach in which the electronic structure is propagated along with the molecular structure. Strong field chemistry is a rapidly growing area, and to improve the accuracy of molecular dynamics in intense laser fields, we included the time-varying electric field in a novel predictor-corrector trajectory integration algorithm. Since intense laser fields can excite and ionize molecules, we extended our studies to include electron dynamics. Specifically, we developed code for time-dependent configuration interaction electron dynamics to simulate strong field ionization by intense laser pulses. Our initial application of ab initio direct dynamics in 1994 was to CH 2 O → H 2 + CO; the calculated vibrational distributions in the products were in very good agreement with experiment. In the intervening years, we have used direct dynamics to explore energy partitioning in various dissociation reactions, unimolecular dissociations yielding three fragments, reactions with branching after the transition state, nonstatistical dynamics of chemically activated molecules, dynamics of molecular fragmentation by intense infrared laser pulses, selective activation of specific dissociation channels by aligned intense infrared laser fields, angular dependence of strong field ionization, and simulation of sequential double ionization.

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

Understanding the chemical bonding of ground and excited states of HfO and HfB with correlated wavefunction theory and density functional approximations

Knowledge of the chemical bonding of HfO and HfB ground and low-lying electronic states provides essential insights into a range of catalysts and materials that contain Hf–O or Hf–B moieties. Here, we carry out high-level multi-reference configuration interaction theory and coupled cluster quantum chemical calculations on these systems. We compute full potential energy curves, excitation energies, ionization energies, electronic configurations, and spectroscopic parameters with large quadruple-ζ and quintuple-ζ quality correlation consistent basis sets. We also investigate equilibrium chemical bonding patterns and effects of correlating core electrons on property predictions. Differences in the ground state electron configuration of HfB(X 4 Σ - ) and HfO(X 1 Σ + ) lead to a significantly stronger bond in HfO than HfB, as judged by both dissociation energies and equilibrium bond distances. We extend our analysis to the chemical bonding patterns of the isovalent HfX (X = O, S, Se, Te, and Po) series and observe similar trends. We also note a linear trend between the decreasing value of the dissociation energy (D e ) from HfO to HfPo and the singlet–triplet energy gap (ΔE S–T ) of the molecule. Finally, we compare these benchmark results to those obtained using density functional theory (DFT) with 23 exchange–correlation functionals spanning multiple rungs of “Jacob’s ladder.” When comparing DFT errors to coupled cluster reference values on dissociation energies, excitation energies, and ionization energies of HfB and HfO, we observe semi-local generalized gradient approximations to significantly outperform more complex and high-cost functionals.

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

On the Finite Nuclear Effect and Gaussian Basis Sets for Four-Component Dirac Hartree-Fock Calculations

In this perspective we describe an overview of Dirac Hartree-Fock theory, which is the foundation of relativistic quantum chemistry methods vital in the molecular modeling of heavy element chemistry and physics. A special emphasis was placed upon the physical connection between the finite nuclear model and the kinetic-balance condition between the large and small components. Additionally, we review the consequences of basis set choices for Dirac Hartree-Fock calculations in the restricted kinetic balance formalism, and show that the basis set limit does not have a hard lower bound due to the nature of the Hamiltonian, referred to in the literature as prolapse. Standard Gaussian basis set contraction schemes can break down in Dirac Hartree-Fock calculations, and special care must be taken for the description of core orbitals near the nucleus. The phenomena described in this perspective are vital for theoreticians and practitioners to understand when selecting a basis set for heavy element quantum chemical calculations using Dirac Hartree-Fock theory.

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

Electronic Structure Calculations of the Fundamental Interactions in Solvent Extraction Desalination

Solvent extraction desalination (SED) is a promising membrane-free desalination technology, but new solvents are still needed that can improve the molecular-level water selectivity, while minimizing solvent contamination in the aqueous phase. Improvements in solvent design for SED require a fundamental understanding of the molecular interactions governing solvation behavior. Herein, quantum chemical calculations are used to delineate different fundamental interactions in water-solvent systems involving 38 different solvents, with a particular focus on the role of the electrostatic potential characteristics. Furthermore, these interactions are rigorously analyzed using the independent gradient model based on Hirshfeld partition (IGMH), Symmetry-Adapted Perturbation Theory (SAPT), as well as the energy decomposition analysis based on force field (EDA-FF) to compare against a standard force field model. The solvation free energies and partition coefficients are also calculated to estimate water extraction performance in a bulk solvent system.

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

Advancing molecular machine learning representations with stereoelectronics-infused molecular graphs

Molecular representation is a critical element in our understanding of the physical world and the foundation for modern molecular machine learning. Previous molecular machine learning models have used strings, fingerprints, global features and simple molecular graphs that are inherently information-sparse representations. However, as the complexity of prediction tasks increases, the molecular representation needs to encode higher fidelity information. This work introduces a new approach to infusing quantum-chemical-rich information into molecular graphs via stereoelectronic effects, enhancing expressivity and interpretability. Learning to predict the stereoelectronics-infused representation with a tailored double graph neural network workflow enables its application to any downstream molecular machine learning task without expensive quantum-chemical calculations. We show that the explicit addition of stereoelectronic information substantially improves the performance of message-passing two-dimensional machine learning models for molecular property prediction. We show that the learned representations trained on small molecules can accurately extrapolate to much larger molecular structures, yielding chemical insight into orbital interactions for previously intractable systems, such as entire proteins, opening new avenues of molecular design. Finally, we have developed a web application (simg.cheme.cmu.edu) where users can rapidly explore stereoelectronic information for their own molecular systems.

Boiko, Daniil A↗
DOE OSTI2024

Inhibition of asphaltene aggregation using deep eutectic solvents: COSMO-RS calculations and experimental validation

Asphaltene precipitation/deposition adversely affects various oil and gas processes, such as oil recovery, transportation, and petroleum processing. The resulting increase in viscosity of crude oil decreases distillate yields creating unstable phase separation. Deep eutectic solvents (DESs) have recently gained interest as inhibitors of asphaltene aggregation. The goal of the present study is to derive the mechanism of inhibition of petroleum asphaltene aggregation by a novel screening of DESs. An archipelago-based chemical structure of asphaltene was adopted for performing quantum chemical calculations. The structure was used in the conductor-like screening model for real solvents (COSMO-RS) model to screen potential DESs for asphaltene precipitation. It was found that DESs containing thymol were the most promising of the 153 DES combinations screened. The COSMO-RS predictions were validated experimentally using solubility data of asphaltene in DESs. Among the studied DESs, thymol-diphenyl ether provided the highest solubility for asphaltene, which was further validated using the experimental and the COSMO-RS predicted data. In addition, to characterize the structural interactions between asphaltene and DESs, Fourier transform infrared (FTIR) spectroscopy and nuclear magnetic resonance (NMR) measurements were performed. It was found that strong interaction between asphaltene and the DESs is responsible for the higher asphaltene dispersion. Finally, the present approach opens pathways to rationally design and understand the impact of structural variation of DESs based on their interactions with asphaltene.

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

Activation of Semicore Electrons in Alkali Metals and Their Role in the B1–B2 Phase Transition under Pressure

Alkali fluorides are often thought of as archetypical ionic compounds whose structures can be understood in terms of the packing of rigid spheres. At ambient pressure, they assume the rocksalt (B1) structure, while under only a few GPa of pressure, the cesium chloride (B2) structure, with higher coordination numbers, is assumed for KF, RbF, and CsF. NaF requires almost 10 times more pressure to undergo this same phase transition, which has not been observed for LiF. Herein, we provide a detailed analysis, based upon quantum chemical calculations, explaining this behavior. We show that for the heavier alkali metals, the semicore p orbitals engage in metal-metal bonding in the B2 phase, facilitating the pressure-induced B1 → B2 structural transition. Furthermore, these findings suggest that the semicore orbitals of heavy alkali metals can be activated without the need of strong oxidants at very mild levels of compression, resulting in the formation of chemical bonds, challenging both traditional and modern core-valence distinctions. In addition, we argue that Cs 5p-5p bonding in CsCl occurs already at ambient pressure, stabilizing the B2 phase, and suggest experiments that may be able to detect signatures of such bonding.

Anions↗