Pulsar Timing Arrays Current and Future Instrumentation Joseph Lazio Jet Propulsion Laboratory, California Institute of Technology © 2012 California Institute of Technology. Government sponsorship acknowledged. Pulsar Timing Arrays Current and Future Instrumentation • What are the requirements? Timing precision and instrumental performance • What exists? Telescopes Feeds, receivers, and backends IT infrastructure • What is on the horizon? Pulsar Timing Arrays Current Results J1713+0747 30 ns RMS J1012+5307 280 ns RMS J1643-1224 1500 ns RMS Demorest et al. 2012, submitted Timing Precision • Assume σ ≤ 100 ns TOA SNR May be necessary, but not sufficient W P condition W σ = TOA, n SNR • W pulse width • SNR signal-to-noise ratio • Assumes sufficient calibration of telescope system Other contributions can include pulse-phase jitter, uncorrected propagation effects Radio Telescope Radio Telescope Data Data Radio Acquisition Storage Frequency System Intermediate Frequency Baseband Local Local Oscillator Oscillator Radio Telescope Noise Temperature • Replace entire telescope T system by resistor in a heat bath • Output voltage equivalent • Not necessarily physical temperature • Tsys = Tsky + Tspill + TRx + … • P = kBT∆ν Radiometer Equation For TPSR << Tsys T 2k T ∆T = sys ∆S = B sys 2∆ν∆t Aeff 2∆ν∆t Where can improvements be made? • T portion determined by telescope S spectral flux sys density [W/m2/Hz] T = T + T + T + … sys sky spill Rx 1 Jy = 10−26 W/m2/Hz • ∆ν – processed bandwidth • Aeff – effective area of telescope • ∆t – observation time within limits imposed by pulsar or ISM Improvements benefit both timing and survey programs Pulsar Observations Observational frequency scattering determined by balancing Pulsar • pulsar spectrum vs. spectrum •sky spectrum vs. •scattering Typically about 1 GHz sky spectrum Pulsar Timing Arrays Current Instrumentation Pulsar Timing Arrays Current Instrumentation Discovery of millisecond Robert C. Byrd Green Bank pulsars (MSPs) Telescope (GBT) dedicated Pulsar Timing Historical Context • 1974: Hulse-Taylor binary pulsar – Arecibo telescope @ 430 MHz – ≤ 8 MHz bandwidth – 175 K system temperature Tsky = 25 K @ 430 MHz • 1982: first millisecond pulsar (B1937+21) – Arecibo telescope @ 1400 MHz – 16 MHz bandwidth – 40 K system temperature Tsky ~ 7 K @ 1400 MHz Pulsar Timing Arrays Current Instrumentation Typical Parameters Current Publications • Yardley et al. 2011, “On detection of the • Aeff ~ 100 m stochastic gravitational-wave background using the Parkes pulsar timing array,” Mon. Not. R. • ∆ν ~ 100 MHz Astron. Soc., 414, 1777 • van Haasteren et al. 2011, “Placing limits on the • Tsys ~ 30 K stochastic gravitational-wave background using European Pulsar Timing Array data,” Mon. Not. • ∆t ~ 15 min. R. Astron. Soc., 414, 3117 • Demorest et al. 2012, “Limits on the Stochastic Gravitational Wave Background from the North American Nanohertz Observatory for Gravitational Waves,” ApJ, in press • Hobbs et al. 2010, “The International Pulsar Timing Array project: using pulsars as a gravitational wave detector,” Class. Quant. Grav., 27, 084013 Pulsar Timing Arrays Near Future Instrumentation New Pulsars Receiver performance • Major radio pulsar reaching fundamental surveys worldwide limits, e.g., – HTRU (Parkes) • Tsys ≈ 20 K – GBNCC (GBT) • ∆ν ~ 1 GHz for – GBT drift scan (GBT) observations near 1 GHz – P-ALFA (Arecibo) • Multi-wavelength MSPs a.k.a. Fermi Pulsar Surveys • New pulsars – Add more arms to PTA, increase sensitivity; or – For fixed amount of observing time, improve quality of PTA • Survey algorithm 1. Observe position on sky 2. Search resulting time a “position” on the sky series for periodic signal at period P with dispersion DM, with orbital parameters – Loop Pulsar Timing Arrays Near Future Pulsar Timing Arrays Digital Backends Green Bank Ultimate Pulsar Processing Instrument (GUPPI) BEE2 feeds 8 gaming systems w/NVIDIA GPUs 100, 200, or 800 MHz bandwidth Large improvement in timing precision ~1 TFLOP in real time operational Pulsar Timing Arrays Digital Backends Wide-bandwidth (1+ GHz), real-time systems implementing RFI mitigation and folding (in construction) • Parkes HIPSR • Effelsberg ASTERIX and follow-on • Arecibo PUPPI FPGA-based Reconfigurable Open Architecture Computing Clone of GBT GUPPI Hardware (ROACH), developed by Center for Astronomical Signal Processing and Electronics Research (CASPER, at UC Berkeley) Pulsar Timing Arrays New Feed & Rx Systems • GBT Wide-Band Pulsar System • Arecibo Wide-Band System • Effelsberg Ultra-broad Band • Arecibo AO40 • Five-hundred metre Aperture Spherical Telescope (FAST) multi-feed system 19 or 100 feed horns Tsys ~ 29K (5K sky + 9K spill + 4K coax jct + 7K dewar jct + 3K LNA) Pulsar Timing Arrays New Feed & Rx Systems • GBT Wide-Band Pulsar System • Arecibo Wide-Band System • Effelsberg Ultra-broad Band • Arecibo AO40 • Five-hundred metre Aperture Spherical Telescope (FAST) multi-feed system Quasi Self-Complementary Feed (developed under U.S. SKA auspices) 19 or 100 feed horns Pulsar Timing Arrays New Feed & Rx Systems • GBT Wide-Band Pulsar System • Arecibo Wide-Band System • Effelsberg Ultra-broad Band • Arecibo AO40 • Five-hundred metre Aperture Spherical Telescope (FAST) multi-feed system 19 or 100 feed horns continuous coverage from 600-3000 MHz cryogenically cooled receiver design by Weinreb (JPL/TDP) + MPIfR Pulsar Timing Arrays New Feed & Rx Systems • GBT Wide-Band Pulsar System • Arecibo Wide-Band System • Effelsberg Ultra-broad Band • Arecibo AO40 ?Five-hundred metre Aperture Spherical Telescope (FAST) multi-feed system ? 19 or 100 feed horns Pulsar Timing Arrays New Approaches SRT, Sardinia, Italy The European Pulsar Timing Array (EPTA): 100-m class telescopes Effelsberg 100-m, Germany Lovell, Jodrell Bank, NRT, Nancay, France WSRT, Westerbork, NL UK Ultimately forming the Large European Array for Pulsars (LEAP) Pulsar Timing Arrays New Approaches Large European Array for Pulsars = LEAP! Current status: - Hardware to record >128 MHz BW, 8 bits, baseband data in place at all telescopes: ASTERIX-like systems at Effelsberg and Jodrell Bank PUMA-II at Westerbork BON at Nançay DFB (in APSR-mode, to be tested) at Sardinia - 24 hr observations at L-band, once per month (in addition to regular EPTA observations – 30 TB/site/session) - Data currently shipped by disk, internet tested Almost there: 160-m dish! - Successful addition of EFF-JB-WSRT, Only 10-sec of data: Nançay added within days/week, SRT in Q4/2012 Pulsar Timing Arrays New Approaches For a set of N telescopes and M pulsars, what is the optimal scheduling that maximizes the probability of detecting gravitational waves? Telescope Infrastructure Pulsar Sample Pulsar Timing Arrays Future Telescopes Five-hundred metre Aperture Karl G. Jansky Very Large Array Spherical Telescope (FAST) • Aeff ~ 130 m • Deff ~ 300 m • YUPPI backend • 19 or 100-beam system Clone of GUPPI, implemented in existing • ~ 500 new MSPs discovered (projected) correlator • Primarily a timing instrument Pulsar Timing Arrays Future Telescopes Australian Square Kilometre Karoo Array Telescope Array Pathfinder (ASKAP) (MeerKAT) • Aeff ~ 60 m • Phased array feed ~ 30 deg2 field of view • Primarily a search instrument • Aeff ~ 100 m • Primarily a timing instrument Square Kilometre Array The Global Radio Wavelength Observatory • Originally: “Hydrogen telescope” Detect H I 21-cm emission from Milky Way-like galaxy at z ~ 1 • SKA science much broader ⇒ Multi-wavelength, multi- messenger • On-going technical development • International involvement Radio Telescopes and Astrophysical Sources Search for additional dual SMBHs, progenitors of GW-emitting binary SMBHs • With VLBA now, SKA in the future • Relevant for future space-based GW mission 7 pc B0402+679 Very Long Baseline Array (VLBA) Pulsar Timing Arrays Current and Future Instrumentation • What are the requirements? Instrumental performance improving steadily • What exists? Powerful set of telescopes • What is on the horizon? Increasing capability with new telescopes .