Density anomaly of charged hard spheres of different diameters in a mixture with core-softened model solvent. Monte Carlo simulation results Author(s):

Very recently the effect of equisized charged hard sphere solutes in a mixture with core-softened fluid model on the structural and thermodynamic anomalies of the system has been explored in detail by using Monte Carlo simulations and integral equations theory (J. Chem. Phys., 137, 244502 (2012)). O...

Ausführliche Beschreibung

Gespeichert in:

Bibliographische Detailangaben
Datum:	2013
Hauptverfasser:	Hribar-Lee, B., Pizio, O.
Format:	Artikel
Sprache:	English
Veröffentlicht:	Інститут фізики конденсованих систем НАН України 2013
Schriftenreihe:	Condensed Matter Physics
Online Zugang:	http://dspace.nbuv.gov.ua/handle/123456789/120858
Tags:	Tag hinzufügen Keine Tags, Fügen Sie den ersten Tag hinzu!
Назва журналу:	Digital Library of Periodicals of National Academy of Sciences of Ukraine
Zitieren:	Density anomaly of charged hard spheres of different diameters in a mixture with core-softened model solvent. Monte Carlo simulation results / B. Hribar-Lee, O. Pizio// Condensed Matter Physics. — 2013. — Т. 16, № 4. — С. 43607:1-10. — Бібліогр.: 42 назв. — англ.

Institution

Digital Library of Periodicals of National Academy of Sciences of Ukraine

id	irk-123456789-120858
record_format	dspace
spelling	irk-123456789-1208582017-06-14T03:04:31Z Density anomaly of charged hard spheres of different diameters in a mixture with core-softened model solvent. Monte Carlo simulation results Author(s): Hribar-Lee, B. Pizio, O. Very recently the effect of equisized charged hard sphere solutes in a mixture with core-softened fluid model on the structural and thermodynamic anomalies of the system has been explored in detail by using Monte Carlo simulations and integral equations theory (J. Chem. Phys., 137, 244502 (2012)). Our objective of the present short work is to complement this study by considering univalent ions of unequal diameters in a mixture with the same soft-core fluid model. Specifically, we are interested in the analysis of changes of the temperature of maximum density (TMD) lines with ion concentration for three model salt solutes, namely sodium chloride, potassium chloride and rubidium chloride models. We resort to Monte Carlo simulations for this purpose. Our discussion also involves the dependences of the pair contribution to excess entropy and of constant volume heat capacity on the temperature of maximum density line. Some examples of the microscopic structure of mixtures in question in terms of pair distributions functions are given in addition. Недавно було дослiджено вплив заряджених твердих сфер однакового розмiру як розчиненої речовини на структурнi i термодинамiчнi аномалiї системи, в якiй застосовано модель плину з потенцiалом з м’яким кором як розчинник. Застосовувався метод моделювання Монте Карло та теорiя iнтегральних рiвнянь [J. Chem. Phys., 2012, 137, 244502]. Метою даної роботи є доповнити попереднє досл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лу приведено на додаток. 2013 Article Density anomaly of charged hard spheres of different diameters in a mixture with core-softened model solvent. Monte Carlo simulation results / B. Hribar-Lee, O. Pizio// Condensed Matter Physics. — 2013. — Т. 16, № 4. — С. 43607:1-10. — Бібліогр.: 42 назв. — англ. 1607-324X PACS: 61.20.-p (61.20.Gy, 6120.Ja, 61.20.Ne, 61.20.Qg), 65.20.-w DOI:10.5488/CMP.16.43607 arXiv:1312.4566 http://dspace.nbuv.gov.ua/handle/123456789/120858 en Condensed Matter Physics Інститут фізики конденсованих систем НАН України
institution	Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection	DSpace DC
language	English
description	Very recently the effect of equisized charged hard sphere solutes in a mixture with core-softened fluid model on the structural and thermodynamic anomalies of the system has been explored in detail by using Monte Carlo simulations and integral equations theory (J. Chem. Phys., 137, 244502 (2012)). Our objective of the present short work is to complement this study by considering univalent ions of unequal diameters in a mixture with the same soft-core fluid model. Specifically, we are interested in the analysis of changes of the temperature of maximum density (TMD) lines with ion concentration for three model salt solutes, namely sodium chloride, potassium chloride and rubidium chloride models. We resort to Monte Carlo simulations for this purpose. Our discussion also involves the dependences of the pair contribution to excess entropy and of constant volume heat capacity on the temperature of maximum density line. Some examples of the microscopic structure of mixtures in question in terms of pair distributions functions are given in addition.
format	Article
author	Hribar-Lee, B. Pizio, O.
spellingShingle	Hribar-Lee, B. Pizio, O. Density anomaly of charged hard spheres of different diameters in a mixture with core-softened model solvent. Monte Carlo simulation results Author(s): Condensed Matter Physics
author_facet	Hribar-Lee, B. Pizio, O.
author_sort	Hribar-Lee, B.
title	Density anomaly of charged hard spheres of different diameters in a mixture with core-softened model solvent. Monte Carlo simulation results Author(s):
title_short	Density anomaly of charged hard spheres of different diameters in a mixture with core-softened model solvent. Monte Carlo simulation results Author(s):
title_full	Density anomaly of charged hard spheres of different diameters in a mixture with core-softened model solvent. Monte Carlo simulation results Author(s):
title_fullStr	Density anomaly of charged hard spheres of different diameters in a mixture with core-softened model solvent. Monte Carlo simulation results Author(s):
title_full_unstemmed	Density anomaly of charged hard spheres of different diameters in a mixture with core-softened model solvent. Monte Carlo simulation results Author(s):
title_sort	density anomaly of charged hard spheres of different diameters in a mixture with core-softened model solvent. monte carlo simulation results author(s):
publisher	Інститут фізики конденсованих систем НАН України
publishDate	2013
url	http://dspace.nbuv.gov.ua/handle/123456789/120858
citation_txt	Density anomaly of charged hard spheres of different diameters in a mixture with core-softened model solvent. Monte Carlo simulation results / B. Hribar-Lee, O. Pizio// Condensed Matter Physics. — 2013. — Т. 16, № 4. — С. 43607:1-10. — Бібліогр.: 42 назв. — англ.
series	Condensed Matter Physics
work_keys_str_mv	AT hribarleeb densityanomalyofchargedhardspheresofdifferentdiametersinamixturewithcoresoftenedmodelsolventmontecarlosimulationresultsauthors AT pizioo densityanomalyofchargedhardspheresofdifferentdiametersinamixturewithcoresoftenedmodelsolventmontecarlosimulationresultsauthors
first_indexed	2025-07-08T18:44:39Z
last_indexed	2025-07-08T18:44:39Z
_version_	1837105449339453440
fulltext	Condensed Matter Physics, 2013, Vol. 16, No 4, 43607: 1–10 DOI: 10.5488/CMP.16.43607 http://www.icmp.lviv.ua/journal Density anomaly of charged hard spheres of different diameters in a mixture with core-softened model solvent. Monte Carlo simulation results B. Hribar-Lee1, O. Pizio2¤ 1 Faculty of Chemistry and Chemical Technology, University of Ljubljana, Aškerčeva c. 5, SI-1000 Ljubljana, Slovenia, 2 Instituto de Química, Universidad Nacional Autónoma de México, Circuito Exterior, 04510, México, D.F., México Received August 7, 2013, in final form September 12, 2013 Very recently the effect of equisized charged hard sphere solutes in a mixture with core-softened fluid model on the structural and thermodynamic anomalies of the system has been explored in detail by using Monte Carlo simulations and integral equations theory [J. Chem. Phys., 2012, 137, 244502]. Our objective of the present short work is to complement this study by considering univalent ions of unequal diameters in a mixture with the same soft-core fluid model. Specifically, we are interested in the analysis of changes of the temperature of maximum density (TMD) lines with ion concentration for three model salt solutes, namely sodium chloride, potassium chloride and rubidium chloride models. We resort to Monte Carlo simulations for this purpose. Our discussion also involves the dependences of the pair contribution to excess entropy and of constant volume heat capacity on the temperature of maximum density line. Some examples of the microscopic structure of mixtures in question in terms of pair distribution functions are given in addition. Key words: liquid theory, soft-core model fluid, electrolyte solution, canonical Monte Carlo simulations, density anomaly PACS: 61.20.-p, 61.20.Gy, 6120.Ja, 61.20.Ne, 61.20.Qg, 65.20.-w 1. Introduction It is our honor and pleasure to make contribution to this special issue dedicated to brilliant scien- tist Prof. Myroslav Holovko on the occasion of his 70th birthday. During last years we were influenced by his developments and profound results along different lines of research in the theory of liquids and mixtures, interfacial phenomena and phase transitions. We appreciate his mentorship in some projects, many useful pieces of advice and friendly discussions concerning topics of our common interest. During many years Myroslav Holovko explored electrolyte solutions in very detail, see e.g., the classical mono- graph [1] written by him together with his teacher Prof. I. Yukhnovskyi that summarizes their effort prior to 1980. Later results can be found in many articles and reviews. The principal objective of the present work is to complete one extension of our very recent study of solutions comprising solute ions and a soft-core model fluid as a solvent [2]. During last decade many studies in the theory of fluids have been focused on single-component, isotropic core-softened models in which the repulsive part of the interparticle interaction potential is peculiar (“softer”) at “intermediate” interparticle separations while a hard-core repulsion acts at short distances. This region of peculiar de- pendence of the pair interparticle interaction is often described by linear or nonlinear ramp [3–8], by a shoulder [9, 10] or Gaussian term [11–13], or another more sophisticated functional dependence, see e.g., recent review [14] and references therein. A more complete list of references has been given in [2] as well. ¤E-mail: pizio@unam.mx © B. Hribar-Lee, O. Pizio, 2013 43607-1 http://dx.doi.org/10.5488/CMP.16.43607 http://www.icmp.lviv.ua/journal B. Hribar-Lee, O. Pizio Unfortunately, the complex shape of the interaction in this type of models requires application of solely numerical procedures for their solution. However, peculiar structural, thermodynamic and dy- namic properties relevant to certain real liquids can be reached. Models of single-component fluids with isotropic core-softened potential have been studied by computer simulations in numerous works, see the list of references in [2]. For purposes of the present work it is necessary to mention that the fluids with spherically symmetric core-softened potentials exhibit anomalies in thermodynamic and dynamic prop- erties in close similarity to those observed in liquids with directional inter-particle interactions, such as water, silica, phosphorous, and some liquid metals. One of the examples of such anomalous behavior is the presence of the fluid density maximum as a function of temperature. Moreover, some of the soft- core fluid models exhibit puzzling liquid-liquid coexistence envelope terminating at the critical point, in addition to common liquid-vapor phase transition. Thermodynamic and other anomalies can be related to this additional phase transition, see e.g., references [15–17]. The soft-core models by design miss the orientational interparticle correlations. Thus, their correspondence to, for example, water or aqueous solutions is restricted to thermodynamic properties. Nevertheless, the soft-core models are simple and permit to capture the peculiarities of thermodynamic properties quite adequately [14]. Theoretical studies of solutions in which a model solvent exhibits thermodynamic properties with anomalies have been attempted recently by different authors, see e.g., [18–24]. Several important par- ticular issues were explored, namely the solubility of apolar solutes in such solvents, the effect of solute species on the temperature of maximum density line, on heat capacity extrema and others. On the other hand, following our previous theoretical and simulation developments [25, 26], we applied similar mod- elling, but involved ionic solutes rather than nonpolar species [2]. We examined the effect of ionic solutes on the microscopic structure and thermodynamics of solutions with a soft-core solvent that possesses anomalies in one-phase region. Experiments [27–30] and computer simulations [31–33] have shown that certain anomalies of thermodynamic properties are preserved in ionic aqueous solutions for low and up to moderate ion concentrations. They become weaker and finally disappear, if the concentration of ions increases but remains below the solubility limit. In contrast to the initial stage of this project [2], we consider a mixture of singly charged cations and anions modelled as charged hard spheres and solvent particles described by a soft-core model, all of them having nonequal diameters. The system is studied by the canonical Monte Carlo simulations. The following discussion of the model and methods is brief in order to avoid unnecessary repetitions of [2], only the essential elements are given to benefit the reader. 2. Model and method The system in question consists of ionic solute (salt) and solvent. The ionic subsystem is a mixture of positively and negatively charged hard spheres with arbitrary diameters,¾ Å and¾ ¡ , the charge is located at the particle center. Ions and solvent particles are immersed in a dielectric background characterized by dielectric constant ". The interaction potential between ions i and j is, U i j (r )Æ 8 < : 1 , r Ç (¾ i Å¾ j )/2, e 2 z i z j 4¼" 0 "r , r Ê (¾ i Å¾ j )/2, (1) where e is the elementary charge, z i Æ jz ¡ j Æ 1 are the valencies of ions, " 0 is the permittivity of free space, and r is the inter-particle distance. As in [2], we use the core-softened model by Barros de Oliveira et al. [11–13] that possesses thermodynamic, structural, and dynamic anomalies. The potential is a sum of the Lennard-Jones and a Gaussian terms, U ss (r )Æ 4" ss · ³ ¾ s r ´ 12 ¡ ³ ¾ s r ´ 6 ¸ Åa" s exp · ¡ 1 2 µ r ¡ r 0 ¾ s ¶ 2 ¸ , (2) with " ss — the solvent-solvent interaction energy and ¾ s — the diameter of a solvent particle. Parameters a, , and r 0 are taken as in the original work [11] (a Æ 5, Æ 1, and r 0 /¾ s Æ 0.7). With this choice of parameters, the inter-particle attraction is strongly depressed and thus there is no liquid-gas transition in the solvent subsystem [34]. 43607-2 Density anomaly. Monte Carlo simulations The pressure-temperature curves calculated for the pure solvent at certain densities have aminimum (i.e., (�P/�T ) ½ Æ 0) and consequently (�½/�T ) P vanishes for some thermodynamic states. The region of density anomaly corresponds to the state points for which (�½/�T ) P È 0 and is bounded by the locus of points for which the thermal expansion coefficient equals zero. There is some room for the choice of ion-solvent interaction. However, to be consistent with our previous work [2], the interaction between an ion (i ), and a solvent (s) particle is taken simply in the form of Lennard-Jones (12-6) potential, U is (r )Æ 4" is · ³ ¾ is r ´ 12 ¡ ³ ¾ is r ´ 6 ¸ , (3) where " is is the ion-solvent interaction energy and ¾ is Æ (¾ i Å¾ s )/2. The reduced temperature for the solvent subsystem is defined as T ¤ Æ k B T /" s , where k B is the Boltzmann constant and T is the absolute temperature. In order to describe the energetic aspects of interactions between all species, it is necessary to define the parameter ®Æ T ¤ /T ¤ el as the ratio between the reduced temperature for the solvent subsystem (T ¤), and the reduced temperature of the ionic subsystem, T ¤ el Æ ¾ ¡ /¸ B , where ¸ B Æ e 2 /(4¼" 0 "k B T ) is the Bjerrum length, ¸ B Æ 7.14 Å. Similar to the previous study [2], we choose ® Æ 0.25. Thus, if the solvent is even at a rather low temperature, the ionic system remains far from the possible vapor – liquid type criticality. The mixing rule for ion-solvent potential is applied such that: " is Æ (" i i " j j ) 1/2, where " ss Æ 1/T ¤, " ¡¡ Æ 1/T ¤ el and " ÅÅ Æ (¾ ¡ /¾ Å )/T ¤ el . We performed Monte Carlo computer simulations in the canonical ensemble and in the cubic sim- ulation box according to common algorithms [35, 36]. The simulation was started from a configura- tion of randomly distributed particles. The number of ions in the box, N ion (N ion Æ N Å ÅN ¡ ) was cho- sen to provide their reasonable statistics and along with the desired ion concentration, x ion [x ion Æ N ion /(N ion ÅN s )], where N s is the number of solvent particles resulting from N ion and x ion . In addition, we picked up the desired total dimensionless density of the system, ½¤, ½¤ Æ (N Å ¾ 3 Å ÅN ¡ ¾ 3 ¡ ÅN s ¾ 3 s )/L 3, where L is the length of the box edge resulting from previously given parameters. Coulomb interactions were taken into account by using the Ewald summation [35, 36] with 337 wave vectors and the inverse length equal to 5/L). In this work, we explore three model salts in the soft-core solvent, namely the NaCl (¾ Å Æ 2.21 Å, ¾ ¡ Æ 4.42 Å), KCl (¾ Å Æ 4.74 Å, ¾ ¡ Æ 4.42 Å) and RbCl (¾ Å Æ 5.27 Å, ¾ ¡ Æ 4.42 Å) see e.g., [37, 38]. In all cases, we assume ¾ s Æ 3.166 Å in close similarity to the SPCE model of water. Each of the systems was first equilibrated over 10 8 attempted configurations followed by a series of production runs of at least 5£10 7 attempted particle displacements to obtain the relevant average properties of the system under investigation (excess internal energy, compressibility factor, heat capacity at constant volume, pair excess entropy, and radial distribution functions). The number of ions present in the simulation box was in the range between 150 and 400 (dependent on the desired density). The excess internal energy of the system was obtained directly as an ensemble average. The equation of state was calculated by [39] PV Nk B T Æ 1¡ 1 3N ¿ r �¯U tot �r À Å 2¼ 3 X i X j x i N j V ¾ 3 i j g i j (¾ i j ) , (i , j ÆÅ,¡), (4) whereU tot is the total excess internal energy of the system and g i j (¾ i j ) is the contact value of the appro- priate ion-ion radial distribution function. Pair excess entropy was determined from the expression [40] S 2,ex k B Æ¡2¼½ X i ,i x i x j Z £ g i j (r ) lng i j (r )¡ g i j (r )Å1 ¤ r 2 dr , (5) where x i is the number fraction of species i , and g i j (r ) are the pair distribution functions of species i and j . The heat capacity at constant volume was determined from the energy fluctuations as C V Æ ¡ hU 2 tot i¡hU tot i 2 ¢ /k B T 2 . (6) 43607-3 B. Hribar-Lee, O. Pizio 3. Results and discussion First series of our calculations concern the NaCl salt model (¾ Å Æ 2.21 Å, ¾ ¡ Æ 4.42 Å) in the soft- core solvent. The cation is rather small and its diameter is smaller compared to the anion and to the diameter of solvent particles. The results for P¤ (P¤ Æ P¾ 3 s /" s ) as a function of T ¤ for a set of densities of the mixture at fixed ion concentration x ion Æ 0.1 are given by short-dashed lines in figure 1 (a). Each line has a minimum (it is difficult to see it on the scale of such combined figure) marked by a circle. By joining the circles, one can construct the temperature of maximum density line at this particular ion concentration, x ion Æ 0.1. The P ¤ (T ¤ ) lines at a given ½ ¤ are built by joining points obtained as an average from several runs. Actually, it is quite a tedious procedure. We perform a series of calculations to construct the TMD lines for different x ion , up to x ion Æ 0.3. Each of the lines has an extremal point that shifts to a higher pressure and lower temperature with an increasing x ion . Moreover, the TMD lines shrink with an increasing ion concentration. These trends were already discussed by us for the model in which all the species involved are of equal diameter [2]. Slightly different view on the changes of the TMD lines is given in figure 1 (b). The shift of the TMD lines to higher ½¤ with a growing x ion can be explained by augmenting the excluded volume effects in the system, if ions are added to the solvent. On the other hand, the extremum point of the TMD along the temperature axis goes down due to the overall growth of attraction, principally between oppositely charged ions. Previously, we discussed the evolution of the microscopic structure of the model characterized by equal diameters of all species on the density at a fixed T ¤ and x ion , see figures 1 and 2 of [2]. Trends of behavior of the pair distribution functions are similar to the model with non-equal diameters. Two examples of the distribution of particles at thermodynamic states on the relevant TMD line are given here in figures 2 and 3. In the case of a lower density and lower ion concentration (figure 2), we can attribute the shape of the ion-ion distribution functions to the formation of Na–S–Cl, Cl–S–Cl and Na– S–Na structures or complexes, i.e., of ions separated by a single solvent molecule, corresponding to the local maxima at (¾ Na Å¾ Cl )/2Åw 1 , ¾ Cl Åw 1 and ¾ Na /2Åw 1 distances, respectively. Here, w 1 is the distance corresponding to the first maximum of the distribution of solvent species, g ww (r ) at 3.63 Å, approximately. Similarly, one can observe trends of formation of clusters at a distance (¾ Na Å¾ Cl )/2Åw 2 , where w 2 is the distance corresponding to the position of the principal maximum of the distribution of solvent species at 7 Å, approximately. Due to a stronger repulsion between small cations, the formation of Na–S–Na structure is less favorable compared to Cl–S–Cl structure at these conditions. On the other hand, if one considers the model at x ion Æ 0.3 at a point of the TMD line with slightly higher total density but substantially lower temperature (figure 3) compared to the case described in figure 2, the inter-ion electrostatic correlations are expected to be stronger. Nevertheless, screening of 0.1 0.15 0.2 T * 0.5 0.6 0.7 0.8 0.9 1 P * x ion = 0.1 0.2 0.25 x ion = 0.3; 0.29; 0.27 a 0.125 0.15 0.175 0.2 0.225 ρ* 0.1 0.125 0.15 0.175 0.2 T * TMD lines NaCl x ion = 0.1 x ion = 0.2 x ion = 0.25 x ion = 0.27 x ion = 0.28 x ion = 0.29 x ion = 0.3 b Figure 1. (Color online) Temperature of maximum density lines of the NaCl model in the P¤ ¡T¤ plane (panel a) T¤¡½¤ plane (panel b). 43607-4 Density anomaly. Monte Carlo simulations 5 10 15 20 25 r(Å) 0 1 2 3 g ij (r) NaCl a T * = 0.175 ρ* = 0.16 x ion = 0.1 Na - Cl Cl - Cl Na - Na 5 10 15 20 25 r(Å) 0 1 2 3 4 5 g ij (r) b Na - w Cl - w w - w Figure 2. (Color online) Pair distribution functions of the NaCl model. Parameters of the system are given in the figure. interactions between ions is expected to be more pronounced. Under the conditions given in figure 3, one can see that the contact value of g NaCl (r ) is practically of the same order as in figure 2. However, the height of the shoulder of this function at r ¼ 6.9 Å describing the correlation between oppositely charged ions separated by a single solvent particle is lower than in figure 2. On the other hand, the height of the first maximum of g ClCl is higher as an evidence of weaker repulsion between equally charged anions. Possible formation of Na–w 1 –Na is also more pronounced, cf. figure 2. The cation-solvent correlations are enhanced in the present case (at this low temperature) compared to anion-solvent correlation and compared to the case discussed in figure 2. We now again return to thermodynamic manifestations of the existence of the TMD lines. A quan- titative cumulative measure of the distribution of particles of different species is commonly given in terms of the pair contribution to the excess internal energy per particle as a function of the total density for different compositions of the mixture and at different reduced temperatures. In previous studies of the single-component soft-core model fluids it was suggested that the relation between the density and diffusion anomalies and the microscopic structure of the system can be explored in terms of the excess- entropy-based framework, see e.g., [41]. It has been shown that the presence of the density and the dif- 5 10 15 20 25 r(Å) 0 1 2 3 g ij (r) a Na - Cl Cl - Cl Na - Na NaCl x ion = 0.3 T * = 0.1325 ρ* = 0.21 5 10 15 20 25 r(Å) 0 2 4 6 g ij (r) b Na - w Cl - w w - w Figure 3. (Color online) Pair distribution functions of the NaCl model. Parameters of the system are given in the figure. 43607-5 B. Hribar-Lee, O. Pizio 0.12 0.15 0.18 0.21 T * tmd -2.5 -2 -1.5 -1 S 2, ex * a x ion = 0.1 0.2 0.25 0.3 0.12 0.15 0.18 0.21 T * tmd 1.5 1.8 2.1 2.4 C V * b x ion = 0.1 0.2 0.25 0.3 Figure 4. (Color online) The behavior of the pair contribution to excess entropy (panel a) and of reduced heat capacity (panel b) on the temperature of maximum density of the NaCl model at different concen- tration of ions. fusion anomalies is related to the effect of excess entropy on the fluid density. Our discussion concerns S 2,ex defined for a mixture by equation (5). However, in contrast to previous works, we have calculated the dependence of S 2,ex as a function of the temperature of maximum density rather than the total density of the model. In this respect, our results given in figure 4 (a) correspond to the TMD lines presented in figure 1. For example, the left hand part of the TMD line for x ion Æ 0.1 in figure 1 (b) that is characterized by dT ¤ /d½ ¤ È 0 corresponds to the upper part (lower absolute values) of the corresponding line for S 2,ex on T ¤ tmd . What can we learn from these functions? First observation is that the S 2,ex (T ¤ tmd ) curves shift to lower values of S 2,ex and lower values for T ¤ tmd with an increasing concentration of ionic solutes. Moreover, these curves shrink with increasing ion concentration. At low values of x ion (and at vanishing x ion as well) the curve S 2,ex (T ¤ tmd ) has a maximum such that dS 2,ex /d(T ¤ tmd )Ç 0 in the small region of T ¤ tmd after the point dS 2,ex /d(T ¤ tmd )Æ 0. This point does not coincide with the maximum temperature of the TMD line. Discontinuity of this derivative is expected, however, around the point at which T ¤ tmd reaches its maximum. This behavior is much less pronounced at high values for ion concentration in a mixture. The existence of such type of peculiar behavior seems to be worth to investigate in more detail for other soft-core solvent models and for say nonpolar solutes in various solvents for the sake of comparison with the case of ionic solutes. A rather similar discussion can be presented while considering the dependence of the reduced con- stant volume heat capacity, C¤ V on T ¤ tmd , figure 4 (b) [C¤ V Æ C V /(N ion ÅN s )k B ]. The heat capacity was “measured” in simulations through the fluctuations of internal energy according to equation (6). It is quite difficult to evaluate the heat capacity precisely due to fluctuations, especially at low temperatures. One way to suppress the fluctuations of the C¤ V values is to consider bigger boxes with a bigger number of particles, but then the simulations become quite lengthy in time. Shift and shrinking of the C¤ V (T ¤ tmd ) curves with an increasing ion concentration are well pronounced. However, the behavior of C¤ V (T ¤ tmd ) around the extremum point of the TMD line must be studied in more detail. After rather comprehensively examining the properties of the TMD lines of the model NaCl salt in the soft-core solvent in question at different ion concentrations, we would like to proceed to the results obtained for other two salts with bigger cations, namely the KCl and RbCl. The assumed diameters of cations were mentioned in the previous section, but we repeat the numbers for the sake of convenience of the reader KCl (¾ Å Æ 4.74 Å) and RbCl (¾ Å Æ 5.27 Å). In the former case, the diameter of cation does not differ much from the diameter of the anion ¾ Cl Æ 4.42 Å. From the analysis of the T ¤¡½¤ projections of the TMD lines given in figure 5 (a), we would like to conclude that in the present formulation of the interactions in the model, the principal factor determining the location of the TMD lines seems to be the excluded volume effects. Namely, at x ion Æ 0.1, the difference between TMD lines is reasonably big, if one 43607-6 Density anomaly. Monte Carlo simulations 0.12 0.15 0.18 0.21 0.24 0.27 ρ* 0.1 0.15 0.2 T * NaCl (x ion = 0.1;0.2;0.25) KCl (x ion = 0.1;0.2;0.25) RbCl (x ion = 0.1;0.2;0.25) TMD lines a 0.09 0.12 0.15 0.18 0.21 T * tmd -4 -3.5 -3 -2.5 -2 -1.5 -1 S * 2, ex c NaCl (x ion =0.25;0.1) KCl (x ion =0.25;0.1) RbCl (x ion =0.25;0.1) b 0.09 0.12 0.15 0.18 0.21 T * tmd 1.5 2 2.5 C V * NaCl (x ion = 0.1;0.25) KCl (x ion = 0.1;0.25) RbCl (x ion = 0.1;0.25) c Figure 5. (Color online) Temperature of maximum density lines of the NaCl, KCl and RbCl ion models with the soft-core solvent in the T¤¡½¤ plane at different concentration of ions (panel a). The behavior of the pair contribution to excess entropy (panel b) and of reduced heat capacity (panel c) on the temperature of maximum density of the NaCl, KCl and RbCl models at different concentration of ions. changes from NaCl salt to KCl salt. The difference between KCl and RbCl cases is small as expected, be- cause the cation diameter does not changemuch. It seems that the shift of the TMD line to higher densities is determined by the amount of excluded volume added to the system. On the other hand, the enhanced excluded volume effects also lead to a lower extremal point along the temperature axis of the TMD lines. However, this shift of TMD lines on temperature axis is to a great extent effected by the electrostatic inter- actions between ions. If the ion concentration increases for each salt model, the lines substantially shift to lower temperatures. Unfortunately, the solvent in question (due to a set of parameters assumed) does not exhibit vapor – liquid criticality. Therefore, we do not have a natural lower boundary for the TMD lines in, for example, the T ¤ ¡½ ¤ projection. Our preliminary calculations have shown that quite short TMD lines still exist for x ion Æ 0.31 and for 0.32 for NaCl model. The shifts of the TMD lines along the density and temperature axes given in figure 5 (a) are also reflected by the behavior of excess pair entropy and of the reduced heat capacity, see figures 5 (b) and 5 (c), respectively. Shrinking of these curves with an increasing ion concentration is quite evident. Trends of the behavior of S¤ 2,ex and of C¤ V are similar to those discussed for the model NaCl salt in the soft-core solvent just above. To conclude, in the present short communication we considered one extension of our very recent work [2] concerning a simple model mixture of ions and soft-core solvent fluid by using Monte Carlo 43607-7 B. Hribar-Lee, O. Pizio computer simulations. Our extension consists in considering ions of non-equal diameters and different from the diameter of solvent particles. We were interested to study how the presence of ion solutes af- fects thermodynamic properties of the system compared to the single-component fluid, specifically the temperature of maximum density line. For this purpose, we studied three model solute salts, namely the model for NaCl, KCl and RbCl. Our results confirm that the mixtures based on a soft-core solvent in question preserve this anomalous behavior at x ion Æ 0.1 0.2, 0.25, and 0.3, if the cation is small, e.g., for NaCl model. However, for the model with a big cation, e.g., RbCl the density anomaly disappears at x ion around 0.27. The TMD lines were constructed from the simulation data for a set of isochores. They shift to lower temperatures and higher densities with an increasing ion concentration in the T ¤ ¡½ ¤ plane. The TMD line shrinks and moves to lower values of the reduced temperature with an increasing ion concentration for each model salt. These observations agree with the results of molecular dynamics simulations [32] for much more elaborate model for aqueous solutions. We have seemingly made interesting observa- tions about the dependence of the excess pair entropy and constant volume heat capacity along the TMD line at different ion concentrations. Still it is necessary to complement our study by simulations of un- charged solutes with different characteristics mixed with applied or another soft-core model for the sol- vent subsystem. It is known that in the presence of uncharged, non-polar solutes, the curves describing the anomalies shift as well, see e.g. [18]. However, for example, the TMD line changes essentially depend on the solute concentration and strength of solute-solvent interactions. In order to compare our present observations with the models of uncharged solutes in soft-core solvents, additional simulations are nec- essary. This project is in progress in our laboratory at present. Intuitively, one should expect the most pronounced differences between the effects of ions and non-polar solutes if the ion concentration is low or, in other words, when the screening of ion-ion interaction is not well pronounced. There remainsmuch room for extensions of the present work. In order to validate the results concern- ing the trends for anomalies with increasing ion concentrations, it is worth to extend the calculations and establish the solubility limits. It is necessary to explore how the modification of the cross solute-solvent interaction potentials would affect the TMD lines. Another important issue is to investigate various prop- erties of ionic solutions and their peculiarities by using a solvent model characterized by the presence of liquid – vapor phase transition and possibly liquid-liquid separation, e.g., [42]. Theoretical developments are in fact desirable in many of these aspects. Acknowledgements O. P. has been supported in part by European Union IRSES Grant STCSCMBS 268498. B. H.-L. acknowl- edges financial support of the Slovenian Research Agency (ARRS) under grants P1-0201, and J1-4148. O. P. is grateful to David Vazquez for technical support of this work. References 1. Yukhnoskii I.R., Holovko M.F., Microscopic theory of classical equilibrium systems, Naukova Dumka, Kiev, 1980 (in Russian). 2. Lukšič M., Hribar-Lee B., Vlachy V., Pizio O., J. Chem. Phys., 2012, 137, 244502; doi:10.1063/1.4772582. 3. Hemmer P.C., Stell G., Phys. Rev. Lett., 1970, 24, 1284; doi:10.1103/PhysRevLett.24.1284. 4. Jagla E.A., Phys. Rev. E, 1998, 58, 1478; doi:10.1103/PhysRevE.58.1478. 5. Jagla E.A., J. Chem. Phys., 1999, 110, 451; doi:10.1063/1.478105. 6. Jagla E.A., J. Chem. Phys., 1999, 111, 8980; doi:10.1063/1.480241. 7. Jagla E.A., Phys. Rev. E, 2001, 63, 061501; doi:10.1103/PhysRevE.63.061501. 8. Jagla E.A., Phys. Rev. E, 2001, 63, 061509; doi:10.1103/PhysRevE.63.061509. 9. Franzese G., J. Molec. Liq., 2007, 136, 267; doi:10.1016/j.molliq.2007.08.021. 10. Vilaseca P., Franzese G., J. Chem. Phys., 2010, 133, 084507; doi:10.1063/1.3463424. 11. Barros de Oliveira A., Netz P.A., Colla T., Barbosa M.C., J. Chem. Phys., 2006, 124, 084505; doi:10.1063/1.2168458. 12. Barros de Oliveira A., Netz P.A., Colla T., Barbosa M.C., J. Chem. Phys., 2006, 125, 124503; doi:10.1063/1.2357119. 13. Barros de Oliveira A., Barbosa M.C., Netz P.A., Physica A, 2007, 386, 744; doi:10.1016/j.physa.2007.07.015. 43607-8 http://dx.doi.org/10.1063/1.4772582 http://dx.doi.org/10.1103/PhysRevLett.24.1284 http://dx.doi.org/10.1103/PhysRevE.58.1478 http://dx.doi.org/10.1063/1.478105 http://dx.doi.org/10.1063/1.480241 http://dx.doi.org/10.1103/PhysRevE.63.061501 http://dx.doi.org/10.1103/PhysRevE.63.061509 http://dx.doi.org/10.1016/j.molliq.2007.08.021 http://dx.doi.org/10.1063/1.3463424 http://dx.doi.org/10.1063/1.2168458 http://dx.doi.org/10.1063/1.2357119 http://dx.doi.org/10.1016/j.physa.2007.07.015 Density anomaly. Monte Carlo simulations 14. Chaimovich A., Shell M.S., Phys. Chem. Chem. Phys., 2009, 11, 1901; doi:10.1039/B818512C. 15. Barraz N.M., Salcedo E., Barbosa M.C., J. Chem. Phys., 2009, 131, 094504; doi:10.1063/1.3213615. 16. Barraz N.M., Salcedo E., Barbosa M.C., J. Chem. Phys., 2011, 135, 104507; doi:10.1063/1.3630941. 17. Salcedo E., Barros de Oliveira A., Barraz N.M., Chakravarty C., Barbosa M.C., J. Chem. Phys., 2011, 135, 044517; doi:10.1063/1.3613669. 18. Chatterjee S., Ashbaugh H.S., Debenedetti P.G., J. Chem. Phys., 2005, 123, 164503; doi:10.1063/1.2075127. 19. Chatterjee S., Debenedetti P.G., J. Chem. Phys., 2006, 124, 154503; doi:10.1063/1.2188402. 20. Buldyrev S.V., Kumar P., Debenedetti P.G., Rossky P.J., Stanley H.E., Proc. Nat. Acad. Sci. USA, 2007, 104, 20177; doi:10.1073/pnas.0708427104. 21. Maiti M., Weiner S., Buldyrev S.V., Stanley H.E., Sastry S., J. Chem. Phys., 2012, 136, 044512; doi:10.1063/1.3677187. 22. Corradini D., Gallo P., Buldyrev S.V., Stanley H.E., Phys. Rev. E, 2012, 85, 051503, doi:10.1103/PhysRevE.85.051503. 23. Su Z., Buldyrev S.V., Debenedetti P.G., Rossky P.J., Stanley H.E., J. Chem. Phys., 2012, 136, 044511; doi:10.1063/1.3677185. 24. Egorov S.A., J. Chem. Phys., 2011, 134, 234509; doi:10.1063/1.3602217. 25. Pizio O., Dominguez H., Duda Yu., Sokołowski S., J. Chem. Phys., 2009, 130, 174504; doi:10.1063/1.3125930. 26. Lukšič M., Hribar-Lee B., Pizio O., Molec. Phys., 2011, 109, 893; doi:10.1080/00268976.2011.558029. 27. Archer D.G., Carter R.W., J. Phys. Chem. B, 2000, 104, 8563, doi:10.1021/jp0003914. 28. Carter R.W., Archer D.G., Phys. Chem. Chem. Phys., 2000, 2, 5138, doi:10.1039/B006232O. 29. Mishima O., J. Chem. Phys., 2005, 123, 154506; doi:10.1063/1.2085144. 30. Mishima O., J. Chem. Phys., 2007, 126, 244507, doi:10.1063/1.2743434. 31. Corradini D., Gallo P., Rovere M., J. Phys.: Conf. Ser., 2009, 177, 012003; doi:10.1088/1742-6596/177/1/012003. 32. Corradini D., Gallo P., Rovere M., J. Phys.: Condens. Matter, 2010, 22, 284104; doi:10.1088/0953-8984/22/28/284104. 33. Holzmann J., Ludwig R., Geiger A., Paschek D., Angew. Chem. Int. Ed., 2007, 46, 8907, doi:10.1002/anie.200702736. 34. Nunes da Silva J., Salcedo E., Barros de Oliveira A., Barbosa M.C., J. Chem. Phys., 2010, 133, 244506, doi:10.1063/1.3511704. 35. Allen M.P., Tildesley D.J., Computer Simulation of Liquids, Clarendon Press, Oxford, 1987. 36. Frenkel D., Smit B., Understanding Molecular Simulation: From Algorithms to Applications, Academic Press, San Diego, 2002. 37. Balbuena P., Johnston K.P., Rossky P.J., J. Phys. Chem., 1996, 100, 2706, doi:10.1021/jp952194o. 38. Mountain R.D., Int. J. Thermophys., 2007, 28, 536; doi:10.1007/s10765-007-0154-6. 39. Hansen J.-P., McDonald I.R., Theory of Simple Liquids, Elsevier, London, 2006. 40. Pond M.J., Krekelberg W.P., Shen V.K., Errington J.R., Truskett T.M., J. Chem. Phys., 2009, 131, 161101; doi:10.1063/1.3256235. 41. Errington J.R., Truskett T.M., Mittal J., J. Chem. Phys., 2006, 125, 244502; doi:10.1063/1.2409932. 42. Lomba E., Almarza N.G., Martín C., McBride C., J. Chem. Phys., 2007, 126, 244510; doi:10.1063/1.2748043. 43607-9 http://dx.doi.org/10.1039/B818512C http://dx.doi.org/10.1063/1.3213615 http://dx.doi.org/10.1063/1.3630941 http://dx.doi.org/10.1063/1.3613669 http://dx.doi.org/10.1063/1.2075127 http://dx.doi.org/10.1063/1.2188402 http://dx.doi.org/10.1073/pnas.0708427104 http://dx.doi.org/10.1063/1.3677187 http://dx.doi.org/10.1103/PhysRevE.85.051503 http://dx.doi.org/10.1063/1.3677185 http://dx.doi.org/10.1063/1.3602217 http://dx.doi.org/10.1063/1.3125930 http://dx.doi.org/10.1080/00268976.2011.558029 http://dx.doi.org/10.1021/jp0003914 http://dx.doi.org/10.1039/B006232O http://dx.doi.org/10.1063/1.2085144 http://dx.doi.org/10.1063/1.2743434 http://dx.doi.org/10.1088/1742-6596/177/1/012003 http://dx.doi.org/10.1088/0953-8984/22/28/284104 http://dx.doi.org/10.1002/anie.200702736 http://dx.doi.org/10.1063/1.3511704 http://dx.doi.org/10.1021/jp952194o http://dx.doi.org/10.1007/s10765-007-0154-6 http://dx.doi.org/10.1063/1.3256235 http://dx.doi.org/10.1063/1.2409932 http://dx.doi.org/10.1063/1.2748043 B. Hribar-Lee, O. Pizio Аномалiя густини у системi заряджених твердих сфер рiзних розмiрiв у сумiшi з розчинником у моделi з м’яким кором. Результати моделювання Монте Карло Б. Грiбар-Лi1, О. Пiзiо 2 1 Факультет хiмiї i хiмiчної технологiї, Унiверситет м. Любляна, 1000 Любляна, Словенiя 2 Iнститут хiмiї, Нацiональний автономний унiверситет м. Мехiко, Мехiко, Мексика Недавно було дослiджено вплив заряджених твердих сфер однакового розмiру як розчиненої речовини на структурнi i термодинамiчнi аномалiї системи, в якiй застосовано модель плину з потенцiалом з м’яким кором як розчинник. Застосовувався метод моделювання Монте Карло та теорiя iнтегральних рiвнянь [J. Chem. Phys., 2012, 137, 244502]. Метою даної роботи є доповнити попереднє досл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я густини 43607-10 Introduction Model and method Results and discussion

Density anomaly of charged hard spheres of different diameters in a mixture with core-softened model solvent. Monte Carlo simulation results Author(s):

Institution

Ähnliche Einträge