Semiclassical approach to model quantum fluids using the statistical associating fluid theory for systems with potentials of variable range

Thermodynamic properties of quantum fluids are described using an extended version of the statistical associating fluid theory for potentials of variable range (SAFT-VR) that takes into account quantum corrections to the Helmholtz free energy A, based on the Wentzel-Kramers-Brillouin approximation. We present the theoretical background of this approach (SAFT-VRQ), considering two different cases depending on the continuous or discontinuous nature of the particles pair interaction. For the case of continuous potentials, we demonstrate that the standard Wigner-Kirkwood theory for quantum fluids can be derived from the de Broglie-Bohm formalism for quantum mechanics that can be incorporated within the Barker and Henderson perturbation theory for liquids in a straightforward way. When the particles interact via a discontinuous pair potential, the SAFT-VR method can be combined with the perturbation theory developed by Singh and Sinha [J. Chem. Phys. 67, 3645 (1977); and ibid. 68, 562 (1978)]. We present an analytical expression for the first-order quantum perturbation term for a square-well potential, and the theory is applied to model thermodynamic properties of hydrogen, deuterium, neon, and helium-4. Vapor-liquid equilibrium, liquid and vapor densities, isochoric and isobaric heat capacities, Joule-Thomson coefficients and inversion curves are predicted accurately with respect to experimental data. We find that quantum corrections are important for the global behavior of properties of these fluids and not only for the low-temperature regime. Predictions obtained for hydrogen compare very favorably with respect to cubic equations of state.     

Computer simulation of molecular organization in a nematogen. The role of thermodynamic parameters

Computer simulations of molecular organization in a nematogen, alkenyl bicyclohexylnitrile (ALBCHNL), have been carried out with respect to translational and orientational motions. The atomic net charge and dipole moment at each atomic center has been evaluated using the complete neglect of differential overlap (CNDO/2) method. The modified Rayleigh-Schroedinger perturbation theory, along with multicentered-multipole expansion method, has been employed to evaluate long-range intermolecular interactions, while a “6-exp” potential function has been assumed for short-range interactions. The total interaction energy values obtained through these computations have been used as input to calculate the thermodynamical parameters such as entropy and Helmholtz free energy of each configuration at room temperature (300 K), nematic-isotropic transition (364.7 K) and above transition temperature (450 K). An attempt has been made to understand the molecular organization, to develop a new and interesting model for nematogen based on the thermodynamic parameters introduced in this article.     

Thermodynamic perturbation theory of simple liquids

It was proven that after averaging over the canonical Gibbs ensemble, the mean perturbation energy was singled out of the classical partition function before the expansion in a series of perturbation theory. Therefore, the term that formally coincides with first order perturbation theory in a decomposition of the Helmholtz free energy bears no relationship to perturbation theory. Then the proper series of the thermodynamic perturbation theory always starts with a second order infinitesimal. Therefore, the wellknown condition of applicability of the thermodynamic perturbation theory, “...the requirement that the perturbation energy per particle be small compared with T...” (L. D. Landau and E. M. Livshits, Statistical Physics, Vol. V, Pt. I), can be substantially weakened. The most important factor for applicability of thermodynamic perturbation theory is the value of many-particle correlations in an unperturbed system, but not the smallness of the perturbation potential.     

Thermodynamic properties of polar fluids from a perturbed-dipolar-hard-sphere equation of state

Based on theoretical results for a system of hard spheres with dipoles, a new equation of state is applied to the correlation of thermodynamic properties for four fluids: argon, ammonia, water and acetonitrile. The reference system has the same dependence on density as that given by the Carnahan-Starling equation, but the coefficients are now functions of temperature through the reduced dipole moment. These coefficients are chosen to match the Padé approximant developed by Rushbrooke, Stell and Hoye for the Helmholtz energy of dipolar hard spheres. The reference system proposed here shows a phase transition for reduced dipole moments greater than 1.9. A simple, empirical perturbation term is added to the reference system to account for induction and dispersion forces. For polar fluids, the equation gives results significantly better than those obtained from conventional cubic equations of state, when using the same limited experimental data for determining equation-of-state parameters.     

The changes in free energies and entropies from analytical radial distribution functions

The changes in Helmholtz free energies and entropies in dense fluids have been evaluated using three known analytical expressions for radial distribution functions (RDFs) of Lennard–Jones (L-J) fluid. This method provides a simpler and a more expeditious way for the calculation of free energy and entropy in L-J dense fluids through statistical mechanics. Previously, integral equations or perturbation theories were used for this purpose. Such approach not only tests the power of analytical distribution functions in predicting the changes in Helmholtz free energies and entropies, but also specifies better expressions in determining these properties. The results are compared with experimental data and an accurate analytic equation of state for the L-J fluid. It is shown if an expression properly presents RDFs as a function of interparticle distance, density and temperature, it is possible to calculate the changes in Helmholtz free energies and entropies from analytical distribution functions.     

A perturbation method for the Ornstein-Zernike equation and the generic van der Waals equation of state for a square well potential model

We calculate the generic van der Waals parameters A and B for a square well model by means of a perturbation theory. To calculate the pair distribution function or the cavity function necessary for the calculation of A and B, we have used the Percus-Yevick integral equation, which is put into an equivalent form by means of the Wiener-Hopf method. This latter method produces a pair of integral equations, which are solved by a perturbation method treating the Mayer function or the well width or the functions in the square well region exterior to the hard core as the perturbation. In the end, the Mayer function times the well width is identified as the perturbation parameter in the present method. In this sense, the present perturbation method is distinct from the existing thermodynamic perturbation theory, which expands the Helmholtz free energy in a perturbation series with the inverse temperature treated as an expansion parameter. The generic van der Waals parameters are explicitly calculated in analytic form as functions of reduced temperature and density. The van der Waals parameters are recovered from them in the limits of vanishing density and high temperature. The equation of state thus obtained is tested against Monte Carlo simulation results and found reliable, provided that the temperature is in the supercritical regime. By scaling the packing fraction with a temperature-dependent hard core, it is suggested to construct an equation of state for fluids with a temperature-dependent hard core that mimicks a soft core repulsive force on the basis of the equation of state derived for the square well model.     

Thermodynamic perturbation theory of simple liquids: Monte Carlo simulation of a hard sphere system and the Helmholtz free energy of SW fluids

The Monte Carlo method is used to simulate similar sized hard sphere systems in a wide range of densities (from η 0.005 to 0.530 with a step of η = 0.005). The models are used to calculate the coefficients of the thermodynamic perturbation theory series for SW fluids up to the fourth order. The width of the attraction zone of the SW potential λ is varied from 1 to 2.5 sphere diameters. The analytical expressions approximating the obtained coefficients by polynomials with respect to the variables η and λ are determined. The absolute accuracy of the approximation is estimated to be better than ±0.001. All the necessary data for the calculation of the Helmholtz free energy of SW fluids up to the fourth-order perturbation theory are given.     

Resummed thermodynamic perturbation theory for central force associating potential. Multi-patch models

A resummed thermodynamic perturbation theory for associating fluids with multiply bondable central force associating potential is extended for the fluid with multiple number of multiply bondable associating sites. We consider a multi-patch hard-sphere model for associating fluids. The model is represented by the hard-sphere fluid system with several spherical attractive patches on the surface of each hard sphere. Resummation is carried out to account for blocking effects, i.e., when the bonding of a particle restricts (blocks) its ability to bond with other particles. Closed form analytical expressions for thermodynamical properties (Helmholtz free energy, pressure, internal energy, and chemical potential) of the models with arbitrary number of doubly bondable patches at all degrees of the blockage are presented. In the limiting case of total blockage, when the patches become only singly bondable, our theory reduces to Wertheim's thermodynamic perturbation theory (TPT) for polymerizing fluids. To validate the accuracy of the theory we compare to exact values, for the thermodynamical properties of the system, as determined by Monte Carlo computer simulations. In addition we compare the fraction of multiply bonded particles at different values of the density and temperature. In general, predictions of the present theory are in good agreement with values for the model calculated using Monte Carlo simulations, i.e., the accuracy of our theory in the case of the models with multiply bondable sites is similar to that of Wertheim's TPT in the case of the models with singly bondable sites.     

Theoretical investigation on the structural and thermodynamic properties of FeSe at high pressure and high temperature

A theoretical investigation on structural and thermodynamic properties of 11-type iron-based superconductor FeSe at high pressure and high temperature was performed by employing the first-principles method based on the density functional theory. Some structural parameters of FeSe in both tetragonal and hexagonal phases are reported. According to the fourth-order Birch-Murnaghan equation of states, the transition pressure P(t) of FeSe from the PbO-type phase to the NiAs-type phase was determined. The calculated results are found to be in good agreement with the available experimental data. Based on the quasi-harmonic Debye model, the pressure and temperature dependence of the thermodynamic properties for hexagonal phase FeSe were investigated. Our theoretical calculations suggest that the pressure and temperature have significant effects on the heat capacity, vibrational internal energy, vibrational entropy, vibrational Helmholtz free energy, thermal expansion coefficient and Debye temperature. Even though few theoretical reports on the structural properties of FeSe are found in the current literature, to our knowledge, this is a novel theoretical investigation on the structural and thermodynamic properties of FeSe at high temperature. We hope that the theoretical results reported here can give more insight into the structural and thermodynamic properties of other iron-based superconductors at high temperature.     

Hard-sphere perturbation theory for a model of liquid Ga

Investigating thermodynamic properties of a model for liquid Ga, we have extended the application of the hard-sphere (HS) perturbation theory to an interatomic pair potential that possesses a soft repulsive core and a long-range oscillatory part. The model is interesting for displaying a discontinuous jump on the main-peak position of the radial distribution function at some critical density. At densities less than this critical value, the effective HS diameter of the model, estimated by the variational HS perturbation theory, has a substantial reduction with increasing density. Thus, the density dependence of the packing fraction of the HS reference fluid has an anomalous behavior, with a negative slope, within a density region below the critical density. By adding a correction term originally proposed by Mon to remedy the inherent deficiency of the HS perturbation theory, the extended Mansoori-Canfield/Rasaiah-Stell theory [J. Chem. Phys. 120, 4844 (2004)] very accurately predicts the Helmholtz free energy and entropy of the model, including an excess entropy anomaly. Almost occurring in the same density region, the excess entropy anomaly is found to be associated with the anomalous packing faction of the HS fluid.     

Toward a statistical mechanical theory for water: analytical theory for a short-ranged reference system

Starting from a realistic Hamiltonian and making use of recent findings that the properties of associating fluids are determined primarily by short-ranged interactions, this methodology has been implemented using statistical mechanical approaches and thermodynamic perturbation theory for the TIP4P model of water. We focus on the short-range reference system for which an analytic expression for the Helmholtz free energy is derived. It is found that the model (reference system) exhibits, in addition to a faithful representation of the structure of water, the same features that are characteristic for real water, namely, (i) the temperature of the density maximum and its pressure dependence, including the inflection point at high pressures and (ii) the temperature minima of the constant pressure heat capacity and the coefficient of isothermal compressibility.     

用密度泛函理论研究Lennard-Jones 流体在狭缝中的相平衡

用改进的基础度量理论(modified fundamental measure theory, MFMT)和密度Taylor展开分别表达过剩自由能中的短程作用和色散作用. 流体分子与狭缝壁之间的相互作用以10-4-3势能函数表达. 由巨势最小原理确定Lennard-Jones (LJ)流体在狭缝中的密度分布和过剩吸附量, 所得结果与分子模拟数据吻合良好. 根据平衡时两相温度, 化学势及巨势相等, 计算了LJ流体在狭缝中的相平衡.     

Finite-temperature full configuration interaction

The exact basis-set values of various thermodynamic potentials of a molecule are evaluated by the finite-temperature full configuration-interaction (FCI) method using ab initio molecular integrals over Gaussian-type orbitals. The thermodynamic potentials considered are the grand partition function, grand potential, internal energy, entropy, and chemical potential in the grand canonical ensemble as well as the partition function, Helmholtz energy, internal energy, and entropy in canonical ensemble. Approximations to FCI that are accurate at low and high temperatures are proposed, implemented, and tested. The results of finite-temperature FCI and its approximations are compared with one another as well as with the results of finite-temperature zeroth-order many-body perturbation theory, in which the Fermi–Dirac statistics is exact. Analytical asymptotic properties in the low- or high-temperature limits of some of these thermodynamic potentials are also given.     

Phase equilibria and plate-fluid interfacial tensions for associating hard sphere fluids confined in slit pores

The excess Helmholtz free energy functional for associating hard sphere fluid is formulated by using a modified fundamental measure theory [Y. X. Yu and J. Z. Wu, J. Chem. Phys. 117, 10156 (2002)]. Within the framework of density functional theory, the thermodynamic properties including phase equilibria for both molecules and monomers, equilibrium plate-fluid interfacial tensions and isotherms of excess adsorption, average molecule density, average monomer density, and plate-fluid interfacial tension for four-site associating hard sphere fluids confined in slit pores are investigated. The phase equilibria inside the hard slit pores and attractive slit pores are determined according to the requirement that temperature, chemical potential, and grand potential in coexistence phases should be equal and the plate-fluid interfacial tensions at equilibrium states are predicted consequently. The influences of association energy, fluid-solid interaction, and pore width on phase equilibria and equilibrium plate-fluid interfacial tensions are discussed.     

Evaluating perturbation contributions in SAFT models by comparing to molecular simulation of n-alkanes

SAFT models are generally written as a perturbation series of the Helmholtz energy with reciprocal temperature as the argument. The perturbation coefficients are then functions of density and molecular size. The variation of the perturbation coefficients with molecular size is given primarily by Wertheim's theory [6], [7], [8] and [9], but there may be additional variations as in the PC-SAFT model. In the present work, we compare the characterization of perturbation coefficients inferred from PC-SAFT to those derived from molecular simulations.The molecular simulations are based on Discontinuous Molecular Dynamics (DMD) and second order Thermodynamic Perturbation Theory (TPT). DMD simulation is applied to the repulsive part of the potential model with molecular details like fused hard spheres for the interaction sites and 110° bond angles. The thermodynamic effects of disperse attractions are treated by rigorous application of TPT. The present work re-examines the related work of Elliott and Gray [35] in the low density and critical regions, focusing on n-alkanes with carbon numbers ranging from 3 to 80.We find that SAFT theory overestimates the repulsive contribution (A₀) and underestimates the first order contribution (A₁) of Helmholtz energy relative to simulation. Nevertheless, the correlations are qualitatively reasonable. Significant inconsistencies arise when considering the second order contribution (A₂). For example, the PC-SAFT characterization of A₂ becomes larger than A₁ in the low density, long chain limit, raising concerns about the convergence of the series. Furthermore, fluctuations are underestimated in the critical region and overestimated in the liquid region. In each case, we can suggest improved characterizations. Altogether, these results suggest ways to modify the SAFT formalism to achieve greater consistency between atomistic and coarse-grained models.     

Density functional study on the structural and thermodynamic properties of aqueous DNA-electrolyte solution in the framework of cell model

A density functional theory (DFT) in the framework of cell model is proposed to calculate the structural and thermodynamic properties of aqueous DNA-electrolyte solution with finite DNA concentrations. The hard-sphere contribution to the excess Helmholtz energy functional is derived from the modified fundamental measure theory, and the electrostatic interaction is evaluated through a quadratic functional Taylor expansion around a uniform fluid. The electroneutrality in the cell leads to a variational equation with a constraint. Since the reference fluid is selected to be a bulk phase, the Lagrange multiplier proves to be the potential drop across the cell boundary (Donnan potential). The ion profiles and electrostatic potential profiles in the cell are calculated from the present DFT-cell model. Our DFT-cell model gives better prediction of ion profiles than the Poisson-Boltzmann (PB)- or modified PB-cell models when compared to the molecular simulation data. The effects of polyelectrolyte concentration, ion size, and added-salt concentration on the electrostatic potential difference between the DNA surface and the cell boundary are investigated. The expression of osmotic coefficient is derived from the general formula of grand potential. The osmotic coefficients predicted by the DFT are lower than the PB results and are closer to the simulation results and experimental data.     

Density-functional theory and Monte Carlo simulation for the surface structure and correlation functions of freely jointed Lennard-Jones polymeric fluids

We present a nonlocal density-functional theory of polymeric fluids consisting of freely jointed Lennard-Jones chains with explicit consideration of the segment size, van der Waals attraction, and structural correlations due to chain connectivity. The excess Helmholtz energy functional is derived from a modified fundamental measure theory for the short-ranged repulsion and the first-order thermodynamic perturbation theory for chain connectivity. The contribution of the long-ranged attraction to the Helmholtz energy functional is taken into account using a quadratic density expansion with the direct correlation function obtained from the first-order mean-spherical approximation. The numerical performance of the density-functional theory is compared well with the simulation results from this work as well as those from the literature for the segment-level density profiles and correlation functions of Lennard-Jones chains in slit pores, near isolated nanoparticles, or in bulk.     

用重整化群理论研究胶体模型体系的相行为

用yukawa势能函数表达胶体颗粒之间的吸引作用。用Duh-Mier-Y-Teran状态方 程表达液相Helmholtz自由能。用一阶微扰理论、固体硬球径向分布函数解析式和 改进的胞腔模型建立固相状态方程，结合建立的状态方程和重整化群理论。研究了 胶体模型体系的液-液相平衡和液-固相平衡。研究表明，颗粒之间色散作用量程参 数的变化对胶休到本世纪末茶杯 系的相行为有特殊需要影响。所得结果与分子模 拟数据吻合良好。