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.     

基于密度泛函理论研究二元排斥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.     

Surface tension of associating fluids by Monte Carlo simulations

Canonical Monte Carlo (NVT-MC) simulations were performed to obtain surface tension and coexistence densities at the liquid-vapor interface of one-site associating Lennard-Jones and hard-core Yukawa fluids, as functions of association strength and temperature. The method to obtain the components of the pressure tensor from NVT-MC simulations was validated by comparing the equation of state of the associative hard sphere system with that coming from isothermal-isobaric Monte Carlo simulations. Surface tension of the associative Lennard-Jones fluid determined from NVT-MC is compared with previously reported results obtained by molecular dynamics simulations of a pseudomixture model of monomers and dimers. A good agreement was found between both methods. Values of surface tension of associative hard-core Yukawa fluids are presented here for the first time.     

Phase diagram of a binary symmetric hard-core Yukawa mixture

We assess the accuracy of the self-consistent Ornstein-Zernike approximation for a binary symmetric hard-core Yukawa mixture by comparison with Monte Carlo simulations of the phase diagrams obtained for different choices of the ratio alpha of the unlike-to-like interactions. In particular, from the results obtained at alpha=0.75 we find evidence for a critical endpoint in contrast to recent studies based on integral equation and hierarchical reference theories. The variation of the phase diagrams with range of the Yukawa potential is investigated.     

Thermodynamic and structural properties of repulsive hard-core Yukawa fluid: integral equation theory, perturbation theory and Monte Carlo simulations

The thermodynamic and structural properties of purely repulsive hard-core Yukawa particles in the fluid state are determined through Monte Carlo simulation and modeled using perturbation theory and integral equation theory in the mean spherical approximation (MSA). Systems of particles with Yukawa screening lengths of 1.8, 3.0, and 5.0 are examined with results compared to variations of MSA and perturbation theory. Thermodynamic properties were predicted well by both theories in the fluid region up to the fluid-solid phase boundary. Further, we found that a simplified exponential version of the MSA is the most accurate at predicting radial distribution function at contact. Radial distribution function of repulsive hard-core Yukawa particles are also reported. The results show that methods based on MSA and perturbation theory that are typically applied to the attractive hard-core Yukawa potential can also be extended to the purely repulsive hard-core Yukawa potential.     

Storage of hydrogen at 303 K in graphite slitlike pores from grand canonical Monte Carlo simulation

Grand canonical Monte Carlo (GCMC) simulations were used for the modeling of the hydrogen adsorption in idealized graphite slitlike pores. In all simulations, quantum effects were included through the Feynman and Hibbs second-order effective potential. The simulated surface excess isotherms of hydrogen were used for the determination of the total hydrogen storage, density of hydrogen in graphite slitlike pores, distribution of pore sizes and volumes, enthalpy of adsorption per mole, total surface area, total pore volume, and average pore size of pitch-based activated carbon fibers. Combining experimental results with simulations reveals that the density of hydrogen in graphite slitlike pores at 303 K does not exceed 0.014 g/cm(3), that is, 21% of the liquid-hydrogen density at the triple point. The optimal pore size for the storage of hydrogen at 303 K in the considered pore geometry depends on the pressure of storage. For lower storage pressures, p < 30MPa, the optimal pore width is equal to a 2.2 collision diameter of hydrogen (i.e., 0.65 nm), whereas, for p congruent with 50MPa, the pore width is equal to an approximately 7.2 collision diameter of hydrogen (i.e., 2.13 nm). For the wider pores, that is, the pore width exceeds a 7.2 collision diameter of hydrogen, the surface excess of hydrogen adsorption is constant. The importance of quantum effects is recognized in narrow graphite slitlike pores in the whole range of the hydrogen pressure as well as in wider ones at high pressures of bulk hydrogen. The enthalpies of adsorption per mole for the considered carbonaceous materials are practically constant with hydrogen loading and vary within the narrow range q(st) congruent with 7.28-7.85 kJ/mol. Our systematic study of hydrogen adsorption at 303 K in graphite slitlike pores gives deep insight into the timely problem of hydrogen storage as the most promising source of clean energy. The calculated maximum storage of hydrogen is equal to approximately 1.4 wt %, which is far from the United States Department of Energy (DOE) target (i.e., 6.5 wt %), thus concluding that the total storage amount of hydrogen obtained at 303 K in graphite slitlike pores of carbon fibers is not sufficient yet.     

Density functional theory study on the structure and capillary phase transition of a polymer melt in a slitlike pore: effect of attraction

A density functional theory is proposed to investigate the effects of polymer monomer-monomer and monomer-wall attractions on the density profile, chain configuration, and equilibrium capillary phase transition of a freely jointed multi-Yukawa fluid confined in a slitlike pore. The excess Helmholtz energy functional is constructed by using the modified fundamental measure theory, Wertheim's first-order thermodynamic perturbation theory, and Rosenfeld's perturbative method, in which the bulk radial distribution function and direct correlation function of hard-core multi-Yukawa monomers are obtained from the first-order mean spherical approximation. Comparisons of density profiles and bond orientation correlation functions of inhomogeneous chain fluids predicted from the present theory with the simulation data show that the present theory is very accurate, superior to the previous theory. The present theory predicts that the polymer monomer-monomer attraction lowers the strength of oscillations for density profiles and bond orientation correlation functions and makes the excess adsorption more negative. It is interesting to find that the equilibrium capillary phase transition of the polymeric fluid in the hard slitlike pore occurs at a higher chemical potential than in bulk condition, but as the attraction of the pore wall is increased sufficiently, the chemical potential for equilibrium capillary phase transition becomes lower than that for bulk vapor-liquid equilibrium.     

Excluded volume effects in macromolecular forces and ion-interface interactions

A charged Yukawa liquid confined in a slit nanopore is studied in order to understand excluded volume effects in the interaction force between the pore walls. A previously developed self-consistent scheme [S. Buyukdagli, C. V. Achim, and T. Ala-Nissila, J. Stat. Mech. 2011, P05033] and a new simpler variational procedure that self-consistently couple image forces, surface charge induced electric field, and pore modified core interactions are used to this aim. For neutral pores, it is shown that with increasing pore size, the theory predicts a transition of the interplate pressure from an attractive to a strongly repulsive regime associated with an ionic packing state, an effect observed in previous Monte Carlo simulations for hard core charges. We also establish the mean-field theory of the model and show that for dielectrically homogeneous pores, the mean-field regime of the interaction between the walls corresponds to large pores of size d > 4 ?. The role of the range of core interactions in the ionic rejection and interplate pressure is thoroughly analyzed. We show that the physics of the system can be split into two screening regimes. The ionic packing effect takes place in the regime of moderately screened core interactions characterized with the bare screening parameter of the Yukawa potential b ? 3/l(B), where l(B) is the Bjerrum length. In the second regime of strongly screened core interactions b ? 3/l(B), solvation forces associated with these interactions positively contribute to the ionic rejection driven by electrostatic forces and enhance the magnitude of the attractive pressure. For weakly charged pores without a dielectric discontinuity, core interactions make a net repulsive contribution to the interplate force and also result in oscillatory pressure curves, whereas for intermediate surface charges, these interactions exclusively strengthen the external pressure, thereby reducing the magnitude of the net repulsive interplate force. The pronounced dependence of the interplate pressure and ionic partition coefficients on the magnitude and the range of core interactions indicates excluded volume effects as an important ion specificity and a non-negligible ingredient for the stability of macromolecules in electrolyte solutions.     

Adsorption of short heteropolymers in slitlike pores

Adsorption of short linear heteropolymers in slitlike pores is studied using the density functional theory and Monte Carlo simulations. The molecules are assumed to be freely jointed tangent hard spheres. The segments have different affinity with regard to the walls. Each molecule contains one surface-binding segment that interacts with the walls via Lennard-Jones (3,9) potential and a number of segments interacting with surfaces via the hard-wall potentials. A position of the surface-binding segment in the chain can be arbitrarily chosen. We have studied the influence of the pore width, the chain length and the chemical structure of molecules on adsorption and the microscopic structure of the confined fluid. The theoretical predictions are compared with Monte Carlo simulations carried out for different 'isomeric' pentamers.     

Gas adsorption in active carbons and the slit-pore model 2: Mixture adsorption prediction with DFT and IAST

We use a fast density functional theory (a "slab-DFT") and the polydisperse independent ideal slit-pore model to predict gas mixture adsorption in active carbons. The DFT is parametrized by fitting to pure gas isotherms generated by Monte Carlo simulation of adsorption in model graphitic slit-pores. Accurate gas molecular models are used in our Monte Carlo simulations with gas-surface interactions calibrated to a high surface area carbon, rather than a low surface area carbon as in all previous work of this type, as described in part 1 of this work. We predict the adsorption of binary mixtures of carbon dioxide, methane, and nitrogen on two active carbons up to about 30 bar at near-ambient temperatures. We compare two sets of results; one set obtained using only the pure carbon dioxide adsorption isotherm as input to our pore characterization process, and the other obtained using both pure gas isotherms as input. We also compare these results with ideal adsorbed solution theory (IAST). We find that our methods are at least as accurate as IAST for these relatively simple gas mixtures and have the advantage of much greater versatility. We expect similar results for other active carbons and further performance gains for less ideal mixtures.     

On the corresponding states law of the Yukawa fluid

We have analyzed the currently available simulation results as well as performed some additional Monte Carlo simulation for the hard-core attractive Yukawa fluid in order to study its corresponding state behavior. We show that the values of reduced surface tension map onto the master curve and a universal equation of state can be obtained in the wide range of the attractive Yukawa tail length after a certain rescaling of the number density. Some comparisons with other nonconformal potentials are presented and discussed.     

Density functional theory for Yukawa fluids

We develop an approximate field theory for particles interacting with a generalized Yukawa potential. This theory improves and extends a previous splitting field theory, originally developed for counterions around a fixed charge distribution. The resulting theory bridges between the second virial approximation, which is accurate at low particle densities, and the mean-field approximation, accurate at high densities. We apply this theory to charged, screened ions in bulk solution, modeled to interact with a Yukawa potential; the theory is able to accurately reproduce the thermodynamic properties of the system over a broad range of conditions. The theory is also applied to "dressed counterions," interacting with a screened electrostatic potential, contained between charged plates. It is found to work well from the weak coupling to the strong coupling limits. The theory is able to reproduce the counterion profiles and force curves for closed and open systems obtained from Monte Carlo simulations.     

Phase diagram and surface tension of the hard-core attractive Yukawa model of variable range: Monte Carlo simulations

The liquid-vapor phase diagram and surface tension for hard-core Yukawa potential with 4相似文献   

Polymorph selection in the crystallization of hard-core Yukawa system

Colloid-colloid interactions in charge-stabilized dispersions can to some extent be represented by the hard-core Yukawa model. The crystallization process and polymorph selection of hard-core Yukawa model are studied by means of smart Monte Carlo simulations in the region of face-centered-cubic (fcc) phase. The contact value of hard-core Yukawa potential and the volume fraction of the colloids are fixed, while the Debye screening length can be varied. In the early stage of the crystallization, the precursors with relatively ordered liquid structure have been observed. Although the crystal structure of thermodynamically stable phase is fcc, the system crystallizes into a mixture of fcc and hexagonal close-packed (hcp) structures under small Debye screening length since the colloidal particles act as effective hard spheres. In the intermediate range of Debye screening length, the system crystallizes into a mixture of fcc, hcp, and body-centered-cubic (bcc). The existence of metastable hcp and bcc structures can be interpreted as a manifestation of the Ostwald’s step rule. Until the Debye screening length is large enough, the crystal structure obtained is almost a complete fcc suggesting the system eventually reaches to a thermodynamically stable state.     

Determination of surface tension in binary mixtures using transition-matrix Monte Carlo

We present a methodology based on grand-canonical transition-matrix Monte Carlo and finite-size scaling analysis to calculate surface tensions in binary mixtures. In particular, mixture transition-matrix Monte Carlo is first used to calculate apparent, system-size-dependent free-energy barriers separating coexisting fluid phases. Finite-size scaling is then used to extrapolate these values to the infinitely large system limit to determine the true thermodynamic surface tension. A key distinction of the methodology is that it yields the entire isothermal surface-tension curve for a binary mixture in a relatively small number of simulations. We demonstrate the utility of the method by calculating surface-tension curves for three binary Lennard-Jones mixtures. While we have only examined the surface tension of simple fluids in this work, the method is general and can be extended to molecular fluids as well as to determine interfacial tensions of liquid-liquid interfaces.     

Thermodynamic and structural properties of mixed colloids represented by a hard-core two-Yukawa mixture model fluid: Monte Carlo simulations and an analytical theory

The interaction between colloidal particles is well represented by a hard-core two-Yukawa potential. In order to assess the accuracy of theoretical predictions for the thermodynamic and structural properties of mixed colloids, standard Monte Carlo simulations are carried out for the hard-core two-Yukawa mixtures. In the simulations, one range parameter in the two-Yukawa potential is taken as 1.8 or 2.8647, and another is taken as 4, 8, or 13.5485. Both attractive and repulsive dominant cases of the potential outside the hard core are considered. The effects of temperature, density, composition, size and energy parameter ratios on internal energy, compressibility factor, and radial distribution function are investigated extensively. Theoretical calculations are performed in the framework of analytical solution for the Ornstein-Zernike equation with the first-order mean spherical approximation (FMSA). Our analysis shows that the FMSA is very accurate for the prediction of the compressibility factor of the hard-core two-Yukawa mixtures at all conditions studied. The FMSA generally predicts accurate internal energy, but overestimates the internal energy of the systems at lower temperatures. Furthermore, we found that a simplified exponential version of the FMSA predicts fairly good radial distribution function at contact for the mixed two-Yukawa fluids. The comparison of the theoretical compressibility factor with that from the Monte Carlo simulations suggests that the FMSA can be used to investigate the fluid-fluid equilibria of hard-core two-Yukawa mixtures.     

Effect of quenched size polydispersity on the fluid-solid transition in charged colloidal suspensions

We study the effect of quenched size polydispersity on the phase behavior of charged colloidal suspensions using free-energy calculations in Monte Carlo simulations. The colloids are assumed to interact with a hard-core repulsive Yukawa (screened-Coulomb) interaction with constant surface potential, so that the particles are polydisperse both in size and charge. In addition, we take the size distribution to be fixed in both the fluid and crystal phase (no size fractionation is allowed). We study the fluid-solid transition for various screening lengths and surface potentials, finding that upon increasing the size polydispersity the freezing transition shifts toward higher packing fractions and the density discontinuity between the two coexisting phases diminishes. Our results provide support for a terminal polydispersity above which the freezing transition disappears.     

Phase separation of the Widom-Rowlinson mixture in confined geometry

The behavior of binary Widom-Rowlinson [J. Chem. Phys. 52, 1670 (1970)] mixture in confined geometry is investigated. The influence of confinement on phase separation is examined. Coexistence curves for the mixture in slitlike pores and pores of complex geometry are calculated in Monte Carlo simulations. Finite size scaling analysis is used to determine precisely the location of critical point for the bulk mixture.