Search

Most low-mass stars form in stellar clusters that also contain massive stars, which are sources of far-ultraviolet (FUV) radiation. Theoretical models predict that this FUV radiation produces photo-dissociation regions (PDRs) on the surfaces of protoplanetary disks around low-mass stars, impacting planet formation within the disks. We report JWST and Atacama Large Millimetere Array observations of a FUV-irradiated protoplanetary disk in the Orion Nebula. Emission lines are detected from the PDR; modelling their kinematics and excitation allows us to constrain the physical conditions within the gas. We quantify the mass-loss rate induced by the FUV irradiation, finding it is sufficient to remove gas from the disk in less than a million years. This is rapid enough to affect giant planet formation in the disk.

The synergy between high-resolution rotational spectroscopy and quantum-chemical calculations is essential for exploring future detection of molecules, especially when spectroscopy parameters are not available yet. By using highly correlated ab initio quartic force fields (QFFs) from explicitly correlated coupled cluster theory, a complete set of rotational constants and centrifugal distortion constants for D$_2$CS and cis/trans-DCSD isomers have been produced. Comparing our new ab initio results for D$_2$CS with new rotational spectroscopy laboratory data for the same species, the accuracy of the computed B and C rotational constants is within 0.1% while the A constant is only slightly higher. Additionally, quantum chemical vibrational frequencies are also provided, and these spectral reference data and new experimental rotational lines will provide additional references for potential observation of these deuterated sulfur species with either ground-based radio telescopes or space-based infrared observatories., Comment: 15 pages, 2 figures and 1 appendix table; published in Molecular Physics (Tim Lee Memorial Issue)

(Abridged) We investigate the impact of radiative feedback from massive stars on their natal cloud and focus on the transition from the HII region to the atomic PDR (crossing the ionisation front (IF)), and the subsequent transition to the molecular PDR (crossing the dissociation front (DF)). We use high-resolution near-IR integral field spectroscopic data from NIRSpec on JWST to observe the Orion Bar PDR as part of the PDRs4All JWST Early Release Science Program. The NIRSpec data reveal a forest of lines including, but not limited to, HeI, HI, and CI recombination lines, ionic lines, OI and NI fluorescence lines, Aromatic Infrared Bands (AIBs including aromatic CH, aliphatic CH, and their CD counterparts), CO2 ice, pure rotational and ro-vibrational lines from H2, and ro-vibrational lines HD, CO, and CH+, most of them detected for the first time towards a PDR. Their spatial distribution resolves the H and He ionisation structure in the Huygens region, gives insight into the geometry of the Bar, and confirms the large-scale stratification of PDRs. We observe numerous smaller scale structures whose typical size decreases with distance from Ori C and IR lines from CI, if solely arising from radiative recombination and cascade, reveal very high gas temperatures consistent with the hot irradiated surface of small-scale dense clumps deep inside the PDR. The H2 lines reveal multiple, prominent filaments which exhibit different characteristics. This leaves the impression of a "terraced" transition from the predominantly atomic surface region to the CO-rich molecular zone deeper in. This study showcases the discovery space created by JWST to further our understanding of the impact radiation from young stars has on their natal molecular cloud and proto-planetary disk, which touches on star- and planet formation as well as galaxy evolution., Comment: 52 pages, 30 figures, submitted to A&A

The JWST has captured the most detailed and sharpest infrared images ever taken of the inner region of the Orion Nebula, the nearest massive star formation region, and a prototypical highly irradiated dense photo-dissociation region (PDR). We investigate the fundamental interaction of far-ultraviolet photons with molecular clouds. The transitions across the ionization front (IF), dissociation front (DF), and the molecular cloud are studied at high-angular resolution. These transitions are relevant to understanding the effects of radiative feedback from massive stars and the dominant physical and chemical processes that lead to the IR emission that JWST will detect in many Galactic and extragalactic environments. Due to the proximity of the Orion Nebula and the unprecedented angular resolution of JWST, these data reveal that the molecular cloud borders are hyper structured at small angular scales of 0.1-1" (0.0002-0.002 pc or 40-400 au at 414 pc). A diverse set of features are observed such as ridges, waves, globules and photoevaporated protoplanetary disks. At the PDR atomic to molecular transition, several bright features are detected that are associated with the highly irradiated surroundings of the dense molecular condensations and embedded young star. Toward the Orion Bar PDR, a highly sculpted interface is detected with sharp edges and density increases near the IF and DF. This was predicted by previous modeling studies, but the fronts were unresolved in most tracers. A complex, structured, and folded DF surface was traced by the H2 lines. This dataset was used to revisit the commonly adopted 2D PDR structure of the Orion Bar. JWST provides us with a complete view of the PDR, all the way from the PDR edge to the substructured dense region, and this allowed us to determine, in detail, where the emission of the atomic and molecular lines, aromatic bands, and dust originate.

(Abridged) Mid-infrared observations of photodissociation regions (PDRs) are dominated by strong emission features called aromatic infrared bands (AIBs). The most prominent AIBs are found at 3.3, 6.2, 7.7, 8.6, and 11.2 $\mu$m. The most sensitive, highest-resolution infrared spectral imaging data ever taken of the prototypical PDR, the Orion Bar, have been captured by JWST. We provide an inventory of the AIBs found in the Orion Bar, along with mid-IR template spectra from five distinct regions in the Bar: the molecular PDR, the atomic PDR, and the HII region. We use JWST NIRSpec IFU and MIRI MRS observations of the Orion Bar from the JWST Early Release Science Program, PDRs4All (ID: 1288). We extract five template spectra to represent the morphology and environment of the Orion Bar PDR. The superb sensitivity and the spectral and spatial resolution of these JWST observations reveal many details of the AIB emission and enable an improved characterization of their detailed profile shapes and sub-components. While the spectra are dominated by the well-known AIBs at 3.3, 6.2, 7.7, 8.6, 11.2, and 12.7 $\mu$m, a wealth of weaker features and sub-components are present. We report trends in the widths and relative strengths of AIBs across the five template spectra. These trends yield valuable insight into the photochemical evolution of PAHs, such as the evolution responsible for the shift of 11.2 $\mu$m AIB emission from class B$_{11.2}$ in the molecular PDR to class A$_{11.2}$ in the PDR surface layers. This photochemical evolution is driven by the increased importance of FUV processing in the PDR surface layers, resulting in a "weeding out" of the weakest links of the PAH family in these layers. For now, these JWST observations are consistent with a model in which the underlying PAH family is composed of a few species: the so-called 'grandPAHs'., Comment: 25 pages, 10 figures, to appear in A&A

We present the first high spectral resolution mid-infrared survey in the Orion BN/KL region, covering 7.2 to 28.3 micron. With SOFIA/EXES we target the enigmatic source Orion IRc2. While this is in the most prolifically studied massive star-forming region, longer wavelengths and molecular emission lines dominated previous spectral surveys. The mid-infrared observations in this work access different components and molecular species in unprecedented detail. We unambiguously identify two new kinematic components, both chemically rich with multiple molecular absorption lines. The "blue clump" has vLSR = -7.1 \pm 0.7 km/s and the "red clump" 1.4 \pm 0.5 km/s. While the blue and red clumps have similar temperatures and line widths, molecular species in the blue clump have higher column densities. They are both likely linked to pure rotational H2 emission also covered by this survey. This work provides evidence for the scenario that the blue and red clumps are distinct components unrelated to the classic components in the Orion BN/KL region. Comparison to spectroscopic surveys towards other infrared targets in the region show that the blue clump is clearly extended. We analyze, compare, and present in depth findings on the physical conditions of C2H2, 13CCH2, CH4, CS, H2O, HCN, H13CN, HNC, NH3, and SO2 absorption lines and an H2 emission line associated with the blue and red clumps. We also provide limited analysis of H2O and SiO molecular emission lines towards Orion IRc2 and the atomic forbidden transitions [FeII], [SI], [SIII], and [NeII]., Comment: Accepted to ApJ; 44 pages, 14 figures, 13 tables

Massive stars disrupt their natal molecular cloud material through radiative and mechanical feedback processes. These processes have profound effects on the evolution of interstellar matter in our Galaxy and throughout the Universe, from the era of vigorous star formation at redshifts of 1-3 to the present day. The dominant feedback processes can be probed by observations of the Photo-Dissociation Regions (PDRs) where the far-ultraviolet photons of massive stars create warm regions of gas and dust in the neutral atomic and molecular gas. PDR emission provides a unique tool to study in detail the physical and chemical processes that are relevant for most of the mass in inter- and circumstellar media including diffuse clouds, proto-planetary disks and molecular cloud surfaces, globules, planetary nebulae, and star-forming regions. PDR emission dominates the infrared (IR) spectra of star-forming galaxies. Most of the Galactic and extragalactic observations obtained with the James Webb Space Telescope (JWST) will therefore arise in PDR emission. In this paper we present an Early Release Science program using the MIRI, NIRSpec, and NIRCam instruments dedicated to the observations of an emblematic and nearby PDR: the Orion Bar. These early JWST observations will provide template datasets designed to identify key PDR characteristics in JWST observations. These data will serve to benchmark PDR models and extend them into the JWST era. We also present the Science-Enabling products that we will provide to the community. These template datasets and Science-Enabling products will guide the preparation of future proposals on star-forming regions in our Galaxy and beyond and will facilitate data analysis and interpretation of forthcoming JWST observations., Comment: Submitted to PASP

We present the first mid-infrared (MIR) detections of HNC and H13CN in the interstellar medium, and numerous, resolved HCN rovibrational transitions. Our observations span 12.8 to 22.9 micron towards the hot core Orion IRc2, obtained with the Echelon-Cross-Echelle Spectrograph aboard the Stratospheric Observatory for Infrared Astronomy (SOFIA). Exceptional, ~5 km/s, resolution distinguishes individual rovibrational transitions of the HNC and HCN P, Q, and R branches; and the H13CN R branch. This allows direct measurement of the species' excitation temperatures, column densities, and relative abundances. HNC and H13CN exhibit a local standard rest velocity of -7 km/s that may be associated with an outflow from nearby radio source I and an excitation temperature of about 100 K. We resolve two velocity components for HCN, the primary component also being at -7 km/s with temperature 165 K. The hottest component, which had never before been observed, is at 1 km/s with temperature 309 K. This is the closest component to the hot core's centre measured to date. The derived 12C/13C=13 is below expectation for Orion's Galactocentric distance, but the derived HCN/HNC=72 is expected for this extreme environment. Compared to previous sub-mm and mm observations, our SOFIA line survey of this region shows that the resolved MIR molecular transitions are probing a distinct physical component and isolating the chemistry closest to the hot core., Comment: 20 pages, 6 figures, submitted to ApJ, comments welcome

We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize their atmospheric structure, composition, and circulation, from gas giants to rocky planets. However, exoplanet atmospheric models capable of interpreting the upcoming observations are often limited by insufficiencies in the laboratory and theoretical data that serve as critical inputs to atmospheric physical and chemical tools. Here we provide an up-to-date and condensed description of areas where laboratory and/or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry, building off a larger 2016 white paper, and endorsed by the NAS Exoplanet Science Strategy report. Now is the ideal time for progress in these areas, but this progress requires better access to, understanding of, and training in the production of spectroscopic data as well as a better insight into chemical reaction kinetics both thermal and radiation-induced at a broad range of temperatures. Given that most published efforts have emphasized relatively Earth-like conditions, we can expect significant and enlightening discoveries as emphasis moves to the exotic atmospheres of exoplanets., Comment: Submitted as an Astro2020 Science White Paper

The infrared (IR) emission of polycyclic aromatic hydrocarbons (PAHs) permeates our universe; astronomers have detected the IR signatures of PAHs around many interstellar objects. The IR emission of interstellar PAHs differs from their emission as seen under conditions on Earth, as they emit through a collisionless cascade down through their excited vibrational states from high internal energies. The difficulty in reproducing interstellar conditions in the laboratory results in a reliance on theoretical techniques. However, the size and complexity of PAHs requires careful consideration when producing the theoretical spectra. In this work we outline the theoretical methods necessary to lead to a fully theoretical IR cascade spectra of PAHs including: an anharmonic second order vibrational perturbation theory (VPT2) treatment; the inclusion of Fermi resonances through polyads; and the calculation of anharmonic temperature band shifts and broadenings (including resonances) through a Wang--Landau approach. We also suggest a simplified scheme to calculate vibrational emission spectra that retains the essential characteristics of the full IR cascade treatment and can directly transform low temperature absorption spectra in IR cascade spectra. Additionally we show that past astronomical models were in error in assuming a 15 cm$^{-1}$ correction was needed to account for anharmonic emission effects.

Aims. Infrared (IR) spectroscopy is a powerful tool to study molecules in space. A key issue in such analyses is understanding the effect that temperature and anharmonicity have on different vibrational bands, and thus interpreting the IR spectra for molecules under various conditions. Methods. We combined second order vibrational perturbation theory and the Wang-Landau random walk technique to produce accurate IR spectra of highly excited PAHs. We fully incorporated anharmonic effects, such as resonances, overtones, combination bands, and temperature effects. Results. The results are validated against experimental results for the pyrene molecule (C16H10). In terms of positions, widths, and relative intensities of the vibrational bands, our calculated spectra are in excellent agreement with gas-phase experimental data.

We present high-spectral resolution observations from 12.96 - 13.33 microns towards Orion IRc2 using the mid-infrared spectrograph, EXES, on SOFIA. These observations probe the physical and chemical conditions of the Orion hot core, which is sampled by a bright, compact, mid-infrared background continuum source in the region, IRc2. All ten of the rovibrational C2H2 transitions expected in our spectral coverage, are detected with high S/N, yielding continuous coverage of the R-branch lines from J=9-8 to J=18-17, including both ortho and para species. Eight of these rovibrational transitions are newly reported detections. The isotopologue, 13CCH2, is clearly detected with high signal-to-noise. This enabled a direct measurement of the 12C/13C isotopic ratio for the Orion hot core of 14 +/- 1 and an estimated maximum value of 21. We also detected several HCN rovibrational lines. The ortho and para C2H2 ladders are clearly separate and tracing two different temperatures, 226 K and 164 K, respectively, with a non-equilibrium ortho to para ratio (OPR) of 1.7 +\- 0.1. Additionally, the ortho and para V_LSR values differ by about 1.8 +/- 0.2 km/s, while, the mean line widths differ by 0.7 +/- 0.2 km/s, suggesting that these species are not uniformly mixed along the line of sight to IRc2. We propose that the abnormally low C2H2 OPR could be a remnant from an earlier, colder phase, before the density enhancement (now the hot core) was impacted by shocks generated from an explosive event 500 yrs ago., Comment: accepted for publication in The Astrophysical Journal

In this work we report on high-resolution IR absorption studies that provide a detailed view on how the peripheral structure of irregular polycyclic aromatic hydrocarbons (PAHs) affects the shape and position of their 3 {\mu}m absorption band. To this purpose we present mass-selected, high-resolution absorption spectra of cold and isolated phenanthrene, pyrene, benz[a]antracene, chrysene, triphenylene, and perylene molecules in the 2950-3150 cm-1 range. The experimental spectra are compared with standard harmonic calculations, and anharmonic calculations using a modified version of the SPECTRO program that incorporates a Fermi resonance treatment utilizing intensity redistribution. We show that the 3 {\mu}m region is dominated by the effects of anharmonicity, resulting in many more bands than would have been expected in a purely harmonic approximation. Importantly, we find that anharmonic spectra as calculated by SPECTRO are in good agreement with the experimental spectra. Together with previously reported high-resolution spectra of linear acenes, the present spectra provide us with an extensive dataset of spectra of PAHs with a varying number of aromatic rings, with geometries that range from open to highly-condensed structures, and featuring CH groups in all possible edge configurations. We discuss the astrophysical implications of the comparison of these spectra on the interpretation of the appearance of the aromatic infrared 3 {\mu}m band, and on features such as the two-component emission character of this band and the 3 {\mu}m emission plateau., Comment: Accepted for publication on The Astrophysical Journal

Sulphur dioxide is well-known in the atmospheres of planets and satellites, where its presence is often associated with volcanism, and in circumstellar envelopes of young and evolved stars as well as the interstellar medium. This work presents a line list of 1.3 billion $^{32}$S$^{16}$O$_2$ vibration-rotation transitions computed using an empirically-adjusted potential energy surface and an ab initio dipole moment surface. The list gives complete coverage up to 8000 cm$^{-1}$ (wavelengths longer than 1.25 $\mu$m) for temperatures below 2000 K. Infrared absorption cross sections are recorded at 300 and 500 C are used to validated the resulting ExoAmes line list. The line list is made available in electronic form as supplementary data to this article and at www.exomol.com., Comment: Submitted to MNRAS

We report on an experimental and theoretical investigation of the importance of anharmonicity in the 3 micron CH stretching region of Polycyclic Aromatic Hydrocarbon (PAH) molecules. We present mass-resolved, high-resolution spectra of the gas-phase cold (~4K) linear PAH molecules naphthalene, anthracene, and tetracene. The measured IR spectra show a surprisingly high number of strong vibrational bands. For naphthalene, the observed bands are well separated and limited by the rotational contour, revealing the band symmetries. Comparisons are made to the harmonic and anharmonic approaches of the widely used Gaussian software. We also present calculated spectra of these acenes using the computational program SPECTRO, providing anharmonic predictions enhanced with a Fermi-resonance treatment that utilises intensity redistribution. We demonstrate that the anharmonicity of the investigated acenes is strong, dominated by Fermi resonances between the fundamental and double combination modes, with triple combination bands as possible candidates to resolve remaining discrepancies. The anharmonic spectra as calculated with SPECTRO lead to predictions of the main modes that fall within 0.5% of the experimental frequencies. The implications for the Aromatic Infrared Bands, specifically the 3 micron band are discussed., Comment: v.2 (02/11/15) Updated: references, figures and affiliation

It has been shown that rotational lines observed in the Horsehead nebula PDR are probably not caused by l-C3H+, as was originally suggested. In the search for viable alternative candidate carriers, quartic force fields are employed here to provide highly accurate rotational constants, as well as fundamental vibrational frequencies, for another candidate carrier: 1 ^1A' C3H-. The ab initio computed spectroscopic constants provided in this work are, compared to those necessary to define the observed lines, as accurate as the computed spectroscopic constants for many of the known interstellar anions. Additionally, the computed Deff for C3H- is three times closer to the D deduced from the observed Horsehead nebula lines relative to l-C3H+. As a result, 1 ^1A' C3H- is a more viable candidate for these observed rotational transitions. It has been previously proposed that at least C6H- may be present in the Horsehead nebular PDR formed by way of radiative attachment through its dipole-bound excited state. C3H- could form in a similar way through its dipole-bound state, but its valence excited state increases the number of relaxation pathways possible to reach the ground electronic state. In turn, the rate of formation for C3H- could be greater than the rate of its destruction. C3H- would be the seventh confirmed interstellar anion detected within the past decade and the first CnH- molecular anion with an odd n., Comment: 26 pages, 1 figure, 4 tables

It is hypothesized that the A ^1B_1 <- X ^1A' excitation into the dipole-bound state of the cyanomethyl anion (CH2CN-) is proposed as the carrier for one diffuse interstellar band. However, this particular molecular system has not been detected in the interstellar medium even though the related cyanomethyl radical and the isoelectronic ketenimine molecule have been found. In this study we are employing the use of proven quartic force fields and second-order vibrational perturbation theory to compute accurate spectroscopic constants and fundamental vibrational frequencies for ^1A' CH2CN- in order to assist in laboratory studies and astronomical observations.

Very recently, molecular rotational transitions observed in the photon-dominated region of the Horsehead nebula have been attributed to l-C3H+. In an effort to corroborate this finding, we employed state-of-the art and proven high-accuracy quantum chemical techniques to compute spectroscopic constants for this cation and its isotopologues. Even though the B rotational constant from the fit of the observed spectrum and our computations agree to within 20 MHz, a typical level of accuracy, the D rotational constant differs by more than 40%, while the H rotational constant differs by three orders of magnitude. With the likely errors in the rotational transition energies resulting from this difference in D on the order of 1 MHz for the lowest observed transition (J = 4 <- 3) and growing as J increases, the assignment of the observed rotational lines from the Horsehead nebula to l-C3H+ is questionable.

From large basis set coupled cluster calculations and a minor empirical adjustment, an anharmonic force field for silane has been derived that is consistently of spectroscopic quality ($\pm 1 cm^{-1}$ on vibrational fundamentals) for all isotopomers of silane studied. Inner-shell polarization functions have an appreciable effect on computed properties and even on anharmonic corrections. From large basis set coupled cluster calculations and extrapolations to the infinite-basis set limit, we obtain TAE_0=303.80 \pm 0.18 kcal/mol, which includes an anharmonic zero-point energy (19.59 kcal/mol), inner-shell correlation (-0.36 kcal/mol), scalar relativistic corrections (-0.70 kcal/mol), and atomic spin-orbit corrections (-0.43 kcal/mol). In combination with the recently revised \HVSI{0}, we obtain $\Delta H^{\circ}_{f,0}[SiH_4(g)]=9.9 \pm 0.4 kcal/mol$, in between the two established experimental values., Comment: Mol. Phys., in press

Numerous laboratory studies of astrophysical ice analogues have shown that their exposure to ionizing radiation leads to the production of large numbers of new, more complex compounds, many of which are of astrobiological interest. We show here that the irradiation of astrophysical ice analogues containing H2O, CH3OH, CO, and NH3 yields quantities of hexamethylenetetramine-methanol (hereafter HMT-methanol; C7N4H14O) that are easily detectible in the resulting organic residues. This molecule differs from simple HMT, which is known to be abundant in similar ice photolysis residues, by the replacement of a peripheral H atom with a CH2OH group. As with HMT, HMT-methanol is likely to be an amino acid precursor. HMT has tetrahedral (T(d)) symmetry, whereas HMT-methanol has C1 symmetry. We report the computed expected infrared spectra for HMT and HMT-methanol obtained using ab initio quantum chemistry methods and show that there is a good match between the observed and computed spectra for regular HMT. Since HMT-methanol lacks the high symmetry of HMT, it produces rotational transitions that could be observed at longer wavelengths, although establishing the exact positions of these transitions may be challenging. It is likely that HMT-methanol represents an abundant member of a larger family of functionalized HMT molecules that may be present in cold astrophysical environments.

Over the last several years, our group has been involved in developing approaches to compute highly accurate spectroscopic constants and vibrational frequencies for small transient molecules that may be used in the interpretation and assignment of high-resolution laboratory experiments as well as high-resolution astronomical spectra. Additionally, we have used the computed spectroscopic constants to simulate purely rotational and rovibrational spectra so that these may be compared directly with high- resolution astronomical observations, and we have worked on developing approaches that can be applied to much larger molecules, such as PAHs, where we can explicitly determine an harmonic corrections to vibrational modes as well as take into account intensity sharing due to resonances. Another part of our work in spectroscopic signatures involves computing highly accurate line lists for common molecules, such as CO2, SO2, and NH3, which occur in many astrophysical environments, including the atmospheres of exoplanets, and often need to have their lines identified in high- resolution observations in order to determine which lines are due to other molecules. In order to characterize the atmospheres of hot exoplanets, these line lists need to be very accurate and extend to very high energies. I will discuss our latest work in these areas of astrochemical spectroscopy research.

Carbon-carbon multiply bonded systems are improperly described with standard correlation methods and basis sets. For computations of vibrational modes, the out-of- plane bends can be reported as imaginary at worst or simply too low at best. Utilizing the simplest of aromatic structures (cyclopropenylidene) and various levels of theory, this work diagnoses this known behavior for the first time. A combined 1-particle and n-particle basis set effect conspire to produce these non-physical results. When moving from sp2 to sp3 hybridization in the carbon atoms, the larger number of basis functions overcorrects the energy. This is exacerbated by correlation methods. These allow for occupation of the π and π∗ orbitals in the expanded wave function that combine with the hydrogen s orbitals. As a result, the improperly described space can be further and non-physically stabilized by post-Hartree-Fock correlation. This represents a fundamental problem with at least Hartree-Fock based methods of all flavors in describing carbon. Beyond being a flaw in quantum chemical theory, other repercussions will be present in computations regarding spectroscopy as well as energy and environmental studies where highly-accurate hydrocabon vibrational transitions or thermochemical data are needed.

In this work we report on high-resolution IR absorption studies that provide a detailed view on how the peripheral structure of irregular polycyclic aromatic hydrocarbons (PAHs) affects the shape and position of their 3-micrometers absorption band. To this purpose we present mass-selected, high-resolution absorption spectra of cold and isolated phenanthrene, pyrene, benz[a]antracene, chrysene, triphenylene, and perylene molecules in the 2950-3150 per cm range. The experimental spectra are compared with standard harmonic calculations, and anharmonic calculations using a modified version of the SPECTRO program that incorporates a Fermi resonance treatment utilizing intensity redistribution. We show that the 3-micrometers region is dominated by the effects of anharmonicity, resulting in many more bands than would have been expected in a purely harmonic approximation. Importantly, we find that anharmonic spectra as calculated by SPECTRO are in good agreement with the experimental spectra. Together with previously reported high-resolution spectra of linear acenes, the present spectra provide us with an extensive dataset of spectra of PAHs with a varying number of aromatic rings, with geometries that range from open to highly-condensed structures, and featuring CH groups in all possible edge configurations. We discuss the astrophysical implications of the comparison of these spectra on the interpretation of the appearance of the aromatic infrared 3-micrometers band, and on features such as the two-component emission character of this band and the 3-micrometers emission plateau.

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

Nucleobases are the carriers of the genetic information in ribonucleic acid and deoxyribonucleic acid (DNA) for all life on Earth. Their presence in meteorites clearly indicates that compounds of biological importance can form via non-biological processes in extraterrestrial environments. Recent experimental studies have shown that the pyrimidine-based nucleobases uracil and cytosine can be easily formed from the ultraviolet irradiation of pyrimidine in H2O-rich ice mixtures that simulate astrophysical processes. In contrast, thymine, which is found only in DNA, is more difficult to form under the same experimental conditions, as its formation usually requires a higher photon dose. Earlier quantum chemical studies confirmed that the reaction pathways were favorable provided that several H2O molecules surrounded the reactants. However, the present quantum chemical study shows that the formation of thymine is limited because of the inefficiency of the methylation of pyrimidine and its oxidized derivatives in an H2O ice, as supported by the laboratory studies. Our results constrain the formation of thymine in astrophysical environments and thus the inventory of organic molecules delivered to the early Earth and have implications for the role of thymine and DNA in the origin of life.

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.

We report on an experimental and theoretical investigation of the importance of anharmonicity in the 3 micrometers CH stretching region of Polycyclic Aromatic Hydrocarbon (PAH) molecules. We present mass-resolved, high-resolution spectra of the gas-phase cold ((is) approximately 4K) linear PAH molecules naphthalene, anthracene, and tetracene. The measured IR spectra show a surprisingly high number of strong vibrational bands. For naphthalene, the observed bands are well separated and limited by the rotational contour, revealing the band symmetries. Comparisons are made to the harmonic and anharmonic approaches of the widely used Gaussian software. We also present calculated spectra of these acenes using the computational program SPECTRO, providing anharmonic predictions enhanced with a Fermi-resonance treatment that utilizes intensity redistribution. We demonstrate that the anharmonicity of the investigated acenes is strong, dominated by Fermi resonances between the fundamental and double combination modes, with triple combination bands as possible candidates to resolve remaining discrepancies. The anharmonic spectra as calculated with SPECTRO lead to predictions of the main modes that fall within 0.5% of the experimental frequencies. The implications for the Aromatic Infrared Bands, specifically the 3-m band are discussed.

Polycyclic Aromatic Hydrocarbons (PAHs) play a significant role in the chemistry of the interstellar medium (ISM) as well as in hydrocarbon combustion. These molecules can have high levels of diversity with the inclusion of heteroatoms and the addition or removal of hydrogens to form charged or radical species. There is an abundance of data on the cationic forms of these molecules, but there have been many fewer studies on the anionic species. The present study focuses on the anionic forms of deprotonated PAHs. It has been shown in previous work that PAHs containing nitrogen heteroatoms (PANHs) have the ability to form valence excited states giving anions electronic absorption features. This work analyzes how the isoelectronic pure PAHs behave under similar structural constructions. Singly-deprotonated forms of benzene, naphthalene, anthracene, and teteracene classes are examined. None of the neutral-radicals possess dipole moments large enough to support dipole-bound excited states in their corresponding closed-shell anions. Even though the PANH anion derivatives support valence excited states for three-ringed structures, it is not until four-ringed structures of the pure PAH anion derivatives that valence excited states are exhibited. However, anisotropically-extended PAHs larger than tetracene will likely exhibit valence excited states. The relative energies for the anion isomers are very small for all of the systems in this study.

A semi-global potential energy surface (PES) and quartic force field (QFF) based on fitting high-level electronic structure energies are presented to describe the structures and spectroscopic properties of NNHNN+. The equilibrium structure of NNHNN+ is linear with the proton equidistant between the two nitrogen groups and thus of D(sub ∞h) symmetry. Vibrational second-order perturbation theory (VPT2) calculations based on the QFF fails to describe the proton “rattle” motion, i.e., the antisymmetric proton stretch, due to the very flat nature of PES around the global minimum, but performs properly for other modes with sharper potential wells. Vibrational self-consistent field/virtual state configuration interaction (VSCF/VCI) calculations using a version of MULTIMODE without angular momentum terms successfully describe this motion and predict the fundamental to be at 759 cm(exp -1). This is in good agreement with the value of 746 cm(exp -1) from a fixed-node diffusion Monte Carlo calculation and the experimental Ar-tagged result of 743 cm(exp -1). Other VSCF/VCI energies are in good agreement with other experimentally reported ones. Both double-harmonic intensity and rigorous MULTIMODE intensity calculations show the proton transfer fundamental has a very strong intensity.

Silacyclopropynylidene, SiC2, is a known and highly abundant circumstellar molecule. Its spectrum has been established as a major component of lines observed toward the carbon-rich star IRC +10216 (CW Leonis). It has been detected in its low-lying v(sub 3) = 1 and 2 vibrational states as well as in various isotopic compositions. Increasing sensitivity and spatial resolution will enable many more emission or absorption lines to be detected. In order to detect new molecular species, unassigned lines of known species must be identified. This work uses established ab initio quartic force fields to produce data necessary for this classification of lines related to SiC2. Agreement between the theoretical vibrational frequencies and known rotational and spectroscopic constants is quite good, as good as 5 cm(exp -1) and 3 MHz, respectively in some cases. Additionally, experimentally unknown vibrational frequencies and rotational constants are provided for the first overtones and combination bands in addition to 3(sub v3), the second overtone of the low-lying antisymmetric stretch/carbide rotation mode. Frequencies of v(sub 3) = 3 low-J rotational transitions of the main isotopic species are also estimated from published data for v(sub 3) ≤ 2. Further, we determine rotational and centrifugal distortion parameters for which in most cases vibrational effects due to the v(sub 3) mode were reduced to first, and in several cases also to second order. These values may approximate equilibrium values better than the ground state values. The data produced herein will aid in the experimental and observational characterization of this known astromolecule in order to identify some of the unassigned lines for a known entity.

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.

A purely ab initio potential energy surface (PES) was refined with selected (32)S(16)O2 HITRAN data. Compared to HITRAN, the root-mean-squares error (RMS) error for all J=0-80 rovibrational energy levels computed on the refined PES (denoted Ames-1) is 0.013 cm(exp -1). Combined with a CCSD(T)/aug-cc-pV(Q+d)Z dipole moment surface (DMS), an infrared (IR) line list (denoted Ames-296K) has been computed at 296K and covers up to 8,000 cm(exp -1). Compared to the HITRAN and CDMS databases, the intensity agreement for most vibrational bands is better than 85-90%. Our predictions for (34)S(16)O2 band origins, higher energy (32)S(16)O2 band origins and missing (32)S(16)O2 IR bands have been verified by most recent experiments and available HITRAN data. We conclude that the Ames-1 PES is able to predict (32/34)S(16)O2 band origins below 5500 cm(exp -1) with 0.01-0.03 cm(exp -1) uncertainties, and the Ames-296K line list provides continuous, reliable and accurate IR simulations. The Ka-dependence of both line position and line intensity errors is discussed. The line list will greatly facilitate SO2 IR spectral experimental analysis, as well as elimination of SO2 lines in high-resolution astronomical observations.

The spectroscopic constants and vibrational frequencies for the 1(sup 3)A' states of HNC, DNC, HOC+, and DOC+ are computed and discussed in this work. The reliable CcCR quartic force field based on high-level coupled cluster ab initio quantum chemical computations is exclusively utilized to provide the anharmonic potential. Then, second order vibrational perturbation theory and vibrational configuration interaction methods are employed to treat the nuclear Schroedinger equation. Second-order perturbation theory is also employed to provide spectroscopic data for all molecules examined. The relationship between these molecules and the corresponding 1(sup 3)A' HCN and HCO+ isomers is further developed here. These data are applicable to laboratory studies involving formation of HNC and HOC+ as well as astronomical observations of chemically active astrophysical environments.

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.

Through established, highly-accurate ab initio quartic force fields (QFFs), a complete set of fundamental vibrational frequencies, rotational constants, and rovibrational coupling and centrifugal distortion constants have been determined for both the cyclic 1(sup 1) 1A' and bent 2(sup 1)A' DCCN, H(C13)CCN, HC(C-13)N, and HCC(N-15) isotopologues of HCCN. Spectroscopic constants are computed for all isotopologues using second-order vibrational perturbation theory (VPT2), and the fundamental vibrational frequencies are computed with VPT2 and vibrational configuration interaction (VCI) theory. Agreement between VPT2 and VCI results is quite good with the fundamental vibrational frequencies of the bent isomer isotopologues in accord to within a 0.1 to 3.2 / cm range. Similar accuracies are present for the cyclic isomer isotopologues. The data generated here serve as a reference for astronomical observations of these closed-shell, highly-dipolar molecules using new, high-resolution telescopes and as reference for laboratory studies where isotopic labeling may lead to elucidation of the formation mechanism for the known interstellar molecule: X 3A0 HCCN.

The interstellar presence of protonated nitrous oxide has been suspected for some time. Using established high-accuracy quantum chemical techniques, spectroscopic constants and fundamental vibrational frequencies are provided for the lower energy O-protonated isomer of this cation and its deuterated isotopologue. The vibrationally-averaged B0 and C0 rotational constants are within 6 MHz of their experimental values and the D(subJ) quartic distortion constants agree with experiment to within 3%. The known gas phase O−H stretch of NNOH(+) is 3330.91 cm(exp−1), and the vibrational configuration interaction computed result is 3330.9 cm(exp−1). Other spectroscopic constants are also provided, as are the rest of the fundamental vibrational frequencies for NNOH(+) and its deuterated isotopologue. This high-accuracy data should serve to better inform future observational or experimental studies of the rovibrational bands of protonated nitrous oxide in the ISM and the laboratory.

The ion-molecule association mechanism of acetylene and ethylene with their cations is investigated by ab initio quantum chemical methods to understand the structures, association energies, and the vibrational and electronic spectra of the products. Stable puckered cyclic isomers are found as the result of first forming less stable linear and bridge isomers. The puckered cyclic complexes are calculated to be strongly bound, by 87, 35 and 56 kcal/mol for acetylene-acetylene cation, ethylene-ethylene cation and acetylene-ethylene cation, respectively. These stable complexes may be intermediates that participate in further association reactions. There are no association barriers, and no significant inter-conversion barriers, so the initial linear and bridge encounter complexes are unlikely to be observable. However, the energy gap between the bridged and cyclic puckered isomers greatly differs from complex to complex: it is 44 kcal/mol in C4H4 +, but only 6 kcal/mol in C4H8 +. The accurate CCSD(T) calculations summarized above are also compared against less computationally expensive MP2 and density functional theory (DFT) calculations for structures, relative energies, and vibrational spectra. Calculated vibrational spectra are compared against available experiments for cyclobutadiene cation. Electronic spectra are also calculated using time-dependent DFT.

It has been shown that rotational lines observed in the Horsehead nebula photon-dominated-region (PDR) are probably not caused by l-C3H+, as was originally suggested. In the search for viable alternative candidate carriers, quartic force fields are employed here to provide highly accurate rotational constants, as well as fundamental vibrational frequencies, for another candidate carrier: 1 (sup 1)A' C3H(−). The ab initio computed spectroscopic constants provided in this work are, compared to those necessary to define the observed lines, as accurate as the computed spectroscopic constants for many of the known interstellar anions. Additionally, the computed D-eff for C3H(−) is three times closer to the D deduced from the observed Horsehead nebula lines relative to l-C3H(+). As a result, 1 (sup 1)A' C3H(−). is a more viable candidate for these observed rotational transitions and would be the seventh confirmed interstellar anion detected within the past decade and the first C(sub n)H(−) molecular anion with an odd n.

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

The A 1B1 <-1A0 excitation into the dipole-bound state of the cyanomethyl anion (CH2CN􀀀) has been hypothesized as the carrier for one di use interstellar band. However, this particular molecular system has not been detected in the interstellar medium even though the related cyanomethyl radical and the isoelectronic ketenimine molecule have been found. In this study we are employing the use of proven quartic force elds and second-order vibrational perturbation theory to compute accurate spectroscopic constants and fundamental vibrational frequencies for ~X 1A0 CH2CN􀀀 in order to assist in laboratory studies and astronomical observations. Keywords: Astrochemistry, ISM: molecular anions, Quartic force elds, Rotational constants, Vibrational frequencies

Highly correlated ab initio quartic force field (QFFs) are used to calculate the equilibrium structures and predict the spectroscopic parameters of three HC2N isomers. Specifically, the ground state quasilinear triplet and the lowest cyclic and bent singlet isomers are included in the present study. Extensive treatment of correlation effects were included using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T). Dunning s correlation-consistent basis sets cc-pVXZ, X=3,4,5, were used, and a three-point formula for extrapolation to the one-particle basis set limit was used. Core-correlation and scalar relativistic corrections were also included to yield highly accurate QFFs. The QFFs were used together with second-order perturbation theory (with proper treatment of Fermi resonances) and variational methods to solve the nuclear Schr dinger equation. The quasilinear nature of the triplet isomer is problematic, and it is concluded that a QFF is not adequate to describe properly all of the fundamental vibrational frequencies and spectroscopic constants (though some constants not dependent on the bending motion are well reproduced by perturbation theory). On the other hand, this procedure (a QFF together with either perturbation theory or variational methods) leads to highly accurate fundamental vibrational frequencies and spectroscopic constants for the cyclic and bent singlet isomers of HC2N. All three isomers possess significant dipole moments, 3.05D, 3.06D, and 1.71D, for the quasilinear triplet, the cyclic singlet, and the bent singlet isomers, respectively. It is concluded that the spectroscopic constants determined for the cyclic and bent singlet isomers are the most accurate available, and it is hoped that these will be useful in the interpretation of high-resolution astronomical observations or laboratory experiments.

A molecule's absorption in the atmospheric infrared (IR) window (IRW) is an indicator of its efficiency as a greenhouse gas. A model for estimating the absorption of a fluorinated molecule within the IRW was developed to assess its radiative impact. This model will be useful in comparing different hydrofluorocarbons and hydrofluoroethers contribution to global warming. The absorption of radiation by greenhouse gases, in particular hydrofluoroethers and hydrofluorocarbons, was investigated using ab initio quantum mechanical methods. Least squares regression techniques were used to create a model based on this data. The placement and number of fluorines in the molecule were found to affect the absorption in the IR window and were incorporated into the model. Several group increment models are discussed. An additive model based on one-carbon groups is found to work satisfactorily in predicting the ab initio calculated vibrational intensities.

High levels of theory have been used to compute quartic force fields (QFFs) for the cyclic and linear forms of the C H + molecular cation, referred to as c-C H + and I-C H +. Specifically the 33 3333 singles and doubles coupled-cluster method that includes a perturbational estimate of connected triple excitations, CCSD(T), has been used in conjunction with extrapolation to the one-particle basis set limit and corrections for scalar relativity and core correlation have been included. The QFFs have been used to compute highly accurate fundamental vibrational frequencies and other spectroscopic constants using both vibrational 2nd-order perturbation theory and variational methods to solve the nuclear Schroedinger equation. Agreement between our best computed fundamental vibrational frequencies and recent infrared photodissociation experiments is reasonable for most bands, but there are a few exceptions. Possible sources for the discrepancies are discussed. We determine the energy difference between the cyclic and linear forms of C H +, 33 obtaining 27.9 kcal/mol at 0 K, which should be the most reliable available. It is expected that the fundamental vibrational frequencies and spectroscopic constants presented here for c-C H + 33 and I-C H + are the most reliable available for the free gas-phase species and it is hoped that 33 these will be useful in the assignment of future high-resolution laboratory experiments or astronomical observations.

A feasible initiation step for particle growth in the interstellar medium (ISM) is simulated by means of ab quantum chemistry methods. The systems studied are dimer ions formed by pairing nitrogen containing small molecules known to exist in the ISM with ions of unsaturated hydrocarbons or vice versa. Complexation energies, structures of ensuing complexes and electronic excitation spectra of the encounter complexes are estimated using various quantum chemistry methods. Moller-Plesset perturbation theory (MP2, Z-averaged perturbation theory (ZAP2), coupled cluster singles and doubles with perturbative triples corrections (CCSD(T)), and density functional theory (DFT) methods (B3LYP) were employed along with the correlation consistent cc-pVTZ and aug-cc-pVTZ basis sets. Two types of complexes are predicted. One type of complex has electrostatic binding with moderate (7-20 kcal per mol) binding energies, that are nonetheless significantly stronger than typical van der Waals interactions between molecules of this size. The other type of complex develops strong covalent bonds between the fragments. Cyclic isomers of the nitrogen containing complexes are produced very easily by ion-molecule reactions. Some of these complexes show intense ultraviolet visible spectra for electronic transitions with large oscillator strengths at the B3LYP, omegaB97, and equations of motion coupled cluster (EOM-CCSD) levels. The open shell nitrogen containing carbonaceous complexes especially exhibit a large oscillator strength electronic transition in the visible region of the electromagnetic spectrum.

We have investigated the physical characteristics of greenhouse gases (GHGs) to assess which properties are most important in determining the efficiency of a GHG. Chlorofluorcarbons (CFCs), hydrofluorocarbons (HFCs), perfluorocarbons (PFCs), nitrogen fluorides, and various other known atmospheric trace molecules have been included in this study. Compounds containing the halogens F or Cl have in common very polar X-F or X-Cl bonds, particularly the X-F bonds. It is shown that as more F atoms bond to the same central atom, the bond dipoles become larger as a result of the central atom becoming more positive. This leads to a linear increase in the total or integrated XF bond dipole derivatives for the molecule, which leads to a non-linear (quadratic) increase in infrared (IR) intensity. Moreover, virtually all of the X-F bond stretches occur in the atmospheric IR window as opposed to X-H stretches, which do not occur in the atmospheric window. It is concluded that molecules possessing several F atoms will always have a large radiative forcing parameter in the calculation of their global warming potential. Some of the implications for global warming and climate change are discussed.

This report will discuss the use of advanced simulation techniques to optimize the performance of the proposed Orion Crew Module airbag landing system design. The Boeing Company and the National Aeronautic and Space Administration s Langley Research Center collaborated in the analysis of the proposed airbag landing system for the next generation space shuttle replacement, the Orion spacecraft. Using LS-DYNA to simulate the Crew Module landing impacts, two main objectives were established and achieved: the investigation of potential methods of optimizing the airbag performance in order to reduce rebound on the anti-bottoming bags, lower overall landing loads, and increase overall Crew Module stability; and the determination of the Crew Module stability and load boundaries using the optimized airbag design, based on the potential Crew Module landing pitch angles and ground slopes in both the center of gravity forward and aft configurations. This paper describes the optimization and stability and load boundary studies and presents a summary of the results obtained and key lessons learned from this analysis.

A paper [Chem. Phys. Lett. 345, 295 (2001)] describes theoretical studies of excited electronic states of nitrogen molecules, with a view toward utilizing those states in synthesizing tetrahedral N4, or Td N4 a metastable substance under consideration as a high-energy-density rocket fuel. Several ab initio theoretical approaches were followed in these studies, including complete active space self-consistent field (CASSCF), state-averaged CASSCF (SA-CASSCF), singles configuration interaction (CIS), CIS with second-order and third-order correlation corrections [CIS(D) and CIS(3)], and linear response singles and doubles coupled-cluster (LRCCSD). Standard double zeta polarized and triple zeta double polarized one-particle basis sets were used. The CASSCF calculations overestimated the excitation energies, while SACASSCF calculations partly corrected these overestimates. The accuracy of the CIS calculations varied, depending on the particular state, while the CIS(D), CIS(3), and LRCCSD results were in generally good agreement. The energies of the lowest six excited singlet states of Td N4 as calculated by the LRCCSD were compared with the energies of possible excited states of N2 + N2 fragments, leading to the conclusion that the most likely route for synthesis of Td N4 would involve a combination of two bound quintet states of N2.