Explicit description of complexation between oppositely charged polyelectrolytes as an advantage of the random phase approximation over the scaling approach.

A polyelectrolyte complex (PEC) of oppositely charged linear chains is considered within the Random Phase Approximation (RPA). We study the salt-free case and use the continuous model assuming a homogeneous distribution of the charges throughout the polyions. The RPA correction to the PEC free energy is renormalized via subtraction of polyion self-energy in order to find the correlation free energy of the complex. An analogous procedure is usually carried out in the case of the Debye-Hückel (DH) plasma (a gas of point-like ions), where the infinite self-energy of point-like charges is subtracted from the diverging RPA correction. The only distinction is that in the PEC both the RPA correction and chain self-energy of connected like charges are convergent. This renormalization allows us to demonstrate that the correlation free energy of the PEC is negative, as could be expected, while the scaling approach postulates rather than proving the negative sign of the energy of interactions between the blobs. We also demonstrate that the increasing concentration of oppositely charged polyions in the solution first results in the formation of neutral globules of the PEC consisting of two polyions as soon as the concentration reaches a certain threshold value, cgl, whereas solution macroscopic phase separation (precipitation of globules) occurs at a much higher concentration, ccoac, ccoac ≫ cgl. Partitioning of polyions between different states is calculated and analytical dependencies of cgl and ccoac on the polyion length, degree of ionization and solvent polarity are found.


Introduction
Correlation effects are of major importance in many charged systems, especially if they comprise macro-and/or multivalent co-and counterions, 1 which control the strength and character of electrostatic interactions. In addition to the valency of the ions, phase behavior of charged liquids (solutions of colloidal particles, polyelectrolytes (PEs), etc.) depends on relative ''asymmetry'' of co-and counterions. For example, bulky macroions (like charged colloidal particles, microgels 2,3 or linear PEs 4 ) with small (pointlike) monovalent counterions can be attributed to strongly ''asymmetric'' charged solutions where high entropic mobility of the counterions protects complete local neutralization of the system and formation of a condensed phase. Therefore, the presence of charged groups in such asymmetric systems can promote solubility of the macroions. In contrast, in charge ''symmetric'' systems, like solutions of oppositely charged PEs, [5][6][7][8][9] condensation of the macroions occurs in the isoelectric point (total charge of the macrions is equal to zero). Therefore, the electrostatic attraction decreases the solubility of the macroions in the ''symmetric'' system. [10][11][12] In the intermediate cases of the asymmetry or beyond the isoelectric point, the interplay between attractive and repulsive electrostatic interactions is responsible for the clustering of the macroions in soluble aggregates. In many cases, solubility of the aggregates is provided by the effect of overcharging (charge inversion): [13][14][15][16] the total charge of the complex becomes opposite to the bare charge of the macroions.
The spatial charge distribution in the ''asymmetric'' systems is essentially inhomogeneous. It allows developing mean-field theories (Poison-Boltzmann or multi-zone approximation) like for dilute solutions of PE micelles [17][18][19][20] or microgels. 21 On the other hand, the inhomogeneous charge distribution means strong spatial charge fluctuations which are equivalent to the strong charge fluctuations at a certain spatial coordinate evolving over time as long as the ergodic hypothesis is valid. Therefore, calculations of the fluctuation contribution to the free energy in such systems are sophisticated. In contrast, charge fluctuations are much smaller in the ''symmetric'' solutions and the mean-field treatment is not able to describe the Coulomb interactions properly (they are equal to zero at this level). The electrostatics due to the small charge fluctuations can be taken into account within the framework of the Random Phase Approximation (RPA), which is a linear response approximation analogous to the one-loop approximation used in quantum field theory.
Over 90 years ago, Debye and Hückel (DH) have considered the problem of electrolyte solutions (the simplest ''symmetric'' system) demonstrating deviations from the ideal gas behaviour. 22 Using the Poison-Boltzmann equation linearization, they have calculated corrections to the osmotic pressure and to the free energy of the ideal gas, both are negative. However, independently of the valency of the ions, their condensation or ion pair formation cannot be described within such an approach. The condition for linearization is that the energy of the electrostatic interactions between the ions is much smaller than the energy of their thermal motion. Three decades later, the first theoretical attempt to describe saline solutions of oppositely charged PEs was undertaken by Overbeek and Voorn. 23 They have treated non-electrostatic interactions within the Flory-Huggins lattice theory, while Coulomb interactions between all charges of the system (mobile lowmolecular-weight ions and charged groups of the chains) have been described using the DH expression. Thus, connectivity of the charged groups into the chains was neglected.
In the late 1980s in pioneer works of Borue and Erukhimovich, the RPA was applied to semidilute PE solutions and complexes, where the charge connectivity into the chain was taken into consideration by the addition of the polyion conformational entropy term to the free energy functional. 24,25 The beauty of the RPA for the oppositely charged PEs is that the condensation of the chains (complex formation) can properly be described at this level because polymer chains are known to be ''poor in entropy'': even small attraction between charged groups is able to overcome the entropic penalty due to the complex formation. Possibility of microphase separation in semidilute PE solutions and mixtures of oppositely charged polyions was predicted within the framework of this approach. 24,[26][27][28] Later on, similar theory was applied for the description of PE multilayers, 29 polyampholytes 30 and overcharged interpolyelectrolyte complexes. 31,32 Further development of the RPA included consideration of ion pairing, 28,33,34 short wavelength integration cutoff, 33,35 hard-core and finite-size ion effects. 33,34,36,37 In our group, fluctuations of polymer and charge densities were examined independently (discrete charge distribution), which allowed getting a general expression for the free energy unifying different limiting cases. 38 The RPA theory for rod-like polymer chains capable of liquid-crystalline ordering was developed in ref. [39][40][41][42]. Treatment of single PE chain trapping own counterions as one-component plasma (i.e. point-like charges immersed into homogeneously and oppositely charged background) is also based on the RPA. 35,[43][44][45][46] Castelnovo and Joanny, 47 and then Fredrickson and co-workers 48 extended the fluctuation theory of complexation beyond the RPA.
Recently account for fluctuations by means of liquid state theory approaches allowed describing solutions 49 and complexes 50 of highly charged PEs. In the framework of this method it was theoretically shown that low-molecular-weight salt ions prefer a supernatant rather than a coacervate, 50 in contrast to predictions of Vroon-Overbeek theory that neglects both connectivity of PE charges and excluded volume interactions. 37,51,52 This finding on salt ion partitioning being in principal agreement with experiment was later confirmed using computer simulations. 53,54 In the present paper, we propose a procedure of physically motivated renormalization of the RPA results for the symmetric PEC allowing us to distinguish the correlation (free) energy of the PEC and self-energy of polyions in the salt-free case. Proper renormalization makes possible to find partitioning of polyions between different states (non-aggregated chains, neutral PEC globules and coacervates) in a stoichiometric mixture of polyanions and polycations. Threshold concentrations of globule formation and macroscopic phase separation are found analytically.
The model with a charge continuously smeared throughout the flexible chain is used for the theoretical treatment of the PEC. 24,55 The RPA correction to the free energy is known to be convergent within the framework of this model. [29][30][31][32] However, the knowledge of this positive expression, which allows us to find the electrostatic contribution to the osmotic pressure and the polymer volume fraction in the complex properly, is not enough for calculation of the correlation (free) energy. Nonzero self-energy of the chains caused by connectivity of the charged groups has to be calculated and subtracted from the RPA results. The self-energy takes into account repulsive interactions of the charges in each individual chain. Due to the subtraction procedure, the correlation energy, which corresponds to the attraction between oppositely charged chains, occurs to be negative as one would expect. In contrast to the scaling approach, where the interaction energy between the blobs was postulated to be negative, 56 the RPA proposes a straightforward way for calculation of the correlation free energy defining the numerical factor as soon as one deals with the polymers in a Y-solvent. In particular, the knowledge of the self-energy allows quantifying the association of single chains into a complex (calculation of the critical concentration of association).
Here it is necessary to mention that initially a renormalization procedure of the RPA results for PE solutions based on the subtraction of chains self-energy was proposed by Mahdi and Olvera de la Cruz who motivated it primarily by mathematical convenience. 36 Lately Shen and Wang have explicitly distinguished self-energy of polyions in PE solution and found solution correlation energy via subtraction of the chain selfenergy from the formal field-theoretic result. 57 In the present paper, we apply the same renormalization procedure to the case of the PEC. It makes possible to argue for the first time the negative sign of the complex correlation energy within the RPA.
In order to make the calculation of the correlation free energy in the case of PECs more lucid and clear, we use analogy with wellknown DH plasma wherein subtraction of infinite self-energy of point-like ions from the diverging RPA result is necessary to get a finite negative expression for the correlation free energy. An analogy between correlation functions in the case of the PE complex and the DH plasma also occurs to be very useful.
The paper is organized as follows. In the next section the DH plasma is discussed in detail, first without and then with volume interactions, and well-known results for ion-ion correlation functions and correlation free energy of plasma are reproduced. In Section 3, we move on to the case of the PEC and propose a renormalization procedure of the RPA free energy correction. It allows us to obtain the correlation free energy of the complex. Section 4 deals with the calculation of the critical concentrations of PEC formation and coacervation. Closed-form expressions for these concentrations are found.

Debye-Hückel plasma
We consider a neutral plasma consisting of oppositely charged ions, each of the size a and volume a 3 . Let the average volume fractions of positively and negatively charged ions are equal: hn + i = hn À i = n 0 . Here n 0 is dimensionless, so that the concentration of each type of the ions equals n 0 /a 3 . Neglecting short-range interactions between the ions under the condition of a small concentration, n 0 { 1, one gets the classical case of the DH plasma of point-like charges. The total free energy is a sum of two terms: the ideal gas entropy and the energy of long-range Coulomb interactions between the ions ð n þ ðrÞ À n À ðrÞ ð Þn þ ðr 0 Þ À n À ðr 0 Þ ð Þ jr À r 0 j d 3 rd 3 r 0 : Here, l b = e 2 /(ek B T) is the Bjerrum length. Correlated electrostatic attraction makes the osmotic pressure of the DH plasma lower than that of the ideal gas. The RPA can be used in order to find the free energy correction to the ideal gas term and to calculate the density-density correlation functions of the DH plasma. The concentration fluctuations are defined as dn k (r) = n k (r) À n 0 , and Fourier transforms are given by where k = +, À. Deviation of the free energy functional from the mean-field value, F MF /k B T = 2Vn 0 ln(n 0 /e)/a 3 , in the Gaussian (square) approximation reads: where V is the total system volume. Owing to the following representation of the free energy fluctuations and introducing the Debye radius r D 2 = a 3 /(8pl b n 0 ), one can obtain the direct matrix of the plasma correlation functions in the Fourier space: The spherically symmetric density-density correlation functions G kl (r) = hdn k (0)Ádn l (r)i are calculated using the inverse Fourier transforms and take the form The second set of the correlation functions, G +À = G À+ , is caused only by the electrostatic attraction between the oppositely charged ions: it is equal to zero at l b = 0. For this reason it can be designated as G +À = G corr . The first set of correlation functions, G ++ = G ÀÀ , should be divided into two parts: G ++ = G self À G corr . The first nonvanishing at the l b = 0 part, G self = n 0 a 3 d(r) is known as a correlation function of the neutral ideal gas while the second part, ÀG corr , should be ascribed to the repulsion of like charges. Note that both G corr and G self are positive at any values of r.
Owing to the expression hdn k (r)dn l (r 0 )i = G kl (|r À r 0 |) and denoting R = |R| = |r À r 0 |, one can calculate Thus, the self-energy of point-like ions and the energy of interaction between the ions are given by It is well-known that the self-energy of point-like charge diverges: where the delta function in the spherical coordinates was substituted, d(R) = d(R)/(4pR 2 ). On the other hand, the correlation internal energy is finite: The energy is negative and decreases with the increase of the concentration n 0 . It means that the electrostatic interactions are attractive. To pass to the Helmholtz free energy correction, one can apply the relation between thermodynamic potentials As a result, the free energy density correction reads and constant A 1 = 0 because the condition F corr = (F tot À F id ) -0 at T -N should be satisfied, F id being the free energy of neutral ideal gas equal to F MF . 58 Another way to get the correlation free energy correction is based on the calculation of Gaussian functional integrals In order to calculate the integrals, quadratic form should be diagonalized. After transformations, one gets the following expression diverging at q -N, i.e. at r -0. This divergence is a consequence of infinite self-energy of point-like charges. F RPA includes both the self-and the interaction energies, The self-energy is given by eqn (15) and can be rewritten using Fourier transform: Finally, the converging result for the correlation free energy coincides with one deduced above. It is also worth noting that considering ions of non-zero radius and finite charge density function instead of point-like ions with d-function type charge distribution, one can avoid all energy divergences. 59 As we mentioned in Introduction, the above calculations of the correlation (free) energy are valid when this correction is small as compared to the energy of thermal motion, |F corr | { |F tr |. This condition is fulfilled at a low concentration of the ions, n 0 { (a/l b ) 3 .

Debye-Hückel plasma with volume interactions
In order to take into account volume interactions of the plasma ions, we can use the virial expansion and the following term should be added to the total free energy Here B and C are the second and the third dimensionless virial coefficients. After expansion of the total free energy functional into the series in powers of the density fluctuations, the Fourier transforms of the correlation functions take the form: were m = 4n 0 (B + 6Cn 0 ). The value m = 0 corresponds to the case without volume interactions, and the matrix coincides with the matrix 8. The correlation functions read Due to the RPA linearity, the correlation functions can be represented as a sum of (i) the ideal gas self-correlation function, G self = n 0 a 3 d(r), (ii) the term caused by the volume interactions, G vol = Àmn 0 a 3 d(r)/(2 + 2m), and (iii) the term G corr = l b n 0 2 e Àr/rD /r induced by the electrostatic interactions: The terms G self and G corr coincide with those defined above for the DH plasma without volume interactions. G vol term vanishes when the triple repulsion compensates pairwise attraction, B = À6Cn 0 . Following eqn (19), one can calculate the RPA correction to the free energy caused by the Coulomb interaction. Despite the presence of the volume interactions, the result for F RPA coincides with eqn (21). Thus, the self-energy given by eqn (13) should be subtracted from F RPA to find the correlation energy, and the latter coincides with correlation energy for the DH plasma without volume interactions, F corr /k B TV = À1/(12pr D   3 ). Independence of F RPA and F corr of the volume interactions should be attributed to the linearity of the RPA.
Note that eqn (11)- (14) remain valid in the case of non-zero volume interactions as well, so that definitions of the self-energy and the correlation energy are the same. Indeed, relationship G ++ À G +À = G self À 2G corr is satisfied despite G vol a 0. In the general case, the term caused purely by the electrostatic interactions should be defined as In the next section we apply the proposed approach to describe complex formation between oppositely charged PE chains.

Interpolyelectrolyte complex
Let us consider a complex of two oppositely charged linear PEs in a salt-free, Y-solvent. For the sake of simplicity, we study the case of a symmetric complex, i.e. f + = f À = f are the fractions of charged groups in these polyions, and F + and F À are their average polymer volume fractions in complex, respectively. The complex electric neutrality is provided by condition F + = F À = F tot /2, where F tot is the total polymer volume fraction within the complex. The lengths of polycations and polyanions are also assumed to be equal, N + = N À = N (see ref. 60 for the asymmetric case). We study the classical model proposed by Borue and Erukhimovich which assumes that electric charge is continuously smeared throughout the chains. 24,25 The total free energy in k B T units takes the form: The first term is responsible for the translational entropy of polymer chains: The second term is the conformational entropy of polyions 61,62 with a being the length of the statistical segment. Deduced independently by Edwards and later by Lifshitz, it takes into account the entropy cost of the polymer non-uniform distribution. The case of flexible polymer chains is under consideration, so that a 3 is of the order of the segment volume. The Coulomb energy is given by We restrict our consideration to the case of weakly charged PEs, f { 1, and a rather polar solvent, u B 1, so that effects of crosschain ion pairing [63][64][65] and attraction of resultant dipoles [66][67][68] can be neglected. Finally, the volume interactions, namely triple repulsive interactions between the segments in the Y-solvent, are taken into account within the last term of the total free energy with C being the dimensionless third virial coefficient. The mean-field value of the total free energy of the neutral complex takes the form: where the second term is negligible at high N, i.e. the volume interactions dominate over the chain translational entropy. Thus, the repulsive interactions make the complex unstable. The complex formation cannot be described at the mean-field level because the latter neglects the electrostatic attractions between charged chains. At the same time, the complex formed by oppositely charged PEs is globular, and weak fluctuation within the globule provides applicability of the RPA (see ref. 25 for detailed discussion). Expansion of the free energy functional into the series in powers of the polymer density fluctuations allows calculating the inverse matrix of correlation functions being the correlation radius of the ideal polymer chain (i.e. the Gaussian chain size), l vol ¼ a being the correlation radius of infinite neutral chain caused solely by the volume interactions and x = a/(48puf 2 F tot ) 1/4 being the screening radius of the electrostatic interactions in PE complex; u = l b /a is the dimensionless Bjerrum length. The RPA correction to the free energy of the PEC can be found in a similar way to that for the DH plasma (see eqn (19) and (21)): and this integral converges at any r c and any finite x (infinite x corresponds to the case of neutral polymers, f = 0). Thus, one gets F RPA is always positive as soon as the complex is formed (x o r c ). In order to get the correlation free energy of the PEC, the selfenergy of polyions should be subtracted from the RPA results.
Following an analogy between the PEC and the DH plasma, the correlation functions should be divided into the parts related to connectivity of monomer units into the chain, G self , electrostatic, G corr , and volume, G vol , interactions: The term G self appears only in the correlation functions of like charges as a consequence of their connectivity into a chain. Indeed, the charges of opposite signs belong to different chains and are not linked with each other by the covalent bonds. In order to find the correlation free energy and the selfenergy, one should explicitly calculate these components of the correlation functions: G vol ðrÞ ¼ 3aF tot 4pr e Àr=lcorr À e Àr=rc Â Ã (43) Here l corr is known to be the correlation radius of neutral polymer solution with the radii x AE ¼ x . ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi 1 Ç x 2 =2r c 2 p and dimensionless d = arctan(x + x À /2r c 2 ). Note that the electrostatic interactions vanish at l b = 0, and G corr = 0. Similarly, G vol = 0 in the absence of volume interactions, i.e. at C = 0. G self is known to be the Ornstein-Zernike type correlation function of a single Gaussian chain and is equal to the correlation function in the solution of the ideal chains. 69 Using relationship similar to eqn (13) and introducing additional multiplier f 2 for the transition from the polymerpolymer to charge-charge correlation function, one obtains where N chains = VF tot /(a 3 N) is the total number of polymer chains (both polycations and polyanions) in the system. Since u B 1/T, the free and the internal electrostatic energies coincide, F self = U self , and the density of the self-energy can be represented as follows: It is necessary to emphasize that the self-energy term is not a mean-field energy of the charged chains with Gaussian statistics but the correction term caused by Gaussian correlations in the location of charges (see Appendix A for a detailed discussion). The value of the self-energy term is defined by the choice of the reference polymer system: it is the system without volume and electrostatic interactions, C = 0 and l b = 0. Therefore, the chains in the reference system possess Gaussian statistics at any length scales owing to the conformational entropy term F conf in the total free energy.
Thus, similar to the case of the DH plasma (eqn (23)), the correlation energy can be written as follows: The correlation energy (48) is always negative owing to the negative integrand, and this result fully coincides with common comprehension of correlation energy as negative energy taking into account the correlation attraction of opposite charges. In the case of an infinitely long chain, N -N, the second term in square brackets diverges at q -0, i.e. at large distances. Indeed, the Coulomb correlation energy of an infinitely long PE chain with the Gaussian statistics at any length scales is infinite and grows with increasing chain length proportional to ffiffiffiffi N p , see eqn (47). For this reason the term accounting for the chain translational entropy cannot be neglected in the total free energy.
Result similar to eqn (48) has been used earlier by Mahdi and Olvera de la Cruz who calculated the phase diagram of semidilute PE solution in the presence of salt. 36 The authors justified subtraction of the second term in square brackets by mathematical convenience: this ''. . . irrelevant subtracted term only facilitates the algebra and gives an electrostatic contribution that reduces to the simple Debye-Hückel electrostatic free energy in the correct limit.'' 36 In fact, subtraction of the selfenergy does not influence results for PE solution since it is a constant term and can be omitted: the self free energy density is proportional to the polymer volume fraction F and does not contribute to the osmotic pressure. 70 However, if we would like to calculate the correlation energy of the PEC or to compare two states of the PE chains -aggregated into PEC and free -it is necessary to properly calculate and comprehend calculation of the free energy of the system. Subtraction of the self-energy of the chains is very important for distinguishing the reference state of the system and, hence, proper calculation of energy of Coulomb interactions. Application of the renormalization procedure to the case of the PEC and the proof of the negative sign of the correlation free energy are the novelties of the above part of the present work.
The self-energy is proportional to the number of PE subchains in the system N chains , like the self-energy is proportional to the number of charges in the DH plasma. The difference between these two systems is the following: point-like charges are structureless and their self-energy is unchanged under any conditions, while the PE chain can adopt different conformations resulting in different mutual spatial locations of charges and different self-energy.
The same identification of the PE chain self-energy as the energy necessary to connect charges onto a single chain has been recently used by Shen and Wang, and our choice of the reference state coincides with theirs: ''the zero energy of the electrostatics is taken to be the state where charges are dispersed into infinitesimal bits at infinity.'' 57 The dependence of the chain self-energy on its conformation was discussed in ref. 57 as well. The authors introduced chain self-energy in order to find correlation free energy of PE solution, and we use very similar renormalization with the same chain self-energy definition. The only difference is that we apply it to the case of the interpolyelectrolyte complexes rather than solution of similarly charged chains.
The positive value of the PRA correction (39) means that the number of like charges surrounding any selected charge on the chain in the PEC is more than the number of opposite charges, despite repulsion of like charges and attraction of opposite charges. This fact is caused by the binding of like charges into the chain. Since these bonds exist before complex formation, their impact on the charge-charge correlations was separated from the correlations induced by purely Coulomb attraction. The Coulomb energies caused by these correlations were treated as the self-energy and the correlation energy, respectively.

Equilibrium density of the PEC
In order to find the equilibrium polymer volume fraction within the complex, the sum of the mean-field and the correlation free energies should be minimized. The minimization is equivalent to the equating of the osmotic pressure of the complex to zero: Here we neglect the translational entropy of polyions within the complex and the nonzero concentration of free polyions in the outer solution since their concentration is exponentially small (see Section 4). The RPA correction F RPA can be used instead of F corr in order to find the PEC density because the self-energy does not contribute to the osmotic pressure. At that, F RPA represents the total electrostatic energy of the system, i.e. the energy with respect to the state when all infinitesimal (not elementary) charges are moved away from each other at infinite distances. This choice of the reference state assumes that polymer chains are broken, and the energy required to gather them back into a single chain and provide the chain with the Gaussian conformation U G,ch el-st equals mean-field electrostatic chain energy U G,ch MF plus the self-energy of a chain U G,ch corr , see Appendix A. In the case of infinitely long chains (N -N), we get , so that the polymer volume fraction within the complex reads (50) and the free energy of the PEC per unit volume equals Thus, the free energy per chain within the macroscopic PEC phase takes the form: Scaling analysis of the chain free energy within the PEC performed in ref. 56 predicts F coac E ÀNu 2/3 f 4/3 . The sign minus in the scaling result was justified by attraction of oppositely charged blobs within the complex. 56 However, the negative energy of correlation attraction between oppositely charged blobs is of the order of the blob self-energy. So the scaling consideration does not allow us to define the sign (plus or minus) of the free energy within the complex, F coac . It is necessary to mention here that properly defined correlation energy of the PEC given by eqn (48) is indeed negative, but the free energy of the complex being the sum of the correlation free energy and the self-energy of the chains is positive. Moreover, scaling treats polyions as a set of disrupted blobs resulting in a linear dependence of the chain correlation energy on its length N, while RPA accounts for correlations within the chain on all length scales and yields asymptotically N 3/2 .

Equilibrium between nonaggregated charged chains and their complex
In order to calculate the critical concentration of PEC formation properly, one should take into account correlation corrections to the mean-field energy of both neutral stoichiometric complex and single charged chains. Below we propose a way of calculation of the self-energy of free PE chains in the solution.

Self-energy of a single polyelectrolyte chain
It is well-known that the PE chain in dilute solution adopts elongated (swollen) rather than Gaussian coil conformation,  monomer units in the blob g B (x blob /a) 2 and the average distance between i-th and j-th units 1 r ij ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi 6=pa 2 ji À jj p . The electrostatic energy of a single blob is given by The number of the blobs in the chain N blob = N/g = Na 2 /x blob 2 . Since the PE chain is stretched, the average distance between two blobs, k-th and l-th, is given by hr kl i = x blob |k À l|, 1 r k, l r N blob . The charge of each blob equals q = efg and the total electrostatic energy of interactions between all blobs reads Thus, the total electrostatic energy is given by while the elastic energy of deformation reads and L = x blob N blob = Na 2 /x blob . The total free energy of the chain can be expressed as a function of x blob which has to be minimized with respect to x blob : The exact value of x blob is determined by the equation: Taking into account Na 2 c x blob 2 and N c 1, one gets 71 and the total free energy can be estimated as This final result, which takes into account not only the meanfield free energy of a chain 31,32 but also contribution due to the correlation of charges (see Appendix B), should be used in calculations of the critical concentration of PEC formation. The scaling estimations for the single PE chain free energy coincide with the above results up to the logarithmic and numerical factors, which are usually omitted in scaling estimations. 56 However, it will be shown below that account for these corrections, 2.16(ln N) 2/3 , is crucial for equilibrium between the free chains and the PEC. A more common procedure for the self-energy calculation of a swollen PE chain in both dilute and semidilute salt solutions can be found in the recent work of Shen and Wang, where the authors took into account coupling between the chain conformation and the screening of Coulomb interactions in the solution. 57 Owing to the low polyion ionization degree, f { 1, we have also neglected chain entanglements (knots), which are known to be localized, adopt a tight configuration and relax very slowly (or even remain frozen) in the case of highly charged polyelectrolytes. 72

Critical concentration of PEC formation
In this section, we discuss the state of oppositely charged polyions in 1 : 1 solution. One can expect that in extremely dilute solutions polyions exist as single chains, while their increasing concentration should result in complexation and formation of either single globules on the basis of one polycation and one polyanion or a macroscopic condensed phase (coacervate). For the sake of simplicity, the equilibrium between the three states of the polyions, namely (i) single chains, (ii) neutral globules and (iii) macroscopic liquid coacervates will be considered, though formation of a solid precipitate is also possible in some cases, 73 e.g. for oppositely charged chiral polypeptides extra capable of hydrogen bonding. 73,74 Here we neglect formation of globules consisting of several couples of oppositely charged polyions assuming that their fraction will be smaller in comparison with the globules on the basis of one couple. In other words, we assume that the size distribution function of the globules has a narrow maximum for the one couple.
Let us assume that the chains are pretty long, N(uf 2 ln N/3) 2/3 c 1, i.e. the number of the electrostatic blobs within the chain exceeds unity by far. Free energy of the chains in each state can be divided into concentration-dependent entropic and concentration-independent energy contributions. The latter has the form: F coac = 2.17(uf 2 ) 2/3 NC 1/9 (62) for the single swollen and complexed chains, respectively. Here we assume that the single PE chains do not retain counterions since energy gain provided by the trapping of one counterions is N times lower than that as a consequence of PE complex formation. We also neglect the DH-type correlation attraction of these chains at an extremely small polymer concentration, which is negligible in comparison with their translational entropy, This estimation follows from the treatment of the solution as the DH plasma with ion charges equal to Q = ef N. Since the critical concentration of complex formation depends on the chain length N exponentially (see below), this condition is fulfilled for long enough chains. It is also assumed that inequality x { r c is fulfilled owing to the long enough PE chains. The energy contribution for the case of the globule includes the energy per chain of the complex coacervate and the excess surface energy per chain, 2pR gl 2 g. Here ga 2 /(k B T) = 0.11(uf 2 ) 2/3 C À7/18 is the surface tension coefficient of the globule. 32,75 This result takes into consideration connectivity of PE charges 32 and for this reason differs from that found in ref. 76, where combination of Vroon-Overbeek and Cahn-Hilliard theories neglecting polymer specificity of the system was used. Thus, free energy per chain in globule reads F gl = 2.17(uf 2 ) 2/3 NC 1/9 + 0.68(uf 2 ) 4/9 N 2/3 C À5/54 . (63) Denoted by M tot the total number of polycations (polyanions) in the solution of volume V. Let the fraction (1 À s) of polyions forms coacervate, and the total number of single polycations and polycations in the globules (i.e. in supernatant) equals M = sM tot . If all polyions in the supernatant are single chains, their free energy reads Formation of the globules can be treated as a reversible process of dimerization. Let us assume that fraction p of the supernatant polyions are involved in the globules. The free energy of the supernatant can be calculated via calculation of the partition function of the system, following the approach proposed in ref. 77. The state when all chains in the supernatant are free is the reference state with the corresponding partition function Z 0 , F 0 = Àk B T ln Z 0 . The partition function of the supernatant Z with pM globules reads W vol is the probability to find Mp polycations in the vicinity of corresponding Mp polyanions, namely, the center of mass of polyanions and polycations located within the volume of the order of the globule volume v gl : with v glob = 2Na 3 /F tot . The combinatorial factor P comb is the total number of ways to select Mp cations, and assemble Mp globules: Thus, the free energy of the supernatant F sn = Àk B T ln Z reads Finally, free energy of chains in the coacervate can be written as follows: The equilibrium distribution of the polyions between three states can be found via minimization of the total free energy F tot (s, p) = F sn + F coac with respect to p and s. It has to be noted that minimization with respect to p results in the law of mass action connecting concentrations of single chains and globules: The fractions of polyions in the coacervate, globules and non-aggregated state are shown in Fig. 2. In Fig. 3 we plot corresponding concentrations of polyions and globules in the supernatant. It is seen that increasing the polyion concentration first results in the formation of single globules. The single chains do not aggregate only at extremely small polymer concentrations. Their coexistence with the globules is observed also at extremely small polymer concentrations. The single globules dominate in the solution in a very wide concentration range, c = 10 À20 -10 À6 , above which the macroscopic coacervate is formed. These numerical results coincide well with analytical estimations of the critical concentrations of the globules c gl and coacervate c coac formation, which can be found from eqn (68)- (70):  Indeed, for long enough chains c gl { c coac as long as 2F gl À F coac À F sw o 0, and the formation of the globules starts at much lower concentrations than their precipitation in the solution. At chosen values of the numerical parameters, all assumptions adopted above are valid. Namely, polyions are long enough to contain a large number of the electrostatic blobs and they exist in the non-aggregated state only at extremely small concentrations. In the concentration range shown in Fig. 2 and 3, the coacervate volume is much smaller than the total volume of the system, so that disregard of the supernatant volume diminution in the course of precipitation is justified. Our principal result on the formation of single neutral interpolyelectrolyte globules prior to the macroscopic phase separation at increasing concentration of polyions in solution coincides with that recently reported by Delaney and Fredrickson. 78 They proved this fact by means of both scaling analysis and theoretical-field methods. At that, analytical expressions for threshold concentrations, c gl and c coac , found in the present paper seem to be more precise as compared to scaling estimations 78 where both the ln N correction as well as the surface effects were neglected, while strict theoretical-field analysis performed in their work did not allow finding analytical formulas.

Conclusion
Attraction of oppositely charged PEs within the complex is caused by fluctuation of charges and local deviation from electric neutrality. This correlation-induced attraction is opposed by repulsive volume interactions, and the equilibrium polymer volume fraction within the complex is defined from their balance. It is natural to expect that the correlation electrostatic free energy of the complex, F corr , should be negative, similar to the case of the DH plasma. However, direct calculations of the RPA contribution, F RPA , within the framework of a continuous model of charge distribution lead to the positive expression. This result is explained as follows. The formal calculations of the RPA Gaussian integrals take into account not only electrostatically induced correlations but also correlations caused by initial connectivity of like charges into the chain. The latter are caused by the structure of the polymer chain and should be included not into the correlation energy but into the self-energy of PE chains, F self .
The self-energy of polyions depends on the conformation of the chain, and different chain conformations can be used as a reference state. However, in the framework of the Borue-Erukhimovich model 25 which was widely used in the past few decades, the reference state is the Gaussian coil. This choice of the reference state is caused not only by mathematical (the structure of the free energy functional containing the conformational entropy term) but also by physical reasons. Indeed, a PEC in a Y-solvent is known to be a polymer globule, so that polyions within the globular complex/in coacervate adopt the conformation of a Gaussian coil on any length scales.
The self-energy of a PE chain within the complex equal to ffiffi ffi 3 p uf 2 N 3=2 is the Coulomb energy due to the correlations of charges along the polyion in ideal coil conformation. It is important to stress that the mean-field free energy of charged polyions with Gaussian coil statistics should not be included into their self-energy because of the overall complex electric neutrality (see Appendix A). The physical meaning of the chain self-energy can be defined as follows. It is the work which is necessary to simultaneously assemble an equal number of polycations and polyanions in the same volume (volume of PEC) providing them conformations of Gaussian coil. This work is required not to get non-zero overall charge of this volume (as in the case of single polyions, eqn (60)), since the complex is electrically neutral, but to achieve connectivity of charges into Gaussian coils. In turn, energy of fluctuationinduced attraction between these already assembled polyions is the true correlation energy, F corr .
For the first time within Borue-Erukhimovich polyelectrolyte complex model we split the RPA result for the energy of Coulomb interactions within the complex, F RPA , into two parts: self-energy of the chains, F self , and correlation energy of the complex, F corr . We demonstrate that the correlation free energy of the complex within the framework of this model considering continuous charge distribution, F corr = F RPA À F self , is negative, in accordance with general expectations. The proposed procedure for the RPA result renormalization (subtraction of the chain selfenergy) is based on physical arguments and uniquely defined. This procedure and PE chain self-energy definition are akin to those used formerly by Mahdi and Olvera de la Cruz 36 and later by Shen and Wang in the case of PE solutions. 57 Accounting for not only the mean-field term but also the correlation electrostatic free energy of the swollen single chain allows us to find the critical concentration c gl of the singleglobule PEC formation and concentration of precipitation of these globules, c coac . Analytical expressions for these concentrations are found for the case of pretty long chains providing a large number of electrostatic blobs. It is shown that the single globules of the PEC exist in a wide range of polyion concentrations, c gl { c coac , when complexation is already favourable owing to a high energy gain in the course of polyion neutralization, while precipitation is still unfavourable because of entropic reasons.

Appendix A. Mean-field and correlation energies of polyions with Gaussian statistics
Consider a bare charged PE chain with Gaussian statistics on any length scales assuming that all counterions leave the chain. Electrostatic energy of the chain can be calculated as follows: After substitution hF(r)F(r 0 )i = hF(r)ihF(r 0 )i + g chain (|r À r 0 |), the chain electrostatic energy can be represented as a sum of two terms: U G,ch el-st = U G,ch MF + U G,ch corr (74) Note that the correlation function of a single Gaussian chain g G,ch (|r À r 0 |) = g G,ch (R), R = |r À r 0 | coincides with G self (R) obtained earlier in eqn (42). This fact is a consequence of identical statistics of the single Gaussian chain and the chain in the solution of ideal (i.e. non-interacting) chains. The meanfield electrostatic energy written in the conventional form reads and omitting numerical factor one gets uQ 2 /r c with Q = Nf being the chain charge in the elementary charge e units and r c ¼ a ffiffiffiffi N p ffiffiffiffiffi 12 p being the characteristic Gaussian size of the chain. On the other hand, U G,ch el-st,MF = 0 within the complex since the excess charge of the PE molecule is compensated by the oppositely charged counterpart. The correlation term can be easily calculated: This result coincides with the RPA self-energy calculations, eqn (46), and supports the proposed interpretation for the self-energy being excess Coulomb energy caused by Gaussian correlations of linked chain charges. Account for the correlations is of great importance since both the mean-field self-energy of the Gaussian chain, U G,ch MF , and the self-energy correction caused by correlations, U G,ch corr , depend on the parameters in similar ways, both Buf 2 N 3/2 . They differ only by numerical coefficients.
Appendix B. Mean-field and correlation energies of the swollen polyion Let F 0 be the average polymer volume fraction within the sphere of the diameter BL circumscribed around the swollen PE chain (see Fig. 1), F 0 = Na 3 /L 3 = x el-st 3 /(Na 3 ) = u À1 f À2 N À2 where x el-st B au À1/3 f À2/3 is the electrostatic blob. To calculate the Coulomb energy of the swollen PE chain, we use the method similar to that presented in Appendix A in the case of the Gaussian charged chain. The electrostatic energy of the chain is given by and F(r) = F 0 + dF(r). Owing to hF(r)F(r 0 )i = hF(r)ihF(r 0 )i + hdF(r)dF(r 0 )i = F 0 2 + G sw,ch (r, r 0 ), one can find since hF(r)i = F 0 within the sphere circumscribed around the chain. The correlation electrostatic energy of the chain, U sw,ch corr , can be found owing to its known (within the scaling consideration) correlation function, G sw,ch (R). The intra-chain correlation function is given by since the chain is Gaussian at R o x el-st and linear at R 4 x el-st , and f B u 1/3 f 2/3 is the polymer volume fraction inside the blob. Then and e is the Euler's number. Thus, correlation electrostatic energy U sw,ch corr exceeds the mean-field term U sw,ch MF by a factor of ln N. For this reason account for the correlation term is requisite when charge of the chain is considered to be smeared throughout the sphere of diameter L.

Conflicts of interest
There are no conflicts to declare.