Edelsys Codorniu-Hernandez, Adrian Daniel Boese, Peter G. Kusalik

• Despite the critical importance of the hydroxyl radical in major scientific fields, there are still open questions on the behavior of this species in the aqueous phase. In particular, there has been much debate on the existence of a hemibonded interaction between the hydroxyl radical and water molecules. While some reports indicate that the hemibonded radical might explain some experimental data, others have claimed that this interaction is simply a density functional theory (DFT) artifact. Here, we provide results from high level (basis set limit of coupled-cluster levels up to single, double, triple excitations (CCSD(T)) and beyond) ab initio calculations of different OH center dot(H2O)(n) clusters in the gas phase to accurately explore the existence of the hemibonded interaction and its energy difference with respect to other well-defined hydrogen bond interactions. Additional comparisons with second order perturbation theory (MP2) and DFT are also presented. Constrained molecular dynamics was applied to determine the free energyDespite the critical importance of the hydroxyl radical in major scientific fields, there are still open questions on the behavior of this species in the aqueous phase. In particular, there has been much debate on the existence of a hemibonded interaction between the hydroxyl radical and water molecules. While some reports indicate that the hemibonded radical might explain some experimental data, others have claimed that this interaction is simply a density functional theory (DFT) artifact. Here, we provide results from high level (basis set limit of coupled-cluster levels up to single, double, triple excitations (CCSD(T)) and beyond) ab initio calculations of different OH center dot(H2O)(n) clusters in the gas phase to accurately explore the existence of the hemibonded interaction and its energy difference with respect to other well-defined hydrogen bond interactions. Additional comparisons with second order perturbation theory (MP2) and DFT are also presented. Constrained molecular dynamics was applied to determine the free energy for the formation/disruption and ice systems. Overall, our findings confirm that the hemibond can be an alternative structure for the hydroxyl radical in the condensed phase when the formation of hydrogen bonds is impeded. These results will aid the understanding of theoretical and experimental data and help future experimental designs for the detection of this important species.…