A systematic investigation of monatomic spherical sorbates in the supercages of zeolites Y and A by molecular dynamics technique is presented. Rates of intercage diffusion, rates of cage visits, and the diffusion coefficients have been calculated as a function of the sorbatezeolite interaction strength. These properties exhibit markedly different dependences on interaction strength for the two zeolites. The observed behavior is shown to be a consequence of the two principal mechanisms of intercage diffusion and the energetic barrier associated with them. The diffusion coefficient and other properties associated with intercage diffusion are found to be directly proportional to the reciprocal of the square of the sorbate diameter when the sorbate diameter is significantly smaller than the window diameter. As the sorbate diameter increases, a peak is observed in all the transport properties investigated including the diffusion coefficient. We call this surprising effect as the ring or levitation effect and it explains several anomalous results reported in the literature and suggests a breakdown of the geometrical criterion for diffusion of sorbates. It shows that under certain conditions nongeometrical factors play a major role and geometrical factors become secondary in the determination of the molecular sieve property. A generalized parameter has been proposed which suggests conditions under which one can expect the ring or levitation effect in any porous medium. Inverse size selectivity becomes operative under these conditions.

The nature of the adsorption sites and the mobility of benzene in sodium-Y zeolite around room temperature is investigated by using molecular dynamics calculations carried out at a loading of two molecules per cage. The motion of benzene is shown to be highly anisotropic, being confined largely to the cage walls. The simulation results correlate well with recent neutron diffraction and infrared absorption data.

Solid state ionic conductors are important from an industrial viewpoint. A variety of such conductors have been found. In order to understand the reasons for high ionic conductivity in these solids, there have been a number of experimental, theoretical and computational studies in the literature. We provide here a survey of these investigations with focus on what is known and elaborate on issues that still remain unresolved. Conductivity depends on a number of factors such as presence of interstitial sites, ion size, temperature, crystal structure etc. We discuss the recent results from atomistic computer simulations on the dependence of conductivity in NASICONs as a function of composition, temperature, phase change and cation among others. A new potential for modelling of NASICON structure that has been proposed is also discussed.

We report a molecular dynamics study of a binary mixture consisting of a large (host) particle and a smaller (guest) particle whose radius is varied over a range. These simulations investigate the possible existence of a diffusion anomaly or levitation effect in dense fluids, previously seen for guest molecules diffusing within porous solids. The voids in the larger component have been characterized in terms of void and neck distributions by means of Voronoi polyhedral analysis. Four different mixtures with differing ratios of guest to host diffusivities (D) have been studied. The results suggest that the diffusion anomaly is seen in both close-packed solids with disorder and dense fluids. In the latter, the void network is constantly and dynamically changing and possesses a considerable degree of disorder. The two regimes, viz., the linear regime (LR) and the anomalous regime (AR), found for porous solids are shown to exist for a dense medium as well. The linear regime is characterized by D(g) proportional to 1/sigma(gg)(2), where sigma(gg) is the diameter of the guest. The anomalous regime exhibits a maximum in D up to rather high temperatures (T = 1.663), even though in porous solids the maximum disappears at higher temperatures. In agreement with previous studies on porous solids, a particle in the AR is associated with lower activation energy, lower friction, and less backscattering in the velocity autocorrelation function when compared to a particle in the LR. Wavevector dependent self-diffusivity, Delta, and decay of the intermediate scattering function, F(s)(k, t), exhibit contrasting behaviors for the LR and AR. For LR, Delta exhibits a minimum at values of k at which there are spatial correlations in S(k) while a smooth decrease with k is seen for AR. For LR, F(s)(k, t) shows a biexponential decay corresponding to two different time scales of motion. Probably, the fast decay is associated with motion within the first shell of solvent neighbors and the slow decay with motion past these shells. For AR, a single-exponential decay is seen. The results indicate a breakdown of the Stokes-Einstein (SE) relationship. The relevant quantity that determines the validity of the SE relationship is the levitation parameter which is indirectly related to the solute/solvent radius ratio and not either the size of the solute or the solvent alone.

Long molecular dynamics simulations (120 ns) of guests confined in zeolite NaCaA as well as that of pure fluid have been carried out in the microcanonical ensemble to obtain an estimate of the error associated with the diffusion coefficient (D). It is found that the error is about 50% for a 1 ns run for argon parameters in NaCaA zeolite. It is found that shorter runs lead to an overestimation of the value of D. It is also found that a linear fit to a region of the msd (mean squared displacement) below ∼200 ps leads to an overestimation of D. For runs that are sufficiently long, the initial configuration has no effect on D. Further, the values obtained from the msd and the velocity autocorrelation function converge for the runs reported here. It is shown that the error in D depends on the nature of the system and its underlying potential energy surface. The calculated statistical inefficiency suggests that the averages of properties such as the total energy will be uncorrelated over blocks longer than 150 ps. Dependence of the error in D on the simulation length and size is reported.

Self-diffusivity, D, of diffusants in widely differing mediums such as liquids (e.g., solution), porous solids (e.g., guests in zeolites), or ions in polar solvents exhibit strong size dependence. We discuss the nature of the size dependence observed in these systems. Altogether, different theoretical approaches have been proposed to understand the nature of size dependence of D not only across these widely differing systems but even in just one medium or class of systems such as, for example, ions in polar solvents. But molecular dynamics investigations in the past decade have shown that the size dependence of self-diffusion in guest-porous solids could have origins in the mutual cancellation of forces that occurs when the size of the diffusant is comparable to the size of the void. The effect leading to the maximum in D is known as the levitation effect (LE). Such a cancellation is a consequence of symmetry. This effect exists in all porous solids irrespective of the geometrical and topological details of the pore network provided by the solid. Recent studies show that the levitation effect and size-dependent diffusivity maximum exists for uncharged solutes in solvents. One of the consequences of this is the breakdown in the Stokes-Einstein relationship over a certain range of solute-solvent size ratio. Experimental measurements of ionic conductivity over the past hundred years have found the existence of a size-dependent diffusivity maximum leading to violation of the Walden's rule for ions in polar solvents. Molecular dynamics simulations and experimental data suggest that even this maximum has its origin in LE. Simulation studies of impurity atom diffusion in close-packed solids as well as ions in superionic and other solids suggest the existence of a size-dependent diffusivity maximum in these materials as well. The levitation effect is a universal effect leading to a maximum in diffusivity of a diffusant in a variety of condensed matter phases. The only condition for its existence appears to be the presence of van der Waals or electrostatic interactions.

No abstract

We propose an alternative explanation of the size-dependent maximum in ion mobility in water in terms of the levitation effect, which accounts for the observed size-dependent maximum in the mobility of guest diffusion in porous media. In this explanation, the size at which the maximum occurs is related to the structure of the void space of the water; at the mobility maximum, the diffusant passes smoothly through necks connecting voids, and its potential energy shows minimum fluctuations. Molecular dynamics simulations of charged spheres of varying sizes are used to support this hypothesis. As in the levitation effect, the friction coefficient, the potential energy fluctuations, and the activation energy are found to be minima for particles with maximum self diffusivities similar to the guest diffusion in zeolites. Wavelength-dependent self diffusivities indicate a monotonic and oscillatory dependence, respectively, on wavenumber k for anomalous (AR) and linear regimes (LR). These are associated with single and biexponential decay of the incoherent intermediate scattering function.