Open Access
Issue
A&A
Volume 669, January 2023
Article Number A66
Number of page(s) 10
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/202245015
Published online 10 January 2023

© F. Haberl et al. 2023

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.

Open access funding provided by Max Planck Society.

1. Introduction

Type-I X-ray bursts are observed from low-mass X-ray binaries (LMXBs) and they uniquely identify the compact object in the system as a neutron star. The bursts occur when a sufficiently large amount of matter rich in H and He, which is accreted from the companion star, is gathered on the surface of the neutron star to ignite a thermonuclear runaway. Depending on the accretion rate, most bursts repeat on timescales of hours to days and typically last from ∼10 s to minutes. A few bursts were observed with much longer durations between several hours and a day. These very rare superbursts with a long recurrence time are thought to burn carbon, rather than the H-He fuel for short bursts (for reviews, see Galloway & Keek 2021; Parikh et al. 2013; Strohmayer & Bildsten 2006; Lewin et al. 1993). Temporal and spectral properties of a large sample of more than 7000 bursts from 85 bursting sources are described in Galloway et al. (2020).

The first type-I X-ray bursts outside of the Milky Way were detected in the Local Group galaxy M 31 using XMM-Newton data (Pietsch & Haberl 2005). The authors found two bursts, which were significantly detected in at least two of the EPIC instruments (Strüder et al. 2001; Turner et al. 2001). In a systematic search for fast transients in XMM-Newton EPIC-pn observations, Pastor-Marazuela et al. (2020) reported the identification of four bursts in the direction of M 31. The first burst in their Table 6 is identical to the second burst reported by Pietsch & Haberl (2005)1. The other three bursts are new, thus increasing the number of X-ray bursters known in M 31 to five.

In the course of our X-ray source population studies of the Large and Small Magellanic Cloud (LMC and SMC) with eROSITA, we also monitored the Magellanic Bridge between the two satellite galaxies of the Milky Way. Together with the Magellanic Clouds, the Magellanic Bridge is part of a gaseous structure around the Magellanic Clouds spanning 200 degrees on the sky in total, also including a leading arm and a trailing stream (Nidever et al. 2010). The Magellanic Bridge consists of neutral gas and a young stellar component (Irwin et al. 1990). Moreover, Bagheri et al. (2013) found evidence for an older stellar population. Such a population is expected if the bridge was formed by tidal interactions between the LMC and SMC, which stripped gas and stars from the clouds (Mathewson 1985).

During the third 6 month long all-sky survey (eRASS3), the eROSITA instrument (Predehl et al. 2021) on board the Russian/German Spektrum-Roentgen-Gamma (SRG) mission (Sunyaev et al. 2021) discovered a new hard X-ray transient, designated eRASSt J040515.6 − 745202 (Rau et al. 2021). The source is located in the direction of the Magellanic Bridge, was not detected in eRASS1 nor eRASS2, and had faded again below the detection limit in eRASS4. Here, we report on X-ray and UV follow-up observations of eRASSt J040515.6 − 745202 with XMM-Newton, Swift, and NICER. An optical identification is provided from Dark Energy Camera (DECam) archival exposures.

2. X-ray observations

2.1. eROSITA

After its launch in July 2019, eROSITA started to scan the sky in December 2019 and completed four full sky surveys (eRASS1 to eRASS4) by December 2021. eROSITA consists of seven co-aligned Wolter type-I telescopes, each equipped with a camera based on a charge coupled device (CCD) of pn type with an integration time of 50 ms (Predehl et al. 2021).

eROSITA scanned eRASSt J040515.6 − 745202 eleven times during eRASS3 between 2021 May 1, 01:51 (MJD 59335.07738) and 2021 May 2, 17:52 (MJD 59336.74441; see Fig. 1). For the data analysis, we used the eROSITA Standard Analysis Software System (eSASS version eSASSusers_211214_0_3; Brunner et al. 2022). To extract source and background events for light curves and spectra, we used the eSASS task srctool (see e.g. Maitra et al. 2021; Haberl et al. 2022). We used circular regions with radii of 90″ and 120″ around the position of the source and a nearby source-free region, respectively, and selected all valid pixel patterns (PATTERN = 15). For the light curves, we combined the data from all cameras (telescope modules TM 1–7) and applied a cut in the fractional exposure of 0.15 (FRACEXP > 0.15) to exclude data from the edge of the detectors. We created a combined spectrum from the data of TMs 1–4 & 6, the five cameras with an on-chip optical blocking filter, with an average exposure time of 294 s. TM5 and TM7 suffer from a light leak (Predehl et al. 2021) and no reliable energy calibration is available yet.

thumbnail Fig. 1.

Background-subtracted eROSITA light curve of eRASSt J040515.6 − 745202 during eRASS3 in the energy band 0.2–8.0 keV. Each data point represents one scan with a typical duration of 40 s and combines data from all seven cameras. The count rates were corrected for vignetting and point-spread function losses. The time in hours starts at the beginning of the first scan.

2.2. Swift

Following the eROSITA discovery of eRASSt J040515.6−745202 (Rau et al. 2021), we initiated follow-up target of opportunity observations with the Neil Gehrels Swift Observatory (Gehrels et al. 2004). We used the data obtained from the X-ray telescope (XRT) in photon counting mode with a time resolution of 2.5 s. The observations were performed on 2021 May 12 and a total XRT exposure time of 1682 s was obtained. The XRT data were analysed using standard procedures (Evans et al. 2007, 2009).

2.3. XMM-Newton

To investigate the nature of eRASSt J040515.6 − 745202, we triggered one of our XMM-Newton anticipated target of opportunity observations (PI FH) to follow up on new hard transient sources in the Magellanic system. For this purpose, we utilised data from the European Photon Imaging Camera (EPIC), sensitive to the 0.15–12 keV band. Two of the three XMM-Newton telescopes are equipped with Metal Oxide Semi-conductor (MOS) CCD arrays (Turner et al. 2001) and the third with a pn-CCD (a predecessor of the eROSITA CCDs; Strüder et al. 2001). The observation was performed on 2021 May 14 (observation ID 0860800401, start MJD = 59348.0497). We used the cameras with medium optical blocking filters, EPIC-pn in large-window readout mode (48 ms time resolution), and EPIC-MOS in full-frame mode (2.6 s). We obtained net exposure times of 25.7 ks and 27.9 ks for EPIC-pn and EPIC-MOS, respectively, after we removed intervals of high background flaring activity, which occurred only at the end of the observation.

XMM-Newton/EPIC data were processed using the XMM-Newton data analysis software SAS, version 19.1.02. We used the SAS task evselect for event extraction from circular regions around the source (45″ radius) and nearby background (60″ radius). Events with PATTERN 1–4 for EPIC-pn (Strüder et al. 2001) and PATTERN 1–12 for EPIC-MOS (Turner et al. 2001) were selected. We created light curves in the energy bands 0.2–2.2 keV, 2.2–8.0 keV, and 0.2–8.0 keV, applying the standard filtering with flags #XMMEA_EP and #XMMEA_EM for EPIC-pn and EPIC-MOS, respectively. For spectra we used the conservative FLAG = 0 filtering and created response files using the SAS tasks arfgen and rmfgen. As the best X-ray source position, we used the one derived by the XMM-Newton pipeline of and with a 1σ statistical uncertainty of and a remaining systematic error of after astrometric correction.

2.4. NICER

The NICER X-ray Timing Instrument (XTI, Gendreau et al. 2012, 2016) is a non-imaging, soft X-ray telescope on board the International Space Station (ISS). The XTI consists of an array of 56 co-aligned concentrator optics, each associated with a silicon drift detector (Prigozhin et al. 2012), operating in the 0.2–12 keV band. The XTI provides high time resolution (∼100 ns) and spectral resolution of ∼85 eV at 1 keV. It has a field of view of ∼30 arcmin2 in the sky and an effective area of ∼1900 cm2 at 1.5 keV (with 52 active detectors).

For the current study, we used NICER data obtained between MJD 59339 and 59346 (i.e. before the XMM-Newton observation) which were grouped into six observation intervals (observation IDs 4570010101–6). Data were analysed with HEASOFT version 6.29c, NICER DAS version 2020-04-23_V007a, and the calibration database (CALDB) version 20200722 (redistribution matrix and ancillary response files). For event selection and background filtering, we used standard screening criteria and intervals with magnetic cut-off rigidity between 1–15 GeV/c. Since NICER is not an imaging instrument, the X-ray background must be calculated from a database of data collected from empty background fields. The ‘space weather’ method (Gendreau et al., in prep.) contains functions which can be used to create the background spectrum based on estimates of the space weather environment. On the other hand, the 3C50 tool (Remillard et al. 2022) uses a number of background proxies from each NICER observation to define the basis states of the background database. We opted to use the latter since, for observations in the direction of the Magellanic Clouds, it has been shown to produce better estimates of the background (Treiber et al. 2021). After filtering the data, we were left with an exposure of ∼12.8 ks.

3. X-ray data analysis

3.1. Temporal analysis

The background-subtracted 0.2–8.0 keV light curves of eRASSt J040515.6 − 745202 obtained from the three EPIC instruments are shown in the upper part of Fig. 2. A short flare-like feature is consistently seen in all three light curves, with a fast rise and exponential decay as it is demonstrated in the zoom-in in the lower part of Fig. 2. A hardness ratio (HR) derived from the light curves in the soft (0.2–2.2 keV) and hard (2.2–8.0 keV) energy bands exhibits a similar rise and decay as the flux of the flare (Fig. 3). An investigation of the burst light curve (0.2–8 keV, total EPIC count rate) at higher time resolutions down to 1 s allowed us to estimate the rise time to ∼7–10 s. The flux rise appears linear without significant dips, indicating the absence of photospheric radius expansion. Also, no significant HR variations are seen around the burst peak at 10 s time resolution and the statistics does not allow for this to be resolved better.

thumbnail Fig. 2.

Top: background-subtracted EPIC-MOS1, MOS2, and pn light curves of eRASSt J040515.6 − 745202 in the energy band 0.2–8.0 keV with time bins of 50 s. Time zero corresponds to 2021 May 14 00:39:53 UT. The red dashed line marks the average count rate obtained from the full EPIC-pn light curve. Bottom: zoom-in showing the X-ray burst with a time binning of 4 s. The light curves for MOS1 and MOS2 were summed up to gain statistics. Time zero corresponds to the burst onset (start of first bin with a significantly increased count rate).

thumbnail Fig. 3.

Background-subtracted EPIC light curves (here MOS1, MOS2, and pn count rates are summed up) of eRASSt J040515.6 − 745202 in two energy bands with time bins of 10 s, together with the hardness ratio (count rate in the hard band divided by the count rate in the soft band) in the bottom panel. Time zero corresponds to the burst onset, as in Fig. 2

The temporal and spectral evolution (seen through the HR) of the flare is consistent with that of a type-I X-ray burst seen from LMXBs. These bursts are mainly detected from LMXBs in the Galaxy with various observatories, which allows for their spectral evolution with high time resolution to be investigated (e.g. Keek et al. 2018; Degenaar et al. 2016; Chelovekov et al. 2005; Haberl et al. 1987; Swank et al. 1977).

Apart from the burst, the EPIC light curves show a dipping behaviour during the last ∼1 h of the observation. The dips typically last 150–200 s and reach minima, which are consistent with zero flux. The dips are more pronounced in the 0.2–2.2 keV band, which results in increased HR. This is indicative of absorption as the origin for the dips.

We created power density spectra (PDS) for both the burst and quiescent intervals from the above light curves which were divided into segments of 4096 s with a time resolution of 10 s. The PDS obtained from all these segments were rebinned with a geometric rebinning factor and averaged to produce the final PDS. The average PDS was also normalised such that its integral gives the squared rms fractional variability and the expected white noise level was subtracted. The PDS from the quiescent interval shows two putative quasi-periodic oscillations (QPOs) at ∼3.8 mHz and its harmonic at ∼7.6 mHz as seen in Fig. 4. We fitted the PDS with a power-law component corresponding to the continuum and fitted the QPO features with two Lorentzian functions. The detection significance of the fundamental frequency is 2.4σ, a quality factor (ν/FWHM) of ∼14, and an rms fractional variability of ∼3% in the 0.2–8 keV range. The phenomenon of millihertz oscillations has been reported from a small subset of bursters at low frequencies (≃7–10 mHz) with an rms variability of ∼2–3% and it can be associated with the transition from stable to unstable nuclear burning, or with the regime of marginally stable nuclear burning (see e.g. Revnivtsev et al. 2001; Heger et al. 2007). Millihertz QPOs are observed to be suppressed after an X-ray burst before they reappear later (Lyu et al. 2014). Hence, we created separate PDS from the persistent emission before and after the burst. Unfortunately, the statistics in the two parts is insufficient to detect the QPOs.

thumbnail Fig. 4.

PDS of eRASSt J040515.6 − 745202 produced from the EPIC-pn light curve in the energy range of 0.2–8 keV as a histogram with error bars. The solid line represents a model fit with a power law and two Lorentzian line profiles.

In Fig. 5 we show the NICER X-ray light curve that is used to further investigate the variability of the source. One can clearly recognise a few low-flux states that sporadically occur. The source seems to transition between a high and a low state as it can be seen by the histogram distribution of the count rate. The low state may be defined by the 0.5–1.2 cts s−1 intensity range.

thumbnail Fig. 5.

NICER X-ray light curve in the 0.2–8.0 keV energy band using 64 s bins. The distribution of fluxes (shown in the inset) is clearly bimodal, with six low-flux states identified as shown in the shaded region (i.e. 0.5–1.2 cts s−1).

We searched for a possible periodic nature of the dips on short (< 20 ks) and longer (20–200 ks) timescales. For short timescales, we used the NICER light curve with photons grouped into 10 s bins, and we applied a method based on the phase dispersion minimisation test described by Stellingwerf (1978), which has been effectively applied to light curves of other sources (e.g. Vasilopoulos et al. 2021). We found hints of a periodic nature at a period of ∼6 ks and its harmonic. However, the low states could not be completely disentangled from the high (normal) state as seen in Fig. 6.

thumbnail Fig. 6.

NICER light curve (10 s bins) folded for a period of 6 ks. The shaded region marks the low-flux states as in Fig. 5.

Power spectra of the EPIC light curves, which cover the full frequency range between ∼4 × 10−5 Hz (limited by the observation length) and ∼10 Hz (constrained by the EPIC-pn time resolution), do not show evidence for any significant signal near 6 ks (1.7 × 10−4 Hz) nor at the highest frequencies. For longer timescales, we found that the Lomb-Scargle periodogram (Lomb 1976; Scargle 1982) was more effective. We applied the method following VanderPlas (2018) to the NICER light curve. We experimented with various bin sizes between 10 s and 300 s, which yielded consistent results. The periodogram revealed peaks at ∼78.5 ks and its harmonic (Fig. 7, top). We folded the NICER light curve for periods of 78 ks and 79 ks as shown in the bottom of Fig. 7. The phases covered by the XMM-Newton observation are indicated at the top of the plots. For the case of a period of 78 ks, dips are expected at the end of the XMM-Newton observation as they are observed (Fig. 2). In particular, from the EPIC light curves, it is evident that the associated decrease in flux does not follow the pattern expected from total or partial eclipses, but it more so resembles a dipping behaviour.

thumbnail Fig. 7.

Top: Lomb-Scargle periodogram of NICER data, a peak is visible at a period of 79 ks. Bottom: NICER light curve (300 s bin) folded for a period of 78 ks (left) and 79 ks (right). With red arrows, we mark the duration of the XMM-Newton observation. The shaded region in phase outlines the low-flux states as shown in Fig. 5.

Finally, we investigated the NICER data for possible X-ray bursts using 5 s and 10 s time bins, but no bursts were identified. This is not unexpected since the combined duration of the NICER exposures is shorter than the XMM-Newton observation.

3.2. Spectral analysis

Analysis of the X-ray spectra was performed using XSPEC v12.10.1f (Arnaud 1996). For the spectral analysis of the burst, we extracted EPIC spectra from the quiescent emission of the full observation, excluding a ∼3100 s interval around the burst (source spectra binned to a minimum of 20 counts per bin to use Gaussian statistics), and spectra from the burst starting with burst rise and exposure times of 166 s (pn) and 174 s (MOS) (binned to a minimum of 1 count to use Cash statistics, Cash 1979). The number of extracted source counts during the burst is ∼1200 counts from all three instruments together, which did not allow us to monitor spectral evolution during the burst. To investigate the burst spectra, we first fitted the (background-subtracted) quiescent emission using a two-component model comprised of a power law and a disk blackbody, both attenuated by photo-electric absorption along the line of sight. The spectra with the best-fit model are presented in Fig. 8 with the corresponding model parameters listed in Table 1. For the fit to the burst spectra, the persistent emission components were included in the model with fixed parameters as derived above, that is to say we assume that the persistent emission did not change during the burst (however, see Worpel et al. 2013; Degenaar et al. 2016). For simplicity, we also adopted that the burst and persistent emission were attenuated by the same absorption. The derived NH value of ∼11.5 × 1020 cm−2 is somewhat larger than the Galactic foreground H I value of 8.4× 1020 cm−2 (Dickey & Lockman 1990) and consistent with 10.5× 1020 cm−2 as derived from the HI4PI map in this direction (HI4PI Collaboration et al. 2016) or 13.9× 1020 cm−2 when accounting for the contribution of H2 (Willingale et al. 2013). The latter values include the H I emission of the Galactic foreground and the Magellanic Bridge (see also Sasaki et al. 2022). If the absorption mainly originates in the interstellar medium, the column density determined from the X-ray spectra suggests a location of the source in the Magellanic Bridge.

thumbnail Fig. 8.

EPIC spectra of the persistent emission from eRASSt J040515.6 − 745202 (MOS1, black; MOS2, red; and pn, green) extracted from the full observation, excluding the burst interval. The best-fit model consisting of a disk-blackbody and a power-law component is plotted as histogram. The two emission components are indicated as dotted lines (for better visibility for the pn spectrum only). The bottom panel shows the residuals.

Table 1.

Results of spectral fits to the X-ray spectra of the persistent emission.

For the burst emission, we started with a single blackbody component, which formally yielded an acceptable fit (C-statistic 758 for 951 d.o.f.), but it showed systematic residuals. Therefore, we added a second blackbody component to account for the cooling during the burst decay, as indicated by the hardness ratio (Fig. 3). The burst spectra with the best-fit model (C-statistic 741 for 949 d.o.f.) are presented in Fig. 9. For the temperatures, values of kT1 = 1.55 keV and kT2 = 0.23 keV were obtained. These two emission components dominate the spectrum at high and low energies, respectively (see Fig. 9). The observed average burst flux in the 0.2–10 keV band was 1.65 × 10−11 erg cm−2 s−1, which corresponds to an absorption-corrected X-ray luminosity of 6.2 × 1036 erg s−1, assuming a distance of 55 kpc (between that of the SMC and LMC). The bolometric luminosities inferred from the two blackbody normalisations are 6.54 × 1036 erg s−1 and 0.33 × 1036 erg s−1 for the high- and low-temperature component, respectively, which yield a total average bolometric burst luminosity of 6.87 × 1036 erg s−1. The average EPIC count rate over the 170 s burst interval inferred from the light curve with a 4 s binning (Fig. 2) was 6.86 cts s−1, while the maximum reached 16.1±1.9 cts s−1. This implies a bolometric luminosity at a burst peak of at least 1.4 × 1037 erg s−1, but probably higher because the burst peak is not resolved with 4 s bins. The apparent emission radii for the two blackbody components, assuming spherical emission, are R1 = 3.0±0.5 km and R2 = 31 km. However, it should be noted that these can only be treated as a rough estimate because they were derived from a spectrum averaged over the burst which is expected to show large temperature variations during its time evolution.

thumbnail Fig. 9.

EPIC spectra of the X-ray burst seen from eRASSt J040515.6 − 745202 (MOS1, black; MOS2, red; and pn, green). The best-fit model consists of two blackbody emission components in addition to the persistent emission. The best-fit model is plotted as a histogram. The two blackbody burst emission components are indicated as green dotted lines and the two components of the persistent emission as blue dotted lines (for better visibility for the pn spectrum only). The bottom panel shows the residuals.

To approximate a multi-temperature spectrum, we used a disk-blackbody model with a power-law distribution of the temperature (diskpbb in XSPEC). Although this model is not physical for an application to a temperature decrease with time (which is likely more exponential), it probably allows us to better estimate the maximum temperature. In the two-blackbody model, the temperatures only provide averages over a certain time interval during the burst decay. The diskpbb model yields a similar fit quality (C-statistic of 747 for 950 d.o.f., with one free parameter less than for the two-blackbody model) with a maximum keV. Peak temperatures above 2 keV are typically observed from known bursters (70% of the bursts have a peak kT > 2 keV, Galloway et al. 2020).

We fitted the NICER spectra with a similar persistent-emission model to look for possible spectral evolution. Given the nature of the NICER background, we performed the fitting using PG statistics within XSPEC, which implements Cash statistics (Cash 1979), with a non-Poisson background model3. To maximise signal-to-noise, we only fitted NICER data of the high (normal) state (i.e. exposure of 8 ks). Data were initially fitted by an absorbed power law while a soft disk blackbody was added to account for residuals at low energies. Best-fit (PG statistics of 880 for 764 d.o.f.) parameter values are given in Table 1, while the best-fit model is plotted in Fig. 10.

thumbnail Fig. 10.

NICER spectrum of the high (normal) state (i.e. count rate > 1.5 cts s−1). As for Fig. 8, the contributions of disk blackbody (dash-dotted) and power law (dotted) are shown separately.

We then compared the intensity of the low states seen in the NICER light curve with the NICER background level. For that, we extracted source and background spectra using the respective time intervals. This resulted in a net count rate of 0.66 cts s−1 and a background of 0.57 cts s−1, using the 3C50 model. Consequently, the low states are marginally brighter than the background level, and the source almost fades away during these intervals.

We also fitted the eROSITA spectrum with the two-component model of a power law and disk blackbody. Because the statistical quality of the spectrum is too low to justify a fit with so many degrees of freedom, we fixed the NH at the average value obtained from the EPIC and NICER spectra. The spectrum with the best-fit model is presented in Fig. 11. For the model parameters, readers can refer to Table 1. The disk-blackbody temperature derived from the eROSITA spectrum is lower than that found from the EPIC and NICER spectra, which may be caused by the higher low-energy contribution of the steeper power law. Fixing either the power-law index or the temperature at the EPIC/NICER average values, however, does not bring the corresponding parameter in agreement, indicating a somewhat lower contribution of the disk component during the eROSITA observations. We note that the spectral model used by Rau et al. (2021) includes a blackbody component instead of a disk blackbody. This has the effect that in their spectral fit, the steeper power law accounts for the soft part of the spectrum, while the hot blackbody represents the harder emission. In the model of this work, the role of the two emission components is reversed, which leads to the different best-fit parameters.

thumbnail Fig. 11.

eROSITA spectrum obtained from eRASS3. The best-fit model consists of a disk-blackbody (dash-dotted) and a power-law (dotted) component as used for the persistent emission seen during the XMM-Newton observation.

The source was detected with Swift/XRT with a net count rate of (0.094±0.008) cts s−1. Given the small number of source counts (∼150), spectra were binned to a minimum of 1 count per bin, while spectral fitting was performed based on Cash statistics. As for the fit to the eROSITA spectrum, the NH was fixed to 11.5× 1020 cm−2. A simple absorbed power-law model was adequate to explain the data (C-Statistics of 86.8 for 120 degrees of freedom, with a corresponding χ2 of 106.6, Table 1). The power-law photon index had a value of ∼1.6, consistent with what was derived from the other X-ray spectra, while the absorption-corrected X-ray luminosity was found to be ∼2.3 × 1036 erg s−1 (0.2–10.0 keV), which is also in the range seen during the other observations.

4. UV, optical, and near-infrared data

As described by Rau et al. (2021), Swift observed the region around eRASSt J040515.6 − 745202 with the Ultra-Violet/Optical telescope (UVOT, Roming et al. 2005) and the UVW1 filter4 with an exposure time of 1754 s. We downloaded the Swift/UVOT data from the UK Swift Science Data Centre5, which includes the UVW1 image (shown in Fig. 12, left) and a source detection list produced by the data centre pipeline. Source magnitudes are given in the Vega system. The two UVOT sources listed in Rau et al. (2021) are not contained in this source list and are fainter than the faintest source (21.48 mag) in the list, indicating their lower significance.

thumbnail Fig. 12.

UVW1-band images from Swift/UVOT (left) and XMM-Newton/OM. The large white circle (radius 5″) marks the eROSITA error circle from Rau et al. (2021), and the three concentric white circles indicate the 1, 2, and 3σ position uncertainties obtained from the XMM-Newton/EPIC X-ray images. The magenta circles mark the positions of the UVOT sources given by Rau et al. (2021), and the green circles (1σ errors) indicate the sources detected in the XMM-Newton/OM UVW1 image.

The source region was also observed with the XMM-Newton optical/UV monitor telescope (OM; Mason et al. 2001), using the UVW1 filter6 and the detector in fast imaging mode. The total exposure time was 20.98 ks, which is significantly higher than the UVOT exposure. We downloaded the pipeline products generated from the OM data from the XMM-Newton science archive7. The OM/UVW1 image is shown in Fig. 12 (right). UVOT source 1 is also detected in the OM image (source ID 44 in the OM detection list). Source 1 is ruled out in any case as the UV counterpart of eRASSt J040515.6 − 745202 by the improved XMM-Newton position. A clear source (OM 42) is detected inside the 2σXMM-Newton/EPIC error radius with a UVW1 magnitude of 21.68±0.08. It is located at RA = 04h05m14928 and Dec. = −74°52′0046 (J2000) with a 1σ error of 0.19″.

Comparing OM and UVOT images (Fig. 12), UVOT source 2 might suffer from a cross-like feature, which is visible in the UVOT image. This could be caused by statistical background fluctuations due to the short exposure time or an independent transient source (at the position of UVOT 2) blended with OM 42. We regard OM source 42 as more reliable as its shape is mostly circular, which is consistent with the point spread function of the instrument. We found no counterpart at any wavelength within 5″ to OM source 42 using the VizieR catalogue access tool8.

An optical counterpart was identified in archival data from DECam (Flaugher et al. 2015), a wide-field CCD camera on the Blanco 4-m telescope at the Cerro Tololo Inter-American Observatory. Twenty-four exposures in g, r, i, and z filters from three programmes were used (including the DECam eROSITA Survey, DeROSITAS, an optical companion survey to the eROSITA survey). The exposures were instrumentally calibrated by the DECam Community Pipeline (Valdes 2021) and catalogued for the Legacy Surveys’ (Dey et al. 2019) data release (DR10 in preparation). The analysis uses the Tractor (Lang et al. 2016) forward modelling method. The method fits models to the observations, taking the point spread functions from the exposures into account. The Tractor fitting agrees with both a very compact galaxy or a point source; that is to say, the data are consistent with either a slightly extended source or a star. The measured AB magnitudes, an absolute spectral flux density photometric system, are [24.08 ± 0.11, 23.21 ± 0.07, 23.05 ± 0.05, 22.27 ± 0.09] in [g, r, i, z], respectively. The photometric calibration is tied to Pan-STARS DR1 (Chambers et al. 2016) through an uber-cal self-calibration method (Schlafly et al. 2012). The images centred on the position obtained for the UV counterpart of eRASSt J040515.6 − 745202 are shown in Fig. 13.

thumbnail Fig. 13.

Cutouts of the DECam imaging at the same size and position as Fig. 12. Each image is a coadd of many exposures in each filter and also rendered as a colour composite. The positions of OM source 42 and UVOT source 2 are overlaid in a copy of the composite for reference. The exposures used are public archival data included in the Legacy Surveys data release 10 (DR10).

5. Discussion

During our eROSITA project to identify hard X-ray transients in the Magellanic system, an alert was generated by the eROSITA Near-Real Time Analysis system (Predehl et al. 2021). The transient, eRASSt J040515.6 − 745202, was only detected during eRASS3 and is located in the direction of the Magellanic Bridge (Rau et al. 2021).

A follow-up observation with XMM-Newton revealed an X-ray burst with energetics that suggest the source is located at the distance of the Magellanic Bridge. Assuming a distance of 55 kpc (between that of the SMC and LMC), a burst peak bolometric luminosity of more than 1.4 × 1037 erg s−1 is derived, which is typical of most X-ray bursts. Given the large distance, the statistics is low; however, we can infer certain characteristics of the burst.

According to theory, type-I X-ray bursts, which are caused by unstable burning of the accreted matter, occur at certain mass accretion rates. Three different regimes are predicted which produce different types of X-ray bursts (Bildsten et al. 2000; Fujimoto et al. 1981).

  1. Mixed H/He burning triggered by thermally unstable H ignition at low accretion rates of 10−14 M yr−1 ≲ 2 × 10−10 M yr−1 resulting in a type I X-ray burst in a H-rich environment.

  2. Pure He shell ignition after steady H burning at intermediate accretion rates of 2 × 10−10 M yr−1 ≲ 4–11 × 10−10 M yr−1.

  3. Mixed H/He burning triggered by thermally unstable He ignition at high accretion rates of 4–11 × 10−10 M yr−1 ≲ 2 × 10−8 M yr−1.

Moreover, sedimentation of heavy elements was suggested to play an important role for the burst environment. This can lead to additional burst regimes (Peng et al. 2007).

The critical accretion rates depend on metallicity, and in particular on the mass fraction of CNO (ZCNO ∼ 0.01 for solar abundances). Between regime 1 and 2 the transition is modified by a factor of (ZCNO/0.01)1/2. For Magellanic Cloud metallicities between 0.2 (SMC) and 0.5 solar (LMC), the critical accretion rates are, therefore, lowered by about a factor of 0.45–0.71.

To convert the persistent luminosity of eRASSt J040515.6 − 745202 (2.6 × 1036 erg s−1) to a mass accretion rate, we assumed standard neutron star parameters (1.4 M, a radius of 12 km, and a 100% efficiency for the conversion of potential energy into X-rays). This results in 2.6 × 10−10 M yr−1, which is somewhat above the critical rate between regime 1 and 2 where a He burst is expected for sub-solar metallicity.

Following in’t Zand et al. (2017), we fitted an exponential decay to the burst light curve (total EPIC count rate in the 0.2–8.0 keV band, Fig. 2). We find an exponential decay time of 70.6 s (1σ uncertainties). From the analysis of more than 1200 X-ray bursts, in’t Zand et al. (2017) find that the most common decay time is 5 s and only 1% have decay times longer than 70 s. For bursts driven by He burning, short (∼10 s) decay times are expected (Parikh et al. 2013), suggesting that the observed burst from eRASSt J040515.6 − 745202 ignited in a H-rich environment. This is consistent with the recurrence time, which must be longer than 5.06 h given by the constraint of the observation start. For He bursts, shorter recurrence times may be expected. For the α value, the ratio between integrated fluxes of persistent emission between the burst and the previous one, and the burst emission, we can infer a lower limit of 38. In the simplest model, assuming isotropic emission and all accreted matter is processed during the bursts; α values of ∼30 and ∼120 are expected for pure H and He burning, respectively. However, the picture is more complicated and α values between 10 and 1000 were observed (Sect. 3.7 of Lewin et al. 1993).

eRASSt J040515.6 − 745202 was scanned by eROSITA for about two days each during four epochs separated by 0.5 year intervals. The only eROSITA detection on May 1–2, 2021 and during follow-up observations with NICER (May 5–12), Swift (May 12), and XMM-Newton (May 14) suggest that the source was in a moderate outburst, lasting at least two weeks. During this outburst, the accretion rate was sufficiently high to gather enough hydrogen for the ignition of an X-ray burst. On longer timescales, the source accretes at a much lower rate, which makes the detection of bursts in a typical X-ray follow-up observation very unlikely.

During the NICER observations, eRASSt J040515.6−745202 alternated between two flux states. The high-flux level was similar to the persistent emission seen during the XMM-Newton observation, while during low states the flux was typically a factor of 4–5 lower. These dips in the NICER light curve could be due to absorption or even eclipses similar to those found in other LMXBs observed with NICER (e.g. Homan et al. 2021; Buisson et al. 2021). The source also exhibited dipping behaviour for the last ∼1 h of the XMM-Newton observation. The dips are resolved in the EPIC light curves, and they typically lasted for 150–200 s and reached minimum count rates consistent with zero. During the dips, the hardness ratio increased, which suggests they are the result of absorption and obscuration by optically thick material. As was proposed for other dipping LMXBs, this material is likely located in a thickened region of the accretion disk where the accretion flow from the secondary star impacts on the disk (Díaz Trigo et al. 2006). In systems with sufficiently high inclination, the neutron star then gets partially occulted by this region. From the NICER data, we find some evidence for a periodic nature of the dips with a period of around 78.5 ks (21.8 h). Since the dips are not resolved in the NICER light curves, the period remains uncertain, but the dip phases derived from the XMM-Newton observation are consistent with the NICER-derived period of 78–79 ks. With an orbital period of ∼21.8 h, eRASSt J040515.6 − 745202 would be very similar to Swift J1858.6−0814 (orbital period 21.34 h), which also exhibits X-ray bursts and dips and, in addition, eclipses with a 1.14 h duration (Buisson et al. 2021; Knight et al. 2022). Eclipses could also be present in the NICER light curves of eRASSt J040515.6 − 745202, but the low count rate and the duration of NICER snapshots combined with the low number of visits makes it difficult to detect the ingress and egress of eclipses and differentiate them from dips. Considering the similarities with Swift J1858.6−0814, the dipping behaviour could also be due to disk outflows (Castro Segura et al. 2022).

The detection of a type-I X-ray burst uniquely characterises eRASSt J040515.6 − 745202 as an LMXB with the compact object identified as a neutron star. From UV, optical, and IR images, we identified the counterpart to the X-ray source with a spectral energy distribution suggesting two emission components. The flux measurements in the g, r, i, and z bands increase with wavelength, consistent with radiation from a late-type M star (effective temperature of around 3000 K). The contribution of such a star in the UV is negligible and the bright UV emission detected by the XMM-Newton/OM most likely originates from the accretion disk around the neutron star. An early-type OB star to account for the UV emission can definitely be ruled out as it would be expected to be brighter than ∼16th magnitude in the optical B and V bands. This picture is fully consistent with the LMXB nature of the system and confirms the identification of the multi-wavelength counterpart to the X-ray source.

The only other LMXB known in the Magellanic System is LMC X-2, a persistently bright system in the southern part of the LMC (e.g. Agrawal & Nandi 2020; Smale & Kuulkers 2000; Bonnet-Bidaud et al. 1989). The location of LMC X-2 is outside the area where younger stellar populations are found, which are traced in X-rays by high-mass X-ray binaries (Haberl & Sturm 2016) and core-collapse supernova remnants (Maggi et al. 2016). eRASSt J040515.6 − 745202 is located even further out in the Magellanic Bridge, about one-third of the angular distance to the SMC. The discovery of an LMXB in the Magellanic Bridge adds support to the existence of an old population (Bagheri et al. 2013) and hence the theory of tidal interaction.


1

This was probably not recognised because a wrong observation ID was given by Pietsch & Haberl (2005) in the caption of their Fig. 4.

2

Science Analysis Software (SAS): https://www.cosmos.esa.int/web/xmm-newton/sas

4

Bandwidth (full width at half maximum) of 69.3 nm with a central wavelength of 260 nm.

6

Band width of 83 nm at an effective wavelength of 291 nm.

Acknowledgments

This work is based on data from eROSITA, the soft X-ray instrument aboard SRG, a joint Russian-German science mission supported by the Russian Space Agency (Roskosmos), in the interests of the Russian Academy of Sciences represented by its Space Research Institute (IKI), and the Deutsches Zentrum für Luft- und Raumfahrt (DLR). The SRG spacecraft was built by Lavochkin Association (NPOL) and its subcontractors, and is operated by NPOL with support from the Max Planck Institute for Extraterrestrial Physics (MPE). The development and construction of the eROSITA X-ray instrument was led by MPE, with contributions from the Dr. Karl Remeis Observatory Bamberg & ECAP (FAU Erlangen-Nürnberg), the University of Hamburg Observatory, the Leibniz Institute for Astrophysics Potsdam (AIP), and the Institute for Astronomy and Astrophysics of the University of Tübingen, with the support of DLR and the Max Planck Society. The Argelander Institute for Astronomy of the University of Bonn and the Ludwig Maximilians Universität Munich also participated in the science preparation for eROSITA. The eROSITA data shown here were processed using the eSASS software system developed by the German eROSITA consortium. This work used observations obtained with XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA Member States and NASA. The XMM-Newton project is supported by the DLR and the Max Planck Society. Work on NICER science at NRL is funded by NASA and the Office of Naval Research. This research has made use of the VizieR catalogue access tool, CDS, Strasbourg, France (DOI: http://dx.doi.org/10.26093/cds/vizier). The original description of the VizieR service was published in Ochsenbein et al. (2000). We thank the anonymous referee for the very useful suggestions, which helped to improve the paper.

References

  1. Agrawal, V. K., & Nandi, A. 2020, MNRAS, 497, 3726 [NASA ADS] [CrossRef] [Google Scholar]
  2. Arnaud, K. A. 1996, ASP Conf. Ser., 101, 17 [Google Scholar]
  3. Bagheri, G., Cioni, M. R. L., & Napiwotzki, R. 2013, A&A, 551, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Bildsten, L. 2000, AIP Conf. Ser., 522, 359 [NASA ADS] [CrossRef] [Google Scholar]
  5. Bonnet-Bidaud, J. M., Motch, C., Beuermann, K., et al. 1989, A&A, 213, 97 [NASA ADS] [Google Scholar]
  6. Brunner, H., Liu, T., Lamer, G., et al. 2022, A&A, 661, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Buisson, D. J. K., Altamirano, D., Armas Padilla, M., et al. 2021, MNRAS, 503, 5600 [NASA ADS] [CrossRef] [Google Scholar]
  8. Cash, W. 1979, ApJ, 228, 939 [Google Scholar]
  9. Castro Segura, N., Knigge, C., Long, K. S., et al. 2022, Nature, 603, 52 [NASA ADS] [CrossRef] [Google Scholar]
  10. Chambers, K. C., Magnier, E. A., Metcalfe, N., et al. 2016, arXiv e-prints [arXiv:1612.05560] [Google Scholar]
  11. Chelovekov, I. V., Lutovinov, A. A., Grebenev, S. A., & Sunyaev, R. A. 2005, Astron. Lett., 31, 681 [NASA ADS] [CrossRef] [Google Scholar]
  12. Degenaar, N., Koljonen, K. I. I., Chakrabarty, D., et al. 2016, MNRAS, 456, 4256 [NASA ADS] [CrossRef] [Google Scholar]
  13. Dey, A., Schlegel, D. J., Lang, D., et al. 2019, AJ, 157, 168 [Google Scholar]
  14. Díaz Trigo, M., Parmar, A. N., Boirin, L., Méndez, M., & Kaastra, J. S. 2006, A&A, 445, 179 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Dickey, J. M., & Lockman, F. J. 1990, ARA&A, 28, 215 [Google Scholar]
  16. Evans, P. A., Beardmore, A. P., Page, K. L., et al. 2007, A&A, 469, 379 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Evans, P. A., Beardmore, A. P., Page, K. L., et al. 2009, MNRAS, 397, 1177 [Google Scholar]
  18. Flaugher, B., Diehl, H. T., Honscheid, K., et al. 2015, AJ, 150, 150 [Google Scholar]
  19. Fujimoto, M. Y., Hanawa, T., & Miyaji, S. 1981, ApJ, 247, 267 [NASA ADS] [CrossRef] [Google Scholar]
  20. Galloway, D. K., & Keek, L. 2021, Astrophys. Space Sci. Lib., 461, 209 [NASA ADS] [CrossRef] [Google Scholar]
  21. Galloway, D. K., in’t Zand, J., Chenevez, J., et al. 2020, ApJS, 249, 32 [NASA ADS] [CrossRef] [Google Scholar]
  22. Gehrels, N., Chincarini, G., Giommi, P., et al. 2004, ApJ, 611, 1005 [Google Scholar]
  23. Gendreau, K. C., Arzoumanian, Z., & Okajima, T. 2012, SPIE Conf. Ser., 8443, 844313 [Google Scholar]
  24. Gendreau, K. C., Arzoumanian, Z., Adkins, P. W., et al. 2016, SPIE Conf. Ser., 9905, 99051H [NASA ADS] [Google Scholar]
  25. Haberl, F., & Sturm, R. 2016, A&A, 586, A81 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Haberl, F., Stella, L., White, N. E., Priedhorsky, W. C., & Gottwald, M. 1987, ApJ, 314, 266 [NASA ADS] [CrossRef] [Google Scholar]
  27. Haberl, F., Maitra, C., Carpano, S., et al. 2022, A&A, 661, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Heger, A., Cumming, A., & Woosley, S. E. 2007, ApJ, 665, 1311 [NASA ADS] [CrossRef] [Google Scholar]
  29. HI4PI Collaboration (Ben Bekhti, N., et al.) 2016, A&A, 594, A116 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Homan, J., Gendreau, K. C., Sanna, A., et al. 2021, ATel, 14606, 1 [NASA ADS] [Google Scholar]
  31. in’t Zand, J. J. M., Visser, M. E. B., Galloway, D. K., et al. 2017, A&A, 606, A130 [CrossRef] [EDP Sciences] [Google Scholar]
  32. Irwin, M. J., Demers, S., & Kunkel, W. E. 1990, AJ, 99, 191 [NASA ADS] [CrossRef] [Google Scholar]
  33. Keek, L., Arzoumanian, Z., Chakrabarty, D., et al. 2018, ApJ, 856, L37 [NASA ADS] [CrossRef] [Google Scholar]
  34. Knight, A. H., Ingram, A., & Middleton, M. 2022, MNRAS, 514, 1908 [NASA ADS] [CrossRef] [Google Scholar]
  35. Lang, D., Hogg, D. W., & Mykytyn, D. 2016, Astrophysics Source Code Library [record ascl:1604.008] [Google Scholar]
  36. Lewin, W. H. G., van Paradijs, J., & Taam, R. E. 1993, Space Sci. Rev., 62, 223 [CrossRef] [Google Scholar]
  37. Lomb, N. R. 1976, Ap&SS, 39, 447 [Google Scholar]
  38. Lyu, M., Méndez, M., & Altamirano, D. 2014, MNRAS, 445, 3659 [NASA ADS] [CrossRef] [Google Scholar]
  39. Maggi, P., Haberl, F., Kavanagh, P. J., et al. 2016, A&A, 585, A162 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Maitra, C., Haberl, F., Vasilopoulos, G., et al. 2021, A&A, 647, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Mason, K. O., Breeveld, A., Much, R., et al. 2001, A&A, 365, L36 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Mathewson, D. S. 1985, PASA, 6, 104 [NASA ADS] [Google Scholar]
  43. Nidever, D. L., Majewski, S. R., Butler Burton, W., & Nigra, L. 2010, ApJ, 723, 1618 [NASA ADS] [CrossRef] [Google Scholar]
  44. Ochsenbein, F., Bauer, P., & Marcout, J. 2000, A&AS, 143, 23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Parikh, A., José, J., Sala, G., & Iliadis, C. 2013, Prog. Part. Nucl. Phys., 69, 225 [NASA ADS] [CrossRef] [Google Scholar]
  46. Pastor-Marazuela, I., Webb, N. A., Wojtowicz, D. T., & van Leeuwen, J. 2020, A&A, 640, A124 [EDP Sciences] [Google Scholar]
  47. Peng, F., Brown, E. F., & Truran, J. W. 2007, ApJ, 654, 1022 [NASA ADS] [CrossRef] [Google Scholar]
  48. Pietsch, W., & Haberl, F. 2005, A&A, 430, L45 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Predehl, P., Andritschke, R., Arefiev, V., et al. 2021, A&A, 647, A1 [EDP Sciences] [Google Scholar]
  50. Prigozhin, G., Gendreau, K., Foster, R., et al. 2012, SPIE Conf. Ser., 8453, 845318 [NASA ADS] [Google Scholar]
  51. Rau, A., Grotova, I., Gokus, A., et al. 2021, ATel, 14646, 1 [NASA ADS] [Google Scholar]
  52. Remillard, R. A., Loewenstein, M., Steiner, J. F., et al. 2022, AJ, 163, 130 [NASA ADS] [CrossRef] [Google Scholar]
  53. Revnivtsev, M., Churazov, E., Gilfanov, M., & Sunyaev, R. 2001, A&A, 372, 138 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Roming, P. W. A., Kennedy, T. E., Mason, K. O., et al. 2005, Space Sci. Rev., 120, 95 [Google Scholar]
  55. Sasaki, M., Knies, J., Haberl, F., et al. 2022, A&A, 661, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Scargle, J. D. 1982, ApJ, 263, 835 [Google Scholar]
  57. Schlafly, E. F., Finkbeiner, D. P., Jurić, M., et al. 2012, ApJ, 756, 158 [NASA ADS] [CrossRef] [Google Scholar]
  58. Smale, A. P., & Kuulkers, E. 2000, ApJ, 528, 702 [NASA ADS] [CrossRef] [Google Scholar]
  59. Stellingwerf, R. F. 1978, ApJ, 224, 953 [Google Scholar]
  60. Strohmayer, T., & Bildsten, L. 2006, in Compact Stellar X-ray Sources, eds. W. Lewin, & M. van der Klis (Cambridge, UK: Cambridge Univ. Press), 39, 113 [NASA ADS] [CrossRef] [Google Scholar]
  61. Strüder, L., Briel, U., Dennerl, K., et al. 2001, A&A, 365, L18 [Google Scholar]
  62. Sunyaev, R., Arefiev, V., Babyshkin, V., et al. 2021, A&A, 656, A132 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Swank, J. H., Becker, R. H., Boldt, E. A., et al. 1977, ApJ, 212, L73 [NASA ADS] [CrossRef] [Google Scholar]
  64. Treiber, H., Vasilopoulos, G., Bailyn, C. D., et al. 2021, MNRAS, 503, 6187 [CrossRef] [Google Scholar]
  65. Turner, M. J. L., Abbey, A., Arnaud, M., et al. 2001, A&A, 365, L27 [CrossRef] [EDP Sciences] [Google Scholar]
  66. Valdes, F. 2021, The DECam Community Pipeline Instrumental Signature Removal [https://legacy.noirlab.edu/noao/staff/fvaldes/Pipelines/Docs/PL206/] [Google Scholar]
  67. VanderPlas, J. T. 2018, ApJS, 236, 16 [Google Scholar]
  68. Vasilopoulos, G., Koliopanos, F., Haberl, F., et al. 2021, ApJ, 909, 50 [Google Scholar]
  69. Verner, D. A., Ferland, G. J., Korista, K. T., & Yakovlev, D. G. 1996, ApJ, 465, 487 [Google Scholar]
  70. Willingale, R., Starling, R. L. C., Beardmore, A. P., Tanvir, N. R., & O’Brien, P. T. 2013, MNRAS, 431, 394 [Google Scholar]
  71. Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914 [Google Scholar]
  72. Worpel, H., Galloway, D. K., & Price, D. J. 2013, ApJ, 772, 94 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1.

Results of spectral fits to the X-ray spectra of the persistent emission.

All Figures

thumbnail Fig. 1.

Background-subtracted eROSITA light curve of eRASSt J040515.6 − 745202 during eRASS3 in the energy band 0.2–8.0 keV. Each data point represents one scan with a typical duration of 40 s and combines data from all seven cameras. The count rates were corrected for vignetting and point-spread function losses. The time in hours starts at the beginning of the first scan.

In the text
thumbnail Fig. 2.

Top: background-subtracted EPIC-MOS1, MOS2, and pn light curves of eRASSt J040515.6 − 745202 in the energy band 0.2–8.0 keV with time bins of 50 s. Time zero corresponds to 2021 May 14 00:39:53 UT. The red dashed line marks the average count rate obtained from the full EPIC-pn light curve. Bottom: zoom-in showing the X-ray burst with a time binning of 4 s. The light curves for MOS1 and MOS2 were summed up to gain statistics. Time zero corresponds to the burst onset (start of first bin with a significantly increased count rate).

In the text
thumbnail Fig. 3.

Background-subtracted EPIC light curves (here MOS1, MOS2, and pn count rates are summed up) of eRASSt J040515.6 − 745202 in two energy bands with time bins of 10 s, together with the hardness ratio (count rate in the hard band divided by the count rate in the soft band) in the bottom panel. Time zero corresponds to the burst onset, as in Fig. 2

In the text
thumbnail Fig. 4.

PDS of eRASSt J040515.6 − 745202 produced from the EPIC-pn light curve in the energy range of 0.2–8 keV as a histogram with error bars. The solid line represents a model fit with a power law and two Lorentzian line profiles.

In the text
thumbnail Fig. 5.

NICER X-ray light curve in the 0.2–8.0 keV energy band using 64 s bins. The distribution of fluxes (shown in the inset) is clearly bimodal, with six low-flux states identified as shown in the shaded region (i.e. 0.5–1.2 cts s−1).

In the text
thumbnail Fig. 6.

NICER light curve (10 s bins) folded for a period of 6 ks. The shaded region marks the low-flux states as in Fig. 5.

In the text
thumbnail Fig. 7.

Top: Lomb-Scargle periodogram of NICER data, a peak is visible at a period of 79 ks. Bottom: NICER light curve (300 s bin) folded for a period of 78 ks (left) and 79 ks (right). With red arrows, we mark the duration of the XMM-Newton observation. The shaded region in phase outlines the low-flux states as shown in Fig. 5.

In the text
thumbnail Fig. 8.

EPIC spectra of the persistent emission from eRASSt J040515.6 − 745202 (MOS1, black; MOS2, red; and pn, green) extracted from the full observation, excluding the burst interval. The best-fit model consisting of a disk-blackbody and a power-law component is plotted as histogram. The two emission components are indicated as dotted lines (for better visibility for the pn spectrum only). The bottom panel shows the residuals.

In the text
thumbnail Fig. 9.

EPIC spectra of the X-ray burst seen from eRASSt J040515.6 − 745202 (MOS1, black; MOS2, red; and pn, green). The best-fit model consists of two blackbody emission components in addition to the persistent emission. The best-fit model is plotted as a histogram. The two blackbody burst emission components are indicated as green dotted lines and the two components of the persistent emission as blue dotted lines (for better visibility for the pn spectrum only). The bottom panel shows the residuals.

In the text
thumbnail Fig. 10.

NICER spectrum of the high (normal) state (i.e. count rate > 1.5 cts s−1). As for Fig. 8, the contributions of disk blackbody (dash-dotted) and power law (dotted) are shown separately.

In the text
thumbnail Fig. 11.

eROSITA spectrum obtained from eRASS3. The best-fit model consists of a disk-blackbody (dash-dotted) and a power-law (dotted) component as used for the persistent emission seen during the XMM-Newton observation.

In the text
thumbnail Fig. 12.

UVW1-band images from Swift/UVOT (left) and XMM-Newton/OM. The large white circle (radius 5″) marks the eROSITA error circle from Rau et al. (2021), and the three concentric white circles indicate the 1, 2, and 3σ position uncertainties obtained from the XMM-Newton/EPIC X-ray images. The magenta circles mark the positions of the UVOT sources given by Rau et al. (2021), and the green circles (1σ errors) indicate the sources detected in the XMM-Newton/OM UVW1 image.

In the text
thumbnail Fig. 13.

Cutouts of the DECam imaging at the same size and position as Fig. 12. Each image is a coadd of many exposures in each filter and also rendered as a colour composite. The positions of OM source 42 and UVOT source 2 are overlaid in a copy of the composite for reference. The exposures used are public archival data included in the Legacy Surveys data release 10 (DR10).

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.