Diffusion of sorbates in zeolites Y and A novel dependence on sorbate size and strength of sor ba t e-Zeoli t e interaction

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 sorbate-zeolite 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.

Dynamics of zeolite cage and its effect on the diffusion properties of sorbate: persistence of diffusion anomaly in NaA zeolite

Recent computer simulations on zeolites Y and A have found that the diffusion coefficient and the rate of intercage diffusion exhibit, apart from a linear dependence on the reciprocal of the square of the sorbate diameter, an anomalous peak as sorbate diameter approaches the window diameter. Here we report molecular dynamics simulations of zeolite NaA incorporating framework flexibility as a function of sorbate diameter in order to verify the existence of anomalous diffusion. Results suggest persistence of anomalous diffusion or ring effect. This suggests that the anomalous behavior is a general effect characteristic of zeolites Y and A. The barrier for diffusion across the eight-ring window is seen to be negative and is found to decrease with sorbate size. The effect of sorbate on the cage motion has also been investigated. Results suggest that the window expands during intercage migration only if the sorbate size is comparable to the window diameter. Flexible cage simulations yield a higher value for the diffusion coefficient and also the rate of intercage diffusion. This increase has been shown to be due to an increase in the intercage diffusions via the centralized diffusion mode rather than the surface-mediated mode. It is shown that this increase arises from an increase in the single particle density distribution in the region near the cage center.

Diffusion Anomaly in Silicalite and VPI-5 from Molecular Dynamics Simulations

Molecular dynamics (MD) simulations on rigid and flexible framework models of silicalite and a rigid framework model of the aluminophosphate VPI-5 for different sorbate diameters are reported. The sorbate-host interactions are modeled in terms of simple atom-atom Lennard-Jones interactions. The results suggest that the diffusion coefficient exhibits an anomaly as gamma approaches unity. The MD results confirm the existence of a linear regime for sorbate diameters significantly smaller than the channel diameter and an anomalous regime observed for sorbate diameters comparable to the channel diameter. The power spectra obtained by Fourier transformation of the velocity autocorrelation function indicate that there is an increase in the intensity of the low-frequency component for the velocity component parallel to the direction of motion for the sorbate diameter in the anomalous regime. The present results suggest that the diffusion anomaly is observed irrespective of (1) the geometry and topology of the pore structure and (2) the nature of the host material. The results are compared with the work of Derouane and co-workers, who have suggested the existence of ''floating molecules'' on the basis of earlier theoretical and computational approaches.

Decoupling of diffusion from viscosity: Difference scenario for translational and rotational motions

Recent experiments have indicated a dramatically different viscosity dependence of the translational and the rotational diffusion coefficients in a supercooled liquid as the glass transition temperature is approached from above. While the translational motion seems to be decoupled from the rising viscosity (eta), the rotational motion seems to remain firmly coupled to eta. In order to understand the microscopic origin of this behavior, we have carried nut detailed theoretical calculations of both the quantities by using a self-consistent mode-coupling theory (MCT). it is found that when the size of the solute is same as that of the solvent molecules, the conventional MCT fails to predict the observed decoupling. The solvent inhomogeneity is found to play a decisive role in determining the decoupling. The difference in the viscosity dependence between rotation and translational diffusion coefficient is discussed.

Robust solver based on modified particle swarm optimization for improved solution of diffusion transport through containment facilities

Accurate estimation of mass transport parameters is necessary for overall design and evaluation processes of the waste disposal facilities. The mass transport parameters, such as effective diffusion coefficient, retardation factor and diffusion accessible porosity, are estimated from observed diffusion data by inverse analysis. Recently, particle swarm optimization (PSO) algorithm has been used to develop inverse model for estimating these parameters that alleviated existing limitations in the inverse analysis. However, PSO solver yields different solutions in successive runs because of the stochastic nature of the algorithm and also because of the presence of multiple optimum solutions. Thus the estimated mean solution from independent runs is significantly different from the best solution. In this paper, two variants of the PSO algorithms are proposed to improve the performance of the inverse analysis. The proposed algorithms use perturbation equation for the gbest particle to gain information around gbest region on the search space and catfish particles in alternative iterations to improve exploration capabilities. Performance comparison of developed solvers on synthetic test data for two different diffusion problems reveals that one of the proposed solvers, CPPSO, significantly improves overall performance with improved best, worst and mean fitness values. The developed solver is further used to estimate transport parameters from 12 sets of experimentally observed diffusion data obtained from three diffusion problems and compared with published values from the literature. The proposed solver is quick, simple and robust on different diffusion problems. (C) 2012 Elsevier Ltd. All rights reserved.

Reactive diffusion in the Ti-Si system and the significance of the parabolic growth constant

Solid diffusion couple experiments are conducted to analyse the growth mechanism of the phases and the diffusion mechanism of the components in the Ti-Si system. The calculation of the parabolic growth constants and the integrated diffusion coefficients substantiates that the analysis is intrinsically prone to erroneous conclusions if it is based on just the parabolic growth constants determined for a multiphase interdiffusion zone. The location of the marker plane is detected based on the uniform grain morphology in the TiSi2 phase, which indicates that this phase grows mainly because of Si diffusion. The growth mechanism of the phases and morphological evolution in the interdiffusion zone are explained with the help of imaginary diffusion couples. The activation enthalpies for the integrated diffusion coefficient of TiSi2 and the Si tracer diffusion are calculated as 190 +/- 9 and 197 +/- 8 kJ/mol, respectively. The crystal structure, details on the nearest neighbours of the components, and their relative mobilities indicate that the vacancies are mainly present on the Si sublattice.

Diffusion on a rugged energy landscape with spatial correlations

Rugged energy landscapes find wide applications in diverse fields ranging from astrophysics to protein folding. We study the dependence of diffusion coefficient (D) of a Brownian particle on the distribution width (epsilon) of randomness in a Gaussian random landscape by simulations and theoretical analysis. We first show that the elegant expression of Zwanzig Proc. Natl. Acad. Sci. U.S.A. 85, 2029 (1988)] for D(epsilon) can be reproduced exactly by using the Rosenfeld diffusion-entropy scaling relation. Our simulations show that Zwanzig's expression overestimates D in an uncorrelated Gaussian random lattice - differing by almost an order of magnitude at moderately high ruggedness. The disparity originates from the presence of ``three-site traps'' (TST) on the landscape - which are formed by the presence of deep minima flanked by high barriers on either side. Using mean first passage time formalism, we derive a general expression for the effective diffusion coefficient in the presence of TST, that quantitatively reproduces the simulation results and which reduces to Zwanzig's form only in the limit of infinite spatial correlation. We construct a continuous Gaussian field with inherent correlation to establish the effect of spatial correlation on random walk. The presence of TSTs at large ruggedness (epsilon >> k(B)T) gives rise to an apparent breakdown of ergodicity of the type often encountered in glassy liquids. (C) 2014 AIP Publishing LLC.

Mode coupling theory analysis of electrolyte solutions: Time dependent diffusion, intermediate scattering function, and ion solvation dynamics

A self-consistent mode coupling theory (MCT) with microscopic inputs of equilibrium pair correlation functions is developed to analyze electrolyte dynamics. We apply the theory to calculate concentration dependence of (i) time dependent ion diffusion, (ii) intermediate scattering function of the constituent ions, and (iii) ion solvation dynamics in electrolyte solution. Brownian dynamics with implicit water molecules and molecular dynamics method with explicit water are used to check the theoretical predictions. The time dependence of ionic self-diffusion coefficient and the corresponding intermediate scattering function evaluated from our MCT approach show quantitative agreement with early experimental and present Brownian dynamic simulation results. With increasing concentration, the dispersion of electrolyte friction is found to occur at increasingly higher frequency, due to the faster relaxation of the ion atmosphere. The wave number dependence of intermediate scattering function, F(k, t), exhibits markedly different relaxation dynamics at different length scales. At small wave numbers, we find the emergence of a step-like relaxation, indicating the presence of both fast and slow time scales in the system. Such behavior allows an intriguing analogy with temperature dependent relaxation dynamics of supercooled liquids. We find that solvation dynamics of a tagged ion exhibits a power law decay at long times-the decay can also be fitted to a stretched exponential form. The emergence of the power law in solvation dynamics has been tested by carrying out long Brownian dynamics simulations with varying ionic concentrations. The solvation time correlation and ion-ion intermediate scattering function indeed exhibit highly interesting, non-trivial dynamical behavior at intermediate to longer times that require further experimental and theoretical studies. (c) 2015 AIP Publishing LLC.

Role of the Molar Volume on Estimated Diffusion Coefficients

The role of the molar volume on the estimated diffusion parameters has been speculated for decades. The Matano-Boltzmann method was the first to be developed for the estimation of the variation of the interdiffusion coefficients with composition. However, this could be used only when the molar volume varies ideally or remains constant. Although there are no such systems, this method is still being used to consider the ideal variation. More efficient methods were developed by Sauer-Freise, Den Broeder, and Wagner to tackle this problem. However, there is a lack of research indicating the most efficient method. We have shown that Wagner's method is the most suitable one when the molar volume deviates from the ideal value. Similarly, there are two methods for the estimation of the ratio of intrinsic diffusion coefficients at the Kirkendall marker plane proposed by Heumann and van Loo. The Heumann method, like the Matano-Boltzmann method, is suitable to use only when the molar volume varies more or less ideally or remains constant. In most of the real systems, where molar volume deviates from the ideality, it is safe to use the van Loo method. We have shown that the Heumann method introduces large errors even for a very small deviation of the molar volume from the ideal value. On the other hand, the van Loo method is relatively less sensitive to it. Overall, the estimation of the intrinsic diffusion coefficient is more sensitive than the interdiffusion coefficient.

Experimental Analysis of Water Vapor Diffusion Through Porous Membranes in a Proton Exchange Membrane Fuel Cell

We report the diffusion characteristics of water vapor through two different porous media, viz., membrane electrode assembly (MEA) and gas diffusion layer (GDL) in a nonoperational fuel cell. Tunable diode laser absorption spectroscopy (TDLAS) was employed for measuring water vapor concentration in the test channel. Effects of the membrane pore size and the inlet humidity on the water vapor transport are quantified through mass flux and diffusion coefficient. Water vapor transport rate is found to be higher for GDL than for MEA. The flexibility and wide range of application of TDLAS in a fuel cell setup is demonstrated through experiments with a stagnant flow field on the dry side.

Relationship between entropy and diffusion: A statistical mechanical derivation of Rosenfeld expression for a rugged energy landscape

Diffusion-a measure of dynamics, and entropy-a measure of disorder in the system are found to be intimately correlated in many systems, and the correlation is often strongly non-linear. We explore the origin of this complex dependence by studying diffusion of a point Brownian particle on a model potential energy surface characterized by ruggedness. If we assume that the ruggedness has a Gaussian distribution, then for this model, one can obtain the excess entropy exactly for any dimension. By using the expression for the mean first passage time, we present a statistical mechanical derivation of the well-known and well-tested scaling relation proposed by Rosenfeld between diffusion and excess entropy. In anticipation that Rosenfeld diffusion-entropy scaling (RDES) relation may continue to be valid in higher dimensions (where the mean first passage time approach is not available), we carry out an effective medium approximation (EMA) based analysis of the effective transition rate and hence of the effective diffusion coefficient. We show that the EMA expression can be used to derive the RDES scaling relation for any dimension higher than unity. However, RDES is shown to break down in the presence of spatial correlation among the energy landscape values. (C) 2015 AIP Publishing LLC.

Observation Of Diffusion Process Of A Water Droplet Immersed In Protein Solution In Microgravity

Pure liquid - liquid diffusion driven by concentration gradients is hard to study in a normal gravity environment since convection and sedimentation also contribute to the mass transfer process. We employ a Mach - Zehnder interferometer to monitor the mass transfer process of a water droplet in EAFP protein solution under microgravity condition provided by the Satellite Shi Jian No 8. A series of the evolution charts of mass distribution during the diffusion process of the liquid droplet are presented and the relevant diffusion coefficient is determined.

Creep and oxygen diffusion in magnetite

The creep rate of polycrystalline Fe3O4 has been measured as a fonction of stress and oxygen partial pressure in the temperature range 480-1100°C. A regime of power law creep is found at high stress, with a stress exponent of ≈- 3.1 and an activation energy of 264 kJ/mol. A regime of diffusional flow is found at lower stresses and is interpreted as Nabarro-Herring creep. The data for the two regimes are combined to deduce an oxygen diffusion coefficient of ≈-10-5 exp(-264 kJ/mol/RT) m2s-1, with oxygen vacancies suggested as the mobile species. © 1990.

Effects of convection and solid wall on the diffusion in microscale convection flows

The diffusive transport properties in microscale convection flows are studied by using the direct simulation Monte Carlo method. The effective diffusion coefficient D is computed from the mean square displacements of simulated molecules based on the Einstein diffusion equation D = x2 t /2t. Two typical convection flows, namely, thermal creep convection and Rayleigh– Bénard convection, are investigated. The thermal creep convection in our simulation is in the noncontinuum regime, with the characteristic scale of the vortex varying from 1 to 100 molecular mean free paths. The diffusion is shown to be enhanced only when the vortex scale exceeds a certain critical value, while the diffusion is reduced when the vortex scale is less than the critical value. The reason for phenomenon of diffusion reduction in the noncontinuum regime is that the reduction effect due to solid wall is dominant while the enhancement effect due to convection is negligible. A molecule will lose its memory of macroscopic velocity when it collides with the walls, and thus molecules are hard to diffuse away if they are confined between very close walls. The Rayleigh– Bénard convection in our simulation is in the continuum regime, with the characteristic length of 1000 molecular mean free paths. Under such condition, the effect of solid wall on diffusion is negligible. The diffusion enhancement due to convection is shown to scale as the square root of the Péclet number in the steady convection regime, which is in agreement with previous theoretical and experimental results. In the oscillation convection regime, the diffusion is more strongly enhanced because the molecules can easily advect from one roll to its neighbor due to an oscillation mechanism. © 2010 American Institute of Physics. doi:10.1063/1.3528310

Configuration-dependent diffusion can shift the kinetic transition state and barrier height of protein folding

We show that diffusion can play an important role in protein-folding kinetics. We explicitly calculate the diffusion coefficient of protein folding in a lattice model. We found that diffusion typically is configuration- or reaction coordinate-dependent. The diffusion coefficient is found to be decreasing with respect to the progression of folding toward the native state, which is caused by the collapse to a compact state constraining the configurational space for exploration. The configuration- or position-dependent diffusion coefficient has a significant contribution to the kinetics in addition to the thermodynamic free-energy barrier. It effectively changes (increases in this case) the kinetic barrier height as well as the position of the corresponding transition state and therefore modifies the folding kinetic rates as well as the kinetic routes. The resulting folding time, by considering both kinetic diffusion and the thermodynamic folding free-energy profile, thus is slower than the estimation from the thermodynamic free-energy barrier with constant diffusion but is consistent with the results from kinetic simulations. The configuration- or coordinate-dependent diffusion is especially important with respect to fast folding, when there is a small or no free-energy barrier and kinetics is controlled by diffusion.Including the configurational dependence will challenge the transition state theory of protein folding.