electrostatic interactions of charged macromolecules and colloids are often governed by small oppositely charged ions ( counterions ) that maintain global electroneutrality .
these counterion - mediated interactions play a fundamental role in classical charged ( coulomb ) fluids that are abundant in biological and soft matter context @xcite and include many charged macromolecules ( such as nucleic acids dna and rna , actin filaments , microtubules and globular proteins ) , affecting their functional , structural and dynamical behavior . in spite of the importance of electrostatic interactions , there is no general method that would allow for an accurate prediction of electrostatic effects in all regions of the parameter space , defined by the surface charge density of macroions , charge valency of counterions , dielectric mismatches between the often hydrophobic core of the macromolecule and the surrounding aqueous medium , etc .
often the electrostatic interactions are treated on the poisson - boltzmann ( pb ) level leading to effective interactions which turn out to be always repulsive between like - charged macromolecules .
conceptually , the pb approach corresponds to a mean - field treatment of electrostatic interactions and is asymptotically valid only for sufficiently large separations between macromolecules , low enough surface charge densities and low counterion valency @xcite .
it characterizes a situation where electrostatic fluctuations and correlations due to the counterions are negligible .
there are other regions in the parameter space of charged macromolecules where one expects the mean - field framework to break down leading to the emergence of a completely different non - pb - type physics .
a notorious example is the phenomenon of like - charge attraction , which emerges between highly charged macroions or in the presence of high valency counterions and has been at the focus of both experimental @xcite and theoretical investigations in recent years ( see refs .
@xcite for an extensive reference list ) .
it appears to us that among the most important recent advances in this field has been the systematization of non - pb effects based on the notions of _ weak coupling _ ( wc ) and _ strong coupling _ ( sc ) approximations .
these terms refer to the strength of electrostatic coupling in the system and may be understood conceptually in terms of the electrostatic interactions of mobile counterions with fixed external charges ( macroions ) in the system when compared with direct electrostatic interactions between the counterions themselves .
this latter contribution may be characterized in terms of the _ bjerrum length _ , @xmath0 which corresponds to the separation at which two unit charges , @xmath1 , interact with thermal energy @xmath2 in a medium of dielectric constant @xmath3 ( in water and at room temperature , the value is @xmath4 nm ) . if the charge of the counterions is @xmath5 then the bjerrum length scales as @xmath6 .
the interaction of counterions with macroion charges ( of surface charge density @xmath7 ) can be characterized in terms of the so - called _ gouy - chapman length _ ,
@xmath8 which gives the separation at which the counterion - surface interaction energy equals @xmath2 .
the ratio of these two fundamental length scales introduces a dimensionless parameter @xmath9 which is known as the ( netz - moreira ) electrostatic coupling parameter @xmath10 @xcite and quantifies the strength of electrostatic coupling in the system .
this parameter is closely related to the plasma parameter of ionic systems @xcite and may be written also in terms of the typical lateral spacing , @xmath11 , between counterions in the proximity of a charged surface , i.e. , @xmath12 as set by the local electroneutrality condition @xmath13 .
it then follows that in the weak coupling regime , defined by @xmath14 , the width of the gouy - chapman layer @xmath15 is much larger than the separation between two neighbouring counterions in the counterion layer and thus it behaves basically as a three - dimensional gas .
each counterion in this case interacts with many others and the collective mean - field approach of the poisson - boltzmann type is thoroughly justified . on the other hand , in the strong coupling regime , defined by @xmath16 ,
the mean distance between counterions , @xmath11 , exceeds the layer width and thus the counterion layer behaves essentially as a two - dimensional layer @xcite . in this case
the mean - field approach breaks down as each counterion is isolated laterally in a relatively large correlation hole of size @xmath11 .
in fact , as each counterion can move almost independently from the others along the direction perpendicular to the charged surface , the properties of the system are dominated by single - particle contributions on the leading order , which is in stark contrast with the collective mean - field picture and emerges as a direct consequence of strong electrostatic correlations in the system .
the two dychotomous limits are characterised by a low ( high ) valency of the counterions , a small ( large ) value of the surface charge density and/or large ( small ) medium dielectric constant .
conceptually the study of the sc regime has been pioneered in several recent works @xcite using various analytical methods .
it was shown @xcite that both the wc and the sc limits may be described analytically as two exact _ limiting laws _ from a single general theory for classical coulomb fluids : while the pb theory follows in the limit of @xmath17 , a limiting single - particle sc theory follows in the limit @xmath18 , which thus allow for an exact statistical mechanical treatment of charged systems at two opposed limiting conditions
. the parameter space in between can be analyzed by approximate methods @xcite , being accessible effectively only _ via _ computer simulations @xcite .
exact solutions for the whole range of coupling parameters are available only in one dimension @xcite .
the wc - sc paradigm has been tested extensively @xcite and was found to describe computer simulations quantitatively correctly in the respective regimes of validity , thus providing a unifying view of the behavior of coulomb fluids .
the main facets of these results are retained even when the model system is generalized in order to include more realistic features describing the bathing solution or the nature of the fixed or mobile charges in the system .
though , for instance , multipolar charge distribution of mobile counterions @xcite or statistically disordered distribution of fixed charges @xcite unavoidably introduce novel features in the counterion - mediated electrostatic interaction , the application of the same general philosophy embodied in the weak and strong coupling limits remains sound and valid .
the only case where it needs to be amended in an essential manner is when the bathing solution contains a mixture of univalent as well as polyvalent salts , which incidentally are also the most common experimental conditions . in that case
a more sophisticated mixed weak - strong coupling analysis is in order leading to qualitatively different results @xcite .    though originally formulated in the context of planar macroions , these advances have transpired also in the dna - like models of polyelectrolytes which deal with electrostatic interactions in the context of macroions with cylindrical or indeed helical fixed charge distribution .
indeed this particular variant of counterion - mediated interactions has a venerable history @xcite . during the last two decades several distinct analytical approaches aimed in different directions improved the classical pb results for simple dna - like models and revealed the importance of correlation effects as well as several other factors including the discrete or helical charge distributions , chain flexibility , finite molecular size , and dielectric inhomogeneities @xcite .
guided by these developments we set ourselves the goal of systematically analyzing the interactions between cylindrical macroions mediated by mobile counterions in the presence of dielectric inhomogeneities .
contrary to the case of planar macroions which can be completely characterized by a single dimensionless coupling parameter , cylindrical macroions require in general two dimensionless coupling parameters that consistently describe the range of validity of the strong and weak coupling approximations . the existence of two coupling parameters is due to the simple fact that , if compared to the planar case , the cylinder has a finite radius that needs to enter the fundamental description of the problem .
this other dimensionless parameter brought fourth by cylindrical geometry of the macroion is nothing but the so - called _ manning parameter _ @xcite . for a cylinder of radius @xmath19 and linear charge density @xmath20 ,
it is given by @xmath21 the manning parameter thus represents the dimensionless linear charge density or , on the other hand , also the rescaled radius of the cylindrical charge distribution .
note that the two parameters ( i.e. , @xmath10 and @xmath22 ) in the cylindrical geometry are independent and can be set separately . for double - stranded dna ( @xmath23 ) , the manning parameter and the ( netz - moreira ) electrostatic coupling parameter are given by @xmath24 and @xmath25 ( at room temperature in water ) in terms of the counterion valency @xmath26 .
in this paper we will analyze systematically the interactions between two cylindrical macroions characterized by a fixed surface charge density as well as a finite dielectric jump between the external dielectric background ( corresponding to a continuum solvent ) and the hydrophobic cores of cylindrical macroions .
analytically derived results in both limits of strong and weak coupling will be compared directly with monte - carlo simulation results .
the organization of this paper is as follows : we first introduce our model and study the image - charge effects within the problem of a single charged cylinder with neutralizing counterions in section [ sec : one_cyl ] .
we then focus on the interaction between two identical charged dielectric cylinders within the wc and the sc theory as well as mc simulations in section [ sec : two_cylinder ] .
the results in the presence and absence of dielectric discontinuity at the cylindrical boundaries are analyzed in detail in order to bring up the effects due to the image charges in the two - cylinder geometry .
let us focus first on the problem of a single infinitely long uniformly charged cylindrical macroion of radius @xmath19 .
the charge of the cylinder is assumed to be distributed uniformly on its surface according to the charge distribution function @xmath27 with @xmath20 being the absolute _ linear _ charge density and @xmath28 the radial coordinate from the cylinder axis .
the cylinder is arbitrarily chosen to be negatively charged and its electrical charge is thus neutralized by positively charged mobile counterions of valency @xmath29 , which are present in the region @xmath30 only .    in general ,
the interior of the cylindrical macromolecule can be characterized by a different dielectric constant @xmath31 than the surrounding continuum solvent medium , @xmath3 .
this is certainly the case for dna that has a hydrophobic inner core composed mostly of stacked nitrogen bases that have a vastly different dielectric response from an aqueous solution @xcite .
other charged polymers that usually do not share structural features with dna can nevertheless also posses hydrophobic inner cores with a local dielectric response that differs from the one of the solvent .
the hydrophobic apolar core of the macromolecular backbone would then have a static dielectric constant @xmath32 ( hydrocarbon ) , while the surrounding continuum dielectric medium which is usually water , a polar associated liquid , would be characterized by @xmath33 .
one can thus introduce the dielectric discontinuity parameter @xmath34 which measures the relative dielectric mismatch at the interface of the two materials . in a dielectrically homogeneous system
, we have @xmath35 and no image charge effects are present , while , for instance , in water - hydrocarbon systems , one has @xmath36 , which suggests strong image charge effects ( note that @xmath37 and the largest value for @xmath38 is obtained when one medium is ideally polarizable , i.e. , is an ideal metal ) .
therefore , a single cylinder can be described by three different dimensionless parameters , namely , the electrostatic coupling parameter @xmath10 , the manning parameter @xmath22 and the dielectric discontinuity parameter @xmath38 as defined in eqs .
( [ eq : xicouple ] ) , ( [ eq : xi_m ] ) and ( [ eq : jump ] ) , respectively .
the presence of a dielectric inhomogeneity across the boundary of the cylinder , see fig .
[ fig : model_1cyln ] , influences the electrostatic potential that can be thus described by the green s function connecting two points @xmath39 outside the cylinder as @xmath40 where @xmath41 is the direct standard coulomb kernel in the absence of dielectric inhomogeneities , @xmath42 and @xmath43 is the  image correction " due to the dielectric jump just as in the case of a planar discontinuity .
unfortunately in cylindrical geometry the concept of kelvin image charges , that plays such a fundamental role for planar dielectric discontinuities , is a bit diluted since the image correction can not be in general formulated in a way that would entail a summation of properly positioned point image charges .
we nevertheless consistently refer to the effects of dielectric inhomogeneities in this system as dielectric image effects .
if explicitly stated we believe this inconsistency in nomenclature can not be a source of confusion .    instead of using the concept of point image charges one must in fact
explicitly solve the poisson equation in cylindrical geometry specified by cylindrical coordinates @xmath44 with appropriate boundary conditions at the surface of the cylinder , which requires the normal component of the electric field to fulfill the following relation @xmath45 where @xmath46 is the cylinder surface charge density . solving the poisson equation with the proper eigenfunction expansion @xcite
one can express the image part of the green s function as @xmath47 with the following definition @xmath48 where @xmath49 and @xmath50 are the angle and height differences between the position vectors @xmath51 and @xmath52 , respectively .
@xmath53 and @xmath54 are the modified bessel functions of the first and the second kind that enter the eigenfunction expansion @xcite .    having derived the appropriate green s function for the cylindrical geometry , we now investigate the behavior of counterions in proximity to a charged dielectric cylinder in the two limits specified by the weak and strong coupling approximations that can be treated analytically .
the behavior at intermediate coupling strengths will be analyzed by using extensive mc simulations .
as noted before , the wc regime is characterized by a small coupling parameter @xmath14 , which is adequate for low valency counterions , small surface charge , high temperature or high solvent dielectric constant . in the strict limit of @xmath17 , ignoring the usually small fluctuations around the mean field configuration , the system is described exactly by the mean - field pb theory @xcite .
formally the pb equation corresponds to the saddle - point condition imposed on the action field - functional in this limit @xcite .
for the present system , the pb equation governing the mean electrostatic potential takes the form @xmath55 where the right - hand side is obviously the number density of the counterions , i.e. @xmath56 note that due to axial symmetry and translational invariance in the @xmath57 direction , the problem is reduced to a one - dimensional formulation with @xmath58 and the solution of the pb equation depends only on the radial coordinate @xmath28 .
the analytical solution of the pb equation in this case is well known and was first discussed in refs .
@xcite within the so - called _ cell model _ where the cylindrical macroion is confined ( coaxially ) in an outer cylindrical boundary of radius @xmath59 .
the cell model guarantees a finite solution for the counterion density despite the extremely slowly varying ( logarithmic ) electrostatic potential in two dimensions .
the solution of the pb equation for the counterion density takes different forms depending on whether the manning parameter is larger or smaller than a threshold value given by @xmath60 $ ] @xcite .
here we are interested mainly in the situations with sufficiently large manning parameters @xmath61 , where the normalized counterion density can be written as @xmath62 ,   \label{npb}\ ] ] with @xmath63 determined from the transcendental equation @xmath64}.\ ] ] note that the density profile is renormalized such that we have @xmath65 the above density distribution is a monotonically decaying function with very slowly convergent asymptotics @xcite .
it is also independent of the dielectric jump parameter , @xmath38 , due to the axial symmetry which implies that the electrostatic field vanishes inside the cylinder .
the absence of dielectric discontinuity effects is specific to the mean - field limit where no fluctuations are taken into account and can be derived also in the case of planar slabs @xcite .
when the manning parameter is decreased , the system exhibits a continuous counterion condensation transition @xcite at a critical manning parameter @xmath66 .
the nature of this transition has been analyzed throughly by means of analytical and numerical methods @xcite .
it was in particular shown that the behavior close to the transition point is described by the mean - field theory and fluctuation and correlation effects ( that will be important eventually at sufficiently large manning parameters ) play no role at the transition point . here
we shall not consider the behavior of counterions close to the transition point but do note that since the image - charge effects are absent in the mean - field limit , such effects are expected to have no influence on the counterion condensation transition itself .    in
what follows , we shall also take into account the hard - core repulsion between counterions and the cylinder by assuming that the counterions have a finite radius @xmath67 . within the pb theory ,
this can not be done explicitly unless full corrections due to excluded - volume repulsions between counterions are taken into account @xcite . here
we shall consider a simplified version of the ion size effects by taking an effective hard - core radius for the cylinder , i.e. , by setting @xmath68 .
though approximate on the wc level , this procedure turns out to be exact in the sc limit @xcite .
( a )    0.25 cm     ( b )    0.25 cm     ( c )    0.25 cm      in the regime when @xmath69 , the mean - field approximation breaks down , and a different kind of approach is needed .
we shall focus on the limit of @xmath18 , where the system can be treated via a systematic strong coupling virial expansion @xcite leading to an exact analytical theory in the leading order in @xmath10 . by construction
the sc theory is a single particle theory and contains contributions only from single - particle interactions .
we do not elaborate on the derivation of the sc theory here as further details can be found in previous publications @xcite .    in this section
, we shall be interested in the counterion density profile around the cylinder in the sc limit .
the general form of the counterion distribution within the sc theory turns out to be given in terms of a boltzmann factor containing the interaction energy of an individual counterion with a fixed charged macroion @xcite . in a dielectrically inhomogeneous system , such as a charged cylindrical core considered here ,
the sc counterion density may be written as @xcite @xmath70 ,   \label{eq : density_ci}\ ] ] where @xmath71 is the dielectric image _ self - energy _ of a single counterion , i.e. , the electrostatic energy of a charged point particle in the vicinity of a neutral dielectric cylinder .
this energy of course represents the contribution from the interaction of a counterion with its image charge and thus depends crucially on the value of the dielectric discontinuity at the macroion surface .
furthermore , @xmath72 where @xmath73 is the bare electrostatic interaction energy of a single counterion located at a given position @xmath74 with the cylinder charge ( up to an irrelevant additive constant ) .
it goes to zero when the charge density on the macroion surface goes to zero .
the image - dependent part of the interaction energy between a point charge and the surface charge distribution on the dielectric cylinder , @xmath75 , can be shown to be nil due to the axial symmetry , i.e. @xmath76 if the surface charge distribution on the macroion is non - uniform , such as in the case of charged helical stripes @xcite , the image dependent part of the interaction energy is non - zero and should be considered explicitly .    note that physically the dielectric image - dependent part of the interaction energy corresponds to the interaction between the dielectrically
induced image charge of the surface charge on the cylinder with the counterion as well as the image charge of the counterion with the surface charge on the cylinder , as follows from the general definition of the green s function ( [ eq : u_0u_im ] ) .
using the above equations , we obtain the sc counterion number density profile ( @xmath30 ) as @xmath77 ,   \label{nsc}\ ] ] where we have introduced the dimensionless integral @xmath78 and the numerical prefactor @xmath79 is determined from the normalization condition , eq .
( [ normrho ] ) .
we can extract two limiting behaviors for the function @xmath80 , viz . @xmath81
the limiting form of @xmath80 for @xmath82 implies that at very small separations between a single counterion and the cylindrical surface , the dielectric self - energy has the same form as in the case of a counterion next to a planar wall , where the polarization effects can be described by an image charge inside the wall @xcite .
the sc density for @xmath35 reduces to the well - known form for a homogeneous system @xmath83 @xcite .
if the dielectric core has a smaller dielectric constant than the medium ( @xmath84 ) , as one typically encounters in the case of biological macromolecules , the image charges have the same sign as counterions and thus lead to a depletion of counterions in the vicinity of the cylinder .
this behavior can be seen from the sc density profiles plotted in fig .
[ fig : dens]a as the dielectric mismatch parameter @xmath38 or the coupling parameter @xmath10 is increased .
they show a reduced density close to the cylinder ( depletion zone ) and a peak some distance away from the cylinder surface .
one can estimate the location of the density peak by setting the density ( [ nsc ] ) derivative to zero , @xmath85 . using the approximation ( [ eqi ] ) for @xmath86 , where @xmath87 is the peak distance from the cylinder surface
, we obtain to the leading order for small @xmath87 @xmath88 which shows good agreement with the peak location values of the sc density profile as seen in fig .
[ fig : dens]a .
as noted before , effects of a finite counterion radius @xmath67 within the sc approximation can be taken into account exactly ( due to the single particle nature of the theory ) by increasing the hard - core radius of the cylinder to @xmath89 .
the counterion - counterion excluded - volume repulsion effects are absent within the leading - order sc theory and enter only in the subleading terms .      we have performed extensive mc simulations in order to verify the validities of the weak- and strong - coupling approaches .
the detailed description of mc simulation is given in appendix [ app : mc ] .    in fig .
[ fig : dens]b we show the radial density profile of counterions around the cylinder without dielectric image effect , @xmath35 . apparently in this case the simulation data are always bracketed by the two analytical results obtained from the wc and sc theories , eqs .
( [ npb ] ) and ( [ nsc ] ) , respectively .
the wc theory is found to be valid at sufficiently small coupling parameter or sufficiently large radial distances from the cylinder .
the sc theory is valid in the opposite regime , i.e. , for sufficiently large coupling parameters or sufficiently small radial distances .
in fact , in agreement with the general trend obtained for planar surfaces @xcite , the validity regime of the sc ( wc ) theory expands to large ( smaller ) separations as the coupling becomes larger ( smaller ) . in the case of a planar charged wall
, the validity regime of the sc theory can be estimated systematically as @xmath90 , where @xmath57 is the distance from the charged plane @xcite . in the case of a charged cylinder ,
a similar criterion is yet to be obtained .    in fig .
[ fig : dens]c , we show the results for the case where the dielectric constant of the cylinder is different from that of the medium such that @xmath91 .
as already discussed the dielectric discontinuities have no effects in the wc limit as confirmed also by the agreement obtained between the simulation data and the wc ( mean - field ) result for small @xmath10 .
by contrast , the image effects become quite significant at large coupling parameters as the density profile deviates qualitatively from those obtained with @xmath35 ( fig .  [
fig : dens]b ) .
the simulation data again show good agreement with the sc approximation at high enough couplings and small enough radial distances and in particular exhibit a depletion zone and a peak at small distances from the cylinder surface as predicted within the sc theory .
note that in a dielectrically inhomogeneous system ( @xmath92 ) the sc theory has an explicit dependence on the coupling parameter @xmath10 because of the self - interaction of counterions with their image charges .
we now turn our attention to the problem of the interaction between two like - charged cylinders , which is commonly used as a model for the interaction between rigid polyelectrolyte chains @xcite .
we consider two identical parallel cylinders of radius @xmath19 oriented along the @xmath57 axis at an interaxial separation of @xmath93 ( see fig .
[ fig : model_2cyln ] ) .
the electric charge is assumed to be uniformly distributed on the surface for both cylinders with equal linear charge density @xmath20 .
the cylinders are assumed to be infinitely long of length @xmath94 and are standardly confined within a confining square box of lateral size @xmath95 ( we shall assume @xmath96 in section [ subsec : wc_2cyln ] and @xmath97 in sections [ subsec : sc_2cyln]-[subsec : sc_2cyln_im ] ) , but we should emphasize that in the regime of manning parameters considered here , the lateral box size plays no significant role @xcite .    ( a )    ( b )    we follow the same approach as in the case of a single cylinder and thus first analyze the problem in the weak - coupling limit and then in the strong - coupling limit .
we then compare the results from these two analytical theories with the mc simulations .      the weak coupling regime ( @xmath98 )
is again characterized by the mean - field pb equation . in this case however , no closed - form solution is available for the nonlinear pb theory .
we thus take full recourse to the numerical methods appropriate to the problem . in the two - cylinder geometry ,
the axial symmetry is broken but the translational symmetry along the @xmath57 axis remains intact ; hence one has to deal with a two - dimensional pb equation in a finite bounding square @xmath99 , which is then solved numerically @xcite .
we shall focus on the interaction force between the cylinders , which can be evaluated via the electrostatic stress tensor in a standard manner @xcite .
note that the wc scheme used here based on the numerical solution of the pb equation allows for a full analysis of the dielectric discontinuity effects .
( later on in our sc analysis , which will be based on analytical methods , we shall restrict our discussion to only the first - order image approximation , as a full analysis of dielectric discontinuity effects in the sc limit is not yet available . ) here again the counterion size effect will be taken into account approximately by setting the distance of closest approach to each cylinder as @xmath89 .
fluctuation contribution around the mean - field solution in cylindrical geometry is non - trivial to calculate @xcite but is by construction small and we thus skip its detailed analysis
. it does however depend crucially on the dielectric discontinuity at the cylinder surface .    in general , the wc force obtained between two like - charged cylinders
is repulsive and decays monotonically with their interaxial separation @xmath93 .
the mean - field prediction for the rescaled force per unit length defined via @xmath100 is shown in fig .
[ fig : forcedelta_wc ] as a function of the interaxial separation @xmath93 and for different values of the dielectric discontinuity parameter @xmath38 ( see eq .
( [ eq : jump ] ) ) .
note that in contrast to the case of a single cylinder studied in the previous section , the dielectric discontinuity is expected to matter in the two - cylinder geometry even within the mean - field approximation .
this is due to the axial symmetry breaking which leads to electric field penetrating into the cylindrical cores . as seen in the figure
, the wc mean - field results compare very well with the mc data ( to be discussed later ) at sufficiently large interaxial separations .
the dielectric discontinuity has a relatively small effect on the interaction force .
it leads to an increased repulsion ( and also larger deviations from the mean - field prediction ) at small interaxial separations but its effects diminish at larger separations .
as noted before , the sc theory follows systematically from the leading order contribution to the partition function in the limit of infinite coupling parameter @xmath18 @xcite .
the corresponding sc free energy involves only interactions between single counterions and the charged cylinders as well as the direct electrostatic interaction between the cylinders themselves ( see refs .
@xcite ) . in the presence of the dielectric discontinuity effects ,
the sc free energy per counterion can be written as @xmath101\rmd{{\mathbf r}}. \label{f1}\ ] ] here @xmath102 is the total number of counterions and @xmath103 is the electrostatic energy due to the interaction between the two cylinders in the absence of any counterions , which contains both the direct coulomb interaction between their surface charges as well as the contribution due to their images that results from the polarization of their dielectric cores .
furthermore , @xmath104 is the energy due to the interaction between a single counterion and the surface charges on both cylinders and includes contributions from the corresponding image charges as well .
finally , @xmath105 is the image self - energy of counterions , i.e. , the contribution from the interaction of an individual counterion with its own image charges in both cylindrical cores .
note that the volume integration should be again taken over the space available to counterions , i.e. , inside the confining square box of lateral size @xmath95 excluding the two cylindrical cores .
note also that on the sc level , the counterion excluded - volume repulsions are absent ( as implied by the fact that strongly coupled counterions are highly isolated within large correlation holes in the sc limit ) and only the excluded - volume interaction between the counterions and the cylinders will be present .
in fact , the counterion size effects can be accounted for exactly within the sc theory via a hard - core repulsion which simply amounts to setting the effective ( hard core ) cylinder radius as @xmath106 @xcite .    as in the single cylinder case
all charges are coupled with the interaction kernel , composed of direct and image part , but with the difference that now the image part of the kernel @xmath107 corresponds to two dielectric cylinders separated by the distance @xmath93 .
this kernel should satisfy the electrostatic boundary conditions on both cylindrical surfaces ( [ bc ] ) , which leads to very complicated numerical expressions and is unfortunately not available in a closed analytical form .
our aim is nevertheless to give a simple approximate analytical expression for the final results , so we need to make an approximation at this step . in what follows we focus on the first - order - image approximation , where we neglect higher orders of inter - cylindrical image interaction .
this approximation implies that the image green s functions of the two cylinders can be written as the sum of the green s functions of the isolated cylinders , eq .  ( [ vim ] ) , @xmath108 where @xmath109 is the one - cylinder image kernel , eq .
( [ vim ] ) . here ,
@xmath110 , @xmath111 are distances centered at the first and @xmath112 , @xmath113 at the second cylinder , for instance , @xmath114 with @xmath115 a unit vector pointing in the @xmath116 direction , i.e. , horizontally on fig .
[ fig : model_2cyln ] .
recall that the surface charge is uniformly distributed on both cylinders so that the fixed charge distribution is composed of two single cylinder distributions @xmath117 where @xmath118 is the charge distribution of a single cylinder , eq .
( [ n0 ] ) . with these assumptions
one can then evaluate all interaction terms in eq .
( [ f1 ] ) analytically as follows .
the cylinder - cylinder interaction energy per counterion can be obtained as @xmath119 furthermore the image self - energy of a single counterion interacting with the dielectric cores of both cylinders can be derived in the form @xmath120 , \label{wself}\end{aligned}\ ] ] where @xmath80 is defined in eq .
( [ eqw ] ) .
finally , the single particle cylinder - counterion contribution can be written as the sum of the interaction with each cylinder separately and assumes the form @xmath121 the cylinder - counterion energy contribution from the @xmath15-th cylinder ( with @xmath122 ) is itself composed of three parts , @xmath123 which correspond respectively to three parts of the green s function , i.e. , @xmath124 , within the approximation implied by eq .
( [ vim2 ] ) and will be evaluated as follows .
the first contribution in eq .
( [ eq : v_contrib ] ) , @xmath125 , is the direct interaction of the bare surface charge of each cylinder with a single counterion , i.e. , for the the @xmath15-th cylinder , @xmath126 where @xmath127 is the radial distance from the center of the @xmath15-th cylinder to the counterion , and @xmath22 is the manning parameter as defined in eq .
( [ eq : xi_m ] ) . the second contribution , @xmath128 corresponds to the interaction of a counterion image charge with the surface charge on the _ same _ cylinder , which is thus _ zero _ by symmetry reasons as already discussed in the one - cylinder case , eq .
( [ vself ] ) .
this contribution can be also considered as the interaction between the image of the surface charge on the same cylinder with the counterion .
finally , we have the cross contributions due to the interaction of a counterion image in one cylinder with the surface charge on the _ other _ cylinder .
formally , it follows for each cylinder as ( see appendix [ app : cross ] ) @xmath129 for @xmath130 , where @xmath131 is the radial distance from an axis shifted by @xmath132 from the center of the @xmath15-th cylinder toward the axis of the other cylinder , see fig .
[ fig : model_2cyln]b , i.e. , @xmath133 or in cartesian coordinates , @xmath134    the sc free energy thus follows by inserting the above expressions into eq .
( [ f1 ] ) in the form @xmath135- \nonumber \\      & -&\xi(1-\delta)\ln\,\rho_1 ^ 2\rho_2 ^ 2-      \xi      \delta\,\ln\,{\rho_1^{*}}^2{\rho_2^{*}}^2\biggr\}\rmd{{\mathbf r}}. \label{free}\end{aligned}\ ] ] the corresponding sc force between the two cylinders is obtained simply by differentiating the free energy with respect to @xmath93 @xmath136 which may be expressed in dimensionless units according to eq .
( [ forcedef_wc ] ) .
the sc density profile follows from the general expression ( [ eq : density_ci ] ) , which is the same as the integrand in ( [ free ] ) , viz .
@xmath137- \label{dens2}\\      & -&\,\xi(1-\delta)\ln\,\rho_1 ^ 2\rho_2 ^ 2-      \xi      \delta\,\ln\,{\rho_1^{*}}^2{\rho_2^{*}}^2\biggr\}\nonumber.\end{aligned}\ ] ] where the normalization prefactor is determined from the electroneutrality condition .
note that the forms of the sc free energy , eq .
( [ free ] ) , and the sc counterion density profile , eq .
( [ dens2 ] ) , reflect the fact that , within the first - order image approximation ( and using the appropriate volume constraints for counterions ) , the two charged dielectric cylindrical cores can be equivalently replaced by four parallel lines of charges , two of which are located along the central axis for each of the cylinders with ( renormalized ) linear charge density @xmath138 and the two other lines of charges with linear charge density @xmath139 are placed along two axes shifted by @xmath132 from the central axis for each of the cylinders , see fig .
[ fig : model_2cyln ] .
note that the image charge language is applicable for two charged dielectric cylinders in the absence of pointlike charges ( counterions ) @xcite , and can be exactly solved by applying higher orders of images .
but in order to be consistent with the first - order image counterion - cylinder approximation , we have used the same image kernel ( [ vim2 ] ) , which accounts for the first - order image approximation between the two cylinders as well ( [ f1 ] ) .     at interaxial separtaion @xmath140 and in the absence of image charges ( @xmath35 ) .
the density is shown across an arbitrary plane perpendicular to the two cylinders in a color - coded fashion as indicated on the graph .
the small depletion zone around each cylinder is due to the counterion volume exclusion as counterions are assumed to have a finite radius of @xmath141 . ]
0.25 cm     at interaxial separtaion @xmath140 and in the absence of image charges ( @xmath35 ) .
the density is shown across an arbitrary plane perpendicular to the two cylinders in a color - coded fashion as indicated on the graph .
the small depletion zone around each cylinder is due to the counterion volume exclusion as counterions are assumed to have a finite radius of @xmath141 . ]
0.25 cm     ( a )    0.25 cm     ( b )    0.25 cm     ( c )    0.25 cm      in order to proceed , let us first consider the case where there is no dielectric mismatch between the cylinders and the surrounding medium , i.e. , @xmath35 ; the image charges will thus be absent .
the sc interaction in this case has been considered also in previous works @xcite .
we shall reproduce some of these previous results for the sake of completeness , but will also provide several new results , including a direct comparison with mc simulations for the force dependence as a function of separation at different values of the coupling parameter as well as a global interaction phase diagram , which have not been considered previously .    in the absence of image charges , the sc free energy ( [ free ] ) reduces to a simple form as @xmath142 where the first term is again the bare repulsion between the two cylinders and the second term contains energetic and entropic contributions from counterions .
the results for the sc force that follow from eq .
( [ free0 ] ) are shown in fig .
[ fig : force0 ] for a few different manning parameters , where we have accounted for the finite counterion radius by choosing @xmath143 .
the sc force is repulsive at small separations and becomes attractive beyond the equilibrium interaxial distance @xmath144 where the force vanishes .
the sc attraction is mediated by counterions that are strongly coupled to both cylinders in this limit and are thus accumulated mainly in the region between the two cylinders as can be seen from the counterion distribution ( [ dens2 ] ) in fig .
[ fig : densdelta0 ] .
this is a direct consequence of the energy contributions included in the second term in eq .
( [ free0 ] ) @xcite .    by decreasing the separation between the cylinders ,
the entropic osmotic pressure from counterions sandwiched between the two cylinders becomes increasingly important and the effective interaction becomes repulsive .
this repulsion is reduced as the surface - surface separation , @xmath145 , becomes smaller than the counterion diameter , @xmath146 as the counterions are depleted from the intervening region due to excluded - volume effects .
this behavior due to the presence of the counterion depletion interaction has been investigated thoroughly in the case of cylinders with helical charge pattern and we shall not delve on it any further here @xcite .    as seen in fig .  [
fig : force0 ] the entropic repulsion effects at small separations decrease when the manning parameter @xmath22 is increased .
for very large manning parameters , the free energy is dominated by purely energetic contributions and does not contain any temperature effects .
in fact the limit @xmath147 formally corresponds to the zero temperature limit as we have already taken the limit @xmath18 appropriate within the sc approximation .
the connection between the sc limit and the zero temperature limit has been analyzed in detail in refs . @xcite .
the asymptotic form of the sc force in the limit of large @xmath22 can be obtained as @xmath148 which is shown as a black solid line in fig .
[ fig : force0 ] .
the equilibrium hard - core surface - surface separation @xmath149 tends to zero as @xmath150 as @xmath147 @xcite .
note that even for moderate values of the manning parameters , @xmath151 , the sc theory predicts a closely packed _ bound state _ with a relatively small surface - surface separation @xmath152 .
let us now consider the situation where the coupling parameter has a finite value . in this case we shall study the system using mc simulations as detailed in appendices [ app : mc ] and [ app : mc_force ] .
the mc results for the interaction force between like - charged cylinders ( with no dielectric mismatch ) are shown in fig .
[ fig : force0mc ] along with the wc and sc predictions . as
seen in fig .
[ fig : force0mc]a , the wc results agree with the simulation data for a sufficiently small coupling parameters @xmath153 as expected . furthermore , by increasing the coupling parameter , the mc results start deviating strongly from the wc prediction and tend to the sc prediction , where a reasonable agreement is obtained already at coupling parameters of the order @xmath154 .
note that the wc and sc predictions bracket the simulation data and thus establish upper and lower bounds for the interaction force at any realistic value of the coupling parameter .
( a )    0.25 cm     ( b )     ( c )    in figs .
[ fig : force0mc]b and c we show the results for larger manning parameters of @xmath155 and 100 , respectively . as manning parameter becomes larger , a larger @xmath10 is required in order to achieve the same level of agreement with the sc prediction at a given separation and thus the convergence to the sc limit becomes weaker .
this effect is more pronounced at larger interaxial separations .
in fact , as is generally known @xcite , the sc effects become more dominant and the sc theory becomes increasingly more accurate at smaller surface - surface separations . in general , the applicability regime of the sc theory ( obtained in the limit of @xmath18 ) to the situations where the coupling parameter @xmath10 is finite can be specified via simple validity criteria as applied successfully in several previous studies @xcite . in what follows
we shall derive the particular form of the sc validity criterion for the case of two like - charged cylinders .    it should be noted first that the sc theory ( obtained as an exact asymptotic theory in the limit of @xmath18 ) contains only contributions stemming from the interaction of individual counterions with charged objects and thus the counterion - counterion repulsions and other higher order many - body effects are absent .
these effects however become increasingly more important as the coupling parameter is decreased , depending crucially on the precise value of all the system parameters . in the case of two cylinders ,
as noted before , counterions are mostly accumulated in the intervening region between the two cylinders ( fig .
[ fig : densdelta0 ] ) where the typical counterion spacing can be estimated by stipulating the local electroneutrality condition as @xmath156 or in rescaled units as @xmath157 .
it is thus evident that as the manning parameter becomes larger at a fixed coupling parameter the spacing between counterions tends to become small .
hence , a larger effect due to counterion - counterion repulsions and thus larger deviations from the sc theory would be expected to be observed at larger @xmath22 . in other words ,
the sc theory is expected to overestimate the counterion density in the intervening region at any _ finite _ value of the coupling parameter , which is consistent with the general result that the sc theory gives the upper bound for the density profile of the counterions @xcite .
one can thus identify the validity regime of the sc theory by estimating the counterion - counterion contributions and comparing them with the sinlge - particle counterion - cylinder contribution at a given interaxial separation .
this can be done by noting that the ground state configuration ( obtained for @xmath158 and @xmath18 ) is predicted within the sc theory to be the configuration where all counterions are localized between the two cylinders and are lined - up along the @xmath57 axis as obviously favored by the electrostatic interaction energy .    when the coupling parameter is reduced to a finite value ,
repulsions between counterions become more important and some counterions tend to escape from the intervening region to the external side of the cylinders ( shown by an arrow in fig .
[ fig : diagrams0]a ) . in order to prevent this from happening
, the coupling parameter @xmath10 has to be large enough . by comparing the energy of the ground state with the energy of a single defect obtained by allowing one counterion to move to the external side of one of the cylinders
, it follows analytically that for a large @xmath22 and in the absence of image charges , one should have @xmath159 these criteria indeed capture the trends observed in fig .
[ fig : force0mc ] when viewed in terms of the three main parameters describing the system , i.e. , @xmath10 , @xmath22 and @xmath93 .
they automatically cover the one proposed before based on the surface - surface separation between the cylinders , where one requires that the surface - surface separation between the cylinders must be smaller than the spacing between counterions , @xmath160 , or in rescaled units @xcite @xmath161 this latter criterion is inspired by the observations in the case of charged planar walls where it can indeed be derived from systematic analysis of higher - order corrections to the sc theory @xcite .
it can describe the validity regime of the sc predictions for the equilibrium separation between strongly coupled cylinders that form a closely packed bound state @xcite , but when the applicability of the sc theory for the interaction force is considered , the criteria ( [ sccriterion ] ) should be used instead . the difference between the new criteria ( [ sccriterion ] ) and the one in eq .
( [ eq : crit_old ] ) is in fact related to the fundamental difference between the topology of the space available to counterions in the case of two cylinders as compared to two planar walls : in the latter case counterions can not  escape " from the intersurface gap and the ground state of the system has a universal two - dimensional configuration @xcite , while in the former the ground state corresponds to a one dimensional arrangement of counterions and the whole space is available for thermal excitations from the ground state .
the criteria ( [ sccriterion ] ) are more stringent and cover also the situation where the cylinders are placed at large interaxial separations .
the ground state of neutralizing counterions between two charged cylinders ( without the effects of the dielectric mismatch ) has been investigated extensively by arnold and holm in ref .
@xcite , where a validity criterion for the sc theory has been proposed as @xmath162 ( when expressed in rescaled units ) based on computer simulations .
this agrees with our estimate in eq .
( [ sccriterion ] ) within 16% of the numerical prefactor .
let us now turn our attention to the behavior of the equilibrium interaxial separation or the so - called bound - state separation , @xmath144 , where the interaction force between the two cylinders vanishes . as noted
before , the sc theory is generally more accurate at smaller separations . a closer inspection of fig .
[ fig : force0mc ] shows that the sc prediction for the bound - state separation agrees with the mc data even outside the regime set by the sc criteria ( [ sccriterion ] ) , which reiterates the point discussed above that a less stringent criterion such as eq .  ( [ eq : crit_old ] ) would be sufficient to describe the validity of the sc prediction for the bound - state separation . in fig .
[ fig : diagrams0]b the bound - state interaxial separation is plotted as a function of the manning parameter , which shows that a reasonable agreement between mc data and the sc theory for this quantity can be achieved for a coupling parameter as small as @xmath163 , in agreement with previous results in ref .
as seen in the figure , the bound - state separation diverges ( i.e. , the cylinders unbind ) for manning parameter approaching a minimum value of @xmath164 @xcite below which the two cylinders merely repel each other .
figure [ fig : diagrams0]b can be also thought of as a phase diagram identifying the attraction ( the region above the lines ) and repulsion ( the region below the lines ) regimes between two like - charged cylinders in terms of the parameters @xmath144 and @xmath22 for a given coupling parameter @xmath10 .    in fig .
[ fig : diagrams0]c we show a global phase diagram in terms of the parameters @xmath10 and @xmath22 .
the boundary lines are determined from mc simulations ( open symbols ) and separate the region ( above the lines ) where the interaction force becomes attractive at some finite interaxial separation between the cylinders and the region ( under the lines ) where the force never becomes attractive at _ any _ interaxial separation . in other words ,
the boundary lines themselves correspond to the parameter values where the force - distance diagrams exhibit a global minimum at @xmath165 , i.e. the force curve only touches the abscissa .
they were determined by bisection procedure for a series of mc runs .
if the cylinder parameter values are chosen to describe dna , we have @xmath166 and @xmath167 for @xmath26 valency counterions , which are shown as solid circles in fig .
[ fig : diagrams0]c .
this suggests that dna - dna interaction in the case of monovalent counterions falls well within the repulsive region , while with divalent counterions and beyond one can expect to observe an attractive force , in qualitative agreement with recent experiments @xcite .     and the dielectric discontinuity parameter @xmath36 at interaxial separtaion @xmath140 .
the density is shown across an arbitrary plane perpendicular to the two cylinders in a color - coded fashion as indicated on the graph .
the small depletion zone around each cylinder is due to the counterion volume exclusion as counterions are assumed to have a finite radius of @xmath141 .
a large depletion zone develops due to the influence of image charges as the coupling parameter increases from @xmath168 to 5000 ( a - d ) .
, title="fig : " ] ( a ) @xmath168   and the dielectric discontinuity parameter @xmath36 at interaxial separtaion @xmath140 .
the density is shown across an arbitrary plane perpendicular to the two cylinders in a color - coded fashion as indicated on the graph .
the small depletion zone around each cylinder is due to the counterion volume exclusion as counterions are assumed to have a finite radius of @xmath141 .
a large depletion zone develops due to the influence of image charges as the coupling parameter increases from @xmath168 to 5000 ( a - d ) .
, title="fig : " ] ( c ) @xmath169    0.25 cm     and the dielectric discontinuity parameter @xmath36 at interaxial separtaion @xmath140 .
the density is shown across an arbitrary plane perpendicular to the two cylinders in a color - coded fashion as indicated on the graph .
the small depletion zone around each cylinder is due to the counterion volume exclusion as counterions are assumed to have a finite radius of @xmath141 .
a large depletion zone develops due to the influence of image charges as the coupling parameter increases from @xmath168 to 5000 ( a - d ) .
, title="fig : " ] ( b ) @xmath170   and the dielectric discontinuity parameter @xmath36 at interaxial separtaion @xmath140 .
the density is shown across an arbitrary plane perpendicular to the two cylinders in a color - coded fashion as indicated on the graph .
the small depletion zone around each cylinder is due to the counterion volume exclusion as counterions are assumed to have a finite radius of @xmath141 .
a large depletion zone develops due to the influence of image charges as the coupling parameter increases from @xmath168 to 5000 ( a - d ) .
, title="fig : " ] ( d ) @xmath171     ( a )    0.25 cm     ( b )     ( c )      we now consider the case of two charged cylinders with a dielectric core whose dielectric constant may be in general different from that of the surrounding medium . in this case , the induced image charges ( due to the polarization of the dielectric cores ) can play an important role as well .
however , as we discussed before ( fig .  [
fig : forcedelta_wc ] ) , the image - charge effects turn out to be small in the wc regime where the coupling parameter is small ( e.g. , when counterions are monovalent ) .
this may be the reason why such dielectric effects have not yet been throughly investigated in the particular case of two charged cylinders .
the situation turns out to be quite different in the opposite limit of strong coupling as we will see in this section .
this is tied in with a fundamental difference between counterion distribution in the sc limit as opposed to the wc limit : while in the former limit one deals with highly isolated single counterion close to a charged surface , in the latter limit individual nature of counterions fades in the wake of dominant collective many - body effects ( as counterions form a diffuse and uncorrelated ionic cloud around the cylinders ) and thus the polarization effects are highly reduced .
if the dielectric cores have a smaller dielectric constant than the medium ( @xmath84 ) , as relevant to most macromolecules in water , the induced image charges will have the same sign as the counterions and thus tend to cause depletion of counterions from the vicinity of the dielectric cores .
this behavior is shown in fig .
[ fig : densdelta ] , where the sc density is plotted across an arbitrary plane perpendicular to the two cylinders in a color - coded fashion . as can be seen , the depletion effect is enhanced as the coupling parameter , @xmath10 , is increased such that a large  depletion zone " develops that eventually encloses both cylinders .
( note that the sc free energy in the presence of a dielectric discontinuity , eq .
( [ free ] ) , has an explicit dependence on the coupling parameter @xmath10 as already known from the case of planar charged dielectric slabs @xcite ) .
in fact , the image depletion at larger couplings turns out to have a similar effect qualitatively as the counterion volume exclusion , i.e. , as if the counterion radius @xmath67 is renormalized to a large effective value .
note that the region with highest concentration ( shown in red color ) is  squeezed " and eventually splits in two regions shifted away from the common @xmath172 plane that passes through the axes of the two cylinders .
the splitting appears when the second derivative of the density ( [ dens2 ] ) with respect to the @xmath173 coordinate becomes positive at the midpoint , @xmath174 . using the small - distance approximation @xmath175 from eq .
( [ eqi ] ) , we can estimate the distance @xmath93 at which the density splitting happens as @xmath176 thus for the manning parameter @xmath155 and at separation @xmath140 , the splitting is estimated to occur at @xmath177 as is the case in fig .
[ fig : densdelta ] .    in figs .
[ fig : forceqdelta]a and b , we show the sc force per unit length ( solid lines ) as obtained from eqs .
( [ free ] ) and ( [ free_to_force ] ) along with the mc simulation results ( symbols ) for two different manning parameters and several different values of the coupling parameter . as seen upon increasing the coupling parameter , the repulsive peak at small separations shifts towards larger values of the interaxial separation @xmath93 , which reflects the emergence of a depletion zone due to repulsions between counterions and their image charges as discussed before . as a result ,
the main contribution to the inter - cylinder interaction at small separations comes from the bare interaction between the two cylinders and the force is approximately given by the bare repulsion force @xmath178 as given by eq .
( [ f00 ] ) in the appendix ( dashed line , fig .
[ fig : forceqdelta]b ) . at larger separations ,
we still find strong sc - type attraction even in the presence of image charges .
one should thus note that in general the repulsive forces are enhanced in the effective interaction in the sc limit especially at small separations as compared to a homogeneous system with no image charges .
therefore , unlike in the wc limit ( see fig .  [
fig : forcedelta_wc ] ) , the introduction of image charges may lead to qualitative changes in the behavior for a strongly coupled system .
in general , the sc results for the interaction force show better agreement with the simulations for larger @xmath10 and smaller separation distances @xmath93 . especially the results for the bound - state
interaxial separation show very good agreement between the sc results ( solid lines ) and the mc data ( symbols ) as shown in fig .
[ fig : forceqdelta]c ( compare with fig .
[ fig : diagrams0]b ) .
in plan - parallel geometry the repulsive self - image interaction compresses counterions toward the mid - plane and hence expands the range of the sc regime @xcite . here in the cylindrical geometry such mechanism can not work , since counterions can escape from inter - cylindrical region .
therefore , it seems that dielectric image effects do not have a noticable influence on the range of validity of the sc approximation .
note that just as in the case with @xmath35 ( blue line in fig .
[ fig : diagrams0]c ) , the simulation results show no attraction between the cylinders at too low electrostatic coupling parameters .
the attraction can only appear once the coupling parameter exceeds a threshold value .
this threshold value is shown as a red line in the global phase diagram in fig .
[ fig : diagrams0]c for the dielectric discontinuity parameter @xmath179 .
note that bare cylinder - cylinder repulsion is enhanced in the case of images ( @xmath92 ) , eq.([f1 ] ) .
therefore , larger coupling parameter @xmath10 is needed to get attraction at smaller manning parameters compared to no - image case ( red line stays above blue for @xmath180 and @xmath36 in fig .
[ fig : diagrams0]c ) . for larger
manning parameters , the attractive interaction between the counterion image and the cylindrical charges ( represented by the term @xmath181 , eq .
( [ vcross ] ) ) becomes more important and , hence , attraction can appear even at lower coupling parameters ( red line is below the blue line for @xmath182 on fig .  [
fig : diagrams0]c ) .
we have analyzed the electrostatic interaction between two like - charged cylindrical macromolecules surrounded by counterions .
we have in particular examined the role of image charges when the cylinders have a dielectric constant which is different from that of the surrounding medium within the the weak- and strong - coupling frameworks elaborated first by netz and coworkers @xcite .
the two limits are defined in terms of a single coupling parameter @xmath10 .
the weak coupling limit , or equivalently the mean - field theory , is relevant when the valency of counterions and/or the macromolecule surface charge density is small and/or when the dielectric constant of the medium is large enough .
it is based on the poisson - boltzmann equation which turns out to give purely repulsive interactions between two - like charged cylinders .
the image - charge effects in this case turn out to be small even at small separations and quickly diminish as the inter - cylinder separation is increased .    on the contrary ,
the strong coupling limit is relevant when the valency of counterions is large and/or the macromolecule surface charge density is large and/or when the dielectric constant of the medium is small enough .
the sc approach is based on a single - particle description which is obtained as an exact limiting result for large coupling parameters @xmath18 @xcite . while some aspects of the sc theory for cylindrical macromolecules were studied in previous works @xcite , other aspects have remained unadressed .
our work is aimed specifically at addressing the effects due to image charges , which indeed turn out to be quite significant in the sc limit . using a generalized sc theory and extensive mc simulations
we have shown that the counterions are strongly depleted away from the cylindrical cores due to repulsion from image charges leading to stronger effective repulsions between the two cylinders at smaller separations .
the counterion - mediated attraction will be present at intermediate to large separation when the coupling parameter is sufficiently large .
these features are in marked contrast with those obtained with the wc theory .
we have mapped a global phase diagram where attractive and repulsive interactions emerge between two like - charged cylinders with or without the image - charge effects .
here we have employed a first - order image approximation to deal with the image - charge effects within the sc theory ( note , however , that the numerical scheme used within the wc limit allows to account for the full effect of the dielectric discontinuity ) . generally , the treatment of the image charges in nontrivial geometries becomes very complicated and the two - cylinder model is no exception .
the first - order image approximation ( [ vim2 ] ) has the obvious advantage that it allows for an analytical treatment of the system .
the study of the full effect due to the dielectric discontinuity in the two - cylinder model requires more advanced numerical and analytical developments that might become available in the future .
another interesting effect which can be investigated in the present context is that of the additional salt that may be present besides the neutralizing counterions .
recent generalizations @xcite should allow for a systematic study of salt screening effects within the wc - sc framework .
r.p . would like to acknowledge the financial support by the slovenian research agency under contract nr .
p1 - 0055 ( biophysics of polymers , membranes , gels , colloids and cells ) .
m.k . would like to acknowledge the financial support by the slovenian research agency under the young researcher grant .
is a newton international fellow .
in this section , we explain how the cross energy term ( [ vcross ] ) may be evaluated in the case of two charged cylinders ( section [ subsec : sc_2cyln ] ) .
this requires employing some mathematical identities as we shall explain below .
note that here the position vectors centered at two different cylinders appear simultaneously , @xmath130 , @xmath184 inserting the image kernel @xmath109 , eq .
( [ vim ] ) , and the cylinder charge density @xmath118 , eq .
( [ n0 ] ) , we first transform the coordinates @xmath127 and @xmath185 centered at the @xmath15-th cylinder to the those centered at the @xmath186-th cylinder using the graf s addition theorems @xcite ,
@xmath187 with @xmath188 and @xmath189 being integers and @xmath190i_n(k\rho_\nu)\hspace{3ex}\\ c_{mn}(k\rho_\nu)&=&[k_{m+n}(kr)+k_{m - n}(kr)]i_n(k\rho_\nu)\nonumber\\      & & \times\ > ( 1-\frac{1}{2}\delta_{n0}).\end{aligned}\ ] ] we can then carry out the integration in ( [ app : vcross_def ] ) over @xmath191 , @xmath192 and @xmath193 .
integration over @xmath193 produces a @xmath194 term , where @xmath195 is the upper limit and thus we obtain @xmath196 in the next step we integrate over the wave - vector @xmath197 , where the integrand contains the rapidly oscillating factor @xmath198 , which supresses contributions of the integrand in whole the range , except at @xmath199 .
formally , we use the following mathematical identity @xmath200 where @xmath201 is the residue of function @xmath202 at @xmath203 .
this integral is finite only if @xmath202 behaves as @xmath204 for @xmath205 and goes to 0 for @xmath206 .
this statement can be proven by considering the split of function @xmath202 into a @xmath204 term and a remainder , i.e. , @xmath207 .
the first part can be integrated straightforwardly , whearas the remainder yields 0 according to the riemann - lebesgue lemma @xcite .
this finally leads us to @xmath208 which may be simplified further using the identity @xmath209\nonumber\end{aligned}\ ] ] and yields the result given in eq .
( [ vcross ] ) .
we performed monte - carlo simulations in order to study the system of one and two charged dielectric cylinders at arbitrary coupling parameters , which thus go beyond the weak and strong coupling theories .
all simultations were performed in the canonical ensemble ( @xmath210 ) using standard metropolis algorithm @xcite . in the case of one charged cylinder
the system ( including the neutralizing counterions ) is enclosed in a cylindrical simulation box of outer radius @xmath211 whereas in the case of two cylinders the system is enclosed in a square box of lateral size @xmath96 ( wc ) and 100 ( sc ) ( note that in this work we focus on the regime of large manning parameters where the lateral size and shape of the confining box becomes unimportant @xcite ) .
the simulation box height is assumed to have a finite value @xmath212 , which in terms of other physical parameters by invoking the global electroneutrality condition : in the case of a single cylinder it is given by @xmath213 , and in the case of two cylinders @xmath214 .
we use periodic boundary conditions in the @xmath57 direction by replicating the main simulation box infinitely many times in that direction .
the energy of the system for a given configuration of the cylinder(s ) is composed of two main parts @xmath215 the first term is the counterion - cylinder interaction energy as given by eqs .
( [ eq_w0c ] ) and ( [ eq : w0c_two ] ) in the text .
the second term in ( [ total_en ] ) is counterion - counterion interaction energy which includes also the contribution from interactions with counterion image charges .
due to the periodic boundary conditions used in the simulations , this latter term involves infinite summation series @xcite .
this is because the counterions and image charges in the main simulation box interact also with their periodic  copies " as obtained by the replication of the main simulation box to an infinite number of simulation box copies .
these summations can be evaluated as explained in the forthcoming sections .
the interaction energy between two given counterions @xmath216 and @xmath217 can be written as @xmath218 the first equation corresponds to the hard - core repulsion between two counterions overlapping counterions .
the second equation gives the image self - energy of the @xmath216-th counterion , i.e. , it takes into account the interaction of a given counterion with its dielectrically induced images in the polarizable dielectric core . the third equation in ( [ ionion ] ) gives the interaction between two different counterions @xmath216 and @xmath217 .
it is written as the sum of two different contributions which will be calculated below : i ) the coulomb interaction , @xmath219 , between the @xmath216-th ( located at @xmath220 ) and the @xmath217-th counterions and all its periodic copies ( located at @xmath221 for integer @xmath222 ) , @xmath223 and ii ) the contribution involving interactions with the dielectrically induced image charges @xmath224 that will be defined later in this appendix .      the first contribution @xmath219 can be written as @xmath225 where we have introduced two different summation schemes depending on the radial separation @xmath226 between the two position vectors @xmath220 and @xmath227 in order to increase the convergence of the sum over the long - range electrostatic interactions . if radial distance @xmath228 is smaller than a threshold @xmath229 , we first sum up @xmath230 terms directly and then use an approximation to estimate the remaining part of the series ( from @xmath231 to infinity ) .
we refer to this as the  direct " summation scheme .
this method is different from the so - called sperb @xcite summation scheme which involves calculating various special functions . for our purposes of relative precision up to @xmath232 we found it more efficient and
can be made in principle as accurate as necessary .
it follows as @xmath233 where the last sum over @xmath234 is equal to @xmath235 , where @xmath236 is digamma function , and needs to be calculated only once in the course of the simulations . the term @xmath237 is defined as @xmath238 this scheme gives an error of the order @xmath239 . for larger radial separations ( @xmath240 ) ,
we use the so - called lekner summation scheme @xcite as it converges more rapidly .
it is given by @xmath241 where @xmath242 accounts for the offset generated at the threshold , i.e. , @xmath243    we choose the threshold value as @xmath244 and use a cut off of @xmath245 in the direct scheme and sum up to @xmath246 $ ] terms in the lekner scheme , which give a relative error smaller than @xmath232 .
let us now consider the second contribution that enters eq .
( [ ionion ] ) , i.e. , @xmath247 , which is the interaction energy obtained by summing up all contributions that involve dielectrically induced image charges throughout different periodic copies of the main simulation box along the @xmath57 axis .
it can be written in the case of one cylinder as @xmath248 where @xmath249 is a unit vector pointing in @xmath250 direction . inserting the image kernel ( [ vim ] ) into the above equation , one ends up with an infinite summation of terms including @xmath251 , which can be expressed in terms of @xmath252 functions , i.e. , @xmath253 this enables us to simply integrate over the wave - vector variable @xmath197 and thus obtain @xmath254 where @xmath255 .
the @xmath256 term needs to be handled carefully by taking the limit @xmath257 as @xmath258 using the summation identity as used in eq .
( [ app : vcross_2 ] ) , we find the final expression @xmath259+\nonumber\\ & & + \frac{4}{l_z}\sum_{m=0}^\infty\sum_{n=1}^\infty\xi_{m}(k_n a)k_m(k_n\rho_i)k_m(k_n\rho_j)\nonumber\\ & & \times\>\cos\ , k_n\delta z\ , \cos\ , m\delta \varphi .
\label{app : imagesum}\end{aligned}\ ] ] this expression is the most computationally expensive part of the simulations , and its efficiency decreases for @xmath260 , i.e. , for counterions near the cylindrical core . in the simulations we truncate these summations in such a way that the relative cut - off error is less than @xmath261 .    in the case of two cylinders ,
we use the first - order image approximation in the simulations as well , and thus the image - charge contributions are simply the sum of the contributions for each cylinder separately , i.e. , @xmath262 where indices 1 and 2 represent coordinates centered at the first and the second cylinder , respectively .
in this section , we explain the method used to directly evaluate the force between two charged cylinders within mc simulations . the force acting on a given cylinder ( say cylinder 1 )
is composed of three contributions , @xmath263 where @xmath178 is the bare electrostatic force due to the interaction with the other cylinder , @xmath264 , is the electrostatic force due to the interaction of the cylinder with the @xmath216-th counterion ( summed over all counterions @xmath265 ) , and @xmath266 is the force due to the osmotic pressure from the @xmath216-th counterion .
the first contribution follows from the interaction energy @xmath103 , eq .
( [ w00_delta ] ) , as @xmath267 which gives @xmath268 the second contribution follows by differentiating @xmath104 , eq .
( [ eq : w0c_two ] ) , with respect to @xmath93 by assuming that the cylinder 1 and the @xmath269-th counterion are fixed , i.e. , @xmath270 which after some algebra gives @xmath271.\qquad\ ] ] here @xmath272 and @xmath273 are the radial distance and the azimuthal angle of the @xmath216-th counterion position with respect to cylinder 1 , and @xmath274 and @xmath275 are the corresponding coordinates with respect to the shifted axis by @xmath132 from the axis of cylinder 1 toward the axis of the cylinder 2 , fig .
[ fig : model_2cyln ] .
similar definitions apply to @xmath276 , @xmath277 , @xmath278 and @xmath279 .
the osmotic ( third ) term in eq .
( [ force ] ) results from the collisions of counterions with the cylinder .
this contribution amounts to a pressure of counterions exerted on the cylinder surface , which is proportional to the contact density of counterions at the cylinder surface @xmath280 .
thus the infinitesimal force acting on the cylinder is given by @xmath281 .
note that the contact counterion density @xmath282 is given by the number of counterions inside a small box of dimensions @xmath283 next to the cylindrical surface .
the total osmotic force of all ions in simulation cell is then @xmath284 where @xmath285 is the heaviside step function which is @xmath286 only if the @xmath216-th counterion falls inside a shell thickness @xmath287 from the cylinder surface .
the above procedure in principle gives the exact value of the osmotic force when @xmath288 .
it also gives a convenient method to calculate the osmotic contribution within mc simulations by evaluating the force for different values of the shell thickness and estimating the limiting value by extrapolation .
r. r. netz , eur .
j. e * 5 * , 557 ( 2001 ) ; a. g. moreira and r. r. netz , _ ibid _ * 8 * , 33 ( 2002 ) ; phys .
lett * 87 * , 078301 ( 2001 ) ; europhys . lett . *
52 * , 705 ( 2000 ) ; _ ibid _ * 57 * , 911 ( 2002 ) .
r. h. french , v. a. parsegian , r. podgornik _ et al .
( 2010)in press .
l. guldbrand , b. jnsson , h. wennerstrm , and p. linse , j. chem . phys . * 80 * , 2221 ( 1984 ) ; d. bratko , b. jnsson , and h. wennerstrm , chem .
lett . * 128 * , 449 ( 1986 ) ; j. p. valleau , r. ivkov , and g. m. torrie , j. chem . phys . * 95 * , 520 ( 1991 ) ; r. kjellander , t. kesson , b. jnsson , and s. marelja , j. chem . phys . * 97 * , 1424 ( 1992 ) .                    a. naji and r. podgornik , phys .
e , 2005 , * 72 * 041402 ; r. podgornik and a. naji , europhys .
letts . , 2006 , * 74 * , 712 ; y. s. mamasakhlisov , a. naji , and r. podgornik , j. stat . phys .
* 133 * , 659 ( 2008 ) ; a. naji _ et al .
_ , phys . rev .
* 104 * , 060601 ( 2010 ) .
stevens , m.o .
robbins , europhys . lett .
* 12 * , 81 ( 1990 ) .
a. diehl , m.n .
tamashiro , m.c .
barbosa , y. levin , physica a * 274 * , 433 ( 1999 ) ; m.c .
barbosa , m. deserno , c. holm , europhys .
* 52 * , 80 ( 2000 ) .
kornyshev , s. leikin , j. chem . phys . * 107 * , 3656 ( 1997 ) ; phys .
lett . * 82 * , 4138 ( 1999 ) .
arenzon , j.f .
stilck , y. levin , eur .
j. b * 12 * , 79 ( 1999 ) ; y. levin , j.j .
arenzon , j.f .
stilck , phys .
. lett . * 83 * , 2680 ( 1999 ) ; a. diehl , h.a .
carmona , y. levin , phys .
e * 64 * , 011804 ( 2001 ) .
henle , c.d .
santangelo , d.m .
patel , p.a .
pincus , europhys . lett . * 66 * , 284 ( 2004 ) . m. manghi and r. r. netz , eur .
j. e * 14 * , 67 ( 2004 ) s. liu , m. muthukumar , j. chem .
phys . * 116 * , 9975 ( 2002 ) ; s. liu , k. ghosh , m. muthukumar , j. chem . phys . * 119 * , 1813 ( 2003 ) .
m. muthukumar , j. chem . phys . *
120 * , 9343 ( 2004 ) .
d. chapot , l. bocquet , e.trizac , j. chem . phys . * 120 * , 3969 ( 2004 ) a. g. cherstvy , j. phys . chem .
b * 111 * , 12933 ( 2007 ) ; j. phys . : condens .
matter * 17 * , 1363 ( 2005 ) .
we make use of the commercial software comsol multiphysics which employes a finite element mesh approach in order to solve the nonlinear pb equation in the two - cylinder geometry ( http://www.comsol.com ) .