Servicios
Descargas
Buscar
Idiomas
P. Completa
Extensive photometry of V1838 Aql during the 2013 superoutburst
J. Echevarría; E. de Miguel; J. V. Hernández Santisteban;
J. Echevarría; E. de Miguel; J. V. Hernández Santisteban; R. Michel; R. Costero; L. J. Sánchez; A. Ruelas-Mayorga; J. Olivares; D. González-Buitrago; J. L. Jones; A. Oskanen; W. Goff; J. Ulowetz; G. Bolt; R. Sabo; F.-J Hambsch; D. Slauson; W. Stein
Extensive photometry of V1838 Aql during the 2013 superoutburst
Revista mexicana de astronomía y astrofísica, vol. 55, no. 1, pp. 21-30, April , 2019
Instituto de Astronomía, UNAM
resúmenes
secciones
referencias
imágenes

ABSTRACT: We present an in-depth photometric study of the 2013 superoutburst of the recently discovered cataclysmic variable V1838 Aql and subsequent photometry near its quiescent state. A careful examination of the development of the superhumps is presented. Our best determination of the orbital period is PORB = 0.05698(9) days, based on the periodicity of early superhumps. Comparing the superhump periods at stages A and B with the early superhump value we derive a period excess of = 0.024(2) and a mass ratio of q = 0.10(1). We suggest that V1838 Aql is approaching the orbital period minimum and thus has a low-mass star as a donor instead of a substellar object.

Key Words: techniques, photometric, stars, novae, cataclysmic variables, stars, dwarf novae, stars, individual, V1838 Aql.

RESUMEN: Presentamos un estudio fotométrico detallado de la super-erupción de V1838 Aql, una variable cataclísmica recientemente descubierta, desde el máximo en 2013 hasta su regreso al mínimo. Examinamos en detalle la evolución de las superjorobas (superhumps). Determinamos el período orbital PORB = 0.05698(9) días a partir de la periodicidad de las superjorobas tempranas. Comparando los períodos de las superjorobas en las etapas A y B con el valor del período orbital, derivamos un valor del cambio en el período orbital de = 0.024(2) y un cociente de masa para el sistema de q = 0.10(1). Sugerimos que V1838 Aql se está acercando al mínimo período orbital, por lo que la secundaria sería una estrella de baja masa y no un objeto sub-estelar.

Carátula del artículo

Articles

Extensive photometry of V1838 Aql during the 2013 superoutburst

J. Echevarría
Universidad Nacional Autónoma de México, Mexico
E. de Miguel
Universidad de Huelva, Spain
J. V. Hernández Santisteban
University of Amsterdam, Netherlands
R. Michel
Universidad Nacional Autónoma de México, Mexico
R. Costero
Universidad Nacional Autónoma de México, Mexico
L. J. Sánchez
Universidad Nacional Autónoma de México, Mexico
A. Ruelas-Mayorga
Universidad Nacional Autónoma de México, Mexico
J. Olivares
Université of Bordeaux 1, France
D. González-Buitrago
University of California, USA
J. L. Jones
CBA-Oregon, USA
A. Oskanen
CBA-Finland, Finland
W. Goff
CBA, USA
J. Ulowetz
CBA, USA
G. Bolt
CBA, Australia
R. Sabo
CBA, USA
F.-J Hambsch
CBA, Belgium
D. Slauson
CBA, USA
W. Stein
CBA, USA
Revista mexicana de astronomía y astrofísica, vol. 55, no. 1, pp. 21-30, 2019
Instituto de Astronomía, UNAM

Received: 08 August 2018

Accepted: 24 October 2018

1. Introduction

Cataclysmic variables (CVs) are close binary systems in which a white dwarf (WD) accretes from a low-mass star via Rochelobe overflow, often creating an accretion disc (for a review see Warner 1995). A large fraction of CVs belong to the subclass of dwarf novae (DNe). They undergo recurrent outbursts with typical amplitudes of   2   -   6 mag in the optical, which are commonly accepted to be caused by a thermal-viscous instability in the disc (Osaki 1974). In addition, SU UMa-type DNe (the subclass of DNe systems with short orbital periods, Porb < 2.5 hr) exhibit occasional eruptions that are less frequent, longer lasting, and slightly brighter (by   0.5   -   1.0 mag) than the normal outbursts. The key feature during these so-called superoutbursts is the presence of superhumps a modulation in the light curve with an underlying periodicity, Psh, a few percent longer than the orbital period. They are thought to arise from a precessing non-axisymmetric disc (Vogt 1982), with the eccentricity being produced by the tidal instability developed at the radius of the 3:1 resonance (Whitehurst 1988). The analysis of the timing and evolution of such light oscillations provides estimates of the system’s parameters using empirical (Patterson et al. 2005) and theoretical (Kato & Osaki 2013) relationships between the superhump/orbital period excess, ϵ (Psh − Porb)/Porb, and the binary mass ratio qM2/M1. Several ϵ (q) relations have been proposed by these authors based on different stages of the superhumps, although Otulakowska-Hypka et al. (2016) show that the scatter in the ϵ - q diagram is considerable for very short orbital periods. This is due to a lack of objects with dynamically confirmed small values of q (Patterson 2011; Kato & Osaki 2013).

Among the SU UMa-type systems, there is a large group that accumulates around the minimum of the orbital period distribution of CVs (Porb ≈ 78 min) (Paczynski & Sienkiewicz 1981; Gänsicke et al. 2009; Knigge et al. 2011). These are binaries with extremely low mass-transfer rates, named WZ Sge-type objects, which are characterised by rare (commonly detected every ≈ 10 years), and large amplitude superoutbursts of duration of ≈ 30 days, caused by an instability in low viscosity accretion discs with α ≈ 0.01 − 0.001 (Smak 1993; Osaki 1994). Some of them are systems currently evolving towards longer periods and are collectively known as period bouncers (e.g. Littlefair, Dhillon & Martin 2003). These binaries are expected to harbour a substellar secondary companion (Howell et al. 1997), i.e. brown dwarfs (e.g. Littlefair et al. 2006; Harrison 2016; Hernández Santisteban et al. 2016; Neustroev et al. 2017).

Worth noting is that WZ Sge-type objects are characterised not only by a long superoutburst recurrence time in comparison to typical SU UMa stars, but also by the presence of early superhumps (double-wave modulation) during the first few days of the eruption, with a periodicity (Pesh) essentially equal to the orbital period of the binary (details in O’Donoghue et al. 1991; Kato 2015). The advent of all-sky surveys (e.g. Breedt et al. 2014) and worldwide citizen-telescope networks has contributed to the discovery of a large population of faint DNe. Among these discoveries, the elusive population of short-period systems, in particular period bouncers, has been found and investigated (Patterson 2011; Coppejans et al. 2016; Otulakowska-Hypka et al. 2016).

The discovery of a new transient, initially proposed as possible nova, was reported by Itagaki on 2013 May 31. Henden16 pointed out that the colour indices of the object and the un-reddened field suggested a DN rather than a nova. As pointed out by Hurst17, Kojima reported a pre-discovery image on 2013 May 30.721 UT, when the magnitude was at about 9.8 mag (un-filtered). A CBET report of this new DNe in Aquila, can be found in Itagaki (2013).

Although a preliminary analysis of the behaviour of V1838 Aql (originally designated as PNV J19150199+0719471) was published by Kato et al. (2014), we present here a full analysis of the superhump behaviour based on our extensive data.

In § 2 we present the observations and their reduction methods. The photometric data and the period analysis are presented in § 3, while in § 4 we address the discussion of our results. We present our conclusions in § 5.

2. Observations and reduction

Photometric observations in the V band were obtained in 2013 during the nights of June 3, 4, 5, 6, 17, 18, 28 and September 2 and 25 at the 0.84 m telescope of the Observatorio Astronómico Nacional at San Pedro Mártir (SPM). We used the Blue-ESOPO CCD detector18 on a 2 × 2 binning configuration. The exposure time of the SPM observations varied between 10 and 30 s. In addition, time-series photometry of the superoutburst was obtained from 10 observatories of the Center for Backyard Astrophysics (CBA) -a network of small (0.2 − 0.4 m) telescopes that covers a wide range in terrestrial longitude. Skillman & Patterson (1993) and de Miguel et al. (2016) describe the methods and observing stations of the CBA network. These observations amounted to 162 separate time-series during 58 nights from June 1 to August 2, 2013, and the typical exposure time ranged from 20 to 120 seconds, depending on the brightness of V1838 Aql. Nearly half of these observations were obtained in V light, while the rest (mainly during the post-outburst regime) were unfiltered. We did not attempt any absolute calibration of the data during the eruption, but the magnitude scale is expected to resemble closely V magnitudes with a zero-point uncertainty of   0.05 mag. Further observations in the R band using the 2.1 m telescope at SPM during 2018, July 18 were conducted. Unfortunately, the weather was unstable and we only managed to obtain differential photometry over three orbital cycles. In the following, we report times and refer to specific dates in a truncated form defined as HJD − 2, 456, 000.

3. Photometry and period analysis
3.1. Photometric Observations

Most of our photometric observations come from the CBA network, with additional V -band observations obtained with the 0.84 m telescope at SPM at some critical stages of the outburst and during the late decline.19 We observed the typical pattern seen in SU UMa stars during superoutburst: a plateau phase -lasting ≈ 25 days, from HJD 444 to 469- where the mean brightness varied smoothly from 10.5 to 13.0 mag, followed by a rapid decline (≈ 3 mag in 2 days) at the end of the main eruption. The subsequent fading towards quiescence occurred at a rather low rate (≈ 0.035 mag d-1), and even 3.5 years after the end of the eruption, the system was found to be ≈ 0.5 mag above the pre-eruption quiescent brightness. However, based on the observations obtained with the NTT telescope (La Palma, Spain), we confirmed that the object had reached the pre-outburst level by June 2017. These observations and the general spectral distribution at quiescence have no further relevance here and will be discussed in a future publication.

3.2. Full Analysis of the Different Stages of the Superhump

Our primary tool for studying periodic signals was the Period 04 package (Lenz & Breger 2005). First, we subtracted the mean and (linear) trend from each individual light curve and formed nightly-spliced light curves. Then, after combining light curves from adjacent nights, a search for periodic signals was done. This approach allows us to improve the frequency resolution, but it has to be implemented with caution, since variations in the amplitude and/or the period of the modulation -both effects known to affict erupting DNe- can distort the outcome of the frequency analysis.

A general view of the superhump transitions can be looked up by identifying the different stages of the superhumps: early superhumps, Stages A, B and C as well as the post-outburst stage (see Kato et al. 2009, 2014, for this terminology in our general discussion). Thus, we looked at the time variations of the superhump period and its amplitude by examining the variation in time with respect to a well-defined feature of the superhump signal.

After this general analysis, a detailed explanation of these stages was made. First, we derived the timings of superhump maxima. A total of 310 times of superhump maxima was identified in the light curves in the interval HJD 449.6-498.6 days. These maxima are shown in Table 1. The early stage of the eruption was not considered, since the signal there was of very low amplitude and individual maxima were not well defined. A linear regression to these timings provides the following test ephemeris:

T m a x H J D = 2,456,452.8035 23 + 0.058191 6   E . (1)

Table 1
Times of superhump maxima of V1838 Aql during the 2013 superoutburst.a

a Individual errors in these timings are not explicitly included here. b E (cycle number). c Superhump maxima expressed as HJD - 2,456,000. d O - C value (in cycles) according to the ephemeris. T_max (HJD) = 2,456,452.8035 + 0.058191 E.

The top panel of Figure 1 displays the general photometric behaviour of the system during our campaign. Next, on the middle panel of Figure 1 is presented the variation of the amplitude of the superhump modulation, defined as the semi-amplitude of the sine wave that best fits the nightly photometric data.


Fig. 1
Photometric behaviour of V1838 Aql during its superoutburst in 2013. Top frame: global light curve of the observational campaign. Middle frame: amplitude variations of the superhumps along the superoutburst (see text). Bottom frame: O-C diagram for the superhump maxima in the HJD 449.6-498.6 day interval with respect to the ephemeris given in equation (1). The arrows indicate the (approximate) location of transitions between different regimes of superhump period variations (see text).

The O-C residuals of the times of maximum light relative to the ephemeris given by equation 1 are shown in the lower panel of Figure 1. The resulting O-C diagram is complex, but it displays a number of features that are usually observed in other SU UMa-type systems (Kato 2015). Among the most relevant features visible in this diagram we point out the following:

  1. 1. During the first four days of the outburst a weak modulation (early superhumps) with a period Pesh ≈ 0.057 d was visible in the light curve.
  2. 2. The onset of fully-grown (Stage A superhumps) took place in a short time-scale (≈ 2 d) and involved an increase in the amplitude of the modulations. Their (mean) period, Psh(A) ≈ 0.059 d, was longer than Pesh.
  3. 3. Once the superhump modulation reached full amplitude, the system entered Stage B where the amplitude of the superhump decreased slowly, and the mean period became shorter (Psh(B) ≈ 0.058 d). The upward curvature of the residuals during this stage (days HJD 449-466) signifies that the period of the superhumps was not constant, but increased over time. From a quadratic fit of the residuals in this interval, we find an increase rate of d P s h ( B ) / d t = 5.8 ( 4 ) ×   10 - 5 .
  4. 4. Before the end of the main eruption, the amplitude of the superhumps was found to grow larger (≈ 0.10 mag). The system entered Stage C, extending from day HJD 466 to the end of the main plateau (around day HJD 470), where the period of the superhump remained essentially constant (Psh(C) ≈ 0.0582 d).
  5. 5. Worth noting is the increase of the amplitude variations as the system dropped by about 3 mag between Stage C and the post-outburst stage.
  6. 6. After the end of the main eruption (day ≥ HJD 473) the superhumps were still visible with a significantly larger amplitude ≈ 0.2 mag), and with a period which remained constant for at least the subsequent ≈ 25 days. The period of the post-outburst modulation was shorter than Psh.

A summary of the main periodicities along the eruption is given in Table 2, as found in the next subsections.

Table 2
Mean photometric periods and frequencies of the superhump modulation.*

* Values correspond to the 2013 superoutburst of V1838 Aql. The errors in parentheses correspond to the last significant figures.

3.2.1. Early Superhumps

From the beginning of our campaign, a weak modulation of about 0.010 mag full amplitude was observed in the light curves. This signal persisted over days HJD 445-448. The power spectrum of the spliced light curve covering this 4-day segment is shown in the upper frame of Figure 2. It is dominated by two broad peaks centred at frequencies 35.10(3) and 17.55(3) cycles d-1. These signals were weak, with amplitudes of 0.0036 and 0.0030 mag, respectively. Although they were barely detected above the noise, we interpret them as a likely manifestation of early superhumps (Kato et al. 2014).


Fig. 2
Upper frame: Power spectrum during days HJD 445-448 (early superhumps), showing broad peaks centred at 17.55 cycles d−1 and its first harmonic. Lower frame: Waveform of the early superhumps obtained after folding the data with Pesh = 0.05698 d. The zero phase is arbitrary.

This photometric feature is known to be typical of WZ Sge-type stars, and is not shown by any other type of dwarf nova. Although its physical origin is still under debate, there is increasing observational evidence that its period (Pesh) is essentially equal to Porb (Patterson et al. 1996; Kato 2015). A folded curve of the spliced light curve with Pesh = 0.05698(9) d is also shown in the lower frame of Figure 2, which shows the double-humped pattern characteristic of early superhumps. The value of Pesh obtained in this paper is slightly different from, but consistent with, the value of 0.05706(2) d reported in Kato et al. (2014).

3.2.2. Common Superhumps

The double-humped pattern of the early superhumps turned into single-peaked humps on day HJD 449. Over days HJD 449-451, the mean amplitude was around ≈ 0.007 mag and the period was ≈ 4% longer than the period found for early superhumps. The dominant signal occurred at 16.83(8) cycles d−1 corresponding to a period of 0.0594(3) d. Since the modulation was better defined in this 2-day interval, we were able to determine the times of maximum in the signal. We identified a total of 9 maxima in the HJD 449.6-451.6 day interval, and obtained a period of 0.05934(11) d (corresponding to a frequency of 16.85(3) cycles d−1) from a linear regression. This value is fully consistent with the one found from the Fourier analysis. The modulation over this 2-day segment is interpreted as Stage-A superhumps. The period we find is close to, but slightly different from, the value of 0.05883(6) d reported in Kato et al. (2014).

Fully-grown, large-amplitude superhumps were finally observed on day HJD 452 (amplitude of 0.10 mag). As a representative example, we show in the upper frame of Figure 3 the light curve from day HJD 455. As the eruption proceeded, the mean amplitude of the superhumps decreased (as shown in the middle panel of Figure 3). The variation in amplitude was smooth in the HJD 452-466 day interval. We formed a spliced light curve in this interval, and obtained the power spectrum shown in the middle frame of Figure 3. The strongest signals occurred at f1 = 17.128(3) and f2 = 34.265(3) cycles d−1.


Fig. 3
Upper frame: A 12-hour spliced light curve obtained on day HJD 455, dominated by large-amplitude common superhumps. The zero level in the figure corresponds to V ≈ 12.1 mag. Middle frame: Power spectrum during the common superhump era (days HJD 442- 466), with main peaks centred at frequencies 17.128 and 34.265 cycles d−1. Lower frame: Mean waveform of the common superhump obtained after folding the data on P = 0.058384 d. The zero phase is arbitrary.

They were interpreted as the frequency of Stage-B superhumps (period of Psh = 0.058384(10) d) and its first harmonic, respectively. Other (weaker) peaks, not shown in Figure 3, were found at f3 = 51.381(6) and f4 = 68.403(6) cycles d−1. The mean waveform of the superhump modulation during this interval is shown in the lower panel of Figure 3.

We note that after subtracting the superhump signal and its harmonics, the power spectrum of the residual light curve showed peaks at 17.20 and 17.28 cycles d−1. But we do not give any physical significance to these detections, and interpret them as the result of period and amplitude variations of the superhump wave during the eruption.

The amplitude of the superhumps increased around day HJD 466, and decreased thereafter until the end of the main eruption (Stage C). The strongest signal in the power spectrum in the HJD 466-470 day interval occurred at 17.191(10) cycles d−1, corresponding to a period of 0.05817(3) d, with additional peaks at higher harmonics.

3.2.3. Post-Outburst Stage

Once the main eruption was over, the light curve was still dominated by superhumps, but now with significantly larger amplitude (≈ 0.15 mag) which decreased slowly (≈ 0.013 mag d−1). This behaviour remained essentially unchanged for nearly 20 days of our observations after the main fading.

We formed a spliced light curve including all the observations from day HJD ≥ 473, and found a power spectrum (shown in the upper frame of Figure 4), with a peak at 17.244(1) cycles d−1. This was interpreted as the frequency of the post-outburst superhump, and dominated the spectrum. Higher-order harmonics were also found, but their amplitude was very low (< 0.0065 mag). This signifies that the waveform of the post-outburst superhump was nearly sinusoidal. The lower frame of Figure 4 shows that this was indeed the case.


Fig. 4
Upper frame: Power spectrum after the main eruption (day HJD ≥ 473) showing a strong peak centred at 17.244(1) cycles d−1 (post-outburst superhump). Lower frame: Mean waveform of the post-outburst superhump, obtained after folding the data on P = 0.057991 d. The zero phase is arbitrary.

3.3. Photometry Near and at Minimum Light

As detailed in § 2, we took two runs near minimum light covering around one orbital period each. The light curves are shown in the upper panel of Figure 5. When folded with the orbital period (0.05698 d), the light curves seem to be out of phase. But this is not surprising: over the 23 days (nearly 400 orbital cycles) elapsed between both runs, an uncertainty of 0.0001 days in Porb involves an uncertainty of 0.7 in phase. We carried out a period analysis of both nights using the Phase Dispersion Minimization (PDM) technique (Stellingwerf 1978) in the Peranso package (Paunzen & Vanmunster 2016). This technique is frequently used to detect variations of superhumps in SU UMa systems (e.g. Kato et al. 2014). The lower frame in Figure 5 shows the results of combining the two nights with the best period estimate (0.0576 d) determined from the PDM technique. Assuming that the observed light comes from the accretion disc, the zero point obtained in this case is HJD 2456537.6946 (time of inferior conjunction of the secondary). The period found using the PDM method yielded a value which is still close to the post-outburst state, but the sinusoidal shape is gone. There was only a small peak around phase 0.25. No double modulation with orbital period was found as would be expected in a bounce-back object. Further observations in the R band were obtained on 2018, July 18 covering three orbital cycles. Since the night was not photometric, we were unable to make absolute calibrations and only differential photometry is shown in Figure 6. No obvious orbital modulation was detected within the individual errors, which are rather large (≈ 0.03 mag).


Fig. 5
Upper panel: V light curves obtained near minimum light in 2013. Open dots are from September 2 and filled dots from September 25. The size of the points corresponds to the mean individual errors. The orbital phases have been taken from the ephemeris obtained in this paper. Lower panel: The two nights folded with the ephemeris obtained by using the PDM technique (see more details in text).


Fig. 6
Differential R light curve obtained at minimum light in July 18, 2018. The size of the individual errors (≈ 0.03 mag) is rather large. The orbital phases have been taken from the ephemeris obtained in this paper. No obvious orbital modulation, within the errors, is detected at this level.

4. Discussion

The values derived for the superhump period in Stages A and B allowed us to estimate the mass ratio of the system through the known “superhump excess” - q relations (Patterson 1998, 2011; Kato & Osaki 2013). Considering the lack of a reliable determination of the orbital period from spectroscopic observations, we assume here that Porb is equal to the period of early superhumps. We find A = 0.042(5) and B = 0.024(2). Thus, our estimates for the mass ratio are qA = 0.12(2) and qB = 0.10(1), respectively. Comparing these two values with the − q relation shown in Bakowska et al. (2017, Figure 19) we can see clearly that B is well within the expected value, while ǫA is not. This is further supported by using the updated Stolz & Schoembs (1984) relation in Otulakowska-Hypka et al. (2016, equation 4), which for our assumed orbital period gives = 0.019(10). Although we are inclined to use the Stage B results, we point out that both values are suggestive of a low-mass donor, although as pointed out in § 1 (e.g. Otulakowska-Hypka et al. 2016, see their Figure 7), empirical relations -q at low q values may carry large systematic uncertainties. For a typical white dwarf with mass 0.8 M (Zorotovic, Schreiber & Gänsicke 2011) and mass ratio q 0.1 , the mass of the secondary is very close to the substellar limit i.e. 0.072 M (Chabrier & Baraffe 2000). Further characterisation of systems like V1838 Aql will allow us to discern empirically where this limit lies for mass-losing donors.

It has been noted by many authors, both theoretically and observationally, that the CV orbital period distribution should present a sharp cut-off at about ≈ 80 min, usually termed as the minimum period (e.g. Rappaport et al. 1982; Ritter & Kolb 1998; Gänsicke et al. 2009). V1838 Aql has an orbital period of about 82 min, very close to the minimum period, which makes it diffcult to discern whether it is approaching to, or receding from, this minimum orbital period. Before asserting its true nature, we could look at some observational features in those CVs systems around the minimum orbital period. Most of these systems possess WZ Sge-like features. Their optical spectra are mostly dominated by the white dwarf and accretion disc itself, with no visible features from the donor. Since the rate of accretion is an order of magnitude smaller than that for systems before reaching the minimum period ( m ˙ 10 - 11 M ≈ yr−1), the accretion discs become very faint, and the broad absorption lines of the white dwarf become visible below   5000   Å ) e.g. WZ Sge (Howell et al. 2008). However, the donor’s observed properties should vary significantly for systems with the same orbital period but evolving towards or away from the period minimum. This is a consequence of the donor’s temperature steep relationship as a function of orbital period (e.g. Knigge et al. 2011). From the superhump analysis presented here, we suspect that the system is approaching the period minimum and thus, the donor is probably a late-M dwarf with an observed effective temperature of ≈ 2400 K (Knigge et al. 2011).

The donor of V1838 Aql is therefore an ideal candidate for NIR time-resolved spectroscopy (e.g. SDSS J143317.78+101123.3, Hernández Santisteban et al. 2016), which would render a fully independent measurement of the orbital period and the mass ratio, to confirm or reject its substellar nature. This is particularly important since few low-q systems have been observed in outburst and for which a dynamical measurement of their components is feasible. (Figure 3 in Kato & Osaki 2013). Thus, V1838 Aql could be a system to calibrate the empirical superhump relations in this poorly-explored region of parameter space.

5. Conclusions

We have presented a long-term study of the 2013 superoutburst of V1838 Aql from its peak to quiescence. Our main results are as follows:

  • The observed early superhumps suggest an orbital period of Porb = 0.05698(9) d, which locates V1838 Aql close to the minimum of the orbital period distribution in CVs.

  • From Stages A and B and early superhump periods, we found the mass ratio to be qA = 0.12(2) and qB = 0.10(1), respectively.

  • Based on the obtained values of the mass ratio, we claim that the donor in V1838 Aql is a low-mass star rather than a substellar object, and that the system is approaching the period minimum.

Given the long interval between outbursts in low q systems, it is of paramount importance to confirm by dynamical methods the orbital parameters of such systems. This would indicate which systems may be used in order to calibrate the empirical superhump excess relations.

The authors are indebted to DGAPA (Universidad Nacional Autónoma de México) support, PAPIIT projects IN111713, IN122409, IN100617, IN102517, IN102617, IN108316 and IN114917. GT acknowledges CONACyT grant 166376. JE acknowledges support from a LKBF travel grant to visit the API at UvA. JVHS is supported by a Vidi grant awarded to N. Degenaar by the Netherlands Organization for Scientific Research (NWO) and acknowledges travel support from DGAPA/UNAM. E. de la F. wishes to thank CGCI-UdeG staff for mobility support. We thank the day and night-time support staff at the OAN-SPM for facilitating and helping obtain our observations.

Supplementary material
References
Bąkowska, K., Olech, A., Pospieszyński, R., et al. 2017, A&A, 603, 72
Breedt, E., Gänsicke, B. T., Drake, A. J., et al. 2014, MNRAS, 443, 3174
Chabrier, G. & Baraffe, I. 2000, ARA&A, 38, 337
Coppejans, D. L., Körding, E. G., Knigge, C., et al. 2016, MNRAS, 456, 4441
de Miguel, E., Patterson, J., Cejudo, D., et al. 2016, MNRAS, 457, 1447
Gänsicke, B. T., Dillon, M., Southworth, J., et al. 2009, MNRAS, 397, 2170
Harrison, T. E. 2016, ApJ, 816, 4
Hernández Santisteban, J. V., Knigge, C., Littlefair, S. P., et al. 2016, Natur, 533, 366
Howell, S. B., Rappaport, S., & Politano, M. 1997, MNRAS, 287, 929
Howell, S. B., Hoard, D. W., Brinkworth, C., et al. 2008, ApJ, 685, 418
Itagaki, K. et al. 2013, CBET, 3554, 1
Kato, T. 2015, PASJ, 67, 108
Kato, T., et al. 2009, PASJ, 61, 395
Kato, T. et al. 2014, PASJ, 66, 90
Kato, T. & Osaki, Y. 2013, PASJ, 65, 115
Knigge, C., Baraffe, I., & Patterson, J. 2011, ApJS, 194, 28
Lenz P., & Breger M. 2005, CoAst, 146, 53
Littlefair, S. P., Dhillon, V. S. & Martin, E. L. 2003, MNRAS, 340, 264
Littlefair, S. P., Dhillon, V. S., Marsh, T. R., et al. 2006, Sci, 314, 1578
Neustroev, V. V., Marsh, T. R., Zharikov, S. V., et al. 2017, MNRAS, 467, 597
O’Donoghue, D., Chen, A., Marang, F., Mittaz, J. P. D., Winkler, H. & Warner, B. 1991, MNRAS, 250, 363
Osaki, Y. 1974, PASJ, 26, 429
______. 1994, ASIC, 417, 93
Otulakowska-Hypka, M., Olech, A., & Patterson, J. 2016, MNRAS, 460, 2526
Paczynski, B. & Sienkiewicz, R. 1981, ApJ, 248, 27
Patterson, J. 1998, PASP, 110, 1132
______. 2011, MNRAS, 411, 2695
Patterson, J. et al. 1996, PASP, 108, 748
Patterson, J. et al. 2005, PASP, 117, 1204
Paunzen, E. & Vanmunster, T. 2016, AN, 337, 239
Rappaport, S., Joss, P. C., & Webbink, R. F. 1982, ApJ, 254, 616
Ritter, H. & Kolb, U. 1998, A&A, 129, 83
Skillman D. R. & Patterson J. 1993, ApJ, 417, 298
Smak, J. 1993, AcA, 43, 101
Stellingwerf, R. F. 1978, ApJ, 224, 953
Stolz, B. & Schoembs, R. 1984, A&A 132, 187
Vogt, N. 1982, ApJ, 252, 653
Warner, B. 1995, Cataclysmic Variable Stars (Cambridge, MA: CUP)
Whitehurst, R. 1988, MNRAS, 232, 35
Zorotovic, M., Schreiber, M. R., & Gänsicke, B. T., 2011, A&A, 536, A42
Notes
Author notes
G. Bolt: CBA-Australia, 295 Camberwarra Drive, Craigie, Western Australia 6025, Australia.
R. Costero, C. Echevarría, A. Ruelas-Mayorga, and L. J. Sánchez: Instituto de Astronomía, Apartado Postal 70-264, México, CDMX, C.P. 04510 (jer@astro.unam.mx).
E. de Miguel: Departamento de Ciencias Integradas, Universidad de Huelva, E-21071 Huelva, Spain.
W. Goff: CBA-California, 13508 Monitor Lane, Sutter Creek, CA 95685, USA.
D. González-Buitrago: Departament of Physics and Astronomy 4129, Frederick Reines Hall, University of California, Irvine, CA 92697-4575, USA.
F.-J Hambsch: CBA-Mol, Andromeda Observatory, Oude Bleken 12, B-2400 Mol, Belgium.
J. V. Hernández Santisteban: Anton Pannekoek Institute for Astronomy, University of Amsterdam, Science Park 904, 1098XH Amsterdam, The Netherlands.
J. L. Jones: CBA-Oregon, Jack Jones Observatory, 22665 Bents Road NE, Aurora, OR, USA.
R. Michel: Instituto de Astronomía, Universidad Nacional Autónoma de México, Apartado Postal 877, Ensenada, Baja California, C.P. 22830, México.
J. Olivares: Laboratoire d’astrophysique de Bordeaux, Univ. Bordeaux, CNRS, B18N, allée Geoffroy Saint-Hilaire, 33615 Pessac, France.
A. Oskanen: CBA-Finland, Hankasalmi Observatory, Verkkoniementie 30, FI-40950 Muurame, Finland.
R. Sabo: CBA-Montana, 2336 Tailcrest Dr., Bozeman, MT 59718, USA.
D. Slauson: CBA-Iowa, Owl Ridge Observatory, 73 Summit Avenue NE, Swisher, IA 52338, USA.
W. Stein: CBA-Las Cruces, 6025 Calle Paraiso, Las Cruces, NM 88012, USA.
J. Ulowetz: CBA-Illinois, Northbrook Meadow Observatory, 855 Fair Lane, Northbrook, IL 60062, USA.
Table 1
Times of superhump maxima of V1838 Aql during the 2013 superoutburst.a

a Individual errors in these timings are not explicitly included here. b E (cycle number). c Superhump maxima expressed as HJD - 2,456,000. d O - C value (in cycles) according to the ephemeris. T_max (HJD) = 2,456,452.8035 + 0.058191 E.

Fig. 1
Photometric behaviour of V1838 Aql during its superoutburst in 2013. Top frame: global light curve of the observational campaign. Middle frame: amplitude variations of the superhumps along the superoutburst (see text). Bottom frame: O-C diagram for the superhump maxima in the HJD 449.6-498.6 day interval with respect to the ephemeris given in equation (1). The arrows indicate the (approximate) location of transitions between different regimes of superhump period variations (see text).
Table 2
Mean photometric periods and frequencies of the superhump modulation.*

* Values correspond to the 2013 superoutburst of V1838 Aql. The errors in parentheses correspond to the last significant figures.

Fig. 2
Upper frame: Power spectrum during days HJD 445-448 (early superhumps), showing broad peaks centred at 17.55 cycles d−1 and its first harmonic. Lower frame: Waveform of the early superhumps obtained after folding the data with Pesh = 0.05698 d. The zero phase is arbitrary.

Fig. 3
Upper frame: A 12-hour spliced light curve obtained on day HJD 455, dominated by large-amplitude common superhumps. The zero level in the figure corresponds to V ≈ 12.1 mag. Middle frame: Power spectrum during the common superhump era (days HJD 442- 466), with main peaks centred at frequencies 17.128 and 34.265 cycles d−1. Lower frame: Mean waveform of the common superhump obtained after folding the data on P = 0.058384 d. The zero phase is arbitrary.

Fig. 4
Upper frame: Power spectrum after the main eruption (day HJD ≥ 473) showing a strong peak centred at 17.244(1) cycles d−1 (post-outburst superhump). Lower frame: Mean waveform of the post-outburst superhump, obtained after folding the data on P = 0.057991 d. The zero phase is arbitrary.

Fig. 5
Upper panel: V light curves obtained near minimum light in 2013. Open dots are from September 2 and filled dots from September 25. The size of the points corresponds to the mean individual errors. The orbital phases have been taken from the ephemeris obtained in this paper. Lower panel: The two nights folded with the ephemeris obtained by using the PDM technique (see more details in text).

Fig. 6
Differential R light curve obtained at minimum light in July 18, 2018. The size of the individual errors (≈ 0.03 mag) is rather large. The orbital phases have been taken from the ephemeris obtained in this paper. No obvious orbital modulation, within the errors, is detected at this level.
Buscar:
Contexto
Descargar
Todas
Imágenes
Scientific article viewer generated from XML JATS4R by Redalyc