the currently best - studied theories of gamma - ray bursts ( grbs ) and their afterglows ( ags ) are the _ fireball _ models ( see , e.g. , zhang & meszaros 2003 for a recent review ) and the _ cannonball _ ( cb ) model ( see , e.g. , dar & de rjula 2003a ; dado , dar & de rjula , 2002a ; 2003a and references therein ) .
the first set of models is often considered to be _ the standard model _ of grbs . in spite of their similarly - sounding names , these two models are ( or were initially ) completely different in their basic hypothesis , in their description of the data , and in their predictions . in this note
we concentrate on a cb - model prediction which is not ( to date ) a standard - model one , the apparently superluminal motion of the source of grbs and their afterglows : the  cannonballs " . in quite exceptional cases
relatively close - by grbs with sufficiently bright ( radio ) ags it may be possible to observe this superluminal motion of cbs relative to `` fixed stars '' .
prior to grb 030329 , the case in which a possible superluminal motion ( dar & de rjula 2000 ) came closest to being observable was that of grb 980425 , which could `` almost '' be classified ( dar & de rjula 2003a ) as an x - ray flash ( xrf ) .
for xrfs , the observation of a superluminal motion may be simpler than for grbs , for the source s apparent displacement in the sky is proportional to the ( small ) observer s viewing angle , and we interpret xrfs as jetted grbs viewed at larger angles ( dar & de rjula 2003a ; dado , dar & de rjula 2003f ) . in a different model ,
this interpretation of xrfs has also been advocated by yamazaki , yonetoku & nakamura ( 2003 ) .
the cb - model is a very explicit elaboration of the original proposal by shaviv and dar ( 1995 ) : that the @xmath0-rays of a grb would be generated by inverse compton scattering ( ics ) of stellar light by the electron constituents of transient narrow jets , emitted in stellar processes leading to gravitational collapse .    in the cb model , long - duration grbs , as well as xrfs ( dar & de rjula 2003a ;
dado , dar & de rjula 2003f ) are produced in the explosions of _ ordinary _ core - collapse sne , akin to sn1998bw ( dar & de rjula 2000 ) , the first sn to be observed in `` association '' with a grb ( grb 980425 , galama et al .
two opposite jets of cbs are emitted in the process , travelling with initially large lorentz factors : @xmath1 .
the cbs initially expand ( in their rest system ) at a velocity comparable to , or smaller than , the speed of sound in a relativistic plasma ( @xmath2 ) , so that the jet opening angle ( subtended by a cb s radius as observed from its emission point ) is @xmath3 .
the cbs highly relativistic motion collimates their emitted radiation the grb and its ag within a forward beam of characteristic opening angle @xmath4 .
an observer sees the `` doppler - favoured '' jet , travelling at a small angle @xmath5 relative to the line of sight .
typically @xmath6 , so that the jet s opening angle can be neglected and the observer s angle is the _ only _ relevant one .    in the cb model , during the ag , the jet opening angle _ diminishes _ with time , while the beaming angle @xmath7 increases , both evolutions being due to the cb s interaction with the matter of the interstellar medium ( dado , dar & de rjula , 2002a ) . as a consequence ,
the ag s source becomes increasingly `` pointlike '' , and its motion in the sky can in principle be followed . since a cb s motion is relativistic for days or months of ( observer s ) time ,
its apparent displacement in the sky is _ superluminal _ ( courdec 1939 , rees 1967 ) , as we argued in dar & de rjula ( 2000 ) .
the closest - by grb observed so far grb 980425 , at a redshift @xmath8 came very close to having an observable superluminal motion ( see dado , dar & de rjula 2003a for a detailed discussion ) . in a gcn note ( dar & de rjula 2003b )
we argued that this motion may be observable in the next - closest grb ( 030329 , at @xmath9 ) . in this note
we sharpen the predictions for this putative motion , given the current availability of precise optical data at early times and times later than the first @xmath10 days ( lipkin et al .
2003 and references therein ) , as well as sparse x - ray data ( marshall & swank 2003 ; marshall , markwardt & swank 2003 ; tiengo et al .
2003 ) and abundant radio data ( sheth et al .  2003 ; berger et al .  2003 ; pooley 2003 ; kuno et al .
2004 ) .    a very relevant new information ( lipkin et al .
2003 and references therein ) on the ag of grb 030329 is the abundance of multiple deviations of the optical light curves relative to a smoothly declining behaviour .
we shall refer to these deviations as _ features_. in dado et al .
( 2003c ) we attributed the most obvious optical - ag feature a  shoulder " starting at @xmath11 1 day ( after burst) to a transition between a first to a second dominant cb , a choice supported by the fact that the @xmath0-ray light - curve of this grb has a very clear two - pulse structure ( vanderspek et al .
2003 ; http:// space.mit.edu/ hete/ bursts/ grb030329/ ; see also vanderspeck et al .  2004 ) , as shown in fig .
( [ grb ] ) . with the emergence of a handful of other similarly - significant features in the first week fast ups and downs of the optical fluences , @xmath12 , by some 20 to 40% a more elaborate interpretation is required . in the ag model
@xmath12 is a direct and _ quasi - local _ tracer of the density of the ism through which a cb travels , a fact to which we have attributed previous similar observations , e.g.  the  humps " in the optical ags of grbs 000301c and 970508 ( dado et al .
it is therefore necessary to investigate the effect of local ism density - inhomogeneities on the expected superluminal motion .
this is what we do in this note for the case of grb 030329 .
we conclude that the sources motion results in a displacement of @xmath13 ( 0.6 ) mas from day 3 to day 30 ( 100 ) after burst . detecting such a motion may not be out of the question .    for the benefit of readers not familiar with the cb model , we present in an appendix a brief overview of the model and its current confrontation with data .
we also offer there some commentary on the evolution of the standard model .
the lorentz factor @xmath14 of a cb diminishes with the distance @xmath15 from the parent sn , as its motion is decelerated by collisions with the constituents of the interstellar medium ( ism ) . in an approximately hydrogenic ism of local number density @xmath16 , the evolution of @xmath17 is determined by energy - momentum conservation to be : @xmath18 with @xmath19 the cb s baryon number and @xmath20 its calculable asymptotic radius , reached within minutes of observer s time . for typical parameters , by the time @xmath21 , the function @xmath17 is still extremely well approximated by its initial value @xmath22 ( dado et al .
2002a ) .    for an object travelling at nearly the speed of light the relation between @xmath23 and the element of observer s time , @xmath24 ,
is : @xmath25 where @xmath26 is the cosmological redshift and @xmath27 is the doppler factor by which the energy of photons emitted by the cb is locally boosted by its relativistic motion . for the lorentz factors and viewing angles relevant to our discussion , @xmath28 and @xmath29 , and @xmath30 to an excellent approximation .
we may integrate eqs .
( [ dgamma1],[dgamma2 ] ) to obtain : @xmath31 which , for any given @xmath16 , yields the explicit @xmath15-dependent lorentz factor @xmath17 .
substitute this result , with use of eq .
( [ doppler ] ) , into eq .
( [ dxdt ] ) and integrate : @xmath32 to obtain the explicit relation between observer s time @xmath33 and travelled distance @xmath34 the last two equations are a parametric description of the function of actual interest , @xmath35 .
for some simple density profiles the above formal exercise can be carried explicitly to the end . for a constant - density ism ,
for instance , we may define : @xmath36 and conclude that @xmath35 is the real root of the cubic : @xmath37= { 2\,c\ , t\over 3\ , ( 1+z)\ , x_\infty}\ ; , \label{cubic0}\ ] ] that is : @xmath38\ , , \nonumber\\   c&\equiv & \left[{2/ \left(b^2 + 2\,\theta^6+b\,\sqrt{b^2 + 4\,\theta^6}\right)}\right]^{1/3}\ , ,
\nonumber\\   b&\equiv & { 1/ \gamma_0 ^ 3}+{3\,\theta^2/\gamma_0}+ { 6\,c\ , t/ [ ( 1+z)\ , x_\infty]}\ , , \label{cubic}\end{aligned}\ ] ] while the distance travelled by a cb at a given observer s time is : @xmath39\ ; , \label{xoft}\ ] ] so that it takes a distance @xmath40 for the cb to half its original lorentz factor . as we shall see in detail in the case of grb 030329 , by the time
the distance @xmath41 is reached , the ag fluence , proportional to a high power of @xmath35 , has decreased by more than one order of magnitude .
we have performed cb - model fits to the ag measurements of all the grbs of known redshift , and found that the values of @xmath41 are spread over more than an order of magnitude , centering at @xmath42 pc ( dado et al 2002a ) .
these are comparable to the radii of the _ superbubbles _ wherein most sn explosions take place , created by the stellar winds and explosions of many previous sne in dense star - formation regions .
the approximation of a constant density @xmath43 , employed in the quoted fits , may be a fair first try for the density within the superbubble and in the galactic - halo regions close - to but without the superbubble . in the case of grb 030329
we shall see that the density profile of the superbubble s surface , as traced by the cb s x - ray and optical fluences , is quite rich and interesting .
the transverse projected velocity in the sky of a cb relative to its parent sn can be obtained from eq .
( [ doppler ] ) : @xmath44 which , for typical parameters , is extremely superluminal .
the resulting angular separation at time @xmath33 is : @xmath45\ ; , \label{deltaalpha}\ ] ] where @xmath46 is the angular distance to the sn / cb system and @xmath47 is the luminosity distance ( we use throughout a cosmology with @xmath48 and @xmath49 ) .
the last expression in the r.h.s .  of eq .
( [ deltaalpha ] ) is valid for a constant - density ism .    at late time ,
when @xmath50 ^ 2\ll 1 $ ] , eq .
( [ doppler ] ) implies that @xmath51 , while eq .
( [ cubic ] ) implies that when @xmath52 ^ 2\ll   1 $ ] , @xmath35 approaches its asymptotic behaviour , @xmath53^{-1/3}$ ] .
the corresponding asymptotic behaviour of @xmath54 is , for a constant @xmath43 : @xmath55^{1/3}\ , , \label{latealpha}\ ] ] which , in practice , is a fair approximation for @xmath33 larger than a few days . since xrfs are grbs seen at a relatively large @xmath56 , eq .
( [ latealpha ] ) implies that the superluminal motion in the former systems is relatively easier to observe .
moreover , the larger the value of @xmath56 , the smaller the values of @xmath57 favoured by selection effects ( detection thresholds ) .
this also privileges xrfs over grbs as potential targets for a search of a superluminal cb motion .
in the cb model , the ags of grbs and xrfs consist of three contributions , from the cbs themselves , the concomitant sn , and the host galaxy : @xmath58 the latter contribution is usually determined by late - time observations , when the cb and sn contributions become negligible , or from measurements with sufficient angular resolution to tell apart @xmath59 from @xmath60 .    let the unattenuated energy flux density of sn1998bw at redshift @xmath61 ( galama et al .
1998 ) be @xmath62 $ ] . for a similar sn placed at a redshift @xmath26 ( dar 1999 , dar & de rjula 2000 , dado et al .
2002a ) : @xmath63 = { 1+z \over 1+z_{bw}}\ ; { d_l^2(z_{bw})\over d_l^2(z)}\ , a(\nu , z)\ , f_{bw}[\nu',t']\ , , \label{bw}\ ] ] where @xmath64 is the attenuation along the line of sight , @xmath65 , and @xmath66 .
the simple ansatz that _ all _ long - duration grbs would be associated with sn1998bw - like sne ( dar & plaga 1999 ; dar & de rjula 2000 ; dado et al .
2002a ) has proven to be unexpectedly precise and successful . for the most precise test so far , that of the grb 030329/sn2003dh pair , see , e.g. , dado et al .
( 2003c ) , stanek et al .
( 2003 ) , matheson et al .
( 2003 ) and hjorth et al .
( 2003 ) .
the ag of the cbs is mainly due to synchrotron radiation from accelerated electrons in the cb s chaotic magnetic field . at optical and higher frequencies ,
the ag has the approximate form ( dado et al .  2003a , c ) : @xmath67\propto
n_p^{(1+\hat\alpha)/2}\ , r_\infty^2\ ,   \gamma^{3\hat\alpha-1}\ , \delta^{3+\hat\alpha}\ , a(\nu , t)\ , \nu^{-\hat\alpha}\ , , \label{afterglow}\ ] ] with @xmath68 changing from @xmath69 to @xmath70 as each given observed frequency exceeds the time - dependent _ `` injection bend '' : _ @xmath71 ^ 3\ , \left[{n_p\over 10^{-3}\;cm^3}\right]^{1/2}\ , { \rm hz}\ , , \label{nubend}\ ] ] where @xmath43 is the baryon density of the interstellar medium changes from @xmath72 to @xmath73 ,
is interpolated by @xmath74^{1.1}}$ ] ( dado et al .
2003a ) . ] . in the same approximation
in which eq .
( [ latealpha ] ) was derived , the ags of grbs and xrfs have the asymptotic behaviour @xmath75 ( dado et al .
2002a ) . at radio frequencies ,
the ag spectrum is affected by self - absorption in the cbs themselves , characterizable by a single parameter per cb : a `` free - free '' absorption frequency , @xmath76 ( dado et al .
the attenuation @xmath77 is a product of the attenuation in the host galaxy , in the intergalactic medium , and in our galaxy .
the attenuation in our galaxy in the direction of the grb or xrf is usually estimated from the galactic maps of selective extinction , @xmath78 , of schlegel , finkbeiner & davis ( 1998 ) , using the extinction functions of cardelli et al .
the extinction in the host galaxy and the intergalactic medium , @xmath79 in eq .
( [ bw ] ) , can be estimated from the difference between the observed spectral index _ at very early time when the cbs are still near the sn _ and that expected in the absence of extinction .
indeed , the cb model predicts and the data confirm with precision the gradual evolution of the effective optical spectral index towards the constant value @xmath80 observed in all `` late '' ags ( dado et al .  2002a ; 2003a ) .
the `` late '' index is independent of the attenuation in the host galaxy , since at @xmath81 ( observer s ) days after the explosion , the cbs are typically already moving in the low - column - density , optically - transparent halo of the host galaxy .
the cb - model description of the r - band ags of six representative grbs of interest to the analysis of grb 030329 are shown in fig .
( [ tiriri ] ) , and discussed anon . in dado
et al .  ( 2002a ) we showed that , in the cb model , the ags of the grbs of known redshift measured at the time could be successfully explained by a single dominant cb launched in the explosion of a sn akin to sn1998bw and moving , after a couple of observer hours , in an approximately constant - density ism .
a typical example is that of grb 990510 .
two exceptions are grb 990123 and grb 021211 , whose ags were measured very early , at a time when their cbs may be piercing the expected @xmath82 density profile generated by the pre - sn `` wind '' of the parent star . in the cb model ,
the early fluence of optical ags is proportional to @xmath83 , in excellent agreement with the observations ( dado et al .
2003b ) .    prior to grb 030329
, there were three cases ( grbs 970508 , 000301c and 021004 ) for which the data showed clear evidence for deviations from a smooth ag behaviour . in the case of grb 000301c , we attributed the `` residua '' of the observational data relative to its smooth cb - model fit to moderate deviations of the ism density ( far from the parent sn ) from a constant value . in the case of grb 970508 , we have actually fit the ag to a jump from one to another constant - density value ( dado et al .
2002a ) . for grb 021004
we found it more natural to fit the ag to the contribution of two cbs , since in this case the @xmath0-ray count - rate as a function of time has a very clear two - pulse ( i.e.  two - cb ) structure ( dado et al .
2003d ) . moreover , unlike for grb 970508 , our attempts to describe the ag of grb 021004 in terms of density variations failed , while the two - cb description fits the general trend of the ag very well .
grb 030329 also has a two - cb @xmath0-ray structure , but an unprecedentedly non - smooth ag , presumably for the simple reason that the precision and continuity of the data are also unprecedented . in the case of this grb
, the question will arise whether or not we should interpret its peculiar ag shape to two cbs , to density fluctuations , or to a combination of both effects .
in dado et al .  ( 2003c ) we made a
first - round " cb - model fit to the nir - optical observations of the ag of this grb , then extending up to day @xmath10 after burst . since experience with
all previously - measured ags had given us confidence in the quality of the model , we extrapolated these fits to later times , and we predicted the presence of a sn akin to sn1998bw , luminous enough to compete with the ag and be discovered ten days after burst they have also reached the predicted asymptotic behaviour @xmath84 , and in the cb model there is no need to hypothesize whether `` breaks '' in the ag fluence have occurred by a given time , there being no breaks . ] .
sn2003dh was discovered 9.6 days after burst and turned out to be surprisingly similar to sn1998bw ( stanek el al .
2003 ; matheson et al .  2003 ; hjorth et al .
2003 ) .    in a second round
, we have now extended our original fit to include the x - ray data of rxte ( marshall & swank , 2003 ; marshall , markwardt & swank , 2003 ) and xmm - newton ( tiengo et al .
2003 ) , as well as the radio data of sheth et al .
( 2003 ) and berger et al .
( 2003 ) and many more nir - optical measurements ( lipkin et al .
2003 and references therein ) . the optical and x - ray data and the broad - band cb - model fit are shown in fig .
( [ figone ] ) .
the inclusion of the new data modifies the parameters of the original data - poor fit only at the 10% level , a satisfactory result .
the photon - number light - curve of the @xmath0 rays of grb 030329 consisted of two clear pulses , corresponding in the cb model to two dominant cbs , see fig .
( [ grb ] ) .
consequently , the fit to the ag was again performed with the additive contribution of two separate cbs . in the optical domain
, the two contributions correspond to the two shoulders , as illustrated in fig .
( [ figtwo ] ) , wherein we show the r - band results and the separate contributions of the two cbs .
the data variations relative to the predicted smooth ag light - curve , ups and downs of @xmath85 magnitude , are , as for grbs 970508 , 000301c ( dado et al .
2002a ) and 021004 ( dado et al .
2003b ) , to be expected : they trace moderate deviations from a constant - density interstellar medium , as implied by eqs .
( [ afterglow],[nubend ] ) .
these density variations are discussed in detail in the next section .
the radio data at different frequencies and their comparison with the cb - model fit are shown in fig .
( [ figthree ] ) .
the fit to the data at @xmath86 ghz is shown fig .
( [ figfour ] ) , wherein the contributions of the two cbs are separately shown .
as for the optical data , fig .
( [ figtwo ] ) , we conclude that the  late " ( @xmath87 d ) radio data are also dominated by one of the cbs .
the rich structure of the ag of this grb is shown in fig .
( [ figseven ] ) , in which lipkin et al .  ( 2003 ) have compared the r - band data to a  reference " ( a double power - law of indices @xmath88 and @xmath89 , with the break at day @xmath90 ) by plotting the magnitude difference between the data and the reference . the cb - model description of the optical and radio ags in figs .  ( [ figone ] ) and ( [ figthree ] ) satisfactorily reproduces the observed general trends .
but the description is unsatisfactory in that its residua are also significant , as we show in fig .
( [ figten ] ) .
these residua are not unlike those in fig .
( [ figseven ] ) , but for the absence of a prominent  @xmath91 " feature , which we view as the result of the comparison of a power law unmotivated at early times with the smoothly - varying data .
the remaining features , particularly the more prominent ones occurring after @xmath92 day , require an explanation .
the cb - model s fluence of eq .
( [ afterglow ] ) is proportional to a power of the ism s electron density @xmath93 and , in that sense , it
 traces " its local , instantaneous value
. this tracing is not perfectly _ local _ , because @xmath94 depends on @xmath35 , and this function reflects the integral effect of the ism density that a cb has swept through its prior voyage . to reproduce the features still prominent in fig .
( [ figten ] ) , we must assume a given non - trivial density profile and solve eqs .
( [ dgamma2 ] ) to ( [ tint ] ) explicitly to obtain the function @xmath35 to be used as input in the expression for the fluence , eq .
( [ afterglow ] ) .
we have in the past done this for the ag of grb 970508 , as reproduced in fig .
( [ tiriri ] ) .
the cbs of grbs 970508 and 030329 , as we shall see explicitly for the second one , are at a distance @xmath95 pc from the sne that emitted them , at an observer s time @xmath92 day . at a distance of that order
, we expect the cbs to be exiting the superbubble where the sn explosion is very likely to have taken place .
as they exit it , they may encounter successive density inhomogeneities produced by a succession of past sn explosions and stellar mass ejections , which created the bubble in the first place .
we have considered the _ structured density profile _ shown ( as a function of observer s time ) in fig .
( [ figeight ] ) , in which successive  onion peals " appear as abrupt density increases , followed by a decline proportional to @xmath96 ( or , approximately , @xmath97 ) , analogous to the closer - by profile of a
wind - fed " circumburst material .
with such a profile of overdensities we obtain a description of the r - band ag whose residua ( relative to the reference broken power - law introduced by the observers ) are shown in fig .
( [ figseven ] ) as the ( red ) line .
the parameters of this fit are only slightly different from the ones for a constant @xmath98 , but for a somewhat larger initial value of @xmath99 for the cb dominating the ag at late times .
the @xmath91 feature is fairly well reproduced with the density profile of fig .
( [ figeight ] ) , which is constant up to day @xmath100 .
the feature is , in this sense , a fake .
the description of the later - time features , however , requires the structured density profile of fig .
( [ figeight ] ) .
clearly the many - parameter exercise we just described is not a  fit " : it could easily be improved to obtain an even better description , and it is certainly not as  predictive " as our prior fits to the early part of the ags of grbs 990123 and 021211 , from which we claimed to have successfully traced the magnitude and @xmath96 profile of the close - by circumburst density ( dado et al .
yet , the description is phenomenologically satisfactory , in the sense that we can once again claim that , in the cb model , the shapes of ags provide interesting and consistent _ quasi - local _ tracers of the ism density through which the cbs travel .
the actual distance @xmath101 from the sn at which the late - ag - dominating cb in grb 030329 was , as seen at a given observer s time @xmath33 , is shown in fig .
( [ figkpc ] ) , constructed with use of eq .
( [ tint ] ) . indeed , at the time at which the density inhomogeneities are apparent ( @xmath33 between @xmath102 and @xmath90 days ) the cb is 100 to 200 pc away from the sn , a reasonable radius for the overdensities surrounding a superbubble .
given the need to introduce density variations to describe the detailed behaviour of the ag , the question arises whether or not two cbs as opposed to just one are needed .
our attempts to describe the ag with just one cb and a variety of density profiles failed .
the reason is simple and unavoidable : the density increase required to produce the ag s  shoulder " at @xmath87 day inescapably increases the rate at which @xmath35 decreases with time .
this entails a fast fall - down of the fluence , following soon after the initial increase produced by the local density increase .
the resulting light - curves rise and fall fast like the one shown in fig .
( [ tiriri ] ) does in the case of grb 970508 and they fail in the description of the ag of grb 030329 .
thus , two cbs appear to be necessary to understand the ag of grb 030329 , as they are for the description of its two - pulse grb structure , shown in fig .
( [ grb ] ) .
we are primarily concerned in this note with the putative superluminal motion of the cbs of grb 030329 , and not with a complete description of its wide - band ag .
thus , we have limited our detailed analysis with a structured density profile to the r - band ag , as in figs .
( [ figseven ] ) , the complete broad - bend analysis with density variations being very laborious .
the outcome of the discussion that follows is that the structure seen in the r - band ag should be _ almost _ achromatic for data ranging from x - ray to nir frequencies , while the structures should progressively disappear at lower and lower radio frequencies .
the x - ray and optical ags are not perfectly achromatic because the injection - bend frequency of eq .
( [ nubend ] ) `` crosses '' a given frequency at different times .
but the crossing of the optical frequencies typically occurs at @xmath92 day and the light curves become increasingly achromatic thereafter .
the radio ags , on the other hand , are predicted to be strongly chromatic ( dado et al .
2003a ) . to be more specific
, we quote from the mentioned reference :    _  electrons that enter a cb with an injection lorentz factor @xmath35 are rapidly fermi accelerated ... on a longer time scale , they lose energy by synchrotron radiation , and their energy distribution evolves ... electrons with a large @xmath103 $ ] emit synchrotron radiation , with no significant time - delay , at the observer s optical and x - ray wavelengths . but the emission of radio is delayed by the time it takes the electrons to `` descend '' to an energy at which their characteristic emission is in the observer s radio band ... thus , the optical and x - ray ag emission starts ... a few observer s seconds after the corresponding @xmath0-ray pulse .
the radio signal , on the other hand , must await a time @xmath104 _
_ for the cumulated electrons to cool down . "
@xmath105 was explicitly estimated in dado et al .
( 2003a ) to be of @xmath106 day .    applied to the case at hand of a varying ism density , what this means is that the features observed in the optical ags should be smoothed in the radio over intervals of @xmath106 day .
these intervals are frequency - dependent : longer for the lower radio frequencies .
the expected trend of features that progressively disappear as the frequency diminishes is precisely the trend observed in the radio data reproduced in fig .
( [ figthree ] ) .
the `` second - round '' fit shown in the figure systematically overestimates the late results , also a consequence of having ignored the density enhancements .
the location of the source of the ag of this grb has been followed in detail by vlba observers from day @xmath107 3 to day @xmath107 84 after burst ( taylor et al .
, to be published ) .
the associated sn2003dh dominates the optical ag after day @xmath108 , but is negligible in the radio ag at all measured times .
thus , the location of the radio source(s ) is the location of the cb(s ) , expected to be extremely superluminal in the cb model .
the values of the fit observer s angles are very similar for the two cbs : @xmath109=2.2 $ ] mrad , @xmath110=2.3 $ ] mrad .
the values of the initial lorentz factors are not so different : @xmath111=1037 $ ] , @xmath112=1606 $ ] .
but the values of the deceleration parameter are very different : @xmath113=0.033 $ ] mpc , @xmath114=0.37 $ ] mpc ( these numbers refer to the initial @xmath99 as in eq .
( [ xinf ] ) , still unaffected by the density variations occurring at @xmath87 day ) .
this is the main reason why the contribution of cb1 to the ag decays much faster with time than that of cb2 , as implied by eq .
( [ cubic0 ] ) and illustrated in figs .
( [ figtwo],[figfour ] ) .
it is also the reason why the superluminal motion of cb1 is much slower than that of cb2 , as implied by eq .
( [ deltaalpha ] ) .
we are therefore interested in the fastest - moving and late - ag - dominating cannonball : cb2 .    in the absence of the density fluctuations of fig .
( [ figeight ] ) , the predicted angular displacement @xmath115 of cb2 , as given by the integrated form of eq .
( [ deltaalpha ] ) , would be that of the upper panel of fig .
( [ tororo ] ) .
it is on the basis of our  first - round " constant - density fit that we predicted such a displacement in dar & de rjula ( 2003b ) . in the presence of the observed density fluctuations ,
the predicted @xmath115 is shown in the lower panel of fig .
( [ tororo ] ) , the result of the integration of the velocity @xmath116 in eqs .
( [ supervelocity],[deltaalpha ] ) , obtained from the full - fledged determination of @xmath35 from eqs .
( [ dgamma1 ] ) to ( [ tint ] ) .
the density fluctuations are dominantly density enhancements , implying a faster deceleration of the cb .
thus the significantly reduction in the predicted apparent superluminal motion .
alas , this prediction of the cb - model is , in the case of grb 030329 , harder to test than we originally thought .
in its current form , the theory of ags in the cb model is based on an analogy with high - resolution observations of relativistic jets emitted by systems such as the microquasars ss 433 ( http://chandra.harvard.edu/ photo/2002/ 0214/ index.html ; http:// chandra .
edu/ press/04@xmath117releases/ press@xmath117010504.html ) , grs 1915 + 105 ( mirabel & rodriguez 1999 ; dhawan et al .  2000 ) and xte j1550 - 564 ( corbel et al .
2002 ) ; and active galactic nuclei , e.g.  the quasar pictor a , whose emitted cbs , as seen at x - ray wavelengths , seem to stop expanding laterally shortly after ejection , and to travel for hundreds of kpc , before they finally stop and blow up ( wilson , young & shopbell 2001 ; grandi et al .
the cb - model s simple explanation for this surprising fact is the following ( dado et al .
the ambient protons intercepted by a cb in its voyage encounter its inner chaotic magnetic field and are reemitted roughly isotropically ( in the cb s rest system ) .
this implies an inwards pressure that stops the cb s expansion and that is asymptotically equal to the pressure of the cb s inner magnetic field , which is thereby calculable .
the corresponding deceleration law ( at constant radius @xmath20 and ism density @xmath43 ) is that of eq .
( [ cubic0 ] ) .
the magnitude and time - dependence of the magnetic field determined as in the previous paragraph plays a crucial role in the prediction of the broad - band ag spectra , which in the cb model is extremely simple , parameter - thrifty , and successful ( dado et al .
2003a , 2003e ) , particularly in comparison with its standard counterpart ( see , e.g.  granot & sari 2002 ) . in spite of their phenomenological success , our assumptions leading to calculable cb radii and magnetic fields are very bold , and are no doubt oversimplifications . to investigate their solidity , in dado et al .
( 2002a ) we also studied a different ansatz : that the cb s radius @xmath118 would keep on increasing constantly with @xmath23 , with the ism density still kept approximately constant .
this ansatz , which corresponds to a much faster cb deceleration than that of eq .
( [ cubic0 ] ) , failed miserably .
but in between the two extremes ( a non - increasing and a steadily increasing cb s radius ) , it is clear that there should be choices of @xmath119 and @xmath16 that are phenomenologically satisfactory and yet , correspond to a _
faster _ cb deceleration than the simple one dictated by eq .  ( [ cubic0 ] ) at constant @xmath120 . for such choices the simplest predictions of superluminal velocities
that we have discussed would result in overestimates .
but it would be premature to modify the simple predictions of the cb model , which so far are succesful .
in the internal / external shock model of grbs , the variability of ag light curves the set of  features " that we have discussed is attributed to patchy shells and refreshed shocks ( see , e.g.  piran , nakar & granot 2003 ; granot , nakar & piran 2003 and references therein ) . in the case of grb 030329
, the features are attributed to delayed collisions between late , relatively low-@xmath0 shells and the earlier higher-@xmath0 shells having decelerated in their interactions with the ism . the timing and the lorentz factors
are chosen so that the early fast shells collide in a matter of ( observer ) seconds to produce the grb peaks , while the slow shells catch up with the decelerated earlier ones at times of order days : some five order of magnitudes longer .
the required amount of fine - tuning is considerable .
the late shells bring in an energy injection an order of magnitude larger than that of the original blast wave
. it would be interesting to know whether their collisions with the slowed - down shells are expected to give rise to grb pulses as well , in which case grbs should be  repeaters " .    in the cb model , as we have seen , all that is required to explain the ag  features " are ism density inhomogeneities occurring at a natural distance from the parent sn : that corresponding to the radius of a superbubble .
no new phenomena must be invoked and no grb parameters must be fine - tuned with special care .
the data on grb 030329 are now sufficiently complete to allow for a detailed prediction of the motion of the source of its ag allegedly superluminal in the cb model .
the @xmath0-ray light - curve and the optical ag require the presence of two cbs , one of which dominates the ag at late times .
the parameters needed to predict the motion of the two cbs in the sky are determined by the optical data , so that the individual motion of each cb is predicted .
the main result of this paper is the prediction of the sky - motion of the fastest - moving cb , which dominates the late ag , and is shown in the lower panel of fig .
( [ tororo ] ) .
the main caveat concerning a putative superluminal signature ( dar & de rjula 2000 , 2003b ) concerns the input to the estimate of its magnitude .
indeed , we have shown in this paper that the predictions are very sensitive to the details of the density profile of the ism .
but we are encouraged by the fact that in the only case in which the superluminal jets of cbs made by a core - collapse sn could be seen , they were seen .
indeed , observations of sn 1987a ( ninenson and papaliolios , 1999 ) showed two sources , moving in opposite directions along the sn s axis at ( real ) velocities compatible with the speed of light and at an apparent superluminal velocity for the approaching source . mercifully , the jet of that sn was not pointing in our direction ( dar , laor & shaviv 1998 ; dar & de rjula 2001b ) .    regarding the search for a superluminal motion , we learned by reading the e - version of nyt 030529 ( the new york times of that date , in grb s parlance ) that , according to dale frail _  [ our observations ] are sufficient to rule out predictions of the cannonball model"_. we have shown that , indeed , the observations of complicated features in the optical ag of grb 030329 imply that our earlier results ( dar & de rjula 2003b ) which ignored the presence of these features constituted an overestimate of the predicted superluminal displacement . in this sense ,
frail was right in stating that the observations ruled out the _ predictions _ , as opposed to the model itself .    in a setting more scientific than the nyt , bloom et al .
( 2003 ) state : _  owing to the proximity and bright radio emission , high - resolution ( @xmath121 pc ) _ very long baseline array _ imaging of the compact afterglow was used by frail ( 2003 ) to _ * unequivocally disprove * _ the cannonball model for the origin of grbs . "
_ the emphasis is ours . we have seen that these news of the death of the cb model may have been premature . even though mark twain eventually died for sure
, the cb model though probably not immortal is still in an excellent shape . yet
, trying to disprove the best available model(s ) , or even the proof of a difficult theorem , is the acceptable standard attitude in many a realm of the exact sciences . in this sense ,
the apparently strong motivation of the quoted observers to disprove the cb model is in our opinion , and regarding this particular model the healthiest of all possible attitudes .
* acknowledgements .
* one of us ( a.  de r. ) is indebted to the physics department and space research institute of technion for its hospitality .
this research was supported in part by the helen asher space research fund for research at the technion .
in this model , the _ engines _ generating long - duration grbs and xrfs are ordinary core - collapse sne . following the collapse of the stellar core into a neutron star or a black hole , and given the characteristically large specific angular momentum of stars , an accretion disk or torus is hypothesized to be produced around the newly - formed compact object , either by stellar material originally close to the surface of the imploding core and left behind by the explosion - generating outgoing shock , or by more distant stellar matter falling back after its passage ( de rjula 1987 ) .
a cb is emitted , as observed in microquasars , when part of the accretion disk falls abruptly onto the compact object ( e.g.  mirabel & rodrigez 1999 ; rodriguez & mirabel 1999 and references therein ) .
the cbs of a grb are assumed to me made of ordinary hydrogenic matter .
this is also in analogy to microquasar cbs , observed to be made of ordinary matter , as opposed to the standard - model contrived mixture of baryons and @xmath122 pairs .
indeed , blueshifted and redshifted h , he , metal and heavy - element optical and uv lines were detected from the approaching and receding cbs emitted by the microquasar ss433 ; see , e.g.  eikenberry et al .
( 2001 ) , gies et al .
( 2002 ) .    in the cb model , sn1998bw , associated with grb 980425 , is an ordinary core - collapse sn ( dar and plaga 1999 ; dar and de rjula 2000 ) : its `` peculiar '' x - ray and radio emissions were not emitted by the sn , but were part of the grb s ag ( dado et al .
2002a , 2003a ) .
the grb appeared to be underluminous because it was viewed further off axis than others ( e.g. , dar and de rjula 2000 ) .
thus , it makes sense in this model to expect an association between ( long ) grbs and sn1998bw - like engines , transported to the grb location ( dar 1999 ) .
the a - posteriori analysis of dado et al .
( 2002a ) of all the available data supported the conclusion that _ all _ ( long ) grbs are indeed associated with such sne .
in several more recent cases we used the cb model to _ predict _ how the associated sn would compare with the rest of the ag ( grb 011121 , 020405 , 021211 ; dado et al .  2002b , 2002c , 2003d ) .
the most interesting case is that of grb 030329/sn2003dh , for which we also foretold the date when the sn would be dominant enough to be discovered ( dado et al .
2003c ) .
the fraction of visible grbs , relative to the sne that produce them , is @xmath123 , where @xmath124 is some effective value of the maximum observer s angle for which grbs have been detectable . in practice
@xmath124 is a few milliradians ( dado , dar & de rjula 2002a ) .
this brings the _ total _ rate of ( observed or unobserved ) grbs to a value close to that of core collapse sne .
this implies that a good fraction of such sne may produce long - duration grbs ( dar & plaga 1999 ; dar & de rjula 2000 ) .
in dar & de rjula ( 2000 ) the initial lorentz factor of the cbs was argued to be @xmath1 , a choice corroborated by all our subsequent analyses .
the high - energy photons of a single pulse in a grb or an xrf are produced as a cb coasts through the `` ambient light '' permeating the surroundings of the parent sn .
the electrons enclosed in the cb compton up - scatter photons to energies that , close to the cbs direction of motion , correspond to the @xmath0-rays of a grb and less close to it , to the x - rays of an xrf .
each pulse of a grb or an xrf corresponds to one cb .
the timing sequence of emission of the successive individual pulses ( or cbs ) reflects the chaotic accretion process and its properties are not predictable , but those of the single pulses are ( dar & de rjula 2003a and references therein ) . in practice grbs are observable for @xmath125 , xrfs are the same phenomenon , observed from somewhat larger angles ( dar & de rjula 2003a ; dado , dar & de rjula 2003f )
.    grbs are notorious for their variety .
yet , they have some two dozen common properties .
the characteristic @xmath0-ray energy is surprisingly narrowly distributed around @xmath126 kev .
the energy spectra are well fit by a specific  band " function .
the total `` equivalent spherical energies '' are distributed over a wide range , that becomes very narrow for the `` true '' ( beaming - corrected ) energy values .
the photons arrive in successive pulses with fairly similar time - profiles which , in  long " grbs , have a median width of 0.5 seconds .
the @xmath0-ray polarization is nearly maximal .
a long list of phenomenologically - observed correlations between many different grb observables appear to be well satisfied . _
all _ of these properties are very simply and successfully explained by the cb model , in a manner that does not require the introduction of _ any _ ad - hoc parameters ( dar & de rjula 2003a and references therein ) .      in the ag phase , and within a few observer minutes from the end of a grb s intense @xmath0 emission , the emissivity of a cb is dominated by synchrotron emission from the electrons that penetrate in it as it propagates in the ism . integrated over frequency ,
this synchrotron emissivity is proportional to the energy - deposition rate of the ism electrons in the cb .
these electrons are  fermi - accelerated " in the cb s tangled magnetic maze to a broken power - law energy distribution with a `` bend '' energy equal to their incident energy in the cbs rest frame , @xmath127 .
their synchrotron radiation the afterglow is collimated and doppler - boosted by the relativistic motion of the cbs . the radiation is also redshifted by the cosmological expansion and attenuated on its way to an earthly observer , during its passage through the cb itself , the host galaxy , the intergalactic space and our own galaxy .
the approximate cb - model analysis of the fluences of grb ags , as functions of time and frequency , is fairly simple . besides @xmath26 and @xmath56 ( which are not parameters specific to the model ) it involves only three quantities : @xmath128 , @xmath99 , an overall normalization , and @xmath76 ( a _ single _ free - parameter frequency pertinent to the radio - absorption within a cb ) . in spite of its economy of parameters ( notorious in the case of the wide - band spectral shapes )
the cb model satisfactorily describes in a unified manner all measured grb ags , including that of grb 980425 ( e.g. , dado et al .
2002a , 2003a , 2003e ) .
our impression is that no such claims can be made in the realm of standard grb models .
berger , e. , et al .
2003 , nature , 426 , 154 blake , c. & bloom , j. s. 2003 , gcn circ .
2011 bloom , j. s. , et al .
2003 , astro - ph/0308034 cardelli , j. a. , clayton , g. c. & mathis , j. s. 1989 , apj , 345 , 245 corbel , s. , et al . 2002 ,
science , 298 , 196 courdec , p. 1939
, annales dastrophysique 2 , 271 dado , s. , dar a. & de rjula , a. 2002a , a&a , 388 , 1079 dado , s. , dar a. & de rjula a. 2002b , apj , 572 , l143 dado , s. , dar a. & de rjula a. 2002c , a&a , 393 , l25 dado , s. , dar a. & de rjula a. 2003a , a&a , 401 , 243 dado , s. , dar a. & de rjula a. 2003b , apj , 585 , l15 dado , s. , dar a. & de rjula a. 2003c , apj , 594 , l89 dado , s. , dar a. & de rjula a. 2003d , apj , 593 , 961 dado , s. , dar a. & de rjula a. 2003e , phys . lett .
b562 , 161 dado s. , dar a. & de rjula a. 2003f , astro - ph/0309294 dahawan , v. , mirabel , f. & rodriguez , l. f. 2000 , apj , 543 , 373 dar , a. & de rjula , a. 2000 , astro - ph/0008474 dar , a. & de rjula , a. 2001 , in _ astrophysics and gamma ray physics in space _ ; eds .
a. morselli and p. picozza ; frascatti physics series vol .
xxiv ( 2002 ) , pp .
513 - 523 ; astro - ph/0110162 dar , a. & de rjula , a. 2003a , astro - ph/0308248 dar , a. & de rjula , a. 2003b , gcn circ
2133 dar , a. , laor , a. & shaviv , n. j. 1998 , phys .
80 , 5813 dar , a. , & plaga , r. 1999 , a&a , 349 , 259 dar , a. 1999 , a&as , 138 , 505 dar , a. 2003 , astro - ph/0301389 de rjula , a. 1987 , phys .
193 , 514 de rjula , a. 2002 , astro - ph/0207033 donaghy , t. q. , lamb , d. q. & graziani .
c , 2003 , astro - ph/031496 eikenberry , s. s. , et al .
2001 , apj , 561 , 1027 esin , a. a. & blandford , r. 2000 , apj , 534 , l151 frail , d. a. , et al . , 2001 ,
apj , 562 , l55 frail , d. a. 2003 , american astronomical society meeting , 2003 frederiksen , j. t. , hedelal , c. b. , haugbolle , t. & nordlund , a. astro - ph/0303360 galama , t. j. , et al .
1998 , nature 395 , 670 ghisellini , g. & lazzati , d. 1999 , astro - ph/9906471 gies , d. r. , et al .
2002 , apj , 566 , 1069 grandi , p. , et al .
2003 , 2003 , apj , 586 , 123 granot , j. , nakar , e. , piran , t. 2003 , nature 426 , 138 granot , j. & sari , r. 2002 , apj .
568 , 820 granot , j. , panaitescu , a. , kumar , p. & woosley , s. e. 2002 , apj , 570 , l61 greiner , j. , et al . 2003a , gcn circ .
2020 greiner , j. , et al .
2003b , astro - ph/0311282 hjorth , j. , et al .
nature 2003 , 423 , 847 hurley , k. , sari , r. & djorgovski , s. g. 2002 , astro - ph/0211620 henden , a. a. , et al .
2003 , gcn circ .
2082 kuno , n. , et al .
2004 , astro - ph/0401258 lamb , d. q. , donaghy , t. q. & graziani , c. 2003a , astro - ph/0309456 lamb , d. q. , donaghy , t. q. & graziani , c. 2003b , astro - ph/0309463 lamb , d. q. , donaghy , t. q. & graziani , c. 2003c , astro - ph/030504 lazzati , d. , rossi , e. , ghisellini , g. & rees , m. 2003 , astro - ph/0309038 , lazzati , d. , ghisellini , g. , celotti , a. & rees , m. 2000 , apj , 529 , l17 lazzati , d. 2003 , astro - ph/0312331 li , w. , filippenko , a.v .
, chornock , r. & jha , s. 2003 , apj , 586 , l9 lipkin , h. y. , et al .
2003 , astro - ph/0312594 margon , b. a. 1984 , ara&a , 22 , 507 marshall , f. e. & swank , j. h. 2003 , gcn circ .
1996 marshall , f. e. , markwardt , c. & swank , j. h. 2003 , gcn circ .
2052 matheson , t. , et al .
apj , 2003 , 599 , 394 meszaros , p. & rees , m. j. 1999 , mnras , 306 , l39 mirabel , f. & rodriguez , l. f. 1999 , ara&a , 37 , 409 nisenson , p. & papaliolios , c. 1999 , apj , 518 , l29 patat , f. , et al .
2001 , apj , 555 , 900 peterson , b. a. & price , p. a. 2003 , gcn circ .
1985 piran , t. , nakar , e. & granot , j. 2003 , astro - ph/0312138 pooley , g. 2003 , gcn circs . 2043 , 2072 ramirez - ruiz , e. , dray , l. m. , madau , p. & tout , c. a. 2001 , mnras , 327 , 829 rees , m. j. 1967 , mnras , 135 , 345 rees , m. j. & meszaros , p. 1992 , mnras , 258 , l41 reichart , d. e. , 2001 , apj , 554 , 643 rodriguez , l. f. & mirabel , i. f. 1999 , apj , 511 , 398 rykoff , e. s. , et al .
2003 , gcn circ .
1995 sari , r. & piran , t. 1999 , a&as , 138 , 537 sari , r. 1999 , apj , 524 , l43 shaviv , n. j. & dar , a. 1995 , apj , 447 , 863 sheth , k. , et al .
2003 , apj , 595 , l33 schlegel , d. j. , finkbeiner d. p. & davis m. 1998 , apj , 500 , 525 stanek , k. z. , et al .
2003 , apj , 591l , 17 tiengo , a. , mereghetti , s. , ghisellini , g. , rossi , e. , ghirlanda , g. & schartel , n. 2003 , astro - ph/0305564 torii , k. 2003 , gcn circ .
1986 vanderspek , v. , et al .
2003 , gcn circ .
1997 vanderspeck , r. , et al .
2004 , astro - ph/0401311 waxman , e. 2003 , astro - ph/030351 waxman , e. & draine , b. t. 2000 , apj , 537 , 796 wilson , a. s. , young , a. j. & shopbell , p. l. 2001 , apj , 547 , 740 yamazaki , r. , yonetoku , d. & nakamura , t. apj , 2003 , 594 , l79 zeh , a. , klose , s. & hartmann , d. h. astro - ph/0311610 zhang , b. & meszaros , p. 2003 , astro - ph/0311321


recent advances in the construction of low cost , low power , and mass produced micro sensors and micro - electro - mechanical ( mem ) systems have ushered in a new era in system design using distributed sensor networks @xcite .
the advent of sensor network technology provides a variety of applications that have been considered unrealistic in the past .
one such application is in the area of space communications : with ad hoc sensor networks , a number of sensor nodes randomly placed on a planet can collaboratively collect information and then , also collaboratively , send the information back to earth . in this scenario
, the sensors must have an ability to transmit information over very long distances with high energy efficiency . in this kind of point - to - point communication ,
directional antennas are a preferred means to avoid interference .    in general ,
this can be achieved by adaptive beamforming . given a number of well - designed antenna elements at the transmitting /
receiving sensor nodes , each node could in principle autonomously transmit / receive the information to / from any desired direction .
the advantages and applications of beamforming with antenna arrays are well known ; in wireless communications , this enables space - division multiplex access ( sdma ) , a technology which has the potential to significantly increase the capacity of multiple access channels .
one of the most important constraints on wireless sensors is energy efficiency .
since the sensor nodes are often distributed in places where manual maintenance is costly , such as remote locations , on top of buildings and so on , it should be possible to operate these for several months without battery replacement . considering the fact that each antenna element requires analog circuitry ( and thus leads to costly hardware ) , in practice each distributed sensor is likely to be equipped with only a single antenna and this precludes the use of autonomous beamforming in scenarios of very energy efficient communication .
nevertheless , if the sensors in the cluster share the information _ a priori _ and synchronously transmit the data collaboratively as sketched in fig.[fig : adhoc_network_model ] , it may be possible to beamform when transmitting ( and also receiving ) the data in a distributed manner
. the resultant overhead due to intra - cluster information sharing may be relatively small as this can be done by low - cost short distance broadcasting - type communication among nodes .
thus , with distributed collaborative beamforming , the nodes can send the collected information to the far - end receiver over long distances with high efficiency .
also , only the sensor cluster in the specified target direction receives the data with high signal power and no significant interference occurs for clusters in other directions .
overall there is thus a potential to increase the capacity of the multiple access channel significantly despite the additional overhead for information sharing .
the obvious question is whether one can form a nice beampattern with a narrow mainbeam , or achieve a reasonable directional gain ( directivity ) .
the sensor nodes in ad hoc networks are located randomly by nature , and the resultant beampattern depends on the particular realization of the sensor node locations .
therefore , it may be quite natural to treat the beampattern with probabilistic arguments . in this paper
, assuming idealized channel model conditions and antenna properties , we analyze the achievable performance of collaborative beamforming based on distributed sensor nodes in a probabilistic sense . specifically , the statistical properties of the achievable beampattern of the random sensor arrays are studied based on the following assumptions .
the sensors form an ad hoc network and the geometry of the cluster is given by a two - dimensional disk of a given radius over which all sensor nodes are distributed uniformly as illustrated in fig.[fig : adhoc_network_model ] .
since the corresponding far - field beampattern depends on the particular realization of the random array of nodes , the probability distribution of the far - field beampattern is of particular interest .    to the best of the authors knowledge ,
the beampattern aspects of collaborative beamforming using random arrays have not been analyzed before in the framework of wireless ad hoc sensor networks .
nevertheless , in the antenna design literature , probabilistic analysis of random arrays is not new . in the framework of linear array design with a large number of sensors ,
lo @xcite has developed a comprehensive theory of random arrays in the late 1960s .
it has been shown that randomly generated linear arrays with a large number of nodes can in fact form a good beampattern with high probability and that with linear random arrays of @xmath0 sensors , the directivity approaches @xmath0 asymptotically .
although our scenario is quite different in that our main goal is not to design array geometry but to _ exploit _ the randomness of the distributed sensor network , it turns out that the results we shall develop in this paper can be seen as an extension of the theory of linear random arrays , @xcite , to random arrays on a disk .
thus , the same conclusion will be reached : with @xmath0 collaborative sensor nodes , one can achieve a directivity of order @xmath0 asymptotically .    the major difference between classical beamforming by antenna arrays and distributed beamforming
is that whereas the geometry of the former is usually known _ a priori _ , the exact location of the sensor nodes in ad hoc network is not , and it should be acquired dynamically . even if their relative location is estimated by some adaptive algorithm ( e.g. , @xcite for receiver beamforming ) , considering the low snr operation of the sensor nodes , it is almost certain that the acquired geometric information has some inaccuracy . also , since all nodes are operated with physically different local oscillators , each node may suffer from statistically independent phase offsets . in order to model and
elucidate the effect of these impairments , we consider the following two scenarios : _ closed - loop _ and _ open - loop_. the closed - loop scenario may be described as follows .
each node independently synchronizes itself to the beacon sent from the destination node ( such as a base station ) and adjusts its initial phase to it .
thus , the beam will be formed in the direction of arrival of the beacon .
this kind of system is often referred to as a _ self - phasing _ array in the literature , and may be effective for systems operating in time - division duplex ( tdd ) mode .
the residual phase jitter due to synchronization and phase offset estimation among sensor nodes is then often the dominant impairment . on the other hand , in the open - loop scenario
we assume that all nodes within the cluster acquire their relative locations from the beacon of a nearby reference point or cluster head .
the beam will then be steered toward an arbitrary direction .
thus , the destination need not transmit a beacon , but each node requires knowledge of its relative position from a predetermined reference point within the cluster .
this case may occur in ad hoc sensor networks where sensor nodes do not have sufficient knowledge of the destination direction _ a priori_. in this scenario ,
since the acquisition of precise knowledge is not realistic , the effects of location estimation ambiguity among sensors upon the beampattern may be of particular interest .    throughout the paper ,
the nodes and channel are assumed to be static over the communication period , and for simplicity the information rate is assumed to be sufficiently low that inter - symbol interference ( isi ) , due to residual timing offset , is negligible .
it will also be assumed that all nodes share the same transmitting information _ a priori _ , as the main focus of the paper is on the beampattern , rather than the front - end communication performance .
the paper is organized as follows .
section [ sec : model ] describes the assumptions , model , and main parameters that describe beam characteristics associated with the framework of wireless ad hoc sensor networks . in section [ sec : average_property ] , the average properties of the beampattern are derived .
the average beampattern of linear random arrays has been derived in @xcite , and our results can be seen as its extension to our sensor network model . for analytical purpose , we also introduce the concept of 3 db sidelobe region . in section [ sec : distribution_beam ] , the statistical distribution of the beampattern in a specific direction is derived .
lo @xcite has derived the distribution of the beampattern in linear arrays based on a gaussian approximation of the array factor , which is a common assumption in the random array literature .
in contrast , we shall develop a numerical method to calculate the exact distribution of the beampattern and also examine the accuracy of the gaussian approximation in detail .
the distribution of the maximum of the sidelobe region is discussed in section [ sec : distribution_maximum ] .
this aspect of beampattern was analyzed by steinberg @xcite , agrawal and lo @xcite , and donvito and kassam @xcite in the framework of linear random arrays . in this paper
, we derive an upper bound on the distribution of the maximum sidelobe in our framework of collaborative beamforming based on the approach of @xcite .
the effect of phase jitter or location estimation errors on the resultant beampattern , associated with the closed - loop and open - loop scenarios , is analyzed in section [ sec : imperfect ] .
steinberg @xcite has analyzed the effects of phase estimation errors in linear arrays , and based on a similar approach we analyze the effects of the average beampattern with phase estimation errors .
finally , section [ sec : conclusion ] concludes the paper .
[ bl][bl][.8]@xmath1 [ bl][bl][.8]@xmath2 [ bl][bl][.8]@xmath3 [ bl][bl][.8]@xmath4 [ bl][bl][.8]@xmath5 [ bl][bl][.8]@xmath6 [ bl][bl][.8]@xmath7 [ bl][bl][.8]@xmath8 [ bl][bl][.8]@xmath9     the geometrical configuration of the distributed nodes and destination ( or target ) is illustrated in fig.[fig : notaion ] where , without loss of generality , all the collaborative sensor nodes are assumed to be located on the @xmath2-@xmath3 plane .
the @xmath10th node location is thus denoted in polar coordinates by @xmath6 .
the location of the destination is given in spherical coordinates by @xmath7 .
following the standard notation in antenna theory @xcite , the angle @xmath11 $ ] denotes the elevation direction , whereas the angle @xmath12 $ ] represents the azimuth direction . in order to simplify the analysis ,
the following assumptions are made :    1 .
the location of each node is chosen randomly , following a uniform distribution within a disk of radius @xmath4 .
each node is equipped with a single ideal isotropic antenna .
all sensor nodes transmit identical energies , and the path losses of all nodes are also identical .
thus the underlying model falls within the framework of phased arrays .
there is no reflection or scattering of the signal .
thus , there is no multipath fading or shadowing .
5 .   the nodes are sufficiently separated that any mutual coupling effects @xcite among the antennas of different sensor nodes are negligible .
furthermore , we also assume that all the nodes are perfectly synchronized so that no frequency offset or phase jitter occurs .
the effects of phase ambiguities on the beampattern will be discussed in section[sec : imperfect ] .
let @xmath13 denote the euclidean distance between the @xmath10th node and the reference location @xmath14 , which is written as @xmath15 if the initial phase of node @xmath16 is set to @xmath17 the corresponding array factor , given the realization of node locations @xmath18 \in [ 0 , r]^n$ ] and @xmath19 \in [ -\pi , \pi]^n$ ] , is written as @xmath20 } , \label{far_field_f0}\end{aligned}\ ] ] where @xmath0 is the number of sensor nodes and @xmath21 is the wavelength of the radio frequency ( rf ) carrier .    in this paper , we are interested in the radiation pattern in the far - field region , and we assume that the far - field condition @xmath22 holds .
the far - field distance @xmath13 in can then be approximated as @xmath23 the far - field beampattern is thus approximated by @xmath24 }   \triangleq \tilde{f}(\phi , \theta | { \boldsymbol r},{\boldsymbol \psi } )   .
\label{far_field_f1}\end{aligned}\ ] ]    alternatively , instead of applying @xmath25 as in , if we choose @xmath26 then we obtain the array factor as @xmath27 }   \nonumber \\ & = e^{j \frac{2 \pi}{\lambda } a } \frac{1}{n }   \sum_{k=1}^{n }   e^{j \frac{2 \pi}{\lambda } r_k \left [ \sin \theta_0 \cos(\phi_0 - \psi_k ) -   \sin \theta \cos(\phi - \psi_k ) \right ] }    \triangleq \tilde{f}^{\dag } ( \phi , \theta    \label{far_field_f2}\end{aligned}\ ] ] the only difference between @xmath28 in and @xmath29 .
the far - field beampattern is thus identical for both systems , and the received signal exhibits no difference as long as the base station compensates for this phase rotation
.    therefore , there are two ways of forming a beam .
one way is to use , but this approach requires accurate knowledge of the distance , relative to the wavelength @xmath21 , between each node and the destination .
thus , this applies to the closed - loop case such as self - phasing arrays .
alternatively , the use of requires knowledge of the node positions relative to some common reference ( such as the origin in this example ) , and thus corresponds to the open - loop case .
knowledge of the elevation direction @xmath8 is also required , but this may be assumed to be known _ a priori _ in many applications . in both cases , the synchronization among sensors is critical , which may be achieved by the use of reference signals such as those of the global positioning system ( gps ) .    of particular interest in practice is the case where @xmath30 , i.e. , the destination node is in the same plane as the collaborative sensor nodes .
therefore , we will consider the beampattern in this plane and thus assume that @xmath31 for the rest of the paper .
we then denote @xmath32 in by @xmath33 and @xmath34 for simplicity .    by assumption , the node locations @xmath6 follow the uniform distribution over the disk of radius @xmath4 .
thus , the probability density functions ( pdfs ) of @xmath35 and @xmath36 are given by @xmath37 from , we have ( with @xmath38 ) @xmath39 where @xmath40 and @xmath41 .
the compound random variable @xmath42 has the following pdf : @xmath43 note that since the above model is symmetric with respect to the azimuth direction @xmath44 , any particular choice of @xmath9 does not change the results in the following .
therefore , without loss of generality , we assume that @xmath45 , and the parameter @xmath44 simply corresponds to the difference angle between the target direction and the reference . also , note that @xmath46 .
the array factor of can then be rewritten as @xmath47 where @xmath48 is the radius of the disk normalized by the wavelength .
finally , the far - field beampattern can be defined as @xmath49 where @xmath50
we start by investigating the average beampattern of the random array resulting from the distributed sensor network model in the previous section . here
, the average is taken over all realizations of @xmath51 , and from the average beampattern is expressed as @xmath52 where @xmath53 denotes expectation with respect to the random variables @xmath54 . from and , it can be readily shown that @xmath55 where @xmath56 is the @xmath57th order bessel function of the first kind . although the function @xmath58 is oscillatory , the local maxima of oscillation tend to decrease with increasing @xmath2 . in , the first term represents the average power level of the sidelobe , which does not depend on the node location , whereas the second term is the contribution of the mainlobe factor . since , conditioned on @xmath44 , the array factor of the form is an average of bounded independent and identically distributed ( i.i.d
) complex random variables , by the weak law of large numbers the beampattern converges to the ensemble average in probability as @xmath59 .
the average beampattern of is plotted in fig.[fig : avrage_beam ] for several values of @xmath60 with @xmath61 and 256 . as can be observed ,
the sidelobe approaches @xmath62 as the beam angle moves away from the target direction .
[ bl][bl][.8]average power [ db ] [ bl][bl][.8]angle [ degree ] [ bl][bl][.8]@xmath63 [ bl][bl][.8]@xmath64 [ bl][bl][.8]@xmath65   and @xmath66.,title="fig:",width=472 ]    to gain further insight , consider the asymptotic expansion of the bessel function @xmath67 for @xmath68 as @xmath69 we then have @xmath70 and becomes , for @xmath71 , @xmath72    the @xmath57th peak of the average sidelobe will appear around @xmath73 , @xmath74 , and its corresponding value becomes @xmath75 ^ 3 , \label{peak_nth}\ ] ] which does not depend on @xmath60 .
the @xmath57th peak and @xmath57th zero positions ( in the sense of the second term in ) can then be expressed asymptotically as @xmath76 since the peak sidelobe value does not depend on @xmath60 and is less sensitive to the value of @xmath0 , it is apparent that the only way one can avoid high peaks in the sidelobe region is to increase @xmath60 such that most of the major peaks are relatively concentrated around the mainlobe .
this phenomenon will be further examined in the following subsections .
one of the important figures of merit in directional antenna design is the 3db beamwidth . in the deterministic antenna , the 3db beamwidth is the threshold angle at which the power of the beampattern drops 3db below that in the target direction @xmath9 . in our scenario ,
the 3db beamwidth itself is a random variable and it is not easy to characterize analytically .
thus , as an alternative measure , we may define the 3db beamwidth of the _ average _ beampattern denoted by @xmath77 as the angle @xmath44 that satisfies @xmath78 in the limit as @xmath79 , one may obtain @xmath80 by numerically solving . for @xmath81 , can be approximated as @xmath82 .
therefore , the beamwidth is asymptotically independent of @xmath0 and is mainly determined by the inverse of the disk radius of the cluster .
consequently , sparsely distributed sensors form a narrow beam on average provided that the cluster radius is sufficiently large .
this sharp mainbeam property may be desirable , but if the far - field destination node has mobility , it should be designed carefully ; the calibration should take place before the mobile node moves out of the mainbeam , but the mainbeam width is inversely proportional to the normalized radius @xmath60 as observed in fig.[fig : avrage_beam ] .
therefore , calibration should be performed more frequently if the destination node moves rapidly or when @xmath60 is increased .
similar to the 3db beamwidth concept , it may also be convenient for our subsequent analysis to define the region within which the average of the sidelobe beampattern falls below some threshold level .
as we have seen , for large @xmath60 , the sidelobe of the average beampattern is given by @xmath62 asymptotically .
therefore , we shall define the _
3db sidelobe region _ as the region in which neither neighboring sidelobe peak in the average beampattern exceeds 3db above @xmath62 . let @xmath83 denote the minimum index of the peak position such that the corresponding peak value satisfies this 3db condition . specifically , from , @xmath83 is the minimum integer @xmath57 that satisfies @xmath84 ^ 3 & \leq 2 , \label{3db_region_def}\end{aligned}\ ] ] and it can be bounded by @xmath85 let @xmath86 denote the angle corresponding to the zero point next to the @xmath83th peak sidelobe which can be obtained from with @xmath87 .
consequently , in this paper , the 3db sidelobe region is defined as @xmath88 fig.[fig:3dbdefs ] illustrates the definitions of the 3db sidelobe region together with that of the 3db beamwidth .
[ r][r][1.0]1 [ r][r][1.0]@xmath89 [ r][r][1.0]@xmath90 [ b][b][.8]1st peak [ b][b][.8]@xmath83th peak [ b][b][.8]1st zero [ b][b][.8]@xmath83th zero [ t][t][.8]3db beamwidth [ b][b][.8]mainbeam [ b][b][.8]3db sidelobe region     as will be shown in section [ sec : zero_mean ] , the idea behind the introduction of 3db sidelobes is that in this region one may assume that the mean value of the random array factor of sampled at @xmath91 becomes a random variable with approximately zero mean .
thus we may reasonably assume that the array factor has zero mean in this region , and this assumption significantly simplifies the analysis of the statistical distribution in the following sections .
fig.[fig : three_side ] shows the threshold angle above which the 3db sidelobe region begins .
the asymptotic 3db beamwidth is also shown for reference . as can be observed , whereas the 3db beamwidth is less sensitive to the number of nodes @xmath0 , the 3db sidelobe region will be considerably reduced as @xmath0 increases .
this means that as @xmath0 increases the dominant non - negligible sidelobe peak may occur with high probability unless @xmath60 is also increased .
this trade - off will be clarified by the study of directivity in the following subsection .
[ bl][bl][.8]@xmath60 [ bl][bl][.8]threshold angle [ degree ]   and number of nodes @xmath0.,title="fig:",width=472 ]      the directivity or directional gain is the parameter that characterizes how much radiated energy is concentrated in the desired direction , relative to a single isotropic antenna .
specifically , it may be defined as @xmath92 where @xmath93 is the radiated power density in the direction of @xmath44 . in the scenario of this paper , since @xmath93 depends on the particular realization of @xmath51 , the corresponding gain may be expressed , by substituting @xmath94 of into the above , as @xmath95^{-1 } .
\label{directivity}\ ] ] the mean value of is given by @xmath96 unfortunately , direct calculation of does not result in a closed - form or insightful expression .
thus , we shall consider the following as an alternative measure : @xmath97 note that by jensen s inequality , we have @xmath98 which means that @xmath99 in is a lower bound on @xmath100 .
however , since by the law of large numbers the denominator of @xmath101 may approach its average value with high probability as @xmath0 increases , the above bound is expected to become tight as @xmath0 increases .
this will be verified in the numerical results below . substituting into
, we obtain @xmath102 where @xmath103 is a generalized hypergeometric function which monotonically decreases with increasing @xmath2 and converges to 0 as @xmath104 .
therefore , unlike well - designed deterministic linear arrays , the gain of a given realization is very likely to be less than @xmath0 , and the limit @xmath0 can be approached only by increasing @xmath60 .
this agrees with the previous observation that the average mainbeam becomes narrow as @xmath60 increases and thus improves the directivity .
it should be noted that although has a simple form and offers some insight into the asymptotic behavior of directivity , the calculation of the generalized hypergeometric function involved in becomes numerically unstable as @xmath60 increases , and it is much easier to numerically integrate the denominator of directly .
[ bl][bl][.8]normalized directivity @xmath105 [ bl][bl][.8]normalized radius @xmath106   and normalized radius @xmath60.,title="fig:",width=472 ]    fig.[fig : directivity ] shows the relationship between the normalized directivity bound @xmath107 and @xmath60 . also shown in the figure as diamond - shaped points are the corresponding exact average directivities @xmath108 obtained by the simulation of 1000 realizations . as can be observed , the bound is very tight compared to the exact performance .
thus , it follows that in order to achieve high normalized directivity ( i.e. directivity close to @xmath0 ) with @xmath0 nodes , the distribution of the nodes should be as sparse as possible .
in fact , we have the following theorem :    [ theorem_directivity ] for large @xmath60 and @xmath0 , @xmath109 is lower bounded by @xmath110 where @xmath111 is a positive constant independent of @xmath0 and @xmath60 ( @xmath112 ) .
see appendix [ app : directive ] .
note that the factor @xmath113 which appears in can be seen as a _ one - dimensional node density_. to verify the above theorem , fig.[fig : directivity_density ] shows the relationship between @xmath107 and the node density @xmath113 , as well as the lower bound in .
[ bl][bl][.8]normalized directivity @xmath114 [ bl][bl][.8]exact ( n=16 ) [ bl][bl][.8]exact ( n=256 ) [ bl][bl][.8]lower bound ( theorem 1 ) [ bl][bl][.8]node density @xmath115   and node density @xmath113.,title="fig:",width=472 ]    the above theorem indicates that there is a simple relationship between directivity and node density .
it can be seen that the node density almost uniquely determines the normalized directivity @xmath109 .
it is important to note that in order to achieve a certain normalized directivity with large numbers of nodes @xmath0 , the node density should be maintained to the desired value by spreading the nodes as sparsely as possible .
alternatively , if the normalized region @xmath60 is fixed , it is not efficient in terms of normalized directivity to increase the number of sensor nodes .
the above theorem also indicates that if the sensor nodes are uniformly distributed and if we are to choose @xmath0 nodes out of them , in terms of normalized directivity it may be better to choose them as sparsely as possible .
in the previous section , we have seen that random arrays have nice average beampatterns with low sidelobes .
however , the average behavior does not necessarily approximate a beampattern of any given realization unless @xmath116 .
in fact , even though the average beampattern has a sharp mainbeam and sidelobes always close to @xmath62 , there is a large dynamic range of sidelobes among randomly generated beampatterns . as an example , the average beampattern and one particular realization of randomly generated sensor locations is shown in fig.[fig : av_vs_instant_beam ] .
the mainbeam of the realization closely matches the average , but sidelobes may fluctuate with a large dynamic range and easily exceed the average level .     and @xmath61.,width=472 ]    therefore , in practice , the statistical _ distribution _ of beampatterns and sidelobes in particular , is of interest . by approximating the beampattern sidelobes as a complex gaussian process ,
lo @xcite has derived the distribution of the beampattern in the case of linear random arrays .    in the following ,
we first derive a numerical method that allows calculation of the exact distribution of the beampattern without applying gaussian approximations .
we then derive a convenient asymptotic form of the sidelobe distribution using a gaussian approximation similar to @xcite , and evaluate its validity in our framework .      since the array factor is a sum of i.i.d .
random variables , its distribution can be computed numerically by the characteristic function method . to this end , from let @xmath117 where @xmath118 and @xmath119 is defined in .
the joint characteristic function of @xmath120 and @xmath121 is written as @xmath122 } \right\ } .
\label{joint_chf1}\end{aligned}\ ] ] for a given pair of @xmath123 and @xmath124 , the above expectation is a single integral of a well - behaved function and can be calculated numerically .    since @xmath125 is a sum of @xmath0 i.i.d .
complex random variables , the joint probability density of @xmath126 and @xmath127 in is given by @xmath128^n   e^{-j(\omega x + \nu y ) } \,d\omega \ , d\nu .\ ] ] the above integral can be computed efficiently using the two - dimensional fast fourier transform ( fft ) .
finally , the complementary cumulative distribution function ( ccdf ) of the beampattern , i.e. , the probability that the instantaneous power of a given realization in the direction @xmath44 exceeds a threshold power , @xmath129 , is given by @xmath130 & = \pr \left [ \frac{\tilde{x}^2 + \tilde{y}^2}{n^2 } > p_0   \right ]   = \int \!\!\ !   \int_{x^2 + y^2 > n^2 p_0 } f_{\tilde{x } , \tilde{y } }   ( x , y ) \,dx\ , dy .
\label{exact_dist_num}\end{aligned}\ ] ]      the exact evaluation of the ccdf outlined above is computationally demanding , especially if the desired numerical precision is high .
considering that the array factor consists of a sum of @xmath0 statistically independent random variables , as @xmath0 increases , by the central limit theorem we may expect that the array factor with any given direction , except at the deterministic angle @xmath131 , approaches a complex gaussian distribution .
this approximation may typically result in a simpler distribution formula . to this end
, we write as @xmath132 where @xmath133 since the @xmath134 s are i.i.d .
random variables , as @xmath0 increases the distribution of @xmath135 and @xmath136 at the direction @xmath137 may approach that of a gaussian random variable with @xmath138 ^ 2 \triangleq \sigma^2_x \label{gauss_var_x }   \\ e \left\ { y \right\ }   & =   0   \label{gauss_mean_y } \\
\text{var } \left ( y \right )   & = \frac{1}{2}\left ( 1 -   \frac{j_1\left ( 2 \alpha ( \phi )   \right ) } { \alpha ( \phi ) }   \right ) \triangleq \sigma^2_y   .
\label{gauss_var_y}\end{aligned}\ ] ] note that @xmath139 , i.e. , @xmath135 and @xmath136 are orthogonal and thus statistically uncorrelated .
the joint pdf of @xmath135 and @xmath136 is then given by @xmath140 the ccdf of @xmath141 can be expressed as @xmath142   & = \pr \left [    \frac{{x}^2 + { y}^2}{n } > { p}_0 \right ]   = \pr \left [    \sqrt { { { x}^2 + { y}^2 } } >
\sqrt{n { p}_0 } \right ]   \nonumber \\ & = \int_{\sqrt{n{p}_0}}^{\infty }   \int_{-\pi}^{\pi } \frac{r}{2 \pi \sigma_x \sigma_y } \exp \left ( -\frac{| r \cos\omega - m_x |^2   } { 2\sigma^2_x } -\frac { r^2 \sin^2 \omega } { 2 \sigma^2_y } \right ) d\omega \ , dr \nonumber \\ & =   \int_{-\pi}^{\pi } \frac{1}{4 \pi \sigma_x \sigma_y u^2_{\omega } } e^{v^2_{\omega } - \frac{m_x^2}{2 \sigma_x^2 } } \left [ \sqrt{\pi } v_{\omega } \text { erfc}\left ( w_{\omega } - v_{\omega } \right ) + e^{-(w_{\omega } -v_{\omega } ) ^2 } \right ] d\omega , \label{exact_ccdf_gauss}\end{aligned}\ ] ] where @xmath143    for @xmath144 , the terms @xmath145 and @xmath146 in the variance expressions and rapidly decrease and their contribution to the resulting variances becomes minor .
therefore , it is very likely that both variances are approximately equal in the sidelobe region . when this is the case , i.e. , if @xmath147 , the distribution of the complex envelope becomes a nakagami - rice distribution .
consequently , the resulting integral can be expressed in terms of the first - order marcum - q function @xmath148    =   q \left ( \frac{m_x}{\sigma_x } , \frac{\sqrt{n { p}_0}}{\sigma_x } \right )    =   q \left ( \sqrt{2 } { m_x } ,   \sqrt{2 n { p}_0 } \right ) .
\label{marcum_formula}\ ] ]    furthermore , if the mean @xmath149 is zero , the envelope follows a rayleigh distribution and we simply have @xmath150   & =   e^{- n { p}_0 } .
\label{full_gaussian}\end{aligned}\ ] ]      as we have seen , if the mean value of the array factor can be assumed to be zero , the distribution can be significantly simplified and thus analysis becomes readily tractable . from it is apparent that under the constant variance constraint the mean value increases as @xmath0 increases .
therefore , when @xmath0 is large , the zero mean assumption may not be guaranteed in general . in section [ sec:3dbsidelobe_region ] , we have introduced the 3db sidelobe region , and in the following we derive properties of the mean value of the array factor in this region .    from the definition of ,
the sidelobe in the 3db region satisfies @xmath151 it follows that @xmath152 from , , and , we have @xmath153 and thus we get @xmath154 therefore , in the 3db sidelobe region , the square of the mean is bounded by unity in the large-@xmath0 asymptote and thus the mean does not grow unbounded with the number of nodes @xmath0 .      , and @xmath155.,width=472 ]    in fig.[fig : distribution_r2 ] , the ccdf s computed with various formulae are shown with @xmath156 and @xmath157 , which corresponds to the sidelobe region .
the exact formula of , the equal variance gaussian approximation of , and the zero - mean gaussian approximation of are shown in the figure .
note that the precise gaussian case of was also calculated by numerical integration , but it is almost identical to the marcum - q function approximation in for this case and thus is not shown . as observed from fig.[fig : distribution_r2 ] , even the zero - mean gaussian approximation may be valid for this sidelobe region , but for @xmath158 the gaussian approximation will have some noticeable discrepancy with the exact value .
this is due to the fact that the zero - mean approximation does not hold for this case .
in fact , fig.[fig : three_side ] indicates that for this value of @xmath0 , the angle falls between the 3db sidelobe region and the mainlobe region and thus the zero mean assumption may not be accurate .
fig.[fig : three_main ] shows the distribution at 3db beamwidth of the average beampattern defined by . in this case , the exact form , the precise gaussian , and the marcum - q approximation show different results even for relatively large @xmath0 , since at this angle the variance of the array factor is small and a large number of random variables must be summed in for a non - zero mean gaussian random variable to adequately approximate the sum . as observed ,
as the number of nodes increases , the mainbeam variance becomes small and approaches the mean value of -3db , as expected .
therefore , for large @xmath0 the mainbeam can be made stable .
this observation agrees with a similar result for linear arrays in @xcite .
[ l][l][.8]exact [ l][l][.8]precise gaussian [ l][l][.8]marcum - q   with exact , precise gaussian , and marcum - q formulas.,title="fig:",width=472 ]
it is well known that unlike periodic or equally - spaced antenna arrays , many arrays with unequal spacing will yield no grating lobes for a large number of elements .
this property is also preserved for random arrays @xcite , but in order to verify this , one may need to find the distribution of the maximum power of the sidelobes . in this section , we develop an approximate upper bound on the distribution of the peak sidelobes for random sensor networks .    in the previous section
, we have seen that the distribution of the beampattern samples within the sidelobe region can be characterized by a zero mean gaussian random variable if the zero - mean condition is satisfied . in the following ,
we further assume that the beampattern is a gaussian random _
process_. in this case , any two samples taken from the beampattern should be characterized by jointly gaussian random variables . in the linear random array framework ,
the distribution of peak sidelobes has been studied in @xcite , assuming the array factor is a gaussian process . for simplicity ,
only the 3db sidelobe region is considered and it will be assumed that the process is stationary with zero mean .
the extension to the non - stationary case is studied in @xcite .    in the following ,
the ccdf of the maximum peak sidelobe , which is the probability that the maximum peak sidelobe exceeds a given power level , will also be referred to as _ outage probability _ and denoted by @xmath159 .
let @xmath160 denote the random variable representing the number of upward crossings at a given level @xmath161 per interval in the 3db sidelobe region @xmath162 .
as shown in appendix [ app : level_crossing ] , assuming that the array factor in this region can be approximated as a zero - mean gaussian process , the mean of @xmath160 is given by @xmath163 note that the above function monotonically decreases with increasing @xmath161 only for @xmath164 , and thus is meaningful only in this region .
finally , noticing that the outage probability is the probability that at least one peak exceeds level @xmath161 and is given by @xcite @xmath165 & = \sum_{k=1}^{\infty } \pr [ \nu(a ) = k ]   \leq \sum_{k=1}^{\infty } k \pr [ \nu(a ) = k ] = e \left\ {   \nu(a ) \right\},\end{aligned}\ ] ] then serves as an upper bound for the outage probability for the maximum sidelobe peak for @xmath164 .
thus , we obtain the ccdf upper bound as @xmath166 & \leq 4 \left ( 1 -   \sin \frac{\phi^{\text{zero}}_{n_0}}{2 } \right ) \sqrt{\pi } \tilde{r } \sqrt{n p_0 } e^{-n p_0 } , \qquad   \text { for } \quad n p_0   > \frac{1}{2 } .
\label{upperbound_sidelobe}\end{aligned}\ ] ]      fig.[fig : ccdf_sidelobepeak ] shows a comparison between simulation results and the upper bound . for the simulation , the outage probability
is calculated based on 10000 randomly generated realizations with the node density @xmath167 and only the peaks within the 3db sidelobe region are examined .
also , in order to capture peak values accurately , the entire 3db sidelobe region of @xmath44 is sampled at a rate as large as @xmath168 . as can be observed ,
the bound is in good agreement with simulation for large @xmath0 .
[ bl][bl][.8]@xmath167 .,title="fig:",width=472 ]    let @xmath169 denote the threshold of the maximum peak value ( @xmath129 ) normalized by the average sidelobe level ( @xmath62 ) .
since from @xmath170 approaches zero as @xmath60 increases , reduces to @xmath171 the above inequality illuminates the relationship between the outage probability and @xmath60 ( assuming that @xmath170 is negligibly small ) .
fig.[fig : outage_sidelobepeak ] shows the maximum possible value of @xmath172 for a given outage probability and @xmath60 .
as can be observed , the maximum sidelobe may grow as @xmath60 increases , but the amount is below 12db for many cases of interest .
consequently , the maximum sidelobe level ( in the 3db region ) of randomly generated arrays may be written as @xmath173 , where the required margin @xmath174 depends on the parameters @xmath60 and @xmath175 , but not on @xmath0 .
thus , increasing @xmath0 always results in a reduction of maximum sidelobe level .
[ bl][bl][.8]@xmath176 [ bl][bl][.8]@xmath177 [ bl][bl][.8]@xmath178 [ bl][bl][.8]normalized radius@xmath60 [ b][b][.8]normalized sidelobe level @xmath172 .,title="fig:",width=472 ]
so far , we have evaluated the beampattern assuming perfect knowledge of the initial phase for each node . in this section ,
we analyze the effect of the phase ambiguities in the closed - loop scenario as well as location estimation errors in the open - loop scenario . for each of the two scenarios ,
we derive the average beampattern and calculate the amount of mainbeam degradation .      in the closed - loop case ,
the effects of imperfect phase may be easily derived , following the approach developed by steinberg @xcite .
the initial phase of node @xmath10 in will now be given by @xmath179 where @xmath180 corresponds to the phase offset due to the phase ambiguity caused by carrier phase jitter or offset between the transmitter and receiver nodes . in the following ,
the phase offset @xmath180 s are assumed to be i.i.d . random variables .
then , from , , , and , the far - field array factor ( with @xmath181 ) will be given by @xmath182 the average beampattern of will be replaced by @xmath183 similar to , direct calculation of results in @xmath184 where @xmath185 thus , as @xmath79 , the average beampattern will simply become a version of the original scaled by a factor of @xmath186 .
let us now assume that the phase offset follows a tikhonov distribution , a typical phase jitter model for phase - locked loop ( pll ) circuits given in @xcite , @xmath187 where @xmath188 is the variance of the phase noise and @xmath189 is the @xmath57th order modified bessel function of the first kind .
the corresponding attenuation factor is given by @xmath190 the variance of the phase noise @xmath188 is related to the loop snr of the pll by @xmath191 fig .
[ fig : degradation_close ] shows the degradation factor @xmath192 with respect to the loop snr .
as observed from the figure , a loop snr of at least 3db may be necessary for each node in order to reduce the overall beampattern degradation to less than 3db .          in the open - loop case , our model of the initial phase
is given in with @xmath30 , and if there are estimation errors in the location parameters @xmath35 and @xmath36 , the initial phase will be replaced by @xmath193 where @xmath194 and @xmath195 are the corresponding error random variables , each set assumed to be i.i.d . and also independent of @xmath35 and @xmath36 for simplicity .
with the far - field approximation , we have @xmath196   + \delta { r}_k \cos(\phi_0 - ( { \psi}_k   + \delta \psi_k ) ) \right\ } \nonumber \\ & =   \frac{2\pi}{\lambda } a +   \frac{4 \pi}{\lambda } { r}_k \left [ \sin\left ( \psi_k -   \frac{\phi_0 + \phi - \delta \psi_k}{2 } \right ) \sin\left ( \frac{\phi_0 - \phi - \delta \psi_k}{2 } \right ) \right ] \nonumber \\ & \quad + \frac{2\pi}{\lambda } \delta { r}_k   \cos ( { \psi}_k    -   ( \phi_0 - \delta { \psi}_k ) )   .
\label{dist_phase2}\end{aligned}\ ] ] let @xmath197 .
then , the right - hand side ( rhs ) of is given by @xmath198 the resulting far - field array factor of will then be given by @xmath199 and the beampattern is expressed as @xmath200 where @xmath201 conditioned on @xmath202 and @xmath195 , the angle @xmath203 can be seen as a uniformly distributed random variable , and thus the pdf of @xmath134 is given by .
considering the fact that @xmath35 and @xmath194 are assumed to be statistically independent , we further assume for analytical purposes that @xmath134 and @xmath204 are statistically independent .
then , again , the beampattern does not depend on the particular choice of @xmath9 .
furthermore , on modeling @xmath205 as being uniformly distributed over @xmath206 $ ] and assuming the phase term of @xmath204 to be uniformly distributed over @xmath207 $ ] , the probability density function of @xmath204 will be given by @xmath208 , \qquad |v| \leq r_{\text{max } } .\ ] ] consequently , the average beampattern can be written as @xmath209 where @xmath210 and without loss of generality @xmath211 was assumed . in , @xmath212 denotes a generalized hypergeometric function which has an oscillatory tail but converges to zero as @xmath2 increases .
also , assuming that the @xmath195 are uniformly distributed over @xmath213 $ ] and using the approximation @xmath214 which is valid for the beampattern around the mainbeam , we obtain @xmath215 since the hypergeometric function @xmath216 has a maximum peak value of 1 at @xmath217 , the above expression indicates that regardless of the value of @xmath60 , there may be two symmetric peaks around the mainbeam at @xmath218 resulting in a _ pointing error_. therefore , the mainbeam may spread over by a factor of @xmath219 . at the center of the mainbeam , we have @xmath220    fig . [
fig : degradation_open ] shows the degradation factor @xmath221 and @xmath222 for a given @xmath223 and @xmath224 .
as observed from the figure and discussion above , the angle estimation error has two effects , i.e. , pointing error and mainbeam degradation .
in particular , if we wish to suppress the mainbeam degradation below 3db , from the figure , we should choose @xmath225 .
this means that the maximum angle estimation error should satisfy @xmath226 and as @xmath60 becomes large , the requirement of minimum angle ambiguity from becomes severe .
[ bl][bl][.8]@xmath227 @xmath228
in this paper , we have analyzed the stochastic performance of random arrays for distributed collaborative beamforming , in the framework of wireless ad hoc sensor networks .
it has been shown that under ideal channel and system assumptions , directivity of order @xmath0 can be achieved asymptotically with @xmath0 sensor nodes , as long as the sensor nodes are located sparsely enough .
we have studied the average and the distribution of the beampattern as well as the distribution of the sidelobe peaks .
several forms of the ccdf of the beampattern have been derived and compared , with particular emphasis on the gaussian approximation of the array factor .
we have considered two scenarios of distributed beamforming and investigated the effects of phase ambiguity and location estimation error upon the resultant average beampatterns .
our main conclusion is that , given a number of nodes randomly distributed over a large disk , one may form a nice beampattern with narrow mainlobe and sidelobes as low as @xmath62 plus some margin for maximum sidelobe peaks . also , the directivity approaches @xmath0 if the nodes are located as sparsely as possible
however , our analysis is based on a number of ideal assumptions on the system and channel model . in practice ,
a number of open issues remain , such as applicability of beamforming when the destination or nodes in the cluster are in rapid motion or the channel suffers severe multipath fading .
also , specific algorithms should be developed for frequency offset correction of each node as well as methods for initial phase or location estimation . finally , efficient protocols for sharing the transmit as well as calibration information among nodes are required .
we first prove the following lemma :    [ lemma_hyper ] a generalized hypergeometric function @xmath103 with @xmath68 can be bounded as @xmath229 where @xmath230 is a constant ( @xmath231 ) .
we start with the integral form @xmath232 since the asymptotic form of @xmath233 given by is valid for @xmath68 , we have the following inequalities @xmath234 @xmath235 for some threshold value @xmath236 which should be determined numerically , and for some interval @xmath237 .
the parameter @xmath238 is chosen to be the smallest non - negative value such that @xmath239 should hold , and this guarantees a continuity of the function at the threshold @xmath240 .
fig.[fig : function_approximation ] illustrates the relationship of and with @xmath236 chosen as a cross point of @xmath241 between the functions @xmath242 and @xmath243 , yielding @xmath244 .
the corresponding value of @xmath238 is 0.4664 .
[ bl][bl][1.0]function of @xmath241 [ bl][bl][1.0]@xmath241 [ bl][bl][.8]@xmath245 [ bl][bl][.8]@xmath246 [ bl][bl][.8]@xmath247   and its upper bound with @xmath248.,title="fig:",width=472 ]    substituting and into , we get for @xmath249 @xmath250 the first term on the rhs of is given by @xmath251 the second term on the rhs of is given by @xmath252 \nonumber \\ & = \frac{8}{\pi^2 x^2_0 } \frac{1}{x } \left ( 1 - \left ( \frac{x_0}{x } \right)^2 + o\left(1/x^4 \right ) \right ) + o\left ( \ln(x)/x^3 \right ) \nonumber \\ & = \frac{1}{x } \frac{8}{\pi^2 x^2_0 } + o\left ( \ln(x)/x^3 \right ) .
\label{fx3_2}\end{aligned}\ ] ] consequently , we may write @xmath253 and the second term on the rhs of drops as @xmath2 becomes large . with @xmath244 and @xmath254 ,
the coefficient of @xmath255 can be calculated to be @xmath256 .    from , , and lemma [ lemma_hyper ] , we have @xmath257 for large @xmath0 , the rhs of converges to with @xmath258 .
in this appendix , we obtain the mean number of upward crossings of a given level of the zero mean gaussian process based on the approach of rice @xcite .
assume that @xmath135 and @xmath136 are uncorrelated zero - mean gaussian processes with variance @xmath259 .
let @xmath260 and @xmath261 and @xmath262 denote the corresponding processes differentiated by @xmath263 . by assumption , @xmath261 and
@xmath262 become zero mean gaussian processes . in order to calculate the variance , first consider the autocorrelation function of @xmath135 at instants @xmath264 and @xmath265 given by @xmath266 where the other terms become zero by the zero mean assumption .
also for the same reason , the first term of the rhs of may be approximated by zero .
therefore , letting @xmath267 , we obtain @xmath268 differentiating the above with respect to @xmath269 twice , setting @xmath270 , and carrying out the statistical average with respect to @xmath1 , the variance of @xmath261 is given by @xcite @xmath271 likewise , one may obtain @xmath272 , and the joint pdf of @xmath273 is given by @xmath274 on changing the random variables in the polar coordinates via @xmath275 , @xmath276 and integrating out @xmath277 and @xmath278 , we obtain @xmath279 the number of positive ( upward ) crossings of the process @xmath123 at level @xmath161 per interval @xmath280 is given by @xcite @xmath281 consequently , the mean number of upward crossings for the interval @xmath162 is given by @xmath282 which results in .
j.  c. chen , l.  yip , j.  elson , h.  wang , d.  maniezzo , r.  e. hudson , k.  yao , and d.  estrin , `` coherent acoustic array processing and localization on wireless sensor networks , '' _ proc .
ieee _ , vol .
91 , pp .  11541162 , aug .