Methanethiol (M) and water (W) clusters like dimers (M1W1, M2, and W2), trimers
(M1W2, M2W1, M3,
and W3), and tetramers (M1W3, M2W2, M3W1, M4,
and W4) were studied to assess the strength of sulfur-centered
hydrogen bonding using different levels of theories, viz, HF, MP2,
MP3, MP4, B3LYP, B3LYP-D3, CCSD, CCSD(T)-F12, and CCSD(T) along with
aug-cc-pVNZ (where N = D, T, and Q) basis sets. Interaction energies
were found to be in the range of −3.3 to −5.3 kcal/mol
for the dimers, −8.0 to −16.7 kcal/mol for the trimers,
and −13.5 to −29.5 kcal/mol for the tetramers at the
B3LYP-D3/CBS limit level of theory. Normal modes of vibrations computed
at the B3LYP/cc-pVDZ level of theory were seen to be in good agreement
with the experimental values. Local energy decomposition calculations
using the DLPNO-CCSD(T) level of theory indicated the domination of
electrostatic interactions’ contribution to the interaction
energy in all cluster systems. Furthermore, atoms in molecules and
natural bond orbital calculations both carried out at the B3LYP-D3/aug-cc-pVQZ
level of theory aided in visualizing the hydrogen bonds besides proving
a rationale for the strength of the hydrogen bonds and thereby the
stability of these cluster systems.