Microscopic structure and thermodynamics of a core-softened model fluid from the second-order integral equations theory

We have studied the structure and thermodynamic properties of isotropic three dimensional core-softened fluid by using the second-order Ornstein-Zernike integral equations completed by the hypernetted chain and Percus-Yevick closures. The radial distribution functions are compared with those from si...

Повний опис

Збережено в:

Бібліографічні деталі
Дата:	2011
Автори:	Pizio, O., Sokołowska, Z., Sokołowski, S.
Формат:	Стаття
Мова:	English
Опубліковано:	Інститут фізики конденсованих систем НАН України 2011
Назва видання:	Condensed Matter Physics
Онлайн доступ:	http://dspace.nbuv.gov.ua/handle/123456789/119939
Теги:	Додати тег Немає тегів, Будьте першим, хто поставить тег для цього запису!
Назва журналу:	Digital Library of Periodicals of National Academy of Sciences of Ukraine
Цитувати:	Microscopic structure and thermodynamics of a core-softened model fluid from the second-order integral equations theory / O. Pizio, Z. Sokołowska, S. Sokołowski // Condensed Matter Physics. — 2011. — Т. 14, № 1. — С. 13601: 1–12. — Бібліогр.: 62 назв. — англ.

Репозитарії

Digital Library of Periodicals of National Academy of Sciences of Ukraine

id	irk-123456789-119939
record_format	dspace
spelling	irk-123456789-1199392017-06-11T03:04:21Z Microscopic structure and thermodynamics of a core-softened model fluid from the second-order integral equations theory Pizio, O. Sokołowska, Z. Sokołowski, S. We have studied the structure and thermodynamic properties of isotropic three dimensional core-softened fluid by using the second-order Ornstein-Zernike integral equations completed by the hypernetted chain and Percus-Yevick closures. The radial distribution functions are compared with those from singlet integral equations and with computer simulation data. The limits of the region of density anomaly resulting from different approximate theories are established. The obtained results show that the second-order hypernetted chain approximation can be used to describe both the structure and the density anomaly of this model fluid. Moreover, we present the results of calculations of the bridge functions. Ми дослiдили структурнi i термодинамiчнi властивостi однорiдного тривимiрного плину з пом’якшеним кором, використовуючи iнтегральнi рiвняння Орнштейна - Цернiке другого порядку i згiперланцюговим замиканням та замиканням Перкуса - Євiка. Зроблено порiвняння радiальних функцiй розподiлу з вiдповiдними функцiями, отриманими iз синглетнихiнтегральних рiвнянь, а також з даними комп’ютерного моделювання. З рiзних наближених теорiй встановлено границi областi аномалiї густини. Отриманi результати показують, що гiперланцюгове наближення iнтегральних рiвнянь другого порядку може бути використане для опису як структури, так i аномалiї густини цього модельного плину. Крiм того, ми представляємо результати обчислень мiсткових функцiй. 2011 Article Microscopic structure and thermodynamics of a core-softened model fluid from the second-order integral equations theory / O. Pizio, Z. Sokołowska, S. Sokołowski // Condensed Matter Physics. — 2011. — Т. 14, № 1. — С. 13601: 1–12. — Бібліогр.: 62 назв. — англ. 1607-324X PACS: 61.20.Ja, 61.25.-f, 64.70.pm, 65.20.Jk DOI:10.5488/CMP.14.13601 arXiv:1106.3231 http://dspace.nbuv.gov.ua/handle/123456789/119939 en Condensed Matter Physics Інститут фізики конденсованих систем НАН України
institution	Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection	DSpace DC
language	English
description	We have studied the structure and thermodynamic properties of isotropic three dimensional core-softened fluid by using the second-order Ornstein-Zernike integral equations completed by the hypernetted chain and Percus-Yevick closures. The radial distribution functions are compared with those from singlet integral equations and with computer simulation data. The limits of the region of density anomaly resulting from different approximate theories are established. The obtained results show that the second-order hypernetted chain approximation can be used to describe both the structure and the density anomaly of this model fluid. Moreover, we present the results of calculations of the bridge functions.
format	Article
author	Pizio, O. Sokołowska, Z. Sokołowski, S.
spellingShingle	Pizio, O. Sokołowska, Z. Sokołowski, S. Microscopic structure and thermodynamics of a core-softened model fluid from the second-order integral equations theory Condensed Matter Physics
author_facet	Pizio, O. Sokołowska, Z. Sokołowski, S.
author_sort	Pizio, O.
title	Microscopic structure and thermodynamics of a core-softened model fluid from the second-order integral equations theory
title_short	Microscopic structure and thermodynamics of a core-softened model fluid from the second-order integral equations theory
title_full	Microscopic structure and thermodynamics of a core-softened model fluid from the second-order integral equations theory
title_fullStr	Microscopic structure and thermodynamics of a core-softened model fluid from the second-order integral equations theory
title_full_unstemmed	Microscopic structure and thermodynamics of a core-softened model fluid from the second-order integral equations theory
title_sort	microscopic structure and thermodynamics of a core-softened model fluid from the second-order integral equations theory
publisher	Інститут фізики конденсованих систем НАН України
publishDate	2011
url	http://dspace.nbuv.gov.ua/handle/123456789/119939
citation_txt	Microscopic structure and thermodynamics of a core-softened model fluid from the second-order integral equations theory / O. Pizio, Z. Sokołowska, S. Sokołowski // Condensed Matter Physics. — 2011. — Т. 14, № 1. — С. 13601: 1–12. — Бібліогр.: 62 назв. — англ.
series	Condensed Matter Physics
work_keys_str_mv	AT pizioo microscopicstructureandthermodynamicsofacoresoftenedmodelfluidfromthesecondorderintegralequationstheory AT sokołowskaz microscopicstructureandthermodynamicsofacoresoftenedmodelfluidfromthesecondorderintegralequationstheory AT sokołowskis microscopicstructureandthermodynamicsofacoresoftenedmodelfluidfromthesecondorderintegralequationstheory
first_indexed	2025-07-08T16:56:39Z
last_indexed	2025-07-08T16:56:39Z
_version_	1837098651354136576
fulltext	Condensed Matter Physics 2011, Vol. 14, No 1, 13601: 1–12 DOI:10.5488/CMP.14.13601 http://www.icmp.lviv.ua/journal Microscopic structure and thermodynamics of a core-softened model fluid from the second-order integral equations theory O. Pizio1, Z. Sokołowska2, S. Sokołowski3 1 Instituto de Quı́mica de la UNAM, Coyoacán 04510, México 2 Department of Physical Chemistry of Porous Materials, Institute of Agrophysics Polish Academy of Sciences, Doświadczalna 4, 20–290 Lublin, Poland 3 Department for the Modeling of Physico-Chemical Processes, Maria Curie-Skłodowska University, 20–031 Lublin, Poland Received April 18, 2010, in final form July 27, 2010 We have studied the structure and thermodynamic properties of isotropic three-dimensional core-softened fluid by using the second-order Ornstein-Zernike integral equations completed by the hypernetted chain and Percus-Yevick closures. The radial distribution functions are compared with those from singlet integral equa- tions and with computer simulation data. The limits of the region of density anomaly resulting from different approximate theories are established. The obtained results show that the second-order hypernetted chain ap- proximation can be used to describe both the structure and the density anomaly of this model fluid. Moreover, we present the results of calculations of the bridge functions. Key words: liquid theory, second-order integral equations, density anomaly, bridge functions PACS: 61.20.Ja, 61.25.-f, 64.70.pm, 65.20.Jk 1. Introduction Recently, significant research activity has been focused on the so-called core-softened models in which the repulsive part of the intermolecular interaction potential exhibits a softening region at “intermediate” interparticle separations, in addition to usual hard or soft (as in the case of Lennard-Jones (12,6) potential) repulsive branch at short distances. The softening region can be described by a linear or nonlinear ramp, a shoulder, a single or multiple attractive well, or a combination of all these features [1–11]. Model fluids with spherically symmetric core-softened potentials exhibit anomalous thermodynamic and dynamic behaviours [12–14] that usually occur in real fluids with directional interparticle interactions, e.g. water, silica, some liquid metals and phosphorus. Examples of the anomalous behavior include the existence of a density maximum as a function of temperature, an increase of the diffusion coefficient upon compression, and, for some systems, the existence of multiple fluid-fluid transitions. The reason why the soft-core model potentials can reproduce the density anomaly is that at low densities and temperatures, the neighboring molecules are separated by distances coinciding with the range of attractive potential well to minimize the energy. At higher temperatures, the particles can penetrate into the energetically less favorable softened core to gain more entropy, thus giving rise to an anomalous contraction upon heating, see e.g., [15, 16]. The systems with core-softened potential have been studied using computer simulations [4–21] as well as integral equations methods [6–9, 15, 18–20]. Among different models of the core-softened potentials, the potential used by Barros de Oliveira et al. [22–24] that yields both the liquid-gas phase transition and thermodynamic, structural and dynamic anomalies seems to be interesting. The potential is the sum of a Lennard-Jones potential c© O. Pizio, Z. Sokołowska, S. Sokołowski 13601-1 http://dx.doi.org/10.5488/CMP.14.13601 http://www.icmp.lviv.ua/journal O. Pizio, Z. Sokołowska, S. Sokołowski of the well depth εf and a Gaussian barrier centered on radius r = r0 with the height of aεf , u(r) = 4εf [ (σf r )12 − (σf r )6 ] + aεf exp [ − 1 c2 ( r − r0 σf )2 ] . (1) Depending on the choice of parameters, equation (1) represents the whole family of two length scales intermolecular interactions, from a deep double-well potential [3, 25, 26] to a repulsive shoulder [4, 27, 28]. Moreover, due to its continuity, the potential (1) is convenient both for the Monte Carlo and molecular dynamics simulations. For specific choices of the parameters, the potential (1) can exhibit a double well behavior, similar to the potential studied by Cho et al. [3, 25]. However, the attractive double well may bring both the liquid-gas phase transition and the anomalies to higher temperatures into an unstable region of the phase diagram in the pressure-temperature plane [26]. In order to circumvent this behavior, several studies [22–24, 29–32] were carried out for a specific set of parameters that give rise to a potential with a repulsive ramp followed by a small attractive well. In particular, the following values were often used: a = 5, c = 1 and r0 = 0.7σf . It was shown [22–24] that for the above set of parameters, the pressure-temperature curves have a minimum if calculated at certain densities (i.e. (∂P/∂T )ρ = 0). Consequently, the derivative (∂ρ/∂T )P = 0 for some thermodynamic states. The region of density anomaly corresponds to state points for which (∂ρ/∂T )P > 0 and is bounded by the locus of points for which the thermal expansion coefficient is equal to zero. In our recent work [31] we applied grand canonical Monte Carlo simulation and integral equa- tions with hypernetted chain (HNC), as well as Rogers-Young closures to study the above defined model. We found still another anomaly in the system, namely we demonstrated that for some tem- peratures the derivative of the density with respect to the chemical potential exhibits a minimum at a certain density, followed by a maximum. Also, a peculiarity in the dependence of the specific heat upon density was found. We compared the pair distribution functions resulting from integral equations and from computer simulations and established that common HNC approximation was not successful in capturing the region of anomalies in contrast to Rogers-Young approximation that imposes thermodynamic self-consistency by construction. However, the HNC was very accurate at high fluid densities. It is worth mentioning that a popular mean spherical approximation has not been used to study the model in question because this approach intrinsically requires an adequate splitting of the potential into a short- and long-ranged parts. This approximation does not seem to be beneficial for some soft-core models with two length scales and of this study in particular. The results of integral equations reported in the literature for the potential (1) were based on the singlet OZ equation [33]. A more sophisticated approach, but one that is also more demanding of computational resources, is based on the inhomogeneous OZ equation [34]. The inhomogeneous OZ equation has usually been used not only to calculate the structure and thermodynamic properties of fluids in contact with surfaces [35–39], but also to determine the structure and the surface tension at the liquid-vapor interface [40]. However, one can also assume that for a single-component fluid, the source of an external potential field is just a single distinguished particle identical to all remaining molecules of the system. This is the so-called Percus’ trick [41, 42] or the “test particle” method. The method was extended later by Attard [43, 44] in the framework of the second-order OZ approach for bulk fluids, see e.g. [45–47] for the discussion of related important issues. Moreover, this approach was successfully applied to several model fluids with spherically symmetric intermolecular potentials [48–54]. Better description of the desired properties or in some cases even qualitatively new findings in comparison to the singlet level theory were obtained. However, to our best knowledge, this approach has never been tested for the core-softened models of the interparticle potentials. Therefore, the motivation of this work is to use the second-order integral equations to the fluid of particles interacting via the core-softened potential (1). The second-order equations are applied to the study of the microscopic structure and thermodynamic properties of the fluid, in particular, in view of its anomalous behavior in a certain region of thermodynamic states. The results are compared with computer simulation data and with the predictions of the singlet integral equations [31]. We explore here the second-order Percus-Yevick and hypernetted closures to the nonuniform OZ equation. 13601-2 Core-softened model fluid 2. Theory The OZ equation for an inhomogeneous fluid, wherein the density ρ(r) is not constant reads h2(r1, r2) = c2(r1, r2) + ∫ dr3h2(r1, r3)c2(r3, r2)ρ(r3), (2) where h2(r1, r2) and c2(r1, r2) are the total and direct correlation functions of an nonuniform fluid. We have used here the subscript “2” in order to distinguish these functions from the common uniform fluid correlation functions. Generally, equation (2) applies when the inhomogeneity is due to an external potential field. However, one distinguished molecule can be also considered as a source of the external potential. In such cases equation (2) can be solved using any of the common approximations, as HNC h2(r1, r2) = exp[h2(r1, r2)− c2(r1, r2)− βu(\|r1 − r2\|)]− 1, (3) or Percus-Yevick (PY) approximation c2(r1, r2) = {1− exp[βu(\|r1 − r2\|)]}[h2(r1, r2) + 1]. (4) In the above u(\|r1−r2\|) is the pair potential and β = 1/kT is the inverse temperature. However, we stress that the above closures do not relate uniform, but rather nonuniform correlation functions and therefore we refer to the results of equation (2) with the closure given by equation (3), or by equation (4), as to the HNC2 and PY2 results, respectively, in contrast to the usual HNC and PY theories [33]. To have a set of equations completed we also need a relation between ρ(r) and the pair correla- tion functions. The exact equation developed by Lovett, Mou, Buff and Wertheim (LMBW) [55, 56] reads ∇ ln[y(r1)] = ∫ dr2c2(r1, r2)∇ρ(r2), (5) where y(r1) is the one-particle background (cavity) function, ρ(r1) = exp[−βv(r1)]y(r1) and v(r1) is the potential due to the particle which is regarded as the source of the inhomogeneity. In our case, v(r1) is just the pair potential, v(r1) ≡ u(\|r1 − 0\|), where the distinguished particle is set at the origin, r2 = 0. One should mention here an important difference between the theories used to study systems in contact with a surface [35–39] (or the systems involving a gas-liquid interface [40]) and theories based on the method of Attard. In the former case the two-particle correlation functions reduce to the functions for bulk fluids, providing the both particles to be located far away from an inhomogeneous region. In the Attard’s approach, the function h2 corresponds rather to the third- order (conditional) correlation function of a bulk fluid, while the “usual” pair correlation function is related to the local density. Indeed, because the bulk pair distribution function, g(r), gives the probability density of finding a pair of fluid molecules at a separation, r, the g(r) is related to the local density, ρ(r), via g(r) = ρ(\|r\|)/ρ, (6) where \|r\| is the distance from the distinguished particles that is the source of the system inhomo- geneity and ρ = lim\|r\|→∞ ρ(\|r\|) is the bulk density. Of course, the cavity function y(r) that enters equation (5) satisfies the relation g(r) = exp[−βu(r)]y(r). To solve equations (2) and (5) with the closure given by equation (3) or equation (4) we use the numerical algorithm Attard’s [43] that relies on the expansion of the two-particle functions (h2 and c2) in series of Legendre polynomials, for details see [43]. Note that the closure equations (3) and (4) constitute an approximation, similar to the common bulk theory, whereas the OZ equation (2) and the equation for the profile, (5) are exact. Equation (5) satisfies the relation g(r) = exp[−βu(r)]y(r). The bulk direct correlation function, c(r), can be recovered from g(r) using the bulk Ornstein-Zernike (OZ) equation h(\|r12\|)− c(\|r12\|) = ρ ∫ dr3h(\|r13\|)c(\|r23\|), (7) 13601-3 O. Pizio, Z. Sokołowska, S. Sokołowski where h(r) = g(r) − 1. The knowledge of the total, h(r), and the cavity, y(r), functions enable us to determines the bridge function, B(r). The bridge function is defined as, see e.g. [33, 57] B(r) = ln y(r)− γ(r), (8) where γ(r) = h(r)− c(r). The bridge function plays a key role in theories based on the singlet OZ equation (7). It is given by the smallest set of diagrams in the graphical expansion, and all other functions can be calculated from it via the OZ equation [33]. The main difference between the first- and the second-order integral equation theories for bulk fluids is that in the latter approach, the bridge function results from the calculations performed, whereas in the case of first-order theories the bridge function is imposed as a closure. Previous calculations carried out for hard spheres, as well as for Lennard-Jones fluids, gave insight into the course of the bridge functions. However, to our best knowledge no results of calculations of the bridge function for core-softened potential models have been presented so far. 3. Results and discussion The second-order integral equations were solved using the Attard’s algorithm [43]. Usually, 80 Legendre polynomials were used and the grid in r was 0.04σ. To test the accuracy we increased the number of Legendre polynomials to 120 and decreased the grid size to 0.025σ. Moreover, we also solved the singlet integral equations and carried out grand canonical ensemble Monte Carlo simulations. The details on the two latter methods are presented in our previous work [31]. We also note that the pressures given below were calculated from the virial equation of state (see equation (6) of [23]). Similarly to our previous work[31], all the calculations have been carried out assuming that the parameters of the potential (1) are a = 5, c = 1 and r0 = 0.7σf . First, we concentrate on a comparison of the pair distribution functions obtained from the second-order OZ equation, the singlet OZ equation and the Grand Canonical ensemble (GCMC) computer simulations [31]. Comparisons are for thermodynamic conditions that cover the density anomaly region found in molecular dynamics studies [23]. In figures 1 (a)–(c) we display the functions g(r) evaluated at the density ρ∗ = ρσ3 f = 0.1 and at three temperatures, T ∗ = kT/εf = 0.5, 0.3 and 0.2. At the highest considered temperature the results of the singlet theories (PY, HNC and Rogers-Young, (RY)), the second-order theories, PY2 and HNC2 are compared with GCMC data. Before discussing the results shown in figure 1, we recall that the singlet RY closure contains one adjustable parameter and this parameter is evaluated requiring the equality of pressures from the virial and compressibility routes, for details see [31]. In figure 1 (a) we see that the singlet PY approximation yields the radial distribution function that significantly deviates from all remaining results. The second-order PY2 approximation is def- initely more accurate than the singlet PY theory, but its results still differ much from computer simulations. In contrast to the singlet PY theory, the singlet HNC closure yields the results that qualitatively agree with computer simulations, but among all the singlet theories the thermody- namically self-consistent RY closure gives the best predictions. For the thermodynamic state from figure 1 (a) the second-order hypernetted chain approximation reproduces the GCMC data with the highest precision. We stress that in contrast to the singlet RY theory the HNC2 approach contains no adjustable parameter. The accuracy of the second-order HNC2 theory is also confirmed by the results presented in figure 1 (b). HNC2 theory well reproduces the location and the height of the first maximum of g(r) and the shoulder at distances r/σf < 1.4. However, the first minimum of g(r) at r/σf ≈ 3.32 is a bit better predicted by the RY approximation. Similar conclusions can be drawn from the inspection of the curves presented in figure 1 (c). Since the PY2 approximation actually fails to reproduce the simulation data with a reasonable accuracy, the PY and PY2 results are omitted in figure 1 (c) and in the following figures 2 and 3. The singlet PY approximation is quite accurate when applied to a hard-sphere system [33]. On the contrary, in the case of soft potentials, e.g., 13601-4 Core-softened model fluid 2 4 6 r/σ f 0 0.5 1 1.5 g( r) PY HNC RY PY2 HNC2 CGMC a 2 4 6 r/σ f 0 0.5 1 1.5 g( r) b 2 4 6 r/σ f 0 0.5 1 1.5 2 g( r) c Figure 1. (Color on-line) Radial distribution functions evaluated from the singlet (PY, HNC and RY), second-order (PY2, HNC2) theories and from GCMC simulation at ρ∗ = 0.1 and at T ∗ = 0.5 (part a), 0.3 (part b) and 0.2 (part c). The nomenclature of the lines is given in part a. Note that consecutive parts b and c contain the results of selected theories. Gaussian-like potential, the singlet HNC approximation works quite satisfactory [58]. It is thus not surprising that for core-softened potential, the hypernetted chain approximation (at singlet and second-order levels) performs better than the the Percus-Yevick closures. However, the observed big differences between the predictions of these approximations are surprising, because in the case of the systems studied so far [48–53] the differences between the PY2 and HNC2 approximations were rather small [59]. Further tests of the accuracy of the HNC2 approximation are shown in figures 2 and 3. Figure 2 shows the results at the density ρ∗ = 0.14, consecutive parts a – d are for the temperatures T ∗ = 0.5, 0.3, 0.2 and 0.15, respectively. Again, the HNC2 approximation is superior over the best singlet RY theory. In particular, it better describes the formation of the first peak of g(r) within the repulsive ramp of the potential (1) (the plot of the potential (1) gives figure 1 of [23]), though both RY and HNC2 theories underestimate the height of this maximum. By contrast, the singlet HNC approach describes much higher first peak of g(r) compared to GCMC data. Finally, we calculated the radial distribution functions at higher densities, ρ∗ = 0.16 and 0.20 (figure 3). For these densities, the second-order HNC2 equation performs well. In particular, its predictions of the height of consecutive maxima and minima of g(r) are more accurate than the singlet RY approximation. However, for some state points, the accuracy of the singlet HNC equation 13601-5 O. Pizio, Z. Sokołowska, S. Sokołowski 2 4 6 r/σ f 0 0.5 1 g( r) a 2 4 6 r/σ f 0 0.5 1 1.5 g( r) b 2 4 6 r/σ f 0 0.5 1 1.5 g( r) c 2 4 6 r/σ f 0 0.5 1 1.5 g( r) d Figure 2. (Color on-line) Radial distribution functions evaluated from the singlet (HNC and RY), second -order (HNC2) theories and from GCMC simulation at ρ∗ = 0.14 and at T ∗ = 0.5 (part a), 0.3 (part b), 0.2 (part c) and 0.15 (part d). The nomenclature of the lines is given in figure 1 (a). is quite good and even better than the HNC2 theory, cf. figures 3 (b) and 3 (c). Before discussing an anomalous behavior of the system in the pressure-temperature plane, that results from the second-order integral equations, we briefly recall previous findings. The line of temperatures of maximum density (TMD) was determined from NPT molecular dynamics [23] and singlet PY [23] and RY [23, 31] integral equations. It was also demonstrated [23, 31] that the singlet HNC approach fails to predict this anomaly. Molecular dynamic simulations revealed that for densities 0.12 < ρ∗ < 0.14, the pressure- temperature curves at constant density have a minimum (see also our comment in the introductory section) which implies density anomaly. The maximum temperature at which the density anomaly still exists is T ∗ max ≈ 0.25 [23]. The predictions of the singlet RY approximation depend on the method according to which its adjustable parameter was determined [31]. Barros de Oliveira et al. [23] obtained the adjustable parameter of the RY equation by checking the consistency between the compressibilities, calculated from the virial and compressibility equations of state. However, in our work [31] we imposed the so-called global coexistence criterion, according to which the pressures from two above mentioned thermodynamic routes were compared. Barros de Oliveira et al. [23] found that the density anomaly exists for the densities from the range of [0.12, 0.14] and that T ∗ max ≈ 0.23. Our calculations, however, led to a wider interval of the densities, [0.12, 0.19]. We should also note that the singlet PY approximation predicts [23] much wider region of densities, at 13601-6 Core-softened model fluid 2 4 r/σ f 0 0.5 1 g( r) a 2 4 r/σ f 0 0.5 1 g( r) b 2 4 r/σ f 0 0.5 1 g( r) c 2 4 r/σ f 0 0.5 1 g( r) d Figure 3. (Color on-line) Radial distribution functions evaluated from the singlet (HNC and RY), second -order (HNC2) theories and from GCMC simulation at ρ∗ = 0.16 and at T ∗ = 0.5 (part a), ρ∗ = 0.16 and at T ∗ = 0.3 (part b), ρ∗ = 0.2 and at T ∗ = 0.5 (part c), and ρ∗ = 0.2 and at T ∗ = 0.3 (part d). The nomenclature of the lines is given in figure 1 (a). which the density anomaly occurs, [0.13, 0.3]. Instantaneously, this region is shifted towards much higher temperatures, with T ∗ max = 0.86. Figure 4 shows the isochores evaluated from different theories (the reduced pressure is defined as P ∗ = Pσ3 f εf). At low densities (ρ∗ < 0.07; the relevant curves are not shown for the sake of brevity) none of the approximations predicts the existence of the density anomaly. For ρ∗ = 0.08 (part a) only the PY2 theory yields a weak minimum of the pressure, all the remaining approximations do not show any anomaly. However, it is interesting to note that the results of the singlet PY and HNC2 approximations almost coincide, except for low temperatures. For ρ∗ = 0.1 (part b) only the PY2 isochore exhibits anomaly, but now the minimum of the pressure is preceded by a maximum. This behavior is completely unexpected, but, it is rather an artifact of the theory, because no computer simulations [23] provide its confirmation. For ρ∗ = 0.14 (part c) all the approximations, but the HNC1 predict the density anomaly, whereas at higher bulk density, ρ∗ = 0.2, (part d) only PY1 and PY2 predict anomalies. The range of the densities for which the HNC2 leads to anomalous behavior is rather narrow, 0.13 6 ρ∗ 6 0.155. As we have shown in our previous work [31], the RY approximation with the global coexistence criterion predicts anomalous behavior up to ρ∗ ≈ 0.19. Surprisingly, the PY2 approximation yields two density anomaly regions: the first one is for the densities [0.07, 0.11] and the second one is for ρ∗ > 0.13. We must recall, however, that the structure predicted from PY2 13601-7 O. Pizio, Z. Sokołowska, S. Sokołowski 0 0.5 1 T * 0.9 1 P * 0 0.5 T * 0.2 0.3 P * HNC PY RY HNC2 PY2 0.5 0.4 0.5 P * 0 0.5 1 T * 1.8 1.9 2 2.1 P * ρ∗ =0.08a ρ∗ =0.14 ρ∗ =0.1 ρ∗ =0.2 b c d Figure 4. (Color on-line.) Isochores resulting from different approximate theories for different densities, given in consecutive parts. The nomenclature of the lines is displayed in part a. 0.2 0.4 0.6 0.8 T * 0.1 0.2 0.3 ρ∗ 0.15 0.3 T * 0.15 0.2 ρ∗ RY-gcc HNC2 RY MD PY Figure 5. (Color on-line) The TMD lines from different approximate theories and from MD simulations [23] (open circles). Abbreviations are as follows: HNC2 denotes the second-order hypernetted chain theory results; PY is the singlet Percus-Yevick theory [23], RY denote the Roger-Young results obtained by Barros de Oliveira et al. [23] and RY-gcc are the Roger-Young results calculated using the global coexistence criterion. The inset magnifies a part of the main figure. Dotted lines here indicate the thermodynamic conditions at which the bridge functions displayed in figure 6 were calculated. 13601-8 Core-softened model fluid approximation differs very much from the computer simulation results. Summary of our calculations is presented in figure 5. We show here TMD lines resulting from different approximate theories. Since the PY2 predictions quantitatively differ from the computer simulation results, the relevant curve has been omitted. The TMD line from the RY approach evaluated using global coexistence criterion differs from that obtained by Barros de Oliveira et al. [23] None of the approximate theories is capable of reproducing the simulation results at quantitative level. Finally, we show the results of the bridge function calculations. The bridge function for the model in question has never been investigated so far. The function γ(r), entering equation (8) was calculated from γ(r) = 1 2π2r ∫ ∞ 0 dkk sin(kr) ρh̃(k)2 1 + ρh̃(k) , (9) where h̃(k) is the Fourier transform of the function h(r). We should stress that the calculations of the Fourier transform h̃(k) and the inverse transform in equation (9) should be carried out with a special care, as described by Kolafa et al. [60]. Our calculations were performed at a constant density equal to ρ∗ = 0.152 for a set of temperatures across the TMD line, as well as at a constant temperature, T ∗ = 0.175 and for several densities, again across the TMD line (cf. figure 5). The functions B(r) significantly differ from the functions for hard-spheres [48, 60]. In particular, the decay of the bridge functions for the core-softened potential is much slower than for hard-spheres. Even for r/σf ≈ 10, the bridge functions exhibit well pronounced oscillations. These oscillations almost vanish at a distance as large as r/σf ≈ 16, whereas for hard-spheres at very high densities, the oscillations already vanish at r/σf ≈ 5, cf. figure 3 of [60]. Crossing the TMD line seems to have no effect on the shape of the bridge functions. The evolution of the functions B(r) along the temperature and the density branches (figure 6) is smooth and we do not observe any peculiarities connected with specific thermodynamic behavior of the system. However, the range of oscillations of a bridge function is worth stressing once again. 0 1 2 r/σ f -40 -30 -20 -10 0 B (r ) 2 4 6 8 10 12 r/σ f -0.2 -0.1 0 0.1 B (r ) 0.10 0.11 0.125 0.14 0.175 T * = a 0 1 2 r/σ f -30 -25 -20 -15 -10 -5 0 B (r ) 2 4 6 8 10 12 r/σ f -0.2 -0.1 0 0.1 B (r ) 0.13 0.14 0.16 ρ* = b Figure 6. (Color on-line) Changes of the bridge functions across the TMD line. Part a is at a fixed density, ρ∗ = 0.152, and for different temperatures given in the figure, whereas part b is at a fixed temperature, T ∗ = 0.175, and for three different densities given in the figure. The paths along which the calculations were performed are shown in the inset to figure 5 as dotted lines. Let us briefly summarize our findings. We have demonstrated that the radial distribution func- tions predicted by the second-order HNC2 approach are more accurate than the predictions of any singlet theory. In contrast to the singlet RY approach, the HNC2 approximation does not involve any adjustable parameter. However, the PY2 approximation is inaccurate and the distribution func- tions resulting from it significantly differ from those of all remaining theories. The discrepancies between PY2 and HNC2 approximations are puzzling and we cannot offer any physical interpre- tation of these trends. The second-order HNC2 approach yields the density anomaly in the system for densities from the interval [0.13, 0.15]. The maximum temperature at which this phenomenon is 13601-9 O. Pizio, Z. Sokołowska, S. Sokołowski observed is T ∗ max ≈ 0.2. Moreover, for the first time we have evaluated the bridge functions for the core softened fluid. Unlike the hard-sphere or Lennard-Jones fluids, the oscillations of the bridge functions extend over much larger interparticle separations. A successful parametrization of the bridge function similar to hard-sphere systems, see e.g. [60–62], would be helpful in developing new closures for singlet level theory for a class of core-softened potentials. References 1. Hemmer P.C., Stell G., Phys. Rev. Lett., 1970, 24, 1284; doi:10.1103/PhysRevLett.24.1284. 2. Stillinger F.H., Head-Gordon T., Phys. Rev. E, 1993, 47, 2484; doi:10.1103/PhysRevE.47.2484. 3. Cho C.H., Singh S., Robinson G.W., Phys. Rev. Lett., 1996, 76, 1651; doi:10.1103/PhysRevLett.76.1651. 4. Jagla E.A., Phys. Rev. E, 1998, 58, 1478; doi:10.1103/PhysRevE.58.1478. 5. Sadr-Lahijany M.R., Scala A., Buldyrev S.V., Stanley H.E., Phys. Rev. Lett., 1998, 81, 4895; doi:10.1103/PhysRevLett.81.4895. 6. Egorov S.A., J. Chem. Phys., 2008, 129, 024514; doi:10.1063/1.2953332. 7. Malescio G., Pellicane G., Phys. Rev. E, 2001, 63, 020501(R); doi:10.1103/PhysRevE.63.020501. 8. Kumar P. et al., Phys. Rev. E, 2005, 72, 021501; doi:10.1103/PhysRevE.72.021501. 9. Chakraborty S.N., Chakravarty C., Phys. Rev. E, 2007, 76, 011201; doi:10.1103/PhysRevE.76.011201. 10. Zhou S., Lajovic A., Jamnik A., J. Chem. Phys., 2008, 129, 124503; doi:10.1063/1.2982162. 11. Chaimovich A., Shell M.S., Phys. Chem. Chem. Phys., 2009, 11, 1901; doi:10.1039/b818512c. 12. Wilding N.B., Magee J.E., Phys. Rev. E, 2002, 66, 031509; doi:10.1103/PhysRevE.66.031509. 13. Gibson H.M., Wilding N.B., Phys. Rev. E, 2006, 73, 061507; doi:10.1103/PhysRevE.74.019903. 14. Yan Z.Y. et al., Phys. Rev. E, 2006, 73, 051204; doi:10.1103/PhysRevE.73.051204. 15. Krekelberg W.P., Mittal J., Ganesan V., Truskett T.M., Phys. Rev. E, 2008, 77 041201; doi:10.1103/PhysRevE.77.041201. 16. Barros de Oliveira A., Netz P.A., Barbosa M.C., Eur. Phys. J. B, 2008, 64, 481; doi:10.1140/epjb/e2008-00101-6. 17. Gribova N.V., Fomin Yu.D., Frenkel D., Ryzhov V.N., Phys. Rev. E, 2009, 79, 051202; doi:10.1103/PhysRevE.79.051202. 18. Lomba E., Almarza N.G., Mart́ın C., McBride C., J. Chem. Phys., 2007, 126, 244510; doi:10.1063/1.2748043. 19. Franzese G. et al., Phys. Rev. E, 2002, 66, 051206; doi:10.1103/PhysRevE.66.051206. 20. Malescio G. et al., Phys. Rev. E, 2005, 71, 061504; doi:10.1103/PhysRevE.71.061504. 21. Barraz N.M., Salcedo E., Barbosa M.C., J. Chem. Phys., 2009, 131, 094504; doi:10.1063/1.3213615. 22. Barros de Oliveira A., Netz P.A., Colla T., Barbosa M.C., J. Chem. Phys., 2006, 125, 124503; doi:10.1063/1.2357119. 23. Barros de Oliveira A., Netz P.A., Colla T., Barbosa M.C., J. Chem. Phys., 2006, 124, 084505; doi:10.1063/1.2168458. 24. Barros de Oliveira A., Barbosa M.C., Netz P.A., Physica A, 2007, 386, 744; doi:10.1016/j.physa.2007.07.015 . 25. Cho C.H., Singh S., Robinson G.W., J. Chem. Phys., 1997, 107, 7979; doi:10.1063/1.475060. 26. Netz P.A., Raymundi J.F., Camera A.S., Barbosa M.C., Physica A, 2004, 342, 48; doi:10.1016/j.physa.2004.04.058. 27. Jagla E.A., J. Chem. Phys., 1999, 110, 451; 1999, 111, 8980; doi:10.1063/1.478105. 28. Jagla E.A., Phys. Rev. E, 2001, 63, 061501; 2001, 63, 061509; doi:10.1103/PhysRevE.63.061501. 29. Dominguez H., Pizio O., Pusztai L., Soko lowski S., Adsorpt. Sci. Technol., 2007, 25, 479. 30. Pizio O., Dominguez H., Pusztai L., Soko lowski S., Physica A, 2009, 388, 2278; doi:10.1016/j.physa.2009.03.014. 31. Pizio O., Dominguez H., Duda Yu., Soko lowski S., J. Chem. Phys., 2009, 130 174504; doi:10.1063/1.3125930. 32. Barros de Oliveira A., Salcedo E., Chakravarty C., Barbosa M.C., J. Chem. Phys., 2010, 132, 234509; doi:10.1063/1.3429254. 33. Hansen J.P., McDonald I.R., Theory of Simple Liquids. Academic Press, 3rd edition, New York, 2006. 34. Henderson D., Soko lowski S., Borówko M., in: Computational Mathods in Surface and Colloid Science, ed. M. Borówko. Marcel Dekker, New York, 2000. 35. Soko lowski S., J. Chem. Phys., 1980, 73, 3507; doi:10.1063/1.440507. 13601-10 http://dx.doi.org/10.1103/PhysRevLett.24.1284 http://dx.doi.org/10.1103/PhysRevE.47.2484 http://dx.doi.org/10.1103/PhysRevLett.76.1651 http://dx.doi.org/10.1103/PhysRevE.58.1478 http://dx.doi.org/10.1103/PhysRevLett.81.4895 http://dx.doi.org/10.1063/1.2953332 http://dx.doi.org/10.1103/PhysRevE.63.020501 http://dx.doi.org/10.1103/PhysRevE.72.021501 http://dx.doi.org/10.1103/PhysRevE.76.011201 http://dx.doi.org/10.1063/1.2982162 http://dx.doi.org/10.1039/b818512c http://dx.doi.org/10.1103/PhysRevE.66.031509 http://dx.doi.org/10.1103/PhysRevE.74.019903 http://dx.doi.org/10.1103/PhysRevE.73.051204 http://dx.doi.org/10.1103/PhysRevE.77.041201 http://dx.doi.org/10.1140/epjb/e2008-00101-6 http://dx.doi.org/10.1103/PhysRevE.79.051202 http://dx.doi.org/10.1063/1.2748043 http://dx.doi.org/10.1103/PhysRevE.66.051206 http://dx.doi.org/10.1103/PhysRevE.71.061504 http://dx.doi.org/10.1063/1.3213615 http://dx.doi.org/10.1063/1.2357119 http://dx.doi.org/10.1063/1.2168458 http://dx.doi.org/10.1016/j.physa.2007.07.015 http://dx.doi.org/10.1063/1.475060 http://dx.doi.org/10.1016/j.physa.2004.04.058 http://dx.doi.org/10.1063/1.478105 http://dx.doi.org/10.1103/PhysRevE.63.061501 http://dx.doi.org/10.1016/j.physa.2009.03.014 http://dx.doi.org/10.1063/1.3125930 http://dx.doi.org/10.1063/1.3429254 http://dx.doi.org/10.1063/1.440507 Core-softened model fluid 36. Soko lowski S., Mol. Phys, 1983, 49, 1481; doi:10.1080/00268978300102091. 37. Plischke M., Henderson D., J. Phys. Chem., 1984, 88, 6544; doi:10.1021/j150670a015. 38. Henderson D., Soko lowski S., Trokhymchuk A., J. Chem. Phys., 1995, 103, 4693; doi:10.1063/1.470656. 39. Henderson D., Pizio O., Trokhymchuk A., Soko lowski S., Physica A, 1997, 244, 147; doi:10.1016/S0378-4371(97)00237-9. 40. Omelyan I., Hirata F., Kovalenko A., Chem. Phys. Phys. Chem., 2005, 7, 4132; doi:10.1039/b507761c. 41. Percus J.K., Phys. Rev. Lett., 1962, 8, 462; doi:10.1103/PhysRevLett.8.462. 42. Percus J.K. – In: The Equlibrium Theory of Classical Fluids, Eds. H.L. Frisch and J.L. Lebowitz. Benjamin, New York, 1964, 33. 43. Attard P., J. Chem. Phys., 1989, 91, 3072; doi:10.1063/1.456930. 44. Attard P., J. Chem. Phys., 1989, 91, 3083; doi:10.1063/1.456931. 45. Henderson D., Condens. Matter Phys., 1993, 2, 31. 46. Holovko M.F., Pizio O. Preprint ICMP–91–2U, 1991 (in Ukrainian). 47. Holovko M.F., Pizio O., Haymet A.D.J., Henderson D. Preprint ITP–91–55E, 1991. 48. Henderson D., Soko lowski S., J. Chem. Phys., 1995, 103, 7541; doi:10.1063/1.470322. 49. Henderson D., Soko lowski S., J. Chem. Phys., 1996, 104, 2971; doi:10.1063/1.471118. 50. Henderson D., Soko lowski S., Condens. Matter Phys., 1996, 7, 51. 51. Duh D.M., Henderson D., Mier y Teran L., Soko lowski S., Mol. Phys., 1997, 90, 563; doi:10.1080/00268979709482638. 52. Henderson D., Soko lowski S., Mol. Phys., 1997, 90, 85; doi:10.1080/002689797172895. 53. Henderson D., Soko lowski S., Wasan D., Mol. Phys., 1998, 93, 294; doi:10.1080/002689798169302. 54. Bryk P., Pizio O., Soko lowski S., J. Chem. Soc. Faraday Trans., 1997, 93, 2367; doi:10.1039/a608565b. 55. Lovett R., Mou C.Y., Buff F.P., J. Chem. Phys., 1976, 65, 570; doi:10.1063/1.433110. 56. Wertheim M.S., J. Chem. Phys., 1976, 65, 2377; doi:10.1063/1.433352. 57. Yukhnovski I.R., Holovko M.F., Statistical Theory of Classical Equilibrium Systems. Naukova Dumka, Kiev, 1980. 58. Louis A.A., Bolhuis P.G., Hansen J.P., Phys. Rev. E, 2000, 62, 7961; doi:10.1103/PhysRevE.62.7961. 59. Henderson D. – In: Fundamentals of Inhomogeneous Fluids, ed. D. Henderson. Marcel Dekker, New York, 1992, Chapter 4. 60. Kolafa J., Lab́ık S., Malijevský A., Mol. Phys., 2002, 100, 2629; doi:10.1080/00268970210136357. 61. Grundke E.W., Henderson D., Mol. Phys., 1972, 24, 269; doi:10.1080/00268977200101431; Henderson D., Grundke E.W., J. Chem. Phys., 1975, 63, 601; doi:10.1063/1.431378. 62. Lab́ık S., Malijevský A., Mol. Phys., 1984, 53, 381; doi:10.1080/00268978400102381; Malijevský A., Lab́ık S., Mol. Phys., 1987, 60, 663; doi:10.1080/00268978700100441. 13601-11 http://dx.doi.org/10.1080/00268978300102091 http://dx.doi.org/10.1021/j150670a015 http://dx.doi.org/10.1063/1.470656 http://dx.doi.org/10.1016/S0378-4371(97)00237-9 http://dx.doi.org/10.1039/b507761c http://dx.doi.org/10.1103/PhysRevLett.8.462 http://dx.doi.org/10.1063/1.456930 http://dx.doi.org/10.1063/1.456931 http://dx.doi.org/10.1063/1.470322 http://dx.doi.org/10.1063/1.471118 http://dx.doi.org/10.1080/00268979709482638 http://dx.doi.org/10.1080/002689797172895 http://dx.doi.org/10.1080/002689798169302 http://dx.doi.org/10.1039/a608565b http://dx.doi.org/10.1063/1.433110 http://dx.doi.org/10.1063/1.433352 http://dx.doi.org/10.1103/PhysRevE.62.7961 http://dx.doi.org/10.1080/00268970210136357 http://dx.doi.org/10.1080/00268977200101431 http://dx.doi.org/10.1063/1.431378 http://dx.doi.org/10.1080/00268978400102381 http://dx.doi.org/10.1080/00268978700100441 O. Pizio, Z. Sokołowska, S. Sokołowski Мiкроскопiчна структура i термодинамiка моделi плину з пом’якшеним кором з теорiї iнтегральних рiвнянь другого порядку O. Пiзiо1, З. Соколовська2, С. Соколовський3 1 Iнститут хiмiї УГАМ, Койокан, Мексика 2 Iнститут агрофiзики Польської академiї наук, Люблiн, Республiка Польща 3 Унiверситет iм. Марiї Складовської-Кюрi,Люблiн, Республiка Польща Ми дослiдили структурнi i термодинамiчнi властивостi однорiдного тривимiрного плину з пом’як- шеним кором, використовуючи iнтегральнi рiвняння Орнштейна-Цернiке другого порядку iз гiперланцюговим замиканням та замиканням Перкуса-Євiка. Зроблено порiвняння радiальних функцiй розподiлу з вiдповiдними функцiями, отриманими iз синглетних iнтегральних рiвнянь, а також з даними комп’ютерного моделювання. З рiзних наближених теорiй встановлено границi областi аномалiї густини. Отриманi результати показують, що гiперланцюгове наближення iнтегральних рiвнянь другого порядку може бути використане для опису як структури, так i аномалiї густини цього модельного плину. Крiм того, ми представляємо результати обчислень мiсткових функцiй. Ключовi слова: теорiя рiдин, iнтегральнi рiвняння другого порядку, аномалiя густини, мiстковi функцiї 13601-12 Introduction Theory Results and discussion

Microscopic structure and thermodynamics of a core-softened model fluid from the second-order integral equations theory

Репозитарії

Схожі ресурси