DRAFT VRSIO DCMBR 12, 2017 Typeset using L A TX twocolumn style in AASTeX61 GW170817 MOST LIKLY MAD A BLACK HOL DAVID POOLY, 1, 2 PAWA KUMAR, 3 AD J. CRAIG WHLR 3 1 Department of Physics and Astronomy, Trinity University, San Antonio, Texas arxiv:1712.03240v1 [astro-ph.h] 8 Dec 2017 2 ureka Scientific, Inc. 3 Department of Astronomy, University of Texas at Austin, Austin, Texas (Received 2017 December 08) ABSTRACT The detection of the neutron-star merger kilonova and short hard gamma-ray burst event GW170817 has been widely described, but one outstanding issue remains: the nature of the remnant of the event. Within the initial uncertainties, the remnant could be either a massive, rotating, magnetic neutron star or a black hole. One of the ways to distinguish these possibilities is with sensitive X-ray observations. We report here Chandra X-ray Observatory Director s Discretionary Time observations made on 2017 Dec 03 and Dec 06 and conclude that X-ray data is consistent with synchrotron radiation in the external shock, and the most likely remnant is a black hole. Keywords: radiation mechanisms: non-thermal - methods: analytical - gamma-rays: bursts, theory - gravitational radiation Corresponding author: David Pooley dpooley@trinity.edu
2 POOLY, KUMAR, & WHLR 1. ITRODUCTIO The discovery of gravitational waves from binary neutron star merger by LIGO (Abbott et al. 2017a) and the associated electromagnetic signal that arrived with a delay of 1.74s (Abbott et al. 2017b) has opened a new and exciting frontier that should provide answers to long-standing questions regarding the synthesis of r-process elements (eg. Tanvir et al. 2017; Kasen et al. 2017), and the generation of relativistic jets and γ-ray photons (cf. Kasliwal et al. 2017) in these events and possibly other high energy sources. The published data for has not enabled a determination as to whether the merged product is a black hole or a massive neutron star. The X-ray flux a couple of weeks after the merger event when the system has become transparent can help answer this question. The Director s Discretionary Time observations with the Chandra X-ray Observatory obtained on - is described in 2. In section 3 we discuss the interpretation of the X-ray data. The main conclusion is summarized in 4. 2. CHADRA OBSRVATIO AD RSULTS Chandra observed the field of starting at 01:39:54 UTC on 2017 Dec 03 for an exposure time of 74.1 ks (ObsID 20860) and again starting at 10:44:40 UTC on 2017 Dec 06 for an exposure time of 24.7 ks (ObsID 20861). Both observations were taken with the telescope aimpoint on the Advanced CCD Imaging Spectrometer (ACIS) S3 chip. Data reduction was performed with the chandra repro script, part of the Chandra Interactive Analysis of Observations (CIAO) software. We used CIAO version 4.9 and calibration database CALDB version 4.7.6. Images of the field of are shown in Figure 1. In each observation, we extracted spectra from a 1. 0 radius region centered on and a 1 source-free background region to the northwest. In the 0.3 8 kev band, we detect 105 counts in the extraction region in the first observation with an estimated contribution of 0.50 background counts and 33 counts in the second observation with an estimated background contribution of 0.16 counts. Combining both observations we report 137.3±12.0 net counts from the two observations combined, for a count-rate of 1.39 ± 0.12 counts / ks. We performed a simultaneous fit of the unbinned source spectra in the 0.5 8 kev band with Sherpa using the modified Cash (1979) statistic cstat and the simplex optimization method. We fit the data with a power-law model with two absorption components. We use the Tuebingen-Boulder ISM absorption model (Wilms, Allen, & McCray 2000) and fix one absorbing column to the Galactic value of n H = 7.20 10 20 cm 2 calculated from the ffelsberg-bonn HI Survey (Winkel et al. 2016) using the online tool at the Argelander-Institut für Astronomie 1. We let the column density of the other absorption component vary to allow for any absorption local to the event. 1 https://www.astro.uni-bonn.de/hisurvey/allsky profiles/ In addition to the simultaneous fit of both the Dec 03 and Dec 06 spectra, we fit each spectrum separately. For consistency, we also fit this model to the spectra we extracted from publicly available Chandra observations from Aug 26 (ObsID 19294) and Sep 01 (ObsID 20728), which we reduced in the manner described above. Our fits agree with the results reported in Troja et al. (2017). We summarize the observations and fit results in Table 1. All uncertainties are the 1σ confidence intervals. The reported fluxes are integrated from the unabsorbed model. Uncertainties on those fluxes were calculated as the standard deviation of integrated, unabsorbed fluxes from Monte Carlo realizations (1000 samples) of the best-fit models, taking into account the uncertainties in the best-fit parameters. The X-ray light curve with uncertainties is shown in Figure 2. We note that, although the spectrum of Dec 03 and the joint fit do prefer a non-zero column for the additional absorption component, the 3σ confidence intervals go down to zero column density so we do not claim that any additional absorption is required by the data. This confidence intervals on the flux are relatively insensitive to this additional component. The large uncertainties on the power-law indices make it difficult to determine whether the X-ray spectrum has changed slope. To address this, we perform a Kolmogorov- Smirnov test on the detected photon energies and find marginal evidence for spectral evolution. The six pairwise KS tests of the four observations show probabilities of the detected photon energies being drawn from the same parent distribution in the 10% 50% range, which indicates consistency of the observed spectra. The smallest probability for being drawn from the same parent distribution is for the first and last observations (ObsIDs 19294 and 20861), which have a 6% probability of being drawn from the same parent distribution. To sum up the early December Chandra observations, the best fit power-law index is 2.2±0.3 (1σ), and the unabsorbed 0.5-8 kev flux is 3.8 10 14 erg cm 2 s 1, which is a factor 5 larger than the X-ray flux 10 days after the merger event. At a distance of 40 Mpc, this corresponds to a luminosity of L x = 7.3 10 39 erg s 1. This luminosity is larger by a factor 20 larger than the ddington luminosity of the merged object (mass 2.7M Abbott et al. 2017a). 3. ITRPRTATIO OF CHADRA DATA The X-ray spectrum f ν ν 1.2±0.3 is consistent with the radio spectrum at a similar epoch (Mooley et al. 2017). Moreover, the 1.6 GHz flux of 1 mjy at 93 days after the merger, extrapolated to the X-ray band with the spectrum ν 0.6, consistent within the errors with both the radio and X-ray spectra, yields flux at 1 kev of 12 njy which is consistent with the Chandra flux at 108 days. The most conservative scenario, therefore, is that the radio and the X-ray photons are produced in the same source via the synchrotron radiation mechanism. The synchrotron cooling frequency in a relativistic shock wave can be shown to be
RMAT 3 Table 1. Power-law Fits to Available Chandra Data on Date ObsID xp. (ks) n H (10 21 cm 2 ) PL Index Unabs. F 0.5 8 kev (erg cm 2 s 1 ) 2017 Aug 26 19294 49.4 3.0 +7.9 3.0 1.0 +0.9 0.8 7.6 +3.5 2.1 10 15 2017 Sep 01 20728 46.7 0.0 +4.0 0.0 1.4 +0.8 0.4 7.4 +2.5 1.7 10 15 2017 Dec 03 20860 74.1 4.3 +2.2 2.0 2.1 +0.4 0.4 4.0 +0.8 0.9 10 14 2017 Dec 06 20861 24.7 3.0 +4.5 3.0 2.4 +0.8 0.7 2.6 +0.6 0.5 10 14 2017 Dec 03/06 joint fit 4.0 +2.0 1.8 2.2 +0.3 0.3 3.8 +0.8 0.8 10 14 2017 Aug 26.4 49.4 ksec 2017 Sep 01.6 46.7 ksec 2017 Dec 03.1 74.1 ksec 0 2 4 6 8 10 12 Figure 1. Chandra images of the field of GW170817. ach image is 30 30 and made with counts in the range 0.3 8 kev. The left and middle images show data presented in Troja et al. (2017), and the right image shows one of the observations presented and discussed here. The colorbar indicates the number of counts in each pixel. F 0.5-8 kev (erg cm -2 s -1 ) 5 10-14 4 10-14 3 10-14 2 10-14 1 10-14 0 0 20 40 60 80 100 120 Days since Figure 2. Chandra light curve of the unabsorbed 0.5 8 kev flux from the available observations of GW170817. Uncertainties in flux take account of uncertainties in the spectral fits; see text for details. (e.g. Kumar & Zhang, 2015) hν c 4.5 kev 1/2 50 n 1 ɛ 3/2 B, 4 t 1/2 7 (1 + Y ) 2, (1) where is the energy in the shocked plasma, n is the number density of particles in the medium in the vicinity of the neutron star merger, ɛ B is the fraction of the thermal energy of the shocked plasma in magnetic fields, t is the time since the merger event, and Y is the Compton Y -parameter which is less than 1 based on radio and X-ray data; integer subscript for a variable X, i.e. X m, is a convenient short hand notation for X/10 m in cgs units. The density n is expected to be of order 0.1 cm 3 or less in the interstellar medium of the host galaxy of. We conclude that hν c > 10 kev, beyond the observed band even a year after the merger. The observed synchrotron flux at 1 kev is given by (e.g. Kumar & Zhang, 2015) f ν 70 njy 1.3 50 ɛ 1.2 e n 1/2 ɛ 0.8 B, 4t 0.9 7, (2) where ɛ e is the fraction of shock energy given to electrons; in deriving the above expression we have taken the electron energy distribution to be dn/de e p with p = 2.2 as suggested by the X-ray and radio spectra f ν ν 0.6. Since the flux at 1 kev at 108 days (12 njy) was larger than the flux at 9 days by factor 5, we conclude from quation 2 that the energy in the shocked plasma increased by a factor 20 between 9 days and 108 days 2. There are two ways 2 The scenario where the rising X-ray light-curve is due to an off-axis relativistic outflow can be ruled out (see Mooley et al. 2017).
4 POOLY, KUMAR, & WHLR that the energy in the external shock could increase with time. One is that the compact remnant left over in the merger continues to pump energy into the outward moving shock front. The other is that there is more energy in the slower moving ejecta that catches up with the decelerating shock front and supplies energy to it. The first possibility is unlikely to work. This is because there is 3 10 3 M of debris between the central object and the forward shock front. Therefore, the speed at which this energy can be supplied is less than 0.5 c (even in the extreme case of 10 52 erg of rotational energy generated by the central engine). This speed is smaller than the outward velocity of the relativistic shock-front, and hence the energy produced by the central engine would not reach the external shock in the time of interest to us. Any such energy will be added to the shock on a much longer time scale when its speed has fallen below 0.5 c, and observations carried out a few years from now would be able to constrain the total energy in the outflow. We are thus left with the second possibility, i.e. that at 108 days there is a lot more energy in slower moving ejecta (d/d(βγ) [βγ] 7 ), which is roughly consistent with the conclusion of Mooley et al. (2017). The data do not provide good constraint in the total energy in the shock front because of the unknown parameters n, ɛ B and ɛ e. 4. DISCUSSIO The Chandra X-ray data for the S merger event GW 170817 obtained 108 days after the merger shows that lightcurve increase in the X-ray band is essentially the same as the radio flux increase reported by Mooley et al. (2017). The spectra at these very different frequencies (1.6 GHz and 2.4x10 8 GHz) are also consistent with each other. These results suggest that radio and X-ray photons are produced in the same source, which is most likely the shocked inter-stellar medium. The rising lightcurve between 10 and 108 days requires the energy in the shocked plasma to increase by a factor 20 during this period, and that is likely due to more energy in slower moving ejecta where the energy scales with shock-front velocity βγ as (βγ) 7, which is roughly consistent with radio data as reported by Mooley et al. (2017). The merger of two neutron stars with mass 1.48±0.12M and 1.26±0.1M where the merged object has a mass of 2.74 +0.04 0.01 M (Abbott et al. 2017a) could result in either a neutron star or a black hole. There might also be a debris disk that gets accreted onto the central object over a period of time, and which could be source of kev X-rays. The observed X-ray luminosity of 7.3 10 39 erg s 1 at 108 days is about 20 times the ddington luminosity for a 2.7 M remnant of the binary S merger. This suggests that X-ray photons are not coming from the debris disk, or fall back disk, left over in the merger. X-rays also could not be produced in a long lived relativistic jet emanating from the compact remnant since the X-ray luminosity in that case is expected to fall off as t 5/3 and the flux at 10 2 days would be much smaller than the observed value. We show next that if the merged object were a hypermassive neutron star endowed with a strong magnetic field, then the X-ray luminosity associated with the dipole radiation would be larger than the observed luminosity 10 days after the event, but much smaller than the observed flux at t 100 days. This argues against the formation of a hypermassive neutron star in this merger. The bolometric spindown luminosity of a neutron star due to dipole radiation is L d (t) (6 10 43 erg s 1 )B 2 12 P 4 3 (1 + t/t SD ) 2, (3) where t SD (5 10 8 s) B 2 12 P 2 3, (4) is the spin-down time, B 12 is the dipole magnetic field in units of 10 12 G, and P 3 is the rotation period in units of milliseconds. The fraction of the spin-down luminosity that appears in the 1-10 kev X-ray band is small. We can use the data for pulsars in our Galaxy to estimate the expected X-ray flux if the merger remnant of GW170817 were a hyper-massive S. The periodic pulses from the Crab are produced relatively close to the neutron star (perhaps near the light cylinder or at a distance of 10 8 cm). The bolometric M luminosity of the pulsed radiation from the Crab is about 0.1% of the spin-down luminosity. Of this flux, only about 10% comes out as 1-10 kev photons (Buhler & Blandford, 2013; Durant et al. 2011). Thus, the fraction of spindown rate of energy loss released as X-ray photons is 10 4 (Durant et al. 2011). Moreover, the fraction of spindown energy going into X-ray photons does not seem to be dependent on pulsar age or surface magnetic field strength (see Table 8 of Durant et al. 2011). Geminga at 3 10 5 years has an X-ray efficiency within a factor 3 of the Crab that is only 10 3 yrs of age. Therefore, a crude estimate is that 10 4 of the spin-down luminosity of the merged neutron star in GW170817 should be expected as 1-10 kev X-rays. Thus, if the merged object is a neutron star, the X-ray flux from GW170817 would be: 3 10 14 erg cm 2 s 1 B12 2 P 3 4 B 12 < 22t 1/2 6 P 3 f x 7 10 9 erg cm 2 s 1 B12 2 t 2 6 B 12 > 22t 1/2 6 P 3 (5) where we took the distance to the object to be 40 Mpc. For typical magnetar magnetic field strength larger than 10 14 G the pulsed X-ray flux at 108 days after the merger is an order of magnitude smaller than the observed flux; however, at t = 10 days the X-ray flux from a putative would at least a factor a few larger than the observed flux for B 2 10 15 G. The system becomes optically thin to X-ray photons roughly a few days after the merger. Further check of the possibility of a magnetar remnant could be provided by continued monitoring of the system in X-rays. If the remnant is really a young magnetar then we should expect to see X-ray outbursts before too long.
RMAT 5 The current X-ray and radio observations 100 days after the merger are best explained by continued emission of the merger-induced shock continuing to propagate into the interstellar medium surrounding GW170817. The early X-ray data suggest that the remnant is not a hyper-massive neutron star with magnetic field between 10 11 G and 2 10 15 G. This suggests that the merged object was most likely a black hole. and. L. Robinson, U. akar and C. Froning for very useful discussions. JCW is supported in part by the Samuel T. and Fern Yanagisawa Regents Professorship. Facility: CXO Software: CIAO (Fruscione et al. 2006), Sherpa (Freeman et al. 2001) We thank Belinda Wilkes and Chandra team for carrying out the Chandra Director s Discretionary Time observations RFRCS Abbott et al, 2017a, PRL 119, 161101 Abbott et al, 2017b, ApJL 848, L12 Buhler, R. & Blandford, R., 2013, arxiv:1309.7046 Cash, 1979, ApJ, 228, 939 Durant, M., Kargaltsev, O., and Pavlov, G. G., 2011, ApJ 743, 38 Freeman, P., Doe, S., & Siemiginowska, A. 2001, Proc. SPI, 4477, 76 Fruscione, A., McDowell, J. C., Allen, G.., et al. 2006, Proc. SPI, 6270, 62701V Kasen, D., Metzger, B., Barnes, J., Quataert,. and Ramirez-Ruiz,., 2017, ature Kasliwal, M.M. et al. 2017, Science Kumar, P. and Zhang, B., 2015, Physics Reports 561, 1 Mooley, K.P. et al., 2017, arxiv:1711.11573 Mösta, P., Ott, C. D., Radice, D., et al. 2015, ature 528, 376 Tanvir et al., 2017, ApJL 848, L27 Troja,. et al., 2017, ature Wilms, Allen and McCray, 2000, ApJ, 542, 914 Winkel, B., Kerp, J., Flöer, L., et al. 2016, A&A, 585, A41