the interaction of a charged particles beam with a plasma is an important subject of relevance for many fields of physics , such as inertial confinement fusion ( icf ) driven by ion or electron beams @xcite , high energy density physics @xcite , and related astrophysical phenomena @xcite .
this interaction is also relevant , among others , for the fast ignition scenario ( fis ) @xcite , where the precompressed deuterium - tritium ( dt ) core of a fusion target is to be ignited by a laser - generated relativistic electron beam .
in addition , a promising icf scheme has been recently proposed @xcite , in which the plasma target is irradiated simultaneously by intense laser and ion beams . within this scheme several experiments @xcite have been performed to investigate the interactions of heavy ion and laser beams with plasma or solid targets .
an important aspect of these experiments is the energy loss measurements for the ions in a wide - range of plasma parameters .
it is expected in such experiments that the ion propagation would be essentially affected by the parametric excitation of the plasma target by means of laser irradiation .
this effect has been supported recently by particle - in - cell numerical simulations @xcite .    motivated by the experimental achievements mentioned above , in this paper we present a study of the effects of intense radiation field ( rf ) on the interaction of projectile particles with an electron plasma .
previously this has been a subject of great activity , starting with the pioneering work of tavdgiridze , aliev , and gorbunov @xcite . more recently
the need of the comprehensive analysis of the complexity of the beam - matter interaction in the presence of strong laser fields has stimulated a number of theoretical studies .
@xcite and @xcite have developed a time - dependent hamiltonian formulation to describe the effect of a strong laser field on energy losses of swift ions moving in a degenerate electron gas .
it has been shown that the energy loss is reduced for ions at intermediate velocities , but is increased for slow ions , due to a resonance process of plasmon excitation in the target assisted by a photon absorption @xcite . for plasma targets , on the other hand
, it has been obtained that the energy loss is reduced in the presence of a laser field when the projectile velocity is smaller than the electron thermal velocity @xcite .
moreover , the reduction of the energy loss rate increases with decreasing the angle between the projectile velocity and the polarization vector of the laser field @xcite .
next , an interesting effect has been obtained by @xcite considering the oscillatory motion of the light projectile particles in the laser field . in this case
the energy loss depends on the projectile charge sign and mass and , therefore , results in a different stopping powers for electrons and positrons . in recent years
special attention has been also paid to the regimes of high intensities of the laser field when the quiver velocity of the plasma electron exceeds the projectile velocity @xcite .
it has been shown that the projectile particles can be even accelerated in a plasma due to the parametric effects excited by an intense laser field .
furthermore , @xcite have studied the influence of a high - intensity laser field on the coulomb explosion and stopping power for a swift h@xmath3 cluster ion in a plasma target .
they have demonstrated that the laser field affects the correlation between the ions and contributes to weaken the wake effect as compared to the laser - free case .
all papers mentioned above deal with a nonrelativistic projectile particles interacting with a laser irradiated plasma . in the present paper we study more general configurations thus considering the interaction of relativistic particles with an electron plasma .
our objective is then to study some aspects of the energy loss of relativistic particles which have not been considered previously . as an important motivation of the present paper is the research on the topic of the fis for inertial confinement fusion
@xcite which involves the interaction of a laser - generated relativistic electron beam ( reb ) with a plasma .
although the concept of the fis implies an overdense plasma and the propagation of a relativistic electron beam from the border of a precompressed target to the dense core occurs without crossing the laser beam , the target plasma is assumed to be parametrically excited by the rf through high - harmonics generation .
in this situation , the presence of the rf as well as the relativistic effects of the beam can dramatically change the main features of the standard stopping capabilities of the plasma .
the plan of the paper is as follows . in section  [ sec:2 ]
, we outline the theoretical formalism for the electromagnetic response of a plasma due to the motion of a relativistic charged particles beam in a plasma in the presence of an intense rf .
the stopping of the relativistic beam is formulated in section  [ sec:3 ] . in this section
we calculate the effects of the rf on the mean energy loss ( stopping power ) of the test particle considering two somewhat distinct cases with a weak ( sec .
[ sec:3 - 1 ] ) and intense ( sec .
[ sec:3 - 2 ] ) radiation field as well as a wide range for the test particle energies , extending from nonrelativistic up to the ultrarelativistic energies .
the dynamics of the electron plasma is treated by a system of fluid equations which results in a high - frequency approximation for the dielectric response function .
the results are summarized in section  [ sec:4 ] .
the whole interaction process of the relativistic charged particles beam with plasma involves the energy loss and the charge states of the projectile particle and  as an additional aspect  the ionization and recombination of the projectile driven by the rf and the collisions with the plasma particles .
a complete description of the interaction of the beam requires a simultaneous treatment of all these effects including , in particular , in the case of the structured particles the effect of the charge equilibration on the energy loss process . in this paper
, we do not discuss the charge state evolution of the projectiles under study , but concentrate on the rf effects on the energy loss process assuming an equilibrium charge state of the projectile particle with an effective charge @xmath4 .
this is motivated by the fact that the charge equilibration occurs in time scales , which are usually much smaller than the time of passage of the beam through target .
in addition , the effects mentioned above are not important in the case of the relativistic electron beams .
the problem is formulated using the hydrodynamic model of the plasma and includes the effects of the rf in a self - consistent way .
the external rf is treated in the long - wavelength limit , and the plasma electrons are considered nonrelativistic .
these are good approximations provided that ( i ) the wavelength of the rf ( @xmath5 ) is much larger than the debye screening length @xmath6 ( with @xmath7 the thermal velocity of the electrons and @xmath8 the plasma frequency ) , and ( ii ) the quiver velocity  of the electrons in the rf @xmath9 ( where @xmath10 is the amplitude of the rf ) is much smaller than the speed of light @xmath11 .
these conditions can be alternatively written as ( i ) @xmath12 , and ( ii ) @xmath13 , where @xmath14 is the intensity of the rf . as an estimate in the case of a dense plasma , with electron density @xmath15  @xmath16 , we get @xmath17  w/@xmath18 . thus the conditions ( i ) and ( ii ) are well above the values obtained with currently available high - power sources of the rf , and these approximations are well justified .
the response of the plasma to the incoming relativistic beam is governed by the hydrodynamic equations for the density @xmath19 and the velocity @xmath20 of the plasma electrons as well as by the maxwell equations for the electromagnetic fields .
thus , @xmath21 \right ] ,   \label{eq:2}\end{aligned}\]]where @xmath22 is the rf , and @xmath23 and @xmath24 are the self - consistent electromagnetic fields which are determined by the maxwell equations , @xmath25here @xmath26 is the induced current density , @xmath27 and @xmath28 are the charge and current densities for the beam , respectively , @xmath29 is the equilibrium plasma density in an unperturbed state .
as mentioned above we consider a relativistic beam of charged particles and , therefore , the influence of the rf @xmath30 on the beam is ignored .
consider now the solution of eqs .  - for the beam - plasma system in the presence of the high - frequency rf . in an unperturbed state ( i.e. , neglecting the self - consistent electromagnetic fields @xmath23 and @xmath24 in eq .   and assuming the homogeneous initial state ) the plasma velocity satisfies the equation @xmath31 which yields the equilibrium velocity for the plasma electrons , @xmath32 . here
@xmath33 and @xmath34 are the quiver velocity and the oscillation amplitude of the plasma electrons , respectively , driven by the rf .
next we consider the linearized hydrodynamic and maxwell s equations for the plasma and electromagnetic fields , respectively . for sufficiently small perturbations ,
we assume @xmath35 , @xmath36 , ( with @xmath37 and @xmath38 ) .
thus , introducing the fourier transforms @xmath39 , @xmath40 , @xmath41 , and @xmath42 with respect to the coordinate @xmath43 , the linearized hydrodynamic equations read @xmath44 \delta n(\mathbf{k},t ) = -in_{e}\left ( \mathbf{k}\cdot \delta \mathbf{v}(\mathbf{k},t)\right ) ,   \label{eq:5 } \\ & & \left [ \frac{\partial } { \partial t}+i\left ( \mathbf{k}\cdot \mathbf{u}_{e}(t)\right ) \right ] \delta \mathbf{v}(\mathbf{k},t ) = -\frac{e}{m}\left [ \mathbf{e}(\mathbf{k},t)+\frac{1}{c}\left [ \mathbf{u}_{e}(t)\times \mathbf{b}(\mathbf{k},t)\right ] \right ] .
\label{eq:6}\end{aligned}\]]in order to solve eqs .   and
it is convenient to introduce instead of the fourier transform @xmath45 a new unknown function @xmath46 via the relation @xcite
@xmath47where @xmath48 , @xmath45 represents one of the quantities @xmath49 , @xmath50 , @xmath41 , and @xmath42 . substituting this relation for the quantities @xmath39 , @xmath50 , @xmath41 , and @xmath42 into eqs .   and
it is easy to see that the obtained equation for the unknown function @xmath46 constitutes an equation with periodic coefficients @xcite .
therefore , we introduce the decompositions @xmath51where @xmath52 and @xmath53 are the corresponding amplitudes of the @xmath54th harmonic .
it is also useful to find the connection between these amplitudes .
this can be done using the fourier series representation of the exponential function @xmath55 as well as the summation formula @xmath56 for the bessel functions @xcite which yields @xmath57here @xmath58 is the bessel function of the first kind and of the @xmath59th order . using these results from eqs .   and
it is straightforward to obtain @xmath60where @xmath61 and we have neglected the term of the order of @xmath62 ( see the last term in eq .  ) according to the approximation ( ii ) .
using the transformation formulas and from the system of equations and it is easy to obtain @xmath63 next let us evaluate the induced current density which is given by @xmath64 $ ] in the coordinate space .
employing eqs .  , , and
for the fourier transform of the induced current density we obtain @xmath65 \right\ }   \label{eq:16 } \\ & & = \frac{i\omega _ { p}^{2}}{4\pi } \sum_{s , r=-\infty } ^{\infty } \frac{1}{\omega _ { r}}j_{\ell -r}(\zeta ) j_{s - r}(\zeta ) \left [ \mathbf{e}_{s}(\mathbf{k},\omega ) + \frac{(\ell -r)\omega _ { 0}}{\omega _ { r}}\frac{\mathbf{a}}{\zeta } \left ( \mathbf{k}\cdot \mathbf{e}_{s}(\mathbf{k},\omega ) \right ) \right ] , \nonumber\end{aligned}\]]where @xmath66 is the plasma frequency .    from the maxwell equations we express the magnetic field through the electric field .
in terms of the amplitudes of the @xmath54th harmonics [ see definition given by eqs .   and
] the electric field is determined by the system of equations @xmath67 ,   \label{eq:17 } \\ & & \left ( \mathbf{k}\cdot \mathbf{e}_{\ell } ( \mathbf{k},\omega ) \right ) = 4\pi i\left [ e\delta n_{\ell } ( \mathbf{k},\omega ) -\delta _ { \ell 0}\rho _ { 0}(\mathbf{k},\omega ) \right ] .
\label{eq:18}\end{aligned}\]]here @xmath68 and @xmath69 are the ordinary fourier transforms of the current and charge densities of the beam , respectively .
note that the @xmath59th harmonics of these quantities are given by @xmath70 and @xmath71 .
also the longitudinal part of the electric field is determined by the poisson equation .
insertion of eq .   into eq .
results in @xmath72 \mathbf{e}_{\ell } ( \mathbf{k},\omega ) -\mathbf{k } ( \mathbf{k}\cdot \mathbf{e}_{\ell } ( \mathbf{k},\omega ) )   \nonumber \\ & & = -\frac{\omega _ { p}^{2}}{c^{2}}\sum_{s , r=-\infty } ^{\infty } \frac{(\ell -r)\omega _ { 0}}{\omega _ { r}}j_{\ell -r}(\zeta ) j_{s - r}(\zeta ) \left [ \mathbf{e}_{s}(\mathbf{k},\omega ) + \frac{\omega _ { \ell } } { \omega _ { r}}\frac{\mathbf{a}}{\zeta } \left ( \mathbf{k}\cdot \mathbf{e}_{s}(\mathbf{k},\omega ) \right ) \right ]   \label{eq:19 } \\ & & + \frac{4\pi i\omega _ { \ell } } { c^{2}}\delta _ { \ell 0}\mathbf{j}_{0}(\mathbf{k},\omega ) ,   \nonumber\end{aligned}\]]where @xmath73 is the longitudinal dielectric function of a plasma .
the obtained expression completely determines the electromagnetic response in the beam - plasma system in the presence of the rf .
the resulting equation represents a coupled and infinite system of linear equations for the quantities @xmath74 ( with @xmath75 ) .
the ( infinite ) determinant of this system determines the dispersion equation for the electromagnetic plasma modes in the presence of the rf .
it should be also emphasized that the hydrodynamic description of a plasma is justified at high frequencies when @xmath76 . in particular
, this condition requires that the velocity of the projectile particle should exceed the thermal velocity @xmath7 of the plasma electrons .
equation   can be further simplified excluding there the longitudinal part @xmath77 of the electric field by means of the poisson equation together with the induced charge density , eq .  .
in the first step we multiply both sides of eq .   by @xmath78 and sum up over the harmonic numbers @xmath59 and @xmath79 . using the summation formula involving two bessel functions ( see above eq .  )
this results in @xmath80 similarly repeating the first step for the latter expression for the longitudinal electric field one finally obtains @xmath81 this expression has been previously obtained and studied by many authors @xcite assuming nonrelativistic beam of charged particles .
thus , inserting eq .   into eq .
we finally obtain @xmath72 \mathbf{e}_{\ell } ( \mathbf{k},\omega ) = \frac{4\pi i\omega } { c^{2}}\delta _ { \ell 0}\mathbf{j}_{0}\left ( \mathbf{k},\omega \right )   \nonumber \\ & & -\frac{\omega _ { p}^{2}}{c^{2}}\sum_{s , r=-\infty } ^{\infty } \frac{(\ell -r)\omega _ { 0}}{\omega _ { r}}j_{\ell -r}(\zeta ) j_{s - r}(\zeta ) \mathbf{e}_{s}(\mathbf{k},\omega )   \label{eq:22 } \\ & & -4\pi i\rho _ { 0}(\mathbf{k},\omega ) \sum_{r=-\infty } ^{\infty } ( -1)^{r}\frac{\boldsymbol{\chi } _ { \ell r}(\mathbf{k},\omega ) } { \varepsilon ( \omega _ { r})}j_{\ell -r}(\zeta ) j_{r}(\zeta ) ,   \nonumber\end{aligned}\]]where we have introduced the vector quantity @xmath82 \frac{\mathbf{a}}{c^{2}\zeta } .
\label{eq:23}\ ] ]    as mentioned above eq .   as well as its equivalent eq .
represent an infinite system of the coupled linear equations for the amplitudes @xmath83 which , in general , do not allow an analytic solution .
next , we will make an _ ad hoc _ assumption neglecting the second term in the right - hand side of eq .   which contains an infinite sum over the harmonics @xmath84 .
in this case the solution of eq .   is trivial and it is given by @xmath85 .
\nonumber\end{aligned}\]]now let us examine the physical consequences of the neglect of the second term in eq .  .
for this purpose we insert the last expression into the second term of the right - hand side of eq .  .
then it is straightforward to see that this term corresponds to the erenkov radiation ( or absorbtion ) in a plasma by a moving charged particles beam .
note that this effect is absent in the field - free case since the physical conditions of the erenkov radiation can not be fulfilled in this case .
thus the rf essentially changes the dispersion properties of the plasma and the erenkov effect becomes now possible .
in particular , an inspection of a general expression shows that in a long wavelength limit the emission of the transverse electromagnetic waves occurs at the frequencies @xmath86 for @xmath87 ( emission ) , and at @xmath88 for @xmath89 ( absorbtion ) , where @xmath90 and @xmath91 are the relativistic factor and the velocity of the beam , respectively .
moreover , the ratio of the neglected term , related to the erenkov effect , to the contribution given by eq .
is of the order of @xmath92 and the neglect of the second term in eq .
is justified at the high - frequencies of the rf ( @xmath93 ) . in the next sections
the further calculations will be done using the approximate expression .    before closing this section
let us consider briefly some limiting cases of eq .  .
for instance , in the electrostatic limit ( assuming that @xmath94 ) from eq .   one obtains @xmath95which confirms eq .  .
in the limit of the vanishing radiation field with @xmath96 , eq .
yields @xmath97 , where @xmath98 .
\label{eq:26}\]]it is seen that in the last expression @xmath99 is the fourier transform of the electric field created in a medium by an external charge with the density @xmath100 and current @xmath68 @xcite .
with the basic results presented in the previous sec .
[ sec:2 ] , we now take up the main topic of this paper , namely the investigation of the stopping power of a relativistic charged particle beam in a plasma irradiated by an intense laser field . for further progress the charge and current densities of the beam in eq .   should be specified .
furthermore , we will assume that these quantities are given functions of the coordinates and time .
in particular , for a beam of charged particles having a total charge @xmath4 and moving with a constant velocity @xmath101 , the charge and current densities are determined by @xmath102 and @xmath103 , respectively .
here @xmath104 is the spatial distribution of the charge in the beam with the normalization @xmath105 .
accordingly , the fourier transforms of these quantities are @xmath106q(\mathbf{k})\delta ( \omega -\mathbf{k}\cdot \mathbf{u})$ ] , @xmath107 , and @xmath108 is the fourier transform of the distribution @xmath104 ( form - factor of the beam ) .
note that @xmath108 is normalized in a such way that @xmath109 .
in particular , for the gaussian beam with a charge distribution @xmath110 , we obtain @xmath111 . here
@xmath112 and @xmath113 are the characteristic sizes of the beam along and transverse to the direction of motion , respectively .    from eqs .
and it is straightforward to calculate the time average ( with respect to the period @xmath114 of the laser field ) of the stopping field acting on the beam .
then , the averaged stopping power ( sp ) of the beam becomes @xmath115 \frac{(\mathbf{u}\cdot \mathbf{a})}{c^{2}\zeta } \right ] \right\ } ,   \nonumber\end{aligned}\]]where @xmath116 and @xmath90 is the relativistic factor of the beam . in eq .
the symbol @xmath117 denotes an average with respect to the period of the laser field .
hence , the sp depends on the beam velocity @xmath101 , the frequency @xmath118 and the intensity @xmath119 of the rf ( the intensity dependence is given through the quiver amplitude @xmath120 ) . moreover , since the vector @xmath121 in eq .
is integrated over the angular variables , @xmath122 becomes also a function of the angle @xmath123 between the velocity @xmath101 , and the direction of rf polarization , represented by @xmath120 .
in addition , it is well known that within classical description an upper cut - off parameter @xmath124 ( where @xmath125 is the effective minimum impact parameter or the distance of closest approach ) must be introduced in eq .   to avoid the logarithmic divergence at large @xmath126 .
this divergence corresponds to the incapability of the classical perturbation theory to treat close encounters between the projectile particles and the plasma electrons properly . for @xmath125 ,
we use the effective minimum impact parameter excluding hard coulomb collisions with a scattering angle larger than @xmath127 .
the resulting cut - off parameter @xmath128 is well known for energy loss calculations ( see , e.g. , @xcite ; @xcite and references therein ) .
here @xmath7 is the thermal velocity of the electrons .
in particular , at high - velocities this cut - off parameter reads as @xmath129 .
when , however , @xmath130 , the de broglie wavelength exceeds the classical distance of closest approach . under these circumstances
we choose @xmath131 .
equation   contains two types of the contributions .
the first one is proportional to the imaginary part of the inverse of the longitudinal dielectric function @xmath132 and is the energy loss of the beam due to the excitations of the longitudinal plasma modes .
the second contribution comes from the imaginary part of @xmath133^{-1}$ ] , where @xmath134 has been introduced above and determines the energy loss of the beam due to the excitations of the transverse plasma modes .
therefore , due to the resonant nature of the integral in eq .
, we may replace the imaginary parts of the functions @xmath135^{-1}$ ] and @xmath133^{-1}$ ] using the property of the dirac @xmath136 function . however , as discussed in the preceding section , the terms containing the factor @xmath137 with zero harmonic number ( @xmath138 ) correspond to the usual erenkov effect which is absent here .
neglecting such terms we finally arrive at @xmath139 \delta \left ( \varepsilon ( \omega _ { \ell } ) \right ) .
\nonumber\end{aligned}\]]by comparison , the sp in the absence of the rf is given by @xcite @xmath140which after straightforward integration for the point - like particles ( with @xmath141 ) yields @xmath142 in the presence of the rf , the sp @xmath143 is modified and is given by the @xmath138 term in eq .
( no photon  sp ) .
it is evident that the rf decreases the sp @xmath143 by a factor of @xmath144 .
finally , we would like to mention that the relativistic effects in eq .
are given by the second term in the square brackets and by the dispersion function @xmath134 . in a nonrelativistic limit ( @xmath145 )
the former vanishes while the latter is replaced by @xmath146 .    to illustrate the features of the stopping of the relativistic charged particles in a laser - irradiated plasma
, we consider below two examples when the intensity of the rf is small but the angle between the vectors @xmath101 and @xmath120 is arbitrary ( sec .
[ sec:3 - 1 ] ) .
in the second example we assume that the polarization vector @xmath10 of the rf is parallel to the beam velocity @xmath101 ( sec .
[ sec:3 - 2 ] ) .      in this section
we consider the case of a weak radiation field ( @xmath147 ) at arbitrary angle @xmath123 between the velocity @xmath101 of the beam and the direction of polarization of the rf @xmath120 . assuming a point - like test particle ( with @xmath141 ) , in eq .
we keep only the quadratic terms with respect to the quantity @xmath120 ( or @xmath148 ) and for the stopping power @xmath122 we obtain @xmath149 \right .
\label{eq:31 } \\ & & + \frac{|\omega _ { -1}|}{\omega _ { -1}}\frac{\delta \left ( \varepsilon ( \omega _ { -1})\right ) } { d(k,\omega ) } \left [ ( \mathbf{k}\cdot \mathbf{a})-\frac{\omega _ { p}^{2}}{\omega _ { -1}^{2}}\frac{(\mathbf{u}\cdot \mathbf{v}_{e})}{c^{2}\zeta } \right ]   \nonumber \\ & & -\frac{2|\omega |}{\omega } \frac{(\mathbf{k}\cdot \mathbf{a})}{k^{2}}\delta \left ( \varepsilon ( \omega ) \right)\bigg\ } ,   \nonumber\end{aligned}\]]where @xmath150 and @xmath143 is the field - free sp given by eq .  .
the remaining terms are proportional to the intensity of the radiation field ( @xmath151 ) .
note that due to the isotropy of the dielectric function @xmath152 as well as the dispersion function @xmath134 the angular integrations in eq .   can be easily done .
this results in @xmath153 , \label{eq:32}\]]where @xmath154 is the scaled quiver velocity and two functions have been introduced as follows @xmath155with @xmath156 , @xmath157 , @xmath158 .
also @xmath159 is the relativistic factor of the beam with @xmath160 . in particular , in the case of the ultrarelativistic beam with @xmath161 from eqs .
and we obtain @xmath162 ,   \label{eq:33a } \\ & & \phi _ { 2}(\nu ) = 2\ln ( 1 + \eta ^{2 } ) ,   \label{eq:34a}\end{aligned}\]]and in eq .   the additional energy loss of the beam ( i.e. the term of eq .
proportional to @xmath163 ) is mainly determined by the function @xmath164 .
consider next the angular distribution of the sp at low - intensities of the rf .
figure  [ fig:1 ] shows the dimensionless quantity @xmath165/s_{0}$ ] ( the relative deviation of @xmath122 from @xmath143 ) as a function of the angle @xmath166 for the scaled laser intensity @xmath167 and for @xmath168 , @xmath169 , @xmath170 , where @xmath171 is the electron classical radius .
the quantity @xmath172 corresponds to the density of the plasma @xmath15 @xmath16 .
in addition , left and right panels of fig .
[ fig:1 ] demonstrate @xmath173 for the stopping of the nonrelativistic ( with @xmath174 and @xmath175 ) and relativistic ( with @xmath176 , @xmath177 and @xmath178 ) charged particles , respectively . at nonrelativistic velocities the second term in eq .
can be neglected and the angular distribution of the quantity @xmath173 has a quadrupole nature ( the angular average of @xmath173 vanishes ) . at @xmath179 the excitation of the waves with the frequencies @xmath180 leads to additional energy loss . at @xmath181
the energy loss changes the sign and the total energy loss decreases .
when the particle moves at @xmath182 with respect to the polarization vector @xmath120 the radiation field has no any influence on the sp . at the relativistic velocities ( fig .
[ fig:1 ] , right panel ) the field - free sp @xmath143 as well as the relative deviation @xmath173 are strongly reduced and the critical angle @xmath183 , which now depends on the factor @xmath184 , is shifted towards higher values .
finally , at ultrarelativistic velocities ( fig .  [ fig:1 ] , dotted line in the right panel ) , @xmath173 is positive for arbitrary @xmath123 and the radiation field systematically increases the energy loss of the particle ( see also eq .   with eqs .   and ) .
let us now investigate the influence of the radiation field on the stopping process of the relativistic charged particle when its velocity @xmath101 is parallel to @xmath120 .
it is expected that the effect of the rf is maximal in this case . in the case of the point - like particles after straightforward integrations in eq .
we arrive at @xmath186 .
\label{eq:35}\ ] ] we have used here the same notations as in sec .  [ sec:3 - 1 ] . from eq .
it is seen that  no photon
sp ( the term with @xmath138 ) oscillates with the intensity of the laser field .
however , the radiation field suppresses the excitation of the plasma modes and this term by a factor of @xmath187 is smaller than the field - free sp @xmath143 . as follows from eq .   at high - intensities of the rf
the no photon  sp vanishes when @xmath188 with @xmath189 , where @xmath190 are the zeros of the bessel function @xmath191 ( @xmath192 , @xmath193 ,  ) .
then the energy loss of the charged particle is mainly determined by other terms in eq .   and is stipulated by excitation of plasma waves with the combinational frequencies @xmath194 .
in the case of the ultrarelativistic charged particle ( with @xmath161 ) eq .   gets essentially simplified . in this equation
neglecting the small corrections proportional to @xmath195 we obtain @xmath196with @xmath197 .
equation   can be further simplified neglecting the small quantity @xmath198 in the argument of the bessel function .
the remaining infinite sum is evaluated in appendix  [ sec : app1 ] [ see eq .  ] which results in @xmath199thus , in the case of the ultrarelativistic charged particles the energy loss of the beam essentially increases with the intensity of the rf .
let us recall , however , that eq .
is only valid below relativistic intensities of the rf when @xmath200 . for treating precisely the resonance in eq .
at @xmath201 we should go beyond the present long - wavelength approximation of the laser field .
figure  [ fig : inf ] demonstrates the quantity @xmath202 for three distinct values of the laser frequency @xmath1 as a function of the rf intensity @xmath197 , where the sp is given by eq .   for an ultrarelativistic projectile particle .
note that the quantity @xmath203 is independent on the plasma density @xmath29 .
it is seen that the sp is indeed not sensitive to the variation of the laser frequency and the approximation   is justified at @xmath200 .
next let us examine an interesting regime of the high - intensities of the rf when @xmath204 .
it is clear that such a regime requires nonrelativistic velocities of the charged particles beam .
thus , at @xmath205 we consider here the relativistic effects as small corrections to the corresponding nonrelativistic expressions .
using the asymptotic behavior of the bessel function @xcite in a general eq .
it is straightforward to obtain @xmath206 \right\ } ,   \nonumber\end{aligned}\]]where we have introduced the following quantities @xmath207 in appendix  [ sec : app1 ] we have found some approximate analytical expressions for these functions .
recalling that @xmath208 the relativistic corrections ( @xmath209 ) can be easily deduced from eq .  .
let us consider briefly two particular cases . when @xmath210 ( where @xmath211 ) all terms in eq
.   have a positive sign and the function @xmath212 in eq .
is maximal in this case due to the constructive interference of the excited waves with the frequencies @xmath194 .
note that the relativistic correction given by @xmath213 vanishes in this case . from eq .
it follows that @xmath214 and eq .
can be evaluated explicitly .
the result reads @xmath215}{\sinh [ \pi \gamma a(0)]}\right\ } , \label{eq:41a}\]]where @xmath216 , @xmath217 .    in the second regime
when @xmath218 , @xmath219 and the relativistic corrections given by @xmath220 again vanishes , @xmath221 .
equation   then yields a negative contribution to the energy loss with @xmath222}{\tanh [ \frac{1}{2}\pi \gamma a(\eta ) ] } \right\ } < 0 .
\label{eq:41b}\]]the function @xmath223 is minimal in this case since the different harmonics cancel effectively each other .
however , it should be emphasized that the contribution of the term containing the function @xmath223 to the sp can either be positive or negative depending on the dimensionless quantity @xmath224 .
moreover , the absolute value of this contribution decreases with the parameter @xmath224 .
the results of the numerical evaluation of the sp given by a general eq .   are shown in figs .
[ fig:2][fig:4 ] both for nonrelativistic and relativistic projectiles , where the ratios @xmath225 and @xmath226 are plotted as the functions of the rf intensity @xmath227 ( figs .
[ fig:2 ] and [ fig:3 ] ) and the relativistic factor @xmath228 of the electron beam ( fig .
[ fig:4 ] ) , respectively . here
@xmath229 and @xmath230 are measured in units @xmath231
w/@xmath18 and @xmath232 m , respectively . also in figs .
[ fig:2 ] and [ fig:3 ] the parameter @xmath233 varies in the interval @xmath234 to ensure the applicability of the long - wavelength approximation of the laser field . from fig .
[ fig:2 ] it is seen that the sp of a nonrelativistic particle exhibits a strong oscillations .
furthermore , it may exceed the field - free sp and change the sign due to plasma irradiation by intense ( @xmath235 ) laser field .
these properties of a nonrelativistic sp have been obtained previously for a classical plasma @xcite as well as for a fully degenerated plasma @xcite .
the effect of the acceleration of the projectile particle shown in fig .
[ fig:2 ] occurs at @xmath236 ( with @xmath237 ) when the  no photon
sp nearly vanishes .
it should be noted that in the laser irradiated plasma a parametric instability is expected @xcite with an increment increasing with the intensity of the radiation field .
this restricts the possible acceleration time of the incident particle .
the similar behavior of the sp is demonstrated in the right panel of fig .
[ fig:2 ] , where however the plasma density is @xmath238  @xmath16 and two orders of magnitude exceeds the value adopted in the left panel of fig .
[ fig:2 ] .
we note that although the ratio @xmath203 is reduced with a plasma density , the field - free sp given by eq .
is proportional to @xmath29 which in turn results in a strong enhancement of the sp @xmath239 with a plasma density .
next , in fig .
[ fig:3 ] , we consider weakly relativistic regimes of the electron beam with @xmath240 ( left panel ) and @xmath241 ( right panel ) . at smaller intensities of the rf , @xmath242 , the sp
is described qualitatively by eq .   and reaches its maximum at @xmath243 ( or at @xmath244 ) .
this is because the electrons are driven by the laser field in the direction of the beam with nearly same velocity @xmath91 effectively enhancing the amplitude of the excited waves .
the minimum or two maxima of the sp at @xmath245 in the left panel of fig .
[ fig:3 ] are formed due to the constructive interference of the excited harmonics as discussed above ( see eq .  ) .
the other minima of the sp at @xmath246 in the left panel of fig .
[ fig:3 ] are formed due to the destructive interference of the excited harmonics ( see eq .  ) .
let us note that at @xmath247 the effect of the enhancement of the sp of the test particle moving in a laser irradiated plasma is intensified at smaller frequency of the laser field @xmath248 ( see figs .
[ fig:2 ] and [ fig:3 ] ) while at higher intensities , @xmath235 , the quantity @xmath249 increases with the frequency @xmath1 .
finally , in fig .  [ fig:4 ] , the quantity @xmath250 is shown as a function of the electron beam energy at various intensities of the laser field given here by a dimensionless parameter @xmath62 .
it is seen that at the moderate energies of the electron beam the sp strongly increases with the laser intensity as has been also demonstrated in figs .
[ fig:2 ] and [ fig:3 ] . comparing two panels of fig .
[ fig:4 ] we may conclude that the position of the maximum of the sp depends only on the intensity of the rf , @xmath251 , while the maximum value of the sp is scaled as @xmath252 and is sensitive to both the intensity and frequency of the rf .
furthermore , with increasing the relativistic factor @xmath228 of the electron beam one asymptotically arrives at the ultrarelativistic regime described by eqs .   and .
interestingly this occurs at the moderate energies of the electron beam with @xmath253 which , in particular , pertain to the fis relativistic electron beam in the typical 1- to 2-mev energy range of practical interest .
in this paper , we have investigated the energy loss of a relativistic charged particles beam moving in a laser irradiated plasma , where the laser field has been treated in the long - wavelength approximation .
the dynamics of the beam - plasma system in the presence of the rf is studied by the linearized fluid and maxwell equations for the plasma particles and the electromagnetic fields , respectively .
the full electromagnetic response of the system is derived involving all harmonics of the rf .
it has been shown that , in general , the excited longitudinal and transversal modes are parametrically coupled due to the presence of the rf . as a result , the erenkov radiation by the charged particles beam , which is absent in a plasma in the field - free case , becomes possible .
however , in contrast to the usual erenkov effect in a laser - free medium , the spectrum of the radiation in the present case is discrete and the emitted frequencies are localized around the harmonics of the rf depending essentially on the energy of the incoming beam . in particular , at ultrarelativistic energies of the beam emitted frequencies are shifted towards higher harmonics of the rf .    in the course of our study
, we have derived a general expression for the sp , which has been also simplified in the limit of a low - intensity laser field . as in the field - free case ,
the sp in a laser irradiated plasma is completely determined by the dielectric function of the plasma which in the present context is given by the hydrodynamic approximation . explicit calculations have been done for two particular cases assuming ( i ) a low - intensity laser field but arbitrary orientation of the laser polarization vector @xmath10 with respect to the charged particles beam velocity @xmath101 , and ( ii ) an arbitrary ( but nonrelativistic ) intensity of the laser field when @xmath10 is parallel to @xmath101 . in the case
( i ) it has been shown that the rf increases the sp for small angles between the velocity of the beam and @xmath10 while decreasing it at large angles .
furthermore , the relative deviation of the sp from the field - free value is strongly reduced with increasing the beam energy . in the case ( ii ) and at high - intensities of the rf
two extrem regimes have been distinguished when the excited harmonics cancel effectively each other reducing strongly the energy loss or increasing it due to the constructive interference . as in the previous nonrelativistic treatments @xcite an acceleration of the projectile particle
is expected at high - intensities of the rf when the quiver velocity of the plasma electrons exceeds the beam velocity @xmath91 .
special attention has been paid to the relativistic effects of the beam .
we have demonstrated that in general the energy loss increases with the beam energy forming a maximum at moderate relativistic factors @xmath254 which , however , are shifted towards higher energies of the beam with the laser intensity .
in addition , the enhancement of the sp is more pronounced when the laser frequency approaches the plasma frequency in agreement with pic simulations @xcite .
finally , it has been also shown that the sp of the ultrarelativistic electron beam increases systematically with the intensity of the rf exceeding essentially the field - free sp .
going beyond the present model , which is based on several approximations , we can envisage a number of improvements .
in particular , a simple possibility qualitatively involving relativistic intensities of the laser field is to replace the quiver velocity @xmath255 by its relativistic counterpart @xmath256 , where @xmath257^{-1/2}$ ] .
the latter approximation can be easily obtained from a single - particle dynamics in a plane rf .
in the course of our study we have also neglected the erenkov radiation which might be important in the total balance of the energy loss process . however , in a general eq .
this effect is well separated by its spectral characteristics and can be treated independently .
we intend to address this issue in our forthcoming investigations .
the work of h.b.n . has been supported by the state committee of science of the armenian ministry of higher education and science ( project no .
13 - 1c200 ) .
in this appendix we evaluate the infinite sum containing the square of the bessel functions and involved in eq .  .
also we will derive some analytical expressions for the functions @xmath223 and @xmath213 introduced in sec .
[ sec:3 - 2 ] [ see eqs .   and ] .
our starting point is the neumann s formula @xcite @xmath258using this expression as well as the known relation @xcite @xmath259the desired summation can be performed easily @xmath260the latter formula is valid at @xmath261 .
next we derive approximate but very accurate expressions for the functions @xmath223 and @xmath213 which determine the sp in eq .  .
for that purpose consider the following two functions @xmath262 where @xmath263 and @xmath264 are some arbitrary positive variables .
note that @xmath265 at @xmath266 and @xmath267 .
consider the derivatives of these functions with respect to the variable @xmath268 .
using the known summation formulas involving the trigonometric functions @xcite we arrive at @xmath269}{\sinh(\pi u)}-\frac{1}{\pi u}\right\ } , \label{eq : ap4a } \\ & & \frac{\partial}{\partial a}g_{2}(a , b ) = 2a\sum_{\ell=1}^{\infty}\frac{(-1)^{\ell}}{\ell}\frac{\sin(2\alpha\ell)}{\ell^{2}+u^{2 } } = \frac{\pi a}{u^{2}}\cos\left [   \pi q(\alpha)\right ]   \left\ {   \frac { \sinh\left [   2\pi up(\alpha)\right ]   } { \sinh(\pi u)}-2p(\alpha)\right\ }   \label{eq : ap5a}\end{aligned}\ ] ] with @xmath270 , @xmath271 $ ] , and @xmath272where @xmath273 and @xmath274 are the fractional and integer parts of @xmath275 , respectively .
consequently , from eqs .
and we obtain @xmath276}{\sinh ( x ) } dx -\ln\frac{u}{u_{0 } } , \label{eq : ap7 } \\ & & g_{2}(a , b ) = g_{2}(a_{0},b)+\pi\cos [ \pi q(\alpha ) ] \left\{\int_{\pi u_{0}}^{\pi u } \frac{\sinh [ 2xp(\alpha ) ] } { \sinh ( x)}\frac{dx}{x } -2p(\alpha ) \ln\frac{u}{u_{0}}\right\ } .
\label{eq : ap8}\end{aligned}\ ] ] here @xmath277 and @xmath278 is an arbitrary initial value of @xmath268 .
thus , the infinite sums in eqs .   and have been represented in the equivalent integral forms with finite integration limits . assuming that @xmath279 let us apply eqs .   and in a particular case with @xmath280 .
taking into account the initial conditions ( see eqs .   and )
@xmath281 from eqs .
and we obtain @xmath282}{\sinh(x)}dx , \label{eq : ap9a }   \\ & & g_{2}(a , b ) = \pi\cos[\pi q(\alpha)]\left\ {   \int_{\pi b}^{\pi u}\frac{\sinh [ 2xp(\alpha ) ] } { \sinh ( x)}\frac{dx}{x } - 2p(\alpha)\ln\frac{u}{b } \right\ } .
\label{eq : ap10a}\end{aligned}\ ] ] for further progress we note that at @xmath279 the hyperbolic functions in eqs .   and can be replaced by an exponential function , @xmath283 , which yield the following approximations for the functions @xmath284 and @xmath285 @xmath286 -\exp\left [   -2\pi u q_{-}(\alpha)\right ]   \right\ }   \label{eq : ap9 }   \\ & & + \frac{1}{2q_{+}(\alpha)}\left\{\exp\left [   -2\pi bq_{+}(\alpha)\right ] -\exp\left[-2\pi u q_{+}(\alpha)\right ]   \right\ } ,   \nonumber   \\ & & g_{2}(a , b ) \simeq\pi\cos[\pi q(\alpha)]\big\ { e_{1}\left[2\pi u q_{+}(\alpha ) \right ] -e_{1}\left[2\pi u q_{-}(\alpha ) \right ]   \label{eq : ap10 }   \\ & & + e_{1}\left [   2\pi bq_{-}(\alpha)\right ]   -e_{1}\left [   2\pi bq_{+}(\alpha)\right ] -2p(\alpha)\ln\frac{u}{b } \big\ } , \nonumber\end{aligned}\ ] ] where @xmath287 is the exponential integral function and @xmath288 .
let us now consider the opposite case of small @xmath289 . then at @xmath290 we may neglect the variable @xmath291 in eqs .   and and represent these expressions in the approximate form @xmath292 here @xmath293 are the functions @xmath294 at @xmath295 and are evaluated in a similar way as approximate eqs .   and .
for that purpose consider again eqs .   and with @xmath296 and @xmath295 , that is @xmath297}{\sinh ( x)}dx -\ln a , \label{eq : ap11 }   \\ & & g_{2}(a,0 ) = g_{2}(1,0)+\pi\cos [ \pi q(\alpha ) ] \left\{\int_{\pi}^{\pi a}\frac{\sinh [ 2xp(\alpha ) ] } { \sinh ( x ) } \frac{dx}{x } -2p(\alpha ) \ln a\right\ } .
\label{eq : ap12}\end{aligned}\ ] ] in the latter expressions replacing the hyperbolic sine functions by the exponential functions we arrive at the following approximate formulas for the quantities @xmath298 and @xmath299 @xmath300 -\exp\left[-2\pi aq_{-}(\alpha ) \right ]   \right\ }    \label{eq : ap13 }    \\ & & + \frac{1}{2q_{+}(\alpha ) } \left\{\exp\left[-2\pi q_{+}(\alpha ) \right ] -\exp\left[-2\pi aq_{+}(\alpha ) \right ]   \right\ }   -\ln a ,    \nonumber   \\ & & g_{2}(a,0 ) \simeq g_{2}(1,0)+\pi\cos\left[\pi q(\alpha ) \right ] \left\{e_{1}\left[2\pi q_{-}(\alpha ) \right ] -e_{1}\left[2\pi q_{+}(\alpha ) \right ]   \right .
\label{eq : ap14 }   \\ & & \left .
+ e_{1}\left[2\pi aq_{+}(\alpha ) \right ] -e_{1}\left[2\pi aq_{-}(\alpha ) \right ] -2p(\alpha ) \ln a\right\ } .
\nonumber\end{aligned}\ ] ]    finally , for the quantities @xmath301 and @xmath302 it follows from eqs .   and that @xcite @xmath303 ,   \label{eq : ap15 } \\ & & g_{2}(1,0 ) \simeq c_{1 } \sin ( 2\alpha ) -\frac{1}{2 } c_{2}\sin ( 4\alpha ) + \frac{2}{3}\pi ^{3}\cos \left [ \pi q(\alpha ) \right ] p(\alpha ) \left [ p^{2}(\alpha ) -\frac{1}{4}\right ] ,   \label{eq : ap16}\end{aligned}\]]where @xmath304 , @xmath305 . for approximate evaluation of the quantities @xmath301 and @xmath302
we have used the expansion @xmath306 at @xmath307 for the logarithmic function
. then at @xmath289 the final analytical expressions for @xmath284 and @xmath285 are obtained from eqs .  , , and , , , respectively .
the relative accuracy of the derived approximations is less than @xmath308 in a wide range of the parameters
.    99 akopyan , e.a . ,
nersisyan , h.b . & matevosyan , h.h .
energy losses of a charged particle in a plasma in an external field allowing for the field action on plasma and particle motion .
_ radiophys .
_ * 40 * , 823826 .
couillaud , c. , deicas , r. , nardin , ph . ,
beuve , m.a . ,
guihaume , j.m . , renaud , m. , cukier , m. , deutsch , c. & maynard , g. ( 1994 ) .
ionization and stopping of heavy ions in dense laser  ablated plasmas .
e _ * 49 * , 15451562 .
frank , a. , blaevi , a. , grande , p.l . , harres , k. , heling , th . , hoffmann , d.h.h .
, knobloch - maas , r. , kuznetsov , p.g . ,
nrnberg , f. , pelka , a. , schaumann , g. , schiwietz , g. , schkel , a. , schollmeier , m. , schumacher , d. , schtrumpf , j. , vatulin , v.v . ,
vinokurov , o.a . &
roth , m. ( 2010 ) .
energy loss of argon in a laser - generated carbon plasma .
e _ * 81 * , 026401 ( 16 ) .
frank , a. , blaevi , a. , bagnoud , v. , basko , m.m .
, brner , m. , cayzak , w. , kraus , d. , heling , th . , hoffmann , d.h.h . ,
ortner , a. , otten , a. , pelka , a. , pepler , d. , schumacher , d. , tauschwitz , an .
& roth , m. ( 2013 ) .
energy loss and charge transfer of argon in a laser - generated carbon plasma .
lett . _ * 110 * , 115001 ( 15 ) .
hoffmann , d.h.h . ,
tahir , n.a . ,
udrea , s. , rosmej , o. , meister , c.v . , varentsov , d. , roth , m. , schaumann , g. , frank , a. , blaevi , a. , ling , j. , hug , a. , menzel , j. , hessling , th . , harres , k. , gnther , m. , el - moussati , s. , schumacher , d. & imran , m. ( 2010 ) .
high energy density physics with heavy ion beams and related interaction phenomena .
_ contrib .
plasma phys . _
* 50 * , 715 .
oguri , y. , tsubuku , k. , sakumi , a. , shibata , k. , sato , r. , nishigori , k. , hasegawa , j. & ogawa , m. ( 2000 ) .
heavy ion stripping by a highly - ionized laser plasma .
methods phys .
b _ * 161 - 163 * , 155158 .
roth , m. , cowan , t.e . , key , m.h . ,
hatchett , s.p . ,
brown , c. , fountain , w. , johnson , j. , pennington , d.m . , snavely , r.a . ,
wilks , s.c .
, yasuike , k. , ruhl , h. , pegoraro , f. , bulanov , s.v . , campbell , e.m . ,
perry , m.d .
& powell , h. ( 2001 ) .
fast ignition by intense laser  accelerated proton beams .
lett . _ * 86 * , 436439 .
roth , m. , stckll , c. , s , w. , iwase , o. , gericke , d.o . ,
bock , r. , hoffmann , d.h.h . ,
geissel , m. & seelig , w. ( 2000 ) .
energy loss of heavy ions in laser - produced plasmas .
_ europhys .
lett . _ * 50 * , 2834 .
stckl , c. , frankenheim , o.b . ,
roth , m. , su , w. , wetzler , h. , seelig , w. , kulish , m. , dornik , m. , laux , w. , spiller , p. , stetter , m. , stwe , s. , jacoby , j. & hoffmann , d.h.h . ( 1996 ) .
interaction of heavy ion beams with dense plasmas .
_ laser part .
beams _ * 14 * , 561574 .
tabak , m. , hammer , j. , glinsky , m.e .
, kruer , w.l . , wilks , s.c . , woodworth , j. , campbell , e.m .
, perry , m.d . & mason , r.j .
ignition and high gain with ultrapowerful lasers .
plasmas _ * 1 * , 16261634 .
tahir , n.a . ,
deutsch , c. , fortov , v.e .
, gryaznov , v. , hoffmann , d.h.h . , kulish , m. , lomonosov , i.v . ,
mintsev , v. , ni , p. , nikolaev , d. , piriz , a.r . , shilkin , n. , spiller , p. , shutov , a. , temporal , m. , ternovoi , v. , udrea , s. & varentsov , d. ( 2005 ) .
proposal for the study of thermophysical properties of high - energy - density matter using current and future heavy - ion accelerator facilities at gsi darmstadt .
* 95 * , 035001 ( 14 ) .
tavdgiridze , t.l . & tsintsadze , n.l .
energy losses by a charged particle in an isotropic plasma located in an external high frequency electric field .
_ zh . eksp .
* 58 * , 975978 [ english translation : _ sov .
. jetp _ * 31 * , 524525 ( 1970 ) ] .
wang , g .- q . , song , y .- h . ,
wang , y .- n . & mikovi , z.l . ( 2002 ) .
influence of a laser field on coulomb explosions and stopping power for swift molecular ions interacting with solids .
a _ * 66 * , 042901 ( 111 ) .
wang , g .- q . , e , p. , wang , y .-
hu , z .- h . , gao , h. , wang , y .- c . , yao , l. , zhong , h .- y . , cheng , l .- h . , yang , k. , liu , w. & xu , d .-
influence of a strong laser field on coulomb explosion and stopping power of energetic h@xmath3 clusters in plasmas .
plasmas _ * 19 * , 093116 ( 15 ) .


the magnitude of finite metric spaces is a partially defined numerical invariant which was introduced by leinster via a category theoretic analogue of the euler characteristic of spaces , group and posets .
if a finite metric space is scaled up sufficiently then the magnitude tends upwards to the number of points in the metric space and , typically , as the metric space is scaled down the magnitude tends to one ; this allows one to think of the magnitude as the ` effective number of points ' in the metric space .
interestingly , magnitude has applications in measuring biodiversity ( see @xcite and @xcite ) .    in previous papers @xcite and @xcite we speculatively studied the magnitude of _ infinite _ metric spaces using finite subsets as approximations .
whilst we do not know that the magnitude can always be defined in that way it seems to give meaningful answers for many subsets of euclidean space , both numerically and theoretically .
for instance , in  @xcite it is shown using that approach that a closed , bounded subinterval of the real line can be assigned a magnitude equal to one plus half its length ; similarly magnitudes can be assigned to circles and to cantor sets .
( these results are corroborated here in section  1 using the measure theoretic approach of this paper ) . in  @xcite
the magnitudes of various simple subsets of @xmath1 and @xmath2 were calculated numerically and for various spaces considered it was observed that if the metric space was scaled up then the magnitude seemed to asymptotically depend on classical intrinsic volumes of the space such as volume and euler characteristic . for _ convex
_ sets it appeared that this dependence held even non - asymptotically , so for a square it seems , numerically at least , that the magnitude is given by @xmath3 .    in this paper
we look at the magnitude of homogeneous riemannian manifolds , using a measure theoretic approach rather than the finite subset approach mentioned above .
the magnitude is calculated explicitly for the @xmath0-sphere of radius @xmath4 with its intrinsic metric , leading to the following answer .
r \bigl(\bigl(\frac{r}{n-1}\bigr)^2 +   1\bigr))\bigl(\bigl(\frac{r}{n-3}\bigr)^2 + 1\bigr ) \dots \bigl(\bigl(\frac{r}{2}\bigr)^2 + 1\bigr ) } { 1-e^{-\pi r}}\quad &   \text{for $ n$ odd}. \end{cases}\ ] ] for a general homogeneous riemannian manifold @xmath6 we calculate the leading three terms in the magnitude as the manifold is scaled up : the first term is proportional to the volume , the second term vanishes and the third term is proportional to the total scalar curvature ( @xmath7 ) .
writing @xmath8 for @xmath6 scaled up by a factor of @xmath9 and writing @xmath10 for the volume of the unit @xmath0-ball , we have @xmath11 this suffices to determine the asymptotics of the magnitude of a homogeneous surface , this asymptotic expression involving just the area and euler characteristic of the surface @xmath12 : @xmath13 this expression in the area and euler characteristic is precisely the penguin valuation of @xcite .
there are a few things to note .
the first thing to note is that although we use the measure theoretic approach rather than the finite approximation approach , we show that all the results of @xcite obtained using the finite approximations are all also obtained using the measure theoretic approach .
the second thing to note is that in the cases that are known , namely the @xmath0-spheres and surfaces , the constant term is precisely the euler characteristic .
this is particularly interesting as the magnitude arose from a category theoretic generalization of the notion of euler characteristic of various mathematical objects .
indeed we might have named magnitude as `` the euler chacteristic of metric spaces '' if that were nt such a potentially confusing term .
the final thing to note is that the intrinsic volumes , which are the terms so far identified in the asymptotic expansion all appear in the asymptotic expansion of the heat kernel @xcite and , as the definition of is slightly reminiscent of the heat kernel , this looks intriguing . before giving a more detailed synopsis of the contents of the paper ,
it is worth making some comments about the metric that we are considering on a riemannian manifold .
there is potential for confusion here in that we are considering metric spaces and riemannian metrics , so the word ` metric ' is being used with two different senses . to reduce the likelihood of confusion , in this paragraph and the next paragraph we will talk about
` distance metric ' for the metric on a metric space and ` riemannian metric ' for an inner product on each tangent space of a smooth manifold : a distance metric is a global object , a riemannian
metric is an infinitessimal object . in order to define the magnitude of a riemannian manifold
we need a distance metric on the manifold , we will be considering the intrinsic distance metric which arises in the following way .
a riemannian metric gives rise to an arc - length function and this gives rise to a distance metric by defining the distance between two points to be the infimum of the lengths of the paths between the two points .
if there is a riemannian embedding of a riemannian manifold in euclidean space then there are two obvious distance metrics , the intrinsic metric and the subspace metric : in this paper we are primarily considering the intrinsic metric whereas in @xcite and @xcite we considered mainly the subspace metric .
so , for instance , on the earth , the intrinsic distance metric is realized by travelling over the surface of the earth , but the subspace distance metric is realized by tunnelling _ through _ the earth .
we compare the asymptotics for the magnitude of the @xmath0-sphere with the subspace metric and to that with the intrinsic metric ; they are seen to have the same leading term but different lower order terms .    there
now follows a more detailed synopsis of what is in this paper .      in section  1
we define the magnitude of a metric space using a measure on the space .
the key observation here is that if @xmath6 is a homogeneous metric space and @xmath14 is a non - zero homogeneous measure on @xmath6 , then the magnitude of @xmath6 is given by the following formula : for any @xmath15 , @xmath16 provided that the denominator is finite and non - zero .
it is also shown in this section that the measure theoretic definition gives the same magnitude as the finite subset approach did for the circle , line segment and cantor set in  @xcite .    in section
2 we give a brief exposition on the notion of intrinsic volumes of riemannian manifolds .
intrinsic volumes are defined for a wide class of subsets of euclidean space and include such invariants as the volume , the surface area of the boundary and the euler characteristic . in @xcite and @xcite we mainly considered these for polyconvex sets  that is unions of convex sets  following klain and rota s treatment @xcite , but here we consider a differential treatment , utilizing @xcite and @xcite .
one key point for us is that for @xmath17 a closed submanifold of @xmath18 these intrinsic volumes are intrinsic in the sense of differential geometry  that is to say they are independent of the embedding and depend only on the riemannian metric of @xmath17 .
so for such a riemannian manifold there are @xmath19 intrinsic volumes @xmath20 , but it turns out that about half of these vanish because @xmath17 is closed , @xmath21 if @xmath22 is odd : the fact that @xmath23 can be explained by the fact that @xmath24 measures the surface of the boundary , which is of course trivial when @xmath17 is closed .
the non - zero intrinsic volumes are also known in this case as killing - lipschitz curvatures and in a few cases are well - known invariants : @xmath25 is the volume of @xmath17 ; the next non - trivial intrinsic volume @xmath26 is the total scalar curvature ( tsc ) of @xmath17 ; and @xmath27 is the euler characteristic of @xmath17 . in general
these can be obtained by integrating over @xmath17 a certain form defined in terms of the curvature .
in section  3 we get to one of the main results which is an explicit formula for @xmath28 the magnitude of the @xmath0-sphere of radius @xmath4 equipped with its intrinsic metric .
( here it is slightly misleading to talk of the radius of the @xmath0-sphere as it is being thought of as an abstract manifold with the arclength distance metric and not as a subset of @xmath29 with the subspace distance metric ; however it is perhaps slightly less comfortable to talk of the sphere of curvature @xmath30 . )
the formula has many of the properties you might expect , such as @xmath31 as @xmath32 .
asymptotically as @xmath33 , writing @xmath10 for the volume of the unit @xmath0-ball , we get @xmath34 this is certainly consistent with the asymptotics being described by the intrinsic volumes : for instance the alternate terms vanish and we pick up the volume , the total scalar curvature and the euler characteristic .    in section  4
we see the other main result that for all homogeneous riemannian manifolds the first three terms in the asymptotic expansion are like those in the above expansion , further work will be needed to get more terms in the expansion .
however , this is sufficient to show that for a homogeneous riemannian surface asymptotically the magnitude has the form conjectured for the asymptotic magnitude of metric subsets of euclidean space in @xcite and @xcite .    in section  5 the magnitude of the @xmath0-sphere with its subspace metric
is considered .
a closed form is given for @xmath35 , but not for general @xmath0 ; however , for general @xmath0 it is shown that the leading term in the asymptotic expansion is the same as for the intrinsic metric , namely @xmath36 , but that the next non - trivial term is different .
it is a pleasure to thank jonathan jordan , tom leinster and david speyer for their contributions to section  1 , as well as for other informative discussions .
i would also like to thank mark meckes for comments on an earlier draft .
in this section we first recall the definition of the magnitude of finite metric spaces in terms of weights and generalize it to infinite metric spaces using weight measures . we show that speyer s theorem for homogeneous spaces immediately generalizes and that this definition is consistent with previous calculations for intervals , circles and cantor sets .
leinster introduced the notion of the euler characteristic of a finite category which is not always defined , but which generalizes many pre - existing notions of euler characteristic , such as the euler characteristic of a group or a poset or the cardinality of a set . then using lawvere s observation that a metric space can be viewed as an enriched category he defined the notion of the euler characteristic or cardinality of a finite metric space
 it was decided that these terms were already too overburdened , so we renamed the concept the magnitude of a metric space .
this has the following definition .    for @xmath6 a finite metric space , a _ weighting _ on @xmath6
is an assignment of a number @xmath37 to each point @xmath38 such that the _ weight equation _ is satisfied for each point @xmath15 : @xmath39 if a weighting exists for @xmath6 then the magnitude @xmath40 is defined to be the sum of the weights : @xmath41    note the following two things about the definition .
firstly , the weights can be negative .
secondly , if more than one weighting exists then the magnitude is independent of the choice of weighting ( a proof appears below in the more general measure theoretic setting ) .
if you require a little intuition then one way of thinking about the weight equations is as follows .
each point can be thought of as an organism , with the weight being the amount of ` heat ' being generated by that point ( a negative weight denoting some sort of refrigeration ) .
each organism experiences the heat generated by the other organisms in a fashion which decays exponentially with the distance
. the organisms each wish to experience unit heat and this is what is expressed by the weight equation .    we can try to extend the idea of magnitude to non - finite metric spaces . two approaches present themselves .
* use finite approximation to the metric space .
* use a measure in the definition of a weighting .
the first of these approaches was explored in  @xcite and  @xcite . in this approach
a metric space @xmath6 is approximated by taking a sequence @xmath42 of finite subsets of @xmath6 , which converges to @xmath6 in the hausdorff topology , so @xmath43 as @xmath44 .
one can then take the corresponding sequence of magnitudes @xmath45 and hope that they tend to a limit , and furthermore hope that the limit is independent of the choice of approximating sequence .
this seems to work in many cases such as compact subsets of euclidean space .
this approach has the advantage that it can be implemented on a computer to obtain approximations to the magnitude of some spaces as in  @xcite .
the second approach is the one that will be used in this paper .
this was independently suggested to me by tom leinster and david speyer ( my own thoughts on the matter lay more in realm of currents and differential forms ) .
the idea is to generalize the definition of magnitude using measures in the following way .
if @xmath6 is a metric space then a _ weight measure _ on @xmath6 is a signed measure @xmath46 on @xmath6 such that for all @xmath15 @xmath47 if a weight measure @xmath46 exists then @xmath40 , the _ magnitude _ of @xmath6 , is defined to be the mass of @xmath6 with respect to @xmath46 : @xmath48    note that the magnitude , if it exists , is independent of the choice of weight measure .
suppose that @xmath46 and @xmath49 are weight measures then @xmath50    it is worth giving a brief mention here to some connections between measure theory and riemannian geometry ( for more details see , for example ,  @xcite ) .
if @xmath17 is an orientable riemannian manifold then there is a volume form @xmath51 , well - defined up to sign  the choice of sign corresponding to the choice of orientation .
the volume form is defined by smoothly picking an ordered orthonormal basis @xmath52 for each tangent space , taking the dual basis @xmath53 and defining @xmath54 .
integrating with respect to this volume form does not quite give a measure because of the sign ambiguity  think how on the real line we have @xmath55 .
however , there is an associated measure @xmath56 , called the density , which obviates the sign problem . in the case of euclidean space
this density is the usual lebesgue measure . as no confusion should occur
, we will not distinguish between the volume form and the associated density .
if @xmath17 is an orientable @xmath0-dimensional submanifold of euclidean space @xmath18 then there are two obvious measures that one might consider .
firstly the riemannian metric on @xmath18 can be pulled back to @xmath17 to make it into an orientable riemannian manifold so we can use the associated density .
alternatively , one could use the @xmath0-dimensional hausdorff measure on @xmath18 restricted to @xmath17 .
provided that the hausdorff measure is normalized so that the @xmath0-dimensional cube @xmath57^n$ ] has unit measure , these two measures on @xmath17 coincide .
one of the key tools of this paper is the following generalization of speyer s theorem on the magnitude of finite homogeneous spaces given in @xcite . recall that a homogeneous space is one for which there is a transitive group of isometries , or , in other words , one in which all of the points look the same .
[ thm : speyer ] suppose that @xmath6 is a homogeneous metric space and @xmath14 is an invariant metric on @xmath6 then @xmath58 is independent of @xmath59 .
if this quantity is non - zero and finite then a weight measure on @xmath6 is given by @xmath60 thus the magnitude is given by @xmath61    this follows immediately from the definition of magnitude .
it is worth pointing out here that every _ compact _ homogeneous space has a unique - up - to - scaling , invariant , positive measure , thus has a well - defined magnitude .    aside from the calculations for the interval and cantor
set given below , this theorem is the only way i currently know for calculating the magnitude of non - finite metric spaces .
we can now look at the spaces whose cardinality was defined using finite approximations in  @xcite , and show that the same results are obtained using weight measures .
these spaces are closed straight line segments , circles and cantor sets .
we consider first @xmath62 the length @xmath63 closed line segment .
the following was essentially told to me by tom leinster and david speyer independently .
[ thm : intervalweightmeasure ] let @xmath14 be the lebesgue measure on @xmath62 the length @xmath63 line segment and let @xmath64 and @xmath65 be the dirac delta measures supported at the two respective end points .
then a weight measure on @xmath62 is given by @xmath66 .
hence the magnitude is simply :
@xmath67    it is just necessary to show that the given measure satisfies the weight equation for every point @xmath68 . we can consider @xmath62 as the subset @xmath69\subset { \mathbb{r}}$ ] where @xmath70 , then for every @xmath68 @xmath71^y_a-[e^{y - x}]_y^b\right )       = 1.\end{aligned}\ ] ] thus the magnitude is given by @xmath72 as required .
now turn to the ternary cantor set .
the following lemma , pointed out to me by jonathan jordan , will be rather useful in calculating the magnitude .
it says that if we have a weight measure on a subset of the line which on some subinterval is a multiple of the lebesgue measure then if the subinterval is removed then a weight measure is obtained for the resulting space by putting simple dirac delta measures at the end points of the removed subinterval .
[ lemma : measureonhole ] if @xmath6 is a subset of @xmath73 with @xmath46 a weight measure on @xmath6 with @xmath69\subset x$ ] and @xmath74}=\frac12\mu|_{[a , b]}$ ] , then a weight measure on @xmath75 is given by @xmath76 so the magnitude satisfies @xmath77    we need the easy - to - verify identity @xmath78 to obtain @xmath79 .
thus for @xmath80 @xmath81\\    & =     \int_{x\in x}e^{-d(x , y)}{\text{\rm d}}\nu-\int_{x\in   [ a , b]}e^{-d(x , y)}\tfrac12{\text{\rm d}}\mu\\     & \qquad+\tfrac12\tanh\bigl(\tfrac{b - a}2\bigr ) \int_{x\in x\setminus ( a , b)}e^{-d(x , y)}({\text{\rm d}}\delta_a+{\text{\rm d}}\delta_b)\\ & = 1+\tfrac12[e^{y - b}-e^{y - a}]+\tfrac12\tanh\bigl(\tfrac{b - a}2\bigr)\bigl(e^{y - a}+e^{y - b}\bigr)\\ & = 1.\end{aligned}\ ] ] the expression for the weight follows immediately
.    this can now be repeatedly to obtain the magnitude of the cantor set , and this agrees with the answer in @xcite .
note that there the cantor set was approximated below via finite subsets , and here it is approximated above by subsets of @xmath82 $ ] with non - zero measure .
the length @xmath63 ternary cantor set @xmath83 has magnitude given by @xmath84    start by defining @xmath85 to be a length @xmath63 line interval which is parametrize by the interval @xmath57 $ ] with the points in the unit interval being identified with their ternary exansion , so that a point a third of the way along @xmath85 has expansion @xmath86 . the interval @xmath85 is given the standard weight measure as in theorem  [ thm : intervalweightmeasure ] , @xmath87 , which will be called @xmath88 .
define @xmath89 and its weight measure @xmath90 inductively by first removing from @xmath91 the open middle third of each maximal subinterval ; this means removing points having expansions @xmath92 of the form @xmath93 the length of each removed open interval will be @xmath94 . by lemma  [ lemma : measureonhole ] above we obtain a weight measure @xmath90 by restricting the weight measure @xmath95 on @xmath91 , except on the newly pruned end points where we have to put a dirac delta measure of mass @xmath96 .    taking the limit , we obtain we obtain a measure @xmath46 on the ternary cantor set @xmath83 , this measure is the restriction of half the lebesgue measure plus a dirac delta on each point in the cantor set which has a finite ternary expansion  the mass at a point with @xmath97 terms in the ternary expansion is @xmath96 , and there are @xmath98 such points , so the total mass of @xmath46 , i.e. , the magnitude is precisely @xmath99 as required .
now finish this section by considering the circle .
the cardinality of @xmath100 the circle of circumference @xmath63 is given by @xmath101    this follows immediately from the homogeneity theorem .    in conclusion ,
the confidence expressed in  @xcite when defining the magnitude in terms of subsets is now justified . in  @xcite
it was unclear that the answers given for the cantor set and circle were independent of the choice of approximating subsets ; but they agree with the answers given here which are independent of the choice of weight measure .
this section is a brief exposition on intrinsic volumes for riemannian manifolds ; these invariants are also known as curvature measures , lipshitz - killing curvatures , generalized curvatures and quermassintegrale .
i will mention weyl s tube formula and the connection with total scalar curvature and volumes of geodesic speres .
this treatment is mainly obtained from the books  @xcite and  @xcite .
this section also serves to fix some notation and conventions .    in  @xcite , following   @xcite , we considered the intrinsic volumes @xmath102 as functions on polyconvex subsets of euclidean space @xmath103 .
in fact the intrinsic volumes are defined for a much larger class of subsets of @xmath103 via the theory of curvature measures ( see @xcite ) .
the intrinsic volumes can be defined for subsets which have an appropriate notion of unit normal bundle  a so - called normal cycle .
this class of subsets includes smooth submanifold with boundary and more generally sets of positive reach . in general
, the basic properties of the intrinsic volumes are the following .
additivity : :    @xmath104 .
homogeneity : :    @xmath105 .
normalization : :    @xmath106^i)=1 $ ] .
then @xmath107 is the usual @xmath108-dimensional volume of @xmath6 and @xmath109 is the euler characteristic of @xmath6 .    for a closed smooth submanifold of @xmath103
there are many formulas for the intrinsic volumes .
the simplest formulas are for when the submanifold in question is a hypersurface @xmath110 . in this case , at each point in @xmath111 there are @xmath0 principal curvatures @xmath112 , and the @xmath113th symmetrized curvatures @xmath114 is defined to be the @xmath113th - elementary symmetric function in the principal curvatures : @xmath115 .
the intrinsic volumes are given , up to factor , by integrating these symmetrized curvatures over the submanifold .
more precisely , recalling that @xmath116 denotes the volume of the unit @xmath113-sphere , for @xmath117 , @xmath118 from this one obtains , for instance , @xmath119 .
this allows us to calculate the intrinsic volumes of @xmath120 the @xmath0-sphere of radius @xmath4 sitting in @xmath121 in the usual way . here
all of the principal curvatures are the reciprocal of the radius , @xmath122 , so the symmetrized curvatures are given by @xmath123 and so , for @xmath117 , @xmath124    in the case of a submanifold of codimension greater than @xmath125 there are not principal curvatures , but the intrinsic volumes can be expressed in terms of the marginally more complicated lipschitz - killing curvatures
. however it is still true that for @xmath126 a closed submanifold @xmath127 whenever @xmath22 is odd . from these curvature expressions , as weyl notes , a good undergraduate could derive the following tube formula . as an easy exercise ,
you may wish to verify this for the @xmath0-sphere using the formulas given above for the intrinsic volumes .
suppose that @xmath128 is a closed , smooth submanifold of euclidean space , that @xmath129 is the set of points of distance at most @xmath130 from @xmath131 , and that @xmath132 is the volume of the unit @xmath97-ball , then @xmath133    weyl s key observation , however , was that the intrinsic volumes of @xmath6 can be expressed purely in terms of its _ intrinsic _ geometry  we have already observed that @xmath134 and @xmath135 .
so the intrinsic volumes only depend on the riemannian metric of @xmath6 and not on the way it is embedded , nor on the codimension in which it is embedded .
for instance , the intrinsic volumes can be written as polynomial riemann curvature tensor . however , in this paper there is just one of these intrinsic volumes that we are interested in .
after the @xmath0-volume , the next non - trivial intrinsic volume is proportional to the _ total scalar curvature _ : @xmath136 here @xmath137 is the _ scalar curvature _ which can be thought of as measuring the `` lack of stuff '' near @xmath138 . for instance , if you draw a circle on the surface of the earth then it has a circumference than a euclidean circle with the same radius  this deficiency in the circumference is measured by the scalar curvature ( in the case of a surface the scalar curvature is twice the gaussian curvature ) . to make this precise , for @xmath38 and @xmath139 ,
let @xmath140 be the set of points at a distance of @xmath141 from @xmath138 ; for @xmath141 sufficiently small this will be an @xmath142-sphere , known as a _
geodesic sphere_. the scalar curvature @xmath143 at @xmath138 is uniquely determined by the second order term in the power series expansion of the volume of geodesic spheres centered at @xmath138 , in the following sense : @xmath144 it is in this way that the scalar curvature can be viewed as a measure of the `` lack of stuff '' near @xmath138 .
the above formula will be used later in the calculation of sub - leading term in the asymptotics of the magnitude of homogeneous metric spaces . as an example , note that in an @xmath0-sphere of radius @xmath4 we have , for all @xmath138 , @xmath145 .
so the scalar curvature is given by @xmath146 which you can check is consistent with the value calculated for @xmath147 above .
this section consists of one of the main results of this paper which is an explicit formula for the magnitude of the @xmath0-sphere .
the asymptotic behaviour is also considered and the euler characteristic is seen to be the constant term in the asymptotic expansion .
we take @xmath148 the @xmath0-sphere with the standard metric of curvature @xmath149 ( which means that it isometrically embeds in @xmath29 with radius @xmath4 ) thought of as a homogeneous space . from speyer
s homogeneous magnitude theorem ( theorem  [ thm : speyer ] ) the magnitude is given by @xmath150 this can be calculated explicitly .
[ thm : magnitudeofspheres ] for @xmath151 , the magnitude of the homogeneous @xmath0-sphere of radius @xmath4 is given by @xmath5 \dfrac{\pi r \bigl(\bigl(\frac{r}{n-1}\bigr)^2 +   1\bigr))\bigl(\bigl(\frac{r}{n-3}\bigr)^2 + 1\bigr ) \dots \bigl(\bigl(\frac{r}{2}\bigr)^2 + 1\bigr ) } { 1-e^{-\pi r}}\quad &   \text{for $ n$ odd } , \end{cases}\ ] ] where for @xmath152 the numerator is interpreted as being @xmath153 .
the idea is to find some explicit integrals to calculate magnitude , evaluate them to obtain @xmath154 and @xmath155 , show that the induction formula holds , @xmath156 and then apply induction .
we first fix some point @xmath157 and define the coordinate @xmath158 on the @xmath0-sphere to be the distance from the point @xmath157 , so @xmath159 .
the points at a distance @xmath158 from @xmath157 form an @xmath142-sphere @xmath160 of euclidean radius @xmath161 , see figure  [ fig : localball ] ;    @xmath162    thus , writing @xmath163 for the volume of the unit @xmath142-sphere , the volume of @xmath160 is given by @xmath164 .
locally , as @xmath158 parametrizes distance , the volume form on @xmath165 is the product of @xmath166 and the volume form on @xmath160 .
so @xmath167 where we made the substitution @xmath168 .
this can be evaluated for the @xmath125-sphere and @xmath169-sphere : @xmath170_0^\pi }       = \frac{\pi r}{1-e^{-\pi r}}\;;\\    \left|s^2_r\right|       =       \frac{\int_{0}^{\pi}\sin(r)\,{\text{\rm d}}r }             { \int_{0}^{\pi } e^{-rr}\sin(r)\ , { \text{\rm d}}r }       = \frac{2}{\text{im}\int_{0}^{\pi } e^{-rr}e^{ir}\ ,
{ \text{\rm d}}r }      = \frac{2(r^2 + 1)}{1+e^{-\pi r}}\;.\end{gathered}\ ] ] these will form the base cases for the induction .
for the inductive step , write @xmath171 so that @xmath172
. we will show that @xmath173 from which it is immediate that @xmath174 , and the required result follows by induction .    to prove the relation for @xmath175 we observe @xmath176^\pi_0+\int_{0}^{\pi}\left(-\sin(r)\right)\left(\tfrac1n \sin^{n}(r)\right){\text{\rm d}}r\\      & = k_n-0-\tfrac1n k_{n+2},\end{aligned}\ ] ] from which we deduce @xmath177 , as we wanted .    to prove the relation for @xmath178 we perform integration by parts twice : @xmath179^\pi_0     -\int_{0}^{\pi }
-\tfrac 1re^{-rr}(n+1)\cos(r)\sin^{n}(r)\ , { \text{\rm d}}r\\     & =     0 +     \tfrac{(n+1)}{r}\left[-\tfrac 1r e^{-rr}\cos(r)\sin^{n}(r)\right]^\pi_0\\       & \qquad -\tfrac{(n+1)}{r}\int_{0}^{\pi } -\tfrac   1re^{-rr}\left(-\sin(r)\sin^{n}(r )       + n\cos(r)\cos(r)\sin^{n-1}(r)\right ) { \text{\rm d}}r\\     & =     0 + 0     + \tfrac{(n+1)}{r^2}\int_{0}^{\pi } e^{-rr}\left(-\sin^{n+1}(r )          + n\left(1-\sin^2(r)\right)\sin^{n-1}(r)\right ) { \text{\rm d}}r\\     & = -\tfrac{(n+1)^2}{r^2}i_{n+2 } + \tfrac{(n+1)n}{r^2}i_n.\end{aligned}\ ] ] from this we deduce that @xmath180 which suffices to complete the inductive step and hence we complete the proof .    the first thing to note from this is that the formula for @xmath181 in the above theorem also works for @xmath182 provided that this notation is interpreted correctly .
if @xmath182 is taken to mean the two points on the equator of @xmath183 the radius @xmath4 circle , then @xmath184 is the metric space consisting of two points a distance @xmath153 apart , and then it is immediate that the magnitude is given by @xmath185    the next thing to note is the @xmath31 as @xmath32 , so effectively , the @xmath0-sphere looks more and more like a point as it gets smaller and smaller .    a third thing to note is that the two formulae can be combined into a single formula which clearly demonstrates the roots of the numerator lying at alternate integral points on the imaginary axis : @xmath186 where @xmath187 and @xmath188 is a constant which ensures that @xmath31 as @xmath32 .
the roots of the numerator are also roots of the denominator , but i have no idea of the relevance of this .
a fourth things worthy of note is that the denominators very quickly converge to @xmath125 .
this means that the behaviour for @xmath4 away from zero , say for @xmath189 , is dominated by the numerator .
we can define @xmath190 , for @xmath191 , to be the numerator : @xmath192 where the empty product is interpreted as @xmath125 for @xmath193 and @xmath152 .
this polynomial in @xmath194 gives us the asymptotic behaviour of @xmath181 in the sense that @xmath195 as @xmath196 .
previously for the closure of an open subset of euclidean space the asymptotics of the magnitude were conjectured to be controlled by certain intrinsic volumes . in the case of the spheres
the leading term and constant term can be seen to behave in the same way .
the leading term of @xmath197 is proportional to the volume of the @xmath0-sphere @xmath165 and the constant term is equal to its euler characteristic : writing @xmath10 for the volume of the unit @xmath0-ball , @xmath198    to deal with the euler characteristic first , note that the euler characteristic of a sphere just depends on the parity of the dimension : @xmath199 and @xmath200 . the constant term in the polynomial @xmath190 is immediately seen to have exactly this property .    that the leading term has the required form can be proved inductively using the inductive formula @xmath201 together with the following relations ( see , for example ,  @xcite ) for @xmath10 the volume of the unit @xmath0-ball and @xmath202 the volume of the unit @xmath0-sphere : @xmath203 assuming that the result is true for @xmath204 and writing @xmath205 for the top coefficient of a polynomial , we have @xmath206 the result follows by induction after checking the base cases , @xmath152 and @xmath35 , using @xmath207 , @xmath208 , @xmath209 and @xmath210 .
we will see in the next section that the leading term in the asymptotic expansion of the magnitude of any homogeneous is of this form . in  @xcite and  @xcite
we gave a precise conjecture for the asymptotic behaviour of the magnitude of ` nice ' subsets of euclidean space and this involved the being asymptotically the same as the ` penguin valuation ' defined for @xmath6 a subset of euclidean space as @xmath211 we are here considering the sphere with its intrinsic metric so it is not a metric subspace of euclidean space , but the leading term and constant term we obtain are the same as for in the penguin valuation .
however , the subdominant term is different
. we would expect @xmath212 but we get the following : it can be proved inductively as the formula for the leading term was in the above theorem , it is also a consequence of the theorem in the next section on general homogeneous riemannian manifolds .
the first non - trivial , subdominant term of @xmath197 is given by @xmath213    these results are compared with analogous ones for the magnitude of the @xmath0-sphere with the subspace metric in section  5 .
in this section we will look at the three leading terms in the asymptotics of the magnitude of a homogeneous riemannian manifold , as the manifold is scaled up .
the leading term will be seen to be proportional to the volume of the manifold , the second - leading term will be seen to be zero and the third - leading term will be seen to be proportional to its total scalar curvature , or in other words the top two leading terms are proportional to the top two intrinsic volumes . in particular , for a homogeneous surface @xmath214 , the magnitude is asymptotically @xmath215 , and this is precisely the evaluation of the functional from @xcite .    by speyer s homogeneous magnitude theorem ( theorem  [ thm : speyer ] ) we know that the magnitude of a homogeneous riemannian manifold @xmath6 is given by @xmath216 if we look at this as we scale up the space we should consider @xmath217 we wish to consider the asymptotics of the denominator @xmath218 . the relevant piece of asymptotic analysis we will require is a version of watson s lemma ( see for example  @xcite ) .    if @xmath219\to [ 0,\infty)$ ] is a bounded function which , for some @xmath220 , has an expansion @xmath221 then @xmath222    the idea behind watson s lemma is that as a function of @xmath223 $ ] , when @xmath224 gets larger and larger @xmath225 looks more and more like a function supported on a smaller and smaller neighbourhood of @xmath226 , so the integral will , up to exponentially small corrections , only depend on the germ of @xmath227 near @xmath226 .
this can now be used to identify the leading asymptotics of the magnitude for homogeneous manifolds .
[ thm : homogeneousasymptotics ] if @xmath6 is an @xmath0-dimensional homogeneous riemannian manifold and @xmath7 denotes total scalar curvature then as @xmath6 is scaled up the asymptotics of the magnitude are as follows : @xmath228    the idea is to pick a point @xmath229 and define the coordinate @xmath230 .
let @xmath231 be the set of points at a distance @xmath141 from @xmath157 , then if @xmath232 denotes the scalar curvature of @xmath6 , which will be the same at every point , then we know @xmath233 because @xmath141 parametrizes distance , locally the volume form splits as @xmath234 .
we can write @xmath235 for the diameter of @xmath6
. then we find , using watson s lemma appropriately , @xmath236 where in the last equality we have used the fact that @xmath237 .
thus we get @xmath238 this gives the first description of the asymptotics , to get the second it is just necessary to recall that @xmath239 that @xmath240 and that @xmath241 .
this theorem immediately gives the asymptotics of the magnitude of a homogeneous surface , by taking @xmath35 in the above and noting that @xmath242 is precisely the euler characteristic @xmath243 .
for @xmath12 a homogeneous riemannian surface , the magnitude is asymptotically given in terms of the area and the euler characteristic by @xmath244    it is interesting to note that in the examples we know of , namely for spheres and surfaces , the constant term in the asymptotic expansion of the magnitude is precisely the euler characteristic .
in this section we look at the metric on the @xmath0-sphere obtained by thinking of it as sub - metric space of @xmath29 via the standard embedding .
an explicit formula is given for the magnitude of the @xmath169-sphere with this metric , whereas for @xmath245 just the leading few terms of the asymptotic expansion are given .
this is sufficient to see that the leading term is the same as for the magnitude of the sphere with the intrinsic metric but that the lower order terms are not the same .
in @xcite we showed that equipping the radius @xmath4 _ circle _ with the intrinsic metric and the subspace metric gave rise to different magnitudes , but that asymptotically the magnitudes differed by a vanishingly small piece : @xmath246 with the techniques so far developed in this paper it is not difficult to calculate the magnitude of the @xmath169 sphere with these two metrics and see that only the leading term in the asymptotic expansion are the same .    for the @xmath169-sphere of radius @xmath4 with
the subspace metric the magnitude is given by @xmath247 thus the sub - leading terms of the asymptotic behaviour are different for the two different metrics .
the intrinsic metric result is in theorem  [ thm : magnitudeofspheres ] , so it is just necessary to do the calculation for the subspace metric .
this is analogous to the calculation for the intrinsic metric . by the homogeneous magnitude theorem @xmath248
using the substitution @xmath249 gives the following evaluation of the integral : @xmath250 the result follows immediately .
whilst i do nt know how to get a general closed form for the magnitude of the sphere with the subspace metric , using the techniques of the previous section it is straight forward to calculate the first few terms in the asymptotic expansion of the magnitude and compare them with those for the intrinsic metric .
for the @xmath0-sphere of radius @xmath4 with the subspace metric the magnitude is given asymptotically by @xmath251 for the @xmath0-sphere of radius @xmath4 with the intrinsic metric the magnitude is given asymptotically by @xmath252    the intrinsic metric result follows from theorem  [ thm : homogeneousasymptotics ] together with the fact from section  [ section : intrinsicvolumes ] that the scalar curvature of the @xmath0-sphere with radius @xmath4 is @xmath253 .
for the subspace metric result , proceed as usual by using speyer s homogeneous magnitude theorem and writing @xmath202 for the volume of the unit @xmath0-sphere .
@xmath254 again making the substitution @xmath249 whilst noting that @xmath255 and @xmath256 the integral can be manipulated as follows : @xmath257 substituting this back into the above formula gives the required result on recalling that @xmath258 .
one can think of this result as saying that just as the intrinsic and subspace metrics are infinitesimally the same in a certain sense  though i do nt know the correct way to formalize that  so the corresponding magnitudes are asymptotically the same in a certain sense , namely @xmath259 as @xmath33 , although they are not the same in a stronger sense because , in general , @xmath260 as @xmath33 .
h.  donnelly , _
heat equation and the volume of tubes _ , inventiones mathematicae ( 1975 ) * 29 * ( 3 ) pp .  239 - 243 .
a.  gray , tubes , birkauser ( 2004 ) .
d.  klain and g .- c .
rota , introduction to geometric probability , cambridge university press ( 1997 ) .
t.  leinster , _ a maximum entropy theorem with applications to the measurement of biodiversity _ , preprint http://arxiv.org/abs/0910.0906[arxiv:0910.0906 ]
. t.  leinster , _ an adventure in analysis _
, blog post at the @xmath0-category caf , november 2009 , http://golem.ph.utexas.edu/category/2009/11/an_adventure_in_analysis.html t.  leinster and s.  willerton , _ on the asymptotic magnitude of metric spaces _ , preprint http://arxiv.org/abs/0908.1582[arxiv:0908.1582 ] .
p.  d.  miller , _ applied asymptotic analysis _ , graduate studies in mathematics 75 , american mathematical society , 2006 .
morvan , generalized curvatures , springer - verlag ( 2008 ) .
s.  schanuel , _ what is the length of a potato ?
an introduction to geometric measure theory _ , in categories in continuum physics ( buffalo , n.y . , 1982 ) , 118126 , lecture notes in mathematics , * 1174 * , springer , berlin , 1986 .
a.  solow and s.  polasky , _ measuring biological diversity _ , environmental and ecological statistics ( 1994 ) * 1 * ( 2 ) 95103 .
s.  willerton , _ heuristic and computer calculations of the magnitude of metric spaces _ , preprint http://arxiv.org/abs/0910.5500[arxiv:0910.5500 ] .
a.  weyl , _ on the volume of tubes _ , american journal of mathematics ( 1939 ) * 61 * ( 2 ) pp .  461472