Research Article  Open Access
On Water Density Fluctuations with Helices of Hydrogen Bonds
Abstract
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 twostructural fluid by dynamic forming the two topological kinds of clusters as a consequence of condensedmatter density fluctuations. Another feature of the densewaterpart 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 wellestablished 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 Xray emission spectroscopy (XES) and Xray Raman spectroscopy (XRS, BL62), 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 welldefined H bonds (one donating and one accepting) is in a strong contrast to the accepted picture of four quasitetrahedral H bonds.
Very recently, also Xray 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 lonepair 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 smallangle Xray scattering (SAXS) studies, they have found evidence for density nonhomogeneity on a length scale of 12โ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 โringsandchainsโ 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 initioalmost uniformly support the standard model.
At the same time, it is known that a topological densityfluctuations 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 VoronoiDelaunay 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 densepart structure of liquids, they have selected those Delaunay simplexes which have a tetrahedral factor more than or equal to the some given one, (convention criterion). Then, it is possible to build atomic ensembles in the MD model of liquids and investigate their densepart structure by choosing assorted simplexes of the regular subdivision in snapshots of the MD cell [10โ13].
Opposite, a topological criterion offered in [13โ16] for finding simplexes of the atomicensemble dense part exactly (rigorous criterion) defines densepart 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 HBordered clusters of water molecules are formed and broken uninterruptedly [19].
MD studying density fluctuations of the liquid water have shown [18] that the lowdensity tetrahedralordered ice regions are divided by higher density tetrahedral clusters with ramified structure (see Figure 1). It is obvious that these densepart 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 (, ) of secondclass particles relative to coordinates (, ) of firstclass particles which are defined as 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 firstclass particle in the same units; is the angular correction of firstclass particle shift in the plane (, ); is the angular shift of the secondclass particle projection on the plane (, ) relative to the firstclass (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: = = and angles between them:
From the equality of angles, , , and , and noting , one can obtain a system of four independent equations for defining the parameters: , , :
Solutions of this task for and as functions of and are and . Then, for , we have , which has a single real solution, . From here, we obtain , , , and .
Using these parameters, it is easy to find the value of which gives .
Now, one can estimate a packing factor of water molecules into the helix provided that all particles are inside the cylinder in radius of , that is, when particles forming the helix belong to it entirely. Then, we have the bottom estimation of this factor as .
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 per the unit of intermolecular distance. From here, we will find.
One can estimate a volume portion of helical clusters as the dense part of liquid water from , taking into account that the density of ice is equal to 0.92โg/cm^{3} and the water one is equal to 1.0โg/cm^{3} at 4ยฐC. From here, we obtain an upper bound of this portion as helical clusters .
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 hydrogenbonded chains due to a collective proton motion [20โ22]. Another basis is electromagnetic inducing the collective proton motion [23].
For describing such coherent vibrations, we can use equation [24]: modified for oscillatory tunneling of heavy quantum particles (protons) through the potential barrier [25] with a permeability, : Here, and are, respectively, proton and electron masses is the light velocity is the selfinduction of closed circuit is its electrical capacity and is the dielectric water permeability. A solution of (5) presents electromagnetic waves with frequency of For defining its value, we shall estimate the parameters in (5). It is known that ~ 1836, ~ 80, ~ 3ยท10^{10}โcm/s. From defining and in [24], one can present them for given task as Here, a = 2.82 , , is the average radius of water molecule, , is the number of water molecules in helical cluster (see Figure 4), and is the number of coin coils which is approximately equal to . An expression of is obtained in [26] as where is equal to the length difference of covalent and hydrogen bonds, that is, = 0.84 , J, and 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 .
Finally, from (6), we find a spectral series of coherent vibrations of protons in helical clusters of the dense part of liquid water: for . 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 crystallineice clusters as a nondense part of liquid water. Therefore, one can tell about water as a twostructural 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 fluctuationinduced 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 twostructural 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.
Acknowledgments
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. 100800217a).
References
 C. Huang, K. T. Wikfeldt, T. Tokushima et al., โThe inhomogeneous structure of water at ambient conditions,โ Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 36, pp. 15214โ15218, 2009. View at: Publisher Site  Google Scholar
 P. Wernet, D. Nordlund, U. Bergmann et al., โThe structure of the first coordination shell in liquid water,โ Science, vol. 304, no. 5673, pp. 995โ999, 2004. View at: Publisher Site  Google Scholar
 E. Cartlidge, โThe strangest liquid,โ New Scientist, February 2010. View at: Google Scholar
 K. Tuttle, โResearchers rediscover the structure of water,โ SLAC Today, February 2010. View at: Google Scholar
 T. Tokushima, Y. Harada, O. Takahashi et al., โHigh resolution Xray emission spectroscopy of liquid water: the observation of two structural motifs,โ Chemical Physics Letters, vol. 460, no. 4–6, pp. 387โ400, 2008. View at: Publisher Site  Google Scholar
 J. D. Smith, C. D. Cappa, B. M. Messer, W. S. Drisdell, R. C. Cohen, and R. J. Saykally, โProbing the local structure of liquid water by Xray absorption spectroscopy,โ Journal of Physical Chemistry B, vol. 110, no. 40, pp. 20038โ20045, 2006. View at: Publisher Site  Google Scholar
 H. A. Stern and B. J. Berne, โQuantum effects in liquid water: pathintegral simulations of a flexible and polarizable ab initio model,โ Journal of Chemical Physics, vol. 115, no. 16, pp. 7622โ7628, 2001. View at: Publisher Site  Google Scholar
 M. J. McGrath, J. I. Siepmann, I.F. W. Kuo et al., โSimulating fluidphase equilibria of water from first principles,โ Journal of Physical Chemistry A, vol. 110, no. 2, pp. 640โ646, 2006. View at: Publisher Site  Google Scholar
 J. D. Bernal, โThe structure of liquids,โ Proceedings of the Royal Society of London. Series A, vol. 280, no. 1382, pp. 299โ322, 1964. View at: Google Scholar
 Y. Waseda, The structure of noncrystalline materials; liquids and amorphous solids, McGrawHill, New York, NY, USA, 1980.
 N. N. Medvedev, VoronoyDelaunay method in research of structure of noncrystalline systems Novosibirsk, Publishing House of the Siberian Branch of the Russian Academy of Science, 2000.
 N. P. Smolin, B. R. Gelchinski, A. A. Mirzoev, and E. V. Dyuldina, โThe analysis of the shortrangeorder atomic structure of liquid metals by the Voronoi polyhedron method and a check of the adequacy of the results gained by the RMC method,โ Journal of NonCrystalline Solids, vol. 312–314, pp. 90โ94, 2002. View at: Publisher Site  Google Scholar
 A. S. Kolokol et al., โMDsimulation of liquid metals in the substantiation of a concept for designing liquid metal coolants,โ Preprint IAE6281/11. Moscow, Russia, 2003. View at: Google Scholar
 A. S. Kolokol and A. L. Shimkevich, โTopological structure of liquid metals,โ Atomic Energy, vol. 98, no. 3, pp. 187โ190, 2005. View at: Publisher Site  Google Scholar
 A. L. Shimkevich, A. S. Kolokol, and I. Yu. Shimkevich, โTwostructure model for simple metals,โ Journal of NonCrystalline Solids, vol. 353, no. 3240, pp. 3472โ3474, 2007. View at: Publisher Site  Google Scholar
 A. S. Kolokol and A. L. Shimkevich, โTopological structure of density fluctuations in condensed matter,โ Journal of NonCrystalline Solids, vol. 356, no. 45, pp. 220โ223, 2010. View at: Publisher Site  Google Scholar
 A. L. Shimkevich, โTetrahedralchaincluster model for thermodynamic description of fluids,โ in Proceedings of the 16th International Conference on Nuclear Engineering (ICONE '08), Orlando, Fla, USA, 2008, Paper ICONE1648566. View at: Google Scholar
 P. N. Alekseev et al., โOn basic principles for modifying water as a coolant of PWR,โ in Transactions of European Nuclear Conference, vol. 4, p. 9, 2010. View at: Google Scholar
 R. Roy, W. A. Tiller, I. Bell, and M. R. Hoover, โThe structure of liquid water; novel insights from materials research; potential relevance to homeopathy,โ Materials Research Innovations, vol. 9, no. 4, pp. 98โ103, 2005. View at: Google Scholar
 G. Zundel and K. Leberle, โProton polarizability of poly(Ltyrosine)hydrogen phosphate hydrogen bonds as a function of alkali cations,โ Biopolymers, vol. 23, no. 4, pp. 695โ705, 1984. View at: Publisher Site  Google Scholar
 G. Zundel, โHydrogen bonds with large proton polarizability and proton transfer processes in electrochemistry and biology,โ Advances in Chemical Physics, vol. 111, pp. 1โ217, 2000. View at: Google Scholar
 M. Shaikh, J. Mohanty, D. K. Maity, S. K. Nayak, and H. Pal, โCollective proton motion in the intramolecular hydrogen bonding network and the consequent enhancement in the acidity of hydroxycalixarenes,โ Journal of Photochemistry and Photobiology A, vol. 195, no. 1, pp. 116โ126, 2008. View at: Publisher Site  Google Scholar
 G. Careri, โCorrelations & conductivity: geometric aspects of physics chemistry and biology,โ in Proceedings of the NATO Advanced Study Institute Held at Cargese, Corsica, France, H. Stanley and N. Ostrowsky, Eds., Kluwer Academic Publishers, Dordrecht, The Netherlands, 1990. View at: Google Scholar
 L. D. Landau and E. M. Lifshits, Electric Dynamics of Condensed Matter, GosTechLit, Moscow, Russia, 1957.
 K. HyeonDeuk and K. Ando, โQuantum effects of hydrogen atoms on the dynamical rearrangement of hydrogenbond networks in liquid water,โ Journal of Chemical Physics, vol. 132, no. 16, Article ID 164507, 2010. View at: Publisher Site  Google Scholar
 A. L. Shimkevich, The Composition Principles for Designing NuclearReactor materials, N. N. PonomarevStepnoi, Ed., IzdAt, Moscow, Russia, 2008.
Copyright
Copyright © 2011 Alexander Shimkevich and Inessa Shimkevich. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.