Search

The azirinyl cation is an aromatic cyclic molecule that is isoelectronic with cyclopropenylidene, c-C3H2, and c-C3H3+. Cyclopropenylidene has been shown to be ubiquitous, existing in many different astrophysical environments. Given the similar chemistry between C and N, and the relative abundances between C and N in astrophysical environments, it is expected that there should be aromatic ringed molecules that incorporate N in the ring, but as yet, no such molecule has been identified. To address this issue, the present study uses high levels of electronic structure theory to compute a highly accurate quartic force field (QFF) for the azirinyl cation and its two lowest lying isomers, the cyanomethyl and isocyanomethyl cations. The theoretical approach uses the singles and doubles coupled-cluster method that includes a perturbative correction for connected triple excitations, CCSD(T), together with extrapolation to the one-particle basis set limit and corrections for scalar relativity and core-correlation. The QFF is then used in a second-order vibrational perturbation theory analysis (VPT2) to compute the fundamental vibrational frequencies and rovibrational spectroscopic constants for all three C2NH2+ isomers. The reliability of the VPT2 vibrational frequencies is tested by comparison to vibrational configuration interaction (VCI) calculations, and excellent agreement is found between the two approaches. Fundamental vibrational frequencies and rovibrational spectroscopic constants for all singly substituted 13C, 15N, and D isotopologues are also reported. It is expected that the highly accurate spectroscopic data reported herein will be useful in the identification of these cations in high-resolution experimental or astronomical observations.

The simple aromatic hydrocarbon, cyclopropenylidene (c-C3H2), is a known, naturally-occurring molecule. The question remains as to whether its isoelectronic, cyclic, fellow aromatics of c-N3+, c-CNN, HCNN+, and c-CNC− are as well. Each of these are exciting objects for observation of Titan, and the rotational constants and vibrational frequencies produced here will allow for remote sensing of Titan's atmosphere or other astrophysical or terrestrial sources. None of these four aromatic species are vibrationally strong absorbers/emitters, but the two ions, HCNN+ and c-CNC−, have dipole moments of greater than 3 D and 1 D, respectively, making them good targets for rotational spectroscopic observation. Each of these molecules is shown here to exhibit its own, unique vibrational properties, but the general trends put the vibrational behavior for corresponding fundamental modes within close ranges of one another, even producing nearly the same heavy atom, symmetric stretching frequencies for HCNN+ and c-C3H2 at 1600 cm−1. The c-N3+ cation is confirmed to be fairly unstable and has almost no intensity in its ν2 fundamental. Hence, it will likely remain difficult to characterize experimentally.

The vibrational spectra for the HOCO radical in both of its conformers and deuterated isotopologues are shown here for the first time. Building on previous work with coupled cluster quartic force fields (QFFs) in the computation of the fundamental vibrational frequencies for both cis- and trans-HOCO, coupled cluster dipole surfaces are now provided for both HOCO conformers and their corresponding deuterated isotopologues. These surfaces and subsequent vibrational configuration interaction (VCI) computations produce the intensities of transitions into vibrational states including the fundamentals, overtones, and first few combination bands of less than 4000 cm(-1), slightly beyond the O-H stretch. Simulated spectra with an artificial full width at half-maximum broadening of 10 cm(-1) are also provided in order to aid in the characterization of HOCO's vibrational frequencies and to assist detection in various laboratory or astronomical observations.

Building on previous studies involving coupled cluster quartic force fields for the description of spectroscopic constants and vibrational frequencies of astronomically relevant molecules, this work applies the same techniques to the elucidation of such properties for the bent 1 (3)A' state of HCN and the isoelectronic 1 (3)A' HCO(+). Core correlation is treated both by explicit means and as a correction. Each approach gives closely comparable spectroscopic constants and vibrational frequencies once more, indicating that the composite method is a viable and less costly alternative. We are providing fundamental vibrational frequencies for these systems where agreement with experiment in previous studies has been within 4 cm(-1) or better. Frequencies for the first overtones and combination bands as well as various spectroscopic constants are also reported.

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

Besides the ν(1) O-H stretching mode at 3435 cm(-1) for HOCS(+), the fundamental vibrational frequencies for this cation and its HSCO(+) isomer have not been determined experimentally. Because these systems are analogues to HOCO(+), a detected interstellar molecule, and are believed to play an important role in reactions of OCS, which has also been detected in the interstellar medium, these cations are of importance to interstellar chemistry and reaction surface studies. This work provides the fundamental vibrational frequencies and spectroscopic constants computed with vibrational perturbation theory (VPT) at second order and the vibrational configuration interaction (VCI) method conjoined with the most accurate quartic force field (QFF) applied to date for these systems. Our computations match experiment to better than 2 cm(-1) for the known O-H stretch. Additionally, there is strong agreement in the prediction of the fundamentals across methods and choices of QFFs. The consistency in the computations and the correspondence for the known mode should give accurate reference data for the rovibrational spectra of these cations and their singly substituted isotopologues for D, (18)O, and (34)S.

The Born−Oppenheimer potential energy surface(s) underlies theoretical and computational chemistry (whether one considers a single or multiply coupled surfaces). The recent progress in representing these surfaces, rigorously obtained from electronic structure calculations, is the focus of this Perspective. Examples of potentials of complex molecules, namely, CH3CHO, CH5+, and H5+, and molecular complexes, namely, water clusters, are given.

Aims. We aim to elucidate the spectral changes in the 3 micron region that result from chemical changes in the molecular periphery of polycyclic aromatic hydrocarbons (PAHs) with extra hydrogens (H-PAHs) and methyl groups (Me-PAHs). Methods. Advanced laser spectroscopic techniques combined with mass spectrometry were applied on supersonically cooled 1,2,3,4-tetrahydronaphthalene, 9,10-dihydroanthracene, 9,10-dihydrophenathrene, 1,2,3,6,7,8-hexahydropyrene, 9-methylanthracene, and 9,10-dimethylanthracene, allowing us to record mass-selective and conformationally selective absorption spectra of the aromatic, aliphatic, and alkyl CH-stretches in the 3.175-3.636 micron region with laser-limited resolution. We compared the experimental absorption spectra with standard harmonic calculations and with second-order vibrational perturbation theory anharmonic calculations that use the SPECTRO program for treating resonances. Results. We show that anharmonicity plays an important if not dominant role, affecting not only aromatic, but also aliphatic and alkyl CH-stretch vibrations. The experimental high-resolution data lead to the conclusion that the variation in Me- and H-PAHs composition might well account for the observed variations in the 3 micron emission spectra of carbon-rich and star-forming regions. Our laboratory studies also suggest that heavily hydrogenated PAHs form a significant fraction of the carriers of IR emission in regions in which an anomalously strong 3 micron plateau is observed., Comment: 13 pages, 8 Figures, 1 table, accepted for publication in A&A

The effect of approximating the three- and four-virtual molecular orbital integrals in single and double coupled-cluster theory including a perturbational correction for connected triple excitations [CCSD(T)] is investigated for the calculation of higher-order properties, specifically the calculation of a molecular quartic force field and spectroscopic constants. The approach was proposed previously, but investigated for only second- and lower-order properties. It is shown that the conclusions reached previously are essentially unchanged on moving to higher-order properties. That is, approximating the selected integrals has essentially no effect on the accuracy of CCSD(T) calculations, and the error due to approximating integrals is much smaller than the residual error due to one-particle basis set deficiencies. The advantage of this approach is that it significantly reduces the amount of data needed to perform CCSD(T) calculations, thereby reducing computational requirements associated with input/output ...

We present a comparison of the infra-red spectra of Cl-(H2O) and H+(H2O)(2) obtained with classical and ring-polymer molecular dynamics with previous quantum calculations. Full-dimensional ab initio-based potential and dipole-moment surfaces are used in these calculations. (C) 2007 Elsevier B.V. All rights reserved.

Current efforts to characterize and study interstellar polycyclic aromatic hydrocarbons(PAHs) rely heavily on theoretically predicted infrared (IR) spectra. Generally, such studies use the scaled harmonic frequencies for band positions and double harmonic approximation for intensities of species, and then compare these calculated spectra with experimental spectra obtained under matrix isolation conditions. High-resolution gas-phase experimental spectroscopic studies have recently revealed that the double harmonic approximation is not sufficient for reliable spectra prediction. In this paper, we present the anharmonic theoretical spectra of three PAHs: naphthalene, anthracene, and tetracene, computed with a locally modified version of the SPECTRO program using Cartesian derivatives transformed from Gaussian 09 normal coordinate force constants. Proper treatments of Fermi resonances lead to an impressive improvement on the agreement between the observed and theoretical spectra, especially in the C–H stretching region. All major IR absorption features in the full-scale matrix-isolated spectra, the high-temperature gas-phase spectra, and the most recent high-resolution gas-phase spectra obtained under supersonically cooled molecular beam conditions in the CH-stretching region are assigned.

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

Extensions to the code MULTIMODE to obtain rovibrational wave functions and properties are described. An application of these new capabilities is made to a calculation of the Franck–Condon factors for photoionization of CF3 to CF . These calculations make use of a new, full-dimensional ab initio potential energy surface, which is also described here.

We present two quantum calculations of the infrared spectrum of protonated methane (CH 5 + ) using full-dimensional, ab initio–based potential energy and dipole moment surfaces. The calculated spectra compare well with a low-resolution experimental spectrum except below 1000 cm –1 , where the experimental spectrum shows no absorption. The present calculations find substantial absorption features below 1000 cm –1 , in qualitative agreement with earlier classical calculations of the spectrum. The major spectral bands are analyzed in terms of the molecular motions. Of particular interest is an intense feature at 200 cm –1 , which is due to an isomerization mode that connects two equivalent minima. Very recent high-resolution jet-cooled spectra in the CH stretch region (2825 to 3050 cm –1 ) are also reported, and assignments of the band origins are made, based on the present quantum calculations.

A full-dimensional ab initio potential energy surface (PES) and dipole moment surface (DMS) are reported for the water dimer, (H2O)2. The CCSD(T)-PES is a very precise fit to 19,805 ab initio energies obtained with the coupled-cluster (CCSD(T)) method, using an aug-cc-pVTZ basis. The standard counterpoise correction was applied to approximately eliminate basis set superposition errors. The fit is based on an approach that incorporates the permutational symmetry of identical atoms [Huang, X.; Braams, B.; Bowman, J. M. J. Chem.Phys. 2005, 122, 044308]. The DMS is a fit to the dipole moment obtained with Møller-Plesset (MP2) theory, using an aug-cc-pVTZ basis. The PES has an RMS fitting error of 31 cm(-1) for energies below 20,000 cm(-1), relative to the global minimum. This surface can describe various internal floppy motions, including various monomer inversions, and isomerization pathways. Ten characteristic stationary points have been located on the surface, four of which are transition structures and the rest are higher order saddle points. Their geometrical and vibrational properties are presented and compared with available previous theoretical work. The CCSD(T)-PES and MP2-DMS dissociate correctly (and symmetrically) to two H2O monomers, with D(e) = 1665.7 cm(-1) (19.93 kJ/mol). Accurate quantum calculations of the zero-point energy of the dimer (using diffusion Monte Carlo) and the monomers (using a vibrational configuration interaction approach) are reported, and these together with D(e) give a value of D0 of 1042 cm(-1) (12.44 kJ/mol). A best estimated value is 1130 cm(-1) (13.5 kJ/mol).

We report an ab initio calculation of the potential surface, quantum structures, and zero-point energies of CH5+ and CH2D3+ in full dimensionality. This potential energy surface is a very precise f...

The vibrational mode coupling in the protonated water dimer is investigated by performing two types of quantum calculations of the vibrational levels of H5O2+ and D5O2+, utilizing the OSS3(p) potential energy surface by Ojamae et al. [L. Ojamae, I. Shavitt, and S. J. Singer, J. Chem. Phys. 109, 5547 (1998)]. One is four-dimensional (4D), treating only the central O⋯H(D)+⋯O moiety. Three of the four modes considered, the asymmetric stretch and the two bends, are largely the vibrations of the central proton, while the fourth mode is essentially the O⋯O stretching vibration. The vibrational levels of O⋯H(D)+⋯O are calculated rigorously, as fully coupled (FC), and also in an adiabatic (3+1)D approximation, where the proton asymmetric stretch is treated as adiabatically separated from the other three degrees of freedom. The second set of calculations, designated VCI, is full-dimensional, 15D; it is performed by the code MULTIMODE, which does configuration interaction (CI) calculations using a basis determined ...

This review focuses on the calculation of rovibrational energies of polyatomic molecules using the code MULTIMODE. This code, which uses normal coordinates and a hierarchical n-mode representation of the potential, aims to be applicable to a wide class of molecules and molecular complexes. The theoretical and computational methods used in this code are described, followed by a review of selected applications. These applications illustrate various features of the code and also point out some limitations of the current version of the code. The review concludes with some ideas about possible future directions in this area of research.

A new potential energy surface, based on high quality ab initio electronic structure calculations, is presented for the hydronium ion (H3O+). The new potential surface is used in rigorous calculations of vibrational energies of H3O+, D3O+, H2DO+, and HD2O+. Comparison with experiment shows significant improvement over our previous calculations using an earlier potential [X. Huang, S. C. Carter, and J. M. Bowman, J. Phys. Chem. B 106, 8182 (2002)]. Vibrational calculations are also presented with a new version of the code MULTIMODE. In this version the maximum number of coupled modes in the potential in any grouping of modes is increased from four (the previous maximum) to five. The importance of five-mode terms in the potential is demonstrated for several vibrational states in H3O+ and H2DO+. Also, in the new version of MULTIMODE the number of coupled modes in the Coriolis term can be varied independently from the number of coupled modes in the potential. Rovibrational calculations for J=1 are also presen...

We report a new full dimensional potential energy surface for the hydronium ion (H 3 O + ). This surface is constructed by a least-squares fit of ab initio electronic energies obtained using the CCSD(T) method with an aug-cc-pVTZ basis set, augmented by some calculations using an aug-cc-pVQZ basis set. The calculated inversion barrier is 693 cm - 1 . Full dimensional vibrational calculations are reported using the new surface for H 3 O + , D 3 O + , H 2 DO + , and HD 2 O + using the codes RVIB4 and MULTIMODE. Comparison with all available experimental data on both splittings and vibrational energies shows significant improvement over our previous calculations using the OSS3(p) potential. A number of vibrational band energies and splittings not measured experimentally are reported for D 3 O + , HD 2 O + , and H 2 DO + .

Protonated acetylene, C2H3(+), is among the simplest carbocations. Comprehensive experimental or highly accurate computational spectroscopic data is lacking for this system due to its inherent complexities. Utilizing state-of-the-art quartic force fields (QFFs), the spectroscopic constants and fundamental vibrational frequencies are provided in this work for the nonclassical, bridged, cyclic global minimum. The rotational constants match experiment to better than 0.1%, and the computed ν2 antisymmetric HCCH stretch is less than 3.0 cm(–1) different from experiment. Hence, the rovibrational spectroscopic data provided herein for c-C2H3(+) and its deuterated isotopologues enrich the chemical understanding of this system. Unfortunately, the same rovibrational spectroscopic data is not as trustworthy for the classical, linear form of protonated acetylene due to the shallow well in which it resides on the potential energy surface. However, spectroscopic data are provided for this isomer in the Supporting Information to enhance future studies.

Highly accurate, coupled-cluster-based quartic force fields (QFFs) have been employed recently to provide spectroscopic reference for a myriad of molecules. Here, we are extending the same approach to provide vibrational and rotational spectroscopic reference data for the sulfur analogues of HOCO, HSCO, and HOCS, in both the cis and trans conformations as well as the D and (34)S isotopologues of each system. The resulting energies corroborate previous computations showing that trans-HSCO is the lowest-energy isomer for this system. The vibrational frequencies are computed with both second-order vibrational perturbation theory (VPT2) and vibrational configuration interaction (VCI) methods. The VPT2 and VCI QFF frequencies largely agree with one another to better than 5.0 cm(-1) (often better than 1.0 cm(-1)) and are also consistent with the type of behavior exhibited in previous studies. As such, the reference data provided here should assist in analysis of environments in which these sulfur systems may be found, including the interstellar medium, combustion flames, or laboratory simulations of either.

Contains fulltext : 166299.pdf (Author’s version preprint ) (Open Access) 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 μ m absorption band. For 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 μ 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 data set 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 μ m band, and on features such as the two-component emission character of this band and the 3 μ m emission plateau. 11 p.

In this work, computational procedures are employed to compute the rotational and rovibrational spectra and line lists for H2O, CO2, and SO2. Building on the established use of quartic force fields, MP2 and CCSD(T) Dipole Moment Surfaces (DMSs) are computed for each system of study in order to produce line intensities as well as the transition energies. The computed results exhibit a clear correlation to reference data available in the HITRAN database. Additionally, even though CCSD(T) DMSs produce more accurate intensities as compared to experiment, the use of MP2 DMSs results in reliable line lists that are still comparable to experiment. The use of the less computationally costly MP2 method is beneficial in the study of larger systems where use of CCSD(T) would be more costly.

The study of interstellar polycyclic aromatic hydrocarbons (PAHs) relies heavily on theoretically predicted infrared spectra. Most earlier studies use scaled harmonic frequencies for band positions and the double harmonic approximation for intensities. However, recent high-resolution gas-phase experimental spectroscopic studies have shown that the harmonic approximation is not sufficient to reproduce experimental results. In our previous work, we presented the anharmonic theoretical spectra of three linear PAHs, showing the importance of including anharmonicities into the theoretical calculations. In this paper, we continue this work by extending the study to include five non-linear PAHs (benz[a]anthracene, chrysene, phenanthrene, pyrene, and triphenylene), thereby allowing us to make a full assessment of how edge structure, symmetry, and size influence the effects of anharmonicities. The theoretical anharmonic spectra are compared to spectra obtained under matrix isolation low-temperature conditions, low-resolution, high-temperature gas-phase conditions, and high-resolution, low-temperature gas-phase conditions. Overall, excellent agreement is observed between the theoretical and experimental spectra although the experimental spectra show subtle but significant differences.

An isotopic-independent, highly accurate potential energy surface (PES) has been determined for CO(2) by refining a purely ab initio PES with selected, purely experimentally determined rovibrational energy levels. The purely ab initio PES is denoted Ames-0, while the refined PES is denoted Ames-1. Detailed tests are performed to demonstrate the spectroscopic accuracy of the Ames-1 PES. It is shown that Ames-1 yields σ(rms) (root-mean-squares error) = 0.0156 cm(-1) for 6873 J = 0-117 (12)C(16)O(2) experimental energy levels, even though less than 500 (12)C(16)O(2) energy levels were included in the refinement procedure. It is also demonstrated that, without any additional refinement, Ames-1 yields very good agreement for isotopologues. Specifically, for the (12)C(16)O(2) and (13)C(16)O(2) isotopologues, spectroscopic constants G(v) computed from Ames-1 are within ±0.01 and 0.02 cm(-1) of reliable experimentally derived values, while for the (16)O(12)C(18)O, (16)O(12)C(17)O, (16)O(13)C(18)O, (16)O(13)C(17)O, (12)C(18)O(2), (17)O(12)C(18)O, (12)C(17)O(2), (13)C(18)O(2), (13)C(17)O(2), (17)O(13)C(18)O, and (14)C(16)O(2) isotopologues, the differences are between ±0.10 and 0.15 cm(-1). To our knowledge, this is the first time a polyatomic PES has been refined using such high J values, and this has led to new challenges in the refinement procedure. An initial high quality, purely ab initio dipole moment surface (DMS) is constructed and used to generate a 296 K line list. For most bands, experimental IR intensities are well reproduced for (12)C(16)O(2) using Ames-1 and the DMS. For more than 80% of the bands, the experimental intensities are reproduced with σ(rms)(ΔI) < 20% or σ(rms)(ΔI∕δ(obs)) < 5. A few exceptions are analyzed and discussed. Directions for future improvements are discussed, though it is concluded that the current Ames-1 and the DMS should be useful in analyzing and assigning high-resolution laboratory or astronomical spectra.

Beginning with the seminal paper of Born and Oppenheimer (BO) [1] in 1927, the concept of the potential energy surface (PES) plays a critical role in the description, simulation, and modeling of molecular systems. It provides the basis [2] for understanding the processes associated with the nuclear motions in molecules. By going beyond the characteristic stationary points and barriers, full dimensional, accurate potential energy surfaces have a very broad range of applications in many areas of physical chemistry; for example, they provide insight into structure, reactivity, and spectroscopy of molecules.While the majority of stable structures on PES are associated with the covalent or ionic bonding [3], the regions of PES dominated by van derWaals interactions are essential for low-temperature phenomena and molecular stacking which are critical for understanding biomolecular structures involving DNA and RNA molecules [4]. Furthermore, the advances in “cold chemistry” [5] make it possible to test the theoretical predictions (see, e.g., [6]) involving very small barriers of

The use of accurate quartic force fields together with vibrational configuration interaction recently predicted gas phase fundamental vibrational frequencies of the trans-HOCO radical to within 4 cm(-1) of experimental results for the two highest frequency modes. Utilizing the same approach, we are providing a full list of fundamental vibrational frequencies and spectroscopic constants for the cis-HOCO system in both radical and anionic forms. Our predicted geometrical parameters of the cis-HOCO radical match experiment and previous computation to better than 1% deviation, and previous theoretical work agrees equally well for the anion. Correspondence between vibrational perturbation theory and variational vibrational configuration interaction for prediction of the frequencies of each mode is strong, better than 5 cm(-1), except for the torsional motion, similar to what has been previously identified in the trans-HOCO radical. Among other considerations, our results are immediately applicable to dissociative photodetachment experiments which initially draw on the cis-HOCO anion since it is the most stable conformer of the anion and is used to gain insight into the portion of the OH + CO potential surface where the HOCO radical is believed to form, and we are also providing highly accurate electron binding energies relevant to these experiments.

In the search for a full mechanism creating CO(2) from OH + CO, it has been suggested that creation of the hydroxyformyl or HOCO radical may be a necessary step. This reaction and its transient intermediate may also be responsible for the regeneration of CO(2) in such high quantities in the atmosphere of Mars. Past spectroscopic observations of this radical have been limited and a full gas phase set of the fundamental vibrational frequencies of the HOCO radical has not been reported. Using established, highly accurate quantum chemical coupled cluster techniques and quartic force fields, we are able to compute all six fundamental vibrational frequencies and other spectroscopic constants for trans-HOCO in the gas phase. These methods have yielded rotational constants that are within 0.01 cm(-1) for A(0) and 10(-4) cm(-1) for B(0) and C(0) compared with experiment as well as fundamental vibrational frequencies within 4 cm(-1) of the known gas phase experimental ν(1) and ν(2) modes. Such results lead us to conclude that our prediction of the other four fundamental modes of trans-HOCO are also quite reliable for comparison to future experimental observation, though the discrepancy for the torsional mode may be larger since it is fairly anharmonic. With the upcoming European Space Agency/NASA ExoMars Trace Gas Orbiter, these data may help to establish whether HOCO is present in the Martian sky and what role it may play in the retention of a CO(2)-rich atmosphere. Furthermore, these data may also help to clear up questions built around the fundamental chemical process of how exactly the OH + CO reaction progresses.

We report full-dimensional, ab initio potential energy and dipole moment surfaces, denoted PES and DMS, respectively, for arbitrary numbers of water monomers. The PES is a sum of 1-, 2-, and 3-body potentials which can also be augmented by semiempirical long-range higher-body interactions. The 1-body potential is a spectroscopically accurate monomer potential, and the 2- and 3-body potentials are permutationally invariant fits to tens of thousands of CCSD(T)/aug-cc-pVTZ and MP2/aug-cc-pVTZ electronic energies, respectively. The DMS is a sum of 1- and 2-body DMS, which are covariant fits to tens of thousands MP2/aug-cc-pVTZ dipole moment data. We present the details of these new 2- and 3-body potentials and then extensive applications and tests of this PES are made to the structures, classical binding energies, and harmonic frequencies of water clusters up to the 22-mer. In addition, we report the dipole moment for these clusters at various minima and compare the results against available and new ab initio calculations.

A new approach is proposed and investigated for approximately including the effects of diffuse functions in one-particle basis sets when high accuracy is desired. The method is cost-effective for use in computing quartic force fields (QFFs), global potential energy surfaces (PESs), or other situations when a large part of the PES is needed. It is conservatively estimated that the use of this approximation leads to a computational savings of a factor of five, and it is argued that this could be significantly larger if input/output wait times are considered. It can be used when extrapolation to the one-particle basis set limit is performed, or it can be used simply to approximate the effect of diffuse functions for a larger basis set. The new approach is based on scaling the diffuse function effect for a smaller basis set to approximate the effect for a larger basis or an extrapolated energy in which larger basis set(s) are used. The scale factor is written as a function of the geometrical coordinates of the molecule and thus it includes a geometry dependence. We report results where the scale factor is a constant, includes through gradient terms, includes through second derivative terms, and includes through diagonal second and third derivative terms. The method has been tested in the calculation of accurate QFFs, equilibrium structures, and harmonic and fundamental vibrational frequencies for NH(2)(-), OH(-), H(2)O, and CH(3)OH. It is found that including up through diagonal second derivative terms leads to reliable fundamental vibrational frequencies and is cost-effective. It is also concluded that the use of a 5Z-quality basis set is essential if high accuracy is desired for these properties, even with extrapolation to the one-particle basis set limit.

A procedure for the calculation of molecular quartic force fields (QFFs) is proposed and investigated. The goal is to generate highly accurate ab initio QFFs that include many of the so-called "small" effects that are necessary to achieve high accuracy. The small effects investigated in the present study include correlation of the core electrons (core correlation), extrapolation to the one-particle basis set limit, correction for scalar relativistic contributions, correction for higher-order correlation effects, and inclusion of diffuse functions in the one-particle basis set. The procedure is flexible enough to allow for some effects to be computed directly, while others may be added as corrections. A single grid of points is used and is centered about an initial reference geometry that is designed to be as close as possible to the final ab initio equilibrium structure (with all effects included). It is shown that the least-squares fit of the QFF is not compromised by the added corrections, and the balance between elimination of contamination from higher-order force constants while retaining energy differences large enough to yield meaningful quartic force constants is essentially unchanged from the standard procedures we have used for many years. The initial QFF determined from the least-squares fit is transformed to the exact minimum in order to eliminate gradient terms and allow for the use of second-order perturbation theory for evaluation of spectroscopic constants. It is shown that this step has essentially no effect on the quality of the QFF largely because the initial reference structure is, by design, very close to the final ab initio equilibrium structure. The procedure is used to compute an accurate, purely ab initio QFF for the H(2)O molecule, which is used as a benchmark test case. The procedure is then applied to the ground and first excited electronic states of the HO(2) (+) molecular cation. Fundamental vibrational frequencies and spectroscopic constants from these highly accurate QFFs are compared and contrasted with previous experiment and theory. It is concluded that the spectroscopic constants determined for the X(3)A(") and A(1)A(') states of HO(2) (+) are the most accurately available.

An accurate potential energy surface has been determined for methanol from ab initio potential data at the CCSD(T) level of theory with an aug-cc-pVTZ basis. The resulting potential function is valid over all twelve vibrational degrees of freedom for all near-equilibrium and torsional configurations. A torsional reaction path has been derived for this potential, from which the low-lying vibrational levels of methanol have been calculated by the reaction path version of MULTIMODE. Comparisons with experiment and other calculations are made.

We present argon predissociation vibrational spectra of the OH(-).H(2)O and Cl(-).H(2)O complexes in the 1000-1900 cm(-1) energy range, far below the OH stretching region reported in previous studies. This extension allows us to explore the fundamental transitions of the intramolecular bending vibrations associated with the water molecule, as well as that of the shared proton inferred from previous assignments of overtones in the higher energy region. Although the water bending fundamental in the Cl(-).H(2)O spectrum is in very good agreement with expectations, the OH(-).H(2)O spectrum is quite different than anticipated, being dominated by a strong feature at 1090 cm(-1). New full-dimensionality calculations of the OH(-).H(2)O vibrational level structure using diffusion Monte Carlo and the VSCF/CI methods indicate this band arises from excitation of the shared proton.

The CH5+ molecular ion is well-known for its large amplitude motions that lead to complete scrambling of the hydrogen atoms, even in the vibrational ground state. Experiments have been reported that probe the consequences of these large amplitude motions. We recently reported that quantum zero-point effects partially quench the scrambling when CH5+ is partially deuterated. Here, the consequences of this quantum localization are investigated through calculations of the low-resolution spectra of CH4D+, CHD4+, and CD5+. The spectra are obtained by convoluting stick spectra, evaluated for individual stationary points on an ab initio potential surface, multiplying them by Diffusion Monte Carlo ground state density at that stationary point, and taking the sum. The CH/D stretch regions of CH4D+ and CD5+ are red-shifted relative to CH5+, while the overall shape of the envelope remains unaffected. In contrast, for CHD4+, the shape of the spectral envelope in the CH/D stretch region differs from the other three isotopologs. These signatures of the quantum localization of the deuterium on the spectra are discussed.

The vibrations of H3O2- and D3O2- are investigated using diffusion Monte Carlo (DMC) and vibrational configuration-interaction approaches, as implemented in the program MULTIMODE. These studies use the potential surface recently developed by Huang et al. [ J. Am. Chem. Soc. 126, 5042 (2004)]. The focus of this work is on the vibrational ground state and fundamentals which occur between 100 and 3700 cm(-1). In most cases, excellent agreement is obtained between the fundamental frequencies calculated by the two approaches. This serves to demonstrate the power of both methods for treating this very anharmonic system. Based on the results of the MULTIMODE and DMC treatments, the extent and nature of the couplings in H3O2- and D3O2- are investigated.

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.

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 configuration-interaction (VCI) approach. Agreement is within 10 cm(-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 C2v-symmetry, single (13)C-substituted isotopologues are presented, which may help identification of cyclic C4 in future experimental analyses or astronomical observations.

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(J) quartic distortion constants agree with experiment to within 3%. The known gas phase O-H stretch of NNOH(+) is 3330.91 cm(-1), and the vibrational configuration interaction computed result is 3330.9 cm(-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 interstellar medium and the laboratory.

We report a full-dimensional potential energy surface for H3O2-, based on fitting 66,965 ab initio electronic energies. A major feature of this potential is a barrier of roughly 200 cm-1 to internal rotation of the two hydroxyl groups about a line connecting the two oxygen atoms and the bridging hydrogen atom. The potential is used in calculations of vibrational energies, performed with the "Reaction Path" version of the code "MULTIMODE". The results are compared to recent infrared messenger experiments and are used to propose interpretations of the experimental results.

Several aspects of ammonia rovibrational spectra have been investigated using the new HSL-2 potential energy surface that includes an approximate correction for nonadiabatic effects. The unprecedented accuracy of rovibrational energy levels and transition energies computed using HSL-2 was demonstrated in Part I of this study. For (14)NH(3), new assignments for a few ν(3) + ν(4) band transitions and energy levels are suggested, and discrepancies between computed and HITRAN energy levels in the 2ν(4) band are analyzed (2ν(4) is the most difficult band below 5000 cm(-1)). New assignments are suggested for existing or missing 2ν(4) levels. Several new vibrational bands are identified from existing, unassigned HITRAN data, including 2ν(2) + ν(4), (ν(3) + ν(4)) -A(')∕A("), ν(1) + 2ν(2), and 2ν(2) + 2ν(4). The strong mixing between the 2ν(4) and 2ν(2) + ν(4) bands is carefully examined and found to be the source of the difficulties in the experimental modeling of 2ν(4). Discussion is presented for preliminary J = 10 results, where the overall root-mean-square error is estimated to be less than 0.039 cm(-1). The analysis of the 4ν(2) band demonstrates both the reliability and the accuracy of predictions from HSL-2. The full list of computed J = 0 band origins (with assignments) and the inversion splittings up to 7000-8000 cm(-1) above the zero-point energy are presented. J = 0-2 levels are reported for those bands below 5100 cm(-1) that are missing from the HITRAN database. For (15)NH(3), excellent agreement is found for the available ν(2) and ν(3) + ν(4)(E) transition energies, but significant deficiencies are shown for HITRAN levels and several corrections are suggested. The (15)N isotopic effects are presented for the J = 0-6 levels of 13 HITRAN bands. For (14)ND(3), we reproduce the pure rotational inversion spectra line frequencies with an accuracy similar to that for (14)NH(3). However, it is not possible to reproduce simultaneously all four pairs of inversion-split vibrational fundamentals to better than 0.05 cm(-1) uncertainty. It is suggested that a reanalysis of some suspicious (14)ND(3) fundamental bands is required. The analyses presented here and in Part I show that rovibrational energy levels and transition frequencies computed with HSL-2 (with nonadiabatic corrections) remain highly accurate well beyond the experimental data used in the refinement procedure. Calculations using HSL-2 are capable of revealing many deficiencies in experimental analyses of ammonia spectra and provide reliable predictions with similar accuracy. It is expected that the results of this study will be useful in the future interpretation of high-resolution spectra from laboratory experiments or from astronomical observations. The present work represents a very significant advance in the state of our knowledge of the spectroscopy of ammonia and its isotopologues.

In this work, we build upon our previous work on the theoretical spectroscopy of ammonia, NH(3). Compared to our 2008 study, we include more physics in our rovibrational calculations and more experimental data in the refinement procedure, and these enable us to produce a potential energy surface (PES) of unprecedented accuracy. We call this the HSL-2 PES. The additional physics we include is a second-order correction for the breakdown of the Born-Oppenheimer approximation, and we find it to be critical for improved results. By including experimental data for higher rotational levels in the refinement procedure, we were able to greatly reduce our systematic errors for the rotational dependence of our predictions. These additions together lead to a significantly improved total angular momentum (J) dependence in our computed rovibrational energies. The root-mean-square error between our predictions using the HSL-2 PES and the reliable energy levels from the HITRAN database for J = 0-6 and J = 7∕8 for (14)NH(3) is only 0.015 cm(-1) and 0.020∕0.023 cm(-1), respectively. The root-mean-square errors for the characteristic inversion splittings are approximately 1∕3 smaller than those for energy levels. The root-mean-square error for the 6002 J = 0-8 transition energies is 0.020 cm(-1). Overall, for J = 0-8, the spectroscopic data computed with HSL-2 is roughly an order of magnitude more accurate relative to our previous best ammonia PES (denoted HSL-1). These impressive numbers are eclipsed only by the root-mean-square error between our predictions for purely rotational transition energies of (15)NH(3) and the highly accurate Cologne database (CDMS): 0.00034 cm(-1) (10 MHz), in other words, 2 orders of magnitude smaller. In addition, we identify a deficiency in the (15)NH(3) energy levels determined from a model of the experimental data.

The full dimensional rate constant reported by Huarte-Larranaga and Manthe for the H+CH4→H2+CH3 reaction [Huarte-Larranaga and Manthe, J. Chem. Phys. 113, 5115 (2000)] is corrected by using an accurate vibrational partition function for CH4 instead of the harmonic normal-mode one used by them. The correction is shown to be substantial over the temperature range considered by Huarte-Larranaga and Manthe.

One-particle basis set extrapolation is compared with one of the new R12 methods for computing highly accurate quartic force fields (QFFs) and spectroscopic data, including molecular structures, rotational constants, and vibrational frequencies for the H2O, N2H+, NO2+, and C2H2 molecules. In general, agreement between the spectroscopic data computed from the best R12 and basis set extrapolation methods is very good with the exception of a few parameters for N2H+ where it is concluded that basis set extrapolation is still preferred. The differences for H2O and NO2+ are small and it is concluded that the QFFs from both approaches are more or less equivalent in accuracy. For C2H2, however, a known one-particle basis set deficiency for C–C multiple bonds significantly degrades the quality of results obtained from basis set extrapolation and in this case the R12 approach is clearly preferred over one-particle basis set extrapolation. The R12 approach used in the present study was modified in order to obtain hi...

We report the properties of two novel transition states of the bimolecular hydrogen exchange reaction in the water dimer, based on an ab initio water dimer potential [A. Shank et al., J. Chem. Phys. 130, 144314 (2009)]. The realism of the two transition states is assessed by comparing structures, energies, and harmonic frequencies obtained from the potential energy surface and new high-level ab initio calculations. The rate constant for the exchange is obtained using conventional transition state theory with a tunneling correction. We employ a one-dimensional approach for the tunneling calculations using a relaxed potential from the full-dimensional potential in the imaginary-frequency normal mode of the saddle point, Q(im). The accuracy of this one-dimensional approach has been shown for the ground-state tunneling splittings for H and D-transfer in malonaldehyde and for the D+H(2) reaction [Y. Wang and J. M. Bowman, J. Chem. Phys. 129, 121103 (2008)]. This approach is applied to calculate the rate constant for the H(2)O+H(2)O exchange and also for H(2)O+D(2)O→2HOD. The local zero-point energy is also obtained using diffusion Monte Carlo calculations in the space of real-frequency-saddle-point normal modes, as a function of Q(im).

A series of high-quality, purely ab initio, quartic force fields (QFFs), computed using a procedure we recently proposed, is reported for NH2− and CCH−. The singles and doubles coupled-cluster method with a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T), was used with TZ, QZ, and 5Z quality basis sets and was combined with extrapolation to the one-particle basis-set limit, core-correlation effects, scalar relativistic effects, and higher-order correlation effects to yield accurate QFFs. A “best-guess” reference geometry was determined at the CCSD(T)/5Z level of theory. Analytical transformation removes nonzero gradients to facilitate a second-order perturbation theory spectroscopic analysis. The QFF is transformed into Morse/cosine coordinates in order to perform exact vibrational configuration interaction computations. Equilibrium structures, vibrational frequencies, rotational constants, and selected spectroscopic constants are reported in comparison with experim...