Structure and rotational dynamics of water around hydrogen peroxide

Abstract

The density functional theory-based molecular dynamics simulations were induced to study various properties of water around different solvation pockets around hydrogen peroxide(H2O2). Our study aims to understand water behavior in the vicinity of H2O2 molecule through analyzing structure, hydrogen bond statistics, orientational and non-diffusive jump profile. A water-like peroxide molecule can form four hydrogen bonds with nearby water molecules: two as a proton donor and two as a proton acceptor, which must disturb the local structure and dynamics of the surrounding water. The presence of two different peak maxima in radial distribution function (RDF) allows us to resolve the solvation environment into three solvation sites: solvation shell I, solvation shell II, and bulk. The probability distribution of the average distance of central water molecule from 4 nearest water molecules shows a shift in the peak for solvation shell I water molecules towards a ∼0.3Å higher value compared to two other solvation shells. This further affects tetrahedral order parameters as water molecules shift from the regular tetrahedron. For the first solvation layer, the water molecules lose their tetrahedrality to a significant extent as qtet probability distribution peaks at 0.5 and broaden in the higher values of qtet without any shoulder. For the second solvation layer, we find the qtet peak at 0.87 with a minor shoulder appearing at qtet = 0.55, a sign of slight loss of tetrahedrality in the unstructured water compared to undisturbed water. Loss of tetrahedrality hints at the change of orientation of H2O molecules in the vicinity of H2O2 and affects the water-water hydrogen bond (HB) formation. A noteworthy reduction in the HB donor/acceptor account and per water HB count shows interference of H2O2. The water molecules prefer H2O2 compared to neighboring water molecules for HB formation. The orientation autocorrelation function provides the rotational pace of water molecules in the different solvation layers. The reorientation time constant difference in the solvation shells I and II fluctuate by mere 0.20 ps, showcasing the trivial nature of hydrogen peroxide and water HBs. Using the HB descriptors, we also report the non-diffusive jump profile involving vicinal water molecules of H2O2. The jump angle profile of solvation shell water molecules is in line with bulk water.