Search

A new 626-based Ames-2 PES refinement and Ames-2016 line lists for 13 CO2 isotopologues are reported. A consistent σRMS = ±0.02 cm−1 is established for hundreds of isotopologue band origins using the Ames-2 PES. Ames-2016 line lists are computed at 296 K, 1000 K and 4000 K using the Ames-2 PES and the same DMS-N2 dipole surface used previously, with J up to 150, E′ up to 24,000 cm−1 or 18,000 cm−1 and appropriate intensity cutoffs. The lists are compared to the CDSD-296, CDSD-4000 databases, UCL line lists, and a few recent highly accurate CO2 intensity measurements. Both agreements and discrepancies are discussed. Compared to the old Ames CO2 lists, the Ames-2016 line lists have line position deviations reduced by 50% or more, which consequently leads to more reliable intensities. The line shape parameters in the Ames-2016 line lists are predicted using the newly assigned conventional vibrational polyad quantum numbers for rovibrational levels below 12,000 cm−1 so the quality of the line shape parameters is similar to that of CDSD or HITRAN. This study further proves that a semi-empirically refined PES (Ames-1 and Ames-2) coupled with a high quality ab initio DMS (DMS-N2 and UCL) may generate IR predictions with consistent accuracy and is thus helpful in the analysis of laboratory spectra and simulations of various isotopologues. The Ames-2016 lists based on DMS-N2 have reached the ∼1% intensity prediction accuracy level for the recent 626 30013-00001 and 20013-00001 bands, but further quantification and improvements require sub-percent or sub-half-percent accurate experimental intensities. The inter-isotopologue consistency of the intensity prediction accuracies should have reached better than 1–3% for regular bands not affected by resonances. Since the Effective Dipole Models (EDM) in CDSD and HITRAN have 1–20% or even larger uncertainties, we show that the Ames lists can provide better alternative IR data for many hard-to-determine isotopologue bands. Comparison at 4000 K suggests that the Ames-4000 K 12C16O2 line list is reliable and consistent within the current cutoffs of J ≤ 150 and E′ ≤ 24,000 cm−1, but intensity contributions involving higher energy levels should not be omitted and future computations need to be converged up to at least 32,000 cm−1 or higher. The remaining issues are discussed regarding the source of energy level discrepancies, intensity underestimations by ∼50% for some weak bands, etc. and also future work.

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

New semi-empirical Infrared (IR) line lists for 12C16O2, Ames-296 K and Ames-1000 K, have been computed using a newly updated ab initio CCSD(T)/aug-cc-pVQZ dipole moment surface (denoted DMS-N2) and an empirically refined potential energy surface (Ames-1). J=0–150 rovibrational levels are computed up to 30,000 cm−1, and related transitions are cut off at 1E−42 cm molecule−1 (296 K) and 1E−36 cm molecule−1 (1000 K). These are the first line lists available to cover reliably the energy region as high as ~20,000 cm−1. Recent experimental data at 1.1 μm has confirmed the predicted intensities for the 50013-00001 and 50014-00001 band transitions have better than 90% agreement. Comparisons are made against the Wattson 750 K line list and HITEMP/HITRAN at 300 K, 500 K, 725 K, 1000 K, 1500 K, 2000 K and 3000 K. The temperature dependence and accuracy of the new Ames-296 K/1000 K line lists are investigated and we claim both line lists are capable of providing reliable opacities up to 18,000–23,000 cm−1, while the highest applicable wavenumber range drops as T rises. We suggest caution is used for T>1000 K simulations. Comparison to recent experiments at 1000 K, 1550 K and 1773 K shows that the Ames-1000 K line list and HITEMP perform similarly in the 3200–3800 cm−1 and 4600–5200 cm−1 ranges. In the 2000–2100 cm−1 range, Ames-1000 K yields better agreement relative to experiment. Existing problems and possible future solutions in the new Ames-296 K/1000 K line lists for line positions and intensities are also discussed.

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.

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.

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

Model infrared spectra for non-rotating H 2 O are calculated at 0 K, based on exact quantum, standard classical and semiclassical calculations. An accurate potential energy surface is used along with a realistic dipole function. An analysis of the classical and quantum spectrum in the harmonic approximation is presented at 0 K. This clearly reveals that the magnitude of the classical intensities is essentially arbitrary, depending on the total energy. Thus, the intensity of classical harmonic spectrum disagrees with the corresponding quantum one. A very simple correction to the classical spectrum is suggested that largely restores agreement with the harmonic quantum spectrum. A second, more general classical correction is also suggested, which, however, requires knowledge of the normal modes.

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

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.

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.

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.

The anion chemistry of the interstellar medium (ISM) has almost exclusively been limited to linear hydrocarbons and cyanocarbons. Of the hydrocarbons, only the even n C n H– chains have been detected in the ISM, and lines hypothesized to originate with b-C3H– have been conclusively linked to the corresponding cation, as originally claimed. However, no reason has yet been provided as to why other anions cannot form, and the cyclic form of C3H– is actually the lowest-energy isomer on the anion's potential energy surface. As such, this work provides the necessary rovibrational reference data for the potential detection of this anion in the ISM or related laboratory experiments. Improvements over previously calculated rovibrational spectroscopic constants are contained herein along with graphical depictions of the pure rotational spectra at 100 K, 40 K, 20 K, and 2.7 K.

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.

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.

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 global potential energy surface (PES) that includes short and long range terms has been determined for the NH(3) molecule. The singles and doubles coupled-cluster method that includes a perturbational estimate of connected triple excitations and the internally contracted averaged coupled-pair functional electronic structure methods have been used in conjunction with very large correlation-consistent basis sets, including diffuse functions. Extrapolation to the one-particle basis set limit was performed and core correlation and scalar relativistic contributions were included directly, while the diagonal Born-Oppenheimer correction was added. Our best purely ab initio PES, denoted "mixed," is constructed from two PESs which differ in whether the ic-ACPF higher-order correlation correction was added or not. Rovibrational transition energies computed from the mixed PES agree well with experiment and the best previous theoretical studies, but most importantly the quality does not deteriorate even up to 10 300 cm(-1) above the zero-point energy (ZPE). The mixed PES was improved further by empirical refinement using the most reliable J=0-2 rovibrational transitions in the HITRAN 2004 database. Agreement between high-resolution experiment and rovibrational transition energies computed from our refined PES for J=0-6 is excellent. Indeed, the root mean square (rms) error for 13 HITRAN 2004 bands for J=0-2 is 0.023 cm(-1) and that for each band is always