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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Computational Study of Field Initiated Surface Reactions for Synthesis of Diamond and Silicon

This project involves using quantum chemistry to simulate surface chemical reactions in the presence of an electric field for nanofabrication of diamond and silicon. A field delivered by a scanning tunneling microscope (STM) to a nanometer scale region of a surface affects chemical reaction potential energy surfaces (PES) to direct atomic scale surface modification to fabricate sub-nanometer structures. Our original hypothesis is that the applied voltage polarizes the charge distribution of the valence electrons and that these distorted molecular orbitals can be manipulated with the STM so as to change the relative stabilities of the electronic configurations over the reaction coordinates and thus the topology of the PES and reaction kinetics. Our objective is to investigate the effect of applied bias on surface reactions and the extent to which STM delivered fields can be used to direct surface chemical reactions on an atomic scale on diamond and silicon. To analyze the fundamentals of field induced chemistry and to investigate the application of this technique for the fabrication of nanostructures, we have employed methods capable of accurately describing molecular electronic structure. The methods we employ are density functional theory (DFT) quantum chemical (QC) methods. To determine the effect of applied bias on surface reactions we have calculated the QC PESs in various applied external fields for various reaction steps for depositing or etching diamond and silicon. We have chosen reactions which are thought to play a role in etching and the chemical vapor deposition growth of Si and diamond. The PESs of the elementary reaction steps involved are then calculated under the applied fields, which we vary in magnitude and configuration. We pay special attention to the change in the reaction barriers, and transition state locations, and search for low energy reaction channels which were inaccessible without the applied bias.

Musgrave, Charles Bruce↗
NASA NTRS2017

On the Development of a New Nonequilibrium Chemistry Model for Mars Entry

This paper represents a summary of results to date of an on-going effort at NASA Ames Research Center to develop a physics-based non-equilibrium model for hypersonic entry into the Martian atmosphere. Our approach is to first compute potential energy surfaces based on accurate solutions of the electronic Schroedinger equation and then use quasiclassical trajectory calculations to obtain reaction cross sections and rate coefficients based on these potentials. We have presented new rate coefficients for N2 dissociation and CO dissociation and exchange reactions. These results illustrate shortcomings with some of the rate coefficients in Parks original T-Tv model for Mars entries and with some of the 30-45 year old shock tube data. We observe that the shock tube experiments of CO + O dissociation did not adequately account for the exchange reaction that leads to formation of C + O2. This reaction is actually the primary channel for CO removal in the shock layer at temperatures below 10,000 K, because the reaction enthalpy for exchange is considerably lower than the comparable value for dissociation.

quasiclassical trajectory calculation↗
NASA NTRS2016

Computing Highly Accurate Spectroscopic Line Lists for Characterization of Planetary Atmospheres: CO2 and SO2 Line Lists Needed for Modeling Venus

Over the last decade, it has become apparent that the most effective approach for determining highly accurate rotational and rovibrational line lists for molecules of interest in planetary atmospheres and other astrophysical environments is through a combination of highresolution laboratory experiments coupled with state-of-the art ab initio quantum chemistry methods. The approach involves computing the most accurate potential energy surface (PES) possible using state-of-the art electronic structure methods, followed by computing rotational and rovibrational energy levels using an exact variational method to solve the nuclear Schrödinger equation. Then, reliable experimental data from high-resolution experiments is used to refine the ab initio PES in order to improve the accuracy of the computed energy levels and transition energies. From the refinement step, we have been able to achieve an accuracy of approximately 0.015 cm-1 for rovibrational transition energies, and even better for purely rotational transitions. This combined "experiment / theory" approach allows for determination of essentially a complete line list, with hundreds of millions of transitions, and having the transition energies and intensities be highly accurate. Our group has successfully applied this approach to determine highly accurate line lists for NH3, CO2 and isotopologues, and SO2 and isotopologues. Here I will report our latest results for CO2 and SO2 including all isotopologues. Comparisons to the available data in HITRAN2012 and other available databases will be shown, though we note that our line lists for SO2 are significantly more complete than any other databases. Since it is important to span a large temperature range in order to model the spectral signature of Venus as well as exoplanets, we will demonstrate how the spectra change on going from low temperatures (100 K) to higher temperatures (500 K to 1500 K).

potential energy surface (PES)5405 PLANETARY SCIEN↗
NASA NTRS2024

Accurate Infrared Line Lists for 20 Isotopologues of CS 2 at Room Temperature

To facilitate atmospheric and spectroscopic studies of carbon disulfide, or CS 2 , in both planetary and exoplanetary atmospheres, we adopt the “Best Theory + Reliable High-resolution Experiment” algorithm to generate semiempirical IR line lists for the 20 most abundant CS 2 isotopologues, denoted as Ames-296K. The IR lists are computed using the Ames-1 potential energy surface, refined using the experimental transition set and an ab initio dipole moment surface fitted from CCSD(T)/aug-cc-pV(T/Q/5+d)Z dipoles extrapolated to a one-particle basis set limit. The IR lists cover the range of 0–10,000 cm -1 , with an S 296K cutoff at 10 -31 cm -1 /molecule·cm -2 (abundance included). A “natural” IR line list at 296 K includes about 10 million lines of the 20 isotopologues, with their intensities scaled by the corresponding abundances. The zero-point energy, partition functions, and abundances are reported for each isotopologue. The energy levels in the global effective Hamiltonian model for 12 C 32 S 2 are adopted to improve the line position accuracy. This new IR list for the main isotopologue is denoted as A+I.296K. Reliable HITRAN2020 line positions are also utilized to improve the accuracy of the 32 S 12 C 34 S, 32 S 12 C 33 S, and 32 S 13 C 32 S isotopologue line lists. The final composite line list is validated against Pacific Northwest National Laboratory experimental cross sections, showing excellent agreement. The agreement supports the quality of the composite line list and the power of synergy between experiment and theory. The new data are proposed for use in updating and expanding the CS 2 data in HITRAN and other high-resolution IR databases. Supplementary files are available in Zenodo and AHED.

Exoplanet atmospheric composition↗
NASA NTRS2023

AI-3000K Infrared Line List for Hot CO 2

AI-3000K is a semi-empirical IR line list constructed for hot CO 2 spectra analysis and simulation up to 3000 – 4000 K. Compared to previously published Ames-1000K and UCL-4000, it represents a major upgrade, utilizing a new algorithm for optimization and including the latest improvements in potential energy surface (PES), dipole moment surface (DMS), and room temperature IR line list (Ames-2021 296K). To maximize the success of introducing experimental based accurate line positions, a new PES (X01d) was refined with respect to >800 selected CDSD2019 [Tashkun et al. JQSRT (2019) 228, 124] energy levels in the range of 0 – 24,000 cm -1 , with σ rms = 0.5 – 0.7 cm -1 . Most differences between the X01d PES based levels and CDSD2019 energies are within ± 2 cm -1 . A new DMS is fitted from extrapolated CCSD(T)/aug-cc-pV(T,Q,5)Z dipole calculations, with σ rms = 5.1 × 10 -6 au for 11,155 geometries up to 40,000 cm -1 , denoted Ames-2021–40K. Compared to the best available Ames-2021 DMS and room temperature IR line list [Huang et al J. Phys. Chem. A (2022) 126, 5940], the relative intensity differences are expected to be ∼ 1 ‰. The line position accuracy of ''X01d + Ames-2021–40K'' IR line list is significantly improved by adopting CDSD2019 energy levels up to 24,000 cm -1 (J ≤ 150). The Einstein A21 coefficients for E' < 15,000 cm -1 transitions are replaced by more accurate values from the Ames-2021 296K IR line list. In short, the AI-3000K is the X01d PES and Ames-2021–40K DMS based line list enhanced with the A 21 of Ames-2021 296K line list and CDSD energy levels. It provides continuous coverage from 0 to 20,000 cm -1 for the four most abundant isotopologues: 12 C 16 O 2 (626), 13 C 16 O 2 (636), 16 O 12 C 18 O (628), and 16 O 12 C 17 O (627). The impacts of isotopologue and E'/E'' cutoffs have been examined. Intensity convergence (not accuracy) of AI-3000K line list is quantitatively estimated in 1 cm -1 bins. It is better than 99% in the whole range of 0 – 20,000 (or 10,000) cm -1 at 1000 K (or 2000 K), or better than 90% in the whole range of 0 – 15,000 (or 9000) cm -1 at 2000 K (or 3000 K), respectively. Convergence beyond 3000 K will require a new PES and DMS for E' > 40,000 cm -1 . The AI-3000K and HITEMP based IR simulations are compared to high resolution shock tube experiments for CO 2 -Ar mixture up to 2000 K. With line position accuracy comparable to that of HITEMP, AI-3000K IR line list yields better agreements at more frequencies. Potential sources of discrepancies with experiment are discussed.

AI-3000K↗
NASA NTRS1983

A theoretical study on the mechanism of electronic to vibrational energy transfer in Hg/3P/ + CO

The mechanism of electronic-to-vibrational (E-V) energy transfer in Hg(3P) + CO collisions has been studied theoretically. The configuration interaction (CI) method was employed to calculate potential energy surfaces of the collision system. A simplified theoretical model, based on the reaction coordinate concept and the calculated potential energy characteristics, was used to discuss the mechanism of the singlet-triplet transition and the energy disposal in the collision. The results obtained were that: (a) the quenching process processed via a collision complex mechanism; and that (b) the triplet-singlet transition occurs near the collinear geometry. A model classical trajectory calculation gives a product CO vibrational distribution in good agreement with the experimental result.

Kato, S.↗
NASA NTRS1992

Equilibrium properties of transition-metal ion-argon clusters via simulated annealing

The geometrical structures of M(+) (Ar)n ions, with n = 1-14, have been studied by the minimization of a many-body potential surface with a simulated annealing procedure. The minimization method is justified for finite systems through the use of an information theory approach. It is carried out for eight potential-energy surfaces constructed with two- and three-body terms parametrized from experimental data and ab initio results. The potentials should be representative of clusters of argon atoms with first-row transition-metal monocations of varying size. The calculated geometries for M(+) = Co(+) and V(+) possess radial shells with small (ca. 4-8) first-shell coordination number. The inclusion of an ion-induced-dipole-ion-induced-dipole interaction between argon atoms raises the energy and generally lowers the symmetry of the cluster by promoting incomplete shell closure. Rotational constants as well as electric dipole and quadrupole moments are quoted for the Co(+) (Ar)n and V(+) (Ar)n predicted structures.

Asher, Robert L.↗
NASA NTRS1978

Calculated rate constants for the reaction ClO + O yields Cl + O2 between 220 and 1000 deg K

Classical trajectory calculations are presented for the reaction ClO + O yields Cl + O2, a reaction which is an important step in the chlorine-catalyzed destruction of ozone which is thought to occur in the 220 and 1000 K. The calculated rate constant is 4.36 x 10 to the minus 11th power exp (-191/T)cu cm molecule (-1)s(-1) and its value at 300 K is 2.3 plus or minus 10 to the 11th power cu cm molecule (-1)s(-1), about a factor of 2 lower than recent experimental data. The empirical potential energy surface used in the calculations was constructed to fit experimental data for ClO, O2 and ClOO molecules. Other important features of this potential surface, such as the barrier to reaction, were varied systematically and calculations were performed for a range of conditions to determine the best theoretical rate constants. Results demonstrate the utility of classical trajectory methods for determining activation energies and other kinetic data for important atmospheric reactions.

Jaffee, R. L.↗
NASA NTRS2001

Effective Potential Energies and Transport Cross Sections for Atom-Molecule Interactions of Nitrogen and Nitrogen

The potential energy surfaces for H2-N and N2-N interactions are calculated by accurate ab initio methods and applied to determine transport data. The results confirm that an effective potential energy for accurately determining transport properties can be calculated using a single orientation. A simple method is developed to determine the dispersion coefficients of effective potential energies Effective potential energies required for O2-O collisions are determ=ined. The H2-N, N2-N, O2-H, and O2-O collision integrals are calculated and tabulated for a large range of temperatures. The theoretical values of the N2-N and O2-O diffusion coefficients compare well with measured data available at room temperature.

Stallcop, James R.↗
NASA NTRS1993

Calculations on the competition between association and reaction for C3H(+) + H2

The ab initio results presently obtained for the potential energy surface of competing associative and reactive ion-molecule processes of the reactants H2 and C3H(+) show that the linear ion is able to directly access the deep potential well of the propargyl ion H2CCCH(+). Using the minimum energy potential pathway and properties of the stationary point structures determined via the ab initio methods, calculation results are obtained for both the association-rate coefficient for C3H3(+) production as a function of density, and the branching ratio between the propargyl and cyclic structures of the ion.

Maluendes, Sergio A.↗
NASA NTRS1998

Energetics and Dynamics of GaAs Epitaxial Growth via Quantum Wave Packet Studies

The dynamics of As(sub 2) molecule incorporation into the flat Ga-terminated GaAs(100) surface is studied computationally. The time-dependent Schrodinger equation is solved on a two-dimensional potential energy surface obtained using density functional theory calculations. The probabilities of trapping and subsequent dissociation of the molecular As(sub 2) bond are calculated as a function of beam translational energy and vibrational quantum number of As(sub 2).

Dzegilenko, Fedor N.↗
NASA NTRS2001

Ab Initio Calculations of Anharmonic Vibrational Spectroscopy for Hydrogen Fluoride (HF)n (n=3,4) and Mixed Hydrogen Fluoride/Water (HF)n(H20)n (n=1,2,4) Clusters

Anharmonic vibrational frequencies and intensities are computed for hydrogen fluoride clusters (HF)n with n=3,4 and mixed clusters of hydrogen fluoride with water (HF)n(H2O)n where n=1,2. For the (HF)4(H2O)4 complex, the vibrational spectra are calculated at the harmonic level, and anharmonic effects are estimated. Potential energy surfaces for these systems are obtained at the MP2/TZP level of electronic structure theory. Vibrational states are calculated from the potential surface points using the correlation-corrected vibrational self-consistent field (CC-VSCF) method. The method accounts for the anharmonicities and couplings between all vibrational modes and provides fairly accurate anharmonic vibrational spectra that can be directly compared with experimental results without a need for empirical scaling. For (HF)n, good agreement is found with experimental data. This agreement shows that the MP2 potential surfaces for these systems are reasonably reliable. The accuracy is best for the stiff intramolecular modes, which indicates the validity of MP2 in describing coupling between intramolecular and intermolecular degrees of freedom. For (HF)n(H2O)n experimental results are unavailable. The computed intramolecular frequencies show a strong dependence on cluster size. Intensity features are predicted for future experiments.

Chaban, Galina M.↗
NASA NTRS2002

Ab initio calculations of anharmonic vibrational spectroscopy for hydrogen fluoride (HF)n (n = 3, 4) and mixed hydrogen fluoride/water (HF)n(H2O)n (n = 1, 2, 4) clusters

Anharmonic vibrational frequencies and intensities are computed for hydrogen fluoride clusters (HF)n, with n = 3, 4 and mixed clusters of hydrogen fluoride with water (HF)n(H2O)n where n = 1, 2. For the (HF)4(H2O)4 complex, the vibrational spectra are calculated at the harmonic level, and anharmonic effects are estimated. Potential energy surfaces for these systems are obtained at the MP2/TZP level of electronic structure theory. Vibrational states are calculated from the potential surface points using the correlation-corrected vibrational self-consistent field method. The method accounts for the anharmonicities and couplings between all vibrational modes and provides fairly accurate anharmonic vibrational spectra that can be directly compared with experimental results without a need for empirical scaling. For (HF)n, good agreement is found with experimental data. This agreement shows that the Moller-Plesset (MP2) potential surfaces for these systems are reasonably reliable. The accuracy is best for the stiff intramolecular modes, which indicates the validity of MP2 in describing coupling between intramolecular and intermolecular degrees of freedom. For (HF)n(H2O)n experimental results are unavailable. The computed intramolecular frequencies show a strong dependence on cluster size. Intensity features are predicted for future experiments.

Water/chemistry↗
NASA NTRS1982

Raindrop oscillations

A model of the change in shape of a raindrop is presented. Raindrops measured by two orthogonal cameras were classified by shape and orientation to determine the nature of the oscillation. A physical model based on potential energy was then developed to study the amplitude variation of oscillating drops. The model results show that oscillations occur about the equilibrium axis ratio, but the time average axis ratio if significantly more spherical for large amplitudes because of asymmetry in the surface potential energy. A generalization of the model to oscillations produced by turbulence yields average axis ratios that are consistent with the camera measurements. The model results for average axis ratios were applied to rainfall studies with a dual polarized radar.

Beard, K. V.↗
NASA NTRS2001

Ab Initio and Improved Empirical Potentials for the Calculation of the Anharmonic Vibrational States and Intramolecular Mode Coupling of N-Methylacetamide

The second-order Moller-Plesset ab initio electronic structure method is used to compute points for the anharmonic mode-coupled potential energy surface of N-methylacetamide (NMA) in the trans(sub ct) configuration, including all degrees of freedom. The vibrational states and the spectroscopy are directly computed from this potential surface using the Correlation Corrected Vibrational Self-Consistent Field (CC-VSCF) method. The results are compared with CC-VSCF calculations using both the standard and improved empirical Amber-like force fields and available low temperature experimental matrix data. Analysis of our calculated spectroscopic results show that: (1) The excellent agreement between the ab initio CC-VSCF calculated frequencies and the experimental data suggest that the computed anharmonic potentials for N-methylacetamide are of a very high quality; (2) For most transitions, the vibrational frequencies obtained from the ab initio CC-VSCF method are superior to those obtained using the empirical CC-VSCF methods, when compared with experimental data. However, the improved empirical force field yields better agreement with the experimental frequencies as compared with a standard AMBER-type force field; (3) The empirical force field in particular overestimates anharmonic couplings for the amide-2 mode, the methyl asymmetric bending modes, the out-of-plane methyl bending modes, and the methyl distortions; (4) Disagreement between the ab initio and empirical anharmonic couplings is greater than the disagreement between the frequencies, and thus the anharmonic part of the empirical potential seems to be less accurate than the harmonic contribution;and (5) Both the empirical and ab initio CC-VSCF calculations predict a negligible anharmonic coupling between the amide-1 and other internal modes. The implication of this is that the intramolecular energy flow between the amide-1 and the other internal modes may be smaller than anticipated. These results may have important implications for the anharmonic force fields of peptides, for which N-methylacetamide is a model.

Gregurick, Susan K.↗
NASA NTRS1974

Rotational excitation of HCN by collisions

Rate constants for the rotational excitation of HCN by collisions with He atoms at temperatures below 100 K were computed from first principles and are presented in tabular form. The potential energy surface was obtained by using the uniform electron gas model of Gordon and Kim (1972) and then joined smoothly to the asymptotic long-range perturbation theory potential valid at large separations. Quantum close-coupling theory was used to analyze the collision dynamics. Individual rates are believed to be accurate to within 50% above 30 K and within a factor of two below 20 K. The results should be extendable to excitation by collision with H2 and may therefore be of value in the study of interstellar clouds.

Green, S.↗