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, КT<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 [10–13].
Opposite, a topological criterion offered in [13–16] 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.
The dense-part cluster in MD model of water at 300 K and its frame (broken red line); blue points are the molecules and red points are centers of the cluster tetrahedrons.
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 (x2, y2) of second-class particles relative to coordinates (x1, y1) of first-class particles which are defined as
r⃗k1:{xk1=aρcos[(4π5-β)k]yk1=aρsin[(4π5-β)k]zk1=a∆k,r⃗k2:{xk2=aρcos[(4π5-β)k-γ]yk2=aρsin[(4π5-β)k-γ]zk2=a∆k+a1-4ρ2sin2(γ2).
We obtain a helix of hydrogen bonds with the radius, ρ, in units of HB distance, a; ∆ is the shift value (on the helix axis “z”) of the following first-class particle in the same units; β is the angular correction of first-class particle shift in the plane (x, y); γ is the angular shift of the second-class particle projection on the plane (x, y) relative to the first-class (see Figure 2) with the same index, k.
The diagram of angular shift, γ, for coordinates (x2, y2) of second-class particles relative to coordinates (x1, y1) of first-class particles.
These parameters are defined under the assumption of helix edges equality as the distance between oxygen atoms of nearest water molecules: a = |r⃗k1-r⃗k-22| = |r⃗k+11-r⃗k-22| and angles between them:cosα1=(1a2)(r⃗k1-r⃗k-22)(r⃗k+11-r⃗k-22),cosα2=(1a2)(r⃗k2-r⃗k1)(r⃗k-22-r⃗k1),cosα3=(1a2)(r⃗k+12-r⃗k+11)(r⃗k-22-r⃗k+11).
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=∆2-∆1-4ρ2sin2(γ2)+ρ2sinαsin(α-γ),0=9∆2-6∆1-4ρ2sin2(γ2)+ρ2[2cosγ+2cos(3α-γ)],0=6∆2-3∆1-4ρ2sin2(γ2)+ρ2[-cosα-cos2α+cosγ+cos(3α-γ)],0=∆1-4ρ2sin2(γ2)+ρ2[cos2α+cos3α-cosγ+2sinαsin(α-γ)-cos(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.189685°≈4π/59, cosγ=58/63, ∆=7/190, and ρ=13/190.
Using these parameters, it is easy to find the value of cosαi=-11/38 which gives αi=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≥(2k+4)/k∆π(ρ+0.5)2k→∞=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 cubic model of ice.
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.
The dense part of water density fluctuations with the coherent exchange of protons (red balls) between oxygen atoms (blue balls) along the helical hydrogen bonds (black lines) and HB bridges (pink balls); the projections of oxygen atoms in the (x, y) plane are denoted by numbers.
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 [20–22]. Another basis is electromagnetic inducing the collective proton motion [23].
For describing such coherent vibrations, we can use equation [24]: 1c2Ld2qdt2+qεC=0,
modified for oscillatory tunneling of heavy quantum particles (protons) through the potential barrier [25] with a permeability, D: mpmec2Ld2qdt2+DqεC=0.
Here, mp and me are, respectively, proton and electron masses c is the light velocity L is the self-induction of closed circuit C is its electrical capacity and ε is the dielectric water permeability. A solution of (5) presents electromagnetic waves with frequency of ν=cmeD/εmpLC/2π.
For defining its value, we shall estimate the parameters in (5). It is known that mp/me ~ 1836, ε ~ 80, c ~ 3·1010 cm/s. From defining L and C in [24], one can present them for given task as L=πa2ρ2N2/l,C=rwl/2(l-rw).
Here, a = 2.82 Å, rw=1.41Å, is the average radius of water molecule, l=naΔ, n is the number of water molecules in helical cluster (see Figure 4), and N is the number of coin coils which is approximately equal to n/6. An expression of D is obtained in [26] as D=[1+sh2(Λk2)(k12+k22)24k12k22]-1,
where k1=2mpE/ℏ,k2=2mp(Uo-E)/ℏ,Λ is equal to the length difference of covalent and hydrogen bonds, that is, Λ = 0.84 Å, E=kBT=0.42·10-20 J, and Uo=1.66·10-20 J, as one fourth energy difference of covalent and hydrogen bonds, kB is Boltzmann constant, and T is Kelvin temperature.
In substituting these values in formula (8), we shall obtain D~1.1·10-4.
Finally, from (6), we find a spectral series of coherent vibrations of protons in helical clusters of the dense part of liquid water: ν~22n-1n[THz],
for n≥6. Then, we will receive the spectral absorption modes of coherent vibrating protons (Table 1).
The spectral absorption modes of the coherent proton vibrations as a function of the water-molecules number in the helix, n.
n
Spectral modes
n
Spectral modes
THz
cm−1
THz
cm−1
6
8.2
273
20
4.8
160
10
6.6
220
25
4.3
144
15
5.5
183
30
3.9
132
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.
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. 10-08-00217a).
HuangC.WikfeldtK. T.TokushimaT.NordlundD.HaradaY.BergmannU.NiebuhrM.WeissT. M.HorikawaY.LeetmaaM.LjungbergM. P.TakahashiO.LenzA.OjamäeL.LyubartsevA. P.ShinS.PetterssonL. G. M.NilssonA.The inhomogeneous structure of water at ambient conditions20091063615214152182-s2.0-7034932616610.1073/pnas.0904743106WernetP.NordlundD.BergmannU.CavalleriM.OdeliusN.OgasawaraH.NäslundL. Å.HirschT. K.OjamäeL.GlatzelP.PetterssonL. G. M.NilssonA.The structure of the first coordination shell in liquid water200430456739959992-s2.0-234255805310.1126/science.1096205CartlidgeE.The strangest liquidFebruary 2010TuttleK.Researchers rediscover the structure of waterFebruary 2010TokushimaT.HaradaY.TakahashiO.SenbaY.OhashiH.PetterssonL. G. M.NilssonA.ShinS.High resolution X-ray emission spectroscopy of liquid water: the observation of two structural motifs20084604–63874002-s2.0-4894911611610.1016/j.cplett.2008.04.077SmithJ. D.CappaC. D.MesserB. M.DrisdellW. S.CohenR. C.SaykallyR. J.Probing the local structure of liquid water by X-ray absorption spectroscopy20061104020038200452-s2.0-3375036662210.1021/jp063661cSternH. A.BerneB. J.Quantum effects in liquid water: path-integral simulations of a flexible and polarizable ab initio model2001115167622762810.1063/1.1407287McGrathM. J.SiepmannJ. I.siepmann@chem.umn.eduKuoI.-F. W.MundyC. J.VandevondeleJ.HutterJ.MohamedF.KrackM.Simulating fluid-phase equilibria of water from first principles2006110264064610.1021/jp0535947BernalJ. D.The structure of liquids19642801382299322WasedaY.1980New York, NY, USAMcGraw-HillMedvedevN. N.2000Publishing House of the Siberian Branch of the Russian Academy of ScienceSmolinN. P.GelchinskiB. R.MirzoevA. A.DyuldinaE. V.The analysis of the short-range-order atomic structure of liquid metals by the Voronoi polyhedron method and a check of the adequacy of the results gained by the RMC method2002312–31490942-s2.0-003678873410.1016/S0022-3093(02)01655-1KolokolA. S.MD-simulation of liquid metals in the substantiation of a concept for designing liquid metal coolantsPreprint IAE-6281/11. Moscow, Russia, 2003KolokolA. S.ShimkevichA. L.Topological structure of liquid metals200598318719010.1007/s10512-005-0191-9ShimkevichA. L.shall@dhtp.kiae.ruKolokolA. S.ShimkevichI. Yu.Two-structure model for simple metals200735332-403472347410.1016/j.jnoncrysol.2007.05.098KolokolA. S.ShimkevichA. L.Topological structure of density fluctuations in condensed matter20103564-52202232-s2.0-7454920475910.1016/j.jnoncrysol.2009.11.020ShimkevichA. L.Tetrahedral-chain-cluster model for thermodynamic description of fluidsProceedings of the 16th International Conference on Nuclear Engineering (ICONE '08)2008Orlando, Fla, USAPaper ICONE16-48566AlekseevP. N.On basic principles for modifying water as a coolant of PWR4Transactions of European Nuclear Conference20109RoyR.TillerW. A.BellI.HooverM. R.The structure of liquid water; novel insights from materials research; potential relevance to homeopathy200594981032-s2.0-30344462987ZundelG.LeberleK.Proton polarizability of poly(L-tyrosine)-hydrogen phosphate hydrogen bonds as a function of alkali cations198423469570510.1002/bip.360230409ZundelG.Hydrogen bonds with large proton polarizability and proton transfer processes in electrochemistry and biology200011112172-s2.0-22844455879ShaikhM.MohantyJ.MaityD. K.NayakS. K.PalH.Collective proton motion in the intramolecular hydrogen bonding network and the consequent enhancement in the acidity of hydroxycalixarenes200819511161262-s2.0-3914910642010.1016/j.jphotochem.2007.09.013CareriG.StanleyH.OstrowskyN.Correlations & conductivity: geometric aspects of physics chemistry and biology1990Dordrecht, The NetherlandsKluwer Academic PublishersLandauL. D.LifshitsE. M.1957Moscow, RussiaGosTechLitHyeon-DeukK.AndoK.Quantum effects of hydrogen atoms on the dynamical rearrangement of hydrogen-bond networks in liquid water2010132162-s2.0-7795235775510.1063/1.3397809164507ShimkevichA. L.2008Moscow, RussiaIzdAtN. N. Ponomarev-Stepnoi, Ed.