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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS1994

New calculations on the ion-molecule processes C2H2(+) + H2 --> C2H3(+) + H and C2H2(+) + H2 --> C2H4+

New high-level quantum chemical calculations have been undertaken to understand the rates and mechanisms of the reactive and associative channels for the reactants C2H2(+) + H2. The reactive channel, which produces C2H3(+) + H, has been shown to be slightly endothermic, confirming earlier calculations at a somewhat lower level and in agreement with some recent experimental work. The associative channel, leading to C2H4+, has been shown to proceed via a transition state with negative energy relative to the reactants, so that association is predicted to be efficient. This result is in conflict with an earlier theoretical study but in agreement with low-temperature experimental measurements.

NASA Discipline Exobiology↗
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↗
NASA NTRS2022

Design of Hypervelocity Flow Generator (HFG) and Its Diagnostics

Ground facilities for hypersonic research are a key element for successful development of aerodynamically proven hypervelocity vehicles. Design concepts and diagnostics of a hypervelocity flow generator (HFG) were made as a test platform for hypersonic/hypervelocity spacecraft models at the NASA Langley Research Center. The HFG is a hypersonic flow field generator using optically heated gas which is blown into an 80 cu m vacuum chamber. The vacuum chamber is kept at a stable vacuum pressure with a combination of three large vacuum pumps, while the HFG is in the test mode. The HFG provides a relatively small test section with approximately a 20 cm window. This facility was designed to generate 2.45 km/sec of flow speed, and potentially generate a continuous flow with the nozzle and vacuum system. The window in test section provides a direct view of the shock wave around a model in order to measure temperature, pressure, and density profile within the shock layer. One of the key test goals under this project is to understand why the emission spectra from the standing shockwave plasma predicted by Lora-Loran codes are significantly different from the measured emission spectra from the Flight Investigation Reentry Environment (FIRE II) Flight. The correct estimate of the thermal loading on the leading edge of hypersonic vehicles greatly affects the aerodynamic design, the material selection for the vehicle, and the cooling requirement and can be obtained by the precise modeling of emission spectra from shockwave plasma. However, the estimation of thermal loading is not an easy task due to complex non-equilibrium radiative process within high temperature shock layers that still falls into a category of cold plasma. Direct flight experiments are the most desirable, but not a cost-effective approach. Analysis by computational fluid dynamics (CFD) offers many test flexibilities. However, the CFD codes must be fully tested and validated with experimental data before the codes are effectively used for practical design. Large discrepancies between experiments and codes appear in hypersonic/hypervelocity flow regimes at high altitudes of 60 km ~ 90 km. This HFG facility offers some important parameters for CFD code validation, such as collision cross-sections, relaxation times, reaction rate coefficients and transportation coefficients. The HFG test facility is based on the ejection flow of high temperature gas heated over to 3500 K through a nozzle. The tungsten gas chamber of the HFG is heated up to a desired temperature by a 60 kW optical power beam source. This system consists of an optical power source, a thermal chamber, an expansion nozzle, a test section, and an 80 cu m vacuum tank. 60 kW optical input power is obtained from the 150 kW Vortek arc lamp system (by Vortek Industries, Vancouver, Canada). This optical beam is focused to heat the gas chamber within which a flow media is heated. The maximum achievable temperature of the flow medium reached approximately 3500 K or even higher but is limited by the melting point of the chamber material used. The exhaust velocity through the nozzle was determined by the stagnation temperature and the molecular weight of the working medium at the test section. To provide design parameters, a NASA Chemical Equilibrium with Application (CEA) computer program is used for the simulation of aerothermal data. This CEA program can calculate chemical equilibrium and properties of complex mixtures using shock tube parameters. For nitrogen gas at 2666 K stagnation temperature, the maximum achievable velocity at the test section is approximately 2.45 km/sec which is within the range of the thermal velocity of 8000 K shock layer. Based on the calculation through the CEA program, the design parameters of a HFG were determined and implemented for the test section that includes an expansion nozzle. The installed test section of the experimental facility can sustain a condition of a re-entry vehicle from the Space at an altitude of 60-90 km. After installation of the HFG, the system was fully tested and its operational parameters were measured. An 80 cu m vacuum chamber of HFG was set at 1 torr level to keep a stable downstream condition. The pumping time to reach the minimum vacuum pressure (~ 1 torr) at the test chamber from atmospheric pressure was approximately an hour using Kinney (MBV-14000/MB -1600/KT-300) and Stokes (Model 1772 and 412) pumps. Such a setup condition allowed a continued stable operation of the HFG experiment with flowrates through 1-, 3-, or 5-mm diameters of nozzle throat. The flow characteristics of the HFG for various operating conditions were performed using a focused Schlieren method. At 0.8 torr chamber pressure, a barrel shock was observed at the test section. With a sphere obstacle of 12 mm diameter in the flow, a bow shock (~ 2-mm thickness) was observed by a focused Schlieren visualization method.

Hypervelocity flow↗
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↗
NASA NTRS1992

Effect of inert propellant injection on Mars ascent vehicle performance

A Mars ascent vehicle is limited in performance by the propellant which can be brought from Earth. In some cases the vehicle performance can be improved by injecting inert gas into the engine, if the inert gas is available as an in-situ resource and does not have to be brought from Earth. Carbon dioxide, nitrogen, and argon are constituents of the Martian atmosphere which could be separated by compressing the atmosphere, without any chemical processing step. The effect of inert gas injection on rocket engine performance was analyzed with a numerical combustion code that calculated chemical equilibrium for engines of varying combustion chamber pressure, expansion ratio, oxidizer/fuel ratio, and inert injection fraction. Results of this analysis were applied to several candidate missions to determine how the required mass of return propellant needed in low Earth orbit could be decreased using inert propellant injection.

Colvin, James E.↗
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↗
NASA NTRS1985

A finite volume method for the calculation of compressible chemically reacting flows

Several efficient pseudo time techniques have been developed for calculating steady state chemically reacting flows. The techniques include the implicit treatment of the chemical source term, point implicit multiple grid accelerator and a constant CFL condition. It turns out that these methods can be viewed as ways of rescaling the equations in time such that all chemical and convective phenomena evolve at comparable pseudo time scales. Consequently the number of iterations needed to solve reacting problems is approximately the same as for non-reacting problems. The techniques are demonstrated for a simple dissociation model and a nontrivial H2 - Air combustion model.

Bussing, T. R. A.↗
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↗
NASA NTRS1994

Computer program for calculation of complex chemical equilibrium compositions and applications. Part 1: Analysis

This report presents the latest in a number of versions of chemical equilibrium and applications programs developed at the NASA Lewis Research Center over more than 40 years. These programs have changed over the years to include additional features and improved calculation techniques and to take advantage of constantly improving computer capabilities. The minimization-of-free-energy approach to chemical equilibrium calculations has been used in all versions of the program since 1967. The two principal purposes of this report are presented in two parts. The first purpose, which is accomplished here in part 1, is to present in detail a number of topics of general interest in complex equilibrium calculations. These topics include mathematical analyses and techniques for obtaining chemical equilibrium; formulas for obtaining thermodynamic and transport mixture properties and thermodynamic derivatives; criteria for inclusion of condensed phases; calculations at a triple point; inclusion of ionized species; and various applications, such as constant-pressure or constant-volume combustion, rocket performance based on either a finite- or infinite-chamber-area model, shock wave calculations, and Chapman-Jouguet detonations. The second purpose of this report, to facilitate the use of the computer code, is accomplished in part 2, entitled 'Users Manual and Program Description'. Various aspects of the computer code are discussed, and a number of examples are given to illustrate its versatility.

Gordon, Sanford↗