arXiv:1804.00453v2 [astro-ph.HE] 23 May 2018 aiyo eetr ol epae a wyi space, in away if far these wish placed from would signals be One then could and detectors. detectors limited of the be family of a will nature in sources very detectors the the near-earth of by our location of the very ability triangulating occur the to Secondly, likely power- are the far. waves of gravitational most of as fac- factor sources main such ful one two is from location Our arise based tors. will earth detectors these on space-based and wavelengths limitation and of ultimate range The wide strengths. with gravita- sources for wave search the tional complement detect further to space-based will Future source detectors ability tremendously. the improve our will of sky globe, localization the in the good detec- with around wave waves up gravitational gravitational set more the and are of more tors identification As clear with [2]. events source such opportu- with probe the to along us neutron nity provided waves, binary has gravitational radiation, of via electromagnetic observation merger spec- recent (BNS) electromagnetic most star the The oth- even in which where signatures trum. detected universe no be can leave the mergers erwise to hole window black of new events a opened has illa otasetdfrain fteN changing NS the which of NS deformations the transient through detector to passing a its lead wave as will for gravitational act probe a should (NS) sensitive of This star a neutron changes. providing the configurational of observable (MI) be inertia may of moment changes pul- the minute very of in stars) timings neutron rotating of (rapidly measurements sars accurate with that fact paper. this in possibility detectors). a earth such near propose conventional We (possibly with source the conjunction of in location accu- the and of detection determination GW rate for precision high with collected bevto fgaiainlwvs(W yLG [1] LIGO by (GW) waves gravitational of Observation anpyisudryn u rpsli ae nthe on based is proposal our underlying physics Main 2 oiBah ainlIsiue riigSho ope,A Complex, School Training Institute, National Bhabha Homi ra Das, Arpan ffr a emd odtc n osbecagsi h signal the in changes possible any be detect d can to events this made pulsar transmit be future will can these look and effort If should events, future. one these foreseeable so that by any for gravita suggests disturbed the importance proposal been For our have crucial mergers. VIRGO), hole of black (and be as LIGO such can radiation in This of additional form provide location. detected can are source GW detectors the if Even pulsar earth. such on monitored detectors, be can which of all ae hs inlcnb oioe nerh ihpsil g p possible can With neighborhoods their earth. in pulsars on s universe, monitored remotely the a be in as can around acts sense, signal this dete whose in be waves de pulsar, should wobbling The and induced resonance signal. at to pulsar pronounced (due most profile extrem be pulse the will its affect effect as will well This as rotation. pulsar its affecting pulsar the rvttoa aepsigtruhapla illa oav a to lead will pulsar a through passing wave gravitational A ,2, 1, usr sWbrgaiainlwv detectors wave gravitational Weber as Pulsars ∗ a away far heas .Dave, S. Shreyansh 1 nttt fPyis hbnsa 515 India 750115, Bhubaneswar Physics, of Institute eetr ol be could detectors ,2, 1, † idiaGanguly, Oindrila hr xetG,n te om frdainwssup- was radiation LIGO of forms holes first other black the no of for GW, merger except possible from where originating be GW not of the could detection of This identification ray precise gamma source. for with observations as supplemented optical such was and merger radiation BNS LIGO recent of by the form example, detection For other mergers. any hole sources emit black for not importance crucial do of which is This the be of location. source determination source will the of accuracy source, about the improving that information direction, additional near detec- provide lie to pulsar they able such when then, based especially earth/space Even tors, by then detected are detectors. directly These conventional be nearby. be to de- to likely the likely sources very of are GW most them galaxy, by as our GW tectable present, within of At are detection pulsars space. the pul- observed in of the reach out of our further inertia extend sources can of This moment astrophysical sar.) any the for disturbs matter, which that event for (Or, waves. tional aay,wihaeacrtl oioe,tu c sa as act thus monitored, of accurately family our are within its (most which in im- pulsars galaxy), variations of get via Thousands which detected pattern. star, being rotation pulses, the its of on oscillations printed with sense, [3–5], induced this detector GW in wave star, gravitational the Weber neutron a The natu- like NS. a acts the matches of frequency at frequency GW pronounced ral the most rota- when be NS resonance will the effects the w.r.t the GW propagation Clearly, allow from of axis. will amplitude direction tion its profile, strain as the pulse well of as the strength in a the changes Hence, rate, probing rotation the the induced. with in be axis, changes along of also rotation detection may pulsar of pulsar combination the the between and of angle direction the wobbling source on GW Depending the inertia. of moment its omto o cuaedtriainof determination accurate for formation aindWbrdtco fgravitational of detector Weber tationed l cuaeymaue pnrt fthe of rate spin measured accurately ely rdce ihacrc hnafocused a then accuracy with predicted inlwv vnsaraydtce by detected already events wave tional oieu aiyof family a us rovide edn ntesuc ieto) The direction). source the on pending uhkiNgr ubi408,India 400085, Mumbai Nagar, nushakti tbeb cuaeosrain fthe of observations accurate by ctable eoeystationed remotely subnevaterplesgasin signals pulse their via isturbance ftoeseicpulsars. specific those of s re hc ontei n other any emit not do which urces ietyb at ae conventional based earth by directly raini h oeto nri of inertia of moment the in ariation 1, aiainlwv ore spread sources wave ravitational o pcfi usr hc would which pulsars specific for ‡ n jtM Srivastava M. Ajit and remote ee eetr fgravita- of detectors Weber detectors ,2, 1, § 2 posedly emitted. Even with a future large network of these being the typical frequencies for BNS merger and earth/space based conventional GW detectors, it is un- also where NS resonant modes are expected). likely to achieve high enough angular resolution for pre- We can write Eij = Ri0j0 in terms of the GW strain cise determination of the location of such distant pure amplitude h for a specific polarization. In transverse GW-sources. In such a case, any other remotely stationed traceless (TT) gauge, detector (such as these pulsar detectors), can provide ad- ditional crucial information. Our proposal also provides −1 T T Rµ0ν0 = ∂0∂0hµν (3) an opportunity to, in some sense, revisit the gravitational 2 wave events which have already been detected by LIGO where (and VIRGO). For this one should look for specific pul- σ T T ikσ x sars which would have been disturbed by these events, hµν = Cµν e (4) and will transmit this disturbance via their pulse signals in any foreseeable future. If these future pulsar events and can be predicted with accuracy then a focused effort can 00 0 0 be made to detect any possible changes in the signals of   0 h+ h× 0 those specific pulsars. This will be a direct test of any Cµν = (5) 0 h× −h+ 0 feasibility of the present proposal. Any possible such de-   00 0 0 tection will give valuable information about the neutron   star parameters (such as equation of state). The suffixes ‘+’ and ‘×’ denote two different polarizations We mention that our proposal is fundamentally differ- of GW. Thus, for the + polarization, the amplitude of re- ent from the conventional method of using pulsar timing sulting Eij is given by (in the transverse-traceless gauge) array for the detection of gravitational waves [6]. In the [8]: latter case, pulse arrival times are monitored for gradual changes due to passage of very long wavelength gravita- 2π2hc2 E = −E = , (6) tional waves in the intervening region. Thus, it is limited xx yy λ2 to extremely low frequency sources with typical frequen- cies less than 10−6 Hz arising from rare sources such as λ being the wavelength of GW. Here we use h to denote supermassive black hole mergers or exotic proposed ob- h+. Resulting changes in the quadrupole moment tensor jects such as cosmic strings. Further, with pulsar timing are given by eq. (1). For simplicity, we take the initial array, to get any information about the direction of GW configuration of NS to be spherically symmetric and the source(s) a network of pulsars will be required. In con- deformation to be ellipsoidal (with the dimension in the trast, in principle, a single pulsar detector as proposed direction of propagation of GW remaining unchanged). here, through changes in its spin rate and pulse profile, Using values Qij obtained from eq. (1), we can then cal- can give some information about the GW source direc- culate changes in the moment of inertia tensor Iij of the tion. neutron star as follows (assuming ellipsoidal deformation A gravitational wave represents periodic variations in of an initially spherical object) the metric of spacetime. Associated changes in the Rie- mann curvature tensor Rµνλρ lead to deformations of a ∆I ∆I k R3c2 xx = − yy ≃ 2 20h (7) matter body through which the GW is passing. For a I I 3 GMλ2 neutron star, changes in its quadrupole moment Qij re- sulting from a GW can be written in the following form where M is the mass of the neutron star. We take sam- [7]: ple values with M = 1.0M⊙ and R = 10 km. Highest sensitivity will be reached for smallest values of λ. (Still Qij = −λdEij , (1) retaining the validity of Eqn.(1) requiring λ much larger than NS radius R, as we will see range of frequencies where Eij is the external tidal field and λd is the tidal we consider are below kHz). As a typical astrophysical deformability given by source of GW, we can take binary NS merger (such as the one recently observed by LIGO and Virgo) with the 2 R5 λd = k2 . (2) highest value of GW frequency being about 1 kHz Then, 3 G we get (taking k2 = 0.1 as a sample value)

Here, k2 is the second Love number with its value for a ∆Ixx NS lying in the range k ≃ 0.05 − 0.15 [2] and R denotes ≃ 10−2h . (8) 2 I the radius of the neutron star. Eqn.(1) is valid for static gravitational fields, hence applicable only for wavelengths The BNS merger recently detected by LIGO-Virgo had larger than the NS radius R. This will hold true for the peak strength of the signal h ≃ 10−19 [2]. Source distance range of frequencies we are interested in (kHz or smaller, has been estimated to be about 130 million light years. 3

We can then imagine that our prototype pulsar detector depinning model of glitches). is in the vicinity of such a GW source with an optimistic Note that the above estimates of changes of spin rate of separation between the source and the pulsar detector NS have not accounted for resonances which can dramat- being about one light year (see, footnote [15]). Resulting ically increase the effects of GW. (For example, resonant value of h at the pulsar location will then be enhanced tidal deformations from orbiting BNS leading to rupture to h ≃ 10−11. Fractional changes in the spin rate ν of of NS crust have been discussed in the literature [10]. For the pulsar will be the same as the fractional change in specific modes, the resonant frequencies for NS can be in the MI. Thus, we estimate changes in the spin rate of the the range of 100 Hz to 1 kHz, which will be relevant for pulsar to be a typical BNS merger GW source.) Under the influence of a gravitational wave the amplitude for a specific mode ∆ν ∆I = ≃ 10−13 . (9) with natural frequency ω can be written as [11], ν I Here, we have written ∆I for the change in the relevant ∝ 1 2i∆Ωt component of MI. Changes in the fractional spin rate of ξ 2 2 e (10) ω − 4∆Ω − 2i∆Ω/tvisc the pulsar of this order should be detectable by preci- sion measurements of the pulses. This is certainly true where ∆Ω = Ω − Ωs, Ω being the driving (GW) fre- for certain millisecond pulsars whose pulse timings have quency and Ωs the spin angular frequency of the pulsar. −15 been measured to an accuracy of about 10 seconds (or tvisc is the time-scale for dissipation. It has been argued better, thereby increasing the sensitivity of these pulsar in [11] that tvisc is much larger than the orbital decay detectors). One may be concerned whether this accuracy time scale tD of the binary neutron star system (even will apply to the changes induced by the GW pulse as much before merger). As the time scale for GW pulse we typically one needs to fold a very large number of pulses are considering will be much smaller, it is possible that to get high accuracy. For this issue it is important to tvisc can be assumed to be large here as well, thereby note that the neutron star acting as a Weber detector leading to significant enhancement of the amplitude at at resonance will exhibit ringing effect. Thus, even for a resonance when ω = 2∆Ω. (For certain alignment be- GW pulse, its effect will be present in the pulsar signal tween the spin axis and the direction of the gravitational for a much longer duration. wave, there may be an averaging effect, which will sup- Clearly, with higher precision measurements of pul- press the effect.) sar timings, weaker, or more distant sources (from the The expression in eq.(10) is not applicable for a GW detector-pulsar) may be detected, increasing the fre- pulse which does not last for many cycles. Here we recall quency of such detections. We have discussed BNS that for Weber detector use of materials with high qual- merger as one particular example of a source. Gravi- ity factor Q was important even for a short GW pulse. tational waves from any other source, with suitable am- For a sustained signal, in principle, one could get enor- plitude may be detected in this manner. Note that the mous amplitude enhancement at resonance with Q factor duration of the changes in the pulse amplitude (and its of about 106. However, even for a short pulse, there is profile) depends on the source and can be large, vary- ringing effect for Weber detector operating at resonance, ing from a millisecond to hundreds of seconds (as for the which was helpful in enhancing signal to noise ratio. Due recent BNS merger). The pulsar timings will then show to this ringing effect, the detector continues to vibrate in variations which may be a sharp variation, or a slow mod- the resonant mode for some time even after the passing of ulation extending over many pulse periods (depending on the pulse through the detector (due to energy absorbed the pulsar spin rate). Further note that for a generic di- from the pulse in the resonant mode). For example, for rection of propagation of the GW, it will induce wobbling a GW burst lasting a few ms, the resonant bar can con- of the pulsar (on top of any already present). This will tinue to ring for a long time, of order 10 min [3, 5]. This lead to a modulation of the pulse intensity profile. De- particular mode having a definite frequency allows for tailed simulations of these changes can give information separation between the random noise and the signal. In about relative direction of GW propagation w.r.t the ro- the same way, if the pulsar continues to ring for some tation axis of NS. Such signals could be searched for in time after the GW pulse has passed through it, then the the recorded data of pulsar timings, looking for sudden radio pulses will continue to retain this definite frequency transient changes, as well as relatively longer modula- signal hidden within. Folding of many pulses may be able tions of the pulses. It may be noted that the GW induced to separate this ringing signal. If there is any such signif- sudden changes should in general be different from the icant gain at resonance then there may be a possibility standard glitches where recovery time scales are much of these pulsar detectors competing with earth based de- longer. We also note that GW induced transient change tectors even for sources in our galaxy. in the spin rate may be positive as well as negative, and We also mention that, in a previous work [12], some so could also appear as an anti-glitch [9] (which cannot be of us have proposed that phase transition dynamics in- explained within the framework of the standard vortex side a pulsar core can lead to density fluctuation changes 4 which can be detected by measuring changes in the pulse deformations by gravitational waves (as proposed here) profiles (and may even account for certain glitches/anti has not been made. glitches). Further, we argued that rapidly evolving den- Let us estimate what level of sensitivity will be re- sity fluctuations can lead to changes in the quadrupole quired for this pulsar-detection of GW if the source lies moment providing new powerful GW sources, possibly in Virgo cluster. With about 2000 galaxies in the Virgo at very high frequencies of order MHz. Combining with cluster, and assuming the rate of BNS merger in a typ- our present results, we can infer that phase transitions ical galaxy ranging from 1 in 10,000 years to 1 in 3000 happening inside any NS core (not necessarily a pulsar) years, one can expect one event in few years originating may provide additional GW sources at very high frequen- from Virgo cluster. With typical distance of source being cies which may be detectable by our presently proposed about 50-60 Mly, we can expect the strain amplitude to remote pulsar GW detectors. be h ∼ 10−18 − 10−19, similar to the one observed by the As globular clusters have a very rich population of neu- recent BNS merger. With eqns.(8)-(9), we expect frac- tron stars, they may be an ideal place to look for such tional change in moment of inertia, and hence the frac- events. In particular, collapsed core globular clusters can tional change in the spin rate to be about 10−20 − 10−21. have inter-star distance as small as 1/10 ly in the core, This is too small to be detected by pulsar timing measure- with larger probability of BNS mergers, or black hole ments (assuming best accuracy for the measurement of mergers, and at the same time, a possibility of small pul- fractional change in pulse timing to be about 10−13) un- sar detector - GW source separation. We may make a less there is a strong effect of resonance. As we mentioned very rough estimate for the expected rates as follows. above, resonance is very likely as typical frequency of GW Note that for eqn.(9), the separation between pulsar and from BNS merger overlaps with the natural frequencies source was taken to be about 1 ly when resonance effects of some of the modes of neutron star oscillations. Taking were not incorporated. If the resonance effects allow the a clue from the Weber bar detector, if the sensitivity can sensitivity to be increased even by a factor of about hun- be enhanced by a factor of 106 or more then pulsars in dred (with ringing of pulsar as we mentioned like for We- Milky way can detect these signals. An enhancement in ber bar) then the distance can easily be allowed to be of the accuracy of measurement of pulsar timings will also the size of a typical globular cluster, that is about 100 help in this direction. Important thing to realize is that, parsec. Most of the pulsars (more that 50%) and BNS are essentially all the pulsars in Milky way will oscillate with inside these globular clusters. With about 2000 pulsars the GW hitting them at widely different times. So the in these globular clusters (which number about 200) in measurements of pulsar timing fluctuations due to a sin- the Milky Way, one can expect that for any BNS merger gle GW event can be carried out over period of many in a globular cluster, there will be few pulsars there which years. A systematic analysis of just the arrival time of will be able to detect the emitted GWs (again, assuming GW at different pulsars can help pin point the location of some resonance effect). The expected rate of BNS merger the GW source. Of course additional information about in the Milky Way is about 1 per 10,000 years and may the source direction can be obtained from any wobbling even be 1 per 3000 years [2] and this estimate suggests induced in the pulsar rotation due to GW induced de- that most of these should be detectable by pulsars. This formations. Clearly, this method will be most useful if estimate does not look optimistic from the point of view one could use pulsars very far away. As globular clus- of observations, especially when we know that very few of ters seem very common in other galaxies, in that sense these pulsars are stable millisecond pulsars with the de- this technique should provide motivation for detailed and sired accuracy. But a further increase in sensitivity factor precision studies of pulsars very far away. We mention by few orders of magnitude, either by stronger resonance that for far away pulsars measurements of quantities like effect, or by a higher precision for pulsar timing mea- timing residuals will not be vary accurate due to lack surement, or simply having a much stronger source, may of detailed knowledge about changes in the intergalactic make such pulsars in the milky way sensitive to events medium etc. (apart from the fact that one needs larger in other nearby galaxies. That may bring the expected telescopes just to detect weak signals from very far away rates of such pulsar-detections within observable reach. pulsars). However, pulsar detection method we have pro- We mention that Dyson had proposed long ago [13] that posed uses changes in pulser timings over very short pe- by using an array of seismometers, earth itself may be riod, and changes in intergalactic medium etc. may not used as gravitational wave detector possibly for gravi- be much relevant for this. tational waves emitted by pulsars. More recently, pro- We conclude by pointing out the main features of our posals have been made for the detection of gravitational proposal that pulsars spread out in space can be used as waves by their effects on nearby stars, especially on the a family of (Weber) gravitational wave detectors. GW is solar acoustic modes (helioseismology and astroseismol- detected by its effect on the moment of inertia of the NS, ogy) [14]. These proposals are very similar in spirit to and resulting signal, in terms of variations in the pulsar our proposals, though to our knowledge, the monitoring pulse timing and its profile (resulting from induced wob- of transient changes in the pulsar signals due to pulsar bling depending on the source direction w.r.t the pulsar 5 rotation axis) is transmitted to earth based pulsar ob- § Electronic address: [email protected] servatories. The effect will be most pronounced at reso- [1] B.P. Abbott et al. (LIGO Scientific Collaboration and nance. It has the possibility of extending our reach fur- Virgo Collaboration), Phys. Rev. Lett. 116, 061102 ther out in space for this new era of gravitational wave (2016). [2] B.P. Abbott et al. (LIGO Scientific Collaboration and astronomy for the detection of GW. Even for relatively Virgo Collaboration), Phys. Rev. Lett. 119, 161101 nearby GW sources, these pulsar detectors, combined (2017). with other conventional GW detectors, can help in im- [3] J. Weber, Phys. Rev. Lett. 18, 498 (1967); ibid Phys. proving the accuracy for determination of precise location Rev. Lett.20, 1307 (1968); ibid Phys. Rev. Lett. 22, 1320 of the GW source. This can be of crucial importance for (1969). sources which do not emit any other form of radiation [4] O. D. Aguiar, Res.Astron.Astrophys. 11, 1 (2011); apart from GW, such as black hole mergers. We also arXiv:1009.1138 [5] M. Maggiore, Gravitational Waves: Volume 1: The- point out that our proposal provides an opportunity to, ory and Experiments Oxford publications, 2007; DOI: in some sense, revisit the gravitational wave events which 10.1093/acprof:oso/9780198570745.001.0001. have already been detected by LIGO (and VIRGO). For [6] G. Hobbs and S. Dai, arXiv:1707.01615 this one should look for specific pulsars which would have [7] T. Hinderer, Astrophys.J. 677, 1216-1220 (2008). been disturbed by these events, and will transmit this dis- [8] S. M. Carroll, Spacetime and Geometry: An Introduction turbance via their pulse signals in any foreseeable future. to General Relativity, Addison Wesley, 2004. 497 If these future pulsar events can be predicted with ac- [9] R. F. Archibald et al., Nature , 591 (2013). [10] D. Tsang, J. S. Read, T. Hinderer, A. L. Piro, R. Bon- curacy then a focused effort can be made to detect any darescu, Phys. Rev. Lett. 108, 011102 (2012). possible changes in the signals of those specific pulsars. [11] D. Lai, Mon. Not. R. Astron. Soc. 270, 611 (1994). This will be a direct test of any feasibility of the present [12] P. Bagchi, A. Das, B. Layek, A. M. Srivastava, Phys.Lett. proposal. Any possible such detection will give valuable B747, 120-124 (2015). information about the neutron star parameters (such as [13] F. Dyson, Astrophys. J. 156, 529, 1969. equation of state). [14] L. Lopes and J. Silk, Astrophys. J. 794, 32 (2014); ibid We thank Kandaswamy Subramanian, Rajaram 807, 135 (2015); ibid, 844, 39 (2017). Nityananda, Balasubramanian Iyer, Joseph Samuel, [15] We note that this is not likely a realistic possibility with the present stage of observations of pulsars as most ob- Raghunathan Srianand, Raghavan Rangarajan and served pulsars are within our galaxy, the most distant be- Bharat Kumar for very useful suggestions and comments. ing about 50 Mly away. Thus, these particular numbers may be taken as a future possibility. We may say that our suggestion should provide motivation for detailed and precision studies of pulsars very far away. For nearby sources, the strength of this technique lies in its abil- ∗ Electronic address: [email protected] ity to give additional information about the GW source † Electronic address: [email protected] direction as mentioned above. ‡ Electronic address: [email protected]