Source: https://pubs.rsc.org/en/content/articlehtml/2019/sc/c9sc01299k
Timestamp: 2019-04-22 14:10:56+00:00

Document:
We describe a general model to explore responsive adsorption processes in flexible porous materials. This model combines mean field formalism of the osmotic potential, classical density functional theory of adsorption in slit pore models and generic potential functions which represent the Helmholtz free energy landscape of a porous system. Using this model, we focus on recreating flexible adsorption phenomena observed in prototypical metal–organic frameworks, especially the recently discovered effect of negative gas adsorption (NGA). We identify the key characteristics required for the model to generate unusual adsorption processes and subsequently employ an extensive parametric study to outline conditions under which gate-opening and NGA are observed. This powerful approach will guide the design of responsive porous materials and the discovery of entirely new adsorption processes.
SPC materials that are responsive to the specific stimulus of gas adsorption typically exhibit one or more distinct steps in their gas adsorption isotherms, attributed to structural transformation of the host material.10 Examples of these responsive adsorption processes include gate-opening and breathing phenomena demonstrated by metal–organic framework (MOF) materials.11,12 Gate-opening refers to a particular (and sometimes subtle) change in the structure of the porous network at a specific adsorptive pressure.13 This change in the structure “opens” previously inaccessible pore volume to the adsorbing molecules producing a transition from low porosity narrow pore (np) state to more porous open pore (op) state, leading to a stepwise change in the adsorption isotherm. Recently, it was demonstrated that this effect can be exploited to maximize the deliverable amount of methane in storage applications.9 Breathing processes, however, begin with an open and high porosity conformation (op) at zero or low pressure of adsorbate. Upon adsorption, the structure spontaneously collapses into a less porous state (np). The deformation to a np state arises from the strong adsorption enthalpy attributed to smaller pore sizes.14 Subsequently, an increase in pressure and the amount adsorbed leads to a second transition back to the op configuration. The specific pressure of these transitions depends on the adsorbing gas and temperature.15 In many systems the op configuration is only a stable empty host configuration above a certain threshold temperature, which can result in the absence of the transitions detailed above.16 Coudert and coworkers developed a systematic description of the thermodynamics of these transitions, which has enabled the production of detailed phase diagrams and improved understand of coadsorption behavior.14,17,18 However, previous work has not considered a diverse assortment of host free energy profiles, which are representative of different SPCs.
Of particular interest in this work is the possibility for a significant energy barrier between op and np phases. This may result in the adsorption process populating metastable states and, upon transition to the np form, some of the adsorbate is actually released back to the gas phase, leading to negative gas adsorption (NGA).19 NGA phenomena were recently discovered for methane and n-butane adsorption in DUT-49.20 This has since spurred extensive experimental investigation to elucidate the effect of crystallite size and the collective nature of the process.21,22 Evans, Bocquet and Coudert also provided an initial theoretical description of NGA highlighting the ligand buckling motif that produces a large energy barrier between the op and np structures of DUT-49.23 Importantly, much of the investigation of NGA has centered on one material (DUT-49). However, in the present work we aim to demonstrate this counterintuitive phenomena in general models of porous materials.
Here we specifically develop a general, minimal theoretical model capable of reproducing responsive adsorption processes that include metastable transitions. The model is inspired by several recent studies on adsorption phenomena in flexible materials. Coudert and coworkers elaborated the thermodynamics of adsorption in flexible porous materials using the osmotic ensemble framework.24 In particular, they proposed a mean field simplification of the osmotic configurational partition function, which leads to an elegant separation of variables within the configurational integral—a simple formulation of the osmotic potential in terms of the free energy of the solid phase and the adsorption isotherm of fluid inside that phase. Ravikovitch and Neimark employed classical non-local density functional theory (NLDFT) to propose a model of elastic deformations of the solid matrix describing contraction or swelling upon adsorption.25 These ideas were later extended by Kowalczyk and coworkers to explain adsorption induced deformation in microporous carbons.26 While in this study they used Monte Carlo simulations of adsorption instead of NLDFT, a common approach begun to emerge where the flexible porous structure is treated as a slit pore with a spring attached to the walls. Several MOFs have a structure and deformation modes that can be adequately represented as a system of stacked layers separated by springs. In particular, Numaguchi and coworkers developed a theoretical description of structural transitions in ELM-11.27 Recently, Siderius and coworkers also used a combination of flat-histogram sampling Monte Carlo methods and slit pore models to study capillary phase transitions in deformable adsorbent materials.28 They proposed how these approaches can be used to construct a pore size distribution for materials with deformable structure.
In this study we build on all these previous developments to produce the general method illustrated in Fig. 1. Flexible porous materials are treated as a simple slit pore with a spring attached to its walls, and various functional forms are used to describe this spring potential. A density functional theory (DFT) description of adsorption in slit pore models provides a route for direct calculation of the grand potential of the adsorbed fluid in a pore of fixed width (a much more computationally efficient technique than Monte Carlo, for this purpose). The mean field formulation of the osmotic potential by Coudert et al. subsequently allows for the combination of grand potential from the DFT calculations with the Helmholtz free energy of the solid, represented by a potential. This simple model allows for the development of the minimum number of parameters required to formulate a new model capable of capturing the NGA effect. Furthermore, we investigate the general Helmholtz free energy profile of the host material and the external thermodynamic criteria required for responsive adsorption phenomena.
Fig. 1 Illustration of the methodology used to chart representative gas adsorption isotherms of responsive porous materials. (a) A series of isotherms a different pore widths, identified by different shades of blue. (b) The free energy profile describing the pore width. (c) A responsive isotherm calculated from the series of isotherms, as the minimum of Ωos, and identified by an orange path.
To generate a library of states (of different pore width) of the system and corresponding grand potential (Ω(T,P,w)) we use classical density functional theory (DFT) for simple fluids confined in slit pore structures.29 Slit pore models have been employed to study materials with a broad pore size distribution26 or in stacked-layer materials, such as ELM-11.27 While these models may prove difficult to provide quantitative agreement to experimental results of other porous materials they allow for the exploration of adsorption phenomena in a simplified system. Additionally, the described methodology can be readily improved by using a more complex kernel, if required.
We define the grand potential functional of the average one-body density ρ(r), eqn (2), where F[ρ(r)] is the intrinsic Helmholtz free-energy functional, Vext is the external potential and μ(ρB) is the reservoir fluid chemical potential, which is determined by the reservoir bulk density (ρB).
This methodology allows for a direct comparison of the grand potential for a given pressure and temperature at different pore widths, demonstrated in Fig. S4.† Subsequently, the complete adsorption profile for a flexible pore is generated by following the minimum of Ωos as a function of pressure.
In this study, we considered two athermal functions of Helmholtz free energy to describe the flexibility, or phases, of the slit pore systems of width w. This includes a harmonic-type potential for a system with a single minimum (eqn (3)) and a bistable-type potential (eqn (4)), which describes a system with two minima separated by barrier.
The variables (k and c) of eqn (3) are chosen to generate a potential surface centered on pore width (wmin) with a given stiffness, or looseness, defined by (k). Similarly, the coefficients (ci) of eqn (4) are selected to produce a potential surface with minima at chosen width values (wthin, wwide), relative energetic stability (Ethin, Ewide) and barrier height, measured from the lowest energetic minimum (Ebarrier).
Importantly, gas adsorption does not only follow the minimum of Ωos. We impose an additional kinetic criteria, a threshold activation energy (Ωcrit), to our model. As discussed by Numaguchi et al., structural transitions should only occur when the activation energy (or energetic barrier between the current state and the Ωos minimum state) becomes equal or less than a characteristic system energy fluctuation.27 The inclusion of this vital condition allows for the simulation of metastable transitions and also reflects the finite observation time of experimental adsorption. We set energy fluctuations of the system (Ωcrit) to be 6kT and we note this is used here as a free parameter. The ESI† includes a detailed description of Ωcrit, as it relates to these systems.
An example of the code used to generate the responsive isotherms can be found at the data repository of J. D. Evans at https://github.com/jackevansadl/supp-data. We encourage researchers to explore this methodology by combining alternative isotherm kernels and complex host free energy potentials.
This generalized theoretical framework can reproduce a rich variety of gas adsorption processes. Free energy profiles, describing the potential energy surface of the slit pore width, when combined with the energetics associated with gas adsorption, allow for flexible and responsive gas adsorption mechanisms to arise. As an initial example, we demonstrate the result of three different potential profiles dictating the behavior of the slit pore system exemplified in Fig. 2.
Fig. 2 (a) Three cases of host free energy profiles and their respective critical parameters. (b) Representative gas adsorption isotherms, where N* is the adsorbed density and λ* is configurational activity, at T* = 0.80 for each of the three cases.
Firstly, a tight harmonic potential with a stiffness k* = 10, case 1, produces a typical type-1 gas adsorption representative of a rigid porous material. During the process of adsorption, the osmotic potential does not display new minima, thus no changes in pore width can occur (Fig. S5†).
The free energy surface described as case 2 includes two minima with the lowest energy minimum at smaller width Ethin < Ewide. At low external pressure, or configurational activity λ*, the system behaves as a rigid material at pore width at wthin. When the system reaches a characteristic gate-opening pressure the system expands to wwide, allowing for increased amount of gas uptake producing a gate-opening adsorption profile. The osmotic potential during this adsorption process (Fig. S6†) clearly shows that the smaller width minimum is further stabilized at low pressure, but at high pressure the larger width phase becomes lower in energy, allowing for an np → op transition.
Case 3 defines another free energy surface with two minima. However, for this case the lowest energy minimum corresponds to the larger width, Ewide < Ethin. This potential produces a breathing-type isotherm profile where the system begins at large width, which, with increasing pressure transforms to a smaller pore state. Subsequently, at higher pressures the system will undergo a second transformation to the original pore state. These transformations op → np and np → op result from significant differences in grand potential of adsorption for different pore widths and different pressures (Fig. S7†). The reason for these differences in grand potential correspond to favorable interactions following adsorption of a single layer or multiple layers of adsorbate molecules.31–33 For example, the minima of the grand potential in pores of smaller size (below w* = 5.00) is associated with layering effects for pores completely filled with a liquid-like adsorbate.
These illustrative examples highlight thermodynamic insights into responsive gas adsorption processes. For instance, the adsorption-induced deformations of the case 3 system do not simply follow the pore-width minima defined by the host free energy profile. The osmotic potential at external pressure range λ* ≈ 1 × 10−3 shows two distinct minima at w* = 2.175 and w* = 2.875, which correspond to favorable layering effects previously discussed (Fig. S7†). Moreover, we also observe further expansion of the pore width at higher external pressure. This is because the system can accommodate more adsorbate at higher activities. In the current model this is captured by the osmotic potential achieving lower values because of smaller values of the grand potential in fully filled larger pores, even if it comes with a penalty of increasing the slit pore width (due to the host free energy potential).
The above examples are representative cases that demonstrate the ability of this model to generate interesting and important gas adsorption processes. However, it is not correct to assume that any bistable system with Ewide < Ethin (case 3) to result in NGA. To explore the result of different forms of the free energy surface we conducted parameter screening resulting in responsive adsorption maps illustrated in Fig. 3. These responsive adsorption maps provide a rich understanding of the requirements of the host pore structure to result in flexibility or NGA. The amount of flexibility here is illustrated as the sum of the absolute value of differential widths (sdw) eqn (5). The magnitude of this value will be relative to the change in width during any transition and will be even larger for breathing processes, in which there are two changes in pore width.
Fig. 3 Responsive adsorption maps for three cases of free energy surfaces at T* = 0.80, (a) harmonic potential, (b) bistable potential with Ethin = 0 and (c) bistable potential with Ewide = 0. Top row: an illustration of the free energy surface and critical parameters. Middle row: amount of adsorption-induced flexibility demonstrated by the sum of differential widths (sdw). Bottom row: density of gas released (N*) by a negative gas adsorption step, if one occurs.
The maps were produced by generating gas adsorption isotherms for a wide grid search of conditions ranging between: 0.1 ≤ k* ≤ 10 and and and for harmonic, bistable (gate-opening) and bistable (breathing) potentials, respectively.
For specific values of k* and we find that a harmonic potential can show significant flexibility. Small values of k* produce a flexible system with a broad minimum, which allows for the system of adapt to an optimal size for adsorbate multilayer formation, which as described previously is extremely favorable for slit pore systems (Fig. S4†).35 This produces a continuous breathing response (Fig. S9†), which has been reported for a number of porous materials.36,37 While under the conditions for metastability we imposed this system does not demonstrate NGA, however, the responsive adsorption exhibited by this potential may open new avenues for the exploration of materials that exhibit novel adsorption processes.
Fig. 3 also clearly demonstrates a bistable potential with Ethin < Ethin is crucial to produce NGA. NGA transitions occur only for a narrow subsection of this flexibility range. Only flexible systems with large barrier heights (Ebarrier) can permit a metastable state and produce an NGA step. Importantly, exploring the gas adsorption profiles for a wide parameter space for the other potentials did not result in NGA.
We have demonstrated the profound effects the host free energy profile can have on phase changes during adsorption including the population of metastable states. Finally, we have used the above methodology to explore the effect of increased temperature on possible NGA processes, Fig. 4.
Fig. 4 The relative frequency of systems that show changes of phase (flex.) and NGA steps (NGA) during adsorption for three temperatures compared to the temperature ranges for N2, CH4, Xe and n-butane, for which NGA is experimentally observed in DUT-49.
The above temperature range acts as a guide for the characterization of new materials to enable the discovery of new responsive behaviors in porous materials.
We have described an approach to investigate responsive gas adsorption processes using a general simulation methodology combining the osmotic potential and classical DFT of slit pore models. This has enabled the systematic and efficient simulation and mapping of gas adsorption isotherms in soft porous materials by allowing for the pore width of the system to change size during gas adsorption according to the osmotic potential.
Employing specific free energy profiles of the host system, this model can reproduce adsorption phenomena such as gate-opening, breathing and NGA. Moreover, we have conducted comprehensive scans of the critical parameters of the host free energy profiles to produce responsive adsorption maps. These maps provide crucial insight into the specific thermodynamics responsible for phase changes during adsorption and NGA.
We subsequently used this model to highlight important design and discovery guidelines for NGA materials. Firstly, bistable systems with large energetic barriers combined with large pore sizes are required for NGA. Secondly, subcritical temperature is observed to produce more unique systems that exhibit breathing and NGA.
Our described approach enables new investigations in this burgeoning research area, allowing for both the understanding of existing material properties and the imaginative design of new responsive adsorption phenomena.
This project has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement no. 742743) in addition to the Deutsche Forschungsgemeinschaft (grant agreement no. FOR2433). J. D. E acknowledges HPC platforms provided by a GENCI grant (A0050807069) and Center for Information Services and High Performance Computing (ZIH) at TU Dresden.
M. E. Davis, Nature, 2002, 417, 813–821 CrossRef CAS PubMed.
C. Perego and R. Millini, Chem. Soc. Rev., 2013, 42, 3956–3976 RSC.
D. L. Gin and R. D. Noble, Science, 2011, 332, 674–676 CrossRef CAS PubMed.
S. Horike, S. Shimomura and S. Kitagawa, Nat. Chem., 2009, 1, 695–704 CrossRef CAS PubMed.
F.-X. Coudert, Chem. Mater., 2015, 27, 1905–1916 CrossRef CAS.
T. D. Bennett, A. K. Cheetham, A. H. Fuchs and F.-X. Coudert, Nat. Chem., 2016, 9, 11–16 CrossRef PubMed.
M. D. Allendorf, R. J. T. Houk, L. Andruszkiewicz, A. A. Talin, J. Pikarsky, A. Choudhury, K. A. Gall and P. J. Hesketh, J. Am. Chem. Soc., 2008, 130, 14404–14405 CrossRef CAS PubMed.
P. Horcajada, C. Serre, G. Maurin, N. A. Ramsahye, F. Balas, M. Vallet-Regí, M. Sebban, F. Taulelle and G. Férey, J. Am. Chem. Soc., 2008, 130, 6774–6780 CrossRef CAS PubMed.
J. A. Mason, J. Oktawiec, M. K. Taylor, M. R. Hudson, J. Rodriguez, J. E. Bachman, M. I. Gonzalez, A. Cervellino, A. Guagliardi, C. M. Brown, P. L. Llewellyn, N. Masciocchi and J. R. Long, Nature, 2015, 527, 357–361 CrossRef CAS PubMed.
L. Vanduyfhuys, S. M. J. Rogge, J. Wieme, S. Vandenbrande, G. Maurin, M. Waroquier and V. Van Speybroeck, Nat. Commun., 2018, 9, 204 CrossRef CAS PubMed.
A. Schneemann, V. Bon, I. Schwedler, I. Senkovska, S. Kaskel and R. A. Fischer, Chem. Soc. Rev., 2014, 43, 6062–6096 RSC.
H. Furukawa, K. E. Cordova, M. O'Keeffe and O. M. Yaghi, Science, 2013, 341, 1230444 CrossRef PubMed.
M. Pera-Titus and D. Farrusseng, J. Phys. Chem. C, 2012, 116, 1638–1649 CrossRef CAS.
A. Boutin, M.-A. Springuel-Huet, A. Nossov, A. Gédéon, T. Loiseau, C. Volkringer, G. Férey, F.-X. Coudert and A. Fuchs, Angew. Chem., Int. Ed., 2009, 48, 8314–8317 CrossRef CAS PubMed.
J. Wieme, K. Lejaeghere, G. Kresse and V. Van Speybroeck, Nat. Commun., 2018, 9, 4899 CrossRef CAS PubMed.
A. Boutin, F.-X. Coudert, M.-A. Springuel-Huet, A. V. Neimark, G. Férey and A. H. Fuchs, J. Phys. Chem. C, 2010, 114, 22237–22244 CrossRef CAS.
F.-X. Coudert, A. Boutin, M. Jeffroy, C. Mellot- Draznieks and A. H. Fuchs, ChemPhysChem, 2011, 12, 247–258 CrossRef CAS PubMed.
A. U. Ortiz, M.-A. Springuel-Huet, F.-X. Coudert, A. H. Fuchs and A. Boutin, Langmuir, 2011, 28, 494–498 CrossRef PubMed.
S. Krause, V. Bon, I. Senkovska, U. Stoeck, D. Wallacher, D. M. Többens, S. Zander, R. S. Pillai, G. Maurin, F.-X. Coudert and S. Kaskel, Nature, 2016, 532, 348–352 CrossRef CAS PubMed.
U. Stoeck, S. Krause, V. Bon, I. Senkovska and S. Kaskel, Chem. Commun., 2012, 48, 10841 RSC.
S. Krause, V. Bon, I. Senkovska, D. M. Többens, D. Wallacher, R. S. Pillai, G. Maurin and S. Kaskel, Nat. Commun., 2018, 9, 1573 CrossRef PubMed.
J. Schaber, S. Krause, S. Paasch, I. Senkovska, V. Bon, D. M. Többens, D. Wallacher, S. Kaskel and E. Brunner, J. Phys. Chem. C, 2017, 121, 5195–5200 CrossRef CAS.
J. D. Evans, L. Bocquet and F.-X. Coudert, Chem, 2016, 1, 873–886 CAS.
F.-X. Coudert, M. Jeffroy, A. H. Fuchs, A. Boutin and C. Mellot-Draznieks, J. Am. Chem. Soc., 2008, 130, 14294–14302 CrossRef CAS PubMed.
P. I. Ravikovitch and A. V. Neimark, Langmuir, 2006, 22, 10864–10868 CrossRef CAS PubMed.
P. Kowalczyk, A. Ciach and A. V. Neimark, Langmuir, 2008, 24, 6603–6608 CrossRef CAS PubMed.
R. Numaguchi, H. Tanaka, S. Watanabe and M. T. Miyahara, J. Chem. Phys., 2013, 138, 054708 CrossRef PubMed.
V. K. Shen and D. W. Siderius, J. Chem. Phys., 2014, 140, 244106 CrossRef PubMed.
M. B. Sweatman and N. Quirke, Langmuir, 2002, 18, 10443–10454 CrossRef CAS.
M. B. Sweatman, Mol. Phys., 2000, 98, 573–581 CrossRef CAS.
P. B. Balbuena and K. E. Gubbins, Langmuir, 1993, 9, 1801–1814 CrossRef CAS.
C. Lastoskie, K. E. Gubbins and N. Quirke, J. Phys. Chem., 1993, 97, 4786–4796 CrossRef CAS.
C. Lastoskie, K. E. Gubbins and N. Quirke, Langmuir, 1993, 9, 2693–2702 CrossRef CAS.
G. Manos and L. Dunne, Nanomaterials, 2018, 8, 818 CrossRef PubMed.
T. X. Nguyen, S. K. Bhatia and D. Nicholson, J. Chem. Phys., 2002, 117, 10827–10836 CrossRef CAS.
H. J. Choi, M. Dincă and J. R. Long, J. Am. Chem. Soc., 2008, 130, 7848–7850 CrossRef CAS PubMed.
E. J. Carrington, C. A. McAnally, A. J. Fletcher, S. P. Thompson, M. Warren and L. Brammer, Nat. Chem., 2017, 9, 882–889 CrossRef CAS PubMed.
F. Nouar, T. Devic, H. Chevreau, N. Guillou, E. Gibson, G. Clet, M. Daturi, A. Vimont, J. M. Grenèche, M. I. Breeze, R. I. Walton, P. L. Llewellyn and C. Serre, Chem. Commun., 2012, 48, 10237 RSC.
G. Férey and C. Serre, Chem. Soc. Rev., 2009, 38, 1380 RSC.

References: V. 
 V. 
 V. 
 V. 
 V. 
 V. 
 V. 
 V. 
 V. 
 V. 

V.