New analysis of the 2ν4, ν4+ν6, 2ν6, ν3+ν4, ν3+ν6, ν1, ν5, ν2+ν4, 2ν3, ν2+ν6 and ν2+ν3 bands of formaldehyde H212C16O: Line positions and intensities in the 3.5μm spectral region

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.

On the separability of large-amplitude motions in anharmonic frequency calculations

Nuclear vibrational theories based upon the Watson Hamiltonian are ubiquitous in quantum chemistry, but are generally unable to model systems in which the wavefunction can delocalise over multiple energy minima, i.e. molecules that have low-energy torsion and inversion barriers. In a 2019 Chemical Reviews article, Puzzarini et al. note that a common workaround is to simply decouple these problematic modes from all other vibrations in the system during anharmonic frequency calculations. They also point out that this approximation can be "ill-suited", but do not quantify the errors introduced. In this work, we present the first systematic investigation into how separating out or constraining torsion and inversion vibrations within potential energy surface (PES) expansions affects the accuracy of computed fundamental wavenumbers for the remaining vibrational modes, using a test set of 19 tetratomic molecules for which high quality analytic potential energy surfaces and fully-coupled anharmonic reference fundamental frequencies are available. We find that the most effective and efficient strategy is to remove the mode in question from the PES expansion entirely. This introduces errors of up to +10 cm-1 in stretching fundamentals that would otherwise couple to the dropped mode, and ±5 cm-1 in all other fundamentals. These errors are approximately commensurate with, but not necessarily additional to, errors due to the choice of electronic structure model used in constructing spectroscopically accurate PES.

Structures of the dehydrogenation products of methane activation by 5d transition metal cations revisited: Deuterium labeling and rotational contours

A previous infrared multiple photon dissociation (IRMPD) action spectroscopy and density functional theory (DFT) study explored the structures of the [M,C,2H]⁺ products formed by dehydrogenation of methane by four, gas-phase 5d transition metal cations (M⁺ = Ta⁺, W⁺, Ir⁺, and Pt⁺). Complicating the analysis of these spectra for Ir and Pt was observation of an extra band in both spectra, not readily identified as a fundamental vibration. In an attempt to validate the assignment of these additional peaks, the present work examines the gas phase [M,C,2D]⁺ products of the same four metal ions formed by reaction with perdeuterated methane (CD₄). As before, metal cations are formed in a laser ablation source and react with methane pulsed into a reaction channel downstream, and the resulting products are spectroscopically characterized through photofragmentation using the free-electron laser for intracavity experiments in the 350-1800 cm^-1 range. Photofragmentation was monitored by the loss of D for [Ta,C,2D]⁺ and [W,C,2D]⁺ and of D₂ in the case of [Pt,C,2D]⁺ and [Ir,C,2D]⁺. Comparison of the experimental spectra and DFT calculated spectra leads to structural assignments for all [M,C,2H/2D]⁺ systems that are consistent with previous identifications and allows a full description of the systematic spectroscopic shifts observed for deuterium labeling of these complexes, some of the smallest systems to be studied using IRMPD action spectroscopy. Further, full rotational contours are simulated for each vibrational band and explain several observations in the present spectra, such as doublet structures in several bands as well as the observed linewidths. The prominent extra bands in the [Pt,C,2D/2H]⁺ spectra appear to be most consistent with an overtone of the out-of-plane bending vibration of the metal carbene cation structure.

A tunable diode laser absorption spectrometer for formaldehyde atmospheric measurements validated by simulation chamber instrumentation

A tunable diode laser absorption spectrometer (TDLAS) for formaldehyde atmospheric measurements has been set up and validated through comparison experiments with a Fourier transform infrared spectrometer (FT-IR) in a simulation chamber. Formaldehyde was generated in situ in the chamber from reaction of ethene with ozone. Three HCHO ro-vibrational line intensities (at 2909.71, 2912.09 and 2914.46 cm(-1)) possibly used by TDLAS were calibrated by FT-IR spectra simultaneously recorded in the 1600-3200 cm(-1) domain during ethene ozonolysis, enabling the on-line deduction of the varying concentration for HCHO in formation. The experimental line intensities values inferred confirmed the calculated ones from the updated HITRAN database. In addition, the feasibility of stratospheric in situ HCHO measurements using the 2912.09 cm(-1) line was demonstrated. The TDLAS performances were also assessed, leading to a 2sigma detection limit of 88 ppt in volume mixing ratio with a response time of 60 sec at 30 Torr and 294 K for 112 m optical path. As part of this work, the room-temperature rate constant of this reaction and the HCHO formation yield were found to be in excellent agreement with the compiled literature data.

Calculations of rovibrational energies and dipole transition intensities for polyatomic molecules using MULTIMODE

We report rigorous calculations of rovibrational energies and dipole transition intensities for three molecules using a new version of the code MULTIMODE. The key features of this code which permit, for the first time, such calculations for moderately sized but otherwise general polyatomic molecules are briefly described. Calculations for the triatomic molecule BF(2) are done to validate the code. New calculations for H(2)CO and H(2)CS are reported; these make use of semiempirical potentials but ab initio dipole moment surfaces. The new dipole surface for H(2)CO is a full-dimensional fit to the dipole moment obtained with the coupled-cluster with single and double excitations and a perturbative treatment of triple excitations method with the augmented correlation consistent triple zeta basis set. Detailed comparisons are made with experimental results from a fit to relative data for H(2)CS and absolute intensities from the HITRAN database for H(2)CO.