基于密度泛函理论研究二元排斥Yukawa流体的表面结构性质

基于自由能密度泛函理论(DFT)考察了二元排斥Yukawa (HCRY)流体在不同外场下的密度分布. 基于微扰理论, 体系的Helmholtz自由能泛函采用硬球排斥部分和长程色散部分贡献之和, 其中Kierlik和Rosinberg的加权密度近似(WDA)被用来计算硬球排斥部分, 而色散部分采用平均场理论(MFT)进行描述. 为了验证DFT计算结果的合理性, 研究中采用巨正则Monte Carlo(GCMC)模拟计算了在不同主体相密度、硬核直径和位能参数比的条件下二元HCRY混合流体的密度分布. 结果表明, 该DFT计算结果与GCMC模拟值吻合良好.     

Balancing local order and long-ranged interactions in the molecular theory of liquid water

A molecular theory of liquid water is identified and studied on the basis of computer simulation of the TIP3P model of liquid water. This theory would be exact for models of liquid water in which the intermolecular interactions vanish outside a finite spatial range, and therefore provides a precise analysis tool for investigating the effects of longer-ranged intermolecular interactions. We show how local order can be introduced through quasichemical theory. Long-ranged interactions are characterized generally by a conditional distribution of binding energies, and this formulation is interpreted as a regularization of the primitive statistical thermodynamic problem. These binding-energy distributions for liquid water are observed to be unimodal. The Gaussian approximation proposed is remarkably successful in predicting the Gibbs free energy and the molar entropy of liquid water, as judged by comparison with numerically exact results. The remaining discrepancies are subtle quantitative problems that do have significant consequences for the thermodynamic properties that distinguish water from many other liquids. The basic subtlety of liquid water is found then in the competition of several effects which must be quantitatively balanced for realistic results.     

Benchmarking the thermodynamic analysis of water molecules around a model beta sheet

Water molecules play a vital role in biological and engineered systems by controlling intermolecular interactions in the aqueous phase. Inhomogeneous fluid solvation theory provides a method to quantify solvent thermodynamics from molecular dynamics or Monte Carlo simulations and provides an insight into intermolecular interactions. In this study, simulations of TIP4P‐2005 and TIP5P‐Ewald water molecules around a model beta sheet are used to investigate the orientational correlations and predicted thermodynamic properties of water molecules at a protein surface. This allows the method to be benchmarked and provides information about the effect of a protein on the thermodynamics of nearby water molecules. The results show that the enthalpy converges with relatively little sampling, but the entropy and thus the free energy require considerably more sampling to converge. The two water models yield a very similar pattern of hydration sites, and these hydration sites have very similar thermodynamic properties, despite notable differences in their orientational preferences. The results also predict that a protein surface affects the free energy of water molecules to a distance of approximately 4.0 Å, which is in line with previous work. In addition, all hydration sites have a favorable free energy with respect to bulk water, but only when the water–water entropy term is included. A new technique for calculating this term is presented and its use is expected to be very important in accurately calculating solvent thermodynamics for quantitative application. © 2012 Wiley Periodicals, Inc.     

Thermodynamic perturbation theory of simple liquids. A hierarchy of equations for expansion coefficients

The properties of an expansion of the statistical sum of a simple liquid with respect to the potential in thermodynamic perturbation theory are analyzed. The coefficients of this expansion are determined by the unperturbed potential, depend on temperature and density, and can be calculated by means of mathematical modeling. It is shown here that the derivatives of these coefficients with respect to temperature and density are expressed through the higher expansion coefficient (these relations are usually called a hierarchy of equations). These coefficients determine the expansion of the Helmholtz free energy and RDF with respect to the perturbation potential. The thermodynamic characteristics of the system (entropy, internal energy, pressure) are expressed through both the differential relations for the Helmholtz free energy and the integral expressions containing RDF. It is found that the hierarchy of equations obtained in this work makes these different methods equivalent. This is important for the application of thermodynamic perturbation theory because it becomes unnecessary to model any other equilibrium properties of the system apart from the expansion coefficients.     

First‐principles modeling of hafnia‐based nanotubes

Hybrid density functional theory calculations were performed for the first time on structure, stability, phonon frequencies, and thermodynamic functions of hafnia‐based single‐wall nanotubes. The nanotubes were rolled up from the thin free layers of cubic and tetragonal phases of HfO₂. It was shown that the most stable HfO₂ single‐wall nanotubes can be obtained from hexagonal (111) layer of the cubic phase. Phonon frequencies have been calculated for different HfO₂ nanolayers and nanotubes to prove the local stability and to find the thermal contributions to their thermodynamic functions. The role of phonons in stability of nanotubes seems to be negligible for the internal energy and noticeable for the Helmholtz free energy. Zone folding approach has been applied to estimate the connection between phonon modes of the layer and nanotubes and to approximate the nanotube thermodynamic properties. It is found that the zone‐folding approximation is sufficiently accurate for heat capacity, but less accurate for entropy. The comparison has been done between the properties of TiO₂, ZrO₂, and HfO₂. © 2017 Wiley Periodicals, Inc.     

Thermodynamic properties of nitrogen from a perturbation theory for ISM fluids

A first-order Barker-Henderson perturbation theory for interaction-site model (ISM) fluids has been applied to calculate the Helmholtz free energy, entropy and internal energy of liquid nitrogen. Comparison with experiment reinforces the idea that the theory is accurate over a wide range of temperatures and densities corresponding to the liquid state, except for the critical region.     

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.     

New parametrization method for dissipative particle dynamics

We introduce an improved method of parametrizing the Groot-Warren version of dissipative particle dynamics (DPD) by exploiting a correspondence between DPD and Scatchard-Hildebrand regular solution theory. The new parametrization scheme widens the realm of applicability of DPD by first removing the restriction of equal repulsive interactions between like beads, and second, by relating all conservative interactions between beads directly to cohesive energy densities. We establish the correspondence by deriving an expression for the Helmoltz free energy of mixing, obtaining a heat of mixing which is exactly the same form as that for a regular mixture (quadratic in the volume fraction) and an entropy of mixing which reduces to the ideal entropy of mixing for equal molar volumes. We equate the conservative interaction parameters in the DPD force law to the cohesive energy densities of the pure fluids, providing an alternative method of calculating the self-interaction parameters as well as a route to the cross interaction parameter. We validate the new parametrization by modeling the binary system SnI(4)SiCl(4), which displays liquid-liquid coexistence below an upper critical solution temperature around 140 degrees C. A series of DPD simulations were conducted at a set of temperatures ranging from 0 degrees C to above the experimental upper critical solution temperature using conservative parameters based on extrapolated experimental data. These simulations can be regarded as being equivalent to a quench from a high temperature to a lower one at constant volume. Our simulations recover the expected phase behavior ranging from solid-liquid coexistence to liquid-liquid coexistence and eventually leading to a homogeneous single phase system. The results yield a binodal curve in close agreement with the one predicted using regular solution theory, but, significantly, in closer agreement with actual solubility measurements.     

Internal and free energy in a pair of like-charged colloids: Monte Carlo simulations

The effective interaction between two colloidal particles in a bath of monovalent co- and counterions is studied by means of lattice Monte Carlo simulations with the primitive model. The internal electrostatic energy as a function of the colloid distance is studied fixing the position of the colloids. The free energy of the whole system is obtained introducing a bias parabolic potential, that allows us to sample efficiently small separations between the colloidal particles. For small charges, both the internal and free energy increase when the colloids approach each other, resulting in an effective repulsion driven by the electrostatic repulsion. When the colloidal charge is large enough, on the other hand, the colloid-ion coupling is strong enough to form double layers. The internal energy in this case decreases upon approaching the colloids because more ions enter the double layer. This attractive contribution to the interaction between the colloids is stronger for larger charges and larger ionic concentrations. However, the total free energy increases due to the loss of ionic entropy, and resulting finally in a repulsive interaction potential driven by the entropic contributions. The loss of ionic entropy can be almost quantitatively reproduced with the ideal contribution, the same level of approximation as the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory. The overall behavior is captured by the DLVO theory qualitatively, and a comparison is made with the functional form predicted by the theory, showing moderate agreement.     

Calculating thermodynamic properties from perturbation theory: Part III. The mixtures of square-well chain fluid

A completely analytic perturbation theory has been developed to calculate the Helmholtz energy, compressibility factor, internal energy and constant-volume heat capacity for square-well chain fluid mixtures. This theory is based on the improved Barker–Henderson macroscopic compressibility (mc) approximation proposed by Zhang, the first-order perturbation theory of Wertheim in which Zhang’s analytic monomer radial distribution function as the function of temperature and monomer density is used, and a simple mixing rule similar to that of Hino–Prausnitz. The validity of the perturbation theory is evaluated by comparing the calculated compressibility factor, internal energy and constant-volume heat capacity for the freely jointed square-well chain mixtures from the theory to MC simulation data. The results show that the theory predicts results in good agreement with simulation results.     

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.     

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.     

Effect of Core-Ion Potentials on Thermodynamic Properties of Liquid Alkali Metals and Alloys

Abstract

A pseudopotential perturbation scheme based on Gibbs-Bogoliubov variational technique is considered to investigate the structure and volume dependent contributions to the internal energies of liquid alkali metals. This has paved the way to analyse the effect of core-ion potentials on Helmholtz free energy, enthalpy, pairwise interactions and on entropy of liquids Na, K, Rb and Cs. Excess entropy of mixing for Na-K. Na-Cs and K-Rb are also calculated and are found very sensitive to the core-ion potentials.     

Sensitivity analysis of thermodynamic properties of liquid water: a general approach to improve empirical potentials

A sensitivity analysis of bulk water thermodynamics is presented in an effort to understand the relation between qualitative features of molecular potentials and properties that they predict. The analysis is incorporated in molecular dynamics simulations and investigates the sensitivity of the Helmholtz free energy, internal energy, entropy, heat capacity, pressure, thermal pressure coefficient, and static dielectric constant to components of the potential rather than the parameters of a given functional form. The sensitivities of the properties are calculated with respect to the van der Waals repulsive and the attractive parts, plus short- and long-range Coulomb parts of three four site empirical water potentials: TIP4P, Dang-Chang and TTM2R. The polarization sensitivity is calculated for the polarizable Dang-Chang and TTM2R potentials. This new type of analysis allows direct comparisons of the sensitivities for different potentials that use different functional forms. The analysis indicates that all investigated properties are most sensitive to the van der Waals repulsive, the short-range Coulomb and the polarization components of the potentials. When polarization is included in the potentials, the magnitude of the sensitivity of the Helmholtz free energy, internal energy, and entropy with respect to this part of the potential is comparable in magnitude to the other electrostatic components. In addition similarities in trends of observed sensitivities for nonpolarizable and polarizable potentials lead to the conclusion that the complexity of the model is not of critical importance for the calculation of these thermodynamic properties for bulk water. The van der Waals attractive and the long-range Coulomb sensitivities are relatively small for the entropy, heat capacity, thermal pressure coefficient and the static dielectric constant, while small changes in any of the potential contributions will significantly affect the pressure. The analysis suggests a procedure for modification of the potentials to improve predictions of thermodynamic properties and we demonstrate this general approach for modifying potentials for one of the potentials.     

A Thermodynamic Model for Square-well Chain Fluid: Theory and Monte Carlo Simulation

A thermodynamic model for the freely jointed square-well chain fluids was developed based on the thermodynamic perturbation theory of Barker-Henderson, Zhang and Wertheim. In this derivation Zhang's expressions for square-well monomers improved from Barker-Henderson compressibility approximation were adopted as the reference fluid, and Wertheim's polymerization method was used to obtain the free energy term due to the bond connectivity. An analytic expression for the Helmholtz free energy of the square-well chain fluids was obtained. The expression without adjustable parameters leads to the thermodynamic consistent predictions of the compressibility factors, residual internal energy and constant-volume heat capacity for dimer, 4-mer, 8-mer and 16-mer square-well fluids. The results are in good agreement with the Monte Carlo simulation. To obtain the MC data of residual internal energy and the constant-volume heat capacity needed, NVT MC simulations were performed for these square-well chain fluids.     

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.     

Ion-water clusters, bulk medium effects, and ion hydration

Thermochemistry of gas-phase ion-water clusters together with estimates of the hydration free energy of the clusters and the water ligands are used to calculate the hydration free energy of the ion. Often the hydration calculations use a continuum model of the solvent. The primitive quasichemical approximation to the quasichemical theory provides a transparent framework to anchor such efforts. Here we evaluate the approximations inherent in the primitive quasichemical approach and elucidate the different roles of the bulk medium. We find that the bulk medium can stabilize configurations of the cluster that are usually not observed in the gas phase, while also simultaneously lowering the excess chemical potential of the ion. This effect is more pronounced for soft ions. Since the coordination number that minimizes the excess chemical potential of the ion is identified as the optimal or most probable coordination number, for such soft ions the optimum cluster size and the hydration thermodynamics obtained with and without account of the bulk medium on the ion-water clustering reaction can be different. The ideas presented in this work are expected to be relevant to experimental studies that translate thermochemistry of ion-water clusters to the thermodynamics of the hydrated ion and to evolving theoretical approaches that combine high-level calculations on clusters with coarse-grained models of the medium.