Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Potential energy surfaces”

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 NTRS1987

Potential energy surfaces for atomic oxygen reactions: Formation of singlet and triplet biradicals as primary reaction products with unsaturated organic molecules

The experimental study of the interaction of atomic oxygen with organic polymer films under LEO conditions has been hampered by the inability to conduct detailed experiments in situ. As a result, studies of the mechanism of oxygen atom reactions have relied on laboratory O-atom sources that do not fully reproduce the orbital environment. For example, it is well established that only ground electronic state O atoms are present at LEO, yet most ground-based sources are known to produce singlet O atoms and molecules and ions in addition to O(3P). Engineers should not rely on such facilities unless it can be demonstrated either that these different O species are inert or that they react in the same fashion as ground state atoms. Ab initio quantum chemical calculations have been aimed at elucidating the biradical intermediates formed during the electrophilic addition of ground and excited-state O atoms to carbon-carbon double bonds in small olefins and aromatic molecules. These biradicals are critical intermediates in any possible insertion, addition and elimination reaction mechanisms. Through these calculations, we will be able to comment on the relative importance of these pathways for O(3P) and O(1D) reactions. The reactions of O atoms with ethylene and benzene are used to illustrate the important features of the mechanisms of atomic oxygen reaction with unsaturated organic compounds and polymeric materials.

Jaffe, Richard L.↗
NASA NTRS1992

Computed potential energy surfaces for chemical reactions

The work on the NH + NO system which was described in the last progress report was written up and a draft of the manuscript is included in the appendix. The appendix also contains a draft of a manuscript on an Ar + H + H surface. New work which was completed in the last six months includes the following: (1) calculations on the (1)CH2 + H2O, H2 + HCOH, and H2 + H2CO product channels in the CH3 + OH reaction; (2) calculations for the NH2 + O reaction; (3) calculations for the CH3 + O2 reaction; and (4) calculations for CH3O and the two decomposition channels--CH2OH and H + H2CO. Detailed descriptions of this work will be given in manuscripts; however, brief descriptions of the CH3 + OH and CH3 + O2 projects are given.

Heinemann, K.↗
NASA NTRS1993

Computed Potential Energy Surfaces for Chemical Reactions

A manuscript describing the calculations on the (1)CH2 + H2O, H2 + HCOH, and H2 + H2CO product channels in the CH3 + OH reaction, which were described in the last progress report, has been accepted for publication in J. Chem. Phys., and a copy of the manuscript is included in the appendix. The production of (1)CH2 in this reaction is important in hydrocarbon combustion since (1)CH2 is highly reactive and would be expected to insert into N2, possibly leading to a new source for prompt NO(x) (vide infra). During the last six months new calculations have been carried out for the NH2 + NO system, which is important in the thermal de-NO(x) process.

Heinemann, K.↗
NASA NTRS1998

Molecular Processes in Comets

The research focused on molecular hydrogen and its response to ultraviolet radiation, photoelectron impact excitation and X-ray radiation and on the interpretation of the ultraviolet spectra of the Jupiter dayglow and auroras. A systematic effort was made to obtain reliable rate coefficients for rovibrational energy transfer of H2, particularly in collisions with hydrogen atoms. We carried out elaborate quantum-mechanical calculations of the scattering of H and H2 on what had been characterized as a reliable potential energy surface of the H3 molecular system. These calculations took into account reactive channels and rate coefficients for ortho-para transitions were obtained. Quantal calculations are too complex to be applied to all the possible rovibrational transitions and we turned to a semi-classical method. With it we calculated rate coefficients for transitions, reactive and non-reactive, for all the rovibrational levels. We carried out the calculations for three of the available H3 potential energy surfaces. We discovered an unexpected sensitivity of the rate coefficients for the non-reactive channels to the potential energy surface. This discovery stimulated more thorough investigations of the potential energy surface elsewhere and to the construction of a new surface. We have used it in further semi-classical calculations in work that will shortly be 2 completed and which, together with new quantum-mechanical calculations should comprise a set of reliable rate coefficients that can be used in discussions of H2 on the Jovian planets. We carried out a detailed study of the Jovian ultraviolet dayglow. There has been a long-running argument about the dayglow on Jupiter. There are two sources of excitation: fluorescence and photoelectron impact excitation. It had been argued that a third source "the electron glow" was needed to bring theory and observation into agreement. We believe we have shown conclusively that the third source is unnecessary. We have achieved a close quantitative agreement between the predicted spectrum arising from fluorescence and photoelectrons and the measured spectrum. We also demonstrated a method by which the presence of HD could be established observationary. Similar calculations were carried out of the ultraviolet spectra of the Jovian auroras. We again found close agreement. Indeed, the agreement was so detailed we were able to derive the temperature of the atmosphere. We found, contrary to the standard model of the time, a high temperature between 400 K and 600 K and established the presence of a significant temperature gradient. A large temperature gradient was indeed found by direct sampling with the Jupiter probe.

Dalgarno, A.↗
NASA NTRS1969

Semiempirical energy surface for H sub 3

Semiempirical potential energy surface generation for triatomic hydrogen from London equation by evaluating Coulomb and exchange integrals taking into account effective orbital overlap

Salomon, M.↗
NASA NTRS1993

The quantum dynamics of electronically nonadiabatic chemical reactions

Considerable progress was achieved on the quantum mechanical treatment of electronically nonadiabatic collisions involving energy transfer and chemical reaction in the collision of an electronically excited atom with a molecule. In the first step, a new diabatic representation for the coupled potential energy surfaces was created. A two-state diabatic representation was developed which was designed to realistically reproduce the two lowest adiabatic states of the valence bond model and also to have the following three desirable features: (1) it is more economical to evaluate; (2) it is more portable; and (3) all spline fits are replaced by analytic functions. The new representation consists of a set of two coupled diabatic potential energy surfaces plus a coupling surface. It is suitable for dynamics calculations on both the electronic quenching and reaction processes in collisions of Na(3p2p) with H2. The new two-state representation was obtained by a three-step process from a modified eight-state diatomics-in-molecules (DIM) representation of Blais. The second step required the development of new dynamical methods. A formalism was developed for treating reactions with very general basis functions including electronically excited states. Our formalism is based on the generalized Newton, scattered wave, and outgoing wave variational principles that were used previously for reactive collisions on a single potential energy surface, and it incorporates three new features: (1) the basis functions include electronic degrees of freedom, as required to treat reactions involving electronic excitation and two or more coupled potential energy surfaces; (2) the primitive electronic basis is assumed to be diabatic, and it is not assumed that it diagonalizes the electronic Hamiltonian even asymptotically; and (3) contracted basis functions for vibrational-rotational-orbital degrees of freedom are included in a very general way, similar to previous prescriptions for locally adiabatic functions in various quantum scattering algorithms.

Truhlar, Donald G.↗
DOE OSTI2026

Pressure-dependent kinetics of cyclopentene consumption

Recent experiments have uncovered deficiencies in the ability of current kinetic models for cyclopentane combustion to predict species profiles as a function of O2 concentration. One likely source of these discrepancies flows from the limited availability of relevant kinetic data for key functionalized intermediates, requiring the use of distant analogies to estimate their rates. This work fills this gap for cyclopentene, the most prominent C5 intermediate in low-temperature cyclopentane combustion. Ab initio transition-state-theory based master equation methods are used to calculate pressure-dependent rate constants across four potential energy surfaces (C5H8 + OH, C5H8 + HO2, C5H7, and C5H7 + O2). Potential energy surfaces are characterized at the CCSD(T)-F12/cc-pVDZ-F12//revDSD-PBEP86-D3BJ/def2-TZVPP level of theory. Pressure-dependent rate constants and branching fractions are determined via 1D master equation calculations. The results are used to analyze the competition for radicals between cyclopentene and cyclopentane in the intermediate temperature region, and to identify the dominant pathways flowing from these initiation channels. On the OH initiation surface, we find good agreement with recently published experimental rate constants at high temperature, and we elucidate the complex temperature dependence below 800 K, where 𝜋𝜋 addition becomes prominent. On the HO2 initiation surface, we find an interesting alternative route to bicyclic ether formation via well-skipping over the HO2 adduct well. On the cyclopentenyl radical surface, the major product is cyclopentadiene plus H. On the cyclopentenyl + O2 surface, we find that resonance stabilization largely prevents the allylic radical from undergoing oxidation, whereas the alkylic radical has a somewhat elevated propensity for oxidation pathways relative to cyclopentyl.

ab initio calculations↗
NASA NTRS2017

The Effect of the Spin-Forbidden Co((sup 1) Sigma plus) plus O((sup 3) P) Yields CO2 (1 Sigma (sub G) plus) Recombination Reaction on Afterbody Heating of Mars Entry Vehicles

Vibrationally excited CO2, formed by two-body recombination from CO((sup 1) sigma plus) and O((sup 3) P) in the wake behind spacecraft entering the Martian atmosphere reaction, is potentially responsible for the higher than anticipated radiative heating of the backshell, compared to pre-flight predictions. This process involves a spin-forbidden transition of the transient triplet CO2 molecule to the longer-lived singlet. To accurately predict the singlet-triplet transition probability and estimate the thermal rate coefficient of the recombination reaction, ab initio methods were used to compute the first singlet and three lowest triplet CO2 potential energy surfaces and the spin-orbit coupling matrix elements between these states. Analytical fits to these four potential energy surfaces were generated for surface hopping trajectory calculations, using Tully's fewest switches surface hopping algorithm. Preliminary results for the trajectory calculations are presented. The calculated probability of a CO((sup 1) sigma plus) and O((sup 3) P) collision leading to singlet CO2 formation is on the order of 10 (sup -4). The predicted flowfield conditions for various Mars entry scenarios predict temperatures in the range of 1000 degrees Kelvin - 4000 degrees Kelvin and pressures in the range of 300-2500 pascals at the shoulder and in the wake, which is consistent with a heavy-particle collision frequency of 10 (sup 6) to 10 (sup 7) per second. Owing to this low collision frequency, it is likely that CO((sup 1) sigma plus) molecules formed by this mechanism will mostly be frozen in a highly nonequilibrium rovibrational energy state until they relax by photoemission.

CO2 recombination↗
NASA NTRS2002

Theoretical Study of Tautomerization Reactions for the Ground and First Excited Electronic States of Adenine

Geometrical structures and energetic properties for different tautomers of adenine are calculated in this study, using multi-configurational wave functions. Both the ground and the lowest singlet excited state potential energy surfaces are studied. Four tautomeric forms are considered, and their energetic order is found to be different on the ground and the excited state potential energy surfaces. Minimum energy reaction paths are obtained for hydrogen atom transfer (tautomerization) reactions in the ground and the lowest excited electronic states. It is found that the barrier heights and the shapes of the reaction paths are different for the ground and the excited electronic states, suggesting that the probability of such tautomerization reaction is higher on the excited state potential energy surface. This tautomerization process should become possible in the presence of water or other polar solvent molecules and should play an important role in the photochemistry of adenine.

Salter, Latasha M.↗
NASA NTRS1995

The Heat of Formation of HNO

The HNO molecule is of interest in both combustion and atmospheric chemistry. For example, Guadagnini et al. have recently presented ab initio potential energy surfaces for the three lowest lying electronic states of HNO and then used these in examining several chemical reactions that take place in the combustion of nitrogen containing fuels and in the oxidation of atmospheric nitrogen. We have previously studied the ground state potential energy surface (i.e., stationary points along the HNO rev. reaction HON path), vibrational spectrum (using an accurate quartic force field), zero-point energy, and bonding of HNO using coupled-cluster ab initio methods. HNO is also very interesting because of the unique nature of its bonding characteristics. That is, the potential energy surface is very flat along the H-N bonding coordinate thereby giving unusual harmonic and fundamental vibrational frequencies, and the H-N bond energy is rather weak in comparison to other H-N bond energies. In fact, using experimental heats of formation for HO, H, and NO, the H- bond energy is computed to be only 49.9 kcal/ mol (298 K). However, ab initio calculations of isodesmic reaction energies involving HNO, FNO, ClNO, and several other molecules have shown that there is an inconsistency in the experimental heats of formation of the XNO (X=H, F, and Cl) species. Hence the motivation for this study was to determine a very accurate(DELTA)H(sup o)(sub f) value for HNO using state of-the-art ab initio methods. Based on many recent studies it is evident that the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T), in conjunction with large one-particle basis sets should be reliable to better than +0.8 kcal/mol for this quantity. The computational methodology is described in the next section followed by our results and discussion. Conclusions are presented in the final section.

Lee, Timothy J.↗
NASA NTRS1995

The Heat of Formation of HNO

The HNO molecule is of interest in both combustion and atmospheric chemistry. For example, Guadagnini et al. have recently presented ab initio potential energy surfaces for the three lowest lying electronic states of HNO and then used these in examining several chemical reactions that take place in the combustion of nitrogen containing fuels and in the oxidation of atmospheric nitrogen. We have previously studied the ground state potential energy surface (i.e., stationary points along the HNO reversible reaction HON path), vibrational spectrum (using an accurate quartic force field), zero-point energy, and bonding of HNO using coupled-cluster ab initio methods. HNO is also very interesting because of the unique nature of its bonding characteristics. That is, the potential energy surface is very flat along the H-N bonding coordinate thereby giving unusual harmonic and fundamental vibrational frequencies, and the H-N bond energy is rather weak in comparison to other H-N bond energies. In fact, using experimental heats of formation for HNO, H, and NO, the H- bond energy is computed to be only 49.9 kcal/ mol (298 K). However, ab initio calculations of isodesmic reaction energies involving HNO, FNO, ClNO, and several other molecules have shown that there is an inconsistency in the experimental heats of formation of the XNO (X double bond H, F, and Cl) species. Hence the motivation for this study was to determine a very accurate (Delta)H(sub f, sup o) value for HNO using state-of-the-art ab initio methods. Based on many recent studies it is evident that the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T), in conjunction with large one-particle basis sets should be reliable to better than +/- 0.8 kcal/mol for this quantity. The computational methodology is described in the next section followed by our results and discussion. Conclusions are presented in the final section.

Lee, Timothy J.↗
DOE OSTI2024

Full‐dimensional coupled‐channel statistical approach to atom‐triatom systems and applications to H/D + O 3 reaction

Abstract The statistical quantum model (SQM), which assumes that the reactivity is controlled by entrance/exit channel quantum capture probabilities, is well suited for chemical reactions with a long‐lived intermediate complex. In this work, a time‐independent coupled‐channel implementation of the SQM approach is developed for atom‐triatom systems in full dimensionality. As SQM treats the capture dynamics quantum mechanically, it is capable of handling quantum effects such as tunneling. A detailed study of the H/D + O 3 capture dynamics was performed by applying the newly developed SQM method on an accurate global potential energy surface. Agreement with previous ring polymer molecular dynamics (RPMD) results on the same potential energy surface is excellent except for very low temperatures. The SQM results are also in reasonably good agreement with available experimental rate coefficients. The strong H/D kinetic isotope effect underscores the dominant role of quantum tunneling under an entrance channel barrier at low temperatures.

Yang, Dongzheng↗
DOE OSTI2025

Ab-initio simulation of spin-vibronic spectra of methoxy radical

Despite the fact that experimental and theoretical work on the spectrum of methoxy has stretched from the microwave to the ultraviolet and proceeded for nearly 50 years, parts of the spectrum have remained a challenge to simulate theoretically and make reliable line-by-line assignments. The spectral complexity arises because the radical has a non-zero electron spin and significant vibronic coupling between the two elec- tronic components of the ground state due to the presence of a conical intersection. This work describes a completely ab initio effort to understand and assign the spin- vibronic levels of the X 2E state from 0 to above 3000 cm−1, a region that includes the fundamental transitions of the C-H symmetric and asymmetric stretches that have not previously been identified uniquely. A potential energy surface for methoxy was calculated at the EOM-CCSDT/ANO1 level of theory. Subsequently this potential energy surface was fit to a quartic power series expansion of all nine vibrational nor- mal coordinates (as determined at the minimum of the conical intersection) by the use of a machine-learning-based algorithm. After the addition of spin-orbit coupling, the spin-vibronic problem was solved using both the Krylov-Schur and Lanczos algorithms with the SOCJT3 software to converge eigenvalues up to 3500 cm−1 and their eigen- vectors. The latter were used, in conjunction with the calculated dipole moment and its derivatives (calculated using finite differences at EOM-CCSDT/ANO1 level), to determine spectral intensities for the spin-vibronic spectra. The calculated transition frequencies and intensities were used to simulate and assign the observed transitions of the spin-vibronic spectra of the radical. The credibility of the assignments and their significance is discussed in detail.

Sharma, Ketan [University of Florida, Gainesville,↗
NASA NTRS2002

Effect of Electronic Excitation on Hydrogen Atom Transfer (Tautomerization) Reactions for the DNA Base Adenine

Geometrical structures and energetic properties for four different tautomers of adenine are calculated in this study, using multi-configurational wave functions. Both the ground and the lowest single excited state potential energy surface are studied. The energetic order of the tautomers on the ground state potential surface is 9H less than 7H less than 3H less than 1H, while on the excited state surface this order is found to be different: 3H less than 1H less than 9H less than 7H. Minimum energy reaction paths are obtained for hydrogen atom transfer (9 yields 3 tautomerization) reactions in the ground and the lowest excited electronic state. It is found that the barrier heights and the shapes of the reaction paths are different for the ground and the excited electronic state, suggesting that the probability of such tautomerization reaction is higher on the excited state potential energy surface. The barrier for this reaction in the excited state may become very low in the presence of water or other polar solvent molecules, and therefore such tautomerization reaction may play an important role in the solution phase photochemistry of adenine.

Chaban, Galina M.↗