Supercooled water confined in a metal-organic framework

Within the so-called “no-man’s land” between about 150 and 235 K, crystallization of bulk water is inevitable. The glass-like freezing and a liquid-to-liquid transition of water, predicted to occur in this region, can be investigated by confining water in nanometer-sized pores. Here, we report the molecular dynamics of water within the pores of a metal-organic framework using dielectric spectroscopy. The detected temperature-dependent dynamics of supercooled water matches that of bulk water as reported outside the borders of the no-man’s land. In confinement, a different type of water is formed, nevertheless still undergoing a glass transition with considerable molecular cooperativity. Two different length scales seem to exist in water: a smaller one, of the order of 2 nm, being the cooperativity length scale governing glassy freezing, and a larger one (> 2 nm), characterizing the minimum size of the hydrogen-bonded network needed to create “real” water with its unique dynamic properties. When water is held in a porous material and cooled below its freezing point the typical molecular interactions leading to crystallisation can be inhibited leading to a supercooled liquid state. Here, the authors use dielectric spectroscopy to analyse a metal organic framework with a symmetrical nanometer sized porous structure and reveal the temperature dependent processes and how they compare with those of bulk water.

T he confinement of liquids in spaces of nanometer size has proven a useful tool for the investigation of the glass transition [1][2][3][4][5][6] . Such experiments can provide valuable information about the increasingly cooperative nature of molecular motion assumed to arise during the gradual transition from a liquid into a glass under cooling 1,6 . Moreover, confinement studies enable to investigate the glassy and supercooled states of liquids that cannot be easily supercooled in their bulk form, avoiding crystallization when encaged in pores holding a limited number of molecules only. One of the most prominent examples in this respect is water 7-13 , a liquid of obvious relevance in numerous fields. Its glassy freezing and a liquid-to-liquid transition were suggested to occur within the so-called "no-man's land" between about 150 and 235 K, where crystallization of its bulk form is inevitable [14][15][16][17] . However, these suggestions are highly controversial [18][19][20][21][22][23] . Consequently, there are numerous investigations trying to solve this issue by studying the molecular dynamics and glassy freezing of supercooled water in confinement. Many different host materials were employed, e.g., porous glasses and molecular sieves 8,[24][25][26][27][28][29] , clay 7,30 , cement-related gels 31 , graphite oxide 32 , and even white bread 30 , to name just a few (see Swenson et al. 12 and Cerveny et al. 13 for an overview). Unfortunately, many of these materials exhibit highly hydrophilic inner pore surfaces, a considerable pore-size distribution, and/or asymmetric pore geometries, which may hamper the straightforward interpretation of the results.
Different scenarios were proposed to explain the numerous, partly contradicting experimental findings for confined supercooled water, speaking both in favor or against the suggested liquid-to-liquid transition (see Cerveny et al. 13 , Capaccioli et al. 33 , and Swenson et al. 34 for examples). Obviously, despite the vast database on confined water, drawing unequivocal conclusions that would lead to a consensus in the community, is difficult. This certainly also is due to the fact that the results considerably depend on the chosen host material 12,13 . Among them, in Cerveny et al. 13 water confined in the porous amorphous silica MCM-41 (Mobil Composition of Matter No. 41), which contains very regular cylindrical pores, was suggested to represent an "ideal model system of confined water". In a very recent work 33 , relating the dynamics of confined water for various host systems to a secondary relaxation, the molecular water dynamics in MCM-41 also was considered to effectively represent the intrinsic bulk-water dynamics.
Metal-organic frameworks (MOFs) are an alternative material class of confinement hosts. In general, MOFs comprise threedimensional networks with significant porosity which are formed by metal ions or clusters, joined together by bridging organic linkers [35][36][37] . Their highly ordered crystal structure ensures welldefined pore sizes. MOFs have outstanding design versatility: For example, by using different organic linkers the pore sizes and the apertures between them can be varied. In two recent reports 6,38 , we have demonstrated that MOFs are well-suited host materials for confinement studies of the glass transition. These works concentrated on the cooperativity length scale of molecular motions in glycerol, a good glass former that, in contrast to water, can be easily supercooled even in its bulk form. Water confined in MOFs has been previously investigated numerous times, however, mostly not addressing the question of glassy freezing within the no-man's land. To the best of our knowledge, there are only three studies of MOF-confined water using dielectric spectroscopy that touch this topic [39][40][41] . Unfortunately, the characteristic times of water dynamics, detected there, were rather far off those reported for the assumed "ideal" behavior in MCM-41. This fact is probably related to the limited frequency and temperature ranges of these studies and/or to the rather small pore sizes of the used host systems, ranging between about 0.4 and 1.1 nm (compared to about 2 nm for the most common form of MCM-41).
The MOF MFU-4l-HCOO (where MFU stands for "Metal-Organic Framework Ulm-University") seems to be a reasonable alternative to MCM-41. Just as for MCM-41, its pores are highly regular and of comparable size. In addition, they are close to spherical, leading to a 3D confinement of the water molecules, in contrast to the highly asymmetric confinement in the channellike pores of MCM-41 4 . In Pelster et al. 42 it was demonstrated that the dimensionality of the pores can influence the properties of confined liquids and, thus, confinement in 3D pores seems preferable for meaningful conclusions concerning the bulk properties. Finally, MFU-4l-HCOO exhibits less hydrophilic inner pore surfaces than MCM-41, whose internal surfaces are covered with numerous hydroxyl groups 43 . All these benefits make this MOF an even better-suited host system for the investigation of supercooled water.
In the present work, we find that water confined in the 1.9 nm pores of MFU-4l-HCOO can be well supercooled and exhibits relaxation times close to those found in MCM-41. Thus, this MOF indeed represents a well-suited host for the investigation of the glass transition of confined water. Much clearer than for other confinement hosts, a pronounced non-Arrhenius temperature dependence of the confined-water dynamics is revealed above a crossover temperature T cr ≈ 175 K. The high data precision, the broad frequency range, the performed evaluation of real and imaginary part including all processes, and the ideal host system available to us allow for a clear proof of this non-Arrhenius behavior. This crossover indicates cooperative molecular motions of supercooled confined water with growing cooperativity length for T > T cr . Below T cr , the cooperativity length is restricted by the pore size of about 1.9 nm resulting in Arrhenius behavior. Our data provide clear hints at two different length scales associated with glassy cooperativity and an extended hydrogen-bonded network, characterizing the short-and meso-scale molecular order in supercooled bulk water. Such a scenario was only rarely considered previously and never such clear experimental indications were provided.
The MFU-4l-related MOFs comprise two types of pores of different size. MFU-4l-HCOO was carefully chosen to achieve a diameter of the bigger pores (1.86 nm) that comes close to the channel diameter of the cylindrical pores in the "ideal" confinement host MCM-41 (2.1 nm for its most common variant). The smaller pores have a diameter of about 1.15 nm. As revealed by Fig. 1, in contrast to MCM-41, the pores have a highly symmetric, three-dimensional shape. Substituting HCOO − for the Cl − ligands, pointing inside the pores, reduces the sizes of the apertures between adjacent pores, leading to more effective confinement. The pore diameters of MFU-4l-HCOO were derived from the atomic coordinates in its crystalline lattice as reported by Denysenko et al. 47 , based on X-ray measurements. As these types of MOFs represent rigid covalent networks, the pore sizes should be unaffected by the water filling.
Compared to MCM-41, revealing numerous hydroxyl groups in its inner surfaces 43 , MFU-4l-HCOO represents an almost completely hydrophobic network. The only significant interaction of its pore walls with water may arise from the eight negatively charged formate groups pointing into its pores. Compared to the Cl − ligands of MFU-4l, the more delocalized charge of the HCOO − ligand should further reduce these interactions. Thus, overall in MFU-4l-HCOO the pore-wall interactions of the confined water molecules should be significantly weaker than in other host systems like MCM-41.
To estimate the theoretical maximum filling grade of MFU-4l-HCOO with water molecules, simulations were performed (see Methods section) leading to about 720 H 2 O molecules per unit cell. The right part of Fig. 1 shows the filling with water molecules as resulting from the simulations.
Permittivity spectra. Figure 2 presents spectra of the real and imaginary parts of the dielectric permittivity (ε′ and ε″, respectively), as measured at various temperatures for water in MFU-4l-HCOO. As commonly found for confined supercooled liquids, several relaxational processes are detected, signified by a steplike decrease in ε′(ν) and peaks or shoulders in ε″(ν), which continuously shift to lower frequencies for decreasing temperature. There is no trace of any abrupt qualitative change of the spectra when temperature is decreased that would point to a crystallization of confined water. This is confirmed by differentialscanning-calorimetry measurements between 100 and 200 K, which reveal heat-flow curves without anomalies, making crystallization of any significant fractions of water in MFU-4l-HCOO unlikely (see Supplementary Fig. 1).
The found four spectral contributions are numbered 1-4 as indicated in Fig. 2. For comparison, we have also measured waterfree MFU-4l-HCOO with empty pores. As an example, the plusses in Fig. 2 show the results for 189 K, which are featureless and, in case of ε, close to the instrument resolution. This proves that the processes detected for the MFU-4l-HCOO sample with water indeed are related to the water filling. Process 1 shows up as a strong power-law increase of both, ε′ and ε″ at low frequencies, which is most pronounced for higher temperatures. It leads to extremely large values of ε′ at low frequencies (not completely shown in Fig. 2), typical for non-intrinsic effects caused by electrode polarization 48 . For the present samples, it most likely is caused by the unavoidable ionic charge transport within hydrous samples, which essentially becomes blocked at the electrodes, grain boundaries, and/or pore walls. The relaxational nature of process 2 is clearly revealed by the step in ε′ (ν) [e.g., at about 10 4 Hz for the 210 K curve; Fig. 2(a)]. Due to the overlap with process 1, the corresponding loss peak is only visible as a shoulder [ Fig. 2(b)]. In contrast, the relaxational process 3 is well  The lines are simultaneous fits of the real and imaginary parts employing Cole-Cole functions for the three observed relaxation processes and a power law for a phenomenological description of the strong low-frequency increase. This power law is defined by ε″ = a ν s−1 and ε′ = a tan(sπ/2) ν s−1 with an exponent parameter s < 1 and a denoting a prefactor. The numbers in the figure indicate the different processes (1-3 are indicated for the 210 K spectrum; four denotes the fastest process, detected at low temperatures only). For comparison, the plusses show results for the pure metal-organic framework without water ("empty") at 188 K. pronounced both in ε′(ν) and ε″(ν). Finally, in the loss spectra of Fig. 2(b), at low temperatures another fast process is evidenced by an excess contribution showing up at high-frequencies at the right flanks of the peaks associated with relaxation 3. Its loss peaks are only incompletely seen in the covered frequency range and the corresponding ε′(ν) step is too weak to be detectable. Process 2 clearly is much slower than the relaxation dynamics of supercooled water that is unaffected by pore-wall interactions. This can be concluded from a comparison with results on supercooled bulk water (however, only available at temperatures beyond the borders of the no-man's land), aqueous water solutions, and water confined in different host systems as MCM-41 13,33,49 . It may be ascribed to water molecules interacting with the pore walls (probably by forming hydrogen bonds with the HCOO groups), which considerably slows down the molecular motions as often found in confined systems. 1,4,25,26 However, other explanations also are possible, e.g., by a partial solvation of the formate ion or in terms of a non-intrinsic Maxwell-Wagner relaxation, arising from the heterogeneous nature of the guest-host system 28 . In any case, this process does not reflect any dynamics close to that of bulk water 13,33 and therefore will be not be treated here in further detail.
In contrast, the characteristic time scale of process 3, estimated by its loss-peak frequency [ Fig. 2b], lies well within the region often considered for the "true" water dynamics of supercooled water within the no-man's land 13,33 . For example, the peak frequencies are located at similar frequencies as for the "ideal" confinement host 13,33 MCM-41 with 2.1 nm pore size (cf., e.g., Fig. 3 of Sjöström et al. 28 ). As both relaxational features, the ε′ step and ε″ peak, are well resolved in the spectra, reliable information on the temperature-dependence of its relaxation time can be obtained by fitting the data as described below.
Relaxation 4 obviously is significantly faster than the time scale of the bulk-water dynamics supposed for the no-man's land 13,33 . Various explanations of this process are possible. For example, it may be ascribed to water within the smaller pores of about 1.2 nm diameter existing in MFU-4l-HCOO. There the stronger confinement could lead to a considerable acceleration of molecular dynamics because the development of intermolecular cooperativity (which slows down the dynamics) can be suppressed in small pores. Such behavior is often found in supercooled liquids and can be used to estimate the size of the cooperativity length scale in dependence of temperature (see, e.g., these references 1,6,42 ). Within the framework of a cooperativity-driven glass transition [50][51][52] , this would imply that, for the temperature range where this process is detected (about 136-180 K), the cooperativity length in bulk water always exceeds the pore size of about 1.2 nm. As an alternative, process 4 can be interpreted as a secondary relaxation associated to the main water relaxation given by process 3. The occurrence of secondary relaxations, also termed β relaxations, which are faster than the α relaxation that reflects the main molecular motion governing the glass transition, is a common finding for supercooled liquids [53][54][55] . To distinguish such secondary relaxations from those sometimes arising from molecule-specific degrees of freedom like side-chain motions in polymers, they are commonly termed Johari-Goldstein (JG) relaxation 53 . Indeed, very recently a rather fast water relaxation (with relaxation times around 10 −5 s at 140 K), found in various confinements, was proposed to be of JG type 34 . However, one should be aware that in several works the much slower main confined-water relaxation, corresponding to process 3 of the present work, was also suggested to reflect the JG relaxation of bulk water 13,33,56 . In Swenson et al. 34 this apparent discrepancy was resolved by proposing two different β relaxations for confined water. It is clear that, based on the present data that do not fully resolve this process, we cannot clarify its true origin. Especially, it is not possible to deduce reliable values for the relaxation times from the fits of the spectra, described below. While the present work does not aim at the clarification of the origin of this process, nevertheless it is worth to mention that the lowtemperature loss spectra of Fig. 2(b) closely resemble those of supercooled liquids with a β relaxation that is submerged under the dominating α peak and leads to a so-called excess wing 57 . This also includes aqueous solutions like H 2 O:LiCl 58 .
Relaxation time. Of the detected processes, obviously relaxation 3 is the most relevant one to be compared to water dynamics in pure water and in other confinement systems as MCM-41. To determine its most important parameter, the relaxation time τ, we have fitted our broadband spectra (for a discussion of the other parameters of relaxation 3, see supplementary discussion). To account for process 1 in the fits, we phenomenologically described it by power laws ε" = a ν s−1 and ε′ = a tan(sπ/2) ν s−1 with an exponent parameter s < 1 and a denoting a prefactor (the additional tangent prefactor for ε′ ensures compliance with the Kramers-Kronig relation). The three relaxational processes were described by a sum of Cole-Cole functions, an empirical formula often used to parameterize dielectric relaxation processes 59 . One should be aware that, depending on temperature, only part of these contributions had to be included in the overall fit function, limiting the total number of parameters. For example, at high temperatures, process 4 and at low temperatures, process 1 and 2 have shifted out of the frequency window. The lines in Fig. 2 are fits, simultaneously performed for the real and imaginary part of the permittivity. Obviously, a reasonable agreement of fits and experimental spectra is achieved in this way. Figure 3 shows an Arrhenius representation of the temperature-dependent relaxation times of process 3 as obtained from these fits (closed circles). For comparison, data on MCM-41 with 2.1 nm pore size 13,28 are included in the figure (diamonds) as well as results on bulk water at high temperatures 49,58 (triangles) and on bulk low-density water well below the noman's land (crosses) 60 . At high temperatures, both τ(1/T) traces come close to that of liquid and supercooled bulk water, shown by the triangles 49,58 . This supports the notion 13,33 that water in MCM-41 represents a reasonable model system of confined water, which obviously also is the case for hydrous MFU-4l-HCOO. A closer look, however, reveals that both data sets do not seem to provide a perfect low-temperature continuation of the high-temperature bulk results measured outside the no-man's land. This is consistent with (highly disputed) [18][19][20][21][22][23] ideas that the true bulk-water data would exhibit a very strong upturn just after entering the no-man's land under (super)cooling (dash-dotted line in Fig. 3) 61 and that the MCM-41-confined water dynamics essentially reflects the β relaxation of water 13,34 .
While the τ(1/T) data for the MFU-4l-HCOO host at T > 175 K (1000/T < 5.7 K −1 ) are of similar order of magnitude as the corresponding data in MCM-41, they exhibit a remarkable difference: In contrast to MCM-41, they reveal clear deviations from simple thermally activated behavior, τ = τ 0 exp [E/(k B T)] (where E is an energy barrier and τ 0 an inverse attempt frequency). Instead they can be reasonably described by the Vogel-Fulcher-Tammann (VFT) law (solid line in Fig. 3), written in its modified form as proposed by Angell 62 : Here T VF denotes the Vogel-Fulcher temperature, where τ diverges and the strength parameter D characterizes the deviation from Arrhenius behavior 62 . VFT behavior is commonly found for the α relaxation of supercooled liquids and regarded as typical for glassy freezing 59,62,63 It can be ascribed to increasingly cooperative molecular motions when approaching the glass transition [50][51][52]64 . Non-Arrhenius behavior of confined water at temperatures above about 180 K, based on dielectric spectroscopy, was also discussed by some earlier works 13,28,32,56 . However, either at high temperatures only Arrhenius behavior was detected in these works, or the VFT temperature dependence was based on rather few data point only and, to our knowledge never, revealed as clearly as in Fig. 3. The effective 3D confinement, reduced pore-wall interactions, the performed evaluation of real and imaginary part including all processes and/or the extension of our spectra up to 1 GHz, leading to better precision of the high-temperature values of τ, may explain this fact. In any case, the steeper high-temperature Arrhenius region, revealed by hydrous MCM-41 (diamonds in Fig. 3 at 1000/T < 5.7 K −1 ) 13,28 , leads to an unphysical attempt frequency, 1/(2πτ 0 ), which should be of the order of a typical phonon frequency but is many decades smaller (about 10 −27 s) 28 . This indicates that deviations from the Arrhenius law must also be present for this host material 28 . The well-defined VFT law for water in MFU-4l-HCOO allows to determine its glass temperature using the usual definition τ(T g ) ≈ 100 s, which leads to T g ≈ 159 K. Its strength parameter D ≈ 11.3, corresponding to a fragility parameter of m ≈ 68 65 , characterizes this water phase in confinement as moderately fragile within the strong/fragile classification scheme proposed by Angell 62 . The notion of a liquidlike non-Arrhenius behavior of water confined in MCM-41 is consistent with Tombari et al. 66 , where the authors point out that the found specific-heat anomaly at roughly 225 K, observed in MCM-41 with 1.8 nm pores, is not due to a glass transition. Notably, the τ(1/T) curves of both water in MCM-41 and in MFU-4l-HCOO exhibit a change from stronger to weaker temperature dependence below a crossover temperature T cr ≈ 175 K (1000/T > 5.7 K −1 ; Fig. 3). In both cases, at T < T cr Arrhenius behavior is found with essentially identical activation energies of ≈0.48 eV, as evidenced by the similar slope of their τ (1/T) curves in Fig. 3. In this region, however, water in MFU-4l-HCOO relaxes by about a factor of 10 slower than in MCM-41. Interestingly, at low temperatures our data rather well match those reported for the low-density liquid (LDL) phase of bulk water (also termed low-density amorphous ice), which were measured beyond the low-temperature border of the non-man's land 60 .

Discussion
For the change in the temperature dependence of τ(T), as observed in MCM-41 at about 180 K and as also found in several other confinement systems 12,13,32 , various explanations were proposed. This includes the interaction of the water molecules with the pore walls 33 or an artifact of the analysis 27 that can arise when the β and α relaxations merge 67 . Moreover, various explanations in terms of a β relaxation were discussed, either assuming a transition from α-like to β dynamics 12,13,32,34 or in terms of a change of the temperature dependence of the β relaxation alone 11,68 . The latter would correspond to the crossover of the β-relaxation time τ β (T) from Arrhenius to VFT behavior, as often seen in bulk supercooled liquids when crossing T g under heating 54,63,69 . This scenario thus would imply a glass transition of T g = T cr ≈ 175 K for water confined in MFU-4l-HCOO.
However, an alternative explanation of our τ(T) results arises when considering their qualitative similarity to findings for confined glassforming liquids 1,6,38,42 . The motivation of such confinement studies originally was the verification of increasing cooperativity length scales ξ arising when approaching T g . As mentioned above, such a scenario was suggested for the explanation of the common non-Arrhenius behavior of the α dynamics of bulk systems, assuming that the glass transition is due to an underlying phase transition [50][51][52] . Indeed, when overcoming obstacles as pore-wall interactions or pore sizes that are much smaller or larger than ξ, a dynamical crossover of the α relaxation time from bulklike VFT temperature dependence to Arrhenius behavior could be observed under cooling for various glassforming liquids 1,6,38,42 . It occurs when ξ exceeds the pore diameter and, thus, cannot grow any further. Using hosts with different pore diameters then can provide information about the temperature-dependent cooperativity length 6 .
Therefore it seems natural to explain our results in confined water in the same way: Process 3 is the α relaxation of confined supercooled water (however, differing significantly from bulk water as discussed below) with T g ≈ 159 K. With decreasing temperature, at T cr ≈ 175 K its cooperativity length has grown to about 1.9 nm and, without confinement, would exceed the pore size under further cooling. Below T cr , the detected relaxation still reflects cooperative α dynamics, but with temperature-independent cooperativity because ξ is limited by the pore diameter d. This leads to the observed  13 . In addition, relaxation times for liquid 58 and supercooled bulk water 49 are included at high temperatures (upright and inverted triangles, respectively) and for bulk low-density liquid (LDL) water at low temperatures (crosses) 60 . The dashdotted line is a critical law with T c = 228 K, proposed for bulk water 61 . The solid line is a fit of the high-temperature (T > 175 K, i.e., 1000/T < 5.7 K −1 ) data of process 3 with the Vogel-Fulcher-Tammann (VFT) function (Eq. 1; parameters: T VF = 122 K, D = 11.3, τ 0 = 6.2 × 10 −15 s). The dashed line at lower temperatures represents an Arrhenius law with its energy barrier of 0.48 eV fixed to the same value as for the low-temperature data of water in MCM-41 (dotted line). The approximate extent of the non-man's land is indicated by the double-arrow on top of the figure. The inset demonstrates that τ(T) of process 3 can be described by the VFT (above T cr ) and Arrhenius laws (below T cr ) using identical prefactors τ 0 , which can be read off at 1000/T = 0. COMMUNICATIONS PHYSICS | https://doi.org/10.1038/s42005-020-0363-x ARTICLE COMMUNICATIONS PHYSICS | (2020) 3:95 | https://doi.org/10.1038/s42005-020-0363-x | www.nature.com/commsphys Arrhenius behavior, in agreement with model considerations 70 based on the Adam-Gibbs theory 50 .
Interestingly, within this scenario the approximate agreement of τ(1/T) at low temperatures with the relaxation times reported for bulk LDL water (Fig. 3) are consistent with ideas that the latter do not reflect the α dynamics of supercooled bulk water with freely unfolded cooperativity. By Cerveny et al. 13 it was suggested that the methods used to produce LDL water 60 may not be able to create the same type of water as could be reached by slow supercooling if crystallization would not intervene. For hyperquenched LDL water, this is supported by theoretical computations revealing characteristic length scales of less than about 3 nm, consistent with the present results 71 . In some respect, the detected confined water dynamics at T < T cr reminds of a β relaxation, which is often assumed to have no or only minor cooperativity. This view is, e.g., advocated in Capaccioli et al. 33 , based on an estimate of the β-relaxation time of bulk water using the coupling model 72 , which leads to values consistent with those experimentally found for hydrous MCM-41 at low temperatures, T < T cr 33 .
The above rationalization of the crossover at T cr is supported by the fact that both the VFT and the Arrhenius regions of the τ (1/T) curve can be consistently described assuming the same attempt frequency, 1/(2πτ 0 ), as demonstrated in the inset of Fig. 3. Further support for this cooperativity-related explanation is provided by Cerveny et al. 32 where such a scenario for confined water was already previously considered. There a systematic variation of the crossover temperatures in dependence of the pore sizes of various confinement hosts was found, just as expected: Under cooling, for larger pores τ(1/T) should adhere to VFT behavior down to lower temperatures, because ξ can grow further before exceeding the pore diameter at T cr . One should be aware that the absolute values of ξ deduced from such studies can only be estimates. For example, NMR experiments on water in MCM-41 have revealed the formation of crystalline water fractions near the pore center 11,68 and pore-size dependent partial crystallization was also considered in other works [73][74][75] . This implies a reduced confinement volume. In addition, a layer of slowed-down molecules interacting with the pore walls, which should depend on the used host material, may also reduce the effective pore size. Interestingly, T cr is somewhat higher for MCM-41 than for MFU-4l-HCOO (see the crossing points of the solid and dashed lines in Fig. 3) despite the nominally larger pore size of MCM-41 (2.1 instead of 1.9 nm). This is consistent with the suggested stronger pore-wall interactions in MCM-41 and the partial water crystallization in its pores, which both should reduce the effective pore volume available for cooperative motions.
One should note that, most likely, the confined α dynamics at T > T cr does not reflect that of the bulk if it could be detected within the no man's land. This is evidenced by the mentioned one-decade offset of the (extrapolated) τ(1/T) traces of bulk and confined water showing up at about 250 K (Fig. 3). To explain this finding, one may assume that the α-relaxation time of bulk water exhibits an extreme upturn of τ(1/T), just after entering the noman's land, as shown by the dash-dotted line in Fig. 3. The latter represents a critical law with T c = 228 K as proposed by Speedy et al. 61 Such a scenario was considered, e.g., in two papers by Swenson et al. 34,56 , where the confined water dynamics, even at T > T cr , consequently was denoted as β relaxation. As discussed in these works, indeed in some respect the confined dynamics detected below T c has β-like characteristics. However, one should be aware that, within this scenario, the glass temperature of bulk water [usually defined by τ(T g ) ≈ 100 s] would be only slightly above T c = 228 K (this becomes obvious when extrapolating the critical law in Fig. 3 to 100 s). This would imply that nearly all the confinement data of Fig. 3 were collected below T g . However, as mentioned above, below the glass transition, τ(T) of a β relaxation should exhibit Arrhenius behavior instead of the VFT temperature dependence observed in the present work. Even the αrelaxation time only follows VFT temperature dependence at T > T g 63,76 , except for extremely slow cooling or if τ is deduced from aging experiments 77 . Thus, even within the scenario of a critical water temperature of 228 K, we suggest that the observed relaxation process should be denoted as the α relaxation of confined water. It is this cooperative molecular motion that governs the glassy freezing of this special supercooled water phase in confinement and there is no slower relaxation that could assume the role of the α relaxation in confinement.
Our results imply that the MOF MFU-4l-HCOO represents an ideal confinement host with superior properties as effective 3D pore geometry and highly hydrophobic pore walls. Water confined in this system in several respects behaves similar as in MCM-41, previously considered as ideal confinement system 13,33 . However, in marked contrast to the latter and much clearer than in any other host system, our measurements up to 1 GHz reveal significant deviations from simple thermally activated dynamics above a crossover temperature T cr ≈ 175 K, where τ(T) follows the VFT law. This evidences cooperative, α-like relaxation dynamics of confined water, whose cooperativity length increases under cooling in this temperature region. Below T cr , τ(T) still reflects cooperative dynamics but with temperature-independent cooperativity length ξ, limited by the pore size of 1.9 nm as also found in supercooled liquids like glycerol 1,6,38,42 . As discussed above, the observed confined-water dynamics, both above and below T cr , in some respects exhibits the characteristics of a secondary relaxation as pointed out in several previous confinement studies of water (e.g., in these references 13,33,34 ). However, it can be considered as the α process of confined water (in accord with the term α conf in Cerveny et al. 13 ) as it governs the cooperativitydriven glassy freezing of this water phase.
As revealed by Fig. 3, the deduced τ(T) data of water confined in MFU-4l-HCOO at first glance provide a reasonable continuation of the high-temperature bulk data into the no-man's land (however, with some offset) and closely match the relaxation times of bulk LDL water at low temperatures. Thus, it is tempting to simply claim that the present results provide information on the properties of hypothetical supercooled bulk water within the no-man's land, which would imply a liquid-liquid transition at 175 K. However, in fact the situation is not so straightforward and the detected confined-water relaxation represents the dynamics of a different kind of water, whose behavior is significantly modified due to the confinement. Interestingly, at least at T > T cr this modification does not arise from the usual restriction of the molecular cooperativity associated with the glassy freezing as found in other confined liquids. This is signified by the detected non-Arrhenius behavior evidencing unimpeded growth of ξ(T) under cooling. These findings thus point to the presence of two different length scales in water: The smaller one, also found in other supercooled liquids, is the cooperativity length scale ξ (< 1.9 nm for T > T cr ) governing glassy freezing of confined water in terms of a hidden "amorphous-order" phase transition underlying the glass transition of all supercooled liquids 51,64 . The longer one (significantly larger than 1.9 nm) represents the characteristic length scale of an extended network of hydrogen-bonded molecules, which is essential for the specific molecular dynamics of bulk water. This length scale may, e.g., be related to the "patches" of fourbonded molecules considered within the percolation theory of water by Stanley and Teixeira 78 . As such an extended network cannot form in nm-sized pores, markedly different but still cooperative relaxation dynamics is revealed for water in confinement. The glass temperature of the resulting confined water phase is 159 K, it behaves moderately fragile (m ≈ 68), and, for infinitely slow cooling rate, its cooperativity length should diverge close to the Vogel-Fulcher temperature T VF = 122 K. Similar considerations involving such a network were also made 9,34 to explain the absence of a bulklike α relaxation in confined water which now is confirmed by the experimental results of the present work. By Tombari et al. 66 it also was suggested that the restricted pore size affects the hydrogen bonding and the resulting arrangement of the water molecules, leading to properties that differ from those of bulk water. There the confined-water state in MCM-41 was described as distorted ice-like structure in equilibrium with the water melt. As shown in the present work, water in MFU-4l-HCOO at least exhibits significant liquidlike fractions with relaxation times considerably faster than in ice 79 .
To corroborate the scenario promoted in the previous paragraph, measurements in MOFs with different pore diameters should be performed. The use of pores that are sufficiently large to allow for the unfolding of the extended hydrogen-bonded network, thus leading to bulklike properties, may be limited by the increasing importance of crystallization 73,75 . However, pore sizes only somewhat smaller and larger than 2 nm should be able to reveal the temperature dependence of the cooperativity length ξ governing the glass transition of confined water.

Methods
Host material synthesis. MFU-4l was synthesized according to the procedure described in Denysenko et al. 44 . MFU-4l-HCOO, used as confinement host in the present work, was prepared by postsynthetically exchanging the chloride sideligands in MFU-4 with HCOO-anions, similar to the procedure described in Ettlinger et al. 46 . The synthesis and structural parameters of MFU-4l-HCOO have been reported in the Supporting Information file of Denysenko et al. 47 .
Water content calculation. The theoretical water content of MFU-4l-HCOO was calculated as follows: In a volume of 1000 Å 3 , a content of 34 H 2 O molecules leads to a density of ca. 1.0 g/cm 3 . The void volume of MFU-4l is about 23575 Å 3 (ca. 78.7%) which would correspond to about 800 H 2 O molecules per unit cell. For the structural models, the initial positions of the water molecules were obtained from sorption calculations (Sorption Module in the Materials Studio 2019 software) employing force-field (Dreiding) derived energy values and partial charges (Qeq). In a first run, an empty box was filled with H 2 O molecules and the final fugacity was varied until a water density of about 1.0 g/cm 3 was reached. The same fugacity parameter (20,000 kPa) was then employed to calculate the position and numbers of water molecules in MFU-4l-HCOO. (NB: The very high fugacity required for this is due to the Dreiding force-field parameter sets who have not been optimized for condensed phases). These simulations led to a final value of 720 H 2 O molecules per unit cell.
Sample preparation. The host material was first degassed for 24 h at 150°C in vacuum. The loading with water was directly done within the dielectric measurement setup by exposing the degassed samples to a water-vapor atmosphere at room temperature. After 24 h the permittivity stopped increasing, indicating full saturation of the pores, whereupon the temperature-dependent dielectric measurements were started. The permittivity is a well-suited quantity to monitor the water content because water has a high dipolar moment and also introduces considerable ionic conduction, thus leading to strong contributions in both the real and imaginary part of the permittivity.
Dielectric measurements. For the dielectric measurements, the powder samples were filled into parallel-plate capacitors with a plate distance of 70 µm. Using powder avoids possible pressure-induced deterioration of the host material that could arise when pressing pellets. Only slight pressure was applied to the capacitor plates to enhance the packing density. Nevertheless, the detected absolute values of the dielectric permittivity may be somewhat reduced due to a packing density less than one. The dielectric measurements of the complex permittivity at frequencies from 1 mHz up to 1 GHz were performed using two experimental techniques: At frequencies up to about 3 MHz, a frequency-response analyzer (Novocontrol Alpha) was employed. High-frequency measurements up to about 1 GHz were carried out using a coaxial reflectometric setup 80 with an impedance analyzer (Keysight E4991A). For cooling, a closed-cycle refrigerator was used.

Data availability
The data that support the findings of this study are available from the corresponding author upon reasonable request. Received: 7 January 2020; Accepted: 30 April 2020;