Inverse compton scattered merger-nova: late X-ray counterpart of gravitational wave signals from NS-NS/BH mergers
The recent observations of GW170817 and its electromagnetic (EM) counterparts show that double neutron star mergers could lead to rich and bright EM emissions. Recent numerical simulations suggest that neutron star and neutron star/black hole (NS-NS/BH) mergers would leave behind a central remnant surrounded by a mildly isotropic ejecta. The central remnant could launch a collimated jet and when the jet propagating through the ejecta, a mildly relativistic cocoon would be formed and the interaction between the cocoon and the ambient medium would accelerate electrons via external shock in a wide angle. So that the merger-nova photons (i.e., thermal emission from the ejecta) would be scattered into higher frequency via inverse compton (IC) process when they propagating through the cocoon shocked region. We find that the IC scattered component peaks at X-ray band and it will reach its peak luminosity in order of days (simultaneously with the merger-nova emission). With current X-ray detectors, such a late X-ray component could be detected out to 200 Mpc, depending on the merger remnant properties. It could serve as an important electromagnetic counterpart of gravitational wave signals from NS-NS/BH mergers. Nevertheless, simultaneous detection of such a late X-ray signal and the merger-nova signal could shed light on the cocoon properties and the concrete structure of the jet.
Subject headings:gamma-ray burst: general - radiation mechanisms: non-thermal - gravitational waves
Department of Astronomy, Beijing Normal University, Beijing 100875, China; email@example.com
LIGO-VIRGO cooperation has reported five gravitational wave (GW) detection (GW150914, GW151226, GW170104, GW170608 and GW170814) from black hole-black hole (BH-BH) mergers, and one GW detection (GW170817) from neutron star-neutron star (NS-NS) merger, which means a new era of GW astronomy is coming (Abbott et al., 2016a, b, 2017a, 2017b, 2017c, 2017d). Nevertheless, the electromagnetic (EM) counterparts of GW170817 have also been detected, including a weak short-duration gamma-ray burst (SGRB) (Goldstein et al., 2017), an optical/IR transient in the galaxy NGC 4993 (Coulter et al., 2017), a radio counterpart (Hallinan et al., 2017) and a late re-brightening X-ray signal(Troja et al., 2017). A large number of teams across the world contributed to the observation of this source using ground- and space-based telescopes (Abbott et al., 2017e, for details).
Before the detection of GW170817, many associated EM counterparts have already been proposed for NS-NS/BH mergers, and their relative brightness is essentially determined by the properties of the post-merger central remnant object. In general, the merger remnant for NS-BH merger would be a BH. But NS-NS merger could lead to a BH or a supra-massive NS, depending on the total mass of the NS-NS system and the NS equation of state (Dai et al., 2006; Gao & Fan, 2006; Zhang, 2013; Lasky et al., 2014; Gao et al., 2016). It is generally believed that, for both cases, short duration gamma-ray bursts (sGRBs) and their afterglow emission are expected as one of the major EM counterparts of NS-NS/BH mergers (Eichler et al., 1989; Narayan et al., 1992; Berger, 2014). On the other hand, for both cases, the merger remnant would be surrounded by a mildly isotropic, sub-relativistic ejecta (which is composed of the tidally ripped and dynamically launched materials during the merger and the matter launched from the neutrino-driven wind from the accretion disk or neutron star surface (Rezzolla et al., 2011; Rosswog et al., 2013; Bauswein et al., 2013; Hotokezaka et al., 2013; Lei et al., 2013; Fernández et al., 2015; Song & Liu, 2017)). These ejecta are mostly composed of neutron-rich materials, and the radioactivity of these materials and the decay of r-process nuclei would heat the ejecta and then power an optical/IR transient (Li & Paczyński, 1998; Metzger et al., 2010). When the merger remnant being a BH, r-process related radioactivity would be the only heating source, so that the luminosity of the optical/IR transient would be roughly times of the nova luminosity (Metzger et al., 2010). But if the merger remnant being a supra-massive NS, its magnetic dipole radiation could serve as an additional heating sources to the ejecta (Zhang, 2013; Gao et al., 2013), which could easily exceed the r-process power, so that the thermal emission from the ejecta would be significantly enhanced (Yu et al., 2013; Metzger & Piro, 2014). The luminosity of the optical transient, in this case, would be systematically brighter by more than one order of magnitude than the r-process dominating cases (Gao et al., 2017). Since the thermal emission from the ejecta are essentially isotropic and also non-relativistic or mildly relativistic (due to the heavy mass loading) and therefore can be detected from any direction if the flux is high enough (Metzger, 2017, for a review).
All observed GW170817 EM counterparts are predicted, but some of them show unexpected behavior (Metzger, 2017b, for a review). For instance, the fluence () and spectral peak energy ( keV) of the associated SGRB 170817A fall into the lower portion of the distributions of known sGRBs, but its peak isotropic luminosity () is abnormally low comparing with other SGRBs(Goldstein et al., 2017). Considering the relatively large upper limit of binary inclination angle relative to our line of site from the GW signal analysis (Abbott et al., 2017d) and the self-consistency between -ray, X-ray, optical/IR and radio observations, a cocoon emission has been proposed as one of the most concordant model (Kasliwal et al., 2017; Piro & Kollmeier, 2017; Gottlieb et al., 2017b; Xiao et al., 2017). Specifically, when a relativistic jet propagating through the surrounded ejecta, a mildly relativistic cocoon would be formed embracing the main jet, so that some relatively weak emission would be expected from a widen-angle structure located in the peripheral of the jet to explain the observed -ray emission (Gottlieb et al., 2017; Kathirgammaraju et al., 2017; Lazzati et al., 2017a, b). Within this scenario, the X-ray and radio observations could be interpreted as the off-axis afterglow from the main jet or from the on-axis afterglow emission of the cocoon, depending on the concrete properties of the jet and the ejecta (Kasliwal et al., 2017; Piro & Kollmeier, 2017; Gottlieb et al., 2017b; Troja et al., 2017; Guidorzi et al., 2017).
If this interpretation is correct, probably all NS-NS mergers would generate a widen-angle mildly relativistic cocoon. The interaction between the relativistic cocoon and the ambient medium could generate an external shock, where particles are believed to be accelerated, giving rise to broad-band synchrotron radiation (cocoon afterglow emission). Considering the large opening angle of the cocoon, photons from other isotropic emission component, such as the thermal component from the ejecta, would be scattered into higher frequency via Inverse Compton (IC) process when they propagating through the cocoon external shock region. In this work, we will estimate the peak frequency and peak flux for this new emission component, and show the dependability of such emission by the currently available X-ray telescopes.
2. Seed photons from ejecta thermal emission
2.1. Numerical model
Considering that one NS-NS merger event leaves behind a central remnant (either BH or NS), surrounded by a neutron rich ejecta with mass and initial dimensionless speed . In any case, the ejecta would receive heating from the radioactive decay of the heavy nuclei synthesized in the ejecta via the r-process. If the central remnant is a NS, its magnetic dipole radiation could serve as an additional heating source to the ejecta. Moreover, the NS wind would continuously push from behind and accelerate the ejecta. Considering the energy dissipation through sweeping up the ambient medium, the dynamical evolution of the ejecta can be determined by (Yu et al., 2013)
where is the bulk Lorentz factor, is time in the observer frame, is the Doppler factor, is the internal energy in the comoving frame, is time in the comoving frame and is the shock swept mass from the interstellar medium (with density ), where is the radius of the ejecta in the lab frame.
With energy conservation, we have
where is the injected energy from the central engine, is the comoving radioactive power and is the comoving radiated bolometric luminosity. When the central engine is a BH, we normally expect (Ma et al., 2017, for other opinion). When the central engine is a NS, normally a fraction of the dipole radiation luminosity is assumed to be injected into the ejecta (Yu et al., 2013, for details), i,e., , where , with being the spin down luminosity and being the spin down timescale, where , and are the initial spin period, the dipole magnetic strength and the radius of the NS. Throughout the paper, the convention is used for cgs units, except for the ejecta mass , which is in unit of solar mass .
ere we adopt the empirical expression for the comoving radioactive power proposed by Korobkin et al. (2012)
where s and s. The radiated bolometric luminosity could be expressed as
where is the radiation dominated pressure. The comoving volume evolution can be fully addressed by together with .
With being solved, one can easily estimate the evolution of the effective temperature of the ejecta in the comoving frame . If a blackbody spectrum is assumed for the thermal emission from the ejecta, one can calculate the observed flux for a given frequency
where is the Boltzmann constant, is the radiation constant, is the Planck constant and is the luminosity distance.
2.2. Analytical estimation
In principle, one can apply above numerical modeling to calculate the observed flux of the thermal emission from the ejecta (henceforth we call it merger-nova emission) at any time and any frequency. The merger-nova photons could serve as the seed photons to be scattered by the cocoon-medium shock into X-ray band, which will be discussed later in detail. In order to better understand the features of these seed photons, we present some analytic approximation for the results of the merger-nova emission, such as the peaking time of the merger-nova, the peak frequency and peak flux at that time.
Black hole as merger remnant
When the merger remnant is a BH, following the analytical estimation from Metzger et al. (2010), the bolometric luminosity of the r-process-powered merger-nova would reach its peak when the photon diffusion time-scale equals the expansion time-scale, where the radius of ejecta being as
where , and are the mass, velocity and opacity of the ejecta, and for spherical outflow (e.g., Padmanabhan, 2000). Under the assumption of free expansion, would be reached on a time-scale
The peak bolometric luminosity is given as
Using Stefan-Boltzmann law, the effective temperature is roughly estimated as
Assuming a blackbody spectrum, we can estimate the peak frequency of the thermal emission as
The according peak flux could be given as
Massive neutron star as merger remnant
If the equation of state of nuclear matter is stiff enough, the central product for a binary neutron star merger could be a stable or a supra-massive NS rather than a black hole. This newborn massive NS would be rotating with a rotation period on the order of milliseconds, and may also contain a strong magnetic field G similar to “magnetars”. In this case, the magnetar dipole radiation could easily dominate the heating and accelerating process for the ejecta, so that the dynamics of the ejecta could be defined by energy conservation (Gao et al., 2013, for details)
Considering that the number density of ambient medium in the NS-NS merger scenario should be usually low, for most situations, the deceleration time for the ejecta should be larger than the spin down timescale of the NS. At this stage, we have , so that . With the energy injection from the millisecond magnetar, the ejecta could be accelerated into a mildly relativistic or even relativistic speed (Gao et al., 2013). In this case, we have the approximation of . With , we have .
The merger-nova emission peaks at when , where is defined as the optical depth of the ejecta, and is the bulk Lorentz factor of the ejecta, and are the radius and volume of the ejecta in the lab frame. The radius of the ejecta at could be simply estimated as
When , we can estimate the peak time in the observer frame as
where is the radius of the ejecta at spin-down timescale and is the Lorenzt factor of the ejecta at the spin-down timescale and (with for a massive neutron star) is the total spin energy of the millisecond magnetar.
At , the effective temperature of the ejecta in the observer frame would be
where we assume a constant fraction () of the NS dipole luminosity would be thermalized in the outflow. With the numerical modeling, it is found that around the transparent time (i.e., ), the radiated bolometric luminosity of the merger-nova is close to the injection energy power (Yu et al., 2013).
Assuming a blackbody spectrum, we can estimate the peak frequency of the thermal emission in the observer frame as
The according peak flux could be given by
Similarly, for case, we have and when and and when (Gao et al., 2013), where is taken when . Then we can estimate the peak time of mergernova as
In this case, the temperature of the ejecta at would be
The peak frequency of the merger-nova emission in the observer frame is
and the corresponding peak flux is
2.3. Inverse Compton scattering
During the propagation of the cocoon, an external shock would form upon interaction with the ambient medium. The shock-accelerated electrons behind the blast wave are usually assumed to be distributed with a power-law function of electron energy, with a minimum Lorentz factor : . Assuming that a constant fraction of the shock energy is distributed to electrons, the minimum injected electron Lorentz factor can be estimated as
where is the Lorentz factor of the cocoon, and is adopted as commonly used in GRB afterglow modeling (Kumar & Zhang, 2015).
When the seed photons (with frequency ) from the merger-nova propagate through the cocoon external shock region, they would be scattered into higher frequency via IC process. The typical photon frequency of the IC scattered merger-nova would be estimated as
When merger remnant is a BH, we have
When the merger remnant is a NS, we have
The peak flux for the IC scattering component could be estimated as
where is the optical depth for IC scatterings for a constant density medium (Sari & Esin, 2001), and is Thompson scattering cross section. We thus have
for case and case respectively.
We use the numerical method described in section 2.1 to calculate the light curve of the inverse compton scattered merger-nova emission. Given some fiducial parameters, such as the ejecta mass, velocity, opacity being as , , and , and the magnetar parameters being as , , and , we find that the IC scattered merger-nova peaks at X-ray band. In figure 1, we compared the peak flux of the IC scattered merger-nova with the sensitivity of current X-ray facilities, such as Swift/XRT, Chandra and XMM-Newton (Burrows et al., 2005; Weisskopf et al., 2002; Jansen et al., 2001). We find that for NS-BH mergers or NS-NS mergers with BH being the merger remnant, such a X-ray component is only detectable out to 2 Mpc with current available facilities. However, for NS-NS mergers with massive NS being the merger remnant, the IC scattered X-ray component would become much brighter, and it will be detectable out to 200 Mpc (the designed horizon of aLIGO for NS-NS mergers (Abbott et al., 2009)). In figure 1, we also plot the expected IC scattered merger-nova for GW170817, and we find that it is too dim to account for the late Chandra X-ray observations.
4. Conclusion and discussion
The recent observations of GW170817 and its EM counterparts have proven the prediction that NS-NS/BH mergers could lead to rich and bright EM emissions, invoking several emission components. In general, the jet component would give rise to a SGRB and its afterglow, and the isotropic component would give rise to a merger-nova emission and its afterglow. Recent studies suggest that the jet component is structured, with a relativistic jet surrounded by a mild relativistic cocoon. In this case, the interaction between the cocoon and the ambient medium would accelerate electrons via external shock in a wide angle. So that the merger-nova photons would be scattered into higher frequency via IC process when they propagating through the cocoon external shock region.
In this work, we find that the IC scattered component peaks at X-ray band and it will reach its peak luminosity simultaneously with the merger-nova. For NS-BH mergers or NS-NS mergers with BH being the merger remnant, the X-ray component is detectable out to 2 Mpc with current facilities (such as Chandra and XMM-Newton). On the other hand, if the total mass of binary neutron star system is small enough and the equation of state of nuclear matter is stiff enough, the merger of two NSs could leave behind a supra-massive NS. In this case, the merger-nova emission could be significantly enhanced, so that the IC scattered X-ray component also becomes brighter. It will be detectable out to 200 Mpc with current facilities. Note that even for BH remnant case, the magnetic wind driven by Blandford-Payne process (Blandford & Payne, 1982) from new-born BH accretion disk or fallback accretion disk would significantly enhance the merger-nova emission (Chen et al., 2017; Ma et al., 2017, for details), in this case, the IC scattered X-ray component could become as bright as the magnetar remnant case.
Our newly proposed late X-ray emission could serve as an important EM counterpart of GW signals. Simultaneous detection of such X-ray signal and the merger-nova signal could help to investigate the cocoon properties and the concrete structure of the jet.
It is worth noticing that some other mechanism could also generate late X-ray emission, which may outshine the proposed signal here. For instance, if the merger remnant of NS-NS is a supra-massive NS, the X-rays powered by NS wind dissipation would diffuse out at late time when the ejecta becomes (or be close to) optically thin, a late X-ray re-brightening would be expected (Metzger & Piro, 2014; Gao et al., 2015, 2017). But if the supra-massive NS have collapsed into a black hole before the surrounding ejecta becomes transparent, such signal would disappear. On the other hand, the afterglow emission from the structure jet could also provide X-ray photons, but its strength is sensitively depending on the viewing angle and the energy distribution within the jet. For most proper viewing angles, the corresponding jet energy is usually small, so that a relatively weak X-ray afterglow emission is expected (Lazzati et al., 2017b).
- In principle, one needs to firstly transform the seed photon energy from the observer frame to the cocoon-medium shock front frame, then calculate the inverse Compton scattering in that frame, and finally transform the scattered photon energy back to the observer frame. For simplicity, we assume that the seed photon injection direction is aline with the moving direction of the cocoon-medium shock front, in which case the two step relativistic transformations could be canceled out.
- Abbott, B. P., Abbott, R., Adhikari, R., et al. 2009, Reports on Progress in Physics, 72, 076901
- Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2016,PhRvL.116f1102A
- Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2016,PhRvL.116x1103A
- Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017,PhRvL.,118v1101A
- Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017, Arxiv: 1709.09660
- Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017, arXiv:1711.05578
- Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017, PhRvL., 119v161101A
- Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017, ApJL., 848,13
- Andreoni, I., Ackley, K., Cooke, J., et al. 2017, arXiv:1710.05846
- Bauswein, A., Goriely, S., & Janka, H.-T. 2013, ApJ, 773, 78
- Berger, E. 2014, ARA&A, 52, 43
- Blandford, R. D., & Payne, D. G., 1982, MNRAS, 199, 883
- Burrows, D. N., Hill, J. E., Nousek, J. A., et al. 2005, Space Sci. Rev., 120, 165
- Chen, W., Xie, W., Lei, W. H., et al., 2017, accepted for publication in ApJ, arXiv:1709.08285
- Chu, Q., Howell, E. J., Rowlinson, A., et al. 2016, MNRAS, 459, 121
- Coulter, D. A., Foley, R. J., Kilpatrick, C. D., et al. 2017, arXiv:1710.05452
- Dai, Z. G., Wang, X. Y., Wu, X. F., & Zhang, B. 2006, Science, 311, 1127
- Eichler, D., Livio, M., Piran, T., & Schramm, D. N. 1989, Nature, 340, 126
- Fernández, R., Kasen, D., Metzger, B. D., & Quataert, E. 2015, MNRAS, 446, 750
- Jansen, F., Lumb, D., Altieri, B., et al. 2001, A&A, 365, L1
- Gao, H. and Ding, X. and Wu, X.-F. and Zhang, B. and Dai, Z.-G., 2013,ApJ,771,86G
- Gao, H. and Ding, X. and Wu, X.-F. and Dai, Z.-G. and Zhang, B., 2015,ApJ,807,163G
- Gao, H., Zhang, B., Lü, H.-J. 2016, Phys. Rev. D, 93, 044065
- Gao, H., Zhang, B., Lü, H.-J., & Li, Y. 2017, ApJ, 837, 50
- Gao, W.-H., & Fan, Y.-Z. 2006, Chinese Journal of Astronomy and Astrophysics, 6, 513
- Giacconi, R., Rosati, P., Tozzi, P., et al., 2001, ApJ, 551,624G
- Goldstein, A., Veres, P., Burns, E., et al. 2017, arXiv:1710.05446
- Gottlieb, O., Nakar, E., & Piran, T. 2017, arXiv:1705.10797
- Gottlieb, O., Nakar, E., Piran, T., & Hotokezaka, K. 2017, arXiv:1710.05896
- Guidorzi, C., Margutti, R., Brout, D., et al. 2017, arXiv:1710.06426
- Hallinan, G., Corsi, A., Mooley, K. P., et al. 2017, arXiv:1710.05435
- Hotokezaka, K., Kiuchi, K., Kyutoku, K., et al. 2013, Phys. Rev. D, 87, 024001
- Kasen, D., & Bildsten, L. 2010, ApJ, 717, 245
- Kasliwal, M. M., Nakar, E., Singer, L. P., et al. 2017, arXiv:1710.05436
- Kathirgamaraju, A., Barniol Duran, R., Giannios, D., et al. 2017,arXiv:1708:07488K
- Korobkin, O., Rosswog, S., Arcones, A., & Winteler, C. 2012, MNRAS, 426, 1940
- Kotera, K., Phinney, E. S., & Olinto, A. V. 2013, MNRAS, 432, 3228
- Kumar, P., & Zhang, B. 2015, Phys. Rep., 561, 1
- Lasky, P. D., Haskell, B., Ravi, V., Howell, E. J., & Coward, D. M. 2014, Phys. Rev. D, 89, 047302
- Lazzati, D. and Deich, A. and Morsony, B. J. and Workman, J. C. , 2013,ApJ,776L,40Y
- Lazzati, D. and Lopez-Camara, D. and Cantiello, M. and Morsony, B. J. and Perna, R. and Workman, J. C., 2017,arXiv:1709.01468L
- Lei, W. H., Zhang, B. & Liang, E. W. 2013, ApJ, 756, 125
- Li, L.-X. and Paczyński, B., 1998,ApJ,507L,59L
- Ma, S. B., Lei, W. H., Gao, H., et al., 2017, arXiv:1710.06318
- Metzger, B. D. and Martínez-Pinedo, G. and Darbha, S. and Quataert, E. and Arcones, A. and Kasen, D. and Thomas, R. and Nugent, P. and Panov, I. V. and Zinner, N. T., 2010, MNRAS,406.2650M
- Metzger, B. D., & Piro, A. L. 2014a, MNRAS, 439, 3916
- Metzger, B. D. 2017, Living Reviews in Relativity, 20, 3
- Metzger, B. D. 2017, arXiv:1710.05931
- Narayan, R., Paczynski, B., & Piran, T. 1992, ApJ, 395, L83
- Padmanabhan, T., 2013,Theoretical Astrophysics - Volume 1, Astrophysical Processes, Cambridge University Press,22Z
- Pierre, M., Pacaud, F., Adami, C., et al., 2016, A&A, 592,1P
- Piro, A. L., & Kollmeier, J. A. 2017, arXiv:1710.05822
- Rezzolla, L., Giacomazzo, B., Baiotti, L., et al. 2011, ApJ, 732, L6
- Rosswog, S., Piran, T., & Nakar, E. 2013, MNRAS, 430, 2585
- Sari, R., & Esin, A. A. 2001, ApJ, 548, 787
- Song, C.-Y., & Liu, T. 2017, arXiv:1710.00142
- Troja, E., Piro, L., van Eerten, H., et al. 2017, Nature, 551,71T
- Weisskopf, M. C., Brinkman, B., Canizares, C., et al. 2002, PASP, 114, 1
- Xiao, D., Liu, L.-D., Dai, Z.-G., & Wu, X.-F. 2017, arXiv:1710.05910
- Yu, Y.-W. and Zhang, B. and Gao, H., 2013,ApJ,776L,40Y
- Zhang, B. and Mészáros, P., 2001,ApJ,552L,35Z
- Zhang, B., 2013,ApJ,763L,22Z