the skyrme model is a classical field theory which describes atomic nuclei and their low energy interactions @xcite .
it is a non - linear theory admitting topological solitons , so - called skyrmions , labelled by a topological charge @xmath4 which can be identified with the baryon number .
when the theory is quantized the skyrmions give a good description of nuclei , @xmath5 resonance @xcite and even bound states of nuclei @xcite . a striking feature
is that skyrmions can be quantized as fermions @xcite .
when the theory is coupled to a fermion field two different ways of describing fermions appear in the same theory .
physically , the fermion field can be thought of as light quarks in the presence of atomic nuclei @xcite .
hiller calculated the spectrum of the fermions in the background of a skyrmion @xcite .
however , the approximation for the skyrmion configuration was very crude .
the spectrum shows a spectral flow of the eigenvalues depending on the value of the coupling constant .
in particular one mode crosses zero .
this general behaviour has been verified by other authors @xcite . to gain a better understanding we consider a generalization of the model .
since static field configurations in the original skyrme model in flat space are topologically equivalent to maps @xmath6 , the model can be generalized to the base space being a sphere of radius @xmath7 .
the original model is recovered in the limit @xmath8 .
it turns out that for @xmath9 and baryon number @xmath10 the faddeev - bogomolny bound is satisfied and the minimal energy solution can be worked out explicitly @xcite . in this case
the map @xmath6 is an isometry . for larger @xmath7 the skyrmions
are localized and tend to the usual flat space skyrmions in the limit @xmath8 .
monopoles and skyrmions have many properties in common , and it is therefore instructive to recall the results for monopoles .
there is an index theorem making use of the bogomolny equations which states that for @xmath11 monopoles there is a @xmath11 dimensional vector space of fermionic zero modes @xcite .
furthermore , the monopole and its fermionic zero mode can be calculated explicitly for @xmath12 .
we will examine how far we can extend the analogy between monopoles and skyrmions .
in this paper we study fermions coupled to skyrmions on @xmath13 . the dirac equation on @xmath13
is derived in sect .
[ diracequation ] .
we also explain the skyrmion - fermion interaction . in sect .
[ diracsolution ] , we derive the ordinary differential equation which describes the solutions with total angular momentum @xmath2 for general shape functions . in sect .
[ differentialequations ] , we give a brief discussion of ordinary differential equations and their regular singular points . in sect .
[ constantskyrmion ] , the equation of sect .
[ diracsolution ] is solved for a constant skyrme configuration , which is equivalent to a mass term for the fermions .
we also discuss the boundary conditions which are necessary to derive the spectrum . in sect .
[ hedgehogskyrmion ] , we derive the eigenfunctions and the spectrum for the hedgehog configuration with the special shape function @xmath14 . in sect .
[ numerics ] the results of sect .
[ constantskyrmion ] and [ hedgehogskyrmion ] are verified numerically .
the spectrum for more general shape functions is also evaluated numerically and compared with results in the literature .
we find an interesting relationship between the topological charge of the skyrmion and certain singular points .
we end with a conclusion .
in this section we derive the dirac equation on a @xmath0-sphere of radius @xmath9 .
first we define an atlas of two charts @xmath15 and @xmath16 which are the stereographic projections from the north pole @xmath17 and the south pole @xmath18 , respectively .
next we define vierbeins on both charts which uniquely determine a connection one - form using the metric compatibility condition and the torsion - free condition .
the connection one - form induces the spin connection so that we can finally write down the dirac equation for each chart . at the end of this section
we discuss the interaction term which couples fermions and skyrmions .    in order to construct the stereographic projection
let @xmath13 be embedded in euclidean @xmath19 with coordinates @xmath20 .
the points of @xmath13 can now be labelled by the projection from the north pole @xmath17 onto the @xmath21 plane through the equator introducing the coordinates @xmath15 .
alternatively , one can project from the south pole @xmath18 introducing the coordinates @xmath16 .
these two charts are well defined everywhere apart from their projection point @xmath17 and @xmath18 , respectively . in terms of the @xmath19 cartesian
coordinates we obtain @xmath22    in the following , we derive the equations of motion for both charts . as we shall see in sect .
[ constantskyrmion ] , this enables us to discuss the regularity conditions at the north and south pole in a simple way .
we will also use the two charts to discuss a nontrivial symmetry of the energy spectrum in sect .
[ generalization ] .
the induced metric on @xmath13 can be calculated from the euclidean metric in @xmath19 .
taking into account that our base space is a lorentzian space time @xmath23 , the metric can be written as @xmath24 the other coordinate system gives rise to the same metric replacing @xmath15 by @xmath16 .
on the overlap the two charts are related via the transition functions @xmath25 therefore , the coordinate vectors @xmath26 and @xmath27 are related by @xmath28 note that this transformation has the determinant @xmath29 .
we now choose non - coordinate bases @xmath30 and @xmath31 .
it is convenient to choose diagonal vierbeins @xmath32 and @xmath33 such that @xmath34 where all other components vanish .
the minus sign occurs because @xmath15 and @xmath16 are related by inversion which changes the orientation .
the vierbeins satisfy @xmath35 where @xmath36 is the flat minkowski metric .
the same relations are true for the vierbeins @xmath33 .
the non - coordinate basis vectors are related by an @xmath37-transformation @xmath38 this is in fact a @xmath39-rotation about the vector @xmath40 . with this choice of vierbeins
the connection one - form @xmath41 can be calculated from the metric compatibility condition @xmath42 and the torsion - free condition @xmath43 where @xmath44 is the dual basis of @xmath45 . after a short calculation
we obtain @xmath46 and the same relation holds for @xmath47 expressed in @xmath16-coordinates .
the spin connection @xmath48 is now given by @xmath49 where @xmath50 $ ] , and the @xmath51-matrices satisfy the standard anticommutation relations @xmath52 .
the lagrangian for fermions in curved space time is given by @xmath53 with our choice of coordinates and vierbeins we obtain @xmath54 using the coordinates @xmath15 and @xmath55 for the coordinates @xmath16 . equation ( [ l2 ] ) can be transformed into ( [ l1 ] ) by using @xmath56 and the induced spinor transformation @xmath57 of the fermion field @xmath58 which we calculate in the following . the spinor transformation @xmath57 is defined as @xmath59 the transformation @xmath57 is a lifting of the local lorentz transformation @xmath60 of the vierbeins in ( [ erotation ] ) .
note that the transformation @xmath61 can be rewritten as @xmath62 for @xmath63 . here
@xmath64 . since we are only concerned with the spatial part we work with the euclidean @xmath65 rather than with @xmath66
infinitesimally , the rotation is given by @xmath67 .
this can be lifted to a spinor transformation by @xmath68 .
for @xmath69 we obtain @xmath70 note that there is a sign choice .
as the spinor transformation @xmath57 is an element of the double cover of the rotation group , the transformation @xmath71 is also a lifting of the rotation @xmath72 .    in the following we discuss fermions coupled to skyrmions .
the full lagrangian @xmath73 of this model is the sum of the fermion lagrangian @xmath74 , the skyrmion lagrangian @xmath75 and the interaction lagrangian @xmath76 . in this paper
we consider fermions in the background of a static skyrme configuration and neglect the backreaction . therefore , we do not discuss the skyrmion lagrangian any further , and the interested reader is referred to @xcite .
the interaction lagrangian @xmath76 given by @xmath77 where @xmath78 is a convenient parameterization of the skyrme field .
equation ( [ lint ] ) is the standard parity - invariant manner in which fermions couple to a skyrme field and is also called the yukawa interaction term .
this lagrangian was first introduced by gell - mann and levy ( without a fourth order skyrme term ) , @xcite .
also see @xcite for a discussion of various different chiral models .
recently , the lagrangian ( [ lint ] ) has been discussed by balachandran _
_ @xcite .
note that @xmath79 is a spin and an isospin spinor . in our notation
the pauli matrices @xmath80 act on the isospin indices , whereas the @xmath51-matrices act on the spin indices .
it is convenient to split the dirac spinor into two @xmath81 spin isospin matrices @xmath82 and @xmath83 such that @xmath84 in this notation the @xmath51-matrices can be written in terms of quaternionic @xmath81 matrices using @xmath85 times the spin pauli matrices @xmath86 and the two dimensional identity matrix @xmath87 as a basis for the quaternions . since any complex @xmath81 matrix can also be expressed in terms of quaternions with complex coefficients it is convenient to choose the same basis for the quaternions as before
. then the spin isospin matrices @xmath88 can be written as @xmath89 . with this notation
the spin operators act on @xmath79 by left - multiplication , @xmath90 , whereas the isospin matrices act on @xmath79 by right multiplication @xmath91 where the transpose in ( [ quat1 ] ) guarantees the group structure , and ( [ quat2 ] ) gives a useful trick to evaluate @xmath92 .
as mentioned above , we make the approximation that the skyrmion is a fixed static background field , _ i.e. _ we ignore the backreaction of the fermion fields .
furthermore , we only consider spherically symmetric skyrmions . for @xmath10 ,
skyrmions are believed to be spherically symmetric , but for @xmath93 this is no longer true .
however , we expect that some generic features of the skyrmion - fermion interaction are not sensitive to the precise shape of the skyrmion
. for more details see sect .
[ generalization ] .
spherically symmetric skyrme fields are best written in polar coordinates @xmath94 where @xmath95 is the `` radial '' shape function and @xmath96 is the normal vector in @xmath97-direction , see ( [ unitvectors ] ) .
the shape function @xmath95 can be calculated by solving the euler - lagrange equation of the skyrme lagrangian numerically . however , instead of using a numerical solution it is more convenient in this context to approximate @xmath95 by an explicit function .
manton showed in @xcite that a good ansatz for @xmath7 not too large is the conformal ansatz given by the ansatz tends to @xmath98 .
qualitatively , this is still a good ansatz . ]
@xmath99 here @xmath100 is a constant which has to be determined by minimizing the energy of the skyrmion .
this ansatz describes a skyrme configuration which is either localized around one of the poles or uniformly distributed .
the value @xmath101 corresponds to the @xmath102-symmetric ansatz of the previous section , whereas for @xmath103 , or @xmath104 , the skyrme field is localized at the south pole , or the north pole , respectively .
manton also showed that for @xmath105 the stable solution is @xmath101 .
finally , we can write down the dirac equation for fermions coupled to a spherically symmetric background skyrmion . it is given by @xmath106 using the coordinates @xmath15 and @xmath107 using the coordinates @xmath16 . here
the skyrmion gives rise to    @xmath108    and the same equation holds in the coordinate system @xmath16 .
in order to solve the dirac equation we introduce spherical polar coordinates for equation ( [ dirac1 ] ) and ( [ dirac2 ] ) . in terms of the coordinates
@xmath109 in @xmath19 , the polar coordinates for the @xmath0-sphere of radius @xmath110 can be written as @xmath111 the angles @xmath112 can now be expressed in terms of @xmath15 as @xmath113 in terms of the @xmath114 , the relations for @xmath115 and @xmath116 are identical to the ones above , and for the @xmath97 coordinate we obtain @xmath117 equation ( [ dirac1 ] ) gives rise to the following dirac equation in polar coordinates : @xmath118 where @xmath119 and the unit vectors @xmath120 are given by @xmath121 @xmath122 can be written as @xmath123 .
similarly , equation ( [ dirac2 ] ) give rise to @xmath124 the spin connection term could be removed by setting @xmath125 and @xmath126 , respectively .
however , it is easier to discuss the regularity of the solution in terms of the original fields .
we adopt the standard representation of the @xmath51-matrices where @xmath127 is diagonal .
this is convenient because the dirac equation is invariant under parity .
@xmath128 now , the spinor transformation @xmath57 in ( [ rhospin ] ) can be evaluated explicitly : @xmath129 with this choice of @xmath51-matrices , we can derive the differential equation for the stationary solutions of the dirac equation by setting @xmath130 : @xmath131 where @xmath132 the components of the matrix in equation ( [ eq1 ] ) commute with the total angular momentum operator @xmath133 where @xmath134 is the orbital angular momentum operator , the spin operator is @xmath135 and the isospin operator is @xmath136 . in polar coordinates ,
the orbital angular momentum @xmath134 is given by @xmath137 equation ( [ eq1 ] ) is also invariant under parity @xmath138 with @xmath139 it is worth discussing the parity operator @xmath138 in more detail . @xmath138
acts on the @xmath51-matrices and the partial derivatives as follows : @xmath140 using the spinor transformation @xmath57 in ( [ rhospin ] ) we can show that @xmath141 therefore , the parity operator @xmath142 which is defined for the coordinates @xmath16 differs from @xmath138 by a minus sign .    the @xmath143 positive parity spin isospin spinors are given by @xmath144 this gives rise to the following system of first order differential equations : @xmath145 where @xmath146 , @xmath147 and @xmath148 are functions of @xmath97 .
note that the @xmath149 spinor @xmath150 gives rise to the same equation with @xmath151 replaced by @xmath152 .    in @xmath16-coordinates
positive and negative parity are interchanged such that the @xmath153 is given by @xmath154 this gives rise to @xmath155 again the @xmath156 spinor @xmath157 gives rise to the same equation with @xmath151 replaced by @xmath152 .
equation ( [ general1 ] ) can be transformed into a second order differential equation for @xmath158 where @xmath159 .
note that since @xmath160 $ ] it is consistent to set @xmath161 .
@xmath162 @xmath163 and similarly equation ( [ general2 ] ) can be transformed into @xmath164 @xmath165 it is not easy to transform equation ( [ od1 ] ) and ( [ od2 ] ) into each other .
in fact , due to the spinor transformation ( [ trafo ] ) @xmath166 and @xmath167 are related via @xmath168 also note that @xmath167 is given by @xmath169 and similarly @xmath170 can be written as @xmath171
in this section we solve the dirac equation which has been derived in the previous section for various shape functions @xmath95 . however , first we review some basic facts about ordinary differential equations and regular singular points , following @xcite .
any linear second order differential equation can be brought into the standard form @xmath172 if @xmath173 and @xmath174 are analytic functions in the neighbourhood of @xmath175 then @xmath176 is called a regular point . in this case
there are locally two linearly independent solutions @xmath177 and @xmath178 such that any solution of ( [ diffeq ] ) is a linear combination of @xmath177 and @xmath178 .
moreover , @xmath179 is an analytic function , _
i.e. _ it can be written as @xmath180 if either @xmath181 or @xmath182 is singular then @xmath175 is called a singular point .
in the special case that @xmath183 where @xmath184 and @xmath185 are analytic functions , _
i.e. _ @xmath186 the point @xmath175 is called a _
regular singular point_. equations which only have regular or regular singular points are known as fuchsian differential equations .
equation ( [ diffeq ] ) can then be solved by the following ansatz : @xmath187 the exponent @xmath57 can be calculated by solving the so - called _ indicial equation _
@xmath188 which is derived by inserting the ansatz ( [ eqansatz ] ) into ( [ diffeq ] )
. the two solutions @xmath189 of the indicial equation ( [ indicial ] ) are called the _ exponents _ of the regular singular point @xmath190 .
our convention will be that @xmath191 .
if the difference of the exponents @xmath192 is not a positive integer or zero , then the ansatz ( [ eqansatz ] ) gives rise to two linearly independent solutions of ( [ diffeq ] ) which are uniformly convergent in some neighbourhood of @xmath190 .
note , however , that the solution diverges at @xmath190 if its exponent @xmath57 is negative .
let @xmath192 be a positive integer or zero .
this situation is slightly more subtle because the two solutions can become linearly dependent , and this could lead to further logarithmic terms .
the solution @xmath177 with exponent @xmath193 is always uniformly convergent in some neighbourhood of @xmath190 .
therefore , we reduce equation ( [ diffeq ] ) to a differential equation of first order by setting @xmath194 .
this gives rise to the following equation for @xmath195 : @xmath196 this equation can be solved by integration and separation of variables , and we obtain @xmath197 where @xmath198 is an analytic function in a suitable domain , and @xmath199 and @xmath200 are constants .
a short calculation shows that @xmath201 let @xmath202 then for @xmath203 the most general solution is @xmath204 where @xmath199 and @xmath4 are arbitrary constants , and @xmath177 is the solution of the form ( [ eqansatz ] ) .
the constant @xmath205 and the function @xmath206 are completely determined by the differential equation ( [ difeqz ] ) .
the function @xmath207 is given by @xmath208 the coefficients @xmath209 being constants . in the special case that @xmath210 no logarithms occur , _ i.e. _ the two solutions @xmath177 and @xmath178 are linearly independent .
if @xmath211 the corresponding form of the solution is @xmath212    it is also possible to discuss the singularity at infinity . the transformation @xmath213 maps the point at infinity to the origin .
if we apply this transformation to equation ( [ diffeq ] ) we obtain a new differential equation in terms of @xmath214 .
if the new equation has a regular or a regular singular point at the origin @xmath215 then ( [ diffeq ] ) has a regular or regular singular point at infinity .
if a fuchsian differential equation has @xmath11 finite regular singular points @xmath216 with exponents @xmath217 and a regular singular point at infinity with exponents @xmath218 then the following equation holds : @xmath219 this is a useful consistency check for the exponents .
we are interested in the behaviour of fermions in the presence of a background skyrmion .
the skyrmion introduces a space - dependent complex mass for the fermions .
therefore , it is instructive to consider first a massive fermion on @xmath13 .
this corresponds to @xmath220 so that the fermion mass is @xmath151 .
this example can be solved explicitly in terms of hypergeometric functions . in the following sections we proceed to more complicated configurations . in sect .
[ hedgehogskyrmion ] , we consider the spherically symmetric shape function @xmath14 which can also be solved explicitly , but in a less standard way . in sect .
[ generalization ] , the conformal shape function ( [ conformal ] ) is discussed . in this case
, the differential equation can no longer be solved explicitly .
we present numerical algorithms in sect .
[ numericalalgorithm ] .    setting @xmath221 , equation ( [ od1 ] ) can be written as @xmath222 where @xmath223 .
equation ( [ od1f0 ] ) has three regular singular points , one at @xmath224 , one at @xmath225 and one at @xmath226 with exponents @xmath227 , @xmath228 and @xmath229 respectively .
note that relation ( [ erelation ] ) between the exponents is satisfied , _ i.e. _ the sum of the exponents is equal to @xmath110 . in this example
we will show that the solution is non - singular on @xmath13 if its exponents are @xmath230 and @xmath231 .
equation ( [ od1f0 ] ) can be solved explicitly : @xmath232 where @xmath233 is a hypergeometric function , _
e.g. _ see @xcite .
the hypergeometric function is non - singular for @xmath234 and is normalized at @xmath235 , _
i.e. _ @xmath236 .
it may be singular at @xmath237 , corresponding to @xmath238 , which is the north pole .
as @xmath158 has to be non - singular at the south pole , @xmath239 , we have to set @xmath240 .
since @xmath158 is not well defined at @xmath224 we have to express @xmath158 in terms of @xmath166 and calculate @xmath241 . using the transformation ( [ trafo ] ) from @xmath242 to @xmath79 we
obtain @xmath243 which is the analogue of equation ( [ fgtilde ] ) .
setting @xmath244 , in equation ( [ od2 ] ) we can derive a similar equation for @xmath245 which has exponents @xmath246 , @xmath247 and @xmath248 . the solution of this equation can again be expressed in terms of hypergeometric functions : @xmath249 this is only non - singular at @xmath224 if we set @xmath250 . assuming that @xmath251 , @xmath170 can be calculated from equation ( [ tildef ] ) : @xmath252
therefore , @xmath158 can be written as @xmath253 @xmath166 has to be non - singular at @xmath224 .
equation ( [ g(1)nonsing ] then implies that @xmath158 is also non - singular at @xmath224 .
therefore , the regular solution @xmath158 has the exponents @xmath230 and @xmath231 .
note that this argument only relies on the local behaviour and explicit solutions are not needed .
as we shall see , in the generic case , the exponents at the singular points do not depend on the chosen shape function , so that our analysis is also valid for general shape functions .
problems can occur if additional singular points coalesce with the singular points at @xmath254 because this changes the exponents .. ] in the generic case , the regular solution can be calculated by imposing that the solution can be expanded as a power series at the north and at the south pole . later , this will be the basis for our numerical algorithm .
now , we impose this condition in order to calculate the spectrum of equation ( [ od1f0 ] ) .
using formula ( 15.3.6 ) in @xcite we can express @xmath158 in terms of hypergeometric functions at the second singular point . ) only holds if @xmath255 .
however , this inequality is satisfied because @xmath256 . ]
@xmath257 the last term has to vanish , otherwise @xmath170 would be singular at the north pole . therefore , either @xmath258 or @xmath259 has to have a pole . since @xmath260
, we obtain @xmath261 , and the spectrum is therefore @xmath262 with this condition the eigenfunctions @xmath263 reduce to jacobi polynomials ( @xcite , formula 15.4.6 ) .
the spectrum is displayed in figure [ fermion1 ] .
it is symmetric with respect to @xmath264 .
the importance of this result lies in the fact that for @xmath265 the energy @xmath266 and the corresponding solutions @xmath263 are the same for all possible shape functions @xmath95 .
therefore , these solutions can be used as starting values for the relaxation method , which will be described in sect .
[ numerics ] .
as the radius of the @xmath0-sphere tends to infinity the spectrum tends to the usual spectrum of the dirac equation , which is continuous for @xmath267 . also note that if we assume that @xmath158 and its first and second derivative remain finite as @xmath151 tends to infinity then the energy @xmath268 has to be of order @xmath151 in this limit , so @xmath269 is the leading order asymptotic for @xmath268
. therefore , we displayed the lines @xmath270 in figure [ fermion1 ] . as we shall see in sect .
[ numerics ] the region @xmath271 also plays a special role for our numerical algorithm .     as a function of the coupling constant @xmath151 for @xmath244 .
@xmath272 is shown as dashed lines .
[ fermion1],height=472 ]    another way of deriving the spectrum ( [ specf=0 ] ) is to note that if the exponents at @xmath273 both vanish then @xmath158 is an integral function , that is it has only one singular point at infinity .
the only integral functions which do not possess an essential singular point at infinity are the polynomials , _
e.g. _ @xcite .
therefore , one of the exponents at infinity is equal to @xmath274 for @xmath275 .
since the exponents at infinity are @xmath276 we obtain the same spectrum as before .
note however , that since the exponents @xmath277 differ by an integer , we have to check in this approach that the solution is indeed a polynomial , _
i.e. _ @xmath205 in equation ( [ gs ] ) vanishes .
this spectrum has also been calculated by other authors , _ e.g. _ @xcite .      in this section
we construct the solutions of equation ( [ od1 ] ) for @xmath14 .
configurations with @xmath278 are called `` hedgehog configurations '' because at each point the vector in isospin space is proportional to the coordinate vector . for @xmath279
the map from the base @xmath13 to the target @xmath13 is even an isometry and possesses @xmath102-symmetry . setting @xmath14 in equation ( [ od1 ] )
we obtain @xmath280 @xmath281 this is a homogeneous differential equation of fuchsian type with four singular points .
they are located at @xmath273 , @xmath282 and @xmath283 .
assuming @xmath284 the exponents are @xmath285 for @xmath286 the singular points at @xmath287 and @xmath288 coalesce , and therefore , the exponents at @xmath287 become @xmath289 and the exponent at infinity is @xmath290 .
the relation for the exponents ( [ erelation ] ) is satisfied .
for @xmath291 the singular points at @xmath224 and @xmath292 coalesce .
this changes the exponents at @xmath224 to @xmath293 . in this case
, there is an essential singular point at infinity and the equation is no longer fuchsian .
since we are interested in solutions that can be expanded as a power series around @xmath273 , it is natural to look for polynomial solutions .
the singular point at @xmath294 might impose a constraint . in order to verify whether there is a logarithmic term we have to calculate the coefficient @xmath205 in equation ( [ gs ] ) . here
after a standard calculation we obtain @xmath296 . therefore , there are two independent power series at the point @xmath297 .
this imposes no further restriction , since we do not have to exclude any of these two solutions .
this implies that @xmath158 is again an integral function and we can calculate the energy spectrum using the exponents at infinity . again , the exponents @xmath277 differ by an integer , and we have to verify that the solutions are indeed polynomials .
it turns out that the only problem occurs for @xmath298 .    in the following ,
we derive the recurrence relation for these polynomials .
yet , it is instructive to consider first the solution where @xmath158 is a constant @xmath299 , _ i.e. _ @xmath298
. then we can solve equation ( [ exk=1 ] ) for @xmath268 and obtain @xmath300 and @xmath301 . setting @xmath302 in equation ( [ g(u)k=1 ] ) , we obtain @xmath303 it is easy to see that this is not a solution for @xmath304 .
@xmath302 is only a solution for arbitrary @xmath151 if @xmath305 .
this solution is very interesting because its spectrum crosses @xmath264 for @xmath306 .
now we derive the general solution for equation ( [ g(u)k=1 ] ) .
since we are interested in the solutions for which the relevant exponent vanishes for @xmath287 , we expand @xmath158 around the south pole @xmath287 : @xmath307 assuming that @xmath308 , we derive the following recurrence relations .
the initial condition is @xmath309 and the recurrence relation is given by @xmath310 since we have an initial condition we expect to be able to reduce this second order recurrence relation to a recurrence of first order . therefore , we make the ansatz @xmath311 substituting this ansatz into the recurrence relation ( [ recurrence ] ) we obtain a non - linear first order recurrence for the @xmath312 .
however , this can be solved by setting @xmath313 by construction the series ( [ ansatza ] ) has the right behaviour at the singular point @xmath287 . to find the behaviour at @xmath224 we have to discuss the radius of convergence @xmath314 of the series .
the radius @xmath314 can be calculated as the following limit : @xmath315 therefore , the series converges for @xmath316 . for @xmath224 , equation ( [ ansatza ] ) can be written as @xmath317 for large @xmath11 , @xmath312 behaves like @xmath318 .
therefore , @xmath319 is proportional to @xmath320 so that the terms in the sum ( [ g(1 ) ] ) are of order @xmath110 and do not vanish .
hence , the series is divergent , unless it terminates .
but because of the special form of @xmath319 and @xmath312 it is easy to see when the series terminates .
in fact , for @xmath321 we can now write @xmath319 as the following product : @xmath322 where @xmath323 is the product of all odd numbers less than or equal to @xmath324 .
all @xmath319 vanish for @xmath325 which is the constant solution that we have derived previously .
if the term @xmath326 vanishes then the coefficient for @xmath319 vanishes individually .
however , if one of the terms in the product vanishes , say for @xmath327 , then not only @xmath328 vanishes but also all @xmath329 for @xmath330 , _
i.e. _ the series terminates .     as functions of the coupling constant @xmath151 for @xmath14
@xmath272 is shown as dashed lines.[fermion2],height=472 ]    from this condition the complete energy spectrum can be calculated .
we obtain the special level @xmath331 which crosses from positive to negative energies , and two series of energy levels @xmath332 and @xmath333 which are given by @xmath334 the spectrum is displayed in figure [ fermion2 ] .
if we only consider the modes with energy @xmath335 the spectrum is symmetric under reflection at the axes @xmath336 and @xmath337 .
however , the mode @xmath338 with energy @xmath339 is only symmetric under combined reflections @xmath340 and @xmath341 , which corresponds to a @xmath39-rotation around the point @xmath342 .
this spectrum agrees with our previous calculation where we used the fact that @xmath158 is an integral function .
in this section we first describe numerical algorithms to calculate the spectrum for more general shape functions @xmath95 .
then we discuss the spectrum for the conformal ansatz ( [ conformal ] ) .
we also discuss the simple ansatz @xmath343 which gives rise to skyrmions with @xmath93 .
finally , we compare our results with the literature .
the energy spectrum can be calculated for various shape functions using a standard relaxation method @xcite .
the differential equation ( [ od1 ] ) is rewritten as a system of first order differential equations .
@xmath344 the energy @xmath268 is considered as an additional dependent variable whose derivative vanishes , @xmath345 , so that the complete system becomes @xmath346 for simplicity , let us discuss the numerical algorithm for @xmath102-symmetric shape function @xmath347 .
equation ( [ g(u)k=1 ] ) has a regular singular point at @xmath287 which implies that at this point we can not choose value and derivative of the solution independently .
however , once the value of @xmath348 is chosen we can calculate @xmath349 .
in fact , inserting @xmath287 into equation ( [ g(u)k=1 ] ) we obtain @xmath350 similarly we obtain for @xmath224 , @xmath351 the system ( [ system ] ) is discretized such that the boundary conditions at the singular points @xmath273 are satisfied . starting from an initial configuration
the relaxation algorithm calculates an improved configuration by linearizing the equation around the solution and minimizing an error function which indicates how well the discretized equations are satisfied .    the secret of a successful relaxation method is a good guess for the starting configuration .
fortunately , for @xmath265 the analytic solution of ( [ od1 ] ) is known for all possible shape functions @xmath95 .
therefore , we start with this configuration at @xmath265 and increase @xmath151 by a small amount @xmath352 .
the relaxed solution is a good initial configuration for @xmath353 , and so on .    as a test
, the spectrum for @xmath244 is calculated using the relaxation method .
this reproduces the analytic result . for @xmath14
the relaxation method converges well as long as @xmath354 is not inside the interval @xmath355 $ ] .
however , if the regular singular point @xmath354 is inside the interval , the relaxation algorithm no longer converges and has to be modified .
this is one of the main reasons for including the lines @xmath272 in the figures .    in this case
a specially chosen shooting method proves more suitable .
it deals with the additional singularities in @xmath355 $ ] by simply avoiding them . as mentioned before
, our differential equation can be defined on the complex plane . due to the boundary condition
the solution @xmath158 is an analytic function .
the original solution is recovered by restricting @xmath158 to be a real function of a real variable @xmath356 $ ] . in order to avoid further singular points on the real line
we can set @xmath357 where @xmath358 $ ]
. then equation ( [ od1 ] ) induces a differential equation for @xmath116 , that is the equation is solved on a semi - circular contour .
the resulting differential equation can be written as a system of 5 real first order differential equations , one of which corresponds to the energy @xmath268 .
now , the spectrum can be obtained by a standard shooting method , @xcite .
the boundary conditions for the shooting method are that at the singular points @xmath359 the solution is non - singular and furthermore that the equation is real at @xmath239 .
since the equation has real coefficients , the fact that @xmath158 is real for @xmath287 implies that @xmath158 is real for @xmath356 $ ] .
a successful shooting method needs a good initial guess . for @xmath360 we can rely on our explicit solution and impose @xmath361 and @xmath362 where equation ( [ initial1a ] ) follows from equation ( [ hypid ] ) and equation ( [ initial1b ] ) follows from the reflection relation for the gamma function ( see ( 6.1.17 ) in @xcite ) .
as in the relaxation method the initial guess for the coupling constant @xmath363 is the value of the solution for coupling constant @xmath151 .
it is also possible to reconstruct the solution @xmath158 for @xmath356 $ ] using the solution @xmath364 for @xmath365 $ ] .
cauchy s formula states that if @xmath366 is an analytic function and @xmath367 is a contour enclosing the point @xmath368 then @xmath369 see for example @xcite . in our case @xmath158 is a real function on the interval @xmath355 $ ] , so @xmath370 where the coefficients @xmath371 are real .
it follows that @xmath158 can be written as @xmath372 and a similar formula can be obtained for its derivative @xmath373 . for @xmath374 both the relaxation method and the shooting method give the same result .      with the conformal shape function
( [ conformal ] ) , equation ( [ od1 ] ) gives rise to a rather lengthy differential equation with five regular singular points , see @xcite for more details .
assuming that the singular points do not coalesce we obtain the following exponents : @xmath375 where @xmath376 and @xmath377 and @xmath378 . the relation ( [ erelation ] ) is again satisfied .
moreover , the exponents at infinity are the same as for @xmath244 . in the limit @xmath379 the exponents at infinity change
discontinuously ( see ( [ exk=1 ] ) ) .
the exponents change if one or more singular points coalesce .
however , in general , there are five singular points which all have different characteristics .
the regular singular points at @xmath359 provide the boundary conditions for the regular solution , and are not very sensitive to the exact form of the shape function .
they imply that the solution of this equation locally still is an analytic function of @xmath380 .
one might expect that the regular singular point at infinity could be used to calculate the degree of the polynomial as in the previous two sections .
however , for general @xmath100 , the solutions are not regular on the entire complex plane , so that the exponents at infinity do not contain any information about the spectrum .
this is because there is a further singular point at @xmath381 with exponents @xmath382 .
generically , the solution diverges logarithmically at this singular point such that @xmath383 is a branch point , see equation ( [ log ] ) .
the existence of polynomial solutions might be related to the spherical @xmath102 symmetry of the problem , which is broken for @xmath384 .
note that the singular point at @xmath385 is independent of @xmath268 , and it can never be inside the interval @xmath386 because @xmath387 for @xmath388 .
finally , the singular point at @xmath389 depends on the energy @xmath268 .
moreover , for certain values of @xmath268 , @xmath100 and @xmath151 this point can be inside the interval @xmath386 .
this singular point could force a potential solution to diverge , such that it would no longer be well - defined .
therefore , we are interested in the local behaviour near @xmath390 . after a straightforward but long calculation using maple , we obtain @xmath391 .
therefore , locally , both solutions are convergent power series , and no logarithmic divergencies occur .    ) for the shape function @xmath392 where @xmath393 .
@xmath272 is shown as dashed lines.[fermion3],height=472 ]    in the following we calculate spectra as functions of the coupling constant @xmath151 for different values of @xmath100 numerically .
for @xmath101 we reproduce the spectrum for @xmath14 of sect .
[ hedgehogskyrmion ] .
also note that @xmath394 corresponds to the spectrum for @xmath395 which agrees with the spectrum for @xmath244 in sect .
[ constantskyrmion ] . as @xmath100 increases the spectrum moves towards the @xmath396 spectrum .
at the same time , the energy @xmath339 intersects @xmath264 for increasing values of the coupling constant @xmath151 , such that for @xmath397 the intersection point also appears to move to @xmath398 . in figure [ fermion3 ]
we display the spectrum for @xmath393 .    in the following
we discuss the eigenfunction which crosses zero in more detail .
we restrict out attention to the following lorentz scalar @xmath399 which is defined as @xmath400 a short calculation shows that @xmath401 is a function of @xmath97 only .
therefore , we can integrate out the angular coordinates @xmath115 and @xmath116 and obtain : @xmath402 the density @xmath403 measures where the fermion field is localised . to obtain equation ( [ qtilde2 ] ) we used equations ( [ 0+spinor ] ) and ( [ fgtilde ] ) . @xmath158 and @xmath166 can be calculated from ( [ od1 ] ) and ( [ od2 ] ) , respectively .
recall @xmath404 .
expressing @xmath167 in terms of @xmath166 avoids having to worry about the denominator in ( [ f ] ) which is rather difficult to handle numerically .
note however , that @xmath167 is non - singular for @xmath405 .
this follows from the behaviour at the corresponding regular singular point of equation ( [ od1 ] ) .
the function @xmath166 is related to @xmath158 via equations ( [ fgtilde ] ) and ( [ f ] ) .
this implies @xmath406 in the following , we calculate @xmath403 for @xmath101 , @xmath407 and @xmath408 and compare it to the baryon density @xmath409 which measures where the skyrme field is localized .
is defined with a minus sign , but this is of no relevance here . ]    the expression for @xmath79 can be calculated explicitly for @xmath101 . in this case @xmath302 .
since we are interested in the zero crossing mode we have to set @xmath410 and obtain @xmath411 surprisingly , this is independent of @xmath151 .
also note that in this case @xmath403 is proportional to @xmath412 . in the following ,
we normalize @xmath403 so that @xmath413 .    for @xmath414 ,
the skyrme configuration is localised near the north pole @xmath415 ( see figure [ fermion4 ] ) . in figure [ fermion4 ] ( a ) and
( b ) we show @xmath412 for @xmath416 and @xmath393 , respectively . in figure [ fermion4 ] ( c ) and ( d )
we show @xmath403 for several values of the coupling constant @xmath151 . for @xmath265 , @xmath403
can also be calculated explicitly and is independent of @xmath100 .
for @xmath378 and @xmath417 , the fermion configuration is also localised near the north pole . for @xmath416 , the energy @xmath268 vanishes when @xmath418 . for @xmath393 ,
the energy @xmath268 vanishes when @xmath419 . in figure [ fermion4 ] ( e ) and
( f ) we display @xmath403 for these values of @xmath151 and compare it to the baryon density @xmath412 .
there is quite a good agreement between these two functions .
therefore , the @xmath264 states can be interpreted as a fermion - skyrmion bound states .    the simplest method for obtaining a skyrmion with baryon number
@xmath93 is the ansatz @xmath420 . for @xmath93
this skyrme configuration has very high energy
. however , we are mainly interested in the interplay of topology and the spectrum of the dirac equation which should be independent of the precise form of the skyrmion .    .
@xmath272 is shown as dashed lines .
[ fermion5],height=472 ]    for @xmath421 the spectrum is displayed in figure [ fermion5 ] .
it is symmetric with respect to the @xmath268-axis .
note that the symmetry of the spectrum can be derived from equation ( [ od1 ] ) and ( [ od2 ] ) .
for even @xmath4 the function @xmath422 is an even function of @xmath380 and @xmath423 is an odd function of @xmath380 . under @xmath424 and @xmath425 equation ( [ od1 ] ) transforms into equation ( [ od2 ] ) .
however , the equations are just expressed in different coordinate systems , and hence have the same energy eigenvalues .
therefore , we have shown that the spectrum is symmetric under @xmath424 if @xmath4 is even , provided that @xmath148 is an odd function of @xmath380 .    for @xmath265
all the eigenmodes have the same value , independent of the background configuration , and parity . as @xmath151 increases , the modes increasingly feel the influence of the background configuration , in particular its topology .
it appears that the winding number @xmath421 of the skyrmion forces @xmath407 modes to cross zero , one for @xmath417 and one for @xmath426 . in order to follow the spectrum for @xmath271 the shooting method has to be used .    .
here we only show the spectrum for @xmath427 , see text .
@xmath272 is shown as dashed lines.[fermion6],height=472 ]    this pattern can also be verified for @xmath428 , see figure [ fermion6 ] .
now there are three modes which cross zero , two of positive parity ( @xmath417 ) and one of negative parity @xmath429 .
it seems probable that the number zero - crossing modes of positive parity differs at most by one from the number of zero - crossing modes of negative parity .
further calculations for @xmath430 and @xmath408 confirm this pattern . for @xmath430
we find two zero - crossing modes of positive parity and two of negative parity , whereas for @xmath431 there are three positive parity modes and two negative parity modes .
this numerical evidence shows that the singular points of equation ( [ od1 ] ) play an important role . in the following
we will assume that @xmath95 is an analytic function , and in particular that it can be expanded in a taylor series with respect to the variable @xmath159 , for @xmath432 .
it follows that @xmath422 and @xmath423 are also analytic functions of @xmath380 .
therefore , the only singular points of ( [ od1 ] ) inside the interval @xmath386 occur for @xmath433 .
if such a singular point is non - degenerate its exponents are @xmath338 and @xmath407 .
furthermore , it can be shown be shown by expanding @xmath148 as a power series in @xmath380 that the coefficient @xmath434 always vanishes . in our previous examples @xmath95 was always monotonic .
let @xmath95 be monotonic with boundary conditions @xmath435 and @xmath436 . if @xmath271 then @xmath433 has @xmath4 nondegenerate solutions for @xmath437 . in this case
, these singularities of equation ( [ od1 ] ) all have exponents @xmath438 and the coefficient @xmath434 vanishes , so that the solution of ( [ od1 ] ) is regular at these points .
we conjecture that the existence of @xmath4 regular singular points inside the interval @xmath386 indicates that there are @xmath4 modes crossing from one side of the spectrum to the other .
this conjecture is equivalent to the statement that for topological charge @xmath4 there are @xmath4 modes that cross zero .
if we allow for more general skyrme configurations than the hedgehog , such a statement becomes even more difficult to prove , because now we have to solve partial differential equations .
however , both the topological number and the number of singular points are stable under small perturbations .
therefore , there seems to be a deep interplay between the topology of the configuration and the singular points of a differential operator .
this can be compared to monopoles where we have an index theorem which connects the topology of the monopole configuration with the index of the dirac operator .    in the following we compare our results with the literature .
kahana and ripka considered the problem of coupling a hedgehog configuration to fermions in flat space using an approximation for the skyrme configuration with a size parameter . varying this parameter they computed the spectral flow of the eigenvalues .
they found a @xmath439 zero mode for a special value of the `` size '' of the hedgehog . since the size is given by minimizing the energy density for the skyrmion their result corresponds to varying the coupling constant @xmath151 .
in @xcite a linear approximation for the shape function , @xmath440 , was used . here
@xmath441 so that the dirac equation for the fermions is independent of @xmath151 . however , we defined the skyrmion as a given background configuration which minimizes the energy . therefore , @xmath442 where @xmath443 is the size of the skyrmion .
@xmath443 can be calculated by minimizing the energy @xcite : @xmath444 .
kahana and ripka showed by numerical calculation that there is a zero mode for @xmath445 which corresponds to @xmath446 . in @xcite an exponential function was used to approximate the shape function : @xmath447 . in this case minimizing the energy results in @xmath448 and the zero mode occurs for @xmath449 which corresponds to @xmath450 . the fact that there is a mode which crosses zero seems to be independent of the special ansatz for the shape function . this is consistent with our results .
in @xcite it has been argued that a skyrmion in flat space is reasonably well approximated by the stereographic projection of a skyrmion on a @xmath0-sphere of optimal radius .
the optimal radius for @xmath10 is @xmath9 so we can use our analytic result for the mode @xmath339 in sect .
[ hedgehogskyrmion ] . for @xmath306
this is a zero mode , and the value of @xmath151 is in the same order of magnitude as in the literature mentioned above .
we have derived the dirac equation for fermions on the @xmath0-sphere which are chirally coupled to a skyrmion . in this paper
, we only considered spherically symmetric skyrme configurations and fermionic states with the total angular momentum @xmath451 .
we were particularly interested in the shape functions @xmath244 and @xmath14 and the conformal ansatz ( [ conformal ] ) .
the shape function @xmath221 leads to a hypergeometric equation and can therefore be solved explicitly in terms of hypergeometric functions . to derive the spectrum we have to impose that the solutions are non - singular at the north and south pole using the correct chart .
we derived the spectrum @xmath452 which agrees with the spectrum calculated by sen @xcite using maurer cartan forms .
moreover , the eigenfunctions @xmath263 are given by jacobi polynomials .
the shape function @xmath14 is more complicated .
it leads to a fuchsian equation with four regular singular points , two at the boundary @xmath359 , one at infinity and one depending on @xmath268 and @xmath151 .
however , the equation could still be solved in terms of polynomials .
we were able to solve the recurrence relation for the coefficients of the polynomial and derived the corresponding spectrum : @xmath453 the eigenfunctions have been evaluated explicitly : @xmath454 where @xmath455 we also considered more general shape functions numerically .
for @xmath10 we discussed the conformal ansatz ( [ conformal ] ) and calculated the spectrum .
for @xmath101 the spectrum agrees with the analytic solution .
in particular there is a mode @xmath338 whose energy @xmath339 crosses zero for a certain value of the coupling @xmath151 .
the graph of @xmath339 as a function of @xmath151 changes as the parameter @xmath100 in the ansatz is varied .
however , for finite @xmath456 , this mode always stays present . according to our numerical evidence
the eigenfunction for @xmath264 can be interpreted as a fermion - skyrmion bound state .    for larger @xmath4 we considered the ansatz @xmath420 .
this leads to @xmath4 modes crossing zero .
we conjectured a connection between zero modes and certain singular points in the equation which in turn are connected to the winding number of the configuration .
the author wants to thank n s manton , h pfeiffer and y shnir for many discussions during various stages of this project .
the author also would like to thank j m speight for many helpful suggestions and continuous support .
t.  h.  r. skyrme , _ a nonlinear field theory _ , proc .
. lond . * a260 * : 127 ( * 1961 * ) , g.  s. adkins , c.  r. nappi and e.  witten , _ static properties of nucleons in the skyrme model _ , nucl
. phys . * b228 * : 552 ( * 1983 * ) , r.  a. leese , n.  s. manton and b.  j. schroers , _ attractive channel skyrmions and the deuteron _ ,
. phys . * b442 * : 228 ( * 1995 * ) , , p.  irwin , _
zero mode quantization of multi - skyrmions _ ,
phys . rev .
* d61 * : 114024 ( * 2000 * ) , , s.  krusch , _ homotopy of rational maps and the quantization of skyrmions _ ( * 2002 * ) , , d.  finkelstein and j.  rubinstein , _ connection between spin , statistics , and kinks _ , j. math . phys .
* 9 * : 1762 ( * 1968 * ) , e.  witten , _ global aspects of current algebra _
, nucl . phys . * b223 * : 422 ( * 1983 * ) , a.  p. balachandran and s.  vaidya , _ skyrmions , spectral flow and parity doubles _ , int
a14 * : 445 ( * 1999 * ) , , j.  r. hiller and t.  f. jordan , _ solutions of the dirac equation for fermions in skyrme fields _ , phys
* d34 * : 1176 ( * 1986 * ) , s.  kahana and g.  ripka , _ baryon density of quarks coupled to a chiral field _ , nucl .
phys . * a429 * : 462 ( * 1984 * ) , g.  ripka and s.  kahana , _ the stability of a chiral soliton in the fermion one loop approximation _ , phys
. lett . * b155 * : 327 ( * 1985 * ) , g.  ripka , _ quarks bound by chiral fields _ , oxford science publications ( * 1997 * ) .    n.  s. manton and p.  j. ruback , _ skyrmions in flat space and curved space _ , phys . lett . * b181 * : 137 ( * 1986 * ) , n.  s. manton and b.  j. schroers , _ bundles over moduli spaces and the quantization of @xmath457 monopoles _ , ann
225 * : 290 ( * 1993 * ) , m.  nakahara , _ geometry , topology and physics _ , graduate student series in physics , iop publishing ltd ( * 1990 * ) .
t.  eguchi , p.  b. gilkey and a.  j. hanson , _ gravitation , gauge theories and differential geometry _ ,
. rept . * 66 * : 213 ( * 1980 * ) , s.  krusch , _
@xmath458 skyrmions and the rational map ansatz _ ,
nonlinearity * 13 * : 2163 ( * 2000 * ) , , m.  gell - mann and m.  levy , _ the axial vector current in beta decay _ , nuovo cim . * 16 * : 705 ( * 1960 * ) , n.  s. manton , _ geometry of skyrmions _ ,
commun . math .
phys . * 111 * : 469 ( * 1987 * ) , e.  t. whittaker and g.  n. watson , _ a course of modern analysis _ , cambridge mathematical library , cambridge university press , fourth edition ( * 1927 * ) .
d.  sen , _ fermions in the space - time @xmath459 _ , j. math .
* 27 * : 472 ( * 1986 * ) , m.  carmeli and s.  malin , _ field theory on @xmath459 topology .
3 . the dirac equation _ , found . phys . * 15 * : 1019 ( * 1985 * ) , w.  h. press , b.  p. flannery , s.  a. teukolsky and w.  t. vetterling , _ numerical recipes _ , cambridge university press ( * 1992 * ) .