A new chapter in the never ending story of cycloadditions: The puzzling case of SO2 and acetylene

A comprehensive study of the different classes of cycloaddition reactions ([3+2], [2+2], and [2+1]) of SO2 to acetylene and ethylene has been performed using density functional theory (DFT) and composite wavefunction methods. The [3+2] cycloaddition reaction, that was previously explored in the context of the cycloaddition of thioformaldehyde S-methylide (TSM) to ethylene and acetylene, proceeds in a concerted way to the formation of stable heterocycles. In this paper, we extend our study to the [2+2] and [2+1] cycloadditions of SO2 to acetylene, which would produce 1,1-oxathiete-2-oxide and thiirene-1,1-dioxide, respectively. One of the main conclusions is that cyclic 1,1-oxathiete-2-oxide can open through a relatively easy breaking of the SO single bond and rearrange toward sulfinyl acetaldehyde (SA). The SA molecule can easily undergo several internal rearrangements, which eventually lead to sulfenic acid and sulfoxide derivatives of ethenone, 1,2,3-dioxathiole, and CO plus sulfinylmethane. The most probable path, however, produces 2-thioxoacetic acid, whose derivatives (or those of the corresponding acetate) are usually obtained by Willgerodt-Kindler-type sulfuration of acetates. This product can in turn decompose, leading to the final products CO2 and H2CS. Comparison of this decomposition path with that of 2-amino-2-thioxoacetic acid shows that the process occurs through different H-transfer processes.

Toward Spectroscopic Accuracy for the Structures of Large Molecules at DFT Cost: Refinement and Extension of the Nano-LEGO Approach

The SE100 database collecting accurate equilibrium geometries of medium size molecules obtained by the semiexperimental (SE) approach has been extended to species containing Br and I atoms. This has allowed the determination of accurate linear regressions between DFT and SE values for all the main bonds and angles involving H, B, C, N, O, F, P, S, Cl, Br, and I atoms. An improved Nano-LEGO tool has been developed, which is based on suitable hybrid and double hybrid functionals and combines in a fully coherent way the templating molecule and linear regression approaches. A number of case studies show that the new Nano LEGO tool provides geometrical parameters on par with state-of-the-art composite wave function methods, but can be routinely applied to medium- to large-size molecules. The accuracy reached for structural parameters is mirrored on rotational constants that can be predicted with an average error within 0.2%.

Reliable Gas Phase Reaction Rates at Affordable Cost by Means of the Parameter-Free JunChS-F12 Model Chemistry

A recently developed strategy for the computation at affordable cost of reliable barrier heights ruling reactions in the gas phase (junChS, [Barone, V.; J. Chem. Theory Comput. 2021, 17, 4913-4928]) has been extended to the employment of explicitly correlated (F12) methods. A thorough benchmark based on a wide range of prototypical reactions shows that the new model (referred to as junChS-F12), which employs cost-effective revDSD-PBEP86-D3(BJ) reference geometries, has an improved performance with respect to its conventional counterpart and outperforms the most well-known model chemistries without the need of any empirical parameter and at an affordable computational cost. Several benchmarks show that revDSD-PBEP86-D3(BJ) structures and force fields provide zero point energies and thermal contributions, which can be confidently used, together with junChS-F12 electronic energies, for obtaining accurate reaction rates in the framework of the master equation approach based on the ab initio transition-state theory.

Correcting the Experimental Enthalpies of Formation of Some Members of the Biologically Significant Sulfenic Acids Family

Sulfenic acids are important intermediates in the oxidation of cysteine thiol groups in proteins by reactive oxygen species. The mechanism is influenced heavily by the presence of polar groups, other thiol groups, and solvent, all of which determines the need to compute precisely the energies involved in the process. Surprisingly, very scarce experimental information exists about a very basic property of sulfenic acids, the enthalpies of formation. In this Article, we use high level quantum chemical methods to derive the enthalpy of formation at 298.15 K of methane-, ethene-, ethyne-, and benzenesulfenic acids, the only ones for which some experimental information exists. The methods employed were tested against well-known experimental data of related species and extensive CCSD(T) calculations. Our best results consistently point out to a much lower enthalpy of formation of methanesulfenic acid, CH₃SOH (Δ_fH⁰(298.15K) = -35.1 ± 0.4 kcal mol^-1), than the one reported in the NIST thermochemical data tables. The enthalpies of formation derived for ethynesulfenic acid, HC≡CSOH, +32.9 ± 1.0 kcal/mol, and benzenesulfenic acid, C₆H₅SOH, -2.6 ± 0.6 kcal mol^-1, also differ markedly from the experimental values, while the enthalpy of formation of ethenesulfenic acid CH₂CHSOH, not available experimentally, was calculated as -11.2 ± 0.7 kcal mol^-1.

In Vitro and In Silico Vibrational-Rotational Spectroscopic Characterization of the Next-Generation Refrigerant HFO-1123

Very short-lived substances have recently been proposed as replacements for hydrofluorocarbons (HFCs), in turn being used in place of ozone-depleting substances, in refrigerant applications. In this respect, hydro-fluoro-olefins (HFOs) are attracting particular interest because, due to their reduced global warming potential, they are supposed to be environmentally friendlier. Notwithstanding this feature, they represent a new class of compounds whose spectroscopic properties and reactivity need to be characterized to allow their atmospheric monitoring and to understand their environmental fate. In the present work, the structural, vibrational, and ro-vibrational properties of trifluorothene (HFO-1123, F₂C = CHF) are studied by state-of-the-art quantum chemical calculations. The equilibrium molecular structure has an expected error within 2 mÅ and 0.2° for bond lengths and angles, respectively. This represents the first step toward the computation of highly accurate rotational constants for both the ground and first excited fundamental vibrational levels, which reproduce the available experimental data well within 0.1%. Centrifugal distortion parameters and vibrational–rotational coupling terms are computed as well and used to solve some conflicting experimental results. Simulation of the vibrational transition frequencies and intensities beyond the double harmonic approximation and up to three quanta of vibrational excitation provides insights into the couplings ruling the vibrational dynamics and guides the characterization of the gas-phase infrared spectrum experimentally recorded in the range of 200–5000 cm^–1. The full characterization of the IR features is completed with the experimental determination of the absorption cross sections over the 400–5000 cm^–1 region from which the radiative forcing and global warming potential of HFO-1123 are derived.

Dipolar 1,3-cycloaddition of thioformaldehyde S-methylide (CH2 SCH2 ) to ethylene and acetylene. A comparison with (valence) isoelectronic O3 , SO2 , CH2 OO and CH2 SO

Methods rooted in the density functional theory and in the coupled cluster ansatz were employed to investigate the cycloaddition reactions to ethylene and acetylene of 1,3-dipolar species including ozone and the derivatives issued from replacement of the central oxygen atom by the valence-isoelectronic sulfur atom, and/or of one or both terminal oxygen atoms by the isoelectronic CH₂ group. This gives rise to five different 1,3-dipolar compounds, namely ozone itself (O₃ ), sulfur dioxide (SO₂ ), the simplest Criegee intermediate (CH₂ OO), sulfine (CH₂ SO), and thioformaldehyde S-methylide (CH₂ SCH₂ , TSM). The experimental and accurate theoretical data available for some of those molecules were employed to assess the accuracy of two last-generation composite methods employing conventional or explicitly correlated post-Hartree-Fock contributions (jun-Cheap and SVECV-f12, respectively), which were then applied to investigate the reactivity of TSM. The energy barriers provided by both composite methods are very close (the average values for the two composite methods are 7.1 and 8.3 kcal mol^-1 for the addition to ethylene and acetylene, respectively) and comparable to those ruling the corresponding additions of ozone (4.0 and 7.7 kcal mol^-1 , respectively). These and other evidences strongly suggest that, at least in the case of cycloadditions, the reactivity of TSM is similar to that of O₃ and very different from that of SO₂ .

Gliding on Ice in Search of Accurate and Cost-Effective Computational Methods for Astrochemistry on Grains: The Puzzling Case of the HCN Isomerization

The isomerization of hydrogen cyanide to hydrogen isocyanide on icy grain surfaces is investigated by an accurate composite method (jun-Cheap) rooted in the coupled cluster ansatz and by density functional approaches. After benchmarking density functional predictions of both geometries and reaction energies against jun-Cheap results for the relatively small model system HCN···(H2O)2, the best performing DFT methods are selected. A large cluster containing 20 water molecules is then employed within a QM/QM' approach to include a realistic environment mimicking the surface of icy grains. Our results indicate that four water molecules are directly involved in a proton relay mechanism, which strongly reduces the activation energy with respect to the direct hydrogen transfer occurring in the isolated molecule. Further extension of the size of the cluster up to 192 water molecules in the framework of a three-layer QM/QM'/MM model has a negligible effect on the energy barrier ruling the isomerization. Computation of reaction rates by the transition state theory indicates that on icy surfaces, the isomerization of HNC to HCN could occur quite easily even at low temperatures thanks to the reduced activation energy that can be effectively overcome by tunneling.

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 (CH₂OHCOOH), 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.

A Computational Journey across Nitroxide Radicals: From Structure to Spectroscopic Properties and Beyond

Nitroxide radicals are characterized by a long-lived open-shell electronic ground state and are strongly sensitive to the chemical environment, thus representing ideal spin probes and spin labels for paramagnetic biomolecules and materials. However, the interpretation of spectroscopic parameters in structural and dynamic terms requires the aid of accurate quantum chemical computations. In this paper we validate a computational model rooted into double-hybrid functionals and second order vibrational perturbation theory. Then, we provide reference quantum chemical results for the structures, vibrational frequencies and other spectroscopic features of a large panel of nitroxides of current biological and/or technological interest.

Accurate Biomolecular Structures by the Nano-LEGO Approach: Pick the Bricks and Build Your Geometry

The determination of accurate equilibrium molecular structures plays a fundamental role for understanding many physical–chemical properties of molecules, ranging from the precise evaluation of the electronic structure to the analysis of the role played by dynamical and environmental effects in tuning their overall behavior. For small semi-rigid systems in the gas phase, state-of-the-art quantum chemical computations rival the most sophisticated experimental (from, for example, high-resolution spectroscopy) results. For larger molecules, more effective computational approaches must be devised. To this end, we have further enlarged the compilation of available semi-experimental (SE) equilibrium structures, now covering the most important fragments containing H, B, C, N, O, F, P, S, and Cl atoms collected in the new SE100 database. Next, comparison with geometries optimized by methods rooted in the density functional theory showed that the already remarkable results delivered by PW6B95 and, especially, rev-DSDPBEP86 functionals can be further improved by a linear regression (LR) approach. Use of template fragments (taken from the SE100 library) together with LR estimates for the missing interfragment parameters paves the route toward accurate structures of large molecules, as witnessed by the very small deviations between computed and experimental rotational constants. The whole approach has been implemented in a user-friendly tool, termed nano-LEGO, and applied to a number of demanding case studies.

Development and Validation of a Parameter-Free Model Chemistry for the Computation of Reliable Reaction Rates

A recently developed model chemistry (jun-Cheap) has been slightly modified and proposed as an effective, reliable, and parameter-free scheme for the computation of accurate reaction rates with special reference to astrochemical and atmospheric processes. Benchmarks with different sets of state-of-the-art energy barriers spanning a wide range of values show that, in the absence of strong multireference contributions, the proposed model outperforms the most well-known model chemistries, reaching a subchemical accuracy without any empirical parameter and with affordable computer times. Some test cases show that geometries, energy barriers, zero point energies, and thermal contributions computed at this level can be used in the framework of the master equation approach based on the ab initio transition-state theory for obtaining accurate reaction rates.

A computational insight into the relationship between side chain IR line shapes and local environment in fibril-like structures

Infrared spectroscopy is a widely used technique to characterize protein structures and protein mediated processes. While the amide I band provides information on proteins' secondary structure, amino acid side chains are used as infrared probes for the investigation of protein reactions and local properties. In this paper, we use a hybrid quantum mechanical/classical molecular dynamical approach based on the perturbed matrix method to compute the infrared band due to the C=O stretching mode of amide-containing side chains. We calculate, at first, the infrared band of zwitterionic glutamine in water and obtain results in very good agreement with the experimental data. Then, we compute the signal arising from glutamine side chains in a microcrystal of the yeast prion Sup35-derived peptide, GNNQQNY, with a fibrillar structure. The infrared bands obtained by selective isotopic labeling of the two glutamine residues, Q4 and Q5, of each peptide were experimentally used to investigate the local hydration in the fibrillar microcrystal. The experimental spectra of the two glutamine residues, which experience different hydration environments, feature different spectral signals that are well reproduced by the corresponding calculated spectra. In addition, the analysis of the simulated spectra clarifies the molecular origin of the experimentally observed spectroscopic differences that arise from the different local electric field experienced by the two glutamine residues, which is, in turn, determined by a different hydrogen bonding pattern.

Looking for the Elusive Imine Tautomer of Creatinine: Different States of Aggregation Studied by Quantum Chemistry and Molecular Spectroscopy

New spectroscopic experiments and state-of-the-art quantum-chemical computations of creatinine in different aggregation states unequivocally unveiled a significant tuning of tautomeric equilibrium by the environment: from the exclusive presence of the amine tautomer in the solid state and aqueous solution to a mixture of amine and imine tautomers in the gas phase. Quantum-chemical calculations predict the amine species as the most stable tautomer by about 30 kJ mol^-1 in condensed phases. On the contrary, moving to the isolated forms, both Z and E imine isomers become more stable by about 7 kJ mol^-1 . Since the imine isomers and one amine tautomer are separated by significant energy barriers, all of them should be present in the gas phase. This prediction has indeed been confirmed by high-resolution rotational spectroscopy, which provides the first experimental characterization of the elusive imine tautomer. The interpretation of the complicated hyperfine structure of the rotational spectrum, originated by three ¹⁴ N nuclei, makes it possible to use the spectral signatures as a sort of fingerprint for each individual tautomer in the complex sample.

Isomerization and Fragmentation Reactions on the [C2SH4] Potential Energy Surface: The Metastable Thione S-Methylide Isomer

Thione S-methylide, parent species of the thiocarbonyl ylide family, is a 1,3-dipolar species on the [C₂SH₄] potential energy surface, not so much studied as its isomers, thiirane, vinyl thiol, and thioacetaldehyde. The conrotatory ring-closure reaction toward thiirane was studied in the 90s, but no complete analysis of the potential energy surface is available. In this paper, we report a computational study of the reaction scheme linking all species. We employed several computational methods (density functional theory, CCSD(T) composite schemes, and CASSCF/CASPT2 multireference procedures) to find the best description of thione S-methylide, its isomers, and transition states. The barrier from thiirane to thione S-methylide amounts to 52.2 kcal mol^-1 (against 17.6 kcal mol^-1 for the direct one), explaining why thiocarbonyl ylides cannot be prepared from thiiranes. Conversion of thiirane to vinyl thiol implies a large barrier, supporting why the reaction has been observed only at high temperatures. Fragmentations of thiirane to S(³P) or S(¹D) and ethylene as well as decomposition to hydrogen sulfide plus acetylene were also explored. Triplet and singlet open-shell species were identified as intermediates in the fragmentations, with energies lower than the transition state between thiirane and vinyl thiol, explaining the preference of the latter at low temperatures.

Accuracy Meets Interpretability for Computational Spectroscopy by Means of Hybrid and Double-Hybrid Functionals

Accuracy and interpretability are often seen as the devil and holy grail in computational spectroscopy and their reconciliation remains a primary research goal. In the last few decades, density functional theory has revolutionized the situation, paving the way to reliable yet effective models for medium size molecules, which could also be profitably used by non-specialists. In this contribution we will compare the results of some widely used hybrid and double hybrid functionals with the aim of defining the most suitable recipe for all the spectroscopic parameters of interest in rotational and vibrational spectroscopy, going beyond the rigid rotor/harmonic oscillator model. We will show that last-generation hybrid and double hybrid functionals in conjunction with partially augmented double- and triple-zeta basis sets can offer, in the framework of second order vibrational perturbation theory, a general, robust, and user-friendly tool with unprecedented accuracy for medium-size semi-rigid molecules.

Reinvestigation of the Deceptively Simple Reaction of Toluene with OH and the Fate of the Benzyl Radical: The "Hidden" Routes to Cresols and Benzaldehyde

In a previous work, we have investigated the initial steps of the reaction of toluene with the hydroxyl radical using several quantum chemical approaches including density functional and composite post-Hartree-Fock models. Comparison of H-abstraction from the methyl group and additions at different positions of the phenyl ring showed that the former reaction channel is favored at room temperature. This conclusion appears at first sight incompatible with the experimental observation of a lower abundance of the product obtained from abstraction (benzaldehyde) with respect to those originating from addition (cresols). Further reactions of the intermediate radicals with oxygen, water, and additional OH radicals are explored in this paper through theoretical calculations on more than 120 species on the corresponding potential energy surface. The study of the addition reactions, to obtain the cresols through hydroxy methylcyclodienyl intermediate radicals, showed that only in the case of o-cresol the reaction proceeds by addition of O₂ to the ring, internal H-transfer, and hydroperoxyl abstraction and not through direct H-abstraction. For both p- and m-cresol, instead, the reaction occurs through a higher-energy direct H-abstraction, thus explaining in part the observed larger concentration of the ortho isomer in the final products. It was also found that the benzyl radical, formed by H-abstraction from the methyl group, is able to react further if additional OH is present. Two reaction paths leading to o-cresol, two leading to p-cresol, and one leading to m-cresol were determined. Moreover, in this situation, the benzyl radical is predicted to produce benzyl alcohol, as was found in some experiments. The commonly accepted route to benzaldehyde was found to be not the energetically favored one. Instead, a route leading to the benzoyl radical (and ultimately to benzoic acid) with the participation of one water molecule was clearly more favorable, both thermodynamically and kinetically.

The Role of State-of-the-Art Quantum-Chemical Calculations in Astrochemistry: Formation Route and Spectroscopy of Ethanimine as a Paradigmatic Case

The gas-phase formation and spectroscopic characteristics of ethanimine have been re-investigated as a paradigmatic case illustrating the accuracy of state-of-the-art quantum-chemical (QC) methodologies in the field of astrochemistry. According to our computations, the reaction between the amidogen, NH, and ethyl, C2H5, radicals is very fast, close to the gas-kinetics limit. Although the main reaction channel under conditions typical of the interstellar medium leads to methanimine and the methyl radical, the predicted amount of the two E,Z stereoisomers of ethanimine is around 10%. State-of-the-art QC and kinetic models lead to a [E-CH3CHNH]/[Z-CH3CHNH] ratio of ca. 1.4, slightly higher than the previous computations, but still far from the value determined from astronomical observations (ca. 3). An accurate computational characterization of the molecular structure, energetics, and spectroscopic properties of the E and Z isomers of ethanimine combined with millimeter-wave measurements up to 300 GHz, allows for predicting the rotational spectrum of both isomers up to 500 GHz, thus opening the way toward new astronomical observations.

Modeling amino-acid side chain infrared spectra: the case of carboxylic residues

Infrared (IR) spectroscopy is commonly utilized for the investigation of protein structures and protein-mediated processes.

State-of-the-art computation of the rotational and IR spectra of the methyl-cyclopropyl cation: hints on its detection in space

Recent measurements by the Cassini Ion Neutral Mass Spectrometer demonstrated the presence of numerous carbocations in Titan's upper atmosphere. In [Ali et al., Planet. Space Sci., 2013, 87, 96], an analysis of these measurements revealed the formation of the three-membered cyclopropenyl cation and its methyl derivatives. As a starting point of a future coordinated effort of laboratory experiments, quantum-chemical calculations, and astronomical observations, in the present work the molecular structure and spectroscopic properties of the methyl-cyclopropenyl cation have been investigated by means of state-of-the-art computational approaches in order to simulate its rotational and infrared spectra. Rotational parameters have been predicted with an expected accuracy better than 0.1% for rotational constants and on the order of 1-2% for centrifugal-distortion terms. As for the infrared spectrum, despite the challenge of a large amplitude motion, fundamental transitions have been computed to a good accuracy, i.e., the uncertainties are expected to be smaller than 5-10 wavenumbers.

Molecular synthons for accurate structural determinations: the equilibrium geometry of 1-chloro-1-fluoroethene

The equilibrium structure for 1-chloro-1-fluoroethene is obtained by high-level quantum chemical calculations and the semi-experimental approach.

Unveiling the Sulfur-Sulfur Bridge: Accurate Structural and Energetic Characterization of a Homochalcogen Intermolecular Bond

By combining rotational spectroscopy in supersonic expansion with the capability of state-of-the-art quantum-chemical computations in accurately determining structural and energetic properties, the genuine nature of a sulfur-sulfur chalcogen bond between dimethyl sulfide and sulfur dioxide has been unveiled in a gas-jet environment free from collision, solvent and matrix perturbations. A SAPT analysis pointed out that electrostatic S⋅⋅⋅S interactions play the dominant role in determining the stability of the complex, largely overcoming dispersion and C-H⋅⋅⋅O hydrogen-bond contributions. Indeed, in agreement with the analysis of the quadrupole-coupling constants and of the methyl internal rotation barrier, the NBO and NOCV/CD approaches show a marked charge transfer between the sulfur atoms. Based on the assignment of the rotational spectra for 7 isotopologues, an accurate semi-experimental equilibrium structure for the heavy-atom backbone of the molecular complex has been determined, which is characterized by a S⋅⋅⋅S distance (2.947(3) Å) well below the sum of van der Waals radii.

Towards the SMART workflow system for computational spectroscopy

Is it possible to convert highly specialized research in the field of computational spectroscopy into robust and user-friendly aids to experiments and industrial applications?

Correct Modeling of Cisplatin: a Paradigmatic Case

Quantum chemistry is a useful tool in modern approaches to drug and material design, but only when the adopted model reflects a correct physical picture. Paradigmatic is the case of cis‐diaminodichloroplatinum(II), cis‐[Pt(NH₃)₂Cl₂], for which the correct simulation of the structural and vibrational properties measured experimentally still remains an open question. By using this molecule as a proof of concept, it is shown that state‐of‐the‐art quantum chemical calculations and a simple model, capturing the basic physical flavors, a cis‐[Pt(NH₃)₂Cl₂] dimer, can provide the accuracy required for interpretative purposes. The present outcomes have fundamental implications for benchmark studies aiming at assessing the accuracy of a given computational protocol.

VMS-ROT: A New Module of the Virtual Multifrequency Spectrometer for Simulation, Interpretation, and Fitting of Rotational Spectra

The Virtual Multifrequency Spectrometer (VMS) is a tool that aims at integrating a wide range of computational and experimental spectroscopic techniques with the final goal of disclosing the static and dynamic physical-chemical properties "hidden" in molecular spectra. VMS is composed of two parts, namely, VMS-Comp, which provides access to the latest developments in the field of computational spectroscopy, and VMS-Draw, which provides a powerful graphical user interface (GUI) for an intuitive interpretation of theoretical outcomes and a direct comparison to experiment. In the present work, we introduce VMS-ROT, a new module of VMS that has been specifically designed to deal with rotational spectroscopy. This module offers an integrated environment for the analysis of rotational spectra: from the assignment of spectral transitions to the refinement of spectroscopic parameters and the simulation of the spectrum. While bridging theoretical and experimental rotational spectroscopy, VMS-ROT is strongly integrated with quantum-chemical calculations, and it is composed of four independent, yet interacting units: (1) the computational engine for the calculation of the spectroscopic parameters that are employed as a starting point for guiding experiments and for the spectral interpretation, (2) the fitting-prediction engine for the refinement of the molecular parameters on the basis of the assigned transitions and the prediction of the rotational spectrum of the target molecule, (3) the GUI module that offers a powerful set of tools for a vis-à-vis comparison between experimental and simulated spectra, and (4) the new assignment tool for the assignment of experimental transitions in terms of quantum numbers upon comparison with the simulated ones. The implementation and the main features of VMS-ROT are presented, and the software is validated by means of selected test cases ranging from isolated molecules of different sizes to molecular complexes. VMS-ROT therefore offers an integrated environment for the analysis of the rotational spectra, with the innovative perspective of an intimate connection to quantum-chemical calculations that can be exploited at different levels of refinement, as an invaluable support and complement for experimental studies.

On the competition between weak O-H···F and C-H···F hydrogen bonds, in cooperation with C-H···O contacts, in the difluoromethane - tert-butyl alcohol cluster

The 1:1 complex of tert-butyl alcohol with difluoromethane has been characterized by means of a joint experimental-computational investigation. Its rotational spectrum has been recorded by using a pulsed-jet Fourier-Transform microwave spectrometer. The experimental work has been guided and supported by accurate quantum-chemical calculations. In particular, the computed potential energy landscape pointed out the formation of three stable isomers. However, the very low interconversion barriers explain why only one isomer, showing one O-H···F and two C-H···O weak hydrogen bonds, has been experimentally characterized. The effect of the H → tert-butyl- group substitution has been analyzed from the comparison to the difluoromethane-water adduct.

Development and Implementation of Advanced Fitting Methods for the Calculation of Accurate Molecular Structures

The determination of accurate equilibrium molecular structures plays a fundamental role for understanding many physical-chemical properties of molecules, ranging from the precise evaluation of the electronic structure to the analysis of dynamical and environmental effects in tuning their overall behavior. For this purpose the so-called semiexperimental approach, based on a nonlinear least-squares fit of the moments of inertia associated with a set of available isotopologues, allows one to obtain very accurate results, without the unfavorable computational cost characterizing high-level quantum chemical methods. In the present work the MSR (Molecular Structure Refinement) software for the determination of equilibrium structures by means of the semiexperimental approach is presented, and its implementation is discussed in some detail. The software, which is interfaced with a powerful graphical user interface, includes different optimization algorithms, an extended error analysis, and a number of advanced features, the most remarkable ones concerning the choice of internal coordinates and the method of predicate observations. In particular, a new black-box scheme for defining automatically a suitable set of nonredundant internal coordinates of A₁ symmetry in place of the customary Z-matrix has been designed and tested. Finally, the implementation of the method of the predicate observations is discussed and validated for a set of test molecules. As an original application, the method is employed for the determination of the semiexperimental structure for the most stable conformer of glycine.

Accurate Vibrational-Rotational Parameters and Infrared Intensities of 1-Bromo-1-fluoroethene: A Joint Experimental Analysis and Ab Initio Study

The medium-resolution gas-phase infrared (IR) spectra of 1-bromo-1-fluoroethene (BrFC═CH₂, 1,1-C₂H₂BrF) were investigated in the range 300-6500 cm^-1, and the vibrational analysis led to the assignment of all fundamentals as well as many overtone and combination bands up to three quanta, thus giving an accurate description of its vibrational structure. Integrated band intensity data were determined with high precision from the measurements of their corresponding absorption cross sections. The vibrational analysis was supported by high-level ab initio investigations. CCSD(T) computations accounting for extrapolation to the complete basis set and core correlation effects were employed to accurately determine the molecular structure and harmonic force field. The latter was then coupled to B2PLYP and MP2 computations in order to account for mechanical and electrical anharmonicities. Second-order perturbative vibrational theory was then applied to the thus obtained hybrid force fields to support the experimental assignment of the IR spectra.

Structural features of the carbon-sulfur chemical bond: a semi-experimental perspective

In this work semi-experimental and theoretical equilibrium geometries of 10 sulfur-containing organic molecules, as well as 4 oxygenated ones, are determined by means of a computational protocol based on density functional theory. The results collected in the present paper further enhance our online database of accurate semi-experimental equilibrium molecular geometries, adding 13 new molecules containing up to 8 atoms, for 12 of which the first semi-experimental equilibrium structure is reported, to the best of our knowledge. We focus in particular on sulfur-containing compounds, aiming both to provide new accurate data on some rather important chemical moieties, only marginally represented in the literature of the field, and to examine the structural features of carbon-sulfur bonds in the light of the previously presented linear regression approach. The structural changes issuing from substitution of oxygen by sulfur are discussed to get deeper insights on how modifications in electronic structure and nuclear potential can affect equilibrium geometries. With respect to our previous works, we perform non-linear constrained optimizations of equilibrium SE structures with a new general and user-friendly software under development in our group with updated definition of useful statistical indicators.

Probing the existence of phase transitions in one-dimensional fluids of penetrable particles

Phase transitions in one-dimensional classical fluids are usually ruled out by using van Hove's theorem. A way to circumvent the conclusions of the theorem is to consider an interparticle potential that is everywhere bounded. Such is the case of, e.g., the generalized exponential model of index 4 (GEM-4 potential), which in three dimensions gives a reasonable description of the effective repulsion between flexible dendrimers in a solution. An extensive Monte Carlo simulation of the one-dimensional GEM-4 model [S. Prestipino, Phys. Rev. E 90, 042306 (2014)] has recently provided evidence of an infinite sequence of low-temperature cluster phases, however, also suggesting that upon pushing the simulation forward what seemed a true transition may eventually prove to be only a sharp crossover. We hereby investigate this problem theoretically by use of three different and increasingly sophisticated approaches (i.e., a mean-field theory, the transfer matrix of a lattice model of clusters, and the exact treatment of a system of point clusters in the continuum) to conclude that the alleged transitions of the one-dimensional GEM-4 system are likely just crossovers.

Insights into the interaction between CH2F2 and titanium dioxide: DRIFT spectroscopy and DFT analysis of the adsorption energetics

Difluoromethane (CH2F2, HFC-32) has been proposed as a valid replacement for both CFCs and HCFCs (in particular HCFC-22), and nowadays it is widely used in refrigerant mixtures. Due to its commercial use, in the last years, the atmospheric concentration of HFC-32 has increased significantly. However, this molecule presents strong absorptions within the 8-12μm atmospheric window, and hence it is a greenhouse gas which contributes to global warming. Heterogeneous photocatalysis over TiO2 surface is an interesting technology for removing atmospheric pollutants since it leads to the decomposition of organic compounds into simpler molecules. In the present work, the adsorbate-substrate interaction between CH2F2 and TiO2 is investigated by coupling experimental measurements using DRIFT spectroscopy to first-principle simulations at DFT/B3LYP level. The experimental results confirm that CH2F2 interacts with the TiO2 surface (∼80% rutile, 20% anatase) through both F and H atoms and show that the DRIFT technique is well suited to study the adsorption of halogenated methanes over semiconductor surfaces. DFT calculations are carried out by considering different periodicities and surface coverages, according to a structure involving an acid-base interaction between the F and Ti(4+) atoms as well as an H-bond between the CH2 group and an O(2-) ion. Lateral effects and energetics are analyzed in the limit of low coverage according to a procedure taking into account the binding, interaction, and distortion energies. The simulation at the different surface coverages and periodicities suggests similar decomposition pathways for the different investigated ensemble configurations.

Unveiling the non-covalent interactions of molecular homodimers by dispersion-corrected DFT calculations and collision-induced broadening of ro-vibrational transitions: application to (CH2F2)2 and (SO2)2

(CH₂F₂)₂ and (SO₂)₂ are investigated using DFT-D3 computations, and experimental dissociation energies are determined by TDL-IR spectroscopy. DFT-D3 dramatically improves over uncorrected DFT.

N2-, O2- and He-collision-induced broadening of sulfur dioxide ro-vibrational lines in the 9.2 μm atmospheric window

Sulfur dioxide (SO2) is a molecule of considerable interest for both atmospheric chemistry and astrophysics. In the Earth's atmosphere, it enters in the sulfur cycle and it is ubiquitous present in polluted atmospheres, where it is responsible for acid rains. It is also of astrophysical and planetological importance, being present on Venus and in interstellar clouds. In this work the collisional broadening of a number of ν1 ro-vibrational lines of SO2 perturbed by N2, O2 and He are investigated at room temperature in the 9 μm atmospheric region by means of high resolution tunable diode laser (TDL) infrared spectroscopy. From N2- and O2-broadening coefficients, the broadening parameters of sulfur dioxide in air, useful for atmospheric applications, are derived as well. From the present measurements some conclusions on the quantum number dependence of the N2-, O2- and He-broadening coefficients are drawn. While the J dependence is weak for all the perturbers investigated, different trends with Ka are reported. N2-broadening coefficients show a slight decrease with increasing values of Ka, whereas O2 and He broadening cross sections first increase up to Ka(″)≈6 and then they keep a nearly constant value. A comparison and a brief discussion on the efficiency of self-, N2-, O2- and He-collisional dynamics are given. The data obtained represent a significant analysis on foreign broadening of SO2 useful for atmospheric remote sensing and astrophysical applications.

An integrated experimental and quantum-chemical investigation on the vibrational spectra of chlorofluoromethane

The vibrational analysis of the gas-phase infrared spectra of chlorofluoromethane (CH2ClF, HCFC-31) was carried out in the range 200-6200 cm(-1). The assignment of the absorption features in terms of fundamental, overtone, combination, and hot bands was performed on the medium-resolution (up to 0.2 cm(-1)) Fourier transform infrared spectra. From the absorption cross section spectra accurate values of the integrated band intensities were derived and the global warming potential of this compound was estimated, thus obtaining values of 323, 83, and 42 on a 20-, 100-, and 500-year horizon, respectively. The set of spectroscopic parameters here presented provides the basic data to model the atmospheric behavior of this greenhouse gas. In addition, the obtained vibrational properties were used to benchmark the predictions of state-of-the-art quantum-chemical computational strategies. Extrapolated complete basis set limit values for the equilibrium geometry and harmonic force field were obtained at the coupled-cluster singles and doubles level of theory augmented by a perturbative treatment of triple excitations, CCSD(T), in conjunction with a hierarchical series of correlation-consistent basis sets (cc-pVnZ, with n = T, Q, and 5), taking also into account the core-valence correlation effects and the corrections due to diffuse (aug) functions. To obtain the cubic and quartic semi-diagonal force constants, calculations employing second-order Møller-Plesset perturbation (MP2) theory, the double-hybrid density functional B2PLYP as well as CCSD(T) were performed. For all anharmonic force fields the performances of two different perturbative approaches in computing the vibrational energy levels (i.e., the generalized second order vibrational treatment, GVPT2, and the recently proposed hybrid degeneracy corrected model, HDCPT2) were evaluated and the obtained results allowed us to validate the spectroscopic predictions yielded by the HDCPT2 approach. The predictions of the deperturbed second-order perturbation approach, DVPT2, applied to the computation of infrared intensities beyond the double-harmonic approximation were compared to the accurate experimental values here determined. Anharmonic DFT and MP2 corrections to CCSD(T) intensities led to a very good agreement with the absorption cross section measurements over the whole spectral range here analysed.

Anharmonic theoretical simulations of infrared spectra of halogenated organic compounds

The recent implementation of the computation of infrared (IR) intensities beyond the double-harmonic approximation [J. Bloino and V. Barone, J. Chem. Phys. 136, 124108 (2012)] paved the route to routine calculations of infrared spectra for a wide set of molecular systems. Halogenated organic compounds represent an interesting class of molecules, from both an atmospheric and computational point of view, due to the peculiar chemical features related to the halogen atoms. In this work, we simulate the IR spectra of eight halogenated molecules (CH2F2, CHBrF2, CH2DBr, CF3Br, CH2CHF, CF2CFCl, cis-CHFCHBr, cis-CHFCHI), using two common hybrid and double-hybrid density functionals in conjunction with both double- and triple-ζ quality basis sets (SNSD and cc-pVTZ) as well as employing the coupled-cluster theory with basis sets of at least triple-ζ quality. Finally, we compare our results with available experimental spectra, with the aim of checking the accuracy and the performances of the computational approaches.