SciELO - Scientific Electronic Library Online

 
vol.58 número1Off-Axis Exact Ray Tracing Algorithm for Zero Coma Point Determination in Classical and Non-Classical Reflective TelescopesMorphological Study of a Sample of Dwarf Tidal Galaxies Using the C-A Plane índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados

Journal

Artigo

Indicadores

Links relacionados

  • Não possue artigos similaresSimilares em SciELO

Compartilhar


Revista mexicana de astronomía y astrofísica

versão impressa ISSN 0185-1101

Rev. mex. astron. astrofis vol.58 no.1 Ciudad de México Abr. 2022  Epub 23-Jan-2023

https://doi.org/10.22201/ia.01851101p.2022.58.01.04 

Articles

New Photometry and Spectroscopy of DW Cancri

S. H. Ramírez1 

O. Segura Montero1 

R. Michel2 

J. Echevarría1 

1Instituto de Astronomía, Universidad Nacional Autónoma de México, Ciudad Universitaria, Ciudad de México, México.

2Instituto de Astronomía, Universidad Nacional Autónoma de México, Ensenada, Baja California, México.


ABSTRACT

We present new observations of the cataclysmic variable DW Cancri, after the system recovered from a low state. We performed a power spectrum analysis that reveals a clear signal of the 38 min spin period in our photometric data. Our spectroscopic power spectrum search was consistent with studies performed before the low state, showing the orbital and spin modulations. Our Doppler Tomography study exhibits a disc structure and an enhanced asymmetric region, possibly related to a hot spot component. Furthermore, a wavelet transform analysis reveals the 70 min spin-orbit beat period. We interpret these results as an indication of a partial recovery of the system. However, DW Cnc does not yet show all the original photometric modulations reported before the low state. Thus, we propose further observations to elucidate if such original signals require more time to reactivate.

Key Words: binaries; spectroscopic; novae; cataclysmic variables; stars; individual; DW Cnc; techniques; photometric; techniques; spectroscopic

RESUMEN

Presentamos nuevas observaciones de la variable cataclísmica DW Cancri después de que se ha recobrado de un estado bajo. Un análisis de espectro de potencias revela clara señal del período de espín de 38 minutos en nuestra fotometría. Nuestro espectro de potencias espectroscópico concuerda con estudios anteriores al estado bajo, mostrando los períodos orbital y de espín. Nuestro estudio de Tomografía Doppler exhibe una estructura de disco y una asimetría intensa, posiblemente causada por una zona de impacto. Un análisis de ondículas revela la presencia localizada del período de batimiento de 70 minutos. Interpretamos estos resultados como una recuperación parcial del sistema. Sin embargo, DW Cnc no muestra todas las señales reportadas, previo al estado bajo, en la fotometría. Proponemos observaciones adicionales para dilucidar si dichas señales requieren mayor tiempo para reactivarse.

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.

TABLE 1 LOG OF PHOTOMETRIC OBSERVATIONS FOR DW CNC 

Date Julian Date
(2450000 +)
No. of exposures
V R
08 March 2020 8916 201 210
09 March 2020 8917 201 219
15 March 2020 8923 197 191
16 March 2020 8924 131 133

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.

TABLE 2 LOG OF SPECTROSCOPIC OBSERVATIONS FOR DW CNC 

Date Julian Date
(2450000 +)
No. of spectra
15 March 2020 8923 28
16 March 2020 8924 51

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:

V(t)=γ+K1sin(2πt-t0Porb), (1)

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 χν2 as our goodness-of-fit parameter. Note that we have fixed the orbital period, as derived in Segura Montero et al. (2020), and therefore we only fit the other three parameters.

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 χν2=1. The parameters for the optimal orbital fit of both emission lines are shown in Table 3, with the respective estimated standard errors for the best-fit values.

Fig. 1 Diagnostic diagram of the Hα emission line. The vertical blue dashed line indicates an optimal separation of 103.9 Å (62 pixels). The used width for the Gaussians was 16.8 Å (10 pixels). See text for further discussion. The colour figure can be viewed online. 

Fig. 2 Radial velocity curve for the best solution of the Hα emission line. The best fit is shown as the blue line, and the 1σ error bars have been scaled so that the goodness-of-fit parameter χν2=1. The colour figure can be viewed online. 

Fig. 3 Diagnostic diagram of the He I λ 5876 Å emission line. The vertical blue dashed line indicates an optimal separation of 78.7 Å (47 pixels). The used width for the Gaussians was 11.7 Å (7 pixels). See text for further discussion. The colour figure can be viewed online. 

Fig. 4 Radial velocity curve for the best solution of the He I emission line. The best fit is shown as the blue line, and the 1σ error bars have been scaled so that the goodness-of-fit parameter χν2=1. The colour figure can be viewed online. 

TABLE 3 ORBITAL PARAMETERS OBTAINED FROM Hα AND He I λ 58761 

Parameter 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.

Fig. 5 Power spectra of the photometric observations. From top to bottom the panels depict the power spectra of the V band, the R band and of the combined data. See text for further discussion. The colour figure can be viewed online. 

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.

Fig. 6 Power spectrum of the residuals, after removing the spin period signal from the photometric data. The colour figure can be viewed online. 

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.

Fig. 7 Photometric data folded by the period signals found in the power spectrum search. The data have been averaged into 60 phase bins. The HJD 0 was selected manually to begin each modulation at its highest value, except for the orbital period for which we used the HJD 0 found in Segura Montero et al. (2020). The colour figure can be viewed online. 

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.

Fig. 8 Top: Power spectra of Hα using GAU2. Bottom: Power spectra of the He I λ5876 Å emission line using GAU2. See text for further discussion. The colour figure can be viewed online. 

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.

Fig. 9 The upper panel shows the power spectrum of the residuals of the GAU2 Hα radial velocity data, after removing the orbital period signal. The lower panel shows the same for He I λ5876 Å. On each panel, the solid orange line represents the power search performed on the residuals. For comparison, we superposed the power spectrum of the original data, plotted as the faint solid blue line. The colour figure can be viewed online. 

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.

Fig. 10 Folded spectroscopic data. The top panel is folded by the orbital period. The bottom panel is folded by the spin period. The colour figure can be viewed online. 

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.

Fig. 11 Power spectra of the Hα (top) and He I λ 5876 Å (bottom) emission lines, using the DGAU option with a Gaussian width of 50 pixels (see text). The colour figure can be viewed online. 

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.

Fig. 12 The upper panel shows the power spectrum of the residuals of the DGAU Hα radial velocity data, after removing the spin period and orbital period signals. The lower panel shows the same for He I λ5876 Å, after subtracting the orbital period. On each panel, the solid orange line represents the power search performed on the residuals. For comparison, we superposed the power spectrum of the original data, plotted as the faint solid blue line. The colour figure can be viewed online. 

In Figure 13 we folded the data by the spin and orbital periods, where both signals show clear modulations of the DGAU radial velocities.

Fig. 13 Radial Velocity data of Hα and and He I λ 5876 Å, folded by the spin (top panel) and orbital (bottom panel) periods found in the power spectrum search. The colour figure can be viewed online. 

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.

Fig. 14 Scalogram of the V -band photometric data. The global spectrum is exhibited in the utmost right panel, while the local spectrum appears in the panels below the data of each night. An orbital cycle signal (≈16 cycles/day) is present throughout all 4 nights of observations. See text for further discussion. The colour figure can be viewed online. 

Fig. 15 Scalogram of the R-band photometric data. The global spectrum is exhibited in the utmost right panel, while the local spectrum appears in the panels below the data of each night. The second night shows the presence of the spin period and a ≈21 cycles/day signal (consistent with the spin-orbit beat period). The rest of the nights are dominated by the orbital cycle. See text for further discussion. The colour figure can be viewed online. 

Fig. 16 Scalogram of the Hα gau2 spectroscopic data. The global spectrum is exhibited in the utmost right panel, while the local spectrum appears in the panels below the data of each night. The first night shows prominent signals at ≈37 cycles/day and ≈22 cycles/day. The latter signal also appears during the second night but is surpassed in intensity by a ≈33 cycles/day signature. See text for further discussion. The colour figure can be viewed online. 

Fig. 17 Scalogram of the He I λ5876Å GAU2 spectroscopic data. The global spectrum is exhibited in the utmost right panel, while the local spectrum appears in the panels below the data of each night. The first night shows a dominant signal at ≈22 cycles/day and a second one at ≈40 cycles/day. The second night is dominated by a ≈33 cycles/day signature, accompanied by a persistent secondary signal at ≈16 cycles/day. See text for further discussion. The colour figure can be viewed online. 

Fig. 18 Scalogram of the Hα DGAU spectroscopic data. The global spectrum is exhibited in the utmost right panel, while the local spectrum appears in the panels below the data of each night. The first night shows a strong signal, consistent with the spin cycle, at ≈38 cycles/day. During the second night the spin cycle signal is also present although with a broader structure; a secondary signature also appears this night at ≈20 cycles/day, which is consistent with the spin-orbit beat period. See text for further discussion. The colour figure can be viewed online. 

Fig. 19 Scalogram of the He I λ5876Å DGAU spectroscopic data. The global spectrum is exhibited in the utmost right panel, while the local spectrum appears in the panels below the data of each night. The first night shows a strong and persistent detection at the ≈20 cycles/day beat signal. Another signature is evident at the ≈38 cycles/day spin signal. The latter is also present during the second night, although its structure considerably broadens. The second night shows also the appearance of a persistent signal of the orbital cycle at ≈16 cycles/day. See text for further discussion. The colour figure can be viewed online. 

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.

Fig. 20 Trail spectra and Doppler tomography of the Hα emission line. The relative emission intensity is shown in a scale of colours, where the strongest intensity is represented by black, followed by red, then blue, and finally yellow. The cross marks represent (from top to bottom) the position of the secondary, the centre of mass and the primary component. The Roche lobe of the secondary is depicted around its cross. The Keplerian and ballistic trajectories of the gas stream are marked as the upper and lower curves, respectively. The colour figure can be viewed online. 

Fig. 21 Trail spectra and Doppler tomography of the He I λ 5876 Å emission line. The relative emission intensity is shown in a scale of colours, where the strongest intensity is represented by black, followed by red, then blue, and finally yellow. The cross marks represent (from top to bottom) the position of the secondary, the centre of mass and the primary component. The Roche lobe of the secondary is depicted around its cross. The Keplerian and ballistic trajectories of the gas stream are marked as the upper and lower curves, respectively. The colour figure can be viewed online. 

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.

Acknowledgements

The authors are indebted to DGAPA (Universidad Nacional Autónoma de México) support, PAPIIT projects IN114917 and IN103120. We thank the staff at the OAN-SPM for facilitating and helping us to obtain our observations. This research made extensive use of astropy, a community-developed core Python package for Astronomy (Astropy Collaboration et al. 2013); Python’s SciPy signal processing library (Virtanen et al. 2020) and matplotlib (Hunter 2007). We also thank the anonymous referee, whose comments substantially helped to improve the content of this paper.

REFERENCES

Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, 33, https://doi.org/10.1051/0004-6361/201322068 [ Links ]

Benítez, R., Bolós, V. J., & Ramírez, M. E. 2010, CMwA, 60, 634, https://doi.org/10.1016/j.camwa.2010.05.010 [ Links ]

Bravo, J. P., Roque, S., Estrela, R., Leão, I. C., & De Medeiros, J. R. 2014, A&A, 568, 34, https://doi.org/10.1051/0004-6361/201323032 [ Links ]

de Lira, S. R., Bravo, J. P., Leão, I. C., et al. 2019, MNRAS, 484, 3935, https://doi.org/10.1093/mnras/stz122 [ Links ]

Echevarria, J. 1983, RMxAA, 8, 109 [ Links ]

Echevarría, J., de la Fuente, E., & Costero, R. 2007, AJ, 134, 262, https://doi.org/10.1086/518562 [ Links ]

Figueira, P., Faria, J. P., Adibekyan, V. Zh., Oshagh, M., & Santos, N. C. 2016, OLEB, 46, 385, https://doi.org/10.1007/s11084-016-9490-5 [ Links ]

Hameury, J.-M. & Lasota, J.-P. 2017, A&A, 602, 102, https://doi.org/10.1051/0004-6361/201730760 [ Links ]

Hellier, C. 2001, Cataclysmic Variable Stars - How and Why They Vary (London, UK: Springer) [ Links ]

Hernández Santisteban, J. V. 2021, Astrophysics Source Code Library, ascl:2106.003 [ Links ]

Horne, K. & Marsh, T. R. 1986, MNRAS , 218, 761, https://doi.org/10.1093/mnras/218.4.761 [ Links ]

Hunter, J. D. 2007, CSE, 9, 90, https://doi.org/10.1109/MCSE.2007.55 [ Links ]

Knigge, C. 2006, MNRAS , 373, 484, https://doi.org/10.1111/j.1365-2966.2006.11096.x [ Links ]

Krischer, L., Megies, T., Barsch, R., et al. 2015, CS&D, 8, 014003, https://doi.org/10.1088/1749-4699/8/1/014003 [ Links ]

Marsh, T. R., Horne, K. 1988, MNRAS , 235, 269, https://doi.org/10.1093/mnras/235.1.269 [ Links ]

Newville, M., Stensitzki, T., Allen, D. B., & Ingargiola, A. 2014, zndo, 11813, https://doi.org/10.5281/zenodo.11813 [ Links ]

Nucita, A. A., Conversi, L., & Licchelli, D. 2019, MNRAS , 484, 3119, https://doi.org/10.1093/mnras/stz118 [ Links ]

Patterson, J., Thorstensen, J. R., Vanmunster, T., et al. 2004, PASP, 116, 516, https://doi.org/10.1086/421034 [ Links ]

Rodríguez-Gil, P., Gänsicke, B. T., Araujo-Betancor, S., & Casares, J. 2004, MNRAS , 349, 367, https://doi.org/10.1111/j.1365-2966.2004.07512.x [ Links ]

Scargle, J. D. 1982, ApJ, 263, 835, https://doi.org/10.1086/160554 [ Links ]

Schneider, D. P. & Young, P. 1980, ApJ, 238, 946, https://doi.org/10.1086/158059 [ Links ]

Segura Montero, O., Ramírez, S. H., & Echevarría, J. 2020, MNRAS , 494, 4110, https://doi.org/10.1093/mnras/staa856 [ Links ]

Shafter, A. W., Szkody, P., & Thorstensen, J. R. 1986, ApJ, 308, 765, https://doi.org/10.1086/164549 [ Links ]

Spruit, H. C. 1998, arXiv:astro-ph/9806141 [ Links ]

Stepanyan, D. A. 1982, PZ, 21, 691 [ Links ]

VanderPlas, J. T. 2018, ApJS, 236, 16, https://doi.org/10.3847/1538-4365/aab766 [ Links ]

Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, NatMe, 17, 261, https://doi.org/10.1038/s41592-019-0686-2 [ Links ]

Warner, B. 1995, Cataclysmic, Variable Stars (Cambridge, UK: CUP) [ Links ]

Warner, B. & Nather, R. E. 1971, MNRAS , 152, 219, https://doi.org/10.1093/mnras/152.2.219 [ Links ]

1IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc. (AURA), under cooperative agreement with the National Science Foundation (NSF).

Received: June 30, 2021; Accepted: November 17, 2021

J. Echevarría, S. H. Ramírez, and O. Segura Montero: Instituto de Astronomía, Universidad Nacional Autónoma de México, Apartado Postal 70-264, Ciudad Universitaria, Ciudad de México, C. P. 04510, México.

R. Michel: Instituto de Astronomía, Universidad Nacional Autónoma de México, Apartado Postal 877, Ensenada, Baja California, C. P. 22830, México.

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License