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 37 records · Page 2

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↗

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↗

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↗

Multi-Scale Modeling Framework for Mercury Biogeochemistry

Multi-Scale modeling of mercury (Hg) geochemical speciation and reactions has been performed by integrating atomistic quantum chemical calculations with continuum scale speciation models. Major progress has been made in the improvement of quantum chemical models to calculate critical thermodynamic data for Hg complexes in aquatic environments. Rapid and reliable quantum chemical approaches have been developed for calculating acid dissociation constants (pK a ) and stability constants (log K), with calculated mean unsigned errors of 0.5 and 1.5 log units, respectively for ligand molecules and Hg complexes. At the continuum scale, systematic analysis of uncertainty propagation in mercury (Hg) speciation modeling has been conducted and was used to identify environmental conditions under which thermodynamic constant uncertainties are significant and recommended to be accounted for. The integrated framework for multi-scale modeling of mercury geochemistry is open to the research community through the web-based multiscale modeling aqueous speciation resource, AQUA-MER. The improved quantum chemical approaches for thermodynamic constant calculations are accessible through AQUA-MER and can be used to provide the missing constants in the continuum scale speciation calculations. In addition to low molecular mass Hg complex speciation, modeling natural aquatic environments also involve the transport of high molecular weight dissolved organic matter (DOM) in reactive flows simultaneously with equilibrium and kinetic reactions. To this end, atomistic MD simulations were performed to capture the details of aggregation, mechanisms and distribution of functional groups in DOM at the molecular level. The elemental composition and calculated bulk properties of the DOM models are in close agreement with experimental measurements. A travel-time based reactive transport model in the hyporheic zone of stream corridors was established for the multicomponent Hg-DOM-S system and implemented through PFLOTRAN.

54 ENVIRONMENTAL SCIENCES↗

Local and global approaches to treat the torsional barriers of 4-methylacetophenone using microwave spectroscopy

The Fourier transform microwave spectrum of 4-methylacetophenone recorded from 8 GHz to 18 GHz under jet-cooled conditions has revealed large tunneling splittings arising from a low barrier to internal rotation of the ring methyl group and small splittings from a high torsional barrier of the acetyl methyl group. The large splittings are especially challenging to model, while the small splittings are difficult to analyze due to the resolution limit of 120 kHz. The combination of two methyl groups undergoing internal rotations caused each rotational transition to split into five torsional species, which were resolved and fitted using a modified version of the XIAM code and the newly devel- oped ntop code to a root-mean-square deviation close to measurement accuracy, providing an estimate of the V 3 potential barriers of about 22 cm -1 and 584–588 cm -1 for the ring and the acetyl methyl groups, respectively. The assignment was aided by separately fitting the five torsional species using odd-power order operators. Only one conformer in which all heavy atoms are located on a symmetry plane could be identified in the spectrum, in agreement with results from conformation analysis using quantum chemical calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Insights into Molecular Magnetism in Metal–Metal Bonded Systems as Revealed by a Spectroscopic and Computational Analysis of Diiron Complexes

A pair of bimetallic compounds featuring Fe–Fe bonds, [Fe( i PrNPPh 2 ) 3 FeR] (R = PMe 3 , ≡N t Bu), have been investigated using High-Frequency Electron Paramagnetic Resonance (HFEPR) as well as field- and temperature-dependent 57 Fe nuclear γ resonance (Mössbauer) spectroscopy. To gain insight into the local site electronic structure, we have concurrently studied a compound containing a single Fe(II) in a geometry analogous to that of one of the dimer sites. Our spectroscopic studies have allowed for the assessment of the electronic structure via the determination of the zero-field splitting and 57 Fe hyperfine parameters for the entire series. We also report on our efforts to correlate structure with physical properties in metal–metal bonded systems using ligand field theory guided by quantum chemical calculations. Through the insight gained in this study, we discuss strategies for the design of single-molecule magnets based on polymetallic compounds linked via direct metal–metal bonds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Conformational diversity of 1-phenylpiperidin-4-one in the gas phase

Conformational preferences of substituted, saturated cyclic ketones may deviate from those of their parent analogues. 1-Phenylpiperidin-4-one was found to exist in several forms differing by not only position (axial or equatorial) of the aryl-substituent but also by a configuration of the piperidine cycle itself: chair-Eq, chair-Ax and twist-Tw. Here, the ratios of the co-existing species contributions in gas phase predicted by various quantum chemical calculations were confirmed by synchronous gas-phase electron diffraction and mass spectrometric experiments carried out for the vapor at 337 K demonstrating a best fit with the B3LYP-D3/cc-pVTZ combination: Eq:Ax:Tw=55(13):22(9):23(10) vs. 40:35:25.

1-phenylpiperidin-4-one↗

Electronic structure of NdO via slow photoelectron velocity-map imaging spectroscopy of NdO –

Electronically excited NdO is a possible product of the chemistry associated with the release of Nd into the ionosphere, and emission from these states may contribute to the observations following such experiments. To better characterize the energetics and spectroscopy of NdO, we report a combined experimental and theoretical study using slow photoelectron velocity-map imaging spectroscopy of cryogenically cooled NdO – anions (cryo-SEVI) supplemented by wave function-based quantum-chemical calculations. Using cryo-SEVI, we measure the electron affinity of NdO to be 1.0091(7) eV and resolve numerous transitions to low-lying electronic and vibrational states of NdO that are assigned with the aid of the electronic structure calculations. Additionally, temperature-dependent data suggest contributions from the (2)4.5 state of NdO – residing 2350 cm –1 above the ground anion state. Finally, photodetachment to higher-lying excited states of NdO is also reported, which may help to clarify observations from prior release experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-throughput ab initio reaction mechanism exploration in the cloud with automated multi-reference validation

Quantum chemical calculations on atomistic systems have evolved into a standard approach to studying molecular matter. These calculations often involve a significant amount of manual input and expertise, although most of this effort could be automated, which would alleviate the need for expertise in software and hardware accessibility. Here, we present the AutoRXN workflow, an automated workflow for exploratory high-throughput electronic structure calculations of molecular systems, in which (i) density functional theory methods are exploited to deliver minimum and transition-state structures and corresponding energies and properties, (ii) coupled cluster calculations are then launched for optimized structures to provide more accurate energy and property estimates, and (iii) multi-reference diagnostics are evaluated to back check the coupled cluster results and subject them to automated multi-configurational calculations for potential multi-configurational cases. All calculations are carried out in a cloud environment and support massive computational campaigns. Key features of all components of the AutoRXN workflow are autonomy, stability, and minimum operator interference. We highlight the AutoRXN workflow with the example of an autonomous reaction mechanism exploration of the mode of action of a homogeneous catalyst for the asymmetric reduction of ketones.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theoretical Study of the Phenoxy Radical Recombination with the O( 3 P) Atom, Phenyl plus Molecular Oxygen Revisited

Quantum chemical calculations of the C 6 H 5 O 2 potential energy surface (PES) were carried out to study the mechanism of the phenoxy + O( 3 P) and phenyl + O 2 reactions. CASPT2(15e,13o)/CBS//CASSCF(15e,13o)/DZP multi-reference calculations were utilized to map out the minimum energy path for the entrance channels of the phenoxy + O( 3 P) reaction. Stationary points on the C 6 H 5 O 2 PES were explored at the CCSD(T)-F12/cc-pVTZ-f12//B3LYP/6-311++G** level for the species with a single-reference character of the wave function and at the CASPT2(15e,13o)/CBS//B3LYP/6-311++G** level of theory for the species with a multi-reference character of the wave function. Conventional, variational, and variable reaction coordinate transition state theories were employed in Rice-Ramsperger-Kassel-Marcus Master Equation calculations to assess temperature- and pressure-dependent phenomenological rate constants and product branching ratios. Here, the main bimolecular product channels of the phenoxy + O( 3 P) reaction are concluded to be para/ortho-benzoquinone + H, 2,4-cyclopentadienone + HCO, and, at high temperatures, also phenyl + O 2 . The main bimolecular product channels of the phenyl + O 2 reaction include 2,4-cyclopentadienone + HCO at lower temperatures, and phenoxy + O( 3 P) at higher temperatures. Both for the phenoxy + O( 3 P) and phenyl + O 2 reactions, the collisional stabilization of peroxybenzene at low temperatures and high pressures competes with the bimolecular product channels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Elaborated thermochemical treatment of HF, CO, N 2 , and H 2 O: Insight into HEAT and its extensions

Here, empirical, highly accurate non-relativistic electronic total atomization energies (eTAEs) are established by combining experimental or computationally converged treatments of the nuclear motion and relativistic contributions with the total atomization energies of HF, CO, N 2 , and H 2 O obtained from the Active Thermochemical Tables. These eTAEs, which have estimated (2σ) uncertainties of less than 10 cm -1 (0.12 kJ mol -1 ), form the basis for an analysis of high-level ab initio quantum chemical calculations that aim at reproducing these eTAEs for the title molecules. The results are then employed to analyze the performance of the high-accuracy extrapolated ab initio thermochemistry, or High-Accuracy Extrapolated Ab Initio Thermochemistry (HEAT), family of theoretical methods. The method known as HEAT-345(Q), in particular, is found to benefit from fortuitous error cancellation between its treatment of the zero-point energy, extrapolation errors in the Hartree-Fock and coupled cluster contributions, neglect of post-(T) core-correlation, and the basis-set error involved in higher-level correlation corrections. In addition to shedding light on a longstanding curiosity of the HEAT protocol—where the cheapest HEAT-345(Q) performs comparably to the theoretically more complete HEAT-456QP procedure—this study lays the foundation for extended HEAT variants that offer substantial improvements in accuracy relative to the established approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Advanced Method Optimization for Sampling and Analysis Instrumentation

This work presents a generalized approach for analytical method optimization that branches the gap between techniques historically employed and accurate modern optimization techniques suitable for various applications. The novelty of the described strategy is the utilization of multivariate, multiobjective optimization with Karush-Kuhn-Tucker conditions to bound the optimization space to solutions within the physical limitations of instrumentation. Briefly, the basic steps outlined in this paper are to (1) determine the objective(s) that should be maximized or minimized based on the goals of the analytical application, (2) conduct a screening experiment, (3) perform ANOVA to determine the parameters which have a statistically significant effect on the objective, (4) conduct an experiment (e.g., Box-Behnken design) to collect data for fitting the objective equation, and (5) determine the physical constraints of the parameters and solve the Lagrangian to determine the optimal method parameters. A broad approach to optimization target selection allows for robust method tuning to develop improved data sets amenable for chemometrics and machine learning algorithm development. Gas chromatography-mass spectrometry was selected as a use case due to its broad use across scientific fields and time-consuming method development involving numerous parameters. In conclusion, this strategy can reduce the cost of research, improve data quality, and enable the rapid development of new analytical technique.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗