In this thesis, molecular simulations of porous materials are performed, ranging from testing theories using simple models to predicting thermodynamic properties of realistic complex systems. Porous materials such as zeolites, MOFs, and cyclodextrins (CDs) are studied for waste water treatment and energy related applications. The purpose of the presented studies is twofold: (1) identify the best performing materials for a specific application from a large pool of structures, (2) obtain molecular insight into the molecular level phenomena involved in a specific application. These studies aim to provide guidance for experimentalist in finding new, better performing materials, and a more detailed understanding of the involved phenomena...

Cyclodextrins (CDs) are widely used in drug delivery, catalysis, food and separation processes. In this work, a comprehensive simulation study on the diffusion of the native α-, β- and γ-CDs in aqueous solutions is carried out using Molecular Dynamics simulations. The effect of the system size on the computed self-diffusivity is investigated and it is found that the required correction can be as much as 75% of the final value. The effect of the water force field is examined and it is shown that the q4md-CD/TIP4P/2005 force field combination predicts the experimentally measured self-diffusion coefficients of CDs very accurately. The self-diffusion coefficients of the three native CDs were also computed in aqueous-NaCl solutions using the Joung and Cheatham (JC) and the Madrid-2019 force fields. It is found that Na⁺ ions have higher affinity towards the CDs when the JC force field is used and for this reason the predicted diffusivity of CDs is lower compared to simulations using the Madrid-2019 force field. As a model system for drug delivery and waste-water treatment applications, the diffusion of the β-CD:Ibuprofen inclusion complex in water is studied. In agreement with experiments for similar components, it is shown that the inclusion complex and the free β-CD have almost equal self-diffusion coefficients. Our analysis revealed that this is most likely caused by the almost full inclusion of the ibuprofen in the cavity of the β-CD. Our findings show that Molecular Dynamics simulation can be used to provide reasonable diffusivity predictions, and to obtain molecular-level understanding useful for industrial applications of CDs.

According to the ISO 14687-2:2019 standard, the water content of H2 fuel for transportation and stationary applications should not exceed 5 ppm (molar). To achieve this water content, zeolites can be used as a selective adsorbent for water. In this work, a computational screening study is carried out for the first time to identify potential zeolite frameworks for the drying of high-pressure H2 gas using Monte Carlo (MC) simulations. We show that the Si/Al ratio and adsorption selectivity have a negative correlation. 218 zeolites available in the database of the International Zeolite Association are considered in the screening. We computed the adsorption selectivity of each zeolite for water from the high-pressure H2 gas having water content relevant to vehicular applications and near saturation. It is shown that due to the formation of water clusters, the water content in the H2 gas has a significant effect on the selectivity of zeolites with a helium void fraction larger than 0.1. Under each operating condition, five most promising zeolites are identified based on the adsorption selectivity, the pore limiting diameter, and the volume of H2 gas that can be dried by 1 dm3 of zeolite. It is shown that at 12.3 ppm (molar) water content, structures with helium void fractions smaller than 0.07 are preferred. The structures identified for 478 ppm (molar) water content have helium void fractions larger than 0.26. The proposed zeolites can be used to dry 400-8000 times their own volume of H2 gas depending on the operating conditions. Our findings strongly indicate that zeolites are potential candidates for the drying of high-pressure H2 gas.

The accurate description of the behavior of fluids in nanoporous materials is of great importance for numerous industrial applications. Recently, a new approach was reported to calculate the pressure of nanoconfined fluids. In this approach, two different pressures are defined to take into account the smallness of the system: the so-called differential and the integral pressures. Here, the effect of several factors contributing to the confinement of fluids in nanopores are investigated using the definitions of the differential and integral pressures. Monte Carlo (MC) simulations are performed in a variation of the Gibbs ensemble to study the effect of the pore geometry, fluid-wall interactions, and differential pressure of the bulk fluid phase. It is shown that the differential and integral pressure are different for small pores and become equal as the pore size increases. The ratio of the driving forces for mass transport in the bulk and in the confined fluid is also studied. It is found that, for small pore sizes (i.e., < 5σfluid ), the ratio of the two driving forces considerably deviates from 1.

It is known that thermodynamic properties of a system change upon confinement. To know how, is important for modelling of porous media. We propose to use Hill’s systematic thermodynamic analysis of confined systems to describe two-phase equilibrium in a nanopore. The integral pressure, as defined by the compression energy of a small volume, is then central. We show that the integral pressure is constant along a slit pore with a liquid and vapor in equilibrium, when Young and Young–Laplace’s laws apply. The integral pressure of a bulk fluid in a slit pore at mechanical equilibrium can be understood as the average tangential pressure inside the pore. The pressure at mechanical equilibrium, now named differential pressure, is the average of the trace of the mechanical pressure tensor divided by three as before. Using molecular dynamics simulations, we computed the integral and differential pressures, ̂p and p, respectively, analysing the data with a growing-core methodology. The value of the bulk pressure was confirmed by Gibbs ensemble Monte Carlo simulations. The pressure difference times the volume, V, is the subdivision potential of Hill, (p − ̂p)V = ɛ. The combined simulation results confirm that the integral pressure is constant along the pore, and that ɛ/V scales with the inverse pore width. This scaling law will be useful for prediction of thermodynamic properties of confined systems in more complicated geometries.

Recently, β-cyclodextrin (βCD)-based polymers with enhanced adsorption kinetics and high removal capacity of organic micropollutants (OMPs) and uptake rates have been synthesized and tested experimentally. Although the exact physical-chemical mechanisms via which these polymers capture the various types of OMPs are not yet fully understood, it is suggested that the inclusion complex formation of OMPs with βCD is very important. In this study, the inclusion complex formation of OMPs with βCD in an aqueous solution is investigated by using the well-established attach-pull-release method in force field-based molecular dynamics simulations. A representative set of OMPs is selected based on the measured occurrences in surface and ground waters and the directives published by the European Union. To characterize the formation of the inclusion complex, the binding free energies, enthalpies, and entropies are computed and compared to experimental values. It is shown that computations using the q4md-CD/GAFF/Bind3P force field combination yield binding free energies that are in reasonable agreement with the experimental results for all OMPs studied. The binding enthalpies are decomposed into the main contributing interaction types. It is shown that, for all studied OMPs, the van der Waals interactions are favorable for the inclusion complexion and the hydrogen bond formation of the guest with the solvent and βCD plays a crucial role in the binding mechanism. Our findings show that MD simulations can adequately describe the inclusion complex formation of βCD with OMPs, which is the first step toward understanding the underlying mechanisms via which the βCD-based polymers capture OMPs.

The adsorption of organic micropollutants (OMPs) on high-silica zeolites is characterized by adsorption isotherms with various shapes. The occurrence of an S-shaped adsorption isotherm indicates the lack of adsorption affinity for OMPs at low, environmentally relevant equilibrium concentrations. In this study, S-shaped isotherms were observed during batch experiments with 2,4,6-trichlorophenol (TCP) and FAU zeolites. This is the first time that an S-shaped isotherm is reported for the adsorption of OMPs on high-silica zeolites. Monte Carlo (MC) simulations in the grand-canonical ensemble were used to obtain a better understanding of the mechanism of the S-shaped adsorption isotherms. From the MC simulation results, it was observed that multiple TCP molecules were adsorbed in the supercages of the FAU zeolites. It was found that the π-π interactions between TCP molecules give rise to the adsorption of multiple TCP molecules per supercage, and thus causing an S-shaped adsorption isotherm. Simulations also revealed that water molecules were preferentially adsorbed in the supercages and sodalite cages of the FAU zeolites. FAU zeolites with a higher Al content adsorbed a higher amount of water molecules and a lower amount of TCP, and showed less pronounced S-shaped isotherms.

We present a new molecular simulation code, Brick-CFCMC, for performing Monte Carlo simulations using state-of-the-art simulation techniques. The Continuous Fractional Component (CFC) method is implemented for simulations in the NVT/NPT ensembles, the Gibbs Ensemble, the Grand-Canonical Ensemble, and the Reaction Ensemble. Molecule transfers are facilitated by the use of fractional molecules which significantly improve the efficiency of the simulations. With the CFC method, one can obtain phase equilibria and properties such as chemical potentials and partial molar enthalpies/volumes directly from a single simulation. It is possible to combine trial moves from different ensembles. This enables simulations of phase equilibria in a system where also a chemical reaction takes place. We demonstrate the applicability of our software by investigating the esterification of methanol with acetic acid in a two-phase system.

We present a detailed overview of classical molecular simulation studies examining the self-diffusion coefficient of water. The self-diffusion coefficient is directly associated with the calculations of tracer or mutual diffusion coefficient of mixtures and, therefore, is a fundamental transport property, essential for an accurate description of mass transfer processes in biological, geological (i.e. energy or environmentally related), and chemical systems. In the current review we explore two distinct research areas. Namely, we discuss the self-diffusion of water in the bulk phase and under confinement. Different aspects that affect the diffusion process, including the molecular models, the system-size effects, the temperature and pressure conditions and the type of confinement are discussed. Finally, possible directions for future research are outlined.

A computational screening of 2930 experimentally synthesized metal-organic frameworks (MOFs) is carried out to find the best-performing structures for adsorption-driven cooling (AC) applications with methanol and ethanol as working fluids. The screening methodology consists of four subsequent screening steps for each adsorbate. At the end of each step, the most promising MOFs for AC application are selected for further investigation. In the first step, the structures are selected on the basis of physical properties (pore limiting diameter). In each following step, points of the adsorption isotherms of the selected structures are calculated from Monte Carlo simulations in the grand-canonical ensemble. The most promising MOFs are selected on the basis of the working capacity of the structures and the location of the adsorption step (if present), which can be related to the applicable operational conditions in AC. Because of the possibility of reversible pore condensation (first-order phase transition), the mid-density scheme is used to efficiently and accurately determine the location of the adsorption step. At the end of the screening procedure, six MOFs with high deliverable working capacities (∼0.6 mL working fluid in 1 mL structure) and diverse adsorption step locations are selected for both adsorbates from the original 2930 structures. Because the highest experimentally measured deliverable working capacity to date for MOFs with methanol is ca. 0.45 mL mL^-1, the selected six structures show the potential to improve the efficiency of ACs.