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

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.↗

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.↗

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.↗

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↗

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.↗

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.↗

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.↗

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.↗

Ca + HF - The anatomy of a chemical insertion reaction

A comprehensive first-principles theoretical investigation of the gas phase reaction Ca + HF - CaF + H is reported. Ab initio potential energy calculations are first discussed, along with characteristics of the computed potential energy surface. Next, the fitting of the computed potential energy points to a suitable analytical functional form is described, and maps of the fitted potential surface are displayed. The methodology and results of a classical trajectory calculation utilizing the fitted potential surface are presented. Finally, the significance of the trajectory study results is discussed, and generalizations concerning dynamical aspects of Ca + HF scattering are drawn.

Jaffe, R. L.↗

Accurate Energetics for the Unimolecular Decomposition of HN2

We report internally contracted configuration interaction (ICCI) calculations along the constrained energy minimum (CEM) path for the addition of H to N2 with both CC-pVQZ and augmented CC-pVQZ basis sets. The results agree to within 0.5 kcal/mol with the larger basis set results of Walch, which were used to calibrate the HN2 potential energy surface used by Koizumi, Schatz, and Walch in coupled channel calculations to determine the lifetimes for unimolecular decay of the lowest vibrational levels of HN2. These small changes in the potential energy surface should have no significant effect on the computed lifetimes.

Walch, Stephen P.↗

Ab Initio Calculation of Accurate Vibrational Frequencies for Molecules of Interest in Atmospheric Chemistry

Due to advances in quantum mechanical methods over the last few years, it is now possible to determine ab initio potential energy surfaces in which fundamental vibrational frequencies are accurate to within +/- 8 cm(sup -1) on average, and molecular bond distances are accurate to within +/- 0.001-0.003 A, depending on the nature of the bond. That is, the potential energy surfaces have not been scaled or empirically adjusted in any way, showing that theoretical methods have progressed to the point of being useful in analyzing spectra that are not from a tightly controlled laboratory environment, such as rovibrational spectra from the interstellar medium. Some recent examples demonstrating this accuracy win be presented and discussed. These include the HNO, CH4, C2H4, and ClCN molecules. The HNO molecule is interesting due to the very large H-N anharmonicity, while ClCN has a very large Fermi resonance. The ab initio studies for the CH4 and C2H4 molecules present the first accurate full quartic force fields of any kind (i.e., whether theoretical or empirical) for a five-atom and six-atom system, respectively.

Lee, Timothy J.↗

Ab initio study of C + H3+ reactions

The reaction C + H3+ --> CH(+) + H2 is frequently used in models of dense interstellar cloud chemistry with the assumption that it is fast, i.e. there are no potential energy barriers inhibiting it. Ab initio molecular orbital study of the triplet CH3+ potential energy surface (triplet because the reactant carbon atom is a ground state triplet) supports this hypothesis. The reaction product is 3 pi CH+; the reaction is to exothermic even though the product is not in its electronic ground state. No path has been found on the potential energy surface for C + H3+ --> CH2(+) + H reaction.

NASA Discipline Exobiology↗

The Role of Conical Intersections in Electronic Quenching of a State OH by H2 and N2

Lester and coworkers have experimentally characterized complexes of OH in the X and A states with H2 and N2. Recently, we have carried out ab initio calculations of relevant portions of the ground state and excited state potential energy surfaces for these systems, including the conical intersection regions, which are responsible for electronic quenching of the A state of OH by H2 and N2. Both of these systems have weakly bound complexes in the X state and strongly bound complexes in the A state. The OH-H2 complex is T-shaped, while the OH-N2 complex is collinear. In both cases the H end of OH is oriented toward the H2 or N2 molecule, respectively. Rotation of the OH about its center of mass involves only a small barrier and rotation by 1800 so that the 0 end of OH is oriented toward the H2 or N2 molecule leads to conical intersections with the ground state surface. Since there is about 95 kcal/mol of available energy after crossing to the ground state surface, chemical reactions on the ground state surface are possible, in addition to electronic quenching. In the case of OH-H2, the conical intersection is much lower in energy than the OH (A state) + H2 asymptotic energy and the A state complex can be characterized as a hanging well on the upper cone of the conical intersection. Passage through the conical intersection places the system on the ground state potential energy surface with the possibility of going to OH (X state) + H2 (i.e. electronic quenching) or to H2O + H (reaction). For OH-N2, the A state complex is also a hanging well on the upper cone of the conical intersection, but the conical intersection is only slightly below the OH (A state) + H2 asymptotic energy. Passage through the conical intersection can lead to OH (X state) + N2 (i.e. electronic quenching), but so far reactive pathways have not been found.

Walch, Stephen P.↗

Theoretical Studies of the HO/HO2 Catalytic Cycle for Ozone Destruction

Recently it has been determined that the HO/HO2 catalytic cycle accounts for nearly one-half of the total ozone depletion in the lower stratosphere. The catalytic cycle is: (1) HO + O3 yields HO2 + O2; (2) HO2 + O3 yields HO + O2 + O2. The net reaction is 2O3 yields 3O2. The rate limiting step in this process is the reaction of HO2 with ozone. There is a problem extending the experimental measurement of the rate of this reaction over the range 233-400 K down to stratospheric temperatures of 210-220 K. Therefore we have undertaken a project to determine the temperature dependence of the rate constant for this reaction in the low temperature region. The first step in this project, which is described in this poster, is the determination of the relevant potential energy surfaces. The calculations use CASSCF/derivative methods to define the pathways followed by CASSCF/ACPF to determine the energetics. The HO + O3 reaction is found to proceed through an HO4 complex, which is unstable with respect to HO2 + O2. The HO2 +O3 reaction is more complex. One pathway, which has been characterized, is the formation of an HO5 complex which decomposes to HO3 + O2 and subsequently to HO + O2 + O2. Another pathway, which is believed to also play a role, is hydrogen abstraction to give O2 + HO3 and subsequent decomposition of HO3 to HO + O2. Isotopic labeling experiments indicate that the later pathway is dominant. However, so far attempts to locate the saddle point for this pathway have not been successful. We have also characterized the potential energy surfaces for a number of species involved in these reactions, including HO3 and triplet O4. The triplet O4 species is probably involved in the reaction of vibrationally excited O2 with ground state O2 leading to O3 + O. The latter reaction is believed to be important as an additional source of stratospheric ozone.

Walch, Stephen P.↗

Reactive Resonances in N+N2 Exchange Reaction

Rich reactive resonances are found in a 3D quantum dynamics study of the N + N2 exchange reaction using a recently developed ab initio potential energy surface. This surface is characterized by a feature in the interaction region called Lake Eyring , that is, two symmetric transition states with a shallow minimum between them. An L2 analysis of the quasibound states associated with the shallow minimum confirms that the quasibound states associated with oscillations in all three degrees of freedom in Lake Eyring are responsible for the reactive resonances in the state-to-state reaction probabilities. The quasibound states, mostly the bending motions, give rise to strong reasonance peaks, whereas other motions contribute to the bumps and shoulders in the resonance structure. The initial state reaction probability further proves that the bending motions are the dominating factors of the reaction probability and have longer life times than the stretching motions. This is the first observation of reactive resonances from a "Lake Eyring" feature in a potential energy surface.

Wang, Dunyou↗

Half-projected Hartree-Fock calculations on several small molecules

The half-projected Hartree-Fock (HPHF) method is examined with respect to its ability to obtain molecular correlation, describe molecular potential energy surfaces, and provide a one-particle basis for more elaborate treatments. The equivalence, aside from questions of efficiency, of two different HPHF algorithms is demonstrated. The results of calculations on H2O, C2, N2, and CH2 indicate that the performance of the HPHF method, with spin projection in appropriate cases, is roughly equivalent to limited MCSCF treatments. In particular, a small but important fraction of the correlation energy, qualitatively correct potential energy surfaces, and good one-particle orbital bases are obtained.

Lengsfield, B. H., III↗