arXiv:astro-ph/0701181v1 8 Jan 2007 2 .A ehad,N .R ht .Catre,J .Hn .M G H Bac M. Faucher-Gigu`e C. T. C.-A. D. B. W. Lommen, and Han, J. N. DAmico L. McLaughlin, A. N. J. Crawford, Possenti, 1 F. Chatterjee, A. A. Kaspi, Hobbs, S. M. B. Bhat, V. M. G. Lazio, R. Manchester, W. Faulkn D. J. J. N. T. N. A. Deshpande, Ramachandran, A. R. R. A. Arzoumanian, Leeuwen, van J. Z. Champion, J. D. et,alwn st efr opeesv essof census population. comprehensive a instru- Galactic perform current the to with us detectable allowing detecting not ments, of capability the of have will thousands and op- become 2014 to in the is with erational SKA survey The (SKA). pulsar Array projected Kilometer Square a of impact scientific the remotely. analysis of data users allow perform will to which archive developed the being the are at services tools hosted Web and facility. via database supercomputing community Center and Theory raw the archive Cornell of to tape (PB) available Petabyte sophisticated made 1 a be of will total which a data accumulate stars: to neutron going radio-loud survey PALFA is of The Transients). class RAdio defined (Rotating 3.98 RRATs newly a of character- a in bursts period of sporadic pulsar istic orbital exhibits 144ms an object young Another with a hours. 24 binary is found which relativistic has of highly one processing PALFA pulsars, data the new ex- preliminary of present is far discoveries We and recent So pulsars. survey. years and new 3-5 progress, of next goals, hundreds the the discover for ALFA to continue new to pected the slated commissioned. when is was receiver It 2004 Array) in Feed L-band began (Arecibo plane Galactic the Abstract. .S Deneva S. PALFA Survey J. Arecibo SKA The Projected Future: and and Survey Present Surveys Pulsar rcb eecp.I ea n20 n ssae ocon- to discovering slated is years, and 3-5 the 2004 at for in receiver tinue began ALFA It 7-beam by telescope. new enabled Arecibo the was of plane installation Galactic the the of survey PALFA The parameters Survey 1.1. Survey Pulsar-ALFA Arecibo The 1. h uvycvr aatclatitudes Galactic covers survey The evl flniueacsil oteAeiotelescope, Arecibo the to accessible 32 longitude of tervals hsketu a onf emn,May. Germany, Honnef, on: Bad Seminar Physikzentrum WE-Heraeus 363. the of Proceedings onl nvriy 1 pc cecs taa Y183 U 14853, NY Ithaca, Sciences, Space 516 University, Cornell http://www.naic.edu/alfa/pulsar ◦ eas ics aaees xetddsoeisand discoveries expected parameters, discuss also We ≤ l ≤ 77 h rcb usrAF PLA uvyof survey (PALFA) Pulsar-ALFA Arecibo The ◦ 1 n 168 and .M ods .R oie,P .C rie .Cml,I .S H. I. Camilo, F. Freire, C. C. P. Lorimer, R. D. Cordes, M. J. , ◦ ≤ l ≤ 214 ∼ ◦ 00plas(Fig.1). pulsars 1000 h LAreceiver ALFA The . 14 | b ≤ | − NurnSasadPlas Psesadcnrbtdtalks contributed and (Posters Pulsars” and Stars “Neutron 19 06 d.W ekr ..Hag P eot21 p 2 55 52- pp. 291, Report MPE Huang, H.H. Becker, W. eds. 2006, , 5 ◦ ntein- the in SA, eta em h afpwrba imtra .2GHz a 1.42 around at diameter configuration is beam hexagonal half-power The a beam. in central beams 6 features Galaxy the in density 2002). electron Lazio the & of (Cordes on model based estimated NE2001 were to the Distances accessible plane telescope. Galactic Arecibo the the of regions delimit PALFA Slant- lines discoveries. ing simulated PALFA actual denote denote stars triangles discoveries, pulsars, known note 1. Fig. h eta emand beam central the oe h k ihgi qa oo xedn h half- closely the be exceeding must or pointings to 7-beam equal three gain gain, with maximum sky the cover plane Galactic the of is out looking temperature system The [email protected] ∼ ∼ rmnts n h nai ani 04KJy K 10.4 is gain on-axis the and arcminutes, 3 4K(els20,Cre ta.20) nodrto order In 2005). al. et Cordes 2004, (Heiles K 24 usrlctosi h aatcpae rse de- Crosses plane. Galactic the in locations Pulsar re r .G ye .M Ransom, M. S. Lyne, G. A. er, e,M rmr .W Stappers, W. B. Kramer, M. ker, ese,L ain .Reid, B. Kasian, L. aensler, ∼ 8 ses .Vlemmings, W. essels, . Jy K 2 ) ar,D .Nice, J. D. tairs, − 1 o h te beams. other the for − 1 for Deneva, J. S. et al.: The Arecibo PALFA Survey and Future SKA Survey 53 tiled. Rather than proceed sequentially with tiling the sky, Table 1. PSR J1906+0746 parameters (Lorimer et al. we do sparse sampling, completing one pointing per triple 2005). with the intention to fill in the gaps after all pointings from this first sparse grid are completed. The sparse sam- Parameter Value pling scheme makes use of the large sidelobes of the 6 Right ascension (J2000) 19h06m48.s673(6) outer beams which have ∼ 16% of the on-axis gain, thus ◦ surveying more sky area per unit time (Cordes et al. 2005 Declination (J2000) 07 46’28.6(3)” and references therein). Spin period, P (ms) 144.071929982(3) Spin period derivative, P˙ 2.0280(2) × 10−14 Currently the PALFA survey uses the four Wideband Epoch (MJD) 53590 Arecibo Pulsar Processor (WAPP) backends (Dowd et al. Orbital period, Pb (days) 0.165993045(8) 2000) with a bandwidth of 100 MHz centered on 1.42 GHz, Projected semi-major axis, x (lt s) 1.420198(2) 256 channels, and a sampling time of 64 µs. In early 2007, Orbital eccentricity, e 0.085303(2) we will begin using a new spectrometer capable of pro- Epoch of periastron, T0 (MJD) 53553.9126685(6) cessing the full 300 MHz receiver bandwidth with 1024 Periastron advance rate,ω ˙ (◦ yr−1) 7.57(3) 3 channels. The data is decimated and processed in quasi- Dispersion measure, DM (pc cm− ) 217.780(2) −2 real time during the observations using publicly available Rotation measure, RM (rad m ) +150(10) “quicklook” software (Lorimer 2001). In this pipeline some Flux density at 1.4 GHz, S1.4 (mJy) 0.55(15) Main pulse widths at 50% and 10% (ms) 0.6 and 1.7 resolution is traded for speed, and the data is decimated by P Characteristic age, τc = (kyr) 112 a factor of 8 in frequency and 16 in time. This means that 2P˙ Inferred distance, d (kpc) ∼ 5.4 initial data processing is sensitive mainly to non-recycled Spectral index, α -1.3(2) pulsars with periods greater than ∼ 30 ms. Decimated Mass function, f (M⊙) 0.1116222(6) datasets have a time resolution of 1024 µs and are dedis- Total system mass, Mtot (M⊙) 2.61(2) persed with 96 trial DM values. The resulting time se- Grav. wave coalescence time, τg (Myr) ∼ 300 ries are then searched for periodic signals using FFTs and harmonic summing, and for single pulses using a matched filtering algorithm (Cordes & McLaughlin 2003). Process- ing full-resolution data will be performed on dedicated clusters at the home institutions of PALFA consortium members and will include acceleration searches, which will greatly increase sensitivity to millisecond pulsars and pul- searches for pulsed signals from the companion have not sars in short-period binaries. found any plausible candidates. From these results, we can place an upper limit of 0.1 mJy kpc2 on the luminosity of the hypothetical pulsar companion. Another possibility is 1.2. Results that the companion may be a pulsar which is not beamed As full-resolution processing pipelines are still being towards us. J1906+0746 shows a pulse profile variation tested, the 28 pulsars discovered so far were found in deci- between 1998 and 2005 (Fig.2) which can be explained mated data only. The highlight among these new objects is by geodetic precession. With time, the precessing beam J1906+0746, a young, relativistic binary pulsar with a pe- of the pulsar changes its orientation with respect to our riod of 144 ms and DM of 217.78 pc cm−3 (Table.1). The line of sight and that causes the observed pulse shape to discovery observation in September 2004 detected it with evolve. In this case, an interpulse is not visible in the 1998 a S/N ∼ 11. The short integration time of the PALFA observation, but appears in the 2005 data. pointing (134 seconds) did not allow for its binary na- Among the 28 pulsars discovered by PALFA so far, ture to be immediately seen, but subsequent processing two objects fit the characteristics of a newly defined class of an archive 35-minute Parkes observation of the same of radio-loud neutron stars, Rotating RAdio Transients region of the sky revealed that the 144 ms signal was ac- (RRATs, McLaughlin et al. 2006). J0628+09 was found celerated and follow-up observations revealed its orbital only by the single pulse search algorithm; it exhibits spo- period to be 3.98 hours (Lorimer et al. 2005). Without radic bursts with a peak S/N of up to 100. Only 3 pulses accounting for acceleration, the pulsar was below the de- were found in the 67 seconds of data of the discovery obser- tection threshold of 8 − 9σ for the Parkes Multibeam Sur- vation (Fig.3), and the period was estimated to be 2.48 s. vey and was rejected as a candidate because of known Longer follow-up observations allowed the pulsar to also radio interference close to its period. J1906+0746 has a be detected with the periodicity search algorithm and nar- characteristic age of 112 kyr, and the total system mass row down its period to 1.24 s. J1928+15 was also discov- is 2.61 ± 0.02 M⊙, as indicated by the measured rate of ered only by its single pulses. In that case, only 2 bursts advance of periastron. If we assume the pulsar mass to be were detected during the 134 second observation, which close to 1.4 M⊙, this leaves leeway for the companion to was not enough to constrain the period. Follow-up obser- be either a or a massive white dwarf. Deep vations of this object are still forthcoming. 54 Deneva, J. S. et al.: The Arecibo PALFA Survey and Future SKA Survey

1.3. Data Archive All raw data and data processing products from the PALFA survey will be archived permanently at a state- of-the-art tape library at Cornell Theory Center. This archive is expected to grow to ∼ 1 PB by the time the survey is completed, and presents a challenge both with respect to database management and insuring data in- tegrity and developing appropriate tools to make such a large archive publicly available online. Currently effort is under way to develop web services, web applications and a comprehensive online portal to the archive which will provide PALFA consortium members and external users with a means to request, download, and remotely manip- ulate raw data and data processing products. Through Fig. 2. J1906+0746 pulse profiles at 1.374 GHz from a standardized web services and protocols, the PALFA sur- 1998 (top) and 2005 (middle) Parkes observation. The vey portal will also be interfaced with and made available bottom panel shows the difference of the two profiles after through the National Virtual Observatory. scaling them for the area of the main pulse. The dashed horizontal lines delimit ±3σ calculated from the pulse pro- 2. Future SKA Pulsar Survey file baseline noise (Lorimer et al.2005). With its unprecedented collecting area and sensitivity, Sep 6 09:42:50 2004 the Square Kilometer Array (SKA) promises to revolu- tionize several areas of radio astronomy, including pul- sar astrophysics. The frequency range of the SKA will be 0.1-25 GHz, with a bandwidth equal to 25% of the ob- serving frequency, to a maximum of 4 GHz bandwidth for frequencies above 16 GHz. It will be possible to ob- serve simultaneously at two independent frequency bands with the same field-of-view center and 2 polarizations per band. 75% of the collecting area will be within a 150 km diameter and the maximum baseline length will be up to 3000 km, allowing for an angular resolution of 0.1 arcseconds at 1.4 GHz (for full specification, see http://www.skatelescope.org). The sensitivity of the SKA makes it especially attrac- tive for pulsar searching. For a pulsar with a period P and pulse width W the minimum detectable flux density is

mσTsys W Smin = , (1) G nptint∆ν P − W  p where m is the number of σ corresponding to the detec- tion threshold, Tsys is the system temperature, G is the gain, np is the number of polarizations summed, tint is the integration time and ∆ν is the bandwidth. For the Fig. 3. Discovery observation of RRAT J0628+09. Panels SKA, the factor Tsys/G is ∼ 10 times smaller than for the on top show total number of events vs. S/N (left), num- and ∼ 100 times smaller than for the ber of events vs. DM channel (middle), and DM channel . Assuming m = 8, Tsys = 25 K, vs. S/N (right). The panel on the bottom shows the loca- tint = 60 s and ∆ν = 0.5ν gives Smin ≈ 1.4µJy (Kramer tion of events with S/N ≥ 5 in the DM-time plane. The 2003). This corresponds to a minimum detectable lumi- two dispersed pulses occur at approximately t = 42 and nosity of 0.1 mJy kpc2 for the Galactic Center region and t = 49 s, and have a distinct signature from radio inter- 0.8 mJy kpc2 for the farther side of our Galaxy. A simu- ference, which spans an unrealistically large range of trial lated SKA pulsar survey assuming an integration time of DM values (e.g. at t = 4 and t = 62 s). 10 minutes and all-sky coverage yields ∼ 10000 new pul- sar discoveries (Cordes, priv.comm.), amounting to a com- plete Galactic pulsar census. Among the most eagerly an- Deneva, J. S. et al.: The Arecibo PALFA Survey and Future SKA Survey 55 ticipated discoveries are pulsars closely orbiting the mas- sive black hole in the Galactic Center. The high stellar density in this region means that similar to globular clus- ters, it will be a prime target for searching for binary and milliseconds pulsars, and through timing of these objects as they move through the potential of the black hole, we can hope to measure the black hole spin. The large num- ber of projected SKA pulsar discoveries means that we will have probes of electron number density along many lines of sight through the Galaxy and will be able to refine models of ionized gas structure like the NE2001 model of Cordes & Lazio (2002). In a similar manner, with the SKA we will be able to detect the brightest pulsars in nearby galaxies and use them as probes of the ionized intergalac- tic medium. Even though pulsars in other galaxies may not be detectable by periodic searches, young Crab-like pulsars can be detected by their giant pulses out to inter- galactic distances. Currently our sample of known pulsars outside our own Galaxy is limited to ∼ 10 objects in the Magellanic Clouds. With projected extragalactic pulsar discoveries by the SKA we can study the pulsar popula- Fig. 4. Comparison of expected sources of gravitational tion statistics in different galaxy types. waves and sensitivity limits for various methods of their Another use of pulsars as probes is for gravitational detection. wave detection. Multiple milliseconds pulsars found by the SKA can serve as the endpoints of arms of a huge gravi- the National Astronomy and Ionospheric Center, and the Aus- tational wave detector: the Pulsar Timing Array (PTA). tralia Telescope National Facility. We are grateful to Cornell The concept for the PTA hinges on the fact that pul- Theory Center staff for their work in designing the PALFA sars are extremely accurate natural clocks. A gravitational database and porting of data processing code. The PALFA raw wave passing in the vicinity of the Earth will stretch and data archive and database work is supported by a National Sci- compress the space through which the pulsar signals have ence Foundation RI grant, a Microsoft E-Science grant, and by to travel in order to reach us. Correlated changes in the the Unisys Corporation. arrival times of pulses from pulsars in different areas of The SKA development project is supported by scientific organizations and observatories in the participating countries the sky will indicate space-time distortion due to gravi- (http://www.skatelescope.org). tational waves. The dimensionless strain quantity hc (f) used to measure the magnitude of the distortion can be related to the observed scatter in arrival times and the References integration time: hc (f) ∼ σTOA/tint. Gravitational waves Cordes, J. M. & Lazio, T. J. W. 2002, astro-ph/0207156 with frequencies on the order of nHz (nanoHertz) would Cordes J. M. et al., 2005, ApJ, 637, 446 be accessible only with the PTA (Fig.4). Cordes, J. M. & McLaughlin, M. A. 2003, ApJ, 596, 1142 Last but not least, the large number of expected SKA Dowd, A., Sisk, W., & Hagen, J. 2000, Astronomical Society pulsar discoveries means that it is statistically likely the of the Pacific Conference Series, 202, 275 sample will contain a corresponding number of exotic Heiles, C. 2004, ALFA Memo, http://alfa.naic.edu/memos objects: relativistic binaries, double pulsar binaries, per- Kramer, M. et al. 2004, New Astronomy Reviews, 48, 993 haps also the long sought-after sub-millisecond pulsars Lorimer, D. R. 2001, Arecibo Technical Memo No. 200101 and pulsar-black hole binaries. The sample of binary pul- Lorimer D. R. et al., 2006, ApJ, 640, 428 sars whose masses can be determined can constrain the McLaughlin, M. A. et al. 2006, Nature, 439, 817 neutron star mass distribution. We do not yet know what are the physical limits on pulsar masses or periods. The large SKA pulsar sample would allow us to refine estimates of the likelihood of the slowest, fastest, least or most mas- sive pulsar likely to exist, and to constrain the equation of state of matter under extreme conditions. Acknowledgements. The PALFA survey is supported by the National Science Foundation through a cooperative agreement with to operate the . The survey has benefited from the expertise of staff at Arecibo,