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 541 records · Page 30

Experimental and Ab Initio Studies of the HDO Absorption Spectrum in the 13165-13500 1/cm Spectral Region

The HDO absorption spectrum has been recorded in the 13165 - 13500 cm(exp-1) spectral region by Intracavity Laser Absorption Spectroscopy. The spectrum (615 lines), dominated by the 2n2 + 3n3 and n1+3n3 bands was assigned and modeled leading to the derivation of 196 accurate energy levels of the (103) and (023) vibrational states. Finally, 150 of these levels have been reproduced by an effective Hamiltonian involving two vibrational dark states interacting with the (023) and ( 103) bright states. The rms deviation achieved by variation of 28 parameters is 0.05-1 cm, compared to an averaged experimental uncertainty of 0.007-1 cm, indicating the limit of validity of the effective Hamiltonian approach for HDO at high vibrational excitation. The predictions of previous ab initio calculations of the HDO spectrum were extensively used in the assignment process. The particular spectral region under consideration has been used to test and discuss the improvements of new ab initio calculations recently performed on the basis of the same potential energy surface but with an improved dipole moment surface. The improvements concern both the energy levels and the line intensities. In particular, the strong hybrid character of the n1+3n3 band is very well accounted for by the the new ab initio calculations.

Schwenke, David↗

AB Initio Study of the Structure and Spectroscopic Properties of Halogenated Thioperoxy Radicals

Thioperoxy (XSO or XOS) radicals exist in a variety of chemical environments, and they have as a consequence drawn some interest. HSO, an important species in the chemistry of the troposphere, has been examined both experimentally. The halogenated (X = F, Cl or Br) peroxy species and isovalent thioperoxy species have been studied less, but they too are potentially interesting because oxidized sulfur species and halogen sources are present in the atmosphere. Learning the fate of XSO and XOS radicals is important to understanding the atmospheric oxidation chemistry of sulfur compounds. Of these, FSO and ClSO are particularly interesting because they have been directly detected spectroscopically. Recent studies in our laboratory on the photochemistry of thionyl halides (X2SO; where X = F or Cl) have suggested new ways to generate XSO species. The laser-induced photodissociation of thionyl fluoride, F2SO, at 193 nm and thionyl chloride, ClSO, at 248 nm is characterized by a radical mechanism, X2SO -> XSO + X. The structure of FSO has been characterized experimentally by Endo et cd. employing microwave spectroscopy. Using the unrestricted Hartree-Fock (UHF) self-consistent field (SCF) method, Sakai and Morokuma computed the electronic structure of the ground (sup 2)A" and the first excited (sup 2)A' states of FSO. Electron correlation was not taken into account in their study. In a laser photodissociation experiment, Huber et al. identified ClSO mass spectromctrically. ClSO has also been detected in low temperature matrices by EPR and in the gas phase by far IR laser magnetic resonance. Although the structure of FSO is known in detail, the only study, experimental or theoretical, of CISO has been an ab initio HFSCF study by Hinchliffe. Electron correlation corrections were also excluded from this study. In order to better understand the isomerization and dissociation dynamics of the radical species, we have performed ab initio correlated studies of the potential energy surfaces (PES) of ClSO and its isomer ClOS at the QCISD(T)/6-31 G* level of theory. For FSO and FOS, more extensive QCISD/6-31 1G(2df) calculations have been possible, and the results are summarized here.

Munoz, Luis A.↗

Ab Initio Vibrational Levels For HO2 and Vibrational Splittings for Hydrogen Atom Transfer

We calculate vibrational levels and wave functions for HO2 using the recently reported ab initio potential energy surface of Walch and Duchovic. There is intramolecular hydrogen atom transfer when the hydrogen atom tunnels through a T-shaped saddle point separating two equivalent equilibrium geometries, and correspondingly, the energy levels are split. We focus on vibrational levels and wave functions with significant splitting. The first three vibrational levels with splitting greater than 2/cm are (15 0), (0 7 1) and (0 8 0) where V(sub 2) is the O-O-H bend quantum number. We discuss the dynamics of hydrogen atom transfer; in particular, the O-O distances at which hydrogen atom transfer is most probable for these vibrational levels. The material of the proposed presentation was reviewed and the technical content will not reveal any information not already in the public domain and will not give any foreign industry or government a competitive advantage.

Barclay, V. J.↗

Computed Minimum Energy Pathway for Isomerization in Ketene

A minimum energy pathway for interchange of the CH groups in ketene via a C2v structure has been obtained using complete active space self consistent field (CASSCF) derivative methods with a polarized valence double zeta basis set to define the reaction pathway followed by multi-reference internally contracted configuration interaction (ICCI) calculations with a [3s3p2d/3s2p] basis set to determine the energetics. Qualitatively, the C2v structure is found to be a shallow minimum on the potential energy surface separated from ketene by a small barrier (0.2 kcal/mol), a second minimum, and a larger barrier (3.0 kcal/mol). The minimum energy pathway leading from the C2v minimum to ketene starts by simultaneous rotation of the farther CH group out of the plane and away from the oxygen followed by increase of the CCO angle and subsequent 1,2-migration of the H of the nearer CH group toward the carbon of the farther CH group.

Walch, Stephen P.↗

MP2 Studies of Relativistic Effects on the Linear Stationary Points of the H+Cl2 Yields HCl+Cl and Cl+HCl Yields ClH+Cl Reactions

The influence of relativistic effects on the linear stationary points of the potential energy surfaces of the ClHCl and HCl2 systems is studied. Scalar relativistic effects have little influence on the geometries and the energies of these points. Spin-orbit effects have no influence on the barrier of the forward reaction of H with Cl2, but increase the reaction energy. The total relativistic effect is an increase of the reaction energy by 2.3 kJ /mol. The barrier for the iso-energetic exchange reaction of HCl with Cl is decreased by 0.5 kJ /mol due to the scalar relativistic effects. This small decrease is cancelled by the larger increase of 2.5 kJ /mol due to spin-orbit effects, giving a total relativistic effect of 2.0 kJ /mol.

Visscher, L.↗

Towards the Synthesis and Identification of Tetrahedral N4

Tetrahedral (Td)-N4 has attracted attention as a high energy density material candidate. Although computational studies have shown Td-N4 to be a stable structure on the ground-state potential energy surface, it has yet to be experimentally observed. This study investigates possible synthetic pathways via excited state surfaces. Preliminary results include a comparison of vertical excitation energies from Td-N4 using CIS, CIS(D), CIS(3), LRCCSD, CASSCF, CASPT2, and MRCI ab initio methods. Previous studies have predicted weak IR vibrational intensities. Therefore, alternative detection schemes of Td-N4 are also investigated, such as UV and vibrational Raman spectroscopy.

Dateo, Christopher E.↗

Anomalous Symmetries Of The Rovibrational Levels of HO2: Consequences Of A Conical Intersection

We show that the geometric phase arising from a conical intersection of the lowest potential energy surfaces of HO2 causes its bending vibrational wavefunctions to be double valued, which allows them to be locally symmetric on one side of the intersection and locally antisymmetric on the other. The material of the proposed publication was reviewed and the technical content will not reveal any information not already in the public domain and will not give any foreign industry or government a competitive advantage.

Barclay, V. J.↗

Minimum Energy Pathways for Chemical Reactions

Computed potential energy surfaces are often required for computation of such parameters as rate constants as a function of temperature, product branching ratios, and other detailed properties. We have found that computation of the stationary points/reaction pathways using CASSCF/derivative methods, followed by use of the internally contracted CI method to obtain accurate energetics, gives useful results for a number of chemically important systems. The talk will focus on a number of applications to reactions leading to NOx and soot formation in hydrocarbon combustion.

Walch, S. P.↗

Spectrum Synthesis of Hot Water in Sunspots and Selected Cool Stars

Very recently, Partridge and Schwenke completed an elaborate theoretical computation of the potential energy surface and dipole moment function for H2O. They have used their results to predict the positions and strengths of nearly 308 million lines. This line tabulation is the most complete now available. It extends to sufficiently high excitations that the spectra of M-stars may be modelled with greater accuracy than ever before provided the predicted line parameters of Partridge and Schwenke are themselves accurate. We have computed synthetic sunspot spectra using the Partridge and Schwenke line list and the sunspot umbral models of Maltby et al. In this display, we compare these synthetic spectra with published high resolution sunspot atlases. We demonstrate the extent to which the new line list successfully predicts the sunspot spectrum and suggest where improvements are necessary. Using the new tabulation, we also illustrate the extent to which hot stellar blankets the H, K and L passbands for select K and M star model atmospheres.

Carbon, D. F.↗

Electron-Impact Ionization of Methane

We report a study of the total ionization of CH_4 by electron impact and its dissociative ionization from the ^2T_2 channel. The calculation of the total ionization cross section uses the improved Binary-Encounter-Dipole model (iBED).(W. A Huo, Phys. Rev. A (submitted for publication).) The dipole Born cross section in the model is expressed in terms of a three-term representation and the optical oscillator strengths are taken from Backx and Van der Wiel.(C. Backx and M. J. Van der Wiel, I Phys. B 18) 3020 (1975). The nuclear dynamics for the dissociation of the ^2T_2 channel is studied using the statistical model. A search of the potential energy surface of the ^2T_2 state of CH_4^+ shows two minima, of C_2v and C_3v symmetries, in agreement with earlier calculations. ((a) K. Takeshita, J. Chem. Phys. 86), 329 (1987). (b) R. F. Frey and E. R. Davidson, J. Chem. Phys. 88, 1775 (1988). The dissociation of the CH_4^+ to CH_3^+ + H goes through a saddle point. Comparison with recent experimental data will be presented and the role of Jahn-Teller effect discussed.

Huo, Winifred M.↗

Dissociative Ionization of Benzene by Electron Impact

We report a theoretical study of the dissociative ionization (DI) of benzene from the low-lying ionization channels. Our approach makes use of the fact that electron motion is much faster than nuclear motion and DI is treated as a two-step process. The first step is electron-impact ionization resulting in an ion with the same nuclear geometry as the neutral molecule. In the second step the nuclei relax from the initial geometry and undergo unimolecular dissociation. For the ionization process we use the improved binary-encounter dipole (iBED) model. For the unimolecular dissociation step, we study the steepest descent reaction path to the minimum of the ion potential energy surface. The path is used to analyze the probability of unimolecular dissociation and to determine the product distributions. Our analysis of the dissociation products and the thresholds of the productions are compared with the result dissociative photoionization measurements of Feng et al. The partial oscillator strengths from Feng et al. are then used in the iBED cross section calculations.

Huo, Winifred↗

Modeling of Switching and Hysteresis in Molecular Transport

The conventional way of modeling current transport in two and three terminal molecular devices could be inadequate for certain cases involving switching and hysteresis. Here we present an alternate approach. Contrary to the regular way where applied bias directly modulates the conducting energy levels of the molecule, our method introduces a nonlinear potential energy surface varying with the applied bias as a control parameter. A time-dynamics is also introduced properly accounting for switching and hysteresis behavior. Although the model is phenomenological at this stage, we believe any detailed model would contain similar descriptions at its core.

Samanta, Manoj P.↗

First Principle Predictions of Isotopic Shifts in H2O

We compute isotope independent first and second order corrections to the Born-Oppenheimer approximation for water and use them to predict isotopic shifts. For the diagonal correction, we use icMRCI wavefunctions and derivatives with respect to mass dependent, internal coordinates to generate the mass independent correction functions. For the non-adiabatic correction, we use scaled SCF/CIS wave functions and a generalization of the Handy method to obtain mass independent correction functions. We find that including the non-adiabatic correction gives significantly improved results compared to just including the diagonal correction when the Born-Oppenheimer potential energy surface is optimized for H2O-16. The agreement with experimental results for deuterium and tritium containing isotopes is nearly as good as our best empirical correction, however, the present correction is expected to be more reliable for higher, uncharacterized levels.

Schwenke, David W.↗

Anharmonic Vibrational Spectroscopy of the F-(H20)n, complexes, n=1,2

We report anharmonic vibrational spectra (fundamentals, first overtones) for the F-(H(sub 2)O) and F-(H(sub 2)O)2 clusters computed at the MP2 and CCSD(T) levels of theory with basis sets of triple zeta quality. Anharmonic corrections were estimated via the correlation-corrected vibrational self-consistent field (CC-VSCF) method. The CC-VSCF anharmonic spectra obtained on the potential energy surfaces evaluated at the CCSD(T) level of theory are the first ones reported at a correlated level beyond MP2. We have found that the average basis set effect (TZP vs. aug-cc-pVTZ) is on the order of 30-40 cm(exp -1), whereas the effects of different levels of electron correlation [MP2 vs. CCSD(T)] are smaller, 20-30 cm(exp -1). However, the basis set effect is much larger in the case of the H-bonded O-H stretch of the F-(H(sub 2)O) cluster amounting to 100 cm(exp -1) for the fundamentals and 200 cm (exp -1) for the first overtones. Our calculations are in agreement with the limited available set of experimental data for the F-(H(sub 2)O) and F-(H(sub 2)O)2 systems and provide additional information that can guide further experimental studies.

Chaban, Galina M.↗

Ab Initio Calculations of Water Line Strengths

We report on the determination of a high quality ab initiu potential energy surface (PES) and dipole moment function for water. This PES is empirically adjusted to improve the agreement between the computed line positions and those from the HITRAN 92 data base with J less than 6 for H2O. The changes in the PES are small, nonetheless including an estimate of core (oxygen 1s) electron correlation greatly improves the agreement with experiment. Using this adjusted PES, we can match 30,092 of the 30,117 transitions in the HITRAN 96 data base for H2O with theoretical lines. The 10,25,50,75, and 90 percentiles of the difference between the calculated and tabulated line positions are -0.11, -0.04, -0.01, 0.02, and 0.07 l/cm. Non-adiabatic effects are not explicitly included. About 3% of the tabulated line positions appear to be incorrect. Similar agreement using this adjusted PES is obtained for the oxygen 17 and oxygen 18 isotopes. For HDO, the agreement is not as good, with root-mean-square error of 0.25 l/cm for lines with J less than 6. This error is reduced to 0.02 l/cm by including a small asymmetric correction to the PES, which is parameterized by simultaneously fitting to HDO md D2O data. Scaling this correction by mass factors yields good results for T2O and HTO. The intensities summed over vibrational bands are usually in good agreement between the calculations and the tabulated results, but individual lines strengths can differ greatly. A high temperature list consisting of 307,721,352 lines is generated for H2O using our PES and dipole moment function.

Schwenke, David W.↗

Quantum Scattering Study of Ro-Vibrational Excitations in N+N(sub 2) Collisions under Re-entry Conditions

A three-dimensional time-dependent quantum dynamics approach using a recently developed ab initio potential energy surface is applied to study ro-vibrational excitation in N+N2 exchange scattering for collision energies in the range 2.1- 3.2 eV. State-to-state integral exchange cross sections are examined to determine the distribution of excited rotational states of N(sub 2). The results demonstrate that highly-excited rotational states are produced by exchange scattering and furthermore, that the maximum value of (Delta)j increases rapidly with increasing collision energies. Integral exchange cross sections and exchange rate constants for excitation to the lower (upsilon = 0-3) vibrational energy levels are presented as a function of the collision energy. Excited-vibrational-state distributions for temperatures at 2,000 K and 10,000 K are included.

Wang, Dunyou↗

The C4H7+ cation. A theoretical investigation

The potential energy surface of the C4H7+ cation has been investigated with ab initio quantum chemical theory. Extended basis set calculations, including electronic correlation, show that cyclobutyl and cyclopropylcarbinyl cation are equally stable isomers. The saddle point connecting these isomers lies 0.6 kcal/mol above the minima. The global C4H7+ minimum corresponds to the 1-methylallyl cation, which is 9.0 kcal/mol more stable than the cyclobutyl and the cyclopropylcarbinyl cation and 9.5 kcal/mol below the 2-methylallyl cation. These results are in excellent agreement with experimental data.

NASA Discipline Number 52-10↗

Vibrational energy levels for CH4 from an ab initio potential

Many areas of astronomy and astrophysics require an accurate high temperature spectrum of methane (CH4). The goal of the present research is to determine an accurate ab initio potential energy surface (PES) for CH4. As a first step towards this goal, we have determined a PES including up to octic terms. We compare our results with experiment and to a PES based on a quartic expansion. Our octic PES gives good agreement with experiment for all levels, while the quartic PES only for the lower levels.

Methane/chemistry↗