The structured winds of single massive stars can be classified into two broad groups: stochastic structure and organized structure. While the former is typically identified with clumping, the latter is typically associated with rotational modulations, particularly the paradigm of Co-rotating Interaction Regions (CIRs). While CIRs have been explored extensively in the UV band, and moderately in the X-ray and optical, here we evaluate radio variability from CIR structures assuming free-free opacity in a dense wind. Our goal is to conduct a broad parameter study to assess the observational feasibility, and to this end, we adopt a phenomenological model for a CIR that threads an otherwise spherical wind. We find that under reasonable assumptions, it is possible to obtain radio variability at the 10% level. The detailed structure of the folded light curve depends not only on the curvature of the CIR, the density contrast of the CIR relative to the wind, and viewing inclination, but also on wavelength. Comparing light curves at different wavelengths, we find that the amplitude can change, that there can be phase shifts in the waveform, and the the entire waveform itself can change. These characterstics could be exploited to detect the presence of CIRs in dense, hot winds.

Wolf-Rayet stars have strong, hot winds, with mass-loss rates at least a factor of ten greater than their O-star progenitors, although their terminal wind speeds are similar. In this paper we use the technique of multiband linear polarimetry to extract information on the global asymmetry of the wind in a sample of 47 bright Galactic WR stars. Our observations also include time-dependent observations of 17 stars in the sample. The path to our goal includes removing the dominating component of wavelength-dependent interstellar polarization (ISP), which normally follows the well-known Serkowski law. We include a wavelength-dependent ISP position angle parameter in our ISP law and find that 15 stars show significant results for this parameter. We detect a significant component of wavelength-independent polarization due to electron scattering in the wind for 10 cases, with most WR stars showing none at the $\sim$0.05% level precision of our data. The intrinsically polarized stars can be explained with binary interaction, large-scale wind structure, and clumping. We also found that 5 stars out of 19 observed with the Str\"omgren $b$ filter (probing the complex $\lambda$4600--4700 emission line region) have significant residuals from the ISP law and propose that this is due to binary illumination or wind clumping. We provide a useful catalogue of ISP for 47 bright Galactic WR stars and upper limits on the possible level of intrinsic polarization.

PolStar is an Explorer-class far ultraviolet (FUV) spectropolarimetry mission designed to target massive stars and their environments. PolStar will take advantage of resonance lines only available in the FUV to measure for the first time the magnetic and wind environment around massive stars to constrain models of rotation and mass loss.

Among Wolf-Rayet stars, those of subtype WN8 are the intrinsically most variable. We have explored the long-term photometric variability of the brightest known WN8 star, WR 40, through four contiguous months of time-resolved, single-passband optical photometry with the BRIght Target Explorer (BRITE) nanosatellite mission. The Fourier transform of the observed light-curve reveals that the strong light variability exhibited by WR 40 is dominated by many randomly-triggered, transient, low-frequency signals. We establish a model in which the whole wind consists of stochastic clumps following an outflow visibility promptly rising to peak brightness upon clump emergence from the optically thick pseudo-photosphere in the wind, followed by a gradual decay according to the right-half of a Gaussian. Free electrons in each clump scatter continuum light from the star. We explore a scenario where the clump size follows a power-law distribution, and another one with an ensemble of clumps of constant size. Both scenarios yield simulated light curves morphologically resembling the observed light curve remarkably well, indicating that one cannot uniquely constrain the details of clump size distribution with only a photometric light curve. Nevertheless, independent evidence favours a negative-index power law, as seen in many other astrophysical turbulent media.

We use a 3D Monte Carlo radiative transfer code to study the polarimetric and photometric variability from stationary corotating interaction regions (CIR) in the wind of massive stars. Our CIRs are approximated by Archimedean spirals of higher (or lower) density formed in a spherical wind originating from the star and we also made allowance for a bright gaussian spot at the base of the CIR. Comparing results from our code to previous analytical calculations in the optically thin case, we find differences which we attribute mainly to a better estimation of the total unpolarized flux reaching the observer. In the optically thick case, the differences with the analytical calculations are much larger, as multiple scattering introduces extra complexities including occultation effects. The addition of a gaussian spot does not alter the shape of the polarization curve significantly but does create a small excess in polarization. On the other hand, the effect can be larger on the light curve and can become dominant over the resulting CIR, depending on the spot parameters and density of the wind.

Massive Wolf-Rayet (WR) stars dominate the radiative and mechanical energy budget of galaxies and probe a critical phase in the evolution of massive stars prior to core-collapse. It is not known whether core He-burning WR stars (classical WR, cWR) form predominantly through wind-stripping (w-WR) or binary stripping (b-WR). With spectroscopy of WR binaries so-far largely avoided due to its complexity, our study focuses on the 44 WR binaries / binary candidates of the Large Magellanic Cloud (LMC, metallicity Z~0.5 Zsun), identified on the basis of radial velocity variations, composite spectra, or high X-ray luminosities. Relying on a diverse spectroscopic database, we aim to derive the physical and orbital parameters of our targets, confronting evolution models of evolved massive stars at sub-solar metallicity, and constraining the impact of binary interaction in forming them. Spectroscopy is performed using the Potsdam Wolf-Rayet (PoWR) code and cross-correlation techniques. Disentanglement is performed using the code Spectangular or the shift-and-add algorithm. Evolutionary status is interpreted using the Binary Population and Spectral Synthesis (BPASS) code, exploring binary interaction and chemically-homogeneous evolution. No obvious dichotomy in the locations of apparently-single and binary WN stars on the Hertzsprung-Russell diagram is apparent. According to commonly used stellar evolution models (BPASS, Geneva), most apparently-single WN stars could not have formed as single stars, implying that they were stripped by an undetected companion. Otherwise, it must follow that pre-WR mass-loss/mixing (e.g., during the red supergiant phase) are strongly underestimated in standard stellar evolution models.

The hot WN star WR2 (HD6327) has been claimed to have many singular characteristics. To explain its unusually rounded and relatively weak emission line profiles, it has been proposed that WR2 is rotating close to break-up with a magnetically confined wind. Alternatively, the line profiles could be explained by the dilution of WR2's spectrum by that of a companion. In this paper, we present a study of WR2 using near-infrared AO imaging and optical spectroscopy and polarimetry. Our spectra reveal the presence of weak photospheric absorption lines from a ~B2.5-4V companion, which however contributes only ~5-10% to the total light, suggesting that the companion is a background object. Therefore, its flux cannot be causing any significant dilution of the WR star's emission lines. The absence of intrinsic linear continuum polarization from WR2 does not support the proposed fast rotation. Our Stokes V spectrum was not of sufficient quality to test the presence of a moderately strong organized magnetic field but our new modelling indicates that to confine the wind the putative magnetic field must be significantly stronger than was previously suggested sufficiently strong as to make its presence implausible.

The most violent strange mode pulsations (SMPs) are expected in classical WR stars, the bare, compact helium-burning cores of evolved massive stars, where SMPs are predicted to manifest themselves in cyclic photometric variability with periods ranging from minutes to hours and amplitudes in the several to tens of mmag range. Since WR\,2 is one of the hottest WR stars known, it should make it an excellent candidate for SMPs. However, our high-precision near infrared photometric campaign shows no periodic variations with an amplitude higher than 2.5 mmag on a timescale of ~6 hours, which is more than an order of magnitude smaller than the theoretical predictions.

The intriguing WN4b star WR6 has been known to display epoch-dependent spectroscopic, photometric and polarimetric variability for several decades. In this paper, we set out to verify if a simplified analytical model in which Corotating Interaction Regions (CIRs) threading an otherwise spherical wind is able to reproduce the many broadband continuum light curves from the literature with a reasonable set of parameters. We modified the optically thin model we developed in Ignace, St-Louis & Proulx-Giraldeau (2015) to approximately account for multiple scattering and used it to fit 13 separate datasets of this star. By including two CIRs in the wind, we obtained reasonable fits for all datasets with coherent values for the inclination of the rotation axis ($i_0=166^{\circ}$) and for its orientation in the plane of the sky, although in the latter case we obtained two equally acceptable values ($\psi=63^{\circ}$ and $\psi=152^{\circ}$) from the polarimetry. Additional line profile variation simulations using the Sobolev approximation for the line transfer allowed us to eliminate the $\psi=152^{\circ}$ solution. With the adopted configuration ($i_0=166^{\circ}$ and $\psi=63^{\circ}$), we were able to reproduce all datasets relatively well with two CIRs located near the stellar equator and always separated by $\sim 90^{\circ}$ in longitude. The epoch-dependency comes from the fact that these CIRs migrate along the surface of the star. Density contrasts smaller than a factor of two and large opening angles for the CIR ($\beta \geq 35^{\circ}$) were found to best reproduce the type of spectroscopic variability reported in the literature.

From $5.5$ months of dual-band optical photometric monitoring at the $1$ mmag level, BRITE-Constellation has revealed two simultaneous types of variability in the O4I(n)fp star $\zeta$ Puppis: one single periodic non-sinusoidal component superimposed on a stochastic component. The monoperiodic component is the $1.78$ d signal previously detected by Coriolis/SMEI, but this time along with a prominent first harmonic. The shape of this signal changes over time, a behaviour that is incompatible with stellar oscillations but consistent with rotational modulation arising from evolving bright surface inhomogeneities. By means of a constrained non-linear light curve inversion algorithm we mapped the locations of the bright surface spots and traced their evolution. Our simultaneous ground-based multi-site spectroscopic monitoring of the star unveiled cyclical modulation of its He II $\lambda4686$ wind emission line with the $1.78$-day rotation period, showing signatures of Corotating Interaction Regions (CIRs) that turn out to be driven by the bright photospheric spots observed by BRITE. Traces of wind clumps are also observed in the He II $\lambda4686$ line and are correlated with the amplitudes of the stochastic component of the light variations probed by BRITE at the photosphere, suggesting that the BRITE observations additionally unveiled the photospheric drivers of wind clumps in $\zeta$ Pup and that the clumping phenomenon starts at the very base of the wind. The origins of both the bright surface inhomogeneities and the stochastic light variations remain unknown, but a subsurface convective zone might play an important role in the generation of these two types of photospheric variability.

We report on the first multi-color precision light curve of the bright Wolf-Rayet binary $\gamma^2$ Velorum, obtained over six months with the nanosatellites in the BRITE- Constellation fleet. In parallel, we obtained 488 high-resolution optical spectra of the system. In this first report on the datasets, we revise the spectroscopic orbit and report on the bulk properties of the colliding winds. We find a dependence of both the light curve and excess emission properties that scales with the inverse of the binary separation. When analyzing the spectroscopic properties in combination with the photometry, we find that the phase dependence is caused only by excess emission in the lines, and not from a changing continuum. We also detect a narrow, high-velocity absorption component from the He I $\lambda$5876 transition, which appears twice in the orbit. We calculate smoothed-particle hydrodynamical simulations of the colliding winds and can accurately associate the absorption from He I to the leading and trailing arms of the wind shock cone passing tangentially through our line of sight. The simulations also explain the general strength and kinematics of the emission excess observed in wind lines such as C III $\lambda$5696 of the system. These results represent the first in a series of investigations into the winds and properties of $\gamma^2$ Velorum through multi-technique and multi-wavelength observational campaigns.

WR 148 (HD 197406) is an extreme runaway system considered to be a potential candidate for a short-period (4.3173 d) rare WR + compact object binary. Provided with new high resolution, high signal-to-noise spectra from the Keck observatory, we determine the orbital parameters for both the primary WR and the secondary, yielding respective projected orbital velocity amplitudes of $88.1\pm3.8$ km s$^{-1}$ and $79.2\pm3.1$ km s$^{-1}$ and implying a mass ratio of $1.1\pm0.1$. We then apply the shift-and-add technique to disentangle the spectra and obtain spectra compatible with a WN7ha and an O4-6 star. Considering an orbital inclination of $\sim67^\circ$, derived from previous polarimetry observations, the system's total mass would be a mere 2-3 M$_{\odot}$ , an unprecedented result for a putative massive binary system. However, a system comprising a 37 M$_{\odot}$ secondary (typical mass of an O5V star) and a 33 M$_{\odot}$ primary (given the mass ratio) would infer an inclination of $\sim18^\circ$. We therefore reconsider the previous methods of deriving the orbital inclination based on time-dependent polarimetry and photometry. While the polarimetric results are inconclusive requiring better data, the photometric results favour low inclinations. Finally, we compute WR 148's space velocity and retrace the runaway's trajectory back to the Galactic plane (GP). With an ejection velocity of $198\pm27$ km s$^{-1}$ and a travel time of $4.7\pm0.8$ Myr to reach its current location, WR 148 was most likely ejected via dynamical interactions in a young cluster.

We report on interferometric observations with the CHARA Array of two classical Wolf-Rayet stars in suspected binary systems, namely WR 137 and WR 138. In both cases, we resolve the component stars to be separated by a few milliarcseconds. The data were collected in the H-band, and provide a measure of the fractional flux for both stars in each system. We find that the WR star is the dominant H-band light source in both systems ($f_{\rm WR, 137} = 0.59\pm0.04$; $f_{\rm WR, 138} = 0.67\pm0.01$), which is confirmed through both comparisons with estimated fundamental parameters for WR stars and O dwarfs, as well as through spectral modeling of each system. Our spectral modeling also provides fundamental parameters for the stars and winds in these systems. The results on WR 138 provide evidence that it is a binary system which may have gone through a previous mass-transfer episode to create the WR star. The separation and position of the stars in the WR 137 system together with previous results from the IOTA interferometer provides evidence that the binary is seen nearly edge-on. The possible edge-on orbit of WR 137 aligns well with the dust production site imaged by the Hubble Space Telescope during a previous periastron passage, showing that the dust production may be concentrated in the orbital plane.

During the summer of 2013, a 4-month spectroscopic campaign took place to observe the variabilities in three Wolf-Rayet stars. The spectroscopic data have been analyzed for WR 134 (WN6b), to better understand its behaviour and long-term periodicity, which we interpret as arising from corotating interaction regions (CIRs) in the wind. By analyzing the variability of the He II $\lambda$5411 emission line, the previously identified period was refined to P = 2.255 $\pm$ 0.008 (s.d.) days. The coherency time of the variability, which we associate with the lifetime of the CIRs in the wind, was deduced to be 40 $\pm$ 6 days, or $\sim$ 18 cycles, by cross-correlating the variability patterns as a function of time. When comparing the phased observational grayscale difference images with theoretical grayscales previously calculated from models including CIRs in an optically thin stellar wind, we find that two CIRs were likely present. A separation in longitude of $\Delta \phi \simeq$ 90$^{\circ}$ was determined between the two CIRs and we suggest that the different maximum velocities that they reach indicate that they emerge from different latitudes. We have also been able to detect observational signatures of the CIRs in other spectral lines (C IV $\lambda\lambda$5802,5812 and He I $\lambda$5876). Furthermore, a DAC was found to be present simultaneously with the CIR signatures detected in the He I $\lambda$5876 emission line which is consistent with the proposed geometry of the large-scale structures in the wind. Small-scale structures also show a presence in the wind, simultaneously with the larger scale structures, showing that they do in fact co-exist.

MWC 314 is a bright candidate luminous blue variable that resides in a fairly close binary system, with an orbital period of 60.753$\pm$0.003 d. We observed MWC 314 with a combination of optical spectroscopy, broad-band ground- and space-based photometry, as well as with long baseline, near-infrared interferometry. We have revised the single-lined spectroscopic orbit and explored the photometric variability. The orbital light curve displays two minima each orbit that can be partially explained in terms of the tidal distortion of the primary that occurs around the time of periastron. The emission lines in the system are often double-peaked and stationary in their kinematics, indicative of a circumbinary disc. We find that the stellar wind or circumbinary disc is partially resolved in the K\prime-band with the longest baselines of the CHARA Array. From this analysis, we provide a simple, qualitative model in an attempt to explain the observations. From the assumption of Roche Lobe overflow and tidal synchronisation at periastron, we estimate the component masses to be M1 $\approx 5$ M$_\odot$ and M2$\approx 15$ M$_\odot$, which indicates a mass of the LBV that is extremely low. In addition to the orbital modulation, we discovered two pulsational modes with the MOST satellite. These modes are easily supported by a low-mass hydrogen-poor star, but cannot be easily supported by a star with the parameters of an LBV. The combination of these results provides evidence that the primary star was likely never a normal LBV, but rather is the product of binary interactions. As such, this system presents opportunities for studying mass-transfer and binary evolution with many observational techniques.

Massive star winds are complex radiation-hydrodynamic (sometimes magnetohydrodynamic) outflows that are propelled by their enormously strong luminosities. The winds are often found to be structured and variable, but can also display periodic or quasi-periodic behavior in a variety of wind diagnostics. The regular variations observed in putatively single stars, especially in UV wind lines, have often been attributed to corotating interaction regions (CIRs) like those seen in the solar wind. We present light curves for variable polarization from winds with CIR structures. We develop a model for a time-independent CIR based on a kinematical description. Assuming optically thin electron scattering, we explore the range of polarimetric light curves that result as the curvature, latitude, and number of CIRs are varied. We find that a diverse array of variable polarizations result from an exploration of cases. The net polarization from an unresolved source is weighted more toward the inner radii of the wind. Given that most massive stars have relatively fast winds compared to their rotation speeds, CIRs tend to be conical at inner radii, transitioning to a spiral shape at a few to several stellar radii in the wind. Winds with a single CIR structure lead to easily identifiable polarization signatures. By contrast allowing for multiple CIRs, all emerging from a range of azimuth and latitude positions at the star, can yield complex polarimetric behavior. Although our model is based on some simplifying assumptions, it produces qualitative behavior that we expect to be robust, and this has allowed us to explore a wide range of CIR configurations that will prove useful for interpreting polarimetric data.

In this paper, high quality spectropolarimetric observations of the Wolf-Rayet (WR) Star WR1 (HD4004) obtained with ESPaDOnS at CFHT are presented. All major emission lines present in the spectrum show depolarization in the relative Stokes parameters Q/I and U/I. From the behaviour of the amount of line depolarization as a function of line strength, the intrinsic continuum light polarization of WR1 is estimated to be P/I=0.443 \pm 0.028 % with an angle of \theta=-26.2^o. Although such a level of polarization could in principle be caused by a wind flattened by fast rotation, the scenario in which it is a consequence of the presence in the wind of Corotating Interaction Regions (CIRs) is preferred. This is supported by previous photometric and spectroscopic observations showing periodic variations with a period of 16.9 days. This is now the third WR star thought to exhibit CIRs in its wind that is found to have line depolarization. Previous authors have found a strong correlation between line depolarization and the presence of an ejected nebula, which they interpret as a sign that the star has reached the WR phase relatively recently since the nebula are thought to dissipate very fast. In the cases where the presence of CIRs in the wind is favoured to explain the depolarization across spectral lines, the above mentioned correlation may be indicating that those massive stars have only very recently transited from the previous evolutionary phase to the WR phase.

On the final day of the Stellar Polarimetry conference, participants split up into three "breakout sessions" to discuss the future of the field in the areas of instrumentation, upcoming opportunities, and community priorities. This contribution compiles the major recommendations arising from each breakout session. We hope that the polarimetric community will find these ideas useful as we consider how to maintain the vitality of polarimetry in the coming years.

This study is the second part of a survey searching for large-scale spectroscopic variability in apparently single Wolf-Rayet (WR) stars. In a previous paper (Paper I), we described and characterized the spectroscopic variability level of 25 WR stars observable from the northern hemisphere and found 3 new candidates presenting large-scale wind variability, potentially originating from large-scale structures named Co-rotating Interaction Regions (CIRs). In this second paper, we discuss an additional 39 stars observable from the southern hemisphere. For each star in our sample, we obtained 4-5 high-resolution spectra with a signal-to-noise ratio of ~100 and determined its variability level using the approach described in Paper I. In total, 10 new stars are found to show large-scale spectral variability of which 7 present CIR-type changes (WR 8, WR 44, WR 55, WR 58, WR 61, WR 63, WR 100). Of the remaining stars, 20 were found to show small-amplitude changes and 9 were found to show no spectral variability as far as can be concluded from the data in hand. Also, we discuss the spectroscopic variability level of all single galactic WR stars that are brighter than v~12.5, and some WR stars with 12.5 < v <= 13.5; i.e. all the stars presented in our two papers and 4 more stars for which spectra have already been published in the literature. We find that 23/68 stars (33.8 %) present large-scale variability, but only 12/54 stars (~22.1 %) are potentially of CIR-type. Also, we find 31/68 stars (45.6 %) that only show small-scale variability, most likely due to clumping in the wind. Finally, no spectral variability is detected based on the data in hand for 14/68 (20.6 %) stars. Interestingly, the variability with the highest amplitude also have the widest mean velocity dispersion.

A 30-day contiguous photometric run with the MOST satellite on the WN5-6b star WR 110 (HD 165688) reveals a fundamental periodicity of P = 4.08 +/- 0.55 days along with a number of harmonics at periods P/n, with n ~ 2,3,4,5 and 6, and a few other possible stray periodicities and/or stochastic variability on timescales longer than about a day. Spectroscopic RV studies fail to reveal any plausible companion with a period in this range. Therefore, we conjecture that the observed light-curve cusps of amplitude ~ 0.01 mag that recur at a 4.08 day timescale may arise in the inner parts, or at the base of, a corotating interaction region (CIR) seen in emission as it rotates around with the star at constant angular velocity. The hard X-ray component seen in WR 110 could then be a result of a high velocity component of the CIR shock interacting with the ambient wind at several stellar radii. Given that most hot, luminous stars showing CIRs have two CIR arms, it is possible that either the fundamental period is 8.2 days or, more likely in the case of WR 110, there is indeed a second weaker CIR arm for P = 4.08 days, that occurs ~ two thirds of a rotation period after the main CIR. If this interpretation is correct, WR 110 therefore joins the ranks with three other single WR stars, all WN, with confirmed CIR rotation periods (WR 1, WR 6, and WR 134), albeit with WR 110 having by far the lowest amplitude photometric modulation. This illustrates the power of being able to secure intense, continuous high-precision photometry from space-based platforms such as MOST. It also opens the door to revealing low-amplitude photometric variations in other WN stars, where previous attempts have failed. If all WN stars have CIRs at some level, this could be important for revealing sources of magnetism or pulsation in addition to rotation periods.

Aims. What is the origin of the large-amplitude variability in Wolf-Rayet WN8 stars in general and WR123 in particular? A dedicated spectroscopic campaign targets the ten-hour period previously found in the high-precision photometric data obtained by the MOST satellite. Methods. In June-August 2003 we obtained a series of high signal-to-noise, mid-resolution spectra from several sites in the {\lambda}{\lambda} 4000 - 6940 A^{\circ} domain. We also followed the star with occasional broadband (Johnson V) photometry. The acquired spectroscopy allowed a detailed study of spectral variability on timescales from \sim 5 minutes to months. Results. We find that all observed spectral lines of a given chemical element tend to show similar variations and that there is a good correlation between the lines of different elements, without any significant time delays, save the strong absorption components of the Hei lines, which tend to vary differently from the emission parts. We find a single sustained periodicity, P \sim 9.8 h, which is likely related to the relatively stable pulsations found in MOST photometry obtained one year later. In addition, seemingly stochastic, large-amplitude variations are also seen in all spectral lines on timescales of several hours to several days.

The Wolf-Rayet star WR 46 is known to exhibit a very complex variability pattern on relatively short time scales of a few hours. Periodic but intermittent radial velocity shifts of optical lines as well as multiple photometric periods have been found in the past. Non-radial pulsations, rapid rotational modulation or the presence of a putative low-mass companion have been proposed to explain the short-term behaviour. In an effort to unveil its true nature, we observed WR 46 with FUSE (Far Ultraviolet Spectroscopic Explorer) over several short-term variability cycles. We found significant variations on a time scale of ~8 hours in the far-ultraviolet (FUV) continuum, in the blue edge of the absorption trough of the OVI {\lambda}{\lambda}1032, 1038 doublet P Cygni profile and in the SVI {\lambda}{\lambda}933, 944 P Cygni absorption profile. We complemented these observations with X-ray and UV light-curves and an X-ray spectrum from archival XMM-Newton (X-ray Multi-Mirror Mission - Newton Space Telescope) data. The X-ray and UV light-curves show variations on a time scale similar to the variability found in the FUV. We discuss our results in the context of the different scenarios suggested to explain the short-term variability of this object and reiterate that non-radial pulsations is the most likely to occur.

We present the results of an intensive photometric and spectroscopic monitoring campaign of the WN4 Wolf-Rayet (WR) star WR1=HD4004. Our broadband V photometry covering a timespan of 91 days shows variability with a period of P=16.9$^{+0.6}_{-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 estimate 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 WR1 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 WR137 is of the same nature. Finally, assuming that the period of CIRs is related to the rotational period, we estimate the rotation rate of the four stars for which sufficient monitoring has been carried out; i.e. v$_{rot}$=6.5, 40, 70 and 275 km/s for WR1, WR6, WR134 and WR137, respectively.

The Onfp class of rotationally broadened, hot spectra was defined some time ago in the Galaxy, where its membership to date numbers only eight. The principal defining characteristic is a broad, centrally reversed He II $\lambda$4686 emission profile; other emission and absorption lines are also rotationally broadened. Recent surveys in the Magellanic Clouds (MCs) have brought the class membership there, including some related spectra, to 28. We present a survey of the spectral morphology and rotational velocities, as a first step toward elucidating the nature of this class. Evolved, rapidly rotating hot stars are not expected theoretically, because the stellar winds should brake the rotation. Luminosity classification of these spectra is not possible, because the principal criterion (He II $\lambda$4686) is peculiar; however, the MCs provide reliable absolute magnitudes, which show that they span the entire range from dwarfs to supergiants. The Onfp line-broadening distribution is distinct and shifted toward larger values from those of normal O dwarfs and supergiants with >99.99% confidence. All cases with multiple observations show line-profile variations, which even remove some objects from the class temporarily. Some of them are spectroscopic binaries; it is possible that the peculiar profiles may have multiple causes among different objects. The origin and future of these stars are intriguing; for instance, they could be stellar mergers and/or gamma-ray-burst progenitors.

We report the detection of a strong, organised magnetic field in the O9IV star HD 57682, using spectropolarimetric observations obtained with ESPaDOnS at the 3.6-m Canada-France-Hawaii Telescope within the context of the Magnetism in Massive Stars (MiMeS) Large Program. From the fitting of our spectra using NLTE model atmospheres we determined that HD 57682 is a $17^{+19}_{-9}$ M$_{\odot}$ star with a radius of $7.0^{+2.4}_{-1.8}$ R$_\odot$, and a relatively low mass-loss rate of $1.4^{+3.1}_{-0.95}\times10^{-9}$ M$_{\odot}$ yr$^{-1}$. The photospheric absorption lines are narrow, and we use the Fourier transform technique to infer $v\sin i=15\pm3$ km s$^{-1}$. This $v\sin i$ implies a maximum rotational period of 31.5 d, a value qualitatively consistent with the observed variability of the optical absorption and emission lines, as well as the Stokes $V$ profiles and longitudinal field. Using a Bayesian analysis of the velocity-resolved Stokes $V$ profiles to infer the magnetic field characteristics, we tentatively derive a dipole field strength of $1680^{+134}_{-356}$ G. The derived field strength and wind characteristics imply a wind that is strongly confined by the magnetic field.

We present the results of a systematic search for large-scale spectroscopic variability in apparently single Wolf-Rayet stars brighter than ~12.5. In this first paper we characterize the various forms of variability detected and distinguish several separate groups. For each star in our sample, we obtained 4-5 high-resolution spectra with a signal-to-noise ratio ~100. Our ultimate goal is to identify new candidates presenting variability that potentially comes from Co-rotating Interaction Regions (CIR). Out of a sample of 25 stars, 10 were found to display large-scale changes of which 4 are of CIR-type (WR1, WR115, WR120 and WR134). The star WR134 was already known to show such changes from previous studies. Three WN8 stars present a different type of large-scale variability and we believe deserve a group of their own. Also, all three WC9d stars in our sample present large-scale variability, but it remains to be checked if these are binaries, as many dust-making WR stars are double. Finally, of the remaining stars, 10 were found to show small-amplitude spectral changes which we attribute to normal line-profile variability due to inhomogeneities in the wind, and 5 were found to show no spectral variability, as far as can be concluded from the data in hand. Follow-up studies are required to identify potential periods for our candidates showing CIR-type changes and eventually estimate a rotation rate for these WR stars.

Using the Very Large Telescope's Spectrograph for INtegral Field Observation in the Near-Infrared (VLT/SINFONI), we have obtained repeated AO-assisted, NIR spectroscopy of the six central luminous, Wolf-Rayet (WR) stars in the core of the very young (~1 Myr), massive and dense cluster R136, in the Large Magellanic Cloud (LMC). We also de-archived available images that were obtained with the Hubble Space Telescope's Space Telescope Imaging Spectrograph (HST/STIS), and extracted high-quality, differential photometry of our target stars to check for any variability related to binary motion. Previous studies, relying on spatially unresolved, integrated, optical spectroscopy, had reported that one of these stars was likely to be a 4.377-day binary. Our study set out to identify the culprit and any other short-period system among our targets. However, none displays significant photometric variability, and only one star, BAT99-112 (R136c), located on the outer fringe of R136, displays a marginal variability in its radial velocities; we tentatively report an 8.2-day period. The binary status of BAT99-112 is supported by the fact that it is one of the brightest X-ray sources among all known WR stars in the LMC, consistent with it being a colliding-wind system. Follow-up observations have been proposed to confirm the orbital period of this potentially very massive system.

We report the results of a spectroscopic and polarimetric study of the massive, hydrogen-rich WN6h stars R144 (HD 38282 = BAT99-118 = Brey 89) and R145 (HDE 269928 = BAT99-119 = Brey 90) in the LMC. Both stars have been suspected to be binaries by previous studies (R144: Schnurr et al. 2008b; R145: Moffat 1989). We have combined radial-velocity (RV) data from these two studies with previously unpublished polarimetric data. For R145, we were able to establish, for the first time, an orbital period of 158.8 days, along with the full set of orbital parameters, including the inclination angle i, which was found to be i = 38 \pm 9 deg. By applying a modified version of the shift-and-add method developed by Demers et al. (2002), we were able to isolate the spectral signature of the very faint-line companion star. With the RV amplitudes of both components in R145, we were thus able to estimate their absolute masses. We find minimum masses M_WR sin^{3}i = (116 \pm 33) M_sol and M_O sin^{3}i = (48 \pm 20)$ M_sol for the WR and the O component, respectively. Thus, if the low inclination angle were correct, resulting absolute masses of the components would be at least 300 and 125 M_sol, respectively. However, such high masses are not supported by brightness considerations when R145 is compared to systems with known, very high masses such as NGC3603-A1 or WR20a. An inclination angle close to 90 degrees would remedy the situation, but is excluded by the currently available data. More and better data are thus required to firmly establish the nature of this puzzling, yet potentially very massive and important system. As to R144, however, the combined data sets are not sufficient to find any periodicity.

We report the results of an intense, spectroscopic survey of all 41 late-type, nitrogen-rich Wolf-Rayet (WR) stars in the Large Magellanic Cloud (LMC) observable with ground-based telescopes. This survey concludes the decade-long effort of the Montr\'eal Massive Star Group to monitor every known WR star in the Magellanic Clouds except for the 6 crowded WNL stars in R136, which will be discussed elsewhere. The focus of our survey was to monitor the so-called WNL stars for radial-velocity (RV) variability in order to identify the short- to intermediate-period ($P \la 200$ days) binaries among them. Our results are in line with results of previous studies of other WR subtypes, and show that the binary frequency among LMC WNL stars is statistically consistent with that of WNL stars in the Milky Way. We have identified four previously unknown binaries, bringing the total number of known WNL binaries in the LMC to nine. Since it is very likely that none but one of the binaries are classical, helium-burning WNL stars, but rather superluminous, hence extremely massive, hydrogen-burning objects, our study has dramatically increased the number of known binaries harbouring such objects, and thus paved the way to determine their masses through model-independent, Keplerian orbits. It is expected that some of the stars in our binaries will be among the most massive known. With the binary status of each WR star now known, we also studied the photometric and X-ray properties of our program stars using archival MACHO photometry as well as Chandra and ROSAT data. We find that one of our presumably single WNL stars is among the X-ray brightest WR sources known. We also identify a binary candidate from its RV variability and X-ray luminosity which harbours the most luminous WR star known in the Local Group.

Using VLT/SINFONI, we have obtained repeated AO-assisted, NIR spectroscopy of the three central WN6ha stars in the core of the very young (~1 Myr), massive and dense Galactic cluster NGC3603. One of these stars, NGC3603-A1, is a known 3.77-day, double-eclipsing binary, while another one, NGC3603-C, is one of the brightest X-ray sources among all known Galactic WR stars, which usually is a strong indication for binarity. Our study reveals that star C is indeed an 8.9-day binary, although only the WN6ha component is visible in our spectra; therefore we temporarily classify star C as an SB1 system. A1, on the other hand, is found to consist of two emission-line stars of similar, but not necessarily of identical spectral type, which can be followed over most the orbit. Using radial velocities for both components and the previously known inclination angle of the system, we are able to derive absolute masses for both stars in A1. We find M_1 = (116 \pm 31) Mo for the primary and M_2 = (89 \pm 16) Mo for the secondary component of A1. While uncertainties are large, A1 is intrinsically half a magnitude brighter than WR20a, the current record holder with 83 and 82 Mo, respectively; therefore, it is likely that the primary in A1 is indeed the most massive star weighed so far.

We report the detection with FUSE of strong, highly blue-shifted absorption features appearing in the absorption troughs of practically all major P Cygni profiles in the presumably single Wolf-Rayet star WR135. These features also appear in the shock-sensitive OVI 1032/38 \AA doublet, coincident both in time and velocity space with the rest of the lower-ionization species. Choosing between two alternative interpretations: large-scale, coherent structures vs. localized, random shocks, we favor the latter. The absolute value of velocity, as well as velocity dispersion in the shocked region, the density of the shocked gas, and the time scales of the observed variability allow us to relate the observed shocks to the incidence of numerous over-dense clumps (blobs) in the wind of a hot, massive star.

CX Cep (WR 151) is the WR+O binary (WN5+O5V) with the second shortest period known in our Galaxy. To examine the circumstellar matter distribution and to better constraint the orbital parameters and mass-loss rate of the WR star, we obtained broadband and multi-band (i.e. UBVRI) linear polarization observations of the system. Our analysis of the phase-locked polarimetric modulation confirms the high orbital inclination of the system (i.e. $i=65^o$). Using the orbital solution of Lewis et al. (1993) we obtain masses of $33.9 M_{\odot}$ and $23.9 M_{\odot}$ for the O and WR stars respectively, which agree with their spectral types. A simple polarimetric model accounting for finite stellar size effects allowed us to derive a mass-loss rate for the WR star of $0.3-0.5\times10^{-5} M_{\odot}/yr$. This result was remarkably independent of the model's input parameters and favors an earlier spectral type for the WR component (i.e. WN4). Finally, using our multi-band observations, we fitted and subtracted from our data the interstellar polarization. The resulting constant intrinsic polarization of $3-4%$ is misaligned in relation to the orbital plane (i.e. $\theta_{CIP}=26^o$ vs. $\Omega=75^o$) and is the highest intrinsic polarization ever observed for a WR star. This misalignment points towards a rotational (or magnetic) origin for the asymmetry and contradicts the most recent evolutionary models for massive stars (Meynet & Maeder 2003) which predict spherically symmetric winds during the WR phase (i.e. $CIP=0%$).

In an effort to better calibrate the supernova rate of starburst galaxies as determined from near-IR [Fe II] features, we report on a [Fe II] 1.644 microns line-imaging survey of a sample of 42 optically-selected SNRs in M33. A wide range of [Fe II] luminosities are observed within our sample (from less than 6 to 695 L_sun). Our data suggest that the bright [Fe II] SNRs are entering the radiative phase and that the density of the local ISM largely controls the amount of [Fe II] emission. We derive the following relation between the [Fe II] 1.644 microns line luminosity of radiative SNRs and the electronic density of the postshock gas, n_e: L_[Fe II] (L_sun) ~ 1.1 n_e (cm^-3). We also find a correlation in our data between L_[Fe II] and the metallicity of the shock-heated gas, but the physical interpretation of this result remains inconclusive, as our data also show a correlation between the metallicity and n_e. The dramatically higher level of [Fe II] emission from SNRs in the central regions of starburst galaxies is most likely due to their dense environments, although metallicity effects might also be important. The typical [Fe II]-emitting lifetime of a SNR in the central regions of starburst galaxies is found to be of the order of 10^4 yr. On the basis of these results, we provide a new empirical relation allowing the determination of the current supernova rate of starburst galaxies from their integrated near-IR [Fe II] luminosity.

A new high-quality set of orbital parameters for the O-type spectroscopic binary HD93205 has been obtained combining echelle and coude CCD observations. The radial velocity orbits derived from the HeII 4686 A (primary component) and HeI 4471 A (secondary component) absorption lines yield semiamplitudes of 133 +/- 2 and 314 +/- 2 km/s for each binary component, resulting in minimum masses of 31 and 13 Msol (q = 0.42). We also confirm for the binary components the spectral classification of O3V+O8V previously assigned. Assuming for the O8V component a ``normal'' mass of 22 - 25 Msol we would derive for the primary O3V a mass of ``only'' 52 - 60 Msol and an inclination of about 55 degrees for the orbital plane. We have also determined for the first time a period of apsidal motion for this system, namely 185 +/- 16 years using all available radial velocity data-sets of HD 93205 (from 1975 to 1999). Phase-locked variations of the X-ray emission of HD 93205 consisting of a rise of the observed X-ray flux near periastron passage, are also discussed.

A numerical code aiming at the automatic detection of spherical expanding HI shells in radio data-cubes is presented. The following five shell parameters are allowed to vary within a specified range: angular radius, expansion velocity and three center coordinates (galactic coordinates l and b, and systemic radial velocity V). We discuss several factors which can reduce the sensitivity of the shell detection: the presence of noise (both instrumental and "structural"), the fragmentation and asphericity of the shell and the inhomogeneity of the background and/or foreground emission. The code is tested on four objects: two HII regions and two early B stars. We use HI data from the Canadian Galactic Plane Survey. In all four cases the evidence for an expanding HI shell is found.

We report on the first near-IR [Fe II] line imaging survey of extragalactic SNRs. Observations of a sample of 42 objects drawn from an optically-selected catalogue of SNRs in M33 provide evidence for a wide range in the [Fe II] 1.644 microns luminosities. This can be understood as being primarily due to variations in the chemical abundances and density of the local ISM, although shock conditions may also play a significant role. We briefly discuss how these results may be used to better calibrate the supernova rate of star-forming galaxies.

In recent years, much studies have focused on determining the origin of the large-scale line-profile and/or photometric patterns of variability displayed by some apparently single Wolf-Rayet stars, with the existence of an unseen (collapsed?) companion or of spatially extended wind structures as potential candidates. We present observations of WR 1 which highlight the unusual character of the variations in this object. Our narrowband photometric observations reveal a gradual increase of the stellar continuum flux amounting to Delta v = 0.09 mag followed by a decline on about the same timescale (3-4 days). Only marginal evidence for variability is found during the 11 following nights. Strong, daily line-profile variations are also observed but they cannot be easily linked to the photometric variations. Similarly to the continuum flux variations, coherent time-dependent changes are observed in 1996 in the centroid, equivalent width, and skewness of He II 4686. Despite the generally coherent nature of the variations, we do not find evidence in our data for the periods claimed in previous studies. While the issue of a cyclical pattern of variability in WR 1 is still controversial, it is clear that this object might constitute in the future a cornerstone for our understanding of the mechanisms leading to the formation of largely anisotropic outflows in Wolf-Rayet stars.

We present the results of an intensive campaign of spectroscopic and photometric monitoring of the peculiar Wolf-Rayet star WR 134 from 1989 to 1997. This unprecedentedly large data set allows us to confirm unambiguously the existence of a coherent 2.25 +/- 0.05 day periodicity in the line-profile changes of He II 4686, although the global pattern of variability is different from one epoch to another. This period is only marginally detected in the photometric data set. Assuming the 2.25 day periodic variability to be induced by orbital motion of a collapsed companion, we develop a simple model aiming at investigating (i) the effect of this strongly ionizing, accreting companion on the Wolf-Rayet wind structure, and (ii) the expected emergent X-ray luminosity. We argue that the predicted and observed X-ray fluxes can only be matched if the accretion on the collapsed star is significantly inhibited. Additionally, we performed simulations of line-profile variations caused by the orbital revolution of a localized, strongly ionized wind cavity surrounding the X-ray source. A reasonable fit is achieved between the observed and modeled phase-dependent line profiles of He II 4686. However, the derived size of the photoionized zone substantially exceeds our expectations, given the observed low-level X-ray flux. Alternatively, we explore rotational modulation of a persistent, largely anisotropic outflow as the origin of the observed cyclical variability. Although qualitative, this hypothesis leads to greater consistency with the observations.