## 1 Introduction

Swimming microorganisms do not live in an infinite, unbounded fluid domain, but instead inhabit complex geometries confined by fluid interfaces and solid boundaries, and populated by other organisms and passive particles. Much, if not most, of the rich variety of behaviour that is seen (Lauga & Powers Reference Lauga and Powers2009; Marchetti *et al.*
Reference Marchetti, Joanny, Ramaswamy, Liverpool, Prost, Rao and Simha2013) cannot be explained outside of the context of confinement. The most basic interaction is of a single swimmer with a boundary or object; even in such cases we see striking behaviour, such as the ‘waltzing’ of a pair of Volvox colonies (Drescher *et al.*
Reference Drescher, Leptos, Tuval, Ishikawa, Pedley and Goldstein2009). Flagellated bacteria such as *Escherichia coli* and *Vibrio alginolyticus* are known to trace out circles (Berg Reference Berg2000; DiLuzio Reference DiLuzio, Turner, Mayer, Garstecki, Weibel, Berg and Whitesides2005; Magariyama *et al.*
Reference Magariyama, Ichiba, Nakata, Baba, Ohtani, Kudo and Goto2005) near solid boundaries and, remarkably, if the boundary is replaced by a free surface the direction of rotation changes (Lauga *et al.*
Reference Lauga, Diluzio, Whitesides and Stone2006; Di Leonardo *et al.*
Reference Di Leonardo, Dell’Arciprete, Angelani and Iebba2011). An important effect is the attraction of swimmers to boundaries, noted by Rothschild (Reference Rothschild1963) for bull spermatozoa but subsequently also seen in bacteria (Frymier *et al.*
Reference Frymier, Ford, Berg and Cummings1995; Berke *et al.*
Reference Berke, Turner, Berg and Lauga2008), which is thought to contribute to the navigation of sperm cells in the female reproductive system (Denissenko *et al.*
Reference Denissenko, Kantsler, Smith and Kirkman-Brown2012) and plays a fundamental role in biofilm formation at surfaces (O’Toole, Kaplan & Kolter Reference O’Toole, Kaplan and Kolter2000). Swimmers adhere to surfaces due to strong lubrication forces, causing catalytic self-propelled rods to orbit colloidal spheres (Takagi *et al.*
Reference Takagi, Palacci, Braunschweig, Shelley and Zhang2014) and the microalga Chlamydomonas reinhardtii to orbit cylindrical posts until flagellar beating detaches it (Contino *et al.*
Reference Contino, Lushi, Tuval, Kantsler and Polin2015). As the number of interacting components increases so does the complexity of the behaviour, resulting in the phase behaviours seen in dense suspensions of bacteria, including long-range orientational order (Cisneros *et al.*
Reference Cisneros, Kessler, Ganguly and Goldstein2011), and the formation of large-scale turbulent structures (Dombrowski *et al.*
Reference Dombrowski, Cisneros, Chatkaew, Goldstein and Kessler2004; Dunkel *et al.*
Reference Dunkel, Heidenreich, Drescher, Wensink, Bär and Goldstein2013) and stable spiral vortices (Wioland *et al.*
Reference Wioland, Woodhouse, Dunkel, Kessler and Goldstein2013). The latter can only occur in confinement. An understanding of these phenomena is important for the design of microfluidic systems, such as devices to direct swimmers (Denissenko *et al.*
Reference Denissenko, Kantsler, Smith and Kirkman-Brown2012) and harness them for mass transport (Weibel *et al.*
Reference Weibel, Garstecki, Ryan, Diluzio, Mayer, Seto and Whitesides2005; Koumakis *et al.*
Reference Koumakis, Lepore, Maggi and Di Leonardo2013) and to extract mechanical work from their activity (Di Leonardo *et al.*
Reference Di Leonardo, Angelani, Dell’Arciprete, Ruocco, Iebba, Schippa, Conte, Mecarini, De Angelis and Di Fabrizio2010); and to construct biomimetic artificial swimmers (Dreyfus *et al.*
Reference Dreyfus, Baudry, Roper, Fermigier, Stone and Bibette2005; Golestanian, Liverpool & Adjari Reference Golestanian, Liverpool and Adjari2005; Paxton, Sen & Mallouk Reference Paxton, Sen and Mallouk2005; Howse *et al.*
Reference Howse, Jones, Ryan, Gough, Vafabakhsh and Golestanian2007) which may fulfil a number of nanotechnological and medical roles.

The interaction of swimming microorganisms with each other and their environment is a complex combination of several factors, including biology, such as the taxis which allows them to move in search of food and tolerable living conditions and to aggregate and form patterns through chemical signalling and quorum sensing (Pedley & Kessler Reference Pedley and Kessler1992; Cates *et al.*
Reference Cates, Marenduzzo, Pagonabarraga and Tailleur2010); and hydrodynamics and physical contact. It has been shown that the scattering of swimmers from planes (Kantsler *et al.*
Reference Kantsler, Dunkel, Polin and Goldstein2013) and posts (Contino *et al.*
Reference Contino, Lushi, Tuval, Kantsler and Polin2015) arises from physical contact of the flagella with the surface, so that hydrodynamics is not the dominant contributor to the phenomena in these cases. Nonetheless, the fact that there is physical contact with the surface emphasises that any hydrodynamic effects have to be considered in this contact regime. Other cases are less clear-cut; for instance, the typical density profile of a suspension of swimmers close to a wall has been reproduced both by considering hydrodynamics (Berke *et al.*
Reference Berke, Turner, Berg and Lauga2008) and Brownian motion combined with collisions (Li & Tang Reference Li and Tang2009).

The importance of hydrodynamics to the interactions of swimmers means exact solutions to the Stokes equations are desirable. In the case of a single swimmer in an unbounded domain several such solutions exist, notably for the motion of a single axisymmetric squirmer (Lighthill Reference Lighthill1952; Blake Reference Blake1971*b*
), later generalised to non-axisymmetric slip velocities (Pak & Lauga Reference Pak and Lauga2014), and for the motion of a ‘treadmilling’ spheroidal (Leshansky *et al.*
Reference Leshansky, Kenneth, Gat and Avron2007) or toroidal swimmers (Taylor Reference Taylor1952; Purcell Reference Purcell1977; Leshansky & Kenneth Reference Leshansky and Kenneth2008), as well as a two-dimensional analogue for a squirming cylinder (Blake Reference Blake1971*a*
) or waving sheet (Taylor Reference Taylor1951). The squirmer solutions have been used to find the advection of tracer particles due to a squirmer (Pushkin, Shum & Yeomans Reference Pushkin, Shum and Yeomans2013). Dropping instead to two dimensions, a number of additional solutions in confinement become available using conformal mapping techniques, such as the motion of an active cylinder near a planar or concave boundary (Crowdy Reference Crowdy2011, Reference Crowdy2013; Papavassiliou & Alexander Reference Papavassiliou and Alexander2015) or under a free surface (Crowdy *et al.*
Reference Crowdy, Lee, Samson, Lauga and Hosoi2011). However, in three dimensions, hydrodynamic interactions have only been calculated by approximate methods. For instance, Ishikawa, Simmonds & Pedley (Reference Ishikawa, Simmonds and Pedley2006) find the far-field interactions between two squirmers by considering multipole expansions of stresslets, and the near-contact interactions using lubrication theory. In many cases point-singularity methods are valuable. The image systems calculated by Blake & Chwang (Reference Blake and Chwang1974) allow the interaction with walls to be found (Zargar, Najafi & Miri Reference Zargar, Najafi and Miri2009; Spagnolie & Lauga Reference Spagnolie and Lauga2012), recovering the experimentally observed attraction to walls (Berke *et al.*
Reference Berke, Turner, Berg and Lauga2008) and the swimming in circles close to surfaces (Lauga *et al.*
Reference Lauga, Diluzio, Whitesides and Stone2006; Papavassiliou & Alexander Reference Papavassiliou and Alexander2015). Other phenomena relevant to low Reynolds number swimming that have been successfully treated using point singularities include flagellar beating for feeding (Higdon Reference Higdon1979) and resulting in synchronisation (Brumley *et al.*
Reference Brumley, Wan, Polin and Goldstein2014). A comparison of real flow fields with point-singularity approximations shows them to be in good agreement even close to the organism in several cases (Drescher *et al.*
Reference Drescher, Goldstein, Michel, Polin and Tuval2010, Reference Drescher, Dunkel, Cisneros, Ganguly and Goldstein2011).

We find exact solutions for the axisymmetric translation and rotation of a squirmer in the presence of a spherical or planar boundary. These are valid at any separation, both in the far field where point singularity solutions are accurate and in the contact limit of vanishing separation, where such approximate solutions are not accurate. They also account for any type of squirming motion and not simply the lowest-order modes considered by Ishikawa *et al.* (Reference Ishikawa, Simmonds and Pedley2006) and that point-singularity descriptions are restricted to. These solutions are obtained using the Lorentz reciprocal theorem for Stokes flows (Happel & Brenner Reference Happel and Brenner1983), first applied to calculate the motion of individual swimmers by Stone & Samuel (Reference Stone and Samuel1996) and recently extended to a many-body setting (Papavassiliou & Alexander Reference Papavassiliou and Alexander2015). In this form of the reciprocal theorem, the stress tensor associated with the Stokes drag on an object of the same shape as the swimmer serves as the integration kernel to extract the speed and angular frequency of the swimming from the slip velocity. This can be viewed as a specialisation of the boundary-element method (Pozrikidis Reference Pozrikidis1992), simplified by the requirement that a swimmer be free of net forces and torques. The simplicity of the Stokes drag solution on a sphere means this calculation is straightforward for a single spherical microorganism, although since the full hydrodynamic solution of a squirming sphere exists (Lighthill Reference Lighthill1952; Blake Reference Blake1971*b*
; Pak & Lauga Reference Pak and Lauga2014) the only advantage of the reciprocal theorem is computational convenience. Indeed, the swimming speed found for an active sphere self-propelling by means of a metachronal wave on its surface by Stone & Samuel (Reference Stone and Samuel1996) had been derived by other means not long previously by Ehlers *et al.* (Reference Ehlers, Samuel, Berg and Montgomery1996). Nevertheless the simplicity of the calculation means it has become a standard tool in the active matter literature (Squires & Bazant Reference Squires and Bazant2004; Golestanian, Liverpool & Adjari Reference Golestanian, Liverpool and Adjari2007). More recently there have been interesting extensions of the reciprocal theorem and other related integral theorems to cases such as propulsion by the Marangoni effect (Masoud & Stone Reference Masoud and Stone2014) and self-propulsion through viscoelastic and non-Newtonian fluids (Lauga Reference Lauga2014) where direct solutions are not so readily available.

The reciprocal theorem is immediately applicable to swimmers in confined fluid domains. This has been exploited in two dimensions to find the motion of squirming (Crowdy Reference Crowdy2011) and self-diffusiophoretic cylinders close to walls (Crowdy Reference Crowdy2013) using as a conjugate solution the Stokes drag on a cylinder in the half-space (Jeffrey & Onishi Reference Jeffrey and Onishi1981), and it has been noted that the reciprocal theorem may be used to find the motion of any number of swimmers (Papavassiliou & Alexander Reference Papavassiliou and Alexander2015). More generally, the reciprocal theorem may be used to find the full hydrodynamics for a given active problem, provided the existence of an appropriate conjugate solution: if the Green’s function for the Stokes equations in a particular confined geometry is known, the reciprocal theorem extracts the flow due to activity on the boundaries. Following such an approach, Michelin & Lauga (Reference Michelin and Lauga2015*b*
) found the fluid flux through an active pipe by relation to the flow solution for a no-slip channel. Furthermore, an approximate integration kernel for a swimmer in a given geometry may be constructed using one of the many existing flows for point forces; relevant examples include the solution for point singularities near walls (Blake & Chwang Reference Blake and Chwang1974), outside spheres (Higdon Reference Higdon1979), between two plates (Liron & Mochon Reference Liron and Mochon1976) or in a cylindrical pipe (Liron & Shahar Reference Liron and Shahar1978).

By following this approach, we are able to find exact solutions for swimmer interactions by using exact solutions for the conjugate Stokes drag problem. Such solutions are available for the Stokes drag on a pair of spheres, or of a single sphere close to a planar wall or fluid interface. The symmetries of the geometry mean there are two independent directions, namely the common diameter of the two spheres and any axis perpendicular to this, and the solution consists of translation and rotation in each of these directions, so that the general motion separates into four components that can be treated individually. The axisymmetric rotation was solved exactly by Jeffery (Reference Jeffery1915), and has since been supplemented by the closely related solution for rotation of a sphere beneath a planar fluid interface (Brenner Reference Brenner1964; O’Neill & Ranger Reference O’Neill and Ranger1979). The solution for axisymmetric translation was given by Stimson & Jeffery (Reference Stimson and Jeffery1926) and found to be in remarkably good agreement with experiment (Happel & Pfeffer Reference Happel and Pfeffer1960). The special case of sedimentation of a sphere towards a solid plane was subsequently given a more detailed analysis both in the limit of large separation (Brenner Reference Brenner1961) and of contact (Cox & Brenner Reference Cox and Brenner1967), with the latter giving a comparison to results obtained from lubrication theory. Several attempts to find the non-axisymmetric motions and rotations have been unable to give the solution in a closed form; the problem is reduced to a system of difference equations, of which an analytic solution has not been found. Nevertheless it is possible to compute the flow to any degree of accuracy (Dean & O’Neill Reference Dean and O’Neill1963; O’Neill Reference O’Neill1964; Goldman, Cox & Brenner Reference Goldman, Cox and Brenner1966, Reference Goldman, Cox and Brenner1967*a*
,Reference Goldman, Cox and Brenner
*b*
; O’Neill & Stewartson Reference O’Neill and Stewartson1967; O’Neill & Majumdar Reference O’Neill and Majumdar1970*a*
).

All of these results rely on the use of bispherical coordinates, in which any configuration of two convex or concave spherical boundaries, as well as the intermediate limit of a plane, is an isosurface. This coordinate system greatly simplifies the imposition of boundary conditions; furthermore, since it is conformally equivalent to spherical coordinates, Laplace’s equation is separable (Jeffery Reference Jeffery1912), allowing for a general solution to the Stokes equations (Jeffery Reference Jeffery1922). Another notable application of Stimson & Jeffery’s solution to swimmer problems is to study the hydrodynamics of catalytic dimers. Catalytic dimers are artificial self-propelled particles composed of a pair of chemically active colloidal beads powered by self-diffusiophoresis (Rückner & Kapral Reference Rückner and Kapral2007); their simplicity facilitates manufacture and allows experiments involving many interacting units (Valadares *et al.*
Reference Valadares, Tao, Zacharia, Kitaev, Galembeck, Kapral and Ozin2010). Using the reciprocal theorem together with bispherical coordinates, Popescu, Tasinkevych & Dietrich (Reference Popescu, Tasinkevych and Dietrich2011) and Michelin & Lauga (Reference Michelin and Lauga2015*a*
) calculated exact expressions for the propulsion speed of catalytic dimers, and were able to discuss optimisation of their swimming speed through changes to the relative sizes and separation of the two beads. Bispherical coordinates also provide a way to calculate hydrodynamic interactions of two spherical objects, such as a sphere sedimenting against a plane (Cox & Brenner Reference Cox and Brenner1967); here, we investigate their use in calculating interactions driven by force-free swimming, for which the reciprocal theorem is ideally suited (Papavassiliou & Alexander Reference Papavassiliou and Alexander2015). In this way Mozaffari *et al.* (Reference Mozaffari, Sharifi-Mood, Koplik and Maldarelli2016) and Sharifi-Mood, Mozaffari & Córdova-Figueroa (Reference Sharifi-Mood, Mozaffari and Córdova-Figueroa2016) found the interaction of spherical self-diffusiophoretic particles with each other and with planar boundaries, finding that the chemical interaction is dominant over the hydrodynamic and usually results in repulsion, except where coverage of the chemically active site over the swimmers is large. Their consideration of non-axisymmetric components of motion necessitated numerical solution.

We outline the equations of viscous flow and discuss the use of the reciprocal theorem to obtain exact solutions for interactions in § 2. In § 3 we review the Stokes drag solutions of Jeffery (Reference Jeffery1915) and Stimson & Jeffery (Reference Stimson and Jeffery1926) that we use with the reciprocal theorem. This section is included for reference and may be skipped if desired. The main results of this paper are contained in § 4, where we find the motion of a squirmer interacting with a passive spherical boundary. The contribution to the motion from the azimuthal squirming coefficients is found explicitly for all orders and is shown for a model organism driven by a rotating cap, while a simple extension to a case of interaction with a planar free surface is discussed in § 4.4. The meridional and radial squirming coefficients do not, at present, have a general form for the interaction valid at all orders; in § 4.5 we calculate the interaction due to the lowest two orders of these modes. Finally we discuss the results obtained and possible extensions to the work presented here in § 5.

## 2 Stokes flows and the reciprocal theorem

The motions of a collection of swimmers or active particles can be determined from the Lorentz reciprocal theorem for Stokes flows. Specifically, for a collection of $N$ force and torque-free swimmers generating motion through active surface slip velocities $\boldsymbol{u}_{i}^{s}$ on their boundaries $\unicode[STIX]{x2202}D_{i}$ , $i=1,\ldots ,N$ , their translational speeds $\tilde{\boldsymbol{U}}_{i}$ and rotations $\tilde{\unicode[STIX]{x1D734}}_{i}$ are given by (Stone & Samuel Reference Stone and Samuel1996; Papavassiliou & Alexander Reference Papavassiliou and Alexander2015)

Here, $\hat{\boldsymbol{n}}$ is the unit outward normal to the fluid domain and $\boldsymbol{\unicode[STIX]{x1D70E}}$ is the stress tensor of a conjugate Stokes flow solution for the same set of particles acted upon by forces $\boldsymbol{F}_{i}$ and torques $\boldsymbol{T}_{i}$ , and with no-slip boundary conditions. Thus the fundamental object in application of the reciprocal theorem to swimmer problems is the normal stress of the Stokes drag problem,

where $p$ is the pressure and $\unicode[STIX]{x1D707}$ the viscosity, whose integral against the slip velocities yields the swimmer motion. For instance, in the case of a single spherical swimmer, the classic flow solution for the Stokes drag on a single sphere may be used to calculate the swimming speed of a squirmer (Stone & Samuel Reference Stone and Samuel1996).

Of course, for a spherical squirmer the full flow field, in addition to the swimmer motion, may be calculated directly (Lighthill Reference Lighthill1952; Blake Reference Blake1971*b*
; Pak & Lauga Reference Pak and Lauga2014), without significant additional effort, and gives more information. In the case of hydrodynamic interactions between two objects, even as simple as two squirming spheres, a direct solution for the full flow is not available.

However, full solutions are available for the Stokes drag of two spheres under conditions of axisymmetry, both for rotational and translational motion. These may be used in the reciprocal theorem to deduce the corresponding axisymmetric interactions of an arbitrary pair of squirming spheres, or of a single squirming sphere with a spherical, or planar, boundary. The solution is founded upon having an expression for the stress tensor for a conjugate Stokes drag problem. The work of Jeffery (Reference Jeffery1915) and Stimson & Jeffery (Reference Stimson and Jeffery1926), as well as subsequent extensions and generalisations (Payne & Pell Reference Payne and Pell1960; Brenner Reference Brenner1961; Kanwal Reference Kanwal1961; Dean & O’Neill Reference Dean and O’Neill1963; Brenner Reference Brenner1964; O’Neill Reference O’Neill1964; Cox & Brenner Reference Cox and Brenner1967; O’Neill & Stewartson Reference O’Neill and Stewartson1967; O’Neill & Majumdar Reference O’Neill and Majumdar1970*a*
,Reference O’Neill and Majumdar
*b*
; Majumdar & O’Neill Reference Majumdar and O’Neill1977), gives the flow for these problems, from which the stress can be computed directly. However, to keep our work self-contained and in a consistent notation, we rederive the Stokes drag solutions ourselves. This requires relatively little work for the rotation, while for the translation we feel that our solution offers some modest improvements on the original calculation of Stimson & Jeffery.

It should be noted that the reciprocal theorem may also be used to calculate approximate expressions for the motion of a squirmer near a wall (Davis & Crowdy Reference Davis and Crowdy2015; Papavassiliou & Alexander Reference Papavassiliou and Alexander2015): since the leading-order flow about a sedimenting sphere is a Stokeslet, the well-known solution for a Stokeslet and for a rotlet near a wall (Blake & Chwang Reference Blake and Chwang1974) may be used to construct stress tensors which can obtain the translational and rotational motion respectively, to third and fourth order in the swimmer’s size; the approximation is improved by also including the stress tensor derived from Blake’s solution for the source dipole near a wall. The results given are identical to those found by matched asymptotics (Davis & Crowdy Reference Davis and Crowdy2015), but by approximating the integration kernel rather than the swimming stroke, the slip velocity may be kept completely general (Papavassiliou & Alexander Reference Papavassiliou and Alexander2015). This approach is obviously extensible to other geometries: for instance Higdon (Reference Higdon1979) has given the solution for a Stokeslet outside a sphere, which would allow the interactions between two swimmers to be determined approximately.

## 3 Bispherical coordinates and conjugate solutions

We record in this section the solutions to Stokes drag problems involving two spheres that we will use with the reciprocal theorem to obtain exact swimmer hydrodynamics. The reader who is primarily interested in these applications may safely skip to § 4 and only refer back as necessary.

### 3.1 Bispherical coordinates

Problems involving two spheres, such as we consider here, are naturally treated by employing a bispherical coordinate system. If $z+\text{i}\unicode[STIX]{x1D70C}$ is a complex coordinate on a Cartesian grid, the bipolar coordinate grid $\unicode[STIX]{x1D709}+\text{i}\unicode[STIX]{x1D702}$ is defined by

*a*,

*b*) $$\begin{eqnarray}\unicode[STIX]{x1D709}+\text{i}\unicode[STIX]{x1D702}=\ln \left[{\displaystyle \frac{z+\text{i}\unicode[STIX]{x1D70C}-R}{z+\text{i}\unicode[STIX]{x1D70C}+R}}\right],\quad z+\text{i}\unicode[STIX]{x1D70C}=-R{\displaystyle \frac{(\sinh \unicode[STIX]{x1D709}-\text{i}\sin \unicode[STIX]{x1D702})}{(\cosh \unicode[STIX]{x1D709}-\cos \unicode[STIX]{x1D702})}},\end{eqnarray}$$

where
$R$
is a positive real number. This can be thought of as a stereographic projection of the lines of latitude and longitude on a sphere about a point on the equator, as demonstrated in figure 1(*a*), with the poles mapping to two symmetric points,
$z=\pm R$
. Finally, a rotation about the
$z$
axis gives an azimuthal coordinate
$\unicode[STIX]{x1D719}$
, which coincides for bispherical and ordinary cylindrical coordinates. Surfaces of constant
$\unicode[STIX]{x1D709}$
are non-intersecting spheres centred on
$-R\,\text{coth}\,\unicode[STIX]{x1D709}$
, with radius
$r=R\,|\text{cosech}\,\unicode[STIX]{x1D709}|$
. We consider the fluid to be the region
$\unicode[STIX]{x1D709}_{2}<\unicode[STIX]{x1D709}<\unicode[STIX]{x1D709}_{1}$
, where
$\unicode[STIX]{x1D709}_{1}$
is taken to be positive. The choice of
$\unicode[STIX]{x1D709}_{2}$
then defines the geometry: if it is positive the fluid is the region between two nested spherical boundaries and if it is negative the fluid is external to the two spheres, while the intermediate case
$\unicode[STIX]{x1D709}_{2}=0$
represents the half-space, as shown in figure 1(*c*). In what follows we will use both cylindrical coordinates
$(z,\unicode[STIX]{x1D70C},\unicode[STIX]{x1D719})$
and bispherical coordinates
$(\unicode[STIX]{x1D709},\unicode[STIX]{x1D702},\unicode[STIX]{x1D719})$
, and denote by
$W$
the conformal factor
$R/(\cosh \unicode[STIX]{x1D709}-\cos \unicode[STIX]{x1D702})$
that appears frequently.

### 3.2 Coaxial rotation

The general solution for axisymmetric azimuthal flows was given by Jeffery (Reference Jeffery1915), together with a number of specific examples, including the coaxial rotation of two spheres, and was subsequently expanded upon by Kanwal (Reference Kanwal1961) to consider objects which do not intersect the axis of symmetry. For an axisymmetric, purely azimuthal Stokes flow, the fluid velocity takes the form $\boldsymbol{u}=u_{\unicode[STIX]{x1D719}}(z,\unicode[STIX]{x1D70C})\boldsymbol{e}_{\unicode[STIX]{x1D719}}$ . The pressure is constant everywhere and may be taken to be equal to zero without loss of generality, so the flow satisfies the scalar equation $(\unicode[STIX]{x1D6FB}^{2}-\unicode[STIX]{x1D70C}^{-2})u_{\unicode[STIX]{x1D719}}=0$ , of which the general solution in bispherical coordinates is

The constants $c_{l}$ and $d_{l}$ are determined from the boundary conditions of solid-body rotation of the sphere $\unicode[STIX]{x1D709}=\unicode[STIX]{x1D709}_{1,2}$ with angular velocity $\unicode[STIX]{x1D714}_{1,2}$ about their common diameter

The force per area on the spheres’ surfaces is purely azimuthal, and integrating it gives the torques

These expressions are uniformly convergent since they are bounded independently of $\unicode[STIX]{x1D709}_{1,2}$ : taking the first term in (3.5) as an example, using that $\sinh n\unicode[STIX]{x1D709}>n\sinh \unicode[STIX]{x1D709}$ for $n>0$ and $\unicode[STIX]{x1D709}>0$ we have that

for $\unicode[STIX]{x1D709}_{2}<0,\,\unicode[STIX]{x1D709}_{1}>0$ . The uniform convergence of the remaining terms is demonstrated similarly, although note that (3.6) diverges as $\unicode[STIX]{x1D709}_{2}\rightarrow \unicode[STIX]{x1D709}_{1}$ .

A limit that is of particular interest for us in considering the near-field hydrodynamics of swimmers is that of vanishing separation where the spheres touch. In this limit $\unicode[STIX]{x1D709}_{1}$ and $\unicode[STIX]{x1D709}_{2}$ both tend to zero in a way that preserves the ratio $r_{1}/r_{2}$ ,

Then, since (3.5) and (3.6) are uniformly convergent we may interchange the limit taking with the summation, finding that the torques converge to the finite values

where $\unicode[STIX]{x1D701}(s,q)\equiv \sum _{n=0}^{\infty }(n+q)^{-s}$ is the Hurwitz zeta function and $\unicode[STIX]{x1D701}(s,1)\equiv \unicode[STIX]{x1D701}(s)$ the Riemann zeta function. When one sphere encloses the other ( $\unicode[STIX]{x1D709}_{2}>0$ ) the torques on the two spheres are equal and opposite, meaning that only relative motion can be deduced from the reciprocal theorem, the left-hand side of (2.1) reducing to $T_{1}(\tilde{\unicode[STIX]{x1D6FA}}_{1}-\tilde{\unicode[STIX]{x1D6FA}}_{2})$ . It is natural to take the concave boundary to set the frame of reference. In the intermediate limit of a plane ( $\unicode[STIX]{x1D709}_{2}=0$ ) the two solutions (3.9) and (3.10) coincide and as with the enclosed system the torque on the wall is equal and opposite to the torque on the finite-sized sphere.

### 3.3 Translation along the common axis

The translational motion of two spheres along their common diameter was first studied by Stimson & Jeffery (Reference Stimson and Jeffery1926) and subsequently adapted for the special case of a sphere sedimenting towards a planar surface (Brenner Reference Brenner1961; Cox & Brenner Reference Cox and Brenner1967). The approach adopted involves the introduction of a streamfunction to solve the continuity equation by construction and then the Stokes equations reduce to a fourth-order operator acting on the streamfunction. Later studies (Dean & O’Neill Reference Dean and O’Neill1963; O’Neill Reference O’Neill1964; O’Neill & Stewartson Reference O’Neill and Stewartson1967; O’Neill & Majumdar Reference O’Neill and Majumdar1970*a*
,Reference O’Neill and Majumdar
*b*
; Majumdar & O’Neill Reference Majumdar and O’Neill1977), motivated in part by a desire to extend to translations perpendicular to the common axis, follow the opposite approach: the Stokes equation is first solved by constructing a harmonic vector from an appropriate combination of the flow and the pressure (
$\boldsymbol{u}-(1/2\unicode[STIX]{x1D707})p\boldsymbol{x}$
), the coefficients of which are then to be determined from boundary conditions and the imposition of incompressibility, resulting in a set of second-order difference equations that unfortunately proves analytically intractable. We quote the solution here in a form that is a minor variation of that given by Stimson & Jeffery (Reference Stimson and Jeffery1926). We also generalise to arbitrary translational speeds
$V_{1}$
and
$V_{2}$
and to the full range of geometries allowed by the bispherical coordinate system.

The domain has non-trivial cohomology in degree 2, associated with expansions or contractions of each of the two spherical boundary surfaces $\unicode[STIX]{x1D709}=\unicode[STIX]{x1D709}_{1,2}$ , such as might occur for two small gas bubbles. Neglecting these motions, the fluid velocity can be written as the curl of a vector $\unicode[STIX]{x1D74D}$ that satisfies the biharmonic equation, $\unicode[STIX]{x1D6FB}^{4}\unicode[STIX]{x1D74D}=\pmb{{\it 0}}$ . For an axisymmetric flow, $\unicode[STIX]{x1D74D}$ may be chosen to be purely azimuthal, $\unicode[STIX]{x1D74D}=\unicode[STIX]{x1D713}\boldsymbol{e}_{\unicode[STIX]{x1D719}}$ and then the general solution may be written in the form

The boundary conditions that determine the real constants ${\mathcal{A}}_{l},{\mathcal{B}}_{l},{\mathcal{C}}_{l},{\mathcal{D}}_{l}$ are that the axial flow should equal the constant translation speeds $V_{1,2}\,\boldsymbol{e}_{z}$ of the two spheres $\unicode[STIX]{x1D709}=\unicode[STIX]{x1D709}_{1,2}$ and that the radial flow $u_{\unicode[STIX]{x1D70C}}$ should vanish on their surfaces. These conditions may be combined into the statement $\unicode[STIX]{x1D735}(\unicode[STIX]{x1D70C}\unicode[STIX]{x1D713}-(1/2)\unicode[STIX]{x1D70C}^{2}V_{j})_{\unicode[STIX]{x1D709}=\unicode[STIX]{x1D709}_{j}}=\pmb{{\it 0}}$ for $j=1,2$ (Stimson & Jeffery Reference Stimson and Jeffery1926), giving four equations to determine the four unknown coefficients, that reduce to the $2\times 2$ block-diagonal form

and

Inversion is straightforward, and the explicit forms of the coefficients ${\mathcal{A}}_{l},{\mathcal{B}}_{l},{\mathcal{C}}_{l},{\mathcal{D}}_{l}$ are shown in appendix A.

In bispherical coordinates, the flow is given by

which allows the pressure at a point $\boldsymbol{x}$ to be calculated as

where $p_{\infty }$ is its asymptotic value. There is no dependence on the path of integration as the domain is simply connected.

The solution we have presented is equivalent to those given previously (Stimson & Jeffery Reference Stimson and Jeffery1926; Brenner Reference Brenner1961), although it is not identical because the manner in which we have solved the continuity equation differs slightly. A consequence is the expansion of the vector potential in associated Legendre polynomials $P_{l}^{1}$ , rather than in Gegenbauer polynomials $C_{n+1}^{-1/2}$ . Furthermore the final scheme we arrive at for determining the coefficients ${\mathcal{A}}_{l},{\mathcal{B}}_{l},{\mathcal{C}}_{l},{\mathcal{D}}_{l}$ presents the $4\times 4$ problem in block diagonal form, which we have not seen in the previous literature.

Finally, using (3.14) and (3.15) a stress tensor may be constructed and integrated by parts over the spheres to give the hydrodynamic drag force

When the fluid has a finite volume or in the limiting case of the half-space ( $\unicode[STIX]{x1D709}_{2}\geqslant 0$ ), the net force on the fluid is zero since the contributions from the two boundaries are equal and opposite. This is also seen for a point force in the half-space (Blake & Chwang Reference Blake and Chwang1974), which is obtained from our solution in the limit $\unicode[STIX]{x1D709}_{1}\rightarrow \infty$ , $\unicode[STIX]{x1D709}_{2}\rightarrow 0$ , with $R$ held constant.

Since the coefficients ${\mathcal{A}}_{l},{\mathcal{B}}_{l},{\mathcal{C}}_{l},{\mathcal{D}}_{l}$ have exponential decay in $\unicode[STIX]{x1D709}_{1}$ and $\unicode[STIX]{x1D709}_{2}$ the sums converge rapidly for large separation; however, as the separation between the spheres vanishes the forces diverge. This limit has been treated in detail by Cox & Brenner (Reference Cox and Brenner1967) and we adapt their method here.

## 4 Swimmer interactions

### 4.1 Squirming

The results of the previous section for Stokes drag of two spheres allow a variety of axisymmetric swimmer motions to be determined via the reciprocal theorem. Despite the absence of expressions for non-axisymmetric motions this is enough to, for instance, give an exact description of the circular motion of microorganisms such as *E. coli* close to planar boundaries (Berg Reference Berg2000; Lauga *et al.*
Reference Lauga, Diluzio, Whitesides and Stone2006), and can also shed light on the hydrodynamics of a daughter colony of *Volvox* inside its parent (Drescher *et al.*
Reference Drescher, Leptos, Tuval, Ishikawa, Pedley and Goldstein2009), or the contact interaction of swimmers with passive particles (Wu & Libchaber Reference Wu and Libchaber2000). Many of these have been studied asymptotically using leading-order point-singularity descriptions (Berke *et al.*
Reference Berke, Turner, Berg and Lauga2008; Spagnolie & Lauga Reference Spagnolie and Lauga2012), but using the exact solutions for Stokes drag we are able to describe the behaviour for arbitrarily small separation and arbitrary squirming motions. Since the reciprocal theorem and the Stokes equations are linear, it suffices to calculate the interaction of a swimmer and passive sphere (Ishikawa *et al.*
Reference Ishikawa, Simmonds and Pedley2006), whose motion is given by

where the stress tensor is the corresponding to the Stokes drag problems given in § 3.

We consider a swimmer of radius $r_{1}$ centred a perpendicular distance $d$ away from the surface of a passive sphere of radius $r_{2}$ , which may be convex, flat or concave and we respectively term the ‘tracer’, ‘wall’ or ‘shell’. The swimmer approaches the surface at an angle $\unicode[STIX]{x1D6FC}$ and its squirming motion is described in terms of a local orthonormal basis $\{\boldsymbol{s}_{r},\boldsymbol{s}_{\unicode[STIX]{x1D703}},\boldsymbol{s}_{\unicode[STIX]{x1D719}}\}$ and polar coordinate system $(\unicode[STIX]{x1D703}_{s},\unicode[STIX]{x1D719}_{s})$ relative to this direction, as shown in figure 2. Its slip velocity may be decomposed into squirming modes (Lighthill Reference Lighthill1952) as

where
$V_{n}(x)\equiv -2P_{n}^{1}(x)/n(n+1)$
and
$A_{n}$
,
$B_{n}$
and
$C_{n}$
are real coefficients with the units of velocity. The free swimming speed, asymptotically far from the surface, is
$U_{free}=(2B_{1}-A_{1})/3$
(Lighthill Reference Lighthill1952; Blake Reference Blake1971*b*
). The addition of the azimuthal modes
$C_{n}$
(Pak & Lauga Reference Pak and Lauga2014) allows for axial rotation of the swimmer, as seen in several real microorganisms, with rotation speed
$\unicode[STIX]{x1D6FA}_{free}=-C_{1}$
about the axisymmetry axis.

To perform the integral in (4.1) it is convenient to express the swimmer’s slip velocity in bispherical coordinates. Defining the angle $\unicode[STIX]{x1D6FD}$ as the rotation angle about $\boldsymbol{e}_{\unicode[STIX]{x1D719}}$ between the cylindrical ( $\boldsymbol{e}_{z},\boldsymbol{e}_{\unicode[STIX]{x1D70C}}$ ) and bispherical ( $\boldsymbol{e}_{\unicode[STIX]{x1D709}},\boldsymbol{e}_{\unicode[STIX]{x1D702}}$ ) basis vectors,

the swimmer polar angle is found using the spherical cosine law to be,

as illustrated in figure 2, and the surface unit vectors are given by the transformation

The slip velocity involves Legendre polynomials in $\cos \unicode[STIX]{x1D703}_{s}$ , which are expanded using the addition theorem for Legendre functions (Sneddon Reference Sneddon1956; Maleček & Nádeník Reference Maleček and Nádeník2001),

As the stress tensor in (4.1) is axisymmetric the integral over $\unicode[STIX]{x1D719}$ only affects the slip velocity components. Hence it is convenient to perform the $\unicode[STIX]{x1D719}$ -integral first and define a vector of the resulting azimuthally averaged slip velocity components,

which is contracted against the stress tensor and integrated over $\unicode[STIX]{x1D702}$ to give the motion. The radial and meridional modes $A_{n}$ and $B_{n}$ cannot drive axisymmetric rotation, since the normal stress corresponding to axisymmetric rotation is purely azimuthal; similarly, the azimuthal modes $C_{n}$ cannot drive axisymmetric translation. The dependence of the motion on the swimmer’s orientation $\unicode[STIX]{x1D6FC}$ is a purely geometric factor for each order of squirming mode, and at large separations where higher-order modes may be neglected the orientation dependence is simply $P_{2}(\cos \unicode[STIX]{x1D6FC})$ , as found using point-singularity models of swimmer interactions with walls (Spagnolie & Lauga Reference Spagnolie and Lauga2012; Davis & Crowdy Reference Davis and Crowdy2015; Papavassiliou & Alexander Reference Papavassiliou and Alexander2015).

The contributions from the tangential modes, $B_{n},C_{n}$ , are evaluated straightforwardly (albeit tediously) using orthogonality of Legendre polynomials. The radial modes, $A_{n}$ , pick up a contribution from the pressure, which may be rewritten in terms of the flow by integrating by parts using the identity

### 4.2 Rotation

In this section we calculate explicitly the rotational motion of a squirmer close to a surface. Combined with self-propulsion parallel to the surface this rotation results in circling behaviour, which has been observed experimentally for flagellated bacteria such as *E. coli* and *Vibrio alginolyticus* in close proximity to a planar boundary (Berg Reference Berg2000; Lauga *et al.*
Reference Lauga, Diluzio, Whitesides and Stone2006). The effect is highly local, with the gap between the bacterium and the wall typically much smaller than the size of the bacterium itself. While point-singularity methods predict such behaviour just as a result of the
$C_{2}$
mode and indeed agree that it should be strongly localised close to the surface, with an inverse-fourth dependence on the separation (Lauga *et al.*
Reference Lauga, Diluzio, Whitesides and Stone2006; Lopez & Lauga Reference Lopez and Lauga2014; Davis & Crowdy Reference Davis and Crowdy2015; Papavassiliou & Alexander Reference Papavassiliou and Alexander2015), higher-order modes can be expected to play an important role at such small gap widths.

The induced rotation is calculated by performing the integral (4.1) using the slip velocity (4.7) and the stress corresponding to axisymmetric rotation,

where $c_{i},d_{i}$ are as given in (3.3) and (3.4). The factor of $V_{l}(\cos \unicode[STIX]{x1D6FD})$ may be written as a polynomial of order $l-1$ in $W$ ,

where the coefficients $w_{n}(\unicode[STIX]{x1D709})$ are determined using any of the various series representations of Legendre polynomials (Whittaker & Watson Reference Whittaker and Watson1996), so that (4.9) becomes

$(W/R)^{1/2}$ is the generating function for Legendre polynomials (Whittaker & Watson Reference Whittaker and Watson1996). Successive differentiations of the generating function give the identity

which reduces (4.9) to a pair of integrals over orthogonal associated Legendre polynomials. Performing these integrals and resumming the results we find that near a concave shell or wall the motion is

while for interaction with a tracer it is

where the torques are given by (3.5) and (3.6). To find $\tilde{\unicode[STIX]{x1D6FA}}_{1}$ , $\unicode[STIX]{x1D714}_{2}$ must be chosen so that $T_{2}=0$ . Conversely, choosing $\unicode[STIX]{x1D714}_{2}$ such that $T_{1}=0$ allows the tracer motion $\tilde{\unicode[STIX]{x1D6FA}}_{2}$ to be found. These expressions are exact, for any separation, any axisymmetric slip velocity and any of the geometries covered by bispherical coordinates. We will show in § 4.3 that this reproduces the rotational hydrodynamic interactions that have been determined previously using asymptotic methods, such as minimal reflections of point singularities. First, however, we examine the limit of small separation, where the swimmers approach contact and the hydrodynamic interactions are strongest.

The limit of vanishing separation is treated using dominated convergence as described in § 3.2. For a shell this gives

and for a tracer

It can be readily verified that these coincide for $r_{1}/r_{2}\rightarrow 0$ with the value

representing the rotation of a squirmer touching a no-slip wall.

To illustrate the near-field behaviour that can be found exactly using the reciprocal theorem, we give a specific example of a swimmer whose slip velocity is an azimuthal circulation within a polar-cap region of opening angle $\unicode[STIX]{x1D703}_{0}$ . Although crude, this provides a squirmer representation of a rotating flagellar bundle, and counter-rotating cell body. Explicitly, we take the slip velocity to be

as depicted schematically in figure 4(*a*). The slip velocity within the cap region is
$\unicode[STIX]{x1D6FA}_{c}$
, which is balanced by a counter-rotation of the body,
$\unicode[STIX]{x1D6FA}_{b}$
, chosen so that the coefficient
$C_{1}=0$
to remove any free rotation and focus on the effects of interactions. The squirming coefficients are given by

and the counter-rotation $\unicode[STIX]{x1D6FA}_{b}$ required to cancel out the free rotation is

When
$\unicode[STIX]{x1D703}_{0}=\unicode[STIX]{x03C0}/2$
we have that
$\unicode[STIX]{x1D6FA}_{b}=\unicode[STIX]{x1D6FA}_{c}$
, as expected on symmetry grounds. *E. coli* has a body counter-rotation measured to be of the order of one-tenth the rotation of its flagellar bundle, with large variation between specimens (Magariyama, Sugiyama & Kudo Reference Magariyama, Sugiyama and Kudo2001), and inversion of (4.20) gives an appropriate value of approximately
$0.28\unicode[STIX]{x03C0}$
for
$\unicode[STIX]{x1D703}_{0}$
, which we idealise as
$\unicode[STIX]{x03C0}/4$
.

The dependence of the swimmer’s rotation on its orientation at large distances is given by the slowest-decaying squirming mode,
$C_{2}$
, and hence by
$P_{2}(\cos \unicode[STIX]{x1D6FC})$
, which is head–tail symmetric. However in the near field there may be significant asymmetry in the orientation dependence which could persist for relatively large separations. Figure 4(*c*) shows how the orientation dependence changes for distances up to 100 times the swimmer radius, for a model *E. coli* interacting with a no-slip wall. A comparison between the interaction with a no-slip wall of a spherical-cap swimmer calculated using all modes up to
$C_{100}$
, and an equivalent squirming sphere with only the dominant far-field
$C_{2}$
mode, is shown in figure 4(*b*) and further illustrates the importance of including higher-order modes in calculating near-field interactions: at separations of the order of the swimmer’s size the effect of including the higher-order modes can be dramatic. In the case of
$\unicode[STIX]{x1D6FC}=\unicode[STIX]{x03C0}/4$
the swimmer’s rotation changes sense as it approaches the wall. When
$\cos \unicode[STIX]{x1D6FC}=3^{-1/2}$
the contribution of the
$C_{2}$
mode is identically zero since
$P_{2}(3^{-1/2})=0$
; however there is still motion driven by higher-order modes of non-negligible magnitude.

### 4.3 Asymptotics at large separation

The tendency of flagellated bacteria such as *E. coli* and *Vibrio alginolyticus* to follow distinctive circular trajectories near boundaries (DiLuzio *et al.*
Reference DiLuzio, Turner, Mayer, Garstecki, Weibel, Berg and Whitesides2005) has resulted in considerable theoretical work on the rotation of a swimmer normal to a planar surface (Lauga *et al.*
Reference Lauga, Diluzio, Whitesides and Stone2006; Di Leonardo *et al.*
Reference Di Leonardo, Dell’Arciprete, Angelani and Iebba2011; Lopez & Lauga Reference Lopez and Lauga2014; Papavassiliou & Alexander Reference Papavassiliou and Alexander2015). This phenomenon has been explained by noting that the rotation of the body and counter-rotation of the tail gives a flow field resembling a rotlet dipole and, when aligned parallel to a wall (as extensile swimmers generically do (Spagnolie & Lauga Reference Spagnolie and Lauga2012)), results in rotation by a mechanism analogous to the turning of a tank. The availability of the exact solutions, (4.13) and (4.14), for this component of motion allow extension to further geometries than those appearing in the literature.

The flow field generated by the $C_{l}$ contribution to the slip velocity has an asymptotic decay of $d^{-(l+2)}$ in an unbounded domain (Pak & Lauga Reference Pak and Lauga2014). Focusing on the $l=1$ contribution we recognise the sum in (4.13) and (4.14) as the torque, (3.6) and (3.5) respectively. Hence the contribution to the rotation from this squirming mode is $-C_{1}\cos \unicode[STIX]{x1D6FC}$ for any separation and in any configuration. This is precisely the same as the rotation found asymptotically (Pak & Lauga Reference Pak and Lauga2014; Davis & Crowdy Reference Davis and Crowdy2015), and demonstrates that this mode corresponds only to self-rotation and does not result in interaction. Therefore, the slowest-decaying contribution to the normal rotation of a squirmer due to interactions is from $C_{2}$ , which represents a rotlet dipole. Here we discuss the interaction of this squirming mode with a passive sphere as a leading-order behaviour which is generic for all swimmers.

Using dimensional analysis, Lopez & Lauga (Reference Lopez and Lauga2014) argued that a squirmer circling parallel to a wall has an asymptotic angular frequency decaying no slower than
$d^{-4}$
; this was confirmed by considering the flow induced by a rotlet dipole near a wall, and indeed has been found to be the leading-order behaviour of a swimmer with arbitrary azimuthal slip velocity near a wall, with
$\tilde{\unicode[STIX]{x1D6FA}}_{1}=C_{2}(r_{1}/2d)^{4}P_{2}(\cos \unicode[STIX]{x1D6FC})/5$
(Papavassiliou & Alexander Reference Papavassiliou and Alexander2015). Figure 3(*a*) shows that this behaviour agrees with our exact solution (4.13) up to a separation of approximately a squirmer diameter. An explicit form for the rotation near a wall is obtained from equation (4.13) by setting
$\unicode[STIX]{x1D709}_{2}=0$
,
$\tilde{\unicode[STIX]{x1D6FA}}_{2}=0$
and
$\unicode[STIX]{x1D709}_{1}=\log (d/r_{1}+\sqrt{d^{2}/r_{1}^{2}-1})$
, and it is found that the rotation induced by interaction with the wall goes as
$d^{-(l+2)}$
for
$C_{l}$
. The behaviour in a shell, shown in figure 3(*b*), also has this form since the separation is always smaller than the radius of curvature of the shell. Note that the rotation of the swimmer when it is precisely in the centre is zero by symmetry and changes sense as the swimmer crosses between hemispheres.

The asymptotic rotation of a swimmer in the presence of a tracer may be calculated using the leading-order forms $\unicode[STIX]{x1D709}_{1}\sim \log (r_{1}/d)$ and $\unicode[STIX]{x1D709}_{2}\sim \log (r_{2}/d)$ , giving a decay of

This may be understood in terms of multipole reflections (Kim & Karrila Reference Kim and Karrila2013). At large separation the swimmer’s motion is driven by the flow reflected in the tracer, which has the leading behaviour of a stresslet since the tracer must remain force free. Dimensional analysis suggests that the reflected flow at the swimmer should have a strength going as
$d^{-6}$
, and therefore a vorticity of
$d^{-7}$
, but for this case of an axisymmetric, azimuthal flow the leading reflected flow is identically zero, so the rotation is driven by a vorticity of
$d^{-9}$
. Figure 3(*c*) shows a cross-over to this behaviour when the separation exceeds the radius of the tracer. In the near field the passive sphere resembles as a wall and we see a
$d^{-4}$
dependence of the swimmer’s rotational speed. For the passive sphere (not shown) the cross-over is not seen, and the asymptotic interaction is
$d^{-4}$
, equal to the asymptotic vorticity generated by a rotlet dipole; thus, the dominant effect of the
$C_{2}$
squirming mode is the motion of the tracer, and by superposition two squirmers with this slip velocity would tend to move each other more than themselves.

### 4.4 Rotation near a free surface

An interesting application of the exact solution presented above is to find the rotation of a squirmer close to a free surface. It has been hypothesised (Lauga *et al.*
Reference Lauga, Diluzio, Whitesides and Stone2006) and subsequently observed experimentally (Di Leonardo *et al.*
Reference Di Leonardo, Dell’Arciprete, Angelani and Iebba2011) that the circular trajectories of *E. coli* near a free surface have the opposite sense to those near a no-slip wall. Using the known hydrodynamic solution for the rotation of a sphere beneath the interface between two fluid phases we find that both cases of rotation near a wall and a free surface may be described as image systems, using the two-sphere solution presented previously. This allows the swimming close to such boundaries to be found and compared without further calculation, and we find that the change of direction depending on the type of boundary is generic and explained by these image systems.

If a sphere rotates beneath the flat interface between the fluid that contains it, and another fluid of viscosity $\tilde{\unicode[STIX]{x1D707}}$ (O’Neill & Ranger Reference O’Neill and Ranger1979), the flow may be found explicitly by supposing an ansatz of the form (3.2) in each phase and matching flow and stress across the boundary. Then the torque on the sphere is

where $\unicode[STIX]{x1D6EC}=(\unicode[STIX]{x1D707}-\tilde{\unicode[STIX]{x1D707}})/(\unicode[STIX]{x1D707}+\tilde{\unicode[STIX]{x1D707}})$ . When $\unicode[STIX]{x1D6EC}=-1$ the empty phase is infinitely viscous and corresponds a no-slip wall; instead, when $\unicode[STIX]{x1D6EC}=+1$ the boundary is a free surface. Since the torque (3.5) corresponding to the two-sphere solution when $r_{1}=r_{2}=r$ is

it can be seen that the result for a free surface is recovered when $\unicode[STIX]{x1D714}_{2}=\unicode[STIX]{x1D714}_{1}$ (Brenner Reference Brenner1964), while the result for a no-slip wall is given by $\unicode[STIX]{x1D714}_{2}=-\unicode[STIX]{x1D714}_{1}$ . Hence a rotating sphere near a free surface has as its image system a corotating sphere which decreases the torque compared to the free-space value, while near a wall the image system is an antirotating sphere which increases the torque.

The rotation near a free surface may then be calculated exactly using the reciprocal theorem and compared to our expressions for squirming near a wall, (4.13). Although the activity of the squirmer generates tangential flows on the interface, since the stress in the conjugate problem is zero there is no contribution to the reciprocal theorem from an integral over the free surface and an expression for the rotation is obtained immediately from (4.14) by substituting $\unicode[STIX]{x1D709}_{2}=-\unicode[STIX]{x1D709}_{1}$ and $\unicode[STIX]{x1D714}_{2}=\unicode[STIX]{x1D714}_{1}$ , giving

In both cases of a wall and a free surface the leading far-field contribution from the $l$ -th squirming mode is equal and opposite, with a strength $\tilde{\unicode[STIX]{x1D6FA}}_{1}\propto d^{-(l+2)}$ . Hence the nature of the boundary determines the sense of rotation generically in the asymptotic limit. In the contact limit, which for a wall is given by (4.17) and for a free surface

we find that the contribution to the rotation from each squirming mode is smaller at a free surface than at a wall, and the ratio of these contributions has a faster-than-exponential decay with increasing
$l$
, indicating that higher-order effects due to microscopic details of a swimmer are less important at a free surface than at a wall. This can be seen in figure 4(*d*), which gives the orientation dependence of the same spherical-cap swimmer as considered before near a free surface; compared to the analogous trace for rotation near a wall, figure 4(*c*), the one for a free surface has the opposite sign and is smoother.

The reciprocal theorem relies on the swimmer problem and the conjugate Stokes drag being defined in the same region. Here we have assumed that the free surface does not deform in either solution, so that this region is the half-space with an embedded sphere in both cases; however, there is no reason not to expect deformation, particularly in the close proximity regime and furthermore, it cannot be assumed that this deformation of the surface will be the same for a swimmer and a dragged sphere. Nevertheless, if the deformation of the surface is sufficiently small it may be treated in a linear fashion by projecting onto the plane and expressing as a slip velocity in both solutions (much as Lighthill’s deforming squirmer has its activity projected onto the surface of a sphere for determination of the swimming speed (Lighthill Reference Lighthill1952)).

### 4.5 Translation

Phenomena such as the tendency of swimming microorganisms to aggregate at surfaces (Berke *et al.*
Reference Berke, Turner, Berg and Lauga2008) have been neatly explained by modelling a swimmer as a collection of point singularities. The sign of a swimmer’s stresslet, which characterises it as extensile or contractile (Ramaswamy Reference Ramaswamy2010), causes it to align parallel or normal to the wall, respectively (Spagnolie & Lauga Reference Spagnolie and Lauga2012), while the inclusion of a source dipole ensures self-propulsion (Drescher *et al.*
Reference Drescher, Goldstein, Michel, Polin and Tuval2010). By adopting the reciprocal theorem the behaviour due to any slip velocity may, in principle, be found (Papavassiliou & Alexander Reference Papavassiliou and Alexander2015); while here the conjugate solution used restricts us to axisymmetric motions, we have the freedom to generalise to curved surfaces. The actual calculation for the translational motion is analogous to the calculation for rotation shown in § 4.2, but is rather more involved and will not be shown explicitly; explicit expressions are given in appendix B. The general result for an arbitrary squirming mode has not yet been found and each contribution must be calculated separately. Instead, we will attempt to describe the behaviour using a few illustrative examples.

We consider the first few translational squirming modes,
$A_{1}$
,
$B_{1}$
,
$A_{2}$
and
$B_{2}$
. The first two of these set the self-propulsive speed in free space and asymptotically resemble source dipoles.
$A_{2}$
and
$B_{2}$
give the asymptotic stresslet of the swimmer (Ishikawa *et al.*
Reference Ishikawa, Simmonds and Pedley2006) and while they generate no motion in an unbounded domain they are of fundamental importance in the interactions of the swimmer with boundaries, since both the swimmer and the boundaries must remain force free and the lowest-order image singularity will be a stresslet. In the far field these point-singularity descriptions are sufficient to fully characterise the generic behaviour (Spagnolie & Lauga Reference Spagnolie and Lauga2012). For the special case of interaction with a wall an explicit asymptotic estimate (Davis & Crowdy Reference Davis and Crowdy2015; Papavassiliou & Alexander Reference Papavassiliou and Alexander2015) of the swimming speed is available as

thus, asymptotically, $A_{1}$ and $B_{1}$ contribute behaviour that differs only in a numerical factor, while behaviour due to $A_{2}$ and $B_{2}$ is distinguishable only by a sign change.

Figures 5 and 6 show, respectively, the normal translation speed of a swimmer induced by the
$A_{1}$
and
$B_{1}$
, and the
$A_{2}$
and
$B_{2}$
modes respectively, in interaction with a passive concave or convex sphere. It can be seen that the far-field equivalence of
$A_{1}$
and
$B_{1}$
, and
$A_{2}$
and
$B_{2}$
, also holds for the concave and convex geometries. These figures indicate that the cross-over to far-field behaviour that is well described by point-singularity models (Ishikawa *et al.*
Reference Ishikawa, Simmonds and Pedley2006; Spagnolie & Lauga Reference Spagnolie and Lauga2012; Davis & Crowdy Reference Davis and Crowdy2015; Papavassiliou & Alexander Reference Papavassiliou and Alexander2015) occurs at very small separations, of the order of a few swimmer diameters.

The
$A_{2}$
and
$B_{2}$
squirming modes generate an asymptotic flow field of
$d^{-2}$
, while the propulsive modes
$A_{1}$
and
$B_{1}$
give a flow field decaying as
$d^{-3}$
. Hence mixing is dominated by the swimmers’ dipoles, and the speed of a passive tracer has a dependence of
$d^{-2}$
, by Fáxen’s law, until the separation becomes small and higher-order effects become important. Since
$A_{2}$
and
$B_{2}$
do not drive self-propulsion, the motion of the swimmer resulting from these modes is due to reflected flow in the boundary of the passive sphere. At separations smaller than the tracer’s radius of curvature the leading order of the reflected flow is equal to that of the flow itself, and gives rise to a local
$d^{-2}$
behaviour. As the separation increases the finite size of the tracer becomes important. Higdon (Reference Higdon1979) gives the image system for a force dipole in a fixed, finite-sized sphere as the sum of a Stokeslet with leading-order strength proportional to
$d^{-2}$
and a dipole with strength
${\sim}d^{-3}$
. Thus, if the passive sphere were fixed the leading-order reflection would go as
$d^{-3}$
, but as it is free to move in such a way as to cancel any force acting on it, we see
$d^{-5}$
. This dependence may also be calculated using a second-order multipole expansion, in which case the leading-order motion of the swimmer is driven by the reflected stresslet inside the tracer (Kim & Karrila Reference Kim and Karrila2013). The cross-over between the two types of behaviour is shown in figure 6(*c*).

The availability of exact solutions for the motion due to these squirming modes means it is possible to calculate explicit trajectories in time, albeit only for motion along the common diameter of two spheres. Hence we consider a head-on collision of the swimmer with a wall, which allows comparison to analogous trajectories calculated by integrating the approximate results (4.26). This is shown in figure 7, and it can be seen that the trajectories differ very little between radial and tangential modes (red and black lines respectively), and the corresponding asymptotic approximation (grey dashed line). These trajectories become distinct only at very small separation, again of the order of a swimmer diameter.

In contrast, the near-field behaviour due to radial and tangential slip is rather different. It can be seen from figures 5 and 6 that the tangential modes $B_{1}$ and $B_{2}$ result in a swimming speed which goes to zero as contact with the surface is approached; this results in collision taking a longer time than predicted by a point-singularity. The radial modes $A_{1}$ and $A_{2}$ result in acceleration to a finite speed as contact is approached; this results from the incompatibility of the boundary conditions of no-slip and radial flow on two touching surfaces. The consequences of this can be seen in figure 7, where the radial modes $A_{1}$ and $A_{2}$ result in collisions in finite time while the tangential modes $B_{1}$ and $B_{2}$ cause deceleration close to contact. Although we have been unable to explicitly integrate the expressions of the swimming speed to determine whether physical contact occurs within finite time or not, we note that the exact solution for a swimming disc with tangential slip collides with a wall in infinite time, as may be verified using the equations of motion calculated by Crowdy (Reference Crowdy2011).

The hydrodynamic force in the conjugate problem is divergent in the near field. This divergence has a leading part going as $F_{i}\sim \unicode[STIX]{x1D709}_{i}^{-2}$ , but the representation of the force as an infinite sum contains a harmonically divergent subleading term. By approximating the sum as an integral Cox & Brenner (Reference Cox and Brenner1967) found that this subleading divergence may be expressed as ${\sim}\log \unicode[STIX]{x1D709}_{i}$ as $\unicode[STIX]{x1D709}_{i}\rightarrow 0$ . When calculating swimmer motions using the reciprocal theorem the force appears as a denominator, with the numerator given by the integral of the slip velocity against the conjugate stress tensor. This numerator also diverges, although no faster than the force; specifically, the divergence is the same as for the force for the integrals involving $A_{1}$ and $A_{2}$ , and one power of $\unicode[STIX]{x1D709}$ slower for those with $B_{1}$ and $B_{2}$ . This is enough to explain the behaviour shown in figures 5–7 and a detailed analysis of the subleading terms in the vein of Cox & Brenner (Reference Cox and Brenner1967) is unnecessary.

An interesting consequence of the near-field distinction between radial and tangential slip is the behaviour of a swimmer inside a small shell, with a radius smaller than the threshold for cross-over to asymptotic behaviour. When a
$B_{1}$
swimmer is inside such a shell the swimming speed is attenuated by the presence of the boundaries, and by symmetry attains a maximum value,
$U_{c}$
, in the centre of the shell. Conversely, a swimmer with
$A_{1}$
activity has an increased speed due to the interactions, as a result of the divergent interaction of the radial modes near boundaries. A comparison is shown for a variety of shell sizes in figure 5(*a*). The speed at the centre of the shell,
$U_{c}$
, may be calculated analytically and depends on the relative sizes as

Thus as the shell becomes large,
$U_{c}$
approaches the free swimming speed in proportion to the volume of fluid displaced by the swimmer, see figure 5(*b*). The same occurs in two dimensions where an exact result is available for the swimming of an active disc inside a circular boundary (Papavassiliou & Alexander Reference Papavassiliou and Alexander2015) and where the maximum speed of a self-propulsive swimmer approaches its free-space value in proportion to the excluded area.

## 5 Discussion

We have found exact expressions for the axisymmetric translation and rotation of a spherical squirmer close to a convex, planar or concave no-slip boundary, as well as the axisymmetric rotation beneath a free surface, by making use of the Lorentz reciprocal theorem and the known Stokes drag solutions in these geometries. This covers the hydrodynamics at all separations, including at contact in the case of rotation, and for arbitrary squirming motion. The near-field regime, where separations are comparable to the swimmer size or smaller, is the regime of greatest relevance to many experimental settings and our exact solution provides rigorous, generic insight. In particular, while the radial and meridional squirming modes show the same asymptotic behaviour, in the near field, at separations smaller than a couple of swimmer diameters, they are markedly different, with the former giving a divergent interaction strength and the latter a hard repulsion. Azimuthal squirming results in the circling behaviour near boundaries seen in flagellated bacteria and our results describe this situation in some detail. The experimentally reported reversal of orbit direction at a free surface is found to be a generic effect. At large separations, the exact solution reproduces results found previously from asymptotic calculations using point-singularity approximations of swimmers and also generalises these to interactions of squirmers with spherical boundaries and tracer particles.

Our solution is founded upon the reciprocal theorem for swimmer problems (Stone & Samuel Reference Stone and Samuel1996) and appears to be the first application of this method to deduce exact solutions that are not currently available by any other method. Given the widespread significance of hydrodynamic interactions, with confining surfaces and with other organisms, to swimmer motion, there are obvious merits to developing applications of this technique in other settings. For instance, we have only been able to provide a partial solution to the interaction of two swimmers, as the non-axisymmetric components of the motion have not been determined. This is because the solution is founded upon the reciprocal theorem and requires the corresponding Stokes drag problem to be solved. For the non-axisymmetric Stokes drag of two spheres, there is, at present, no exact closed-form solution, although there is a scheme in terms of a set of difference equations that could be solved numerically to any desired degree of accuracy. Such an approach would allow the full hydrodynamic interaction of an arbitrary pair of squirmers to be computed, although not in closed form. Furthermore the large range of validity of the approximate far-field solutions here compared with our exact results indicates that asymptotic estimates are valuable and there is merit to pursuing an approximate approach to find the non-axisymmetric behaviour. This may be done, for instance, by constructing an approximate stress tensor using the solution for a Stokeslet outside a sphere (Higdon Reference Higdon1979).

## Acknowledgements

We are grateful to T. Machon, M. Polin and especially G. Rowlands for fruitful discussions. This work was partially supported by the UK EPSRC through grant no. A.MACX.0002.

## Appendix A. The Stimson–Jeffery drag force

The real coefficients ${\mathcal{A}}_{l},{\mathcal{B}}_{l},{\mathcal{C}}_{l},{\mathcal{D}}_{l}$ appearing in (3.11) are found by inverting the solution schemes (3.12)–(3.13).

It is easiest to consider the concave ( $\unicode[STIX]{x1D709}_{1}>\unicode[STIX]{x1D709}_{2}>0$ ) and convex ( $\unicode[STIX]{x1D709}_{1}>0>\unicode[STIX]{x1D709}_{2}$ ) geometries separately. In the former we find

*a*-

*d*) $$\begin{eqnarray}\displaystyle {\mathcal{A}}_{l} & = & \displaystyle 2L(V_{1}-V_{2})\text{e}^{-(2l+1)(\unicode[STIX]{x1D709}_{1}+\unicode[STIX]{x1D709}_{2})}\left[{\displaystyle \frac{(2l+1)}{(2l+3)}}(\text{e}^{(2l-1)\unicode[STIX]{x1D709}_{1}}-\text{e}^{(2l-1)\unicode[STIX]{x1D709}_{2}})+\text{e}^{(2l+1)\unicode[STIX]{x1D709}_{2}}-\text{e}^{(2l+1)\unicode[STIX]{x1D709}_{1}}\right]\nonumber\\ \displaystyle {\mathcal{B}}_{l} & = & \displaystyle 2L(V_{1}-V_{2})\text{e}^{-(2l+1)(\unicode[STIX]{x1D709}_{1}+\unicode[STIX]{x1D709}_{2})}\left[{\displaystyle \frac{(2l+1)}{(2l-1)}}\left(\text{e}^{(2l+3)\unicode[STIX]{x1D709}_{1}}-\text{e}^{(2l+3)\unicode[STIX]{x1D709}_{2}}\right)+\text{e}^{(2l+1)\unicode[STIX]{x1D709}_{2}}-\text{e}^{(2l+1)\unicode[STIX]{x1D709}_{1}}\right]\nonumber\\ \displaystyle {\mathcal{C}}_{l} & = & \displaystyle L\bigg[(2l+3)(V_{1}+V_{2})+(2l+1)(\text{e}^{-2\unicode[STIX]{x1D709}_{2}}-\text{e}^{-2\unicode[STIX]{x1D709}_{1}})(V_{1}-V_{2})\nonumber\\ \displaystyle & & \displaystyle +\,{\displaystyle \frac{4(\text{e}^{-(2l+1)(\unicode[STIX]{x1D709}_{1}-\unicode[STIX]{x1D709}_{2})}V_{1}+\text{e}^{(2l+1)(\unicode[STIX]{x1D709}_{1}-\unicode[STIX]{x1D709}_{2})}V_{2})-(2l+1)^{2}(2(\unicode[STIX]{x1D709}_{1}-\unicode[STIX]{x1D709}_{2})V_{1}+\text{e}^{-2(\unicode[STIX]{x1D709}_{1}-\unicode[STIX]{x1D709}_{2})}V_{2})}{(2l-1)}}\bigg]\nonumber\\ \displaystyle {\mathcal{D}}_{l} & = & \displaystyle -L\bigg[(2l-1)(V_{1}+V_{2})+(2l+1)(\text{e}^{2\unicode[STIX]{x1D709}_{2}}-\text{e}^{2\unicode[STIX]{x1D709}_{1}})(V_{1}-V_{2})\nonumber\\ \displaystyle & & \displaystyle +\,{\displaystyle \frac{4(\text{e}^{-(2l+1)(\unicode[STIX]{x1D709}_{1}-\unicode[STIX]{x1D709}_{2})}V_{1}+\text{e}^{(2l+1)(\unicode[STIX]{x1D709}_{1}-\unicode[STIX]{x1D709}_{2})}V_{2})-(2l+1)^{2}(-2(\unicode[STIX]{x1D709}_{1}-\unicode[STIX]{x1D709}_{2})V_{1}+\text{e}^{2(\unicode[STIX]{x1D709}_{1}-\unicode[STIX]{x1D709}_{2})}V_{2})}{(2l+3)}}\bigg],\nonumber\\ \displaystyle & & \displaystyle\end{eqnarray}$$

where

while in the convex geometry we have

*a*-

*d*) $$\begin{eqnarray}\displaystyle {\mathcal{A}}_{l} & = & \displaystyle L\left[2V_{1}\text{e}^{-(2l+1)(\unicode[STIX]{x1D709}_{1}+\unicode[STIX]{x1D709}_{2})}\left({\displaystyle \frac{(2l+1)}{(2l+3)}}(\text{e}^{(2l-1)\unicode[STIX]{x1D709}_{1}}-\text{e}^{(2l-1)\unicode[STIX]{x1D709}_{2}})-\text{e}^{(2l+1)\unicode[STIX]{x1D709}_{1}}+\text{e}^{(2l+1)\unicode[STIX]{x1D709}_{2}}\right)\right.\nonumber\\ \displaystyle & & \displaystyle \left.+\,V_{2}\left({\displaystyle \frac{(2l+1)^{2}\text{e}^{-2(\unicode[STIX]{x1D709}_{1}-\unicode[STIX]{x1D709}_{2})}-4\text{e}^{-(2l+1)(\unicode[STIX]{x1D709}_{1}-\unicode[STIX]{x1D709}_{2})}}{(2l+3)}}-(2l-1)-(2l+1)(\text{e}^{-2\unicode[STIX]{x1D709}_{1}}-\text{e}^{-2\unicode[STIX]{x1D709}_{2}})\right)\right]\nonumber\\ \displaystyle {\mathcal{B}}_{l} & = & \displaystyle L\left[2V_{1}\text{e}^{-(2l+1)(\unicode[STIX]{x1D709}_{1}+\unicode[STIX]{x1D709}_{2})}\left({\displaystyle \frac{(2l+1)}{(2l-1)}}(\text{e}^{(2l+3)\unicode[STIX]{x1D709}_{1}}-\text{e}^{(2l+3)\unicode[STIX]{x1D709}_{2}})-\text{e}^{(2l+1)\unicode[STIX]{x1D709}_{1}}+\text{e}^{(2l+1)\unicode[STIX]{x1D709}_{2}}\right)\right.\nonumber\\ \displaystyle & & \displaystyle \left.+\,V_{2}\left({\displaystyle \frac{4\text{e}^{-(2l+1)(\unicode[STIX]{x1D709}_{1}-\unicode[STIX]{x1D709}_{2})}-(2l+1)^{2}\text{e}^{2(\unicode[STIX]{x1D709}_{1}-\unicode[STIX]{x1D709}_{2})}}{(2l-1)}}+(2l+3)+(2l+1)(\text{e}^{2\unicode[STIX]{x1D709}_{1}}-\text{e}^{2\unicode[STIX]{x1D709}_{2}})\right)\right]\nonumber\\ \displaystyle {\mathcal{C}}_{l} & = & \displaystyle L\left[2V_{2}\left({\displaystyle \frac{(2l+1)}{(2l-1)}}\text{e}^{-2(\unicode[STIX]{x1D709}_{1}+\unicode[STIX]{x1D709}_{2})}(\text{e}^{(2l+3)\unicode[STIX]{x1D709}_{1}}-\text{e}^{(2l+3)\unicode[STIX]{x1D709}_{2}})-\text{e}^{(2l+1)\unicode[STIX]{x1D709}_{1}}+\text{e}^{(2l+1)\unicode[STIX]{x1D709}_{2}}\right)\right.\nonumber\\ \displaystyle & & \displaystyle \left.+\,V_{1}\left({\displaystyle \frac{4\text{e}^{-(2l+1)(\unicode[STIX]{x1D709}_{1}-\unicode[STIX]{x1D709}_{2})}-(2l+1)^{2}\text{e}^{2(\unicode[STIX]{x1D709}_{1}-\unicode[STIX]{x1D709}_{2})}}{(2l-1)}}+(2l+3)-(2l+1)(\text{e}^{-2\unicode[STIX]{x1D709}_{1}}-\text{e}^{-2\unicode[STIX]{x1D709}_{2}})\right)\right]\nonumber\\ \displaystyle {\mathcal{D}}_{l} & = & \displaystyle L\left[2V_{2}\left({\displaystyle \frac{(2l+1)}{(2l+3)}}\text{e}^{2(\unicode[STIX]{x1D709}_{1}+\unicode[STIX]{x1D709}_{2})}(\text{e}^{(2l-1)\unicode[STIX]{x1D709}_{1}}-\text{e}^{(2l-1)\unicode[STIX]{x1D709}_{2}})-\text{e}^{(2l+1)\unicode[STIX]{x1D709}_{1}}+\text{e}^{(2l+1)\unicode[STIX]{x1D709}_{2}}\right)\right.\nonumber\\ \displaystyle & & \displaystyle \left.+\,V_{1}\left({\displaystyle \frac{(2l+1)^{2}\text{e}^{-2(\unicode[STIX]{x1D709}_{1}-\unicode[STIX]{x1D709}_{2})}-4\text{e}^{-(2l+1)(\unicode[STIX]{x1D709}_{1}-\unicode[STIX]{x1D709}_{2})}}{(2l+3)}}-(2l-1)+(2l+1)(\text{e}^{2\unicode[STIX]{x1D709}_{1}}-\text{e}^{2\unicode[STIX]{x1D709}_{2}})\right)\right].\nonumber\\ \displaystyle & & \displaystyle\end{eqnarray}$$

Using these expressions the explicit form of the force on two spheres may be found using (3.16).

## Appendix B. Expressions for the translational motion

In what follows we give explicit expressions for the coaxial translation of two spheres, as driven by the squirming modes $A_{1},A_{2},B_{1},B_{2}$ on the surface $\unicode[STIX]{x1D709}=\unicode[STIX]{x1D709}_{1}$ (for squirming on $\unicode[STIX]{x1D709}=\unicode[STIX]{x1D709}_{2}$ , $\unicode[STIX]{x1D709}_{1}$ and $\unicode[STIX]{x1D709}_{2}$ should be interchanged throughout). The reciprocal theorem, (4.1), results in an expression for $\tilde{U} _{1}F_{1}+\tilde{U} _{2}F_{2}$ , where the forces $F_{1}$ and $F_{2}$ are given by (3.16).

To isolate $\tilde{U} _{1}$ and $\tilde{U} _{2}$ separately appropriate choices of the coefficients should be made such that one of the two forces is zero, or so that the forces are equal and opposite, in which case the relative, rather than absolute, motions are found.

### B.1 Radial modes

The reciprocal theorem integral corresponding to $A_{1}$ gives the expression

while that for $A_{2}$ gives

### B.2 Tangential modes

The reciprocal theorem integral corresponding to $B_{1}$ gives the expression

while that for $B_{2}$ gives