## A UNIFIED BINARY NEUTRON STAR MERGER MAGNETAR MODEL FOR THE CHANDRA X-RAY TRANSIENTS CDF-S XT1 AND XT2

HUI SUN 1 ,YE LI 2 , BIN-BIN ZHANG 3,4 , BING ZHANG 5 , FRANZ E. BAUER 6,7,8 , YONGQUAN XUE 9,10 , WEIMIN YUAN 1,11 Draft version December 3, 2019

## ABSTRACT

Two bright X-ray transients were reported from the Chandra Deep Field South (CDF-S) archival data, namely CDF-S XT1 and XT2. Whereas the nature of the former is not identified, the latter was suggested as an excellent candidate for a rapidly spinning magnetar born from a binary neutron star (BNS) merger. Here we propose a unified model to interpret both transients within the framework of the BNS merger magnetar model. According to our picture, CDF-S XT2 is observed from the 'free zone' where the magnetar spindown powered X-ray emission escapes freely, whereas CDF-S XT1 originates from the 'trapped zone' where the X-ray emission is initially blocked by the dynamical ejecta and becomes transparent after the ejecta is pushed to a distance where Thomson optical depth drops below unity. We fit the magnetar model to the light curves of both transients and derived consistent parameters for the two events, with magnetic field, initial spin period and X-ray emission efficiency being ( Bp = 10 16 G Pi , = 1 2ms, . η = 0 001) and ( . Bp = 10 15 8 . G Pi , = 4 4ms, . η = 0 001) for XT1 and . XT2, respectively. The 'isotropic equivalent' ejecta mass of XT1 is M ej ∼ 10 -3 M glyph[circledot] , while it is not constrained for XT2. Our results suggest that more extreme magnetar parameters are required to have XT1 detected from the trapped zone. The model parameters for both events are generally consistent with those derived from SGRB X-ray plateau observations. The host-galaxy properties of both transients are also consistent with those of SGRBs. The event rate densities of both XT1 and XT2 are consistent with that of BNS mergers. Subject headings:

## 1. INTRODUCTION

The discovery of the first gravitational-wave event GW170817 from binary neutron star (BNS) merger (Abbott et al. 2017a) and its broadband electromagnetic (EM) counterparts (Abbott et al. 2017b) ushered in the multimessenger era of astrophysics. The EM counterparts, including the short gamma-ray burst (SGRB) 170817A, kilonova AT2017gfo, and broadband afterglow, have confirmed the related theoretical models and posed interesting constraints on model parameters (e.g., see Metzger 2017 for a summary). The merger product of this event is not well constrained. Several groups argued that the merger product is a black hole (BH), probably formed after a brief hypermassive neutron star (HMNS) phase (e.g. Margalit &amp; Metzger 2017; Pooley et al. 2018; Rezzolla et al. 2018; Ruiz et al. 2018). However, the possibility of a long-lived supramassive neutron star (SMNS) or stable neutron star (SNS) with a low dipolar magnetic field cannot be ruled out from the data (Ai et al. 2018), and the existence of such a long-lived remnant is helpful to interpret some of the observations (Yu et al. 2018; Li et al. 2018; Piro et al. 2019).

- 1 Key Laboratory of Space Astronomy and Technology, National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, China, hsun@nao.cas.cn;
- 2 Kavli Institute for Astronomy and Astrophysics, Peking University, Beijing 100871, China;
- 3 School of Astronomy and Space Science, Nanjing University, Nanjing 210093, China;
- 4 Key Laboratory of Modern Astronomy and Astrophysics (Nanjing University), Ministry of Education, China
- 5 Department of Physics and Astronomy, University of Nevada, Las Vegas, NV 89154, USA, zhang@physics.unlv.edu;
- 6 Instituto de Astrofísica and Centro de Astroingeniería, Facultad de Física, Pontificia Universidad Católica de Chile, Casilla 306, Santiago 22, Chile;
- 7 Millennium Institute of Astrophysics (MAS), Nuncio Monseñor Sótero Sanz 100, Providencia, Santiago, Chile;
- 8 Space Science Institute, 4750 Walnut Street, Suite 205, Boulder, CO 80301;
- 9 CAS Key Laboratory for Research in Galaxies and Cosmology, Department of Astronomy, University of Science and Technology of China, Hefei, 230026, China;
- 10 School of Astronomy and Space Science, University of Science and Technology of China, Hefei, 230026, China;
- 11 School of Astronomy and Space Science, University of Chinese Academy of Sciences, 19A Yuquan Road, Beijing, 100049, China

The possibility of BNS mergers producing a long-lived neutron star (NS) has been suggested in the literature to interpret some of the X-ray features in SGRB afterglows, including Xray flares (Dai et al. 2006; Gao &amp; Fan 2006), extended emission (Metzger et al. 2008), and especially the so-called internal plateau observed in a good fraction of SGRBs (Rowlinson et al. 2010, Rowlinson et al. 2013, hereafter R13, Lü et al. 2015, hereafter L15, Gao et al. 2016, Li et al. 2016a, 2017, cf. Rezzolla &amp; Kumar 2015). The existence of such long-lived neutron star requires that the equation of state of the neutron star is relatively stiff, which is indeed supported by some recent observations of Galactic massive neutron stars (Lattimer, &amp; Prakash 2010; Cromartie et al. 2019). The X-ray plateau can be powered by the dissipation of a magnetar wind which is essentially isotropic for a rapidly spinning magnetar. In such a case, one would expect BNS-merger powered X-ray transients with no SGRB detection if the line of sight misses the bright jet zone of the event (Zhang 2013; Sun et al. 2017).

Sun et al. (2017) studied the possible X-ray light curves of a long-lived magnetar generated from BNS mergers for different observer's viewing angles. In particular, they defined three geometric zones where the observer would view different light curves:

- · The jet zone. In this zone, the observer would see a bright SGRB. SGRBs with extended emission or internal plateau belong to this configuration;
- · The free zone. In this zone, the observer can see X-rays freely but not γ -rays. There could still be a

GRB 170817A-like weak GRB in that viewing direction which clears a funnel to allow X-rays to escape, but such a SGRB is not detectable if the distance is large enough, e.g. glyph[greaterorsimilar] 80 Mpc in the case of GRB 170817A (Zhang et al. 2018);

- · The trapped zone. In this zone, X-rays are initially trapped behind the ejecta from the BNS merger system, but eventually become free when the ejecta is pushed to a large enough radius.

Assuming that the X-ray radiation efficiency does not sensitively depend on the spindown luminosity, the X-ray light curve in the jet/free zone would generally follow the dipole spindown law, and thus appears as a bright plateau followed by a decay with a temporal index between -1 and -2 if the magnetar does not collapse, or steeper than -3 if it collapses. In the trapped zone, the X-ray emission is initially absorbed until the ejecta reaches the transparent radius, i.e. when the photosphere radius has traversed the ejecta from larger radii to smaller radii. The X-ray light curve, therefore, should show a rapid rise before the transparent time (as the optical depth τ drops below unity) and then follow the dipolar radiation law afterward. For typical magnetar parameters, the spindown timescale is shorter than the transparent time. As a result, the light curves in the trapped zone should follow a decay segment after the rising phase. The transparent time depends on the surface magnetic field and the initial spin period of the magnetar, and the opacity and the mass of the ejecta. Through Monte Carlo simulations, Sun et al. (2017) found that for typical parameters of magnetars and ejecta, the peak luminosity of the X-ray counterparts of BNS mergers should be around 10 46 4 . and 10 49 6 . ergs -1 for the line of sight in the trapped zone and free zone, respectively. One important question is whether such events have been detected by current telescopes and how their event rate density compares with those of other X-ray transients in the universe (Sun et al. 2015).

The 7 Ms Chandra Deep Field South (CDF-S) archive data (Luo et al. 2017) is an excellent resource to search for such transients (Zheng et al. 2017; Yang et al. 2019). Recently, a peculiar X-ray transient, CDF-S XT2, was discovered to be associated with a host galaxy at z ∼ 0 738 (Xue et al. 2019). . Its light curve tracks well the prediction of the spindown luminosity evolution of a millisecond magnetar, with a peak rest-frame 0.3-10 keV luminosity of L peak XT2 , = 2 7 . + 6 3 . -2 3 . × 10 45 ergs -1 . The light curve shows a shallow plateau phase lasting for about 2 ks followed by a ∼ t -2 decay. The source was located in the outskirts of a host galaxy with a low starformation rate. Both the type of the host galaxy and its offset from the center of the galaxy are typical for known SGRBs (Xue et al. 2019). All the data strongly suggest that XT2 originates from a rapidly spinning magnetar formed from a BNS merger.

Besides XT2, there was another bright transient, CDF-S XT1, also discovered in the Chandra Deep Field South survey (Bauer et al. 2017). The light curve shows a quick rise at around 110 s to a peak flux (0.3-10 keV) of 5 1 . × 10 -12 ergs -1 followed by a power-law decay with a slope -1 53 . ± 0 27 (Bauer et al. 2017). . The average spectral slope is Γ = 1 43 . + 0 23 . -0 13 . . No additional X-rays above the background rate have been detected in coincidence with this position from Chandra and XMM Newton archives. The photometric redshift of the host galaxy is z ph = 2 23 (0.39-3.21 at 2 . σ confidence), which leads to a peak luminosity (2-10 keV)

6 8 . × 10 46 ergs -1 ( ∼ 1 -140 × 10 45 ergs -1 over the redshift range). The faint host galaxy ( mR = 27 5 mag) was identified . from the CANDELs survey. Bauer et al. (2017) discussed several possible physical origins of this transient but no conclusive result was claimed. The BNS merger magnetar model was considered as one of the possible scenarios to interpret the data, but other possibilities including a tidal disruption event, a supernova shock breakout and a GRB orphan afterglow cannot be completely ruled out, even though they are disfavored for various reasons.

Motivated by the discovery of XT2, in this paper we reinvestigate XT1 and propose a unified model to interpret both transients within the framework of a BNS merger magnetar model. The detailed description of the model and the light curve fitting are presented in Sections 2 and 3, respectively. In Sections 4 and 5, we compare the properties of these two transients with SGRBs in terms of the X-ray plateau and host properties and estimate the event rate density of these transients. The results are summarized in Section 6.

## 2. MODEL

## 2.1. Magnetar wind emission

Consider a rapidly spinning magnetar produced from a BNS merger. It loses its angular momentum through both magnetic dipole radiation and quadrupole gravitational wave (GW) radiation, with the energy loss rate described as (Shapiro &amp; Teukolsky 1983; Zhang, &amp; Mészáros 2001)

$$\dot { E } = I \Omega \dot { \Omega } = - \frac { B _ { p } ^ { 2 } R ^ { 6 } \Omega ^ { 4 } } { 6 c ^ { 3 } } - \frac { 3 2 G I ^ { 2 } \epsilon ^ { 2 } \Omega ^ { 6 } } { 5 c ^ { 5 } }, \quad \quad ( 1 )$$

where Ω = 2 π/ P is the angular frequency and ˙ Ω its time derivative, I is the moment of inertia, Bp is the dipolar field strength at the magnetic poles on the NS surface, R is the radius of the NS, and glyph[epsilon1] is the ellipticity of the NS. The first term is the magnetic dipole radiation or magnetar wind spindown term, and the second term describes the GW radiation energy loss rate.

Similar to Sun et al. (2017), we parameterize the X-ray emission luminosity due to magnetar wind dissipation as being proportional to the dipole spindown luminosity, i.e.,

$$L _ { X, j \text{et/free} } ( t ) = \eta L _ { \text{sd} } \equiv \eta \frac { B _ { p } ^ { 2 } R ^ { 6 } \Omega ^ { 4 } } { 6 c ^ { 3 } }, \text{ \quad \ \ } ( 2 )$$

where η is the efficiency of converting the dipole spindown luminosity to the observed X-ray luminosity. By solving the Ω -evolution using Eq.(1) and plugging it into Eq.(2), one can obtain the X-ray light curve given a quantified η evolution. For simplicity, we take η as a constant.

In the jet zone or free zone, the light curve starts with a relatively flat plateau, which is followed by a decay of ∝ t -2 (dipole spindown dominated) or ∝ t -1 (GW spindown dominated).

In the trapped zone, the X-ray photons are initially trapped in the ejecta and energize the 'merger-nova' (which is the kilonova with extra energy injection from a magnetar, Yu et al. 2013). They become free after the ejecta reaches the transparent radius. The corresponding X-ray light curve is expected to show a fast rise, followed by an evolution defined by magnetar wind dissipation. There are three critical timescales to determine the shape of light curves: the magnetar spindown time t sd , the transparent time t τ , and the collapse time t c . For t τ &lt; t sd glyph[lessmuch] t c , the X-ray plateau phase (followed by a decay) will emerge during the rising phase. On the other hand, for

t sd &lt; t τ glyph[lessmuch] t c a decaying light curve should immediately follow the rising phase. Depending on which component dominates the spindown, the temporal slope in the decay slope varies in the range of [-2,-1] (Sun et al. 2017).

## 2.2. Light curve in trapped zone

The bolometric luminosity of merger-nova peaks in the ultraviolet/optical/infrared band. The merger-nova is usually dim in X-rays but can be bright under extreme conditions, e.g. an extremely large magnetic field, a very small initial spin period, or a low ejecta mass (Sun et al. 2017, cf. Siegel &amp; Ciolfi 2016a,b). We calculate the merger-nova light curve following Yu et al. (2013) and Sun et al. (2017) (see those papers for detailed treatments). The Lorentz factor of the ejecta evolves as

$$\frac { d \Gamma } { d t } = \frac { L _ { \text{sd} } + L _ { \text{ra} } - L _ { e } - \Gamma \mathcal { D } ( d E ^ { \prime } _ { \text{int} } / d t ^ { \prime } ) } { M _ { \text{ej} } c ^ { 2 } + E ^ { \prime } _ { \text{int} } }, \quad \ \ ( 3 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where L sd , L ra , Le are the spindown luminosity from magnetar, radioactive heating luminosity and bolometric emission of merger-nova, respectively, M ej is the ejecta mass, and E ′ int is the internal energy of the ejecta in comoving frame. D = 1 / [ Γ (1 -β cos θ )] is the Doppler factor, β = √ 1 -Γ -2 is the dimensionless velocity, Γ is the Lorentz factor, and θ is the viewing angle ( θ = 0 for an on-beam observer).

By tracking the comoving temperature of the ejecta T ′ , one can obtain the merger-nova luminosity at a given frequency ν as

$$( \nu L _ { \nu } ) _ { \text{bb} } = \frac { 8 \pi ^ { 2 } \mathcal { D } ^ { 2 } R ^ { 2 } } { h ^ { 3 } c ^ { 2 } } \frac { ( h \nu / \mathcal { D } ) ^ { 4 } } { e x p ( h \nu / \mathcal { D } k T ^ { \prime } ) - 1 }, \quad ( 4 )$$

In general, the total X-ray luminosity from the trapped zone can be calculated as

$$L _ { X, \text{trapped} } ( t ) = e ^ { - \tau } \frac { \eta B _ { p } ^ { 2 } R ^ { 6 } \Omega ^ { 4 } ( t ) } { 6 c ^ { 3 } } + ( \nu _ { X } L _ { \nu, X } ) _ { \text{bb} }, \quad ( 5 ) \quad \text{res} [ \text{mul} ]$$

where

$$\tau = \kappa ( M _ { e j } / V ^ { \prime } ) ( R / \Gamma ). \text{ \quad \ \ } ( 6 ) \text{ \quad e j }$$

is the optical depth of the ejecta, κ is the opacity of the ejecta for X-rays, and V ′ is the comoving volume. The thermal component (the second term) is usually negligibly dim in X-rays. Therefore the luminosity in the trapped zone is dominated by the nonthermal component (the first term of Eq. 5).

## 2.3. Pair production, ionization state, and opacity of the ejecta

One complication in calculating the trapped-zone light curves is the optical depth for nonthermal X-rays. Pair production in the magnetar wind may increase the optical depth of X-ray photons, and the ionization state of the ejecta would significantly affect the opacity (Metzger &amp; Piro 2014).

The effect of pair production can be evaluated by assuming the spectrum of nonthermal X-ray emission from magnetar wind dissipation, which is not well constrained from the observational data of SGRB X-ray plateau emission. Based on GRB phenomenology, we assume a broken-power-law spectrum similar to the Band function of GRB prompt emission (Band et al. 1993). This is because the X-ray plateau emission of SGRBs usually does not extend to gamma-rays, and if they do (detected by Swift BAT as the so-called 'extended emission', the spectrum is quite soft. We assume an E peak =10 keV, β = -2 5 and . α = -Γ p = -1 43 for the band function and .

extend the spectrum to 0.511-10 MeV in the rest frame to estimate the pair production optical depth. This gives

$$\mu _ { \sqrt { L _ { X } } } ^ { \mu } & \simeq \frac { \int _ { 0. 5 1 1 \text{MeV} } ^ { 1 0 \text{MeV} } E N ( E ) d E } { \int _ { 0. 2 \text{keV} } ^ { 1 0 \text{keV} } E N ( E ) d E } \simeq 0. 1 8 \quad \quad ( 7 ) \\ \mu _ { \sqrt { L _ { X } } } ^ { \mu } & \simeq \frac { \mu } { L _ { X } } \times \frac { \mu } { L _ { X } } \times \frac { \mu } { L _ { X } } \times \frac { \mu } { L _ { X } } \times \frac { \mu } { L _ { X } } \times \frac { \mu } { L _ { X } } \times \frac { \mu } { L _ { X } } \times \frac { \mu } { L _ { X } } \times \frac { \mu } { L _ { X } } \times \frac { \mu } { L _ { X } } \times \frac { \mu } { L _ { X } } \times \frac { \mu } { L _ { X } } \times \frac { \mu } { L _ { X } } \times \frac { \mu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \frac { \nu } { L _ { X } } \times \sinh \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \min \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \ max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \text{ }$$

With the total energy E γ glyph[similarequal] L γ ∆ t and ∆ t glyph[similarequal] 100 s, the optical depth for pair production is estimated as

$$\mathfrak { e u } _ { \gamma \gamma } & = n _ { \gamma } \sigma _ { T } R = \frac { \sigma _ { T } L _ { \gamma } \Delta t R } { \epsilon _ { \gamma } ( 4 / 3 ) \pi R ^ { 3 } } = \frac { 0. 1 8 \sigma _ { T } L _ { X } } { \epsilon _ { \gamma } ( 4 / 3 ) \pi c ^ { 2 } \Delta t } \\ \mathfrak { e u } _ { \gamma } & = n _ { \gamma } \sigma _ { T } R = \frac { \sigma _ { T } L _ { \gamma } \Delta t R } { \epsilon _ { \gamma } ( 4 / 3 ) \pi c ^ { 2 } \Delta t } \\ \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { a u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { ele u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { \ e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak{ e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak'e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak \ e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfra { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e n } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e s u _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u  _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u. } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak { e u } _ { \mathfrak$$

where we have adopted the following parameters: R = c ∆ t , ∆ t ∼ 100 s, LX is as adopted as the observed peak luminosity, and glyph[epsilon1] γ = mec 2 = 511 keV is the pair production threshold. Notice that this is a very conservative estimate. With the extreme magnetar parameters needed to fit the data of XT1, the ejecta would reach mild relativistic speeds. This would increase R and decrease the photon energy in the comoving frame. Both effects lead to a decrease in τ γγ , which drops below unity when the bulk Lorentz factor Γ &gt; . 4 3.

Let us conservatively consider that pair production is important in the ejecta. One can compare the total number of produced pairs, N ± , and the total number of electrons already in the ejecta, Ne . They can be estimated as 12

$$\begin{array} { c } ^ { \prime } \\ ( 4 ) & & N _ { \pm } \simeq \frac { \lambda } { \epsilon _ { \gamma } } \lambda ^ { \prime } \lambda ^ { \prime } \lambda ^ { \prime } \lambda ^ { \prime } \\ \epsilon _ { \gamma } & & \text{and} \end{array} \begin{array} { c } ^ { \prime } \\ ( \frac { L _ { X } } { 7 \times 1 0 ^ { 4 6 } \, \text{erg} \, s ^ { - 1 } } \end{array} \begin{array} { c } \Delta t \\ ( \frac { L _ { X } } { 1 0 0 \, s } \end{array} \begin{array} { c }, \\ ( 9 ) \end{array}$$

and

$$N _ { e } \simeq Y _ { e } \frac { M _ { e j } } { m _ { p } } = 4. 8 \times 1 0 ^ { 5 3 } \left ( \frac { M _ { e j } } { 1 0 ^ { - 3 } M _ { \odot } } \right ) \left ( \frac { Y _ { e } } { 0. 4 } \right ), \quad ( 1 0 )$$

respectively, where Ye is the electron fraction. Since the pair multiplicity M ≡ e N ± / Ne glyph[similarequal] 3, one can see that pair production would moderately increase the opacity of X-rays in the ejecta.

Next, one can check the ionization state of the ejecta. With a millisecond magnetar as the central engine, the ejecta is continuously illuminated by the X-ray flux from the magnetar. The ionization energy for the innermost electronic state for Fe is ∼ 9 3 keV. For . Ep ∼ 10 keV, the X-ray number flux at 10 keV is

$$\dot { N _ { X } } & = L _ { X } / ( 1 0 \, \text{keV} ) = 4. 4 \times 1 0 ^ { 5 4 } \, s ^ { - 1 } \, \left ( \frac { L _ { X } } { 7 \times 1 0 ^ { 4 6 } \, \text{erg} \, s ^ { - 1 } } \right ). \\ \tilde { \tilde { n } } & = \Omega \, \text{an can} \, \text{can} \, \text{that} \, \text{aiac} \, \text{can} \, \text{has fill} \, \text{inized} \, \text{by this cross}$$

One can see that the ejecta can be fully ionized by this strong X-ray flux within 0.1 s for M ej = 10 -3 M glyph[circledot] , and in around 1 s even for M ej = 10 -2 M glyph[circledot] . In the following, we assume that the ejecta is fully ionized. 13

12 To calculate N ± , we have assumed that each 511 keV photon is converted to a lepton via γγ → e + e -. Photons with energy greater than 511 keV will interact with lower energy photons to produce pairs. Since the pairproducing γ -rays are in the N ( glyph[epsilon1] ) d glyph[epsilon1] ∝ glyph[epsilon1] -2 5 . d glyph[epsilon1] regime, the pairs are predominantly produced at the energy glyph[epsilon1] γ =511 keV.

13 In order to access whether the ejecta can be fully ionized, one also needs to consider the recombination time scale of the ejecta. The time scale of recombination depends on density and recombination coefficients, which further depends on the ionization state of Fe and temperature (Woods et al. 1981). The ionization rate also depends on the ionization state with the ionization energy ranging from 7.9 eV(Fe I ) to 9.3 keV (Fe XXVI). Detailed calculations using numerical tools such as CLOUDY (Ferland et al. 2017) are needed to further justify the full ionization hypothesis.

The opacity of X-rays is dominated by electron Thomson scattering for a fully ionized ejecta. The X-ray opacity can be estimated as

$$\kappa = \frac { \sigma _ { T } n _ { e } } { \rho } = \frac { \sigma _ { T } \xi _ { e } } { m _ { p } } \simeq ( 0. 6 \, \text{cm} ^ { 2 } \mathfrak { g } ^ { - 1 } ) \frac { \xi _ { e } } { 1. 6 }, \quad \ \ ( 1 2 )$$

where ξ e = (1 + M e ) Ye ∼ 1 6 . for our nominal parameters. In our following modeling, κ = 1 cm 2 g -1 is approximately adopted.

## 3. THE LIGHT-CURVE FIT

Comparing the observed light curves with the model above, one can infer that XT1 likely originates from the trapped zone while XT2 originates from the free zone. We perform the light-curve fitting with the least square method. The free parameters include Bp , Pi and η (for both events) and additional parameters of M ej and the zero time point T 0 for XT1. In Bauer et al. (2017), T 0 is poorly constrained and was arbitrarily set as 10 s prior to the arrival of the first photon and ∼ 150 s before the peak (Bauer et al. 2017). This gap could be larger if the early emission is obscured. It is better constrained for XT2 as no photon was detected 10 s prior to the peak luminosity (Xue et al. 2019).

The light-curve fitting is presented in Figure 1, with the results summarized in Table 1. The rest-frame 0.3-10 keV luminosities are presented for both transients, with XT1 derived from k -correction assuming a spectral slope of 1.43 and XT2 directly from Xue et al. (2019). For XT1 (upper panel of Figure 1), the light-curve data points can be generally reproduced with the trapped zone geometry (red line). An example good fit (with χ 2 /dof close to unity) gives Bp = 10 16 G Pi , = 1 2 ms, . M ej = 0 001 . M glyph[circledot] , η = 0 001 and . T 0 = -140 s. The fast-rising phase corresponds to the emergence of the X-ray emission produced from magnetar wind dissipation as the optical depth drops with time, with the peak at the epoch τ ∼ 1. The decline phase of the XT1 light curve can be well fitted by the spindown luminosity after the spindown time scale up to 10 ks. The ejecta mass is fitted to be 10 -3 M glyph[circledot] .

For XT2 (the lower panel of Figure 1), the light curve can be well fitted by the free-zone model (magenta line) with Bp = 10 15 8 . G Pi , = 4 4 ms, . η = 0 001, which is consistent with . the modeling of this event (Xue et al. 2019; Xiao et al. 2019; Lü et al. 2019). The χ / 2 dof is slightly over 2 but still leads to a relatively small p-value ( p ∼ 0 012 at the significance of . 0.05). The low efficiency of ∼ 10 -3 for both XT1 and XT2 is theoretically expected within the slow magnetic wind dissipation model with a saturation Lorentz factor Γ sat = 10 3 -10 4 (Xiao et al. 2019).

It is interesting to note that the fitting results of XT1 and XT2 suggest comparable magnetar parameters (e.g., magnetic field, initial spin period and transfer efficiency) for the two transients, implying a unified origin for the two transients with different viewing geometries (summarized in Figure 2).

Similar to Xue et al. (2019), who placed an upper limit on the gamma-ray luminosity of a putative SGRB L 8 -100keV = 3 5 . × 10 49 ergs -1 for XT2, we also conducted a search for a possible gamma-ray signal associated with XT1 using the Fermi -GBM data. We searched the signal in a time interval from -10 4 to 2 × 10 4 s around T 0 (from Bauer et al. (2017)) and found no significant transient above 10 keV. The 1σ flux upper limit in 8-100 keV range between -50 s and 50 s is 4 73 . + 8 19 . -4 67 . × 10 -9 ergcm -2 s -1 assuming a power-law spectrum. The 1σ flux upper limit in the same energy range between 200 s and 500 s is 1 2 . + 0 35 . -0 35 . × 10 -8 ergcm -2 s -1 . Both flux upper limits lead to an upper limit luminosity of ∼ 10 50 ergs -1 , which is higher than the luminosity of GRB 170817A but lower than those of the majority of on-axis SGRBs.

FIG. 1.- Light-curve fitting of luminosities in the source rest frame for the two transients, CDF-S XT1 (a) and XT2 (b). The black dots with errors are the binned data with time dilation corrected. (a) For XT1, the data are taken from Bauer et al. (2017). The red curve is an example good fit with the trapped zone model, with the unabsorbed trapped zone luminosity ( η L sd tz , ) marked in gray. (b) For XT2, the data are taken from Xue et al. (2019), and an example free zone fit is shown as the magenta curve.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01107v3-figure-09280f08-78ed-4762-8876-5776499a97ab.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01107v3-figure-f1950f4e-6fa3-45d5-b86a-d250b0c7c1f1.jpg)

The ejecta mass of XT1 in our example fit, M ej = 10 -3 M glyph[circledot] falls into the range of the ejecta mass from numerical simulations of binary neutron star mergers (e.g. Hotokezaka et al. 2013). On the other hand, recent simulations (mostly prompted by the observation of the bright kilonova associated with GW170817) show that the ejecta mass could be as high as 10 -2 M glyph[circledot] , especially when the strong neutrino-driven disk wind is taken into account (Just et al. 2015; Siegel, &amp; Metzger 2017; Fernández et al. 2019b). Our example M ej seems lower than this. We notice that the modeling presented here, unlike the kilonova modeling, does not constrain the total M ej , but only constrains the 'isotropic equivalent' value of M ej along the line of sight. This is because the non-thermal X-ray emitter (the magnetar wind) is moving relativistically. Once the line-of-sight ejecta becomes transparent, the X-ray flux would rise significantly. Since the ejecta mass is not dis-

TABLE 1 SUMMARY OF PROPERTIES OF XT1 AND XT2 AND FITTING RESULTS.

|     | Viewing Direction   | Redshift     | L p / (ergs - 1 )   |   Γ p | B p / (G)   |   P i / (ms) |     η | M ej / ( M glyph[circledot] )   | T 0 / ( s )   | χ 2 / dof   |
|-----|---------------------|--------------|---------------------|-------|-------------|--------------|-------|---------------------------------|---------------|-------------|
| XT1 | Trapped Zone        | 2.23( z ph ) | 6 . 8 × 10 46       |  1.43 | 10 16       |          1.2 | 0.001 | 0.001                           | -140          | 4 . 15 / 7  |
| XT2 | Free Zone           | 0.738        | 2 . 7 × 10 45       |  1.93 | 10 15 . 8   |          4.4 | 0.001 | -                               | -             | 19 . 48 / 8 |

FIG. 2.- Schematic diagram of the geometry highlighting three emission (jet, free and trapped) zones and the possible viewing geometries of XT1 and XT2. A structured jet is shown with gradients in yellow. XT2 (magenta) is likely observed from the free zone with a viewing angle similar to GRB 170817A (green) . XT1 (red) is more likely from the trapped zone.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01107v3-figure-541ddd34-e5da-4fd4-9545-6aff00631e5a.jpg)

tributed isotropically, it is possible to derive an effective M ej smaller than the true ejecta mass if the line of sight is not too deep into the trapped zone.

## 4. COMPARISON WITH SGRBS

Since SGRBs are believed to originate from BNS mergers, as manifested in the GW170817/GRB 170817A association, we perform a comparison of the two X-ray transients (XT1 and XT2) with the previously observed SGRBs.

## 4.1. Magnetic fields and initial spin periods

A fraction (1/3 to 1/2) of SGRBs detected by Swift are followed by extended emission or an X-ray plateau, which can be interpreted as the emission of the dissipating wind of a post-merger magnetar (Rowlinson et al. 2013; Lü et al. 2015). Since the young magnetar wind is essentially isotropic, it is expected that SGRB-less X-ray transients similar to XT1 and XT2 exist and should share similar properties as the X-ray plateaus (Zhang 2013).

We first compare the derived magnetar parameters ( Bp and Pi ) of the SGRB X-ray plateau population and those of XT1 and XT2. For the SGRBs, we adopt the sample of L15 and R13 and their derived magnetar parameters. Since the L15 sample includes many plateaus followed by a steep decay segment (the so-called internal plateaus and best interpreted as the collapse of an SMNS at the end of the plateau), the derived Bp and Pi are only upper limits for these events. Only two candidates in the sample are SNSs with the true fitted values of Bp and Pi . As shown in Figure 3, the SGRB magnetar sample typically has a Bp distribution between 10 15 G and 10 17 G and a Pi distribution between 1 ms and 10 ms. The values for XT1 and XT2 both fall into these ranges, suggesting a similar origin.

## 4.2. Photon indices

In Figure 4, we compare the photon indices of XT1 and XT2 with those of the SGRB X-ray plateaus. Both the R13 sample ( Γ X , 2 for GRBs with two or more breaks and Γ X in the segment with a shallower α for GRBs with one break) and the L15 sample (BAT 15-150 keV photon index Γ γ for the extended emission sample and XRT-band photon index for the X-ray plateau sample) are adopted. One can see that the photon index of XT2 is typical for SGRBs. Even though that of XT1 is not typical, it nonetheless falls within the SGRB photon index distribution.

FIG. 3.- Comparison of the derived magnetar parameters ( Bp and Pi ) of XT1/XT2 with SGRB plateaus. Blue squares mark the magnetar sample from R13. The orange dots mark the magnetar sample from L15, among which those with arrows represent the supramassive population which collapse to black holes so that the data only give the upper limits of Bp and Pi . The upper and right panels are the histograms for Bp and Pi , respectively, with the red and magenta dashed lines marking the transients XT1 and XT2, respectively. The upper limits from L15 are directly taken as face values for the orange histograms.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01107v3-figure-c7b3cb60-9881-48b6-aa67-73dbb0d7760c.jpg)

## 4.3. Host-galaxy properties

We compare the host-galaxy properties of XT1 and XT2 with those of redshift-binned SGRBs and long GRBs (LGRBs) in Figure 5. The host galaxy data of GRBs are adopted from Li et al. (2016b) and references therein (e.g., Berger 2009; Fong et al. 2013; Blanchard et al. 2016). They

FIG. 4.- Comparison of the average photon indices of XT1/XT2 with SGRB plateaus. Red and magenta dashed lines mark the XT1 (Bauer et al. 2017) and XT2 (Xue et al. 2019), respectively. The blue and orange histograms are for the samples of R13 and L15, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01107v3-figure-f74717c2-93d0-4ab7-98f6-c064b9ed6d28.jpg)

are presented in both the specific star formation rate (sSFR) stellar mass ( M ∗ ) plane and the offset - half-light radius R 50 plane. We only include SGRBs with SED-estimated stellar masses and emission-line estimated SFRs in the M ∗ vs. sSFR scatter plot. The host properties of XT2 resemble those of SGRBs at similar redshift ( z ∼ 0.7) in all four quantities. For XT1, the host-galaxy mass appears to be smaller than any among the known SGRBs, while the SFR is comparable to them. Since the redshifts of all known SGRBs are (much) smaller than photometric redshift zph = 2 23 of XT1 (Bauer . et al. 2017), and since galaxies are typically smaller at high redshifts (Bouwens et al. 2004; Mosleh et al. 2012), the XT1 host may not be regarded as abnormally small compared with the SGRB sample. The galaxy size R 50 of XT1 is estimated from the Kron radius r kron as R 50 = r kron / 1 19 by assuming . the Sersic index n = 1, which is generally consistent with latetype galaxies. The R 50 and the offset of XT1, in units of kpc, are estimated with zph = 2 23 (Bauer et al. 2017). . The offset of XT1 belongs to the lower end of the offset distribution of SGRBs. The lower panel of Figure 5 describes O(II : I)host, representing the 'odds' or probabilities that the sources belong to the LGRB (massive-star core collapse type, or Type II) vs. SGRB (compact-star merger type, or Type I) populations based on the statistical properties of the host-galaxy data of the two types (Li et al. 2016b). As pointed out by Xue et al. (2019), the O(II : I)host of XT2 falls right on the peak of the distribution of SGRBs. We similarly calculate O(II : I)host for XT1 and find that it can in principle belong to the SGRB population (high end of distribution) even though it is also consistent with the LGRB population.

## 5. EVENT RATE DENSITY

Based on GW170817 detection, the local event rate density of NS-NS mergers is estimated as 1 5 . + 3 2 . -1 2 . × 10 3 Gpc -3 yr -1 (Abbott et al. 2017a; Zhang et al. 2018). It would be interesting to check how the event rate densities of XT1 and XT2 compares with this rate.

The event rate density in general evolves with redshift as ρ ( z ) = ρ 0 f ( z ). Sun et al. (2015) derived the redshift evolution function f ( z ) of NS-NS mergers considering the merger delay distribution with respect to the star formation history. The local event rate density ρ 0 can be estimated via N = ρ 0 V max Ω T / (4 π ), where V max represents the maximum volume a transient like XT1 or XT2 can be detected by CDF-S, and Ω

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01107v3-figure-2dd64d3e-d5eb-41b9-8daa-acecff632b30.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01107v3-figure-4a5cf25d-4e62-4bea-a961-2a3817b2a51d.jpg)

FIG. 5.- Comparison of the host-galaxy properties of XT1/XT2 with SGRBs (squares) and LGRBs (dots). Different colors mark different redshift ranges. The red/magenta stars represent the host galaxies of XT1/XT2. The upper panel shows the distributions in the sSFR vs. M ∗ domain. The middle panel shows the offset vs. R 50 , both in units of kpc. The low panel shows the histograms of O(II : I)host .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01107v3-figure-147dde9a-eb54-4e08-8d55-215124875897.jpg)

and T represent the field of view (FOV) and on-sky exposure time, respectively. Even though both transients were discovered from CDF-S, different search strategies were performed by the two discovery teams. Different criteria lead to different limits of sky coverage and exposure time.

For the peak flux of XT1, Chandra can detect similar events even if they are fainter by a factor of 10. Therefore, given the peak luminosity of L 2 -10keV = 6 8 . × 10 46 ergs -1 , Chandra can detect similar events to a much higher redshift ( z max = 3 2), . with a maximum redshift-corrected comoving volume 14 ∼ 3000 Gpc 3 . The total on-sky exposure time given by Bauer et al. (2017) is the combination of 46.6 Ms for four of the six ACIS-I detectors, with a FOV of 289arcmin , 62.1 Ms for 2 three ACIS-S detectors with a FOV of 217arcmin 2 , and 3.7 Ms of central 100arcmin 2 for HRC. This gives ρ , 0 XT1 = 54 + 124 -45 Gpc -3 yr -1 (with 1 σ errors hearafter). This is much smaller than ρ , 0 BNS.

For the peak flux of XT2, Chandra can detect similar events up to ∼ z max = 1 9, . with a maximum redshift-corrected comoving volume ∼ 2200 Gpc 3 . The total on-sky exposure time given by Xue et al. (2019) is 7 Ms, with a FOV of 0.05 square degrees. This gives 15 ρ , 0 XT2 = 1 4 . + 3 3 . -1 2 . × 10 3 Gpc -3 yr -1 . The XT1 was also discovered during the search for XT2. The corresponding event rate density of XT1 following this strategy is ρ , 0 XT1 = 1 1 . + 2 5 . -0 9 . × 10 3 Gpc -3 yr -1 . Both the event rate densities of the XT1 and XT2 are consistent with ρ , 0 BNS.

In order to make transients XT1 and XT2, a BNS merger needs to leave behind either an SMNS or an SNS. Based on various constraints, this fraction is generally above 60% (Gao et al. 2016; Margalit &amp; Metzger 2019). Our example fits suggest that the last observed data point for both cases is after the characteristic spindown time scale. Since there is no evidence of magnetar collapse from the light curve (which would appear as a much steeper decay), there is a high probability that the remnant for both cases is an SNS. There is a disagreement on the fraction of SNS remnants for BNS mergers. Margalit &amp; Metzger (2019) suggested that it is less than 3% based on the assumption that the merger remnant of GW170817 is a BH followed by a brief HMNS phase and that the maximum mass of a non-spinning neutron star ( M TOV) is 2.17 M glyph[circledot] (Margalit &amp; Metzger 2017; Ruiz et al. 2018; Rezzolla et al. 2018). If this is the case, our interpretation of XT1 and XT2 as BNS merger remnants would be challenged. On the other hand, this low M TOV is inconsistent with the short GRB X-ray plateau data, which suggest M TOV &gt; . 2 3 M glyph[circledot] and SNS fraction ∼ 30% (Gao et al. 2016). The merger remnant of GW170817 is not well settled (cf. Ai et al. 2018; Yu et al. 2018; Li et al. 2018; Piro et al. 2019). Even if the merger product is a BH, Shibata et al. (2019) recently showed that M TOV upper limit should be 2 3 . M glyph[circledot] rather than 2.17 M glyph[circledot] . This is close to the value from the short GRB X-ray plateau constraint. If so, the derived event rate densities from XT1 and XT2 are still broadly consistent with BNS mergers with an SNS remnant.

The fact that the event rate density of XT1 is lower than that of XT2 may seem unexpected, since the trapped zone solid angle is typically much larger than that of the free zone (Sun et al. 2017). On the other hand, the majority of the trapped zone X-ray transients should have much lower luminosities than XT1. In order to be still detectable in the trapped zone,

14 Here we take the SGRB redshift distribution considering a Gaussian merger model, i.e., Equation (20) in Sun et al. 2015.

15 It is consistent with the event rate density derived by Xue et al. (2019), i.e. ρ 0 XT2 , = 1 8 . + 4 1 . -1 6 . × 10 3 Gpc -3 yr -1 .

the magnetar parameters should be even more extreme so that the plateau luminosity is higher. Indeed XT1 has a stronger Bp and a shorter Pi than XT2, suggesting more extreme conditions. The lower event rate density of XT1 is therefore understandable.

## 6. CONCLUSIONS

In this paper, we have proposed a unified magnetar model of BNS merger to explain both X-ray transients discovered from the Chandra Deep Field-South Survey, i.e. CDF-S XT1 and CDF-S XT2. The model can explain well the observed light curves of the two events, with XT1 from the trapped zone (fitting parameters Bp = 10 16 G Pi , = 1 2ms, . η = 0 001, . M ej = 0 001 . M glyph[circledot] , T 0 = -140 s) and XT2 from the free zone (fitting parameters Bp = 10 15 8 . G Pi , = 4 4ms, . η = 0 001). . This suggestion is supported by the consistency of the properties of these two events with those of the SGRB X-ray plateaus in terms of light curves and spectra as well as the host-galaxy properties. The estimated event rate densities of these transients are consistent with that of BNS mergers (Abbott et al. 2017a; Zhang et al. 2018), suggesting that most of the BNS mergers may have left behind long-lived massive neutron stars. The X-ray opacity which was assumed to be 1 cm 2 g -1 in this work needs to be proven by future simulations.

It is predicted that future multimessenger observations of BNS merger events may catch more such X-ray transients associated with BNS gravitational wave events. Whereas SGRBs are beamed, the magnetar-powered X-ray transients have much wider solid angles, so that most BNS mergers may be associated with a SGRB-less X-ray transient (Zhang 2013; Sun et al. 2017). Such transients are not easy to detect with the current wide-field GRB detectors but could be ideal targets for future wide field X-ray telescopes such as the Einstein Probe (Yuan et al. 2016). The detections of these transients in the future will play a key role in identifying the BNS merger remnant, constraining the NS equation of state, and probing the energy power of kilonova/merger-nova.

Finally, even though we proposed a unified model for CDFXT1 and XT2, it is still possible that the two have distinct origins. Furthermore, it may be possible to interpret both events with scenarios other than BNS origins. For example, Fernández et al. (2019a) argued that XT2 can also be interpreted as the merger product of a white dwarf and a neutron star or black hole. Peng et al. (2019) suggested that tidal disruption of white dwarfs by an intermediate black hole may be possible to interpret both XT1 and XT2. Our proposed model can be eventually confirmed via joint GW / X-ray detections in the future.

We thank the referee for important comments. H. Sun thanks helpful discussions on the magnetar sample with Houjun Lü and on lightcurve fitting with He Gao. H.S. and W.M.Y. acknowledge support by the Strategic Pioneer Program on Space Science, Chinese Academy of Sciences, Grant No.XDA15052100 and the Strategic Priority Research Program of the Chinese Academy of Sciences Grant No. XDB23040100. Y.L. is supported by the KIAA-CAS Fellowship, which is jointly supported by Peking University and Chinese Academy of Sciences. This work is also partially supported by the China Post-doctoral Science Foundation (No. 2018M631242). B.-B.Z. acknowledges support from National Thousand Young Talents program of China, the National Key Research and Development Program of China

(2018YFA0404204), and NSFC-11833003. F.E.B. acknowledges support from CONICYT-Chile Basal AFB-170002 and the Ministry of Economy, Development, and Tourism's Millennium Science Initiative through grant IC120009, awarded to The Millennium Institute of Astrophysics, MAS. Y.Q.X. acknowledges support from NSFC grants (11890693 and 11421303).

## REFERENCES

Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017a, Physical Review Letters,

119, 161101 Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017b, ApJ, 848, L12 Ai, S., Gao, H., Dai, Z.-G., et al. 2018, ApJ, 860, 57 Band, D., Matteson, J., Ford, L., et al. 1993, ApJ, 413, 281 Bauer, F. E., Treister, E., Schawinski, K., et al. 2017, MNRAS, 467, 4841 Berger, E. 2009, ApJ, 690, 231 Blanchard, P. K., Berger, E., &amp; Fong, W.-f. 2016, ApJ, 817, 144 Bouwens, R. J., Illingworth, G. D., Blakeslee, J. P., et al. 2004, ApJ, 611, L1 Cromartie, H. T., Fonseca, E., Ransom, S. M., et al. 2019, Nature Astronomy, 439 Dai, Z. G., Wang, X. Y., Wu, X. F., &amp; Zhang, B. 2006, Science, 311, 1127 Ferland, G. J., Chatzikos, M., Guzmán, F., et al. 2017, Rev. Mexicana Astron. Astrofis., 53, 385 Fernández, R., Margalit, B., &amp; Metzger, B. D. 2019a, MNRAS, 488, 259 Fernández, R., Tchekhovskoy, A., Quataert, E., et al. 2019b, MNRAS, 482, 3373 Fong, W., Berger, E., Chornock, R., et al. 2013, ApJ, 769, 56 Gao, H., Zhang, B., Lü, H.-J. 2016, Phys. Rev. D, 93, 044065 Gao, W.-H., &amp; Fan, Y.-Z. 2006, Chinese J. Astron. Astrophys., 6, 513 Hotokezaka, K., Kiuchi, K., Kyutoku, K., et al. 2013, Phys. Rev. D, 87, 024001 Just, O., Bauswein, A., Ardevol Pulpillo, R., et al. 2015, MNRAS, 448, 541 Lattimer, J. M., &amp; Prakash, M. 2010, arXiv e-prints, arXiv:1012.3208 Li, A., Zhang, B., Zhang, N.-B., et al. 2016a, Phys. Rev. D, 94, 083010 Li, A., Zhu, Z.-Y., &amp; Zhou, X. 2017, ApJ, 844, 41 Li, Y., Zhang, B., &amp; Lü, H.-J. 2016b, ApJS, 227, 7 Li, S.-Z., Liu, L.-D., Yu, Y .-W., et al. 2018, ApJ, 861, L12 Luo, B., Brandt, W. N., Xue, Y. Q., et al. 2017, ApJS, 228, 2 Lü, H.-J., Zhang, B., Lei, W.-H., Li, Y., &amp; Lasky, P. D. 2015, ApJ, 805, 89 Lü, H.-J., Yuan, Y., Lan, L., et al. 2019, arXiv:1904.06664 Margalit, B., &amp; Metzger, B. D. 2017, ApJ, 850, L19 --. 2019, ApJ, 880, L15

- Metzger, B. D. 2017, arXiv:1710.05931
- Metzger, B. D., &amp; Piro, A. L. 2014, MNRAS, 439, 3916 Metzger, B. D., Quataert, E., &amp; Thompson, T. A. 2008, MNRAS, 385, 1455 Mosleh, M., Williams, R. J., Franx, M., et al. 2012, ApJ, 756, L12 Peng, Z.-K., Yang, Y.-S., Shen, R.-F., et al. 2019, ApJ, 884, L34 Piro, L., Troja, E., Zhang, B., et al. 2019, MNRAS, 483, 1912 Pooley, D., Kumar, P., Wheeler, J. C., &amp; Grossan, B. 2018, ApJ, 859, L23 Rezzolla, L., &amp; Kumar, P. 2015, ApJ, 802, 95 Rezzolla, L., Most, E. R., &amp; Weih, L. R. 2018, ApJ, 852, L25 Rowlinson, A., O'Brien, P. T., Tanvir, N. R., et al. 2010, MNRAS, 409, 531 Rowlinson, A., O'Brien, P. T., Metzger, B. D., Tanvir, N. R., &amp; Levan, A. J. 2013, MNRAS, 430, 1061 Ruiz, M., Shapiro, S. L., &amp; Tsokaros, A. 2018, Phys. Rev. D, 97, 021501 Shapiro, S. L., &amp; Teukolsky, S. A. 1983, Research supported by the National Science Foundation. New York, Wiley-Interscience, 1983, 663 p., Shibata, M., Zhou, E., Kiuchi, K., &amp; Fujibayashi, S. 2019, PRD, 100, 023015 Siegel, D. M., &amp; Ciolfi, R. 2016a, ApJ, 819, 14 Siegel, D. M., &amp; Ciolfi, R. 2016b, ApJ, 819, 15 Siegel, D. M., &amp; Metzger, B. D. 2017, Phys. Rev. Lett., 119, 231102 Sun, H., Zhang, B., &amp; Li, Z. 2015, ApJ, 812, 33 Sun, H., Zhang, B., &amp; Gao, H. 2017, ApJ, 835, 7 Woods, D. T., Shull, J. M., &amp; Sarazin, C. L. 1981, ApJ, 249, 399 Xiao, D., Zhang, B.-B., &amp; Dai, Z.-G. 2019, ApJ, 879, L7 Xue, Y. Q., Zheng, X. C., Li, Y., et al. 2019, Nature, 568, 198 Yang, G., Brandt, W. N., Zhu, S. F., et al. 2019, MNRAS, 487, 4721 Yu, Y.-W., Liu, L.-D., &amp; Dai, Z.-G. 2018, ApJ, 861, 114 Yu, Y.-W., Zhang, B., &amp; Gao, H. 2013, ApJ, 776, L40 Yuan, W., Amati, L., Cannizzo, J. K., et al. 2016, Space Sci. Rev., 202, 235 Zhang, B. 2013, ApJ, 763, L22 Zhang, B., &amp; Mészáros, P. 2001, ApJ, 552, L35 Zhang, B.-B., Zhang, B., Sun, H., et al. 2018, Nature Communications, 9, 447

Zheng, X. C., Xue, Y. Q., Brandt, W. N., et al. 2017, ApJ, 849, 127
Hui Sun Ye Li Binbin Zhang Bing Zhang Franz E. Bauer Yongquan Xue Weimin Yuan
A unified binary neutron star merger magnetar model for the Chandra X-ray transients CDF-S XT1 and XT2
Two bright X-ray transients were reported from the Chandra Deep Field South
archival data, namely CDF-S XT1 and XT2. Whereas the nature of the former is
not identified, the latter was suggested as an excellent candidate for a
rapidly spinning magnetar born from a binary neutron star (BNS) merger. Here we
propose a unified model to interpret both transients within the framework of
the BNS merger magnetar model. According to our picture, CDF-S XT2 is observed
from the "free zone" where the magnetar spindown powered X-ray emission escapes
freely, whereas CDF-S XT1 originates from the "trapped zone" where the X-ray
emission is initially blocked by the dynamical ejecta and becomes transparent
after the ejecta is pushed to a distance where Thomson optical depth drops
below unity. We fit the magnetar model to the light curves of both transients
and derived consistent parameters for the two events, with magnetic field,
initial spin period and X-ray emission efficiency being ($B_p=10^{16}\,G$,
$P=1.2\,\rm ms$, $\eta = 0.001$) and ($B_p=10^{15.8}\,G$, $P=4.4\, \rm ms$,
$\eta = 0.001$) for XT1 and XT2, respectively. The "isotropic equivalent"
ejecta mass of XT1 is $M_{\rm ej} \sim 10^{-3}$ $M_{\odot}$, while it is not
constrained for XT2. Our results suggest that more extreme magnetar parameters
are required to have XT1 detected from the trapped zone. The model parameters
for both events are generally consistent with those derived from SGRB X-ray
plateau observations. The host galaxy properties of both transients are also
consistent with those of SGRBs. The event rate densities of both XT1 and XT2
are consistent with that of BNS mergers.