A&A 478, 815-822 (2008)
DOI: 10.1051/0004-6361:20077886
A. Dobrotka1 - A. Retter2 - A. Liu3
1 - Departement of Physics,
Faculty of Materials Science and Technology,
Slovak University of Technology in Bratislava,
Jána Bottu 25,
91774 Trnava,
The Slovak Republic
2 -
Departement of Astronomy and Astrophysics,
Penn State University,
525 Davey Lab, University Park,
PA, 16802-6305,
USA
3 -
Norcape Observatory,
PO Box 300,
Exmounth, 6707,
Australia
Received 15 May 2007 / Accepted 29 October 2007
Abstract
Aims. We present the period analysis of unfiltered photometric observations of V5116 Sgr (Nova Sgr 2005 #2) and search for superhump candidates in novae remnants.
Methods. The PDM method for period analysis was used. The masses of the novae componets we are estimated from the secondary mass - orbital period and primary mass - decline time relations.
Results. We find that 13 nights of V5116 Sgr observations in the year 2006 are modulated with a period of
d (
h). Following the shape of the phased light curves and no apparent change in the value of the periodicity in different subsamples of the data, we interpreted the period as orbital in nature. The binary system then falls within the period gap of the orbital period distribution of cataclysmic variables. From the maximum magnitude - rate of decline relation, we estimated a maximum absolute visual magnitude of
mag using the measured value of decline
d. The mass-period relation for cataclysmic variables yields a secondary mass estimate of about
.
We propose that V5116 Sgr is a high inclination system showing an irradiation effect of the secondary star. No fully developed accretion disc up to the tidal radius with the value lower than
cm is probable. The mass ratio was estimated in a few novae and the presence or absence of superhumps in these systems was compared with the mass ratio limit for superhumps of about 0.35. We find that, in the majority of novae with expected superhumps, this variability has not been found yet. Therefore, more observations of these systems is encouraged.
Key words: stars: novae, cataclysmic variables - stars: individual: V5116 Sgr - accretion, accretion disks
Novae are a subclass of cataclysmic variable stars. In these interracting binaries, the white dwarf is accreting the matter transferred from the secondary star. The accretion disc may be formed in the non magnetic case. The intermediate polar systems have a truncated disc and polar systems have magnetic field strong enough to prevent disc formation (see Warner 1995, for a review). The accumulation of a critical amount of accreted material onto the white dwarf surface results in a nova explosion. The distribution of orbital periods in cataclysmic variables shows a period gap between about 2 and 3 h, but novae do not show this lack of objects in the mentioned interval.
The source V5116 Sgr (Nova Sgr 2005 #2) was discovered by Liller (2005) on
2005 July 4.049 UT. The nova had a magnitude 8.0 on two red
photographs. An unfiltered CCD image from 2005 July 5.085 UT showed the object
at mag 7.2. The spectrum from 2005 July 5.099 UT shows H
with
the FWHM of
970 km s-1. The expansion velocity derived from the
sharp P Cyg profile was
1300 km s-1. The position of the nova was
measured by Gilmore & Kilmartin (2005) and Jacques (2005). Gilmore &
Kilmartin (2005) searched for the nova precursor, but no convincing candidate
has been found. Russell et al. (2005) performed a 0.8-2.5
m
spectroscopy of the nova on 2005 July 15. The object showed emission lines of
H I, He I, C I, N I, Ca II, and O I with an FWHM
2200 km s-1, and
He I showed a P Cyg profile at 1.0830 and 2.0581
m. No thermal dust
emission was observed.
After the nova explosion, the accretions disc is destroyed. The new disc is reformed by the stream of matter flowing from the secondary, interracting with itself and forming a ring in the circularisation radius. The disc starts to form by the viscous shearing. Matter losing angular momentum is moving inward and the excess of angular momentum is transported by the matter flowing outward. The inward moving matter touches the white dwarf and the disc is reformed (Pringle 1981). If the white dwarf is magnetic, the matter interacts with the magnetosphere and is then conducted by the magnetic field to the magnetic poles. The interaction of the flow with the poles of the rotating star is observed as a periodic signal with the spin period of the white dwarf. In the case of intermediate polars, the spin period is usually much shorter than the orbital period (Patterson 1994; Hellier 1996). The polar systems without a disc are synchronous rotators, so the spin period is equal or close to the orbital period (see e.g. Schmidt & Stockman 1991). Nova V1500 Cyg (polar system without a disc) changed the period from 0.141 d to 0.138 d and then stabilised at 0.140 d (Patterson 1978, 1979a). The difference of 1.8% between the rotation period of the white dwarf and the orbital period is ascribed to the effects of the nova explosion in 1975. The synchronisation can be then corrupted by the nova explosion, but the spin period remain still very close to the orbital period. This is in contrast with intermediate polars with disc presence.
Searching for periods in novae allows study of the orbital distrubution and evolution of these systems. Currently, there are about 50 novae with known orbital periods (Warner 2002) with typical values ranging from 2 to 9 h. The existence of the accretion disc or its renovation after the nova explosion is confirmed by the detection of the superhump period (see e.g. Retter et al. 1997; Kang et al. 2006a) or the spin period of a magnetic white dwarf in the case of intermediate polars (see e.g. Retter et al. 1998).
Superhumps are caused by a precessing accretion disc generally in systems with
mass ratio
(see Patterson et al. 2005,
for review), in which the disc radius reaches the 3:1 resonance. The mass
ratio indicates that systems with massive primaries and low-mass secondaries
are probable superhumpers. Novae in general possess high mass white dwarfs
(see e.g. Warner 1995; Smith & Dhillon 1998; Webbink 1990) and short orbital
periods sugest low mass secondaries. It is therefore meaningfull to search for
superhump variability in novae with short orbital period.
We have an ongoing program with small telescopes to observe novae and to
search for periodicities in their optical light curves. In this paper we
report detection of one periodicity in our photometric data (
d) of V5116 Sgr and discuss the presence of superhumps in nova
remnants. In Sect. 2 we present our observational material. The
long-term light curve with the period analysis of the data is presented in
Sect. 3 and in Sect. 4 we discuss the long-term
light-curve behaviour (Sect. 4.1), the results of the period
analysis (Sect. 4.2), and superhump search in nova remnants
(Sect. 4.3).
V5116 Sgr was observed during 13 nights in 2006. The observations included
45.2 h and 1256 data points in total. The list of observations is
presented in Table 1. The photometry was carried out with a 0.3-m
f/6.3 telescope coupled to an SBIG ST7E CCD camera. The telescope is located
in Exmounth, Western Australia. The pixel size of the CCD is
microns. This camera is attached to an Optec f5 focal reducer giving an image
field of view of
arcmin. The range of seeing for the data was 2.5-3 arcsec. The exposure times were between 30 and 90 s every 120 s, and no
filter was used. Aperture photometry was used in the reduction, with an
aperture size of 12 pixels (radius). As comparison star, we used GSC 7392 856
with V = 12 mag and as check star GSC 7392 292 with
V = 11.6 mag. The
instrumental magnitudes of the stars were derived from the SBIG CCDOPS
software and the magnitudes of the comparison and check stars from the GSC.
Figure 1 displays our light curve and the long-term light curve of the nova from the American Association of Variable Star Observers (AAVSO). In Fig. 2 we show the 4 best examples of our observational runs with the sinusoidal fit (with first harmonic) to the data using the period derived in this paper. These variations are seen clearly.
Table 1: The observational log.
![]() |
Figure 1: The AAVSO light curve ( upper panel) and the unfiltered light curve of our data ( bottom panel). |
Open with DEXTER |
![]() |
Figure 2: A selection of our observations (4 detrended best runs; 1, 2, 11, and 13). The solid curve is the sinusoidal fit (plus first harmonic) to the data using the period derived in this paper. |
Open with DEXTER |
By using the maximum visual magnitude of V =7.15, we measured the time
required for a decline of two magnitudes from maximum
d and three magnitudes from maximum
d
(Fig. 1 - upper panel). This makes V5116 Sgr a very fast nova
according to the classification given in Table 5.4 of Warner (1995). We
calculated the visual absolute magnitude at maximum brightness using Eq. (2) of Della Valle & Livio (1995) and obtained
mag.
The bottom panel of the Fig. 1 shows our unfiltered photometry for
the nova. The long-term variations with a time scale 30-40 days
with amplitude
1 mag
450 days after the outburst is
clear. Similar large-scale quasi-periodic brightness oscillations with
amplitude 1-1.5 mag months after the maximum were observed, for example, in
GK Per, V603 aql, DK Lac (see e.g. Warner 1995) or V4745 Sgr (Csák et al. 2005), with time scales of 5, 12, 25, and
20 days, respectively.
For the period analysis we used the PDM method (Stellingwerf 1978). First, the
data were detrended by substracting the linear fit from each night. A simple
visual inspection of the runs suggests that the data are modulated with a
periodicity of 1.5 h. This value was then identified as the first
harmonic of the real periodicity, making the exact value of the frequency
d-1, which corresponds to the periodicity
d (
h). We calculated the errors in the
frequency from the half width at the half minimum of the peaks in the PDM power spectrum.
The results from the period analysis of all the data are depicted in
Fig. 3. The suspected periodicity (
d-1,
d) is seen with its 1-d aliases, and the dominant peak
is the first harmonic (
d-1). For testing
which value is the real signal and which the alias, we performed an O-C
diagram with a linear fit to the data (Fig. 4; upper panel - O-C for
the frequency 8.08 d-1 and lower panel - O-C for the frequency
9.08 d-1). The lower panel shows a larger scatter of the O-C values
compared with the upper panel. This indicates that the real signal is
d-1 and the dominant peak at
16 d-1 is the
first harmonic of the real signal.
To test the reality of the peaks, we performed a power spectrum of the
synthetic data using a sinusoids with the derived periods ( including
fundamental f, first harmonic
and
described later)
and phases with amplitudes found from the fits to the data sampled as the
original data with introduced Gaussian noise (Fig. 3e). The
similarity with the original periodogram is clear. Furthermore, we performed a
period analysis of two different subsamples of the data (first subsample -
runs 3-8; second subsample - runs 9-13) depicted in
Fig. 5. The frequency is present together with its first harmonic
in both datasets. The values are
d-1 and
d-1 (first and second subsample, respectively), which is within
the errors in both cases. We also performed an analysis of the check star and
did not identify any periodicity associated with its data.
![]() |
Figure 3:
Power spectra of all the data. Panel a) - raw detrended data (night
means substracted); panel b) - data after the substraction of the harmonic
frequency
![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 4:
O-C diagram. Upper panel - for frequency
![]() ![]() ![]() ![]() |
Open with DEXTER |
There is a strong peak near the value 8.58 d-1. This signal is not
real as seen in Figs. 3b,c. After the substraction of
(Fig. 3b), this peak disappeared, and after a further substraction of f (Fig. 3c), no peak is dominant in this frequency region. This
feature also appears in the synthetic power spectrum (Fig. 3e), so
it is definitely an alias of f.
Another strong feature is a peak at 35 d-1. Substracting the
harmonic frequency
from the data, after the substraction
described in the previous paragraph, this signal dissapared
(Fig. 3d). It is another artefact of the real signal f. No other
dominant peaks are visible with frequencies higher than 40 d-1. No
other frequency was then found.
![]() |
Figure 5: Power spectra of two different subsamples of the data. Upper panel - runs 3-8; lower panel - runs 9-13. The dashed line represents f = 8.08 d-1. |
Open with DEXTER |
In Fig. 6 we show the light curve of V5116 Sgr folded on the
0.1238 d period. The points are averaged magnitudes in 50 equal bins that
cover the 0-1 phase interval. The full amplitude of the mean variation is
mag for primary minimum and
mag for secondary
minimum. Using the sinusoidal fit to determine the amplitude is inadequate,
because of the non sinusoidal shape of the light curve, we therefore estimated
the amplitude by measuring the difference between extrema. The bars are the
errors of the mean value. The light curves folded on the periodicity found
show typical features of an eclipsing system.
We determided the epoch of minima by fitting low-order polynomials to selected
parts of the light curve. The best-fitted (using sinusoidal fits) ephemeris of
the periodicity is
![]() |
Figure 6: The light curve of V5116 Sgr folded on the 0.1238 d period and binned into 50 equal bins. Panel a) - all the data; panel b) - runs 3-8; panel c) - runs 9-13. |
Open with DEXTER |
The long-term light curve of the nova V5116 Sgr (Fig. 1 - bottom panel) shows a transition from smooth decline to probable oscillations. Several models have been suggested for this ``transient'' phase. One of them explains the transition as the time when the accretion disc is re-established (Retter 2002). The author proposes a connection between this phase and intermediate polars. There is no indication of the magnetic nature of the white dwarf in V5116 Sgr yet. A spin period of the white dwarf indicating the intermediate polar type of this ``transient'' phase in novae was detected, for example, in V4745 Sgr by optical observations (Dobrotka et al. 2006a), V1494 Aql (Drake et al. 2003), and V4743 Sgr by X-ray observations (Ness et al. 2003). If the ``transient'' phase interpretation of Retter (2002) is applicable to V5116 Sgr, the accretion must be then re-established and the accretion disc should be present.
We have identified one periodicity in the light curve of the nova V5116 Sgr
about 15 months after the maximum brightness. Its value is
d (
h). The upper limits of the period
difference between two subset of data with
16 days of mean distance
analysed in Fig. 5 is
d, which gives
.
Following Patterson et al. (1993), the mean variation in the
superhump period in SU UMa superhumping sytems is
(see their Table. 1). For the recurrent nova VY Aqr, the authors
derived a variation of
.
Our value for
probably comes from the errors in the period measurements rather
than from real period changes. The period then seems to be constant during our
observations, and the shape of the folded light curve suggests a primary and a
secondary eclipse. We therefore propose that this periodicity is the orbital
period of the binary system. Such a period is at the lower end of the mostly
populated region of orbital periods in novae (Warner 2002). The dominant first
harmonic frequency in the power spectra is a result of the clear structure
primary - secondary eclipse in the folded light curve, which suggests a high
inclination angle of the binary system. Using the orbital period of the system
and Eq. (9) from Smith & Dhillon (1998), we obtain a rough estimate for the
secondary star mass of
.
Using a mean
white-dwarf mass of
from Smith & Dhillon
(1998), we find a mass ratio (secondary/primary mass) of
.
After the nova explosion, the hot white dwarf may heat up and irradiate the
cooler companion. The observed orbital light curve of the nova can
be the result of the aspect variations or eclipses of the secondary due to
heating from the hot primary, and the asymmetry in the pulse profiles could be
produced once the shape of the secondary approaches that of a tear drop
model. The irradiation effects in classical novae can also be detected long
after the outburst stage (e.g., V1500 Cyg; Sommers & Naylor 1999, DN Gem;
Retter et al. 1999). Two eclipse-like features are present in the folded light
curve of V5116 Sgr (Fig. 6). The shape is similar to the light
curve of V2540 Oph in 2003 (Ak et al. 2005) with a dip at phase
0.5 (
amplitude of the primary minimum for V2540 Oph and
for V5116 Sgr). The authors conclude that V2540 Oph is likely a high
inclination system either showing an irradiation effect or having a spiral
structure in its accretion disc. Woudt & Warner (2003a) note that one of the
following requirements must be fulfilled for the large-amplitude orbital
modulations to be seen in the light curve of a recent nova in which the
accretion disc does not dominate the luminosity of the system: 1) the disc is
foreshortened but the secondary is seen (a high inclination angle), 2) the
disc has small dimensions, 3) there is no disc (case of polars). In the case
of V5116 Sgr, the polar interpretation is possible because of the orbital
period distribution of polars that peaks below 5 h (Warner 1995), but
nothing else indicates this option. The small dimension of the disc is
supported by 1) the short orbital period, and 2) the post nova stage when the
disc is reforming after explosion.
A comparison to two novae within the period gap IM Nor and DD Cir with present irradiation effect (Woudt & Warner 2003b, c) can be made. The light curve of V5116 Sgr is different from those of IM Nor and DD Cir but similar to V2540 Oph. IM Nor and DD Cir showed very small dips at phase 0.5 interpreted as partial eclipses of the irradiated secondary by the disc or matter stream. The light curve shape depends on the disc radius and on the inclination angle. The deep secondary eclipses in V5116 Sgr require a large disc or a high inclination. The strength of the irradiation effect (white dwarf post-outburst temperature) can also play a significant role. The differences in the phase of the maxima in Figs. 6b and c can be due to spiral structures in the disc as mentioned before for the nova V2540 Oph.
The reconstruction of an accretion disc after the nova eruption is indicated
by the discovery of the white dwarf spin or by the superhump period. A
probable spin period of the white dwarf was detected 1 yr after the
outburst in V4745 Sgr (Dobrotka et al. 2006a),
2.75 yr in V4743 Sgr
(Kang et al. 2006b), and
15 months after the maximum in V1425 Aql
(Retter et al. 1998). Several systems show superhumps as early as two and a
half months after the outburst like V4633 Sgr (a spin period is another option
in this case, Lipkin et al. 2001) or two years after the outburst like
V1974 Cyg (Retter et al. 1997). In our V5116 Sgr light curves extending
470 days after outburst, we did not find any photometric variations
consistent with white dwarf spin modulation or superhump properties. We cannot
say anything about the magnetic nature of the white dwarf, but the component's
mass ratio indicates that the superhump existence is likely. The mass ratio
using the primary mass average is not decisive. Using the derived secondary
mass 0.26
(Sect. 4.2), the mass ratio is
lower than 0.35 for a primary mass higher than 0.74
.
Taking
all known primary masses in cataclysmic variables (Ritter & Kolb 2003), 63%
have higher or equal mass than 0.74
.
This probability value
is not enough to make certain that superhumps are expected inV5116 Sgr, but
the search for such variability could be fruitful; however V5116 Sgr is a very
fast nova. According to Eq. (13) from Livio (1992) and taking
d derived in this paper (Sect. 3.1), we obtain a mass
estimate of the white dwarf of
and thus
.
The presence of superhumps is then
expected if a disc is fully developed up to tidal radius.
The superhump period is a few percent longer or shorter than the orbital period. The possibility that the periodicity found in this paper is a superhump is rejected following the stability discussion and the shape of the folded light curve (Fig. 6). The mean shape of superhumps is typically an asymetric sinusoid and our data show typical eclipse like features.
The presence of the disc in cataclysmic variables depends on the mass loss from the secondary. The material supplied from the secondary within the period gap depends on the strength of the magnetic braking driving the secondary out of thermal equilibrium (mass loss time scale shorter than the thermal time scale). The stronger the braking, the wider the gap and the higher the upper end of the period gap. If the magnetic braking is low enough, the mass loss time scale may never become shorter than the thermal time scale. In the case of novae, another condition is important. The nova explosion heats the secondary that leads to an enhanced mass transfer. Therefore the complete absence of accretion discs in novae within the period gap is not expected.
Table 2: Selected nova systems.
The orbital period 2.462 h of the mentioned IM Nor (Woudt &
Warner 2003b) yields a secondary mass estimate
of 0.20
following Smith & Dhillon (1998). The decay time
d (Kato et al. 2002) yields a primary mass estimate of
0.86
(Livio 1992), giving a mass ratio
0.23. Woudt & Warner (2003b) conclude that they had not observed the
superhump modulation. This conclusion with the mass ratio safely lower than 0.35 indicates that the existence of a fully developed disc is not
probable. DD Cir (Woudt & Warner 2003c) has a period of 2.339 h yielding
a secondary mass estimate of 0.18
.
The decay time
d (Liller 1999; no other estimates are available) places the
object in the class of very fast novae (similar to V5116 Sgr). The decay time
d following the equation
from Warner (1995) gives a rough white dwarf mass estimate of
1.16
(mass ratio <0.2). Woudt & Warner (2003b)
interpreted the deep eclipse as obscuration of the disc and did not detect any
periodicity near the orbital period in the period analysis. Therefore,
following the tidal instability model, it is again probable that the disc is
not fully developed up to 3:1 resonance radius 3 years after the outburst;
however, the authors argue that the disc radius is 47% of the orbital
separation. Using the third Kepler law, we obtain a disc radius of
cm. The disc in DD Cir is then large enough to reach the 3:1
resonance radius calculated from the Eq. (3.39) in Warner (1995);
cm. The result marginally contrasts with the absence of
superhumps. With the derived values for V5116 Sgr, we get
cm as an upper limit of the disc radius, when considering the
absence of superhumps into consideration. The boundaries of the eclipse in our
light curve are not as clear as in the DD Cir case, therefore not suitable to
estimate the disc radius.
Comparison to other systems follows the same way as in the case of IM Nor and
DD Cir. We took component masses from the Ritter & Kolb's catalogue (2003),
and other masses we are estimated from the orbital period (using Smith &
Dhillon 1998) and
time (using Livio 1992). We rejected systems
with unknown or non measurable
or
time (peculiar light
curve, strong oscillations in the decline) without catalogue mass values (no
information on mass ratio, ex: RS Car), magnetic novae (AM Her systems without
disc, ex: V1500 Cyg), systems with orbital period >10 h (Smith and Dhillon
fitting not adequate, systems with evolved secondary are then excluded too,
ex: DI Lac, V841 Oph, V723 Cas, GK Per), systems with insufficient photometric
observations to detailed period analysis (ex: V500 Aql, V446 Her, HZ Pup,
DY Pup, CT Ser). The final list of systems is in Table 2.
![]() |
Figure 7:
Selected nova remnants in the mass ratio (q) - orbital period
(
![]() ![]() ![]() |
Open with DEXTER |
The results are depicted in Fig. 7. The lower panel only shows the
systems with known catalogue mass values and the upper panel covers the
addition of estimated mass values in this work. The critical mass ratio of
0.35 is also shown. The results presented in the lower panel are as expected
from the mass ratio limit for superhumps, but the statistical set is
small. Adding other mass estimates changes the situation. Seven systems below
or close to the limit 0.35 show superhumps but 14 systems do not. The
estimated mass ratio has a mean difference of 0.2 from the catalogue
values. The differences are distributed randomly. Therefore there is no reason
to suspect that all mass ratios are systematically underestimated. The seven systems with detected superhumps are safely located in comparing to the 0.35 critical value, as expected by theory, and are only a third of all systems
falling below or close to the mass ratio of 0.35. The majority of all systems
occupy the orbital period interval of 2-4 h. Following Patterson
(2005), this interval has a 40-90% probability of observing
superhumps. Seven systems versus 14 from our investigation yields only
30%.
Superhumps were observed during the outbursts of short orbital period dwarf novae (SU UMa stars) where the disc reaches the 3:1 resonance radius (superoutburst case). The disc during this active stage (outburst - active and hot stage, quiescence - not active and cold stage) reaches larger radius than in quiescence (see Lasota 2001, for review). The irradiation of the secondary by the central white dwarf in nova remnants can produce an enhanced mass transfer rate that is strong enough to keep the disc in this stable hot active stage (as in the nova like stars) with the larger diameter (Hameury & Lasota 2005). In addition, superoutbursts in SU UMa stars, when superhumps are observed, can be explained by the irradiation and the enhanced mass transfer (Smak 1995, 2004; Hameury et al. 2000; Schreiber et al. 2004). Therefore, permanent superhumpers are possible after a nova outburst. Some of the possible interpretations of the superhump lack can be: 1) not systematically applicable tidal instability theory (tidal torques are probably not the main and only condition of superhump existence as concluded by Hameury & Lasota 2005), 2) disc radius not developed up to the tidal radius, 3) superhumps present but not detected because of insufficient sets of data.
By using the maximum visual magnitude of V = 7.15, we measured the decay
time
d from the long-term light curve. This makes
V5116 Sgr a very fast nova. We estimated a visual absolute magnitude in
maximum of
mag. We found a period of
d (
h), which we interpret as the orbital
period of the binary system. We propose that V5116 Sgr is a high inclination
system showing a strong irradiation effect of the secondary star. No fully
developed accretion disc up to the tidal radius is possible with a radius
lower than
cm.
Searching for superhumps in nova remnants using the limiting value for mass ratio 0.35 and estimated mass values shows that, in the majority of systems with expected superhumps, this variability has not been found yet. Looking for superhumps in these systems could be fruitful.
Acknowledgements
This work was supported by the Slovak Academy of Sciences Grant No. 2/7011/7. We acknowledge with thanks the variable star observations from the AAVSO International Database contributed by observers worldwide and used in this research. We also acknowledge Michael Friedjung for his valuable comments.