Ab initio study of water. II. Liquid structure, electronic and thermodynamic properties over a wide range of temperature and density

The electronic and liquid structures of water and its thermodynamic propert ies are studied over a wide range of temperature (0 degrees-600 degrees C) and density (0.6-1.4 g/cm(3)) based on the ab initio molecular orbital theo ry combined with the integral equation method of liquid. Unlike standard tr eatments of water by means of the classical statistical mechanics including molecular simulations, the effective charges on oxygen and hydrogen atoms in water molecules are not "input parameters,'' but naturally predicted fro m the theory in the course of self-consistent determination of the electron ic structure and the density pair correlation functions in liquids. It is f ound that the molecular dipole moments and electronic polarization energies decrease with increasing temperature and/or density. The theoretical resul ts for dipole moments are in quantitative accord with the experimental data , which has been determined based on the NMR chemical shift coupled with th e molecular dynamics simulation [N. Matsubayashi, C. Wakai, and M. Nakahara , J. Chem. Phys. 110, 8000 (1999)]. The state dependence of the electronic structure is discussed in terms of the thermal activation of molecules and intermolecular interactions including the hydrogen bonds. The liquid struct ure of water is discussed in the wide range of thermodynamic states in term s of S(r), an average response of the pair correlation functions to tempera ture change which is introduced in the present study in order to make struc tural characteristics of water more distinctive. It is concluded from the b ehavior of the function that the short range structure of water retains the characteristics to ice, or the tetrahedral configuration, over relatively wide range of temperature in the normal density (1.0 g/cm(3)). The ice-like characteristics of water disappear to large extent both at high (1.4 g/cm( 3)) and low (0.6 g/cm(3)) densities for different reasons: in the high dens ity due to the packing effect, while in the low density due to essentially entropic cause, or increased configuration space available to a molecule. T he distance between the nearest-neighbor molecules in water are insensitive to the density change compared with those corresponding to the Lennard-Jon es fluid. The difference in the behaviors between the two fluids is explain ed in terms of the intermolecular interactions and liquid structures. The n umber of hydrogen bonds is calculated from the oxygen-hydrogen pair correla tion function using a new definition based on S( r), which enables us to di stinguish the hydrogen-bonded OH pairs from those just in contact due to pa cking effect. The temperature and density dependence of the quantity is dis cussed in terms of the liquid structure of water. (C) 1999 American Institu te of Physics. [S0021-9606(99)50142-6].