1. INTRODUCTION
As described by the standard model (e.g. Warner & Nather 1971) cataclysmic variables (CVs) are semi-detached binary systems consisting of a late- type secondary star filling its critical Roche surface, that transfers matter into a white dwarf (WD) via an accretion disk (see Warner 1995, for a comprehensive review on CVs). CVs can be divided into systems with outbursts, like dwarf novae, and non-outburst systems, like intermediate and polar systems (e.g. Hellier 2001). Intermediate polar systems (IPs) are a class of CVs whose primary component is a magnetic white dwarf. The presence of a moderate magnetic field (0.1-10 MG) in these systems is strong enough to inhibit the formation of the innermost regions of the accretion disk, but is not sufficiently strong to synchronise the rotation of the white dwarf with the orbit (e.g. Hameury & Lasota 2017). The number of intermediate polars is small (about one percent) compared with the total population of CVs (Warner 1995).
DW Cancri (hereafter DW Cnc) was identified as a CV by Stepanyan (1982). Rodríguez-Gil et al. (2004) performed its first time-resolved spectroscopic study. Their analysis led them to suggest that the system is a short-period intermediate polar below the period-gap, whose photometric data resemble the behaviour of VY Scl stars. They found the Balmer and He I lines to be modulated with two periods: 86.10±0.05 min, associated with the orbital period; and 38.58±0.02 min, likely corresponding to the WD spin period. Patterson et al. (2004) confirmed the intermediate polar nature of DW Cnc. A radial velocity search performed by these authors detected the 86 and 38 min periods. Additionally, in a photometric power search they found a 70 min signal, consistent with the beat period of the binary (compelling evidence of its IP nature), and also a weak periodic signal at 110 min, which was left as an unsolved problem. Nucita et al. (2019) reported positive XMM-Newton observations in 2012 in the range 0.3 - 10 Kev; their light curves show evidence of a period around 38 min and also a signature around 75 min, both consistent within the errors to the spin and beat periods, respectively.
In a previous publication (Segura Montero et al. 2020), a radial velocity study of DW Cnc was presented, with observations performed during a low state in 2018-2019. Through a power spectrum analysis they found the 86 min signal associated with the orbital period and two much weaker modulations associated with the 70 min beat period and the 38 min spin period. Particularly, the 38 min signal was significantly weaker than that previously published by Rodríguez-Gil et al. (2004) and Patterson et al. (2004). To explain this substantial change Segura Montero et al. (2020) suggested that the sudden drop into a low state-caused by an episode of low mass transfer from the companion-inhibited the lighthouse effect produced by the rebound emission, thus rendering the spin period of the WD difficult to detect. Such variability in the behaviour of DW Cnc exemplifies the importance of pursuing follow-up observations of the system.
Hence, in § 2 of this paper we present new photometric and spectroscopic observations of DW Cnc. In § 3 we show a radial velocity study of the Hα and He I 5876 Å emission lines. We performed a power spectrum search of our photometric and spectroscopic data, which is shown in § 4, followed by a wavelet transform analysis in § 5. We also carried out a Doppler tomography study in § 6. We close the article with a discussion of our results and our conclusions in § 7 and § 8, respectively.
2. OBSERVATIONS AND REDUCTION
2.1. Photometry
CCD photometry was obtained on 2020 March 8-9 and 15-16, with the 0.84m telescope at the Observatorio Astronómico Nacional at San Pedro Mártir, located in Baja California, Mexico. V and R images were obtained sequentially during the four nights with an e2vm E2V-4240 2048×2048 CCD using 2×2 binning. The exposure times were of 30s and 20s for the V and R filters, respectively. Data reduction was carried out with the IRAF1 software system. After bias and flat field corrections, aperture photometry of DW Cnc and some field stars was obtained with the phot routine. The same comparison star (RA=7:58:58, DEC=+16:15:07) used by Patterson et al. (2004) was adopted, assuming the reported magnitudes of B=15.89, V =15.21, and R=14.82. The log of photometric observations is shown in Table 1.
2.2. Spectroscopy
Spectra were obtained with the 2.1m telescope of the Observatorio Astronómico Nacional at San Pedro Mártir, using the Boller and Chivens spectro-graph and a Spectral Instrument CCD detector in the 5500 - 6500 Å range (resolution ≈ 1200), on the nights of 2020 March 15 and 16. The exposure time for each spectrum was 300 s. Standard IRAF procedures were used to reduce the data. The log of spectroscopic observations is shown in Table 2. The spectra show strong Hα λ6563 Å and He I λ5876 Å emission lines. The typical S/N ratio of the individual spectrum is of ≈20 for the emission lines. The spectra were not normalized.
3. RADIAL VELOCITIES
The radial velocity of the emission lines in each spectrum were computed using the RVSAO package in IRAF, with the CONVRV function, constructed by J. Thorstensen (2008, private communication). This routine follows the algorithm described by Schneider & Young (1980), convolving the emission line with an antisymmetric function, and assigning the zero value of this convolution as the midpoint of the line profile. As in Segura Montero et al. (2020), we initially used the GAU2 option, available in the routine, which uses a negative and a positive Gaussian to convolve the emission line, and needs the input of the width and separation of the Gaussians. This method traces the emission of the wings of the line profile, presumably arising from the inner parts of the accretion disc.
Following the methodology described by Shafter et al. (1986), we made a diagnostic diagram to find the optimal Gaussian separation, by performing a non-linear least-squares fit (Newville et al. 2014) of a simple circular orbit to each trial:
where y is the systemic velocity, K
1 the semi-amplitude, t
0 the time of inferior conjunction of the donor and P
orb is the orbital period. We employed
In particular, a control parameter is defined in this diagnostic, σ
K
/K, whose minimum is a very good indicator of the optimal fit. The diagnostic diagram for Hα is displayed in Figure 1, while the orbital fit for its best solution is exhibited in Figure 2. The diagnostic diagram of the He I λ 5876 Å emission line is shown in Figure 3, and its orbital fit appears in Figure 4. The 1σ error bars of the radial velocity fits were scaled so that the goodness-of-fit
Parameter | Hα | He I λ5876 |
γ (km s−1) | 10 ±2 | 40 ±2 |
K1 (km s−1) | 21.5 ±2.5 | 27.3 ±2.9 |
HJD0* | 0.803 ±0.001 | 0.833 ±0.001 |
P orb (min) | Fixed** | Fixed** |
1Using the wings of the lines.
*(24558923+ days).
**86.10169 ± 0.00031 min.
Furthermore, as explained in § 4.2, following Patterson et al. (2004), we also implemented the DGAU convolution option to perform an additional power spectrum. This option uses the derivative of a single Gaussian, and only requires the input of the Gaussian width.
4. POWER SPECTRUM SEARCH
We made a power spectrum analysis of the V and R photometric bands, and a power search of the measured radial velocities of the Hα and the He I λ5876 Å emission lines, using a Lomb-Scargle algorithm (Scargle 1982) in both cases.
For each periodogram, we computed the false- alarm-probability (hereafter FAP), using a function included in Astropy Collaboration et al. (2013). As explained by VanderPlas (2018), the FAP quantifies the significance of a peak by calculating the probability that the random variations in the data lead to a peak of similar magnitude, conditioned on the assumption of the null hypothesis of having no periodic signal present in the data. Following de Lira et al. (2019), we identified as significant peaks those whose FAP was less than 0.01, i.e., with a significance level greater than 99 percent.
4.1. Photometric Data
In Figure 5 we show the results of the photometric frequency analysis. The results from the V band are in the upper panel, where we observe a prominent peak at a frequency of 37.47 cycles/day (associated with the 38 min spin period). Two other weak signals are also present: 13.26 and 23.25 cycles/day, which correspond to periods of 108 min and 62 min, respectively. We do not find the orbital period, and contrary to Patterson et al. (2004) and Nucita et al. (2019), our periodogram does not detect the 70 min spin-orbit beat period.
The results from the R band (middle panel) show the same signals as the V band. In the bottom panel, we show the combined analysis using both the V and R bands. The results are the same as before, where no significant power signals are found for the orbital and beat periods. We note here that the 23.25 signal did not yield a FAP below the 0.01 threshold. To probe for its legitimacy and make an attempt to unmask new signals, we followed Patterson et al. (2004), by fitting a sine wave with a periodicity corresponding to the predominant 37.47 cycles/day signal, and then subtracting the fitted sinusoid from the photometric time series. We proceeded to make a new Lomb-Scargle search of the residuals as exhibited in Figure 6. It can be observed that after subtracting the spin cycle, both of the 13.26 and 23.25 signals remain present, with their recalculated FAPs below the 0.01 cutoff value. This result is different from that obtained by (Patterson et al. 2004), whose power search yields a weak bump at 16 cycles/day after subtracting the spin and beat modulations; however, this detection weakens when they study the long-term behaviour of their data.
Before folding the photometric time series by the modulations found in the power spectra, we made an analysis to check whether the V and R data are correlated. For this purpose we used a code2 developed by Figueira et al. (2016), which implements a Bayesian approach to produce the probability distribution of the correlation coefficient ρ. The V and R data yielded a distribution with a mean value of 0.367, a standard deviation of 0.034, and with a 95% credible interval of [0.298,0.431]. The lower limit of the 95% credible interval is well above ρ = 0, which establishes sufficient confidence in the correlation.
In Figure 7 we show, from top to bottom, the photometric data folded by the associated spin period, the 62-min period, the associated orbital period and the 108-min period. To reduce the influence of noise we averaged the data into 60 phase bins. Folding by the spin period clearly depicts a strong modulation, as expected from the power spectrum analysis. The 62-min and 108-min folded data show a slightly noisier, yet clearly perceptible, sinusoidal oscillation. It is worth noting that the 62 min signal corresponds to the beat period between the spin and 108-min periods. Finally, as expected from the lack of signal in the power search, folding by the orbital period yields no sinusoidal modulation for the photometry.
4.2. Spectroscopic Data
4.2.1. GAU2 Option
As described in § 3, we performed the measurement of the radial velocities using the GAU2 convolution option. We then proceeded to implement the Lomb-Scargle power search on the convolved data, which is shown in Figure 8. The upper panel is the analysis of Hα, which shows a peak frequency at 16.82 cycles/day, corresponding to the 86 min orbital period.
No significant power signals for the spin and beat periods were found. In the lower panel we repeated the analysis for He I λ5876 Å, which shows a peak frequency at 16.76 cycles/day, associated also with the 86 min orbital period. Again, the spin and beat period signals are not present. The mean value of the orbital period for both lines is 85.76 ± 0.15 min. This value is, within the errors, compatible with that measured by Segura Montero et al. (2020). Since we observed a smaller number of spectra than the previous authors, we have a larger error. Therefore we will adopt their value of P orb = 86.10169 ± 0.00031 min, throughout this paper.
Following the methodology used for the photometric data described in § 4.1, we subtracted the conspicuous orbital signal from both the Hα and He I λ 5876 Å data sets and performed a power search on the residuals (shown in Figure 9). The subtracted data of Hα (top panel), as with the DGAU option (see § 4.2.2), show the appearance of a weak picket fence around 20.70 cycles/day, with an overly high FAP value of 0.92. On the other hand, the residuals from the He I λ 5876 Å data (bottom panel) show a signal at 36.39 cycles/day, a frequency comparable to the spin modulation, but its 0.97 FAP value puts its legitimacy in doubt.
We folded the data by the orbital and spin periods, as shown in Figure 10. The upper panel exhibits a clear modulation when the data is folded by the orbital period. However, folding by the spin period (bottom panel) shows no evident indication of a periodic signal.
4.2.2. DGAU Option
Following Patterson et al. (2004) we measured the radial velocities by convolving the emission lines with the derivative of a Gaussian (DGAU option) of 83.8 Å(50 pixels) of width. As can be observed in Figure 11, performing a power search on the data yielded the spin and orbital periods for Hα (top panel). The He I λ5876 Å emission line (bottom panel) shows a clear signal of the orbital period and a weak power peak at the spin period. The computed FAP for this weaker spin peak yielded 0.94, rendering it as a questionable signal for He I λ5876 Å. This results are in agreement with those obtained by Patterson et al. (2004) in their spectroscopic power search, where they found the orbital and spin periods.
As before, we subtracted the strong modulations from the DGAU radial velocity data sets of both emission lines, and performed a new Lomb-Scargle search on the residuals. After subtracting the orbital signal from the He I λ 5876 Å data (bottom panel of Figure 12), the spin signal peak in the power search of its residuals shows a considerable increase, showing a new FAP value of 0.01 and establishing a better reliability upon this modulation. After subtracting both the orbital and spin signals from the Hα data (top panel of Figure 12), its residuals hint at a weak surge of the 20.70 cycles/day signal, consistent with the GAU2 option in § 4.2.1. However, the FAP of this signal yields 0.93, a value far too high. Nonetheless, we also detected this signature in the wavelet analysis in § 5, indicating that its presence could be real.
In Figure 13 we folded the data by the spin and orbital periods, where both signals show clear modulations of the DGAU radial velocities.
5. WAVELET TRANSFORM ANALYSIS
The wavelet transform is a method that applies the convolution of the signal with a set of wavelets to map the variations occurring in both the time and frequency domains (See Bravo et al. 2014, and references therein for a detailed formulation of the of the scales (or frequencies) that contribute most to the total energy of the signal (e.g. Benítez et al. 2010).
Following Bravo et al. (2014) and de Lira et al. (2019), we performed a wavelet transform analysis of the photometric lightcurve and the spectroscopic radial velocity data. We applied the continuous 6th order Morlet wavelet transform, implemented from the Python ObsPy package3 (Krischer et al. 2015). In Figures 14-19 we show the local and global power spectra of our data sets. As explained by Bravo et al. (2014), the local spectrum depicts the energy distribution in time-frequency space, and the global spectrum is obtained by time integration of the local map.
5.1. Photometry
The scalogram of the V -band (Figure 14) shows a predominant signal at 15.4 cycles/day. This signal is detected in all four nights, and its relative intensity increases on the third (HJD +2458923) and fourth night (HJD +2458924). This signature is consistent with the orbital cycle, which was absent in the Lomb-Scargle power spectrum.
The R-Band analysis (Figure 15) shows a similar structure to that exhibited by the V -band scalogram, except for the second night (HJD +2458917), which yields a strong signal at ≈21 cycles/day; this frequency is comparable to the 20.60 cycles/day modulation found by Patterson et al. (2004), corresponding to the beat period between the spin and orbital cycles. We note that we did not detect this signal in the photometric power Lomb-Scargle analysis (see § 4.1), but it did appear after subtracting the stronger signals in the spectroscopic data in § 4.2. A secondary signal peaking at ≈35 cycles/day (possibly associated with the spin period) is also visible during the second night (HJD +2458917).
5.2. Spectroscopy
We performed the wavelet analysis on the radial velocity data sets obtained both from the GAU2 and DGAU convolution methods. We now proceed to describe this results.
5.2.1. GAU2 Option
The Hα GAU2 scalogram in Figure 16, shows a broad power peak extending from ≈20 to ≈39 cycles/day at the outset of the first night (HJD +2458923.0). Such broad power peak narrows down into two localized peaks, of which the most prominent and persistent shows a midpoint at a frequency of ≈22 cycles/day; a signal consistent with that also found in the photometric Lomb-Scargle analysis of ≈23 cycles/day. The second night of this scalogram (HJD +245894) displays a strong signal at ≈33 cycles/day, which we consider to be a possible alias of the spin modulation. The ≈22 cycles/day signature is also present during the second night but with a decrease in relative intensity.
The He I λ 5876 Å GAU2 scalogram, exhibited in Figure 17, shows a prominent signal that persists throughout the first night (HJD +2458923) at ≈22 cycles/day, in agreement with the Hα gau2 data. During this first night a secondary signal appears at ≈40 cycles/day, which gradually shifts to ≈37 cycles/day (related to the spin period) as the night progresses. During the second night (HJD+2458924) we see an overall shift towards smaller frequency values, displaying a strong signal at ≈33 cycles/day and a slightly milder yet very persistent signal at ≈16 cycles/day (consistent with the orbital period).
5.2.2. DGAU Option
Figure 18 depicts the Hα DGAU scalogram, showing during the first night (HJD +2458923), a conspicuous signal at the expected spin period of ≈38 cycles/day. In the second night (HJD +2458924), the signal observed the previous night is still visible, but its structure considerably broadens in frequency, eventually adopting a two-pronged shape, whose contribution reflects on the loss of the 38 cycles/day signal in the global spectrum. The second night also shows a clear signal at ≈20 cycles/day, related to the 70 min beat period.
The He I λ 5876 Å DGAU scalogram, in Figure 19, displays a strong signal at ≈20 cycles/day that persists all throughout the first night (HJD +2458923), and a weaker signature at ≈38 cycles/day, consistent with the spin modulation. On the second night (HJD +2458924) the spin cycle signal becomes enhanced and its structure considerably broadens. This night also shows the appearance of a secondary detection at ≈16 cycles/day, related to the spin modulation, which persists over the whole night.
6. DOPPLER TOMOGRAPHY
Doppler tomography is an indirect imaging technique developed by Marsh & Horne (1988). It produces two-dimensional mappings of the emission intensity in velocity space of the accretion disc, using the phase-resolved profiles of the spectral emission lines. We obtained the Doppler tomography of the Hα and of the He I λ 5876 Å emission lines, using a Python wrapper4 (Hernandez Santisteban 2021) of the original FORTRAN routines published by Spruit (1998) within an IDL environment. In the top left panel of Figure 20 we show the observed trailed spectra of Hα, while the reconstructed trailed spectra appears in the top right panel; the tomogram is displayed in the bottom panel. With the same layout, the trailed spectra and tomography of He I λ 5876 Å are exhibited in Figure 21. The parameters used to plot the features in the tomograms are as follows: an inclination of i = 50o; a value of the mass of the primary star of M w =0.75 Mʘ consistent with the average mass for white dwarfs in CVs below the period gap (Knigge 2006); a mass ratio q = 0.2, estimated following Echevarria (1983); and an orbital period of P orb = 86.10 minutes (Patterson et al. 2004; Rodríguez-Gil et al. 2004; Segura Montero et al. 2020). We now proceed to describe the results obtained for each emission line.
6.1. Hα
The Hα observed trailed spectra display a complex behaviour. From orbital phase 0.0 to ≈0.10 they show a single peaked structure. In the interval from 0.10 to 0.30, the trailed spectra display a double-peaked profile. After this interval, the profile briefly becomes single peaked, and from 0.35 to 0.55 the blue shifted peak becomes more intense than the red shifted peak. From 0.55 onward the line profile again displays a symmetric double-peaked structure with a brief single-peaked intrusion at ≈ 0.8.
The tomography shows a disc signal in red colour (Marsh & Horne 1988) with a superimposed intense region (in black) at the position of the Roche Lobe of the secondary, which could be caused by emission from a hot spot component (e.g. Echevarría et al. 2007). The disc structure was not detected in the tomography by Segura Montero et al. (2020), but this finding is in good agreement with the dominant double-peaked structure, characteristic of discs in systems of high inclination (Horne & Marsh 1986), exhibited by the Hα line profiles reported by Rodríguez-Gil et al. (2004).
6.2. He I 5876 Å
The He I λ 5876 Å trailed spectra show the oscillation of a broad single-peaked profile, consistent with the line profiles of this emission line put forward by Rodríguez-Gil et al. (2004). The tomography shows a blob-like region of high intensity in the upper quad- rants that overlays the position of the Roche Lobe of the secondary, and further extends towards negative velocities, covering the position where the emission coming from the hotspot is expected in velocity space. Note that the position of the blob in the He I λ 5876 Å tomography is consistent with that of the region of maximum intensity observed for Hα.
7. DISCUSSION
We performed a study of new photometric and spectroscopic observations of DW Cnc, after the system recovered from a low state, presumably caused by an episode of low mass transfer from the secondary that inhibited the lighthouse effect (Segura Montero et al. 2020). Photometry from the AAVSO shows that at the time of our observations, DW Cnc had already reached a state comparable to that from 1999-2003 reported by Patterson et al. (2004) and Rodríguez-Gil et al. (2004). With this in mind, we performed various analyses to compare the behaviour displayed by DW Cnc before and after experiencing the low state.
Our photometric power search analysis (see § 4.1) shows a clear spin-cycle modulation, in agreement with Patterson et al. (2004) and Rodríguez-Gil et al. (2004). Our photometry also exhibits a moderate signal at 108 min. Patterson et al. (2004) found evidence of a similar weak signal at 110.85(9) min; they left this finding as an unsolved problem, which as they put it, does not seem related to any other clocks in the binary. Furthermore, we have also detected a new 62 min period signal which, although weak, corresponds to the beat period of the spin and the 108 min signatures. Finding this beat signal suggests that the 108 min period is not caused by spurious effects. Still, the origin of these modulations remains unknown. However, we did not find the 70 min beat period signal, nor the 86 min orbital period after the subtraction of the main spin modulation; the lack of signals is the main contrast in our photometric results from those by Patterson et al. (2004). We find the same contrast with the modulations found in the X-ray observations of Nucita et al. (2019), who also detected the spin-orbit beat period and a spin modulation signal.
We also conducted two different power spectrum analyses of the radial velocities of the Hα and He I λ 5876 Å emission lines (see § 4.2): the first one by convolving the lines with two antisymmetric Gaussians (GAU2), and the second one by employing the derivative of a Gaussian (DGAU) as the convolution function. When employing the same method as Patterson et al. (2004), i.e. DGAU, we obtain signals consistent with these authors, detecting both the orbital and spin modulations. On the other hand, the power search in the GAU2 option (which traces the inner regions of the disc) yielded only a strong signal for the orbital period. This, in a way, is also similar to the results reported by Rodríguez-Gil et al. (2004), who find both signals when using a broad Gaussian but obtain exclusively the orbital signature when the correlation of the emission line is made with a narrow Gaussian.
In § 5, we conducted a wavelet transform analysis which confirmed the presence of some of the signals found in the Lomb-Scargle periodograms. However, and perhaps more remarkably, this analysis also detected a localized signature of the 70 min beat period in the DGAU scalograms. We note that a weak hint of the beat modulation also appeared in the Hα periodograms after subtracting the dominant signals from the spectroscopic data, but we regard this detection with caution, given the high FAP yielded by the signature (see § 4.2).
We implemented a Doppler tomography study of the binary in § 6. The trailed spectrum of Hα exhibits a double-peaked structure that becomes single-peaked in short intervals, and with a changing relative intensity of the peaks. On the other hand, the He I λ 5876 Å trailed spectrum shows a broad single-peaked profile. Both of these trailed spectra are consistent with the structure of the profiles of the same emission lines reported by Rodríguez-Gil et al. (2004). Moreover, our Hα Doppler tomogram dis- plays the presence of the accretion disc, a structure not detected in the tomograms by Segura Montero et al. (2020), hinting at a possible replenishing of the disc. Our tomography also shows what appears to be a hot spot component for both emission lines, in good agreement with the S-wave, presumably originated at the location of the bright spot, exhibited in the trailed spectra diagrams from Rodríguez-Gil et al. (2004).
The similarities with previous studies, listed above, lead us to believe that the system has under- gone at least a partial recovery from the low state. In particular, finding a clear indication of the spin period both in our photometry and in the DGAU option, a signal which avoided detection in Segura Montero et al. (2020), suggests that the outer disc has undergone enough replenishing to provoke a detectable lighthouse beacon. This is further supported by the tomography showing a clear indication of the accretion disc and a hot spot emission.
Nonetheless, the main differences regarding the signals found (and not found) in the Lomb-Scargle power search of the photometry prevent us from declaring complete recovery of the previous state of DW Cnc. It remains to be seen if the system will fully recover to the condition reported in Patterson et al. (2004). The results in the present article show some progress, but further observations are still required to see whether if the mechanisms causing the behaviour reported in 2004 require more time to get kick-started and allow detection.
8. CONCLUSIONS
Photometric and spectroscopic observations of DW Cnc show, to some extent, a behaviour that agrees with that exhibited before experiencing a low state; pointing at a partial recovery of the system. Namely, our photometry yields a strong modulation of the spin period, which eluded detection by Segura Montero et al. (2020), indicating the reactivation of the lighthouse effect; our analysis also showed a weak unresolved signal at 108 min. Furthermore, when implementing the same methodology as Patterson et al. (2004) to measure the radial velocities of the emission lines, we obtained a result consistent with these authors in the periodogram, where we detect both the orbital and spin periods. Furthermore, in agreement with Rodríguez-Gil et al. (2004), we find evidence of the disc structure and hot spot emission in our Doppler tomography study. Finally, our wavelet transform analysis displays a localized detection of the 70 min beat period.
However, the signatures in the photometric periodograms do not completely match those reported before the low state by Nucita et al. (2019) and Patterson et al. (2004), who find not only the spin cycle, but also a signal corresponding to the spin-orbit beat period, and even a weak detection of the orbital signal after further treatment of their data. Instead we found a new 62 min period which corresponds to the beat between the spin and 108 min periods.
We could not find public data that would enable us to replicate our analyses. Therefore, we propose additional observations of the system to assess if it is possible that the mechanisms that gave rise to the signatures exhibited before the low state, require more time to completely rekindle.