Anharmonic Vibrational Analysis of the Infrared and Raman Gas-Phase Spectra of s-trans- and s-gauche-1,3-Butadiene

VIBFREQ1295: A New Database for Vibrational Frequency Calculations

High-throughput approaches for producing approximate vibrational spectral data for molecules of astrochemistry interest rely on harmonic frequency calculations using computational quantum chemistry. However, model chemistry recommendations (i.e., a level of theory and basis set pair) for these calculations are not yet available and, thus, thorough benchmarking against comprehensive benchmark databases is needed. Here, we present a new database for vibrational frequency calculations (VIBFREQ1295) storing 1295 experimental fundamental frequencies and CCSD(T)(F12*)/cc-pVDZ-F12 ab initio harmonic frequencies from 141 molecules. VIBFREQ1295's experimental data was complied through a comprehensive review of contemporary experimental data, while the ab initio data was computed here. The chemical space spanned by the molecules chosen is considered in-depth and is shown to have good representation of common organic functional groups and vibrational modes. Scaling factors are routinely used to approximate the effect of anharmonicity and convert computed harmonic frequencies to predicted fundamental frequencies. With our experimental and high-level ab initio data, we find that a single global uniform scaling factor of 0.9617(3) results in median differences of 15.9(5) cm^-1. A far superior performance with a median difference of 7.5(5) cm^-1 can be obtained, however, by using separate scaling factors (SFs) for three regions: frequencies less than 1000 cm^-1 (SF = 0.987(1)), between 1000 and 2000 cm^-1 (SF = 0.9727(6)), and above 2000 cm^-1 (SF = 0.9564(4)). This sets a lower bound for the performance that could be reliably obtained using scaling of harmonic frequency calculations to predict experimental fundamental frequencies. VIBFREQ1295's most important purpose is to provide a robust database for benchmarking the performance of any vibrational frequency calculations. VIBFREQ1295 data could also be used to train machine-learning models for the prediction of vibrational spectra and as a reference and data starting point for more detailed spectroscopic modeling of particular molecules. The database can be found as part of the Supporting Information for this paper or in the Harvard DataVerse at https://doi.org/10.7910/DVN/VLVNU7.

Calculation of absolute Raman scattering cross-sections using vibrational self-consistent field/vibrational configuration interaction wave functions

In the present study, the differential scattering cross-sections, depolarization ratios and Raman shifts of small molecular systems are obtained from configuration iteration wave functions of vibrational self-consistent field (VSCF) states. The transition polarizabilities were modeled using the Placzek approximation, neglecting those contributions not arising from the electric dipole mechanism. This theoretical approach is considered a good approximation for samples that absorb in the UV range if the excitation radiation falls in the visible region, as is the case of the molecules selected for the present study, namely: water, methane, and acetylene. Potential energy and electronic polarizability surfaces are calculated by the CCSD(T) and CC3 methods with aug-cc-p(C)V(T,Q,5)Z basis sets. The vibrational Hamiltonian includes the vibrational angular momentum contribution of the Watson kinetic energy operator. As expected, due to the variational nature of the VSCF and vibrational configuration interaction (VCI) methods, the Raman transition wavenumbers are substantially improved over the harmonic predictions. Surprisingly, the scattering cross-sections obtained using the harmonic approximation or the VSCF method better agrees with the experimental values than those cross-sections predicted using VCI wave functions. The more significant deviations of the VCI results from the experimental reference may be related to the significant uncertainties of the measured cross-sections. Still, it may also indicate that the VCI Raman transition moments may require a more accurate description of the electronic polarizability surface. Finally, the depolarization ratios calculated for H₂ O and C₂ D₂ using harmonic and VCI wave functions have similar accuracy, whereas, for C₂ H₂ and C₂ HD, the VCI results are more accurate.

Fundamental studies of vibrational resonance phenomena by multivalued resummation of the divergent Rayleigh-Schrödinger perturbation theory series: deciphering polyad structures of three H216O isotopologues

A fundamental quantitative study of the vibrational resonances of three H₂¹⁶O H,D-isotopologues with a quartic Watson Hamiltonian was carried out using the resummation of the high-order (∼λⁿ, n ≤ 203) divergent Rayleigh-Schrödinger perturbation theory (RSPT) series by quartic Padé-Hermite multivalued diagonal approximants PH^[m,m,m,m,m], m ≤ 40. The resonance condition between a pair of states is formulated as the existence of a common complex energy solution branch point inside the unit circle: |E(λ_j)|,  Re(λ_j)² + Im(λ_j)² ≤ 1. For the matrix formulation of the vibrational problem (VCI), the existence of common branch points is governed by the Katz theorem and they can be found as roots of discriminant polynomials. The main branches of the Padé-Hermite approximants typically reproduce VCI energies with high accuracy while alternative branches often fit nearby resonant states. The resummation of the RSPT series for H₂O and D₂O (up to the tenth polyad) revealed not only Fermi and Darling-Dennison resonances, but also unusual (0,2,-5) and (5,-2,0) resonance effects matching the (5,2,5) polyad structure, while the (3,2,1) structure was rigorously confirmed for HDO. It is demonstrated that the (5,2,5) polyad structure ensures good organization of high-energy resonating states and breaks down the classic (2,1,2) structure. The advocated methodology of quantitative description of resonance phenomena and revealing polyad structures is suitable for larger molecules and can be adapted to linear molecules and symmetric tops. Its application ensures rigorous classification of vibrational states and can be used in quantitative vibration-rotation spectroscopy.

On the impact of side methyl groups on the structure and vibrational properties of β-carotenoids. The case of butadiene and isoprene

Theoretical consideration about the impact of methyl groups on the structure and vibrational properties of β-carotenoids, using medium size molecules of trans-butadiene and trans-isoprene, are reported. Density functional theory (DFT) calculations with correlation-consistent and polarization-consistent basis sets were applied to trans-1,3-butadiene and trans-isoprene as the smallest building bricks of β-carotenoids. Their structure and harmonic vibrations were estimated in the complete basis set limit (CBS) using the non-linear least square fit. Optimized geometries and harmonic frequencies, obtained with B3LYP and BLYP density functionals and large basis sets, were favorably reproduced by a significantly faster approach, using a recently modified STO(1M)-3G Slater-type basis set. Selected density functionals with STO(1M)-3G and 6-311++G** basis sets were also successful in predicting β-carotene structures and harmonic vibrations. This work demonstrates the potential applicability of the proposed level of theory for larger molecules, including β-carotenoids, present in numerous natural food sources. The proposed scheme of molecular modeling, applied to biologically active compounds in food, could provide a deeper insight into their function in vivo, which is directly related to their structure and spectroscopic properties. It could also support the experimental qualitative analysis, based on peak assignment of β-carotenoids in various food sources.

Overtone spectroscopy of v(C=O) stretching vibration of hexafluoroacetone: Experimental and ab initio determination of peak positions, absolute intensities, and band shapes

Infrared overtone spectra of the ν(C = O) stretching vibration (ν₁) of gaseous hexafluoroacetone ((CF₃)₂C = O, HFA) were recorded in the spectral range of 7450-3300 cm^-1 with a resolution of 0.1 cm^-1. Experimental absolute IR intensities and vibrational band centers of the overtones 2ν₁, 3ν₁, 4ν₁ of HFA were measured and compared with their ab initio counterparts, calculated by the second-order canonical vibrational perturbation theory (CVPT2). A hybrid potential energy surface (PES) was evaluated using the quantum-mechanical models MP2/cc-pVQZ for harmonic and MP2/cc-pVTZ for anharmonic parts. Cubic surfaces of dipole moment components were determined using the MP2/cc-pVTZ model. The predicted IR intensities for the first and second overtones reproduced the experimental values with a discrepancy of 6% and 9%, respectively. A weak Fermi resonance between the first overtone 2ν₁ and combination tone ν₁ + ν₂ + ν₈ was predicted. The appropriate model was employed for simulating the bands studied using: (i) the asymmetric-top vibration-rotational structure of the ν₁ mode, (ii) the inhomogeneous band structure due to contributions of the vibrational states populated at the room temperature, and (iii) the homogeneous broadening of each vibration-rotation transition due to the intramolecular vibrational redistribution (IVR). The rotational and vibrational anharmonic constants were taken from the ab initio calculations, whereas the IVR data were obtained from the experimental data of Chekalin et al. (JPCA 2014, 118:955) with a time resolution of ≈100 fs. A good level of agreement of the predicted shapes of bands studied with experiment is achieved.

Computational studies, NMR, Raman and infrared spectral analysis of centrosymmetric (2Z,4Z)-Hexa-2,4-dienedinitrile

The Raman (50–3500[Formula: see text]cm[Formula: see text]) and infrared (200–3500[Formula: see text]cm[Formula: see text]) spectra of (2Z,4Z)-Hexa-2,4-dienedinitrile (C6H4N2; cis,cis-HDDN) have been recorded. Initially, three conformers were proposed based on following point groups: C2h, C2v and [Formula: see text] (gauche). For comparison purposes, M05-2X, M06-2X method were used in addition to B3LYP, MP2[Formula: see text]full and MP4[Formula: see text]full quantum mechanical calculations employing 6-31G(d) and 6-311+G(d,p) basis sets using Gaussian 09 program. Aided by potential energy surface scan, the C2h conformer is verified as a global minimum and the gauche ([Formula: see text]) being a local minimum with an energy barrier of 2.82[Formula: see text]kcal/mol. The energy difference ranged from 4.87 to 13.70[Formula: see text]kcal/mol favors the C2h conformer in good agreement with Raman and infrared spectral analysis. These results were also supported by the observed [Formula: see text] value (H9-H[Formula: see text]) at 10.8[Formula: see text]Hz which is consistent to 10.19[Formula: see text]Hz estimated for centrosymmetric conformer (C2h) rather than 4.81[Formula: see text]Hz predicted for gauche ([Formula: see text]). Complete vibrational assignments are proposed herein based on normal coordinate analysis (NCA) and potential energy distributions (PEDs) combined with theoretical vibrational frequencies and force constants in internal coordinates. It is crucial that NCA using VEDA 4 program with mixing shows large deviations from ours due to neglecting symmetry elements while mixing molecular vibrational motions. The results of these spectroscopic and theoretical studies are reported herein and compared with similar molecules whenever appropriate.

The Hunt for Elusive Molecules: Insights from Joint Theoretical and Experimental Investigations

Rotational spectroscopy is an invaluable tool to unambiguously determine the molecular structure of a species, and sometimes even to establish its very existence. This article illustrates how experimental and theoretical state-of-the-art tools can be used in tandem to investigate the rotational structure of molecules, with particular emphasis on those that have long remained elusive. The examples of three emblematic species-gauche-butadiene, disilicon carbide, and germanium dicarbide-highlight the close, mutually beneficial interaction between high-level theoretical calculations and sensitive microwave measurements. Prospects to detect other elusive molecules of chemical and astronomical interest are discussed.

Anharmonic vibrational analysis of s-trans and s-cis conformers of acryloyl fluoride using numerical-analytic Van Vleck operator perturbation theory

A new gas-phase infrared (IR) spectrum of acryloyl fluoride (ACRF, CH₂CHCFO) with a resolution of 0.1cm^-1 in the range 4000-450cm^-1 was measured. Theoretical ab initio molecular structures, full quartic potential energy surfaces (PES), and cubic surfaces of dipole moments and polarizability tensor components (electro-optical properties, EOP) of the s-trans and s-cis conformers of the ACRF were calculated by the second-order Møller-Plesset electronic perturbation theory with a correlation consistent Dunning triple-ζ basis set. The numerical-analytic implementation of the second-order operator canonical Van Vleck perturbation theory was employed for predicting anharmonic IR and Raman scattering (RS) spectra of ACRF. To improve the anharmonic predictions, harmonic frequencies were replaced by their counterparts evaluated with the higher-level CCSD(T)/cc-pVTZ model, to form a "hybrid" PES. The original operator representation of the Hamiltonian is analytically reduced to a quasi-diagonal form, integrated in the harmonic oscillator basis and diagonalized to account for strong resonance couplings. Double canonical transformations of EOP expansions enabled prediction of integral intensities of both fundamental and multi-quanta transitions in IR/RS spectra. Enhanced band shape analysis reinforced the assignments. A thorough interpretation of the new IR experimental spectra and existing matrix-isolation literature data for the mixture of two conformers of ACRF was accomplished, and a number of assignments clarified.