Rotational Spectroscopy and Conformational Flexibility of 2-Phenylethanethiol: The Dominant S-H⋠⋠⋠π Intramolecular Hydrogen Bond.

We present a rotational-computational investigation of the aromatic mercaptan 2-phenylethanethiol, addressing its potential energy surface, conformational equilibrium, internal dynamics and intramolecular interactions. The experiment used broadband chirped-pulse Fourier transform microwave spectroscopy in a supersonic jet expansion, recording the rotational spectrum in the 2-8 GHz frequency region. Two different conformers were detected in the spectrum. The most intense transitions correspond to a skew (gauche-gauche) conformation, identified as the global minimum. The spectra of ten different isotopologues were assigned for this species, leading to accurate effective and substitution structures. The weaker spectrum presents small tunnelling doublings caused by the torsional motion of the thiol group, which are only compatible with an antiperiplanar skeleton and a gauche thiol. The larger stability of the global minimum is attributed to an intramolecular S-H⋅⋅⋅π weak hydrogen bond. A comparison of the intramolecular interactions in the title molecule and 2-phenylethanol, similarly stabilized by a O-H⋅⋅⋅π hydrogen bond, shows the different strength of these interactions. Density functional (B3LYP-D3, B2PLYP-D3) and ab initio (MP2) calculations were conducted for the molecule.

Iron Complexes as Potential Carriers of Diffuse Interstellar Bands: The Photodissociation Spectrum of Fe+(H2O) at Optical Wavelengths.

The fullerene ion C60+ is the only carrier of diffuse interstellar bands (DIBs) identified so far. Transition-metal compounds feature electronic transitions in the visible and near-infrared regions, making them potential DIB carriers. Since iron is the most abundant transition metal in the cosmos, we here test this idea with Fe+(H2O). Laboratory spectra were obtained by photodissociation spectroscopy at 80 K. Spectra were modeled with the reflection principle. A high-resolution spectrum of the DIB standard star HD 183143 served as an observational reference. Two broad bands were observed from 4120 to 6800 Å. The 4120-4800 Å band has sharp features emerging from the background, which have the width of DIBs but do not match the band positions of the reference spectrum. Calculations show that the spectrum arises from a d-d transition at the iron center. While no match was found for Fe+(H2O) with known DIBs, the observation of structured bands with line widths typical for DIBs shows that small molecules or molecular ions containing iron are promising candidates for DIB carriers.

Sulfur-arene interactions: the S⋯π and S-H⋯π interactions in the dimers of benzofuran⋯sulfur dioxide and benzofuran⋯hydrogen sulfide.

Non-covalent interactions between sulfur centers and aromatic rings play important roles in biological chemistry. We examined here the sulfur-arene interactions between the fused aromatic heterocycle benzofuran and two prototype sulfur divalent triatomics (sulfur dioxide and hydrogen sulfide). The weakly-bound adducts were generated in a supersonic jet expansion and characterized with broadband (chirped-pulsed) time-domain microwave spectroscopy. The rotational spectrum confirmed the detection of a single isomer for both heterodimers, consistent with the computational predictions for the global minima. The benzofuran⋯sulfur dioxide dimer exhibits a stacked structure with sulfur closer to benzofuran, while in benzofuran⋯hydrogen sulfide the two S-H bonds are oriented towards the bicycle. These binding topologies are similar to the corresponding benzene adducts, but offer increased interaction energies. The stabilizing interactions are described as S⋯π or S-H⋯π, respectively, using a combination of density-functional theory calculations (dispersion corrected B3LYP and B2PLYP), natural bond orbital theory, energy decomposition and electronic density analysis methods. The two heterodimers present a larger dispersion component, but nearly balanced by electrostatic contributions.

Equilibrium structures of selenium compounds: The torsionally flexible molecule of selenophenol.

The equilibrium structure of selenophenol has been investigated using rotational spectroscopy and high-level quantum mechanical calculations, offering electronic and structural insight into the scarcely studied selenium compounds. The jet-cooled broadband microwave spectrum was measured in the 2-8 GHz cm-wave region using broadband (chirped-pulse) fast-passage techniques. Additional measurements up to 18 GHz used narrow-band impulse excitation. Spectral signatures were obtained for six isotopic species of selenium (80Se, 78Se, 76Se, 82Se, 77Se, and 74Se), together with different monosubstituted 13C species. The (unsplit) rotational transitions associated with the non-inverting µa-dipole selection rules could be partially reproduced with a semirigid rotor model. However, the internal rotation barrier of the selenol group splits the vibrational ground state into two subtorsional levels, doubling the dipole-inverting µb transitions. The simulation of the double-minimum internal rotation gives a very low barrier height (B3PW91: 42 cm-1), much smaller than for thiophenol (277 cm-1). A monodimensional Hamiltonian then predicts a huge vibrational separation of 72.2 GHz, justifying the non-observation of µb transitions in our frequency range. The experimental rotational parameters were compared with different MP2 and density functional theory calculations. The equilibrium structure was determined using several high-level ab initio calculations. A final Born-Oppenheimer (reBO) structure was obtained at the coupled-cluster CCSD(T)_ae/cc-wCVTZ level of theory, including small corrections for the wCVTZ → wCVQZ basis set enlargement calculated at the MP2 level. The mass-dependent method with predicates was used to produce an alternative rm(2) structure. The comparison between the two methods confirms the high accuracy of the reBO structure and offers information on other chalcogen-containing molecules.

Microhydration of Phenyl Formate: Gas-Phase Laser Spectroscopy, Microwave Spectroscopy, and Quantum Chemistry Calculations.

Herein, we have investigated the structure of phenyl formate⋅⋅⋅water (PhOF⋅⋅⋅H2 O) dimer and various non-covalent interactions present there using gas-phase laser spectroscopy and microwave spectroscopy combined with quantum chemistry calculations. Two conformers of PhOF⋅⋅⋅H2 O (C1 and T1), built on the two cis/trans conformers of the bare molecule, have been observed in the experiment. In cis-PhOF, there is an nCO → π A r * ${{{\rm \pi }}_{{\rm A}{\rm r}}^{{\rm {^\ast}}}}$ interaction between the lone-pair orbital of the carbonyl oxygen atom and the π* orbital of the phenyl ring, which persists in the monohydrated C1 conformer of PhOF⋅⋅⋅H2 O according to the NBO and NCI analyses. On the other hand, this interaction is absent in the trans-PhOF conformer as the C=O group is away from the phenyl ring. The C1 conformer is primarily stabilized by an interplay between O-H⋅⋅⋅O=C hydrogen bond and O-H⋅⋅⋅π interactions, while the stability of the T1 conformer is primarily governed by the O-H⋅⋅⋅O=C hydrogen bond. The most important finding of the present work is that the conformational preference of the PhOF monomer is retained in its monohydrated complex.

Torsional chirality and molecular recognition: the homo and heterochiral dimers of thenyl and furfuryl alcohol.

Furfuryl alcohol and thenyl alcohol contain a labile torsional chiral center, producing transiently chiral enantiomers interconverting in the nanosecond time-scale. We explored chiral molecular recognition using the weakly-bound intermolecular dimers of both alcohols, freezing stereomutation. Supersonic jet broadband microwave spectroscopy revealed homo and heterochiral diastereoisomers for each alcohol dimer and the structural characteristics of the clusters. All dimers are primarily stabilized by a moderately intense O-H⋯O hydrogen bond, but differ in the secondary interactions, which introduce additional hydrogen bonds either to the ring oxygen in furfuryl alcohol or to the π ring system in thenyl alcohol. Density-functional calculations (B2PLYP-D3(BJ)/def2-TZVP) show no clear preferences for a particular stereochemistry in the dimers, with relative energies of the order 1-2 kJ mol-1. The study suggests opportunities for the investigation of chiral recognition in molecules with torsional barriers in between transient and permanent interconversion regimes.

The evolution towards cyclic structures in the aggregation of aromatic alcohols: the dimer, trimer and tetramer of 2-phenylethanol.

Gas-phase spectroscopic studies of alcohol clusters offer accurate information on the influence of non-covalent interactions on molecular recognition, and are of paramount importance to model supramolecular and biological chemical processes. Here, we examine the role of the aliphatic side chain in the self-aggregation of aromatic alcohols, using a multi-methodological gas-phase approach which combines microwave spectroscopy and mass-resolved electronic and vibrational laser spectroscopy. Spectroscopic and electronic structure computations were carried out for the dimer, trimer and tetramer of 2-phenylethanol, extending previous investigations on smaller aromatic alcohols. While the conformational flexibility of the ethyl group anticipates a variety of torsional isomers, the intra- and inter-molecular interactions restrict molecular conformations and favour particularly stable isomers. The conformational landscape of the clusters is very shallow and multiple competing isomers were rotationally and/or vibrationally detected, including three dimer species, two trimers and two tetramers. Cluster growth is associated with a tendency to form cyclic hydrogen bond structures.

Hydrogen Bonding in the Dimer and Monohydrate of 2-Adamantanol: A Test Case for Dispersion-Corrected Density Functional Methods.

Weakly-bound intermolecular clusters constitute reductionist physical models for non-covalent interactions. Here we report the observation of the monomer, the dimer and the monohydrate of 2-adamantanol, a secondary alcohol with a bulky ten-carbon aliphatic skeleton. The molecular species were generated in a supersonic jet expansion and characterized using broadband chirped-pulse microwave spectroscopy in the 2-8 GHz frequency region. Two different gauche-gauche O-H···O hydrogen-bonded isomers were observed for the dimer of 2-adamantanol, while a single isomer was observed for the monomer and the monohydrate. The experimental rotational parameters were compared with molecular orbital calculations using density functional theory (B3LYP-D3(BJ), B2PLYP-D3(BJ), CAM-B3LYP-D3(BJ), ωB97XD), additionally providing energetic and electron density characterization. The shallow potential energy surface makes the dimer an interesting case study to benchmark dispersion-corrected computational methods and conformational search procedures.

Molecular Rotation Spectrum of Tetracyclic Quinolizidines: Observation of trans-Matrine and the Elusive cis-Matrine.

We characterized the bis-quinolizidine tetracyclic alkaloid (5S, 6S, 7R, 11R)-matrine in a supersonic jet expansion, using chirped-pulsed broadband microwave spectroscopy. Previous crystal diffraction analyses suggested 16 diastereoisomers associated with matrine's four carbon stereocenters but were inconclusive whether the lactamic nitrogen atom would additionally produce separated trans-/cis- diastereoisomers or if both species may interconvert through low potential barriers. Our experiment simultaneously detected trans- and cis-matrine through their rotational spectrum, confirming the possibility of conformational rearrangement in matrine alkaloids. The two matrine conformers mainly differ in the envelope or half-chair lactamic ring, as evidenced by the experimental rotational and nuclear quadrupole coupling parameters. Molecular orbital calculations with ab initio (MP2) and density functional methods (B3LYP-D3(BJ) and MN15) were tested against the experiment, additionally offering an estimation of the cis-/trans- barrier of 24.9-26.9 kJ mol-1. The experiment illustrates the structural potential of chirped-pulsed broadband microwave spectroscopy for high-resolution rotational studies of biomolecules in the range of 20-40 atoms.

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.

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 N2O⋯CO. In order to assure a large initial dataset of rotational parameters, we first measured the microwave spectra of the 15N2O⋯12CO and 15N2O⋯13CO 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.

Sulfur hydrogen bonding and internal dynamics in the monohydrates of thenyl mercaptan and thenyl alcohol.

The monohydrates of thenyl alcohol and thenyl mercaptan have been probed in a supersonic jet expansion using chirped-pulse and Fabry-Perot Fourier-transform microwave spectroscopy. The rotational spectra revealed a single isomer for each of the dimers. The thenyl alcohol hydrate is stabilized by an O-HOw hydrogen bond between the alcohol and water, with water acting as a proton acceptor and additionally engaging in an Ow-Hπ interaction with the thenyl ring. Conversely, water behaves as a proton donor in the thenyl mercaptan hydrate, linking to the thiol group though an Ow-HS hydrogen bond and secondary Ow-Hπ interactions with the ring. In both dimers water retains internal mobility, as tunneling doublings in the spectrum confirm an internal rotation motion of water inside the cluster. The experimental results have been complemented with density-functional-theory molecular orbital calculations, binding energy decomposition and a topological analysis of the electronic density, providing a comparative description of the effects of hydrogen bonding of water to the alcohol and thiol groups in the dimers, relevant to understand hydrogen bonding to sulfur centers.

The Six Isomers of the Cyclohexanol Dimer: A Delicate Test for Dispersion Models.

The cyclohexanol homodimer acts as a delicate test model of the role of dispersion forces in intermolecular association. Whereas phenol produces a single dimer, the suppression of π interactions and the larger conformational flexibility in cyclohexanol results in multiple isomerism, with six competing dimers of the free molecule being observed in a supersonic jet expansion. Rotational spectroscopy reveals accurate structural data, specifically the formation of homo- and heterochiral diastereoisomers and the presence of both equatorial and axial forms in the dimers. Four dispersion-corrected density-functional molecular orbital calculations were tested against the experiment, with B3LYP-D3(BJ) offering good structural reproducibility with an Alrich's triple-ζ basis set. However, the prediction of the dimer energetics is largely model-dependent, thus offering a testbed for the validation of dispersion-corrected computational models.

The Hydrogen Bond and Beyond: Perspectives for Rotational Investigations of Non-Covalent Interactions.

In the last decade, experiment and theory have expanded our vision of non-covalent interactions (NCIs), shifting the focus from the conventional hydrogen bond to new bridging interactions involving a variety of weak donor/acceptor partners. Whereas most experimental data originate from condensed phases, the introduction of broadband (chirped-pulse) microwave fast-passage techniques has revolutionized the field of rotational spectroscopy, offering unexplored avenues for high-resolution studies in the gas phase. We present an outlook of hot topics for rotational investigations on isolated intermolecular clusters generated in supersonic jet expansions. Rotational spectra offer very detailed structural data, easily discriminating the isomeric or isotopic composition and effectively cancelling any solvent, crystal, or matrix bias. The direct comparison with quantum mechanical predictions provides insight into the origin of the inter- and intramolecular interactions with much greater precision than any other spectroscopic technique, simultaneously serving as test-bed for fine-tuning of theoretical methods. We present recent examples of rotational investigations around three topics: oligomer formation, chiral recognition, and identification of halogen, chalcogen, pnicogen, or tetrel bonds. The selected examples illustrate the benefits of rotational spectroscopy for the structural and energetic assessment of inter-/intramolecular interactions, which may help to move from fundamental research to applications in supramolecular chemistry and crystal engineering.

The S-S Bridge: A Mixed Experimental-Computational Estimation of the Equilibrium Structure of Diphenyl Disulfide.

The disulfide bridge (-S-S-) is an important structural motif in organic and protein chemistry, but only a few accurate equilibrium structures are documented. We report the results of supersonic-jet microwave spectroscopy experiments on the rotational spectra of diphenyl disulfide, C6 H5 -S-S-C6 H5 (including all 13 C and 34 S monosubstituted isotopologues), and the determination of the equilibrium structure by the mixed estimation (ME) method. A single conformation of C2 symmetry was observed in the gas phase. This disulfide is a challenging target since its structure is determined by 34 independent parameters. Additionally, ab initio calculations revealed the presence of three low-frequency vibrations (<50 cm-1 ) associated to phenyl torsions which would prevent the calculation of an accurate force field. For this reason, instead of the semiexperimental method, we used the mass-dependent (rm ) method to fit the structural parameters concurrently to moments of inertia and predicate parameters, affected with appropriate uncertainties. The predicates were obtained by high-level quantum-chemical computations. A careful analysis of the results of different fits and a comparison with the ab initio optimizations confirms the validity of the used methods, providing detailed structural information on the title compound and the disulfide bridge.

Rotational spectroscopy of organophosphorous chemical agents: cresyl and phenyl saligenin phosphates.

Cresyl and phenyl saligenin phosphates have been probed in a jet expansion by broadband chirp-excitation microwave spectroscopy, revealing the most stable conformations and their structural properties. The rotational parameters offer a high-resolution univocal route for characterization of organophosphorous agents and a testbed for computational models.

Internal dynamics of cyclohexanol and the cyclohexanol-water adduct.

Two conformers of cyclohexanol and the cyclohexanol-water adduct have been characterized in a jet expansion using rotational spectroscopy. In the gas phase, cyclohexanol adopts an equatorial position for the hydroxyl group, with the two conformers differing in the orientation of the hydroxylic hydrogen, either gauche or trans with respect to the aliphatic hydrogen at C(1). Axial cyclohexanol was not detected in the jet. The transitions of the gauche conformer are split into two component lines due to the tunneling effect of the O-H internal rotation, which connects two equivalent gauche minima. The tunneling splitting in the vibrational ground state has been determined to be ΔE0+0- = 52(2) GHz. From this splitting, the inversion barriers connecting the two equivalent gauche conformers have been determined using a flexible model to be B2 = 377 cm-1. A single isomer is detected for the cyclohexanol-water dimer, in which the water molecule acts as a proton donor to the equatorial gauche ring. The presence of torsional tunneling in the adduct suggests a concerted large-amplitude-motion in which the internal rotation in the ring is accompanied by a torsion of the water molecule, to produce an equivalent enantiomer. The torsional tunneling in the adduct is reduced to ΔE0+0- = 32.7(4) GHz and the potential barrier in the complex increases to B2 = 494 cm-1.

How flexible is the disulfide linker? A combined rotational-computational investigation of diallyl disulfide.

The symmetrically substituted diallyl disulfide adopts a non-symmetric conformation in the gas-phase, as observed with supersonic-jet rotational spectroscopy. The determination of the equilibrium structure with a predicate mixed regression illustrates both the benefits of the mass-dependent method for moderately large molecules and the structural peculiarities of the disulfide bridge.

Rotational Spectrum, Tunneling Motions, and Intramolecular Potential Barriers in Benzyl Mercaptan.

The rotational spectrum of benzyl mercaptan (parent and four isotopologues) has been assigned in a supersonic jet expansion using chirped-pulse Fourier transform microwave spectroscopy. The spectrum is characterized by torsional tunneling doublings, strongly perturbed by Coriolis interactions. The experimental rotational constants reveal that the sulfur atom is located above the ring plane in a gauche conformation. The torsion dihedral θ0 = φ (SCα-C1C2) is approximately 74°, according to a flexible molecular model calculation reproducing the energy separation (ΔE01 ∼ 2180.4 MHz) between the first two torsional substates. The global minimum configuration is 4-fold degenerate, corresponding to potential minima with θ0 ≈ ±74° and ±(180-74)°. The four equivalent minima are separated by potential barriers at θ = ±90°, 0°, or 180°. The tunneling splittings are caused by the potential barrier at θ = ± 90°, while the barriers at torsions of 0° or 180° are too large to generate detectable splittings. The tunnelling barrier has been determined as 248 cm-1, similar to the value obtained with high-level MP2 ab initio calculations (259 cm-1), but smaller than in benzyl alcohol (280 cm-1).

Sulfur Hydrogen Bonding in Isolated Monohydrates: Furfuryl Mercaptan versus Furfuryl Alcohol.

The hydrogen bonds involving sulfur in the furfuryl mercaptan monohydrate are compared with the interactions originating from the hydroxyl group in furfuryl alcohol. The dimers with water were created in a supersonic jet expansion and characterized using microwave spectroscopy and supporting molecular orbital calculations. In furfuryl alcohol-water, a single isomer is observed, in which the water molecule forms an insertion complex with two simultaneous hydrogen bonds to the alcohol (O-H⋅⋅⋅Ow ) and the ring oxygen (Ow -H⋅⋅⋅Or ). When the alcohol is replaced by a thiol group in furfuryl mercaptan-water, two isomers are observed, with the thiol group preferentially behaving as proton donor to water. The first isomer is topologically equivalent to the alcohol analog but the stronger hydrogen bond is now established by water and the ring oxygen, assisted by a thiol S-H⋅⋅⋅Ow hydrogen bond. In the second isomer the sulfur group accepts a proton from water, forming a Ow -H⋅⋅⋅S hydrogen bond. Binding energies for the mercaptan-water dimer are predicted around 12 kJ mol-1 weaker than in the alcohol hydrate (B3LYP-D3(BJ)). The non-covalent interactions in the furfuryl dimers are dominantly electrostatic according to a SAPT(0) energy decomposition, but with increasing dispersion components in the mercaptan dimers, which are larger for the isomer with the weaker Ow -H⋅⋅⋅S interaction.