Slightly more water than the solubility limit (0.57% water volume) was added to benzene-d6 in a sealed NMR tube at 298 K, and the tube was then warmed at 343 K to make the homogeneous sample solution fully saturated with dissolved water. The typical
1H-NMR spectra before and after warming are shown in Fig.
1a,b, respectively, where the proton signals are assigned to bulk water and dissolved water.
1H-NMR chemical shift at a higher magnetic field indicates that dissolved water exists as single molecules
15,21. The solution was then cooled to 298 K. The
1H-NMR spectrum of the supersaturated solution contains an additional sharp signal at a significantly lower magnetic field than that of trace bulk water (Fig.
1c). The proton signal with a very large chemical shift of 5.25 ppm suggests the presence of water with a very strong hydrogen-bonding structure, which we call a “water cluster”, compared with that of bulk water
15. The additional sharp signal with the chemical shift of 5.25 ppm appears also through sonication of the same mixture of benzene-d6 with a small amount of water (Supplementary Fig.
S1b). We also measured the D
2O concentration-dependent chemical shift of water clusters in Fig.
S1c,d. These results validate that the observed peak is from a water entity. The dissolved water (single molecules) is stable and entropic favourable below the solubility limit. When the dissolved water becomes supersaturated, the single water molecules quickly and spontaneously form water clusters.
It should be noted that the instantaneous appearance of the
1H-NMR signal is not observed for
water mixtures of polar solvents, where the proton signal of dissolved water continuously shifts downfield with increasing water concentration (Supplementary Fig.
S2a). The proton signal assigned to the water cluster is also observed for benzene derivatives, such as toluene and xylene (chemical shifts of 5.19 and 5.17 ppm, respectively, Supplementary Table
S1), and even for other hydrophobic solvents, such as chloroform (Supplementary Fig.
S2b). The appearance of a new very sharp proton signal of water with a large chemical shift is specific to water supersaturated in hydrophobic solvents. The solvent effect to the
1H-NMR signals of water clusters at different solvent ratios is shown in Fig.
S2c. This suggests that there is only one type of water cluster structures to be generated in hydrophobic solvents such as benzene, toluene, xylene, chlorobenzene, dichlorobenzene, trichlorobenzene, cyclohexane, carbon tetrachloride and chloroform.
As a control experiment, FT-IR spectroscopy was performed for the same sample solutions (Fig.
2). The IR absorption peaks assigned to dissolved water (and bulk water) are overlapped in a broad peak (probably ascribed to a water cluster). For the IR measurements under ambient conditions, it is not surprising that the absorption peak of bulk water (and also the water cluster) broadens because of dynamic hydrogen bonding, which will be discussed later.
Diffusion-ordered and nuclear Overhauser effect (DOSY and NOESY) spectroscopy and spin-lattice
T1 and spin-spin
T2 relaxation time measurement was performed to characterize the water clusters. DOSY spectroscopy gave diffusion coefficients for the protons ascribed to dissolved and bulk water of 5.0 and 2.3 × 10
−9 m
2 s
−1 (Supplementary Fig.
S3), respectively, which agree well with previously reported values
22,23,24. However, the diffusion coefficient of the water cluster (0.5 × 10
−9 m
2 s
−1) is surprisingly low (about 1/10 and 1/5 that of dissolved and bulk water, respectively). This means that the water molecules in the cluster have very restricted mobility. We showed the picture of the NMR tube (Fig.
S1b inset) to show the good homogeneity of the sample. The good homogeneity indicating that the measured effect is not related to artefacts on the glass surface. The
T1 and
T2 measurement gave the correlation time
TC which characterized the interaction induced by molecular motions and indicated that the protons in water clusters had a longer correlation time with the nearest neighbour water molecules than that in bulk water (Supplementary Table
S2)
14.
NOESY spectroscopy (Supplementary Fig.
S4) shows proton exchange (a negative nuclear Overhauser effect) between the dissolved water and the water cluster, which means that both the dissolved water and the water cluster are solutes in benzene and coexist (are in equilibrium) with mutual interaction. In addition, the chemical exchange build-up curve by different mixing time indicates that the exchange rate between water clusters and dissolved water is faster than that of the bulk water and dissolved water.
The
1H-NMR signal assigned to the water cluster slowly decreases with time (Fig.
1d), although the signal is still detectable after 3 days even at 298 K (Supplementary Fig.
S5). The water cluster is in a metastable state, and the data measured within 15 min do not cause a significant difference (error <2%).
The temperature dependences of the NMR signal intensities assigned to dissolved water and the water cluster are shown in Supplementary Fig.
S3a. The concentration of the water cluster increases with decreasing temperature accompanied by a quantitative and complementary decrease in the dissolved water amount, and the water cluster concentration reversibly decreases by increasing the temperature.
To investigate the equilibrium between dissolved water and the water cluster, the reversibility, or quasi-thermodynamical stability, of water cluster formation was analysed by the classical can’t Hoff plot (Inset of Fig.
3a). The Δ
H, Δ
S, and Δ
G values for water cluster formation from dissolved water were determined from the straight line in the van Hoff plot to give the energy diagram (Fig.
3b). The large enthalpy gains for the formation of the water cluster from dissolved water (about 34 kJ/mol) could be the driving force for cluster formation, which can be ascribed to the formation of multiple hydrogen bonds. However, cluster formation is accompanied by a large entropy loss (about −110 J/mol K). This can be explained by cluster formation or the formation of a more ordered structure. Δ
G to form the water cluster from dissolved water, which is the sum of the large enthalpy gain and entropy loss, is slightly negative, and the water cluster is a thermodynamically metastable state.
The transformation from the water cluster to bulk water has a negative Δ
G value (Fig.
3b), and the water cluster is indeed finally converted to bulk water. An example of the decay of the water cluster with time is shown in Fig.
1d. During a separate series of
1H-NMR measurements, spinning of the sample tube in the instrument accelerates the transformation of the water cluster to bulk water, where second-order kinetics is dominant. The temperature dependency of the decay rate in the Arrhenius plots gives apparent activation energy of decay of 38 kJ/mol (Supplementary Fig.
S6). The large activation energy of cluster decay supports that the water cluster is a metastable state.
The size-specific small water cluster number
n has been discussed using the experimentally determined continuous shift of the IR frequency with
n and by ab initio calculations
6,12,28. Here, we calculated the
1H-NMR chemical shifts of the protons of water for different cluster number
n (Fig.
4a).
The experimentally observed
1H-NMR chemical shift of 5.25 ppm (Fig.
1c and Supplementary Table
S1) is closest to the predicted chemical shift of 5.36 ppm for the cubic octamer D2d structure in Fig.
4b. This assignment does not conflict with the previous experimental and theoretical studies of cubic octamer
4,5,31.
The calculation also suggests that the dipole moment of the cubic octamer is zero (Fig.
4a), meaning that the water cluster is a nonpolar entity and differs from “normal” polar water. This helps to explain the long life of the water cluster in hydrophobic solvents. The dipole moment calculation also supports that the observed cluster species is the octamer by excluding polar water clusters with
n = 5–7, 9, and 10. The suggested cage-like octameric configuration of the water cluster also explains the very low diffusion coefficient of water measured by DOSY spectroscopy, because the diffusivity of the water molecules in the highly ordered water cluster is significantly reduced.
Shields
et al.
4 calculated Δ
H for the formation of an octamer with a static stable hydrogen-bonded structure to be −241 kJ/mol, whereas our experimental value is Δ
H = −34 kJ/mol (Fig.
3b). The difference between the two values suggests that the water cluster formed in this study could be composed of dynamic rather than static hydrogen bonding. This agrees with the experimental
1H-NMR signal with a large chemical shift being observed as a sharp peak rather than split peaks owing to individual protons. Dynamic hydrogen bonding in the size-specific cluster (
n = 8) is one of the features of the thermodynamically metastable water cluster formed in hydrophobic solvents.
Conventional
1H-NMR is very effective to investigate the water clusters that easily form in hydrophobic solvents under ambient conditions because of its superior signal resolution to IR and the additional information provided by DOSY, NOESY and relaxation time analysis. A thermodynamic study was performed to characterize the “metastable” state of the water cluster and cluster formation from dissolved water. The dynamic hydrogen bonding in the size-specific cluster is also described, in addition to the nonpolar property of the water cluster. We concluded that water clusters are a cubic octamer, which does not conflict with the previous papers, such as the “square-ice” observed in the hydrophobic conditions recently shown by Algara-Siller, G.
et al.
32.
Water clusters have attracted significant interest in many biological and chemical systems, for example, in bioinspired materials and devices where water molecules are in contact with or incorporated in organic hydrophobic materials
33,34,35,36,37,38,39. Investigation of water in organic entities as a matrix by
1H-NMR is expected to assist the formation and stabilization of water clusters and reveal the unique properties of water clusters.