- Split View
-
Views
-
Cite
Cite
M. Del Santo, A. A. Nucita, G. Lodato, L. Manni, F. De Paolis, J. Farihi, G. De Cesare, A. Segreto, The puzzling source IGR J17361–4441 in NGC 6388: a possible planetary tidal disruption event, Monthly Notices of the Royal Astronomical Society, Volume 444, Issue 1, 11 October 2014, Pages 93–101, https://doi.org/10.1093/mnras/stu1436
- Share Icon Share
Abstract
On 2011 August 11, INTEGRAL discovered the hard X-ray source IGR J17361–4441 near the centre of the globular cluster NGC 6388. Follow-up observations with Chandra showed the position of the transient was inconsistent with the cluster dynamical centre, and thus not related to its possible intermediate mass black hole. The source showed a peculiar hard spectrum (Γ ≈ 0.8) and no evidence of QPOs, pulsations, type-I bursts, or radio emission. Based on its peak luminosity, IGR J17361–4441 was classified as a very faint X-ray transient, and most likely a low-mass X-ray binary. We re-analysed 200 d of Swift/XRT observations, covering the whole outburst of IGR J17361–4441 and find a t−5/3 trend evident in the light curve, and a thermal emission component that does not evolve significantly with time. We investigate whether this source could be a tidal disruption event, and for certain assumptions find an accretion efficiency ϵ ≈ 3.5 × 10−4(MCh/M) consistent with a massive white dwarf, and a disrupted minor body mass Mmb ≈ 1.9 × 1027(M/MCh) g in the terrestrial-icy planet regime. These numbers yield an inner disc temperature of the order kTin ≈ 0.04 keV, consistent with the blackbody temperature of kTin ≈ 0.08 keV estimated by spectral fitting. Although the density of white dwarfs and the number of free-floating planets are uncertain, we estimate the rate of planetary tidal disruptions in NGC 6388 to be in the range 3 × 10−6–3 × 10−4 yr−1. Averaged over the Milky Way globular clusters, the upper limit value corresponds to 0.05 yr−1, consistent with the observation of a single event by INTEGRAL and Swift.
INTRODUCTION
The globular cluster (GC) NGC 6388 is considered one of the best candidates to harbour an intermediate mass black hole (hereafter IMBH; Baumgardt et al. 2005) at its centre. Optical observations allowed a first estimate of the IMBH mass of ≃5700 M⊙ (Lanzoni et al. 2007). Under the hypothesis that the IMBH is accreting the surrounding matter, it is natural to expect signatures in the X-ray band. XMM–Newton and Chandra observations towards NGC 6388 showed that several X-ray sources exist in the GC centre (see Nucita et al. 2008, Cseh et al. 2010) so that a unique identification of the putative IMBH is missing.
Based on the correlation between the X-ray and radio flux from black holes (BH; Merloni et al. 2003), Maccarone (2004) realized that radio observations may be very useful to pinpoint faint objects in GC. Following this suggestion, Cseh et al. (2010) combined radio observations of the central region of NGC 6388 with the Chandra X-ray flux of the IMBH candidate, in order to put an upper limit (3σ) on the BH mass of ≃1500 M⊙. A similar result was obtained combining high spatial resolution and wide-field spectroscopy of more than 600 stars in the direction of NGC 6388 (Lanzoni et al. 2013).
A general excitement started on 2011 August 11 when the IBIS telescope (Ubertini et al. 2003) on-board the INTEGRAL satellite (Winkler et al. 2003) identified a new hard X-ray source, labelled IGR J17361–4441, close to the NGC 6388 centre (Gibaud et al. 2011). The discovery of this transient opened the possibility that we were witnessing the awakening of the IMBH.
In order to explore the nature of IGR J17361–4441, several observing campaigns, in particular in the X and radio domains (Bozzo et al. 2011), have been performed. The X-ray position obtained with Swift/XRT was consistent with the centre of NGC 6388 (Ferrigno et al. 2011a). However, Chandra follow-up revealed that IGR J17361–4441 was located ∼ 2.7 arcsec from the cluster centre and it was not consistent with the position of any known X-ray source (Pooley et al. 2011). The transient source in NGC 6388 was also the target of two XMM–Newton slew observations (see for details Nucita et al. 2012) which revealed a spectrum softer than that initially observed by XRT.
Based on the X-ray luminosity of 6–9 × 1035 erg s−1 estimated with XRT, Wijnands et al. (2011) classified IGR J17361–4441 as very faint X-ray transient (VFXT; see Wijnands et al. 2006 for a discussion about VFXTs). Considering the location in a GC, these authors proposed the system being a neutron star (NS) low-mass X-ray binary (LMXB). However, they noted that the observed photon index of Γ ∼ 0.6–1.0 is atypical for LMXB transients in this luminosity range, where normally the photon index is around 1.6–2.2 (e.g. Del Santo et al. 2007; Degenaar & Wijnands 2009).
After the discovery, Swift/XRT followed the evolution of the IGR J17361–4441 outburst for ≃200 d until 2011 November 5 when the source was no longer observable (Bozzo et al. 2012). We have reanalysed these XRT data (obtaining consistent results) and performed a deeper analysis of the light curve and spectral evolution.
We noted that a ∝ t−5/3 trend is clearly evident in the light curve and we started to investigate whether this source might be a tidal disruption event (TDE).
TDE have been observed over the years as X-ray (Komossa & Bade 1999; Esquej et al. 2007; Bloom et al. 2011; Cenko et al. 2012) and UV (Gezari et al. 2009, 2012) transients. In most cases, such outbursts follow from the disruption of a star by a supermassive black hole (SMBH), which results in a powerful super-Eddington flare with a duration of several weeks. Recently, it has also been shown that such TDEs can also produce jets, resulting also in radio emission (Bloom et al. 2011; Levan 2011; Zauderer et al. 2011). In a few cases, the luminosity of the outburst appeared to imply that the disrupted object was a relatively low-mass object, such as a planet or brown dwarf. This is the case, for example, of the giant X-Ray flare in NGC 5905 (Li, Narayan & Menou 2002), or, more recently, of IGR J12580+0134, a flare discovered by INTEGRAL in 2011 and associated with NGC 4845 (Nikołajuk & Walter 2013).
In both cases, the best fit to the light curve implied the tidal disruption of a giant planet by a SMBH. In one other case, Campana et al. (2011) had interpreted the GRB101225A as a Galactic TDE, where an asteroid had been accreted by an NS. Such interpretation has, however, been ruled out by a redshift measurement of the source (Levan et al. 2014), that demonstrates its extragalactic origin.
Still, the sudden accretion of rocky objects by compact objects in the Galaxy is probably not uncommon. For example, there is substantial evidence for tidal disruption of asteroids by white dwarfs (WDs). To date, there are 30 published cases where debris is seen via thermal emission in the infrared (e.g. Farihi et al. 2010), and sometimes also gaseous metallic emission in the optical (e.g. Gänsicke et al. 2006), and via accreted metals in the stellar atmosphere. These circumstellar discs lie completely within the Roche limit of the stars, approximately 1 R⊙, and are thus consistent with tidally disrupted planetesimals (Jura 2003). The planetary nature of the disc material is given both by the abundance ratios in the polluted stellar atmospheres (e.g. Gänsicke et al. 2012), and via infrared emission of evolved solids associated with planet formation (Jura, Farihi & Zuckerman 2009).
Among these polluted WDs, there appears to be evidence in support of intense accretion episodes of planetary debris (Farihi, Gänsicke & Koester 2013). Briefly, trace metals in an H-rich WD provide instantaneous measures of accretion rates, which are typically in the range 105–108 g s−1 (Koester & Wilken 2006). On the other hand, atmospheric metals in a comparable He-rich WD give a historical average of accretion over time-scales up to 106 yr (Farihi, Jura & Zuckerman 2009). Remarkably, while there are examples of stars with time-averaged accretion rates up to 1011 g s−1, no instantaneous accretion rate above 109 g s−1 has yet been witnessed among many dozens of metal polluted WDs. Analysis of available data suggest that there are short-lived bursts (thigh < 103 yr) of high-rate accretion in these polluted WD systems, such that the product of |$t_{\rm high}\dot{M}_{\rm high}$| is around 1024 g, or 1/6000 the mass of Earth (Farihi et al. 2012).
The structure of the paper is as follows. In Section 2, we briefly discuss the INTEGRAL and Swift data analysis. In Section 3, we analyse the spectral evolution and the light curve of the system, demonstrating a remarkable t−5/3 decline, typical of TDEs. In Section 4, we discuss the possible interpretation of the source: after ruling out a possible X-ray binary scenario, we concentrate on a TDE interpretation, showing that the parameter set that best describe the system indicates the disruption of a rocky object by a WD. We also discuss and rule out a possible interpretation as an extragalactic TDE in a background galaxy. In Section 5, we draw our conclusions.
OBSERVATION AND DATA ANALYSIS
Swift/XRT
We have analysed the Swift/XRT observations performed on the target IGR J17361–4441 (Id. numbers from 00032072001 to 00032072038; see also Bozzo et al. 2012) for a total exposure time of ≃105 ks. Swift data have been analysed by using the standard procedures described in Burrows et al. (2005) and with the latest calibration files.1 We processed the XRT data with the xrtpipeline (v.0.12.6) task and we applied standard screening criteria by using ftools (heasoft v.6.13.0). The source spectra have been extracted from a circular region centred on the target nominal coordinates (with radius of ≃ 40 arcsec) by using the xselect routine. The background spectra were obtained from circular regions far enough from the IGR J17361–4441 position. When necessary, we corrected observations affected by pile-up (i.e. with source count rates above ≃0.5 ct s−1). We then used the XRTMKARF task to create the ancillary response files. Of the 38 XRT spectra, we avoided to use the last six because of the low statistics. Each of the 32 spectrum has been rebinned in order to have at minimum 25 counts per channel with GRPPHA. All fits have been performed with xspec (Arnaud, Dorman & Gordon 2007) and errors calculated at 90 per cent confidence level.
INTEGRAL/IBIS
We have analysed all data collected by the low-energy IBIS detector, i.e. ISGRI (Lebrun et al. 2003), when IGR J17361–4441 was present in the total IBIS field of view (14| $_{.}^{\circ}$|5×14| $_{.}^{\circ}$|5). Observations have been performed in the period 2011 August 11–October 22, i.e. INTEGRAL revolutions 1078–1101 for a total of 454 science windows (SCW) and ≃1 Ms time exposure. Data analysis was carried out with the Off-line Scientific Analysis (osa) v. 10.0 distributed by the ISDC (Courvoisier et al. 2003). We processed data with the ibis_science_analysis task and extracted images in three energy bands: 18–40, 40–80, 80–150 keV. We did not get any detection at SCW level (∼2 ks), so that we mosaicked the 454 images in groups of 38, 39, 61, 63, 130 and 123. A significant signal to noise (above 5σ) has been obtained only in the first energy band (18–40 keV). The errors on the IBIS/ISGRI count rate have been obtained via the significance (signal-to-noise ratio) in output from the osa software. We extracted an averaged ISGRI spectrum (53 ks) corresponding to the flat part of the XRT light curve (see Section 3.2) by using the osa task MOSAIC_SPEC (useful for faint sources).
Swift/BAT
The Swift/BAT survey data collected from 2004 November up to end of 2013 retrieved from the HEASARC public archive were processed using batimager (Segreto et al. 2010). This code is dedicated to the processing of coded mask instrument data. We found that IGR J17361–4441 is detected in only one time interval, since MJD ∼ 55770 up to MJD ∼ 55828, with a statistical significance of 15 standard deviations in the 20–85 keV all-sky map. Light curve has been extracted with 5 d temporal bins.
RESULTS
Temporal analysis
We have estimated the 1-10 keV unabsorbed fluxes by fitting the 32 Swift/XRT spectra of IGR J17361–4441 with a simple power law. Plotting the flux evolution, we noted that it could be well represented by a broken power law. Thus, we fit the light curve with a constant plus an exponential law (see Fig. 1). It results a plateau flux Fmax = (4.1 ± 0.1) × 10−11 erg cm−2 s−1 and a knee time of tk = 36 ± 1 d from the INTEGRAL trigger MJDIBIS = 55784.9479, i.e. roughly 5 d before the first XRT follow-up observation. Then, it was very interesting to observe that after the knee, the data were well fitted with a ∝ t−5/3 law (|$\chi _{\nu }^{2}$| = 1.14 with 30 d.o.f.).
In Fig. 2, the INTEGRAL/IBIS (up) and Swift/BAT (bottom) hard X-ray light curves are shown. Thanks to BAT, we found that the IGR J17361–4441 outburst was started (in hard X-rays) roughly 14 d (MJDBAT = 55770.4480) before the INTEGRAL trigger, thus implying a total event (observed) duration of 99 d, when the last XRT point (MJDstop = 55869.8738 d) is considered as the end.
Spectral analysis
Extending the XRT spectral energy range down to 0.3 keV, an additional disc blackbody component to the power law is requested by the data. Then, the 32 XRT spectra of IGR J17361–4441 have been fit with an absorbed power-law plus DISKBB (in xspec). The column density is constant at a mean value of ≃0.8 × 1022 cm−2 (also consistent with Bozzo et al. 2011), so that it has been fixed in our fit procedure. The evolutions of the inner disc temperature (kTin) and power-law slope (Γ) are shown in Fig. 3. It is interesting to note that in spite of the variation of the power-law component, i.e. an abrupt change of the power-law slope at roughly tn ≃ 30 d from the INTEGRAL trigger is observed, the disc temperature results constant at a value around 0.08 keV (Fig. 3).
To obtain the hard X-ray spectrum up to ≃100 keV, we averaged the IBIS/ISGRI data corresponding to the first four bins of the light curve (Fig. 2, top). These data are simultaneous to the plateau period in the XRT light curve, when IGR J17361–4441 did not show any spectral variation (Fig. 3). Then, we have combined the IBIS/ISGRI spectrum with one of the brighter XRT spectra and fit the broad-band spectrum with an absorbed cut-off power-law (CUTOFFPL) plus DISKBB, obtaining parameters (see Table 1) consistent with those reported in Bozzo et al. (2011).
In addition, we fit this spectrum with a physical model (Fig. 4) by using a thermal Comptonization model (COMPTT in xspec) in place of the cut-off power law obtaining an electron population temperature of roughly 10 keV, an inner disc radius Rin ≃ 12 600 km (from the normalization of the DISKBB model; Mitsuda et al. 1984) and a bolometric (0.1–100 keV) flux of about 2 × 10−9 erg cm−2 s−1 (see Table 1). Assuming a distance of 13.2 kpc (Dalessandro et al. 2008), we have estimated the bolometric peak luminosity as Lbol ≃ 3.5 × 1037 erg s−1.
It is known that in accreting compact objects LX corresponds (as first approximation) to the accretion luminosity |$L_{\rm acc} = \epsilon \, \dot{M} \, c^2$|. Thus, a lower limit to the accreted mass on to the central object can be evaluated.
We have integrated the bolometric luminosity over the time for the whole outburst duration and obtained Macc ≃ 3.4 × 1023 ϵ− 1 g.
DISCUSSION: WHAT KIND OF SOURCE CAN IT BE?
X-ray binary scenario
Based on the peak luminosity of |$L_{\rm 2{\rm -}10\, keV}$| = 8.5 × 1035 erg s−1, IGR J17361–4441 has been classified as VFXT (Wijnands et al. 2011). This is a non-homogeneous class of sources characterized by a peak luminosity (in 2–10 keV) within the range 1034-36 erg s−1 (Wijnands et al. 2006). Very faint transients in GC are most likely LMXBs with NS or BH as accreting compact object.
As discussed in details in Kulkarni et al. (1993), Portegies Zwart et al. (2000), and Strader et al. (2012), one expects that even if hundreds of stellar-mass black holes form in a typical globular cluster, only a few of them survive since the others are likely ejected through dynamical interactions. According to Kalogera, King & Rasio (2004), they are likely to have extremely low duty cycles (below about 10−3), consistent with the absence of such an object in Galactic GC. Observational support for such picture was provided by the lack of X-ray binaries with BH in Milky Way GC until 2011, when Strader et al. (2012) proposed that two radio sources in the GC M22 are stellar-mass BH.
The possibility that IGR J17361–4441 is a binary system with a BH can be ruled out because of the low radio-to-X-ray flux ratio (see Bozzo et al. 2011; Chomiuk et al. 2013). Ferrigno et al. (2011b) reported on the ATCA non-detection of IGR J17361–4441; they inferred an upper limit of 19 μJy at 9 GHz which translates into a luminosity at 8.5 GHz of 4 × 1028 erg s−1. We note that this value is at least a factor of 100 fainter than the typical radio luminosity observed in BH binaries with LX ≃ 8 × 1035 erg s−1(see fig. 3 in Migliari & Fender 2006). Moreover, BH transients in the hard state usually show quasi-periodic oscillation (see Homan & Belloni 2005; McClintock & Remillard 2006). These features were also observed in the very faint system Swift J1357.2–0933 (Casares et al. 2011; Armas Padilla et al. 2014). However, the power density spectrum constructed after accumulating RXTE/PCA light curves of IGR J17361–4441 did not show the presence of any coherent signal (against white noise) in the frequency range 10−3–103 Hz at 3σ confidence level (see fig. 2 in Bozzo et al. 2011).
The other possibility is that IGR J17361–4441 is an LMXB with an NS accreting at very low |$\dot{M}$|, either with low magnetic field (possibly showing suddenly type-I X-ray bursts, i.e. burst-only sources; Cornelisse et al. 2004, Del Santo et al. 2007, Campana 2009) or with high magnetic field, i.e. an accreting millisecond X-ray pulsars (AMSP; Wijnands et al. 2008; Wijnands 2008).
However, the hard X-ray spectrum (Γ ≃ 0.8; Table 1) observed at the beginning of the outburst of IGR J17361–4441 is not observed in VFXT with low magnetic field NS. Usually, these objects show Γ ∼ 1.5–2.2 in this luminosity range (see Degenaar & Wijnands 2010 and references therein) and blackbody temperatures (when present) consistently higher than that observed in IGR J17361–4441 (i.e. 0.3–0.5 keV; see Armas Padilla, Degenaar & Wijnands 2013). Moreover, the soft spectra observed at the end of the outburst are typical of the soft spectral state, which is observed in NS burster sources only when in high accretion rate regime (i.e. luminosities; Barret et al. 2000).
Model . | NH . | kTin . | Ec . | Γ . | kTe . | τ . | |$\chi ^2_{\nu }$|(dof) . | Fbol . |
---|---|---|---|---|---|---|---|---|
. | 1022 (cm−2) . | (keV) . | (keV) . | . | (keV) . | . | . | (erg cm−2 s−1 ) . |
DISKBB+CUTOFFPL | 0.7 ± 0.1 | |$0.089^{+0.004}_{-0.005}$| | |$41^{+18}_{-12}$| | 0.8 ± 0.1 | - | - | 0.80(62) | 1.0 × 10−9 |
DISKBB+COMPTT | 0.7 ± 0.1 | |$0.087^{+0.004}_{-0.006}$| | - | - | 12 ± 2 | 8 ± 1 | 1.01(62) | 1.7 × 10−9 |
Model . | NH . | kTin . | Ec . | Γ . | kTe . | τ . | |$\chi ^2_{\nu }$|(dof) . | Fbol . |
---|---|---|---|---|---|---|---|---|
. | 1022 (cm−2) . | (keV) . | (keV) . | . | (keV) . | . | . | (erg cm−2 s−1 ) . |
DISKBB+CUTOFFPL | 0.7 ± 0.1 | |$0.089^{+0.004}_{-0.005}$| | |$41^{+18}_{-12}$| | 0.8 ± 0.1 | - | - | 0.80(62) | 1.0 × 10−9 |
DISKBB+COMPTT | 0.7 ± 0.1 | |$0.087^{+0.004}_{-0.006}$| | - | - | 12 ± 2 | 8 ± 1 | 1.01(62) | 1.7 × 10−9 |
Model . | NH . | kTin . | Ec . | Γ . | kTe . | τ . | |$\chi ^2_{\nu }$|(dof) . | Fbol . |
---|---|---|---|---|---|---|---|---|
. | 1022 (cm−2) . | (keV) . | (keV) . | . | (keV) . | . | . | (erg cm−2 s−1 ) . |
DISKBB+CUTOFFPL | 0.7 ± 0.1 | |$0.089^{+0.004}_{-0.005}$| | |$41^{+18}_{-12}$| | 0.8 ± 0.1 | - | - | 0.80(62) | 1.0 × 10−9 |
DISKBB+COMPTT | 0.7 ± 0.1 | |$0.087^{+0.004}_{-0.006}$| | - | - | 12 ± 2 | 8 ± 1 | 1.01(62) | 1.7 × 10−9 |
Model . | NH . | kTin . | Ec . | Γ . | kTe . | τ . | |$\chi ^2_{\nu }$|(dof) . | Fbol . |
---|---|---|---|---|---|---|---|---|
. | 1022 (cm−2) . | (keV) . | (keV) . | . | (keV) . | . | . | (erg cm−2 s−1 ) . |
DISKBB+CUTOFFPL | 0.7 ± 0.1 | |$0.089^{+0.004}_{-0.005}$| | |$41^{+18}_{-12}$| | 0.8 ± 0.1 | - | - | 0.80(62) | 1.0 × 10−9 |
DISKBB+COMPTT | 0.7 ± 0.1 | |$0.087^{+0.004}_{-0.006}$| | - | - | 12 ± 2 | 8 ± 1 | 1.01(62) | 1.7 × 10−9 |
The possibility that the system might harbour an AMSP was proposed at the onset of the Swift observations, hypothesis which was immediately in contrast with the lack of pulsations in the PCA data (Bozzo et al. 2011). Even considering that pulsation in LMXBs might be transient as observed in Aql X–1 (Casella et al. 2008), we note that AMSPs have similar power-law slopes than the non-pulsating NS and blackbody temperatures even higher, i.e. up to 1 keV (see Papitto et al. 2010; Patruno & Watts 2012). Moreover, they usually do not show variations of the spectral shape at different flux levels (see Falanga et al. 2012 and references therein).
We conclude that the hypothesis either of low magnetic (burster) or high magnetic (AMSP) NS for the compact object of IGR J17361–4441 is unlikely. In addition, this is supported by the fact that NS transients have been observed to show (usually) quiescent luminosities of about 1032-33 erg s−1(see Rea et al. 2011, and references therein), i.e. higher than the upper limit of few 1031 erg s−1 derived for IGR J17361–4441.
Hints for a TDE?
Based on the considerations reported in Section 4.1, we suggest IGR J17361–4441 being a TDE. There are two empirical evidences that point in the direction of interpreting this flare of IGR J17361–4441 as a TDE. First and foremost, the apparent decline of the light curve as t−5/3, that is typical of such events, even though we can fit the expected decline only to the late light curve, while the early light curve shows a plateau which is generally not expected in TDEs. We tried to fit the light curve with a pure t−5/3 type law and we did not obtain a good fit (|$\chi ^{2}_{\nu }$| = 2.3, 30 d.o.f., see Fig. 5).
Secondly, several TDEs produced by the disruption of stars by an SMBH appear to be characterized by a thermal emission component that does not evolve significantly with time, as in the present case. Such component has been observed in TDEs discovered in the optical (van Velzen et al. 2011; Gezari et al. 2012; Chornock et al. 2014) and in X-rays (Maksym et al. 2013; Donato et al. 2014). A slow evolution of the thermal component is not expected from theoretical models of TDE, that imply a cooling thermal component (Lodato & Rossi 2011), unless the fallback rate is strongly super-Eddington (Loeb & Ulmer 1997), or unless some reprocessing of the high-energy emission by an opaque shroud of debris is invoked (Guillochon, Manukian & Ramirez-Ruiz 2014, which would only be able to explain the events observed in the optical).
In a TDE, the accreted mass is half of the mass of the disrupted object. According to the numbers quoted above (end of Section 3), we thus estimate that the disrupted minor body would have a mass of Mmb ≈ 7 × 1023ϵ−1 g. Even if the efficiency parameter is low, if the event is caused by accretion on to a compact object it would put the minor body in the asteroid – terrestrial planet mass range, and would thus most likely be of rocky composition.
A critical question to be addressed is whether the disrupted material, that would initially orbit relatively far from the WD surface, would be able to accrete fast enough so that the accretion rate on to the WD matches the fallback rate given in equation (1). For this to happen, we need to require that the viscous time in the disc is much smaller than the fallback time tmin. Two recent studies have addressed this question. Metzger, Rafikov & Bochkarev (2012) argue that most of the accretion occurs at a late stage after the disc formation, and estimate viscous times of the order of thousands of years. In this case, the WD would accrete at a rate much smaller than the fallback rate and the disc would shine mostly in the infrared. Conversely, Bear & Soker (2013) conclude that the energy liberated during the formation of the disc is high enough to lead to a rapid accretion of the material, which would lead to an initial spike of intense accretion at rates comparable to the fallback rate. In this case, we do expect to see the system to emit in X-rays. Thus, in order to explain this source as the tidal disruption of a terrestrial planet by a WD, we need to invoke the mechanism suggested by Bear & Soker (2013).
The hard X-ray emission is interpreted to originate from inverse-Compton scattering of lower energy photons (kTin ≃ 0.09 keV) by high-energy electrons in a hot corona (kTe ≃ 12 keV; Table 1) forming around the accretion flow (see also Nikołajuk & Walter 2013). This is the same mechanism which is thought to occur in all compact objects accreting matter via an accretion disc. In TDE, we do not have detailed models for the formation and evolution of the corona and the associated time-scales. On the other hand, hard X-ray emission is not uncommon in TDE candidates (Burrows et al. 2011; Cenko et al. 2012; Nikołajuk & Walter 2013), but it is difficult to make a comparison with our IGR J17361–4441 event since they are all thought to have occurred on an SMBH.
Estimating the rate of planetary tidal disruption
In this section, we estimate the rate of tidal disruption of planets by WD in the core of a GC. It should be noted that these estimates rely on a number of assumptions that are either poorly constrained observationally or poorly understood theoretically, so they should be treated with caution.
Numerical simulations have demonstrated that dissolution of planetary systems occurs on time-scales that are shorter than the age of a typical open or GC (Spurzem et al. 2008). We thus expect that the disruption does not follow from the deflection of a planet that is initially bound to the WD.
However, the interacting object may be a free-floating planet (FFP), resulting from the dissolution of planetary systems. The retention of a population of FFP in GC is not obvious. Simple estimates based on equipartition arguments would imply that a planetary population would acquire an extremely large velocity dispersion, well in excess of the escape speed from the cluster, and would thus evaporate within a few relaxation times (note that the relaxation time for NGC 6388 is ≈5.2 × 107 yr at the core and ≈8 × 108 yr at the half-mass radius; Harris 1996). However, numerical simulations (Fregeau et al. 2002; Hurlay & Shara 2002) show that a significant population of FFP is retained in relaxed systems. A possible explanation for this behaviour is due to a lack of equipartition. Indeed, again through numerical simulations including a mass spectrum, Trenti & van der Marel (2013) have shown that, as a consequence of Spitzer instability, equipartition does not occur even for relaxed stellar systems, and that the velocity dispersion depends very weakly on mass, such that σ ∝ m−0.08. In our particular case, with a mass of a third Earth mass and an average stellar mass of 0.1 M⊙, we estimate that σpl ≈ 2.5σstar.
In order to estimate the rate of planetary disruptions, we need to know the number density of FFP and WD in a typical GC, both of which are very uncertain. Here, we provide a pessimistic estimate of the rate of disruptions by considering the lowest values for both parameters, and an optimistic estimate considering the highest values. For the WD number density, we assume either a low value of ≈104 pc−3 (Raskin et al. 2009), or a high value – more appropriate for collapsed GC, such as NGC 6388 – of 105 pc−3, in agreement with simulations performed by Ivanova et al. (2006). Estimating the number density of FFP of rocky composition is even more difficult. It has been estimated that the density of FFP in GC could exceed NFFP = 106 pc−3 (Soker et al. 2001), i.e. approximatively one per star in NGC 6388 considering the star density reported in Pryor & Meylan (1993). Indeed, N-body simulations suggested that FFPs could evolve to the current epoch with a population that exceeds the stellar population at the cluster centre by a factor of ∼100 (Fregeau et al. 2002; Hurlay & Shara 2002). This is also supported by N-body simulation of planet formation performed by Ida, Kokubo & Kominami (2003). We thus assume a low value for NFFP = 107 pc−3 (i.e. ten per star) and a high value of 108 pc−3 (i.e. 100 per star).
In the HST/ACS sources catalogue of NGC 6388 used by Lanzoni et al. (2007), there are 31 point sources within the Chandra error box of IGR J17361–4441. No one is a WD (Lanzoni, private communication). A more detailed analysis of the data available in the HST archive will be presented in a forthcoming paper. However, it is worth to note that only very bright WDs can be observed by HST in the redden core of NGC 6388.
Extragalactic TDE
We considered the possibility that IGR J17361–4441 has an extragalactic origin, i.e. it is due by a sudden increase in the accretion rate on an SMBH at the centre of an otherwise quiescent galaxy (see Komossa et al. 2004; Burrows et al. 2011; Saxton et al. 2012; Nikołajuk & Walter 2013).
First of all, we have estimated the number of background AGNs expected towards the target within a circle of radius ≃ 0.6 arcsec through the log N–log S diagram (Hasinger, Miyaji & Schmidt 2005) and by using the Chandra flux upper limit ≃4.8 × 10−16 erg cm−2 s−1 (Pooley et al. 2011). Assuming a power-law model with Γ = 1.7 and an absorption column density NH = 3 × 1021 cm−2 (Dickey & Lockmann 1990), we have obtained an unabsorbed flux in the 0.5–2 keV band of 2.9 × 10−16 erg cm−2 s−1 (by using webpimms v3.92). The expected number of background AGNs is then ≃1 × 10−4. Thus, the probability of the existence of a background AGN is unlikely although not negligible.
On the other hand, we have estimated the likelihood that the event is an extragalactic TDE occurring in a background galaxy, either with an AGN activity or not. A decay time-scale of ≈70 d is in agreement with the expected time-scale for the disruption of a solar-type star by an SMBH. In particular, we can reproduce the observed inner disc temperature at peak of kTin ≈ 0.08 keV and a decay time-scale of 68 d if the disrupted object has a mass M* = 0.48 M⊙ and the SMBH has a mass of M ≈ 5.8 × 106 M⊙.
We can then use the Marconi & Hunt (2003) LK, bul-MBH relation to obtain the K-band luminosity of the putative host galaxy for the assumed BH mass of ∼5.8 × 106 M⊙. At a distance of ∼150 Mpc, the resulting apparent luminosity magnitude would be of the order of 15. Now, such a bright and relatively nearby galaxy would be seen as extended object by HST/ACS. Indeed, the limiting magnitude for bright galaxies (above which they are observed as point sources) is around 23.3 (see Blanton & Moustakas 2009; Bono et al. 2010 and references therein). No extended source has been observed in the field of IGR J17361–4441 by HST and we can thus rule out the possibility of an extragalactic TDE.
CONCLUSION
The nature of the puzzling new INTEGRAL source IGR J17361–4441 is still unknown after a number of follow-up observations. We noted two empirical evidences typical of TDEs: a ∝ t−5/3 trend of the light curve and a thermal disc blackbody emission almost constant with time. Thus, we started to investigate the possibility that a minor body was tidally disrupted by a compact object. Our observational analysis and theoretical calculations indicate a free-floating terrestrial planet tidally disrupted by a WD. In particular, the mass of the disrupted body is Mmb ≈ 1.9 × 1027(M/MCh) g, while the accreting object is consistent with a WD close to the Chandrasekhar limit.
Finally, we have estimated the rate of a planetary disruption event by a WD in the GC NGC 6388: although many of the relevant parameters are highly uncertain, we estimate the rate |$\dot{N}_{\rm TE}$| to be in the range ≃10−6–10−4 yr−1. The lower limit suggests that the detection of such an event in the Milky Way is unlikely. However, if we consider the optimistic estimate for the rate, and considering that the total number of GC in the Galaxy is roughly 150, we obtain a total rate of events of about 0.05 yr−1, i.e. one every ∼20 yr, comparable with the lifetime of INTEGRAL and Swift.
We acknowledge interesting discussions with Giuseppe Bertin, Piergiorgio Casella, Domitilla De Martino, Alessandra De Rosa, Immacolata Donnarumma, Giuliana Fiorentino, Barbara Lanzoni, Cristina Pallanca, Patrizia Romano, Alexander Tschekoskoy, Dimitri Veras, Andrea Tiengo, Alice Zocchi. MDS and GDC acknowledge financial support from the agreement ASI-INAF I/009/10/0 and from PRIN-INAF 2009 (PI: L. Sidoli). GL acknowledges financial support from PRIN MIUR 2010–2011, ‘The Chemical and Dynamical Evolution of the Milky Way and Local Group Galaxies’, prot. 2010LY5N2T. MDS thanks the Mathematics and Physics Department at the University of Salento for hospitality. AN thanks INAF/IAPS for hospitality.