\section{Introduction}\label{section-intro}
The exactness and the integrality of morphisms of log schemes,
introduced by Kazuya Kato \cite{KK1}, have been widely accepted as
practically useful concepts in log geometry.
For instance, a log smooth morphism is underlain by a flat morphism if
it is integral, while the underlying morphism of a log smooth morphism
in general might be fairly pathological.

In this note we will prove that, under a few minor conditions, any
morphism of fs log schemes can be modified by so-called log
modifications to be an exact, and moreover an integral morphism.
Our main theorem is the following (Theorem \ref{thm-blowupbasechange2}
and Theorem \ref{thm-blowupbasechange3}):

\vspace{1ex}
\noindent
{\bf Theorem.} {\it Let $f\colon X\rightarrow Y$ be a morphism of
locally noetherian fs log schemes.
Suppose the underlying structure of $f$ and $Y$ are quasi-compact.}
\begin{enumerate}\vspace{-1ex}
\item {\it There exists a log-ideal $\mathcal{K}$ of $Y$ such that the base
      change of $f$ by the log blow-up of $Y$ along $\mathcal{K}$ is exact.}
      \vspace{-1.7ex}
\item {\it There exists a commutative diagram}
      $$
      \begin{array}{ccc}
      \til{X}&\longrightarrow&X\\
      \llap{$\scriptstyle{\til{f}}$}\vphantom{\bigg|}\Big\downarrow&&\vphantom{\bigg|}\Big\downarrow
      \rlap{$\scriptstyle{f}$}\\
      \til{Y}&\longrightarrow&Y
      \end{array}
      $$
      {\it of fs log schemes such that the horizontal arrows are log
      modifications and $\til{f}$ is integral.}\vspace{-1ex}
\end{enumerate}

The log blow-up is also the concept which was already studied by Kazuya
Kato in his unpublished works (cf.\ \cite{KK2}\cite{KK3}). Roughly
speaking, a log blow-up is the composite of a usual blow-up centered
at an ideal which comes from an ideal of the log strtucture
(log-ideal) followed by normalization. The log modification is a
certain class of universally surjective log \'etale morphisms, which of
course includes log blow-up's.

Notice that the second part of the theorem above (especially for a log
smooth morphism $f$) can be compared to a weak log-version of the
famous flattening theorem by Raynaud and Gruson \cite{RG1}.

\vspace{1ex}
{\it Acknowledgments.}
This work owes much to conversations with Professor Richard Pink. The
author expresses gratitude to him. The author thanks Max-Planck-Institut
f\"ur Mathematik Bonn for the nice hospitality.

\vspace{1ex}
{\it Notation and conventions.}
All sheaves are considered with respect to the \'etale topology. Log
schemes are denoted simply by capital letters. For a log scheme $X$, we
denote by $\under{X}$, $\O_X$, $\mathcal{M}_X$, and $\ovl{\mathcal{M}}_X$ the
underlying scheme, the structure sheaf of $\under{X}$, the log structure
of $X$, and the characteristic $(=\mathcal{M}_X/\Oi_X)$ of $X$, respectively.
For a commutative and unitary monoid $P$ the associated group (resp.\
the subgroup of invertible elements) is denoted by $\gp{P}$ (resp.\
$\inv{P}$).

\setcounter{equation}{0}
\section{Exact morphism and integral morphism}\label{section-imtegral}
In this section we recall the definition of exactness and integrality,
and basic properties of them, especially in relation to charts; we also
recall the notion of neat and good charts which is useful in the next
section.
\begin{dfn}\label{dfn-exacthomo}
\begin{enumerate}\vspace{-1ex}
\item
A homomorphism $Q\rightarrow P$ of monoids
is said to be {\it exact} if the diagram
$$
\begin{array}{ccc}
Q&\longrightarrow&P\\
\vphantom{\bigg|}\Big\downarrow&&\vphantom{\bigg|}\Big\downarrow\\
\gp{Q}&\longrightarrow&\gp{P}
\end{array}
$$
is Cartesian in the category of sets.
\vspace{-1.7ex}
\item
A morphism $f\colon X\rightarrow Y$ of log schemes is said to be {\it exact
at} $x\in X$ if the homomorphism $\ovl{\mathcal{M}}_{Y,\ovl{y}}\rightarrow
\ovl{\mathcal{M}}_{X,\ovl{x}}$ is exact, where $y=f(x)$.
({\it Here and hereafter $\ovl{\cdot}$ means the separable closure of points.})
The morphism $f$ is said to be {\it exact} if it is exact at any point of
$X$.
\vspace{-1ex}
\end{enumerate}
\end{dfn}

If a homomorphism $h\colon Q\rightarrow P$ of monoids is exact and
$\inv{Q}=1$, then $h$ is injective.
In particular, if a morphism $f\colon X\rightarrow Y$ is exact at $x\in X$,
then $\ovl{\mathcal{M}}_{Y,\ovl{y}}\rightarrow\ovl{\mathcal{M}}_{X,\ovl{x}}$ is injective.

Let $P$ be a monoid and $A$ a ring.
We denote by $\mathop{\mathrm{S}}\nolimits(P)_A$ the log scheme underlain by the scheme $\Spec A[P]$
endowed with the log structure induced by $P\rightarrow A[P]$.
We write, in particular, $\mathop{\mathrm{S}}\nolimits(P)=\mathop{\mathrm{S}}\nolimits(P)_{\Z}$.
The log structure of $\mathop{\mathrm{S}}\nolimits(P)$ is denoted by $\mathcal{M}_P$.
(Caution: Even if $Q\rightarrow P$ is exact, $\mathop{\mathrm{S}}\nolimits(P)\rightarrow\mathop{\mathrm{S}}\nolimits(Q)$ is
not necessarily exact.)

Let $X$ be a log scheme.
A chart of $X$ modeled on a monoid $P$ is a strict morphism $\la\colon
X\rightarrow\mathop{\mathrm{S}}\nolimits(P)$ of log schemes.
The chart $\la$ is said to be good at a point $x\in X$ if $P\rightarrow
\ovl{\mathcal{M}}_{X,\ovl{x}}$ is an isomorphism.
Good charts always exist at any point, provided $X$ is fs.

Let $f\colon X\rightarrow Y$ be a morphism of log schemes.
A chart of $f$ modeled on a homomorphism $h\colon
Q\rightarrow P$ of monoids is a commutative diagram
$$
\begin{array}{ccc}
X&\longrightarrow&\mathop{\mathrm{S}}\nolimits(P)\\
\llap{$\scriptstyle{f}$}\vphantom{\bigg|}\Big\downarrow&&\vphantom{\bigg|}\Big\downarrow
\rlap{$\scriptstyle{\mathop{\mathrm{S}}\nolimits(h)}$}\\
Y&\longrightarrow&\mathop{\mathrm{S}}\nolimits(Q)
\end{array}
$$
with the horizontal arrows strict.
This chart of $f$ is said to be {\it neat} at $x\in X$ if $\gp{h}$ is
injective, $\gp{P}/\gp{Q}\rightarrow\ovl{\mathcal{M}}_{X/Y,\ovl{x}}$ is an
isomorphism, and $P\rightarrow\ovl{\mathcal{M}}_{X,\ovl{x}}$ is exact.
This chart of $f$ is said to be {\it good} at $x\in X$ if the charts
$X\rightarrow \mathop{\mathrm{S}}\nolimits(P)$ and $Y\rightarrow\mathop{\mathrm{S}}\nolimits(Q)$ are good at $x$ and $f(x)$,
respectively.
\begin{lem}\label{lem-neatchart}
(A.\ Ogus \cite[Theorem 2.13]{AG1})\
{\it Let $f\colon X\rightarrow Y$ be a morphism of fine log schemes and
$x\in X$ a point.
Suppose there exists a chart} $Y\rightarrow\mathop{\mathrm{S}}\nolimits(Q)$ {\it of $Y$ modeled on
a fine monoid $Q$.}
{\it Then $f$ has fppf locally a neat chart at $x$ modeled on a homomorphism
of fine monoids extending the given chart} $Y\rightarrow\mathop{\mathrm{S}}\nolimits(Q)$.
{\it Furthermore, if the order of the torsion part of}
$\gp{\ovl{\mathcal{M}}}_{X/Y,\ovl{x}}$ {\it is invertible in the residue field
$k(\ovl{x})$, then the neat chart as above can be found \'etale locally.}
\end{lem}

\begin{cor}\label{cor-goodchart}
{\it If $f\colon X\rightarrow Y$ is a morphism of fs log schemes exact at
$x\in $X, then $f$ allows fppf locally a good chart at $x$.}
\end{cor}

\noindent{\bf Proof.}\hspace{2mm}
Take \'etale locally a good chart $Y\rightarrow\mathop{\mathrm{S}}\nolimits(Q)$ at $y=f(x)$ and
extend it fppf locally to a neat chart $(X\rightarrow\mathop{\mathrm{S}}\nolimits(P),
Y\rightarrow\mathop{\mathrm{S}}\nolimits(Q),Q\rightarrow P)$ at $x$.
Since $\ovl{\mathcal{M}}_{Y,\ovl{y}}\rightarrow\ovl{\mathcal{M}}_{X,\ovl{x}}$ is injective,
the chart $X\rightarrow\mathop{\mathrm{S}}\nolimits(P)$ is automatically good at $x$.
\hfill$\square$
\begin{dfn}\label{dfn-integralhomo}
A homomorphism $h\colon Q\rightarrow P$ of integral monoids is said to be
{\it integral} if for any integral monoid $Q'$ and a homomorphism
$Q'\rightarrow Q$, the push-out $P\oplus_{Q}Q'$ in the category of monoids
is integral.
\end{dfn}

It is known that the integrality of $h\colon Q\rightarrow P$ is equivalent
to the following condition:
If $a_1,a_2\in Q$, $b_1,b_2\in P$ and $h(a_1)b_1=h(a_2)b_2$,
there exist $a_3,a_4\in Q$ and $b\in P$ such that $b_1=h(a_3)b$,
$b_2=h(a_4)b$, and $a_1a_3=a_2a_4$ (cf.\ \cite[(4.1)]{KK1}).
It is also known that $h\colon Q\rightarrow P$ is integral and injective
if and only if the map $\Z[h]\colon\Z[Q]\rightarrow\Z[P]$ is flat.
\begin{lem}\label{lem-integralexact}
\begin{enumerate}\vspace{-1ex}
\item {\it Let $h\colon Q\rightarrow P$ be an integral
homomorphism of integral monoids such that} $h^{-1}(\inv{P})=\inv{Q}$.
{\it Then $h$ is exact.}
\vspace{-1.7ex}
\item {\it Consider the commutative diagram}
$$
\begin{array}{ccc}
P'&\stackrel{a}{\longrightarrow}&P\\
\llap{$\scriptstyle{h'}$}\vphantom{\bigg|}\Big\uparrow&&\vphantom{\bigg|}\Big\uparrow
\rlap{$\scriptstyle{h}$}\\
Q'&\longrightarrow&Q
\end{array}
$$
{\it of integral monoids such that the horizontal arrow are surjective.}
{\it Then}:
\begin{enumerate}\vspace{-1ex}
\item {\it If $h'$ is integral and there exists a submonoid $P''\subseteq P'$
such that} $P'/P''\cong P$, {\it then $h$ is integral.}
\vspace{-1ex}
\item {\it If $h$ is integral and the horizontal arrows are exact, then
$h'$ is integral.}
\vspace{-1ex}
\end{enumerate}
({\it Here the quotient of a monoid $M$ by a submonoid $N$ is defined by the
coset space $M/\sim$ with the induced monoid structure, where $a\sim b$ if
and only if there exist $c,d\in N$ such that $ac=bd$.})
\vspace{-1ex}
\end{enumerate}
\end{lem}

\noindent{\bf Proof.}\hspace{2mm}
Straightforward.
\hfill$\square$
\begin{dfn}\label{dfn-integralmorphism}
A morphism $f\colon X\rightarrow Y$ of integral log schemes is said to be
{\it integral at} $x\in X$ if $\ovl{\mathcal{M}}_{Y,\ovl{y}}\rightarrow
\ovl{\mathcal{M}}_{X,\ovl{x}}$ is integral, where $y=f(x)$.
We say $f$ is {\it integral} if it is integral at any point of $X$.
\end{dfn}

Note that the integrality of $\ovl{\mathcal{M}}_{Y,\ovl{y}}\rightarrow
\ovl{\mathcal{M}}_{X,\ovl{x}}$ is equivalent to that of
$\mathcal{M}_{Y,\ovl{y}}\rightarrow\mathcal{M}_{X,\ovl{x}}$.
By Lemma \ref{lem-integralexact} (1), integral morphisms are exact.
If $h\colon Q\rightarrow P$ is integral, then by Lemma \ref{lem-integralexact}
(2) (a), the morphism $\mathop{\mathrm{S}}\nolimits(h)\colon\mathop{\mathrm{S}}\nolimits(P)\rightarrow\mathop{\mathrm{S}}\nolimits(Q)$ is integral.

In what follows, all push-outs are taken in the category of fs monoids,
unless otherwise explicitely specified.
Accordingly, all fiber products of fs log schemes are taken in the category
fs log schemes, unless otherwise specified.

Let $P\leftarrow Q\rightarrow Q'$ be homomorphisms of fs monoids and
$\gp{P}\leftarrow\gp{Q}\rightarrow\gp{Q'}$ the induced morphisms, which
gives rise to the push-out $\gp{P}\oplus_{\gp{Q}}\gp{Q'}$ in the
category of finitely generated abelian groups.
\begin{lem}\label{lem-pushoutfs}
{\it The push-out of $P\leftarrow Q\rightarrow Q'$ in the category of fs
monoids is given by the smallest
fs submonoid in} $\gp{P}\oplus_{\gp{Q}}\gp{Q'}$ {\it containing
those images of $P$ and $Q'$, or what amounts to the same, the saturation of
the submonoid in} $\gp{P}\oplus_{\gp{Q}}\gp{Q'}$ {\it generated by those
images of $P$ and $Q'$.}
\end{lem}

\noindent{\bf Proof.}\hspace{2mm}
It is enough to prove that the push-out in the category of fine monoids is
the submonoid in $\gp{P}\oplus_{\gp{Q}}\gp{Q'}$ generated by the
images of $P$ and $Q'$, which is easy to verify.
\hfill$\square$
\begin{cor}\label{cor-pushoutfs}
{\it For a homomorphism $Q\rightarrow P$ of fs monoids such that $\gp{Q}\cong
\gp{P}$, the push-out of
$P\leftarrow Q\rightarrow P$ is given by $P$.}
\end{cor}

We denote the push-out of $P\leftarrow Q\rightarrow Q'$ by $P_{Q'}$ if
the monoid $Q$ is evident in the context.
\begin{lem}\label{lem-exactstable}
{\it For an exact homomorphism $Q\rightarrow P$ of fs monoids
and a homomorphism $Q\rightarrow Q'$ of fs monoids,
the induced homomorphism $Q'\rightarrow P_{Q'}$ is exact.}
\end{lem}

\noindent{\bf Proof.}\hspace{2mm}
It is straightfoward to see that the induced homomorphism $Q'\rightarrow
P'$ is exact, where $P'$ is the submonoid in $\gp{P}\oplus_{\gp{Q}}\gp{{Q'}}$
generated by the images of $P$ and $Q'$.
The monoid $P_{Q'}$ is the saturation of $P'$.
Since $Q'$ is saturated, we see that $Q'\rightarrow P_{Q'}$ is also exact.
\hfill$\square$
\begin{cor}\label{cor-exactstable}
{\it Exactness is stable under base change.}
\end{cor}
\setcounter{equation}{0}
\section{Log blow-up and log modification}\label{section-logblowup}
\begin{dfn}\label{dfn-idealmon}
Let $Q$ be a monoid. A subset $K\subseteq Q$ (possibly empty) is said to
be an ideal of $Q$ if $ax\in K$ for any $a\in Q$ and $x\in K$.
\end{dfn}

Any ideal of a finitely generated monoid is finitely generated as an ideal.
Products and intersections of ideals are obviously ideals.
The correspondence $K\mapsto\ovl{K}\colon=\pi(K)\cdot\ovl{Q}$ gives rise to
the natural bijection between the set of all ideals in $Q$ and the set of all
ideals in $\ovl{Q}$, where $\pi\colon Q\rightarrow\ovl{Q}=Q/\inv{Q}$ is the
natural projection.
\begin{ntn}
Let $Q$ be an fs monoid.
For a subset $E\subseteq\gp{Q}$ we denote by $Q\langle E\rangle$ the
smallest fs submonoid in $\gp{Q}$ which contains $Q$ and $E$.
\end{ntn}
\begin{dfn}
Let $X$ be a log scheme.
A {\it log-ideal} of $X$ is a sheaf of ideals $\mathcal{K}$ of $\mathcal{M}_X$.
\end{dfn}
\begin{ntn}
\begin{enumerate}\vspace{-1ex}
\item For an ideal $K\subseteq Q$ we denote by $\til{K}$ the log-ideal of
      $\mathop{\mathrm{S}}\nolimits(Q)$ generated by the image of $K$.
\vspace{-1.7ex}
\item Let $f\colon X\rightarrow Y$ be a morphism of log schemes and
      $\mathcal{K}$ a log-ideal of $Y$. Then we denote by $f^{\bullet}\mathcal{K}$
      the log-ideal of $X$ generated by the image of $f^{-1}\mathcal{K}$.
\vspace{-1ex}
\end{enumerate}
\end{ntn}

Clearly, for $f\colon X\rightarrow Y$ and $g\colon Y\rightarrow Z$ and a
log-ideal $\mathcal{K}$ of $Z$, we have $(g\circ f)^{\bullet}\mathcal{K}=f^{\bullet}
g^{\bullet}\mathcal{K}$.

A log scheme $X$ is called {\it locally noetherian} if the underlying scheme
$\under{X}$ is locally noetherian.
\begin{lem}\label{lem-coherentideal}
(A.\ Ogus \cite[Proposition 2.6]{AG1})\
{\it Let $X$ be a locally noetherian fine log scheme and $\mathcal{K}$ a log-ideal.
Then the following conditions are equivalent}:
\begin{enumerate}\vspace{-1ex}
\item {\it For any $x,y\in X$ such that $x$ belongs to the topological closure
of $\{y\}$, the natural map $\mathcal{K}_{\ovl{x}}\rightarrow\mathcal{K}_{\ovl{y}}$ is
surjective.}
\vspace{-1.7ex}
\item {\it For any point $x\in X$ and any \'etale neighborhood $\iota\colon
U\rightarrow X$ of $x$ with chart} $\la\colon U\rightarrow\mathop{\mathrm{S}}\nolimits(P)$ {\it set}
$K=\beta_{\ovl{x}}^{-1}(\mathcal{K}_{\ovl{x}})$, {\it where} $\beta\colon
P\rightarrow\mathcal{M}_X|_U$.
{\it Then replacing $U$ by a sufficiently small \'etale neighborhood of $x$,
we have}
$\iota^{\bullet}\mathcal{K}=\la^{\bullet}\til{K}$.
\vspace{-1ex}
\end{enumerate}
\end{lem}

\begin{dfn}\label{dfn-coherentideal}
Let $X$ be a locally noetherian fine log scheme.
Then a log-ideal $\mathcal{K}$ of $X$ is said to be {\it coherent} if it
satisfies the conditions in Lemma \ref{lem-coherentideal}.
\end{dfn}
\begin{dfn}\label{dfn-locallyprincipal}
A coherent log-ideal $\mathcal{K}$ of $X$ is said to be {\it locally principal}
if it is \'etale locally generated by a single global section.
\end{dfn}
\begin{dfn}\label{dfn-logblowupuniversal}
Let $X$ be a locally noetherian fs log scheme and $\mathcal{K}$ a coherent log-ideal
of $X$.
Then the {\it log blow-up} $\rho_{\mathcal{K}}\colon\mathop{\mathrm{B}}\nolimits_{\mathcal{K}}(X)\rightarrow X$
of $X$ along $\mathcal{K}$ is a morphism of locally noetherian fs log schemes
satisfying the following conditions:
\begin{enumerate}\vspace{-1ex}
\item The log ideal $\rho_{\mathcal{K}}^{\bullet}\mathcal{K}$ of $\mathop{\mathrm{B}}\nolimits_{\mathcal{K}}(X)$ is
locally principal.\vspace{-1.7ex}
\item For any morphism $f\colon X'\rightarrow X$ of locally
noetherian fs log schemes such that $f^{\bullet}\mathcal{K}$ is locally principal,
there exists a unique morphism $\til{f}\colon X'\rightarrow\mathop{\mathrm{B}}\nolimits_{\mathcal{K}}(X)$
such that $f=\rho_{\mathcal{K}}\circ\til{f}$.
\vspace{-1ex}
\end{enumerate}
\end{dfn}

The log blow-up along $\mathcal{K}$ is unique up to isomorphisms if it exists.
We are going to construct $\rho_{\mathcal{K}}\colon\mathop{\mathrm{B}}\nolimits_{\mathcal{K}}(X)\rightarrow X$.
We should first start with the case $X=\mathop{\mathrm{S}}\nolimits(P)$ for an fs monoid $P$.

The following construction is due to Kazuya Kato \cite{KK3}:
Let $K\subseteq P$ be an ideal. Then the log blow-up
$\mathop{\mathrm{B}}\nolimits_K(P)\rightarrow\mathop{\mathrm{S}}\nolimits(P)$ of $\mathop{\mathrm{S}}\nolimits(P)$ along $K$ is constructed as follows:
Let $I(K)$ be the ideal of $\Z[P]$ generated by $K$.
Consider the natural morphism $\mathop{\mathrm{Proj}}\nolimits\oplus_nI(K)^n\rightarrow
\Spec\Z[P]$.
$\mathop{\mathrm{Proj}}\nolimits\oplus_nI(K)^n$ has the affine open covering
$$
\mathop{\mathrm{Proj}}\nolimits\oplus_nI(K)^n=\bigcup_{p\in K}\Spec\Z[P\langle p^{-1}K
\rangle].
$$
(It is enough for $p$ to run over the generators of $K$.)
The canonical log structures given by $P\langle p^{-1}K\rangle\rightarrow
\Z[P\langle p^{-1}K\rangle]$ obviously glue to a fine log structure on
$\mathop{\mathrm{Proj}}\nolimits\oplus_nI(K)^n$.
The fs log scheme $\mathop{\mathrm{B}}\nolimits_K(P)$ is defined by the saturation of this resulting
log scheme.
The morphisms $\mathop{\mathrm{S}}\nolimits(P\langle p^{-1}K\rangle)\rightarrow\mathop{\mathrm{S}}\nolimits(P)$ give rise to
the desired morphism $\mathop{\mathrm{B}}\nolimits_K(P)\rightarrow\mathop{\mathrm{S}}\nolimits(P)$, which turns out to be the
log blow-up in this case.

Let $X$ be an fs log scheme underlain by a locally noetherian scheme and
$\mathcal{K}$ a coherent log-ideal.
Suppose $X$ has a chart $\la\colon X\rightarrow\mathop{\mathrm{S}}\nolimits(P)$ modeled on an fs monoid
$P$, and that the log-ideal $\mathcal{K}$ is generated by global sections.
Take a point $x\in X$ and set $K=\beta_{\ovl{x}}^{-1}(\mathcal{K}_{\ovl{x}})$, where
$\beta\colon P\rightarrow\mathcal{M}_X$.
Since $\mathcal{K}$ is generated by global sections, the ideal $K$ does not depend
on the choice of $x$.
In this situation we define the log blow-up $\mathop{\mathrm{B}}\nolimits_{\mathcal{K}}(X)\rightarrow X$ by
$$
\mathop{\mathrm{B}}\nolimits_{\mathcal{K}}(X)\colon=X\times_{\mathop{\mathrm{S}}\nolimits(P)}\mathop{\mathrm{B}}\nolimits_K(P)\longrightarrow X.
$$
It can be checked without so much pain that this morphism satisfies the
conditions in Definition \ref{dfn-logblowupuniversal}.
In particular, it does not depend on the choice of the chart $\la$.

In general we take an \'etale covering $\bigcup_{i\in I}U_i\rightarrow X$
such that each $U_i$ has allows a chart $U_i\rightarrow\mathop{\mathrm{S}}\nolimits(P_i)$
and $\mathcal{K}|_{U_i}$ is generated by
global sections.
Then from what we have seen above the local log blow-up's
$\mathop{\mathrm{B}}\nolimits_{\mathcal{K}|_{U_i}}(U_i)\rightarrow U_i$ glue to a morphism of log schemes
$\mathop{\mathrm{B}}\nolimits_{\mathcal{K}}(X)\rightarrow X$ which is nothing but the log
blow-up of $X$ along the coherent log-ideal $\mathcal{K}$.
Note that log blow-up's are log \'etale.

The following lemmas are easy to see by the universality of log blow-up's:
\begin{lem}\label{lem-logblowupselfint}
{\it Let $X$ be a locally noetherian fs log scheme and $\mathcal{K}$ a log-ideal
of $X$. Then we have}
$\mathop{\mathrm{B}}\nolimits_{\mathcal{K}}(X)\times_X\mathop{\mathrm{B}}\nolimits_{\mathcal{K}}(X)\cong\mathop{\mathrm{B}}\nolimits_{\mathcal{K}}(X)$.
\end{lem}
\begin{lem}\label{lem-logblowupcommutative}
{\it Let $\mathcal{K}$ and $\mathcal{K}'$ be coherent log-ideals of $X$. Then we have
a canonical isomorphisms}
$$
\mathop{\mathrm{B}}\nolimits_{\mathcal{K}\ast\mathcal{K}'}(X)\cong\mathop{\mathrm{B}}\nolimits_{\rho_{\mathcal{K}}^{\bullet}\mathcal{K}'}
(\mathop{\mathrm{B}}\nolimits_{\mathcal{K}}(X))\cong\mathop{\mathrm{B}}\nolimits_{\rho_{\mathcal{K}'}^{\bullet}\mathcal{K}}
(\mathop{\mathrm{B}}\nolimits_{\mathcal{K}'}(X))
$$
{\it over $X$, where $\mathcal{K}\ast\mathcal{K}'$ denote the log-ideal of $X$
generated by $\mathcal{K}$ and $\mathcal{K}'$.}
\end{lem}
\begin{exa}\label{exa-logblowuptoric}
It is worth while describing log blow-up's in the toric case:
Let $\Sigma$ be a fan of $N\cong\Z^r$ and $T_{\Sigma}$ the assiciated
toric variety over a ring $A$.
For any cone $\sigma\in\Sigma$, the affine patch $T_{\sigma}=
\Spec A[\sigma^{\vee}\cap M]$ (where $M=\Hom_{\Z}(N,\Z)$) has the canonical
log structure which gives rise to $\mathop{\mathrm{S}}\nolimits(\sigma^{\vee}\cap M)_A$.
These patches glue to an fs log scheme which is also denoted by $T_{\Sigma}$.
Let $T=\Spec A[M]$ be the torus included in $T_{\Sigma}$ as a Zariski open
set and $i\colon T\hookrightarrow T_{\Sigma}$ the open immersion.
Then $T$ is exactly the locus on which the log structure is trivial.
Furthermore, the log structure $\mathcal{M}_{T_{\Sigma}}$ on $T_{\Sigma}$ is
isomorphic to
$$
\O_{T_{\Sigma}}\cap i_{\ast}\Oi_T\lhook\joinrel\longrightarrow\O_{T_{\Sigma}}.
$$

Log blow-up's of toric varieties coincides with those described
in \cite[I, \S 2]{KKM}:
For a log-ideal $\mathcal{K}$ of $T_{\Sigma}$, the ideal
$\mathcal{I}(\mathcal{K})$ of
$\O_{T_{\Sigma}}$ generated by $\mathcal{K}$ is invariant under the
action of $T$.
Conversely, to a $T$-invariant ideal of $\O_{T_{\Sigma}}$ we can associate
the log ideal of $T_{\Sigma}$ by taking the intersection with
$\mathcal{M}_{T_{\Sigma}}$.
This correspondence gives rise to the bijection between the set of all
log-ideals of $T_{\Sigma}$ and that of all $T$-invariant ideals of
$\O_{T_{\Sigma}}$.
If $\mathcal{I}$ is a $T$-invariant ideal of $\O_{T_{\Sigma}}$, the log blow-up
along $\mathcal{I}$ is given by a certain subdivision
$\til{\Sigma}\rightarrow\Sigma$
of the fan $\Sigma$; the fan $\til{\Sigma}$ is the subdivision consisting of
cones on which the associated piecewise linear function
 $\mathrm{ord}\mathcal{I}$.
\end{exa}

Any log blow-up of a toric variety therefore comes from subdivision of
the fan. But a morphism of toric varieties by subdivision is not
necessarily a log blow-up.
\begin{lem}\label{lem-blowupbasechange1}
{\it Let $Q\rightarrow P$ be a homomorphism of fs monoids such that}
$\gp{Q}\cong\gp{P}$.
{\it Then there exists an ideal $K\subseteq Q$ and a strict open immersion}
$\mathop{\mathrm{S}}\nolimits(P)\hookrightarrow\mathop{\mathrm{B}}\nolimits_K(Q)$ {\it such that}
$$
\begin{array}{ccc}
\mathop{\mathrm{S}}\nolimits(P)&=&\mathop{\mathrm{S}}\nolimits(P)\\
\vphantom{\bigg|}\Big\downarrow&&\vphantom{\bigg|}\Big\downarrow\\
\mathop{\mathrm{S}}\nolimits(Q)&\longleftarrow&\mathop{\mathrm{B}}\nolimits_K(Q)
\end{array}
$$
{\it is a Cartesian diagram of fs log schemes.}
\end{lem}

\noindent{\bf Proof.}\hspace{2mm}
Write $T=\mathop{\mathrm{S}}\nolimits(P)$ and $S=\mathop{\mathrm{S}}\nolimits(Q)$.
We can set $P=Q\langle\frac{a_1}{b},\ldots,\frac{a_r}{b}\rangle$ in
$\gp{Q}\cong\gp{P}$, where $b,a_1,\ldots,a_r\in Q$.
Let $K=(b,a_1,\ldots,a_r)$ be the ideal of $Q$.
Then the log blow-up $\mathop{\mathrm{B}}\nolimits_K(Q)$ has the covering
$$
\mathop{\mathrm{B}}\nolimits_K(Q)=\bigcup_{i=0,\ldots,r}V_i,
$$
where
$$
V_0=\mathop{\mathrm{S}}\nolimits(Q\langle\frac{a_1}{b},\ldots,\frac{a_r}{b}\rangle),
$$
which is isomorphic to $T$, and
$$
V_j=\mathop{\mathrm{S}}\nolimits(Q\langle\frac{b}{a_j},\frac{a_1}{a_j},\ldots,\frac{a_r}{a_j}
\rangle)
$$
for $j>0$.
By Lemma \ref{lem-pushoutfs} we have $T\times_SV_0\cong T$ and
$T\times_SV_j\cong\mathop{\mathrm{S}}\nolimits(P\langle\frac{b}{a_j}\rangle)$ for $j>0$.
Since we can check $T\times_S(V_0\cap V_j)=T\times_SV_j$, we have
$T\times_S\mathop{\mathrm{B}}\nolimits_K(Q)\cong T$.
\hfill$\square$
\begin{thm}\label{thm-blowupbasechange2}
{\it Let $f\colon X\rightarrow Y$ be a morphism of locally noetherian fs log
schemes.
Suppose $\under{f}$ and $\under{Y}$ are quasi-compact.
Then there exists a log-ideal $\mathcal{K}$ of $Y$
such that the base change} $X_{\mathop{\mathrm{B}}\nolimits_{\mathcal{K}}(Y)}\rightarrow\mathop{\mathrm{B}}\nolimits_{\mathcal{K}}(Y)$
{\it by the log blow-up of $Y$ along $\mathcal{K}$ is exact.}
\end{thm}

\noindent{\bf Proof.}\hspace{2mm}
First we assume that $\gp{\ovl{\mathcal{M}}}_{X/Y}$ is constant and that there
exists a chart $Y\rightarrow\mathop{\mathrm{S}}\nolimits(Q)$.
Then by Lemma \ref{lem-neatchart}
there exists an fppf covering $\{X_i\rightarrow X\}$ of $X$ consisting of
finitely many members such that each $X_i\rightarrow Y$ allows a chart
$(\la_i\colon X_i\rightarrow\mathop{\mathrm{S}}\nolimits(P_i),\mu\colon Y\rightarrow\mathop{\mathrm{S}}\nolimits(Q),
h_i\colon Q\rightarrow P_i)$
such that $h_i$ is injective and $\gp{\ovl{\mathcal{M}}}_{X/Y}\cong\gp{P}_i/\gp{Q}$.
Set $Q_i=P_i\cap\gp{Q}$. Then $Q\rightarrow Q_i$ yields $\gp{Q}\cong\gp{Q}_i$.
Set $Y_i=Y\times_{\mathop{\mathrm{S}}\nolimits(Q)}\mathop{\mathrm{S}}\nolimits(Q_i)$ for each $i$.
By Lemma \ref{lem-blowupbasechange1} we have an ideal $K_i$ of $Q$ such
that $\mathop{\mathrm{S}}\nolimits(Q_i)\hookrightarrow\mathop{\mathrm{B}}\nolimits_{K_i}(Q)$ gives the base change.
Set $\mathcal{K}_i=\mu^{\bullet}\til{K}_i$ and consider the log blow-up
$\mathop{\mathrm{B}}\nolimits_{\mathcal{K}_i}(Y)\rightarrow Y$.
Then we have $X_i\times_Y\mathop{\mathrm{B}}\nolimits_{\mathcal{K}_i}(Y)\cong X_i$, since
$Y_i\times_Y\mathop{\mathrm{B}}\nolimits_{\mathcal{K}_i}(Y)\cong Y_i$.
Since $Y_i\rightarrow\mathop{\mathrm{B}}\nolimits_{\mathcal{K}_i}(Y)$ is strict, we have
$\gp{\ovl{\mathcal{M}}}_{X_i/\mathop{\mathrm{B}}\nolimits_{\mathcal{K}_i}(Y)}\cong\gp{\ovl{\mathcal{M}}}_{X_i/Y}$,
which is constant and is isomorphic to $\gp{P}_i/\gp{Q}_i$.
Hence it follows from the exactness of $Q_i\rightarrow P_i$ that
$X_i\rightarrow Y_i$ is exact, and hence $X_i\rightarrow\mathop{\mathrm{B}}\nolimits_{\mathcal{K}_i}(Y)$ is
exact.
Let $\mathcal{K}$ be the log-ideal of $Y$ generated by $\mathcal{K}_i$'s.
Since the exactness is stable under base changes, the base change
$X_{\mathop{\mathrm{B}}\nolimits_{\mathcal{K}}(Y)}\rightarrow\mathop{\mathrm{B}}\nolimits_{\mathcal{K}}(Y)$ restricted on each $X_i$
is exact, and hence the whole
$X_{\mathop{\mathrm{B}}\nolimits_{\mathcal{K}}(Y)}\rightarrow\mathop{\mathrm{B}}\nolimits_{\mathcal{K}}(Y)$ is exact.
Hence we have proved the theorem in this particular case.

In general we may still assume $Y$ allows a chart $Y\rightarrow\mathop{\mathrm{S}}\nolimits(Q)$.
The sheaf $\gp{\ovl{\mathcal{M}}}_{X/Y}$ is constructible.
We look at each locally closed set of $\under{X}$ on which
$\gp{\ovl{\mathcal{M}}}_{X/Y}$ is constant, and take a log blow-up to make $f$ exact
on this locus.
Finally, we gather up the resulting log-ideal of $Y$ and then take
sufficiently large log blow-up.
\hfill$\square$
\begin{dfn}\label{dfn-logmodification}
A morphism $\rho\colon\til{X}\rightarrow X$ of locally noetherian
fs log schemes is called a {\it log modification} if for any point
$x\in X$ there exists an \'etale neighborhood $U\rightarrow X$ of
$x$, and a log blow-up $Z\rightarrow\til{X}_U$ such that the composite
$Z\rightarrow\til{X}_U\rightarrow U$ is also a log blow-up.
\end{dfn}

Note that if $X$ is quasi-compact the morphism
$\rho\colon\til{X}\rightarrow X$ is a log modification if and only if
there exists a log blow-up $Z\rightarrow\til{X}$ such that the composite
$Z\rightarrow\til{X}\rightarrow X$ is also a log blow-up.

The above definition of log modifications is motivated by the following
facts:
\begin{itemize}\vspace{-1ex}
\item A morphism $\rho\colon\til{X}\rightarrow X$ of quasi-projective
varieties over a field is a modification (i.e.\ proper and birational
morphism) if and only if there exists a blow-up $Z\rightarrow\til{X}$
such that the composite $Z\rightarrow X$ is a blow-up; this follows
easily from \cite[Theorem 7.17]{Ha1} and Chow's lemma.
\vspace{-1.7ex}
\item Any morphism of toric varieties coming from subdivision of fans
are log modifications due to De Concini-Procesi's Theorem (cf.\ \cite[\S
      1.7]{TO1}).
\vspace{-1ex}
\end{itemize}

Note that log modifications are log \'etale.

\begin{lem}\label{lem-logmodificationcomposite}
\begin{enumerate}\vspace{-1ex}
\item {\it Log modifications are universally surjective. $($A morphism
      $X\rightarrow Y$ of fs log schemes is said to be universally
      surjective if for any
morphism $Y'\rightarrow Y$ of fs log schemes the underlying morphism of
the base change $X_{Y'}\rightarrow Y'$ is surjective.$)$}
\vspace{-1.7ex}
\item {\it If $X\rightarrow Y$ and $Y\rightarrow Z$ are log modifications,
then the composite $X\rightarrow Z$ is a log modification.}
\vspace{-1.7ex}
\item {\it If $\til{X}\rightarrow X$ is a log modification and $Y\rightarrow X$
is a morphism, then the base change $\til{X}_Y\rightarrow Y$ is a log
modification.}
\vspace{-1ex}
\end{enumerate}
\end{lem}

\noindent{\bf Proof.}\hspace{2mm}
First we prove (1).
It is enough to prove that log blow-up's are universally surjective.
By the universality of log blow-up's it suffices to show that any blow-up
is surjective.
By an easy observation we reduce to verify whether $\mathop{\mathrm{B}}\nolimits_{K}(P)\rightarrow
\mathop{\mathrm{S}}\nolimits(P)$ is surjective for a monoid $P$ and an ideal $K$.
But this is trivial since the underlying morphism of $\mathop{\mathrm{B}}\nolimits_{K}(P)\rightarrow
\mathop{\mathrm{S}}\nolimits(P)$ is a normalized blow-up in the usual sense.
To prove (2) we may assume there exists a log blow-up $\mathop{\mathrm{B}}\nolimits(Z)\rightarrow Z$
which factors through $Y$ by a log blow-up.
We can take a log blow-up $\mathop{\mathrm{B}}\nolimits(Y)\rightarrow Y$ which factors through
$X$ by a log blow-up, since $Y$ is easily seen to be quasi-compact.
Take the log blow-up $\mathop{\mathrm{B}}\nolimits(\mathop{\mathrm{B}}\nolimits(Z))\rightarrow\mathop{\mathrm{B}}\nolimits(Z)$ along the pull-back of
the blown-up center of $\mathop{\mathrm{B}}\nolimits(Y)\rightarrow Y$.
Then we have the unique morphism $\mathop{\mathrm{B}}\nolimits(\mathop{\mathrm{B}}\nolimits(Z))\rightarrow\mathop{\mathrm{B}}\nolimits(Y)$, and then
(2) is proved.
(3) is straightforward.
\hfill$\square$

\vspace{1ex}
\begin{thm}\label{thm-blowupbasechange3}
{\it Let $f\colon X\rightarrow Y$ be a morphism of locally noetherian fs log
schemes.
Suppose $\under{f}$ and $\under{Y}$ are quasi-compact.
Then there exists a commutative diagram}
$$
\begin{array}{ccc}
\til{X}&\longrightarrow&X\\
\llap{$\scriptstyle{\til{f}}$}\vphantom{\bigg|}\Big\downarrow&&\vphantom{\bigg|}\Big\downarrow
\rlap{$\scriptstyle{f}$}\\
\til{Y}&\longrightarrow&Y
\end{array}
$$
{\it of fs log schemes such that the horizontal arrows are log modifications
and $\til{f}$ is integral.}
\end{thm}

\noindent{\bf Proof.}\hspace{2mm}
By Theorem \ref{thm-blowupbasechange2} we may assume $f$ is exact, since
base changes by log blow-up's are log modifications.
Since we have assumed $\under{Y}$ is quasi-compact, we may assume $Y$ allows
a chart $\mu\colon Y\rightarrow\mathop{\mathrm{S}}\nolimits(Q)$ which is good at a point.
By Corollary \ref{cor-goodchart} we can take fppf locally on $X$ a chart
$(\la\colon X\rightarrow\mathop{\mathrm{S}}\nolimits(P),\mu\colon Y\rightarrow\mathop{\mathrm{S}}\nolimits(Q),h\colon
Q\rightarrow P)$ of $f$ which is good at a point.
Then $\mathop{\mathrm{S}}\nolimits(h)\colon\mathop{\mathrm{S}}\nolimits(P)\rightarrow\mathop{\mathrm{S}}\nolimits(Q)$ is a morphism of toric varieties
which comes from a morphism of fans $\Sigma'\rightarrow\Sigma$ (both of which
are defined by single cones) underlain by $\varphi\colon
N'\rightarrow N$.

For any ray $\rho\in\Sigma'_1$ the image of $\rho$ in $N$ is either a
ray or a point.
If it is a ray, then take the primitive base $n_1\in N$ of $\varphi(\rho)$,
and extend it to a $\Z$-basis $n_1,\ldots,n_r$ of $N$.
The $r+1$-rays spanned by $n_1,\ldots,n_r,-(n_1+\cdots+n_r)$ give
the projective space
$\mathbf{P}^{r}_{\Z}$, and the very ample divisor $\O(1)$ gives rise to
the support function, which we denote by $s_{\rho}$.
In case $\rho$ is mapped to the point we set $s_{\rho}=0$.
Then
$$
s=\sum_{\rho\in\Sigma'_1}s_{\rho}
$$
gives a piecewise linear function on $N\times_{\Z}\R$ which yields a
log-ideal $\mathcal{K}$ on $\mathop{\mathrm{S}}\nolimits(Q)$
(cf.\ \cite[I, Theorem 9]{KKM}).
Set $\mathcal{K}'=\mathop{\mathrm{S}}\nolimits(h)^{\bullet}\mathcal{K}$.
The induced morphism $\mathop{\mathrm{B}}\nolimits_{K'}(P)\rightarrow\mathop{\mathrm{B}}\nolimits_K(Q)$ is also an
equivariant morphism of toric varieties which comes from $\til{\Sigma}'
\rightarrow\til{\Sigma}$, where $\til{\Sigma}$ (resp.\ $\til{\Sigma}'$) is
the subdivision of $\Sigma$ (resp.\ $\Sigma'$) consisting of biggest cones
on which $s$ (resp.\ $\varphi_{\R}^{\ast}s$) is linear.
Since any cones in $\til{\Sigma}'$ is mapped surjectively onto some cone in
$\til{\Sigma}$, that the morphism of associated toric varieties
$\mathop{\mathrm{B}}\nolimits_{K'}(P)\rightarrow\mathop{\mathrm{B}}\nolimits_K(Q)$ has constant fiber dimension.
Taking more auxiliary log blow-up's of $\mathop{\mathrm{B}}\nolimits_K(Q)$ and $\mathop{\mathrm{B}}\nolimits_{K'}(P)$, we may
assume $\mathop{\mathrm{B}}\nolimits_K(Q)$ is non-singular.
In this case the morphism $\mathop{\mathrm{B}}\nolimits_{K'}(P)\rightarrow\mathop{\mathrm{B}}\nolimits_K(Q)$ is flat (since
toric varieties are Cohen-Macaulay), or what amounts to the same, integral
(as a morphism of fs log schemes) by \cite[(4,1)]{KK1}.

We gather up these resulting local log-ideals on $Y$ to make the log-ideal
$\mathcal{K}$.
We set $\mathcal{K}'=f^{\bullet}\mathcal{K}$ and then consider $\mathop{\mathrm{B}}\nolimits_{\mathcal{K}'}(X)
\rightarrow\mathop{\mathrm{B}}\nolimits_{\mathcal{K}}(Y)$.
Then by what we have seen above it is fppf locally integral, and hence
is integral.
\hfill$\square$
\begin{small}


\section{Introduction}

Effective field theories have found increasing interest
as powerful tools for describing the dynamics of physical systems
where global symmetries are spontaneously broken and continuous order
parameter fields represent the relevant low-energy degrees of freedom.
Depending on space dimensionality and the manifold on which the fields
live classical localized static solutions may fall into topologically
distinct classes characterized by integer winding numbers,
carrying energy, momentum and internal properties which suggest their
interpretation as particle-like excitations of the uniform ground
state. Their spatial size is determined by scaling
properties of different competing terms in the effective lagrangian.

Quantization of effective field theories allows to assign proper
quantum numbers to quasi-particle properties, identify
excited states of quasi-particles, and calculate loop
corrections to the classical results for observable quantities,
which may be crucial for experimental verification. The evaluation of
loop corrections necessarily brings about the need for renormalization.
Although, generally effective field theories might be
non--renormalizable in the strict sense, they still may be
renormalized order by order in terms of a gradient expansion.
A well-known example is chiral perturbation theory (ChPT), applied
successfully in hadron physics. In the vacuum sector such an expansion
is truncated by allowing only for external momenta small compared to the
underlying scale of the theory. However in the soliton sector,
the soliton itself constitutes "external" fields with gradients
comparable to the scale of the theory which cannot be made small
by definition. Therefore, in dealing with solitons, the problem
of truncating the expansion
can only be solved by the ad hoc assumption that the renormalized
couplings of the higher gradient terms are small.

A prominent example for the application of this program are 3D-$SU(N_f)$
skyr\-mions in $N_f$-flavor meson fields, where the topological charge
is identified with baryon number, with impressive success for baryonic
properties, resonances, and meson-baryon dynamics.
Similarly, the conjecture to consider 2D-$O(3)$ spin textures as charged
quasi-particles in ferromagnetic quantum Hall systems~\cite{SKK93}
and antiferromagnetic high--temperature superconductors~\cite{dpw88},
with topological winding density identified with the deviation of the
electron density from its uniform background value, suggests a
corresponding investigation.
For magnons and vortices with unity charge such attempts have
been presented for ferromagnets~\cite{a97} and
antiferromagnets~\cite{r89}. The field theoretical approach
forwarded here is rather
related to the antiferromagnet due to its preserved time--reversal
invariance which makes the analysis fairly similar to what applies
to relativistic systems. The ferromagnet where time--reversal
invariance is broken would require the consideration of the
Landau--Lifshitz dynamics.
The evaluation of the Casimir energies involves bound--state
energies and
a sum over scattering phase--shifts~\cite{dhn74}.
This provides also complete
information about resonant excited states in the continuum
which in the 3D-$SU(N_f)$ case successfully describe well-established
baryon resonances.

For the general outline of the necessary steps and technique we discuss
the case where
the second-order exchange energy is taken in the time--reversal
invariant form of the non--linear sigma model as it naturally appears in
relativistic field theories. In two dimensions it is scale invariant
and therefore irrelevant for the spatial extent of the static
localized solution. To fix the soliton size two more terms are needed.
We use the standard fourth-order Skyrme term, and a symmetry-breaking
coupling of the $O(3)$-order parameter to an anisotropy field.
The static part of the Skyrme term represents a local
approximation to the Coulomb energy of the charged excitations; in our
field-theoretic example we keep also the time-derivative parts of it.
So our discussion will mainly serve a demonstrative purpose as a model
field theory.

We shall specifically address the question of the binding energy of
doubly charged quasi particles.
For 3D-$SU(N_f)$ skyrmions this question has a long history since it was
discovered that in the winding number $B=2$ sector the static
lowest-energy solution is bound with respect to decay into two separate
$B=1$ skyrmions and displays only axial symmetry in contrast to the
radially symmetric hedgehog skyrmions. There has been much discussion
about the physical relevance (in the deuteron) of such a torus
configuration,
and it was speculated that quantum corrections might reverse the sign
of the binding energy $E(B=2)-2E(B=1)$. Evaluation of loop
corrections to the 3D-$SU(2)$ $B=2$ torus is a formidable task which to our
know\-ledge has not yet been achieved. It is interesting that the same
situation occurs for the 2D-$O(3)$ skyrmions: the classical $n=2$
solution shows a ring-like density distribution and it is bound with
respect to decay into two individual $B=1$ skyrmions. In this case,
however, the evaluation of loop corrections for both configurations is
of comparable complexity, and we will show that they in fact reverse
the sign of the binding energy.

Actually, there is an ongoing discussion~\cite{LLKS97,NK97}
about the binding of 2-skyrmions also in condensed matter applications.
However, it should be stressed that for ferromagnetic
systems the time-dependent part of the effective lagrangian has to be
replaced by the T-violating Landau-Lifshitz form with only one time
derivative.
For antiferromagnets
the time derivative coupling of the staggered
spin to the magnetic field, which does not contribute to the static
stabilization should be included~\cite{b98,h98,rs99}.
And, in any case, for quantitative
conclusions, the non--local character of the Coulomb energy should be
respected~\cite{mm97}.

\section{Static solutions}

The Lagrangian of the $O(3)$ model in $2+1$ dimensions
is conveniently written in terms of the $3$-component field
$\mbox{\boldmath $\Phi$}$ satisfying the constraint $\mbox{\boldmath $\Phi$} \cdot \mbox{\boldmath $\Phi$} =1$,
\begin{equation}\label{lag}
{\cal{L}} = \frac{f^2}{2} \partial_\mu \mbox{\boldmath $\Phi$} \partial^\mu \mbox{\boldmath $\Phi$}
- \frac{1}{4e^2} ( \partial_\mu \mbox{\boldmath $\Phi$} \times \partial_\nu \mbox{\boldmath $\Phi$} )^2
- f^2 m^2 (1-\Phi_3)    \, .
\end{equation}
The three terms represent the non--linear sigma
(N$\ell \sigma$) model, the Skyrme and the
potential term. There exists a conserved topological current
\begin{equation}\label{top}
T^\mu = \frac{1}{8 \pi} \epsilon^{\mu \nu \rho} \mbox{\boldmath $\Phi$}
( \partial_\nu \mbox{\boldmath $\Phi$} \times \partial_\rho \mbox{\boldmath $\Phi$} ) \, , \qquad
\partial_\mu T^\mu = 0
\, .
\end{equation}
The radially symmetric hedgehog ansatz, written
in polar coordinates $(r,\varphi)$
\begin{equation}\label{hedgehog}
\mbox{\boldmath $\Phi$} = \left( \begin{array}{c}
\sin F(r) \cos n \varphi   \\
\sin F(r) \sin n \varphi     \\
\cos F(r)
\end{array} \right) \, ,
\end{equation}
corresponds to winding number
\begin{equation}\label{winding}
\int d^2 r T^0 = \frac{n}{2}[\cos F(\infty) - \cos F(0)] = n
\, , \qquad
T^0 = - \frac{n F^\prime \sin F}{4 \pi r}
\end{equation}
and the soliton's size may be defined according to
\begin{equation}\label{radius}
<r^2>_n  \, = \frac{1}{n} \int d^2 r \, r^2 T^0 =
-\frac{1}{4\pi} \int d^2 r \, r F^\prime \sin F
\, .
\end{equation}
It is convenient to absorb the length $1/\sqrt{fem}$
which sets the scale for the size of localized structures
into the dimensionless spatial coordinate $x=\sqrt{fem} \, r$.
The static energy functional connected with the lagrangian
(\ref{lag}) is then given by
\begin{equation}\label{static}
E^{class}_n  =  \frac{f^2}{2} \int d^2 x \left[ F^{\prime 2}
+ \frac{n^2 s^2}{x^2} + a \frac{n^2 F^{\prime 2} s^2 }{x^2}
+ 2 a (1-c) \right]
 =  4\pi f^2 E_n^0 (a) \, ,
\end{equation}
with the abbreviations $s=\sin F$ and $c= \cos F$.
Technically, $a=m/fe$ is the only nontrivial parameter,
while $4 \pi f^2$ sets the overall energy scale. We shall therefore present
results for different values of $a$. The limit
$a \to 0$ describes the pure Belavin--Polyakov (BP) solution~\cite{bp75}.
The opposite limit, $a \to \infty$, is technically also well
defined; it describes a system without spin-spin aligning force.

The Euler--Lagrange equation following from the variation
of the energy functional
\begin{equation}\label{stability}
\frac{1}{x} \left(x F^\prime \right)^\prime - \frac{n^2 sc}{x^2}
+ a \frac{n^2 s}{x} \left( \frac{F^\prime s}{x} \right)^\prime
- a s = 0
\end{equation}
is a second order non--linear differential equation which is solved
subject to the boundary conditions
\begin{equation}\label{boundary}
F (0) = \pi \, , \qquad  F (x \to \infty) = 0
\, .
\end{equation}
In principle the boundary condition at the origin
according to (\ref{winding}) could be
$F(0) = \nu \pi$ with $\nu=1,3,\dots$ an odd multiple of $\pi$,
however it turns out that the solution with $\nu=1$ has the
smallest static energy.
\begin{figure}[htbp]
\begin{center}
\begin{tabular}{ll}
\epsfig{figure=baby1a.ps,width=5cm,angle=270}
&
\epsfig{figure=baby1b.ps,width=5cm,angle=270}
\end{tabular}
\protect\caption{Hedgehog profiles and energy densities for the
skyrmions with $n=1$ and $n=2$. The parameter $a^2=0.1$ was
used.
}
\end{center}
\end{figure}

The hedgehog profiles and the corresponding energy densities for
$a^2=0.1$ are shown in Fig. 1 for the 1- and 2-skyrmions.
It is noticed that the maximum energy density for the $n=2$ soliton
is not located at the
origin but in a spherical shell of radius $x \simeq 1$.

The classical energies and square radii are given in Table \ref{tab_1}.
\begin{table}
\begin{center}
\parbox{10cm}{\caption{\label{tab_1}
Classical soliton energies and square radii
for the stable 1- and 2-skyrmions
for various parameters $a$.
}}
\begin{tabular}{|c|cccccc|}
\hline
$a$    & $0.0$  &  $0.01$  &  $0.0316$  & $0.1$ & $0.316$ & $1.0$ \\
\hline
$E_1^0$     & $1.00$ & $1.03$  & $1.08$   & $1.21$ & $1.56$ & $2.57$   \\
$E_2^0$     & $2.00$ & $2.03$  & $2.10$   & $2.30$ & $2.94$ & $4.85$   \\
\hline
$<x^2>_1$   & $\infty$ & $2.44$ & $2.11$  & $1.80$   & $1.57$ & $1.42$   \\
$<x^2>_2$   & $ \pi  $ & $3.11$ & $3.06$  & $2.99$   & $2.78$ & $2.76$   \\
\hline
\end{tabular}
\end{center}
\end{table}
In the two limiting cases, $a \to 0$ and $a \to \infty$, the
differential equation (\ref{stability}) may be solved
analytically (see appendix). The corresponding static energies
and radii are
\begin{eqnarray}\label{limit}
E_n^0=n \, , \qquad \quad & \, \qquad
<x^2>_n=\sqrt{\frac{2n}{3} (n^2-1)}
\cdot \frac{\pi/n}{\sin (\pi/n)}
\qquad \qquad & a<<1 \nonumber \\
E_n^0=4n/3 \cdot a \, , \, & <x^2>_n=4n/3
\qquad \qquad \qquad \qquad  & a>>1
\, .
\end{eqnarray}
It is noticed that for the BP solution with $a \to 0$
the 1-skyrmion's radius diverges which is due to the particular
choice of the potential term in (\ref{lag}). This divergence
is perceptible only for extremely small values of $a$, for
$a \stackrel{>}{_\sim} 0.001$ the radius for the 1-soliton is
still smaller than that of the 2-soliton which stays finite
in the limit $a \to 0$ (for a more thorough discussion of this
problem we refer to the appendix).
The dependence of the energies on the dimensionless parameter
$a$ is shown in Fig. 2. Independently from this parameter
$E_2^0 \le 2E_1^0$ always holds, the 2-skyrmion is classically
stable against decay into two 1-skyrmions.
\begin{figure}[htbp]
\centerline{\epsfig{figure=baby2.ps,width=7cm,angle=270}}
\protect\caption{
Classical soliton energies for $n=1$ and $n=2$ as functions
of the dimensionless parameter $a$. For better comparison
the energies plotted are divided by the corresponding topological
charge. The 2-skyrmion is classically stable against decay
into two 1-skyrmions.
}
\end{figure}
Hedgehog solutions with higher topological
charge, $n \ge 3$, are not stable (see subsection 3.3
for their stability)
and do not represent the minimum energy configuration of the
corresponding sector. The 3-soliton's minimum energy configuration
is not rotationally symmetric, it is a linear molecule which
consists of three distorted 1-solitons. Higher multi--solitons
are then supposedly molecules of stable 2-skyrmions (and
1-skyrmions in case of $n$ odd)~\cite{psz95}.  This
picture may change fundamentally when the Casimir energies
are taken into account.

\section{Small amplitude fluctuations}

The evaluation of quantum corrections is based on the normal modes
of the classical configurations. With the appropriate boundary
conditions these normal modes describe excited localized (bound) states,
and the scattering of the vacuum fluctuations (mesons or magnons)
off the soliton backgound. (The S-matrix for magnon-vortex scattering
has been considered (without the stabilizing Skyrme and potential
terms) in~\cite{ikw96} for antiferromagnets and
in~\cite{ismw98} for ferromagnets).
Because of the constraint $\mbox{\boldmath $\Phi$} \cdot \mbox{\boldmath $\Phi$} =1$ the model
effectively possesses two independent fields.
We therefore introduce two--component small amplitude fluctuations
\begin{equation}\label{fluct}
\mbox{\boldmath $\eta$} = \mbox{\boldmath $\hat r$}_n \eta_L
+ \mbox{\boldmath $\hat \varphi$}_n \eta_T \, , \qquad
\mbox{\boldmath $\hat r$}_n = \left( \begin{array}{c}
\cos n \varphi   \\ \sin n \varphi  \end{array} \right) \, , \quad
\mbox{\boldmath $\hat \varphi$}_n = \left( \begin{array}{c}
-\sin n \varphi   \\ \cos n \varphi  \end{array} \right)
\end{equation}
which we decompose into components parallel and perpendicular
to the soliton configuration with winding number $n$. Of course,
there are many different ways to parametrize the total
time-dependent field $\mbox{\boldmath $\Phi$}$ which lead to different equations of
motion for the corresponding fluctuations. However,
the bound--states and the scattering matrix,
in particular the phase--shifts are independent of the chosen
parametrization. A very convenient choice
different from the common Polyakov parametrization~\cite{p75},
which only at first sight may seem complicated, is given by
\begin{equation}\label{ck}
\mbox{\boldmath $\Phi$} = \left( \begin{array}{c}
\mbox{\boldmath $\hat r$}_n \sin F (1-\mbox{\boldmath $\eta$}^2/2)
+ \mbox{\boldmath $\hat r$}_n \cos F \eta_L
+ \mbox{\boldmath $\hat \varphi$}_n \eta_T  \\
\cos F (1-\mbox{\boldmath $\eta$}^2/2) - \sin F \eta_L
\end{array} \right)
\stackrel{r \longrightarrow \infty}{\longrightarrow}
\left( \begin{array}{c} \mbox{\boldmath $\eta$} \\
1-\mbox{\boldmath $\eta$}^2/2 \end{array} \right)  \, .
\end{equation}
The main advantage of this parametrization is that it leads to a flat
metric for the non--linear sigma model part of the lagrangian.
For 3D-$SU(N)$ skyrmions this parametrization is well--known
as Callan--Klebanov ansatz~\cite{ck85}.

It is now straightforward to write down the e.o.m. for the
fluctuating components $\eta_L$ and $\eta_T$. The Lagrangian
(\ref{lag}) has to be expanded to quadratic order
in the fluctuations, then the e.o.m. can be read off. We give these
coupled linear differential equations explicitely, where we have
already exploited the time--dependence as well as the angular
dependence of the fluctuations
\begin{equation}\label{phonon}
\eta_L = \sum_M f_M (x) e^{iM\varphi}  e^{-i \omega t} \, , \qquad
\eta_T = i \sum_M g_M (x) e^{iM\varphi}  e^{-i \omega t}
\, .
\end{equation}
The e.o.m. decouple in the magnetic quantum number $M$, which is the
analogon of the phonon- (or grand-) spin familiar from scattering
calculations for 3D-$SU(N)$ skyrmions
\begin{eqnarray}\label{eom}
&&-\frac{1}{x} \left( x b_L f_M^\prime \right)^\prime
+ \frac{M^2 }{x^2} f_M + \frac{n^2}{x^2} ( b_T-2s^2 ) f_M
-a \frac{2n^2c}{x} \left( \frac{F^\prime s}{x} \right)^\prime f_M
+ a c f_M \qquad
\nonumber \\
&& \qquad \qquad -\frac{nMc}{x^2} ( 1+b_T ) g_M
+ a \frac{nM F^\prime s}{x^2} g_M^\prime
+ a \frac{2nM}{x} \left( \frac{F^\prime s}{x} \right)^\prime g_M
= \omega^2 b_L f_M  \nonumber \\
&& -\frac{1}{x} \left( x g_M^\prime \right)^\prime
+ \frac{M^2 }{x^2} b_T g_M
- ( F^{\prime 2} - \frac{n^2 c^2}{x^2} ) g_M
-a \frac{n^2c}{x} \left( \frac{F^\prime s}{x} \right)^\prime g_M
+ a c g_M  \qquad \\
&& \qquad \qquad -\frac{nMc}{x^2} ( 1+b_T) f_M
- a \frac{nM F^\prime s}{x^2} f_M^\prime
+ a \frac{nM}{x} \left( \frac{F^\prime s}{x} \right)^\prime f_M
= \omega^2 b_T g_M \nonumber
\, .
\end{eqnarray}
Here we introduced the metric functions $b_L(x)=1+a n^2 s^2/x^2$
and $b_T(x)=1+a F^{\prime 2}$. The energies $\omega$ are
understood in natural units of $\sqrt{fem}$ such that the threshold occurs
at $\omega = \sqrt{a}$.
The equations are invariant with
respect to the simultaneous replacements $M \to -M$ and $g_M \to -g_M$.
These equations contain
all the information about zero--modes, bound--states, the scattering
matrix and the stability of the hedgehog solutions.

\subsection{Zero-modes}

The soliton's energy (\ref{static}) is invariant under spatial
rotations around the $z$-axis (which for the hedgehog
is equivalent to an iso--rotation around the internal $3$-axis)
as well as under a translation in the
$x$-$y$ plane. The infinitesimal rotation (iso--rotation)
and translation give rise to zero--modes which are solutions
of the e.o.m. (\ref{eom}) for $\omega^2=0$.

The rotational zero--mode is obtained by a shift of the azimuthal
angle $\varphi \to \varphi + \alpha$ in the hedgehog ansatz
(\ref{hedgehog}). Comparison with (\ref{ck}) and (\ref{phonon})
shows then that this zero--mode
\begin{equation}\label{rot}
f_0(x)=0 \, , \quad g_0(x)=s \qquad \mbox{(rotational zero-mode)}
\end{equation}
is located in the $M=0$ partial--wave.
Similarly, the translational zero--mode
$\mbox{\boldmath $r$} \to \mbox{\boldmath $r$} + \mbox{\boldmath $a$}$
\begin{equation}\label{trans}
f_1(x)=F^\prime \, , \quad g_1(x)=\frac{ns}{x}
\qquad \mbox{(translational zero-mode)}
\end{equation}
sits in the $M=1$ partial--wave.

With the stability condition (\ref{stability})
it is easily checked that the above modes satisfy the e.o.m. for
$\omega^2=0$. As both zero--modes possess a finite norm, they will act
as bound--states in the scattering calculation and influence the
phase--shifts via Levinson's theorem. Besides these
zero energy bound--states there are also bound--states at finite
energy.

\subsection{Bound-states}

In both systems, $n=1$ and $n=2$, we observe bound--states
in the $M=n$ partial--wave. With decreasing parameter $a$ the
energies of these bound--states move towards the threshold,
in particular for the $n=1$ system this happens very quickly
(cf. Table \ref{tab_2}).
There is a simple explanation of this phenomenon. As is noticed
from the e.o.m. the iso--rotation with respect to the $1$- and
$2$-axes
\begin{equation}\label{bound}
f_n(x)=1 \, , \quad g_n(x)=c
\qquad (\mbox{zero-energy solution for} \quad a \to 0)
\end{equation}
is a zero--energy solution of the $M=n$ partial--wave for $a \to 0$.
The potential connected with a finite $a \not= 0$ distorts this
continuum solution into a bound--state, which in the limit $a \to 0$
is shifted towards the threshold.
\begin{table}[htbp]
\begin{center}
\parbox{10.5 cm}{\caption{\label{tab_2}
Bound-state energies relative to the threshold
in the $M=n$ partial--wave. With decreasing
parameter $a$ the bound--states move towards threshold.
}}
\begin{tabular}{|c|cccccc|}
\hline
$a$   & $1.0$   & $0.707$   & $0.316$       & $0.224$      &
                               $0.1$         & $0.071$      \\
\hline
$n=1$  & $0.943$ & $0.99994$ &  $\simeq 1.0$ & $\simeq 1.0$ &
                                $\simeq 1.0$ & $\simeq 1.0$ \\
$n=2$  & $0.330$ & $0.394$   &  $0.622$      & $0.690$      &
                                $0.866$      & $0.943$      \\
\hline
\end{tabular}
\end{center}
\end{table}

In the $n=1$ system this is the only bound--state apart from the
zero--modes already discussed. However in the $n=2$ case some of the
low-lying resonances may also become bound for sufficiently strong $a$,
e.g. for $a=1$ we find additional bound--states in the $M=0$ and
$M=1$ partial--waves at energies $\omega_0=0.923$ and
$\omega_1=0.793$ respectively. All these bound--states must be
considered for the Casimir energy (compare also the discussion
of the phase--shifts in subsection 3.4).

\subsection{Stability of hedgehog solitons}

Formally, eqs. (\ref{eom}) are the e.o.m. for small amplitude fluctuations
around a hedgehog soliton with arbitrary topological quantum number
$n$. For this reason they contain information about the stability of
the $n$-skyrmion. For the $n$-skyrmion to be stable there must not
exist solutions to (\ref{eom}) with negative $\omega^2$.

Indeed, we do find solutions of (\ref{eom}) with negative $\omega^2$
for topological quantum numbers $n \ge 3$. For example for $n=3$
with $a^2=0.1$ there exist such solutions in the
$M=2$ ($\omega^2=-0.14$), $M=3$ ($\omega^2=-0.30$) and
all higher partial waves. In the $n=1$ and $n=2$ sectors
there exist no negative $\omega^2$ solutions.

We conclude, that the $n=1$ and $n=2$ hedgehog solitons are stable
while those with higher topological charges are unstable. This
agrees with the findings in ref.~\cite{psz95}. In what follows we
consider the scattering off the stable skyrmions with $n=1$ and
$n=2$.

\subsection{Phase--shifts}

For the $2$-dimensional scattering problem in spherical
coordinates~\cite{mf} the incident plane wave
is decomposed into partial waves
carrying magnetic quantum numbers
\begin{equation}\label{plane}
e^{i \mbox{\boldmath $pr$}} = \sum_{M=-\infty}^{\infty} i^M
J_M(pr) e^{iM \varphi} =
\sum_{M=0}^{\infty} \epsilon_M i^M J_M(pr) \cos (M \varphi) \, .
\end{equation}
The latter transformation with the multiplicities $\epsilon_0=1$
and $\epsilon_M=2$ for $M \ge 1$
follows from the properties of the regular
and irregular Bessel functions of the first kind, $J_M(pr)$ and
$N_M(pr)$, when $M \to -M$. Thus, it suffices to consider the
partial waves $M \ge 0$. Similarly, for the scattered wave
(plane wave plus outgoing radial wave) in the case of a single
channel we have
\begin{eqnarray}\label{scatt}
\Psi &=& \sum_{M=0}^{\infty} \epsilon_M i^M \psi_M (p,r)
\cos (M \varphi)
\, .
\end{eqnarray}
The partial--wave projected scattering waves $\psi_M (p,r)$ with
appropriate boundary
conditions at the origin are integrated according to the e.o.m. and
the phase--shifts then are obtained from the asymptotic form
\begin{equation}\label{scatt1}
\psi_M (p,r) \to [ J_M(pr) \cos \delta_M(p) + N_M(pr) \sin \delta_M(p) ]
e^{i \delta_M(p)}
\, .
\end{equation}
The phase--shifts are related to the cross--section in the usual way,
and the property $\delta_M(p)=\delta_{-M}(p)$
follows from the corresponding symmetry of the e.o.m.. The
generalization to two coupled channels now is straightforward.
It is noticed that our scattering equations (\ref{eom})
decouple asymptotically when the functions
$f_{\pm}= (f_M \mp g_M) / \sqrt{2} $ are introduced
\begin{equation}\label{asym}
-\frac{1}{r} \left(r f_\pm^\prime \right)^\prime
+ \frac{(M \pm n)^2}{r^2} f_\pm = p^2 f_\pm \, , \qquad
\omega^2=p^2+m^2 \, ,
\end{equation}
with the solution
\begin{equation}\label{asym1}
f_{\pm} (p,r) = A_{\pm} (p) J_{|M \pm n|}(pr)
+ B_{\pm} (p) N_{|M \pm n|}(pr)
\, .
\end{equation}
From the coefficients $A_{\pm} (p)$ and $B_{\pm} (p)$ of the
asymptotical solution, the $2 \times 2$ scattering matrix $S_M$
is obtained for every partial wave $M$. Again, as the e.o.m. are
invariant with respect to $M \to -M$ so is the scattering matrix
and it suffices to consider $M \ge 0$. It is convenient to
diagonalize the scattering matrix
\begin{equation}\label{smatrix}
S_M = U_M \left( \begin{array}{cc}
e^{i\delta_M^{(1)}}     &    0   \\
0        &    e^{i\delta_M^{(2)}}      \end{array} \right)
(U_M)^\dagger \, ,
\end{equation}
in order to obtain the intrinsic eigen phase--shifts $\delta_M^{(1)}$
and $\delta_M^{(2)}$. Their sum $\delta_M = \delta_M^{(1)} +
\delta_M^{(2)}$ is plotted for the partial waves with $M \le 4$
in Fig. 3, where the $n=1$ and $n=2$ systems are considered separately.

The picture that emerges looks qualitatively quite similar to what
has been obtained for 3D-$SU(N)$ skyrmions long ago~\cite{we84}.
For topological charge $n=1$ we obtain in detail:
The phases for $M=0$ and $M=1$, where the rotational and translational
zero modes are located together with the $M=1$ bound--state,
start at $\pi$ resp. $2\pi$ according to Levinson's theorem.
For the smaller parameters $a^2=0.01$ and $0.1$ the
wave--function of the $M=1$
bound--state close to threshold (Table \ref{tab_2})
is already quite similar to the infinitesimal translation
(cf. subsection 3.5 for the limit $a \to 0$).
As a consequence the phase--shift at threshold
reacts with a sudden drop from $2\pi$ to $\pi$ which in the
figures is only noticed because the phase--shift falls below
$\pi$ before it bends to the right. A weakly pronounced breathing
mode is observed at low energies for $M=0$. The
partial--waves $M=2,3,\dots$ then contain a band of resonances
which are flattened out for higher $M$ and also with decreasing
parameter $a$.

Similarly for topological charge $n=2$: For the smaller parameters
$a^2=0.01$ and $0.1$ the $M=0$ and $M=1$ phase--shifts start
at $\pi$ because of the zero modes. In all but the $M=2$
partial--wave with its bound--state (cf. Table \ref{tab_2})
sharp resonances occur, which are much more pronounced compared
to the $n=1$ case. For strong potentials then, the $M=0$ and $M=1$
resonances become bound (e.g. for $a=1$ at $\omega_0=0.923$
and $\omega_1=0.793$). Also less pronounced secondary resonances
do appear.
\begin{figure}
\begin{center}
\begin{tabular}{cc}
\epsfig{figure=baby3a.ps,width=5cm,angle=270} &
\epsfig{figure=baby3b.ps,width=5cm,angle=270} \\
\epsfig{figure=baby3c.ps,width=5cm,angle=270} &
\epsfig{figure=baby3d.ps,width=5cm,angle=270} \\
\epsfig{figure=baby3e.ps,width=5cm,angle=270} &
\epsfig{figure=baby3f.ps,width=5cm,angle=270}
\end{tabular}
\protect\caption{Phase shifts for the lowest partial waves
scattered off the stable $n=1$ and $n=2$ skyrmions.
The relevant parameter was fixed at $a^2=0.01$, $0.1$ and $1.0$.
The momentum $p$ is plotted in natural units.
According to Levinson's theorem the phase--shifts start at
multiples of $\pi$. Many resonances are observed in the various
channels.
}
\end{center}
\end{figure}

\subsection{Belavin-Polyakov soliton}

Finally we would like to add a few comments on the pure
Belavin--Polyakov solution~\cite{bp75}.
As mentioned we recover this solution
for small parameters $a$ with its size $R$ still fixed by the
balance of the Skyrme and potential terms (cf. appendix). In the
pure N$\ell \sigma$ model without additional stabilizing terms
the size $R$ of the BP soliton is an undetermined parameter.
In contrast to the full lagrangian (\ref{lag}) the pure
N$\ell \sigma$ model does not break $O(3)$--symmetry and is also
conformally invariant. Due to these additional symmetries a larger
number of zero--modes is expected.

The e.o.m. for the fluctuations obtained in that case from (\ref{eom})
with $a=0$ and $x=r/R$ decouple for $f_{\pm}= (f_M \mp g_M) / \sqrt{2}$
(not only asymptotically) and are easily solved. Still the counting of the
zero--modes is a bit tricky. In addition to the rotation
(\ref{rot}), the breathing--mode (conformal invariance) becomes a
zero--mode in the $M=0$ partial--wave as expected. But for the
1--soliton the iso--rotation (\ref{bound}) coincides with the
translation (\ref{trans}) giving rise to only $2 \times 1$
zero--modes in the $M=1$ partial--wave. This makes altogether
$4$ zero--modes for the 1--soliton. Similar counting for the
2--soliton yields $8$ zero--modes.

Of course, this discussion is somewhat academic, because in reality
there should always be some stabilizing term which fixes the size
of the soliton and which at least breaks conformal invariance.
Therefore we do not show the corresponding phase--shifts.

\section{Casimir energy}

With the bound--state energies and the phase--shifts provided in
the previous section we are in the position to evaluate the
1-loop contribution to the soliton energy, i.e. the Casimir energy.
The UV singularities contained in the loop are related to the
high momentum behaviour of the phase--shifts
\begin{equation}\label{psum}
\delta (p) = \sum_{M=0}^\infty \epsilon_M \delta_M (p)
\stackrel{p \longrightarrow \infty}{\longrightarrow}
a_0 p^2 + a_1
\, .
\end{equation}
In the case of a vanishing Skyrme term the coefficients
\begin{equation}\label{coeff}
a_0=0 \, , \qquad a_1=\frac{1}{4} \int d^2 r \left[ F^{\prime 2}
+ \frac{n^2 s^2}{r^2} + 2 m^2 (1-c) \right]
\end{equation}
are known analytically from the Born terms.
For the full model they have to
be extracted numerically from the phase--shift sum (\ref{psum}),
which for that purpose has to be calculated to a high precision
with some $100$ partial waves taken into account.
Typical values obtained, e.g. for $a^2=0.1$, are $a_0=1.0 (fem)^{-1}$
and $a_1=7.2$ for the 1-soliton, and $a_0=2.0 (fem)^{-1}$
and $a_1=13.4$ for the 2-soliton respectively.
According to (\ref{coeff})
the N$\ell \sigma$ plus potential term contributions to these
coefficients are $a_1=8.3$ for $n=1$ and $a_1=15.6$ for $n=2$.
With this established we may use the phase--shift formula~\cite{dhn74},
subtract the troublesome high momentum behaviour from the phase--shifts,
and separately add the corresponding counterterms
\begin{eqnarray}\label{1-loop}
E^{1-loop} \! &=& \! \frac{1}{2\pi}
\left[ -\int_0^\infty d p \frac{p}{\omega}
\delta (p) - m \delta (0) \right] + \frac{1}{2} \sum_B \omega_B
\\
\! &=& \! -\frac{1}{2\pi} \int_0^\infty d p \frac{p}{\omega}
\left( \delta (p) - a_0 p^2 - a_1 \right)
+ \frac{1}{2} \sum_B (\omega_B - m)
- \int \frac{d^2p}{(2\pi)^2} \, \frac{a_0p^2+a_1}{\omega}
\, . \nonumber
\end{eqnarray}
This
procedure is closely related to what has been employed in the $3+1$
dimensional Skyrme model for the calculation of the Casimir energy to the
nucleon mass~\cite{m93,mw97}.
The counter terms could e.g. be evaluated using a 2-momentum cutoff
$\Lambda >> m$
\begin{equation}\label{cut}
\int \frac{d^2p}{(2\pi)^2} \, \frac{1}{\omega} =
\frac{1}{2\pi} [\Lambda -m] \, , \qquad
\int \frac{d^2p}{(2\pi)^2} \, \frac{p^2}{\omega} =
\frac{1}{6\pi} [\Lambda^3 - \frac{3}{2} m^2 \Lambda + 2m^3] \, ,
\end{equation}
which makes the linear and cubic divergencies explicit
\begin{equation}\label{loop}
E^{1-loop} =
-\frac{a_0}{6\pi} \left[ \Lambda^3 - \frac{3}{2} m^2 \Lambda \right]
- \frac{a_1}{2\pi} \Lambda + E^{cas}
\, .
\end{equation}
While the N$\ell \sigma$ and the potential term contributions
located in $a_1$ may be absorbed in a renormalized constant
\begin{equation}\label{fren}
f^2 = \bar f^2 - \frac{\Lambda}{4\pi}
\end{equation}
($\bar f$ denotes the bare constant), the divergencies stemming
from the Skyrme term renormalize the Skyrme parameter $e$ as well
as the couplings of all the higher gradient terms not listed in
the lagrangian (\ref{lag}). The latter renormalization we do not
carry out explicitely, instead we simply assume that the renormalized
Skyrme parameter be $e$ and the renormalized couplings of the
higher gradient terms be zero. To which extent this is a consistent
assumption will be discussed in the following section in connection
with the scale--dependence of our results. Now all the divergencies
residing in (\ref{loop}) are taken care of, leaving the finite
Casimir energy
\begin{equation}\label{casimir}
E^{cas}= -\frac{1}{2\pi} \int_0^\infty d p \frac{p}{\omega}
\left( \delta (p) - a_0 p^2 - a_1 \right)
+ \frac{1}{2} \sum_B (\omega_B - m)
- \frac{a_0 m^3}{3\pi} + \frac{a_1 m}{2\pi}
\, .\end{equation}
This result is independent of the employed regularization scheme.
For instance, using dimensional regularization the counter
terms (\ref{cut}) in $2+1$ dimensions become finite and equal
to the last terms in the brackets, respectively,
which yields the same expression for the Casimir energy.
However, although dimensional regularization in odd
dimensions has been used right from the beginning~\cite{bg72} and
is now commonly applied in $\Phi^4_3$, QCD$_3$ and
Chern-Simons theories~\cite{bg86,m90},
we hesitate to apply this scheme for an odd number of dimensions
because the fate of the UV singularities remains obscure.
As the divergencies do not show up in odd dimensional
regularization there seems to be no possibility to introduce
a renormalization scale. For strictly (super) renormalizable theories
this may be of no further importance, however in theories which are
renormalized order by order in terms of a gradient expansion a scale
is desirable in order to control the convergence of the series
(cf. section 5).

\begin{figure}[htbp]
\centerline{\epsfig{figure=baby4.ps,width=7cm,angle=270}}
\protect\caption{Subtracted phase--shift sum
as a function of the momentum $p$ in natural units
for the $n=1$ and $n=2$ systems calculated with $a^2=0.1$.
}
\end{figure}
The Casimir energy (\ref{casimir}) consists
of three parts, i.e. the phase--shift integral, the
bound--state contribution and the contribution from the
counter--terms.

The subtracted phase--shift sum $\delta (p) - a_0 p^2 - a_1$ which
enters the phase--shift integral is plotted in Fig. 4 for the
$n=1$ and $n=2$ systems ($a^2=0.1$). The maximum momentum to
which this sum has to be integrated depends sensitively on the
parameter $a$.
The various contributions to the Casimir energy are given in
Table \ref{tab_3} for several values of $a$.
\begin{table}[htbp]
\begin{center}
\parbox{12cm}{\caption{\label{tab_3}
Individual contributions to the Casimir energy according
to (\ref{casimir}) for the $n=1$ and $n=2$ systems. All
energies are given in natural units.
}}
\begin{tabular}{|c|rrr|rrr|}
\hline
   &&  $n=1$ &&& $n=2$ & \\
\hline
$a$  & $0.1$   &  $0.316$  &  $1.0$ & $0.1$   &  $0.316$  &  $1.0$ \\
\hline
phase-shifts
   & $-0.29$  & $-0.08$ & $0.00$ & $-0.63$  & $-0.26$ & $0.05$  \\
bound-states
   & $-0.47$ & $-0.84$ & $-1.56$ & $0.52$ & $-1.06$ & $-2.47$  \\
counter terms
   & $0.31$  & $0.62$  & $1.85$  & $0.61$  & $1.16$  & $3.36$  \\
\hline
total
   & $-0.45$  & $-0.30$ & $0.29$ & $-0.54$  & $-0.15$ & $1.00$  \\
\hline
\end{tabular}
\end{center}
\end{table}
The dependence of the Casimir energies
on this parameter is shown in Fig. 5. Note that for an easier
comparison again only half of the 2-skyrmion's Casimir energy is
plotted.
\begin{figure}[htbp]
\centerline{\epsfig{figure=baby5.ps,width=7cm,angle=270}}
\protect\caption{Casimir energies in natural units
for the 1- and 2-skyrmion
devided by their topological charge.
}
\end{figure}
It is noticed that the Casimir energy of the 1-skyrmion stays always
smaller compared to half of that of the 2-skyrmion. Therefore, in
contrast to the classical energy the Casimir energy favors single
skyrmions. These competing effects may be summarized in the formula
for the total energy cast into the form
\begin{eqnarray}\label{total}
E_n^{total}&=&4\pi f^2 E_n^0(a) +  \sqrt{fem}  E_n^1(a)
\nonumber \\
&=& 4\pi f^2 \left[ E_n^0(a) +   y \,  E_n^1(a) \right]
\, , \qquad \qquad  y=\frac{\sqrt{fem}}{4\pi f^2}
\, ,
\end{eqnarray}
where $E_n^1(a)$ represents the Casimir energy (\ref{casimir})
in natural units $\sqrt{fem}$.
Thus, apart from the overall scale $4\pi f^2$ the model is
characterized by the two dimensionless parameters $a$ and $y$.
If this ratio $y$ exceeds a critical value (which
increases slowly with $a$,
e.g. $y=0.3, 0.4, 0.7$ for $a^2=0.01, 0.1, 1.0$), then
$E_2^{total}-2E_1^{total}$ becomes positive and the 2-skyrmion
decays into two single 1--skyrmions.

\subsection{Numerical example from hadron physics}

Although different dimensionalities may cause qualitative
differences, let us illustrate these results by a numerical example
using the scale of hadron physics:
Take the value $a=0.1$, fix the size scale at
$1/\sqrt{fem}=0.2 fm$  ($\sqrt{fem} = 1$ GeV) and make use
of the numbers given in Table 1 and Table 3.
With the above parameters the 1-soliton gets a
topological radius $0.27$ fm and a Casimir energy
$E_1^{cas}=-0.45$GeV, a situation very similar to what has been
obtained in the hadronic Skyrme model~\cite{mw97}, if in addition
the classical soliton mass is fixed at $E_1^{class}=1.21 (4\pi f^2)
=1.39$ GeV in order to obtain the nucleon mass at
$E_1^{total}=0.94$ GeV. Then the $n=2$ soliton's total energy
$E_2^{total}=(2.64 - 0.53)$ GeV $= 2.11$ GeV
turns out to be larger than that of two separated 1-skyrmions.
The parameter in this example,
$y = 0.9$, lies far above the critical value $0.3$.
Of course, the $2+1$ dimensional model discussed
here cannot be taken literally for hadron
physics, but exactly this mechanism may shift the
undesired torus configuration obtained in hadronic soliton models
with baryon number $B=2$ to higher energies.

\subsection{Casimir energy of the pure Belavin-Polyakov soliton}

Finally we would like to discuss the Casimir energy of the pure
BP soliton with free size parameter $R$ in order
to make contact to the quantum corrections as obtained in
ref.~\cite{r89}. With $m=0$ and the coefficients $a_0=0 \, ,
\, a_1=2\pi n$ and no additional bound--states present (cf. the
discussion about the zero--modes in subsection 3.5) eq.(\ref{1-loop})
for the 1--loop contribution simplifies
\begin{eqnarray}\label{bp}
E^{1-loop} \! &\stackrel{BP}{=}& \!
- \frac{1}{2\pi} \int_0^{(\Lambda)} d p \delta (p)
\nonumber \\
\! &=& \! -\frac{1}{2\pi} \int_0^{(\Lambda)} d p
\left( \delta (p) - a_1 \right)   - \frac{a_1}{2\pi} \Lambda
\\
\! &=& \! -\frac{1}{2\pi R} \int_0^{(\Lambda R)} d (pR)
\left( \delta (pR) - 2 \pi n \right)   - n \Lambda
\, . \nonumber
\end{eqnarray}
We indicate here the possibility to limit the momentum integration
by a finite Debye momentum $\Lambda=p_D$ related to the lattice constant.
Because the soliton should at least cover several lattice
sites, $\Lambda R > 1$ is large enough to extend
the integration over the subtracted phase--shifts
to infinity (see Fig. 4) just as in our field theoretical treatment.
Through the substitution in the last step the dependence of the subtracted
phase--shift integral on the soliton size $R$ becomes explicit:
the magnitude of the Casimir energy decreases with increasing
soliton size.

Numerically, we find
\begin{equation}
E^{cas}_1=-\frac{0.5}{R_1} \, , \qquad \qquad
E^{cas}_2=-\frac{1.7}{R_2}
\end{equation}
for the pure BP solitons with topological charges $n=1$ and $n=2$
respectively. Their sizes $R_1$ and $R_2$ are arbitrary and generally
different. If the sizes are
determined by the stabilizing terms in the limit $a \to 0$,
we obtain $R_n=b_n/\sqrt{fem}$ with $b_1 \to 0$ and $b_2 \to \sqrt{2}$
(appendix). For that reason
the Casimir energy of the 1-soliton plotted in Fig. 5 bends to negative
values and finally diverges as $a$ approaches zero. The Casimir energy
of the 2-soliton stays finite,
$E^1_2 = -1.7/\sqrt{2} + {\cal O} ( \sqrt{a} )$, but with an infinite
slope at $a=0$.

In ref.~\cite{r89} the subtracted phase--shift integral was neglected for
$p_D R >> 1$ and the last
term in (\ref{bp}) then lead to the result
$E^{1-loop} \stackrel{BP}{=} -n p_D$. Because this is exactly the
term we absorb into a renormalized $f^2$ (\ref{loop},\ref{fren}),
this result corresponds
to $E^{cas}=0$ in our notation. The opposite limit $p_D R << 1$
where the 1-loop contribution (\ref{bp}) vanishes
corresponds to soliton sizes much smaller than the lattice spacing
and does not make sense. The conclusion, that the quantum corrections
lower the soliton energy as soliton size increases, is wrong.

\section{Scale dependence}

In this section we introduce a scale $\mu$ which allows to shift finite
pieces from the tree to the 1-loop contribution and vice versa.
It will be introduced in such a way
that we recover the results of the previous section for
$\mu = \mu_0$. Tentatively, $\mu_0=4\pi f^2$ may be identified with
the underlying energy scale of the model in
analogy to the chiral scale $\mu_0=4\pi f_\pi \simeq 1$GeV of
ChPT~\cite{e96}. Our results will however be presented in such a way
that no fixation of $\mu_0$ is required.
For this purpose we write the cutoff characterizing the divergencies
of the 1-loop contribution (\ref{loop}) as a sum of two parts
$\Lambda = [\Lambda - (\mu-\mu_0)] +(\mu-\mu_0)$. The first part
is then renormalized into a scale-dependent strength of the
N$\ell \sigma$ and potential term
\begin{equation}\label{fren1}
f^2 (\mu) = \bar f^2 - \frac{\Lambda - (\mu-\mu_0)}{4\pi}
= f^2 + \frac{\mu-\mu_0}{4\pi} \, .
\end{equation}
According to this formula $\mu \not= \mu_0$ shifts
the energy scale $4 \pi f^2(\mu)$ away from its original value
$4 \pi f^2$. In order to exhibit the scale--dependence it is
now convenient to present our results as function of
the ratio $f^2(\mu)/f^2$ rather than $\mu$ itself.
The presentation
is then independent of $\mu_0$ and also robust
against changes in the regu\-larisation scheme
which may introduce numerical factors to the scale (e.g. in the
3-momentum instead of the 2-momentum cutoff scheme
the scale is changed by a factor of $\pi/2$).
The second part of the above decomposition, $(\mu-\mu_0)$,
remains explicit in the finite Casimir energy

\begin{eqnarray}\label{casimir1}
E^{cas} (\mu) &=& -\frac{1}{2\pi} \int_0^\infty d p \frac{p}{\omega}
\left( \delta (p) - a_0 p^2 - a_1 \right)
+ \frac{1}{2} \sum_B (\omega_B - m) \\
&& + \frac{a_0}{6\pi} \left[ (\mu-\mu_0)^3 -
\frac{3m^2}{2}(\mu-\mu_0)-2m^3 \right]
- \frac{a_1}{2\pi} \left[ \mu - \mu_0 - m \right]
\, . \nonumber
\end{eqnarray}
Also for the $a_0$ - term with its cubic divergence
we kept all contributions from the second part
in the Casimir energy although other prescriptions are possible.
However this ambiguity will not influence the results strongly
because of the smallness of $a_0$ (cf. numbers in the previous
section). In this way both, the tree  and the 1-loop contributions
become scale--dependent. While for the N$\ell \sigma$ and potential
terms the scale--dependence in the total soliton energy is exactly
compensated, the 1-loop contribution which arises from
the Skyrme term would not only yield a scale--dependent Skyrme
parameter $e(\mu)$, but would also switch on all higher gradient
terms assumed to be zero for $\mu=\mu_0$. Having assumed all
these couplings to be independent of the scale implies
that the tree + 1-loop contribution to the soliton energy
cannot be strictly scale--invariant. In fact, the resulting
scale--dependence measures the magnitude of the higher gradient terms not
accounted for. All the more it comes as a surprise, just as in the
hadronic case~\cite{mw97}, that for specific parameter choices an
almost scale--independent soliton mass may be obtained. Such a case
($a^2=0.1, y=1$) is depicted in Fig.6 for the 1-soliton. The rather
strong scale--dependence of the tree contribution is nicely
compensated over a wide range $f^2 (\mu)=0.5 f^2, \dots , 2.0 f^2$
when the 1-loop contribution is added.
Of course, this statement has to depend on the parameters used.
Therefore in Fig. 7 we plotted the total 1-soliton mass for $a^2=0.1$
and various values of the ratio $y$. It is noticed that weak
scale--dependence requires a ratio close to $y \simeq 1$,
the value used
in the previous figure.
\begin{figure}[htbp]
\centerline{\epsfig{figure=baby6.ps,width=7cm,angle=270}}
\protect\caption{Scale-dependence of the 1-soliton's energy
in tree and tree + 1-loop. The dimensionless parameters are
$a=a(\mu_0)=.316$ and $y=y(\mu_0)=1.0$. The energy is
in units $4\pi f^2$ and the scale is given by the ratio
$f^2(\mu)/f^2$.
}
\centerline{\epsfig{figure=baby7.ps,width=7cm,angle=270}}
\protect\caption{Scale-dependence of the 1-soliton's energy
in tree + 1-loop for $a=a(\mu_0)=.316$ and
various parameters $y=y(\mu_0)$.
Energy and scale as in Fig. 6.
}
\end{figure}
Lower and higher values of $y$ enhance the
scale--dependence. We do not show the corresponding plots for the
2-soliton because they look quite similar, with the scale--dependence
being even somewhat weaker in that case.

We may now pose the question for which parameter combinations $a$ and
$y$ we may expect a modest scale--dependence in accordance with the
assumption that the couplings of all higher gradient terms are zero.
The answer is illustrated in Fig.8.
\begin{figure}[htbp]
\centerline{\epsfig{figure=baby8.ps,width=7cm,angle=270}}
\protect\caption{Parameter space $(a,y)$ for 1-skyrmions.
Inside the inner contour the scale--dependence is less than $3\%$,
inside the outer contour less than $6\%$.
}
\end{figure}
Inside the inner contour the average scale--dependence in the intervall
$f^2(\mu) / f^2 \in [0.5,2]$ is less than $3\%$. For comparison
the $6\%$ contour is also plotted.

All the parameter combinations lying inside the contours
yield negative Casimir energies whose absolute values though sizeable
do not exceed $50\%$ of the classical soliton mass. They also
lead to an
unstable 2-soliton which decays into two individual 1-solitons.
In particular the parameters $(a=0.1, y=0.9)$ of our numerical example
from hadron physics (section 4.1) lie in the
center of the almost scale--independent region.

The weak scale dependence obtained
for parameter combinations lying inside the contours in Fig. 8
implies that loop--corrections
can be reliably calculated within the framework
of the lagrangian (\ref{lag}) with all higher couplings set
equal to zero.

\section{Conclusions}

We studied the magnon--vortex system in the $2+1$ dimensional
$O(3)$ model in a field theoretical approach. For that purpose the
N$\ell \sigma$ model was augmented by stabilizing standard fourth--order
Skyrme and potential terms.

Complete information on bound and scattering states in all
partial--waves was established. We find an extremely rich excitation
spectrum with pronounced resonances as in the 3D-$SU(N_f)$ case with
its baryon resonances. In principle these resonances should be
accessible by measuring the excitation spectrum of spin--waves
in the presence of skyrmions.

Furthermore, the quantum fluctuations will allow to study their
influence on the shape
of the soliton~\cite{a97}.
Note that in tree approximation the
third component of the order parameter field is always
tied to $-1$ at the origin for topological reasons,
in variance with microscopic Hartree-Fock calculations
for the ferromagnet~\cite{HF97}
which suggest a vanishing value in case of small solitons.

Finally with the bound state energies and the scattering
phase--shifts the Casimir energies were evaluated. The
appearing UV singularities were renormalized under the
assumption that the (renormalized) couplings of the higher
gradient terms be small. A criterium for the consistence of this
assumption is the approximate scale--independence
of the results such that the effective action represents
an almost 1-loop renormalizable theory. This criterium requires a detailed
balance of tree and 1-loop contributions
which limits the parameter space where this requirement is met.
The following results apply for this
restricted parameter space. Parameters lying outside
that range may lead to other conclusions, but there the
1-loop contribution cannot be reliably calculated within the model.
\begin{itemize}
\item  The Casimir energy is negative and generally large leading
       to a considerable reduction of the total soliton energy.
       However, it does not exceed $50\%$ of
       the tree contribution and thus may still be considered
       a (sizeable) correction.
\item  The magnitude of the Casimir energy decreases with increasing
       soliton size. But it also decreases as the stabilizing terms
       become more important relative to the N$\ell \sigma$ model
       term.
\item  With the 1-loop corrections included the $n=2$ soliton becomes
       unstable and decays into two individual $n=1$ solitons (which
       may still be weakly bound by a dipole force). We conjecture
       that the same situation might occur for the 3D-$SU(2)$
       $B=2$ torus in hadron physics.
\end{itemize}
It is obvious that the program outlined in this paper has to be
tailored for specific applications concerning ferromagnets and
antiferromagnets. Most importantly for ferromagnets the
time derivative part of the lagrangian has to be replaced by the
$T$ violating Landau-Lifshitz term with one time derivative only.
This replacement may change the results for the Casimir energies
appreciably. Also, for antiferromagnets an external magnetic field
should couple through the time component of the
covariant derivative. Finally, proper consideration of the
non--local Coulomb interaction complicates the situation: instead
of coupled differential equations coupled integro--differential
equations have to be solved for the fluctuations.

The evaluation of loop corrections as presented here naturally suggest
the inclusion of temperature into the formalism. This opens the
possibility to study the properties of 2D spin--textures at finite
temperature, which is of particular interest near the
symmetry restoring phase--transition.

\bigskip
\noindent
This work is supported in parts by funds provided by the FCT, Portugal
(Contract PRAXIS/4/4.1/BCC/2753).

\section*{Appendix}

In this appendix we present analytical soliton solutions for
the two limiting cases $a \to 0$ and $a \to \infty$.

\subsubsection*{Belavin-Polyakov solution for $a \to 0$}

The analytical solution of the Euler-Lagrange equation
(\ref{stability}) in the limit $a \to 0$
$$
\frac{1}{x} \left( x F^\prime \right)^\prime - \frac{n^2 s c}{x^2} =0
\eqno{({\rm A}.1)}
$$
with boundary conditions (\ref{boundary}) is the well--known
Belavin-Polyakov soliton
$$
\tan \frac{F}{2} = \left( \frac{b_n}{x} \right)^n \, .
\eqno{({\rm A}.2)}
$$
The size parameter $b_n^2 = \sqrt{ 2n(n^2-1)/3 }$
is determined by the interplay of the Skyrme- and potential terms
contributing to the static energy functional
$$
E_n^0 = n + a \, \frac{\pi}{3b_n^2} \, \frac{n^2-1}{\sin \pi/n }
      + a \, \frac{b_n^2}{2} \, \frac{\pi/n}{\sin \pi/n }
\, , \qquad n \not= 1 \, .
\eqno{({\rm A}.3)}
$$
The corresponding square radius follows from its definition
(\ref{radius})
$$
< x^2 >_n = \sqrt{\frac{2n}{3} (n^2-1) } \, \frac{\pi/n}{\sin \pi/n }
\, \qquad n \not= 1 \, .
\eqno{({\rm A}.4)}
$$
It is noticed that in case of the 1-soliton the contribution
from the potential term to ({\rm A}.3) diverges. For that reason
we modify ({\rm A}.2),
$$
\tan \frac{F}{2} = b_n \,
\sqrt{ \frac{\epsilon}{e^{\epsilon x^2}-1}} \, ,
\eqno{({\rm A}.5)}
$$
such that for $\epsilon \to 0$ the Belavin--Polyakov solution is
recovered. The variation of the static energy functional in lowest
order $\epsilon$
$$
E_1^0 = 1 + \frac{\epsilon}{2} + \frac{2a}{3b_1^2}
- \frac{ab_1^2}{2} \ell n \epsilon
\eqno{({\rm A}.6)}
$$
gives then $\epsilon = ab_1^2$ and $\ell n \epsilon =  - 4 / 3b_1^4$
with the result
$$
E_1^0 = 1 + \epsilon ( \frac{1}{2} - \ell n \epsilon ) \, , \qquad \qquad
a^2 = - \frac{3}{4} \epsilon^2 \ell n \epsilon \, .
\eqno{({\rm A}.8)}
$$
Indeed with $a$ also $\epsilon$ tends to zero. The square radius
behaves like
$$
<x^2>_1 = \sqrt{-\frac{4}{3} \ell n \epsilon }
\eqno{({\rm A}.9)}
$$
and diverges weakly in the limit $a \to 0$. For
$a \stackrel{<}{_\sim} 0.001$ the radius of the 1-soliton exeeds
that of the 2-soliton. Exactly this behaviour may be traced in the
numerical solution if the differential equation (\ref{stability})
is solved for very small $a$ with great precision. Because the
square radius ({\rm A}.9) equals the derivative of the soliton
energy with respect to the parameter $a$, this also explains the
infinite slope at the origin indicated in the 1-soliton's curve
plotted in Fig. 2.

\subsubsection*{Solution for $a \to \infty$}

The Euler-Lagrange equation (\ref{stability}) in the opposite
limit $a \to \infty$
$$
\frac{n^2 s}{x} \left( \frac{F^\prime s}{x} \right)^\prime - s = 0
\eqno{({\rm A}.10)}
$$
may also be solved analytically
$$
\cos F(x) = \frac{1}{8n^2} (b_n^2 - x^2) x^2 + 2 \frac{x^2}{b_n^2} - 1
\, , \qquad x \le b_n \, .
\eqno{({\rm A}.11)}
$$
The integration constants are fixed by the condition $F(b_n)=0$
connected with topological charge $n$. The size parameter $b_n^2=4n$
again minimizes the static energy functional
$$
E_1^0 = na \left[ \frac{b_n^2}{4n} +
\frac{1}{2} \left( \frac{4n}{b_n^2} \right)
- \frac{1}{6} \left( \frac{b_n^2}{4n} \right)^3 \right]
\eqno{({\rm A}.12)}
$$
with the results quoted in (\ref{limit})
$$
E_n^0 = \frac{4n}{3} a \, , \qquad \qquad
<x^2>_n =  \frac{4n}{3} \, .
\eqno{({\rm A}.13)}
$$
Needless to say, that numerical solutions of the
Euler-Lagrange equation (\ref{stability}) reproduce these
results for large enough $a$.