Issue |
A&A
Volume 521, October 2010
|
|
---|---|---|
Article Number | A84 | |
Number of page(s) | 11 | |
Section | Interstellar and circumstellar matter | |
DOI | https://doi.org/10.1051/0004-6361/200913936 | |
Published online | 25 October 2010 |
A long term spectroscopic and photometric study of the old nova HR Delphini![[*]](/icons/foot_motif.png)
M. Friedjung1 - M. Dennefeld1, -
I. Voloshina2,
1 - Institut d'Astrophysique de Paris -UMR 7095,
CNRS/Université Pierre et Marie Curie, 98bis boulevard Arago, 75014
Paris, France
2 - Sternberg Astronomical Institute, 119992 Moscow, Russia
Received 21 December 2009 / Accepted 17 June 2010
Abstract
Context. The Nova HR Del, discovered in 1967, was found to
be exceptionally bright in the optical and UV during the whole lifetime
of the IUE satellite (ending in 1996) and appears to still be extremely
luminous today. The reason for this continuing activity is not clear,
but continuing weak thermonuclear burning might be involved.
Aims. We therefore need to better understand the processes at play in HR Del.
Methods. HR Del was thus monitored over several years, both in broad band photometry and spectroscopically in the H
spectral region.
Results. The profile of the H
line shows two components: a narrow, central component and broader
wings. The former is most easily understood as the result of an
accretion disk whose geometry might lead it to partly occult certain
regions of itself. That component shows something like an S wave with
an orbital phase dependence, suggesting that it could be due to a spot
bright in H
.
The wide component must come from another region, with a probably
non-negligible contribution from the material ejected during the 1967
outburst. Non-orbital variations of the H
equivalent width were found both on long and short time scales. Similar
variations were found in the photometry, showing a component with a
clear dependence on the orbital phase, however with no obvious relation
with the H
variations.
Conclusions. The orbital part of the photometric variations can be explained by irradiation of the companion, while the properties of H
are explicable by the presence of an accretion disk and a spot bright in H
.
Key words: novae, cataclysmic variables - circumstellar matter - stars: individual: HR Del - stars: individual: V723 Cas - stars: individual: V603 Aql
1 Introduction
HR Delphini (or Nova Del 1967) was a bright classical nova with, however,
unusual properties. It was discovered on July 8 1967 by Alcock (see
Candy 1967),
when it had brightened to a magnitude near 5.5 from a pre-outburst
magnitude of about 12 at the beginning of June. After having remained 5
months near this
magnitude, it brightened again to a short peak at mv of 3.5 in December 1967,
which probably corresponds to what is considered to be the maximum for classical novae.
This was followed by an extremely slow irregular decline, as can be seen for instance
in Sanyal (1974).
During the stage before the brightest peak, the velocity of the main
absorption component continuously decreased, unlike the usually
observed behaviour of velocities in classical novae after maximum. This
reinforces the idea that the short peak in December 1967 was indeed the
``real'' maximum. The unusual behaviour
of this nova in such an exceptionally long pre-maximum stage might,
according to
Friedjung (1992), be explainable by the presence of a wind,
which was optically thin in the continuum, and with a decreasing velocity during that
pre-maximum stage, unlike the velocities of the optically thick winds normally seen in
classical novae after maximum. HR Del may, indeed, only have
marginally satisfied the conditions needed for thermonuclear runaway. This could be due
to the very low mass (0.55-0.75 )
of the white dwarf component of the binary
(Selvelli & Friedjung 2003). A few other
classical novae such as V723 Cas (1995) may have a similar nature.
In addition, this nova showed signs of unusual activity both before and long
after its outburst (Selvelli & Friedjung 2003). This can be seen by
its exceptional luminosity, compared with that of other old novae both from IUE
observations in the UV and from its optical apparent brightness. The time
variable P Cygni profile of the CIV resonance doublet also indicates a very
high wind velocity of the post-nova of the order of 5000 km s-1. The deduced values of distance of 850 pc and
E(B-V) = 0.16 led to a UV luminosity
of 56
and a visual absolute magnitude of
Mv = 1.85. The latter
value has to be corrected for the inclination of the system: assuming that the
radiation came from a gravitationally heated accretion disk (then having an
accretion rate of
yr-1 onto a white dwarf with a radius of 0.0125
), the correction leads to an absolute magnitude of
Mv = 2.30. Such observations suggested unusual activity; a possible explanation
was continuing very weak thermonuclear burning, with the radiation of the white dwarf being
reprocessed by an accretion disk (and the companion star),
a theoretically possible situation when the white dwarf component's
mass is low (Sion & Starrfield 1994). In this connection we can
note that V723 Cas was still an active X-ray source more than 12 years after the outburst (Ness et al. 2008), but in the case of HR Del, only weak X-ray
emission was reported at early times when instruments were much less sensitive
than now (Hutchings 1980, with
erg s-1 after correction for our adopted distance), and nothing more recently.
It is possible to obtain a better clue to the luminosity by estimating
the ionisation of the nebular remnant by the central source. The
remnant of HR Del
has been studied in a number of papers, e.g. Harman & O'Brien (2003).
Very recently, Moraes & Diaz (2009) have studied the clumpy structure of the nebula. Detailed photoionisation modelling suggests that there is a
disk shaped ionising source, with a luminosity of 1036 erg s-1 and a temperature of 65 000 K. The two sides of a disk, emitting like a black body at such a temperature, would then
have a radius of
cm or 0.2
.
Such a relatively low temperature might explain why the X-ray emission is weak . A total shell mass
of
was also found by these authors, at least half being contained in the neutral clumps.
In the light of these results, time variations of HR Del need to be followed
by ground-based observations, as they could reveal some signs of exceptional
activity. It was with this in mind that we observed the region of Hwith
a high spectral resolution over several years, since this line could
give particularly significant results. Parallel broad band photometric
monitoring was also performed. Some preliminary results of this
observing campaign have already been presented by Friedjung et al.
(2005).
2 Observations
2.1 Photometry
Photometric observations of HR Del were made in order to detect possible
long-term variations, and also to verify the orbital phases for the
spectroscopic observations. They were carried out with two different
instruments: firstly, the single-channel photoelectric photometer on
the 60-cm telescope of the Sternberg Astronomical Institute in Crimea was
used from 2002 to 2009, with continuously changing
Johnson UBV filters and a 10 s integration in each filter. The usual
method of differential measurements was applied. A local standard was
used permanently during observations and
calibrated against the well-known standard star SAO 106418
(magnitudes
).
More than 1600 measurements of HR Del were obtained in
3 bands during the total time of our observations, with usually
2-3 measurements per night. These data were used to determine the
long-term behaviour of HR Del, which is shown in Fig. 1. They are accurate to 1% in the V and the B bands and
in U band.
![]() |
Figure 1: Overall behaviour of HR Del in the V band during 2002-2008. Photoelectric V magnitudes are shown versus Julian Day. The arrows indicate the epochs of the spectroscopic observations. |
Open with DEXTER |
Secondly, CCD observations of HR Del were obtained with the goal of studying
variability on time scales from minutes to days. Several runs of a few nights
were made in August and October 2005, during the period of July-October 2006,
then during October-November 2008, and finally during three nights in October
2009. Most of these observations were carried out in the V band with the
60 cm telescope of the Crimean station of the Sternberg Astronomical Institute,
using an Apogee 47 CCD detector with
m pixels, and a field
of view of
arcmin. During 2008, in addition to the above instrument used both in the V and R bands, the 38cm telescope of the Crimean Astrophysical
Observatory was also used from November 10 to 13, but only in the R band,
with a CCD-ST7 camera (
m pixels). In 2009, a Pictor
MEADE-416 camera was also used on the 50 cm Sternberg telescope in the V band.
All these data were obtained in a
binning mode to reduce the read-out time. The duration of individual
observational series varied from 2.5 up to 4.5 h. A star close to
HR Del in the field was always used as a local standard (a few suitable
stars were used, among them the one used for the
photoelectric observations), and repeatedly calibrated using well-known
comparison stars during a few nights with perfect weather. Two other
stars, one with coordinates
and
and another with
and
,
were used as control stars.
More than 4200 measurements of HR Del were obtained this way in the V band during 2005-2008, covering more than 30 nights, while about 700 datapoints were obtained in R over 9 nights.
The CCD measurement errors do not exceed
% in the V band, and slightly less for the R band.
The details of the observations are presented in Table 1. Examples of daily light curves are presented in Fig. 2. A more detailed description of these observations will be presented in a companion paper (Voloshina et al., in prep.).
![]() |
Figure 2: Examples of CCD observations for different dates: CCD-V magnitudes versus phase. |
Open with DEXTER |
2.2 Spectroscopy
The spectroscopic monitoring was done at the Observatoire de Haute-Provence (OHP) in
France, using the 1.52 m telescope, equipped with the Aurélie high-resolution
spectrograph (Gillet et al. 1994). Light passing through a 3'' circular entrance
aperture then encounters a modified Bowen-Walraven image slicer (Walraven
1972), whose 5 slices of 0''6 each are
imaged onto a
EEV CCD with 13.5
m pixels.
To gain in sensitivity,
all the corresponding lines of the CCD were summed before read-out,
producing a single-line spectrum as output. However as a consequence,
the possibility of identifying and cleaning cosmic rays hits before
final reduction disappears. The spectra were reduced following standard
procedures: bias and flat-field corrections, and wavelength calibration
using a thorium-argon internal lamp. In view of the narrow spectral
range covered and used, no spectral response correction was performed.
Many of the observations (but not all) were done in service modes, and
used a 1200 l/mm grating giving a dispersion of 7.6 per mm and a spectral resolution of 0.32
(0.10 per pixel). Typical individual exposure times are 30 min,
but the sequencing of observations and phase coverage are irregular, a tribute to the service mode operation. The details of
observations are given in Tables 2 and 3 .
An example of a high-resolution spectrum is given in Fig. 3.
![]() |
Figure 3:
Example of a high dispersion spectrum, obtained on November 1, 2003 at the OHP. The core and the low-intensity wings of H |
Open with DEXTER |
A few low-dispersion spectra were obtained towards the end of the
campaign to check the general aspect of the object. We used the OHP
1.93m telescope in December 2007, the Carelec spectrograph and a
300 l/mm grating giving 1.8
per pixel and a 6.5
resolution with the 2'' slit used.
The same setting was used again in December 2008.
In June 2008, we used the IDS at the INT 2.5 m telescope in La Palma. There, a 630 l/mm grating gave 0.9
per pixel and 2.7
resolution with
a 1''5 slit. Those spectra were reduced with standard procedures, including
the use of standard star observations to correct for the wavelength-dependent spectral response. In
the case of the INT, however, only the central part of the spectrum can be
exploited due to strong vignetting. Summation was done over the full H
extent along the slit and the spectra are displayed in Fig. 4.
![]() |
Figure 4:
Low-dispersion spectra. Lower part: OHP, Dec. 2007; Middle: INT, June 2008. Upper part: OHP, Dec. 2008. Small offsets were made between spectra (see text). The units in the ordinates are 10-16 erg/cm2/s/ |
Open with DEXTER |
3 Results and analysis
3.1 Analysis of the photometry
As can be seen from Fig. 2, daily light curves of HR Del show
night-to-night variations and aperiodic stochastic variations, that is
so-callled flickering.
Because they contain long time series of observations at similar
epochs, the CCD data were used to search for the orbital period of the
light curve. Although covering a longer total timespan we may note that
the photo-electric data are too scarcely
sampled for that purpose. The analysis was made primarily in the V band,
where more data are available, but the R band data were also used to check
for possible wavelength-dependent effects. They were analysed with the
``Period04'' PC code from Lenz & Breger (2004) and with a programme
based on the Deeming (1975) method, in the frequency range
0-50 d-1. The resulting periodogram is displayed in Fig. 5:
the most prominent peak is for a period of 0.214164 days, corresponding,
within the uncertainties, to the period obtained by Kürster and Barwig
from their spectroscopic analysis (0.214165), the other peaks being daily
aliases. Our observations folded with this period are shown in
Fig. 6. The scatter that remains in this figure comes from
night to night variations and possibly also from flickering.
Attempts to fold data with periods corresponding to other peaks in the
periodogram lead to unacceptable fits, thus confirming the real period.
If we approximate this light curve with a sine function (thick, red line), the resulting semi-amplitude is
.
![]() |
Figure 5: Power spectrum of the 2008 CCD observations. |
Open with DEXTER |
![]() |
Figure 6: Light curve of HR Del versus photometric phase, in the V band, from observations obtained during 2005-2006. The averaged orbital sine wave is superposed onto the observed data. |
Open with DEXTER |
No other short-term variations were found in these observations besides the orbital ones. To better detect long time-scale variations, such as those found by Bianchini (1987) in some old novae, our data were corrected for orbital variations. There is a suspicion of long-term variations with a possible period of around 1500 days (that is about 4 years), already apparent in Fig. 1, but several more years of observations would be needed to confirm this.
![]() |
Figure 7: Equivalent width variation with orbital phase in 2003 November. Filled circles (above) represent the total equivalent width, triangles (middle) the central component and the squares (lower) the wings alone. |
Open with DEXTER |
![]() |
Figure 8: Equivalent width variation in 2005 mid September. Symbols as above. |
Open with DEXTER |
3.2 Analysis of the H
spectral region
The H
profiles show a well-defined central component, with
wider wings of lower intensity (Fig. 3),
suggesting formation in different regions. In view of this,
measurements were made of both the total equivalent width and the
line centre alone.
The time of each spectrum (at mid-integration), the
corresponding orbital phase according to the ephemeris of Kürster & Barwig
(1988) and the measured equivalent widths of H
(both total, and central component only) are given in Tables 2 and 3. Examples of the variation of the equivalent
width of the different H
line components with Julian date and with Kürster and
Barwig's orbital phase are shown for two epochs in Figs. 7 and 8, over a time interval not lasting more than a few days.
Clear differences are seen between these two examples: on one hand, the average value
of the H
equivalent width is much lower at the first epoch; on the other, the
scatter of the variations with phases is more important during the
second epoch. Such changes are seen over the whole duration of our
observing campaign.
One concern is a possible, time variable, contamination by emission from the ejected
nebula in the H
line due to the central region (plus a much smaller
contribution from the fainter, adjacent [N II] emission).
To estimate the possible contribution from the nebula, several tests have been made.
A detailed map of the nebular emission was obtained by Harman & O'Brien
(2003), combining HST images with ground-based, high-resolution
spectroscopy. They show that the nebular emission has an elliptical shape, with
major and minor axes of roughly
arcsec, with a large
H
contribution coming from a bright rim. These dimensions are much
larger than the 3'' entrance aperture, centred on the star, used for our
high-dispersion spectroscopy, so that the contribution from the nebula
should not be significant in the central H
component.
The structure of the nebula, as seen in the HST images of Harman & O'Brien
(2003), is however clumpy, and several bright knots, located inside
the ellipsoid, had been identified. To measure their possible contribution to
the H
flux from the central object (integrated over a circular aperture,
like the 3'' used), it would be sufficient to measure the intensity of those blobs
relative to the one of the central star. Unfortunately, the central object is
saturated in Harman's data, so we cannot estimate this ratio. We
instead used several of our ground-based images, obviously of lesser spatial
resolution, but where some of the brightest knots can be identified, and found
that the intensity ratio between the brightest knots and the central star is typically 1 to 50 in H
.
It would therefore require several knots within the aperture to make a 10%
contribution, and only a few are seen in the whole ellipsoid in the HST map.
![]() |
Figure 9: Equivalent width variation of the central component with phase, for four different observing periods. Triangles are for 2003 November, crosses for 2005 September, open stars for 2005 October and filled polygons for 2008 August. |
Open with DEXTER |
We also used the spatial resolution along the slit in our low-dispersion
spectra to estimate the change in H
equivalent width when including more and more outer regions, and found that 10% was a conservative
upper limit for a variable contribution from the nebula (see the
discussion below of the low-dispersion spectroscopy). Finally, the
measured velocities do not correspond either, at least for the central
component, to what is expected from the nebula (see discussion later).
In view of all this, we consider H
emission line equivalent
width variations of more than 10% as intrinsic, as can be seen in
the figures. The low level of the contamination is understandable,
following the expansion and fading of the nebula, as even our earliest
observations occurred about 35 years after the outburst of the
nova.
As the line wings are always faint, the value of the total H
equivalent
width is dominated by the value of the narrow central component, so the
variations in these two components are similar;
however, it is not clear to what extent any of the orbital phase variations
are significant. In some cases, like 2003 November (Fig. 7) or
2005 October, the variations are small: for instance, on November 1-2,
2003, the ratio between the highest and lowest values of equivalent width is only 1.13. This can
give a limit to deviations from circular symmetry at
that time. In other cases, e.g. 2005 September (Fig. 8),
or 2008 August, the variations, and the scatter are much larger. As more
observing points are available in those longer series, these larger variations are
probably indicative of rapid variations from night to night.
![]() |
Figure 10: Series of profiles obtained during 2002-2004, averaged in ten bins of 0.1 in phases. The phases are going from 0 at the bottom (bin 0 to 0.1) to 1 at the top (bin 0.9 to 1. |
Open with DEXTER |
![]() |
Figure 11: Series of averaged profiles obtained in 2005 (binning and phases as above). |
Open with DEXTER |
There are, however, significant differences in average equivalent width from one epoch to the next, as is illustrated in Fig. 9 for the central component: lower values are found for the 2003 epoch, with comparatively stable values over different phases, while higher values are found in September 2005, and a much larger change with phase in 2008.
Quite large, longer time-scale variations of the total and line centre equivalent widths are seen (e.g. Fig. 13, for the central component). The values are clearly lower in 2003 and 2004 than earlier or later. Such variations would be hard to explain by contamination from varying contributions of flux from the ejected nebula. For the wings, no clear systematic variation with time of their equivalent width can be established, as the scatter is large, but changes in their profiles and/or terminal velocity are definitely seen.
The central component of the H
profile is assymmetric and varies with
orbital phase, suggesting something like the well known S wave of cataclysmic
variables. Comparison of profiles from successive spectra, obtained at some
observational epochs when fairly long time series of observations could be
made, reveals this effect.
Three examples are displayed, from the beginning, middle and end of the campaign,
in Figs. 10-12, respectively, where series
of spectra from the same epoch have been binned over 0.1 intervals of phases.
They show the permanency and the stability of the phenomenon. The detailed profiles
may not look exactly the same, for a given bin, at each of the three epochs displayed,
but these small differences are due to the averaging process; although
the total number of observed spectra was quite large, the number of spectra available
for a given bin at a given epoch is
small, varying between 1 and 5, and furthermore the phase distribution within a
given bin is quite inhomogeneous, so that the ``average'' aspect can be strongly
influenced by the quality or phase of a single spectrum in the average.
Nevertheless, it is clear that the main features, like the phase of the maximum blue
or red peak in the central component are quite clear and reproducible.
This also confirms that the spectroscopically determined period of
Kürster & Barwig (1988) is the most appropriate for our analysis.
![]() |
Figure 12: Series of averaged profiles obtained in 2008 (binning and phases as above). |
Open with DEXTER |
![]() |
Figure 13: Long time-scale equivalent width variation of only the central component. |
Open with DEXTER |
3.3 Low-dispersion spectroscopy
The three spectra shown in Fig. 4 are rather similar: the lower
one was obtained at the OHP on 2007 Dec. 10 through thin clouds
(exposure of 10 mn), while the upper one was obtained at the OHP on
2008 Dec. 7, in good conditions (exposure of 5 mn); the middle spectrum is from
the INT on June 25, 2008 (exposure of 200 s). The extremities
of this INT spectrum (bluewards of 5000
and redwards of
7000
)
are not usable for determining the fluxes because of the
onset of a strong vignetting in the camera, which is difficult to correct
for in the reduction (the bump seen in the continuum around 5000
is also artificial). All three
spectra are calibrated, and the differences in ordinates mainly reflect the
difference in atmospheric transparency at the moment of the observations (except
that the upper spectrum has been offset upwards by 300 units for clarity).
The emission features detected comprise the Balmer lines of hydrogen
(H
,
H
,
H
,
H
), lines of neutral helium (7065,
6678, 5876, 5015, 4921
), the HeII lines at 5411 and 4686
,
as
well as the broad CIII/NIII
4640 and CIV
5805
features. The slightly higher spectral resolution of the INT (middle) spectrum
clearly shows [N II] 6548, 6584
components contributing to H
.
There is no indication of a significant variation, from one spectrum to another,
of the high-excitation lines, which could be related to the variability of the exciting source.
When comparing with the latest available spectrum in the literature, obtained
in 1990 (see Ringwald et al. 1996), we can note that the high
excitation lines (HeII, CIV) are still strong, but that the [OIII]
4959, 5007
lines, present in addition to the classical
accretion disk lines, are much weaker now, as expected
because of the expansion of the nebula. In practice, we
detect only one, faint [OIII] component, blueshifted by 500 km s-1, which
may correspond to one of the nebular spots detected by Harman & O'Brien
(2003).
The equivalent width of H
has also been measured on our
low-dispersion spectra; the comparison of the long-slit,
integrated values with the high spectral resolution values found with the 3"
circular aperture, can indicate the contribution of the nebula
to the measurements. When integrating all emission along the slit
(extension about 13'') in our 2008 December spectrum, for instance, we find a
total equivalent width of H
of
,
which is larger than any
value found on high-dispersion spectra at the same epoch.
But if we restrict the integration in the long-slit spectrum
to only the central 3'', we only find
,
in good agreement with the
high-resolution value. There is therefore at most a 20% contribution of
something else to the central value extending to large radii from the
central object.
The INT spectrum, obtained with a 1
5 slit in good seeing conditions
(
1
0) similarly yields an equivalent width for the central core of the
H
line of
,
in excellent agreement with the high-resolution
spectra core values. We can thus safely consider that most of the H
emission comes from the very central parts and that
10% is a conservative upper limit for the possible contribution of the
nebula (or the nebular blobs) to the equivalent width of the central
component of H
in the high-resolution spectra.
4 Discussion
The first point to emphasise is that the real spectroscopic phases of
HR Del are uncertain. Kürster & Barwig (1988)
estimate an error of
days for a spectroscopic period of 0.214165 days.
Their observations were obtained in 1978-1980, of the order of 9000 days before our
observations, so the phase uncertainty since then could be not far
from one orbital period. We can thus only be reasonably
certain about the relative spectroscopic phases at each epoch of our
observations, but we can seek information about
the real phases using another approach.
It should be kept in mind that Kürster & Barwig (1988)
used a convention for spectroscopic binaries, with zero phase
at maximum radial velocity, which is different from the one generally
used for cataclysmic binaries. The photometric data used the convention
for eclipsing variables, with zero phase at the minimum of the light
curve (as in Fig. 6,
the zero point being given in the
figure caption). In fact, by comparing their zero point with our
photometric zero point, we find that the Kürster and Barwig phase (NOT
the true spectroscopic phase) corresponding to our photometric zero
point, zero phase (as in our Fig. 6) is 0.42.
A particularly striking feature of our spectra is that the H
profile
seems to be a sum of the already mentioned central and wing components.
It is hard to make a comparison with the lower resolution spectra of
Kürster & Barwig (1988) with resolutions of 1 and 2
,
taken from
1978 to 1980, because those spectra still had serious interpretation problems due to
contamination by nebular emission.
We can compare the radial velocities corresponding to the widths of the
different components of our high-resolution H
profiles with the
velocities of the emission produced by the expanding nebula. Harman
& O'Brien (2003) find a maximum of about 600 km s-1 for the
radial velocity of major nebular H
emission, such a low velocity
being expected for this slow nova. Moraes & Diaz
(2009) also find that the highest radial velocity of H
emitting material is associated with the polar caps, and has a radial
velocity of 630 km s-1. If a nebular contribution were important, our H
profile
would be expected to show a component with about twice this width in
velocity units and yet the dominant central component is much narrower.
The limits in radial velocity of the line centre and line wing
components are easily measured in our spectra, as seen for instance
in Fig. 3: the violet edges are at about -250 and -650 km s-1 for the central component and the blue wing,
respectively and similar values are found on the red side (although
there they could be slightly modified owing to a possible contamination
by [N II]). We thus see that a major nebular contribution is
implausible for the
central component, but cannot be easily excluded in the wings.
Weak 6548.1 and 6583.6 Å [N II] lines are also present in the nebular emission (Harman & O'Brien 2003; Moraes & Diaz 2009). There is no evidence in our spectra (e.g. Fig. 3)
of a central peak due to these lines, that is, neither for the fainter,
6548.1 Å line nor for the line that is three times stronger at
6583.5 .
The nebular [N II]
contribution can thus be deemed undetectible. Separation of the nebular
component may be easier in the future with space-resolved spectroscopic
observations of the centre of the HR Del image.
With an inclination of 42,
a white dwarf mass of 0.65
and a
radial velocity of 250 km s-1 (measured from our spectra), the radius of the
outermost circular orbit of emission for the central component of the
profile would be at
cm or 0.89
from the centre of the
white dwarf stellar component. This is significantly larger than the radius of the
photo-ionising disk found from Moraes & Diaz (2009). The reason
for this difference is not clear at present; however, parts of the disk at
intermediate radii may not produce much ionising radiation and/or may only contribute to
the profile wings. The radius of the emitting region appears somewhat
too large when compared with the equatorial radius for an accretion disk around
the white dwarf according to the two possible system parameters of Kürster &
Barwig (1988), but the distance from the centre of the white dwarf to the
inner Lagrangian point is close to the calculated value of the outer
radius of the disk, using the expressions on page 33 of Warner (1995). The simplest
explanation of these results is that some radiation is emitted by regions of
higher radial velocity at the outer edge of the disk than is directly
deduced from the observations, but that it is occulted because of
a particular geometry of the line-emitting region.
High-velocity components have been observed in HR Del in the
past, in particular in the high velocity C IV resonance
line P Cygni absorption component in the ultraviolet, which had an edge radial velocity of
-5000 km s-1 (Selvelli & Friedjung 2003).
Such a component could have been produced in a high-velocity jet.
No evidence for it is seen today in our H
profiles nor in older optical profiles such as in
the 1990 spectrum (Ringwald et al. 1996).
One can thus wonder to what extent the ultraviolet spectrum of HR Del
has changed between the lifetime of IUE and the epoch of our optical spectra, but
no UV spectra could be obtained since the end of IUE in 1996. The only
recent UV data available are two photometric points from the Galex all-sky survey,
which give the following AB magnitudes (observations of 2006 September): 12.78 in
the near-UV (central wavelength 2315
)
and 12.96 in the far-UV (
1539
). If taken at face value, these magnitudes would indicate a small fading
in the near-UV, and a stronger one (about a factor of two) in the far-UV with respect to
the average IUE spectrum from Selvelli & Friedjung (2003). As however
HR Del is slightly brighter than the saturation limit of the Galex detectors,
those magnitudes have to be corrected for dead-time effects. When this is done according to the
precepts described in Morrissey et al. (2007), the corrected AB magnitudes become
12.36 and 12.38 in the NUV and FUV, respectively.
![]() |
Figure 14: The two Galex photometric points (lower, black points) superposed on the IUE average spectrum from Selvelli & Friedjung (2003). The horizontal dashed lines represent the bandpass of the two (NUV and FUV) filters. The upper, red points, linked to the corresponding observed (lower) values, represent the real, estimated flux after correction for deadtime effects in the detectors. |
Open with DEXTER |
This means that the NUV magnitude is still about the same as in 1988,
within the uncertainties of the correction and that in the FUV the
fading has not exceeded 30
since 1988, which is quite remarkable.
This is illustrated in Fig. 14, which is Fig. 1 of Selvelli &
Friedjung (2003), where we have plotted the Galex values over the spectrum
(uncorrected for reddening). Of course, we still need a good UV spectrum
to assess what are the real line profiles today.
In connection with non-orbital variations of H,
Selvelli & Friedjung (2003) suspected quite rapid non-orbital
variations of the C IV and N V ultraviolet resonance line
absorption in HR Del, similar to the rapid variations of
C IV seen for another old nova (V603 Aql in Friedjung et al.
1997, and in Prinja et al. 2000). Such variations
in the ultraviolet resonance absorption components of P Cygni profiles may be
expected for the winds of accretion disks of cataclysmic binaries and have been
observed (Froning 2005; Proga 2005).
The ``S wave type'' pattern of H,
displayed in
Figs. 10 to 12, is quite clear on the blue side, but somewhat less so
on the red side of the profiles. It does, however, not appear to be a classical
cataclysmic binary S wave. However, since it moves over the
central component of the line profile like the classical wave, it can be
understood as near the outer edge of a rotating disk in the
neighbourhood of a fixed point in the rotating frame of the binary.
The profile shows a maximum assymmetry towards the violet around
Kürster and Barwig phase 0.75 or at photometric phase 0.33, which
is only 0.08 of an orbital period after photometric phase
0.25. If we interpret the photometic variations as due to varying
visibility of the regions of the surface of the mass loser,
heated by radiation from the disk and white dwarf (see below),
the photometric zero phase is when the mass loser is nearest
the observer in inferior conjunction. The phase of maximum
assymetry of the profile around photometric phase 0.33 is then not
far from quadrature when the spot is near a line perpendicular to the
line of sight. An accretion disk should rotate in the same direction
as the binary (conservation of angular momentum), so the H
bright spot should be on the side of the disk opposite to the mass
loser. Such a situation is not classical. That could possibly be
related to the large size of the accretion disk found above.
We can emphasise however that the stability of the S wave over the six years of our spectroscopic observations (about 10 000 orbital cycles) indicates that it moves with a period which is quite close to the orbital period of Kürster and Barwig. Indeed, using their estimated error for the period, we calculate that the phase shift over the time of our spectroscopic observations would be of the order of 0.23. Positive and negative superhump periods, which can be permanent, would, on the other hand, differ from the orbital period by a factor of a few percent when the orbital period is that of HR Del (Olech et al. 2009).
The variations in the V magnitude cannot be caused by the aspect changes of
the distorted Roch lobe-filling cool component in a normal state. The visual
magnitude was around 12.5 during our observations, which gives an absolute
magnitude around 2.35, using a distance of 850 pc and a visual extinction
of
(Selvelli & Friedjung 2003). An ordinary
star not heated by its companion, with a mass of 0.53 solar masses
(mass of the companion, according to Selvelli & Friedjung
2003) should have an absolute magnitude of the order of 8.5 (assuming the star is on
the main sequence). The contribution of an unheated
cool component to the V magnitude would thus be around 1/300,
and variations in its visibility could not produce the observed photometric changes in V
with a full amplitude of 0.08 mag. It would moreover be hard to
produce the observed symmetrical light curve, by supposing extra light
produced by an assymmetrical bright spot.
If the ellipticity (or, more precisely, the shape of the Roche lobe) of
the cool component produced the photometric variations, this component would
then be strongly heated by the hot component, consisting of a white
dwarf plus an accretion disk. As we have seen, this interpretation poses
a problem in understanding the phasing of the H
bright spot.
Such an interpretation disagrees with the relative
phases of the S wave and the photometric variations unless the S wave
is not produced by a classical bright spot. A more detailed discussion
of this aspect is beyond the scope of this paper and will be
given in a forthcoming one (Voloshina et al., in prep.).
![]() |
Figure 15:
Combined longer time-scale photometric and spectroscopic behaviour. The upper part shows the V magnitude (corrected for orbital variations) versus Julian date,
and the lower part the total equivalent width of H |
Open with DEXTER |
The large ``activity'' of HR Del found by Selvelli & Friedjung (2003) and confirmed by the high luminosty determined by Moraes & Diaz (2009) still appears to be unusual, regardless of the explanation given for it. The high luminosity leads to heating of the mass loser.
Long time-scale variations are apparent both in the total magnitude of the object and in the equivalent width of H
(see Fig. 15),
but the lack of obvious relation between the two does not give us any
clues as to their origin. Those long-time scale variations cannot be
similar to those of dwarf novae, which, when not in outburst, are
fainter than recent old
novae and not brighter like HR Del. The behaviour of dwarf
novae is
understood as the result of a lower accretion rate of the white dwarf
component outside outburst,
compared with that of old novae. Present theories of dwarf nova
outbursts explain them by
accretion disk instabilities, which can occur at low accretion rates.
Small
outbursts, somewhat similar to those of dwarf novae, have been observed
for the classical novae GK Per and V446 Her some
decades
after the outburst. This is understood as caused by a drop of the
accretion rate at such times, associated with the appearence of
instability. The situation is less clear for
the low-amplitude (less than 1 mag), ``stunted'' outbursts seen
for
other old novae and novalike systems (Honeycutt et al. 1998).
Solar type cycles of the mass loser, leading to variations in the accretion rate on time scales of about a decade, have also been suggested as an explanation of long-term changes in the luminosity of cataclysmic variables (Bianchini 1987; Warner 1988). It remains to be seen whether this is possible when there is strong irradiation of the mass loser, if one wishes to explain the long timescale variations of HR Del in this way.
5 Conclusions
We have performed photometric and high-resolution H
spectroscopic
observations of HR Del from 2002 to 2008. The photometry shows
variations depending on orbital phase, plus apparently random variations
from night to night and perhaps variations on a longer-time scale of
several years. The orbital phase-dependent variations may be interpreted by irradiation
of the mass losing companion star by a still very bright, expanded white dwarf and its accretion
disk. It is not clear at present whether the long time scale variations are
interpretable in terms of solar type cycles of the mass loser.
Ultraviolet, broad-band photometric data from the GALEX archive suggest only a weak fading since the epochs of the IUE, but no recent spectra are available to estimate the present activity from the UV spectral lines.
The H
profile shows both a narrow and a wide component. The
former appears to be due to emission by an accretion disk, while the latter
probably contains a non-negligible contribution from the expanding nebula. The central
component does not show a very clear dependence of
equivalent width with orbital phase. Its profile, however, shows a clear
orbital variation, which looks similar to the S wave generally seen in spectra of
cataclysmic binaries, but which has the wrong phase dependence to be
produced by a classical bright spot. The stability of the S wave with respect to
the orbital period over 8 years indicates no connection with the period of
any positive or negative superhump. The exact nature of our ``S wave''
still remains to be elucidated.
In view of the exceptional brightness and activity of this old nova about 35 years after its eruption, continuing low-level thermo-nuclear combustion in the envelope of a low-mass white dwarf may be considered a possible explanation. The evolution of this interesting stellar system should therefore be closely followed over a large wavelength range in the future.
AcknowledgementsI. Voloshina acknowledges support from the CNRS-EARA, which sponsored her visits to the IAP, where part of this work was completed. M. Dennefeld and M. Friedjung thank D. Briot who kindly cooperated in obtaining a long series of spectra in August 2008, and the OHP night assistants (most notably R. Giraud and J.C. Mevolhon) for their valuable contribution to the service observations. The authors are grateful to an anonymous referee for useful suggestions.
Galex is a NASA Small Explorer, developed in cooperation with the French space agency CNES and the Korean Ministry of Science and Technology.
References
- Baldwin, J. A., Philipps, M. M., & Terlevich, R. 1981, PASP, 93, 5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bianchini, A. 1987, Mem. Soc. Astr. It., 58, 245 [Google Scholar]
- Candy, M. P. 1967, IAU Circ., 2022 [Google Scholar]
- Deeming, T. J. 1975, Ap&SS, 36, 137 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
- Friedjung, M. 1992, A&A, 262, 487 [NASA ADS] [Google Scholar]
- Friedjung, M., Selvelli, P., & Cassatella, A. 1997, A&A, 318, 204 [NASA ADS] [Google Scholar]
- Friedjung, M., Bruch, A., Voloshina, I., et al. 2005, in The Astrophysics of Cataclysmic Variables and Related Objects, ed. J. M. Hameury, & J. P. Lasota, ASP Conf. Ser., 330, 437 [Google Scholar]
- Froning, C. S. 2005, in The Astrophysics of Cataclysmic Variables and Related Objects, ed. J. M. Hameury, & J. P., Lasota, ASP Conf. Ser., 330, 81 [Google Scholar]
- Gillet, D., Burnage, R., Kohler, D., et al. 1994, A&AS, 108, 181 [Google Scholar]
- Harman, D. J., & O'Brien, T. J. 2003, MNRAS, 344, 1219 [NASA ADS] [CrossRef] [Google Scholar]
- Honeycutt, R. K., Robertson, J. W., & Turner, G. W. 1998, AJ, 115, 2527 [NASA ADS] [CrossRef] [Google Scholar]
- Hutchings, J. B. 1980, PASP, 92, 458 [NASA ADS] [CrossRef] [Google Scholar]
- Kürster, M., & Barwig, H. 1988, A&A, 199, 201 [NASA ADS] [Google Scholar]
- Lamers, H. J. G. L. M., Snow, T. P., & Lindholm, D. M. 1995, ApJ, 455, 269 [NASA ADS] [CrossRef] [Google Scholar]
- Lenz, P., & Breger, M. 2004, IAUS, 224, 786 [Google Scholar]
- Moraes, M., & Diaz, M. 2009, AJ, 138, 1541 [NASA ADS] [CrossRef] [Google Scholar]
- Morrissey, M. P., Conrow, T., Barlow, T. A., et al. 2007, ApJS, 173, 682 [NASA ADS] [CrossRef] [Google Scholar]
- Ness, J. U., Schwarz, G., Starrfiels, S., et al. 2008, ApJ, 135, 1228 [Google Scholar]
- Olech, A., Rutowski, A., & Schwarzenberg-Czerny, A. 2009, MNRAS, 399, 465 [NASA ADS] [CrossRef] [Google Scholar]
- Prinja, R. K., Knigge, C., Ringwald, F. A., & Wade, R. A. 2000, MNRAS, 318, 368 [Google Scholar]
- Proga, D. 2005, in The Astrophysics of Cataclysmic Variables and Related Objects, ed. J. M. Hameury, & J. P. Lasota, ASP Conf. Ser., 330, 103 [Google Scholar]
- Ringwald, F. A., Naylor, T., & Mukai, K. 1996, MNRAS, 281, 192 [NASA ADS] [Google Scholar]
- Sanyal, A. 1974, ApJS, 28, 115 [NASA ADS] [CrossRef] [Google Scholar]
- Selvelli, P. L., & Friedjung, M. 2003, A&A, 401, 297 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sion, E. M., & Starrfield, S. G. 1994, ApJ, 421, 261 [NASA ADS] [CrossRef] [Google Scholar]
- Walraven, Th., & Walraven, J. H. 1972, in Auxiliary Instrumentation for Large Telescope, ESO/CERN conference, ed. S. Laustsen, & A. Reiz, 175 [Google Scholar]
- Warner, B. 1988, Nature, 336, 129 [NASA ADS] [CrossRef] [Google Scholar]
- Warner, B. 1995, Cataclysmic Variable Stars (Cambridge University Press) [Google Scholar]
Online Material
Table 1: Journal of the CCD observations.
Table 2: Details of spectroscopic observations.
Table 3: Details of spectroscopic observations, continued.
Footnotes
- ... Delphini
- A journal of the CCD photometric observations is presented in Table 1. Observing data for individual spectra are presented in Tables 2 and 3. All these tables are only available in electronic form at http://www.aanda.org
- ...
- Visiting Astronomer, Observatoire de Haute-Provence, CNRS, France.
- ...
- Visiting Astronomer, Crimean station of the Sternberg Astronomical Observatory.
All Tables
Table 1: Journal of the CCD observations.
Table 2: Details of spectroscopic observations.
Table 3: Details of spectroscopic observations, continued.
All Figures
![]() |
Figure 1: Overall behaviour of HR Del in the V band during 2002-2008. Photoelectric V magnitudes are shown versus Julian Day. The arrows indicate the epochs of the spectroscopic observations. |
Open with DEXTER | |
In the text |
![]() |
Figure 2: Examples of CCD observations for different dates: CCD-V magnitudes versus phase. |
Open with DEXTER | |
In the text |
![]() |
Figure 3:
Example of a high dispersion spectrum, obtained on November 1, 2003 at the OHP. The core and the low-intensity wings of H |
Open with DEXTER | |
In the text |
![]() |
Figure 4:
Low-dispersion spectra. Lower part: OHP, Dec. 2007; Middle: INT, June 2008. Upper part: OHP, Dec. 2008. Small offsets were made between spectra (see text). The units in the ordinates are 10-16 erg/cm2/s/ |
Open with DEXTER | |
In the text |
![]() |
Figure 5: Power spectrum of the 2008 CCD observations. |
Open with DEXTER | |
In the text |
![]() |
Figure 6: Light curve of HR Del versus photometric phase, in the V band, from observations obtained during 2005-2006. The averaged orbital sine wave is superposed onto the observed data. |
Open with DEXTER | |
In the text |
![]() |
Figure 7: Equivalent width variation with orbital phase in 2003 November. Filled circles (above) represent the total equivalent width, triangles (middle) the central component and the squares (lower) the wings alone. |
Open with DEXTER | |
In the text |
![]() |
Figure 8: Equivalent width variation in 2005 mid September. Symbols as above. |
Open with DEXTER | |
In the text |
![]() |
Figure 9: Equivalent width variation of the central component with phase, for four different observing periods. Triangles are for 2003 November, crosses for 2005 September, open stars for 2005 October and filled polygons for 2008 August. |
Open with DEXTER | |
In the text |
![]() |
Figure 10: Series of profiles obtained during 2002-2004, averaged in ten bins of 0.1 in phases. The phases are going from 0 at the bottom (bin 0 to 0.1) to 1 at the top (bin 0.9 to 1. |
Open with DEXTER | |
In the text |
![]() |
Figure 11: Series of averaged profiles obtained in 2005 (binning and phases as above). |
Open with DEXTER | |
In the text |
![]() |
Figure 12: Series of averaged profiles obtained in 2008 (binning and phases as above). |
Open with DEXTER | |
In the text |
![]() |
Figure 13: Long time-scale equivalent width variation of only the central component. |
Open with DEXTER | |
In the text |
![]() |
Figure 14: The two Galex photometric points (lower, black points) superposed on the IUE average spectrum from Selvelli & Friedjung (2003). The horizontal dashed lines represent the bandpass of the two (NUV and FUV) filters. The upper, red points, linked to the corresponding observed (lower) values, represent the real, estimated flux after correction for deadtime effects in the detectors. |
Open with DEXTER | |
In the text |
![]() |
Figure 15:
Combined longer time-scale photometric and spectroscopic behaviour. The upper part shows the V magnitude (corrected for orbital variations) versus Julian date,
and the lower part the total equivalent width of H |
Open with DEXTER | |
In the text |
Copyright ESO 2010
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.