the large - order behavior of a perturbative expansion in gauge theories is inevitably dominated by the factorial growth of renormalon diagrams [ 1 - 4 ] . in the case of quantum chromodynamics ( qcd ) , the coefficients of perturbative expansions in the qcd coupling @xmath1 can increase dramatically even at low orders .
this fact , together with the apparent freedom in the choice of renormalization scheme and renormalization scales , limits the predictive power of perturbative calculations , even in applications involving large momentum transfers , where @xmath1 is effectively small .
a number of theoretical approaches have been developed to reorganize the perturbative expansions in an effort to improve the predictability of the perturbative qcd ( pqcd ) .
for example , optimized scale and scheme choices have been proposed , such as the method of effective charges ( ech ) [ 5 ] , the principle of minimal sensitivity ( pms ) [ 6 ] , and the brodsky - lepage - mackenize ( blm ) scale - setting prescription [ 7 ] and its generalizations [ 8 - 20 ] . in ref.[4 ] , the resummation of the formally divergent renormalon series and the parametrization of related higher - twist power - suppressed contributions has been given .
in general , a factorially divergent renormalon series arises when one integrates over the logarithmically running coupling @xmath2 in a loop diagram .
such contributions do not occur in conformally invariant theories which have a constant coupling .
of course , in the physical theory , the qcd coupling does run .    among the fundamental predictions of qcd are asymptotic scaling laws for large - angle exclusive processes [ 21 - 28 ] .
qcd counting rules were formalized in refs.[22,23 ] .
these reactions probe hadronic constituents at large relative momenta , or equivalently , the hadronic wave function at short distances .
in particular , the meson wave function determines the leading higher - twist contribution to meson production at high @xmath0 .
the hadronic wave functions in terms of quark and gluon degrees of freedom play an important role in the quantum chromodynamics predictions for hadronic processes . in the perturbative qcd theory ,
the hadronic distribution amplitudes and structure functions which enter exclusive and inclusive processes via the factorization theorems at high momentum transfers can be determined by the hadronic wave functions , and therefore they are the underlying links between hadronic phenomena in qcd at large ( nonperturbative ) and small distances ( perturbative ) .
if the hadronic wave functions were accurately known , then we could calculate the hadronic distribution amplitude and structure functions for exclusive and inclusive processes in qcd .
the study of hadron production at large transverse momenta in hadronic interactions is a valuable testing ground of the perturbative regime of qcd , providing information on both the parton distribution functions ( pdfs ) in the proton , and the parton - to - hadron fragmentation functions ( ffs)[29 ] .
the frozen coupling constant approach can be applied for investigation , not only for exclusive processes , but also for the calculation of higher - twist contributions to some inclusive processes , for example as large -@xmath3 meson photoproduction [ 30 ] , two - jet+meson production in the electron - positron annihilation [ 31 ] . in the frozen coupling constant approach
was used refs.[30,31 ] for calculation of integrals , such as @xmath4 according to ref.[7 ] , in pqcd calculations , the argument @xmath5 of the running coupling constant should be taken equal to the square of the momentum transfer of a hard gluon in a corresponding feynman diagram , in both the renormalization and factorization scale . but defined in this way , @xmath6 suffers from infrared singularities .
for example in ref.[32 ] , @xmath5 equals to @xmath7 and @xmath8 , where @xmath9 , @xmath10 are the subprocess s mandelstam invariants . therefore , in the soft regions
@xmath11 , @xmath12 , the integrals in ( 1.1 ) diverge and we need some regularization methods of @xmath13 in these regions for their calculation . in ref.[33 ] , the authors investigated the phenomenology of infrared renormalons in inclusive processes .
the dispersive approach has been devised to extend properly modified perturbation theory calculations towards the low - energy region [ 34 ] .
connections between power corrections for the three deep inelastic scattering sum rules have also been explored in ref.[35 ] .
investigation of the infrared renormalon effects in various inclusive and exclusive processes is one of the most important and interesting problems in the perturbative qcd .
it is known that infrared renormalons are responsible for factorial growth of coefficients in perturbative series for the physical quantities .
but , these divergent series can be resummed by means of the borel transformation [ 1 ] and the principal value prescription [ 36 ] , and effects of infrared renormalons can be taken into account by a scale - setting procedure @xmath14 at the one - loop order results .
technically , all - order resummation of infrared renormalons corresponds to the calculation of the one - loop feynman diagrams with the running coupling constant @xmath15 at the vertices or , alternatively , to calculation of the same diagrams with nonzero gluon mass .
studies of infrared renormalon problems have also opened new prospects for evaluation of power - suppressed corrections to processes characteristics [ 37 ] .
power corrections can also be obtained by means of the landau - pole free expression for the qcd coupling constant .
the most simple and elaborated variant of the dispersive approach , the shirkov and solovtsov analytic perturbation theory , was formulated in ref.[38 ] .    by taking these points into account
, it may be argued that the analysis of the next - to - leading order higher - twist effects on the dependence of the pion wave function in pion production at proton - proton collisions by the running coupling approach , are significant from both theoretical and experimental[39 ] points of view .    in this work
we apply the running coupling approach[40 ] in order to compute the effects of the infrared renormalons on the meson production in proton - proton collisions in next - to - leading order .
this approach was also employed previously[41,42 ] to calculate the inclusive meson production in proton - proton and photon - photon collisions . the running coupling method in next - to - leading order pion and kaon electromagnetic form factor was computed in [ 43 ] .    as we know ,
power - suppressed contributions to exclusive processes in qcd , which are commonly referred to as higher - twist corrections .
the higher - twist approximation describes the multiple scattering of a parton as power corrections to the leading - twist cross section .
we will show that higher - twist terms contribute substantially to the inclusive meson cross section at moderate transverse momenta .
in addition , we shall demonstrate that higher - twist reactions necessarily dominate in the kinematic limit where the transverse momentum approaches the phase - spase boundary .
a precise measurement of the inclusive charged pion production cross section at @xmath16 and @xmath17 is important for the proton - proton collisions program at the relativistic heavy ion collider ( rhic ) at the brookhaven national laboratory .
another important aspect of this study is the choice of the meson model wave functions . in this respect ,
the contribution of the higher - twist feynman diagrams to a pion production cross section in proton - proton collisions is been computed by using various pion wave functions .
also , higher - twist contributions which are calculated by the running coupling constant and frozen coupling constant approaches are been estimated and compared to each other . within this context
, this paper is organized as follows : in sec.[ht ] , we provide formulas for the calculation of the contribution of the high twist diagrams . in sec .
[ ir ] we present formulas and an analysis of the next - to - leading order higher - twist effects on the dependence of the pion wave function by the running coupling constant approach . in sec .
[ lt ] , we provide formulas for the calculation of the contribution of the leading - twist diagrams . in sec .
[ results ] , we give the numerical results for the cross section and discuss the dependence of the cross section on the pion wave functions .
we present our conclusions in sec .
the higher - twist feynman diagrams , which describe the subprocess @xmath18 for the pion production in the proton - proton collision are shown in fig.1 . in the higher - twist diagrams , the pion of a proton quark is directly observed .
their @xmath19 power suppression is caused by a hard gluon exchange between pion constituents .
the amplitude for this subprocess can be found by means of the brodsky - lepage formula [ 27 ] :    @xmath20    in eq.(2.1 ) , @xmath21 is the sum of the graphs contributing to the hard - scattering part of the subprocess .
the hard - scattering part for the subprocess under consideration is @xmath22 , where a quark and antiquark form a pseudoscalar , color - singlet state @xmath23 . here
@xmath24 is the pion wave function , i.e. , the probability amplitude for finding the valence @xmath25 fock state in the meson carry fractions @xmath26 and @xmath27 , @xmath28 . remarkably , this factorization is gauge invariant and only requires that the momentum transfers in @xmath21 be large compared to the intrinsic mass scales of qcd . since the distribution amplitude and the hard - scattering amplitude are defined without reference to the perturbation theory , the factorization is valid to leading order in @xmath29 , independent of the convergence of perturbative expansions .
the hard - scattering amplitude @xmath21 can be calculated in perturbation theory and represented as a series in the qcd running coupling constant @xmath30 .
the function @xmath31 is intrinsically nonperturbative , but its evolution can be calculated perturbatively . in our calculation
, we have neglected the pion and the proton masses . turning to extracting the contributions of the higher - twist subprocesses
, there are many kinds of leading - twist subprocesses in @xmath32 collisions as the background of the higher - twist subprocess @xmath33 , such as @xmath34 , @xmath35 , @xmath36 etc .
the contributions from these leading - twist subprocesses strongly depend on some phenomenological factors , for example , quark and gluon distribution functions in the proton and fragmentation functions of various constituents , etc .
most of these factors have not been well determined , neither theoretically nor experimentally .
thus they cause very large uncertainty in the computation of the cross section of process @xmath37 . in general ,
the magnitude of this uncertainty is much larger than the sum of all the higher - twist contributions , so it is very difficult to extract the higher - twist contributions .
the production of a hadron at large transverse momentum , @xmath0 , in a hadronic collisions is conventionally analyzed within the framework of perturbative qcd by convoluting the leading - twist @xmath38 hard subprocess cross sections with evolved structure and fragmentation functions .
the most important discriminant of the twist of a perturbative qcd subprocess in a hard hadronic collision is the scaling of the inclusive invariant cross section [ 22 - 25,44 ] ,    @xmath39    at fixed @xmath40 and center - of - mass angle @xmath41 .
as we know in the original parton model [ 45 ] the power fall - off is simply @xmath42 since the @xmath43 subprocess amplitude for point - like partons is scale invariant , and there is no dimensional parameter as in a conformal theory .
the detected hadron @xmath44 can be produced directly in the hard subprocess reaction as in an exclusive reaction , then such direct higher - twist processes can give a significant contribution since there is no suppression from jet fragmentation at large momentum fraction carried by the hadron , @xmath45 , and the trigger hadron is produced without any waste of energy[44 ] .
the mandelstam invariant variables for subprocesses @xmath18 are defined as @xmath46    in our calculation , we have also neglected the quark masses . we have aimed to calculate the pion production cross section and to fix the differences due to the use of various pion model functions . we have used five different wave functions : the asymptotic wave function ( asy ) , the chernyak - zhitnitsky wave function [ 28,46 ] , the cleo wave function [ 47 ] , the braun - filyanov pion wave functions [ 48 ] and the bakulev - mikhailov - stefanis pion wave function [ 49 ] .
it should be noted that the wave functions of pions also are developed in refs.[50 - 52 ] by the dubna group :    @xmath47,\ ] ] @xmath48,\ ] ] @xmath49,\ ] ] @xmath50,\ ] ] @xmath51 @xmath52 @xmath53 where @xmath54 is the pion decay constant .
here , we have denoted by @xmath55 , the longitudinal fractional momentum carried by the quark within the meson .
then , @xmath56 and @xmath57 .
the pion wave function is symmetric under the replacement @xmath58 .
several important nonperturbative tools have been developed which allow specific predictions for the hadronic wave functions directly from theory and experiments . the qcd sum - rule technique and lattice gauge theory
provide constraints on the moments of the hadronic distribution amplitude .
however , the correct pion wave function is still an open problem in qcd .
it is known that the pion wave function can be expanded over the eigenfunctions of the one - loop brodsky - lepage equation , _
i.e. _ , in terms of the gegenbauer polynomials @xmath59    @xmath60.\ ] ]    the evolution of the wave function on the factorization scale @xmath61 is governed by the functions @xmath62 , @xmath63^{\gamma_n/\beta_0},\ ] ] @xmath64    in eq.(2.6 ) , @xmath65 are anomalous dimensions defined by the expression ,    @xmath66.\ ] ]    the constants @xmath67 are input parameters that form the shape of the wave functions and which can be extracted from experimental data or obtained from the nonperturbative qcd computations at the normalization point @xmath68 .
the qcd coupling constant @xmath13 at the one - loop approximation is given by the expression ,    @xmath69    here , @xmath70 is the fundamental qcd scale parameter , @xmath71 is the qcd beta function one - loop coefficient , respectively , @xmath72 the cross section for the higher - twist subprocess @xmath73 is given by the expression @xmath74 ^ 2}{{\hat s}^3}\left[\frac{1}{{\hat u}^2}+\frac{1}{{\hat t}^2}\right],\ ] ] where @xmath75+e_2\hat u\int_{0}^{1}dx\left[\frac{\alpha_{s}(\hat q_2 ^ 2)\phi_{\pi}(x,\hat q_2 ^ 2)}{1-x}\right].\ ] ] here @xmath76and @xmath77 , represent the momentum squared carried by the hard gluon in fig.1 , @xmath78 is the charge of @xmath79 and @xmath80 .
the higher - twist contribution to the large-@xmath3 pion production cross section in the process @xmath81 is [ 53 ] : @xmath82 @xmath83 @xmath84 @xmath85 @xmath86 @xmath87 @xmath88 @xmath89 @xmath90    where @xmath91  is the transverse mass of pion , which is given by @xmath92
let us first consider the frozen coupling approach . in this approach
, we take the four - momentum square @xmath93 of the hard gluon to be equal the pion s transverse momentum square @xmath94 . in this case , the qcd coupling constant @xmath95 in the integral ( 2.10 ) does not depend on the integration variable . after
this substitution calculation of integral ( 2.10 ) becomes easy .
hence , the effective cross section obtained after substitution of the integral ( 2.10 ) into the expression ( 2.9 ) is referred as the frozen coupling effective cross section .
we will denote the higher - twist cross section obtained using the frozen coupling constant approximation by @xmath96 .
we have extracted the following higher - twist subprocesses contributing to the two covariant cross sections in eq.(2.11 ) @xmath97 by charge conjugation invariance , we also have @xmath98
in this section , we shall calculate the integral ( 2.10 ) using the running coupling constant method and also discuss the problem of normalization of the higher - twist process cross section in the context of the same approach .    as is seen from ( 2.10 ) , in general , one has to take into account not only the dependence of @xmath99 on the scale @xmath100 , but also an evolution of @xmath101 with @xmath100 .
the meson wave function evolves in accordance with a bethe - salpeter - type equation .
therefore , it is worth noting that , the renormalization scale ( argument of @xmath95 ) should be equal to @xmath102 , @xmath103 , whereas the factorization scale [ @xmath61 in @xmath104 is taken independent from @xmath105 , we assume @xmath106 . such a approximation does not considerably change the numerical results , but the phenomenon considered in this article ( effect of infrared renormalons ) becomes transparent .
the main problem in our investigation is the calculation of the integral in ( 2.10 ) by the running coupling constant approach .
this integral in the framework of the running coupling approach takes the form    @xmath107    the @xmath108 has the infrared singularity at @xmath109 , if @xmath110 or @xmath111 , if @xmath112 and as a result integral @xmath113 diverges ( the pole associated with the denominator of the integrand is fictitious , because @xmath114 , and therefore , the singularity of the integrand at @xmath115 is caused only by @xmath116 ) . for the regularization of the integral , we express the running coupling at scaling variable @xmath116 with the aid of the renormalization group equation in terms of the fixed one @xmath13 .
the renormalization group equation for the running coupling @xmath117 has the form [ 36 ]    @xmath118 ^ 2-\frac{\beta_{1}}{16}[\alpha(\lambda q^2)]^3\ ] ]    where @xmath119 the solution of eq.(3.2 ) , with the initial condition @xmath120 is [ 35 ]    @xmath121^{-1}.\ ] ]    this transcendental equation can be solved iteratively by keeping the leading @xmath122 and next - to - leading @xmath123 powers . for @xmath124
these terms are given by @xmath125}{[1+\ln\lambda / t]^2}\ ] ]    the first term in eq.(3.4 ) is the solution of the renormalization group eq.(3.2 ) with leading power accuracy , whereas the whole expression ( 3.4 ) is the solution of eq.(3.2 ) with next - to - leading power accuracy . after substituting eq.(3.4 ) into eq.(2.10 ) we get : @xmath126 @xmath127 @xmath128 @xmath129^{\gamma_{n}/\beta_{0}}c_{n}^{3/2}(2x-1)\right]}{(1-x ) ( 1+\ln\lambda / t_{1})}-\ ] ] @xmath130^{\gamma_{n}/\beta_{0}}c_{n}^{3/2}(2x-1)\right ] \ln(1+\ln\lambda / t_{1})}{(1-x ) ( 1+\ln\lambda / t_{1})^2}+\ ] ] @xmath131^{\gamma_{n}/\beta_{0}}c_{n}^{3/2}(2x-1)\right]}{(1-x ) ( 1+\ln\lambda / t_{2})}-\ ] ] @xmath132^{\gamma_{n}/\beta_{0}}c_{n}^{3/2}(2x-1)\right]\ln(1+\ln\lambda / t_{2})}{(1-x ) ( 1+\ln\lambda / t_{2})^2}\ ] ]    where @xmath133 and @xmath134 .
it should be note that in [ 41 ] , we used @xmath135 .
the integral ( 3.5 ) is common and , of course , still divergent , but now it is recast into a form , which is suitable for calculation . using the running coupling constant approach ,
this integral may be found as a perturbative series in @xmath13 @xmath136 the expression coefficients @xmath137 can be written as power series in the number of light quark flavors or , equivalently , as a series in power of @xmath71 : @xmath138 the coefficients @xmath139 of this series demonstrate factorial growth @xmath140 , which might indicate an infrared renormalon nature of divergences in the integral ( 3.5 ) and corresponding series ( 3.6 ) .
the procedure for dealing with such ill - defined series is well known : one has to perform the borel transform of the series [ 15 ] @xmath141(u)=\sum_{n=0}^{\infty}\frac{d_n}{n ! } u^n,\ ] ] then invert @xmath142(u)$ ] to obtain the resummed expression ( the borel sum ) @xmath143 .
after this we can find directly the resummed expression for @xmath144 .
the change of the variable @xmath105 to @xmath145 , as @xmath146 .
then @xmath147 @xmath148    the basic theoretical problem is how to define the borel sum @xmath142(u)$ ] of the integral in eq.(3.7 ) for the quantities we are interested in . in qcd
this problem is usually solved using perturbative methods and calculating the corresponding multiloop feynman diagrams with a one - gluon line , dressed by the chains of fermion bubbles .    for the calculation
the expression ( 3.7 ) we will apply the integral representations of @xmath149 and @xmath150 [ 54,55 ] . after this operation , formula ( 3.7 ) is simplified and we can extract the borel sum of the perturbative series ( 3.6 ) and the corresponding borel transform in dependence from the wave functions of the meson , respectively . also after such manipulations the obtained expression can be used for numerical computations .
it is convenient to use the following integral representation for @xmath149 and @xmath151 @xmath152 and @xmath153 where @xmath154 is the euler - mascheroni constant .
after inserting eq.(3.8 ) and ( 3.9 ) into ( 3.7 ) .
then , we obtain @xmath155 @xmath156 @xmath157 @xmath158 in the case of @xmath159 for @xmath143 , we get @xmath160-\frac{4\sqrt{3 } \pi f_{\pi } e_{1 } \beta_{1 } \hat t}{\beta_{0}^3 } \cdot\ ] ] @xmath161(1-c-\ln u-\ln { t_1})u+\ ] ] @xmath162-\frac{4\sqrt{3 } \pi f_{\pi } e_{2}\beta_{1 } \hat u}{\beta_{0}^3}\cdot\ ] ] @xmath163(1-c-\ln u-\ln { t_2})u.\ ] ] in the case of the @xmath164 wave function , we find @xmath165 @xmath166^{50/81 } \left[\frac{4}{1-u}-\frac{24}{2-u}+\frac{40}{3-u}- \left.\frac{20}{4-u}\right]\right]-\ ] ] @xmath167 @xmath166^{50/81 } \left[\frac{4}{1-u}-\frac{24}{2-u}+\frac{40}{3-u}- \left.\frac{20}{4-u}\right]\right](1-c-\ln u-\ln { t_1})u+\ ] ] @xmath168 @xmath166^{50/81 } \left[\frac{4}{1-u}-\frac{24}{2-u}+\frac{40}{3-u}- \left.\frac{20}{4-u}\right]\right]-\ ] ] @xmath169 @xmath166^{50/81 } \left[\frac{4}{1-u}-\frac{24}{2-u}+\frac{40}{3-u}- \left.\frac{20}{4-u}\right]\right](1-c-\ln u-\ln { t_2})u\ ] ]    in the case of the @xmath170 wave function , we get @xmath171^{50/81}\cdot\ ] ] @xmath172 - 0.4125\left[\frac{\alpha_{s}(p_{t}^2)}{\alpha_{s}(\mu_{0}^2)}\right]^{364/405}\cdot\ ] ] @xmath173\right]-\ ] ] @xmath174^{50/81}\cdot\ ] ] @xmath172 - 0.4125\left[\frac{\alpha_{s}(q^2)}{\alpha_{s}(\mu_{0}^2)}\right]^{364/405}\cdot\ ] ] @xmath173\right ] ( 1-c-\ln u-\ln { t_1})u+\ ] ] @xmath175^{50/81}\cdot\ ] ] @xmath172 - 0.4125\left[\frac{\alpha_{s}(p_{t}^2)}{\alpha_{s}(\mu_{0}^2)}\right]^{364/405}\cdot\ ] ] @xmath173\right]-\ ] ] @xmath176^{50/81}\cdot\ ] ] @xmath172 - 0.4125\left[\frac{\alpha_{s}(p_{t}^2)}{\alpha_{s}(\mu_{0}^2)}\right]^{364/405}\cdot\ ] ] @xmath173\right ] ( 1-c-\ln u-\ln { t_2})u\ ] ]    also , in the case of the @xmath177 wave function , we get @xmath178^{50/81}\cdot\ ] ] @xmath172 - 0.244\left[\frac{\alpha_{s}(p_{t}^2)}{\alpha_{s}(\mu_{0}^2)}\right]^{364/405}\cdot\ ] ] @xmath173\right]-\ ] ] @xmath179^{50/81}\cdot\ ] ] @xmath172 - 0.244\left[\frac{\alpha_{s}(p_{t}^2)}{\alpha_{s}(\mu_{0}^2)}\right]^{364/405}\cdot\ ] ] @xmath173\right](1-c-\ln u-\ln { t_1})u+\ ] ] @xmath180^{50/81}\cdot\ ] ] @xmath172 - 0.244\left[\frac{\alpha_{s}(p_{t}^2)}{\alpha_{s}(\mu_{0}^2)}\right]^{364/405}\cdot\ ] ] @xmath173\right]-\ ] ] @xmath181^{50/81}\cdot\ ] ] @xmath172 - 0.244\left[\frac{\alpha_{s}(p_{t}^2)}{\alpha_{s}(\mu_{0}^2)}\right]^{364/405}\cdot\ ] ] @xmath173\right ] ( 1-c-\ln u-\ln { t_2})u\ ] ]    equations(3.11)-(3.14 ) is nothing more than the borel sum of the perturbative series ( 3.6 ) , and the corresponding borel transform in the case @xmath159 is @xmath141(u)=\frac{1}{1-u}-\frac{1}{2-u},\ ] ] in the case @xmath164 is @xmath141(u)=\frac{1}{1-u}-\frac{1}{2-u}+0.84\left(\frac{\alpha_{s}(p_{t}^2)}{\alpha_{s}(\mu_{0}^2)}\right)^{50/81 } \left(\frac{4}{1-u}-\frac{24}{2-u}+\frac{40}{3-u}-\frac{20}{4-u}\right),\ ] ]    in the case @xmath170 is @xmath141(u)=\frac{1}{1-u}-\frac{1}{2-u}+0.405\left(\frac{\alpha_{s}(p_{t}^2)}{\alpha_{s}(\mu_{0}^2)}\right)^{50/81 } \left(\frac{4}{1-u}-\frac{24}{2-u}+\frac{40}{3-u}-\frac{20}{4-u}\right)-\ ] ] @xmath182 and in the case @xmath177 is @xmath141(u)=\frac{1}{1-u}-\frac{1}{2-u}+0.282\left(\frac{\alpha_{s}(p_{t}^2)}{\alpha_{s}(\mu_{0}^2)}\right)^{50/81 } \left(\frac{4}{1-u}-\frac{24}{2-u}+\frac{40}{3-u}-\frac{20}{4-u}\right)-\ ] ] @xmath183 the series ( 3.6 ) can be recovered by means of the following formula @xmath184(u)\mid_{u=0}\ ] ] the borel transform @xmath142(u)$ ] has poles on the real @xmath185 axis at @xmath186 which confirms our conclusion concerning the infrared renormalon nature of divergences in ( 3.6 ) . to remove them from eqs.(3.11)-(3.14 )
we applied the principal value prescription .
hence , the effective cross section obtained after substitution of the expressions ( 3.10 - 3.14 ) into the expression ( 2.10 ) is referred as the running coupling effective cross section .
we will denote the higher - twist cross section obtained using the running coupling constant approach by @xmath187
regarding the higher - twist corrections to the pion production cross section , a comparison of our results with leading - twist contributions is crucial .
we take two leading - twist subprocesses for the pion production:(1 ) quark - antiquark annihilation @xmath188 , in which the @xmath189 meson is indirectly emitted from the gluon , @xmath190 and ( 2 ) quark - gluon fusion , @xmath191 , with subsequent fragmentation of the final quark into a meson , @xmath192 .
the corresponding cross sections are obtained in @xmath193 @xmath194 for the leading - twist contribution , we find @xmath195 @xmath196 where @xmath197 @xmath198 and @xmath199 represents gluon and quark fragmentation functions into a meson containing gluon and quark of the same flavor . in the leading - twist subprocess , the @xmath189 meson is indirectly emitted from the gluon and quark with the fractional momentum @xmath45 .
the @xmath200 function can be expressed in terms of the parton kinematic variables , and the @xmath45 integration can then be done .
the final form for the cross section is @xmath201 @xmath202 @xmath203 @xmath204
in this section , we discuss the numerical results for next - to - leading order higher - twist effects with higher - twist contributions calculated in the context of the running coupling and frozen coupling approaches on the dependence of the chosen meson wave functions in the process @xmath205 . in the calculations , we use the asymptotic wave function @xmath206 , the chernyak - zhitnitsky @xmath207 , the cleo pion wave function [ 41 ] , the braun - filyanov pion wave functions [ 42 ] , and the bakulev - mikhailov - stefanis pion wave function [ 43 ] . for the higher - twist subprocess , we take @xmath208 and we have extracted the following four higher - twist subprocesses contributing to @xmath209 cross sections : @xmath210 @xmath211 , @xmath212 , @xmath213 contributing to cross sections . for the dominant leading - twist subprocess for the pion production , we take the quark - antiquark annihilation @xmath188 , in which the @xmath189 meson is indirectly emitted from the gluon and quark - gluon fusion , @xmath214 , with subsequent fragmentation of the final quark into a meson , @xmath215 .
as an example for the quark distribution function inside the proton , the mrst2003c package [ 56 ] has been used . the higher twist subprocesses probe the meson wave functions over a large range of @xmath61 squared momentum transfer , carried by the gluon .
therefore , in the diagram given in fig.1 we take @xmath216 , @xmath103 , which we have obtained directly from the higher - twist subprocesses diagrams .
the same @xmath61 has been used as an argument of @xmath30 in the calculation of each diagram .
the results of our numerical calculations are plotted in figs.2 - 15 .
first of all , it is very interesting to compare the resummed higher - twist cross sections with the ones obtained in the framework of the frozen coupling approach . in figs.2 - 4
we show the dependence of next - to - leading order higher - twist cross sections @xmath217 calculated in the context of the running coupling constant approach and the ratios @xmath218/@xmath219 , @xmath217/@xmath220 as a function of the pion transverse momentum @xmath3 for different pion wave functions at @xmath221 .
it is seen that the values of cross sections @xmath217 , and ratios for fixed @xmath222 and @xmath223 depend on the choice of the pion wave function .
as seen from fig.2 the next - to - leading order higher - twist differential cross section is monotically decreasing with an increase in the transverse momentum of the pion . in figs.3 - 5
, we show the dependence of the ratios @xmath218/@xmath219 , @xmath217/@xmath220 , @xmath217/@xmath224 as a function of the pion transverse momentum @xmath3 for different pion wave functions . here
@xmath217 , @xmath225 , @xmath224 are the higher - twist cross sections calculated in the context of the running coupling method , in the framework of the frozen coupling approach and is the leading - twist cross section , respectively . as seen from fig.5 , in the region @xmath226
next - to - leading order higher - twist cross section calculated in the context of the running coupling method is suppressed by about 1 - 2 orders of magnitude relative to the leading - twist cross section , but in the region @xmath227 is comparable with the cross section of leading - twist . in figs.6 - 8 , we have depicted next - to - leading order higher - twist cross sections@xmath217 , ratios @xmath228 , @xmath229 as a function of the rapidity @xmath222 of the pion at @xmath16 and @xmath230 . at @xmath16 and @xmath230 ,
the pion rapidity lies in the region @xmath231 .
as seen from figs.6 - 8 , next - to - leading order higher - twist cross section and ratios have a different distinctive . in the region ( @xmath232 )
the next - to - leading order higher - twist cross section @xmath217 , is suppressed by about 1 - 2 order of magnitude relative to the leading - twist cross section . in the region ( @xmath233 )
, the ratio for all wave functions increase with an increase of the @xmath222 rapidity of the pion and has a maximum approximately at the point @xmath234 .
besides that , the ratio decreases with an increase in the @xmath222 rapidity of the pion . as is seen from figs.7 - 8 , the ratio @xmath235 very sensitive to the choice of the meson wave functions .
but , as seen from fig.8 , the ratio @xmath217/@xmath224 for all wave functions has a minimum approximately at the point @xmath234 .
also , the distinction between @xmath236 with @xmath237 , @xmath238 , @xmath239 and @xmath240 have been calculated .
for example , in the case of @xmath16 , @xmath221 , the distinction between @xmath236 with @xmath241(i = cleo , cz , bf , bms ) as a function of the pion transverse momentum @xmath3 is shown in table [ table1 ] .
thus , the distinction between @xmath236 and @xmath242 is maximum at @xmath243 , with @xmath244 at @xmath245 but the distinction between @xmath236 with @xmath246 is minimum at @xmath245 , with @xmath244 at @xmath243 and increase with an increase in @xmath0 .
such a behavior of @xmath235 may be explained by reducing all moments of the pion model wave functions to those of @xmath159 for high @xmath61 .
also , we have calculated the distinction between @xmath236 with @xmath237 , @xmath238 , @xmath239 and @xmath240 as a function of the rapidity @xmath222 of the pion .
for example , in the case of @xmath247 , @xmath248 the distinction between @xmath236 with @xmath241(i = cleo , cz , bf , bms ) as a function of the rapidity @xmath222 of the pion is presented in table [ table2 ] .
we have also carried out comparative calculations in the center - of - mass energy @xmath17 .
the results of our numerical calculations in the center - of - mass energies @xmath249 are plotted in figs.9 - 15 .
analysis of our calculations at the center - of - mass energies @xmath16 and @xmath17 , show that with the increase in beam energy values of the cross sections , ratio @xmath250 , and contributions of next - to - leading order higher - twist to the cross section decrease by about 1 - 2 order .
therefore the experimental investigation of higher - twist effects include renormalon effects conveniently in low energy .
on the other hand , the higher - twist corrections and ratios @xmath235 a very sensitive to the choice of the pion wave function . analysis of our calculations at the center - of - mass energies @xmath16 and @xmath17 , show that values of the next - to - leading order cross sections decrease by about @xmath251 percent of magnitude relative to the leading - order cross sections .
also , the distinction between @xmath236 with @xmath237 , @xmath238 , @xmath239 and @xmath240 have been calculated .
for example , in the case of @xmath17 , @xmath221 , the distinction between @xmath236 with @xmath241(i = cleo , cz , bf , bms ) as a function of the pion transverse momentum @xmath3 is shown in table [ table3 ] .
thus , the distinction between @xmath236 with @xmath241,(i= cleo , bf , ) is maximum at @xmath252 , with @xmath244 at @xmath253 , but the distinction between @xmath236 with @xmath238 , @xmath237 , @xmath239 is minimum at @xmath253 , with @xmath244 at @xmath254 and increase with an increase in @xmath0 .
also , we have calculated the distinction between @xmath236 with @xmath237 , @xmath238 , @xmath239 and @xmath240as a function of the rapidity @xmath222 of the pion .
for example , in the case of @xmath255 , @xmath256 the distinction between @xmath236 with @xmath241 as a function of the rapidity @xmath222 of the pion is presented in table [ table4 ] .
the calculations show that the ratio @xmath241/@xmath236 , ( i = cleo , cz , bf , bms ) for all values of the transverse momentum @xmath0 of the pion identically equivalent to ratio @xmath257/@xmath258 .    in our calculations of the next - to - leading order higher - twist cross section of the process the dependence of the transverse momentum of
meson appears in the range of @xmath259 .
therefore , higher - twist cross section obtained in our paper should be observable at rhic
in this work we have calculated the inclusive meson production via higher - twist mechanism and obtained the expressions for the subprocess @xmath260 cross section for mesons with symmetric wave functions . for calculation of the next - to - leading order cross section
we have applied the running coupling constant method and revealed infrared renormalon poles in the cross section expression .
infrared renormalon induced divergences have been regularized by the means of the principal value prescripton and the resummed expression ( the borel sum ) for the next - to - leading order higher - twist cross section has been found . in the resummed higher - twist cross section
differs considerably from that found using the frozen coupling approach , some regions .
also we have demonstrated that next - to - leading order higher - twist contributions to meson production cross section in the proton - proton collisions have important phenomenological consequences .
our investigation enables us to conclude that the higher - twist pion production cross section in the proton - proton collisions depends on the form of the pion model wave functions and may be used for their study .
analysis of our calculations shows that the magnitude of next - to - leading order cross sections calculated in the running coupling approach in some regions is larger than the leading - twist cross sections in 1 - 2 order .
further investigations are needed in order to clarify the role of high twist effects in this process .
we have demonstrated that the resummed result depends on the pion model wave functions used in calculations .
the proton - proton collisions provide us with a new opportunity to probe a proton s internal structure .
in particular , meson production in proton - proton collisions takes into account infrared renormalon effects : this opens a window toward new types of parton distributions which can not be measured by the deep inelastic lepton - proton scatterings . finally , we discuss the phenomenological consequences of possible higher - twist contributions to meson production in proton - proton collisions in next - to - leading order at rhic .
future rhic measurements will provide further tests of the dynamics of large @xmath0 hadron production beyond leading twist .
the work presented in this paper was completed while one of the authors , a.i.ahmadov , was visiting the hecap section of the abdus salam ictp , trieste , italy .
he would like to express his gratitude to the members of the section especially to the head prof .
s. randjbar - daemi for their hospitality
. financial support by ictp is also gratefully acknowledged .
hooft , in the why s of subnuclear physics , erice , 1977 , edited by a.  zichichi ( plenum , new york , 1979 ) , p.94 a. h.  mueller , nucl . phys . *
b250 * , 327 ( 1985 ) ; phys . lett . *
b308 * , 355 ( 1993 ) .
v. i.  zakharov , nucl .
b385 * , 452 ( 1992 ) .
m.  beneke , phys .
rep . * 317 * , 1 ( 1999 ) .
g.  grunberg , phys . lett . *
b95 * , 70 ( 1980 ) ; * b110 * , 501(e ) , ( 1982 ) ; phys . rev . *
d29 * , 2315 ( 1984 ) .
p. m.  stevenson , phys . rev . *
d23 * , 2916 ( 1981 ) .
s. j.  brodsky , g. p.  lepage and p. b.  mackenzie , phys .
* d28 * , 228 ( 1983 ) ; g. p.  lepage and p. b.  mackenzie , ibid . * 48 * , 2250 ( 1993 ) .
g.  grunberg and a. l.  kataev , phys .
lett * b279 * , 352 ( 1992 ) . s. j.  brodsky and j.  rathsman , hep - ph/9906339 .
h. j.  lu , and c. a. r. sa de  melo , phys .
lett * b273 * , 260 ( 1991 ) ; * b285 * , 399(e ) ( 1992 ) ; h. j.  lu , phys . rev . *
d45 * , 1217 ( 1992 ) . c. p.  lepage and p. b.  mackenzie phys .
* d48 * , 2250 ( 1993 ) . s. j.  brodsky and h. j.  lu , phys . rev .
* d51 * , 3652 ( 1995 ) .
j.  rathsman , phys . rev . *
d54 * , 3420 ( 1996 ) .
m.  neubert , phys . rev . *
d51 * , 5924 ( 1995 ) .
m.  beneke and v. m.  braun , phys . lett . * b348 * , 513 ( 1995 ) ; p. ball , m.  beneke and v. m.  braun , nucl . phys . *
b452 * , 424 ( 1995 ) ; m.  beneke , nucl . phys . *
b405 * , 563 ( 1993 ) . c. n.  lovett - turner and c. j. maxwell , nucl . phys .
b432 * , 147 ( 1994 ) . s. j.  brodsky , j.  ellis , e.  gardi , m.  karliner , and m. a.  samuel , phys .
rev . * d56 * , 6980 ( 1997 ) .
s. j.  brodsky , g. t.  gabadadze , a. l.  kataev , and h. j.  lu , phys .
b372 * , 133 ( 1996 ) .
s. j.  brodsky , m.  melles , and j.  rathsman , phys .
rev . * d60 * , 096006 ( 1999 ) .
s. j.  brodsky , m. s.  gill , m.  melles , and j.  rathsman , phys .
rev . * d58 * , 116006 ( 1998 ) .
s. v.  mikhailov , jhep 0706:009 ( 2007 ) .
v. a.  matveev , r. m.  muradyan , and a. n.  tavkhelidze , lett .
nuovo cimento 7 , 719 ( 1973 ) .
s. j.  brodsky and g. r.  farrar , phys .
lett . * 31 * , 1153 ( 1973 ) . j. f.  gunion , s. j.  brodsky , r.  blankenbecler , phys.rev .
* d6 * , 2652 ( 1972 ) d. w. sivers , s. j. brodsky , r. blankenbecler , phys . rept .
* 23 * , 1 ( 1976 ) .
v. a.  matveev , l. a.  slepchenko , and a. n.  tavkhelidze , phys . lett . *
b100 * , 75 ( 1981 ) .
g. l.  lepage and s. j.  brodsky , phys . rev . * d22 * , 2157 ( 1980 ) . v. l.  chernyak and a. r.  zhitnitsky , phys
. rep . * 112 * , 173 ( 1984 ) .
franois  arleo and david  denterria , phys
* d78 * : 094004 ( 2008 ) .
j. a.  bagger and j. f.  gunion , phys
* d25 * , 2287 ( 1982 ) . v. n.  baier and a. grozin , phys . lett . *
b96 * , 181 ( 1980 ) ; s.  gupta , phys .
rev . * d24 * , 1169 ( 1981 ) .
a. i.  ahmadov , i.  boztosun , r. kh .
muradov , a.  soylu and e. a.  dadashov , int . j. mod
e15 * , 1209 ( 2006 ) .
m.  maul , e.  stein , a.  schfer , l.  mankiewich , phys .
lett.*b401*,100 ( 1997 ) .
y. l.  dokshitzer , v. a.  khoze and s. i.  troyan , phys . rev . *
d53 * , 89 ( 1996 ) .
a. l.  kataev , mod .
* a20 * , 2007 ( 2005 ) .
h.  contopanagos and g.  sterman , nucl . phys . * b419 * , 77 ( 1994 ) .
w.  greiner , s.  schramm and e.  stein , quantum chromodynamics , 2nd edn.(berlin , springer , 2002 ) , pp.551 .
d. v.  shirkov and i. l.  solovtsov , phys .
lett . * 79 * , 1209 ( 1997 ) . star collaboration ( j.adams et.al ) , phys
b637 * , 161 ( 2006 ) .
agaev , phys . lett . *
b360 * , 117 ( 1995 ) ; * b369 * , 379(e ) ( 1996 ) . a. i.  ahmadov , coskun  aydin , sh . m.  nagiyev , yilmaz a.  hakan , and e. a.  dadashov , phys
* d80 * , 016003 ( 2009 ) .
a. i.  ahmadov , coskun  aydin , e. a.  dadashov and sh .
m.  nagiyev , phys . rev . *
d81 * , 054016 ( 2010 ) .
agaev , modern physics letters * a13 * , 2637 ( 1998 ) .
franois  arleo , stanley j.  brodsky , dae suang hwang , and anne m.sickles , phys
* 105 * , 062002 ( 2010 ) . s. m.  berman , j. d.  bjorken and j. b.  kogut , phys . rev . *
d18 * , 900 ( 1978 ) .
v. l.  chernyak and a. r.  zhitnitsky , nucl . phys . *
b201 * , 492 ( 1982 ) ; v. l.  chernyak , a. r.  zhitnitsky and i. r.  zhitnitsky , nucl . phys . *
b204 * , 477 ( 1982 ) .
a.  schmedding and o.  yakovlev , phys
d62 * , 116002 ( 2000 ) . v. m.  braun and i. f.  filyanov , z. phys , * c44 * , 157 ( 1989 ) .
a. p.  bakulev , s. v.  mikhailov and n. g.  stefanis , phys .
lett , * b587 * , 91 ( 2004 ) . s. v.  mikhailov and a. v.  radyushkin , phys . rev . *
d45 * , 1754 ( 1992 ) .
bakulev , s. v.  mikhailov , z. phys . *
c68 * , 451 ( 1995 ) .
a. e.  dorokhov , jetp lett .
* 77 * , 63 ( 2003 ) .
j. f.  owens , rev .
phys . * 59 * , 465 ( 1987 ) .
j. zinn - justin , phys .
rept . * 70 * , 109 ( 1981 ) .
a.  erdelyi , higher transcendental functions ( mcgrow - hill book company , new york , 1953 ) , vol.2 .
mrst2003c.f can be obtained from http://durpdg.dur.ac.uk/hepdata/pdf.html .
see also , a.  d.  martin , r.  g.  roberts , w.  j.  stirling and r. s.thorne , hep - ph/0307262 ; r. s.thorne , hep - ph/0309343 .


post - main - sequence mass loss affects both present and future generations of stars .
it determines the remnant state of stars of intermediate mass and below , as well as affecting the state and contents of the interstellar medium out of which new stars are formed . however , the rate and amount of mass lost by red giant branch ( rgb ) and asymptotic giant branch ( agb ) stars is a major uncertainty in stellar evolutionary theory .
as @xcite and @xcite point out , existing formulations rely on empirical laws based on observations of population  i stars .
there is little theoretical guidance on how mass loss should be included in stellar evolutionary models , or on how it affects population synthesis models based on the stellar models .
globular clusters ( gcs ) have long been the prime testing grounds for stellar evolution theory . with large populations of low - mass stars of the same age , metallicity , and distance
, they should be ideal places to test ideas about mass loss .
there are several observational confirmations of the theoretical prediction that stars more massive than @xmath2 lose a significant fraction of their mass after leaving the main sequence .
these include initial - final mass estimates of white dwarfs ( e.g. * ? ? ?
* ) and evidence from line profiles of red giants ( e.g. , * ? ? ?
to explain the horizontal branch it was argued ( e.g. , * ? ? ?
* ) that much of this mass is shed on the rgb , but as agb stars , post - agb stars , and pne are also observed in gcs @xcite some mass must be lost on the agb as well .
there is also ample observational evidence that many gcs have stars which are currently losing mass , at rates of @xmath3  yr@xmath4 @xcite . from these observations
it is difficult to obtain an overall picture of time - averaged mass loss rates and duty cycles and their relation to cluster properties such as metallicity , since only a few stars at a time lose mass .
a more global picture of mass loss within gcs might be provided by the intra - cluster medium ( icm ) .
mass lost from gc stars is expected to accumulate in the icm until the cluster passes through the galactic disk , when the icm will be stripped by ram pressure .
the icm is expected to be mostly gaseous ( although the ionization state is unknown ) , with a small amount of dust formed in the extended atmospheres of red giants .
the dust should be optically thin and heated by starlight to 5080  k @xcite ; its thermal emission should peak in the mid - to - far - infrared ( ir ) and could be detected as an ir excess at these wavelengths .
although gc stars are clearly losing mass , detections of the icm have been surprisingly difficult .
searches for intra - cluster gas have been conducted for several decades @xcite ; searches for dust @xcite are more recent but have yielded similar results , mostly upper limits .
only one galactic gc , ( m15 ) , shows clear evidence for an infrared ( ir ) excess @xcite .
diffuse infrared emission near may be associated with the cluster or with galactic cirrus , while has an infrared point source which may be intracluster dust @xcite . an additional cluster , , has a tentative detection of dust at submillimeter wavelengths @xcite .
all of the detections imply dust masses well below the expected levels .
this paper presents results of a _
@xcite survey of eight milky way gcs using 24 and 70  @xmath0 m  observations with the multiband imaging photometer for _ spitzer _ ( mips ; * ? ? ?
we describe our observational strategy and data analysis techniques in  2 .
section 3 introduces our results , including the detection of excess ir emission in several of the target clusters and upper limits on the remaining clusters .
section 4 explores the implication of these new mid - ir results , including a comparison with predicted dust masses .
we conclude in  5 with a summary and suggested follow - up observations .
the target list ( table  [ tab : targets ] ) comprises the milky way globular clusters that we believed to be the best candidates for icm dust detection .
these clusters are reasonably nearby , at high galactic latitude ( to reduce confusing fore- and back - grounds ) , and have high escape velocities .
they are among the most massive galactic gcs .
more massive clusters are more likely to harbor a larger number of stars in the short - lived superwind phase : if the icm is being removed from the cluster on short time scales , then we are more likely to see it in a more massive cluster where the icm material is being replenished more quickly .
clusters with escape velocities that are much larger than agb wind speeds ( @xmath5  km  s@xmath4 ) will also be more likely to retain icm material . at the time this program was planned , all had existing or scheduled _ spitzer _
mid - ir observations with the infrared array camera ( irac ; * ? ? ?
observing clusters with a range of times since last disk passage is desirable , but we did not use this as a selection criterion because orbital information is only available for a fraction of milky way globular clusters ( e.g. , from * ? ? ?
the clusters span a factor of thirty in metallicity ( 1.5  dex ) , important for testing the effects of metallicity on dust production and retention .
all of the clusters are detected with iras @xcite , and both ngc  104 ( 47 tuc ) and ngc  362 contain stars with known circumstellar dust @xcite .
irac photometry traces the rayleigh - jeans tail of the spectral energy distribution of globular cluster stars .
although cold intracluster dust is not expected to radiate at irac wavelengths , understanding the near - ir sed is very useful for quantifying any ir excess from the dust .
irac observations of the clusters were made over the period 2005 july  2006
september as part of _ spitzer _ program identification ( pid ) 20298 ( pi r. rood ) , except for ngc  5904 which was part of pid 3256 ( pi p. goudfrooij ) .
dataset identifiers are given in table  [ tab : iracphot ] .
all but one of the clusters were observed in a single irac field of view ( @xmath6 ) in all four channels , with 90252 frames ( 12  s each ) per position ( 30 frames , 30  s each for ngc  5904 ) .
ngc  104 was observed in a @xmath7 map .
ngc  104 and ngc  6752 were observed in the high - dynamic - range ( hdr ) mode which includes short ( 0.6  s ) frames taken before the 12  s frames .
the half - light radius of ngc  104 ( 190 ) is slightly too large to be fully covered by the pid 20298 observations .
fortunately , observations of this cluster covering a wider field of view were available , taken as part of instrument commissioning in 2003 november ( pid 623 ) .
these observations were taken in 12-second hdr mode with 3 frames per position , and cover a @xmath8 area .
for all clusters , the _ spitzer _ science center ( ssc ) pipeline version s14 post - bcd mosaics were used for integrated photometry .
apphot _ was used to measure integrated magnitudes within the half - light radii given in table  [ tab : targets ] . for each cluster
the background was measured in an annulus with inner and outer radii 108 and 117 , respectively , and subtracted from the integrated magnitudes . for ngc  104 and
ngc  6752 , the mosaics made from the short - exposure frames were used to avoid bright star saturation ; the other clusters had no saturation problems .
the magnitudes measured in the half - light radii were converted to total flux densities with multiplication by 2 , and corrected with application of the ` extended source calibration ' given by the ssc .
the results are given in table  [ tab : iracphot ] , with electronic links to the observational data .
the uncertainties in the flux densities due to poisson noise from the cluster and sky background noise are very small ( @xmath9% ) ; uncertainties in the absolute ( a few percent ; * ? ? ?
* ) and extended source calibrations ( 10% ) dominate the total flux measurements .
the mips photometry astronomical observing template ( aot ) was used to produce @xmath6 images of the clusters at 24 and 70  @xmath0 m .
this field size is comparable to the cluster sizes ( half - light radii from * ? ? ?
* @xmath10 ) .
although a somewhat larger field of view would have been desirable , this would have substantially increased the observing time . to set exposure times , we predicted the 24 and 70  @xmath0 m  stellar fluxes from the clusters @xmath11-band fluxes ( measured from 2mass images ; * ? ? ?
* ) assuming a rayleigh - jeans @xmath12 spectrum .
the 24  @xmath0 m  observations were required to characterize the cluster sed , so high - significance detections of individual stars were not required .
the estimated 24  @xmath0 m  flux was used to compute the expected average 24  @xmath0 m  surface brightness inside the clusters half - light radii .
the resulting values , 0.84.2  mjy  sr@xmath4 , were expected to be detectable in 1 cycle of the photometry aot with 10 second frame times .
the intra - cluster dust is expected to emit most strongly at 70  @xmath0 m ; however , even at this long wavelength , the integrated flux of the cluster stars is much brighter than the mips detection limits .
we therefore targeted the exposure times to have strong detections of the integrated stellar photospheric flux ; any dust signal would appear as an ir excess at 70  @xmath0 m .
the resulting predicted 70  @xmath0 m  fluxes ranged from 6150 mjy , requiring 110 cycles of the mips photometry aot with 10 second frame times .
the observations were made over the period 2006 june
2006 november under _
pid 30031 , a guaranteed time program of g.  fazio .
the data were processed with the s14 version of the mips pipeline .
post - pipeline processing for the 24  @xmath0 m  bcd images included ` self - calibration ' ( dividing by the normalized median image of ` off - source ' frames ) and mosaicking ; the 70  @xmath0 m  bcd images were time- and column - filtered and then mosaicked .
all mosaicking used the mopex program @xcite with the standard parameters recommended in ssc documentation .
the 24  @xmath0 m  mosaics have a pixel scale of 2.5 , and the 70  @xmath0 m  mosaics have a scale of 5.0 .
figures  [ fig : images24 ] and [ fig : images70 ] show the cluster mosaic images .
all of the clusters were detected at high signal - to - noise in the 24  @xmath0 m  images .
the same is not true at 70  @xmath0 m , and the major result of this paper is the lack of any strong dust signal from any of the eight clusters .    as for irac ,
flux densities in the mips bands were measured with aperture photometry .
inspection of figure  [ fig : images24 ] shows that the mid - infrared light for several clusters ( ngc  5272 and ngc  6205 ) appears to have slightly different centers from the optical centers given by @xcite .
we computed new centers by smoothing the 24  @xmath0 m  images and computing the flux centroid of the smoothed images in a box of width 4 times the smoothed fwhm , centered on the optical center .
the mid - ir centers are given in table  [ tab : mipsphot ] .
compared to the irac data , the choice of aperture for the mips images is less obvious : the intracluster medium seems more likely to sink to the center of the cluster potential well than to follow the same spatial distribution as the cluster stars .
[ however , we note that both the extended emission detected in ngc  7078 by @xcite and the possible intracluster dust in ngc  6024 reported by @xcite are offset from the cluster centers . ]
we decided to use two apertures for each cluster : the optical half - light radius as for the irac data , and a smaller ` core ' radius , to maximize signal - to - noise for detection of dust at the cluster center .
the half - light radii are given in table  [ tab : targets ] ; for these apertures , background annuli with inner and outer radii 150 and 175  were used .
following @xcite , we chose a core aperture of 16  radius ( diameter twice the fwhm of the psf ) , with background annuli 120140  in radius . to allow for subtraction of stellar contribution to the 70  @xmath0 m  signal
( see  [ sec : nsflux ] ) , the 24  @xmath0 m  photometry in this aperture was measured on versions of the images convolved to the 70  @xmath0 m  resolution using a kernel described in @xcite .
the results of the aperture photometry are listed in table  [ tab : mipsphot ] , with electronic links to the observational data ; no color corrections have been made to these values .
the only aperture correction made is a multiplicative factor of 1.13 applied to the large - aperture 70  @xmath0 m  photometry .
again following @xcite , uncertainties in the mips flux density were calculated as @xmath13 where @xmath14 accounts for the noise correlation between pixels due to re - sampling during mosaicing ( @xmath15 for 24  @xmath0 m  and 2 for 70  @xmath0 m ) , @xmath16 is discussed below , @xmath17 is the standard deviation of the pixels in the background annulus , and @xmath18 and @xmath19 are the number of pixels in the photometry and background apertures , respectively .
( poisson noise from the clusters is negligable compared to the background noise . ) after examining noise in their 70  @xmath0 m  background and photometry apertures , @xcite concluded that ` excess noise ' corresponding to @xmath20 was present in their photometry apertures , a conclusion validated by their distribution of signal - to - noise ratios .
although this precise value may not apply to our observations , which were made in a slightly different manner , we adopt it here to be conservative .
the uncertainties in table  [ tab : mipsphot ] do not include uncertainties in the absolute calibration ( 2% at 24  @xmath0 m  and 5% at 70  @xmath0 m ; * ? ? ?
* ; * ? ? ?
confusion noise due to extragalactic sources is not expected to affect these observations : @xcite derived the 70  @xmath0 m  confusion limit due to extragalactic sources as 0.3  mjy , well below the flux regime of the observations presented here .
confusion due to galactic cirrus , however , is an issue : @xcite computed the mips 70  @xmath0 m  point source detection limits due to galactic cirrus confusion at the latitudes of our clusters ( @xmath21 ) as @xmath22  mjy , again below our detection limits . for larger apertures , however , the confusion noise and detection limits scale as @xmath23 where @xmath24 is the aperture size and @xmath25 is the cirrus power spectrum index .
the scaling factors range from about 6 to several hundred for the half - light radius apertures used , and hence cirrus confusion will affect our ability to detect extended emission for clusters with large angular sizes .
figure  [ fig : seds ] shows the cluster spectral energy distributions ( seds ) from 0.3 to 100  @xmath0 m . the new irac measurements from table  [ tab : iracphot ] and half - light radius mips measurements from table  [ tab : mipsphot ] ( multiplied by 2 to convert to total fluxes , even though this may not be appropriate for non - stellar emission ) are shown as filled symbols .
also plotted are extinction - corrected total optical and near - ir magnitudes given by @xcite and @xcite .
iras measurements from the faint source catalog @xcite are also shown ; the 12 and 25  @xmath0 m  measurements have been corrected for the cluster extent following the procedure used by @xcite and flux densities in all iras bands have been color - corrected using the values for a 5000  k blackbody given in the _ iras explanatory supplement_. upper limits to the 90  @xmath0 m  flux densities given by @xcite are shown for the five clusters in common between the two samples .
overplotted are black - body seds derived from fits to the flux densities ( equally - weighted ) at @xmath26  @xmath0 m .
all are consistent with the expected average temperatures for globular cluster stars ( 45005500  k ) .
a first check of our measurements is to compare them with previous work .
spitzer _ and iras results are mutually consistent except for ngc  6341 , where the 70  @xmath0 m  data are slightly above the iras 60  @xmath0 m  upper limits .
we suggest that the iras limit might be slightly too low if there is extended emission from the cluster . for ngc  104 ,
both the 60 and 70  @xmath0 m  observations are below the predictions based on the black - body fit , possibly because of the large spatial extent of the cluster .
two of our target clusters were observed with the phot instrument on iso by @xcite : using those authors method to convert surface brightness to flux density , their derived 70  @xmath0 m  flux densities are @xmath27  mjy for ngc  104 and @xmath28  mjy for ngc  7252 .
our measurements for these clusters are consistent with these values .
the 70  @xmath0 m  measurements are generally consistent with the 90  @xmath0 m  upper limits given by @xcite , except for ngc  6341 , where the 70  @xmath0 m  value appears to be too large .
the detection for this cluster is at @xmath29 significance and may be spurious .    figure  [ fig : mips_compare ] compares the predicted ( from the black - body fit ) and measured flux densities ( the same values as in the previous figure ) at 24 and 70  @xmath0 m .
we expected that any icm dust would likely be too cool to be detectable at 24  @xmath0 m ; however , circumstellar dust should be warmer and possibly detected at this wavelength .
all of the 24  @xmath0 m  flux densities are within 20% of the prediction ; six clusters are within 10% .
this is reasonable accuracy considering the difficulties in deriving integrated measurements of these resolved objects .
the non - detection of 24  @xmath0 m  excesses in integrated photometry may simply reflect the small number of mass - losing agb stars compared to the bulk of the population .
point - spread - fitting photometry of the 24  @xmath0 m  images , while beyond the scope of this work , may reveal which stars are surrounded by circumstellar dust .
seven of the eight clusters in the sample have measurements or upper limits on 70  @xmath0 m  flux that are more than a factor of 2 above the black - body prediction .
the exception is ngc  104 , for which the 70  @xmath0 m  upper limit is _ lower _ than the predicted photospheric value .
the problem could be background subtraction : the half - light radius is almost as large as the 70  @xmath0 m  image and a `` background '' region is difficult to define .      figures  [ fig : seds ] and [ fig : mips_compare ] clearly demonstrate that only one of the clusters in our sample
ngc  6341has a significant infrared excess . to quantify this
, we computed ( limits on ) the 70  @xmath0 m  non - stellar emission in both the half - light and core apertures . to compute the non - stellar emission in the half - light aperture
, we assumed that any photospheric emission at 70  @xmath0 m  follows the optical emission , i.e. , half of the total emission is contained within @xmath30 . as an estimate of the stellar emission at 70  @xmath0 m  we use the predicted value from the black - body fit described above .
thus the non - stellar emission at 70  @xmath0 m  in the half - light radius is @xmath31 the fraction of stellar emission contained within the core aperture is not as straightforward to compute ; to estimate it , we scale the flux measured at 24  @xmath0 m  in the core aperture ( recall that the image used to measure this flux was convolved to the 70  @xmath0 m  resolution ) by the ratios of the predicted total fluxes at 70 and 24  @xmath0 m
. therefore the non - stellar emission at 70  @xmath0 m  in the core aperture is : @xmath32\ ] ] where the factor of 1.76 is the aperture correction for a point source measured in a 16  aperture .
table  [ tab : ns_dust ] gives the derived non - stellar fluxes for each cluster , together with the dust mass predictions and observational limits derived in
[ sec : dust_limit ] .
only one cluster shows even tentative evidence for the detection of excess ir emission at 70  @xmath0 m : ngc  6341 .
figure  [ fig : images70 ] shows that the apparent 70  @xmath0 m  emission in this cluster comes from an extended region offset from the cluster core by roughly 43  arcsec , or about @xmath33 . within the emission region
there are three separate peaks which contain about 60% of the total flux ; it is unclear whether these peaks are separate point sources ( none has a co - located 24  @xmath0 m  source ) or merely represent structure in the emission . in comparison ,
the dust emission detected in ngc  7078 by @xcite is much rounder and is located closer to the center of that cluster ( @xmath34 ) ; however , the possible dust cloud in ngc  5024 detected by @xcite is also far from the cluster center , at @xmath35 .
the iris 60  @xmath0 m  sky maps do show some nebulosity near ngc  6341 , so it is possible that the detected flux is in fact not associated with the cluster ; however ngc  1851 , at a similar galactic latitude , also has some nearby nebulosity and does not have a 70  @xmath0 m  detection .
future data from the akari all - sky survey should have sufficient spatial resolution to address the issue of whether there is diffuse emission near ngc  6341 .
the lack of any excess infrared emission in the 90  @xmath0 m  akari observations of ngc  6341 by @xcite ( see fig .  [ fig : seds ] ) casts doubt on the reality of our low - significance detection ; however , we retain it for analysis in the remainder of this paper .
while only some of the clusters have spatially diffuse ir emission within @xmath30 , all of the 70  @xmath0 m  images except those of ngc  104 and ngc  5904 do contain some faint point sources which are not concentrated toward the cluster centers .
[ most of the point sources identified by @xcite are off the fields of view of our images . ]
what is the nature of these sources ?
the 70  @xmath0 m  data provide only marginal detections and poor spatial localization , but the 24  @xmath0 m  data provide additional clues to the nature of the sources .
initially , visual examination was used to generate lists of possible point sources in each cluster image . using an estimated 70  @xmath0 m  position as a starting point , we performed aperture photometry on the 24 and 70  @xmath0 m  images using iraf / phot , with the photometry routine allowed to re - center the sources .
several 70  @xmath0 m  sources were off the edges of the 24  @xmath0 m  images , or not visible at the shorter wavelength : the latter group could either be very red sources , or noise fluctuations in the 70  @xmath0 m  images .
standard photometric apertures as given in the mips data handbook were used ( at 70  @xmath0 m , aperture radius 18  and background annulus @xmath36 ; at 24  @xmath0 m , aperture radius 6  and background annulus @xmath37 ) and appropriate aperture corrections were applied .
we did not attempt irac photometry for the 70  @xmath0 m  sources : the irac images of the gcs are significantly more crowded , and we were doubtful about the reliability of measurements in these regions .
the 70  @xmath0 m  sources in the globular cluster fields could be stars within the cluster , unrelated galactic stars , or background galaxies . to discriminate between these possibilities ,
a color - magnitude diagram was made from the mips photometry , shown in figure  [ fig : mips_cmd ] .
for comparison , the agb and post - agb models ( restricted to oxygen - rich models with @xmath38 , which are most appropriate for comparison with gcs ) of @xcite are also shown .
the model fluxes have been scaled down by a factor of 5 to reflect a dust - to - gas ratio appropriate for gc metallicities .
the model fluxes are scaled to a common fiducial distance ( 8.5  kpc ) so the same practice was followed for the cluster sources .
in addition to the evolved star models , we show the ( unscaled ) flux densities and ratios for 70  @xmath0 m  sources in the elais - n2 field of the swire legacy survey @xcite ; these sources are predominantly extragalactic .
all but one of the cluster field 70  @xmath0 m  sources are very red , with typical values of @xmath39 , compared to the rayleigh - jeans value of @xmath40 .
the cluster source colors are more similar to those of the galaxies than to most of the evolved - star models , and the stellar models with similar colors are significantly brighter .    from the colors , we tentatively conclude that the 70  @xmath0 m  sources in the cluster fields are likely background galaxies .
the high galactic latitude of the clusters ( @xmath41 ) supports this conclusion , as do the results of @xcite from their study of 70  @xmath0 m  sources near @xmath42  cen and those of @xcite .
the number density of sources is also consistent with this conclusion : there are 23 sources in a total of 346 arcmin@xmath43 for a density of 0.066 arcmin@xmath44 .
the number counts given by @xcite for the goods - north field show that there are 0.07 galaxies  arcmin@xmath44 to flux limits of 11  mjy ( appropriate for our average field ) , comparable to the number in the cluster fields .
@xcite observed ngc  104 and 362 at mid - infrared wavelengths from 2.424  @xmath0 m .
those authors found eight red sources ( @xmath45 ) in ngc  362 ; no such sources were found in ngc  104 . comparing our images of ngc  362 to those of @xcite
, we find that their red sources a , d and f are within our field of view .
only source d is detected at 70  @xmath0 m  ; it has @xmath46 , consistent with the colors of a galaxy .
the non - stellar 70  @xmath0 m  fluxes computed in  [ sec : nsflux ] allow us to estimate the amount of dust present in these clusters . for the clusters with negative non - stellar fluxes , we computed the dust mass upper limits as @xmath47 .
we follow @xcite in assuming any dust to be optically thin , so the dust mass is given by : @xmath48 } \frac{d^2_{\rm kpc}}{{\kappa}_{\nu } b_\nu(t_d)}\ ] ] where @xmath49 is the dust absorption coefficient and @xmath50 the planck function , both in cgs units . our data
do not constrain the dust temperature , so we use the value found by @xcite for ngc  7078 ( @xmath51  k ) , and the same @xmath52  cm@xmath43  g@xmath4 .
different assumptions about these two quantities can change the derived dust masses by several orders of magnitude ; for example , @xcite estimated the dust mass in ngc  7078 by using @xmath54 for amorphous fayalite at 50  k ; at 70  @xmath0 m  this substance has @xmath55
cm@xmath43  g@xmath4 , or five times the value used by @xcite . on the opposite extreme ,
@xcite assume silicate grains with @xmath56  k ; with @xmath57 , the absorption coefficient at 70  @xmath0 m  is 36  cm@xmath43  g@xmath4 .
table  5 of @xcite shows the effects of different assumptions about the dust temperature on derived dust masses : changing from @xmath51  k to @xmath58  k , for the same composition , increases the inferred @xmath59 by about a factor of ten , while increasing from @xmath51  k to @xmath60  k decreases @xmath59 by about a factor of 2.5 .
uncertainties in the dust composition and temperature therefore result in uncertainties of @xmath61 orders of magnitude in derived dust masses .
table  [ tab : ns_dust ] gives the dust mass upper limits derived from the 70  @xmath0 m  measurements in both the half - light and core apertures .
all values reported are @xmath62 upper limits , except for the possible ( @xmath63 ) detection of total dust mass in ngc  6341 .
this is our only observation resulting in even a tentative detection ; we show it as a value rather than an upper limit in the table and following plots in order to give a visual estimate of the uncertainties .
because the cluster emission is generally very faint , we have not attempted to subtract any point source emission from the total fluxes ; if point sources within the apertures are unrelated to the clusters then the true upper limits would be even lower . for all of the clusters in our sample ,
the limits on the total dust masses are @xmath64 ; the limits on core masses are an order of magnitude lower , @xmath65 .
@xcite gave upper limits to the dust masses for ngc  104 and ngc  5272 ; our total value for ngc  104 is slightly above their limits while our core value is slightly below .
our total value for ngc  5272 is comparable to that of @xcite for 70  k dust .
our value for ngc  5272 is also below that given by @xcite for 80  k dust based on their observations of this cluster at millimeter wavelengths .
our upper limits are above those of @xcite for the five clusters in common , due to the greater sensitivity of their observations .
the dust mass limits we infer are consistent with measurements and upper limits on gas masses . converting dust masses to gas masses by multiplying by @xmath66}}$ ] , we derive gas masses or upper limits for our sample clusters in the ranges @xmath67 ( core ) and @xmath68 ( total ) .
the only confirmed h  i detection , @xmath69 in ngc  7078 @xcite , is nicely within this range .
our inferred limits on the gas mass are consistent with the results of previous searches for intra - cluster neutral and ionized hydrogen , as summarized in @xcite .
two clusters have reported detections of ionized gas : ngc  104 @xcite and ngc  1851 @xcite .
@xcite derived a total gas mass of @xmath70 and stated that they expected essentially all the gas in the cluster to be ionized ; @xcite estimated the total mass of ionized hydrogen in ngc  1851 as @xmath71 .
our gas mass upper limit for ngc  104 is about the same as the @xcite detection while our upper limit for ngc  1851 is above the @xcite level .      to understand the significance of the dust masses and upper limits derived from the _ spitzer _ observations , we need theoretical estimates of the expected dust content of the clusters in the sample .
various estimates of the dust mass within a globular cluster are possible , but little detailed modeling has been done with the exception of the work of @xcite .
those authors combined dynamical models of star clusters with radiative transfer models to predict dust grain temperatures and luminosities for a few individual clusters , given the presence of @xmath72 of dust . while it would be valuable to repeat such modeling with modern radiative transfer codes and understanding of dust properties , for present purposes a simpler estimate of globular cluster dust content will suffice .
following @xcite , the dust mass may be estimated by : @xmath73}}{100 } \label{eq : dustpred}\ ] ] where @xmath74  yr is the horizontal branch lifetime , @xmath75 is the time since last crossing of the galactic plane , @xmath76 is the number of stars on the horizontal branch , @xmath77 is the gas mass lost from each star , and @xmath78}$ ] is the cluster metallicity .
the last factor scales for the change in gas - to - dust ratio with metallicity .
values for @xmath75 are listed in table  [ tab : targets ] ; for the two clusters not listed by @xcite we assume a lower limit of @xmath79  yr . to estimate @xmath76 we use the ` specific evolutionary flux ' method of @xcite
( see also * ? ? ?
* ) : @xmath80 assuming @xmath81  stars  yr@xmath4  l@xmath82 as appropriate for globular clusters , and @xmath83 @xcite using @xmath84 and @xmath85 given by @xcite .
for one cluster , we can compare this dust mass estimate to the current dust production rate : @xcite have recently used _
data to identify the mass - losing stars in ngc  104 . summing their mass loss rates for the brightest stars ( @xmath86 ) and converting to a dust production rate with a gas - to - dust ratio of 200
, we estimate @xmath87  yr@xmath4 .
when multiplied by @xmath88  yr , this yields a total dust content of @xmath89 , about twice that predicted from eq .
[ eq : dustpred ] .
the observations of @xcite are likely to be affected by blending ; correcting for this would lower the inferred dust mass to a value closer to the prediction .
even without such a correction , the comparison gives us confidence that the dust mass predictions are of the correct order of magnitude .
figure  [ fig : dustpred ] compares the predicted and observed dust masses and upper limits . to augment our sample we analyse additional clusters observed by other groups : @xcite , @xcite , and ngc  7078 @xcite , and the six other clusters observed by @xcite : , , , , and .
where several values for dust mass or upper limits are given in published works , only the one with assumed dust temperature closest to 70  k is used .
all of our measurements are well below the predicted values ; for most clusters the disagreement is more than an order of magnitude .
the same is true for clusters observed by other groups , in that only two clusters have the inferred dust mass close to the predicted value : ngc  6356 @xcite ( only for specific assumptions about the dust distribution ) , and ngc  5634 @xcite , one of the lowest - mass , lowest - metallicity clusters shown .
while assuming a lower dust temperature or higher absorption coefficient would increase the dust masses inferred from the observations , it seems unlikely that this can explain the three order of magnitude discrepancy for clusters such as ngc  5272 .
the other explanation is that one or more of the assumptions involved in eq .
[ eq : dustpred ] is incorrect : what if globular cluster stars lose much less than @xmath90 on average , the number of horizontal branch stars does not scale with @xmath91 , the dust - to - gas ratio does not scale with metallicity , or the dust is removed from clusters on a timescale much shorter than the disk passage time @xmath75 ? direct observations of red giant mass loss and color - magnitude diagram modeling are consistent with mass losses of approximately @xmath90 @xcite , and the scaling of @xmath76 with @xmath91 is well - grounded in stellar evolutionary theory .
we address the remaining factors below .    with a small sample of clusters , and an even smaller number of low - significance detections ,
deriving conclusions about the effects of gc properties on mass loss is difficult .
figure  [ fig : masslim ] shows the dust masses and limits plotted as a function of @xmath78}$ ] , @xmath91 , @xmath75 , and @xmath92 .
( in this and the following figure , we plot the more stringent dust mass limits of @xcite for the five clusters in common with our sample . )
there are no clear correlations : while the clusters with large inferred dust masses tend to have large @xmath75 , consistent with the expectation that intracluster dust builds up between cluster passage through the disk , other clusters with similar values of @xmath75 have much smaller inferred masses .
this agrees with the finding of @xcite that the icm dust in ngc  7078 must have been collecting for a time @xmath93 , and with the computation by @xcite that cluster dust lifetimes @xmath94 .
the clusters with dust detections do not have the highest central escape velocities , suggesting that possession of a deep potential well is not in itself enough to guarantee the presence of an icm .
the clusters cover a relatively small range in bolometric luminosity ; this parameter does not seem to have a strong effect on inferred dust mass .
how does metallicity relate to intracluster dust ? the two clusters with the largest inferred dust masses ( ngc  7078 and ngc  6356 ) have the lowest and highest metallicities in the sample .
ngc  5024 , with a tentative dust detection , also has a low metallicity .
the globular cluster data do not help to clarify the controversy over the connection between metallicity and dust production : while some studies show that metallicity and mass loss are not correlated in carbon - rich agb stars in the magellanic clouds @xcite , others find that metal - poor carbon stars , like oxygen - rich giants , form less dust @xcite .
the clusters considered here are almost all more metal - poor than the magellanic clouds , so their agb stars are expected to be oxygen , rather than carbon - rich ; however the effects ( if any ) of metallicity on dust mass loss are not evident from the present data .
observations of mass - losing stars within gcs show that mass loss is occurring at roughly the expected rate , so the non - detection of the icm means that it must somehow be removed .
one possible mechanism is ram - pressure stripping of the cluster icm by the galactic halo .
support for this mechanism is provided by detections of diffuse x  ray emission near several galactic globular clusters @xcite ; @xcite show that the temperature and luminosity of the diffuse emission is consistent with it arising from stripped intracluster gas , and that the energetics require nearly all of the cluster gas to be removed .
those authors also state , to show diffuse x  ray emission , a cluster should have a velocity relative to the galactic halo @xmath95  km  s@xmath4 .
we computed @xmath96 for our sample clusters from the space velocities given in @xcite , @xcite , or @xcite assuming a non - rotating halo ; these are plotted against dust masses and limits in the left panel of figure  [ fig : dust_remove ] .
the high velocities of both ngc  7078 and ngc  5024 appear to be problematic for halo stripping ; a measurement of the space velocity of ngc  6356 might help to clarify the situation .
there is certainly no evidence that clusters with low @xmath96 have more dust ; however , the correlation of cluster kinematics with metallicity combined with any dependence on metallicity of the dust mass limits may also confuse the issue .
another possible icm removal mechanism is kinetic energy injected by stellar collisions , as proposed by @xcite for m15 .
those authors use an n - body model to estimate the time between collisions in the cluster core ; their result is of the same order of magnitude as a less - sophisticated estimate based on eq .
8 - 125 of @xcite : @xmath97 \frac{\sigma_c}{\rho_0 ^ 2 r_c^3}\ ] ] where @xmath98 is the number of stars within the cluster core ( of radius @xmath99  pc ) , @xmath100 is the central velocity dispersion in km  s@xmath4 , and @xmath101 is the central mass density in @xmath102  pc@xmath103 . using the values for @xmath99 , @xmath101 , and @xmath104 , from from king - model fits given in @xcite , we computed @xmath105 for each cluster in the sample ; it is plotted against @xmath59 in the right panel of fig .
[ fig : dust_remove ] .
again , there is no clear correlation of the time between collisions and dust masses : while both ngc  6356 and ngc  5024 have long @xmath106 , so do many other clusters .
ngc  7078 is core - collapsed and should have a short central collision time . we conclude that neither ram pressure stripping and stellar collisions as the cause of icm removal can be ruled out by the current data .
the intracluster medium of globular clusters remains elusive .
our sample comprised eight galactic globular clusters thought to have the best chance for retaining an intracluster medium .
no significant extended emission is detected from seven of the eight target clusters : any such emission may well be hidden underneath confusion noise from galactic cirrus .
one cluster , ngc  6341 , has tentative evidence for excess ir emission from cool dust , although the lack of a 90  @xmath0 m  excess casts doubt on the detection .
all of our inferred dust masses and upper limits are @xmath107 , compared to predicted dust masses ranging from @xmath108 for the same clusters .
searches for intracluster gas have been similarly unsuccessful ( e.g. , * ? ? ?
the lack of intracluster dust and gas implies that either ( 1 ) evolved stars in globular clusters make less dust than predicted , or ( 2 ) the icm escapes from the clusters more quickly than expected .
the first possibility seems unlikely because of the evidence to the contrary in both ngc  7078 and ngc  5139 , where two of the three dustiest , most mass - losing stars have @xmath78 } < -2.0 $ ] @xcite .
it is important to point out that icm escape is not just a trivial occurrence in the life of a globular cluster : substantial mass loss over a gc s history could affect its dynamical evolution and perhaps even its survival as a bound system @xcite .
several mechanisms for icm escape have been discussed above , with energy injected by stellar collisions perhaps being favored over ram pressure from galactic halo gas .
stellar winds provide another possible mechanism ; giant stars with wind velocities above the escape velocity of a typical cluster have been observed @xcite . while there is observational evidence that metallicity does not strongly affect the mass loss rate of individual stars @xcite , stellar wind velocities do appear to scale with metallicity @xcite .
faster winds in agb stars in metal - rich clusters could therefore increase the rate at which the icm is driven out of the cluster .
the existing small sample of clusters does not permit us to distinguish between the possibilities for icm removal .
decades of searching have yielded few detections of the elusive globular cluster icm . a possible strategy for future observations hoping to detect
the icm is to study the removal mechanisms in detail and identify those clusters most likely to retain some icm
. observations of distant clusters which have had few interactions with the milky way could place limits on the importance of halo stripping , although the larger distance of such clusters reduces the effective sensitivity . while infrared facilities available in the near future will not greatly improve on the sensitivity of _ spitzer
_ , observing a larger sample of clusters would allow a more thorough exploration of the dust content as a function of cluster properties .
observations with future radio facilities such as the square kilometer array will place more stringent limits on the gas content . together with more sophisticated modeling efforts and detailed examination of mass loss in individual cluster stars , a multi - wavelength observational approach provides the best chance for finding and explaining the elusive intracluster medium in globular clusters .      we thank the anonymous referee for a careful reading of the manuscript and helpful comments and r.t .
rood for early access to the irac data on ngc  6205 and ngc  6341 .
this work is based on observations made with the _
spitzer space telescope _ , which is operated by the jet propulsion laboratory , california institute of technology under nasa contract 1407 .
p.  b. was partially supported by a discovery grant from the natural sciences and engineering research council of canada .
m.  l.  b. , c.  e.  w. , e.  p. , and r.  d.  g. are supported in part by nasa through _ spitzer _ contracts 1276760 , 1256406 , and 1215746 issued by jpl / caltech to the university of minnesota .
lrlrrrr ngc 104 ( 47 tuc ) & 4.5 & @xmath109 & 65.9 & 190.1 & @xmath110&@xmath111 + ngc 362 & 12.1 & @xmath112 & 46.7 & 49.1 & @xmath113&@xmath114 + ngc 1851 & 8.5 & @xmath115 & 48.7 & 30.4 & @xmath116 & + ngc 5272 ( m3 ) & 10.4 & @xmath117 & 27.6 & 138.4 & @xmath118&@xmath119 + ngc 5904 ( m5 ) & 7.5 & @xmath120 & 37.2 & 105.9&@xmath121&@xmath119 + ngc 6205 ( m13 ) & 7.7 & @xmath122 & 32.1 & 101.4 & @xmath123&@xmath124 + ngc 6341 ( m92 ) & 8.2 & @xmath125 & 35.9 & 61.4 & @xmath126&@xmath111 + ngc 6752 & 4.0 & @xmath127 & & 114.8 & @xmath128 & +    lrrrrl ngc 104 & 107 & 67.9 & 46.4 & 24.7 & , + ngc 362 & 7.10 & 4.82 & 3.44 & 1.73 & + ngc 1851 & 3.57 & 2.31 & 1.64 & 0.82 & + ngc 5272 & 8.67 & 5.68 & 3.63 & 1.85 & + ngc 5904 & 9.25 & 6.56 & 4.81 & 2.42 & + ngc 6205 & 9.66 & 6.43 & 4.35 & 2.17 & + ngc 6341 & 5.05 & 3.38 & 2.31 & 1.16 & + ngc 6752 & 13.25 & 8.41 & 6.32 & 3.42 & +    llrrrrl ngc 104 & 00:24:07.4 @xmath12972:04:48&@xmath130 & @xmath131&@xmath132 & @xmath133 & , + ngc 362 & 01:03:14.8 @xmath12970:50:51&@xmath134 & @xmath135 & @xmath136 & @xmath137 & , + ngc 1851 & 05:14:06.8 @xmath12940:02:49&@xmath138 & @xmath139 & @xmath140 & @xmath141 & , + ngc 5272 & 13:42:11.6 + 28:22:47&@xmath142 & @xmath143 & @xmath144 & @xmath145 & , + ngc 5904 & 15:18:33.0 + 02:04:50&@xmath146 & @xmath147 & @xmath148 & @xmath149 & , + ngc 6205 & 16:41:39.9 + 36:27:28&@xmath150 & @xmath151 & @xmath152 & @xmath153 & , + ngc 6341 & 17:17:07.8 + 43:08:11&@xmath154 & @xmath155 & @xmath156 & @xmath157 & , + ngc 6752 & 19:10:52.0 @xmath12959:59:02&@xmath158 & @xmath159 & @xmath156 & @xmath160 & , +    lrrrrr ngc 104 & @xmath161 & @xmath162 & @xmath163&@xmath164&@xmath165 + ngc 362 & @xmath166 & @xmath167 & @xmath168&@xmath169 & @xmath170 + ngc 1851 & @xmath171 & @xmath172 & @xmath173&@xmath174 & @xmath175 + ngc 5272 & @xmath176 & @xmath177 & @xmath178&@xmath179 & @xmath180 + ngc 5904 & @xmath181 & @xmath182&@xmath183&@xmath184 & @xmath185 + ngc 6205 & @xmath186 & @xmath187 & @xmath188&@xmath189 & @xmath190 + ngc 6341 & @xmath191 & @xmath192 & @xmath193&@xmath194 & @xmath195 + ngc 6752 & @xmath196 & @xmath197&@xmath198&@xmath199 & @xmath200 +