@article{AldorettaStLouisRichardsonetal.2016,
  author    = {Aldoretta, E. J. and St-Louis, N. and Richardson, N. D. and Moffat, Anthony F. J. and Eversberg, T. and Hill, G. M. and Shenar, Tomer and Artigau, E. and Gauza, B. and Knapen, J. H. and Kubat, Jiř{\´i} and Kubatova, Brankica and Maltais-Tariant, R. and Munoz, M. and Pablo, H. and Ramiaramanantsoa, T. and Richard-Laferriere, A. and Sablowski, D. P. and Simon-Diaz, S. and St-Jean, L. and Bolduan, F. and Dias, F. M. and Dubreuil, P. and Fuchs, D. and Garrel, T. and Grutzeck, G. and Hunger, T. and Kuesters, D. and Langenbrink, M. and Leadbeater, R. and Li, D. and Lopez, A. and Mauclaire, B. and Moldenhawer, T. and Potter, M. and dos Santos, E. M. and Schanne, L. and Schmidt, J. and Sieske, H. and Strachan, J. and Stinner, E. and Stinner, P. and Stober, B. and Strandbaek, K. and Syder, T. and Verilhac, D. and Waldschlaeger, U. and Weiss, D. and Wendt, A.},
  title     = {An extensive spectroscopic time series of three Wolf-Rayet stars - I. The lifetime of large-scale structures in the wind of WR 134},
  series = {Monthly notices of the Royal Astronomical Society},
  volume    = {460},
  journal   = {Monthly notices of the Royal Astronomical Society},
  publisher = {Oxford Univ. Press},
  address   = {Oxford},
  issn      = {0035-8711},
  doi       = {10.1093/mnras/stw1188},
  pages     = {3407 -- 3417},
  year      = {2016},
  abstract  = {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 analysed 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 analysing the variability of the He ii lambda 5411 emission line, the previously identified period was refined to P = 2.255 +/- 0.008 (s.d.) d. The coherency time of the variability, which we associate with the lifetime of the CIRs in the wind, was deduced to be 40 +/- 6 d, or similar to 18 cycles, by cross-correlating the variability patterns as a function of time. When comparing the phased observational grey-scale difference images with theoretical grey-scales 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 I center dot a parts per thousand integral 90A degrees 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.},
  language  = {en}
}
@article{CorcoranNicholsPabloetal.2015,
  author    = {Corcoran, Michael F. and Nichols, Joy S. and Pablo, Herbert and Shenar, Tomer and Pollock, Andy M. T. and Waldron, Wayne L. and Moffat, Anthony F. J. and Richardson, Noel D. and Russell, Christopher M. P. and Hamaguchi, Kenji and Huenemoerder, David P. and Oskinova, Lida and Hamann, Wolf-Rainer and Naze, Yael and Ignace, Richard and Evans, Nancy Remage and Lomax, Jamie R. and Hoffman, Jennifer L. and Gayley, Kenneth and Owocki, Stanley P. and Leutenegger, Maurice and Gull, Theodore R. and Hole, Karen Tabetha and Lauer, Jennifer and Iping, Rosina C.},
  title     = {A coordinated X-Ray and optical campaign of the nearest massive eclipsing binary, delta ORIONIS Aa. I. Overview of thr X-Ray spectrum},
  series = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  volume    = {809},
  journal   = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  number    = {2},
  publisher = {IOP Publ. Ltd.},
  address   = {Bristol},
  issn      = {0004-637X},
  doi       = {10.1088/0004-637X/809/2/132},
  pages     = {15},
  year      = {2015},
  abstract  = {We present an overview of four deep phase-constrained Chandra HETGS X-ray observations of delta Ori A. Delta Ori A is actually a triple system that includes the nearest massive eclipsing spectroscopic binary, delta Ori Aa, the only such object that can be observed with little phase-smearing with the Chandra gratings. Since the fainter star, delta Ori Aa2, has a much lower X-ray luminosity than the brighter primary (delta Ori Aa1), delta Ori Aa provides a unique system with which to test the spatial distribution of the X-ray emitting gas around delta Ori Aa1 via occultation by the photosphere of, and wind cavity around, the X-ray dark secondary. Here we discuss the X-ray spectrum and X-ray line profiles for the combined observation, having an exposure time of nearly 500 ks and covering nearly the entire binary orbit. The companion papers discuss the X-ray variability seen in the Chandra spectra, present new space-based photometry and ground-based radial velocities obtained simultaneously with the X-ray data to better constrain the system parameters, and model the effects of X-rays on the optical and UV spectra. We find that the X-ray emission is dominated by embedded wind shock emission from star Aa1, with little contribution from the tertiary star Ab or the shocked gas produced by the collision of the wind of Aa1 against the surface of Aa2. We find a similar temperature distribution to previous X-ray spectrum analyses. We also show that the line half-widths are about 0.3-0.5 times the terminal velocity of the wind of star Aa1. We find a strong anti-correlation between line widths and the line excitation energy, which suggests that longer-wavelength, lower-temperature lines form farther out in the wind. Our analysis also indicates that the ratio of the intensities of the strong and weak lines of Fe XVII and Ne X are inconsistent with model predictions, which may be an effect of resonance scattering.},
  language  = {en}
}
@article{EvansDeGioiaEastwoodGagneetal.2011,
  author    = {Evans, Nancy Remage and DeGioia-Eastwood, Kathleen and Gagne, Marc and Townsley, Leisa and Broos, Patrick S. and Wolk, Scott J. and Naze, Yael and Corcoran, Michael and Oskinova, Lida and Moffat, Anthony F. J. and Wang, Junfeng and Walborn, Nolan R.},
  title     = {The search for low-mass companions of b stars in the Carina Nebula cluster trumpler 16},
  series = {The astrophysical journal : an international review of spectroscopy and astronomical physics ; Supplement series},
  volume    = {194},
  journal   = {The astrophysical journal : an international review of spectroscopy and astronomical physics ; Supplement series},
  number    = {1},
  publisher = {IOP Publ. Ltd.},
  address   = {Bristol},
  issn      = {0067-0049},
  doi       = {10.1088/0067-0049/194/1/13},
  pages     = {9},
  year      = {2011},
  abstract  = {We have developed lists of likely B3-A0 stars (called "late B" stars) in the young cluster Trumpler 16. The following criteria were used: location within 3' of eta Car, an appropriate V and B - V combination, and proper motion (where available). Color and magnitude cuts have been made assuming an E(B - V) = 0.55 mag +/- 0.1, which is a good approximation close to the center of Trumpler 16. These lists have been cross-correlated with X-ray sources found in the Chandra Carina Complex Project. Previous studies have shown that only very rarely (if at all) do late main-sequence B stars produce X-rays. We present evidence that the X-ray-detected sources are binaries with low-mass companions, since stars less massive than 1.4 M-circle dot are strong X-ray sources at the age of the cluster. Both the median X-ray energies and X-ray luminosities of these sources are in good agreement with values for typical low-mass coronal X-ray sources. We find that 39\% of the late B stars based on a list with proper motions have low-mass companions. Similarly, 32\% of a sample without proper motions have low-mass companions. We discuss the X-ray detection completeness. These results on low-mass companions of intermediate-mass stars are complementary to spectroscopic and interferometric results and probe new parameter space of low-mass companions at all separations. They do not support a steeply rising distribution of mass ratios to low masses for intermediate-mass (5 M-circle dot) primaries, such as would be found by random pairing from the initial mass function.},
  language  = {en}
}
@article{GagneFehonSavoyetal.2011,
  author    = {Gagne, Marc and Fehon, Garrett and Savoy, Michael R. and Cohen, David H. and Townsley, Leisa K. and Broos, Patrick S. and Povich, Matthew S. and Corcoran, Michael F. and Walborn, Nolan R. and Evans, Nancy Remage and Moffat, Anthony F. J. and Naze, Yael and Oskinova, Lida},
  title     = {Carina ob stars: x-ray signatures of wind shocks and magnetic FIELDS},
  series = {The astrophysical journal : an international review of spectroscopy and astronomical physics ; Supplement series},
  volume    = {194},
  journal   = {The astrophysical journal : an international review of spectroscopy and astronomical physics ; Supplement series},
  number    = {1},
  publisher = {IOP Publ. Ltd.},
  address   = {Bristol},
  issn      = {0067-0049},
  doi       = {10.1088/0067-0049/194/1/5},
  pages     = {26},
  year      = {2011},
  abstract  = {The Chandra Carina Complex contains 200 known O- and B-type stars. The Chandra survey detected 68 of the 70 O stars and 61 of 127 known B0-B3 stars. We have assembled a publicly available optical/X-ray database to identify OB stars that depart from the canonical L-X/L-bol relation or whose average X-ray temperatures exceed 1 keV. Among the single O stars with high kT we identify two candidate magnetically confined wind shock sources: Tr16-22, O8.5 V, and LS 1865, O8.5 V((f)). The O4 III(fc) star HD 93250 exhibits strong, hard, variable X-rays, suggesting that it may be a massive binary with a period of > 30 days. The visual O2 If* binary HD 93129A shows soft 0.6 keV and hard 1.9 keV emission components, suggesting embedded wind shocks close to the O2 If* Aa primary and colliding wind shocks between Aa and Ab. Of the 11 known O-type spectroscopic binaries, the long orbital-period systems HD 93343, HD 93403, and QZ Car have higher shock temperatures than short-period systems such as HD 93205 and FO 15. Although the X-rays from most B stars may be produced in the coronae of unseen, low-mass pre-main-sequence companions, a dozen B stars with high L-X cannot be explained by a distribution of unseen companions. One of these, SS73 24 in the Treasure Chest cluster, is a new candidate Herbig Be star.},
  language  = {en}
}
@article{MunozMoffatHilletal.2017,
  author    = {Munoz, Melissa and Moffat, Anthony F. J. and Hill, Grant M. and Shenar, Tomer and Richardson, Noel D. and Pablo, Herbert and St-Louis, Nicole and Ramiaramanantsoa, Tahina},
  title     = {WR 148: identifying the companion of an extreme runaway massive binary},
  series = {Monthly notices of the Royal Astronomical Society},
  volume    = {467},
  journal   = {Monthly notices of the Royal Astronomical Society},
  publisher = {Oxford Univ. Press},
  address   = {Oxford},
  issn      = {0035-8711},
  doi       = {10.1093/mnras/stw2283},
  pages     = {3105 -- 3121},
  year      = {2017},
  abstract  = {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 ± 3.8\&\#8201;km\&\#8201;s\&\#8722;1 and 79.2 ± 3.1\&\#8201;km\&\#8201;s\&\#8722;1 and implying a mass ratio of 1.1 ± 0.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 \&\#8764;67°, derived from previous polarimetry observations, the system's total mass would be a mere 2-3M\&\#8857;\&\#8288;, an unprecedented result for a putative massive binary system. However, a system comprising a 37M\&\#8857; secondary (typical mass of an O5V star) and a 33M\&\#8857; primary (given the mass ratio) would infer an inclination of \&\#8764;18°. 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 ± 27\&\#8201;km\&\#8201;s\&\#8722;1 and a travel time of 4.7 ± 0.8 Myr to reach its current location, WR 148 was most likely ejected via dynamical interactions in a young cluster.},
  language  = {en}
}
@article{NazeBroosOskinovaetal.2011,
  author    = {Naze, Y. and Broos, Patrick S. and Oskinova, Lida and Townsley, L. K. and Cohen, David H. and Corcoran, M. F. and Evans, N. R. and Gagne, M. and Moffat, Anthony F. J. and Pittard, J. M. and Rauw, G. and Ud-Doula, A. and Walborn, N. R.},
  title     = {GLOBAL X-RAY PROPERTIES OF THE O AND B STARS IN CARINA},
  series = {ASTROPHYSICAL JOURNAL SUPPLEMENT SERIES},
  volume    = {194},
  journal   = {ASTROPHYSICAL JOURNAL SUPPLEMENT SERIES},
  number    = {1},
  publisher = {IOP PUBLISHING LTD},
  address   = {BRISTOL},
  issn      = {0067-0049},
  doi       = {10.1088/0067-0049/194/1/7},
  pages     = {20},
  year      = {2011},
  abstract  = {The key empirical property of the X-ray emission from O stars is a strong correlation between the bolometric and X-ray luminosities. In the framework of the Chandra Carina Complex Project, 129 O and B stars have been detected as X-ray sources; 78 of those, all with spectral type earlier than B3, have enough counts for at least a rough X-ray spectral characterization. This leads to an estimate of the L-X-L-BOL ratio for an exceptional number of 60 O stars belonging to the same region and triples the number of Carina massive stars studied spectroscopically in X-rays. The derived log(L-X/L-BOL) is -7.26 for single objects, with a dispersion of only 0.21 dex. Using the properties of hot massive stars listed in the literature, we compare the X-ray luminosities of different types of objects. In the case of O stars, the L-X-L-BOL ratios are similar for bright and faint objects, as well as for stars of different luminosity classes or spectral types. Binaries appear only slightly harder and slightly more luminous in X-rays than single objects; the differences are not formally significant (at the 1\% level), except for the L-X-L-BOL ratio in the medium (1.0-2.5 keV) energy band. Weak-wind objects have similar X-ray luminosities but they display slightly softer spectra compared with "normal" O stars with the same bolometric luminosity. Discarding three overluminous objects, we find a very shallow trend of harder emission in brighter objects. The properties of the few B stars bright enough to yield some spectral information appear to be different overall (constant X-ray luminosities, harder spectra), hinting that another mechanism for producing X-rays, besides wind shocks, might be at work. However, it must be stressed that the earliest and X-ray brightest among these few detected objects are similar to the latest O stars, suggesting a possibly smooth transition between the two processes.},
  language  = {en}
}
@article{NicholsHuenemoerderCorcoranetal.2015,
  author    = {Nichols, Joy and Huenemoerder, David P. and Corcoran, Michael F. and Waldron, Wayne and Naze, Yael and Pollock, Andy M. T. and Moffat, Anthony F. J. and Lauer, Jennifer and Shenar, Tomer and Russell, Christopher M. P. and Richardson, Noel D. and Pablo, Herbert and Evans, Nancy Remage and Hamaguchi, Kenji and Gull, Theodore and Hamann, Wolf-Rainer and Oskinova, Lida and Ignace, Rosina and Hoffman, Jennifer L. and Hole, Karen Tabetha and Lomax, Jamie R.},
  title     = {A COORDINATED X-RAY AND OPTICAL CAMPAIGN OF THE NEAREST MASSIVE ECLIPSING BINARY, delta ORIONIS Aa. II. X-RAY VARIABILITY},
  series = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  volume    = {809},
  journal   = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  number    = {2},
  publisher = {IOP Publ. Ltd.},
  address   = {Bristol},
  issn      = {0004-637X},
  doi       = {10.1088/0004-637X/809/2/133},
  pages     = {21},
  year      = {2015},
  abstract  = {We present time-resolved and phase-resolved variability studies of an extensive X-ray high-resolution spectral data set of the delta Ori Aa binary system. The four observations, obtained with Chandra ACIS HETGS, have a total exposure time of approximate to 479 ks and provide nearly complete binary phase coverage. Variability of the total X-ray flux in the range of 5-25 is is confirmed, with a maximum amplitude of about +/- 15\% within a single approximate to 125 ks observation. Periods of 4.76 and 2.04 days are found in the total X-ray flux, as well as an apparent overall increase in the flux level throughout the nine-day observational campaign. Using 40 ks contiguous spectra derived from the original observations, we investigate the variability of emission line parameters and ratios. Several emission lines are shown to be variable, including S XV, Si XIII, and Ne IX. For the first time, variations of the X-ray emission line widths as a function of the binary phase are found in a binary system, with the smallest widths at phi = 0.0 when the secondary delta Ori Aa2 is at the inferior conjunction. Using 3D hydrodynamic modeling of the interacting winds, we relate the emission line width variability to the presence of a wind cavity created by a wind-wind collision, which is effectively void of embedded wind shocks and is carved out of the X-ray-producing primary wind, thus producing phase-locked X-ray variability.},
  language  = {en}
}
@article{PabloRichardsonMoffatetal.2015,
  author    = {Pablo, Herbert and Richardson, Noel D. and Moffat, Anthony F. J. and Corcoran, Michael and Shenar, Tomer and Benvenuto, Omar and Fuller, Jim and Naze, Yael and Hoffman, Jennifer L. and Miroshnichenko, Anatoly and Apellaniz, Jesus Maiz and Evans, Nancy and Eversberg, Thomas and Gayley, Ken and Gull, Ted and Hamaguchi, Kenji and Hamann, Wolf-Rainer and Henrichs, Huib and Hole, Tabetha and Ignace, Richard and Iping, Rosina and Lauer, Jennifer and Leutenegger, Maurice and Lomax, Jamie and Nichols, Joy and Oskinova, Lida and Owocki, Stan and Pollock, Andy and Russell, Christopher M. P. and Waldron, Wayne and Buil, Christian and Garrel, Thierry and Graham, Keith and Heathcote, Bernard and Lemoult, Thierry and Li, Dong and Mauclaire, Benjamin and Potter, Mike and Ribeiro, Jose and Matthews, Jaymie and Cameron, Chris and Guenther, David and Kuschnig, Rainer and Rowe, Jason and Rucinski, Slavek and Sasselov, Dimitar and Weiss, Werner},
  title     = {A coordinated X-Ray and optical campaign of the nearest massive eclipsing binary, delta ORIONIS Aa. III. Analysis of optical photometric (most) and spectroscopic (ground based) variations},
  series = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  volume    = {809},
  journal   = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  number    = {2},
  publisher = {IOP Publ. Ltd.},
  address   = {Bristol},
  issn      = {0004-637X},
  doi       = {10.1088/0004-637X/809/2/134},
  pages     = {11},
  year      = {2015},
  abstract  = {We report on both high-precision photometry from the Microvariability and Oscillations of Stars (MOST) space telescope and ground-based spectroscopy of the triple system delta Ori A, consisting of a binary O9.5II+early-B (Aa1 and Aa2) with P = 5.7 days, and a more distant tertiary (O9 IV P > 400 years). This data was collected in concert with X-ray spectroscopy from the Chandra X-ray Observatory. Thanks to continuous coverage for three weeks, the MOST light curve reveals clear eclipses between Aa1 and Aa2 for the first time in non-phased data. From the spectroscopy, we have a well-constrained radial velocity (RV) curve of Aa1. While we are unable to recover RV variations of the secondary star, we are able to constrain several fundamental parameters of this system and determine an approximate mass of the primary using apsidal motion. We also detected second order modulations at 12 separate frequencies with spacings indicative of tidally influenced oscillations. These spacings have never been seen in a massive binary, making this system one of only a handful of such binaries that show evidence for tidally induced pulsations.},
  language  = {en}
}
@article{ParkinBroosTownsleyetal.2011,
  author    = {Parkin, E. R. and Broos, Patrick S. and Townsley, L. K. and Pittard, J. M. and Moffat, Anthony F. J. and Naze, Y. and Rauw, G. and Oskinova, Lida and Waldron, W. L.},
  title     = {X-RAY EMISSION FROM THE DOUBLE-BINARY OB-STAR SYSTEM QZ CAR (HD 93206)},
  series = {ASTROPHYSICAL JOURNAL SUPPLEMENT SERIES},
  volume    = {194},
  journal   = {ASTROPHYSICAL JOURNAL SUPPLEMENT SERIES},
  number    = {1},
  publisher = {IOP PUBLISHING LTD},
  address   = {BRISTOL},
  issn      = {0067-0049},
  doi       = {10.1088/0067-0049/194/1/8},
  pages     = {11},
  year      = {2011},
  abstract  = {X-ray observations of the double-binary OB-star system QZ Car (HD 93206) obtained with the Chandra X-ray Observatory over a period of roughly 2 years are presented. The respective orbits of systems A (O9.7 I+b2 v, P-A = 21 days) and B (O8 III+o9 v, P-B = 6 days) are reasonably well sampled by the observations, allowing the origin of the X-ray emission to be examined in detail. The X-ray spectra can be well fitted by an attenuated three-temperature thermal plasma model, characterized by cool, moderate, and hot plasma components at kT similar or equal to 0.2, 0.7, and 2 keV, respectively, and a circumstellar absorption of similar or equal to 0.2 x 10(22) cm(-2). Although the hot plasma component could be indicating the presence of wind-wind collision shocks in the system, the model fluxes calculated from spectral fits, with an average value of similar or equal to 7x10(-13) erg s(-1) cm(-2), do not show a clear correlation with the orbits of the two constituent binaries. A semi-analytical model of QZ Car reveals that a stable momentum balance may not be established in either system A or B. Yet, despite this, system B is expected to produce an observed X-ray flux well in excess of the observations. If one considers the wind of the O8 III star to be disrupted by mass transfer, the model and observations are in far better agreement, which lends support to the previous suggestion of mass transfer in the O8 III+o9 v binary. We conclude that the X-ray emission from QZ Car can be reasonably well accounted for by a combination of contributions mainly from the single stars and the mutual wind-wind collision between systems A and B.},
  language  = {en}
}
@article{RamiaramanantsoaMoffatHarmonetal.2018,
  author    = {Ramiaramanantsoa, Tahina and Moffat, Anthony F. J. and Harmon, Robert and Ignace, R. and St-Louis, Nicole and Vanbeveren, Dany and Shenar, Tomer and Pablo, Herbert and Richardson, Noel D. and Howarth, Ian D. and Stevens, Ian R. and Piaulet, Caroline and St-Jean, Lucas and Eversberg, Thomas and Pigulski, Andrzej and Popowicz, Adam and Kuschnig, Rainer and Zoclonska, Elzbieta and Buysschaert, Bram and Handler, Gerald and Weiss, Werner W. and Wade, Gregg A. and Rucinski, Slavek M. and Zwintz, Konstanze and Luckas, Paul and Heathcote, Bernard and Cacella, Paulo and Powles, Jonathan and Locke, Malcolm and Bohlsen, Terry and Chen{\´e}, Andr{\´e}-Nicolas and Miszalski, Brent and Waldron, Wayne L. and Kotze, Marissa M. and Kotze, Enrico J. and B{\"o}hm, Torsten},
  title     = {BRITE-Constellation high-precision time-dependent photometry of the early O-type supergiant zeta Puppis unveils the photospheric drivers of its small- and large-scale wind structures},
  series = {Monthly notices of the Royal Astronomical Society},
  volume    = {473},
  journal   = {Monthly notices of the Royal Astronomical Society},
  number    = {4},
  publisher = {Oxford Univ. Press},
  address   = {Oxford},
  issn      = {0035-8711},
  doi       = {10.1093/mnras/stx2671},
  pages     = {5532 -- 5569},
  year      = {2018},
  abstract  = {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 ζ 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/Solar Mass Ejection Imager, 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 multisite spectroscopic monitoring of the star unveiled cyclical modulation of its He ii λ4686 wind emission line with the 1.78-d rotation period, showing signatures of corotating interaction regions 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 λ4686 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 ζ 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.},
  language  = {en}
}