cosmological gamma - ray bursts are one of the most interesting and exotic phenomena in astrophysics . in standard burst models ( e.g. see piran 1999 ; mszros 2001 for reviews ) , the rapidly expanding fireball decelerates , converting the internal energy of the hot plasma into kinetic energy of the beamed , relativistically moving ejecta and electron - positron pairs .
at the point where the fireball becomes optically thin and the gamma - ray burst ( grb ) we see is emitted , the matter will not emit non - thermal gamma - rays unless some mechanism can efficiently re - convert the kinetic energy back into internal energy , i.e. , unless some particle acceleration process takes place .
diffusive shock acceleration is widely believed to be this mechanism ( e.g. , rees & mszros 1992 ; piran 1999 ) for the prompt emission ; it is also likely to be a key element for the underlying physics of x - ray and optical flashes and burst afterglows .
an important question concerning burst models is whether or not their prompt emission can be accurately described by a detailed shock acceleration analysis , and if so , how is the pertinent parameter space of models limited by the observations ?
tavani ( 1996a , b ) proposed a synchrotron shock emission model that invoked such diffusive acceleration as the means for generating non - thermal distributions .
this work has been an interpretative driver in the field , due in part to the impressive model fits of data he obtained .
yet , tavani s model and many subsequent works do not treat the critical involvement the acceleration process has on shaping the relationship between the thermal and non - thermal portions of the electron distribution .
the normal approach is to simplify the distribution to a truncated or broken power - law , occasionally with an additional thermal component .
the viability of the synchrotron model , or emission scenarios that invoke any other radiation mechanism , hinges on such subtleties of the particle distribution , which are inextricably linked to the nature of shock acceleration , diffusive or coherent .
in addition , they must provide compatibility with the hard x - ray spectral index , an issue raised by preece et al .
( 1998 ; 2000 ) , who observed that around 1/3 of all batse bursts exhibited spectra that were too flat below the peak to accommodate a synchrotron interpretation .    this paper explores the general expectations for burst emission spectra from shock acceleration theory .
parametric descriptions for a linear combination of thermal and non - thermal distributions are used to obtain fits to burst spectra for the cases of the synchrotron , inverse compton and synchrotron self - compton ( ssc ) mechanisms , and these distributions are interpreted in the context of known results from acceleration simulations .
cases of particles cooling much slower than their acceleration are treated in this exposition ; justification for these are discussed in section  [ sec : cooling ] .
the focus is on grbs observed by the compton gamma - ray observatory , specifically those with egret detections that provide broad spectral coverage and therefore more constrained fits around the spectral peak .
it is found that acceptable fits are only possible with a marked dominance of non - thermal electrons , contrasting the particle distributions obtained in acceleration simulations .
this poses a general problem for any burst acceleration paradigm where the non - thermal population is drawn probabilistically from a thermal gas .
in addition , fits are found to be more or less equally acceptable for the synchrotron and inverse compton processes ; it is anticipated that swift will provide a number of bright bursts where spectral fitting down into the x - ray band might enable discrimination between such emission mechanisms . however , at this juncture , spectroscopic discrimination against unabsorbed synchrotron self - compton scenarios seems probable , given that the characteristically broad ssc spectra in space are difficult to reconcile with cgro data . in service of this analysis ,
a compact formalism for ssc emissivities from thermal and truncated power - law electron distributions is present in the appendix .
it is anticipated that glast s broad spectral coverage from kev to gev , supplied by the gamma - ray burst monitor in the batse band and the large area telescope in and above the egret band , will render glast a powerful tool in constraining and interpreting burst spectral properties in the future .
this section explores the relationship of synchrotron spectra to the electron energy and angular distributions , inferring lorentz factor distributions using data from the bright burst grb 910503 as a representative case study .
the formalism for synchrotron radiation in optically thin environs is standard and readily available in many textbooks , including the treatments of bekefi ( 1966 ) , jackson ( 1975 ) and rybicki & lightman ( 1979 ) . throughout this paper ,
the standard convention for the labelling of the photon polarizations will be adopted , namely that refers to the state with the photon s _ electric _ field vector parallel to the plane containing the magnetic field and the photon s momentum vector , while denotes the photon s electric field vector being normal to this plane .
the classical ( photon ) angle - integrated emissivities for the and polarization for monoenergetic electrons of lorentz factor moving at angle to the field are ( e.g. westfold 1959 ; jackson 1975 ; rybicki and lightman 1979 ) @xmath0 in units of , where the are modified bessel functions of the second kind .
the generally higher probability of emitting photons than ones reflects the intrinsic dipolar radiative nature of accelerating / oscillating charges ; the polarization is elliptical for photon propagation along the field .
the characteristic synchrotron rate factor that scales these emissivities is @xmath1 here , is the fine structure constant .
the photon energies that appear in these functions that form the scaling for the energies of emission are [ with the electron compton wavelength over being ] : @xmath2 where is the particle s larmor radius , and gauss is the quantum critical field .
note that since , the photon angles to the field form a narrow distribution about .
furthermore , the expressions in eq .
( [ eq : synchrate ] ) are appropriate for cases where the particle momenta are not closely aligned with the magnetic field ; situations of such near alignment are discussed in section  [ sec : pitch_synch ] below .
the total synchrotron emissivity is obtained by integrating over the energy and angular distribution of the particles . for electron distributions in units of , summing over polarizations forms the expression for the emissivity that will be used mostly throughout this paper ( in units of ): @xmath3
this is the differential photon spectrum ; frequently in this work , reference will be made to the so - called spectrum , a spectral energy distribution that is proportional to .
physical diagnostics on the grb environment are obtained by comparing the emission spectrum in eq .
( [ eq : synch_emiss ] ) with observations .
the obvious database of choice comprises batse observations from the compton gamma - ray observatory ( cgro ) , specifically the batse spectroscopy catalog of preece et al .
( 2000 ) with its pubicly - accessible electronic compilation .
yet this data possesses a principal limitation , namely that imposed by the restricted energy band of batse .
inferences of particle distributions are generally on a more secure footing when using broader burst spectra , particularly in domains of significant spectra structure , which is the intrinsic nature of the batse pool .
thus , the two decades of energy afforded by batse often are insufficient to describe the nature of the non - thermal particle distribution invoked in models .
data from the egret experiment on cgro suitably augment this situation and extend the dynamic range by 24 orders of magnitude in energy .
accordingly , a focus on egret bursts is insightful and offers more secure determinations of particle distributions .
use is therefore made of spectral compilations of batse , comptel and egret data that are offered in schaefer et al .
specifically , a case study is made here of grb 910503 , observing that the inferences to be made generally extend to the handful of other egret bursts , and also to a large number of bright batse bursts that possess no confirmed egret detections .
note that the focus here is primarily on time - integrated spectra .
given the complexity and diversity of burst time profiles , and associated variations in spectral shape and the well - known hard - to - soft evolution , detailed spectral fits will depend on the time window chosen . a classic example for this
is provided by the gradual emergence of a hard gamma - ray component in the spectrum of grb 941017 ( gonzalez , et al .
the main conclusions of this paper apply to the mean properties of the emitting particle distributions , and are striking enough to indicate their applicability for most or nearly all of the duration of the bursts studied herein .
a variety of particle distributions are possible in modeling burst spectra .
clearly , all should contain a non - thermal component , since the data are clearly non - thermal : purely isothermal models do not work ( pacynski 1986 ; goodman 1986 ) , though constructed convolutions of thermal spectra of different temperatures are possible .
the emphasis here is motivated by physical interpretation , and the leading contender for non - thermal particle energization in bursts is diffusive acceleration at relativistic shocks in the burst outflow .
such acceleration is well studied in non - relativistic flows , and is widely believed to originate with thermal particles that are diffusively transported when interacting with field turbulence in the shock environs .
such transport effects many crossings of the shock layer ( e.g. see drury 1983 ; blandford & eichler 1987 ; jones & ellison 1991 for reviews ) , that lead to a friction term in the momentum evolution , i.e. acceleration .
while other acceleration models such as reconnection and coherent electrodynamic energization exist and are appropriate in a variety of astrophysical environments ( e.g. the solar corona , pulsar magnetospheres ) , diffusive shock acceleration is the most widely invoked in astrophysical models ; it forms the centerpiece of discussion here as well as for many investigations of bursts .
tavani ( 1996a , b ) proposed a shock emission model that invoked such diffusive acceleration as the means for generating non - thermal distributions .
this work has been an interpretative driver in the field , in no small part due to the impressive model fits of data obtained .
the principal purpose of this section is to delve deeper into the physical implications of such data fitting .
accordingly , it is appropriate to use an electron distribution similar to tavani s ( 1996a , b ) to add insight to the discussion . here , the quasi - isotropic electron distribution @xmath4 is adopted , where is a step function with for and zero otherwise , and is the system or post - shock temperature .
the factor directly relates to the efficiency of acceleration , and is a parameter in the fitting protocol , as are ( assumed much greater then unity ) , the dimensionless , the non - thermal index , and the normalization factor .
the ensuing discussion will provide physical grounds justifying this choice of distribution in shocked environs ; other possibilities with different mandates include the broken power - law adopted by lloyd & petrosian ( 2000 ) .
note that is independent of for true isotropy . in many cases in this paper , will be chosen for simplicity , since this corresponds to the dominant radiative contribution from the synchrotron mechanism .
the parameter is introduced to generalize from tavani s ( 1996a , b ) fits , permitting the non - thermal distribution to emerge perhaps at truly super - thermal energies , rather than exactly in the thermal peak ( tavani s restricted case ) .
typical acceleration spectra in non - relativistic shocks assume forms similar in general appearance to eq .
( [ eq : elec_dist ] ) . specifically , they correspond to the portion of parameter space ( at least for ions ) , so that dissipation in the shock layer leads to a heating plus non - thermal acceleration that ensues only at super - thermal energies in the maxwell - boltzmann tail .
evidence for such a structured distribution can be found both in heliospheric measurements of accelerated ions near shocks ( e.g. earth s bow shock : see ellison , mbius & paschmann 1990 ; interplanetary shocks : see baring et al .
1997 ) , and also a variety of acceleration simulations and theoretical analyses ( e.g. see scholer , trattner & kucharek 1992 ; giacalone et al
. 1993 , baring , ellison & jones 1993 ; ellison , baring & jones 1996 ; kang & jones 1997 ) .
the situation for relativistic shocks is less certain , due in no small part to the absence of _ in situ _
spacecraft determinations of particle spectra ; simulation expectations for such environments are discussed below . in this paper , values of appear ,
being in general only weakly constrained by the fits and largely motivated by the above physical and theoretical considerations .    while numerical evaluations of eq .
( [ eq : synch_emiss ] ) given the electron distribution in eq .
( [ eq : elec_dist ] ) were performed to obtain approximations to burst spectra , standard asymptotic limits are readily attainable in analytic form .
the non - thermal contribution to the synchrotron spectrum ( from the truncated power - law ) , dominant whenever , assumes the following power - law forms ( derived , for example , from eq .
[ [ eq : synchfinal ] ] ) : @xmath5 here , is the component of the field perpendicular to the instantaneous electron motion .
the thermal component , sampled when , necessarily declines exponentially at photon energies well above the thermal peak , but possesses the same power - law dependence as in the non - thermal case for low photon energies : @xmath6 this spectral similarity at low energies follows from the narrow nature of the distributions being convolved with the synchrotron emissivity functions in eq .
( [ eq : synchrate ] ) .
accordingly , it is not possible to distinguish between emission from thermal or power - law distributions if spectrum measurements exist only at energies below the peak ; some other diagnostic is required .
the exponential decline above the thermal peak possesses a somewhat complex two - domain structure ; the reader is referred to pavlov & golenetskii ( 1986 ) , brainerd & petrosian ( 1987 ) and baring ( 1988a ) for the presentation of appropriate asymptotic forms in cases .
it must be noted that throughout this paper , photon spectra are generated in the shock rest frame ( srf ) and then boosted in energy to the observer s frame ( of ) by a fixed doppler factor , independent of photon angle in the srf .
this approximation , essentially adopted in many models of grb spectra , is expedient but does not accommodate the fact that photon distributions in the srf are inherently anisotropic due both to the instrinsic nature of the synchrotron and inverse compton processes , and also to the fact that relativistic shocks naturally generate anisotropic ions and electrons . notwithstanding , this boost approximation does not dramatically distort the angle - integrated of spectral shape for typical srf angular distributions and suffices to describe the continuum for the purposes of the discussions here .      before presentation of results
, it is appropriate to address the modeling procedure used in this paper .
approximate spectra , as opposed to data fits , were derived visually using both log - log and semi - log representations of spectra to increase confidence in the results .
least squares trials were performed in some cases , but with no addition of insight or perceived accuracy .
this choice of a visual approach was one of expediency , and did not compromise the scientific conclusions presented here ; the essential points are unambiguously made without the need of statistical analysis .
this latitude is afforded by the limited energy range of the data in the batse catalog : slight adjustments in parameters used to fit the extant data for a given burst would almost certainly have been dominated by those needed to accommodate a broader spectral range . for the specific cases of grb 910503 and grb 910601 , where egret data are incorporated in the study ,
the uncertainty in the mev data points is a major contributor to the determination of model parameters . in bursts where no egret or comptel data were used ,
the `` fit '' was actually to the band ( 1993 ) model spectral fit presented in the preece et al .
( 2000 ) spectroscopy catalog , not to the actual data . in such cases ,
the high energy band model index was adopted , and other distribution parameters adjusted . in this sense , there was in implicit second order incorporation of the statistics embodied in the band model fit .
it must be emphasized that there is no uniquely preferred method for statistically fitting public data since it is already a convolution of model and instrumental response . the correct approach ( as in lloyd & petrosian 2000 ) would be to fold spectral models here through the batse or egret detector response function , an enterprise that rapidly becomes more ambitious when different instrumental responses must be blended .
such complexity seems premature for the explorations of this paper , and is deferred until the more refined data of swift and glast becomes available .
results for the spectral approximation of the grb 910503 data are presented in figure  [ fig : synch_spec_dist ] , together with the multi - instrument , time - integrated cgro data as provided by schaefer ( private communication ) prior to the public release in schaefer et al .
we remark that the spectral deconvolutions published in schaefer et al .
( 1998 ) were made independently for each data set from the different cgro instruments , applying an artificial normalization to reconcile offsets in the different data amplitudes .
such a procedure is less accurate than a joint spectral deconvolution of the entire spectral data set ( see for example the analysis of grb 990123 data in briggs et al .
notwithstanding , the major conclusions of this paper are robust and are not affected by such subtleties of data analysis .
figure  [ fig : synch_spec_dist ] displays both synchrotron emission spectra for the case of pitch angles ( the so - called _ quasi - isotropic _ case ) , and the inferred electron distribution that generates the continuum . for the solid curve approximate case ,
the overall normalization is obviously a free parameter and was offset slightly from the data for visual clarity .
the electron distributions were normalized to unit area in space .
the solid curve depicts the case that would be required to explain the photon spectrum .
the non - thermal particle distribution clearly dominates the thermal component in this `` fit . ''
its index differs from the value of 3.8 that tavani ( 1996a ) obtained in his analysis of this burst , a discrepancy that is actually insignificant given the error bars in the egret data .
adjustment of the parameter to improve the appearance of the model approximation around and above the spectral peak led to the assumed value of , accurate to around 15% .
notwithstanding , there appears to be some spectral structure around the peak that can not be explained by the simple particle distribution in eq .
( [ eq : elec_dist ] ) .
the positioning of the peak in energy produced the fit value , where is the lorentz factor associated with the bulk motion of the grb expansion . here
, is interpreted as the magnetic field in the burst expansion frame , i.e. , not in the ism or shock rest frames .
the parameter describing the relative normalization of the non - thermal and thermal components was found to be ( ) , a strikingly large value that is not addressed explicitly in detail by tavani ( 1996a ) .
a discussion of the implications of an inferred particle distribution with such a large follows just below .
a similar property is obtained for synchrotron model fits for batse bursts grb 940619 ( trigger number 3035_6 ) and grb 940817 ( trigger number 3128_5 ) from the preece et al .
( 2000 ) catalog , and also for grb 910601 and grb 910814 whose composite spectra are presented in schaefer et al .
the spectra for bursts grb 930131 and grb 940217 that possess egret detections also suggest a similar result , though the lack of a composite broad - band multi - instrument spectrum for each in the literature precluded a detailed analysis .
the parameter regime of is expected to be pervasive for batse bursts , due to the smoothness of their circumpeak continua . for comparative purposes , another distribution and its resultant synchrotron
is illustrated in figure  [ fig : synch_spec_dist ] as dashed curves .
the distribution was deliberately chosen to possess a somewhat dominant thermal component , motivated by expectations from shock acceleration theory , and clearly yields a spectrum that grossly deviates from the observations both well above and well below the peak .
while these deductions are obtained using a time - integrated spectrum , such a conclusion must apply either for most of the burst duration , or at sub - intervals when the burst was at its brightest .    for the synchrotron mechanism to be a viable explanation of grb spectra ,
the underlying particle distribution must be physically realistic . here , expectations of particle distributions from shock acceleration theory are the predominant focus , partly because this was the preferred acceleration mechanism in tavani ( 1996a , b ) , and largely because this is the underlying energization assumed in most burst models . in a nutshell ,
both theory and observation pertaining to non - relativistic shocks strongly indicate and cases .
note that is essentially the parameter that demarcates the relative normalization of the non - thermal and thermal components in equation  ( [ eq : elec_dist ] ) .
applications of shock acceleration in the heliosphere and at supernova remnant shells more or less mandate these parameter regimes , since the seed pool of particles that can be subjected to acceleration are thermal and moreover plentiful . in heliospheric environs ,
the source is the cold solar wind , and there is little evidence of significant pre - heating that is not connected to the shock system under study .
the _ in situ _ spacecraft measurements exhibit strong heating at shocks and the generation of dominant thermal ions in the downstream regions , corresponding to cases .
data taken at the high sonic mach number earth s bow shock ( e.g. see ellison , mbius & paschmann 1990 ; scholer , trattner & kucharek 1992 ) and low sonic and afvnic mach number interplanetary shocks ( e.g. see gosling et al . 1981 ; baring et al .
1997 ) clearly support this assertion . in supernova remnants ,
direct particle measurements are , or course , impossible , and inferences can only be made from radiative signatures . while electron temperatures and densities can be estimated from x - ray signals and molecular emission in the ir / optical bands , the data can not yet clearly discriminate and situations .    from a theoretical perspective
, diffusive shock acceleration naturally draws particles from the thermal pool via transport in shock - associated turbulence . whether this be first or second order in momentum transport , thermal and non - thermal particles are subjected to similar diffusion , at least as far as ions are concerned , so that circumstances would require virtually complete suppression of thermalization and a 100% efficiency of acceleration out of the thermal pool .
such a scenario would require unusually anomalous spatial diffusion coefficients ( i.e. extremely small ) at thermal momenta , or perhaps strange and non - dissipative electrodynamic properties in the shock layer .
note that is more likely to increase at low energies due to the possibility of wave damping by thermal populations ( e.g. see forman , jokipii & owens 1974 , for the relationship between and wave power spectra within the confines of quasi - linear theory ) . without the assumption of anomalous forms for , monte carlo transport simulations
( e.g. see baring , ellison & jones 1993 ; ellison , baring & jones 1996 ) clearly indicate a downstream thermal population that dominates the extended non - thermal distribution in number , corresponding to .
the `` injection '' of ions from the thermal pool is extremely efficient , but not incredibly so .
hybrid plasma simulations ( e.g. see trattner & scholer 1991 , 1993 ; scholer , trattner & kucharek 1992 ; giacalone et al .
1992 , 1993 ; liewer , goldstein & omidi 1993 ) of wave generation and ion transport near non - relativistic shocks , which treat the electron component as a background fluid and therefore address only some of the electrodynamic aspects of the plasma in the shock layer , produce similar distributions .
these indications are strongly suggestive , but not decisive for the gamma - ray burst problem , since they largely do nt focus on electron properties and equally importantly do not explore acceleration at relativistic shocks . in situ measurements of non - thermal electrons at shocks are limited ; examples include recent geotail observations of interplanetary shocks ( shimada , et al .
1999 ; terasawa et al . , 2001 ) .
accordingly , discussion is dependent on simulational information , as is the case for relativistic shocks , and the literature is sparse on such a subject .
monte carlo simulations of ion acceleration at relativistic shocks ( e.g. ellison , jones & reynolds 1990 ; baring 1999 ; ellison & double 2002 ) generate the same injection properties as do non - relativistic shocks , i.e. .
plasma simulations of relativistic shocks are similarly sparsely studied .
gallant et al . ( 1992 ) and hoshino et al . (
1992 ) presented results from one - dimensional ( 1d ) particle - in - cell ( pic ) simulations of electron - positron and electron - positron - proton plasmas in ultrarelativistic flows with perpendicular shocks , where the mean field was orthogonal to the shock normal .
such simulations treat the full electrodynamics of such plasmas , but are severely limited by cpu constraints .
accordingly , they tend to use unrealistically small proton to electron mass ratios , and model relatively small spatial scales .
these constraints tend to mask the production of some wave modes and limit consideration to particles with energies at most a few to ten times thermal .
the context for these investigations was pulsar wind flows , and they revealed no acceleration in pure plasmas . with the addition of protons in significant abundances , hoshino et al .
( 1992 ) demonstrated that non - thermal positrons ( and not electrons ) could be created , as the ions generated elliptically - polarized magnetosonic modes that resonantly interacted with the positrons .
mace and jones ( 1994 , unpublished ; f. jones , private communication ) explored a 1d pic simulation similar to gallant et al .
( 1992 ) for an ultrarelativistic plasma shock of arbitrary field obliquity , and found that electron / positron acceleration arose in quasi - parallel scenarios , but disappeared when the field obliquity to the shock normal increased above around degrees , concurring with zero acceleration in the earlier results for perpendicular shocks . a scaled version of their parallel shock distribution ( binned in energy ) is depicted in fig .  [ fig : synch_spec_dist ] , and is fairly proximate to the , case .
a notable feature of the simulation data is the very limited range of energies accessible to the pic technique , imposed by the severe cpu limitations for such complex codes .
some interesting recent pic simulations results have been presented by shimada & hoshino ( 2000 ) and hoshino & shimada ( 2002 ) for electron - proton plasmas containing mildly - relativistic shocks , where suprathermal power - law electrons are generated and the distribution closely resembles the pic simulation data depicted in fig .  [
fig : synch_spec_dist ] . yet , as with the mace and jones results , the complete distribution reveals that the power - law emerges only well into the maxwellian tail , corresponding to the situation elicited in the monte carlo simulations .
it is noted in passing that the one - dimensional nature of these simulations potentially imposes an important absence of critical physics to the acceleration problem : jokipii , kta & giacalone ( 1993 ) and jones , jokipii and baring ( 1998 ) demonstrate that particle diffusion across field lines is eliminated when simulations restrict the spatial dimensions to less than three .
the implication of this result for the shocks is a cessation of acceleration when thermal particles can not be transported across fields to return to the shock from the downstream side ; baring , ellison & jones ( 1993 ) demonstrated that for high mach number , non - relativistic ion shocks , a quenching of acceleration through this effect would arise at around , a result that matched the behavior seen in the relativistic pic simulations of mace and jones .
moreover this result is consistent with the absence of acceleration in non - relativistic 1d hybrid plasma simulations applied to the quasi - perpendicular heliospheric termination shock ( e.g. liewer , goldstein & omidi 1993 ; kucharek & scholer 1995 ) , a site that is believed to generate the non - thermal anomalous cosmic rays .
the essential loss of physics in simulations of restricted dimensionality is clearly crucial to the correct modeling and interpretation of particle acceleration in gamma - ray bursts .
the principal conclusion of this discussion is that there is little observational or simulational evidence for quasi - isotropic electron distributions in the environs of shocks that would render a synchrotron spectrum commensurate with the grb910503 data in particular and typical burst spectra in general .
an issue that has already been substantially discussed in the literature is whether the synchrotron model is compatible with spectral slopes below the burst peak . for grb 910503 , fig .
[ fig : synch_spec_dist ] clearly indicates that there is no conflict between the data and the model .
the fixed value of for the low energy index in eqs .
[ eq : ntsynch_limits ] and  [ eq : thsynch_limits ] suggests that this is not the situation for all bursts , a fact that was demonstrated by preece et al .
( 1998 ; see also preece et al .
2000 ; 2002 ) .
synchrotron theory distinctly predicts this canonical index at low energies for any isotropic electron distribution that has a characteristic lorentz factor as its lower scale .
correlations between lorentz factors and pitch angles can alter this property , as discussed in section  [ sec : pitch_synch ] below .
any index greater than this value can not be accommodated by quasi - isotropic distributions ; indices less than this are easily described by a steepening particle distribution above .
preece et al.s ( 2000 ) display of the distribution of low energy indices clearly indicated a conflict , as illustrated in fig .
[ fig : batse_low_index ] .
there the low energy indices for various empirical fitting models used in the batse spectroscopy catalog are binned in a histogram , with the permitted range of clearly marked .
the most common fitting model used in the catalog was the popular band ( 1993 ) grb model for the observed flux in photons/@xmath7/sec : @xmath8 most bursts possess values and . in folding the band model , broken power - law models and others through the batse detector response , preece et al .
( 2000 ) found that the band model performed best in roughly 2/3 of the burst sample .
note that the index recorded in the catalog is the asymptotic low energy index ( preece et al .
2000 ) , and does not exactly match the index local to the lowest energy kev in the batse window .
the inferred values of are subject to instrumental selection effects , particularly when the peak energy is low .
hence , there is the possibility that the distribution might move to slightly higher when the detection band moves to lower energy , the situation for swift .
note that such an effect was observed by preece et al .
( 1998 ) when applying a correction to accommodate the energy of the batse window .
a depiction of the histogram for time - integrated spectra in the batse spectroscopy catalog is given in lloyd & petrosian ( 2000 ) ; it exhibits a similar distribution , as expected .
the quasi - isotropic synchrotron model clearly can not account for nearly 25% of sample burst spectra .
this is true even for pathological choices of electron energy distributions : from the form in eq .
( [ eq : synch_emiss ] ) , it is easily shown that the spectral index is a monotonically increasing function of energy for any electron distribution , due to the broad and smooth nature of the synchrotron emissivities in eq .
( [ eq : synchrate ] ) . while this is disturbing , it may be more of an indication of problems with assumptions underpinning the bound than a suggestion that synchrotron radiation is not responsible for grb prompt emission .
for instance , relaxing an adherence to isotropy or distributions can potentially alter the spectral slope .
it is demonstrated in section  [ sec : anis_synch ] below that for an independent of , the spectrum below the peak can not be flattened beyond cases , and in fact steepens . however , introducing a -dependent angular distribution can suitably ease the bound ; this is exemplified in the focus on pitch angle synchrotron emission in section  [ sec : pitch_synch ] .
an alternative recourse is synchrotron self - absorption , i.e. relaxation of an optically thin assumption .
given the theoretical expectation of extremely flat spectra below the break ( for power - law electrons and a raleigh - jeans tail for thermal : see , e.g. , rybicki & lightman 1979 ) , different groups have invoked such a mechanism to account for cases in the batse catalog ( e.g. crider & liang 1999 ; granot , piran & sari 2000 ) . pushing the absorption break to gamma - ray energies requires extreme physical conditions ( @xmath9 , gauss ) , as noted by lloyd & petrosian ( 2000 ) , and can be inferred from eq .
( [ eq : selfabs_energy ] ) .
these parameter regimes would be more extreme in the absence of the strong doppler boosting associated with grb outflows .
equipartition is usually invoked to justify such high fields , however a widely - accepted and well - understood mechanism for the generation of such fields remains elusive . a convenient resolution of this issue is afforded by the proposition of panaitescu & meszaros ( 2000 ) , namely that the prompt emission is an inverse compton image of a self - absorbed synchrotron spectrum .
the concomitant moving of the absorption break energy to below the optical band lowers the required values of and to less extraordinary ranges .
the x - ray spectral indices in the scenario of panaitescu & meszaros ( 2000 ; see also liang , boettcher & kocevski 2003 ) can comfortably accommodate the domain , as discussed in section  [ sec : alphalow_comp ] .
an additional concern is that the turnover energy ( in the grb expansion frame ) satisfying is a function of a number of environmental parameters .
equation  ( 6.53 ) of rybicki & lightman ( 1979 ) for the optical depth for self - absorption can be manipulated to arrive at the equation for .
this result can be found in many papers ; here eq .
( 7 ) of baring ( 1988b ) is adapted to accommodate the specific form for the non - thermal portion of the distribution in eq .
( [ eq : elec_dist ] ) : @xmath10 for @xmath11 here , is the thomson optical depth of the emission region , and , as before . note that represents the gamma function in eq .
( [ eq : phidef ] ) . the density and size of the absorbing medium
are incorporated in .
these parameters possess implicit dependence on the bulk lorentz factor of the flow , with different evolutionary dependences for adiabatic and radiative expansion cases ( e.g. , see dermer , chiang & bttcher 1999 ) . if eq .
( [ eq : selfabs_energy ] ) is interpreted in the frame of the burst outflow , then is in the euv / soft x - ray band for typical burst lorentz factors . given the number of parameters ( five in all ) contributing to the computation of the self - absorption optical depth ,
it is difficult to comprehend how can be fine - tuned to the batse band in the observer s frame ; such a conspiracy is an issue also for optically thin synchrotron emission , whose peak depends on three parameters , , and .    a strong test of the self - absorption proposition could be offered by swift , which should determine whether or not batse was fully capable of determining the `` asymptotic '' values of , especially for those bursts with peaks around 100 kev or lower .
the gamma - ray burst monitor ( gbm ) on glast will provide an additional broad - band probe .
should the batse distribution be mirrored by the swift and gbm ones , synchrotron self - absorption in the gamma - ray band could then only be operating in a small minority of bursts , rendering it less germane to the spectral shape discussion , except perhaps for synchrotron self - compton models .
the foregoing discussion has focused on the case of synchrotron emission from quasi - isotropic or large pitch angle electron distributions .
distributing the angles of electrons in a non - uniform manner can alter the spectral index of the emitted radiation .
a prime example of this was the flattening incurred in distributing magnetic pair creation turnovers in neutron star environments , probed by baring ( 1990 ) . here , obviously
a distribution skewed towards small pitch angles is a potential candidate for modifying .
such a distribution is not at all contrived : synchrotron cooling can produce it quite naturally ( e.g. see brainerd & lamb 1987 ) .
however , as alluded to above , it is now demonstrated that such distributions , if uncorrelated with lorentz factor , yield only spectral steepening , and so can not resolve the low energy index problem .    the spectral behavior can adequately be demonstrated using the non - thermal portion of eq .
( [ eq : elec_dist ] ) , which for the purposes of this discussion will be written in the form for and is zero otherwise . here . if is defined as the dimensionless cyclotron energy , then the differential rate of production of synchrotron photons in eq .
( [ eq : synch_emiss ] ) can be quickly manipulated by reversing the integrations over and .
the result is @xmath12 where @xmath13 is the synchrotron emissivity scale , and @xmath14 this form for the synchrotron rate will be used later in the discussions on synchrotron self - compton emission . by adding and subtracting suitable forms like eq .
( [ eq : synchfinal ] ) with different choices of , a general expression for the synchrotron emission from a piece - wise continuous power - law distribution of electrons can be obtained . defining the photon spectral index to be ,
it is easily shown that is a monotonically increasing function of energy for any electron distribution with a low cutoff , using piecewise continuous power - laws to successively refine approximations to a smooth distribution .
( [ eq : synchfinal ] ) possesses two readily identifiable asymptotic limits , namely when and .
these lead to the forms given in eq .
( [ eq : ntsynch_limits ] ) . now consider an angular distribution defined over the range for .
the superposition of spectral forms in eq .
( [ eq : ntsynch_limits ] ) then immediately indicates that such an angular distribution must generate a spectrum with asymptotic forms at , and when . the spectrum in between depends on how rapidly varies with angle , and can be either of slope or steeper , but never flatter than and never steeper than .
this follows from the convolution over angles of spectra that are monotonically declining in energy necessarily being also monotonically declining in .
this assertion holds even for that is strongly peaked towards , and monotonically declining with .
hence , angular distributions with uncorrelated with can not ameliorate the low energy index issue due to the spectral shape implied by eq .
( [ eq : synchfinal ] ) .
a more viable possibility for flattening the low energy index is when there is a correlation between the pitch angle and the lorentz factor of charged particles .
this provides a distinct kinematic difference from the foregoing considerations of synchrotron radiation from quasi - isotropic electron distributions .
the archetypical example for consideration here is so - called _ pitch - angle synchrotron _ radiation , a mechanism explored in detail by epstein ( 1973 ; for applications , see epstein & petrosian 1973 ) .
this mechanism corresponds to synchrotron emission from particles possessing very small pitch angles .
the polarization - averaged emissivity derived by epstein ( 1973 ) in this limit , for electrons in a uniform magnetic field , can be written in the form @xmath15 this form is simply obtained from the flux emissivity in equation  ( 10 ) of epstein ( 1973 ) , and is applicable for ; at higher energies , the usual synchrotron formula in eq .
( [ eq : synch_emiss ] ) is operable if .
note that this process exhibits strong circular polarization due to the near coalignment of particles with the field , though there is significant linear polarization in narrow bands of photon energy .
observe that defines the intrinsic energy scale of pitch angle synchrotron ( pas ) emission .
( [ eq : synch_pas_rate ] ) is virtually independent of for domains , thereby defining an extremely flat spectrum .
this property is a direct consequence of the kinematic coupling , yielding a dependence that contrasts the isotropic synchrotron dependence .
this weaker dependence on immediately implies an inference of stronger source fields in order to move the characteristic pas photon energy into the batse window .
the spectrum drops dramatically above if , otherwise it merges into the synchrotron form if .
these shapes imply ( i ) that prompt grb spectra that are flatter than below the peak ( i.e. around 98% of bursts : see figure  [ fig : batse_low_index ] ) can be easily accommodated by a pitch angle synchrotron scenario ( just like the inverse compton scattering case considered in section  [ sec : icformalism ] below ) , and ( ii ) accordingly that a variety of grb spectral shapes can be accurately modeled due to the characteristically narrow peaking of the pas emissivity in eq .
( [ eq : synch_pas_rate ] ) . these were conclusions of lloyd - ronning and petrosian ( 2002 ) , who explored fits to batse bursts for pitch angle synchrotron models , and addressed correlations between spectral parameters and temporal properties .
hence , pitch angle synchrotron is an attractively versatile spectroscopic possibility for bursts .
note that the mechanism of `` jitter '' radiation in turbulent magnetic fields ( medvedev 2000 ) possesses similar spectral properties due to the intrinsic kinematic coupling , thereby offering an alternative possibility of similar promise .
pitch angle synchrotron might arise when electrons stochastically increase their pitch angle from virtually zero during an acceleration process . for it to provide the major contribution ,
a prevalence of small pitch angles must be strongly favored , i.e. the pitch angles must not be significantly populated , so that normal synchrotron emission does not dominate .
it is not clear how this might come about , though rapid cooling and inefficient diffusion transverse to the field is an interesting possibility ( e.g. see lloyd - ronning and petrosian 2002 ) .
in addition , small pitch angles must be prepared in the first instance , which suggests a coherent electrodynamic mode of acceleration as opposed to diffusive one .
while pas may afford a means of enhancing the non - thermal radiative signal relative to the thermal component , a detailed exploration of pas spectral shapes and polarization for thermal and non - thermal distributions in the context of diffusive acceleration will be the subject of future work .      while detailed exploration of cooling effects is beyond the scope of this paper ,
a brief discussion is appropriate .
cooling is known to strongly influence spectral shape in quasi - isotropic synchrotron scenarios ( e.g. ghisellini , celotti & lazzati 2000 ) , introducing a second `` line of death '' with a value of .
this presents problems for models invoking strong cooling in generating spectra like those seen by cgro .
such models are most applicable to grb environments where the emission regions are somewhat remote from the acceleration ones .
this is because even though synchrotron radiation is a remarkably efficient process , diffusive acceleration at shocks is even more so .
accordingly , coincident acceleration and emission zones of the same physical scale can exhibit prolific acceleration up to some maximum lorentz factor , before radiative cooling in the magnetic field causes cessation of acceleration .
the spectral turnover that marks this cessation is quasi - exponential in nature , and have not been observed in burst spectra to date .
this combination acceleration / radiation scenario is widely envisaged in the contexts of x - ray and gamma - ray emission from supernova remnants and blazars .
whenever multiple spatial or temporal scales are sampled , convolution of spectral contributions can blur the turnover into a broad distribution , essentially the case explored by ghisellini , celotti & lazzati ( 2000 ) .
the rate of acceleration scales as the inverse gyrofrequency , which has a weaker dependence on electron lorentz factor than the synchrotron cooling rate .
specifically , the acceleration timescale ( in seconds ) can be approximated by , where is the energy of the electrons in units of tev , is the shock speed in units of , is the mean magnetic field strength ( measured in gauss ) , and is a factor of order unity that accounts for corrections to the diffusive timescale appropriate for relativistic shocks . the synchrotron cooling timescale ( also in seconds ) is .
various versions of these formulae are widely available , though a specific discussion of them in the context of blazars is in baring ( 2002 ) ; they are effectively applicable in the comoving frame of the burst outflow , as opposed to the observer s frame .
the electron energy and the magnetic field are coupled by placing the synchrotron peak in the soft gamma - ray band , i.e. setting gauss . here
is the bulk lorentz factor of the outflow .
these quickly lead to the ratio estimate @xmath16 acceleration is clearly much more rapid than synchrotron cooling for typical , a property that changes only when considering electrons that emit above about 10 gev , i.e. for gauss . in the special case of , the result in eq .
( [ eq : cooling_ratio ] ) corresponds to the well - known synchrotron cooling cutoff property of the crab pulsar wind nebula ( e.g. see de jager et al . , 1996 ) ,
namely that the spectral turnover occurs at around 50 mev , independent of the magnetic field strength .
since significant questions surround the viability of the synchrotron mechanism in explaining burst spectra , consideration of alternative emission processes is warranted . here
, the focus is on inverse compton scattering , another very efficient radiative mechanism . for the purposes of the discussion here
, it comes in two relevant varieties : scattering off a narrow - band soft photon distribution , and upscattering of synchrotron photons ( synchrotron self - compton , or ssc ) that provide an inherently broad seed spectrum .
it shall be seen that the resulting predictions for prompt emission spectra for these two cases are manifestly different .
the inverse compton ( ic ) mechanism possesses two obvious advantages over synchrotron radiation : ( i ) high values of can be avoided , since there is no requirement that the synchrotron peak appear in the gamma - ray band , and ( ii ) the low energy index can , in cases to be discerned in sections  [ sec : icformalism ] and  [ sec : alphalow_comp ] , be as flat as . at the same time , inverse compton scenarios ( particularly the ssc case ) have the disadvantage of requiring that synchrotron radiation or other components be unobservable at optical and lower wavebands during the prompt and early afterglow phases .
the parameter constraints thereby imposed are discussed in section  [ sec : constraint ] .    the well - known criterion for ic to dominate synchrotron as an electron cooling mechanism ( i.e. also in total bolometric luminosity for a spatially uniform medium )
is simply that the energy density of the soft photons exceeds that of the ambient magnetic field , ( e.g. rybicki & lightman , 1979 ) . in general
, the value of can vary enormously depending on the particular invocation for soft photon generation . for the specific case of a ssc model , couples to the field and can be determined by integrating in eq .
( [ eq : synch_emiss ] ) over and multiplying by some fiducial photon escape timescale , , that reflects the physical size of the emission region . for monoenergetic electrons of lorentz factor , the familiar result ( e.g.
see p.  168 of rybicki & lightman , 1979 ) for the synchrotron energy loss rate can be used to write .
hence , the ssc emission can dominate the bolometric luminosity for significant but not particularly large values of the product .
it was precisely this property that led gonzalez et al .
( 2003 ) to argue that the intriguing hard gamma - ray component seen by the egret tasc in grb 941017 was not due to ssc emission , but perhaps due to inverse compton scattering seeded by an external soft photon source .
the criterion for an ssc component to fall in the soft gamma - ray band is , which can easily accommodate fields in the gauss range for and .
a more thorough exploration of ssc parameter spaces for burst prompt and afterglow emission can be found , for example , in mszros , rees , & papathanassiou ( 1994 ) , sari & esin ( 2001 ) and zhang & mszros ( 2001 ) .
the formalism for the inverse compton ( ic ) emissivity in optically thin environs is readily available in many astrophysical textbooks , such as rybicki & lightman ( 1979 ) .
since the description of ic polarization effects is more subtle than for synchrotron radiation , being dependent on anisotropies in the soft photon population , the developments here are restricted to cases where photon polarizations are summed over .
the rate of emission , integrated over the angles of the upscattered photons , is independent of the direction of momentum of the seed photons .
hence , no integration over is required , contrasting eq .
( [ eq : synch_emiss ] ) for synchrotron radiation .
however , the rate is dependent on the soft photon energy , which must be integrated over , yielding the form ( in units of ) @xmath17 here is the thomson cross section , is the energy distribution of soft photons , and @xmath18\ ;    \theta ( z ) \quad , \quad    z\ ; = \ ; \dover{\erg}{4\gamma^2\erg_{\rm s}}\;\ ; ,   \label{eq : fdef}\ ] ] for a heaviside step function when and otherwise .    mirroring the treatment of synchrotron radiation , while numerical evaluations of eq .
( [ eq : ic_emiss ] ) were performed using the electron distribution in eq .
( [ eq : elec_dist ] ) to obtain approximations to burst spectra in section  [ sec : part_comp ] , asymptotic spectral indices are readily attainable . for monoenergetic soft photons , the non - thermal contribution to the inverse compton spectrum ( from the truncated power - law ) is dominant whenever , and assumes the following power - law forms : @xmath19 the index at high energies , , is identical to that for synchrotron radiation , since both mechanisms emit photons with the same dependence on the electron lorentz factor : . at low energies ,
the differential spectrum is approximately independent of the upscattered energy of the photon , due to the limit of and the fact that electrons near the cutoff dominate the low energy photon production .
this index of has been known since the early days of application of inverse compton to astrophysical problems : see , for example , figure  3 of the comprehensive treatment of jones ( 1968 ) that also illustrates that the compton spectrum flattens further to in the domain of downscattering .    as in the case of synchrotron radiation , the thermal component , sampled when , necessarily declines exponentially at photon energies well above the thermal peak , but possesses the same power - law dependence as in the non - thermal case for low photon energies : @xmath20 this follows because the low energy photons are produced by electrons mostly near the peak of the thermal distribution , again sampling const .      from the perspective of spectral shapes at and below the gamma - ray peak , it is important to explore the properties of synchrotron self - compton ( ssc ) radiation , to distinguish its character from ic emission that is seeded by quasi - monoenergetic soft photons . here
the analysis will first be confined to the case of a purely non - thermal electron distribution , i.e. regimes , where the distribution can be written in the form for and is zero otherwise .
hereafter , when making the identification with the non - thermal portion of the distribution function in eq .
( [ eq : elec_dist ] ) , , and one can set and to establish a correspondence .
for such a truncated power - law , a compact analytic form for the ssc spectrum was obtained , essentially reducing the quadruple integration to one over a single variable . defining a fiducial photon escape timescale , for a one - zone emission region of size , the average density of synchrotron photons that seed the ssc emission can be written as , so that the form in eq .
( [ eq : synchfinal ] ) can be directly inserted into eq .
( [ eq : ic_emiss ] ) .
the analytic reduction is expounded in appendix a , with the result being @xmath21   \label{eq : sscfinal}\\[-5.5pt ]     & \times & \biggl\ {       \int_{0}^{\chi } dy \ , k_{5/3}(y )          p_{\delta } \bigl(\dover{y}{\chi}\bigr ) +       \int_{\chi}^{\infty } dy \ ,
k_{5/3}(y )          q_{\delta } \bigl(\dover{\chi}{y}\bigr ) \ ; \biggr\}\;\ ; , \nonumber\end{aligned}\ ] ] where @xmath22 is the effective scale for the rate of the synchrotron - self - compton process . here , is the thomson optical depth of the emission region .
the emergent photon energy is encapsulated in a single dimensionless parameter @xmath23 that denotes the energy scale of the peak for ssc emission , which putatively could be in the gamma - ray band , and is the only spectral structure in the emissivity .
the functions and are relatively compact and are given in eqs .
( [ eq : pdef ] ) and  ( [ eq : qdef ] ) in appendix a in terms of elementary functions .
( [ eq : sscfinal ] ) can be routinely integrated numerically .
asymptotic energy dependences of the ssc rate in eq .
( [ eq : sscfinal ] ) are readily obtained . in cases
where , the term never contributes to the integration since it is exponentially suppressed by the behavior of the bessel function .
accordingly , it is quickly deduced that the ssc spectrum approximately obeys at energies , as expected : a convolution of an extended synchrotron spectrum with a truncated power - law electron distribution with index is a power - law of index .
there is a slight logarithmic modulation as outlined in eq .
( [ eq : sscnt_highchi ] ) .
the limit can also be promptly obtained .
the term samples only the asymptotic limit of the bessel function , i.e. the form . from this
, it is quickly verified that the integration is proportional to .
it is simple to show using the eqs .
( [ eq : pdef ] ) and  ( [ eq : qdef ] ) in appendix a that the term possesses an identical dependence .
the net result is , as embodied in eq .
( [ eq : sscnt_lowchi ] ) , a property that has been evident in various ssc models of bursts ( e.g. see sari & esin , 2001 ) .
these two limits can be summarized via @xmath24 the low energy index images the synchrotron low energy tail , a result that is not surprising since such a tail is steeper ( and therefore its constituent photons are more populous ) than the tail resulting from ic upscattering of monoenergetic soft photons in eq .
( [ eq : ntic_limits ] ) .
accordingly , the bound is placed on ssc spectra so that this mechanism can not aid with the low energy spectral index issues discussed in section  [ sec : alphalow_synch ] .        to facilitate the fitting of grb spectra with ssc spectra , the case of ultra - relativistic thermal particles
must also be considered , for which the distribution can be written in the form . when making the identification with the thermal part of the distribution function in eq .
( [ eq : elec_dist ] ) , one can set and to establish a correspondence with correct normalization . in appendix
b , an expedient derivation of the emission spectrum is offered , making use of the result already presented in eq .
( [ eq : sscfinal ] ) for a truncated power - law electron distribution .
basically , the limit is taken to generate the spectrum for monoenergetic electrons , and this is then convolved with a relativistic maxwell - boltzmann distribution . setting @xmath25 as the effective scale for the rate of the synchrotron - self - compton process for the hot thermal particles , one arrives at @xmath26^{1/4}\bigr)\quad ,   \label{eq : sscfinal_th_too}\ ] ] where the function is defined in eq .
( [ eq : caljdef ] ) , and is steeply declining as increases . here , .
asymptotic forms for the spectrum in eq .
( [ eq : sscfinal_th_too ] ) are routinely derived , being presented in eq .
( [ eq : sscthlimits ] ) .
as expected , an form arises at energies below the characteristic energy , and the spectrum is a weakly declining exponential far above this .
results for the spectral fitting of the grb 910503 data for the case of inverse compton scattering are presented in figure  [ fig : ic_spec_dist ] .
the fits were performed using the formalism presented in section  [ sec : icformalism ] , specifically with numerical evaluations of eq .
( [ eq : ic_emiss ] ) .
again , the multi - instrument , time - integrated cgro data as published in schaefer et al .
( 1998 ) was used .
the figure displays both inverse compton emission spectra for the case of an isotropic , monoenergetic soft photon distribution , and the inferred electron distribution that generates the continuum .
the solid curve depicts the case that would be required to explain the photon spectrum , with the overall normalization again being a free parameter ; this curve was offset slightly from the data for visual clarity .
the electron distributions were normalized to unit area in space .
the non - thermal particle distribution clearly dominates the thermal component in this `` fit , '' as is the case for synchrotron emission .
the index was controlled by the egret data , and so is identical to that in fig .
[ fig : synch_spec_dist ] .
adjustment of the parameter to improve the appearance of the model approximation around and above the spectral peak led to the assumed value of , again accurate to around 15% .
the positioning of the peak in energy produced the fit value , where is the lorentz factor associated with the bulk motion of the grb expansion .
the parameter describing the relative normalization of the non - thermal and thermal components was found to be ( % ) , still a large value though only a tenth of that obtained for the synchrotron fit .
the conclusion is essentially identical to that in section  [ sec : part_synch ] : the observational data can only be accommodated by an electron distribution that has a preponderance of non - thermal electrons as fig .
( [ fig : ic_spec_dist]b ) .
such a determination was also obtained for batse bursts grb 940619 ( trigger number 3035_6 ) and grb 940817 ( trigger number 3128_5 ) from the preece et al .
( 2000 ) catalog .
moreover , the fits in figs .
[ fig : synch_spec_dist ] and  [ fig : ic_spec_dist ] are of comparable quality , at least below 1 mev and above 5 mev ; i.e. there is little potential in the cgro data to discriminate between the two emission mechanisms .
there is a small bump / excess in the theoretical spectrum in fig .
[ fig : ic_spec_dist ] just above 1 mev that is a consequence of the sharper inverse compton scattering kernel ; its prominence can be muted somewhat by slightly reducing both and the index .
such a liberty may not be afforded in the glast era when burst spectral indices will be better constrained because of its broad energy range , sampled by the gamma - ray burst monitor ( gbm ) in the batse band and the large area telescope ( lat ) in the egret band .
the last point to note is that the inverse compton fit is tolerated here because the low energy asymptotic index is possibly not fully realized in the batse window .
the data possess a slight curvature that may be consistent with a continued flattening of the spectrum below 30 kev , as would be expected for the inverse compton mechanism . yet
the data match the synchrotron emissivity also , indicating that batse can not discriminate between these two mechanisms on the basis of the soft gamma - ray portion of the bandwidth .
this impasse should be addressed when swift is launched in the coming year .
swift possesses the sensitivity in the critical 1 -100 kev band to probe the spectrum below the peak and afford diagnostics on the emission mechanism .
furthermore , the spectral coverage afforded by the gbm and the lat on glast will also provide powerful observational constraints .
the low energy asymptotic forms in eqs .
( [ eq : ntic_limits ] ) and  ( [ eq : thic_limits ] ) raise the possibility that inverse compton emission can eliminate the so - called `` line of death '' problem that is pervasive in synchrotron models .
the ic mechanism can accommodate low energy asymptotic indices as flat as , with steeper indices easily being generated by appropriately distributed electrons .
the domain of compatibility is indicated in fig .
[ fig : batse_low_index ] , where around 98% of batse burst spectra in the spectroscopy catalogue of preece et al .
( 2000 ) can be accounted for by the inverse compton process .
this is an enticing property , but is not a definitive statement of the operation of this mechanism in grb prompt emission regions .
it must be emphasized that the establishment of an domain is contingent upon the soft photons that seed the upscattering assuming a narrow energy ( i.e. quasi - monoenergetic ) distribution .
this caveat is exemplified by the ssc value of being identical to that for synchrotron radiation due to the breadth of the synchrotron seed spectrum . hence invoking a synchrotron self - compton model
only dramatically impacts the operating range of and not the spectral shape in the very soft gamma - rays .
given the discussion in section  [ sec : anis_synch ] pertaining to synchrotron radiation , it is natural to ask whether anisotropies in either the soft photon or electron populations can dramatically alter the spectrum .
the inverse compton photon production rate in eq .
( [ eq : ic_emiss ] ) implicitly assumes isotropy of the soft photons constituting .
furthermore , it assumes that the angular distribution of electrons is not seriously deficient along the line of sight to the observer ; i.e. that the viewing perspective is not towards the side of an electron jet .
calculating ic spectra for anisotropic distributions is quite involved , and moreover sensitive to the choice of geometry . yet
spectral behavior depends mostly on the interplay between such geometry and the kinematics of ic scattering .
if the electron distribution is highly collimated and the observer s direction lies outside the cone of jet collimation , then the kinematics of scattering requires a smaller angle of scattering in the electron rest frame to send the photon towards the observer .
this small yields a lower final photon energy in the observer s frame , so that the emitted inverse compton spectrum is softer outside the jet core .
the spectrum then flattens at lower energies ( e.g. see dermer & schlickeiser 1993 ; baring 1994 for blazar contexts ) due to the kinematic depletion of photons in the observer s direction , and can be even slightly flatter than .
a problem with this scenario immediately arises : the energetics are dominated by emission in the jet cone so that assuming `` off - axis '' ic emission places more severe constraints on the grb radiative efficiency .
note that fluctuations in the jet orientation or structure would lead to the emergent inverse compton spectra being highly variable both in flux ( observed in bursts and therefore not a problem , in principle ) and in spectral shape  any structure such as the peak would be expected to vary in time rapidly in energy due to the scattering kinematics .
while modest fluctuations ( as opposed to hard - to - soft evolution ) are present in gamma - ray burst prompt emission , detailed modeling would be required to discern whether they could match those produced by anisotropic scattering scenarios .
any anisotropy in the soft photon population is not sampled until they become more collimated than the electrons ; only then does the emissivity change appreciably from eq .
( [ eq : ic_emiss ] ) .
then , electrons within the photon `` cone '' of beaming experience an almost full range of incident photon angles in their rest frames , so that they contribute an inverse compton spectrum similar to that for isotropic soft photons .
however , the majority of electrons lie outside this cone , and for these , collisions with photons are predominantly head - on , leading to a depletion of the softer photons from the overall inverse compton spectrum just as with electron collimation . hence
, one expects that such soft photon anisotropies could in principal also flatten the spectrum below the peak .      to provide a meaningful comparison with the preceding mechanisms , grb 910503
is again used as a focal point . yet , so as to provide a variation , here a case study for another egret burst , grb 910601 is also presented .
results for the spectral fitting for synchrotron self - compton scattering are presented in figure  [ fig : ssc_spec_dist ] .
the `` fits '' were performed using the formalism presented in section  [ sec : sscformalism ] , specifically with numerical evaluations of eqs .
( [ eq : sscfinal ] ) and  ( [ eq : sscfinal_th_too ] ) .
the figure displays the ssc emission spectra for the case of an isotropic electron distribution , for each burst .
the solid curve depicts the case that would most closely explain the photon spectrum , with the overall normalization again being a free parameter . due to the strong dependence of the characteristic energy on electron lorentz factor ,
the values of required to move the ssc peak into the gamma - ray band lead to inferences of sub - gauss fields if .
consider first grb 910503 .
the most obvious features are that the `` fit '' is more difficult than for the emission mechanisms previously considered , with poor results for the lower energies much below the peak , and that the thermal contribution must be dramatically suppressed .
the cause is easily identified : the ssc spectrum is intrinsically much broader than either the ic or synchrotron spectra .
this is due to the characteristic energy for ssc embodied in the dimensionless parameter ( see equation  [ [ eq : chidef ] ] ) scaling as ( or for the thermal case ) .
consequently , small changes in the electron near the peak of the maxwell - boltzmann distribution incur large variations in the ssc photon energy so that the spectral peak is broad .
this effect is manifested in the non - thermal ssc emissivity in eq .
( [ eq : sscfinal ] ) , which leads to the excesses evident around 30 kev in the `` fits '' of figure  [ fig : ssc_spec_dist ] , and more profoundly for the thermal ssc contribution , which is extremely broad . from eq .  ( [ eq : sscthlimits ] ) a notable conclusion can be drawn : that the thermal ssc spectrum in the depicted representation only starts to fall above around 30 mev , and then only slowly , i.e. it is blueward of the peak from the non - thermal ssc contribution .
this drives the requisite dramatic suppression of the normalization of a thermal component to the underlying electron distribution , leading to the same issues ( but in this case more strikingly emphasized ) for acceleration from thermal energies that were discussed above .
note that for values of , a comparably satisfactory ssc fit can be obtained with .
observe that for such a large , the tail of the thermal maxwellian emerges above that of the non - thermal distribution for a small range of lorentz factors .
the situation for grb 910601 is even more constrained , due to the markedly steeper gamma - ray spectrum .
the thermal electrons must be even less populous to accommodate the comptel and egret data , and the non - thermal electrons must possess a very steep index , considerably steeper than the required for a synchrotron fit .
in fact , both the synchrotron and ic models had much greater success fitting the grb 910601 data , though are not depicted .
steepening the index to produces only marginal improvement in the fit .
the inevitable conclusion of this analysis is that the breadth of the ssc emissivity for even monoenergetic electrons ( derived in equation  ( [ eq : sscfinal_mono ] ) , and closely approximated by the case in figure  [ fig : ssc_spec_dist ] ) is large enough to provide significant or severe problems for the ssc interpretation of prompt emission from gamma - ray bursts .
this contention is not based solely on the grb 910503 and in particular the grb 910601 data ; the similarity of the spectra of grb 910814 ( schaefer et al .
1998 ) and grb 990123 ( briggs et al .
1999 ) to that of grb 910601 strongly suggests that a spectroscopic incompatibility of the ssc mechanism for prompt emission may be commonplace , specifically for differential burst spectra steeper than around above the peak .
bursts such as the superbowl one ( grb 930131 ) with high energy index may be much more compatible with an ssc scenario .
while this assessment may achieve closure in the swift / glast era , a caveat is possibly provided if the synchrotron continuum is strongly self - absorbed .
this is the scenario envisaged by panaitescu & meszaros ( 2000 ) and liang , boettcher & kocevski ( 2003 ) , with the latter group suggesting that ssc with moderate optical depths (  ) could explain x - ray excesses seen in some burst spectra . with such strong absorption ,
the supply of soft photons for inverse compton upscattering is more sharply peaked than even the synchrotron continuum from electrons possessing a truncated power - law distribution , thereby approximately mimicking the pure inverse compton spectral shapes explored in section  [ sec : part_comp ] .
the requirement that self - absorption be effective in narrowing the synchrotron self - compton peak is that the turnover energy exceeds that of the synchrotron peak , which is below the x - ray band . using eq .
( [ eq : selfabs_energy ] ) , one quickly arrives at the requisite condition : @xmath27 for @xmath28 here is the energy of the peak in the prompt emission . for , one deduces that  , so that optical depths can be tolerated for typical without pushing too low , i.e. into the mildly - relativistic regime .
detailed computations of self - absorbed ssc spectra are deferred to future work .
in addition to the soft gamma - ray spectral shape , a probing diagnostic on the emission mechanisms is whether their extrapolation down to the optical band can be accommodated by optical monitoring data for bursts .
such broadband constraints turn out not to be profoundly limiting for synchrotron or pure inverse compton emission .
however , they do impact ssc models , since these can potentially have a prominent synchrotron component in the optical that might violate upper bounds or , in the case of grb 990123 , prompt optical detections .
it is straightforward to assess how optical band data can constrain ssc model parameters using the formalism derived here .
first , the focus is on observational data .
grb 990123 ( with redshift ) is to date the only burst with a prompt optical detection ( akerlof et al . 1999 ) , peaking at around magnitude 8.8 some 47 seconds after batse trigger and contemporaneous with the batse signal .
these prompt observations were made by the robotic optical transient search experiment - i ( rotse - i ) .
the recent bright burst grb 030329 ( redshift ) was only detected optically by rotse - iii in the afterglow epoch at visual magnitude 12.5 over an hour after the burst trigger on hete-2 ( smith et al .
kehoe et al .
( 2001 ) report upper bounds of around magnitude 1315 to half a dozen bursts , in two cases starting just 15 - 20 seconds after trigger .
these low limiting magnitudes can potentially be severely constraining on models with a high optical to gamma - ray flux expectation , though contemporaneous multiwavelength detections are requisite for unambiguous diagnostics .
taking magnitude 9 as a benchmark optical signal for the purposes of exploring broadband constraints , a number that must be considered modulo the source distance and intrinsic luminosity , this corresponds to a differential photon flux of the order of photons sec@xmath29
@xmath30 kev@xmath29 at around 2ev .
this can be compared with the batse spectrum of grb 910503 , a burst brighter than most of those in the sample of kehoe et al .
( 2001 ) , which has a time - integrated differential flux of photons sec@xmath29 @xmath30 kev@xmath29 at around 100 kev .
the interpolation between such an optical magnitude and gamma - ray flux is a roughly spectrum .
immediately , one infers that synchrotron and inverse compton models for the gamma - rays can accommodate such an optical magnitude ; only limiting magnitudes of around 19 or higher can provide conflicts for these mechanisms .
note that a similar conclusion would apply for pitch angle synchrotron models , which possess very flat low energy asymptotic spectra like inverse compton scenarios .
the situation for ssc broadband spectra is very different .
the ratio of the energy of the ssc peak to that of the synchrotron peak is simply of the order , which would set the synchrotron peak in the infra - red for lorentz factors of , independent of the bulk flow .
lower can move the synchrotron peak into the optical and even the uv .
the ratio of emissivities at the synchrotron and ssc peaks is quickly determined from equations  ( [ eq : ndotsyn_def ] ) and  ( [ eq : ndotsscdef ] ) to be in the notations of this paper , i.e. with .
hence , for , the peaks of the two components would lie roughly on an spectrum .
moreover , opting for to position a synchrotron peak in the optical window , normalizing to the batse spectrum of grb 910503 would yield a count rate an order of magnitude below the magnitude 9 benchmark posited above . at face value
, this would render the ssc scenario compatible with limiting magnitudes of around 11 for this source .
yet comptonization of the spectrum is not evident in the gamma - rays , so it is inferred that , thereby increasing the ratio and lowering the permissible limiting magnitude for a viable ssc model .
accordingly , for the only source with a prompt optical detection , grb 990123 , for the ssc model not to overestimate the observed optical flux , the synchrotron peak has to move somewhat out of the optical , though not much .
this amounts to a modest constraint in phase space that can be quickly deduced , though at this juncture is not particularly enlightening ; such an exercise will be more timely once swift is providing a number of simultaneous hard x - ray and optical detections of bursts .
the interpretative modeling presented here provides a strong indication that the synchrotron and inverse compton processes may provide a more viable explanation of prompt burst spectra than the ssc mechanism .
the spectroscopic bias against ssc is due to its inherently broad character around the peak , and is only marginally constrained by broad - band considerations , i.e. limiting magnitudes in the optical band .
however , strong self - absorption of the synchrotron component can provide an escape clause for ssc scenarios , though this imposes significant fine - tuning of parameter space .
inverse compton provides an attractive possibility for the vast majority of bursts , as does pitch angle synchrotron emission , though quasi - isotropic synchrotron radiation may only be reconciled with around 2/3 of bursts at the lowest energies in the batse window .
discrimination between these processes may prove possible after swift launches and in the glast era , though polarimetric measurements would enhance the diagnostic capability of observations . in analyzing ssc spectral models ,
attractively compact analytic forms for the synchrotron self - compton emissivity from thermal and truncated power - law electron distributions were developed ; these will prove useful for applications of ssc spectral models to a variety of astrophysical environments .
a distinctive conclusion of this work is that agreeable `` fits '' with either synchrotron or inverse compton emission scenarios are only attainable if the underlying electron distribution is profoundly dominated by a non - thermal component .
this result is in contradistinction with standard understanding of acceleration mechanisms such as the diffusive ( fermi ) process at shocks , where the non - thermal particles are drawn probabilistically from the shocked thermal population , and therefore are subdominant in number .
this constraint has profound implications for any proposed grb model , providing a challenge for theorists in an era when computational resources are permitting new avenues for exploration using intensive numerical simulations of particle acceleration at shocks .
this work was supported by the nsf extragalactic astrophysics program .
matthew braby acknowledges the support of a shell foundation scholarship towards this research during his undergraduate degree at rice university .
we thank the referee , rob preece , for insightful and constructive suggestions for polishing the paper , and nicole lloyd - ronning , jerry fishman and demos kazanas for reading through the manuscript and providing helpful comments .
we also thank brad schaefer for supplying the grb 910503 compilation data prior to its publication in 1998 .
in this appendix , developments leading to the form in eq .
( [ eq : sscfinal ] ) for the rate of synchrotron self - compton emission from a truncated power - law distribution of electrons are expounded .
the starting point is the insertion of times the synchrotron rate in eq .
( [ eq : synchfinal ] ) , serving as the soft photon density , directly into the expression in eq .
( [ eq : ic_emiss ] ) for the inverse compton emissivity .
this presents a triple integral that can be reduced to a compact single integration by appropriate manipulations .
the lorentz factor of the electron that upscatters the synchrotron photons proves to be a less convenient integration variable than , which is defined in eq .
( [ eq : fdef ] ) . likewise , to proves expedient to use @xmath31 as an integration variable , rather than the synchrotron photon energy . performing these two changes of variables leads quickly to an ssc emissivity of @xmath32 where is the ssc rate scale factor given in eq .
( [ eq : ndotsscdef ] ) , and is the dimensionless ssc energy variable defined in eq .
( [ eq : chidef ] ) . here , and
the maximum value of can be written @xmath33 the -integration over the inverse compton scattering kernel leads to a definition of : @xmath34   \label{eq : gzdef}\\[-5.5pt ]     g_{\delta}(z ) & = & \dover{2z}{\delta+3}\ ;     \biggl ( 2\log_ez + \dover{\delta-1}{\delta+3 } \biggr )     + \dover{2}{\delta+1 } - \dover{4z^2}{\delta+5}\quad , \nonumber\end{aligned}\ ] ] the integration can be reduced by noting that the definitions of and establish two critical values of , namely and . which delineate distinct dependences of the integrand on .
the relative ordering of and depends on the value of assumed : and .
this naturally divides the integration into two ranges , , and .
the developments are straightforward , and are facilitated using an function defined via @xmath35   \label{eq : hzdef}\\[-5.5pt ]     h_{\delta}(z ) & = & \dover{4z}{(\delta+3)^2}\ ;     \biggl ( 2\log_ez + \dover{\delta-5}{\delta+3 } \biggr )     + \dover{4}{(\delta+1)^2 } - \dover{8z^2}{(\delta+5)^2}\quad .\nonumber\end{aligned}\ ] ] the case requires an integration by parts in one contributing interval , and one quickly arrives at the compact form @xmath36 i.e. eq .  ( [ eq : sscfinal ] ) , using the following functional forms ( with and defined as just above ) : @xmath37\ , g_{\delta}(1 ) + h_{\delta}(1)\,\biggr\ } \;\ ; ,   \label{eq : pdef}\ ] ] for , and for @xmath38 + h_{\delta}(z)\,\biggr\ } \quad , \nonumber\\[-5.5pt ]   \label{eq : qdef}\\[-5.5pt ]         & = & - 2z\,\dover { \delta + 1}{(\delta + 3)^2}\biggl [ 2 \ ,             \log_ez - \dover{\delta + 11}{\delta + 3 } \biggr ]        - 4 \dover{\delta -1}{(\delta + 1)^2 } -2 \dover{\log_ez}{\delta + 1 }         + 2z^2\ , \dover{\delta + 1}{(\delta + 5)^2}\quad .\nonumber\end{aligned}\ ] ] observe that as , , as demanded by continuity of the spectrum .
asymptotic dependences of the ssc rate on the upscattered photon energy , encapsulated in the proportionalities of eq .
( [ eq : ntssc_limits ] ) , can be simply derived . in the limit , both
the and terms contribute to eq .
( [ eq : sscfinal_app ] ) , with the integrals being dominated by the domain .
the bessel function can be replaced by its power - law asymptotic form , resulting in elementary integrals . for the integration
converges to @xmath39 the case is less simple , though the integration over does not contribute significantly since it is exponentially suppressed by the dependence of the bessel function .
the integration can be developed using identity 6.561.16 of gradshteyn and ryzhik ( 1980 ) , which can be specialized to the case appropriate for this analysis : @xmath40 the derivative of this with respect to can also be used to expedite the developments , and can be written in the form @xmath41 where is the digamma function . setting , since the upper limit to
the can be replaced by infinity , one readily arrives at the result @xmath42   \label{eq : sscnt_highchi}\\[-5.5pt ]      & & \hbox{\hskip 0pt }           - \dover{1}{2 }   \psi \biggl ( \dover{\delta}{4 } + \dover{19}{12}\biggr )          - \dover{1}{2 } \psi \biggl ( \dover{\delta}{4 } - \dover{1}{12}\biggr)\ , \biggr ]     \ , g_{\delta}(1)+ h_{\delta}(1 ) \biggr\}\ ;     \dover{\ndotssc}{\chi^{(\delta + 1)/2}}\quad , \quad \chi\,\gg\ , 1\;\ ;   .\nonumber\end{aligned}\ ] ] the logarithmic term will clearly dominate for extremely large values of .
the comparison of synchrotron self - compton ( ssc ) models invoking shock acceleration with burst spectra requires also the ssc rate for relativistically hot thermal particles , which implies using a four - dimensional integral defined by inserting the thermal contribution to eq .
( [ eq : elec_dist ] ) into eq .
( [ eq : synch_emiss ] ) and then the result subsequently into eq .
( [ eq : ic_emiss ] ) . rather than laboriously reducing the integral along the lines of the developments of appendix a , the result in eq .
( [ eq : sscfinal_app ] ) can be used to advantage . taking the limit generates the ssc emissivity for monoenergetic electrons with lorentz factor .
the result can then be convolved with the thermal electron distribution to yield the required emissivity .
as becomes very large , approaches zero at all positive less than unity , so that the integration over in eq .
( [ eq : sscfinal_app ] ) does not contribute . using the result @xmath43
the integration over the thermal electron distribution produces a double integral that can be evaluated by reversing the order of integrations , yielding @xmath26^{1/4}\bigr)\quad ,   \label{eq : sscfinal_th}\ ] ] where is given in eq .
( [ eq : ndotsscdef ] ) , , and @xmath44   \label{eq : caljdef}\\[-5.5pt ]     & \equiv & 2\dover{e^{-\kappa}}{\kappa } \int_0^{\infty }      \dover{e^{-v}s(v)}{\kappa + v}\ , dv       \ ; = \ ; 2\dover{e^{-\kappa}}{\kappa^3 } \int_0^{\infty }      \dover{e^{-v } \ , v^2\ , s(v)}{\kappa + v}\ , dv \quad , \nonumber\end{aligned}\ ] ] where @xmath45 the first integral defining can be expressed as a sum of incomplete gamma functions and one of their derivatives , using identities 3.381.3 and 4.358.1 of gradshteyn and ryzhik ( 1980 ) , though this does not expedite the evaluation of . however , the integral representation of the incomplete gamma function in 8.353.3 of gradshteyn and ryzhik ( 1980 ) proves useful , leading to the alternative integral forms posited on the second line of eq .
( [ eq : caljdef ] ) .
these latter forms are readily amenable to numerical computation in a fashion similar to evaluation of euler s integral form for the gamma function .
note that the first form on the second line is expedient for cases , while the second form on this line , which was obtained using the fact that the integrals of and on are identically zero , is advantageous for evaluation when .
asymptotic forms for are simply obtained : for , and for .
the appropriate asymptotic limits of eq .
( [ eq : sscfinal_th ] ) can then be quickly determined : @xmath46 ^ 2\;\chit^{-2/3}\;\ ; , &         $ \quad \chit\ll 1\quad , \vphantom{\biggl(}$\cr            \dover{8\pi}{\sqrt{5}}\ ; \dover{2^{9/10}}{\chit^{4/5}}\ ;            \exp \biggl\ { -\dover{5}{2^{8/5}}\ ; \chit^{1/5 } \biggr\}\;\ ; , &          $ \quad \chit\gg 1\quad , \vphantom{\biggl(}$\cr }   \label{eq : sscthlimits}\ ] ] where the low argument power - law dependence for the bessel function is used when , and the method of steepest descents is employed to derive the case , for which the exponential portion of is sampled .    note that the ssc emissivity for a monoenergetic electron distribution can be determined by omitting the convolution with the thermal distribution , so that eq .
( [ eq : qlimit ] ) can be used efficaciously to yield @xmath47 as is a monotonically decreasing function of on , for it flattens the ssc spectrum relative to that for synchrotron emission , and so effects a broadening of the peak .
this monoenergetic formula can be convolved with arbitrary electron distributions to obtain resultant synchrotron emissivities . by the same token , eq .
( [ eq : sscfinal_th ] ) can be used to facilitate a laplace transform approach for special cases for the electron distribution .
akerlof , c. , et al .
1999 , @xmath48398 , 400 .
band , d. , et al .
1993 , @xmath48413 , 281 .
baring , m.  g. 1988a , @xmath48235 , 51 . baring , m.  g. 1988b , @xmath48235 , 79 . baring , m.  g. 1990 , @xmath48244 , 49 . baring , m.  g. 1994 , in _ proc . of the second compton symposium _ , eds .
gehrels , n. & norris , j. , ( aip conf .
304 , new york ) p.  639 . baring , m.  g. 1995 , proc .
of the erice nato asi _ currents in high energy astrophysics , eds .
shapiro , m. , wefel , j. and silberberg , r. , ( kluwer academic , dordrecht ) p.  21 .
baring , m. g. 1999 , in _ proc .
of the 26th international cosmic ray conference _ , vol .
iv , p.  5
, held in salt lake city , utah . baring , m.  g. 2002 , pub .
@xmath4819 , 60 .
baring , m.  g. , ellison , d.  c. & jones , f.  c. 1993 , @xmath48409 , 327 . baring , m.  g. , ogilvie , k.  w. , ellison , d.  c. , & forsyth , r.  j. 1997 , @xmath48476 , 889 .
bekefi , g. 1966 , _
radiation processes in plasmas , _ ( wiley and sons , new york ) blandford , r.  d. & eichler , d. 1987 , phys .
@xmath48154 , 1 .
brainerd , j.  j. & lamb , d.  q. 1987 , @xmath48313 , 231 .
brainerd , j.  j. & petrosian , v. 1987 , @xmath48320 , 703 .
briggs , m.  s. , et al .
1999 , @xmath48524 , 82 .
crider , a. & liang , e.  p. 1999
, @xmath48138 , 405 . de jager , o.  c. , et al .
1996 , @xmath48457 , 253 .
dermer , c.  d. , chiang , j. & bttcher , m. 1999 , @xmath48513 , 656 .
dermer , c.  d. & schlickeiser , r. 1993 @xmath48416 , 209 .
drury , l.  oc .
@xmath4846 , 973 .
ellison , d.  c. , baring , m.  g. & jones , f.  c. 1996 , @xmath48473 , 1029 .
ellison , d.  c. & double , g.  p. 2002
, @xmath4818 , 213 .
ellison , d.  c. , jones , f.  c. , & reynolds , s.  p. 1990
, @xmath48360 , 702 .
ellison , d.  c. , mbius , e. , & paschmann , g. 1990 , @xmath48352 , 376 .
epstein , r.  i. 1973 , @xmath48183 , 593 .
epstein , r.  i. & petrosian , v. 1973 , @xmath48183 , 611 .
forman , m.  a. , jokipii , j.  r. & owens , a.  j. 1974 , @xmath48192 , 535 .
gallant , y. a. , hoshino , m. , langdon , a. b. , arons , j. & max , c. e. 1992 , @xmath48391 , 73 .
ghisellini , g. , celotti , a. & lazzati , d. 2000 , @xmath48313 , l1 .
giacalone , j. , burgess , d. , schwartz , s.  j. & ellison , d.  c. 1992 , @xmath4819 , 433 .
giacalone , j. , burgess , d. , schwartz , s.  j. & ellison , d.  c. 1993 , @xmath48402 , 550 .
goodman , j. 1986 , @xmath48308 , l47 .
gonzalez , m.  m. , et al .
2003 , @xmath48424 , 749 .
gosling , j.  t. , et al .
1981 , @xmath4886 , 547 .
gradshteyn , i.  s. and ryzhik , i.  m. 1980 , _ table of integrals , series and products _ , ( academic press , new york ) granot , j. , piran , t. & sari , r. 2000 , @xmath48534 , l163 .
hoshino , m. , arons , j. , gallant , y.  a. & langdon , a.  b. 1992 , @xmath48390 , 454 .
hoshino , m. & shimada , n. 2002 , @xmath48572 , 880 .
jackson , j.  d. 1975 , _ classical electrodynamics , _ ( wiley and sons , new york ) jokipii , j.  r. , kta , j. , & giacalone , j. 1993 , @xmath4820 , 1759 .
jones , f.  c. 1968 , @xmath48167 , 1159 .
jones , f.  c. & ellison , d.  c. 1991 , @xmath4858 , 259 .
jones , f.  c. , jokipii , j.  r. & baring , m.  g. 1998 , @xmath48509 , 238 .
kang , h. & jones , t.  w. 1997 , @xmath48476 , 875 .
kehoe , r. , et al .
2001 , @xmath48554 , l159 .
kucharek , h. , & scholer , m. 1995 , @xmath48100 , 1745 .
liang , e.  p. , boettcher , m. & kocevski , d. 2003 , in _ gamma - ray burst and afterglow astronomy 2001 _ , eds .
g.  r. ricker & r.  k. vanderspek ( aip conf .
662 ) , p.  295 .
liewer , p.  c. , goldstein , b.  e. & omidi , n. 1993 , @xmath4898 , 15,211 .
lloyd , n.  m. & petrosian , v. 2000 , @xmath48543 , 722 .
medvedev , m.  v. 2000 , @xmath48540 , 704 .
mszros , p. 2001 , science @xmath48291 , 79 .
mszros , p. , rees , m.  j. & papathanassiou , h. 1994 , @xmath48432 , 181 .
pacyski , b. 1986 , @xmath48308 , l43 .
panaitescu , a. & mszros , p. 2000 , @xmath48544 , l17 .
pavlov , g.  g. & golenetskii , s.  v. 1986 , @xmath48128 , 341 .
piran , t. 1999 , phys .
@xmath48314 , 575 .
preece , r.  d. , et al .
1998 , @xmath48506 , l23 .
preece , r.  d. , et al .
2000 , @xmath48126 , 19 .
preece , r.  d. , et al .
2002 , @xmath48581 , 1248 .
rees , m. j. , & mszros 1992 , @xmath48258 , 41 .
rybicki , g.  b. & lightman , a.  p. 1979 , _ radiative processes in astrophysics , _
( wiley and sons , new york ) sari , r. & esin , a. 2001 , @xmath48548 , 787 .
schaefer , b.  e. , et al .
1998 , @xmath48492 , 696 .
scholer , m. , trattner , k.  j. & kucharek , h. 1992 , @xmath48395 , 675 .
shimada , n. & hoshino , m. 2000 , @xmath48543 , l67 .
shimada , n. , et al .
1999 , @xmath48264 , 481 .
smith , d.  a. et al .
2003 , @xmath48596 , l151 .
tavani , m. 1996a , @xmath4876 , 3478 .
tavani , m. 1996b , @xmath48466 , 768 .
terasawa , t. , et al .
2001 , in proc .
27th international cosmic ray conference , ( copernicus gesellschaft , germany ) , sh2.3/sh119 [ see _ http://www.copernicus.org/icrc/program.htm_ ] trattner , k.  j. & scholer , m. 1991 , @xmath4818 , 1817 .
trattner , k.  j. & scholer , m. 1993 , ann .
geophysicae @xmath4811 , 774 .
westfold , k.  c. 1959 , @xmath48130 , 241 .
zhang , b. & mszros , p. 2001
, @xmath48559 , 110 . _