Multiscale Modeling of Physical Properties of Nanoporous Frameworks: Predicting Mechanical, Thermal, and Adsorption Behavior.

ConspectusNanoporous frameworks are a large and diverse family of supramolecular materials, whose chemical building units (organic, inorganic, or both) are assembled into a 3D architecture with well-defined connectivity and topology, featuring intrinsic porosity. These materials play a key role in various industrial processes and applications, such as energy production and conversion, fluid separation, gas storage, water harvesting, and many more. The performance and suitability of nanoporous materials for each specific application are directly related to both their physical and chemical properties, and their determination is crucial for process engineering and optimization of performances. In this Account, we focus on some recent developments in the multiscale modeling of physical properties of nanoporous frameworks, highlighting the latest advances in three specific areas: mechanical properties, thermal properties, and adsorption.In the study of the mechanical behavior of nanoporous materials, the past few years have seen a rapid acceleration of research. For example, computational resources have been pooled to create a public large-scale database of elastic constants as part of the Materials Project initiative to accelerate innovation in materials research: those can serve as a basis for data-based discovery of materials with targeted properties, as well as the training of machine learning predictor models.The large-scale prediction of thermal behavior, in comparison, is not yet routinely performed at such a large scale. Tentative databases have been assembled at the DFT level on specific families of materials, such as zeolites, but prediction at larger scale currently requires the use of transferable classical force fields, whose accuracy can be limited.Finally, adsorption is naturally one of the most studied physical properties of nanoporous frameworks, as fluid separation or storage is often the primary target for these materials. We highlight the recent achievements and open challenges for adsorption prediction at a large scale, focusing in particular on the accuracy of computational models and the reliability of comparisons with experimental data available. We detail some recent methodological improvements in the prediction of adsorption-related properties: in particular, we describe the recent research efforts to go beyond the study of thermodynamic quantities (uptake, adsorption enthalpy, and thermodynamic selectivity) and predict transport properties using data-based methods and high-throughput computational schemes. Finally, we stress the importance of data-based methods of addressing all sources of uncertainty.The Account concludes with some perspectives about the latest developments and open questions in data-based approaches and the integration of computational and experimental data together in the materials discovery loop.

Molecular Simulation of the Impact of Defects on Electrolyte Intrusion in Zeolites.

We have investigated through molecular simulation the intrusion of electrolytes in two representative pure-silica zeolites, silicalite-1 and chabazite, in which point defects were introduced in varying amounts. We distinguish between two types of defects, considering either "weak" or "strong" silanol nest defects, resulting in different hydration behaviors. In the presence of weak defects, the hydration process occurs through a homogeneous nucleation process, while with strong defects, we observe an initial adsorption followed by a filling of the nanoporous volume at a higher pressure. However, we show that electrolytes do not penetrate the zeolites, and these defects appear to have only marginal influence on the thermodynamics of electrolyte intrusion. While replacing pure water by the electrolyte solution shifts the intrusion pressure toward higher values because of the drop of water saturation vapor pressure, an increase in hydrophilicity of the framework due to point defects has the opposite effect, showing that controlling the amount of defects in zeolites is crucial for storage energy applications.

Defective Nature of CdSe Quantum Dots Embedded in Inorganic Matrices.

Quantum dots (QDs) embedded in inorganic matrices have been extensively studied for their potential applications in lighting, displays, and solar cells. While a significant amount of research studies focused on their experimental fabrication, the origin of their relatively low photoluminescence quantum yield has not been investigated yet, although it severely hinders practical applications. In this study, we use time-dependent density functional theory (TDDFT) to pinpoint the nature of excited states of CdSe QDs embedded in various inorganic matrices. The formation of undercoordinated Se atoms and nonbridging oxygen atoms at the QD/glass interface is responsible for the localization of a hole wave function, leading to the formation of low-energy excited states with weak oscillator strength. These states provide pathways for nonradiative processes and compete with radiative emission. The photoluminescence performance is predicted for CdSe QDs in different matrices and validated by experiments. The results of this study have significant implications for understanding the underlying photophysics of CdSe QDs embedded in inorganic matrices that would facilitate the fabrication of highly luminescent glasses.

Chiral Lanthanum Metal-Organic Framework with Gated CO2 Sorption and Concerted Framework Flexibility.

A metal-organic framework (MOF) CTH-17 based on lanthanum(III) and the conformationally chiral linker 1,2,3,4,5,6-hexakis(4-carboxyphenyl)benzene, cpb6-: [La2(cpb)]·1.5dmf was prepared by the solvothermal method in dimethylformamide (dmf) and characterized by variable-temperature X-ray powder diffraction (VTPXRD), variable-temperature X-ray single-crystal diffraction (SCXRD), and thermogravimetric analysis (TGA). CTH-17 is a rod-MOF with new topology och. It has high-temperature stability with Sohncke space groups P6122/P6522 at 90 K and P622 at 300 and 500 K, all phases characterized with SCXRD and at 293 K also with three-dimensional (3D) electron diffraction. VTPXRD indicates a third phase appearing after 620 K and stable up to 770 K. Gas sorption isotherms with N2 indicate a modest surface area of 231 m2 g-1 for CTH-17, roughly in agreement with the crystal structure. Carbon dioxide sorption reveals a gate-opening effect of CTH-17 where the structure opens up when the loading of CO2 reaches approximately ∼0.45 mmol g-1 or 1 molecule per unit cell. Based on the SCXRD data, this is interpreted as flexibility based on the concerted movements of the propeller-like hexatopic cpb linkers, the movement intramolecularly transmitted by the π-π stacking of the cpb linkers and helped by the fluidity of the LaO6 coordination sphere. This was corroborated by density functional theory (DFT) calculations yielding the chiral phase (P622) as the energy minimum and a completely racemic phase (P6/mmm), with symmetric cpb linkers representing a saddle point in a racemization process.

The changing state of porous materials.

Porous materials contain regions of empty space into which guest molecules can be selectively adsorbed and sometimes chemically transformed. This has made them useful in both industrial and domestic applications, ranging from gas separation, energy storage and ion exchange to heterogeneous catalysis and green chemistry. Porous materials are often ordered (crystalline) solids. Order-or uniformity-is frequently held to be advantageous, or even pivotal, to our ability to engineer useful properties in a rational way. Here we highlight the growing evidence that topological disorder can be useful in creating alternative properties in porous materials. In particular, we highlight here several concepts for the creation of novel porous liquids, rationalize routes to porous glasses and provide perspectives on applications for porous liquids and glasses.

A pressure-amplifying framework material with negative gas adsorption transitions.

Adsorption-based phenomena are important in gas separations, such as the treatment of greenhouse-gas and toxic-gas pollutants, and in water-adsorption-based heat pumps for solar cooling systems. The ability to tune the pore size, shape and functionality of crystalline porous coordination polymers--or metal-organic frameworks (MOFs)--has made them attractive materials for such adsorption-based applications. The flexibility and guest-molecule-dependent response of MOFs give rise to unexpected and often desirable adsorption phenomena. Common to all isothermal gas adsorption phenomena, however, is increased gas uptake with increased pressure. Here we report adsorption transitions in the isotherms of a MOF (DUT-49) that exhibits a negative gas adsorption; that is, spontaneous desorption of gas (methane and n-butane) occurs during pressure increase in a defined temperature and pressure range. A combination of in situ powder X-ray diffraction, gas adsorption experiments and simulations shows that this adsorption behaviour is controlled by a sudden hysteretic structural deformation and pore contraction of the MOF, which releases guest molecules. These findings may enable technologies using frameworks capable of negative gas adsorption for pressure amplification in micro- and macroscopic system engineering. Negative gas adsorption extends the series of counterintuitive phenomena such as negative thermal expansion and negative refractive indices and may be interpreted as an adsorptive analogue of force-amplifying negative compressibility transitions proposed for metamaterials.

Emergence of Coupled Rotor Dynamics in Metal-Organic Frameworks via Tuned Steric Interactions.

The organic components in metal-organic frameworks (MOFs) are unique: they are embedded in a crystalline lattice, yet, as they are separated from each other by tunable free space, a large variety of dynamic behavior can emerge. These rotational dynamics of the organic linkers are especially important due to their influence over properties such as gas adsorption and kinetics of guest release. To fully exploit linker rotation, such as in the form of molecular machines, it is necessary to engineer correlated linker dynamics to achieve their cooperative functional motion. Here, we show that for MIL-53, a topology with closely spaced rotors, the phenylene functionalization allows researchers to tune the rotors' steric environment, shifting linker rotation from completely static to rapid motions at frequencies above 100 MHz. For steric interactions that start to inhibit independent rotor motion, we identify for the first time the emergence of coupled rotation modes in linker dynamics. These findings pave the way for function-specific engineering of gear-like cooperative motion in MOFs.

Transient Catenation in a Zirconium-Based Metal-Organic Framework and Its Effect on Mechanical Stability and Sorption Properties.

Interpenetration of two or more sublattices is common among many metal-organic frameworks (MOFs). Herein, we study the evolution of one zirconium cluster-based, 3,8-connected MOF from its non-interpenetrated (NU-1200) to interpenetrated (STA-26) isomer. We observe this transient catenation process indirectly using ensemble methods, such as nitrogen porosimetry and X-ray diffraction, and directly, using high-resolution transmission electron microscopy. The approach detailed here will serve as a template for other researchers to monitor the interpenetration of their MOF samples at the bulk and single-particle limits. We investigate the mechanical stability of both lattices experimentally by pressurized in situ X-ray diffraction and nanoindentation as well as computationally with density functional theory calculations. Both lines of study reveal that STA-26 is considerably more mechanically stable than NU-1200. We conclude this study by demonstrating the potential of these MOFs and their mixed phases for the capture of gaseous n-hexane, used as a structural mimic for the chemical warfare agent sulfur mustard gas.

Water Adsorption in Soft and Heterogeneous Nanopores.

ConspectusLiquids under confinement differ in behavior from their bulk counterparts and can acquire properties that are specific to the confined phase and linked to the nature and structure of the host matrix. While confined liquid water is not a new topic of research, the past few years have seen a series of intriguing novel features for water inside nanoscale pores. These unusual properties arise from the very specific nature of nanoporous materials, termed "soft porous crystals"; they combine large-scale flexibility with a heterogeneous internal surface. This creates a rich diversity of behavior for the adsorbed water, and the combination of different experimental characterization techniques along with computational chemistry at various scales is necessary to understand the phenomena observed and their microscopic origins. The range of systems of interest span the whole chemical range, from the inorganic (zeolites, imogolites) to the organic (microporous carbons, graphene, and its derivatives), and even encompass the hybrid organic-inorganic systems (such as metal-organic frameworks).The combination of large scale flexibility with the strong physisorption (or even chemisorption) of water can lead to unusual properties (belonging to the "metamaterials" category) and to novel phenomena. One striking example is the recent elucidation of the mechanism of negative hydration expansion in ZrW2O8, by which adsorption of ∼10 wt % water in the inorganic nonporous framework leads to large shrinkage of its volume. Another eye-catching case is the occurrence of multiple water adsorption-driven structural transitions in the MIL-53 family of materials: the specific interactions between water guest molecules and the host framework create behavior that has not been observed with any other adsorbate. Both are counterintuitive phenomena that have been elucidated by a combination of experimental in situ techniques and molecular simulation.Another important direction of research is the shift in the systems and phenomena studied, from physical adsorption toward studies of reactivity, hydrothermal stability, and the effect of confinement on aqueous phases more complex than pure water. There have been examples of water adsorption in highly flexible metal-organic frameworks being able to compete with the materials' coordination bonds, thereby limiting its hydrothermal stability, while tweaking the functional groups of the same framework can lead to increased stability while retaining the flexibility of the material. However, this additional complexity and tunability in the macroscopic behavior can occur from changes in the confined fluid rather than the material. Very recent studies have shown that aqueous solutions of high concentration (such as LiCl up to 20 mol L-1) confined in flexible nanoporous materials can have specific properties different from pure water and not entirely explained by osmotic effects. There, the strong ordering of the confined electrolyte competes with the structural flexibility of the framework to create an entirely new behavior for the {host, guest} system.

Thermodynamic exploration of xenon/krypton separation based on a high-throughput screening.

Nanoporous framework materials are a promising class of materials for energy-efficient technology of xenon/krypton separation by physisorption. Many studies on Xe/Kr separation by adsorption have focused on the determination of structure/property relationships, the description of theoretical limits of performance, and the identification of top-performing materials. Here, we provide a study based on a high-throughput screening of the adsorption of Xe, Kr, and Xe/Kr mixtures in 12 020 experimental MOF materials, to provide a better comprehension of the thermodynamics behind Xe/Kr separation in nanoporous materials and the microscopic origins of Xe/Kr selectivity at both low and ambient pressure.

The role of temperature and adsorbate on negative gas adsorption transitions of the mesoporous metal-organic framework DUT-49.

Unusual adsorption phenomena, such as breathing and negative gas adsorption (NGA), are rare and challenge our thermodynamic understanding of adsorption in deformable porous solids. In particular, NGA appears to break the rules of thermodynamics in these materials by exhibiting a spontaneous release of gas accompanying an increase in pressure. This anomaly relies on long-lived metastable states. A fundamental understanding of this process is desperately required for the discovery of new materials with this exotic property. Interestingly, NGA was initially observed upon adsorption of methane at relatively low temperature, close to the respective standard boiling point of the adsorptive, and no NGA was observed at elevated temperatures. In this contribution, we present an extensive investigation of adsorption of an array of gases at various temperatures on DUT-49, a material which features an NGA transition. Experiments, featuring a wide range of gases and vapors at temperatures ranging from 21-308 K, were used to identify for each guest a critical temperature range in which NGA can be detected. The experimental results were complemented by molecular simulations that help to rationalize the absence of NGA at elevated temperatures, and the non-monotonic behavior present upon temperature decrease. Furthermore, this in-depth analysis highlights the crucial thermodynamic and kinetic conditions for NGA, which are unique to each guest and potentially other solids with similar effects. We expect this exploration to provide detailed guidelines for experimentally discovering NGA and related "rule breaking" phenomena in novel and already known materials, and provide the conditions required for the application of this effect, for example as pressure amplifying materials.

Ab Initio Molecular Dynamics of CdSe Quantum-Dot-Doped Glasses.

We have probed the local atomic structure of the interface between a CdSe quantum dot (QD) and a sodium silicate glass matrix. Using ab initio molecular dynamics simulations, we determined the structural properties and bond lengths, in excellent agreement with previous experimental observations. On the basis of an analysis of radial distribution functions, coordination environment, and ring structures, we demonstrate that an important structural reconstruction occurs at the interface between the CdSe QD and the glass matrix. The incorporation of the CdSe QD disrupts the Na-O bonds, while stronger SiO4 tetrahedra are reformed. The existence of the glass matrix breaks the stable 4-membered (4MR) and 6-membered (6MR) Cd-Se rings, and we observe a disassociated Cd atom migrated in the glass matrix. Besides, the formation of Se-Na and Cd-O linkages is observed at the CdSe QD/glass interface. These results significantly extend our understanding of the interfacial structure of CdSe QD-doped glasses and provide physical and chemical insight into the possible defect structure origin of CdSe QD, of interest to the fabrication of the highly luminescent CdSe QD-doped glasses.

Pressure promoted low-temperature melting of metal-organic frameworks.

Metal-organic frameworks (MOFs) are microporous materials with huge potential for chemical processes. Structural collapse at high pressure, and transitions to liquid states at high temperature, have recently been observed in the zeolitic imidazolate framework (ZIF) family of MOFs. Here, we show that simultaneous high-pressure and high-temperature conditions result in complex behaviour in ZIF-62 and ZIF-4, with distinct high- and low-density amorphous phases occurring over different regions of the pressure-temperature phase diagram. In situ powder X-ray diffraction, Raman spectroscopy and optical microscopy reveal that the stability of the liquid MOF state expands substantially towards lower temperatures at intermediate, industrially achievable pressures and first-principles molecular dynamics show that softening of the framework coordination with pressure makes melting thermodynamically easier. Furthermore, the MOF glass formed by melt quenching the high-temperature liquid possesses permanent, accessible porosity. Our results thus imply a route to the synthesis of functional MOF glasses at low temperatures, avoiding decomposition on heating at ambient pressure.

Mixed-metal metal-organic frameworks.

Mixed-metal MOFs are metal-organic frameworks that contain at least 2 different metal ions as nodes of their frameworks. They are prepared relatively easily by either a one-pot synthesis with a synthesis mixture containing the different metals, or by a post-synthetic ion-exchange method by soaking a monometallic MOF in a concentrated solution of a different (but compatible) metal-ion. More difficult is the accurate characterization of these materials. Is the formed product a mixture of monometallic MOFs or indeed a MOF with different metallic nodes? Are the metals randomly distributed or do they form domains? What is the oxidation state of the metals? How do the metals mutually influence each other, and impact the material's performance? Advanced characterization techniques are required e.g. X-ray absorption spectroscopy, magnetic resonance and electron microscopy. Computational tools at multiple scales are also often applied. In almost every case, a judicious choice of several techniques is required to unambiguously characterize the mixed-metal MOF. Although still in their infancy, several applications are emerging for mixed-metal MOFs, that improve on conventional monometallic MOFs. In the field of gas sorption and storage, especially the stability and affinity towards the target gases can be largely improved by introducing a second metal ion. In the case of flexible MOFs, the breathing behavior, and in particular the pressure at which the MOF opens, can be tailored. In heterogeneous catalysis, new cascade and tandem reactions become possible, with particular focus on reactions where the two metals in close proximity truly form a mixed-metal transition state. The bimetallic MOF should have a clear benefit over a mixture of the respective monometallic MOFs, and bimetallic enzymes can be a huge source of inspiration in this field. Another very promising application lies in the fields of luminescence and sensing. By tuning the lanthanide metals in mixed-metal lanthanide MOFs and by using the organic linkers as antennae, novel smart materials can be developed, acting as sensors and as thermochromic thermometers. Of course there are also still open challenges, as also mixed-metal MOFs do not escape the typical drawbacks of MOFs, such as low stability in moisture and possible metal leaching in liquids. The ease of synthesis of mixed-metal MOFs is a large bonus. In this critical review, we discuss in detail the synthesis, characterization, computational work and applications of mixed-metal MOFs.

Rich Polymorphism of a Metal-Organic Framework in Pressure-Temperature Space.

We present an in situ powder X-ray diffraction study on the phase stability and polymorphism of the metal-organic framework ZIF-4, Zn(imidazolate)2, at simultaneous high pressure and high temperature, up to 8 GPa and 600 °C. The resulting pressure-temperature phase diagram reveals four, previously unknown, high-pressure-high-temperature ZIF phases. The crystal structures of two new phases-ZIF-4-cp-II and ZIF-hPT-II-were solved by powder diffraction methods. The total energy of ZIF-4-cp-II was evaluated using density functional theory calculations and was found to lie in between that of ZIF-4 and the most thermodynamically stable polymorph, ZIF- zni. ZIF-hPT-II was found to possess a doubly interpenetrated diamondoid topology and is isostructural with previously reported Cd(Imidazolate)2 and Hg(Imidazolate)2 phases. This phase exhibited extreme resistance to both temperature and pressure. The other two new phases could be assigned with a unit cell and space group, although their structures remain unknown. The pressure-temperature phase diagram of ZIF-4 is strikingly complicated when compared with that of the previously investigated, closely related ZIF-62 and demonstrates the ability to traverse complex energy landscapes of metal-organic systems using the combined application of pressure and temperature.

MOF Decomposition and Introduction of Repairable Defects Using a Photodegradable Strut.

Photoswitchable components can modulate the properties of metal organic frameworks (MOFs); however, photolabile building blocks remain underexplored. A new strut NPDAC (2-nitro-1,4-phenylenediacetic acid) that undergoes photodecarboxylation has been prepared and incorporated into a MOF, using post-synthetic linker exchange (PSLE) from the structural analogue containing PDAC (p-phenylenediacetic acid). Irradiation of NPDAC-MOF leads to MOF decomposition and concomitant formation of amorphous material. In addition to complete linker exchange, MOFs containing a mixture of PDAC and NPDAC can be obtained through partial linker exchange. In NPDAC30-MOF, which contains approximately 30 % NPDAC, the MOF retains crystallinity after irradiation, but the MOF contains defect sites consistent with loss of decarboxylated NPDAC linkers. The defect sites can be repaired by exposure to additional PDAC or NPDAC linkers at a much faster rate than the initial exchange process. The photoremoval and replacement process may lead to a more general approach to customizable MOF structures.

Modelling of framework materials at multiple scales: current practices and open questions.

The last decade has seen an explosion of the family of framework materials and their study, from both the experimental and computational points of view. We propose here a short highlight of the current state of methodologies for modelling framework materials at multiple scales, putting together a brief review of new methods and recent endeavours in this area, as well as outlining some of the open challenges in this field. We will detail advances in atomistic simulation methods, the development of material databases and the growing use of machine learning for the prediction of properties. This article is part of the theme issue 'Mineralomimesis: natural and synthetic frameworks in science and technology'.

Liquid metal-organic frameworks.

Metal-organic frameworks (MOFs) are a family of chemically diverse materials, with applications in a wide range of fields, covering engineering, physics, chemistry, biology and medicine. Until recently, research has focused almost entirely on crystalline structures, yet now a clear trend is emerging, shifting the emphasis onto disordered states, including 'defective by design' crystals, as well as amorphous phases such as glasses and gels. Here we introduce a strongly associated MOF liquid, obtained by melting a zeolitic imidazolate framework. We combine in situ variable temperature X-ray, ex situ neutron pair distribution function experiments, and first-principles molecular dynamics simulations to study the melting phenomenon and the nature of the liquid obtained. We demonstrate from structural, dynamical, and thermodynamical information that the chemical configuration, coordinative bonding, and porosity of the parent crystalline framework survive upon formation of the MOF liquid.

Negative Hydration Expansion in ZrW_{2}O_{8}: Microscopic Mechanism, Spaghetti Dynamics, and Negative Thermal Expansion.

We use a combination of x-ray diffraction, total scattering, and quantum mechanical calculations to determine the mechanism responsible for hydration-driven contraction in ZrW_{2}O_{8}. The inclusion of H_{2}O molecules within the ZrW_{2}O_{8} network drives the concerted formation of new W─O bonds to give one-dimensional (─W─O─)_{n} strings. The topology of the ZrW_{2}O_{8} network is such that there is no unique choice for the string trajectories: the same local changes in coordination can propagate with a large number of different periodicities. Consequently, ZrW_{2}O_{8}·H_{2}O is heavily disordered, with each configuration of strings forming a dense aperiodic "spaghetti." This new connectivity contracts the unit cell via large shifts in the Zr and W atom positions. Fluctuations of the undistorted parent structure towards this spaghetti phase emerge as the key negative thermal expansion (NTE) phonon modes in ZrW_{2}O_{8} itself. The large relative density of NTE phonon modes in ZrW_{2}O_{8} actually reflects the degeneracy of volume-contracting spaghetti excitations, itself a function of the particular topology of this remarkable material.

Structure and Dynamics of Water Confined in Imogolite Nanotubes.

We have studied the properties of water adsorbed inside nanotubes of hydrophilic imogolite, an aluminum silicate clay mineral, by means of molecular simulations. We used a classical force field to describe the water and the flexible imogolite nanotube and validated it against the data obtained from first-principles molecular dynamics. With it, we observe a strong structuration of the water confined in the nanotube, with specific adsorption sites and a distribution of hydrogen bond patterns. The combination of number of adsorption sites, their geometry, and the preferential tetrahedral hydrogen bonding pattern of water leads to frustration and disorder. We further characterize the dynamics of the water, as well as the hydrogen bonds formed between water molecules and the nanotube, which is found to be more than 1 order of magnitude longer than water-water hydrogen bonds.