Transferable Potentials for Phase Equilibria. 9. Explicit Hydrogen Description of Benzene and Five-Membered and Six-Membered Heterocyclic Aromatic Compounds

Quantifying the Uncertainty of Force Field Selection on Adsorption Predictions in MOFs

Comparisons between simulated and experimental adsorption isotherms in MOFs are fraught with challenges. On the experimental side, there is significant variation between isotherms measured on the same system, with a significant percentage (∼20%) of published data being considered outliers. On the simulation side, force fields are often chosen "off-the-shelf" with little or no validation. The effect of this choice on the reliability of simulated adsorption predictions has not yet been rigorously quantified. In this work, we fill this gap by systematically quantifying the uncertainty arising from force field selection on adsorption isotherm predictions. We choose methane adsorption, where electrostatic interactions are negligible, to independently study the effect of the framework Lennard-Jones parameters on a series of prototypical materials that represent the most widely studied MOF "families". Using this information, we compute an adsorption "consensus isotherm" from simulations, including a quantification of uncertainty, and compare it against a manually curated set of experimental data from the literature. By considering many experimental isotherms measured by different groups and eliminating outliers in the data using statistical analysis, we conduct a rigorous comparison that avoids the pitfalls of the standard approach of comparing simulation predictions to a single experimental data set. Our results show that (1) the uncertainty in simulated isotherms can be as large as 15% and (2) standard force fields can provide reliable predictions for some systems but can fail dramatically for others, highlighting systematic shortcomings in those models. Based on this, we offer recommendations for future simulation studies of adsorption, including high-throughput computational screening of MOFs.

Elucidating the structure of donor-acceptor conjugated polymer aggregates in liquid solution

High-spin donor-acceptor conjugated polymers are extensively studied for their potential applications in magnetic and spintronic devices. Inter-chain charge transfer among these high-spin polymers mainly depends on the nature of the local structure of the thin film and π-stacking between the polymer chains. However, the microscopic structural details of high-spin polymeric materials are rarely studied with an atomistic force field, and the molecular-level local structure in the liquid phase remains ambiguous. Here, we have examined the effects of oligomer chain length, side chain, and processing temperature on the organization of the high-spin cyclopentadithiophene-benzobisthiadiazole donor-acceptor conjugated polymer in chloroform solvent. We find that the oligomers display ordered aggregates whose structure depends on their chain length, with an average π-stacking distance of 3.38 ± 0.03 Å (at T = 298 K) in good agreement with the experiment. Also, the oligomers with longer alkyl side chains show better solvation and a shorter π-stacking distance. Furthermore, the clusters grow faster at higher temperature with more ordered aggregation between the oligomer chains.

Effects of Fluorinated Functionalization of Linker on Quercetin Encapsulation, Release and Hela Cell Cytotoxicity of Cu-Based MOFs as Smart pH-Stimuli Nanocarriers

Controlled delivery of target molecules is required in many medical and chemical applications. For such purposes, metal-organic frameworks (MOFs), which possess desirable features such as high porosity, large surface area, and adjustable functionalities, hold great potential as drug carriers. Herein, Quercetin (QU), as an anticancer drug, was loaded on Cu2 (BDC)2 (DABCO) and Cu2 (F4 BDC)2 )DABCO) MOFs (BDC=1,4-benzenedicarboxylate and DABCO=1,4-diazabicyclo[2.2.2]octane). As these Cu-MOFs have a high surface area, an appropriate pore size, and biocompatible ingredients, they can be utilized to deliver QU. The loading efficiency of QU in these MOFs was 49.5 % and 41.3 %, respectively. The drug-loaded compounds displayed sustained drug release over 15 days, remarkably high drug loading capacities and pH-controlled release behavior. The prepared nanostructures were characterized by different characterization technics including FT-IR, PXRD, ZP, TEM, FE-SEM, UV-vis, and BET. In addition, MTT assays were carried out on the HEK-293 and HeLa cell lines to investigate cytotoxicity. Cellular apoptosis analysis was performed to investigate the cell death mechanisms. Grand Canonical Monte Carlo simulations were conducted to analyze the interactions between MOFs and QU. Moreover, the stability of MOFs was also investigated during and after the drug release process. Ultimately, kinetic models of drug release were evaluated.

Effect of Ligand Aromaticity on Cyclohexane and Benzene Sorption in IRMOFs: A Computational Study

A series of four isoreticular MOFs (IRMOF-1, -10, -14, and -16) were selected for a computational investigation of the effect of ligand aromaticity on the adsorption capacity of an aromatic VOC (benzene) compared to its nonaromatic analog (cyclohexane). The affinity of the adsorbates was evaluated by calculating Henry's constants and adsorption enthalpies. It has been evidenced that while KH values decrease with ligand elongation (IRMOF-10 and -16), inserting a pyrene core into the MOF structure (IRMOF-14) increases both the cyclohexane and benzene adsorption efficiency by ∼290 and 54%, respectively. To elucidate host-guest interactions, we sought to locate preferential adsorption sites in MOF structures for the two VOCs studied by using the GCMC method. It appears that benzene interacts with the metal center (Zn4O clusters) and most of the ligand while cyclohexane tends to localize preferentially only near the Zn4O clusters. Coadsorption isotherms (equimolar mixture of benzene and cyclohexane) demonstrated the preferential adsorption of cyclohexane due to the stronger affinity for the MOF structure. On the other hand, for other isoreticular structures, the ligand elongation leads to a shift of the adsorption curve of cyclohexane caused by pore size increase and therefore less interactions with the walls. This phenomenon is counterbalanced in the case of IRMOF-14 due to stronger interactions between the cyclohexane and pyrene groups. The present results thus open perspectives in the design of promising MOF candidates for high-performing separation and sorption/detection of hydrocarbon VOCs.

Extension of the MolMod Database to Transferable Force Fields

MolMod, a web-based database for classical force fields for molecular simulations of fluids [Mol. Sim. 45, 10 (2019), 806-814], was extended to transferable force fields. Eight transferable force fields, including all-atom and united-atom type force fields, were implemented in the MolMod database: OPLS-UA, OPLS-AA, COMPASS, CHARMM, GROMOS, TraPPE, Potoff, and TAMie. These transferable force fields cover a large variety of chemical substance classes. The system is designed such that new transferable force fields can be readily integrated. A graphical user interface was implemented that enables the construction of molecules. The MolMod database compiles the force field for the specified component and force field type and provides the corresponding data and meta data as well as ready-to-use input files for the molecule for different simulation engines. This helps the user to flexibly choose molecular models and integrate them swiftly in their individual workflows, reducing risks of input errors in molecular simulations.

Predicting Critical Properties and Acentric Factors of Fluids Using Multitask Machine Learning

Knowledge of critical properties, such as critical temperature, pressure, density, as well as acentric factor, is essential to calculate thermo-physical properties of chemical compounds. Experiments to determine critical properties and acentric factors are expensive and time intensive; therefore, we developed a machine learning (ML) model that can predict these molecular properties given the SMILES representation of a chemical species. We explored directed message passing neural network (D-MPNN) and graph attention network as ML architecture choices. Additionally, we investigated featurization with additional atomic and molecular features, multitask training, and pretraining using estimated data to optimize model performance. Our final model utilizes a D-MPNN layer to learn the molecular representation and is supplemented by Abraham parameters. A multitask training scheme was used to train a single model to predict all the critical properties and acentric factors along with boiling point, melting point, enthalpy of vaporization, and enthalpy of fusion. The model was evaluated on both random and scaffold splits where it shows state-of-the-art accuracies. The extensive data set of critical properties and acentric factors contains 1144 chemical compounds and is made available in the public domain together with the source code that can be used for further exploration.

A Palladium-Based MOF for the Preferential Sorption of Benzene

Selective sorption of volatile aromatic compounds is a challenging issue for their total abatement. Despite the well-known affinity of palladium toward rich π systems, studies dedicated to volatile organic compound (VOC) capture with Pd(II)-based metal-organic frameworks (MOFs) are still very scarce. Intending to shed more light on this complex topic, this work compares the adsorption properties of two isostructural MOFs [Cu(2-pymo)₂]_n and [Pd(2-pymo)₂]_n and their selectivity for the sorption of linear, cyclic, or aromatic VOCs. The combination of both experimental and computational investigations highlights an increasing aromatic affinity over saturated hydrocarbons when palladium is chosen as a metal center (n_Benzene/n_n-hexane = 1.8 at 0.5 p/p₀) in the MOF instead of copper (n_Benzene/n_n-hexane = 0.7 at 0.5 p/p₀). Furthermore, [Pd(2-pymo)₂]_n clearly exhibits preferential adsorption of benzene over toluene (n_Benzene/n_Toluene = 1.7 at 0.5 p/p₀), due to the steric hindrance effects of the latter. The present results clearly underline the attractiveness of Pd-based MOFs for the design of selective aromatic adsorbents. Moreover, they also highlight the [Pd(2-pymo)₂]_n MOF as a relevant candidate for the selective capture of benzene, by a synergistic combination of both charge interactions and steric hindrance effects.

Modeling Alkyl Aromatic Hydrocarbons with Dissipative Particle Dynamics

Building on previous work studying alkanes, we develop a dissipative particle dynamics (DPD) model to capture the behavior of the alkyl aromatic hydrocarbon family under ambient conditions of 298 K and 1 atmosphere. Such materials are of significant worldwide industrial importance in applications such as solvents, chemical intermediates, surfactants, lubricating oils, hydraulic fluids, and greases. We model both liquids and waxy solids for molecules up to 36 carbons in size and demonstrate that we can correctly capture both the freezing transition and liquid-phase densities in pure substances and mixtures. We also demonstrate the importance of including specialized bead types into the DPD model (rather than solely relying on generic bead types) to capture specific local geometrical constructs such as the benzene ring found in the benzyl chemical group; this can be thought of as representing subtle real-world many-body effects via customized pairwise non-bonded potentials.

Predicting the Solubility of Nonelectrolyte Solids Using a Combination of Molecular Simulation with the Solubility Parameter Method MOSCED: Application to the Wastewater Contaminants Monuron, Diuron, Atrazine and Atenolol

Methods to predict the equilibrium solubility of nonelectrolyte solids are indispensable for early-stage process development, design, and feasibility studies. Conventional analytic methods typically require reference data to regress parameters, which may not be available or limited for novel systems. Molecular simulation is a promising alternative, but is computationally intensive. Here, we demonstrate the ability to use a small number of molecular simulation free energy calculations to generate reference data to regress model parameters for the analytical MOSCED (modified separation of cohesive energy density) model. The result is an efficient analytical method to predict the equilibrium solubility of nonelectrolyte solids. The method is demonstrated for the wastewater contaminants monuron, diuron, atrazine and atenolol. Predictions for monuron, diuron and atrazine are in reasonable agreement with MOSCED parameters regressed using experimental solubility data. Predictions for atenolol are inferior, suggesting a potential limitation in the adopted molecular models, or the solvents selected to generate the necessary reference data.

Excess chemical potential of thiophene in [C4MIM] [BF4, Cl, Br, CH3COO] ionic liquids, determined by molecular simulations

The excess chemical potential of thiophene in imidazolium-based ionic liquids [C₄mim][BF₄], [C₄mim][Cl], [C₄mim][Br], and [C₄mim][CH₃COO] were determined by means of molecular dynamics in conjunction with free energy perturbation techniques employing non-polarizable force fields at 300 K and 343.15 K. In addition, energetic and structural analysis were performed such as: interaction energies, averaged noncovalent interactions, radial, and combined distribution functions. The results from this work revealed that the ionic liquids (ILs) presenting the most favorable excess chemical potentials ([C₄mim][BF₄], [C₄mim][CH₃COO]) are associated with the strongest energetic interaction between the thiophene molecule and the ionic liquid anion, and with the weakest energetic interaction between the thiophene molecule and the ionic liquid cation.

Excess chemical potential of thiophene in imidazolium-based ionic liquids [C₄mim][BF₄], [C₄mim][Cl], [C₄mim][Br], and [C₄mim][CH₃COO] determined by molecular simulations.

Simple corrections for the static dielectric constant of liquid mixtures from model force fields

Pair-wise additive force fields provide fairly accurate predictions, through classical molecular simulations, for a wide range of structural, thermodynamic, and dynamical properties of many materials. However, one key property that has not been well captured is the static dielectric constant, which characterizes the response of a system to an applied electric field and is important in determining the screening of electrostatic interactions through a system. A simple correction has been found to provide a relatively robust method to improve the estimate of the static dielectric constant from molecular simulations for a broad range of compounds. This approach accounts for the electronic contribution to molecular polarizability and assumes that the charges that couple a molecule to an applied electric field are proportional to the effective force field charges. In this work, we examine how this correction performs for systems at different temperatures and for binary mixtures. Using a value for the electronic polarizability, based on the experimental index of refraction, and a charge scaling factor, determined at a single temperature, we find that the static dielectric constant can be predicted remarkably well, in comparison to the experimentally measured values. This provides good evidence that the effective charges that appear in pair-wise additive force fields developed to reproduce the potential energy surface of a system are not the same as those that determine the static dielectric constant; however, they can be captured in a relatively simple manner, which is dependent on the particular force field.

Thermodynamic, structural and dynamic properties of selected non-associative neat liquids

Non-associative neat liquids benzene, acetone and carbon tetrachloride have been examined via molecular dynamics simulations. Several models of each neat liquid have been simulated, and selected thermodynamic and structural results are presented. However, the models have been compared mainly in terms of their dynamic quantities. Since models are rarely parametrized with the dynamic properties in mind, the principal goal of this work is to present quantities such as the power spectra, rotational correlation functions and relaxation times, diffusion coefficients and self and distinct parts of the van Hove functions in relation to available experimental data. The general trends of the calculated data provide a benchmark for the behavior of neat simple liquids which will be built upon in the cases of mixtures with associative liquids.

Molecular Mechanism for Azeotrope Formation in Ethanol/Benzene Binary Mixtures through Gibbs Ensemble Monte Carlo Simulation

Azeotropes have been studied for decades due to the challenges they impose on separation processes but fundamental understanding at the molecular level remains limited. Although molecular simulation has demonstrated its capability of predicting mixture vapor-liquid equilibrium (VLE) behaviors, including azeotropes, its potential for mechanistic investigation has not been fully exploited. In this study, we use the united atom transferable potentials for phase equilibria (TraPPE-UA) force field to model the ethanol/benzene mixture, which displays a positive azeotrope. Gibbs ensemble Monte Carlo (GEMC) simulation is performed to predict the VLE phase diagram, including an azeotrope point. The results accurately agree with experimental measurements. We argue that the molecular mechanism of azeotrope formation cannot be fully understood by studying the mixture liquid-state stability at the azeotrope point alone. Rather, azeotrope occurrence is only a reflection of the changing relative volatility between the two components over a much wider composition range. A thermodynamic criterion is thus proposed on the basis of the comparison of partial excess Gibbs energy between the components. In the ethanol/benzene system, molecular energetics shows that with increasing ethanol mole fraction, its volatility initially decreases but later plateaus, while benzene volatility is initially nearly constant and only starts to decrease when its mole fraction is low. Analysis of the mixture liquid structure, including a detailed investigation of ethanol hydrogen-bonding configurations at different composition levels, reveals the underlying molecular mechanism for the changing volatilities responsible for the azeotrope.

A metal-organic framework for efficient water-based ultra-low-temperature-driven cooling

Efficient use of energy for cooling applications is a very important and challenging field in science. Ultra-low temperature actuated (Tdriving < 80 °C) adsorption-driven chillers (ADCs) with water as the cooling agent are one environmentally benign option. The nanoscale metal-organic framework [Al(OH)(C6H2O4S)] denoted CAU-23 was discovered that possess favorable properties, including water adsorption capacity of 0.37 gH2O/gsorbent around p/p0 = 0.3 and cycling stability of at least 5000 cycles. Most importantly the material has a driving temperature down to 60 °C, which allows for the exploitation of yet mostly unused temperature sources and a more efficient use of energy. These exceptional properties are due to its unique crystal structure, which was unequivocally elucidated by single crystal electron diffraction. Monte Carlo simulations were performed to reveal the water adsorption mechanism at the atomic level. With its green synthesis, CAU-23 is an ideal material to realize ultra-low temperature driven ADC devices.

Curvature Dependence of the Mass Accommodation Coefficient

Mass accommodation coefficient, a parameter that captures molecular transport phenomena at liquid-vapor interfaces, is essential for predicting the growth of liquid droplets during condensation processes but is difficult to obtain experimentally. Molecular simulations have been widely used to obtain accommodation coefficients for planar interfaces, but the applicability of planar accommodation coefficients to the high-curvature interfaces present in very small droplets is not clear. In this work, molecular dynamics simulations are used to compute equilibrium mass accommodation coefficients at different temperatures for small droplets of various fluids, including Lennard-Jones and Buckingham fluids, benzene, butane, methane, methanol, and water. For all fluids studied, the mass accommodation coefficient increases with droplet size to a constant limiting value and decreases with temperature. Furthermore, the accommodation coefficient curvature dependence collapses onto a universal curve when appropriately scaled.

Predicting vapor liquid equilibria using density functional theory: A case study of argon

Predicting vapor liquid equilibria (VLE) of molecules governed by weak van der Waals (vdW) interactions using the first principles approach is a significant challenge. Due to the poor scaling of the post Hartree-Fock wave function theory with system size/basis functions, the Kohn-Sham density functional theory (DFT) is preferred for systems with a large number of molecules. However, traditional DFT cannot adequately account for medium to long range correlations which are necessary for modeling vdW interactions. Recent developments in DFT such as dispersion corrected models and nonlocal van der Waals functionals have attempted to address this weakness with a varying degree of success. In this work, we predict the VLE of argon and assess the performance of several density functionals and the second order Møller-Plesset perturbation theory (MP2) by determining critical and structural properties via first principles Monte Carlo simulations. PBE-D3, BLYP-D3, and rVV10 functionals were used to compute vapor liquid coexistence curves, while PBE0-D3, M06-2X-D3, and MP2 were used for computing liquid density at a single state point. The performance of the PBE-D3 functional for VLE is superior to other functionals (BLYP-D3 and rVV10). At T = 85 K and P = 1 bar, MP2 performs well for the density and structural features of the first solvation shell in the liquid phase.

Transferable Potentials for Chloroethenes: Insights into Nonideal Solution Behavior of Environmental Contaminants

Predicting the nonideal phase behavior of binary and multicomponent systems remains a significant challenge for particle-based simulations. Here, we develop a transferable force field for chloroethenes, common environmental contaminants, that can accurately model the vapor liquid phase equilibria including azeotrope formation. The new all-atom force field reproduces saturated liquid densities, saturated vapor pressures, boiling points, and critical properties within 1, 10, 1, and 1% of the experiment data, respectively. Furthermore, the vapor liquid equilibria of trichloroethylene and 1-propanol binary mixture, which forms a minimum boiling point azeotrope, is predicted with a reasonable accuracy. The microstructure of neat and binary systems is explored using pair correlation functions and spatial distribution functions. As the new force field is consistent with transferable potentials for phase equilibria (TraPPE) force field, it expands the applicability of TraPPE force field to chloroethenes.

Ionic-Functionalized Polymers of Intrinsic Microporosity for Gas Separation Applications

Ionic-functionalized microporous materials are attractive for energy-efficient gas adsorption and separation processes and have shown promising results in gas mixtures at pressure ranges and compositions that are relevant for industrial applications. In this work, we studied the influence of different counterions (Li⁺, Na⁺, K⁺, Rb⁺, and Mg²⁺) on the porosity, carbon dioxide (CO₂) gas adsorption, and selectivity in ionic-functionalized PIM-1 (IonomIMs), a polymer belonging to the class of linear and amorphous microporous polymers known as polymers of intrinsic microporosity (PIMs). It was found that an increase in the concentration of ionic groups led to a decrease in the free volume, resulting in a less porous polymer framework, and Mg²⁺-functionalized IonomIMs exhibited a relatively larger porosity compared to other IonomIMs. The CO₂ adsorption capacity was affected by the different counterions for IonomIM-1, and a higher loading capacity for pure CO₂ was observed for Mg²⁺. Furthermore, the IonomIMs showed an enhanced CO₂ selectivity in CO₂/CH₄ and CO₂/N₂ gas mixtures at conditions used in pressure swing adsorption and vacuum swing adsorption applications. It was also observed that the concentration of ionic groups plays a vital role in changing the CO₂ gas adsorption and selectivity.

Performance of density functionals for modeling vapor liquid equilibria of CO2 and SO2

Vapor liquid equilibria (VLE) and condensed phase properties of carbon dioxide and sulfur dioxide are calculated using first principles Monte Carlo (FPMC) simulations to assess the performance of several density functionals, notably PBE-D3, BLYP-D3, PBE0-D3, M062X-D3, and rVV10. GGA functionals were used to compute complete vapor liquid coexistence curves (VLCCs) to estimate critical properties, while the hybrid and nonlocal van der Waals functionals were used only for computing density at a single state point due to the high computational cost. Our results show that the BLYP-D3 functional performs well in predicting VLE properties for both molecules when compared with other functionals. In the liquid phase, pair correlation functions reveal that there is not a significant difference in the location of the peak for the first solvation shell while the peak heights are different for different functionals. Overall, the BLYP-D3 functional is a good choice for modeling VLE of acidic gases with significant environmental implications such as CO₂ and SO₂ . © 2017 Wiley Periodicals, Inc.

Revealing at the molecular level the role of the surfactant in the enhancement of the thermal properties of the gold nanofluid system used for concentrating solar power

A surfactant acts as a kind of net around gold that enhances at the molecular level the gold nanofluid thermal properties.

Extension of the SAFT-VR Mie EoS To Model Homonuclear Rings and Its Parametrization Based on the Principle of Corresponding States

The statistical associating fluid theory of variable range employing a Mie potential (SAFT-VR-Mie) proposed by Lafitte et al. (J. Chem Phys. 2013, 139, 154504) is one of the latest versions of the SAFT family. This particular version has been shown to have a remarkable capability to connect experimental determinations, theoretical calculations, and molecular simulations results. However, the theoretical development restricts the model to chains of beads connected in a linear fashion. In this work, the capabilities of the SAFT-VR Mie equation of state for modeling phase equilibria are extended for the case of planar ring compounds. This modification proposed replaces the Helmholtz energy of chain formation by an empirical contribution based on a parallelism to the second-order thermodynamic perturbation theory for hard sphere trimers. The proposed expression is given in terms of an extra parameter, χ, that depends on the number of beads, m_s, and the geometry of the ring. The model is used to describe the phase equilibrium for planar ring compounds formed of Mie isotropic segments for the cases of m_s equals to 3, 4, 5 (two configurations), and 7 (two configurations). The resulting molecular model is further parametrized, invoking a corresponding states principle resulting in sets of parameters that can be used indistinctively in theoretical calculations or in molecular simulations without any further refinements. The extent and performance of the methodology has been exemplified by predicting the phase equilibria and vapor pressure curves for aromatic hydrocarbons (benzene, hexafluorobenzene, toluene), heterocyclic molecules (2,5-dimethylfuran, sulfolane, tetrahydro-2H-pyran, tetrahydrofuran), and polycyclic aromatic hydrocarbons (naphthalene, pyrene, anthracene, pentacene, and coronene). An important aspect of the theory is that the parameters of the model can be used directly in molecular dynamics (MD) simulations to calculate equilibrium phase properties and interfacial tensions with an accuracy that rivals other coarse grained and united atom models, for example, liquid densities, are predicted, with a maximum absolute average deviation of 3% from both the theory and the MD simulations, while the interfacial tension is predicted, with a maximum absolute average of 8%. The extension to mixtures is exemplified by considering a binary system of hexane (chain fluid) and tetrahydro-2H-pyran (ring fluid).

Molecular Insights into Gelation of Di-Fmoc-l-Lysine in Organic Solvent-Water Mixtures

Despite significant interest in molecular gels due to their intriguing structure formation through self-assembly and their stimuli-responsive behavior, our understanding of the gel formation mechanism of a low-molecular-weight gelator (LMWG) is incomplete. Here, we report a combined experimental and computational study on a LMWG, di-Fmoc-l-lysine, that has two aromatic moieties and multiple hydrogen bond donors and acceptors. Gelation in various organic solvent-water mixtures was obtained through the solvent-triggered technique. We show that an approach based on approximate cohesive energy density derived from density functional theory (DFT) calculations can capture the experimental solubility trend of LMWGs in different organic solvents. Furthermore, DFT calculations indicate parallel and helical structures to be the preferred structural motifs for gelator dimers. We believe that these motifs can potentially lead to fiber formation as observed with microscopy. Our work provides a relatively simple yet effective approach to quantify interactions between solvents and complex gelators that can help rationalize solubility and gelation behavior.

Adsorption and molecular siting of CO2, water, and other gases in the superhydrophobic, flexible pores of FMOF-1 from experiment and simulation

CO₂ isotherms for FMOF-1 reveal 11.0 mol L^–1 max uptake and suggest framework expansion, substantiated by in situ neutron diffraction and GCMC simulations.

FMOF-1 is a flexible, superhydrophobic metal–organic framework with a network of channels and side pockets decorated with –CF₃ groups. CO₂ adsorption isotherms measured between 278 and 313 K and up to 55 bar reveal a maximum uptake of ca. 6.16 mol kg^–1 (11.0 mol L^–1) and unusual isotherm shapes at the higher temperatures, suggesting framework expansion. We used neutron diffraction and molecular simulations to investigate the framework expansion behaviour and the accessibility of the small pockets to N₂, O₂, and CO₂. Neutron diffraction in situ experiments on the crystalline powder show that CO₂ molecules are favourably adsorbed at three distinct adsorption sites in the large channels of FMOF-1 and cannot access the small pockets in FMOF-1 at 290 K and oversaturated pressure at 61 bar. Stepped adsorption isotherms for N₂ and O₂ at 77 K can be explained by combining Monte Carlo simulations in several different crystal structures of FMOF-1 obtained from neutron and X-ray diffraction under different conditions. A similar analysis is successful for CO₂ adsorption at 278 and 283 K up to ca. 30 bar; however, at 298 K and pressures above 30 bar, the results suggest even more substantial expansion of the FMOF-1 framework. The measured contact angle for water on an FMOF-1 pellet is 158°, demonstrating superhydrophobicity. Simulations and adsorption measurements also show that FMOF-1 is hydrophobic and water is not adsorbed in FMOF-1 at room temperature. Simulated mixture isotherms of CO₂ in the presence of 80% relative humidity predict that water does not influence the CO₂ adsorption in FMOF-1, suggesting that hydrophobic MOFs could hold promise for CO₂ capture from humid gas streams.

Metal Ion Modeling Using Classical Mechanics

Metal ions play significant roles in numerous fields including chemistry, geochemistry, biochemistry, and materials science. With computational tools increasingly becoming important in chemical research, methods have emerged to effectively face the challenge of modeling metal ions in the gas, aqueous, and solid phases. Herein, we review both quantum and classical modeling strategies for metal ion-containing systems that have been developed over the past few decades. This Review focuses on classical metal ion modeling based on unpolarized models (including the nonbonded, bonded, cationic dummy atom, and combined models), polarizable models (e.g., the fluctuating charge, Drude oscillator, and the induced dipole models), the angular overlap model, and valence bond-based models. Quantum mechanical studies of metal ion-containing systems at the semiempirical, ab initio, and density functional levels of theory are reviewed as well with a particular focus on how these methods inform classical modeling efforts. Finally, conclusions and future prospects and directions are offered that will further enhance the classical modeling of metal ion-containing systems.