Effect of confinement in nano-porous materials on the solubility of a supercritical gas

By combining Gibbs Ensemble Monte Carlo simulations and density functional theory, we investigate the influence of confinement in a slit-shaped carbon pore on the solubility of a supercritical solute gas in a liquid solvent. In the cases studied here, competing adsorption of the solvent and solute determines whether the solubility is enhanced or suppressed for larger pores. We find that the solubility in the confined system is strongly dependent on pore width, and that molecular packing effects are important for small pore widths. In addition, the solubility decreases on increase in the temperature, as for the bulk mixture, but the rate of decrease is greater in the pore due to a decrease in the partial molar enthalpy of the solute in the pore; this effect becomes greater as pore width is decreased. The solubility is increased on increasing the bulk pressure of the gas in equilibrium with the pore, and obeys Henry's law at lower pressures. However, the Henry constant differs significantly from that for the bulk mixture, and the range of pressure over which Henry's law applies is reduced relative to that for the bulk mixture. The latter observation indicates that solute–solute interactions become more important in the pore than for the bulk at a given bulk pressure. Finally, we note that different authors use different definitions of the solubility in pores, leading to some confusion over the reported phenomenon of ‘oversolubility’. We recommend that solubility be defined as the overall mole fraction of solute in the pores, since it takes into account the increase in density of the solvent in the pores, and avoids ambiguity in the definition of the pore volume.     

Instability of the Randall-Sundrum Model and Exact Bulk Solutions

Five dimensional geodesic equation is used to study the gravitational force acted on a test particle in the bulk of the Randall-Sundrum two-brane model. This force could be interpreted as the gravitational attraction from matters on the two branes and may cause the model to be unstable. By analogy with star models in astrophysics, a fluid RS model is proposed in which the bulk is filled with a fluid and this fluid has an anisotropic pressure to balance the gravity from the two branes. Thus a class of exact bulk solutions is obtained which shows that any 4D Einstein solution with a perfect fluid source can be embedded in y = constant hypersurfaces in the bulk to form an equilibrium state of the brane model. By requiring a 4D effective curvature to have a minimum, the compactification size of the extra dimension is discussed.     

方解石纳米孔隙内二氧化碳毛细凝聚的分子模拟

在二氧化碳地质封存、增产非常规油气以及孔隙材料表征测量方面,纳米孔隙中二氧化碳相态的准确预测具有重要意义。然而,由于纳米尺度下毛细力、分离压等作用力占据主导因素,流体在孔隙中的相行为与体相流体存在根本不同。已有实验和模拟表明,Kelvin毛细凝聚理论无法预测特征尺度10nm下的,孔隙内流体凝聚压力与体相饱和蒸气压的偏离程度。本文利用分子模拟方法,研究了孔径范围为0.83-8.0nm的方解石纳米孔隙中二氧化碳毛细凝聚。结果表明,微孔(小于2nm)中二氧化碳受吸附层影响,凝聚压力远低于体相饱和蒸气压,介孔(2?50nm)中二氧化碳受壁面影响,凝聚压力仍低于体相饱和蒸气压,且受孔径影响显著,并获得了纳米孔隙中二氧化碳相态随压力的变化规律。     

Probing interfacial dynamics of water in confined nanoporous systems by NMRD

The confined dynamics of water molecules inside a pore involves an intermittence between adsorption steps near the interface and surface diffusion and excursions in the pore network. Depending on the strength of the interaction in the layer(s) close to the surface and the dynamical confinement of the distal bulk liquid, exchange dynamics can vary significantly. The average time spent in the surface proximal region (also called the adsorption layer) between a first entry and a consecutive exit allows estimating the level of ‘nanowettablity’ of water. As shown in several seminal works, NMRD is an efficient experimental method to follow such intermittent dynamics close to an interface. In this paper, the intermittent dynamics of a confined fluid inside nanoporous materials is discussed. Special attention is devoted to the interplay between bulk diffusion, adsorption and surface diffusion on curved pore interfaces. Considering the nano or meso length scale confinement of the pore network, an analytical model for calculating the inter-dipolar spin–lattice relaxation dispersion curves is proposed. In the low-frequency regime (50?KHz–100?MHz), this model is successfully compared with numerical simulations performed using a 3D-off lattice reconstruction of Vycor glass. Comparison with experimental data available in the literature is finally discussed.     

Demixing of a binary mixture in slit-like pores at high temperatures

We investigate the fluid—fluid demixing transitions in the case of adsorption of so-called symmetric binary mixtures in slit-like pores at temperatures higher than the bulk gas—liquid critical temperature. The aim of the study is to determine how the demixing of such mixtures in the pores depends on the bulk phase composition and on the parameters characterizing the pore. The calculations have been carried out by means of a density functional theory. In the case of an equimolar bulk mixture, the demixing transition inside the pore occurs only when the adsorption potentials of both species are identical. The occurrence of this transition is manifested by a cusp in the adsorption isotherm. For nonequimolar bulk phase compositions, the transition can also take place if the adsorption energies of both components are different. However, the difference in the adsorption energies should be small enough, otherwise a continuous demixing takes place. For non-equimolar compositions two branches of the grand canonical potential intersect, whereas for equimolar bulk composition they meet tangentially. We have determined phase behaviour for several model systems.     

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.     

The phase coexistence method to obtain surface free energies and nucleation barriers: a brief review

ABSTRACT

A recently developed method where one analyses the finite size effects associated with liquid–solid phase equilibria including vapour–crystal coexistence is briefly reviewed. It is shown that the estimation of the chemical potential of the vapour surrounding the crystal as function of the crystal volume yields information on the bulk coexistence conditions, when an extrapolation to the thermodynamic limit is performed. Estimating the pressure of the fluid surrounding the crystal nucleus in the finite simulation box and the volume of this nucleus that coexists with the fluid in thermal equilibrium, an estimate for the total surface excess free energy can be obtained, which to a very good approximation is independent of the size of the simulation box. The free energy barrier against homogeneous nucleation of crystals thus can be estimated as a function of the nucleus volume. Monte Carlo simulations for the soft effective Asakura–Oosawa model of colloid-polymer mixtures which form face-centered cubic colloidal crystals are used to exemplify this method, computing the surface excess free energy of these crystals over a wide range of crystal volumes, without the need to characterise the non-spherical crystal shape. A possible extension of these concepts to heterogeneous nucleation is also briefly discussed.     

A solvable model for localized adsorption in a coulomb system

A model for an interface with localized adsorption is presented, in which the surface has a distribution of sticky adhesive sites in contact with a Coulomb fluid. Contrary to the current literature on the electrical double layer the surface charge is in dynamic equilibrium with the bulk fluid. The sum rules obeyed by the one- and two-body correlation functions are investigated. Explicit results are obtained for a solvable model, the two-dimensional one-component plasma at reduced temperature 2. The effect of the granularity of the adsorbed charge on the adsorption isotherm is discussed.     

气泡体积分数对沙质沉积物低频声学特性的影响

由于有机物质分解等原因,实际的海底沉积物中存在气泡,气泡的存在会显著影响沉积物低频段的声学特性,因此研究气泡对沉积物低频段声速的影响机理具有重要意义.考虑到外场环境的不可控性,在室内水池中搭建了大尺度含气非饱和沙质沉积物声学特性获取平台,在有界空间中应用多水听器反演方法首次获取了含气非饱和沙质沉积物300—3000 Hz频段内的声速数据(79—142 m/s),并同时利用双水听器法获取了同一频段的数据(112—121 m/s).在声波频率远低于沉积物中最大气泡的共振频率时,根据等效介质理论,将孔隙水和气泡等效为一种均匀流体,改进了水饱和等效密度流体近似模型.模型揭示了气泡对沉积物低频段声学特性的影响规律,理论上解释了沉积物中声速的降低.通过分析模型预报声速对模型参数的敏感性,根据测量得到的声速分布反演得到了沉积物不同区域的气泡体积分数,气泡体积分数从1.07%变化到2.81%.改进的模型为沉积物中气泡体积分数估计提供了一种新方法.     

FP-LMTO calculations of the structural,elastic, thermodynamic,and electronic properties of the ideal-cubic perovskite BiGaO3

Using the first-principles full-potential linear muffin-tin orbital method within the local density approximation, we have studied the structural, elastic, thermodynamic, and electronic properties of the ideal-cubic perovskite BiGaO₃. It is found that this compound has an indirect band gap. The valence band maximum (VBM) is located at Γ-point, whereas the conduction band minimum (CBM) is located at X-point. The pressure and volume dependences of the energy band gaps have been calculated. The elastic constants at equilibrium are also determined. We derived the bulk and shear moduli, Young’s modulus, and Poisson’s ratio. The thermodynamic properties are predicted through the quasi-harmonic Debye model, in which the lattice vibrations are taken into account. The variation of the bulk modulus, heat capacities, and Debye temperature with pressure and temperature are successfully obtained.     

Measuring nanopore size from the spin-lattice relaxation of CF4 gas

The NMR 19F spin-lattice relaxation time constant T1 for CF4 gas is dominated by spin-rotation interaction, which is mediated by the molecular collision frequency. When confined to pores of approximately the same size or smaller than the bulk gas mean free path, additional collisions of molecules with the pore walls should substantially change T1. To develop a method for measuring the surface/volume ratio S/V by measuring how T1 changes with confinement, we prepared samples of known S/V from fumed silica of known mass-specific surface area and compressed to varying degrees into cylinders of known volume. We then measured T1 for CF4 in these samples at varying pressures, and developed mathematical models for the change in T1 to fit the data. Even though CF4 has a critical temperature below room temperature, we found that its density in pores was greater than that of the bulk gas and that it was necessary to take this absorption into account. We modeled adsorption in two ways, by assuming that the gas condenses on the pore walls, and by assuming that gas in a region near the wall is denser than the bulk gas because of a simplified attractive potential. Both models suggested the same two-parameter formula, to which we added a third parameter to successfully fit the data and thus achieved a rapid, precise way to measure S/V from the increase in T1 due to confinement in pores.     

Covariance for cone and wedge complete filling

Interfacial phenomena associated with fluid adsorption in two dimensional systems have recently been shown to exhibit hidden symmetries, or covariances, which precisely relate local adsorption properties in different confining geometries. We show that covariance also occurs in three-dimensional systems and is likely to be verifiable experimentally and in Ising model simulations studies. Specifically, we study complete wetting in wedge (W) and cone (C) geometries as bulk coexistence is approached and show that the equilibrium midpoint heights satisfy l(c)(h,alpha)=l(w)(h / 2,alpha), where h measures the partial pressure and alpha is the tilt angle. This covariance is valid for both short-ranged and long-ranged intermolecular forces and identifies both leading and next-to-leading-order critical exponents and amplitudes in the confining geometries.     

Adsorption of argon and nitrogen in cylindrical pores of MCM-41 materials: application of density functional theory

Adsorption of argon and nitrogen at their respective boiling points in cylindrical pores of MCM-41 type silica-like adsorbents is studied by means of a non-local density functional theory (NLDFT), which is modified to deal with amorphous solids. By matching the theoretical results of the pore filling pressure versus pore diameter against the experimental data, we arrive at a conclusion that the adsorption branch (rather than desorption) corresponds to the true thermodynamic equilibrium. If this is accepted, we derive the optimal values for the solid-fluid molecular parameters for the system amorphous silica-Ar and amorphous silica-N₂, and at the same time we could derive reliably the specific surface area of non-porous and mesoporous silica-like adsorbents, without a recourse to the BET method. This method is then logically extended to describe the local adsorption isotherms of argon and nitrogen in silica-like pores, which are then used as the bases (kernel) to determine the pore size distribution. We test this with a number of adsorption isotherms on the MCM-41 samples, and the results are quite realistic and in excellent agreement with the XRD results, justifying the approach adopted in this paper.     

Further observations on the polarized absorption spectrum of benzene crystals at 2000 Å

The structure of a hard sphere fluid confined by model slit and cylindrical pores is investigated. Results from grand canonical Monte Carlo (GCMC) simulations and from the hypernetted chain/mean spherical approximation (HNC/MSA) equation are reported. GCMC results are compared with those from the HNC/MSA equation, and agreement is good. The effect of confinement on liquids at the same chemical potentials is that the absorption of the hard sphere fluid into the pores decreases with increasing confinement, i.e., when going from planar to cylindrical geometry or by narrowing the pores. The adsorption on the pore walls has, in general, the opposite behaviour. For high bulk concentrations and certain values of cylindrical pore diameter the concentration profile is higher at the centre of the pore than next to the pore wall. A very strong, but continuous, transition occurs in the concentration profile, as a function of the cylinder's diameter. These results could be of some interest in catalysis studies.     

First-Principles Calculations of Elastic Properties of LaNi5 Compound

The equilibrium lattice constants, temperature dependence of bulk modulus, the pressure dependence of the normalized volume V/V0, elastic constants Cij and bulk modulus of LaNi5 crystal are obtained using the firstprincipies piane-wave pseudopotential method in the GGA-PBE generalized gradient approximation as well as the quasi-harmonic Debye model. We analyse the relationship between bulk modulus and temperature up to 2000 K and obtain the relationship between bulk modulus B and pressure at diFFerent temperatures. It is found that the bulk modulus B increases monotonously with increasing pressure. Moreover, the pressure dependences of Debye temperatures and the pressure derivatives of lattice constants are also successfully obtained. The calculated results are in agreement with the experimental data and the other theoretical results.     

Adsorption and interfacial phenomena of a Lennard-Jones fluid adsorbed in slit pores: DFT and GCMC simulations

ABSTRACT

Confinement of fluids in porous media leads to the presence of solid–fluid (SF) interfaces that play a key role in many different fields. The experimental characterisation of SF interfacial properties, in particular the surface tension, is challenging or not accessible. In this work, we apply mean-field density functional theory (DFT) to determine the surface tension and also density profile of a Lennard-Jones fluid in slit-shaped pores for realistic amounts of adsorbed molecules. We consider the pore walls to interact with fluid molecules through the well-known 10-4-3 Steele potential. The results are compared with those obtained from Monte Carlo simulations in the Grand Canonical Ensemble (GCMC) using the test-area method. We analyse the effect on the adsorption and interfacial phenomena of volume and energy factors, in particular, the pore diameter and the ratio between SF and fluid–fluid dispersive energy parameters, respectively. Results from DFT and GCMC simulations were found to be comparable, which points to their reliability.     

Theoretical studies of argon adsorption in MCM-41 mesoporous systems

Adsorption isotherms are predicted for spherical adsorbates in cylindrical channels of MCM-41 mesoporous materials over a wide range of temperatures by using the “fragment method”. This prediction shows that an equilibrium capillary condensation is impossible for pores with diameters smaller than 2.5 nm. The adsorbate distribution in relatively large pore channels was described by the quasi-chemical approximation (QCA) that takes into account direct pair correlations between interacting molecules. In order to improve the lattice-gas model in the vicinity of the critical point, a calibration function that takes into account information from the fragment method, was introduced into the QCA equations. The influence of the size factor of pores on argon adsorption isotherms was demonstrated.     

含少量气泡流体饱和孔隙介质中的弹性波

考虑孔隙流体中含有少量气泡,且气泡在声波作用下线性振动,研究声波在这种孔隙介质中的传播特性.本文先由流体质量守恒方程和孔隙度微分与流体压力微分的关系推导出了含有气泡形式的渗流连续性方程;在处理渗流连续性方程中的气体体积分数时间导数时,应用Commander气泡线性振动理论导出气体体积分数时间导数与流体压强时间导数的关系,进而得到了修正的Biot形式的渗流连续性方程;最后结合Biot动力学方程求得了含气泡形式的位移场方程,便可得到两类纵波及一类横波的声学特性.通过对快、慢纵波的频散、衰减及两类波引起的流体位移与固体位移关系的考察,发现少量气泡的存在对快纵波和慢纵波的传播特性影响较大.     

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.     

MD模拟方法研究圆柱形纳米微孔的吸附平衡

本文用MD（分子动力学）模拟方法对圆柱形纳米微孔中的Lennard-Jones流体的物理吸附平衡及界面现象进行了研究。模拟中引入新型势模型,采用了计数法计算平衡压力,预测了在不同孔径下（直径为 3～5nm）的吸附等温线,分析了毛细凝聚理论用于纳米微孔吸附的局限性。并比较了同一纳米尺度的窄缝形微孔和圆柱形微孔吸附的差异。