微孔中简单流体扩散行为的分子动力学模拟研究

用分子动力学模拟方法研究了受限在微孔中的简单流体氩的扩散行为,考察了微孔类型、孔径、温度和密度对微孔中流体扩散系数的影响.研究发现,微孔中流体的扩散系数均小于体相流体,并且随孔径的减小而减小,同时沿孔道或狭缝方向的扩散系数分量远大于沿孔径方向的分量,并且流体在通道型微孔中的扩散系数小于在狭缝型微孔中的扩散系数,温度和密度也是影响微孔中扩散的重要因素.     

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.     

疏水性微孔中水的结构和扩散性质的分子模拟

用分子动力学（MD）方法模拟了受限在疏水性微孔中的水的结构与动力学行为.分别考察了孔径、温度和压力对水在孔道方向的密度分布和自扩散系数的影响,计算了不同温度下水的径向分布函数.发现在小孔径的微孔中,随着温度的降低,水分子沿孔道的分布逐渐变得不均匀,最终导致气-液相分离,微孔孔道内有明显的分段现象.受限在小孔径微孔中水的自扩散系数大约为体相流体水的20％～30％,并且随着孔径的减小,自扩散系数也减小.同时还发现沿孔道方向的自扩散系数分量大约为孔径方向的4～5倍.提出了微孔中水自扩散系数的关联模型.     

Transport properties and distribution of water molecules confined in hydrophobic nanopores and nanoslits

The transport properties, including the diffusivity and viscosity, of water confined in hydrophobic nanopores and nanoslits were studied by molecular dynamics simulations. The results show that the diffusion coefficient in nanopores and nanoslits is markedly lower than that in the bulk. But the viscosity is much larger than that in bulk. The parallel diffusion coefficient is obviously larger than the perpendicular ones. The diffusion coefficient in the channel pore is ever less than that in the slit pore at the same pore width, but the viscosity is larger. The temperature and density affect significantly the diffusivity and viscosity in nanopores and nanoslits. Lower density water exhibits some special characteristics on density profiles in nanopores and nanoslits at lower temperatures, and the density profiles show a change from homogeneous to inhomogeneous as the pore width is reduced. Even clusters occurred in micropores.     

Transport Properties of Fluids in Micropores by Molecular Dynamics Simulation

The transport properties of fluid argon in micropores, i.e. diffusivity and viscosity, were studied by molecular dynamics simulations. The effects of pore width, temperature and density on diffusivity and viscosity were analyzed in micropores with pore widths from 0.8 to 4.0 nm. The results show that the diffusivity in micropores is much lower than the bulk diffusivity, and it decreases as the pore width decreases; but the viscosity in micropores is significantly larger than the bulk one, and it increases sharply in narrow micropores. The diffusivity in channel parallel direction is obviously larger than that in channel perpendicular direction. The temperature and density are important factors that obviously affect diffusivity and viscosity in micropores.     

Molecular Simulation on Transport Properties of Confined Water

The diffusivity and viscosity of water confined in micropores were studied by molecular dynamics simulations. The effects of pore width and density were analyzed at pore widths from 0.9 to 2.6nm. The diffusivity in micropores is lower than that of the bulk, and it decreases as pore width decreases and as density increases. But the viscosity in micropores is much larger than that of the bulk, and it increases as pore width decreases and as density increases. The diffusivity in channel parallel direction is obviously larger than that in channel perpendicular directions.     

Grand canonical-like molecular dynamics simulations: application to anisotropic mass diffusion in a nanoporous medium

In this work, we describe two grand canonical-like molecular dynamics approaches to investigate mass diffusion phenomenon of a simple Lennard-Jones fluid confined between solid surfaces and in direct contact with reservoirs. In the first method, the density is used as the control variable in the reservoir whereas it is the pressure in the second method. Both methods provide consistent results, however, the constant density approach is the most efficient with respect to the computational time and implementation. Then, employing the constant density approach, we have studied the transient behavior of the diffusion process associated with the migration of one fluid into another one confined between parallel solid walls. Results have shown that the evolution of molar fraction of the invading fluid follows roughly a 1D diffusion model when the solid phase is weakly or moderately adsorbent with a characteristic time increasing when the pore width decreases. However, when the adsorption is high and the pore width small (i.e., below ten molecular sizes), the apparent mass diffusion in the adsorbed layer is reduced compared to that in the center of the slit pore. Hence, this mass diffusion process becomes a two-dimension phenomenon that must take into account an effective mass diffusion coefficient varying locally.     

Structural and dynamical properties for confined polymers undergoing planar Poiseuille flow

The authors present the results from nonequilibrium molecular dynamics simulations for the structural and dynamical properties of highly confined linear polymer fluids undergoing planar Poiseuille flow. They study systems confined within pores of several atomic diameters in width and investigate the dependence of the density profiles, the mean squared radius of gyration, the mean squared end-to-end distance, streaming velocity, strain rate, shear stress, and streaming angular velocity as functions of average fluid density and chain length. Their simulation results show that, sufficiently far from the walls, the radius of gyration for molecules under shear in the middle of the pore follows the power law Rg=ANbnu, where Nb is the number of bonds and the exponent has a value of 0.5 which resembles the value for a homogeneous equilibrium fluid. Under the conditions simulated, the authors find the onset of flat velocity profiles but with very little wall slippage. These flat profiles are most likely due to the restricted layering of the fluid into just one or two molecular layers for narrow pore widths compared to chain length, rather than typical plug-flow conditions. The angular velocity is shown to be proportional to half the strain rate in the pore interior when the chain length is sufficiently small compared to the pore width, consistent with the behavior for homogeneous fluids in the linear regime.     

Relaxation and short time dynamics of bulk liquids and fluids confined in spherical cavities and slit pores

The density of states for bulk and confined fluids have been modeled using a recently proposed gamma distribution (Krishnan, S. H.; Ayappa, K. G. J. Chem. Phys. 2004, 121, 3197). The gamma distribution results in a closed form analytical expression for the velocity autocorrelation function and the relaxation time of the fluid. The two parameters of the gamma distribution are related analytically to the second and fourth frequency moments of the fluid using short time expansions. The predictions by the proposed gamma model are compared with the velocity autocorrelation functions obtained using the theory of instantaneous normal modes (INMs) and from molecular dynamics simulations. The model is applied to a bulk soft sphere liquid and fluids confined in a spherical cavity and slit-shaped pores. The gamma model is able to capture the resulting changes in relaxation time due to changes in density and temperature extremely well for both the bulk liquid and confined inhomogeneous fluid situations. In all cases, the predictions by the gamma model are superior to those obtained from the INM theory. In the case of the fluid confined in a slit pore, the loadings were obtained from a grand canonical Monte Carlo simulation where the pore is equilibrated with a bulk fluid. This is similar to a confinement situation in a surface force apparatus. The predicted relaxation times vs pore widths from the gamma model are seen to accurately capture the oscillations due to formation and disruption of layers within the slit pore.     

The effect of corrugation of pore wall on the thermal diffusion in nanopores by molecular dynamics simulations

In this work we have studied the effect of corrugation on the thermal diffusion (soret effect) in isotopic and non-isotopic fluid mixtures confined in a slit pore. We used a boundary driven non-equilibrium molecular dynamics to simulate thermal diffusion in Lennard–Jones (LJ) binary mixtures confined in structureless Steele 10-4-3 and atomistic Lennard–Jones pore walls. The results showed that for the isotopic mixture thermal diffusion factor for both wall types agrees and the corrugation of the LJ wall has no effect in isotopic mixture. However, for non-isotopic mixture confined in atomistic LJ pore the component with stronger attraction adsorbs more to the wall than the structureless Steele wall. The effect of corrugation of pore wall on the thermal diffusion is noticeable in narrow slit pore and mixture with large difference in molecular attraction parameter of components.     

Mean field kinetic theory for a lattice gas model of fluids confined in porous materials

We consider the mean field kinetic equations describing the relaxation dynamics of a lattice model of a fluid confined in a porous material. The dynamical theory embodied in these equations can be viewed as a mean field approximation to a Kawasaki dynamics Monte Carlo simulation of the system, as a theory of diffusion, or as a dynamical density functional theory. The solutions of the kinetic equations for long times coincide with the solutions of the static mean field equations for the inhomogeneous lattice gas. The approach is applied to a lattice gas model of a fluid confined in a finite length slit pore open at both ends and is in contact with the bulk fluid at a temperature where capillary condensation and hysteresis occur. The states emerging dynamically during irreversible changes in the chemical potential are compared with those obtained from the static mean field equations for states associated with a quasistatic progression up and down the adsorption/desorption isotherm. In the capillary transition region, the dynamics involves the appearance of undulates (adsorption) and liquid bridges (adsorption and desorption) which are unstable in the static mean field theory in the grand ensemble for the open pore but which are stable in the static mean field theory in the canonical ensemble for an infinite pore.     

Effect of solvation and confinement on the trans-gauche isomerization reaction in n-butane

The effect of solvation and confinement on the conformational equilibria and kinetics of n-butane is examined using molecular dynamics simulations of the bulk and confined fluids and compared to appropriately chosen reference states. Clear evidence for a solvent shift of the preferred conformation in bulk n-butane is found. At a temperature of 292 K and a density of 6.05 nm-3 a small solvent shift in favor of gauche is observed (similar to previously reported values), and the shift increases substantially with an increase in density to 8.28 nm-3. The rate of torsional interconversion from the trans to the gauche state, calculated using the relaxation function method, was found to increase with increasing temperature and density. The rate constants kTG and kGT have an Arrhenius temperature dependence yielding activation energies significantly lower than the trans-gauche and gauche-trans barrier heights in the torsional potential for a free molecule, depending on the density. In the confined phase, we considered the same densities as simulated in the bulk phase, and for four different values of the physical pore width (approximately 1.5-4.0 nm). At the high density, we find that the position of the trans-gauche equilibrium is displaced towards excess trans compared with the bulk phase, reflecting the confinement and interactions of the molecules with the pore wall. The isomerization rate is found to decrease with decreasing pore width. Again, we find that the kinetics obeys an Arrhenius rate law and the activation energy for the trans-gauche and gauche-trans interconversions is slightly smaller than that of the bulk fluid at the same density.     

Diffusion enhancement in core-softened fluid confined in nanotubes

We study the effect of confinement in the dynamical behavior of a core-softened fluid. The fluid is modeled as a two length scales potential. This potential in the bulk reproduces the anomalous behavior observed in the density and in the diffusion of liquid water. A series of NpT molecular dynamics simulations for this two length scales fluid confined in a nanotube were performed. We obtain that the diffusion coefficient increases with the increase of the nanotube radius for wide channels as expected for normal fluids. However, for narrow channels, the confinement shows an enhancement in the diffusion coefficient when the nanotube radius decreases. This behavior, observed for water, is explained in the framework of the two length scales potential.     

Phase Equilibria of Hydrogen Bonding Fluid in a Slit Pore with Broken Symmetry

Phase equilibria of hydrogen bonding (HB) fluid confined in a slit pore with broken symmetry were investigated by the density functional theory incorporated with modified fundamental measure theory, where the symmetry breaking originated from the distinct interactions between fluid molecules and two walls of the slit pore. In terms of adsorption-desorption isotherms and the corresponding grand potentials, phase diagrams of HB fluid under various conditions are presented. Furthermore, through phase coexistences of laying transition and capillary condensation, the effects of HB interaction, pore width, fluid-pore interaction and the broken symmetry on the phase equilibrium properties are addressed. It is shown that these factors can give rise to apparent influences on the phase equilibria of confined HB fluid because of the competition between intermolecular interaction and fluid-pore interaction. Interestingly, a significant influence of broken symmetry of the slit pore is found, and thus the symmetry breaking can provide a new way to regulate the phase behavior of various confined fluids.     

Anisotropic dynamics of dipolar liquids in narrow slit pores

We report molecular dynamics simulation results for Stockmayer fluids confined to narrow slitlike pores with structureless, nonconducting walls. The translational and rotational dynamics of the dipolar particles have been investigated by calculating autocorrelation functions, diffusion coefficients, and relaxation times for various pore widths (five or less particle diameters) and directions parallel and perpendicular to the walls. The dynamic properties of the confined systems are compared to bulk properties, where corresponding bulk and pore states at the same temperature and chemical potential are determined in parallel grand canonical Monte Carlo simulations. We find that the dynamic behavior inside the pore depends on the distance from the walls and can be strongly anisotropic even in globally isotropic systems. This concerns especially the particles in the surface layers close to the walls, where the single particle and collective dipolar relaxation resemble that of true two-dimensional dipolar fluids with different in-plane and out-of-plane relaxations. On the other hand, bulklike relaxation is observed in the pore center of sufficiently wide pores.     

微孔中简单流体粘度的分子动力学模拟及关联模型

用分子动力学模拟计算了微孔介质中流体氩在不同温度、不同密度和不同孔径下的剪切粘度.并根据Chapman-Enskog关于硬球流体传递性质的理论以及Heyes的关于Lennard-Jones流体粘度的表达式,提出了两个描述微孔介质中流体粘度的模型,该模型可以计算微孔中流体氩在不同状态下的粘度值.通过与计算机模拟值的比较,证明这两个微孔流体粘度模型是可用的.     

Tractable molecular theory of transport of Lennard-Jones fluids in nanopores

We present here a tractable theory of transport of simple fluids in cylindrical nanopores, which is applicable over a wide range of densities and pore sizes. In the Henry law low-density region the theory considers the trajectories of molecules oscillating between diffuse wall collisions, while at higher densities beyond this region the contribution from viscous flow becomes significant and is included through our recent approach utilizing a local average density model. The model is validated by means of equilibrium as well nonequilibrium molecular dynamics simulations of supercritical methane transport in cylindrical silica pores over a wide range of temperature, density, and pore size. The model for the Henry law region is exact and found to yield an excellent match with simulations at all conditions, including the single-file region of very small pore size where it is shown to provide the density-independent collective transport coefficient. It is also shown that in the absence of dispersive interactions the model reduces to the classical Knudsen result, but in the presence of such interactions the latter model drastically overpredicts the transport coefficient. For larger micropores beyond the single-file region the transport coefficient is reduced at high density because of intermolecular interactions and hindrance to particle crossings leading to a large decrease in surface slip that is not well represented by the model. However, for mesopores the transport coefficient increases monotonically with density, over the range studied, and is very well predicted by the theory, though at very high density the contribution from surface slip is slightly overpredicted. It is also seen that the concept of activated diffusion, commonly associated with diffusion in small pores, is fundamentally invalid for smooth pores, and the apparent activation energy is not simply related to the minimum pore potential or the adsorption energy as generally assumed.     

Capillary condensation of short-chain molecules

A density-functional study of capillary condensation of fluids of short-chain molecules confined to slitlike pores is presented. The molecules are modeled as freely jointed tangent spherical segments with a hard core and with short-range attractive interaction between all the segments. We investigate how the critical parameters of capillary condensation of the fluid change when the pore width decreases and eventually becomes smaller than the nominal linear dimension of the single-chain molecule. We find that the dependence of critical parameters for a fluid of dimers and of tetramers on pore width is similar to that of the monomer fluid. On the other hand, for a fluid of chains consisting of a larger number of segments we observe an inversion effect. Namely, the critical temperature of capillary condensation decreases with increasing pore width for a certain interval of values of the pore width. This anomalous behavior is also influenced by the interaction between molecules and pore walls. We attribute this behavior to the effect of conformational changes of molecules upon confinement.     

Influence of a Surface Microheterogeneity on the Disjoining Pressure of Adsorbed Lennard–Jones Fluid. Computer Simulation

Computer simulation of the adsorption of Lennard–Jones fluid in slitlike micropores with structured walls was performed by the Monte Carlo and molecular dynamic methods. The influence of pore width, temperature, pressure, and surface structure on the disjoining pressure was considered. Structural and energy characteristics of the surface affect significantly the value of the disjoining pressure. Among the surface structural characteristics, the important role is played by the number and localization of vacancy-type defects.     

Nanostructures in a binary mixture confined in slit-like pores with walls decorated with tethered polymer brushes in the form of stripes: dissipative particle dynamics study

Using dissipative particle dynamics, we investigate the behavior of a binary mixture, exhibiting demixing in a bulk phase, confined in slit-like pores with walls modified by the stripes of tethered brush of chains. Our main interest is to determine possible morphologies that can be formed inside the pore, depending on the geometrical parameters characterizing the system (the size of the pore and the width of the stripes). In order to describe the observed morphologies we calculate several characteristics, as the density and local temperature profiles, the radii of gyration for the attached polymers, and the minimum polymer-polymer distances in the direction parallel and perpendicular to the pore walls. The summary of our findings is presented as a sketch of the diagram of morphologies.