An adaptive model is developed here for the liquid water density fluctuations as momentary dense clusters with helices of hydrogen bonds and nondense tetrahedral clusters of ice. This model can be useful for explanation of liquid water structural anomalies including the high quantity of hydrogen bonds with quasitetrahedral orientation in the nonordered liquid water. The topology of such clusters is essentially differed from the one of the crystalline ice. From this and only this point of view, the liquid water can be considered as a two-structural fluid by dynamic forming the two topological kinds of clusters as a consequence of condensed-matter density fluctuations. Another feature of the dense-water-part clusters is helical ordering of protons which can realize coherent vibrations. A spectral series of such vibrations is determined as a function of the number of molecules into the helical cluster.

1. Introduction

Although the liquid water is all around us, we still do not understand how it behaves on a microscopic level. Its dynamic hydrogen bonds (HB) structure has been the subject of intense debate for decades. Ice with well-established HB structure forms a tight tetrahedral lattice of water molecules each binding to four others. The prevailing model of the liquid water holds that the molecules loosen their grip but remain generally arranged in the same tetrahedral groups [1].

However, the majority of water molecules in the liquid are found in higher density regions [2]. The structure and states of the molecules in these regions are found through X-ray emission spectroscopy (XES) and X-ray Raman spectroscopy (XRS, BL6-2), techniques which allows to analyze electronic states of molecules. The researchers have found that the molecules in the high density regions form an asymmetrical disordered structure where some islands of tetrahedral order are floated. This new research is requiring some sort of secondary structure. The greater density of liquid water implies that the molecules are more closely packed than the simple tetrahedrons seen in the ice [3, 4].

The conclusion that a dominant fraction of molecules in the liquid water is very asymmetrically hydrogen bonded with only two well-defined H bonds (one donating and one accepting) is in a strong contrast to the accepted picture of four quasitetrahedral H bonds.

Very recently, also X-ray emission spectroscopy (XES) is applied at high resolution to the liquid water [5], arriving at the similar conclusions as [1, 2]. They interpret the observed splitting of the lone-pair state into two sharp features as evidence of a bimodal distribution of local structures in the liquid: either strongly tetrahedral or highly disordered, which interconvert discontinuously with temperature. From small-angle X-ray scattering (SAXS) studies, they have found evidence for density nonhomogeneity on a length scale of 1-2 nm indicating that the two components are spatially separated on the time scale of the experiment [3, 4].

The controversy about the microstructure of liquid water pits a new model involving water molecules in a relatively stable “rings-and-chains” structure against the standard model that considers water molecules in the distorted tetrahedral coordination. An analysis [6] and molecular dynamics (MD) simulations [7, 8]—both classical and ab initio-almost uniformly support the standard model.

At the same time, it is known that a topological density-fluctuations structure of liquids is ramified clusters of almost regular Delaunay simplexes (tetrahedrons) built on the fours of densely packed particles and connected in pairs by faces in tetrahedral Bernal chains [9, 10]. The Voronoi-Delaunay method is used for analyzing atomic configurations in disordered systems. The review of publications on this subject is presented in the monograph [11].

For identifying the dense-part structure of liquids, they have selected those Delaunay simplexes which have a tetrahedral factor more than or equal to the some given one, К𝑇<1 (convention criterion). Then, it is possible to build atomic ensembles in the MD model of liquids and investigate their dense-part structure by choosing assorted simplexes of the regular subdivision in snapshots of the MD cell [1013].

Opposite, a topological criterion offered in [1316] for finding simplexes of the atomic-ensemble dense part exactly (rigorous criterion) defines dense-part simplexes on the overall length of their edges in maximum of number of tetrahedral clusters in the MD cell. Then, one can get a statistics of these clusters for different liquids.

These data may be useful for creating an analytical model of the cluster statistics for calculating the configuration entropy of any condensed matter and estimating its thermodynamic properties using only a pair potential of particles [17] and for analyzing the microstructure features of water [18].

2. Same Structural Features of the Liquid Water

Precise experimental techniques for determining the local structure of the liquid water are lacking since each water molecule undergoes a quick rearrangement (during femtoseconds). The need for better understanding the liquid water at the microscopic level has forced developing computational methods that describe the individual and cooperative dynamics and structure of water molecules. Many investigations are carried out using these techniques. They show that locally HB-ordered clusters of water molecules are formed and broken uninterruptedly [19].

MD studying density fluctuations of the liquid water have shown [18] that the low-density tetrahedral-ordered ice regions are divided by higher density tetrahedral clusters with ramified structure (see Figure 1). It is obvious that these dense-part clusters of water density fluctuations are more complicated due to hydrogen bonds.

3. Helical Model for the Dense Part of Liquid Water

One can expect that water molecules have enough time for rebuilding hydrogen bonds in the dense part of water density fluctuations due to fast librations and rotations of them. From this, we will build a model for tetrahedral clusters of the dense part of liquid water provided that each water molecule in such cluster is hydrogen bonded with three molecules of this cluster and all angles between these bonds are equal to each other.

It turned out that such cluster represents “nanotube” of water molecules and geometric parameters of this tube are defined strictly. For illustration, Figure 2 shows an angular shift, 𝛾, for coordinates (𝑥2, 𝑦2) of second-class particles relative to coordinates (𝑥1, 𝑦1) of first-class particles which are defined as 𝑟1𝑘𝑥1𝑘=𝑎𝜌cos4𝜋5𝑘𝑦𝛽1𝑘=𝑎𝜌sin4𝜋5𝑘𝑧𝛽1𝑘=𝑎𝑘,𝑟2𝑘𝑥2𝑘=𝑎𝜌cos4𝜋5𝑦𝛽𝑘𝛾2𝑘=𝑎𝜌sin4𝜋5𝑧𝛽𝑘𝛾2𝑘=𝑎𝑘+𝑎14𝜌2sin2𝛾2.(1) We obtain a helix of hydrogen bonds with the radius, 𝜌, in units of HB distance, 𝑎; is the shift value (on the helix axis “𝑧”) of the following first-class particle in the same units; 𝛽 is the angular correction of first-class particle shift in the plane (𝑥, 𝑦); 𝛾 is the angular shift of the second-class particle projection on the plane (𝑥, 𝑦) relative to the first-class (see Figure 2) with the same index, 𝑘.

These parameters are defined under the assumption of helix edges equality as the distance between oxygen atoms of nearest water molecules: 𝑎 = |𝑟1𝑘𝑟2𝑘2| = |𝑟1𝑘+1𝑟2𝑘2| and angles between them:cos𝛼1=1𝑎2𝑟1𝑘𝑟2𝑘2𝑟1𝑘+1𝑟2𝑘2,cos𝛼2=1𝑎2𝑟2𝑘𝑟1𝑘𝑟2𝑘2𝑟1𝑘,cos𝛼3=1𝑎2𝑟2𝑘+1𝑟1𝑘+1𝑟2𝑘2𝑟1𝑘+1.(2)

From the equality of angles, 𝛼1, 𝛼2, and 𝛼3, and noting 𝛼𝜋/5+𝛽, one can obtain a system of four independent equations for defining the parameters: 𝛼, 𝛾, и𝜌:0=214𝜌2sin2𝛾2+𝜌2sin𝛼sin(𝛼𝛾),0=92614𝜌2sin2𝛾2+𝜌2[],2cos𝛾+2cos(3𝛼𝛾)0=62314𝜌2sin2𝛾2+𝜌2[],cos𝛼cos2𝛼+cos𝛾+cos(3𝛼𝛾)0=14𝜌2sin2𝛾2+𝜌2[].cos2𝛼+cos3𝛼cos𝛾+2sin𝛼sin(𝛼𝛾)cos(3𝛼𝛾)(3)

Solutions of this task for 𝛾 and 𝜌 as functions of 𝛼 and Δ are tg[(𝛼𝛾)/2]=(cos𝛼+cos2𝛼)/(sin2𝛼+2sin𝛼) and 𝜌2=(3/2)2/(cos𝛼+cos2𝛼). Then, for 𝛼, we have 0=12cos3𝛼+16cos2𝛼4cos𝛼8, which has a single real solution, cos𝛼=2/3. From here, we obtain 𝛽=12.1896854𝜋/59, cos𝛾=58/63, =7/190, and 𝜌=13/190.

Using these parameters, it is easy to find the value of cos𝛼𝑖=11/38 which gives 𝛼𝑖=106.83.

Now, one can estimate a packing factor of water molecules into the helix provided that all particles are inside the cylinder in radius of 𝜌+0.5, that is, when particles forming the helix belong to it entirely. Then, we have the bottom estimation of this factor as 𝜂helix(2𝑘+4)/𝑘𝜋(𝜌+0.5)2𝑘=2/𝜋(𝜌+0.5)2=0.70.

In using the cubic model of ice shown in Figure 3, it is easy to estimate its packing factor for water molecules.

The elementary ice cell contains 7.5 molecules of water. The edge length of such cell is equal to 4/3 per the unit of intermolecular distance. From here, we will find𝜂ice=7.5(3/4)3=0.61.

One can estimate a volume portion Ω of helical clusters as the dense part of liquid water from 1.0=(0.92/0.61)[0.70Ω+0.61(1Ω)], taking into account that the density of ice is equal to 0.92 g/cm3 and the water one is equal to 1.0 g/cm3 at 4°C. From here, we obtain an upper bound of this portion as helical clusters Ω(0.08/0.92)(0.61/0.09)=0.59.

These clusters having higher molecular density provide, in principle, four hydrogen bonds for each molecule of water: three of them are internal and one is external for connecting ice crystallites with less density. However, if protons are ordered in the helix (see Figure 4), there can arise resonant vibrations of proton current in a coin of the helical cluster. One can estimate frequencies of these vibrations.

4. Evaluating Resonance Frequencies of the Coherent Proton Vibrations

One assumes that the motion of all protons in helical circuit (coin), as shown in Figure 4, is synchronous; that is, the protons vibrate coherently in the helical cluster of hydrogen bonds for the dense part of liquid water. This hypothesis is based on anomalous proton polarizability in hydrogen-bonded chains due to a collective proton motion [2022]. Another basis is electromagnetic inducing the collective proton motion [23].

For describing such coherent vibrations, we can use equation [24]: 1𝑐2𝐿𝑑2𝑞𝑑𝑡2+𝑞𝜀𝐶=0,(4) modified for oscillatory tunneling of heavy quantum particles (protons) through the potential barrier [25] with a permeability, 𝐷: 𝑚𝑝𝑚𝑒𝑐2𝐿𝑑2𝑞𝑑𝑡2+𝐷𝑞𝜀𝐶=0.(5) Here, 𝑚𝑝 and 𝑚𝑒 are, respectively, proton and electron masses 𝑐 is the light velocity 𝐿 is the self-induction of closed circuit 𝐶 is its electrical capacity and 𝜀 is the dielectric water permeability. A solution of (5) presents electromagnetic waves with frequency of 𝜈=𝑐𝑚e𝐷/𝜀𝑚𝑝𝐿𝐶/2𝜋.(6) For defining its value, we shall estimate the parameters in (5). It is known that 𝑚𝑝/𝑚𝑒 ~ 1836, 𝜀 ~ 80, 𝑐 ~ 3·1010 cm/s. From defining 𝐿 and 𝐶 in [24], one can present them for given task as 𝐿=𝜋a2𝜌2𝑁2/𝑙,𝐶=𝑟w𝑙/2𝑙𝑟w.(7) Here, a = 2.82 Å, 𝑟w=1.41Å, is the average radius of water molecule, 𝑙=𝑛aΔ, 𝑛 is the number of water molecules in helical cluster (see Figure 4), and 𝑁 is the number of coin coils which is approximately equal to 𝑛/6. An expression of 𝐷 is obtained in [26] as 𝐷=1+sh2Λ𝑘2𝑘21+𝑘2224𝑘21𝑘221,(8) where 𝑘1=2𝑚𝑝E/,𝑘2=2𝑚𝑝(U𝑜E)/,Λ is equal to the length difference of covalent and hydrogen bonds, that is, Λ = 0.84 Å, 𝐸=𝑘𝐵𝑇=0.421020 J, and 𝑈𝑜=1.661020 J, as one fourth energy difference of covalent and hydrogen bonds, 𝑘𝐵 is Boltzmann constant, and 𝑇 is Kelvin temperature.

In substituting these values in formula (8), we shall obtain 𝐷1.1104.

Finally, from (6), we find a spectral series of coherent vibrations of protons in helical clusters of the dense part of liquid water: 𝜈22𝑛1𝑛[]THz,(9) for 𝑛6. Then, we will receive the spectral absorption modes of coherent vibrating protons (Table 1).

5. Discussion of Results

The topology of the dense liquid water part as helices of hydrogen bonds differs from the one for tetrahedral crystalline-ice clusters as a nondense part of liquid water. Therefore, one can tell about water as a two-structural liquid though this is a dynamic distinction having only statistical sense because such clusters arise and disappear as density fluctuations of the liquid water.

At the same time, possible coherent proton vibrations by external impact at their resonant frequencies can selectively amplify some modes and, thus, strengthen water microheterogeneity far from its equilibrium. In this connection, studying the electromagnetic absorption spectra in this range of frequencies can be interested both for revealing spectral lines (9) and for searching for a possible technique for controlling the microstructure of the liquid water.

Moreover, any additives can be concentrated in the expected ramified border between two dynamic structures of water solution that, in turn, can cause fluctuation-induced clustering and impurity [26].

Thus, the helices of hydrogen bonds in water density fluctuations and the high resonant mobility of proton in them as a catalytic stimulator for chemical reactions can create templates for building macromolecules and structural blocks of biosystems.

6. Conclusions

The model of hydrogen bond helix is built for the dense part of the liquid water density fluctuations. This model allows providing the high quantity of tetrahedral hydrogen bonds in the nonordered liquid matrix and explaining the density anomaly of liquid water.

The topology of such helical clusters is essentially differed from the one of the crystalline ice. From this and only this point of view, the liquid water can be considered as the two-structural liquid because the formation and decay of such clusters has dynamic character and is natural consequence of water density fluctuations.

The other interesting feature of these clusters is possible helical ordering protons in them which can have coherent vibrations. A spectral series of such resonances is defined.


The authors are pleased to acknowledge Dr. A. S. Kolokol for giving some data on molecular dynamic simulation of water structure and discussing this work which is supported by the Russian Foundation of Basic Research (RFBR) (Grant no. 10-08-00217a).