Theoretical structural and spectroscopic characterization of peroxyacetic acid (CH3-CO-OOH): study of the far infrared region.

Peroxyacetic acid, a non-rigid oxygenated organic molecule which acts in the atmosphere as a reservoir of HOX and ROX radicals, is studied using highly correlated ab initio methods with the aim of its spectroscopic characterization in the gas phase. The study focuses on the far infrared region providing reliable rovibrational parameters such as energy levels and splittings. The molecule presents three conformers that inter-convert by internal rotation, drawing a potential energy surface of 12 minima. One of them shows prominent stability due to the formation of one weak intramolecular bond between the hydrogen atom of the hydroperoxy group and the oxygen atom of the carbonyl group. For the three minimum energy structures, rotational constants and centrifugal distortion constants are provided. It may be expected that the most stable conformer is the only one contributing to the spectral features in further measurements at low temperature. In this structure, the methyl torsional barrier has been found to be very low, V3 = 88.6 cm-1 producing a splitting of 2.262 cm-1 for the ground vibrational state. The study confirms that the ν20 torsional mode interacts strongly with the other two torsional modes ν19 and ν21, but slightly with the remaining vibrations. Then, a variational procedure in three dimensions allows the exploration of the low-frequency modes. The methyl torsional fundamental ν21 was found to be 49.1 cm-1 (Ai) and 33.4 cm-1 (E). The fundamentals of ν20 (C-O bond torsion) and ν19 (OH torsion) have been computed to be 216.7 cm-1 (A2) and 218.5 cm-1 (E) and 393.6 cm-1 (A2) and 394.1 cm-1. Since non-rigidity can have effects on the reactivity due to the conformer interconversion, and transitions involving low-lying levels can be observed with many spectroscopic techniques, this work can help kinetic studies and assignments of further spectroscopic studies needed for the detection in the gas phase of trace molecules.

High resolution far-infrared synchrotron spectroscopy of 2-furfural conformers: Fundamental and hot bands.

In the continuity of a previous jet-cooled rovibrational study of trans and cis conformers of 2-furfural in the mid-infrared region (700-1750 cm-1) [Chawananon et al., Molecules 28 (10), 4165 (2023)], the present work investigates the far-infrared spectroscopy of 2-furfural using a long path absorption cell coupled to a high-resolution Fourier transform spectrometer and synchrotron radiation at the AILES beamline of the SOLEIL synchrotron. Guided by anharmonic calculations, vibrational energy levels and excited-state rotational constants are sufficiently predictive for a complete assignment of all fundamental and combination bands up to 700 cm-1, as well as the rovibrational analysis of 4 (1) low-frequency modes of trans-(cis-)2-furfural. A global rovibrational simulation, including far-infrared rovibrational lines and microwave and millimeter-wave rotational lines assigned in a previous study [Motiyenko et al., J. Mol. Spectrosc., 244, 9 (2007)] provides a reliable set of ground- and excited-state rotational parameters involving ring torsion, bending, and ring puckering modes of 2-furfural. In a second step, a rovibrational analysis of several hot band sequences, mainly involving the lowest frequency ring CHO torsion mode, is carried out. Reliable values of some anharmonic coefficients are obtained experimentally and could serve as a benchmark for validating advanced anharmonic calculations related to these large amplitude motions of flexible molecules.

The Role of Water in the Adsorption of Nitro-Organic Pollutants on Activated Carbon.

The density functional theory (DFT) is applied to theoretically study the capture and storage of three different nitro polycyclic aromatic hydrocarbons, 4-nitrophenol, 2-nitrophenol, and 9-nitroanthracene by activated carbon, with and without the presence of water. These species are pollutants derived from vehicle and industry emissions. The modeling of adsorption is carried out at the molecular level using a high-level density functional theory with the B3LYP-GD(BJ)/6-31+G(d,p) level of theory. The adsorption energies of polluting gases considered isolated and in a humid environment are compared to better understand the role of water. The calculations reveal different possible pathways involving the formation of chemical bonds between adsorbent and adsorbate on the formation of intermolecular van der Waals interactions. The negative adsorption energy on AC for the three species is obtained when they are treated individually and in mixture with H2O. The basis-set superposition error, estimated using the counterpoise correction, varies the adsorption energies by 2-13%. Dispersion effects were also taken into account. The adsorption energy ranges from -10 to -414 kJ/mol suggesting a diversity of pathways. The resulting analysis suggests three preferred pathways for capture. The main pathway is physical interaction due to π-π stacking. Other means are capture due to the formation of hydrogen bonds resulting from water adsorbed on the surface and the simultaneous adsorption of pollutant and water where water can act as a link that promotes adsorption. The thermodynamic properties give a clue to the most eco-friendly approaches for molecular adsorption.

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.

Collision excitation of c-C3H-(X1A1) by He.

Accurate modeling of anionic abundances in the interstellar and circumstellar media requires calculations of collisional data with the most abundant species that are usually He atoms and H2 molecules. In this paper, we focus on smaller cyclic molecular anion, c-C3H-, an astrophysical candidate, following the detection of larger CnH- carbon chains. From a new three-dimensional potential energy surface, the rotational (de-)excitation of the c-C3H-(X1A1) anion by collision with He is investigated. The surface is obtained in the supermolecular approach at the CCSD(T)-F12/aug-cc-pVTZ level of theory. Fully quantum close-coupling calculations of inelastic integral cross sections are performed on a grid of collisional energies large enough to ensure the convergence of the state-to-state rate coefficients for the 34 first rotational levels up to jKa,Kc = 77,0 of c-C3H- and temperatures ranging from 5 to 100 K. For this collisional system, rate coefficients exhibit a strong dominance in favor of 21,2 → l1,1 downward transition. This transition was previously used for the detection of the cyclic parent c-C3H. The c-C3H--He rate coefficients (∼10-11 cm3 s-1) are of the same order of magnitude as those of the detected anions CnH- (as C2H-, C4H-, and C6H-) in collision with He and one order of magnitude smaller than those with H2. The critical densities of H2 were also estimated, and a discussion on the validity of the local thermodynamic equilibrium conditions is carried out. This work represents the contribution to understanding and modeling abundances and chemistry of hydrocarbon radicals, CnH, in astrophysical media.

Large Amplitude Motions of Pyruvic Acid (CH3-CO-COOH).

Torsional and rotational spectroscopic properties of pyruvic acid are determined using highly correlated ab initio methods and combining two different theoretical approaches: Second order perturbation theory and a variational procedure in three-dimensions. Four equilibrium geometries of pyruvic acid, Tc, Tt, Ct, and CC, outcome from a search with CCSD(T)-F12. All of them can be classified in the Cs point group. The variational calculations are performed considering the three internal rotation modes responsible for the non-rigidity as independent coordinates. More than 50 torsional energy levels (including torsional subcomponents) are localized in the 406-986 cm-1 region and represent excitations of the ν24 (skeletal torsion) and the ν23 (methyl torsion) modes. The third independent variable, the OH torsion, interacts strongly with ν23. The A1/E splitting of the ground vibrational state has been evaluated to be 0.024 cm-1 as it was expected given the high of the methyl torsional barrier (338 cm-1). A very good agreement with respect to previous experimental data concerning fundamental frequencies (νCAL - νEXP ~ 1 cm-1), and rotational parameters (B0CAL - B0EXP < 5 MHz), is obtained.

Computational analysis of the far infrared spectral region of various deuterated varieties of ethylene glycol.

The far infrared spectra of three deuterated isotopologues of ethylene glycol, CH2OD-CH2OD, CH2OD-CH2OH and CH2OH-CH2OD, where the latter two species differ in their intramolecular hydrogen-bonding arrangement, are studied using highly correlated ab initio methods, vibrational second order perturbation theory and a variational procedure of reduced dimensionality. New subroutines suitable for the study of large systems with more than two interacting large amplitude motions were implemented and applied. The molecular symmetry of ethylene glycol decays by the formation of weak intramolecular bonds producing very asymmetrical stable structures. Three internal rotations contribute to the formation of a very anisotropic potential energy surface and to the puzzling distribution of the rovibrational energy levels. The ground vibrational state rotational constants and the centrifugal distortion constants (S-reduction, Ir representation) corresponding to the aGg' (G1) and gGg' (G2) forms are provided for the studied isotopologues. The low-lying vibrational levels up to 550 cm-1 are obtained variationally for J = 0. Two series of sublevels of the ground vibrational state are obtained: eight components localized in G1 lying between 0.0 and 0.3 cm-1 and eight sublevels localized in G2 lying between 138.1 and 138.4 cm-1. The gap between both sets is lower in CH2OD-CH2OD and more dispersed in the monodeuterated variety.

Understanding propyl cyanide and its isomers formation: ab initio study of the spectroscopy and reaction mechanisms.

Recent detection of propyl cyanide (C3H7CN) with both linear and branched structures has stimulated many experimental and theoretical studies. In this theoretical work, we present the spectroscopic properties of the far infrared spectra of these species and we investigate their different paths of formation in the gas phase. Our spectroscopic study concerns the far infrared spectra of iso, anti and gauche propyl cyanide isomers. The equilibrium structures and the potential energy surfaces are calculated using explicitly correlated cluster ab initio methods (CCSD(T)-F12) and a variational procedure designed for non-rigid species and large amplitude motions. Accurate rotational constants, centrifugal distortion constants, potential energy barriers and surfaces are provided. The rovibrational parameters in the ground vibrational states compare very well with experimental data. The low energy vibrational levels correspond to torsional modes. Far infrared energies are calculated up to 500 cm-1 using the variational approach and the vibrational second order theory (VPT2), and a good agreement with previous experimental values is found. We have also investigated the gas phase formation of the different C3H7CN isomers. After several trials of reacting gaseous species, we considered that a possible formation route of the C3H7CN isomers can be from the bimolecular reaction of HCN with propene. At the UMP2(full)/aug-cc-pVTZ level of theory, this reaction involves two steps for each isomer; the first one corresponds to the association of the two radicals while the second one corresponds to H transfer. From highly correlated ab initio calculations by means of CCSD(T)/aug-cc-pVTZ//UMP2(full)/aug-cc-pVTZ, the geometries, energetics and minimum energy paths of the reactions are obtained. Also, the first step's transition state disappears, as the diradical minimum is much less stabilized with UCCSD(T) than by UMP2(full). We employ the zero curvature tunneling and canonical variational (CVT/ZCT) semiclassical method to predict rate constants for propyl cyanide isomers formation in the gas phase. However, due to the presence of a significant barrier, this reaction leads to very small rate constants. Very recently, we probed reaction mechanisms involving radical addition and we found that these reactions are barrierless and highly exothermic leading to expected fast reactive processes for the formation of propyl cyanide products. The kinetics of such diradical reactions is under study.

A new ab initio potential energy surface for the collisional excitation of N2H(+) by H2.

We compute a new potential energy surface (PES) for the study of the inelastic collisions between N2H(+) and H2 molecules. A preliminary study of the reactivity of N2H(+) with H2 shows that neglecting reactive channels in collisional excitation studies is certainly valid at low temperatures. The four dimensional (4D) N2H(+)-H2 PES is obtained from electronic structure calculations using the coupled cluster with single, double, and perturbative triple excitation level of theory. The atoms are described by the augmented correlation consistent triple zeta basis set. Both molecules were treated as rigid rotors. The potential energy surface exhibits a well depth of ≃2530 cm(-1). Considering this very deep well, it appears that converged scattering calculations that take into account the rotational structure of both N2H(+) and H2 should be very difficult to carry out. To overcome this difficulty, the "adiabatic-hindered-rotor" treatment, which allows para-H2(j = 0) to be treated as if it were spherical, was used in order to reduce the scattering calculations to a 2D problem. The validity of this approach is checked and we find that cross sections and rate coefficients computed from the adiabatic reduced surface are in very good agreement with the full 4D calculations.

The Far-Infrared Spectrum of Methoxymethanol (CH3-O-CH2OH): A Theoretical Study.

Methoxymethanol (CH3OCH2OH), an oxygenated volatile organic compound of low stability detected in the interstellar medium, represents an example of nonrigid organic molecules showing various interacting and inseparable large-amplitude motions. The species discloses a relevant coupling among torsional modes, strong enough to prevent complete assignments using effective Hamiltonians of reduced dimensionality. Theoretical models for rotational spectroscopy can improve if they treat three vibrational coordinates together. In this paper, the nonrigid properties and the far-infrared region are analyzed using highly correlated ab initio methods and a three-dimensional vibrational model. The molecule displays two gauche-gauche (CGcg and CGcg') and one trans-gauche (Tcg) conformers, whose relative energies are very small (CGcg/CGcg'/Tcg = 0.0:641.5:792.7 cm-1). The minima are separated by relatively low barriers (1200-1500 cm-1), and the corresponding methyl torsional barriers V 3 are estimated to be 595.7, 829.0, and 683.7 cm-1, respectively. The ground vibrational state rotational constants of the most stable geometry have been computed to be A 0 = 17233.99 MHz, B 0 = 5572.58 MHz, and C 0 = 4815.55 MHz, at ΔA 0 = -3.96 MHz, ΔB 0 = 4.76 MHz, and ΔC 0 = 2.51 MHz from previous experimental data. Low-energy levels and their tunneling splittings are determined variationally up to 700 cm-1. The A/E splitting of the ground vibrational state was computed to be 0.003 cm-1, as was expected given the methyl torsional barrier (∼600 cm-1). The fundamental levels (100), (010), and (001) are predicted at 132.133 and 132.086 cm-1 (methyl torsion), 186.507 and 186.467 cm-1 (O-CH3 torsion), and 371.925 and 371.950 cm-1 (OH torsion), respectively.

Ab Initio Study of the Large Amplitude Motions of Various Monosubstituted Isotopologues of Methylamine (CH3-NH2).

CCSD(T)-F12 theory is applied to determine electronic ground state spectroscopic parameters of various isotopologues of methylamine (CH3-NH2) containing cosmological abundant elements, such as D, 13C and 15N. Special attention is given to the far infrared region. The studied isotopologues can be classified in the G12, G6 and G4 molecular symmetry groups. The rotational and centrifugal distortion constants and the anharmonic fundamentals are determined using second order perturbation theory. Fermi displacements of the vibrational bands are predicted. The low vibrational energy levels corresponding to the large amplitude motions are determine variationally using a flexible three-dimensional model depending on the NH2 bending and wagging and the CH3 torsional coordinates. The model has been defined assuming that, in the amine group, the bending and the wagging modes interact strongly. The vibrational levels split into six components corresponding to the six minima of the potential energy surface. The accuracy of the kinetic energy parameters has an important effect on the energies. Strong interactions among the large amplitude motions are observed. Isotopic effects are relevant for the deuterated species.