Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: stability and interactions in biological and soft - matter systems often depends on the underlying properties of electrostatic interactions @xcite . charges on macromolecular surfaces in aqueous environments , as in the case of membranes , self - assembled micelles , globular proteins and fibrous polysaccharides , affect a wealth of functional , structural and dynamical properties @xcite . the traditional approach to charged ( bio)colloidal systems has been the mean - field poisson - boltzmann ( pb ) formalism applicable at weak surface charges , low counter - ion valency and high temperature @xcite . the limitations of this approach become practically important in highly - charged systems where counterion - mediated interactions between charged bodies start to deviate substantially from the mean - field accepted wisdom @xcite . one of the fundamental recent advances in this field has been the systematization of these non - pb effects based on the notions of _ weak _ and _ strong _ coupling approximations . the latter approach has been pioneered by rouzina and bloomfield @xcite , elaborated later by shklovskii _ et al . _ @xcite , levin _ et al . _ @xcite , and brought into final form by netz _ these two approximations allow for an explicit and exact treatment of charged systems at two disjoint _ limiting conditions _ whereas the parameter space in between can be analyzed only approximately @xcite and is mostly accessible solely _ via _ computer simulations @xcite .    in the absence of a general approach that would cover thoroughly all the regions of the parameter space one has to take recourse to various partial formulations that take into account only this or that facet of the problem . in this respect the counterion - only or the one - component coulomb fluid model system has proved to be of substantial value @xcite . heuristically as well as numerically . a proper understanding of the behavior of charged systems would thus start with the analysis of counter - ion distribution around charged macromolecular surfaces , neglecting completely the effects of salt . both the weak and the strong coupling approximations are based on a functional integral or field - theoretic representation @xcite of the grand canonical partition function of a system composed of fixed surface charges with intervening mobile counterions , and depend on the value of a single dimensionless coupling parameter @xmath0 @xcite . the distance at which two unit charges interact with thermal energy @xmath1 is known as the bjerrum length @xmath2 ( in water at room temperature , one has @xmath3 nm ) . if the charge valency of the counterions is @xmath4 then the aforementioned distance scales as @xmath5 . similarly , the distance at which a counterion interacts with a macromolecular surface ( of surface charge density @xmath6 ) with an energy equal to @xmath1 is called the gouy - chapman length , defined as @xmath7 . a competition between ion - ion and ion - surface interactions can be quantitatively measured with a ratio of these characteristic lengths , that is @xmath8 , which is known as the ( netz - moreira ) electrostatic coupling parameter @xcite . the weak coupling ( wc ) regime @xmath9 ( appropriate for low valency counterions and/or weakly charged surfaces ) , is characterized by the fact that the width of the counterion layer @xmath10 is much larger than the separation between two neighboring counterions in solution and thus the counterion layer 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 ( pb ) type is completely justified . on the other hand in the strong coupling ( sc ) regime @xmath11 ( appropriate for high valency counterions and/or highly charged surfaces ) , the mean distance between counterions , @xmath12 , is much larger than the layer width ( _ i.e. _ , @xmath13 ) , indicating that the counterions are highly localized laterally and form a strongly correlated quasi - two - dimensional layer next to a charged surface . in this case , the weak - coupling approach breaks down due to strong counterion - surface and counterion - counterion correlations . since counterions can move almost independently from the others along the direction perpendicular to the surface , the collective many - body effects that enable a mean - field description are absent , necessitating a complementary sc description @xcite . the range of validity of both limiting theories at intermediate values of the coupling parameter has been explored thoroughly in the literature @xcite . formally the weak coupling limit can be straightforwardly identified with the saddle - point approximation of the field theoretic representation of the grand canonical partition function , and is reduced to the mean - field pb theory in the lowest order for @xmath14 . the quadratic fluctuations around the mean field provide a second - order correction to the mean - field solution for small finite @xmath15 @xcite . the strong coupling approximation has no pb - like correlates @xcite since it is formally equivalent to a single particle description obtained from a systematic @xmath16 expansion in the limit @xmath17 , and corresponds to two lowest order terms in the virial expansion of the grand canonical partition function . the consequences and the formalism of these two limits of the coulomb fluid description have been explored widely and in detail ( for reviews , see refs . @xcite ) .    considering the inhomogeneity of charged surfaces in various biological contexts it has always been of interest to investigate not just electrostatic interactions between symmetrical charged surfaces , _ i.e. _ those bearing equal charges of the same sign , but also interactions between surfaces bearing unequal charges or even charges of opposite sign @xcite . this problem has a venerable history starting from the seminal work of parsegian and gingell @xcite who formulated a linearized pb theory of the interactions in the presence of salt . the linearization _ ansatz _ was later generalized in the work of lau and pincus @xcite and ben - yaakov _ et al . _ @xcite who formulated the appropriate non - linear mean - field theory of non - symmetric electrostatic interactions . it is thus our goal in this contribution to show how and to what extent the asymmetry in the distribution of charges on two apposed planar surfaces affects the interactions between macromolecular surfaces carrying them . below we shall present a complete analysis of the asymmetric case in the weak coupling limit , _ i.e. _ the mean - field poisson - boltzmann theory supplemented with a quadratic - fluctuations analysis , as well as in the strong coupling limit _ via _ the asymptotic strong - coupling theory and evaluate how these analytical results compare with extensive numerical simulations . we will show that in their respective regimes of validity ( _ i.e. _ small / large couplings ) both approximations present a very accurate quantitative statistical description of the system . in our model system we consider uniform surface charge distributions on two plane - parallel surfaces ( located at @xmath18 ) given by the surface charge density of the form @xmath19 we may interchangeably use the half - separation @xmath20 , or @xmath21 to identify the surface - surface distance . we assume furthermore that the charge of both bounding surfaces is compensated by mobile counterions of charge valency @xmath4 immersed in an aqueous medium of dielectric constant @xmath22 and distributed in between the two surfaces ( see fig . [ fig : geometry ] ) . we thus neglect all coions . this approximation is relevant for low salt concentrations where the debye screening length is much larger than the scales of interest @xcite . we consider the surfaces as impenetrable to counterions and neglect the dielectric discontinuity across the bounding surfaces which was addressed at various levels of approximation in @xcite .    without loss of generality we can assume here that @xmath23 and @xmath24 it will be helpful for our later developments to introduce an _ asymmetry parameter _ @xmath25 that will allow us to quantify the dissimilarity between the two bounding surfaces as @xmath26 furthermore , by suitably normalizing the results one can concentrate exclusively on the interval @xmath27 . all other cases can be mapped onto this interval with appropriate rescaling of the parameters . the values @xmath28 and @xmath29 represent exceptional points in the parameter space : @xmath28 is the standard symmetric case ( @xmath30 ) already amply treated in the literature , and @xmath29 represents the antisymmetric case ( @xmath31 ) with no counterions between surfaces that reduces to the trivial case of a planar capacitor . these two well - understood limiting cases will be thus omitted from our discussion . counterions between the surfaces satisfy electroneutrality condition that can be written in the form @xmath32 where @xmath33 denotes the ( infinite ) area of each surface . because of the asymmetry present in the system , we have two length scales describing the interaction of the counterions with each of the bounding surfaces . these two length scales are given by the corresponding gouy - chapman lengths associated with the two surfaces as @xmath34 for the same reason we can thus define two different coupling parameters @xmath35 each one being defined by the ratio between the bjerrum length and the corresponding gouy - chapman length . in what follows , we rescale the surface separation as @xmath36 ( or the rescaled half - distance as @xmath37 ) with respect to plate 1 . with an appropriate rescaling one could also equivalently define all dimensionless lengths with respect to plate 2 . thus the minimal set of dimensionless parameters that fully characterize the system in the thermodynamic limit is given by @xmath38 . other physical quantities such as the mean electrostatic potential @xmath39 , number density of counterions , @xmath40 , and the pressure , @xmath41 , acting on each surface can be rescaled as well . we shall use the standard rescaled electrostatic potential @xmath42as well as the rescaled density and pressure @xmath43 where @xmath44 , and all the other quantities have been defined above . in the weak coupling regime , the leading contribution to the partition function comes from the saddle - point configuration of the local fluctuating electrostatic potential , @xmath45 @xcite . the saddle - point configuration can be straightforwardly translated into a solution of the pb equation and corresponds to an exact asymptotic result in the limit @xmath14 @xcite . for a system containing only counterions , the pb equation for the dimensionless potential , @xmath46 , can be written in the standard form @xcite @xmath47 with boundary conditions @xmath48 integration of the pb equation gives rise to the first integral of the system of the form @xmath49 where the constant @xmath50 is nothing but the mean - field pb pressure acting between the bounding surfaces @xcite and @xmath51 is the pb number density profile of counterions between the surfaces . the normalization factor @xmath52 follows from the electroneutrality condition ( [ eq : electroneutral ] ) as @xmath53    the nature of the solution @xmath54 obviously crucially depends on the sign of the pressure @xmath50 @xcite . different forms are obtained for positive and negative pressures , corresponding to repulsion and attraction between the bounding surfaces respectively . we review these different cases separately .      in the case of repulsive pressure the appropriate solution of eq . ( [ eq : pbeq ] ) can be written as @xmath56 where the constants @xmath57 and @xmath58 are obtained from the boundary conditions ( [ bc - rep ] ) and satisfy the set of two equations @xmath59 eliminating @xmath57 we obtain an equation for @xmath60 of the form @xmath61 the solution of this equation provides the final result for the repulsive pb pressure @xmath50 . note in particular that in rescaled units and by definition one has @xmath62 where @xmath63 . once @xmath60 is known , the parameter @xmath57 can be simply obtained from eqs . ( [ z01 ] ) or ( [ z02 ] ) and thus the potential @xmath64 , eq . ( [ psi0 ] ) , is fully determined . the density profile of counterions then follows from @xmath65    a positive ( repulsive ) solution for the pressure as considered in this section is always possible for any given asymmetry parameter @xmath25 ( excluding the trivial case of @xmath66 ) . in particular , it easily follows that within the mean - field theory two surfaces of equal sign ( @xmath67 ) always repel , that is at all separation distances @xmath68 . when the surfaces bear charges of opposite sign @xmath69 , they attract at large separations ( see below ) and a repulsion emerges only at sufficiently small separations .    at small separations @xmath70 , we can obtain the limiting solution for @xmath60 and thus the limiting small - distance pressure @xmath50 as @xmath71 for arbitrary @xmath72 as noted in scetion [ sec : geo ] . this is of course nothing but the ideal - gas osmotic pressure of counterion confined between the two plates ( _ i.e. _ , @xmath73 in actual units ) , which dominates over the energetic contributions at small separations . at large separations @xmath74 and for @xmath67 , we obtain the asymptotic expansion @xmath75 , \label{rep - limit - zeta}\ ] ] which is valid for @xmath76 . thus in the limit @xmath77 , the pressure behaves as @xmath78 which agrees with the asymptotic pressure between two equally charged surfaces ( @xmath79 ) . for smaller @xmath25 than determined above , _ i.e. _ for @xmath80 , one needs to invoke a different asymptotic expansion and specifically for @xmath81 ( one surface being neutral ) , one obtains @xmath82 this latter asymptotic result may be obtained from the one in eq . ( [ rep - limit ] ) by redefining @xmath83 . this may be understood simply by noting that because of symmetry a system with @xmath67 may be decomposed into two halves each with an effective asymmetry parameter @xmath84 . an attractive pressure on the mean - field level is possible only if the surfaces are oppositely charged @xmath69 . the appropriate solution in this case is given by @xmath86 where the constants @xmath57 and @xmath87 can again be obtained from boundary conditions , this time in the form @xmath88 eliminating @xmath57 we obtain an equation for @xmath60 as @xmath89 in this case we have in rescaled units @xmath90 and for the density profile of counterions @xmath91    the asymptotic form of the attractive pressure at large separations @xmath92 can be derived as @xmath93 where @xmath69 as noted above . for infinite separations , this pressure does _ not _ vanish and exponentially approaches @xmath94 since for @xmath95 the system behaves partially as a simple capacitor .      in the case of charged surfaces with opposite sign ( @xmath69 ) , the large - distance attraction regime and the short - distance repulsion regime merge at the point of zero pressure , @xmath97 , where the surfaces are at equilibrium . in this case , the pb solution for the potential reads @xmath98 and the density profile of counterions is given by @xmath99 where @xmath57 is found from the boundary conditions as @xmath100 , and the _ bound - state separation _ , @xmath101 , follows in rescaled units as @xmath102 the surfaces attract for @xmath103 and repel for @xmath104 . in the vicinity of @xmath105 , that is for @xmath106 , the pressure behaves as @xmath107 this concludes the calculation of the inter - surface pressure on the mean - field pb level strictly valid for @xmath14 . the preceding results may be summarized in a phase diagram shown in fig . [ fig : bound - sc ] in terms of @xmath101 and the asymmetry parameter @xmath25 displaying the mean - field attraction and repulsion regimes separated by the boundary line ( [ eq : pb_a * ] ) . the forms of the pressure here are completely consistent with those derived by lau and pincus @xcite _ via _ a different route . the first non - zero correction to the saddle point is second order in the fluctuations of the local electrostatic potential around the mean - field pb solution , @xmath108 . our goal here is to calculate the corrections in pressure , @xmath109 , stemming from these quadratic fluctuations , which leads then to the total wc pressure @xmath110 this approach has correlates in many diverse areas of physics where fluctuations around a mean - field solution are important @xcite and goes under different names , though the physics is always the same . we may conventionally refer to the mean - field pb term , @xmath111 , and the fluctuations contribution , @xmath109 , as the _ zeroth - order _ and the _ second - order _ correction terms on the wc level , respectively . this procedure formally also corresponds to a series expansion in powers of @xmath0 ( loop expansion ) around the asymptotic mean - field solution ( @xmath14 ) @xcite and is thus expected to be valid for sufficiently small coupling parameters as will be determined later . note also that in this latter sense the second - order pressure turns out to be proportional to @xmath0 , that is @xmath112 , and thus corresponds to a first - loop correction @xcite .    in order to proceed , one needs to evaluate the appropriate hessian of the field action in the partition function and study its fluctuation spectrum ( see refs . @xcite for more details ) . the hessian of the field action can be derived in the form @xmath113 where @xmath114 is the inverse coulomb operator and @xmath115 is the zeroth - order pb density as derived in the previous section . hence , @xmath116    the corresponding correction , @xmath117 , to the free energy of the system is then given by the trace - log of the hessian . it can be written equivalently in the following form @xcite @xmath118 this form can be derived rather straightforwardly by using the argument principle @xcite and converting the discrete sum of eigenvalues of the hessian operator into an integral over the transverse wave - vector @xmath119 , with density of modes @xmath120 of the logarithm of the secular determinant @xmath121 of the same operator . the index @xmath122 in the secular determinant refers to the eigenvalue equation that can be derived in the form @xmath123 by simply writing @xmath124 for the quotient @xmath125 , and noting that the secular determinant depends explicitly also on the value of the inter - surface spacing , @xmath20 , the free energy contribution from the quadratic fluctuations can be equivalently expressed exactly in a dimensionless form as @xmath126 where @xmath127 , @xmath128 and the rescaled area @xmath129 . here the secular determinant of the hessian for homogeneous transverse modes has been written as a function of dimensionless quantities @xmath130 , @xmath128 . this determinant has to be standardly regularized so that all irrelevant constants , _ i.e. _ all the terms not depending on the separation between the bounding surfaces , are dropped , amounting to a rescaling @xmath131 this corresponds to a subtraction of the part of the free energy for two separate interfaces at infinite separation from the total free energy .    in the next step one has to calculate the secular determinant @xmath132 for each of the pressure regimes separately , since the appropriate eigenfunctions of the hessian depend on the mean - field solution that in its turn depends on the sign of the interaction pressure , see eq . ( [ sign - eigen ] ) . in what follows we shall follow closely the derivations in refs . @xcite . the total pressure in the weak - coupling limit is thus the sum of the pb pressure and the quadratic fluctuations correction and can be written as @xmath133      in this regime the secular determinant of the hessian , eq . ( [ eigf ] ) , can be obtained by solving @xmath134 with appropriate boundary conditions implying continuity of the solution and its derivative across the bounding surfaces at @xmath18 . the general solution of eq . ( [ sec - rep ] ) for various regions in the perpendicular direction can be written in the form @xmath135 where @xmath136,\\      & & y_{2}=\cfrac{q^2}{q^2+\alpha^2}\>e^{-qz}\left[1- \cfrac{\alpha}{q}\,\tg\,\alpha ( z - z_0)\right].\end{aligned}\ ] ] taking into account the continuity of the solution and its derivatives we get a set of four homogeneous equations for the coefficients @xmath137 and @xmath138 . the solution exists only if the ( secular ) determinant of this system equals zero . thus we derive the secular determinant of the hessian operator in this case in the rescaled form @xmath139 where @xmath140 the regularized form of the secular determinant is obtained by taking the quotient as indicated in eq . ( [ renorm ] ) . while doing this , it is important to realize that @xmath141 also depends on the inter - surface distance . in fact from eq . ( [ rep - limit ] ) it follows that the appropriate limit of @xmath141 is @xmath142 in the regularization of the secular determinant this limiting behavior should be consistently taken into account .    finally , the dimensionless quadratic fluctuations free energy @xmath143 can be calculated numerically _ via _ eqs . ( [ second - f ] ) and ( [ eq : seculard_final ] ) . the fluctuations contribution to the pressure , @xmath144 , then follows from eq . ( [ secondord ] ) . the asymptotic form of the second - order dimensionless pressure can be obtained analytically . note that at large separations @xmath145 , a repulsive mean - field pressure @xmath55 , as considered in this section , is possible only for non - negative @xmath25 . for not too small @xmath146 , _ i.e. _ when @xmath147 , we find @xmath148 while for @xmath81 , _ i.e. _ when @xmath80 , we get @xmath149 again the difference in the two cases above is due to the symmetry of the problem in the latter case , that can be described by redefining @xmath83 and discarding the sub - dominant terms . the second - order pressure is obviously _ attractive _ and in this regime leads to a reduction of the total pressure from the mean - field value @xmath50 . this clearly shows that electrostatic correlations favor attraction between two repelling asymmetrically charged plates . however , the total pressure never becomes negative as the fluctuations are assumed to be small within the second - order weak - coupling analysis .      in this case the secular determinant of the hessian , eq . ( [ eigf ] ) , is obtained by solving @xmath150 the general solutions for particular regions in the @xmath151 direction are @xmath152 where @xmath153,\\      & & y_{2}=\cfrac{q^2}{q^2-\alpha^2}\,e^{-qz}\left[1+\cfrac{\alpha}{q}\,\cth\,\alpha ( z - z_0)\right].\end{aligned}\ ] ] again the solution exists only if the determinant of the system of equations , which connect coefficients @xmath154 and stems from the application of the boundary conditions at @xmath155 , is identically zero . this again defines the secular determinant @xmath124 appropriate for this case . it is easy to show that the secular determinant @xmath124 can be obtained from the @xmath55 result , eq . ( [ eq : seculard_final ] ) , simply by substituting @xmath156 and so we can write @xmath157 where @xmath158 here we can again regularize the secular determinant to discard divergences , eq . ( [ renorm ] ) . again one has to be careful by taking the correct limit for @xmath141 in the above regularization scheme . in this regime , the appropriate limit is given by @xmath159 as follows straightforwardly from eq . ( [ att - limit ] ) . the fluctuations contribution to the pressure , @xmath144 , can then be evaluated numerically from eq . ( [ secondord ] ) . the asymptotic form of @xmath144 for @xmath74 can be derived analytically as @xmath160 which is applicable only for charged surfaces of opposite sign , @xmath161 , which can attract ( @xmath85 ) at large separations . the function @xmath162 is defined as @xmath163 for @xmath164 , and @xmath165 for @xmath166 . the second - order pressure thus asymptotically decays exponentially and can be only attractive . for not too small values of @xmath25 , it is thus qualitatively very different from the case @xmath167 .    the total weak - coupling pressure , @xmath168 , is shown in fig . [ fig : wc_p ] for a few different asymmetry parameters and a relatively small value of the coupling parameter @xmath0 . as noted above the foregoing weak - coupling analysis is valid as long as the quadratic corrections are sufficiently small so that the series expansion around the mean - field solution does not diverge @xcite . as an approximate measure for the validity regime of this scheme , one can require that the second - order correction term is smaller than the leading order term , _ i.e. _ @xmath169 this leads to a useful criterion identifying the regime of coupling parameters and distances in which the weak - coupling theory is applicable . for @xmath55 and by employing the closed - form expressions obtained for large separations @xmath92 , we find the validity condition @xmath170 this indicates that at a given non - vanishing @xmath0 , the weak - coupling scheme becomes increasingly more accurate at larger separations , while as the surfaces get closer a smaller coupling parameter needs to be chosen . on the other hand , for @xmath85 ( which occurs for @xmath161 ) and at large separations @xmath92 , we obtain @xmath171 the right hand side here is exponentially large meaning that for charged surfaces of opposite sign , the weak - coupling analysis performs far better at finite coupling parameters and smaller inter - surface separations than for the surfaces of equal sign ( @xmath67 ) . finally , note that for @xmath96 that corresponds to the equilibrium phase boundary line in fig . [ fig : bound - sc ] for @xmath161 , we deal with a situation where the leading order term is zero and the fluctuations are _ dominant _ at any finite value of @xmath0 . the convergence of the loop expansion has to be determined in this case by evaluating the higher order terms which we shall not consider in this paper . the strong - coupling approximation coincides with the lowest order non - trivial expansion of the partition function in terms of the fugacities of the counterions . this expansion may be expressed as a @xmath16 series expansion @xcite , whose leading order term ( @xmath17 ) corresponds to the so - called sc theory . we will not delve into the strong - coupling expansion in more detail since it has been exhaustively reviewed in the literature @xcite . on the leading order , the free energy is obtained as @xmath172 where @xmath173 is electrostatic interaction energy of charged surfaces @xmath174 with @xmath33 representing surface area of each plate , and @xmath175 and @xmath176 are electrostatic interaction energies between a single counterion and individual charged surfaces , _ i.e. _ @xmath177 since in the strong - coupling regime the free energy is given _ via _ simple quadratures , it is much simpler to evaluate it than on the weak - coupling level . defining the rescaled free energy @xmath178 we obtain @xmath179.\ ] ] differentiating the free energy with respect to the surface - surface distance @xmath180 we get the corresponding pressure acting between the bounding surfaces @xmath181 . \label{eq : p_sc}\ ] ] the dependence of this dimensionless pressure on the separation for different values of @xmath25 is presented in fig . [ fig : sc_p ] . note that the sc pressure can become attractive for _ both _ like - charged and oppositely charged surfaces which contrasts with the mean - field theory that does not allow attraction between like - charged surfaces . this is because of the strong electrostatic correlations mediated by counterions between the charged surfaces for @xmath182 and has been investigated throughly before for equally charged surfaces @xcite . our results show that a similar attraction mechanism holds for asymmetrically charged surfaces in the sc limit . the pressure exhibits two well - defined limiting laws obtainable for small and large inter - surface separations . for small separations @xmath70 , we have @xmath183 for arbitrary @xmath72 as noted in section [ sec : geo ] . for large separations @xmath74 , we obtain @xmath184 note that in the limit @xmath185 the sc pressure coincides with the pb result on the leading order ( compare eqs . ( [ rep - limit1 ] ) and ( [ psc_limit ] ) ) and represents the ideal - gas osmotic pressure of counterions which dominates over the electrostatic contributions . the pb and sc forms for @xmath161 coincide also in the limit of @xmath74 ( compare eqs . ( [ att - limit ] ) and ( [ newlimit ] ) ) and reduce to the pressure in the capacitor . the dependence of the pressure on the inter - surface separation points to the existence of a bound state defined _ via _ @xmath186 . the sc bound - state separation @xmath101 can be expressed analytically as @xmath187 and is presented in fig . [ fig : bound - sc ] as a function of @xmath25 . obviously both like - charged and oppositely charged surfaces can form bound states at small surface - surface separations . the bound - state separation approaches infinity and the surfaces _ unbind _ asymptotically as @xmath188 , that is when one plate becomes electroneutral . finally one can also derive the explicit form of the counterion density as a function of the normal coordinate @xmath151 . this can be read off simply from the integrand in eq . ( [ deff ] ) , that is @xmath189 , where @xmath190 is a normalization factor @xcite . according to the electroneutrality condition , we normalize the density as @xmath191 , or in rescaled units @xmath192 where we have defined @xmath193 with @xmath194 being the gouy - chapman length with respect to plate 1 . from here , the density profile is obtained as @xmath195 } , \label{denssc}\ ] ] as a function of @xmath151 and @xmath25 . the regime of applicability of the leading order sc theory follows from a simple criterion that has been discussed and confirmed previously in the case of equally charged surfaces by both mc simulations and higher - order calculations @xcite . the generalization to asymmetrically charged surfaces is straightforward .    for large couplings , counterions are strongly attracted to an oppositely charged surface as the counterion - surface interaction becomes large and equivalently , the gouy - chapman length , @xmath10 , and the thickness of the counterionic layer at the surface become small . the layer thickness has to be compared with the typical lateral spacing between counterions @xmath196 . for counterions sandwiched between two asymmetrically charged surfaces , this latter quantity follows from the local electroneutrality condition as @xmath197 up to a factor of the order unity and assuming that the surfaces are sufficiently close so that they may be strongly coupled _ via _ the counterions as will be determined consistently here . in rescaled units , one gets @xmath198 where @xmath199 . obviously , @xmath200 becomes large relative to the layer thickness as @xmath0 grows . note that on the other hand the lateral coulomb repulsion between counterions in this quasi - two - dimensional layer becomes much larger than the thermal energy , _ i.e. _ @xmath201 , indicating that counterions form a strongly correlated liquid in which they are highly localized within correlation holes of lateral size @xmath202 @xcite . thus for surface - surface separations , @xmath68 , smaller than the correlation hole size , _ i.e. _ @xmath203 counterions can move almost independently from each other in the direction normal to the surface and one can safely assume that the effective surface - surface interaction as well as the counterionic density profile follow only from the interactions of _ individual _ counterions with the bounding charged surfaces . the counterion - counterion interactions contribute on the sub - leading order and matter at larger separations . this picture is of course confirmed on a systematic level by the sc expansion analysis @xcite , and the above equation sets a criterion for the validity regime of the single - particle leading order sc theory ( @xmath204 ) , when applied to finite coupling parameters . we performed monte - carlo ( mc ) simulations in order to study the system of two asymmetrically charged surfaces beyond the analytical limits of weak and strong coupling discussed above . all simulations were performed in the canonical ensemble ( nvt ) using the standard metropolis algorithm @xcite . the mobile counterions were modeled as point charges @xcite enclosed in a simulation box bounded in @xmath151 direction by two charged surfaces of distance @xmath138 and surface charge densities @xmath205 and @xmath206 ( compare fig . [ fig : geometry ] ) . periodic boundary conditions were applied in the lateral directions parallel to the bounding surfaces . the lateral size of the charged surfaces , @xmath207 , which is set equal to the lateral size of the simulation box , was held fixed throughout all the simulations . the number of counterions , @xmath208 , was varied between 60 and 1800 , depending on the system parameters in order to fulfill the electroneutrality condition . the counterions interact through the coulombic potential @xmath209 with @xmath210 being the charge of the @xmath211-th counterion and @xmath212 the separation distance between the @xmath211-th and the @xmath213-th counterions ( with @xmath214 ) . the interaction energy of the @xmath211-th counterion with the charged surfaces in the simulation box are given by @xmath215 ^ 2}\bigg )          + \frac{\pi}{2}\bigg\}\bigg ] , \label{coulomb_wallterm}\end{aligned}\ ] ] where @xmath216 and @xmath206 is surface charge density of the @xmath217-th surface ( with @xmath218 ) and @xmath219 is the normal distance between the @xmath211-th counterion and the @xmath217-th surface . the long - ranged coulomb interactions in this system were accounted for _ via _ a charged sheet scheme similar to that proposed by torrie and valleau @xcite . this scheme makes use of the counterion profile in the simulation box in order to calculate an external field , stemming from the long - ranged interactions . this external field is iteratively updated and self - consistency is achieved normally in a few iterations . in the course of simulations , new configurations were created by trial displacements of the counterions and equilibration was accomplished by running through @xmath220 configurations . the @xmath221 following configurations were then used for the production runs . the pressure @xmath41 was calculated in the production runs according to the contact - value theorem as @xmath222 where @xmath223 is the density of counterions at contact with the @xmath217-th surface . in thermodynamic equilibrium , the pressure does not depend on which surface ( @xmath224 or @xmath225 ) is chosen in order to calculate the pressure from the above equation , and the contact condition at both surfaces leads to precisely the same value for the pressure . all simulations were conducted at fixed temperature @xmath226  k , lateral simulation box size @xmath227   and dielectric constant @xmath228 , which is assumed to be the same throughout the system . the simulations were performed at different values of the coupling parameter @xmath0 and the asymmetry parameter @xmath25 . we explored the @xmath0 parameter space extensively by using @xmath229 , @xmath230 and @xmath231 in order to cover exhaustively both the weak coupling and the strong coupling regimes . the concurrent values of the asymmetry parameter were always taken as @xmath232 at each value of the coupling parameter . the results are plotted in the form of rescaled density and pressure as previously defined in this paper ( section [ sec : rescal ] ) . in order to asses the validity of the weak and strong coupling results presented above for asymmetrically charged surfaces , we performed extensive mc simulations and compared them to analytical results in both limits . it transpires from this comparison that the simulation results corresponding to an exact evaluation of the partition function are always bracketed by the wc and the sc limiting forms , smoothly approaching them in the appropriate limits of the coupling parameter @xmath0 . first we compare the density profiles of simulations with theoretical results given by eqs . ( [ denswc1 ] ) and ( [ denswc2 ] ) for the pb limit ( solid lines ) and eq . ( [ denssc ] ) for the sc limit ( dashed lines ) . as seen from fig . [ dens ] the theoretical pb and sc rescaled density profiles represent two extremal cases and all mc simulations results with finite values of @xmath0 are located consistently between these two limits . the mc results for small @xmath0 are almost exactly spot on the pb prediction , while larger discrepancies are observed as @xmath0 grows . for large enough @xmath233 , the mc results slowly converge to the sc result . this is especially clear for surfaces with charges of equal sign , @xmath234 , whereas for surfaces with opposite sign , @xmath235 , there is no big difference between pb and sc profiles . next we consider the inter - surface pressure as obtained from the simulations ( symbols in fig . [ fig : sim_pressure ] ) as well as the pb theory , eqs . ( [ pres - pb - rep ] ) and ( [ pres - pb - att ] ) , and the sc theory , eq . ( [ eq : p_sc ] ) ( solid and dashed lines , respectively ) . the pb result is expected to be valid for separations @xmath236 @xcite ( see also section [ sec : wc_crit ] ) . therefore , for @xmath237 the pb line expectedly agrees nicely with the simulation data ( open squares ) in the whole range of separations shown in the figure . upon closer inspection , however , we find small deviations from the pb result as shown in the insets in fig . [ fig : sim_pressure ] for all three values of the asymmetry parameter @xmath25 . in this case , the fluctuation correction to the mean - field pressure accurately compensates for these deviations and adding the second - order correction to the pb pressure leads to a total wc pressure , eq . ( [ secondord ] ) ( shown as a dotted line in the inset ) that matches the simulation data perfectly . note that the pressure changes can be drastic even on the wc level as @xmath25 assumes different values . for @xmath238 and 0 the wc pressure is strictly repulsive , while it turns attractive and leads to a bound state ( zero pressure point ) for @xmath239 .    in the intermediate regime of coupling parameters , the simulation results for the pressure are clearly bracketed by the two limiting analytical forms , given by the pb plus the second order correction and the sc expressions of the interaction pressure . the sc prediction is expected to be valid for separations @xmath240 as discussed in section [ sec : sc_crit ] . consistently , the interaction pressure starts off close to the strong coupling limit at small separations and then smoothly converges to the weak coupling limit for larger separations . this is strictly true for @xmath238 and 0 . in the case of @xmath241 the difference between the strong and weak coupling results for the rescaled interaction pressure is marginal and the simulation data and the analytical results nearly coincide for all rescaled separations @xmath242 . we emphasize that the pressures and the density profiles are plotted here in _ rescaled representation _ ; in actual units , fig . [ fig : sim_pressure ] corresponds to different ranges of separation , @xmath138 , for the wc and sc regimes as the gouy - chapman length , @xmath10 , is typically very different between the two limits ( small at high couplings and large at small couplings as may be realized , _ e.g. _ , by changing the counterion valency at fixed surface charge densities and bjerrum length ) .    for large values of the coupling parameter the simulation results for the interaction pressure expectedly follow very closely the strong coupling prediction for a wider range of inter - surface separations . the correspondence between the sc theory and simulations is better for @xmath241 than for @xmath238 , which can be again traced back to the fact that the strong and the weak coupling results are very close to one another for charged surfaces of opposite sign in the whole range of rescaled separations , whereas they differ significantly in the case of surfaces of equal sign . for intermediate and large couplings , we have not attempted to compare our data with the second - order wc approximation as this approximation breaks down at the range of distances shown in the figures ( section [ sec : wc_crit ] ) . note also that in all cases considered here the theoretical and the simulated values of the interaction pressure converge for very small surface - surface separations . in fact , the sc and pb results coincide in the leading order as the rescaled distance , @xmath68 , tends to zero , eqs . ( [ rep - limit1 ] ) and ( [ psc_limit ] ) , since both are dominated by the osmotic pressure of counterions . the sub - leading corrections for very small @xmath68 are different in the pb and sc limits and on this level the simulation data with finite @xmath0 are generally expected to agree better with the sc prediction at small separations @xcite . we now consider the simulated bound - state separation , @xmath101 , as a function of the coupling parameter in fig . [ fig : bound_xi ] ( symbols ) . for @xmath234 , the pb theory ( @xmath14 ) gives only repulsion and predicts no bound state . while the sc theory predicts a closely packed bound state with @xmath243 given by eq . ( [ eq : d*_sc ] ) , explicitly here @xmath244 . as seen , by increasing the coupling parameter the simulation data for @xmath101 ( filled squares ) decrease monotonically and rather slowly converge to the sc prediction ( dotted line ) . the comparison is again worse for the @xmath238 than for the @xmath241 case ( dashed line and filled circles ) . for @xmath234 , even at @xmath245 the difference between the simulations and the analytical result is still close to 10% . the deviations are quite pronounced for smaller values of the coupling parameter . the opposite is true for @xmath241 . here , the simulation results are close to the strong - coupling analytical limit in the whole range of @xmath0 values . for @xmath246 we can discern , see inset in fig . [ fig : bound_xi ] , weak coupling behavior that starts off with the pb - predicted value @xmath247 in the limit of @xmath248 ( eq . ( [ eq : pb_a * ] ) , solid line ) , that later follows the pb plus second - order corrections line ( dot - dashed line ) and then rapidly approaches the strong - coupling result @xmath249 ( dashed line ) . as already noted the differences between strong and weak coupling in this case are marginal in the rescaled representation .    the dependence of @xmath105 on @xmath25 in fig . [ fig : bound_zeta ] complements the above observations . the bound - state separation diverges for @xmath250 both in simulations as well as in the analytical limits . here we reproduce the simulation results at @xmath234 and @xmath251 , which again clearly show the convergence to the sc result for @xmath234 and the crossover from the pb result ( solid line ) to the sc result ( dashed line ) for @xmath235 upon increasing the coupling parameter . the quantitative agreement between the simulations and the analytical results in the two limiting cases of pb and sc is excellent . note here again that the second - order wc expansion ( dotted line in the inset ) begins to fail at small separations as it deviates from the limiting pb results as well as from the simulation data . since the pb pressure is zero on the pb phase boundary line , the validity of the second order wc correction can be assessed analytically only by performing a two - loop calculation which goes beyond the scope of this paper . to summarize , we have derived theoretical forms for the interaction pressure as well as the counterionic density profiles of asymmetrically charged planar surfaces with neutralizing counterions in between . based on the field - theoretical methods we analyzed two different regimes of weak and strong coupling as defined by the electrostatic coupling parameter @xmath0 . the crossover between these two regimes is studied _ via _ monte - carlo simulations .    for small values of @xmath0 , the system is described very well by the _ weak coupling _ ( wc ) theory , that in the lowest ( zeroth ) order coincides with the mean - field poisson - boltzmann ( pb ) result . the second order of wc corresponds to a first - order loop expansion and represents the contribution from correlated quadratic fluctuations around the mean - field or saddle - point solution . this second - order correction , which is proportional to @xmath0 , always lowers the interaction pressure between the surfaces and thus leads to an attractive contribution to the total interaction pressure . since it corresponds to an expansion of the partition function around the mean - field saddle point it has to be smaller , in absolute terms , than the mean - field result . net attraction given by second - order fluctuations term ( for interacting surfaces of equal sign ) is therefore inconsistent with the nature of the wc approximation .    for large values of the coupling parameter @xmath0 , the weak coupling approach breaks down and the virial expansion amounting to the _ strong coupling _ ( sc ) approximation must be used . the sc approach is effectively a one - particle theory and takes properly into account the strong correlation and interaction of the counterions with external surface fields on the leading order @xcite . following standard procedures , we derived an analytical expression for interaction pressure in the sc limit . the interaction pressure in this case is always lower than the pb result and can be negative ( corresponding to a net attractive force ) even for charged surfaces of equal sign . we compared both our theories , _ i.e. _ wc with second - order corrections and sc , with monte - carlo simulations . we found very good agreement for both theories in their expected regime of validity . as expected , the wc approach is valid for separations @xmath252 . the second - order correction improves the small discrepancies between pb and mc results at large separations but it tends to fail for smaller distances when discrepancies get more pronounced . on the other hand , the sc theory describes the behavior perfectly at small separations @xmath253 . for small enough coupling parameter @xmath0 , the validity of the pb approximation spreads to smaller separations , @xmath68 , where pb and sc results nearly coincide ( in the rescaled representation ) . therefore , we may conclude that for sufficiently small @xmath0 , the pb result is valid on the whole interval @xmath68 . note that the second - order wc correction term consistently diverges ( toward large negative values ) for small inter - surface separations , @xmath254 , irrespective of @xmath25 , which makes it in general inapplicable in this limit . the reason for this is simple . for small inter - surface separations the mean - field solution becomes more and more homogeneous , almost a constant , and the interaction free energy approaches its standard zero - frequency van der waals form that diverges for small separations . in the case of surface charges with equal sign , @xmath67 , the wc theory predicts no attraction and hence no bound state . the attraction and the corresponding bound state appear only for coupling parameters @xmath0 that are large enough , as predicted by the sc theory . in the case of charged surfaces of opposite sign , @xmath161 , the attraction appears also in wc limit above a threshold value @xmath255 that represents the equilibrium surface - surface separation . it is notable that for charged surfaces of opposite sign , the wc analysis in general performs much better than for the surfaces of equal sign and that the sc and the wc results are very close to one another for charged surfaces of opposite sign in the whole range of rescaled separations , whereas they differ significantly in the case of surfaces of equal sign . there is also only a marginal difference between the pb and sc counterion density profiles in the rescaled representation for surfaces of opposite sign . a reasonable explanation for this would be in our opinion that for oppositely charged surfaces the counterions mostly feel the effect of the strong uniform external field provided by the surface charges , which acts similarly in the strong as well as the weak coupling limit . thus the mean - field and the strong - coupling approaches should converge . in the case of similarly charged surfaces , the mean - field theory depends more on the local counterion density whereas the strongly coupled counterions still feel mostly the external field . thus the difference between the wc and the sc frameworks in the @xmath161 and @xmath67 cases . our results support an emerging new paradigm , according to which the wc and the sc limit bracket the exact results for the interaction pressure between charged surfaces neutralized by mobile counterions . they indeed provide quantitatively correct results for the interaction pressure in the limit of small and large inter - surface separations , while at intermeditae separations the exact results are always located between the two limits . it thus seems advisable that in analyzing the electrostatic interactions in colloidal systems one always calculates both analytic limits , the wc as well as the sc , in order to get a good handle on the range of values that the interaction can assume for any value of the electrostatic coupling parameter . in future we intend to study the same system in the presence of added salt and dielectric discontinuities . m.k . and r.p . would like to acknowledge the financial support by the agency for research and development of slovenia ( grants p1 - 0055(c ) , z1 - 7171 , l2 - 7080 ) . this study was supported by the intramural research program of the nih , national institute of child health and human development . this research was supported in part by the national science foundation under grant no . phy05 - 51164 .                                                                    in our model excluded - volume counterion - counterion repulsions are neglected . these two - body interactions will be irrelevant in the leading order within the strong - coupling limit @xcite . however , the counterion size ( or the counterion - wall exluded - volume interactions ) is incorporated exactly by taking @xmath256 as the closest - approach surface - surface distance leaving out the counterion diameter from the actual distance . on a more fundamental level it is indeed possible to incorporate excluded- volume effects within a modifed poisson - boltzmann formulation as in _ e.g. _ : i. borukhov , d. andelman , h. orland , phys . lett . * 79 * 435 ( 1997 ) ; electrochimica acta * 46 * 221 ( 2000 ) . r. r. netz , eur . j. e * 3 * 131 ( 2000 ) ; r. menes , p. pincus , b. stein , phys . e * 62 * 2981 ( 2000 ) ; j. schmit , r. menes , p. a. pincus , phys e * 66 * 061502 ( 2002 ) ; m. m. hatlo , r. a. curtis , l. lue , j. chem . phys . * 128 * 164717 ( 2008 ) .
Abstract:

Answer: we compare weak and strong coupling theory of counterion - mediated electrostatic interactions between two asymmetrically charged plates with extensive monte - carlo simulations . analytical results in both weak and strong coupling limits compare excellently with simulations in their respective regimes of validity . the system shows a surprisingly rich structure in terms of interactions between the surfaces as well as fundamental qualitative differences in behavior in the weak and the strong coupling limits .