约束条件下的硬球流体

利用密度泛函理论和分子动力学方法,对处于两平行硬墙之间的硬球流体的密度分布进行了计算通过比较两种方法的结果,发现在墙之间距离较大时,Rosenfeld密度泛函理论的结果与分子动力学模拟的结果符合很好;当两堵墙间的距离很小时,这两个结果之间存在明显的不一致.另外,还研究了约束条件下密度分布的结构.     

AB2型聚合物流体的表面结构性质

在密度泛函理论(DFT)框架下, 应用改进的基本度量理论(MFMT)表达硬球作用对自由能泛函的贡献, 根据统计力学理论结合加权密度近似(WDA)表达聚合作用对自由能泛函的贡献, 建立了描述AB2型聚合物流体的化学势, 得到了聚合物流体在硬球颗粒表面的密度分布表达式, 计算了聚合物流体在硬球颗粒表面附近的密度分布, 并探讨了体积分数、聚合程度和硬球颗粒尺度对体系密度分布的影响. 此外, 通过体系密度分布, 进一步分析了体积分数、聚合程度和硬球颗粒尺度与剩余吸附的关系.     

平行板间超支化聚合物流体的密度分布和溶剂化力

在密度泛函理论(DFT)框架下,应用改进的基本度量理论(MFMT)和统计力学理论结合加权密度近似(WDA)分别表示硬球作用和聚合作用对自由能泛函的贡献,建立了描述ABg型超支化聚合物流体的化学势,得到了聚合物流体在平行板间的密度分布表达式,计算了聚合物流体在两平行板间的密度分布,并探讨了体积分数、反应程度和单体中B类官能团数对体系密度分布的影响.此外,通过体系密度分布,进一步分析了反应程度、板间宽度与溶剂化力的关系.     

氢气在碳纳米管内吸附的密度泛函理论研究

本文应用密度泛函理论研究纯流体氢气在单壁碳纳米管内吸附过程,采用硬球状态方程改进的基本测量理论表征硬球的斥力作用,离子吸引项的贡献则用微扰理论描述.在温度为300k,氢气本体对比密度范围为0.2～0.7条件下,计算了三种不同尺寸的碳纳米管氢气吸附的密度分布,其密度泛函计算结果与计算机分子模拟数据完全一致.     

基于密度泛函理论研究二元排斥Yukawa流体的表面结构性质

基于自由能密度泛函理论(DFT)考察了二元排斥Yukawa (HCRY)流体在不同外场下的密度分布. 基于微扰理论, 体系的Helmholtz自由能泛函采用硬球排斥部分和长程色散部分贡献之和, 其中Kierlik和Rosinberg的加权密度近似(WDA)被用来计算硬球排斥部分, 而色散部分采用平均场理论(MFT)进行描述. 为了验证DFT计算结果的合理性, 研究中采用巨正则Monte Carlo(GCMC)模拟计算了在不同主体相密度、硬核直径和位能参数比的条件下二元HCRY混合流体的密度分布. 结果表明, 该DFT计算结果与GCMC模拟值吻合良好.     

一个新的桥泛函及其在非均一流体密度泛函理论中的应用

基于对OZ 方程的渐近行为与Taylor级数展开的分析，提出了一个新的桥泛函，桥泛函被表达为间接相关函数的函数，Taylor级数展开的重整化导致了一个可调参数，通过将所提出的桥泛函与一个最近提出的密度泛函理论方法学，以及单个硬墙的sum 规则结合，可以确定可调参数.所提出的桥泛函能预言如下非均一流体的密度分布：硬球流体接近一个硬墙与在球形空隙内，Lennard-Jones 流体与缔合硬球流体在两个硬墙之内.理论预言与文献所报导的模拟数据符合很好.     

受限于纳米狭缝中的四缔合Lennard-Lones流体的相平衡和界面张力研究

分别用改进的基础测量理论和平均球近似理论表达短程作用和长程作用对四缔合Lennard-Jones流体的过剩自由能的贡献. 在密度函泛理论的框架下, 研究了平均密度等温线, 密度分布, 未缔合分子在平衡汽相和液相中的分布, 相平衡以及平衡时的界面张力等热力学性质. 分析了缔合能量, 流体-固体作用和孔宽对受限于纳米狭缝中的四缔合Lennard-Jones流体相行为的影响.     

电场中B2分子特性研究

利用密度泛函理论(DFT)B3LYP/6-3l1＋G(2d)方法研究在不同方向电场(0～＋0.02 a.u.)作用下的B2分子的基态键长、总能量、偶极矩、最高占据轨道(HOMO)能量、最低空轨道(LUMO)能量、能隙及势能曲线的变化规律. 结果表明: 在一定外加电场范围内, 随电场强度的增大, 分子键长变大; 总能量降低; 偶极矩增大; HOMO能级、LUMO能级均降低; 能隙依赖外电场方向, 平行分子轴(Z)方向电场使能隙递减, 垂直分子轴(X)方向电场使能隙递增; 分子势能降低, 平行分子轴线(Z)方向电场对分子势能的影响随着核间距的增大而增大, 原有的“势能平台”遭到破坏.     

用密度泛函理论研究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流体在狭缝中的相平衡.     

Structure of spherical electric double layers: a density functional approach

A density functional theory is presented for the structure of spherical electric double layers within the restricted primitive model, where the macroion is considered as a hard sphere having uniform surface charge density, the small ions as charged hard spheres, and the solvent is taken as a dielectric continuum. The theory is partially perturbative as the hard-sphere contribution to the one-particle correlation function is evaluated using suitably averaged weighted density and the ionic part is obtained through a second-order functional Taylor expansion around the uniform fluid. The theory is in quantitative agreement with Monte Carlo simulation for the density profiles and the zeta potentials over a wide range of macroion sizes and electrolyte concentrations. The theory is able to provide interesting insights about the layering and the charge inversion phenomena occurring at the interface.     

A novel method to describe the interaction pressure between charged plates with application of the weighted correlation approach

Based on the Euler-Lagrange equation for ion density distribution in an inhomogeneous, charged, and hard-sphere fluid, a novel method is proposed to determine the interaction pressure between charged plates. The resulting expression is a sum of distinct physical contributions to the pressure, which involves different contributions to the single-particle direct correlation function. It can, therefore, be conveniently used in any density functional approach to facilitate analysis of the pressure components. In this study, the so-called fundamental measure theory (FMT)∕weighted correlation approach (WCA) approach is applied to estimate both the hard-sphere and the electric residual contributions to the single-particle direct correlation function, upon the calculation of the ionic density profiles between charged plates. The results, against the Monte Carlo simulations, show that the FMT∕WCA approach is superior to the typical FMT∕mean spherical approximation approach of the density functional theory in predicting the interaction pressure between charged plates immersed in an electrolyte solution upon various conditions in the primitive model. The FMT∕WCA approach can well capture the fine features of the pressure-separation dependence, to reproduce not only the shoulder shape and the weak attractions in monovalent electrolytes but also the strongly oscillatory behavior of pressure in divalent electrolytes where pronounced attractions are observed. In addition, it is found that the FMT∕WCA approach even has an advantage over the anisotropic, hyper-netted chain approach in that it agrees with the Monte Carlo results to a very good extent with, however, much less computational effort.     

Polarization relaxation in an ionic liquid confined between electrified walls

The response of a room temperature molten salt to an external electric field when it is confined to a nanoslit is studied by molecular dynamics simulations. The fluid is confined between two parallel and oppositely charged walls, emulating two electrified solid-liquid interfaces. Attention is focused on structural, electrostatic, and dynamical properties, which are compared with those of the nonpolarized fluid. It is found that the relaxation of the electrostatic potential, after switching the electric field off, occurs in two stages. A first, subpicosecond process accounts for 80% of the decay and is followed by a second subdiffusive process with a time constant of 8 ps. Diffusion is not involved in the relaxation, which is mostly driven by small anion translations. The relaxation of the polarization in the confined system is discussed in terms of the spectrum of charge density fluctuations in the bulk.     

Molecular dynamics study of screening at ionic surfaces

Molecular dynamics simulations of NaCl fluid are used to understand the behavior of ionic fluid to screen the field generated by charges on the ionic crystal surfaces in absence of any external electric field. The NaCl fluid in the strongly coupled regime (corresponding to the melt) in contact with the charged octopolar (111) NaCl surface shows that the spatial correlations decay in an oscillatory manner, with a screening length lambdaQ given by the envelope of the damped oscillations. By contrast to the Debye-Huckel theory, in the strongly coupled regime, lambdaQ increases with increasing coupling strength (also seen in bulk ionic simulations). The NaCl fluid confined between neutral (100) NaCl surfaces also shows weak oscillatory charge decay near the surface. Similar oscillatory exponential decay was seen when the NaCl fluid was confined between two analytically smooth neutral walls. The origin of these oscillations was due to the difference in ion sizes. NaCl fluid confined between neutral octopolar (110) and dipolar (110) surface show stronger density oscillations than (100) surface but comparatively very weak charge oscillations. This paper shows that the strength of the charges on the crystal surfaces is enough to induce a characteristic spatial distribution of charges in the contacting fluid and the extent of distribution depends on the type of surface.     

Does confining the hard-sphere fluid between hard walls change its average properties?

We use grand canonical transition-matrix Monte Carlo and discontinuous molecular dynamics simulations to generate precise thermodynamic and kinetic data for the equilibrium hard-sphere fluid confined between smooth hard walls. These simulations show that the pronounced inhomogeneous structuring of the fluid normal to the confining walls, often the primary focus of density functional theory studies, has a negligible effect on many of its average properties over a surprisingly broad range of conditions. We present one consequence of this insensitivity to confinement: a simple analytical equation relating the average density of the confined fluid to that of the bulk fluid with equal activity. Nontrivial implications of confinement for average fluid properties do emerge in this system, but only when the fluid is both (i) dense and (ii) confined to a gap smaller than approximately three particle diameters. For this limited set of conditions, we find that "in-phase" oscillatory deviations in excess entropy and self-diffusivity (relative to the behavior of the bulk fluid at the same average density) occur as a function of gap size. These paired thermodynamic/kinetic deviations from bulk behavior appear to reflect the geometric packing frustration that arises when the confined space cannot naturally accommodate an integer number of particle layers.     

Density functional theory for copolymers confined in a nanoslit

A density functional theory is developed for copolymers confined in a nanoslit on the basis of our previous work for homopolymers. The theory accurately captures the structural characteristics for diblock and alternating copolymers composed of hard-sphere or square-well segments. Satisfactory agreement is obtained between the theoretical predictions and simulation results in segment density profiles, segment fractions, and partition coefficients. Structures under confinement strongly depend on the substituent segment sizes for the hard-sphere copolymers and also on the segment-wall attractions for the square-well copolymers. Alternating copolymers are found to behave as homopolymers with effective segment size, and effective segment-segment and segment-wall interactions.     

Structures and correlation functions of multicomponent and polydisperse hard-sphere mixtures from a density functional theory

The structures of nonuniform binary hard-sphere mixtures and the correlation functions of uniform ternary hard-sphere mixtures were studied using a modified fundamental-measure theory based on the weight functions of Rosenfeld [Rosenfeld, Phys. Rev. Lett. 63, 980 (1989)] and Boublik-Mansoori-Carnahan-Starling-Leland equation of state [Boublik, J. Chem. Phys. 53, 471 (1970); Mansoori et al., J. Chem. Phys. 54, 1523 (1971)]. The theoretical predictions agreed very well with the molecular simulations for the overall density profiles, the local compositions, and the radial distribution functions of uniform as well as inhomogeneous hard-sphere mixtures. The density functional theory was further extended to represent the structure of a polydisperse hard-sphere fluid near a hard wall. Excellent agreement was also achieved between theory and Monte Carlo simulations. The density functional theory predicted oscillatory size segregations near a hard wall for a polydisperse hard-sphere fluid of a uniform size distribution.     

Density Functional Theory and the Capillary Evaporation of a Liquid in a Slit

Density functional theory is applied to a Lennard-Jones fluid near a single hard wall and in a slit formed by two walls. We use some simplified versions of the Weeks-Chandler-Andersen (WCA) and the Barker-Henderson (BH) theories. Only the most crude mean field version of the WCA theory, in which the hard-sphere correlation function is set equal to unity for all distances, seems useful. Use of the full WCA approximation is impractical because the effective hard-sphere diameter is density dependent. Generally, the best results are obtained using the BH macroscopic compressibility approximation. Our earlier study of "evaporation" of Lennard-Jones molecules in a slit is extended to other densities using the mean field theory. Copyright 2000 Academic Press.     

Electro‐osmotic flow in nanochannels with voltage‐controlled polyelectrolyte brushes: Dependence on grafting density and normal electric field

Molecular dynamics simulations were performed for electro‐osmotic flow (EOF) confined in a polyelectrolyte‐grafted nanochannel under variable grafting density and normal electric field. With decreasing the value of the normal electric field, the brush undergoes a collapse transition, and the ion distribution is changed significantly. The brush thickness increases on increasing the grafting density at positive and weak negative electric fields, whereas a reduced brush thickness is observed at strong negative electric field. Our results further reveal that the flow velocity is not only dependent on conformational transition of the brush but also related to the cation and anion distributions. At low grafting density, the EOF is almost completely quenched at high electric field strength due to strong surface friction between ions and walls. For the case of very dense grafting, the flow velocity is influenced weakly within the brush when varying the grafting density. Additionally, a bidirectional flow occurs at an intermediate electric field. The investigation on fluid flux indicates that the fluid flux is insensitive to the grafting density, when the normal electric field is removed. For nonzero normal electric fields, a significant change in the fluid flux is observed at low grafting densities. © 2012 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys, 2012     

Response of Bi-disperse Polyelectrolyte Brushes to External Electric Fields——A Numerical Self-consistent Field Theory Study

The self-consistent field theory has been employed to numerically study the response of bi-disperse flexible polyelectrolyte (PE) brushes grafted on an electrode to electric fields generated by opposite surface charges on the PE-grafted electrode and a second parallel electrode. The numerical study reveals that, under a positive external electric field, the shorter and negatively charged PE chains are more responsive than the longer PE chains in terms of the relative changes in their respective brush heights. Whereas under a negative external electric field, the opposite was observed. The total electric force on the grafted PE chains was calculated and it was found that, under a positive external electric field, the magnitude of the total electric force acting on one shorter PE chain is larger than that on one longer PE chain, or vice versa. The underlying mechanism was unraveled through analyzing the total electric field across the two oppositely charged electrodes.