1
|
Munguía-Valadez J, Chávez-Rojo MA, Sambriski EJ, Moreno-Razo JA. The generalized continuous multiple step (GCMS) potential: model systems and benchmarks. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2022; 34:184002. [PMID: 35090143 DOI: 10.1088/1361-648x/ac4fe8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/05/2021] [Accepted: 01/28/2022] [Indexed: 06/14/2023]
Abstract
The generalized continuous multiple step (GCMS) potential is presented in this work. Its flexible form allows forrepulsiveand/orattractivecontributions to be encoded through adjustable energy and length scales. The GCMS interaction provides a continuous representation of square-well, square-shoulder potentials and their variants for implementation in computer simulations. A continuous and differentiable energy representation is required to derive forces in conventional simulation algorithms. Molecular dynamics simulations are of particular interest when considering the dynamic properties of a system. The GCMS potential can mimic other interactions with a judicious choice of parameters due to the versatile sigmoid form. In this study, our benchmarks for the GCMS representation include triangular, Yukawa, Franzese, and Lennard-Jones potentials. Comparisons made with published data on volumetric phase diagrams, liquid structure, and diffusivity from model systems are in excellent agreement.
Collapse
Affiliation(s)
- Jorge Munguía-Valadez
- Departamento de Física, Universidad Autónoma Metropolitana-Iztapalapa, Avenida San Rafael Atlixco No. 186, Colonia Vicentina, Delegación Iztapalapa, Mexico City 09340 Mexico
| | - Marco Antonio Chávez-Rojo
- Facultad de Ciencias Químicas, Universidad Autónoma de Chihuahua, Circuito Universitario s/n, Campus II, Chihuahua, Chihuahua 31125, Mexico
| | - Edward John Sambriski
- Department of Chemistry, Delaware Valley University, 700 East Butler Avenue, Doylestown, PA 18901 United States of America
| | - José Antonio Moreno-Razo
- Departamento de Física, Universidad Autónoma Metropolitana-Iztapalapa, Avenida San Rafael Atlixco No. 186, Colonia Vicentina, Delegación Iztapalapa, Mexico City 09340 Mexico
| |
Collapse
|
2
|
Santos A, Yuste SB, López de Haro M. Structural and thermodynamic properties of hard-sphere fluids. J Chem Phys 2020; 153:120901. [PMID: 33003724 DOI: 10.1063/5.0023903] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
This Perspective article provides an overview of some of our analytical approaches to the computation of the structural and thermodynamic properties of single-component and multicomponent hard-sphere fluids. For the structural properties, they yield a thermodynamically consistent formulation, thus improving and extending the known analytical results of the Percus-Yevick theory. Approximate expressions linking the equation of state of the single-component fluid to the one of the multicomponent mixtures are also discussed.
Collapse
Affiliation(s)
- Andrés Santos
- Departamento de Física and Instituto de Computación Científica Avanzada (ICCAEx), Universidad de Extremadura, Badajoz E-06006, Spain
| | - Santos B Yuste
- Departamento de Física and Instituto de Computación Científica Avanzada (ICCAEx), Universidad de Extremadura, Badajoz E-06006, Spain
| | - Mariano López de Haro
- Instituto de Energías Renovables, Universidad Nacional Autónoma de México (U.N.A.M.), Temixco, Morelos 62580, Mexico
| |
Collapse
|
3
|
Zerón IM, Vega C, Benavides AL. Continuous version of a square-well potential of variable range and its application in molecular dynamics simulations. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1481232] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Affiliation(s)
- I. M. Zerón
- División de Ciencias e Ingenierías, Universidad de Guanajuato, Guanajuato, México
| | - C. Vega
- Depto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040, Madrid, Spain
| | - A. L. Benavides
- División de Ciencias e Ingenierías, Universidad de Guanajuato, Guanajuato, México
| |
Collapse
|
4
|
Paganini IE, Pastorino C, Urrutia I. Structure, thermodynamic properties, and phase diagrams of few colloids confined in a spherical pore. J Chem Phys 2015; 142:244707. [PMID: 26133449 DOI: 10.1063/1.4923164] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
Abstract
We study a system of few colloids confined in a small spherical cavity with event driven molecular dynamics simulations in the canonical ensemble. The colloidal particles interact through a short range square-well potential that takes into account the basic elements of attraction and excluded-volume repulsion of the interaction among colloids. We analyze the structural and thermodynamic properties of this few-body confined system in the framework of inhomogeneous fluids theory. Pair correlation function and density profile are used to determine the structure and the spatial characteristics of the system. Pressure on the walls, internal energy, and surface quantities such as surface tension and adsorption are also analyzed for a wide range of densities and temperatures. We have characterized systems from 2 to 6 confined particles, identifying distinctive qualitative behavior over the thermodynamic plane T - ρ, in a few-particle equivalent to phase diagrams of macroscopic systems. Applying the extended law of corresponding states, the square well interaction is mapped to the Asakura-Oosawa model for colloid-polymer mixtures. We link explicitly the temperature of the confined square-well fluid to the equivalent packing fraction of polymers in the Asakura-Oosawa model. Using this approach, we study the confined system of few colloids in a colloid-polymer mixture.
Collapse
Affiliation(s)
- Iván E Paganini
- Departamento de Física de la Materia Condensada, Centro Atómico Constituyentes, CNEA, Av.Gral. Paz 1499, 1650 Pcia. de Buenos Aires, Argentina
| | - Claudio Pastorino
- Departamento de Física de la Materia Condensada, Centro Atómico Constituyentes, CNEA, Av.Gral. Paz 1499, 1650 Pcia. de Buenos Aires, Argentina
| | - Ignacio Urrutia
- Departamento de Física de la Materia Condensada, Centro Atómico Constituyentes, CNEA, Av.Gral. Paz 1499, 1650 Pcia. de Buenos Aires, Argentina
| |
Collapse
|
5
|
Fantoni R, Giacometti A, Santos A. Bridging and depletion mechanisms in colloid-colloid effective interactions: A reentrant phase diagram. J Chem Phys 2015; 142:224905. [DOI: 10.1063/1.4922263] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Affiliation(s)
- Riccardo Fantoni
- Dipartimento di Scienze Molecolari e Nanosistemi, Università Ca’ Foscari Venezia, Calle Larga S. Marta DD2137, I-30123 Venezia, Italy
| | - Achille Giacometti
- Departamento de Física and Instituto de Computación Científica Avanzada (ICCAEx), Universidad de Extremadura, E-06071 Badajoz, Spain
| | - Andrés Santos
- Departamento de Física and Instituto de Computación Científica Avanzada (ICCAEx), Universidad de Extremadura, E-06071 Badajoz, Spain
| |
Collapse
|
6
|
Santos A, Yuste SB, de Haro ML, Bárcenas M, Orea P. Structural properties of fluids interacting via piece-wise constant potentials with a hard core. J Chem Phys 2013; 139:074505. [DOI: 10.1063/1.4818601] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
|
7
|
Bakhti B, Müller G, Maass P. Interacting hard rods on a lattice: distribution of microstates and density functionals. J Chem Phys 2013; 139:054113. [PMID: 23927249 DOI: 10.1063/1.4816379] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We derive exact density functionals for systems of hard rods with first-neighbor interactions of arbitrary shape but limited range on a one-dimensional lattice. The size of all rods is the same integer unit of the lattice constant. The derivation, constructed from conditional probabilities in a Markov chain approach, yields the exact joint probability distribution for the positions of the rods as a functional of their density profile. For contact interaction ("sticky core model") between rods, we give a lattice fundamental measure form of the density functional and present explicit results for contact correlators, entropy, free energy, and chemical potential. Our treatment includes inhomogeneous couplings and external potentials.
Collapse
Affiliation(s)
- Benaoumeur Bakhti
- Fachbereich Physik, Universität Osnabrück, Barbarastrasse 7, 49076 Osnabrück, Germany
| | | | | |
Collapse
|
8
|
Piasecki J, Szymczak P, Kozak JJ. Stability of phases of a square-well fluid within superposition approximation. J Chem Phys 2013; 138:164506. [PMID: 23635156 DOI: 10.1063/1.4801329] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The analytic and numerical methods introduced previously to study the phase behavior of hard sphere fluids starting from the Yvon-Born-Green (YBG) equation under the Kirkwood superposition approximation (KSA) are adapted to the square-well fluid. We are able to show conclusively that the YBG equation under the KSA closure when applied to the square-well fluid: (i) predicts the existence of an absolute stability limit corresponding to freezing where undamped oscillations appear in the long-distance behavior of correlations, (ii) in accordance with earlier studies reveals the existence of a liquid-vapor transition by the appearance of a "near-critical region" where monotonically decaying correlations acquire very long range, although the system never loses stability.
Collapse
Affiliation(s)
- Jarosław Piasecki
- Institute of Theoretical Physics, Faculty of Physics, University of Warsaw, Warsaw, Poland
| | | | | |
Collapse
|
9
|
Valadez-Pérez NE, Benavides AL, Schöll-Paschinger E, Castañeda-Priego R. Phase behavior of colloids and proteins in aqueous suspensions: Theory and computer simulations. J Chem Phys 2012; 137:084905. [PMID: 22938263 DOI: 10.1063/1.4747193] [Citation(s) in RCA: 69] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Affiliation(s)
- Néstor E Valadez-Pérez
- División de Ciencias e Ingenierías, Campus León, Universidad de Guanajuato, Loma del Bosque 103, Lomas del Campestre, 37150 León, Guanajuato, Mexico
| | | | | | | |
Collapse
|
10
|
Pavlyukhin YT. Thermodynamic perturbation theory of simple liquids: Monte Carlo simulation of a hard sphere system and the Helmholtz free energy of SW fluids. J STRUCT CHEM+ 2012. [DOI: 10.1134/s0022476612030092] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
11
|
|
12
|
Rohrmann RD, Santos A. Multicomponent fluids of hard hyperspheres in odd dimensions. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 83:011201. [PMID: 21405686 DOI: 10.1103/physreve.83.011201] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/21/2010] [Indexed: 05/30/2023]
Abstract
Mixtures of hard hyperspheres in odd-space dimensionalities are studied with an analytical approximation method. This technique is based on the so-called rational function approximation and provides a procedure for evaluating equations of state, structure factors, radial distribution functions, and direct correlation functions of additive mixtures of hard hyperspheres with any number of components and in arbitrary odd-dimension space. The method gives the exact solution of the Ornstein-Zernike equation coupled with the Percus-Yevick closure, thus, extending the solution for hard-sphere mixtures [J. L. Lebowitz, Phys. Rev. 133, A895 (1964)] to arbitrary odd dimensions. Explicit evaluations for binary mixtures in five dimensions are performed. The results are compared with computer simulations, and a good agreement is found.
Collapse
Affiliation(s)
- René D Rohrmann
- Instituto de Ciencias Astronómicas, de la Tierra y del Espacio (ICATE-CONICET), San Juan, Argentina.
| | | |
Collapse
|
13
|
Espíndola-Heredia R, del Río F, Malijevsky A. Optimized equation of the state of the square-well fluid of variable range based on a fourth-order free-energy expansion. J Chem Phys 2009; 130:024509. [DOI: 10.1063/1.3054361] [Citation(s) in RCA: 59] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
14
|
Zhou S, Solana JR. Third-order thermodynamic perturbation theory for effective potentials that model complex fluids. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2008; 78:021503. [PMID: 18850837 DOI: 10.1103/physreve.78.021503] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/03/2008] [Revised: 06/23/2008] [Indexed: 05/26/2023]
Abstract
We have performed Monte Carlo simulations to obtain the thermodynamic properties of fluids with two kinds of hard-core plus attractive-tail or oscillatory potentials. One of them is the square-well potential with small well width. The other is a model potential with oscillatory and decaying tail. Both model potentials are suitable for modeling the effective potential arising in complex fluids and fluid mixtures with extremely-large-size asymmetry, as is the case of the solvent-induced depletion interactions in colloidal dispersions. For the former potential, the compressibility factor, the excess energy, the constant-volume excess heat capacity, and the chemical potential have been obtained. For the second model potential only the first two of these quantities have been obtained. The simulations cover the whole density range for the fluid phase and several temperatures. These simulation data have been used to test the performance of a third-order thermodynamic perturbation theory (TPT) recently developed by one of us [S. Zhou, Phys. Rev. E 74, 031119 (2006)] as compared with the well-known second-order TPT based on the macroscopic compressibility approximation due to Barker and Henderson. It is found that the first of these theories provides much better accuracy than the second one for all thermodynamic properties analyzed for the two effective potential models.
Collapse
Affiliation(s)
- Shiqi Zhou
- School of Physics Science and Technology, Central South University, Changsha, Hunan, 410083, China.
| | | |
Collapse
|
15
|
|
16
|
Rohrmann RD, Santos A. Structure of hard-hypersphere fluids in odd dimensions. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2007; 76:051202. [PMID: 18233644 DOI: 10.1103/physreve.76.051202] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/21/2007] [Indexed: 05/25/2023]
Abstract
The structural properties of single component fluids of hard hyperspheres in odd space dimensionalities d are studied with an analytical approximation method that generalizes the rational function approximation earlier introduced in the study of hard-sphere fluids [S. B. Yuste and A. Santos, Phys. Rev. A 43, 5418 (1991)]. The theory makes use of the exact form of the radial distribution function to first order in density and extends it to finite density by assuming a rational form for a function defined in Laplace space, the coefficients being determined by simple physical requirements. Fourier transform in terms of reverse Bessel polynomials constitute the mathematical framework of this approximation, from which an analytical expression for the static structure factor is obtained. In its most elementary form, the method recovers the solution of the Percus-Yevick closure to the Ornstein-Zernike equation for hyperspheres at odd dimensions. The present formalism allows one to go beyond by yielding solutions with thermodynamic consistency between the virial and compressibility routes to any desired equation of state. Excellent agreement with available computer simulation data at d=5 and d=7 is obtained.
Collapse
Affiliation(s)
- René D Rohrmann
- Departamento de Física, Universidad de Extremadura, E-06071 Badajoz, Spain.
| | | |
Collapse
|
17
|
Zhou S. Performance Evaluation of Third-Order Thermodynamic Perturbation Theory and Comparison with Existing Liquid State Theories. J Phys Chem B 2007; 111:10736-44. [PMID: 17713938 DOI: 10.1021/jp071044y] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
To evaluate the performance of a recently proposed third-order thermodynamic perturbation theory (TPT), we employ the third TPT for calculation of thermodynamic properties such as compressibility factor, internal energy, excess chemical potential, gas-liquid coexistence curve, and critical properties of several fluids. By comparing the third-order TPT results with corresponding simulation data available in literature and supplied in the present report and theoretical results from several other theoretical approaches, one concludes that the third-order TPT is, in general, more accurate than other approaches such as Barker-Henderson second-order TPT using a macroscopic compressibility approximation (MCA-TPT), self-consistent Ornstein-Zernike approach, Monte Carlo perturbation theory, and a specially devised equation of state. Specifically, the third-order TPT can predict quantitatively a double critical phenomena of gas-liquid transition and a low-density liquid (LDL)-high-density liquid (HDL) transition associated with a soft core (SC) potential fluid very satisfactorily, but the predictions for the LDL-HDL transition based on the second-order MCA-TPT are quantitatively very bad or qualitatively incorrect. The failure of the second-order MCA-TPT for the SC fluid can be ascribed to the facts that for the SC potential the second-order and third-order terms of the perturbation expansion are not small quantities and that the second-order term is underestimated by the MCA. It is concluded that the present third-order version of the TPT is reliable for varying model fluids.
Collapse
Affiliation(s)
- Shiqi Zhou
- School of Physics Science and Technology, Central South University, Changsha, Hunan, China.
| |
Collapse
|
18
|
Zhou S. Improvement on macroscopic compressibility approximation and beyond. J Chem Phys 2006; 125:144518. [PMID: 17042620 DOI: 10.1063/1.2353834] [Citation(s) in RCA: 32] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A numerical procedure is proposed to extend the thermodynamic perturbation expansion (TPE) to a higher order. It is shown that the present second order term is superior to that due to a macroscopic compressibility approximation (MCA), a local compressibility approximation, and a superposition approximation by Barker and Henderson [Rev. Mod. Phys. 48, 587 (1976)]. Extensive model calculation and comparison with simulation data available in literature and supplied in the present report indicate that the present third order TPE is superior to a previous second order TPE based on the MCA, two previous perturbation theories, which are respectively based on an analytical mean spherical approximation for an Ornstein-Zernike equation, and an assumed explicit functional form for the Laplace transform of radial distribution function multiplied by radial distance, and a recent generalized van der Waals theory. The present critical temperature for a hard core attractive Yukawa fluid of varying range is in very good agreement with that due to a hierarchical reference theory. The present third order TPE is computationally far more modest than the self-consistent integral equation theory, and therefore is a viable alternative to use of the latter.
Collapse
Affiliation(s)
- Shiqi Zhou
- Institute of Modern Statistical Mechanics, Hunan University of Technology, Wenhua Road, Zhuzhou 412008, People's Republic of China.
| |
Collapse
|
19
|
López-Rendón R, Reyes Y, Orea P. Thermodynamic properties of short-range square well fluid. J Chem Phys 2006; 125:084508. [PMID: 16965030 DOI: 10.1063/1.2338307] [Citation(s) in RCA: 48] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The interfacial properties of short-range square well fluid with lambda=1.15, 1.25, and 1.375 were determined by using single canonical Monte Carlo simulations. Simulations were carried out in the vapor-liquid region. The coexistence curves of these models were calculated and compared to those previously reported in the literature and good agreement was found among them. We found that the surface tension curves for any potential model of short range form a single master curve when we plot gamma* vs TT(c). It is demonstrated that the critical reduced second virial coefficient B(2)* as a function of interaction range or T(c)* is not constant.
Collapse
Affiliation(s)
- R López-Rendón
- Departamento de Química UAM-Iztapalapa, Avenida San Rafael Atlixco 186, 09340 México D.F., México
| | | | | |
Collapse
|
20
|
Malijevský A, Yuste SB, Santos A. How “sticky” are short-range square-well fluids? J Chem Phys 2006; 125:074507. [PMID: 16942351 DOI: 10.1063/1.2244549] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The aim of this work is to investigate to what extent the structural properties of a short-range square-well (SW) fluid of range lambda at a given packing fraction eta and reduced temperature T* = kBT/epsilon can be represented by those of a sticky-hard-sphere (SHS) fluid at the same packing fraction and an effective stickiness parameter tau(T*,lambda). Such an equivalence cannot hold for the radial distribution function g(r) since this function has a delta singularity at contact (r = sigma) in the SHS case, while it has a jump discontinuity at r = lambda sigma in the SW case. Therefore, the equivalence is explored with the cavity function y(r), i.e., we assume that [formula: see text]. Optimization of the agreement between y(SW) and y(SHS) to first order in density suggests the choice tau(T*,lambda) = [12(e(1/T* - 1)(lambda - 1)](-1). We have performed Monte Carlo (MC) simulations of the SW fluid for lambda = 1.05, 1.02, and 1.01 at several densities and temperatures T* such that tau(T*,lambda) = 0.13, 0.2, and 0.5. The resulting cavity functions have been compared with MC data of SHS fluids obtained by Miller and Frenkel[J. Phys.: Condens. Matter 16, S4901 (2004)]. Although, at given values of eta and tau, some local discrepancies between y(SW) and y(SHS) exist (especially for lambda = 1.05), the SW data converge smoothly toward the SHS values as lambda-1 decreases. In fact, precursors of the singularities of y(SHS) at certain distances due to geometrical arrangements are clearly observed in y(SW). The approximate mapping y(SW)-->y(SHS) is exploited to estimate the internal energy and structure factor of the SW fluid from those of the SHS fluid. Taking for y(SHS) the solution of the Percus-Yevick equation as well as the rational-function approximation, the radial distribution function g(r) of the SW fluid is theoretically estimated and a good agreement with our MC simulations is found. Finally, a similar study is carried out for short-range SW fluid mixtures.
Collapse
Affiliation(s)
- Alexandr Malijevský
- E. Hála Laboratory of Thermodynamics, Academy of Science of the Czech Republic, 16502 Prague 6, Czech Republic.
| | | | | |
Collapse
|
21
|
López de Haro M, Santos A, Yuste SB. On the radial distribution function of a hard-sphere fluid. J Chem Phys 2006; 124:236102. [PMID: 16821959 DOI: 10.1063/1.2201699] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
22
|
Zhou S. Isostructural solid–solid transitions in binary asymmetrical hard sphere system: Based on solvent-mediated potential. J Colloid Interface Sci 2005; 288:308-12. [PMID: 15927593 DOI: 10.1016/j.jcis.2005.02.070] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2005] [Revised: 02/16/2005] [Accepted: 02/18/2005] [Indexed: 11/27/2022]
Abstract
Size dependence is imparted onto a modified bridge functional, adopted for a recently proposed semi-analytical hard sphere reference system theory for calculation of solvent-mediated potential (SMP). The SMP for two large hard sphere particles immersed in a small hard sphere solvent bath predicted by the present improved version is in satisfactory agreement with the prediction from a theoretically based fitting formula. Isostructural solid-solid transitions in the binary asymmetrical hard sphere system are investigated based on a single-component macrofluid approximation combined with the improved version. It is found that the isostructural solid-solid transition appears when size asymmetry increases. The limiting asymmetry size ratio is near 1/8. As the size asymmetry increases, critical density for both large and small hard sphere components for the fcc isostructural solid-solid transition increases and decreases, respectively.
Collapse
Affiliation(s)
- Shiqi Zhou
- Research Institute of Modern Statistical Mechanics, Zhuzhou Institute of Technology, Zhuzhou City 412008, People's Republic of China.
| |
Collapse
|
23
|
Largo J, Solana JR, Yuste SB, Santos A. Pair correlation function of short-ranged square-well fluids. J Chem Phys 2005; 122:84510. [PMID: 15836066 DOI: 10.1063/1.1855312] [Citation(s) in RCA: 52] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We have performed extensive Monte Carlo simulations in the canonical (NVT) ensemble of the pair correlation function for square-well fluids with well widths lambda-1 ranging from 0.1 to 1.0, in units of the diameter sigma of the particles. For each one of these widths, several densities rho and temperatures T in the ranges 0.1< or =rhosigma(3)< or =0.8 and T(c)(lambda) less or approximately T less or approximately 3T(c)(lambda), where T(c)(lambda) is the critical temperature, have been considered. The simulation data are used to examine the performance of two analytical theories in predicting the structure of these fluids: the perturbation theory proposed by Tang and Lu [Y. Tang and B. C.-Y. Lu, J. Chem. Phys. 100, 3079 (1994); 100, 6665 (1994)] and the nonperturbative model proposed by two of us [S. B. Yuste and A. Santos, J. Chem. Phys. 101 2355 (1994)]. It is observed that both theories complement each other, as the latter theory works well for short ranges and/or moderate densities, while the former theory works for long ranges and high densities.
Collapse
Affiliation(s)
- J Largo
- Departamento de Física Aplicada, Universidad de Cantabria, E-39005 Santander, Spain.
| | | | | | | |
Collapse
|
24
|
Skibinsky A, Buldyrev SV, Franzese G, Malescio G, Stanley HE. Liquid-liquid phase transitions for soft-core attractive potentials. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2004; 69:061206. [PMID: 15244553 DOI: 10.1103/physreve.69.061206] [Citation(s) in RCA: 67] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/23/2003] [Revised: 02/02/2004] [Indexed: 05/24/2023]
Abstract
Using event-driven molecular dynamics simulations, we study a three-dimensional one-component system of spherical particles interacting via a discontinuous potential combining a repulsive square soft core and an attractive square well. In the case of a narrow attractive well, it has been shown that this potential has two metastable gas-liquid critical points. Here we systematically investigate how the changes of the parameters of this potential affect the phase diagram of the system. We find a broad range of potential parameters for which the system has both a gas-liquid critical point C1 and a liquid-liquid critical point C2. For the liquid-gas critical point we find that the derivatives of the critical temperature and pressure, with respect to the parameters of the potential, have the same signs: they are positive for increasing width of the attractive well and negative for increasing width and repulsive energy of the soft core. This result resembles the behavior of the liquid-gas critical point for standard liquids. In contrast, for the liquid-liquid critical point the critical pressure decreases as the critical temperature increases. As a consequence, the liquid-liquid critical point exists at positive pressures only in a finite range of parameters. We present a modified van der Waals equation which qualitatively reproduces the behavior of both critical points within some range of parameters, and gives us insight on the mechanisms ruling the dependence of the two critical points on the potential's parameters. The soft-core potential studied here resembles model potentials used for colloids, proteins, and potentials that have been related to liquid metals, raising an interesting possibility that a liquid-liquid phase transition may be present in some systems where it has not yet been observed.
Collapse
Affiliation(s)
- A Skibinsky
- Center for Polymer Studies and Department of Physics, Boston University, Boston, Massachusetts 02215, USA
| | | | | | | | | |
Collapse
|
25
|
|
26
|
Reiner A, Kahl G. The hierarchical reference theory as applied to square well fluids of variable range. J Chem Phys 2002. [DOI: 10.1063/1.1483258] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
27
|
|