Waterlike structural and excess entropy anomalies in liquid beryllium fluoride

The relationship between structural order metrics and the excess entropy is studied using the transferable rigid ion model (TRIM) of beryllium fluoride melt, which is known to display waterlike thermodynamic anomalies. The order map for liquid BeF2, plotted between translational and tetrahedral order metrics, shows a structurally anomalous regime, similar to that seen in water and silica melt, corresponding to a band of state points for which average tetrahedral (q(tet)) and translational (tau) order are strongly correlated. The tetrahedral order parameter distributions further substantiate the analogous structural properties of BeF2, SiO2, and H2O. A region of excess entropy anomaly can be defined within which the pair correlation contribution to the excess entropy (S2) shows an anomalous rise with isothermal compression. Within this region of anomalous entropy behavior, q(tet) and S2 display a strong negative correlation, indicating the connection between the thermodynamic and the structural anomalies. The existence of this region of excess entropy anomaly must play an important role in determining the existence of diffusional and mobility anomalies, given the excess entropy scaling of transport properties observed in many liquids.     

Relationships between self-diffusivity, packing fraction, and excess entropy in simple bulk and confined fluids

Static measures such as density and entropy, which are intimately connected to structure, have featured prominently in modern thinking about the dynamics of the liquid state. Here, we explore the connections between self-diffusivity, density, and excess entropy for two of the most widely used model "simple" liquids, the equilibrium Lennard-Jones and square-well fluids, in both bulk and confined environments. We find that the self-diffusivity data of the Lennard-Jones fluid can be approximately collapsed onto a single curve (i) versus effective packing fraction and (ii) in appropriately reduced form versus excess entropy, as suggested by two well-known scaling laws. Similar data collapse does not occur for the square-well fluid, a fact that can be understood on the basis of the nontrivial effects that temperature has on its static structure. Nonetheless, we show that the implications of confinement for the self-diffusivity of both of these model fluids, over a broad range of equilibrium conditions, can be predicted on the basis of knowledge of the bulk fluid behavior and either the effective packing fraction or the excess entropy of the confined fluid. Excess entropy is perhaps the most preferable route due to its superior predictive ability and because it is a standard, unambiguous thermodynamic quantity that can be readily predicted via classical density functional theories of inhomogeneous fluids.     

Core-softened fluids, water-like anomalies, and the liquid-liquid critical points

Molecular dynamics simulations are used to examine the relationship between water-like anomalies and the liquid-liquid critical point in a family of model fluids with multi-Gaussian, core-softened pair interactions. The core-softened pair interactions have two length scales, such that the longer length scale associated with a shallow, attractive well is kept constant while the shorter length scale associated with the repulsive shoulder is varied from an inflection point to a minimum of progressively increasing depth. The maximum depth of the shoulder well is chosen so that the resulting potential reproduces the oxygen-oxygen radial distribution function of the ST4 model of water. As the shoulder well depth increases, the pressure required to form the high density liquid decreases and the temperature up to which the high-density liquid is stable increases, resulting in the shift of the liquid-liquid critical point to much lower pressures and higher temperatures. To understand the entropic effects associated with the changes in the interaction potential, the pair correlation entropy is computed to show that the excess entropy anomaly diminishes when the shoulder well depth increases. Excess entropy scaling of diffusivity in this class of fluids is demonstrated, showing that decreasing strength of the excess entropy anomaly with increasing shoulder depth results in the progressive loss of water-like thermodynamic, structural and transport anomalies. Instantaneous normal mode analysis was used to index the overall curvature distribution of the fluid and the fraction of imaginary frequency modes was shown to correlate well with the anomalous behavior of the diffusivity and the pair correlation entropy. The results suggest in the case of core-softened potentials, in addition to the presence of two length scales, energetic, and entropic effects associated with local minima and curvatures of the pair interaction play an important role in determining the presence of water-like anomalies and the liquid-liquid phase transition.     

Waterlike hierarchy of anomalies in a continuous spherical shouldered potential

We investigate by molecular dynamics simulations a continuous isotropic core-softened potential with attractive well in three dimensions, introduced by Franzese [J. Mol. Liq. 136, 267 (2007)], that displays liquid-liquid coexistence with a critical point and waterlike density anomaly. Besides the thermodynamic anomalies, here we find diffusion and structural anomalies. The anomalies, not observed in the discrete version of this model, occur with the same hierarchy that characterizes water. We discuss the differences in the anomalous behavior of the continuous and the discrete model in the framework of the excess entropy, calculated within the pair correlation approximation.     

Non-hard sphere thermodynamic perturbation theory

A non-hard sphere (HS) perturbation scheme, recently advanced by the present author, is elaborated for several technical matters, which are key mathematical details for implementation of the non-HS perturbation scheme in a coupling parameter expansion (CPE) thermodynamic perturbation framework. NVT-Monte Carlo simulation is carried out for a generalized Lennard-Jones (LJ) 2n-n potential to obtain routine thermodynamic quantities such as excess internal energy, pressure, excess chemical potential, excess Helmholtz free energy, and excess constant volume heat capacity. Then, these new simulation data, and available simulation data in literatures about a hard core attractive Yukawa fluid and a Sutherland fluid, are used to test the non-HS CPE 3rd-order thermodynamic perturbation theory (TPT) and give a comparison between the non-HS CPE 3rd-order TPT and other theoretical approaches. It is indicated that the non-HS CPE 3rd-order TPT is superior to other traditional TPT such as van der Waals/HS (vdW/HS), perturbation theory 2 (PT2)/HS, and vdW/Yukawa (vdW/Y) theory or analytical equation of state such as mean spherical approximation (MSA)-equation of state and is at least comparable to several currently the most accurate Ornstein-Zernike integral equation theories. It is discovered that three technical issues, i.e., opening up new bridge function approximation for the reference potential, choosing proper reference potential, and/or using proper thermodynamic route for calculation of f(ex-ref), chiefly decide the quality of the non-HS CPE TPT. Considering that the non-HS perturbation scheme applies for a wide variety of model fluids, and its implementation in the CPE thermodynamic perturbation framework is amenable to high-order truncation, the non-HS CPE 3rd-order or higher order TPT will be more promising once the above-mentioned three technological advances are established.     

Confinement, entropy, and single-particle dynamics of equilibrium hard-sphere mixtures

We use discontinuous molecular dynamics and grand-canonical transition-matrix Monte Carlo simulations to explore how confinement between parallel hard walls modifies the relationships between packing fraction, self-diffusivity, partial molar excess entropy, and total excess entropy for binary hard-sphere mixtures. To accomplish this, we introduce an efficient algorithm to calculate partial molar excess entropies from the transition-matrix Monte Carlo simulation data. We find that the species-dependent self-diffusivities of confined fluids are very similar to those of the bulk mixture if compared at the same, appropriately defined, packing fraction up to intermediate values, but then deviate negatively from the bulk behavior at higher packing fractions. On the other hand, the relationships between self-diffusivity and partial molar excess entropy (or total excess entropy) observed in the bulk fluid are preserved under confinement even at relatively high packing fractions and for different mixture compositions. This suggests that the excess entropy, calculable from classical density functional theories of inhomogeneous fluids, can be used to predict some of the nontrivial dynamical behaviors of fluid mixtures in confined environments.     

Application of a renormalization-group treatment to the statistical associating fluid theory for potentials of variable range (SAFT-VR)

An accurate prediction of phase behavior at conditions far and close to criticality cannot be accomplished by mean-field based theories that do not incorporate long-range density fluctuations. A treatment based on renormalization-group (RG) theory as developed by White and co-workers has proven to be very successful in improving the predictions of the critical region with different equations of state. The basis of the method is an iterative procedure to account for contributions to the free energy of density fluctuations of increasing wavelengths. The RG method has been combined with a number of versions of the statistical associating fluid theory (SAFT), by implementing White's earliest ideas with the improvements of Prausnitz and co-workers. Typically, this treatment involves two adjustable parameters: a cutoff wavelength L for density fluctuations and an average gradient of the wavelet function Φ. In this work, the SAFT-VR (variable range) equation of state is extended with a similar crossover treatment which, however, follows closely the most recent improvements introduced by White. The interpretation of White's latter developments allows us to establish a straightforward method which enables Φ to be evaluated; only the cutoff wavelength L then needs to be adjusted. The approach used here begins with an initial free energy incorporating only contributions from short-wavelength fluctuations, which are treated locally. The contribution from long-wavelength fluctuations is incorporated through an iterative procedure based on attractive interactions which incorporate the structure of the fluid following the ideas of perturbation theories and using a mapping that allows integration of the radial distribution function. Good agreement close and far from the critical region is obtained using a unique fitted parameter L that can be easily related to the range of the potential. In this way the thermodynamic properties of a square-well (SW) fluid are given by the same number of independent intermolecular model parameters as in the classical equation. Far from the critical region the approach provides the correct limiting behavior reducing to the classical equation (SAFT-VR). In the critical region the β critical exponent is calculated and is found to take values close to the universal value. In SAFT-VR the free energy of an associating chain fluid is obtained following the thermodynamic perturbation theory of Wertheim from the knowledge of the free energy and radial distribution function of a reference monomer fluid. By determining L for SW fluids of varying well width a unique equation of state is obtained for chain and associating systems without further adjustment of critical parameters. We use computer simulation data of the phase behavior of chain and associating SW fluids to test the accuracy of the new equation.     

Density dependence of the entropy and the solvation shell structure in supercritical water via molecular dynamics simulation

We perform molecular dynamics simulations of supercritical water (SCW) with a wide range of densities along a near critical isotherm using the simple point charge extended (SPC/E) pair potential in order to study the entropy and the solvation shell structure around a central water molecule. It is shown that both the translational and orientational two-particle correlation entropy terms can serve as the metrics of the translational-orientational structural orders in water and it is revealed that the translational structural order is very sensitive to the density variation in the gas-like and liquid-like region, while the orientational structural order is much more dependent upon compression in the medium-density SCW region. The comparison of the magnitudes of the full thermodynamic excess entropy and two-particle correlation entropy confirms the recent findings that the many-body terms other than two-body ones also present significant and non-neglectable contributions to the full excess entropy for the highly anomalous fluids like water. The analysis of entropy terms as a function of intermolecular distance and the orientational distribution functions as well as the three-dimensional spatial distribution functions indicate that the structural order occurs only in a much more diffused first solvation shell due to the elongated hydrogen bonds under supercritical conditions. It is revealed that no obvious second or higher neighbor shells occur in SCW, in contrast with the feature of normal liquid water that the anomalous decrease of translational order upon compression occurs mainly in the second shell.     

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.     

How critical fluctuations influence adsorption properties of a van der Waals fluid onto a spherical colloidal particle

A recently proposed 3rd-order thermodynamic perturbation theory (TPT) is extended to its 5th-order version and non-uniform counterpart by supplementing with density functional theory (DFT) and a number of ansatzs for a bulk 2nd-order direct correlation function (DCF). Employment of the ansatzs DCF enables the resultant non-uniform formalism devoid of any adjustable parameter and free from numerically solving an Ornstein–Zernike integral equation theory. Density profiles calculated by the present non-uniform formalism for a hard core attractive Yukawa (HCAY) fluid near a spherical geometry are favorably compared with corresponding simulation data available in literature, and are more accurate than those based on a previous 3rd + 2nd-order perturbation DFT. The non-uniform 5th-order TPT is employed to investigate adsorption of the HCAY fluid onto a colloidal particle; it is disclosed that a depletion adsorption can be induced when the coexistence bulk fluid is situated in neighborhood of a critical point or near a bulk vapor–liquid coexistence gaseous phase or liquid phase density. A physical interpretation is given for such depletion adsorption and for its connection with parameters of the potential under consideration, which is ascribed to critical density fluctuations existing within a wide region of the bulk diagram. For a large spherical external potential inducing wetting transition, it is found that only round wetting transition is found instead of 1st-order pre-wetting transition in the case of a planar wall external potential, and the wetting transition temperature increases relative to that for the planar wall external potential. The present theoretical results for wetting transitions are supported by previous investigation based on thermodynamic considerations and a phenomenological Landau mean field theory, and are also in conformity with the present qualitative physical interpretation.     

Does confining the hard-sphere fluid between hard walls change its average properties?

We use grand canonical transition-matrix Monte Carlo and discontinuous molecular dynamics simulations to generate precise thermodynamic and kinetic data for the equilibrium hard-sphere fluid confined between smooth hard walls. These simulations show that the pronounced inhomogeneous structuring of the fluid normal to the confining walls, often the primary focus of density functional theory studies, has a negligible effect on many of its average properties over a surprisingly broad range of conditions. We present one consequence of this insensitivity to confinement: a simple analytical equation relating the average density of the confined fluid to that of the bulk fluid with equal activity. Nontrivial implications of confinement for average fluid properties do emerge in this system, but only when the fluid is both (i) dense and (ii) confined to a gap smaller than approximately three particle diameters. For this limited set of conditions, we find that "in-phase" oscillatory deviations in excess entropy and self-diffusivity (relative to the behavior of the bulk fluid at the same average density) occur as a function of gap size. These paired thermodynamic/kinetic deviations from bulk behavior appear to reflect the geometric packing frustration that arises when the confined space cannot naturally accommodate an integer number of particle layers.     

Thermodynamic Properties of Multicomponent Liquid Mixtures

Considering a ternary liquid mixture to be made up of three binary mixtures, by means of cell model using Sutherland type potential function for pair interaction between molecules, a statistical theory for binary liquid mixtures has been extended for ternary systems. In the light of above extension, excess volume ( V E ), excess energy ( E E ) and excess entropy ( TS E ) have been computed for three binary (benzene + cyclohexane, benzene + chlorobenzene and cyclohexane + chlorobenzene) and the resultant ternary system (benzene + cyclohexane + chlorobenzene) at 298.15 K. All the above mentioned excess properties have been computed from the data of ultrasonic velocity and density only.     

Structural anomalies for a three dimensional isotropic core-softened potential

Using molecular dynamics simulations we investigate the structure of a system of particles interacting through a continuous core-softened interparticle potential. We found for the translational order parameter t a local maximum at a density rho(t-max) and a local minimum at rho(t-min)>rho(t-max). Between rho(t-max) and rho(t-min), the t parameter anomalously decreases upon increasing pressure. For the orientational order parameter Q(6) a maximum was observed at a density rho(t-max)相似文献   

Core-softened system with attraction: trajectory dependence of anomalous behavior

In the present article we carry out a molecular dynamics study of the core-softened system and show that the existence of the water-like anomalies in this system depends on the trajectory in P-ρ-T space along which the behavior of the system is studied. For example, diffusion and structural anomalies are visible along isotherms as a function of density, but disappears along the isochores and isobars as a function of temperature. On the other hand, the diffusion anomaly may be seen along adiabats as a function of temperature, density, and pressure. It should be noted that it may be no signature of a particular anomaly along a particular trajectory, but the anomalous region for that particular anomaly can be defined when all possible trajectories in the same space are examined (for example, signature of diffusion anomaly is evident through the crossing of different isochors. However, there is no signature of diffusion anomaly along a particular isochor). We also analyze the applicability of the Rosenfeld entropy scaling relations to this system in the regions with the water-like anomalies. It is shown that the validity of the Rosenfeld scaling relation for the diffusion coefficient also depends on the trajectory in the P-ρ-T space along which the kinetic coefficients and the excess entropy are calculated.     

Extending the simple weighted density approximation for a hard-sphere fluid to a Lennard-Jones fluid II. Application

A simple weighted density approximation (SWDA) was extended to nonuniform Lennard-Jones fluids by following the spirit of a partitioned density function theory [S. Zhou, Phys. Rev. E 68 (2003) 061201] and mapping the hard-core part onto an effective hard-sphere fluid whose higher order terms beyond the second order of the functional perturbation expansion are treated by the SWDA. The resultant DFT formalism performs well for Lennard-Jones fluids under the influence of diverse external fields. With the present DFT formalism, we investigate in detail the structure and adsorption properties of a low-density LJ gas in a spherical cavity with a wall consisting of hard-sphere or LJ particles. It was found that when the cavity wall exerts an attractive external potential on the LJ particles in the cavity, the excess adsorption decreases as the temperature increases, while when the cavity wall exerts a hard repulsive external potential on the LJ particles in the cavity, the excess adsorption increases as the temperature increases.     

Simplified crossover droplet model for adsorption of pure fluids in slit pores

We present a generalized crossover (GC) model for the excess adsorption of pure fluids at a flat solid-liquid interface, which reproduces scaling behavior of the excess adsorption in the critical region and is reduced to the classical, van der Waals-type analytical model far away from the bulk critical point. In developing this model, we used the density-functional theory (DFT) approach for the order parameter profile calculations with a generalized corresponding states model for the local free-energy density. The GC DFT model well represents the available experimental adsorption data for Kr/graphite, C2H4/graphite, C3H8/graphite, CO2/silica, and SF6/graphite systems in the entire density range 0 < rho < or = 3rhoc and temperatures up to 1.7Tc. In the critical region 0.5 rhoc < r < or = 1.5rhoc and T < or = 1.15Tc, the GC DFT model is consistent with the predictions of the asymptotic renormalization-group crossover model for the critical adsorption in a semi-infinite system developed earlier. For the excess adsorption on the critical isochore, both theories predict a scaling-law behavior Gamma proportional tau(-nu+beta), but fail to reproduce a "critical depletion" of the excess adsorption along the critical isochore of the SF6/graphite system near Tc. We show that an anomalous decrease of adsorption observed in this system at tau = T/Tc - 1 < 10(-2) can be explained by finite-size effect and develop a simplified crossover droplet (SCD) model for the excess adsorption in a slit pore. With the effective size of the pore of L = 50 nm, the SCD model reproduces all available experimental data for SF6/graphite, including the critical isochore data where tau-->0, within experimental accuracy. At L > xib (where xib is a bulk correlation length) the SCD model is transformed into the GC DFT model for semi-infinite systems. Application of the SCD model to the excess adsorption of carbon dioxide on the silica gel is also discussed.     

Packing density of HS(CH2)(n)COOH self-assembled monolayers

Self-assembled monolayers (SAMs) of HS(CH(2))(n)COOH, n = 5, 10, 15 deposited from ethanol solution onto gold are prepared by five approaches, and their packing densities are evaluated by X-ray photoelectron spectroscopy (XPS) measurements. The five approaches are: (1) direct deposition; (2) acetic-acid-assisted deposition; (3) butyl-amine-assisted deposition; (4) displacement of a preformed HS(CH(2))(n)CH(3) (n = 5, 10, 15) SAMs; and (5) co-deposition with HS(CH(2))(n)CH(3) (n = 5, 10, 15). Packing density metrics are calculated from measurements of SAM and substrate photoemission intensities and their attenuations by two methods. In one case the attenuated photoemissions are expressed as a ratio relative to comparable measurements on an experimental HS(CH(2))(n)CH(3) model system. In the other case a new method is introduced where a calculated attenuation based on theoretical random coil and extended chain models is used as the reference to determine a packing density fraction. Packing densities are also correlated with the S2p(Au-bonded):Au4f peak area ratios and with shifts in the C1s binding energies. SAMs prepared by the direct deposition are a partial multilayer where a second molecular layer is physisorbed onto the SAM and not removable by solvent washing. The addition of acetic acid to the deposition solution disrupts dimer associations of HS(CH(2))(n)COOH in solution and at the surface of the monolayer and yields the most ordered monolayer with the highest density of -COOH groups. The addition of butyl amine results in a labile ammonium carbonate ion pair formation but results in a lower packing density in the SAM. The displacement of the preformed HS(CH(2))(n)CH(3) SAM and the co-deposition of HS(CH(2))(n)CH(3) with HS(CH(2))(n)COOH result in SAMs with little incorporation of the -COOH component.     

Vapor-liquid equilibrium and critical behavior of the square-well fluid of variable range: a theoretical study

The vapor-liquid phase behavior and the critical behavior of the square-well (SW) fluid are investigated as a function of the interaction range, lambdain [1.25, 3], by means of the self-consistent Ornstein-Zernike approximation (SCOZA) and analytical equations of state based on a perturbation theory [A. L. Benavides and F. del Rio, Mol. Phys. 68, 983 (1989); A. Gil-Villegas, F. del Rio, and A. L. Benavides, Fluid Phase Equilib. 119, 97 (1996)]. For this purpose the SCOZA, which has been restricted up to now to a few model systems, has been generalized to hard-core systems with arbitrary interaction potentials requiring a fully numerical solution of an integro-partial differential equation. Both approaches, in general, describe well the liquid-vapor phase diagram of the square-well fluid when compared with simulation data. SCOZA yields very precise predictions for the coexistence curves in the case of long ranged SW interaction (lambda>1.5), and the perturbation theory is able to predict the binodal curves and the saturated pressures, for all interaction ranges considered if one stays away from the critical region. In all cases, the SCOZA gives very good predictions for the critical temperatures and the critical pressures, while the perturbation theory approach tends to slightly overestimate these quantities. Furthermore, we propose analytical expressions for the critical temperatures and pressures as a function of the square-well range.     

Negative linear compressibility in confined dilatating systems

The role of a matrix response to a fluid insertion is analyzed in terms of a perturbation theory and Monte Carlo simulations applied to a hard sphere fluid in a slit of fluctuating density-dependent width. It is demonstrated that a coupling of the fluid-slit repulsion, spatial confinement, and the matrix dilatation acts as an effective fluid-fluid attraction, inducing a pseudocritical state with divergent linear compressibility and noncritical density fluctuations. An appropriate combination of the dilatation rate, fluid density, and the slit size leads to the fluid states with negative linear compressibility. It is shown that the switching from positive to negative compressibility is accompanied by an abrupt change in the packing mechanism.     

Effect of critical fluctuations on adsorption of van der Waals fluid in a spherical cavity

A recently proposed non-uniform fifth-order thermodynamic perturbation theory (TPT) is employed to investigate the adsorption of a hard core attractive Yukawa (HCAY) fluid in a spherical cavity. Extensive comparison with available simulation data indicate that the non-uniform fifth-order TPT is sufficiently reliable in calculating the density profiles of the HCAY fluid in the highly confining geometry, and generally is more accurate than a previous third-order?+?second-order perturbation density functional theory. The non-uniform fifth-order TPT is free from numerically solving an Ornstein–Zernike integral equation, and also free of any adjustable parameter; consequently, it can be applied to both supercritical and subcritical temperature regions. The non-uniform fifth-order TPT is employed to investigate critical adsorption of the HCYA fluid in a single spherical cavity – it is disclosed that the critical fluctuations near the critical point induce depletion adsorption – quantitative theoretical calculation on relationship between the critical depletion adsorption, parameters of coexistence bulk phase and the responsible external field is in agreement with qualitative physical analysis.