© ESO, 2014

1. Introduction

The massive binary nature of HD 47129 was first reported by Plaskett (1922), who inferred a total mass of 138.9 M_⊙, the highest ever observed at that time. He described the binary as consisting of two early-type stars: an Oe5 primary and a fainter secondary with weak and broad lines. For almost a century, HD 47129, also known as Plaskett’s star, has been the target of many studies, with tremendous progress in the past twenty years (e.g. Bagnuolo et al. 1992; Wiggs & Gies 1992; Linder et al. 2006, 2008; Mahy et al. 2011; Grunhut et al. 2013a). The masses of the components have been revised downwards to m₁ sin³i = 45.4 for the primary and m₂ sin³i = 47.3M_⊙ for the secondary, leading to a mass ratio of about 1.05 ± 0.05 (Linder et al. 2008). Unfortunately, Plaskett’s star is a non-eclipsing binary, preventing a precise determination of its orbital inclination. The latter was estimated from polarimetry to 71 ± 9° (Rudy & Herman 1978). Plaskett’s star is usually considered a member of the Mon OB2 association. However, Linder et al. (2008) find a discrepancy between the luminosity and the dynamical masses of the two components of the system, which could be solved by assuming a greater distance for the star.

In many respects, HD 47129 offers a textbook example of evolutionary effects and the effects of interactions in massive binaries. Linder et al. (2008) used a disentangling technique, based on the method of González & Levato (2006), on high-resolution optical spectra to derive the spectral types O8 III/I and O7.5 V/III for the primary and secondary stars, respectively. The broad and shallow absorption lines of the secondary star suggest that this star rotates rapidly (vsini ~ 300 km s^-1), whereas the much sharper absorption lines of the primary indicate a projected rotational velocity of vsini ~ 75 km s^-1 (Linder et al. 2008). Using the CMFGEN model atmosphere code (Hillier & Miller 1998) to analyse the disentangled spectra, Linder et al. (2008) show that the atmosphere of the primary has an enhanced N and He abundance and a depletion of C. The secondary atmosphere is possibly depleted in N. These results indicate that the binary system is in a post-Roche lobe-overflow evolutionary stage where matter and angular momentum have been transferred from the primary to the secondary.

Because of its high rotation speed, the wind of the secondary star could be rotationally flattened. Such a situation would affect the properties of the wind interaction zone in this binary. This was confirmed by the studies of the Hα emission region by Wiggs & Gies (1992) and Linder et al. (2008), as well as by the study of the X-ray emission of the system by Linder et al. (2006). An alternative explanation for the origin of the equatorial wind of the secondary could be magnetic confinement. Indeed, Grunhut et al. (2013a) have recently reported the presence of a magnetic field in Plaskett’s star. In the framework of the magnetism in massive stars (MiMeS) survey, they detected Zeeman signatures from the rapidly rotating secondary in high-resolution Stokes V spectra. The strength of the highly organized field has to be at least 2850 ± 500 G, and variations compatible with rotational modulation of an oblique field were also found (Grunhut et al. 2013a).

Mahy et al. (2011) studied the CoRoT (Convection, Rotation, and planetary Transits, Baglin et al. 2006; Auvergne et al. 2009) light curve of Plaskett’s star and extracted 43 significant frequencies. Among these frequencies there are three major groups: 0.823 d^-1 and six harmonics, the orbital frequency 0.069 d^-1 and two harmonics, as well as two frequencies at 0.368 d^-1 and 0.399 d^-1. The latter two could be related to the rotation period of the secondary.

The discovery of the 0.823 d^-1 frequency by Mahy et al. (2011) prompted us to perform a spectroscopic monitoring of the system, to constrain the origin of this variation. In this paper, we report the results of this study. The paper is organized as follows. Section 2 describes the observations, data reduction and disentangling treatment. In Sect. 3, the line profile variations are discussed and the results of the Fourier analyses are given. Finally, in Sect. 4, we discuss the strengths and weaknesses of three possible scenarios and summarize our results and conclusions.

2. Observations and data reduction

2.1. OHP spectroscopy

Spectroscopic time series of HD 47129 were obtained during two six-nights observing campaigns in December 2009 and December 2010 at the Observatoire de Haute Provence (OHP, France). We used the Aurélie spectrograph fed by the 1.52 m telescope. The spectrograph was equipped with a 1200 lines mm^-1 grating blazed at 5000 Å and a CCD EEV42-20 detector with 2048 × 1024 pixels of 13.5 μm². Our set-up covered the wavelength domain from 4460 to 4670 Å with a resolving power of 20 000. Typical exposure times were 20−30 min, depending on the weather conditions, and a total of 118 spectra were collected. To achieve the most accurate wavelength calibration, Th–Ar lamp exposures were taken regularly over each observing night (typically once every 90 min).

The data were reduced using the MIDAS software provided by ESO. The normalization was done self-consistently using a series of continuum windows. Table 1 provides a summary of our observing campaigns and the characteristics of the sampling.

2.2. Disentangling

The goal of the present paper is to study and interpret the short-term line profile variability of Plaskett’s star. Because of the Doppler shifts and line blending induced by the orbital motion in such a binary system, it is not trivial to know which star triggers the variations. The first step in our analysis was thus to disentangle the spectra of the primary and secondary components. For this purpose, we applied our code based on the González & Levato (2006) technique to the data of the two observing campaigns simultaneously. This method allows to derive the radial velocities (RVs) as well as the individual spectra of the components. However, because of the very broad and shallow lines of the secondary, its RVs obtained with the disentangling were sometimes not well constrained. To avoid artefacts, we thus fixed the RVs of the secondary to the theoretical values derived from the orbital solution of Linder et al. (2008). The situation for the primary was better, as we recovered RVs that were in good agreement with the orbital solution of Linder et al. (2008): the dispersions of the RVs about the latter orbital solution were 11 and 6 km s^-1 for the He i λ 4471 and He ii λ 4542 lines, respectively. There were very little differences between the disentangled spectrum of the primary obtained with either the primary RVs being fixed to the orbital solution or allowing them to vary. In the subsequent analysis of the line profile variability, we thus adopted the mean spectra as reconstructed keeping the RVs of the primary fixed to the orbital solution of Linder et al. (2008). The reconstructed spectra of the individual components are displayed in Fig. 1. These individual spectra are important to correctly interpret the results of our line profile variability analysis. For instance, Fig. 1 reveals that the only secondary line that has at least some parts of its profile that are not blended with primary lines at some orbital phases is He i λ 4471. The latter line displays weak emission humps that were already noted by Linder et al. (2008).

Fig. 1 Mean disentangled spectra of the primary and secondary (shifted upwards by 0.1 continuum units). Important spectral lines are identified by the labels.

Fig. 2 Line-profile variations in the He i λ 4471 line in the 2009 (left) and 2010 (right) campaigns. The top panels show the dynamic spectra of the observed data in the heliocentric frame of reference. The bottom panels yield the secondary dynamic spectra. The labels on the vertical axis indicate the date of the first and last observation of each night in the format HJD − 2 455 000.

3. Variability analysis

3.1. He I λ 4471

The top panels of Fig. 2 display the dynamic spectra of the observed data in the heliocentric frame of reference (hereafter called the raw dynamic spectra), whilst the bottom panels illustrate the dynamic spectra after subtracting the mean disentangled primary spectrum and shifting the data into the secondary’s frame of reference (hereafter called the secondary dynamic spectra).

The raw dynamic spectra are dominated by the orbital motion of the relatively sharp primary line. These figures show which parts of the core of the secondary absorption are most affected by blends with the primary. We note that the observed position of the primary line slightly deviates from the position expected from the orbital solution. Another feature that can be seen on the first and third night of the 2009 campaign as well as on the third night of the 2010 data is a discrete depression that moves from the blue to the red wing of the broad secondary absorption.

Fig. 3 Fourier power spectra of the 2009 (left) and 2010 (right) time series of the secondary’s He i λ 4471 line. The top panels of each row show the full 2D periodograms (between 0 and 2 d-1), whilst the second and third panels illustrate the mean periodogram over the emission humps, respectively before and after prewhitening (see text). The spectral window is shown in the lower panel.

The secondary dynamic spectra reveal another feature that is more difficult to distinguish on the raw dynamic spectra. The strengths of the blue (4462−4465 Å) and red (4475−4480 Å) emission humps that flank the secondary’s absorption (see Fig. 1) undergo a strong modulation. The latter seems correlated with the presence of the discrete depression pointed out hereabove. Indeed, when the redwards moving discrete depression is observed, the emission humps are at their low level. The deviations between the observed and the theoretical positions and shapes of the primary profile lead to residuals when subtracting the mean primary spectrum that contaminate the secondary dynamic spectra between about 4466 and 4475 Å.

The dynamic spectra reveal that the variations extend over a wider wavelength range than the widths of the primary lines. We thus conclude that the secondary very probably displays variations, although, at this stage, we cannot exclude that the primary is also variable, especially in view of the residuals after subtracting the mean primary profile.

To further characterize the variations in the He i λ 4471 line, we performed a 2D Fourier analysis (Rauw et al. 2008) based on the periodogram method of Heck et al. (1985) refined by Gosset et al. (2001). This method is explicitly designed to account for the uneven sampling of the time series.

Applying our Fourier analysis to the time series of the observed spectra, without shifting them into the frame of reference of either star and/or subtracting the spectrum of the primary, results in power spectra dominated by low frequencies due to the orbital motion of the primary and, to a lesser extent, secondary lines. To get rid of these orbital frequencies and analyse the genuine short-term line profile variations, we have subtracted the mean disentangled spectrum of the primary shifted to the appropriate RV from the observed spectra. The resulting difference spectra were then further shifted into the secondary’s frame of reference. Therefore, we obtained the “individual spectra” of the secondary in the frame of reference of this star for each observation. The resulting 2D periodograms of the corresponding 2009 and 2010 time series are shown in Fig. 3.

To start, we focus on the emission humps (4462−4465 and 4475−4480 Å). We find that the mean periodogram of each observing campaign is dominated by a frequency near 1.65 d^-1 and its aliases. The error on the centre of a peak in the periodogram can be estimated as Δν = 0.1/ΔT¹. Here ΔT is the time interval between the first and the last observations of the time-series (see Table 1) and we assume that the uncertainty on the peak frequency amounts to 10% of the peak width. The uncertainties are of 0.02 d^-1 for the campaigns separately and of 3 × 10^-4 d^-1 when both campaigns are considered simultaneously². However, in the present case, these estimates are quite optimistic because they neither account for the signal to noise ratio of the observations nor for the impact of the treatment of the data prior to the Fourier analysis. Therefore, we estimate an overall error of ± 0.05 d^-1 on the peak values.

Because of this rather large error on the position of the peak, some caution is needed when comparing the results of our present analysis with those of the CoRoT light-curve analysis of Mahy et al. (2011). The latter authors reported a large number of frequencies and some of the peaks in our periodograms could be blends of two or more frequencies. Indeed, our spectroscopic time series are not sufficiently long to distinguish closely spaced frequencies. Moreover, our time series are heavily affected by 1 day aliasing, which was absent in the CoRoT data. For instance, the 1 day aliasing introduces some coupling between frequencies near 0.4 d^-1 and those near 1.6 d^-1. Nevertheless, we note that our 1.65 d^-1 peak nicely fits the 1.646 d^-1 frequency of the strongest signal found in the CoRoT data. Table 2 lists the frequencies from the CoRoT data that are potentially consistent with the frequencies found in the current analysis.

To check whether the photometric frequencies can account for the periodograms of our data, we used a prewhitening technique (Rauw et al. 2008). In this way, we found that the 2009 periodogram can be prewhitened with a single frequency (1.64 d^-1), whilst efficient prewhitening of the 2010 data requires three frequencies (1.65, 0.82, and 0.37 d^-1), which are the three strongest, non-orbital, frequencies reported by Mahy et al. (2011). The prewhitened periodograms are shown in Fig. 3. We thus conclude that the variations in the emission humps of the secondary’s He i λ 4471 line can be explained by the two dominant frequencies of the CoRoT photometry: 0.82 d^-1 and its first harmonic at 1.65 d^-1, with some possible contribution of the 0.37 d^-1 frequency.

The power spectra over the absorption core are affected by the residuals from the subtraction of the primary spectrum (see Fig. 3). For the 2010 data, the corresponding mean power spectrum features essentially the same set of frequencies as for the emission humps, although with larger residuals after prewhitening. For the 2009 campaign, the power spectrum is more complex, showing considerable power at low frequencies. This is especially the case at wavelengths between 4466 and 4471 Å because the primary line remains in this wavelength domain over the entire duration of the 2009 campaign (see Fig. 2). At longer wavelengths where the moving discrete depression is seen in the dynamic spectra, the dominant frequency is near 0.4 d^-1.

As pointed out above, the observed wavelengths of the primary line deviate from those expected from the orbital solution of Linder et al. (2008). We thus performed a Fourier analysis of the He i λ 4471 line in the primary’s frame of reference (adopting the Linder et al. 2008 orbital solution). We restrict our 2D Fourier analysis to the core of the primary’s line between 4470 and 4473 Å, which is roughly equivalent to the velocity range [− v sini,v sini]. For both observing campaigns, the power spectrum is completely explained by the 0.82 and 1.65 d^-1 frequencies and their aliases. Since these are the same frequencies as those found in the variations in the secondary’s emission humps, it seems quite possible that the variations are in fact due to the secondary star.

3.2. Other lines

The variability of the He ii λ 4542 and N iii λλ 4511−4518 lines is less prominent than in the case of He i λ 4471. The raw dynamic spectra of the He ii and N iii lines are dominated by the primary’s orbital motion and do not reveal clear redwards-moving discrete depressions (see Fig. 4). The secondary dynamic spectra are qualitatively similar to those of He i λ 4471, but the variability is more subtle.

The power spectra averaged over the He ii λ 4542 line in the frame of reference of the secondary star can be described by the combination of the 0.82 and 1.65 d^-1 frequencies and their aliases (see Fig. 5). For the 2010 data, a slightly cleaner prewhitening is achieved with 0.37 instead of 1.65 d^-1, and for the 2009 time series, a more efficient result is achieved including another harmonic frequency (2.47 d^-1). Shifting the spectra into the frame of reference of the primary star essentially yields the same results.

Fig. 4 Same as Fig. 2, but for the He ii λ 4542 line.

Fig. 5 Same as Fig. 3 but for the He ii λ 4542 secondary’s line. The prewhitening in the middle panels was done with the 0.82 and 1.65 d-1 frequencies.

Given the nitrogen overabundance of the primary star and the possible underabundance of the secondary reported by Linder et al. (2008), we expect a priori that the N iii lines should mainly be formed in the atmosphere of the primary star³. We have analysed the two strongest N iii lines (λλ 4511 and 4515) in the frame of reference of the primary. The resulting power spectra are again explained by the 0.82 and 1.65 d^-1 frequencies and (in the 2009 data) some contribution at 0.37 d^-1.

4. Discussion

The analysis of our spectroscopic time series of Plaskett’s star in the previous section revealed three main frequencies 0.82, 1.65, and 0.37 d^-1, although the properties of the latter two are interdependent because of the aliasing problem. Before we discuss the possible interpretation of these results, let us first recall the ingredients of Plaskett’s star as (we think) we know this system. The primary star is a rather slow rotator presenting chemical enrichment indicating that it is most probably in a post-RLOF status. The secondary is a fast rotator probably surrounded by a flattened wind. The orbit is circular and the stars currently do not fill their Roche lobes. A magnetic field, probably associated with the secondary, was detected. This magnetic field could be at the origin of the flattened wind structure.

With these ingredients in mind, three scenarios could possibly explain the observed line profile variations: pulsations, tidal interactions, and an oblique magnetic rotator.

4.1. Pulsations

Mahy et al. (2011) tentatively assigned the 0.82 d^-1 frequency and its six harmonics to a low-order (2 ≤ l ≤ 4) non-radial pulsation (NRP). Our Fourier analyses provide us with the amplitude (A) and the phase constant (ϕ₀) of the variation as a function of wavelength. These quantities can help us test the NRP scenario. Indeed, for NRPs with moderate degrees l the phase constant usually changes monotonically across the line profile and the difference between the phase constant in the red wing and in the blue wing is directly related to the values of l and m of the pulsation (e.g. Schrijvers & Telting 1999; Zima et al. 2004).

The amplitude and phase distributions of the 0.82 and 1.64 d^-1 frequencies of the He i λ 4471 secondary line are shown in Fig. 6. The uncertainties were evaluated via Monte Carlo simulations (see Rauw et al. 2008). Whilst these distributions are possibly affected by the residuals of the primary line, they are quite erratic and do not ressemble what one would expect for typical NRPs.

Fig. 6 Amplitude and phase variations for the He i λ 4471 line in the frame of reference of the secondary during the 2009 (left) and 2010 (right) campaigns. Two frequencies have been considered simultaneously: ν1 = 0.82 and ν2 = 1.65 d-1. The top panels yield the mean secondary profile, whilst the amplitudes associated with the two frequencies are shown in the second panels. The two lower panels indicate the phase constants of these frequencies.

Fig. 7 Same as Fig. 6, but in the frame of reference of the primary.

If we consider the amplitudes and phase constants in the frame of reference of the primary, the variations are less erratic and, at least in the case of the 0.82 d^-1 frequency, are more reminiscent of genuine NRPs. This would imply that the primary is responsible for the 0.82 and 1.65 d^-1 frequencies. In the frame of reference of the primary, significant amplitude of variability is found far away from the core of the primary line, out to −500 (−6.6 vsini) and + 400 km s^-1 (5.3 v sini) away from the line centre. Similar results are obtained for the He ii λ 4542 line. Variability over such an extended wavelength range around the line core is however rather unexpected for conventional NRPs in a slow rotator, therefore suggesting that the observed variability is most likely associated with the secondary star. The apparent variability of the primary could then be induced by travelling features of the much broader secondary line when they cross the position of the primary line.

4.2. Tidal interactions

The analysis of the line profile variations in the slightly eccentric binary system Spica (Palate et al. 2013a and references therein) has shown that tidal interactions could explain Spica’s line profile variations. Although Plaskett’s star is not an eccentric system, the non-synchronicity of the secondary rotation could lead to tidally induced variations such as those reported by Moreno et al. (2005). Therefore, we can wonder whether tidal interactions could be responsible for the variations reported hereabove.

We have thus tried to model Plaskett’s star with the TIDES code (see Moreno et al. 1999; 2005; 2011) combined to CoMBiSpeC⁴ (see Palate & Rauw 2012 and Palate et al. 2013b). The TIDES code computes the time-dependent shape of the stellar surface for eccentric and/or asynchronous systems accounting for centrifugal and Coriolis forces, gas pressure, viscous effects and gravitational interactions. From that, CoMBiSpeC computes the synthetic spectra of the stars at several orbital phases.

We stress that our goal here is not to fine tune the parameters that reproduce perfectly the spectra of the components of Plaskett’s star and their observed variations. Here, we rather wish to test whether or not tidal interactions can produce detectable variations in the spectra of this system. Indeed, fitting the spectra rigorously is a challenging task because of a large number of free/unknown parameters and the non-solar composition of the stars (see Linder et al. 2008); CoMBiSpeC currently only works with atmosphere models that have a solar composition. In the present paper, we have only adjusted the temperatures to reproduce the observed strength of the He i λ 4471 and He ii λ 4542 lines. The parameters used are listed in Table 3.

Mahy et al. (2011) suggested that the frequencies near 0.4 d^-1 could correspond to the rotational frequency of the secondary. The projected equatorial rotational velocity of the secondary has been evaluated to ~300 km s^-1 by Linder et al. 2008. These authors also derived a log g equal to 3.5 ± 0.1. Adopting an inclination of 67° as proposed by Mahy et al. (2011), and considering that the rotation frequency is 0.4 d^-1, the projected rotational velocity of 300 km s^-1 yields a radius of 16.1 R_⊙ for the secondary. The log g and minimum mass of the secondary (m_S sin³i = 47.3 M_⊙, Linder et al. 2008) along with the adopted inclination of 67° (Mahy et al. 2011) rather suggest a radius⁵ of R_⊙. In this latter case, the projected rotational velocity would be equal to ~385 km s^-1. The discrepancy between these rotational velocities can be explained by several factors: the uncertainties on the inclination (Mahy et al. 2011 indicated a range from 30 to 80°), the error on the rotational velocity determination, on the log g, and on the frequency.

Our simulations indicate that tidal interactions can produce visible variations in the lines of the secondary. Because the primary is thought to be in synchronous rotation in a circular orbit, it is expected that there is no variability of the primary due to tidal interactions. The strength of the simulated variations is comparable to the observed ones. Reducing the secondary radius to 18 R_⊙ slightly lowers the amplitude of the tidally induced variations, although they remain well visible in the synthetic spectra. Figure 8 displays the simulated line profile variations in the secondary during the orbital cycle. Whilst our code cannot simulate the impact of the tidal interactions on a co-rotating confined wind, it seems very plausible that the tidal interactions also affect the latter and thus extend into the emission humps. This result brings up the interesting possibility that breaking tidal waves near the secondary’s equator could produce an enhanced mass-loss in this region (see Osaki 1999). Combined with wind confinement due to the magnetic field (Grunhut et al. 2013a; 2013b), this effect could produce the equatorial torrus of circumstellar material that was advocated by Wiggs & Gies (1992) and Linder et al. (2006; 2008). The much smaller radius (10.7 R_⊙) proposed by Grunhut et al. (2013a) not only fails to reproduce the log g value inferred from the model atmosphere code fits (see below), but would also considerably reduce the amplitude of surface oscillations due to tidal interactions.

We have performed a Fourier analysis of the individual simulated spectra of the secondary stars. The power spectrum displays a peak at the secondary’s rotation frequency imposed in the model, but does not show peaks at the harmonic frequencies. Tidal interactions can thus induce profile variations, provided the stellar radii are sufficiently large. However, the current model does not explain the full frequency content of HD 47129’s spectroscopic and photometric variations.

Fig. 8 Line profile variations in the synthetic spectra of the secondary during the orbital cycle. The spectra (separated by 0.05 in phase) have been shifted vertically by 0.1 continuum unit for clarity.

4.3. Magnetically confined winds

Grunhut et al. (2013a) reported on the discovery of a magnetic field in Plaskett’s star with a minimum surface dipolar strength of 2.8 kG at the magnetic poles. Such a strong magnetic field could confine the secondary’s stellar wind, thereby producing the emission humps seen around the He i λ 4471 line as well as the double-peaked Hα emission. Using additional spectropolarimetric observations, Grunhut et al. (2013b) inferred a tilt angle of the magnetic axis of (80 ± 5)° with respect to the rotation axis. This would result in a magnetic configuration very similar to HD 57682 (Grunhut et al. 2012). One would then expect a double-wave modulation of the strength of the emission features. Based on a periodicity search in the variations in the Hα equivalent width and the longitudinal magnetic field, Grunhut et al. (2013b) inferred a rotational period of 1.215 d for the secondary, corresponding to the 0.82 d^-1 frequency.

This scenario could thus explain our detection of the 0.82 and 1.65 d^-1 frequencies quite naturally in the emission humps. However, adopting 0.82 d^-1 as the secondary’s rotational frequency would lead to a radius estimate from the projected rotational velocity that is twice smaller than what we have obtained above, thus enhancing the discrepancy between the radius estimated from v sini and the one estimated from log g.

Concerning the latter discrepancy, we note that Wiggs & Gies (1992) previously reported a 2.78 d period in the equivalent width of the Hα emission wings⁶.

4.4. Summary and conclusion

Our analysis of the He i λ 4471, N iii λλ 4510-4518, and He iiλ 4542 lines of Plaskett’s star revealed variations at several frequencies: 0.37, 0.82, and 1.65d^-1. The strongest variations are found for the He i line.

Considering the current knowledge of the system, it is difficult to assign the spectroscopic variability to either the primary or the secondary star. Indeed, our analysis yielded conflicting results in this respect: whilst the velocity range over which variations are detected clearly favours the secondary star, the pattern of the amplitude and phase of the variations makes more sense if they originate in the primary star. As a result, none of the current scenarios (pulsations, tidal interactions, magnetically confined winds) accounts for all the aspects of the observed variations.

Acknowledgments

We would like to thank Pr. G. Koenigsberger who provided us with the TIDES code and helped us to use it. We are grateful to the referee, Dr. Jason Grunhut who helped us improve this paper. We acknowledge support through the XMM/INTEGRAL PRODEX contract (Belspo), from the Fonds de Recherche Scientifique (FRS/FNRS), as well as by the Communauté Française de Belgique − Action de recherche concertée − Académie Wallonie − Europe.

References

All Tables

All Figures

	Fig. 1 Mean disentangled spectra of the primary and secondary (shifted upwards by 0.1 continuum units). Important spectral lines are identified by the labels.
In the text

	Fig. 2 Line-profile variations in the He i λ 4471 line in the 2009 (left) and 2010 (right) campaigns. The top panels show the dynamic spectra of the observed data in the heliocentric frame of reference. The bottom panels yield the secondary dynamic spectra. The labels on the vertical axis indicate the date of the first and last observation of each night in the format HJD − 2 455 000.
In the text

	Fig. 3 Fourier power spectra of the 2009 (left) and 2010 (right) time series of the secondary’s He i λ 4471 line. The top panels of each row show the full 2D periodograms (between 0 and 2 d-1), whilst the second and third panels illustrate the mean periodogram over the emission humps, respectively before and after prewhitening (see text). The spectral window is shown in the lower panel.
In the text

	Fig. 4 Same as Fig. 2, but for the He ii λ 4542 line.
In the text

	Fig. 5 Same as Fig. 3 but for the He ii λ 4542 secondary’s line. The prewhitening in the middle panels was done with the 0.82 and 1.65 d-1 frequencies.
In the text

Fig. 6 Amplitude and phase variations for the He i λ 4471 line in the frame of reference of the secondary during the 2009 (left) and 2010 (right) campaigns. Two frequencies have been considered simultaneously: ν1 = 0.82 and ν2 = 1.65 d-1. The top panels yield the mean secondary profile, whilst the amplitudes associated with the two frequencies are shown in the second panels. The two lower panels indicate the phase constants of these frequencies.

In the text

	Fig. 7 Same as Fig. 6, but in the frame of reference of the primary.
In the text

	Fig. 8 Line profile variations in the synthetic spectra of the secondary during the orbital cycle. The spectra (separated by 0.05 in phase) have been shifted vertically by 0.1 continuum unit for clarity.
In the text