The liquid-vapor interface of the restricted primitive model of ionic fluids from a density functional approach

We investigate the liquid-vapor interface of the restricted primitive model for an ionic fluid using a density functional approach. The applied theory includes the elect rostatic contribution to the free energy functional arising from the bulk energy equation of state and the mean spherical appro...

Повний опис

Збережено в:

Бібліографічні деталі
Дата:	2011
Автори:	Patrykiejew, A., Sokołowski, S., Pizio, O.
Формат:	Стаття
Мова:	English
Опубліковано:	Інститут фізики конденсованих систем НАН України 2011
Назва видання:	Condensed Matter Physics
Онлайн доступ:	http://dspace.nbuv.gov.ua/handle/123456789/119941
Теги:	Додати тег Немає тегів, Будьте першим, хто поставить тег для цього запису!
Назва журналу:	Digital Library of Periodicals of National Academy of Sciences of Ukraine
Цитувати:	The liquid-vapor interface of the restricted primitive model of ionic fluids from a density functional approach / A. Patrykiejew, S. Sokołowski, O. Pizio // Condensed Matter Physics. — 2011. — Т. 14, № 1. — С. 13603: 1-12. — Бібліогр.: 48 назв. — англ.

Репозитарії

Digital Library of Periodicals of National Academy of Sciences of Ukraine

id	irk-123456789-119941
record_format	dspace
spelling	irk-123456789-1199412017-06-11T03:04:37Z The liquid-vapor interface of the restricted primitive model of ionic fluids from a density functional approach Patrykiejew, A. Sokołowski, S. Pizio, O. We investigate the liquid-vapor interface of the restricted primitive model for an ionic fluid using a density functional approach. The applied theory includes the elect rostatic contribution to the free energy functional arising from the bulk energy equation of state and the mean spherical approximation for a restricted primitive model, as well as the associative contribution, due to the formation of pairs of ions. We compare the density profiles and the values of the surface tension with previous theoretical approaches. Ми досл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 The liquid-vapor interface of the restricted primitive model of ionic fluids from a density functional approach / A. Patrykiejew, S. Sokołowski, O. Pizio // Condensed Matter Physics. — 2011. — Т. 14, № 1. — С. 13603: 1-12. — Бібліогр.: 48 назв. — англ. 1607-324X PACS: 68.08.-p, 68.43.Fg, 82.35.Gh, 68.43.-h DOI:10.5488/CMP.14.13603 arXiv:1106.3236 http://dspace.nbuv.gov.ua/handle/123456789/119941 en Condensed Matter Physics Інститут фізики конденсованих систем НАН України
institution	Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection	DSpace DC
language	English
description	We investigate the liquid-vapor interface of the restricted primitive model for an ionic fluid using a density functional approach. The applied theory includes the elect rostatic contribution to the free energy functional arising from the bulk energy equation of state and the mean spherical approximation for a restricted primitive model, as well as the associative contribution, due to the formation of pairs of ions. We compare the density profiles and the values of the surface tension with previous theoretical approaches.
format	Article
author	Patrykiejew, A. Sokołowski, S. Pizio, O.
spellingShingle	Patrykiejew, A. Sokołowski, S. Pizio, O. The liquid-vapor interface of the restricted primitive model of ionic fluids from a density functional approach Condensed Matter Physics
author_facet	Patrykiejew, A. Sokołowski, S. Pizio, O.
author_sort	Patrykiejew, A.
title	The liquid-vapor interface of the restricted primitive model of ionic fluids from a density functional approach
title_short	The liquid-vapor interface of the restricted primitive model of ionic fluids from a density functional approach
title_full	The liquid-vapor interface of the restricted primitive model of ionic fluids from a density functional approach
title_fullStr	The liquid-vapor interface of the restricted primitive model of ionic fluids from a density functional approach
title_full_unstemmed	The liquid-vapor interface of the restricted primitive model of ionic fluids from a density functional approach
title_sort	liquid-vapor interface of the restricted primitive model of ionic fluids from a density functional approach
publisher	Інститут фізики конденсованих систем НАН України
publishDate	2011
url	http://dspace.nbuv.gov.ua/handle/123456789/119941
citation_txt	The liquid-vapor interface of the restricted primitive model of ionic fluids from a density functional approach / A. Patrykiejew, S. Sokołowski, O. Pizio // Condensed Matter Physics. — 2011. — Т. 14, № 1. — С. 13603: 1-12. — Бібліогр.: 48 назв. — англ.
series	Condensed Matter Physics
work_keys_str_mv	AT patrykiejewa theliquidvaporinterfaceoftherestrictedprimitivemodelofionicfluidsfromadensityfunctionalapproach AT sokołowskis theliquidvaporinterfaceoftherestrictedprimitivemodelofionicfluidsfromadensityfunctionalapproach AT pizioo theliquidvaporinterfaceoftherestrictedprimitivemodelofionicfluidsfromadensityfunctionalapproach AT patrykiejewa liquidvaporinterfaceoftherestrictedprimitivemodelofionicfluidsfromadensityfunctionalapproach AT sokołowskis liquidvaporinterfaceoftherestrictedprimitivemodelofionicfluidsfromadensityfunctionalapproach AT pizioo liquidvaporinterfaceoftherestrictedprimitivemodelofionicfluidsfromadensityfunctionalapproach
first_indexed	2025-07-08T16:56:50Z
last_indexed	2025-07-08T16:56:50Z
_version_	1837098667019862016
fulltext	Condensed Matter Physics 2011, Vol. 14, No 1, 13603: 1–12 DOI:10.5488/CMP.14.13603 http://www.icmp.lviv.ua/journal The liquid-vapor interface of the restricted primitive model of ionic fluids from a density functional approach A. Patrykiejew1, S. Sokołowski1, O. Pizio2∗ 1 Department for the Modelling of Physico-Chemical Processes, Maria Curie-Skłodowska University, 20–031 Lublin, Poland, 2 Instituto de Quı́mica de la UNAM, Coyoacán 04510, México, D.F., Mexico Received September 13, 2010, in final form November 23, 2010 We investigate the liquid-vapor interface of the restricted primitive model for an ionic fluid using a density functional approach. The applied theory includes the electrostatic contribution to the free energy functional arising from the bulk energy equation of state and the mean spherical approximation for a restricted primitive model, as well as the associative contribution, due to the formation of pairs of ions. We compare the density profiles and the values of the surface tension with previous theoretical approaches. Key words: density functional, adsorption, chains, crystals PACS: 68.08.-p, 68.43.Fg, 82.35.Gh, 68.43.-h 1. Introduction According to the simplest model of ionic fluids, called the restricted primitive model (RPM), the fluid consists of charged hard spheres (at total density ρ) of equal diameter, d, half of which carry a charge of +q, while the other half carry a charge of −q. These spheres are immersed in a dielectric medium of dielectric constant ǫ. Reduced thermodynamic quantities, temperature and density, appropriate to this model are defined as follows T ∗ = kT ǫd/q2 and ρ∗ = ρd3. Perhaps, the first successful statistical-mechanical description of this model was proposed by Debye and Hückel [1] about ninety years ago. The Debye-Hückel approach predicted the existence of the first-order transition for RPM and was subsequently applied in several works, e.g., to describe phase behavior of three-component ionic fluids [2]. More recently, the bulk RPM was studied by means of the Ornstein-Zernike (OZ) equation supplemented by the mean spherical approximation (MSA) [3–5], which permits to obtain an analytical solution for the pair correlation functions. The MSA correlation functions can be used to evaluate an equation of state [6, 7]. In spite of quite reasonable results obtained for the structure and thermodynamics in a general sense, neither virial nor compressibility equations of state permit to obtain gas-liquid transition for the RPM, which contradicts the computer simulations [8–15]. By contrast, the equation of state resulting from the energy route predicts the first-order phase transition between an ionic vapor and a dense ionic liquid [6, 7]. However, there is a big discrepancy between the values of critical temperature and critical density of MSA theory and the critical parameters resulting from computer simulations. Moreover, a more sophisticated generalized mean-spherical approximation [16], which was designed to reconcile thermodynamic inconsistency of the MSA, does not lead to a better agreement of the critical parameters with simulations as well. Since both MSA and generalized MSA theories underestimate the critical density and overes- timate the critical temperature, attempts were made to include ad hoc the effect of ion pairing (association) into the RPM, which is most pronounced along the vapor branch of the coexistence envelope and near the critical point [13, 17]. The inclusion of the association into the theory improves the agreement of theoretical predictions for the critical constants with simulations (cf. ∗E-mail: pizio@servidor.unam.mx c© A. Patrykiejew, S. Sokołowski, O. Pizio 13603-1 http://dx.doi.org/10.5488/CMP.14.13603 http://www.icmp.lviv.ua/journal A. Patrykiejew, S. Sokołowski, O. Pizio table I of [13]). However, the overall shape of the coexistence envelope coming from several theoreti- cal approaches is not well described compared to simulation data. In particular, the liquid densities along coexistence are poorly described for different temperatures. We should bear in mind that the Hamiltonian of the model with association used in the theoretical developments differs from the genuine RPM used in simulations. Recent developments in the modelling of different types of associations and clustering in liquids and solutions have been comprehensively and critically reviewed by Holovko in [18]. The liquid-vapor interface is one of the simplest examples of nonuniform systems. However, compared to the case of Lennard-Jones fluids [19], much less is known about the liquid-vapor interface in the RPM. Theoretical work on the gas-liquid interface for RPM was pioneered by Telo da Gama et al. [20], that used the gradient expansion method [21]. The problem of describing an ionic liquid-ionic vapor interface has been also studied by Groh et al. [22], who proposed a density functional approach to evaluate the surface tension and the density profiles through the interface. Their approach involved a local density approximation for the hard-sphere part of the free-energy functional and a nonlocal treatment of Coulombic contributions. The latter term was evaluated approximating the inhomogeneous pair correlation functions by their bulk counterparts. There are significant differences between the theory used by Groh et al. [22] and the density functional approaches that were proved successful in the description of electrical double layers at charged walls [23–31]. Namely, the theories of [23–27] treat the electrostatic part of the free energy by means of a second-order density expansion about the density of a reference fluid, which was taken to be the homogeneous bulk fluid far from the surface. Although such an approach can be appropriate for many purposes, but it is problematical when it comes to a liquid-vapor interface, where two bulk fluids are involved. Moreover, it is known that in the case of Lennard-Jones fluids, the corresponding second-order expansion with respect to a homogeneous fluid fails to account for the liquid-vapor coexistence. More recently, the description of a liquid-vapor interface was the subject of investigations of Weiss and Schröer [32]. Using a square-gradient type theory they computed the density profiles and the interfacial tension at different temperatures using Debye- Hückel theory and its extension to ion-pair formation, as well as adding the free energy term describing correlations between ion pairs as entities and free ions, as developed by Fisher and Levin [17]. Any density functional theory determines thermodynamic properties of an inhomogeneous fluid from the Helmholtz free energy F and its functional dependence on the local densities {ρi(r)}. The free energy functional is commonly decomposed into the sum of three contributions, namely into the ideal, the hard-sphere, and the electrostatic terms. Of course, one can also include here an additional term, due to possible association of ions. Various formulations of the DFT have been discussed in the literature. As we have already stressed, a majority of the DFTs have followed perturbational second-order expansion of the electrostatic free energy functional with respect to a bulk homogeneous fluid [23–27]. Recently, Gillespie et al. [33, 34] proposed a version of the electrostatic free energy functional that replaces a uniform reference system with a suitably chosen position-dependent reference fluid. The inhomogeneous reference fluid densities are then computed from the local densities by a self- consistent iteration procedure. Actually, Gillespie et al. [33, 34] proposed a reference fluid density functional which permits to construct a reference model that locally satisfies electroneutrality con- dition and has the same ionic strength at every point as the inhomogeneous fluid in question. Such a construction of the reference fluid permits to apply the expression for the electrostatic contribu- tion to the free energy, which results from the equation of state for a bulk ionic system and makes it unnecessary to employ the direct correlation functions as input quantities. This kind of approach was proposed by us to study liquid-vapor transitions in RPMs confined to slit-like pores [31, 35]. The results of our approach reasonably well reproduce the structure of ionic fluids at a charged wall at sufficiently high temperatures and correctly predict [36] the temperature dependence of the double layer capacitance, in agreement with experiments and simulation data [37]. Quite re- cently, the theory was also extended to include the effects of association between unlike ions [38]. However, the applications of the approach described above would not be complete without explor- 13603-2 The liquid-vapor interface of the restricted primitive model ing liquid-vapor interface of ionic fluid. Therefore, in this communication we intend to apply the theory of [35, 36, 38] to the study of liquid-vapor interface of the RPM. We calculate the density profiles and the interfacial tension. Moreover, we would like to investigate the effect of association leading to the formation of ion pairs and whether this type of effects improves the description of the interfacial properties similarly to the bulk structure and thermodynamics. 2. Theory We consider a binary mixture of ionic species. The symbols d, Zi = qi/e and µi denote, respec- tively, the hard-sphere diameter, valence of ions and the chemical potential of species i = 1, 2, and e is the electron charge. The interaction between the ions is uij(r) = { ∞, r < 1, e2ZiZj/εr, r > 1, (1) where r is the distance between a pair of ions. We also assume that Z1 = −Z2 = 1 and that the dielectric constant ε is uniform throughout the entire system. The DFT we use in this work is identical to that described in [35, 38] and therefore we present only its brief description. If there is no external potential field the grand potential of the system is written in the form, Ω = ∫ drf [{ρi}]− ∑ i=1,2 ∫ ρi(z)µidr. (2) According to the usual density functional treatment, the free energy density functional, f [{ρi}], is decomposed into ideal, hard sphere, electrostatic and associative terms f [{ρi}] = fid[{ρi}] + fhs[{ρi}] + fel[{ρi}] + fas[{ρi}]. The ideal term is fid[{ρi}] = ∑ i=1,2[ρi(z) ln ρi(z)− ρi(z)], whereas for the hard-sphere contribution we apply an expression resulting from a recent version [39] of the Fundamental Measure Theory [40, 41], with the free energy density consisting of terms dependent on scalar and on vector weighted densities, nα(r) (α = 0, 1, 2, 3) and nα(r) (α = V 1, V 2) [39], f(hs) = −n0 ln(1− n3) + n1n2 − nV 1 · nV 2 1− n3 + n3 2(1− 3ξ2) n3 + (1− n3) 2 ln(1 − n3) 36πn2 3(1 − n3)2 , (3) where ξ(r) = \|nV 2(r)\|/n2(r). The definitions of weighted densities, nα(r), α = 0, 1, 2, 3, V 1, V 2, are given in [40, 41]. The electrostatic contribution is evaluated from the approach described in [35, 38], according to which fel[{ρi}]/kT = − e2 T ∗ [Z2 1 ρ̄1 + Z2 2 ρ̄2] Γ({ρ̄i}) 1 + Γ({ρ̄i})d + Γ3({ρ̄i}) 3π , (4) where ρ̄i(z) denote suitably defined inhomogeneous average densities of a reference fluid. The form for fel[{ρ̄i(z)}] results from the MSA equation of state evaluated via the energy route [3, 4, 6– 8, 13]. In the above Γ({ρ̄i}) = ( √ 1 + 2κ({ρ̄i})d− 1)/2d, where κ({ρ̄i}) denotes the inverse Debye screening length, κ2({ρ̄i}) = (4πe2/T ∗)[Z2 1 ρ̄1 + Z2 2 ρ̄2]. The reference fluid densities {ρ̄i} are evaluated by employing the approach of Gillespie et al. [33, 34]. In the case of gas-liquid interface, the spatial symmetry of the RPM implies that the density profiles of the two ionic species should be the same. Therefore, in the presently considered situation, the approach of Gillespie et al. [33, 34] is simplified and the reference fluid densities are given by ρ̄i(z) = ∫ ρi(z)Wi(\|r− r ′\|)dr′, (5) where the weight function Wi(r) is just a normalized step function Wi(\|r− r ′\|) = θ(\|r − r ′\| −Rf (r ′)) (4π/3)R3 f (r ′) . (6) 13603-3 A. Patrykiejew, S. Sokołowski, O. Pizio The radius of the sphere over which the averaging is performed, Rf , is approximated by the “capacitance” radius, i.e., by the ion radius plus the screening length Rf (r) = d 2 + 1 2Γ({ρ̄i(r)}) . (7) The evaluation of Rf requires an iteration procedure. This iteration loop should be carried out in addition to the main iteration procedure used to evaluate the density profiles, see for details [33, 34]. Finally, the associative contribution, fas, is formulated at the level of the first-order thermody- namic perturbation theory [38, 42] fas/kT = ∑ i=1,2 ρ̄i [ lnα({ρ̄i}) + 1 2 − α({ρ̄i}) 2 ] , (8) where α is the degree of dissociation according to the mass action law, 2α = 1−K(ρ̄1 + ρ̄2)α 2. The association constant K ≡ K({ρ̄i}, T ) is a product of solely temperature dependent term, K0 and Kγ , K = K0Kγ . The constant K0, is [43] K0 = 8πd3 ∞ ∑ m=2 (T ∗)−2m (2m)!(2m− 3) , (9) whereas Kγ follows from the so-called simple interpolation scheme [13, 38, 44] Kγ = (1− η̄/2) (1− η̄)3 exp [ − Γd T ∗ (2 + Γd) (1 + Γd)2 ] . (10) The theory presented above uses different weighted densities to evaluate the hard sphere and electrostatic free energy contributions. Note that the idea of using different weighted densities to evaluate these contributions has also been employed by Patra et al. [25, 26]. The density profile is obtained by minimizing the excess grand potential functional ∆Ω = Ω−Ωb, δ∆Ω δρi(z) = 0, for = 1, 2, (11) where Ωb is the grand potential of a bulk uniform system of density {ρb,i}. The evaluation of the density profiles requires the knowledge of the densities of both coexisting phases. Therefore, we have precisely evaluated the bulk phase diagram prior to the density profile calculations. Next, the local density calculations were carried out assuming that for z 6 −Lz and for > Lz the density of the fluid is equal to the bulk densities of gaseous and liquid phases, respectively. The value of Lz was evaluated at each temperature by testing the convergence of the density profiles towards the final solutions. We have started with Lz = 30d and the consecutive runs were carried out doubling the value of Lz for each new run. All the integrations were performed using Simpson method with the grid size of 0.01d. The convergence criterion definitly states that the maximum difference between two consecutive iterations must be smaller than 1E-8 percent. The knowledge of the density profiles of ions permits to calculate the excess grand thermodynamics potential per unit surface area, A, i.e., the surface tension γ = ∆Ω/A. 3. Results and discussion Figure 1 shows the phase diagram in the density-temperature plane. The phase diagram result- ing from the MSA energy equation of state has been presented in several previous works, cf. [13, 22]. Nevertheless, we have decided to include it here for completeness of the study. The solid line corre- sponds to the system without association (i.e. the associative free energy term is neglected in the grand thermodynamics potential functional), whereas the dashed line describes phase behavior of the model with the association effects included. Inclusion of chemical association into the theory 13603-4 The liquid-vapor interface of the restricted primitive model provides a mechanism for the formation of ionic pairs. The fraction of pairs depends on density, on temperature and on the association constant. From physical point of view, the formation of pairs of oppositely charged ions, alters the effective interactions between all the particles. In effect, the critical temperature of liquid-vapor transition, T ∗ c ≈ 0.0745, becomes lower compared to the criti- cal temperature of the model without association, T ∗ c ≈ 0.0786. The liquid-vapor phase diagrams are strongly asymmetric. At low temperatures the vapor phase is very dilute. For example, for the system without association at T ∗ = 0.05 the vapor density is ρ∗b ≈ 1.094E− 5, whereas the density of the coexisting liquid phase is ρ∗b ≈ 0.2135. 0.0001 0.01 ρ b * 0.06 0.07 0.08 T * Figure 1. Phase diagram resulting from MSA approach and the energy equation of state. Solid and dashed lines denote the results for the model without and with association. -8 -4 0 4 8 z/d 0 0.1 0.2 ρ∗ (z ) T/T c =0.709 T/T c =0.827 T/T c =0.9456 Figure 2. A comparison of the total density profiles, ρ∗(z) = (ρ1(z)+ρ1(z))d3 resulting from the DFT without association theory (solid lines) and from MSA1 theory of Groh et al. [22] (dashed lines). The temperatures T/Tc are given in the figure. First, we consider the system without association. The evaluation of the liquid-vapor density profiles was carried out assuming that the limiting values ρi(−Lz) and ρi(Lz) are equal to the MSA densities of coexisting phases. Of course, the agreement of our results with the results published 13603-5 A. Patrykiejew, S. Sokołowski, O. Pizio previously [22, 32] essentially depends on the agreement between the bulk phase diagrams. One should keep in mind the above remark while analyzing figure 2, where we show a comparison of our results with those of Groh et al.[22]. The latter profiles were obtained from a different bulk theory (called by them “the MSA1 approach” [22]), which yields the value of the critical temperature, T ∗ c ≈ 0.0846, different from that obtained by MSA theory. Note that a comparison of the entire phase diagrams resulting from the MSA and MSA1 approaches in given in figure 1 of [22]. 0 2 4 6 z/d -10 -5 ln [ρ b, il -ρ i(z )] y = -3.9231-0.41785 (z/d) y = -3.0234-1.015(z/d) y = -2.7294-1.3354(z/d) T * = 0.074 T * =0.065 T * =0.06 Figure 3. The convergence of the liquid parts of the single species density profiles towards the bulk liquid equilibrium density at three different temperatures given in the figure. Figure also contains the equations for straight lines (marked by solid lines) approximating the data, computed from density functional theory. Here y = ln[ρ)b, il−ρi(z)] and ρlb,i is the single species density of an ionic liquid at coexistence. The calculations are for the model without association. The way in which the density profile approaches the asymptotic bulk values has been discussed by Groh et al. [22]. We have checked that the decay of the function ρlb,i − ρi(z) (where ρlb,i is the liquid-phase bulk density) with z → Lz is exponential, i.e. ln[ρi(z) − ρlb,i] ∝ z/ξ. This point is illustrated in figure 3, where we have plotted ln[ρlb,i − ρi(z)] versus z/d. At larger distances this dependence is fairly well approximated by a straight line with the correlation coefficient, R2, very close to unity. We now consider the effect of the association, cf. figure 4. Part a shows a comparison of the total density profiles for the systems with (symbols) and without association (lines) at different temperatures. At low temperatures the liquid density for the system with the association included is higher than for the system with the association effects switched off, cf. figure 1. The profiles for the model with association are more “smeared out”, but this is the effect of lowering the critical temperature due to association effects. Similarly to the case of a system without association we have found exponential decay of the function ρlb,i − ρi(z) with z → Lz (the corresponding plot was omitted for the sake of brevity). Figure 4b shows the profiles of α(z) (the phase diagram in the α-T plane has been already presented in [38]) We stress that the interfacial region of α(z) is shifted towards rarefied phase compared to the density profiles. It means that the position of the “pseudo-Gibbs” dividing surface, which can be introduced for the dissociation degree in analogy to the usual liquid-vapor dividing surface is also shifted compared to the position of the usual dividing surface (which is located at z = 0). Note that similar shift of the dissociation degree was also observed in the case of liquid-vapor interfaces of non-ionic fluids [45]. The surface tension obtained from the density functional theory and from the MSA1 and MSA approaches, developed by Groh et al. [22], as well as from the square gradient theory [20] are shown in figure 5. We have plotted here the ratio γd2/kTc versus the temperature reduced by the critical temperature, T/Tc. Except for the temperatures very close to the critical temperature, the present 13603-6 The liquid-vapor interface of the restricted primitive model -9 -6 -3 0 3 6 9 z/d 0 0.1 0.2 ρ* (z ) T=0.055 T=0.060 T=0.065 T=0.070 T=0.075 a -9 -6 -3 0 3 6 9 z/d 0 0.2 0.4 α( z) T=0.070 T=0.065 T=0.060 T=0.055 b Figure 4. Part a. A comparison of the total density profiles for the models without (solid lines) and with association (symbols). Part b. Dissociation degrees across the vapor-liquid interface. The temperatures are given in the figure. 0.4 0.6 0.8 1 T/T c 0 0.05 0.1 0.15 0.2 γd 2 /( kT c) MSA1 MSA square gradient theory DFT DFT+ASS Figure 5. The dependence of the surface tension on the temperature. Curves labelled MSA, MSA1 were evaluated by Groh et al. [22], the curve labelled “square gradient theory” was evaluated by Telo da Gama et al. [20], the curve abbreviated as DFT was obtained for the model without association, whereas black circles – for the model with the association effects included. theory predicts lower values of the surface tension than all the rest approaches. The associative free energy terms still lower the value of γd2/kTc. Note that the MSA theory grossly underestimates the densities of the coexisting liquid phase. This can suggest that the values of the surface tension are also underestimated. On the other hand, the value of the critical temperature is overestimated by the MSA approach, and hence can lead to too high values of the interfacial tension. Near critical temperature, both models (i.e., the model with and without association) yield the standard mean-field critical behavior, i.e., γd2/kTc ∝ (1 − T/Tc) (3/2) (see figure 6). As the temperature approaches the critical point, the theory yields a nearly linear variation of ln[γd2/kTc] with ln[(1 − T/Tc)] and the slope of the straight line approximating the obtained data is almost perfectly equal to (3/2). 13603-7 A. Patrykiejew, S. Sokołowski, O. Pizio -6 -5 -4 -3 -2 -1 0 ln(1-T/T c ) -10 -8 -6 -4 -2 ln (γ d2 /k T c) γd 2 /kT c ∝ (1-T/T c ) (3/2) Figure 6. The decay of the surface tension on approaching the critical temperature. Solid line and empty circles are for the model without association, whereas dashed line and filled circles – for the model with association. The straight lines slope is 3/2. -6 -5 -4 -3 -2 -1 0 ln(1-T/T c ) 0 1 2 3 4 ln (W /d ) W/d ∝ (1−Τ/Τ c ) (-0.5) Figure 7. The divergence of the interfacial width W , on approaching the critical temperature. The straight lines slope is −1/2. Abbreviations are the same as in figure 6. The width of the interfacial zone can be characterized by the parameter W , defined as [46] W = −[ρ(z = Lz)− ρ(z = −Lz)] [ dρ(z) dz ] −1 z=z0 , (12) where z0 is given by ρ(z0) = (1/2)[ρs(z = Lz)+ρ(z = −Lz)] and ρ(z = Lz) and ρ(z = −Lz) are the total densities of the coexisting liquid and vapor phases. Obviously, the interface becomes wider as the temperature increases. The effects of association slightly increases the interface width. As the temperature approaches the critical temperature, the interfacial width diverges. To investigate the character of this divergence, we have plotted the values of lnW/d versus ln(1 − T/Tc), see figure 7. For T approaching the critical temperature, the dependence of lnW/d versus ln(1−T/Tc) is linear. We have found the following scaling: W ∗ ∝ (1 − T/Tc(M))−0.5. The obtained value of 13603-8 The liquid-vapor interface of the restricted primitive model the exponent (0.5) is characteristic of mean-field type theories [22]. The scaling of W is consistent with the scaling of the surface tension. -10 0 10 z/d 0 5 10 ρ (z )/ ρ c T * =0.0391 T * =0.035 DFT DFT+A Figure 8. A comparison of the total density profiles for the models without association (solid and dotted lines) and with association (dashed and dash-dotted lines) with computer simulation data. The simulation results for the profiles (see details in the text) are shown by symbols [47, 48]. The temperatures are given in the figure. Our discussion involved the results of different, but solely theoretical approaches. However, it must be clarified how accurate our approaches are (with and without association) compared to the available computer simulation data. In particular, Alejandre and co-workers investigated the liquid- vapor interface for the RPM using molecular dynamics simulation method supplemented by certain specific technical peculiarities dealing specifically with the RPM [47, 48]. These authors presented the density profiles, ρ∗(z) at four reduced temperatures, T ∗=0.035, 0.038, 0.040 and 0.043 [47]. On the other hand, a single density profile at T ∗ = 0.391 relevant to our study was given in [48] (this profile is for the soft primitive model but with parameters permitting comparison with the RPM). The authors have not performed estimates of the critical parameters coming from their method of study of the RPM. Thus, they somewhat arbitrarily used the critical temperature obtained in [11], T ∗ c = 0.0492, for their purposes. We have rescaled the temperature and density with respect to the critical parameters of both theoretical approaches (with and without association), respectively, and present here a comparison of the density profiles with computer simulation data in figure 8. Note that the critical density given in [11] was used to rescale the simulation data. We observe that the rescaled liquid phase density at coexistence coming from the DFT for the model with ion pairing is in reasonable agreement with the simulation data. The DFT for the model without association greatly overestimates these liquid densities. However, the width of the interfacial region is not described well. Theoretical approaches yield narrower interface width for two reduced temperatures in question. It is difficult to attribute this specific inaccuracy to the particular term of the free energy functional. However, it seems that the mean field consideration requires improvement in order to better describe interparticle correlations and hopefully to obtain a more adequate interface width. As concerns the relationship between the simulation data for the surface tension and the results of theoretical approaches we would like to comment on the following points. A set of computer simulation results for the surface tension were put together with the predictions of some theoretical approaches [20, 22, 32] (that were discussed in the introductory part of our communication) on arbitrary temperature scale, see figure 2 of [47]. Apparently, the theories yield the surface tension of the same magnitude as computer simulation data. This permitted to conclude that the surface tension from this set of theories [20, 22] is overestimated because the effect of ion association is not 13603-9 A. Patrykiejew, S. Sokołowski, O. Pizio taken into account and solely the theory that accounts for association [32] is in good agreement with simulations. Then, it would seem that our approach taking ion pairing into account is the best, see figure 5. Moreover, a comparison of the theories, simulations and experimental data was performed [48] using an arbitrary scale. Such an idyllic picture is, however, entirely destroyed, if one rescales the reduced temperature scale by the critical temperature of each theory and of simulations. This is actually an adequate comparison. According to [47] the simulations yield, e.g. γd2/kTc=0.3 at T/Tc=0.6. Such a value is very much higher than the results of any theory, c.f. figure 5, at this particular temperature. Moreover, this means that the inclusion of association effects makes things worse than from any other theory. This is difficult to accept, however. According to the more recent simulation work of the same authors [48] it is difficult to draw a definite conclusion about the dependence of the surface tension on temperature for the model in question due to the statistical error associated with the value for surface tension. Moreover, it has been discussed [48] that the inflection point at T ∗ = 0.04 (again it is a reduced temperature of simulation) can be present on the dependence of surface tension on temperature as a result of chemical association of ions in the RPM in close similarity to hydrogen bonded fluids. In our opinion, additional work, both theoretical and simulational, is necessary to obtain definite answers about the behavior of surface tension on temperature for the model in question and the related models. From theoretical viewpoint one needs a theory that yields a more adequate shape of the bulk coexistence envelope. Then, the free energy expression coming from this approach, if available, would contribute to the development of better density functional approaches for inhomogeneous ionic fluids. To summarize briefly, in this work we have applied density functional theory to the study of the liquid-vapor interface of a RPM fluid. Of course, the theory we use here is ad hoc. However, our previous calculations [35, 36] have shown that it reproduces reasonably well the structure of a fluid at a charged and uncharged wall, predicts the dependence of the capacitance of the double layer on the temperature and yields “capillary evaporation” phase diagrams for confined ionic systems [35, 36, 38]. We have shown that the inclusion of ion pairing effect leads to a better agreement of the density profiles of ions across the liquid-vapor interface. A more sophisticated approaches [13], compared to the present one, including the effects of ionic association can lead to even better agreement of the coexistence envelope with simulations and of the density profiles across the liquid-vapor interface. We plan to extend our study in this respect in the nearest future. However, the inaccuracies of the values for the liquid phase density at coexistence at different temeperatures prevent us from obtaining a reasonable agreement with the available simulation data for the surface tension. It seems, however, that a more ample set of simulation data using different techniques is necessary. Nevertheless, the theories of this study yield a decreasing surface tension with increasing temperature as one intuitively would expect. Still we are not able to conclude that the theoretical approach is entirely successful both for the fluid-solid and fluid-fluid interfaces and look forward to improving it in future work. References 1. Debye P.W., Hückel E., Phys. Z., 1923, 24, 185. 2. Moreira A.G., Netz R.R., European Phys. J. D, 2001, 13, 61; doi:10.1007/s100530170287. 3. Blum L., J. Chem. Phys., 1974, 61, 2129; doi:10.1063/1.1682224. 4. Blum L., Molec. Phys., 1975, 30, 1529; doi:10.1080/00268977500103051. 5. Sánchez-Diaz L.E., Vizcarra-Rendón A., Medina-Noyola M., J. Chem. Phys., 2010, 132, 234506; doi:10.1063/1.3455336. 6. Waisman E., Lebowitz J.L., J. Chem. Phys., 1970, 52, 4307; doi:10.1063/1.1673642. 7. Waisman E., Lebowitz J.L., J. Chem. Phys., 1972, 56, 3093; doi:10.1063/1.1677645. 8. Panagiotopoulos A.Z., Fluid Phase Equilibria, 1992, 76, 97; doi:10.1016/0378-3812(92)85080-R. 9. Caillol J.N., J. Chem. Phys., 1994, 100, 2161; doi:10.1063/1.466513. 10. Romero-Enrique J.M., Orkoulas G., Panagiotopoulos A.Z., Fisher M.E., Phys. Rev. Lett., 2000, 85, 4558; doi:10.1103/PhysRevLett.85.4558. 11. Yan Q.L., de Pablo J.J., J. Chem. Phys., 2001, 114, 1727; doi:10.1063/1.1335653. 12. Orkoulas G., Panagiotopouloas A.Z., J. Chem. Phys., 1994, 101, 1452; doi:10.1063/1.467770. 13603-10 http://dx.doi.org/10.1007/s100530170287 http://dx.doi.org/10.1063/1.1682224 http://dx.doi.org/10.1080/00268977500103051 http://dx.doi.org/10.1063/1.3455336 http://dx.doi.org/10.1063/1.1673642 http://dx.doi.org/10.1063/1.1677645 http://dx.doi.org/10.1016/0378-3812(92)85080-R http://dx.doi.org/10.1063/1.466513 http://dx.doi.org/10.1103/PhysRevLett.85.4558 http://dx.doi.org/10.1063/1.1335653 http://dx.doi.org/10.1063/1.467770 The liquid-vapor interface of the restricted primitive model 13. Jiang J., Blum L., Bernard O., Prausnitz J.M., Sandler S.I., J. Chem. Phys., 2002, 116, 7977; doi:10.1063/1.1468638. 14. Valeriani C., Camp P.J., Zwanikken J.W., van Roij R., Dijkstra M., J. Phys. Condens. Mat., 2010, 22, 104122; doi:10.1088/0953-8984/22/10/104122. 15. Diehl A., Panagiotopoulos, A.Z., J. Chem. Phys., 2006, 124, 194509; doi:10.1063/1.2192498. 16. Stell G., J. Phys.: Condens. Mat., 1996, 8, 9329; doi:10.1088/0953-8984/8/47/024. 17. Fisher M., Levin Y., Phys. Rev. Lett., 1993, 71, 3826; doi:10.1103/PhysRevLett.71.3826. 18. Holovko M.F., J. Mol. Liq., 2002, 96–97, 65; doi:10.1016/S0167-7322(01)00327-0. 19. Rowlinson J.S., Widom B., Molecular Theory of Capillarity, Claredon, Oxford, 1982. 20. Telo da Gama M.M., Evans R., Sluckin T.J., Molec. Phys., 1980, 41, 1355; doi:10.1080/00268978000103591. 21. Evans R., Sluckin T.J., Molec. Phys., 1980, 40, 413; doi:10.1080/00268978000101581. 22. Groh B., Evans R., Dietrich S., Phys. Rev. E, 1998, 57, 6944; doi:doi:10.1103/PhysRevE.57.6944. 23. Mier y Teran L., Suh S.H., White H.S., Davis H.T., J. Chem. Phys., 1990, 92, 5087; doi:10.1063/1.458542. 24. Kierlik E., Rosinberg M.L., Phys. Rev. A, 1991, 44, 5025; doi:10.1103/PhysRevA.44.5025. 25. Patra C.N., Ghosh S.K., Phys. Rev. E, 1993, 47, 4088; doi:10.1103/PhysRevE.47.4088. 26. Patra C.N., J. Chem. Phys., 1999, 111, 9832; doi:10.1063/1.480319 . 27. Boda D., Henderson D., Mier y Teran L., Soko lowski S., J. Phys. Condens. Mat., 2002, 14, 11945; doi:10.1088/0953-8984/14/46/305. 28. Pizio O., Bucior K., Patrykiejew A., Soko lowski S., J. Chem. Phys., 2005, 123, 214902; doi:10.1063/1.2128701. 29. Goel T., Patra C.N., Ghosh S.K., Mukherjee T., Molec. Phys., 2009, 107, 19; doi:10.1080/00268970802680497. 30. Goel T., Patra C.N., Ghosh S.K., Mukherjee T., J. Chem. Phys., 2008, 129, 154906; doi:10.1063/1.2992525. 31. Pizio O., Patrykiejew A., Sokolowski S., Condens. Matter Phys., 2004, 7, 779. 32. Weiss C.V., Schröer W., J. Phys. Condens. Mat., 2000, 12, 2637; doi:10.1088/0953-8984/12/12/306. 33. Gillespie D., Nonner W., Eisenberg R.S., J. Phys. Condens. Mat., 2002, 14, 12129; doi:10.1088/0953-8984/14/46/317. 34. Gillespie D., Nonner W., Eisenberg R.S., Phys. Rev. E, 2003, 68, 031503; doi:10.1103/PhysRevE.68.031503. 35. Pizio O., Patrykiejew A., Soko lowski S., J. Chem. Phys., 2004, 121, 11957; doi:10.1063/1.1818677. 36. Reszko-Zygmunt J., Soko lowski S., Henderson D., Boda D., J. Chem Phys., 2005, 122, 084504; doi:10.1063/1.1850453. 37. Holovko M., Kapko V., Henderson D., Boda D., Chem. Phys. Lett., 2001, 341, 363; doi:10.1016/S0009-2614(01)00505-X. 38. Pizio O., Soko lowski S., J. Chem. Phys., 2005, 122, 144707; doi:10.1063/1.1883165. 39. Yu X.Y., Wu J.Z., J. Chem. Phys., 2002, 117, 10156; doi:10.1063/1.1520530. 40. Rosenfeld Ya., Phys. Rev. Lett., 1989, 63, 980; doi:10.1103/PhysRevLett.63.980. 41. Rosenfeld Ya., J. Chem. Phys., 1993, 98, 8126; doi:10.1063/1.464569. 42. Wertheim M.S., J. Stat. Phys., 1986, 42, 459; doi:10.1007/BF01127721; doi:477.10.1007/BF01127722. 43. Ebeling W., Z. phys. Chem. (Leipzig), 1968, 238, 400. 44. Zhou Q., Yeh S., Stell G., J. Chem. Phys., 1995, 102, 5785; doi:10.1063/1.469310. 45. Borowko M., Pizio O., Soko lowski S. – In: Computational Methods in Surface and Colloid Science, Edited by Borowko M. New York, Marcel Dekker, 2000, Chapter 4. 46. Fischer J., Methfessel M., Phys. Rev. A, 1980, 22, 2836; doi:10.1103/PhysRevA.22.2836. 47. González-Melchor M., Alejandre J., Bresme F., Phys. Rev. Lett., 2003, 90, 135506; doi:10.1103/PhysRevLett.90.135506. 48. González-Melchor M., Alejandre J., Bresme F., J. Chem. Phys., 2005, 122, 104710; doi:10.1063/1.1861878. 13603-11 http://dx.doi.org/10.1063/1.1468638 http://dx.doi.org/10.1088/0953-8984/22/10/104122 http://dx.doi.org/10.1063/1.2192498 http://dx.doi.org/10.1088/0953-8984/8/47/024 http://dx.doi.org/10.1103/PhysRevLett.71.3826 http://dx.doi.org/10.1016/S0167-7322(01)00327-0 http://dx.doi.org/10.1080/00268978000103591 http://dx.doi.org/10.1080/00268978000101581 http://dx.doi.org/doi:10.1103/PhysRevE.57.6944 http://dx.doi.org/10.1063/1.458542 http://dx.doi.org/10.1103/PhysRevA.44.5025 http://dx.doi.org/10.1103/PhysRevE.47.4088 http://dx.doi.org/10.1063/1.480319 http://dx.doi.org/10.1088/0953-8984/14/46/305 http://dx.doi.org/10.1063/1.2128701 http://dx.doi.org/10.1080/00268970802680497 http://dx.doi.org/10.1063/1.2992525 http://dx.doi.org/10.1088/0953-8984/12/12/306 http://dx.doi.org/10.1088/0953-8984/14/46/317 http://dx.doi.org/10.1103/PhysRevE.68.031503 http://dx.doi.org/10.1063/1.1818677 http://dx.doi.org/10.1063/1.1850453 http://dx.doi.org/10.1016/S0009-2614(01)00505-X http://dx.doi.org/10.1063/1.1883165 http://dx.doi.org/10.1063/1.1520530 http://dx.doi.org/10.1103/PhysRevLett.63.980 http://dx.doi.org/10.1063/1.464569 http://dx.doi.org/10.1007/BF01127721 http://dx.doi.org/477.10.1007/BF01127722 http://dx.doi.org/10.1063/1.469310 http://dx.doi.org/10.1103/PhysRevA.22.2836 http://dx.doi.org/10.1103/PhysRevLett.90.135506 http://dx.doi.org/10.1063/1.1861878 A. Patrykiejew, S. Sokołowski, O. Pizio Мiжфазна границя рiдина-пара обмеженої примiтивної моделi iонних плинiв з теорiї функцiоналу густини А. Патрикеєв1, С. Соколовскi1, О. Пiзiо2 1 Унiверситет iм. Марiї Складовської-Кюрi, Люблiн, Республiка Польща, 2 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я, ланцюжки, кристали 13603-12 Introduction Theory Results and discussion

The liquid-vapor interface of the restricted primitive model of ionic fluids from a density functional approach

Репозитарії

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