Journal of Chemistry

Volume 2013 (2013), Article ID 545367, 7 pages

http://dx.doi.org/10.1155/2013/545367

## Molecular-Dynamics Simulation of Self-Diffusion of Molecular Hydrogen in X-Type Zeolite

School of Materials Science and Engineering, Shenyang Ligong University, Shenyang 110159, China

Received 3 May 2013; Revised 15 July 2013; Accepted 17 July 2013

Academic Editor: KRS Chandrakumar

Copyright © 2013 Xiaoming Du. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The self-diffusion of hydrogen in NaX zeolite has been studied by molecular-dynamics simulations for various temperatures and pressures. The results indicate that in the temperature range of 77–293 K and the pressure range of 10–2700 kPa, the self-diffusion coefficients are found to range from 1.61 × 10^{−9} m^{2}*·*s^{−1} to 3.66 × 10^{−8} m^{2}*·*s^{−1} which are in good agreement with the experimental values from the quasielastic neutron scattering (QENS) and pulse field gradients nuclear magnetic resonance (PFG NMR) measurements. The self-diffusion coefficients decrease with increasing pressure due to packing of sorbate-sorbate molecules which causes frequent collusion among hydrogen molecules in pores and increase with increasing temperature because increasing the kinetic energy of the gas molecules enlarges the mean free path of gas molecule. The activated energy for hydrogen diffusion determined from the simulation is pressure-dependent.

#### 1. Introduction

Zeolite having the faujasite (FAU) framework structure type is a microporous material with regular crystalline structure. It is widely employed as molecular sieves in various industrial processes, such as adsorption, catalysis, purification, and separation of mixtures [1, 2]. The synthetic low silica FAU zeolites, namely, X-type zeolites, have the property that the diameter of the pores can be controlled by changing the size and charge of the exchangeable cations. This offers a variable internal surface area, which is a useful attribute for the physisorption of gases.

Recently, zeolites are attractive materials for the safe and efficient storage of hydrogen [3–5]. One of the potential advantages of using zeolites as hydrogen storage materials is their thermal and chemical stability, which has been shown to be useful in catalytic applications. This is exemplified by their lack of combustibility both in air and, in particular, in a hydrogen atmosphere. The hydrogen storage capacities of X-type zeolites have been reported experimentally and theoretically [6–9]. However, little attention was previously paid to the transport properties of adsorbed hydrogen molecules which play a central role in catalytic and separation processes [10] that take place within zeolite cavities [11]. Understanding the host-guest interactions that control molecular diffusion may suggest new materials with advanced performance. We have recently calculated [12] self-diffusion coefficients for hydrogen in all-silica ZSM-5 over the temperature range of 77–293 K and the pressure range of 14–3360 kPa. The self-diffusion coefficients are found to range from 1.2 × 10^{−9} m^{2}/s to 3.8 × 10^{−7} m^{2}/s. The diffusion of the hydrogen in all-silica ZSM-5 zeolite channels is anisotropic. Fu et al. [13] have studied the diffusion of molecular hydrogen adsorbed in zeolite 13X at high coverage using quasielastic neutron scattering (QENS) for temperatures ranging from 10 to 60 K. The diffusive motion of the adsorbed H_{2} could be described by a liquid-like jump diffusion model above 35 K. DeWall et al. [14] also obtained similar results to Fu et al. Bar et al. [15] reported the results that the self-diffusion coefficients of hydrogen in NaX, NaA, and CHA zeolites decreased with the decrease of the pore diameter for these zeolites using QENS and NMR techniques. Van Den Berg et al. [16] investigated the influence of the cation distribution in sodalite zeolite on the hydrogen self-diffusion using molecular-dynamics simulation; they found that the hydrogen diffusion coefficient was much larger for the ordered cation distribution than for the disordered distribution. However, very little work on diffusion behavior of hydrogen in FAU zeolites as a function of loading was done so far.

In this paper, we investigate the effect of temperature and pressure on self-diffusivities determined using molecular-dynamics method for hydrogen molecules adsorbed in X-type zeolite. All results are for diffusion over the temperature range of 77–293 K and the pressure range of 10–2700 kPa. Computational data on the diffusivities of hydrogen in zeolites are presented and compared with the experimental data.

#### 2. Model and Methods

Structure of initial unit cell of NaX zeolite (Na_{88}Si_{104}Al_{88}O_{384}) was constructed according to the result of X-ray crystallographic studies [17]. The crystal is cubic with unit cell parameters Å and space group (number 203). The pore topology of the NaX zeolite is shown schematically in Figure 1. NaX is a 3-dimensional zeolite and has a 12-membered ring pore opening of 7.4 Å in diameter providing access to a supercage of 12.4 Å in diameter. In dehydrated NaX zeolite, there are four Na^{+} cation sites. 2.7, 26.6, 29.3, and 29.3 Na^{+} cations are located at the site I in (the hexagonal prism), I′ (in the -cage, close to the hexagonal window to the hexagonal prism), site II (located in the sixth ring of the -cage), and site III (in the supercage, close to a square window between two other square windows), respectively [17].

All MD simulations have been performed using the simulation program code DL_POLY [18]. DL_POLY is a package of molecular simulation routines written by Smith and Forester, copyright CCLRC for the Central Laboratory of the Research Councils, Daresbury Laboratory at Daresbury, Nr. Warrington (1996). Lennard-Jones potentials are used for adsorbate-adsorbate and adsorbate-framework interactions. The interactions between hydrogen and the oxygen, silicon, and aluminum atoms of the framework and the extraframework cations are derived and used in [19–21]. The values of the potential parameters are summarized in Table 1. The Lorenz-Berthelot mixing rules are applied to obtain the cross potentials and .

The hydrogen molecule is treated as a centrosymmetrical Lennard-Jones particle, which has proven to be a valid and accurate approximation with respect to extended representation in other studies [22, 23]. To define the zeolite structure by molecular mechanics, most simulation studies are performed using the method proposed by Bezus et al. [24]. In a Kiselev-type potential the atoms from the zeolite framework are fixed at their crystallographic position. In our calculations, the zeolite framework was kept rigid and the nonframework cations (Na^{+}) were allowed to move freely in the zeolite. In order to get better starting structures for the barrier optimizations, we performed energy minimization calculations to determine the equilibrium structures using density function theory (DFT) methods that were successfully used to explore hypothetical FAU and LTA with dense topologies. DFT calculations were carried out with the code Quantum ESPRESSO [25] which uses atom centered basis functions that are particularly efficient for total energy studies of very low density materials with very large unit cells. Periodic DFT calculations were carried out using the PBE exchange correlation functional within the GGA approximation. The minimized zeolites are used as adsorbent for hydrogen adsorption.

For the simulations, a simulation box with one-unit cell was used. Periodic boundary conditions were applied in three dimensions in order to simulate an infinite system. The canonical (NVT) ensemble was used, with a step size of 0.001 ps, at 77 K, 160 K, 195 K, 230 K, and 293 K, using a Nosé-Hoover thermostat [26]. The initial positions for the hydrogen adsorbate were determined from grand canonical Monte Carlo (GCMC) adsorption simulations. The loading of each starting configuration was determined from adsorption isotherms at 77 K, 160 K, 195 K, 230 K, and 293 K and 10–3000 kPa [27] to calculate the self-diffusion coefficients of hydrogen inside NaX zeolite. The starting configuration was equilibrated for 50 ps before collecting averages. Molecular-dynamics simulations were run for 500 ps. The longer simulation runs yielded results that were not significantly different from the reported results. The nonbond interaction cut-off distance was chosen to be 13 Å which is less than half of the length of the simulation box. Test runs with a larger cut-off distance did not yield significantly different results.

The primary means of tracking the movement of adsorbate molecules over the course of the simulation was via the mean-squared displacement. The overall mean-squared displacement is where is the position of molecule at time and is an arbitrary starting point.

The mean-squared displacement results were used to calculate self-diffusion coefficients for hydrogen in each simulation using the Einstein equation. The Einstein equation for the self-diffusion coefficient in three dimensions is [28] where is the diffusion coefficient of the molecules and is the number of the diffusing molecules. Furthermore, the diffusion coefficient components , , and in -, -, and -direction can be calculated by replacing with , , and in (2), respectively.

The radial distribution function gives a measure of the probability that, given the presence of an atom (molecule, colloid, etc.) at the origin of an arbitrary reference frame, there will be an atom with its center located in a spherical shell of infinitesimal thickness at a distance, , from the reference atom. This concept also embraces the idea that the atom at the origin and the atom at distance may be of different chemical types, say and . The resulting function is then commonly given the symbol and is defined by Hansen and McDonald [29]: where and represent the numbers of and particles, respectively, is the model volume, and represents the number of particles in the range from to at the origin of th particle.

#### 3. Results and Discussion

Figure 2 displays the mean-squared displacement (MSD) of the hydrogen molecule of diffusion in NaX zeolite at 77 K for different pressures which are presented as a function of time. Figure 3 displays the same data at 10 kPa for different temperatures. The relationship between the observed MSD and time is linear to a very good approximation. It indicates that the normal diffusion occurs at this time scale. The self-diffusion coefficients were obtained from the slope of the curves via a linear regression fit to linear portion of MSD curves to (2). The calculated self-diffusion coefficient for hydrogen at 77 K for pressures of 10 and 100 kPa is 5.87 × 10^{−}9 and 1.90 × 10^{−}9 m^{2}·s^{−1}, respectively (Figure 2). These values are in good agreement with reported experimental diffusion coefficients obtained by Bar et al. [15] from the quasielastic neutron scattering (QENS) and pulse field gradients nuclear magnetic resonance (PFG NMR) methods. They are also in good agreement with those obtained from previous MD simulations [30]. It can be seen from Figure 2 that as hydrogen pressure increases there is an overall decrease in the diffusion coefficients for all temperatures. That is, the diffusivity decreases with increasing loading. Thus with each successive increase in pressure there is less available space for the fast moving hydrogen molecules; hence, the movement becomes hindered. At low pressure, the molecules have the highest mobility. NaX zeolite has a 12-membered ring pore opening of 7.4 Å in diameter that provides access to a supercage of 12.4 Å in diameter. Therefore, molecules have the freedom to move in all directions from one supercage to another cage (intercage). Due to their small size, hydrogen molecules have the good ability to move from one cage to another. At high pressure, we observe reduced translational mobility due to a reduced volume for the individual molecule. The diffusivity also decreases with increasing pressure. The sorbate-sorbate interaction further contributes to the decreased mobility with increasing pressure. Similar results are also found for alkane diffusion in mordenite zeolite [31].

In general as the simulation temperature increases there is an overall increase in diffusion coefficients for the hydrogen molecules. With comparison of the values of diffusion coefficients at a certain pressure and different temperatures as shown in Figure 4, the diffusion coefficients increase with the increase in temperature. This shows that the diffusion of hydrogen in NaX zeolite is a thermally activated process. Figure 5 shows the relationship of diffusion coefficients changing with the inverse of the temperature at different pressure. It appears that the diffusion coefficients only for 10 and 100 kPa versus follow a linear relation. The activation energy for the diffusion of hydrogen in NaX is calculated by fitting diffusion coefficients for 10 and 100 kPa versus (Figure 5) to the Arrhenius equation [32]. The values obtained are 2.01 and 1.95 kJ/mol for 10 and 100 kPa, respectively. The activation energy of hydrogen in zeolite NaX was found experimentally to be 1–4.0 kJ/mol [15], and it was 2 kJ/mol for NaA zeolite [33]. Schuring et al. [31] considered that the activated energy is possibly dependent on pressure of gas molecules due to the increase of diffusion coefficient with the increase in temperature and the decrease in pressure in studying the diffusion of alkene in MFI zeolite. To solve the question, they proposed that the activated energy should be calculated at low pressure.

In order to investigate the effects of pressure and temperature on the distance among hydrogen molecules, the radial distribution functions (RDFs) of the centers of mass for hydrogen molecules in zeolite pores at different pressures and temperatures are shown in Figures 6 and 7, respectively. Only one peak is observed which is centered around 2.5 Å–3.3 Å in RDF plots of the hydrogen molecules, indicating the adsorbed hydrogen in zeolite pores is still unordered. For 77 K, at the pressures of 10 kPa and 100 kPa the distances are both about 2.9 Å which is approximate to the dynamics diameter of hydrogen molecule, 2.958 Å. The probability of the collusion among hydrogen molecules in pores is less than that of the collusion between hydrogen molecules and pore walls. The diffusion process of hydrogen molecules is controlled from the collusion between hydrogen molecules and pore walls. The distance decreases to about 2.7 Å and 2.5 Å for 600 kPa and above 1700 kPa, respectively, as shown in Figure 6. It is indicated that at high pressure hydrogen molecules are induced to change orientation or conformation to achieve a more compact packing. The collusion among hydrogen molecules in pores becomes frequent. It can be argued that this is also the reason for the very slow diffusion in the case of higher pressure. For the effect of temperature, in the case of 10 kPa, as shown in Figure 7, the distances among hydrogen molecules are about 2.9 Å at the temperatures of 77 K, 160 K, and 195 K. This increases to about 3.1 Å for 230 K and 3.3 Å for 293 K, respectively. The reason is that the kinetic energy of the gas molecules increases with increasing temperature. This can increase the mean free path of gas molecule, showing the very quick diffusion in the case of higher temperature.

To obtain a more detailed investigation of the information on the short-time dynamics of the system, the velocity autocorrelation function (VACF) of the diffusion of hydrogen in zeolite is analyzed. The velocity autocorrelation function is defined as where and denote the velocity of particle at time and time zero (i.e., ), respectively. Again, averaging is performed over all time origins in the trajectory.

Figures 8 and 9 show the velocity autocorrelation functions of the diffusion of hydrogen in NaX zeolite for different pressures and temperatures, respectively. In general, the structure of the velocity autocorrelation function can provide information on the short and intermediate time dynamics of the system. To quantify this short-time dynamical behavior of hydrogen molecules, we consider two quantities: (i) The time, , at which first turns negative represents the average time at which the sorbate first encounters a repulsive barrier and is often closely related to the crossover time from ballistic to diffusion motion; (ii) , the position of the first minimum in indicates the time scale on which a particle reverses its velocity as a result of repulsive encounters. As collisions with the wall and other sorbates increase, the correlation function decays to zero. The values of and in Figures 8 and 9 are listed in Tables 2 and 3. The magnitude of provides an estimate of the probable time, on average, that a molecule first encounters a repulsive barrier. From Table 2, we see that as pressure increases, continuously decreases. This may be attributed to the increasing frequency of intermolecular collision and collision between molecules and the walls of the pore. This trend in is reversed with increasing temperature (Table 3). The reason is possibly that as temperature increases, the dynamical energy of a molecule increases, which shortens the crossover time from ballistic to diffusion motion. Also seen from Tables 2 and 3, the values of are always larger than that of . However, values do not exhibit the good consistence with the variation of pressure and temperature.

#### 4. Conclusions

The diffusion behavior of molecular hydrogen in NaX zeolite was investigated by using molecular-dynamics (MD) simulations. Our molecular-dynamics results show that the self-diffusion coefficient of hydrogen in NaX decreases as the pressure increases; this is in agreement with the earlier quasielastic neutron scattering (QENS) and pulse field gradients nuclear magnetic resonance (PFG NMR) intracrystalline self-diffusivity measurements of the other group [15]. The reason is found from RDF of hydrogen molecules that at high pressure hydrogen molecules are induced to change orientation or conformation to achieve a more compact packing. The collusion among hydrogen molecules in pores becomes frequent, showing the very slow diffusion in the case of the higher pressure. The temperature dependence of the self-diffusion coefficient shows an increase with increasing temperature between 77 K and 293 K. An explanation for the effect is that the kinetic energy of the gas molecules increases with increasing temperature. This can increase the mean free path of gas molecule, showing the very quick diffusion in the case of higher temperature. The activated energy for hydrogen diffusion determined from the simulation is pressure-dependent. The value obtained below 100 kPa is about 2 kJ/mol. The results and data of hydrogen diffusion properties obtained from the simulations are theoretically significant for understanding of the mechanism of hydrogen storage on microporous zeolites.

#### Acknowledgment

The authors gratefully acknowledge the financial support from the program for Liaoning Excellent Talents in University (LNET), China (LJQ2012016).

#### References

- R. M. Barrer,
*Zeolites and Clay Minerals as Sorbents and Molecular Sieves*, Academic Press, New York, NY, USA, 1978. - J. Kärger and D. M. Ruthven,
*Diffusion in Zeolites and other Microporous Solids*, Wiley Interscience, New York, NY, USA, 1992. - A. Züttel, P. Sudan, P. Mauron, T. Kiyobayashi, C. Emmenegger, and L. Schlapbach, “Hydrogen storage in carbon nanostructures,”
*International Journal of Hydrogen Energy*, vol. 27, no. 2, pp. 203–212, 2002. View at Publisher · View at Google Scholar · View at Scopus - A. W. C. van den Berg and C. O. Areán, “Materials for hydrogen storage: current research trends and perspectives,”
*Chemical Communications*, no. 6, pp. 668–681, 2008. View at Publisher · View at Google Scholar - X. Lin, J. Jia, X. Zhao et al., “High H
_{2}adsorption by coordination-framework materials,”*Angewandte Chemie*, vol. 45, no. 44, pp. 7358–7364, 2006. View at Publisher · View at Google Scholar · View at Scopus - H. W. Langmi, A. Walton, M. M. Al-Mamouri et al., “Hydrogen adsorption in zeolites A, X, Y and RHO,”
*Journal of Alloys and Compounds*, vol. 356-357, pp. 710–715, 2003. View at Publisher · View at Google Scholar · View at Scopus - Y. Li and R. T. Yang, “Hydrogen storage in low silica type X zeolites,”
*Journal of Physical Chemistry B*, vol. 110, no. 34, pp. 17175–17181, 2006. View at Publisher · View at Google Scholar · View at Scopus - K. P. Prasanth, R. S. Pillai, S. A. Peter et al., “Hydrogen uptake in palladium and ruthenium exchanged zeolite X,”
*Journal of Alloys and Compounds*, vol. 466, no. 1-2, pp. 439–446, 2008. View at Publisher · View at Google Scholar · View at Scopus - L. Kang, W. Deng, K. Han, T. Zhang, and Z. Liu, “A DFT study of adsorption hydrogen on the Li-FAU zeolite,”
*International Journal of Hydrogen Energy*, vol. 33, no. 1, pp. 105–110, 2008. View at Publisher · View at Google Scholar · View at Scopus - J. Weitkamp, “New directions in zeolite Catalysis,” in
*Catalysis and Adsorption By Zeolites*, G. Olhmann, J. C. Vedrine, and P. A. Jacobs, Eds., pp. 21–45, Elsevier, Amsterdam, The Netherlands, 1991. View at Google Scholar - J. M. Newsam, “Zeolites,” in
*Solid State Chemistry: Compounds*, A. K. Cheetham and P. Day, Eds., pp. 234–280, Oxford University Press, Oxford, UK, 1992. View at Google Scholar - X. M. Du, Y. Huang, and E. D. Wu, “Molecular-dynamics studies of the diffusion of H
_{2}in all-silica ZSM-5,”*Materials Science Forum*, vol. 704-705, pp. 401–406, 2012. View at Google Scholar - H. Fu, F. Trouw, and P. E. Sokol, “A quasi-elastic and inelastic neutron scattering study of H
_{2}in zeolite,”*Journal of Low Temperature Physics*, vol. 116, no. 3-4, pp. 149–165, 1999. View at Google Scholar · View at Scopus - J. DeWall, R. M. Dimeo, and P. E. Sokol, “Slow diffusion of molecular hydrogen in zeolite 13X,”
*Journal of Low Temperature Physics*, vol. 129, no. 3-4, pp. 171–184, 2002. View at Publisher · View at Google Scholar · View at Scopus - N. K. Bar, H. Ernst, H. Jobic, and J. Karger, “Combined quasi-elastic neutron scattering and NMR study of hydrogen diffusion in zeolites,”
*Magnetic Resonance in Chemistry*, vol. 37, pp. S79–S83, 1999. View at Google Scholar - A. W. C. Van Den Berg, S. T. Bromley, E. Flikkema, and J. C. Jansen, “Effect of cation distribution on self-diffusion of molecular hydrogen in Na
_{3}Al_{3}Si_{3}O_{12}sodalite: a molecular dynamics study,”*Journal of Chemical Physics*, vol. 121, no. 20, pp. 10209–10216, 2004. View at Publisher · View at Google Scholar · View at Scopus - D. H. Olson, “The crystal structure of dehydrated NaX,”
*Zeolites*, vol. 15, no. 5, pp. 439–443, 1995. View at Google Scholar · View at Scopus - T. R. Forester and W. Smith,
*The DL_POLY_2 User Manual*, CCLRC Daresbury Laboratory, 1997. - A. Michels, W. de Graaff, and C. A. Ten Seldam, “Virial coefficients of hydrogen and deuterium at temperatures between −175°C and +150°C. Conclusions from the second virial coefficient with regards to the intermolecular potential,”
*Physica*, vol. 26, no. 6, pp. 393–408, 1960. View at Google Scholar · View at Scopus - K. Watanabe, N. Austin, and M. R. Stapleton, “Investigation of the air separation properties of zeolites types A, X and Y by monte carlo simulations,”
*Molecular Simulation*, vol. 15, no. 4, pp. 197–221, 1995. View at Publisher · View at Google Scholar - G. Maurin, P. Llewellyn, T. Poyet, and B. Kuchta, “Influence of extra-framework cations on the adsorption properties of X-faujasite systems: microcalorimetry and molecular simulations,”
*Journal of Physical Chemistry B*, vol. 109, no. 1, pp. 125–129, 2005. View at Publisher · View at Google Scholar · View at Scopus - A. W. C. Van Den Berg, S. T. Bromley, N. Ramsahye, and T. Maschmeyer, “Diffusion of molecular hydrogen through porous materials: the importance of framework flexibility,”
*Journal of Physical Chemistry B*, vol. 108, no. 16, pp. 5088–5094, 2004. View at Google Scholar · View at Scopus - A. I. Skoulidas, D. M. Ackerman, J. K. Johnson, and D. S. Sholl, “Rapid transport of gases in carbon nanotubes,”
*Physical Review Letters*, vol. 89, no. 18, pp. 185901–185904, 2002. View at Google Scholar · View at Scopus - A. G. Bezus, A. V. Kiselev, A. A. Lopatkin, and P. Q. Du, “Molecular statistical calculation of the thermodynamic adsorption characteristics of zeolites using the atom-atom approximation. Part 1.-Adsorption of methane by zeolite NaX,”
*Journal of the Chemical Society Faraday Transactions 2*, vol. 74, pp. 367–379, 1978. View at Google Scholar - http://www.pwscf.org.
- D. Frenkel and B. Smit,
*Understanding Molecular Simulation*, Academic Press, San Diego, Calif, USA, 1996. - X. M. Du, Y. H, Q. Zhang, J. Li, and E. D. Wu,
*Acta Petrolei Sinica Suppl*, vol. 1, pp. 37–40, 2012. - D. Hofmann, L. Fritz, C. Schepers, and M. Bohning, “Detailed-atomistic molecular modeling of small molecule diffusion and solution processes in polymeric membrane materials,”
*Macromolecular Theory and Simulations*, vol. 9, no. 6, pp. 293–327, 2000. View at Publisher · View at Google Scholar - J. P. Hansen and I. R. McDonald,
*Theory of Simple Liquids*, Academic Press, London, UK, 2nd edition, 1990. - G. K. Papadopoulos and D. N. Theodorou, “Simulation studies of methane, carbon dioxide, hydrogen and deuterium in ITQ-1 and NaX zeolites,”
*Molecular Simulation*, vol. 35, no. 1-2, pp. 79–89, 2009. View at Publisher · View at Google Scholar · View at Scopus - D. Schuring, A. P. J. Jansen, and R. A. Van Santen, “Concentration and chainlength dependence of the diffusivity of alkanes in zeolites studied with MD simulations,”
*Journal of Physical Chemistry B*, vol. 104, no. 5, pp. 941–948, 2000. View at Google Scholar · View at Scopus - D. Do,
*Adsorption Analysis: Equilibrium and Kinetics*, Imperial College Press, London, UK, 1998. - R. Kahn, E. Cohen De Lara, and E. Viennet, “Diffusivity of the hydrogen molecule sorbed in NaA zeolite by a neutron scattering experiment,”
*The Journal of Chemical Physics*, vol. 91, no. 8, pp. 5097–5102, 1989. View at Google Scholar · View at Scopus