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

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

微孔中简单流体混合物扩散系数的分子动力学模拟与关联

为了了解简单流体混合物在微孔介质中的流动和传递性质,对微孔中氩和氪流体混合物的扩散系数进行了计算机模拟和关联模型研究.运用平衡分子动力学方法模拟了宏量条件下饱和氩流体的扩散系数和恒温氪流体的扩散系数,模拟值与文献实验值符合良好,从而程序的正确性得到验证.然后,采用类似Bitsanis等人的方法模拟了平板湿壁微孔中氩和氪等摩尔流体混合物在不同对比温度、不同对比密度以及不同对比孔径条件下的扩散系数,发现孔径很小的时候扩散系数会急剧的增大.同时基于这些模拟值,参考CE理论和Heyes关系式,以对比温度、对比密度以及对比孔径为变量,关联出两个简单流体等摩尔混合物在微孔中扩散系数的计算模型.模型的计算结果与计算机模拟值能够较好地吻合.     

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

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

用密度泛函理论研究Lennard-Jones 流体在狭缝中的相平衡

用改进的基础度量理论(modified fundamental measure theory, MFMT)和密度Taylor展开分别表达过剩自由能中的短程作用和色散作用. 流体分子与狭缝壁之间的相互作用以10-4-3势能函数表达. 由巨势最小原理确定Lennard-Jones (LJ)流体在狭缝中的密度分布和过剩吸附量, 所得结果与分子模拟数据吻合良好. 根据平衡时两相温度, 化学势及巨势相等, 计算了LJ流体在狭缝中的相平衡.     

用密度泛函理论研究Lennard-Jones 流体在狭缝中的相平衡

用改进的基础度量理论(modified fundamental measure theory, MFMT)和密度Taylor展开分别表达过剩自由能中的短程作用和色散作用. 流体分子与狭缝壁之间的相互作用以10-4-3势能函数表达. 由巨势最小原理确定Lennard-Jones (LJ)流体在狭缝中的密度分布和过剩吸附量, 所得结果与分子模拟数据吻合良好. 根据平衡时两相温度, 化学势及巨势相等, 计算了LJ流体在狭缝中的相平衡.     

分子在纯硅β分子筛内扩散的随机行走模型

建立了一个β分子筛上分子扩散的模型. 该模型中, 分子在β分子筛中运动是在不同吸附点位上作无规行走. β孔道的拓扑结构和在两种孔道吸附位上不同的跃迁几率导致分子沿两个主轴方向扩散, 扩散系数存在一个关联关系; 分子动力学对不同温度下苯分子在β分子筛上扩散模拟证实了这一关联关系. 氩原子在不同作用半径下的动力学模拟表明, 分子作用半径大小是满足随机行走假设的重要条件, 即该模型要求扩散分子作用半径足够大, 与分子筛孔径相近.     

ITQ-3分子筛中二元混合分子扩散性质的模拟

采用平衡分子动力学方法先是模拟研究了纯组分Ar、SF_6以及CF4在ITQ-3分子筛中的扩散行为,结果表明,在窄孔道中Ar的扩散系数随负载的增加呈现先增加再减小的变化趋势,在宽孔道中则随负载逐渐减小。纯组分的大分子SF_6以及CF4只在y方向扩散,扩散系数与小分子在z方向的变化规律一致。随后,又模拟了二元混合组分SF_6和CF4在ITQ-3分子筛中的扩散行为,模拟结果与各自纯组分SF_6和CF4进行了对比,发现二元混合的SF_6和CF4整体变化趋势与其单组分保持一致,在y方向的扩散系数都是随着负载的增加而先增大后减小,但是两者之间还存在一种相互作用,使得在整个负载范围内扩散较快的CF4的扩散系数比纯组分的扩散系数相对小一些,而纯组分时扩散相对缓慢的SF_6在二元混合中的扩散系数则相对增大。     

分子在ITQ-3分子筛窄孔道内扩散的过渡态理论模型

建立了一个基于过渡态理论的分子在ITQ-3窄孔道方向扩散的模型. 该模型中, 由于分子从空腔中的一个吸附位越过势垒到相邻的另一个空腔中的吸附位需要克服很大的势垒能, 因而分子在势垒处可以简化处理为只存在排斥势, 可得到扩散系数依赖温度和Lennard-Jones作用参数的解析关系. 用分子动力学方法对CF₄在ITQ-3上扩散进行了模拟并验证了解析关系的合理性. 分子动力学模拟计算得到的扩散活化能、势垒能和吸附位势能与实际值相吻合. 模拟结果也显示了扩散系数依赖于附载量, 表现为先增大后减小的变化模式. 扩散活化能的计算证实了这一变化机理, 即当附载量增加时, 由于处于空腔中的分子彼此抬高了势能, 降低了扩散活化能, 使得扩散系数随附载量的增加而增大, 之后由于堵塞效应, 扩散系数随附载量的增加而逐步减小.     

巨正则系综Monte Carlo模拟方法确定活性炭的微孔尺寸

根据299 K下甲烷在活性炭中的吸附实验数据，通过调节狭缝微孔的孔宽参数，利用巨正则系综Monte Carlo（GCEMC）方法得到不同孔宽下流体的微观结构以及吸附等温线．比较并拟合模拟结果和实验数据，确定了活性炭微孔的平均孔宽，为下一步求解微孔尺寸分布以及为预测吸附剂在不同温度下吸附不同吸附质分子时的吸附性能提供了基础与指导．模拟中，甲烷分子采用单点Lennard-Jones球型分子模型，活性炭用狭缝孔来近似表征，流体分子与单个狭缝墙的相互作用采用著名的Steele的10－4-3势能模型．模拟表明，此方法为考察介孔材料的微孔分布以及微孔平均孔宽提供了新的思路．     

巨正则系综Monte Carlo模拟方法研究活性炭的微孔尺寸

根据299K下甲烷在活性炭中的吸附实验数据，通过调节狭缝微孔的孔宽参数，利用巨正则系综Monte Carlo（GCEMC）方法得到不同也宽下流体的微观结构以及吸附等温线，比较并拟合模拟结果和实验数据，确定了活性炭微孔的平均孔宽，为下一步求解微孔尺寸分布以及为预测吸附剂在不同温度下吸附不同吸附质分子的吸附性能提供了基础与指导，模拟，甲烷分子采用单点Lennard-Jones球型分子模型，活性炭用狭缝孔来近似表征，流体分子与单个狭缝墙的相互作用采用著名的Steele的10－4－3势能模型，模拟表明，此方法为考察介孔材料的微孔分布以及微孔平均孔宽提供了新的思路。     

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.     

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.     

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.     

Surface and porosity of nanocrystalline boehmite xerogels

Boehmite xerogels are prepared by hydrolysis of Al(OC4H9)3 followed by peptization with HNO3 (H+/Al = 0, 0.07, 0.2). XRD and TEM show that these gels are made of nanosized crystals (5-9 nm in width and 3 nm thick). According to the amount of acid, no significant differences are found in size and shape, but only in the spatial arrangement of the crystallites. Nitrogen adsorption-desorption isotherms of nonpeptized gels are of type IV, whereas isotherms of peptized gels are of type I. These isotherms are analyzed by the t-plot method. The majority of pore volume results from intercrystalline mesopores, but the peptized gels also contain intercrystalline micropores. The particle packing is very dense for the gel peptized with H+/Al = 0.2 (porosity = 0.26), but it is less dense in non-peptized gel (porosity = 0.44). Heating these gels under vacuum creates, from 250 degrees C onwards, an intracrystalline microporosity resulting from the conversion of boehmite into transition alumina. But heating also causes intercrystalline micropores collapsing. The specific surface area increases up to a limit temperature (300 degrees C for nonpeptized gels and 400 degrees C for peptized) beyond which sintering of the particles begins and the surface decreases. The PSD are calculated assuming a cylindrical pore geometry and using the corrected Kelvin equation proposed by Kruk et al. Peptized xerogels give a monomodal distribution with a maximum near 2 nm and no pores are larger than 6 nm. Nonpeptized gels have a bimodal distribution with a narrow peak near to 2 nm and a broad unsymmetrical peak with a maximum at 4 nm. Heating in air above 400 degrees C has a strong effect on the porosity. As the temperature increases, there is a broadening of the distribution and a marked decrease of small pores (below 3 nm). However, even after treatment at 800 degrees C, micropores are still present.     

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.     

Effects of fiber diameter and CO2 activation temperature on the pore characteristics of polyacrylonitrile based activated carbon nanofibers

The effects of fiber diameter and activation temperature on the pore characteristics of polyacrylonitrile based activated carbon nanofibers are investigated. It was found that lower fiber diameters as well as higher activation temperatures lead to a higher weight loss, specific surface area and total pore volume. The nitrogen adsorption capacity of activated carbon nanofibers is almost three times that of activated carbon fiber with a diameter of 10 µm. As far as the size of pores in activated carbon nanofibers is concerned, it is basically the micropores that dominate the scene. Moreover, tailoring the pore characteristics by adjusting the activation temperature and fiber diameter is plausible. Copyright © 2009 John Wiley & Sons, Ltd.     

Tailoring of nanoscale porosity in carbide-derived carbons for hydrogen storage

The poor performance of hydrogen storage materials continues to hinder development of fuel cell-powered automobiles. Nanoscale carbons, in particular (activated carbon, exfoliated graphite, fullerenes, nanotubes, nanofibers, and nanohorns), have not fulfilled their initial promise. Here we show that carbon materials can be rationally designed for H2 storage. Carbide-derived carbons (CDC), a largely unknown class of porous carbons, are produced by high-temperature chlorination of carbides. Metals and metalloids are removed as chlorides, leaving behind a collapsed noncrystalline carbon with up to 80% open pore volume. The detailed nature of the porosity-average size and size distribution, shape, and total specific surface area (SSA)-can be tuned with high sensitivity by selection of precursor carbide (composition, lattice type) and chlorination temperature. The optimum temperature is bounded from below by thermodynamics and kinetics of chlorination reactions and from above by graphitization, which decreases SSA and introduces H2-sorbing surfaces with binding energies too low to be useful. Intuitively, pores of different size and shape should not contribute equally to hydrogen storage. By correlating pore properties with 77 K H2 isotherms from a wide variety of CDCs, we experimentally confirm that gravimetric hydrogen storage capacity normalized to total pore volume is optimized in materials with primarily micropores ( approximately 1 nm) rather than mesopores. Thus, in agreement with theoretical predictions, a narrow size distribution of small pores is desirable for storing hydrogen, while large pores merely degrade the volumetric storage capacity.     

A diffusion model for the fluids confined in micropores

The self-diffusion coefficients were calculated by molecular dynamics simulations and the effects of pore width, temperature, and fluid density on diffusion behavior of simple fluid argon and polar fluid water confined in micropores were analyzed and studied. A mathematical model describing diffusion behavior of fluids confined in micropores was proposed from the theories of molecular dynamics and molecular kinematics, and validated on the basis of the simulation results at various conditions. The model indicates that the diffusion coefficient is proportional to the square root of the pore width and to the temperature divided by the density squared. It is applicable to either liquid or gas states and only two parameters are required.