Constant-NTμ simulations: Free energy difference method for excess adsorption

We describe a new method to simulate confined fluids in equilibrium with a bulk. The equilibrium is first located at a low density, where conventional methods relying upon real or virtual particle exchanges are reliable. Thereafter, the chemical potential of the fluid is increased by the same amount in both systems using a variation of a recently developed free energy difference method in an isobaric (isotension) ensemble. The method, illustrated here for a simple fluid in a planar slit, is shown to be reliable up to high densities. As the method does not rely upon particle exchanges it is ideally suited to the simulation of equilibria in complex fluids, e.g., molecular liquids and polymers. © 1996 by John Wiley & Sons, Inc.     

Simulation of chemical potentials and phase equilibria in two- and three-dimensional square-well fluids: finite size effects

We study the simulation cell size dependence of chemical potential isotherms in subcritical square-well fluids by means of series of canonical ensemble Monte Carlo simulations with increasing numbers of particles, for both three-dimensional bulk systems and two-dimensional planar layers, using Widom-like particle insertion methods. By estimating the corresponding vapor/liquid coexistence densities using a Maxwell-like equal area rule for the subcritical chemical potential isotherms, we are able to study the influence of system size not only on chemical potentials but also on the coexistence properties. The chemical potential versus density isotherms show van der Waals-like loops in the subcritical vapor/liquid coexistence range that exhibit distinct finite size effects for both two- and three-dimensional fluids. Generally, in agreement with recent findings for related studies of Lennard-Jones fluids, the loops shrink with increasing number of particles. In contrast to the subcritical isotherms themselves, the equilibrium vapor/liquid densities show only a weak system size dependence and agree quantitatively with the best-known literature values for three-dimensional fluids. This allows our approach to be used to accurately predict the phase coexistence properties. Our resulting phase equilibrium results for two-dimensional square-well fluids are new. Knowledge concerning finite size effects of square-well systems is important not only for the simulation of thermodynamic properties of simple fluids, but also for the simulation of models of more complex fluids (such as aqueous or polymer fluids) involving square-well interactions.     

石墨烯/聚合物纳米复合电流变材料

电流变智能流体在外电场刺激下能快速可逆地改变自身流变性能,具有重要技术应用价值.传统的基于微米颗粒的电流变流体易于沉降并且电致屈服强度不高限制了技术应用,最近基于纳米颗粒的非传统电流变材料研究受到重视,特别是具有各向异性形貌的纳米颗粒悬浮液被发现具有明显增强的电/磁流变效应.本文介绍了最近基于石墨烯的二维纳米复合电流变材料的研究进展,主要包括石墨烯/半导聚合物、石墨烯/极性聚合物、石墨烯/碳等几种典型的电流变材料的制备、结构和电流变行为.研究表明利用石墨烯独特的二维纳米结构、优异的电学和热学性质可能为制备新颖的高性能纳米电流变材料提供途径。     

Coarse-grained molecular dynamics modeling of associating fluids: thermodynamics, liquid structure, and dynamics in the limit of zero association strength

A continuous coarse-grained potential model for associating fluids, consisting of an off-center specific site bonded with a harmonic potential to a center particle, has been developed and used in canonical molecular dynamics simulations. The thermodynamic, structural, and dynamic properties of the limiting nonassociating reference coarse-grained fluid are investigated as functions of the mass distribution and bond strength between center and site particles. It is theoretically shown and confirmed by simulation that in this limit variations in these potential parameters do not alter the equation of state of the reference coarse-grained fluid but that they have profound influences on both the translational and the rotational dynamics. From the simulation results we arrive at some guidelines that should be kept in mind in the selection of appropriate values for the model parameters. This work provides the precursory knowledge for the study of coarse-grained associating fluids using the conventional molecular dynamics method.     

A consistent calculation of the chemical potential for dense simple fluids

A general method to calculate the excess chemical potential betamuex, that is based on the Kirkwood coupling parameter's dependence of the correlation functions, is presented. The expression for the one particle bridge function B(1)r is derived for simple fluids with spherical interactions. Only the knowledge of the bridge function B(2)r is required. The accuracy of our approach is illustrated for a dense hard sphere fluid. As far as B(2)r is considered as exact, B(1)r is found to be, at high densities, the normalized bridge function -B(2)rB(2)(r=0). This expression ensures a consistent calculation of the excess chemical potential by satisfying implicitly the Gibbs-Duhem constraint. Only the pressure-consistency condition is necessary to calculate the structural and thermodynamic properties of the fluid.     

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

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

Optimized theory for simple and molecular fluids

An optimized closure approximation for both simple and molecular fluids is presented. A smooth interpolation between Perkus-Yevick and hypernetted chain closures is optimized by minimizing the free energy self-consistently with respect to the interpolation parameter(s). The molecular version is derived from a refinement of the method for simple fluids. In doing so, a method is proposed which appropriately couples an optimized closure with the variant of the diagrammatically proper integral equation recently introduced by this laboratory [K. M. Dyer et al., J. Chem. Phys. 123, 204512 (2005)]. The simplicity of the expressions involved in this proposed theory has allowed the authors to obtain an analytic expression for the approximate excess chemical potential. This is shown to be an efficient tool to estimate, from first principles, the numerical value of the interpolation parameters defining the aforementioned closure. As a preliminary test, representative models for simple fluids and homonuclear diatomic Lennard-Jones fluids were analyzed, obtaining site-site correlation functions in excellent agreement with simulation data.     

Chemical potential perturbation: a method to predict chemical potentials in periodic molecular simulations

A new method, called chemical potential perturbation (CPP), has been developed to predict the chemical potential as a function of density in periodic molecular simulations. The CPP method applies a spatially varying external force field to the simulation, causing the density to depend upon position in the simulation cell. Following equilibration the homogeneous (uniform or bulk) chemical potential as a function of density can be determined relative to some reference state after correcting for the effects of the inhomogeneity of the system. We compare three different methods of approximating this correction. The first method uses the van der Waals density gradient theory to approximate the inhomogeneous Helmholtz free energy density. The second method uses the local pressure tensor to approximate the homogeneous pressure. The third method uses the Triezenberg-Zwanzig definition of surface tension to approximate the inhomogeneous free energy density. If desired, the homogeneous pressure and Helmholtz free energy can also be predicted by the new method, as well as binodal and spinodal densities of a two-phase fluid region. The CPP method is tested using a Lennard-Jones (LJ) fluid at vapor, liquid, two-phase, and supercritical conditions. Satisfactory agreement is found between the CPP method and an LJ equation of state. The efficiency of the CPP method is compared to that for Widom's method under the tested conditions. In particular, the new method works well for dense fluids where Widom's method starts to fail.     

Communication: Generalizing Rosenfeld's excess-entropy scaling to predict long-time diffusivity in dense fluids of Brownian particles: from hard to ultrasoft interactions

Computer simulations are used to test whether a recently introduced generalization of Rosenfeld's excess-entropy scaling method for estimating transport coefficients in systems obeying molecular dynamics can be extended to predict long-time diffusivities in fluids of particles undergoing Brownian dynamics in the absence of interparticle hydrodynamic forces. Model fluids with inverse-power-law, Gaussian-core, and Hertzian pair interactions are considered. Within the generalized Rosenfeld scaling method, long-time diffusivities of ultrasoft Gaussian-core and Hertzian particle fluids, which display anomalous trends with increasing density, are predicted (to within 20%) based on knowledge of interparticle interactions, excess entropy, and scaling behavior of simpler inverse-power-law fluids.     

Anisotropic dynamics of dipolar liquids in narrow slit pores

We report molecular dynamics simulation results for Stockmayer fluids confined to narrow slitlike pores with structureless, nonconducting walls. The translational and rotational dynamics of the dipolar particles have been investigated by calculating autocorrelation functions, diffusion coefficients, and relaxation times for various pore widths (five or less particle diameters) and directions parallel and perpendicular to the walls. The dynamic properties of the confined systems are compared to bulk properties, where corresponding bulk and pore states at the same temperature and chemical potential are determined in parallel grand canonical Monte Carlo simulations. We find that the dynamic behavior inside the pore depends on the distance from the walls and can be strongly anisotropic even in globally isotropic systems. This concerns especially the particles in the surface layers close to the walls, where the single particle and collective dipolar relaxation resemble that of true two-dimensional dipolar fluids with different in-plane and out-of-plane relaxations. On the other hand, bulklike relaxation is observed in the pore center of sufficiently wide pores.     

气体混合物各组分化学势的Monte Carlo模拟

在则系综测试粒子Ｍｏｎｔｅ Ｃａｒｌｏ （ＧＣＭＣ）方法模拟常温下空气（以氮气为代表）及其污染物微量有机物（以苯为例）的混合物中各组分的化学势。模拟中,氮气和苯分子采用ＬＪ球型分子势能模型,采用Ｍｅｔｒｏｐｏｌｉｓ抽样及周期边界条件。通过模拟并拟合得到了３００．２Ｋ、苯的摩尔分数为０．０６２５,氮气及苯化学势与压力的关联式,以用于狭缝碳孔中该混合物体系的选择性吸附。     

Test of athermal terms of activity coefficient models by Monte Carlo simulation with hard-core models

The infinite dilution activity coefficients of exactly athermal fluids were calculated by Monte Carlo simulation with hard-core models. The hard-core models used in this work were hard-sphere and hard-spherocylinder models. The Widom test particle method was adopted to calculate the residual chemical potentials of solutes in pure solvent and in pure solute solutions. The infinite dilution activity coefficients of solutes were obtained from the residual chemical potentials of solutes. The infinite dilution activity coefficients calculated by Monte Carlo simulation were compared with those of athermal terms in activity coefficient equations. Staverman–Guggenheim equation overestimates the activity coefficients. The deviations of activity coefficients increase with increasing the hard-core volume of solute. Flory–Huggins equation based on molar volume gives good results for the hard-spherocylinder systems. Elbro-FV equation gives good results for both the hard-sphere and hard-spherocylinder systems.     

Equilibrium properties of dense hydrogen isotope gases based on the theory of simple fluids

We present a new method for the prediction of the equilibrium properties of dense gases containing hydrogen isotopes. The proposed approach combines the Feynman-Hibbs effective potential method and a deconvolution scheme introduced by Weeks et al. The resulting equations of state and the chemical potentials as functions of pressure for each of the hydrogen isotope gases depend on a single set of Lennard-Jones parameters. In addition to its simplicity, the proposed method with optimized Lennard-Jones potential parameters accurately describes the equilibrium properties of hydrogen isotope fluids in the regime of moderate temperatures and pressures. The present approach should find applications in the nonlocal density functional theory of inhomogeneous quantum fluids and should also be of particular relevance to hydrogen (clean energy) storage and to the separation of quantum isotopes by novel nanomaterials.     

Calculation Methods of Solution Chemical Potential and Application in Emulsion Microencapsulation

Several new biased sampling methods were summarized for solution chemical potential calculation methods in the field of emulsion microencapsulation. The principles, features, and calculation efficiencies of various biased Widom insertion sampling methods were introduced, including volume detection bias, simulation ensemble bias, and particle insertion bias. The proper matches between various types of solution in emulsion and biased Widom methods were suggested, following detailed analyses on the biased insertion techniques. The volume detection bias methods effectively improved the accuracy of the data and the calculation efficiency by inserting detection particles and were suggested to be used for the calculation of solvent chemical potential for the homogeneous aqueous phase of the emulsion. The chemical potential of water, argon, and fluorobenzene (a typical solvent of the oil phase in double emulsion) was calculated by a new, optimized volume detection bias proposed by this work. The recently developed Well-Tempered(WT)-Metadynamics method skillfully constructed low-density regions for particle insertion and dynamically adjusted the system configuration according to the potential energy around the detection point, and hence, could be used for the oil-polymer mixtures of microencapsulation emulsion. For the macromolecule solutes in the oil or aqueous phase of the emulsion, the particle insertion bias could be applied to greatly increase the success rate of Widom insertions. Readers were expected to choose appropriate biased Widom methods to carry out their calculations on chemical potential, fugacity, and solubility of solutions based on the system molecular properties, inspired by this paper.     

Synthesis and properties of magnetic fluids produced on the basis of magnetite particles

Magnetic fluids based on magnetite synthesized by the chemical condensation method at temperatures of 25, 40, 60, and 80°C were obtained and studied. Magnetite particles were examined by X-ray phase and X-ray fluorescence analyses and electron microscopy. The average size of the coherent scattering region of magnetite particles was 13–17 nm, depending on the synthesis temperature. Magnetic fluids were synthesized from magnetite particles obtained at 25 and 80°C, with water and octane serving as carrier fluids. The NMR method was used to determine the saturation magnetization and average magnetic moment of the particles: for water-based magnetic fluids, 2100 A m^–1 and 5.7 × 10^–19 A m² at magnetite particle synthesis temperature of 25°C and 3670 A m^–1 and 4.6 × 10^–19 A m² at magnetite particle synthesis temperature of 80°C; for octane-based magnetic fluids, 2250 A m^–1 and 4.1 × 10^–19 A m² at magnetite particle synthesis temperature of 25°C.     

Toward efficient chemical potential calculations by expanded ensemble simulations; to make the free energy pathway fairly level

A scheme is suggested of how to construct good bias potentials ("balancing factors") to be used in expanded ensemble (EE) calculations of chemical potentials of solutions. A combination of two strategies are used: (i) to use a pathway for particle insertions that avoids large variations in free energy and (ii) to use calculated free energy derivatives to construct a bias potential that makes the pathway fairly level. Only a few very short simulations are needed to accomplish the latter, and then, a full EE simulation is done to obtain the chemical potential. By practical calculations of the chemical potential of benzene, cyclohexane, and benzylamine in water, it is shown that this method is at least equally efficient to the recent adaptive EE (AEE) method by Aberg et al. (J. Chem. Phys. 2004, 120, 3370). Furthermore, the new method provides an alternative strategy that complements existing EE methods.     

Applications of Wang-Landau sampling to determine phase equilibria in complex fluids

Applications of the Wang-Landau algorithm for simulating phase coexistence at fixed temperature are presented. The number density is sampled using either volume scaling or particle insertion/deletion. The resulting algorithms, while being conceptually easy, are of comparable efficiency to existing multicanonical methods but with the advantage that neither the chemical potential nor the pressure at phase coexistence has to be estimated in advance of the simulation. First, we benchmark the algorithm against literature results for the vapor-liquid transition in the Lennard-Jones fluid. We then demonstrate the general applicability of the algorithm by studying vapor-liquid coexistence in two examples of complex fluids: charged soft spheres, which exhibit a transition similar to that in the restricted primitive model of ionic fluids, being characterized by strong ion pairing in the vapor phase; and Stockmayer fluids with high dipole strengths, in which the constituent particles aggregate to form chains, and for which the very existence of a transition has been widely debated. Finally, we show that the algorithm can be used to locate a weak isotropic-nematic transition in a fluid of Gay-Berne mesogens.     

Investigation of viscoelastic focusing of particles and cells in a zigzag microchannel

Microfluidic particle focusing has been a vital prerequisite step in sample preparation for downstream particle separation, counting, detection, or analysis, and has attracted broad applications in biomedical and chemical areas. Besides all the active and passive focusing methods in Newtonian fluids, particle focusing in viscoelastic fluids has been attracting increasing interest because of its advantages induced by intrinsic fluid property. However, to achieve a well-defined focusing position, there is a need to extend channel lengths when focusing micrometer-sized or sub-microsized particles, which would result in the size increase of the microfluidic devices. This work investigated the sheathless viscoelastic focusing of particles and cells in a zigzag microfluidic channel. Benefit from the zigzag structure of the channel, the channel length and the footprint of the device can be reduced without sacrificing the focusing performance. In this work, the viscoelastic focusing, including the focusing of 10 μm polystyrene particles, 5 μm polystyrene particles, 5 μm magnetic particles, white blood cells (WBCs), red blood cells (RBCs), and cancer cells, were all demonstrated. Moreover, magnetophoretic separation of magnetic and nonmagnetic particles after viscoelastic pre-focusing was shown. This focusing technique has the potential to be used in a range of biomedical applications.     

Molecular simulations of aqueous electrolyte solubility: 1. The expanded-ensemble osmotic molecular dynamics method for the solution phase

We have developed a molecular-level simulation technique called the expanded-ensemble osmotic molecular dynamics (EEOMD) method, for studying electrolyte solution systems. The EEOMD method performs simulations at a fixed number of solvent molecules, pressure, temperature, and overall electrolyte chemical potential. The method combines elements of constant pressure-constant temperature molecular dynamics and expanded-ensemble grand canonical Monte Carlo. The simulated electrolyte solution systems contain, in addition to solvent molecules, full and fractional ions and undissociated electrolyte molecular units. The fractional particles are coupled to the system via a coupling parameter that varies between 0 (no interaction between the fractional particle and the other particles in the system) and 1 (full interaction between the fractional particle and the other particles in the system). The time evolution of the system is governed by the constant pressure-constant temperature equations of motion and accompanied by random changes in the coupling parameter. The coupling-parameter changes are accepted with a probability derived from the expanded-ensemble osmotic partition function corresponding to the prescribed electrolyte chemical potential. The coupling-parameter changes mimic insertion/deletion of particles as in a crude grand canonical Monte Carlo simulation; if the coupling parameter becomes 0, the fractional particles disappear from the system, and as the coupling parameter reaches unity, the fractional particles become full particles. The method is demonstrated for a model of NaCl in water at ambient conditions. To test our approach, we first determine the chemical potential of NaCl in water by the thermodynamic integration technique and by the expanded-ensemble method. Then, we carry out EEOMD simulations for different specified values of the overall NaCl chemical potential and measure the concentration of ions resulting from the simulations. Both computations give consistent results, validating the EEOMD methodology.     

纳米Fe3O4磁性液体稳定性的研究

The effect of interaction（attraction effect and repulsion effect）between magnetic particles on the stability of nanometer magnetic fluids was analyzed theoretically. It was pointed that the size and concentration of magnetic particles and the coating of surfactant were important factors affecting the stability of magnetic fluids. Nanosized Fe3O4 magnetic fluids were prepared by the chemical co-precipitation method. Characterized by means of TEM,zeta potential and so on,the magnetic fluids had the mean size of 10 nm and exellent stability. The stability of magnetic fluids was researched experimentally as to the factors such as the method of adding starting materials,the speed of stirring,dispersion effect,surfactants coating occasion,coating time and adding amount,and pH. The affecting mechanism of these factors was analyzed further.