Thermodynamic properties for the triangular-well fluid

The thermodynamic properties of the triangular-well fluid with a well range of up to twice the hard sphere diameter were studied by means of a new developed equation of state and molecular simulation. This EoS is based on the perturbation theory of Barker and Henderson with the first and second-order perturbation terms evaluated by molecular simulation and then a fit with a simple function based on the radial distribution function of the reference fluid. The thermodynamic properties for the triangular-well fluid were also obtained directly by Gibbs ensemble and NPT Monte Carlo simulations. Good agreement is observed between the proposed EoS and the molecular simulation results. A model for the triangular-well solid is also presented; this has been used to calculate the solid–liquid transition line. Very good agreement is obtained with previously report values for this line and for the triple point temperature and pressure.     

A new study of associating inhomogeneous fluids with classical density functional theory

ABSTRACT

A new density functional for the study of associating inhomogeneous fluids based on Wertheim's first-order thermodynamic perturbation theory is presented and compared to the most currently used associating density functionals. This functional is developed using the weighted density approximation in the range of association of hard spheres. We implement this functional within the framework of classical density functional theory together with modified fundamental measure theory to account for volume exclusion of hard spheres. This approach is tested against molecular simulations from literature of pure associating hard spheres and mixtures of non-associationg and associating hard spheres with different number of bonding sites close to a hard uniform wall. Furthermore, we compare and review our results with the performance of associating functionals from literature, one based on fundamental measure theory and the inhomogeneous version of Wertheim's perturbation theory. Results obtained with classical DFT and the three functionals show excellent agreement with molecular simulations in systems with one hard wall. For the cases of small pores where only one or two layers of fluid are allowed discrepancies between results with classical DFT and molecular simulations were found.     

Monte Carlo study of dipolar and quadrupolar hard prolate spherocylinders

Monte Carlo simulations of hard prolate spherocylinders (HPSs) with embedded dipole or quadrupole moment are reported for two elongations (L?=?0.5 and 1) and several values of the packing fraction. The MC values of the residual internal and Helmholtz energy and compressibility factor were determined. Our work represents the first simulation study of dipolar HPSs focused on the determination of the thermodynamic properties. In the case of quadrupolar HPSs, our results enlarge the range of state conditions for which the simulation data are available. The obtained MC data were used for a test of the perturbation theory of polar non-spherical molecule fluids. In order to evaluate the perturbation contributions containing the two-particle integrals, the values of the shape integrals (evaluated recently for dipolar and quadrupolar hard prolate spherocylinders) were employed and we were allowed to avoid the use of the similarity between Kihara and Gaussian overlap models. Fair agreement between the simulation data and the theoretical predictions was reached.     

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.     

Effect of the attractive forces on the density change melting of 2D LJ fluids

The Weeks-Chandler-Andersen perturbation theory of fluids can be considered as a modern version of van der Waals'. Whereas in the old van der Waals theory the effect of the attractive forces on the thermodynamic properties is introduced through a constant, a, in Weeks-Chandler-Andersen theory this constant is replaced by a function of temperature and density α(β, ϱ). This function has been determined by means of Molecular Dynamics simulation of a two-dimensional Lennard-Jones fluid, and here is used to analyse the effect of attractive forces on the density change in melting. The influence is found to be slightly greater than the predictions of van der Waals theory.     

Equation of state for the Lennard-Jones truncated and shifted fluid with a cut-off radius of 2.5 σ based on perturbation theory and its applications to interfacial thermodynamics

ABSTRACT

An equation of state is presented for describing thermodynamic properties of the Lennard-Jones truncated and shifted (LJTS) potential with a cut-off radius of 2.5 σ. It is developed using perturbation theory with a hard-sphere reference term and labelled with the acronym PeTS (perturbed truncated and shifted). The PeTS equation of state describes the properties of the bulk liquid and vapour and the corresponding equilibrium of the LJTS fluid well. Furthermore, it is developed so that it can be used safely in the entire metastable and unstable region, which is an advantage compared to existing LJTS equations of state. This makes the PeTS equation of state an interesting candidate for studies of interfacial properties. The PeTS equation of state is applied here in two theories of interfaces, namely density gradient theory (DGT) and density functional theory (DFT). The influence parameter of DGT as well as the interaction averaging diameter of DFT are fitted to data of the surface tension of the LJTS fluid obtained from molecular simulation. The results from both theories agree very well with those from the molecular simulations.     

Generalized Debye-Falkenhagen energy in condensed matter

The thermodynamic perturbation theory of Stell, Rasaiah and Narang is generalized to the case of multi-component Stockmayer fluids. The effects of polarization by the permanent dipoles is included to first order and the perturbation series is summed approximately to all orders in the dipole moments by means of two Padé approximants. The two-body and triplet terms in the expansion have been evaluated by Monte Carlo integration for a specific choice of the reference Lennard-Jones mixture and the results are used to study systematically the effect of dipolar interactions on excess thermodynamic properties. The utility of a one-fluid theory based on the van der Waals model is investigated and it is shown that the approach is useful for obtaining semi-quantitative information on mixtures of weakly polar fluids. The special case of mixtures of hydrogen chloride and xenon is considered and it is shown that the potential model used is inadequate to account for the large deviations from ideality which are observed for this system.     

Thermodynamics of a long-range triangle-well fluid

The long-range triangle-well fluid has been studied using three different approaches: firstly, by an analytical equation of state obtained by a perturbation theory, secondly via a self-consistent integral equation theory, the so-called self-consistent Ornstein–Zernike approach (SCOZA) which is presently one of the most accurate liquid-state theories, and finally by Monte Carlo simulations. We present vapour–liquid phase diagrams and thermodynamic properties such as the internal energy and the pressure as a function of the density at different temperatures and for several values of the potential range. We assess the accuracy of the theoretical approaches by comparison with Monte Carlo simulations: the SCOZA method accurately predicts the thermodynamics of these systems and the first-order perturbation theory reproduces the overall thermodynamic behaviour for ranges greater than two molecular diameters except that it overestimates the critical point. The simplicity of the equation of state and the fact that it is analytical in the potential range makes it a good candidate to be used for calculating other thermodynamic properties and as an ingredient in more complex theoretical approaches.     

Gibbs ensemble simulation of the vapour-liquid equilibrium of square well spherocylinders

Gibbs ensemble Monte Carlo simulations have been performed for systems of square-well spherocylinders of different length-to-breadth ratio. The results are used to test a recent perturbation theory proposed for this kind of system. In addition, the results are compared to similar simulations performed for a Kihara fluid of elongated molecules. An unexpected good agreement is found for the coexistence thermodynamic and structural properties of both model fluids, hence suggesting that the hard spherocylinder plus square-well interaction should be considered as a reference potential for a perturbative treatment of more complex fluid models.     

Simple analytic equations of state for the generalized hard-core Mie(α, β) and Mie(α, β) fluids from perturbation theory

New, simple and analytic perturbation theory equations of state for generalized hard-core Mie HCMie(α, β) and Mie(α, β) fluids are proposed. They are based on the second-order Barker-Henderson perturbation theory in the macroscopic compressibility approximation and the new analytical expression of the radial distribution function of hard spheres, g^HS(r), developed by Sun in terms of a polynomial expansion of base functions adapted to the square-well and Sutherland potentials [Can. J. Phys. 83 (2005) 55], the combination of which yields the HCMie(α, β) and Mie(α, β) functions. The compressibility factors, the residual internal energies and the radial distribution function at contact with the hard core are then obtained from this equation of state for the HCLJ(12, 6) potential, which is a particular case of the HCMie(α, β) potentials with α = 12 and β = 6. The results are in good agreement with the existing Monte Carlo (MC) simulation data, and compare favorably with those obtained from five other equations of state, three of which contain numerical coefficients fitted to the Monte Carlo results. For the Mie(α, 6) (α = 8, 10, 12), fluids, the present equation of state is a good representation of recent molecular dynamics (MD) simulations of the pressure and internal energy. It is more accurate than the statistical associating fluid theory of variable range (SAFT-VR Mie(n, 6)) theory for n = 8, and 10, while for n = 12 the SAFT-VR theory is best. For the Mie(14, 7) fluid, which is outside the range of application of the SAFT-VR theory, the results for the pressure are in good agreement with the analytical equation of state obtained from the MC simulation data.     

Studies of a primitive model of mixtures of nonpolar molecules with water

The paper presents calculations of the properties of binary mixtures of hard spheres and directionally associating hard spheres, a simple model for mixtures of nonpolar molecules with water that was developed by Nezbeda and his coworkers. Extensive results from Monte Carlo simulations in the isobaric, isothermal ensemble are presented for the density, configurational energy and chemical potentials in the mixtures for fluid states over a range of temperatures, pressures and compositions. A species exchange technique is used to compute the chemical potential difference between components in the mixtures. The results obtained are compared with the predictions of first-order thermodynamic perturbation theory (TPT). It is found that this theory provides an accurate picture of the system over most of the conditions considered. Calculations are also made of vapour–liquid coexistence for the model using TPT and calculations of solid–fluid coexistence for the model using TPT and existing results for the free energy of the pure component solids. It is found that the vapour–liquid coexistence for the model is pre-empted by the solid–fluid coexistence, as had previously been found for the pure component directionally associating hard sphere system.     

Equation of state of a supercritical fluid based on the Ornstein-Zernike equation

A numerical modeling of the thermodynamic properties of a fluid is performed using the method of integral equations. The predictions are compared with the results of MC and MD simulations. The problem of stability of the numerical solution is examined. The methods for correcting the correlation functions and for estimating their uncertainties are proposed.     

Quasi-Particle Theory of Alfven Soliton Interaction in Plasmas

The collision of solitons due to Alfven waves in plasmas is studied in this paper by the aid of quasi-particle theory. The suppression of the interaction of solitons, in presence of the perturbation terms, is acheived by means of this theory. The perturbation terms that are considered in this paper are nonlinear damping, finite conductivity and Landau damping. The numerical simulations support the theory that was developed. PACS Codes: 02.30.Ik, 02.30.Jr, 52.35.Sb.     

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.     

Thermodynamic properties of the Lennard-Jones FCC solid: perturbation theory parameterisation and Monte Carlo simulation

This work is concerned with a valid representation of the solid-phase equation of state (EOS), the validity of which is evaluated by comparing to Monte Carlo (MC) simulation results. The proposed EOS has been developed by employing an optimal division of the Lannard-Jones (LJ) potential and an effective temperature- and density-dependent diameter into the framework of the simplified perturbation theory. Then, with the aim of extending to the chain systems, the conventional chain contribution (i.e. TPT1) is added to the proposed model (i.e. the atomic LJ system). Finally, the solid-state EOS based on Helmholtz free energy will be introduced for low temperature and high density conditions. To verify the accuracy of the proposed model, its performance is compared with the results of MC simulation. The comparison between the obtained results from the proposed model and the MC simulations shows that the EOS can satisfactorily predict the properties of the solid LJ system, both for the atomic system and for the chains.     

广义Morse流体解析状态方程及对氮的应用

 根据Ross变分微扰理论以及硬球流体Percus-Yevick（PY） 径向分布函数表达式,建立了广义Morse势流体的解析状态方程。与模拟结果的比较一方面证实了广义Morse 势模型的合理性;另一方面表明了解析Ross变分微扰理论的精度相当或略好于非解析的Weeks-Chandler-Anderson （mWCA）理论,而优于复杂的优化超网络链积分方程理论（RHNC）。该解析状态方程被应用于拟合处于环境温度和压强小于1 GPa情形流体氮的实验数据,所得到的势能参数被用于预测高温高密度情形氮流体的压强,预测结果证实,该解析状态方程可以很好地适用于较宽的压强和温度范围。     

Resummation in hot field theories

There has been significant progress in our understanding of finite-temperature field theory over the past decade. In this paper, we review the progress in perturbative thermal field theory focusing on thermodynamic quantities. We first discuss the breakdown of naive perturbation theory at finite temperature and the need for an effective expansion that resums an infinite class of diagrams in the perturbative expansion. This effective expansion which is due to Braaten and Pisarski, can be used to systematically calculate various static and dynamical quantities as a weak-coupling expansion in powers of g. However, it turns out that the weak-coupling expansion for thermodynamic quantities are useless unless the coupling constant is very small. We critically discuss various ways of reorganizing the perturbative series for thermal field theories in order to improve its convergence. These include screened perturbation theory (SPT), hard-thermal-loop perturbation theory, the Φ-derivable approach, dimensionally reduced (DR) SPT, and the DR Φ-derivable approach.     

Simple optimized Brenner potential for thermodynamic properties of diamond

We have examined the commonly used Brenner potentials in the context of the thermodynamic properties of diamond. A simple optimized Brenner potential is proposed that provides very good predictions of the thermodynamic properties of diamond. It is shown that, compared to the experimental data, the lattice wave theory of molecular dynamics (LWT) with this optimized Brenner potential can accurately predict the temperature dependence of specific heat, lattice constant, Grüneisen parameters and coefficient of thermal expansion (CTE) of diamond.     

Perturbative study of thermodynamic properties for the two-leg spin ladder

We apply finite-temperature perturbation theory to study thermodynamic properties of the two-leg antiferromagnetic spin ladder in the strong interchain coupling limit. The internal energy, specific heat and uniform susceptibility are calculated analytically by third-order perturbation expansions. At zero temperature, the present method results in the same ground state energy as that obtained by the strong coupling expansion without temperature. At finite-temperature, we obtain a peak in the specific heat and a broad maximum in the uniform susceptibility. The results agree quite well with experimental data for the material Cu₂(C₅H₁₂N₂)₂Cl₄ and the numerical data of 8-order series expansion theory.