==== Front Nano Lett Nano Lett nl nalefd Nano Letters 1530-6984 1530-6992 American Chemical Society 37294683 10.1021/acs.nanolett.3c00235 Letter Mechanism of Water Intrusion into Flexible ZIF-8: Liquid Is Not Vapor https://orcid.org/0000-0002-9090-2408 Amayuelas Eder † https://orcid.org/0000-0002-3197-2780 Tortora Marco ‡ https://orcid.org/0000-0001-9649-1470 Bartolomé Luis † Littlefair Josh David § Paulo Gonçalo ‡ Le Donne Andrea § Trump Benjamin ∥ Yakovenko Andrey Andreevich ⊥ https://orcid.org/0000-0002-8912-9024 Chorążewski Mirosław ∇ https://orcid.org/0000-0003-2735-6982 Giacomello Alberto *‡ https://orcid.org/0000-0003-1220-5866 Zajdel Paweł *¶ https://orcid.org/0000-0002-3925-3799 Meloni Simone *§ https://orcid.org/0000-0001-6523-1780 Grosu Yaroslav *†∇ † Centre for Cooperative Research on Alternative Energies (CIC energiGUNE), Basque Research and Technology Alliance (BRTA), Alava Technology Park, Albert Einstein 48, 01510 Vitoria-Gasteiz, Spain ‡ Dipartimento di Ingegneria Meccanica e Aerospaziale, Sapienza Università di Roma, via Eudossiana 18, 00184 Rome, Italy § Dipartimento di Scienze Chimiche e Farmaceutiche (DipSCF), Università degli Studi di Ferrara (Unife), Via Luigi Borsari 46, I-44121 Ferrara, Italy ∥ NIST Center for Neutron Research, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA ⊥ X-Ray Science Division, Advanced Photon Source, Argonne National Laboratory, Argonne, Illinois 60439, USA ∇ Institute of Chemistry, University of Silesia, Szkolna 9, 40-006 Katowice, Poland ¶ Institute of Physics, University of Silesia, 75 Pulku Piechoty 1, 41-500 Chorzow, Poland * alberto.giacomello@uniroma1.it * pawel.zajdel@us.edu.pl * simone.meloni@unife.it * ygrosu@cicenergigune.com 09 06 2023 28 06 2023 23 12 54305436 18 01 2023 21 03 2023 © 2023 The Authors. Published by American Chemical Society 2023 The Authors https://creativecommons.org/licenses/by/4.0/ Permits the broadest form of re-use including for commercial purposes, provided that author attribution and integrity are maintained (https://creativecommons.org/licenses/by/4.0/). Zeolitic Imidazolate Frameworks (ZIF) find application in storage and dissipation of mechanical energy. Their distinctive properties linked to their (sub)nanometer size and hydrophobicity allow for water intrusion only under high hydrostatic pressure. Here we focus on the popular ZIF-8 material investigating the intrusion mechanism in its nanoscale cages, which is the key to its rational exploitation in target applications. In this work, we used a joint experimental/theoretical approach combining in operando synchrotron experiments during high-pressure intrusion experiments, molecular dynamics simulations, and stochastic models to reveal that water intrusion into ZIF-8 occurs by a cascade filling of connected cages rather than a condensation process as previously assumed. The reported results allowed us to establish structure/function relations in this prototypical microporous material, representing an important step to devise design rules to synthesize porous media. synchrotron radiation MOFs intrusion mechanism hydrophobic surfaces molecular dynamics H2020 European Research Council 10.13039/100010663 803213 H2020 Future and Emerging Technologies 10.13039/100010664 101017858 document-id-old-9nl3c00235 document-id-new-14nl3c00235 ccc-price This paper originally published ASAP on June 9, 2023. A change was made to the Acknowledgment, and a new version reposted on June 12, 2023. ==== Body pmcWetting and dewetting of microporous media, i.e., materials with pores diameter smaller than 2 nm, is of great importance for broad range of technological and natural systems.1,2 Among others, energy-absorbing and energy-storing systems have attracted considerable attention in recent years as they are expected to play a crucial role in a future carbon free society.3 Many of these systems transfer mechanical energy through pressurized intrusion–extrusion cycles of a liquid into/out of a porous material.4,5 The wetting and dewetting of a microporous material is a truly multiscale process, in which the emerging characteristics of the intrusion/extrusion cycle exhibit a nontrivial dependence on physio-chemical details at the molecular level: the type of material, its structure and morphology and its chemical composition.6,7 Thus, a microscopic understanding of the intrusion and extrusion processes is key to design novel, high-performance materials. In this context, metal–organic frameworks (MOFs) have arisen as a promising class of materials to achieve breakthrough performance in energy absorption during mechanical impact.8−11 The large inner surface area (up to 10,000 m2g–1), the highly tunable framework architecture, and the chemical composition of MOFs12 make these materials attractive candidates for energy-dissipation systems. Among the vast number of MOFs, only a handful have been studied for this application,13−16 owing to their good chemical, thermal, and mechanical stability, and the relatively inexpensive chemicals and synthesis protocol for their production; they predominantly belong to the hydrophobic microporous family of zeolitic imidazolate frameworks (ZIFs).17 Arguably, the most extensively studied microporous material to date is ZIF-8, which has become a reference for microporous materials in the area of water intrusion–extrusion and related applications.8,15,18−20 Traditionally, research in the field focused on the empirical evaluation of energy absorption performance of ZIF-8 and other microporous media. However, attention has recently turned toward the investigation of the microscopic mechanism of hydrophobic ZIF wetting.21 Deepening the understanding of hydrophobic wetting has been the main goal of several groups in the scientific community in recent years, exploring the effects of confining water in hydrophobic nanopores, the interactions of water with different hydrophobic surfaces, and the liquid and vapor interfaces in nanoconfinement under hydrostatic pressure.22−27 In a recent work, Sun et al.,28 proposed the intrusion of ZIF-8 to be governed by the condensation of vapor present in the nanoscale cages, postulating that the kinetics of the intrusion process is determined by the intrinsic length (single nanometer) and time scales (nanoseconds) necessary for critical water clusters to nucleate inside individual cages. Understanding the key rate-limiting step of intrusion is of paramount importance to design improved materials for energy dissipation and the key contribution of Sun et al.28 paved the way to such investigation. Indeed, the hysteresis, i.e., the difference of thermodynamic conditions between intrusion and extrusion of a nanoporous material, determines how much energy per intrusion-extrusion cycle is dissipated. Hysteresis, in turn, depends on how quickly the system reaches the equilibrium state corresponding to the given thermodynamic conditions, whose microscopic determinants are investigated here. In this work, combining intrusion experiments, in operando synchrotron radiation, atomistic simulations, classical continuum theories, and stochastic models, we report some important insights into the multiscale mechanism of ZIF-8 MOF wetting by H2O. We reveal that the intrusion of water in ZIF-8 proceeds by the cascade penetration of water in connected cages, forming coherent domains of ZIF-8 cages with alike structural characteristics. This mechanism is fundamentally different from the previously proposed condensation-driven intrusion in individual ZIF-8 cages. We show that the proposed cascade intrusion mechanism crucially depends on the hydrogen bonds which form across neighboring cages, which would instead be irrelevant in the condensation scenario. On a larger scale, water intrusion progresses by the formation and growth of coherent domains of wet cages, which progressively and cooperatively advance through the MOF crystallite. Our conclusions support new design criteria for MOFs for energy applications, hinging on pore connectivity rather than on individual cage properties. First, we focus on structural changes in ZIF-8 (flexible sodalite, SOD, topology; Figure 1c)29 during water intrusion at 30 °C by in operando synchrotron radiation and intrusion porosimetry. Figure 1a shows the evolution of the unit cell parameter a during the experiment, namely the expansion of a during H2O intrusion (Figure S1). We have already reported the negative compressibility of ZIF-8 upon intrusion in a recent article.21 The analysis of the shape of the diffraction pattern peaks allowed us to identify the strain (Figure 1b) produced in ZIF-8 (Basolite Z1200) under hydrostatic pressure scans (further experimental details in the Supporting Information). Additionally, we were able to identify the coherent crystalline domain size, i.e., the largest domain of crystallites with the same structural characteristics. In fact, wet cavities are larger than dry ones (Figure 1d), and if, during intrusion, domains of fully wet and completely dry cages are formed (Figure 1e), these can be detected by synchrotron radiation and the size of the largest one determined. Figure 1 a) Evolution of lattice parameter during H2O intrusion at 30 °C. b) Evolution of apparent strain and coherent crystalline domain size versus pressure for ZIF-8 + water system at 30 °C. c) SOD topology of ZIF-8 showing the interconnected cavities. d) Schematic simplification of the wetting of external cavities in a ZIF-8 crystal, becoming larger than the dry ones (blue color stands for water. e) Wetting process of an extended ZIF-8 structure, with different domain sizes: the wet ones larger (blue) and the dry ones smaller (white). We remark that the cartoon shown in this panel is purely illustrative: synchrotrons do not allow us to identify the geometry of wet/dry domains, whether cubic-like, spherical-like or more complex, nor if a single or multiple wet/dry domains are formed during intrusion. Both the strain and the size of coherent domains present a sharp variation in correspondence of the intrusion process (compare Figure 1a,b, and Figure S1). Considering that wet ZIF-8 cages undergo an expansion relative to dry ones,21 intrusion introduces some strain in partly filled ZIF-8 crystallites due to the size mismatch between wet and dry cages (Figure 1e). In other words, expansion of wet cages during intrusion breaks the local symmetry introducing a strain. This strain continuously increases with the number of wet cages up to a maximum. We hypothesize that this maximum is attained when the number of wet cages equals that of dry ones: in this condition one should have the maximum mismatch and strain. Simultaneously to the increase of strain one observes a decrease in the size of the largest coherent domains. Our interpretation of this empirical evidence is that wet cages form continuous, coherent domains—in a broad sense “droplets”—within the initially empty ZIF-8 crystallites. These droplets grow up to the point that the coherent domain of dry cages (“bubbles”) becomes the smallest coherent domain in the ZIF-8 crystallite; from this point on, the largest coherent domain becomes the wet one, which steadily grows (Figure 1b) while the dry domain shrinks. Summarizing, synchrotron data suggest that ZIF-8 is intruded by water following a cascade mechanism by progressive filling of connected cages, which account for the observed shrinkage and growth of coherent domains. To clarify the molecular mechanism of wetting of complex microporous systems, we performed molecular dynamics atomistic calculations in which a z-oriented ZIF-8 slab is immersed in bulk water, Figure 2a; a 25 MPa pressure is applied by two pistons parallel to the slab, according to the method proposed by Marchio et al.30 Time scales accessible by atomistic simulations, tens of nanoseconds, are too short to observe the wetting of the slab; thus, the process was accelerated using the restrained molecular dynamics (RMD). RMD forces the system to visit the microscopic states corresponding to an overall level of wetting of the ZIF-8 computational sample as measured by the number of water molecules within the MOF slab. It is worth remarking that these microscopic states are sampled with a statistical distribution consistent with the constant number of particles, pressure, and temperature ensemble and in quasi-static conditions, i.e., in conditions consistent with the experimental ones (see details of the method and the computational setup in the Supporting Information). RMD also allows us to compute the free energy profiles characterizing the wetting process.21 At a variance with previous simulations focusing on a three-periodic ZIF-8 bulk sample, with water condensation investigated by Grand Canonical Monte Carlo simulations,28 in which new water atoms are inserted at random locations within the material, our approach allows ZIF-8 wetting to take place according to most probable mechanism without excluding a priori neither condensation nor cascade penetration. Figure 2 a) Computational ZIF-8 sample, made by a (100) oriented slab comprising 7 layers of ZIF-8 cages. Each layer contains 4 cages, for a total of 28 complete cages. Cages are numbered in ascending order starting from the left layer of the slab, as indicated in the sketch. b) Number of water molecules per cage during the intrusion process at three levels of overall filling. These histograms show that intrusion starts from cages in contact with bulk water (cages 1–4 and 25–28) and proceeds toward the interior through a cascade, cage-by-cage filling mechanism. c) Free energy profile of wetting for the same ZIF-8 cage (denoted by the yellow circle in panel d) when (i) water is allowed to enter by all 5 wet cages surrounding the wetted one (blue) or (ii) water is allowed to enter only from one of them (red). The shadowed regions denote the (statistical) error, standard deviation. The difference between the two curves, well beyond the error, confirms that wetting from neighboring cages is favored as compared to water condensation, the latter depending only on the number of water molecules in the cage. d–g) Snapshots of the single cage wetting process in blue. Panel f corresponds to the transition state (TS) of the single cage intrusion process, the state corresponding to the maximum of the free energy (see also Figure S2 in the Supporting Information). To identify the actual wetting mechanism, we computed the number of molecules in each of the 28 cages comprising the ZIF-8 slab (Figure 2a). If intrusion took place via condensation in bulk cages, one should observe water molecules in the gas phase permeating the slab and condensing at random cages within the material. On the contrary, one notices that the wetting of the slab starts by the complete filling of specific cages belonging to the outermost layers, those in contact with bulk water, and proceeds by filling cages connected to the wet ones through 6MR apertures, which is consistent with the cascade intrusion mechanism proposed above. To further clarify the ZIF-8 wetting mechanism, we simulated water intrusion in a single cage. We started from the partly filled ZIF-8 slab in the state containing ∼960 water molecules and ran additional RMD simulations to refine the free energy profile for the filling of an individual cage (blue line in Figure 2c). Before the intrusion of this cage starts, i.e., when the number of water molecules #H2O/cage is zero, some water molecules inside the surrounding wet cages are present near the connecting 6MR apertures. When wetting starts, i.e., when #H2O/cage grows, these molecules penetrate the cage and the free energy grows. This can be explained considering that in this early stage the penetrating water molecules are undercoordinated; i.e., they form fewer hydrogen bonds than in the initial state.31 This energetically unfavorable process is driven by the external pressure. The free energy ceases to increase at the transition state (free energy maximum) which occurs when water molecules penetrating from neighboring wet cages can form hydrogen bonds between them (Figure 2f), which suddenly reduce the energy penalty associated with H2O undercoordination. This event requires that enough water molecules penetrate the cage, ca. 8 water molecules according to our simulations. From this point onward, the free energy decreases until a minimum is reached at #H2O/cage ∼38, which is consistent with the typical water molecule occupancy obtained from the histograms of Figure 2b for fully filled cages. To confirm that the wetting mechanism is due to cascade penetration rather than condensation, we performed additional free energy calculations for the wetting of a single cage but allowing water to enter from only one of the surrounding wetting cages (red line in Figure 2c). If wetting is by condensation, with the number of water molecules in the liquid-like water nucleus as the “reaction coordinate”, there should be no significant difference between the free energy profiles in the two cases. Instead, Figure 2c shows that allowing water to enter in a cage from one or multiple surrounding cages completely changes the free energy profile. Specifically, it is energetically more convenient to form a water bridge across several 6MR apertures with already filled ZIF-8 cavities. This finding confirms our hypothesis that the kinetically favored mechanism to wet ZIF-8 is through cascade penetration, with water penetrating from cage to cage, with the transition state corresponding to the formation of hydrogen bonds among water molecules entering from neighboring cages (Figure 2f, see also Figure S2). This brings us to propose that the key structural parameter controlling hysteresis is the pore connectivity, including the number and distance of connecting apertures through which water molecules can propagate capillary penetration. To some extent, these key characteristics are related to the cage size of the MOF, but the latter characteristic is not as important for determining hysteresis via limiting the size of critical liquid-like water nuclei, as previously proposed.28 The results of the single cage wetting simulation clearly show that the microscopic ZIF-8 wetting mechanism is not condensation but that the transition state is determined by the formation of a hydrogen bond bridge between water molecules penetrating from neighboring wet cages, rather than by the competition of free energy gain corresponding to larger liquid-like water clusters and the free energy penalty, which characterize condensation processes (see Supporting Information for a summary of classical models of condensation processes). The simulated microscopic penetration mechanism suggests that there is a correlation between the wetting state of neighboring cages, which can favor the formation of coherent domains of wet or dry cages, causing in turn different structural and mechanical characteristics (strain). It is helpful to explain the formation of these coherent domains in terms of the macroscopic theory of capillarity. Let us consider the case of n fully wet cages (∼40 H2O per cage) and focus on two cases, (i) one in which the wet cages are randomly distributed in the ZIF-8 crystallite and (ii) the other in which they form a single, droplet-like, coherent domain. According to our recent finding that water molecules can form hydrogen bonds across secondary apertures of porous systems,32 the droplet-like state is energetically favored. MD simulations reveal that hydrogen bonding across 6MR apertures stabilizes the system by 6.5 kBT (see the Supporting Information) per pair of adjacent wet cages connected by the aperture. Thus, the droplet-like configuration is energetically favored because it minimizes the number of wet cages in contact with dry ones. Indeed, the number of “interface” wet–dry 6MR apertures is proportional to the surface of the coherent domain of wet cages which tends to be minimized for energetic reasons. It is worth remarking that the argument developed here is analogous to the one at the basis of the capillarity, which predicts that isolated droplets are spherical to minimize surface energy; similarly to classical capillarity,33 intrusion in ZIF-8 is controlled by a competition between bulk-like volume terms favoring intrusion—mainly the liquid pressure—and the surface cost related to the interface of wet–dry cages which introduces a (free) energy penalty that disfavors the formation of many small wet domains. Atomistic simulations, with their (relatively) small computational sample preclude a direct empirical confirmation of the wetting mechanism on the crystallite scale. Thus, we developed a stochastic model in which cages undergo wetting or drying with a probability depending on the number of wet neighbors, coherently with the results in Figure 2c, with a ∼2.5 kBT reduction of the wetting barrier per surrounding wet cage (see Supporting Information for the details). The simulated crystallite has the same type of connectivity of ZIF-8. Mimicking experiments, in our model, we linearly increase the pressure in a ramp, although for computational reasons we use compression rate orders of magnitude faster than the typical experimental ones (106 MPa/s vs 3.16 × 10–6 MPa/s). The results of our model, shown in Figure 3, confirm the conclusion that the stabilization provided by hydrogen bonding across 6MR apertures results in the formation of one large coherent domain (see Supporting Information for a full movie of the process) while the random filling of cages by condensation is negligible. Because of the high rate at which the pressure is changed, the system forms straight “liquid fronts” connected by curved edges which bear memory of the square geometry and of the initial conditions. In an infinitely long calculation, one would reach the surface with minimal energy—a single spherical empty domain in the center. Additionally, we remark that the advancing flat front of Figure 3, or the single spherical domain mentioned in the previous sentence, is the result of some simplification adopted in the stochastic model, which is based on a defect-free, perfectly cubic, ZIF-8 crystallite, with all wet cages containing the same amount of liquid. This might change the fine details of the mechanism, whether the advancing from is regular or not, but not the overall, cage-by-cage intrusion mechanism. Figure 3 Stochastic model of water intrusion in a ZIF-8 cubic crystallite. Gray and blue squares represent dry and wet cages, respectively. One notices that, as pressure increases, water starts to intrude in the crystallite. Here, intrusion takes place according to a cascade penetration mechanism forming a small number of growing “droplets” of wet cages. These droplets eventually merge, and water completely fills the crystallite. In summary, via a combination of experimental and theoretical approaches we have clarified the intrusion mechanism of ZIF-8 under hydrostatic pressure, which proceeds by cascade filling of neighboring cages. In operando synchrotron radiation tracked the evolution of domains of wet/dry cages during intrusion, indicating the formation coherent domains of wet cages which grow during the process. RMD simulations gave microscopic evidence of this process, bringing to the conclusion that coherent domains are formed due to penetration of water in connected ZIF-8 cages, rather than water condensation in individual ones, as previously proposed. The detailed wetting mechanism of a single cage has been analyzed, showing that the kinetic bottleneck of the process is the formation of hydrogen bonds bridging water molecules across neighboring apertures. This suggests that tuning the intrusion or extrusion characteristics of hydrophobic MOFs can be obtained by optimizing the pore connectivity; in particular, the geometrical and chemical characteristics of the apertures are crucial to favor or disfavor the formation of hydrogen bonds, which, in turn, drive the cascade filling of neighboring cavities. Finally, on the crystallite scale, the shape of the domains of wet cages seems compatible with the classical capillarity concept of surface energy because they tend to minimize the surface area. Indeed, a stochastic model informed by MD simulations confirms the prevalence of cascade penetration in coherent domains as opposed to a condensation scenario. The critical characteristics responsible for the water intrusion mechanism if ZIF-8 is the presence of interconnected cages. Thus, we expect a similar mechanism to hold also for other ZIF MOFs. For example, for ZIF-67, which is an isomorph to ZIF-8, simulations show cage-by-cage intrusion (Figure S5). Other ZIFs, e.g., ZIF-12, with different morphologies and cage and aperture sizes, might show some differences with respect to the mechanism discussed above, which, however, we expect to be quantitative rather than qualitative. We plan to investigate this in the future. Supporting Information Available The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.nanolett.3c00235.Experimental porosimetry and synchrotron details and data, PV isotherm, atomistic simulations details, density maps, average volume of a ZIF-8 cage along intrusion, schematic model of the mechanism of water condensation, number of water molecules per cage during the intrusion process, and stochastic model additional information (PDF) Supplementary Material nl3c00235_si_001.pdf Author Contributions The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript. The authors declare no competing financial interest. Acknowledgments This project leading to this application has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No 101017858. This research is part of a project that has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No. 803213). The authors acknowledge PRACE for awarding us access to Marconi100 at CINECA, Italy. This article is part of the grant RYC2021-032445-I funded by MICIN/AEI/10.13039/501100011033 and by the European Union NextGenerationEU/PRTR. Synchrotron powder diffraction data were collected at beamline 17-BM at the Advanced Photon Source, Argonne National Laboratory. Use of the Advanced Photon Source was supported by the U. S. Department of Energy, Office of Science, Office of Basic Energy Sciences, under Contract No. DE-AC02-06CH11357. ==== Refs References Le Donne A. ; Tinti A. ; Amayuelas E. ; Kashyap H. K. ; Camisasca G. ; Remsing R. C. ; Roth R. ; Grosu Y. ; Meloni S. Intrusion and Extrusion of Liquids in Highly Confining Media: Bridging Fundamental Research to Applications. Adv. Phys. X 2022, 7 (1 ), 2052353 10.1080/23746149.2022.2052353. Lynch C. I. ; Rao S. ; Sansom M. S. P. Water in Nanopores and Biological Channels: A Molecular Simulation Perspective. Chem. Rev. 2020, 120 (18 ), 10298–10335. 10.1021/acs.chemrev.9b00830.32841020 Clough E. C. ; Plaisted T. A. ; Eckel Z. C. ; Cante K. ; Hundley J. M. ; Schaedler T. A. Elastomeric Microlattice Impact Attenuators. Matter 2019, 1 (6 ), 1519–1531. 10.1016/j.matt.2019.10.004. Fraux G. ; Coudert F.-X. ; Boutin A. ; Fuchs A. H. Forced Intrusion of Water and Aqueous Solutions in Microporous Materials: From Fundamental Thermodynamics to Energy Storage Devices. Chem. Soc. Rev. 2017, 46 , 7421–7437. 10.1039/C7CS00478H.29051934 Tinti A. ; Giacomello A. ; Grosu Y. ; Casciola C. M. Intrusion and Extrusion of Water in Hydrophobic Nanopores. Proc. Natl. Acad. Sci. U. S. A. 2017, 114 (48 ), E10266–E10273. 10.1073/pnas.1714796114.29138311 Wu H. ; Qian X. ; Zhu H. ; Ma S. ; Zhu G. ; Long Y. Controlled Synthesis of Highly Stable Zeolitic Imidazolate Framework-67 Dodecahedra and Their Use towards the Templated Formation of a Hollow Co3O4 Catalyst for CO Oxidation. RSC Adv. 2016, 6 (9 ), 6915–6920. 10.1039/C5RA18557B. Zhang Y. ; Jia Y. ; Li M. ; Hou L. Influence of the 2-Methylimidazole/Zinc Nitrate Hexahydrate Molar Ratio on the Synthesis of Zeolitic Imidazolate Framework-8 Crystals at Room Temperature. Scientific Reports 2018 8:1 2018, 8 (1 ), 1–7. 10.1038/s41598-018-28015-7. Khay I. ; Chaplais G. ; Nouali H. ; Marichal C. ; Patarin J. Water Intrusion-Extrusion Experiments in ZIF-8: Impacts of the Shape and Particle Size on the Energetic Performances. RSC Adv. 2015, 5 (40 ), 31514–31518. 10.1039/C5RA02636A. Ortiz G. ; Nouali H. ; Marichal C. ; Chaplais G. ; Patarin J. Versatile Energetic Behavior of ZIF-8 upon High Pressure Intrusion–Extrusion of Aqueous Electrolyte Solutions. J. Phys. Chem. C 2014, 118 , 7321–7328. 10.1021/jp412354f. Ortiz G. ; Nouali H. ; Marichal C. ; Chaplais G. ; Patarin J. Energetic Performances of the Metal–Organic Framework ZIF-8 Obtained Using High Pressure Water Intrusion–Extrusion Experiments. Phys. Chem. Chem. Phys. 2013, 15 (14 ), 4888–4891. 10.1039/c3cp00142c.23443335 Grosu Y. ; Li M. ; Peng Y.-L. ; Luo D. ; Li D. ; Faik A. ; Nedelec J.-M. ; Grolier J.-P. A Highly Stable Nonhysteretic {Cu2(Tebpz) MOF+water} Molecular Spring. ChemPhysChem 2016, 17 (21 ), 3359–3364. 10.1002/cphc.201600567.27442186 Furukawa H. ; Cordova K. E. ; O’Keeffe M. ; Yaghi O. M. The Chemistry and Applications of Metal-Organic Frameworks. Science (1979) 2013, 341 (6149 ), 1230444 10.1126/science.1230444. Ortiz G. ; Nouali H. ; Marichal C. ; Chaplais G. ; Patarin J. Energetic Performances of the Metal–Organic Framework ZIF-8 Obtained Using High Pressure Water Intrusion–Extrusion Experiments. Phys. Chem. Chem. Phys. 2013, 15 (14 ), 4888–4891. 10.1039/c3cp00142c.23443335 Ortiz G. ; Nouali H. ; Marichal C. ; Chaplais G. ; Patarin J. Energetic Performances of “ZIF-71-Aqueous Solution” Systems: A Perfect Shock-Absorber with Water. J. Phys. Chem. C 2014, 118 (37 ), 21316–21322. 10.1021/jp505484x. Grosu Y. ; Gomes S. ; Renaudin G. ; Grolier J. P. E. ; Eroshenko V. ; Nedelec J. M. Stability of Zeolitic Imidazolate Frameworks: Effect of Forced Water Intrusion and Framework Flexibility Dynamics. RSC Adv. 2015, 5 (109 ), 89498–89502. 10.1039/C5RA19879H. Grosu Y. ; Li M. ; Peng Y. L. ; Luo D. ; Li D. ; Faik A. ; Nedelec J. M. ; Grolier J. P. A Highly Stable Nonhysteretic {Cu2(Tebpz) MOF+water} Molecular Spring. ChemPhysChem 2016, 17 (21 ), 3359–3364. 10.1002/cphc.201600567.27442186 Park K. S. ; Ni Z. ; Côté A. P. ; Choi J. Y. ; Huang R. ; Uribe-Romo F. J. ; Chae H. K. ; O’Keeffe M. ; Yaghi O. M. Exceptional Chemical and Thermal Stability of Zeolitic Imidazolate Frameworks. Proc. Natl. Acad. Sci. U. S. A. 2006, 103 (27 ), 10186–10191. 10.1073/pnas.0602439103.16798880 Fraux G. ; Boutin A. ; Fuchs A. H. ; Coudert F.-X. Structure, Dynamics, and Thermodynamics of Intruded Electrolytes in ZIF-8. J. Phys. Chem. C 2019, 123 , 15589–15598. 10.1021/acs.jpcc.9b02718. Grosu Y. ; Mierzwa M. ; Eroshenko V. A. ; Pawlus S. ; Chorazewski M. ; Nedelec J. M. ; Grolier J. P. E. Mechanical, Thermal, and Electrical Energy Storage in a Single Working Body: Electrification and Thermal Effects upon Pressure-Induced Water Intrusion-Extrusion in Nanoporous Solids. ACS Appl. Mater. Interfaces 2017, 9 (8 ), 7044–7049. 10.1021/acsami.6b14422.28177602 Grosu Ya. ; Renaudin G. ; Eroshenko V. ; Nedelec J.-M. ; Grolier J.-P. E. Synergetic Effect of Temperature and Pressure on Energetic and Structural Characteristics of {ZIF-8 + Water} Molecular Spring. Nanoscale 2015, 7 (19 ), 8803–8810. 10.1039/C5NR01340B.25907279 Tortora M. ; Zajdel P. ; Lowe A. R. ; Chorążewski M. ; Leão J. B. ; Jensen G. v. ; Bleuel M. ; Giacomello A. ; Casciola C. M. ; Meloni S. ; Grosu Y. Giant Negative Compressibility by Liquid Intrusion into Superhydrophobic Flexible Nanoporous Frameworks. Nano Lett. 2021, 21 (7 ), 2848–2853. 10.1021/acs.nanolett.0c04941.33759533 Amabili M. ; Grosu Y. ; Giacomello A. ; Meloni S. ; Zaki A. ; Bonilla F. ; Faik A. ; Casciola C. M. Pore Morphology Determines Spontaneous Liquid Extrusion from Nanopores. ACS Nano 2019, 13 (2 ), 1728–1738. 10.1021/acsnano.8b07818.30653291 Beckstein O. ; Sansom M. S. P. Liquid-Vapor Oscillations of Water in Hydrophobic Nanopores. Proc. Natl. Acad. Sci. U. S. A. 2003, 100 (12 ), 7063–7068. 10.1073/pnas.1136844100.12740433 Powell M. R. ; Cleary L. ; Davenport M. ; Shea K. J. ; Siwy Z. S. Electric-Field-Induced Wetting and Dewetting in Single Hydrophobic Nanopores. Nat. Nanotechnol 2011, 6 (12 ), 798–802. 10.1038/nnano.2011.189.22036811 Smirnov S. N. ; Vlassiouk I. v. ; Lavrik N. v. Voltage-Gated Hydrophobic Nanopores. ACS Nano 2011, 5 (9 ), 7453–7461. 10.1021/nn202392d.21838311 Smirnov S. ; Vlassiouk I. ; Takmakov P. ; Rios F. Water Confinement in Hydrophobic Nanopores. Pressure-Induced Wetting and Drying. ACS Nano 2010, 4 (9 ), 5069–5075. 10.1021/nn101080k.20690599 Tinti A. ; Camisasca G. ; Giacomello A. Structure and Dynamics of Water Confined in Cylindrical Nanopores with Varying Hydrophobicity. Philosophical Transactions of the Royal Society A 2021, 379 (2208 ), 20200403 10.1098/rsta.2020.0403. Sun Y. ; J Rogge S. M. ; Lamaire A. ; Vandenbrande S. ; Wieme J. ; Siviour C. R. ; van Speybroeck V. ; Tan J.-C. High-Rate Nanofluidic Energy Absorption in Porous Zeolitic Frameworks. Nat. Mater. 2021, 20 , 1015–1023. 10.1038/s41563-021-00977-6.33888902 Fairen-Jimenez D. ; Moggach S. A. ; Wharmby M. T. ; Wright P. A. ; Parsons S. ; Düren T. Opening the Gate: Framework Flexibility in ZIF-8 Explored by Experiments and Simulations. J. Am. Chem. Soc. 2011, 133 (23 ), 8900–8902. 10.1021/ja202154j.21553843 Marchio S. ; Meloni S. ; Giacomello A. ; Casciola C. M. Wetting and Recovery of Nano-Patterned Surfaces beyond the Classical Picture. Nanoscale 2019, 11 (44 ), 21458–21470. 10.1039/C9NR05105H.31686077 Paulo G. ; Gubbiotti A. ; Grosu Y. ; Meloni S. ; Giacomello A. The Impact of Secondary Channels on the Wetting Properties of Interconnected Hydrophobic Nanopores. Commun. Phys. 2023, 6 (21 ), 21 10.1038/s42005-023-01140-0. Bushuev Y. G. ; Grosu Y. ; Chorążewski M. A. ; Meloni S. Subnanometer Topological Tuning of the Liquid Intrusion/Extrusion Characteristics of Hydrophobic Micropores. Nano Lett. 2022, 22 (6 ), 2164–2169. 10.1021/acs.nanolett.1c02140.35258978 Giacomello A. ; Casciola C. M. ; Grosu Y. ; Meloni S. Liquid Intrusion in and Extrusion from Non-Wettable Nanopores for Technological Applications. European Physical Journal B 2021 94:8 2021, 94 (8 ), 1–24. 10.1140/epjb/s10051-021-00170-3.