Simulation Studies of the Dynamics and the Connectivity Patterns of Hydrogen Bonds in Water from Ambient to Supercritical Conditions

Pressurized high-temperature water attracts attention as a promising medium for chemical synthesis, biomass processing or destruction of hazardous waste. Adjustment to the desired solvent properties requires knowledge on the behavior of populations of hydrogen-bonded molecules. In this work, the interconnection between the hydrogen bond (HB) dynamics and the structural rearrangements of HB networks have been studied by molecular dynamics simulation using the modified central force flexible potential and the HB definition controlling pair interaction energy, HB length and HB angle. Time autocorrelation functions for molecular pairs bonded continuously and intermittently and the corresponding mean lifetimes have been calculated for conditions ranging from ambient to supercritical. A significant reduction in the continuous and intermittent lifetimes has been found between (293 K, 0.1 MPa) and (373 K, 25 MPa) and attributed to the decreasing size of patches embedded in the continuous HB network. The loss of global HB connectivity at ca. (573 K, 10 MPa) and the investigated supercritical conditions do not noticeably affect the HB dynamics. Over the whole temperature range studied, the reciprocal intermittent lifetime follows the transition state theory dependence on temperature with the activation energy of 10.4 kJ/mol. Calculations of the lifetime of molecules that do not form hydrogen bonds indicate that at supercritical temperatures, the role of reactions involving an unbound H₂O molecule as a reactant can be enhanced by lowering system density.