Adsorption-Induced Deformation of Zeolites 4A and 13X: Experimental and Molecular Simulation Study.

Gas adsorption in zeolites leads to adsorption-induced deformation, which can significantly affect the adsorption and diffusive properties of the system. In this study, we conducted both experimental investigations and molecular simulations to understand the deformation of zeolites 13X and 4A during carbon dioxide adsorption at 273 K. To measure the sample's adsorption isotherm and strain simultaneously, we used a commercial sorption instrument with a custom-made sample holder equipped with a dilatometer. Our experimental data showed that while the zeolites 13X and 4A exhibited similar adsorption isotherms, their strain isotherms differed significantly. To gain more insight into the adsorption process and adsorption-induced deformation of these zeolites, we employed coupled Monte Carlo and molecular dynamics simulations with atomistically detailed models of the frameworks. Our modeling results were consistent with the experimental data and helped us identify the reasons behind the different deformation behaviors of the considered structures. Our study also revealed the sensitivity of the strain isotherm of zeolites to pore size and other structural and energetic features, suggesting that measuring adsorption-induced deformation could serve as a complementary method for material characterization and provide guidelines for related technical applications.

Closer Look at Adsorption of Sarin and Simulants on Metal-Organic Frameworks.

The development of effective protection against exposure to chemical warfare agents (CWAs), such as sarin, relies on studies of its adsorption on the capturing materials and seeking candidates capable of adsorbing large amounts of sarin gas. Many metal-organic frameworks (MOFs) are promising materials for the effective capture and degradation of sarin and simulant substances. Among the simulants capable of mimicking thermodynamic properties of the agent, not all of them have been investigated on the ability to act similarly in the adsorption process, in particular, whether the agent and a simulant have similar mechanisms of binding to the MOF surface. Molecular simulation studies not only provide a safe way to investigate the aforementioned processes but can also help reveal the mechanisms of interactions between the adsorbents and the adsorbing compounds at the molecular level. We performed Monte Carlo simulations of the adsorption of sarin and three simulants, dimethyl methylphosphonate (DMMP), diisopropyl methylphosphonate (DIMP), and diisopropyl fluorophosphate (DIFP), on selected MOFs that have previously shown strong capabilities to adsorb sarin. On the basis of the calculated adsorption isotherms, enthalpy of adsorption, and radial distribution functions, we revealed common mechanisms among the particularly efficient adsorbents as well as the ability of simulants to mimic them. The findings can help in selecting a suitable simulant compound to study CWA adsorption on MOFs and guide further synthesis of efficient MOFs for the capture of organophosphorus compounds.

Molecular Simulations of Vapor-Liquid Equilibrium of Isocyanates.

The wide range of applications of the isocyanates across multiple industries sparks the interest in the study of their phase behavior. A molecular simulation is a powerful tool that can go beyond experimental investigations relying on a molecular structure of a chemical. The success of a molecular simulation relies on a description of the system, namely, force field, and its parameterization on reproducing properties of interest. In this work, we propose a united-atom force field based on the transferable potentials for phase equilibria (TraPPE) to model the vapor-liquid phase behavior of isocyanates. With Monte Carlo and molecular dynamics simulation methods and the introduced force field, we modeled vapor-liquid equilibrium for a family of linear mono-isocyanates, from methyl isocyanate to hexyl isocyanate, and hexamethylene diisocyanate. Additionally, we performed similar calculations for methyl, ethyl, and butyl isocyanates based on the all-atom GAFF-IC force field available in the literature for modeling isocyanate viscosities. We showed that the developed TraPPE-based force field generally overperformed the GAFF-IC force field and overall showed excellent performance in modeling phase behavior of isocyanates. Based on the simulated vapor pressures for the considered compounds, we estimated the Antoine equation parameters to calculate the vapor pressure in a range of temperatures. The predictions are of particular use in the investigation of thermodynamic properties for those isocyanates lacking experimental vapor pressure data. Results can also be employed in modeling the phase behavior of isocyanate mixtures to investigate their sensing and capturing. Furthermore, from the vapor-liquid equilibrium binodals, we predicted the critical properties of isocyanates which can be used in thermodynamic models based on an equation of state.

Force Fields for Molecular Modeling of Sarin and its Simulants: DMMP and DIMP.

Even three decades after signing the Chemical Weapons Convention, organophosphorus chemical warfare agents (CWAs), such as sarin, remain a threat. The development of novel methods for the detection of CWAs, protection from CWAs, and CWA decontamination motivates research on their physicochemical properties. Due to the extreme toxicity of sarin, most of the experimental studies are carried out using less toxic simulant compounds. In addition to experimental studies of sarin simulants, both sarin and simulants can be studied using in silico experiments-molecular simulations. The results of classical molecular modeling of the compounds and their agreement with experimental data rely on the force field used to describe the system. In recent years, there have been several force fields proposed for sarin and its most common simulant dimethyl methylphosphonate (DMMP). However, other simulants frequently used in experiments received less attention from the molecular simulation perspective, for example, to date, there is no force field and no simulation data for diisopropyl methylphosphonate (DIMP). Here, we compare the literature force fields for sarin and DMMP, focusing specifically on the vapor-liquid equilibrium for the pure species. We carried out Monte Carlo and molecular dynamics simulations using the existing literature force fields from which we predicted the liquid densities and vapor pressures developing the entire binodal curves. We compared the predictions to the experimental data and showed that the TraPPE-UA force field outperformed the other force fields. Thus, we extended TraPPE-UA for DIMP, utilized this force field in molecular simulations, and predicted the liquid densities and vapor pressures for a range of temperatures (binodal curve), which agreed well with the published experimental data. From the binodal, we calculated the critical properties of DIMP and demonstrated that these parameters can be used in the Peng-Robinson equation of state for this compound.

Porous Structure of Silica Colloidal Crystals.

We prepared silica colloidal crystals with different pore sizes using isothermal heating evaporation-induced self-assembly in quantities suitable for nitrogen porosimetry and studied their porous structure. We observed pores of two types in agreement with the description of silica colloidal crystals as face-centered cubic packed structures containing octahedral and tetrahedral voids. We calculated the sizes of these pores using the Derjaguin-Broekhoff-de Boer theory of capillary condensation for spherical pores. We also described the pore geometry mathematically and showed that the octahedral pore radii measured experimentally matches closely the radii of the spheres of the same volume. In the case of the tetrahedral pores, the proposed approach underestimated the pore radius by ca. 40%. Overall, this simple geometrical description provides a good representation of the porous system in silica colloidal crystals.