radio pulsars show a wide variety of modulations of their pulse amplitudes , including bursts and nulls , that affect their detectability in surveys .
phenomena seen in some pulsars include short - period nulling , in which a pulsar is not detected for several pulse periods , only to reappear with full strength ( @xcite ) ; eclipses , in which a companion star or its wind or magnetosphere absorbs or disperses the pulsar signal ( @xcite,@xcite,@xcite,@xcite ) ; long - term nulling or intermittent behavior , in which a pulsar is quiescent for days or weeks ( @xcite ) ; and rotating radio transients ( rrats , @xcite ) , pulsar - like objects from which only occasional radio bursts are detected .
this paper describes analysis of a large - scale survey using the arecibo telescope that is sensitive to both periodic and aperiodic signals .
rrats were first discovered in archive parkes multibeam survey data ( @xcite ) .
eleven objects , with periods ranging from 0.7 to 7  s and pulse widths of @xmath0  ms , were found using a single pulse search algorithm ( @xcite )
. the longer periods of rrats compared with the general pulsar population suggest similarities with the x - ray populations of x - ray dim isolated neutron stars ( xdinss ) and magnetars .
rrat j1819@xmath11458 has been detected at x - ray energies ( @xcite ) with properties that are similar to those of xdinss and high magnetic field radio pulsars .
a different type of pulse modulation is observed in the case of pulsars emitting giant pulses .
such pulses are tens to thousands of times brighter and an order of magnitude or more narrower than the average pulse ( see @xcite for an overview ) .
giant pulses from the crab pulsar have substructure on timescales of 2  ns ( @xcite ) , and psr  b1937 + 21 emits giant pulses as narrow as 16  ns ( @xcite ) .
giant micro - pulses from the vela pulsar have widths @xmath2s ( @xcite ) , and the slowly rotating pulsars psr  b1112 + 50 , psr  b0031@xmath107 , and psr  j1752 + 2359 occasionally emit bright pulses which are @xmath3  ms wide , @xmath4 times narrower than the average pulse ( @xcite ) .
the detection of giant pulses is a potentially powerful method for finding extragalactic pulsars too distant for their normal emission to be detectable by periodicity searches ( @xcite )
.    a variety of energetic phenomena other than pulsar emission can give rise to fast transients potentially detectable in radio pulsar surveys .
within the solar system , transient radio events may be generated by energetic particles impacting the earth s atmosphere , solar flares , and decameter radio flares originating in jupiter s atmosphere .
analogously to the latter , extrasolar planets with strong magnetic fields are expected to be detectable in the 10 - 1000  mhz range ( @xcite , @xcite , @xcite ) .
magnetic activity on the surfaces of brown dwarfs and particle acceleration in the magnetic fields of flare stars are also known radio flare progenitors ( @xcite , @xcite , @xcite , @xcite ) .
gamma - ray bursts are predicted to have detectable radio emission at @xmath5  mhz ( @xcite , @xcite ) , and radio flares have been observed from some x - ray binaries ( @xcite , @xcite ) . among the most energetic and exotic events in the universe ,
supernovae , merging neutron stars and coalescing black holes may produce wide - band radio bursts detectable at extragalactic distances ( @xcite ) .    in this paper
we describe an ongoing survey for pulsars and transient radio sources with the arecibo telescope .
the survey addresses outstanding questions about the nature and emission mechanisms of intermittent radio sources . in   [ section_surveyparams ]
we present the palfa survey parameters , and in   [ section_spmethods ]
we describe the single pulse search and radio frequency interference excision algorithms which are part of the survey data processing pipeline .
section  [ section_palfastats ] contrasts palfa detection statistics on known pulsars and new discoveries , and
[ section_biases ] examines selection effects influencing the detection and classification of transient sources . in
  [ section_intermittency ] , we apply an intermittency measure method for comparing the efficiency of periodicity and single pulse pulsar searches . in
  [ section_objects ] , we discuss the properties of individual intermittent objects discovered by palfa . in   [ section_faroutmodel ] and
  [ section_constraints ] , we apply constraints derived from the survey sensitivity and results to the detectability of various energetic phenomena expected to emit radio bursts . finally , in   [ section_conclusions ] , we present our main conclusions .
the palfa survey started in 2004 , shortly after the installation of the seven - beam alfa receiver on the arecibo telescope .
the survey searches for pulsars and transients in the inner and outer galactic plane regions accessible to arecibo ( see below ) .
the alfa receiver is well - suited for survey observations , allowing simultaneous data collection from seven fields , each @xmath6 ( fwhm ) across . taking into account the hexagonal arrangement of the beams on the sky and the near sidelobes ,
the combined power - pattern is approximately @xmath7 ( @xcite ) .
we observe a 100 mhz passband centered on 1440  mhz in each of the seven telescope beams .
wideband arecibo pulsar processors ( wapps ; @xcite ) are used to synthesize 256-channel filterbanks spanning these bands at intervals of 64  @xmath8s . during observations ,
full - resolution data are recorded to disk and , in parallel , decimated down to a time resolution of 1  ms and searched for periodic signals and single pulses by a quick - look processing pipeline running in real time at the arecibo observatory ( @xcite ) .
this approach allows for immediate discovery of relatively bright pulsars with periods longer than a few milliseconds .
searching full - resolution data allows detection of millisecond pulsars and narrower single pulses and is done offline at participating palfa institutions as the processing is much more computationally intensive .
table  [ pointing_stats ] lists various alfa system and survey parameters , including the sky area corresponding to processed and inspected data reported on in the present paper and the total sky area observed to date ( see @xcite for a detailed explanation of other parameters ) .
standard observation times are 268  s for inner galaxy pointings ( @xmath9 ) and 134  s for outer galaxy pointings ( @xmath10 ) .
some early observations had a duration of 134  s and 67  s for inner and outer galaxy pointings , respectively .
the quoted system temperature of 30  k is measured looking out of the galactic plane .
the initial threshold of @xmath11 for the single pulse search is used when selecting events from dedispersed time series based on their signal - to - noise ratios only , before any filtering is applied . while there are a significant number of events due to random noise above this threshold , identification of an event as a genuine pulse takes into account not only its signal - to - noise ratio but also the fact that it is detected at a contiguous range of trial dispersion measures , which is in general not true of spurious events .
thus weak pulses can be correctly identified , while they would be excluded if the threshold was set according to gaussian noise statistics .
lccc center frequency ( ghz ) & 1.440 + total bandwidth ( mhz ) & 100 + channel bandwidth ( mhz ) & 0.39 + sampling time ( @xmath8s ) & 64 + nominal system temperature ( k ) & 30 + gain ( k / jy ) & + center pixel & 10.4 + ring pixels & 8.2 + beam width (  ) & + main beam ( 1 pixel ) & 3.5 + main beam plus near sidelobes ( 1 pixel ) & 18 + main beam plus near sidelobes ( 7 pixels ) & @xmath12 + inner galaxy & + observation time per pointing ( s ) & 268 , 134 + observed ( deg@xmath13 ) & 156 + processed ( deg@xmath13 ) & 99 + anticenter & + observation time per pointing ( s ) & 134 , 67 + observed ( deg@xmath13 ) & 119 + processed ( deg@xmath13 ) & 87 + detection threshold ( @xmath14 ) & + single pulse search ( initial threshold ) & 5 + fourier transform periodicity search & 7.5      here we compute the maximum distance at which sources of a given luminosity are be detectable by the palfa survey , @xmath15 .
we then compare the sensitivity of the palfa survey to previous work done with the parkes multibeam system ( @xcite , @xcite ) .
the rms noise in a radio transient search , where the effective observation time is equal to the pulse width , is _
n = , where @xmath16 is the system - equivalent flux density , @xmath17 is the number of polarization channels summed and @xmath18 is the bandwidth .
a pulse s observed width @xmath19 may be broadened compared to its intrinsic width @xmath20 by several effects .
after dedispersing the raw data and obtaining a dedispersed time series , there is residual dispersive broadening due to the finite width of a frequency channel and the error of the trial dispersion measure ( dm ) used compared to the actual pulsar dm .
scatter broadening is not correctable and will have a contribution that depends on observing frequency and varies with direction on the sky . in general , w ( w_i^2 + t_dm , ch^2 + t_dm ,
err^2 + t_sc^2)^1/2 , where @xmath21 is the dispersive broadening across a frequency channel of width @xmath22 ( mhz ) for an observing frequency @xmath23 ( ghz ) , @xmath24 is the dispersive broadening due to the difference between the trial and actual dm of the source , and @xmath25 is the scattering broadening . broadening conserves pulse area ,
so that the intrinsic and observed peak flux densities are related through @xmath26 .
if @xmath27 is the detection threshold , the minimum detectable _ intrinsic _ peak flux density is s_p , i , min = ( ) . for a one steradian pulsar radio beam , a source of intrinsic peak luminosity @xmath28 can be detected out to a maximum distance of d_max = ( ) ^1/2 = l_p , i^1/2 ( ) ^1/2 . for a steadily emitting pulsar with period - averaged luminosity @xmath29 and duty cycle @xmath30
, we have @xmath31 .
the amount of pulse broadening depends on system parameters as well as dispersion and scattering , which vary with direction on the sky so that @xmath32 .
we use the ne2001 model of galactic ionized electron density ( @xcite ) to calculate representative results for @xmath15 in the direction @xmath33 , a region of overlap between palfa and the parkes multibeam survey .
[ fig : dmax_vs_lp ] shows @xmath15 vs. @xmath28 detection curves using a threshold @xmath34 for both surveys . for lower luminosities ,
sources are not visible to large enough distances for scattering to affect detectability and the inverse square law dominates the detection curve so that @xmath35 . for larger distances and smaller intrinsic pulse widths ,
scattering and ( residual ) dispersion smearing make pulses increasingly harder to detect and @xmath15 increases more slowly with @xmath28 .
is about twice as large for palfa than for the parkes multibeam survey because of the larger sensitivity of the arecibo telescope .
breaks in the curves correspond to a transition from luminosity - limited to scattering - limited detection for palfa and to a dispersion - limited regime for parkes .
the difference is due to the smaller channel width of palfa ( 0.4  mhz ) compared to parkes ( 3  mhz ) . ]
this section describes processing methods employed by the cornell pulsar search pipeline at the cornell center for advanced computing and the swinburne university of technology .
the presto search code is run independently at west virginia university , university of british columbia , mcgill university and the university of texas at brownsville .
presto uses a similar matched filtering algorithm to the one described below , but a different rfi excision scheme and trial dm list .
we dedisperse raw data with 1272 trial dms in the range @xmath36  pc  @xmath37 . in order to find individual pulses from intermittent sources we operate on dedispersed time series with two time domain algorithms : matched filtering , which has been the standard in single pulse searching so far , and a friends - of - friends algorithm . after single pulse candidates
have been identified we use a stacking method in the time - frequency plane to verify that the pulses are dispersed and the sweep observed across frequencies follows the dispersion relation , as expected for non - terrestrial sources .
certain types of terrestrial signals , e.g. swept - frequency radars , have pulses whose appearance in the time - frequency plane may approximate dispersion by ionized gas , but on closer inspection such pulses generally deviate significantly from the cold plasma dispersion relation .      matched filtering detection of broadened pulses relies on the convolution of a pulse template with the dedispersed time series .
ideally , a pulse template would consist of one or several superimposed gaussian templates , reflecting the diversity of pulsar pulse profiles , some of which exhibit multiple peaks .
in addition , multipath propagation due to scattering adds an exponential tail to the pulse profile .
we approximate true matched filtering by smoothing the time series by adding up to @xmath38 neighboring samples , where @xmath39 , and selecting events above a threshold after each smoothing iteration .
the smoothing is done in pairs of samples at each @xmath40-stage , so that the resulting template is a boxcar of length @xmath38 samples ( @xcite ) .
the sampling time used in palfa observations is 64  @xmath8s and therefore our matched filtering search is most sensitive to pulse widths of 64  @xmath8s to 8.2  ms .
this search strategy is not optimal for single pulses from heavily scattered pulsars because the templates are symmetric , while the scattered pulse shape with an exponential tail is not , and significantly scattered pulses can be wider than @xmath41  ms .
two issues make a complement to matched filtering necessary .
the pulse templates described above have discrete widths of @xmath42 samples by algorithm design and there is decreased sensitivity to pulses with widths that are significantly different .
in addition , the matched filter search output may be dominated by bright , wide pulses from radio frequency interference ( rfi ) . in that case ,
a single rfi burst is detected as an overwhelming number of individual events instead of a single event .
the cluster algorithm is similar to the friends - of - friends search algorithm used to find galaxies in optical images ( @xcite ) , and complements matched filtering by not restricting the width of expected pulse detections .
the dedispersed time series is processed sequentially .
each event above a threshold which is found is designated as the first of a cluster .
a cluster of events is augmented by broadening it to include any adjacent samples above a threshold .
gaps of @xmath43 samples are allowed within a cluster , and palfa data is processed with @xmath44 .
the brightest sample of a cluster is recorded as the event amplitude , and the total number of samples in the cluster as its width .
this approach is less sensitive to weak , narrow pulses but results in significantly fewer spurious events due to rfi .    a limiting factor for the largest pulse width detectable by both the matched filtering and friends - of - friends search algorithms is the fact that data are analysed in blocks much shorter than the complete time series length .
the mean and standard deviation of a block are used for thresholding events found within the block .
this approach minimizes the effect of baseline variations with time scales much longer than typical pulsar pulse widths of a few to a few tens of ms .
the disadvantage is that only pulses with @xmath45 can be detected . in processing palfa data ,
we use blocks of length @xmath46s  @xmath47  s. according to the ne2001 model of ionized gas in the galaxy ( @xcite ) , a scattering broadening time of that magnitude in the inner - galaxy part of the palfa survey corresponds to @xmath48  pc  @xmath37 and a maximum search distance well outside of the milky way for most directions we survey .
however , lines of sight intersecting hii regions can result in large dispersion measures and scattering times and are therefore selected against . any intergalactic scattering that broadens pulses beyond about 0.1
s will also cause events to be selected against .
[ j0627best ] shows standard single pulse search output for the discovery observation of psr  j0627 + 16 (   8.1 ) .
the main panel shows signal - to - noise ratio ( s / n ) vs. dm and time for events with @xmath49 .
the panels on top illustrate event statistics for the observation : from left to right , number of events vs. s / n and dm , and dm vs. s / n . the pulsar detection is manifested as a peak in the histogram of number of events vs. dm . there is a corresponding peak in the dm vs. s / n plot , since the several detected pulses are significantly brighter than background noise in the dedispersed time series .     vs. time and dm with larger circles denoting brighter bursts .
panels on top from left to right show histograms of the number of events vs. s / n and dm , and a scatter plot of event dm vs. s / n .
a fit to the narrow peak of event dm vs. s / n indicates a pulse width @xmath50  ms ( @xcite ) . ]
[ rratquick ] shows events in dm - time space for another palfa single pulse discovery , psr  j1928 + 15 (   8.5 ) . in this case
, three closely spaced bursts were found at @xmath51  s , @xmath52  pc  @xmath37 .
the pulsar is detected at a range of dms , with the signal to noise ratio increasing as trial dm values approach the actual pulsar dm and decreasing as trial dms further on recede from the pulsar dm .
in contrast , events due to interference from terrestrial signals at @xmath53  s span the entire range of trial dms and their signal to noise ratios do not show a significant variation with dm .
if an excess of candidate pulses is identified in the dedispersed time series for a particular trial dm , we use the expected dispersion sweep across the frequency band in order to test if the pulses are from non - terrestrial origin . for a pulse with @xmath54 in the time series , we extract a chunk of raw data centered on its time of arrival and look for a sweep across frequency that follows the dispersion relation ( fig .
[ 4dynspecs ] a , b ) .          in the case of multiple weaker pulses detected in the same beam , we extract a raw data chunk centered on the arrival time of each pulse , stack the chunks and look for a dispersion sweep in the resulting cumulative dynamic spectrum ( fig .
[ 4dynspecs ] c , d ) .
stacking can induce spurious dispersion sweeps in the time - frequency plane even when the data consist of only gaussian - distributed noise .
single pulse events identified in dedispersed time series by definition correspond to sums over a given dispersion path that are above average . when these are used to select chunks for summing
, the statistical fluctuations will build up to show a dispersed pulse in the time - frequency plane that follows the cold - plasma dispersion law perfectly .
we simulated this effect by generating a time - frequency plane of gaussian noise , dedispersing it , and selecting events above threshold from the resulting time series .
we extracted chunks from the fake time - frequency plane centered on each event and stacked them , thus reproducing the procedure used on real data .    selecting events with signal - to - noise above a threshold @xmath55 in the dedispersed time series yields a biased set of noise - only samples in the time - frequency plane . on the other hand ,
when stacking dynamic spectra , points in the stacked dynamic spectrum must have a minimum signal - to - noise @xmath56 for a dispersion sweep to be discernible by eye .
the rms noise in the dedispersed time series after summing @xmath57 channels is a factor of @xmath58 larger than the rms noise in the time - frequency plane .
the rms noise in the stacked dynamic spectrum after stacking @xmath59 chunks is a factor of @xmath60 larger than the rms noise in the time - frequency plane .
we calculate the average deviation from the mean in the time - frequency plane implied by the two thresholds @xmath55 and @xmath61 , equate them , and obtain n_c ( m_c / m_d)^2n_ch .    for the palfa survey and our processing pipeline , @xmath62 and @xmath63
. therefore if @xmath64 data chunks are added there will be a spurious , dispersed pulse in the resulting stacked dynamic spectrum even if the chunks contain only noise .
since the noise decorrelates over @xmath65 samples , the pulse will have a width on the order of @xmath66s . to avoid the induced spurious - event effect ,
the width and s / n of stacked dynamic spectra should be compared with those that would result from noise only .
all of our detections that use stacking satisfy the criteria @xmath67s and @xmath68 .
typically we sum chunks centered on the brightest @xmath69 pulses for a single pulse candidate , which is safely below the @xmath70 limit . in the case of psr
j1909 + 06 ( fig .
[ 4dynspecs]d ) only 2 chunks were added and the dispersed swath in the time - frequency plane is 1  ms wide .      the rfi environment at the arecibo telescope and the 7-beam configuration of the alfa receiver present both challenges and opportunities for rfi mitigation to facilitate searching for single pulses .
palfa survey data are dedispersed with trial dispersion measures of @xmath36  pc  @xmath37 .
rfi pulse intensity typically peaks at @xmath71  pc  @xmath37 , and incidental low - intensity pulses whose s / n peaks at dm = 0  pc  @xmath37 tend to be smeared below the detection threshold for dedispersed time series with trial @xmath72  pc  @xmath37 .
a more complex signature is observed for federal aviation administration ( san juan airport ) radar pulses , which are unfortunately common in arecibo observations .
the radar rotation period is @xmath73  s , and each pulse has an envelope that is @xmath74  s wide and consists of sub - pulses with variable period on the order of @xmath75  ms . depending on telescope orientation ,
radar pulses may be detected in all , some , or none of the alfa beams due to reflections off the telescope structure .
radar pulses are up to two orders of magnitude brighter than pulsar pulses and , without mitigation , can completely dominate single pulse search results .
in addition , the modulation of the radar signal is manifested as detections with @xmath76  pc  @xmath37 so that unlike other non - radar rfi , their s /
n vs. dm signature can not be used for excision .
we exploit the known radar characteristics as well as the pattern of pulse detection in the 7 alfa beams to excise both radar and non - radar rfi .
the first part of our excision algorithm targets radar pulses . after a list of single pulse events
is generated for all trial dms , we bin the events for trial @xmath77  pc  @xmath37 in time ( by 0.1  s ) and record the number of events in each time bin .
then we treat the histogram as a time series and perform a discrete fourier transform .
a peak near the radar rotation frequency indicates a significant number of radar pulses in the data . from the fourier components we extract the phase and
find the arrival time of the earliest radar pulse , @xmath78 .
since the envelope width is @xmath74  s , events within 0.5  s of that location are excised , and the procedure is repeated for events near @xmath79 , where @xmath80 is an integer .
the more radar pulses are present , the better the performance of this technique because the radar peak in the dft is more prominent and the pulses arrival times are determined more precisely .
however , if only a couple of strong radar pulses are present within the typical 268  s palfa integration time , they may be bright enough to dominate event statistics , yet the dft method does not excise non - periodic incidental rfi .
consequently , after applying the dft - based method we use an additional rfi filter that handles aperiodic cases .
the second filter uses the number and proximity of beams in which an event is detected in order to determine if it is due to rfi . again
events for trial @xmath77  pc  @xmath37 from each beam are binned in time . after detecting peaks indicating an excess of events for
the respective time bins , the algorithm cross - checks between results for all seven beams , and each event falling within a histogram peak receives a penalty grade based on how many beams histograms exhibit a peak and how close to each other they are on the sky .
most pulsars detected blindly via a single pulse search appear in one beam or two adjacent beams , and very bright pulsars may be detected in three or four adjacent beams .
we set the excision penalty threshold just below the value corresponding to the latter configurations and excise events accordingly .
this method complements the dft cleaning scheme and effectively excises sparse radar blasts as well as non - periodic rfi detected in multiple beams .
the application of the two excision algorithms makes a marked difference in the final single pulse search output for pointings contaminated with rfi ( fig .
[ spcleaned ] ) .
the figure shows some false positives , for example a clump of pulsar pulses in beam 4 around @xmath81  s are excised due to low - level rfi at low dms occurring in non - adjacent beams in the same 0.1  s time bins as the pulses .
lowering the threshold according to which an excess of events is defined and time bins are marked for excision reduces false positives but also diminishes the effectiveness of rfi excision . while tens of pulses are detected within 134  s from the bright pulsar b2020 + 28 ( fig .
[ spcleaned ] ) , for sources discovered via a single pulse search the number is an order of magnitude smaller ( see   7 ) .
therefore the chance of rfi occurring simultaneously with pulses of intermittent sources and causing them to be excised is much lower .     vs. dm channel and time , number of events vs. dm channel , and event s / n vs dm channel .
the pulsar detection in beams 3 and 4 is evident after rfi events are excised . ]     vs. dm channel and time , number of events vs. dm channel , and event s / n vs dm channel .
the pulsar detection in beams 3 and 4 is evident after rfi events are excised . ]
the palfa survey has made a total of 354 blind detections of 172 pulsars up to late 2008 . following quick - look processing of the data with degraded time and frequency resolutions at the time that data are acquired ( @xcite ) ,
data are shipped to external sites for processing through the cornell pulsar search code and presto .
table  [ detection_stats ] shows a breakdown of cornell detections , including pulsars originally discovered by the quick - look pipeline , by the periodicity ( fft ) and single pulse ( sp ) search algorithms for known pulsars and new discoveries , and table  [ psr_params ] lists parameters for palfa single pulse discoveries that are discussed in more detail in   [ section_objects ] .
columns in table  [ psr_params ] are ( 1 ) pulsar name ; ( 2)-(3 ) : equatorial coordinates ; ( 4 ) period ; ( 5 ) pulse width ( fwhm ) ; ( 6 ) dispersion measure ; ( 7 ) distance ; ( 8) peak flux density of the brightest detected pulse ; ( 9 ) total observation time ; ( 10 ) total number of pulses detected ; and ( 11 ) implied average pulse rate . the last column ( 12 ) points to the paper section that discusses the object .
ccccc known & 48 ( 38% ) & 5 ( 4% ) & 73 ( 58% ) & 126 + new & 28 ( 61% ) & 6 ( 13% ) & 12 ( 26% ) & 46    lllcccccrrrc j0627 + 16 & 06:27:13(7 ) & 16:12(2 ) & 2.180 & 0.3 & 113 & 3.2 & 150 & 7454 & 48 & 23 &   [ 0627sec ] + j0628 + 09 & 06:28:33(7 ) & 09:09(2 ) & 1.241 & 10 & 88 & 2.5 & 85 & 1072 & 42 & 141 &   [ 0628sec ] + j1854 + 03 & 18:54:09(7 ) & 03:04(2 ) & 4.559 & 50 & 216 & 5.5 & 14 & 388 & 9 & 84 &   [ 1854sec ] + j1909 + 06 & 19:09:24(7 ) & 06:40(2 ) & 0.741 & 1.5 & 35 & 2.2 & 82 & 536 & 10 & 67 &   [ 1909sec ] + j1919 + 17 & 19:19:47(7 ) & 17:44(2 ) & 2.081 & 100 & 148 & 5.3 & 12 & 393 & 35 & 320 &   [ 1919sec ] + j1928 + 15 & 19:28:20(7 ) & 15:13(2 ) & 0.403 & 5 & 242 & 7.4 & 180 & 2880 & 3 & 4 &   [ 1928sec ] +
j1946 + 24 & 19:46:00(7 ) & 23:58(2 ) & 4.729 & 4 & 96 & 4.3 & 101 & 268 & 4 & 54 &   [ 1946sec ]
+    more than half of the detected known pulsars and only a quarter of the newly discovered objects were seen by both the periodicity and single pulse search algorithms .
there is also a significant difference in the fraction of pulsars detected only via a single pulse search . in total ,
74% of the new pulsars as opposed to 42% of the known pulsars were detected either only by fft search or only by single pulse search .
the much higher percentage of new pulsars to be detected only by one algorithm means that most of the palfa discoveries are either periodic emitters too weak to show up in single pulse searches or intermittent objects whose emission is too heavily modulated be detected by a periodicity search .
considering that bright pulsars are more likely to be seen by both the fft and sp search , and they are also more likely to have already been found by previous , less sensitive surveys , this indicates that palfa is probing deeper and finding pulsars that are farther away or less bright than the known objects in the same region of the galactic plane .
single pulse searches have not been routinely carried out on survey data in the past , with the exception of @xcite , @xcite , who discovered psr  j1918 + 08 via a single pulse search , and @xcite , who found 11 rrats .
therefore , the proportion of single - pulse - only detections is higher for the palfa discoveries than the known pulsars .
the five known pulsars detected only by single pulse search are known steady emitters which were seen away from the beam center and therefore with reduced sensitivity .
the sample of palfa single pulse discoveries (   [ section_objects ] ) includes one object which has not been detected again despite reobservations ( j1928 + 15 ) , two sources detectable in re - observations through their time - averaged emission ( psr  j0628 + 09 and psr  j1909 + 06 ) , and four long - period objects , one of which has the smallest known duty cycle for any radio pulsar ( psr  j0627 + 16 ) .    among the 11 parkes rrats ,
10 were successfully redetected in subsequent parkes observations , and one ( j1839@xmath101 ) has not been redetected despite multiple attempts ( @xcite ) .
psr  j0848@xmath143 and psr  j1754@xmath130 can sometimes be detected through their time - averaged emission with sensitive , low - frequency observations ( @xcite ) .
six of the sources
have periods greater than 2  s. in this section we discuss selection effects which may account for some of the observational characteristics of intermittent sources . characteristics of this population which remain unexplained by observational biases may indicate underlying intrinsic differences between these objects and conventional radio pulsars .
they may belong to a genuinely intermittent and physically different class of neutron stars .
we can define a general model for observed signal intensity in terms of the time and frequency - dependent flux density @xmath82 as i(t , ) = g_t(t , ) g_iss(t , ) s_i(t , ) + n(t , ) + rfi(t , ) , where @xmath83 is the telescope gain , @xmath84 is the variation factor due to interstellar scintillation , @xmath85 is radiometer noise , and @xmath86 is the contribution from terrestrial radio frequency interference .    variations in telescope gain across the beams of the multi - beam receiver or within the power pattern of a single beam impact the detectability of a pulsar as either a single pulse or periodic source . as shown by the reobservations of psr  j0627 + 16 , psr  j1909 + 06 , psr  j0848@xmath143 and psr  j1754@xmath130 , classifying a source as intermittent may be due to insufficient sensitivity to detect its periodic emission during the discovery observation .
the effectiveness of single pulse vs. fft - based periodic search depends on the pulse amplitude distribution , and the observed pulse amplitude distribution in turn depends on the pulsar location with respect to the telescope beam power pattern .
[ pulsepdf ] shows the probability density function ( pdf ) of event amplitude for psr  b1933 + 16 , a strong pulsar detected in two beams of the same palfa pointing .
one beam of the alfa receiver is pointing directly at the pulsar and the low - intensity tail of the distribution is clearly visible since all pulses are above the detection threshold .
the pulsar is also detected in the first sidelobe of an adjacent beam , in which case only the brightest pulses are above the detection threshold and only the high - intensity tail of the distribution is visible .
thus an off - axis detection of a steadily emitting pulsar may misrepresent it as an intermittent source .
a different set of effects which affect detectability in single pulse vs. periodic search is due to intervening ionized gas .
turbulence in the ionized interstellar gas and the motion of the pulsar with respect to the gas introduce phase modulations in the pulsar emission which cause the observed intensity to change over time and frequency . in general , the scintillation gain factor has contributions from both diffractive ( diss ) and refractive ( riss ) interstellar scintillation , such that @xmath87 @xcite .
the detectability of nearby pulsars with low dms can be significantly affected by diss , which can have time scales comparable to typical survey pointing times ( minutes ) and frequency scales @xmath88 . in terms of the observing frequency @xmath23 and the pulsar distance @xmath89 , the diss time scale @xmath90 .
the scintillation bandwidth @xmath91 and scattering timescale @xmath92 are related through @xmath93 . an empirical fit from @xcite for @xmath23 in ghz and @xmath92 in ms gives an estimate for @xmath92 based on dm : _ s = -6.46 + 0.154  dm + 1.07 (  dm)^2 - 3.86  _ f _ [ eqn_taus ] , but there is significant scatter about this relationship . when the scintillation bandwidth and time scale are small , averaging over the observation time and bandwidth can quench diss .    for refractive scintillation , @xmath94 and @xmath95 .
unlike for diss , the time scale of riss increases with distance and therefore it can become very long for distant pulsars with relatively high dms . both riss and diss can either help or hinder pulsar detection due to the time - variable constructive or destructive interference of wave fronts they cause . for a single - pass survey
, this means that some weak pulsars may be detected because they are modulated above the detection threshold , and some pulsars that are otherwise bright enough may be missed because of scintillation modulation .
@xcite discuss in detail the detection probability for scintillating pulsars by single - pass and multi - pass surveys .
since palfa is a single - pass survey , it may have missed pulsars that scintillate on time scales comparable to or longer than the observation time .
none of our discoveries in table  [ psr_params ] appear to be affected by diss , though riss over long time scales may affect the objects with higher dms .
the detectability of a pulsar as a periodic source depends on the integrated pulse flux within the observation , which is determined by the number of periods within the observation time , @xmath96 , along with the pulse amplitude distribution . for the same integrated flux per pulse ,
fewer periods within a fixed @xmath97 mean that long - period pulsars are less likely to be detected than short - period pulsars by any method relying on detecting integrated flux such as an fft - based search or a continuum imaging survey .
for fft searches in particular , this effect is compounded by the fact that when @xmath98 is small , the harmonics of interest are in the low - frequency part of the power spectrum , which is often dominated by non - gaussian red noise and rfi .
hereafter we use the term `` small-@xmath98 bias '' to refer to the combined influence of these effects on pulsar detectability .
both the periodicity and single pulse searches perform with varying efficiency depending on an object s degree of intermittency . in this section
we present an intermittency measure method to quantify the relative performance of the two search algorithms , apply it to results from the parkes multibeam and palfa surveys , and discuss general implications for surveys .
we can compare the detectability of objects by periodicity and single pulse search and attempt to narrow down different classes of intermittent emitters by calculating the intermittency ratio r = ( s / n)_sp/(s / n)_fft for each object in our sample of palfa and parkes detections . @xcite
derive r = ( ) [ eqn_r ] , where @xmath99 and @xmath100 for a gaussian pulse shape , @xmath101 is the maximum expected pulse intensity within @xmath102 periods and @xmath103 is a modified average pulse peak intensity .
the intermittency ratio is a system - independent measure of the efficiency of the two types of search for objects with different degrees of intermittency . in fig .
[ r_vs_np ] , expected values of @xmath104 vs. @xmath105 are shown for an exponential pulse amplitude pdf and power - law distributions with various indices .
lower limits on @xmath104 are given for intermittent sources which were not detected via periodicity search .
in addition , we show @xmath104 values for pulsars that were detected using both the periodicity and single - pulse search algorithms in the palfa and parkes multibeam surveys . a total of 283 palfa pulsar detections and 255 parkes multibeam detections are shown in the plot , including multiple detections of some objects .
vs @xmath96 for 283 palfa ( filled hexagons ) and 255 parkes multibeam ( squares , open hexagons ) pulsar detections .
squares denote the 10 parkes rrats with period measurements , including two which can sometimes be detected through their time - averaged emission .
the criterion for inclusion of canonical pulsars was a simultaneous periodic and single pulse detection .
palfa points include both blind survey detections and targeted test observations of known pulsars .
arrows denote lower limits for palfa and parkes intermittent sources .
points with bars show the average @xmath104 for a set of observations of the same object , with the bars denoting the maximum and minimum @xmath104 values .
millisecond pulsar detections appear at lower right .
dotted lines show @xmath104 for power - law pulse amplitude pdfs with indices from 0.5 to 2.0 .
a ratio of @xmath106 was assumed for the cutoff intensities @xmath107 and @xmath108 ( @xcite ) .
the dashed line shows @xmath104 for an exponential pdf . a single pulse search performs better than an fft periodicity search for @xmath109 . ]      in fig .
[ r_vs_np ] , long - period pulsars are predominantly found in the upper left ( small @xmath105 , large @xmath104 ) and millisecond pulsars in the lower right ( large @xmath105 , small @xmath104).psr  b0525 + 21 , with @xmath110 , has @xmath111  s and @xmath112  pc  @xmath37 , which make it susceptible to both the small-@xmath105 bias and diffractive scintillation and therefore it is detected more readily as a single pulse source . at the other extreme is psr  b1937 + 21 , a millisecond pulsar which emits giant pulses but is nevertheless detected with a much higher s / n in the periodicity search than in the single pulse search ( @xmath113 ) because its period is small compared to the palfa observation time and the normal pulse amplitude is exceptionally steady ( @xcite ) .
this indicates that the effect of @xmath105 on detectability can overshadow individual properties , and in survey mode millisecond pulsars will show up overwhelmingly as periodic sources . in that case , processing the data with a single pulse search algorithm can reveal the presence of giant pulses in an otherwise steady emitter .    in - between psr  b0525 + 21 and psr  b1937 + 21 , from upper left to lower right
is a relatively smooth distribution of pulsars whose intermittency ratios are predominantly determined by the influence of @xmath105 on the detectability of all types of periodic emitters .
we look for unusual objects among the outliers from that trend .
there are a number of pulsars detected both in periodicity and single pulse search in the region where @xmath114 , including psr  b0656 + 14 , which occasionally emits bursts much brighter than the average pulse and would arguably be classified as a rrat if located farther away ( @xcite ) .
interspersed among them in the region are most of the palfa and parkes intermittent sources .
psr  j1909 + 06 , which was not detected in a periodicity search but is visible as a periodic source when folded with its known period , has @xmath115 , and the long - period palfa intermittent sources psr  j1854 + 03 and psr  j1946 + 24 have @xmath116 and 1.8 , respectively . of the parkes
intermittent sources , seven have @xmath117 .
psr  j1754@xmath130 ( @xmath118 ) and psr  j0848@xmath143 ( @xmath119 ) were detected through their time - averaged emission in low frequency follow - up observations . in the @xmath120 region
we find only five objects : psr  b0525 + 21 , the crab pulsar , parkes intermittent sources psr  j1317@xmath15759 and psr  j1819@xmath11458 , and palfa intermittent source psr  j1928 + 15 .
the crab pulsar has a high intermittency ratio due to its giant pulses but is otherwise a relatively steady emitter as periodic emission is consistently detected even if the giant pulses are ignored . of the remaining three intermittent sources none
are young pulsars like the crab .
the pulse amplitude distributions of psr  j1317@xmath15759 and psr  j1819@xmath11458 are described by power laws with indices @xmath74 ( @xcite ) , smaller than measured giant pulse indices of @xmath121 ( e.g. @xcite , @xcite ) .
we conclude from fig .
[ r_vs_np ] is that periodicity and single pulse searches should be used in combination by pulsar surveys since there are types of sources which are best detected by either method .
as more intermittent objects are found , placing them in @xmath122 space may help us identify physically distinct populations .
in addition , repeated sky coverage can probe intermittency at time scales of days to months .
table  [ psr_params ] lists parameters of the seven pulsars discovered by palfa via a single pulse search . in this section
we discuss each object s properties in detail and describe the steps taken in addition to the palfa processing pipeline in order to verify the signals celestial origin and calculate the pulsar period from times of arrival of single pulses .
psr j0627 + 16 was discovered when nine single pulses were detected at a trial @xmath123  pc  @xmath37 ( fig .
[ j0627best ] ) .
the period estimated from the pulse arrival times in the discovery observation is @xmath124  s. what is unusual about this object is its very narrow peak in s / n vs. dm space , indicating a narrow pulse ( @xcite ) . fig .
[ 4dynspecs]c shows the cumulative dynamic spectrum after stacking raw data chunks aligned around the five brightest pulses .
a fit to the dispersion sweep in the time - frequency plane gives a best dm of 113  pc  @xmath37 and when the raw data chunk around the brightest pulse is dedispersed with this dm , the fwhm pulse width is only 0.3  ms .
two follow - up observations at 0.33 and 1.4  ghz yielded 17 pulses in 50 minutes and 22 pulses in 72 minutes , respectively . after removing bright single pulses from the dedispersed time series and replacing the respective time samples with a running average , periodic emission with @xmath124  s was detected throughout the 0.33  ghz observation , confirming the period estimate from the discovery observation at 1.4  ghz and the presence of underlying normal emission .
many of the rrats detected by @xcite have short duty cycles , calculated as the ratio between the average width of single pulses and the period .
psr  j0627 + 16 has the smallest known duty cycle , @xmath125 by this definition .
however , for steady emitters the duty cycle is typically defined as the ratio of the folded pulse profile fwhm and the period . individual pulses may occur in a phase window wider than any individual pulse and corresponding to the folded profile width . at 1.4  ghz
, the individual 0.3  ms pulses of j0627 + 16 occur in a 8  ms window , while at 0.33  ghz individual pulses have widths of @xmath126  ms and occur in a 15  ms window .
the folded pulse profile at 0.33  ghz has fwhm of 60  ms , suggesting that bright pulses may be confined to a narrower window than pulses comprising the underlying normal emission .
the detection of bright single pulses and a weak periodic signal from psr  j0627 + 16 when data are folded with an accurate period estimate , along with the absence of a periodic detection in search mode weighs in favor of the proposition of @xcite that some rrats may have periodic emission with highly variable pulse amplitudes .
psr j0628 + 09 was discovered by detecting three pulses at @xmath127  pc  @xmath37 .
the period estimated from the pulse times of arrival was 2.48  s. in follow - up observations the pulsar was also detected in periodicity searches as it emits on average several bursts per minute .
the periodicity detections and the larger number of single pulses in subsequent observations allowed the actual pulsar period of 1.241  s to be determined ( @xcite ) .
psr j1854 + 03 was discovered in 2008 via a single pulse search performed on full - resolution survey data .
four pulses were detected at @xmath128  pc  @xmath37 during the 268  s observation . the period @xmath129  s was estimated by taking the smallest difference between times of arrival ( toas ) of two consecutive pulses and verifying that intervals between all four toas are integer multiples of it . a confirmation observation with the more sensitive central alfa beam yielded within 120  s five pulses whose arrival times match the estimated period .
psr j1909 + 06 was observed in 2006 but not identified as a candidate until 2007 when the full - resolution data were searched for single pulses .
two pulses with a width of @xmath74  ms and signal to noise ratio of 6 and 9 were detected at @xmath130  pc  @xmath37 .
the stacked dynamic spectrum clearly shows a dispersion sweep ( fig .  [ 4dynspecs]d ) .
psr  j1909 + 06 was discovered in data from the off - axis beam 2 of the multi - beam alfa receiver . since the on - axis gain of the center beam is 10.4  k / jy compared to 8.2  k / jy for the other six beams ( @xcite ) , we aimed the center beam at the discovery coordinates for a confirmation observation .
for the same integration time of 268  s , 8 pulses with @xmath131 were detected in the confirmation observation .
we used the pulse arrival times to determine a period and arrived at a best estimate of @xmath132  s.      during the 2007 discovery observation of psr  j1919 + 17 , multiple pulses at @xmath133  pc  @xmath37 were detected in beam 4 of the alfa receiver , with no corresponding periodic detection . a fast folding analysis ( @xcite , @xcite )
was used to narrow down the period to @xmath134  s , which was confirmed in 2008 , when the pulsar was detected as a normal periodic emitter with the more sensitive central alfa beam .
the pulsar has a double - peak profile with the two peaks @xmath5  ms apart , which most likely made it difficult to find its period from pulse arrival times alone .
psr j1928 + 15 was discovered in 2005 by detection of what looked like a single bright pulse at @xmath135  pc  @xmath37 in a 120  s observation ( fig .
[ rratquick ] , top ) .
more detailed analysis revealed that the event was in fact composed of 3 separate pulses occurring at intervals of 0.403  s , with the middle pulse being brighter by an order of magnitude than the other two ( fig .
[ rratquick ] , bottom ) . in fig .
[ 4dynspecs]a the dispersion of the brightest pulse by ionized interstellar gas is shown in the time - frequency plane , evidence of the non - terrestrial origin of the pulses .
a fit to the pulse signal in the time - frequency plane resulted in a refined estimate of @xmath136  pc  @xmath37 . despite multiple follow - up observations ,
the source has not been detected again .    given the dm of this source , 242  pc  @xmath37 , it is unlikely that the non - detection can be attributed to diffractive scintillation .
since the three pulses are equally spaced , they can be interpreted as a single event seen in successive rotations of a neutron star .
this signature might be accounted for by an object that is dormant or not generally beamed toward the earth and whose magnetosphere is perturbed sporadically by accretion of material from an asteroid belt @xcite .
psr j1946 + 24 has @xmath137  pc  @xmath37 and @xmath138  s and was discovered by detecting 4 individual pulses .
the intervals between detected pulses in this case were all @xmath139  s. the brightest detected pulse of psr  j1946 + 24 has @xmath140 and its dispersion sweep is visible in the time - frequency plane without any stacking ( fig .  [ 4dynspecs]b ) .
single pulse arrival times collected over multiple follow - up observations were used to obtain partial timing solutions for psr  j0628 + 09 , psr  j1909 + 06 , psr  j1919 + 17 , psr  j1854 + 03 and psr  j1946 + 24 and verify their estimated periods .
the periods of psr  j0627 + 16 and psr  j1909 + 06 were verified by manually folding raw data and detecting a periodic signal .
finally , despite the fact that psr  j1928 + 15 was not detected after the discovery observation , the @xmath141  s intervals between the three pulses detected from that source differ by only @xmath50  ms , which is approximately half the fwhm pulse width and therefore provides a period estimate to that precision .
we are currently timing the six sources which were successfully redetected .
obtaining phase - connected timing solutions will enable us to compare them to other neutron star populations and will yield positions which will facilitate observations at higher wavelengths .
the palfa survey was designed to detect pulsars and pulsar - like objects by identifying either ( or both of ) their periodic or single - pulse emission , as discussed earlier in this paper .
the same data may also be used to detect transient events from other galactic and extragalactic sources . at minimum ,
the palfa survey can place limits on the rate and amplitude distribution of radio transients that are shorter than about 1 second in duration .
possible source classes for short - duration transient events , to name a few , include flare stars , magnetar bursts , gamma - ray bursts ( grbs ) , and evaporating black holes .
grbs have diverse time scales spanning @xmath41  ms to 1000  s. mergers of double neutron - star ( dns ) binaries and neutron star - black hole ( ns - bh ) binaries are proposed sources for short bursts , and they may also emit contemporaneous radio pulses ( @xcite , @xcite ) .
merger rates in the galaxy are estimated to be @xmath74 to 150  myr@xmath142 for dns binaries from pulsar surveys and a factor of 20 - 30 less for ns - bh binaries from population synthesis studies ( e.g. , @xcite ) .
event rates of detectable events depend critically on the luminosity of radio bursts but are likely to be comparable to the rate of short - duraction grbs , and thus smaller than @xmath74  event  day@xmath142  hemisphere@xmath142 .
another phenomenon that may produce rare , bright pulses is the annihilation of primordial black holes ( @xcite ) .
@xcite estimate the galactic rate for such events to be @xmath143  kpc@xmath144  yr@xmath142 .
other possibilities are discussed in @xcite .    motivated by these considerations , we discuss the constraints that can be made on transient events from palfa survey data obtained to date .
first , we note that most  but not all  of the aperiodic events detected to date in palfa data ( as reported in   [ section_objects ] ) and in the parkes multibeam survey by @xcite , appear to be periodic when they are reobserved sufficiently to establish a periodicity .
additionally , _ all _ events detected so far in both surveys have dispersion measures that can be accounted for by ionized gas in the galaxy ( using the ne2001 model ) , suggesting that the emitting sources are closely related to the standard galactic pulsar population . by constrast
, @xcite recently detected a strong ( 30 jy ) , isolated burst with duration @xmath145  ms that shows dm @xmath146  pc  @xmath37 , too large to be accounted for by modeled foreground material in the galaxy or by plasma in the small magellanic cloud , which is somewhat near the direction the event was found .
in addition to detecting radio sources on axis in data streams from each of the seven beams of the alfa system , it is possible to detect strong events off axis with low but non - zero sensitivity over almost the entire hemisphere centered on the zenith .
consequently , it is possible to detect relatively low - rate but very bright transients off axis as well as the much weaker but higher rate single pulses that we have detected in the on - axis parts of the beam pattern .    to establish that a given event is celestial in nature , we need to distinguish it from terrestrial interference . for on - axis events ,
this classification process is aided by the directionality of the seven beams and the expectation that most real events are likely to be seen in three or fewer contiguous beams rather than in all seven beams .
bright events may violate this expectation , however , because near - in sidelobes are fairly large ( c.f .
figure 1 of @xcite ) .
celestial and terrestrial radio signals can also enter the seven beams of the alfa receiver indirectly  far off the nominal pointing axis  through reflection and scattering off telescope support structures , as with any antenna .
extremely bright signals ( celestial or terrestrial ) may be detected after scattering into one or more beams , and the effective gain is approximately that of an isotropic radiator with small gain for events incident far off axis from the pointing direction .
for the arecibo telescope , which has an intricate support structure with many possible scattering surfaces , such wide - angle incidence is likely to provide multiple paths for radiation to enter the feed optics .
radiation arriving along multiple paths will constructively and destructively interfere across alfa s focal plane , so that the event s strength will be nonuniform in the different data streams .    when radiation enters the telescope optics from wide angles , it is difficult to distinguish celestial from terrestrial signals based on the nominal pointing direction of the telescope and on whether they are detected in a single beam , a subset of beams , or all beams
. conceivably , extremely bright events from a single radio source could be emitted at very rare intervals that would be detected when the telescope is pointed nominally at directions separated by tens of degrees ( but within the same hemisphere ) .
such sources and events will be exceedingly difficult to disentangle from terrestrial rfi .
of course the most powerful confirmation of a particular event is the detection of similar events in reobservations of the same sky position . for most of the events
reported in   [ section_objects ] and by @xcite ,
redetections have been made in this manner .
they have also relied on establishing that their signatures in the frequency - time plane conform to what is expected .
the simplest celestial signals are narrow pulses that show differential arrival times in accord with the cold - plasma dispersion law ( e.g. @xcite ) .
celestial objects may also show the effects of multipath propagation through the interstellar medium ( ism ) either in the frequency structure that is produced in the spectra of nearby pulsars or in the asymmetric broadening of the pulses from more distant objects . for cases where spatial information associated with the nominal telescope pointing direction is insufficient to constrain the celestial nature of a particular event
, we must rely especially on a detailed study of the event s structure in the time - frequency plane .
indeed the time - frequency structure of the event discussed by @xcite was a key part of the argument for its classification as a celestial event .
a given reflector and feed antenna have a net antenna power pattern with a main lobe and near - in sidelobes that is similar to an airy function with a main - lobe angular width @xmath147 , where @xmath148 is the effective reflector diameter .
the wide - angle part of the antenna pattern has far - out sidelobes that are independent of @xmath148 and have average gains @xmath149 , similar to that of an isotropic antenna .
however peaks in the far - out sidelobe pattern can have gains significantly larger than @xmath150 , corresponding to special directions where scattering from the telescope support structure is especially efficient .
for the arecibo telescope , we expect sidelobes to have the same general properties , with the caveat that support structures will introduce further complexity in the far - out sidelobe pattern that varies as the telescope is used to track a source .
tracking a sky position is done by rotating the azimuth arm and moving the gregorian dome along the azimuth arm , thus changing the scattering geometry .
we therefore expect the sidelobe structure to change significantly with azimuth and elevation of the targeted position on the sky .
for the seven - beam alfa system , the main beams by design sample different parts of the sky , as do the near - in sidelobes . at wider angles , the power pattern of each feed overlaps the others but with considerable variation in gain from feed to feed .    to derive simple constraints on burst amplitudes and rates , we expand the gain for each feed into inner beam and wide - angle terms : g ( , ) = g_max p_n ( , ) + ( 1 - _ b ) , where @xmath151 and @xmath152 are polar and azimuthal angles , respectively ; @xmath153 is the power pattern for the main beam and near - in sidelobes ; @xmath154 is the main - beam efficiency , and @xmath155 is the average level of the far - out sidelobes .
the boresight gain is g_max = = = , [ eqn_gmax ] where @xmath156 is the effective telescope area , and @xmath157 and @xmath158 are the solid angles of the antenna power pattern and of the main beam , respectively . the system equivalent flux density is s_sys = = , [ eqn_ssys ] , where @xmath159 is the boltzmann constant and @xmath160 is the system temperature .
the minimum detectable flux density is s_min = , [ eqn_smin ] which is a function of the event duration , @xmath19 , as well as system parameters .
it is useful to define @xmath161 as the minimum detectable flux density for unit gain using eq .
[ eqn_ssys ] with @xmath162 . substituting into eq .
[ eqn_smin ] yields , for two polarizations and @xmath163  mhz , s_min , 1 = 10 ^ 4.9  jy   ( ) ( ) ^1/2 ( ) ; we note that for some directions , the system temperature can differ from the fiducial 30  k used .
we now derive limits that can be placed on the rate and amplitudes of transient events using our simple model for the antenna power pattern .
let @xmath164 be the event rate per unit solid angle from a population of sources , which implies a total number @xmath165 in a total observation time @xmath97 .
if no events are detected , we can place an upper limit on the rate @xmath166 for minimum flux density , @xmath167 .
we treat the response of the antenna power pattern as azimuthually symmetric about the bore axis , so @xmath168 , and calculate the gain as a function of polar angle from the bore axis . for each annulus of solid angle @xmath169 ,
the upper bound on the rate is _ s , max , for flux densities greater than s_min = . the multiplier @xmath170 accounts for whether a given alfa feed and receiver combination samples the same patch of sky or not .
for the main lobes , @xmath171 while for the far - out sidelobes , @xmath172 .
we evaluate fiducial values by integrating respectively over the main beams and over the remainder of the power patterns . for a main - lobe beam width of 3.5 arc min , @xmath173  deg@xmath13 and @xmath174 . using the average gain over the main beam out to the half - power point , @xmath175
, we have s_min = 11  mjy   ( ) ( ) ^1/2 ( ) and _ s , max = 0.12  hr^-1  deg^-2 ( ) . integrating over wide angles , and using @xmath176 , we estimate s_min 10 ^ 5.4  jy   ( ) ( ) ^1/2 ( ) , where the far - out sidelobes have an average gain of @xmath177 , and _ s ,
max =  hr^-1  deg^-2 ( ) . where @xmath178 is the fraction of the hemisphere actually covered by the far - out sidelobes .
results for palfa data are shown in fig .
[ farout_palfa ] . within the half - power beam width
the constraints are strong on the flux density but weak on the event rate while the opposite is true for the far - out sidelobes , represented by the part of the curve below @xmath179  h@xmath142  deg@xmath180 .
the asympotic value of the curve corresponding to the maximum , boresight gain differs from the fiducial value calculated above because the latter is an average over the main beam . in our simple model the near - in sidelobes
are not included explicitly so our curves overestimate @xmath167 at the corresponding polar angles . however , the limits change rapidly , so this does not change the salient features of the figure , which are the seven order - of - magnitude difference in constrained amplitudes and the eight order of magnitudes on the event rates between the main lobs and the far - out sidelobes .
h and an assumed pulse width of @xmath181  ms .
the three points , from left to right , correspond to the far - out sidelobes , the break between the main beam and far - out sidelobes ( at @xmath182  h@xmath142  deg@xmath180 ) , and the boresight gain .
the shaded region is excluded by the palfa data .
[ farout_palfa ] ]      a brief summary of our results is that ( a ) a small number of events has been detected from sources that are consistent with their membership in the radio pulsar population , but with a much greater deal of modulation in some cases ; ( b ) no very strong pulses have been identified in the far - out sidelobes of the arecibo telescope and ( c ) no bursts have been detected that are extragalactic in origin .
here we compare our results with those of @xcite , who surveyed high - galactic - latitude regions that included the magellanic clouds and who discovered a single reported pulse , with amplitude 30  jy .
they surveyed @xmath183  deg@xmath13 with a detection threshold @xmath184  mjy , about 100 times fainter than the detected pulse amplitude .
to compare any two surveys , we initially assume a population of radio sources that is homogeneously distributed in euclidean space and we ignore propagation effects that might limit the detectability of bursts from distant sources . the number of events detected @xmath185 for a survey of duration @xmath186 that samples instantaneously a solid angle @xmath187 . for multibeam systems with @xmath188 pixels , @xmath189 , where @xmath190 is the solid angle of a single beam .
we assume events have the same peak luminosity and pulse width @xmath19 and are detectable to a distance @xmath15 .
the ratio of number of detected events in two surveys with total observation times @xmath191 and @xmath192 , and bandwidths @xmath193 and @xmath194 is & = & ( ) ( ) ( ) ^3 + & = & ( ) ( ) ( ) ^3/2 ( ) ^3/4 .
[ eqn_nbna ]    as pointed out by @xcite , the same survey that detected the 30-jy pulse should have yielded many more detections if the assumptions of a volume limited sample apply . using eq .
[ eqn_nbna ] , the lone 30  jy pulse implies there should have been @xmath195 additional detections above a 100-times fainter threshold . in the 90  h of follow up observations ,
there should have been @xmath196 additional detections .
comparing the parkes and palfa surveys , a system - equivalent flux density of 4  jy , and 100 mhz bandwidth .
for the parkes survey we use 13 beams of 14  arcmin diameter , 480  h , a threshold of 600 , @xmath197  jy , and 288  mhz bandwidth . ] , we estimate that the palfa survey should have identified @xmath198 pulses above threshold .
the palfa survey s null result on extragalactic events is therefore in accord with lorimer et al.s non - detection of events weaker than the single strong event they reported .
with the assumptions made , these results suggest that the non - detections of pulses weaker than the 30-jy pulse are inconsistent with the source population having a distribution that is homogeneous and isotropic .
there are caveats on these results related to the assumptions about the source population and to the possible role of propagation effects in limiting detectability .
the minimum detectable flux density for the palfa survey is 30  mjy for a 5-ms pulse , implying that a pulse of 30  jy strength emitted at a distance @xmath199 is detectable to a distance @xmath15 given by @xmath200 .
if the lone event detected by lorimer et al . is from a source at the attributed distance @xmath201  gpc ( based on assigning about half of the dm to the intergalactic medium ) , the palfa survey would detect fewer pulses than estimated above ( for constant luminosity ) because the universe is not old enough and from redshifting of the spectrum if it is steep in frequency .
nonetheless , a cosmological population of constant - luminosity sources would extend to @xmath202 times the attributed distance , implying from the first form of eq .
[ eqn_nbna ] that @xmath203 pulses should have been detected in the palfa survey , on average , while @xmath204 weaker pulses should have been seen in the discovery survey with parkes .
the difference in yield under this alternative scaling results from the fact that the palfa survey can detect sources well beyond the cosmological population while the parkes survey is comparatively shallower ( by a factor of ten in flux density , i.e. 300 compared to 30  mjy ) , but covers a factor of 30 greater solid angle .
consideration of a broad luminosity function does not alter these conclusions qualitatively .
however , we would generally expect the pulse amplitude distribution to extend to _ lower _ flux densities ( from lower luminosities ) given that the sole reported event was much larger than the threshold and that , therefore , many more pulses are expected compared to the constant - luminosity assumption .
pulse broadening from multipath propagation in either the ism or the intergalactic medium ( igm ) can also limit the numbers of events detected in a survey because the broadening increases with source distance .
the maximum detectable distance is therefore diminished , as demonstrated for the ism in figure  [ fig : dmax_vs_lp ] .
lorimer et al .
argued that the strong frequency dependence of the width of the 30-jy pulse was consistent with multipath propagation through a turbulent , ionized igm .
if so , sources from a greater distance would be less - easily detectable . in the simplest scattering geometries ,
scattering conserves the area of the pulse , so matched filtering yields a test statistic with s / n @xmath205 when the pulse broadening time @xmath206 is much larger than the intrinsic pulse width ( @xcite ) .
the 30-jy pulse therefore could have been broadened by an additional factor of @xmath207 and the s / n still would have been above threshold .
the much longer pulses ( @xmath208  ms = 50  s ) would have been difficult if not impossible to identify in the time series owing to high - pass filtering used to mitigate baseline fluctuations .
the high - pass filtering was explicit in the parkes multibeam data acquisition but is part of the post - processing in the palfa analysis .
such very broadened pulses could be detected by imaging surveys with relatively short integration times like the nrao vla sky survey ( nvss ) and future surveys with the australian square kilometer array pathfinder or the allen telescope array .    in between this extreme case of 50-s broadening and the reported pulse width of @xmath145
ms is a substantial search volume in which more distant sources could have produced detectable events . broadening from sources at redshifts greater than one
is increased by the larger electron density but is lessened by the frequency redshift , so the net scaling of the pulse broadening with redshift depends on how the scattering material is distributed along the line of sight .
we have examined the current state of knowledge about various classes of intermittent radio - loud neutron stars , summarized selection effects that may affect the classification of intermittent sources , and presented our data processing methods and results from single pulse searches performed on palfa survey data .
one of our single - pulse discoveries is a relatively persistent emitter that , like many other pulsars , shows a broad amplitude distribution , from which a few bright pulses were detected .
two sources were detected as normal periodic emitters when observed at a lower frequency ( psr  j0627 + 16 ) and with the more sensitive central beam of the alfa receiver ( psr  j1919 + 17 ) .
four objects were most likely not detected via periodicity search because of their long periods ( @xmath209  s ) compared to the palfa observation time of @xmath210  s , and one intermittent object ( psr  j1928 + 15 ) was discovered by detecting three pulses emitted on successive rotations but it was not detected again despite multiple reobservations .
most of the palfa sources that were discovered solely through the single - pulse analysis have subsequently been found to be _ sporadic _ but with measureable and consistent rates of detectable pulses .
the parkes survey had similar results . however , two sources ( the parkes source psr j1839@xmath101 and the palfa source psr j1928 + 15 ) are much more irregular even though an underlying periodicity has been found in each case .
we therefore distinguish between cases where the apparent sporadic behavior is due to a detection threshold that selects only the strongest pulses from a broad distribution ( e.g. @xcite ) and those where emission is truly _
intermittent_.    a major issue is related to the question , if we see transients from a given sky direction only once , how can we figure out what phenomena produce them .
we must archive radio survey data and characterize the wide variety of detected radio transients outside the framework of rotating neutron stars if necessary .
archived data may yield new discoveries if reprocessed with improved search algorithms and a database of transient signals would be an invaluable reference as theoretical understanding of energetic radio events and prediction of their signatures improves .
the palfa survey is taking steps in both directions : raw data are stored in a tape archive and data processing products are indexed in a database at the cornell center for advanced computing and currently work is under way to make them publicly available via a web portal .
we thank the staff at naic and atnf for developing the alfa receiver and the associated backend systems . in particular , we thank jeff hagen , bill sisk and steve torchinsky at naic and graham carrad at the atnf .
we are also grateful to the parkes multibeam survey team for providing the results incorporated in fig .
[ r_vs_np ] .
this work was supported by the nsf through a cooperative agreement with cornell university to operate the arecibo observatory .
nsf also supported this research through grants ast-02 - 05853 and ast-05 - 07376 ( columbia university ) , ast-02 - 06035 and ast-05 - 07747 ( cornell university ) and ast-06 - 47820 ( bryn mawr college ) .
, d.r.l . and p.c.c.f .
are supported by a wvepscor research challenge grant .
is an alfred p. sloan research fellow .
f.  crawford is supported by grants from research corporation and the mount cuba astronomical foundation .
j.w.t.h . is a nwo veni fellow .
pulsar research at ubc is supported by nserc and the canada foundation for innovation .
i.h.s . acknowledges support from the atnf distinguished visitor program and the swinburne university of technology visiting distinguished researcher scheme .
l.e.k . held an nserc canada graduate scholarship while most of this work was performed .
v.m.k . holds a canada research chair and the lorne trottier chair , and is supported by nserc , fqrnt , cifar , and by the canada foundation for innovation .
the parkes telescope is part of the australia telescope , which is funded by the commonwealth government for operation as a national facility managed by csiro .
nrao is a facility of the nsf operated under cooperative agreement by associated universities , inc .
research in radio astronomy at the nrl is supported by the office of naval research .
mclaughlin , m. a. et al .
2006 , nature , 439 , 817 mclaughlin , m. a. , rea , n. , gaensler , b. m. , chatterjee , s. , camilo , f. , kramer , m. , lorimer , d. r. , lyne , a. g. , israel , g. l. , possenti , a. 2007 , apj , 670 , 1307