NMR Studies on the Temperature-Dependent Dynamics of Confined Water

We use $^2$H NMR to study the rotational motion of supercooled water in silica pores of various diameters, specifically, in the MCM-41 materials C10, C12, and C14. Combination of spin-lattice relaxation, line-shape, and stimulated-echo analyses allows us to determine correlation times in very broad time and temperature ranges. For the studied pore diameters, 2.1-2.9 nm, we find two crossovers in the temperature-dependent correlation times of liquid water upon cooling. At 220-230 K, a first kink in the temperature dependence is accompanied by a solidification of a fraction of the confined water, implying that the observed crossover is due to a change from bulk-like to interface-dominated water dynamics, rather than to a liquid-liquid phase transition. Moreover, the results provide evidence that $\alpha$ process-like dynamics is probed above the crossover temperature, whereas $\beta$ process-like dynamics is observed below. At 180-190 K, we find a second change of the temperature dependence, which resembles that reported for the $\beta$ process of supercooled liquids during the glass transition, suggesting a value of $T_g\!\approx\!185$ K for interface-affected liquid water. In the high-temperature range, $T\!>\!225$ K, the temperature dependence of water reorientation is weaker in the smaller C10 pores than in the larger C12 and C14 pores, where it is more bulk-like, indicating a significant effect of the silica confinement on the $\alpha$ process of water in the former 2.1 nm confinement. By contrast, the temperature dependence of water reorientation is largely independent of the confinement size and described by an Arrhenius law with an activation energy of $E_a\!\approx\!0.5\ $eV in the low-temperature range, $T\!<\!180 $K, revealing that the confinement size plays a minor role for the $\beta$ process of water.


Introduction
Water exhibits many intriguing properties as a bulk and a confined liquid, which are of enormous relevance in nature and technology. 1,2 Bulk water shows several well-known anomalies, including the important density maximum. It was proposed that these anomalies originate from the existence of a second critical point in the deeply supercooled regime, which terminates a liquid-liquid (LL) phase transition between low-density and high-density water phases. 3,4 Experimental validation of the conjectured scenario is, however, not straightforward because bulk water crystallizes in the relevant temperature range, the no-man's land: 150-235 K. The existence of the LL phase transition thus still remains a subject of controversial scientific discussions. [5][6][7][8] The properties of confined water can depend on the features of the confining matrix and, hence, deviate from the bulk behavior.
It is well known that the melting temperature is reduced until regular crystallization is suppressed when decreasing the size of nanoscopic confinements. [9][10][11][12][13] Such suppression of crystallization allows one to keep confined water in the liquid state in the no-man's land and, hence, to explore the possible existence of a LL phase transition. Yet, it remains a crucial question to what extent the nature of water in severe confinements reflects that of water in the bulk liquid. When tackling this issue, it is important to consider that the properties of confined waters can vary across the confining geometry. 12,14,15 In particular, one may expect that confinement has stronger effects on the behavior of interfacial water near the matrix surface than on that of internal water in the confinement center.
MCM-41 compounds are ideal matrices to confine water on nanoscales because these silica materials exhibit nanopores with defined and tunable diameters. 16 Suppression of water crystallization was reported for MCM-41 with pore diameters d p r 2.1 nm. [10][11][12] Previous studies, which exploited this effect to characterize the properties of liquid water in the no-man's land, arrived at different conclusions about the existence of a LL phase transition. On the one hand, neutron scattering (NS) works found a sharp kink in temperature-dependent correlation times of water dynamics at B225 K, which was interpreted in terms of a fragile-to-strong transition, related to a LL phase transition. 10,17 On the other hand, dielectric spectroscopy (DS) studies did not observe such transitions, but rather a gentle crossover in the range of 180-190 K, which was attributed to a subtle interplay of structural a and local b relaxations during a glass transition. [18][19][20] Our recent 2 H NMR approach rationalized the observation of two crossovers in temperature-dependent correlation times at 220-230 K and 180-190 K, respectively, based on a two-step solidification scenario for water. 21 It was proposed that a fraction of water becomes solid at the higher of the two temperatures, leading to a change from bulk-like to interface-dominated dynamics for the other fraction, which stays liquid until it undergoes a confinementaffected glass transition at the lower of the two temperatures, causing another variation of the dynamical behavior.
Here, we investigate the role of the cavity size on water behaviors, in particular, on the proposed relation between crossovers in temperature-dependent correlation times and solidification events. For this purpose, the 2 H NMR studies on D 2 O in MCM-41 C10, 21 featuring 2.1 nm confinements, are extended to D 2 O in MCM-41 materials C12 and C14 with pore diameters up to 2.9 nm. A combination of spin-lattice relaxation (SLR), lineshape, and stimulated-echo (STE) analyses enables insights into both rates and mechanisms for water reorientation in broad dynamic and temperature ranges, as was also exploited in 2 H NMR work on water in zeolites. 22 In addition, we perform 1 H SLR measurements for H 2 O in MCM-41 C12.

Theoretical background
While our 1 H NMR studies are sensitive to the dipolar interaction between the nuclear magnetic moments of the protons, our 2 H NMR approaches probe the quadrupolar interaction between the nuclear quadrupole moment of the deuterons and an electric field gradient at the nuclear site, as resulting from anisotropic charge distributions within chemical bonds. The associated quadrupolar frequency is given by Here, the anisotropy parameter d and the asymmetry parameter Z describe the shape of the electric field gradient tensor and the angles y and f specify the orientation of this tensor and, hence, of the water molecule, with respect to the external magnetic field B 0 . The AE signs correspond to the two allowed transitions between the three Zeeman levels of the deuteron. 23 For D 2 O, the anisotropy parameter amounts to d = 2p Â 161 kHz, corresponding to a quadrupole coupling constant of 2p Â 215 kHz. Moreover, the asymmetry parameter is small, Z = 0.10, 21,24 so that o Q is approximately proportional to the second order Legendre polynomial P 2 (cos y), where y is the angle between O-D bond and B 0 field. Consequently, the fluctuations of the quadrupolar frequency o Q provide access to the rotational correlation function of the water molecules. Throughout this contribution, pointed brackets denote the ensemble average. 2 H SLR analysis yields information about the spectral density J 2 (o), which is related to the correlation function F 2 (t) via Fourier transformation. For isotropic reorientation, the relaxation time T 1 and the spectral density J 2 are linked according to 25 1 with a Larmor frequency o 0 of 2p Â 46.1 MHz in our case. When the correlation function of the observed molecular dynamics is a single exponential, F 2 (t) = exp(Àt/t), the spectral density takes the Debye form leading to a T 1 minimum for a correlation time t E 1/o 0 E 1 ns. Thus, upon cooling, the value of T 1 decreases when o 0 t { 1 at high temperatures, while it increases when o 0 t c 1 at low temperatures.
For confined water, the dynamical behavior is more complex and the correlation function is not a single exponential, but distributions of correlation times G(log t) exist. The Cole-Cole (CC) and Cole-Davidson (CD) spectral densities proved useful to describe such complex molecular dynamics: 26 They are characterized by the time constants t CC and t CD and by the width parameters b CC and b CD . The CC and CD spectral densities correspond to symmetric and asymmetric distributions of correlation times G(log t), respectively. 27 When observing the 2 H NMR line shape of confined liquids, a crossover from a broad Pake spectrum S P (o) to a narrow Lorentzian spectrum S L (o) occurs when isotropic molecular reorientation becomes faster than the time scale of the experiment, t E 1/d E 1 ms, upon increasing the temperature. 27 In the solid-echo experiments of the present approach, the transition between the broad and narrow spectra is accompanied by a reduction R(log t) of the signal intensity. Such reduced echo intensity for the intermediate motional regime, t E 1 ms, results because molecular dynamics during the dephasing and rephasing periods of the echo experiment interferes with echo formation. 24,28,29 For broad or bimodal distributions of correlation times, fast (t { 1/d) and slow (t c 1/d) molecules coexist in a certain temperature range. Then, the observed spectra S(o) are approximately described by weighted superpositions of Pake and Gðlog tÞ½1 À Rðlog tÞdðlog tÞ (7) Here, we take into account that the outcome of the experiment is determined by an effective distribution of correlation times, G(log t)[1 À R(log t)] due to attenuated echo signals of molecules with t E 1/d E 1 ms. Moreover, we use normalized line-shape components S L (o) and S P (o). 2 H STE measurements provide access to slow water reorientation, 10 À5 s o t o 10 0 s. In the STE approach, three pulses divide the experimental time into two short evolution times t p { t, which are separated by a mixing time t m E t. Using appropriate pulse lengths and pulse phases, it is possible to measure the rotational correlation functions 23,31-33 As the angular resolution of STE experiments increases with an extension of the evolution times, measurements of F cc 2 (t m ) and F ss 2 (t m ) for various values of t p provide access to the motional mechanism. 31,32,34 The angles of elementary rotational jumps can be determined when the time constant of STE decays is analyzed as a function of t p . While the decay time is virtually independent of the value of t p for large-angle jumps, e.g., tetrahedral jumps, it strongly decreases with increasing evolution time for small-angle jumps, e.g., rotational diffusion. The overall geometry of the reorientation is available from the evolution-time dependence of the residual correlation at long mixing times, 34 e.g., from F cc N (t p ) F cc 2 (t m ct, t p ). In particular, strongly anisotropic reorientation does not cause a complete loss of correlation, but leaves substantial residual correlation. The correlation function F 2 (t m ) is obtained from F ss 2 (t m ) when the evolution time is short and the anisotropy parameter, as in our case, is negligible so that sin(o Q t p ) E o Q t p p P 2 (cos y), see eqn (2).
For the analysis of experimental data, it is important to consider that, in addition to molecular reorientation, spin relaxation damps the STE signal. Therefore, we fit the normalized decays of F cc 2 (t m ) and F ss 2 (t m ) to the function Thus, we use a stretched exponential to describe the signal decay due to water dynamics and introduce F N to consider a residual correlation resulting, e.g., in the case of anisotropic reorientation. Moreover, we employ F(t m ) to take into account SLR damping. While the SLR damping of F cc 2 (t m ) can be determined in independent SLR measurements, allowing us to fix F(t m ) in the STE analysis, such independent determination of SLR effects is not possible when fitting F ss 2 (t m ). 35 Therefore, most of our studies involve the former correlation function.
Owing to broadly diversified dynamics of confined liquids, we use the mean logarithmic correlation time t m to characterize the time scale of water reorientation. In the SLR analysis, t m = t CC for the symmetric CC distribution, while this average value is obtained from hln ti ln t m = ln t CD + c(b CD ) + Eu (11) for the asymmetric CD distribution. 36 Here, Eu E 0.58 is Euler's constant and c(x) denotes the Digamma function. In the STE analysis, the mean logarithmic correlation time can be calculated according to 36 3 Experimental section

Sample preparation
The synthesis of the MCM-41 materials was performed according to the standard protocol of Grün et al., 37 which was adapted by Grünberg and coworkers. 38 Tetraethoxy silane, C 12 TAB (dodecyltrimethyl-ammonium-bromide) or C 14 TAB (tetradecyl-trimethylammonium-bromide), distilled water, and 25% NH 3 aq. were used in a ratio of 1 : 0.1349 : 138.8 : 2.8. The template was removed by calcination of the as-synthesized materials at 650 1C for 16 hours resulting in pure mesoporous silica materials. Tetraethoxy silane was acquired from Acros Organics. C x TABs were purchased from ABCR. All chemicals were used as received without further purification. Following the nomenclature in the literature, we denote the prepared mesoporous silica as C12 and C14, specifying the C x TAB chain length. In addition, we use mesoporous silica C10, which was synthesized by Kittaka and coworkers. 9,10,18,39,40 The synthesized MCM-41 samples were characterized employing N 2 gas adsorption together with multi-point BET to determine the sample surface area S p and the BJH method to obtain the specific pore volume V p and the pore diameter d p . The properties of the used MCM-41 C12 and C14 materials are shown together with that of MCM-41 C10 9,10,18,39-41 in Table 1. Prior to use, all MCM-41 matrices were dried by heating to 180 1C for at least 24 hours.
Filling with H 2 O or D 2 O employed different routes. The C12 and C14 samples were filled by carefully dripping water onto the dried mesoporous silica materials. The used amounts of water were calculated to obtain fillings of (100 AE 10)% of the pore volumes. For C10, we prepared two samples exhibiting different amounts of excess water outside the pores. Soaking the mesoporous silica in water and subsequent dabbing with filter paper  39 2.14 -1096 resulted in a sample with a substantial amount of excess water. Afterwards, the material was air-dried, resulting in a sample with a negligible amount of excess water, see below. All samples were sealed in a NMR tube.  2 H NMR spectra were recorded with the solid-echo sequence 90 x 1 À D À 90 y 1 where the delay D was set to 20 ms. In 2 H STE experiments, we used the pulse sequences 90 x 1 À t p À 90 Àx 1 À t m À 90 x 1 À t p and 90 x 1 À t p À 45 y 1 À t m À 45 y 1 À t p to obtain F cc 2 (t m ,t p ) and F ss 2 (t m ,t p ), respectively. When studying the evolution-time dependence, a fourth 901 pulse was added after a delay of 15 ms to overcome the dead time of the receiver for small values of t p . Appropriate phase cycles were used to cancel out unwanted single-quantum and double-quantum coherences. 33 4 Results and discussion

Spin-lattice relaxation
First, we analyze the buildup of 2 H magnetization, M(t), to investigate water dynamics in mesoporous silica. Fig. 1 shows results for C10 and C12 samples at several characteristic temperatures. For the studied samples, the buildup of magnetization is monoexponential at ambient temperatures, while it can involve several steps at lower temperatures. To rationalize the relaxation scenario in the latter range, it is useful to recall the freezing behavior of water in the studied MCM-41 materials. While crystallization of excess water outside the pores occurs near 273 K, freezing of confined water inside the pores depends on the pore size. [9][10][11][12] In C12, the freezing temperature of water amounts to about 220 K, whereas regular crystallization of water was found to be largely suppressed in C10. Nevertheless, calorimetric signals were observed in the latter confinement at cryogenic temperatures and attributed to a formation of glassy water or distorted ice. 9,13,42,43 Exploiting this knowledge about the freezing behavior, comparison of buildup curves M(t) for various water contents and confinement sizes enables an assignment of the observed relaxation steps to specific water species. First, we inspect results for C10 samples with substantial and negligible fractions of excess water, see Section 3.1, in Fig. 1(a). When excess water exists at 236 K, a relaxation step at short times, which continues the high-temperature relaxation, is accompanied by a relaxation step at long times, which sets in slightly below 273 K, but is absent in samples without excess water, as obtained from a more detailed analysis, see below. Therefore, we attribute the short-time and long-time steps to confined water and excess water, respectively. As was reported in our previous work, 21 a further change of the relaxation behavior occurs near 225 K, where yet another relaxation starts to split off from the short-time process. To demonstrate the effect, we compare buildup curves for the C10 samples with substantial and negligible fractions of excess water at 167 K. We observe a third relaxation step at times intermediate between the short-time and long-time relaxations. While this medium-time step, like the short-time step, does not depend on the amount of excess water, the amplitude of the long-time step is again smaller for smaller amounts of excess water. In Fig. 1(b), we see that the relaxation behavior is qualitatively similar for the C12 sample, in particular, three relaxation steps are distinguishable at sufficiently low temperatures. We conclude that water inside the pores causes not only the shorttime relaxation, but also the medium-time relaxation, while ice outside the pores contributes the long-time relaxation.
A quantitative analysis of the 2 H SLR behavior enables further insights into the nature of the different water species. Therefore, we fit the buildup curves of the 2 H magnetization to Here, M N denotes the equilibrium magnetization and the index n refers to the different relaxation steps. For reasons to be specified, we use n = l (liquid) for the short-time step, n = s i (solid inside) for the medium-time step, and n = s o (solid outside) for the longtime step. To allow for the possibility of nonexponential 2 H SLR, we utilize stretched exponentials for the interpolation of the individual contributions. Employing the Gamma function G(x), the mean SLR times hT 1,n i are obtained from the fit parameters according to Fig. 2 shows results of the SLR analysis for D 2 O in C10, C12, and C14. In panel (a), we see for C10 that, when a substantial amount of excess water exists and freezes slightly below 273 K, monomodal SLR turns into bimodal SLR. Upon further cooling, hT 1,l i decreases, as anticipated for fast dynamics (o 0 t { 1) in liquid water inside the pores, while hT 1,s o i rises, as expected for slow dynamics (o 0 t c 1) in an ice phase outside the pores. A comparison with hT 1 i values for bulk ice (I h ) 44 supports the latter assignment. Near 225 K, hT 1,s o i becomes too long for a reliable determination in a reasonable amount of time, whereas hT 1,l i passes a minimum, indicating that confined water exhibits correlation times t E 1/o 0 E 1 ns. In the same temperature range, the medium-time process appears as a new phenomenon and, hence, a third water species with a distinguishable dynamical behavior emerges. The findings that hT 1,s i i ranges well between the other relaxation times and increases when reducing the temperature mean that the dynamics of this water species is intermediate between the slow dynamics of the regular ice and the fast dynamics of the confined liquid. For the C10 sample with a negligible amount of excess water, a splitting into bimodal SLR in the vicinity of 273 K does not occur, but the medium-time relaxation step again appears near 225 K. Thereby, the value of h T 1,s i i is independent of the amount of excess water, supporting our conjecture that this water species, like the liquid species, resides inside the pores. In panel (b), we see that, qualitatively, both confined water species exhibit comparable SLR behaviors in all studied mesoporous silica. Quantitatively, the hT 1,l i minimum is shifted to a higher temperature in C12 and C14 as compared to C10 and, on the high-temperature side of the minimum, hT 1,l i is slightly longer in C12 and C14 than in C10, implying that water dynamics is faster in the larger confinements, at least in the weakly supercooled regime.
Additional knowledge can be obtained from the shape of the relaxation steps, as described by the stretching parameters b n . We find that, if existent, the medium and slow relaxations are nonexponential. They are characterized by stretching parameters b s i E 0.6 and b s o E 0.6, essentially independent of temperature and sample. By contrast, for all pore sizes, the fast relaxation is exponential above B185 K and nonexponential below this temperature. In the latter range, the value of b l continuously decreases when the temperature is decreased, resulting in b l E 0.6-0.7 at 150-160 K. In Fig. 2(b), we use dark and light symbols for hT 1,l i to indicate this crossover from exponential to nonexponential 2 H SLR near 185 K.
The observations for b l can be rationalized when we consider that water dynamics strongly varies across the pore volume, as found in previous studies on confined water, 15,45-49 including our work on D 2 O in C10. 21 In general, in 2 H NMR, such distribution of correlation times t results in a distribution of relaxation times T 1 and, hence, in nonexponential SLR. However, it is crucial for this argument that the correlation time of a molecule does not change during the buildup of the magnetization, which usually occurs on a much longer time scale than the dynamics of the molecule itself. Thus, our finding that the fast relaxation step is exponential above B185 K indicates that the correlation times of the associated water fraction are time dependent, implying that these water molecules sample a substantial part of the pore volume on the milliseconds time scale of the buildup process so that the resulting exchange of t values averages over any distribution of T 1 times, reconstituting exponential relaxation, as expected for a liquid. Vice versa, the continuous development of nonexponentiality found for the fast relaxation step below B185 K implies that the molecular diffusion becomes too slow to complete the exploration of a relevant part of the pore volume on the time scale of h T 1,l i, i.e., the system becomes nonergodic. Likewise, the nonexponentiality of the medium-time and long-time steps (b s i ,b s o o 1) shows that the associated fractions of water molecules do not scout different local environments on the time scale of the SLR process and, thus, they are solid. With the same arguments, the existence of several relaxation steps at sufficiently low temperatures provides strong evidence that there are dynamically distinguishable water fractions that do not exchange molecules until the buildup of magnetization is complete. An absence of this exchange is unlikely between two liquid phases. Therefore, our SLR results indicate that liquid and solid water species coexist inside C10 pores below B225 K, at least down to B185 K, where the fast relaxation step becomes nonexponential, too.
As for the liquid water fraction, the observed crossover from exponential to nonexponential SLR resembles the situation for supercooled liquids undergoing a glass transition. 32,50 At the glass transition temperature T g , the structural relaxation ceases to ensure ergodicity, resulting in this change of SLR behavior. It is important to notice that, at such temperatures, structural arrangements of supercooled liquids are too slow to be effective for SLR so that the T 1 values are determined by faster b relaxations rather than by this a relaxation. 51 However, changes of b dynamics in response to the freezing of a dynamics usually lead to a kink in the temperature-dependent SLR times hT 1 (T)i at T g . 32,50 For the liquid water fraction in all studied mesoporous silica, such a kink of h T 1,l (T)i is observed at B185 K, supporting our conclusion that the structural reorganization of interface-affected water vanishes near 185 K in a way resembling a glass transition.
So far, the 2 H SLR results imply that, in the studied silica matrices, a fraction of confined water solidifies at 220-230 K, while another fraction of confined water remains liquid below these temperatures. The liquid confined water diffuses in a subvolume of the pores restricted by the silica walls and the solid confined water, at least down to ca. 185 K. The solid water forming at 220-230 K inside the pores (s i ) is different from and may not be mistaken with the frozen water outside the pores (s o ) since the latter gives rise to yet another SLR step at even longer times.
A straightforward relation between relaxation times T 1 and correlation times t exists for exponential SLR, see eqn (3). Therefore, we restrict the analysis to the short-time relaxation step at T 4 185 K where b l = 1 and, consequently, hT 1,l i T 1,l . In Fig. 2, we see for all studied samples that the minimum value of T 1,l is larger than expected for a Debye process, indicating that the correlation function F 2 (t) of the water reorientation is not a single exponential. 32 This result is not surprising as dynamical heterogeneities are a characteristic feature of viscous liquids, in particular of confined liquids showing diverse molecular mobilities in various pore regions. To consider the existence of distributions of correlation times G(log t), we utilize the CC and CD spectral densities. Then, the width parameters of these spectral densities, b CC and b CD , can be determined from the minimum values of T 1,l . 32 The results are compiled in Table 2. Assuming temperatureindependent width parameters b CC and b CD and inserting the corresponding spectral densities J CC (o) and J CD into eqn (3), we determine mean logarithmic correlation times t m .
In Fig. 3, we present the results for liquid confined water in C10, C12, and C14. In the inserted graph, mean logarithmic correlations times obtained from the CC and CD spectral densities are compared with each other for the example of C10. For temperatures in the vicinity of the T 1,l minimum, both spectral densities yield consistent results. Differences become significant above 240 K and below 210 K. Thus, in the temperature range of the proposed LL transition of water, the 2 H SLR results hardly depend on the choice of the spectral density. Since application of DS to H 2 O dynamics in MCM-41 pores revealed a CC shape, 18 we focus on the mean logarithmic correlation times t m resulting from this spectral density in the following. In the main graphs, we see a crossover in the temperature dependence of these time constants at T = 220-230 K, which is prominent for C12 and C14, but exists also for C10. 21 At higher temperatures, there are deviations from an Arrhenius (ARR) law, which are strong for the former larger pores and weak for the latter smaller pores. In this range, a Vogel-Fulcher-Tammann (VFT) behavior, may describe the data, as proposed in NS work. 52 At lower temperatures, the temperature dependence of t m is weaker and, at least for C10, consistent with DS data, 18 which were described by an ARR law. Interestingly, the crossover between VFT-like and ARR-like behaviors for liquid confined water and the emergence of solid confined water occur in the same temperature range 220-230 K. Therefore, the kink of t m (T) does not necessarily provide evidence for the existence of a LL transition, as will be discussed in more detail below.
In the following 2 H NMR studies, we will exploit the different 2 H SLR times of the various water species to single out the contribution from the liquid fraction and to suppress that from the solid fractions. In detail, we destroy the magnetization and  start the acquisition after delays hT 1,l i { t d { hT 1,s i i, hT 1,s o i, ensuring that the magnetizations form deuterons in the liquid and solid fractions have recovered to major and minor extents, respectively, when the data are recorded. Hence, such partially relaxed (PR) experiments allow us to ascertain the dynamics of liquid water, essentially unperturbed by contributions from solid water inside or outside the pores. In all experiments, the delay t d was chosen to maximize the ratio of the signal contributions of the liquid and solid water species. Before continuing our 2 H NMR approach, we analyze 1 H SLR to investigate H 2 O dynamics in C12. We find that, when excess water is absent, the buildup of magnetization is exponential at all studied temperatures in 1 H NMR, different from the situation in 2 H NMR. This difference is a consequence of the fact that spin diffusion, i.e., a transfer of magnetization due to flip-flop processes of spins, is faster for 1 H than for 2 H mainly due to large differences of the gyromagnetic ratios of these nuclei. Specifically, when solid and liquid fractions are intimately mixed in the pore volumes, spin diffusion is sufficiently fast to exchange magnetization between both species and, hence, to establish a common T 1 in 1 H NMR, while a common 2 H SLR behavior is not necessarily established, as is observed in our measurements. In Fig. 2, we see that 1 H T 1 exhibits a clear minimum at B240 K, indicating a correlation time of t E 1/o 0 = 0.44 ns at this temperature. Inspection of Fig. 3 reveals that this correlation time from 1 H SLR well agrees with that from 2 H SLR. Hence, differences between H 2 O and D 2 O dynamics in C12 are small.  { 1 ms) and isotropic, a broad Pake spectrum is found at sufficiently low temperatures, revealing that the reorientation dynamics of water is slow (t c 1 ms). In the intermediate temperature range, the spectra can be described as a weighted superposition of Lorentzian and Pake components, suggesting that fast (t { 1 ms) and slow (t c 1 ms) molecules coexist. Comparing the spectra for the different pore sizes, we observe that the crossover between the Lorentzian and Pake shapes occurs at lower temperatures for C10 than for C12, indicating that water dynamics is faster in the smaller than in the larger confinement in the vicinity of 200 K.

Line-shape analysis
On first glance, it may be tempting to assign the Lorentzian and Pake lines to the liquid and solid water species, respectively. However, one has to bear in mind that PR experiments are designed to suppress the latter contribution. In previous work, 21 it was confirmed that a Pake contribution resulting from solid water is eliminated in the PR spectrum of D 2 O in C10. Therefore, the assignment of the Lorentzian and Pake components to the liquid and solid fractions of confined water was ruled out. Rather, it was argued that the liquid water by itself gives rise to a superposition of narrow and broad lines. Specifically, a coexistence of Lorentzian and Pake contributions can result from not only bimodal dynamics, but also a broad distribution G(log t), see eqn (7), as found for the liquid water species in our SLR analysis. Then, fast (t { 1 ms) and slow (t c 1 ms) molecules from the distribution are at the origin of the narrow and broad lines, respectively, while molecules with correlation times t E 1 ms can be neglected because their signals are poorly refocused during the used solid-echo pulse sequence, as described by the reduction factor R(log t), see eqn (7). In this scenario, the weighting factor of the Lorentzian line, W(T), i.e., the contribution of this line-shape component to the total spectral intensity, is expected to continuously decrease upon cooling as the distribution G(log t) shifts to longer times and, hence, the fraction of fast molecules decreases. Such a scenario was also proposed for benzene, 53 isobutyric acid, 54 bipyridine, 55 and naphthalene 56 confined in mesoporous silica. By contrast, in 2 H NMR on D 2 O in zeolites, 22 the observation of a comparable coexistence of narrow and broad line-shape features was not attributed to a broad distribution G(log t), but to a specific motional mechanism, which involves tetrahedral jumps and p flips with a single correlation time t.
To ascertain the role of dynamical heterogeneities for the line shape, we calculate solid-echo spectra from random-walk simulations. The methodology of these calculations was described in previous works. 57,58 In Fig. 5(a) and (b), we show simulated spectra for tetrahedral jumps with single and distributed correlation times. In both cases, narrow and broad spectral patterns can coexist during the line-shape transition. However, such a coexistence is found in a narrower dynamic range and, hence, a narrower temperature range for single than for distributed time constants. When using the above findings to map correlation times onto experimental temperatures, we observe that the coexistence range from the simulations for a single correlation time does not agree with that found in the measurements, while the expected and observed Lorentzian-Pake superposition regimes between about 175 and 200 K are in reasonable agreement for a Gaussian distribution G(log t) with s = 2.1, which resembles the CC distribution from the above SLR data. Thus, line-shape analysis provides evidence that confined water exhibits pronounced dynamical heterogeneities. On the other hand, this approach does not allow us to distinguish between tetrahedral jumps and isotropic reorientation since both mechanisms lead to comparable spectra when the same distribution G(log t) is used in the simulations, see Fig. 5

(b) and (c).
To further explore the conjecture that a broad distribution G(log t) governs the line shape of liquid confined water, we determine the weighting factor W(T) from Lorentzian fits of the central region of PR spectra. In Fig. 6, a continuous temperature evolution of W(T) is obvious for D 2 O in C10, C12, and C14. A similar increase of the weighting factor with increasing temperature is observed for C12 and C14, while the rise of W(T) occurs at lower temperatures and in a narrower range for C10. These observations are consistent with the outcome of our SLR analysis. More precisely, when we use the criterion W(T) = 0.5 to extract a correlation time t m = 1 ms, we find good agreement with the time constants from SLR analyses, see Fig. 3. Moreover, higher values of the width parameters b CC and b CD for C10 than for C12 and C14, see Table 2, imply a narrower distribution G(log t) for the former confinement, in harmony with a narrower temperature range of Lorentzian-Pake coexistence.
For a final check of the hypothesis that dynamical heterogeneities associated with the liquid fraction are responsible for the line shape, we calculate expectations for the Lorentzian contribution using the distribution G(log t) reported in a DS study on H 2 O in C10. 18 First, we neglect the reduction factor R(log t) associated with solid-echo experiments and obtain the Lorentzian contribution directly from WðTÞ ¼ Ð À6 À1 Gðlog tÞdðlog tÞ, see eqn (7). In Fig. 6, it is evident that the calculated and measured weighting factors roughly agree, but the transition range is wider for the former than for the latter. Therefore, it is necessary to consider the reduction factor of the experiment and to employ an effective distribution of correlation times for the calculations, explicitly, WðTÞ ¼ Ð À6 À1 Gðlog tÞ½1 À Rðlog tÞdðlog tÞ, see eqn (7). The exact shape of R(log t) depends on the mechanism for molecular dynamics, see Fig. 5(b) and (c). For isotropic random jumps, the full width at half maximum was found to amount to about two orders of magnitude. 24 Therefore, we use a logarithmic Gaussian distribution with a width parameter s = 0.85, which is centered at log(t/s) = À6. In Fig. 6, we see that this approach reasonably well describes the experimental data. Hence, we can conclude that the liquid fraction of confined water exhibits pronounced dynamical heterogeneities, leading to a superposition of Lorentzian and Pake lines in the NMR spectra at intermediate temperatures.

Stimulated-echo experiments
Finally, we use 2 H STE experiments to study slow water dynamics at low temperatures. Again, we perform PR measurements to focus on the behavior of liquid confined water. Fig. 7 shows F cc 2 (t m ) for D 2 O in C10 and C12 at various temperatures. For both samples, the decays are not exponential and shift to longer times when the temperature is decreased, in harmony with water dynamics that is governed by strong heterogeneity and slows down upon cooling. At given temperatures T r 185 K, the loss of correlation is comparable for the C10 and C12 samples, while relaxation-time and line-shape analyses indicated that water dynamics is faster in smaller pores at higher temperatures. Furthermore, the STE decays exhibit a weak two-step character with a short-time decay due to molecular dynamics and a longtime decay due to spin relaxation. Specifically, good interpolations of the experimental data result when fitting F cc 2 (t m ) to Logarithmic Gaussian distributions G(log t) with a width parameter of s = 2.1 were used, which closely resemble the CC distribution for C10 from the present 2 H SLR analysis. Specifically, both distributions have the same second moment. The characteristic time constants of the reorientation processes and the corresponding temperatures, as obtained for C10, see Fig. 3, are indicated on the left hand side. More precisely, the jump correlation times t j characterizing the exponential distribution of waiting times between consecutive reorientation events are indicated. For isotropic random jumps, the jump correlation time t j can be identified with the rotational correlation time t 2 of F 2 , while t j = 4 3 t 2 for tetrahedral jumps, see Fig. 9. For the distributions, the indicated values are mean logarithmic time constants. In addition, the reduced solid-echo intensities 1 À R are specified for the various dynamical scenarios. Fig. 6 Weighting factor W(T) describing the contribution of the Lorentzian line to the total intensity of partially relaxed spectra for D 2 O in C10, C12, and C14. The solid lines are guides to the eye. The other lines are predictions for C10 calculated from distributions of correlation times reported in a dielectric spectroscopy study on H 2 O in C10. 18 The dashed and dotted results were obtained when neglecting and considering the reduction factor R(log t) inherent to solid-echo experiments, respectively. See text for details. eqn (10) and fixing the parameters of the SLR damping F(t m ) at the values obtained from the SLR analysis in Section 4.1.
For a characterization of the decay owing to water reorientation, we calculate mean logarithmic correlation times t m using the fit results for t K and b K together with eqn (12). The results are included in Fig. 3. It is evident that the time constants hardly depend on the pore size. For D 2 O in both C10 and C12, an ARR law with an activation energy of E a = 0.5 eV describes the STE data. Furthermore, we see that the temperature dependence resulting from the STE analysis at lower temperatures T r 185 K is somewhat weaker than that obtained from the SLR data at higher temperatures. In both temperature ranges, the present NMR data are fully consistent with prior DS results for H 2 O in C10, [18][19][20] which reported a mild change of the temperature dependence at 180-190 K.
To obtain insights into the mechanism for water reorientation at low temperatures T o 185 K, we exploit that the angular resolution of STE experiments is enhanced when the evolution time is extended. In Fig. 8, we show correlation functions F cc 2 (t m ) and F ss 2 (t m ) at 160 K for various lengths of the evolution time t p . It can be seen for D 2 O in C10 and C12 that the dependence of the correlation functions on the evolution time is weak. For a quantitative analysis, the STE decays are again interpolated with eqn (10). Fig. 9 shows the resulting time constants and residual correlations. In Fig. 9(a), we see that the time constants t cc K (t p ) and t ss K (t p ) are basically independent from the value of the evolution time and the type of the correlation function. This independence indicates that overall water reorientation results from large-angle rather than small-angle elementary jumps, as is confirmed by comparison with curves calculated for different motional models. 24 Specifically, the experimental findings are in reasonable agreement with expectations for jumps about the tetrahedral angle, while they are in contrast to simulations for isotropic jumps about an angle of 101. In passing, we note that the similarity of the F cc 2 (t m ) and F ss 2 (t m ) data and the weak evolution-time dependence of the results indicate that the choice of the experimental parameters is not critical when we use STE experiments to determine the temperature dependence of water dynamics, as was tacitly assumed in the above analysis.
As for the residual correlation, determination of reliable values is possible for F cc N , but not for F ss N since the long-time plateau of the correlation functions is affected by relaxation effects, which are known from independent experiments in the case of F cc 2 , but not in the case of F ss 2 . 35 In Fig. 9(b), residual correlations F cc N (t p ) E 0.2 AE 0.1 can be observed for the studied evolution times with a mild tendency to decrease when this time interval is extended. Thus, the fit results confirm the above conclusion from visual inspection that the probed low-temperature dynamics of liquid water does not cause complete decays of the correlation functions. Further information is available when we contrast the experimental data with expectations for isotropic reorientation, tetrahedral jumps, and 1801 jumps around the molecular symmetry axis. The observed residual correlation is significantly higher and lower than the expectations for isotropic In these fits, we exploit that the spin-lattice relaxation damping F(t m ) is known from independent analysis in Section 4.1. We obtain stretching parameters in the range b K = 0.25-0.39 for the studied samples and temperatures.  Parameters characterizing the dependence of the correlation functions F cc 2 (t m ,t p ) and F ss 2 (t m ,t p ) on the length of the evolution time: (a) time constants t cc K (t p ) and t ss K (t p ) and (b) residual correlation F cc N (t p ). The results were obtained by fitting stimulated-echo decays of D 2 O in C10 and C12 at T = 160 K. In panel (a), the experimental data t cc,ss K (t p )/t ss K (t p -0) are compared with simulation results for tetrahedral jumps (solid line) and isotropic 101 jumps (dashed line). In panel (b), the measured data are contrasted with calculated expectations for isotropic reorientation (solid line), 1801 jumps about the molecular symmetry axis (dotted line), and tetrahedral jumps. In the latter case, we distinguish between exact tetrahedral jumps (dashed line) and distorted tetrahedral jumps (dash-dotted line), where we randomly choose orientations on the edges of cones with semi-opening angles of 31 around the tetrahedral positions. reorientation and 1801 jumps, respectively. Hence, water reorientation does exhibit significant anisotropy at the studied temperatures, but it is not given by a simple two-site jump related to the molecular symmetry. Rather, the observations for F cc N (t p ) are consistent with a tetrahedral jump, in particular, when we allow for some distortions (AE31) of the molecular orientations, as expected in disordered hydrogen-bond networks. However, a unique determination of the motional geometry is not possible since similar behaviors F cc N (t p ) are produced by other anisotropic motions 31 and uncertainties in the determination of the plateau values result from imperfections in the suppression of the water species with the longer relaxation time T 1,s i in the partially relaxed measurements.

Conclusion
We employed a combination of 1 H and 2 H NMR methods to investigate H 2 O and D 2 O in MCM-41 C10, C12, and C14. For the studied matrices, we found that all confined water is in the liquid state above 220-230 K, whereas fractions of liquid water and solid water coexist inside the pores below this temperature range, as indicated by an observation of separate SLR steps. Most probably, the liquid and solid fractions can be identified with interfacial and internal waters, respectively.
For liquid confined water, the temperature-dependent rotational correlation times show two crossovers. A crossover at 220-230 K accompanies the appearance of solid confined water. While a VFT behavior roughly describes water reorientation in the confined liquid above the crossover, clear deviations exist below. A crossover at 180-190 K is indicated by a mild weakening of the temperature dependence, which is better seen in prior DS results 18-20 than in the present NMR data. The latter effect comes along with a change of the SLR behavior. While an exponentiality of SLR together with an observation of Lorentzian line shapes revealed that liquid water molecules show an isotropic reorientation associated with an exploration of a substantial part of the pore volume above the latter crossover, a nonexponentiality of SLR and an existence of finite residual correlation provided evidence that these molecules exhibit an anisotropic and localized reorientation below. These observations imply that the reorganization of the water network ceases at 180-190 K.
Based on these results, we propose the following scenario. Above 220-230 K, our NMR studies probe the a process of water molecules in all pore regions. In this temperature range, the temperature dependence becomes stronger and, hence, the bulk behavior is approached when increasing the pore radius, as expected. Two findings show that an observation of a kink in temperature-dependent correlation times of confined water at B225 K is not sufficient to conclude on the existence of a LL phase transition at this temperature. First, an appearance of solid water, most likely, freezing or vitrification 19 of water in the inner area of the pore, accompanies the alteration of the temperature dependence found for liquid water at 220-230 K, implying that the space available to the liquid fraction changes from 3D-like to 2D-like, i.e., to the outer area of the pore. Second, liquid water exhibits a prominent b process below 220-230 K, implying that a splitting of the b process from the a process contributes to changes in temperature-dependent correlation times in this range. Notwithstanding the foregoing, the present results do not allow us to rule out a LL phase transition for bulk water.
In the temperature range from B225 K to B185 K, the NMR results related to liquid water start to be governed by the b process, but these water molecules still explore a substantial part of the pore volume on a longer time scale, restoring ergodicity and, thereby, exponentiality of SLR. Thus, structural relaxation still occurs. While the strong confinement provided by the pore walls and water nanosolids may alter the nature of the a process, the scenario, in a sense, resembles that in bulk supercooled liquids exhibiting a-b splitting upon cooling. Compared to other liquids, 58-60 the b process of water is, however, related to molecular reorientation with much larger amplitude and, hence, it dominates the loss of orientational correlation, rendering an observation of the a process of water difficult. Because of this difficulty, there is still some room for discussions about the existence of the a process of confined and bulk waters below 220-230 K. 19 The observations for liquid water at 180-190 K resemble changes reported for the b process of supercooled liquids in response to a glass transition. In particular, the SLR of liquid water changes from exponential to nonexponential, implying that the a process ceases to restore ergodicity in this temperature range, as observed for bulk liquids at T g . 32,50 Also, the temperature dependence of the SLR times of the liquid fraction, hT 1,l i, becomes weaker. For bulk liquids, such effect results from changes of the b process owing to the freezing of the a process during vitrification. 32 Therefore, we propose that interfacial water exhibits a glass transition at 180-190 K. This conjecture is supported by the observation of calorimetric signals for confined water at such temperatures. 13,40 Due to the strong geometrical restrictions, it is, however, questionable whether the present and previous results for confined water are of any relevance for the glass transition of bulk water.
Below 180-190 K, the b process of water is associated with anisotropic and localized motion, yet it continues to involve reorientations about large angles, in agreement with findings for water dynamics at other interfaces. 24,[61][62][63] At these lower temperatures, the correlation times of water reorientation follow an ARR law with an activation energy of E a E 0.5 eV. They hardly depend on the pore size, at variance with our observations at higher temperatures, providing further evidence that distinct relaxation processes are probed in different temperature regimes.
There is a certain similarity between the present findings for water dynamics below B225 K and previous results for small mobile molecules close to large immobile molecules in binary mixtures, which were interpreted in terms of a type A glass transition, 64 as predicted within mode-coupling theory. 65 For such type A glass transition, a nonergodicity level continuously develops upon cooling, in contrast to a discontinuous jump of the nonergodicity parameter at the critical temperature T c for a type B glass transition. Thus, in future studies, the temperature evolution of the long-time plateau of the correlation function should be ascertained in more detail to distinguish between possible descriptions of water dynamics in the low-temperature regime.
Altogether, the present 2 H NMR studies provide evidence that, upon cooling, the dynamical behavior of liquid water in silica pores with diameters of 2-3 nm changes at 220-230 K and 180-190 K. These changes are closely related to a formation of crystalline or glassy water species in combination with an onset of a pronounced b relaxation, which exhibits a number of common features for water in various confinements and mixtures.