The correspondence between the infrared vibrational spectra and the underlying hydrogen bonding network in aqueous clusters: caveats and tactics*

The structural -- spectral correspondence relates the observed infrared (IR) vibrational spectra to the underlying molecular structure. In the case of hydrogen bonded clusters the IR ``fingerprint'' region in the 3,000 -- 4,000 cm$^{-1}$ spectral range provides a direct probe of the connectivity and dynamics of the cluster's hydrogen bonding network. For medium size ($n $ 20) aqueous neutral and ionic clusters, the presence of several closely lying isomers that differ substantially in the oxygen atom network complicates both the sampling of the respective potential energy surfaces as well as the accurate determination of their energy order. Traditionally, a hierarchical procedure based on initial sampling with classical potentials and subsequent refinement with electronic structure methods has been used. We will highlight representative examples for the (H$_{2}$O)$_{20}$, (H$_{2}$O)$_{25}$ and H$_{3}$O$^{+}$(H$_{2}$O)$_{20}$ clusters [1-3] where sampling with classical potentials fails to produce the most stable minima and outline approaches and strategies that are based on a combination of enhanced sampling of configurations in conjunction with electronic structure theory to obtain realistic cluster configurations that are consistent with the measured IR spectra. \\[4pt] [1] S. S. Xantheas, \textit{Can. J. Chem. Eng.} \textbf{90}, 843 (2012).\\[0pt] [2] N. Sahu, S. R. Gadre, P. Bandyopadhyay, E. Miliordos and S. S. Xantheas, \textit{J. Chem. Phys.} \textbf{141}, 164304 (2014).\\[0pt] [3] C. C. Pradzynski, C. W. Dierking, F. Zurheide, R. M. Forck, T. Zeuch, U. Buck and S. S. Xantheas, \textit{Phys. Chem. Chem. Phys.} \textbf{16}, 26691 (2014).