quantum chromodynamics is a perturbative field theory for @xmath10 , corresponding to energies above @xmath11 mev .
however , setting up a good perturbation expansion is more difficult than simply working order by order in @xmath12 .
the difficulty is that when computing a quantity with more than one scale , logarithms of ratios of those scales appear which invalidate the naive perturbation expansion .
for example , in @xmath0 collisions at center - of - mass energy @xmath4 , we might look for the distribution of jets as a function of the invariant mass @xmath13 of the jet .
then the differential cross section will have a correction of the form @xmath14 .
even for @xmath10 these large logarithms can dominate if @xmath13 is sufficiently small .
luckily , these logarithms appear only in certain combinations even at higher order , so that all the terms of the form @xmath15 can be ( re)summed at once .
however , understanding which terms will appear , how to resum them , and how to combine the resummed result with fixed order results can be quite complicated .
it is the goal of this paper to show how it can be done using effective field theory techniques .    in this paper , resummation and matching of event shapes
are studied using soft - collinear effective theory ( scet )  @xcite .
event shapes are observables which are sensitive to the overall distribution of final state particles , and therefore involve both short and long distance physics .
we will consider mainly the event shape @xmath16 where the thrust @xmath17 is defined by @xmath18 summing over all momentum @xmath19-vectors @xmath20 in the event and maximizing over unit 3-vectors @xmath21 . in the threshold region near @xmath22 large logarithms of the form @xmath23 appear at fixed order in perturbation theory .
the resummed result will be valid even if @xmath23 is large , as long as @xmath12 and @xmath24 are small .
effective field theories provide a systematic approach to resummation .
they separate out physics at a hard underlying scale @xmath4 from physics associated with a scale of interest @xmath25 and from even lower scales . at each scale
a separate matching calculation is done which is independent of physics at asymptotically lower or higher energy .
then the scale dependence is calculated using the renormalization group . in this way
, large logarithms can not appear because no two largely separated scales are accessible to the theory at the same time .
the soft - collinear effective theory works by separating the degrees of freedom of qcd into soft modes and collinear modes in different directions .
the relevant scales are then associated with the observable of interest , such as @xmath26 , or with the large components of the collinear modes .
the resummation of event shapes in @xmath0 is not currently of extreme importance phenomenologically .
it nevertheless provides a clean arena ( as compared to hadron collisions ) to explore resummation and matching .
the scet techniques and most of the formulas we discuss here were originally developed for @xmath27-physics , such as resummation in @xmath28 decays  @xcite .
they have also been applied to the study of deep - inelastic - scattering near @xmath29  @xcite and to the production of massive jets initiated by top quark decays  @xcite .
one convenient feature of effective field theories is factorization , which allows us to use objects , such as soft and jet functions , calculated in one process to study another .
thus , most of the hard work required for calculation of the distributions we describe here can be extracted from the literature .
nevertheless , there are certain features , in particular the nlo matching step , for which event shapes are uniquely illuminating .
the breakdown of naive perturbation theory due to the appearance of large logarithms is independent of @xmath12 blowing up and of non - perturbative effects . to emphasize this point , and to simplify the resummed expressions , the running of @xmath12
will be simply turned off ( setting @xmath30 ) for most of this paper . working in the conformal limit
removes one scale ( @xmath31 ) from the problem and thus clarifies which large logs are being resummed .
it is not hard to turn @xmath32 back on , as will be shown in section  [ sec : gen ] .
also , we will be including all one - loop results but no two - loop results .
thus , our expressions will not contain a complete next - to - leading log resummation , which should also resum the two - loop double logs .
there are three steps involved in the calculation of event shapes in scet , associated with three separate scales . at the hard scale @xmath33 ,
qcd is matched onto scet by demanding parton level matrix elements in the two theories be the same .
this can be done at leading order by matching onto an operator @xmath34 with two collinear fields , or at next - to - leading order by matching in addition to an operator @xmath35 with collinear fields in three directions  @xcite .
the threshold resummation only involves @xmath34 as the matrix elements of @xmath35 vanish near @xmath36 .
( an alternative approach to @xmath37 matching is described in  @xcite ) . at the scale @xmath38 , the collinear fields freeze and
can be removed from the theory by integrating them out .
this results in the a jet function @xmath39 .
finally , even though we are interested in a distribution at the scale @xmath40 , the soft degrees of freedom remain relevant down to a scale @xmath41 , after which they too can be integrated out of the theory .
the fixed order result will have large logarithms of @xmath42 and @xmath43 , but in the resummed result all these logs are exponentiated and innocuous .
before we present the factorization formula and calculate the thrust distribution in scet , we will review the way resummation of thrust is traditionally handled .
there are a number of ways to resum event shapes  @xcite .
since the focus of this work is on comparisons to scet , we will confine our attention to the original approach of catani , trentadue , turnock and webber  @xcite , which will be referred to as cttw throughout .
moreover , most of the other approaches reduce to  @xcite at next - to - leading fixed order in @xmath12 ( nlo ) and to leading log , so no significant loss of generality is sustained .
we will find that in the two - jet limit scet also agrees with cttw to the order we are working , although the resummed expressions are not exactly the same .
a more significant difference is in the nlo matching .
a critical advantage of the effective theory approach is that the finite parts of loop amplitudes are automatically incorporated into the perturbative expressions .
for example , the total nlo cross section for @xmath0 is reproduced by combining the finite parts of the hard , jet , and soft functions and a contribution from a finite integral over higher - order operators in scet .
this does not necessarily entail less work than in calculating the total cross section through the traditional combination of real and virtual contributions .
however , due to factorization , the infrared divergent contributions which are absorbed into the jet and soft functions are universal and thus the could potentially be used for many processes .
in this section , we will review some basic results from qcd on thrust , and the resummed expressions presented in  @xcite .    to begin , consider the parton model description of @xmath0 annihilation .
at order @xmath44 , the only process which contributes is @xmath45 .
these two jets have no structure and hence the cross section is simply @xmath46 , where @xmath47 is leading order total @xmath48 annihilation cross section .    at order
@xmath49 , there are two @xmath50 diagrams which contribute @xmath51 where the photon line on the bottom is the @xmath52 current coming in . the differential cross section is @xmath53_{\mathrm{parton } } = \delta ( { s } )    \delta ( t ) + { \bar{\alpha}}\frac{{s}^2 + { t}^2 + 2 { u}}{{s}\,{t } }   \label{sigparton}\ ] ] where we have defined @xmath54 and the reduced mandlestam variables are @xmath55 and @xmath56 with @xmath57 .
now , for 3-parton events the thrust variable @xmath58 reduces to @xmath59 so that @xmath60_{\mathrm{parton } } & = & \delta ( \tau ) +    { \bar{\alpha}}\left [ \frac{2 ( 3 \tau^2 - 3 \tau + 2)}{\tau ( 1 - \tau ) } \log    \frac{1 - 2\tau}{\tau } - \frac{3 ( 1 - 3 \tau ) ( 1 + \tau)}{\tau } \right]\\    & = & \delta ( \tau ) + { \bar{\alpha}}\left [ \frac{- 4 \log \tau - 3}{\tau } \right ] +    { \bar{\alpha}}d_{\mathrm{fin } } ( \tau ) .
\label{partondc}\end{aligned}\ ] ] the second line is written to manifest the singularity structure .
the remainder @xmath61 is finite as @xmath62 .    instead of the differential thrust distribution , it is useful to work directly with the integrated quantity @xmath63 for small @xmath64 , @xmath65 here we see explicitly the large logarithms @xmath23 and @xmath66 which demand resummation .
a useful theoretical trick , used in  @xcite and in  @xcite , is to employ a hemisphere mass definition .
this greatly simplifies the factorization formula in the two - jet limit . the hemisphere momentum @xmath67 is defined to be the sum of the 4-momenta of all the particles in the hemisphere of the left ( right ) jet .
then thrust reduces to @xmath68 , and we can calculate it from @xmath69 which we will see has a closed form expression .
the traditional approach to resummation which we review here is due to catani et al .
the basic idea is that terms of the form @xmath70 come from multiple real collinear and soft emissions .
the kinematics of multiple collinear emissions can be modeled by parton branchings , with branching probabilities proportional to the altarelli - parisi splitting functions @xmath71 in the two - jet limit , the corrections due to soft effects can be modeled by imposing angular ordering .
all of this follows from the coherent - branching algorithm developed in  @xcite .
in fact , more general methods are available for handling soft emissions ( for example ,  @xcite ) .
but , to nlo for event shapes like thrust , the methods set up in  @xcite are sufficient to compare to scet .    in the two - jet limit , the factorized expression for the differential cross section in cttw is @xmath72_\mathrm{cttw } = { j_{\mathrm c } } ( q^2 , p_l^2 ) { j_{\mathrm c } } ( q^2 , p_r^2 ) .
\label{factqcd}\ ] ] the subscript on @xmath73 refers to the 2-jet contribution .
@xmath74 has a simple physical interpretation : it is the probability of finding a final state jet with invariant mass @xmath13 .
thus , it satisfies @xmath75 @xmath4 is a scale associated with the hard process .
we will see later that the equivalent jet function in scet has a different normalization .    at leading order the jet is a massless parton and so @xmath76 . at next - to - leading order ,
the coherent branching algorithm allows one - angular ordered emission .
so @xmath77 here @xmath78 is the energy fraction in the emission and @xmath79 .
the angular ordering constraint for the coherence of soft emissions is implicit in the restriction @xmath80 ; for @xmath81 it cuts off the infrared divergences in the integral as @xmath82 .
evaluating   gives @xmath83_{\star}^{[p^2 , q^2 ] } + \cdots .
\label{jfirst}\ ] ] thus , to first order in @xmath12 , @xmath84_\mathrm{cttw } = \delta ( \tau ) + { \bar{\alpha}}\left [ \frac{- 4 \log \tau    - 3}{\tau } \right]_{\star}^{[\tau , 1 ] } ,   \label{fixedr}\ ] ] which reproduces the divergent as @xmath85 part of the the parton model result  .
here we have introduced the @xmath86 distribution ( alternatively called the r- or @xmath87-distribution ) , which is a generalization of the + -distribution .
it is uniquely defined by the two conditions @xmath88^{[x , a]}_{\star } = f ( x)$ for $ x > 0$}\ ] ] @xmath89^{[x , a]}_{\star } g ( x ) = \int_0^a d x f ( x ) \left [ g    ( x ) - g ( 0 ) \right]\ ] ] for clarity , we have added to the notation an explicit instance of the dependent variable @xmath90 in @xmath91_{\star}^{[x , a]}$ ] .
a useful relation is @xmath92^{[x , a]}_{\star } = [ f ( x)]^{[x , b]}_{\star }   + \delta ( x ) \int_a^b d   x ' f ( x ' ) .\ ] ] this will be used extensively in the next section .
the leading order resummation is preformed by iterating the angular - ordered emissions .
this leads to an integro - differential equation for @xmath93 , similar to the evolution equation for parton - distribution functions .
the equation is solved in an integrated form , and the resummed expression for the integrated thrust in the two jet limit is @xmath94_\mathrm{cttw } = \exp \left [ - 2 { \bar{\alpha}}\log^2 \tau - 3{\bar{\alpha}}\log \tau   \right ]   \frac{e^{- 2 \gamma_e \eta}}{\gamma [ 2 \eta+1 ] } ,   \label{qcdresum}\ ] ] where @xmath95 this equation is a combination of expressions in  @xcite taken with @xmath30 .
expanding @xmath96 to order @xmath12 the fixed order result   is reproduced .
this expression is resummed in the sense that @xmath97 as @xmath85 , in contrast to the fixed order result   for @xmath98 or the parton model expression   which diverge as @xmath62 .
having reviewed the way the resummed thrust distribution is traditionally calculated , we now turn to the equivalent calculation in scet .
we will see that there are a number of advantages of this effective field theory treatment . instead of a classical treatment , where multiple real emissions are summed at the level of the cross section
, the entire resummation is done in scet through the renormalization group .
this makes explicit the various scales in the problem , and allows greater freedom to choose those scales to minimize the large logarithms in the distributions of interest .
the idea behind scet is to separate out the quarks and gluons of qcd into soft and collinear degrees of freedom .
a collinear field is associated with a light - like direction @xmath99 .
the component of its momentum in that direction , @xmath100 , must be much larger than any of the other momentum components .
all the qcd degrees of freedom which could change this momentum have been integrated out , so @xmath101 appears as a fixed label .
for example , a collinear quark is written as @xmath102 with a fixed @xmath101 .
in addition to collinear quarks and gluons , there are soft quarks and gluons .
these soft fields can only interact with each other or transfer momentum to the soft components @xmath103 of a collinear field .
the lagrangian of scet has a separate gauge invariance associated with soft gluons and gluons in each collinear direction .
it is useful to maintain this gauge invariance explicitly in the operators of the theory with the use of wilson lines .
for example , a two - jet operator is @xmath104 where @xmath105 @xmath106 the collinear wilson lines @xmath107 maintain collinear gauge invariance and the soft wilson lines @xmath108 maintain soft gauge invariance .
the starting point in the effective theory approach to event shapes is again factorization in the two - jet limit . in scet ,
the event - shape distributions near the endpoint come from matrix elements of @xmath34  @xcite . in terms of the hemisphere masses @xmath109 and @xmath110 defined above , factorization implies  @xcite @xmath111_\mathrm{scet}= |c_h ( \mu)|^2 \int { \mathrm{d}}k_l { \mathrm{d}}k_r j    ( p_l^2 - q k_l , \mu ) j ( p_r^2 - q k_r , \mu ) s ( k_l , k_r , \mu ) .
\label{factorform}\ ] ] the hard function @xmath112 , the jet functions @xmath113 and the soft function @xmath114 will be defined below .
the form of this factorized expression has a physical explanation .
each jet function @xmath115 comes from one of the collinear quarks .
it represents , like @xmath74 , the probability for producing a jet of invariant mass @xmath13 from that collinear field ( the precise relation to @xmath74 is explored below ) .
recall that the large component of the momentum of a collinear field @xmath116 is fixed .
since , in the two - jet limit , the jets are back to back with center of mass energy @xmath4 , we must have @xmath117
. the small component @xmath118 can vary .
if there are no soft effects , then the hemisphere mass is simply @xmath119 .
however , as the factorization formula implies , the collinear jet can give up some soft momentum to the soft qcd background , leading to @xmath120 .
the hemisphere mass is unchanged by this emission , but now @xmath121 , which explains the form of .
the function @xmath122 is a hard function , which comes from integrating out hard modes of qcd in matching to scet .
demanding @xmath123 for states with two external quarks lets us calculate @xmath122 order by order in perturbation theory .
the computation entails taking the difference between qcd and scet graphs , such as @xmath124 where the @xmath125 refers to an insertion of @xmath126 , and the right diagram is only a representative contribution ( for details , see  @xcite ) .
the difference is finite because the infrared divergences in qcd and scet are the same and the uv divergences are removed with counterterms . at one - loop
the matching gives  @xcite @xmath127 with @xmath128 , @xmath129 and @xmath130 .
the anomalous dimension and renormalization group evolution of @xmath112 can be found from the uv divergences of the same 1-loop graphs .
the rg equation is @xmath131 with solution @xmath132 \left ( \frac{q^2}{\mu_h^2}\right)^{-\gamma_h   \log\frac{\mu^2}{\mu_h^2 } } .\ ] ] to first order in @xmath12 , the @xmath133 dependence drops out and is reproduced .
the natural matching scale @xmath134 is manifest in this expression .
the jet function @xmath115 is the imaginary part of the propagator of a collinear quark .
it is defined by  @xcite @xmath135 .\ ] ] it can be thought of as the spectral density for a jet of collimated particles interacting with a soft qcd background .
the jet function can be calculated order - by - order in perturbation theory through the discontinuity of conventional feynman diagrams , such as @xmath136 some of the infrared divergences in these graphs are cut off because the collinear quark is taken to have invariant mass @xmath13 .
these calculations have been done in the context of @xmath28  @xcite and for deep - inelastic scattering  @xcite and the jet function is known to two loops  @xcite . due to factorization , the same jet function applies in these processes and in @xmath0 annihilation . to first order in @xmath12 ,
the result is @xmath137 + \left [ \frac{\gamma_j    \log \frac{p^2}{\mu^2 } + \gamma_j}{p^2 } \right]^{[p^2 , \mu^2]}_{\star } ,   \label{jetfixed}\ ] ] with @xmath138 , @xmath139 and @xmath140 .
the renormalization group evolution of the jet function , in contrast to that of the hard function , is non - local in @xmath13 @xmath141 j ( p^2 , \mu ) + 2\gamma_j \int_0^{p^2 } d    q^2 \frac{j ( p^2 , \mu ) - j ( q^2 , \mu)}{p^2 - q^2 } .\ ] ] this is similar to the altarelli - parisi equation for the evolution of the parton distribution functions ( pdfs ) .
in contrast to pdfs , however , the jet function @xmath115 is perturbative as long as @xmath142 .
simplification is achieved through use of the laplace transform  @xcite , @xmath143 whereby the evolution becomes local in @xmath144 @xmath145 now the rge can be solved like that of @xmath122 : @xmath146   \left ( \nu e^{\gamma_e } \mu_j^2    \right)^{-\eta_j } \tilde{j}(\nu,\mu_j)\ ] ] with @xmath147 finally , the inverse laplace transform produces  @xcite @xmath148 \tilde{j } \left (    \partial_{\eta_j } \right ) \left [ \frac{1}{p^2 } \left ( \frac{p^2}{\mu_j^2 }    \right)^{\eta_j } \right]_{\star}^{[p^2 , \mu_j^2 ] } \frac{e^{- \gamma_e    \eta_j}}{\gamma [ \eta_j]},\ ] ] where @xmath149 the functional dependence on @xmath150 in @xmath151 comes from a functional dependence on @xmath152 in the fixed order expression for the jet function   at the matching scale @xmath153 .    finally , we have to calculate the soft function @xmath154 .
it is defined though matrix elements of soft wilson lines  @xcite @xmath155 unlike the jet function , the soft function is often evaluated at scales @xmath156 , where it is non - perturbative .
it can , in principle , be measured experimentally and then be evolved perturbatively , like the pdfs . or it can be modeled  @xcite .
however , in our case , since we are taking @xmath12 to be small and fixed , it can be calculated in perturbation theory .
even if @xmath157 is non - perturbative , the perturbative calculation is useful because the uv divergences dictate the anomalous dimensions and hence the evolution equation .
the soft function , being completely defined in terms of wilson lines , can be calculated in qcd through diagrams such as @xmath158    where the kinked lines are wilson lines in the @xmath159 and @xmath160 directions . to order @xmath12 ,
@xmath161 can be written as a product @xmath162 where  @xcite @xmath163   + \left [   \frac{\gamma_s \log \frac{k}{\mu } + \gamma_s}{k } \right]^{[k , \mu]}_\star \label{sfirst } , \ ] ] with @xmath164 , @xmath165 and @xmath166 .
the evolution of the soft function involves the same kind of non - local equation as for the jet function .
this can be seen from examining its anomalous dimension directly , or from rg invariance of the convolution appearing in the factorization formula  .
the solution is therefore similar .
explicitly , @xmath167 \tilde{s } ( \partial_{\eta_s } ) \left [    \frac{1}{k } \left ( \frac{k}{\mu_s } \right)^{\eta_s } \right]_{\star}^{[k ,    \mu_s ] } \frac{e^{- \gamma_e \eta_s}}{\gamma [ \eta_s ] } , \ ] ] where @xmath168 and @xmath169 expanding to first order in @xmath12 reproduces  .
now we have all the ingredients appearing in the the factorization formula  . at order @xmath12 , the result is @xmath170_\mathrm{scet } = \delta ( p^2_l ) \delta ( p_r^2 ) \left [    1 + { \bar{\alpha}}\left ( - 1 + \frac{\pi^2}{3 } \right )   \right ] \nonumber\ ] ] @xmath171^{[p^2_l , q^2]}_\star   + { \bar{\alpha}}\delta ( p_l^2 )   \left
[ \frac{2 \log \frac{q^2}{p^2_r }   -   \frac{3}{2}}{p^2_r } \right]^{[p^2_r , q^2]}_\star .
\label{dplpr}\ ] ] note that the dependence on all the scales @xmath172 and @xmath173 has completely canceled .
this is a non - trivial result which requires three relations among the six anomalous dimensions : @xmath174 and is a strong consistency check on the entire formalism .
s can be understood on more general grounds by relating the @xmath175s to a universal cusp anomalous dimension  @xcite . ]    already scet can be compared to cttw .
the differential distribution in scet , eq .
, matches the cttw distribution , eq .
, when the jet functions are expanded to first order using eq .  .
the scet expression has an additional finite piece ( the @xmath176 term ) which comes from loop graphs which do not enter into the traditional formulation .    still working at leading order in @xmath12
, scet produces a simple form for the thrust distribution near @xmath177 : @xmath178_\mathrm{scet } = \delta ( \tau ) \left [ 1 + { \bar{\alpha}}\left ( - 1 + \frac{\pi^2}{3 } \right ) \right ]   + { \bar{\alpha}}\left [ \frac{- 4 \log \tau - 3}{\tau } \right]^{[\tau , 1]}_\star . \label{scetsig2}\ ] ] this reproduces the leading behavior for small @xmath64 of the both the parton model expression , eq .  , and the resummed expression with the functions @xmath74 expanded to first order , eq .  .
to get the resummed thrust distribution from scet , we need to calculate a couple of convolution integrals .
first , the soft and jet functions must be combined into @xmath179 second , the @xmath180 functions must be integrated against the event shape . for thrust , we need @xmath181 both of these convolutions can be evaluated by performing the laplace transform .
note that the function names can be misleading . here ,
@xmath182 ( and not just the scet jet function @xmath115 ) is playing the role of @xmath183 from the cttw formulation .
thus , the cttw jet function is a combination of the soft and jet functions in scet , as expected because the coherent branching algorithm used to derive it incorporates both soft and collinear effects .    solving the @xmath180 function with the laplace transform techniques gives @xmath184 \exp \left [ \frac{\gamma_s}{2 }    \log^2 \frac{\mu}{\mu_s } - \gamma_s \log \frac{\mu}{\mu_s } \right ] \left (    \frac{q \mu_s}{\mu_j^2 } \right)^{-\eta_s}\ ] ] @xmath185 \left [ \frac{1}{p^2 } \left ( \frac{p^2}{\mu_j^2 }    \right)^{\eta_k } \right]_{\star}^{[p^2 , \mu_j^2 ] } \frac{e^{- \gamma_e    \eta_k}}{\gamma [ \eta_k ] } ,   \label{resid}\ ] ] where @xmath186 and @xmath187 note the residual logarithms of @xmath188 in line which are not resummed .
these imply that we can not choose @xmath189 .
instead , the natural scales ( those which remove the large logs ) should satisfy @xmath190 .
then , we can evolve either the jet function from @xmath191 down to @xmath173 or the soft function from @xmath173 up to @xmath191 but there is not a single scale which can minimize all the large logs .    finally , to evaluate the convolution for the thrust integral , we use the laplace transform of the @xmath192-function @xmath193 combining this with the expression for @xmath194 above , the scet prediction for thrust in the 2-jet limit is @xmath94_\mathrm{scet }   = |c_h ( \mu_h ) |^2 \exp \left [ 2 { \bar{\alpha}}\log^2    \frac{\mu_h^2}{\mu_j^2 } + 3 { \bar{\alpha}}\log \frac{\mu_h^2}{\mu_j^2 } - { \bar{\alpha}}\log^2 \frac{\mu_h^2}{\mu_s^2 } \right ] \left ( \frac{q \mu_s}{\mu_j^2 }    \right)^{- 4 { \bar{\alpha}}\log \frac{\mu_h^2}{\mu_s^2 } } \nonumber\ ] ] @xmath195 ^ 2 \left (    \frac{q^2 \tau}{\mu_j^2 } \right)^{2 \eta_k } \frac{e^{- 2 \gamma_e    \eta_k}}{\gamma [ 2 \eta_k + 1 ] } .\ ] ] note that this expression is explicitly independent of @xmath87 .
now , choosing the natural scales @xmath196 , @xmath197 , and @xmath198 to remove the large logs , the thrust distribution becomes @xmath94_\mathrm{scet }    = \exp \left [ - 2 { \bar{\alpha}}\log^2 \tau - 3 { \bar{\alpha}}\log \tau    \right ] \tilde{r } ( \partial_{\eta } ) \frac{e^{- 2 \gamma_e \eta}}{\gamma    [ 2 \eta + 1]}\ ] ] with @xmath199 and @xmath200 as in eq .  .
this expression would be identical to the cttw expression , eq .
, if @xmath201 .
recall that the function @xmath202 comes from 1-loop matching in scet , which turns into boundary conditions for the renormalization group evolution . in the approach of cttw
the boundary condition is simply that @xmath203 . nevertheless , the effect of @xmath204 to order @xmath12 is only to provide a finite constant , as can been seen from comparing equations   and
that is , there is no contribution of order @xmath205 in the difference , only of order @xmath206 and higher . beyond order @xmath12 ,
the @xmath207 does change the @xmath205 dependence .
but since @xmath208 these terms are subleading to the dominant @xmath209 and can get corrections from higher - loop effects .
for example , at two - loops , a term of the form @xmath210 is generated .
thus , scet and cttw agree to first order in @xmath12 and for the resummation of the leading large logarithms , which is the order to we have been working .
now let us look at how fixed order results and resummation are combined .
there are two elements to this : ( 1 ) matching to hard emissions to get the differential distribution correct away from the two - jet region ( i.e. away from @xmath22 ) ; and ( 2 ) including finite parts of loops to reproduce fixed order inclusive results .
let us begin with the matching procedure described in  @xcite . to match to the hard emissions at order @xmath12 we need the parton model differential cross section from equation ( [ partondc ] ) .
the divergent part is already contained in the two jet contribution , so the remainder is @xmath211 @xmath212 - 2 \log^2 ( 1 - \tau )     + 3 ( 1 - 2 \tau ) \log ( 1 - 2 \tau )   - 4 \mathrm{li}_2 (   \frac{\tau}{1 - \tau } )   + 6\tau +   \frac{9}{2 } \tau^2 .\ ] ] the final result for the matched integrated thrust distribution from  @xcite is @xmath213_{\mathrm{cttw } }   = ( 1 + \sigma_1 )   \left\ {   \left ( 1 + c_1 \right ) \exp \left [ - 2 { \bar{\alpha}}\log^2 \tau -    3{\bar{\alpha}}\log \tau \right ] \frac{e^{- 2\gamma_e \eta}}{\gamma [ 2\eta+1 ] }    + { \bar{\alpha}}{d_\mathrm{fin}}(\tau ) \right\ } ,   \label{matchedqcd}\ ] ] where @xmath214 is the nlo contribution to the total cross section and @xmath215 is chosen so that @xmath216 ( recall that @xmath217 at order @xmath12 ) . in the qcd approach the factor @xmath218
must be be determined by a independent calculation of the total rate .
this involves combining the real and virtual contributions at order @xmath12 using an appropriate infrared regulator . as emphasized in  @xcite
there is arbitrariness in the matching because the hard emissions ( in @xmath219 ) are fixed order but the two - jet contribution is resummed .
for example , to the same accuracy @xmath220 could be multiplied by the exponential .
now , let us turn to the matching in scet . to include thrust distributions away from the endpoint , we can either attempt to add power corrections to scet , or we can match to higher order operators as described in  @xcite . matching is much simpler . to perform the matching , we add new operators @xmath221 where the @xmath222 are the additional terms coming from wilson lines necessary for gauge invariance .
@xmath223 is a collinear gluon in direction @xmath224 .
the matching demands that @xmath225 where the subscript means the matching is done at the hard scale @xmath133 .
there is some arbitrariness in the definition of the matrix elements in scet due to reparameterization invariance .
the matrix element of a collinear quark @xmath102 on a qcd quark state @xmath226 is only defined up to its soft momentum component @xmath227 . moreover
, a basis for summing over directions @xmath159 for @xmath126 and @xmath228 , @xmath224 , and @xmath229 for @xmath37 must be chosen as well .
a certain convention was described in  @xcite for resolving these ambiguities and others are possible . in any case , while different conventions may shift the contributions from @xmath34 and @xmath35 in  , the sum is parameterization invariant .
thus , independently of the convention we have @xmath230_{\mathrm{scet}}^{\mu_h } = [ { \mathrm{d}}\sigma]_{\mathrm{parton}},\ ] ] where @xmath231_{\mathrm{parton}}\sim \langle \mathrm{qcd } \rangle^2 $ ] is the tree - level parton model cross section , as shown in eq .  .
now , we already know that @xmath232 gives   to first order in @xmath12 and that the parton model distribution at order @xmath12 is  .
thus , with obvious implicit phase space factors , @xmath233 so , at leading order , the contribution from the sum of @xmath234 and the interference between @xmath34 and @xmath35 to the differential cross section is unambiguous .
the running of @xmath35 could also be included even though it does not resum any large logs for the event shapes under consideration . with running , the matrix elements of @xmath34 and @xmath35 would appear with separate renormalization kernels , and so the final differential cross section would end up depending on the conventions chosen .
the ambiguity could be resolved by a careful higher order treatment , but for the purposes of comparing to the cttw prediction for thrust , we will simply not renormalize the finite terms .
then , @xmath213_{\mathrm{scet } }   = \exp \left [ - 2    { \bar{\alpha}}\log^2 \tau - 3 { \bar{\alpha}}\log \tau \right ] \tilde{r } ( \partial_{\eta } )    \frac{e^{-2 \gamma_e \eta}}{\gamma [ 1 + 2\eta ] } + { \bar{\alpha}}{d_\mathrm{fin}}(\tau ) .
\label{matchedscet}\ ] ] with @xmath235 and @xmath205 as in eq .   and
@xmath219 in  .
the total cross section in scet is given @xmath236_\mathrm{scet } =   \sigma_0\left[1 + { \bar{\alpha}}\left ( - 1 + \frac{\pi^2}{3 } \right ) + { \bar{\alpha}}{d_\mathrm{fin } } ( \frac{1}{3 } ) \right ] = \sigma_0\left [ 1 + \frac{\alpha_s}{\pi}\right].\ ] ] this is the correct total @xmath0 total cross section to first order in @xmath12 !
let us review the contributions that go into the cross section .
first , at the hard scale @xmath134 , there is the finite part of the loop matching to qcd , @xmath237 .
next , at the scale @xmath238 where we integrate out the collinear fields , the jet functions give @xmath239 .
at the seesaw scale @xmath240 , the soft function gives @xmath241 , and finally the finite part of the real emission integral , away from @xmath22 gives @xmath242 . in producing the total cross section ,
only the soft and jet contributions are infrared divergent .
however , their convolution , which appears in the function @xmath243 is infrared finite . the hard matching and the hard emissions
are ir finite by themselves .
thus , the total cross section can be understood as a combination of a process dependent ir finite hard part and universal but ir regulator dependent contributions from soft and collinear emissions .
the scet thrust distribution   is compared to the cttw thrust distribution   in figure  [ fig : dt ] .
the plot is of @xmath244 with @xmath12 fixed at @xmath245 .
a more careful rendition of the differential thrust distribution would take the derivative of @xmath246 before assigning the matching scale @xmath247 , as is done below .
however , the effect is higher order , and so we just plot @xmath248 directly .
matched resummed differential thrust distribution in scet and in the standard approach of cttw with fixed coupling @xmath249,scaledwidth=95.0% ]
in this section , some simple generalizations are described .
the above results were derived assuming @xmath12 to be constant in order to emphasize the resummation of sudakov logarithms in contrast to large logarithms associated with the scale @xmath31 .
now , it is shown how the results can be modified with running @xmath12 . also , the scet prediction for another event shape , the jet mass @xmath250 , is given .
it is straightforward to allow @xmath12 to run . including 1-loop running , the effect is to modify the single and double logs in the following way  @xcite : @xmath251\ ] ] @xmath252 for example
, the differential thrust distribution for @xmath253 with running @xmath12 becomes @xmath254_{\mathrm{scet } } \!\!\!\!&=&\ !
\frac{1}{\tau } \exp\left [ 4
s(q , q \tau ) + 6 a(q , q \tau )   - 8 s(q\sqrt{\tau},q \tau ) -6 a(q \sqrt{\tau},q \tau)\right ]   \tilde{r } ( \partial_{\eta } ) \frac{e^{- 2 \gamma_e \eta}}{\gamma [ 2 \eta ] } \nonumber \\ & &   + \frac{2 \alpha_s}{3 \pi } d_{\mathrm{fin } } ( \tau ) ,   \label{scetfinal}\end{aligned}\ ] ] where @xmath255 @xmath256 is given in eq .   and @xmath257 @xmath258 \partial_{\eta}^2   - \frac{3}{2 } \alpha_s(q\sqrt{\tau } )    \partial_{\eta } \label{r2 } \right\ } .\ ] ]
this is the same function @xmath202 as in eq .  , but with the @xmath12 factors evaluated at the appropriate matching scales .
at this point , one would hope to compare to data .
however , besides the obvious shortcoming of not containing the full nll resummation ( it does not include effects of the two - loop cusp anomalous dimension ) , this parton - level expression does not include important non - perturbative effects .
due to the running of @xmath12 , the perturbative expression breaks down when the soft scale is of order @xmath31 , that is , when @xmath259 , as can be seen explicitly in  . in fact , even for significantly larger values of thrust power corrections of order @xmath260 become quantitatively important , at least at lep energies .
this problem has been approached elsewhere using scet  @xcite and with other techniques  @xcite .
other event shapes can be studied the same way as thrust .
for example , consider the heavy jet mass @xmath250 defined by @xmath261 in this case , the matching scales are @xmath134 , @xmath262 and @xmath263 and scet gives for @xmath264 @xmath265_\mathrm{scet } & = & \frac{2}{\rho } |1 + c_h |^2   \exp\left [ 4 s(q , q \rho ) + 6 a(q , q \rho )   - 8 s(q\sqrt{\rho},q \rho ) - 6 a(q \sqrt{\rho},q \rho ) \right ] \nonumber\\ & &   \times \left [ \tilde{k }    ( \partial_{\eta } ) \frac{e^{- \gamma_e \eta}}{\gamma [ \eta ] } \right ] \left [    \tilde{k } ( \partial_{\eta } ) \frac{e^{- \gamma_e \eta}}{\gamma [ \eta + 1 ] }    \right ] + \frac{2 \alpha_s}{3 \pi } d_{\mathrm{fin } } ( \rho ) , \end{aligned}\ ] ] where @xmath266 this formula agrees with the jet mass distribution in  @xcite to leading log and first order in @xmath12 .
the same function @xmath267 appears for jet mass and for thrust because to order @xmath12 in the parton model , @xmath268 .
we have investigated how to combine resummation with next - to - leading order matching of event shapes in the original approach of  @xcite ( cttw ) and using scet . in the cttw formulation ,
real emissions from collinear splitting functions are used and various kinematical features associated with soft emission are combined to derive a differential cross section .
the cross section factorizes into the product of two jet functions .
resummation is done by solving a differential equation for the jet functions in terms of the physical scales @xmath13 and @xmath269 of the event .
in contrast , scet factorizes the event shape distribution into a contribution from hard , jet , and soft functions .
these functions are matched at the scales @xmath196 , @xmath270 and @xmath271 respectively and resummation is done through renormalization group evolution .
the resummed thrust distribution in scet and cttw have been compared , and found to agree to next - to - leading order in @xmath12 and for leading - log resummation .
the resummation of thrust in scet brings to light a number of features of the process not obvious in cttw formulation .
for example , the appearance of the seesaw scale @xmath272 as the natural matching scale for the soft function follows from the kinematics of the scet factorization theorem .
of course , the existence of this scale has been known for a long time from qcd , but in the effective field theory approach this scale just drops out of the factorized expression .
thus , for more complicated processes , it is reasonable to expect a similar transparency for the matching scales , which may facilitate subleading log resummation .
in fact , two - loop , three - loop , and some four - loop anomalous dimensions for various quantities are already available  @xcite , and so subleading log resummation appears quite feasible . the biggest impediment to using these more accurate resummed results in a comparison to data
is that power corrections of order @xmath260 have an important quantitative effect on event shapes .
however , these corrections should modify only the soft function while higher order resummed expressions for the the hard and jet functions will remain valid .
thus , the effective theory should be able to weave together the perturbative and non - perturbative calculations .
a new result of this paper is the demonstration that inclusive quantities , such as the total cross section for @xmath0 can be calculated in a new way using scet . instead of taking the full differential @xmath159 + 1 body cross section and combining with the one - loop @xmath159-body cross section
, scet combines finite parts of loops of soft and jet functions with a hard matching calculation and a finite integral over hard emissions .
the soft and jet functions depend on the infrared regulator , but their convolution , and the hard function , do not . for @xmath0 annihilation at nlo
, this may not be so impressive , but the procedure promises to apply to more complicated processes , perhaps even some for which nlo results are not available .
it would also be very interesting to explore this mechanism at nnlo or to work with hadronic processes where the singularity structure is more complicated .
i would like to thank t.  becher and s.  fleming for many useful discussions .
i also benefited from the hospitality of the les houches ecole de physique and discussions with many of its participants .
this work was supported in part by the national science foundation under grant nsf - phy-0401513 and by the johns hopkins theoretical interdisciplinary physics and astronomy ceneter .      c.  w.  bauer , s.  fleming and m.  e.  luke , phys .
d * 63 * , 014006 ( 2001 ) [ arxiv : hep - ph/0005275 ] . c.  w.  bauer , s.  fleming , d.  pirjol and i.  w.  stewart , phys .
d * 63 * , 114020 ( 2001 ) [ arxiv : hep - ph/0011336 ] . c.  w.  bauer , d.  pirjol and i.  w.  stewart , phys .
d * 65 * , 054022 ( 2002 ) [ arxiv : hep - ph/0109045 ] .
m.  beneke , a.  p.  chapovsky , m.  diehl and t.  feldmann , nucl .
b * 643 * , 431 ( 2002 ) [ arxiv : hep - ph/0206152 ] .      c.  w.  bauer and a.  v.  manohar , phys .
d * 70 * , 034024 ( 2004 ) [ arxiv : hep - ph/0312109 ] .
s.  w.  bosch , b.  o.  lange , m.  neubert and g.  paz , nucl .
b * 699 * , 335 ( 2004 ) [ arxiv : hep - ph/0402094 ] .
m.  neubert , eur .
j.   c * 40 * , 165 ( 2005 ) [ arxiv : hep - ph/0408179 ] .
a.  v.  manohar , phys .
114019 ( 2003 ) [ arxiv : hep - ph/0309176 ] . t.  becher and m.  neubert , phys .
* 97 * , 082001 ( 2006 ) [ arxiv : hep - ph/0605050 ] . t.  becher , m.  neubert and b.  d.  pecjak , jhep * 0701 * , 076 ( 2007 ) [ arxiv : hep - ph/0607228 ] .
s.  fleming , a.  h.  hoang , s.  mantry and i.  w.  stewart , arxiv : hep - ph/0703207 .
c.  w.  bauer and m.  d.  schwartz , arxiv : hep - ph/0607296 . c.  w.  bauer and m.  d.  schwartz , phys .
lett .   * 97 * , 142001 ( 2006 ) [ arxiv : hep - ph/0604065 ] .
m.  trott , phys .
d * 75 * , 054011 ( 2007 ) [ arxiv : hep - ph/0608300 ] .
s.  catani , g.  turnock , b.  r.  webber and l.  trentadue , phys .
b * 263 * , 491 ( 1991 ) .
s.  catani , m.  l.  mangano , p.  nason and l.  trentadue , nucl .
b * 478 * , 273 ( 1996 ) [ arxiv : hep - ph/9604351 ] .
y.  l.  dokshitzer , a.  lucenti , g.  marchesini and g.  p.  salam , jhep * 9801 * , 011 ( 1998 ) [ arxiv : hep - ph/9801324 ] .
s.  catani , b.  r.  webber and g.  marchesini , nucl .
b * 349 * , 635 ( 1991 ) .
s.  catani and m.  h.  seymour , nucl .
b * 485 * , 291 ( 1997 ) [ erratum - ibid .
b * 510 * , 503 ( 1998 ) ] [ arxiv : hep - ph/9605323 ] . c.  w.  bauer , a.  v.  manohar and m.  b.  wise , phys .
lett .   * 91 * , 122001 ( 2003 ) [ arxiv : hep - ph/0212255 ] . c.  lee and g.  sterman , phys .
d * 75 * , 014022 ( 2007 ) [ arxiv : hep - ph/0611061 ] .
t.  becher and m.  neubert , phys .
b * 637 * , 251 ( 2006 ) [ arxiv : hep - ph/0603140 ] .
g.  p.  korchemsky and g.  marchesini , phys .
b * 313 * , 433 ( 1993 ) .
g.  p.  korchemsky and s.  tafat , jhep * 0010 * , 010 ( 2000 ) [ arxiv : hep - ph/0007005 ] . g.  p.  korchemsky and g.  sterman , nucl .
b * 555 * , 335 ( 1999 ) [ arxiv : hep - ph/9902341 ] .
g.  p.  korchemsky and a.  v.  radyushkin , nucl .
b * 283 * , 342 ( 1987 ) .
a.  v.  belitsky , phys .
b * 442 * , 307 ( 1998 ) [ arxiv : hep - ph/9808389 ] .
s.  moch , j.  a.  m.  vermaseren and a.  vogt , nucl .
b * 688 * , 101 ( 2004 ) [ arxiv : hep - ph/0403192 ] .
t.  becher and m.  neubert , phys .
b * 633 * , 739 ( 2006 ) [ arxiv : hep - ph/0512208 ] .


it is well established that the diffuse x - ray emission from clusters of galaxies implies the existence of a hot ( 10@xmath1 k ) , thermally radiating intracluster medium ( icm ) . in the dense centers of many clusters
this gas is expected to cool radiatively in less than a hubble time .
as pressure support decreases the gas will gravitate towards the bottom of the cluster potential well , a phenomenon that has come to be known as a `` cooling flow '' ( cf .
fabian 1994 ) .
such cooling flows can be detected by direct x - ray imaging , where they are seen as a surface brightness enhancement in the x - ray emission from a region typically @xmath0 100 - 200 kpc in radius , and by spectroscopic observations which reveal gas at cooler temperatures than the surrounding icm .
the excess x - ray luminosity in the cooling flow region is typically 10@xmath2 ergs s@xmath3 , corresponding to roughly 10 - 50 % of the total cluster x - ray emission ( peres et al .
the mass deposition rate is @xmath4 , where @xmath5 is the gas temperature at r@xmath6 , the cooling radius .
estimated cooling flow rates in the centers of rich clusters vary from 30 - 300 m@xmath7 yr@xmath3 ( white , jones , & forman 1997 ) and in some cases reach mass deposition rates of 1000 m@xmath7 yr@xmath3 or more ( allen 2000 ) .
a long - standing problem with the cooling flow scenario has been the fact that the large quantities of cool gas expected to accumulate from the flow have rarely been detected .
searches for cool gas in atomic or molecular form have generally proven unsuccessful ( e.g. , heckman et al .
1989 ; jaffe 1990 ; mcnamara , bregman , & oconnell 1990 ; odea et al .
1994 ; voit & donahue 1995 ; jaffe & bremer 1997 ; donahue et al .
2000 ) , with reported upper limits two orders of magnitude smaller than the masses implied by x - ray mass deposition rates ( although see edge 2001 ) . while some support for the cooling flow picture is provided by low energy x - ray spectra of clusters which appear to require larger absorption than can be accounted for by foreground gas alone ( white et al .
1991 ; allen et al . 1993
; irwin & sarazin 1995 ) , it now appears that the fraction of clusters that actually possess such excess absorption may be much less than previously thought ( arabadjis & bregman 2000 , hereafter ab2000 ; bhringer et al .
furthermore , recent high - resolution x - ray spectra obtained with the xmm - newton and chandra satellites have failed to show expected emission lines from cooled gas at temperatures below a few kev ( fabian et al . 2001 ; peterson et al . 2001 ; tamura et al . 2001 ) and it has also been found that local isothermality matches the spectral data of `` cooling flow '' regions better than the continually declining temperature profile predicted by the standard cooling flow model ( e.g. molendi & pizzolato 2001 ; bhringer et al .
a variety of mechanisms have been proposed to account for the lack of detectable cool gas in cooling flows ; these include reheating the gas , mixing , differential absorption , efficient conversion of cooling flow gas into low mass stars , inhomogeneous metallicity distributions , or disruption of cooling flows by recent subcluster mergers ( see recent discussions in fabian et al .
2001 and peterson et al .
2001 ) . as of
yet , it is unclear whether any of these mechanisms are capable of resolving the aforementioned problems with the cooling flow model .
one characteristic that is common to the mechanisms mentioned above is that they have mainly been introduced in an `` after the fact '' manner .
that is , they assume the observed mass deposition rates are basically correct but that something is happening to the gas after ( or while ) it cools , be it reheating or conversion into low mass stars or some other process .
perhaps a simpler way of addressing the problems of the cooling flow model are to look for ways to reduce the mass deposition rates from the beginning .
one way to achieve this is to reduce the excess emission ( @xmath8 ) associated with cooling flows by attributing some ( or all ) of it to x - ray sources other than thermal bremsstrahlung or line emission from the icm . but are there sources luminous enough to reproduce @xmath8 , can they be found within galaxy clusters , and , most importantly , can they mimic hot cluster gas ?
soon after x - ray emission was first detected from galaxy clusters , katz ( 1976 ) proposed that much of this emission might come from individual compact sources such as those known to exist in our own galaxy .
however , the subsequent discovery of iron emission lines in x - ray spectra ( forman & jones 1982 , and references therein ) clearly demonstrated that thermal bremsstrahlung from a diffuse intracluster gas is the dominant source of x - ray emission in clusters , and interest in compact sources quickly waned .
nevertheless , compact sources must surely contribute some fraction of the total x - ray emission from galaxy clusters .
a growing body of evidence indicates that a substantial intergalactic population of stars , star clusters and stellar remnants is present in the cores of rich clusters , and some of these objects will be x - ray sources .
for example , ferguson , tanvir , & von hippel ( 1998 ) showed the existence of intergalactic stars in the virgo cluster , and planetary nebulae associated with this population have also been detected ( e.g. , arnaboldi et al .
1996 ; ciardullo et al . 1998 ) .
similarly , west et al . ( 1995 ) hypothesized the existence of a large population of intergalactic globular clusters in the cores of galaxy clusters based on observed galaxy - to - galaxy variations in globular cluster populations .
furthermore , there is mounting evidence that many cluster galaxies have been destroyed over a hubble time , victims of collisions with other galaxies or disruption by the mean tidal field of the cluster potential in which they reside ( e.g. , thompson & gregory 1993 ; lopez - cruz et al .
1997 ; gregg & west 1998 ; moore et al .
1999 ; calcaneo - rodin et al . 2000 ) .
as galaxies are disrupted they will spill their contents throughout intracluster space , including x - ray sources such as x - ray binaries , globular clusters , and black holes .
such disruption is likely to be most frequent in the dense cluster core , and hence would be expected to produce a substantial pool of intracluster material , including compact x - ray sources , at the cluster center .
any significant contribution by these compact sources to the x - ray flux will have to be taken into account when estimating the properties of icm , including temperatures , densities , and cooling flow mass deposition rates .
an intracluster population of compact objects would have a different spectral signature than thermal emission from intracluster gas .
for instance , the spectra of large spiral and elliptical galaxies reveal the presence of hard x - ray emission which is thought to be produced by either high mass x - ray binaries or accretion onto massive black holes and is adequately represented by power - laws having @xmath9 @xmath0 0.6 - 3.0 ( where @xmath10 ; colbert & mushotzky 1999 ; allen , di matteo , & fabian 2000 ) .
luminous agn , on the other hand , are often fit by power - laws with @xmath9 @xmath0 1.6 - 2.0 ( nandra et al .
indeed , the spectral signatures of such objects have already been identified in several x - ray clusters by allen et al .
( 2001 ) . using asca data
, these authors found very luminous sources with power - law spectral indices characteristic of compact sources .
the recent discovery of luminous _ extended _ non - thermal x - ray emission in several weak ( or non- ) cooling flow clusters ( rephaeli , gruber , & blanco 1999 ; fusco - femiano et al .
1999 ; 2000 ; 2001 ) suggests that one should also look for such emission in massive cooling flow clusters . at present , the origin of the diffuse non - thermal emission is still a matter of some debate .
suggested mechanisms include inverse - compton scattering of cmb photons up to x - ray energies via a population of highly relativistic electrons , and non - thermal bremsstrahlung , likewise , by a population of relativistic electrons ( sarazin 2001 ; rephaeli 2001 ) .
those highly relativistic electrons may have been accelerated by recent subcluster merger events which shocked the icm or may have been ejected from a central massive black hole .
the discovery of x - ray holes at the positions of radio lobes of a central compact object in several massive cooling flow clusters ( e.g. bhringer et al .
1993 , mcnamara et al .
2000 ; fabian et al .
2000 ) certainly lends credence to the second hypothesis .
recently , harris et al . ( 2000 ) ( chandra data ) and bhringer et al . ( 2001b ) ( xmm - newton data ) have even identified extended non - thermal emission directly associated with radio jets / lobes in 3c295 and m87 , respectively .
the spectral signature of extended non - thermal emission is also different from that of thermal emission . to date , the extended non - thermal components of x - ray spectra of clusters have adequately been modeled by power - laws having @xmath11 ( fusco - femiano et al . 1999 ; 2000 ; 2001 ) . as with luminous compact sources ,
a proper deduction of cluster properties should explicitly include the effects of extended non - thermal emission in the analysis .    in the present work we assess the case for a power - law component to the x - ray spectra of cluster centers , in particular those with cooling flows .
we show that simultaneous fits of rosat pspcb and asca gis / sis spectra , which have been used to infer the presence of cooling flows in many clusters to date , are in fact consistent with a significant contribution from compact sources and/or extended non - thermal emission . using rosat hri data , we have placed limits on the contribution of any bright point sources within the central regions of these clusters and conclude that the power - law components are most likely produced by either a collection of many low - luminosity compact sources ( with @xmath12 ergs s@xmath3 ) or by diffuse non - thermal emission or by some combination of the two .
we shall assume @xmath13 = 50 km s@xmath3 mpc@xmath3 and @xmath14 throughout .
ten well - studied x - ray luminous galaxy clusters ( table 1 ) were chosen to span a wide range of cooling flow mass deposition rates .
five of these have essentially no cooling flow ( ncf ) , while the other five cooling flow ( cf ) clusters have calculated mass deposition rates ranging from 189 to 645 solar masses per year .
observations by the rosat position sensitive proportional counter ( pspcb ) and high resolution imager ( hri ) and the asca gas imaging spectrometer ( gis ) and solid - state imaging spectrometer ( sis ) were obtained from the heasarc public archives and reduced using ftools v5.01 .
the combination of rosat and asca data provides good sensitivity to both soft and hard x - ray emission , spanning the total energy range 0.14 - 10.0 kev , and should provide a strong test of power - law components in the x - ray spectra of clusters .
the rosat pspcb data were corrected for spatial and temporal gain fluctuations in the detector .
x - ray spectra , spanning the energy range 0.14 - 2.04 kev , of all ten clusters were extracted from circular regions centered on the peak of the emission with a constant extraction radius ( r ) of 3 arcminutes .
background spectra were selected from source - free circular regions of radius r. all spectra were corrected for the effects of cosmic rays in the detector ( snowden et al .
1992 ) , vignetting , spatial variations in detector efficiency , and deadtime .
the background - subtracted spectra were then grouped to ensure a minimum of 20 counts per radial bin so that reliable @xmath15 statistics could be computed .
the rosat hri data was reduced in a manner similar to the pspcb data .
the hri data was used to place limits on the contribution of any point sources within the central regions of these clusters .
unfortunately , the hri lacks spectral resolution and was not used in the main spectral analysis .
a description of the analysis of hri data is given in section 5 .
no hri data were available for a2063 or a3532 .      the screened events from `` revision 2 '' processing of asca gis / sis data were used .
additional cleaning procedures ( discussed in the asca abc guide ) were implemented , including gain and deadtime corrections . for the gis detectors ( gis2 and gis3 ) , x - ray spectra spanning the energy range 0.7 - 10.0 kev were extracted from the same region of sky as that of the rosat pspcb source spectra . for the sis detectors ( sis0 and sis1 ) , x - ray spectra spanning the range 2.0 - 10.0 kev
were extracted , thus avoiding complications with the degrading efficiency of the sis detectors at low energies .
background spectra for both the sis and gis detectors were extracted from the `` blank sky '' observations made during the performance testing stage of the mission .
spectra were extracted from the same regions of the detectors as that of the source spectra ( thus minimizing position - dependent effects associated with background subtraction ) .
the background - subtracted spectra were also grouped ( or binned ) to ensure @xmath15 statistics could be used . for the clusters a2142 and a3667 sis1 data were ignored because of an excess of hot and flickering pixels .
reduced rosat pspcb and asca gis / sis spectra were fit simultaneously using models within the software package xspec v11 ( arnaud 1996 ) .
each cluster was fit with a variety of models based on mekal thermal plasma routine ( mewe & gronenschild 1985 ; kaastra 1992 ; liedahl et al .
1995 ) , and/or redshifted power - laws .
in addition , we used the popular cooling flow model comprising a thermal plasma ( mekal ) and cooling flow component implemented by the mkcflow routine within xspec ( mushotzky & szymkowiak 1988 ) .
fits were obtained by varying the temperature ( @xmath16 ) and metallicity ( @xmath17 ) of thermal plasma component(s ) and by varying the photon index ( @xmath9 ) of any power - law components .
photoelectric absorption ( balucinska - church & mccammon 1992 ) was applied to the spectral models through two separate methods : i ) by applying only one absorption component to the entire spectral model and fixing the hi column density ( @xmath18 ) at the known galactic value ( see table 1 ) and then letting it vary to see if the fit could be improved ; ii ) by applying one absorption component to the entire spectral model ( with @xmath18 fixed at the galactic value ) and a separate redshifted variable absorption component ( @xmath19 ) to any power - law components via the zphabs model in xspec .
the second method allows for _ intrinsic _ absorption to be associated with any power - law components .
spectra of nearby high mass x - ray binaries and accreting massive black holes provide strong evidence for heavy intrinsic absorption , with @xmath20 @xmath21 above the galactic column ( colbert & mushotzky 1999 ; iwasawa et al . 2000 ; akylas et al . 2001 ) .
we therefore regard the addition of a free intrinsic absorption parameter as an important feature of our spectral analysis .
henceforth we will use tp , pl , and cf to refer , respectively , to thermal plasma ( mekal ) , power - law , and mkcflow components .
for example , in this notation tppl identifies a model having thermal plasma and power - law components .
goodness of fit was judged by the reduced chi - square ( @xmath22 ) and by visual inspection of residual plots .
when additional components were added to a spectral model , for example adding a power - law component to a thermal plasma component to make tppl , the corresponding f - statistic was calculated in order to judge whether the addition was reasonable . unless stated otherwise , the derived relative contributions of each of the components to the total flux refer to the energy range 0.14 - 10.0 kev stated , whereas the luminosities presented refer to the asca gis band only ( 0.7 - 10.0 kev ) , since the spectra were fit simultaneously and not combined .
all luminosities presented were corrected for galactic absorption only .    to check our data reduction procedures , spectra from each of the detectors were also fit individually and no significant deviations were found .
we have also attempted to duplicate published results by arabadjis & bregman ( 2000 ) for three clusters derived using the same data ( rosat pspcb ) and the same software .
the results , which refer to the inner 3 arcminute radius , are presented in table 2 . here , the hotter tp component is fixed at the white , jones , & forman ( 1997 ) value , @xmath23 , @xmath18 is fixed at the stark et al .
( 1992 ) value , and the energy range is 0.14 - 2.4 kev and 0.14 - 2.04 kev for ab2000 and our study , respectively .
the comparison shows excellent agreement , with derived temperatures and mass deposition rates having the same values to within one standard deviation .
we find that single - temperature tp models provide acceptable fits to the spectra of a119 , a2063 , a3266 , a3532 , and a3667 ( 1.00 @xmath24 @xmath22 @xmath24 1.17 ) with temperatures in very good agreement with previous studies ( white , jones , & forman 1997 ; markevitch et al .
significant improvement is not usually obtained by varying the hydrogen column density from the known galactic value reported by dickey & lockman ( 1990 ) or stark et al .
the derived metallicities range from @xmath25 which is typical of the intracluster medium ( allen & fabian 1997 ) .    in general , both pure pl models with or without variable intrinsic absorption
do not provide acceptable fits to the ncf cluster spectra with @xmath26 @xmath27 1.4 .
this result is not unexpected since it leaves no room for thermal emission from the icm . with intrinsic absorption
included , an acceptable fit is obtained for a3532 with @xmath28 , however .
this cluster is located at the center of the shapley concentration and is thought to be undergoing a merger with a3528 .
the diffuse non - thermal emission generated from such a merger might possibly account for the ability of the simple pl model to fit the spectrum of a3532 .
unfortunately , no hri data exists for this cluster so we can not place limits on the contribution of any point sources within the central 3 arcminutes .    0.1 in    0.1 in    in the best fitting tppl models , the pl component typically contributes 1 - 14% ( mean of 7 % ) when absorbed by the galactic column only or 10 - 25 % ( mean of 20 % ) when intrinsic absorption is included ( excluding a3532 which has an anomalously large pl component ) .
we note that in virtually every case the inclusion of intrinsic absorption noticeably improves the quality of the fit , making it superior to even the simple tp models ( see table 3 ) .
a mean f - statistic of @xmath29 ( @xmath27 95 % confidence ) was found when the pl component with intrinsic absorption was included in the fit indicating that its addition is , indeed , reasonable .
the corresponding luminosities of intrinsically absorbed pl components range from @xmath30 erg s@xmath3 with a mean value of @xmath31 erg s@xmath3 ( corrected for galactic absorption only ) .
photon indices for these components are typical of compact accretion - driven x - ray sources and extended non - thermal emission .
interestingly enough , the amount intrinsic absorption ranges from @xmath32 @xmath21 which is typical of optically thick compact sources .
we also note that the temperature of the icm tends to decrease significantly when a power - law component is included .
this trend was previously reported by allen et al .
( 2001 ) when including a power - law component in their spectral models .
figure 1 presents the best fit tppl model to the spectrum of a3667 . for clarity ,
only the rosat pspcb data are plotted .
as expected , standard cooling flow ( tpcf ) models of our non - cooling flow sample yield best fits with @xmath33 = 0 m@xmath7 yr@xmath3 to within one sigma .
the mass deposition rate for the modest cooling flow cluster a2063 , however , is 25 @xmath34 12 m@xmath7 yr@xmath3 .
the cooler component in our tptp models typically contributes zero flux to within one sigma ( @xmath0 10 % for a2063 & a3667 ) .
neither of these models provided fits significantly better than those of the standard tp models .    in summary
, we confirm earlier investigations which find that tp models provide satisfactory fits to rosat and asca spectra of the centers of non - cooling flow clusters .
however , tppl models which include intrinsic absorption generally provide fits of superior quality , and predict that the power - law component contributes 20 percent of the total flux . the photon indices of this component are typical of accreting x - ray sources and/or diffuse non - thermal x - ray emission which usually range from @xmath35 .
we have fit the spectra of the cf clusters with the same variety of models used in the previous section , but also with an additional model ; tptppl .
this model represents the emission from the observed region as the sum of the hot icm , a diminished cooling flow ( second tp component ) , and intermixed compact sources and/or extended non - thermal emission characterized by the pl component .
simple tp models generally provide fairly reasonable fits to the spectra with 1.1 @xmath36 @xmath22 @xmath36 1.25 .
derived temperatures and metallicities are in good agreement with previous studies ( white , jones , & forman 1997 , allen & fabian 1997 ; sarazin & mcnamara 1997 ; markevitch et al .
1998 ; allen 2000 ; ab2000 , allen et al . 2001 ) .
none of the fits were significantly improved by letting the value of @xmath18 vary .
hence , we find no evidence for excess absorption in these clusters , which is also consistent with several other investigations ( sarazin 1996 ; allen & fabian 1997 ; sarazin & mcnamara 1997 ; allen 2000 ; ab2000 ) .    as anticipated , pure pl models ( with or without intrinsic absorption ) do not provide acceptable fits to the x - ray spectra of any of the cf clusters . typically , @xmath26 @xmath27 1.5 for this model .
our tppl models fit the spectra of all 5 cf clusters very well , with 1.0 @xmath36 @xmath22 @xmath36 1.2 ( see table 4 ) . with the exception of a2029 , which is best fit by simple tp models
, the pl component typically contributes 5 - 30 % ( mean of 20 % ) of the total flux when absorbed only by the galactic column .
the inclusion of intrinsic absorption significantly improves the quality of the fit and raises the contribution of the pl component to 9 - 47 % ( mean of 27 % ) .
the mean f - statistic value for the addition of this power - law component is @xmath37 which suggests it is significant ( @xmath38 99 % confidence ) .
the luminosity ( 0.7 - 10.0 kev ) of pl components in these four clusters span @xmath39 ergs s@xmath3 , with a mean value of @xmath40 ergs s@xmath3 .
the photon indices of pl components are reasonably well constrained and range from 1.0 @xmath41 2.4 ( mean of 1.8 ) ( see table 4 ) .
the implied intrinsic column densities acting on the pl component are rather high for three of the cf clusters ( a644 , a1689 , & a2597 ) and are consistent with high - mass compact sources . as with the ncf clusters ,
a lower temperature is found for the thermal component than is found from the simple tp models .
tptp models were also implemented .
the quality of the fits provided by these models is generally comparable to that of tppl models .
calculation of the f - statistic suggests that the addition of the extra tp component is significant .
thus , we confirm previous reports that multi - phase thermal plasmas fit the spectra of cf clusters much better than that of isothermal models . the temperatures and metallicites of the two thermal components are consistent with the results of previous studies ( white , jones , & forman 1997 , markevitch et al .
1998 ; allen 2000 ; ab2000 ; allen et al .
2001 ) . for the cooler component in particular , we find @xmath42 and a mean contribution of roughly 22 % of the total flux , which is typical of cooling flow clusters ( peres et al .
1998 ) .
the standard cooling flow model ( tpcf ) provides adequate fits to the spectra of the cf clusters , as expected .
derived mass deposition rates , temperatures , and metallicities are generally consistent with previous studies ( white , jones , & forman 1997 ; peres et al .
1998 ; allen et al . 1999
; allen 2000 ; ab2000 ) . the best fitting tpcf model to the spectrum of a644
is presented in figure 2 . only the asca sis0 and gis2 data are plotted . a comparison of the tpcf and tppl models is presented in table 4 .
it is clear that the tppl models fit the spectra of the 5 cf clusters just as well as standard tpcf models and without any additional complexity .
0.1 in    0.1 in    finally , the tptppl model also fits the x - ray spectra of these clusters well , with 1.0 @xmath36 @xmath22 @xmath36 1.2 .
when absorbed by the galactic column only , the pl component typically contributes 2 - 5 % of the total flux inside 3 arcminutes .
the addition of intrinsic absorption on the pl component , again , improves the fit and the mean contribution of the pl component is raised to @xmath43 32 % . when compared to the tppl model , the mean f - statistic value is about 5 ( @xmath43 68 % confidence ) which suggests the addition of a second tp component only improves the fit marginally . the luminosity ( 0.7 - 10.0 kev ) of the pl component ranges from @xmath44 ergs s@xmath3 with a mean value of @xmath45 ergs s@xmath3 .
the power - law components have photon indices of 1.8 @xmath41 3.3 with a mean of 2.3 . with the exception of a2029 ( which has virtually no contribution from the pl component ) , the intrinsic column densities on the pl component are @xmath46 @xmath21 .
the temperature of the cooler tp component did not change significantly from the value derived from tptp models .
however , it s contribution to the flux was generally lower than that of tptp models ( mean @xmath0 10 % ) .    in summary , we find that in all cases the tppl and tptppl models fit the rosat and asca spectra of the five cf clusters just as well as the standard cooling flow model ( tpcf ) or the two - temperature thermal plasma ( tptp ) .
the pl components of these models are quite luminous ( contributing between 27 - 32 % of the total flux ) and characterized by photon indices typical of compact accreting x - ray sources and/or extended non - thermal emission .
this , then , is the main conclusion of our work : _ rosat and asca spectra can be satisfactorily modeled by attributing a significant part of the luminosity within the central regions of cf clusters to sources other than cooling gas . _
our analysis shows that luminous ( l@xmath47 @xmath0 10@xmath48 ergs s@xmath3 ) power - law components with @xmath49 might be present in the centers of rich clusters . including intrinsic absorption dramatically improves the quality of the fit ; in these circumstances , there is an apparent dichotomy in the importance of the pl components between ncf and cf clusters .
not only is the typical pl component in the cf cluster more luminous by about a factor of 10 than that of its ncf counterpart , but it also contributes a larger fraction of the total x - ray flux within the central 3 arcminutes of its cluster ( between 27 - 32 % for cf clusters as opposed to 20 % for ncf clusters ) .
this may be consistent with the fact that the _ cores _ of cf clusters are more often observed to be active ( i.e. radio loud ) , perhaps implying the presence of a central massive black hole(s ) and a surrounding diffuse plasma of relativistic electrons , than the cores of ncf clusters ( burns 1990 ) .
it is also worth noting that the improvement of the fit to the spectra of cf clusters far exceeds the improvement of the fit to the spectra of ncf clusters when a power - law component is included in the spectral model .
the most attractive feature of the tppl models is that one does not expect large quantities of atomic or molecular gas to be present in the centers of clusters , nor does one expect to see emission lines below a few kev , because cooling is not significant .
this is consistent with current observations . on the other hand , evidence which apparently supports the existence of cooling flows ( e.g. excess absorption , detection of co line emission in extreme cf clusters ) may contradict this model .
however , the success of the single - temperature plasma + power - law model in fitting the spectra of the central regions of rich clusters , regardless of cooling flow status , has also been reported by a number of authors using a variety of x - ray instruments , including markevitch & vikhlinin ( 1997 ) ( asca and rosat data ) , markevitch et al .
( 1998 ) ( asca data ) , rephaeli , gruber & blanco ( 1999 ) ( rxte data ) , guainazzi & molendi ( 1999 ) ( beppo - sax data ) , fusco - femiano et al .
( 1999 ; 2000 ; 2001 ) ( beppo - sax data ) , and recently by sambruna et al .
( 2000 ) who used chandra to uncover a previously _ hidden _ agn in hydra a with @xmath50 ergs s@xmath3 and @xmath51 . however , our study is the first that systematically addresses the problems of the standard cooling flow model by implementing the single - temperature plasma + power - law model .
the tptppl models , which might represent a compromise between the standard cooling flow model and our favoured tppl models , also provide good fits to the spectra of the cf clusters .
however , the addition of the second , cooler tp component only marginally improves the fit and at the expense of adding much more complexity ( and uncertainty ) to the model .
furthermore , the predicted temperatures and relative contributions of the cooler tp component in the tptppl models suggest that emission lines below a few kev should be present , however they have failed to turn up in chandra and xmm spectra .
finally , because the mean contribution of the cooler component is @xmath0 10 % , as compared with 20 % in the tptp models , this implies that mass deposition rates are only roughly reduced by a factor of two . using models similar to our tptppl models , allen et al .
( 2001 ) reached the same conclusion for a number of clusters in their sample .
these reduced deposition rates are still far too large to be consistent with the small amounts of atomic and molecular gas found in cooling flow clusters .
we expect that chandra and xmm - newton will place much tighter constraints on the parameters of the complex tptppl models and it will be interesting to see if the models can be reconciled with problems of the standard cooling flow model .
assuming that a power - law component is , indeed , present in the spectra of these clusters , some of the possible candidates are :    * seyfert - like agn could contribute most , if not all , of the pl component luminosity found in our study .
these objects can be quite luminous and are often observed to have @xmath52 @xmath53 ergs s@xmath3 with power - law indices of 1.6 - 2.0 ( nandra et al . 1997 ) .
in fact , allen et al .
2001 have suggested that agn are the likely candidate for the observed pl component in their sample .
likewise , the luminosities and steepness of the pl components in our sample are generally consistent with emission from a central agn(s ) . * diffuse non - thermal emission , which has possibly been seen in the central regions of coma , a2256 , and a2199
. searches of those clusters using beppo - sax and rxte have yielded very luminous ( @xmath43 10@xmath48 ergs s@xmath3 ) sources which are described by power - laws with @xmath54 and may possibly be linked to power - law components observed in the asca band ( rephaeli , gruber , & blanco 1999 ; fusco - femiano et al .
1999 ; 2000 ; 2001 ) .
unfortunately , x - ray spectra usually have to be extracted from regions _ much _ larger than that of the central galaxy and , thus , it is difficult to draw conclusions about the source of the power - law component .
a vla survey of more than 200 nearby clusters found that only 29 had extended radio emission on mpc scales ( giovannini , tordi , & feretti 1999 ) perhaps implying that extended non - thermal emission is not that common in galaxy clusters .
it should be noted , however , that the detection of true intracluster radio emission is exceptionally difficult , as one must remove bright radio galaxies from the analysis ( rephaeli 2001 ) .
furthermore , the non - thermal emission may be associated with the radio galaxies themselves ( e.g. in the radio lobes ) . * massive black holes associated with elliptical and lenticular galaxies .
the regions sampled in our study generally include a single giant elliptical or cd galaxy plus a number of large galaxies , probably e or s0 galaxies .
colbert & mushotzky ( 1999 ) and allen , di matteo , & fabian ( 2000 ) have recently discovered hard power - law components in several nearby ge galaxies with @xmath55 ergs s@xmath3 with photon indices of 0.5 - 3.0 .
such sources are typically less luminous than seyfert nuclei , perhaps as a result of low radiative efficiency from an advection - dominated accretion flow ( adaf ; narayan & yi ) . while several of these sources could , in principle , account for the pl component in the ncf clusters they would have to be scaled up by a factor of 100 to account for the emission in the cf clusters .
faint galaxies , on the other hand , are much more numerous and many of them might harbor less massive central black holes . the local group elliptical m32 , whose central object has @xmath56 erg s@xmath3 ( loewenstein et al .
1998 ) offers a possible guide to the x - ray luminosity of faint cluster members .
hence @xmath57 such galaxies , each with a supermassive black hole , would be required to account for the power - law component of the clusters studied here . on this basis ,
faint galaxies seem incapable of producing significant power - law emission inside the regions we have observed .
we note , however , that if massive black holes originally belonging to bright or faint galaxies are now scattered throughout intracluster space ( i.e. no longer associated with optically visible galaxies ) as the result of past mergers or stripping both of these limits may be relaxed . *
galactic x - ray binaries , having @xmath58 ( white et al .
1988 , verbunt et al .
1995 ) , and which in our galaxy are most commonly associated with globular clusters .
intergalactic globular clusters , if similar to globular clusters in the milky way , would have to be present in large numbers to collectively account for @xmath59 erg s@xmath3 attributable to the pl component in our models . assuming a typical intergalactic globular contains a binary x - ray source with @xmath60 @xmath61 ergs s@xmath3 ( hertz & wood 1985 ) and the specific frequency of globulars in the central ge galaxy is @xmath62 ( harris 1991 ) then even a very bright central galaxy with m@xmath63 = -24 would have a total x - ray luminosity of @xmath64 ergs s@xmath3 , too low by a factor of @xmath65 .
even the most optimistic models of intracluster globular clusters of west et al .
( 1995 ) predict @xmath66 such clusters , significantly less than is required to account for the entire pl component .
hence , it appears that there are probably two or three sources capable of producing the luminosity of the pl components of our models : massive black holes associated ( or once associated ) with large elliptical galaxies , diffuse non - thermal emission , and agn .
does existing data favor one type of source over the others ?    to determine the luminosity of any point sources within the central 3 arcminutes we used rosat hri data .
two different procedures were implemented : i ) spectra of clearly visible point sources were extracted from a 5  region centered on the source and grouped into one large energy bin ( 0.5 - 2.0 kev ) .
the spectra were then fit within the xspec package using a simple power - law model with @xmath67 and absorbed by the galactic column only . the 0.5 - 2.0 kev flux was then converted into a 0.7 - 10.0 kev unabsorbed flux using the program w3pimms .
if more than one point source was present the individual fluxes were summed .
as an example , figure 3 is the hri observation of a3667 , smoothed with a 5  gaussian beam .
there are at least two bright compact sources within the central 3 arcminutes of this cluster . ;
ii ) to place limits on the contribution of point sources which are hidden in the bright diffuse background ( what we refer to as `` barely visible '' point sources ) we calculated the mean standard deviation of counts per pixel ( within 3 ) and then scaled this to a standard deviation per circular region of diameter = 5  ( comparable to the hri psf ) .
the latter value was expressed as unabsorbed flux ( 0.7 - 10.0 kev ) using the program w3pimms and assuming a simple power - law spectral model with @xmath67 absorbed only by the galactic column listed in table 1 .
we regard the result as a 1-sigma limit to the contribution of any point sources .
0.1 in    0.1 in    the results of hri data analysis are presented in table 5 .
one or more point sources are easily visible within the central 3 arcminutes of all three ncf clusters having hri observations .
it is clear from a comparison with the luminosities of the observed pl component from the best fitting tppl models ( with intrinsic absorption ) , however , that the bright point sources can not be solely responsible for the power - law emission , as they contribute only @xmath43 5 % of the pl component luminosity .
significant contributions are required from low - luminosity point sources and/or diffuse non - thermal emission .
the third column in table 5 ( method ii ) presents the 3-sigma limit of `` barely visible '' point sources , suggesting that at least @xmath68 of them are needed ( in addition to the visible point sources ) to account for the entire power - law component .
there are no visible point sources within the central 3 of the 5 cf clusters studied .
this is almost certainly the result of the cf clusters being much more centrally peaked than the ncf clusters .
at least @xmath68 `` barely visible '' point sources are required to account for the observed pl component in these clusters .
perhaps another way to differentiate between compact sources and diffuse non - thermal emission is through the presence or absence of flux variability , since most luminous accreting x - ray sources are known to vary significantly on short time scales ( e.g. nandra et al . 1997 ) .
however , because the pl component typically contributes only @xmath69 % this might make detection of variability difficult in clusters ( especially in the cf clusters ) .
furthermore , if the pl component is produced by @xmath68 point sources , and the variations in each are random , it will make the variation in the collective luminosity quite small , perhaps rendering this method unusable .
if compact or diffuse non - thermal sources provide significant emission from the centers of rich clusters , it will be necessary to revise a number of conclusions derived from x - ray data . fully exploring the consequences
is beyond the scope of the present paper ; however , we will mention a few important issues .
allen ( 2000 ) reports that a deprojection analysis of rosat hri data yields cooling flow mass deposition rates that agree reasonably well with those derived from asca spectra .
a consistent re - analysis would require first subtracting the power - law contribution from the x - ray images , but of course the spatial distribution of the power - law emission must first be understood .
the presence of emission lines from heavy ions originally led astronomers to accept thermal emission as the dominant mechanism in rich cluster icms .
accretion - driven x - ray sources do exhibit line emission ( kallman , vrtilek , & kahn 1989 ; narayan & raymond 1999 ) , but the lines are produced under quite different physical conditions than those within the canonical hot icm .
the implications for icm abundance estimates and enrichment scenarios will require careful study .
likewise , previously derived temperatures of the icm will also be affected ( as demonstrated in table 3 ) .
might emission from compact sources account for the excess hard x - radiation recently discovered ( xu et al .
1998 ; ikebe et al .
1997 ) in the centers of several clusters ?
the next generation of x - ray detectors and satellites ( i.e. xmm and chandra ) will be needed to map the x - ray emission with adequate spatial resolution to address this and other questions raised by the present work .
we have shown that rosat pspcb and asca gis2/gis3 spectra of ten x - ray luminous clusters permit a significant contribution from compact sources and/or diffuse non - thermal emission .
although this conclusion holds irrespective of the presence of a cooling flow , we find that power - law components are especially luminous in clusters with massive cooling flows .
our study indicates that inferred cooling flow mass deposition rates may be lowered substantially , thus ameliorating , or possibly resolving , the long - standing problem of the missing cool gas .
existing x - ray observations of large elliptical galaxies , agn , and extended non - thermal emission suggest that any of these sources could account for the observed luminosity of the power - law component . however , the absence of extremely bright point sources in the rosat hri observations rules out the possibility that the power - law emission comes from a few extremely bright compact sources .
more likely is that most of the power - law emission is produced by some combination of diffuse and low - luminosity compact sources , such as massive black holes .
the deep imaging capabilities and superior spectral resolution of chandra and xmm should be able to place much tighter constraints on the relative contributions of power - law sources and thermal bremsstrahlung to cluster x - ray emission and , perhaps , even determine the source of the power - law emission .
we thank the rosat and asca help teams , in particular jane turner , for advice with data reduction techniques .
we also thank paul nandra , peter thomas , arif babul , ann gower , pat henry , and alastair edge for useful comments and suggestions .
igm acknowledges support from nserc of canada ( usra & pgsa ) .
gaw and mjw were supported by nserc research grants .
mjw also acknowledges support from nsf grant ast-0071149 .
akylas , a. , georgantopoulos , i. , & comastri , a. , 2001 , mnras , 324 , 521 allen , s.w . , 2000 ,
mnras , 315 , 269 allen , s.w . , di matteo , t. , & fabian , a.c . 2000 , mnras , 311 , 493 allen , s.w . , et al . , 1993 , mnras , 262 , 901 allen , s.w . , et al . , 2001 , mnras , 322 , 589 allen , s.w . , & fabian , a.c . , 1997 , mnras , 286 , 583 arabadjis , j.s . , & bregman , j.n . , 2000 , apj , 536 , 144 ( ab2000 ) arnaboldi , m. , et al . , 1996 , apj , 472 , 145 arnaud , k.a . , 1996 , in asp conf .
ser . 101 , astronomical data analysis software and systems v , ed .
jacoby , g. , & barnes , j. , ( san francisco : asp ) , 17 balucinska - church , m. , & mccammon , d. , 1992 , apj , 400 , 699 bardelli , s. , zucca , e. , & baldi , a. , 2001 , mnras , 320 , 387 bauer , f. , & sarazin , c.l . , 2000 , apj , 530 , 222 bazzano , a. , et al .
, 1984 , apj , 279 , 515 bhringer , h. , et al . , 1993 , mnras , 264 , l25 bhringer , h. , et al .
, 2001a , a&a , in press ( astro - ph/0111112 ) bhringer , h. , et al . , 2001b , a&a , 365 , l181 burns , j.0 . , 1990 , aj , 99 , 14 calcneo
roldn , c. , et al . , 2000 , mnras , 314 , 324 ciardullo , r. , et al . , 1998 ,
aj , 118 , 488 colbert , e.j.m . , & mushotzky , r.f . , 1999 , apj , 519 , 89 dickey , j.m & lockman , f.j .
1990 , ara&a , 28 , 215 donahue , et al . , 2000 , apj , 545 , 670 edge , a.c . , 2001 , mnras , in press ( astro - ph/0106225 ) fabian , a.c . , 1994 , ara&a , 32 , 277 fabian , a.c . , et al
mnras , 318 , 65 fabian , a.c . , mushotzky , r.f . , nulsen , p.e.j . , & peterson , j.r . , 2001 ,
mnras , 321 , l33 feretti , l. , et al . , 1999 ,
a&a , 344 , 472 ferguson , h.c . , tanvir , n.r . ,
& von hippel , t. , 1998 , nature , 391 , 461 flores , r.a . , quintana , h. , & way , m.j .
, 2000 , apj , 532 , 206 forman , w. , & jones , c. , 1982 , ara&a , 20 , 547 fusco - femiano , r. , et al .
, 1999 , apj , 513 , l21 fusco - femiano , r. , et al . , 2000 , apj , 534 , l7 fusco - femiano , r. , et al . , 2001 ,
apj , 552 , l97 gregg , m.d . , & west , m.j . , 1998 , nature , 396 , 549 giovannini g. , tordi , m. , & feretti , l. , 1999 , new ast . , 4 , 141 guainazzi , m. , & molendi , s. , 1999 , a&a , 351 , l19 harris , w. e. , 1991 , ara&a , 29 , 543 harris , d.e .
et al . , 2000 ,
apj , 530 , l81 heckman , t.m . ,
baum , s.a . ,
van breugel , w.j.m . , & mccarthy , p.j .
, 1989 , apj , 338 , 48 hertz , p. , & wood , k.s . , 1985 , apj , 290 , 171 ikebe , y. , et al . , 1997 , apj , 481 , 660 irwin , j.a .
, & sarazin , c.l . , 1995 ,
apj , 455 , 497 iwasawa , k. , et al . , 2000 ,
mnras , 318 , 879 jaffe , w. , 1990 , a&a , 240 , 254 jaffe , w. , & bremer , m.n .
, 1997 , mnras , 284 , l1 kaastra , j.s . , 1992 , an x - ray spectral code for optically thin plasmas ( internal sron - leiden report , updated version 2.0 ) kallman , t.r . , vrtilek , s.d . , & kahn , s.m .
, 1989 , apj , 345 , 498    katz , j.i . , 1976 ,
apj , 207 , 25 liedahl , d.a . ,
osterheld , a.l . , & goldstein , w.h . , 1995 , apjl , 438 , 115 loewenstein , m. , hayashida , k. , toneri , t. , & davis , d.s .
, 1998 , apj , 497 , 681 lopez - cruz , o. , et al .
, 1997 , apj , 475 , l97 markevitch , m. , forman , w.r . ,
sarazin , c.l . , & vikhlinin , a. , 1998 , apj , 503 , 77 markevitch , m. , et al . , 2000 , apj , 541 , 542 markevitch , m. & vikhlinin , a. , 1997 , apj , 491 , 467    mcnamara , b.r . , bregman , j.n . ,
& oconnell , r.w . , 1990 , apj , 360 , 20 mcnamara , b.r . , et al . , 2000 , apj , 534 , l135 mewe , r. , gronenschild , e.h.b.m . ,
& van den oord , g.h.j . , 1985 ,
a&as , 62 , 197 molendi , s. , & de grandi , s. 1999 , a&a , 351 , l41 molendi , s. , & pizzolato , f. , 2001 , apj , in press ( astro - ph/0106552 ) moore , b. , lake , g. , quinn , t. , & stadel , j. , 1999 , mnras , 304 , 465 mushotzky , r.f . , &
szymkowiak , a.e . , 1988 , in cooling flows in clusters and galaxies , ed .
fabian , ( dordrecht : kluwer ) , 53 nandra , k. , et al .
, apj , 477 , 602 narayan , r. , & raymond , j. , 1999 , apj , 515 , l69 narayan , r. , & yi , i. , 1995 , apj , 452 , 710 odea , c.p . , et al .
, 1994 , apj , 422 , 467 oegerle , w.r . , & hill , j.m . , 1994 ,
aj , 107 , 857 peres , c.b . , et al . , 1998 , mnras , 298 , 416 peterson , j.r . , et al . , 2001 , a&a , 365 , l104 rephaeli , y. 2001 , astro - ph/0101363 rephaeli , y. , gruber , d.e . , & blanco , p. , 1999 , 511
, l21 sarazin , c.l . ,
1996 , in galactic and cluster cooling flows , ed .
n. soker , ( san francisco : publ .
pacific ) , 115 , 172 sarazin , c.l .
, 2001 , in merging processes in clusters of galaxies , ed
. l. feretti , m. gioia , and g. giovannini ( dordrecht : kluwer ) , in press ( astro - ph/0105418 ) sambruna , r.m . , et al .
, 2000 , apj , 532 , 91 .
sarazin , c. l. , burns , j.o .
, roettiger , k. , & mcnamara , b.r . , 1995 , apj , 447 , 559 .
sarazin , c.l . , & mcnamara , b.r . , 1997 ,
apj , 480 , 203 snowden , s.l . , et al . , 1992 , apj , 393 , 819 stark , a.a . ,
et al . , 1992 , apjs , 79 , 77 tamura , t. , et al . , 2001 , a&a , 365 , l86 thompson , l.a . , &
gregory , s.a .
, 1993 , aj , 106 2197 verbunt , f. , bunk , w. , hasinger , g. , johnston , h.m . , 1995 , a&a , 300 , 732 voit , g.m . &
donahue , m. , 1995 , apj , 452 , 164 west , m. , ct , p. , jones , c. , forman , w. , & marzke , r. , 1995 , apj , 453 , l77 white , d.a . , et al .
, 1991 , mnras , 252 , 72 white , d.a . ,
jones , c. , & forman , w. , 1997 , mnras , 292 , 419 white , n.e . , stella , l. , & parmar , a.n . , 1988 , apj , 324 , 363 xu , h. , et al . , 1998 , apj , 500 , 738    ccccccll a119 & 0.0442 & 3.17 & ncf & 0 & 5.3 & 3 bright radio galaxies near center@xmath70 + a644 & 0.0704 & 8.00@xmath71 & cf & 189 & 6.5 & minor merger event@xmath72
+ a1689 & 0.1832 & 1.81 & cf & 645 & 8.7 & + a2029 & 0.0773 & 3.14 & cf & 556 & 7.4 & non - thermal emission not detected@xmath73 + a2063 & 0.0353 & 2.90@xmath71 & ncf & 37 & 4.1 & subcluster merger?@xmath74 + a2142 & 0.0909 & 4.16 & cf & 350 & 11.4 & merger@xmath75 & hard x - ray emission detected@xmath76 + a2597 & 0.0852 & 2.49 & cf & 271 & 9.1 & powerful cd radio galaxy@xmath77 + a3266 & 0.0589 & 3.00@xmath71 & ncf & 0 & 6.2 & major merger@xmath78 + a3532 & 0.0554 & 5.98 & ncf & 0 & 4.7 & merger@xmath79 + a3667 & 0.0556 & 4.00@xmath71 & ncf & 0 & 7.1 & hard x - ray emission found@xmath80 +      a119 & tp & - & - & 0.987 & & - & - & 1.011 + & tpcf & 0.08 @xmath34 @xmath81 & 0 @xmath34 3 & 0.994 & & 0.01 @xmath34 151 & 0.6 @xmath34 5 & 1.026 + a2029 & tptp & 1.29 @xmath34 0.15 & - & 1.592 & & 1.19 @xmath34 0.31 & - & 1.499 + & tpcf & 0.08 @xmath34 6.97 & 132 @xmath34 31 & 1.911 & & 0.03 @xmath34 4.82 & 170 @xmath34 43 & 1.772 + a2142 & tpcf & 0.08 @xmath34 7.41 & 188 @xmath34 43 & 1.190 & & 0.10 @xmath34 9.16 & 230 @xmath34 73 & 1.130 +      a119 & tp & @xmath82 & & & & 467.4/446 + & tppl & @xmath83 & @xmath84 & 25 % & @xmath85 & 454.0/443 + a2063 & tp & @xmath86 & & & & 945.6/807 + & tppl & @xmath87 & @xmath88 & 8.5 % & @xmath89 & 933.3/804 + a3266 & tp & @xmath90 & & & & 1160.5/1117 + & tppl & @xmath91 & @xmath92 & 23 % & @xmath93 & 1158.0/1114 + a3532 & tp & @xmath94 & & & & 776.4/746 + & tppl & @xmath95 & @xmath96 & 87.5 % & @xmath97 & 756.0/743 + a3667 & tp & @xmath98 & & & & 809.0/791 + & tppl & @xmath99 & @xmath100 & 22 % & @xmath101 & 797.6/788 +      a644 & tppl & & @xmath102 & 35 % & @xmath103 & 1855.7/1613 + & tpcf & @xmath104 & & & & 1841.2/1613 + a1689 & tppl & & @xmath105 & 9 % & @xmath106 & 1643.6/1409 + & tpcf & @xmath107 & & & & 1653.9/1409 + a2029 & tppl & & @xmath108 & 5 % & @xmath109 & 1460.3/1149 + & tpcf & @xmath110 & & & & 1462.1/1149 + a2142 & tppl & & @xmath111 & 47 % & @xmath112 & 1010.4/868 + & tpcf & @xmath113 & & & & 1013.2/868 + a2597 & tppl & & @xmath114 & 17 % & @xmath115 & 1305.1/1281 + & tpcf & @xmath116 & & & & 1302.3/1281 +      a119 & @xmath117 & @xmath118 & @xmath119 + a644 & & @xmath120 & @xmath121 + a1689 & & @xmath122 & @xmath123 + a2029 & & @xmath124 & @xmath125 + a2142 & & @xmath126 & @xmath127 + a2597 & & @xmath128 & @xmath129 + a3266 & @xmath130 & @xmath131 & @xmath132 + a3667 & @xmath133 & @xmath134 & @xmath135 +