Calculating free energies for diffusion in tight-fitting zeolite-guest systems: Local normal-mode Monte Carlo

Simulating adsorptive expansion of zeolites: application to biomass-derived solutions in contact with silicalite

We have constructed and applied an algorithm to simulate the behavior of zeolite frameworks during liquid adsorption. We applied this approach to compute the adsorption isotherms of furfural-water and hydroxymethyl furfural (HMF)-water mixtures adsorbing in silicalite zeolite at 300 K for comparison with experimental data. We modeled these adsorption processes under two different statistical mechanical ensembles: the grand canonical (V-Nz-μg-T or GC) ensemble keeping volume fixed, and the P-Nz-μg-T (osmotic) ensemble allowing volume to fluctuate. To optimize accuracy and efficiency, we compared pure Monte Carlo (MC) sampling to hybrid MC-molecular dynamics (MD) simulations. For the external furfural-water and HMF-water phases, we assumed the ideal solution approximation and employed a combination of tabulated data and extended ensemble simulations for computing solvation free energies. We found that MC sampling in the V-Nz-μg-T ensemble (i.e., standard GCMC) does a poor job of reproducing both the Henry's law regime and the saturation loadings of these systems. Hybrid MC-MD sampling of the V-Nz-μg-T ensemble, which includes framework vibrations at fixed total volume, provides better results in the Henry's law region, but this approach still does not reproduce experimental saturation loadings. Pure MC sampling of the osmotic ensemble was found to approach experimental saturation loadings more closely, whereas hybrid MC-MD sampling of the osmotic ensemble quantitatively reproduces such loadings because the MC-MD approach naturally allows for locally anisotropic volume changes wherein some pores expand whereas others contract.

Simulating the Relaxation Dynamics of Microwave-Driven Zeolites

We have performed equilibrium and nonequilibrium molecular dynamics simulations to study how microwave (MW)-heated zeolite systems relax to thermal equilibrium. We have simulated the relaxation of both ionic and dipolar phases in FAU-type zeolites, finding biexponential relaxation in all cases studied. Fast-decay times were uniformly below 1 ps, while slow-decay times were found to be as long as 14 ps. Fast-decay times increase with an increase in the initial temperature difference between MW-heated ions/dipoles and the equilibrium system. Slow-decay times were found to be relatively insensitive to the details of the MW-heated nonequilibrium state. Velocity, force, and orientational correlation functions, calculated at equilibrium to explore the natural dynamics of energy transfer, decay well before 1 ps and show little evidence of biexponential decay. In contrast, kinetic energy correlation functions show strong biexponential behavior with slow-decay times as long as 14 ps. We suggest a two-step mechanism involving initial, efficient energy transfer mediated by strongly anharmonic zeolite-guest forces, followed by a slower process mediated by weakly anharmonic couplings among normal modes of the zeolite framework. In addition to elucidating relaxation from MW-heated states, we expect that these studies will shed light on energy transfer in other contexts, such as adsorption and reaction in zeolites, which often involve significant heat release.

A Study of Pore Blockage in Silicalite Zeolite Using Free Energy Perturbation Calculations

Binary systems consisting of large coadsorbed molecules (n-hexane, cyclohexane, and benzene) with smaller penetrant molecules (methane) were simulated to investigate the mechanisms of pore blockage in the zeolite silicalite. Benzene and cyclohexane trap the methane molecules in the zeolite channels on the time scales of molecular dynamics simulations. Minimum energy paths for methane diffusion past the blocking molecules were determined, and free energy perturbation calculations were carried out along the paths to get the rate constants of methane hopping past coadsorbed benzene and cyclohexane molecules, which adsorb in the channel intersections. Three principal diffusion pathways were found in both the methane/benzene and methane/cyclohexane systems. Minima which were connected by low-energy pathways were grouped together into macrostates. Using the calculated hopping rates between macrostates, kinetic Monte Carlo was then used to obtain the diffusivity of methane with a coadsorbate benzene loading such that all channel intersections are filled by benzene - conditions where molecular dynamics simulations fail. Passage of methane across cyclohexane molecules involved pushing the cyclohexane molecules into the channels from their preferred channel intersection positions.