Structure and polarization properties of water: molecular dynamics with a nonadditive intermolecular potential

The temperature and density dependence of the structure and polarization properties of bulk water were systematically investigated using the ab initio MCYna potential [Li et al., J. Chem. Phys. 127, 154509 (2007)], which includes nonadditive contributions to intermolecular interactions. Molecular dynamics simulations were conducted for isochores of 1, 0.8, and 0.6 g/cm^{3} and temperatures from 278 to 750 K. Special attention was paid to the structural change of water in the range from the normal boiling point to supercritical temperatures. At temperatures below the normal boiling temperature, water exhibits a tetrahedral structure along the 0.8 and 0.6 g/cm^{3} isochores. A significant collapse of the hydrogen bonding network was observed at temperatures of 450, 550, and 650 K. The MCYna potential was able to successfully reproduce the experimental dielectric constant. The dielectric constant and average dipole moments decrease with increasing temperature and decreasing density due to weakened polarization. A comparison is also made with SPC-based models.

Percolation Transition in Supercritical Water: A Monte Carlo Simulation Study

Computer simulations of water have been performed on the canonical ensemble at 15 different molecular number densities, ranging from 0.006 to 0.018 A-3, along the supercritical isotherm of 700 K, in order to characterize the percolation transition in the system. It is found that the percolation transition occurs at a somewhat higher density than what is corresponding to the supercritical extension of the boiling line. We have shown that the fractal dimension of the largest cluster and the probability of finding a spanning cluster are the most appropriate properties for the location of the true percolation threshold. Thus, percolation transition occurs when the fractal dimension of the largest cluster reaches 2.53, and the probability of finding a cluster that spans the system in at least one dimension and in all the three dimensions reaches 0.97 and 0.65, respectively. On the other hand, the percolation threshold cannot be accurately located through the cluster size distribution, as it is distorted by appearance of clusters crossing the finite simulated system even far below the percolation threshold. The structure of the largest water cluster is dominated by a linear, chainlike arrangement, which does not change noticeably until the largest cluster becomes infinite.

Local Density Inhomogeneities and Dynamics in Supercritical Water: A Molecular Dynamics Simulation Approach

Molecular dynamics atomistic simulations in the canonical ensemble (NVT-MD) have been used to investigate the "Local Density Inhomogeneities and their Dynamics" in pure supercritical water. The simulations were carried out along a near-critical isotherm (Tr = T/Tc = 1.03) and for a wide range of densities below and above the critical one (0.2 rho(c) - 2.0 rho(c)). The results obtained reveal the existence of significant local density augmentation effects, which are found to be sufficiently larger in comparison to those reported for nonassociated fluids. The time evolution of the local density distribution around each molecule was studied in terms of the appropriate time correlation functions C(Delta)rhol(t). It is found that the shape of these functions changes significantly by increasing the density of the fluid. Finally, the local density reorganization times for the first and second coordination shell derived from these correlations exhibit a decreasing behavior by increasing the density of the system, signifying the density effect upon the dynamics of the local environment around each molecule.

Solvation in supercritical water

Solvation in supercritical water under equilibrium and nonequilibrium conditions is studied via molecular dynamics simulations. The influence of solute charge distributions and solvent density on the solvation structures and dynamics is examined with a diatomic probe solute molecule. It is found that the solvation structure varies dramatically with the solute dipole moment, especially in low-density water, in accord with many previous studies on ion solvation. This electrostrictive effect has important consequences for solvation dynamics. In the case of a nonequilibrium solvent relaxation, if there are sufficiently many water molecules close to the solute at the outset of the relaxation, the solvent response measured as a dynamic Stokes shift is almost completely governed by inertial rotations of these water molecules. By contrast, in the opposite case of a low local solvent density near the solute, not only rotations but also translations of water molecules play an important role in solvent relaxation dynamics. The applicability of a linear response is found to be significantly restricted at low water densities.

Spectroscopic characterization of microscopic hydrogen-bonding disparities in supercritical water

The local hydrogen-bonding environment in supercritical water (380 degrees C, 300 bars, density 0.54 gcm3) was studied by x-ray Raman scattering at the oxygen K edge. The spectra are compared to those of the gas phase, liquid surface, bulk liquid, and bulk ice, as well as to calculated spectra. The experimental model systems are used to assign spectral features and to quantify specific local hydrogen-bonding situations in supercritical water. The first coordination shell of the molecules is characterized in more detail with the aid of the calculations. Our analysis suggests that approximately 65% of the molecules in supercritical water are hydrogen bonded in configurations that are distinctly different from those in liquid water and ice. In contrast to liquid water the bonded molecules in supercritical water have four intact hydrogen bonds and in contrast to ice large variations of bond angles and distances are observed. The remaining approximately 35% of the molecules exhibit two free O-H bonds and are thus either not involved in hydrogen bonding at all or have one or two hydrogen bonds on the oxygen side. We determine an average O-O distance of 3.1+/-0.1 A in supercritical water for the H bonded molecules at the conditions studied here. This and the corresponding hydrogen bond lengths are shown to agree with neutron- and x-ray-diffraction data at similar conditions. Our results on the local hydrogen-bonding environment with mainly two disparate hydrogen-bonding configurations are consistent with an extended structural model of supercritical water as a heterogeneous system with small patches of bonded molecules in various tetrahedral configurations and surrounding nonbonded gas-phase-like molecules.