Direct Pair Correlation Function Research Articles

Indirect Solution of Ornstein-Zernike Equation Using the Hopfield Neural Network Method

Microscopic information, such as the pair distribution and direct correlation functions, can be obtained from experimental data. From these correlation functions, thermodynamical quantities and the potential interaction function can be recovered. Derivations of Ornstein-Zernike equation and Hopfield Neural Network method are given first, as a theoretical background to follow the present work. From these two frameworks, structural information, such as the radial distribution (g(r)) and direct correlation (C(r)) functions, were retrieved from neutron scattering experimental data. The problem was solved considering simple initial conditions, which does not require any previous information about the system, making it clear the robustness of the Hopfield Neural Network method. The pair interaction potential was estimated in the Percus-Yevick (PY) and hypernetted chain (HNC) approximations and a poor agreement, compared with the Lennard-Jones 6-12 potential, was observed for both cases, suggesting the necessity of a more accurate closure relation to describe the system. In this sense, the Hopfield Neural Network together with experimental information provides an alternative approach to solve the Ornstein-Zernike equations, avoiding the limitations imposed by the closure relation.

The solvent mediated interaction potential between solute particles: theory and applications.

In this paper we develop a theory to calculate the solvent mediated interaction potential between solute particles dispersed in a solvent. The potential is a functional of the instantaneous distribution of solute particles and is expressed in terms of the solute-solvent direct pair correlation function and the density-density correlation function of the bulk solvent. The dependence of the direct pair correlation function on multi-point correlations of the solute distribution is simplified with a mean field approximation. A self consistent approach is developed to calculate the effective potential between solute particles, the solute-solvent and the solute-solute correlation functions. The significance of the solvent fluctuations on the range of the effective potential is elucidated. The theory is applied to calculate equilibrium properties of the Asakura-Oosawa (AO) model for several values of solute and solvent densities and for several values of the particles size ratio. The results give a quantitative description of many-body effect on the effective potential and on the pair correlation functions.

On the decay of the pair correlation function and the line of vanishing excess isothermal compressibility in simple fluids.

We revisit the competition between attractive and repulsive interparticle forces in simple fluids and how this governs and connects the macroscopic phase behavior and structural properties, as manifested in pair correlation functions. We focus on the asymptotic decay of the total correlation function h(r) which is, in turn, controlled by the form of the pair direct correlation function c(r). The decay of rh(r) to zero can be exponential (monotonic) if attraction dominates repulsion and exponentially damped oscillatory otherwise. The Fisher-Widom (FW) line separates the phase diagram into two regions characterized by the two different types of asymptotic decays. We show that there is a new and physically intuitive thermodynamic criterion which approximates well the actual FW line. This new criterion defines a line where the isothermal compressibility takes its ideal gas value χT=χT id. We test our hypothesis by considering four commonly used models for simple fluids. In all cases, the new criterion yields a line in the phase diagram that is close to the actual FW line for the thermodynamic state points that are most relevant. We also investigate (Widom) lines of maximal correlation length, emphasizing the importance of distinguishing between the true and Ornstein-Zernike correlation lengths.

The standard mean-field treatment of inter-particle attraction in classical DFT is better than one might expect.

In classical density functional theory (DFT), the part of the Helmholtz free energy functional arising from attractive inter-particle interactions is often treated in a mean-field or van der Waals approximation. On the face of it, this is a somewhat crude treatment as the resulting functional generates the simple random phase approximation (RPA) for the bulk fluid pair direct correlation function. We explain why using standard mean-field DFT to describe inhomogeneous fluid structure and thermodynamics is more accurate than one might expect based on this observation. By considering the pair correlation function g(x) and structure factor S(k) of a one-dimensional model fluid, for which exact results are available, we show that the mean-field DFT, employed within the test-particle procedure, yields results much superior to those from the RPA closure of the bulk Ornstein-Zernike equation. We argue that one should not judge the quality of a DFT based solely on the approximation it generates for the bulk pair direct correlation function.

Density-functional theory for fluid-solid and solid-solid phase transitions.

We develop a theory to describe solid-solid phase transitions. The density functional formalism of classical statistical mechanics is used to find an exact expression for the difference in the grand thermodynamic potentials of the two coexisting phases. The expression involves both the symmetry conserving and the symmetry broken parts of the direct pair correlation function. The theory is used to calculate phase diagram of systems of soft spheres interacting via inverse power potentials u(r)=ε(σ/r)^{n}, where parameter n measures softness of the potential. We find that for 1/n<0.154 systems freeze into the face centered cubic (fcc) structure while for 1/n≥0.154 the body-centred-cubic (bcc) structure is preferred. The bcc structure transforms into the fcc structure upon increasing the density. The calculated phase diagram is in good agreement with the one found from molecular simulations.

Analytic solution of the Ornstein-Zernike relation for inhomogeneous liquids.

The properties of a classical simple liquid are strongly affected by the application of an external potential that supports inhomogeneity. To understand the nature of these property changes, the equilibrium particle distribution functions of the liquid have, typically, been calculated directly using either integral equation or density functional based analyses. In this study, we develop a different approach with a focus on two distribution functions that characterize the inhomogeneous liquid: the pair direct correlation function c(r1,r2) and the pair correlation function g(r1,r2). With g(r1,r2) considered to be an experimental observable, we solve the Ornstein-Zernike equation for the inhomogeneous liquid to obtain c(r1,r2), using information about the well studied and resolved g(0)(r1,r2) and c(0)(r1,r2) for the parent homogeneous ((0)) system. In practical cases, where g(r1,r2) is available from experimental data in a discrete form, the resulting c(r1,r2) is expressed as an explicit function of g(r1,r2) in a discrete form. A weaker continuous form of solution is also obtained, in the form of an integral equation with finite integration limits. The result obtained with our formulation is tested against the exact solutions for the correlation and distribution functions of a one-dimensional inhomogeneous hard rod liquid. Following the success of that test, the formalism is extended to higher dimensional systems with explicit consideration of the two-dimensional liquid.

Fluid-solid transition in simple systems using density functional theory.

A free energy functional for a crystal which contains both the symmetry-conserved and symmetry-broken parts of the direct pair correlation function has been used to investigate the fluid-solid transition in systems interacting via purely repulsive Weeks-Chandler-Anderson Lennard-Jones potential and the full Lennard-Jones potential. The results found for freezing parameters for the fluid-face centred cubic crystal transition are in very good agreement with simulation results. It is shown that although the contribution made by the symmetry broken part to the grand thermodynamic potential at the freezing point is small compared to that of the symmetry conserving part, its role is crucial in stabilizing the crystalline structure and on values of the freezing parameters.

Structure of a planar electric double layer containing size-asymmetric ions: density functional approach

We have studied the structure of the electrolytes with asymmetries in charge and size near a charged planar electric double layer by a density functional theory. In the present theory, the hard-sphere contribution has been approximated as the direct pair correlation function with the coupling parameter, whereas the electronic contribution has been approximated as the mean-spherical approximation in the bulk phase. This theoretical approach for the size-symmetric and size-asymmetric electrolytes displays a good agreement with the simulation results over a wide range of surface charge densities and electrolyte concentrations. However, the accuracy between the present theory and the simulation results slightly deteriorates for the highly size-asymmetric electrolytes and the multivalent electrolytes. In these cases, the performance of the present theory is comparable to those of the simplified extension of the Poisson–Boltzmann theory and the modified Poisson–Boltzmann theory. The calculated result indicates that the surface charge distribution function, which was introduced as an indicator for studying the charge reversal, layering effect, and surface charge amplification in a planar electric double layer, describes the electronic properties of a planar electric double layer well.

Communication: Integral equation theory for pair correlation functions in a crystal.

A method for calculating pair correlation functions in a crystal is developed. The method is based on separating the one- and two-particle correlation functions into the symmetry conserving and the symmetry broken parts. The conserving parts are calculated using the integral equation theory of homogeneous fluids. The symmetry broken part of the direct pair correlation function is calculated from a series written in powers of order parameters and that of the total pair correlation function from the Ornstein-Zernike equation. The results found for a two-dimensional hexagonal lattice show that the method provides accurate and detailed informations about the pair correlation functions in a crystal.

Electric double layer for a size-asymmetric electrolyte around a spherical colloid

We have studied the structure of a size-asymmetric electrolyte on charged colloids by a density functional perturbation theory. The hard-sphere contribution has been approximated as the direct pair correlation function with the coupling parameter, whereas the electronic contribution has been approximated as the mean-spherical approximation in the bulk phase. The calculated results for the ionic density distributions and mean electrostatic potentials are in very good agreement with the computer simulations over a wide range of colloid sizes and electrolyte concentrations. The present theory provides better structural results than the hypernetted-chain equation based on the mean spherical approximation. We have confirmed that the overcharging appears when the counterions are larger than the coions. The overcharging disappears everywhere when the electrostatic repulsion becomes strong enough, while the charge reversal is observed when the coions are larger than the counterions, and the reversal effect appears for a size-asymmetric electrolyte at high surface charge densities. The charge reversal occurs even for the point of zero charge, mainly due to the depletion force between two ions. The present theory is able to provide interesting insights about the charge reversal and overcharging phenomena occurring at the interface.

Correlation functions in liquids and crystals: Free-energy functional and liquid-to-crystal transition

A free-energy functional for a crystal that contains both the symmetry-conserved and symmetry-broken parts of the direct pair-correlation function has been used to investigate the crystallization of fluids in three dimensions. The symmetry-broken part of the direct pair-correlation function has been calculated using a series in ascending powers of the order parameters and which contains three- and higher-body direct correlation functions of the isotropic phase. It is shown that a very accurate description of freezing transitions for a wide class of potentials is found by considering the first two terms of this series. The results found for freezing parameters including the structure of the frozen phase for fluids interacting via the inverse power potential u(r)=ε(σ/r)(n) for n ranging from 4 to ∞ are in very good agreement with simulation results. It is found that for n>6.5 the fluid freezes into a face-centered cubic (fcc) structure while for n≤6 the body-centered cubic (bcc) structure is preferred. The fluid-bcc-fcc triple point is found to be at 1/n=0.158, which is in good agreement with simulation result.

A Site Density Functional Theory for Water: Application to Solvation of Amino Acid Side Chains.

We report a site density functional theory (SDFT) based on the conventional atomistic models of water and the universality ansatz of the bridge functional. The excess Helmholtz energy functional is formulated in terms of a quadratic expansion with respect to the local density deviation from that of a uniform system and a universal functional for all higher-order terms approximated by that of a reference hard-sphere system. With the atomistic pair direct correlation functions of the uniform system calculated from MD simulation and an analytical expression for the bridge functional from the modified fundamental measure theory, the SDFT can be used to predict the structure and thermodynamic properties of water under inhomogeneous conditions with a computational cost negligible in comparison to that of brute-force simulations. The numerical performance of the SDFT has been demonstrated with the predictions of the solvation free energies of 15 molecular analogs of amino acid side chains in water represented by SPC/E, SPC, and TIP3P models. For theTIP3P model, a comparison of the theoretical predictions with MD simulation and experimental data shows agreement within 0.64 and 1.09 kcal/mol on average, respectively.

Freezing of a two-dimensional fluid into a crystalline phase: Density functional approach

A free-energy functional for a crystal proposed by Singh and Singh [Europhys. Lett. 88, 16005 (2009)] which contains both the symmetry conserved and symmetry broken parts of the direct pair correlation function has been used to investigate the crystallization of a two-dimensional fluid. The results found for fluids interacting via the inverse power potential u(r)=ε(σ/r)(n) for n=3,6, and 12 are in good agreement with experimental and simulation results. The contribution made by the symmetry broken part to the grand thermodynamic potential at the freezing point is found to increase with the softness of the potential. Our results explain why the Ramakrishnan-Yussouff [Phys. Rev. B 19, 2775 (1979)] free-energy functional gave good account of freezing transitions of hard-core potentials but failed for potentials that have soft core and/or attractive tail.

Spatial inhomogeneities in ionic liquids, charged proteins, and charge stabilized colloids from collective variables theory

Effects of size and charge asymmetry between oppositely charged ions or particles on spatial inhomogeneities are studied for a large range of charge and size ratios. We perform a stability analysis of the primitive model of ionic systems with respect to periodic ordering using the collective variables-based theory. We extend previous studies [Ciach et al., Phys. Rev. E 75, 051505 (2007)] in several ways. First, we employ a nonlocal approximation for the reference hard-sphere fluid which leads to the Percus-Yevick pair direct correlation functions for the uniform case. Second, we use the Weeks-Chandler-Anderson regularization scheme for the Coulomb potential inside the hard core. We determine the relevant order parameter connected with the periodic ordering and analyze the character of the dominant fluctuations along the λ lines. We show that the above-mentioned modifications produce large quantitative and partly qualitative changes in the phase diagrams obtained previously. We discuss possible scenarios of the periodic ordering for the whole range of size and charge ratios of the two ionic species, covering electrolytes, ionic liquids, charged globular proteins or nanoparticles in aqueous solutions, and charge-stabilized colloids.

Weighted correlation approach: An extended version with applications to the hard-sphere fluid

The purpose of this study is to extend the weighted correlation approach (WCA) for inhomogeneous fluids. It now introduces a generic expression to evaluate the single-particle direct correlation function in terms of a series of pair direct correlation functions weighted by different correlation-weight functions and adjustable weight factors. When applied for practical use, however, approximations of the pair direct correlation functions have to be made, together with appropriate definitions of the weighted densities and the choices of the correlation-weight functions. The WCA approach would, then, not only help us to connect and compare different strategies and their underlying assumptions in the density functional approaches, but also enable us to propose and apply density functional theory methods to predict the density profile of, e.g., the hard-sphere fluid confined between a pair of parallel planar hard walls. Numerical results of the extended WCA approach, against the Monte Carlo (MC) simulations in a range of surface separations and bulk densities, suggest that it is capable of representing the fine features of the hard-sphere density distributions. The WCA results also agree well with the calculations from the fundamental measure theory. In addition, the thermodynamic self-consistency of the WCA approach is confirmed by its fairly good agreement with the MC fitted data for the surface tension of a hard-sphere fluid at a planar hard wall. All these tests show that a pure WCA approach can be constructed to investigate the states of ionic hard-sphere fluids.

Free-energy functional for freezing transitions: Hard-sphere systems freezing into crystalline and amorphous structures

A free-energy functional that contains both the symmetry-conserved and symmetry-broken parts of the direct pair correlation function has been used to investigate the freezing of a system of hard spheres into crystalline and amorphous structures. The freezing parameters for fluid-crystal transition have been found to be in very good agreement with the results found from simulations. We considered amorphous structures found from molecular dynamics simulations at packing fractions η lower than the glass close packing fraction η(J) and investigated their stability compared to that of a homogeneous fluid. The existence of a free-energy minimum corresponding to a density distribution of overlapping Gaussians centered around an amorphous lattice depicts a deeply supercooled state with a heterogeneous density profile.

Density functional theory of liquid crystals and surface anchoring: Hard Gaussian overlap-sphere and hard Gaussian overlap-surface potentials

This article applies the density functional theory to confined liquid crystals, comprised of ellipsoidal shaped particles interacting through the hard Gaussian overlap (HGO) potential. The extended restricted orientation model proposed by Moradi and co-workers [J. Phys.: Condens. Matter 17, 5625 (2005)] is used to study the surface anchoring. The excess free energy is calculated as a functional expansion of density around a reference homogeneous fluid. The pair direct correlation function (DCF) of a homogeneous HGO fluid is approximated, based on the optimized sum of Percus-Yevick and Roth DCF for hard spheres; the anisotropy introduced by means of the closest approach parameter, the expression proposed by Marko [Physica B 392, 242 (2007)] for DCF of HGO, and hard ellipsoids were used. In this study we extend an our previous work [Phys. Rev. E 72, 061706 (2005)] on the anchoring behavior of hard particle liquid crystal model, by studying the effect of changing the particle-substrate contact function instead of hard needle-wall potentials. We use the two particle-surface potentials: the HGO-sphere and the HGO-surface potentials. The average number density and order parameter profiles of a confined HGO fluid are obtained using the two particle-wall potentials. For bulk isotropic liquid, the results are in agreement with the Monte Carlo simulation of Barmes and Cleaver [Phys. Rev. E 71, 021705 (2005)]. Also, for the bulk nematic phase, the theory gives the correct density profile and order parameter between the walls.

Density functional theory approach to the freezing transition in two-dimensional colloid system

The freezing transition of two-dimensional colloid system is studied by the density functional theory (DFT). The liquid direct-pair correlation function (DPCF) is simulated by molecular dynamics simulations. With the DPCF, the freezing transition is determined using DFT with the Fourier components of the density ρ K as the order parameters. The DFT results are consistent with the simulation results if ρ Ks of the lowest three shells of reciprocal lattice vector (RLV) are used in the DFT calculation. The first order freezing transition takes place at DPCF ρ0 C(ka) = 0.788 on the first RLV shell. The laser induced freezing of colloid system is also studied using the DFT.

Interfacial and wetting properties of a binary point Yukawa fluid

We investigate the interfacial phase behavior of a binary fluid mixture composed of repulsive point Yukawa particles. Using a simple approximation for the Helmholtz free energy functional, which yields the random phase approximation for the pair direct correlation functions, we calculate the equilibrium fluid density profiles of the two species of particles adsorbed at a planar wall. We show that for a particular choice (repulsive exponential) of the wall potentials and the fluid pair-potential parameters, the Euler-Lagrange equations for the equilibrium fluid density profiles may be transformed into a single ordinary differential equation and the profiles obtained by a simple quadrature. For certain other choices of the fluid pair-potential parameters fluid-fluid phase separation of the bulk fluid is observed. We find that when such a mixture is exposed to a planar hard wall, the fluid exhibits complete wetting on the species 2 poor side of the binodal, i.e., we observe a thick film of fluid rich in species 2 adsorbed at the hard wall. The thickness of the wetting film grows logarithmically with the concentration difference between the fluid state point and the binodal and is proportional to the bulk correlation length of the intruding (wetting) fluid phase. However, for state points on the binodal that are further from the critical point, we find there is no thick wetting film. We determine the accompanying line of first-order (prewetting) surface phase transitions which separate a thin and thick adsorbed film. We show that for some other choices of repulsive wall potentials the prewetting line is still present, but its location and extent in the phase diagram is strongly dependent on the wall-fluid interaction parameters.

On spatial variations of nematic ordering

We present a theory of orientational order in nematic liquid crystals which interpolates between several distinct approaches based on the director field (Oseen and Frank), order parameter tensor (Landau and de Gennes), and orientation probability density function (Onsager). As in density-functional theories, the suggested free energy is a functional of spatially-dependent orientation distribution, however, the nonlocal effects are taken into account via phenomenological elastic terms rather than by means of a direct pair-correlation function. In illustration of this approach we consider a simplified model with orientation parameter on a circle and reveal its relation to the complex Ginzburg–Landau theory.