Calculation of complex chemical equilibria.
Complex chemical equilibrium calculation, analyzing thermodynamic principles and solution to nonlinear equations
Gordon, S.↗
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.
US EPA Superfund Remediation Program Update: Revisions to Risk and Dose Assessment Models - 20368
The U.S. Environmental Protection Agency (EPA) Superfund remedial program is finishing up a significant number of revisions to its guidance for the risk assessment process at radioactively contaminated Superfund sites. The six Preliminary Remediation Goal (PRG) and Dose Compliance Concentration (DCC) internet-based calculators for radiation risk and dose assessment at Superfund sites are being reformatted to match the chemical calculators and improve the user experience, revised to add a new peak risk and dose output option, and to add some adjustment options for external exposure. A new Radon Vapor Intrusion Screening Level (RVISL) calculator is expected to be finished in 2020, will be similar to the Vapor Intrusion Screening Level (VISL) for chemical calculator that was finished in 2018. (authors)
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
Uncertainty propagation in pore water chemical composition calculation using surrogate models
Performance assessment in deep geological nuclear waste repository systems necessitates an extended knowledge of the pore water chemical conditions prevailing in host-rock formations. In the last two decades, important progress has been made in the experimental characterization and thermodynamic modeling of pore water speciation, but the influence of experimental artifacts and uncertainties of thermodynamic input parameters are seldom evaluated. In this respect, we conducted an uncertainty propagation study in a reference geochemical model describing the pore water chemistry of the Callovian-Oxfordian clay formation. Nineteen model input parameters were perturbed, including those associated to experimental characterization (leached anions, exchanged cations, cation exchange selectivity coefficients) and those associated to generic thermodynamic databases (solubilities). A set of 13 quantities of interest were studied by the use of polynomial chaos expansions built non-intrusively with a least-squares forward stepwise regression approach. Training and validation sets of simulations were carried out using the geochemical speciation code PHREEQC. The statistical results explored the marginal distribution of each quantity of interest, their bivariate correlations as well as their global sensitivity indices. The influence of the assumed distributions for input parameters uncertainties was evaluated by considering two parametric domain sizes.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Adaptive basis sets for practical quantum computing
Electronic structure calculations on small systems such as H 2 , H 2 O, LiH, and BeH 2 with chemical accuracy are still a challenge for the current generation of noisy intermediate-scale quantum (NISQ) devices. One of the reasons is that due to the device limitations, only minimal basis sets are commonly applied in quantum chemical calculations, which allows one to keep the number of qubits employed in the calculations at a minimum. However, the use of minimal basis sets leads to very large errors in the computed molecular energies as well as potential energy surface shapes. One way to increase the accuracy of electronic structure calculations is through the development of small basis sets better suited for quantum computing. In this work, we show that the use of adaptive basis sets, in which exponents and contraction coefficients depend on molecular structure, provides an easy way to dramatically improve the accuracy of quantum chemical calculations without the need to increase the basis set size and thus the number of qubits utilized in quantum circuits. As a proof of principle, we optimize an adaptive minimal basis set for quantum computing calculations on an H 2 molecule, in which exponents and contraction coefficients depend on the H—H distance, and apply it to the generation of H 2 potential energy surface on IBM-Q quantum devices. The adaptive minimal basis set reaches the accuracy of the double-zeta basis sets, thus allowing one to perform double-zeta quality calculations on quantum devices without the need to utilize twice as many qubits in simulations. This approach can be extended to other molecular systems and larger basis sets in a straightforward manner.
97 MATHEMATICS AND COMPUTING↗
Horizontal mixing coefficients for two-dimensional chemical models calculated from National Meteorological Center Data
Calculations of the two-dimensional, species-independent mixing coefficients for two-dimensional chemical models for the troposphere and stratosphere are performed using quasi-geostrophic potential vorticity fluxes and gradients from 4 years of National Meteorological Center data for the four seasons in both hemispheres. Results show that the horizontal mixing coefficient values for the winter lower stratosphere are broadly consistent with those currently employed in two-dimensional models, but the horizontal mixing coefficient values in the northern winter upper stratosphere are much larger than those usually used.
Newman, P. A.↗
Metastable isomers - A new class of interstellar molecules
The abundances of a variety of metastable isomers of small organic molecules, analogous to HNC/HCN, in dense interstellar clouds are considered. These metastable species, some of which are thought to exist as intermediates in laboratory organic chemical reactions, are of considerable interest to chemists. Current ideas of gas-phase, ion-molecule chemistry are utilized to demonstrate that such metastable species should often be present in dense clouds in sufficient abundance to be observed. Unfortunately, the spectral constants of metastable isomers have rarely been determined in the laboratory, and quantum chemical calculations of a varying degree of accuracy must be utilized; results are included of some new quantum chemical calculations. The interstellar chemistry and expected microwave spectra of a representative sample of possibly important interstellar metastable isomers are discussed.
Green, S.↗
Applicability of PM3 to transphosphorylation reaction path: Toward designing a minimal ribozyme
A growing body of evidence shows that RNA can catalyze many of the reactions necessary both for replication of genetic material and the possible transition into the modern protein-based world. However, contemporary ribozymes are too large to have self-assembled from a prebiotic oligonucleotide pool. Still, it is likely that the major features of the earliest ribozymes have been preserved as molecular fossils in the catalytic RNA of today. Therefore, the search for a minimal ribozyme has been aimed at finding the necessary structural features of a modern ribozyme (Beaudry and Joyce, 1990). Both a three-dimensional model and quantum chemical calculations are required to quantitatively determine the effects of structural features of the ribozyme on the reaction it catalyzes. Using this model, quantum chemical calculations must be performed to determine quantitatively the effects of structural features on catalysis. Previous studies of the reaction path have been conducted at the ab initio level, but these methods are limited to small models due to enormous computational requirements. Semiempirical methods have been applied to large systems in the past; however, the accuracy of these methods depends largely on a simple model of the ribozyme-catalyzed reaction, or hydrolysis of phosphoric acid. We find that the results are qualitatively similar to ab initio results using large basis sets. Therefore, PM3 is suitable for studying the reaction path of the ribozyme-catalyzed reaction.
Manchester, John I.↗
Quantum-Chemically Informed Machine Learning: Prediction of Energies of Organic Molecules with 10 to 14 Non-hydrogen Atoms
High-fidelity quantum-chemical calculations can provide accurate predictions of molecular energies, but their high computational costs limit their utility, especially for larger molecules. We have shown in previous work that machine learning models trained on high-level quantum-chemical calculations (G4MP2) for organic molecules with one to nine non-hydrogen atoms can provide accurate predictions for other molecules of comparable size at much lower costs. Here we demonstrate that such models can also be used to effectively predict energies of molecules larger than those in the training set. To implement this strategy, we first established a set of 191 molecules with 10–14 non-hydrogen atoms having reliable experimental enthalpies of formation. We then assessed the accuracy of computed G4MP2 enthalpies of formation for these 191 molecules. The error in the G4MP2 results was somewhat larger than that for smaller molecules, and the reason for this increase is discussed. Two density functional methods, B3LYP and ωB97X-D, were also used on this set of molecules, with ωB97X-D found to perform better than B3LYP at predicting energies. The G4MP2 energies for the 191 molecules were then predicted using these two functionals with two machine learning methods, the FCHL-Δ and SchNet-Δ models, with the learning done on calculated energies of the one to nine non-hydrogen atom molecules. The better-performing model, FCHL-Δ, gave atomization energies of the 191 organic molecules with 10–14 non-hydrogen atoms within 0.4 kcal/mol of their G4MP2 energies. Thus, this work demonstrates that quantum-chemically informed machine learning can be used to successfully predict the energies of large organic molecules whose size is beyond that in the training set.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Quantum Monte Carlo Calculations of Chemical Binding and Reactions
The auxiliary field quantum Monte Carlo method developed by the PIs has been shown to provide the most accurate description of strongly correlated electronic systems, from molecules to solids. Unlike other explicitly many‐body approaches, the quantum Monte Carlo method scales as a low order polynomial of systems size, similar to mean‐field methods such as density functional theory. However, the auxiliary field quantum Monte Carlo algorithm is significantly more expensive than traditional density functional calculations. This creates a bottleneck for applications to extended systems, such as large molecules and solids. One principal objective of this proposal was to develop new auxiliary field quantum Monte Carlo computational strategies to achieve improved scaling with system size, using downfolding and localization schemes, without sacrificing the predictive power of the calculations. A second goal is to extend the reach of auxiliary field quantum Monte Carlo to calculate excited states. This final report summarizes what has been achieved during the course the project toward these goals.
97 MATHEMATICS AND COMPUTING↗
The efficient calculation of chemically reacting flow
A semi-implicit finite volume formulation is used to study flows with chemical reactions. In this formulation the source terms resulting from the chemical reactions are treated implicitly and the resulting system of partial differential equations is solved using two time-stepping schemes. The first is based on the Runge-Kutta method while the second is based on an Adams predictor-corrector method. Results show that improvements in computational efficiency depend to a large extent on the manner in which the source term is treated. Further, analysis and computation indicate that the Runge-Kutta method is more efficient than the Adams methods. Finally, an adaptive time stepping scheme is developed to study problems involving shock ignition. Calculations for a hydrogen-air system agree well with other methods.
Eklund, D. R.↗
Activation of methane by U + studied by guided ion beam tandem mass spectrometry and quantum chemistry
Reaction pathways of all products formed in the U + + CH 4 (CD 4 ) reaction were explored as a function of kinetic energy using guided ion beam tandem mass spectrometry and quantum chemical calculations. UH + , UC + , UCH + , UCH 2 + , and UCH 3 + (and their perdeuterated analogues) are formed in endothermic reactions. In both systems, the UCH 2 + (UCD 2 + ) dehydrogenated product was the dominant product in the low-energy region, whereas the UH + (UD + ) hydride product became predominant at high energies. The kinetic energy behavior of the various products is consistent with a common intermediate of H–U + –CH 3 (D–U + –CD 3 ). Here, the kinetic energy dependence of all product cross sections was modeled to obtain experimental bond dissociation energies at 0 K (in eV): D 0 (U + –H) = 2.42 ± 0.10, D 0 (U + –C) = 3.95 ± 0.12, D 0 (U + –CH) = 4.91 ± 0.09, D 0 (U + –CH 2 ) = 4.11 ± 0.04, and D 0 (U + –CH 3 ) = 2.41 ± 0.09. Quantum chemical calculations using the UCCSD(T) and UB3LYP approaches with the cc-pwCVXZ-PP basis set with MDF-60 pseudopotential for U + and the aug-cc-pCVXZ and aug-cc-pVXZ (X = T, Q) basis set for carbon and hydrogen, respectively, validate the experimental bond dissociation energies and outline the potential energy surface for all reactions observed. In addition, spin–orbit corrections of the bond energies for all products were calculated at a CASSCF-CASPT2-RASSI level.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Triple Emission of 5'-( para -R-Phenylene)vinylene-2-(2'-hydroxyphenyl)benzoxazole (PVHBO). Part I: Dual Emission from the Neutral Species
The effects of 5'-(para-R-phenylene)vinylene (PV) substituents on the emission properties of 2-(2'-hydroxyphenyl)benzoxazole (HBO) are analyzed using steady-state and time-resolved absorption and emission spectroscopies in addition to quantum chemical calculations. All members in the series of PVHBOs are capable of excited-state intramolecular proton transfer (ESIPT) with a solvent sensitivity that is typical of a HBO derivative to produce a normal (aka enol) emission and an excited-state tautomer (aka keto) emission. These two emission bands of the neutral dyes are discussed in the current paper. The intermolecular proton transfer, i.e., the deprotonation, of a PVHBO results in the third band of the triple emission, which is described in the succeeding paper. The placement of an electron-withdrawing substituent R on the PVHBO scaffold increases the intensity of the keto emission relative to the enol emission in hydrogen-bonding solvents. The R substituents do not significantly alter the wavelengths of the enol and keto emission bands, which are located in the blue and green regions, respectively, of the visible spectrum. The ultrafast time-resolved spectroscopies and quantum chemical calculations offer explanations on how the R group and the solvent affect the enol and keto emission properties (i.e., wavelength, lifetime, fluorescence quantum yield, and relative ratio of their emissions). Here, the key findings include the following: (1) the emission energies of both enol and keto forms are not sensitively dependent on the R substituent and (2) the solvent-engaged enol excited state is quenched more efficiently as the R substituent becomes more electron-withdrawing. A PVHBO acts as a fusion of HBO and stilbenoid that intersect at the hydroxyphenyl moiety. Depending on the solvent and other environmental conditions, PVHBOs may exhibit the ESIPT property of HBO or the substituent-dependent emission of stilbenoid. This paper and the succeeding article provide a photophysical model of PVHBOs to explain the wavelengths and relative abundances of the three emission bands (enol, keto, and anion) that these compounds are able to produce. Judicial selection of the environmental factors may drive the emission of a PVHBO into the spectral regions of blue, green, and, in a couple of cases, orange or red.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Polyrate 2023: A computer program for the calculation of chemical reaction rates for polyatomics. New version announcement
Not Available
Meana-Pañeda, Rubén↗
Towards Quantum-Chemical Level Calculations of SARS-CoV-2 Spike Protein Variants of Concern by First Principles Density Functional Theory
The spike protein (S-protein) is a crucial part of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), with its many domains responsible for binding, fusion, and host cell entry. In this review we use the density functional theory (DFT) calculations to analyze the atomic-scale interactions and investigate the consequences of mutations in S-protein domains. We specifically describe the key amino acids and functions of each domain, which are essential for structural stability as well as recognition and fusion processes with the host cell; in addition, we speculate on how mutations affect these properties. Such unprecedented large-scale ab initio calculations, with up to 5000 atoms in the system, are based on the novel concept of amino acid–amino acid-bond pair unit (AABPU) that allows for an alternative description of proteins, providing valuable information on partial charge, interatomic bonding and hydrogen bond (HB) formation. In general, our results show that the S-protein mutations for different variants foster an increased positive partial charge, alter the interatomic interactions, and disrupt the HB networks. We conclude by outlining a roadmap for future computational research of biomolecular virus-related systems.
60 APPLIED LIFE SCIENCES↗
SCF calculation of chemical shifts in polycyclic aromatic hydrocarbons
Explore the source record for details and available documents.
Long, E. R.↗
Chemical model calculations of C2, C3, CH, CN, OH, and NH2 abundances in cometary comae
Explore the source record for details and available documents.
Mitchell, G. F.↗
Efficient calculation of chemically reacting flow
Explore the source record for details and available documents.
Eklund, Dean R.↗
CET93 and CETPC: An interim updated version of the NASA Lewis computer program for calculating complex chemical equilibria with applications
The NASA Lewis chemical equilibrium program with applications continues to be improved and updated. The latest version is CET93. This code, with smaller arrays, has been compiled for use on an IBM or IBM-compatible personal computer and is called CETPC. This report is intended to be primarily a users manual for CET93 and CETPC. It does not repeat the more complete documentation of earlier reports on the equilibrium program. Most of the discussion covers input and output files, two new options (ONLY and comments), example problems, and implementation of CETPC.
Mcbride, Bonnie J.↗