Second-order Vibrational Perturbation Theory Research Articles

Infrared spectroscopy of the syn-methyl-substituted Criegee intermediate: A combined experimental and theoretical study.

An IR-vacuum ultraviolet (VUV) ion-dip spectroscopy method is utilized to examine the IR spectrum of acetaldehyde oxide (CH3CHOO) in the overtone CH stretch (2νCH) spectral region. IR activation creates a depletion of the ground state population that reduces the VUV photoionization signal on the parent mass channel. IR activation of the more stable and populated syn-CH3CHOO conformer results in rapid unimolecular decay to OH + vinoxy products and makes the most significant contribution to the observed spectrum. The resultant IR-VUV ion-dip spectrum of CH3CHOO is similar to that obtained previously for syn-CH3CHOO using IR action spectroscopy with UV laser-induced fluorescence detection of OH products. The prominent IR features at 5984 and 6081cm-1 are also observed using UV + VUV photoionization of OH products. Complementary theoretical calculations utilizing a general implementation of second-order vibrational perturbation theory provide new insights on the vibrational transitions that give rise to the experimental spectrum in the overtone CH stretch region. The introduction of physically motivated small shifts of the harmonic frequencies yields remarkably improved agreement between experiment and theory in the overtone CH stretch region. The prominent features are assigned as highly mixed states with contributions from two quanta of CH stretch and/or a combination of CH stretch with an overtone in mode 4. The generality of this approach is demonstrated by applying it to three different levels of electronic structure theory/basis sets, all of which provide spectra that are virtually indistinguishable despite showing large deviations prior to introducing the shifts to the harmonic frequencies.

Probing Ion-Receptor Interactions in Halide Complexes of Octamethyl Calix[4]Pyrrole.

Ion receptors are molecular hosts that bind ionic guests, often with great selectivity. The interplay of solvation and ion binding in anion host-guest complexes in solution governs the binding efficiency and selectivity of such ion receptors. To gain molecular-level insight into the intrinsic binding properties of octamethyl calix[4]pyrrole (omC4P) host molecules with halide guest ions, we performed cryogenic ion vibrational spectroscopy (CIVS) of omC4P in complexes with fluoride, chloride, and bromide ions. We interpret the spectra using density functional theory, describing the infrared spectra of these complexes with both harmonic and anharmonic second-order vibrational perturbation theory (VPT2) calculations. The NH stretching modes of the pyrrole moieties serve as sensitive probes of the ion binding properties, as their frequencies encode the ion-receptor interactions. While scaled harmonic spectra reproduce the experimental NH stretching modes of the chloride and bromide complexes in broad strokes, the high proton affinity of fluoride introduces strong anharmonic effects. As a result, the spectrum of F-·omC4P is not even qualitatively captured by harmonic calculations, but it is recovered very well by VPT2 calculations. In addition, the VPT2 calculations recover the intricate coupling of the NH stretching modes with overtones and combination bands of CH stretching and NH bending modes and with low-frequency vibrations of the omC4P macrocycle, which are apparent for all of the halide ion complexes investigated here. A comparison of the CIVS spectra with infrared spectra of solutions of the same ion-receptor complexes in d3-acetonitrile and d6-acetone shows how ion solvation changes the ion-receptor interactions for the different halide ions.

Performance of EOM-CCSD(T)(a)*-Based Quartic Force Fields in Computing Fundamental, Anharmonic Vibrational Frequencies of Molecular Electronically Excited States with Application to the Ã1A″ State of :CCH2 (Vinylidene).

Highly accurate anharmonic vibrational frequencies of electronically excited states are not as easily computed as their ground electronic state counterparts, but recently developed approximate triple excited state methods may be changing that. One emerging excited state method is equation of motion coupled cluster theory at the singles and doubles level with perturbative triples computed via the (a)* formalism, EOMEE-CCSD(T)(a)*. One of the most employed means for the ready computation of vibrational anharmonic frequencies for ground electronic states is second-order vibrational perturbation theory (VPT2), a theory based on quartic force fields (QFFs),fourth-order Taylor series expansions of the potential portion of the internuclear Watson Hamiltonian. The combination of these two is herein benchmarked for its performance for use as a means of computing rovibrational spectra of electronically excited states. Specifically, the EOMEE-CCSD(T)(a)* approach employing a complete basis set extrapolation along with core electron inclusion and relativity (the so-called "CcCR" approach) defining the QFF produces anharmonic fundamental vibrational frequencies within 2.83%, on the average, of reported gas-phase experimentally assigned values for the test set including the states of HCF, HCCl, HSiF, HNO, and HPO. However, some states have exceptional accuracy in the fundamentals, most notably for ν2 of HCCl in which the CcCR QFF value is within 1.8 cm-1 at 927.9 cm-1 (or 0.2%) of the experiment. Additionally, this approach produces rotational constants to, on the absolute average, within 0.41% of available experimental data, showcasing notable accuracy in the computation of rovibronic spectral data. Furthermore, utilizing a hybrid approach composed of harmonic CcCR force constants along with a set of simple EOMEE-CCSD(T)(a)*/aug-cc-pVQZ QFF cubic and quartic force constants is faster than using pure CcCR and better represents those modes that suffer from numerical instability in the anharmonic portion of the QFF, implying that this so-called "CcCR + QZ" QFF approach may be the best for future applications. Finally, complete, rovibrational spectral data are provided for :CCH2, a molecule of potential astrochemical interest, in order to aid in its potential future experimental rovibronic characterization.

Development of a universal method for vibrational analysis of the terminal alkyne C≡C stretch.

The terminal alkyne C≡C stretch has a large Raman scattering cross section in the "silent" region for biomolecules. This has led to many Raman tag and probe studies using this moiety to study biomolecular systems. A computational investigation of these systems is vital to aid in the interpretation of these results. In this work, we develop a method for computing terminal alkyne vibrational frequencies and isotropic transition polarizabilities that can easily and accurately be applied to any terminal alkyne molecule. We apply the discrete variable representation method to a localized version of the C≡C stretch normal mode. The errors of (1) vibrational localization to the terminal alkyne moiety, (2) anharmonic normal mode isolation, and (3) discretization of the Born-Oppenheimer potential energy surface are quantified and found to be generally small and cancel each other. This results in a method with low error compared to other anharmonic vibrational methods like second-order vibrational perturbation theory and to experiments. Several density functionals are tested using the method, and TPSS-D3, an inexpensive nonempirical density functional with dispersion corrections, is found to perform surprisingly well. Diffuse basis functions are found to be important for the accuracy of computed frequencies. Finally, the computation of vibrational properties like isotropic transition polarizabilities and the universality of the localized normal mode for terminal alkynes are demonstrated.

Exploring Alternate Methods for the Calculation of High-Level Vibrational Corrections of NMR Spin-Spin Coupling Constants.

Traditional nuclear magnetic resonance (NMR) calculations typically treat systems with a Born-Oppenheimer-derived electronic wave function that is solved for a fixed nuclear geometry. One can numerically account for this neglected nuclear motion by averaging over property values for all nuclear geometries with a vibrational wave function and adding this expectation value as a correction to an equilibrium geometry property value. Presented are benchmark coupled-cluster singles and doubles (CCSD) vibrational corrections to spin-spin coupling constants (SSCCs) computed at the level of vibrational second-order perturbation theory (VPT2) using the vibrational averaging driver of the CFOUR program. As CCSD calculations of vibrational corrections are very costly, cheaper electronic structure methods are explored via a newly developed Python vibrational averaging program within the Dalton Project. Namely, results obtained with the second-order polarization propagator approximation (SOPPA) and density functional theory (DFT) with the B3LYP and PBE0 exchange-correlation functionals are compared to the benchmark CCSD//CCSD(T) and experimental values. CCSD//CCSD(T) corrections are also combined with literature CC3 equilibrium geometry values to form the highest-order vibrationally corrected values available (i.e., CC3//CCSD(T) + CCSD//CCSD(T)). CCSD//CCSD(T) statistics showed favorable statistics in comparison to experimental values, albeit at an unfavorably high computational cost. A cheaper CCSD//CCSD(T) + B3LYP method showed quite similar mean absolute deviation (MAD) values as CCSD//CCSD(T), concluding that CCSD//CCSD(T) + B3LYP is optimal in terms of cost and accuracy. With reference to experimental values, a vibrational correction was not worth the cost for all of the other methods tested. Finally, deviation statistics showed that CC3//CCSD(T) + CCSD//CCSD(T) vibrational-corrected equilibrium values deteriorated in comparison to CCSD//CCSD(T) attributed to the use of a smaller basis set or lack of solvation effects for the CC3 equilibrium calculations.

Scaling-up VPT2: A feasible route to include anharmonic correction on large molecules

Vibrational analysis plays a crucial role in the investigation of molecular systems. Though methodologies like second-order vibrational perturbation theory (VPT2) have paved the way to more accurate simulations, the computational cost remains a difficult barrier to overcome when the molecular size increases. Building upon recent advances in the identification of resonances, we propose an approach making anharmonic simulations possible for large-size systems, typically unreachable by standard means. This relies on the fact that, often, only portions of the whole spectra are of actual interest. Therefore, the anharmonic corrections can be included selectively on subsets of normal modes directly related to the regions of interest. Starting from the VPT2 equations, we evaluate rigorously and systematically the impact of the truncated anharmonic treatment onto simulations. The limit and feasibility of the reduced-dimensionality approach are detailed, starting on a smaller model system. The methodology is then challenged on the IR absorption and vibrational circular dichroism spectra of an organometallic complex in three different spectral ranges.

Influence of fourth-order vibrational corrections on semi-experimental (reSE) structures of linear molecules.

Semi-experimental structures (reSE) are derived from experimental ground state rotational constants combined with theoretical vibrational corrections. They permit a meaningful comparison with equilibrium structures based on high-level abinitio calculations. Typically, the vibrational corrections are evaluated with second-order vibrational perturbation theory (VPT2). The amount of error introduced by this approximation is generally thought to be small; however, it has not been thoroughly quantified. Herein, we assess the accuracy of theoretical vibrational corrections by extending the treatment to fourth order (VPT4) for a series of small linear molecules. Typical corrections to bond distances are on the order of 10-5 Å. Larger corrections, nearly 0.0002 Å, are obtained to the bond lengths of NCCN and CNCN. A borderline case is CCCO, which will likely require variational computations for a satisfactory answer. Treatment of vibrational effects beyond VPT2 will thus be important when one wishes to know bond distances confidently to four decimal places (10-4 Å). Certain molecules with shallow bending potentials, e.g., HOC+, are not amenable to a VPT2 description and are not improved by VPT4.

Fundamental vibrational frequencies of pnictogen (Pn) containing linear triatomic anions: OCPn - and SCPn - where Pn = N, P, As and Sb

The fundamental vibrational frequencies of the isolated OCAs − and SCAs − ions, as well as their heavier antimony analogues, are reported for the first time. A thorough analysis of basis set convergence of the bond lengths and harmonic vibrational frequencies, as well as the examination of anharmonic corrections, is conducted using the MP2 and CCSD(T) ab initio methods with robust basis sets for the OC P n − and SC P n − anions moving down the pnictogen series from Pn = N to Sb. Second-order vibrational perturbation theory (VPT2) and vibrational configuration interaction (VCI) theory are used to compute the fundamental vibrational frequencies of each triatomic anion approaching the MP2 and CCSD(T) complete basis set (CBS) limits. While fundamental vibrational frequencies have been reported for OCN − and SCN − in the gas phase, only experimental vibrational frequencies of the salts for the heavier pnictogen analogues (P and As) have been reported. Experimental Raman vibrational frequencies for OCAs − and SCAs − salts were found to be coincidentally in good agreement with the CCSD(T) VCI frequencies near the CBS limit, differing by at most 13 cm − 1 from the former and 9 cm − 1 from the latter. These modest differences are likely due to the presence of counter-ions and other environmental effects in the solid state, and the overall agreement between the salts and gas-phase frequencies is quite similar to that observed for the OCP − and SCP − ions.

Approximating large-basis coupled-cluster theory vibrational frequencies using focal-point approximations.

The focal-point approximation can be used to estimate a high-accuracy, slow quantum chemistry computation by combining several lower-accuracy, faster computations. We examine the performance of focal-point methods by combining second-order Møller-Plesset perturbation theory (MP2) with coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] for the calculation of harmonic frequencies and that of fundamental frequencies using second-order vibrational perturbation theory (VPT2). In contrast to standard CCSD(T), the focal-point CCSD(T) method approaches the complete basis set (CBS) limit with only triple-ζ basis sets for the coupled-cluster portion of the computation. The predicted harmonic and fundamental frequencies were compared with the experimental values for a set of 20 molecules containing up to six atoms. The focal-point method combining CCSD(T)/aug-cc-pV(T + d)Z with CBS-extrapolated MP2 has mean absolute errors vs experiment of only 7.3 cm-1 for the fundamental frequencies, which are essentially the same as the mean absolute error for CCSD(T) extrapolated to the CBS limit using the aug-cc-pV(Q + d)Z and aug-cc-pV(5 + d)Z basis sets. However, for H2O, the focal-point procedure requires only 3% of the computation time as the extrapolated CCSD(T) result, and the cost savings will grow for larger molecules.

Accuracy Meets Feasibility for the Structures and Rotational Constants of the Molecular Bricks of Life: A Joint Venture of DFT and Wave-Function Methods.

A fully unsupervised computational protocol is proposed with the aim of obtaining reliable structural properties for molecular bricks of life in the gas phase. The results of the new composite scheme approach spectroscopic accuracy at a moderate cost without any empirical parameter in addition to those of the underlying electronic structure method. The whole workflow is fully automated and provides optimized geometries and equilibrium rotational constants. Direct comparison with experimental ground state rotational constants can be performed thanks to the effective computation of vibrational corrections in the framework of second-order vibrational perturbation theory. The results for all the nucleic acid bases and several flexible molecules of biological or medicinal interest show that the accuracy of the new tool is close to that delivered by state-of-the-art composite wave function methods for small semirigid molecules.

Pbqff: Push-Button Quartic Force Fields.

pbqff is an open-source program for fully automating the production of quartic force fields (QFFs) and their corresponding anharmonic spectroscopic data. Rather than being a monolithic piece of code, it consists of several key modules including a generic interface to quantum chemistry codes and, notably, queuing systems; a molecular point group symmetry library; an internal-to-Cartesian coordinate conversion module; a module for the ordinary least-squares fitting of potential energy surfaces; and an improved second-order rotational and vibrational perturbation theory package for asymmetric and symmetric tops that handles type-1 and -2 Fermi resonances, Fermi resonance polyads, and Coriolis resonances. All of these pieces are written in Rust, a modern, safe, and performant programming language that has much to offer for scientific programming. This work introduces pbqff and its surrounding ecosystem, in addition to reporting new anharmonic vibrational data for c-(C)C3H2 and describing how the components of pbqff can be leveraged in other projects.

Gas-phase formation and spectroscopic characterization of the disubstituted cyclopropenylidenes c-C3(C2H)2, c-C3(CN)2, and c-C3(C2H)(CN)

Aims. The detection of c-C3HC2H and possible future detection of c-C3HCN provide new molecules for reaction chemistry in the dense interstellar medium (ISM) where R-C2H and R-CN species are prevalent. Determination of chemically viable c-C3HC2H and c-C3HCN derivatives and their prominent spectral features can accelerate potential astrophysical detection of this chemical family. This work characterizes three such derivatives: c-C3(C2H)2, c-C3(CN)2, and c-C3(C2H)(CN). Methods. Interstellar reaction pathways of small carbonaceous species are well replicated through quantum chemical means. Highly accurate cc-pVXZ-F12/CCSD(T)-F12 (X = D,T) calculations generate the energetics of chemical formation pathways as well as the basis for quartic force field and second-order vibrational perturbation theory rovibrational analysis of the vibrational frequencies and rotational constants of the molecules under study. Results. The formation of c-C3(C2H)2 is as thermodynamically and, likely, as stepwise favorable as the formation of c-C3HC2H, rendering its detectability to be mostly dependent on the concentrations of the reactants. Both c-C3(C2H)2 and c-C3(C2H)(CN) will be detectable through radioastronomical observation with large dipole moments of 2.84 D and 4.26 D, respectively, while c-C3(CN)2 has an exceedingly small and likely unobservable dipole moment of 0.08 D. The most intense frequency for c-C3(C2H)2 is v2 at 3316.9 cm–1 (3.01 μm), with an intensity of 140 km mol–1. The mixed-substituent molecule c-C3(C2H)(CN) has one frequency with a large intensity, v1, at 3321.0 cm–1 (3.01 μm), with an intensity of 82 km mol–1. The molecule c-C3(CN)2 lacks intense vibrational frequencies within the range that current instrumentation can readily observe. Conclusions. Both c-C3(C2H)2 and c-C3(C2H)(CN) are viable candidates for astrophysical observation, with favorable reaction profiles and spectral data produced herein, but c-C3(CN)2 will not be directly observable through any currently available remote sensing means, even if it forms in large abundances.

Изотопные эффекты в спектрах комплексов с водородными связями. Расчет колебательных спектров поглощения димеров (D-=SUB=-2-=/SUB=-CO)-=SUB=-2-=/SUB=- и D-=SUB=-2-=/SUB=-CO·sDF и тримеров D-=SUB=-2-=/SUB=-CO·s(DF)-=SUB=-2-=/SUB=- и (D-=SUB=-2-=/SUB=-CO)_2·sDF

The frequencies and intensities for the vibrational bands of absorption spectra of hydrogen-bonded (D2CO)2 and D2CO∙∙∙DF dimers, two D2CO∙∙∙(DF)2 trimers and four (D2CO)2∙∙∙DF trimers are calculated in the MP2/aug-cc-pVTZ approximation with the basis set superposition error taken into account. Anharmonic values of spectral parameters were obtained using the vibrational second-order perturbation theory. The influence of hydrogen bonds on the spectral parameters was determined from comparison of the values calculated for monomers, dimers, and trimers in the same approximation. The data obtained were compared with the results of previous calculations of (H2CO)2 and H2CO∙∙∙HF dimers and H2CO∙∙∙(HF)2 and (H2CO)2∙∙∙HF trimers. It was shown that one D2CO∙∙∙(DF)2 trimer and two (D2CO)2∙∙∙DF trimers have significant binding energies and strong absorption bands, which makes them promising candidates for detection by spectroscopic methods.

Theoretical Spectroscopic Study of Two Ketones of Atmospheric Interest: Methyl Glyoxal (CH3COCHO) and Methyl Vinyl Ketone (CH3COCH═CH2).

Two ketones of atmospheric interest, methyl glyoxal and methyl vinyl ketone, are studied using explicitly correlated coupled cluster theory and core–valence correlation-consistent basis sets. The work focuses on the far-infrared region. At the employed level of theory, the rotational constants can be determined to within a few megahertz of the experimental data. Both molecules present two conformers, trans/cis and antiperiplanar (Ap)/synperiplanar (Sp), respectively. trans-Methyl glyoxal and Ap-methyl vinyl ketone are the preferred structures. cis-Methyl glyoxal is a secondary minimum of very low stability, which justifies the unavailability of experimental data in this form. In methyl vinyl ketone, the two conformers are almost isoenergetic, but the interconversion implies a relatively high torsional barrier of 1798 cm–1. A very low methyl torsional barrier was estimated for trans-methyl glyoxal (V3 = 273.6 cm–1). Barriers of 429.6 and 380.7 cm–1 were computed for Ap- and Sp-methyl vinyl ketone. Vibrational second-order perturbation theory was applied to determine the rovibrational parameters. The far-infrared region was explored using a variational procedure of reduced dimensionality. For trans-methyl glyoxal, the ground vibrational state was estimated to split by 0.067 cm–1, and the two low excited energy levels (1 0) and (0 1) were found to lie at 89.588 cm–1/88.683 cm–1 (A2/E) and 124.636 cm–1/123.785 cm–1 (A2/E). For Ap- and Sp-methyl vinyl ketone, the ground vibrational state splittings were estimated to be 0.008 and 0.017 cm–1, respectively.

Quantum chemical rovibrational analysis of aminoborane and its isotopologues.

Aminoborane, H2 NBH2 and its isotopologues, H2 N10 BH2 , D2 NBD2 , and D2 N10 BD2 , have been studied by high-level ab initio methods. All calculations rely on multidimensional potential energy surfaces and dipole moment surfaces including high-order mode coupling terms, which have been obtained from electronic structure calculations at the level of explicitly correlated coupled-cluster theory, CCSD(T)-F12, or the distinguishable cluster approximation, DCSD. Subsequent vibrational structure calculations based on second-order vibrational perturbation theory, VPT2, and vibrational configuration interaction theory, VCI, were used to determine rotational constants, centrifugal distortion constants, vibrationally averaged geometrical parameters and (an)harmonic vibrational frequencies. The impact of core-correlation effects is discussed in detail. Rovibrational VCI calculations were used to simulate the gas phase spectra of these species and an in-depth analysis of the ν7 band of aminoborane is provided. Color-coding is used to reveal the identity of the individual progressions of the rovibrational transitions for this particular mode.

Rotational and Vibrational Spectra of the Pyridyl Radicals: A Coupled-Cluster Study.

Pyridyl is a prototypical nitrogen-containing aromatic radical that may be a key intermediate in the formation of nitrogen-containing aromatic molecules under astrophysical conditions. On meteorites, a variety of complex molecules with nitrogen-containing rings have been detected with nonterrestrial isotopic abundances, and larger nitrogen-containing polycyclic aromatic hydrocarbons (PANHs) have been proposed to be responsible for certain unidentified infrared emission bands in the interstellar medium. In this work, the three isomers of pyridyl (2-, 3-, and 4-pyridyl) have been investigated with coupled cluster methods. For each species, structures were optimized at the CCSD(T)/cc-pwCVTZ level of theory and force fields were calculated at the CCSD(T)/ANO0 level of theory. Second-order vibrational perturbation theory (VPT2) was used to derive anharmonic vibrational frequencies and vibrationally corrected rotational constants, and resonances among vibrational states below 3500 cm-1 were treated variationally with the VPT2+K method. The results yield a complete set of spectroscopic parameters needed to simulate the pure rotational spectrum of each isomer, including electron-spin, spin-spin, and nuclear hyperfine interactions, and the calculated hyperfine parameters agree well with the limited available data from electron paramagnetic resonance spectroscopy. For the handful of experimentally measured vibrational frequencies determined from photoelectron spectroscopy and matrix isolation spectroscopy, the typical agreement is comparable to experimental uncertainty. The predicted parameters for rotational spectroscopy reported here can guide new experimental investigations into the yet-unobserved rotational spectra of these radicals.

Accurate Quantum Chemical Spectroscopic Characterization of Glycolic Acid: A Route Toward its Astrophysical Detection.

The first step to shed light on the abiotic synthesis of biochemical building blocks, and their further evolution toward biological systems, is the detection of the relevant species in astronomical environments, including earthlike planets. To this end, the species of interest need to be accurately characterized from structural, energetic, and spectroscopic viewpoints. This task is particularly challenging when dealing with flexible systems, whose spectroscopic signature is ruled by the interplay of small- and large-amplitude motions (SAMs and LAMs, respectively) and is further tuned by the conformational equilibrium. In such instances, quantum chemical (QC) calculations represent an invaluable tool for assisting the interpretation of laboratory measurements or even observations. In the present work, the role of QC results is illustrated with reference to glycolic acid (CH2OHCOOH), a molecule involved in photosynthesis and plant respiration and a precursor of oxalate in humans, which has been detected in the Murchison meteorite but not yet in the interstellar medium or in planetary atmospheres. In particular, the equilibrium structure of the lowest-energy conformer is derived by employing the so-called semiexperimental approach. Then, accurate yet cost-effective QC calculations relying on composite post-Hartree–Fock schemes and hybrid coupled-cluster/density functional theory approaches are used to predict the structural and ro-vibrational spectroscopic properties of the different conformers within the framework of the second-order vibrational perturbation theory. A purposely tailored discrete variable representation anharmonic approach is used to treat the LAMs related to internal rotations. The computed spectroscopic data, particularly those in the infrared region, complement the available experimental investigations, thus enhancing the possibility of an astronomical detection of this molecule.

Isotope Effects in the Spectra of Hydrogen-Bonded Complexes. Calculation of the Structure and Vibrational Absorption Spectra of the H-=SUB=-2-=/SUB=-O...HF, H-=SUB=-2-=/SUB=-O...DF, D-=SUB=-2-=/SUB=-O...HF, and D-=SUB=-2-=/SUB=-O...DF Complexes

The frequencies and intensities of IR absorption bands of the H2O...HF, H2O...DF, D2O...HF, and D2O...DF hydrogen-bonded complexes are calculated using the second-order vibrational perturbation theory. The MP2/aug-cc-pVTZ method with the basis set superposition error taken into account is used to calculate the electronic wave functions in determining the equilibrium configuration, the potential energy and dipole moment surfaces of these complexes, as well as in calculating spectral parameters. It is shown that upon complexation the frequencies and intensities of the stretching vibration of HF (DF) molecules and the intensities of stretching vibrations of H2O (D2O) molecules change most significantly. A variational calculation of librational motion of the water molecule in the two-well potential explained the reason for the higher value of the fundamental transition frequency of the ν1(H-F) mode as compared to the frequency of this mode upon hot transition from the first excited librational state. The dependence of the intermode anharmonic interaction on isotope substitution was analyzed. Keywords: hydrogen bond, calculations of spectra of molecular complexes, anharmonic interactions, isotope effects.

An integrated protocol to study hydrogen abstraction reactions by atomic hydrogen in flexible molecules: application to butanol isomers.

This work presents a protocol designed to study hydrogen abstraction reactions by atomic hydrogen in molecules with multiple conformations. The protocol starts with the search and location of the conformers of the equilibrium structures using the TorsiFlex program. By a simple modification of the starting geometry of reactants, a Python script generates the input for the hydrogen abstraction transition states. Initially, the search of the stationary points (reactants and transition states) is carried out at a low-level employing firstly a preconditioned search and secondly a random search. The low-level conformers were reoptimized using a higher level electronic structure method. This information allows the evaluation of the multistructural harmonic-oscillator partition functions, which are corrected for zero-point energy anharmonicity by the hybrid degeneracy-corrected second-order vibrational perturbation theory and for torsional anharmonicity by the multistructural torsional method, as implemented in the MsTor program. The structural information of the stationary points is used by Pilgrim to evaluate the multipath canonical variational transition state theory thermal rate constants with multidimensional small-curvature corrections for tunneling. Therefore, the thermal rate constants include variational (recrossing) and tunneling effects in addition to the effect of multiple conformations on the thermal rate constants. These features grant the applicability of the method to a wide range of temperatures. The method was applied to each of the hydrogen abstraction sites of the four isomers of butanol. The methodology employed allowed us to calculate the thermal rate constants in the temperature range of 250-2500 K and to accurately fit them to analytical expressions. The variety of abstraction sites shows that the protocol is robust and that it can be employed to study hydrogen abstraction reactions in molecules containing carbon and oxygen as heavy atoms.

Slow monomer vibrations in formic acid dimer: Stepping up the ladder with FTIR and Raman jet spectroscopy.

In an effort to extend the cold gas phase spectroscopic database of the cyclic formic acid dimer (FAD), we present and analyze the jet-cooled vibrational infrared and Raman spectrum of (HCOOH)2 in the monomer fingerprintregion between 600 and 1500 cm-1. The present study bridges the gap between the intermolecular dimerization-induced and the carbonyl stretching fundamentals that have already been reexamined using jet-cooled or high-resolution spectroscopy. This completes the characterization of the jet-cooled vibrational (HCOOH)2 spectrum below the complex OH (CH) stretching fundamentals, and we report resonance-induced FAD combination/overtone transitions that will serve as a valuable reference for a theoretical modeling of its vibrational dynamics. As a by-product, several new formic acid trimer fundamentals are identified in the jet spectra and assigned with the help of second-order vibrational perturbation theory (VPT2). The polar formic acid dimer still eludes detection in a supersonic jet, but we are able to estimate an experimental upper-bound of the polar dimer-to-trimer-to-cyclic dimer intensity ratio to about 1:10:100 under typical expansion conditions. Using VPT2 with resonance treatment (VPT2+K), we reinvestigate the notorious ν22 resonance triad. Generally, we find that VPT2, which is, of course, inadequate for modeling the resonance-rich OH stretching spectrum of FAD, is performing very satisfactorily in predicting fundamental and two-quantum state term values for the slower modes below 1500 cm-1. As these modes are the building blocks for the ultrafast energy dissipation in the OH stretching region, the present work opens the door for its quantitative understanding.