Discovery of a Rare Eclipsing Be/X-ray Binary System, Swift J010902.6-723710 = SXP 182

Abstract

We report on the discovery of Swift J010902.6-723710, a rare eclipsing Be/X-ray Binary system by the Swift SMC Survey (S-CUBED). Swift J010902.6-723710 was discovered via weekly S-CUBED monitoring observations when it was observed to enter a state of X-ray outburst on 10 October 2023. X-ray emission was found to be modulated by a 182s period. Optical spectroscopy is used to confirm the presence of a highly-inclined circumstellar disk surrounding a B0-0.5Ve optical companion. Historical UV and IR photometry are then used to identify strong eclipse-like features re-occurring in both light curves with a 60.623 day period, which is adopted as the orbital period of the system. Eclipsing behavior is found to be the result of a large accretion disk surrounding the neutron star. Eclipses are produced when the disk passes in front of the OBe companion, blocking light from both the stellar surface and circumstellar disk. This is only the third Be/X-ray Binary to have confirmed eclipses. We note that this rare behavior provides an important opportunity to constrain the physical parameters of a Be/X-ray Binary with greater accuracy than is possible in non-eclipsing systems.

1 Introduction

Be/X-ray Binaries (BeXRBs) are a type of interacting High Mass X-ray Binary (HMXB) that contain a main sequence OBe star and a compact object, typically a neutron star (NS). These binaries are characterized by moderately eccentric elliptical orbits ($e\sim 0.3$), orbital periods on the order of $\sim$10-1000 days, and the presence of strong emission lines such as the Balmer series H$\alpha$ line in the optical spectrum of the OBe star (for a comprehensive review of BeXRBs, see Reig (2011)). The optical emission lines are interpreted to be a strong indication of a geometrically thin circumstellar disk (Rivinius, 2019) that surrounds the donor star and is variable in size. Interactions between the NS and disk are capable of producing intermittent X-ray outbursts, which are perhaps the most prominent feature of BeXRB systems.

There are two types of outburst that can be produced by NS interactions with the circumstellar disk (Stella et al., 1986; Okazaki & Negueruela, 2001; Reig, 2011). Type I outbursts are periodic in nature and associated with NS-disk interactions that occur at the periastron passage in each orbit(Stella et al., 1986; Okazaki & Negueruela, 2001). Type II outbursts are much more luminous than Type I outbursts (Okazaki & Negueruela, 2001; Reig, 2011) and can even reach super-Eddington luminosities, such as during the SMC X-3 outburst described by Townsend et al. (2017) and the SMC X-2 outburst that occurred in 2015 (Jaisawal et al., 2023; Roy et al., 2022). These luminous outbursts are independent of the orbital phase of the NS at the time of outburst and can last for multiple orbits (Reig, 2011). Due to both the long duration and phase independence of these events, Type II outbursts are thought to be related to the growth and shape of the Be star’s disk (Martin et al., 2014; Monageng et al., 2017) causing it to interact with the NS for a longer time.

Since X-ray outbursts in BeXRBs are dependent on the disk of their donor star, these systems are prone to experiencing long quiescent states during which they are hard to detect and identify (Coe et al., 2023b). Swift J010902.6-723710 is an example of a system that escaped identification due to a long period of quiescence. This newly-discovered system resides in the Small Magellanic Cloud (SMC), a satellite galaxy of the Milky Way that is well-documented to have a large population of HMXBs (Antoniou et al., 2010) due to an unusually-active period of star formation that occurred in the SMC’s recent past (Harris & Zaritsky, 2004; Rezaeikh et al., 2014). New SMC BeXRBs are still occasionally discovered (Maitra et al., 2023), indicating that a hidden population of previously-undetected systems still exists within the dwarf galaxy.

In an effort identify X-ray outbursts and detect new sources, the Swift SMC Survey (S-CUBED) (Kennea et al., 2018) has operated since 2016 to provide regular X-ray and UV monitoring of the SMC. This survey utilizes the Neil Gehrels Swift Observatory (Gehrels et al., 2004) to tile the galaxy on a weekly cadence with 142 overlapping tiles and $\sim$1 minute spent on each tile. Data recorded by S-CUBED are automatically analyzed to flag new X-ray outburst events that are detected by Swift’s X-ray Telescope (XRT; Burrows et al. 2005). S-CUBED also observes all tiles using the Ultraviolet/Optical Telescope (UVOT; Roming et al. 2005) in the uvw1 band centered at 2600Å. The archive of X-ray and UV photometric data taken by S-CUBED has proved to be an invaluable resource for discovering previously-undetected BeXRBs (Monageng et al., 2020; Kennea et al., 2020; Coe et al., 2021; Kennea et al., 2021; Monageng et al., 2022). It was through an X-ray selected search of the S-CUBED archive that Swift J010902.6-723710 was first identified as a candidate source in the summer of 2023, several months before it entered outburst during October 2023.

In this paper, we report on the recent X-ray outburst that was experienced by Swift J010902.6-723710 which confirms its status as a newly-identified BeXRB system (Coe et al., 2023a). In the following sections, we present the results of observations made both before and during outburst using Swift, the Optical Gravitational Lensing Experiment (OGLE), and the Southern African Large Telescope (SALT). Additionally, all results are combined and analyzed in an effort to obtain information about the orbital and physical properties of both the NS and its mass donor OBe companion.

2 Observations

2.1 Swift - XRT

Figure 1: XRT light curve for Swift J010902.6-723710 from August 2023 to present, combining data from S-CUBED observations and follow-up TOOs. Note source visibility to Swift is indicated by the grey background, white vertical strips indicate periods when the target could not be observed. Vertical dotted lines mark predicted times for the periastron passages.

Swift J010902.6-723710 has remained in a quiescent state for most of the duration of S-CUBED monitoring. XRT first detected this source during weekly S-CUBED monitoring in March of 2020. Two more detections occurred in July and September of 2021. This mini-flare of X-ray photons went un-reported due to its sparse detection frequency. The X-ray luminosity of the flare peaked at an XRT count rate of $0.071_{-0.035}^{+0.051}$ counts s${}^{-1}$ before returning to a quiescent state. The peak XRT count rate implies an approximate 0.3-10 keV band luminosity of $L_{X}=4.76\times 10^{36}$ erg s${}^{-1}$ if the standard distance of 62.44 kpc (Graczyk et al., 2020) to the SMC is assumed.

The source was identified as a possible candidate BeXRB by examining the optical/IR SEDs of all unidentified X-ray sources in the S-CUBED survey (Gaudin et al. in prep.). This identification motivated a deep 10ks Swift observation in August of 2023 to try and constrain the hardness ratio of the source’s X-ray spectrum. However, the source was only weakly detected with a mean luminosity of $L_{X}=9.07\times 10^{35}$ erg s${}^{-1}$.

The results of all XRT observations during the recent outburst are shown in Figure 1. Weekly S-CUBED monitoring detected X-ray emission from the source as it entered outburst on October 10th and again at peak XRT count rate on October 31st. The peak XRT count rate of Swift J010902.6-723710 was $0.22_{-0.058}^{+0.071}$ counts s${}^{-1}$, implying a peak 0.3-10 keV band luminosity of $L_{X}=1.35\times 10^{37}$ erg s${}^{-1}$ at a distance of 62.44 kpc. S-CUBED detection of outbursting behavior triggered deeper follow-up Swift observations of 5 ks per day on November 2nd and 4th in an effort to both monitor the outburst and constrain the spin period of the NS. Deep observations show the XRT count rate start to decline exponentially as expected, reaching a mean XRT count rate value of 0.11 counts s${}^{-1}$ during the early November observations. Additional observations taken on November 22nd and December 2nd confirm the continuing trend of exponential decline in X-ray brightness as the outburst fades. However, weekly S-CUBED detections persisted, with additional detections occurring on December 12th and December 26th. These two detections had 0.3-10 keV luminosities of $L_{X}=3.55\times 10^{36}$ erg s${}^{-1}$ and $L_{X}=3.98\times 10^{36}$ erg s${}^{-1}$, respectively, which do not fit the trend of steadily declining flux values after the peak luminosity was reached.

Figure 2: $Z_{2}^{2}$ Periodogram for four observations with exposures of $\simeq$5 ks of Swift XRT data taken in PC mode. The periodogram shows a peak at 182 seconds in all four observations made after outburst was detected by S-CUBED weekly monitoring. There is evidence for pulsar spin-up over time, however the errors on period measurement are large meaning this result is not statistically significant.

All 5 ks follow-up Swift observations were searched for the presence of periodic pulsations that are indicative of the spin period of the X-ray pulsar. The results of this search, carried out using a $Z_{2}^{2}$ Periodogram (Buccheri et al., 1983), are shown in Figure 2. In all 5 ks observations, the periodogram shows peaks at around $182$s (see Figure 2 for individual period measurements). Pulsar periods show small evidence for spin-up during the outburst, although given large errors on the periodicities due to the low numbers of counts and low (2.5s) PC mode frame time, the spin-up is not strongly detected. Folding the light curve of each individual XRT observation of 182s reveals an asymmetric double-peaked periodic signature, which is the expected shape for an X-ray pulsar light curve. Based on the phase-folded light curve and the strength of the periodogram peaks, we argue that the true spin period of the pulsar is 182s, allowing Swift J010902.6-723710 to be designated as SXP 182.

Spectral fitting is used to derive values for the spectral index of this source and column density along the line of sight. This is done by fitting the time-averaged 0.3-10 keV spectrum to an absorbed power-law model using the methods presented in Evans et al. (2014). Using these methods, best-fit parameters and a 68% confidence interval can be derived for both properties. SXP 182 is shown to have a hard X-ray spectrum, with a derived photon index for the source of $\Gamma=0.52_{-0.15}^{+0.16}$, which is consistent with the photon index of $\Gamma\sim 0-1$ expected of a BeXRB system (Kennea et al., 2018). Additionally, power law fitting indicates that the column density along the line of sight towards this source is $N_{H}=5.1_{-1.9}^{+2.3}\times 10^{21}$ cm${}^{-2}$ which is higher than the average value of the column density of $N_{H}=5.34\times 10^{20}\text{cm}^{-2}$ towards the SMC (Willingale et al., 2013; Kennea et al., 2018).

Figure 3: OGLE IV, Swift UVOT, and Swift XRT light curves for Swift J010902.6-723710.

2.2 Swift - UVOT

Swift J010902.6-723710 has been observed 221 times by UVOT in the uvw1-band since 2016 as part of S-CUBED monitoring, which represents approximately weekly coverage. UVOT light curves are generated for S-CUBED sources using the FTOOLS software package (Blackburn et al., 1999). Photometric data are extracted from a circle with 5 arcsecond radius around the XRT source position for the object using the uvotsource method that is part of FTOOLS. Using this method, the light curve presented in Figure 3 was generated for the entire duration of S-CUBED.

UVOT data identifies Swift J010902.6-723710 as a persistent emitter in the uvw1 band as it is well-detected in every observation with a 14.4 mean magnitude. The most prominent feature of this light curve is the presence of strong UV variability on timescales of 500-1000 days. The source reached its minimum average magnitude of 14.7 in October 2021 after years of gradual dimming from its initial brightness of 14.1 at the start of S-CUBED monitoring. At this point, the source began a period of rapid brightening lasting for over 500 days, reaching a maximum magnitude of 13.7 on June 27, 2023. After the peak magnitude was reached, the source has entered a steep dimming phase that corresponds with the start of the XRT outburst. This dimming phase was observed to be concluded on MJD 60283, when the source was observed to be at its pre-outburst uvw1 magnitude.

The long period of UV brightening that Swift J010902.6-723710 underwent was crucial in identifying the system as a candidate BeXRB before the outburst occurred. This type of UV increase leading to outburst has been observed in other sources such as Swift J004516.6–734703 (Kennea et al., 2020) which demonstrated very little UV variability before experiencing a similar brightening in the lead up to a Type I outburst. In the outburst of Swift J004516.6–734703, the UV brightening was interpreted to be the result of a circumstellar disk forming around the companion Be star in the newly-identified binary. For Swift J010902.6-723710, we can interpret the long-term variability to be an indication of growth and decay in the size of the disk, where the disk expanded over the last 500+ days to a large-enough radius for NS-disk interactions to occur.

2.3 OGLE - the Optical Gravitational Lensing Experiment

The OGLE project (Udalski et al., 2015) provides long term I-band photometry with a cadence of 1-3 days for sources in the Magellanic Clouds. From the X-ray position the optical/IR counterpart to Swift J010902.6-723710 is identified as 2MASS J01090226-7237101. It was observed continuously for nearly 3 decades by OGLE, though the observations were interrupted for approximately 3 years during COVID-19.

The counterpart to the X-ray source is identified in the OGLE catalogue as:

OGLE II (I-band): smc_sc11.107571
OGLE III (I-band): smc110.3.22311
OGLE IV (I band): smc726.26.15515
OGLE IV (V band): smc726.26.v.22358

The I band data from just the OGLE IV project are shown in the top panel of Figure 3 for comparison with the contemporaneous Swift UVOT and XRT data.

Figure 4: Top panel - the whole OGLE data set covering 26 years. Lower panel - the above data detrended with a simple polynomial in preparation for timing analysis.

Figure 5: Generalised Lomb-Scargle analysis of the detrended OGLE data set. The peak is at a period of $60.623\pm 0.001$ days.The horizontal dashed line represents the False Alarm Probabilty of 1%.

Figure 6: Detrended OGLE data divided into 5 consecutive epochs (see Table 1) and each then folded using the ephemeris given in Equation 1. For the purpose of being able to see the separate profiles, those numbered 2-4 have been arbitrarily shifted in the y-axis such as to avoid overlap.

Table 1: OGLE data date ranges used in Figure 6

Data block	Dates	Average I band
number	JD-2450000	magnitude
1	600-2000	16.26
2	2000-3100	16.22
3	3100-5000	16.35
4	5200-6400	16.37
5	6400-9000	16.34

For the purpose of period analysis the whole of the 26.4 years of OGLE data were used and first detrended using a simple polynomial - see Figure 4. This data set was then analysed with a generalised Lomb-Scargle routine and the resulting power spectrum is shown in Figure 5. The peak in the power spectrum is at a period of $60.623\pm 0.001$ days. This peak is driven by the sharp eclipse-like features that can be clearly seen by eye in the top panels of Figures 3 and 4. It is assumed that this represents the orbital period of the system.

If all the detrended OGLE data are folded at the proposed binary period, with a phase 0.0 set to the date of the first OGLE measurement (JD 2450627.9) then the sharp eclipse-like dips are clearly seen in Figure 6. The average FWHM of the dips is 0.1 in phase. It is also clear that the ingress into the eclipse is much less sharp than the egress. To explore the changes in the shape of the eclipse over time the OGLE data were divided up into 5 separate time segments representing the times when visible changes were occurring in the profile. The resulting separate profiles are shown in Figure 6. The 5 epochs chosen are listed in Table 1. Using the phase of the dip in the I-band shown in Figure 6 as a reference point, the ephemeris for the time of the optical eclipses is given by Equation 1:

$$T_{ecl}=2450645.1+N(60.62\pm{0.01})~{}\textrm{~{}JD}$$ (1)

In addition to the regular I-band measurements the OGLE IV project also records V band magnitudes every few days. This provides the opportunity to investigate the overall colour changes seen in the system as function of brightness - a colour-magnitude diagram (CMD). Since the V band measurements are less frequent than the I band, the determination of (V-I) can only occur when the I and V measurements are close enough together in time. In this instance the proximity of the 2 measurements is set to be less than 3 days. This means that occasionally one V band measurement may be partnered with more than one I band measurements. The result is shown in Figure 7 and discussed below.

2.4 SALT - the Southern African Large Telescope

Figure 7: $(V-I)-I$ (left) and $(U-I)-I$ (right) colour-magnitude diagrams of Swift J010902.6-723710

Table 2: The H$\alpha$ equivalent width (EW) and peak separation ($\Delta V$) measured from SALT observations.

Date	MJD	EW (Å)	Grating	$\Delta V$ (km/s)
03-11-2023	60251.88	-7.4$\pm$0.3	PG0900	$-$
09-11-2023	60257.87	-7.3$\pm$0.5	PG2300	373.6$\pm$0.4
02-12-2023	60280.85	-7.4$\pm$0.4	PG2300	383.4$\pm$0.9

Swift J010902.6-723710 was observed with the Southern African Large Telescope (SALT; Buckley et al. 2006) with the Robert Stobie Spectrograph (Burgh et al., 2003; Kobulnicky et al., 2003) using different settings. On 03-11-2023 (MJD60251.88) the PG0900 grating was used (grating angle of 15.125${}^{\circ}$) with an exposure time of 1200 sec covering a wavelength range of $4200-7250$ Å. The PG2300 grating was used on 09-11-2023 (MJD60257.87) and 02-12-2023 (MJD60280.85) (grating angle of 48.875${}^{\circ}$) with an exposure time of 1200 sec covering a wavelength range $6100-6900$ Å. An additional observation was taken on 09-11-2023 (MJD60257.87) with the PG2300 grating (grating angle of 30.5${}^{\circ}$) with an exposure time of 1500 sec covering a wavelength range $3840-4915$ Å. The primary reductions, which include overscan corrections, bias subtraction, gain and amplifier cross-talk corrections, were performed with the SALT science pipeline (Crawford et al., 2012). The remainder of the data reduction steps, which comprise wavelength calibration, background subtraction, and extraction of the one-dimensional spectrum, were done in iraf. All spectra were corrected for the heliocenter and redshift of the SMC of 145.6 km/s (McConnachie, 2012).

Figure 8 shows the H$\alpha$ emission line profiles. The observations obtained on MJD60257.87 and MJD60280.85 with the PG2300 grating exhibit asymmetric double-peak profiles as a result of the Keplerian distribution of matter in the disc when the disc is viewed at non-zero inclination angles. The observation performed on MJD60251.88 with the PG0900 grating shows an asymmetric single-peak profile since the resolution of the grating is insufficient to resolve the two peaks. The H$\alpha$ equivalent width measurements are recorded in Table 2.

3 Discussion

3.1 Corbet Diagram

Combining the spin period with the orbital period allows this new system to be placed on the Corbet diagram. The original Corbet diagram (Corbet, 1984) showed a correlation between the spin and orbital periods of neutron stars in BeXRBs. However, instead of a correlation between spin and orbital periods, a modern version (shown in Figure 9) shows that these properties are well-constrained to a specific region of the Corbet diagram. BeXRB systems are expected to fall above the diagonal of the diagram, which represents a log-linear relationship between the orbital period in days and the spin period in seconds of the NS. Figure 9 shows the location of SXP 182 when placed on the Corbet diagram with all BeXRBs in the SMC that have known spin and orbital periods (Haberl & Pietsch, 2005; Coe & Kirk, 2015; Haberl & Sturm, 2016; McBride et al., 2017; Carpano et al., 2017; Kennea et al., 2018; Lazzarini et al., 2019; Kennea et al., 2020; Carpano et al., 2022; Maitra et al., 2023). The location of SXP 182 is consistent with the trend shown by other SMC BeXRBs. This serves as a check that the periods derived in Sections 2.1 and 2.3 are consistent with the observational properties of other BeXRBs. For Swift J010209.6-723710, the derived orbital period of 60.623 days and pulsar spin period of 182 seconds place the system near the center of the distribution of orbit and spin periods for SMC BeXRBs, providing a strong piece of evidence validating the source as a newly-discovered BeXRB.

3.2 Inclination angle of the Be disc

The peak separation of the double-peak H$\alpha$ emission lines (Figure 8) can be used to estimate the size of the H$\alpha$ emitting region (Huang, 1972):

$$R=\frac{GM_{\ast}\sin^{2}i}{(0.5\Delta V)^{2}}$$ (2)

Using $M_{\ast}=17.8$ M${}_{\odot}$ (Cox, 2000) based on the spectral type of the massive companion, this results in a disc size radius range of $92-97\sin^{2}i$ R${}_{\odot}$. We can estimate the inclination angle of the Be disc, assuming that the disc and orbital planes are aligned and that during the period of X-ray activity the NS was accreting matter from the outermost parts of the Be disc at periastron passage. The semi-major axis of the orbit can be estimated by assuming a companion mass of $M_{\ast}=17.8$ M${}_{\odot}$ and a neutron star mass of $M_{NS}=1.4$ M${}_{\odot}$. Similarly, the periastron passage can be estimated using a conservative value for the eccentricity of $e\sim 0.5$ based on the relationship between eccentricity and orbital period presented in Townsend et al. (2011). This results in a range of disc inclination angles of $72-90^{\circ}$. The suggestion of a high disc inclination angle is corroborated by the H$\alpha$ emission line displaying a double peak morphology with a deep central depression from the high-resolution observations.

Figure 8: The evolution of the H$\alpha$ emission line from SALT observations. The spectra are corrected for the heliocenter and redshift of the SMC.

Figure 7 shows the $(V-I)-I$ and $(U-I)-I$ colour-magnitude diagrams. The general trend of the colour-magnitude plots is a redder-when-brighter pattern, which is indicative of inclination angles below $90^{\circ}$(Harmanec, 1983; Rajoelimanana et al., 2011; Reig & Fabregat, 2015). This trend is more noticeable in the $(V-I)-I$ plot where the range in colour is broader since the simultaneous $V$ and $I$ band observations were taken during a period of substantial optical variability.

Figure 9: Corbet diagram for all BeXRBs in the SMC with known orbital periods and spin periods (Haberl & Pietsch, 2005; Coe & Kirk, 2015; McBride et al., 2017; Carpano et al., 2017; Kennea et al., 2018; Lazzarini et al., 2019; Kennea et al., 2020; Carpano et al., 2022; Maitra et al., 2023). SXP 182 is plotted alongside these data and is shown to demonstrate the same correlation between X-ray pulsar spin period and NS orbital period that is seen in the rest of the SMC BeXRB population.

3.3 Spectral Classification of the massive counterpart

Figure 10 shows the SALT spectrum of Swift J010902.6-723710 covering the blue wavelength range. The spectrum shows clear signatures of an early-type star with several Balmer and helium lines present. The H$\beta$ line is in weak emission with an asymmetric profile, along with the other absorption lines exhibiting infilling due to the presence of the circumstellar disc. According to the criteria presented in Evans et al. (2004), the strong presence of the He I lines at 4026, 4143, 4388 and 4471 Å  as well as the weak presence of the He II lines at 4541 and 4686 Å  constrains the spectral type to B0-0.5. The faintest V-band observation from OGLE measurements during our period of monitoring is $\sim$16.5 magnitudes. Using this and the distance modulus of the SMC of 18.95 (Graczyk et al., 2013), this results in a luminosity class of V (Straizys & Kuriliene, 1981; Pecaut & Mamajek, 2013). In summary, the spectral class of the massive companion in Swift J010902.6-723710 is B0-0.5 V.

Figure 10: The SALT spectrum of Swift J010902.6-723710 covering the blue region with different line species labeled at their expected rest wavelengths. The spectrum is corrected for the heliocenter and redshift of the SMC.

3.4 Eclipsing Behavior

The orbital period of Swift J010902.6-723710 derived from long-term OGLE monitoring can be used to re-examine the UVOT light curve in search of periodic behavior. This was done by de-trending the UVOT data using a 5th-order polynomial and folding the de-trended data at the 60.623 day orbital period of the system. Figure 11 shows the comparison between the binned and time-averaged OGLE and UVOT light curves starting at MJD 57563, which is the date of the first S-CUBED observation. Despite S-CUBED’s weekly cadence providing infrequent sampling of any singular orbit, the eclipsing behavior is clearly visible in the de-trended and folded UVOT light curve. When compared to the OGLE data, it becomes evident that the shape of the eclipse profile is wavelength-dependent. The UVOT data shows increasing emission that peaks just before each eclipse begins, which is similar to the behavior exhibited by the system during the 4th OGLE data block shown in Figure 6. Over the lifetime of S-CUBED, the UVOT eclipse profile is also shown to be broader than the OGLE eclipse profile with similar depth at the time of maximum eclipse. At the maximum depth of the eclipse, the uvw1-band magnitude drops by an average of 0.22 magnitudes and I-band magnitude decreases by an average of 0.25 magnitudes. Converting these magnitude decreases to flux measurements for both telescopes, the implied relative decrease in flux is at both wavelengths $\frac{\Delta F}{F}\sim 0.2$.

Figure 11: A comparison of the OGLE IV and UVOT light curves folded at the proposed binary period of 60.623 days using Equation 1. The OGLE data are offset vertically to avoid overlap.

Equation 2 from Maggi et al. (2013) can be used to estimate the size of the accretion disk based on the relative decrease in flux that is observed during eclipse:

$$\frac{\Delta F}{F}=\left(\frac{R_{X}}{R_{C}}\right)^{2}$$ (3)

where $R_{X}$ is the radius of the eclipsing object and $R_{C}$ is the radius of the optical companion. The value for $R_{X}$ can be estimated due to the constraints placed on $R_{C}$ via the results of Section 3.3. The typical radius for a spectral class B0V star is 7.4 $R_{\odot}$ (Allen, 1976). This value can be used to constrain the upper limit for the radius of the eclipsing object at $R_{X}=R_{C}\sqrt{\frac{\Delta F}{F}}=3.3\,R_{\odot}$. A radius of this size rules out the NS, a large planet, or a Sun-like tertiary star as the cause of the eclipsing behavior. Therefore, it can be concluded that the eclipsing object is likely an extended accretion disk that surrounds the NS. Maggi et al. (2013) derives the theoretical upper limit for the size of an accretion disk in a BeXRB system to be $r_{c}\sim 11.5\,R_{\odot}$, assuming Bondi-Hoyle accretion and a 10 $M_{\odot}$ companion star. The 3.3 $R_{\odot}$ accretion disk size derived for Swift J010902.6-723710 is well within this upper limit value, providing further evidence in favor of an eclipsing NS accretion disk.

Very few BeXRBs have been observed to demonstrate eclipsing behavior. Swift J010902.6-723710 is only the third known eclipsing BeXRB system and the second to be found in the SMC. LXP 168.8 (Maggi et al., 2013) was found to have a eclipsing accretion disk with a 24.3 day orbital period. SXP 5.05 (Coe et al., 2015) was observed to contain a Be star that eclipses its NS companion every 17 days. This newly-discovered eclipsing system thus provides a unique opportunity to further constrain the physical parameters of the system such as the sizes of both disks and the masses of both stars. More observations are needed to further characterize the BeXRB using subsequent eclipses.

3.5 Outburst Type

Date	MJD	Phase	$L_{X}$ (erg s${}^{-1}$)
10-10-2023	60227	0.06	$4.14_{-2.61}^{+4.27}\times 10^{36}$
10-31-2023	60248	0.40	$1.35_{-0.391}^{+0.475}\times 10^{37}$
11-02-2023	60250.0	0.47	$5.66_{-0.275}^{+0.275}\times 10^{36}$
11-02-2023	60250.9	0.48	$7.06_{-0.527}^{+0.527}\times 10^{36}$
11-04-2023	60252.1	0.50	$5.30_{-0.508}^{+0.508}\times 10^{36}$
11-05-2023	60253.0	0.52	$4.96_{-0.253}^{+0.253}\times 10^{36}$
11-05-2023	60253.6	0.53	$4.03_{-0.385}^{+0.385}\times 10^{36}$
11-06-2023	60254	0.54	$4.95_{-0.518}^{+0.518}\times 10^{36}$
11-21-2023	60269.9	0.80	$3.20_{-2.04}^{+3.33}\times 10^{36}$
11-22-2023	60270.6	0.81	$2.25_{-0.190}^{+0.190}\times 10^{36}$
12-02-2023	60280	0.96	$1.83_{-0.177}^{+0.177}\times 10^{36}$
12-08-2023	60286	0.06	$1.43_{-0.192}^{+0.192}\times 10^{36}$
12-12-2023	60290	0.13	$3.55_{-2.28}^{+3.78}\times 10^{36}$
12-14-2023	60292	0.16	$1.17_{-0.397}^{+0.397}\times 10^{36}$
12-26-2023	60304	0.36	$3.98_{-2.26}^{+3.69}\times 10^{36}$

Table 3: S-CUBED and Swift TOO XRT detection dates, orbital phases from Equation 1, and luminosities during outburst.

Figure 12: Truncated OGLE IV, Swift UVOT, and Swift XRT light curves for Swift J010902.6-723710 showing emission from MJD 60100 to MJD 60300. The times of optical eclipse, calculated using Equation 1, are shown as vertical grey lines. Blue arrows show the calculated XRT flux upper limit during S-CUBED observations where Swift J010902.6-723710 is not detected.

The type of outburst observed in Swift J010902.6-723710 has remained uncertain as long-term monitoring of the source continues to produce XRT detections well past the peak of the initial outburst event. In Type I outbursts, X-ray emission is often limited both by the duration of an orbit and the current orbital phase. Emission is expected to be detected at a narrow range of values near the periastron passage of the orbit and is expected to last for no longer than a full orbital period (Stella et al., 1986). Additionally, Type I outbursts typically demonstrate a moderate increase in luminosity, peaking at $L_{X}\sim 10^{36}-10^{37}$ ergs s${}^{-1}$ (Okazaki & Negueruela, 2001; Reig, 2011). Table 3 shows the date of XRT detections of Swift J010902.6-723710 and the orbital phase at which they occurred with respect to Equation 1. The light curve of XRT detections during outburst is plotted with respect to the time of optical eclipses in Figure 12. Both the figure and table show the unusually long duration of this outburst. The first S-CUBED detection of the outburst on October 10th occurred at Phase 0.06, which is approximately the periastron passage of orbit. However, other subsequent detections continue to occur well past the periastron passage of the orbit, demonstrating phase independence of the emission. The last two S-CUBED detections, occurring on December 12th and December 26th, correspond to Phase 0.13 and 0.36 of a new orbit. These subsequent detections are not consistent with the typical behavior demonstrated by a Type I outburst and have more in common with the emission signature of a Type II outburst (Okazaki & Negueruela, 2001; Townsend et al., 2017; Tamang et al., 2022). If this is indeed a Type II outburst, then the peak luminosity is much smaller than is typical for these outbursts. The peak X-ray luminosity of $L_{X}=1.35\times 10^{37}$ erg s${}^{-1}$ is the only detection to occur above $10^{37}$ erg s${}^{-1}$. Based on these features, one of two situations has occurred. One possibility is that the system has produced an abnormally long-duration Type I outburst with a typical peak luminosity (Okazaki & Negueruela, 2001; Reig, 2011). Alternatively, the system has produced a Type II outburst that fails to reach the characteristic large peak luminosity that is expected during these events (Okazaki & Negueruela, 2001; Reig, 2011).

However, it is important to remember that this traditional classification scheme is very quantised - either it is called Type I or Type II. Since originally proposed in 1986 (Stella et al., 1986) there have been many more diverse examples of X-ray outbursts from BeXRB systems observed. In reality there must be a scope for a whole spectrum of outburst types since the outburst duration and the X-ray luminosity seen depend almost entirely on the interaction between the circumstellar disc and the neutron star. That in turn depends upon the characteristics of the circumstellar disc (density, size, inclination, mass outflow rate etc.) and the characteristics of the neutron star orbit (eccentricity, inclination to the circumstellar disk plane, phase of periastron etc.). So, given all those free parameters one should expect a whole range of observational characteristics, and not be driven to simply call them Type I or Type II. SXP 182 is an excellent example of how the original classification scheme can be overly simplistic.

4 Conclusion

This paper reports the detection of a previously-unknown BeXRB via weekly observations of the S-CUBED survey. This new system, Swift J010902.6-723710, was identified via a transient X-ray outburst and followed up via multi-wavelength observations. Deep follow-up X-ray observations identify a proposed spin period of 182s for the NS in this binary system. Historical light curve analysis of both UV and IR emission reveal the presence of strong eclipse-like features that re-appear every 60.623 days in the light curve, which is adopted as the proposed orbital period of the system. Optical spectroscopy reveals a strongly double-peaked H$\alpha$ emission line, indicating a highly-inclined system with an inclination of $i=72-90^{\circ}$. Spectroscopic observations are also used to constrain the spectral class of the optical companion as a B0-0.5 star of spectral class Ve.

The proposed spin and orbital periods place Swift J010902.6-723710 place in the center of the expected distribution for similar BeXRBs on the Corbet Diagram. Eclipsing behavior is found to be caused by a $3.3R_{\odot}$ accretion disk that surrounds the NS, making this the third eclipsing BeXRB to be detected so far. Finally, the type of outburst observed is found to be uncertain, with characteristics of both Type I and Type II outbursts being found in the X-ray emission of the source. More observations of this system are needed, particularly during subsequent eclipses, in order to further constrain the physical parameters of this rare eclipsing system and better understand this long-lasting, moderately luminous X-ray outburst.

5 acknowledgments

This work made use of data supplied by the UK Swift Science Data Centre at the University of Leicester. JAK and TMG acknowledge the support of NASA contract NAS5-00136. We acknowledge the use of public data from the Swift data archive.

The SALT observations reported in this paper were obtained under the SALT Large Science Programme on transients (2021-2-LSP-001; PI: DAHB), which is also supported by Poland under grant no. MNiSW DIR/WK/2016/07.

PAE acknowledges UKSA support. DAHB acknowledges research support from the National Research Foundation. LJT is supported by the SALT Foundation and the NRF.