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: in view of various applications such as the study of two  dimensional melting @xcite , investigations of mesoscale structure formation @xcite or engineering of colloidal crystals on spherical surfaces @xcite , the self  assembly of sub-@xmath0 m colloidal particles at water  air or water  oil interfaces has gained much interest in recent years . these particles are trapped at the interface if the colloid is only partially wetted by both the water and the oil . this configuration is stable against thermal fluctuations and appears to be the global equilibrium state , as it is observed experimentally that the colloids immersed in the bulk phases are attracted to the interface @xcite ( see also sec . [ sec:1coll ] ) . the mutual interaction between the trapped colloids at distances close to contact , i.e. , within the range of molecular forces , is dominated by strong van  der  waals attraction . in order to avoid coagulation due to this attraction , the colloids can be stabilized sterically with polymers or with charges such that the colloids repel each other . variants of charge stabilization may include the coverage with ionizable molecules which dissociate in water , or the labeling of colloids with charged fluorescent markers . for charge  stabilized colloids at large distances , the resulting repulsive force at water interfaces stems from a dipole  dipole interaction as shown theoretically for point charges on surfaces @xcite and verified experimentally for polystyrene ( ps ) spheres on water  oil interfaces @xcite .    nonetheless , charged colloids at interfaces also show attractions far beyond the range of van  der  waals forces . the corresponding experimental evidence can be roughly classified as follows . ( i ) according to refs . @xcite , ps spheres ( radii @xmath1 @xmath0 m ) on flat water  air interfaces using highly deionized water exhibit spontaneous formation of complicated metastable mesostructures . they are consistent with the presence of an attractive , secondary minimum in the intercolloidal potential at distances @xmath2 with a depth of a few @xmath3 . the use of water slightly contaminated by ions seems to move the minimum further out and to reduce its depth @xcite . ( ii ) in ref . @xcite , pmma spheres with radius @xmath4 @xmath0 m were trapped at the interface of water droplets immersed in oil . here , the secondary minimum has been measured at a distance @xmath5 and is reported to be surprisingly steep . the tentative explanation of these findings given in ref . @xcite invokes an analogue of long  ranged flotation or capillary forces which decay @xmath6 . this interpretation was criticized in ref . @xcite ( with which the authors of ref . @xcite agreed @xcite ) and in ref . @xcite which both concluded that possible capillary forces in this system are much shorter ranged , i.e. , @xmath7 , but the authors of these references disagree with respect to the sign of this shorter  ranged force . in yet another twist of the story , after completion of our work we encountered the very recent ref . @xcite in which the authors claim that long  ranged capillary forces @xmath6 caused by the colloidal charges persist for sub-@xmath0 m particles . this conclusion is based on measurements of the meniscus shape around single glass spheres with radii 200  300 @xmath0 m floating at water  oil and water  air interfaces .    motivated by the experimental data summarized above and the still incomplete theoretical understanding , here we undertake a quite general analysis of capillary interactions between two spherical colloids trapped at fluid interfaces . we characterize the system by a general stress field which acts on the interface , e.g. , due to a discontinuous electrostatic field at the interface , and by forces on the colloid of gravitational and/or electrostatic origin . special attention will be given to the role of a restoring force fixing the interface , e.g. , due to gravity or interface pinning . in sec . [ sec:1coll ] we present a free energy model for a single colloid trapped at an interface in the limit of small stresses and forces . the general solution for the interface deformation will be used in sec . [ sec:2coll ] in order to determine the effective potential between two colloids within the superposition approximation . in view of the differing theoretical results in the literature , the derivation of the interface deformation and the resulting effective potential is presented in detail in order to reveal properly the subtleties involved . in sec . [ sec : discussion ] we apply the results to the case of charged polymeric spheres on water interfaces . it turns out that the constraint of approximate mechanical isolation of the experimental systems renders the capillary interaction basically short  ranged . long  ranged attractive forces @xmath6 only arise through a restoring force acting on the interface ( which is , however , weak ) or in the presence of an external electric field . we shall discuss the relation to previous theoretical results , especially in view of the experimental results reported in refs . directions for further research will be pointed out . in sec . [ sec : summary ] we summarize our results . in this section we consider as a first step the equilibrium state of a single colloid of radius @xmath8 at the interface between two fluid phases denoted as 1 and 2 . the contact angle formed by the interface and the colloid surface is given by young s law , @xmath9 where @xmath10 is the surface tension between phases 1 and 2 , and @xmath11 is the surface tension between the colloid and phase @xmath12 . as a reference state , with respect to which changes in free energy will be measured , we take a planar meniscus configuration with the colloid at such a height @xmath13 that young s law is satisfied ( fig . [ fig : ref ] ) . this corresponds to the equilibrium configuration of an uncharged colloid at the interface if its weight can be neglected  which for generic cases is a safe approximation for @xmath14 @xmath0 m @xcite ( see also sec . [ sec : discussion ] ) . we model the colloid as a smooth sphere so that the system is invariant under rotations around the colloid axis perpendicular to the reference planar meniscus . the presence of charges induces a shift of the system ( colloid and interface ) with respect to the reference state . here we neglect corresponding changes in the surface tensions @xmath10 and @xmath11 ; this approximation is expected to be valid provided the concentration of charges is sufficiently small . this shift is characterized by the meniscus profile @xmath15 relative to the planar configuration and by the height @xmath16 of the colloid center . in the reference configuration , the charge distribution is assumed to be already in equilibrium . in the following we do not consider the degree of freedom `` charge density field '' explicitly but take it to be fixed to that of the reference configuration . this amounts to neglecting the feedback of the interface displacement on the charge distribution . it turns out to be useful to introduce the radius @xmath17 of the three  phase contact line and the angle @xmath18 as auxiliary variables ( see fig .  [ fig : pert ] ) .      in this subsection we formulate a free energy functional for the degrees of freedom @xmath15 and @xmath16 . as shown later , for the cases of interest here the deviations from the reference configuration are small enough to justify a perturbative treatment , and the free energy can be expanded up to quadratic order in @xmath19 and @xmath15 . we denote with @xmath20 the vertical force per unit area acting on the meniscus surface in the reference configuration . in the case of a charged colloid , @xmath20 is given by the @xmath21component of the difference of the maxwell stress tensor right above and below the meniscus , plus the pressure difference acting across the meniscus ( incluiding an imbalance in osmotic pressure due to the different concentration of ions just above and below the meniscus ) @xcite . we introduce the following dimensionless parameter to measure the relative strength of this force : @xmath22 where the integral extends over the flat meniscus @xmath23 ( the plane @xmath24 with a circular hole of radius @xmath25 ) . we assume @xmath26 with @xmath27 , so that the integral converges . in the same spirit we introduce the total vertical force @xmath28 acting on the charged colloid in the flat reference configuration , which includes gravity , the electrostatic force , and the total ( i.e. , hydrostatic and osmotic ) pressure exerted by the fluids . this leads to the definition @xmath29 these two dimensionless parameters will appear naturally in the course of the calculations . if @xmath30 and @xmath31 vanish , the reference configuration _ is _ the equilibrium state ; it is the global minimum of the free energy functional given in eqs . ( [ eq : fdef ] ) and ( [ eq : f ] ) below if @xmath32 , @xmath33 . the aforementioned perturbative expansion of the free energy can be rephrased as an expansion in terms of the small parameters @xmath30 and @xmath31 . the free energy @xmath34 of the colloid expressed relative to the aforementioned reference configuration consists of five terms : @xmath35 in the following we discuss each contribution :    * _ fluid contact of the colloid_. with @xmath36 denoting the surface area of the colloid which is in contact with phase @xmath12 , the surface free energy of the colloid due to its exposure to the phases 1 and 2 is @xmath37 in appendix  [ sec : fcont ] we express this contribution as a function of @xmath38 and @xmath15 . the final result , valid up to corrections of at least third order in @xmath30 or @xmath31 , reads @xmath39 ^ 2 + \pi \gamma ( r_0 ^ 2 - r_{0,\rm ref}^2 ) . \ ] ] * _ change of the meniscus area_. the free energy contribution due to variations in the meniscus area relative to the reference state reads @xmath40 where @xmath41 is the surface of the fluid interface projected onto the plane @xmath24 ( in which the reference interface is located ) , and @xmath42 is the gradient operator on the flat reference interface . for small slopes ( @xmath43 ) one obtains : @xmath44 since the @xmath45dependent term is of second order in @xmath45 , we have approximated the integration domain @xmath41 by @xmath46 ; the corrections are at least of third order in the small parameters @xmath30 or @xmath31 . the first term in this expression represents the change in the area of the meniscus which is cut out by the colloid , and in eq . ( [ eq : fdef ] ) it cancels the second term of eq . ( [ eq : fcont_iso ] ) . * _ volume forces on the fluids_. we consider the case that the only volume force acting on the fluid phases is gravity . the electrostatic forces are active only at the surfaces , where a net charge can accumulate . if the spherical colloid is replaced by a long cylinder of radius @xmath25 ( fig . [ fig : cylinder ] ) the change in gravitational potential energy relative to the reference state due to the displacements of volumes of the fluid phases can be calculated easily : @xmath47 assuming that the mass density of phase 2 is larger than the one of phase 1 , i.e. , @xmath48 ; @xmath49 denotes the capillary length @xmath50 in terms of the acceleration @xmath51 of gravity . for a spherical colloid @xmath52 in eq . ( [ eq : fvol ] ) has to be corrected to account for the specific dependence of the displaced fluid volume on the colloid height @xmath16 , yielding a slightly cumbersome expression @xcite . this correction is usually numerically small and we neglect it for reasons of simplicity . as it will be shown in subsec . [ sec : lambda ] this is justified because the results of interest here are insensitive to the precise form of @xmath53 in the limit @xmath54 ( @xmath55 mm , which is much larger than any other relevant length scale of the system ) . * _ force on the fluid interface_. the aforementioned surface force density @xmath20 acts on the fluid interface between phase 1 and phase 2 . the free energy change due to the ensuing displacements of the meniscus is @xmath56 here @xmath20 is the surface force in the _ reference _ configuration ( @xmath57 corresponds to a force pointing upward ) . changes in the force induced by meniscus deformations and colloidal displacements contribute terms of higher orders in @xmath30 or @xmath31 in eq . ( [ eq : finter ] ) . the replacement of @xmath58 by @xmath23 in eq . ( [ eq : finter ] ) introduces terms of higher order , too . * _ contribution from the colloid_. the free energy change due to a vertical displacement of the colloid is @xmath59 where @xmath60 is the vertical force on the colloid in the _ reference _ configuration . like for @xmath61 , changes of @xmath60 due to deviations from the reference configuration contribute to higher order terms .    in conclusion , by adding eqs . ( [ eq : fcont_iso ] ) , ( [ eq : fmen ] ) , ( [ eq : fvol ] ) , ( [ eq : finter ] ) , and ( [ eq : fcoll ] ) we obtain the following approximate expression for the total free energy , which is correct up to second order in @xmath30 or @xmath31 , and which is a function of @xmath38 and a functional of @xmath15 : @xmath62 +    \pi \gamma [ u_0 - \delta h ] ^2 -    f \delta h , \ ] ] where @xmath63 . the equilibrium configuration of the system minimizes the free energy expression  ( [ eq : f ] ) . the minimization proceeds in two stages . first we seek the minimum with respect to @xmath38 at fixed @xmath15 : @xmath64 where we have used the definition  ( [ eq : epsf ] ) . using repeatedly the definitions of @xmath16 and @xmath17 in terms of the auxiliary angle @xmath18 ( fig . [ fig : pert ] ) we can compute the change of the contact radius : @xmath65 plus corrections of second order in @xmath30 , @xmath31 . in the second step , we minimize eq . ( [ eq : f ] ) with respect to @xmath15 at fixed @xmath38 . this is a problem of variations with a free boundary condition at @xmath66 @xcite . variation with respect to @xmath67 yields a second  order ordinary differential equation , @xmath68 while variation with respect to @xmath69 provides a boundary condition : @xmath70 where the last equality follows from eq . ( [ eq : h_eq ] ) . the second boundary condition one has to impose on eq . ( [ eq : younglaplace ] ) follows from the requirement that the meniscus is asymptotically flat far from the colloid ( assuming that @xmath71 ) : @xmath72 equation  ( [ eq : younglaplace ] ) describes mechanical equilibrium of the interface such that the laplace pressure balances the forces acting on the interface . the boundary condition  ( [ eq : u0 ] ) expresses mechanical equilibrium of the colloidal particle : at the contact line ( i.e. , the circle with radius @xmath17 at @xmath73 ) the interface exerts a force onto the colloid which has a non  vanishing contribution only in @xmath74direction with the magnitude @xmath75 \approx 2 \pi \gamma r_{0,\rm ref } u'_0 $ ] . this contact line force is balanced by the total force @xmath60 . the solution of eqs . ( [ eq : younglaplace][eq : u_infty ] ) can be written in terms of the modified bessel functions of zeroth order : @xmath76 , \ ] ] where the integration constant @xmath77 is determined by the boundary condition  ( [ eq : u0 ] ) .      for typical values of the parameters , @xmath49 is of the order of millimeters and therefore much larger than any other length scale occurring for experiments with sub - micrometer colloids . in order to study the intermediate asymptotics ( @xmath79 ) of @xmath15 as given by eq . ( [ eq : solbessel ] ) , we insert the asymptotic expansions of the bessel functions @xcite as @xmath78 and retain those terms which do not vanish in this limit . assuming that @xmath80 decays sufficiently fast , eq . ( [ eq : solbessel ] ) reduces to @xmath81 where @xmath82 and @xmath83 is euler s constant . in eq . ( [ eq : sol1 ] ) , we have expressed the integration constant @xmath77 appearing in eq . ( [ eq : solbessel ] ) in terms of @xmath31 by using the boundary condition  ( [ eq : u0 ] ) . the first term in eq .  ( [ eq : sol1 ] ) is a solution of the homogeneous part of eq . ( [ eq : younglaplace ] ) ( with @xmath84 in the equation ) , demonstrating that the limit @xmath78 is singular as long as @xmath85 . the second term corresponds to a particular solution of the inhomogeneous differential equation . if the surface force @xmath20 decays algebraically , @xmath86 , this term decays like @xmath87 ( since we have assumed @xmath27 ) , so that the logarithmic contribution is dominant . if @xmath88 , the asymptotic behavior is no longer given by eq . ( [ eq : sol1 ] ) but a different dependence on @xmath49 arises .    at distances @xmath89 of the order of @xmath49 , the expression  ( [ eq : sol1 ] ) is not valid and a crossover to the exact solution  ( [ eq : solbessel ] ) takes place in order to satisfy the boundary condition at infinity ( eq .  ( [ eq : u_infty ] ) ) . figure [ fig : u(r ) ] sketches the behavior of the meniscus profile @xmath15 . as can be checked directly in the differential equation  ( [ eq : younglaplace ] ) , one has @xmath90 asymptotically as @xmath91 , which expresses the balance between the surface force @xmath20 on the meniscus and the gravitational force . there is , however , another intermediate regime @xmath92 but not too large , in which @xmath15 decays @xmath93 and which corresponds to the solution of eq . ( [ eq : younglaplace ] ) with @xmath94 set to zero , i.e. , when gravity is balanced by the laplace pressure induced by the meniscus curvature . these conclusions , as well as the functional dependence of @xmath15 on @xmath20 , are robust and independent of the details of the boundary condition at @xmath95 . in order to support this statement , we consider two cases which implement the distant boundary condition differently :    * _ pinned interface_. in the absence of gravity the interface is assumed to be pinned at a finite distance @xmath96 from the colloid . this corresponds to setting @xmath84 in eq . ( [ eq : younglaplace ] ) and replacing eq . ( [ eq : u_infty ] ) by the new boundary condition @xmath97 . in this case the solution of eq . ( [ eq : younglaplace ] ) is given by @xmath98 where @xmath99 in analogy to eq . ( [ eq : epspi ] ) . in the intermediate asymptotic regime ( @xmath100 ) , the meniscus profile is then given by @xmath101 assuming that @xmath20 decays sufficiently fast so that the integrals in eqs . ( [ eq : sol_pin ] ) and ( [ eq : epspi_pin ] ) converge as @xmath102 . equation ( [ eq : sol2 ] ) resembles eq . ( [ eq : sol1 ] ) with @xmath103 replaced by @xmath96 . * _ pinned curved reference interface_. in some experiments the interface between phases 1 and 2 is in fact closed , so that the colloidal particle lies at the surface of a large nonvolatile spherical droplet of phase 2 which is immersed in phase 1 @xcite ( fig . [ fig : droplet ] ) and fixed by certain means ( e.g.  by a glass plate ) . + the free energy functional  ( [ eq : f ] ) has to be modified to account for the curvature of the reference interface as well as of the constraint that the droplet volume remains unchanged under deformation . to determine the interface deformation , we minimize the functional and employ the boundary condition that the droplet is fixed at some point far from the colloid . the mathematical details and the corresponding solution for @xmath15 are presented in appendix  [ sec : droplet ] . here we quote only the intermediate asymptotic behavior ( @xmath104 , where @xmath105 is the radius of the undeformed droplet ) : @xmath106 with @xmath107 a numerical constant given by eq . ( [ eq : tildec ] ) . again eq . ( [ eq : sol3 ] ) closely resembles eqs . ( [ eq : sol1 ] ) and ( [ eq : sol2 ] ) . the physical reason for the occurrence of the singularity in the limit @xmath108 is that a `` restoring force '' far from the particle is required to yield a well  defined unperturbed interface which allows one to determine the deformation @xmath15 unambiguously . for example , if one takes the limit @xmath109 in the young  laplace eq . ( [ eq : younglaplace ] ) , it is inconsistent to impose the boundary condition @xmath110 in the corresponding solution . the special case @xmath111 , however , is not singular ; this will be discussed in sec . [ sec : discussion ] . furthermore , the comparison of eqs . ( [ eq : sol1 ] ) , ( [ eq : sol2 ] ) , and ( [ eq : sol3 ] ) demonstrates that the functional form of the intermediate asymptotic behavior is independent of how the restoring force is implemented . this corresponds to the so  called intermediate asymptotic behavior of the second kind @xcite , characterized by the following features :    1 . there is a length scale ( @xmath49 , @xmath96 , or @xmath105 ) which is much larger than the other length scales of the system under consideration and which seems  at first sight  to be irrelevant . 2 .   nevertheless , this length scale determines the dominant logarithmic ( or more generally , the power  law ) dependence . 3 .   the detailed physical origin of this length scale ( in the examples we have considered , gravity , pinning of a reference flat or curved interface ) does not matter . well known examples of this kind of asymptotic behavior are critical phenomena in phase transitions @xcite . in that case , it is the microscopic length scale given by the amplitude of the correlation length which can not be set to zero although it is much smaller than the correlation length itself . this microscopic length scale is required to formulate the power  law behavior of certain properties of the system , but its detailed physical origin is unimportant for the universal decay exponents of the power laws . in this section we consider the equilibrium state of two identical colloids floating at the interface at a fixed lateral distance @xmath112 and compute the effective interaction potential @xmath113 generated by the meniscus . the free energy can be derived along the same lines leading to eq . ( [ eq : f ] ) , but with due account for the fact that in the presence of two colloids the meniscus slope no longer exhibits rotational symmetry . nonetheless @xmath114 depends only on the distance @xmath112 between the centers of the two spheres . the reference configuration is that of two colloids floating on a planar interface with the corresponding reference free energy being independent of @xmath112 . in this case one has for the free energy relative to that of the reference configuration a contribution @xmath115 from the meniscus and a contribution of the form @xmath116 from each colloid . ( the total free energy includes also the direct interaction between the colloids ; this contribution will be considered in subsec . [ sec : elec ] . ) from eqs . ( [ eq : fcont_aniso ] ) , ( [ eq : fmen ] ) , ( [ eq : fvol ] ) , ( [ eq : finter ] ) , and ( [ eq : fcoll ] ) one obtains : @xmath117 + \sum_{\alpha=1,2 } \left\ {      \frac{\gamma}{2 r_{0 , \rm ref } } \oint_{\partial s_\alpha } \!\!\ ! d\ell \ ; [ \delta \hat{h}_\alpha - \hat{u}]^2       - \hat{f } \delta \hat{h}_\alpha \right\ } .\ ] ] here , @xmath118 is the meniscus profile in the presence of two colloids , @xmath119 are the corresponding heights , @xmath120 is the vertical force per unit area on the meniscus in the reference configuration , and @xmath121 is the force on each colloid . by symmetry , one has @xmath122 and @xmath123 . @xmath124 are the circular disks delimited by the contact lines of the colloids in the reference configuration , @xmath125 are the contact lines , with the convention that we trace them counterclockwise , and @xmath126 ( see fig . [ fig:2coll_top ] ) . the equilibrium configuration is the minimum of the free energy given by eq . ( [ eq : f2 ] ) . the minimization procedure follows closely subsec . [ sec : minf ] . first , minimizing with respect to @xmath127 at fixed meniscus height @xmath118 leads to the height of the colloids , @xmath128 where @xmath129 is the mean height of the contact line . in the next step , we minimize with respect to @xmath118 at fixed @xmath119 . variation in the interior of the domain @xmath130 provides a second  order partial differential equation , @xmath131 while variation at the boundary @xmath132 provides the following _ transversality conditions _ @xcite : @xmath133 where the last equality follows from eq . ( [ eq : h_eq2 ] ) . in this expression , @xmath134 is the derivative in the outward normal direction of @xmath125 . ( in this way , the triad ( @xmath135 ) is right - handed , where @xmath136 is the unit vector in the outward normal direction , @xmath137 is the unit vector in the counterclockwise tangent direction , and @xmath138 is the unit vector in the positive @xmath74direction . ) when applying the transversality condition , in the context of gauss theorem one must keep in mind that the boundary of the region @xmath130 consists of the contours @xmath125 traced in _ clockwise direction _ with the normals directed towards the interior of @xmath124 . ( the boundary at @xmath139 does not contribute due to eq . ( [ eq : u_infty2 ] ) . ) in the special case of a single colloid , rotational invariance reduces eq .  ( [ eq : trans ] ) to eq . ( [ eq : u0 ] ) . finally , one has the additional boundary condition @xmath140 the solution of eq . ( [ eq : younglaplace2 ] ) with the boundary conditions given by eqs . ( [ eq : trans ] ) and ( [ eq : u_infty2 ] ) is a difficult task . ( we are only aware of the  already very involved  solution of the homogeneous eq . ( [ eq : younglaplace2 ] ) , i.e. , @xmath141 , with the simplified boundary condition @xmath142 @xcite . ) for the present purpose , one can use the so - called _ superposition approximation _ @xcite , which yields the correct solution in the asymptotic limit of large separation @xmath143 between the colloids . let @xmath144 denote the equilibrium meniscus profile as if colloid @xmath145 was alone , with @xmath146 and @xmath147 denoting the corresponding forces . the superposition approximation then reads : @xmath148 notice that the fields @xmath149 and @xmath150 are defined in the domain @xmath151 , while the fields @xmath152 and @xmath153 are defined in the smaller domain @xmath130 . equations  ( [ eq : younglaplace2 ] ) and ( [ eq : u_infty2 ] ) are fulfilled by this approximate solution , but the boundary condition  ( [ eq : trans ] ) is violated : using eqs .  ( [ eq : superposition ] ) and the boundary condition in eq .  ( [ eq : u0 ] ) for the single  colloid solution , one obtains @xmath154_{{\bf r } \in \partial s_1 }     \simeq \left [ \frac{\partial { u}_2 ( { \bf r})}{\partial n_1 } -       \frac{1}{r_{0 , \rm ref } } \left\ { { u}_2 ( { \bf r})-        \frac{1}{2 \pi r_{0 , \rm ref } } \oint_{\partial s_1 } \!\!\ ! d\ell \ ; { u}_2 \right\ }     \right]_{{\bf r } \in \partial s_1 } , \ ] ] and a similar expression for the other colloid with the indices 1 and 2 interchanged . in general , this expression does not vanish as required by eq . ( [ eq : trans ] ) . if @xmath112 is large , eq . ( [ eq : trans1 ] ) can be evaluated by expanding @xmath155 into a taylor series around the center of colloid 1 , yielding to lowest order ( @xmath156 is the outward directed normal unit vector of @xmath157 ) @xmath158_{{\bf r } \in \partial s_1 }     \simeq \frac{1}{4 } r_{0 , \rm ref } [ { \bf e}_1 { \bf e}_1 : \nabla \nabla { u}_2 ( d ) +     \nabla^2 { u}_2 ( d ) ]     \ ] ] in dyadic notation . inserting the single  colloid solution given in eq . ( [ eq : sol1 ] ) , one finds that this expression decays like @xmath159 if @xmath160 , and like @xmath161 if @xmath162 . the superposition solution can be used to determine the vertical displacement according to eq . ( [ eq : h_eq2 ] ) : @xmath163 where @xmath164 is the relative vertical displacement of an isolated colloid ( eq .  ( [ eq : h_eq ] ) ) and @xmath165 is the average of the single  colloid meniscus height at the contact line of the other colloid . the meniscus  induced effective potential between the two colloids ( without their direct interaction ) is defined as @xmath166 where @xmath167 is the free energy of the two  colloid equilibrium configuration ( eq .  ( [ eq : f2 ] ) ) while @xmath168 is the single  colloid equilibrium free energy ( eq . ( [ eq : f ] ) ) . as noted before the energy of the reference configuration is independent of the separation @xmath112 and drops out from eq . ( [ eq : veff_a ] ) . we insert eqs . ( [ eq : superposition ] ) , ( [ eq : z_superposition ] ) , and ( [ eq : meanu ] ) into eq . ( [ eq : veff_a ] ) and exploit the invariance of the free energy under exchange of the colloids , i.e. , the symmetry under exchanging indices @xmath169 . after carrying out some analytic manipulations one finds the following expression for the effective potential : @xmath170 \nonumber \\    & & \mbox { } - \int_{s_1 } \!\!\ ! da \ ; \left [       \gamma |\nabla { u}_2|^2 +       \frac{\gamma}{\lambda^2 } { u}_2 ^ 2 -       2 { \pi}_2 { u}_2 \right ] \nonumber \\    & & \mbox { }   +     \frac{\gamma}{r_{0 , \rm ref } } \oint_{\partial s_1 } \!\!\ ! d\ell \ ; [ \bar{{u } } - { u}_2]^2 - 2 f \bar{{u}}\;.\end{aligned}\ ] ] the first integral accounts for the change in surface energy , gravitational potential energy , and surface  stress potential energy of the meniscus due to the overlap of the meniscus deformations caused by the two colloids . the second integral is the corresponding change due to the fact that the interface is reduced by an amount @xmath171 compared to the single  colloid case because of the presence of the second colloid . the third integral is the change in surface free energy of one colloid due to the extra meniscus deformation induced by the second colloid . the last term is the change in energy due to the vertical displacement of one colloid by this extra meniscus deformation .    for the mathematical manipulations to follow , it is suitable to rewrite eq . ( [ eq : vmen1 ] ) by applying gauss theorem to the integrals involving @xmath172 and by using the fact that the functions @xmath173 fulfill eq . ( [ eq : younglaplace ] ) individually , e.g. , @xmath174 = \nonumber \\    & & \mbox { } - \gamma \oint_{\partial s_1 } \!\!\ ! d\ell \ ; \frac{\partial ( { u}_1 { u}_2)}{\partial n_1 } +    \int_{\mathbb{r}^2\backslash(s_1 \bigcup s_2 ) } \!\!\!\!\!\ ! da \ ; \left [ { \pi}_1 { u}_2 -       \frac{\gamma}{\lambda^2 }   { u}_1 { u}_2 \right ] . \nonumber \\ \end{aligned}\ ] ] thus one obtains from eq . ( [ eq : vmen1 ] ) @xmath175 ^ 2     - 2 f \bar{{u } } .\end{aligned}\ ] ] in this form all the integrals , except the first one , are performed over bounded domains . this allows one to carry out a taylor expansion which yields a uniformly valid asymptotic expansions of the terms as @xmath176 . in the following calculations , the integrals are evaluated in polar coordinates with the origin at the center of @xmath124 and @xmath177 is the position vector with respect to this origin . in particular , @xmath178 is the center of colloid 1 with respect to the center of colloid 2 , and @xmath179 is the center of colloid 2 with respect to the center of colloid 1 ( fig . [ fig:2coll_top ] ) so that , e.g. , @xmath180 .    1 .   in order to compute the leading asymptotic behavior of the first term in eq . ( [ eq : vmen2 ] ) as @xmath176 we distinguish two cases : 1 . @xmath85 : in this case , @xmath181 ( eq .  ( [ eq : sol1 ] ) ) while @xmath182 ( @xmath27 ) as @xmath183 . asymptotically the main contribution stems from the region near @xmath171 and @xmath184 employing the definition  ( [ eq : epspi ] ) . it will turn out that there is no need to compute also the next  to  leading term in order to obtain @xmath185 as the leading term will be part of the leading contribution to @xmath114 and it will not be cancelled by other contributions . @xmath186 : in this case @xmath187 if @xmath188 . the regions contributing mainly to the integral are the neighborhoods of @xmath171 and @xmath189 . we employ the taylor expansion of the integrand ( which is valid up to a maximum order depending on how fast @xmath94 decays ) to evaluate the leading and the next  to  leading contributions if @xmath190 ( using dyadic notation ) : @xmath191_{{\bf r}_2 = d { \bf e}_x } \nonumber \\      & & \mbox { } + \int_{\mathbb{r}^2\backslash s_2 } \!\!\!\!\!\ ! da \ ; { u}_2 \left[{\pi}_1 +        \dots \right]_{{\bf r}_1 = - d { \bf e}_x } \nonumber \\      & \simeq & 2 \pi \gamma r_{0 , \rm ref}\ ; \varepsilon_\pi \;{u}(d ) +      \frac{1}{2 } \pi \nabla^2 { u } ( d ) \int_{r_{0,\rm ref}}^\infty dr\ ; r^3\,\pi(r ) \nonumber \\      & & \mbox { } + { \pi}(d ) \int_{\mathbb{r}^2\backslash s_2 } \!\!\!\!\!\ ! da \ ; { u}_2 + \dots \ ; . \end{aligned}\ ] ] we can simplify the result further by evaluating the last integral by repeated partial integration with the explicit solution for @xmath192 given in eq . ( [ eq : sol1 ] ) : @xmath193 .\ ] ] one finally obtains @xmath194 { \pi } ( d ) + \dots \nonumber      \end{aligned}\ ] ] where we have employed eq . ( [ eq : younglaplace2 ] ) ( with @xmath195 for simplicity ) . convergence of the integral of @xmath192 imposes the more stringent constraint @xmath190 on the decay of @xmath196 ( see eq . ( [ eq : int_uda ] ) ) . + the validity of formulae  ( [ eq : pi1u2_a ] ) and ( [ eq : pi1u2_b ] ) has been checked explicitly by comparing their predictions with the numerical evaluation of the integral for a surface stress of the form @xmath197 with @xmath190 . 2 .   the second term in eq .  ( [ eq : vmen2 ] ) can be estimated in the limit @xmath176 by expanding the integrand into a taylor series , which is uniformly valid in the integration domain : @xmath198_{{\bf r}_2 = d { \bf e}_x }               \simeq \pi r_{0 , \rm ref}^2 \ , { \pi } ( d ) \ , { u } ( d ) +               \dots \ ; . \end{aligned}\ ] ] 3 .   the third term in eq . ( [ eq : vmen2 ] ) reads : @xmath199_{r_1=r_{0 , \rm ref } } \\      & \simeq & \int_0^{2 \pi } \!\!\ ! d\varphi \ ; \left\ { ( { \bf r}_1 \cdot \nabla )       \left [ { u}_1 \left (      \phantom { \frac{1}{2 } } \!\!\ !     { u}_2 + { \bf r}_1 \cdot ( \nabla { u})_2 +      \right.\right.\right.\nonumber \\      & & \qquad \qquad \left . \left . \left . \frac{1}{2 } { \bf r}_1 { \bf r}_1 : ( \nabla \nabla { u})_2        + \dots \right)_{{\bf r}_2=d { \bf e}_x } \right ] \right\}_{r_1=r_{0 , \rm ref } } \nonumber \\    & \simeq & 2 \pi r_{0 , \rm ref }   { u}(d ) \left . \frac{d { u}}{d r } \right|_{r = r_{0 , \rm ref } } +     \frac{\pi}{2 } r_{0 , \rm ref } \nabla^2 { u}(d ) \left . \frac{d ( r^2 { u})}{d r } \right|_{r = r_{0 , \rm ref } } +     \dots \nonumber \\    & \simeq & 2 \pi r_{0 , \rm ref } \varepsilon_f { u}(d ) -     \frac{\pi}{2\gamma } r_{0 , \rm ref}^3 \left(\varepsilon_f + \frac{2{u}_0}{r_{0 , \rm ref } } \right ) { \pi}(d ) + \dots      \end{aligned}\ ] ] using the boundary condition for @xmath192 in eq .  ( [ eq : u0 ] ) . 4 .   the fourth term in eq .  ( [ eq : vmen2 ] ) is evaluated analogously : @xmath200_{r_1=r_{0 , \rm ref } } \\      & \simeq & \int_0^{2 \pi } \!\!\ ! d\varphi \ ;       \left\ { { \bf r}_1 \cdot \left [ ( \nabla { u}^2)_2 + { \bf r}_1 \cdot ( \nabla \nabla { u}^2)_2        + \dots \right]_{{\bf r}_2=d { \bf e}_x } \right\}_{r_1=r_{0 , \rm ref } } \nonumber \\      & \simeq & \pi r_{0 , \rm ref}^2 \nabla^2   { u}^2 ( d ) + \dots \ ; .    \end{aligned}\ ] ] 5 .   the fifth term in eq . ( [ eq : vmen2 ] ) is given by @xmath201 ^ 2 & = &      \int_0^{2 \pi } \!\!\ ! d\varphi \ ; { u}_2 ^ 2 - 2 \pi \bar{{u}}^2 \nonumber \\      & \simeq & \frac{1}{2 } \pi r_{0 , \rm ref}^2 \left [ \nabla^2 { u}^2 ( d ) -        2 { u } \nabla^2 { u } ( d ) \right ] + \dots \nonumber \\     & \simeq & \frac{1}{2 } \pi r_{0 , \rm ref}^2 \left [ \nabla^2 { u}^2 ( d ) +        \frac{2}{\gamma } { \pi } ( d ) \ , { u } ( d ) \right ] + \dots   \;.        \end{aligned}\ ] ] 6 . in order to evaluate the sixth term in eq . ( [ eq : vmen2 ] ) , we expand @xmath202 in the definition given by eq . ( [ eq : meanu ] ) : @xmath203    the asymptotic behavior of the effective potential is finally obtained from eq . ( [ eq : vmen2 ] ) by collecting all terms . there are two qualitatively different cases :    * : the limit @xmath204 is singular and the single  colloid meniscus profile exhibits a logarithmic dependence . one finds that the leading contribution is provided by the terms proportional to @xmath205 : @xmath206 which represents a long - ranged attractive effective potential , irrespective of the sign of the forces @xmath60 and @xmath94 acting on the system . * : the single  colloid meniscus profile decays like @xmath207 if @xmath161 . the leading contribution is proportional to @xmath208 , because the terms proportional to @xmath205 cancel each other : @xmath209 this correspond to a shorter  ranged effective interaction , which in principle can be either attractive or repulsive depending on the precise form of the function @xmath210 . in the particular case that @xmath210 decays monotonically to zero , e.g. , @xmath211 , it is easy to check that @xmath114 amounts to a repulsive force , because the sign of @xmath192 is opposite to the sign of @xmath196 . we have seen that the error committed by the superposition approximation in satisfying one of the boundary conditions , eq . ( [ eq : correction_sup ] ) , decays like @xmath208 , too . this suggests that the corrections to the superposition approximation could modify the precise value of the constant factor acting as an amplitude in eq . ( [ eq : shortrangev ] ) , and _ a priori _ it can not be excluded that there are cancellations leading to a vanishing amplitude , and therefore to an even faster decay for large @xmath112 . thus the superposition approximation might not be reliable enough for calculating @xmath114 if @xmath186 . if one traces back the origin of the dominant contributions to @xmath212 , one finds that in both cases only the first , third , and sixth term in eq . ( [ eq : vmen2 ] ) are relevant . they correspond physically to the effect of the overlap of the two single  colloid meniscus profiles and the change of the colloid height ( first integral and the term @xmath213 in eq . ( [ eq : vmen1 ] ) , respectively ) . equations ( [ eq : logv ] ) and ( [ eq : shortrangev ] ) describe the asymptotic behavior of the meniscus  induced effective intercolloidal potential and thus represent a central result of our analysis . they provide the explicit functional dependence on an arbitrary stress field @xmath20 which decays sufficiently fast . the assumptions entering their derivation are ( i ) that the deviations of the meniscus profile from the reference configuration are small , allowing one to confine the analysis to a free energy expression which is quadratic in the deviations , and ( ii ) the superposition approximation , which expresses the two  colloid equilibrium state in terms of the single  colloid state . the analysis shows that the limit @xmath214 is non  singular only in the case @xmath162 . equation ( [ eq : logv ] ) can be used to determine the flotation force . there are no stresses acting at the meniscus , @xmath215 , while the force @xmath60 on the colloids is due to their weight corrected by the buoyancy force . accordingly , eq .  ( [ eq : logv ] ) reduces to the flotation potential @xmath216 where @xmath217 is known as the capillary charge @xcite , by analogy with two  dimensional electrostatics . the asymptotic form for @xmath218 originates from a potential proportional to the modified bessel function @xmath219 ( see eq . ( [ eq : solbessel ] ) ) @xcite . the order of magnitude of the capillary charge is @xmath220 . for a typical value at room temperature and for colloids with a mass density of the order of , the prefactor of the logarithm can be estimated as @xmath221 therefore , compared with the thermal energy , the flotation force is negligible for submicrometer  sized colloids . another application is the case of electrically charged colloids . if one of the liquid phases is water , the charge of the colloid is screened ( the debye length of pure water is @xmath222 @xmath0 m and smaller in the presence of an electrolyte ) , and the effective electric field is that of a dipole oriented perpendicularly to the fluid interface @xcite . the electrostatic field decays as @xmath223 and the stress on the meniscus as @xmath224 . both the electrostatic stress and the osmotic ionic pressure decay in the same manner @xcite . thus the total intercolloidal potential at intermediate distances is @xmath225 if gravity is neglected both as a force on the colloid and as a restoring force for the interface @xmath226 , then one can indeed show that the ensuing condition of mechanical isolation ( no net force on the system ) leads to @xmath186 , i.e. , precisely the situation for which the limit @xmath78 is non  singular . to see this , we consider the total stress tensor @xmath227 which consists of the maxwell stress tensor ( due to the electrostatic field ) and a diagonal osmotic pressure tensor ( due to the electrolytes ) @xcite . at interfaces @xmath227 can be discontinuous . the total volume @xmath228 of the system is divided into volumes @xmath229 , @xmath230 , and @xmath231 ( see fig .  [ fig : forcebalance ] for the explanation of the notation in the following equation ) . the total force reads ( the superscript @xmath232 denotes evaluation on the positive(negative ) side of the oriented surface , i.e. , the side the arrows in fig .  [ fig : forcebalance ] point to ( do not point to ) ) @xmath233 \nonumber \\    & = & \int_{s_{\rm men , ref } } \!\!\!\!\!\ ! da \ ; \pi   \ , { \bf e}_z   + f { \bf e}_z \ ; . \ ] ] in the first line , we have applied gauss theorem with due account of the possible discontinuities of the tensor @xmath227 accross the interfaces . in the second line , @xmath234 in the fluid phases @xmath229 and @xmath230 , because the counterion distribution of the present reference configuration is the equilibrium distribution and thus locally force free . this distribution is considered to be fixed . the second term in the second line is the total force on the meniscus ( which can have only a normal component ) , while the terms in square brackets sum up to the force @xmath60 acting on the colloid . thus the vertical component of the total force is @xmath235 if it vanishes , as it is the case for an isolated system , then @xmath236 . according to eq . ( [ eq : logv ] ) this implies that the long  ranged logarithmic contribution to @xmath114 is absent and thus the limit @xmath78 is regular . physically this means that there is no need for a restoring force acting on the fluid  fluid interface when the deformation is created by localized internal stresses . instead , according to eq . ( [ eq : shortrangev ] ) , one obtains a potential @xmath237 for the present case of a dipolar electric field ( see above ) . this shorter  ranged potential can not counterbalance the direct electrostatic dipolar repulsion @xmath238 . such a counterbalance would be needed for a straightforward explanation of the aforementioned experimentally observed attractions . the line of argument to explain the absence of the logarithm contribution to @xmath114 was already put forward in refs . @xcite and @xcite , where exclusively the case @xmath239 has been addressed . our detailed analysis complements and generalizes these contributions . in ref . @xcite , @xmath240 is estimated by taking into account only the degree of freedom `` meniscus profile '' , @xmath241 , and considering only the change in meniscus area due to the superposition of the dimples . this corresponds to retaining only the term @xmath242 in expression  ( [ eq : vmen1 ] ) . although the type of @xmath112dependence obtained that way is correct , the sign of the force turns out to be wrong ( attraction instead of repulsion ) . the reason is that the contributions @xmath243 and @xmath244 in eq . ( [ eq : vmen1 ] ) are equally important as the retained term . this is taken into account in the more detailed analysis of ref . @xcite , where the limiting case @xmath245 ( point  particle ) is considered from the outset , so that `` height of the colloid '' is not an independent degree of freedom , i.e. , @xmath246 . correspondingly , @xmath60 is set to zero ( @xmath247 ) , and its effect is modelled by a dirac  delta contribution to the stress @xmath20 such that @xmath248 . the potential @xmath249 calculated this way corresponds to keeping the three terms mentioned above which are relevant ( the first and third term in the first integral of eq . ( [ eq : vmen1 ] ) and the last term of eq . ( [ eq : vmen1 ] ) ) and to setting @xmath250 , since this limit is regular when @xmath251 . our analysis has shown that the terms which are dropped in the limit @xmath252 ( second and third integral in eq . ( [ eq : vmen1 ] ) , second , fourth , and fifth terms in eq . ( [ eq : vmen2 ] ) ) yield indeed a subdominant contribution to eq . ( [ eq : shortrangev ] ) . however , the integral appearing as a prefactor in eq . ( [ eq : shortrangev ] ) is divergent for @xmath252 , so that in ref . @xcite a short  distance cutoff @xmath253 has been introduced which is expected to be of the order of @xmath8 . the precise value of @xmath253 depends on the details of the implementation of this cutoff ; for the example used in ref . @xcite , the application of eq . ( [ eq : shortrangev ] ) yields @xmath254 .    in the presence of gravity mechanical isolation is violated and as we have shown @xmath255 . ( the force  balance argument ( eq .  ( [ eq : forcebalance ] ) ) can be easily generalized to include the effect of the gravitational volume force . ) in the case of a curved reference interface ( corresponding to the experiment reported in ref . @xcite using colloids trapped at the interface of water droplets in oil ) , eq . ( [ eq : forcebalance ] ) is replaced by @xmath256 where @xmath257 is the radial unit vector of the unperturbed spherical droplet . force balance in the vertical direction yields a factor @xmath258 ( see fig . [ fig : droplet ] ) in the integral , which can be expanded for small angles in the limit @xmath259 . this leads to a curvature  induced correction @xmath260 even for mechanical isolation . thus we see that independently of the details of the implementation of the distant boundary conditions the logarithmic term in the potential @xmath114 ( eq .  ( [ eq : logv ] ) ) has a strength which is proportional to ( colloid size / system length)@xmath261 , which is nevertheless too weak to explain the reported attractive total interaction . mechanical isolation is violated in the presence of an external homogeneous electric field @xmath262 . for the experimental setup used in ref . @xcite it can not be ruled out that such an external field may have distorted the measurements @xcite . since this is the only experiment which provides quantitative information about the secondary minimum in @xmath263 , we discuss the case of an external field in more detail . in this experiment position  recording measurements were performed on a hexagonal configuration of seven trapped colloids on a water droplet immersed in oil @xmath264 m ) . the latter was confined between two glass plates and the droplet stuck on the upper glass plate with a contact angle close to @xmath265 . residual charges might have resided on the upper plate @xcite . the measurements yielded the position of the secondary minimum , @xmath266 , and the curvature at the minimum , @xmath267 . with regard to the latter value we remark that systematic corrections can be estimated to lower @xmath268 by a factor 2 to 3 . these systematic corrections include ( i ) multiparticle effects and ( ii ) center  of - mass movement of the droplet . we have estimated the effect of ( i ) by carrying out langevin simulations of the seven particle system using the intercolloidal potential eq . ( [ eq : vtot1 ] ) below . as for the effect of ( ii ) , any shape deformation induced by the moving colloids changes the center  of  mass position of the droplet since the droplet is fixed to the upper glass plate . the corresponding change in the gravitational energy of the droplet translates into a weak confining potential for the colloids which limits the stochastic movement of the center  of  mass of the seven colloids . this effect might be part of an explanation for the absence of center  of  mass movement observed in ref . @xcite . according to eqs . ( [ eq : logv ] ) and ( [ eq : vtot ] ) , the total intercolloidal potential in an external field is @xmath269 using the aforementioned experimental data for @xmath270 and @xmath271 , one obtains from eq . ( [ eq : vtot1 ] ) @xmath272 , @xmath273 and @xmath274 which is surprisingly deep , even if reduced by a factor 2 to 3 due to the systematic corrections mentioned above . furthermore , from eq . ( [ eq : logv ] ) we deduce @xmath275 and with @xmath276 n / m we find @xmath277 . the long  ranged meniscus deformation ( see eq . ( [ eq : sol1 ] ) ) is on the scale of nanometers . the short  ranged meniscus deformation near the colloid can only be evaluated with a specific microscopic model for @xmath20 . for a rough estimate of @xmath30 , we consider the colloid charge to be concentrated on the surface . the asymptotic behavior of the stress tensor in this case is given by @xmath278 @xcite . if we assume this form of the stress tensor to hold for all @xmath89 , we find @xmath279 for the values of @xmath10 and @xmath253 given above . for not too small contact angles , this _ a posteriori _ justifies the perturbative approach which we have adopted .    if the system has a net charge @xmath280 , then @xmath281 ) and ( [ eq : logv ] ) ) . using the values for @xmath282 and @xmath10 as given above we find @xmath283 for the value @xmath284 quoted in ref . @xcite , this yields a relatively small electric field @xmath285 , indicating how sensitive the system can be to spurious external fields . alternatively , an electric field @xmath286 is sufficient for the meniscus  induced logarithmic potential to have a depth of the order @xmath287 . thus the external field offers the possibility to tune easily the capillary long  ranged attraction and to manipulate the structures formed by the colloids at the interface .      in ref . @xcite , experimental results for the meniscus deformation around glass particles of radii 200  300 @xmath0 m trapped at water  air or water  oil interfaces are reported . the data for the meniscus slope @xmath288 at the contact circle imply ( eq .  ( [ eq : u0 ] ) ) @xmath289 ( water  oil interface ) , about fifteen times larger than the corresponding @xmath290 caused by gravity alone . furthermore , the reported meniscus shape for one sample contains a logarithmic part which is consistent with @xmath291 ( see eq . ( [ eq : sol1 ] ) ) . as in the experiment analyzed in the previous subsection , the experimental observations could be understood within the framework of the theoretical model we have developed in terms of an external electric field violating mechanical isolation . yet independently of us , the authors of ref . @xcite have developed a theoretical model based on the same physical hypotheses as our approach which , they claim , explains the observed long  ranged meniscus deformation . here we would like to note two important errors which flaw their analysis :    1 . ( 3.15 ) in ref . @xcite , in terms of which they interpret their observations , can be obtained by inserting the large  distance ( @xmath292 ) asymptotic behavior of @xmath20 in our eq . ( [ eq : sol1 ] ) ( equivalent to their eq .  ( 3.14 ) up to an additive constant ) both into the integral term _ and into the definition of @xmath30 _ ( eq . ( [ eq : epspi ] ) ) . since the dominant contribution to the integral defining @xmath30 stems from points @xmath293 , this procedure is clearly inadmissible . as a consequence , they obtain a wrong , nonvanishing prefactor of the logarithm , in spite of their explicit consideration of mechanical isolation ( see their eq . ( 6.6 ) ) . 2 .   in order to calculate the intercolloidal effective potential within the superposition approximation , the formula relating @xmath294 to the meniscus deformation @xmath295 _ in the presence of gravity alone _ ( i.e. , eq .  ( [ eq : flotv ] ) ) is applied even though @xmath296 . thus , an additional term contributing to the leading logarithm is not included ( compare their eq . ( 3.16 ) and our eq . ( [ eq : logv ] ) ) .    in ref . @xcite also a numerical analysis is carried out . a detailed study of the relation between the results of this numerical analysis and our theoretical predictions will be published elsewhere  @xcite . it may be possible that the presence of an external field is consistent with the data from refs . @xcite , as well as with the presence of the secondary potential minimum observed in the experiments using planar troughs , in particular in the cases @xmath297 , which fall into the intermediate asymptotic regime considered here . but this still remains as an open problem . given that already nanometer distortions of the meniscus produce noticeable attractions , the surface topography of colloids might be relevant . in ref . @xcite , colloidal surface roughness is proposed as an explanation of the attraction . the meniscus contact line is assumed to be pinned at defects on the colloid surface caused by surface roughness . this imposes a different boundary condition for the meniscus at contact , which is then deformed even in the absence of electrostatic forces . the corresponding analysis in ref . @xcite is concerned only with the term @xmath242 in the expression  ( [ eq : vmen1 ] ) , leading to the conclusion that @xmath113 decays as @xmath298 and corresponds to an attractive potential with a strength of @xmath299 for meniscus deformations of the order of @xmath300 . this conclusion , however , can not be simply carried over to the case of charged colloids . as we have shown the contributions of the terms @xmath301 and @xmath244 in eq . ( [ eq : vmen1 ] ) are relevant and can change the qualitative behavior of @xmath113 even in the limit of point  colloids . it would be worthwhile to generalize the analysis of ref . @xcite along the line of arguments presented here in order to assess the importance of surface roughness . this should be complemented by more precise experimental information about the actual colloidal topography . recently , an explanation based on a contaminated interface has been advocated  @xcite . the air  water interface would be actually a two  dimensional emulsion consisting of hydrophylic ( water ) patches and hydrophobic ( silicon oil ) patches . the colloidal particles ( hydrophobic in character according to ref . @xcite ) would cluster in the hydrophobic patches . thus , confinement of the colloids by finite  size hydrophobic patches would give the impression of an effective intercolloidal attraction . at present , this explanation is only of qualitative nature . thermal fluctuations of the interface position around its mean value @xmath302 also induce an effective interaction between the colloids which confine these fluctuations ( casimir  like force ) . using a gaussian model of the fluctuating interface in analogy to the procedure employed in ref . @xcite for calculating fluctuation  induced forces between rods in a membrane , one finds for uncharged colloids a fluctuation  induced potential @xmath303 , which is too small and falls off more rapidly than the intercolloidal dipole repulsion . here , the generalization to the charged case might give hints for the effective attractions between _ very _ small particles trapped at interfaces . concerning particle sizes well below the debye length , one should modify our analysis to account for the overlap of the screening ionic clouds ( this would affect , e.g. , the superposition approximation for the stress field @xmath94 , eq . ( [ eq : superposition ] ) ) .    finally , in ref . @xcite the attraction of particles trapped at a nematic  air interface is reported and an explanation in terms of a logarithmic meniscus deformation is proposed which parallels the explanation given in ref . @xcite : in this case , the deformation would be caused by the elastic distortion induced by the particles on the nematic phase  @xcite . our detailed theoretical treatment shows that no long  range ( logarithmic ) meniscus distortion can arise on an interface if the system is mechanically isolated and the excess free energy of the perturbed interface is correctly described by an expression like eq . ( [ eq : f ] ) . thus it appears that the simple explanation of the observed colloidal attractions given in ref . @xcite is not correct . however , it is not clear whether the free energy of a distorted nematic  air interface is equivalent to that of a simple fluid interface due to the long  ranged interactions in the nematic bulk caused by defects . thus , a generalization of our theory to interfaces involving nematic phases would be desirable in order to assess the possibility of long  ranged colloidal attractions in more detail . we have analyzed the effective force induced by capillary deformation between two smooth spherical colloids floating at a fluid interface . the relevant degrees of freedom are the meniscus deformation and the height of the colloids ( fig . [ fig : pert ] ) , whose equilibrium values are given by minimization of a free energy functional . this functional was derived assuming small deviations from a reference configuration ( fig . [ fig : ref ] ) . it incorporates the surface tensions of the three interfaces involved ( `` colloid  fluid phase 1 '' , `` colloid  fluid phase 2 '' , and `` fluid phase 1  fluid phase 2 '' ) , the potential energy of the colloids under the action of a force @xmath60 , the potential energy of the fluid interface in an arbitrary surface stress field @xmath304 , and the potential energy due to a restoring force acting on the interface ( eqs . ( [ eq : fdef])([eq : f ] ) ) . the effective intercolloidal potential ( eq .  ( [ eq : veff_a ] ) ) was calculated in the limit of large separations by using the superposition approximation ( eq .  ( [ eq : superposition ] ) ) . we have shown in this limit that the contribution to the effective potential by the interfaces `` colloid  fluid phases '' is subdominant . if the total force acting on the system , consisting of the two colloids and the meniscus , does not vanish , the presence of the restoring force acting on the fluid interface is essential  although its precise form does not matter ( subsec . [ sec : lambda ] ) . in this case , the effective interaction is long  ranged and attractive ( eq . ( [ eq : logv ] ) ) . if the total force vanishes , the restoring force is irrelevant altogether , the effective interaction is shorter  ranged ( eq .  ( [ eq : shortrangev ] ) ) , and it can not be computed reliably within the superposition approximation .    as an application , we have considered the case of like  charged , micrometer  sized particles when the capillary deformation is due to the ensuing electrostatic field . we have discussed how one can tune the long  ranged attractive interaction by an external electric field , but we conclude that the experimentally observed attraction in an isolated system can not be explained within the present model . possible directions for future research such as colloidal surface roughness and fluctuations of the interface have been discussed . we kindly acknowledge helpful discussions on the subject of the manuscript with m.  nikolaides , a.  bausch , k.  danov , and p.  kralchevsky . in this appendix we determine the contribution to the free energy due to the exposure of the colloid to the two fluid phases . here we consider the general case of no rotational invariance , so that the meniscus height @xmath305 depends on the distance @xmath89 from the axis through the colloid center and on the angle of revolution @xmath306 around this axis ( fig . [ fig : contactline ] ) . accordingly , the auxiliary variables @xmath307 and @xmath308 also depend on the revolution angle and one has @xmath309 and @xmath310 where @xmath311 is the meniscus height at contact . the surface areas of the colloid in contact with phase 1 and phase 2 are @xmath312\ ] ] and @xmath313 respectively . the expression for the free energy in eq . ( [ eq : fcont ] ) is based on these formulae for the special case @xmath314 ( fig . [ fig : ref ] ) and on young s law ( eq .  ( [ eq : young ] ) ) : @xmath315 = \nonumber \\ & = & \frac{1}{2 } \gamma r^2 \int_0^{2 \pi } \!\!\!\ ! d \varphi \ ,   [ \cos \xi(\varphi ) - \cos \theta]^2 +    \frac{1}{2 } \gamma r^2 \int_0^{2 \pi } \!\!\!\ ! d \varphi \ ,   [ \cos^2 \theta - \cos^2 \xi(\varphi ) ] \nonumber \\ & = & \frac{1}{2 } \gamma \int_0^{2 \pi } \!\!\!\ ! d \varphi \ ,            [ u(r_0(\varphi ) , \varphi ) - \delta h]^2 +    \frac{1}{2 } \gamma \int_0^{2 \pi } \!\!\!\ ! d \varphi \ ,   [ r_0 ^ 2 ( \varphi ) - r_{0,\rm ref}^2 ] .\end{aligned}\ ] ] the second term , which arises upon completing the square , represents the change of the meniscus area which is cut out by the colloid . since @xmath45 and @xmath38 are already of first order in @xmath30 or @xmath31 , one can replace @xmath307 by @xmath25 in the first term and obtains @xmath316 ^ 2 + \frac{1}{2 } \gamma \int_0^{2 \pi } \!\!\!\ ! d \varphi \ ,   [ r_0 ^ 2 ( \varphi ) - r_{0,\rm ref}^2 ] , \ ] ] plus corrections of at least third order in @xmath30 or @xmath31 . in the special case of rotational invariance , this expression reduces to eq . ( [ eq : fcont_iso ] ) . for the purpose of sec . [ sec:2coll ] one can rewrite eq . ( [ eq : almostfcont_aniso ] ) as the line integral @xmath317 ^ 2 + [ r_0 ^ 2 - r_{0,\rm ref}^2 ] \right\ } , \ ] ] where @xmath318 is a circular disk of radius @xmath319 , here centered at the origin , and @xmath320 denotes its circumference . in this appendix we calculate the deformation of the interface of a spherical droplet due to a floating colloidal particle . the reference configuration corresponds to a colloid floating at the surface of the spherical droplet such that eq . ( [ eq : young ] ) is fulfilled . here we will consider only axially symmetric configurations ( see fig . [ fig : droplet ] ) . in terms of the polar angle @xmath321 $ ] , the distance @xmath89 measured along the unperturbed spherical droplet surface is @xmath322 . ( note that in this appendix we consider @xmath323 and @xmath324 as functions of @xmath325 instead of @xmath89 . ) we define @xmath326 and we introduce the angle @xmath327 with @xmath328 at which another boundary condition ( to be specified below ) holds reflecting the fact that the droplet is fixed in space . with this additional boundary condition we model closely the actual experimental setup of ref . @xcite , where the droplet was attached to a glass plate with a contact angle less than @xmath265 @xcite .      * in spherical coordinates , @xmath329 reads @xmath330 where @xmath331 is the differential of the solid angle , and @xmath332 is the gradient operator on the unit sphere . as before , this functional is now expanded for small deviations from the reference configuration , and we obtain @xmath333 ,    \end{aligned}\ ] ] where we have neglected terms evaluated at @xmath334 which vanish in the limit @xmath259 . * the term @xmath61 in eq . ( [ eq : finter ] ) has to be amended by the extra work upon deformation which is caused by the pressure difference in the reference configuration between the interior and the exterior of the droplet , @xmath335 : @xmath336        \nonumber \\ \label{eq : finter1 }     & \simeq & - 2 \pi \gamma \int_{\psi_{0}}^{\psi_1 } \!\!\ ! d\psi \ ; \sin\psi \ ,       \left [ \frac{r_{\rm drop}^2}{\gamma } \pi \ , u + ( 2 r_{\rm drop } u + 2 u^2)\right ] . \end{aligned}\ ] ] since @xmath337 is of order @xmath338 , we have to keep also terms @xmath339 in this expression for consistency . for simplicity we have neglected the contribution due to the change of @xmath41 during the deformation . such a term would affect the behavior of the meniscus only near @xmath334 and vanishes in the limit @xmath259 . * the change in volume of the perturbed droplet reads @xmath340 where @xmath341 is a contribution due to the spherical shape of the colloid ; it depends on @xmath38 and @xmath342 but can be safely neglected . ( this is the same kind of term encountered in calculating @xmath53 , see eq . ( [ eq : fvol ] ) . ) expanding @xmath343 for small deformations @xmath45 up to linear order yields @xmath344 physically , the deformation of the interface occurs under the constraint that the droplet volume remains unchanged . this constraint ( @xmath345 ) will be implemented in the free energy functional by means of a lagrange multiplier @xmath0 which itself turns out to be linear in the parameters @xmath30 and @xmath31 . this justifies keeping only the linear term in eq . ( [ eq : deltav ] ) . the contribution @xmath346 ( see eq .  ( [ eq : fcont_iso ] ) ) is concerned only with quantities at @xmath334 , the colloid surface , and thus in the limit of large droplet radii this contribution are the same as in the planar case . the free energy related to the work done upon moving the colloid , @xmath347 ( see eq . ( [ eq : fcoll ] ) ) , remains also unchanged . for simplicity we also neglect gravity and set @xmath84 . in conclusion , the free energy functional to be minimized is the sum @xmath348 plus the constraint term @xmath349 ( where @xmath0 is a dimensionless lagrange multiplier ) : @xmath350 +    \pi \gamma [ u_0 - \delta h ] ^2 -    f \delta h .\ ] ] the quantity @xmath351 can be interpreted as a homogeneous pressure field enforcing the constant  volume constraint . we note that the terms linear in @xmath45 present in eqs . ( [ eq : fmen1 ] ) and ( [ eq : finter1 ] ) have cancelled in the total free energy . minimization with respect to @xmath38 yields the same result as in eq . ( [ eq : h_eq ] ) . subsequent minimization with respect to @xmath352 leads to @xmath353 and ( compare with eq . ( [ eq : u0 ] ) ) @xmath354 the solution must satisfy the constant  volume constraint @xmath355 finally , the second boundary condition mentioned at the beginning expresses the physical requirement that the droplet is fixed in space . ( otherwise , the application of a localized force at the interface would shift the droplet as a whole without deforming it . ) as an example , we assume that the droplet interface is pinned at @xmath327 ( another example , force balance by suitable localized counterstresses , is studied in ref . @xcite ) : @xmath356    the general solution of the inhomogeneous legendre equation ( [ eq : spheryl ] ) is @xmath357 \pi(\sigma ) , \ ] ] where @xmath77 and @xmath358 are integration constants and @xmath359 are solutions of the homogeneous legendre equation . for notational simplicity , evaluation of a function at @xmath334(@xmath327 ) will be denoted by the subscript @xmath360 .    from eq . ( [ eq : u_1 ] ) one obtains @xmath361 from the other boundary condition  ( [ eq : u_0 ] ) it follows that @xmath362 where we have defined the auxiliary function @xmath363 .\ ] ] finally , in order to impose the integral constraint  ( [ eq : constraint ] ) we employ the following identities : @xmath364 \pi(\sigma ) =   \\ \int_{\psi_0}^{\psi_1 } \!\!\!\ ! d\sigma \ , \sin \sigma \ ; \pi(\sigma ) \ ; \int_{\psi_0}^{\sigma } \!\!\!\ ! d\psi \ ; \sin \psi \ ,     [ p(\sigma ) q(\psi ) - p(\psi ) q(\sigma ) ] = \nonumber \\ \mbox { } - \int_{\psi_0}^{\psi_1 } \!\!\!\ ! d\sigma \ , \sin \sigma \ ; \pi(\sigma ) \ ; \frac{1}{2 } \left[1-\frac{w(\sigma)}{w_0}\right ] \;. \nonumber \end{aligned}\ ] ] the first equality follows from interchanging the order of integration . the last equality follows most easily by noticing ( i ) that the @xmath325integral in the second line is , as a function of @xmath365 , a solution of the differential equation  ( [ eq : spheryl ] ) with @xmath45 replaced by this integral , @xmath366 replaced by @xmath367 , @xmath368 , and the analogues of the boundary conditions @xmath369 , and ( ii ) that the function @xmath370 is the solution to the same differential equation with the same boundary conditions . with these identities eq . ( [ eq : constraint ] ) turns into @xmath371 \;,\ ] ] where we have introduced the numerical constants @xmath372 \nonumber \\    { \cal q } & : = & \int_{\psi_0}^{\psi_1 } \!\!\!\ ! d\psi \ , \sin \psi \ ; [ q(\psi ) - q_1 ] .\end{aligned}\ ] ] we note that the constants @xmath77 and @xmath358 are determined by the two linear eqs . ( [ eq : ab1 ] ) and ( [ eq : ab2 ] ) . these results permit us to study the intermediate asymptotic regime @xmath373 . in the expressions derived above , we have to take the limit @xmath259 while keeping fixed @xmath327 , @xmath374 and @xmath375 ( so that @xmath376 ) . in this limit , the constants @xmath377 and @xmath378 are finite , while @xmath379 and @xmath380 . the integral over @xmath94 in eq . ( [ eq : ab1 ] ) simplifies as follows ( provided that the stress field @xmath381 decays sufficiently fast ) : @xmath382 where we have used @xmath383 . the integrals over @xmath94 in eqs . ( [ eq : ab2 ] ) and ( [ eq : sol_drop ] ) can be simplified in the same manner . from eq . ( [ eq : ab1 ] ) we obtain @xmath384 and eq .  ( [ eq : ab2 ] ) yields @xmath385 one finds indeed that the lagrange multiplier @xmath0 determined by eq . ( [ eq : mu ] ) is linear in @xmath386 . finally , the general solution  ( [ eq : sol_drop ] ) simplifies to @xmath387 which renders eq . ( [ eq : sol3 ] ) with the constant @xmath388 } .\ ] ] this constant is finite and non  vanishing provided @xmath327 is not close to @xmath265 . the case @xmath389 is pathological because @xmath390 for a smooth profile , so that the boundary condition  ( [ eq : u_1 ] ) would overdetermine the problem . r.  aveyard , b.  p.  binks , j.  h.  clint , p.  d.  i.  fletcher , t.  s.  horozov , b.  neumann , v.  n.  paunov , j.  annesley , s.  w.  botchway , d.  nees , a.  w.  parker , a.  d.  ward , and a.  n.  burgess , phys . lett . * 88 * , 246102 ( 2002 ) .
Abstract:

Answer: we present a detailed analysis of the effective force between two smooth spherical colloids floating at a fluid interface due to deformations of the interface . the results hold in general and are applicable independently of the source of the deformation provided the capillary deformations are small so that a superposition approximation for the deformations is valid . we conclude that an effective long  ranged attraction is possible if the net force on the system does not vanish . otherwise , the interaction is short  ranged and can not be computed reliably based on the superposition approximation . as an application , we consider the case of like  charged , smooth nanoparticles and electrostatically induced capillary deformation . the resulting long  ranged capillary attraction can be easily tuned by a relatively small external electrostatic field , but it can not explain recent experimental observations of attraction if these experimental systems were indeed isolated .