Phase coexistence and dynamic properties of water in nanopores

The dynamical properties of a confined fluid depend strongly on the (spatially varying) density. Its knowledge is therefore an important prerequisite for molecular-dynamics (MD) simulations and the analysis of experimental data. In a mixed Gibbs ensemble Monte Carlo (GEMC)/MD simulation approach we first apply the GEMC method to find possible phase states of water in hydrophilic and hydrophobic nanopores. The obtained phase diagrams evidence that a two-phase state is the most probable state of a fluid in incompletely filled pores in a wide range of temperature and level of pore filling. Pronounced variations of the average and local densities are observed. Subsequently, we apply constant-volume MD simulations to obtain water diffusion coefficients and to study their spatial variation along the pore radius. In general, water diffusivity slightly decreases in a hydrophilic pore and noticeably increases in a hydrophobic pore (up to about 40% with respect to the bulk value). In the range of gradual density variations the local diffusivity essentially follows the inverse density and the water binding energy. The diffusivity in the quasi-two-dimensional water layers near the hydrophilic wall decreases by 10 to 20% with respect to the bulk value. The average diffusivity of water in incompletely filled pore is discussed on the basis of the water diffusivities in the coexisting phases.Received: 1 January 2003, Published online: 14 October 2003PACS: 61.20.Ja Computer simulation of liquid structure - 64.70.Fx Liquid-vapor transitions     

Molecular simulation study of triangle-well fluids confined in slit pores

Grand canonical Monte Carlo simulations are used to study the behaviour of triangle-well (TW) fluids with variable well widths confined inside slit pores. The effect of individual factors influencing the properties of confined fluids such as fluid–fluid interactions, pore size and pore wall–fluid interactions are obtained using simulations as it is difficult to experimentally determine the same. An interesting observation of this study is that inside the narrow pore of slit height h* = 5 at the high-pressure condition of P* = 0.8, for the TW fluid with long-range attraction or for the fluid at a low temperature for even a short-range attraction, the density profiles show layering such that there is a sticking tendency of the particles at centre, while there is a depletion of particles near the wall (as the layers at the centre have higher density peak heights than near the walls).     

Behavior of a wetting phase near a solid boundary: vapor near a weakly attractive surface

Density profiles of a LJ vapor near a weakly attractive surface with long-range fluid wall potential was studied along the pore coexistence curve. There are two localized density maxima near the pore wall: the first one is caused by localization of the molecules in the minimum of the fluid-wall potential, and the second one reflects adsorption of molecules at the first layer at higher densities. In addition, a third, weak density maximum is observed close to the critical temperature due to the competition between the long-range attractive tail of the fluid-wall potential and the effect of missing neighbors. This maximum separates the region of a gradual density depletion toward the surface due to the missing neighbor effect and the adsorption region further from the surface, where the density gradually increases toward the surface due to the attractive fluid-wall potential. When approaching the bulk critical temperature, this maximum moves away from the surface due to the divergence of the bulk correlation length. Applicability of various equations to describe the vapor density profiles is examined. Excess adsorption of vapor at low temperatures turns into excess depletion at higher temperatures. The crossover temperature increases with increasing pore size and strengthening fluid-wall interaction. The problems of the theory of the surface critical behavior of Ising models in a case of a non vanishing surface field and its mapping on a fluid is discussed.     

Adsorption isotherm of a Lennard-Jones nitrogen in a carbon slitlike pore

A density functional perturbation approximation based both on second-order perturbation theory and on the pore average density has been proposed to study the adsorption hysteresis of nitrogen in a carbon slit pore. The main advantage of the present approximation is that it is computationally much simpler than the original density functional approximation based on the second-order perturbation theory of liquids, and can be applied to several model fluids confined in a strong external field in order to study their structural and thermodynamic properties. The calculated adsorption hysteresis for the confined Lennard-Jones nitrogen is in very good agreement with computer simulation, even if its accuracy slightly deteriorates for the desorption branch. The calculated equilibrium particle density distributions also compare well with computer simulations, and are better than those of a density functional theory based on the so-called mean-field approximation.     

振动驱动颗粒气体体系的局域态本构关系的实验验证

对准二维、水平边界振动驱动的颗粒气体体系的流体力学 参量进行了局域态本构关系的实验研究. 实验观测结果与经典动力学理论预测进行了比较.由于颗粒气体空间分布的不均匀性, 颗粒体系的整体本构关系不成立, 有必要对局域态进行分析. 局域态本构关系是指颗粒系统的局域温度、局域压强和局域数密度之间的关系. 通过颗粒速度的方向变化, 可以得到颗粒的碰撞点. 因此在计算压力张量的对角线项时, 除了动力学部分之外, 我们计入了颗粒碰撞的影响, 得到了一个约为常数的压力张量迹, 即颗粒压强的空间分布, 与流体力学理论预测以及分子动力学模拟结果相符合; 但是颗粒温度和数密度的空间分布, 在振动的正反两个方向的分量出现差异, 并且温度、压强和数密度之间的局域本构关系, 无论在低密度或高密度区域, 实验与理论预测在定性上一致, 但定量上都有较大差别. 因此经典流体力学理论在描述这样的体系时需加以修正. 关键词： 颗粒气体 态方程 流体力学     

Structure and phase transitions in a network-forming associating Lennard-Jones fluid in a slit-like pore: a density functional approach

A study is reported of adsorption of an associating Lennard-Jones fluid with four associative sites per molecule in a slit-like pore. The density distribution of particles in the pore and thermodynamics properties are evaluated by using a density functional method. It is found that at low temperatures the fluid exhibits a set of layering transitions, followed by capillary condensation. Transitions are localized by analysing the grand canonical potential. The density profiles of particles and the distribution of unbound and differently bonded particles demonstrate changes in the structure of the fluid in the pore along the phase coexistence. The critical temperature is lower for a confined fluid, compared with the bulk counterpart. However, an increase in the energy of association increases the critical temperature. The envelope of the capillary condensation is narrower than the bulk liquid-vapour phase diagram. The dependence of the solvation force on the energy of association and on the bulk density is discussed.     

Density profile of water confined in cylindrical pores in MCM-41 silica

Recently, water absorbed in the porous silica material MCM-41-S15 has been used to demonstrate an apparent fragile to strong dynamical crossover on cooling below ～220 K, and also to claim that the density of confined water reaches a minimum at a temperature around 200 K. Both of these behaviours are purported to arise from the crossing of a Widom line above a conjectured liquid-liquid critical point in bulk water. Here it is shown that traditional estimates of the pore diameter in this porous silica material (of order 15 ?) are too small to allow the amount of water that is observed to be absorbed by these materials (around 0.5 g H(2)O/g substrate) to be absorbed only inside the pore. Either the additional water is absorbed on the surface of the silica particles and outside the pores, or else the pores are larger than the traditional estimates. In addition the low Q Bragg intensities from a sample of MCM-41-S15 porous silica under different dry and wet conditions and with different hydrogen isotopes are simulated using a simple model of the water and silica density profile across the pore. It is found the best agreement of these intensities with experimental data is shown by assuming the much larger pore diameter of 25 ? (radius 12.5 ?). Qualitative agreement is found between these simulated density profiles and those found in recent empirical potential structure refinement simulations of the same data, even though the latter data did not specifically include the Bragg peaks in the structure refinement. It is shown that the change in the (100) peak intensity on cooling from 300 to 210 K, which previously has been ascribed to a change in density of the confined water on cooling, can equally be ascribed to a change in density profile at constant average density. It is further pointed out that, independent of whether the pore diameter really is as large as 25 ? or whether a significant amount of water is absorbed outside the pore, the earlier reports of a dynamic crossover in supercooled confined water could in fact be a crystallization transition in the larger pore or surface water.     

Free energies of supercritical solvation from molecular dynamics simulation and integral equation studies

We present here molecular dynamics (MD) simulation and integral equation (IE) studies on free energies of solvation of a non-polar solute in a dilute supercritical solvent to estimate the contribution of inhomogeneities in solvent density to the free energy of solvation. The solvation of a Xe-like solute in an Ne-like solvent as well as that of naphthalene in CO₂ have been investigated. At state points in the compressible region in the neighborhood of the solvent critical point, we have utilized the IE estimates of free energies to model the ideal situation where local density inhomogeneities would be absent. The difference between the free energies in the presence (as derived from MD simulation) and in the absence (from IE) of local density inhomogeneities was studied as a function of density along an isotherm close to the critical point. Although for low density supercritical solvents, a marked difference is observed, a study of the density dependence of this difference across the critical density does not directly reveal any signature of local density enhancement on the thermodynamics of solvation.     

Collective behavior of thermally active colloids

Colloids with patchy metal coating under laser irradiation could act as local heat sources and generate temperature gradients that could induce self-propulsion and interactions between them. The collective behavior of a dilute solution of such thermally active particles is studied using a stochastic formulation. It is found that when the Soret coefficient is positive, the system could be described in a stationary state by the nonlinear Poisson-Boltzmann equation and could adopt density profiles with significant depletion in the middle region when confined. For colloids with a negative Soret coefficient, the system can be described as a dissipative equivalent of a gravitational system. It is shown that in this case the thermally active colloidal solution could undergo an instability at a critical laser intensity, which has similarities to a supernova explosion.     

Molecular dynamics study of confined ionic liquids in Au nanopore

Molecular dynamics simulations were carried to investigate the structure and dynamics of [BMIM][PF₆] ionic liquid (IL) confined inside a slit-like Au metal nanopore with a pore size of 5.0 nm. The calculations show that the mass and number densities of the confined ILs are oscillatory; the solid-like high density layers are formed in the vicinity of the metal surface. The orientational investigation shows that the imidazolium ring of [BMIM] cations prefers to form a small tilt angle with the pore walls. Furthermore, the mean squared displacement (MSD) calculation indicates that the dynamics of confined ILs are remarkably slower than those observed in bulk systems. Our results suggest that the confinement of the Au nanopore can strongly affect the structural and dynamical properties of the confined ILs.     

Computer simulation investigation of diffusion selectivity in graphite slit pores

Equilibrium molecular dynamics simulations are reported of oxygen and nitrogen molecules confined in graphite slit pores. Self- and collective diffusion coefficients have been calculated as a function of pore width, temperature and density for each pure component in the pore space. The aim of this study was to elucidate the mechanism by which oxygen and nitrogen are kinetically separated when air is passed over an adsorbent bed consisting of molecular sieving carbon in the commercial production of oxygen. It was found that a critical pore width exists for each species at which there is a sharp drop in the rate of diffusion (both self- and collective diffusion) of each fluid. The critical pore width is one for which the individual molecules are prevented from rotating freely about one of their axes. The greater length of a nitrogen molecule means that the critical pore width is higher for this species than for oxygen. Consequently, oxygen molecules diffuse substantially faster than nitrogen molecules in the vicinity of the nitrogen critical pore width. From an analysis of correlation functions and their corresponding power spectra it is shown that the restricted rotations, which occur at or below the critical pore width, cause a decoupling of translational and rotational modes, with the net result being a lowering of translational diffusion. The nitrogen critical pore width lies within the range of the mean pore size of most commercial molecular sieving carbons, and so this mechanism may help to explain the high oxygen selectivities reported in the literature.     

Surface critical behavior of fluids: Lennard-Jones fluid near a weakly attractive substrate

The phase behavior of fluids near weakly attractive substrates is studied by computer simulations of the coexistence curve of a Lennard-Jones (LJ) fluid confined in a slitlike pore. The temperature dependence of the density profiles of the LJ fluid was used to study the surface critical behavior. A universal critical behavior of the local order parameter, defined as the difference between the local densities of the coexisting liquid and vapor phases at some distance from the pore walls, , is observed in a wide temperature range and found to be consistent with the surface critical behavior of the Ising model. Near the surface the dependence of the order parameter on the reduced temperature obeys a scaling law ~1 with a critical exponent 1 of about 0.8, corresponding to the surface transition. A crossover from bulk-like to surface-like critical behavior occurs, when the distance to the surface is about twice the correlation length at the given temperature. Relations between the and transitions in Ising systems and the surface critical behavior of fluids are discussed.     

The phase behaviour of a fluid in a slitlike pore with permeable walls

The confinement of a lattice fluid in a set of slitlike pores separated by semipermeable walls with a finite width has been studied. The walls are modelled by a square-well repulsive potential with a finite height. The thermodynamic properties and the phase behaviour of the system are evaluated by means of Monte Carlo simulations. For some states theoretical calculations have been made using a mean-field-type theory. These investigations confirm previous findings for confined Lennard-Jones fluids, obtained from a density functional approach. For intermediate and low potential barriers that separate the pores, the isotherms exhibit two hysteresis loops and the liquid-vapour coexistence curve divides into two branches describing condensation inside the pore and inside the permeable wall. These two branches are separated by a triple point. At temperatures lower than the triple point temperature, the condensation takes place instantaneously in both the pore and inside the permeable wall. It was found that when the temperature is scaled by the bulk critical temperature, the phase diagram emerging from this simple mean-field treatment is close to the phase diagram obtained from simulation.     

Dewetting in associating lattice gas model confined by hydrophobic walls

The phase behavior of a two dimensional fluid confined within hydrophobic walls is obtained by Monte Carlo simulations. The fluid is described by the associating lattice gas model which reproduces the density and diffusion anomalous behavior of water.The confined fluid exhibits a liquid-liquid critical temperature which decreases with the decrease of the distance between the confining walls. In contact with the wall a dewetting is observed. The thickness of this interfacial layer is independent of the distance between the two walls. Even for very small distances between the two walls no total depletion is observed and consequently no drying transition is present.     

Lattice density functional theory for confined Ising fluids: comparison between different functional approximations in slit pore

ABSTRACT

In this paper, Lafuente and Cuesta's cluster density functional theory (CDFT) and lattice mean field approximation (LMFA) are formulated and compared within the framework of lattice density functional theory (LDFT). As a comparison, an LDFT based on our previous work on nonrandom correction to LMFA is also developed, where local density approximation is adopted on the correction. The numerical results of density distributions of an Ising fluid confined in a slit pore obtained from Monte Carlo simulation are used to check these functional approximations. Due to rational treatment on the coupling between site-excluding entropic effect and contact-attracting enthalpic effect by CDFT with Bethe-Peierls approximation (named as BPA-CDFT for short), the improvement of BPA-CDFT beyond LMFA is checked as expected. And it is interesting that our LDFT has a comparative accuracy with BPA-CDFT. Apparent differences between the profiles such as solvation force, excess adsorption quantity and interfacial tension from LMFA and non-LMFAs are found in our calculations. We also discuss some possible theoretical extensions of BPA-CDFT.     

Ionic liquid confined in silica nanopores: molecular dynamics in the isobaric–isothermal ensemble

Molecular dynamics simulations in the isobaric–isothermal ensemble are used to investigate the structure and dynamics of an ionic liquid confined at ambient temperature and pressure in hydroxylated amorphous silica nanopores. The use of the isobaric–isothermal ensemble allows estimating the effect of confinement and surface chemistry on the density of the confined ionic liquid. The structure of the confined ionic liquid is investigated using density profiles and structural order parameters while its dynamics is assessed by determining the mobility and ionic conductivity of the confined phase. Despite the important screening of the electrostatic interactions (owing to the small Debye length in ionic liquids), the local structure of the confined ionic liquid is found to be mostly driven by electrostatic interactions. We show that both the structure and dynamics of the confined ionic liquid can be described as the sum of a surface contribution arising from the ions in contact with the surface and a bulk-like contribution arising from the ions located in the pore centre; as a result, most properties of the confined ionic liquid are a simple function of the surface-to-volume ratio of the host porous material. In contrast, the ionic conductivity of the confined ionic liquid, which is a collective dynamical property, is found to be similar to the bulk. This study sheds light on the complex behaviour of hybrid materials made up of ionic liquid confined in inorganic porous materials.     

Oscillatory structure of confined fluids

The confinement of fluids in pores and wedges is associated with exponentially damped oscillatory packing structure, as observed with the surface forces apparatus. This paper reviews the statistical mechanics of confined fluids and then illustrates the results with density functional data for hard-sphere solvent. The free energy of the pore fluid and its functional derivatives with respect to thermodynamic fields all oscillate, as a function of pore width, with a wavelength close to the solvent diameter. In contrast, the density at the centre of pores oscillates with twice this wavelength, as a function of pore width. The development of a unified physics of confined fluids is considered. Approximations based on one-dimensional physics do extraordinarily well in planar symmetry at three-dimensions.     

Variational Monte Carlo analysis of Bose--Einstein condensation in a two-dimensional trap

The ground-state properties of a system with a small number of interacting bosons over a wide range of densities are investigated. The system is confined in a two-dimensional isotropic harmonic trap, where the interaction between bosons is treated as a hard-core potential. By using variational Monte Carlo method, we diagonalize the one-body density matrix of the system to obtain the ground-state energy, condensate wavefunction and the condensate fraction. We find that in the dilute limit the depletion of central condensate in the 2D system is larger than in a 3D system for the same interaction strength; however as the density increases, the depletion at the centre of 2D trap will be equal to or even lower than that at the centre of 3D trap, which is in agreement with the anticipated in Thomas--Fermi approximation. In addition, in the 2D system the total condensate depletion is still larger than in a 3D system for the same scattering length.     

Diffusivity and conductivity of a primitive model electrolyte in a nanopore

Equilibrium and non-equilibrium molecular dynamics simulations are applied to obtain the diffusion coefficient and electric conductivity of ions in dilute electrolytes confined in neutral cylindrical pores. The electrolyte is described with the restricted primitive model and the wall of the pore is modelled as a soft wall. The equilibrium molecular dynamics simulations show that the axial diffusion coefficient of ions decreases with increasing confinement. For a fixed pore radius the diffusion coefficient decreases with increasing number density of the ions. The current response of the system to an applied electric field is maintained at constant temperature by Gaussian isokinetic equations of motion, and at constant concentration by periodic boundary conditions with recycling of ions in the axial direction. The electric conductivity is calculated from the current density and the electric field applied for different pore sizes. In contrast to the trend in diffusivity, conductivity increases slightly in smaller pores. For a very small pore, however, conductivity is lower than the bulk, because oppositely charged ions moving in opposite directions under the electric field cannot avoid collisions with each other in a narrow channel.     

Grand Canonical Ensemble Monte Carlo Simulation of Depletion Interactions in Colloidal Suspensions

Depletion interactions in colloidal suspensions confined between two parallel plates are investigated by using acceptance ratio method with grand canonical ensemble Monte Carlo simulation. The numerical results show that both the depletion potential and depletion force are affected by the confinement from the two parallel plates. Furthermore, it is found that in the grand canonical ensemble Monte Carlo simulation, the depletion interactions are strongly affected by the generalized chemical potential.