1.Introduction
Cataclysmic variables stars (CVs) are semi-detached binary systems composed by a white dwarf (WD) as primary star and a late-type main sequence secondary star. The later fills its Roche Lobe and transfers matter through the Lagrangian point L1 onto the WD via Roche-lobe overflow.
In non-magnetic systems, the incoming material forms an accretion disk around the WD (Warner 1995). This accretion disk periodically reaches a critical surface density value, leading to the development of thermal instabilities (Meyer & Meyer-Hofmeister 1981), triggering an outburst that increases the system brightness by up to 8 magnitudes. The system returns to quiescence in a few days or weeks. This photometric behavior defines a CV class named Dwarf Nova (Patterson et al. 1981; Howell et al. 1995; Szkody & Mattei 1984). Besides the light curve pattern method to classify CVs, there are other parameters that help to do so, such as the orbital period, the mass transfer rate or the magnetic field strength.
The orbital period defines the separation of the binary components and establishes the evolutionary status of the secondary star. CVs with periods over 3 hours contain secondaries stars with radiative cores, and as the orbital period decreases below 3 hours the companion reaches a mass low enough to become fully convective (Verbunt & Zwaan 1981; Rappaport, Verbunt, & Joss 1983; Paczynski & Sienkiewicz 1983; Spruit & Ritter 1983).
A high mass transfer rate defines the nova-like (NL) CV class. It helps to maintain a hot and stable accretion disk and prevents dwarf-nova-type outbursts. We refer to Warner (1995) for a general review on NLs and CVs in general, or Mizusawa et al. (2010); Balman et al. (2014) where they describe the UV and X-ray properties of several NL systems.
The strength of the magnetic field defines two types of CVs named Polars and
Intermediate Polars (IP). In Polars, the magnetic field is so strong (
Within the CV zoo it is possible to find the SW Sex stars, which are a subclass of high accretion rate CVs; i.e. NL members. They were first identified by Thorstensen et al. (1991) as eclipsing NL CVs with an orbital period around 3 to 4 hours. Radial velocities studies reveal that the emission lines vary periodically; the Balmer series lines lagging behind the expected phase for a WD according to photometric ephemeris. Their spectra show single or double-peaked emission lines, regardless of inclination angle. The single-peaked emission lines are thought to originate from the material encountering a magnetic accretion curtain close to the surface of the WD (Hoard et al. 2003). In contrast, the double-peaked profiles are likely a consequence of phase-dependent absorption components, as revealed by their trailed spectra.
Nowadays, SW Sex stars are thought to be the dominant population of CVs with periods within 3 and 4.5 hours (Rodríguez-Gil et al. 2007). A few cases of confirmed SW Sex stars show circular polarization, reflecting its magnetic nature (Rodríguez-Gil et al. 2001). Additionally, photometric quasi-periodic oscillations (QPOs) are features observed recently in SW Sex stars. It is possible that this QPOs reflect the rotation of an underlying magnetic WD (Patterson et al. 2002). Finally, evidence has been mounting that the SW Sex phenomenon is an evolutionary stage in the evolutionary process of CVs (Schmidtobreick et al. 2012).
SDSS J1224 (12:24:05.58+18:41:02.7) is a relatively bright star with ց= 16.01 mag. The presence of absorption and emission lines in its spectrum point to its being a pre-CV; a possible source for this could be the irradiation from the secondary by a hot WD (Szkody et al. 2011). Szkody et al. (2014) reported 126 min of time-resolved spectra taken in 2011, and concluded that there was no significant radial velocity variation during that time interval, so the system likely has a long orbital period or a low inclination angle.
In this work we present observational evidence that allows to classify SDSS J1224. This paper is organized as follows. In § 2 we present our observations, and in § 3 the data analysis and results. In § 4 a discussion is provided with our main conclusions stated in § 5.
2.Observations and Data Reduction
2.1.Photometric Observations
Differential time-resolved photometry of SDSSJ 1224 was conducted using the direct CCD imaging mode of a 0.84m telescope, located at Observatorio Astronómico Nacional at the Sierra San Pedro Mátir (OAN SPM1) in México. We acquired a long series of photometric data in the V broadband Johnson-Cousins filter, with exposure times ranging from 30 to 60 s. We observed the system on three nights in April 2016, four nights in March 2017 and three nights in April-May 2017 with the same instrument and configuration. Table 1 shows the log of these photometric observations.
Date | HJD Start +245700 | Exp. Time Number of Integrations | Duration h |
01 April 2016 | 480.73269 | 30 s × 387 | 3.2 |
02 April 2016 | 481.64345 | 60 s × 387 | 6.5 |
03 April 2016 | 482.67408 | 60 s × 240 | 4.0 |
23 March 2017 | 836.729960 | 30 s × 200 | 1.7 |
24 March 2017 | 837.682042 | 60 s × 370 | 6.2 |
25 March 2017 | 838.667957 | 60 s × 240 | 4.0 |
26 March 2017 | 839.645395 | 60 s × 210 | 3.5 |
29 April 2017 | 872.667346 | 60 s × 281 | 4.7 |
30 April 2017 | 873.660104 | 60 s × 270 | 4.5 |
01 May 2017 | 874.650873 | 60 s × 210 | 4.0 |
Data reduction was performed using standard IRAF procedures (Tody 1986, 1993). The images were bias-corrected and flat-fielded before aperture photometry was carried out. We used photometric aperture radii of 2.0 times the PSF FWHM. The uncertainty in the differential photometry was estimated to be in the range 0.02 to 0.05 mag, according to the magnitude dispersion observed for the field stars. We used the star TYC1445-830-1 in the field of view as the reference star.
2.2.Spectroscopic Observations
We carried out spectroscopic observations using the 2.12m telescope located at OAN SPM, with the low to intermediate resolution Boller & Chivens (B&Ch) spectrograph. We acquired spectra with a resolution of 5.5 Å using a 600 lines/mm grating to cover a 4300 - 5700 Å range. The observations were made through a 1.5´´ slit, oriented in the east-west direction. CuHeNeAr lamp exposures were taken every 60 min for wavelength calibration and for flux calibration; spectrophotometric standards from Oke (1990) catalog were observed. The exposure time was 900 s per spectrum. The image processing was carried out with standard IRAF procedures (Tody 1986, 1993).
3.Data Analysis and Results
3.1.Light Curve Morphology
We observed CV SDSSJ122405.58+184102.7 over ten nights, spread over two years, covering 42 hr of total photometry. The system was at quiescence during our observations, as we can see in the top panel of Figure 1. However, individual light curves (individual 10 frames in Figure 1 have two conspicuous brightness variations, occurring at two different time scales but with almost the same amplitude: a long photometric signal with ≈ 0.2 mag amplitude and a short one with ≈ 0.1 mag amplitude.
On HJD = 836 (frame N4 on Figure 1) the light curve exhibits a sinusoidal pattern with a mean amplitude of ≈ 0.1 mag and a period of tens of minutes. We can appreciate the later in more detail in the upper panel of Figure 2. However, this modulation nearly disappears when the observation time is longer, like in HJD = 837 (frame N5 in Figure 1). In the bottom panel of Figure 2 the long term photometric signal is highlighted. In that case, the brightness changes mainly because of the orbital motion. Even masked inside the long photometrical signal, the 0.1 mag modulation is present on other nights (darker circles in Figure 1).
3.2.Photometric Orbital Period
In Figure 1, all observations are displayed in different frames. We applied a discrete Fourier transform (DFT) algorithm, using the software Period04 (Lenz & Breger 2005), to our complete photometric data set to search for significant periodic signals. The results obtained in this way are compatible with a Lomb-Scargle periodogram.
Panel (a) of Figure 3 shows a fraction of the total power spectrum for a ten days light curve after subtracting the spectral window. We explored the frequency range from 0 to 600 cycles day-1 (the Nyquist limit). However, the power spectrum is dominated by short frequencies. The maximum frequency is located at Ω = 5.959119 cycles day-1 with a SNR of 7.98. According to empirical results of Breger et al. (1993), and numerical simulations from Kuschnig et al. (1997), a power spectrum with a SNR larger than 4.0 is needed to ensure that the signal is a real feature. We carried out a Monte Carlo simulation to improve the orbital frequency and to calculate its uncertainty. This kind of Monte Carlo simulation is described in Mennickent & Tappert (2001), Mennickent et al. (2002) and Aviles et al. (2018). We generated 1000 data sets with the same HJD as the original time string and the magnitudes were those obtained by the last fit plus Gaussian noise. For every new data set, a least-squares fit was computed. The frequency distribution of this simulation was used to improve the orbital frequency and the uncertainty corresponds to the standard deviation of this distribution. With this procedure, the frequency was improved to 5.959076 cyclesday-1 which corresponds to a period P = 1/𝛺 = 0.167811(1) days (= 4.027464 (3) h) adopted here as the orbital one. In panel (b) of Figure 3 we present a Monte Carlo histogram that shows the frequency improvement compared with that obtained from the FFT power spectrum (inset). Finally, in panel (c) of Figure 3 we show the folded light curve with the period P =4.03 h beside a sinusoidal fit to remark the 0.2 magnitude long term modulation.
3.3.Short Time-Scale Variability in the Light Curve
The power spectrum is the most widely used tool to detect any quasi-periodic oscillation. In such a case, the power spectrum will show a broad peak that covers many frequencies, rather than a narrow peak centered on a given frequency. In other words, these oscillations are not coherent over time. They can change in a time scale of less than a day, like in V442 Oph and RX J1643.7+3402 (Patterson et al. 2002), or can be coherent over 20 cycles like in HS 0728+6738 (Rodríguez-Gil et al. 2004).
To inquire if a short time-scale oscillation is present in other nights, and to gain insight about the HDJ = 836 light curve behavior, we used Period 04 to analyze the data marked with darker circles in Figure 1 (Frames N1, N3, N4 and N10). We selected at least one night from each observational campaign. In some cases it was necessary to detrend the data to be able to carry out the period analysis, as in Kennedy et al. (2016), where they masked the light curve to remove the eclipses in MASTER OTJ192328.22+612413.5. We explored a wide range of frequencies as was done for the orbital period analysis. In Panel (a) of Figure 4 we show a portion of the power spectrum, and the frequency that modulates this data with a high probability ω = 50.423395 day-1, equivalent to a period of P = 28.558 min. We also carried out a Monte Carlo simulation to improve this value. The histogram in Panel (b) of Figure 4 shows the most favored frequency for this variability, ω = 50.423990 days-1. In Panel (c) we present the data folded with a period of 28.6 min, where a wave pattern is observed with an amplitude of ≈ 0.1 magnitudes.
We assume that this short-periodic variation may be related to the spin period of the WD. Therefore, SDSSJ1224 could be classified as an IP. Another evidence in favor of an IP nature is the detection of beat frequencies ω - 2Ω, ω - Ω, ω + Ω, ω + 2 Ω, where ω and Ω are spin and orbital frequencies, respectively. The existence of the spin-orbital sidebands for IPs was pointed out by Warner (1986) and Wynn & King (1992). The optical beat period is thought to arise from the reprocessing of soft X-rays by parts of the system fixed on the binary frame, like the secondary star itself.
3.4.Spectroscopic Features of SDSSJ1224
The average of 36 spectra of SDSSJ1224, acquired on March 19, 2017 (HJD ≈ 834.66), is presented in Figure 5 This spectrum shows single-peaked line profiles dominated by hydrogen (Balmer series) and HeI (λ5015, λ4922 and λ4388) emission lines. High excitation lines like HeII λ4886, the CIII / NIII λ4645 Bowen blend and OII 𝜆4416 are also present. Of particular interest is the presence of the HeII λ4886 line, because it is recognized as evidence of the presence of a source of ionizing photons, typical of magnetic CVs or NL variables (Araujo-Betancor et al. 2003; Rodriguez-Gil et al. 2004). Table 2 provides the equivalent widths (EW) and the full-width half-maximum (FWHM) of the main emission lines detected in the averaged spectrum. The FWHMs were obtained by fitting a single Gaussian to the line profiles. Taking into account the small FWHM of the lines, the inclination angle is expected to be low for SDSSJ1224.
Line | FWHM[Å] | EW[Å] |
H λ | 6.425 | 1.312 |
He IIλ4886 | 9.512 | 0.506 |
Hβ | 7.177 | 2.035 |
He Iλ4922 | 7.032 | 0.234 |
He Iλ5015 | 8.850 | 0.256 |
SDSSJ1224 has an optical spectrum in quiescence similar to RX J1643.7+3402, HS 0728+6738 and HS 1813+6122, as reported by Patterson et al. (2002), Rodríguez-Gil et al. (2004), Rodríguez-Gil et al. (2007), respectively. In all cases, they argue the SW Sex membership according to spectroscopic properties and photometric variabilities. Rodríguez-Gil et al. (2009) detected circular polarization in RX J1643.7+3402, confirming in this way the magnetic nature of SW Sex stars.
In Figure 5 we observe strong evidence for the detection of the WD as absorption in the blue wings of both the H𝛽 and H𝛾 lines. This is in agreement with the results reported by Szkody et al. (2011, 2014). Our photometry and FWHM measurements support the idea of a long orbital period, as well as a low inclination angle for SDSSJ1224, instead of a low accretion rate as the reason for the WD detection.
4.Discussion
CV SDSSJ1224 is a non-eclipsing binary system that has certain photometric features to potentially lead to a classification as a SW Sex star. First, we found an orbital period of 4.03 h. According to Rodríguez-Gil et al. (2007), the orbital periods for such systems are from 3 to 4.5 h. Almost 50% of the CV population in this range has been classified as SW Sex stars. This sample includes eclipsing and non-eclipsing systems.
Second, we detected quasi-periodic oscillations with a period of 28.6 minutes. In terms of short period variabilities (kilo-seconds), SW Sex stars are also characterized by exhibiting quasi-periodic modulations in their light curves. In a compilation of CVs that show fast oscillations, Warner (2004) presented nine SW Sex stars that show quasi-periodic oscillations (QPOs) with a predominant time scale of ≈ 1000 - 2000s (16.7 - 33.3 min).
Third, another fact that supports our conjecture is that assuming an accretion disk-magnetic field interaction model, as proposed by Rodríguez-Gil et al. (2001), the spin period and the orbital one are related by
The previous equation results from considering that the shock between the gas stream and the accretion disk occurs close to the co-rotation radius. Here f represents the co-rotation radius expressed in units of RL1 (the distance between the WD and the inner Lagrangian point, 𝐿 1 ). In the case of SDSSJ1224, the value for f is 0.53, which is in good agreement with values reported for the prototype SW Sex itself (Groot 1999); LS Peg (Rodriguez-Gil 2001); V533 Her (Rodríguez-Gil & Martínez-Pais 2002). As a consequence of this model, SW Sex stars are indeed IP with the highest mass accretion rates.
Fourth, variabilities on time scales of minutes to tens of minutes have been detected in the optical light curve of several IPs (e.g., Patterson et al. 2002), and they are recognized as the spin period of the magnetic WD. We consider that, from its spectral appearance during low state (Figure 5) and the optical variability detected, SDSSJ1224 is indeed a new IP SW Sex class member. We suggest that the detected frequency of 50.423990 day-1 is the WD spin frequency ω(Pspin =1/ω)The quasi-sinusoidal shape of its light curve, which is typical for the spin light curve of IPs (see e.g., the spin light curve of FO Aqr by de Martino et al. 1994), lends support to our hypothesis that this frequency is indeed the spin frequency of the WD.
Fifth, magnetic WDs in spin equilibrium should closely satisfy the ratio Pspin / Porb = 0.1 (King & Lasota 1991), assuming that the accretion process is not primarily through a disk. The QPOs may reflect an underlying rotation of a magnetic WD (Patterson et al. 2002). In the case of SDSSJ1224 this assumption seems quite valid, since if the 28.6 min signal is the spin period of the WD, then Pspin = 0.12 \ Porb it , which is very close to the value around which IPs seem to cluster on hard X-ray surveys (Scaringi et al. 2010).
5.Final Comments and Conclusions
We conclude that, based on the overall light curve features, the reported photometric periods and the spectroscopic properties, the system SDSSJ1224 can be classified as a new non-eclipsing SW Sex. The observed behavior matches all the conditions to be a member of this class of CVs.
It has a photometric orbital period of 4.03 h, which lies in the range of periods where the vast majority of SW Sex stars are grouped. It shows photometric features common to this type of stars and, as such, the number of non-eclipsing systems of this particular type in this range of periods is increased.
We detect a brightness variation amplitude of 0.1 mag and a period of 28.6 min, that we associate to the spin period of the WD. Also, we detect beat frequencies between the WD spin frequency, ω, and the orbital frequency, Ω, of the order of 30 minutes (Porb = 8 Pbeat) like in the system RX J1643.7+3402 (Rodríguez-Gil et.al. 2009).
The photometric behavior fits very well with those observed on the out-of-eclipse light curve of several CVs like HS 0728+6738 (Rodríguez-Gil et al. 2004) or MASTER OTJ192328.22+612413.5 (Kennedy et al. 2016). In particular, these two objects show a QPO amplitude of 0.2 mag with a period around 20 minutes, quite similar to SDSSJ1224. This result may imply that the QPOs are not exclusive to eclipsing systems, as was pointed out by Patterson et al. (2001).
Because quasi-periodic oscillations in the range of kilo-seconds have recently been detected in confirmed SW Sex stars (Patterson 2002; Kennedy et al. 2017; Rodríguez-Gil et al. 2007), we consider that SW Sex stars are actually members of IP CVs stars.
The physical conditions present in this type of binary systems are not well understood and should be addressed both theoretically and observationally. However, such a study is beyond the scope of this work.
This work is based upon observations carried out at the Observatorio Astronómico Nacional at the Sierra San Pedro Mártir (OAN-SPM), Baja California, México. We thank the daytime and night support staff at the OAN-SPM for facilitating and helping us to obtain the observations. The authors also acknowledge PAICyT-UANL grant CE642-18 for resources provided toward this research. We are also thankful to the anonymous referee for a careful reading of the manuscript and for useful comments that improved the content of this work.