A&A 487, 611-618 (2008)
DOI: 10.1051/0004-6361:20078644
C. Papadaki1,2 - H. M. J. Boffin3 - D. Steeghs4 - L. Schmidtobreick1
1 - ESO, Casilla 19001, Santiago19, Chile
2 - Vrije Universiteit Brussel, Pleinlaan 2, 1050 Brussels, Belgium
3 - ESO, Karl-Schwarzschild-Str. 2,
85738 Garching, Germany
4 - Department of Physics, University of Warwick, UK
Received 10 September 2007 / Accepted 1 April 2008
Abstract
Aims. We performed a detailed spectroscopic analysis of the dwarf nova V2051 Oph at the end of its 1999 superoutburst. We studied and interpreted the simultaneous behaviour of various emission lines.
Methods. We obtained high-resolution echelle spectroscopic data at ESO's NTT with EMMI, covering the spectral range of 4000-7500 .
The analysis was performed using standard IRAF tools. The indirect imaging technique of Doppler tomography was applied, in order to map the accretion disc and distinguish between the different emission sources.
Results. The spectra are characterised by strong Balmer emission, together with lines of He I and the iron triplet Fe II 42. All lines are double-peaked, but the blue-to-red peak strength and central absorption depth vary. The primary's velocity was found to be 84.9
.
The spectrograms of the emission lines reveal the prograde rotation of a disc-like emitting region and, for the Balmer and He I lines, an enhancement of the red-wing during eclipse indicates a bright spot origin. The modulation of the double-peak separation shows a highly asymmetric disc with non-uniform emissivity. This is confirmed by the Doppler maps, which apart from the disc and bright spot emission also indicate an additional region of enhanced emission in the 4
quadrant (+Vx, - Vy), which we associate with the superhump light source. Given the behaviour of the iron triplet and its distinct differences from the rest of the lines, we attribute its existence to an extended gas region above the disc. Its origin can be explained through the fluorescence mechanism.
Key words: accretion, accretion disks - novae, cataclysmic variables
The SU UMa stars are cataclysmic variables (CVs) belonging to the sub-category of dwarf novae (DN). Like all CVs they consist of a red dwarf (secondary) that fills its Roche lobe and orbits around a white dwarf (primary). The primary accretes matter from the secondary through an accretion disc (AD). The collision between the accretion stream and the rim of the AD results in the dissipation of energy and the formation of a ``bright spot'' (BS). The DN are characterised by two brightness states, quiescence and outburst of typically 2-5 mag. The main characteristics of the SU UMa stars are that they undergo superoutbursts, where the system's brightness increases even more than in a normal outburst, and have orbital periods below or, in some cases, in the CV period gap of 2-3 h. For an extensive review of DN and SU UMa stars, in particular, see Warner (1995), Osaki (2001) and references therein.
![]() |
Figure 1: V2051 Oph average out-of-eclipse spectrum. All lines are double-peaked with prominent emission from the Balmer and He I lines. He II and Fe II lines are also in emission, albeit weaker. |
Open with DEXTER |
V2051 Oph (RA 17:08:19.09, Dec -25:48:30.8, J2000) was discovered by Sanduleak (1972) and has ever since been an object of great interest and continuous study. Although a lot of questions have been raised on its classification as a low-field polar (Warner & O'Donoghue 1987) or a DN, it is now widely believed to belong to the SU UMa subtype of DN. It has a short period of 1.5 h (Warner & Cropper 1983) and shows recurrent but sparse outbursts (Warner & Cropper 1983; Warner & O'Donoghue 1987) and superoutbursts (Kato et al. 2001). Its visual magnitude is 15 mag in quiescence (Warner & Cropper 1983), 13.2-14.6 mag in outburst and 11.6-12.9 in superoutburst (Kato et al. 2001). A positive superhump, i.e. a modulation originating from the beating between the prograde precession of the AD and the orbital motion of the binary, has also been detected by Kiyota & Kato (1998) and confirmed by Patterson et al. (2003). The light curves of the latter were obtained by G. Garradd
and show a clear hump before each eclipse, indicating the presence of a strong BS. The system's high inclination
and therefore deep eclipses, along with follow-up observations during an unusually low-brightness state (
), allowed Baptista et al. (1998) to accurately define the binary parameters.
Compared to other SU UMa systems, V2051 Oph displays a variety of eclipse morphologies and shows humps at various phases and irregular occurrences (Warner & Cropper 1983; Vrielmann et al. 2002; Warner & O'Donoghue 1987). This behaviour was the cornerstone of the polar scenario, especially since no superoutbursts had been observed up to that point. Another characteristic of the system is the large amplitude flickering activity during quiescence. It reaches 30% of the total light in the optical and occurs at all phases apart from eclipse (Cook & Brunt 1983). The flickering source was studied through eclipse mapping by Baptista & Bortoletto (2004), who found two different dominant contributions: (1) low-frequency flickering, mainly arising from an overflowing gas-stream and (2) high-frequency flickering, produced in the AD.
In the first spectroscopic studies, V2051 Oph clearly showed Balmer, He I and He II lines in emission (Williams 1983). Profile distortion during eclipse, characteristic of a prograde rotation of a disc-like emission line region (Cook & Brunt 1983), was also evident. A study during quiescence, covering 4100-5000 Å (Watts et al. 1986), showed the blue peak to be stronger than the red, in all emission lines, and the peak separation in H
to be larger than that in H
,
pointing to a more central location of the former. UV spectra obtained from the same study were similar to that of DN in quiescence. Warner & O'Donoghue (1987) obtained the first spectra during decline from an outburst (1 mag higher than quiescence), and found the separation and widths of the lines to be similar to the ones published by Watts et al. (1986) during quiescence. However, there was a distinct difference between the set of spectra during the two states: the V/R variations (flux ratio of the blue to the red wings) around the orbit were in almost exact anti-phase, indicating a large change in the distribution of the line-emitting regions (Warner & O'Donoghue 1987). It has to be noted that in both cases no evidence for BS emission nor the resulting S-wave was found. The first Doppler maps of the H
and He I
4471 Å lines, according to the back-projection technique, are presented in Kaitchuck et al. (1994). The trailed spectra of the lines confirmed the double-peaked emission and rotational disturbance, while the maps showed the ring-shaped emission pattern of the disc. Steeghs et al. (2001) were the first to detect and resolve DN oscillations in the emission lines. Coherent oscillations were also detected in the continuum. The same study showed that the disc is highly asymmetric with a stronger contribution from its blue-shifted side. Recently, Saito & Baptista (2006) reported a study of the spectra and the AD while the system was in an unusual low brightness state during 1996. Among others, they proposed a model of a vertically extended disc with the line emission arising from the upper atmospheric layers, gas-stream overflow, an extended gas region above the disc and a non-negligible opening angle of the disc.
The aim of this work is to accomplish a detailed spectroscopic study of V2051 Oph. Albeit as mentioned above, V2051 Oph shows some rather unique behaviour, it is also a very useful representative of the class of dwarf novae, and as such, serves to understand the physics at play in these systems better. In particular, one can investigate the possible presence of stream-disc overflow or of matter outside of the disc's plane. Moreover, our unique set of spectra, with its broad wavelength coverage and high resolution gives us the opportunity to simultaneously study the behaviour of the various emission lines and, through the indirect imaging technique of Doppler tomography, to map for the first time the emission distribution of the AD in detail. This study thus also serves to show the merit of echelle spectroscopy-based Doppler tomography.
This paper is organised as follows: in Sect. 2 we present our observations and data reduction technique, in Sect. 3 we perform the spectral analysis and Sect. 4 deals with the application and results of Doppler tomography. Last, Sect. 5 is devoted to the discussion of our two principal results, while Sect. 6 gives a summary.
The spectra were reduced following the standard reduction procedures in IRAF. The order extraction, wavelength calibration, flux calibration and order merging were performed by the echelle package and the reduction task doecslit in IRAF. The spectrophotometric standard LTT 7379 (Hamuy et al. 1994,1992) was observed during the same night and served for the flux calibration. The spectra were also corrected for cosmic rays.
V2051 Oph was observed 2 days after the end of the 1999 superoutburst (Kato et al. 2001), when its brightness rapidly fell to approximately that of quiescence. Nevertheless, photometry conducted by Patterson et al. (2003) during the same night, revealed clear superhumps, showing therefore that the disc was still large enough to precess and had not fully recovered from outburst.
![]() |
Figure 2: V2051 Oph trailed spectra of the 11 most prominent lines. The spectra have been ordered according to phase. |
Open with DEXTER |
The average out-of-eclipse spectrum is shown in Fig. 1. All lines are in emission and appear double-peaked. The prominent ones are marked on the same figure. The spectrum is characterised by strong Balmer emission on a slightly blue continuum, together with lines of He I and a weak line of He II
.
Lines of the heavier element, Fe II, are also present. All Balmer and He I lines are asymmetric and show a stronger blue peak, although in the Balmer lines the difference in strength is more pronounced. Only in the case of He I
5016 is the red peak stronger but that line is blended with Fe II. As for the unblended Fe II (5169 Å) line, it shows a stronger red-peak. The Balmer decrement from H
to H
is 1:0.67:0.56:0.43, showing a mix of optically thick and thin emission. We also searched for orbital variations of the Balmer decrement, but none were evident. Table 1 shows the full width at zero intensity (FWZI), equivalent width (EW) and integrated flux of the 11 most prominent lines on which we will focus.
The three Fe II lines (4924,
5018 and
5169) compose the iron triplet Fe II 42. Fe II
4924 and
5018 are blended with He I lines, which are fairly strong in our spectra. However, the existence of the Fe II
5169, the deeper central absorption of all Fe II lines (compared to the rest) and the inversion of blue-to-red peak strength in all Fe II lines (except the one at
4924) make us conclude that Fe II has a non-negligible contribution in the blended lines. The only case where Fe II could be less dominant is in the blend with He I
4922. A fourth Fe II emission line is seen at
5317 but is too weak to be included in further analysis.
A red-blue wing asymmetry is observed in the lines on the blue part of the spectrum, in Fig. 1. To be sure that spectra belonging to incomplete orbits are not responsible for the asymmetries, the average spectra of the individual full orbits were calculated as well, but the pattern persisted. We also checked if the merging of the echelle orders could be responsible, but no significant order mismatches were evident. However, the order merging is responsible for the higher noise of the red wing of H
(Figs. 1, 2). The asymmetry could therefore originate either from blended lines or from an enhancement of the continuum. The first scenario could possibly be the cause of the asymmetry for H
and H
,
owing to the proximity of the He I lines at 4121 and 4388 Å, respectively. However, the blend of He I
5016 with Fe II
5018, which shows the same red-wing asymmetry, has no emission line in the proximity interfering with the red wing. The same applies to He I
4471. For longer wavelengths the lines appear symmetric. The same behaviour was indicated by Watts et al. (1986) and Steeghs et al. (2001), giving further credit to its existence.
The trailed spectra, calculated according to the following orbital ephemeris (Baptista et al. 2003),
![]() |
(1) |
The diagnostic diagrams, encompassing the double Gaussian fitting technique described by Shafter (1983), were used to estimate the projected primary velocity K1 and the systemic velocity .
When using the information from all emission lines, K1 was found to be
km s-1. However, if we avoid the weaker lines which give higher discrepancies in the diagnostic diagrams, we obtain a refined value of
km s-1. This value is in good agreement with the one of
km s-1 from Baptista et al. (1998) and 91 km s-1 from Watts et al. (1986). The first one was predicted by a photometric model based on measurements of contact phases in the eclipse light curve, while the second one was derived from the radial velocity curves of H
and H
.
As for
,
even though it varied around 0 for most of the lines, it had greater values for some of them. Discrepancies of
values between the emission lines of V2051 Oph, as well as other CVs, have been reported before (see Watts et al. 1986, and references therein). This anomaly is associated with blendings with neighbouring lines, red-blue wing asymmetries and any deviations from the assumed flatness of the underlying continuum. Any of these reasons can shift
from its correct value and our spectra show that most of them apply in our case. The results from the diagnostic diagrams are given in Table 2. It also includes the phase offset
and the separation of the Gaussians d in Å. Only He I
5876 has not been included in the Table because the diagnostic diagrams showed inconclusive results.
Even though the method developed by Shafter (1983) gives the primary's velocity in close approximation to its real value, Marsh (1988) warned against the risk of biases and the fact that the value obtained might depend on the signal-to-noise ratio. Marsh proposed a different method that is less prone to bias and that takes the increased sensitivity into account. The application of this method to all the emission lines from our data set gives
.
If the weaker lines are not included, as in our previous K1 calculation, we find
.
This last value is in excellent agreement with the one predicted by Baptista et al. (1998). The weights were calculated from the equivalent width of each line, as shown in Table 1.
Taking into account each line's ,
we constructed the spectrograms of Fig. 3. All spectra were folded on phases 0-1 and repeated over two cycles. In order to improve and enhance the features of the spectrograms, the spectra were binned by 4 pixels which gave optimal results. All lines show the characteristic double-peaked structure, as well as the blue wing being eclipsed prior to the red one, characteristic of a prograde rotation of a disc-like emission region. The asymmetry of the emission lines is evident, in some cases though being more prominent than others. In most of the spectrograms, there is a brightness enhancement of the red wing during eclipse, probably indicating a BS origin.
Moreover, as seen from the spectrograms but even more clearly from Fig. 4, the double-peak separation is highly phase dependent. This is a strong signature that the disc is highly asymmetric with a non-uniform emissivity. Through period analysis, it has been found that a periodicity of
best describes the variation of the double-peak separation. This applies to all lines except Fe II
5169, probably due to the high scatter owing to its low strength. Such a variation points to the existence of two emission sources on the AD. The indirect imaging technique of Doppler tomography, as presented in the next section, will address this issue.
Table 1: Characteristics of emission lines.
Table 2: V2051 Oph emission line wing velocity parameters.
![]() |
Figure 3: V2051 Oph trailed spectrograms. The x-axis represents the velocity in km s-1 and the y-axis the orbital phase. The grayscale contrast has been adjusted to emphasise between emission (black) and absorption (white). |
Open with DEXTER |
Assuming all line emitting regions are located on the AD (however this might not be true for all lines and will be discussed in the next sections) and that the velocity field of the AD is Keplerian then the orbital radius
corresponding to a velocity
is given by the equation,
.
Given the fact that the double-peak separation reflects the velocity of the gas on the outer edge of the AD and that the extension of the emission line wings reflects that of the inner disc, then the outer and inner radii of the AD are given by the equations:
![]() |
(2) |
![]() |
(3) |
The indirect imaging technique of Doppler tomography, first introduced by Marsh & Horne (1988), translates orbital phase-dependent line profiles into a two dimensional velocity space map. In this way, the AD emission distribution (ring shaped structure and/or spiral patterns) and any isolated emission sources (secondary irradiation, BS & stream trajectory) can be readily identified. For our set of data, the Doppler maps have been produced using the implementation code of Spruit (1998), which applies the Maximum Entropy Method (MEM).
The resulting Doppler maps of the 11 selected emission lines are shown in Fig. 6. Both binned and unbinned spectra were used as an input but there was no significant change in the resulting doppler maps. Therefore, the unbinned spectra which gave maps of higher resolution, were adopted. The first row of Fig. 6 shows the Balmer series, the second one the He I lines and the third one the Fe II triplet. In all maps the x- and y-axis represents the Vx- and Vy-velocities, in the range of -1500 to 1500
.
The maps of each row have the same intensity scale, while each row has a different one. This was done in order to both highlight the features seen and be able to directly compare lines of the same series.
The Balmer and He I maps show the ring-shaped emission characteristic of an AD - the result of the double-peaked emission lines, in accordance to the two maps (H
and He I
4471) presented by Kaitchuck et al. (1994). However, the higher resolution maps presented in this work have much more to reveal. They show BS emission, albeit of different strength in each case. This emission is located at the lower part of the second quadrant (
-V x, + V y), and was expected since the simultaneous photometric light curve, mentioned in the introduction, shows a prominent hump caused by the BS. The existence of a BS was also seen in the spectrograms of Fig. 3, where there is a very bright red-shifted emission during the eclipse. This occurs when the BS emission is most obvious since most of the disc is eclipsed but the BS which precedes the secondary is fully visible. What is also evident from both Balmer and He I lines, is that the strength of the BS increases when moving to shorter wavelengths. It is most prominent in H
and fades away when moving to lower series. This can be interpreted in terms of optical thickness. The BS appears to be optically thick in H
,
compared to the optically thin disc-emission and progressively less in the rest of the Balmer series. The same is valid for the He I lines. This behaviour has been observed before in other DN as shown for VW Hyi in Smith et al. (2006) and references therein. The map of He I
4471 is not of the same quality as the rest due to its proximity to the continuum, but has been included since it still shows the same features.
![]() |
Figure 4:
Modulation of the double-peaks separation. Superimposed is the dominant periodicity of
![]() |
Open with DEXTER |
But, the most striking feature of all, is the one seen in the
quadrant (
+ V x, - V y) of the maps, with a surprising persistence in all Balmer and He lines. Such a feature does not coincide with the location of a spiral arm, always seen in the first (
+ V x, + V y) and third quadrant (
- V x, -Vy). The existence and persistence of such a rare feature adds a great deal of interest to its origin and so will be discussed in Sect. 5.2. However, a noticeable difference in the feature's location on the different maps has to be mentioned. This is mostly seen in the Balmer series but not in the He lines. A possible trend of the feature moving downwards on the map with increasing wavelength, is visible. Its opening angle seems to differ as well, being greater in H
compared to H
.
The Fe II triplet shows a different behaviour than the rest of the lines. The BS is absent in the 5169 line, and most probably in the other two lines as well. The only reason for a more enhanced emission at the BS region could be the blending with the He I lines. The
quadrant emission seen in the rest of the lines is absent. The difference of these maps compared to the rest is further discussed in Sect. 5.1.
![]() |
Figure 5: Radial extension of all emission lines. Dashed lines correspond to the refined outer radii when the tidal effect is taken into account. |
Open with DEXTER |
![]() |
Figure 6:
V2051 Oph Doppler maps. Top to bottom rows correspond to the Balmer lines, He I lines and the Fe II triplet, respectively. The x and y axis represent the Vx and Vy velocities respectively, spanning from -1500 to 1500
![]() |
Open with DEXTER |
We think that the Fe II emission can be adequately explained as being produced through the fluorescence mechanism, in which UV photons are emitted by a hot source and then absorbed by iron atoms. The absorption of Fe II in the UV triggers its emission in the optical. Absorption bands in the UV, due to Fe II have already been observed for V2051 Oph by Saito & Baptista (2006) and this iron curtain could simply be the UV counterpart of the optical emission. Saito & Baptista (2006) also suggested that this could originate from absorption by an extended gas region above the disc. The same explanation has already been given by Mason & Howell (2005) for the case of a new short orbital period CV, H0242-2802. They inferred an orbital period of 107 min and an inclination of
,
both parameters being similar to those of V2051 Oph. However, the lack of UV spectroscopy has not yet confirmed their interpretation. It has to be noticed, that the Fe II emission lines in their spectra, covering 4200-5600 Å, also showed deeper central absorption than the Balmer lines and a stronger red than blue peak in contrast to the Balmer lines. This is in accordance with what we have observed in our spectra. In their case He I emission was weak and the Fe II 42 observed was considered to be mostly Fe II emission.
Furthermore, V2051 Oph has not shown any Fe II emission during quiescence (Williams 1983) or in decline from a normal outburst (Steeghs et al. 2001). In consequence, its emission is likely related to the shortly terminated superoutburst, whose increased temperature may have excited Fe II in the gas region. The same has also been observed for WZ Sge, which showed no emission during quiescence but after an outburst (for more details see Mason & Howell 2005, and references therein).
Given all the aforementioned facts, the most plausible explanation is that Fe II does not come from the same line-forming region as the Balmer and He I lines, therefore does not form on the AD itself but from an extended gas region above it.
In the Doppler maps of Balmer and He I lines we have seen an arc-shaped region of enhanced emission located on the
quadrant. It should be noted that the existence of two locations of prominent emission in the maps was already evident from the modulation of the double-peaks separation (Fig. 4). Since one of them is the BS and the double-peak separation is modulated with half the orbital period, that means that the second region of emission should be opposite the BS, as seen in the maps.
Such emission has also been reported by Smith et al. (2006) for another SU UMa star, VW Hyi, just prior to the onset of an outburst. However in this case the emission was only seen in the Doppler map of H
and remained absent in the rest of the Balmer and He I lines (they shared a similar spectral coverage to ours). The appearance of this feature in H
maps alone, and their limited phase coverage, left them unconvinced of its existence.
We recall that we have captured the system shortly (2 d) after the end of a superoutburst and that it showed photometric superhumps at the time of our observations (Patterson et al. 2003). It is then quite natural to assume that the enhanced emission region seen in our maps is nothing else but the equivalent of the superhump light source (SLS). Concerning its shape, for the majority of systems, the SLS is an extended region of the disc involving 1/4 of its area (Warner 1995, p. 199), which agrees with our maps. Warner & O'Donoghue (1987) saw a narrow hump in their light curves during outburst which they associated with a possible superhump. They gave phase 0.54 for its emission peak. Such an emission should be seen at the lower half of the
quadrant in our maps. In H
,
H
and He I
5876 this could be true, while it does not appear to hold for H
,
H
,
or He I
4471, for instance. However, in any cases, caution should be drawn to such an interpretation, since Warner & O'Donoghue (1987) did not fully trust the reality of their signal due to rapid variations of brightness at that time.
Additionally, as pointed in Sect. 1 the system at that time had a different emission distribution in its disc. The V/R variations of their spectra were in almost anti-phase with ours and their average spectrum showed a stronger blue-to-red peak for H
in contrast to ours.
Still, the physical origin of the SLS for CVs is unclear. As briefly mentioned before, this emission's origin most probably should not be attributed to stream-disc overflow. SPH simulations by Kunze et al. (2001) offer additional support to this opinion. In their study they used, among others, the parameters of the SU UMa star OY Car, which is considered the ``twin'' of V2051 Oph. Their 3D simulations showed that the re-entry region of the overflowing stream onto the disc is located close to zero velocity, at orbital phase 0.5. Such a region of enhanced energy is best pictured by Doppler tomography and in our maps should be expected on the border of the
and
quadrants. However, the enhanced emission region in all our maps (except perhaps for He I
5876) does not agree with the one predicted by the numerical simulations (see e.g. references in Kunze et al. 2001). We caution however for the fact that our data set gives us only a resolution in the maps of about 33
.
It is thus difficult to determine the boundary between quadrants precisely and reach a definitive conclusion.
However, if this feature is indeed associated with the SLS, it is rather puzzling that it was seen at the onset of the outburst of VW Hyi, when the disc was not large enough for superhumps to develop.
Another possibility is that the arc-shaped emission is due to tidal interactions of the secondary and the disc, like the two spiral arms observed in DN in outburst. But this has to be ruled out since simulations predict and observations verify that the two spirals are seen in the
and
quadrants, which is not the case.
Acknowledgements
We kindly thank C. Tappert for providing the MIDAS adapted version of the Doppler tomography and diagnostic diagrams codes. The authors wish to thank the anonymous referee for pertinent remarks on the manuscript. C.P. wishes to thank A. Bonanos for useful comments on the manuscript and gratefully acknowledges a doctoral research grant by the Belgian Federal Science Policy Office (Belspo).