Self-diffusion of particles interacting through a square-well or square-shoulder potential

The diffusion coefficient and velocity autocorrelation function for a fluid of particles interacting through a square-well or square-shoulder potential are calculated from a kinetic theory similar to the Davis-Rice-Sengers theory and the results are compared to those of computer simulations. At low densities the theory yields too low estimates due to the neglect of correlations between subsequent partial collisions of identical pairs; in particular, the neglect of boundstate effects appears important. At intermediate densities the theory makes reasonable predictions and at high densities it produces too high values, due to the neglect of ring terms and other correlated collision events. The results for the square-shoulder potential generally exhibit better agreement between theory and simulations than do those for the square-well potential.     

On the solvation force of water-like fluid models with square-well attraction and site–site association in slit-like pores: density functional approach

The solvation force of the water-like fluid models with square-well attraction and site–site chemical association confined to slit-like pores has been explored. Theoretical procedure is based on the application of the density functional approach with mean-field approximation for the attractive interparticle interactions. The chemical association effects are treated by using the first-order thermodynamic perturbation theory of Wertheim. Trends of behaviour of the solvation force are put in correspondence with the distribution of molecules in the pores and with the average density of the adsorbate. Moreover, the distribution of non-bonded species on pore width is described. The influence of the width of the square-well and of the gas–solid attraction is discussed. A comparison of theoretical predictions with computer simulations results for water models in slit-like pores is performed.     

The role of higher-order terms in perturbation approaches to the monomer and bonding contributions in a SAFT-type equation of state for square-well chain fluids

A perturbation theory for square-well chain fluids is developed within the scheme of the (generalised) Wertheim thermodynamic perturbation theory. The theory is based on the Pavlyukhin parametrisations [Y. T. Pavlyukhin, J. Struct. Chem. 53, 476 (2012)] of their simulation data for the first four perturbation terms in the high temperature expansion of the Helmholtz free energy of square-well monomer fluids combined with a second-order perturbation theory for the contact value of the radial distribution function of the square-well monomer fluid that enters into bonding contribution. To obtain the latter perturbation terms, we have performed computer simulations in the hard-sphere reference system. The importance of the perturbation terms beyond the second-order one for the monomer fluid and of the approximations of different orders in the bonding contribution for the chain fluids in the predicted equation of state, excess energy and liquid–vapour coexistence densities is analysed.     

Parallelized sampling of the Gibbs ensemble

An expression for the probability distribution of NVT-like sub-ensembles constituting the Gibbs ensemble is derived. Knowledge of this distribution makes it possible to carry out the simulation without the explicit exchange of real particles between the simulation boxes and to evaluate directly any Gibbs ensemble average from a series of independent simultaneous simulations (Monte Carlo or molecular dynamics) performed on a set of NVT-likt sub-ensembles with the fixed distribution of particles. An implementation of the method, which is tailored mainly for complex systems, is exemplified for the square-well fluid, and its efficiency and results are compared with those obtained from the conventional Gibbs ensemble simulations.     

Perturbation theory for the square-well dimer fluid

An analytical equation of state is presented for the square-well dimer fluid of variable well width (1 ≤ λ ≥ 2) based on Barker-Henderson perturbation theory using the recently developed analytical expression for radial distribution function of hard dimers. The integral in the first- and the second-order perturbation terms utilizes the Tang, Y and Lu, B. C.-Y., 1994, J. chem. Phys., 100, 6665 formula for the Hilbert transform. To test the equation of state, NVT and Gibbs ensemble Monte Carlo simulations for square-well dimer fluids are performed for three different well widths (λ = 1.3, 1.5 and 1.8). The prediction of the perturbation theory is also compared with that of thermodynamic perturbation theory in which the equation of state for the square-well dimer is written in terms of that of square-well monomers and the contact value of the radial distribution function.     

A modified soft-core fluid model for the direct correlation function of the square-shoulder and square-well fluids

We propose a simple analytical expression of the direct correlation function for the square-shoulder and square-well fluids. Our approximation is based on an ansatz for the direct correlation function of a modified soft-core fluid, whose parameters are adjusted by fitting the data obtained from Monte Carlo computer simulations. Moreover, it is complemented with a Wertheim-like parametrization to reproduce correctly the direct correlation inside the hard-core. We demonstrate that this approach is in quantitative agreement with the numerical solution of the Ornstein–Zernike equation within the Percus–Yevick approximation. We also show that our results are accurate in a large regime of densities for different interaction ranges and potential strengths. Therefore, this opens up the possibility of introducing the square-shoulder or the square-well potentials as new reference systems in advanced theoretical approximations.     

Monte Carlo simulations in the vicinity of the critical point: Vapour-liquid coexistence curve

A generalized Gibbs-ensemble methodology with a fluctuating particle is used to determine the coexistence vapor-liquid densities for a square-well fluid. It is shown that the presence of the fluctuating particle in sub-states of the Markov chain of states supresses considerably density fluctuations which makes it possible to carry out simulations efficiently even for temperatures very close to the critical temperature.     

Monte Carlo and cell model calculations for the solid—fluid phase behaviour of the triangle-well model

Monte Carlo simulations and cell model calculations are reported for the vapour-liquid and solid-liquid phase behaviour of the triangle-well model system. The behaviour is examined as a function of the range of the triangle-well attraction, from 1.05 to 2.5 times the diameter of the hard core of the potential. Cell model calculations indicate that the stable solid is almost always face-centred cubic (fcc), except for a small set of conditions where hexagonal close-packed (hcp) is favoured. This outcome differs markedly from a much earlier study performed for the square-well model potential, where a much richer phase diagram was observed, with significant regions of stability for hep and body-centred cubic (bcc) phases. Monte Carlo simulations indicate that the cell model calculations represent well the true phase behaviour for this model system. The differing behaviour between the triangle-well and square-well models indicates an important role for the flatness of the potential well in governing the stability of hcp and bcc phases relative to the fcc phase.     

方阱链流体在固液界面分布的密度泛函理论研究

应用Yethiraj的加权密度近似泛函理论研究平板狭缝中方阱链流体的密度分布，系统的Helm holtz自由能泛函分为理想气体的贡献利剩余贡献两部分，其中剩余贡献部分分别采用刘洪 来等人建立的基于空穴相关函数的方阱链流体状态方程和Gil-Villegas等人提出的统计缔合 流体理论状态方程(SAFT-VR)结合简单加权密度近似计算.考察了不同链长、温度、系统密度 和壁面吸引强度下平板狭缝中方阱链流体的密度分布，并与Monte Carlo(MC)模拟结果进行 了比较.结果表明采用不同的状态方程对密度分布的计算有明显的影响，对于受限于硬壁狭 缝中的方阱链流体，温度和密度比较高时，两种状态方程计算的结果均与MC模拟符合得比较 好，在低温和低密度下效果变差，SAFT-VR方程的计算结果更接近于MC模拟结果.对于受限于 方阱壁狭缝中的方阱链流体，由于系统密度分布的非均匀性加强，采用两种状态方程计算的 结果均与MC模拟结果有一定偏差，寻找更合适的权重函数是进一步改进的关键. 关键词： 密度泛函理论 非均匀流体 密度分布 固液界面 方阱链     

Molecular dynamical calculations on the transport properties of a square-well fluid: III. The thermal conductivity

The thermal conductivity of a system consisting of square-well particles has been determined by means of molecular dynamical computer simulations. The calculations were performed for a large number of points in the phase diagram, covering almost the whole fluid region. The effect of the attractive part in the molecular potential is displayed most clearly by looking at the separate contributions of the kinetic, the potential and the cross terms to the coefficient of thermal conductivity. A surprising result is that these terms partially cancel each other, so that the total coefficient of thermal conductivity is rather insensitive to the influence of the attractive well. A comparison with the Davis-Rice-Sengers (DRS) theory shows a discrepancy, be it not as severe as that found in the coefficient of the viscosity.     

Analytical equations for the compressibility factor of a two-dimensional square-well fluid

Analytical approximations for the compressibility factor of a two-dimensional square-well fluid, resulting from perturbational theories are presented.     

Perturbation theory for thermodynamic properties of a ν-dimensional square-well fluid

The Barker-Henderson perturbation theory is used for a ν-dimensional fluid with square-well potential. Analytic expressions are given for the equation of state, excess free energy per particle and internal energy. The numerical results are discussed. A significant feature is the increase of the thermodynamic properties with increasing dimensionality.     

Thermodynamic properties of a two-dimensional square-well fluid

Analytic expressions for the thermodynamic properties of a classical two-dimensional square-well fluid and the first quantum correction to them are derived using the Barker-Henderson perturbation theory. Numerical results are reported. It is found that the quantum effect, which increases with increase of density, is largely determined by the hard-core and the attractive tail has a minor effect at high density.     

Perturbation theory for square-well plus sutherland (∞, 6) fluid

The perturbation approach to the theory of fluids, originally introduced by Zwanzig, is applied to the square-well plus Sutherland (∞, 6) fluid. It is shown that the virial coefficients up to the fourth are in close agreement with those for the LJ (12-6) with a hard-sphere cutoff when the bowl width of the square well used is chosen to be equal to 0.25 σ where σ is the hard-sphere diameter. Also, for this suitable choice of the bowl width, the first and the second-order corrections to the free energy have been found to agree reasonably well with those for the square-well fluid of Barker and Henderson.     

The thermodynamics of molecules with discrete potentials

Fluids formed by molecules interacting with discrete potentials are examined in the context of perturbation theory and the reference hypernetted chain equation (RHNC) solution to the Ornstein—zernike equation. A perturbation theory for discrete-potential fluids (DPT) is presented, which only requires one to know the properties of a square-well fluid of variable range. Several potentials are studied: square-shoulder, a combination of a square-well and square-shoulder, and a discrete representation of a continuous potential model. We have found that the DPT approach reproduces the RHNC predictions in most of the cases.     

Solutions of the Yvon-Born-Green equation for a system of square-well molecules at a temperature below the triple point

The Yvon-Born-Green equation (with superposition approximation) is solved numerically for the pair correlation function for a system of molecules interacting via the square-well potential (with σ₂/σ₁ = 1·85), for an isotherm below the triple point, and over a broad range of densities. The correlation function data and attendant thermodynamics generated for this isotherm are compared with results reported previously by the authors for several supercritical and subcritical isotherms of the square-well fluid. To facilitate the interpretation of these results, particularly in those regions of (T, ρ) space where phase transitions may occur, a geometrical representation of the data is presented (motivated, in part, by recent work by René Thom), and the location of the triple point is discussed in terms of this construction. The differences anticipated between results reported here and those that would be obtained in an exact statistical mechanical analysis, are identified.     

Precise simulation of near-critical fluid coexistence

We present a novel method to derive liquid-gas coexisting densities, rho(+/-)(T), from grand canonical simulations (without knowledge of T(c) or criticality class). The minima of Q(L) identical with (2)(L)/(L) in an LxLxL box with m=rho-(L) are used to generate recursively an unbiased universal finite-size scaling function. Monte Carlo data for a hard-core square-well fluid and for the restricted primitive model electrolyte yield rho(+/-) to +/-1%-2% of rho(c) down to 1 part in 10(4)-10(3) of T(c) (and confirm well Ising character). Pressure mixing in the scaling fields is unequivocally revealed and indicates Yang-Yang ratios R(mu)=-0.04(4) and 0.2(6) for the two models, respectively.     

Quasicrystals in a monodisperse system

We investigate the formation of a two-dimensional quasicrystal in a monodisperse system, using molecular dynamics simulations of hard-sphere particles interacting via a two-dimensional square-well potential. We find that more than one stable crystalline phase can form for certain values of the square-well parameters. Quenching the liquid phase at a very low temperature, we obtain an amorphous phase. By heating this amorphous phase, we obtain a quasicrystalline structure with fivefold symmetry. From estimations of the Helmholtz potentials of the stable crystalline phases and of the quasicrystal, we conclude that the observed quasicrystal phase can be the stable phase in a specific range of temperatures.