modern cosmology faces two outstanding challenges , namely the existence and nature of dark energy and dark matter .
many theories have been put forward to address both issues : for dark matter , several candidates are available , such as weak - interacting particles ( wimps ) arising from extensions to the standard model ( _ e.g. _ axions , neutralinos ) , _ etc . _ ; for dark energy , `` quintessence '' models consider the slow - roll of a scalar field @xcite , amongst other candidates ; others suggest that the averaging of inhomogeneities at a cosmological scale may yield an effective scalar field , thus accounting for the dark energy component of the universe @xcite .
a possible unification of both `` dark '' components has also been suggested , resorting to a scalar field model @xcite or an exotic equation of state , as featured by the so - called modified chaplygin gas @xcite .
a different approach assumes that no extra energy content is needed and that the fundamental laws and tenets of gravitation may be incomplete , perhaps just a low - energy approximation ; as a consequence , modifications of the friedmann equation to include higher order terms in the energy density @xmath0 ( see _ e.g. _ @xcite and references therein ) have been proposed or , at a more fundamental level , changes to the action functional .
a rather straight forward approach lies in replacing the linear scalar curvature term in the einstein - hilbert action by a function of the scalar curvature , @xmath1 ; alternatively , one could resort to other scalar invariants of the theory @xcite .
this has led to some success in replicating the accelerated expansion of the universe , while comparison with its evolution throughout the different ages has yielded some constraints and exclusions to the form of @xmath1 ; perhaps the most well - known proposal of this type is the starobinsky inflationary model , where a quadratic term in the curvature is added to the usual linear form ( plus cosmological constant ) , @xmath2 @xcite .
solar system tests could also bring further insight , mostly arising from the parameterized post - newtonian ( ppn ) metric coefficients derived from this extension of general relativity ( gr ) .
however , some disagreement exists in the community , with some arguing that no changes are predicted at a post - newtonian level ( see _ e_.g . @xcite and references therein ) ; amongst other considerations , this mostly stems from an approach based either in the more usual metric affine connection ( that is , where the affine connection is taken _ a priori _ as depending on the metric ) , or in the so - called palatini approach @xcite ( where both the metric and the affine connection are taken as independent variables ) . as an example of a clear phenomenological consequence of this extension of gr , it has been shown that @xmath3 theories yield a gravitational potential which displays an increasing , repulsive contribution , added to the newtonian term @xcite .
notwithstanding the significant literature on these @xmath1 models , few steps have been taken to address another interesting possibility : not only that the curvature is non - trivial in the einstein - hilbert lagrangian , but also that the coupling between matter and geometry is not minimum ; indeed , these are only implicitly related in the action functional , since one expects that covariantly invariant terms in @xmath4 should be constructed by contraction with the metric ( _ e.g. _ the kinetic term of a real scalar field , @xmath5 ) .
a non - minimum coupling would imply that geometric quantities ( such as the scalar invariants ) would explicitly show in the action ; asides from theoretical elegance , this could have deep phenomenological implications : indeed , in regions where the curvature is high ( which , in gr , are related to regions of high energy density or pressure ) , the implications of such theory could deviate considerably from those predicted by einstein s theory @xcite .
related proposals have been put forward previously to address the problem of the accelerated expansion of the universe @xcite and the existence of a cosmological constant @xcite .    in this sense ,
the immediate question , posed in the title of this work , is : what are the implications for the behaviour of matter under such conditions ?
some work has been put forward concerning this issue , namely changes to geodetic behaviour @xcite , the possibility of modelling dark matter @xcite and the violation of the highly constrained equivalence principle @xcite .
perhaps the most important consequence of the mentioned studies is that energy may no longer be covariantly conserved , that is , @xmath6 , where @xmath7 is the energy - momentum tensor of matter ; this occurs because , due to the presence of extra terms in the equations of motion , the bianchi identities no longer imply in the covariant conservation of the energy - momentum tensor .
this work addresses what we believe is the natural proving ground for a non - minimally coupled gravity model : regions where the density may be high enough , to evidence some deviation from gr , although moderate enough so that effects are still perturbative  a star .
the results rely upon and expand the methodology followed by the authors in previous studies @xcite .
this paper is divided in the following sections : first , we present the model upon which the subsequent work is based ; then , we develop the equations of motion , aiming at the modified tolman - oppenheimer - volkoff ( tov ) equation , with due care taken regarding the validity of the newtonian regime and resulting modified hydrostatic equilibrium equation ; afterwards , we insert the polytropic equation of state into the latter , and compute the necessary observables ; a numerical session follows , where profiles and bounds are computed for the relevant quantities ; finally , a discussion of our results is presented .
following the discussion of the previous section , one postulates the following action for the theory @xcite ,    s = d^4 x [ action ]    where @xmath8 ( with @xmath9 ) are arbitrary functions of the scalar curvature @xmath10 , @xmath4 is the lagrangian density of matter and @xmath11 is the metric determinant . for convenience
, the contribution of the non - minimum coupling of @xmath12 is gauged through the coupling constant @xmath13 ( which has dimensions @xmath14 = [ f_2]^{-1}$ ] ) .
the standard einstein - hilbert action is recovered by taking @xmath15 and @xmath16 , where @xmath17 and @xmath18 is the cosmological constant ( from now on , one works in a unit system where @xmath19 ) .
variation with respect to the metric @xmath20 yields the modified einstein equations of motion , here arranged as    & & [ ee ] ( f_1 + 2 f_2 _ m ) r_- 1 2 f_1 g_= + & & ( _ - g _ ) ( f_1 + 2 f_2 _ m ) + t _  ,    where one defines @xmath21 for convenience , as well as @xmath22 , and omitted the argument .
the matter energy - momentum tensor is , as usually , defined by    t_= -2 ( _ m ) g^  .
as stated before , the bianchi identities , @xmath23 imply the non-(covariant ) conservation law    ^t_= f_2 1 + f_2 ( g__m - t _ ) ^r  , [ non - cons ]    and , as expected , in the gr limit @xmath24 , one recovers the conservation law @xmath25 .
it is our stated purpose to arrive at the modified form of the tov equation , the relativistic version of the hydrostatic equilibrium condition . from eq .
( [ ee ] ) , it is clear that a full treatment of the equations of motion is unattainable unless some specific form for @xmath26 and @xmath27 is provided .
furthermore , the presence of both the pure curvature and non - minimum coupling , respectively , still constitutes a daunting analytical challenge ; hence , and since one is mainly interested in the relevance of the effects within a high curvature and pressure medium , where @xmath12 should overwhelm the modification of the pure curvature term , @xmath28 ( neglecting the contribution of the cosmological constant ) , it appears sensible to discard the latter ; a treatment of the standard @xmath1 scenario with @xmath15 may be found in ref .
@xcite .
mathematically , the chosen approximation reads as    [ ineq1 ] 1 2| f_1 - r| & & |f_2 _
m |  , +        where the second and third inequalities indicate that this regime stems not just from the comparison between contributions to the action functional , but also those involved in the modified einstein equations ( [ ee ] ) ; also , notice that the perturbative condition @xmath29 has not been enforced yet .
the last inequality is satisfied if the following stronger conditions hold ,    & & | d f_1 dr | 2| |  , + & & | d^2 f_1 dr^2 |2 | |  .
this said , one considers the simplest form @xmath30 ; this linear coupling could arise from a first order expansion of a more general @xmath27 function , in the weak field environ of the sun .
also , one takes @xmath31 , a natural choice for the lagrangian density of an ideal fluid @xcite  characterised by the standard energy - momentum tensor    t_= ( + p ) u_u_- p g _  , [ tmn ]    where @xmath32 is the fluid s pressure , @xmath0 its matter - energy density and @xmath33 its 4-velocity vector , with @xmath34 , @xmath35 , so that @xmath36 . inserting the above expressions for @xmath12 and @xmath4 , the inequalities ( [ ineq1 ] ) become    [ ineqa ] |f_1 | & & 2 |(+ p ) r |  , + [ ineqb ] | f_1 | & & 2 | + p |  , + [ ineqc ] | d f_1 dr | & & 2 | p(r ) |  ,  | d^2 f_1 dr^2| 2 | p(r ) |  ,    where the prime denotes derivation with respect to the radial coordinate .
clearly , if @xmath37 , these are trivially satisfied . also , notice that @xmath38 implies that @xmath14 = [ f_2^{-1}]=[r]^{-1 } = m^{-2}$ ] .
one must now ascertain the form of @xmath26 . in ref .
@xcite it is shown that acceptable models with @xmath39 ( @xmath40 , @xmath41 = [ r^m ] = m^{2m}}$ ] and @xmath42 or @xmath43 ( with @xmath44 and @xmath41 = [ r^{-1}]= m^{-2}}$ ] ) are cosmologically viable ; the latter may arise from the renomalizability of the theory near the planck scale , with @xmath45 and @xmath46 @xcite . for the perturbative regime to be valid , one must have , for the @xmath47 case    & & [ ineqf1a ] || |1 r^m p|  , + & & || ( 1-m ) m | 1 r^1+mpd r dr |  , + & & || ( 1-m)m |1 r^1+mp
| | d^2 rdr^2 -(1+m)(d r dr)^2r |  .    while , for the @xmath48 case ( which can be derived from the above inequalities , setting @xmath49 , @xmath50 and ignoring the cosmological constant term ) ,    & & [ ineqf1b ] || |r p|  , + & & || 2 | 1 pd r dr |  , + & & || 2 |1 p
 | | d^2 rdr^2 |
.      with the above choice for @xmath51 and @xmath12 , the equation of motion becomes    & & ( 1 + a p ) r_- 1 2r ( g_+a t _ ) = + & & a ( _ - g_)p + 1 2 t _  ,    where one defines the parameter @xmath52 , with dimension @xmath53 = m^{-4}$ ] ; accordingly , both @xmath54 and @xmath55 are dimensionless quantities . in the above form ,
the physical meaning of the proposed model is more transparent : aside from a pressure - dependent term on the _ r.h.s .
_ , the most interesting modification occurs on the _ l.h.s . _ ; firstly , the contribution of the riemann tensor is modified by a factor @xmath56 ; secondly , the scalar curvature is coupled not only to the metric @xmath20 , but also to the energy - momentum tensor @xmath7 .
one thus gets a clear picture of the matter - geometry interaction , which occurs via scalar - tensor combinations .    by taking the trace of the above equation ,
one obtains    r = - -3a p + t = 3p - + 3 a p  , [ ricci ]    having inserted @xmath57 .
interestingly enough , in the `` strong '' @xmath58 regime it yields the asymptotic result @xmath59 : in this regime , a varying , low density can give origin to extremely high curvatures , while an almost uniform , even if a high density fluid might yield a vanishingly small curvature .
foreseeing a later comparison with solar observables , which one knows are well predicted by gr , it is now assumed that the effect of @xmath12 yields only perturbative corrections : in the approximation @xmath60 , the scalar curvature is given by @xmath61 , resulting in @xmath62 ; anticipating the newtonian approximation @xmath63 , this yields @xmath64 and @xmath65 .
inserting eq .
( [ ricci ] ) in the einstein equation , one obtains , after some algebraic manipulation ,    [ motion ] & & r_= + & & ( 3p - ) g_+2 ( 1-ap ) t_+ a ( 4 _ - g _ ) p  ,    keeping only first order terms in @xmath66
. this shall be the main tool for the derivation of the tov equation .
since one is dealing with an ideal , spherically symmetric system , in which temporal variations are assumed to occur only at the cosmological scale @xmath67 , and hence negligible at an astrophysical time scale , one considers the birkhoff metric ( in its anisotropic form ) , given by the line element    ds^2 = e^(r ) dt^2 - ( e^(r)dr^2 + d ^2 )  , [ line ]    with @xmath68 , so that @xmath69 .
one can resort to the expression of the riemann tensor and eq .
( [ motion ] ) to obtain the following intermediary step ( the full algebraic derivation is shown in appendix a ) :    [ eq1 ] & & = - + & & ap2(+3p ) + a 2 ( 5 e^- _ 00 + 3e^- _ rr + 2 _ r^2 ) p  .
one now defines the parameter @xmath70 , here called the effective mass , to distinguish it from its identification with the gravitational mass of the unperturbed gr scenario , derived from the schwarzschild metric ; it is given by the usual expression    e^- = 1 - 2 gm_e r  , [ massdef ]    which yields    = 2gm_e r^2  ,    where the prime denotes differentiation with respect to the radial coordinate .
inserting the above in eq .
( [ eq1 ] ) and solving for @xmath71 , one gets    & & m_e = 4 r^2 + + & & a r^2 4 g ( 5 e^- _ 00 + 3e^- _ rr + 2 r^2 _ ) p 2 + a ( - 3p )  .    in the perturbative regime
, one may take only first - order terms in @xmath54 and @xmath55 , obtaining    [ mass ] & & m_e 4 r^2 + + & & a r^2 8 g ( 5 e^- _ 00
+ 3e^- _ rr + 2 _ r^2 ) p  ,    using @xmath72 .
the above expression clearly shows the perturbation to the purely gravitational mass , defined by @xmath73 .      before continuing the derivation
, it is opportune to address the issue of the validity of the newtonian regime ; clearly , establishing this limit will enable many fruitful simplifications in the calculations ahead .
in the standard derivation of the hydrostatic equilibrium equation , this arises from a set of simplifications imposed on the relativistic tov equation ; the later reads ,    p(r ) = -g r^2 [ ( r ) + p(r ) ] [ m_e(r ) + 4 p(r ) r^3 ] 1 - 2 gm_e(r)/r  , [ tovgr ]    and the newtonian approximation is valid if the following inequalities are satisfied    r 2gm_e(r )  ,   ( r ) p(r )  ,   m_e(r ) 4 p(r ) r^3  , [ newtongr ]    yielding the non - relativistic hydrostatic equation of state    p(r ) = - g ( r ) m_e(r ) r^2  .
[ hydrogr ]    since is is assumed that the coupling between matter and geometry will only produce a perturbative effect , leading to the redefinition of mass through eq .
( [ mass ] ) , it is clear that no changes occur regarding the validity of the newtonian approximation .
hence , foreseeing the application of the following results to the sun , where such regime is valid , one may simplify the intermediate calculations and , when convenient , insert the inequalities ( [ newtongr ] ) in order to simplify @xmath66-dependent terms . by the same token ,
terms involving the coupling between @xmath66-dependent quantities and any covariant derivatives may be evaluated by taking their newtonian counterparts ,    [ dernewt ] & & _ 00 p = -e^- 2 p 0  , + & & _ rrp = p  + 1 2 p p   , + & & _ p = e^-rp rp  .
notice that this does not imply that one is directly deriving the newtonian hydrostatic equilibrium equation , since leading order terms will not be approximated .      by following a procedure similar to the one leading to eq .
( [ eq1 ] ) ( depicted in appendix a ) , one obtains the following equation :    [ eq2 ] & & = + & & p + ap 2 ( p- ) + a 2 ( 3 e^- _ 00 + 5 e^- _ rr - 2 _ r^2 ) p  .    substituting by the expressions for the covariant derivatives , one gets    & & ( 3 e^- _ 00 + 5 e^- _ rr - 2 _ r^2 ) p = + & & 5 ( 1 - 2 g m_e r ) p
+ + & & p  .
one now introduces the approximations discussed in the previous subsection .
also , the approximation @xmath74 is taken , since the perturbative corrections would produce second order terms in @xmath66 .
one obtains , after a little algebra ,    & & ( 2 + a ) (  r -2 g m_e r^3 ) = + & & p ( 1 - a 2 ) + a 2  .    solving for @xmath75 ,
one gets    = & & g ( 2 + a ) m_e + 4 ( 2 - a ) p r^3 r
( r - 2 g m_e ) + + & & a r^2 r - 2gm_e 5 4 p
+ ( 5 g r - 1 2 r ) p 2 + a ( + 3 4 r p )  .
when @xmath76 , one recovers the standard expression    = g m_e + 4 p r^3 r^2 - 2gm_e r  .
linearizing with respect to @xmath66 yields    = g m_e + 4 p r^3 r^2 - 2gm_e r + a h(p , )  ,    where the function @xmath77 is defined through    & & h(p , ) = + & & ( 5 8p  - 4 g p ) r + ( 5 2 g r^2 -3 8 g m_e r - 1 4 ) p + & & ( 5 8p  - 4 g p ) r -p 4  ,    after considering the inequality @xmath78 and also @xmath79 ; taking the sun s maximum , central , density @xmath80 , and the radius @xmath81 , one gets @xmath82 .    before continuing
, one can rewrite the expression for the gravitational mass , obtained before , but imposing the limit @xmath83 , @xmath84 , which yields    _ 00 p 0  ,   _ rr p p   ,   _ p r p  .
this approximation , together with @xmath63 , implies that    m_e 4 r^2 + ar^2 ( 3 p8 g + p 4 gr - 2 ^2 )  .
similarly , one obtains    r & & -- 3 a ( p  + 2 r p ) ( 2 + a )   + & & - 2 - a 2  .    although the expression for the scalar curvature @xmath10 is somewhat involved , one can approximate its derivative by the unperturbed value @xmath85 , losing only terms of order @xmath86 ; indeed , one could write    r & = & - 2 + a f(,p ) + r & = & - 2 + a f(,p ) = - a 2  + a^2 f(,p )  ,    with    f(,p)= - 1 2   .    to write the modified equation of hydrostatic equilibrium , one now resorts to the non-(covariant ) conservation of the energy - momentum tensor , eq .
( [ non - cons ] ) .
the expression for the covariant derivative is purely geometric ; in order to derive an expression for @xmath87 , one aims at the @xmath88 component of the above equation ; with our choice of @xmath38 and @xmath31 , one gets the modified tov equation    & & p +  2 ( + p ) = 2 1 + r ( 1 - 2gm_e r ) r + & & - 2(+ p ) p- 2 p 1 + r r p + a p   ,    dropping higher order terms and using @xmath85 .
thus , one finally obtains a set of three differential equations for the problem at hand :    [ eqset1 ] p & - & 2 p r = - 2 ( + p )  , + [ eqset2 ] m_e & = & 4 r^2 + ar^2 ( 3 p8 g + p 4 gr - 2 ^2 )  , + [ eqset3 ]  2 & = & g m_e + 4 p r^3 r^2 -2 gm_er + a  .    substituting eq .
( [ eqset1 ] ) into ( [ eqset3 ] ) , after some algebra one gets the modified tov equation ,    [ toveq ] & & p + g ( + p ) m_e + 4 p r^3 r^2 - 2gm_er = + & & a  .
this yields the non - relativistic hydrostatic equilibrium equation ,    [ hydroeq ] p + gm_e r^2 = a  .
realistic stellar models rely on four differential equations , together with appropriate definitions @xcite ; aside from the mass conservation condition ( [ eqset2 ] ) and the hydrostatic equilibrium equation ( [ toveq ] ) or ( [ hydroeq ] ) , these express energy conservation and transport , through    = 4 r^2 ( r )    and    = - l(r ) 4r^2 _ c  ,    respectively .
in the above , @xmath89 is the energy flow across a sphere of radius @xmath90 , @xmath91 is the energy generation rate per mass unit , and @xmath92 is the conductivity coefficient . for given @xmath91 and @xmath92 , which account for the processes ongoing inside the star , one
is left with five unknowns , @xmath71 , @xmath93 , @xmath0 , @xmath94 and @xmath95 ; an additional relation is needed , in the form of a suitable equation of state @xmath96 .
many candidate equations of state and solar models are available , with varying degrees of sofistication , accounting for effects such as chemical composition , solar matter mixing , discontinuities between layers , heat difusion , _ etc .
_ ; two outstanding examples are the mihalas - hummer - dappen and opal equations of state @xcite .
however , solving the above set of differential equations with a realistic equation of state requires heavy - duty numerical integration with complex code designs , and an analysis of the perturbations induced by a non - minimum coupling between matter and curvature is beyond the scope of the present study .
instead , one may resort to a very simplistic assumption , the so - called polytropic equation of state .
this is commonly given by @xmath97 , where @xmath98 is the polytropic constant , @xmath99 is the baryonic mass density and @xmath100 is the polytropic index @xcite .
the polytropic index @xmath100 interpolates between the basic thermodynamical processes : @xmath101 for isobaric , @xmath102 for isometric , and an infinite polytropic index @xmath100 for an isothermal sphere .
adiabatic processes yield @xmath103 , with @xmath104 the adiabatic coefficient .
several crude approximations to relevant astrophysical systems are also obtained : @xmath105 may model the degenerate star cores found in giant ( gaseous ) planets , white or brown dwarfs , and red giants ; a polytropic index @xmath106 yields a boundless system ( that is , with non - vanishing density everywhere ) , which was taken by schuster as the first candidate for a stellar system . finally ,
a polytropic equation of state was used by eddington in his proposal for the first solar model , with @xmath107 ; naturally , it does not offer an accurate description of the solar interior , and has been deprecated by the following advancements .    nonetheless , the use of a polytropic equation of state is still of interest , due to its simplicity and ease of manipulation , which render it a valuable tool in more theoretically driven studies , as is the present case ( as an example , the generalized chaplygin gas is a polytrope with index @xmath108 @xcite ) .
clearly , more realistic assumptions regarding the structure of the sun would improve the final results ; also , the procedure outlined in this work could also be applied to more exotic bodies , either through the use of an adequate polytropic index , or via a more realistic equation of state , possibly yielding a more stringent constraint on the coupling parameter @xmath13 .
recall that the effective mass @xmath70 is defined in terms of the energy density @xmath0 , which appears in the energy - momentum tensor ; these two quantities are related through @xmath109 , yielding the relation    = ( p k)^n / ( n+1 ) +
n p  .
however , since one is interested in probing the newtonian regime which occurs in the sun , the condition @xmath63 holds ; therefore , @xmath110 , and one may take the form @xmath111 for the equation of state @xcite .    in order to transform the modified hydrostatic equilibrium eq .
( [ hydroeq ] ) into a differential equation with a single variable , one first rewrites eq .
( [ hydroeq ] ) as    & & 1 r^2 + & & =
-4 g - a ( 3 p8 + p 4 r - 2 g ^2 )  ,    and inserts the polytropic equation of state , @xmath111 , written as @xmath112 and @xmath113 , with @xmath114 a dimensionless variable and @xmath115 ; as stated before , @xmath80 is the central density , and @xmath116 is the central pressure .
one obtains the perturbed lane - emden equation for the function @xmath117 :    [ le ] 1 ^ 2  = - ^n  ,    where the prime now denotes derivation with respect to the dimensionless radial coordinate @xmath118 , and one has defined @xmath119 and @xmath120 , for convenience . obviously , setting @xmath121
one recovers the unperturbed lane - emden equation @xcite    [ le0 ] 1 ^2 ( ^2  )
= - ^n  .
one may try to solve eq .
( [ le ] ) analytically around @xmath122 , and compare with the solution of the unperturbed equation , given ( in the vicinity of @xmath122 ) by @xcite    ( ) 1 - 1 6 ^2 + n 120 ^4  .
[ eq0 ]    this calculation ( outlined in appendix b ) yields    ( ) 1 -a ^2 + b^4  , [ origin ]    with    a = 1 6 ( 1 - a_c 13 + 25 n 12(n+1 ) )  ,    and    b= n 120 ( 1 - a_c 11 60 39 + 59n n+1 )  ,    neglecting the @xmath123 term and taking the limit @xmath124 .
taking @xmath125 gives back the unperturbed values @xmath126 and @xmath127 , as expected .
for @xmath102 ( a constant density model ) , one finds    ( ) 1- 1 6 ( 1 - 11 6 a_c ) ^2  .    for @xmath128 ( the first proposed model for the sun ) , one obtains    ( ) 1- 1 6 ( 1 - 11 6 a_c ) ^2 + 1 40 ( 1 - 9910 a_c ) ^4  .    finally , for @xmath129 ( an isothermal sphere ) , one gets    ( ) 1- 1 6 ( 1 - 23 12 a_c ) ^2 + 1 40 ( 1 - 1837180 a_c ) ^4
.      inspection of eq .
( [ le ] ) shows that it is a third - order equation , while the unperturbed lane - emden equation is only second - order .
thus , the initial conditions @xmath130 and @xmath131 alone are not sufficient .
indeed , an extra condition regarding the initial behaviour of the second derivative must be provided ; this is derived from the series expansion of @xmath132 around @xmath122 taken before , that is     ( 0 ) = -1 3(1 - a_c 13 + 25n 12(n+1 ) )  .    in order to numerically solve the perturbed lane - emden equation ,
one rewrites eq .
( [ le ] ) ( after a long derivation ) as    & &   = - 8 5 a_c - n ( n - 1 ) ^3 ^2 + & & + 24 5 n_c ^n+1 - 2 5 + + & & 4 5 ^n- n    + & &  .    clearly , when @xmath133 , the first term blows up unless it is compensated by the condition @xmath134 , which is precisely the unperturbed lane - emden equation . by the same token
, the first term expresses the deviation from the unperturbed case , and should be of order @xmath135 , therefore cancelling out the divergence .
however , implementing the above third order differential equation proves too computationally demanding ; thus , one must first approach the perturbed lane - emden equation and , given that one is searching for a perturbation , expand it in terms of @xmath136 , where @xmath137 is the solution to the unperturbed lane - emden equation and @xmath138 is the ( relative ) perturbation , obeying @xmath139 .
thus , one may write , for @xmath140 ,    ^n = _ 0^n ( 1 + a_c ) ^n _ 0^n ( 1 + na_c )  .
if this expansion is introduced in the perturbed lane - emden equation , and considering only terms of order @xmath135 , one obtains a differential equation for @xmath138 , with @xmath137 as source :    & &  + 2 ( _ 0 _ 0 + 1 )  + ( n-1 ) _
0^n-1 = 5n2 _ 0 ^ 2 n-2 _ 0 + & & + ( 2 n + 1 ) n_c _ 0 ^ 2 n-1 _ 0 + 9 n+5 4 ( n+1 ) _ 0 ^ 2 n-1 + 3 n_c _ 0 ^ 2n + & & - 5n ( n-1 ) 8 _ 0^n-3 _ 0^3 + n(3n+7 )
4(n+1 ) _ 0^n-2 _ 0^2 + 1 2_0^n-1 _ 0  .
having eliminated the second derivative of the unperturbed solution through eq .
( [ le0 ] ) . from eq .
( [ origin ] ) , one concludes that this differential equation is supplemented by the initial conditions @xmath141 and @xmath142 .    note
that @xmath138 does not depend on @xmath135 : one must only find the unique @xmath138 ( for each @xmath143 and @xmath100 ) ; one aims at a simultaneous variation of the model s parameter @xmath135 and the polytropic index @xmath100 in the vicinity of the standard solar value @xmath128 , enabling the plotting of an exclusion graph in the @xmath144 plane . for @xmath128 ,
the differential equation for @xmath138 simplifies to    & &  + 2 ( _ 0 _ 0 + 1 )  + 2 _ 0 ^ 2 = 3 n_c _ 0 ^ 6 + 2 _ 0 ^ 5 + + & & 7 n_c _ 0 ^ 5 _ 0+ 152 _ 0 ^ 4 _
0 + 1 2_0 ^ 2 _ 0 + 3 _ 0 _
0^2 - 15 4 _ 0^3  .
one now computes the deviation between the effective mass @xmath70 , defined in eq .
( [ mass ] ) , and the gravitational mass , defined by @xmath73 ; some algebra yields    [ massdif ] m_e - m_g & = & n+1 g a_c n_c
^2 _ 0^n  ,    which is explicitly written in terms of @xmath145 , the model s dimensionless parameter , and @xmath146 , the ratio that measures the validity of the newtonian regime .
one has also to eliminate @xmath147 through eq .
( [ le0 ] ) ; one concludes that not only is the above mass difference small ( due to the perturbative regime , @xmath124 ) , but that it is further suppressed by the factor @xmath123 .
clearly , the model s parameters should be adjusted to the known observables : the star s
radius @xmath148 and mass @xmath149 .
the issue of this identification is rather delicate ; indeed , @xmath149 should not be identified with the gravitational mass @xmath150 , but with the total effective mass @xmath151 .
asides from the physical interpretation of the full energy content of the star affecting geodetic motion around it , one can also resort to the matching conditions of the inner metric with the outer schwarzschild metric ; indeed , writing @xmath152 , @xmath153 , @xmath154 and @xmath155 , one has    _ - ( r _ ) & = & -2 g m_e - m_e r_r_(r_-2gm_e )  , + _ +  & = & - 2gm r_(r_-2gm )  , + _ - ( r _ ) & = & 2 g m_e r_(r_-2gm_e ) + + & & a 4 ( 5 2 p(r _ ) r_- p(r _ ) )  , + _ +  & = & - _ +  = 2gm r_(r_-2gm )  ,    where the @xmath156 and @xmath157 subscripts indicate inner or outer boundary condition .
if one identifies @xmath158 , it follows that @xmath159 , which implies    + p ( r_)r _ = 0  ,    ( after resorting to eq . ( [ mass ] )
, with @xmath160 ) , and    p(r _ ) r_-p(r _ ) = 0  .
clearly , both conditions hold , since    p(r _ ) & & (_1 ) ^n(_1 ) = 0  , + p
( r _ ) & &  ( _ 1)(_1)^n + n(_1)^2 ( _ 1)^n-1 = 0  ,    where @xmath161 signals the star s boundary , through @xmath162 .      given the identification @xmath163
, one can now deduce the perturbation to the central density @xmath164 , which is a model dependent parameter .
however , before using eq .
( [ le ] ) to evaluate the perturbation @xmath138 and extract relevant quantities , it should be noticed that , by inspection , it is clear that @xmath138 diverges when @xmath137 approaches zero .
hence , one can not extend the perturbational approach to the full range of the star , and should instead deal with the full differential equation for the perturbed @xmath132 .    in order to circumvent this issue ,
recall that there is a pronounced deviation between the predictions of the polytropic model and the realistic standard solar model , for @xmath165 ; this reflects the crossing from the radiative zone , where the @xmath128 polytrope is a good approximation for the sun , and the convection zone , where this approach fails .
hence , the irregular behaviour of @xmath138 near @xmath161 may be safely disregarded , since the fundamental equation for @xmath132 is not valid there . instead
, one shall consider only the range @xmath166 or , equivalently , @xmath167 .
in doing so , one is of course neglecting the contribution of the latter for all relevant quantities : however , although the density and pressure are still significant at that point , both the polytropic and the standard solar model show that up to @xmath168 of the sun s mass is located within the radiative zone .    also , it can be numerically shown that the boundary condition @xmath169 does not shift significantly from the unperturbed @xmath170 case , so that no problem arises from neglecting any changes to the scale factor @xmath171 ( which , recall , relates the dimensionless coordinate @xmath118 with the physical distance to the center @xmath90 ) .
furthermore , the matching of the inner and outer derivatives of the metric should not be taken as a realistic condition , but merely a consistency check for the developed model .
one proceeds with the calculation of the total mass of the sun or , more accurately , the mass of the radiative zone : for the unperturbed case , one has ( see appendix c for the derivation of the following results )    m = -4 _ c0 _ 0r _ c0 = -m 4 _ 0r _ r r_r^3  .    where one defines @xmath172 and @xmath173 is the unperturbed central density ( that is , obtained from @xmath149 , @xmath10 and the numerical results for the unperturbed solution @xmath137 ) .    in the perturbed case ,
one obtains    & & m = -4 _ c _ 0r + & &  .
notice that there are two @xmath135-dependent contributions : one arising from @xmath137 and its derivatives , since including the perturbation @xmath138 would only amount to second - order terms @xmath174 , and other involving the integral of @xmath138 .    since the total mass @xmath149 does not change ( only its interpretation as a purely gravitational mass or a sum of gravitational plus `` active '' components ) , one gets    & & 1 - _ c0 _ c = + & & a_c _ r^2 _ 0r _ 0^_r ^2 _
0^n  d  .    clearly , taking @xmath175 yields @xmath176 .
also , notice that @xmath177 is not perturbed : this would reflect a change in @xmath161 ( since @xmath178 ) , which has already been ruled out .
also , notice that the denominator affecting the integral is equal to @xmath179 : thus , one can interpret this integral as the volume averaging of the expression in brackets within the integrand , with the density as weighting function , that is , @xmath180 .    from the polytropic equation of state ,
one gets @xmath181 , which enables writing    & & 1 - ( t_c0 t_c)^n+1 = + & & a_c _ r^2 _ 0r _ 0^_r ^2 _ 0^n  d  .
the above results allow us to construct an exclusion plot for the central temperature on the @xmath182 parameter space , by imposing the constraint @xmath183 , the uncertainty of the central temperature of the sun .
solving numerically the differential equation for @xmath138 for a varying @xmath100 , and computing all relevant quantities for varying @xmath135 , and defining the absolute perturbation @xmath184 , so that @xmath185 , one obtains the results shown in the figures [ absprofiles ] to [ mass ] .
= 5.7 cm    = 5.7 cm    = 5.7 cm    = 5.7 cm    = 5.7 cm    = 5.7 cm    = 5.7 cm    first , notice that , from fig .
[ absprofiles ] , the absolute perturbation @xmath186 is fairly insensitive to the value of @xmath100 ; clearly , given that the profile of @xmath137 for varying @xmath100 differs more sharply as @xmath187 ( as seen in fig .
[ absprofiles ] ) , a approximately constant @xmath186 translates into a relative perturbation @xmath138 that also exhibits this behaviour , as can be seen from fig .
[ relprofiles ] .
furthermore , notice the indication of the divergence of the relative perturbation @xmath138 in fig .
[ relprofiles ] ( which does not depend on @xmath100 , since @xmath188 and @xmath189 yields @xmath190 , as previously discussed ) ; as stated before , this divergence is avoided by dealing only with the radiative region @xmath191 .
also , as can be seen from fig .
[ relprofiles ] , the relative perturbation @xmath138 peaks at @xmath192 ( for @xmath193 , and a smaller value of @xmath194 for @xmath128 ) ; hence , the perturbative condition @xmath195 translates to @xmath124 , leading to the chosen ( extreme ) interval for simulation @xmath196 .
hence , according to figs .
[ reltc ] , [ reltc3d ] and [ reltc2d ] , no relative deviation of the central temperature above the experimentally determined level of @xmath197 is attained . however , the values found , of the order of @xmath198 , indicate that any future refinement of the experimental error of @xmath199 could yield a direct bound on the parameter @xmath135 .
the above results show that @xmath135 is presently unconstrained by solar observables , aside from the perturbative condition @xmath200 ; taking @xmath201 and @xmath202 , this yields @xmath203 or , in natural units , @xmath204 .
however , this study assumes that the effects arising from the non - trivial matter - geometry coupling supersede those from the non - linear curvature term in the modified hilbert - einstein action ( [ action ] ) , as expressed by the inequalities ( [ ineqa ] ) , ( [ ineqf1a ] ) and ( [ ineqf1b ] ) .
having numerically determined the perturbative solution @xmath138 , one can now re - evaluate these conditions . for this , one first replaces @xmath10 by its unperturbed value ( with the newtonian approach @xmath205 ) , @xmath206 , since corrections would amount to second - order terms in @xmath13 .    considering the scaling laws , @xmath112 and @xmath113 , and the inequalities ( [ ineqf1a ] ) and ( [ ineqf1b ] ) , some algebra ( see appendix d ) yields , for the @xmath207 case ,    [ ineqla1 ] || & & ( 2 _ c ^n ) ^m+1 2 n_c  , +        while , for the @xmath208 case ,    [ ineqla2 ] & & 2 n_c  , + & & n n+1 n_c  , + & & n n+1 n_c c_n,-1 ( )  ,    where one defines the quantity    c_n , m ( ) | ( 1 + nm ) ^2 ( ) -  ( ) n ^2()+  ( ) |  .
one can plot the function @xmath209 by evaluating @xmath132 by the unperturbed solution @xmath137 , using eq .
( [ le0 ] ) and taking @xmath210 .
one obtains    c_n , m ( ) | _ 0 + _ 0 ^ 4 _ 0 ( 3_0 + 2 _ 0 ) + _ 0 ^ 4|  .
varying @xmath211 in its domain @xmath212 ( for the @xmath213 case ) and @xmath49 ( for the @xmath214 case ) ; its profile is depicted in figure [ cnmgraph ] .
one concludes that @xmath215 , for both the @xmath216 case ( @xmath212 ) as well as the @xmath217 case ( @xmath49 ) .
this , together with the constraint @xmath218 ( as can be seen from fig .
[ absprofiles ] ) and the quantity @xmath219 yields , for the @xmath47 case ,    [ ineqla3 ] || 2.8 10 ^ 43 ( 8.48 10 ^ 36)^m ( 1-m)m  ev^-2(1+m )  ,    while , for the @xmath48 case ,    [ ineqla4 ] || 6.6 10 ^ 6  .
using the previously discussed value , @xmath220 , which arises from the planck - scale renormalization of the theory , one gets    [ ineqla5 ] || 6.6 10 ^ -36  ev^-2 = ( 3.9 10 ^ 8  gev)^-2  .
also , a recent paper has reported a relation between the coefficient of the quadratic term in @xmath51 and the ppn parameter @xmath221 @xcite ,    = 1 2  .    hence , the current experimental constraint @xmath222 @xcite yields    0.17 m_pl^-2 = ( 3.0 10 ^ 19  gev ) ^-2  ,    where @xmath223 is the planck mass . this yields    [ ineqla6 ] || 7.4 10 ^ -51  ev^-2 = ( 1.2 10 ^ 16  gev)^-2  .
clearly , the @xmath48 case does not impose any significant bound on @xmath13 , with the above value laying many orders of magnitude below the upper bound , @xmath204 , arising from the perturbative treatment .
notice that the above inequalities ( [ ineqla1 ] ) , ( [ ineqla2 ] ) , ( [ ineqla3 ] ) , ( [ ineqla4 ] ) , ( [ ineqla5 ] ) and ( [ ineqla6 ] ) are not to be considered as restrictions on the parameters , but as conditions for the validity of the regime where the effects of the matter - geometry coupling supersede those of the the non - linear curvature term in eq .
( [ action ] ) .
in this work we have examined a model where the usual einstein - hilbert action is modified , not only by allowing non - linear curvature terms to appear , but by also enabling an explicit , non - minimum coupling between the curvature and the lagrangian of matter fields .
we first assume that this matter - geometry coupling is linear in curvature and its effect allows the non - linear curvature terms to be neglected  an assertion that is qualified in the end of the study , for two relevant non - linear @xmath1 models .
we then proceed and perform the necessary calculations to ascertain the effect of the latter in the hydrostatic equilibrium of an @xmath210 polytrope such as the sun .
we assume a perturbative regime to the usual tolman - oppenheimmer - volkoff equation , and take into consideration the validity of the newtonian regime in this modified theory , as well as the redefinition of relevant quantities , which are computed numerically and compared with solar observables .
this goal is achieved through the use of the ( non - relativistic ) polytropic equation of state @xmath97 : as stated before , this is a very simplistic description of the complex behaviour of solar matter , and has been superseded by much more elaborate equations of state ; we use it in our study so to better illustrate the effects of the matter - curvature coupling on an easily understandable , closed model .
however , we should remark that this restrictive treatment , although advantageous from the theoretical point of view , might conceal some of the more intricate phenomenology found in stellar systems , which could affect our results .
clearly , a subsequent study should consider a more realistic solar model .
the results allow us to conclude that no strong constraints on the matter - geometry coupling from the comparison between the model s predictions and current experimental sensitivity , aside from the perturbative approach @xmath224 , which yields @xmath225 .
however , the numerically obtained results show that a slight increase in accuracy would allow an upper bound to be placed on @xmath13 ; this closeness between the prediction of the perturbative model and experiment also seems to validate the latter .                                in order to obtain an approximation to the solution of the lane - emden equation ( [ le ] ) in the vicinity of @xmath228
, one writes @xmath229 . inserting this in eq .
( [ le ] ) yields , after some algebra and expanding up to fourth order in @xmath118 ,                            in the above , @xmath137 is used instead of @xmath132 , since @xmath230 is coupled to @xmath135 , so that including the perturbation @xmath138 would only amount to second - order terms @xmath174 . using the definitions @xmath146 and @xmath231 ,
one gets    & & m = _
0^r_r m_e  dr = + & & 4 _ 0^r_r r^2  dr + n+1 g a_c n_c _ 0^r_r ^2 _ 0^n  dr = + & & 4 r_0 ^ 3 _ c _ 0^_r ^2 ^n  d+n+1 g a_c n_c r_0 _ 0^_r ^2 _ 0^n  d= + & & 4 _ c _ 0^_r ^2 ( _ 0 ( 1 + a_c ) ) ^n  d+ + & & ( 4 r_0 ^ 2 _ c^2 p_c ) a_c n_c r_0 _ 0^_r ^2 _ 0^n  d + & & 4 _ c _ 0^_r ^2 _ 0^n ( 1 + n a_c )  d+ + & & 4 r_0 ^ 3 _ c a_c _ 0^_r ^2 _ 0^n  d= + & & 4 _ c ( _ 0^_r ^2 _ 0^n  d+ a_c _ 0^_r ^2 _ 0^n ( n + )  d ) = + & & 4 _ c ( - _ r^2 _ 0r + a_c _ 0^_r ^2 _ 0^n ( n + )  d ) = + & & -4 _ c _ 0r ( 1 - a_c _ r^2 _ 0r _ 0^_r ^2 _ 0^n ( n + )  d )  .          & & || = ( 2 _ c ) ^m+1 2 n_c ^n(1+m)+1  , + & & || 2^m ( 1-m ) m ^1+m
| (r ) ^1+mp |=
+ & & n n+1 ( 2 _ c ) ^m+1 ( 1-m ) m 2 n_c ^n(1+m ) + 1  , + & & || |  ( r ) -(1+m)^2(r ) | = + & & ( 2 _ c ) ^m+1 ( 1-m)m 2 n_c ^n(1+m)+1 | ( 1 + nm ) ^2 ( ) -  ( ) n ^2()+  ( )        where the prime denotes denotes differentiation with respect to @xmath118 ( the factors @xmath171 arising from changing derivation with respect to @xmath90 to derivation with respect to @xmath114 cancel out ) .
the authors would like to thank f. lobo for fruitful discussions .
the work of j.p . is sponsored by the fct under the grant @xmath232 .
o.b . acknowledges the partial support of the fct project @xmath233 .
s. carroll , v. duvvuri , m. trodden and m. turner , , 043528 ( 2004 ) ; s. capozziello , s. nojiri , s.d .
odintsov and a. troisi , , 135 ( 2006 ) ; s. nojiri and s.d .
odintsov , , 086005 ( 2006 ) , int .
phys . * 4 * , 115 ( 2007 ) .
l. amendola , , 043511 ( 2000 ) ; , 196 ( 2001 ) ; l. amendola and d. tocchini - valentini , , 043509 ( 2001 ) , phys .
d * 66 * , 043528 ( 2002 ) ; a billyard and a. coley , , 083503 ( 2000 ) ; d. holden and d. wands , , 043506 ( 2000 ) ; o. bertolami and p. martins , , 064007 ( 2000 ) ; l. chimento , a. jakubi , d. pavon and w. zimdahl , , 083513 ( 2003 ) ; l. chimento and d. pavon , , 063511 ( 2006 ) ; c. boehmer , t. harko and f. lobo , gr - qc/0709.0046 .
v. faraoni , _ cosmology in scalar - tensor gravity _
( kluwer academis publishers , 2004 ) ; o. bertolami . j. pramos and s. turyshev , gr - qc/0602016 ; o. bertolami , f. pedro and m. le delliou , astro - ph/0703462 ; astro - ph/0705.3118 .
t. multamaki and i. vilja , astro - ph/0612775 ; k. kainulainen , j. piilonen , v. reijonen and d. sunhede , , 024020 ( 2007 ) ; a. bustelo and d. barraco , , 2333 ( 2007 ) ; e. barausse , t. sotiriou and j. miller , gr - qc/0703132 ; s. capozziello , a. stabile and a. troisi , gr - qc/0709.0891 .
`` the fundamentals of stellar astrophysics '' , g. collins ( w.h .
freeman and co. , 1989 ) ; j.n .
bahcall , 47 ( 2000 ) ; `` textbook of astronomy and astrophysics with elements of cosmology '' , v. bhatia ( narosa publishing house , 2001 ) ; `` theoretical astrophysics : stars and stellar systems '' , t. padmanabhan ( cambridge university press , 2001 ) .


bl lac objects are one of the several flavors of radio loud active galactic nuclei . along with peculiar properties at other wavelengths ( lack of strong emission lines , high levels of variability , optical polarization up to 3% , detection at x - ray and @xmath0-ray energies ) , they also show quite extreme behavior in the radio : the parsec - scale structure is dominated by a compact , flat - spectrum core , from which a one - sided jet emerges , showing knots of enhanced brightness . when monitored over long intervals ,
these components are often found to be in motion , sometimes with an apparently superluminal velocity ( see , e.g. , * ? ? ?
the currently accepted explanation for these properties is found in the frame of unified models ( see e.g. * ? ? ?
these models are based on a dust enshrouded super - massive black hole accreting matter in a disk , from which two symmetric collimated jets of relativistic plasma are ejected in opposite directions .
when the angle between the jet axis and the line of sight is small , the resulting doppler boosting accounts for many of the above - mentioned properties .
a source such as markarian 501 ( b1652 + 39 ) is an ideal target to test the assumptions involved in unified models . at its low redshift ( @xmath1)kms@xmath2mpc@xmath2 ] , 1mas  =  0.72pc ; therefore , the milliarcsecond resolution achievable with vlbi techniques makes it possible to investigate in great detail the regions near the central core .
furthermore , mkn  501 is one of the few sources with a clear limb - brightened structure @xcite .
finally , this source is very well studied at other frequencies and it is known for its x@xmath4 and @xmath5ray @xcite activity , which also sets constraints on the parameters describing the physics of the inner jet ( see * ? ? ? * ; * ? ? ?
* and references therein ) .    in this paper we present a large multi - epoch , multi - frequency data - set of vlbi images ; new observations as well as previously published data are used to study the physical properties of this source .
in particular , we will discuss the pc - scale jet morphology , velocity and orientation .
additionally , our data - set allows us to compare images obtained at different frequencies with similar resolutions and thus to perform spectral index such that @xmath6 studies for comparison with the total intensity images .
in  [ sec : observations ] we describe the previous and new observations used in this paper , along with the data reduction methods and procedures for component fitting and spectral index mapping .
results are presented in  [ sec : results ] and discussed in  [ sec : discussion ] .
we present our conclusions in  [ sec : conclusions ] .
table 1 summarizes all observations considered in this paper , listing epoch , frequency , total observing time and array used .
the notes provide references for details .
the present data - set is very similar to the one used by edwards & piner ( 2002 , hereafter @xcite ) to discuss the proper motion in mkn 501 : we did not use the 1995 oct 17 epoch ( due to its low sensitivity ) , or the 1996 jul 10 and 1998 oct 30 epochs ( as better data from other nearby epochs was available ) . on the other hand
, we add 3 more epochs ( see note 4 to table 1 ) at 15 and 22ghz and we consider two new space vlbi ( vsop ) observations obtained on 2001 march 5 and 6 ( see table  2 ) .
images from some of the previous observations have been published or made available to the scientific community on open web sites . however , for the sake of completeness and homogeneity , we requested and obtained the calibrated ( @xmath7)-data for all observations presented here .
after the initial reduction ( see notes to table  1 for references ) , we imported the calibrated data into the nrao astronomical image processing system ( aips ) and inspected them carefully .
then , we produced new images ; one or more self - calibration cycles proved to be useful in some cases , providing a better signal to noise ratio .
we adopted the appropriate ( @xmath7)-range , weights , cell - size and angular resolution to map the spectral index and to look for proper motion of the jet substructures ; we present parameters of final images in table  3 .
we obtained space vlbi images at three different epochs : 1997.8 , 1998.4 , and 2001.3 .
vsop ( vlbi space observatory programme ) observations combine data from ground arrays , such as the vlba , and the 8 m orbiting antenna on board the satellite halca . for the latter two epochs observations at 1.6 and 4.8ghz were performed on successive days to minimize flux density variability problems in spectral index mapping .
the first space vlbi observation was conducted at 1.6ghz only .
further details on these observations are given in table  2 .
the standard vsop observing mode has one polarization ( lcp ) and two ifs of 16 mhz each with two - bit sampled data . the lowest if failed in the 1.6ghz observation in 1997 . in 2001 ,
the `` a priori '' calibration was affected by lack of system temperature measurement for robledo and goldstone .
a mini - cal procedure was done before starting the experiment in robledo , providing a value of 20k that was assumed as constant . for halca
, we used the nominal system temperatures of 75k and 77k for the two ifs at 1.6ghz , and 88k and 92k for those at 4.8ghz @xcite .
all the data were correlated in socorro and then imported and reduced in aips .
global fringe - fitting returned good solutions for all the observations in which large ground telescopes were used , although no solutions were found for the last observation ( 2001 mar 6 ) , unless a value of snr as low as 4.5 was set in the fit . in spite of such a low signal to noise ratio
, we found these solutions reasonable : delays and rates found in this way do not present any discontinuity or suspicious behavior .
thus , we applied these solutions and proceeded with the whole array .
we then performed self calibration using clean - component models , first in phase and , afterward , in both phase and amplitude . in the latter
, we repeated the task several times , increasing the baseline range .
this revealed to be of great importance for those observations where gain information for some antennas was missing .
for example , we could self - calibrate the shortest space - ground baselines using ` uniform weight ' ground - only images ; in fact , we have similar baselines both in the ground array and to halca , as provided by its highly elliptical orbit ( see fig .
[ uv1 ] , [ uv2 ] ) .
figures [ fig3 ]  [ fig9 ] present all the images , in order of decreasing angular resolution .
identification of components is usually the most difficult problem affecting proper motion studies .
long time intervals between epochs can lead to misidentification of components ; moreover , the use of different resolutions can cause confusion , for example by blending components with different properties and velocities ; finally , registering images at different observing frequencies is also necessary : the apparent core shift due to the different spectral index of the core and jet can be significant .    in the present work ,
the time coverage is very dense ( 14 epochs over 4 years , plus the final vsop observation ) and there are no large time gaps . to prevent angular resolution problems , we produced maps at the same resolution at all epochs . even if some observations allowed a better resolution , we convolved all the images obtained at a frequency 5ghz@xmath815ghz with a circular restoring beam of 1.2mas fwhm
however , higher resolution maps have also been made , in order to provide more information when needed .
in particular we obtained images at a resolution of 0.6 @xmath9 0.9 mas ( ra @xmath9 dec ) in all observations with good @xmath10-coverage on long baselines ( see table  3 ) .
a comparison of fit components in convolved and full resolution images revealed good agreement .
the fit to the components of the jet was done within aips and the difmap package .
in aips we used the task ` jmfit ` , which allows one to fit gaussian components to part of an image . up to
four components can be used at a time and the task , starting from given initial conditions , solves in an iterative way for position , flux density and extension of each component . since more than four components are needed to properly fit the images , we had to run the task several times on different regions .
the regions were made to overlap with each other in order to obtain results compatible with previously published models .
we also repeated the fit looking for `` stable '' solutions , i.e. , changing our initial guesses or even leaving them blank .
eventually , very good models were obtained ; the difference in flux density between our fits and the images is always smaller than 5% .    using difmap
, we also carried out model - fitting to the visibilities rather than the image . similarly to the fits made in aips , four to six jet components were required in addition to the core to provide a good representation of the data .
when producing spectral index images , there are three main effects that have to be carefully taken into account .
first , because of flux density variability and possible proper motion of components , the observations must be taken at very close epochs .
next , the angular resolution must be the same at both frequencies and the @xmath10-coverage as similar as possible ; this guarantees the same sensitivity to extended structures in the high and low frequency images . finally , as absolute position
is not available , we need to properly align the images before combining them .    to meet the first two requirements , we produced spectral index images using data separated by less than 24hrs .
we considered 15 and 22ghz data obtained on the same day ( 1997 aug 15 ) and space vlbi observations at 1.6 and 5ghz made on successive days ( 1998 apr 7 and 8  see * ? ? ? * for a preliminary analysis ) .
we note that 1998 space vlbi observations have better @xmath10-coverage than those in 2001 ( see fig .
[ uv1 ] , [ uv2 ] ) . to obtain matching @xmath10-coverages
, we cut the shortest baselines from the low frequency data and the longest baselines in the high frequency observations .
images were then obtained using uniform weights , and the same cell - size and restoring beam .
the problem in registering the images is that the peak of the image , coincident with the core , can correspond to different positions at different frequencies if the core is self - absorbed .
so , we aligned the images using the positions of components along the jet in a region with steep spectral index and no self - absorption effects .
we found that at 15 and 22ghz no shift was necessary , while a shift of @xmath110.2mas was necessary in both ra and dec between the 5 and 1.6ghz images .
see figures [ spix1 ]  [ spix3 ] for the final spectral index maps .
the large scale radio structure has been imaged with the vla @xcite , showing a core dominated source with two - sided diffuse emission oriented at @xmath12 ; this indicates that the jets are non relativistic at a ( projected ) distance of @xmath13kpc from the core .
the pc - scale structure has also been investigated , using vlbi techniques ( see e.g. * ? ? ? * ) , with resolution of @xmath13 mas .
however , the collection of multi - epoch , multi - frequency , high resolution and high sensitivity images will improve our knowledge of the source morphology and evolution .
the present vlbi observations show a strong core and a one - sided jet .
the jet exhibits multiple sharp bends before undergoing a last turn , followed by rapid expansion .
we can distinguish three different regions :    * a first region , extending @xmath14mas from the core , where a high brightness jet structure is present ( fig .
[ fig3 ] , [ fig4 ] , and [ fig5 ] ) .
the jet pa is not constant in this region , being 150@xmath15160@xmath15 near the core ( @xmath161mas ) and moving to @xmath17 ( from 1 to @xmath114mas from the core ) to become @xmath18 from 4 to 10mas .
+ the 22ghz high resolution image ( fig .
[ fig3 ] ) shows that the jet is resolved and limb - brightened starting at @xmath111mas from the core .
this is confirmed by the high sensitivity 5ghz space vlbi image ( fig .
[ fig5 ] ) ; although this image does not have enough angular resolution to resolve the inner jet near to the core , it shows an extended low brightness emission , confirming the transverse extension of the jet in the inner 2mas . in the region from 2 to 10mas from the core ,
the limb brightened jet structure is clearly visible .
+ the brightness profile along the main axis of the jet is not uniform ; high brightness regions are followed by dimmer regions . in particular at 56mas from the core a deep minimum is present just before an extended bright spot with a peak 89mas from the core .
however , we do not refer to these regions as `` knots '' , since they are resolved at higher frequencies . * between @xmath1110mas and @xmath1130mas from the core
, the jet loses its collimation and becomes visible only in lower resolution images ( fig .
[ fig6 ] , [ fig7 ] ) . in this region ,
the jet is oriented at a pa of @xmath11 110@xmath15 and its opening angle is increasing .
a well defined limb brightened structure is visible ( see , e.g. , fig .
[ fig7 ] ) , in agreement with the `` layered '' structure visible in the polarized emission @xcite . as in the first region ,
lower frequency images reveal a wider jet and a more evident external shear layer ; this implies that the spectra in the inner regions are flatter than in the external regions .
* at @xmath1130mas from the core the jet shows strong bending and a large opening angle ( fig .
[ fig8 ] , [ fig9 ] ) . in this third region
the orientation of the diffuse jet is in agreement with the orientation of the kpc - scale emission .
the jet opening angle is @xmath1143@xmath15 .
the radio brightness is uniform and there is no indication of the helical structure suggested by @xcite , which possibly resulted from the poor @xmath19-coverage ( five telescopes only ) of their observations .
a limb - brightened structure is visible in our low resolution images up to @xmath11100mas from the core , as can be seen in fig .
[ fig9 ] .
no radio emission is detected either on the opposite side of the radio core or in the position of the putative counterjet feature f reported by @xcite . however , if we compare our measured total flux densities to single dish data ( e.g. from the umrao database ) , the correlated vlbi flux density accounts only for a part of the total flux density ( @xmath20% at 5ghz and @xmath21% at 15ghz ) ; this is in agreement with the presence of a larger scale structure resolved out by the vlbi data , such as the two sided emission shown by @xcite and @xcite .
single dish monitoring of markarian 501 shows that this source is remarkably stable in the radio band .
@xcite monitored a sample of 23 @xmath5 and x@xmath4 ray loud blazars over a 4 year period ( 19962000 ) with the 32 m medicina radio telescope . from their data ,
markarian 501 is the most quiescent source , with an almost flat light curve at both 5 and 8.4ghz .
moreover , markarian 501 is part of the umrao database , which also shows a fairly steady flux density . during the period of our observation , the average flux densities are @xmath22jy at 4.8ghz , @xmath23jy at 8ghz and @xmath24jy at 14.5ghz .
we consider the vlbi core flux density using our data at 15ghz ( 0.6@xmath90.9 mas fwhm , ra@xmath9dec ) , where we have 8 different epochs ( see table  3 ) . assuming the map peak as the core flux density , we see some dispersion in the values , which range between 360mjy ( 1997 apr 25 ) and 508mjy ( 1995 dec 15 ) .
this suggests possible core variability , which is diluted in single dish observations , since the vlbi core makes up only @xmath25% of the total flux density .
we note , however , that our observations are not homogeneous and were not aimed to monitor flux density variability .
the overall morphology of the source does not present dramatic changes between epochs .
no compact jet component ( knot ) is present in the jet at all frequencies : if we consider the highest resolution data , the local peaks present in the low resolution images become clearly resolved . therefore , we limit our study to the innermost region , where high - resolution images are available ( the jet brightness is too low to be detected at high frequency beyond 10mas from the core ) . in this region ,
the jet is always transversally resolved : it appears centrally peaked at low resolution , but resolved with a limb brightened structure at high resolution .
this implies that we could be mixing together different components with different properties and moving at different velocities .
for all these reasons , we analysed images at the best angular resolution ( typically 0.6@xmath90.9mas , ra@xmath9dec ) and , in model - fitting , we tried to use a number of components that did not produce an oversimplified representation of the data . in particular , to avoid resolution or frequency effects ,
we used only 15ghz data , plus the 8.4ghz observation at epoch 1998.48 .
the results obtained from model - fitting of visibilities in difmap are presented in table  4 and fig .
[ positions ] .
these results are in agreement with those obtained by fitting gaussian components to the images ( using ` jmfit ` in aips ) .
the data are in general well fitted by a core and four jet components , in agreement with @xcite ( not surprisingly , as most data are the same ) .
we labelled components according to the same notation for ease of comparison : c1 to c4 moving from the outermost component inward .
however , for the reasons stated above , we adopted two gaussians for each of the components c2 and c1 at many epochs .
this produced a better fit to the data , in agreement with the resolved limb - brightened jet structure present in these regions .
nominal position errors are very small .
the major uncertainty is related to the complex source structure .
we estimate a position error of @xmath260.1mas for c4 , c3 , c2 and @xmath110.2mas for c1 .    from our analysis we conclude that no proper motion is visible in any of the four components .
consider the case of component c2 : @xcite measured a proper motion of @xmath27 for this component .
however , we note that at the position of c2 two faint external peaks are visible in the 22ghz image ( fig .
[ fig3 ] ) on either side of the lower frequency component .
they could even have a different spectrum ( fig .
[ spix1 ] ) .
the peak position visible at lower resolution in the 15ghz image is therefore the blend of the limb - brightened structure visible at 22ghz . in the space vlbi image at 5ghz , the component c2 is at a different position with respect to the 15ghz images : this could be due to different spectra for the different components .
we conclude that it is not possible to measure a proper motion for such a complex structure and that any pattern velocity measurement from low resolution data is not related to the jet bulk speed in this source .      under the assumption of intrinsically symmetric two - sided jets affected by doppler boosting effects ,
the observational data allow one to estimate the jet bulk velocity . in the following
, we present the results derived from a study of the jet / counterjet ratio , core dominance and jet expansion .      from the lack of detection of a counter - jet in all the images (  [ sec : morphology ] ) , we set a lower limit on the jet / counterjet ratio .
we use the highest ratio between the jet brightness and the counterjet upper limit .
we estimate this ratio excluding the innermost region of the jet ; this avoids contamination from the unresolved central emission . the best image for this analysis is the space vlbi observation obtained in apr 1998 at 1.6ghz . from this image
we can derive a jet / counterjet ratio @xmath28 ( 1  rms ) at @xmath114mas from the core . given the well known relation @xmath29 , and assuming
@xmath30 ( see  [ sec : index ] ) , this value yields @xmath31 . in turn , this implies @xmath320.89 and @xmath3327@xmath15 .
our deep , low resolution images allow us to also put a constraint at large distance from the core . in particular , in the 1.6ghz image from the 2001.3 data , at @xmath34 mas from the core
we still have r @xmath35 70 ; although less extreme , this value is also very interesting , as it implies that the jet remains at least mildly relativistic ( @xmath36 ) even at large distances from the core .
it is not possible to extend this analysis further with our data .
however , we remind the reader of the symmetric morphology of the source on the kpc scale (  [ sec : morphology ] ) .
this implies that the bulk of the jet slows to subrelativistic speed going from the parsec to the kpc scales .
given the general relation between total and core radio power @xcite , we can use low frequency unboosted data to compute an `` expected '' value for the core power . a comparison between the expected and measured 5ghz core radio power allows to constrain the jet orientation and velocity @xcite .
we note that at 5ghz self - absorption effects are small ( see  [ sec : parameters ] ) ; however , to take also into account possible core variability , we allow the observed core flux density to vary by a factor of 2 .
this is a conservative assumption , based on the total flux density monitoring and the range of derived vlbi core flux densities ( see  [ sec : variability ] ) , and , more generally , on statistical flux density variability of agn in the radio .
to estimate the orientation and jet velocity for mkn 501 , we have compared its total unboosted flux density ( 1.81jy ) at 408mhz @xcite , with the nuclear flux density at 5ghz . since @xcite used the arcsecond core radio power to derive their correlation , we adopted the total vlbi correlated flux density at 5ghz ( @xmath11800mjy ) as the mkn 501 core flux density .
any possible underestimate in this assumption is compensated for by the allowed range ( a factor of 2 ) in flux density variability .
we found that mkn 501 has to be oriented at @xmath37 27@xmath15 with @xmath380.88 .
we note that a high jet velocity ( @xmath390.95 ) is allowed only in the range 10@xmath15 @xmath16 @xmath40 @xmath16 27@xmath15 .
high velocities are forbidden at smaller angles ; the resulting doppler factor would require a core much brighter than observed .
the functional dependence of the jet intensity on the jet velocity and radius for an adiabatically expanding jet has been discussed in the case of relativistic motion by @xcite .
since the model represents a simplification of the real situation , the results should be considered with caution and compared with other observational results .
this model has enabled useful constraints to be placed , in agreement with other data , for 3c264 @xcite and other sources , such as ngc315 @xcite and 3c449 @xcite .    to apply this model we used our deep vsop observations at 1.6ghz obtained in march 2001 and in april 1998 . in these observations , the jet emission is visible on scales up to @xmath41mas ( see , e.g. , fig .
[ fig9 ] ) .
first , we derived brightness profiles across the jet ; we used the aips task ` slice ` on an image reconstructed with a circular beam of 7mas diameter , where the longest baselines have been omitted to increase the signal to noise ratio . in the innermost region ( @xmath42mas from the core )
, we have 8 slices , taken every half beam - width at pa 27.2@xmath43 .
after the main bend , we considered 12 more slices , each one beam apart , at pa @xmath44 .
we stopped our analysis @xmath45mas from the core ; further out , the low jet brightness and the presence of extended sub - structures do not allow the data to be fitted well .    using the aips task slfit , we fitted single gaussians to each profile .
this could be done unambiguously in most cases , although in a few slices some deviation from a pure single gaussian profile was present .
however , the difference between the area subtended by the profile and the fit is always smaller than 5% and does not affect the analysis .
we plot the resulting fwhms and peak brightnesses in fig .
[ trends1 ] and [ trends2 ] ; the quantities have been deconvolved from the clean beam , according to the formula given by @xcite .
there are clearly two regimes ; least - squares fits yield a power - law of index @xmath46 in the inner jet ( @xmath47mas from the core ) and @xmath48 in the outer regions . in the following
, we use these best - fit solutions rather then the actual data ; this prevents us from spurious results brought about by random fluctuations .
we recall the equations given in @xcite for the jet surface brightness @xmath49 :    @xmath50 ( predominantly parallel magnetic field )    @xmath51 ( predominantly transverse magnetic field ) .
we assume an injection spectral index @xmath52 , as suggested from the average @xmath53 in the image , excluding the self - absorbed core (  [ sec : index ] ) .
it is not easy to decide which magnetic field regime is preferable : according to the 8.4ghz images of @xcite , the mkn 501 jet is highly polarized at the edges with the magnetic field parallel to the jet flow , while in the jet spine the magnetic field is orthogonal to the jet flow with a lower degree of order .
@xcite report a 5ghz electric vector pa distribution consistent with this result but the lower resolution renders the results less conclusive .
therefore , since the percentage of ordered field and the intrinsic orientation are not well known , we consider the two cases of pure parallel and perpendicular magnetic field orientation ; in a disordered magnetic field , the real result will be in between the two extreme situations .
figure  [ adiabatic ] shows the derived trends of @xmath54 for an initial lorentz factor @xmath55 and @xmath56 ( @xmath57 ) with parallel and perpendicular magnetic fields . in each plot
, we draw five lines , corresponding to angles to the line of sight of 5@xmath43 , 10@xmath43 , ... , 25@xmath43 ( i.e. , in the range of values allowed by the jet sidedness and core dominance ) .
we note that the jet velocity decreases with the core distance , more slowly in the perpendicular case , and on a shorter scale in the parallel case .
the lack of detection of the cj in all our images strongly constrains any possible model since it implies a relativistic jet also at large distance from the core ( @xmath58pc ) .
models with low initial lorentz factor ( @xmath59 ) are definitely ruled out , regardless of the assumed magnetic field orientation , since they disagree with both the observed limb - brightened structure (  [ sec : jvs ] ) and the jet / counterjet ratio ( at least in the case of the parallel magnetic field ) ; moreover , they require a jet deceleration (  [ sec : jvs ] ) between the @xmath5ray region and radio jet region that is too strong . among models with a higher initial speed ( @xmath60 ) ,
the case of the parallel magnetic field requires an orientation angle of 25@xmath15 , at least at @xmath3570mas from the core , to be in agreement with the observed jet / counterjet brightness ratio .
we can not exclude models with a narrow angle to the line of sight at small distances from the core , where solutions are similar .
in the case of the perpendicular magnetic field , the constraints are less severe since we have a fast jet even at large distances from the core .      in fig .
[ spix1 ] , [ spix2 ] and [ spix3 ] , we display the spectral index distribution at various resolutions between 1.6 and 4.8ghz and between 15 and 22ghz .
observational results are summarized in table  5 for high resolution images and described in the following subsections .      at low resolution ( @xmath62 mas , pa=@xmath63 )
, we can study the spectral index distribution at distances from the core larger than 810 mas ( fig .
[ spix1 ] ) .
the limb - brightened jet structure between 8 and 20mas from the core is clearly visible .
the spectrum is always flat in the inner spine ( medium value @xmath53=0.08 ) and steep in the bright shear regions , with values in the range @xmath64 . in the inner region
( @xmath168mas from the core ) , the jet is well resolved with a steep spectrum in the more external regions ( @xmath65 ) but with a puzzling region south of the high brightness region where the spectral index is @xmath66 .
this region is clearly visible also in the lower resolution image in @xcite where it seems connected to the inner spine flat spectrum region , while no such connection is present in the total intensity images .
using the space vlbi baselines , we can study the source details at higher angular resolution ( @xmath67mas , pa=@xmath68 ) even at this frequency ( fig .
[ spix2 ] ) .
the nuclear spectrum is inverted ( @xmath69 ) . except for c4 ( confused with the core and subject to a large uncertainty ) , the jet sub - components show flat to moderately steep spectra : @xmath70
. regions in between the components have steeper spectra .      at these high frequencies ,
we are beyond the self - absorption turnover and the nuclear region shows a flat , non - inverted spectrum ( @xmath53=0.45 ) , generally steepening along the jet .
the resolution of 0.6@xmath90.9mas ( ra@xmath9dec , fig .
[ spix3 ] ) shows that the inner regions c4 , c3 , c2 have @xmath71 in the range @xmath72 , while regions in between are much steeper ( @xmath73 and higher ) .
the large component c1 is too faint and resolved to be properly imaged at these high frequencies .
a marginal suggestion is present of a spectral steepening moving from the inner ( spine ) to the external jet region but at this high angular resolution the jet width is not well defined .
it is interesting to consider the nuclear flux density over the whole range of frequency .
we are aware that the compact core of mkn 501 presents some variability ( see  [ sec : variability ] ) ; however , the large number of observations and frequencies available allow us to carry out such a study based on average values .
we plot in fig .
[ spectrum ] the average core peak flux density , from the images at resolution of 0.6@xmath90.9mas ( see table  3 ) ; at 1.6ghz , we had to consider super - resolved images , in order to avoid a large contamination from the jet .
the spectrum is convex , peaking at @xmath74ghz and steepening toward the high frequency end .
below the turnover frequency , the spectrum is flatter than @xmath75 , as would be expected for a single self - absorbed component .
this suggests that different components with different turnover frequencies may be blended together .
limb brightened jet morphology on parsec scales is present in some fr  i sources but also in a few high power fr  ii sources @xcite . a possible explanation is the presence of velocity structure in the jet .
when viewed at a relatively large angle to the line of sight , this structure yields different doppler factors for the internal spine and the external shear layer .
in particular , an inner high velocity spine could be de - boosted while a slower external layer could be less strongly de - boosted , or even boosted .
transverse jet velocity structure has been found by @xcite in 3c31 , where the jet velocity in outer regions is @xmath110.7 the inner spine velocity .
@xcite invoke such a structure to account for observational discrepancies between fr  i and bl - lac objects .
a possible origin of the velocity structure is the interaction of the jet with the surrounding medium and its consequent entrainment ; the mass loading slows down the jet in the external regions , while the inner spine travels unaffected .
however , the presence of limb brightened jets in both high and low power sources @xcite and the evidence in m87 that the jet structure is present very near to the core emission ( see the image presented by * ? ? ?
* ) , suggest the possibility that the velocity structure may originate at the base of the jet .
@xcite presents a model where the inner high velocity spine is produced in the black hole region , while the external shear moving at lower velocity is produced in a more extended region downstream .    in mkn 501
, it is interesting to note that the jet structure is visible at a distance @xmath161mas from the core , suggesting that it originates very close to the base of the jet .
however , we recall that , at @xmath1 , 1mas corresponds to @xmath76pc , i.e. , @xmath77 schwarzschild radii ( see * ? ? ? * and references therein for a discussion on the central black hole mass in mkn 501 ) and , therefore , we do not have the angular resolution to see this jet inner region .
nevertheless , we note that the presence of velocity structure so near to the nuclear source implies too strong a deceleration if it is completely due to the jet interaction with the ism .
we suggest ( in agreement with m87 data of * ? ? ?
* ) that a transverse velocity gradient is intrinsic to the jet structure ( as discussed by * ? ? ?
* ) and that a further velocity decrease due to mass loading in the jet from the ism is likely to be present and possibly dominates the jet velocity structure at larger distances from the core .
we have strong indications that in mkn 501 the radio jet emission is oriented at a large angle with respect to the line of sight relative to that expected for a blazar :    * from the jet morphology : the jet is limb - brightened in high resolution images at @xmath78mas from the core .
the different brightness can be related to a different doppler factor only with an orientation angle @xmath40 15@xmath15 .
then , we can reproduce the observed ratio between the brightness of the external shear and the inner spine with @xmath79 , assuming @xmath80 ( see tables  6 , 7 ) . with these values , the spine doppler factor ( @xmath81 ) is lower than the external layer doppler factor ( @xmath82 ) . on the contrary , a jet orientation at a smaller angle and/or a lower spine velocity imply @xmath83 and therefore a centrally peaked structure .
* from the core dominance : the core boosting is high but not extreme ; high bulk velocities require , therefore , an orientation angle @xmath40 in the range 10@xmath1527@xmath15 with respect to the line of sight .
* from the fit to the trend of the jet brightness and fwhm : assuming an adiabatic model and consistency with the jet / counterjet brightness ratio , we need a radio jet starting with a high velocity ( @xmath840.998 ) , and in the case of a parallel magnetic field , @xmath40 has to be @xmath1125@xmath15 at @xmath3570mas .    these results seem to be in disagreement with constraints derived from the high frequency ( @xmath0-ray ) emission and variability detected in blazars . studying the rapid burst of tev photons in mkn 421 , @xcite determine that a jet with a bulk velocity @xmath8510 with opening angle and orientation angle @xmath86 is required .
similarly , in mkn  501 the expected doppler factor should be 10 , as obtained from fits of synchrotron self - compton ( ssc ) models @xcite . to reconcile the radio and high frequency results we have to assume that the inner region of the jet ( 0.0010.03pc ) is moving with a lorentz factor 10 and oriented at @xmath87 . from this region
we have high energy radiation and negligible radio emission . over the range from 0.03 to 50pc
( projected distance ) the jet orientation has to change from @xmath115@xmath15 to 25@xmath15 ( or at least @xmath88 in the case of perpendicular magnetic field ) .
we propose that this large change in the jet direction may be gradual and we derive the minimum change requested by observational data at different distances from the core . as shown in table  6 , a progressive change in the jet direction with respect to the line of sight is possible and in agreement with the observed pc - scale jet properties .    at present , there is no obvious explanation for this change in the jet direction in the inner region .
it is possible that it could be related to helical type motion @xcite ( although helical motion is not visible in our high resolution images ) , or to an instability in the jet due to reasons that are thus far unknown .
it is also important to stress that small changes in direction are not uncommon in fr  i radio galaxies and are magnified when projection effects occur . in any case
, we recall that beyond 30mas from the core the jet position angle becomes stable and closely aligned with the kpc - scale structure .
we note that in the case of the perpendicular magnetic field the jet velocity is still very high at large distances from the core , in agreement with the lack of detection of a visible counterjet in all published maps of this source and despite the symmetric structure visible in the vla low resolution image @xcite .
however , we recall that @xmath89 is required to justify the observed limb - brightened structure from doppler boosting alone , even in the case of perpendicular magnetic field .      from the present data
we can derive that :    * the nuclear region shows a self - absorbed spectrum with a peak frequency @xmath90ghz and a peak flux density @xmath91jy .
following @xcite , we can use these data to estimate the magnetic field in the nuclear region ( from 0.03 to 0.15pc ) : + b = @xmath92 + assuming @xmath40 ( the core angular size ) = 0.2mas ( 0.15pc ) and @xmath93=5 ( table  6 ) , we estimate a magnetic field b@xmath110.03  gauss . on the other hand ,
fits of ssc models to a high energy state observed in april 1997 yield values of the magnetic field in the range 0.20.02  gauss in the inner region ( @xmath160.03pc ) @xcite .
this suggests a constant or a not very large decrease of the magnetic field from the gamma ray emitting region to the beginning of the radio jet . *
the radio spectra of jet sub - components show an evident steepening at high frequencies , suggesting the effect of radiative losses in a high magnetic field . * in the region where the jet is well resolved
, the spectrum is flat in the inner spine confirming the presence of a higher transport efficiency , lower radiative losses , or reacceleration mechanisms .
the transverse steepening from the spine to the shear layer regions confirms the larger losses in these regions , probably related to the interaction with the external medium . *
it is not straightforward to determine the equipartition magnetic field @xcite given the complex structure of the jet .
we can obtain a rough estimate by assuming that electrons and heavy particles carry the same amount of energy ( @xmath94 ) , the filling factor @xmath95=1 , and the frequency range is between 10mhz and 100ghz .
this yields @xmath110.015  gauss in the c4 and c3 regions and 0.01  gauss in the c2 and c1 regions . * assuming a viewing angle @xmath40 as in tables  6 and 7 , the observed angular distances correspond to the following de - projected linear distances : +  the inner limb - brightened region is present at @xmath164pc from the core +  the first jet region extends to @xmath1110mas , corresponding to @xmath1130pc +  the vlbi jet is visible up to 170270pc ( assuming @xmath40 = 25@xmath15  15@xmath15 , respectively ) +  the kpc - scale structure is symmetric ( i.e. not relativistic ) at @xmath96kpc from the core +  the large scale symmetric structure extending to @xmath1170@xmath97 corresponds to a linear size of 120200kpc , in agreement with the expected linear size of a fr i radio galaxy with the same low frequency total radio power .
as a result of the many frequencies and resolutions available , we have established that a limb - brightened structure is indeed present in markarian 501 , beginning in the very inner jet .
this has important consequences : ( 1 ) it implies the presence of velocity structure already at very small scale , which jet models will have to take into account ; ( 2 ) it suggests , together with other arguments , that the vlbi jet can not be oriented at too small a viewing angle .
in particular , the detection of a very long one - sided jet and the results of an adiabatic expansion fit yield values of the angle of view in the range @xmath98 , depending also on the magnetic field orientation .
the bulk velocity is strongly relativistic at the beginning ( @xmath99 ) , decreasing to @xmath100 ( @xmath1011.25 ) at the end of the vlbi jet , 100mas from the core ( parallel @xmath102 ) or to @xmath103 ( @xmath104 ) in the case of a perpendicular @xmath102 .
the data considered here are consistent with a pattern speed of zero . previous reports , based mostly on lower frequency observations , have reported low pattern speeds , however the complex limb - brightened structure revealed at 22ghz makes measurement of component speeds problematic at high frequencies .
finally , the study of the spectral indices has revealed a self - absorbed compact core and a complex spectral index distribution in the jet with flat and very steep regions . a clear connection is present between limb - brightening in total intensity and the spectral index images : the inner spine has a flat spectrum while the jet more external regions show steep spectra .
our derived magnetic field strengths are similar to those obtained from data at other wavelengths , yielding values of order @xmath105  gauss in the core region .
the estimated equipartition field in the jet region is in the range 0.0150.01  gauss .
the authors thanks r.  fanti for his critical reading and useful comments on the paper .
mg thanks jive ( joint institute for vlbi in europe ) for support and hospitality and d.  gabuzda for helpful advice in the reduction of the third epoch of vsop data .
we gratefully acknowledge the vsop project , which is led by the japanese institute of space and astronautical science in cooperation with many organizations and radio telescopes around the world .
the national radio astronomy observatory is operated by associated universities , inc .
, under cooperative agreement with the national science foundation .
this research has made use of nasa s astrophysics data system bibliographic services and of the united states naval observatory ( usno ) radio reference frame image database ( rrfid ) .
this research has made use of data from the university of michigan radio astronomy observatory which is supported by funds from the university of michigan .
this material is based in part upon work supported by the italian ministry for university and research ( miur ) under grant cofin 2001 - 02 - 8773 and by the u.s .
national science foundation under grant no .
0098579 .
aaron , s.  1999 , asp conf .
ser .  159 : bl lac phenomenon , 427 baum , s.  a.  et al .  1997 , , 483 , 178 bradbury , s.  m.  et al .
1997 , , 320 , l5 cassaro , p. , stanghellini , c. , bondi , m. , dallacasa , d. , della ceca , r. , & zappal ` a , r.  a.1999 , , 139 , 601 chiaberge , m. , celotti , a. , capetti , a. , & ghisellini , g.  2000 , , 358 , 104 conway , j.  e. & wrobel , j.  m. 1995 , , 439 , 98 de ruiter , h.  r. , parma , p. , fanti , c. , & fanti , r.  1990 , , 227 , 351 cotton , w.  d. , feretti , l. , giovannini , g. , lara , l. , & venturi , t.  1999 , , 519 , 108 edwards , p.  g. , giovannini , g. , cotton , w.  d. , feretti , l. , fujisawa , k. , hirabayashi , h. , lara , l. , & venturi , t.  2000 , , 52 , 1015 edwards , p.  g.  & piner , b.  g.  2002 , , 579 , l67 feretti , l. , giovannini , g. , gregorini , l. , parma , p. , & zamorani , g. 1984 , , 139 , 55 feretti , l. , perley , r. , giovannini , g. , & andernach , h.  1999 , , 341 , 29 fey , a.  l.  & charlot , p.  1997
, , 111 , 95 ficarra , a. , grueff , g. , & tomassetti , g. 1985 , , 59 , 255 giovannini , g. , feretti , l. , gregorini , l. , & parma , p. 1988
, , 199 , 73 giovannini , g. 1999 in _ bl lac phenomenon _ , asp conf .
ser . , 159 , 439 giovannini , g. , cotton , w. d. , feretti , l. , lara , l. , & venturi , t. 2000 , adspr , 26 , 693 giovannini , g. , cotton , w.d . , feretti , l. et al .
2001 , , 552 , 508 giovannini , g. 2002 , newar in press ( astro - ph/0302176 ) hirabayashi , h.  et al.2000 , , 52 , 955 hirabayashi , h.  et al.1998 , science , 281 , 1825 and erratum 282 , 1995 homan , d.  c. , ojha , r. , wardle , j.  f.  c. , roberts , d.  h. , aller , m.  f. , aller , h.  d. , & hughes , p.  a.  2001 , , 549 , 840 junor , w. , biretta , j.  a. , & livio , m.  1999 , , 401 , 891 katarzy ' nski , k. , sol , h. , & kus , a.  2001 , , 367 , 809 kellermann , k.  i. , vermeulen , r.  c. , zensus , j.  a. , & cohen , m.  h.  1998 , , 115 , 1295 killeen , n.  e.  b. , bicknell , g.  v. , & carter , d.  1986 , , 309 , 45 laing , r.  a. , parma , p. , de ruiter , h.  r. , & fanti , r.  1999 , , 306 , 513 laing , r.  a.  & bridle , a.  h.2002a , , 336 , 328 laing , r.  a.  & bridle , a.  h.2002b , , 336 , 1161 ma , c.  et al .
1998 , , 116 , 516 marscher , a.  p.  1987 , superluminal radio sources , 280 marscher , a. p. 1999 , aph , 11 , 19 meier , d.l .
2002 , newar in press pacholczyk , a. g. 1970 , radio astrophysics ( san francisco : freeman & co. ) pollack , l.  k. , taylor , g.  b. , & zavala , r.  t.  2003 , , 589 , 733 quinn , j.  et al .  1996 , , 456 , l83 rieger , f.  m.  & mannheim , k.  2003 , , 397 , 121 salvati , m. , spada , m. , & pacini , f.  1998 , , 495 , l19 tavecchio , f. , maraschi , l. , pian , e. et al . 2001 , apj , 554 , 725 ulvestad , j. s. , johnston , k. j. & weiler , k. w. 1983 , apj , 266 , 18 urry , c. m. & padovani , p. 1995 , , 107
, 803 van breugel , w. & schilizzi r. 1986 , , 301 , 834 venturi , t.  et al .
2001 , , 379 , 755 villata , m.  & raiteri , c.  m.1999 , , 347 , 30      1995 apr 7 & 15 & vlba & 40 + 1995 apr 12 & 8.3 & vlba & 15 + 1995 dec 15 & 15 & vlba & 40 + 1996 apr 23 & 8.3 & vlba & 7.5 + 1996 apr 23 & 15 & vlba & 6 + 1996 jun 7 & 8.1 & vlba & 6 + 1997 mar 13 & 15 & vlba & 40 + 1997 apr 25 & 15 & vlba & 90 + 1997 apr 25 & 22 & vlba & 90 + 1997 may 26 & 15 & vlba & 50 + 1997 may 26 & 22 & vlba & 50 + 1997 aug 4 & 1.6 & vlba@xmath4sc+go+*halca * & 420 + 1997 aug 15 & 15 & vlba & 40 + 1997 aug 15 & 22 & vlba & 40 + 1998 apr 7 & 4.8 & vlba+eb+*halca * & 780 + 1998 apr 8 & 1.6 & vlba+ro+go+*halca * & 600 + 1998 jun 24 & 8.4 & vlba+gc+gn+kk+ & 20 + & & mc+on+wf & + 1999 jul 19 & 15 & vlba@xmath4sc & 40 + 2001 mar 5 & 1.6 & vlba+go+ro+*halca * & 480 + 2001 mar 6 & 4.8 & vlba@xmath4hn+*halca * & 540 +      1997 aug 4 & 1.6 & 7 & nz ( @xmath106 ) & go ( @xmath107 ) + 1998 apr 7 & 4.8 & 13 & rz ( @xmath108+@xmath109 ) , nz ( @xmath110 ) & eb + 1998 apr 8 & 1.6 & 10 & rz ( @xmath111 ) , nz ( @xmath107 ) & go , ro + 2001 mar 5 & 1.6 & 8 & rz ( @xmath112 ) , nz ( @xmath113 ) & ro ( @xmath107 ) , go ( @xmath114 ) + 2001 mar 6 & 4.8 & 9 & rz ( @xmath115 ) , nz ( @xmath112 ) & +      1995 apr 7 & 15 & 0.6 @xmath9 0.9 & 0.20 & 485 + & & 1.2 & 0.25 & 524 + 1995 apr 12 & 8.3 & 1.2 & 0.40 & 524 + 1995 dec 15 & 15 & 0.6 @xmath9 0.9 & 0.35 & 508 + & & 1.2 & 0.30 & 546 + 1996 apr 23 & 8.3 & 1.2 & 0.25 & 521 + & 15 & 0.6 @xmath9 0.9 & 0.30 & 490 + & & 1.2 & 0.35 & 520 + 1996 jun 7 & 8.1 & 1.2 & 0.40 & 535 + 1997 mar 13 & 15 & 0.6 @xmath9 0.9 & 0.15 & 381 + & & 1.2 & 0.20 & 410 + 1997 apr 25 & 15 & 0.6 @xmath9 0.9 & 0.10 & 360 + & & 1.2 & 0.15 & 390 + & 22 & 0.3 @xmath9 0.5 & 0.15 & 327 + 1997 may 26 & 15 & 0.6 @xmath9 0.9 & 0.10 & 408 + & & 1.2 & 0.10 & 438 + & 22 & 0.4 @xmath9 0.6 & 0.20 & 334 + 1997 aug 4 & 1.6 & 1.5 @xmath9 2.9 & 2.3 & 409 + 1997 aug 15 & 15 & 0.6 @xmath9 0.9 & 0.20 & 455 + & & 1.2 & 0.15 & 487 + & 22 & 0.3 @xmath9 0.5 & 0.25 & 350 + 1998 apr 7 & 4.8 & 0.6 @xmath9 0.9 & 0.40 & 465 + & & 1.2 & 0.25 & 509 + 1998 apr 8 & 1.6 & 1.5 @xmath9 2.9 & 0.12 & 421 + 1998 jun 24 & 8.6 & 0.6 @xmath9 0.9 & 0.25 & 477 + & & 1.2 & 0.35 & 508 + 1999 jul 19 & 15 & 0.6 @xmath9 0.9 & 0.20 & 446 + & & 1.2 & 0.20 & 480 + 2001 mar 5 & 1.6 & 1.5 @xmath9 2.9 & 0.50 & 316 + 2001 mar 6 & 4.8 & 1.2 & 0.50 & 428      1995.29 & c4 & 0.6 + & c3 & 2.4 + & c2 & 4.3 + & c1 & 7.6 + 1995.96 & c4 & 0.7 + & c3 & 2.6 + & c2 & 4.0 + & c1a & 7.4 + & c1b & 7.5 + 1996.29 & c4 & 0.5 + & c3 & 2.6 + & c2 & 4.0 + & c1a & 7.7 + & c1b & 7.8 + 1997.21 & c4 & 0.6 + & c3 & 2.6 + & c2a & 3.7 + & c2b & 3.9 + & c1a & 7.7 + & c1b & 7.9 + 1997.29 & c4 & 0.6 + & c3 & 2.3 + & c2a & 3.7 + & c2b & 4.0 + & c1 & 7.8 + 1997.37 & c4 & 0.6 + & c3 & 2.5 + & c2a & 3.8 + & c2b & 4.1 + & c1a & 7.7 + & c1b & 7.9 + 1997.62 & c4 & 0.7 + & c3 & 2.4 + & c2a & 3.6 + & c2b & 4.4 + & c1a & 7.3 + & c1b & 8.0 + 1998.48 & c4 & 0.8 + & c3 & 2.7 + & c2 & 4.7 + & c1 & 7.4 + 1999.55 & c4 & 0.7 + & c3 & 2.6 + & c2 & 4.5 + & c1a & 7.8 + & c1b & 8.9          0.0001  @xmath160.03 & 4 & 15 & 15 & ? & ? & @xmath0-ray region + 0.03  0.15 & 10 & 15 & 4 & 10 & 5 & radio core + 0.15  7 & 15 & 15 & 2 & 3 & 5 & first jet region + 7  20 & 1520 & 15 & 21 & 3 & 43 & before of large bending + 2030 & 25 & 103 & 12 & 2 & 2.5 & after the large bending + 50 & 25 & 1.25 & 1.8 & 1.1 & 1.5 & final vlbi jet region      0.0001  @xmath160.03 & 4 & 15 & 15 & ? & ? & @xmath0-ray region + 0.03  0.15 & 10 & 15 & 4 & 10 & 5 & radio core + 0.15  7 & 15 & 15 & 2 & 3 & 5 & first jet region + 7  20 & 15 & 15 & 2 & 3 & 4 & before of large bending + 2030 & 15 & 10 & 2.5 & 3 & 4 & after the large bending + 50 & 15 & 10 & 2.5 & 3 & 4 & final vlbi jet region