狭缝滞留吸附性质的理论研究

针对Gibbs系综狭缝空间的单原子分子吸附系统, 引入了分子间作用对数的校正因子对格子空间的密度泛函理论(LDFT, lattice density functional theory)进行了修正. 利用连续迭代的方法, 探讨了狭缝滞留吸附行为的微观机理. 对比了修正前后的LDFT理论对滞留吸附行为的预测结果, 发现差别较大. 原LDFT理论由于平均场近似的处理, 对滞留吸附行为预测存在系统误差.     

氢气在A和X型沸石上超临界吸附的格子密度函数模型

采用基于三维Ono-Kondo方程的格子密度函数理论(LDFT)模型模拟了氢气在A和X型微孔沸石上的超临界吸附等温线. 根据沸石孔的尺寸和形状, LDFT模型将氢分子在孔中的吸附位分布近似为简单立方、面心立方和体心立方的团簇结构. 模拟结果表明, LDFT模型可有效地用于描述氢气在A和X型沸石上的单层或多层超临界吸附行为. 模拟所得的吸附等温线与实验测定结果吻合. 特别是, LDFT模型中的氢-沸石作用势能参数的准确性得到了Lennard-Jones(12-6)势能方法的有效验证. 因此, LDFT模型被用于预测了更宽温度和压力范围内氢气在X沸石上的超临界吸附.     

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

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

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

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

氢气在炭狭缝微孔内吸附的预测

通过边界的平均场近似,推导二平板狭缝孔格子理论Ono-Kondo吸附等温方程.利用格子气模型特性和微观物理学理论,计算氢分子在石墨平面的最大吸附容量.比较由氢分子在石墨平面二典型聚集状态标定的Ono-Kondo方程,并用预测精度较高的方程计算了与文献相同条件下的吸附等温线.在比较了计算结果、试验结果和GCMC分子模拟结果后,对Ono-Kondo吸附等温方程的特点、理论基础作了分析,指出了方程的适用范围.     

单壁碳纳米管上毒性苯气体净化的分子模拟(英文)

采用巨正则系综蒙特卡罗(GCMC)方法研究了空气中微量苯组分在单臂碳纳米管(SWNTs)上的吸附净化.模拟表明,具有较大孔径的(20,20)纳米管比较适合吸附纯苯蒸汽,而对于移除空气中的毒性苯物质,苯的吸附选择性分别在(12,12)纳米管及4.0 MPa时和(18,18)纳米管及0.1 MPa时出现最小值和最大值.为了解释这一异常行为,我们进一步分析了N2-O2-C6H6混合物的局部密度分布、吸附分子构型和概率密度分布,发现(18,18)纳米管内外完全被苯分子占据,而对于(12,12)纳米管,由于存在更强的吸附质-吸附剂相互作用,空气分子更倾向于吸附在管与管之间的间隙.此外,吸附分子的空间有序参数表明大多数苯分子采取"平躺"在纳米管表面的定位,而线性的N2和O2分子则多数平行于孔轴方向.最后研究了温度和苯分子主体相浓度对分离效果的影响.我们发现较大孔中的选择性随着温度的增加比小孔下降更加明显.与此对比,主体相苯浓度对小孔中的选择性起到更加重要的作用.     

单壁碳纳米管上毒性苯气体净化的分子模拟

采用巨正则系综蒙特卡罗(GCMC)方法研究了空气中微量苯组分在单臂碳纳米管(SWNTs)上的吸附净化. 模拟表明, 具有较大孔径的(20,20)纳米管比较适合吸附纯苯蒸汽, 而对于移除空气中的毒性苯物质, 苯的吸附选择性分别在(12,12)纳米管及4.0 MPa时和(18,18)纳米管及0.1 MPa时出现最小值和最大值. 为了解释这一异常行为, 我们进一步分析了N₂-O₂-C₆H₆混合物的局部密度分布、吸附分子构型和概率密度分布, 发现(18,18)纳米管内外完全被苯分子占据, 而对于(12,12)纳米管, 由于存在更强的吸附质-吸附剂相互作用, 空气分子更倾向于吸附在管与管之间的间隙. 此外, 吸附分子的空间有序参数表明大多数苯分子采取“平躺”在纳米管表面的定位, 而线性的N₂和O₂分子则多数平行于孔轴方向. 最后研究了温度和苯分子主体相浓度对分离效果的影响. 我们发现较大孔中的选择性随着温度的增加比小孔下降更加明显. 与此对比, 主体相苯浓度对小孔中的选择性起到更加重要的作用.     

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

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

芳香类化合物在ITQ-1分子筛中吸附特性的蒙特卡罗模拟

用巨正则统计系综蒙特卡罗模拟研究了苯、甲苯以及间二甲苯分子在ITQ-1分子筛中的吸附特性.从这3种分子的粒子分布云图上,可发现分子的扩散和吸附主要在十二元环超笼内发生,在十元环通道内的吸附和扩散则相对较难.从一系列不同压力下的蒙特卡罗模拟还预测了3种分子的吸附等温线,预测结果与实验结果相符.这3种分子在一定压力下,都可通过十元环通道或连接十二元环超笼的十元环窗口到达分子筛孔道内部,达到较好的吸附平衡状态.     

应用Polanyi吸附势理论描述超高交联树脂对多环芳烃的吸附

研究了萘、苊、芴三种多环芳烃类化合物在超高交联吸附树脂HC-1上的静态吸附平衡行为,分别采用Langmuir、Freundlich和Poanyi-Dubinin-Manes模型方程对吸附平衡数据进行了拟合分析,并探讨了吸附机理。实验结果表明,基于微孔填充机理的Polanyi-Dubinin-Manes方程能较好地拟合吸附平衡数据,表明树脂的微孔区对吸附起着重要作用。实验结果也表明,吸附质分子尺寸是影响吸附效果的重要因素,吸附质分子越大则吸附量越小;采用吸附质分子尺寸对Polanyi-Dubinin-Manes方程中吸附势密度进行修正,以体积吸附量对修正后吸附势密度作图,可得到一条与吸附质无关的特性曲线,根据该曲线方程可有效预测二、三环的多环芳烃化合物在HC-1上的吸附性能。     

Modeling fluid diffusion using the lattice density functional theory approach: counterdiffusion in an external field

The dependence of the diffusivity on temperature, pressure, and composition is not understood well; consequently, data is preferred significantly over correlations in most practical situations. Even in dilute gases, the contributions of attractions and repulsions to the diffusivity are difficult to understand on a molecular level without performing simulations. We have developed a Lattice Density Functional Theory (LDFT) approach for modeling diffusion to supplement existing methods that are very rigorous but computationally demanding. The LDFT approach is analogous to the van der Waals equation in how it accounts for molecular interactions in that it has first-order corrections to ideal behavior; it is an extension of the Equilibrium LDFT for adsorption and phase behavior. In this work, the LDFT approach is presented and demonstrated by modeling the problem of color counterdiffusion in an externally-applied potential field. This potential field, in combination with the intermolecular potential function, creates a diffusion regime in which repulsions cause oscillations in the density profile. Using the LDFT approach, the oscillations were described and attributed to nearest-neighbor and next nearest-neighbor interactions. The LDFT approach gives qualitative and quantitative agreement with dual control-volume Grand Canonical Molecular Dynamics simulations.     

Modeling of monolayer adsorption of hydrogen on ZSM‐5 zeolite by lattice density functional theory

The adsorption isotherms of hydrogen on microporous zeolite ZSM‐5, at supercritical conditions, have been modeled using the monolayer lattice density functional theory (LDFT) models, where the simple cubic lattice, face‐centered cubic lattice, body‐centered cubic lattice and tetragonal lattice structures are assumed for the arrangements of the adsorption sites inside pores based on the size and shape of the zeolite. The results indicate that the monolayer LDFT models appear to be effective in describing hydrogen adsorption on zeolite ZSM‐5 at supercritical conditions, and the calculated adsorption isotherms agree well with the experimental isotherms measured previously. The layer density of adsorbed phase is presented versus the bulk density and temperature. It is found that the densities of adsorbed phase on adsorbent surface are much higher than the bulk density for temperature range under study. However, in the core region, the layer densities are close to the bulk density. The monolayer adsorption is suitable for hydrogen on ZSM‐5 zeolite. Copyright © 2013 John Wiley & Sons, Ltd.     

Effect of spin polarization for hydrogen adsorbed on Si(111)(1×1) surface: First‐principles calculations

The role of spin polarization on adsorption of atomic and molecular hydrogen on Si(111)(1×1) surface is examined by comparing the results of the local spin density approximation (LSD) and those of the local density approximation (LDA). A large improvement of the adsorption energies (around 0.8 eV/H) was found for the H atom adsorbed on Si(111)(1×1) surface. The inclusion of spin polarization reduces the overbinding between the H atom and the silicon surface and its effect is much more pronounced when the H atom is far away from the surface. Despite of the large changes in the adsorption energies, the main character of the potential energy surface of the H atom on Si(111)(1×1) surface is retained. An opposite effect is found in the charge‐density‐transfer map of LSD results as compared to LDA results for the H atom approaching the surface through the H3 path, in which the H atom loses electrons rather than gains electrons from the surface. The fact that the H atom tends to lose electrons in the silicon bulk has already been reported by the experimental studies for the behavior of the H atom in the p‐type silicon. For the molecular hydrogen on Si(111)(1×1) surface, the effect of the spin polarization is so small that it can be neglected. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 79: 47–55, 2000     

Principles of the theory of adsorption of large molecules in pores with nonuniform walls

Principles of the theory of adsorption of large molecules blocking more than one adsorption center on a surface in slit-like pores are proposed. The theory takes into account lateral adsorbate—adsorbate interactions and nonuniformity of the pore walls. The equations of adsorption isotherms are derived using the cluster approach. The lateral interactions are taken into account in the quasi-chemical approximation, preserving effects of direct correlations, and in the mean field approximation without effects of correlations. The following problems are discussed: 1) distinguishing of partial contributions of nonuniform adsorption centers on the pore walls; 2) exact solution for dimer adsorption in a two-layer pore with uniform walls; 3) basic types of adsorption isotherms, for which the differences are due to various orientations of the adsorbate in micropores with uniform walls; 4) estimates of the pressure responsible for volume filing of micropores; and 5) the effect of nonuniformity of the pore walls on the pressure values. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 8, pp. 1467–1478, August, 1999.     

Effect of critical fluctuations on adsorption of van der Waals fluid in a spherical cavity

A recently proposed non-uniform fifth-order thermodynamic perturbation theory (TPT) is employed to investigate the adsorption of a hard core attractive Yukawa (HCAY) fluid in a spherical cavity. Extensive comparison with available simulation data indicate that the non-uniform fifth-order TPT is sufficiently reliable in calculating the density profiles of the HCAY fluid in the highly confining geometry, and generally is more accurate than a previous third-order?+?second-order perturbation density functional theory. The non-uniform fifth-order TPT is free from numerically solving an Ornstein–Zernike integral equation, and also free of any adjustable parameter; consequently, it can be applied to both supercritical and subcritical temperature regions. The non-uniform fifth-order TPT is employed to investigate critical adsorption of the HCYA fluid in a single spherical cavity – it is disclosed that the critical fluctuations near the critical point induce depletion adsorption – quantitative theoretical calculation on relationship between the critical depletion adsorption, parameters of coexistence bulk phase and the responsible external field is in agreement with qualitative physical analysis.     

过渡金属掺杂的g-GaN吸附Cl2和CO气体分子的第一性原理研究

基于密度泛函理论的第一性原理计算, 系统研究了类石墨烯氮化镓(g-GaN)和掺杂过渡金属原子(TM)的 g-GaN 对 Cl₂和CO气体分子的吸附行为。结果表明, Cl₂和 CO在本征 g-GaN上的吸附均为物理吸附, 2个体系的吸附能均为正值, 表明体系不稳定。相反, Cl₂和 CO在 Fe和 Co掺杂的 g-GaN上吸附时的吸附能为负值, 且吸附能较小, 表明吸附体系稳定。通过分析态密度、电荷密度差和能带结构等性质, 可以得出结论:过渡金属原子的引入能有效增强气体分子与 g-GaN之间的相互作用。     

Density Functional Study of the C Atom Adsorption on the α-Fe2O3 （001） Surface

The adsorption of C atoms on the α-Fe2O3(001) surface was studied based on density function theory(DFT) ,in which the exchange-correlation potential was chosen as the PBE(Perdew,Burke and Ernzerhof) generalized gradient approximation(GGA) with a plane wave basis set. Upon the optimization on different adsorption sites with coverage of 1/20 and 1/5 ML,it was found that the adsorption of C atoms on the α-Fe2O3(001) surface was chemical adsorption. The coverage can affect the adsorption behavior greatly. Under low coverage,the most stable adsorption geometry lied on the bridged site with the adsorption energy of about 3.22 eV;however,under high coverage,it located at the top site with the energy change of 8.79 eV. Strong chemical reaction has occurred between the C and O atoms at this site. The density of states and population analysis showed that the s,p orbitals of C and p orbital of O give the most contribution to the adsorption bonding. During the adsorption process,O atom shares the electrons with C,and C can only affect the outermost and subsurface layers of α-Fe2O3;the third layer can not be affected obviously.