Kinetic theory of gas separation in a nanopore and comparison to molecular dynamics simulation

Kinetic mesoscopic theory derived from an atomistic model is applied to study permeation and separation of gases in a single rectangular pore. The goal is to judge the analytical method against the results of molecular dynamics simulation and to demonstrate the ease and relevance of analytical theories to calculate density profiles, flux, permeance, and separation factors. The permeance is linked to the amount of gas adsorbed in the pore and the effect of the effective gas-wall interaction on adsorption is explored. The effects of pore size, temperature, and the parameters of the pore wall interaction are investigated and reproduce the trends found in the numerical simulation of permeation of a mixture of methane and carbon dioxide in a carbon nanopore.     

Density inhomogeneity and diffusion behavior of fluids in micropores by molecular-dynamics simulation

The density profiles and the diffusion behavior of fluid argon confined in micropores were studied by molecular-dynamics simulations. The effects of pore size (width), temperature and number density on the density profiles and the self-diffusion coefficients in micropores were simulated with pore widths from 0.6 to 4.0 nm. The density profiles are greatly affected by the pore size. Strong inhomogeneities in the channel direction and vapor-liquid phase separation in the micropores were observed when initial conditions were chosen in the coexistence region of the fluid. The self-diffusion coefficient in the channel direction in the pores was found to be much lower than in the bulk, and decreasing with decreasing pore size, decreasing temperature, and increasing density.     

金属-有机骨架材料中吸附气体的扩散速率

采用分子动力学方法,以甲烷为探针分子研究了不同压力条件下气体在具有不同孔道结构的金属-有机骨架材料(MOFs)中的扩散速率.通过计算气体在八种材料中的自扩散系数,并结合气体分子在材料中的质心分布图等,讨论了气体扩散速率与孔道结构之间的关系.研究结果表明:对于同时含有孔笼(pocket)和三维正交孔道(channel)结构的MOF材料(P-C材料),低压时甲烷气体吸附在孔笼结构中,随着压力的升高,气体分子开始进入正交孔道,同时其自扩散系数增加;而对于只含有三维立方孔道结构的IRMOF(isoreticular MOF)系列材料,在中低压范围内,气体分子在其中的自扩散系数随压力变化较小.当压力进一步升高时,气体分子在材料孔道中的吸附逐渐接近饱和,其自扩散系数均降低.因此,在不同MOF材料中气体分子扩散速率的差异主要取决于孔道结构的不同.对P-C材料,中低压下通过控制压力可以控制气体在其中的扩散速率,从而为MOF材料在气体存储、分离等方面的实际应用提供参考信息.     

Surface self-diffusion of organic molecules adsorbed in porous silicon

The pulsed field gradient nuclear magnetic resonance method has been employed to probe self-diffusion of organic guest molecules adsorbed in porous silicon with a 3.6 nm pore size. The molecular self-diffusion coefficient and intrapore adsorption were simultaneously measured as a function of the external vapor pressure. The latter was varied in a broad range to provide pore loading from less than monolayer surface coverage to full pore saturation. The measured diffusivities are found to be well-correlated with the adsorption isotherms. At low molecular concentrations in the pores, corresponding to surface coverages of less than one monolayer, the self-diffusion coefficient strongly increases with increasing concentration. This observation is attributed to the occurrence of activated diffusion on a heterogeneous surface. Additional experiments in a broad temperature range and using binary mixtures confirm this hypothesis.     

Structure and dynamics of micelles and cubic phase structures with ethoxylated phytosterol surfactant in water

The self-assembly of a sterol ethoxylate surfactant with 30 oxyethylene units in water was studied by 1H NMR self-diffusion measurements in a wide concentration range in the micellar region (0-25 wt %). The data showed that the surfactant aggregates do not interact by hard sphere interactions but rather a strong concentration dependence of the diffusion coefficient was noted which was explained by polymer scaling theory. In the cubic phase (30-65 wt %), the self-diffusion data from water, from surfactant, and from free polyoxyethylene suggest spherical micelles, although water diffusion was much restricted due to binding to the surfactant headgroup. From X-ray measurements in the cubic phase, the unit cell size was calculated, and together with surfactant self-diffusion measurements the exchange dynamics between free and aggregated surfactant was obtained.     

Structure-based prediction of the conductance properties of ion channels

The HOLE procedure allows the prediction of the absolute conductance of an ion channel model from its structure. The original prediction method uses an empirically corrected Ohmic method. It is most successful, with predictions being reliable to within a factor of two. A new modification of the procedure is presented in which the self-diffusion coefficients of water molecules from molecular dynamics simulation are used to replace the empirical correction factor. A "prediction" of the conductance for the porin OmpF by the new method is made and shown to be very close to the experimental value. HOLE also allows the prediction of the effect that the addition of non-electrolyte polymers will have on channel conductance. The method has great potential to yield structural information from data provided by single channel recordings but needs further validation by making measurements on channels of known structure. Preliminary results are given of single channel records establishing the effects of non-electrolytes on the conductance of gramicidin D channels. As an example of the potential uses of the procedure application is made to examine the oligomerization of alpha-toxin (alpha-hemolysin) channels. A model for the alpha-toxin hexamer, based on the crystal structure for the heptamer, is generated using molecular mechanics methods. The compatibility of the structures with single channel conductance data is assessed using HOLE.     

Formation of columns of anomalous water by condensation in quartz capillaries

It is shown that these columns are formed in the kinetic mode from unsaturated water vapor. The condensation coefficient for an anomalous column is much less than for ordinary water. The growth rate is dependent on the equilibration time for the concentration of the anomalous component in the column. The diffusion coefficient for the anomalous component in water is less than the self-diffusion coefficient of water by nearly an order of magnitude; this and the low volatility are to be ascribed to large molecular size. It is supposed that the molecules of the anomalous component are strong groups of H₂O molecules that are only slowly converted to monomers.     

Transport coefficients of the Lennard-Jones model fluid. II Self-diffusion

In an extensive computer simulation study, the transport coefficients of the Lennard-Jones model fluid were determined with high accuracy from equilibrium molecular-dynamics simulations. In the frame of time-correlation function theory, the generalized Einstein relations were employed to evaluate the transport coefficients. This second of a series of four papers presents the results for the self-diffusion coefficient, and discusses and interprets the behavior of this transport coefficient in the fluid region of the phase diagram. The uncertainty of the self-diffusion data is estimated to be 1% in the gas region and 0.5% at high-density liquid states. With the very accurate data, even fine details in the shape of the self-diffusion isotherms are resolved, and the previously little-investigated behavior of the self-diffusion coefficient at low-density gaseous states is analyzed in detail. Finally, aspects of the mass transport mechanisms on the molecular scale are explored by an analysis of the velocity autocorrelation functions.     

Study of Self-Diffusion ofn-Decane in NaX Zeolite by the Pulsed Magnetic-Field Gradient NMR Method

Translational mobility of n-decane molecules in a porous space of NaX zeolite was studied within the wide ranges of diffusion times and temperatures. The dependence of the effective self-diffusion coefficient on diffusion time was established. Confined mobility of diffusant molecules inside the crystallite was observed both for complete and partial filling of NaX pores with a liquid, when the adsorption barrier was absent at the interface between intra- and intercrystallite regions. It was suggested that obstacles are present at the surface of NaX crystallites complicating the transfer of liquid molecules from crystallite channels to intercrystallite space. True value of self-diffusion coefficient ofn-decane in the itracrystallite space of NaX was determined and its dependence on the concentration of liquid molecules in zeolite channels was considered. A special attention was paid to the study of molecular exchange between intra- and intercrystallite-confined liquids.     

Quantum diffusion in liquid water from ring polymer molecular dynamics

We have used the ring polymer molecular-dynamics method to study the translational and orientational motions in an extended simple point charge model of liquid water under ambient conditions. We find, in agreement with previous studies, that quantum-mechanical effects increase the self-diffusion coefficient D and decrease the relaxation times around the principal axes of the water molecule by a factor of around 1.5. These results are consistent with a simple Stokes-Einstein picture of the molecular motion and suggest that the main effect of the quantum fluctuations is to decrease the viscosity of the liquid by about a third. We then go on to consider the system-size scaling of the calculated self-diffusion coefficient and show that an appropriate extrapolation to the limit of infinite system size increases D by a further factor of around 1.3 over the value obtained from a simulation of a system containing 216 water molecules. These findings are discussed in light of the widespread use of classical molecular-dynamics simulations of this sort of size to model the dynamics of aqueous systems.     

PGSE-NMR studies of solvent diffusion in poly(N-isopropylacrylamide) colloidal microgels

The solvent self-diffusion coefficient has been studied in thermoshrinking poly(N-isopropyl acrylamide) microgel dispersions by the pulsed-gradient spin-echo PGSE-NMR technique, as a function of temperature and mass fraction. After suitable corrections for the temperature, the H₂O/D₂O ratio and the relative volume fractions, all the self-diffusion data obtained over a temperature range of approximately 40 °C and mass fraction (2–12 % wt/wt) could be superimposed with the volume fraction as the universal factor. The observed reduction in the solvent self-diffusion coefficient with volume fraction was greater than that predicted by simple obstruction theory. After correction for-, and the subsequent removal of the obstruction effect, the diffusion of the solvent through the core of the particle is elucidated. As found for other polymer-solvent systems, there were no specific binding effects. The diffusion of the solvent in these dispersions over such temperature and mass fraction ranges could be rationalised assuming a constant solvent self-diffusion coefficient in the core of the particles.     

Concentration dependences of heat conductivity coefficients of inert gases in narrow pores

Calculation of the transfer of molecules in porous systems requires self-consistent expressions describing the kinetic transfer coefficients for various concentrations and temperatures. The concentration dependences of heat conductivity and self-diffusion coefficients for fluids with different densities, ranging from rarefied gases to liquids, were considered in terms of a unified model. For monoatomic gases (argon), the model takes into account two energy transfer channels, namely, the vacancy mechanism and energy transfer through collisions of molecules. The former channel is characteristic of rarefied gases, while the latter is noted for condensed phases. The energy parameters of the model were determined on the basis of data on the heat conductivity coefficient in the bulk phase. The heat conductivity coefficient follows a linear temperature dependence for low density; in the medium and large density regions, these dependences follow a more complex pattern that changes depending on temperature. The influence of the interaction of atoms with the pore walls on the concentration dependences of the heat conductivity coefficients was investigated for different total amounts of the adsorbate. These coefficients depend appreciably on the distance to the pore wall and on the direction of heat transfer.     

Structure,thermodynamics and diffusion in asymmetric binary mixtures: a molecular dynamics simulation study

Structural and thermodynamic properties as well as diffusion coefficients of binary fluid mixtures with asymmetry in mass, size, charge and their combinations have been studied using classical molecular dynamics simulations. The fluid mixture is modelled as spherical particles interacting via the Weeks–Chandler–Andersen and Coulomb potential. The diameter, charge and mass of the fluid particles are in the range 6–60 Å, 1–10e and 1—500 amu, respectively. Systematic variations in pair-correlation functions, thermodynamic properties as well as the self-diffusion coefficient are found with the size, charge and mass ratio of the particles. The self-diffusion coefficient for systems having more than one type of asymmetry is calculated and expressed in terms of diffusion coefficients of systems with only one type of asymmetry.     

甘油-水-氯化钠三元溶液中甘油浓度对甘油自扩散系数的影响

利用分子动力学模拟方法研究了浓度对甘油-水-氯化钠三元溶液中甘油自扩散系数的影响. 随着甘油浓度的增大, 甘油的自扩散系数逐渐减小. 氢键分析表明, 甘油自扩散系数的减小来源于其参与的甘油-水氢键数目的减少和甘油-甘油氢键数目的增加.     

Computation of some thermodynamic properties of nitrogen using a new intermolecular potential from molecular dynamics simulation

A new pair-potential energy function of nitrogen has been determined via the inversion of reduced viscosity collision integrals and fitted to obtain an analytical potential form. The pair-potential reproduces the second virial coefficient, viscosity, thermal conductivity, self-diffusion coefficient, and thermal diffusion factor of nitrogen in a good accordance with experimental data over wide ranges of temperatures and densities. We have also performed the molecular dynamics simulation to obtain pressure, internal energy, heat capacity at constant volume, and self-diffusion coefficient of nitrogen at different temperatures and densities using our calculated pair-potential and some other potentials. The molecular dynamics of the nitrogen molecules has been also used to determine nitrogen equation of state in two (low and high) pressure ranges. Our results are in a good agreement with experiment and literature values.     

Water transport coefficient distribution through the membrane in a polymer electrolyte fuel cell

The net water transport coefficient through the membrane, defined as the ratio of the net water flux from the anode to cathode to the protonic flux, is used as a quantitative measure of water management in a polymer electrolyte fuel cell (PEFC). In this paper we report on experimental measurements of the net water transport coefficient distribution for the first time. This is accomplished by making simultaneous current and species distribution measurements along the flow channel of an instrumented PEFC via a multi-channel potentiostat and two micro gas chromatographs. The net water transport coefficient profile along the flow channels is then determined by a control-volume analysis under various anode and cathode inlet relative humidity (RH) at 80 °C and 2 atm. It is found that the local current density is dominated by the membrane hydration and that the gas RH has a large effect on water transport through the membrane. Very small or negative water transport coefficients are obtained, indicating strong water back diffusion through the 30 μm Gore-Select^® membrane used in this study.     

Self-diffusion of methane in single-walled carbon nanotubes at sub- and supercritical conditions

The diffusivities of methane in single-walled carbon nanotubes (SWNTs) are investigated at various temperatures and pressures using classical molecular dynamics (MD) simulations complemented with grand canonical Monte Carlo (GCMC) simulations. The carbon atoms at the nanotubes are structured according to the (m, m) armchair arrangement and the interactions between each methane molecule and all atoms of the confining surface are explicitly considered. It is found that the parallel self-diffusion coefficient of methane in an infinitely long, defect-free SWNT decreases dramatically as the temperature falls, especially at subcritical temperatures and high loading of gas molecules when the adsorbed gas forms a solidlike structure. With the increase in pressure, the diffusion coefficient first declines rapidly and then exhibits a nonmonotonic behavior due to the layering transitions of the adsorbed gas molecules as seen in the equilibrium density profiles. At a subcritical temperature, the diffusion of methane in a fully loaded SWNT follows a solidlike behavior, and the value of the diffusion coefficient varies drastically with the nanotube diameter. At a supercritical temperature, however, the diffusion coefficient at high pressure reaches a plateau, with the limiting value essentially independent of the nanotube size. For SWNTs with the radius larger than approximately 2 nm, capillary condensation occurs when the temperature is sufficiently low, following the layer-by-layer adsorption of gas molecules on the nanotube surface. For SWNTs with a diameter less than about 2 nm, no condensation is observed because the system becomes essentially one-dimensional.     

Fluid transport in nanochannels induced by temperature gradients

We investigate the mechanisms of fluid transport driven by temperature gradients in nanochannels through molecular dynamics simulations. It is found that the fluid-wall interaction is critical in determining the flow direction. In channels of very low surface energy, where the fluid-wall binding energy ε(fw) is small, the fluid moves from high to low temperature and the flow is induced by a potential ratchet near the wall. In high surface energy channels, however, the fluid is pumped from low to high temperature and the pressure drop caused by the temperature gradient is the major driving force. In addition, as the fluid-wall interaction is strengthened, the flow flux assumes a maximum, where ε(fw) is close to the lower temperature T(L) of the channel and ε(fw)/kT(L) ≈ 1 is roughly satisfied.     

Pair correlation functions and the self-diffusion coefficient of Lennard-Jones liquid in the modified free volume theory of diffusion

In this paper, we apply the Matteoli-Mansoori empirical formula for the pair correlation function of simple fluids obeying the Lennard-Jones potential to calculate reduced self-diffusion coefficients on the basis of the modified free volume theory. The self-diffusion coefficient thus computed as functions of temperature and density is compared with the molecular dynamics simulation data and the self-diffusion coefficient obtained by the modified free volume theory implemented with the Monte Carlo simulation method for the pair correlation function. We show that the Matteoli-Mansoori empirical formula yields sufficiently accurate self-diffusion coefficients in the supercritical regime, provided that the minimum free volume activating diffusion is estimated with the classical turning point of binary collision at the mean relative kinetic energy 3k(B)T/2, where k(B) is the Boltzmann constant and T is the temperature. In the subcritical regime, the empirical formula yields qualitatively correct, but lower values for the self-diffusion coefficients compared with computer simulation values and those from the modified free volume theory implemented with the Monte Carlo simulations for the pair correlation function. However, with a slightly modified critical free volume, the results can be made quite acceptable.