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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Photoionization in Ultraviolet Processing of Astrophysical Ice Analogs at Cryogenic Temperatures

Two recent experimental studies have demonstrated that amino acids or amino acid precursors are generated when astrophysical ice analogs are subjected to ultraviolet (UV) irradiation at cryogenic temperatures. Understanding the complete phenomenology of photoprocessing is critical to elucidating chemical reaction mechanisms that can function within an ice matrix under very cold conditions. Pushing beyond the much better characterized study of photolytic dissociation of chemical bonds through electronic excitation, this work explored the ability of UV radiation present in the interstellar medium to ionize small molecules embedded in ices. Quantum chemical calculations, including bulk solvation effects, were used to study the ionization of hydrogen (H2), water, and methanol (CH3OH) bound in small clusters of water. Ionization potentials were found to be much smaller in the condensed phase than in the gas phase; even a small cluster can account for large changes in the ionization potentials in ice, as well as the known formation of an OH--H3O+ pair in the case of H2O photoionization. To gauge the impact of photoionization on subsequent grain chemistry, the reaction between OH and CO in the presence of H3O+ was studied and compared with the potential energy surface without hydronium present, which is relevant to chemistry following photolysis. The differences indicate that the reaction is somewhat more likely to proceed to products (H + CO2) in the case of photoionization.

Woon, David E.↗
NASA NTRS2007

Pressure and Temperature Dependence of the Reaction of Vinyl Radical with Ethylene

This work reports measurements of absolute rate coefficients and Rice-Ramsperger-Kassel-Marcus (RRKM) master equation simulations of the C2H3 + C2H4 reaction. Direct kinetic studies were performed over a temperature range of 300-700 K and pressures of 20 and 133 mbar. Vinyl radicals (H2C=CH) were generated by laser photolysis of vinyl iodide (C2H31) at 266 nm, and time-resolved absorption spectroscopy was used to probe vinyl radicals through absorption at 423.2 nm. Measurements at 20 mbar are in good agreement with previous determinations at higher temperature. A weighted three-parameter Arrhenius fit to the experimental rate constant at 133 mbar, with the temperature exponent fixed, gives k = (7 +/- 1) x 10(exp -l4) cu cm/molecule/s (T/298 K)(exp 2) exp[-(1430 +/- 70) K/T]. RRKM master equation simulations, based on G3 calculations of stationary points on the C4H7 potential energy surface, were carried out to predict rate coefficients and product branching fractions. The predicted branching to 1-methylallyl product is relatively small under the conditions of the present experiments but increases as the pressure is lowered. Analysis of end products of 248 nm photolysis of vinyl iodide/ethylene mixtures at total pressures between 27 and 933 mbar provides no direct evidence for participation of I -methylallyl.

Ismail, Huzeifa↗
NASA NTRS2013

On the Use of Quartic Force Fields in Variational Calculations

The use of quartic force fields (QFFs) has been shown to be one of the most effective ways to efficiently compute vibrational frequencies for small molecules. In this paper we outline and discuss how the simple-internal or bond-length bond-angle (BLBA) coordinates can be transformed into Morse-cosine(-sine) coordinates which produce potential energy surfaces from QFFs that possess proper limiting behavior and can effectively describe the vibrational (or rovibrational) energy levels of an arbitrary molecular system. We investigate parameter scaling in the Morse coordinate, symmetry considerations, and examples of transformed QFFs making use of the MULTIMODE, TROVE, and VTET variational vibrational methods. Cases are referenced where variational computations coupled with transformed QFFs produce accuracies compared to experiment for fundamental frequencies on the order of 5 cm(exp −1) and often as good as 1 cm(exp −1).

quartic force field↗
NASA NTRS2014

An Experimental and Theoretical Study of Nitrogen-Broadened Acetylene Lines

We present experimental nitrogen-broadening coefficients derived from Voigt profiles of isotropic Raman Q-lines measured in the 2 band of acetylene (C2H2) at 150 K and 298 K, and compare them to theoretical values obtained through calculations that were carried out specifically for this work. Namely, full classical calculations based on Gordon's approach, two kinds of semi-classical calculations based on Robert Bonamy method as well as full quantum dynamical calculations were performed. All the computations employed exactly the same ab initio potential energy surface for the C2H2N2 system which is, to our knowledge, the most realistic, accurate and up-to-date one. The resulting calculated collisional half-widths are in good agreement with the experimental ones only for the full classical and quantum dynamical methods. In addition, we have performed similar calculations for IR absorption lines and compared the results to bibliographic values. Results obtained with the full classical method are again in good agreement with the available room temperature experimental data. The quantum dynamical close-coupling calculations are too time consuming to provide a complete set of values and therefore have been performed only for the R(0) line of C2H2. The broadening coefficient obtained for this line at 173 K and 297 K also compares quite well with the available experimental data. The traditional Robert Bonamy semi-classical formalism, however, strongly overestimates the values of half-width for both Qand R-lines. The refined semi-classical Robert Bonamy method, first proposed for the calculations of pressure broadening coefficients of isotropic Raman lines, is also used for IR lines. By using this improved model that takes into account effects from line coupling, the calculated semi-classical widths are significantly reduced and closer to the measured ones.

coefficients↗
NASA NTRS2013

Anharmonic Rovibrational Calculations of Singlet Cyclic C4 Using a New Ab Initio Potential and a Quartic Force

We report a CCSD(T)/cc-pCV5Z quartic force field (QFF) and a semi-global CCSD(T)-F12b/aug-cc-pVTZ potential energy surface (PES) for singlet, cyclic C4. Vibrational fundamentals, combinations and overtones are obtained using vibrational second-order perturbation theory (VPT2) and the vibrational configurationinteraction (VCI) approach. Agreement is within 10 cm(exp −1) between the VCI calculated fundamentals on the QFF and PES using the MULTIMODE (MM) program, and VPT2 and VCI results agree for the fundamentals. The agreement between VPT2- QFF and MM-QFF results is also good for the C4 combinations and overtones. The J = 1 and J = 2 rovibrational energies are reported from both VCI (MM) on the PES and VPT2 on the QFF calculations. The spectroscopic constants of (12)C4 and two C(sub 2v)-symmetry, single (13)C-substituted isotopologues are presented, which may help identification of cyclic C4 in future experimental analyses or astronomical observations.

quartic force field↗
NASA NTRS2015

Line Mixing in Parallel and Perpendicular Bands of CO2: A Further Test of the Refined Robert-Bonamy Formalism

Starting from the refined Robert-Bonamy formalism [Q. Ma, C. Boulet, and R. H. Tipping, J. Chem. Phys. 139, 034305 (2013)], we propose here an extension of line mixing studies to infrared absorptions of linear polyatomic molecules having stretching and bending modes. The present formalism does not neglect the internal degrees of freedom of the perturbing molecules, contrary to the energy corrected sudden (ECS) modeling, and enables one to calculate the whole relaxation matrix starting from the potential energy surface. Meanwhile, similar to the ECS modeling, the present formalism properly accounts for roles played by all the internal angular momenta in the coupling process, including the vibrational angular momentum. The formalism has been applied to the important case of CO2 broadened by N2. Applications to two kinds of vibrational bands (sigma yields sigma and sigma yields pi) have shown that the present results are in good agreement with both experimental data and results derived from the ECS model.

Semiclassical theories↗
NASA NTRS2016

Ames S-32 O-16 O-18 Line List for High-Resolution Experimental IR Analysis

By comparing to the most recent experimental data and spectra of the SO2 628 ν1/ν3 bands (see Ulenikov et al., JQSRT 168 (2016) 29-39), this study illustrates the reliability and accuracy of the Ames-296K SO2 line list, which is accurate enough to facilitate such high-resolution spectroscopic analysis. The SO2 628 IR line list is computed on a recently improved potential energy surface (PES) refinement, denoted Ames-Pre2, and the published purely ab initio CCSD(T)/aug-cc-pVQZ dipole moment surface. Progress has been made in both energy level convergence and rovibrational quantum number assignments agreeing with laboratory analysis models. The accuracy of the computed 628 energy levels and line list is similar to what has been achieved and reported for SO2 626 and 646, i.e. 0.01-0.03 cm(exp −1) for bands up to 5500 cm(exp −1). During the comparison, we found some discrepancies in addition to overall good agreements. The three-IR-list based feature-by-feature analysis in a 0.25 cm(exp −1) spectral window clearly demonstrates the power of the current Ames line lists with new assignments, correction of some errors, and intensity contributions from varied sources including other isotopologues. We are inclined to attribute part of detected discrepancies to an incomplete experimental analysis and missing intensity in the model. With complete line position, intensity, and rovibrational quantum numbers determined at 296 K, spectroscopic analysis is significantly facilitated especially for a spectral range exhibiting such an unusually high density of lines. The computed 628 rovibrational levels and line list are accurate enough to provide alternatives for the missing bands or suspicious assignments, as well as helpful to identify these isotopologues in various celestial environments. The next step will be to revisit the SO2 828 and 646 spectral analyses.

Sulfur dioxide↗
NASA NTRS2015

Energy and Spectroscopic Characterization of the Isomers of C4H3-, C6H3-, and C6H5-

Organic and inorganic molecules, neutral and ions have been observed in the interstellar medium. A few anions of organic molecules have also been observed recently. The Cassini spacecraft in the upper atmosphere of Titan has observed anions of large organic molecules. In this project we have studied the physical and spectroscopic properties of C4H3-, C6H3-, and C6H5-. We have optimized the geometrical structures of all low-lying isomers of the anions, calculated rotational, and harmonic vibrational frequencies of the anions mentioned above using the B3LYP density functional along with the augmented correlation consistent polar valence triple zeta (aug-cc-pVTZ) basis set. We have found many low-lying isomers on the potential energy surface of these anions.

Energy↗
NASA NTRS2018

Validating Advanced Thermophysics Models

Aerospace America year in review article for Thermophysics covering the COMARS data from ExoMars, meteor testing in arcj ets and potential energy surface calculations.

PES↗
NASA NTRS2017

Dissociation and Internal Excitation of Molecular Nitrogen Due to N + N2 Collisions Using Direct Molecular Simulation

In this work we present a molecular level study of N2+N collisions, focusing on excitation of internal energy modes and non-equilibrium dissociation. The computation technique used here is the direct molecular simulation (DMS) method and the molecular interactions have been modeled using an ab−initio potential energy surface (PES) developed at NASA's Ames Research Center. We carried out vibrational excitation calculations between 5000K and 30000K and found that the characteristic vibrational excitation time for the N + N2 process was an order of magnitude lower than that predicted by the Millikan and White correlation. It is observed that during vibrational excitation the high energy tail of the vibrational energy distribution gets over populated first and the lower energy levels get populated as the system evolves. It is found that the non-equilibrium dissociation rate coefficients for the N + N2 process are larger than those for the N2 + N2 process. This is attributed to the non-equilibrium vibrational energy distributions for the N + N2 process being less depleted than that for the N2 +N2 process. For an isothermal simulation we find that the probability of dissociation goes as 1/T(sub tr) for molecules with internal energy (epsilon(sub int)) less than approximately 9.9eV, while for molecules with epsilon (sub int) greater than 9.9eV the dissociation probability was weakly dependent on translational temperature of the system. We compared non-equilibrium dissociation rate coefficients and characteristic vibrational excitation times obtained by using the ab-initio PES developed at NASA's Ames Research Center to those obtained by using an ab-initio PES developed at the University of Minnesota. Good agreement was found between the macroscopic properties and molecular level description of the system obtained by using the two PESs.

Grover, Maninder S.↗
NASA NTRS2018

New O+CO2 Scattering Cross Sections and Their Impact on Photochemical Atmospheric Escape of Oxygen from Mars

Dissociative recombination of O2+ and CO2+ molecular ions in the upper atmosphere of Mars is a major source of suprathermal (hot) oxygen atoms. A significant fraction of hot O atoms are sufficiently fast to overcome Mars' gravitational potential and escape to space. This mechanism is known as photochemical escape (PE) and found to be one of the major escape mechanisms presently active on Mars. The Mars Atmosphere and Volatile Evolution Mission (MAVEN) can constrain the PE rate from measurements of electron temperature and density and ion temperature and density. Escape probabilities of hot O can then be constructed from in-situ neutral densities and escaping O flux loss due to collisions with atmospheric background gases. Due to CO2 abundance in the upper atmosphere, O+CO2 scattering is the most important sink of the escaping hot O flux that strongly affects estimates of total atomic oxygen escape rate obtained from MAVEN data, as outlined above. New O+CO2 elastic and inelastic cross sections were constructed from first principles for collision energies from 0.01-5 eV. O(3P)+CO2 collisions were described using new 4D electronic potential energy surfaces (PESs) for restricted geometries. Three lowest triplet states of CO2 were considered for symmetric, bending, and asymmetric vibrational modes. The cross sections were determined from close-coupling calculations performed on individual surfaces using vibrating rotor approximations for CO2 molecule. Differential and momentum transfer cross sections, of use in energy transfer atmospheric models, were also constructed. The obtained cross sections are in excellent agreement with published results. Our new elastic cross sections are significantly smaller than assumed in earlier atmospheric escape models, resulting in larger photochemical escape rate of O from Mars, as estimated from a 1D escape model. Our results agree well with estimated hot O flux (about 9 x 10(sup 25) s(sup -1)) derived from altitude profiles of 130.4 nm O emission measured up to ~4,000 km by the Imaging Ultraviolet Spectrograph on MAVEN orbiter.

Atmospheric↗
NASA NTRS2015

Empirical InfraRed Line Lists for Five SO2 Isotopologues: 32/33/34/36S16O2 and 32S18O2

Using the latest published, empirically refined potential energy surface (PES) Ames-1 and purely ab initio CCSD(T)/aug-cc-pVQZ dipole moment surface (DMS), we have computed Infrared line lists for five symmetric isotopologues of sulfur dioxide: 32S16O2 (626), 33S16O2 (636), 34S16O2 (646), 36S16O2 (666), and 32S18O2 (828). The line lists are based on J=0-80 rovibrational variational calculations with E’ ≤ 8000 cm-1. The 34S16O2 and 33S16O2 line lists are compared to the experiment-based models in the HIgh-resolution TRANsmission molecular absorption database (HITRAN2012, http://www.cfa.harvard.edu/hitran/) and the Cologne Database for Molecular Spectroscopy, CDMS (http://www.astro.uni-koeln.de/cdms/). The accuracy for computed 646 band origins is similar to what has been reported for the main isotopologue, i.e. 0.01 – 0.03 cm-1 for bands up to 5500 cm-1. For rovibrational transitions, the 646 line position and intensity deviation patterns are much simpler and more self-consistent than those of the main isotopologue 626. The discrepancies are mainly found for higher Ka/J transitions. 626 and 646 exhibit comparable line position and intensity agreement for lower Ka/J transitions. The line position deviations for the 636 purely rotational band are parallel to those of 626 and 646, while its line intensity deviations do not show branching patterns as we found in the 626 and 646 cases. Predictions for the other minor isotopologues are expected to exhibit similar accuracy. These line lists are accurate enough to provide alternatives for missing bands of 626 and the minor isotopologues. It may significantly facilitate the laboratory spectroscopic measurement and analysis, as well as to identify these isotopologues in various astrophysical environments.

Sulfur Dioxide↗
NASA NTRS2022

Anisotropy Measurements from the Near-Threshold Photodissociation of the N2−NO Complex

We have used velocity map ion imaging to measure the angular anisotropy of the NO (A) products from the photodissociation of the N 2 −NO complex. Our experiment ranged from 108 to 758 cm −1 above the threshold energy to form NO (A)+N 2 (X) products, and these measurements reveal, for the first time, a strong angular anisotropy from photodissociation. At 108 cm −1 above the photodissociation threshold, we observed NO (A) photoproducts recoil preferentially perpendicular to the laser polarization axis with an average anisotropy parameter, β = −0.25; however, as the available energy was increased, the anisotropy increased, and at 758 cm −1 above the threshold energy, we found an average β = +0.28. The observed changes in the angular anisotropy of the NO (A) photoproduct are qualitatively similar to those observed for the photodissociation of the Ar−NO complex and likely result from changes in the region of the excited state potential energy surface accessed during the electronic excitation. At the lowest available energy, we also noted a large contribution from hot band excitation; however, this contribution decreased as the available energy increased. The outsized contribution at the lowest available energy may result from hot bands having better Franck−Condon overlap with the excited electronic state near threshold. Finally, we contrast the experimental center of mass translational energy distribution with a statistical energy distribution determined from phase space theory. The experimental and statistical distributions show pronounced disagreement, particularly at low kinetic energies, with the experimental one showing less dissociation resulting in high rotational levels of the fragments.

Energy↗
NASA NTRS2021

What it Takes to Compute Highly Accurate Rovibrational Line Lists for Use in Astrochemistry

We review the Best Theory + Reliable High-resolution Experiment (BTRHE) strategy for obtaining highly accurate molecular rovibrational line lists with InfraRed (IR) intensities. The need for highly accurate molecular rovibrational line lists is two-fold: a) assignment of the many rovibrational lines for common stable molecules especially those that exhibit a large amplitude motion, such as NH3, or have a high density of states such as SO2; and b) characterization of the atmospheres of exoplanets which will be one of the main areas of research in astronomy in the coming decades. The first motivation arises due to the need to eliminate lines due to common molecules in an astronomical observation in order to identify lines from new molecules, while the second motivation arises due to the need to obtain accurate molecular opacities in order to characterize the atmosphere of an exoplanet. The BTRHE strategy first consists of using high quality ab initio quantum chemical methods to obtain a global potential energy surface (PES) and dipole moment surface (DMS) that contains the proper physics. The global PES is then refined using a subset of the reliable high-resolution experimental data. The refined PES then gives energy level predictions to an accuracy similar to the reproduction accuracy of the experimental data used in the refinement step in the interpolation region (i.e., within the range of the experimental data used in the refinement step). The accuracy of the energy levels will slowly degrade as they are extrapolated to spectral regions beyond the high-resolution experimental data used in the refinement step. However, because the degradation is slow, the predicted energy levels can be used to assign new high-resolution experiments, and the data from these can then be used in a subsequent refinement step. In this way, the global PES eventually can yield highly accurate energy levels for all desired spectral regions including to very high energies and high J values. We show that IR intensities computed with the BTRHE rovibrational wavefunctions and the DMS can be very accurate provided one has minimized the fitting error of the DMS and tested the completeness of the DMS. Some examples of our work on NH3, CO2, and SO2 are given to highlight the usefulness of the BTRHE strategy and to provide ideas on how to further improve its predictive power in the future. In particular, it is shown how successive refinement steps, once new high-resolution data is available, can lead to PESs that yield highly accurate transition energies to larger spectral regions. The importance of including non-adiabatic corrections to reduce the J-dependence of errors for H-containing molecules is shown with work on NH3. Another very important aspect of the BTRHE approach is the consistency across isotopologues, which allows for highly accurate line lists for any isotopologue once one is obtained for the main isotopologue (which has more high-resolution data available for refinement).

Xinchuan Huang↗
NASA NTRS2020

Exploring the Limits of the Data-Model-Theory Synergy: “Hot” MW Transitions for Rovibrational IR Studies

In order to further improve the accuracy of rovibrational IR line lists generated from the “Best Theory +Reliable High-resolution Experiment” (BTRHE) strategy from 0.01-0.05 cm-1, or 300-1500 MHz, to ~10 MHz, we explore the current limits of the Data-Model-Theory synergy by examining the accuracy and consistency of existing data, then propose that “hot” bands in microwave (MW) spectra is the solution we need for future enhancements. The Ames SO2 J=0-20 rovibrational energy levels computed on the semi-empirically refined Ames-2 potential energy surface (PES) are fit to the Effective Hamiltonian (EH) model regularly used in the experimental infrared (IR) analysis for SO2 isotopologues. In the fitted EH(Ames) model, the rotational constants A/B/C and all 5 quartic centrifugal distortion constants display clear, systematic, and consistent patterns along the vibrational state energy or quanta. Such consistent patterns may facilitate the vibrational assignments for MW hot bands and extract more information from high temperature MW spectra. Some EH(Expt) analyses were carried out with the lowest order Coriolis Coupling term, C1. Their constants should not be directly compared with other EH(Expt) and EH(Ames) results. After excluding them, our  = EH(Ames)- EH(Expt) analyses for 5 isotopologues (626, 636, 646, 628 and 828) indicates some loss of accuracy and consistency starting from vibrational states as low as 22 or 1000 cm-1. Some EH parameters, e.g. K, may have relative deviations as large as 50-100% and totally lose any recognizable patterns. This simply means that current EH(Expt) models do not have the system-wide consistency we need to further refine the EH(Ames) and Ames rovibrational IR line lists. A large part of such defects are probably inherited from the limited precision of experimental line positions, i.e. 1E-3 ~ 1E-4 cm-1, or 3-30 MHz. This is confirmed in a series of truncation tests using the Ames data. Although the EH(Ames) consistency may help identify unreliable rovibrational EH(Expt) parameters, and make reliable predictions for minor isotopologues and unobserved vibrational bands, we believe only the highly accurate “hot” MW transitions can provide real enhancements for EH(Expt) accuracy and consistency. “Hot” MW spectra should play a more significant role in the future synergy of Data, Model, and Theory in the field of rovibrational IR studies.

Xinchuan Huang↗
NASA NTRSDate not supplied

Utilizing Experimental-based Line Positions for Semi-empirical IR Line Lists

Molecular IR line lists computed from semi-empirically refined ab initio potential energy surface and high quality ab initio dipole moment surface may have line position accuracy typically within 0.01 - 0.05 cm−1. For high resolution spectroscopy databases, it is necessary to integrate the computed theoretical IR intensities with the line positions accurately determined from experiments and (or) reliably derived from Effective Hamiltonian models. However, it is not a trivial task, and the choices in reality are heavily contingent upon the coverage, consistency, and accuracy of the data available for a specific molecule or isotopologue. We will present several approaches applied to recent Ames IR line lists of CO 2 1, N 2 O 2 , and OCS to demonstrate how challenging such integration may become, and what we have learned From these projects. Discussions will focus on the relation between each different scenario and the corresponding choice or solution, including their advantages and limits, how the semi-empirical IR line lists can help, and what else may be needed for future improvements. We will emphasize that the “Theory+Experiment” synergy may still play significant role in the determination of the best line positions.

Effective Hamiltonian models↗
NASA NTRSDate not supplied

Reliable IR Line Lists for N2O and OCS Isotopologues for High Temperature Applications up to 2000 K

To address the needs of planetary and exoplanetary atmosphere studies towards JWST, HWO and future IR astronomy missions, theoretical IR line lists of Carbonyl Sulfide (OCS) and Nitrous Oxide (N2O) are generated to provide complete, reliable, and consistent IR reference data in the range of 0 - 15,000 cm-1, or wavelength longer than 0.67 Å, for temperature up to 2000 K. The line lists are quantum mechanically computed using spectroscopically accurate potential energy surfaces empirically refined with reliable high-resolution experimental data, and high-quality dipole moment surfaces constructed at CCSD(T) level with aug-cc-pVT/Q/5(+d)Z basis and extrapolations to complete basis set limit. The Ames N2O line list includes the primary isotopologue 14N2 16O (denoted "446") with J=0-210 and E' up to 25,000 cm-1, and 11 minor isotopologues of 15N and 17/18O with J=0-150 and lower E' cutoffs. For OCS, Ames "natural" 296 K line list includes 12 most abundant isotopologues. The intensity contribution of OCS minor isotopologues at higher temperature is non-negligible, e.g., 8.6 ± 4.9% at 1000 K. The 16O12C32S (93.7%) and 16O12C34S (4.2%) transitions are computed with J=0-240 and E'< 22,000 cm-1, plus J=0-150 and E' < 16,000 cm-1 for five more minor isotopologues. Reliable Effective Hamiltonian models (e.g., NOSL-296 for 14N2 16O), experimental data collections, and reliable HITRAN or MARVEL levels are combined with quantum rovibrational calculations to deliver accurate and consistent empirical corrections for line positions. The reliability and accuracy of our IR intensity predictions are demonstrated by comparison with available databases, line lists, and experiments, at both room temperature and at high temperature. The line lists will be freely available to public online at ZENODO and NAS data portal. We acknowledge the financial support from NASA XRP and APRA programs.

Planetary↗
NASA NTRSDate not supplied

Ab Initio Electronic Structure Calculations of CNN for CN Excitation Studies

The CN molecule is an important contributor to radiative heat flux in shock layers around vehicles entering Titan’s atmosphere. Current data for heavy particle (de)excitation rate coefficients of CN leads to uncertainties in the population of CN in its first and second excited states. This in turn leads to uncertainties in the radiative heat flux predicted by Computational Fluid Dynamics (CFD) simulations of Titan atmospheric entry. This work performs ab initio electronic structure calculations of the CNN complex to create Potential Energy Surfaces (PESs) that correlate to the ground and first and second excited states of CN. Specifically, the state combinations of CN(X,A,B) + N(4S𝑜) correlate to six states of CNN (three Quintet A” and three Triplet A”). Initial calculations of these states suggest that heavy particle (de)excitation of CN by N atoms is likely to proceed through collinear geometries on triplet surfaces. Complete PESs will show all of the reaction pathways in detail, and will be used in nonadiabatic dynamics calculations to evaluate improved rate coefficients and reduce uncertainty in the radiative heat flux during Titan entry.

Eric C Geistfeld↗