Open Access
Issue
A&A
Volume 683, March 2024
Article Number A243
Number of page(s) 27
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/202347629
Published online 22 March 2024

© The Authors 2024

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1. Introduction

From their first three observational runs, the LIGO/VIRGO gravitational wave (GW) detectors discovered 90 binary black holes (BBHs), two binary neutron stars (BNS), and two binary neutron star-black hole (NS-BH) systems (Nitz et al. 2023). Among these, the unique discovery of GW 170817 (Abbott et al. 2017a) and its associated broadband electromagnetic (EM) detections (Abbott et al. 2017b,c) marked the arrival of the multi-messenger astronomy era with the detection of the short gamma-ray burst (SGRB) GRB 170817A (Goldstein et al. 2017; Savchenko et al. 2017; Abbott et al. 2017c), the optical and near-infrared (NIR) kilonova AT2017gfo (Arcavi et al. 2017; Lipunov et al. 2017; Coulter et al. 2017; Tanvir et al. 2017; Soares-Santos et al. 2017; Valenti et al. 2017; Abbott et al. 2017d), and the broadband radio-to-X-ray afterglow (Alexander et al. 2017; Evans et al. 2017; Hallinan et al. 2017; Margutti et al. 2017; Troja et al. 2017; Savchenko et al. 2017; Abbott et al. 2017c; Sugita et al. 2018), which collectively reinforced related theoretical models and provided novel parameter constraints (e.g., Li & Paczyński 1998; Metzger et al. 2010; Goriely et al. 2011; Roberts et al. 2011; Barnes & Kasen 2013; Fernández & Metzger 2016; Smartt et al. 2017; Metzger 2019).

The merger remnant of an NS-BH binary is expected to be a BH, while four different remnant possibilities have been argued to result from a BNS merger: a BH; a differential-rotation-supported hypermassive NS lasting ∼30−300 ms before collapsing to a BH (hereafter HMNS; Sun et al. 2017; Margalit & Metzger 2017; Pooley et al. 2018; Rezzolla et al. 2018; Ruiz et al. 2018; Margalit & Metzger 2019); a rigid-rotation-supported supra-massive neutron star (with a low dipolar magnetic field; hereafter SMNS) that lasts for tens of seconds to ≳104 s (Margalit & Metzger 2019) before collapsing to a BH (Sun et al. 2017; Ai et al. 2018); and a stable NS that persists.

The outcome after the BNS merger depends on the total system mass (typically ≈2.5−2.7 M among Galactic BNS; Kiziltan et al. 2013; Martinez et al. 2015), and the unknown NS equation of state (EoS; Sun et al. 2017). Indeed, recently, some massive NSs have been identified, for instance, PSR J0952−0607 (MNS = 2.35 ± 0.17 M; Romani et al. 2022) and PSR J0740+6620 (MNS = 2.08 ± 0.07 M; Fonseca et al. 2021), providing the most severe constraints on the dense-matter EoS. The possibility of BNS mergers producing a long-lived NS has been suggested in the literature to interpret some of the X-ray features observed in SGRB afterglows, including X-ray flares (Dai et al. 2006; Gao & Fan 2006), GRB extended emission (Metzger et al. 2008; Gompertz et al. 2014), and the so-called internal plateau1 observed in a good fraction of SGRBs (Rowlinson et al. 2010, 2013; Lü et al. 2015; Gao et al. 2016; Li et al. 2016a,b; Yu et al. 2018; Piro et al. 2019). This interpretation has allowed constraints on the maximum NS remnant mass of ≲2.5−2.8 M (Lasky et al. 2014; Gao et al. 2016; Lü et al. 2015, 2017).

A potentially related phenomenon with BNS merger is extragalactic Fast X-ray transients (FXTs), short flashes of X-ray photons with durations from a few minutes to hours, which have been observed in the ∼0.3−10 keV X-ray band by Chandra, XMM-Newton and Swift-XRT (e.g., Soderberg et al. 2008; Jonker et al. 2013; Glennie et al. 2015; Irwin et al. 2016; Bauer et al. 2017; Lin et al. 2018, 2019, 2020; Lin et al. 2021, 2022; Xue et al. 2019; Alp & Larsson 2020; Novara et al. 2020; Ide et al. 2020; Pastor-Marazuela et al. 2020; Eappachen et al. 2023; Quirola-Vásquez et al. 2022, 2023). At present, their nature remains poorly understood in part due to the lack of EM counterparts at other wavebands. Critically, while 34 extragalactic FXTs have been identified to date, both serendipitously and through careful searches, only in one case, XRT 080109/SN 2008D (Soderberg et al. 2008; Mazzali et al. 2008; Modjaz et al. 2009), has there been a detection of a multi-wavelength counterpart after the outburst. The others have only been identified via archival data-mining from days (e.g., Luo et al. 2014; Lin et al. 2021) to years after the outbursts (e.g., Alp & Larsson 2020; Lin et al. 2020, 2022; Quirola-Vásquez et al. 2022, 2023). Moreover, only a fraction of FXTs has clear host-galaxy associations (e.g., 13 FXTs detected by Chandra have been associated with extended objects; Quirola-Vásquez et al. 2022, 2023; Eappachen et al. 2022, 2023), while far fewer have firm distance and energetic constraints (e.g., Soderberg et al. 2008; Irwin et al. 2016; Bauer et al. 2017; Xue et al. 2019; Alp & Larsson 2020; Novara et al. 2020; Lin et al. 2022; Eappachen et al. 2022, 2023; Quirola-Vásquez et al. 2022, 2023).

A variety of theoretical models have been proposed for the origin of extragalactic FXTs, such as: (i) X-ray emission produced from the shock breakout (SBO; LX, peak ∼ 1042 − 1044 erg s−1) of core-collapse supernova (CC-SN) once it crosses the stellar material (e.g., Soderberg et al. 2008; Nakar & Sari 2010; Waxman & Katz 2017; Novara et al. 2020; Alp & Larsson 2020); (ii) the X-ray emission from the off-axis long GRBs (LX, peak ≲ 1045 erg s−1) is produced by a wider, mildly relativistic cocoon jet once it breaks the surface of a massive progenitor star (Ramirez-Ruiz et al. 2002; Zhang et al. 2004; Nakar 2015; Zhang 2018; D’Elia et al. 2018); (iii) tidal disruption events (TDEs; LX, peak ∼ 1042 − 1050 erg s−1) involving a white dwarf (WD) and an intermediate-mass black hole (IMBH), whereby X-ray emission is generated from the accretion of part of the WD by the IMBH (e.g., Jonker et al. 2013; MacLeod et al. 2014; Saxton et al. 2021; Maguire et al. 2020); and (iv) BNS merger (LX, peak ∼ 1044 − 1051 erg s−1, e.g., Dai et al. 2018; Fong et al. 2015; Sun et al. 2017; Bauer et al. 2017; Xue et al. 2019), where the X-ray photons are created by the accretion of fallback material onto the remnant BH, a wider and mildly relativistic cocoon, or the spin-down magnetar emission (e.g., Metzger & Piro 2014; Sun et al. 2017, 2019; Metzger et al. 2018a).

Twenty-two of the above FXTs come from the analysis of Quirola-Vásquez et al. (2022, 2023; hereafter Papers I and II, respectively) using two decades of Chandra data (Paper I was based on the available Chandra Source Catalog 2.0 covering until the end of 2014, while Paper II carried out a systematic reduction of Chandra data from 2014 to April 2022). Two relatively distinct sets of FXTs were distinguished according to their distances: nearby and distant FXTs, with luminosity distances dL ≲ 100 and dL ≳ 100 Mpc, respectively. Based on host-galaxy associations and estimated redshifts, an SN SBO scenario has been ruled out for many distant FXTs, based on inferred X-ray luminosity peak and isotropic energy arguments (Eappachen et al. 2022, 2023; Quirola-Vásquez et al. 2022, 2023). The energetics and host-galaxy properties of the most distant FXTs remain consistent with an origin as a GRB or an IMBH-WD TDE.

It has been predicted that some FXTs could be EM counterparts of BNS mergers (Zhang 2013; Metzger & Piro 2014; Sun et al. 2017, 2019). Luo et al. (2014) and Zheng et al. (2017) identified two new unusual FXTs in the 7 Ms Chandra Deep Field-South (CDF-S) dataset, XRT 141001 and XRT 150322, denoted as “CDF-S XT1” and “CDF-S XT2”, respectively. These two FXTs were later studied in detail by Bauer et al. (2017) and Xue et al. (2019), respectively, and re-identified as FXTs 14 and 16 in Papers I and II, respectively. XRT 150322 is particularly intriguing because it exhibits a flat, extended X-ray light curve that suggests a magnetar-wind origin (Sun et al. 2019; Xiao et al. 2019; Lü et al. 2019), similar to some SGRB X-ray afterglows (e.g., Rowlinson et al. 2010, 2013; Troja et al. 2019). Its host-galaxy properties are consistent with the observed host galaxies of SGRBs (Xue et al. 2019). XRT 141001, on the other hand, is consistent with an “orphan” X-ray afterglow from an off-axis SGRB with weak optical emission (Bauer et al. 2017; Sarin et al. 2021). Sun et al. (2019) proposed a unified model to interpret both transients within the framework of the BNS merger magnetar model by considering different observer viewing angles (following on from earlier models developed by Yu et al. 2013; Sun et al. 2017). According to this new work, CDF-S XT2/FXT 16 is observed from a region where the magnetar spin-down powered X-ray emission escapes freely, whereas CDF-S XT1/FXT 14 originates from a region that is initially opaque to the X-ray emission but becomes optically thin after the ejecta expands and becomes ionised.

In this paper, we explore a possible association of FXTs with massive, rapidly spinning magnetars produced by a BNS merger observed at different viewing angles. The manuscript is organized as follows: we explain the magnetar model in Sect. 2; we describe the data considered in this work and the fitting process in Sects. 3 and 4, respectively; Sects. 57 present the light curve fitting results, the upper limit constraints using optical data, and the potential for GW detections, respectively; and finally, Sect. 8 presents a discussion and conclusions. Throughout the paper, a concordance cosmology with parameters H0 = 70 km s−1 Mpc−1, ΩM = 0.30, and ΩΛ = 0.70 is adopted. Magnitudes are quoted in the AB system.

2. Magnetar model in a nutshell

Some GRBs might have been with millisecond magnetars (see for instance, Rowlinson et al. 2010, 2013; Gompertz et al. 2014; Dainotti et al. 2017). The BNS merger model acknowledges the potential rapid formation of a BH, whose accretion could power a SGRB and its subsequent afterglow. However, the discovery of NSs with masses around ≈2 M (Demorest et al. 2010; Antoniadis et al. 2013; Romani et al. 2022) implies that certain mergers might result in the creation of a transient or stable rapidly-spinning magnetar instead (Duncan & Thompson 1992; Nousek et al. 2006; Zhang et al. 2006; Metzger et al. 2008). As such, once the BNS merger occurs, a magnetar may be formed with enough rotational energy to prevent gravitational collapse and the direct formation of a BH (Zhang & Mészáros 2001). Its rotational energy could be lost via GWs (via the quadrupole moment in the mass distribution) and EM radiation (via the magnetic field dipole distribution), causing the magnetar to spin down. At a certain moment, it may reach a critical point in time (from tens of seconds to > 104 s) at which centrifugal forces can no longer support the degeneracy pressure against gravity, resulting in the formation of a BH (Siegel & Ciolfi 2016a,b; Sun et al. 2017). This scenario can lead to a brief X-ray plateau, as observed in certain SGRB afterglows, resulting from the injection of the spin-down energy of the magnetar and the subsequent collapse into a BH (Dai & Lu 1998; Fan & Xu 2006; Rowlinson et al. 2010, 2013; Lasky et al. 2017). An alternative explanation involving fall-back accretion has also been proposed (Rosswog 2007).

Here, we consider the formalism from Yu et al. (2013) and Sun et al. (2017) to describe the physical model of the magnetar emission. Approximately isotropic X-ray emission may be produced by the internal dissipation of the magnetar wind (i.e., covering a wide solid angle); however, this process is not well understood (e.g., Metzger & Piro 2014). Their luminosity tracks the dipole spin-down luminosity of the magnetar with a certain efficiency. The wind energy injection dissipated from the magnetar would produce a plateau in the X-ray light curve (Zhang & Mészáros 2001) followed by a power-law decay (like FX ∝ tα where α ∼ −1 to −2) due to the magnetar spin down.

We focus on the scenario whereby the BNS merger produces a stable or supra-massive millisecond magnetar, which can be characterized by the initial rotational period (Pi) and magnetic field. In addition to the formation of the magnetar, a fraction of unbound mass is ejected (hereafter, ejecta material, Mej) due to processes that occur on a dynamical timescale (see, e.g., Oechslin et al. 2007; Bauswein et al. 2013; Hotokezaka et al. 2013), and depend on the total binary mass, mass ratio, and EoS (see Metzger 2019, and references therein). This ejecta material, rich in r-process elements, powers a thermal transient called a kilonova (KN; Li & Paczyński 1998). The total dynamical ejecta mass typically ranges from 10−4 − 10−2M for BNS mergers (e.g., Kyutoku et al. 2013, 2015; Foucart et al. 2017) and moves outward with velocities of 0.1−0.3c (some authors propose that outflow velocities greater than ∼0.6c can occur, which through the decay of free neutrons in the outermost layers of the ejecta or prompt shock heating of the ejecta by a relativistic outflow can generate UV precursors of KNe; Dean et al. 2021; Metzger 2019).

Once the merger occurs, a strongly collimated, relativistic jet might be emitted perpendicular to the orbital plane of the system (e.g., Narayan et al. 1992; Metzger et al. 2008; Zhang & Dai 2010; Bucciantini et al. 2012; Berger 2014), and three distinct geometrical zones can be defined (see Fig. 1). The magnetar emission and the ejecta material produce a configuration where the observed phenomenon depends on the observer viewing angle relative to the outflow or jet axis (see Yuan et al. 2015; Sun et al. 2017):

thumbnail Fig. 1.

Cartoon illustrating the jet, free, and trapped zones, according to the magnetar model for FXTs. Left panel: situation at time t < tτ. Right panel: situation at time t > tτ. The insets in both panels show the light curves expected for the free and trapped zones before and after tτ.

1. The jet zone is the on-axis SGRB jet direction. In this direction, the X-ray emission from the magnetar can be observed as a plateau in the SGRB afterglow (e.g., Rowlinson et al. 2010, 2013; Lü et al. 2015; Gompertz et al. 2014).

2. The free zone includes any directions where the ejecta material does not obscure the magnetar emission region from the observer’s line of sight. For instance, along slightly off-axis directions with respect to the jet, resulting in no strong SGRB signal, although the X-ray magnetar emission can escape freely. There could still be a weak (GRB 170817A-like) GRB along this viewing direction, which clears a funnel (such a configuration could be produced by a structured jet-cocoon geometry; Lazzati et al. 2017; Lamb & Kobayashi 2017) to allow X-rays to escape, but such an SGRB-like GRB 170817A would not be detectable with current γ-ray detectors beyond ≳80 Mpc (Zhang et al. 2018; Sun et al. 2019).

3. The trapped zone occurs from directions where the dynamical ejecta obscures the magnetar X-ray emission. In this case, the ejecta material is heated and accelerated (via the Poynting flux) by the magnetar emission, which eventually causes a magnetar-enhanced KN (e.g., Yu et al. 2013; Sun et al. 2017; Metzger & Piro 2014; Metzger 2019, and references therein). In this case, the time when the ejecta transforms from optically thick to thin (i.e., when t = tτ; where tτ denotes the time when the optical depth of the ejecta material is one) plays an important role in the dynamics of the ejecta and eventually obscured the X-ray emission.

Below, in Sects. 2.1 and 2.2, we explain the mechanisms behind the free and trapped zones, respectively. In Appendix A, one can find the mathematical details of this model.

2.1. Free-zone emission

During the spin-down process, a continuous isotropic Poynting-flux-dominated outflow is launched, which could potentially inject energy and modify the ejecta material, and the afterglow properties (Sun et al. 2017; Metzger et al. 2018b; Ai et al. 2018; Yu et al. 2018; Piro et al. 2019; Troja et al. 2020). The millisecond magnetar loses its rotational energy (Erot, see Eq. (A.1)) through EM and GW radiation (from non-axisymmetric deformation Shapiro & Teukolsky 1983; Usov 1992; Zhang & Mészáros 2001; Gao et al. 2016; Lasky & Glampedakis 2016; Sun et al. 2017). The evolution of the NS angular frequency considers the EM and GW losses (see Eq. (A.2)), which both play a role in the spin-down evolution, and we define the dipole spin-down luminosity as (Zhang & Mészáros 2001; Metzger & Piro 2014; Lasky & Glampedakis 2016; Sun et al. 2017; Zhang 2018):

L sd ( t ) = B P 2 R M 6 Ω ( t ) 4 6 c 3 , $$ \begin{aligned} L_{\rm sd}(t)=\frac{B_{\rm P}^2R_{\rm M}^6\Omega (t)^4}{6c^3}, \end{aligned} $$(1)

where Ω(t) is the angular frequency evolution (Ω(t) = 2π/P(t), and P(t) is the rotational period), BP is the dipole component of the magnetic field at the poles2, RM is the magnetar radius, and c is the speed of light. Importantly, a strong magnetic field provides a mechanism for extracting rotational energy from the magnetar via EM spin-down. Indeed, magneto-hydrodynamic simulations show that the magnetic fields in a BNS merger are amplified to values exceeding the field of Galactic magnetars (i.e., BP ∼ 1015 − 1016 G; Price & Rosswog 2006; Zrake & MacFadyen 2013; Kiuchi et al. 2014).

The internal dissipation of the magnetar emission is not perfectly efficient, that is, the magnetar wind emits an X-ray luminosity that is only a fraction of Lsd (labeled as η), converting the dipole spin-down luminosity to the observed X-ray luminosity (ηLsd). In the free zone, without obscuring material, the magnetar wind emission should be

L X , free ( t ) = η L sd = η B p 2 R M 6 Ω 4 ( t ) 6 c 3 · $$ \begin{aligned} L_{\rm X,free}(t)=\eta L_{\rm sd}=\frac{\eta B_{\rm p}^2R_{\rm M}^6\Omega ^4(t)}{6c^3}\cdot \end{aligned} $$(2)

Typically, η is assumed to take values of ≈10−3 and to be constant in time (Sun et al. 2019). According to the geometrical/orientation-based structure (see Fig. 1), the free-zone X-ray emission will be visible regardless of the evolution and dynamics of the ejecta material.

Figure 2 shows the numerical solution for the free zone (using Eqs. (1), (2) and (A.2), and under the assumption of η = 10−3) exploring the parameter space, that is, varying the magnetic field component (BP), the initial rotational period (Pi), magnetar radius (RM), and the ellipticity of the magnetar (ϵ). If the magnetic field or magnetar radius increases, the plateau luminosity is boosted, but its duration declines (see Fig. 2, left panels). Meanwhile, as the initial period increases, the plateau luminosity decreases (see Fig. 2, top-right panel) while the duration of the plateau grows. All these light curves mentioned above follow a decay as LX ∝ t−2 after the plateau, that is, the EM losses dominate the system (see Eqs. (A.4) and (A.6) for more details). Figure 2, bottom-right panel, shows light curves with different ϵ values. As ϵ increases (i.e., more significant magnetar deformation), the GWs losses become more dominant (i.e., ϵ ≳ 10−3; see Eqs. (A.5) and (A.7) for more details), and the light curve decays as LX ∝ t−1 (see Eq. (A.7)), followed by a regime wherein EM losses dominate as LX ∝ t−2 (see Eq. (A.8)).

thumbnail Fig. 2.

Theoretical light curves in the free-zone using Eqs. (1), (2) and (A.2). We explore how the light curve changes as we vary the magnetic field (top-left panel), the initial spin period (top-right panel), the magnetar radius (bottom-left panel), and the ellipticity (bottom-right panel). The text in the upper right of each panel shows the fixed parameters in the numerical solutions, while the legends and colored curves indicate the varied parameters and the range in the values.

2.2. Trapped-zone emission

As we explained above, in the trapped zone, part of the X-ray magnetar emission (Eq. (1)) is trapped behind the ejecta material (see Fig. 1; Yu et al. 2013; Sun et al. 2017). The trapped X-ray radiation heats and accelerates the ejecta material via pdV work. Under these circumstances, the optical depth (τ) of the ejecta is defined as

τ = κ ( M ej V ) ( R ej Γ ) , $$ \begin{aligned} \tau =\kappa \left(\frac{M_{\rm ej}}{V^{\prime }}\right)\left(\frac{R_{\rm ej}}{\Gamma }\right), \end{aligned} $$(3)

where κ and V′ are the opacity and the co-moving volume3, respectively, while Γ and Rej are the Lorentz factor and the radius of the ejecta material, respectively. Initially, the ejecta remains extremely hot and optically thick, that is, τ ≫ 1, immediately after being ejected from the vicinity of the merger. This means that the ejecta material obscures the X-ray emission of the magnetar. At the same time, the magnetar enhances the KN emission, that is, the KN is not powered by radioactive decay only but also by a fraction of the magnetar spin-down energy, rendering it brighter than a KN powered by radioactive decay alone (Yu et al. 2013; Metzger & Piro 2014; Metzger 2019). The KN emission boosted by the X-ray wind from magnetar spin-down is called a merger-nova (a term adopted by Yu et al. 2013). The peak luminosities of the merger-nova candidates were estimated to be above ≳1042 erg s−1, that is, one order of magnitude brighter than a standard r-process radioactivity-powered KN (i.e., ejecta mass and velocity of ∼10−3 − 10−4M and 0.3 − 0.1c, respectively, and opacity between ∼1.0 and 10 cm2 g−1). The merger-nova phenomenon is predicted to cover a wide range of peak luminosities depending on the magnetar parameters (Zhang 2018; Yu et al. 2013). Once the trapped magnetar wind becomes optically thin, the magnetar X-ray emission should rise quickly to a level similar to the free-zone luminosity at that point in time, and the photons can escape freely without being reprocessed by the ejecta (under the condition τ < 1).

Observationally, some KN candidates associated with SGRBs have been reported in the literature, for example, GRB 130603B (Tanvir et al. 2013; Berger et al. 2013), GRB 060614 (Yang et al. 2015), GRB 050709 (Jin et al. 2016) and GRB 191019A (Levan et al. 2023), while a systematic search for KNe in GRB afterglows was done by Rossi et al. (2020); also, three compelling detections have been found: GRB 170817A (e.g., Abbott et al. 2017c,d), GRB 211211A (Yang et al. 2022; Rastinejad et al. 2022; Troja et al. 2022), and recently GRB 230307A (Levan et al. 2024). Additional systematic campaigns to search for merger-nova events have revealed three more candidates: GRBs 050724, 070714B, and 061006 (Gao et al. 2017). In the case of the GRB 170817A, very late-time observations with the Chandra X-ray Observatory show an unabsorbed X-ray flux (luminosity) of ≈1.4 × 10−15 erg cm−2 s−1 (≈2.7 × 1038 erg s−1), which is higher than expected from simple structured-jet model predictions, and could be explained by the energy injection from a long-lived central engine (although other scenarios are not discarded completely, e.g., Troja et al. 2020; Hajela et al. 2022). Sun et al. (2017) and Yu et al. (2013) have modeled the merger-nova complex system by describing the dynamical evolution of the ejecta and how the heating and cooling processes affect it, which we summarize below (for more details, see Appendix A).

From the total energy of the ejecta (Eej, see Eq. (A.10)), excluding the rest-mass energy, and considering the luminosities which boosted or lost energy to the ejecta material (see Eq. (A.11)), such as the luminosity of the magnetar (Lsd), the radioactive decay of the r-process material (Lra), and the bolometric luminosity of the heated electrons (Le) in the system, we can derive the evolution of the Lorentz factor (Γ) of the ejecta as a function of the observer time as

d Γ d t = L sd + L ra L e Γ D ( d E int / d t ) M ej c 2 + E int , $$ \begin{aligned} \frac{\mathrm{d}\Gamma }{\mathrm{d}t}=\frac{L_{\rm sd}+L_{\rm ra}-L_{\rm e}-\Gamma \mathcal{D} (\mathrm{d}E_{\rm int}^{\prime }/\mathrm{d}t^{\prime })}{M_{\rm ej}c^2+E_{\rm int}^{\prime }}, \end{aligned} $$(4)

where E int $ E^\prime_{\mathrm{int}} $ is the internal energy in the co-moving frame, and 𝒟 = 1/[Γ(1 − β cos θ)] is the Doppler factor with β = 1 Γ 2 $ \beta = \sqrt{1-\Gamma^{-2}} $, θ = 0° for an on-beam observer. The internal energy evolution should include heating from the magnetar emission and radioactivity and cooling via the electron radiation and the adiabatic work pdV (Kasen & Bildsten 2010). The evolution in the co-moving frame of the internal energy can be expressed as

d E int d t = ξ L sd + L ra L e P d V d t , $$ \begin{aligned} \frac{\mathrm{d}E_{\rm int}^{\prime }}{\mathrm{d}t^{\prime }}=\xi L_{\rm sd}^{\prime }+L_{\rm ra}^{\prime }-L_{\rm e}^{\prime }-P^{\prime }\frac{\mathrm{d}V^{\prime }}{\mathrm{d}t^{\prime }}, \end{aligned} $$(5)

where the efficiency ξ is defined as the fraction of the magnetar spin-down luminosity used to heat the ejecta, and P′ is the pressure (dominated by radiation) in the co-moving frame (see Eq. (A.14)). The evolution of the internal energy also depends on the co-moving volume evolution (see Eq. (A.15)). The definitions of Lra and Le are provided in Appendix A (see Eqs. (A.13) and (A.17), respectively).

The spectrum of the merger nova should approximately resemble a blackbody under the co-moving temperature, T′,

k b T = { k b ( E int a V τ ) 1 / 4 for τ > 1 k b ( E int a V ) 1 / 4 for τ 1 , $$ \begin{aligned} k_{\rm b}T^{\prime }=\left\{ \begin{matrix} k_{\rm b}\left(\frac{E_{\rm int}^{\prime }}{aV^{\prime }\tau }\right)^{1/4}&\mathrm{for}\ \tau > 1 \\&\\ k_{\rm b}\left(\frac{E_{\rm int}^{\prime }}{aV^{\prime }}\right)^{1/4}&\mathrm{for}\ \tau \lesssim 1 \\ \end{matrix} \right., \end{aligned} $$(6)

where kb and a are the Boltzmann and the blackbody radiation constants, respectively.

For blackbody emission including the co-moving temperature T′, the spectral luminosity at a particular frequency, ν, is given by

( ν L ν ) bb = 8 π 2 D 2 R ej 2 h 3 c 2 ( h ν / D ) 4 exp ( h ν / D k b T ) 1 , $$ \begin{aligned} (\nu L_\nu )_{\rm bb}=\frac{8\pi ^2\mathcal{D} ^2R_{\rm ej}^2}{h^3c^2}\frac{(h\nu /\mathcal{D} )^4}{\exp {(h\nu /\mathcal{D} k_{\rm b}T^{\prime })}-1}, \end{aligned} $$(7)

where h is Planck’s constant. For the trapped-zone case, the observed X-ray luminosity can be written as

L X , trapped = e τ η B p 2 R M 6 Ω ( t ) 4 6 c 3 + ( ν X L ν , X ) bb , $$ \begin{aligned} L_{\rm X,trapped} = e^{-\tau }\frac{\eta B_{\rm p}^2R_{\rm M}^6\Omega (t)^4}{6c^3}+(\nu _{\rm X} L_{\rm \nu ,X})_{\rm bb}, \end{aligned} $$(8)

where the first term comes from the dissipating wind, which is negligible when τ ≫ 1, while the second term is the Wien tail emission of the blackbody merger-nova photosphere calculated from Eq. (7) at X-ray frequencies, νX.

Meanwhile, the specific flux of the merger nova, in the observer frame, is

F ν = 1 4 π d L 2 ( 1 + z ) 8 π 2 D 2 R ej 2 h 3 c 2 ν ( 1 + z ) [ h ν ( 1 + z ) / D ] 4 exp [ h ν ( 1 + z ) / D k b T ] 1 , $$ \begin{aligned} F_\nu =\frac{1}{4\pi d_{\rm L}^2(1+z)}\frac{8\pi ^2\mathcal{D} ^2R_{\rm ej}^2}{h^3c^2\nu (1+z)}\frac{[h\nu (1+z)/\mathcal{D} ]^4}{\exp {[h\nu (1+z)/\mathcal{D} k_{\rm b}T^{\prime }]}-1}, \end{aligned} $$(9)

where dL is the luminosity distance, z is the redshift, and ν(1 + z) is the frequency in the source frame. It is clear that the presence of ejecta material should obscure part of the magnetar wind emission, where the shape of the X-ray light curve should depend on tτ. Nevertheless, events exhibiting a plateau at early times could still be interpreted from a trapped-zone emission perspective.

We analyze different cases for comparison with the free-zone emission, where the plateau has a duration of TX and the trapped-zone emission has an obscuration time tτ:

1. Cases with tτ ≪ TX could be related to sources with a very low ejecta mass (i.e., the line of sight passes through the trapped zone where the ejecta masses are still low), such that the resulting light curves are similar to the free-zone case.

2. For tτ ≫ TX (i.e., the line of sight passes deep into the trapped zone), the resulting light curve (and spectra) will strongly differ from the free zone case, and we expect to observe just (part of) the decay phase.

3. If tτ ∼ TX, we may observe just a small final portion of the plateau emission, followed by the decay phase.

If the X-ray emission is viewed from the free zone, the X-ray plateau approximately begins with the BNS merger coalescence. Then, the multi-wavelength emission of the merger-nova at t < tX, e, where tX, e is the emerging time of the X-rays, is not important because the magnetar emission and the merger nova should be approximately simultaneous. On the other hand, if the X-ray emission is observed from the trapped zone, the merger time should be tX, e − tτ, that is, the multi-wavelength observations of the merger nova at t < tX, e become relevant (Ai & Zhang 2021).

Numerical solutions to the above equations for the free- and trapped-zone emissions are shown in Fig. 3 for a set of ejecta masses from 2 × 10−4 to 5 × 10−3M (color coded), and two different opacities κ = 1.0 (solid lines) and 10 cm2 g−1 (dashed lines) under some fixed initial conditions (see the caption of Fig. 3). Specifically, Fig. 3, left panel, shows the dynamical evolution of the ejecta material (i.e., the Lorentz factor) for different masses. The Lorentz factor increases substantially while t ≲ tτ, beyond which point the rate slows and flattens near some maximum, with a bit of an overshoot observed for higher ejecta masses. The slope of the rise in Γ and the time to reach maximum are both strong functions of the ejecta mass (with lower ejecta masses associated with more highly relativistic ejecta). These trends are expected because the ejecta becomes optically thin, and the magnetar emission does not influence the dynamical evolution at t > tτ. Figure 3, middle panel, depicts the optical depth evolution, which is high at early times and shows a smooth, strong decline with time as the ejecta expands, eventually passing from the optically thick to thin regimes around t ≈ tτ. This transition is an important factor for observing the magnetar emission. In Fig. 3 right panel, we see that the X-ray luminosity of the magnetar exhibits an extremely steep rise at early times, reaching a peak luminosity at t ∼ tτ, after which the ejecta is optically thin, allowing the magnetar X-ray photons to escape freely. In Sects. 2.2.1 and 2.2.2, we discuss the role of opacity in the X-ray light curve and the merger-nova emission, respectively.

thumbnail Fig. 3.

Theoretical results of the trapped-zone emission. Left panel: evolution of Lorentz factor (Γ). Middle panel: logarithmic optical depth (τ). The black horizontal line represents the optical depth limit at τ = 1. Right panel: X-ray luminosity (LX) of the magnetar model (Eq. (8), first term). The solid black line shows the free zone spin-down luminosity (ηLsd). The color code shows the different ejecta masses (Mej = 2 × 10−4 − 5 × 10−3M), while the solid and dashed lines represent opacities of 1.0 and 10.0 cm2 g−1, respectively, expected for the blue and red kilonova (Tanaka et al. 2017; Metzger 2019; Tanaka et al. 2020). The initial conditional assumed are: the initial velocity and internal ejecta energy of the ejecta are βi = 0.1, and E int,i E k =0.5 M ej β i 2 c 2 $ E_{\rm int,i} \simeq E_{\rm k} = 0.5\,M_{\rm ej}\beta_{\rm i}^2c^2 $, respectively. The magnetar parameters are: Bp = 5 × 1015 G, Pi = 10−3 s, I = 1045 g cm2, ϵ = 10−5, and ξ = 0.3.

2.2.1. Opacity

There are several uncertainties associated with the expected opacity adopted in theoretical models and from observations of AT2017gfo (e.g., Tanaka et al. 2017, 2020; Banerjee et al. 2020; Sarin et al. 2021). For instance, Banerjee et al. (2020), from simulations, argued that the Planck mean opacity for lanthanide-free ejecta (so-called blue KN) could be in the range κ ≈ 1 − 0.5 cm2 g−1 at t ∼ 0.1 days, and could increase up to κ ≈ 10 − 5 cm2 g−1 at t ∼ 1 day. In other words, the opacity of the ejecta changes with time. Observationally, lanthanide-free ejecta material (κ ≈ 5 − 3 cm2 g−1) is needed to explain the early blue KN emission from AT2017gfo, while lanthanide-rich material (κ ≈ 30 − 20 cm2 g−1) is necessary to reproduce the more long-lasting NIR emission of AT2017gfo (e.g., Tanaka et al. 2020). Thus, currently, it is accepted that AT2017gfo had at least two components: a red and a blue KN (Arcavi et al. 2017; Cowperthwaite et al. 2017; Chornock et al. 2017; Drout et al. 2017; Evans et al. 2017; Nicholl et al. 2017; Smartt et al. 2017; Kasen et al. 2017; Kilpatrick et al. 2017; Shappee et al. 2017; Tanvir et al. 2017; Gao et al. 2017). Nevertheless, some authors (such as Villar et al. 2017, from semi-analytical models) suggested that the data can be well modeled by including a three-component KN with an intermediate (“purple”) opacity component (κ = 3 cm2 g−1).

Figure 3 explores the influence of the opacity on the dynamical evolution of the ejecta and the X-ray light curves considering opacities of κ = 1.0 (solid lines) and 10.0 cm2 g−1 (dashed lines). We see that an increase in the opacity leads to a more extended overshoot in the Lorentz factor (see Fig. 3, left panel), which is more noticeable for higher ejecta masses and longer timescales for the ejecta to become optically thin (middle panel) and reach a lower X-ray peak luminosity (right panel).

2.2.2. Merger-nova emission

Even if only a modest fraction of the rotational energy liberated in the spin down (see Eq. (A.1)) is converted to EM radiation in the hours to years after the merger, this would substantially enhance the EM luminosity of the merger counterpart (Metzger 2019, and references therein). A merger-nova should emit quasi-blackbody-like emission at a specific temperature (given by Eq. (6)), generating emission (given by Eqs. (7) and (9)) across the EM spectrum. The light curves at X-ray, UV, and optical wavelengths (i.e., at 1 keV, 30 eV, and 1 eV energies, respectively) and the temperature evolution of the merger-nova are shown in Fig. 4. The color code represents different ejecta masses (from 3 × 10−4 to 1 × 10−2M), assuming typical parameters (see caption in Fig. 4).

thumbnail Fig. 4.

Evolution of the thermal component and light curves of the merger-nova emission (Eq. (7)) for different ejecta masses, color-coded from 10−4 to 5 × 10−3M. The top-left panel is the temperature evolution in the observer frame, while the other three panels represent the merger-nova emission at different frequencies (i.e., : 1 keV, 30 eV, and 1 eV). The model fixed parameters are ξ = 0.3, η = 10−3, κ = 1.0 cm2 g−1, Bp = 3 × 1015 G, and Pi = 5 ms.

In general, the temperature of the merger-nova decreases with time (see Fig. 4, top-left panel), from ≈106 (at ∼1 s) to ≈103 K (at ∼106 s), causing the peak of the blackbody emission to shift from X-ray to NIR wavelengths over this time period. The temperature decrease is initially rather shallow up to the pivot point, which coincides with the time when the ejecta changes from optically thick to thin, followed by a faster decay. The effect of the magnetar on the enhancement of the temperature is stronger for low-ejecta masses (see Fig. 4, top-left panel).

In detail, the 1 keV X-ray emission of the merger nova dominates at early times, with a peak luminosity of LX, peak ≲ 1037 erg s−1 at 1 s (see Fig. 4, bottom-right panel). The relatively low (< 106 K) and decreasing temperature as the ejecta cools leads to a strong temporal decline as the Wien tail of the blackbody emission passes from the X-ray to UV to optical bands. The blackbody merger-nova luminosity is highest in the UV band (30 eV), reaching peak luminosities of ≈ 1045 erg s−1 at ≈104 − 105 s (see Fig. 4, bottom-left panel). The optical light curves (1 eV) reach peak luminosities of ≈1043 − 1044 erg s−1 on timescales of 4−5 days (see Fig. 4, top-right panel).

Figure 5 shows the optical/NIR merger-nova light curves for the u- to K-bands (derived using Eq. (7)), considering a magnetic field strength and initial rotational period of Bp = 1015 G and Pi = 5 ms, respectively, and ejecta masses of Mej = 10−4 (solid lines) and Mej = 10−2M (dashed lines). KN light curves with no magnetar spin-down enhancement are shown for comparison (dotted lines). Both low and high values for the ejecta opacity of κ = 0.1 (left panels) and κ = 20 cm2 g−1 (right panels), representing lanthanide-free and lanthanide-rich ejecta scenarios, respectively, are shown.

thumbnail Fig. 5.

Theoretical light curves of the merger-nova emission, across several common optical/NIR filters, considering opacities of κ = 0.1 (left panels) and 10 cm2 g−1 (right panels), representing the lanthanide-free and -rich ejecta cases, respectively. KN emission without any magnetar spin-down enhancement is shown for comparison (dotted lines). The model fixed parameters are Pi = 5 ms, ξ = 0.3, and η = 10−3, along with two adopted ejecta masses of Mej = 10−4 (solid lines) and 10−2M (dashed lines). The upper and lower panels adopt magnetic field strengths of Bp = 1015 G and 1016 G, respectively.

We see that the emission peaks earlier in the bluer filters, with little dependence on ejecta mass (see Fig. 5), similar to AT2017gfo (e.g., Kasen et al. 2017; Villar et al. 2017). For the low-opacity (lanthanide-free) ejecta case, the time to peak only changes by a factor of ∼2−3 between the low (Mej ∼ 10−4M) to high (Mej ∼ 10−2M) ejecta masses, with bluer filters showing slightly smaller shifts. Meanwhile, the luminosities decrease in all filters by a factor of ∼5 for a 100-fold increase in ejecta mass. By contrast, for the higher opacity (lanthanide-rich) case, the time to peak is extended by a factor of ∼10 for the bluest and ∼3 for the reddest filters for the higher ejecta mass scenario (see Fig. 5, right panels), while the peak luminosities do not change considerably.

The lower panels of Fig. 5 consider a higher magnetic field of Bp = 1016 G. An enhanced magnetic field implies a lower EM characteristic spin-down timescale (see Eq. (A.4)), and an increased luminosity (see Eq. (A.6)), as the magnetar should lose rotational energy more quickly, and its contribution to the ejecta emission is enhanced in early epochs, although it is brief. This is shown in Fig. 5, where the merger-nova emission is enhanced at early epochs in blue and red filters. In contrast, the magnetar does not give energy to the ejecta at later epochs, and the emission becomes fainter and decays faster.

Finally, according to Metzger (2019), the contribution of the magnetar emission could have three different effects on the light curves of normal KNe, which are visible in Fig. 5: (i) an increase in the peak luminosity by a few orders of magnitude due to additional heating from the magnetar spin-down; (ii) more rapid evolution, that is, an earlier peak time because the greater ejecta velocity due to the kinetic energy added to the ejecta by the spin-down energy released during the optically thick phase; and (iii) emission shifts toward bluer energy bands because the high luminosity of the transient increases its effective temperature. These effects only consider the optimistic case where the magnetar resists an immediate collapse to form a BH or an essential fraction of its rotational energy does not escape as gamma-rays (e.g., Dall’Osso et al. 2009; Corsi & Mészáros 2009).

3. Fast X-ray transient sample

We consider now the sample of 22 FXT candidates reported in Papers I and II, which can be crudely separated according to their distances: (i) 5 nearby FXTs located at ≲100 Mpc, and (ii) 17 distant FXTs located at ≳100 Mpc. The five nearby FXTs have peak luminosities of LX, peak ≲ 1040 erg s−1. Among the 17 distant FXTs, 8 have been associated with faint spatially resolved galaxies (with spectroscopic or photometric redshifts in the range of ∼0.7−2.2), implying peak luminosities of LX, peak ≳ 1044 erg s−1, and a remaining 9 have no clear host galaxies identified so far.

In this work, we do not consider the nearby sample because their low luminosity peak (Lpeak ≲ 1040 erg s−1) does not match with the BNS emission. In particular, we explore for a subset of the distant FXTs whether they can be described by the spin-down power of a massive magnetar formed in the immediate aftermath of a BNS merger, as modeled by the equations in Sect. 2. Notably, the light curves in the free and trapped zones (see Figs. 2 and 3, respectively) are different, such that some distant FXTs may remain more consistent with free zone viewing-angle models, whereas others may be more consistent with trapped zone viewing-angle models. Among the 17 distant FXTs identified in Chandra data, six (FXTs 7, 10, 12, 16, 19, and 22) show early plateaus in their light curves (potential free-zone cases) while three sources (FXTs 14, 15, and 20) show clear rising and falling behavior (potential trapped-zone cases). The other eight Chandra FXTs do not follow the expected light curves4 of this interpretation, which suggests a different origin, for instance, IMBH TDEs.

Before exploring the parameter space that explains the FXT light curves, we briefly review their current distances. FXT 16 has a firm spectroscopic host redshift of z = 0.738 (Xue et al. 2019), while FXTs 14 and 19 have constrained photometric host redshifts of 2 . 23 1.84 + 0.98 $ 2.23_{-1.84}^{+0.98} $ and 1.44 ± 0.08, respectively (Bauer et al. 2017; Lin et al. 2022; Quirola-Vásquez et al. 2023). FXT 22 has three possible host galaxy associations, denoted as cX, cNE and cW (at 0 . 6 $ 0{{\overset{\prime\prime}{.}}}6 $, 4 . 2 $ 4{{\overset{\prime\prime}{.}}}2 $ and 3 . 6 $ 3{{\overset{\prime\prime}{.}}}6 $ from the transient position, respectively; Eappachen et al. 2023); unfortunately, cX (which is formally inside the X-ray uncertainty region) has only a single detection at g = 25.9 ± 0.1 AB mag, precluding any photometric and spectroscopic redshift determination, while, cNE and cW have spectroscopic and photometric redshifts of zcNE = 1.5082 ± 0.0001 (Andreoni et al. 2021a; Jonker et al. 2021) and z cW = 1 . 04 0.14 + 0.22 $ z_{\mathrm{cW}} = 1.04^{+0.22}_{-0.14} $ (Eappachen et al. 2023), respectively. For FXTs for which no clear host galaxy has been identified so far, we adopt two nominal redshifts of 0.5 and 1.0 (see Papers I and II for details).

4. Fitting data

We fit the observed light curves using the free and trapped models to derive the magnetar and ejecta parameters. For the fitting process, we assume typical values, based on simulations and previous observational constraints, of RM = 1.2 × 106 cm, and η = 10−3 (e.g., Sun et al. 2019, reported values of η ∼ 10−3 for the magnetar emission of CDF-S XT1 and XT25), and the initial velocity, radius and internal energy of the ejecta of βi = 0.1 (i.e., initial Lorentz factor of Γi = 1.005), Rej, i = 2.4 × 107 cm (e.g., Radice et al. 2018a,b), and Eint, i ≃ Ek = 0.5 Mej(βic)2, respectively, where Ek is the kinetic energy. The assumed parameters do not affect the results of the fitting. For instance, a change in the initial velocity and radius of the ejecta of βi = 0.3 and Rej, i = 4.8 × 107 cm (i.e., an increase by a factor of ≈3 − 2), respectively, means a switch in the peak luminosity of ≲1%. Spherical symmetry amongst the ejecta material is assumed, as we expect the ejecta to expand radially over many orders of magnitude from ∼106 to ∼1015 cm (Metzger 2019). We perform the light-curve fitting with the least squares method, using the Python package LMFIT, while to determine the uncertainties of the parameters, we use a Bayesian inference from a Monte Carlo method (Foreman-Mackey et al. 2013; Sharma 2017) and the MULTINEST nested sampling algorithm (Feroz & Hobson 2008; Feroz et al. 2009). We assume uniform logarithmic prior distributions for each free parameter (see Table 1 for more details).

Table 1.

Parameters and uniform priors distributions considered for the free and trapped magnetar models.

The two magnetar parameters BP and Pi are free parameters in both the free-zone and trapped-zone models, while three ejecta parameters (Mej, κ and ξ) in the trapped zone model are also free in the fitting process. However, due to the low number of bins in the light curve of FXT 15, we freeze some parameters during the fitting process to retain a positive number of degrees of freedom (d.o.f.).

Notably, in the trapped-zone scenario, the peak emission depends strongly on both the Mej and the time when the light curve starts (i.e., the BNS merger). As this is not necessarily known, we introduce an additional parameter called T0, which attempts to correct the start point. In the free and trapped models, we assume an ellipticity parameter of ϵ = 10−5 as a consequence of the fact that the light curve decays are not consistent with FX ∝ t−1 (see Eq. (A.7)).

The best-fitting parameters are shown in Table 2, including 5000 simulations followed to convergence to determine errors. The uncertainties in the parameters represent the 16% and 84% percentiles of the Markov chain Monte Carlo (MCMC) distributions. The χ2, d.o.f., Bayesian information criterion (BIC), and Akaike information criteria (AIC) are derived from the best-fitted parameters under the least squares method. In Sects. 4.1 and 4.2, we describe the results obtained for the free- and trapped-zone samples, respectively.

Table 2.

Parameters obtained from the fitting process of the magnetar model.

4.1. Free-zone FXTs

The light curves of known-redshift FXTs 16, 19, and 22 (assuming either zcW = 1.04 or zcNE = 1.5105) are well-fitted with the free-zone model. The best-fit models and their residuals are shown in Fig. 6, while the posterior distributions can be found as Figs. C.4 and C.5; the latter suggest a correlation between BP and Pi with a coefficient between 0.6 and 0.9. The best-fit magnetar parameters cover a narrow region in the parameter space of BP ≈ 3 × 1015 − 7 × 1014 G and Pi ≈ 7.9 − 1.2 ms, and are consistent with literature-derived values in the cases of FXT 16 (i.e., BP ∼ 6 × 1015 G and Pi ∼ 4.4 ms; Xue et al. 2019; Xiao et al. 2019; Lü et al. 2019) and FXT 22 (i.e., BP ≲ 7 × 1014 G; Ai & Zhang 2021).

thumbnail Fig. 6.

Free-zone fitting (solving the system of Eqs. (2) and (A.2)) in the source rest-frame for FXTs with known and nominal redshifts (see Table 2). Time dilation effects correct the data, and the bottom panels show the ratio between the data and the model (D/M).

For the hostless FXTs 7, 10, and 12, we assume two nominal redshifts of znominal = 0.5 and 1.06. These FXTs are well-fitted by the free-zone model, with best-fitting models and residuals shown in Fig. 6 and posterior distributions in Figs. C.1C.3. For FXTs 7 and 12, Lin et al. (2022) also identified that the EM losses are the main contribution to explaining the light curve trends, discarding GW losses as the prime route to dissipating rotational energy. It is important to realize that increasing the redshift from z = 0.5 to z = 1.0 results in slight drops in both BP and Pi (see Table 2). Thus, if these FXTs lie at larger distances, they will likely have lower magnetar parameters, although these are presumably bounded by the NS breakup limit at ∼0.96 ms (Lattimer & Prakash 2004; Rowlinson et al. 2013).

4.2. Trapped-zone FXTs

The light curve of known-redshift FXT 14 is well-fitted with the trapped-zone model. The best fit and residuals are shown in Fig. 7, while the posterior distributions can be found in Fig. C.6. Notably, the posterior distributions of FXT 14 imply that degeneracies exist between the parameters T0, ξ, and κ. Although the ejecta appears to be low mass (Mej ∼ 2 × 10−4M) and lanthanide-free (κ ∼ 0.1 cm2 g−1), the uncertainties for the opacity are significant (∼6−0.03 cm2 g−1), spanning the values adopted or reported in the literature (i.e., κ ≈ 1 cm2 g−1; Sun et al. 2019). The derived magnetar parameters are also consistent with the values obtained in previous works (e.g., BP ≈ 1016 G and Pi ≈ 1.2 ms; Sun et al. 2019). Additionally, we find that only a small fraction of the magnetar emission is transferred to the ejecta material, that is, ξ ∼ 10−3.

thumbnail Fig. 7.

Trapped zone light-curve fitting of luminosities (solving the system of Eqs. (4)–(8)) in the source rest frame for FXTs with known and nominal redshifts (see Table 2). The bottom panels show the ratio between the data and the model values (D/M). The dashed lines represent the unobscured magnetar emission.

The light curves of FXTs 15 and 20 show similar behavior and are well-fitted by plausible magnetar and ejecta parameters (see Table 2). Figures C.7C.10 present their posterior distributions for nominal redshifts of z = 1.0 and 0.5, while Fig. 7 shows their best fits and residuals. Unfortunately, due to the low number of bins for FXT 20 (d.o.f. = 3), parameters such as κ, Mej, ξ, and Pi are not constrained properly. In particular, we find that the association of FXT 20 with a lanthanide-rich or free ejecta depends on the adopted redshift, that is, under the assumptions of z = 0.5 and 1.0 the opacities are κ < 40 and ∼0.1 cm2 g−1, respectively. Similar to FXT 14, the magnetar luminosity contributes only a small fraction to the heating of the ejecta material (because of the low value of ξ ∼ 10−3). In the case of FXT 15, it has even fewer bins, requiring that during the fitting, we fixed the parameters κ and ξ (see Table 2). We consider only the lanthanide-free case based on the results of FXTs 14 and 20. Notably, FXT 15 has the most massive ejecta material (i.e., Mej ∼ 10−3M), and the highest T0 (∼102 s). Both values are a consequence of the longer time to reach the peak luminosity, around ∼103 s, compared to FXTs 14 and 20 (≲102 s).

Within the magnetar scenario, the early rising phases in FXTs 14, 15, and 20 are interpreted as the emergence of the X-ray magnetar emission wind during its optically thick phase, while the declining phases of the light curves are well-fitted as the spin-down luminosity of the magnetar until the end of the detections. The best-fitted parameters suggest comparable magnetar parameters (magnetic field and initial rotational period values of ∼1016 − 1015 G and ∼7.9 − 1.0 ms, respectively).

Finally, we must acknowledge the limitations of our modeling as far as these may influence interpretation. For instance, one weakness of our model is that it considers only a fixed opacity through the numerical solution. Yet, the only known example of a KN, AT2017gfo, exhibited two components with different opacities or multiple opacity compositions in the ejecta material. Clearly, AT2017gfo is a unique event, while a distribution of event properties from future KN detections is necessary to understand better the physics involved in these sources and the behavior of the opacity during and after the merger.

5. Magnetar’s parameters comparison

Around ≈30 − 50% of SGRBs detected by Swift-XRT have light curves that exhibit extended emission or an X-ray plateau compared to the standard afterglow power-law decay, which typically is interpreted as emission from a magnetar remnant (see, for instance, Rowlinson et al. 2013; Lü et al. 2015). However, after the BNS merger, the young magnetar emission should be essentially isotropic, suggesting that a relatively large population of magnetar-driven FXTs without associated beamed gamma-ray emission should exist (e.g., Zhang 2013). As such, we might expect that the magnetar parameters, BP and Pi, for FXTs should be similar to those derived from the X-ray afterglows of SGRBs.

Figure 8 compares the initial rotational period and magnetic field determined for our sample with long GRBs (LGRBs; Lyons et al. 2010; Yi et al. 2014), SGRBs (Rowlinson et al. 2013; Lü et al. 2015), and super-luminous supernovae (SLSNe; Nicholl et al. 2017). In general, SGRBs have the strongest magnetic fields ( B P SGRBs 8 × 10 16 $ B_{\mathrm{P}}^{\mathrm{SGRBs}} \lesssim 8 \times 10^{16} $ G) and a wide range of initial rotational periods ( P i SGRBs 30 $ P_{\mathrm{i}}^{\mathrm{SGRBs}} \lesssim 30 $ ms); however, the majority of the magnetar parameters of SGRBs are just upper limits (e.g., Lü et al. 2015). On the other hand, SLSNe have the weakest magnetic fields ( B P SLSNe 5 × 10 14 $ B_{\mathrm{P}}^{\mathrm{SLSNe}} \lesssim 5 \times 10^{14} $ G) and a narrow range of initial rotational periods ( P i SLSNe 8 1 $ P_{\mathrm{i}}^{\mathrm{SLSNe}} \sim 8{-}1 $ ms). LGRBs fall in parameter space between SGRBs and SLSNe, with P i LGRBs 20 $ P_{\mathrm{i}}^{\mathrm{LGRBs}} \lesssim 20 $ ms and B P LGRBs 5 × 10 15 5 × 10 14 $ B_{\mathrm{P}}^{\mathrm{LGRBs}} \sim 5 \times 10^{15}{-}5 \times 10^{14} $ G, slightly overlapping with both kinds of objects.

thumbnail Fig. 8.

Comparison of the derived magnetar parameters in this work (BP and Pi) with LGRBs (green triangles; Lyons et al. 2010; Yi et al. 2014), SGRBs (black squares; Rowlinson et al. 2010, 2013; Lü et al. 2015) and SLSNe (blue circles; Nicholl et al. 2017). The gray region represents the breakup limit of a neutron star, which is around 0.96 ms (Lattimer & Prakash 2004). FXTs related to the free and trapped models are represented by star and diamond markers, respectively.

Overall, the derived magnetar parameters from FXTs populate a region between LGRBs and SGRBs, overlapping the upper and lower portions of those populations, respectively. By contrast, the sample of FXTs and SLSNe do not overlap at all (refuting a possible association). This reinforces the results inferred from the host-galaxy properties presented in Papers I and II, and from the derived ejecta masses (Mej ∼ 10−4 − 10−2M), which match well with simulations (Hotokezaka et al. 2013; Bauswein et al. 2013; Sekiguchi et al. 2016; Ciolfi et al. 2017; Radice et al. 2018b; Siegel et al. 2019) and observations of GW 170817/AT2017gfo (Abbott et al. 2017d). FXTs related to the free and trapped zones appear to cover a different parameter-space region (see Fig. 8). The magnetar parameters derived for FXTs 15 and 20 carry large uncertainties, and as a result, they are consistent with both the free zone and trapped zone models (even considering the uncertainties related to redshifts and efficiency, which we discuss below). Only the well-constrained magnetar parameters of FXT 14 differ significantly from the free-zone FXTs (see Fig. 8).

In general, this sample of FXTs populates the lower tail of the magnetic field distribution of the magnetars suggested to drive SGRBs (i.e., ≳1015 G; see Fig. 8). Just a few outliers have ≲1015 G (e.g., FXTs 19, and 22). Papers I and II suggest a non-association with on-axis LGRBs based on the lack of gamma-ray detections and the low luminosity of FXTs, although an association with low-luminosity LGRBs remains possible. The whole sample is above the breakup limit (i.e., Pi ≲ 1 ms; Lattimer & Prakash 2004) except for FXTs 15 and 20, although their parameters have large uncertainties or are only lower limits (such as FXT 20).

Because the magnetar model depends on the assumed efficiency, it is important to discuss the implications of the assumed value for η during the comparison with other transients. As we explained above (see Sect. 4), we assumed an efficiency of η ∼ 10−3 based on the previous analysis of FXTs CDF-S XT1 and XT2 (Sun et al. 2019); however, this parameter might take other values such as in previous works. For instance, Rowlinson et al. (2013) assumed an efficiency in the conversion of rotational energy into EM radiation of 100%. In Appendix B, we discuss the changes in the magnetar parameters adopting a range of redshifts (especially relevant for the FXTs without measured redshifts) and efficiencies. Figure B.1 depicts the best-fit magnetar parameters under different redshifts and efficiencies. These higher efficiency values would increase our magnetar parameters by a factor between ≈5 and 7, that is, rescaling our results in the right-top direction of Fig. 8. This displacement would mean that the FXT magnetar parameters overlap with those of SGRBs, implying that the magnetar nature of this sample of FXTs has similar parameters to those derived from the X-ray afterglows of SGRBs (e.g., Rowlinson et al. 2010, 2013; Gompertz et al. 2014; Lü et al. 2015).

Finally, we estimate the volumetric rate (rescaling the values obtained by Quirola-Vásquez et al. 2023 by a factor of 9/17) of the nine FXTs to ≈800−470 Gpc−3 yr−1. This value is consistent with those observed for GW 170817 ( 1 . 5 1.2 + 3.2 × 10 3 $ 1.5_{-1.2}^{+3.2} \times 10^3 $ Gpc−3 yr−1; Abbott et al. 2017c) and GRB 170817A ( 190 160 + 440 $ 190_{-160}^{+440} $ Gpc−3 yr−1; Zhang et al. 2018). However, EM observations of GW 170817/GRB 170817A offer no conclusive evidence to rule out a long-lived or stable NS (e.g., Abbott et al. 2017c; Piro et al. 2019; Troja et al. 2020). The presence of a blue kilonova component related to a large mass of lanthanide-free ejecta and kinetic energy of Ek ≈ 1051 erg, together with a successful relativistic jet, strongly disfavors a prompt collapse to a BH and suggests a first stage HMNS formation that collapsed quickly to a BH after the merger (Granot et al. 2017; Margalit & Metzger 2017; Shibata et al. 2017; Metzger et al. 2018b; Rezzolla et al. 2018; Murguia-Berthier et al. 2021). The uncertain nature of the remnant does not permit a direct comparison with FXTs.

Considering the possible BNS remnant channels of SGRBs, only the SMNS and stable-NS scenarios might be related to the FXT-magnetar sample7. Under both scenarios, the percentage of mergers leading to the formation of SMNS and a stable NS are ∼18–65% and ≲3%8 (Piro et al. 2017; Metzger 2019; Margalit & Metzger 2019; Patricelli & Bernardini 2020), respectively, and the associated rates should be ∼530 and 30 Gpc−3 yr−1, respectively. Comparing both values and the FXT-magnetar sample, an association between the nine FXTs and the SMNS channel seems possible. Similarly, in the case of SGRBs (the rate is ∼(2−0.1) × 103 Gpc−3 yr−1 at z ≲ 1.5, after correcting by a nominal beaming factor of ≈110 and considering a Gaussian merger delay model; Sun et al. 2015; Wanderman & Piran 2015), the rate of the SMNS channel is ∼1300−36 Gpc−3 yr−1, which is also consistent with this sub-sample of FXTs.

Therefore, a possible association with off-axis SGRBs (due to the lack of a gamma-ray counterpart) and the BNS merger progenitor channel seems plausible, reinforcing previous conclusions drawn from the host-galaxy properties, lack of gamma-ray detections and low luminosities (Papers I and II). The consistent range of FXT magnetar parameters also suggests a unified origin, with the distinguishing parameter being the viewing angle.

6. Merger-nova emission in FXTs 14 and 22?

Optical and NIR observations have been taken for FXTs 14 and 22 from tens of minutes to some days after the X-ray trigger, enabling constraints on contemporaneous counterparts. Unfortunately, however, no simultaneous counterparts have been identified (Bauer et al. 2017; Jonker et al. 2021; Andreoni et al. 2021c). Nevertheless, additional parameter space can be constrained from optical and NIR upper limits, provided the magnetar interpretation and merger-nova emission are plausible. Such contemporaneous observations are not available for the other FXTs, which we exclude from this line of discussion. In Sects. 6.1 and 6.2, we compare the most recent optical and NIR upper limits to the merger-nova theoretical models for FXTs 14 and 22, respectively.

6.1. FXT 14

FXT 14 was serendipitously discovered in the Chandra Deep Field South (CDF-S) on 2014 October 1 (2014-10-01 07:04:37 UT; Luo et al. 2014; Bauer et al. 2017). The CDF-S field has been extensively observed across the EM spectrum, and the position of FXT 14 was serendipitously observed ≈80 min after the X-ray trigger by the Visible Multi-Object Spectrograph (VIMOS) instrument mounted in the Very Large Telescope (VLT). No simultaneous optical transient was detected to an estimated magnitude limit of mR ≳ 25.7 mag at 2σ confidence level. Subsequently, the field of the transient was imaged at ∼18, 27, and 111 days after the X-ray trigger by the VLT-FOcal Reducer/low dispersion Spectrograph 2 (FORS2), Gemini South-Gemini Multi-Object Spectrographs (GMOS), and the Hubble Space Telescope (HST)-WFC3 instrument, respectively. Again, no optical or NIR counterpart was detected, yielding upper limits at the position of the transient of mR ≈ 27.0, mr ≈ 26.0, and mF110W ≈ 28.4 AB mag at 2σ confidence level for days ∼18, 27, and 111, respectively.

Using the best-fit parameters obtained for FXT 14 (see Table 2), we compute the numerical solution of the merger-nova emission in different energy bands (from Eq. (9)). Figure 9, left panel, depicts the merger-nova emission model of FXT 14 considering different filters from u- to Ks-bands (solid color lines) and a comparison with the upper-limits of FXT 14 in the R, r and F110W filters (F110W upper limit could be compared with J-band) adopting the redshift z = 2.23. Overall, the upper limits remain consistent with the expected merger-nova theoretical model. Specifically, around the VLT-VIMOS observation (≈80 min or ≈0.055 days after the trigger), we see that the upper limit still accommodates the peak emission in r-band from the numerical model for redshifts z ≳ 1.8. Notably, this provides a more secure lower bound on the redshift of FXT 14 than the existing HST-based photometric redshift based on HST and deep ground-based imaging.

thumbnail Fig. 9.

Numerical solution of the merger-nova emission in different energy bands. Left panel: optical and NIR upper limits of FXT 14 and the expected merger-nova emission (from u- to K-bands) assuming their best fitting parameters (see Table 2) and z = 2.23 (solid lines). As a middle and right panels: optical and NIR upper limits of FXT 22, and the merger-nova emission considering the derived magnetar parameters (see Table 2), two redshifts of z = 1.0 (solid lines) and 0.1 (dashed lines), and assuming lanthanide-free (middle panel; κ = 0.1 cm2 g−1 and Mej = 10−4M) and rich (right panel; κ = 30 cm2 g−1 and Mej = 10−3M) scenarios. The upper limits of FXT 14 were taken from Bauer et al. (2017), while for FXT 22 from Andreoni et al. (2021a,b,c), Xin et al. (2021), Rossi et al. (2021), and Eappachen et al. (2023). The time is measured from the X-ray trigger detection of the FXTs.

6.2. FXT 22

FXT 22 was serendipitously detected by a Chandra calibration observation of Abell 1795 on 2021 April 23 (2021-04-23 22:15:36.63 UT; Lin et al. 2021). No X-ray emission was detected in previous Chandra observations at this location. Three extended optical objects have been identified as host galaxy candidates (denoted cX, cNE, and cW; Eappachen et al. 2023). The transient location was imaged several times from ≈1 day to 45 days after the outburst by different telescopes and instruments.

The Zwicky Transient Facility (ZTF) observed the field of the transient three days before the outburst and one day after the X-ray trigger (Andreoni et al. 2021a). The forced photometry data, using the ZTF point-spread-function (PSF), was subtracted at the location of FXT 22. No source was detected in individual images in the g-, r-, and i-bands, just upper limits of g > 20.5, r > 20.9, and i > 19.6 AB mag (5σ confidence value). Subsequently, the region was imaged by FORS2 on the 8 m VLT telescope on 6 May 2021 (≈13 days after the X-ray trigger) in the R-band filter. No optical counterpart was found within the X-ray uncertainty region to a depth of R = 24.7 mag at 3σ confidence level (Eappachen et al. 2023). The transient location was subsequently imaged with the Wafer-Scale Imager for Prime (WaSP) instrument mounted on the 200-inch Hale Telescope at Palomar Observatory (Andreoni et al. 2021b) on 6 May 2021 at 07:10:00 UT in the i- and r-bands. No optical counterpart was identified, with upper limits of i > 24.8 and r > 25.2 mag. On 7 May 2021, the field was observed by the NAOC 2.16 m telescope at Xinglong Observatory with the I-band filter, yielding a non-detection with a limiting magnitude of I ≳ 20.5 mag in stacked image (Xin et al. 2021). On 8 May 2021 at 07:28:00 UT (i.e., ≈15 days after the transient detection), the field was visited by the 8.4 m Large Binocular Telescope (LBT) using the r-sloan and z-sloan bands with the Large Binocular Cameras (LBC; Rossi et al. 2021). No detection at the position of the transients was identified, with limit magnitudes of r > 26.1 and z > 25.1 mag. Finally, on 10 June 2021 (i.e., ≈45 days after the onset of the outburst), the quintuple-beam imager HiPERCAM instrument mounted on the 10.4 m Gran Telescopio Canarias (GTC) observed the field of the source with simultaneous us, gs, rs, is, and zs-band filters. No detection at the position of the X-ray transients was identified, with limiting magnitudes of us > 26.2, gs > 27.0, rs > 26.1, is > 24.4, and zs > 24.7 mag at 3σ confidence level (Eappachen et al. 2023).

As noted in Sect. 2.2, the difference between the free- and trapped zones is geometrical (see Fig. 1), and events such as FXT 22 (which is related to the free zone emission based on its light curve plateau) could ultimately be associated with a line of sight that passed through the trapped zone where the ejecta mass was relatively low. Under this condition, FXT 22 would follow the relation tτ ≲ TX (i.e., the time at which point the ejecta material becomes optically thin occurs before the plateau appears). In this way, we explore in Fig. 9 whether the merger-nova models, assuming either the lanthanide-free (κ = 0.1 cm2 g−1 and Mej = 10−4M) or -rich (κ = 30 cm2 g−1 and Mej = 10−3M) scenarios, are consistent with the optical upper limits obtained for FXT 22. We find that the upper limits remain consistent with the expected merger-nova theoretical model for redshifts z ≳ 0.1 (see Fig. 9, middle and right panels, dashed lines), and hence remain completely consistent with potential hosts cW at zcW = 1, 04 and cNE at zcNE ≈ 1.5. This limit also implies that the strongest host candidate, the faint source cX, visible in GTC-HiPERCAM g-filter images only (Eappachen et al. 2023), should lie at zcX ≳ 0.1.

7. Rotational energy losses via GWs and detection probability

A fraction of the rotational energy of the magnetar is likely released by GW radiation (Fan et al. 2013; Lasky et al. 2014; Lan et al. 2020), produced via mass quadrupole deformation (related to the ellipticity, ϵ) considering rigid-body rotation (Zhang & Mészáros 2001; Fan et al. 2013; Ho 2016; Lasky & Glampedakis 2016; Lü et al. 2017). There are several possibilities for deforming a newly born magnetar (e.g., a bar-mode instability or inertial quadrupole r-mode; Andersson & Kokkotas 2001; Cutler 2002; Haskell et al. 2008; Corsi & Mészáros 2009; Lasky & Glampedakis 2016). Nevertheless, the most important contribution to ϵ may come from magnetic induction. For a differentially rotating magnetar, its magnetic fields become twisted and stored in a toroidal shape, causing it becomes non-spherical (Cutler 2002; Haskell et al. 2008). The condition to produce significant GW radiation is for the magnetic axis to be strongly misaligned with the rotational axis (Shapiro & Teukolsky 1983; Zhang 2018).

We briefly explore a potential GW emission contribution in the context of the magnetar model. As we explained above, the millisecond magnetar also loses rotational energy via GW emission, where the characteristic spin-down timescale (Eq. (A.5)) of the magnetar can be numerically written as

τ GW = 5 c 5 P i 4 128 G I ϵ 2 ( 2 π ) 4 9.1 × 10 3 ( I 10 45 ) 1 ( ϵ 10 3 ) 2 ( P i 10 3 ) 4 s . $$ \begin{aligned} \tau _{\rm GW}=\frac{5c^5P_{\rm i}^4}{128GI\epsilon ^2(2\pi )^4} \simeq 9.1 \times 10^{3}\left(\frac{I}{10^{45}}\right)^{-1}\left(\frac{\epsilon }{10^{-3}}\right)^{-2} \left(\frac{P_{\rm i}}{10^{-3}}\right)^4\,\mathrm{s}. \end{aligned} $$(10)

Given that the decay phase of the FXTs remains inconsistent with a GW dominant spin-down luminosity trend (see Eq. (A.7)), we expect that τGW > TX (where TX ≳ 1 ks), and combined with Eq. (10), an upper limit on the ellipticity can be expressed as

ϵ < 2.5 × 10 3 ( I 10 45 ) 1 / 2 ( P i 10 3 ) 2 . $$ \begin{aligned} \epsilon < 2.5 \times 10^{-3}\left(\frac{I}{10^{45}}\right)^{-1/2}\left(\frac{P_{\rm i}}{10^{-3}}\right)^2. \end{aligned} $$(11)

To compute the maximum value of ϵ, we should take into account different EoSs (see Table 3). Based on these, the maximum values of ϵ, assuming Pi ≈ 1 ms (see Table 2), lie in the range ϵ ≲ 1.2 × 10−3 − 2.1 × 10−4. This range is lower than the expected values when GW radiation dominates the rotational energy losses (i.e., ϵ ∼ 10−2; Fan et al. 2013; Lasky et al. 2014; Ho 2016; Lan et al. 2020), reinforcing the assumed value of ϵ ∼ 10−5 during the fitting process.

If this new population of magnetars does lose energy via GWs, it should be possible to constrain their GW emission considering a set of EoSs (see Table 3), as well as predict their detectability with the current Advanced LIGO detector (aLIGO) and more sensitive future detectors such as the Einstein Telescope (ET) and Cosmic Explorer (CE). Specifically, if the rotation energy is released via GWs with a frequency f, the GW strain (h(t)) for a rotating NS at luminosity distance dL can be expressed as (Lü et al. 2017, 2019, and references therein)

h ( t ) = 4 G I ϵ d L c 4 Ω ( t ) 2 , $$ \begin{aligned} h(t)=\frac{4GI\epsilon }{d_{\rm L} c^4}\Omega (t)^2, \end{aligned} $$(12)

where the characteristic amplitude of GWs from a rotating NS can be estimated as (Corsi & Mészáros 2009; Hild et al. 2011; Lasky & Glampedakis 2016; Lü et al. 2017)

h c = f h ( t ) d t d f = f d L 5 G I 2 c 3 f 8.22 × 10 24 ( I 10 45 g cm 2 ) 1 / 2 ( f 1000 Hz ) 1 / 2 ( d L 100 Mpc ) 1 , $$ \begin{aligned} h_{\rm c}&= fh(t)\sqrt{\frac{\mathrm{d}t}{\mathrm{d}f}}=\frac{f}{d_{\rm L}} \sqrt{\frac{5GI}{2c^3f}}\nonumber \\&{\simeq }8.22\times 10^{-24}\left(\frac{I}{10^{45}\,\mathrm{g\,cm^{2}}}\right)^{1/2}\left(\frac{f}{1000\,\mathrm{Hz}}\right)^{1/2}\left(\frac{d_{\rm L}}{100\,\mathrm{Mpc}}\right)^{-1}, \end{aligned} $$(13)

where f is the frequency in Hz and dL is in Mpc, and the gravitational wave frequency evolution, df/dt, is derived directly from Eq. (A.2). It is clear that hc does not depend on ϵ, just dL and I. Thus, each EoS will provide a different characteristic amplitude. It is important to emphasize that future GW detectors should also detect the GWs from the BNS merger, which precedes the FXT magnetars.

Figure 10 compares the expected hc considering eight EoSs from Table 3 (shaded regions) for the FXTs analyzed in this work with the sensitivities of aLIGO (cyan line), ET (gray line), and CE (black line). Overall, the maximum strain values, hc, are for the EoS NL3ωρ, which has the highest inertia (see Table 3). All of the FXTs remain firmly beyond the current capabilities of aLIGO (i.e., z ≲ 0.03 or ≲130 Mpc). ET will be sensitive to events such as FXTs 7, 10, 12, 15, and 20 if they lie at redshifts z ≲ 0.5, while CE should push out to z ≲ 0.75 (allowing the inclusion of events like FXT 16). No GW observatory, existing or currently planned, will be able to detect events similar to FXTs 14, 19, and 22 (see Fig. 10).

thumbnail Fig. 10.

GW strain evolution with frequency for the FXTs analyzed in this work. Each shaded region represents the GW strain considering eight EoS (see Table 3) at different distances. The cyan, gray, and black lines are the sensitivity limits for aLIGO, ET, and CE GW telescopes, respectively. For the FXTs without redshifts (FXTs 7, 10, 12, 15, and 20), we assumed z = 0.5 and 1.0 (see Table 2).

8. Conclusions

Papers I and II reported the detection of 22 extragalactic fast X-ray transients (FXTs) hidden in two decades of Chandra data (analyzing ∼259 Ms of data), among which 17 are associated with distant galaxies (≳100 Mpc). Different mechanisms and progenitors have been proposed to explain their properties; nevertheless, their nature remains uncertain. We interpreted a subset of nine distant FXTs within the framework of a two-zone (“free” and “trapped”) BNS merger magnetar model, parametrized by a magnetic field (BP), initial rotational period (Pi), ejecta mass (Mej), opacity (κ), and the fraction of the magnetar energy that is transferred to the heating of the ejecta (ξ).

The model can explain well the observed X-ray light curves of the sources, wherein we associate six FXTs to a free-zone scenario and three FXTs to a trapped-zone scenario. The best-fitted results yield relatively similar BP and Pi values in the ranges 1016 − 1015 G and 7.9 − 1.0 ms, respectively. These values are comparable to those derived from the magnetar interpretation of some X-ray afterglow of SGRBs (for instance Rowlinson et al. 2010, 2013; Lü et al. 2015), reinforcing an association with off-axis SGRBs as suggested previously (Bauer et al. 2017; Xue et al. 2019; Lin et al. 2022; Quirola-Vásquez et al. 2022, 2023).

For two FXTs, XRT 141001 (FXT 14) and XRT 210423 (FXT 22), which have prompt deep optical and near-infrared upper limits available, we explored the constraints on early merger-nova detections. We found that the merger-nova numerical models agree with the optical and near-infrared upper limits, and we identified that the non-detections place lower limits on the distance of z ≳ 1.8 and ≳0.1 for XRT 141001/FXT 14 and XRT 210423/FXT 22, respectively.

Finally, we analyzed the contribution of the gravitational wave (GW) radiation to the magnetar spin-down (related to the ellipticity factor, ϵ, see Eq. (A.2)), and the detectability of GWs associated with our sample of FXTs by current and future GW observatories. The maximum value of ϵ, considering different equations of state models (see Table 3), is ϵ ≲ 10−3 − 10−4, which is below the expected values for GW losses to dominate the system (i.e., ϵ ≳ 10−2; Fan et al. 2013; Lasky et al. 2014; Ho 2016; Lan et al. 2020). We find that it is only possible to detect BNS magnetar remnants similar to the FXTs with the current Advanced LIGO detectors to z ≲ 0.03, while future GW detectors should extend our sensitivity to z ≲ 0.5 − 0.75.


1

An “internal plateau” is a plateau followed by a quick decay (power-law index of ≲ − 3), which is not explained by the evolution of the synchrotron emission from a decelerating forward shock, but could be interpreted within the central engine framework (e.g., Zhang & Mészáros 2001; Tang et al. 2019).

2

We are assuming a pure dipole magnetic field. Nevertheless, it is known that neither purely poloidal nor toroidal (these fields are responsible for causing surface currents in magnetars; Soldateschi et al. 2021) magnetic topologies are stable.

3

The parameters in the co-moving frames are denoted as Q′. Parameters without a prime (′) are in the observer frame.

4

We discarded events that do not follow the expected light curve decays: both XRT 000519 and 110103 have a substructure in their main peaks (Jonker et al. 2013; Glennie et al. 2015; Quirola-Vásquez et al. 2022) which does not keep an explanation under the magnetar model; in contrast, although the others have a power-law decay, their slopes are not consistent with the expected by the magnetar model (see Sects. 2.1 and 2.2).

5

In Appendix B, we discuss the role of the efficiency in deriving the magnetar parameters.

6

In Appendix B, we discuss the role of the nominal redshifts in deriving the magnetar parameters.

7

The rate of BNS is 320 240 + 490 $ 320^{+490}_{-240} $ Gpc−3 yr−1 (Abbott et al. 2021).

8

Assuming the merging extragalactic NS-NS binary population is identical to the known Galactic NS–NS binaries.

9

The parameters in the co-moving frames are denoted as Q′. Parameters without a prime (′) are in the observer frame.

Acknowledgments

The scientific results reported in this article are based on observations made by the Chandra X-ray Observatory. This research has made use of software provided by the Chandra X-ray Center (CXC). We acknowledge support from: CONICYT through Programa de Capital Humano Avanzado, folio #21180886 (JQ-V) Basal AFB-170002 (JQ-V, FEB), and FONDECYT Regular 1141218 (JQ-V, FEB); the Ministry of Economy, Development, and Tourism’s Millennium Science Initiative through grant IC120009, awarded to The Millennium Institute of Astrophysics, MAS (JQ-V, FEB); this project was (partially) funded by NWO under grant number 184.034.002 (P.G.J.); G.Y. acknowledges funding from the Netherlands Research School for Astronomy (NOVA). NSF grant AST-2106990 and Chandra X-ray Center grant GO0-21080X (W.N.B.); the National Natural Science Foundation of China grant 11991053 (B.L.); the NSFC grant (12025303), the Strategic Priority Research Program of the Chinese Academy of Sciences (grant NO. XDB0550300; Y.Q.X.); H.S. is supported by the National Natural Science Foundation of China grant 12103065.

References

  1. Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017a, Phys. Rev. Lett., 119, 161101 [Google Scholar]
  2. Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017b, ApJ, 848, L12 [Google Scholar]
  3. Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017c, ApJ, 848, L13 [CrossRef] [Google Scholar]
  4. Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017d, ApJ, 850, L39 [NASA ADS] [CrossRef] [Google Scholar]
  5. Abbott, R., Abbott, T. D., Abraham, S., et al. 2021, Phys. Rev. X, 11, 021053 [Google Scholar]
  6. Ai, S., & Zhang, B. 2021, ApJ, 915, L11 [NASA ADS] [CrossRef] [Google Scholar]
  7. Ai, S., Gao, H., Dai, Z.-G., et al. 2018, ApJ, 860, 57 [NASA ADS] [CrossRef] [Google Scholar]
  8. Alexander, K. D., Berger, E., Fong, W., et al. 2017, ApJ, 848, L21 [NASA ADS] [CrossRef] [Google Scholar]
  9. Alp, D., & Larsson, J. 2020, ApJ, 896, 39 [NASA ADS] [CrossRef] [Google Scholar]
  10. Andersson, N., & Kokkotas, K. D. 2001, Int. J. Mod. Phys. D, 10, 381 [NASA ADS] [CrossRef] [Google Scholar]
  11. Andreoni, I., Kasliwal, M., & Graham, M. 2021a, ATel, 14600, 1 [Google Scholar]
  12. Andreoni, I., De, K., Kasliwal, M., & Tzanidakis, A. 2021b, ATel, 14640, 1 [Google Scholar]
  13. Andreoni, I., Perley, D., Kasliwal, M., et al. 2021c, ATel, 14641, 1 [NASA ADS] [Google Scholar]
  14. Antoniadis, J., Freire, P. C. C., Wex, N., et al. 2013, Science, 340, 448 [Google Scholar]
  15. Arcavi, I., Hosseinzadeh, G., Howell, D. A., et al. 2017, Nature, 551, 64 [NASA ADS] [CrossRef] [Google Scholar]
  16. Banerjee, S., Tanaka, M., Kawaguchi, K., Kato, D., & Gaigalas, G. 2020, ApJ, 901, 29 [NASA ADS] [CrossRef] [Google Scholar]
  17. Barnes, J., & Kasen, D. 2013, ApJ, 775, 18 [NASA ADS] [CrossRef] [Google Scholar]
  18. Bauer, F. E., Treister, E., Schawinski, K., et al. 2017, MNRAS, 467, 4841 [NASA ADS] [CrossRef] [Google Scholar]
  19. Bauswein, A., Baumgarte, T. W., & Janka, H. T. 2013, Phys. Rev. Lett., 111, 131101 [NASA ADS] [CrossRef] [Google Scholar]
  20. Berger, E. 2014, ARA&A, 52, 43 [CrossRef] [Google Scholar]
  21. Berger, E., Fong, W., & Chornock, R. 2013, ApJ, 774, L23 [NASA ADS] [CrossRef] [Google Scholar]
  22. Bucciantini, N., Metzger, B. D., Thompson, T. A., & Quataert, E. 2012, MNRAS, 419, 1537 [NASA ADS] [CrossRef] [Google Scholar]
  23. Chornock, R., Berger, E., Kasen, D., et al. 2017, ApJ, 848, L19 [NASA ADS] [CrossRef] [Google Scholar]
  24. Ciolfi, R., Kastaun, W., Giacomazzo, B., et al. 2017, Phys. Rev. D, 95, 063016 [Google Scholar]
  25. Corsi, A., & Mészáros, P. 2009, ApJ, 702, 1171 [NASA ADS] [CrossRef] [Google Scholar]
  26. Coulter, D. A., Foley, R. J., Kilpatrick, C. D., et al. 2017, Science, 358, 1556 [NASA ADS] [CrossRef] [Google Scholar]
  27. Cowperthwaite, P. S., Berger, E., Villar, V. A., et al. 2017, ApJ, 848, L17 [CrossRef] [Google Scholar]
  28. Cutler, C. 2002, Phys. Rev. D, 66, 084025 [NASA ADS] [CrossRef] [Google Scholar]
  29. Dai, Z. G., & Lu, T. 1998, A&A, 333, L87 [NASA ADS] [Google Scholar]
  30. Dai, Z. G., Wang, X. Y., Wu, X. F., & Zhang, B. 2006, Science, 311, 1127 [NASA ADS] [CrossRef] [Google Scholar]
  31. Dai, L., McKinney, J. C., Roth, N., Ramirez-Ruiz, E., & Miller, M. C. 2018, ApJ, 859, L20 [Google Scholar]
  32. Dainotti, M. G., Nagataki, S., Maeda, K., Postnikov, S., & Pian, E. 2017, A&A, 600, A98 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Dall’Osso, S., Shore, S. N., & Stella, L. 2009, MNRAS, 398, 1869 [CrossRef] [Google Scholar]
  34. Dean, C., Fernández, R., & Metzger, B. D. 2021, ApJ, 921, 161 [NASA ADS] [CrossRef] [Google Scholar]
  35. D’Elia, V., Campana, S., D’Aì, A., et al. 2018, A&A, 619, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Demorest, P. B., Pennucci, T., Ransom, S. M., Roberts, M. S. E., & Hessels, J. W. T. 2010, Nature, 467, 1081 [Google Scholar]
  37. Drout, M. R., Piro, A. L., Shappee, B. J., et al. 2017, Science, 358, 1570 [NASA ADS] [CrossRef] [Google Scholar]
  38. Duncan, R. C., & Thompson, C. 1992, ApJ, 392, L9 [Google Scholar]
  39. Eappachen, D., Jonker, P. G., Fraser, M., et al. 2022, MNRAS, 514, 302 [NASA ADS] [CrossRef] [Google Scholar]
  40. Eappachen, D., Jonker, P. G., Levan, A. J., et al. 2023, ApJ, 948, 91 [NASA ADS] [CrossRef] [Google Scholar]
  41. Evans, P. A., Cenko, S. B., Kennea, J. A., et al. 2017, Science, 358, 1565 [NASA ADS] [CrossRef] [Google Scholar]
  42. Fan, Y.-Z., & Xu, D. 2006, MNRAS, 372, L19 [NASA ADS] [CrossRef] [Google Scholar]
  43. Fan, Y.-Z., Wu, X.-F., & Wei, D.-M. 2013, Phys. Rev. D, 88, 067304 [NASA ADS] [CrossRef] [Google Scholar]
  44. Fernández, R., & Metzger, B. D. 2016, Annu. Rev. Nucl. Part. Sci., 66, 23 [CrossRef] [Google Scholar]
  45. Feroz, F., & Hobson, M. P. 2008, MNRAS, 384, 449 [NASA ADS] [CrossRef] [Google Scholar]
  46. Feroz, F., Hobson, M. P., & Bridges, M. 2009, MNRAS, 398, 1601 [NASA ADS] [CrossRef] [Google Scholar]
  47. Fong, W., Berger, E., Margutti, R., & Zauderer, B. A. 2015, ApJ, 815, 102 [NASA ADS] [CrossRef] [Google Scholar]
  48. Fonseca, E., Cromartie, H. T., Pennucci, T. T., et al. 2021, ApJ, 915, L12 [NASA ADS] [CrossRef] [Google Scholar]
  49. Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [Google Scholar]
  50. Foucart, F., Desai, D., Brege, W., et al. 2017, Class. Quant. Grav., 34, 044002 [NASA ADS] [CrossRef] [Google Scholar]
  51. Gao, W.-H., & Fan, Y.-Z. 2006, Chin. J. Astron. Astrophys., 6, 513 [NASA ADS] [CrossRef] [Google Scholar]
  52. Gao, H., Zhang, B., & Lü, H.-J. 2016, Phys. Rev. D, 93, 044065 [NASA ADS] [CrossRef] [Google Scholar]
  53. Gao, H., Zhang, B., Lü, H.-J., & Li, Y. 2017, ApJ, 837, 50 [NASA ADS] [CrossRef] [Google Scholar]
  54. Glennie, A., Jonker, P. G., Fender, R. P., Nagayama, T., & Pretorius, M. L. 2015, MNRAS, 450, 3765 [NASA ADS] [CrossRef] [Google Scholar]
  55. Goldstein, A., Veres, P., Burns, E., et al. 2017, ApJ, 848, L14 [CrossRef] [Google Scholar]
  56. Gompertz, B. P., O’Brien, P. T., & Wynn, G. A. 2014, MNRAS, 438, 240 [NASA ADS] [CrossRef] [Google Scholar]
  57. Goriely, S., Bauswein, A., & Janka, H.-T. 2011, ApJ, 738, L32 [Google Scholar]
  58. Granot, J., Guetta, D., & Gill, R. 2017, ApJ, 850, L24 [NASA ADS] [CrossRef] [Google Scholar]
  59. Hajela, A., Margutti, R., Bright, J. S., et al. 2022, ApJ, 927, L17 [NASA ADS] [CrossRef] [Google Scholar]
  60. Hallinan, G., Corsi, A., Mooley, K. P., et al. 2017, Science, 358, 1579 [NASA ADS] [CrossRef] [Google Scholar]
  61. Haskell, B., Samuelsson, L., Glampedakis, K., & Andersson, N. 2008, MNRAS, 385, 531 [NASA ADS] [CrossRef] [Google Scholar]
  62. Hild, S., Abernathy, M., Acernese, F., et al. 2011, Class. Quant. Grav., 28, 094013 [Google Scholar]
  63. Ho, W. C. G. 2016, MNRAS, 463, 489 [NASA ADS] [CrossRef] [Google Scholar]
  64. Hotokezaka, K., Kiuchi, K., Kyutoku, K., et al. 2013, Phys. Rev. D, 87, 024001 [NASA ADS] [CrossRef] [Google Scholar]
  65. Ide, S., Hayashida, K., Noda, H., et al. 2020, PASJ, 72, 40 [NASA ADS] [CrossRef] [Google Scholar]
  66. Irwin, J. A., Maksym, W. P., Sivakoff, G. R., et al. 2016, Nature, 538, 356 [NASA ADS] [CrossRef] [Google Scholar]
  67. Jin, Z.-P., Hotokezaka, K., Li, X., et al. 2016, Nat. Commun., 7, 12898 [NASA ADS] [CrossRef] [Google Scholar]
  68. Jonker, P. G., Glennie, A., Heida, M., et al. 2013, ApJ, 779, 14 [NASA ADS] [CrossRef] [Google Scholar]
  69. Jonker, P., Levan, A., Torres, M., Eappachen, D., & Quirola, J. 2021, Transient Name Server AstroNote, 160, 1 [NASA ADS] [Google Scholar]
  70. Kasen, D., & Bildsten, L. 2010, ApJ, 717, 245 [NASA ADS] [CrossRef] [Google Scholar]
  71. Kasen, D., Metzger, B., Barnes, J., Quataert, E., & Ramirez-Ruiz, E. 2017, Nature, 551, 80 [Google Scholar]
  72. Kilpatrick, C. D., Foley, R. J., Kasen, D., et al. 2017, Science, 358, 1583 [NASA ADS] [CrossRef] [Google Scholar]
  73. Kiuchi, K., Kyutoku, K., Sekiguchi, Y., Shibata, M., & Wada, T. 2014, Phys. Rev. D, 90, 041502 [Google Scholar]
  74. Kiziltan, B., Kottas, A., Yoreo, M. D., & Thorsett, S. E. 2013, ApJ, 778, 66 [NASA ADS] [CrossRef] [Google Scholar]
  75. Korobkin, O., Rosswog, S., Arcones, A., & Winteler, C. 2012, MNRAS, 426, 1940 [NASA ADS] [CrossRef] [Google Scholar]
  76. Kotera, K., Phinney, E. S., & Olinto, A. V. 2013, MNRAS, 432, 3228 [NASA ADS] [CrossRef] [Google Scholar]
  77. Kyutoku, K., Ioka, K., & Shibata, M. 2013, Phys. Rev. D, 88, 041503 [NASA ADS] [CrossRef] [Google Scholar]
  78. Kyutoku, K., Ioka, K., Okawa, H., Shibata, M., & Taniguchi, K. 2015, Phys. Rev. D, 92, 044028 [NASA ADS] [CrossRef] [Google Scholar]
  79. Lü, H.-J., Zhang, H.-M., Zhong, S.-Q., et al. 2017, ApJ, 835, 181 [CrossRef] [Google Scholar]
  80. Lamb, G. P., & Kobayashi, S. 2017, MNRAS, 472, 4953 [NASA ADS] [CrossRef] [Google Scholar]
  81. Lan, L., Lü, H.-J., Rice, J., & Liang, E.-W. 2020, ApJ, 890, 99 [NASA ADS] [CrossRef] [Google Scholar]
  82. Lasky, P. D., & Glampedakis, K. 2016, MNRAS, 458, 1660 [NASA ADS] [CrossRef] [Google Scholar]
  83. Lasky, P. D., Haskell, B., Ravi, V., et al. 2014, Phys. Rev. D, 89, 047302 [NASA ADS] [CrossRef] [Google Scholar]
  84. Lasky, P. D., Leris, C., Rowlinson, A., & Glampedakis, K. 2017, ApJ, 843, L1 [NASA ADS] [CrossRef] [Google Scholar]
  85. Lattimer, J. M., & Prakash, M. 2004, Science, 304, 536 [NASA ADS] [CrossRef] [Google Scholar]
  86. Lazzati, D., Deich, A., Morsony, B. J., & Workman, J. C. 2017, MNRAS, 471, 1652 [NASA ADS] [CrossRef] [Google Scholar]
  87. Levan, A. J., Malesani, D. B., Gompertz, B. P., et al. 2023, Nature, 7, 976 [NASA ADS] [Google Scholar]
  88. Levan, A., Gompertz, B. P., Salafia, O. S., et al. 2024, Nature, 626, 737 [NASA ADS] [CrossRef] [Google Scholar]
  89. Li, L.-X., & Paczyński, B. 1998, ApJ, 507, L59 [NASA ADS] [CrossRef] [Google Scholar]
  90. Li, Y., Zhang, B., & Lü, H.-J. 2016a, ApJS, 227, 7 [NASA ADS] [CrossRef] [Google Scholar]
  91. Li, A., Zhang, B., Zhang, N.-B., et al. 2016b, Phys. Rev. D, 94, 083010 [NASA ADS] [CrossRef] [Google Scholar]
  92. Lin, D., Strader, J., Carrasco, E. R., et al. 2018, Nat. Astron., 2, 656 [Google Scholar]
  93. Lin, D., Irwin, J., & Berger, E. 2019, ATel, 13171, 1 [NASA ADS] [Google Scholar]
  94. Lin, D., Strader, J., Romanowsky, A. J., et al. 2020, ApJ, 892, L25 [Google Scholar]
  95. Lin, D., Irwin, J. A., & Berger, E. 2021, ATel, 14599, 1 [NASA ADS] [Google Scholar]
  96. Lin, D., Irwin, J. A., Berger, E., & Nguyen, R. 2022, ApJ, 927, 211 [NASA ADS] [CrossRef] [Google Scholar]
  97. Lipunov, V. M., Gorbovskoy, E., Kornilov, V. G., et al. 2017, ApJ, 850, L1 [NASA ADS] [CrossRef] [Google Scholar]
  98. Lu, W., & Kumar, P. 2018, ApJ, 865, 128 [NASA ADS] [CrossRef] [Google Scholar]
  99. Lü, H.-J., Zhang, B., Lei, W.-H., Li, Y., & Lasky, P. D. 2015, ApJ, 805, 89 [CrossRef] [Google Scholar]
  100. Lü, H.-J., Yuan, Y., Lan, L., et al. 2019, Res. Astron. Astrophys., 21, 047 [Google Scholar]
  101. Luo, B., Brandt, N., & Bauer, F. 2014, ATel, 6541, 1 [NASA ADS] [Google Scholar]
  102. Lyons, N., O’Brien, P. T., Zhang, B., et al. 2010, MNRAS, 402, 705 [NASA ADS] [CrossRef] [Google Scholar]
  103. MacLeod, M., Goldstein, J., Ramirez-Ruiz, E., Guillochon, J., & Samsing, J. 2014, ApJ, 794, 9 [NASA ADS] [CrossRef] [Google Scholar]
  104. Maguire, K., Eracleous, M., Jonker, P. G., MacLeod, M., & Rosswog, S. 2020, Space Sci. Rev., 216, 39 [NASA ADS] [CrossRef] [Google Scholar]
  105. Margalit, B., & Metzger, B. D. 2017, ApJ, 850, L19 [Google Scholar]
  106. Margalit, B., & Metzger, B. D. 2019, ApJ, 880, L15 [NASA ADS] [CrossRef] [Google Scholar]
  107. Margutti, R., Berger, E., Fong, W., et al. 2017, ApJ, 848, L20 [CrossRef] [Google Scholar]
  108. Martinez, J. G., Stovall, K., Freire, P. C. C., et al. 2015, ApJ, 812, 143 [NASA ADS] [CrossRef] [Google Scholar]
  109. Mazzali, P. A., Valenti, S., Della Valle, M., et al. 2008, Science, 321, 1185 [NASA ADS] [CrossRef] [Google Scholar]
  110. Metzger, B. D. 2019, Liv. Rev. Relat., 23, 1 [Google Scholar]
  111. Metzger, B. D., & Piro, A. L. 2014, MNRAS, 439, 3916 [NASA ADS] [CrossRef] [Google Scholar]
  112. Metzger, B. D., Quataert, E., & Thompson, T. A. 2008, MNRAS, 385, 1455 [CrossRef] [Google Scholar]
  113. Metzger, B. D., Martínez-Pinedo, G., Darbha, S., et al. 2010, MNRAS, 406, 2650 [NASA ADS] [CrossRef] [Google Scholar]
  114. Metzger, B. D., Giannios, D., Thompson, T. A., Bucciantini, N., & Quataert, E. 2011, MNRAS, 413, 2031 [Google Scholar]
  115. Metzger, B. D., Beniamini, P., & Giannios, D. 2018a, ApJ, 857, 95 [NASA ADS] [CrossRef] [Google Scholar]
  116. Metzger, B. D., Thompson, T. A., & Quataert, E. 2018b, ApJ, 856, 101 [NASA ADS] [CrossRef] [Google Scholar]
  117. Modjaz, M., Li, W., Butler, N., et al. 2009, ApJ, 702, 226 [NASA ADS] [CrossRef] [Google Scholar]
  118. Murguia-Berthier, A., Ramirez-Ruiz, E., De Colle, F., et al. 2021, ApJ, 908, 152 [NASA ADS] [CrossRef] [Google Scholar]
  119. Nakar, E. 2015, ApJ, 807, 172 [Google Scholar]
  120. Nakar, E., & Sari, R. 2010, ApJ, 725, 904 [NASA ADS] [CrossRef] [Google Scholar]
  121. Narayan, R., Paczynski, B., & Piran, T. 1992, ApJ, 395, L83 [NASA ADS] [CrossRef] [Google Scholar]
  122. Nicholl, M., Guillochon, J., & Berger, E. 2017, ApJ, 850, 55 [NASA ADS] [CrossRef] [Google Scholar]
  123. Nitz, A. H., Kumar, S., Wang, Y.-F., et al. 2023, ApJ, 946, 59 [NASA ADS] [CrossRef] [Google Scholar]
  124. Nousek, J. A., Kouveliotou, C., Grupe, D., et al. 2006, ApJ, 642, 389 [Google Scholar]
  125. Novara, G., Esposito, P., Tiengo, A., et al. 2020, ApJ, 898, 37 [Google Scholar]
  126. Oechslin, R., Janka, H. T., & Marek, A. 2007, A&A, 467, 395 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  127. Pastor-Marazuela, I., Webb, N. A., Wojtowicz, D. T., & van Leeuwen, J. 2020, A&A, 640, A124 [EDP Sciences] [Google Scholar]
  128. Patricelli, B., & Bernardini, M. G. 2020, MNRAS, 499, L96 [NASA ADS] [CrossRef] [Google Scholar]
  129. Piro, A. L., Giacomazzo, B., & Perna, R. 2017, ApJ, 844, L19 [NASA ADS] [CrossRef] [Google Scholar]
  130. Piro, L., Troja, E., Zhang, B., et al. 2019, MNRAS, 483, 1912 [NASA ADS] [CrossRef] [Google Scholar]
  131. Pooley, D., Kumar, P., Wheeler, J. C., & Grossan, B. 2018, ApJ, 859, L23 [NASA ADS] [CrossRef] [Google Scholar]
  132. Price, D. J., & Rosswog, S. 2006, Science, 312, 719 [NASA ADS] [CrossRef] [Google Scholar]
  133. Quirola-Vásquez, J., Bauer, F. E., Jonker, P. G., et al. 2022, A&A, 663, A168 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  134. Quirola-Vásquez, J., Bauer, F. E., Jonker, P. G., et al. 2023, A&A, 675, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  135. Radice, D., Perego, A., Hotokezaka, K., et al. 2018a, ApJ, 869, 130 [Google Scholar]
  136. Radice, D., Perego, A., Hotokezaka, K., et al. 2018b, ApJ, 869, L35 [NASA ADS] [CrossRef] [Google Scholar]
  137. Ramirez-Ruiz, E., Celotti, A., & Rees, M. J. 2002, MNRAS, 337, 1349 [Google Scholar]
  138. Rastinejad, J. C., Gompertz, B. P., Levan, A. J., et al. 2022, Nature, 612, 223 [NASA ADS] [CrossRef] [Google Scholar]
  139. Ravi, V., & Lasky, P. D. 2014, MNRAS, 441, 2433 [NASA ADS] [CrossRef] [Google Scholar]
  140. Rezzolla, L., Most, E. R., & Weih, L. R. 2018, ApJ, 852, L25 [Google Scholar]
  141. Roberts, L. F., Kasen, D., Lee, W. H., & Ramirez-Ruiz, E. 2011, ApJ, 736, L21 [NASA ADS] [CrossRef] [Google Scholar]
  142. Romani, R. W., Kandel, D., Filippenko, A. V., Brink, T. G., & Zheng, W. 2022, ApJ, 934, L17 [NASA ADS] [CrossRef] [Google Scholar]
  143. Rossi, A., Stratta, G., Maiorano, E., et al. 2020, MNRAS, 493, 3379 [CrossRef] [Google Scholar]
  144. Rossi, A., Cappellaro, E., & D’Avanz, P. 2021, ATel, 14635, 1 [Google Scholar]
  145. Rosswog, S. 2007, MNRAS, 376, L48 [NASA ADS] [Google Scholar]
  146. Rowlinson, A., O’Brien, P. T., Tanvir, N. R., et al. 2010, MNRAS, 409, 531 [NASA ADS] [CrossRef] [Google Scholar]
  147. Rowlinson, A., O’Brien, P. T., Metzger, B. D., Tanvir, N. R., & Levan, A. J. 2013, MNRAS, 430, 1061 [NASA ADS] [CrossRef] [Google Scholar]
  148. Ruiz, M., Shapiro, S. L., & Tsokaros, A. 2018, Phys. Rev. D, 97, 021501 [NASA ADS] [CrossRef] [Google Scholar]
  149. Sarin, N., Ashton, G., Lasky, P. D., et al. 2021, ArXiv e-prints [arXiv:2105.10108] [Google Scholar]
  150. Savchenko, V., Ferrigno, C., Kuulkers, E., et al. 2017, ApJ, 848, L15 [NASA ADS] [CrossRef] [Google Scholar]
  151. Saxton, R., Komossa, S., Auchettl, K., & Jonker, P. G. 2021, Space Sci. Rev., 217, 18 [NASA ADS] [CrossRef] [Google Scholar]
  152. Sekiguchi, Y., Kiuchi, K., Kyutoku, K., Shibata, M., & Taniguchi, K. 2016, Phys. Rev. D, 93, 124046 [NASA ADS] [CrossRef] [Google Scholar]
  153. Shapiro, S. L., & Teukolsky, S. A. 1983, Black Holes, White Dwarfs, and Neutron Stars: The Physics of Compact Objects (New York: Wiley) [CrossRef] [Google Scholar]
  154. Shappee, B. J., Simon, J. D., Drout, M. R., et al. 2017, Science, 358, 1574 [NASA ADS] [CrossRef] [Google Scholar]
  155. Sharma, S. 2017, ARA&A, 55, 213 [Google Scholar]
  156. Shibata, M., Fujibayashi, S., Hotokezaka, K., et al. 2017, Phys. Rev. D, 96, 123012 [Google Scholar]
  157. Siegel, D. M., & Ciolfi, R. 2016a, ApJ, 819, 14 [NASA ADS] [CrossRef] [Google Scholar]
  158. Siegel, D. M., & Ciolfi, R. 2016b, ApJ, 819, 15 [NASA ADS] [CrossRef] [Google Scholar]
  159. Siegel, D. M., Barnes, J., & Metzger, B. D. 2019, Nature, 569, 241 [Google Scholar]
  160. Smartt, S. J., Chen, T.-W., Jerkstrand, A., et al. 2017, Nature, 551, 75 [Google Scholar]
  161. Soares-Santos, M., Holz, D. E., Annis, J., et al. 2017, ApJ, 848, L16 [CrossRef] [Google Scholar]
  162. Soderberg, A. M., Berger, E., Page, K. L., et al. 2008, Nature, 454, 246 [NASA ADS] [CrossRef] [Google Scholar]
  163. Soldateschi, J., Bucciantini, N., & Del Zanna, L. 2021, A&A, 654, A162 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  164. Sugita, S., Kawai, N., Nakahira, S., et al. 2018, PASJ, 70, 81 [NASA ADS] [CrossRef] [Google Scholar]
  165. Sun, H., Zhang, B., & Li, Z. 2015, ApJ, 812, 33 [NASA ADS] [CrossRef] [Google Scholar]
  166. Sun, H., Zhang, B., & Gao, H. 2017, ApJ, 835, 7 [NASA ADS] [CrossRef] [Google Scholar]
  167. Sun, H., Li, Y., Zhang, B.-B., et al. 2019, ApJ, 886, 129 [NASA ADS] [CrossRef] [Google Scholar]
  168. Tanaka, M., Utsumi, Y., Mazzali, P. A., et al. 2017, PASJ, 69, 102 [NASA ADS] [Google Scholar]
  169. Tanaka, M., Kato, D., Gaigalas, G., & Kawaguchi, K. 2020, MNRAS, 496, 1369 [NASA ADS] [CrossRef] [Google Scholar]
  170. Tang, C.-H., Huang, Y.-F., Geng, J.-J., & Zhang, Z.-B. 2019, ApJS, 245, 1 [NASA ADS] [CrossRef] [Google Scholar]
  171. Tanvir, N. R., Levan, A. J., Fruchter, A. S., et al. 2013, Nature, 500, 547 [CrossRef] [Google Scholar]
  172. Tanvir, N. R., Levan, A. J., González-Fernández, C., et al. 2017, ApJ, 848, L27 [CrossRef] [Google Scholar]
  173. Troja, E., Piro, L., van Eerten, H., et al. 2017, Nature, 551, 71 [NASA ADS] [CrossRef] [Google Scholar]
  174. Troja, E., Castro-Tirado, A. J., & Becerra González, J. 2019, MNRAS, 489, 2104 [NASA ADS] [Google Scholar]
  175. Troja, E., van Eerten, H., Zhang, B., et al. 2020, MNRAS, 498, 5643 [NASA ADS] [CrossRef] [Google Scholar]
  176. Troja, E., Fryer, C. L., O’Connor, B., et al. 2022, Nature, 612, 228 [NASA ADS] [CrossRef] [Google Scholar]
  177. Usov, V. V. 1992, Nature, 357, 472 [NASA ADS] [CrossRef] [Google Scholar]
  178. Valenti, S., Sand, D. J., Yang, S., et al. 2017, ApJ, 848, L24 [CrossRef] [Google Scholar]
  179. Villar, V. A., Guillochon, J., Berger, E., et al. 2017, ApJ, 851, L21 [Google Scholar]
  180. Wanderman, D., & Piran, T. 2015, MNRAS, 448, 3026 [NASA ADS] [CrossRef] [Google Scholar]
  181. Waxman, E., & Katz, B. 2017, Handbook of Supernovae (Springer International Publishing AG), 967 [CrossRef] [Google Scholar]
  182. Xiao, D., Zhang, B.-B., & Dai, Z.-G. 2019, ApJ, 879, L7 [NASA ADS] [CrossRef] [Google Scholar]
  183. Xin, L. P., Wang, J., & Wei, J. Y. 2021, ATel, 14610, 1 [Google Scholar]
  184. Xue, Y. Q., Zheng, X. C., Li, Y., et al. 2019, Nature, 568, 198 [NASA ADS] [CrossRef] [Google Scholar]
  185. Yang, B., Jin, Z.-P., Li, X., et al. 2015, Nat. Commun., 6, 7323 [NASA ADS] [CrossRef] [Google Scholar]
  186. Yang, J., Ai, S., Zhang, B.-B., et al. 2022, Nature, 612, 232 [NASA ADS] [CrossRef] [Google Scholar]
  187. Yi, S. X., Dai, Z. G., Wu, X. F., & Wang, F. Y. 2014, arXiv e-prints [arXiv:1401.1601] [Google Scholar]
  188. Yu, Y.-W., Zhang, B., & Gao, H. 2013, ApJ, 776, L40 [CrossRef] [Google Scholar]
  189. Yu, Y.-W., Liu, L.-D., & Dai, Z.-G. 2018, ApJ, 861, 114 [NASA ADS] [CrossRef] [Google Scholar]
  190. Yuan, W., Zhang, C., Feng, H., et al. 2015, arXiv e-prints [arXiv:1506.07735] [Google Scholar]
  191. Zhang, B. 2013, ApJ, 763, L22 [NASA ADS] [CrossRef] [Google Scholar]
  192. Zhang, B. 2018, The Physics of Gamma-Ray Bursts (Cambridge: Cambridge Univeristy Press) [Google Scholar]
  193. Zhang, D., & Dai, Z. G. 2010, ApJ, 718, 841 [NASA ADS] [CrossRef] [Google Scholar]
  194. Zhang, B., & Mészáros, P. 2001, ApJ, 552, L35 [CrossRef] [Google Scholar]
  195. Zhang, W., Woosley, S. E., & Heger, A. 2004, ApJ, 608, 365 [NASA ADS] [CrossRef] [Google Scholar]
  196. Zhang, B., Fan, Y. Z., Dyks, J., et al. 2006, ApJ, 642, 354 [Google Scholar]
  197. Zhang, B.-B., Zhang, B., Sun, H., et al. 2018, Nat. Commun., 9, 447 [NASA ADS] [CrossRef] [Google Scholar]
  198. Zheng, X. C., Xue, Y. Q., Brandt, W. N., et al. 2017, ApJ, 849, 127 [NASA ADS] [CrossRef] [Google Scholar]
  199. Zrake, J., & MacFadyen, A. I. 2013, ApJ, 769, L29 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: the magnetar model in detail

Here, we consider the formalism from Yu et al. (2013) and Sun et al. (2017) to describe the physical model of the magnetar emission. A nutshell description was done in Sect. 2.

A.1. Free-zone emission

The total rotational energy reservoir of a millisecond magnetar is

E rot = 1 2 I Ω 2 , $$ \begin{aligned} E_{\rm rot}=\frac{1}{2}I\Omega ^2, \end{aligned} $$(A.1)

where I is the moment of inertia, Ω = 2π/P is the angular frequency of the magnetar, and P its period. The millisecond magnetar loses its rotational energy through electromagnetic (EM) and gravitational-wave (GW) radiation (from non-axisymmetric deformation; Shapiro & Teukolsky 1983; Usov 1992; Zhang & Mészáros 2001; Gao et al. 2016; Lasky & Glampedakis 2016; Sun et al. 2017) as:

E ˙ rot = I Ω Ω ˙ = E ˙ EM + E ˙ GW = B p 2 R M 6 Ω 4 6 c 3 32 G I 2 ϵ 2 Ω 6 5 c 5 , $$ \begin{aligned} \dot{E}_{\rm rot}=I\Omega \dot{\Omega }=\dot{E}_{\rm EM} +\dot{E}_{\rm GW}=-\frac{B_p^2R_{\rm M}^6\Omega ^4}{6c^3} -\frac{32GI^2\epsilon ^2\Omega ^6}{5c^5}, \end{aligned} $$(A.2)

where Ω ˙ $ \dot{\Omega} $ is the time derivative of the angular frequency, BP is the dipole component of the magnetic field at the poles, ϵ is the ellipticity of the magnetar, RM is the magnetar radius, and c and G are the speed of light and Newton gravitational constant, respectively.

The evolution of Ω(t) considers the EM and GW losses, which both play a role in the spin-down evolution. In the following, we define the dipole spin-down luminosity as (Zhang & Mészáros 2001; Metzger & Piro 2014; Lasky & Glampedakis 2016; Sun et al. 2017; Zhang 2018):

L sd ( t ) = B P 2 R M 6 Ω ( t ) 4 6 c 3 , $$ \begin{aligned} L_{\rm sd}(t)=\frac{B_P^2R_{\rm M}^6\Omega (t)^4}{6c^3}, \end{aligned} $$(A.3)

where Ω(t) is the solution of the previous spin-down Equation A.2. Importantly, a strong magnetic field provides a mechanism for extracting rotational energy from the NS remnant via EM spin-down. Indeed, magneto-hydrodynamic simulations show that the magnetic fields in a BNS merger are amplified to values exceeding the field of Galactic magnetars (i.e., BP ∼ 1015 − 1016 G; Price & Rosswog 2006; Zrake & MacFadyen 2013; Kiuchi et al. 2014). Discarding one of the two terms in Eq. A.2, assuming the other is dominant, we can define the EM and GW characteristic spin-down timescales as τEM and τGW, respectively, (Zhang & Mészáros 2001):

τ EM = 3 c 3 I B P 2 R M 6 Ω i 2 = 3 c 3 I P i 2 B P 2 R M 6 ( 2 π ) 2 s , $$ \begin{aligned}&\tau _{\rm EM}=\frac{3c^3I}{B_P^2R_{\rm M}^6\Omega _i^2}=\frac{3c^3IP_i^2}{B_P^2R_{\rm M}^6(2\pi )^2}\,\mathrm{s}, \end{aligned} $$(A.4)

τ GW = 5 c 5 128 G I ϵ 2 Ω i 4 = 5 c 5 P i 4 128 G I ϵ 2 ( 2 π ) 4 s , $$ \begin{aligned}&\tau _{\rm GW}=\frac{5c^5}{128GI\epsilon ^2\Omega _i^4}=\frac{5c^5P_i^4}{128GI\epsilon ^2(2\pi )^4}\,\mathrm{s}, \end{aligned} $$(A.5)

where Ωi and Pi are the initial angular velocity and period of the magnetar, respectively. On times shorter than the spin-down timescale (called tsd, and defined as tsd = min[τEM, τGW]), it is expected that Ω(t) and Lsd(t) will remain stable. Considering the EM spin-down term only in Eq. A.2, the luminosity evolves as

L sd ( t ) = L 0 EM ( 1 + t τ EM ) 2 erg / s , $$ \begin{aligned} L_{\rm sd}(t)=L_{0}^\mathrm{EM}\left(1+\frac{t}{\tau _{\rm EM}}\right)^{-2}\,\mathrm{erg/s}, \end{aligned} $$(A.6)

where L 0 EM = I Ω i 2 ( 2 τ EM ) 1 $ L_{0}^{\mathrm{EM}}{=}I\Omega_i^2(2\tau_{\mathrm{EM}})^{-1} $ erg s−1. Eq. A.6 implies that the light curve shows a plateau whereby L sd L 0 EM $ L_{\mathrm{sd}}{\propto}L_{0}^{\mathrm{EM}} $ at t ≲ τEM, followed by a decay L sd L 0 EM ( t / τ EM ) 2 $ L_{\mathrm{sd}}{\propto}L_{0}^{\mathrm{EM}}(t/\tau_{\mathrm{EM}})^{-2} $ at t ≫ τEM.

The above equations, however, neglect the enhanced angular momentum losses due to neutrino-driven mass loss, which are important at early times during the magnetar formation (Metzger et al. 2011). On the other hand, considering only the GW spin-down term in Eq. A.2 (i.e., the angular velocity evolution is mainly due to the GW losses), the spin-down luminosity behaves as

L sd ( t ) = L 0 EM ( 1 + t τ GW ) 1 erg / s . $$ \begin{aligned} L_{\rm sd}(t)=L_{0}^\mathrm{EM}\left(1+\frac{t}{\tau _{\rm GW}}\right)^{-1}\,\mathrm{erg/s}. \end{aligned} $$(A.7)

Similar to Eq. A.6, the light curve will show a plateau L sd L 0 EM $ L_{\mathrm{sd}}{\propto}L_{0}^{\mathrm{EM}} $ at t ≲ τGW, followed by a decay as L sd L 0 EM ( t / τ GW ) 1 $ L_{\mathrm{sd}}{\propto}L_{0}^{\mathrm{EM}}(t/\tau_{\mathrm{GW}})^{-1} $ at t ≫ τGW. Such rapid rotation can produce large deformations, which imply significant quadrupole GW radiation, assuming rigid body rotation. It is important to note that the luminosity decays at late times as Lsd(t)∝t−1 (Eq. A.7), rather than Lsd(t)∝t−2 (Eq. A.6) like in the EM dipole spin-down scenario.

In reality, if GW radiation is to dominate spin-down, it must do so early in the evolution, while EM emission becomes important later (Lasky & Glampedakis 2016). The conditions for this come from the transition timescale once both the EM and GW terms are considered, and the spin-down luminosity behaves as

L sd ( t ) { t 0 t < t sd = τ GW t 1 t sd = τ GW < t < τ EM t 2 t > τ EM , $$ \begin{aligned} L_{\rm sd}(t){\propto }\left\{ \begin{matrix} t^0&t{<}t_{\rm sd}=\tau _{\rm GW} \\ t^{-1}&t_{\rm sd}=\tau _{\rm GW}{ < }t{ < }\tau _{\rm EM} \\ t^{-2}&t{>}\tau _{\rm EM} \end{matrix} \right., \end{aligned} $$(A.8)

for τGW < τEM, and

L sd ( t ) { t 0 t < t sd = τ EM t 2 t > t sd , $$ \begin{aligned} L_{\rm sd}(t){\propto }\left\{ \begin{matrix} t^0&t{ < }t_{\rm sd}=\tau _{\rm EM} \\ t^{-2}&t{>}t_{\rm sd} \end{matrix} \right., \end{aligned} $$(A.9)

for τEM < τGW. It remains unclear which is the preferential regime to lose rotational kinetic energy. To acquire evidence, it is necessary to detect a substantial number of new GW events and their associated EM signals using GW detectors and high-energy instruments, respectively, in the future (Lu & Kumar 2018).

A.2. Trapped-zone emission

The total energy of the ejecta (excluding the rest mass energy) can be expressed as

E ej = ( Γ 1 ) M ej c 2 + Γ E int , $$ \begin{aligned} E_{\rm ej}=(\Gamma -1)M_{\rm ej}c^2+\Gamma E_{\rm int}^{\prime }, \end{aligned} $$(A.10)

where E int $ E_{\mathrm{int}}^\prime $ is the internal energy in the co-moving frame9. The first and second terms represent the kinetic minus the rest-mass energy and the internal energy of the ejecta, respectively. At each time step, dt, the ejecta material is boosted by the energy of the magnetar (Lsd) and the radioactive decay of the r-process material (Lra). Moreover, the bolometric luminosity of the heated electrons is a lost energy source (Le) in the system. From energy conservation, we can conclude that in the observer frame

d E ej = ( L sd + L ra L e ) d t . $$ \begin{aligned} dE_{\rm ej}=(L_{\rm sd}+L_{\rm ra}-L_e)dt. \end{aligned} $$(A.11)

The conversion from observer time (t) to co-moving time (t′) is related as dt′=𝒟dt, where 𝒟 = 1/[Γ(1 − β cos θ)] is the Doppler factor with β = 1 Γ 2 $ \beta{=}\sqrt{1-\Gamma^{-2}} $ and θ = 0° for an on-beam observer. Deriving Eq. A.10 in terms of dt and equating to Eq. A.11, we obtain the evolution of the Lorentz factor as a function of the observer time as

d Γ dt = L sd + L ra L e Γ D ( d E int / d t ) M ej c 2 + E int . $$ \begin{aligned} \frac{d\Gamma }{dt}=\frac{L_{\rm sd}+L_{\rm ra}-L_e-\Gamma \mathcal{D} (dE_{\rm int}^{\prime }/dt^{\prime })}{M_{\rm ej}c^2+E_{\rm int}^{\prime }}. \end{aligned} $$(A.12)

In the last equation, we can see that the evolution of Γ also depends on the change of the internal energy of the ejecta material in the co-moving frame ( d E int / d t $ dE_{\mathrm{int}}^\prime/dt^\prime $). To convert the luminosities from the observer to the co-moving frames, we use the relation L′=L/𝒟2 (Zhang 2018).

The radioactive luminosity in the co-moving frame is described as (Korobkin et al. 2012)

L ra = 4 × 10 51 ( M ej M ) × [ 1 2 1 π arctan ( t t 0 t σ ) ] 1.3 erg s 1 , $$ \begin{aligned} L_{\rm ra}^{\prime }=4\times 10^{51}\left(\frac{M_{\rm ej}}{M_\odot }\right)\times \left[\frac{1}{2}-\frac{1}{\pi }\arctan \left(\frac{t^{\prime }-t_0^{\prime }}{t_\sigma ^{\prime }}\right)\right]^{1.3}\,\mathrm{erg}\,\mathrm{s}^{-1}, \end{aligned} $$(A.13)

where t 0 ~1.3 $ t_0^{\prime}{\sim}1.3 $ seconds and t σ ~0.11 $ t_\sigma^{\prime}{\sim}0.11 $ seconds are the best-fitted parameters of the nuclear heating power as a function of time from simulations of heavy element nucleosynthesis in compact binary mergers (Korobkin et al. 2012). In the case of a relativistic gas, the pressure (dominated by radiation) in the co-moving frame can be written as a function of the co-moving internal energy density as

P = 1 3 E int V . $$ \begin{aligned} P^{\prime }=\frac{1}{3}\frac{E_{\rm int}^{\prime }}{V^{\prime }}. \end{aligned} $$(A.14)

The co-moving volume evolution is determined by

d V d t = 4 π R ej 2 β c , $$ \begin{aligned} \frac{dV^{\prime }}{dt^{\prime }}=4\pi R_{\rm ej}^2\beta c, \end{aligned} $$(A.15)

while the radius evolution of the ejecta material in the observer frame can be described as

d R ej dt = β c 1 β · $$ \begin{aligned} \frac{dR_{\rm ej}}{dt}=\frac{\beta c}{1-\beta }\cdot \end{aligned} $$(A.16)

Additionally, as we explain above, the ejecta suffers a cooling via a bolometric electron emission. The co-moving frame radiated bolometric luminosity, which depends on the ejecta optical depth and the diffusion approximation, can be estimated as (Kasen & Bildsten 2010; Kotera et al. 2013)

L e = { E int c Γ τ R ej for t < t τ E int c Γ R ej for t t τ , $$ \begin{aligned} L_e^{\prime }=\left\{ \begin{array}{ll} \frac{E_{\mathrm{int}}^{\prime } c\Gamma }{\tau R_{\mathrm{ej}}}&\mathrm{for}\; t{ < }t_\tau \\ \frac{E_{\mathrm{int}}^{\prime } c\Gamma }{R_{\mathrm{ej}}}&\mathrm{for}\; t{\gtrsim }t_\tau \\ \end{array} \right., \end{aligned} $$(A.17)

where the term for t < tτ takes into account the skin-depth effect of an optically thick emitter (Yu et al. 2013), and tτ is the time when τ = 1.

The spectrum of the merger nova should approximately resemble a blackbody under the co-moving temperature, T′, and the peak energy of the blackbody emission is:

ε γ , p 4 D k b T = { 4 D k b ( E int a V τ ) 1 / 4 for τ > 1 4 D k b ( E int a V ) 1 / 4 for τ 1 , $$ \begin{aligned} \varepsilon _{\gamma ,p}\approx 4\mathcal{D} k_bT^{\prime }=\left\{ \begin{array}{ll} 4\mathcal{D} k_b\left(\frac{E_{\mathrm{int}}^{\prime }}{aV^{\prime }\tau }\right)^{1/4}&\mathrm{for}\; \tau {>}1 \\ 4\mathcal{D} k_b\left(\frac{E_{\mathrm{int}}^{\prime }}{aV^{\prime }}\right)^{1/4}&\mathrm{for}\; \tau {\lesssim }1 \\ \end{array} \right., \end{aligned} $$(A.18)

where kb and a are the Boltzmann and the blackbody radiation constants, respectively.

Appendix B: role of redshift and efficiency factor

As mentioned above, the redshift uncertainties affect the derived magnetar parameters. We recomputed the magnetar parameters adopting a range of redshifts (from z = 0.2 to 1.5) to explore this change. In Figure B.1, we can see how the magnetic field and the initial period move over a wide redshift range. Under a redshift modification from z = 0.2 to 1.0, both the initial period and magnetic field decrease by a factor of ≈7 and 5, respectively (see Fig. B.1).

thumbnail Fig. B.1.

This graph depicts the best-fit magnetar parameters at different redshifts and efficiencies. The number beside each data point is the corresponding redshift, while the colors represent efficiencies of η ∼ 10−3 (red lines), 10−2 (blue lines), and 10−1 (magenta lines). The gray region represents the breakup limit of a neutron star, which is ≈0.96 ms (Lattimer & Prakash 2004).

Moreover, during the fitting process (see Sect. 4), we fix the efficiency factor to η ∼ 10−3 based on literature outcomes (Xue et al. 2019). Here, we analyze its role in the magnetar parameters modifying the efficiency factor in Fig. B.1 (from η = 10−3 to 10−1). It is clear that the magnetic field and the initial rotational period increase by a factor of ≈5 and 7, respectively, where η is increased from 10−3 to 10−1. A similar factor has been found in previous works (e.g., Rowlinson et al. 2010, 2013).

Appendix C: posterior distributions for the free and trapped zones

This appendix shows the posterior distributions (from Fig. C.1 to C.10) of the free and trapped model fitting.

thumbnail Fig. C.1.

MCMC posterior distribution for the free zone fitting result of FXT 7 assuming znominal = 0.5 (left panel) and 1.0 (right panel). The dashed black lines show the 16%, 50% and 84% percentile of the MCMC results, while the red lines depicts the parameters of the best-fitted model.

thumbnail Fig. C.2.

Similar to Fig C.1, but for FXT 10.

thumbnail Fig. C.3.

Similar to Fig C.1, but for FXT 12.

thumbnail Fig. C.4.

MCMC posterior distribution for the free zone fitting result of FXTs 16 (left panel) and 19 (right panel).

thumbnail Fig. C.5.

MCMC posterior distribution for the free zone fitting result of FXT 22 assuming zcW = 1.04 (left panel) and zcNE = 1.5105 (right panel).

thumbnail Fig. C.6.

MCMC posterior distribution for the trapped zone fitting result of FXT 14. The dashed black lines show the 16%, 50% and 84% percentile of the MCMC results, while the red lines depict the parameters of the best-fitted model.

thumbnail Fig. C.7.

Similar to Fig C.6, but for FXT 20 assuming z = 1.0.

thumbnail Fig. C.8.

Similar to Fig C.6, but for FXT 20 assuming z = 0.5.

thumbnail Fig. C.9.

Similar to Fig C.6, but for FXT 15 assuming z = 1.0.

thumbnail Fig. C.10.

Similar to Fig C.6, but for FXT 15 assuming z = 0.5.

All Tables

Table 1.

Parameters and uniform priors distributions considered for the free and trapped magnetar models.

Table 2.

Parameters obtained from the fitting process of the magnetar model.

All Figures

thumbnail Fig. 1.

Cartoon illustrating the jet, free, and trapped zones, according to the magnetar model for FXTs. Left panel: situation at time t < tτ. Right panel: situation at time t > tτ. The insets in both panels show the light curves expected for the free and trapped zones before and after tτ.

In the text
thumbnail Fig. 2.

Theoretical light curves in the free-zone using Eqs. (1), (2) and (A.2). We explore how the light curve changes as we vary the magnetic field (top-left panel), the initial spin period (top-right panel), the magnetar radius (bottom-left panel), and the ellipticity (bottom-right panel). The text in the upper right of each panel shows the fixed parameters in the numerical solutions, while the legends and colored curves indicate the varied parameters and the range in the values.

In the text
thumbnail Fig. 3.

Theoretical results of the trapped-zone emission. Left panel: evolution of Lorentz factor (Γ). Middle panel: logarithmic optical depth (τ). The black horizontal line represents the optical depth limit at τ = 1. Right panel: X-ray luminosity (LX) of the magnetar model (Eq. (8), first term). The solid black line shows the free zone spin-down luminosity (ηLsd). The color code shows the different ejecta masses (Mej = 2 × 10−4 − 5 × 10−3M), while the solid and dashed lines represent opacities of 1.0 and 10.0 cm2 g−1, respectively, expected for the blue and red kilonova (Tanaka et al. 2017; Metzger 2019; Tanaka et al. 2020). The initial conditional assumed are: the initial velocity and internal ejecta energy of the ejecta are βi = 0.1, and E int,i E k =0.5 M ej β i 2 c 2 $ E_{\rm int,i} \simeq E_{\rm k} = 0.5\,M_{\rm ej}\beta_{\rm i}^2c^2 $, respectively. The magnetar parameters are: Bp = 5 × 1015 G, Pi = 10−3 s, I = 1045 g cm2, ϵ = 10−5, and ξ = 0.3.

In the text
thumbnail Fig. 4.

Evolution of the thermal component and light curves of the merger-nova emission (Eq. (7)) for different ejecta masses, color-coded from 10−4 to 5 × 10−3M. The top-left panel is the temperature evolution in the observer frame, while the other three panels represent the merger-nova emission at different frequencies (i.e., : 1 keV, 30 eV, and 1 eV). The model fixed parameters are ξ = 0.3, η = 10−3, κ = 1.0 cm2 g−1, Bp = 3 × 1015 G, and Pi = 5 ms.

In the text
thumbnail Fig. 5.

Theoretical light curves of the merger-nova emission, across several common optical/NIR filters, considering opacities of κ = 0.1 (left panels) and 10 cm2 g−1 (right panels), representing the lanthanide-free and -rich ejecta cases, respectively. KN emission without any magnetar spin-down enhancement is shown for comparison (dotted lines). The model fixed parameters are Pi = 5 ms, ξ = 0.3, and η = 10−3, along with two adopted ejecta masses of Mej = 10−4 (solid lines) and 10−2M (dashed lines). The upper and lower panels adopt magnetic field strengths of Bp = 1015 G and 1016 G, respectively.

In the text
thumbnail Fig. 6.

Free-zone fitting (solving the system of Eqs. (2) and (A.2)) in the source rest-frame for FXTs with known and nominal redshifts (see Table 2). Time dilation effects correct the data, and the bottom panels show the ratio between the data and the model (D/M).

In the text
thumbnail Fig. 7.

Trapped zone light-curve fitting of luminosities (solving the system of Eqs. (4)–(8)) in the source rest frame for FXTs with known and nominal redshifts (see Table 2). The bottom panels show the ratio between the data and the model values (D/M). The dashed lines represent the unobscured magnetar emission.

In the text
thumbnail Fig. 8.

Comparison of the derived magnetar parameters in this work (BP and Pi) with LGRBs (green triangles; Lyons et al. 2010; Yi et al. 2014), SGRBs (black squares; Rowlinson et al. 2010, 2013; Lü et al. 2015) and SLSNe (blue circles; Nicholl et al. 2017). The gray region represents the breakup limit of a neutron star, which is around 0.96 ms (Lattimer & Prakash 2004). FXTs related to the free and trapped models are represented by star and diamond markers, respectively.

In the text
thumbnail Fig. 9.

Numerical solution of the merger-nova emission in different energy bands. Left panel: optical and NIR upper limits of FXT 14 and the expected merger-nova emission (from u- to K-bands) assuming their best fitting parameters (see Table 2) and z = 2.23 (solid lines). As a middle and right panels: optical and NIR upper limits of FXT 22, and the merger-nova emission considering the derived magnetar parameters (see Table 2), two redshifts of z = 1.0 (solid lines) and 0.1 (dashed lines), and assuming lanthanide-free (middle panel; κ = 0.1 cm2 g−1 and Mej = 10−4M) and rich (right panel; κ = 30 cm2 g−1 and Mej = 10−3M) scenarios. The upper limits of FXT 14 were taken from Bauer et al. (2017), while for FXT 22 from Andreoni et al. (2021a,b,c), Xin et al. (2021), Rossi et al. (2021), and Eappachen et al. (2023). The time is measured from the X-ray trigger detection of the FXTs.

In the text
thumbnail Fig. 10.

GW strain evolution with frequency for the FXTs analyzed in this work. Each shaded region represents the GW strain considering eight EoS (see Table 3) at different distances. The cyan, gray, and black lines are the sensitivity limits for aLIGO, ET, and CE GW telescopes, respectively. For the FXTs without redshifts (FXTs 7, 10, 12, 15, and 20), we assumed z = 0.5 and 1.0 (see Table 2).

In the text
thumbnail Fig. B.1.

This graph depicts the best-fit magnetar parameters at different redshifts and efficiencies. The number beside each data point is the corresponding redshift, while the colors represent efficiencies of η ∼ 10−3 (red lines), 10−2 (blue lines), and 10−1 (magenta lines). The gray region represents the breakup limit of a neutron star, which is ≈0.96 ms (Lattimer & Prakash 2004).

In the text
thumbnail Fig. C.1.

MCMC posterior distribution for the free zone fitting result of FXT 7 assuming znominal = 0.5 (left panel) and 1.0 (right panel). The dashed black lines show the 16%, 50% and 84% percentile of the MCMC results, while the red lines depicts the parameters of the best-fitted model.

In the text
thumbnail Fig. C.2.

Similar to Fig C.1, but for FXT 10.

In the text
thumbnail Fig. C.3.

Similar to Fig C.1, but for FXT 12.

In the text
thumbnail Fig. C.4.

MCMC posterior distribution for the free zone fitting result of FXTs 16 (left panel) and 19 (right panel).

In the text
thumbnail Fig. C.5.

MCMC posterior distribution for the free zone fitting result of FXT 22 assuming zcW = 1.04 (left panel) and zcNE = 1.5105 (right panel).

In the text
thumbnail Fig. C.6.

MCMC posterior distribution for the trapped zone fitting result of FXT 14. The dashed black lines show the 16%, 50% and 84% percentile of the MCMC results, while the red lines depict the parameters of the best-fitted model.

In the text
thumbnail Fig. C.7.

Similar to Fig C.6, but for FXT 20 assuming z = 1.0.

In the text
thumbnail Fig. C.8.

Similar to Fig C.6, but for FXT 20 assuming z = 0.5.

In the text
thumbnail Fig. C.9.

Similar to Fig C.6, but for FXT 15 assuming z = 1.0.

In the text
thumbnail Fig. C.10.

Similar to Fig C.6, but for FXT 15 assuming z = 0.5.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.