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

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

Structure of inhomogeneous attractive and repulsive hard-core yukawa fluid: grand canonical Monte Carlo simulation and density functional theory study

A density functional theory is proposed for an inhomogeneous hard-core Yukawa (HCY) fluid based on Rosenfeld's perturbative method. The excess Helmholtz energy functional is derived from a modified fundamental measure theory for the hard-core repulsion and a quadratic functional Taylor expansion for the long-ranged attractive or repulsive interactions. To test the established theory, grand canonical ensemble Monte Carlo simulations are carried out to simulate the density profiles of attractive and repulsive HCY fluid near a wall. Comparison with the results from the Monte Carlo simulations shows that the present density functional theory gives accurate density profiles for both attractive and repulsive HCY fluid near a wall. Both the present theory and simulations suggest that there is depletion for attractive HCY fluid at low temperature, but no depletion is found for repulsive HCY fluid. The calculated results indicate that the present density functional theory is better than those of the modified version of the Lovett-Mou-Buff-Wertheim and other density functional theories. The present theory is simple in form and computationally efficient. It predicts accurate radial distribution functions of both attractive and repulsive HCY fluid except for the repulsive case at high density, where the theory overestimates the radial distribution function in the vicinity of contact.     

Density functional theory of homopolymer mixtures confined in a slit

A density functional theory (DFT) is developed for polymer mixtures with shorted-ranged attractive interparticle interactions confined in a slit. Different weighting functions are used separately for the repulsive part and the attractive part of the excess free energy functional by applying the weighted density approximation. The predicted results by DFT are in good agreement with the corresponding simulation data indicating the reliability of the theory. Furthermore, the center-of-mass profiles and the end-to-end distance distributions are obtained by the single chain simulation; the predictions also agree well with simulation data. The results reveal that both the attraction of the slit wall and the temperature has stronger effect on longer chains than on shorter ones because the intrasegment correlation of chains increases with increasing chain length.     

Structures and adsorption of binary hard-core Yukawa mixtures in a slitlike pore: grand canonical Monte Carlo simulation and density-functional study

The grand canonical ensemble Monte Carlo simulation and density-functional theory are applied to calculate the structures, local mole fractions, and adsorption isotherms of binary hard-core Yukawa mixtures in a slitlike pore as well as the radial distribution functions of bulk mixtures. The excess Helmholtz energy functional is a combination of the modified fundamental measure theory of Yu and Wu [J. Chem. Phys. 117, 10156 (2002)] for the hard-core contribution and a corrected mean-field theory for the attractive contribution. A comparison of the theoretical results with the results from the Monte Carlo simulations shows that the corrected theory improves the density profiles of binary hard-core Yukawa mixtures in the vicinity of contact over the original mean-field theory. Both the present corrected theory and the simulations suggest that depletion and desorption occur at low temperature, and the local segregation can be observed in most cases. For binary mixtures in the hard slitlike pore, the present corrected theory predicts more accurate surface excesses than the original one does, while in the case of the attractive pore, no improvement is found in the prediction of a surface excess of the smaller molecule.     

Model colloidal fluid with competing interactions: bulk and interfacial properties

Using a simple mean field density functional theory (DFT), the authors investigate the structure and phase behavior of a model colloidal fluid composed of particles interacting via a pair potential which has a hard core of diameter sigma, is attractive Yukawa at intermediate separations, and is repulsive Yukawa at large separations. The authors analyze the form of the asymptotic decay of the bulk fluid correlation functions, comparing results from DFT with those from the self-consistent Ornstein-Zernike approximation (SCOZA). In both theories the authors find rich crossover behavior, whereby the ultimate decay of correlation functions changes from monotonic to long wavelength damped oscillatory decay on crossing certain lines in the phase diagram or sometimes from oscillatory to oscillatory with a longer wavelength. For some choices of potential parameters the authors find, within the DFT, a lambda line at which the fluid becomes unstable with respect to periodic density fluctuations. SCOZA fails to yield solutions for state points near such a lambda line. The propensity towards clustering of particles, which is reflected by the presence of a long wavelength (>sigma) slowly decaying oscillatory pair correlation function, and a structure factor that exhibits a very sharp maximum at small but nonzero wave numbers, is enhanced in states near the lambda line. The authors present density profiles for the planar liquid-gas interface and for fluids adsorbed at a planar hard wall. The presence of a nearby lambda transition gives rise to pronounced long wavelength oscillations in the one-body density profiles at both types of interface.     

Density profiles and solvation forces for a Yukawa fluid in a slit pore

The effect of varying wall-particle and particle-particle interactions on the density profiles near a single wall and the solvation forces between two walls immersed in a fluid of particles is investigated by grand canonical Monte Carlo simulations. Attractive and repulsive particle-particle and particle-wall interactions are modeled by a versatile hard-core Yukawa form. These simulation results are compared to theoretical calculations using the hypernetted chain integral equation technique, as well as with fundamental measure density functional theory (DFT), where particle-particle interactions are either treated as a first order perturbation using the radial distribution function or else with a DFT based on the direct-correlation function. All three theoretical approaches reproduce the main trends fairly well, but exhibit inconsistent accuracy, particularly for attractive particle-particle interactions. We show that the wall-particle and particle-particle attractions can couple together to induce a nonlinear enhancement of the adsorption and a related "repulsion through attraction" effect for the effective wall-wall forces. We also investigate the phenomenon of bridging, where an attractive wall-particle interaction induces strongly attractive solvation forces.     

Hybrid density functional theory for homopolymer mixtures confined in a selective nanoslit

By integrating polymer density function theory (DFT) and single-chain molecular simulation, a hybrid DFT is developed for homopolymer mixtures confined in a selective nanoslit. Two weighting functions are adopted separately in the polymer DFT for repulsive and attractive contributions to the excess free energy functional. The theoretical results agree well with simulation data for the density profiles, configurations (tail, loop and train), adsorption amounts, layer thicknesses, and partition coefficients. The polymer-slit interaction is found to have a large effect on the density profiles and partition coefficients but is found to have a small effect on the average sizes and percentages of the configurations. Nearly half of the polymer segments form tails, and the other half form trains. In addition, bridges are observed to form for sufficiently long polymer chains. As the length difference between two polymers increases, the effect of chain connectivity becomes increasingly important.     

Comparison of random-walk density functional theory to simulation for bead-spring homopolymer melts

Density profiles for a homopolymer melt near a surface are calculated using a random-walk polymeric density functional theory, and compared to results from molecular dynamics simulations. All interactions are of a Lennard-Jones form, for both monomer-monomer interactions and surface-monomer interactions, rather than the hard core interactions which have been most investigated in the literature. For repulsive systems, the theory somewhat overpredicts the density oscillations near a surface. Nevertheless, near quantitative agreement with simulation can be obtained with an empirical scaling of the direct correlation function. Use of the random phase approximation to treat attractive interactions between polymer chains gives reasonable agreement with simulation of dense liquids near neutral and attractive surfaces.     

Predicting adsorption enthalpies on silicalite and HZSM‐5: A benchmark study on DFT strategies addressing dispersion interactions

We evaluated the accuracy of periodic density functional calculations for adsorption enthalpies of water, alkanes, and alcohols in silicalite and HZSM‐5 zeolites using a gradient‐corrected density functional with empirical dispersion corrections (PBE‐D) as well as a nonlocal correlation functional (vdW‐DF2). Results of both approaches agree in acceptable fashion with experimental adsorption energies of alcohols in silicalite, but the adsorption energies for n‐alkanes in both zeolite models are overestimated, by 21?46 kJ mol^?1. For PBE‐D calculations, the adsorption of alkanes is exclusively determined by the empirical dispersion term, while the generalized gradient approximation‐DFT part is purely repulsive, preventing the molecule to come too close to the zeolite walls. The vdW‐DF2 results are comparable to those of PBE‐D calculations, but the latter values are slightly closer to the experiment in most cases. Thus, both computational approaches are unable to reproduce available experimental adsorption energies of alkanes in silicalite and HZSM‐5 zeolite with chemical accuracy. © 2014 Wiley Periodicals, Inc.     

Mixing rules for van-der-Waals type equation of state based on thermodynamic perturbation theory

A concept based on the thermodynamic perturbation theory for a `simple fluid' has been applied to the attractive term of a van-der-Waals type equation of state (EOS) to derive a simple mixing rule for the a parameter. The new mixing rule is a small correction to the original one-fluid approximation to account for the influence of particles of j-type on the correlation function of ii-type in a mixture consisting of particles of i and j types. The importance of the correction has been shown by comparison of the calculated results for binary mixtures of Lennard–Jones fluids with the data obtained by numerical method (Monte-Carlo simulation). The new mixing rules can be considered as a flexible generalization of the conventional mixing rules and can be reduced to the original v-d-W mixing rules by defaulting the extra binary parameters to zero. In this way the binary parameters already available in the literature for many systems can be used without any additional regression work. Extension of the new mixing rules to a multicomponent system do not suffer from `Michelsen–Kistenmacher syndrome' and provide the correct limit for the composition dependence of second virial coefficients. Their applicability has been illustrated by various examples of vapor–liquid and liquid–liquid equilibria using a modified Patel–Teja EOS. The new mixing rules can be applied to any EOS of van-der-Waals type, i.e., EOS containing two terms which reflect the contributions of repulsive and attractive intermolecular forces.     

Further test of third order + second-order perturbation DFT approach: hard core repulsive Yukawa fluid subjected to diverse external fields

Grand canonical Monte Carlo simulation is used to investigate density profiles of hard-core repulsive Yukawa (HCRY) model fluid under the influence of various external fields and radial distribution function (RDF) of the bulk HCRY system. The aim of these extensive simulations is to provide exact data for purely repulsive interaction potential against which the validity of a third order + second-order perturbation DFT approach can be tested. It is found that a semiempirical parametrized bridge function due to Malijevsky and Labik performs very well for the RDF of the bulk HCRY fluid. Incorporation of a bulk second-order direct correlation function (DCF) of the HCRY fluid based on the Malijevsky-Labik bridge function into the third order + second-order perturbation DFT approach yields the resulting theoretical predictions for the density profiles of inhomogeneous HCRY fluid that are in a very good agreement with the simulation data, an exception being somewhat larger deviations appearing for the structure of the fluid around the center of a hard spherical cavity. Both theory and simulation predict layering transition and gas-liquid coexistence phenomena occurring with the HCRY model fluid under confined conditions. For the case of an inverse sixth-power repulsive potential under the influence of a flat stationary wall defined by an inverse twelfth-power repulsive potential, the present third order + second-order perturbation DFT approach is found to be superior to several existing weighted density approximations (WDA) and partitioned WDA.     

Structure and thermodynamics of hard-core Yukawa fluids: thermodynamic perturbation approaches

The thermodynamic perturbation theories, which are based on the power series of a coupling constant (λ-expansion), have been proposed for studying the structural and thermodynamic properties of a hard-core Yukawa (HCY) fluid: one (A1-approximation) is the perturbation theory based on the hard-sphere repulsion as a reference system. The other (A2-approximation) is the perturbation theory based on the reference system which incorporates both the repulsive and short-range attractive interactions. The first-order mean-spherical approximation (FMSA) provided by Tang and Lu [J. Chem. Phys. 99, 9828 (1993)] has been employed for investigating the thermodynamic properties of a HCY fluid using the alternative method via the direct correlation function. The calculated results show that (i) the A1 and A2 approximations are in excellent agreements with previous computer simulation results in the literature and compare with the semi-empirical works of Shukla including the higher-order free energy terms, (ii) the A1 and A2 approximations are better than the FMSA and the mean-spherical approximation, (iii) the A2-approximation compares with the A1-approximation, even though the perturbation effect of an A2-approximation is much smaller than that of an A1-approximation, and that (iv) the FMSA study is particularly of advantage in providing the structure and thermodynamics in a simple and analytic manner.     

Solvent mediated interactions close to fluid-fluid phase separation: microscopic treatment of bridging in a soft-core fluid

Using density functional theory we calculate the density profiles of a binary solvent adsorbed around a pair of big solute particles. All species interact via repulsive Gaussian potentials. The solvent exhibits fluid-fluid phase separation, and for thermodynamic states near to coexistence the big particles can be surrounded by a thick adsorbed "wetting" film of the coexisting solvent phase. On reducing the separation between the two big particles we find there can be a "bridging" transition as the wetting films join to form a fluid bridge. The effective (solvent mediated) potential between the two big particles becomes long ranged and strongly attractive in the bridged configuration. Within our mean-field treatment the bridging transition results in a discontinuity in the solvent mediated force. We demonstrate that accounting for the phenomenon of bridging requires the presence of a nonzero bridge function in the correlations between the solute particles when our model fluid is described within a full mixture theory based upon the Ornstein-Zernike equations.     

Density and surface tension of pure ionic liquid 1-butyl-3-methyl-imidazolium l-lactate and its binary mixture with alcohol and water

The density and surface tension of the pure ionic liquid 1-butyl-3-methyl-imidazolium l-lactate were measured from T (293.15 to 343.15) K. The coefficient of thermal expansion, molecular volume, standard entropy, lattice energy, surface entropy, surface enthalpy, and enthalpy of vaporization were calculated from the experimental values. Density and surface tension were also determined for binary mixtures of {1-butyl-3-methyl-imidazolium l-lactate + water/alcohol (methanol, ethanol, and 1-butanol)} systems over the whole composition range from T (298.15 to 318.15) K at atmospheric pressure. The partial molar volume, excess partial molar volume and apparent molar volume of the component IL and alcohol/water in the binary mixtures were discussed as well as limiting properties at infinite dilution and the thermal expansion coefficients of the four binary mixtures. The surface properties of the four binary mixtures were also discussed.     

Application of an improved simplified perturbed hard chain theory equation of state for pure compounds and binary asymmetric mixtures

The simplified perturbed hard chain theory equation of state is modified using recently developed repulsive and attractive terms. Both terms meet the low density and close-packed density boundary conditions and are in reasonable agreement with molecular simulation data. The modified equation of state accurately predicts pure component properties including saturated vapor volume, liquid density, vapor pressure and enthalpy of vaporization. Compared with the original equation, the modified one predicts physical properties more accurately. However, the improvement of predicting enthalpy of vaporization is marginal. The two equations are tested for predicting the phase behavior of binary asymmetric mixtures. It is shown that the difference in predicting the phase behavior is not appreciable.     

Spectroscopic Evidence for Clusters of Like‐Charged Ions in Ionic Liquids Stabilized by Cooperative Hydrogen Bonding

Direct spectroscopic evidence for hydrogen‐bonded clusters of like‐charged ions is reported for ionic liquids. The measured infrared O?H vibrational bands of the hydroxyethyl groups in the cations can be assigned to the dispersion‐corrected DFT calculated frequencies of linear and cyclic clusters. Compensating the like‐charge Coulomb repulsion, these cationic clusters can range up to cyclic tetramers resembling molecular clusters of water and alcohols. These ionic clusters are mainly present at low temperature and show strong cooperative effects in hydrogen bonding. DFT‐D3 calculations of the pure multiply charged clusters suggest that the attractive hydrogen bonds can compete with repulsive Coulomb forces.     

Viscosity,density, interfacial tension and compositional data for near critical mixtures of methane + butane and methane + decane systems at 310.95 K

Single- and two-phase binary mixtures of near critical methane + butane and methane + decane systems have been investigated at 310.95 K. A capillary tube viscometer was used to measure the viscosity and a high-pressure densitometer was employed for density measurements. The mixture was prepared gravimetrically and a direct sampling system was used to measure the composition of equilibrium phases in the two-phase region by gas chromatography. The meniscus height technique was used to generate the interfacial tension data.     

Structure and adsorption of a hard-core multi-Yukawa fluid confined in a slitlike pore: grand canonical Monte Carlo simulation and density functional study

Because of the increasing interest in studying the phenomenon exhibited by charge-stabilized colloidal suspensions in confining geometry, we present a density functional theory (DFT) for a hard-core multi-Yukawa fluid. The excess Helmholtz free-energy functional is constructed by using the modified fundamental measure theory and Rosenfeld's perturbative method, in which the bulk direct correlation function is obtained from the first-order mean spherical approximation. To validate the established theory, grand canonical ensemble Monte Carlo (GCMC) simulations are carried out to determine the density profiles and surface excesses of multi-Yukawa fluid in a slitlike pore. Comparisons of the theoretical results with the GCMC data suggest that the present DFT gives very accurate density profiles and surface excesses of multi-Yukawa fluid in the slitlike pore as well as the radial distribution functions of the bulk fluid. Both the DFT and the GCMC simulations predict the depletion of the multi-Yukawa fluid near a nonattractive wall, while the mean-field theory fails to describe this depletion in some cases. Because the simple form of the direct correlation function is used, the present DFT is computationally as efficient as the mean-field theory, but reproduces the simulation data much better than the mean-field theory.     

Large attractive depletion interactions in soft repulsive-sphere binary mixtures

We consider binary mixtures of soft repulsive spherical particles and calculate the depletion interaction between two big spheres mediated by the fluid of small spheres, using different theoretical and simulation methods. The validity of the theoretical approach, a virial expansion in terms of the density of the small spheres, is checked against simulation results. Attention is given to the approach toward the hard-sphere limit and to the effect of density and temperature on the strength of the depletion potential. Our results indicate, surprisingly, that even a modest degree of softness in the pair potential governing the direct interactions between the particles may lead to a significantly more attractive total effective potential for the big spheres than in the hard-sphere case. This might lead to significant differences in phase behavior, structure, and dynamics of a binary mixture of soft repulsive spheres. In particular, a perturbative scheme is applied to predict the phase diagram of an effective system of big spheres interacting via depletion forces for a size ratio of small and big spheres of 0.2; this diagram includes the usual fluid-solid transition but, in the soft-sphere case, the metastable fluid-fluid transition, which is probably absent in hard-sphere mixtures, is close to being stable with respect to direct fluid-solid coexistence. From these results, the interesting possibility arises that, for sufficiently soft repulsive particles, this phase transition could become stable. Possible implications for the phase behavior of real colloidal dispersions are discussed.