Three Body Fluctuation-Induced Interaction at Fluid Interfaces:A Strong Deviation from the Pairwise Summation

Three Body Fluctuation-Induced Interaction at Fluid Interfaces: A Strong Deviation from the Pairwise Summation


We present a new method based on the scattering technique to investigate fluctuation-induced forces at a fluid interface. The scattering approach, well suited to the study of many body systems of arbitrary geometries, is augmented to include boundary fluctuations. Using this method, we study the deviation of the total fluctuation-induced interaction from the sum of pairwise energies for three colloidal particles. We consider both frozen and fluctuating colloids and obtain a very good agreement between analytical and numerical results. We find a marked difference in the three body fluctuation-induced free energy between the frozen and fluctuating colloids, both in sign and relative size.


With the construction of micro- and nano-eletro-mechanical devices (MEMS and NEMS) (1); (2); (3) and the desire to understand biological processes in greater details, it has become necessary to have a complete understanding of fluctuation-induced forces. These forces were first predicted by Casimir in 1948 who showed that confinement of zero point energy fluctuations of the electromagnetic field between two conducting plates in vacuum gives rise to long ranged attractive interactions (4); (5). Thirty years later, Fisher and de Gennes noted that imposing constraints on the thermal, rather than quantum, fluctuations in fluids in the vicinity of their critical points also results in similar long-ranged forces (6). Since then, the thermal versions of Casimir forces have been investigated in various theoretical and experimental studies (7); (8); (9); (10); (11); (12); (13); (14); (15); (16); (17); (18); (19).

In particular, due to the industrial applications and importance in biological systems, the thermal fluctuation-induced forces between colloidal particles have been studied within the last few years. For example, the experiments of Hertlein et. al. reveal the presence of repulsive and attractive fluctuation-induced forces on colloidal particles immersed in a near-critical binary mixture in the vicinity of a wall (11). The thermal fluctuation-induced force also plays an important role in the effective interaction between colloidal particles trapped between two fluid phases (20); (21).

There have been previous works using various techniques to study fluctuation-induced forces at interfaces with interesting results (12); (13); (14); (15); (16); (19); (17); (18). There are, however, a number of complications associated with each of these methods. The most common is the difficulty of calculating the interaction between colloids at very small separations and matching it with the Derjaguin approximation (22). Moreover, extending the previous methods to various geometries and many particle systems can be tedious.

In this work we study the interaction between colloids trapped at a fluid interface using the scattering formalism (23); (24); (25), which is widely used to obtain quantum fluctuation-induced forces. The scattering formalism presents a straightforward approach for treating various geometries and many body systems and is numerically faster compared to previous methods.

Furthermore, the scattering method simplifies the fluctuation-induced calculations by writing the energy into the translation matrices () and scattering matrices or -matrices (). The translation and -matrices describe how a fluctuation propagates through the field between the objects and how fluctuations interacts with the objects, respectively.

The important difference between the QED Casimir effect and the fluctuation-induced forces between colloids at an interface is that the colloids and their boundaries can fluctuate. This work extends the scattering method to include these effects, by separating out the interface and colloid Hamiltonians. This new extension can simplify the problem by separating the properties of the colloid, as well as the type of fluctuations, from the calculation of the -matrix. Note, in contrast to previous works (12); (13); (14); (15); (16); (19), the calculation of the T- and the translation matrices are completely decoupled.

In order to demonstrate the effectiveness of this new technique, we calculate the fluctuation-induced forces between three spherical Janus colloids at an interface. We study three different fluctuation scenarios: (i) colloids frozen at the interface; (ii) bobbing colloids that fluctuate only vertically; and (iii) bobbing and tilting colloids that both fluctuate vertically and tilt side to side. Surprisingly we find different three body behaviors for fixed and fluctuating colloids. The three body effect for the fixed colloids is repulsive and comparable to the two body interaction, while for fluctuating colloids is attractive and comparable only at short separations.

In the following, we briefly explain the derivation of the scattering formalism for colloids at fluid interfaces. Moreover, we show some numerical and analytical results obtained with this technique. Note that due to the rapid communications format, we only outline the important steps. The details of the derivations and calculations are left for a more complete exposition (26).

We consider an infinite interface between two fluid phases characterized by the surface tension . At equilibrium, the interface is flat and placed at . We suppose that the thermal height fluctuations of the interface are small without overhanges and bubbles. Therefore, we express the interface height profile in the Monge representation, . The capillary wave Hamiltonian (15); (16) then becomes


where is the capillary length with , the difference in the mass density of the two fluid phases.

We now introduce colloidal particles at the interface between two fluid phases. The colloid Hamiltonian will be defined as the energy costs associated with the insertion of a colloid,


where is the height of the center of mass of the colloid with respect to the equilibrium interface, and is the surface tension between the colloid and the fluid phase I(II). The contact line field is the height at which the fluid interface intersects the colloid, and is extended into the interior of the colloid. The cross sections of the colloid with the equilibrium and fluctuating interfaces are denoted by and , respectively, and is the change in the projected area. The surface area of the colloid in contact with the fluid phase I(II) is represented by where is the change in the area of the colloid in the fluid phase I. Note that for Janus particles the total area of the colloid does not change . In addition, the gravitational energy due to the colloid height is negligible and can be ignored in Eq. (2) because .

The first term in Eq. (2) is associated with excluding the projected areas from the interface since the integration domain in Eq. (1) covers the whole space. The second term represents the energy costs related to the change of the projected area for a colloid in the fluctuating vs reference interface, and the third term is the energy costs due to the change in the area of the colloids in the phase I and II.

The total Hamiltonian of the system is the sum of the interface and colloid Hamiltonians, , and the partition function is the sum over all interface and colloid configurations


with representing the constraints imposed by the colloidal particles on the interface height fluctuations. At the surface of the colloid the interface height field is constrained by the contact line field . A Dirac delta functional is used to implement the constraints, , where we have introduced the auxiliary field . Equation (3) can be separated as


where and correspond to the partition functions related to the interface and colloid fluctuations respectively.

The integration over the scalar field in the interface partition function can be performed and results in


where is a constant and


In Eq. (6), is the Green’s function of the capillary wave equation for the free interface and can be expanded in terms of the solutions to the capillary wave equation, i.e. . Using separation of variables in a coordinate system centered on a colloid, the capillary wave equation has two solutions: the solution regular at origin, which corresponds to the incident field in the scattering method ; and a solution that dies off at infinity, which corresponds to the scattered field . The expansion of the Green’s function and the auxiliary field in terms of these solutions are and , respectively, with the expansion coefficient and the multipole coefficient. Plugging these equations into Eq. (6), we find for and for . The translation matrix describes the couplings between the partial waves on distinct colloids. The scattering amplitude matrix relates the incident and scattered fields, and depends only on the shape of the colloid. Inserting the expressions for in Eq. (4), we find


The matrix () in Eq. (7) is a diagonal matrix with the elements that are the coefficient of the Green’s function expansion, and is a vector with elements that are the multipole moments associated with the particle .

We evaluate the colloid partition function following the same steps as the ones we took to calculate the interface partition function. The colloid Hamiltonian in Eq. (2) is expanded to be quadratic in the contact line field , and then the field is expanded in terms of the incident solutions,


After expansion the colloid partition function is


where the exact form of depends on the type of colloid fluctuations and will be presented below. We now insert Eqs. (9) and (7) into Eq. (4) to obtain the total partition function . The resulting partition function is a Gaussian integral over the multipoles and . The fluctuation-induced free energy at temperature can be then obtained through , with corresponding to the partition function with all the colloids placed at infinite distance from each other. The fluctuation-induced energy as a function of the -matrix and the translation matrix then reads as


where and with the modified T-matrix of the colloid given by


Note while depends only on the shape of the colloid, also depends on which in turn depends on the specific fluctuations that the colloid undergoes. Equations (10) and  (11) can be used to obtain the fluctuation-induced interaction between N colloidal particles at a fluid interface.

We now employ Eq. (10) to calculate the fluctuation-induced free energy between three identical spherical colloids with radius trapped at an interface. For simplicity we assume that the colloids are Janus particles, i.e. the contact line is pinned to the colloid surface such that the area of the colloidal particles in two fluid phases does not change, i.e. . In addition, we set the equilibrium contact angle to . Using the solutions to the Helmholtz equation in polar coordinates, the T-matrix for colloidal particles with circular cross section at the flat interface reads


where and are the modified Bessel functions of the first kind and second kind, respectively. In order to obtain the fluctuation-induced energy given by Eq. (10), it is necessary to calculate or more specifically the modified given in Eq. 11. To find we assume that the colloids are: (i) frozen, allowing no fluctuation; (ii) bobbing, allowing the colloids to fluctuate up and down; and (iii) bobbing and tilting, allowing the colloids to fluctuate both up and down as well as tilt side to side.

For a frozen colloid at the interface, the colloidal particles are exposed to an infinite potential barrier to fluctuate, therefore the elements of the matrix tend to infinity and from Eq. (11) . For bobbing colloidal particles, the projected area does not change and therefore in Eq. (2). Substituting Eq. (8) into Eq. (2), we can calculate the matrix . In the limit , we find for and for . From Eq. (11) we immediately find that for , and for , .

For colloids that both bob and tilt, the projected area changes with fluctuations, i.e. . After calculating the change in the projected area for small tilt angles, using Eq. (8) in the limit we find for , and respectively. Therefore, we have , for and otherwise.

To calculate the fluctuation-induced energy given in Eq. (10), we also need to calculate the translation matrix in polar coordinate, which can be obtained by using the Graff’s addition theorem for Bessel functions (27). We find with the angle between the coordinate systems fixed to the colloids and .

Employing the translation and -matrices described above, we write

Figure 1: (Color online) Three body effect for colloids on a line with and on an equilateral triangle with : The ratio of the three body energy to the absolute value of the two body energy vs. the surface-to-surface separation of the colloids 1 and 2, . The dashed and dotted lines represent the asymptotic results of Eqs. (15) and (16) for the colloids in the linear and triangular configurations, respectively. Blue squares (linear configuration) and diamonds (triangular configuration) denote the numerical results for the bobbing and tilting colloids. Red circles (linear configuration) and triangles (triangular configuration) show the numerical results for the bobbing only colloids. Inset: Three body effect for fixed colloids in the linear and triangular configurations. The dashed and dotted lines show the asymptotic result in Eq. (14). The black circles (linear configuration) and red squares (triangular configuration) represent the numerical results for the fixed colloids.

where the first term is sum over the two-body interactions and the second term is the free energy due to the three body effect. For three fixed colloids at the interface, the asymptotic three body effect in the limit reads


with , the radius of the colloid and the center-to-center distance between the colloids and . Note that the energy in Eq. (14) depends only on the distance between the colloids.

For three bobbing colloids at the interface, the asymptotic energy in the limit reads


Equation (15) is in agreement with the result in Refs. (17); (18) and shows that the three body asymptotic energy is much smaller than the two body interaction . Therefore, at large separations, for bobbing colloids the sum of pairwise energies gives a very good approximation to the fluctuation-induced energy.

For three bobbing and tilting colloids, the three body energy in the limit reads


In this case the two body interaction scales with , and as such the pairwise energy summation again provides a good approximation for the total fluctuation-induced energy.

Figure 1 depicts the ratio of the three body effects to the absolute value of the two body energies vs. the surface-to-surface separation, , for two different configurations: colloids sitting on a line () or at the vertices of an equilateral triangle (). As shown in the figure, there is a very good agreement between the asymptotic energies (dashed and dotted lines) given by Eqs. (14)-(16) and the numerical results (symbols) calculated by Eq. (10).

As illustrated in the figure, the three body effect is very small at large separations for bobbing only (red circles and triangles) and bobbing and tilting colloids (blue squares and diamonds) in both configurations. In contrast, the inset in Fig. 1 shows that the three body effect energy for fixed colloids is comparable to their two body interaction with the opposite sign,i.e., there is a strong deviation from the pairwise summation for the case of fixed colloids in both linear and triangular configurations, see the circles and squares in the inset. Quite interestingly, since the three body effect for fixed colloids is positive (repulsive), the total energy is smaller than the sum of pairwise added interaction energies, i.e. .

In summary, we have presented an effective method to study the fluctuation-induced forces at fluid interfaces using a modified scattering method, in which the separation of the interface and colloid free energies enables us to treat various interface and colloid fluctuations. The new approach has many advantages, such as the treatment of many colloid systems, quick numerical calculations, and the treatment of different colloid fluctuations.

As an example, we studied the three body interaction between frozen and fluctuating spherical Janus colloids and found a very interesting difference between them. For both frozen and fluctuating colloids the two body interactions are attractive. However, for frozen colloids the three body interactions are repulsive weakening the interaction compared to the pairwise summation, while for fluctuating colloids the three body interactions are attractive strengthening the interaction. In addition, we found that the three body interaction for the frozen colloids is non-negligible for the complete range of interactions, while the three body interactions for the fluctuating colloids are negligible for larger separations.

Due to the non-trivial behavior of these forces, a better knowledge of them will shed light on structure formation and crystallization phenomena at the fluid interfaces.

Authors would like to thank Mehran Kardar for useful discussions. This work was supported by the NSF through grants DMR-06-45668 (RZ).


  1. F. M. Serry, D. Walliser, and G. J. Maclay, Journal of Applied Physics 84, 2501 (1998)
  2. E. Buks and M. L. Roukes, Phys. Rev. B 63, 033402 (Jan 2001)
  3. H. B. Chan, V. A. Aksyuk, R. N. Kleiman, D. J. Bishop, and F. Capasso, Science 291, 1941 (2001)
  4. H. B. G. Casimir and D. Polder, Phys. Rev. 73, 360 (Feb 1948)
  5. H. B. G. Casimir, Proc. Koninklijke Nederlandse Akad. Wetenschappen B51, 793 (1948)
  6. M. E. Fisher and P. G. de Gennes, C. R. Acad. Sci. Paris B 287, 207 (1978)
  7. R. Zandi, J. Rudnick, and M. Kardar, Phys. Rev. Lett. 93, 155302 (Oct 2004)
  8. R. Garcia and M. H. W. Chan, Phys. Rev. Lett. 83, 1187 (Aug 1999)
  9. M. Goulian, R. Bruinsma, and P. Pincus, EPL (Europhysics Letters) 22, 145 (1993)
  10. M. Krech and S. Dietrich, Phys. Rev. A 46, 1886 (Aug 1992)
  11. C. Hertlein, L. Helden, A. Gambassi, S. Dietrich, and C. Bechinger, Nature 451, 172 (2008)
  12. R. Golestanian, M. Goulian, and M. Kardar, Phys. Rev. E 54, 6725 (Dec 1996)
  13. H. Lehle, E. Noruzifar, and M. Oettel, Eur. Phys. J. E 26, 151 (2008)
  14. E. Noruzifar and M. Oettel, Phys. Rev. E 79, 051401 (May 2009)
  15. H. Lehle, M. Oettel, and S. Dietrich, Europhys. Lett. 75, 174 (2006)
  16. H. Lehle and M. Oettel, Phys. Rev. E 75, 011602 (Jan 2007)
  17. C. Yolcu, I. Z. Rothstein, and M. Deserno, EPL (Europhysics Letters) 96, 20003 (2011)
  18. C. Yolcu, I. Z. Rothstein, and M. Deserno, Phys. Rev. E 85, 011140 (Jan 2012)
  19. H.-K. Lin, R. Zandi, U. Mohideen, and L. P. Pryadko, Phys. Rev. Lett. 107, 228104 (Nov 2011)
  20. F. Bresme and M. Oettel, Journal of Physics: Condensed Matter 19, 413101 (2007)
  21. M. Oettel and S. Dietrich, Langmuir 24, 1425 (2008)
  22. B. V. Derjaguin, Kolloid Z. 69, 155 (1934)
  23. T. Emig, N. Graham, R. L. Jaffe, and M. Kardar, Phys. Rev. Lett. 99, 170403 (Oct 2007)
  24. T. Emig, N. Graham, R. L. Jaffe, and M. Kardar, Phys. Rev. D 77, 025005 (Jan 2008)
  25. S. J. Rahi, T. Emig, N. Graham, R. L. Jaffe, and M. Kardar, Phys. Rev. D 80, 085021 (Oct 2009)
  26. E. Noruzifar, J. Wagner, and R. Zandi, in preparation
  27. “NIST Digital Library of Mathematical Functions,”, Release 1.0.5 of 2012-10-01, online companion to (28)
  28. NIST Handbook of Mathematical Functions, edited by F. W. J. Olver, D. W. Lozier, R. F. Boisvert, and C. W. Clark (Cambridge University Press, New York, NY, 2010) print companion to (27)
This is a comment super asjknd jkasnjk adsnkj
The feedback cannot be empty
Comments 0
The feedback cannot be empty
Add comment

You’re adding your first comment!
How to quickly get a good reply:
  • Offer a constructive comment on the author work.
  • Add helpful links to code implementation or project page.