solutions of highly - charged polyelectrolytes , such as biopolymers ( dna , proteins , _ etc . _ ) and many synthetic polyelectrolytes , exhibit striking phenomena while multivalent salts or charged molecules are added into the solutions  @xcite .
if the concentration of the added salt lies within a window of concentration , the solution is separated in two phases : one is a viscous phase containing condensed or precipitated polyelectrolytes , and the other is a dilute phase in which polyelectrolytes are dissolved .
if it lies outside the window , the solution is in a homogeneous phase .
these phenomena are called `` reentrant condensation ''  @xcite , or commonly referred as `` dna condensation '' in the fields of biology and medicine  @xcite .
this topic has attracted much attention recently due to its potential application in developing nonviral vectors of dna for the purpose of gene therapy  @xcite .
the amount of salt which results in the phase separation is determined by many factors , such as valence of salt , temperature , chain structure , and concentration of polyelectrolytes  @xcite .
it is known that the lower boundary of the window depends linearly on the concentration of polyelectrolytes but the upper one of the window is insensitive to it and roughly a constant .
both flexible polyelectrolytes and stiff polyelectrolytes display similar phase diagrams  @xcite .    despite decades of effort ,
our understanding on polyelectrolytes is still limited and many issues , including the reentrant condensation , are far from being resolved  @xcite .
the aim of this paper is to investigate the properties of polyelectrolytes in salt solutions by means of computer simulations .
there have been many simulations devoted to the study of polyelectrolytes  @xcite . however , most of them investigated the properties in salt - free solutions  @xcite . only recently ,
due to the progress of computing power , large - scale simulations for systems with added salt can be realized .
the works included , for example , flexible chains in monovalent or multivalent salt solutions  @xcite , semiflexible chains with multivalent salt or charged molecules  @xcite , and rigid chains in the presence of multivalent salts  @xcite .
but the salt concentration or the salt valence was not large enough to observe any significant evidence being able to connect to the redissolution of polyelectrolytes . therefore , the reentrant condensation could not be discussed thoroughly in these works .
et al . _
@xcite recently studied the effective interaction between two immobile dnas in the presence of tetravalent counterions .
they have found that the effective interaction yields an attractive minimum at short distances which then disappears in favor of a shallow minimum at large separations , indicating the occurrence of dna condensation and redissolution and a novel stable mesocrystal . in a more recent study  @xcite , we have investigated the properties of mobile polyelectrolytes in tetravalent salt solutions and demonstrated that following a collapsed transition , polyelectrolytes undergo reexpanding transition upon addition of the salt .
this phenomenon is just a single - chain version of the condensation and the redissolution of polyelectrolytes .
moreover , by varying the ion size , we have shown that the excluded volume of ions plays a decisive role on the reentrant condensation .
however , to well understand the phenomena , there is still a far way to go .
for example , is it possible to obtain any information related to the phase boundary of reentrant condensation from simulations ?
there exist theories to explain the phase boundary  @xcite but it has never been studied in simulations . concerning the size of polyelectrolyte in multivalent salts , muthukumar and coworkers have applied single and double screening theories to explain the collapsed transition of chains  @xcite .
however , no existing theory can predict chain reexpansion in the region of high salt concentration . moreover , chain morphology is essential to understand the properties of polyelectrolytes .
it can be described by many quantities , such as asphericity and prolateness parameter  @xcite . to our knowledge
, no such study has been reported yet for strongly - charged polyelectrolytes .
an alternative way to understand polymers is to study scaling behavior .
although the concept of scaling has been successfully applied to neutral polymers  @xcite , application of this concept to polyelectrolytes  @xcite has encountered many difficulties owing to the long - range nature of coulomb interaction . in this paper , we compute the single - chain structure factor .
the swelling exponent is then calculated by studying the power - law - like regime of the structure factor .
since properties of polyelectrolytes depend strongly on ion size , we investigate in this paper the effect of ion size too . in our previous work
@xcite , the sizes of ions were set to identical and varied simultaneously . as a consequence , bjerrum association  @xcite
was strongly intervened , particularly when ions were small .
this setup has a disadvantage that the system was easily jammed by large ions , leading to the impossibility of investigation of the cases of big ions .
therefore , in this study we vary only the size of the multivalent counterions and fix unchanged the one of the other ions , so that the ion size effect is restrictedly investigated .    in the study of polyelectrolytes , it is tempting to discuss the long range nature of coulomb interaction using the concept of electrostatic persistence length .
odijk  @xcite , skolnick and fixman  @xcite predicted that the electrostatic persistence length is proportional to the square of the screening length @xmath0 .
the deduction was based upon debye - hckel theory .
consequently , the effective intramolecular interaction is always repulsive , resulting in a positive electrostatic persistence length  @xcite .
however , experiments have shown that electrostatic persistence length can be negative under certain condition so that electrostatics gives a reversed contribution to chain flexibility  @xcite .
simulations on this issue  @xcite were mainly performed within debye - hckel approximation using a screened coulomb potential @xmath1 to describe the screening effect of salt without involving explicitly salt ions . a negative electrostatic persistence length ,
hence , has never been reproduced in the simulations .
noticeably , it is known that linear screening theory is valid only for a weakly - charged system and breaks down for a highly - charged one , such as dna and many biopolymers , due to strong ion correlations  @xcite . recently , by employing a loop expansion method which goes beyond mean field , ariel and andelman were able to predict the negative regime of the electrostatic persistence length for a strongly - charged rodlike polyelectrolyte  @xcite .
therefore , we expect that to observe a negative electrostatic persistence length , one should incorporate salt ions explicitly and use the bare coulomb potential in simulations .
our study fulfills this requirement and hence offer a good chance to verify the above phenomenon and theory .
charge inversion is a universal phenomenon , occurred when a strongly charged macroion binds so many counterions that its net charge alters sign  @xcite .
_ intended to link charge inversion with the reentrant condensation of polyelectrolytes  @xcite .
the idea was based upon that counterions form strongly correlated liquid on the surface of polyelectrolytes . in an intermediate salt region , the bare charges of the chains are almost neutralized by the condensed counterions ; the correlation - induced attraction results in the condensation of polyelectrolytes . in a high - salt region ,
the counterions overcompensate the chain bare charges ; coulomb repulsion between chains dominates the attraction and , in consequence , the condensates are redissolved into the solution .
solis and olvera de la cruz gave a different point of view using a theory called `` two - state model ''  @xcite .
they predicted that charge inversion is not necessarily occurred with redissolution of polyelectrolytes .
although overcompensated by counterions , the redissolved polyelectrolytes can have either sign due to the association of coions  @xcite .
this association depends strongly on the ion size which was not considered in the study of nguyen _ et al .
_  @xcite .
therefore , solis and olvera de la cruz predicted a charge distribution which oscillates from a polyelectrolyte , in accordance with the profile observed in our previous work  @xcite . in this paper , we investigate more extreme cases , including much smaller and much larger ions .
we also study like - charge attraction between chains by calculating the potential of mean force and determine the conditions under which the attraction takes place .
the rest of the article is organized as follows .
model and simulation details are given in sec .
[ sec_model_method ] .
results are discussed in sec .
[ sec_result_discussion ] .
it is divided into five subsections .
the first subsection ( sec .  [ sec_rg2_cm ] ) is devoted to the study of the radius of gyration of a polyelectrolyte where salt concentration and polyelectrolyte concentration are varied .
the results provide information helpful in understanding the boundaries of the condensation window .
[ sec_shape ] studies the morphology of multiple polyelectrolytes in salt solutions .
we compute several quantities , including shape factors , asphericity , and prolateness parameter .
the effect of size of multivalent ions is discussed .
single - chain structure factor is , then , calculated in sec .
[ sec_sf ] .
by least - square fitting the power - law - like regime , we compute the swelling exponent of a polyelectrolyte , which provides a fundamental vision how strongly - charged polyelectrolytes scale in salt solutions . sec .
[ sec_lp ] deals with persistence length .
the negative regime of the electrostatic persistence length is successfully reproduced in the simulations .
comparison of the results with theoretical predictions is made .
[ sec_pmf ] provides a general picture of the effective interaction between polyelectrolytes by computing the potential of mean force .
as long as knowing the condition to appear like - charge attraction between chains and the integrated charge distribution around a chain ( sec .
[ sec_iondist ] ) , we are able to arbitrate if the phenomena of charge inversion and the redissolution of polyelectrolytes are related .
we give our conclusions in sec .
[ sec_conclusion ] .
we employed a coarse - grained model to simulate a polyelectrolyte system .
the system contains four bead - spring chains , one of which is composed of 48 beads ( monomers ) .
each monomer dissociates a monovalent cation ( counterion ) into a solution and carries a negative unit charge .
there are totally 192 monovalent cations in the solution .
solvent molecules are not modeled explicitly .
their effect is taken into account implicitly by considering them as a medium of constant dielectric constant .
the geometry of the medium is cubic and periodic boundary condition is applied to simulate bulk environment .
salt molecules are added into the system . in the solution , they are dissociated into cations ( counterions ) and anions ( coions ) .
we assume that the cations are tetravalent and the anions are monovalent .
the amount of the cations and anions obeys the condition of charge neutrality .
three types of interactions are included in our simulations .
the first one is the excluded volume interaction described by a purely repulsive lennard - jones ( lj ) potential , @xmath2 + \varepsilon_{\rm lj } & \mbox{\ for $ r \le \root 6 \of 2 \sigma$\;}\\ 0 & \mbox{\ for $ r > \root 6 \of 2 \sigma$\ ; } \end{array } \right.,\ ] ] where @xmath3 is the distance between two particles , @xmath4 is the interaction strength , and @xmath5 is the collision diameter .
we assume that all the particles have identical @xmath4 but the collision diameter may be different .
the collision diameters for monomer , monovalent cation , tetravalent cation , and anion are denoted by @xmath6 , @xmath7 , @xmath8 , and @xmath9 , respectively .
lorentz - berthelot mixing rule  @xcite is applied for the interaction between different kinds of particles .
the second interaction is called finitely extensible nonlinear elastic ( fene ) potential  @xcite and used to describe the bond connection between adjacent monomers on a chain .
it reads as @xmath10 where @xmath11 is the spring constant , @xmath12 is bond length , and @xmath13 is the maximum extension of a bond .
all particles interact via the third interaction , coulomb interaction , which is written as @xmath14 where @xmath15 is the charge of the @xmath16th particle , @xmath17 is the temperature , @xmath18 is the boltzmann constant , and @xmath19 is the bjerrum length defined as the separation distance at which the coulomb potential between two unit charges @xmath20 equals to the thermal energy @xmath21 .
@xmath19 is equal to @xmath22 where @xmath23 is the dielectric constant and @xmath24 is the vacuum permittivity .
we perform molecular dynamics ( md ) simulations in this study  sjplimp / lammps.html . ] .
we set @xmath7 and @xmath9 equal to @xmath6 and vary the size of tetravalent cation @xmath8 from @xmath25 to @xmath26 to investigate the effect of ion size on properties of the polyelectrolyte solution .
simulation parameters are chosen as follows : @xmath27 , @xmath28 , and @xmath29 .
the bjerrum length is set to @xmath30 .
similar systems in salt - free solutions have been investigated by stevens and kremer  @xcite .
we assume that the masses of all the particles are identical and equal to @xmath31 , and choose @xmath6 and @xmath32 to be the length unit and the energy unit , respectively .
therefore , the natural time unit is @xmath33 .
the system is studied in canonical ensemble where the temperature is controlled using nos - hoover thermostat with the frequency of heat transfer with thermostat equal to @xmath34 .
coulomb interaction is calculated by the technique of ewald sum . except in sec .
[ sec_rg2_cm ] , the monomer concentration is fixed at @xmath35 which corresponds to a box size @xmath36 , whereas the salt concentration is varied .
the amount of salt is controlled such that the volume fraction of total particles is smaller than @xmath37 .
equations of motion are integrated applying verlet algorithm where time step is set to @xmath38 .
each simulation starts with an equilibration phase of @xmath39 md steps and is followed by a production run of @xmath40 md steps in which data are collected every 1000 steps . in order to shorten the notation
, we will use @xmath6 as the length unit , @xmath41 as the concentration unit , and @xmath20 as the charge unit in the following text .
we started our study by simulating a single polyelectrolyte with added tetravalent salt at various monomer concentrations .
the chain consists of 48 monomers and the size of the tetravalent cations is set to be identical to the one of the monomers .
we used the radius of gyration to characterize the chain size .
the radius of gyration @xmath42 is calculated by the equation @xmath43 where @xmath44 is the number of monomers on the chain , @xmath45 is the position vector of the @xmath16th monomer , and @xmath46 the center of mass of the chain .
the results are presented in fig .
[ cm_p1n48_rg2.eps ] where the mean square radius of gyration @xmath47 is plotted as a function of salt concentration @xmath48 at four monomer concentrations , @xmath49 , @xmath50 , @xmath51 and @xmath52 .     as a function of @xmath48 at four monomer concentrations , @xmath49 , @xmath50 , @xmath51 , and @xmath52 .
the symbols used to represent the four @xmath53 are indicated at the left - bottom side of the figure .
the value of @xmath47 in a salt - free solution ( @xmath48=0.0 ) is depicted near the left vertical axis .
the four curves in the regime @xmath54 can overlap each other ( shown in the inset ) under a suitable transformation described in the text .
, scaledwidth=40.0% ]    the angular bracket @xmath55 means averaging over chain configurations .
two behaviors were observed .
the first one is that each curve displays a similar trend of behavior against salt concentration .
upon addition of salt , the chain undergoes two structural transitions : a collapsed transition and a followed reexpanding transition , as having been observed in the previous work  @xcite .
the first transition results from the screening of the coulomb interaction by the added salt . the second transition is still not well understood although it is a result of the competition between energetic gain and entropic loose in counterion condensation and in shape transformation of polyelectrolytes to minimize the free energy .
we observed that the minimum of @xmath47 appears near the salt concentration @xmath56 and the more dilute the monomer concentration , the smaller value the minimum of @xmath47 will be .
@xmath57 is the equivalence point at which the amount of charges of the tetravalent counterions exactly neutralizes the charges on polyelectrolytes as if the salts and the polyelectrolytes are in the stoichiometric reaction .
we will call @xmath57 the equivalence concentration . in a salt - free solution ( @xmath58 )
, we found that the chain size increases with decreasing monomer concentration ( refer to the four data points near the left vertical axis of the figure ) .
this effect has been discussed in literatures  @xcite and our results are in accordance with them .
and the second behavior is that @xmath47 for different @xmath53 attains approximately the same value while the salt concentration is larger than the equivalence concentration .
therefore , the curves approximately overlap each other .
we noticed , furthermore , that in the region @xmath54 , the @xmath47 curves for different @xmath53 can also lie over each other under an appropriate transformation .
we demonstrate this overlap in the inset of fig .
[ cm_p1n48_rg2.eps ] where the abscissa is @xmath59 and the ordinate is @xmath60 with @xmath61 the minimum value of @xmath47 .
@xmath62 depends smoothly on @xmath53 and , in the figure , is equal to @xmath63 , @xmath64 , @xmath65 , and @xmath66 for @xmath49 , @xmath50 , @xmath51 , and @xmath52 , respectively .
the size of polyelectrolytes in multivalent salt solutions has been studied by small - angle scattering experiments  @xcite .
the results showed that it diminishes with increasing or decreasing @xmath48 toward the region where the condensation of polyelectrolytes occurs ( cf .
table 2 of ref .
the results of our simulations are in accordance with the experiments .
we make a hypothesis that chain size can be served as a qualitative condition to distinguish between a monophasic solution of polyelectrolyte and a diphasic one .
more precisely , if the chain size is larger than some value , the chains are dissolved in the solution and the system is in a homogeneous phase ; if it is smaller , the chains are precipitated from the solution and a phase separation occurs . under this hypothesis ,
polyelectrolyte solutions of different @xmath53 become monophasic as @xmath67 is increased over a concentration in the high - salt region , since the chain size exceeds the some value nearly at the same @xmath48 due to the similarity of the @xmath47 curves against @xmath48 in this region .
therefore , the upper boundary of the salt concentration window for the condensation of polyelectrolyte is independent of @xmath53 .
similarly , the lower boundary of the window approximately depends linearly on @xmath53 , owing to the similarity of the @xmath60 curves against @xmath59 in the low - salt region ; the weak dependence of @xmath62 and @xmath61 on @xmath53 ensures that the chain size surpasses the some value roughly at the same @xmath59
. therefore , fig .
[ cm_p1n48_rg2.eps ] provides an essential information to understand the boundaries of reentrant condensation .
the aim of this article is to understand the properties of multiple flexible polyelectrolytes with added salt in an aqueous solution .
it is known that the bjerrum length @xmath19 in water is equal to 7.14 at room temperature .
our simulation setup , hence , corresponds to a collision diameter @xmath6 of monomer equal to 2.38 .
this model represents a polyelectrolyte system with linear charge density equal to a unit charge per @xmath68 .
] , and can be employed to simulate a strongly - charged system such as sodium polystyrene sulfonate .
the four monomer concentrations studied in this section , while converted to real unit , are equal to 0.039 m , 0.246 m , 0.986 m , and 1.971 m , respectively .
even the largest @xmath53 is still smaller than the overlap threshold @xmath69 estimated by @xmath70 .
therefore , our systems are dilute polymer solutions .
readers should be aware that the value up to nearly 2 m is a very high monomer concentration for a real polymer system , and our systems are _ dilute _ only in the formal meaning of the world since the chains are very short . in the following sections , we fix @xmath53 at 0.008 and simulate a system comprising multiple chains .
this monomer concentration is one or two order of magnitude higher than a typical monomer concentration used in experiments  @xcite .
however , based on the fact that systems at different monomer concentrations reveal similar evolution of chain size against salt concentration , the results obtained at @xmath71 can be used to understand a system at more dilute @xmath53 .
we mention that this @xmath53 is a choice so that the simulation box is not so big .
therefore , the amount of the added salt is controllable under limited computing resources .
consequently , to study the behavior of polyelectrolytes in salt solutions , covering a broad range of salt concentration , becomes numerically feasible .      in this section and in the followings ,
the studied system contains four polyelectrolytes and the monomer concentration is fixed at @xmath71 .
this section is devoted to the study of chain conformation in tetravalent salt solutions .
the effects of salt concentration and size of tetravalent counterions are investigated .
we vary the size of tetravalent counterion , @xmath8 , from 0.0 to 4.0 and fix the size of the other ions at 1.0 .
a large @xmath8 can represent a large cation , a large ion group , or a large charged colloid .
it can be also a result of the hydration of a small ion .
a small @xmath8 may be less representative for real systems but theorists are generally interested in it because many models are established on the base of vanishing ion size .
morphology of a polymer can be quantified by many quantities .
the first quantity that we studied is the shape factor @xmath72 , defined as the ratio of the mean - square end - to - end distance , @xmath73 , to the mean - square radius of gyration @xmath47 .
@xmath72 attains a value 12 for a rodlike polymer , 6.3 for a flexible chain in a good solvent , and 6 for an ideal chain  @xcite . for a compact structure ,
the value of @xmath72 is small . for example , if the chain has a spherical morphology and the two ends are randomly positioned , @xmath72 is equal to 2 . in fig.[shapep4n48_eta1&xi1.eps](a ) , we show how the shape factor @xmath72 varies with @xmath48 at different @xmath8 .     as a function of @xmath48 , and ( b ) @xmath74 as a function of @xmath48 , for the cases of different size @xmath8 of the tetravalent counterions . the symbol used for each @xmath8
is indicated at the left - bottom side of the figures .
the error bars are smaller than the size of the symbols .
the result for neutral polymers is depicted as reference and marked by `` neutral '' .
, title="fig:",scaledwidth=40.0% ] ( a )   as a function of @xmath48 , and ( b ) @xmath74 as a function of @xmath48 , for the cases of different size @xmath8 of the tetravalent counterions .
the symbol used for each @xmath8 is indicated at the left - bottom side of the figures .
the error bars are smaller than the size of the symbols .
the result for neutral polymers is depicted as reference and marked by `` neutral '' . , title="fig:",scaledwidth=40.0% ] ( b )    in the figure , a reference curve , marked by `` neutral '' , is depicted .
it represents a system composed of four neutral polymers of 48 monomers in tetravalent salt ( @xmath75 ) solutions .
focus firstly on the polyelectrolyte solutions with @xmath76 , 1.4 and 2.0 .
we witnessed that @xmath72 decreases at first and then increases upon addition of salt , in analogy of the radius of gyration studied in the previous section , and shows a v - shaped curve against @xmath48 in the semilog plot . in the low - salt region or in the high - salt region , @xmath72 is larger than or close to the value 6.3 , indicating that the chains have elongated or coil - like structures . in between , it appears a minimum near the equivalence concentration @xmath77 .
the value of @xmath72 tells us that the chains possess compact structures but not so compact as spheres . for the cases with larger @xmath8 ( @xmath78 and @xmath79 )
, @xmath72 shows a shallow minimum and the whole curve stays basically above the value 6.3 . in these cases ,
the chains are in an extended state .
it is worth to notice that for @xmath80 , @xmath72 , in principle , increases with @xmath8 at a given @xmath48 .
now pay attention to the cases with @xmath81 .
@xmath72 displays a similar v - shaped curve against @xmath48 but opposite to the previous cases , its value increases with decreasing @xmath8 .
moreover , the degree of chain reexpansion in the high - salt region is weakened while @xmath8 is decreased . noticeably , for the case @xmath82
, @xmath72 displays a sharp transition from a value 7.2 in the low - salt solution to a value 6.5 in the vicinity of @xmath57 . in the latter region ,
the chains behave similarly to neutral polymers .
the second quantity that we studied is @xmath83 , which is the ratio of the radius of gyration @xmath84 to the hydrodynamic radius @xmath85 computed by @xmath86 .
@xmath83 is equal to @xmath87 for a hard sphere and attains a value 2.25 for a rodlike chain  @xcite .
renormalization group theory has predicted @xmath88 in a good solvent and @xmath89 at @xmath90 point  @xcite , both of which have been shown in agreement with experiments  @xcite .
we present in fig .
[ shapep4n48_eta1&xi1.eps](b ) the variation of @xmath83 against @xmath48 .
the referenced curve for the system of neutral polymers is depicted .
the curve was found very close to the value 1.24 , which seems to indicate an ideal - chain behavior .
however , the neutral chains should behave as random coils and the value of @xmath83 would be larger than 1.24 .
the deviation of @xmath83 from the value of a coil is probably due to the short chain length used in our simulations .
_ have studied the finite size effect and shown that the dependence of @xmath83 on chain length is not negligible  @xcite .
therefore , @xmath83 reported here departs slightly from the asymptotic value of an infinite chain length . nonetheless , the results are still indicative and useful to understand how @xmath83 varies with salt concentration . in the small @xmath48 region ,
the polyelectrolytes are more extended than the neutral polymers . in the middle region of @xmath48 ,
the chains attain compact structures when @xmath91 .
the dependence of @xmath83 upon @xmath8 resembles to that of @xmath72 . for large @xmath8 ( @xmath92 ,
the chains are elongated . for vanishing @xmath8 , we observed that following a sharp decrease , @xmath83 reattains approximately a constant value , close to that for the neutral polymers .    a more fundamental way to study
the shape of a polymer is to investigate the tensor of radius of gyration @xmath93  @xcite , defined by : @xmath94 where @xmath95 , @xmath96 denote the three cartesian components and @xmath97 is the @xmath95 component of the position of the @xmath16th monomer . let @xmath98 , @xmath99 and @xmath100 be the three eigenvalues of @xmath101 .
a quantity called `` asphericity ''  @xcite , which measures the deformation from a spherical geometry , is defined by @xmath102 it takes a value between 0 ( sphere ) and 1 ( rod ) . for a coil polymer , it is 0.431 obtained from simulations  @xcite .
another quantity describing the degree of prolateness is defined by @xmath103 where @xmath104  @xcite .
it is ranged from @xmath105 to @xmath106 and can be used to distinguish between a prolate shape ( @xmath107 ) and an oblate one ( @xmath108 ) . for neutral polymers in the dilute limit , simulations predict a prolate chain with @xmath109  @xcite .
we present in figs .
[ shapep4n48_asph&sp.eps](a ) and [ shapep4n48_asph&sp.eps](b ) the asphericity @xmath110 and the prolateness parameter @xmath111 as a function of @xmath48 .     as a function of @xmath48 , and ( b ) @xmath111 as a function of @xmath48 , for different @xmath8 .
the symbol used for each @xmath8 is indicated at the left - bottom side of the figures .
the error bars are smaller than the size of the symbols .
the result for neutral polymers is depicted as reference and marked by `` neutral '' .
, title="fig:",scaledwidth=40.0% ] ( a )   as a function of @xmath48 , and ( b ) @xmath111 as a function of @xmath48 , for different @xmath8 .
the symbol used for each @xmath8 is indicated at the left - bottom side of the figures .
the error bars are smaller than the size of the symbols .
the result for neutral polymers is depicted as reference and marked by `` neutral '' . , title="fig:",scaledwidth=40.0% ] ( b )    resembling to the behaviors of @xmath72 and @xmath83 , @xmath110 and @xmath111 show a minimum value in the vicinity of @xmath57 .
the curve associated to @xmath75 roughly separates two different behaviors of evolution : @xmath110 and @xmath111 increase with increasing and with decreasing @xmath8 from 1.0 .
the value of @xmath110 lies between 0.275 and 0.55 .
therefore , even in a compact structure ( where @xmath48 is near @xmath57 ) , the chain shape is asymmetric . in a low - salt solution ,
the chains show extended morphology but are not as elongated as rods . on the other hand , @xmath111 is larger than zero for all of the studied cases , indicating that the polyelectrolytes display prolate conformation . at a high salt concentration ,
@xmath110 is close to 0.45 and @xmath111 is close to 0.54 , suggesting once more that the polyelectrolytes behave as the neutral polymers in the high - salt region . for the case with null @xmath8
, chain reexpansion does not occur ; @xmath110 and @xmath111 stay at a constant value when @xmath48 goes beyond @xmath57 .
the results obtained in this section strongly suggest that @xmath75 is the optimal condition to pack polyelectrolytes into the smallest volume .      in this section
, we investigate the single - chain structure factor .
this quantity is experimentally accessible and provides an essential information to theoretical analysis .
the single - chain structure factor is defined by @xmath112 where @xmath113 is the scattering wavevector , @xmath114 is the number of polymer chains , @xmath115 is the number of monomers on a chain , and @xmath116 is the position of the @xmath117th monomer on the @xmath118th chain .
suppose that the system has spherical symmetry .
the dependence of the structure factor on the orientation of @xmath113 can be therefore integrated out : @xmath119 where @xmath120 is the norm of @xmath113 and @xmath121 is the solid angle .
this integration was performed numerically by sampling over a set of randomly - oriented @xmath113 vectors .
we present , in fig .
[ sfp4n48_sxxx&txxx](a ) , @xmath122 at several salt concentrations for the case of @xmath75 . for the reason of clarity ,
the curves in the figure have been multiplied by @xmath123 with @xmath124 being increased by one for each curve from the bottom to the top .
@xmath122 for the neutral polymers is depicted as reference with @xmath125 .     for @xmath75 at different @xmath48 and ( b ) @xmath122 at @xmath126 for different @xmath8 .
the value of @xmath48 or @xmath8 , at which the simulations were running , is indicated at the left side of each corresponding curve .
@xmath122 for the neutral polymers at @xmath126 is depicted as reference and marked by `` neutral '' .
for clarity , the curves , from bottom to top , have been multiplied by a factor @xmath127 with @xmath124 started from zero and increased by one for every curve .
the dotted curve , near each @xmath122 and terminated at @xmath128 , is the guinier function @xmath129 .
two dashed lines , which represent two power - law dependences , @xmath130 and @xmath131 ( marked by `` -1 '' and `` -2 '' , respectively ) , are drawn in the plots.,title="fig:",scaledwidth=40.0% ] + ( a ) +   for @xmath75 at different @xmath48 and ( b ) @xmath122 at @xmath126 for different @xmath8 .
the value of @xmath48 or @xmath8 , at which the simulations were running , is indicated at the left side of each corresponding curve .
@xmath122 for the neutral polymers at @xmath126 is depicted as reference and marked by `` neutral '' . for clarity , the curves , from bottom to top ,
have been multiplied by a factor @xmath127 with @xmath124 started from zero and increased by one for every curve .
the dotted curve , near each @xmath122 and terminated at @xmath128 , is the guinier function @xmath129 .
two dashed lines , which represent two power - law dependences , @xmath130 and @xmath131 ( marked by `` -1 '' and `` -2 '' , respectively ) , are drawn in the plots.,title="fig:",scaledwidth=40.0% ] + ( b ) +    we witnessed that @xmath122 shows typical behavior . in the guinier regime ( @xmath132 )
, @xmath122 behaves similarly to @xmath129 .
it is demonstrated by plotting @xmath129 ( in dotted curve ) near the corresponding @xmath122 with the value of @xmath133 adopted from sec .
[ sec_shape ] .
the dotted curve terminates at @xmath128 in order to indicate the boundary of the guinier regime . in the regime @xmath134
, @xmath122 shows power - law - like behavior , manifested by a linear dependence of @xmath135 on @xmath136 .
the slope @xmath137 of the linear dependence is related to the swelling exponent @xmath138 by the relation @xmath139 , which describes the size extension of a polymer : @xmath140 . in the regime @xmath141 ,
@xmath122 displays non - universal behavior . however , while @xmath142 , self - scattering of monomer is the only contribution , leading , therefore , @xmath122 to a universal value @xmath63 .
we mention that the krathy regime , in which @xmath143 , is not observed in our study because the system is a dilute solution and the chain length is not long enough  @xcite .
for @xmath8 other than 1.0 , we observed analogous trends of behavior .
for example , the variation of @xmath122 for different @xmath8 at the equivalence concentration is presented in fig .
[ sfp4n48_sxxx&txxx](b ) , in which we rediscovered the three regimes : the guinier regime , the power - law - like regime and the non - universal regime .
we noticed that in the last regime , a small bounce appears on the curve and , with increasing the size of tetravalent ions , gradually moves toward small @xmath120 .
the swelling exponent @xmath138 of the polyelectrolytes was calculated by performing least - square fits in the log - log plot of @xmath122 in the power - law - like regime .
the fitting region was chosen to be @xmath144 and the fitting reliability has been verified to be excellent .
we present the results of the swelling exponent @xmath138 in fig .
[ nup4n48.eps ] .     as a function of @xmath48 for different @xmath8 .
the symbol used for each @xmath8 is indicated at the left - bottom side of the figure .
@xmath138 for the referenced neutral - polymer system is depicted and marked by `` neutral '' .
, scaledwidth=40.0% ]    @xmath138 for the referenced neutral polymers was verified at first and found to attain a value close to the flory s one , @xmath145 .
next , we focused on the swelling exponent for the polyelectrolytes .
we observed that @xmath138 exceeds the value for the neutral polymers at a low salt concentration .
it indicates a strong extension of chain size due to coulomb repulsion between monomers . for the cases with @xmath146 , 1,0 , 1,4 , and 2.0
, the swelling exponent shows a v - shaped curve against @xmath48 : @xmath138 decreases and then increases with increasing salt concentration .
a minimum occurs near the equivalence concentration @xmath147 and , at this moment , the chains show compact structures but are less compact than spheres because a spherical structure scales as @xmath148 .
we remark that the swelling exponent obtained here for @xmath75 is consistent with that obtained by direct variation of the chain length of a system containing a single polyelectrolyte  @xcite except in the mid - salt region where @xmath138 in the latter study is slightly smaller and acquires a value @xmath149 at @xmath147 .
this deviation is due to multi - chain aggregation occurred in our system , which deforms the chain shape from a sphere and , consequently , increases @xmath138 . for the case with vanishing @xmath8 , following a drastic decrease
, @xmath138 retains roughly a value 0.6 in the high - salt region and the chains , hence , scale like the neutral polymers . for the cases with large @xmath8 ( @xmath78 and 4.0 ) , @xmath138 is greater than 0.6 but still less than 1.0 .
the chains , hence , do not reveal rigid rod - like structures . in principle , @xmath75 separates two behaviors of evolution at a salt concentration : @xmath138 increases with increasing and with decreasing @xmath8 from 1.0 .
the results demonstrate again that the polyelectrolytes are mostly condensed while the size of the tetravalent counterions is compatible with the one of the monomers .
there is a considerable interest to understand the variation of the persistence length of a polyelectrolyte against salt concentration .
persistence length is a local property , which measures the stiffness of a polymer , and is distinguished from global conformational quantities .
there are many ways to define the persistence length  @xcite . in this section ,
we adopt a microscopic definition to calculate it : @xmath150 where @xmath151 is the @xmath117th bond vector and @xmath152 is the root - mean - square bond length . in this definition ,
the contribution of the finite bond length , which corresponds to the @xmath153 term in eq .
[ persist_lp ] , is considered .
the obtaining @xmath154 as a function of @xmath48 for different @xmath8 is plotted in fig .
[ lpp4n48.eps ] .     as a function of @xmath48 for different @xmath8 .
the symbol used for each @xmath8 is indicated at the left - bottom side of the figure .
the error bars are smaller than the size of the symbols .
the bare persistence length @xmath155 is depicted as reference and marked by `` neutral '' .
, scaledwidth=40.0% ]    resembling to the quantities of chain conformation , @xmath154 decreases in the chain collapse region and increases in the chain reexpansion region .
a minimum value appears in the vicinity of @xmath57 .
the degree of increase diminishes with decreasing @xmath8 . for vanishing @xmath8 ,
it is so weak that @xmath154 is roughly a constant as @xmath156 .
similar to the results in the previous sections , @xmath154 increases with both increasing and decreasing @xmath8 from 1.0 .
@xmath75 , therefore , denotes a condition to separate the two behaviors of evolution of @xmath154 .
the persistence length has generally two contributions : the bare persistence length @xmath155 ( in the absence of electrostatic force ) and the electrostatic persistence length @xmath157 ( owing to the coulomb interaction ) @xcite . in this study ,
the bare persistence length @xmath155 was calculated by performing simulations on the referenced system where four neutral chains and 192 neutral particles are placed in salt solutions .
the result is depicted in fig .
[ lpp4n48.eps ] and the curve is marked by `` neutral '' .
we found that @xmath155 is not sensitive to the salt concentration and roughly attains a value 2.3 .
it slightly decreases in the high - salt region , probably due to the jamming effect by the salt ions . for the polyelectrolyte solutions with large @xmath8 ( @xmath78 and 4.0 ) , the @xmath154 curve thoroughly lies above the @xmath155 curve , indicating a more rigid flexibility than a neutral polymer .
the electrostatic persistence length is , hence , positive . on the other hand , for the polyelectrolyte solutions with @xmath158 , @xmath154 is smaller than @xmath155 in the mid - salt region and shows one or two intersections with the @xmath155 curve .
therefore , @xmath159 is positive in a salt - free solution and becomes eventually negative upon addition of salt .
it increases in the high - salt region and reattains a positive value for @xmath160 or @xmath161 .
the negative contribution of electrostatics to the persistence length has been experimentally observed  @xcite .
recently ariel and andelman ( aa ) pointed out that a negative @xmath157 indicates a mechanical instability ( collapse ) of polyelectrolytes due to the presence of multivalent counterions  @xcite . by taking into account the thermal fluctuations and correlations between bound counterions ,
they derived an expression of @xmath157 for an intrinsically - rigid charged polymer in ( @xmath162:1)-salt solutions in the limit of @xmath163 : @xmath164 where @xmath165 is the coulomb strength parameter between a monomer and a @xmath162-valent counterion , @xmath166 is the inverse debye - hckel screening length , and @xmath167 is the electrostatic persistence length deduced by odijk - skolnick - fixman ( osf ) theory  @xcite . in the derivation , they have assumed point - like charges for the ions .
we present , in fig .
[ lep4n48.eps ] , @xmath157 obtained from the simulations and that predicted by eq .
( [ le_ariel ] ) as a function of @xmath168 .     as a function of @xmath168 for different @xmath8 .
the symbol used for each @xmath8 is indicated at the left - bottom side of the figure .
the error bars are smaller than the size of the symbols .
the dashed curve depicts @xmath169 predicted by eq .
( [ le_ariel ] ) .
, scaledwidth=40.0% ]    we witnessed that @xmath157 depends strongly on the size of the tetravalent ions
. the negative regime of @xmath157 has been successfully captured in the simulations . to make a comparison with the theoretical prediction
, we focused on the case with vanishing @xmath8 .
please notice that this case does not exactly fulfill the assumption of the theory because in the simulations , the monovalent counterions and coions are not point charge - like and , also , the chains are not rigid .
however , according to our experiences  @xcite , ions other than multivalent counterions do not play a decisive role on the properties of polyelectrolytes .
this comparison , therefore , makes sense if the effect of chain rigidness is negligible .
we emphasize that the electrostatic correlations are fully taken into account in the simulations . in our study ,
the region in which eq .
( [ le_ariel ] ) is valid is @xmath170 because @xmath171 , @xmath172 , and @xmath173 is about 1.1 .
notice that eq .
( [ le_ariel ] ) has a singularity at @xmath174 and , hence , the curve is divided into two branches .
the left branch is located in the valid region but fails to predict @xmath157 .
it gives not only a wrong value but also a wrong trend of variation of @xmath157 against @xmath168 .
we noticed that eq .
( [ le_ariel ] ) is a modification of @xmath175 multiplied by a factor but the valid region given by aa theory  @xcite is incompatible with that of osf theory .
osf theory is held when the screening length @xmath176 is much smaller than the bare persistence length  @xcite .
therefore , aa theory may be true in the region @xmath177 .
as shown in the figure , @xmath157 for vanishing @xmath8 is well described by the right branch of the curve of eq .
( [ le_ariel ] ) , which seems to support this conjecture .
this topic deserves a further analysis in the future .
experiments and simulations have shown that like - charged macroions can attract each other in a solution  @xcite .
it is hence relevant to know under which conditions the like - charge attraction takes place between polyelectrolytes . to answer this question , we investigate `` potential of mean force '' in this section .
potential of mean force describes the effective interaction between two particles ( molecules ) in a medium , which includes two contributions : the direct interaction between two particles ( molecules ) and the indirect interaction via other particles ( molecules ) and the medium .
this quantity can be calculated by simulations : @xmath178 where @xmath179 is the radial distribution function between particles of species @xmath95 and @xmath180 .
@xmath179 is computed by @xmath181 where @xmath182 and @xmath183 are the number of @xmath95-particles and that of @xmath180-particles , respectively , @xmath184 is the volume of the simulation box , and @xmath185 with @xmath186 denoting the number of @xmath180-particles lying in a spherical shell , centered at the @xmath95-particle @xmath16 , of radius @xmath3 and shell thickness @xmath187 .    because our system is composed of only four chains , it is not easy to reduce statistical error if we directly calculate the potential of mean force between chains .
we , therefore , calculated the potential of mean force between monomers on different chains , which also provides the information about the chain interaction .
we denote this quantity @xmath188 and the results are presented in fig .
[ pmf11p4n48txxx ] which contains six plots , from top to bottom and from left to right , corresponding , in orders , to the cases @xmath82 , 0.5 , 1.0 , 2.0 , 3.0 and 4.0 .
@xmath188 for the referenced neutral polymers is also depicted in the figure .     for @xmath82 , 0.5 , 1.0 , 2.0 , 3.0 and 4.0 , shown , in orders , from top to bottom and from left to right , in six plots .
each plot includes several @xmath188 at different @xmath48 .
the value of @xmath48 is indicated in the plot or near the corresponding curve .
@xmath188 for the neutral polymer system is depicted as reference and marked by `` neutral '' .
the unit of ordinate axis is @xmath21.,title="fig:",scaledwidth=40.0% ]   for @xmath82 , 0.5 , 1.0 , 2.0 , 3.0 and 4.0 , shown , in orders , from top to bottom and from left to right , in six plots .
each plot includes several @xmath188 at different @xmath48 .
the value of @xmath48 is indicated in the plot or near the corresponding curve .
@xmath188 for the neutral polymer system is depicted as reference and marked by `` neutral '' .
the unit of ordinate axis is @xmath21.,title="fig:",scaledwidth=40.0% ]   for @xmath82 , 0.5 , 1.0 , 2.0 , 3.0 and 4.0 , shown , in orders , from top to bottom and from left to right , in six plots .
each plot includes several @xmath188 at different @xmath48 .
the value of @xmath48 is indicated in the plot or near the corresponding curve .
@xmath188 for the neutral polymer system is depicted as reference and marked by `` neutral '' .
the unit of ordinate axis is @xmath21.,title="fig:",scaledwidth=40.0% ]   for @xmath82 , 0.5 , 1.0 , 2.0 , 3.0 and 4.0 , shown , in orders , from top to bottom and from left to right , in six plots .
each plot includes several @xmath188 at different @xmath48 .
the value of @xmath48 is indicated in the plot or near the corresponding curve .
@xmath188 for the neutral polymer system is depicted as reference and marked by `` neutral '' .
the unit of ordinate axis is @xmath21.,title="fig:",scaledwidth=40.0% ]   for @xmath82 , 0.5 , 1.0 , 2.0 , 3.0 and 4.0 , shown , in orders , from top to bottom and from left to right , in six plots .
each plot includes several @xmath188 at different @xmath48 .
the value of @xmath48 is indicated in the plot or near the corresponding curve .
@xmath188 for the neutral polymer system is depicted as reference and marked by `` neutral '' .
the unit of ordinate axis is @xmath21.,title="fig:",scaledwidth=40.0% ]   for @xmath82 , 0.5 , 1.0 , 2.0 , 3.0 and 4.0 , shown , in orders , from top to bottom and from left to right , in six plots .
each plot includes several @xmath188 at different @xmath48 .
the value of @xmath48 is indicated in the plot or near the corresponding curve .
@xmath188 for the neutral polymer system is depicted as reference and marked by `` neutral '' .
the unit of ordinate axis is @xmath21.,title="fig:",scaledwidth=40.0% ]    in a salt - free solution , @xmath188 displays a repulsive interaction , and is stronger than that for the neutral polymers owing to the electrostatic repulsion between monomers .
this can be seen by comparing the range of repulsion estimated by equaling @xmath188 to the thermal energy @xmath21 .
the range of repulsion is at @xmath189 for the polyelectrolytes and at @xmath190 for the neutral polymers .
upon addition of salt , @xmath188 responses in two distinct ways , which depends on the ion size .
the first one is applied for small or large tetravalent counterions , such as @xmath82 , 0.5 , or 4.0 . in these cases , @xmath188 is purely repulsive in the course of addition of salt .
hence , the monomers on different chains repel between themselves and the system exhibits no like - charge attraction . the second way is happened for an intermediate size such as @xmath75 , 2.0 , or 3.0 . at this moment
, @xmath188 displays an attractive well in a mid - salt region .
therefore , it appears like - charge attraction between chains . in the first case
, the range of repulsion of @xmath188 decreases to a minimum value nearly at @xmath126 .
further increase of @xmath48 either leaves @xmath188 unchanged ( for vanishing @xmath8 and for @xmath191 ) or increases the range of repulsion of @xmath188 ( for @xmath146 ) .
we will call @xmath188 at @xmath126 `` the curve at equivalence '' , which separates two behaviors of evolution of @xmath188 as @xmath48 is increased or decreased from @xmath57 .
we noticed that @xmath188 tends toward the referenced curve of the neutral polymers at high salt concentrations except for the deviation happened at small @xmath3 due to the excluded volume of large ions .
for vanishing @xmath8 , the curve at equivalence is very close to the reference curve and @xmath188 is roughly unchanged while @xmath192 . for @xmath146 ,
the curve at equivalence lies below the referenced curve ; @xmath188 is bounded between the two curves and increases with @xmath48 . and in the second case
, @xmath188 displays a potential well located near @xmath193 .
it suggests an attraction mediated by a tetravalent counterion .
the depth of the well can be more profound than @xmath21 .
therefore , a stable bound state can be formed between monomers on different chains and it is a prerequisite to occur salt - induced phase separation .
the maximum depth of the well is happened while @xmath48 is @xmath57 . in general
, we found that @xmath188 is bounded between the salt - free curve and the curve at equivalence while @xmath194 , and between the latter curve and the referenced neutral curve while @xmath192 .
we have shown that the occurrence of the like - charge attraction depends on both salt concentration and ion size . in fig .
[ pmf11p4n48_boundary.eps ] we present , in the variable space , the region where @xmath188 shows an attractive well .    , @xmath8)-variable space . in region
( i ) , @xmath188 shows an attractive well and in region ( ii ) , it is purely repulsive.,scaledwidth=40.0% ]    in the figure , the data points ( symbolized by star ) denote the salt concentrations for a given @xmath8 , at which a purely repulsive @xmath188 turns to show an attractive well and _ vice versa _ , and the closed solid curve is sketched to help readers to separate the two regions .
it clearly shows that an effective attraction appears only when @xmath8 is neither too small nor too big and salt concentration is intermediate around the equivalence concentration .
we mention that the occurrence of an attractive region in the potential of mean force is not sufficient to have a phase separation of the system .
phase transitions in charged colloid systems have been largely discussed in literatures but there have been relatively few studies of phase transitions in charged chain systems ( see refs .
@xcite and @xcite and references therein ) .    fig .
[ snapshots ] shows snapshots of the simulations for @xmath146 , 1.0 and 4.0 at two salt concentrations @xmath126 and @xmath195 .    ,
( b ) @xmath75 , and ( c ) @xmath191 , at @xmath126 ( left ) and @xmath195 ( right ) . in the pictures , the bead - spring chains are represented in gray color ( red color online ) and the monovalent counterions in light one ( white color online ) .
the tetravalent counterions are represented in dark - colored ( blue - colored online ) particles whereas the coions are not plotted for the clarity of the pictures . in the snapshots
, a chain may appear more than one fragment because of periodic boundary condition.,title="fig:",scaledwidth=40.0% ] , ( b ) @xmath75 , and ( c ) @xmath191 , at @xmath126 ( left ) and @xmath195 ( right ) .
in the pictures , the bead - spring chains are represented in gray color ( red color online ) and the monovalent counterions in light one ( white color online ) .
the tetravalent counterions are represented in dark - colored ( blue - colored online ) particles whereas the coions are not plotted for the clarity of the pictures . in the snapshots
, a chain may appear more than one fragment because of periodic boundary condition.,title="fig:",scaledwidth=40.0% ]    , ( b ) @xmath75 , and ( c ) @xmath191 , at @xmath126 ( left ) and @xmath195 ( right ) .
in the pictures , the bead - spring chains are represented in gray color ( red color online ) and the monovalent counterions in light one ( white color online ) .
the tetravalent counterions are represented in dark - colored ( blue - colored online ) particles whereas the coions are not plotted for the clarity of the pictures . in the snapshots
, a chain may appear more than one fragment because of periodic boundary condition.,title="fig:",scaledwidth=40.0% ] , ( b ) @xmath75 , and ( c ) @xmath191 , at @xmath126 ( left ) and @xmath195 ( right ) .
in the pictures , the bead - spring chains are represented in gray color ( red color online ) and the monovalent counterions in light one ( white color online ) .
the tetravalent counterions are represented in dark - colored ( blue - colored online ) particles whereas the coions are not plotted for the clarity of the pictures . in the snapshots
, a chain may appear more than one fragment because of periodic boundary condition.,title="fig:",scaledwidth=40.0% ]    , ( b ) @xmath75 , and ( c ) @xmath191 , at @xmath126 ( left ) and @xmath195 ( right ) .
in the pictures , the bead - spring chains are represented in gray color ( red color online ) and the monovalent counterions in light one ( white color online ) .
the tetravalent counterions are represented in dark - colored ( blue - colored online ) particles whereas the coions are not plotted for the clarity of the pictures . in the snapshots
, a chain may appear more than one fragment because of periodic boundary condition.,title="fig:",scaledwidth=40.0% ] , ( b ) @xmath75 , and ( c ) @xmath191 , at @xmath126 ( left ) and @xmath195 ( right ) .
in the pictures , the bead - spring chains are represented in gray color ( red color online ) and the monovalent counterions in light one ( white color online ) .
the tetravalent counterions are represented in dark - colored ( blue - colored online ) particles whereas the coions are not plotted for the clarity of the pictures . in the snapshots
, a chain may appear more than one fragment because of periodic boundary condition.,title="fig:",scaledwidth=40.0% ]    at @xmath126 , we saw that the polyelectrolytes favor to be apart from each other for both of the cases @xmath146 and @xmath191 .
however , delicate difference was witnessed that the previous case displays a single - chain collapse but the latter one does not . for @xmath75 , the chains favor to aggregate together , in accordance with the results of fig .
[ pmf11p4n48txxx ] which show the occurrence of an attraction between the negatively - charged polymers .
this attraction is mediated by counterions , mainly the tetravalent counterions , since almost all the tetravalent counterions are condensed on the polyelectrolytes .
we point out that the condensation of tetravalent counterion is not a sufficient condition to occur chain aggregation .
it can be seen from the case @xmath146 in which most of the tetravalent counterions are condensed but the chains does not aggregate . at the high salt concentration @xmath195 ,
the attraction between chains is reduced , in comparison with the cases at @xmath126 .
the chains segregate or show tendency of segregation from each other . in the simulations
, we observed that while @xmath192 , only part of the tetravalent counterions can condense onto the chains .      in this section ,
we investigate the charge distribution around a polyelectrolyte .
it enables to understand how ions distribute around a chain and provides useful information for theorists to develop models .
we define a worm - like tube of radius @xmath3 around a chain to be a union of @xmath44 spheres of radius @xmath3 centered at each monomer on the chain .
we calculated the averaged total charge inside a worm - like tube , including the charge of the ions and the one of the monomers .
this quantity is called `` the integrated charge distribution '' and is denoted by @xmath196 .
we show in fig .
[ iondist ] @xmath196 at various @xmath48 for @xmath8 ranged between 0.0 and 4.0 .     at different @xmath48 for @xmath8=0.0 , 0.5 , 1.0 , 2.0 , 3.0 and 4.0 , shown , in orders , from top to bottom and from left to right , in six figures .
the line arts used for different @xmath48 are indicated at the right - bottom side of each figure .
, title="fig:",scaledwidth=40.0% ]   at different @xmath48 for @xmath8=0.0 , 0.5 , 1.0 , 2.0 , 3.0 and 4.0 , shown , in orders , from top to bottom and from left to right , in six figures .
the line arts used for different @xmath48 are indicated at the right - bottom side of each figure .
, title="fig:",scaledwidth=40.0% ]   at different @xmath48 for @xmath8=0.0 , 0.5 , 1.0 , 2.0 , 3.0 and 4.0 , shown , in orders , from top to bottom and from left to right , in six figures .
the line arts used for different @xmath48 are indicated at the right - bottom side of each figure .
, title="fig:",scaledwidth=40.0% ]   at different @xmath48 for @xmath8=0.0 , 0.5 , 1.0 , 2.0 , 3.0 and 4.0 , shown , in orders , from top to bottom and from left to right , in six figures .
the line arts used for different @xmath48 are indicated at the right - bottom side of each figure .
, title="fig:",scaledwidth=40.0% ]   at different @xmath48 for @xmath8=0.0 , 0.5 , 1.0 , 2.0 , 3.0 and 4.0 , shown , in orders , from top to bottom and from left to right , in six figures .
the line arts used for different @xmath48 are indicated at the right - bottom side of each figure .
, title="fig:",scaledwidth=40.0% ]   at different @xmath48 for @xmath8=0.0 , 0.5 , 1.0 , 2.0 , 3.0 and 4.0 , shown , in orders , from top to bottom and from left to right , in six figures .
the line arts used for different @xmath48 are indicated at the right - bottom side of each figure .
, title="fig:",scaledwidth=40.0% ]    focus firstly on the case with vanishing @xmath8 . in the region
@xmath197 , @xmath196 attains a value equal to the bare charge of a chain , @xmath198 .
a stair - like behavior then appears for @xmath196 in the region @xmath199 while @xmath200 . outside the region , @xmath196 monotonically increases and tends toward zero at large @xmath3 due to electroneutrality .
since @xmath201 is the depletion zone for the particles other than tetravalent counterions , the behavior of @xmath196 in this region is completely determined by the condensation of tetravalent counterions .
the stair - like behavior indicates a _ perfect _ condensation on the surface of the chain , owing to the vanishing size of the tetravalent counterions .
for instance , at @xmath202 , we have 24 tetravalent counterions in the simulation box and , therefore , 6 condensed tetravalent counterions on a chain , in average .
a perfect condensation will lead to a value @xmath203 for @xmath196 and this is exactly what we have observed in the figure . at the equivalence concentration ,
@xmath196 is zero while @xmath204 .
if @xmath48 is further increased ( @xmath192 ) , @xmath196 reveals a positive peak in the region @xmath205 , which indicates an overcharging , and , then , rapidly decays to zero as @xmath206 .
the hight of the peak suggests that only a portion of the tetravalent counterions are condensed on the chain . for the case
@xmath146 , the region in which only the tetravalent counterions can intervene is narrowed to @xmath207 . in the region ,
the step of the stair - like @xmath196 becomes tilted while @xmath194 . if @xmath192 , @xmath196 displays a positive peak , which is more pronounced than that for vanishing @xmath8 .
for the case @xmath75 , an oscillatory behavior occurs at high @xmath48 .
similar phenomena have been observed in other systems such as charged plate  @xcite , colloid system  @xcite , and rigid polyelectrolyte  @xcite .
this oscillatory behavior indicates a multi - layer organization of ions : in the close neighborhood of a polyelectrolyte , tetravalent counterions condense and form a positively - charged layer which overcharges the chain ; exterior to this layer , anionic particles intervene and form a negatively - charged layer which overcompensates the charge inside and , hence , the integrated charge becomes negative , and so forth .
notice that the depletion zone around a chain for tetravalent counterions is @xmath208 .
it becomes larger than that for the other ions while @xmath209 . in these cases
, @xmath196 decreases with increasing @xmath48 in the depletion zone , but shows completely opposed behavior just outside the zone apparently due to the condensation of tetravalent counterions . at very high @xmath48 ( for instance , at @xmath210 for @xmath211 ) , @xmath196 in the zone
can even attain a value smaller than the chain bare charge . in this case ,
negative particles are the major composition of the first ion layer around a chain , which is counterintuitive .
the following three observations are worth to be noticed .
( 1 ) the salt concentration at which @xmath196 becomes ( roughly ) zero in the region @xmath212 decreases with increasing @xmath8 . ( 2 ) the oscillatory behavior appears more markedly for large @xmath8 than for small one . ( 3 ) at a given @xmath48 , the larger the @xmath8 , the higher the maximum peak of @xmath196 will be .    it is tempting to consider that a polyelectrolyte and the ions condensed on the chain form a complex object . under some conditions ,
the effective charge of the complex object reverses its sign .
this phenomenon is called charge inversion  @xcite .
_ intended to link the phenomenon of charge inversion with the reentrant condensation  @xcite .
they proposed that zero effective charge and charge inversion of chains are , respectively , the cause for the condensation and for the redissolution of polyelectrolytes . their theory necessitates further verification either by experiments or by simulations . in the rest of this section ,
we try to calculate the effective charge of a polyelectrolyte , in hope of being able to clarify the relationship between the charge inversion and the reentrant condensation .
we simply define the complex to be the set of particles lying within a worm - like tube of radius @xmath213 around the chain .
this definition has been utilized in many places , for example , in refs .
@xcite and @xcite .
usually , @xmath213 is chosen to be the distance at which the electrostatic attraction between a monovalent counterion and a monomer is equal to the kinetic energy @xmath214 ; in consequence , the counterion inside the tube of radius @xmath213 can not escape from the chain .
we , thus , consider the total charge @xmath215 inside the tube as the effective charge of a chain . in our study
, @xmath213 is equal to 2 .
the results are presented in fig .
[ iondist_p4n48_rc200.eps](a ) .     as a function of @xmath48 with @xmath216 , and ( b ) @xmath217 as a function of @xmath48 , for different @xmath8 .
the symbol used to denote each @xmath8 is indicated in the figures . , title="fig:",scaledwidth=40.0% ] + ( a ) +   as a function of @xmath48 with @xmath216 , and ( b ) @xmath217 as a function of @xmath48 , for different @xmath8 . the symbol used to denote each @xmath8
is indicated in the figures . , title="fig:",scaledwidth=40.0% ] + ( b )    two behaviors were observed : ( 1 ) for @xmath78 and 4.0 , @xmath215 decreases with increasing @xmath48 , ( 2 ) for @xmath218 , @xmath215 increases basically with @xmath48 .
notice that @xmath215 for @xmath160 or 2.0 attains a large positive value in the mid - salt region .
hence , @xmath188 should show strong repulsion in this region , according to the theory of nguyen _
et al_. but , in fact , we observed a strong attractive well for @xmath188 and the chains were observed to aggregate together .
therefore , the results do not support the theory .
we notice that it is too simple to regard @xmath215 as the effective chain charge , because ion valence is not considered in the definition .
for example , the suitable @xmath213 for tetravalent counterions should be 8.0 , based upon the same argument to choice @xmath213 for monovalent counterions .
in addition , @xmath216 becomes too small while @xmath219 , since the depletion zone around a chain for the tetravalent counterions is bigger than the tube .
therefore , we redefine @xmath220 as the novel condensation region for tetravalent counterions , whereas keeping @xmath216 for monovalent counterions and for monomers on the other chains . considering that the condensation of coions is mainly mediated by tetravalent counterions ,
we choose @xmath220 for coions .
the charges inside the different tubs are computed and their sum @xmath217 is then regarded as an alternative definition for the effective chain charge .
the results are presented in fig .
[ iondist_p4n48_rc200.eps](b ) .
we remark that @xmath217 looks quite differently from @xmath215 .
it shows a positive peak in the mid - salt region for @xmath75 , 1.4 , and 2.0 but retains negative value for the other @xmath8 .
the obtained results suggest a disconnection between neutralization of chain charge and condensation of polyelectrolytes .
our simulating system in the high - salt region reveals a picture more close to the one described by solis  @xcite : although _ locally _ overcompensated by condensed counterions , the net charge of a redissolved chain can be positive or negative due to the association of coions .
we must point out that the effective chain charge interpreted by @xmath215 or by @xmath217 in this section can be only served as a reference because both @xmath215 and @xmath217 are sensitive to the choice of @xmath213 . a rigorous way to study it , which can circumvent the ambiguity , is to study the electrophoretic mobility of polyelectrolytes in salt solutions . this part of research is currently undergoing .
we have performed molecular dynamics simulations in canonical ensemble to investigate the properties of highly - charged polyelectrolytes in tetravalent salt solutions , including chain morphology , single - chain structure factor , swelling exponent , persistence length , like - charge attraction between chains and integrated charge distribution around a polyelectrolyte .
we have explored a wide range of salt concentration and studied the effect of size of tetravalent counterions on the above properties .
we have found that a polyelectrolyte exhibits , in turns , two conformational transitions upon addition of salt : a collapsed transition and a reexpanding transition , manifested by a v - shaped curve of the square of radius of gyration @xmath47 against salt concentration @xmath48 .
the systems at different monomer concentration @xmath53 show similar curves .
the curves overlap each other in the chain - reexpansion region . in the chain - collapsed region , we have also observed that @xmath60 against @xmath59 lie over each other for different @xmath53 , where @xmath62 and @xmath221 weakly depends on @xmath53 .
these findings are useful in understanding the phase boundaries of the condensation of polyelectrolytes : on the upper boundary , @xmath48 is roughly independent of @xmath53 and on the lower one , @xmath48 depends linearly on @xmath53 .
the study of chain morphology shows that for the cases with @xmath8 ( the size of the tetravalent counterions ) around 1.0 , the chains reveal compact and prolate structures in the vicinity of the equivalence concentration @xmath57 . in a low - salt region and in a high - salt region , the chains show extended structures . at a very high salt concentration ,
the chains behave similarly to neutral polymers . for the cases with large @xmath8 , the chains are elongated and do not display a collapsed - transition in the mid - salt region . for the case with vanishing @xmath8 ,
the behavior of the chains is analogous to that of neutral polymers as @xmath222 .
noticeably , our results show that @xmath75 is the optimal condition to pack polyelectrolytes into the smallest volume .
we have shown that the single - chain structure factor @xmath122 of the multi - chain system in salt solutions follows the guinier function @xmath129 while @xmath223 .
@xmath122 then exhibits a power - law - like behavior @xmath224 in the region @xmath134 , suggesting that the scaling law , @xmath225 , is held for polyelectrolytes in salt solutions .
the swelling exponent @xmath139 attains a minimum value near @xmath57 and behaves in analogy of the morphological quantities .
it shows once more that the chains are mostly collapsed while @xmath8 is compatible with 1.0 .
our simulations incorporated explicitly the salt ions and , therefore , were able to study the negative regime of the electrostatic persistence length @xmath157 .
this regime appears simultaneously with the mechanical instability of the polyelectrolytes induced by salt .
we have examined the theory of ariel and andelman  @xcite and found that it fails to predict @xmath157 in the _ so - claimed _ valid region @xmath226 but describes @xmath157 well in the unexpected region @xmath227 . since the osf theory is valid in the latter region , we conjectured that the theory of ariel and andelman is applicable in the latter region .
this part deserves a detail investigation in the future .
we have studied like - charge attraction between polyelectrolytes by calculating the potential of mean force @xmath188 between monomers on different chains .
the results show that it depends strongly on salt concentration and ion size .
an attractive region appears in @xmath188 only when the salt concentration is intermediate and the tetravalent counterions possess a size compatible with the monomers .
the most profound attractive well is happened while @xmath75 and @xmath228 .
in the final part , the integrated charge distribution @xmath196 around a polyelectrolyte was investigated .
we have found that the chain charge can be locally overcompensated by condensed ions .
@xmath196 displays an oscillatory behavior at high salt concentrations while @xmath80 .
the oscillatory behavior results from a multi - layer organization of ions around the chain , which is in agreement with the picture predicted by solis  @xcite . we have used two different definitions to calculate the effective chain charge .
the results suggest that charge inversion is not a necessary condition to occur reentrant condensation .
recently wen and tang have provided evidence against the connection between charge inversion and resolubilization of polyelectrolyte bundles  @xcite .
a more rigorous study by electrophoresis performed by molecular dynamics simulations is undergoing .
the author acknowledges c .- l .
lee , y .- f .
wei , and n. pardillos for their comments .
this material is based upon work supported by the national science council , the republic of china , under the contract no .
nsc 94 - 2112-m-007 - 023 .
most of the simulations were running at the national center for high - performance computing .
the author expresses his gratitude to the members and the staffs of the council and the center .