the pion distribution amplitude ( da ) is an important element for applying the qcd factorization theory and qcd light - cone sum rules ( lcsr ) to exclusive processes involving pion . for examples ,
it is important for understanding the semi - leptonic decays @xmath5 and @xmath6 , the pion - photon transition form factor ( tff ) @xmath7 , the exclusive process @xmath8 , and etc .. inversely , one can make use of all those processes to obtain stringent constraints on the pion da  @xcite .
because of its simplest structure , the pion leading - twist da has attracted much attention after the pioneering works done by refs.@xcite .
unfortunately , at present , there is no definite conclusion on the behavior of pion leading - twist da .
for example , the babar measurement for the pion - photon tff  @xcite supports the chernyak - zhitnitsky ( cz)-like  @xcite behavior ; while the corresponding data by the belle collaboration  @xcite supports the asymptotic - like behavior  @xcite .    the pion leading - twist da at the scale @xmath9 can be expanded into a gegenbauer polynomial as  @xcite , @xmath10 , \label{dagegexp}\end{aligned}\ ] ] where the @xmath11 are gegenbauer polynomials and the @xmath12 are gegenbauer moments .
theoretically , the gegenbauer moments have been calculated by qcd sum rules  @xcite and the lattice gauge theory  @xcite .
phenomenologically , the gegenbauer moments have also been estimated from the pion - photon transition form factor  @xcite , the pion electromagnetic form factor  @xcite , the @xmath13 tffs  @xcite , and etc .
most of those estimations are consistent with each other within errors .
it is noted that the central values of them are different from each other , especially , the lattice qcd estimations are bigger than those of qcd sum rules . due to the forthcoming more accurate data ,
it is helpful to provide a more accurate theoretical prediction on the pion da for a better comparison with the data .    since its invention  @xcite
, the shifman - vainshtein - zakharov ( svz ) sum rules has been widely adopted for dealing with the hadron phenomenology .
in addition , the background field theory provides a systematic description for its key components , i.e. the vacuum condensates , from the viewpoint of qcd field theory  @xcite .
it assumes that the quark and gluon fields are composed of the background fields and the quantum fluctuations around them .
thus , to take the background field theory as the foundation for the qcd sum rules , it not only has distinct physical picture , but also greatly simplifies the calculation due to its capability of adopting different gauges for quantum fluctuations and background fields .
the qcd sum rules within the background field theory has already been adopted to study the meson properties , cf.refs.@xcite .
previous sum rule estimations for the pion da moments are usually done up to dimension - four condensates , or including part of the dimension - six condensates contributions . in the present paper
, we shall adopt the svz sum rules under the background field theory to make a detailed study on the pion leading - twist da up to dimension - six operators .
for the purpose , as the first time , we shall expand both the quark propagator and the vertex operator @xmath0 up to a more complex form with full dimension - six operators contributions so as to achieve an accurate qcd sum rule estimation .
we shall show that those dimension - six operators shall result in new dimension - six condensates that do provide sizable contributions and should be taken into consideration for a sound estimation , i.e. they shall bring sizable changes in our previous predictions on the pion da behavior .
the remaining parts of the paper are organized as follows . in section
ii , we first present a brief introduction of the qcd sum rules and then provide the quark propagator and the vertex @xmath0 up to dimension - six operators under the background field theory .
in section iii , we present the calculation technology for deriving the sum rules for pion leading - twist da moments .
numerical results are given in section iv .
section v is reserved for a summary . for convenience , formulas for the propagators and the vacuum matrix elements in @xmath14-dimension space
are presented in the appendix .
within the svz sum rules , the hadrons are represented by their interpolating quark currents taken at large virtualities .
the correlator of those currents can be treated within the operator product expansion , whose short - distance coefficients are calculated using qcd perturbation theory , whereas its long - distance parts are written according to the non - perturbative but universal vacuum condensates , such as the quark condensate @xmath15 , the gluon condensate @xmath16 , and etc .
furthermore , under the background field theory , the quark and gluon fields are composed of background fields and quantum fluctuations around them . because of the influence from background fields , the quark and gluon propagators shall include nonperturbative component inevitably , which however can be treated via a systematic way .    within the background field theory , the gluon field @xmath17 and quark field @xmath18 in the qcd lagrangian or green function should be replaced by @xmath19 where in the right hand side of the arrow , @xmath17 with @xmath20 and @xmath18 indicate the background fields of the gluon and quark , respectively . @xmath21 and @xmath22 stand for the gluon and quark s quantum field , i.e. , the quantum fluctuation on the background fields .
the lagrangian with the background field theory can be found in ref.@xcite .
usually , the background quark and gluon fields satisfy the following equations of motion @xmath23 where @xmath24 and @xmath25 indicate the fundamental and adjoint representations of the gauge covariant derivatives respectively . as an advantage of the background field theory
, one can take different gauges for dealing with the quantum fluctuations and background fields .
more specifically , one can adopt the  background gauge "  @xcite @xmath26 for the gluon quantum field , the schwinger gauge or the  fixed - point gauge "  @xcite @xmath27 for the background field .          within the framework of the background field theory , the quark propagator
would be affected by the background quark and/or gluon fields , which satisfies the equation @xmath28 if taking @xmath29 eq.([quaprofun1 ] ) can be changed as @xmath30 where @xmath31 , and @xmath32    moreover , using the  fixed - point gauge " as shown in eq.([fpgau ] ) , the gluon background field can be expressed by the gauge invariant @xmath33 as    @xmath34    where @xmath33 is a notation for @xmath35 , where the indexes @xmath36 indicate the covariant derivative up to @xmath37-th order .
@xmath38 with @xmath39 stands for the gluon field strength tensor and @xmath40 is the structure constant of @xmath41 group .    substituting eq.([gluexp ] ) into eq.([quaprofun3 ] ) , we obtain the expressions for @xmath42 . by further using eq.([quaprofun2 ] ) , we can obtain the required quark propagator in the background field , which can be expressed as    @xmath43    the quark propagators with various gauge invariant tensors @xmath33 that shall result in up to dimension - six operators are presented in appendix a.
the feynman diagrams for the quark propagators ( [ appa1],@xmath44,[prop2 ] ) that with various gauge invariant tensors are shown in fig.([propagator ] ) , where thirteen figures , i.e. fig.([propagator]a ) , @xmath44 , fig.([propagator]m ) , correspond to @xmath45 , @xmath44 , @xmath46 , respectively .
the cross ( @xmath47 ) attached to the gluon line indicates the tensor of the local gluon background field with
@xmath37 " stands for @xmath37-th order covariant derivative .
because the  fixed - point gauge " , as indicated by eq.([fpgau ] ) , violates the translation invariance , the quark propagator from @xmath48 to @xmath49 , @xmath50 , can not be directly obtained by applying the replacement @xmath51 in eq.([prop4 ] ) .
it can be related with @xmath52 via the relation  @xcite @xmath53 where @xmath54 stands for the charge conjugation matrix and the symbol @xmath55 indicates transpose of both the dirac and the color matrices .
when applying the svz sum rules to calculate the moments of meson , one would encounter the vertex operator , @xmath56 , where @xmath57 indicates some kind of dirac matrices , e.g. , @xmath58 for the pion leading - twist da , @xmath59 or @xmath60 for the pion twist-3 das , and etc .. we have @xmath61 where the @xmath44 stands for the higher - order terms which are irrelevant for our present analysis and @xmath62    we can expand the operator @xmath0 into series with operators @xmath63 and @xmath64 . for the purpose ,
we first expand @xmath0 as @xmath65 \nonumber\\ & & + 3 ( z\cdot \tensor{\partial } ) ( z\cdot \underline{b})^2 + ( z\cdot b)^3 , \nonumber\\ & & \cdots \cdots , \nonumber\end{aligned}\ ] ] where the  underline " below @xmath66 ( or the latter @xmath48 ) indicates that the operation @xmath67 does not act on it . in deriving those equations ,
the following equation has been adopted , @xmath68.\ ] ] by keeping only those terms that shall leads to operators up to dimension - six , we obtain @xmath69 g_{\mu\nu;\rho\sigma } , \nonumber\\ ( z\cdot \tensor{d})^n_{5(1 ) } & = & -\frac{n(n-1)}{3 } ( z\cdot \tensor{\partial})^{n-2 } \nonumber\\ & & \times \underline{x}^\mu \underline{x}^\rho \underline{x}^\lambda z^\nu z^\sigma g_{\mu\nu } g_{\rho\sigma;\lambda } , \nonumber\\ ( z\cdot \tensor{d})^n_{5(2 ) } & = & -\frac{n(n-1)}{3 } ( z\cdot \tensor{\partial})^{n-2 } \nonumber\\ & & \times \underline{x}^\mu \underline{x}^\rho \underline{x}^\lambda z^\nu z^\sigma g_{\mu\nu;\lambda } g_{\rho\sigma } , \nonumber\\ ( z\cdot \tensor{d})^n_{5(3 ) } & = & \left [ -\frac{i}{15 } \times n ( z\cdot \tensor{\partial})^{n-1 } \underline{x}^\mu \underline{x}^\rho \underline{x}^\sigma \underline{x}^\lambda z^\nu \right .
\nonumber\\ & & - \frac{i}{45 } n(n-1)(n-2 ) ( z\cdot \tensor{\partial})^{n-3 } \nonumber\\ & & \left .
\times \underline{x}^\mu \left ( \underline{x}^\rho z^\sigma z^\lambda + z^\rho \underline{x}^\sigma z^\lambda + z^\rho z^\sigma \underline{x}^\lambda \right ) z^\nu \right ] \nonumber\\ & & \times g_{\mu\nu;\rho\sigma\lambda } , \nonumber\\ ( z\cdot \tensor{d})^n_{6(1 ) } & = & \frac{i}{6 } n(n-1)(n-2 ) ( z\cdot \tensor{\partial})^{n-3 } \nonumber\\ & & \times \underline{x}^\mu \underline{x}^\rho \underline{x}^\lambda z^\nu z^\sigma z^\tau g_{\mu\nu } g_{\rho\sigma } g_{\lambda\tau } , \nonumber\\ ( z\cdot \tensor{d})^n_{6(2 ) } & = & \left [ -\frac{n(n-1)}{8 } ( z\cdot \tensor{\partial})^{n-2 } \underline{x}^\mu \underline{x}^\rho \underline{x}^\lambda \underline{x}^\tau z^\nu z^\sigma \right .
\nonumber\\ & & - \frac{n(n-1)(n-2)(n-3)}{12 } ( z\cdot \tensor{\partial})^{n-4 } \nonumber\\ & & \left .
\times \underline{x}^\mu \underline{x}^\rho z^\lambda z^\tau
z^\nu z^\sigma \right ] g_{\mu\nu } g_{\rho\sigma;\lambda\tau } , \nonumber\\ ( z\cdot \tensor{d})^n_{6(3 ) } & = & -\frac{2}{9 } n(n-1 ) ( z\cdot \tensor{\partial})^{n-2 } \nonumber\\ & & \times \underline{x}^\mu \underline{x}^\rho \underline{x}^\lambda \underline{x}^\tau
z^\nu z^\sigma g_{\mu\nu;\lambda } g_{\rho\sigma;\tau } , \nonumber\\ ( z\cdot \tensor{d})^n_{6(4 ) } & = & -\frac{n(n-1)}{8 } ( z\cdot \tensor{\partial})^{n-2 } \nonumber\\ & & \times \underline{x}^\mu \underline{x}^\rho \underline{x}^\lambda \underline{x}^\tau
z^\nu z^\sigma g_{\mu\nu;\lambda\tau } g_{\rho\sigma } , \nonumber\\ ( z\cdot \tensor{d})^n_{6(5 ) } & = & \left [ -\frac{i}{72 } \times n ( z\cdot \tensor{\partial})^{n-1 } \underline{x}^\mu \underline{x}^\rho \underline{x}^\sigma \underline{x}^\lambda \underline{x}^\tau z^\nu \right .
\nonumber\\ & & - \frac{i}{216 } n(n-1)(n-2 ) ( z\cdot \tensor{\partial})^{n-3 } \nonumber\\ & & \times \underline{x}^\mu \left ( z^\rho z^\sigma \underline{x}^\lambda \underline{x}^\tau + z^\rho \underline{x}^\sigma z^\lambda \underline{x}^\tau + z^\rho \underline{x}^\sigma \underline{x}^\lambda z^\tau \right . \nonumber\\ & & \left
. + \underline{x}^\rho z^\sigma z^\lambda \underline{x}^\tau + \underline{x}^\rho z^\sigma \underline{x}^\lambda z^\tau + \underline{x}^\rho \underline{x}^\sigma z^\lambda z^\tau \right ) z^\nu \nonumber\\ & & - \frac{i}{360 } n(n-1)(n-2)(n-3)(n-4 ) \nonumber\\ & & \left .
\times ( z\cdot \tensor{\partial})^{n-5 } \underline{x}^\mu z^\rho z^\sigma z^\lambda z^\tau z^\nu \right ] g_{\mu\nu;\rho\sigma\lambda\tau } , \nonumber\end{aligned}\ ] ] where the subscript @xmath70 with @xmath71 stands for the dimension of the operator , in which @xmath72 stands for the @xmath73-th type of the operator with same dimension . for example , there is two type of dimension four operators , three type of dimension five operators and five type of dimension - six operators .
fig.([vertex ] ) represents the feynman diagrams for the vertex operator @xmath74 under the background field theory , where thirteen figures , i.e. fig.([vertex]a ) , @xmath44 , fig.([vertex]m ) , correspond to @xmath75 , @xmath44 , @xmath76 , respectively .
the cross ( @xmath47 ) attached to the gluon line indicates the tensor of the local gluon background field with
@xmath37 " stands for the @xmath37-th order covariant derivative .
stands for the background quark field . the left big dot and the right big dot stand for the vertex operators @xmath77 and @xmath78 in the correlator , respectively .
the full feynman diagrams can be found in figs.([feyna],[feynb ] ) .
, scaledwidth=45.0% ]            considering the definition @xmath79 where @xmath80 is the pion decay constant and the @xmath37-th order moments of the pion da are defined as @xmath81 especially , the @xmath49-th moment satisfies the normalization condition @xmath82    to derive the sum rules for the pion leading - twist da moments @xmath83 , we introduce the following correlation function ( correlator ) , @xmath84 where @xmath85 , @xmath86 and the currents @xmath87    in physical region , the correlator ( [ correlator ] ) can be treated by inserting a complete set of intermediate hadronic states , which can be simplified with the help of eq.([mom ] ) as @xmath88 where the quark - hadron duality has been adopt and @xmath89 stands for the continuum threshold .    on the other hand , we apply the ope for the correlator ( [ correlator ] ) in deep euclidean region . substituting the replacement rules ( [ bfrep0],[bfrep ] ) into the correlator ( [ correlator ] ) and applying the corresponding feynman rules , we obtain @xmath90    the feynman diagrams for the correlator ( [ ope ] ) are showed in fig.([feyn ] ) , whose details are further presented in figs.([feyna],[feynb ] ) .
the solid and helical lines are for the quark and gluon propagators , respectively .
the cross @xmath91 symbol stands for the background quark field .
the left big dot and the right big dot are the vertex operators @xmath77 and @xmath92 , respectively .
fig.([feyn]a ) corresponds to the first term in eq.([ope ] ) , fig.([feyn]b ) and its permutation - diagrams correspond to the second and third terms of eq.([ope ] ) , respectively .    as for fig.([feyna ] ) , we observe    * figs.([feyna]l,[feyna]p,[feyna]q,[feyna]s ) contribute zero because their quark loop has @xmath93=0 $ ] .
the contribution from fig.([feyna]g ) is negligible due to strong suppression from small @xmath94-quark current masses . *
fig.([feyna]a ) provides the perturbative contribution , figs.([feyna]b,[feyna]r ) provide the contribution proportional to the dimension - four condensate @xmath16 , figs.([feyna]c,[feyna]h,[feyna]m,[feyna]u ) provide the contribution proportional to the dimension - six condensate @xmath95 and figs.([feyna]e,[feyna]j,[feyna]o,[feyna]t ) provide the contribution proportional to the dimension - six condensate @xmath96 . here
@xmath16 and @xmath95 are abbreviations for @xmath97 and @xmath98 . *
the remaining diagrams involves the contributions proportional to the dimension - six condensate either @xmath95 or @xmath99 , where @xmath100 .    as for fig.([feynb ] )
, we observe that    * fig.([feynb]a ) provides the contribution proportional to the dimension three condensate @xmath15 and the dimension five condensate @xmath101 .
fig.([feynb]b ) provides the contribution proportional to either the dimension five condensate @xmath101 or the dimension - six condensate @xmath102 . * the remaining three diagrams ( [ feynb]c,[feynb]d,[feynb]e )
provide the contribution proportional to the dimension - six condensate @xmath102 , where @xmath103 .
we adopt the @xmath104-scheme to deal with the infrared divergences emerged in fig.([feyna ] ) , whose divergent terms shall be absorbed into the poin da following the way suggested by ref.@xcite . during the calculation
, we shall meet with the following matrix elements : @xmath105 the formulas relating the first three vacuum matrix elements with the conventional condensates can be found in ref.@xcite , and the corresponding formulas for the remaining vacuum matrix elements under the @xmath14-dimensional space , @xmath106 , are presented in appendix b.    as a combination the correlator within different regions , the sum rules for the pion da moments can be derived by the dispersion relation @xmath107 where @xmath108 is the borel parameter and the borel transformation operator @xmath109    our final result reads    @xmath110    where @xmath111 , \nonumber\\ & & \mathbb{b } = \frac{8(165n^2 + 323n + 131)}{729(n+1 ) } + \frac{8(51n+44)}{243 } \ln
\left(\frac{m^2}{\mu^2}\right ) + \theta(n-2 ) \left\ { -\frac{4(16n^3 - 89n^2 - 39n-30)}{729n(n+1 ) } + \frac{392n}{243 } \ln \left(\frac{m^2}{\mu^2}\right ) \right . \nonumber\\ & & \quad\quad \left .
+ \sum^{n-2}_{k=0 } ( -1)^k \left [ \frac{64n(k+1)}{243(n+1)(n - k ) } - \frac{8(10k^2 - 10nk+58k-7n+48)}{243(n - k)(n - k-1 ) } + \frac{8(8k-5n-2)}{243(k+1)(n - k+1 ) } \right ] \right\}. \nonumber\end{aligned}\ ] ]    as a final remark , we can obtain the gegenbauer moments of the pion leading - twist da with the help of @xmath83 .
that is , by substituting eq.([dagegexp ] ) into eq.([nmom ] ) , we have @xmath112
to do the numerical calculation , we adopt  @xcite : @xmath113 and @xmath114 . as for the condensates ,    * we adopt @xmath115 which is derived from the current algebra  @xcite and the partially conserved axial current  @xcite . *
we adopt  @xcite , @xmath116 which is derived by using @xmath117 with @xmath118  @xcite . *
as suggested in ref.@xcite , we adopt  @xcite @xmath119 where @xmath120  @xcite . more definitely , we take  @xcite @xmath121 and @xmath122 . and
we take @xmath123 . * combining ( [ relgluqua ] ) with ( [ ratquaglu ] ) , together with the ratio @xmath124  @xcite , we obtain @xmath125 and @xmath126    the leading - order @xmath127 is fixed by @xmath128  @xcite and the renormalization scale is taken as @xmath129 . usually , the continuum threshold @xmath89 is taken as the square of the first exciting state of pion , i.e. , @xmath130 , which however may underestimate some continuum states contributions . at the present , we determine the value of @xmath89 from the sum rules of the @xmath49-th moment @xmath131 together with its normalization condition ( [ normalize0 ] ) , which leads to @xmath132 .
.the moments @xmath83 of the pion leading - twist da at the scale @xmath133 and various condensates contributions over the total dispersion integration . @xmath134 and @xmath135 are abbreviations for @xmath136 and @xmath137 . [
cols="^,^,^,^",options="header " , ]     [ daparameter ]    using those moments ( [ mom21],[mom22],[mom23 ] ) , together with the formulas ( [ gegmom1],[gegmom2],[gegmom3 ] ) , we can obtain the gegenbauer moments @xmath138 , @xmath139 and @xmath140 .
for example , we obtain @xmath141 and @xmath142 .
we present a comparison of @xmath143 under the svz sum rules and the lattice qcd in tab.[taba2 ] .
phenomenologically , a lcsr analysis on the pion electromagnetic form factor gives @xmath144  @xcite and @xmath145  @xcite .
a lcsr analysis on the pion transition form factor gives @xmath144  @xcite , @xmath146  @xcite , @xmath147  @xcite , @xmath148  @xcite and @xmath149  @xcite .
a lcsr analysis on the @xmath150 gives @xmath151  @xcite , @xmath152  @xcite and @xmath153  @xcite .
it shows that our present estimation of @xmath154 agrees with most of the estimations within errors derived in the literature .
especially , our central value is close to the lattice qcd estimations  @xcite .
following the idea of ref.@xcite , by including the correlation effect from the transverse distribution , we can construct a pion da model as    @xmath155 - \textrm{erf } \left [ \sqrt{\frac{m_q^2}{8\beta_\pi^2 x(1-x ) } } \right ] \right\ } , \label{da_model}\end{aligned}\ ] ]    where the light constitute quark @xmath156 gev , the error function @xmath157 , and the longitudinal part can be constructed as @xmath158.\ ] ] the parameters @xmath159 , @xmath160 , @xmath161 and @xmath162 can be fixed by the pion da normalization , the constraint from @xmath163  @xcite and the above determined gegenbauer moments @xmath154 and @xmath164 .
several typical values for those parameters are presented in table [ daparameter ] .
the background field theory provides an unambiguous physical picture for both the qcd physical vacuum and the svz sum rules . under this framework
, the ope of the correlator can be calculated systematically .
for the first time , we provide the quark propagator and the vertex up to dimension - six operators under the background field theory .
it is noted that those newly dimension - six terms for the quark propagator are important for a sound estimation of a physics process up to dimension - six condensates . because the mass terms are kept explicitly in the quark propagator ,
it is applicable for both the massless and massive cases .    in the present paper
, we have studied the pion leading - twist da within the qcd sum rules under the background field theory .
the resultant sum rules for the da moments up to dimension - six condensates have been presented .
it has been found that the dimension - six condensates can provide sizable contributions to the da moments in comparison to contributions from the lower dimension condensates . within errors
, our estimation of @xmath4 agrees with most of the estimations derived in the literature .
its central value is larger than those determined by previous qcd sum rules , which is mainly caused by the contribution from the dimension - six condensate @xmath137 .     by varying the second and four gegenbauer moments within their reasonable region , i.e. @xmath141 and @xmath142 .
the solid line stands for the central values for @xmath165 and the shaded band stands for the theoretical uncertainties caused by the gegenbauer moments .
, scaledwidth=50.0% ]    as an application of the suggested pion da model ( [ da_model ] ) , we redraw the pion - photon tff @xmath7 in fig.([fig_pptff ] ) .
all the formulas for @xmath7 can be found in our previous paper @xcite , only one need to change the pion da used there to be our present model .
fig.([fig_pptff ] ) indicates that in the large @xmath166 region , the pion - photon tff @xmath7 lies in between the belle data and the babar data , which is consistent with the recent @xmath167 factorization estimation with the joint resummation  @xcite .
* acknowledgments * : this work was supported in part by the natural science foundation of china under grant no.11075225 , no.11235005 and no.11275280 .
within the framework of the background field theory , the full quark propagator with various gauge invariant tensors that shall result in up to dimension - six operators can be written as    @xmath168    in which the propagators with various gauge invariant tensors are    @xmath169 \gamma^\nu g_{\mu\nu;\rho } \right\},\\ & & s_f^{4(1)}(x,0 ) = i \int \frac{d^4p}{(2\pi)^4 } e^{-ip\cdot x } \left\ { \left [ \frac{1}{4 } \left ( \frac{\gamma^\mu ( m - \not\ !
p)}{(m^2-p^2)^3 } - 2\frac{p^\mu}{(m^2-p^2)^3 } \right ) \gamma^\nu \gamma^\rho \gamma^\sigma \right.\right .
\nonumber\\ & & \quad\quad\quad\quad\ \ \left.\left .
+ \frac{1}{2 } \left ( \frac{(m + \not\ !
p ) \gamma^\mu}{(m^2-p^2)^3 } g^{\nu\sigma } + 4\frac{\gamma^\mu ( m - \not\
! p)}{(m^2-p^2)^4 }
p^\nu p^\sigma \right ) \gamma^\rho \right ] g_{\mu\nu}g_{\rho\sigma } \right\ } ,   \\ & & s_f^{4(2)}(x,0 ) = i \int \frac{d^4p}{(2\pi)^4 } e^{-ip\cdot x } \left\ { \frac{i}{4 } \left [ \frac{g^{\{\mu\rho}\gamma^{\sigma\ } } ( m - \not\ !
p)}{(m^2-p^2)^3 } - 2\frac{g^{\{\mu\rho}p^{\sigma\}}}{(m^2-p^2)^3 } + 4\frac{\gamma^{\{\mu}p^\rho p^{\sigma\ } } ( m - \not\ !
p)}{(m^2-p^2)^4 } \right ] \gamma^\nu g_{\mu\nu;\rho\sigma } \right\ } , \\ & & s_f^{5(1)}(x,0 ) = i \int \frac{d^4p}{(2\pi)^4 } e^{-ip\cdot x } \left\ { -\frac{i}{3 } \left [ \left ( 3\frac{\gamma^\mu ( m-\not\ !
p ) \gamma^\nu}{(m^2-p^2)^4 } ( p^\lambda \gamma^\rho + p^\rho \gamma^\lambda ) - \frac{\gamma^\nu ( g^{\mu\lambda } \gamma^\rho + g^{\mu\rho } \gamma^\lambda)}{(m^2-p^2)^3 } \right ) \gamma^\sigma \right.\right .
\nonumber\\ & & \quad\quad\quad\quad\ \ \left.\left .
+ 4 \left ( \frac{\gamma^\mu ( m-\not\ ! p)}{(m^2-p^2)^4 } g^{\{\nu\sigma } p^{\lambda\ } } + 2 \frac{p^\mu g^{\{\nu\sigma } p^{\lambda\}}}{(m^2-p^2)^4 } + 6\frac{\gamma^\mu ( m-\not\ ! p)}{(m^2-p^2)^5 }
p^\sigma p^\lambda \right ) \gamma^\rho \right ] g_{\mu\nu}g_{\rho\sigma;\lambda } \right\ } ,   \\ & & s_f^{5(2)}(x,0 ) = i \int \frac{d^4p}{(2\pi)^4 } e^{-ip\cdot x } \left\ { \frac{2i}{3 } \left [ \left ( \frac{g^{\mu\lambda}}{(m^2-p^2)^3 } + \frac{6p^\mu p^\lambda}{(m^2-p^2)^4 }
\right ) \gamma^\nu \gamma^\rho \gamma^\sigma \right.\right .
\nonumber\\ & & \quad\quad\quad\quad\ \ \left.\left .
- 2 \left ( \frac{\gamma^\mu ( m-\not\ ! p)}{(m^2-p^2)^4 } g^{\{\nu\sigma } p^{\lambda\ } } + 2 \frac{p^\mu g^{\{\nu\sigma } p^{\lambda\}}}{(m^2-p^2)^4 } + 6\frac{\gamma^\mu ( m-\not\ ! p)}{(m^2-p^2)^5 }
p^\sigma p^\lambda \right ) \gamma^\rho \right ] g_{\mu\nu;\lambda}g_{\rho\sigma } \right\ } ,   \\ & & s_f^{5(3)}(x,0 ) = i \int \frac{d^4p}{(2\pi)^4 } e^{-ip\cdot x } \left\ { \frac{4}{15 } \left [ \frac{g^{\{\rho\sigma } p^\lambda \gamma^{\mu\ } } ( m-\not\ ! p)}{(m^2-p^2)^4 } - 2 \frac{g^{\{\rho\sigma } p^\lambda p^{\mu\}}}{(m^2-p^2)^4 } + 6 \frac{\gamma^{\{\mu } p^\rho p^\sigma p^{\lambda\ } } ( m-\not\ !
p)}{(m^2-p^2)^5 } \right.\right . \nonumber\\ & & \quad\quad\quad\quad\ \ \left.\left .
- \frac{g^{(\mu\rho\sigma\lambda)}}{(m^2-p^2)^3 } \right ] \gamma^\nu g_{\mu\nu;\rho\sigma\lambda } \right\ } ,   \\ & & s_f^{6(1)}(x,0 ) = i \int \frac{d^4p}{(2\pi)^4 } e^{-ip\cdot x } \frac{i}{8 } \left\ { \left [ \frac{\gamma^\mu ( m - \not\ !
p)}{(m^2-p^2)^4 } - 4\frac{p^\mu}{(m^2-p^2)^4 } \right ] \gamma^\nu \gamma^\rho \gamma^\sigma \gamma^\lambda \gamma^\tau \right . \nonumber\\ & & \quad\quad\quad\quad\ \ \left .
+ 2\left [ 3\frac{\gamma^\mu ( m - \not\
! p)}{(m^2-p^2)^4 } g^{\sigma\tau } + 16\frac{\gamma^\mu ( m - \not\ ! p)}{(m^2-p^2)^5 }
p^\sigma p^\tau - 4\frac{g^{\mu\sigma}p^\tau + g^{\mu\tau}p^\sigma}{(m^2-p^2)^4 } \right ]
\gamma^\nu \gamma^\rho \gamma^\lambda \right\ } g_{\mu\nu}g_{\rho\sigma}g_{\lambda\tau } , \\ & & s_f^{6(2)}(x,0 ) = i \int \frac{d^4p}{(2\pi)^4 } e^{-ip\cdot x } \left(-\frac{1}{8}\right ) \left\ { \left [ 3\frac{\gamma^\mu ( m - \not\ !
p ) \gamma^\nu}{(m^2-p^2)^4 } g^{\{\lambda\tau } \gamma^{\rho\ } } + 16\frac{\gamma^\mu ( m - \not\ !
p ) \gamma^\nu}{(m^2-p^2)^5 } \gamma^{\{\rho } p^\lambda p^{\tau\ } } \right.\right .
\nonumber\\ & & \quad\quad\quad\quad\ \ \left .
- 4\frac{\gamma^\nu}{(m^2-p^2)^4 } g^{\mu\{\lambda } p^\tau \gamma^{\rho\ } } \right ] \gamma^\sigma + 4\left [ \frac{m + \not\ ! p}{(m^2-p^2)^4 } g^{(\nu\sigma\tau\lambda ) } + 6\frac{m+\not\ ! p}{(m^2-p^2)^5 } g^{\{\nu\sigma } p^{\tau } p^{\lambda\ } } \right .
\nonumber\\ & & \quad\quad\quad\quad\ \ \left.\left .
+ 48\frac{m + \not\
! p}{(m^2-p^2)^6 }
p^\nu p^\sigma p^\tau p^\lambda \right ] \gamma^\mu \gamma^\rho \right\ } g_{\mu\nu}g_{\rho\sigma;\lambda\tau } , \\ & & s_f^{6(3)}(x,0 ) = i \int \frac{d^4p}{(2\pi)^4 } e^{-ip\cdot x } \left ( -\frac{2}{9 } \right ) \left\ { 3\left ( \left [ 2\frac{\gamma^\mu ( m - \not\ ! p)}{(m^2-p^2)^5 }
p^\lambda p^\tau - \frac{g^{\{\mu\lambda } p^{\tau\}}}{(m^2-p^2)^4 } - \frac{4p^\mu p^\lambda p^\tau}{(m^2-p^2)^5 } \right ] \gamma^\nu \gamma^\rho \right.\right . \nonumber\\ & & \quad\quad\quad\quad\ \ \left .
+ ( \mu \leftrightarrow \lambda ) + ( \rho \leftrightarrow \tau ) + ( \mu \leftrightarrow \lambda , \rho \leftrightarrow \tau ) \right ) \gamma^\sigma + 4\left [ \frac{m + \not\ ! p}{(m^2-p^2)^4 } g^{(\nu\lambda\sigma\tau ) } \right .
\nonumber\\ & & \quad\quad\quad\quad\ \ \left.\left
. + 6\frac{m+\not\ ! p}{(m^2-p^2)^5 } g^{\{\nu\lambda } p^{\sigma } p^{\tau\ } } + 48\frac{m + \not\ ! p}{(m^2-p^2)^6 }
p^\nu p^\lambda p^\sigma p^\tau \right ] \gamma^\mu \gamma^\rho \right\ } g_{\mu\nu;\lambda}g_{\rho\sigma;\tau } , \\ & & s_f^{6(4)}(x,0 ) = i \int \frac{d^4p}{(2\pi)^4 } e^{-ip\cdot x } \left ( -\frac{1}{2 } \right ) \left\ { \left [ \frac{m + \not\ !
p}{(m^2-p^2)^4 } g^{(\nu\tau\lambda\sigma ) } + 6\frac{m+\not\ !
p}{(m^2-p^2)^5 } g^{\{\nu\tau } p^{\lambda } p^{\sigma\ } } \right.\right.\nonumber\\ & & \quad\quad\quad\quad\ \ \left.\left .
+ 48\frac{m + \not\
! p}{(m^2-p^2)^6 }
p^\nu p^\tau p^\lambda p^\sigma \right ] \gamma^\mu \gamma^\rho - 3\left [ \frac{g^{\{\mu\lambda } p^{\tau\}}}{(m^2-p^2)^4 } + \frac{8p^\mu p^\lambda p^\tau}{(m^2-p^2)^5 } \right ] \gamma^\nu \gamma^\rho \gamma^\sigma \right\ } g_{\mu\nu;\lambda\tau}g_{\rho\sigma } , \\ & & s_f^{6(5)}(x,0 ) = i \int \frac{d^4p}{(2\pi)^4 } e^{-ip\cdot x } \left\ { -\frac{i}{18 } \left [ \frac{g^{[\rho\sigma\lambda\tau } \gamma^{\mu ] } ( m-\not\ ! p)}{(m^2-p^2)^4 } - 4 \frac{g^{[\rho\sigma\lambda\tau } p^{\mu]}}{(m^2-p^2)^4 } + 6 \frac{g^{\{\rho\sigma } p^\lambda p^\tau \gamma^{\mu\ } } ( m-\not\ ! p)}{(m^2-p^2)^5 } \right.\right .
\nonumber\\ & & \quad\quad\quad\quad\ \ \left.\left .
- 12 \frac{g^{\{\rho\sigma } p^\lambda p^\tau p^{\mu\}}}{(m^2-p^2)^5 } + 48 \frac{p^{\{\rho } p^\sigma p^\lambda p^\tau \gamma^{\mu\ } } ( m-\not\!p)}{(m^2-p^2)^6 } \right ] \gamma^\nu g_{\mu\nu;\rho\sigma\lambda\tau } \right\ } , \label{prop2}\end{aligned}\ ] ]    where    @xmath170 } & = & g^{(\mu\nu\rho\sigma ) } p^\lambda + g^{(\lambda\nu\rho\sigma ) } p^\mu + g^{(\mu\lambda\rho\sigma ) } p^\nu + g^{(\mu\nu\lambda\sigma ) } p^\rho + g^{(\mu\nu\rho\lambda ) } p^\sigma , \nonumber\\ g^{\{\mu\nu } p^{\rho\ } } & = & g^{\mu\nu } p^\rho + g^{\mu\rho } p^\nu + g^{\nu\rho } p^\mu , \nonumber\\ g^{\{\mu\nu } p^\rho p^{\sigma\ } } & = & g^{\mu\nu } p^\rho p^\sigma + g^{\mu\rho } p^\nu p^\sigma + g^{\mu\sigma } p^\nu p^\rho + g^{\nu\rho } p^\mu p^\sigma + g^{\nu\sigma } p^\mu p^\rho + g^{\rho\sigma } p^\mu p^\nu , \nonumber\\ g^{\{\mu\nu } p^\rho p^\sigma p^{\lambda\ } } & = & g^{\mu\nu } p^\rho p^\sigma p^\lambda + g^{\mu\rho } p^\nu p^\sigma p^\lambda + g^{\mu\sigma } p^\rho p^\nu p^\lambda + g^{\mu\lambda } p^\rho p^\sigma p^\nu + g^{\nu\rho } p^\mu p^\sigma p^\lambda + g^{\nu\sigma } p^\rho p^\mu p^\lambda \nonumber\\ & & + g^{\nu\lambda } p^\rho p^\sigma p^\mu + g^{\rho\sigma } p^\mu p^\nu p^\lambda + g^{\rho\lambda } p^\mu p^\sigma p^\nu + g^{\sigma\lambda }
p^\rho p^\mu p^\nu \nonumber\\ g^{\{\mu\rho } \gamma^{\sigma\ } } & = & g^{\mu\rho } \gamma^\sigma + g^{\mu\sigma } \gamma^\rho + g^{\rho\sigma } \gamma^{\mu } , \nonumber\\ g^{[\rho\sigma\lambda\tau } \gamma^{\mu ] } & = & \gamma^\mu g^{(\rho\sigma\lambda\tau ) } + ( \mu \leftrightarrow \rho ) + ( \mu \leftrightarrow \sigma ) + ( \mu \leftrightarrow \lambda ) + ( \mu \leftrightarrow \tau ) , \nonumber\\ g^{\{\rho\sigma } p^\lambda \gamma^{\mu\ } } & = & \gamma^\mu g^{\{\rho\sigma } p^{\lambda\ } } + ( \mu \leftrightarrow \rho ) + ( \mu \leftrightarrow \sigma ) + ( \mu \leftrightarrow \lambda ) , \nonumber\\ g^{\{\rho\sigma } p^\lambda p^\tau \gamma^{\mu\ } } & = & \gamma^\mu g^{\{\rho\sigma } p^\lambda p^{\tau\ } } + ( \mu \leftrightarrow \rho ) + ( \mu \leftrightarrow \sigma ) + ( \mu \leftrightarrow \lambda ) + ( \mu \leftrightarrow \tau ) , \nonumber\\ \gamma^{\{\mu } p^\rho p^{\sigma\ } } & = & \gamma^\mu p^\rho p^\sigma + \gamma^\rho p^\mu p^\sigma + \gamma^\sigma p^\mu p^\rho , \nonumber\\ \gamma^{\{\mu } p^\rho p^\sigma p^{\lambda\ } } & = & \gamma^\mu p^\rho p^\sigma p^\lambda + ( \mu \leftrightarrow \rho ) + ( \mu \leftrightarrow \sigma ) + ( \mu \leftrightarrow \lambda ) , \nonumber\\ \gamma^{\{\mu } p^\rho p^\sigma
p^\lambda p^{\tau\ } } & = & \gamma^\mu p^\rho p^\sigma p^\lambda p^\tau + ( \mu \leftrightarrow \rho ) + ( \mu \leftrightarrow \sigma ) + ( \mu \leftrightarrow \lambda ) + ( \mu \leftrightarrow \tau ) , \nonumber\\ g^{\mu\{\lambda } p^\tau \gamma^{\rho\ } } & = & g^{\mu\lambda } p^\tau \gamma^\rho + g^{\mu\tau } p^\lambda \gamma^\rho + g^{\mu\rho } p^\tau \gamma^\lambda + g^{\mu\tau } p^\rho \gamma^\lambda + g^{\mu\rho } p^\lambda \gamma^\tau + g^{\mu\lambda } p^\rho \gamma^\tau . \nonumber\end{aligned}\ ] ]    in those expressions , the mass terms are kept explicitly , which are appropriate not only for the light - quark case but also for the heavy - quark one .
the vacuum matrix elements in the @xmath14-dimension space are @xmath171 , \nonumber\\ & & \textrm{tr}_c \left<0\left| g_{\mu\nu;\lambda } g_{\rho\sigma;\tau } \right|0\right > \nonumber\\ & & \quad\quad = -\frac{8}{27d^2(d-1 ) } g_s^2 \sum_{u , d , s } \left
< g_s\bar{\psi}\psi\right>^2 \nonumber\\ & & \quad\quad \times \left [ 2 g_{\lambda\tau } ( g_{\mu\rho}g_{\nu\sigma } - g_{\mu\sigma}g_{\nu\rho } ) + g_{\lambda\rho } ( g_{\tau\mu}g_{\sigma\nu } - g_{\tau\nu}g_{\sigma\mu } ) \right .
\nonumber\\ & & \quad\quad \left .
- g_{\lambda\sigma } ( g_{\tau\mu}g_{\rho\nu } - g_{\tau\nu}g_{\rho\mu } ) \right ] , \nonumber\\ & & \left<0\left| g^a_{\mu\nu } g^b_{\rho\sigma } g^c_{\lambda\tau } \right|0\right > \nonumber\\ & & \quad\quad = \frac{\left < fg^3\right>}{24 d ( d-1 ) ( d-2 ) } f^{abc } \nonumber\\ & & \quad\quad \times \left\ { \left [ g_{\mu\rho } ( g_{\nu\lambda}g_{\sigma\tau } - g_{\nu\tau}g_{\sigma\lambda } ) - g_{\mu\sigma } ( g_{\nu\lambda}g_{\rho\tau } - g_{\nu\tau}g_{\rho\lambda } ) \right ] \right .
\nonumber\\ & & \quad\quad \left .
- \left [ g_{\nu\rho } ( g_{\mu\lambda}g_{\sigma\tau } - g_{\mu\tau}g_{\sigma\lambda } ) - g_{\nu\sigma } ( g_{\mu\lambda}g_{\rho\tau } - g_{\mu\tau}g_{\rho\lambda } ) \right ] \right\ } , \nonumber\\ & & \textrm{tr}_c \left<0\left| g_{\mu\nu } g_{\rho\sigma } g_{\lambda\tau } \right|0\right > \nonumber\\ & & \quad\quad = \frac{i}{4 d ( d-1 ) ( d-2 ) } \left <
g_s^3fg^3\right > \nonumber\\ & & \quad\quad \times \left\ { \left [ g_{\mu\rho } ( g_{\nu\lambda}g_{\sigma\tau } - g_{\nu\tau}g_{\sigma\lambda } ) - g_{\mu\sigma } ( g_{\nu\lambda}g_{\rho\tau } - g_{\nu\tau}g_{\rho\lambda } ) \right ] \right .
\nonumber\\ & & \quad\quad \left .
- \left [ g_{\nu\rho } ( g_{\mu\lambda}g_{\sigma\tau } - g_{\mu\tau}g_{\sigma\lambda } ) - g_{\nu\sigma } ( g_{\mu\lambda}g_{\rho\tau } - g_{\mu\tau}g_{\rho\lambda } ) \right ] \right\ } , \nonumber\\ & & \left<0\left| g^a_{\mu\nu } g^b_{\rho\sigma;\lambda\tau } \right|0\right > \nonumber\\ & & \quad\quad = \delta^{ab } \left\ { \left ( a \times \sum_{u , d , s } \left
< g_s\bar{\psi}\psi\right>^2 - b \times \left < g_s f g^3 \right > \right ) \right .
\nonumber\\ & & \quad\quad \times \left [ 2 g_{\lambda\tau } ( g_{\mu\sigma}g_{\nu\rho } - g_{\mu\rho}g_{\nu\sigma } ) + g_{\rho\tau } ( g_{\mu\sigma}g_{\nu\lambda } - g_{\mu\lambda}g_{\nu\sigma } ) \right .
\nonumber\\ & & \quad\quad \left .
- g_{\sigma\tau } ( g_{\mu\lambda}g_{\nu\rho } - g_{\mu\rho}g_{\nu\lambda } ) \right ] \nonumber\\ & & \quad\quad + \left ( a \times \sum_{u , d , s } \left < g_s\bar{\psi}\psi\right>^2 + b \times \left < g_s f g^3
\right > \right ) \nonumber\\ & & \quad\quad \left .
\times \left [ g_{\mu\tau } ( g_{\rho\nu}g_{\sigma\lambda } - g_{\rho\lambda}g_{\nu\sigma } ) + g_{\nu\tau } ( g_{\rho\lambda}g_{\sigma\mu } - g_{\rho\mu}g_{\sigma\lambda } ) \right ] \right\ } , \nonumber\\ & & \textrm{tr}_c \left<0\left| g_{\mu\nu } g_{\rho\sigma;\lambda\tau } \right|0\right > \nonumber\\ & & \quad\quad = 4\left ( a \times g_s^2 \sum_{u , d , s } \left
< g_s\bar{\psi}\psi\right>^2 - b \times \left < g_s^3 f g^3
\right > \right ) \nonumber\\ & & \quad\quad \times \left [ 2 g_{\lambda\tau } ( g_{\mu\sigma}g_{\nu\rho } - g_{\mu\rho}g_{\nu\sigma } ) + g_{\rho\tau } ( g_{\mu\sigma}g_{\nu\lambda } - g_{\mu\lambda}g_{\nu\sigma } ) \right .
\nonumber\\ & & \quad\quad \left .
- g_{\sigma\tau } ( g_{\mu\lambda}g_{\nu\rho } - g_{\mu\rho}g_{\nu\lambda } ) \right ] \nonumber\\ & & \quad\quad + 4\left ( a \times g_s^2 \sum_{u , d , s } \left < g_s\bar{\psi}\psi\right>^2 + b \times \left < g_s^3 f g^3
\right > \right ) \nonumber\\ & & \quad\quad \times \left [ g_{\mu\tau } ( g_{\rho\nu}g_{\sigma\lambda } - g_{\rho\lambda}g_{\nu\sigma } ) + g_{\nu\tau } ( g_{\rho\lambda}g_{\sigma\mu } - g_{\rho\mu}g_{\sigma\lambda } ) \right ] , \nonumber\end{aligned}\ ] ] where @xmath172 and the symbol @xmath173 stands for tracing to the colour matrixes . obviously , @xmath174    our results agree with those of refs.@xcite , except for the matrix element @xmath175 .
the result for @xmath175 given ref.@xcite is incorrect .
it is noted that the correct equation for @xmath176 is @xmath177   \nonumber\\ & & + g^a_{\mu\alpha;\alpha\nu}- g^a_{\nu\alpha;\alpha\mu},\end{aligned}\ ] ] while in ref.@xcite , the coefficient of the first term has been taken as @xmath178 .
t. huang , t. zhong and x. g. wu , phys .
d * 88 * , 034013 ( 2013 ) ; t. huang , x. g. wu and t. zhong , chin .
. lett . * 30 * , 041201 ( 2013 ) ; x.g .
wu , t. huang and t. zhong , chin .
c * 37 * , 063105 ( 2013 ) ; x. g. wu and t. huang , phys .
d * 82 * , 034024 ( 2010 ) ; x. g. wu and t. huang , phys .
d * 84 * , 074011 ( 2011 ) .
x. d. xiang , x. n. wang and t. huang , commun .
. phys . * 6 * , 117 ( 1986 ) ; s. v. mikhailov and a. v. radyushkin , phys .
d * 45 * , 1754 ( 1992 ) ; a. p. bakulev , k. passek - kumericki , w. schroers and n. g. stefanis , phys .
d * 70 * , 033014 ( 2004 ) .
s. friot , d. greynat and e. de rafael , jhep * 10 * , 043 ( 2004 ) ; y. chung _ et al_. , phys .
b * * 102 * * , 175 ( 1981 ) ; nucl . phys .
b * 197 * , 55 ( 1982 ) ; h. g. dosch , m. jamin and s. narison , phys .
b * 220 * , 251 ( 1989 ) .
s. j. brodsky , t. huang , and g. p. lepage , in _ particles and fields-2 _ , proceedings of the banff summer institute , ban8 ; alberta , 1981 , edited by a. z. capri and a. n. kamal ( plenum , new york , 1983 ) , p. 143 ; g. p. lepage , s. j. brodsky , t. huang , and p. b.mackenize , _
ibid_. , p. 83 ; t. huang , in _ proceedings of xxth international conference on high energy physics _ , madison , wisconsin , 1980 , edited by l. durand and l. g pondrom , aip conf
69 ( aip , new york , 1981),p . 1000 .


we report the detection of off - pulse emission from two long period pulsars psr b0525 + 21 and psr b204516 using the interferometric mode of the giant metrewave radio telescope ( gmrt ) in the 325 mhz frequency band .
the `` main - pulse '' ( or the equivalent term `` on - pulse '' used in this paper ) is the emission from within the polar - cap region of a pulsar ; and the off - pulse is the emission outside this main - pulse .
off - pulse emission from pulsars has been a subject of interest since the discovery of pulsars four decades ago .    according to the rotating vector model ( rvm ) proposed by radhakrishnan & cooke ( 1969 )
the observed pulses are due to relativistically beamed radiation along the open dipolar field lines .
the plane of linear polarization traces the magnetic field line associated with the emission at every instant .
several statistical studies have revealed that pulsar radio emission is in the form of a circular beam arising from the polar cap ( i.e. the main - pulse ) with an opening angle @xmath2 , and from a height above the magnetic pole equal to 12 per cent of the radius of the light cylinder .
the width of the observed pulse is a geometrical function of @xmath2 ( which itself is a function of the emission height ) , the pulsar period , the angle between the rotation axis and the magnetic axis ( @xmath3 ) and the angle between the magnetic axis and the line of sight to the observer ( @xmath4 ) .
the angles @xmath3 and @xmath4 can be estimated by fitting the rvm to the swing of the polarisation position angle ( ppa ) across the pulse @xcite .
the observed main - pulse covers 5 - 10 per cent of the pulse period for 90 per cent of the pulsars . in several pulsars
this main - pulse consists of multiple components including low - level bridge emission between such components ( rathnasree & rankin 1995 ) . the pulse width scales as 1/sin(@xmath3 ) , from the geometry of the emission model ; in the rare case of a pulsar with closely aligned magnetic and rotation axes ( i.e. @xmath5 or 180 ) the pulse width can be as high as 100 per cent .
several pulsars have low level but detectable bridge emission between two widely spaced pulse components , making them candidates for an aligned rotator geometry .
in rare cases when @xmath3 is close to 90 , emission from the opposite pole ( i.e. the inter - pulse ) may be observed at 180phase from the main - pulse . in all these instances the observed emission is believed to lie within the polar cap ( i.e. main - pulse ) as described earlier .    in a few pulsars , low level emission components known as pre-/post - cursors ( ppc )
have been observed outside the main - pulse ( e.g. see mitra & rankin 2010 for a discussion ) . in these pulsars ,
all geometrical evidence indicate that the main pulse is consistent with emission from open field lines .
the ppc components appear highly polarized and are far from the main - pulse .
the discovery of the ppc component about 60away from the main pulse in psr b0943 + 10 is particularly interesting , as the line of sight almost grazes the emission cone for this pulsar @xcite .
hence the ppc emission originates either from a much larger height , where due to spreading of dipolar field lines the ppc component can lie far away from the main - pulse , or from the regions of closed field lines .
some pulsars occasionally emit giant pulses which are believed to arise close to the light cylinder rather than the polar cap .
interestingly ppc pulsars have low spin - down luminosities ( @xmath6  erg  s@xmath1 ) , whereas the giant pulses arise from pulsars with @xmath7  erg  s@xmath1 .
the ppc components and giant pulses are potential sources of magnetospheric off - pulse emission in pulsars .
they challenge the conventional wisdom of pulsar radio emission arising only from open magnetic field lines close to the pulsar polar cap and raise questions about the origin of these emission components .    the other possible source of off - pulse emission is the interstellar medium ( ism ) around the pulsar .
the pulsar loses most of its rotational energy in the form of a relativistic wind which , when confined by the surroundings , may form a pulsar wind nebula ( pwn ) .
several types of pwne are formed depending on the confinement mechanism .
young and energetic pulsars are often located in their associated supernova remnant ( snr ) .
the pulsar wind streaming into the ambient medium produces standing shocks resulting in a plerionic pwn ( like the one observed in the crab nebula ) . in older pulsars , where the surrounding snr is likely to be dissipated
, the relativistic particles may interact with the ism magnetic field and radiate , creating ghost remnants ( proposed by blandford et al .
1973 , but not observed till date ) .
finally , a pulsar moving through the ism with supersonic speed can produce a bow - shock nebula , where ram pressure balance is established between the pulsar wind and the ambient medium .
several such bow shock nebulae have been detected in h@xmath3 , and some in the radio as plerionic bow shocks , but never in both ( see chatterjee & cordes 2002 for a detailed study ) .
several searches for pwne have been conducted in the past with varying sensitivities and resolutions .
however only in about 10 young pulsars have radio emission been detected outside the main - pulse and all of these are believed to be associated with pwne . in all such cases
the associated pulsars are young ( 10@xmath8 - 10@xmath9yr ) , have high spin - down luminosities ( @xmath10  erg  s@xmath1 ) and all except one are associated with snrs @xcite .
it is believed that the pulsar wind and the environment change as pulsars slow down and age , making them less likely to harbour pwne .
the present study was designed to detect radio emission away from the polar caps and not associated with pwne .
therefore we targeted old and less energetic pulsars unlikely to harbour pwne .
their profiles did not show any components outside the main - pulse .
the targets were neither aligned rotators nor were inter - pulsars .
we attempted to detect emission within a pulse phase range of approximately 80 - 250  from the peak in the main - pulse .
we describe the gated interferometric observations and the data analysis in section 2 , the results in section 3 , the additional tests we carried out to confirm that the off - pulse emission were not instrumental artifacts in section 4 , and discuss the implications of the detections in section 5 .
we chose two pulsars , b0525 + 21 and b204516 , whose pulse profiles did not show any features outside the main - pulse @xcite and were stronger than 50 mjy to be able to put a non - detection upper limit of at most 5 per cent of the pulsed flux .
the properties of the targets are listed in table [ pulsprop ] .
we targeted long period pulsars because : + 1 .
pwne are so far only known from short - period , energetic pulsars in the vicinity of snr .
a long period selects against both high energy pulsars and snr association , and hence pwne .
we wanted at least 8 bins ( of 131 or 262 ms each ) across the pulse period to cleanly separate the off - pulse and the on - pulse regions .
the temporal broadening across the 16 mhz bandwidth due to dispersion by the ism was 183.3 ms for psr b0525 + 21 ( binwidth 262 ms , gate width 1.31 s ) , and 47.8 ms for psr b204516 ( binwidth 131 ms , gate width 0.655 s ) .
therefore , we did not have to dedisperse the signal .
lccccccc pulsar&period & dm & dist&@xmath11 & @xmath12 & @xmath13&@xmath14 + & s & @xmath15 & kpc & year & erg  @xmath16 & km  @xmath16 & mjy + b0525 + 21 & 3.7455 & 50.937 & 2.28 & 1.48@xmath17 & 3.0@xmath18 & 229 & 80.5 + b204516 & 1.9616 & 11.456 & 0.95 & 2.84@xmath17 & 5.7@xmath18 & 511 & 169.3 +      interferometric observations are better than single dish measurements for this study for several reasons : + 1 .
an imaging interferometer is only insensitive to constant flux density background along the spatial axes while a standard pulsar receiver is usually insensitive to the constant background along the time axis ; the detection of off - pulse emission is essentially an attempt to find such a constant background along the time axis .
self - calibration of interferometric data can correct instrumental and atmospheric gain fluctuations on very short time scales .
the corrections are determined by the flux densities of the constant and bright background sources in the field and hence would not be affected by the pulse variation of the relatively weak pulsar .
the higher spatial resolution of an interferometer greatly reduces the coincidence of unrelated sources ; thereby reducing the probability that the off - pulse emission is from an unrelated source within the synthesized beam .
we imaged the pulsars with the giant metrewave radio telescope ( gmrt ) , an aperture - synthesis radio interferometer located near pune , india @xcite .
the 30 antennas of 45 m diameter provide a maximum baseline of 27 km and can be operated at 6 frequency bands between 50 and 1450 mhz .
we observed the pulsars at 325 mhz with a 16 mhz bandwidth split into 128 channels .
the frequency was chosen for its optimal combination of resolution ( 10 ) and sensitivity ( few 100 @xmath19jy in 4 hours of observing ) .
the high frequency - resolution was useful in flagging narrow - band radio - frequency interference ( rfi ) and to avoid bandwidth smearing .
the shortest integration output by the hardware correlator is 0.131 s , though this data is usually averaged and median filtered to provide a standard data output integration of 16.77s .
the two pulsars were first observed on 19 may , 2009 with 0.262s integration to keep the data output rate and size manageable . this time resolution ( of 7 bins )
was found to be inadequate in the case of psr b204516 while the data was of poor quality ( @xmath20 missing antennas ) for both targets .
we reobserved the sources on 19 january , 2010 , with integrations of 0.131s for psr b204516 ( period 1.9616s ) , and 0.262s for psr b0525 + 21 ( period 3.7455s ) .
this allowed us to divide both the pulsar periods into 14 time bins .
several earlier studies used the technique of online gating to separate the on- and off - pulse regions using the timing information of pulsars , and accumulated time - averaged on- and off - pulse data @xcite .
this resulted in relatively smaller size of data files and hence shorter computation times which was a major consideration 1 - 2 decades ago .
the disadvantage was the impossibility of rectifying any errors during the gating process ( see results and discussions in gaensler et al .
we recorded the data at sufficiently high resolution and applied offline gating .      in the absence of accurate absolute time reference in the interferometric mode we used offline gating to isolate the off - pulse region of the pulsar period .
the 30 antennas of gmrt provide 30 self , and 435 cross visibilities at each integration interval ( 0.131s or 0.262s ) .
we folded the self - data to obtain the pulsar profile and identify the on - and off - pulse regions .
this was used to gate the cross - antenna visibilities into on- and off - pulse data sets .
the temporal and frequency gain variations in the self - data were corrected by normalising the instantaneous values by the local median .
the timescale for the temporal median was 30 - 50 times the pulsar period to ensure that the pulse variation was retained .
the data from all frequency channels and antennas were averaged with robust sigma - clipping and processed using the standard folding algorithm for pulsars ( hankins & rickett 1975 ) to determine the pulse profile .
the pulsar profiles and the off- and on - pulse gated are shown in figure [ fig_profilep ] .
the off- and on - pulse gates of psr b0525 + 21 were centred on phases 0.24 and 0.64 , respectively .
the 5 bins nearest to each were averaged to construct the off- and on - pulse data sets .
the corresponding phases for psr b204516 were 0.21 and 0.67 , respectively .
the folding and gating procedures were carried out using software developed for this purpose .
the on- and off - pulse data were cleaned of radio frequency interference ( rfi ) using the rfix algorithm @xcite .
they were separately calibrated , flagged and imaged in a standard manner using the nrao aips package .
the flux - density scale of the images was determined from observations of the calibrator 3c48 and the measurements of baars et al.(1977 ) with the latest corrections of perley et al ( 1999 ; in aips ) .
the data sets were taken through several rounds of phase self - calibration followed by a final round of amplitude and phase self - calibration .
* psr b0525 + 21 : * observations of the calibrator 0521 + 166 were interspersed with the target to correct for the amplitude and phase gain fluctuations .
the presence of the extremely strong and extended crab nebula ( flux @xmath21 1000 jy and angular size @xmath22 10 ) at the edge of the primary beam ( 1.5away ) resulted in enhanced noise and strong ripples in the initial image .
therefore , subsequent analyses were carried out with a lower uv cut - off at 1.5 k@xmath23 which reduced the artifacts and noise .
this eliminated all structures larger than 2.3  from the image ; this did not affect our purpose as our target source was expected to be much smaller than 2 .
* psr b204516 : * the calibrator 2137207 was observed for amplitude and phase calibration .
a pointing error resulted in the pulsar being located 35  away from the field centre ( 60 per cent gain level of the primary beam ) .
the primary beam correction was applied to the image to obtain the correct flux density for the pulsar .    ]
10 mjy ) in the on- and off - pulse images .
the scatter around unity confirmed that the flux - scale was the same for both .
( b0525 + 21 : top ; b204516 : bottom ) [ fig_onvoff],title="fig : " ]   10 mjy ) in the on- and off - pulse images .
the scatter around unity confirmed that the flux - scale was the same for both .
( b0525 + 21 : top ; b204516 : bottom ) [ fig_onvoff],title="fig : " ]
lccccccccc pulsar & frequency & time & beam & size & on rms & on flux & off rms & off flux & int .
flux + & ( mhz ) & ( min ) & ( ) & ( pc ) & ( mjy@xmath24)&(mjy@xmath24 ) & ( mjy@xmath24)&(mjy@xmath24)&(mjy@xmath24 ) + b0525 + 21 & 332.9 & 160 & 9.5@xmath256.5 & 0.088 & 0.55 & 80.2@xmath265.8 & 0.45 & 3.9@xmath260.5 & 30.0@xmath262.1 + b204516 & 317.1 & 180 & 11.9@xmath257.2 & 0.042 & 2.10 & 305@xmath2622 & 0.65 & 4.3@xmath261.1 & 110.5@xmath267.9 +    we detected off - pulse emission in both the pulsars , coincident with the location in the on - pulse image . the results are summarized in table [ results ] .
the on- and the off - pulse emissions appeared to be unresolved for both the pulsars .
we compared the flux densities of sources in the field of view to check the flux - density calibration .
the ratios of on- and off - image flux densities of sources stronger than 10 mjy are plotted in figure [ fig_onvoff ] .
there were 43 and 93 comparison sources in the fields of b0525 + 21 and b204516 , respectively .
the scatter of values around unity , especially for strong sources ( which have smaller fractional flux errors ) confirmed the similarity of the flux scale of the on- and off - images for both the pulsars .
* psr b0525 + 21 : * the pulsar was detected in the off - pulse image at the 8.6@xmath27 level ( flux - density 3.9@xmath260.5 mjy , rms noise 0.45 mjy / beam ) .
the flux density of the pulsar in the on - pulse image was 80.2@xmath265.8 mjy ( rms noise 0.55 mjy / beam ) .
the resolution element at the distance of the pulsar corresponds to 0.09 pc . the data from 19 may , 2009 , yielded an off - pulse detection of 3.6 mjy .
* psr b204516 : * the pulsar was detected in the off - pulse image at the 6.6@xmath27 level ( flux density 4.3@xmath261.1 mjy ; rms noise 0.65 mjy / beam ) .
the flux density of the pulsar in the on - pulse image was 305.2@xmath2621.9 mjy ( rms noise 2.1 mjy / beam ) .
the resolution element at the distance of the pulsar corresponds to 0.04 pc .
the flux - densities averaged over the entire pulsar period is 30.0@xmath262.1 for b0525 + 21 and 110.5@xmath267.9 for b204516 .
these values are smaller than expected ( table [ pulsprop ] ) by factors of 1.5 and 2.7 , which can be caused by effects such as refractive interstellar scintillations @xcite .
we confirmed the correctness of our flux scale by comparing the flux densities of the phase calibrators to their known values ( less than 1 per cent ) and by comparing the flux densities of 3 - 5 strong sources in each field to the interpolated values from nvss ( 1.4 ghz ) and vlss ( 74 mhz ) .
the quoted errors on the flux densities were obtained by adding in quadrature the image rms , calibration errors and the error on the value of the primary flux - density calibrator @xcite .
the signal - to - noise ratio in the published profiles of the pulsars @xcite are insufficient to conclude whether the off - pulse emission detected here is a pedestal ( throughout the period ) or a pulse confined to the off - pulse gate .
firstly , we detected the off - pulse emission in two data sets of b0525 + 21 observed 8 months apart , and within the noise the two measurements were identical .
we can safely rule out the detections being chance occurance of noise peaks .
secondly , b204516 was positioned almost at the half - power point of the primary beam while b0525 + 21 was close to the field centre .
the detection in both cases suggests that they are not the spurious structures occasionally seen at the centre of the field .
we discuss below several effects which could result in spurious sources in the off - pulse data at the location of the pulsar ; and demonstrate that they are unlikely to be responsible for the same .      to determine the probability of finding an unrelated source coinciding with the pulsar we used the vla first catalogue @xcite to determine the following relationship between source counts ( n ) and flux - density ( s ) : @xmath28 where @xmath29 is the flux at 1.4 ghz in mjy , @xmath30 is the flux bin width at 1.4 ghz in mjy and @xmath31 is the area of sky under consideration in @xmath32 .    the area of the synthesized beam is @xmath22100 @xmath33 at 325 mhz .
the flux - density for consideration in the equation is between 3@xmath25image - rms and 5 - 10 per cent of the pulsar flux , translated to 1.4 ghz ( using spectral index @xmath34 ) .
the upper value is related to a realistic upper limit to the flux fraction in the off - pulse from the single - dish pulse profile @xcite .
this yielded a probability of coincident unrelated sources of @xmath35 .
the probability of detecting coincident sources for both pulsars is the square of the above number and therefore highly improbable .
this scenario requires that approximately 1 - 4 per cent of the on - pulse data fall within the off - pulse gate .
however , any monotonic drift in the clock , without recovery , would have smeared out the pulse obtained by folding the self - data ; on the contrary we see the on - pulse detected with a signal- to - noise ratio of 50 - 100 ( figure  [ fig_profilep ] ) .
alternatively , the clock could have maintained an accurate long term average but with large excursions in the values of the individual time stamps , i.e. excursions of the order of 1 - 2 s ( half the pulsar period ) for an output data rate of one visibility per 0.131s or 0.262s .
this sounds somewhat contrived ; if this were the case we should have expected a histogram of the time difference between adjacent data to show a large scatter .
the histograms plotted in figure  [ fig_timehist ] show a tight scatter ( rms 0.007s ) around the expected mean values ( 0.262s and 0.131s ) .    ,
title="fig : " ] , title="fig : " ]    additionally , any artifact due to time - smearing should have resulted in the shorter - period and faster - sampled b204516 having a higher off - pulse fraction , whereas the ratio is much higher for b0525 + 21 .
this scenario requires temporal correlation between data separated by at least half the pulsar period .
the basic sampling interval of the gmrt correlator is 31.25 nanoseconds .
a cross - spectrum is output every 16@xmath19s by fourier transforming 512 such samples .
the data output every 16@xmath19s should be independent from all other such data in the time series .
8188 ( or 16376 ) such independent samples are averaged to obtain the short - term acquisition ( sta ) date output at 0.131 ( or 0.262 s ) .
given the independence at the 16@xmath19s level it is difficult to envisage correlation between one sta and the next , even more so across 6 - 10 stas ( half the pulsar period ) . nevertheless this was a possibility .
observers using the standard interferometric mode of the gmrt would not be affected as none of their targets vary on sub - second timescales , and their integration time is rarely less than 2s ( more typically 17s ) .
the leakage would only redistribute the flux density along the time axis without causing any change in the measured flux density .
observers using the pulsar mode would not have detected such an effect because they are only sensitive to the relative flux in excess of the base level .
we tested this possibility by estimating the temporal correlation in the 0.131s data series . assuming that the off - pulse was entirely due to leaking of the signal from the on - pulse bin to the off - pulse bins we estimated the required leakage for a spurious detection of the observed level .
there was no way of introducing a narrow , clean pulse into the gmrt receiver system to directly measure such a leakage .
so , we investigated the consequence of such a correlation on the noise data .
first , we ran simulations to estimate the auto - correlation function for noise in the presence of leakage sufficient to generate the observed off - pulse , i.e. we generated a noise series and smeared each data point into subsequent data points according to a particular time profile  e.g. constant leakage into the next 13 bins ; or linearly decreasing leakage into the next 13 bins .
the value 13 is related to the period of the two pulsars which covered 14 time bins .
we then recorded 6@xmath2510 minute scans of noise data ( 0.131s integration ) but with front - end termination which sealed the telescope at the antenna feeds .
this ensured that no temporally continuous external source ( cosmic sources and rfi ) , which would be correlated across the entire observing session , was present in the noise data .
the auto - correlation for each baseline - channel ( 435@xmath25120 ) data was separately calculated for each scan . the median autocorrelation profile and the scatter
are shown in figure  [ autocorrfull ] .
the worst data are at the 0.1 per cent level but the median is only about 0.04 per cent .
the profile does not fall to zero at large time lags , indicating non - stationarity of the noise signal .
we obtained a similar behaviour by introducing a time - varying mean level into the noise in our simulations .
this time variation in the mean level may reflect the system gain variations as well , which would have been corrected to a greater or lesser extent by self - calibration during imaging .
therefore the profile is a firm upper limit to the contribution of temporal leaking of signals to an off - pulse detection .
it should be noted that interferometric imaging is an excellent filter of bad data ; a compact source ( like the off - pulse detection ) would require most of the baselines to show high temporal auto - correlation ; a few bad baselines with high auto correlation would only result in noisy ripples across the image and not a localised source .
therefore the median line is a more accurate measure of the leakage for our study than the extreme points in the scatter .
we note that non - stationarity and gain fluctuation only occur at the 0.05 per cent level in the gmrt , and would not have a discernible effect for most imaging exercises .    ]
figure  [ autocorrmag ] shows a magnified representation of the relevant part of the noise auto - correlation plot .
a1 and a2 ( both psr b0525 + 21 ) and b1 ( psr b204516 ) indicate the required noise auto - correlation if the detected off - pulse were due to the leaking of the signal from the on - pulse bin . the analysis recipe for psr b0525 + 21 ( with similar arguments for psr b204516 )
is as follows : + 1 .
the pulsar period includes 14 time bins , each of 0.262 s. the on- and off - pulse fluxes measured over 5 time bins are 80.2 and 3.9 mjy , respectively .
higher time resolution observations @xcite show that the main - pulse narrower than one of our bins .
+ 2 . a leakage of 0.9 per cent of the pulse flux from the first bin into each of the next 13 bins will explain the observed off - pulse .
this corresponds to an auto correlation level of 0.0087 .
the auto - correlation profile calculated for gmrt noise ( the lowest curve n in figure  [ autocorrmag ] ) , is smaller than the required amount by a factor of at least 20 .
a leakage reducing with temporal distance is more likely .
a linearly reducing leakage ( curve a2 ) requires a similar average auto - correlation in the off - pulse region and a higher correlation at short lags .
we conclude that temporal leaking of the signal is incapable of explaining the detected level of off - pulse emission .
the curve b2 represents the corresponding model for psr b204516 : off - pulse 4.3 mjy ; on - pulse flux 305.2 mjy ; required constant leakage 0.43 per cent ; required correlation in the off - pulse region 0.003 ; discrepancy is a factor of @xmath22 6 .    additionally , the argument against time - smearing ( sec .
4.2 , second para ) is also valid against leaking of the signal along the time series .    .
the plots a1 , a2 ( psr b0525 + 21 ) and b1 ( psr b204516 ) are the required ( simulated ) noise auto - correlation for a spurious off - pulse detection due to temporal leaking of the pulse into subsequent bins .
the thick black curve is the observed noise auto - correlation profile ( from figure   [ autocorrfull ] ) .
[ autocorrmag ] ]
a literature survey yielded many reports of searches for `` continuous emission in the direction of pulsars '' ( table [ litrsurvey ] ) . here
we shall discuss only those in which gating techniques were applied to target off - pulse emission , which could have been either from the pulsar itself or from pwne .
perry & lyne ( 1985 ) reported off - pulse emission from 4 pulsars ( psrs b1541 + 09 , b160400 , b1929 + 10 and b2016 + 28 ) using a gated two - element interferometer at 408 mhz .
subsequent studies showed that two of these pulsars , b1541 + 09 and b1229 + 10 , were aligned rotators , undermining the claim of off - pulse origin @xcite .
the off - pulse emission in the other two pulsars , b1604 - 00 and b2016 + 28 , were later identified with unrelated background sources ( storm & van someren greve 1990 ; hankins et al . 1993 ) .
bartel et al .
( 1985 ) reported a non - detection from psr b0329 + 54 using gated vlbi interferometry at 2.3 ghz .
strom & van someren greve ( 1990 ) reported non - detections in 43 pulsars using the wsrt in gated mode at 327 and 609 mhz .
stappers et al .
( 1999 ) reported non - detection in 4 pulsars using atca in gated interferometric mode at 1.3 and 2.2 ghz .
gaensler et al .
( 2000 ) searched for unpulsed emission from 27 energetic and/or high velocity pulsars using gating interferometry with the vla and atca at 1.4 ghz and found emission in two cases .
however they could not conclusively determine whether these detections corresponded to pwne or the pulsar magnetosphere ; they even suggested that the detections could be spurious and a result of improper online gating .
these detections were at the level of 2@xmath27 ( psr b1634 - 45 ) and 13@xmath27 ( psr b1706 - 16 ) .
there have been detections of several radio pulsars with associated plerionic bow shock nebula : b1951 + 32 @xcite , b175724 @xcite , b1853 + 01 @xcite and b164343 @xcite .
all of these are highly energetic pulsars associated with snr @xcite .
gaensler et al . (
1998 ) reported the detection of off - pulse emission from psr b090649 , which appears to be a pwn generated by a slow moving pulsar in the dense ism .
llll non detections/ & detections & tentative & pwne detection + subsequent snr & later refuted & detection & + 1 .
schnhardt 1973 , 1974 & 1 .
gopal - krishna 1978 & 1 .
gaensler et al . 2000
( @xmath36)&1 .
frail & kulkarni 1991 + 2 .
weiler et al .
glushak et al .
1981 & & 2 .
frail et al .
cohen et al .
perry & lyne 1985 ( @xmath36 ) & & 3 .
gaensler et al .
bartel et al .
1985 ( @xmath36 ) & & & 4 .
giacani et al .
strom & van someren & & & + greve 1990 ( @xmath36 ) & & & + 6 .
hankins et al .
1993 ( @xmath36 ) & & & + 7 .
frail & & & & + scharringhausen 1997 & & & + 8
. stappers et al .
1999 ( @xmath36 ) & & & +   +      we discuss here the unsuccessful search for off - pulse emission in psr b0525 + 21 by weiler et al .
( 1974 ) and perry & lyne ( 1985 ) .
weiler et al . ( 1974 ) were looking for extended emission around the pulsar without employing any kind of gating .
they had a synthesized beam size of 47@xmath25 133 with a noise rms of 2.5 mjy / b .
it is apparent that their setup was inadequate for the present level of detection .
perry & lyne ( 1985 ) employed a two - element interferometer with hardware gating at 408 mhz for their studies but could neither image the fields and nor correct for sensitivity variations .
they listed considerable emission in the off - pulse bin for b0525 + 21 ( in fact greater than for b1929 + 10 and b2016 + 28 , which were claimed to be positive detections ) but the noise was also very high precluding a positive detection .
the major advantages of our effort were : + 1 .
low frequency of observation : the steep spectrum pulsar is likely to dominate over emission from any associated nebula with flatter spectrum ; + 2 . a multi - element interferometer : the gmrt is currently the most sensitive low - frequency instrument .
the ability to correct gain variations , filter bad data using imaging residuals and rfi excision tools @xcite allowed us to make sufficiently deep images .
+ 3 . off - line gating of data
: this allowed secure separation of the off- and on - pulse sections .
the confinement of the relativistic wind from pulsars generates pwne which are luminous across the electromagnetic spectrum in synchrotron , inverse compton , and optical line emission from the shocked regions .
all known pwne seen around radio pulsars have spin - down luminosities @xmath7  erg  s@xmath1 .
they all appear to be young and , with the exception of psr b090649 @xcite , are associated with snr .
two classes of pwne can be formed in the absence of associated snr : one , where the pulsar wind is confined by the density of the ism  known as static pwn ; and the other , where the wind is confined by the ram - pressure of motion of the pulsar through the ism  called bow - shock nebula .
our two target pulsars , b0525 + 21 and b204516 , are not associated with snrs .
so , we explored the possibility that the detected off - pulse emission was due to a pwn .
the efficiency factor @xmath37 is defined as the ratio of the radio bolometric luminosity ( @xmath38 ) of a pwn to the spin - down energy ( @xmath12 ) of the pulsar @xmath39 . if we assume a typical pwn spectral index of @xmath40 0.3 its radio luminosity between @xmath41 hz and @xmath42 hz is given by @xmath43 where @xmath44 is the distance to the pwn in kpc and @xmath45 is the integrated flux of pwn at 325 mhz in mjy . using the measured flux - density ( table  [ results ] ) , source distance and the spin - down energies ( table  [ pulsprop ] ) we calculate efficiency factors of 2@xmath46 for b0525 + 21 and 2@xmath47 for b204516 .
these @xmath37 values are 1 - 2 orders of magnitude higher than those of previously known pwne ( typical @xmath48 ; see frail & scharringhausen 1997 , gaensler et al .
2000 for a discussion ) .
[ [ static - pwn ] ] static pwn + + + + + + + + + +    the relativistic particles and poynting flux emanating from the pulsar , at rest relative to the ism , will drive through the ambient medium a shock of radius @xmath49 given by @xcite @xmath50 the velocity of the shock front is given by @xmath51 here @xmath52 , where @xmath53 is the proton mass and @xmath54 is the particle density of the ism . using equation ( [ eqn_rs ] ) the required particle density for a pwn is @xmath55 where @xmath56 is the spindown power in units of @xmath57  erg  s@xmath1 , @xmath58 is the age in units of @xmath59 yr , @xmath60 is the radius of pwn in units of 0.01  pc and @xmath61 is the ism density in units of 0.01 @xmath62 . in the present exercise
we assume that the pwne in our targets are smaller than one synthesized beam width since the off - pulse emission is unresolved in both pulsars . using equation ( [ eqn_nrs ] ) and values in tables [ pulsprop ] and [ results ] we determined the ism density required to drive a pwn to be @xmath63 for b0525 + 21 and @xmath64 for b204516 .
the required particle densities are much higher than typical ism densities of @xmath65 , suggesting that these pulsars are too weak and old to power a static pwn through the ism . since our estimate of the sizes of the nebulae are upper limits the particle density estimates are lower limits .
[ [ bow - shock - pwn ] ] bow shock pwn + + + + + + + + + + + + +    this requires that the shock velocity ( @xmath66 ) be much smaller than the pulsar transverse velocity ( @xmath67 ) .
the typical ism density of @xmath68 in equation  ( [ eqn_rsv ] ) yields @xmath66 @xmath69 for b0525 + 21 and @xmath70 for b204516 , which are indeed much smaller than the @xmath67 listed in table [ pulsprop ] .
the radius of the shock is given by @xcite @xmath71    from equation ( [ eqn_rbs ] ) the particle density required to sustain a bow shock pwn is @xmath72 where @xmath56 is the spin - down power in units of @xmath73 , @xmath60 is the radius of pwn in units of 0.01  pc and @xmath74 is the velocity of the pulsar through the ism in units of 100 @xmath75 .
this requires an ism density of @xmath76 for b0525 + 21 and @xmath77 for b204516 .
these values are 3 - 4 orders of magnitude lower than the typical ism density .
this implies that the pulsars are too weak to drive a bow shock nebula of the size corresponding to the telescope resolution .
however , a bow shock interpretation can be salvaged if the nebulae are @xmath222 orders of magnitude smaller .    in summary , identifying the off - pulse emission detected here with pwne results in unrealistic values of the ism particle density .
this is not surprising because we selected as our targets middle - aged pulsars , which category has hitherto not been known to be associated with pwne .
numerous published studies provide estimates of @xmath3 and @xmath4 for the pulsars studied here .
some use an empirical / geometrical ( e / g ) approach to establish @xmath3 and @xmath4 ( lyne & manchester 1988 ; rankin 1993 ; see everett & weisberg 2001 for a discussion on the e / g approach ) .
one can also obtain @xmath3 and @xmath4 values by fitting the rvm to the ppa traverse @xcite , but these fits yield highly correlated estimates of the two angles .
this issue has been discussed in depth by everett & weisberg ( 2001 ) and mitra & li ( 2004 ) .
everett & weisberg ( 2001 ) list six different studies of psr b0525 + 21 with @xmath3 ranging from 116  to 159 values should span 0 - 180(and not 0 - 90 ) for a consistent definition of the pulsar rotation axis .
everett & weisberg provide a compilation of published @xmath3 and @xmath4 values corrected for this convention . ] .
mitra & li ( 2004 ) estimated @xmath78 for psr b204516 . both rvm
fit and the e / g approach were used to obtain the above estimates of @xmath3 .
these studies provide good evidence that our targets , b0525 + 21 and b204516 , are not aligned rotators and therefore any off - pulse emission must arise far from the magnetic pole .
the pulse profiles of both these pulsars show sharp rising and falling edges .
psr b0525 + 21 is an example of a classical double ( d ) profile , whereas psr b204516 is classified as a triple @xcite .
the main - pulse width for these pulsars span about 5 per cent of the period , and there are no obvious emission components visible outside the main - pulse . based on the pulse widths and the estimates of @xmath3 and @xmath4 , mitra & li ( 2004 ) concluded that the radio emission arises at a distance of 12 per cent of the light cylinder .
hence , the radio on - pulse emission of psr b0525 + 21 and b204516 are classical examples of emission from open dipolar field lines .
the remaining 95 per cent of the period , i.e. the off - pulse region , were hitherto thought to be radio quiet zones of the pulsar magnetosphere @xcite .
our off - pulse phase corresponds to 80 - 208(from the peak of the main - pulse ) for b0525 + 21 and 101 - 229for b204516 . in a few cases pre-/post - cursor emissions
have been detected about 60from the main - pulse , posing questions regarding their origin @xcite .
our detections are even further away from the main - pulse .
it remains to be seen if these off - pulse detections are an extreme example of ppc or if it represents unpulsed emission from throughout the pulse period .
it is to be noted that the opening angle of dipolar field lines scales as the square of the emission height .
if radio emission from pulsars originate from higher up the light cylinder ( @xmath79 1 - 2 per cent of the light cylinder ) , the polar emission could span a larger longitudinal range .
this implies that radio emission at larger distances from the neutron star surface are potential sites of ppc or off - pulse emission in pulsars .
alternatively , if the emission arises close to the neutron star surface , the emitting region will need to encroach upon the closed field lines of the neutron star .
thus , the existence of magnetospheric off - pulse emission should prove to be an important input for pulsar electrodynamic models which try to establish the relation between magnetospheric currents and coherent pulsar radio emission @xcite .
we report gmrt observations at 325 mhz which detected off - pulse emission from two long - period pulsars which have low spin - down rates and are not associated with supernova remnants .
we have adduced evidence that the signals are not artifacts of the observing procedure , nor unrelated background sources which happen to be coincident with the pulsar .
we have argued that explaining these detections as pulsar wind nebula requires ism particle densities which differ from typical measured values by several orders of magnitude .
robust estimates of the geometrical parameters of these two pulsars argue against them being aligned rotators .
this leaves the possibility of emission from the magnetosphere .
if the off - pulse emission arises from a much higher height than the on - pulse then it could arise from open field lines at the edge of the light cylinder . on the other hand
if the off - pulse emission is from the same height as the on - pulse emission then it must be associated with closed field lines .
further studies at multiple frequencies , in polarization and higher resolution are needed to establish the nature of this emission which can impose valuable constraints on the pulsar emission mechanism models .
we thank the director of ncra for discretionary observing time which made this work possible .
we also thank the staff of the gmrt for their help with these non - standard observations possible .
the gmrt is run by the national centre for radio astrophysics of the tata institute of fundamental research .
we also thank the referee , scott ransom , for useful comments .
arons , j. 1983 , , 302 , 301 athreya , r. 2009 , , 696 , 885 baars , j. w. m. , genzel , r. , pauliny - toth , i. i. k. , witzel , a. 1977 , , 61 , 99 backus , i. , mitra , d. , rankin , j. m. 2010 , , 404 , 30 bartel , n. , ratner , m. i. , shapiro , i. i. , cappallo , r. j. , rogers , a. e. e. , whitney , a. r. 1985 , , 90 , 318 becker , r. h. , white , r. l. , helfand , d. j. 1995 , , 450 , 559 blandford , r. d. , ostriker , j. p. , pacini , f. , rees , m. j. 1973 , , 23 , 145 chatterjee , s. , cordes , j. m. 2002 , , 575 , 407 cohen , n. l. , cotton , w. d. , geldzahler , b. j. , marcaide , j. m. 1983 , , 264 , 273 everett , j. e. , weisberg , j. m. 2001 ; , 553 , 341 frail , d. a. , giacani , e. b. , goss , w. m. , dubner , g. 1996 , , 464 , 165 frail , d. a. , kulkarni , s .r .
1991 , , 352 , 785 frail , d. a. , scharringhausen , b. r. 1997 , , 480 , 364 gaensler , b. m. , stappers , b. w. , frail , d. a. , johnston , s. 1998 , , 499 , 69 gaensler , b. m. , stappers , b. w. , frail , d. a. , moffett , d. a. , johnston , s. , chatterjee , s. 2000 , , 318 , 58 giacani , e. b. , frail , d. a. , goss , w. m. , vieytes , m. 2001 , , 121 , 3133 glushak , a. p. , pynzar , a. v. , udaltsov , v. a. 1981 , , 25 , 182 gopal - krishna 1978 , , 185 , 521 goldreich , p. , julian , w. h. 1969 , , 157 , 869 gould , d. m. , lyne , a. g. 1998 , , 301 , 235 hankins , t. h. , moffett , d. a. , novikov , a. , popov , m. 1993 , , 417 , 735 hankins , t. h. , rickett , b. j. 1975 , methods in comp phy . , 14 , 55 hester , j. j. , kulkarni , s. r. 1988 , , 331 , l121 hobbs , g. , lyne , a. g. , kramer , m. , martin , c. e. , jordan , c. 2004 , , 353 , 1311 lorimer , d. r. , yates , j. a. , lyne , a. g. , gould , d. m. 1995 , , 273 , 411 lyne , a. g. , manchester , r. n. 1988 , , 234 , 477 mitra , d. , deshpande , a. a. 1999 , , 346 , 906 mitra , d. , li , x. h. 2004 , , 421 , 215 mitra , d. , rankin , j. m. 2010 , submitted to perry , t. e. , lyne , a. g. 1985 , , 212 , 489 radhakrishnan , v. , cooke , d. j. 1969 ; aph .
lett , 3 , 225 rankin , j. m. 1993 , , 405 , 285 rathnasree , n. , rankin , j. m. 1995 , , 452 , 814 ruderman , m. a. , sutherland , p. g. 1975 , , 196 , 51 schnhardt , r. e. 1973 , nature phy .
sci , 243 , 62 schnhardt , r. e. 1974 , , 35 , 13 spitkovsky , a. 2006 , , 648 , l51 stappers , b. w. , gaensler , b. m. , johnston , s. 1999 , , 308 , 609 stinebring , d. r. , condon , j. j. 1990 , , 352 , 207 strom , r. g. , van someren greve , h. w. 1990 , , 171 , 351 swarup , g. , ananthakrishnan , s. , kapahi , v. k. , rao , a. p. , subrahmanya , c. r. , kulkarni , v. k. 1991 , current sci .
, 60 , 95 taylor , j. h. , manchester , r. n. , lyne , a. g. 1993 , , 88 , 529 weiler , k. w. , goss , w. m. , schwarz , u. j. 1974 , , 35 , 473