Report

The Structure of the First Coordination Shell in Liquid Water

See allHide authors and affiliations

Science  14 May 2004:
Vol. 304, Issue 5673, pp. 995-999
DOI: 10.1126/science.1096205

Abstract

X-ray absorption spectroscopy and x-ray Raman scattering were used to probe the molecular arrangement in the first coordination shell of liquid water. The local structure is characterized by comparison with bulk and surface of ordinary hexagonal ice Ih and with calculated spectra. Most molecules in liquid water are in two hydrogen–bonded configurations with one strong donor and one strong acceptor hydrogen bond in contrast to the four hydrogen–bonded tetrahedral structure in ice. Upon heating from 25°C to 90°C, 5 to 10% of the molecules change from tetrahedral environments to two hydrogen–bonded configurations. Our findings are consistent with neutron and x-ray diffraction data, and combining the results sets a strong limit for possible local structure distributions in liquid water. Serious discrepancies with structures based on current molecular dynamics simulations are observed.

Experimental studies of the hydrogen-bonded network structure in water have mainly relied on neutron and x-ray diffraction and infrared (IR) spectroscopies (1). Diffraction data from noncrystalline materials provide radial distribution functions (RDFs) (24) that do not provide angular correlations needed to uniquely assign local geometries in water (57). A more detailed atomistic picture has been derived theoretically by molecular dynamics (MD) simulations (4, 8) that are consistent with diffraction data. Structural information from IR spectroscopies generally relies on the correlation between the O-H stretching frequency and hydrogen-bond (H-bond) length, which has been shown to be ambiguous for liquid water (9).

Here, we report an independent experimental investigation of local bonding configurations in the first coordination shell of liquid water by using the near-edge fine structure in x-ray absorption spectroscopy (XAS), also denoted XANES and NEXAFS, where a core electron is excited into empty electronic states. The character of these states and, hence, the near-edge fine structure in XAS depends on the chemical environment, bond lengths, and bond angles (10). We also obtained the same spectral information by using nonresonant x-ray Raman scattering (XRS) involving core excitations (11).

XAS and XRS at the oxygen K-edge (12) are sensitive to distortions of H-bonds on the H-sides (donor H-bonds) of the molecules in the condensed phases of water (11, 13, 14). Because the time scale for excitation is much faster than the molecular (vibrational) motions in the liquid, these spectroscopies probe the electronic structure of a distribution of instantaneous configurations and thus allow decomposition in terms of specific H-bond situations (12). We analyzed the near-edge structures in the liquid water XA spectrum (the terminology “XA spectrum” is used for both XAS and XRS) with the aid of experimental model systems and calculated spectra. The XA spectra for water molecules in different H-bonding configurations are depicted in Fig. 1, where ice Ih bulk and surface spectra are compared with spectra of bulk liquid water at two temperatures. Bulk ice Ih is tetrahedrally coordinated, but the exact H-bonding environment at the ice Ih surface still raises questions (15, 16). However, there is consensus that a large fraction (50% or more) of the molecules in the first half bilayer of the ice Ih surface has one free O-H group, whereas the other is H-bonded to the second half bilayer. The liquid water XA spectrum closely resembles that for the ice surface, but it is very different from that of bulk ice. We interpret this finding, and our analysis demonstrates that the molecules in the liquid are not predominantly four-coordinated.

Fig. 1.

(a) Bulk ice Ih (XAS secondary electron yield). (b) Ice Ih surface (topmost surface layer, XAS Auger electron yield). (c) NH3-terminated first half bilayer of the ice Ih surface (XAS Auger electron yield). (d) Liquid water at ambient conditions [XAS fluorescence yield, taken from (13) and additionally corrected for saturation effects]. (e) Bulk liquid water at 25°C (solid line) and 90°C (dashed line) (XRS, spectra normalized to the same area). (f) Difference spectra: 25°C water minus bulk ice (solid curve) and 90°C water minus 25°C water (circles with error bars). The latter difference spectrum has been multiplied by a factor of 10. Note that XRS and XAS essentially give the same information, as indicated by the similarity of the room-temperature water spectra (e) and (d). Their slight differences are due to the poorer energy resolution of XRS. The temperature effects can best be studied with XRS because of the larger probing depth and the resulting easier sample preparation than for XAS. For details on sample preparation and probing depths, see (1113, 15).

The spectra in Fig. 1 can be divided into three main regions: the pre-edge (around 535 eV), the main edge (537 to 538 eV), and the post-edge (540 to 541 eV). The bulk ice spectrum (Fig. 1, curve a) is dominated by intensity in the post-edge region and shows a weak main-edge structure. Both the surface ice (Fig. 1, curve b) and liquid water (Fig. 1, curve d) spectra have a peak in the pre-edge region, a dominant main edge, and less intensity compared with bulk ice in the post-edge region. Termination of the ice surface with NH3 (Fig. 1, curve c) entails a coordination of the free O-H groups and causes the pre-edge peak to vanish and the intensity to shift to the post-edge region. We assign intensities in the pre- and main-edge regions to water molecules with one uncoordinated O-H group, whereas the intensity in the post-edge region is related to fully coordinated molecules. Remarkably, most molecules in bulk liquid water at room temperature exhibit a local coordination comparable to that at the ice surface, with one strong and one non-, or only weakly, H-bonded O-H group. The contribution to the spectrum from molecules with four-fold coordination similar to bulk ice is very small. Performing the measurements with D2O or H2O led to identical spectra within the experimental resolution, and thus tunneling contributions are not decisive.

Comparison of the XRS spectra of room-temperature (25°C) and hot water (90°C) (Fig. 1, curve e) shows that heating increases intensities in the pre- and main-edge regions while decreasing that in the post-edge, but the changes are small compared with the changes observed between ice and the liquid. Figure 1, curve f, shows the difference spectra of 25°C water minus ice (17) (solid curve) and 90°C minus 25°C water (circles with error bars). The latter has been multiplied by a factor of 10 to match the intensities. The similarity of the curves, in particular the same isosbestic point at 538.8 eV, supports a two-component structure in terms of O-H coordination: Configurations with one uncoordinated or weakly H-bonded O-H group replace tetrahedral ones in the ice-liquid phase transition, and heating liquid water causes similar types of changes, but of about one-tenth the magnitude. The contributions of the two different species at room temperature are 80% and 20% for species with only one strong H-bond and tetrahedrally coordinated, respectively; for 90°C, we find 85% and 15% with uncertainties of ±15 to ±20% in both cases (12).

We now compare this phenomenological finding with a detailed theoretical analysis based on density functional theory (DFT). The method used provides highly accurate results for both intensities and energy positions in x-ray spectra of H-bonded systems (12, 18). Here, the calculations provide model spectra for various distributions of bond lengths and angles as expected to occur in the liquid. To allow for calculations of a large number of local configurations, we use a small model cluster of 11 molecules, where a central molecule is surrounded by the first coordination shell (4 molecules) plus part of the second shell (6 molecules saturating the dangling O-H groups of the first shell). Spectra were calculated for the central molecule with systematically varied nearest neighbor H-bond environments.

First, the model cluster was tested by calculating the three known cases of bulk ice, ice surface, and NH3 terminated ice surface (Fig. 2). The calculated spectra based on clusters with 11 water molecules reproduce the general shapes and trends seen in the experimental spectra in Fig. 1: shift of intensity from post to pre and main edges when going from a fully coordinated molecule to one with a free O-H group (compare Fig. 2, curves a and b) and shift of intensity back to the post-edge region when the free O-H group is coordinated upon NH3 adsorption (compare Fig. 2, curves b and c).

Fig. 2.

Calculated spectra for the model systems bulk ice Ih. (a) Fully coordinated molecules. (b) Ice surface (molecules in the topmost surface layer with one free O-H). (c) NH3-terminated ice surface (molecules with an NH3-saturated O-H). Solid lines, small clusters; dashed lines, large clusters. Cluster sizes are 11 molecules for (a) to (c), solid lines; 44 molecules for (a), dashed line; and 27 molecules for (b) and (c), dashed lines. For the bulk ice cases, the calculated molecules are fully coordinated to their four nearest neighbors and surrounded by: (a, solid line) part of the second coordination shell (6 molecules) and (a, dashed line) the full second shell (12 molecules) and 27 further molecules. For the ice-surface cases, an appropriate number of molecules from the bulk ice clusters were removed [3 and 17 molecules for (b, solid line) and (b, dashed line), respectively] to generate configurations where the calculated molecule has one free O-H as found on the ice surface (12, 15, 16). To model the effect of NH3 termination, the same molecule was calculated after saturating the dangling O-H by an NH3 molecule (lone pair electrons of the N oriented toward the free O-H).

These findings did not depend substantially on the cluster size: Calculated spectra for identical nearest neighbor H-bond environments but now with a total of 44 (Fig. 2, curve a, dashed line) and 27 (Fig. 2, curves b and c, dashed lines) water molecules were similar to the spectra for small clusters. Furthermore, thermal fluctuations and/or vibrational motions did not affect the calculations considerably; for example, superposition of spectra with different intramolecular bonding distances was similar to that for the mean bonding configuration. The same holds for intermolecular vibrations (12).

To distinguish possible local configurations in the liquid, spectra were calculated for the central molecule in the 11-molecule cluster with systematically varied nearest neighbor H-bond distances and angles (Fig. 3A, curves a to h) (12). The spectrum for a configuration with tetrahedral coordination [linear H-bonds: all θ= 0° and 2.75 Å for all O-O distances (1, 2), Fig. 3A, curve a] shows a maximum close to 540 eV, similar to Fig. 2, curve a. Donor H-bond distortions were introduced by varying the O-O distance r and the angle θ of the nearest neighbor (and the two attached molecules; see inset of Fig. 3A) on only one H-side (Fig. 3A, curves b to e) and on both H-sides (Fig. 3A, curves f to h). Intermediate elongations smear out the maximum close to 540 eV (Fig. 3A, curve b). Upon further distortion of one donor H-bond, a pre-edge peak emerges at 535 eV and a main-edge structure develops at 537 eV, while post-edge intensities decrease (Fig. 3A, curves c to e). Additional distortion of the second donor H-bond (Fig. 3A, curves f and g) further reduces intensities above 540 eV without considerable change in the pre- and main-edge regions. Note that varying r at constant θ and varying θ at constant r can yield similar spectra (e.g., Fig. 3A, curves c and e).

Fig. 3.

(A) Calculated spectra for a cluster of 11 water molecules with O-O distances ri (in Å) and angles θi (in degrees) of the nearest neighbors i = 1, 2 on the H-sides (nearest neighbor H-bond acceptor molecules) systematically varied as given for each spectrum. Starting from (a) a tetrahedral structure with all four nearest neighbor O-O distances at 2.75 Å and linear H-bonds, for (b, c) only r1, for (d) r1 and θ1, and for (e) only θ1 are varied to introduce H-bond distortions on one H-side. For (f) θ1 and r2 and for (g, h) r1 and r2 are varied to introduce H-bond distortions on both H-sides. Compared with the cluster calculation for ice Ih in Fig. 2, curve a, solid line, relative orientations of the molecules have been changed arbitrarily to better account for the less ordered local configurations in the liquid. In addition, the calculated oscillator strengths are convoluted with a slightly larger width to account for the larger vibrational motion of the molecules in the liquid. This explains the differences between (a) and Fig. 2, curve a, solid line. (B) The systematic spectral changes allow for a definition of zones A and B to denote the positions of the two nearest neighbor H-bond acceptor molecules. Spectra show similar features for positions within one zone [compare, e.g., (a) and (b) or (c) to (e) and as verified by many further calculations], and the cases shown here can be considered to be representative. The boundary between the zones (black line) is used as a cutoff for H-bonding (12) as based on the occurrence of strong pre- and main-edge features (similar to the experimental ice surface spectrum) when the nearest neighbor is outside zone A beyond the black line (e.g., c to e).

Spectra c to e and g to h in Fig. 3A show that donor H-bond distortions of a certain degree manifest as a distinct pre-edge peak and an intense main edge in the oxygen K-edge XA spectrum of water. The analogy of these pre- and main-edge features to the ones observed for the ice-surface molecules with one broken donor H-bond (Figs. 1 and 2) suggests denoting the corresponding distorted H-bonds as “broken.” More specifically, based on spectral calculations of a large number of different configurations in the liquid, we define the cutoff for H-bonding in terms of O-O distances as the boundary between zones A and B (black line in Fig. 3B) (12). For linear H-bonds, we find the same cutoff as conventionally used in the analysis of MD simulations (19). For bent H-bonds, we use a smaller cutoff because of the occurrence of pre- and main-edge intensities for large angular distortions at standard O-O distances (Fig. 3A, curve e). This difference is not important, as will be addressed in detail below.

Zones A and B in Fig. 3B can be used to classify local water configurations: The environment in bulk ice, e.g., is an AA configuration, with both nearest neighbor H-bond acceptor molecules within zone A, whereas that at the ice surface is an AB configuration. We define three groups of configurations in the liquid: Double-donor (DD) configurations with two intact donor H-bonds are characterized by XA spectra without pre-edge peak (Fig. 3A, curves a, b, and f); single-donor (SD) configurations with one intact and one broken donor H-bond exhibit a distinct pre-edge peak and an intense main edge (Fig. 3A, curves c, d, e, and g); and the nondonor (ND) configuration with two broken donor H-bonds gives two distinct peaks (Fig. 3A, curve h) similar to the spectra of gas-phase water and the species at the surface of liquid water (20).

The pre-edge intensities for the SD configurations are due to excitations to unoccupied anti-bonding orbitals that are localized along the internal O-H bond (Fig. 4). The localization is caused by breaking one donor H-bond while keeping the other intact, which entails s-p-rehybridization in the orbitals close to 535 eV (14). The resulting increase of p-character translates to an intensity increase in the pre-edge region of the XA spectra through the dipole selection rule (12). The orbital that contributes to the pre-edge intensity closely matches zone A, as shown in Fig. 4. This comparison gives our H-bond criterion its physical basis and supports our interpretation; the presence of another molecule inside zone A changes this orbital and thereby the corresponding spectral feature, which also explains the sensitivity of XAS (XRS) to H-bond distortions on the H-side.

Fig. 4.

Connecting electronic and geometric structures. The contour plot of a typical unoccupied orbital giving rise to the pre-edge intensity in the XA spectrum of a water molecule with a broken donor H-bond (solid line contour plot) is compared with our cutoff criterion for H-bonding (black line in Fig. 3B) (12). Broken donor H-bonds are characterized by nearest neighbor O's outside the shaded area (12). The orbital contour has been calculated for the central molecule in a cluster with 11 molecules (the corresponding XA spectrum is shown in Fig. 2, curve b, solid line) where the H-bond on one H-side is broken (SD configuration). Only 4 of the 11 molecules are displayed.

The occurrence of the pre-edge feature, furthermore, can be directly correlated with the energetics in the H-bond interaction. For the model clusters considered here, we find that distortions leading to enhanced intensity in the pre-edge correspond to one donor H-bond losing about 40 to 70% of the individual bond energy of bulk ice (12). The SD species that dominate the liquid thus have one strong and one weak or broken donating H-bond. Consistent with commonly used definitions of H-bonds in liquid water (and for simplicity), we use the term “broken” for the distorted donor H-bond of an SD configuration.

In Fig. 3A, only the calculated spectra of SD species in curves c to e and g exhibit all features observed in the experimental liquid water spectrum: a pre-edge peak close to 535 eV and a strong main edge close to 537 eV. This result suggests, consistent with our phenomenological findings, that asymmetric SD species with one strong and one broken donor H-bond dominate the liquid water structure. In a more detailed analysis, our experimental data were fitted with sums of calculated spectra. The results are summarized in Table 1. Different combinations of computed spectra corresponding to different choices of bond lengths and angles within each class give similar spectra. All yield a dominant contribution of SD species both at 25°C and 90°C. We note the close agreement with the previously described analysis based on model experimental spectra (12). For 90°C we find the best fits for a decreased (increased) amount of DD (SD) species as compared with room temperature, indicating the breaking of H-bonds with increasing temperature along with the decrease of fully coordinated molecules (21). In particular, it comes as no surprise that bulk-ice–like configurations are virtually absent in water close to the boiling point. Consistent with the discussion of Fig. 1, we find that DD configurations are replaced by SD configurations both when heating the liquid and when melting ice into liquid water. Computed spectra indicate that XAS is not sensitive to the breaking of H-bonds on the O-side (acceptor H-bond) (13). However, symmetry arguments imply that, with a dominant amount of molecules in SD configurations, most molecules must in addition have a broken acceptor H-bond. This entails that, within the H-bond definition given here [Fig. 3 and (12)], each molecule has on average 2.2 ± 0.5 (2.1 ± 0.5) H-bonds at 25°C (90°C).

Table 1.

Relative amounts of local configurations—double-donor (DD), single-donor (SD), and nondonor (ND) configurations—in liquid water according to different MD simulations at 25°C and 90°C. The values in the first column for each temperature follow from the experimental observations (EXP) in Fig. 1 (12) and from fitted spectra (FIT) such as the ones shown in Fig. 5. The errors are given by the uncertainties related to the experimental observations (Fig. 1) and the calculations (12).

Method
Type EXP + FIT SPC MCYL CPMD
25°C
DD Embedded Image 70 50 79
SD 80 ± 20 27 41 20
ND 5 ± 5 3 9 1
90°C
DD Embedded Image 56 39 63
SD Embedded Image 37 47 34
ND 5 ± 5 7 14 3

As mentioned above, XAS cannot distinguish between SD configurations with a broken bond due to elongation (Fig. 3A, curve c) or to strong bending (Fig. 3A, curve e). However, these configurations exhibit different nearest neighbor O-O and O-H distances. We compare our results with O-O and O-H RDFs from neutron diffraction data taken at ambient conditions (2, 22). Two XA spectra generated with different distributions of H-bond lengths and angles, but both with the same fraction of DD, SD, and ND species, are shown in Fig. 5A, curves a and b. Spectra of this type have been used to estimate the amounts of the different species (Table 1) that are consistent with the measured XA spectrum.

Fig. 5.

Calculated XA spectra and RDFs (solid lines) for three models—(a), (b), and (c)—are compared with (A) experimental XA spectra and RDFs (B and C) at 25°C. Spectra and RDFs have been generated simultaneously for the central molecule in clusters consisting of 11 molecules, for a total of 14 clusters (12). A sampling of the possible DD, SD, and ND configurations for the calculated molecule in the liquid is guaranteed by using a variety of nearest neighbor O-O distances and H-bond angles. The RDFs were calculated for the first coordination shell only. Experimental RDFs [open circles in (B) and (C)] were derived from neutron diffraction (2). The calculated XA spectra were broadened with a 1 eV Gaussian for better comparison with experimental XRS data [open circles in (A)] taken from Fig. 1, curve e. Models (a) and (b) use a balance of DD:SD:ND = 10:85:5. Model a (b) has more (less) angularly distorted H-bonds and less (more) elongated ones. Model (c) (solid line) represents the SPC simulation (Table 1) with a balance DD:SD:ND = 70:27:3, as determined with our criterion for H-bonding (12). For comparison, the dashed lines in B and C, curve c, show the O-O and O-H RDFs calculated from the full SPC simulation. There is good agreement up to 3.2 Å (O-O) and 2.4 Å (O-H), validating our approach to characterize the first shell in the liquid. Only model (a) or similar balances with a dominant fraction of SD species (Table 1) give good fits when simultaneously considering our experimental XRS data and the RDF curves.

The O-O (Fig. 5B) and O-H (Fig. 5C) RDFs are shown for the two model distributions corresponding to the XA spectra in Fig. 5A, curves a and b (12). The 11-molecule model configurations that can adequately describe the local structure probed by XAS are obviously not sufficient for a full description of RDFs beyond the first coordination shell [<3.4 to 3.5 Å in the O-O RDF (2, 4) and <2.5 to 2.6 Å in the O-H RDF (2)]. Molecules from the second shell are not included in the RDFs of models a to c in Fig. 5 and hence we should not expect agreement for larger distances. The “a” model has a larger portion of broken H-bonds as a result of bending, whereas the “b” model has more elongated bonds. Both models reproduce the experimental XA spectrum, but model a clearly appears to provide a better agreement with the RDFs within the first coordination shell. Model b lacks structural contributions in the region 2.9 to 3.3 (1.9 to 2.3) Å and overemphasizes the contributions around 3.5 and 2.5 Å. It is thus inconsistent with water density, as illustrated by the poor agreement with the diffraction data. Combining XAS and diffraction results therefore suggests that the H-bonds are predominantly broken by bending rather than by elongation.

Finally, our results can be compared with two classical MD simulations using the flexible SPC (23) and MCYL (24) pair potential energy models and to an ab initio Car-Parrinello MD (CPMD) simulation (25) [for details of the simulations, see (12)]. The sampling of the distribution of structures is similar to that in the XAS/XRS experiment, i.e., the instantaneous geometries at a sequence of time steps are analyzed. Using our criterion for H-bonding (12), we have determined the amounts of DD, SD, and ND configurations (Table 1). The MCYL potential yields the largest amount of SD, but all of the simulations deviate substantially from experimental results. Note that with our spectroscopically determined electronic structure criterion, SPC (CPMD) at 25°C still gives an average of 3.3 (3.6) H-bonds per molecule, similar to earlier theoretical studies reporting 3.5 H-bonds per molecule (19). For example, with the numbers in Table 1 for SPC at 25°C, we determine that (0.7 × 4) + (0.27 × 2) + (0.03 × 0) = 3.34 H-bonds per molecule. An XA spectrum generated based on the distribution of species obtained from our SPC MD simulation at 25°C is shown in Fig. 5A, curve c. Clearly the agreement with the experimental spectrum is very poor: no pre-edge intensity (SD species) and too much intensity around 540 eV (fully coordinated DD species). The corresponding O-O (O-H) RDFs (12) are depicted in Fig. 5B (C), curve c, as solid lines. We note that the SPC simulation, in comparison with experimental results, exaggerates contributions at r-values characteristic for tetrahedral/near tetrahedral configurations (in the maxima) in both the O-O and O-H RDFs. However, the O-O (O-H) RDFs from SPC and our model a are rather similar at 2.9 to 3.3 (1.9 to 2.3) Å. The poor agreement with the experimental XA spectrum is hence attributed to too few species with large H-bond angles (SD species) in the SPC simulation. Similarly, the other two MD simulations shown in Table 1 cannot create an XA spectrum that fits the data. The MD simulations are, however, consistent with our results in predicting the decrease (increase) of DD (SD) species with increasing temperature.

The large amount of SD species can also explain femtosecond-IR studies indicating two different O-H groups in liquid water: one “strongly” and one “weakly” H-bonded (26). The orientational relaxation dynamics were shown to be directly connected with the strength of the two H-bond groups, where the weak H-bonds relax much faster than the strong H-bonds. According to our results, the number of strong H-bonds in the liquid is substantially smaller than expected, which may seem in contradiction with the small heat of melting compared with the heat of sublimation for ice. However, quantum chemical calculations have shown that each bond in the proposed SD configurations is stronger than the average bond in four-fold coordination because of anticooperativity effects (27, 28). Thus, the large number of weakened/broken H-bonds in the liquid leads to only a small change in energy. A recently developed quantum chemical model (1, 27) that proposes the predominance in the liquid phase of two hydrogen–bonded water molecules in ring conformations is consistent with our results. Water is a dynamic liquid where H-bonds are continuously broken and reformed (29). The present result that water, on the probed subfemtosecond time scale, consists mainly of structures with two strong H-bonds, one donating and one accepting, nonetheless implies that most molecules are arranged in strongly H-bonded chains or rings embedded in a disordered cluster network connected mainly by weak H-bonds.

Supporting Online Material

www.sciencemag.org/cgi/content/full/1096205/DC1

Materials and Methods

SOM Text

Figs. S1 to S11

Tables S1 and S2

References and Notes

References and Notes

View Abstract

Navigate This Article