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

Computational Insight into the Nature and Strength of the π-Hole Type Chalcogen∙∙∙Chalcogen Interactions in the XO2∙∙∙CH3YCH3 Complexes (X = S, Se, Te; Y = O, S, Se, Te)

In recent years, the non-covalent interactions between chalcogen centers have aroused substantial research interest because of their potential applications in organocatalysis, materials science, drug design, biological systems, crystal engineering, and molecular recognition. However, studies on π-hole-type chalcogen∙∙∙chalcogen interactions are scarcely reported in the literature. Herein, the π-hole-type intermolecular chalcogen∙∙∙chalcogen interactions in the model complexes formed between XO2 (X = S, Se, Te) and CH3YCH3 (Y = O, S, Se, Te) were systematically studied by using quantum chemical computations. The model complexes are stabilized via one primary X∙∙∙Y chalcogen bond (ChB) and the secondary C-H∙∙∙O hydrogen bonds. The binding energies of the studied complexes are in the range of -21.6~-60.4 kJ/mol. The X∙∙∙Y distances are significantly smaller than the sum of the van der Waals radii of the corresponding two atoms. The X∙∙∙Y ChBs in all the studied complexes except for the SO2∙∙∙CH3OCH3 complex are strong in strength and display a partial covalent character revealed by conducting the quantum theory of atoms in molecules (QTAIM), a non-covalent interaction plot (NCIplot), and natural bond orbital (NBO) analyses. The symmetry-adapted perturbation theory (SAPT) analysis discloses that the X∙∙∙Y ChBs are primarily dominated by the electrostatic component.

Noncovalent Interactions in the Molecular Geometries of 4-Methylthiazole···H2O and 5-Methylthiazole···H2O Revealed by Microwave Spectroscopy

The pure rotational spectra of 4-methylthiazole···H2O and 5-methylthiazole···H2O were recorded by chirped-pulse Fourier transform microwave (CP-FTMW) spectroscopy. Each complex was generated within the rotationally cold environment of a gas sample undergoing supersonic expansion in the presence of an argon buffer gas. The spectra of five isotopologues of each complex have been measured and analyzed to determine the rotational constants, A0, B0, and C0; centrifugal distortion constants, DJ, DJK, and d1; nuclear quadrupole coupling constants, χaa(N3) and [χbb(N3) - χcc(N3)]; and parameters describing the internal rotation of the CH3 group, V3 and ∠(i,b). The experimentally deduced parameters were obtained using the XIAM and the BELGI-Cs-hyperfine code. For each complex, parameters in the molecular geometry are fitted to experimentally determined moments of inertia. DFT calculations have been performed at the ωB97X-D/aug-cc-pVQZ level in support of the experiments. Each complex contains two hydrogen bonds; a comparatively strong, primary interaction between the N of thiazole and an O-H of H2O, and a weaker, secondary interaction between O and either the hydrogen atom attached to C2 (in 5-methylthiazole···H2O) or the CH3 group attached to C4 (in 4-methylthiazole···H2O). The barrier to internal rotation of the CH3 group, V3, is slightly lower for 4-methylthiazole···H2O (XIAM result is 340.05(56) cm-1) than that for the 4-methylthiazole monomer (357.6 cm-1). This is likely to be a result of internal charge redistribution within the 4-methylthiazole subunit following its coordination by H2O. At the precision of the experiments, V3 of 5-methylthiazole···H2O (XIAM result is 325.16(38) cm-1) is not significantly different from V3 of the 5-methylthiazole monomer (332.0 cm-1).

Hyperfine-resolved spectra of HDS together with a global ro-vibrational analysis

Despite their chemical simplicity, the spectroscopic investigation of light hydrides, such as hydrogen sulfide, is challenging due to strong hyperfine interactions and/or anomalous centrifugal-distortion effects. Several hydrides have already been detected in the interstellar medium, and the list includes H2S and some of its isotopologues. Astronomical observation of isotopic species and, in particular, those bearing deuterium is important to gain insights into the evolutionary stage of astronomical objects and to shed light on interstellar chemistry. These observations require a very accurate knowledge of the rotational spectrum, which is so far limited for mono-deuterated hydrogen sulfide, HDS. To fill this gap, high-level quantum-chemical calculations and sub-Doppler measurements have been combined for the investigation of the hyperfine structure of the rotational spectrum in the millimeter- and submillimeter-wave region. In addition to the determination of accurate hyperfine parameters, these new measurements together with the available literature data allowed us to extend the centrifugal analysis using a Watson-type Hamiltonian and a Hamiltonian-independent approach based on the Measured Active Ro-Vibrational Energy Levels (MARVEL) procedure. The present study thus permits to model the rotational spectrum of HDS from the microwave to far-infrared region with great accuracy, thereby accounting for the effect of the electric and magnetic interactions due to the deuterium and hydrogen nuclei.

Local Charge-Displacement Analysis: Targeting Local Charge-Flows in Complex Intermolecular Interactions

Charge-displacement (CD) analysis has recently proven to be a simple and powerful scheme for quantitatively analyzing the profile of the charge redistribution occurring upon intermolecular interactions along a given interaction axis. However, when two molecular fragments bind through complex interactions involving multiple concurrent charge flows, ordinary CD analysis is capable of providing only an averaged picture of the related charge-flow profiles and no detailed information on each of them. In this article, we combine CD analysis with a Hirshfeld partitioning of the molecular charge redistribution for a local analysis on focused portions of the molecule, allowing for a detailed characterization of one charge flow at a time. The resulting scheme - the local charge-displacement (LCD) analysis - is tested on the intriguing case of the dimethyl sulfide (DMS)-sulfur dioxide (SO2) complex, characterized by concurrent charge flows relating to a sulfur-sulfur homo-chalcogen interaction and a pair of hydrogen bonds. The LCD scheme is then applied to the analysis of multiple hydrogen bonding in the acetic acid dimer, of base-pairing interactions in DNA, and of ambifunctional hydrogen bonding in the ammonia-pyridine complex.

Dipolar spin-spin coupling as an auxiliary tool for the structure determination of small isolated molecules

The "gold standard" for obtaining accurate equilibrium structures is the so-called semi-experimental (SE) approach, which exploits the structural information contained in rotational constants. Within the SE approach, ground-state rotational constants-accurately obtained from high-resolution spectroscopic studies-are computationally corrected in order to remove vibrational effects. The resulting SE equilibrium rotational constants for a significant set of isotopic species allow for retrieving a unique set of equilibrium bond lengths and angles for the molecule under consideration. However, in some cases, the lack of isotopic substitution hampers or even prevents a rigorous and complete structure determination. In this perspective, we introduce the use of dipolar spin-spin coupling constants as an additional source of structural information in support of the standard SE approach. As a proof-of-concept, we tested this new strategy on some prototypical species, such as water, ammonia, phosphine, and their fluorinated counterparts. Our results indicate that-even when the molecular structure can be obtained from a large set of isotopic rotational constants-the use of dipolar spin-spin coupling constants guarantees a better accuracy and reduces the correlations among the geometrical parameters. Moreover, we point out that our approach offers the possibility to fully derive the molecular structure of PF₃, a species for which any isotopic substitution is not possible.

junChS and junChS-F12 Models: Parameter-free Efficient yet Accurate Composite Schemes for Energies and Structures of Noncovalent Complexes

A recently developed model chemistry (denoted as junChS [Alessandrini, S.; et al. J. Chem. Theory Comput. 2020, 16, 988-1006]) has been extended to the employment of explicitly correlated (F12) methods. This led us to propose a family of effective, reliable, and parameter-free schemes for the computation of accurate interaction energies of molecular complexes ruled by noncovalent interactions. A thorough benchmark based on a wide range of interactions showed that the so-called junChS-F12 model, which employs cost-effective revDSD-PBEP86-D3(BJ) reference geometries, has an improved performance with respect to its conventional counterpart and outperforms well-known model chemistries. Without employing any empirical parameter and at an affordable computational cost, junChS-F12 reaches subchemical accuracy. Accurate characterizations of molecular complexes are usually limited to energetics. To take a step forward, the conventional and F12 composite schemes developed for interaction energies have been extended to structural determinations. A benchmark study demonstrated that the most effective option is to add MP2-F12 core-valence correlation corrections to fc-CCSD(T)-F12/jun-cc-pVTZ geometries without the need of recovering the basis set superposition error and the extrapolation to the complete basis set.

Charge-Flow Profiles along Curvilinear Paths: A Flexible Scheme for the Analysis of Charge Displacement upon Intermolecular Interactions

The Charge-Displacement (CD) analysis has proven to be a powerful tool for a quantitative characterization of the electron-density flow occurring upon chemical bonding along a suitably chosen interaction axis. In several classes of interesting intermolecular interactions, however, an interaction axis cannot be straightforwardly defined, and the CD analysis loses consistency and usefulness. In this article, we propose a general, flexible reformulation of the CD analysis capable of providing a quantitative view of the charge displacement along custom curvilinear paths. The new scheme naturally reduces to ordinary CD analysis if the path is chosen to be a straight line. An implementation based on a discrete sampling of the electron densities and a Voronoi space partitioning is described and shown in action on two test cases of a metal-carbonyl and a pyridine-ammonia complex.

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.

How accurate is the determination of equilibrium structures for van der Waals complexes? The dimer N2O⋯CO as an example

Plausible methods for accurate determination of equilibrium structures of intermolecular clusters have been assessed for the van der Waals dimer N₂O⋯CO. In order to assure a large initial dataset of rotational parameters, we first measured the microwave spectra of the ¹⁵N₂O⋯¹²CO and ¹⁵N₂O⋯¹³CO isotopologs, expanding previous measurements. Then, an anharmonic force field was calculated ab initio and a semi-experimental equilibrium structure was determined. The dimer structure was also calculated at the coupled-cluster level of theory using very large basis sets with diffuse functions and counterpoise correction. It was found that the contributions of the diffuse functions and the counterpoise correction are not additive and do not compensate each other although they have almost the same value but opposite signs. The semi-experimental and ab initio structures were found to be in fair agreement, with the equilibrium distance between the centers of mass of both monomers being 3.825(13) Å and the intermolecular bond length r(C⋯O) = 3.300(9) Å. In this case, the mass-dependent method did not permit us to determine reliable intermolecular parameters. The combination of experimental rotational constants and results of ab initio calculations thus proves to be very sensitive to examine the accuracy of structural determinations in intermolecular clusters, offering insight into other aggregates.

Interaction topologies of the SO chalcogen bond: the conformational equilibrium of the cyclohexanolSO2 cluster

The conformational landscape of the cyclohexanolSO2 cluster was revealed in the gas phase using chirped-pulsed broadband rotational spectroscopy and quantum chemical calculations. Four isomers stabilized by a dominant SO chalcogen bond and cooperative C-HO[double bond, length as m-dash]S and O-HO[double bond, length as m-dash]S secondary weak hydrogen bonds were observed, with a near-parallel orientation of the S[double bond, length as m-dash]O and O-H bonds. Isomers formed by equatorial-gauche cyclohexanol are more stable than the isomers containing axial cyclohexanol. The multiple conformations of cyclohexanol and the versatile binding properties of SO2, simultaneously operating as nucleophile and electrophile through its π-holes and non-bonding electrons lead to a complex conformational behavior when the cluster is formed. The long (2.64-2.85 Å) attractive SO interaction between SO2 and cyclohexanol is mainly electrostatic and the contribution of charge transfer is obvious, with an NBO analysis suggesting that the strength of the SO interaction is nearly two orders of magnitude larger than the hydrogen bonds. This study provides molecular insights into the structural and energetic characteristics that determine the formation of pre-nucleation clusters between SO2 and a volatile organic compound like cyclohexanol.

Unveiling Bifunctional Hydrogen Bonding with the Help of Quantum Chemistry: The Imidazole-Water Adduct as Test Case

The ubiquitous role of water and its amphiprotic nature call for a deeper insight into the physical-chemical properties of hydrogen-bonded complexes formed with building blocks of biomolecules. In this work, the semiexperimental (SE) approach combined with the template model (TM) protocol allowed the accurate determination of the equilibrium structure of two isomeric forms of the imidazole-water complex. In this procedure, the integration of experiment (thanks to a recent rotational spectroscopy investigation) and theory is exploited, also providing the means of assessing the reliability and accuracy of different quantum-chemical approaches. Overall, this study demonstrated the robustness of the combined SE-TM approach, which can provide accurate results using affordable quantum-chemical methods. Finally, the structural and energetic characteristics of these complexes have been examined in detail and compared with those of analogous heterocycle-water adducts, also exploiting energy decomposition analyses.

Looking for the bricks of the life in the interstellar medium: The fascinating world of astrochemistry

The discovery in the interstellar medium of molecules showing a certain degree of complexity, and in particular those with a prebiotic character, has attracted great interest. A complex chemistry takes place in space, but the processes that lead to the production of molecular species are a matter of intense discussion, the knowledge still being at a rather primitive stage. Debate on the origins of interstellar molecules has been further stimulated by the identification of biomolecular building blocks, such as nucleobases and amino acids, in meteorites and comets. Since many of the molecules found in space play a role in the chemistry of life, the issue of their molecular genesis and evolution might be related to the profound question of the origin of life itself. Understanding the underlying chemical processes, including the production, reactions and destruction of compounds, requires the concomitant study of spectroscopy, gas-phase reactivity, and heterogeneous processes on dust-grains. The aim of this contribution is to provide a general view of a complex and multifaceted challenge, while focusing on the role played by molecular spectroscopy and quantum-chemical computations. In particular, the derivation of the molecular spectroscopic features and the investigation of gas-phase formation routes of prebiotic species in the interstellar medium are addressed from a computational point of view.

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.

The challenge of non-covalent interactions: theory meets experiment for reconciling accuracy and interpretation

In the past decade, many gas-phase spectroscopic investigations have focused on the understanding of the nature of weak interactions in model systems. Despite the fact that non-covalent interactions play a key role in several biological and technological processes, their characterization and interpretation are still far from being satisfactory. In this connection, integrated experimental and computational investigations can play an invaluable role. Indeed, a number of different issues relevant to unraveling the properties of bulk or solvated systems can be addressed from experimental investigations on molecular complexes. Focusing on the interaction of biological model systems with solvent molecules (e.g., water), since the hydration of the biomolecules controls their structure and mechanism of action, the study of the molecular properties of hydrated systems containing a limited number of water molecules (microsolvation) is the basis for understanding the solvation process and how structure and reactivity vary from gas phase to solution. Although hydrogen bonding is probably the most widespread interaction in nature, other emerging classes, such as halogen, chalcogen and pnicogen interactions, have attracted much attention because of the role they play in different fields. Their understanding requires, first of all, the characterization of the directionality, strength, and nature of such interactions as well as a comprehensive analysis of their competition with other non-covalent bonds. In this review, it is shown how state-of-the-art quantum-chemical computations combined with rotational spectroscopy allow for fully characterizing intermolecular interactions taking place in molecular complexes from both structural and energetic points of view. The transition from bi-molecular complex to microsolvation and then to condensed phase is shortly addressed.

The preferred conformation of the tetrafluoro-1,3-dithietane⋯isopropylamine complex as revealed by rotational spectroscopy

The favored conformation of the C₂F₄S₂–IPA complex is determined by the strength of the S⋯N ChB as revealed by rotational spectroscopy.

Chemical bonding in cuprous complexes with simple nitriles: octet rule and resonance concepts versus quantitative charge-redistribution analysis

Resonance structures for six cuprous complexes with simple nitriles are interpreted by means of a quantitative analysis of charge redistribution upon copper-nitrile bonding.

Theory meets experiment for elucidating the structure and stability of non-covalent complexes: water–amine interaction as a proof of concept

A joint experimental-theoretical spectroscopic investigation has focused on a better understanding of the nature of weak, non-covalent interactions in amine-water model systems.

Xe⋯OCS: relatively straightforward?

Spectroscopy meets theory in a study of Xe⋯OCS complex: accurate near-equilibrium structures, experimental interaction energies, and CCSD(T)/CBS results presented.

Conformational Equilibria and Molecular Structures of Model Sulfur-Sulfur Bridge Systems: Diisopropyl Disulfide

The conformations and molecular structures of diisopropyl disulfide have been studied by high-resolution microwave spectroscopy and quantum chemical calculations. Three conformers, G'GG', G'GT, and GGG', have been observed in the jet expansion. The global minimum, G'GG', adopts a configuration with the G' orientation of H-C-S-S and S-S-C-H and the G orientation of C-S-S-C showing the C₂ symmetry. The rotational spectra of monosubstituted ¹³C and ³⁴S isotopologues have also been recorded for G'GG', leading to an accurate structural determination of this conformer. Two additional ³⁴S isotopologues have also been measured for G'GT. The relative energies of three observed conformers calculated at the MP2/6-311++(d,p) level of theory are within 2 kJ mol^-1, while the relative intensity measurements suggested their population ratio to be N_G'GG'/N_G'GT/N_GGG' ≈ 5:3:2.

Chalcogen bond and internal dynamics of the 2,2,4,4-tetrafluoro-1,3-dithietane⋯water complex

The rotational spectrum of the 2,2,4,4-tetrafluoro-1,3-dithietane⋯water complex has been investigated by high resolution rotational spectroscopy. Inversion of the water around its C₂ axis is hindered by a barrier determined to be 87.4(2) cm⁻¹.

Conformational preference determined by inequivalent n-pairs: rotational studies on acetophenone and its monohydrate

Acetophenone and its complex with water have been investigated by using pulsed jet Fourier transform microwave spectroscopy complemented with quantum chemical calculations.

Rotational characterization of S⋯F chalcogen bonds in the complex of 2,2,4,4-tetrafluoro-1,3-dithietane and difluoromethane

The nature of S⋯F chalcogen bonds and C–H⋯F and C–F⋯F–C contacts was characterized by rotational spectroscopy for the first time.