Search

Dual Coordinate Descent (DCD) is one of the most popular optimization methods. The parallelization of DCD is difficult, as DCD is sequential in nature. As such, simultaneously running multiple DCD threads on batches of data elements causes result inaccuracy and slow convergence, due to the concurrent updates of multiple coordinates. Some parallelization methods adopt separable approximate functions that depend on the degree of parallelism. Such dependencies result in both poor scalability and slow convergence. To address these challenges, in this paper we present a new parallel primal-dual algorithm for DCD, called PPD. In PPD, the block data distribution is utilized to obtain a new approximate function that is independent from the parallelism. Subsequently we design a novel primal-dual acceleration scheme based on PPD to approach the optimal solution closely and quickly. We demonstrate the advantages of PPD in terms of scalability and efficiency through experiments.

A highly accurate CO

ConspectusWe review the Best Theory + Reliable High-Resolution Experiment (BTRHE) strategy for obtaining highly accurate molecular rovibrational line lists with InfraRed (IR) intensities. The need for highly accurate molecular rovibrational line lists is twofold: (a) assignment of the many rovibrational lines for common stable molecules especially those that exhibit a large amplitude motion, such as NH

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 semi-empirical molecular rovibrational IR line lists, such as ExoMol, TheoReTs, and Ames, combine the experimental accuracy and theoretical power to reach better than 0.1 cm−1 accuracy for line positions and better than 80–90% agreement for line intensities. The quality of these existing semi-empirical IR lists allows further improvements of intensity and line positions for those unobserved minor isotopologues. This paper presents our new BTRHE (Best Theory + Reliable High-resolution Experiment) strategy implementation. For line intensity, the isotopologue consistency and the patterns of mass dependence in the Ames-296 K SO2 and CO2 IR lists are quantitatively presented along the mass-inverse coordinates. The consistency and patterns are better than those in existing experimental data. The methodology proposed here can be used to identify inconsistencies, outliers, and mistakes in intensities, and help improve Effective Dipole Model (EDM) and molecular IR databases. We call for an experimental study on the 50006 and 60007 bands of CO2 628. For line position predictions, a simple approach combining the variational IR line lists with Effective Hamiltonian (EH) model may refine the effective rotational constants A0/B0/C0 and quartic centrifugal distortion constants of minor isotopologues. The prediction accuracy may be improved by two orders of magnitude, i.e. reaching 0–5 MHz prediction accuracy in the range of J

The quality of Ames-296K SO2 Infrared (IR) line list intensities is first validated by quantitative exploration of several dipole moment surfaces (DMSs) and partition sum convergence. The DMSs are computed with several of Dunning's correlation-consistent basis sets and their vibrational dependence are compared to the empirical model derived from Stark effect experiments reported by D. Patel, D. Margolese, and T.R. Dyke [J. Chem. Phys. 70, 2740 (1979)]. The effective dipole deviations from the DMS adopted in the Ames IR lists is 0.2–0.4% for vibrational states up to 3ν3. The vibrational dependence of the dipole moment is also in good agreement, except for nν1. Partition sum convergence at 296 K is confirmed by new calculations with rotational quantum number J up to 150 and upper state E’ up to 8000 cm−1. The isotopologue consistency of the Ames IR line lists is superior relative to the regular Effective Hamiltonian (EH) models and Effective Dipole Moment (EDM) models. The ν1 + ν2 and ν2 + ν3 intensity consistency check reveals the recently reported experimental intensities need significant improvement or re-analysis. After the accuracy, convergence, and isotopologue consistency have been confirmed, the theoretical Ames-296K intensities are combined with the experimental line positions or EH models that experimental spectroscopists published after 2009. Three high-resolution IR line sets are reported for the 32/33/34S16O2, 32S18O2 and 16O32S18O isotopologues: (1) the “New Lines Sets” include experimentally measured line positions; (2) the “Expanded Line Sets” include possible transitions among new rovibrational levels assigned in experiments and ground state (GS) levels predicted by reliable EH models; (3) the “Ames + MARVEL Sets” include possible transitions among all those levels reported in a recent MARVEL analysis. [Tobias et al, JQSRT 208, 152 (2018)]. Compared to the limited data in High-resolution TRANsmission molecular absorption database (HITRAN), these line sets have significantly improved the data coverage up to 4000 cm−1. Some missing bands can be traced to the unpublished experimental data. The isotopologue consistency of these line sets will help identify the uncertainties and defects in the experimental EH and EDM models. These line sets are good candidates for the next HITRAN update, if line shape parameters are available. The line sets can be downloaded from supplementary files or from the Ames Molecular Spectroscopic Database at http://huang.seti.org .

Total internal partition sums (TIPS) are reported for 166 isotopologues of 51 molecules important in planetary atmospheres. Molecules 1 to 50 are taken from the HITRAN2016 list, and, in some cases, additional isotopologues are considered for some of the molecules. Molecules 51–53 are C3H4, CH3, and CS2, respectively. TIPS are not reported for the O atom and CF4; thus, while there are 53 species in the list, data are reported for 51 molecules. The TIPS are determined by various methods from 1 K to a Tmax that ensures the TIPS reported have converged. These data are provided with HITRAN2016 and a new version of the TIPS code is available in both FORTRAN and python languages.

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.

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.

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.

Reliable infrared (IR) line lists are reported for the 13 isotopologues of carbon dioxide in HITRAN notation: 626, 636, 628, 627, 828, 727, 827, 638, 637, 737, 838, 738, and 646. Three IR lists are available for each istotopologue: a complete list at 296 K, a reduced-size list at 296 K, plus a reduced-size list at 1000 K. They are denoted Ames-296K, Ames-296K.reduced and Ames-1000K.reduced. With J up to 150, and energy up to 18,000 cm−1 above the zero point energy, these lists are expected to cover the temperature range up to 1500 K. Line shape parameters including temperature dependence are calculated and reported for four temperature ranges: Mars, Earth, Venus, and Hotter (700–2000 K). Comparisons are made against the available transition data in the HITRAN2012 models. Line position accuracy for most transitions up to 10,000–13,000 cm−1 is better than 0.03–0.05 cm−1. Computed transition intensities agree well with most HITRAN data but there exist suspicious exceptions for isotopologues. These line lists will expedite CO2 IR experimental data analysis and provide the scientific community with trustworthy alternatives for unknown IR bands. These line lists may be combined with existing experimental databases to facilitate the analysis of future laboratory experiments or astronomical observations.

The simple aromatic hydrocarbon, cyclopropenylidene (c-C

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.

Quartic force fields (QFFs) have been shown to be one of the most effective ways to efficiently compute vibrational frequencies for small molecules. In this letter we 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 describe the vibrational (or rovibrational) energy levels of an arbitrary molecular system to 5 cm −1 or better compared to experiment. 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.

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.

Nearly 4800 features of ammonia between 6300 and 7000 cm−1 with intensities ≥4×10−24 cm−1/(molecule·cm−2) at 296 K were measured using 16 pure NH3 spectra recorded at various temperatures (296–185 K) with the McMath–Pierce Fourier Transform Spectrometer at Kitt Peak National Observatory, AZ. The line positions and intensities were retrieved by fitting individual spectra based on a Voigt line shape profile and then averaging the values to form the experimental linelist. The integrated intensity of the region was 4.68×10−19 cm−1/(molecule·cm−2) at 296 K. Empirical lower state energies were also estimated for 3567 absorption line features using line intensities retrieved from 10 spectra recorded at gas temperature between 185 and 233 K. Finally, using Ground State Combination Differences (GSCDs) and the empirical lower state energy estimates, the quantum assignments were determined for 1096 transitions in the room temperature linelist, along with empirical upper state energies for 434 levels. The assignments correspond to seven vibrational states, as confirmed from recent ab initio calculations. The resulting composite database of 14NH3 line parameters will provide experimental constraints to ab initio calculations and support remote sensing of gaseous bodies including the atmospheres of Earth, (exo)planets, brown dwarfs, and other astrophysical environments.

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.

Highly accurate IR line lists are ready for 13 CO2 isotopologues and 7 SO2 isotopologues. Ames-296K IR lists carry 0.01 - 0.03 cm−1 accuracy up to 13,000 cm−1 for CO2 and 5500 cm−1 for SO2, and 90–95% intensity agreement for most observable bands. Good for atmospheric modeling on Venus and Exoplanets.

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

We report variational calculations of rovibrational energies of CH 4 using the code MULTIMODE and an ab initio force field of Schwenke and Partridge. The systematic convergence of the energies with respect to the level of mode coupling is presented. Converged vibrational energies calculated using the five-mode representation of the potential for zero total angular momentum are compared with previous, benchmark calculations based on Radau coordinates using this force field for zero total angular momentum and for J = 1. Very good agreement with the previous benchmark calculations is found.

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

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

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

We report full dimensional calculations of vibrational energies of H3O+ and D3O+ using two implementations of the code MULTIMODE. In one implementation, the reference geometry is the minimum of the potential (the standard choice for MULTIMODE). This implementation is not able to readily describe splittings in the vibrational energies due to motion through the inversion barrier. A second implementation, in which the reference geometry is the inversion saddle point, is able to describe the splittings. These full dimensional calculations are done using the realistic, though not spectroscopically accurate, potential of Ojamae, Singer and Shavitt, and the results are compared with experiment.

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.