Modeling and Simulation of Crystallization of Metal–Organic Frameworks

The use of collective variables and enhanced sampling in the simulations of existing and emerging microporous materials

Microporous materials, including zeolites, metal-organic frameworks, and cage compounds, offer diverse functionalities due to their unique dynamics and guest confinement properties. These materials play a significant role in separation, catalysis, and sensing, but their complexity hinders exploration using traditional atomistic simulations. This review explores collective variables (CVs) paired with enhanced sampling as a powerful approach to enable efficient investigation of key features in microporous materials. We highlight successful applications of CVs in studying adsorption, diffusion, phase transitions, and mechanical properties, demonstrating their crucial role in guiding material design and optimisation. The future of CVs lies in integration with techniques like machine learning, allowing for enhanced efficiency and accuracy. By tailoring CVs to specific materials and developing multi-scale approaches we can further unlock the intricacies of these fascinating materials. Simulations are a cornerstone in unravelling the complexities of microporous materials and are crucial for our future understanding.

Microkinetic insights into the role of catalyst and water activity on the nucleation, growth, and dissolution during COF-5 synthesis

The chemical pathway for synthesizing covalent organic frameworks (COFs) involves a complex medley of reaction sequences over a rippling energy landscape that cannot be adequately described using existing theories. Even with the development of state-of-the-art experimental and computational tools, identifying primary mechanisms of nucleation and growth of COFs remains elusive. Other than empirically, little is known about how the catalyst composition and water activity affect the kinetics of the reaction pathway. Here, for the first time, we employ time-resolved in situ Fourier transform infrared spectroscopy (FT-IR) coupled with a six-parameter microkinetic model consisting of ∼10 million reactions and over 20 000 species. The integrated approach elucidates previously unrecognized roles of catalyst pK_a on COF yield and water on growth rate and size distribution. COF crystalline yield increases with decreasing pK_a of the catalysts, whereas the effect of water is to reduce the growth rate of COF and broaden the size distribution. The microkinetic model reproduces the experimental data and quantitatively predicts the role of synthesis conditions such as temperature, catalyst, and precursor concentration on the nucleation and growth rates. Furthermore, the model also validates the second-order reaction mechanism of COF-5 and predicts the activation barriers for classical and non-classical growth of COF-5 crystals. The microkinetic model developed here is generalizable to different COFs and other multicomponent systems.

Brightly Luminescent (TbxLu1-x)2bdc3·nH2O MOFs: Effect of Synthesis Conditions on Structure and Luminescent Properties

Luminescent, heterometallic terbium(III)-lutetium(III) terephthalate metal-organic frameworks (MOFs) were synthesized via direct reaction between aqueous solutions of disodium terephthalate and nitrates of corresponding lanthanides by using two methods: synthesis from diluted and concentrated solutions. For (Tb_xLu_1-x)₂bdc₃·nH₂O MOFs (bdc = 1,4-benzenedicarboxylate) containing more than 30 at. % of Tb³⁺, only one crystalline phase was formed: Ln₂bdc₃·4H₂O. At lower Tb³⁺ concentrations, MOFs crystallized as the mixture of Ln₂bdc₃·4H₂O and Ln₂bdc₃·10H₂O (diluted solutions) or Ln₂bdc₃ (concentrated solutions). All synthesized samples that contained Tb³⁺ ions demonstrated bright green luminescence upon excitation into the ¹ππ* excited state of terephthalate ions. The photoluminescence quantum yields (PLQY) of the compounds corresponding to the Ln₂bdc₃ crystalline phase were significantly larger than for Ln₂bdc₃·4H₂O and Ln₂bdc₃·10H₂O phases due to absence of quenching from water molecules possessing high-energy O-H vibrational modes. One of the synthesized materials, namely, (Tb_0.1Lu_0.9)₂bdc₃·1.4H₂O, had one of the highest PLQY among Tb-based MOFs, 95%.

Autocatalysis and Oriented Attachment Direct the Synthesis of a Metal-Organic Framework

Synthesis of porous, covalent crystals such as zeolites and metal-organic frameworks (MOFs) cannot be described adequately using existing crystallization theories. Even with the development of state-of-the-art experimental and computational tools, the identification of primary mechanisms of nucleation and growth of MOFs remains elusive. Here, using time-resolved in-situ X-ray scattering coupled with a six-parameter microkinetic model consisting of ∼1 billion reactions and up to ∼100 000 metal nodes, we identify autocatalysis and oriented attachment as previously unrecognized mechanisms of nucleation and growth of the MOF UiO-66. The secondary building unit (SBU) formation follows an autocatalytic initiation reaction driven by a self-templating mechanism. The induction time of MOF nucleation is determined by the relative rate of SBU attachment (chain extension) and the initiation reaction, whereas the MOF growth is primarily driven by the oriented attachment of reactive MOF crystals. The average size and polydispersity of MOFs are controlled by surface stabilization. Finally, the microkinetic model developed here is generalizable to different MOFs and other multicomponent systems.

Ultrasound-Assisted Synthesis of Luminescent Micro- and Nanocrystalline Eu-Based MOFs as Luminescent Probes for Heavy Metal Ions

The luminescent coarse-, micro- and nanocrystalline europium(III) terephthalate tetrahydrate (Eu₂bdc₃·4H₂O) metal-organic frameworks were synthesized by the ultrasound-assisted wet-chemical method. Electron micrographs show that the europium(III) terephthalate microparticles are 7 μm long leaf-like plates. According to the dynamic light scattering technique, the average size of the Eu₂bdc₃·4H₂O nanoparticles is equal to about 8 ± 2 nm. Thereby, the reported Eu₂bdc₃·4H₂O nanoparticles are the smallest nanosized rare-earth-based MOF crystals, to the best of our knowledge. The synthesized materials demonstrate red emission due to the ⁵D₀–⁷F_J transitions of Eu³⁺ upon 250 nm excitation into ¹ππ* state of the terephthalate ion. Size reduction results in broadened emission bands, an increase in the non-radiative rate constants and a decrease in both the quantum efficiency of the ⁵D₀ level and Eu³⁺ and the luminescence quantum yields. Cu²⁺, Cr³⁺, and Fe³⁺ ions efficiently and selectively quench the luminescence of nanocrystalline europium(III) terephthalate, which makes it a prospective material for luminescent probes to monitor these ions in waste and drinking water.

Modeling the Layer-by-Layer Growth of HKUST-1 Metal-Organic Framework Thin Films

Metal organic frameworks have emerged as an important new class of materials with many applications, such as sensing, gas separation, drug delivery. In many cases, their performance is limited by structural defects, including vacancies and domain boundaries. In the case of MOF thin films, surface roughness can also have a pronounced influence on MOF-based device properties. Presently, there is little systematic knowledge about optimal growth conditions with regard to optimal morphologies for specific applications. In this work, we simulate the layer-by-layer (LbL) growth of the HKUST-1 MOF as a function of temperature and reactant concentration using a coarse-grained model that permits detailed insights into the growth mechanism. This model helps to understand the morphological features of HKUST-1 grown under different conditions and can be used to predict and optimize the temperature for the purpose of controlling the crystal quality and yield. It was found that reactant concentration affects the mass deposition rate, while its effect on the crystallinity of the generated HKUST-1 film is less pronounced. In addition, the effect of temperature on the surface roughness of the film can be divided into three regimes. Temperatures in the range from 10 to 129 °C allow better control of surface roughness and film thickness, while film growth in the range of 129 to 182 °C is characterized by a lower mass deposition rate per cycle and rougher surfaces. Finally, for T larger than 182 °C, the film grows slower, but in a smooth fashion. Furthermore, the potential effect of temperature on the crystallinity of LbL-grown HKUST-1 was quantified. To obtain high crystallinity, the operating temperature should preferably not exceed 57 °C, with an optimum around 28 °C, which agrees with experimental observations.

Current Research Trends and Perspectives on Solid-State Nanomaterials in Hydrogen Storage

Hydrogen energy, with environment amicable, renewable, efficiency, and cost-effective advantages, is the future mainstream substitution of fossil-based fuel. However, the extremely low volumetric density gives rise to the main challenge in hydrogen storage, and therefore, exploring effective storage techniques is key hurdles that need to be crossed to accomplish the sustainable hydrogen economy. Hydrogen physically or chemically stored into nanomaterials in the solid-state is a desirable prospect for effective large-scale hydrogen storage, which has exhibited great potentials for applications in both reversible onboard storage and regenerable off-board storage applications. Its attractive points include safe, compact, light, reversibility, and efficiently produce sufficient pure hydrogen fuel under the mild condition. This review comprehensively gathers the state-of-art solid-state hydrogen storage technologies using nanostructured materials, involving nanoporous carbon materials, metal-organic frameworks, covalent organic frameworks, porous aromatic frameworks, nanoporous organic polymers, and nanoscale hydrides. It describes significant advances achieved so far, and main barriers need to be surmounted to approach practical applications, as well as offers a perspective for sustainable energy research.

Solvent fluctuations in the solvation shell determine the activation barrier for crystal growth rates

Solution crystallization is a common technique to grow advanced, functional crystalline materials. Supersaturation, temperature, and solvent composition are known to influence the growth rates and thereby properties of crystalline materials; however, a satisfactory explanation of how these factors affect the activation barrier for growth rates has not been developed. We report here that these effects can be attributed to a previously unrecognized consequence of solvent fluctuations in the solvation shell of solute molecules attaching to the crystal surface. With increasing supersaturation, the average hydration number of the glutamic acid molecule decreases and can reach an asymptotic limit corresponding to the number of adsorption sites on the molecule. The hydration number of the glutamic acid molecule also fluctuates due to the rapid exchange of solvent in the solvation shell and local variation in the supersaturation. These rapid fluctuations allow quasi-equilibrium between fully solvated and partially desolvated states of molecules, which can be used to construct a double-well potential and thereby to identify the transition state and the required activation barrier. The partially desolvated molecules are not stable and can attach spontaneously to the crystal surface. The activation barrier versus hydration number follows the Evans-Polanyi relation. The predicted absolute growth rates of the α-glutamic acid crystal at lower supersaturations are in reasonable agreement with the experimental observations.