a century ago , the study of variable stars was a leading area of astronomy . with the increasing availability of large format optical detectors and inexpensive high speed computing this sub - field ,
as witnessed by projects such as asas , ogle , catalina sky survey , the palomar transient factory , panstarss and skymapper , is making a come - back .
radio astronomy appears to be poised for a similar growth . at meter wavelengths ,
commercially available signal processing chips make it feasible to image the entire primary beam of a dipole or a cluster of dipoles .
these technological innovations lie at the heart of lofar @xcite , mwa @xcite and lwa @xcite . at centimeter wavelengths the `` large number , small diameter '' ( lnsd ) array approach ( made possible by inexpensive signal processing , advances in commercial rf technology , innovative ideas in the design of small diameter telescopes and phased array focal planes ) has now been demonstrated to be a cost effective method of building high speed mapping machines @xcite .
the lnsd approach has motivated a new generation of radio facilities : apertif / wsrt @xcite ; meerkat @xcite ; askap @xcite .    lllllll[hbtp ] type ii sne & 0.1 - 1 & 10 & 100mpc & 16 & 0.04 & @xcite + type ib / c sne & 0.1 & 0.3 & 50mpc & 14.5 & @xmath1 & @xcite + sn1998bw - like & 0.1 & 0.1 & 300mpc & 18.4 & @xmath2 & @xcite + swj1644 + 57-like & 0.1 & 1 & @xmath3 & 21.7 & 0.1 & @xcite + orphan afterglows & 1 & 1 & 1gpc & 21.0 & @xmath4 & @xcite + ns - ns mergers & 0.11 & 0.13 & 800mpc & 20.5 & @xmath5 & @xcite + [ tab : listoftrans ]    we divide radio transients into four categories based on two attributes .
the first is the duration of the basic phenomenon ( shorter than or greater than a few seconds ) .
the second is their location ( within the galaxy or extra - galactic ) . roughly speaking the duration maps to coherent versus incoherent emission and the location to repeated versus cataclysmic events .
pulsars and related phenomenon ( giant pulses , nulling pulsars , erratic pulsars , rotating radio transients , and magnetars ) are the dominant category of short duration radio transients at meter and centimeter wavelengths .
there are no secure examples of short duration radio transients that are located beyond the local group . flare stars and associated phenomena are prime examples of long duration radio transients of galactic origin .
the focus of this paper is long duration transients of extra - galactic origin .
known examples in this group are supernovae @xcite and gamma - ray burst afterglows @xcite . in both cases , the radio emission arises as the fast moving debris interacts with the circumstellar matter . in table
[ tab : listoftrans ] we summarize the areal density of radio - emitting supernovae ( including the sub - classes ) and grb afterglows .
note the areal density of `` live transients '' ( transients present at any given instant of time ) of both supernovae and grb afterglows is less than 0.05 per square degree .    in 2007 ,
bower et al .
[ hereafter , b07 ] reported on the analysis of a single field observed every week as a part of the very large array ( vla ) calibration program . the observations were conducted at 4.8ghz and 8.4ghz and lasted 22 years .
the 944 epochs and the weekly cadence makes this data set a most valuable set to probe the decimeter band for long duration transients at the sub - millijansky level . these authors reported the discovery of eight transients found in only one epoch ( hereafter `` single - epoch '' ; duration , 20minutes @xmath6 1week ) and two transients found in rolling 2-month searches ( hereafter , `` multi - epoch '' sources ) .
deep observations towards these sources were undertaken at optical , near - ir and x - ray bands .
the most remarkable feature of the b07 sources is an absence of optical and near - ir counterparts , despite deep searches ( b07 ; @xcite ) . as noted by @xcite all extra - galactic transients (
regardless of the band at which the transient was discovered ) have detectable optical counterparts , namely , their host galaxies .
furthermore , remarkably the areal density of transients live at any given time was estimated to be 1.5deg@xmath7 ( @xmath8mjy ) .
this density exceeds that of all other known radio transient source populations by an order of magnitude ( or more ) ; see table  [ tab : listoftrans ] .
@xcite thus argued that the absence of an optical counterpart means that b07 transients have to be repeating sources of galactic origin , and proposed that b07 transients are old neutron stars ( which naturally satisfy the requirement of being optically almost invisible ) .
given that the search for transient and strong variables is one of the primary motivators for the next generation radio facilities ( described earlier ) it is important to critically investigate the b07 transients since this class nominally dominates over all other known classes of radio transients ( see table  [ tab : listoftrans ] ) .
to this end , here we report on a re - analysis of the original data of b07 (  [ sec : obsreana ] ,  [ sec : findings ] ) and revisit the transient reported by @xcite (  [ sec : othersurveys ] ) . in  [ sec : lognlogs ] we present an update of the expected rate of radio transients . in  [ sec : futuresurveys ] we discuss these rates in relation to future synoptic radio imaging surveys and conclude .
the data used by b07 arose from a calibrator program carried out during the period 19832005 .
all observations were made in the standard continuum mode with 100mhz of total bandwidth in each of two adjacent 50-mhz bands ( ifs ) at center frequencies of 5ghz and 8.4ghz and in both hands of circular polarization .
see b07 for more details about the full data - set .    for the re - analysis we selected , from the archive , only the raw data relevant to the transients reported in b07 .
this means the eight epochs from which the single - epoch transients were first found and the 3 + 8 data sets from which the two multi - epoch transients were found .
data were taken at other radio frequencies in about half of the cases .
some details of the single epoch and two - month transients can be found in table  [ tab : listoffields ] ) .    for the re - analysis we used aips @xcite .
the data reduction and imaging followed the same path used by b07 with a small exception . b07 employed aips for the flagging and analysis of the single - epoch transients , and used the _ miriad _ package @xcite for imaging the two - month averages .
we endeavored to make the calibration and the flagging of uv data ( aips task ` tvflg ` ) data for each epoch as uniform way . following these steps we ran each raw visibility data set through the vla pipeline ( ` vlarun ` ) .
no flux density calibrator was observed during any epoch of these test observations . following b07 ,
the flux density of the phase calibrator ( b1803 + 784 ) was fixed to be 2.2jy ( 5ghz ) and 2.8jy ( 8.4ghz ) .
for those epochs with 22ghz and 1.4ghz observations the flux density of the phase calibrator was taken to be 3jy and 2jy , respectively .
it is evident from the strong variations in the radio light curves for b1803 + 784 that these mean values are only approximate .
our reinvestigation confirm that at least during the period 19811999 the variation was less than 15% .
fortunately , an accurate flux density scale is not crucial for our analysis since we report results in terms of the signal - to - noise .
following b07 we applied a gaussian weighting to the visibility data in order to limit the effects of bandwidth smearing .
this was done by applying a 150-k@xmath9 taper to all visibility data prior to imaging ( ` imagr ` ) . for each field
we required that a source be present in both frequency bands ( ifs ) with similar flux densities and with similar positions .
images were made with extra large fields - of - view .
the wide field - of - view is necessary to reduce the effect of side - lobes that can mimic sources in narrow fields .
these final analysis images had a size of about 40-arcmin at 5ghz data and 27-arcmin at 8.4ghz .
for guidance , the full - width at half power for vla antennas is 45-arcmin/@xmath10(ghz ) , or 9.3-arcmin at 5ghz and 4.3-arcmin for 8.4ghz .
measurements of the vla beam power response beyond the first null are given in @xcite , while the polynomial coefficients needed to correct for the primary beam attenuation can be found in the aips task ` pbcor ` .
we do note that these corrections are uncertain for large angular distances from the beam center .
all data were taken in the b1950 coordinate system .
we stayed in the b1950 system throughout calibration and imaging .
lllllll rt19840502 & se & 4.9 & 4.6 & 0.22 & 6.0 & phase center artifact .
+ rt19840613 & se & 4.9 & 4.6 & 3.3 & 5.7 & side - lobe of bright source . + rt19860115 & se & 4.9 & 4.6 & 1.3 & 14.8 & side - lobe of bright source .
+ rt19860122 & se & 4.9 & 4.6 & 4.6 & 14.5 & artifact .
lower if is bad .
+ rt19870422 & 2 m & 4.9 & 4.6 & 8.0 & 12.8 & artifact .
bad pointing .
+ rt20010331 & 2 m & 8.5 & 2.6 & 4.4 & 1.5 & no detection ( see  [ sec : rt20010331 ] ) .
+   + rt19920826 & se & 4.9 & 4.6 & 2.0 & 21.2 & * d *  snr : _ 6.4 _ ; ` 5.8 ` .
+ rt19970205 & se & 8.5 & 2.6 & 4.4 & 1.4 & * b *   snr : _ 7.7 _ ; ` 5.7 ` ( ca ) + rt19970528 & se & 4.9 & 4.6 & 6.8 & 3.9 & * cnb *  snr : _ 7.5 _ ; ` 5.6 ` .
( ca ) + rt19990504 & se & 4.9 & 4.6 & 8.9 & 18.9 & * d *  snr : _ 7.3 _ ; ` 5.7 ` .
+   + rt19970528 & se & 8.5 & 2.6 & 6.8 & 1.3 & * cnb *   no detection ( ba ) + rt19990504 & se & 8.5 & 2.6 & 8.9 & 8.3 & * d *   no detection ( ba )    [ tab : listoffields ]
below we offer a detailed report for each of the ten sources reported in b07 . summarizing our results for the impatient reader we find four of the eight single - epoch transients [ rt19840502 (  [ sec : rt19840502 ] ) , rt19860122 (  [ sec : rt19860122 ] ) , rt19840613 (  [ sec : rt19840613 ] ) , rt19860115 (  [ sec : rt19860115 ] ) ] and one of the two - month transients ( rt19870422 ;  [ sec : rt19870422 ] ) to be artifacts .
our re - analysis finds that the remaining long - duration transient , rt20010331 (  [ sec : rt20010331 ] ) did not have a significant detection when imaged with aips ( snr of 3.2 ) but the _ miriad _ imaging yields an snr of 7.4 .    for
the remaining single epoch sources [ rt19920826 , rt19970205 , rt19970528 , rt19990504 ] our analysis ( undertaken by daf ; see  [ sec : obsreana ] for a summary of our data reduction ) finds reduced snr .
the discrepancy with respect to the results reported in b07 concerned us and so another author ( gb ) reanalyzed these fields .
some of the discrepancies arise from uv data flagging but even when the same flagging is used different snr algorithms yield measures which differ by 1 to 1.5@xmath11 .
since the sources lie close to threshold of detection even a small shift of a single @xmath11 has an exponential effect in their confidence . in effect
, the reality of this group ( taken as a whole ) apparently depends on details of algorithms in aips and _ miriad _  the investigation of which is beyond the scope of this paper .
[ sec : rt19920826][sec : rt19990504 ] we detail the analyses and a summary of the snrs of these sources can be found in table  [ tab : listoffields ] .
we also re - investigated the lone transient candidate reported in @xcite .
our re - examination which now correctly includes the number of independent beams that were searched for shows that there is a four percent probability that this candidate is due to noise ( see  [ sec : othersurveys ] ) . in the next section (  [ sec : revisedb07 ] )
we synthesize these findings and present our conclusions about the b07 transients .
b07 report finding a transient close to the pointing center ( 13@xmath12 ) with a primary beam - corrected flux density of @xmath13jy , or a snr=6.1 . we imaged the calibrated data set and confirmed the presence of emission at this level at the reported position .
however , we identified a phase center artifact in the visibility data for the left hand polarization of the upper if band .
the effect of this artifact in the image plane is to create strong positive and negative slidelobes , with the positive feature identified as a transient source ( figure  [ fig:19840502 ] ) . when the upper if band data are removed the resulting image is noise - like and the peak flux density at the nominal position of the reported transient is @xmath14jy .
additional observations were taken during the same epoch at 15ghz .
the peak flux density at the same position is @xmath15jy .
, @xmath16 , @xmath170.5 , 0.5 , 1.0 , 1.5mjy , with negative ( positive ) contours given by dashed ( sold ) lines .
the synthesized beam size is shown in the lower left . ]
b07 report finding a transient coincident with a possible host galaxy ( @xmath18 ) and with a primary beam - corrected flux density of @xmath19jy , or a snr=7.0 .
our deconvolved image shows a source at that location .
gaussian fitting suggests that the source is resolved and this conclusion is supported with an integrated flux density ( @xmath20jy ) being clearly larger than the estimate of the peak flux density of @xmath21jy .
imaging the lower and upper ifs separately we find another discrepancy . the peak flux density in the lower if is four times weaker than the upper if band .
the likely source of the problem can be seen by comparing the dirty image with the dirty beam ( figure [ fig:19840613 ] ) .
rt19840613 appears to be an uncleaned side - lobe of j150123 + 781806 , one of the brighter persistent sources detected in 452 images made by b07 at 5ghz .
the putative transient is 56@xmath12 away from j150123 + 781806 to the northwest , close to a local maximum ( 10% of peak ) in the dirty beam at this location .
this side - lobe artifact is stronger in the upper band but it is still present in the lower band .
deconvolution does not fully remove the side - lobe from the image and the effect is to produce a false transient .
we investigated whether the artifact is due to short - timescale ( @xmath22min ) variability of j150123 + 781806 but after dividing the visibility data in half ( by time ) and re - imaging , we found no evidence for variability .
the artifact may be due to some low level interference picked up by some antennas baselines but we were not able to identify the bad data .
it is possible that rt19840613 is a real transient that unfortunately lies on the side - lobe of the dirty beam , but its gaussian fit parameters and its variation in the lower and upper sidebands suggest that it is not a real source .
additional observations were made during this epoch at 1.5ghz .
the data quality is good and the images have no obvious artifacts .
no source is visible at the transient position .
the peak flux density is @xmath23jy .    .
note that the transient candidate rt19840613 lies at the same angle and position as a side - lobe from the bright source j150123 + 781806 located in the bottom right corner of this image .
, title="fig : " ] + .
note that the transient candidate rt19840613 lies at the same angle and position as a side - lobe from the bright source j150123 + 781806 located in the bottom right corner of this image .
, title="fig : " ]      b07 report a transient with a primary beam - corrected flux density of @xmath24jy , or a snr=5.5 .
some of the same issues with the image of the previously discussed rt19840613 were also seen for rt19860115 . in figure  [ fig:19860115 ]
we show the dirty image along with the synthesized beam .
the rt19860115 appears to be an uncleaned side - lobe of j150123 + 781806 and hence not a real transient .
rt19860115 lies at the same angular distance ( 3.3@xmath25 ) and position angle ( 100@xmath26 ccw cc ) of a local maximum in the side - lobe structure ( 25% of the peak beam ) .
we were unable to identify the source of these strong side - lobes .
as in the case of rt19840613 , we were able to rule out that the strong side - lobes originated from short - term variability of j150123 + 781806 .    .
note that the transient candidate rt19860115 lies at the same angle and position as a side - lobe from the bright source j150123 + 781806 located in the upper right corner of this image .
, title="fig : " ] + .
note that the transient candidate rt19860115 lies at the same angle and position as a side - lobe from the bright source j150123 + 781806 located in the upper right corner of this image .
, title="fig : " ]      b07 report a transient with a primary beam - corrected flux density of @xmath27jy , or a snr=6.4 .
images made with all the visibility data do show emission at the position of rt19860122 with the reported peak flux density .
however , there does appear to be some erroneous visibility data resulting in low - level rings in the image centered at the phase center .
the bad data were traced to the lower if . if the data for the lower if band are removed there is no source at the reported transient position ( peak flux density of @xmath28jy ) .
this is one of two transients found by binning individual epochs into two - month averages .
b07 report a 5ghz transient with a primary beam - corrected flux density of @xmath29jy , or a snr of 6.1 .
the source j150050 + 780945.5 is positionally coincident with a blue host galaxy ( @xmath30 ) and has been suggested by @xcite as a candidate afterglow from a binary coalescence event .
three single epochs were used to form the average .
they are in yyyymmdd format : 19870414 , 19870422 and 19870429 .
a close inspection of the data reveals a problem .
the last two of the three epochs used to find rt19870422 were pointed to a different part of the sky , an error that was made when the original observing files were written .
the source is actually j210133.67 + 373528 , and appears only as a transient when added inappropriately to the data from the standard on - the - sky test field .
in hindsight we recognize these to be the same type of errors as seen earlier ( see * ? ? ?
this is the second of the transients which were found by binning individual epochs into two - month averages .
b07 report a 5ghz transient with a primary beam - corrected flux density of @xmath31jy , or a snr=7.4 .
separately , @xcite reported a marginal x - ray source at this position .
eight single epochs were used to form the average .
they are in yyyymmdd format : 20010306 , 20010314 , 20010321 , 20010328 , 20010403 , 20010411 , 20010418 and 20010425 .
we calibrated and imaged the eight epochs of observations from 2001 march 6 to 2001 april 25 ( all in b - configuration ) .
while the four quiescent sources from b07 can be seen in this deep image , we do not identify a significant source at the position of rt20010331 .
the brightest peak within the synthesized beam ( natural weighting ) is @xmath32jy ( or 3.2@xmath11 ) .
an independent reduction using the aips package by gb confirms the absence of this source .
the first problem lies with the primary beam corrections reported in b07 .
rt20010331 lies 4.4@xmath33 from the phase center , close to the 10% response radius of the primary beam .
our estimate for the rms noise of @xmath34jy lies within a few percent of the theoretical value .
the beam - corrected rms noise in this case ( at 10% response of the primary beam ) would then be about @xmath35jy , not the @xmath36jy given in b07 . rather than correcting at the 10% radius
, it appears that the flux density and noise in b07 were mistakenly corrected at the 14% power level .
this multiplicative error has no impact on the signal - to - noise .    in order to investigate this signal - to - noise discrepancy between our image and b07
, we split the data in various ways ( separate epochs , months of march and april epochs , adjacent if bands ) and re - imaged , looking for a bright peak .
none were found .
the discrepancy between b07 and the work reported apparently can be traced back to differences in the _ miriad _ and aips imaging packages . our calibrated visibility data , when processed through _ miriad _ using nearly identical imaging parameters as those in aips , gives a flux density of @xmath37jy ( 7@xmath11 ) .
we have no explanation for the discrepancy between the two results obtained from aips and _ miriad_. it is worth noting that the peak flux densities of the persistent sources identified by b07 agree in these images to @xmath38 . for this paper
we accept the analysis given here .
there was one epoch ( 2001 march 6 ) in which data were also taken at 5ghz . the peak flux density at the position of rt20010331 is @xmath39jy . in summary
, we find no evidence to support that rt20010331 is a significant detection .
b07 report a transient with a primary beam - corrected flux density of @xmath40jy , or a snr=6.4 .
we confirm a source at this position but with a slightly reduced snr .
using natural weighting of the gridded visibilities , our measured peak flux density is @xmath41jy , or 5.8@xmath11 . despite the lower significance
, there is some confidence that rt19920826 is real since it appears in both ifs with comparable flux densities .
further investigation by gb shows that the difference in the snr between b07 and the analysis here can be traced to differences in flagging of the uv data ( of two specific antennas ) .
b07 report a transient with a primary beam - corrected flux density of @xmath42jy , or a snr=7.5 .
these observations were taken during a time when some antennas were being moved to the b configuration and so we applied antenna position corrections ( via aips task ` vlant ` ) before calibration .
rt19970528 is 6.8@xmath25 from the center of the image , where the response of the antennas is only 16% of their peak .
the uncorrected flux density is @xmath43 .
a point source search ( aips task ` sad ` ) of the four million pixels enclosed interior to a radius around this candidate , shows six other un - cataloged candidates with similar signal - to - noise .
the primary - beam corrected flux density @xmath44mjy is similar that from to b07 .
this estimate does not include a correction for temporal smearing due to the sources distance from the phase center , nor for the added uncertainty in the magnitude of the primary beam correction beyond the 20% point . correcting for these terms we get @xmath45mjy .
gb re - investigated this source and found that with the b07 analysis the snr varies with the approach ( peak / rms , ` sad ` , ` jmfit ` ) between 6.2 , 7.4 and 7.1 .
( without any corrections discussed above ) .
the corresponding snrs for the image discussed above is 6.1 , 5.1 and 6.8 .
the source is a bit `` ratty '' and this may explain the variation in snr .
we therefore find this source to be a weak detection .
additional observations were made during this epoch at 8.5ghz .
the source lies close to the first null ( @xmath46 from the phase center ) of the beam at this frequency .
the attenuation by the primary beam is severe and hence the sensitivity is not sufficient to provide any strong spectral index constraints .
b07 report a transient with a primary beam - corrected flux density of @xmath47jy , or a snr=7.3 .
@xcite report an x - ray source in the vicinity of the radio source .
deep multi - wavelength data is consistent with the x - ray source arising from a qso but located five arcseconds from the putative radio transient .
we find the ( uncorrected ) peak flux density is @xmath48jy ( or an snr of 5.7 ) .
rt19990504 lies 8.9@xmath33 from the phase center , close to the first null ( @xmath49 from the phase center ) where we would not expect to find sources .
the polynomial expressions used in aips to correct for the beam attenuation are increasingly inaccurate outside the 20% response radius , and they are not applicable close to the null .
the azimuthally averaged measured value of the beam response is 1.8% , implying a flux density correction factor of 55@xmath50 @xcite or @xmath51mjy .
as with rt19970528 , the analysis of b07 data by gb finds snr of 5.6 ( peak / rms ) , 7.3 ( ` sad ` ) and 6.8 ( ` jmfit ` ) and a similar variation with the image reported here .
it is worth noting that the source may be extended and also that a visual inspection of the annular region reveal a number of similar sources .
we conclude that rt19990504 is not a robust detection .
additional observations were made during this epoch at 8.5ghz .
an image was made from these higher frequency data but , as expected , no source was found .
this is the only single - epoch transient identified at 8.5ghz by b07 .
they report a primary beam - corrected flux density of @xmath52jy , or a snr=7.8 .
rt19970205 is 4.4@xmath25 from the center of the image , where the response of the antennas is only 9.8% of their peak .
there is an elongated source at this position with a peak flux density of @xmath53jy ( using natural weighting ) .
there is some indication that the source is extended since the integrated flux density is about twice the peak flux density .
the peak flux density was likely underestimated by b07 since they did not fully correct for bandwidth or temporal smearing effects .
these data were taken in the bna array with an integration time of 3 1/3 s and a 50-mhz bandwidth . in this observing configuration , the effects of temporal smearing with this dump time will likely reduce the peak flux density by only a few percent .
however , chromatic aberration is expected to be larger with the effect of smearing a point source along the radial direction @xcite .
accounting for the bandwidth effect , we measure a peak flux density of @xmath54mjy ( snr=5.7 ) .
next , due to the uncertainty about the form of the primary response function beyond the 20% point ( aips task pbcor ) , we only know that the peak flux of rt19970205 lies in the range of 2.94.4mjy . as noted above the source is likely extended and this may account for the variation in snr ( from @xmath55 to 7.6 ) .
@xcite presented a survey for 5ghz radio transients at low galactic latitudes .
most of the data were reduced in near - real time ( 2hr delay ) and transient candidates were followed up with radio , visible light and x - ray instruments .
the authors reported a single transient candidate , j@xmath56 with a peak specific flux of @xmath57mjy ( 5.8@xmath11 ) and with no obvious optical counterpart .
we have re - analyzed this data set and confirm that the snr is 5.8 . considering that the search resulted from inspecting @xmath58 independent beams (  [ sec : numberindependentbeams ] ) and assuming gaussian statistics , we find that the probability of the highest event found in these many independent beams is attributable to chance or noise is 3.6% ( see  [ sec : a ] ) .
therefore , we advocate that this candidate is not a real event .
as noted in  [ sec : introduction ] the ten transients reported in b07 were remarkable for their apparent abundance ( 1.5 per square degree , at any given epoch ) _ and _ the lack of quiescent optical counterparts . in
[ sec : findings ] we found that five of the b07 transients are artifacts arising out of ( rare ) data acquisition problems or imaging artifacts ( see table  [ tab : listoffields ] ) .
we find that the 2-month transient rt20010331 is not well detected .
we argue that the single transient reported by @xcite is also not a robust detection .
in  [ sec : findings ] we present our analysis for the remaining four transients ( all of which are single - epoch transients ) : rt19920826 , rt19970205 , rt19970508 and rt19990504 .
we find snrs which are typically 1-@xmath11 below that reported in b07 .
the revised and b07 snrs are given in table  [ tab : listoffields ] .
we note that an independent pipeline written in parseltongue @xcite , a python interface to aips , and run on three of these four sources confirms the lower snr values found here @xcite : rt 19920826 ( snr=6.0 ) , rt 19970528 ( snr=4.4 ) , and rt 19990504 ( snr=4.0 ) .
the lowering of snr ( from between 7 and 8 to between 5 and 6 ) has a pernicious effect when the number of independent beams which were search is included . in  [ sec : numberindependentbeams ] we estimate this number to be @xmath59 . in
[ sec : a ] we derive the probability density function for the highest @xmath60 values of @xmath61 gaussian random numbers . in figure
[ fig : pdistpoints ] we plot the density function for the highest value ( @xmath62 ) and the fourth highest value ( @xmath63 ) . as can be seen from this figure  [ fig : pdistpoints ] , if the snrs reported here
are accepted then the global case for the remaining b07 transients is entirely weakened . if , on the other hand , the snrs reported in b07 are accepted then the four transients reported in b07 do argue for a new class of radio transients .    ) and the fourth highest value ( @xmath63 ) of a population of @xmath64 gaussianly distributed random numbers with zero mean and unit variance . the probability density function plotted here is discussed in  [ sec : a ] and the justification for @xmath61 can be found in  [ sec : numberindependentbeams ] .
the dots represent the snrs of the following single - epoch transients discussed in table  [ tab : listoffields ] .
( from left to right ) : rt19905054 , rt19920826 , rt19970205 and rt19970228 .
( top ) : snrs as reported in this paper .
( bottom ) : snrs as reported in b07 .
note that the vertical location of the dots is arbitrary . for details of these density distributions
see  [ sec : a ] . , title="fig:",width=259 ] ) and
the fourth highest value ( @xmath63 ) of a population of @xmath64 gaussianly distributed random numbers with zero mean and unit variance .
the probability density function plotted here is discussed in  [ sec : a ] and the justification for @xmath61 can be found in  [ sec : numberindependentbeams ] .
the dots represent the snrs of the following single - epoch transients discussed in table  [ tab : listoffields ] .
( from left to right ) : rt19905054 , rt19920826 , rt19970205 and rt19970228 .
( top ) : snrs as reported in this paper .
( bottom ) : snrs as reported in b07 .
note that the vertical location of the dots is arbitrary . for details of these density distributions
see  [ sec : a ] . , title="fig:",width=259 ]    the above approach of using a fixed threshold for all epochs does not result in optimal detection . in particular , the threshold for a low resolution survey is lower than that for a higher resolution survey ( since the latter has a correspondingly larger number of synthesized beams ) .
b07 addressed this problem by requiring that the probability of a false detection ( pfd ) in an individual epoch was constant and less than @xmath65 where @xmath65 is the total number of images . with this approach ,
the expectation number of false detections is 1 for the entire survey . applying the b07 method we find the following pfds : rt19920826 ( log(pfd)=@xmath66 ) ; rt19970205 ( @xmath67 ) ; rt19970528 ( @xmath68 ) ; rt19990504 ( @xmath69 ) . with this more refined approach
only rt19920826 and 19990504 survive .
however , for reasons discussed in  [ sec : rt19990504 ] we have misgivings about rt1990504 .
an entirely different approach ( and in some ways orthogonal to the above snr based approach ) is to look at the angular distribution of the transient sources with respect to the primary axis .
basic interferometry theory informs us that the dirty image is simply the fourier transform of the visibility data .
as such the radiometric noise in the dirty image should be independent of the angular offset from the phase center .
in contrast , the point source sensitivity decreases as one goes away from the pointing center and this is governed by the primary beam response ( assuming that the spectral resolution of the survey is high enough that the delay beam is larger than the primary beam ) .
thus , once the minimum snr for detection is fixed , cosmic sources should be concentrated towards the pointing direction whereas noise spikes ( masquerading as threshold point sources ) should be uniformly distributed .    in  [ sec : primarybeam ] we derive the expected distribution of cosmic sources as a function of the angular offset . in figure
[ fig : primarybeampoints ] we plot the expected cumulative distribution and also the angular offset of the four sources which are not artifacts but whose snr seems to be under dispute , namely rt19920826 (  [ sec : rt19920826 ] ) , rt19970205 (  [ sec : rt19970205 ] ) , rt19970528 (  [ sec : rt19970528 ] ) and rt19990504 (  [ sec : rt19990504 ] ) .
from this figure one can see that only rt19920826 lies in the expected region whereas the remaining three are collectively improbable .    .
the sources are represented by dots and are ( from left to right ) : rt19920826 , rt19970528 , rt19970205 , and rt19990504 .
the points are deliberately placed at cumulative probability of 1% .
rt19920826 is firmly within the region where one naturally expects cosmic sources.,width=259 ]    in summary , two different statistical tests , one based on snr and the other making use of the spatial signature provided by the primary beam , suggest that of the remaining four sources detected at threshold , only one , namely , rt19920826 is a good detection .
thus a simple interpretation of our re - analysis is that the rate of b07 transients is considerably lower than that reported by b07 , perhaps an order of magnitude smaller .
we acknowledge that the discrepancy between the analysis presented in b07 and the analysis presented here ( and re - investigated ) is disturbing . in the previous section (  [ sec : remainingfourtransients ] )
we investigate the reasons for the discrepancy in the snrs and variously find possible and plausible causes : flagging of data , the choice of data reduction package ( aips versus _ miriad _ ) and the specific method used to compute snrs .
however , none of these explanations are satisfactory .
we are continuing this investigation but at the present time we consider this topic to beyond the scope of the paper .
in contrast , the sources which we find as artifacts have ready explanations ( see table  [ tab : listoffields ] ) .
one source is a result of the file header containing a pointing direction of the previous pointing .
another is due to a systematic associated with local signals ( rfi ) .
these signals do not have the natural fringe rate of cosmic sources and appear as candidates close to the phase center .
two are side - lobes of a stronger sources . suffice to say that such rare errors will be found if one inspects sufficiently large number of beams !
in figure  [ fig : rates ] we plot the areal densities of three known transients ( sn1998bw - like , type - ii rsn and swift j1644 + 57-like ) and that of two expected classes ( ns - ns mergers and orphan afterglows of long duration grbs ) of transients .
the areal density of these five classes is also summarized in table  [ tab : listoftrans ] .
we briefly discuss each of these five classes of transients below .
the traditional extra - galactic radio transients are type ii radio supernovae .
the rate that we present here is based on the single radio sn detected in a blind radio survey by @xcite and it agrees with an independent estimate by @xcite . as illustrated by the example presented in @xcite , the main advantage in the search for radio sne in a blind survey is the unique view that it provides to the otherwise hidden population of heavily obscured sne .
the observed rate of other types of radio sne ( e.g. ordinary type ib / c sne ) are considerably lower than that of type ii supernovae and are not discussed further .
radio emission is expected from both classical long duration gamma - ray bursts ( grbs ) as well as the more abundant but less luminous grbs exemplified by grb980425 associated with the energetic sn ic supernova , sn1998bw @xcite .
the radio emission is far brighter than that of ordinary core collapse sn ( ii , ib , ic ) and the increased volume makes up for the intrinsically smaller birth rate .
the recent discovery of the energetic supernova sn2009bb @xcite demonstrates that radio surveys can find such sources without resorting to high energy ( gamma - ray ) missions .
@xcite estimated the number of afterglows from classical grbs and whose explosion axis is directed away from us ( `` orphan '' afterglows ) .
the expected rate depends strongly on the poorly constrained @xmath70-ray beaming . on one hand
this makes any rate prediction uncertain . on the other hand
, even a non - detection by the kind of survey that we discuss below will provide an independent and unique constraint on the average opening angle of long gamma - ray bursts , their true rate and total energy output @xcite . the areal density in figure  [ fig : rates ]
is derived using a typical opening angle of 10 degrees .    a surprising and
apparently an important development in the field of radio transients took place just this year with the discovery of a radio transient associated with the nucleus of a modest size galaxy .
the source , swiftj1644 + 57 was initially detected as a hard x - ray transient @xcite .
subsequent follow up found a bright , compact self - absorbed radio counterpart , localized at the center of a normal galaxy at @xmath71 @xcite .
the current view is that swiftj1644 + 57 arose from a relativistic jet produced when a star was tidally disrupted as it passed too close to an otherwise dormant super - massive nuclear black hole ( e.g. * ? ? ?
shortly , thereafter , a second candidate non - thermal tidal disruption event ( tde ) was recently proposed @xcite .
events such as these give us an opportunity to study the activity of 10@xmath72 - 10@xmath73 m@xmath74 supermassive black holes in otherwise normal galaxies .    the areal density in figure  [ fig : rates ]
is calculated assuming an observed rate of 0.2 @xmath75 swiftj1644 + 57-like events and a gamma - ray beaming factor of @xmath76 @xcite .
nominally , swiftj1644
+ 57-like sources appear to be the most frequent extra - galactic transients that will be found in radio transient searches .
we acknowledge that the uncertainty of both the observed rate and the gamma - ray beaming is high and the true rate may be significantly different .
now we come to the most uncertain as well as potentially the most important extra - galactic radio transient  the merger of two neutron stars ( or a black hole and a neutron star ) .
it is generally accepted ( or expected ) that short hard bursts are on - axis explosions of these mergers @xcite . as in long duration grbs
, radio emission is expected by afterglow ( on - axis or orphan ) .
the rates are highly uncertain because there are very few observations of short hard grbs .
thus there still continues to be a debate about the geometry of these explosions ( `` jetted '' or not ) .
next , while the expected radio emission is straightforward to estimate ( subject to the usual parametric uncertainties of the energy fractions of relativistic electrons and magnetic field ) an additional uncertainty is the density of the ambient gas ( which is necessary for the production of the afterglow emission ) .
regardless of the uncertainty whether neutron star mergers are the sources of short grbs or not , a substantial sub- and mildly relativistic outflow is expected to be ejected during the merger . @xcite
estimate radio emission from these outflows .
the areal density in figure  [ fig : rates ] is calculated based on their estimates for ns - ns mergers @xcite ] , assuming a ns - ns merger rate of 300 @xmath77 and that any merger ejects @xmath78 erg of a mildly relativistic outflow .
we note that @xcite suggested that rt19870422 was the radio emission from the remains of a neutron star merger .
however , as noted in
[ sec : findings ] this source is an artifact .
there are sound reasons to continue the exploration of the dynamic radio sky .
radio searches are an ideal way to discover core - collapse supernovae embedded in or behind dusty regions .
the discovery of sn2009bb shows that large radio searches can find urgently needed additional examples of nearby low luminosity grbs .
next , the many rewards of radio follow up observations of swift1644 + 57 ( accurate localization , energetics , beaming , outflow velocity ) show the tremendous diagnostic power of radio observations of this entirely new class of extragalactic transients .    as exciting as these developments are
, the search for new classes of radio transients has involved several false starts .
the euphoria that followed the discovery of a highly dispersed ( and therefore argued to be of extragalactic origin ) millisecond burst @xcite was rapidly diminished by the discovery of many such bursts , presumably of terrestrial origin @xcite ; but see @xcite .
similarly , a long - duration transient found by @xcite and @xcite was later traced to a glitch in the vla on - line data taking system @xcite , apparently affecting 0.29% of all first survey pointings @xcite .
finally , our re - analysis ( see  [ sec : rt19870422 ] ) shows that the claim of late time radio emission from neutron star coalescence @xcite is premature .
here we focused on the potentially new class of long duration radio transients reported in b07 .
we rule out six of the ten transients and cast doubts on some of the remaining ones .
however , even with one or two survivors the long duration transients of b07 remain of great interest .
first , even with a diminished number of transients the implied areal rate of b07 transients would be comparable to the recently established class of swiftj1644 + 57 transients . however , _ unlike any other long duration radio transient _ ( to wit , supernovae ; active stars ; tidal disruption events ; gamma - ray burst afterglows , beamed or otherwise ) the b07 transients are remarkable for the absence of a quiescent optical counterpart .
the event rt19920826 survives two independent tests . as such it useful to speculate on the origin of this transient .
the absence of an optical ( b07 ) and near ir quiescent source could mean one of two origins .
the event is extra - galactic in origin but the host galaxy is faint enough not to have been detected ( as does happen for a few percent of long duration grb host galaxies ) or that it is offset from a host galaxy ( as in the case for a few short hard bursts ; @xcite ) .
alternatively , the event is a galactic neutron star and we have to then appeal to an optically invisible galactic ( neutron star ) population @xcite .
clearly , a new survey which can net a dozen of such sources ( but brighter ! ) would help resolve the origin .
for the galactic hypothesis one expect no quiescent counterpart even at hst sensitivity whereas detectable galaxies , in the majority , are expected for the extra - galactic hypothesis .
there are great gains in the discovery of new classes of radio transients but at the same time the path to true success is littered with false starts or easy speculations .
the way forward must incorporate the lessons learnt from the false starts .
we elaborate on this conclusion below .    to start
with we believe that the field of radio transients ( at least in the decimeter band ) is sufficiently mature that any new survey which just sets an upper limit relative to the known population is of marginal value .
future surveys have to be sufficiently deep and cover large enough sky that success ( i.e. detection of a few to many transients ) is assured . in our opinion , this means that a survey should be designed to find at least one or more swiftj1644 + 57-like transient ( figure  [ fig : rates ] ) .    next , timely and multi - wavelength follow - up is essential .
for example , the transients reported in @xcite , @xcite and @xcite have plausible origins as supernovae and swiftj1644 + 57-like sources ( see figure  [ fig : rates ] ) .
however , the lack of timely follow - up or deep multi - wavelength followup of these events preclude us from coming to a definitive conclusion .
finally , the search should be restricted to sources with a high level of significance .
this certainly means paying attention to the large number of beams searched .
however , at low thresholds and with large number of beams ( cf  [ sec : findings ] ) it would be prudent to set thresholds beyond mere statistical considerations . a threshold of 9 or even @xmath79 may be appropriate .
alternatively , an immediate verification of a transient by deeper observation or a confirmation by observations at other wavelengths would allow detection of transients closer to threshold .
we start with a discussion of two recent developments .
@xcite undertook an ambitious program similar in spirit to b07 , namely the investigation of fields surrounding vla calibrators .
sources with duration between 4 and 40d were searched for .
the total integration time was 435hr .
no transient source in the ghz range with flux greater than 8mjy was found .
the authors place an upper limit to the areal density of 0.032deg@xmath7 . assuming @xmath80 scaling this areal density is 4.4deg@xmath7 and is well above the b07 rate .
the first survey imaged the sky in an hexagonal grid , in which each position in the survey footprint was observed on average 34 times @xcite .
@xcite used this fact to construct light curves of sources detected in individual first survey snapshots .
they identified 1627 variable candidates with variability exceeding @xmath81 .
this effort is probably the largest variable and transient survey ever carried out .
one disadvantage of such a survey for transient identification is that the co - added images are not much deeper than a single epoch image .
this make it hard to tell if an apparent transient source is really a transient or just a variable source that exceeded the detection threshold in one of the epochs .
the limitations discussed above lead us to suggest a new evla survey specifically tailored to systematically explore the radio sky . following that we review a far more ambitious survey  the vast key project on askap .
for the discussion below we will adopt the rates summarized in table  [ tab : listoftrans ] .
the rates are specified to a flux density of 0.3mjy and are extrapolated to higher flux densities as n(@xmath82 .
a moderately ambitious survey with the evla can result in great progress .
this survey has two virtues .
one , the evla offers excellent spatial resolution .
next , the evla is fully commissioned and is working to specifications .
specifically consider a 100 square degree survey undertaken in the 24ghz band .
an integration time of 85s results in a sensitivity of 0.3mjy ( 10-@xmath11 ) .
a single epoch covering 100 square degrees would require 50 hours .
as can be seen from figure  [ fig : phase ] ( taken from table  [ tab : listoftrans ] ) such a survey would have to explore a variety of time scales to probe the emerging classes of transient .
fifteen epochs could reasonably cover the range of a week to years .
noting the great importance of multi - wavelength imaging data , such a survey would sensibly focus on regions of sky where considerable multi - wavelength data ( including radio ) exists .
one such region is , for example , the sdss equatorial stripe @xcite .
furthermore , the high instantaneous sensitivity of the evla makes rapid follow up of newly transients possible .
after the first four epochs the reference images will be twice as deep as the survey field .
a single new epoch would then yield about ten swiftj1644 + 57-like sources and four supernovae .
ten such images may find a new example of an sn1998bw - like event , several clear examples of orphan afterglows and have an excellent chance of finding the first examples of neutron star mergers .
we note that these different classes of objects have different characteristics , both in duration ( figure  [ fig : phase ] ) and also in host magnitudes and location with respect to host galaxy .
therefore , it is possible to distinguish between different types of objects .    as note
earlier , rapid verification of a transient ( either by additional and deeper radio observations or observations at other wavelengths ) can reduce the requirement for a high detection threshold .
this would then require close rapid reduction
well within the reach of modern computers .
the variable &  slow transient ( vast ) is an approved key project of askap .
the vast - wide survey aims to survey in the 1.2ghz band about 10,000 square degrees every day for 2 years . with 40-s integrations
the expected single - epoch rms is 0.5mjy ( vast memo#1 ) .
since the survey is planned to be undertaken daily the reference image will be built up quite rapidly .
a comparison of data obtained at a new epoch when compared with the built - up reference image can be expected to detect about six swiftj1644 + 57-like sources .
type ii sn and swiftj1644 + 57-like sources are relatively long lived and
so one could consider averaging the daily images ( to say 10 days ) .
the resulting summed data has a sensitivity of 0.16mjy .
the expected number of transients per such summed image is nominally 32 ( type ii sn ) , 82 ( swiftj1644 + 57 ) , 8 ( orphan afterglows ) and 4 ( neutron star coalescences ) .    to sum up , the dynamic radio sky remains a rich area for exploration .
based on what we know about the areal rates for the _ known _ transient sources , future synoptic radio imaging surveys are expected to yield substantial numbers of exotic transients .
such surveys will also provide the definitive test for the b07 population .
daf thanks jim condon and alicia soderberg for important discussions early on in this project .
we thank steve croft for a most careful reading of the paper and j. condon for making several insightful suggestions .
the national radio astronomy observatory is a facility of the national science foundation operated under cooperative agreement by associated universities , inc .
srk thanks the department of astronomy , university of wisconsin at madison for their hospitality .
eoo is supported by an einstein fellowship and nasa grants .
srk s research in part is supported by nasa and nsf .
this research has made use of data from the university of michigan radio astronomy observatory which has been supported by the university of michigan and by a series of grants from the national science foundation , most recently ast-0607523 .
this research has made use of nasa s astrophysics data system .    ,
j. , et  al . 2010 ,
classical and quantum gravity , 27 , 173001    bannister , k.  w. , murphy , t. , gaensler , b.  m. , hunstead , r.  w. , & chatterjee , s.  2011 , , 412 , 634    bannister , k.  w. , murphy , t. , gaensler , b.  m. , hunstead , r.  w. , & chatterjee , s.  2011 , erratum , july 2011    becker , r.  h. , white , r.  l. , & helfand , d.  j.  1995 , , 450 , 559    bell , m.  e.  2011 , `` the low frequency array and the transient and variable radio sky '' , phd thesis , southampton university .
bell , m.  e. , fender , r.  p. , swinbank , j. , et al .
2011 , , 415 , 2    berger , e.  2009 , , 690 , 231    berger , e. , kulkarni , s.  r. , frail , d.  a. , & soderberg , a.  m.  2003 , , 599 , 408    bloom , j.  s. , et al .
2009 , arxiv:0902.1527    bloom , j.  s. , et al .
2011 , science , 333 , 203    booth , r.  s. , de blok , w.  j.  g. , jonas , j.  l. , & fanaroff , b.  2009 , arxiv:0910.2935    bower , g.  c. , saul , d. , bloom , j.  s. , bolatto , a. , filippenko , a.  v. , foley , r.  j. , & perley , d.  2007 , apj , 666 , 346    bower , g.  c.  2011 , , 732 , l12    bower , g.  c. , whysong , d. , blair , s. , croft , s. , keating , g. , law , c. , williams , p.  k.  g. , & wright , m.  c.  h.  2011 , arxiv:1107.1517    burke - spolaor , s. , bailes , m. , ekers , r. , macquart , j .-
p . , & crawford , f. , iii 2011 , , 727 , 18    burrows , d.  n. , kennea , j.  a. , ghisellini , g. , et al .
2011 , , 476 , 421    cenko , s.  b. , fox , d.  b. , moon , d .- s .
, et al .
2006 , , 118 , 1396    cenko , s.  b. , et al .
2011 , arxiv:1107.5307    chatterjee , s. , murphy , t. , & vast collaboration 2010 , bulletin of the american astronomical society , 42 , # 470.12    cotton , w. d. & perley , r. 2010 , obit development memo series no .
17    croft , s. , bower , g.  c. , keating , g. , law , c. , whysong , d. , williams , p.  k.  g. , & wright , m.  2011 , , 731 , 34    croft , s. , tomsick , j.  a. , & bower , g.  c.  2011 , arxiv:1107.5039    dewdney , p.  e. , hall , p.  j. , schilizzi , r.  t. , & lazio , t.  j.  l.  w.  2009 , ieee proceedings , 97 , 1482    ellingson , s.  w. , clarke , t.  e. , cohen , a. , et al .
2009 , ieee proceedings , 97 , 1421    galama , t.  j. , vreeswijk , p.  m. , van paradijs , j. , et al .
1998 , , 395 , 670    gal - yam , a. , et al .  2006 , , 639 , 331    gehrels , n.  1986 , , 303 , 336    gehrels , n. , ramirez - ruiz , e. , & fox , d.  b.  2009 , , 47 , 567    giannios , d. , & metzger , b.  d.  2011 , , 1137    gregory , p.
c. , & taylor , a.  r.  1986 , aj , 92 , 371    gregory , p.  c. , capak , p. , gasson , d. , & scott , w.  k.  2001 , galaxies and their constituents at the highest angular resolutions , 205 , 98    greisen , e.  w.  2003 , information handling in astronomy - historical vistas , 285 , 109    hodge , j.  a. , becker , r.  h. , white , r.  l. , richards , g.  t. , & zeimann , g.  r.  2011 , , 142 , 3    johnston , s. , taylor , r. , bailes , m. , et al .
2008 , experimental astronomy , 22 , 151    jonas , j.  l.  2009 , ieee proceedings , 97 , 1522    , v. , et  al .
2004 , , 601 , l179    keane , e.  f. , kramer , m. , lyne , a.  g. , stappers , b.  w. , & mclaughlin , m.  a.  2011 , , 415 , 3065    kettenis , m. , van langevelde , h.  j. , reynolds , c. , & cotton , b.  2006 , astronomical data analysis software and systems xv , 351 , 497    krding , e. , rupen , m. , knigge , c. , fender , r. , dhawan , v. , templeton , m. , & muxlow , t. 2008 , science , 320 , 1318    kulkarni , s.  r. , frail , d.  a. , wieringa , m.  h. , et al .
1998 , , 395 , 663    law , n.  m. , kulkarni , s.  r. , dekany , r.  g. , et al .
2009 , , 121 , 1395    lazio , j.  2008 , astronomische nachrichten , 329 , 330    levan , a.  j. , et al .
2011 , science , 333 , 199    levinson , a. , ofek , e.  o. , waxman , e. , & gal - yam , a.  2002 , apj , 576 , 923    lien , a. , chakraborty , n. , fields , b.  d. , & kemball , a.  2011 , arxiv:1107.0775    lonsdale , c.  j. , cappallo , r.  j. , morales , m.  f. , et al .
2009 , ieee proceedings , 97 , 1497    lorimer , d.  r. , bailes , m. , mclaughlin , m.  a. , narkevic , d.  j. , & crawford , f.  2007 , science , 318 , 777    metzger , b.  d. , et al.2010 , , 406 , 2650    metzger , b.  d. , & berger , e.   2011 , arxiv:1108.6056    nakar , e.  2007 , , 442 , 166    nakar , e. , piran , t. , & granot , j.  2002 , , 579 , 699    nakar , e. , & piran , t.  2011 , , 478 , 82    , r. , piran , t. , & shemi , a. 1991 , , 379 , l17    nissanke , s.  m. , sievers , j.  l. , dalal , n. , & holz , d.  e.  2011 , arxiv:1105.3184    sault , r.  j. , teuben , p.  j. , & wright , m.  c.  h.  1995 , astronomical data analysis software and systems iv , 77 , 433    thyagarajan , n. , helfand , d.  j. , white , r.  l. , & becker , r.  h.  2011 , arxiv:1107.5901    ofek , e.  o. , breslauer , b. , gal - yam , a. , frail , d. , kasliwal , m.  m. , kulkarni , s.  r. , & waxman , e.  2010 , , 711 , 517    ofek , e.  o. , & frail , d.  a.  2011 , , 737 , 45    ofek , e.  o. , frail , d.  a. , breslauer , b. , kulkarni , s.  r. , chandra , p. , gal - yam , a. , kasliwal , m.  m. , & gehrels , n.  2011 , arxiv:1103.3010    oosterloo , t. , verheijen , m.  a.  w. , van cappellen , w. , et al .
2009 , proceedings of wide field astronomy technology for the square kilometre array ( skads 2009 ) .
november 2009 .
chateau de limelette , belgium .
, e.  s. 1991 , , 380 , l17    perley , r.  a. , schwab , f.  r. , & bridle , a.  h.  1989 , synthesis imaging in radio astronomy , 6 ,    rau , a. , et al .
2009 , , 121 , 1334    rossi , e.  m. , perna , r. , & daigne , f.  2008 , , 390 , 675    rttgering , h. , de bruyn , a.  g. , fender , r.  p. , et al .  2003 , texas in tuscany .
xxi symposium on relativistic astrophysics , 69    scott , w.  k. 1996 , phd thesis , _ compact source variability at centimeter wavelengths _ , university of british columbia    soderberg , a.  m. , et al .
2006 , , 442 , 1014    soderberg , a.  m. , et al .
2010 , , 463 , 513    sylvestre , j.  2003 , , 591 , 1152    thyagarajan , n. , helfand , d.  j. , white , r.  l. , & becker , r.  h.  2011 , arxiv:1107.5901    totani , t. , & panaitescu , a.  2002 , , 576 , 120    de vries , w.  h. , becker , r.  h. , white , r.  l. , & helfand , d.  j.  2004 , , 127 , 2565    van velzen , s. , koerding , e. , & falcke , h.  2011 , arxiv:1104.4105    weiler , k.  w. , panagia , n. , sramek , r.  a. , et al .
2010 , , 81 , 374    welch , j. , backer , d. , blitz , l. , et al .
2009 , ieee proceedings , 97 , 1438    zauderer , b.  a. , et al .
2011 , arxiv:1106.3568    de zotti , g. , massardi , m. , negrello , m. , & wall , j.  2010 , , 18 , 1
in this section , our goal is to compute the probability density function of the @xmath60th highest value of @xmath83 .
let @xmath84 be the probability density function of @xmath85 with @xmath86 being the cumulative function .
we denote the @xmath87 highest value by @xmath88 .
thus the maximum of the series of measurements is @xmath89 and the minimum value is @xmath90 .
let @xmath91 be the probability density function of @xmath92 .
this means that at least one of the measurements lies is in the range @xmath93 $ ] .
the probability density for this event is @xmath94 .
next , then @xmath95 measurements must lie below this range and @xmath96 above this range .
the probability for any value to be smaller than @xmath88 is @xmath97 and the probability for a value to higher than @xmath88 is is @xmath98 .
this is now a binomial distribution with @xmath99 total values .
thus the probability density function for @xmath88 is @xmath100^{m-1}.\          \label{eq : rho}\ ] ] the first combinatorial factor of @xmath61 accounts for the possibility that @xmath88 can occupy any position in the sequence .
the second combinatorial factor , @xmath101 account for the combinations satisfying the condition that @xmath95 values lie below @xmath88 and @xmath96 lie above @xmath88 .
for both the maximum ( @xmath62 ) and minimum ( @xmath102 ) equation  [ eq : rho ] simplifies to that expected from basic considerations .
now let us consider the specific case where @xmath103 follows gaussian statistics : @xmath104 where @xmath103 is normalized in units of @xmath11 .
the probability that an event is extreme or lies within the range @xmath105 is @xmath106 the probability that an event is extreme in only one sense , maximum or minimum , and lies outside the range @xmath107 $ ] ( say ) is thus @xmath108=                   1-\frac{1}{2}{\rm erfc}\big(\frac{h}{\sqrt{2}}\big )          \end{aligned}\ ] ] where @xmath109 .
consider the case where @xmath110 ( say @xmath111 or more ) and @xmath60 is small , say 10
. then we can approximate @xmath112 , @xmath113 , ... , @xmath114 . furthermore using the approximation , @xmath115
, we find in the limiting case where @xmath88 is greater than a few [ so that @xmath116 : @xmath117          \big[\frac{n}{2}{\rm erfc}\big(\frac{h_m}{\sqrt{2}}\big)\big]^{m-1}.          \end{aligned}\ ] ]
provided that there is sufficient spectral resolution , the response of an interferometer to a source follows the antenna response function ( `` primary beam '' ) .
we will assume that the response function is azimuthally symmetric and specified by @xmath118 where @xmath119 is the angular offset from the pointing axis .
let the areal density of sources with flux density greater than @xmath120 be a power law , say , @xmath121 . for euclidean geometry and most reasonable luminosity functions , @xmath122 .
next , we note that detection is really finding sources at a given snr and above .
fortunately the noise distribution for an interferometric image is uniform .
thus a source with a given flux density will have an snr , @xmath123 , that scales with the primary beam response , @xmath124 .
the number of sources above a certain snr and contained outside an angular radius of @xmath125 is @xmath126    for the specific case of a gaussian beam , @xmath127          \label{eq : gaussian}\ ] ] where the traditional `` full width at half maximum '' ( fwhm ) is @xmath128 . substituting equation  [ eq : gaussian ] into equation  [ eq : ns ] we obtain @xmath129 half
the sources will be detected outside the radius @xmath130 .
for @xmath122 we obtain @xmath131 .
the expression @xmath132 is more useful .
a plot of @xmath133 can be found in figure  [ fig : primarybeamdistribution ] .     from the pointing aixs and assuming @xmath122 .
the response curve is normalized by insisting that the integral of the curve ( from @xmath134 to @xmath135 ) is unity .
fifty percent of the sources are within @xmath136 ( marked by a pentagram ) and 97% within @xmath137 ( marked by a cross ; at this radius the primary beam gain is only 0.1 relative to that on - axis ) .
, width=288 ]
it is not unusual to find no source after undertaking a survey .
we wish to determine an upper limit to the number of sources that were being searched .
the number detected is given by poisson statistics .
the probability of finding @xmath138 sources is then given by @xmath139 where @xmath9 is the poisson parameter and equal to @xmath140 .
our goal is to determine the maximum value of @xmath9 given a non - detection . as the value of @xmath9
is increased the probability of detection , by which we mean the probability of detecting one or more events , also increases .
this probability is @xmath141 which we note is @xmath142 .
this probability can be set to the desired confidence level , @xmath143 and thence @xmath144 the reader with a stronger physical bent may find the complement , @xmath145 , more appealing .
regardless , we find @xmath146 $ ] at a confidence level of [ 95% , 99% , 99.9% ] .
here we compute @xmath61 , the number of independent beams for the vla data that went into the analysis of b07 and @xcite . for b07 , a circular region with a radius of two times
the half - power radius was searched for each epoch .
the number of independent beams per epoch is the ratio of that area to the area of the synthesized beam in that particular epoch . for individual epochs ,
this value ranged from as small as @xmath147 to as large as @xmath111 .
the total numbers of independent beams for the 5 and 8.4 ghz data are @xmath148 and @xmath149 , respectively .
the smaller number of independent beams for 8.4ghz are due to the tapering of the visibility data , increasing the typical synthesized beam size at 8.4ghz .
for @xcite the search was made for transients in 4.65radius circular region .
the fwhm of the synthesized beam is @xmath150arc - seconds . though the goal was 16 epochs per pointing we only achieved an average of 15.7 images . with 141 epochs
we derive @xmath151 .