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.

Stochastically-triggered photospheric light variations reaching $\sim$$40$ mmag peak-to-valley amplitudes have been detected in the O8Iaf supergiant V973 Scorpii as the outcome of two months of high-precision time-resolved photometric observations with the BRIght Target Explorer (BRITE) nanosatellites. The amplitude spectrum of the time series photometry exhibits a pronounced broad bump in the low-frequency regime ($\lesssim$$0.9$ d$^{-1}$) where several prominent frequencies are detected. A time-frequency analysis of the observations reveals typical mode lifetimes of the order of $5-10$ days. The overall features of the observed brightness amplitude spectrum of V973 Sco match well with those extrapolated from two-dimensional hydrodynamical simulations of convectively-driven internal gravity waves randomly excited from deep in the convective cores of massive stars. An alternative or additional possible source of excitation from a subsurface convection zone needs to be explored in future theoretical investigations.

Previously, the X-ray emission of zeta Puppis was found to be variable with light curves harbouring "trends" with a typical timescale longer than the exposure length. The origin of these changes was proposed to be linked to large-scale structures in the wind, but further characterization of the variability at high energies was needed. Since then, a number of new X-ray observations have become available. Furthermore, a cyclic behaviour with a 1.78d period was identified in long optical photometric runs, which is thought to be associated with the launching mechanism of large-scale wind structures. We analysed these new X-ray data, revisited the old data, and compared X-ray with optical data, including when simultaneous. We found that the behaviour in X-rays cannot be explained in terms of a perfect clock because the amplitude and shape of its variations change with time. For example, zeta Puppis was much more strongly variable between 2007 and 2011 than before and after this interval. Comparing the X-ray spectra of the star at maximum and minimum brightness yields no compelling difference beyond the overall flux change: the temperatures, absorptions, and line shapes seem to remain constant, well within errors. The only common feature between X-ray datasets is that the variation amplitudes appear maximum in the medium (0.6-1.2keV) energy band. Finally, no clear and coherent correlation can be found between simultaneous X-ray and optical data. Only a subgroup of observations may be combined coherently with the optical period of 1.78d, although the simultaneous optical behaviour is unknown. The currently available data do not reveal any obvious, permanent, and direct correlation between X-ray and optical variations. The origin of the X-ray variability therefore still needs to be ascertained, highlighting the need for long-term monitoring in multiwavelengths, i.e. X-ray, UV, and optical.

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.

The BRITE mission is a pioneering space project aimed at the long-term photometric monitoring of the brightest stars in the sky by means of a constellation of nano-satellites. Its main advantage is high photometric accuracy and time coverage inaccessible from the ground. The main aim of this paper is the presentation of procedures used to obtain high-precision photometry from a series of images acquired by the BRITE satellites in two modes of observing, stare and chopping. We developed two pipelines corresponding to the two modes of observing. The assessment of the performance of both pipelines is presented. It is based on two comparisons, which use data from six runs of the UniBRITE satellite: (i) comparison of photometry obtained by both pipelines on the same data, which were partly affected by charge transfer inefficiency (CTI), (ii) comparison of real scatter with theoretical expectations. It is shown that for CTI-affected observations, the chopping pipeline provides much better photometry than the other pipeline. For other observations, the results are comparable only for data obtained shortly after switching to chopping mode. Starting from about 2.5 years in orbit, the chopping mode of observing provides significantly better photometry for UniBRITE data than the stare mode. This paper shows that high-precision space photometry with low-cost nano-satellites is achievable. The proposed meth- ods, used to obtain photometry from images affected by high impulsive noise, can be applied to data from other space missions or even to data acquired from ground-based observations.

Massive stars play a significant role in the chemical and dynamical evolution of galaxies. However, much of their variability, particularly during their evolved supergiant stage, is poorly understood. To understand the variability of evolved massive stars in more detail, we present a study of the O9.2Ib supergiant $\zeta$ Ori Aa, the only currently confirmed supergiant to host a magnetic field. We have obtained two-color space-based BRIght Target Explorer photometry (BRITE) for $\zeta$ Ori Aa during two observing campaigns, as well as simultaneous ground-based, high-resolution optical CHIRON spectroscopy. We perform a detailed frequency analysis to detect and characterize the star's periodic variability. We detect two significant, independent frequencies, their higher harmonics, and combination frequencies: the stellar rotation period $P_{\mathrm{rot}} = 6.82\pm0.18$ d, most likely related to the presence of the stable magnetic poles, and a variation with a period of $10.0\pm0.3$ d attributed to circumstellar environment, also detected in the H$\alpha$ and several He I lines, yet absent in the purely photospheric lines. We confirm the variability with $P_{\mathrm{rot}}$/4, likely caused by surface inhomogeneities, being the possible photospheric drivers of the discrete absorption components. No stellar pulsations were detected in the data. The level of circumstellar activity clearly differs between the two BRITE observing campaigns. We demonstrate that $\zeta$ Ori Aa is a highly variable star with both periodic and non-periodic variations, as well as episodic events. The rotation period we determined agrees well with the spectropolarimetric value from the literature. The changing activity level observed with BRITE could explain why the rotational modulation of the magnetic measurements was not clearly detected at all epochs.

{\iota} Ori is a well studied massive binary consisting of an O9 III + B1 III/IV star. Due to its high eccentricity (e = 0.764) and short orbital period (P orb = 29.13376 d), it has been considered to be a good candidate to show evidence of tidal effects; however, none have previously been identified. Using photometry from the BRITE-Constellation space photometry mission we have confirmed the existence of tidal distortions through the presence of a heartbeat signal at periastron. We combine spectroscopic and light curve analyses to measure the masses and radii of the components, revealing {\iota} Ori to be the most massive heartbeat system known to date. In addition, using a thorough frequency analysis, we also report the unprecedented discovery of multiple tidally induced oscillations in an O star. The amplitudes of the pulsations allow us to empirically estimate the tidal circularization rate, yielding an effective tidal quality factor Q $\approx 4 \times 10^{4}$ .

We studied the variability of the magnetic O-type supergiant $\zeta$ Ori Aa using multi-colour BRITE photometry. We confirmed the known rotation frequency $f_{\rm rot} = 0.15 \pm 0.02$ c/d, and detected some of its higher harmonics, of which $4f_{\rm rot}$ is compatible with the known DAC recurrence timescale. Thanks to simultaneous high-resolution CHIRON spectroscopy, we could identify another frequency $f_{\rm env} = 0.10 \pm 0.02$ c/d, caused by the circumstellar environment. Variations in the circumstellar environment are believed to cause the observed difference between the BRITE lightcurves.

We report a simultaneous ground and space-based photometric study of the Beta Cephei star Nu Eridani. Half a year of observations have been obtained by four of the five satellites constituting BRITE-Constellation, supplemented with ground-based photoelectric photometry. We show that carefully combining the two data sets virtually eliminates the aliasing problem that often hampers time-series analyses. We detect 40 periodic signals intrinsic to the star in the light curves. Despite a lower detection limit we do not recover all the pressure and mixed modes previously reported in the literature, but we newly detect six additional gravity modes. This behaviour is a consequence of temporal changes in the pulsation amplitudes that we also detected for some of the p modes. We point out that the dependence of theoretically predicted pulsation amplitude on wavelength is steeper in visual passbands than those observationally measured, to the extent that the three dominant pulsation modes of Nu Eridani would be incorrectly identified using data in optical filters only. We discuss possible reasons for this discrepancy.

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.

BRITE (BRIght Target Explorer) Constellation, the first nanosatellite mission applied to astrophysical research, is a collaboration among Austria, Canada and Poland. The fleet of satellites (6 launched, 5 functioning) performs precise optical photometry of the brightest stars in the night sky. A pioneering mission like BRITE - with optics and instruments restricted to small volume, mass and power in several nanosatellites, whose measurements must be coordinated in orbit - poses many unique challenges. We discuss the technical issues, including problems encountered during on-orbit commissioning (especially higher-than expected sensitivity of the CCDs to particle radiation). We describe in detail how the BRITE team has mitigated these problems, and provide a complete overview of mission operations. This paper serves as a template for how to effectively plan, build and operate future low-cost niche-driven space astronomy missions.

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.

This paper aims to precisely determine the masses and detect pulsation modes in the two massive components of Beta Cen with BRITE-Constellation photometry. In addition, seismic models for the components are considered and the effects of fast rotation are discussed. This is done to test the limitations of seismic modeling for this very difficult case. A simultaneous fit of visual and spectroscopic orbits is used to self-consistently derive the orbital parameters, and subsequently the masses, of the components. The derived masses are equal to 12.02 +/- 0.13 and 10.58 +/- 0.18 M_Sun. The parameters of the wider, A - B system, presently approaching periastron passage, are constrained. Analysis of the combined blue- and red-filter BRITE-Constellation photometric data of the system revealed the presence of 19 periodic terms, of which eight are likely g modes, nine are p modes, and the remaining two are combination terms. It cannot be excluded that one or two low-frequency terms are rotational frequencies. It is possible that both components of Beta Cen are Beta Cep/SPB hybrids. An attempt to use the apparent changes of frequency to distinguish which modes originate in which component did not succeed, but there is potential for using this method when more BRITE data become available. Agena seems to be one of very few rapidly rotating massive objects with rich p- and g-mode spectra, and precisely known masses. It can therefore be used to gain a better understanding of the excitation of pulsations in relatively rapidly rotating stars and their seismic modeling. Finally, this case illustrates the potential of BRITE-Constellation data for the detection of rich-frequency spectra of small-amplitude modes in massive pulsating stars.

We report on an analysis of high-precision, multi-colour photometric observations of the rapidly-oscillating Ap (roAp) star $\alpha$ Cir. These observations were obtained with the BRITE-Constellation, which is a coordinated mission of five nanosatellites that collects continuous millimagnitude-precision photometry of dozens of bright stars for up to 180 days at a time in two colours (Johnson B and R). BRITE stands for BRight Target Explorer. The object $\alpha$ Cir is the brightest roAp star and an ideal target for such investigations, facilitating the determination of oscillation frequencies with high resolution. This star is bright enough for complementary interferometry and time-resolved spectroscopy. Four BRITE satellites observed $\alpha$ Cir for 146 d or 33 rotational cycles. Phasing the photometry according to the 4.4790 d rotational period reveals qualitatively different light variations in the two photometric bands. The phased red-band photometry is in good agreement with previously-published WIRE data, showing a light curve symmetric about phase 0.5 with a strong contribution from the first harmonic. The phased blue-band data, in contrast, show an essentially sinusoidal variation. We model both light curves with Bayesian Photometric Imaging, which suggests the presence of two large-scale, photometrically bright (relative to the surrounding photosphere) spots. We also examine the high-frequency pulsation spectrum as encoded in the BRITE photometry. Our analysis establishes the stability of the main pulsation frequency over the last 20 years, confirms the presence of frequency f7, which was not detected (or the mode not excited) prior to 2006, and excludes quadrupolar modes for the main pulsation frequency.

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.

We have used the MOST (Microvariability and Oscillations of STars) microsatellite to obtain four weeks of contiguous high-precision broadband visual photometry of the O7.5III(n)((f)) star {\xi} Persei in November 2011. This star is well known from previous work to show prominent DACs (Discrete Absorption Components) on time-scales of about 2 d from UV spectroscopy and NRP (Non Radial Pulsation) with one (l = 3) p-mode oscillation with a period of 3.5 h from optical spectroscopy. Our MOST-orbit (101.4 min) binned photometry fails to reveal any periodic light variations above the 0.1 mmag 3-sigma noise level for periods of hours, while several prominent Fourier peaks emerge at the 1 mmag level in the two-day period range. These longer-period variations are unlikely due to pulsations, including gravity modes. From our simulations based upon a simple spot model, we deduce that we are seeing the photometric modulation of several co-rotating bright spots on the stellar surface. In our model, the starting times (random) and lifetimes (up to several rotations) vary from one spot to another yet all spots rotate at the same period of 4.18 d, the best-estimated rotation period of the star. This is the first convincing reported case of co-rotating bright spots on an O star, with important implications for drivers of the DACs (resulting from CIRs - Corotating Interaction Regions) with possible bright-spot generation via a breakout at the surface of a global magnetic field generated by a subsurface convection zone.