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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Highly Accurate Potential Energy Surface, Dipole Moment Surface, Rovibrational Energy Levels, and Infrared Line List for (32)S(16)O2 up to 8000 cm(exp -1)

A purely ab initio potential energy surface (PES) was refined with selected (32)S(16)O2 HITRAN data. Compared to HITRAN, the root-mean-squares error (RMS) error for all J=0-80 rovibrational energy levels computed on the refined PES (denoted Ames-1) is 0.013 cm(exp -1). Combined with a CCSD(T)/aug-cc-pV(Q+d)Z dipole moment surface (DMS), an infrared (IR) line list (denoted Ames-296K) has been computed at 296K and covers up to 8,000 cm(exp -1). Compared to the HITRAN and CDMS databases, the intensity agreement for most vibrational bands is better than 85-90%. Our predictions for (34)S(16)O2 band origins, higher energy (32)S(16)O2 band origins and missing (32)S(16)O2 IR bands have been verified by most recent experiments and available HITRAN data. We conclude that the Ames-1 PES is able to predict (32/34)S(16)O2 band origins below 5500 cm(exp -1) with 0.01-0.03 cm(exp -1) uncertainties, and the Ames-296K line list provides continuous, reliable and accurate IR simulations. The Ka-dependence of both line position and line intensity errors is discussed. The line list will greatly facilitate SO2 IR spectral experimental analysis, as well as elimination of SO2 lines in high-resolution astronomical observations.

Sulfur Dioxide (SO2)↗
NASA NTRS1979

F + H2 collisions on two electronic potential energy surfaces - Quantum-mechanical study of the collinear reaction

Collinear quantum calculations are carried out for reactive F + H2 collisions on two electronic potential energy surfaces. The resulting transmission and reflection probabilities exhibit much greater variation with energy than single-surface studies would lead us to anticipate. Transmission to low-lying product channels is increased by orders of magnitude by the presence of the second surface; however, branching ratios among product states are found to be independent of the initial electronic state of the reactants. These apparently contradictory aspects of the calculation are discussed and a tentative explanation put forward to resolve them.

Zimmerman, I. H.↗
DOE OSTI2023

Comprehensive Kinetics on the C 7 H 7 Potential Energy Surface under Combustion Conditions

The automated kinetics workflow code, KinBot, was used to explore and characterize the regions of the C 7 H 7 potential energy surface that are relevant to combustion environments and especially soot inception. We first explored the lowest-energy region, which includes the benzyl, fulvenallene + H, and cyclopentadienyl + acetylene entry points. We then expanded the model to include two higher-energy entry points, vinylpropargyl + acetylene and vinylacetylene + propargyl. The automated search was able to uncover the pathways from the literature. In addition, three important new routes were discovered: a lower-energy pathway connecting benzyl with vinylcyclopentadienyl, a decomposition mechanism from benzyl that results in side-chain hydrogen atom loss to produce fulvenallene + H, and shorter and lower energy routes to the dimethylene-cyclopentenyl intermediates. We systematically reduced the extended model to a chemically relevant domain composed of 63 wells, 10 bimolecular products, 87 barriers, and 1 barrierless channel and constructed a master equation using the CCSD(T)-F12a/cc-pVTZ//ωB97X-D/6-311++G(d,p) level of theory to provide rate coefficients for chemical modeling. Our calculated rate coefficients show excellent agreement with measured ones. We also simulated concentration profiles and calculated branching fractions from the important entry points to provide an interpretation of this important chemical landscape.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

A new diabatic representation of the coupled potential energy surfaces for Na(3p P-2) + H2 yields Na(3s S-2) + H2 or NaH + H

A diabatic representation is presented of the coupled potential-energy surfaces for Na(3p P-2) + H2 yields Na (3s S-2) + H2 or NaH + H. The representation is designed to yield, upon diagonalization, realistic values for the two lowest energy adiabatic states at both asymptotes of the chemical reaction as well as near the conical intersection in the three-body interaction region. It is economical to evaluate and portable. It is suitable for dynamics calculations on both the quenching process and the electronically nonadiabatic chemical reaction.

Halvick, Philippe↗
NASA NTRS1991

Theoretical characterization of the potential energy surface for H + N2 yields HN2. III - Calculations for the excited state surfaces

Additional calculations which characterize potential energy sources (PESs) for the excited 3A-double-prime state, for a bound 2(2A-prime) state, for HN2(+), and for the Rydberg states associated with HN2(+). It is anticipated that these excited state PESs will be important in interpreting and designing experiments to characterize the ground state HN2 species via neutralized ion beam techniques.

Walch, Stephen P.↗
NASA NTRS1991

Converged three-dimensional quantum mechanical reaction probabilities for the F + H2 reaction on a potential energy surface with realistic entrance and exit channels and comparisons to results for three other surfaces

Accurate three-dimensional quantum mechanical reaction probabilities are presented for the reaction F + H2 yields HF + H on the new global potential energy surface 5SEC for total angular momentum J = 0 over a range of translational energies from 0.15 to 4.6 kcal/mol. It is found that the v-prime = 3 HF vibrational product state has a threshold as low as for v-prime = 2.

Lynch, Gillian C.↗
DOE OSTI2020

Theoretical H + O 3 rate coefficients from ring polymer molecular dynamics on an accurate global potential energy surface: assessing experimental uncertainties

Thermal rate coefficients and kinetic isotope effects have been calculated for an important atmospheric reaction H/D + O 3 → OH/OD + O2 based on an accurate permutation invariant polynomial-neural network potential energy surface, using ring polymer molecular dynamics (RPMD), quasi-classical trajectory (QCT) and variational transition-state theory (VTST) with multidimensional tunneling. The RPMD approach yielded results that are generally in better agreement with experimental rate coefficients than the VTST and QCT ones, especially at low temperatures, attributable to its capacity to capture quantum effects such as tunneling and zero-point energy. Additionally, the theoretical results support one group of existing experiments over the other. In addition, rate coefficients for the D + O 3 → OD + O 2 reaction are also reported using the same methods, which will allow a stringent assessment of future experimental measurements, thus helping to reduce the uncertainty in the recommended rate coefficients of this reaction.

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

Exploring the Tl 2H 2 potential energy surface: A comparative analysis with group 13 systems and experiment

Abstract Thallium chemistry is experiencing unprecedented importance. Therefore, it is valuable to characterize some of the simplest thallium compounds. Stationary points along the singlet and triplet TlH potential energy surface have been characterized. Stationary point geometries were optimized with the CCSD(T)/aug‐cc‐pwCVQZ‐PP method. Harmonic vibrational frequencies were computed at the same level of theory while anharmonic vibrational frequencies were computed at the CCSD(T)/aug‐cc‐pwCVTZ‐PP level of theory. Final energetics were obtained with the CCSDT(Q) method. Basis sets up to augmented quintuple‐zeta cardinality (aug‐cc‐pwCV5Z‐PP) were employed to obtain energetics in order to extrapolate to the complete basis set limits using the focal point approach. Zero‐point vibrational energy corrections were appended to the extrapolated energies in order to determine relative energies at 0 K. It was found that the planar dibridged isomer lies lowest in energy while the linear structure lies highest in energy. The results were compared to other group 13 MH (M = B, Al, Ga, In, and Tl) theoretical studies and some interesting variations are found. With respect to experiment, incompatibilities exist.

Chemistry↗
DOE OSTI2025

Two-dimensional ring polymer molecular dynamics determination of the MnO + + H 2 /D 2 reaction rates on a Δ-machine learned potential energy surface

In this work, we investigate the impact of nuclear quantum effects in the kinetics of the MnO + + H 2 reaction, a prototypical system for gas-phase H 2 activation by transition metal oxide ions. The DFT based potential energy surfaces (PESs) for the lowest-lying quintet and septet spin states reported in our previous work are improved by 2953 newly calculated CCSD(T)/AVDZ points using a delta-machine learning (Δ-ML) method. To examine nuclear quantum effects, the rate coefficients are computed using ring-polymer molecular dynamics (RPMD). Due to the mechanistic complexity of the reaction, two reaction coordinates are necessary to map out the free-energy surface and an extended RPMD rate theory is developed. The calculated RPMD rate coefficients on the new PES are in better agreement with experimental data. The calculated kinetic isotope effects (KIEs) range from 1.6 to 1.8, also in good agreement with the experimental values. Additionally, the converged RPMD rate coefficients are 2.3 to 3.1 times higher than their classical counterparts across the studied temperature range, suggesting the presence of moderate nuclear quantum effects in this reaction.

Liu, Yang [University of New Mexico, Albuquerque, ↗
DOE OSTI2023

Comparison of multifidelity machine learning models for potential energy surfaces

Multifidelity modeling is a technique for fusing the information from two or more datasets into one model. It is particularly advantageous when one dataset contains few accurate results and the other contains many less accurate results. Within the context of modeling potential energy surfaces, the low-fidelity dataset can be made up of a large number of inexpensive energy computations that provide adequate coverage of the N-dimensional space spanned by the molecular internal coordinates. The high-fidelity dataset can provide fewer but more accurate electronic energies for the molecule in question. Here, we compare the performance of several neural network-based approaches to multifidelity modeling. We show that the four methods (dual, Δ-learning, weight transfer, and Meng–Karniadakis neural networks) outperform a traditional implementation of a neural network, given the same amount of training data. We also show that the Δ-learning approach is the most practical and tends to provide the most accurate model.

Chemistry↗
DOE OSTI2020

Incorporating Electronic Information into Machine Learning Potential Energy Surfaces via Approaching the Ground-State Electronic Energy as a Function of Atom-Based Electronic Populations

Machine learning (ML) approximations to density functional theory (DFT) potential energy surfaces (PESs) are showing great promise for reducing the computational cost of accurate molecular simulations, but at present, they are not applicable to varying electronic states, and in particular, they are not well suited for molecular systems in which the local electronic structure is sensitive to the medium to long-range electronic environment. With this issue as the focal point, we present a new machine learning approach called "BpopNN" for obtaining efficient approximations to DFT PESs. Conceptually, the methodology is based on approaching the true DFT energy as a function of electron populations on atoms; in practice, this is realized with available density functionals and constrained DFT (CDFT). Additionally, the new approach creates approximations to this function with neural networks. These approximations thereby incorporate electronic information naturally into a ML approach, and optimizing the model energy with respect to populations allows the electronic terms to self-consistently adapt to the environment, as in DFT. We confirm the effectiveness of this approach with a variety of calculations on LinHn clusters.

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

Dissociative electron attachment to 5-bromo-uracil: non-adiabatic dynamics on complex-valued potential energy surfaces

Electron induced dissociation reactions are relevant to many fields, ranging from prebiotic chemistry to cancer treatments. However, the simulation of dissociation electron attachment (DEA) dynamics is very challenging because the auto-ionization widths of the transient negative ions must be accounted for. We propose an adaptation of the ab initio multiple spawning (AIMS) method for complex-valued potential energy surfaces, along the lines of recent developments based on surface hopping dynamics. Our approach combines models for the energy dependence of the auto-ionization widths, obtained from scattering calculations, with survival probabilities computed for the trajectory basis functions employed in the AIMS dynamics. The method is applied to simulate the DEA dynamics of 5-bromouracil in full dimensionality, i.e., taking all the vibrational modes into consideration. The propagation starts on the $π_{2}^{*}$ resonance state and describes the formation of Br - anions mediated by non-adiabatic couplings. The potential energies, gradients and non-adiabatic couplings were computed with the fractional-occupancy molecular orbital complete-active-space configuration-interaction method, and the calculated DEA cross section are consistent with the observed DEA intensities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2016

Comparison of Quantum Mechanical and Empirical Potential Energy Surfaces and Computed Rate Coefficients for N2 Dissociation

Physics-based modeling of hypersonic flows is predicated on the availability of chemical reaction rate coefficients and cross sections for the collisional processes. This approach has been built around the use of quantum mechanical calculations to describe the interaction between the colliding particles. In this approach a potential energy surface (PES) is computed by solving the electronic Schrödinger equation and collision cross sections are determined for that PES using classical, semiclassical or quantum mechanical scattering methods. The rate coefficients are computed by integrating the thermally weighted cross sections. State-to-state rate coefficients are determined by only integrating over a thermal distribution of collisional energies. Finally, thermal rate coefficients are determined by summation of the state-to-state rate coefficients for reactions of molecules in all relevant ro-vibrational energy levels. If the flow is in thermal non-equilibrium, the translational, vibrational and rotational energy modes can be represented in different ways: three unique temperatures can be used to describe the distributions, the populations of individual ro-vibrational energy levels can be determined by solving the Master Equation, or through the use of direct simulation in particle-based Monte Carlo sampling. The PES-to-rate coefficient approach had been proposed and attempted in the early days of digital computing, but it is only in the last 15 years that computer hardware and software have been up to the task of calculating accurate interatomic and intermolecular potentials.

Jaffe, Richard L.↗
DOE OSTI2021

Photodissociation Dynamics of CH 2 OO on Multiple Potential Energy Surfaces: Experiment and Theory

UV excitation of the CH 2 OO Criegee intermediate across most of the broad span of the (B 1A') – (X 1A') spectrum results in prompt dissociation to two energetically accessible asymptotes: O ( 1 D) + H 2 CO (X 1 A 1 ) and O ( 3 P) + H 2 CO (a 3A"). Dissociation proceeds on multiple singlet potential energy surfaces that are coupled by two regions of conical intersection (CoIn). Velocity map imaging (VMI) studies reveal a bimodal total kinetic energy (TKER) distribution for the O ( 1 D) + H 2 CO (X 1 A 1 ) products with the major and minor components accounting for ca. 40% and ca. 20% on average of the available energy (E avl ), respectively. The unexpected low TKER component corresponds to highly internally excited H 2 CO (X 1 A 1 ) products accommodating ca. 80% of E avl . Full dimensional trajectory calculations suggest that the bimodal TKER distribution of the O ( 1 D) + H 2 CO (X 1 A 1 ) products originates from two different dynamical pathways: a primary pathway (69%) evolving through one CoIn region to products and a smaller component (20%) sampling both CoIn regions enroute to products. Those that access both CoIn regions likely give rise to the more highly internally excited H 2 CO (X 1 A 1 ) products. The remaining trajectories (11%) dissociate to O ( 3 P) + H 2 CO (a 3 A") products after traversing through both CoIn regions. Here, the complementary experimental and theoretical investigation provides insight on the photodissociation of CH 2 OO via multiple dissociation pathways through two regions of CoIn that control the branching and energy distributions of products.

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

Trapping Ca + inside a molecular cavity: computational study of the potential energy surfaces for Ca + -[ n ]cycloparaphenylene, n = 5–12

Ion trap quantum computing utilizes electronic states of atomic ions such as Ca + to encode information on to a qubit. To explore the fundamental properties of Ca + inside molecular cavities, we describe here a computational study of Ca + bound inside neutral [n]-cycloparaphenylenes (n = 5–12), often referred to as “nanohoops”. This ab initio study characterizes optimized structures, harmonic vibrational frequencies, potential energy surfaces, and ion molecular orbital distortion as functions of increasing nanohoop size. Here the results of this work provide a first step in guiding experimental studies of the spectroscopy of these ion-molecular cavity complexes.

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

Ab initio quantum scattering calculations for the CO–O 2 system and a new CO–O 2 potential energy surface: O 2 and air broadening of the R(0) line in CO

We present ab initio calculations of the collisional broadening of the R(0) pure rotational line in CO (at 115 GHz) perturbed by O 2 . Our calculations are done in a fully quantum way by solving close-coupling quantum-scattering equations without any approximations. We also report a new, highly accurate CO–O 2 potential energy surface on which we did the quantum-scattering calculations. The calculated collisional broadening agrees with the available experimental data in a wide temperature range. The calculated collisional shift is negligible compared to the broadening, which is also consistent with the experimental data. We combine this result with our previous calculations for the same line in CO perturbed by N 2 [Józwiak et al., J. Chem. Phys. 154, 054314 (2021)]; the obtained air-perturbed broadening of the R(0) pure rotational line in CO and its temperature dependence perfectly agree with the HITRAN database. This result constitutes an important step toward developing a methodology for providing accurate ab initio reference data on spectroscopic collisional line-shape parameters for molecular systems relevant to the Earth’s atmosphere and for populating spectroscopic line-by-line databases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2001

A semiempirical study of the optimized ground and excited state potential energy surfaces of retinal and its protonated Schiff base

The electronic ground and first excited states of retinal and its Schiff base are optimized for the first time using the semiempirical AM1 Hamiltonian. The barrier for rotation about the C(11)-C(12) double bond is characterized by variation of both the twist angle delta(C(10)-C(11)-C(12)-C(13)) and the bond length d(C(11)-C(12)). The potential energy surface is obtained by varying these two parameters. The calculated ground state rotational barrier is equal to 15.6 kcal/mol for retinal and 20.5 kcal/mol for its Schiff base. The all-trans conformation is more stable by 3.7 kcal/mol than the 11-cis geometry. For the first excited state, S(1,) the 90 degrees twisted geometry represents a saddle point for retinal with the rotational barrier of 14.6 kcal/mol. In contrast, this conformation is an energy minimum for the Schiff base. It can be easily reached at room temperature from the planar minima since it is separated from them by a barrier of only 0.6 kcal/mol. The 90 degrees minimum conformation is more stable than the all-trans by 8.6 kcal/mol. We are thus able to present a reaction path on the S(1) surface of the retinal Schiff base with an almost barrier-less geometrical relaxation into a twisted minimum geometry, as observed experimentally. The character of the ground and first excited singlet states underscores the need for the inclusion of double excitations in the calculations.

Retinaldehyde/chemistry↗