Published on November 15, 2007
Recent Results From The LIGO Search For Periodic Gravitational Waves: A Caltech/MIT collaboration supported by the United States National Science Foundation Recent Results From The LIGO Search For Periodic Gravitational Waves Gregory Mendell, LIGO Hanford Observatory on behalf of the LIGO Scientific Collaboration The Laser Interferometer Gravitational-Wave Observatory http://www.ligo.caltech.edu Sources: LIGO-G060177-00-W Sources Mountain on neutron star Precessing neutron star Accreting neutron star Oscillating neutron star A B C D Credits: A. image by Jolien Creighton; LIGO Lab Document G030163-03-Z. B. image by M. Kramer; Press Release PR0003, University of Manchester - Jodrell Bank Observatory, 2 August 2000. C. image by Dana Berry/NASA; NASA News Release posted July 2, 2003 on Spaceflight Now. D. image from a simulation by Chad Hanna and Benjamin Owen; B. J. Owen's research page, Penn State University. Search methods can detect any type of periodic source. Upper limits are set on gravitational-wave amplitude, h0, of rotating triaxial ellipsoid. Searches: Searches 1. Known pulsars (radio & x-ray) (e.g., Crab pulsar) Position & frequency evolution known (including derivatives, timing noise, glitches, orbit). 2. Unknown neutron stars Nothing known, search over sky position, frequency & its derivatives. 3. Accreting neutron stars & LMXBs (e.g., Sco-X1) Position known; some need search over freq. & orbit. 4. Targeted sky position: galactic center, globular clusters, isolated non-pulsing neutron stars (e.g., Cas A) Search over frequency & derivatives. *Searches 2-4 are computationally expensive: e.g., for obs. time T a coherent search over the sky, f, and df/dt scales as T6 while its sensitivity scales as T1/2; orbital params. or higher derivs. add powers of T. Methods: Methods Track Doppler shift and df/dt Semicoherent Methods StackSlide: add the power Hough: add weighted 1 or 0 PowerFlux: add weighted power Coherent Methods Bayesian Param. Estimation Maximum Likelihood & Matched Filtering Frequency Domain Time Domain Weights depend on both noise and antenna patterns: Methods can include multi-detector data and coincidence steps. Hierarchical Methods: combine the above to maximize sensitivity. GWs from triaxial ellipsoid: GWs from triaxial ellipsoid For upper limits have to select a model. (This is not needed for detection!) Ellipticity, , measures asymmetry in triaxially shaped neutron star. Equatorial Ellipticity All results for this talk are 95% confidence ULs on h0 and . f is the GW freq. Astrophysical predictions & payoff: Astrophysical predictions & payoff Neutron, hybrid or quark stars max. ~ 10-6, 10-5, 10-4 respectively. Blandford/LSC statistical estimate: few 10-24 (100 yrs/birthrate)1/2 Age-based limits, e.g., Cas A (see K. Wette’s presentation) Spindown limits (e.g., Crab pulsar) Accreting Stars Torque balanced by GWs or limit cycles Thermo-Elastic mountains Magnetic mountains R-modes For a summary, see: LIGO Scientific Collaboration, gr-qc/0605028, accepted by Phys. Rev. D (2007). For more on indirect limits and astrophysical payoff see B. Owen poster. S5 Known Pulsar Search: Using data from the first thirteen months of S5 Black curve represents one full year of data for all three interferometers running at design sensitivity Blue stars represent pulsars for which we are reasonably confident of having phase coherence with the signal model. Green stars represent pulsars for which there is uncertainty about phase coherence S5 Known Pulsar Search S5 Crab Pulsar Result : S5 Crab Pulsar Result These results give upper limits for the Crab pulsar of < 2.6x10-4, h0 < 5.0x10-25 using S5 data up to the glitch on 23 Aug. 2006 this value of the ellipticity is now in the range of some of the more speculative equations of state (Owen, 2005) These beat the spindown limit of h0 < 1.4x10-24 by a factor of 2.9 – for canonical moment of inertia I = 1038 kgm2 - we even beat Palomba’s limit Start to constrain the amount of spin-down energy in GWs to less than 10% of overall emitted and known spindown (Palomba, 2000, Santostasi) This is significant: the uncertainties on all non-GW contributions add up to 80% of the total! Moment of inertia is uncertain by about a factor of three, but we can plot the result on the moment of inertia – ellipticity plane to give exclusion regions (Pitkin for the LSC, 2005) PRELIMINARY Estimated Joint Sensitivity S3/S4 limit S5 limit Spindown limit Factor of 2.9 Multi-template Crab Search: Multi-template Crab Search UNDERWAY Known pulsar GW searches track phase assuming fGW=2fEM. If the gravitational radiation time evolution is different from that of the electromagnetic radiation it is possible these could miss the gravitational waves. We are considering mechanisms by which any emitted gravitational waves will differ from the electromagnetic. Consideration of free precession or glitches leads to |fGW-2fEM|/(2fEM) 10-4 and corresponding band for time derivatives of the frequency. Need many templates; for the Crab this is underway. S4 One Month All-Sky Search: Hanford, Livingston, and Multi-IFO Results: S4 One Month All-Sky Search: Hanford, Livingston, and Multi-IFO Results See presentation by A. Sintes PRELIMINARY PowerFlux circular-polarization strain H1 S4 upper limits for band with HW Injected Signal.: R.A. Dec. Simulated Pulsar (h0 ~ 8.4 x 10-24 nearly circ. polarized signal) PowerFlux circular-polarization strain H1 S4 upper limits for band with HW Injected Signal. PRELIMINARY [email protected] S3 Final Results: [email protected] S3 Final Results http://einstein.phys.uwm.edu/ PowerFlux S5 Results: PowerFlux S5 Results PRELIMINARY (Using data from 07 Nov. 2005 through 20 July 2006) [email protected] S4 & S5: [email protected] S4 & S5 UPDATE S4 Post-Processing Nearly Complete! S5 Initial Hierarchical Search is underway: computes the fully coherent multi-ifo maximum likelihood statistic for 25 hr segments containing 40 hrs of Hanford and Livingston 4km IFO data. performs Hough transforms of the results. will eventually include automated follow-up of candidates; could include StackSlide option. Summary Preliminary Results and Plans: Summary Preliminary Results and Plans 1. S4 all-sky, 50-1000 Hz, PowerFlux, StackSlide, Hough search: results in preparation (see presentation by A. Sintes) and start of S5 preliminary all-sky PowerFlux search: Best UL: h0 < few 10-24 . 2. Start of S5 preliminary coherent known pulsar search: Best UL h0 < few 10-26; Best UL: a little less than ~ 10-7; Crab limits beat the spindown limit! 3. S5 targeted sources: Cas A (youngest candidate NS) search is underway. (see presentations by K. Wette & B. Owen) 4. S5 all-sky PowerFlux & Multi-IFO initial Hierarchical [email protected] searches are under way. 5. More searches are under development, e.g., LMXBs. End: End Frequency Modulationand S Parameter: Frequency Modulation and S Parameter Relativistic corrections are included in the actual code For analysis < 1 yr sky points with small S have small doppler variation; harder to distinguish GWs from instrument lines at these points. PowerFlux S5 Preliminary Results: PowerFlux S5 Preliminary Results PRELIMINARY (Using data from 07 Nov. 2005 through 20 July 2006) Published Periodic Search Results To June 2007: Published Periodic Search Results To June 2007 (* Just accepted by Phys. Rev. D) (To appear in Phys. Rev. D) * The LIGO/VIRGO Pulsar Search Joint Working Group has started meeting!: The LIGO/VIRGO Pulsar Search Joint Working Group has started meeting! LIGO/Virgo MOU signed Weekly teleconferences Face-to-Face meetings in March & May 2007 Data Sharing started in May 2007 [email protected]:: [email protected]: http://einstein.phys.uwm.edu/ Like [email protected], but for LIGO/GEO matched-filtered search for GWs from rotating compact stars. Support for Windows, Mac OSX, and Linux clients Our own clusters have thousands of CPUs. [email protected] has many times more computing power at low cost. Crab Pulsar Spindown Limit: Crab Pulsar Spindown Limit Spindown limit assumes all the pulsars rotational energy loss is radiated by gravitational wave We know some energy is emitted electromagnetically and is powering the expansion of the Crab nebula This is poorly constrained and allows room for gravitational wave emission Braking index The braking index of the Crab is n=2.5, not n=3 for purely magnetic dipole radiation, and not n=5 for purely gravitational radiation emission Palomba (2000) allows for a combination of mechanisms to account for this braking index and ends up with a GW spin-down limit which is 2.5 times below the n=5 standard limit. Credits: X-ray: NASA/CXC/ASU/J. Hester et al.; Optical: NASA/HST/ASU/J. Hester et al. Slide23: Met with astronomers & astrophysicist in 2006 at MIT A. Melatos: Magnetic Mountains S. Ransom: Longer term, an “Arecibo in the South” would find and time hundreds of new cluster MSPs... (FAST?); Even longer term, the Square Kilometer Array will find thousands of new pulsars C. Palomba/T. Regimbau: population studies Dunc Lorimer: tip of the iceberg => Slide24: Prospects S4 One Month Semicoherent Search Astrophysical Reach: S4 One Month Semicoherent Search Astrophysical Reach PRELIMINARY Other Work:1. X-ray pulsars (for example J0537-6910 glitchiest pulsar) Working with astronomers to get timing data. 2. LMXB Search3. Proposed Unknown Binary Search4. Globular cluster target5. Code speed up6. SN1987A7. Generalized PowerFlux8. LIGO/VIRGO work.: Other Work: 1. X-ray pulsars (for example J0537-6910 glitchiest pulsar) Working with astronomers to get timing data. 2. LMXB Search 3. Proposed Unknown Binary Search 4. Globular cluster target 5. Code speed up 6. SN1987A 7. Generalized PowerFlux 8. LIGO/VIRGO work. Coherent Matched Filtering: Coherent Matched Filtering Jaranowski, Krolak, & Schutz gr-qc/9804014; Schutz & Papa gr-qc/9905018; Williams and Schutz gr-qc/9912029; Berukoff and Papa LAL Documentation F-statistic Results approaching astrophysical interest (also spindown limits/indirect limits; see B. Owen Poster.): Results approaching astrophysical interest (also spindown limits/indirect limits; see B. Owen Poster.) Blandford (1984) as cited by Thorne in 300 Years of Gravitation; see also LIGO Scientific Collaboration, gr-qc/0508065, accepted by Phys. Rev. D (2005); and LIGO Scientific Collaboration, S2 Maximum Likelihood Search. Slide29: Hough (2 of 92 sky patches shown) (North Pole) (Equator) Determines limit (highest-SNR patch) Hough count Sample 0.25-Hz bands Nature of periodic gravitational waves: Nature of periodic gravitational waves The GW signal from a triaxial pulsar can be modelled as The unknown parameters are h0 - amplitude of the gravitational wave signal - polarization angle of signal; embedded in Fx,+ - inclination angle of the pulsar 0 - initial phase of pulsar (0) In the known pulsar searches we currently look for signals at twice the rotation frequency of the pulsars For blind searches the location in the sky and the source’s frequency and its evolution are search parameters.