Timescales for Detecting Magnetized White Dwarfs in Gravitational Wave Astronomy
Total Page:16
File Type:pdf, Size:1020Kb
Article Timescales for detecting magnetized white dwarfs in gravitational wave astronomy Surajit Kalita 1* 1 Department of Physics, Indian Institute of Science, Bangalore 560012, India; [email protected] * Correspondence: [email protected] Version February 20, 2021 submitted to Universe 1 Abstract: Over the past couple of decades, researchers have predicted more than a dozen 2 super-Chandrasekhar white dwarfs from the detections of over-luminous type Ia supernovae. 3 It turns out that magnetic fields and rotation can explain such massive white dwarfs. If these 4 rotating magnetized white dwarfs follow specific conditions, they can efficiently emit continuous 5 gravitational waves and various futuristic detectors, viz. LISA, BBO, DECIGO, and ALIA can detect 6 such gravitational waves with a significant signal-to-noise ratio. Moreover, we discuss various 7 timescales over which these white dwarfs can emit dipole and quadrupole radiations and show that 8 in the future, the gravitational wave detectors can directly detect the super-Chandrasekhar white 9 dwarfs depending on the magnetic field geometry and its strength. 10 Keywords: White dwarfs; gravitational waves; magnetic fields; rotation; luminosity 11 1. Introduction 12 Over the years, the luminosities of type Ia supernovae (SNeIa) are used as one of the standard 13 candles in cosmology to measure distances of various astronomical objects. This is due to the reason that 14 the peak luminosities of all the SNeIa are similar, as they are originated from the white dwarfs (WDs) 15 burst around the Chandrasekhar mass-limit (∼ 1.4M for carbon-oxygen non-rotating non-magnetized 16 WDs [1]). However, the detection of various over-luminous SNeIa, such as SN 2003fg [2], SN 2007if 17 [3], etc., for about a couple of decades, questions the complete validity of the standard candle using 18 SNeIa. Calculations show that such over-luminous SNeIa needs to be originated from WDs with 19 masses ∼ 2.1 − 2.8M . It is well understood that these WDs exceed the standard Chandrasekhar 20 mass-limit, and hence, they are termed as the super-Chandrasekhar WDs. Ostriker and Hartwick first 21 showed that rotation could increase a WD’s mass up-to ∼ 1.8M [4]. Later, Mukhopadhyay and his 22 collaborators, in a series of papers for about a decade, showed that magnetic field could be one of the 23 prominent physics, which can significantly increase the mass of the WDs [5–8]. Nevertheless, various 24 other physics such as modified gravity [9], generalized uncertainty principle [10], noncommutative 25 geometry [11], and many more, can also explain these super-Chandrasekhar WD progenitors. 26 So far, no super-Chandrasekhar WDs have been observed directly in any of the surveys such 27 as GAIA, SDSS, Kepler, etc. as the highly magnetized WDs generally possess very low thermal 28 luminosities [12]. Hence, they are difficult to be observed in the standard electromagnetic (EM) 29 surveys. On the other hand, the detection of gravitational waves (GW) by LIGO/Virgo from the 30 various compact object merger events opens a new era in astronomy. In the future, space-based GW 31 detectors such as LISA, DECIGO, BBO, etc., can detect many other astrophysical objects. In this 32 article, we show that these space-based detectors can also detect the continuous GW emitted from the 33 magnetized super-Chandrasekhar WDs, provided they behave like pulsars. We also show that the 34 timescale for the emission of the electromagnetic (EM) and quadrupole radiations from these WDs is 35 highly dependent on the magnetic field geometry and its strength. Submitted to Universe, pages 1 –5 www.mdpi.com/journal/universe Version February 20, 2021 submitted to Universe 2 of 5 36 2. Structure of rotating magnetized white dwarfs A system can emit GW if it possesses a time-varying quadrupolar moment. Hence, neither a spherically symmetric nor an axially symmetric WD can emit GW. It is well known that a toroidal field can z0 Ω z make a star prolate-shaped, whereas rotation and poloidal magnetic fields can transform it into oblate-shaped. Hence, it is necessary to have a misalignment between the magnetic field and rotation axes to make the WD a tri-axial system so that it can emit gravitational χ x0 radiation. For a WD, as shown in Figure1, rotating with an angular χ frequency W, the two polarizations of GW, h+ and h×, at a time t are given by [13,14] 1 1 + cos2 i h+ = h sin c sin 2i cos c cos Wt − sin c cos 2Wt , 0 4 2 x (1) 1 h× = h0 sin c sin i cos c sin Wt − cos i sin c sin 2Wt , 2 Figure 1. A cartoon diagram of magnetized rotating WD pulsar. with 0 2 Here z and z are respectively the 4G W eIxx h = , (2) rotation and magnetic field axes. 0 c4 d 37 where c is the angle between the rotation and magnetic field axes, i is 38 the angle between the rotation axis and the observers’ line of sight, G 39 is Newton’s gravitational constant, c is the speed of light, d is the distance between the detector and 40 the source, and e = jIzz − Ixxj/Ixx with Ixx and Izz being the moments of inertia of the WD about x− ◦ 41 and z− axes respectively. It is evident that if c is close to 0 , radiation at frequency W is dominated, ◦ 42 whereas if c is near 90 , the same at frequency 2W is dominant. To model the rotating magnetized WDs, 43 we use the XNS code [15], which solves the axisymmetric equilibrium configuration of stellar structure 44 in general relativity. A detailed discussion about all the parameters used to model the WDs is given 45 by Kalita and Mukhopadhyay [7]. Figure2 shows the density isocontours for the WDs with purely 46 (a) toroidal and (b) poloidal magnetic field geometry generated with XNS code. It is evident that the 47 toroidal field makes the WD bigger, whereas the poloidal field reduces its size. The stability of a WD 48 in the presence of rotation and magnetic fields can be understood respectively from the constraints on kinetic-to-gravitational energy ratio [16], and magnetic-to-gravitational energy ratio [17]. 3 3 Density (g cm− ) Density (g cm− ) 9 9 10 10 2000 1000 ) 7 ) 7 10 10 m m k k ( 0 ( 0 R R 5 5 10 10 2000 1000 2000 0 2000 1000 0 1000 (km) (km) R R (a) Toroidal magnetic field (b) Poloidal magnetic field 10 −3 Figure 2. Density isocontours of magnetized WDs with central density rc = 10 g cm . 49 Version February 20, 2021 submitted to Universe 3 of 5 50 3. Gravitational radiation emitting from rotating magnetized white dwarfs We model several rotating WDs using XNS code with different central = 1 s = 10 mins = 1 week T T T densities (rc) considering various -16 LISA magnetic field geometry and strength. 10 Massive binaries The relation between the integrated ) Galactic binaries 2 / -18 signal-to-noise ratio (SNR) and the 1 10 TianQin − observation time T is given by [18] z H -20 ALIA ( 10 s D T S 1 -22 SNR = p h, (3) P Extreme S( ) 10 5 n p mass BBO -24 ratio 10 inspirals 51 where S(n) is the power spectral DECIGO -5 -4 -3 -2 -1 0 1 52 density (PSD) of the detector at a 10 10 10 10 10 10 10 ν (Hz) 53 frequency n and h is the strength of 54 GW at that frequency. Figure3 shows p p p Figure 3. S(n) of different GW detectors along with T /5h for 55 S( ) n for thep various detectors [19,20] the various rotating magnetized WDs with different integration time. 56 along with T /5h for our modeled We assume c = 30◦ and d = 100 pc. The typical masses of the WDs 57 48 49 WDs. The masses of the WDs are are 1.40 − 2.65M with radii ∼ 750 − 3300 km, Ixx ∼ 10 − 10 g 58 2 −5 −3 in the range of 1.40 − 2.65M with cm , and e ∼ 10 − 10 depending on rc, magnetic field geometry −3 −5 59 e ∼ 10 − 10 depending on rc, and its strength. 60 magnetic field configuration and its 61 strength. It is evident that BBO and 62 DECIGO can immediately detect the magnetized WDs, whereas the minimum integrated time of 63 observations for ALIA, TianQin and LISA turns out to be a few minutes, a few days and a few weeks, 64 respectively. Moreover, the rotational frequencies for many of these WDs can be significantly higher 65 than the orbital frequencies of galactic binaries, and hence, they are free from these confusion noise 66 [21]. Other problems of confusion noise coming from the massive binaries can be sorted out with 67 proper source modeling, which is beyond the scope of this work. 68 4. Timescale for emission of gravitational radiation Since the WDs, we considered, as shown in Figure1 are pulsar-like, they can emit both EM (dipole) and quadrupole radiations, which are respectively associated with the EM luminosity LD and quadrupole luminosity LGW. Due to these radiations, the variations of W and c with respect to t are given by [22] B2 R6 3 d(WIz0z0 ) 2G 2 5 2 2 p pW 2 = − (Izz − Ixx) W sin c 1 + 15 sin c − sin c F(x ), (4) dt 5c5 2c3 0 B2 R6 2 dc 12G 2 4 3 p pW I 0 0 = − (Izz − Ixx) W sin c cos c − sin c cos c F(x ), (5) z z dt 5c5 2c3 0 0 where Iz0z0 is the moment of inertia of the WD along z −axis, x0 = R0W/c, Rp is the radius of the pole, R0 is the average radius, Bp is the strength of the magnetic field at the pole of the WD, and the function F(x0) is defined as x4 1 ( ) = 0 + F x0 6 4 2 .