Arxiv:1004.4893V1 [Astro-Ph.SR] 27 Apr 2010 Oni Fcanada
Total Page:16
File Type:pdf, Size:1020Kb
Large-Scale Periodic Variability of the Wind of the Wolf-Rayet Star WR 1 (HD 4004) A.-N. Chen´e∗, ∗∗ Canadian Gemini Office, HIA/NRC of Canada, 5071, West Saanich Road, Victoria (BC), V9E 2E7, Canada [email protected] and N. St-Louis D´epartement de Physique and CRAQ, Universit´ede Montr´eal, C.P. 6128, Succ. Centre-Ville, Montr´eal, Qu´ebec, H3C 3J7, Canada [email protected] ABSTRACT We present the results of an intensive photometric and spectroscopic monitoring campaign of the WN4 Wolf-Rayet (WR) star WR 1=HD 4004. Our broadband V photometry covering a +0.6 timespan of 91 days shows variability with a period of P=16.9−0.3 days. The same period is also found in our spectral data. The light-curve is non-sinusoidal with hints of a gradual change in its shape as a function of time. The photometric variations nevertheless remain coherent over several cycles and we esti- mate that the coherence timescale of the light-curve is of the order of 60 days. The spectroscopy shows large-scale line-profile variability which can be interpreted as excess emission peaks moving from one side of the profile to the other on a timescale of several days. Although we cannot unequivocally exclude the unlikely possibility that WR 1 is a binary, we propose that the nature of the variability we have found strongly suggests that it is due to the presence in the wind of the WR star of large-scale structures, most likely Co-rotating Interaction Regions (CIRs), which are predicted to arise in inherently unstable radiatively driven winds when they are perturbed at their base. We also suggest that variability observed in WR6, WR134 and WR 137 is of the same nature. Finally, assuming that the period of CIRs is related to the arXiv:1004.4893v1 [astro-ph.SR] 27 Apr 2010 rotational period, we estimate the rotation rate of the four stars for which sufficient monitoring −1 has been carried out; i.e. vrot=6.5, 40, 70 and 275 km s for WR1, WR6, WR134 and WR137, respectively. Subject headings: stars: Wolf-Rayet — stars: individual(WR 1, HD 4004) — stars: winds,outflows — stars: rotation *Guest investigator, Dominion Astrophysical Observa- de Recherche Scientifique of France, and the University of tory, Herzberg Institute of Astrophysics, National Research Hawaii. Also based on observations obtained at the Obser- Council of Canada. vatoire du Mont M´egantic with is operated by the Centre **Based on observations obtained at the Canada-France- de Recherche en Astrophysique du Qu´ebec and the Obser- Hawaii Telescope (CFHT) which is operated by the Na- vatoire de Haute-Provence which is operated by the Institut tional Research Council of Canada, the Institut Na- National des Sciences de l’Univers of the Centre National tional des Sciences de l’Univers of the Centre National de la Recherche Scientifique of France. 1 1. Introduction pletely excluded. If the large-scale lpv and the photometric vari- The winds of Wolf-Rayet (WR) stars are ability of single WR stars can be associated with well-known to be non-uniform on small phys- CIRs, it may imply that the period of these vari- ical scales. Because of the inherently unsta- abilities corresponds directly to the rotational ble nature of radiatively-driven hot stellar winds period of the star. Indeed, Cranmer & Owocki (Owocki, Castor & Rybicki 1988), instabilities, (1996) propose that “spots” fixed to the stellar which reveal themselves in the spectrum as narrow surface, caused either by pulsations or magnetic excess emission peaks superposed on the broad field activity, are at the origin of the CIRs. Hence, emission lines, appear stochastically, propagate in once the radius at which the CIR originates is the wind and disappear after several hours (e.g. known, we can determine the rotational velocity Moffat et al. 1988). This clumpiness leads to of the star at that point. It follows that the ro- emission lines that are variable at a level of up to tation rate of at least some WR stars could in ∼ 5 % of the line flux (St-Louis et al. 2009). These principle be determined by carrying out a system- changes are, of course, random and no periodicity atic investigation of the variability of all single is expected. WR stars. As a first step, St-Louis et al. (2009) In certain cases, WR stars also display large- and Chen´e& St-Louis (in prep) set out to identify scale line-profile variability (lpv). For some stars new candidates for CIR-type variability by con- this spectroscopic variability has been shown to be ducting a survey of all apparently single Galactic periodic and to originate from a different physical WR stars brighter than v ∼12.5th magnitude. For processes. Of course, known massive WR+O bi- each star in their sample, they obtained 4–5 spec- naries produce clear periodic radial-velocity (RV) tra which allowed them to establish a list of WR variations from the orbital motion of the stars. It stars showing large-scale lpv. The next step is to has also been realized in the past two decades that observe intensively each of the candidates in order in such systems the winds from the stars collide, to verify and determine the periodicity. forming a shock cone that wraps around the star In the above-mentioned survey, WR 1 was one with the smaller momentum flux (generally the of the most striking cases of large-scale lpv with O star). This interaction induces kinematically- changes reaching 8–10% of the line flux and easily characteristic variations in the line profiles that distinguishable large scale subpeaks superposed on are also periodic. the broad wind emission profiles. Consequently, it However, large-scale lpv has also been found was the first WR star on which we concentrated in WR stars that are not known to be binaries. our efforts. Several previous studies claiming un- The two most studied cases, WR 6 (WN4) and reproducible periods for the variability of this star WR 134 (WN6), have been monitored intensively have been summarized in Morel et al. (1999b). in photometry and spectroscopy by Morel et al. Later, Niedzielski (2000a,b) demonstrated that (1997, 1998, 1999a). In both cases, these authors the spectrum of WR1 varied greatly from night- have observed periodic large-scale lpv and a com- to-night reaching 50% in the equivalent width of plex light-curve with the same periodicity. To the Heiλ5876 line. On the other hand, they found that list of two stars, we may also add WR 137, that line-profile changes during one night are much a WC7pd star in a long-period binary system. smaller with a typical scatter of 0.5 A˚ in equivalent This star shows large-scale lpv with a period of width for the Heiiλ5411, Civλ5808 and Heiλ5876 1.2 d that is unlikely to be related to its O9 lines. The changes, however, were found to be sys- companion nor to the wind-wind collision zone tematic during the course of an entire night and (Lef`evre et al. 2005a). This behaviour is most reach a total of 3-4 A.˚ Moreover, the search for pe- likely explained by the presence of large-scale den- riods smaller than 2 days failed and only indica- sity structures in the wind, such as co-rotating tions of long-term variability could be suggested. interaction regions (CIRs) (Cranmer & Owocki The investigation of photometric variability over 1996; Fullerton et al. 1997), but the possibility more than 16 days carried out by Morel et al. of the presence of a compact companion or of (1999b) did not lead to the identification of a pe- a low-mass main-sequence star cannot be com- riod. Consequently, the period search must be 2 done using data taken over a time range of at least FWHM of the point spread function (PSF) and an twice as long. More recently, Flores et al. (2007) annulus of sky was selected with inner and outer claimed a period of P=7.684 days in lpv of this radii of respectively 4 and 8 times the FWHM of star which they attribute to the ejecta of streams the PSF. or jets from the stellar surface. Once the photometry was obtained for all In this paper, we present the results of an in- frames, we then applied a correction for the differ- tensive monitoring campaign of WR 1 extending ential refraction of the atmosphere. To do so, we over several weeks in photometry and spectroscopy needed to know how the flux of a star varies with with the aim of determining the nature of the airmass and according to the slope of its spec- variability and, eventually, if found to be associ- trum. For each star in our field of view, a relative ated with CIRs as suspected from our survey ob- (B − V ) color was obtained by subtracting the av- servations, the rotation period of WR1. In Sec- erage fluxes measured in the B and V -band frames tion 2, we present our photometric and spectro- observed at the same airmass. We then obtained scopic observations and the data reduction proce- the slope M of the linear relation describing the dures. In Sections 3 and 4 we describe our results, variability of the flux as a function of airmass for and in Section 5 we discuss the possible interpre- each star during the best clear nights. Although a tations. Finally, our conclusions are presented in certain number of stars can present intrinsic flux Section 5.3.