Thermal and Guest-Assisted Structural Transition in the NH 2 -MIL-53(Al) Metal Organic Framework: A Molecular Dynamics Simulation Investigation

Transcription

1 nanomaterials Article Thermal and Guest-Assisted Structural Transition in the NH -MIL-53(Al) Metal Organic Framework: A Molecular Dynamics Simulation Investigation Roald Boulé 1,, Claire Roland, Laurent Le Pollés, Nathalie Audebrand and Aziz Ghoufi 1, * 1 IPR (Institut de Physique de Rennes) UMR 651, CNRS, University Rennes, F-35 Rennes, France; roald.boule@gmail.com ISCR (Institut de Sciences Chimiques de Rennes) UMR 66, CNRS, University Rennes, F-35 Rennes, France; claire.roiland@univ-rennes1.fr (C.R.); laurent.le-polles@ensc-rennes.fr (L.L.P.); nathalie.audebrand@univ-rennes1.fr (N.A.) * Correspondence: aziz.ghoufi@univ-rennes1.fr; Tel.: Received: 1 June 18; Accepted: 1 July 18; Published: 1 July 18 Abstract: Reversible structural transition between the Large (LP) and Narrow Pore (NP) forms (breathing phenomena) of the MIL-53(X, X = Al, Cr, Fe, Ga) Metal Organic Framework (MOF) is probably one of the most amazing physical properties of this class of soft-porous materials. Whereas great attention has been paid to the elucidation of the physical mechanism ruling this reversible transition, the effect of the functionalization on the flexibility has been less explored. Among functionalized MIL-53(Al) materials, the case of NH -MIL-53(Al) is undoubtedly a very intriguing structural transition rarely observed, and the steadier phase corresponds to the narrow pore form. In this work, the flexibility of the NH -MIL-53(Al) metal organic framework was investigated by means of molecular dynamics simulations. Guest (methanol) and thermal breathing of the NH -MIL-53(Al) was thus explored. We show that it is possible to trigger a reversible transition between NP and LP forms upon adsorption, and we highlight the existence of stable intermediate forms and a very large pore phase. Furthermore, the NP form is found thermodynamically stable from to K, which is the result of strong intramolecular hydrogen bonds. Keywords: NH -MIL-53(Al); MOFs; molecular simulation; structural transition 1. Introduction The flexibility of soft-porous Metal-Organic Frameworks (MOFs) is probably one of the most fascinating structural features in the field of nanoporous materials in comparison with the most common materials such as active carbons and zeolites [1 6]. This flexibility is the result of the connection of the organic ligands with inorganic building units. For several MOF materials, this flexibility induces structural transitions triggered by several external stimuli such as light [7 9], electrical fields [1,11], mechanical pressure [1 1] and adsorption of gases [ 6,15 3]. Interestingly, it has been recently shown that this structural modulation (breathing, pore gating, ligand flip, etc.) allowed new adsorption and separation properties involving new areas in adsorption processes [ 6,15 3]. Among the flexible MOFs, the MIL-53 material is probably the most studied in this past decade with evidence of a guest-, thermal- and mechanically-assisted reversible structural transition between Narrow Pore (NP) and Large Pore (LP) forms [,19,,,5]. Indeed, a reversible expansion/contraction of the pore dimension up to % in unit cell volume was observed whatever the external stimuli. At room temperature and atmospheric pressure, the steadier phase of the MIL-53(Al) material is the LP form, i.e., an orthorhombic structure with a large pore of 1 Å, inducing a unit cell volume of 18 Å 3. By applying an external (electrical field, mechanical pressure, etc) or Nanomaterials 18, 8, 531; doi:1.339/nano

2 Nanomaterials 18, 8, 531 of 1 internal (adsorption of gases) constraint, the MIL-53(Cr) is then thermodynamically and mechanically destabilized. To reach a new thermodynamical equilibrium, a structural transition occurs to trend toward a more steadier phase in line with this constraint. The pores of the MIL-53(Al) contract; the symmetry changes; and the volume decreases from 18 Å 3 to 1 Å 3 : the NP form. When the constraint is released, the initial LP form is then recovered, emphasizing the reversible aspect of this structural transition: the breathing effect. Recently, the use of this family of materials was suggested for mechanical energy storage applications [1 1]. Whereas the MIL-53(Al) was intensively studied, the flexibility of the modified MIL-53 was less investigated. The structural behaviors of MIL-53 have been recently analyzed as a function of functionalized organic linkers [6]. It was thus highlighted that these chemical modifications had a deep impact on the breathing phenomena (reversible structural transition). Among functionalized aluminum hydroxo terephthalates (Al(OH)(BDC-X)) (with X -H, -CH 3, -Cl, -Br, -NH, -NO, -(OH), -CO H) crystallizing in the MIL-53-type structure, it was recently exhibited that NH -MIL-53(Al) did not present structural transition, which explains the little attention given to its flexibility [7,8], although Ahmfelt and coworkers have found similar breathing behavior to that unfunctionalized MIL-53 [9]. Herein, we computationally investigate then, for the first time, the effect of the temperature and the adsorption of a polar guest (methanol) on the structure and the flexibility of the NH -MIL-53(Al). Let us mention that, although the comparison of the effect of the adsorption of polar and no polar molecules on the flexibility of NH -MIL-53(Al) can be very interesting, it is not the scope in this work, and that is under investigation.. Materials and Methods The simulation box was built from 3 unit cells of NH -MIL-53(Al) along the x, y and z directions. The initial configuration was built by grafting the NH group in the ortho position of the benzenic cycles of the LP form of the MIL-53(Al) material. Atomic positions of the LP form correspond to the experimental coordinates previously obtained from X-ray powder diffraction []. An illustration of the LP-NH -MIL-53(Al) is provided in Figure 1. NH-MIL-53(Al) in an orthorhombic LP Form (18 Å 3 at 3K) NH-MIL-53(Al) in a monoclinic NP Form (135 Å 3 at 3K) b" c" a" Initial configuration Steadier phase at room temperature Figure 1. Illustration of the initial Narrow Pore (NP) and Large Pore (LP) forms of the NH -MIL-53(Al) metal organic framework. Experimentally speaking, it has been well known that the NH -MIL-53(Al) is thermodynamically stable in the NP form [7,3]. We took the route to use the LP form as the initial configuration in order to ensure that our force field was capable of predicting the steadier phase. The initial box lengths correspond to the 3 unit cells of the LP-MIL-53(Al), L x = 33. Å, L y = 5.8 Å and L z = 6.8 Å. All molecular dynamics simulations were performed by using the DLPOLY software (Version.) [31]. The velocity-verlet algorithm in the NσT statistical ensemble was used. N is the number of molecules, T the temperature and σ the isotropic constraint (1 bar). The Nose Hoover thermostat and barostat was used with a relaxation time of.5 ps. Acquisition phases were performed during 5 ns by using a time step of 1 fs, whereas the equilibration phases lasted.5 ns. The Ewald summation was used for

3 Nanomaterials 18, 8, of 1 calculating the electrostatic interactions. Long- and short-range interactions were truncated at 1 Å. The NH -MIL-53(Al) framework was considered as flexible and was modeled by refining the force field of the MIL-53(Cr) [19]. The force fields of the NH group and Al atom were considered by combining the DREIDING [3] and AMBER [33] models. Partial charges were taken in [3], where calculations were carried out on the cluster represented in Figure. Methanol molecules were modeled by means of the OPLS (Optimized Potentials for Liquid Simulations) force field [35] because it has been shown that the structure and physical properties of liquid and gas were well reproduced [35]. The van der Waals interactions between the NH -MIL-53(Al) material and methanol molecules were calculated by using the mixing rules of Lorentz and Berthelot. Let us mention that the methanol molecules were randomly inserted. All parameters are provided in ForceFIELD.txt (the force field), CONFIG.txt (coordinates of the initial configuration) and CONTROL.txt (thermodynamics parameters and integration algorithms). These files are the input files of the DLPOLY [31] software. OH# HO# Al# OC# HN# Figure. Cluster of atoms used in the charge calculation [3]. Labels of atoms (OH, HO, Al, OC, HN) implicated in the radial distribution functions calculation are also represented. 3. Results and Discussion We report in Figure 3a,b the unit cell volume (V u.c. ) and the unit cell parameters (a, b, c) of the NH -MIL-53(Al) as a function of the temperature. As shown in Figure 3a, V u.c. slightly increases to 3 Å 3 from to K. In this range of volume, the NH -MIL-53(Al) is in the NP form. This result is in good agreement with the experiment where the LP form was not detected whatever the temperature [7]. Interestingly, a deviation of 8% of V u.c. between simulation and experiment was found that was a source of force field validation. Additionally, the fact of finding a stable NP form that is in line with the experiment allowed us to be confident in the so-developed force field. Although the NP phase was systematically recovered, a slight increase in V u.c. was observed that could suggest a possible reopening. As shown in Figure 3a, this evolution was fully reversible, and no hysteresis was observed. Let us note that at 5 K, an NP form of 165 Å 3 was found that bore out the continuous increase in V u.c. and the possibility of a reopening. Actually, the thermal energy (k B T where k B is the Boltzmann constant) is not sufficient to activate the reopening. Probably the internal energy of the NH -MIL-53(Al) has to be too high to be surpassed by the thermal one, which suggests strong intramolecular interactions. As shown in Figure 3a, a strong difference is observed with the MIL-53(Al) material because a structural LP NP transition is observed at 7 K. This difference between both materials can only be explained by the presence of the NH groups and is probably the result of the

4 Nanomaterials 18, 8, 531 of 1 interactions between NH and other atoms such as oxygen and hydrogen atoms of the terephthalate moiety. Furthermore, Figure 3b shows that the increase in V u.c. is due to the increase in the b value, whereas a decreases, and c is almost kept constant. This evolution of the unit cell parameters could suggest a possible reopening because it has been shown that the transition between the NP and LP phase is connected to the decrease of a, the increase of b and the quasi-conservation of c (see Table 1 of [19]). a)# b)# V exp = A a (A ) b & c (A ) 13 NH-MIL-53(Al) from K to K NH-MIL-53(Al) from K to K MIL-53(Al) Figure 3. (a) Unit cell volume as a function of the temperature for the NH -MIL-53(Al) (left axis) and for the MIL-53(Al) (right axis) materials; (b) Unit cell parameters; a (left axis), b (right axis) and c (right axis) as a function of the temperature for the NH -MIL-53(Al) material. On average, unit cell parameters were calculated within a relative error of. Å. We manage in Figure a the radial distribution functions (RDF) between the oxygen atoms of carboxylate groups (labeled as OC in Figure ) and the hydrogen atoms of NH groups (named HN in Figure ). For a central particle, RDF describes how the local density evolves as a function of distance. From statistical physics, RDF corresponds to the probability of finding an atom at a distance of r from a central atom. As shown in Figure a, the first peak is located at. Å from to K, which involves strong intramolecular interactions of the same order of magnitude of a hydrogen bond that is located at.5 Å. This high affinity between OC and HN is probably at the origin of the stability of the NP form from to K. The RDF between OC and HO (see Figure for the definition of labels) is also reported in Figure b. As exhibited in Figure b, the first peak is located around.6 Å, which suggests an interaction weaker than that of OC-HN, but strong enough to contribute to the stability of the NP phase. This OC-HO interaction suggests that the alumina octahedron (Al(OC) (OH) ) could be slightly deformed. We calculated then the average angles of the octahedra, OC-Al-OC, OC-Al-OH and OH-Al-OH. We found OC Al OC = and 9., OC Al OH = 88.8 and OH Al OH = These values are close to the ideal octahedra, 9 and 18, which highlights an undeformed octahedral geometry. Additionally, let us note that the distance of 3. Å between OC and HO is also recovered in the LP of the MIL-53(Al), resulting in that this interaction cannot be at the origin of the stability of the NP form of NH -MIL-53(Al). Therefore, the stable NP form is the result of the strong intramolecular hydrogen bonds between HN and OC. Additionally, as shown in Figure b, the interactions between OH and HN and N and HO seem to not take part in the stabilization of the NP form.

5 Nanomaterials 18, 8, of 1 a) K 3K K 3 RDF (OC-HN) R (A ) b) OC-HN (K) OC-HO (K) OH-HN (K) N-HO (K) 3 RDF R (A ) Figure. (a) Radial distribution functions between OC and HN (see Figure for labels definition) at three temperatures. In the inset on top right an enlargement for the short distances was provided; (b) Radial distribution functions between OC and HN, OC and HO, OH and HN and N and HO at K. These results obviously bear out that the empty NH -MIL-53(Al) is structurally and energetically stable under its NP form from to K. To assess the possibility of a reopening of the evolution of V u.c. upon adsorption of a polar compound, methanol (MeOH) was investigated. MeOH molecules were then progressively inserted in the empty configuration NP-NH -MIL-53(Al) at 3 K. We began by carrying out molecular dynamics (MD) simulation of one molecule by unit cell (1 molec./u.c.), NH -MIL-53(Al)@1 from 3 to K. From the final configuration of NH -MIL-53(Al)@1, additional MeOH molecules were inserted to reach 3 molec./u.c., NH -MIL-53(Al)@3. As shown in Figure 5, this process was repeated to finally fill the material at 8 moles./u.c. We report in Figure 6 the unit cell volume as a function of the temperature for all cases. As shown in Figure 6a, for NH -MIL-53(Al)@1, an increase of 6% of V u.c. with respect to the empty material is observed. From an overall point of view, Figure 6a exhibits a progressive increase of V u.c. from 1 to 8 moles./u.c.. As reported in Figure 6, the unit cell volume evolves from 1 Å 3 (in the empty phase) to 16 Å 3 via intermediate volumes, 11 Å 3, 1 Å 3, 13 Å 3, 15 Å 3. These results highlight that NH -MIL-53(Al) can be opened upon adsorption. Interestingly, Figure 6b,c shows that the reopening does not continue because a transition is observed between 36 K and 37 K for both NH -MIL-53(Al)@1 and NH -MIL-53(Al)@5. The additional reopening is then facilitated by a thermal activation that increases the local intramolecular disorder and triggers an additional increase in unit cell volume. Unexpectedly, Figure 6d shows a large increase in V u.c. for NH -MIL-53(Al)@8. Indeed, a very large pore of a unit cell volume of 16 Å 3 is observed. This over-opening is the result of the combination

6 Nanomaterials 18, 8, of 1 of the strong steric effect of a large amount of methanol molecules and thermal energy involving an increase in the intramolecular motion. Thermodynamically speaking, both of these contributions induce a destabilization of the structure and trigger a transition toward a steadier phase. The domain of existence of these phases (narrow, intermediate, large and very large pore) with respect to the temperature is highlighted in Figure 7a. STARTING Final empty configuration from 3K to K Insertion of 1 molec./u.c. è NH-MIL-53(Al)@1 MD simulation in NσT ensemble from 3K to K Insertion of molec./u.c. from final configuration of NH-MIL-53(Al)@1 è NH-MIL-53(Al)@3 MD simulation in NσT ensemble from 3K to K Insertion of 3 molec./u.c. from final configuration of NH-MIL-53(Al)@3 è NH-MIL-53(Al)@5 MD simulation in NσT ensemble from 3K to K Insertion of 3 molec./u.c. from final configuration of NH-MIL-53(Al)@5 è NH-MIL-53(Al)@8 MD simulation in NσT ensemble from 3K to K END Figure 5. Scheme illustrating our protocol to evaluate the reopening of NH -MIL-53(Al). molec./u.c., molecules by unit cell. a) molec./u.c 1 molec./u.c 3 molec./u.c 5 molec./u.c 8 molec./u.c b) molec./u.c 1 molec./u.c c) d) molec./u.c 5 molec./u.c molec./u.c Figure 6. (a) Unit cell volume as a function of the temperature for, 1, 3, 5 and 8 molec./u.c.; (b) Unit cell volume as a function of the temperature for and 1 molec./u.c.; (c) for 3 and 5 molec./u.c. and (d) for 8 molec./u.c.

7 Nanomaterials 18, 8, of 1 a)# molec./u.c. 1 molec./u.c. 3 molec./u.c. 5 molec./u.c. 8 molec./u.c. Very Large Pore Large Pore Intermediate Pore b (A ) b)# β ( o ) Narrow Pore 3 36 T (K) 38 8 b β 6 Adsorbed amount (molec./u.c.) Figure 7. (a) Unit cell volume of the NH -MIL53(Al) as a function of the temperatures and adsorbed amount; (b) Unit cell parameters (b and β) as a function of the adsorbed amount at 3 K. Actually, the swelling of NH -MIL-53(Al) can be imputed to a steric effect due to the MeOH molecules, which increase the pore size. As shown in Figure 1, the pore size of NH -MIL-53(Al) can be quantified by the the unit cell parameter along the y direction, i.e., b. We report then in Figure 7b the b value as a function of the adsorbed amount. As shown in Figure 7b, an increase in b is observed, which suggests an opening of pore as a function of the adsorbed amount. Interestingly, Figure 7b also depicts that these intermediate phases present an NP symmetry, i.e., monoclinic. However, from 8 moles./u.c., the orthorhombic symmetry, i.e., the LP symmetry, is recovered. To capture the microscopic insight into the reopening, the RDF between HN and OC are reported in Figure 8b for, 1, 3, 5 and 8 molec./u.c. at 3K. As shown in Figure 8b, the location of the first peak that represents the strength of the interaction (short distances involve strong interactions) is slightly modified in the presence of MeOH and during the reopening. Furthermore, the amplitude of the first peak is also un-modified, which suggests that the HN-OC interaction takes place whatever the loading. The reopening is probably due to a steric effect such that the adsorbed MeOH exerted an internal pressure that surpasses the interactions at the origin of the closed structure [16] (which is illustrated in Figure 8a). The kinetic radius of MeOH is around 3.5 Å, whereas the pore radius of the empty material was evaluated at.5 to.7 Å (this value was obtained from the calculation of the Connolly volume of one pore). To obtain an adjustment between the pore radius of NH -MIL53(Al) and the kinetic diameter of MeOH, the swelling is triggered. Additionally, as shown in Figure 8c, the adsorption of MeOH into the NH -MIL53(Al) is energetically favorable given the strong hydrogen bond interaction between the hydroxyl group of MeOH and the OH group of NH -MIL-53(Al), which is in favor of its adsorption. Indeed, Figure 8c shows an interaction located at. Å, i.e., a shorter distance than HN-OC (.8 Å). At the same time, the methyl group of MeOH (labeled Cm) is in a hydrophobic interaction with the carbon atoms of the benzenic cycles (named Cp). Indeed, as shown in Figure 8d, a distance of.5 Å was found as the location of the first peak, which is in line with the typical hydrophobic interaction. Actually, these hydrophobic interactions are in favor of the swelling of NH -MIL-53(Al).

8 Nanomaterials 18, 8, of 1 a) swelling R =.3-.7 Å + R = 3.5 Å Strong hydrogen bond (adsorp6on) Hydrophobic interac6on (swelling) = swelling b) c) r(å) HO-OHm HN-OC RDF (Cp-Cm) RDF (HO-OHm) RDF (HN-OC) RDF (HN-OC) molec./u.c. 1 molec./u.c. 3 molec./u.c. 5 molec./u.c. 8 molec./u.c.. d) r (Å) r (Å) Figure 8. (a) Scheme illustrating the swelling from the adsorption of the MeOH molecule in the empty NH -MIL53(Al); (b) Radial distribution functions between HN and OC for, 1, 3, 5 and 8 molec./u.c. at 3 K; (c) Radial distribution functions between the oxygen atom of the hydroxyl group of MeOH (Om) and HO and between HN and OC for 8 molec./u.c. at 3 K; (d) RDF between the carbon atoms of the benzenic groups of the MOF (Cp) and the carbon atoms of MeOH (Cm) at 3 K.. Conclusions In this work, the thermal structural transition of the NH -MIL-53(Al) metal organic framework was investigated by means of molecular dynamics simulations. The force field of the NH -MIL-53(Al) was developed by refining the force field of the MIL-53(Cr). For the empty material, the NP form of the NH -MIL-53(Al) was found to be steadier, in good agreement with the experiment. This stability and the absence of reopening are the result of the strong affinity between the oxygen atoms of the carboxylate groups and the hydrogen atoms of the amino groups. We showed that the adsorption of methanol triggered a progressive reopening of the structure by generating intermediate phases. We highlighted that this reopening is the result of a combination between hydrophobic interactions and strong hydrogen bonds. In future work, the differences between the influence of methanol and other molecules, on the one hand, water and, on the other hand, molecules that are not capable of hydrogen bond formation, for example dichloromethane or methane, will be investigated. Author Contributions: A.G. wrote the paper and designed the simulations. R.B. carried out the simulations. A.G. analyzed the data. A.G., N.A., L.L.P. and C.R. supervised the project. All authors commented on and approved the paper. Funding: A.G. is grateful to the ANR for its financial support through the project MEACOPA, ANR-17-CE9-3. Conflicts of Interest: The authors declare no conflict of interest. References Ferey, G. Giant flexibility of crystallized organic inorganic porous solids: Facts, reasons, effects and applications. New J. Chem. 16,, [CrossRef] Serre, C.; Férey, G. Large breathing effects in three-dimensional porous hybrid matter: Facts, analyses, rules and consequences. Chem. Soc. Rev. 9, 38, Horike, S.; Shimomora, S.; Kitagawa, S. Soft porous crystals. Nat. Chem. 9, 1, 695. [CrossRef] [PubMed]

9 Nanomaterials 18, 8, of 1. Férey, G.; Serre, C.; Devic, T.; Maurin, G.; Jobic, H.; Llewellyn, P.; Weireld, G.D.; Vimont, A.; Daturi, M.; Chang, J. Why hybrid porous solids capture greenhouse gases? Chem. Soc. Rev. 11,, [CrossRef] [PubMed] 5. Kitawaga, S.; Kitaura, R.; Noro, S. Functional Porous Coordination Polymers. Angew. Chem. Int. Ed., 3, [CrossRef] [PubMed] 6. Schneemann, A.; Bon, V.; Schwedler, I.; Senkovska, I.; Kaskel, S.; Fischer, R. Flexible metal-organicframeworks. Chem. Soc. Rev. 1, 16, [CrossRef] [PubMed] 7. Modrow, A.; Zargarani, D.; Herges, R.; Stock, N. The first porous MOF with photoswitchable linker molecules. Dalton Trans. 11,, 17. [CrossRef] [PubMed] 8. Yanai, N.; Uemura, T.; Inoue, M.; Matsuda, R.; Fukushima, T.; Tsujimoto, T.; Isoda, S.; Kitagawa, S. Guest-to-host transmission of structural changes for stimuli-response adsorption property. J. Am. Chem. Soc. 1, 1, [CrossRef] [PubMed] 9. Brown, J.W.; Henderson, B.L.; Kiesz, M.D.; Whalley, A.C.; Morris, W.; Grunder, S.; Deng, H.; Furukawa, H.; Zink, J.; Stoddart, J.; et al. Photophysical pore control in an azobenzene containing metal-organic framework. Chem. Sci. 13,, [CrossRef] 1. Ghoufi, A.; Benhamed, K.; Bouki-Hacene, L.; Maurin, G. Electrically induced Breathing of the MIL-53(Cr) Metal-Organic Framework. ACS Cent. Sci. 17, 3, [CrossRef] [PubMed] 11. Knebel, A.; Geppert, B.; Volgmann, K.; Kolokolov, D.; Stepanov, A.; Twiefel, J.; Heitjans, P.; Volkmer, D.; Caro, J. Defibrillation of soft porous metal-organic frameworks with electric fields. Science 17, 358, [CrossRef] [PubMed] 1. Yot, P.; Boudene, Z.; Macia, J.; Granier, D.; Vanduyfhuys, L.; Verstraelen, T.; Speybroeck, V.V.; Devic, T.; Serre, C.; Férey, G.; et al. Metal-organic framework as potential shock absorbers, the case of the highly flexible MIL-53(Al). Chem. Commun. 1, 5, [CrossRef] [PubMed] 13. Rodriguez, J.; Beurroies, I.; Loiseau, T.; Denoyel, R.; Llewellyn, P. The direct heat measurement of mechanical energy storage metal-organic frameworks. Angew. Chem. Int. Ed. 15, 5, [CrossRef] [PubMed] 1. Yot, P.; Vanduyfhus, L.; Alvarez, E.; Rodriguez, J.; Itié, J.; Fabry, P.; Guillou, N.; Devic, T.; Beurroies, I.; Llewellyn, P.; et al. Mechanical energy storage performance of an aluminum fumarate metal organic framework. Chem. Sci. 16, 7, 6 5. [CrossRef] [PubMed] 15. Salles, F.; Ghoufi, A.; Maurin, G.; Bell, R.; Mellot-Draznieks, C.; Llewellyn, P.; Serre, C.; Ferey, G. Molecular Dynamics Simulations of Breathing MOFs, Structural transformations of MIL-53(Cr) upon Thermal Activation and CO Adsorption. Angew. Chem. Int. Ed. 8, 7, [CrossRef] [PubMed] 16. Ghoufi, A.; Maurin, G.; Férey, G. Physics Behind the Guest-Assisted Structural Transitions of a Porous Metal-Organic Framework Material. J. Phys. Chem. Lett. 1, 1, [CrossRef] 17. Boutin, A.; Coudert, F.; Springuel-Huet, M.; Neimark, A.; Férey, G.; Fuchs, A. The behavior of Flexible MIL-53(Al) upon CH and CO Adsorption. J. Phys. Chem. C 1, 11, 37. [CrossRef] 18. Gaab, M.; Trukhan, N.; Maurer, S.; Gummaraju, R.; Müller, U. The progression of Al-based metal-organic frameworks From academic research to industrial production and applications. Microporous Mesoporous Mater. 1, 157, [CrossRef] 19. Ghoufi, A.; Subercaze, A.; Ma, Q.; Yot, P.; Ke, Y.; Puente-Orench, I.; Devic, T.; Guillerm, V.; Zhong, C.; Serre, C.; et al. Comparative Guest, Thermal, and Mechanical Breathing of the Porous Metal Organic Framework MIL-53(Cr): A Computational Exploration Supported by Experiments. J. Phys. Chem. C 1, 116, [CrossRef]. Ma, Q.; Yang, Q.; Ghoufi, A.; Férey, G.; Zhong, C.; Maurin, G. Guest dependent pressure behavior of the flexible MIL-53(Cr), A computational approach. Dalton Trans. 1, 1, [CrossRef] [PubMed] 1. Mason, J.; Oktawiec, J.; Taylor, M.; Hudson, R.; Rodriguez, J.; Bachman, J.; Gonzalez, M.; Cervellino, A.; Guagliardi, A.; Brown, C.; et al. Methane storage in flexible metal-organic frameworks with intrinsic thermal management. Nature 16, 57, 357. [CrossRef] [PubMed]. Bon, V.; Kavoosi, N.; Senkovska, I.; Kaskel, S. Tolerance of Flexible MOFs toward Repeated Adsorption Sress. ACS Appl. Mater. Interfaces 15, 7, 9 3. [CrossRef] [PubMed] 3. Krause, S.; Bon, V.; Senkovska, I.; Stoeck, U.; Wallacher, D.; Tobbens, D.; Zander, S.; Pillai, R.; Maurin, G.; Coudert, F.; et al. A pressure-amplifying framework material with negative gas adsorption transitions. Nature 16, 53, 38. [CrossRef] [PubMed]

10 Nanomaterials 18, 8, of 1. Serre, C.; Millange, F.; Thouvenot, C.; Noguès, M.; Marsolier, G.; Louer, D.; Férey, G. Very High Breathing effect in the First Nanoporous Chromium(III)-based Solids: MIL-53 or Cr(OH).OC-C6H- CO.HOC-C6H-COHx.HOy. J. Am. Chem. Soc., 1, [CrossRef] [PubMed] 5. Serre, C.; Bourrelly, S.; Vimont, A.; Ramsahye, N.; Maurin, G.; Llewellyn, P.; Daturi, M.; Filinchuk, Y.; Leynaud, O.; Garnes, P.; et al. An Explanation for the Very Large Breating Effect of a Metal-Organic Framework during CO Adsorption. Adv. Mater. 7, 19, [CrossRef] 6. Devic, T.; Horcajada, P.; Serre, C.; Salles, F.; Maurin, G.; Moulin, B.; Heurtaux, D.; Clet, G.; Vimont, A.; Grenèche, J.; et al. Functionalization in Flexible Porous Solids: Effects on the Pore Opening and the Host-Guest Interactions. J. Am. Chem. Soc. 1, 13, [CrossRef] [PubMed] 7. Munn, A.; Pillai, R.; Biswas, S.; Stock, N.; Maurin, G.; Walton, R. The flexibility of modified-linker MIL-53 materials. Dalton Trans. 16, 5, [CrossRef] [PubMed] 8. Lescouet, T.; Kockrick, E.; Bergeret, G.; Pera-Titus, M.; Farrusseng, D. Engineering MIL-53(Al) flexibility by controlling amino tags. Dalton Trans. 11,, [CrossRef] [PubMed] 9. Ahnfelt, T.; Gunzelmann, D.; Loiseau, T.; Hirsemann, D.; Senker, J.; Ferey, G.; Stock, N. Synthesis and Modification of a Functionalized 3D Open-Framework Structure with MIL-53 Topology. Inorg. Chem. 9, 8, [CrossRef] [PubMed] 3. Stavitski, E.; Pidko, E.A.; Couck, S.; Remy, T.; Hensen, E.J.; Weckhuysen, B.; Denayer, J.; Gascon, J.; Kapteijn, F. Complexity behind CO Capture on NH -MIL-53(Al). Langmuir 11, 7, [CrossRef] [PubMed] 31. Todorov, I.; Smith, W.; Trachenko, K.; Dove, M. DLPOLY3: New dimensions in molecular dynamis simulations via massive parallelism. J. Mater. Chem. 6, 16, [CrossRef] 3. Mayo, S.L.; Olafson, B.D.; Goddard, W. DREIDING: A generic force field for molecular simulations. J. Phys. Chem. 199, 9, [CrossRef] 33. Cornell, W.; Cleplak, P.; Bayly, C.; Gould, I.; Merz, K.; Fergusson, D.; Spellmeyer, D.; Fox, T.; Caldwell, J.; Kollman, P. A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules. J. Am. Chem. Soc. 1995, 117, [CrossRef] 3. Garcia-Pérez, E.; Serra-Crespo, P.; Hamad, S.; Kapteijn, F.; Gascon, J. Molecular simulation of gas adsorption and diffusion in a breathing MOF using a rigid force field. Phys. Chem. Chem. Phys. 1, 16, [CrossRef] [PubMed] 35. Jorgensen, W.; Maxwell, D.; Tirado-Rives, J. Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids. J. Am. Chem. Soc. 1996, 118, [CrossRef] c 18 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (