text:
# 1 Skyrme solitons

## 1 Skyrme solitons

The fundamental field of Skyrme’s theory  is a map

$`U:^3SU(2).`$ (1.1)
We denote points in $`^3`$ by $`x`$ with coordinates $`x_i`$, $`i=1,2,3`$ and Euclidean length $`r=|x|=\sqrt{x_1^2+x_2^2+x_3^2}`$ . Sometimes we write $`\widehat{x}`$ for the unit vector $`x/r`$. It is often useful to parametrise $`U`$ in terms of the Pauli matrices $`\tau _1,\tau _2`$ and $`\tau _3`$ and the triplet of pion fields $`\pi _1,\pi _2`$ and $`\pi _2`$ as

$`U(x)=\sigma (x)+i\pi _a(x)\tau _a,`$ (1.2)
where summation over the repeated index $`a`$ is implied and the field $`\sigma `$ is determined by the constraint $`\sigma ^2+\pi _1^2+\pi _2^2+\pi _3^2=1`$. In this introductory section we do not specify the class of functions to which $`U`$ belongs. It is assumed to be sufficiently smooth for all the following operations to make sense.

The Skyrme energy functional is best written in terms of the Lie-algebra valued currents

$`L_i=U^{}_iU`$ (1.3)
or

$`R_i=U_iU^{},`$ (1.4)
where $`_i=/x_i`$. It is

$`E[U]={\displaystyle d^3x\left(\frac{1}{2}\mathrm{tr}(L_iL_i)+\frac{1}{16}\mathrm{tr}([L_j,L_i][L_j,L_i])\right)}.`$ (1.5)
The Euler-Lagrange equation for stationary points of this functional is conveniently expressed in terms of the modified currents

$`\stackrel{~}{L}_i=L_i{\displaystyle \frac{1}{4}}[L_j,[L_j,L_i]]`$ (1.6)
and

$`\stackrel{~}{R}_i=R_i{\displaystyle \frac{1}{4}}[R_j,[R_j,R_i]],`$ (1.7)
where we again use the convention that repeated indices are summed over. It reads

$`_i\stackrel{~}{L}_i=0`$ (1.8)
or, equivalently,

$`_i\stackrel{~}{R}_i=0.`$ (1.9)
Here we are interested in finite-energy solutions of the Euler-Lagrange equation. It is shown in  that the finite energy requirement implies that the map $`U`$ tends to a constant value at infinity in a weak sense. We choose that constant to be the identity element $`1SU(2)`$ and demand

$`\underset{r\mathrm{}}{lim}U(x)=1.`$ (1.10)
The boundary condition (1.10) means that the domain of $`U`$ is effectively compactified to a three-sphere. Since the target space is also a three-sphere, maps satisfying (1.10) have an associated integer degree. The first rigorous proof that for a finite-energy Skyrme configuration in a very general class of functions the degree

$`\text{deg}[U]={\displaystyle \frac{1}{24\pi ^2}}{\displaystyle d^3xϵ_{ijk}\mathrm{tr}\left(L_iL_jL_k\right)}`$ (1.11)
is an integer was given in . This result means that the space

$`𝒞=\{U:^3SU(2)|E[U]<\mathrm{}\}`$ (1.12)
of finite-energy configurations is a disjoint union of sectors

$`𝒞_k=\{U:^3SU(2)|E[U]<\mathrm{}\text{and \hspace{0.17em} deg}[U]=k\}`$ (1.13)
labelled by the integers $`k`$.

The symmetry group of Skyrme’s theory will play an important role in our discussion. The energy functional (1.5), the boundary condition (1.10) and the degree (1.11) are invariant under the action of the Euclidean group of translations and rotations in $`^3`$ and under rotations of the pion fields

$`\pi _aG_{ab}\pi _b,GSO(3),`$ (1.14)
which we call iso-rotations. Reflections in space $`S:xx`$ and in iso-space $`\pi _a\pi _a`$ both leave the energy invariant but each changes the sign of the degree. The pull-back of Skyrme configurations $`U`$ via $`S`$ provides a map

$`\stackrel{~}{S}:𝒞_k𝒞_k,UUS`$ (1.15)
which preserves the energy.

It was shown in  that the energy in each topological sector is bounded below by a multiple of the degree. It follows from the results of  that the bound cannot be attained for the standard version of the Skyrme model described here, so that we have the strict inequality

$`E[U]>12\pi ^2|k|.`$ (1.16)
The bound ensures that the infima

$`I_k=\text{inf}\{E[U]|U𝒞_k\}`$ (1.17)
are well-defined. The question of whether the infima are attained was first addressed by Esteban in the paper . Amongst other things Esteban proved that, for a suitable class of functions,

$`I_kI_l+I_{kl}`$ (1.18)
for all $`k,l`$. She also showed that infima are attained provided one assumes the strict inequality

$`I_k<I_l+I_{kl}`$ (1.19)
for all $`k\{0,\pm 1\}`$ and $`l\{0,k\}`$ in the range $`|l|+|kl|<\sqrt{2}|k|`$. In  it was shown that the result still holds if one widens the class of allowed functions but the inequality (1.19) remains a necessary assumption in the proof. The strict inequality is also of interest in physics. As we shall see it is related to the question of attractive forces in the Skyrme model.

For low values of $`k`$ the existence and nature of minima of the Skyrme energy functional is understood in more detail. For fields of degree one, the highly symmetric hedgehog ansatz

$`U_H(x)=\mathrm{exp}(if(r)\widehat{x}_a\tau _a),`$ (1.20)
introduced already by Skyrme, leads to an ordinary differential equation for the profile function $`f`$. With the boundary conditions $`f(0)=\pi `$ and $`f(\mathrm{})=0`$ the resulting Skyrme configuration has degree 1 and is called the Skyrmion. It was shown in  that the Skyrmion minimises the Skyrme energy functional amongst all degree one configurations of the hedgehog form. In  the existence of the minimum of the Skyrme energy functional in $`𝒞_1`$ was proved, but it has not been established rigorously that minimising configurations have the symmetry of the hedgehog field (1.20). Note that if $`U_1`$ is a minimal energy configuration in $`𝒞_1`$ then the reflected configuration $`\stackrel{~}{S}\left(U_1\right)`$ has the same energy and minimises the energy in $`𝒞_1`$.

In the following we use the term Skyrme solitons for minimal energy solutions of the Skyrme equation with non-zero degree $`k`$. There is overwhelming numerical evidence that the minimum in $`𝒞_2`$ is attained by a configuration of toroidal symmetry . For higher degree, too, much is known numerically about the minima of Skyrme’s energy functional in the sector $`𝒞_k`$. Numerical searches, assisted by analytical ansätze and investigations of the possible symmetries of Skyrme solitons, suggest the existence of Skyrme solitons in all sectors $`𝒞_k`$ up to $`k=22`$. The energies are sufficiently accurately computed that it appears that the inequality (1.19) is satisfied for all $`k`$ in the range $`2k20`$ and all $`l`$ in the range $`0<l<k`$, see .

The existence of attractive forces between two Skyrmions was already shown by Skyrme, using the product ansatz. In this paper we investigate the existence of attractive forces between general Skyrme solitons. It is perhaps worth stressing that our arguments in fact apply to any finite-energy solution of the Skyrme equation, not just minimal ones. An earlier attempt at proving the existence of attractive forces between Skyrme solitons was made in the unpublished paper . Our approach is partly inspired by ideas in  but also fills important gaps left there. Our main tool is an asymptotic expansion of Skyrme solitons. In Section 2 we show that Skyrme fields have an asymptotic expansion in powers of $`1/r`$ and that for non-trivial Skyrme solitons that expansion always has a non-trivial leading multipole. In Section 3 we study the interaction energy of two Skyrme solitons and show that it is dominated by a certain linear interaction energy provided the orders of the leading multipoles of the Skyrme solitons do not differ by more than 2. In Section 4 we derive an expression for the linear interaction energy of two multipoles, and show that it can always be made negative by suitable relative rotations in space and iso-space. At the end of this paper we briefly comment on the relationship between our results and Esteban’s work, and on the implications for the existence of Skyrme solitons of arbitrary degree.

## 2 The asymptotics of Skyrme solitons

The aim of this section is to show that if $`U`$ and the currents $`L_j`$ are a little better than continuous, then $`U`$ is smooth in $`^3`$, and has a non-trivial asymptotic expansion in powers of $`1/r`$ and $`\mathrm{log}r`$ as $`r\mathrm{}`$. By non-trivial we mean here that if $`U`$ is non-constant, then it cannot happen that all terms in the asymptotic expansion vanish. Put another way, it is not possible for $`U`$ to approach $`1`$ at infinity faster than every power of $`1/r`$ unless $`U`$ is identically equal to $`1`$ in $`^3`$.

To make a precise statement, say that the (possibly matrix-valued) function $`f`$ is in $`C_b^{0,\alpha }(^3)`$, with $`0<\alpha <1`$, if $`f`$ satisfies

$$\underset{x^3}{sup}|f|+\underset{x,x^{}^3,xx^{}}{sup}\frac{(1+r+r^{})^\alpha |f(x)f(x^{})|}{|rr^{}|^\alpha +(1+r+r^{})^\alpha d(\omega ,\omega ^{})^\alpha }<\mathrm{}.$$
(2.1)
Here $`x=r\omega ,x^{}=r^{}\omega ^{}`$, where $`\omega `$ is regarded as an angular variable (unit vector) living on the unit $`2`$-sphere, and $`d(\omega ,\omega ^{})`$ is the $`2`$-sphere distance between $`\omega `$ and $`\omega ^{}`$. If we consider the space $`C_b^{0,\alpha }(K)`$, where $`K`$ is a bounded subset of $`^3`$, then this is precisely the same as the usual space of (bounded) Hölder-continuous functions, with Hölder exponent $`\alpha `$. However, $`C_b^{0,\alpha }(^3)`$ is slightly different from the usual space of bounded Hölder-continuous functions on $`^3`$.

The main result of this section can be summarized as follows:

###### Theorem 2.1

Suppose that for some $`\delta >0`$,

$$(1+r)^\delta (U1)C_b^{0,\alpha },(1+r)^\delta L_jC_b^{0,\alpha }.$$
(2.2)
Suppose further that $`U`$ satisfies the Skyrme differential equations (1.8) in the sense of distributions.

Then $`U`$ is smooth in $`^3`$ and has a complete asymptotic expansion in powers of $`1/r`$ and $`\mathrm{log}r`$, for large $`r`$. If $`U`$ is non-constant, then this expansion has a leading term which is harmonic, hence a multipole.

Note that the hypotheses (2.2) force $`U`$ to approach $`1`$ and the $`L_j`$ to approach $`0`$ like $`r^\delta `$ as $`r\mathrm{}`$. As a technical remark, we point out that the assumption of Hölder-continuous currents implies that $`U`$ will have a Hölder continuous first derivative. The second derivatives $`_j_kU`$ are then defined only in the sense of distributions, but in the Skyrme equation $`_j_kU`$ enters linearly, and is multiplied by continuous functions (smooth functions of the currents) see Section 2.1 below. In particular, the left-hand side of the Skyrme equation $`_j\stackrel{~}{L}_j=0`$ makes sense as a distribution if (2.2) holds.

The assumptions (2.2) of Theorem (2.1) do not follow immediately from the variational analysis used by Esteban in . Her methods only give that the derivatives $`_jU`$ are locally square-integrable (and that the components of $`U`$ are locally bounded). On physical grounds, one expects minimizers of the Skyrme energy functional to satisfy (2.2), but it would be desirable to bridge the gap between the analysis given here and what was proved rigorously in  and . This issue will not be pursued further here.

The proof of Theorem (2.1) proceeds in four steps, each of which takes up one of the following subsections. In the first subsection we re-write the Skyrme equation in order to make explicit the form of the non-linearities. The equation is quasilinear, in the sense that the derivatives of highest order (2) enter linearly. The Skyrme equation can therefore be regarded as a second-order linear elliptic PDE, with Hölder-continuous coefficients. Then standard regularity theorems (Schauder estimates) yield that $`U`$ is smooth.

In Section 2.2 we study the behaviour of the Skyrme equation at spatial infinity, by introducing coordinates $`(s,\omega )`$, where $`s=r^1`$, so that the 2-sphere at infinity becomes a genuine boundary at $`s=0`$. After a rescaling, the Laplacian of $`^3`$ is replaced by

$$\mathrm{\Delta }_b=(s_s)^2s_s+\mathrm{\Delta }_\omega ,$$
(2.3)
where $`\mathrm{\Delta }_\omega `$ is the Laplacian of the unit $`2`$-sphere. The analysis is now guided by the corresponding analysis of a system of ordinary differential equations with regular singular point at $`s=0`$. In particular, one does not expect the solutions to be smooth near $`s=0`$, but one does expect a non-trivial expansion in powers of $`s`$ (and possibly $`\mathrm{log}s`$). The $`b`$-calculus of  enables us to make this kind of argument precise. Thus in Section 2.2 we show that if $`U=1+u`$, then $`u`$ is conormal at $`s=0`$, which is to say that $`u`$ and all derivatives of the form $`(s_s)^m_\omega ^nu`$ are continuous as $`s0.`$ (This condition is strictly weaker than $`u`$ being smooth near $`s=0`$.) In Section 2.3, we show that it is not possible for $`U`$ to approach $`1`$ faster than $`r^\mu `$ for every $`\mu >0`$ unless $`U=1`$ in $`^3`$. Finally in Section 2.4 we combine this fact with another application of the $`b`$-calculus to show the existence of a non-trivial asymptotic expansion in powers of $`r^1`$ and $`\mathrm{log}r`$ (or equivalently in powers of $`s`$ and $`\mathrm{log}s`$).

### 2.1 Rewriting the Skyrme Equation

To begin with we work near a fixed point of $`^3`$, which we may as well take to be the origin $`0`$. By replacing $`U(x)`$ by $`U(0)^1U(x)`$ we can assume that $`U(0)=1`$. Write

$$U(x)=1+u(x),$$
(2.4)
so that $`u(0)=0`$. Because $`U`$ is unitary, the $`2\times 2`$ complex matrix $`u`$ will satisfy the algebraic constraints

$$u+u^{}+uu^{}=0,\mathrm{tr}(u)+detu=0.$$
(2.5)
In particular $`u`$ is neither exactly skew-hermitian nor trace-free. Then

$$L_i=_iu+u^{}_iu$$
(2.6)
and

$$_iL_i=(1+u^{})\mathrm{\Delta }u+_iu^{}_iu.$$
(2.7)
The cubic term in the currents can be written

$$\frac{1}{2}L_j[L_j,L_i]=\frac{1}{2}(_ju+u^{}_ju)[_ju+u^{}_ju,_iu+u^{}_iu].$$
(2.8)
Taking the divergence, and using the notation

$$v_{ij}=_iu^{}_ju,$$
(2.9)
we obtain

$$\frac{1}{2}L_j[L_j,L_i]=(1+u^{})(T+F),$$
(2.10)
where

$$T=T(u,u,^2u)=\frac{1+u}{2}(L_j[L_j,(1+u^{})\mathrm{\Delta }u]+L_j[(1+u^{})_iu_ju,L_i])$$
(2.11)
and

$$F=F(u,u)=\frac{1+u}{2}(v_{ij}[L_i,L_j]+L_j[L_j,v_{ii}]+L_j[(1+u^{})v_{ij},L_i]).$$
(2.12)
The nonlinear terms have been divided here so that $`T`$ is polynomial in the first derivatives of $`u`$ and linear in its second derivatives, while $`F`$ only contains $`u`$ and its first derivatives.

With this notation, the full Skyrme equation can be written

$$P(u,u,^2u)=Q(u,u)+F(u,u)$$
(2.13)
where

$$P(u,u,^2u)=\mathrm{\Delta }u+T(u,u,^2u)\text{ and }Q(u,u)=(1+u)v_{ii}.$$
(2.14)
Notice that $`Q`$ is (approximately) quadratic in $`u`$, $`T`$ is of degree three and $`F`$ is of degree four. Note also that $`T`$ is linear in $`_i_ju`$ (and quadratic in $`_iu`$).

Because of the assumed Hölder continuity of the currents, the coefficients of the differential operator

$$fP(u,u,^2f)$$
(2.15)
are Hölder continuous, and this operator is linear and elliptic in a small neighbourhood $`K`$ of $`0`$. Moreover, the RHS $`Q+F`$ of (2.13) is also in $`C_b^{0,\alpha }(K)`$, so by elliptic regularity, it follows that $`uC_b^{2,\alpha }(K)`$. Then the currents are in $`C_b^{1,\alpha }(K)`$ and the process continues to show that $`u`$ is in $`C_b^{k,\alpha }(K)`$ for every $`k`$. Thus $`u`$ is smooth near $`0`$. Since the point $`0`$ was arbitrary, the argument shows that $`u`$ is smooth in $`^3`$.

### 2.2 Boundary regularity

In order to analyze the asymptotic behaviour of the field $`U`$, we shall make a transformation of the problem which involves passing from $`^3`$ to a ‘compactification’ $`\overline{^3}`$ in which the sphere at infinity becomes a genuine boundary. This is easily achieved by introducing the coordinate $`s=1/r`$ along with angular coordinates $`\theta `$ and $`\phi `$ in $`^3`$. Then in the space $`[0,\mathrm{})\times S^2`$ with coordinates $`(s,\theta ,\phi )`$, the boundary $`s=0`$ corresponds to $`r=\mathrm{}`$ and $`\theta `$ and $`\phi `$ give coordinates on the boundary, which is the 2-sphere ‘at infinity’. $`\overline{^3}`$ is defined to be the union of $`^3`$ with the $`2`$-sphere at infinity attached in this way.

It can be cumbersome to work with explicit coordinates on the $`2`$-sphere, so we again use $`\omega `$ for points on $`S^2`$ and represent any point other than the origin of $`\overline{^3}`$ in the form $`(s,\omega )`$.

Next, introduce rescaled derivatives

$$D_i=r_i=\frac{1}{s}_i.$$
(2.16)
These vector fields have the property that they are linear combinations, with coefficients that are smooth, down to $`s=0`$, of the basic vector fields $`s_s`$, $`_\theta `$ and $`_\phi `$. The Euclidean Laplacian takes the form

$$\mathrm{\Delta }=s^2\mathrm{\Delta }_b$$
(2.17)
with $`\mathrm{\Delta }_b`$ defined in (2.3).

Now we write $`U=1+u`$ for large $`r`$ (that is, for small positive $`s`$) and make the rescalings of (2.13) suggested by (2.16) and (2.17). The result is a ‘$`b`$’ version of the Skyrme equation,

$$P_b(u,Du,D^2u)=Q_b(u,Du)+s^2F_b(u,Du),$$
(2.18)
where

$`P_b(u,Du,D^2u)`$ $`=`$ $`\mathrm{\Delta }_bu+s^2T_b(u,Du,D^2u),`$ (2.19)
$`Q_b(u,Du)`$ $`=`$ $`(1+u)D_iu^{}D_iu,`$ (2.20)
these terms being got by replacing $`_i`$ by the rescaled derivative $`D_i`$ wherever they occur.

The reason for reformulating the equation in this way is that there is a well-established theory, called the $`b`$-calculus, which can be used to analyze equations of this kind . The $`b`$-calculus is concerned with $`b`$-differential operators. In the present situation, a $`b`$-differential operator is just a differential operator of the form

$$P=\underset{a+b+cm}{}C_{abc}(s,\omega )(s_s)^a_\theta ^b_\phi ^c,$$
(2.21)
where the coefficients $`C_{abc}`$ are smooth up to the boundary $`s=0`$. From (2.3) it is clear that the rescaled Laplacian $`\mathrm{\Delta }_b`$ is an example of such an operator. The set of all such operators will be denoted by $`\text{Diff}_b`$, those of order $`k`$ by $`\text{Diff}_b^k`$.

The first aspect of this theory that is needed is the counterpart of the elliptic regularity for Hölder spaces that we used in the previous subsection. For this we need $`b`$-Hölder spaces, already introduced in (2.1). With the new variable $`s`$, we have $`f(s,\omega )C_b^{0,\alpha }`$ if

$$\underset{x^3}{sup}|f|+\underset{(s,\omega )(s^{},\omega ^{})}{sup}\frac{(s+s^{})^\alpha |f(s,\omega )f(s^{},\omega ^{})|}{|ss^{}|^\alpha +(s+s^{})^\alpha d(\omega ,\omega ^{})^\alpha }<\mathrm{}.$$
(2.22)
Now put

$$C_b^{k,\alpha }=\{f:LfC_b^{0,\alpha }\text{ for all }b\text{-differential operators }L\text{Diff}_b^k\}$$

and

$$𝒜=\{u:LuC_b^{0,\alpha }\text{ for all }L\text{Diff}_b\}.$$

In order to force functions to decay as $`s0`$, we introduce weighted versions of these spaces,

$$s^\delta C_b^{k,\alpha }=\{u=s^\delta v:vC_b^{k,\alpha }\},s^\delta 𝒜=\{u=s^\delta v:v𝒜\}.$$

The following is a very special case of elliptic regularity for $`b`$-differential operators .

###### Theorem 2.2

Consider the differential operator $`P=\mathrm{\Delta }_b+s^2E`$, where $`E`$ is a second-order differential operator with coefficients smooth up to the boundary $`s=0`$. Suppose that $`Pu=f`$ near $`s=0`$, with $`u`$ and $`f`$ in $`s^\delta C_b^{0,\alpha }`$. Then if $`\delta `$ is not an integer, it follows that $`us^\delta C_b^{2,\alpha }`$.

We want to apply this to the Skyrme equation, written in the form (2.19). After the last section, we know that the coefficients of the perturbing term $`E`$ are smooth for $`s>0`$, but all we know at $`s=0`$ is the original assumption that the currents are in $`s^\delta C_b^{0,\alpha }`$. The elliptic regularity statement still holds in this case (though this statement does not seem to be available in the literature). Applying this result, for $`\delta `$ a small positive number, we obtain first that $`us^\delta C_b^{2,\alpha }`$, which gives that the coefficients are in $`s^\delta C_b^{1,\alpha }`$, so $`us^\delta C_b^{3,\alpha }`$ and so on. Iterating, we find $`us^\delta 𝒜`$.

This is a major step forward: in particular, the angular variation of $`u`$ is now very well controlled. However $`u`$ may still be very far from being smooth up to the boundary or having an asymptotic expansion there. Indeed, horrors like $`s^\delta \mathrm{sin}\mathrm{log}s`$ lie in $`s^\delta 𝒜`$.

### 2.3 Power-law decay of $`u`$

In this section and the next we establish that a topologically non-trivial solution of the Skyrme equation must have a non-trivial asymptotic expansion in powers of $`1/r`$. We show first that it is not possible for a topologically non-trivial solution to approach $`1`$ faster than every power of $`1/r`$. This result is then fed into an iterative analysis of the equation in the next section. These two sections can, however, be read in either order. The main result of this section is as follows.

###### Theorem 2.3

If the topological charge of the Skyrme field $`U`$ is non-zero and if it satisfies the hypotheses of Theorem 2.1, then there exists some $`\mu ^+`$ such that $`r^\mu u(r,\omega )`$ does not tend to zero as $`r\mathrm{}`$.

To set this result in a more general context, recall that a partial differential equation $`Lu=0`$ is said to have the unique continuation property at a point $`0`$, say, if the following is true:

$$\text{If all derivatives of }u\text{ vanish at }0\text{, then }u=0\text{ in some neighbourhood of }0\text{.}$$
(2.23)
The methods used to establish that a second-order PDE has the unique continuation property establish analogous statements with somewhat weaker hypotheses. For example, a simplified version of Theorem 17.2.6 of  is as follows.

###### Theorem 2.4

Let $`a_{jk}(x)`$ be smooth and positive-definite in a neighbourhood $`X`$ of $`\mathrm{\hspace{0.17em}\hspace{0.17em}0}`$, and suppose that $`a_{jk}(0)=\delta _{jk}`$. Suppose that for $`xX`$, the smooth function $`u`$ satisfies

$$|a_{jk}(x)_j_ku|A\left(|u(x)|+|u(x)|\right)$$
(2.24)
for some constant $`A`$, and

$$|x|^\mu |u(x)|0\text{ as}|x|0\text{for }\text{every}\text{ }\mu ^+.$$
(2.25)
Then $`u=0`$ identically in a neighbourhood of $`0`$.

Using this result we can prove the following:

###### Theorem 2.5

Let $`U`$ satisfy the hypotheses of Theorem (2.1). Then if the topological charge of $`U`$ is non-zero, $`U`$ cannot be constant in any open subset of $`^3`$.

Proof Define

$$W=\{x^3:U(y)=U(x)\text{ for all }y\text{ in some neighbourhood of }x\}.$$
(2.26)
Then $`W`$ is open by definition. By Theorem 2.4, $`W`$ is also closed. To see this, suppose that $`x_nW`$, $`x_nx_0^3`$. Then all derivatives of $`U`$ are zero at $`x_0`$, so (2.25) holds (with $`0`$ replaced by $`x_0`$). The Skyrme equation, in the form (2.13), implies a differential inequality of the form (2.24) in a neighbourhood of $`x_0`$. It follows<sup>1</sup><sup>1</sup>1Unique continuation theorems do not extend wholesale to systems. However, in our case, the leading term is the Laplacian and the other second-order terms $`C(u,u,^2u)`$ are non-scalar but small near $`x_0`$. The proof of Theorem 2.4 goes through in this case. that $`U`$ is identically constant in a neighbourhood of $`x_0`$, so that $`x_0W`$. Since $`^3`$ is connected, $`W=\mathrm{}`$ or $`W=^3`$.

We will use this theorem to give an indirect proof of Theorem 2.3. Suppose $`r^\mu u(x)0`$ as $`r\mathrm{}`$ for every $`\mu `$. Adapating Theorem 2.4 we shall show that then $`u=0`$ for all sufficiently large $`r`$. Thus $`U`$ is constant in an open set, hence by the previous result constant everywhere.

So consider the $`b`$-differential operator

$$P_b=\mathrm{\Delta }_b+s^2T_b$$
(2.27)
on $`\overline{^3}`$, where the coefficients of $`T_b`$ are in $`C_b^{0,\alpha }`$. By rescaling Theorem 17.2.6 in , we obtain

###### Theorem 2.6

Suppose that

$$|P_bu(s,\omega )|As^\delta (|u(s,\omega )|+|Du(s,\omega )|)\text{ for all }0<s<s_0,\omega S^2$$
(2.28)
and that

$$s^\mu |u(s,\omega )|0\text{ as }s0\text{ for all }\mu .$$
(2.29)
Then $`u(s,\omega )=0`$ for $`0s<s_1`$, where $`s_1`$ is some small positive number.

The ‘$`b`$’ version (2.18) of the Skyrme equation implies a differential inequality of the form (2.28), just as before. It follows that if $`u`$ decays faster than any power of $`r`$, then $`u=0`$ for all sufficiently large $`r`$. By the remarks before the statement of Theorem 2.6, the proof of Theorem 2.3 is now complete.

### 2.4 Refined regularity, asymptotic expansions

The main result of this section can be stated as follows:

###### Theorem 2.7

Let $`U=1+u`$ satisfy the hypotheses of Theorem 2.1. Then there is some integer $`M1`$ and an asymptotic expansion

$$u\underset{j=M}{\overset{2M}{}}Y_j(\omega )r^{(j+1)}+\underset{j=2M+1}{\overset{\mathrm{}}{}}r^{(j+1)}w_j(\omega ,\mathrm{log}r)\text{ for large }r.$$
(2.30)
Here the $`Y_j`$ are Lie-algebra valued spherical harmonics,

$$\mathrm{\Delta }_\omega Y_j=j(j+1)Y_j,$$
(2.31)
and $`Y_M0`$, so that the piece $`_{j=M}^{2M}Y_j(\omega )r^{(j+1)}`$ is a non-zero harmonic function. The functions $`w_j`$ are smooth in $`\omega `$ and polynomial in $`\mathrm{log}r`$.

It will follow from the proof that the terms $`w_{2M+1}`$ to $`w_{3M+1}`$ are of at most degree $`1`$ in $`\mathrm{log}r`$, the terms $`w_{3M+2}`$ to $`w_{4M+2}`$ are of at most degree $`2`$ in $`\mathrm{log}r`$ and so on. We remark also that the asymptotic expansion can safely be differentiated term by term to give asymptotic expansions of all derivatives of $`u`$.

In order to motivate the proof, consider the equation

$$\mathrm{\Delta }_bu=f$$
(2.32)
where $`u`$ and $`f`$ are defined for small $`s`$. If $`f`$ has the form $`f=s^\lambda g(\omega )`$, then a solution can be found as follows. Expand $`g`$ as a sum of spherical harmonics,

$$g=\underset{j=0}{\overset{\mathrm{}}{}}g_j,\text{ where }\mathrm{\Delta }_\omega g_j=j(j+1)g_j$$

and seek a solution $`u=u_j(s)g_j`$. Then $`u_j`$ must satisfy

$$[(s_s)^2(s_s)j(j+1)]u_j(s)=s^\lambda .$$

This is solved by

$$u_j(s)=\frac{s^\lambda }{\lambda (\lambda 1)j(j+1)}$$

provided there is no resonance, that is to say

$$\lambda j,j+1.$$

In the resonant case, $`u_j(s)`$ has the form $`s^\lambda (A+B\mathrm{log}s)`$. The general solution is obtained by combining this with an arbitrary solution of the homogeneous equation $`\mathrm{\Delta }_bv=0`$. If we require $`u0`$ as $`s0`$, then $`v`$ must itself go to zero and hence will be a sum of multipoles, $`v=_{j=0}^{\mathrm{}}Y_j(\omega )s^{j+1}`$, where $`Y_j`$ satisfies (2.31).

The $`b`$-calculus extends results of this kind to functions in $`s^\delta 𝒜`$ (which, as we have seen, can be far from having expansions in powers of $`s`$). In order to summarize the needed results, write

$$f=𝒪(s^\delta )\text{ instead of }fs^\delta 𝒜.$$
(2.33)
This conforms to the use of the ‘$`𝒪`$’-notation in the rest of the paper, but has the additional property

$$\text{If }f=𝒪(s^\delta )\text{ then }Lf=𝒪(s^\delta )\text{ for all }L\text{Diff}_b.$$
(2.34)
###### Lemma 2.8

Suppose that $`u`$ and $`f`$ defined near $`s=0`$ satisfy (2.32). If $`u=𝒪(s^\delta )`$ and $`f=𝒪(s^{n+\delta })`$, where $`\delta >0`$ and $`n`$ is a positive integer, then

$$u=\underset{j=0}{\overset{n1}{}}Y_j(\omega )s^{j+1}+𝒪(s^{n+\delta })$$
(2.35)
where $`Y_j`$ satisfies (2.31). In particular the sum on the RHS is harmonic

$$\mathrm{\Delta }_b\underset{j=0}{\overset{n1}{}}Y_j(\omega )s^{j+1}=0.$$
(2.36)
###### Lemma 2.9

Suppose that $`u`$ and $`f`$ defined near $`s=0`$ satisfy (2.32). If $`u=𝒪(s^\delta )`$ and

$$f=\underset{j=0}{\overset{n1}{}}f_j(\omega ,\mathrm{log}s)s^{j+1}+𝒪(s^{n+\delta }),$$
(2.37)
where $`\delta >0`$, $`n`$ is a positive integer and $`f_j`$ is polynomial of degree $`m_j`$ in $`\mathrm{log}s`$, then

$$u=\underset{j=0}{\overset{n1}{}}w_j(\omega ,\mathrm{log}s)s^{j+1}+𝒪(s^{n+\delta })$$
(2.38)
where $`w_j`$ is a polynomial of degree $`m_j+1`$ in $`\mathrm{log}s`$.

These results will be applied to the Skyrme equation, now rewritten as

$$\mathrm{\Delta }_bu=Z_b(u):=Q_b(u)s^2T_b(u)+s^2F_b(u).$$
(2.39)
###### Lemma 2.10

Suppose that $`U=1+u`$ satisfies the Skyrme equation and $`u=𝒪(s^\delta )`$. Then $`u=𝒪(s^2)`$.

Proof: On the RHS of (2.39), $`Q_b`$ is quadratic in $`Du`$ and the other terms in $`Z_b`$ are of even higher degree. Hence $`Z_b=𝒪(s^{2\delta })`$. Applying Lemma 2.36 we obtain that $`u=𝒪(s^{2\delta })+𝒪(s)`$. If $`2\delta <1`$, we can iterate this argument to obtain eventually $`u𝒪(s)`$. In , it is shown that a solution of the Skyrme equation cannot have leading term $`1/r=s`$ in its asymptotic expansion. It follows that $`u=𝒪(s^{1+\delta })`$ (for a possibly smaller $`\delta >0`$). Hence $`Z_b(u)=𝒪(s^{2+2\delta })`$ and so, applying Lemma 2.36 again, $`u=Y_1s^2+𝒪(s^{2+\delta })`$.

Combining Theorem 2.3 with Lemma 2.36, we see that there exists an integer $`M1`$ with the property that

$$u=Y_M(\omega )s^{M+1}+𝒪(s^{M+1+\delta }),$$
(2.40)
where $`Y_M`$ is a non-vanishing spherical harmonic. We can now complete the proof of Theorem 2.7 in the following iterative fashion. From the structure of $`Z_b(u)`$ it follows from (2.40) that

$$Z_b(u)=f_{2M+1}s^{2M+2}+𝒪(s^{2M+2+\delta }).$$
(2.41)
Applying Lemma 2.9,

$$u=\underset{j=M}{\overset{2M}{}}Y_j(\omega )s^{j+1}+w_{2M+1}(\omega ,\mathrm{log}s)s^{2M+2}+𝒪(s^{2M+2+\delta }),$$
(2.42)
where $`w_{2M+1}(\omega ,\mathrm{log}s)`$ is of degree at most $`1`$ in $`\mathrm{log}s`$. Now this expression for $`u`$ is substituted into $`Z_b`$, to give

$$Z_b=\underset{2M+1}{\overset{3M+1}{}}f_j(\omega )s^{j+1}+f_{3M+2}(\omega ,\mathrm{log}s)s^{3M+3}++𝒪(s^{3M+3+\delta })$$
(2.43)
where $`f_{3M+2}`$ is of degree at most $`1`$ in $`\mathrm{log}s`$. Now apply Lemma 2.9 to get

$`u={\displaystyle \underset{j=M}{\overset{2M}{}}}Y_j(\omega )s^{j+1}+{\displaystyle \underset{j=2M+1}{\overset{3M+1}{}}}w_j(\omega ,\mathrm{log}s)s^{j+1}+w_{3M+2}(\omega ,\mathrm{log}s)s^{3M+3}+𝒪(s^{3M+3+\delta }),`$ (2.44)
where the functions $`w_{2M+1},\mathrm{},w_{3M+1}`$ are of degree at most 1 in $`\mathrm{log}s`$ and $`w_{3M+2}`$ is of degree at most 2 in $`\mathrm{log}s`$. Carrying on in this way we obtain a complete asymptotic expansion.

To complete the proof of Theorem 2.7 note finally that the first $`M+1`$ terms in the expansion are actually the pion field. Indeed, if the expansion (2.30) is substituted into the algebraic constraints (2.5) we see that the harmonic piece $`_M^{2M}Y_js^{j+1}`$ is skew-adjoint and trace-free—the quadratic corrections enter at order $`s^{2M+2}`$.

The upshot of this section is that every Skyrme soliton has a leading Lie-algebra valued multipole field (called a $`2^M`$-pole)

$`u_M(x)=i\tau _a{\displaystyle \underset{m=M}{\overset{M}{}}}{\displaystyle \frac{4\pi }{2M+1}}Q_{Mm}^a{\displaystyle \frac{Y_{Mm}(\theta ,\phi )}{r^{M+1}}},`$ (2.45)
where $`Y_{Mm}`$ are the usual spherical harmonics on $`S^2`$, see appendix A. The leading multipole moment $`Q_{Mm}^a`$ is independent of the location of the Skyrme soliton, and is acted on naturally by rotations and iso-rotations. It is a key ingredient in the calculations of the following sections. As already mentioned one can show that Skyrme solitons never have asymptotic monopole fields . The leading multipole field of the $`B=1`$ hedgehog (1.20) is a triplet of dipoles, and dipoles are known to occur as leading multipoles in a number of Skyrme solitons. The highest leading multipole known from numerical work is an octupole, which occurs in a $`B=7`$ configuration with icosahedral symmetry .

## 3 The interaction energy of two Skyrme solitons

Suppose we have Skyrme solitons $`U^{(1)}`$ and $`U^{(2)}`$ of degrees $`k`$ and $`l`$ which minimise the energy in the sectors $`𝒞_k`$ and $`𝒞_l`$. Since the total energies of $`U^{(1)}`$ and $`U^{(2)}`$ are finite there must be balls $`B_1`$ and $`B_2`$ in $`^3`$ so that most of the energy of $`U^{(1)}`$ and $`U^{(2)}`$ is concentrated in, respectively, $`B_1`$ and $`B_2`$. Outside the balls $`B_1`$ and $`B_2`$ the asymptotic analysis of the previous section applies. Suppose that the leading multipole of $`U^{(1)}`$ is a $`2^M`$-pole and the leading multipole of $`U^{(2)}`$ is a $`2^N`$-pole. Denoting the radii of $`B_1`$ and $`B_2`$ by $`D_1`$ and $`D_2`$, and with the abbreviation (2.45) for a generic Lie-algebra valued multipole field we have

$`U^{(1)}(x)1+u_M(x)\text{for}xB_1`$ (3.1)
and

$`U^{(2)}(x)1+v_N(x)\text{for}xB_2.`$ (3.2)
Using the translational invariance of the Skyrme energy functional we can assume without loss of generality that $`B_1`$ is centred at $`X_+=(0,0,R/2)`$ and that $`B_2`$ is centred at $`X_{}=(0,0,R/2)`$, where $`R`$ is so large that $`B_1`$ and $`B_2`$ do not overlap, i.e. $`R>D_1+D_2`$. The parameter $`R`$ will be interpreted as the separation of the Skyrme solitons. There is clearly an ambiguity in the definition of such a separation parameter, but this does not affect our calculation of leading terms in the limit where $`R`$ becomes large. Then we define the following product configuration

$`U_R(x)=U^{(1)}(x)U^{(2)}(x).`$ (3.3)
This configuration has degree $`k+l`$ and we shall see shortly that its energy is finite, so that $`U_R𝒞_{k+l}`$.

Our goal is to study the energy of the product configuration $`U_R`$ perturbatively in the limit of large $`R`$ and to compute the leading terms in powers of $`1/R`$. A similar calculation for moving and spinning Skyrmions was performed in , where some further details are given. Let $`L_i^{(1)}`$ and $`\stackrel{~}{L}_i^{(1)}`$ be the currents (1.3) and (1.6) constructed out of $`U^{(1)}`$, and $`R_i^{(2)}`$ and $`\stackrel{~}{R}_i^{(2)}`$ be the currents (1.4) and (1.7) constructed out of $`U^{(2)}`$. Then one computes

$`E[U_R]`$ $`=`$ $`E[U^{(1)}]+E[U^{(2)}]+W_2+W_4.`$ (3.4)
The energies $`E[U^{(1)}]`$ and $`E[U^{(2)}]`$ are simply the energies of the Skyrme solitons $`U^{(1)}`$ and $`U^{(2)}`$ and therefore independent of $`R`$. The interaction terms $`W_2`$ and $`W_4`$ are given by integrals over $`^3`$

$`W_2={\displaystyle d^3xw_2}\text{and}W_4={\displaystyle d^3xw_4},`$ (3.5)
with integrands

$`w_2=\mathrm{tr}\left(L_i^{(1)}\stackrel{~}{R}_i^{(2)}+\stackrel{~}{L}_i^{(1)}R_i^{(2)}L_i^{(1)}R_i^{(2)}\right)`$ (3.6)
and

$`w_4={\displaystyle \frac{1}{8}}\mathrm{tr}([L_i^{(1)},R_j^{(2)}][L_i^{(1)},R_j^{(2)}]+[L_i^{(1)},R_j^{(2)}][R_i^{(2)},L_j^{(1)}])+[L_i^{(1)},L_j^{(1)}][R_i^{(2)},R_j^{(2)}]).`$ (3.7)
We shall see shortly that the term $`W_2`$ contains the leading contribution to the interaction energy. However, for the presentation of our method of computation it is more convenient to begin with the quartic term $`W_4`$. We split the integration region $`^3`$ into the balls $`B_1`$ and $`B_2`$ and the complement $`C=^3(B_1B_2)`$. To illustrate our method, consider the integral

$`I={\displaystyle \frac{1}{8}}{\displaystyle _{B_1}}d^3x\mathrm{tr}\left([L_i^{(1)},R_j^{(2)}][L_i^{(1)},R_j^{(2)}]\right){\displaystyle \frac{1}{4}}{\displaystyle _{B_1}}d^3x\mathrm{tr}\left(L_i^{(1)}L_i^{(1)}\right)\mathrm{tr}\left(R_j^{(2)}R_j^{(2)}\right).`$ (3.8)
The currents $`L_i^{(1)}`$ are smooth functions and hence bounded on the compact domain $`B_1`$. Therefore

$`|I|K{\displaystyle _{B_1}}d^3x\mathrm{tr}\left(R_j^{(2)}R_j^{(2)}\right).`$ (3.9)
for some positive constant $`K`$. Since $`B_1`$ is far away from the centre of soliton $`U^{(2)}`$, the leading contribution to the integral (3.9) is obtained by replacing $`R_j^{(2)}`$ by the leading multipole component $`i_jv_N`$. Using

$`|_jv_N|(x){\displaystyle \frac{K^{}}{|xX_{}|^{N+2}}}`$ (3.10)
for a further positive constant $`K^{}`$ we conclude that

$`I=𝒪\left({\displaystyle \frac{1}{R^{2N+4}}}\right).`$ (3.11)
A similar calculation for the other terms in $`W_4`$ shows that

$`{\displaystyle _{B_1}}d^3xw_4=𝒪\left({\displaystyle \frac{1}{R^{2N+4}}}\right).`$ (3.12)
Considering the contribution from $`B_2`$ we find by the same argument

$`{\displaystyle _{B_2}}d^3xw_4=𝒪\left({\displaystyle \frac{1}{R^{2M+4}}}\right).`$ (3.13)
Finally the remaining integral over $`C`$ can be estimated as follows. Define

$`F(R)={\displaystyle _C}d^3x{\displaystyle \frac{1}{|xX_+|^{2M+4}}}{\displaystyle \frac{1}{|xX_{}|^{2N+4}}},`$ (3.14)
noting that the integral converges for all values of $`R>D_1+D_2`$. Then there is a positive constant $`K^{\prime \prime }`$ such that

$`{\displaystyle _C}d^3xw_4<K^{\prime \prime }F(R).`$ (3.15)
The large $`R`$ behaviour of $`F`$ can be estimated by a scaling analysis. Fix $`R_0>D_1+D_2`$ and consider $`R>R_0`$. Changing integration variables $`x(R_0/R)x`$ one computes

$`F(R)=\left({\displaystyle \frac{R_0}{R}}\right)^{2M+2N+5}\left[F(R_0)+{\displaystyle _{S_1(R)S_2(R)}}d^3x{\displaystyle \frac{1}{|xX_+|^{2M+4}}}{\displaystyle \frac{1}{|xX_{}|^{2N+4}}}\right],`$ (3.16)
where $`S_1(R)`$ and $`S_2(R)`$ are the thick shells

$`S_1(R)`$ $`=`$ $`\{x^3|(R_0/R)D_1|(x_1,x_2,x_3(R_0/2)|<D_1\}`$
$`S_2(R)`$ $`=`$ $`\{x^3|(R_0/R)D_2|(x_1,x_2,x_3+(R_0/2)|<D_2\}`$ (3.17)
which converge to punctured balls

$`B_1^0`$ $`=`$ $`\{x^3\{(0,0,R_0/2)\}||(x_1,x_2,x_3(R_0/2)|<D_1\}`$
$`B_2^0`$ $`=`$ $`\{x^3\{(0,0,R_0/2)\}||(x_1,x_2,x_3+(R_0/2)|<D_2\}`$ (3.18)
in the limit $`R\mathrm{}`$. In that limit the integral over $`S_1(R)`$ diverges like $`R^{2N+1}`$ and that over $`S_2(R)`$ like $`R^{2M+1}`$. Combining this with with the factor $`R^{(2M+2N+5)}`$ we deduce that $`F(R)`$ decays for large $`R`$ like $`R^{(2N+4)}`$ and $`R^{(2M+4)}`$, just like the contributions (3.12) and (3.13). Combining all the terms, we conclude that the leading terms in $`W_4`$ decay for large $`R`$ according to $`R^{(2N+4)}`$ and $`R^{(2M+4)}`$.

In order to evaluate $`W_2`$ we first divide the region of integration into the half-spaces

$`H^+=\{x^3|x_3>0\}\text{and }H^{}=\{x^3|x_3<0\}.`$ (3.19)
In $`H^+`$ we replace $`U^{(2)}`$ by the leading term $`1+v_N`$ and in $`H^{}`$ we replace $`U^{(1)}`$ by the leading contribution $`1+u_M`$. The result is

$`W_2`$ $``$ $`{\displaystyle _{H^+}}d^3x{\displaystyle \frac{1}{4}}\mathrm{tr}\left(L_i^{(1)}[_jv_N,[_jv_N,_iv_N]]\right)\mathrm{tr}\left(_iv_N\stackrel{~}{L}_i^{(1)}\right)`$ (3.20)
$``$ $`{\displaystyle _H^{}}d^3x{\displaystyle \frac{1}{4}}\mathrm{tr}\left(R_i^{(2)}[_ju_M,[_ju_M,_iu_M]]\right)+\mathrm{tr}\left(_iu_M\stackrel{~}{R}_i^{(2)}\right).`$
Now integrating by parts and using the Euler-Lagrange equations $`_i\stackrel{~}{L}_i^{(1)}=_i\stackrel{~}{R}_i^{(2)}=0`$ for the individual Skyrme solitons we convert two of the terms into an area integral

$`{\displaystyle _{H^+}}d^3x\mathrm{tr}\left(_iv_N\stackrel{~}{L}_i^{(1)}\right)+{\displaystyle _H^{}}d^3x`$ $`\mathrm{tr}\left(_iu_M\stackrel{~}{R}_i^{(2)}\right)`$ (3.21)
$`=`$ $`{\displaystyle _{x_3=0}}𝑑x_1𝑑x_2\mathrm{tr}\left(v_N\stackrel{~}{L}_3^{(1)}+u_M\stackrel{~}{R}_3^{(2)}\right).`$
Since the $`x_1x_2`$-plane is far away from both Skyrme solitons the leading contribution to this area integral can be expressed entirely in terms of the asymptotic fields:

$`\mathrm{\Delta }E=2{\displaystyle \underset{a=1}{\overset{3}{}}}{\displaystyle _{x_3=0}}𝑑x_1𝑑x_2(u_M^a_3v_N^av_N^a_3u_M^a).`$ (3.22)
A simple scaling analysis shows that $`\mathrm{\Delta }E`$ falls off like $`R^{(N+M+1)}`$ for large $`R`$. We skip the details here because we shall show how to evaluate $`\mathrm{\Delta }E`$ exactly in the next section. The remaining terms in (3.20) can be estimated with the techniques used in estimating $`W_2`$. The result is

$`W_2=\mathrm{\Delta }E+𝒪\left({\displaystyle \frac{1}{R^{3N+6}}}\right)+𝒪\left({\displaystyle \frac{1}{R^{3M+6}}}\right)`$ (3.23)
Combining all terms in (3.4) we conclude that

$`E[U_R]=E[U^{(1)}]+E[U^{(2)}]+\mathrm{\Delta }E+𝒪\left({\displaystyle \frac{1}{R^{2N+4}}}\right)+𝒪\left({\displaystyle \frac{1}{R^{2M+4}}}\right).`$ (3.24)
Note that $`\mathrm{\Delta }E`$ is the leading contribution to the interaction energy if $`|NM|2`$ i.e. if the orders of the leading multipoles of the two Skyrme solitons differ by at most two. We will comment on the validity of this assumption at the end of this paper.

## 4 Harmonic functions and their interaction energy

In order to compute the interaction energy $`\mathrm{\Delta }E`$ we need to derive some general results about harmonic functions. We define the regions

$`H_\delta ^{}=\{x^3|x_3<\delta \}\text{and }H_\delta ^+=\{x^3|x_3>\delta \},`$ (4.1)
where the positive parameter $`\delta `$ is introduced for technical reasons. Then we introduce the spaces

$`^{}=\{f:H_\delta ^{}|\mathrm{\Delta }f=0,\underset{r\mathrm{}}{lim}f(x)=0\}`$ (4.2)
and

$`^+=\{g:H_\delta ^+|\mathrm{\Delta }g=0,\underset{r\mathrm{}}{lim}g(x)=0\}.`$ (4.3)
Elements of $`^{}`$ tend to zero at the boundary “at infinity” of $`H_\delta ^{}`$, elements of $`^+`$ tend to zero at the boundary “at infinity” of $`H_\delta ^+`$. No additional restriction is placed on behaviour at the boundaries $`x_3=\pm \delta `$.

For the calculations in this section it is convenient to split $`^3`$ into $`^2\times `$ and denote vectors in $`^2`$ by bold letters, e.g. $`𝒙=(x_1,x_2)`$. We then write $`x=(𝒙,x_3)`$. The most general element of $`^{}`$ can be written as

$`f(x)={\displaystyle \frac{d^2𝒌}{(2\pi )^22k}p(𝒌)\mathrm{exp}(i𝒌𝒙+kx_3)},`$ (4.4)
where $`k=\sqrt{𝒌^2}`$ and the volume element $`d^2𝒌/((2\pi )^22k)`$ arises from the combination of $`d^3k`$ with the delta-function $`\delta (𝒌^2k^2)`$ which ensures that $`f`$ satisfies the Laplace equation. Since $`f`$ is real the Fourier transform $`p`$ satisfies

$`\overline{p}(𝒌)=p(𝒌).`$ (4.5)
Similarly, the most general element of $`^+`$ can be written as

$`g(x)={\displaystyle \frac{d^2𝒍}{(2\pi )^22l}q(𝒍)\mathrm{exp}(i𝒍𝒙lx_3)},`$ (4.6)
with $`l=\sqrt{𝒍^2}`$ and

$`\overline{q}(𝒌)=q(𝒌).`$ (4.7)
There is a natural pairing between elements of the $`^{}`$ and those of $`^+`$

$`f,g={\displaystyle _{x_3=0}}𝑑x_1𝑑x_2(g_3ff_3g)\text{for}f^{},g^+,`$ (4.8)
Using the expansion (4.4) and (4.6) we find, in terms of the Fourier modes,

$`f,g={\displaystyle \frac{d^2𝒌}{(2\pi )^22k}p(𝒌)q(𝒌)}.`$ (4.9)
This pairing is of interest to us since the interaction energy (3.22) is proportional to the sum over the pairings $`u_M^a,v_N^a`$. Therefore, we also refer to the expressions (4.8) and (4.9) as the interaction energy of the harmonic functions $`f`$ and $`g`$.

It is clear that the pairing (4.8) may vanish for some pairs of harmonic functions $`f^{}`$, $`g^+`$. This happens for example if the support of the Fourier transform $`p`$ is complementary to that of the Fourier transform $`q`$. However, we shall now show that the interaction energy of multipoles can always be made non-zero by rotating one of the functions.

### 4.1 Multipole fields

For $`f^{}`$ we have an alternative expansion in spherical harmonics $`Y_{Mm}`$

$`f(x)={\displaystyle \underset{M0}{}}{\displaystyle \underset{m=M}{\overset{M}{}}}{\displaystyle \frac{4\pi }{(2M+1)}}Q_{Mm}{\displaystyle \frac{Y_{Mm}(\theta _+,\phi _+)}{r_+^{M+1}}},`$ (4.10)
where $`r_+=|xX_+|`$ and $`(\theta _+,\phi _+)`$ are spherical coordinates centred at $`X_+=(0,0,R/2)`$. In this section we only need to assume $`R>0`$, but in our applications we will be interested in the large $`R`$ limit. The coefficients $`Q_{Mm}`$ are called the multipole moments of the function $`f`$. Assume that $`f`$ has non-vanishing multipole moments and suppose $`M`$ is the smallest integer such that $`Q_{Mm}0`$ for some $`m=M,\mathrm{},M`$. The function

$`f_M(x)={\displaystyle \frac{4\pi }{(2M+1)}}{\displaystyle \underset{m=M}{\overset{M}{}}}Q_{Mm}{\displaystyle \frac{Y_{Mm}(\theta _+,\phi _+)}{r_+^{M+1}}}.`$ (4.11)
is a $`2^M`$-pole field and $`Q_{Mm}`$ are the leading multipole moments of $`f`$.

It is often convenient to write multipole fields in terms of partial derivatives of the Coulomb potential centred at $`X_+`$:

$`\varphi _+(x)={\displaystyle \frac{1}{r_+}}={\displaystyle \frac{1}{\sqrt{\rho ^2+(x_3R/2)^2}}},`$ (4.12)
where $`\rho ^2=𝒙^2`$. The function $`_3^{m_3}_2^{m_2}_1^{m_1}\varphi _+(x)`$ is a $`2^M`$-pole field if $`m_1+m_2+m_3=M`$. However, not of all of the fields obtained in this way are independent. We introduce the complex derivatives

$`={\displaystyle \frac{1}{2}}(_1i_2)\text{and}\overline{}={\displaystyle \frac{1}{2}}(_1+i_2)`$ (4.13)
and note that $`\mathrm{\Delta }=_3^2+4\overline{}`$. Then, since $`\mathrm{\Delta }\varphi _+(x)=0`$ we have

$`\overline{}\varphi _+(x)={\displaystyle \frac{1}{4}}_3^2\varphi _+(x).`$ (4.14)
Thus a basis for $`2^M`$-pole fields is given by $`_3^{m_3}^n\overline{}^{\overline{n}}\varphi _+(x)`$ where $`M=m_3+n+\overline{n}`$ and either $`n`$ or $`\overline{n}`$ can be taken to be zero. In appendix A we derive the exact relation between the functions $`_3^{M+m}^m\varphi _+(x)`$, $`Mm<0`$, and $`_3^{Mm}\overline{}^m\varphi _+(x)`$, $`0mM`$, on the one hand and the spherical harmonics $`Y_{Mm}`$ centred at $`X_+`$ on the other. The result is that we have the alternative expansion of an $`2^M`$-pole field

$`f_M(x)={\displaystyle \underset{Mm0}{}}A_{Mm}_3^{M+m}^m\varphi _+(x)+{\displaystyle \underset{1mM}{}}A_{Mm}_3^{Mm}\overline{}^m\varphi _+(x),`$ (4.15)
where the coefficients $`A_{Mm}`$, $`m=M,\mathrm{},M`$ are directly proportional to the multipole moments $`Q_{Mm}`$. It follows from the results in appendix A that

$`A_{Mm}=\sqrt{{\displaystyle \frac{4\pi }{2M+1}}}{\displaystyle \frac{(1)^{M+m}2^m}{\sqrt{(Mm)!(M+m)!}}}Q_{Mm}`$ (4.16)
for $`m0`$ and

$`A_{Mm}=\sqrt{{\displaystyle \frac{4\pi }{2M+1}}}{\displaystyle \frac{(1)^M2^{|m|}}{\sqrt{(Mm)!(M+m)!}}}Q_{Mm}`$ (4.17)
for $`m<0`$. Note in particular that the reality of $`f_M`$ is equivalent to

$`A_{M(m)}=\overline{A}_{Mm}.`$ (4.18)
Multipole fields have a remarkably simple Fourier transform, which will be important for us. We use the representation

$`\varphi _+(x)={\displaystyle \frac{d^2𝒌}{2\pi k}e^{k|x_3\frac{R}{2}|}\mathrm{exp}(i𝒌𝒙)},`$ (4.19)
which can be verified as follows. Exploiting the invariance of $`\varphi _+`$ under rotations in the $`x_1x_2`$-plane we may assume that $`𝒙=(\rho ,0)`$. Using polar coordinates $`(k,\psi )`$ for $`𝒌`$ we first carry out the $`dk`$ integration and then the angular integration:

$`{\displaystyle _0^{2\pi }}{\displaystyle \frac{d\psi }{2\pi }}{\displaystyle _0^{\mathrm{}}}𝑑ke^{k|x_3\frac{R}{2}|}\mathrm{exp}(i𝒌𝒙)`$ $`=`$ $`{\displaystyle _0^{2\pi }}{\displaystyle \frac{d\psi }{2\pi }}{\displaystyle \frac{1}{(i\rho \mathrm{cos}\psi |x_3R/2|)}}`$ (4.20)
$`=`$ $`{\displaystyle _{S^1}}{\displaystyle \frac{dw}{2\pi i}}{\displaystyle \frac{2}{(i\rho w^22|x_3R/2|w+i\rho )}},`$ (4.21)
where we changed variables to $`w=e^{i\psi }`$ in the last line. Expanding the integrand in partial fractions and using the residue theorem then yields the expression (4.12). We compute the Fourier transform of the multipole field (4.15) by differentiating (4.19) under the integral sign. Note that

$`\mathrm{exp}(i𝒌𝒙)={\displaystyle \frac{i}{2}}ke^{i\psi }\mathrm{exp}(i𝒌𝒙)\text{and}\overline{}\mathrm{exp}(i𝒌𝒙)={\displaystyle \frac{i}{2}}ke^{i\psi }\mathrm{exp}(i𝒌𝒙),`$ (4.22)
so we find

$`f_M(x)={\displaystyle \frac{1}{2\pi }}{\displaystyle d^2𝒌e^{k(\frac{R}{2}x_3)}k^{M1}\underset{MmM}{}\left(\frac{i}{2}\right)^{|m|}A_{Mm}e^{im\psi }\mathrm{exp}(i𝒌𝒙)}.`$ (4.23)
Here we have used that $`x_3<\delta `$ so that in particular $`x_3<\frac{R}{2}`$. Thus with the normalisation (4.4) we arrive at the following simple expression for the Fourier transform

$`p_M(𝒌)=4\pi e^{k\frac{R}{2}}k^M{\displaystyle \underset{MmM}{}}\left({\displaystyle \frac{i}{2}}\right)^{|m|}A_{Mm}e^{im\psi }.`$ (4.24)
Remarkably, this function factorises into a $`k`$-dependent part and the function

$`\mathrm{\Theta }(\psi )={\displaystyle \underset{MmM}{}}\left({\displaystyle \frac{i}{2}}\right)^{|m|}A_{Mm}e^{im\psi }`$ (4.25)
of the angle $`\psi `$. The $`k`$-dependent part $`e^{k\frac{R}{2}}k^M`$ is non-zero for $`k0`$ and the function $`\mathrm{\Theta }`$ only vanishes identically if $`A_{Mm}=0`$ for all $`m=M,\mathrm{},M`$, i.e. if the $`2^M`$-pole field is trivial.

### 4.2 The interaction energy of two scalar multipoles

The interaction energy of two multipoles can be expressed in a remarkably compact way. Let

$`\varphi _{}(x)={\displaystyle \frac{1}{|xX_{}|}}`$ (4.26)
be the Coulomb potential centred at $`X_{}=(0,0,R/2)`$ and consider the multipole field

$`g_N(x)={\displaystyle \underset{Nn0}{}}B_{Nn}_3^{N+n}^n\varphi _{}(x)+{\displaystyle \underset{1nN}{}}B_{Nn}_3^{Nn}\overline{}^n\varphi _{}(x)`$ (4.27)
with $`B_{N(n)}=\overline{B}_{Nn}`$. By the same calculation as for $`f_M`$ above we find the Fourier transform of $`g_N`$ in the $`x_1x_2`$-plane:

$`q_N(𝒌)=4\pi e^{k\frac{R}{2}}k^M{\displaystyle \underset{NnN}{}}\left({\displaystyle \frac{i}{2}}\right)^{|n|}B_{Nn}e^{in\psi }.`$ (4.28)
The interaction energy of the two multipole fields $`f_M`$ and $`g_N`$

$`V_{MN}=f_M,g_N`$ (4.29)
can now be computed using the formula (4.9). Using the factorisation property of the Fourier transforms $`p_N`$ and $`q_N`$, it is easy to perform the integration over $`𝒌`$. Assuming without loss of generality that $`MN`$ we first carry out the integration over the angle $`\psi `$ to find

$`V_{MN}=4\pi {\displaystyle _0^{\mathrm{}}}𝑑ke^{kR}k^{N+M}{\displaystyle \underset{m=M}{\overset{M}{}}}2^{2|m|}\overline{A}_{Mm}B_{Nm},`$ (4.30)
where we have used the reality condition for the coefficients $`A_{Mm}`$ and $`B_{Nm}`$. Computing the remaining integral we obtain the final result

$`V_{MN}=4\pi {\displaystyle \frac{(M+N)!}{R^{M+N+1}}}{\displaystyle \underset{m=M}{\overset{M}{}}}2^{2|m|}\overline{A}_{Mm}B_{Nm}.`$ (4.31)
This formula has a number of interesting features. The interaction energy depends only on the separation of the multipoles and on the combination $`_{m=M}^M2^{2|m|}\overline{A}_{Mm}B_{Nm}`$ of the multipole components. As explained in appendix A, the multipole moments $`Q_{Nn}`$ of a $`2^N`$-pole can be thought of as elements of the $`(2N+1)`$-dimensional unitary irreducible representation $`W_N`$ of $`SO(3)`$. The vector $`B`$ with $`2N+1`$ components $`B_{Nn}`$, $`NnN`$ is naturally an element of $`W_N`$. Rotations $`GSO(3)`$ about the centre $`X_{}`$ of the multipole field $`g_N`$ act on the multipole components via $`B_{Nn}_{n^{}=N}^NU_{nn^{}}^N(G)B_{Nn^{}}`$, where $`U^N`$ is a $`(2N+1)`$-dimensional irreducible representation of $`SO(3)`$ (because of the rescaling (4.16) and (4.17) this is not the standard unitary representation). With our assumption that $`MN`$ we can use the multipole components $`A_{Mm}`$, $`MmM`$, to define the linear form

$`F_A:W_N,B{\displaystyle \underset{m=M}{\overset{M}{}}}2^{2|m|}\overline{A}_{Mm}B_{Nm}.`$ (4.32)
By assumption, the components $`A_{Mm}`$ are not all zero, and therefore the map $`F_A`$ is non-degenerate. Writing the formula (4.31) in terms of this map as

$`V_{MN}=4\pi {\displaystyle \frac{(M+N)!}{R^{M+N+1}}}F_A(B),`$ (4.33)
we immediately deduce the following result.

###### Theorem 4.1

The interaction energy of an $`2^M`$-pole and a $`2^N`$-pole separated by a distance $`R`$ is always non-vanishing for some relative orientation of the two multipoles. When such an orientation is chosen, the modulus of the interaction energy decreases with the separation as $`R^{(M+N+1)}`$.

Proof: Assuming without loss of generality that the multipoles are separated along the $`x_3`$-axis and that $`MN`$ we have the formula (4.33) for the interaction energy. Since the map (4.32) defined in terms of the (non-vanishing) multipole components $`A_{Mm}`$ of the $`2^M`$ pole at $`X_+`$ is non-degenerate it has a $`2N`$-dimensional kernel. It then follows from the irreducibility of the $`(2N+1)`$-dimensional representation $`W_N`$ that there exists a rotation $`GSO(3)`$ such that $`U^N(G)B`$ is not in the kernel of $`F_A`$ for some $`G`$. For that $`G`$ we thus have $`F_A(U^N(G)B)=\kappa 0`$ and $`V_{MN}=4\pi \kappa (M+N)!R^{(M+N+1)}`$.

## 5 Attractive forces and existence of minima

The arguments of the previous section apply to the asymptotic pion fields of the Skyrme solitons $`U^{(1)}`$ and $`U^{(2)}`$ discussed in Section 3. In particular we note that the interaction energy $`\mathrm{\Delta }E`$ (3.22) for the leading multipole fields $`u_M`$ and $`v_N`$ is just the sum over iso-components of pairings of the form (4.29)

$`\mathrm{\Delta }E_a=2u_M^a,v_N^a=\mathrm{\hspace{0.17em}2}{\displaystyle _{x_3=0}}𝑑x_1𝑑x_2(u_M^a_3v_N^av_N^a_3u_M^a).`$ (5.1)
Now pick one of the iso-indices, say $`a=1`$, and use iso-rotations of the Skyrme solitons to make sure that the first iso-components $`u_M^1`$ and $`v_N^1`$ are non-vanishing. It then follows from theorem (4.1) that we can make the multipole interaction energy $`\mathrm{\Delta }E_1`$ non-zero by spatial rotations of one of the Skyrme solitons. This fact is the main upshot of the calculations in the previous section and a crucial input for the following argument which was missing in .

Now consider the sum

$`\mathrm{\Delta }E=\mathrm{\Delta }E_1+\mathrm{\Delta }E_2+\mathrm{\Delta }E_3.`$ (5.2)
We would like to show that we can always arrange for $`\mathrm{\Delta }E`$ to be negative by a suitable iso-rotation of one of the Skyrme solitons. We may assume that, possibly after re-labelling the pion fields,

$`\mathrm{\Delta }E_1\mathrm{\Delta }E_2\mathrm{\Delta }E_3,`$ (5.3)
If $`\mathrm{\Delta }E<0`$ we are done, so suppose that $`\mathrm{\Delta }E0`$. Since we know that not all $`\mathrm{\Delta }E_a`$ vanish we can conclude that $`\mathrm{\Delta }E_1>0`$. Now perform an iso-rotation of Skyrme soliton $`1`$ by $`180`$ degrees around the third iso-spin axis. This reverses the sign of $`\pi _1^{(1)}`$ and $`\pi _2^{(1)}`$ and hence also the sign of $`\mathrm{\Delta }E_1`$ and $`\mathrm{\Delta }E_2`$. The new value of $`\mathrm{\Delta }E`$ is

$`\mathrm{\Delta }E`$ $`=`$ $`\mathrm{\Delta }E_1\mathrm{\Delta }E_2+\mathrm{\Delta }E_3`$ (5.4)
$`=`$ $`\mathrm{\Delta }E_1(\mathrm{\Delta }E_2\mathrm{\Delta }E_3)<0`$
since $`\mathrm{\Delta }E_1<0`$ and, with our ordering, $`(\mathrm{\Delta }E_2\mathrm{\Delta }E_3)0`$.

Thus, the contribution $`\mathrm{\Delta }E`$ to the interaction energy of two Skyrme solitons $`U^{(1)}`$ and $`U^{(2)}`$ can always be made less than zero by suitable rotations and iso-rotations of Skyrme soliton 1. It follows from the discussion at the end of Section 3 that for $`|NM|2`$ and sufficiently large separation parameter $`R`$,

$`E[U_R]<E[U^{(1)}]+E[U^{(1)}].`$ (5.5)
We conclude with a few comments on the implications of our result for the question of existence of general Skyrme solitons. As explained in Section 1, Esteban proved the existence of Skyrme solitons of arbitrary degree provided the strict inequality (1.19) holds. Our result (5.5) implies the inequality in those cases where minima exist in the sectors $`l`$ and $`kl`$, and where the associated multipoles have orders which do not differ by more than two. Since monopole fields do not arise in Skyrme solitons, the interaction energy $`\mathrm{\Delta }E`$ dominates at large separation if the leading multipole moments in Skyrme solitons are at most octupoles. As explained at the end of Section 2, the $`B=7`$ Skyrme soliton is believed to have octupoles as leading multipoles, but there is no numerical evidence for leading multipoles of higher order. Unfortunately, it seems very difficult to rule out this possibility in general.

Even if one could prove (or circumvent) the assumption concerning multipoles, the existence of attractive forces between Skyrme solitons is not sufficient to establish the inequality (1.19) for infima. Physically it seems reasonable that the existence of attractive forces should imply the existence of minima in every sector. However, we have not been able to develop this observation into a mathematical proof. Further thoughts and speculations in this direction can be found in .

Acknowledgements

BJS acknowledges financial support through an Advanced Research Fellowship of the Engineering and Physical Sciences Research Council. MAS thanks Rafe Mazzeo for directing him to the unique continuation theorem used in Section 2. NSM is grateful to Walter Kohn for drawing his attention to the formula (4.4).

## Appendix A Spherical Harmonics

In this appendix we derive the relation between the standard spherical harmonics and the following functions on $`^3\{0\}`$ used in the multipole expansion in Section 4.1

$`F_{Nn}(x)=\{\begin{array}{cc}\overline{}^n_3^{Nn}\left(\frac{1}{r}\right)\hfill & \text{if}n0\hfill \\ ^n_3^{N+n}\left(\frac{1}{r}\right)\hfill & \text{if}n<0\hfill \end{array}`$ (A.3)
Here $`N0`$ and $`NnN`$ and $`=\frac{1}{2}(_1i_2)`$. These functions are harmonic in their domain:

$`\mathrm{\Delta }F_{Nn}=0.`$ (A.4)
They are also homogeneous of degree $`(N+1)`$ so that they can be written as

$`F_{Nn}={\displaystyle \frac{1}{r^{N+1}}}\mathrm{\Phi }_{Nn}(\theta ,\phi ),`$ (A.5)
where $`(\theta ,\phi )`$ are the usual spherical coordinates on the two-sphere centred at the origin. Since the Laplace operator takes the following form in spherical coordinates

$`\mathrm{\Delta }={\displaystyle \frac{1}{r}}{\displaystyle \frac{^2}{r^2}}r+{\displaystyle \frac{1}{r^2}}\mathrm{\Delta }_\omega ,`$ (A.6)
where $`\mathrm{\Delta }_\omega `$ is the Laplace operator on the 2-sphere of unit radius, it follows from (A.4) that

$`\mathrm{\Delta }_\omega \mathrm{\Phi }_{Nn}=N(N+1)\mathrm{\Phi }_{Nn}.`$ (A.7)
Define the generator of rotations about the $`3`$-axis

$`J_3=i{\displaystyle \frac{}{\phi }}`$ (A.8)
and express it in terms of complex coordinates $`z=x_1+ix_2`$ and complex derivatives in the $`x_1x_2`$ plane:

$`J_3=z\overline{z}\overline{}.`$ (A.9)
The operator $``$ acts as a raising operator and the operator $`\overline{}`$ acts as a lowering operator for $`J_3`$:

$`[J_3,\overline{}]=\overline{}\text{and}[J_3,]=.`$ (A.10)
Thus if $`\varphi _n`$ is a function on $`^3`$ which is an eigenfunction of $`J_3`$ with eigenvalue $`n`$ then $`\overline{}\varphi _n`$ is an eigenfunction with eigenvalue $`n+1`$ provided it is not zero. Similarly $`\varphi _n`$ is an eigenfunction of $`J_3`$ with eigenvalue $`n1`$ provided it is not zero.

It follows from (A.3) that

$`F_{Nn}=\{\begin{array}{cc}\overline{}^nF_{Nn,0}\hfill & \text{if}n0\hfill \\ ^nF_{N+n,0}\hfill & \text{if}n<0.\hfill \end{array}`$ (A.13)
Noting that, by rotational symmetry about the $`3`$-axis,

$`J_3\mathrm{\Phi }_{N0}=0\text{for all}N,`$ (A.14)
we conclude that

$`J_3F_{Nn}=nF_{Nn}.`$ (A.15)
It follows that

$`\mathrm{\Phi }_{Nn}=r^{N+1}F_{Nn}`$ (A.16)
also satisfies

$`J_3\mathrm{\Phi }_{Nn}=n\mathrm{\Phi }_{Nn}.`$ (A.17)
Thus, to sum up, the $`\mathrm{\Phi }_{Nn}`$ are functions on $`S^2`$, which are eigenfunctions of both the Laplace operator (A.7) and the operator $`J_3`$ with eigenvalues respectively $`N(N+1)`$ and $`n`$. It follows from standard harmonic analysis on $`S^2`$ that they must be proportional to the spherical harmonics $`Y_{Nn}`$.

In the case $`n=0`$ we can determine the proportionality constant by evaluating both $`\mathrm{\Phi }_{N0}`$ and $`Y_{N0}`$ on the positive $`3`$-axis. With the usual normalisation  we find

$`Y_{N0}=\sqrt{{\displaystyle \frac{2N+1}{4\pi }}}{\displaystyle \frac{(1)^N}{N!}}\mathrm{\Phi }_{N0}.`$ (A.18)
The relation between $`Y_{Nn}`$ and $`\mathrm{\Phi }_{Nn}`$ for $`n0`$ is harder to compute. Let us assume initially that $`n>0`$. Starting with the standard expression for the associated Legendre function in terms of Legendre polynomials

$`P_N^n(\mathrm{cos}\theta )=(1)^n\mathrm{sin}^n\theta \left({\displaystyle \frac{d}{d\mathrm{cos}\theta }}\right)^nP_N(\mathrm{cos}\theta )`$ (A.19)
and the expression of the spherical harmonic in terms of associated Legendre function (see e.g.  p. 99) we have the relation

$`Y_{Nn}(\theta ,\phi )=(1)^n\sqrt{{\displaystyle \frac{(Nn)!}{(N+n)!}}}\left(\mathrm{sin}\theta e^{i\phi }\right)^n\left({\displaystyle \frac{}{\mathrm{cos}\theta }}\right)^nY_{N0}(\theta ,\phi ).`$ (A.20)
Then using (A.18) and the definition of $`\mathrm{\Phi }_{N0}`$ we deduce

$`Y_{Nn}(\theta ,\phi )=r^{N+1}\sqrt{{\displaystyle \frac{2N+1}{4\pi }}}{\displaystyle \frac{(1)^{N+n}}{N!}}\sqrt{{\displaystyle \frac{(Nn)!}{(N+n)!}}}z^nD_\theta ^n_3^N\left({\displaystyle \frac{1}{r}}\right),`$ (A.21)
where

$`D_\theta ={\displaystyle \frac{1}{r}}{\displaystyle \frac{}{\mathrm{cos}\theta }}={\displaystyle \frac{x_3}{\rho }}{\displaystyle \frac{}{\rho }}+_3`$ (A.22)
and $`z=\rho e^{i\phi }`$ with $`\rho =r\mathrm{sin}\theta `$ as in the main text of the paper. Then we use the commutation relation

$`[D_\theta ,_3]={\displaystyle \frac{1}{\rho }}{\displaystyle \frac{}{\rho }}`$ (A.23)
to move $`D_\theta `$ past $`_3`$ in (A.21). Noting that $`D_\theta r=0`$ we find

$`Y_{Nn}(\theta ,\phi )=r^{N+1}\sqrt{{\displaystyle \frac{2N+1}{4\pi }}}{\displaystyle \frac{(1)^{N+n}}{\sqrt{(Nn)!(N+n)!}}}z^n\left({\displaystyle \frac{1}{\rho }}{\displaystyle \frac{}{\rho }}\right)^n_3^{Nn}\left({\displaystyle \frac{1}{r}}\right).`$ (A.24)
Now exploit that on any function $`f`$ which only depends on $`\rho `$ and $`x_3`$

$`{\displaystyle \frac{1}{\rho }}{\displaystyle \frac{f}{\rho }}(\rho ,x_3)={\displaystyle \frac{2}{z}}\overline{}f(\rho ,x_3)`$ (A.25)
to conclude

$`Y_{Nn}(\theta ,\phi )=r^{N+1}\sqrt{{\displaystyle \frac{2N+1}{4\pi }}}{\displaystyle \frac{(1)^{N+n}\mathrm{\hspace{0.17em}2}^n}{\sqrt{(Nn)!(N+n)!}}}\overline{}^n_3^{Nn}\left({\displaystyle \frac{1}{r}}\right).`$ (A.26)
Thus we finally arrive at promised relationship between $`Y_{Nn}`$ and $`\mathrm{\Phi }_{Nn}`$, valid for $`n0`$:

$`Y_{Nn}=\sqrt{{\displaystyle \frac{2N+1}{4\pi }}}{\displaystyle \frac{(1)^{N+n}\mathrm{\hspace{0.17em}2}^n}{\sqrt{(Nn)!(N+n)!}}}\mathrm{\Phi }_{Nn}.`$ (A.27)
To deduce the corresponding result for $`n<0`$ we note that $`\mathrm{\Phi }_{Nn}=\overline{\mathrm{\Phi }}_{N(n)}`$ and $`Y_{Nn}=(1)^n\overline{Y}_{N(n)}`$ for $`n<0`$. Thus for $`n<0`$:

$`Y_{Nn}=\sqrt{{\displaystyle \frac{2N+1}{4\pi }}}{\displaystyle \frac{(1)^N\mathrm{\hspace{0.17em}2}^{|n|}}{\sqrt{(Nn)!(N+n)!}}}\mathrm{\Phi }_{Nn}.`$ (A.28)