Molecular Recognition, Transient Chirality and Sulfur Hydrogen Bonding in the Benzyl Mercaptan Dimer

: The homodimers of transiently chiral molecules offer physical insight into the process of molecular recognition, the preference for homo or heterochiral aggregation and the nature of the non-covalent interactions stabilizing the adducts. We report the observation of the benzyl mercaptan dimer in the isolation conditions of a supersonic jet expansion, using broadband (chirped-pulse) microwave spectroscopy. A single homochiral isomer was observed for the dimer, stabilized by a cooperative sequence of S-H ··· S and S-H ··· π hydrogen bonds. The structural data, stabilization energies and energy decomposition describe these non-covalent interactions as weak and dispersion-controlled. A comparison is also provided with the benzyl alcohol dimer.


Introduction
Molecular recognition is a subtle chemical process involving a combination of intermolecular interactions and intramolecular factors, largely dependent on the composition of the host and guest molecules. The consequences of molecular recognition, extending to large-scale chemical phenomena like supramolecular [1] and technological applications [2], justify the molecular investigation of the electronic and stereochemical features of the recognition process using high-resolution spectroscopy [3] and quantum mechanical methods [4,5]. In this context, the preparation of mass-selected intermolecular adducts in the gas phase is a well-known procedure for isolation and observation of specific weak interactions between target chemical groups, which otherwise would be blurred in the condensed media. This approach, when combined with rationalizing theoretical calculations has led, in the last decade, to an explosion of new weak intermolecular interactions previously unnoticed [6,7].
Molecular recognition is especially interesting between chiral species, as it may provide insight into biochemical docking, asymmetric synthesis, and chiral analysis. Among chiral molecules, those with very low (5-10 kJ mol −1 ) torsional stereomutation barriers display transient chirality. Transient enantiomers interconvert in nanosecond time scales that would be undetectable with conventional techniques. However, transient chirality can be revealed by the formation of diastereomeric dimers in the gas phase, simultaneously freezing stereomutation and providing information on the structural and energetic factors controlling homo-or heterochiral aggregation [8,9]. The process of chirality synchronization was first observed with electronic and vibrational spectroscopy [10][11][12], but these techniques require computational support for spectral interpretation. Alternatively, rotational spectroscopy [7,13] offers an unequivocal structural identification that complements molecular orbital calculations and may contribute to validating the computational models, which in turn are critical for the determination of the weak chirodiastaltic energies [14][15][16] between homo-and heterochiral aggregates.
In this work, we selected benzyl mercaptan as a dimerization target. Benzyl mercaptan represents an extension of our previous investigation on thiophenol dimerization [36] and will establish if the two rings maintain the π-stacking thiophenol arrangement or tilted geometries similar to the benzyl alcohol dimer [42], recently revisited [43]. Additionally, it will discern whether the homochiral aggregation of benzyl alcohol is respected in the mercaptan and the structural and physical differences in the S-H···S hydrogen bond compared to the canonical O-H···O hydrogen bond.

Materials and Methods
The sample of benzyl mercaptan (>96%, bp 195 • C) was obtained commercially and required no further purification. The sample was vaporized inside the heating reservoir (45-55 • C) of a pulsed solenoid injector and expanded near adiabatically through a 0.8 mm nozzle to form a supersonic jet within an expansion chamber. Neon at stagnation pressures of 0.2 MPa was used as the carrier gas, with typical molecular pulses of 800-900 µs. The expanding jet was probed in the 2-8 GHz cm-wave region with a direct-digital chirpedpulse Fourier transform microwave (CP-FTMW) spectrometer, following Pate's design [44]. The operation sequence is based on short (4 µs, 20 W) chirped pulses, which are broadcasted perpendicularly to the jet. The MW radiation produces a fast-passage broadband transient excitation [45,46], simultaneously covering the full spectral bandwidth. Following the excitation, the molecular ensemble emits a free-induction decay, which is detected in the time-domain (ca. 40 µs) and acquired using a (25 GSamples/s) digital oscilloscope. A Fourier transformation with a Kaiser-Bessel window results in FWHM linewidths of ca. 100 kHz. In this experiment, ca. 1 M averages were acquired at a repetition rate of 5 Hz. The uncertainty of the frequency measurements was estimated below 20 kHz.
Several computational calculations complemented the experimental study. Following an initial conformational screening with molecular mechanics (MMFFs [47]), all further calculations used density-functional theory (DFT) molecular orbital calculations. Two density-functional methods were selected here, including B3LYP [48] and the double hybrid B2PLYP [49] method, in combination with Ahlrichs' polarized triple-zeta basis def2-TZVP [50]. Both methods were supplemented with D3 [51] dispersion corrections and Becke-Johnson damping [52]. Frequency calculations were performed at the same level of theory, using the harmonic approximation. The calculation of complexation energies considered the basis set superposition errors (BSSE) with the counterpoise approximation [53]. All DFT calculations were conducted with Gaussian 16 [54]. The physical contributions to the binding potential of the water clusters were estimated by energy decomposition analysis using second-order symmetry adapted perturbation theory [53,55] (SAPT), implemented in PSI4 [56]. Finally, the presence of non-covalent interactions was analyzed with the NCIPlot methodology, based on a reduced gradient of the electronic density [57].

Benzyl Mercaptan Monomer
Benzyl mercaptan displays a bidimensional potential energy surface, qualitatively similar to benzyl alcohol [58]. The two torsional degrees of freedom are associated with the elevation and orientation of the terminal thiol group, given by dihedrals τ 1 (SC α -C ipso C ortho ) and τ 2 (HS-C α C ipso ). Our previous rotational investigation [59] of benzyl mercaptan confirmed that the monomer presents a single isomer in the gas phase, as in the alcohol. In this conformation the sulfur atom is synclinal to the ring plane (τ 1± 74 • vs. ±55 • in benzyl alcohol) and the thiol hydrogen is synchronously oriented towards the π ring (τ 2~± 74 • vs. ±53 • in benzyl alcohol), denoted gauche-gauche or GG (Tables S1 and S2, Supporting Information, SI). The global minimum is four-fold degenerate, as the thiol group may tunnel between symmetry equivalent conformations either by reflection on the ring plane or in a perpendicular plane bisecting the phenyl ring through carbons C ipso and C para . Inversion through the perpendicular plane creates a detectable barrier (B 2 = 248 cm −1 = 2.97 kJ mol −1 ), which produces characteristic tunneling doublings in the rotational spectrum of Figure S1 (SI) [59], strongly perturbed by Coriolis interactions [60]. The inversion barrier in benzyl mercaptan is ca. 11% lower than in the alcohol (B 2 = 280 cm −1 = 3.35 kJ mol −1 ). The molecule thus changes chirality by inverting between the two gauche-gauche enantiomers GG and G G of Figure 1 (primes denote negative dihedrals). implemented in PSI4 [56]. Finally, the presence of non-covalent interactions was analyzed with the NCIPlot methodology, based on a reduced gradient of the electronic density [57].

Benzyl Mercaptan Monomer
Benzyl mercaptan displays a bidimensional potential energy surface, qualitatively similar to benzyl alcohol [58]. The two torsional degrees of freedom are associated with the elevation and orientation of the terminal thiol group, given by dihedrals τ1(SCα-Cipso-Cortho) and τ2(HS-CαCipso). Our previous rotational investigation [59] of benzyl mercaptan confirmed that the monomer presents a single isomer in the gas phase, as in the alcohol. In this conformation the sulfur atom is synclinal to the ring plane (τ1~ ±74° vs. ±55° in benzyl alcohol) and the thiol hydrogen is synchronously oriented towards the π ring (τ2~ ±74° vs. ±53° in benzyl alcohol), denoted gauche-gauche or GG (Tables S1 and S2, Supporting Information, SI). The global minimum is four-fold degenerate, as the thiol group may tunnel between symmetry equivalent conformations either by reflection on the ring plane or in a perpendicular plane bisecting the phenyl ring through carbons Cipso and Cpara. Inversion through the perpendicular plane creates a detectable barrier (B2 = 248 cm −1 = 2.97 kJ mol −1 ), which produces characteristic tunneling doublings in the rotational spectrum of Figure S1 (SI) [59], strongly perturbed by Coriolis interactions [60]. The inversion barrier in benzyl mercaptan is ca. 11% lower than in the alcohol (B2 = 280 cm −1 = 3.35 kJ mol −1 ). The molecule thus changes chirality by inverting between the two gauche-gauche enantiomers GG and G'G' of Figure 1 (primes denote negative dihedrals).

Benzyl Mercaptan Homodimer
For benzyl mercaptan the dimerization from a single gauche-gauche conformation may in principle produce only homochiral (GG-GG = G'G'-G'G') or heterochiral (GG-G'G' = G'G'-GG) diastereoisomers (proton donor denoted first). The most stable isomers are expected to balance several intermolecular interactions, mainly involving the S-H and C-H groups as proton donors and the sulfur lone pairs and π electrons as acceptors. However, the two diastereotopic acceptor lone pairs (Lp+/−) and the multiple possibilities for interaction between the polar thiol groups and the two ring molecules are expected to produce multiple isomers. The results of a DFT (B3LYP-D3(BJ)) conformational search are summarized in Table S3 and Figure S2 (SI), suggesting 11 isomers within complexation energies below 10 kJ mol −1 . The five most stable isomers are compared to the benzyl alcohol dimer [42] in Figure 2. The two most stable dimers 1 (= GG-GG-Lp−) and 2 (= GG-GG-Lp+) are homochiral and controlled by the two thiol groups (3D Figure S3, SI), which engage in two successive S-H···S and S-H···π interactions. The two rings adopt tilted (C1) unsymmetric orientations to permit further interaction via C(sp 2 )-H···π weak hydrogen bonds, with the alternative use of the two sulfur lone-pairs at the acceptor molecule producing two different ring orientations Lp+ and Lp− ( Figure S4, SI), well separated in energy (ΔEZPE = 1.9 kJ mol −1 , ΔEc = 3.1 kJ mol −1 ). The third isomer (=G'G'-GG-Lp−) is the heterochiral equivalent to isomer 1 ( Figure S5, SI), predicted at ΔEZPE = 3.3 kJ mol −1 (ΔEc = 4.6

Benzyl Mercaptan Homodimer
For benzyl mercaptan the dimerization from a single gauche-gauche conformation may in principle produce only homochiral (GG-GG = G G -G G ) or heterochiral (GG-G G = G G -GG) diastereoisomers (proton donor denoted first). The most stable isomers are expected to balance several intermolecular interactions, mainly involving the S-H and C-H groups as proton donors and the sulfur lone pairs and π electrons as acceptors. However, the two diastereotopic acceptor lone pairs (Lp+/−) and the multiple possibilities for interaction between the polar thiol groups and the two ring molecules are expected to produce multiple isomers. The results of a DFT (B3LYP-D3(BJ)) conformational search are summarized in Table S3 and Figure S2 (SI), suggesting 11 isomers within complexation energies below 10 kJ mol −1 . The five most stable isomers are compared to the benzyl alcohol dimer [42] in Figure 2. The two most stable dimers 1 (= GG-GG-Lp−) and 2 (= GG-GG-Lp+) are homochiral and controlled by the two thiol groups (3D Figure S3, SI), which engage in two successive S-H···S and S-H···π interactions. The two rings adopt tilted (C 1 ) unsymmetric orientations to permit further interaction via C(sp 2 )-H···π weak hydrogen bonds, with the alternative use of the two sulfur lone-pairs at the acceptor molecule producing two different ring orientations Lp+ and Lp− ( Figure S4, SI), well separated in energy (∆E ZPE = 1.9 kJ mol −1 , ∆E c = 3.1 kJ mol −1 ). The third isomer (=G G -GG-Lp−) is the heterochiral equivalent to isomer 1 ( Figure S5, SI), predicted at ∆E ZPE = 3.3 kJ mol −1 (∆E c = 4.6 kJ mol −1 ) The fourth (=GG-GG-ππ) and fifth (=GG-GG-π) isomers reveal different homochiral S-H···π interactions, either alone or combined with a C(sp 3 )-H···π link, which is destabilized ca. 4-5 kJ mol −1 with respect to the thiol S-H···S hydrogen bonding. Isomer GG-GG-ππ uses two identical S-H···π interactions to produce the lowest-lying C 2 -symmetric dimer, characterized by two tilted rings which avoid π-stacking ( Figure S6, SI). Only for the higher energy isomers (>10 kJ mol −1 ) near-parallel π···π or C-H···π interactions are predicted, as in the (C 2 ) homochiral isomer 9 (=GG-GG) and the (C i ) heterochiral isomer 11 (=GG-G G ). The preference for C 2 -symmetric heterochiral S-H···π dimers observed in the crystallographic structure [61] must thus be attributed to matrix effects and does not represent the isolated molecule nor the dimer. The three most stable isomers of the benzyl mercaptan dimer were reoptimized at B2PLYP-D3(BJ) level in Table 1, offering a prediction of the rotational parameters. kJ mol −1 ) The fourth (=GG-GG-ππ) and fifth (=GG-GG-π) isomers reveal different homochiral S-H···π interactions, either alone or combined with a C(sp 3 )-H···π link, which is destabilized ca. 4-5 kJ mol −1 with respect to the thiol S-H···S hydrogen bonding. Isomer GG-GG-ππ uses two identical S-H···π interactions to produce the lowest-lying C2-symmetric dimer, characterized by two tilted rings which avoid π-stacking ( Figure S6, SI). Only for the higher energy isomers (>10 kJ mol −1 ) near-parallel π···π or C-H···π interactions are predicted, as in the (C2) homochiral isomer 9 (=GG-GG) and the (Ci) heterochiral isomer 11 (=GG-G'G'). The preference for C2-symmetric heterochiral S-H···π dimers observed in the crystallographic structure [61] must thus be attributed to matrix effects and does not represent the isolated molecule nor the dimer. The three most stable isomers of the benzyl mercaptan dimer were reoptimized at B2PLYP-D3(BJ) level in Table 1, offering a prediction of the rotational parameters.  Figure S2 for a larger energy window). The right column represents the most stable isomers of the benzyl alcohol dimer.  The jet-cooled microwave spectrum of Figure S1 provided experimental evidence on the nature of the dimerization adduct. Figure 3 shows a spectral section with typical rotational transitions. Noticeably, a single isomer was observed for the dimer of benzyl mercaptan. The spectrum was quite dense and more than 300 rotational transitions were measured experimentally, with all three (µ a , µ b , µ c ) selection rules active. The wide range of angular momentum quantum numbers (J = 3-15 and K −1 < 12) and the diversity of R-branch (J + 1←J) and some Q-branch (J←J) transitions assured a good determination of the rotational parameters. There was no indication of tunneling effects in the spectra, which were reproduced satisfactorily with a Watson's (S-reduced) semirigid rotor Hamiltonian [62] and quartic centrifugal distortion terms. The experimental rotational parameters are compared with B2PLYP-D3(BJ) theoretical predictions in Table 1, with the full listing of observed transitions in Table S4 (SI). Since several low-lying isomers were predicted at small relative energies (<3-5 kJ mol −1 ), the observation of a single dimer species cannot be uniquely attributed to a pre-expansion thermal depopulation and probably suggests a conformational relaxation through low interconversion barriers, previously observed in jet experiments [63,64] and kinetically affecting the jet populations [65].
The identification of the benzyl mercaptan dimer was unequivocal. We confirmed the detection of the predicted global minimum, identified as the homochiral GG-GG-Lp−. The predictions for the equilibrium rotational constants slightly exceed the experimental ground-state values by 6-8 MHz, but the relative differences are below 1.6-2.8%, ensuring a positive identification. The agreement of the (harmonic) centrifugal distortion constants with the experiment values is consistent with the conformational assignment. compared with B2PLYP-D3(BJ) theoretical predictions in Table 1, with the full listing of observed transitions in Table S4 (SI). Since several low-lying isomers were predicted at small relative energies (<3-5 kJ mol −1 ), the observation of a single dimer species cannot be uniquely attributed to a pre-expansion thermal depopulation and probably suggests a conformational relaxation through low interconversion barriers, previously observed in jet experiments [63,64] and kinetically affecting the jet populations [65]. The identification of the benzyl mercaptan dimer was unequivocal. We confirmed the detection of the predicted global minimum, identified as the homochiral GG-GG-Lp−. The predictions for the equilibrium rotational constants slightly exceed the experimental ground-state values by 6-8 MHz, but the relative differences are below 1.6-2.8%, ensuring a positive identification. The agreement of the (harmonic) centrifugal distortion constants with the experiment values is consistent with the conformational assignment.  Table 1. The unmarked transitions belong to different species.

Non-Covalent Interactions
The intermolecular interactions associated to the benzyl mercaptan dimerization were explored using structural, energetic, and electronic density information. Unlike the parallel-displaced hydrogen-bonded dimer of thiophenol [36], the aggregation of benzyl mercaptan follows a pattern of non-stacked cooperative hydrogen bonding. A primary S-H···S hydrogen bond is formed in which the thiol acts as a proton donor to a second sulfur atom, followed by a secondary S-H···π hydrogen bond to the opposed π ring. Weak C-H···π or C-H···S hydrogen bonds are minor necessary contributors to the most stable isomers. The benzyl mercaptan dimers exhibit characteristically long S-H···S hydrogen bonds (B2PLYP:  Table 1 show distances to the ring centroid of r(S-H···centroid) = 2.52 Å (B2PLYP).
Binding energies are collected in Tables 1 and 2, comparing B2PLYP (with BSSE corrections) and a second-order 2+(3) symmetry-adapted perturbation theory (SAPT) energy decomposition (at the optimized geometries of B3LYP-D3(BJ)/def2-TZVP). Not unexpectedly, the predicted binding energy of the benzyl mercaptan dimer (SAPT: −35.3 kJ mol −1 ) is larger than in thiophenol (−25.9 to −27 kJ mol −1 ) [36] but remains ca. 84% smaller than in the benzyl alcohol dimer (−42.1 kJ mol −1 ) [43], reflecting the weaker thiol H-bond. Additional insight into the nature of the thiol hydrogen bond is obtained from the attractive contributions to the energy decomposition [67,68]. The largest attractive contributor in the benzyl mercaptan dimer is dispersive (54%) but in smaller proportion than in the thiophenol dimer (60%) and with larger electrostatic participation (34% vs. 29-31% in thiophenol). This result reverts the calculation for the benzyl alcohol dimer, where the electrostatic component (45%) is larger but nearly matched by dispersion (41%). A comparison with the dispersive pyridine-methane dimer and the water and hydrogen sulfide dimers is also given in Table 2. Similar trends of decreasing strength and dispersion have been observed in other thiol clusters and monohydrates [30,31]. To our knowledge, there are no comparable experimental binding energies involving thiol dimers except for hydrogen sulfide clusters [26,27].
Non-covalent interactions can also be analyzed using the topological properties of the electron density (ρ(r)) with the NCIPlot method [57]. Figure 4 represents a reduced electronic density gradient s = versus the signed electronic density (sign (λ 2 ) ρ) using the second eigenvalue (λ 2 ) of the electron density Hessian, comparing the global minimum of the benzyl mercaptan and benzyl alcohol dimers (B3LYP-D3(BJ)/def2-TZVP geometries). This representation qualitatively identifies the stronger O-H···O hydrogen bond interactions in the alcohol (at more negative abscissas), simultaneously providing a 3D visualization of the most relevant interaction surfaces. The 3D plot clearly identifies the S-H···S hydrogen bond and broad interaction regions between the two rings mainly associated with the S-H···π contact. The reduced electronic density calculations thus confirm the cooperative scheme of intermolecular interactions present in the benzyl mercaptan dimer, complementing the structural and energetic descriptions in Tables 1 and 2. PEER REVIEW Figure 4. NCI plot (left panel) and reduced electronic density gradient (right panel) for the dimer of benzyl mercaptan. In the NCI plot, blue shades indicate attractive interactions (a to the S-H···S hydrogen bond), green colors indicate weak attractive interactions (like S-H red represents repulsive interactions (ring critical points). The reduced electronic density attractive (negative minima) and repulsive (positive minima) interactions in the dimers mercaptan (green) and benzyl alcohol (red), with a stronger O-H···O interaction.

Discussion
The benzyl mercaptan dimer is a model cluster with a primary amphot group and a combination of aliphatic proton donors and a π ring acceptor. The m pivot between the thiol and the aromatic ring permits some conformational flexib The reduced electronic density shows the attractive (negative minima) and repulsive (positive minima) interactions in the dimers of benzyl mercaptan (green) and benzyl alcohol (red), with a stronger O-H···O interaction. Table 2. Binding energy decomposition for the dimer of benzyl mercaptan, and comparison with several thiol and alcohol dimers and the van der Waals pyridine-methane dimer. The interaction energy is decomposed into electrostatic (∆E elec ), inductive (multipole interactions/charge transfer, ∆E ind ), exchange repulsion (∆E exch ) and dispersion (∆E disp ) energy terms, using SAPT(2)+3/aug-cc-pVDZ//B3LYP-D3(BJ)/def2-TZVP (all values in kJ mol −1 ).

Discussion
The benzyl mercaptan dimer is a model cluster with a primary amphoteric thiol group and a combination of aliphatic proton donors and a π ring acceptor. The methylene pivot between the thiol and the aromatic ring permits some conformational flexibility and a variety of intermolecular forces, generally balancing two cooperative interactions. The rotational spectrum confirmed a single isomer in the gas phase, providing data for comparison between experiment and theory. The observed global minimum (GG-GG-Lp−) is homochiral and characterized by a leading thiol-thiol (S-H···S) hydrogen bond, in cooperation with a secondary thiol-aromatic (S-H···π) hydrogen bond. Characteristically, this interaction pattern and stereochemistry is coincident with the most stable (hom-O g π-I) benzyl alcohol dimer, stabilized by stronger O-H···O and O-H···π alcohol hydrogen bonds [42,43]. In consequence, the expected reduction in hydrogen bond strength of the thiol group does not affect the conformational equilibrium observed in the alcohol, and the π ring acceptor maintains its secondary role in both dimers. The smaller, more rigid phenol and thiophenol dimers, also observed in the gas phase [36,69], similarly depend on the S-H···S or O-H···O hydrogen bonds but adopt different hinged or stacked orientations which reflect a different balance with the ring interactions. Since the number of experimental observations is limited additional thiol-alcohol comparisons are presently possible only for a few other aromatic dimers like 2-phenylethanethiol [29,70], also showing similar thiol-thiol interactions in the global minimum.
The second isomer of the mercaptan dimer (GG-GG-Lp+), separated 1.7 kJ mol −1 (B2PLYP)-1.8 kJ mol −1 (B3LYP), differs in the use of a different sulfur lone-pair in the acceptor molecule but maintains its homochiral character. The equivalent isomer in the alcohol (hom-O t π-V: 3.2 kJ mol −1 ) is much higher in energy. Consequently, the benzyl mercaptan dimer shows a reinforced preference for homochirality compared to the alcohol, where the heterochiral isomer (het-O g π-II) is nearly isoenergetic (0.2 kJ mol −1 ) with the (hom-O g π-I) global minimum. Conversely, in the benzyl mercaptan dimer the heterochiral isomer (G G -GG-Lp− ≡ het-O g π-II) is third in energy and more separated ((B2PLYP: 3.4 kJ mol −1 ; B3LYP: 3.3 kJ mol −1 ) from the homochiral partner.
The fourth and fifth isomers of the benzyl mercaptan dimer do not show thiol-thiol hydrogen bonds, which are replaced by separate thiol interactions with different proton acceptors. These isomers introduce the discussion on the strength of the thiol or alcohol groups as proton acceptors in presence of competing groups, which has been mostly studied for alcohol dimers with two acceptor groups. In those cases, the observed alcohol-alcohol hydrogen bonding of (benzyl alcohol) 2 is confirmed with related aromatic acceptors, as in (indanol) 2 [71] or (propargyl alcohol) 2 [72]. However, ethers [21,73], carbonyl groups [74], or carboxylic acids [75] show preference as primary proton acceptors, more intense in amino [76] or cyano [77] groups. Additionally, when the two thiol/alcohol groups bind separately to secondary acceptors both homo and heterochiral symmetric geometries become possible. In the benzyl mercaptan dimer, the fourth isomer (GG-GG-ππ, 4.4 kJ mol −1 ) shows a homochiral C 2 antiparallel geometry stabilized by two S-H···π hydrogen bonds, higher in energy than the equivalent C 2 benzyl alcohol dimer (homππ-IV, 2.0 kJ mol −1 ) and reminiscent of the antiparallel aniline dimer [78]. However, the heterochiral C i benzyl alcohol dimer (hetππ-III, 1.6 kJ mol −1 ) does not appear among the preferred geometries of the benzyl mercaptan dimer, where the fifth unsymmetric isomer (GG-GG-π, 4.6 kJ mol −1 ) instead shows a combination of S-H···π and C-H···π interactions.
Differences in the physical forces operating in the benzyl mercaptan and alcohol dimers are notorious in the topological analysis of the electronic density and the SAPT energy decomposition of Table 2 and Figure 3. The weaker, dispersive character of the thiol dimer, thus, contributes to a better description of non-covalent forces involving low-electronegativity atoms.
Future advances in the understanding of intermolecular forces are expected to extend to different chemical groups, incorporating a synergic combination of empirical data and computational models. In this work, rotational spectroscopy has an increasingly important role thanks to the development of chirped-pulsed broadband techniques.

Conclusions
The combination of rotational data and molecular orbital calculations constitutes a very effective way of interrogating weak non-covalent interactions, simultaneously extending the knowledge base on intermolecular forces and the benchmark capacity of the experiments. The availability of empirical data is expected to contribute to the development of improved density-functional methods, which will be essential for the molecular modelling of large chemical systems.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/sym13112022/s1, Figure S1: Rotational spectrum of benzyl mercaptan in the 3-8 GHz frequency region. The positive trace shows the experimental spectrum; the negative trace is the simulation of fitted rotational constants of the monomer and the dimer (isomer 1, GG-GG-Lp−), Figure S2: The most stable isomers of the benzyl mercaptan dimer, classified by their intermolecular interactions and ordered by complexation energy (∆E c , kJ/mol). Isomers in the dashed rectangle exhibit a sulfur S-H···S hydrogen bond. The only symmetric isomers are 5 (C 2 ), 9 (C 2 ) and 11 (C i ), Figure S3: A rotatable 3D figure of the two most stable isomers of the benzyl mercaptan dimer, Figure S4: A comparison of the three asymmetric (C 1 ) most stable isomers of the benzyl mercaptan dimer. Isomers 1 (G G -G G -Lp+ = GG-GG-Lp−, left) and 2 (G G -G G -Lp− = GG-GG-Lp+, right), are both homochiral but differ in the acceptor lone pair position of sulfur (the dihedral formed by the lone pair or τ(LpS-CαCipso) have different signs for isomer 1 and isomer 2), producing a change of relative orientation between the rings, Figure S5: A comparison of the three asymmetric (C 1 ) most stable isomers of the benzyl mercaptan dimer. Isomers 1 (G G -G G -Lp+ = GG-GG-Lp−, left) and 3 (GG-G G -Lp+ = G G -GG-Lp−, right) differ in the stereochemistry of the donor group, either homochiral in isomer 1 or heterochiral in isomer 3, Figure S6: The most stable symmetric dimer of benzyl mercaptan (isomer 5 GG-GG-ππ) has C 2 symmetry but avoids the parallel arrangement of the two phenyl rings, unlike in the symmetric benzyl alcohol dimer. The Figure shows the (a, b,  c) principal inertial axes, with a vector pointing in the C 2 (≡b) symmetry axis direction, Figure S7: A comparison of the NCI Plots for isomers 1 (G G -G G -Lp+ = GG-GG-Lp−, upper panel) and 2 (G G -G G -Lp− = GG-GG-Lp+, lower panel) of the benzyl mercaptan dimer, Table S1: Rotational parameters of benzyl mercaptan, Table S2: Atomic coordinates of benzyl mercaptan gauche GG, Table S3: Conformational search for the benzyl mercaptan dimer (B3LYP-D3(BJ)) and comparison with the experimental rotational parameters,