[

# [

[
Received December 31, 2011, in final form March 6, 2012
###### Abstract

Equation of state for systems with particles self-assembling into aggregates is derived within a mesoscopic theory combining density functional and field-theoretic approaches. We focus on the effect of mesoscopic fluctuations in the disordered phase. The pressure – volume fraction isotherms are calculated explicitly for two forms of the short-range attraction long-range repulsion potential. Mesoscopic fluctuations lead to an increased pressure in each case, except for very small volume fractions. When large clusters are formed, the mechanical instability of the system is present at much higher temperature than found in mean-field approximation. In this case phase separation competes with the formation of periodic phases (colloidal crystals). In the case of small clusters, no mechanical instability associated with separation into dilute and dense phases appears. \keywordsclusters, self-assembly, equation of state, mesoscopic fluctuations \pacs61.20.Gy, 64.10.+h, 64.60.De, 64.75.Yz

###### Abstract

Р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 густу фази, не виникає.

\keywords

кластери, самоскупчення, рiвняння стану, мезоскопiчнi флуктуацiї

\issue

201215223604 \doinumber10.5488/CMP.15.23604 Equation of state in cluster-forming systems]Effect of mesoscopic fluctuations on equation of state in cluster-forming systems A. Ciach, O. Patsahan]A. Ciach\refaddrlabel1, O. Patsahan\refaddrlabel2 \addresses \addrlabel1 Institute of Physical Chemistry, Polish Academy of Sciences, 01-224 Warszawa, Poland \addrlabel2 Institute for Condensed Matter Physics of the National Academy of Sciences of Ukraine,
1 Svientsitskii Str., 79011 Lviv, Ukraine \authorcopyrightA. Ciach, O. Patsahan, 2012

## 1 Introduction

Recent experimental [1, 3, 2, 4], theoretical [7, 9, 10, 5, 11, 6, 8] and simulation [14, 15, 8, 12, 16, 13] studies reveal that in many systems with competing interactions, clusters or aggregates of different sizes and shapes are formed. These objects in certain thermodynamic states can form ordered structures in space [5, 11, 12]. Notable examples include charged globular proteins in water [1, 17, 3, 13], and mixtures of small nonadsorbing polymers with charged colloids or micelles [1, 2, 18]. Interactions (or in fact effective interactions) in the latter systems can be described by the model potential consisting of short-range attraction, resulting from solvophobic or depletion interactions, and long-range repulsion, resulting from screened electrostatic potential (SALR potential).

Systems containing clusters or aggregates are inhomogeneous on the length scale associated with the average size of the aggregates and average distance between them. The corresponding length scale of the inhomogeneities is significantly larger than the size of the particles. Fluctuations on the mesoscopic length scale corresponding to displacements of the aggregates have an important impact on the grand potential, and thus on the equation of state (EOS). Derivation of an accurate EOS for inhomogeneous systems is less trivial than in the case of homogeneous systems, since it is necessary to perform summation over different spatial distributions of the clusters and over all deformations of them.

Contribution to the grand potential associated with mesoscopic fluctuations can be calculated in the field-theoretic approach [5, 11]. In principle, this contribution can be obtained in the perturbation expansion in terms of Feynman diagrams. In practice, an approximate result can be analytically obtained in the self-consistent Hartree approximation [19, 5, 6]. A formal expression for the fluctuation contribution to the grand potential has been derived in references [19, 5, 6]. However, its explicit form with the chemical potential expressed in terms of temperature and density has not been determined yet. The EOS isotherms for various forms of the SALR potential were not analyzed, and the effect of the mesoscopic fluctuations on pressure remains an open question.

It is important to note that various forms of the SALR potential are associated with different properties of the systems. Depending on the ratios between the strengths and ranges of the attractive and repulsive parts of the potential, separation into uniform phases, formation of clusters of various sizes and shapes (globules, cylinders, slabs) in the so-called microsegregation, or isolated individual particles may occur. In certain conditions, the clusters can be periodically distributed in space in the periodic phases whose densities are smaller than the density of the liquid phase [12, 5, 11]. Possible types of the phase diagram for different SALR potentials are shown in figure 1 (see also references [9, 10, 5]).

Properties of the disordered phase can be influenced by the periodic phases for thermodynamic states close to the stability of the latter. We expect that the disordered phase, although the long-range order is absent, is inhomogeneous on the mesoscopic length scale and resembles ‘molten periodic phases’. In this respect the inhomogeneous disordered phase is similar to microemulsion which can be interpreted as molten lyotropic liquid crystal.

In this work we focus on the stable or metastable disordered inhomogeneous phase in which clusters are formed, but they do not form any ordered structure. We derive the EOS with the contribution from mesoscopic length-scale fluctuations included. We calculate the explicit form of the EOS for two representative examples of the SALR potential within the self-consistent Hartree approximation. The first system corresponds to the formation of small clusters, and the gas-liquid separation is unstable for all temperatures in the mean-field (MF) approximation (figure 1, right panel). In the second system, large clusters are formed. The gas-liquid separation is present in this system as a stable or a metastable transition for low temperatures (figure 1, central panel). We shall compare the effect of mesoscopic fluctuations in these two cases on the isotherms , where is the volume fraction of particles and is pressure.

In the next section we briefly summarize the mesoscopic approach. In section 3 the EOS is obtained by two methods. In section 3.1 we consider mesoscopic fluctuations about the average volume fraction, while in section 3.2 fluctuations about the most probable volume fraction are included. The formulas derived in section 3 are evaluated for the two versions of the SALR potential in section 4. We obtain a completely different effect of the mesoscopic fluctuations in these two cases. The two approaches (sections 3.1 and 3.2) yield very close results provided that the relative fluctuation contribution to the average volume fraction is small. For larger fluctuation-induced shifts of the volume fraction, only qualitative agreement of the two methods is obtained. Short summary is presented in section 5.

## 2 Short summary of the mesoscopic description

We consider a local volume fraction of particles, i.e. the microscopic volume fraction averaged over mesoscopic regions, as an order parameter [6]. The corresponding mesoscopic volume fraction varies on a length scale larger than the size of the particles, and the characteristic size of inhomogeneities is the upper limit for the mesoscopic length scale. A particular form of the mesoscopic volume fraction can be considered as a constraint on the microscopic states. The corresponding mesostate is a subset of microstates compatible with the imposed constraint. The mesoscopic volume fraction (or the mesostate) was defined in references [5, 6]. For a one-component case we fix the mesoscopic length scale and consider spheres of radius and centers at that cover the whole volume of the system. We define the mesoscopic volume fraction at by

 ζ(r):=1VS∫r′∈SR(r)^ζ(r′,M), (1)

where , and the microscopic volume fraction in the microstate is defined by

 ^ζ(r,M):=N∑i=1θ(σ2−|r−ri|), (2)

where is the Heaviside unit step function. The microscopic volume fraction is equal to at points that are inside one of the hard spheres, and zero otherwise. Integrated over the system volume, it yields the volume occupied by the particles. The mesoscopic volume fraction at is equal to the fraction of the volume of the sphere that is occupied by the particles. Note that takes the same value if one particle is entirely included in this sphere, independently of the precise position of its centre. Thus, gives less precise information on the distribution of particles than . In the disordered phase is independent of and equals the fraction of the total volume that is occupied by the particles. The mesostate can be imagined as a fixed distribution of centers of clusters, with arbitrary distribution of particles within the clusters, and small modifications of their shapes. Probability of the mesostate is given by [5, 6]

 p[ζ]=\re−βΩco[ζ]Ξ, (3)

where

 Ξ=∫′Dζ\re−βΩco[ζ]. (4)

The functional integral in (4) is over all mesostates,

 Ωco=U[ζ]−TS[ζ]−μN[ζ], (5)

where are the internal energy, entropy and the number of molecules respectively in the system with the constraint of compatibility with the mesostate imposed on the microscopic volume fractions. is given by the well known expression

 U[ζ]=12∫r1∫r2Vco(r1−r2)ζ(r1)ζ(r2), (6)

where for spherically symmetric interactions

 Vco(r1−r2)=V(r12)gco(r1−r2), (7)

, is the interaction potential, is the volume of the particle, and is the microscopic pair correlation function for the microscopic volume fraction in the system with the constraint of compatibility with the mesostate imposed on the microscopic states. The grand potential can be written in the form

 βΩ[¯ζ]=βΩco[¯ζ]−logΞfluc, (8)

where

 Ξfluc=∫Dϕ\re−βHfluc[¯ζ,ϕ], (9)
 Hfluc[¯ζ,ϕ]=Ωco[¯ζ+ϕ]−Ωco[¯ζ]=∑n=1∫r1…∫rnCcon[¯ζ]n!ϕ(r1)…ϕ(rn). (10)

The average mesoscopic volume fraction, , corresponds to the minimum of , and must satisfy the equation

 δβΩco[¯ζ]δ¯ζ(r)+⟨δ(βHfluc)δ¯ζ(r)⟩=0, (11)

where the averaging is over the fields with the probability . Note that when for odd , then the second term on the LHS in (11) vanishes, and the average volume fraction coincides with the most probable volume fraction given by

 δβΩco[ζ]δζ(r)∣∣∣ζ=ζ0=0. (12)

By contrast, when for odd , then .

In order to evaluate the fluctuation contribution to we decompose into two parts

 Hfluc[¯ζ,ϕ]=HG[¯ζ,ϕ]+ΔH[¯ζ,ϕ], (13)

where in the disordered phase

 HG[¯ζ,ϕ]=12∫k~ϕ(k)~C2(k,¯ζ)~ϕ(−k), (14)

and is the Fourier transform of

 C2(r12,ζ)=δ2βΩ[ζ]δζ(r1)δζ(r2). (15)

The above function calculated for is related to the direct correlation function [20].

Assuming , we obtain [19, 5]

 βΩ[¯ζ]≈βΩco[¯ζ]−log∫Dϕ\re−βHG+⟨βΔH⟩G+O(⟨βΔH⟩2G), (16)

where denotes the averaging with the Gaussian Boltzmann factor .

Since , approximate EOS can be obtained from (16) calculated for satisfying (11), when (see (5)) is known. The chemical potential in (16) should be expressed in terms of and ; its form as a function of and can be determined from equation (11). In order to evaluate the second term in equation (11), we need approximate forms of the correlation functions. In the lowest order approximation, it is necessary to determine

 ⟨~ϕ(k)~ϕ(−k)⟩=~G2(k,¯ζ)=1/~C2(k,¯ζ). (17)

## 3 Approximate results for the fluctuation contributions to the EOS, density and chemical potential

In this section we derive an explicit form of the EOS , the density shift, , and the chemical potential under the following assumptions: (i) local density approximation for and (ii) the lowest-order approximation for the second term in (11). In the local density approximation we have

 −TS[ζ]=∫rfh(ζ(r)), (18)

where is the free-energy density of the hard-sphere system with dimensionless density . We assume the Percus-Yevick approximation

 βfh(ζ)=ρ∗ln(ρ∗)−ρ∗+ρ∗[3ζ(2−ζ)2(1−ζ)2−ln(1−ζ)]. (19)

In the local density approximation are just functions of in the disordered phase, and we can simplify the notation, introducing

 An(ζ)=\rdn[βfh(ζ)]\rdζn. (20)

For we have

 Ccon[ζ]=An(ζ), (21)

whereas for

 ~Cco2(k,ζ)=β~Vco(k)+A2(ζ), (22)

where in the disordered phase is the Fourier transform of the function

 Cco2(r12,ζ)=δ2βΩco[ζ]δζ(r1)δζ(r2). (23)

Equations (19) and (6) define the functional for a given form of (equation (7)).

From equation (3) it follows that the most probable fluctuations correspond to the wavenumbers for which assumes the minimum, and the inhomogeneities on the length scale are energetically favored when . In this work we focus on the effect of the self-assembly into aggregates. Therefore, we restrict our attention to which assumes the minimum for , and  [5, 11]. Since the fluctuations with the wavenumber are most probable, they yield the main fluctuation contribution to the grand potential (16). For such fluctuations we can make the approximation

 β~Vco(k)≈β~Vco(kb)+β~V(2)co(kb)(k−kb)2/2+… (24)

As the energy scale we choose the excess energy associated with the fluctuations having unit amplitude and the wavenumber , and introduce the notation

 β∗=1/T∗=β|~Vco(kb)|, (25)
 v∗2=~V(2)co(kb)2|~Vco(kb)| (26)

and

 v∗0=~Vco(0)~Vco(kb). (27)

### 3.1 Fluctuations around the average volume fraction

In this subsection we consider fluctuations about the average value . We shall first determine the chemical potential as a function of and from (11). In order to calculate the second term in (11), we assume that relevant fluctuations are of small amplitudes, and truncate the expansion in (10) at the fourth-order term. Next we insert the derivative with respect to of the RHS of equation (10) truncated at the quadratic term in , and we obtain from (11) and (5) an approximate equation for the rescaled chemical potential, , of the form

 β¯μ≈β¯μMF(¯ζ)+A3(¯ζ)2G(¯ζ), (28)

where

 (29)

and the last term in (28) is the fluctuation contribution with

 G(ζ)=∫k~G2(k,ζ)=G2(0,ζ). (30)

The same expression can be obtained from , when (13), (14) and the approximation

 \re−βHfluc=\re−βHG[1−βΔH+O(ΔH2)] (31)

are used. The approximation (28) is valid as long as the correction term is not larger than the MF result. When considering particular cases we shall verify if this is the case. The fluctuation contribution to the direct correlation function (15) is obtained by calculating the second derivative of the second term on the RHS of (8) with respect to . In the consistent approximation we insert in the obtained expression the appropriate derivatives of (equation (10)) with the expansion in truncated at the second order. The result is given by [19, 21, 5, 6]

 ~C2(k,ζ)≈~Cco2(k,ζ)+A4(ζ)2G(ζ). (32)

Equations (32), (17) and (30) should be solved self-consistently.

The fluctuation induced shift of the volume fraction, , can be obtained from (11) by expanding the first term on the LHS about ,

 (33)

For small we can truncate the expansion in (33) at the first term. When the RHS in (33) is approximated as in the calculation of from (11), we obtain the result

 Δζ≈−A3(ζ0)2~C2(0,ζ0)G(¯ζ)≈−A3(ζ0)2~C2(0,ζ0)+A3(ζ0)G′(ζ0)G(ζ0). (34)

We used the approximations: , and equation (32).

For the potential given in (24) the approximate form of is [22, 19, 21, 5]

 G(ζ)=2a√T∗Z(ζ), (35)

where

 a=k2b4π√v∗2 (36)

and . The above approximation is valid for [22, 19, 21, 5]. The equation (32) for takes the form

 Z(ζ)3=Z(ζ)~Cco2(kb,ζ)+A4(ζ)a√T∗, (37)

and the explicit expression for is

 Z(ζ)=W(ζ)6+2~Cco2(kb,ζ)W(ζ) (38)

with

 (39)

The fluctuation contribution in equation (16) for the approximations (32)–(37) was calculated in references [22, 19, 21, 5], and has the form

 βΩ[¯ζ]≈βΩco[¯ζ]+2a√T∗Z(¯ζ)V−A4(¯ζ)G(¯ζ)28V. (40)

Taking into account (28) and (19), we obtain from (40) the explicit form of the EOS

 βP(¯ζ)=βPMF(¯ζ)+F(¯ζ,T∗), (41)

where

 βPMF(¯ζ)=−β∗v∗02¯ζ2+¯ζ\rdβfh(¯ζ)\rd¯ζ−βfh(¯ζ)=−β∗v∗02¯ζ2+ρ∗¯ζ2+¯ζ+1(1−¯ζ)3 (42)

and

 F(¯ζ,T∗) = a√T∗A3(¯ζ)¯ζZ(¯ζ)−2a√T∗Z(¯ζ)+a2A4(¯ζ)T∗2Z(¯ζ)2 (43) = a√T∗2[2A3(¯ζ)¯ζ−~Cco2(kb,¯ζ)Z(¯ζ)−3Z(¯ζ)].

The second equality in (42) is valid for the PY approximation for . In order to obtain the last equality in (43), equation (37) was used.

### 3.2 Fluctuations around the most probable volume fraction

In the previous subsection we considered fluctuations about the average value, which in general differs from the most probable value of the volume fraction. In principle, it is possible to consider equations analogous to (8) and (10), but with replaced by . In this new approach equation (12) is satisfied, and thus the expansion in (10) starts with (). On the other hand, . The results obtained in the two approaches — with included fluctuations around the average value or around the most probable value — should be the same in the exact theory. However, when the fluctuation contribution is obtained in an approximate theory, the results may depend on the validity of the assumptions made in the two approaches. In this section we derive an alternative version of the EOS, based on the contribution from the fluctuations around the most probable value. From (12) we obtain for the chemical potential

 ¯μ(ζ0)=¯μMF(ζ0) (44)

with given in equation (29). We consider (13) with replaced by , and the approximation (31). In the above, is given in equation (14) with , where , and

 ΔH = 12∫r′∫r′′ϕ(r′)[Cco2(r′−r′′,ζ0)−C2(r′−r′′,ζ0)]ϕ(r′′)+∫rA3(ζ0)3!ϕ(r)3+∫rA4(ζ0)4!ϕ(r)4+…

Note that the function defined here differs from the correlation function, because in this case . Taking into account that for of the form (13) there holds , we obtain the expression

 ⟨ϕ(r)⟩≈−∫Dϕ13!∫r′∫r′′∫r′′′A3(ζ0)ϕ(r′)ϕ(r′′)ϕ(r′′′)ϕ(r)\re−βHGΞfluc. (46)

Finally, the lowest-order result is

 ¯ζ=ζ0+⟨ϕ⟩≈ζ0−A3(ζ0)2~C2(0,ζ0)G(ζ0). (47)

At the same level of approximation is given in equation (32), except that all quantities are calculated at which satisfies (12) rather than (11). This can be verified by a direct calculation of with the help of (31) and (3.2), in an approximation analogous to (46) (see reference [21]).

The above shift of the volume fraction differs from (34), because instead of , there appears . The dependence of on the average volume fraction is given in equations (44) and (47), with eliminated .

In order to evaluate the EOS, we consider an equation analogous to (40), with calculated at its minimum . The EOS takes the form

 βP(ζ0)=βPMF(ζ0)+F0(ζ0,T∗), (48)

where is defined in (42), satisfies (12), and

 F0(ζ0,T∗)=−2a√T∗Z(ζ0)+a2A4(ζ0)T∗2Z(ζ0)2. (49)

The dependence of on the average volume fraction is given by parametric equations (48) with (49) and (47).

The approximate theory developed in this section is valid for small since we assumed that the relevant fluctuations are small and truncated the expansion in (10) at the term . Moreover, to evaluate we neglected the terms of the order . We may expect that if we obtain large and large discrepancies between the results obtained by the two methods, then the approximate theory is not sufficiently accurate.

### 3.3 Comparison between the two methods

Let us focus on the chemical potential, and compare the two expressions, equations (28) and (44) where and satisfy equations (11) and (12), respectively. We expand the RHS in equation (28) about ,

 ¯μ ≈ β~Vco(0)ζ0+βf′h(ζ0)+Δζ[β~Vco(0)+A2(ζ0)+A4(ζ0)2G(ζ0) (50) +A3(ζ0)2G′(ζ0)]+A3(ζ0)2G(ζ0)+O(Δζ2).

From (34), (32) and (44) we obtain an equality of the two expressions for the chemical potential to the linear order in , when is given in (34).

Similarly, to compare the two expressions for the EOS, equations (41) and (48), we expand the RHS of equation (41) about to the linear order in . Taking into account (34) and (37), we arrive at equation (48), up to the terms proportional to . The latter are disregarded in an approximation consistent with the theory for the fluctuation contribution considered in this work. For relatively large , when the terms beyond the linear order become important, discrepancies between the results obtained by the two methods should be expected.

## 4 Explicit results for two model potentials

In this subsection we shall compare the expressions for the chemical potential and for the pressure obtained by the two approaches for two systems showing a qualitatively different behavior. We shall evaluate the EOS (41) for the representative model SALR potential,

 Vco(r)=[−A1r\re−z1r+A2r\re−z2r]θ(r−1), (51)

where is the inverse range in units. The function is a very crude approximation for the pair distribution function. In Fourier representation, the above SALR potential takes the form

 ~Vco(k)=4π[A2\re−z2z22+k2(z2sinkk+cosk)−A1\re−z1z21+k2(z1sinkk+cosk)]. (52)

We choose two sets of parameters, considered in reference [11] in the context of most probable inhomogeneous structures

 System 1:A1=1,A2=0.05,z1=3,z2=0.5; System 2:A1=1,A2=0.2,  z1=1,z2=0.5. (53)

The relevant parameters, , and , (see (26) and (27)) take the values: and

 System 1:kb≈1.79,81v∗0≈−30.145; System 2:kb≈0.6089,v∗0≈0.035. (54)

The two potentials in Fourier representation are shown in figure 3. In the first system small clusters are formed, since is small. Moreover, , and the clusters repel each other. The gas-liquid separation is entirely suppressed due to the very short range of the attractive part of the potential. In the second system large clusters are formed, and (the clusters attract each other). Therefore, in MF, the metastable separation into disordered low- and high density phases occurs at low temperature. Simulation results show that when large clusters are formed, gas-liquid separation occurs for low temperatures, and periodic phases are stable at higher temperatures [9].

We are interested mainly in the part of the phase diagram where the homogeneous structure is less probable than periodic distribution of particles in space, i.e. when does not assume a minimum for (see (3)). We stress that the most probable structure differs from the average structure due to mesoscopic fluctuations. Cluster formation is associated with the excess volume fraction followed by a depleted volume fraction in mesoscopic regions, and the most probable mesoscopic structure associated with cluster formation is periodic. Displacements of the clusters (i.e., mesoscopic fluctuations) can destroy the long-range order, though. Indeed, when temperature is sufficiently high, the average volume fraction takes the constant value as a result of the averaging over cluster displacements, and the disordered inhomogeneous structure with short-range correlations of the cluster positions is found [22, 23]. On the other hand, for low temperatures, the ordered periodic structures are stable [22, 23]. The phase-space region where the inhomogeneous phases (with either short- or long-range order) are stable is enclosed by the structural line [23, 5] given by and shown in figure 3. The structural line is also referred to as -line in literature [26, 9, 10, 24, 25]. Note that in the reduced units (see (25)) the structural line is universal.

We first compare the change of the average volume fraction induced by mesoscopic fluctuations. The shift calculated from (47) and (33) in System 1 is shown in figure 5. The shift is small for a relevant range of temperatures, and both formulas yield practically the same result — they are indistinguishable on the plot. The shift increases for a decreasing temperature. In System 2, the fluctuation contribution to the volume fraction is much larger than in System 1 (figure 5). As expected, when is not very small, , then the two approaches yield somewhat different results, as shown in figure 6 for System 2.

In the next step, we study the chemical potential. The fluctuation contribution in System 1 is small, except at very small volume fractions (figure 7), whereas in System 2 it is substantial, and increases for decreasing temperature, as shown in figures 9 and 9. The two approaches yield similar results for small , whereas when , significant discrepancy between the two approaches is obtained. We can conclude that on the quantitative level the approximate theory is oversimplified for the range of and for which there are significant discrepancies between the two approaches.

Note that since is large for volume fractions , for very small volume fractions our results are oversimplified.

Finally, we present the isotherms obtained from (41) and (48) for the two systems in figures 1114. In System 1 the pressure is much higher than found in MF, and for all temperatures it monotonously increases with , as shown in figure 11. The increased pressure associated with mesoscopic fluctuations may result from the repulsion between the clusters, because in this case .

Since in System 2 , a mechanical instability develops at the MF spinodal line, with the metastable MF critical point . What is really interesting is that such instability appears at much higher temperatures due to mesoscopic fluctuations. This is in strong contrast to the fluids with purely attractive interactions, where density fluctuations decrease the critical temperature with respect to the mean-field estimate. The present case with dominant fluctuations associated with mesoscopic wavelengths bears some resemblance to the restricted primitive model (RPM) of ionic systems. There is no gas-liquid instability in the RPM at the MF level of a mesoscopic theory analogous to the one considered here, but when the short-wavelength charge-density fluctuations are included, such instability appears [24]. One could imagine that the mesoscopic fluctuations, i.e., displacements of the clusters from their most probable locations lead to their coalescence when , and thus support the phase separation. Note that properties of System 2 are completely different from the previously studied System 1.

In this region of the phase diagram, the shift of the volume fraction is large. Therefore, on the quantitative level, the results are not sufficiently accurate. The inflection point on the isotherm appears at or according to (48) or (41) with (47), respectively. Both temperatures, however, are much higher than in MF. Further studies are required to verify if the separation into disordered inhomogeneous phases, or periodic ordering of clusters occurs. If the ordered phases are formed, the still open question is for which part of the phase diagram such phases are globally stable.

Finally, let us focus on the pressure for very small volume fractions. In the fluctuation correction to pressure (equation (43)) the first term comes from the fluctuation contribution to chemical potential (see (28)). As shown in figure 7 (right panel), for very low volume fractions, our approximation is oversimplified, so the negative pressure is an artifact. For very low volume fractions, we should expect a perfect gas behavior, except that some fraction of particles should form clusters. Pressure should be proportional to the sum of the number densities of monomers and clusters. Since the number of clusters is smaller than the number of particles forming them, pressure should be smaller than in the corresponding perfect gas of isolated particles. Our theory agrees with this expectation (see figure 12, right panel).

## 5 Summary

In this work, the effects of mesoscopic fluctuations on the average volume fraction, chemical potential and pressure as functions of temperature and the average volume fraction were considered within the framework of the mesoscopic theory [5, 6]. We restricted our attention to a stable or metastable disordered phase. The fluctuation contribution to the quantities mentioned above was calculated in two ways. First, we considered the fluctuations about the average volume fraction, and derived equations (33), (28) and (41) for the volume fraction, chemical potential and pressure, respectively. In the second version, we considered the fluctuations about the most probable volume fraction, and obtained equations (47), (44) and (48), with satisfying (12). The chemical potential and the EOS as functions of the average volume fraction are given by parametric equations (47) and (44), and (48), respectively. Our expressions are derived under the assumption that the dominant fluctuations are of small amplitudes. Consistent with the above assumption, the two methods yield the same result to a linear order in the fluctuation contribution to the volume fraction .

The fluctuation contributions to all three quantities were explicitly calculated for two versions of the SALR potential. In System 1 the zeroth moment of the effective interactions is positive and small clusters are formed. In System 2 the zeroth moment of the effective interactions is negative, and the clusters are large. We obtain nearly the same results independently of the method used when the fluctuation-induced shift of the volume fraction is very small (System 1). When (low temperature in System 2), significant discrepancies between the two methods appear for some part of the phase diagram. The largest discrepancies are present when both methods yield the results that strongly deviate from the MF predictions. The larger is the probability of finding inhomogeneous mesoscopic states compared to the homogeneous distribution of particles, the stronger are the discrepancies between the two methods. It is in this part of the phase diagram that the periodic order may appear. We conclude that the first method is superior to the second one because it is easier to implement. Exact results would be necessary to get a comparison between the accuracy achieved by these methods.

We have found that mesoscopic fluctuations play a very important role and lead to a significant change of the chemical potential and pressure. The larger is the probability of finding the inhomogeneities, i.e., the further away from the structural line on the low- side of it (figure 3), the larger is the role of fluctuations. When small clusters are formed (System 1 in section 4), the fluctuation contribution to pressure increases monotonously with an increasing volume fraction. By contrast, for large clusters (System 2 in section 4) the fluctuation contribution to pressure is nonmonotonous; it is negligible for small as well as for large volume fractions, whereas for intermediate volume fractions it is large and increases with a decreasing temperature (figures 1114). Moreover, an inflection point at the pressure – volume fraction isotherm appears at the temperature and volume fraction both much larger than found in MF. Further studies are required to verify if the periodically ordered cluster phases are stable in System 2, or phase separation occurs due to the mechanical instability. Possible scenarios are: (i) phase separation at low , and periodic structures at higher , or (ii) the phase separation is only metastable, and finally (iii) the periodically ordered phases are only metastable. For System 1, the phase separation is not expected.

In addition to the assumptions discussed earlier, we make an approximation concerning the form of the effective potential . Note that in equation (51) we assumed that the pair distribution function vanishes for in units. For the volume fraction, this is a poor approximation, and quantitative results for the structural line depend on the form of the pair distribution function (regularization of the potential [27]).

In the future studies, the EOS for periodically ordered cluster phases should be determined in order to find the phase diagram. Our results indicate that despite the universal properties of the dependence of the most probable structures on the thermodynamic state, the effect of fluctuations on the average distribution of particles may depend on the form of the interaction potential, especially on the sign of the zeroth moment of the effective interactions.

## Acknowledgements

A part of this work was realized within the International PhD Projects Programme of the Foundation for Polish Science, co-financed from European Regional Development Fund within Innovative Economy Operational Programme ‘‘Grants for innovation’’. Partial support by the Ukrainian-Polish joint research project under the Agreement on Scientific Collaboration between the Polish Academy of Sciences and the National Academy of Sciences of Ukraine for years 2009–2011 is also gratefully acknowledged.

## References

• [1] Stradner A., Sedgwick H., Cardinaux F., Poon W.C.K., Egelhaaf S.U., Schurtenberger P., Nature, 2004, 432, 492; \doi10.1038/nature03109.
• [2] Campbell A.I., Anderson V.J., van Duijneveldt J.S., Bartlett P., Phys. Rev. Lett., 2005, 94, 208301;
\doi10.1103/PhysRevLett.94.208301.
• [3] Porcar L., Falus P., Chen W.-R., Faraone A., Fratini E., Hong K., Baglioni P., Liu Y., J. Phys. Chem. Lett., 2010, 1, 126; \doi10.1021/jz900127c.
• [4] Stiakakis E., Petekidis G., Vlassopoulos D., Likos C.N., Iatrou H., Hadjichristidis N., Roovers J., Europhys. Lett., 2005, 72, 664; \doi10.1209/epl/i2005-10283-y.
• [5] Ciach A., Phys. Rev. E, 2008, 78, 061505; \doi10.1103/PhysRevE.78.061505.
• [6] Ciach A., Mol. Phys., 2011, 109, 1101; \doi10.1080/00268976.2011.638329.
• [7] Sear R.P., Gelbart W.M., J. Chem. Phys., 1999, 110, 4582; \doi10.1063/1.478338.
• [8] Archer A.J., Wilding N.B., Phys. Rev. E, 2007, 76, 031501; \doi10.1103/PhysRevE.76.031501.
• [9] Archer A.J., Pini D., Evans R., Reatto L., J. Chem. Phys., 2007, 126, 014104; \doi10.1063/1.2405355.
• [10] Archer A.J., Phys. Rev. E, 2008, 78, 031402; \doi10.1103/PhysRevE.78.031402.
• [11] Ciach A., Góźdź W.T., Condens. Matter Phys., 2010, 13, 23603; \doi10.5488/CMP.13.23603.
• [12] De Candia A., Del Gado E., Fierro A., Sator N., Tarzia M., Coniglio A., Phys. Rev. E, 2006, 74, 010403(R); \doi10.1103/PhysRevE.74.010403.
• [13] Kowalczyk P., Ciach A., Gauden P.A., Terzyk A.P., J. Colloid Interface Sci., 2011, 363, 579;
\doi10.1016/j.jcis.2011.07.043.
• [14] Sciortino F., Mossa S., Zaccarelli E., Tartaglia P., Phys. Rev. Lett., 2004, 93, 055701;
\doi10.1103/PhysRevLett.93.055701.
• [15] Sciortino F., Tartaglia P., Zaccarelli E., J. Phys. Chem. B, 2005. 109, 21942; \doi10.1021/jp052683g.
• [16] Toledano J., Sciortino F., Zaccarelli E., Soft Matter, 2009, 5, 2390; \doi10.1039/b818169a.
• [17] Shukla A., PNAS, 2008, 105, 5075; \doi10.1073/pnas.0711928105.
• [18] Zhang T.H., Groenewold J., Kegel W.K., Phys. Chem. Chem. Phys., 2009, 11, 10827; \doi10.1039/b917254h.
• [19] Ciach A., Góźdź W.T., Stell G., J. Phys.: Condens. Matter, 2006, 18, 1629; \doi10.1088/0953-8984/18/5/016.
• [20] Evans R., Adv. Phys., 1979, 28, 143l \doi10.1080/00018737900101365.
• [21] Patsahan O., Ciach A., J. Phys.: Condens. Matter, 2007, 19, 236203; \doi10.1088/0953-8984/19/23/236203.
• [22] Brazovskii S.A., Sov. Phys. JETP, 1975, 41, 85.
• [23] Ciach A., Patsahan O., Phys. Rev. E, 2006, 74, 021508; \doi10.1103/PhysRevE.74.021508.
• [24] Ciach A., Stell G., J. Mol. Liq., 2000, 87, 255; \doi10.1016/S0167-7322(00)00125-2.
• [25] Ciach A, Góźdź W.T., Evans R, J. Chem. Phys., 2003, 118, 3702; \doi10.1063/1.1539046.
• [26] Stell G., In: New Approaches to Problems in Liquid-State Theory, edited by C. Caccamo, J.-P. Hansen, and G. Stell, Kluwer Academic Publishers, Dordrecht, 1999.
• [27] Patsahan O.V., Mryglod I.M., Condens. Matter Phys., 2004, 7, 755.
\ukrainianpart

Вплив мезоскопiчних флуктуацiй на рiвняння стану кластероутворювальних систем А. Цях\refaddrlabel1, О. Пацаган\refaddrlabel2 \addresses \addrlabel1 Iнститут фiзичної хiмiї, Польська академiя наук, 01-224 Варшава, Польща \addrlabel2 Iнститут фiзики конденсованих систем Нацiональної академiї наук України,
вул. Свєнцiцького, 1, 79011 Львiв, Україна

\makeukrtitle
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters