Molecular Dynamics Studies of Hydrofluorocarbons in Faujasite-type Zeolites:  Modeling Guest-Induced Cation Migration in Dry Zeolites

Abstract

We have developed and applied a new force field for simultaneously modeling the dynamics of hydrofluorocarbons (HFCs) and exchangeable Na cations in faujasite-type zeolites. Our aim is to account for (i) the zeolite's capacity of separating HFC isomers, (ii) the experimentally observed unusual cation migration in Na−Y from the β-cages into the supercages upon the adsorption of HFCs, and (iii) the abnormal trans/gauche ratio in these systems. Energy minimizations and molecular dynamics simulations performed with this force field give excellent agreement with experimental data on heats of adsorption, guest−host distances, infrared spectra, and conformer ratios for different coverages of HFC-134 (CF2H−CF2H) and HFC-134a (CH2F−CF3) in Na−X (Si:Al = 1.2) and Na−Y (Si:Al = 2.4). The force field also accounts partially for the observed cation migration at intermediate loadings and fully at high loadings. The extent of migration is found to be influenced by the competition among Na−O, Na−F, and Na−Na interactions. The Na−O interaction disfavors migration because Na(I‘) sites are found to be energetically more favorable than Na(III‘) sites; the Na−F attraction obviously favors migration; and surprisingly, the Na−Na repulsion also favors migration because moving Na cations into supercages leads to better cation dispersion. This migration occurs in a two-step mechanism that involves first a concerted two-cation jump, S I‘ → S II‘/S II → (S III or S III‘), followed by a S II‘ → S II jump, leading to a net process S I‘ → (S III or S III‘). The preferred binding site in both Na−X and Na−Y involves HFCs anchored by both site II and site III‘ cations. The loading dependence of the heat of HFC adsorption in zeolite Na−X is predicted to be different from that in Na−Y, because of the energetics of cation migration in Na−Y. HFC-134 is generally more strongly bound to both zeolites, because of its ability to make close Na−F and O−H contacts with the zeolites. The binding energy for the gauche conformer of HFC-134 is larger than that for trans at low loadings, but as loading increases, the difference decreases. The highly correlated small-amplitude motion predicted for cations in bare faujasites is quenched upon adsorption of HFCs. Most of the HFCs are too strongly bound to exhibit diffusive behavior during our molecular dynamics simulations.

FT-IR study of the adsorption of carbon monoxide and of some nitriles on Na-faujasites: Additional insight on the formation of complex interactions

“…On the other hand, the possibility of weak interactions of the C-H bonds of nitriles with other oxygen atoms exposed in the cavity (like proposed to occur for hydrocholorocarbons in Nafaujasites [44]), should not cause a decrease in the CN stretching frequency. Taking into account that the cations may displace a little bit upon adsorption of different molecules, the possibility of the formation of adsorbed species where two different cations interact with the N lone pair and with a p-type orbital seems to be possible on NaX.…”

A Numerical Evidence for Nonframework Cation Redistribution Upon Water Adsorption in Faujasite Zeolite

“…In a previous molecular dynamics (MD) study devoted to hydrated zeolite 4A, Faux [16] observed a relatively slow motion of some of the Na + ions (one order of magnitude slower than the water diffusion in the same system). Using the same type of method, Jaramillo et al [17] also observed the diffusion of some of the cations in a model of NaX upon adsorption of two polar hydrofluorocarbons. Standard molecular dynamics is obviously not well-suited for studying equilibrium cation distributions (and redistributions), since cation diffusion takes place on a timescale which is larger (or comparable) to the MD timescale.…”

Adsorption of Water and Aromatics in Faujasite Zeolites: A Molecular Simulation Study

“…The cation force field has been adapted from the work of Jaramillo and Auerbach (Jaramillo et al, 2001). The cation-framework potential consists of an exp-6 repulsion-dispersion term that acts between the cation and the oxygen atoms of the faujasite and a Coulombic term that acts between the cation and both the oxygen and T atoms of the framework.…”