Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “potential energy surface”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Theoretical studies of diatomic and triatomic systems containing the group IB atoms Cu, Ag, and Au

Selected portions of the ground state potential energy surfaces of the Cu3, Ag3, AgCu2, and AuCu2 trimers are studied at the single-reference singles plus doubles configuration interaction and couple pair functional levels correlating 33 electrons. The calculations use the effective core potentials of Hay and Wadt (1985) to replace the deep core levels. The Cu3 and Ag3 molecules are found to have 2B2 obtuse-angled ground states, with low-lying 2A1 acute-angled excited states. The AgCu2 and AuCu2 molecules have 2A1 acute-angled ground states. The Cu3 molecule has a smaller 3d population than Ag3, and Cu3 has a smaller s electron density in the open-shell orbital than does Ag3, in agreement with recent ESR experiments.

Walch, Stephen P.↗
NASA NTRS1988

An improved long range potential for O(1D)+H2

Several dynamics studies have indicated that the product isotopic distributions and vibrational energy distributions in the reaction O(1D) + H2(HD) are very sensitive to the long range interaction of the reactants. In this paper, an improved calculation of the reactant region of the potential energy surface is reported. In agreement with previous work, no barrier is found to edge-on insertion (1 A-prime ground state water surface), but a smaller (less than 0.2 kcal/mol) barrier to collinear addition, 1 Sigma (+) surface, is found than in previous work. The long range potential obtained in the present work most closely resembles the SL3 surface.

Walch, Stephen P.↗
NASA NTRSDate not supplied

Quantum Hardware-Enabled Molecular Dynamics via Transfer Learning

The ability to perform ab initio molecular dynamics simulations using potential energy surfaces provided by quantum computers would open the door to virtually exact dynamics for a variety of chemical and biochemical systems, with impacts on catalysis and biophysics. Nonetheless, performing molecular dynamics on surfaces produced by quantum hardware has been hampered by the noisy energies typically produced by quantum computers and challenges associated with computing gradients and scaling to large systems interest. A recent set of advances in machine learning, known as transfer learning, provides a new path forward for molecular dynamics simulations on quantum hardware. Transfer learning offers a workaround, where one first trains models on larger, less accurate classical datasets and then refines them on smaller, more accurate quantum datasets. We explore this approach by training machine learning models to predict a molecule's potential energy based on its geometric structure using Behler-Parrinello neural networks. When successfully trained, the model enables energy gradient predictions necessary for dynamic simulations. To reduce the quantum resources needed, the model is initially trained with data derived from classical density functional theory and subsequently refined with a smaller dataset obtained from a variational quantum eigensolver optimization of the unitary coupled cluster ansatz. We show that this approach significantly reduces the size of the needed quantum training dataset while capturing the high accuracies needed within quantum chemistry simulations. The success of this two-step training method opens more opportunities to apply machine learning models on quantum data, a significant stride towards efficient quantum-classical hybrid computational models.

quantum computing↗
NASA NTRS2015

The Structure, Anharmonic Vibrational Frequencies, and Intensities of NNHNN+

A semi-global potential energy surface (PES) and quartic force field (QFF) based on fitting high-level electronic structure energies are presented to describe the structures and spectroscopic properties of NNHNN+. The equilibrium structure of NNHNN+ is linear with the proton equidistant between the two nitrogen groups and thus of D(sub ∞h) symmetry. Vibrational second-order perturbation theory (VPT2) calculations based on the QFF fails to describe the proton “rattle” motion, i.e., the antisymmetric proton stretch, due to the very flat nature of PES around the global minimum, but performs properly for other modes with sharper potential wells. Vibrational self-consistent field/virtual state configuration interaction (VSCF/VCI) calculations using a version of MULTIMODE without angular momentum terms successfully describe this motion and predict the fundamental to be at 759 cm(exp -1). This is in good agreement with the value of 746 cm(exp -1) from a fixed-node diffusion Monte Carlo calculation and the experimental Ar-tagged result of 743 cm(exp -1). Other VSCF/VCI energies are in good agreement with other experimentally reported ones. Both double-harmonic intensity and rigorous MULTIMODE intensity calculations show the proton transfer fundamental has a very strong intensity.

Astrochemistry↗
NASA NTRS2022

Ames-2021 CO 2 Dipole Moment Surface and IR Line Lists: Toward 0.1% Uncertainty for CO 2 IR Intensities

A highly accurate CO 2 ab initio dipole moment surface (DMS), Ames-2021, is reported along with 12 C 16 O 2 infrared (IR) intensity comparisons approaching a 1–4‰ level of agreement and uncertainty. The Ames-2021 DMS was accurately fitted from CCSD(T) finite-field dipoles computed with the aug-cc-pVXZ (X = T, Q, 5) basis for C atom and the d-aug-cc-pVXZ (X = T, Q, 5) basis for O atoms, and extrapolated to the one particle basis set limit. Fitting σ(rms) is 3.8 × 10^(–7) au for 4443 geometries below 15 000 per cm . The corresponding IR intensity, S Ames-2021 , are computed using the Ames-2 potential energy surface (PES), which is the best PES available for CO 2 . Compared to high accuracy IR studies for 2001i–00001 and 3001i–00001 bands, S Ames-2021 matches NIST experiment-based intensities [S NIST-HIT16 or S HIT20 ] to −1.0 ± 1.3‰, or matches DLR experiment-based intensities [S DLR-HIT16/UCL/Ames ] to 1.9 ± 3.7‰. This indicates the systematic deviations and uncertainties have been significantly reduced in S Ames-2021 . The S UCL2015 (or S HITRAN2016 ) have larger deviations (vs S DLR ) and uncertainties (vs S DLR , S NIST ) which are attributed to the less accurate Ames-1 PES adopted in UCL-296 line list calculation. The S Ames-2021 intensity of 12 C 16 O 2 and 13 C 16 O 2 is utilized to derive new absolute 13 C/ 12 C ratios for Vienna PeeDee Belemnite (VPDB) with uncertainty reduced by 1/3 or 2/3. Further evaluation of S Ames-2021 are carried out on those CO 2 bands discussed in the HITRAN2020 update paper. Consistent improvements and better accuracies are found in band-by-band analysis, except for those bands strongly affected by Coriolis couplings, or very weak bands measured with relatively larger experimental uncertainties. The Ames-2021 296 K IR line lists are generated for 13 CO 2 isotopologues, with 18,000 per and S 296K > 1E–31) cm/molecule cutoff and then combined with CDSD line positions (except 14 C 16 O 2 ). The Ames-2021 DMS and 296 K IR line lists represent a major improvement over previous CO 2 theoretical IR intensity studies, including Ames-2016, UCL-296, and recent UCL DMS 2021 update. A real 1 permille level of agreement and uncertainty will definitely require both more accurate PES and more accurate DMS.

Chemical calculations↗
NASA NTRS1985

The lower electronic states of ClOO - A computational investigation

Eight doublet and eight quartet states of ClOO were investigated by ab initio CI techniques. The potential energy surfaces of the four lowest energy doublet states of both A-double-prime and A-prime symmetry indicate that only the 1 2A-double-prime state is bound. In contrast to the model provided by the HO2 radical, all of the excited doublet states investigated were repulsive with respect to dissociation to Cl + O2 and metastable or bound with respect to dissociation to ClO + O. The transitions to the excited states investigated span the visible and near UV spectral regions, but the transition moments indicate that they are very weak. Since the photolysis products are the same as those of the rapid thermal dissociation, photolysis is not expected to be an important atmospheric process. The soft bending potential for the 1 2A-prime state and the shape of the 1 4A-double-prime state in the entrance channel of the ClO + O yields Cl + O2 reaction provide a qualitative explanation for the underprediction of the low temperature reaction rate by previous trajectory calculations.

Jafri, J. A.↗
NASA NTRS2022

Recommended Direct Simulation Monte Carlo Collision Model Parameters for Reacting Methane Flows

While hydrocarbon combustion and gas-surface interactions have been simulated for a variety of applications in continuum solvers, such processes have been rarely investigated in Direct Simulation Monte Carlo (DSMC). Thus Variable Hard/Soft Sphere (VHS/VSS)collision parameters are not available for most molecules of interest in the decomposition of methane, a significant pyrolysis gas product and potential atmospheric species for Titan re-entries. As these properties are essential for accurate modeling of gas transport and thermochemical effects, a method must be devised to compute them for use in DSMC. In this study, the methodology outlined in Stephani et al is utilized to compute the VHS/VSS parameters by fitting the Ω(1,1)and Ω(2,2) collision integrals [1]. These collision integrals are typically computed from intermolecular potentials and are used in CFD simulations to calculate transport properties such as viscosity, thermal conductivity, and diffusion. This approach ensures consistency between micro-scale collisions and macro-scale transport properties. To produce the necessary collision integrals needed for the fitting routine, we utilize the method described by Laricchiuta et al, where collision integrals are computed based on the constituent particles’ polarizability and effective electron number [2].Preliminary computations of pure-species Ω(1,1)and Ω(2,2)collision integrals produced from the coupled approach of Stephani et al and Laricchiuta et al for methane show good agreement with those computed via available ab-initio Potential Energy Surface Computations, as seen in Figure 1. Above 300K, no more than 12% deviation is observed between the two results over the temperature range of the fitted PES data [3].Current results of this study include full collision-specificVSS transport properties for a complex-species mixture involving major byproducts of pure methane decomposition. Targeted final results will include a comprehensive list of VSS parameters relevant toDSMC simulations of Titan atmospheric entry, based on the reduced kinetic model of Savajano et al [4].

M R Gosma↗
NASA NTRS1972

Dynamics of the molecular and atomic mechanisms for the hydrogen-iodine exchange reaction.

The molecular and atomic mechanisms for the hydrogen-iodine exchange reaction are treated theoretically by means of extensive classical trajectories calculated on a reasonable potential energy surface on which the single adjustable parameter is the iodine-core effective charge. The analysis shows the molecular mechanism to be dynamically forbidden, but gives an over-all rate constant for the atomic mechanism that is in agreement with the experimental values. It is indicated that the formation of a weak H2I complex plays an important dynamical role if the atomic mechanism is limited to reactions with collision complexes involving no more than two hydrogen atoms and two iodine atoms. Excellent agreement with experiment is obtained for the rate constant for the recombination I+I+H2 yields I2+H2 and its negative temperature coefficient.

Raff, L. M.↗
NASA NTRS1972

Stabilities of nitrogen containing heterocyclic radicals and geometrical influences on non-radiative processes in organic molecules

Stabilities of nitrogen containing heterocyclic radicals were studied to detect radicals of the type R-N-R, and to theoretically rationalize their electronic structure. The computation of simple potential energy surfaces for ground and excited states is discussed along with the photophysical properties of indolizine. Methods of calculation and problems associated with the calculations are presented. Results, tables, diagrams, discussions, and references are included.

Evleth, E. M.↗
NASA NTRS1977

Rotational and vibrational transitions for Li + H2 collisions

Close coupling calculations for integral and differential cross sections have been carried out for Li + H2 collisions with an ab initio Hartree-Fock potential energy surface. Rotational, vibrational, and vib-rotational excitation cross sections are reported at 0.4336 eV, 0.7 eV, and 0.8673 eV in the center of mass system. For pure rotational excitations, which dominate the inelastic scattering, coupling with vibrational states is not very important. For vibrational transitions, the influence of large multiquantum rotational transitions is far less than that found for Li(+) + H2 collisions.

Choi, B. H.↗
NASA NTRS1978

On the possibility of negative activation energies in bimolecular reactions

The temperature dependence of the rate constants for model reacting systems was studied to understand some recent experimental measurements which imply the existence of negative activation energies. A collision theory model and classical trajectory calculations are used to demonstrate that the reaction probability can vary inversely with collision energy for bimolecular reactions occurring on attractive potential energy surfaces. However, this is not a sufficient condition to ensure that the rate constant has a negative temperature dependence. On the basis of these calculations, it seems unlikely that a true bimolecular reaction between neutral molecules will have a negative activation energy.

Jaffe, R. L.↗
NASA NTRS1978

Theoretical study of the photodissociation of HOCl

A study of the UV photodissociation of hypochlorous acid is conducted on the basis of ab initio SCF-Cl calculations. These calculations show that HOCl has only a single peak in the UV photoabsorption spectrum at 220 nm. This result implies that HOCl would have a long lifetime for photodissociation if it were to be formed in the stratosphere. The photodissociation products of HOCl have been identified as Cl+OH, based on an examination of the topographies of the excited electronic state potential energy surfaces. The results of this study indicate that HOCl could be a significant reservoir for stratospheric chlorine.

Jaffe, R. L.↗
NASA NTRS1978

Theoretical study of collinear Be + FH/nu1/ yields BeF/nu2/ + H

The potential energy surface for collinear Be + FH yields BeF + H was studied at various levels of ab initio approximation. A final surface was obtained from a first order configuration interaction wave function, using the iterative orbital method and a medium-sized basis of Slater atomic functions. The exothermicity is computed to be 6 kcal/mole; the barrier height is predicted to be about 28 kcal/mole at a geometry where both internuclear separations are extended by about 0.4 bohr from their asymptotic equilibrium values. This surface differs qualitatively from simple LEPS models.

Schor, H.↗
NASA NTRS1981

A hybrid method for improving MCSCF convergence

It has been found that the convergence problems for many ill conditioned single-configuration SCF calculations arise from mixing among only a small number of orbitals. This orbital set includes the highest closed, the partially filled, and (possibly) a few of the lowest virtual orbitals. The fact that convergence problems can be, in very large measure, linked to a small orbital set is used to design a hybrid MCSCF procedure in which these orbitals are treated using a second-order MCSCF method, while other mixings are treated with a lower-order method which avoids the time consuming integral transformation. Tests on BeO show that the hybrid method yields convergence even when the simple lower-order treatment diverges. The method is expected to facilitate determination of MCSCF wave functions for large basis problems and for the construction of potential energy surfaces.

Bauschlicher, C. W., Jr.↗
NASA NTRS1983

A new strategy for stress analysis using the finite element method

In the paper the authors examine the effectiveness of the Powell-Toint strategy for evaluating the Hessian of the potential energy surface of a finite element model that can be used for linear stress analysis and transient response predictions of structures. Cases for which the Powell-Toint strategy may be cost-effective with the conventional method of stress analysis are identified.

Kamat, M. P.↗
NASA NTRS1984

Metastability of isoformyl ions in collisions with helium and hydrogen

The stability of HOC(+) ions under conditions in interstellar molecular clouds is considered. In particular, the possibility that collisions with helium or hydrogen will induce isomerization to the stable HCO(+) form is examined theoretically. Portions of the electronic potential energy surfaces for interaction with He and H atoms are obtained from standard quantum mechanical calculations. Collisions with He atoms are found to be totally ineffective for inducing isomerization. Collisions with H atoms are found to be ineffective at low interstellar temperatures owing to a small (about 500 K) barrier in the entrance channel; at higher temperatures where this barrier can be overcome, however, collisions with hydrogen atoms do result in conversion to the stable HCO(+) form. Although detailed calculations are not presented, it is argued that low-energy collisions with H2 molecules are also ineffective in destroying the metastable ion.

Green, S.↗
NASA NTRS1984

Hydrogen atom abstraction from aldehydes - OH + H2CO and O + H2CO

The essential features of the potential energy surfaces governing hydrogen abstraction from formaldehyde by oxygen atom and hydroxyl radical have been characterized with ab inito multiconfiguration Hartree-Fock (MCHF) and configuration interaction (CI) wave functions. The results are consistent with a very small activation energy for the OH + H2CO reaction, and an activation energy of a few kcal/mol for the O + H2CO reaction. In the transition state structure of both systems, the attacking oxygen atom is nearly collinear with the attacked CH bond.

Dupuis, M.↗