PHOEBE 2 is a Python package for modeling the observables of eclipsing star systems, but until now has focused entirely on the forward-model -- that is, generating a synthetic model given fixed values of a large number of parameters describing the system and the observations. The inverse problem, obtaining orbital and stellar parameters given observational data, is more complicated and computationally expensive as it requires generating a large set of forward-models to determine which set of parameters and uncertainties best represent the available observational data. The process of determining the best solution and also of obtaining reliable and robust uncertainties on those parameters often requires the use of multiple algorithms, including both optimizers and samplers. Furthermore, the forward-model of PHOEBE has been designed to be as physically robust as possible, but is computationally expensive compared to other codes. It is useful, therefore, to use whichever code is most efficient given the reasonable assumptions for a specific system, but learning the intricacies of multiple codes presents a barrier to doing this in practice. Here we present the 2.3 release of PHOEBE (publicly available from http://phoebe-project.org) which introduces a general framework for defining and handling distributions on parameters, and utilizing multiple different estimation, optimization, and sampling algorithms. The presented framework supports multiple forward-models, including the robust model built into PHOEBE itself.

Results of an analysis of the BRITE-Constellation photometry of the SB1 system and ellipsoidal variable $\pi^5$ Ori (B2\,III) are presented. In addition to the orbital light-variation, which can be represented as a five-term Fourier cosine series with the frequencies $f_{\rm orb}$, $2f_{\rm orb}$, $3f_{\rm orb}$, $4f_{\rm orb}$ and $6f_{\rm orb}$, where $f_{\rm orb}$ is the system's orbital frequency, the star shows five low-amplitude but highly-significant sinusoidal variations with frequencies $f_i$ ($i ={}$2,..,5,7) in the range from 0.16 to 0.92~d$^{-1}$. With an accuracy better than 1$\sigma$, the latter frequencies obey the following relations: $f_2-f_4 = 2f_{\rm orb}$, $f_7 - f_3 = 2f_{\rm orb}$, $f_5 = f_3 - f_4 = f_7 - f_2$. We interpret the first two relations as evidence that two high-order $\ell = 1, m = 0$ gravity modes are self-excited in the system's tidally distorted primary component. The star is thus an ellipsoidal SPB variable. The last relations arise from the existence of the first-order differential combination term between the two modes. Fundamental parameters, derived from photometric data in the literature and the {\em Hipparcos\/} parallax, indicate that the primary component is close to the terminal stages of its main sequence (MS) evolution. Extensive Wilson-Devinney modeling leads to the conclusion that best fits of the theoretical to observed light-curves are obtained for the effective temperature and mass consistent with the primary's position in the HR diagram and suggests that the secondary is in an early MS evolutionary stage.

In this paper, we present results from the spectropolarimetric follow-up of photometrically selected candidate magnetic B stars from the MOBSTER project. Out of four observed targets, one (HD 38170) is found to host a detectable surface magnetic field, with a maximum longitudinal field measurement of 105$\pm$14 G. This star is chemically peculiar and classified as an $\alpha^2$ CVn variable. Its detection validates the use of TESS to perform a photometric selection of magnetic candidates. Furthermore, upper limits are derived for the remaining three stars, and we report the discovery of a previously unknown spectroscopic binary system, HD 25709. Finally, we use our non-detections as case studies to further inform the criteria to be used for the selection of a larger sample of stars to be followed up using high-resolution spectropolarimetry.

Classical novae are thermonuclear explosions that occur on the surfaces of white dwarf stars in interacting binary systems (Bode & Evans 2008). It has long been thought that the luminosity of classical novae is powered by continued nuclear burning on the surface of the white dwarf after the initial runaway (Gallaher & Starrfield 1978). However, recent observations of GeV $\gamma$-rays from classical novae have hinted that shocks internal to the nova ejecta may dominate the nova emission. Shocks have also been suggested to power the luminosity of events as diverse as stellar mergers (Metzger & Pejcha 2017), supernovae (Moriya et al. 2018), and tidal disruption events (Roth et al. 2016), but observational confirmation has been lacking. Here we report simultaneous space-based optical and $\gamma$-ray observations of the 2018 nova V906 Carinae (ASASSN-18fv), revealing a remarkable series of distinct correlated flares in both bands. The optical and $\gamma$-ray flares occur simultaneously, implying a common origin in shocks. During the flares, the nova luminosity doubles, implying that the bulk of the luminosity is shock-powered. Furthermore, we detect concurrent but weak X-ray emission from deeply embedded shocks, confirming that the shock power does not appear in the X-ray band and supporting its emergence at longer wavelengths. Our data, spanning the spectrum from radio to $\gamma$-ray, provide direct evidence that shocks can power substantial luminosity in classical novae and other optical transients.

Traditionally, the effects of interstellar extinction on binary star light curves have been treated as a uniform reduction in the observed brightness of the system that is independent of orbital phase. However, unless the orbital plane of the system coincides with the plane of the sky, or if the two stars are completely identical and present with minimal mutual irradiation and tidal/rotational distortions, then this is unlikely to be an accurate representation of the effect of interstellar extinction. Here, we present an updated treatment of interstellar extinction as incorporated in the PHOEBE 2.2 release (publicly available from http://phoebe-project.org) and assess the importance of using such an approach in the modeling of different types of binary systems. We also present the incorporation of PHOENIX model atmospheres into the PHOEBE 2.2 release, providing increased fidelity on computed observables down to lower temperatures than previously available. The importance of these new code developments is then highlighted via an extincted toy model of the eclipsing white-dwarf-subdwarf binary SDSS~J235524.29+044855.7 -- demonstrating that, in the age of LSST as well as complementary space-based photometric missions, a proper accounting for extinction and as well as the use of realistic model atmospheres will be essential in deriving accurate binary parameters.

We present an updated O-C diagram of the light-time variations of the eclipsing binary (component B) in the system QZ Carinae as it moves in the long-period orbit around the non-eclipsing pair (component A). This includes new Variable Stars South members' measures from 2017 to 2019, BRITE satellite observations in 2017 and 2018, and 100 previously unpublished measures made at Auckland Observatory from 1974 to 1978. We conclude that QZ Carinae has not yet completed one orbit of the two pairs since discovery in 1971. The duration of totality of primary eclipses was measured to be 0.295 +/- 0.02 day (7.08 +/-0.48 hours), rather longer than earlier values from light curve models. Other observational findings include the shape of primary and secondary eclipses and small-scale short-term brightness changes.

Being relatively rare, the properties of magnetic O stars are not fully understood. To date fewer than a dozen of these stars have been confirmed, making any inference of their global properties uncertain due to small number statistics. To better understand these objects it is necessary to increase the known sample. The MOBSTER collaboration aims to do this by identifying candidate magnetic O, B, and A stars from the identification of rotational modulation in high-precision photometry from the Transiting Exoplanet Survey Satellite (TESS). Here we discuss the collaboration's efforts to detect rotational modulation in TESS targets to identify candidate magnetic O stars for future spectropolarimetric observations.

In this contribution, we present the MOBSTER Collaboration, a large community effort to leverage high-precision photometry from the Transiting Exoplanet Survey Satellite (\textit{TESS}) in order to characterize the variability of magnetic massive and intermediate-mass stars. These data can be used to probe the varying column density of magnetospheric plasma along the line of sight for OB stars, thus improving our understanding of the interaction between surface magnetic fields and massive star winds. They can also be used to map out the brightness inhomogeneities present on the surfaces of Ap/Bp stars, informing present models of atomic diffusion in their atmospheres. Finally, we review our current and ongoing studies, which lead to new insights on this topic.

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.

$\varepsilon$ Lupi A is a binary system consisting of two main sequence early B-type stars Aa and Ab in a short period, moderately eccentric orbit. The close binary pair is the only doubly-magnetic massive binary currently known. Using photometric data from the BRITE-Constellation we identify a modest heartbeat variation. Combining the photometry with radial velocities of both components we determine a full orbital solution including empirical masses and radii. These results are compared with stellar evolution models as well as interferometry and the differences discussed. We also find additional photometric variability at several frequencies, finding it unlikely these frequencies can be caused by tidally excited oscillations. We do, however, determine that these signals are consistent with gravity mode pulsations typical for slowly pulsating B stars. Finally we discuss how the evolution of this system will be affected by magnetism, determining that tidal interactions will still be dominant.

Observations of Beta Lyr in four months of 2018 by three BRITE Constellation satellites (the red-filter BTr and BHr, and the blue-filter BLb) permitted a first, limited look into the light-curve variability in two spectral bands. The variations were found to be well correlated outside the innermost primary minima with the blue variations appearing to have smaller amplitudes than the red; this reduction may reflect their presumed origin in the cooler, outer parts of the accretion disk. This result must be confirmed with more extensive material as the current conclusions are based on observations spanning slightly less than three orbital cycles of the binary. The assumption of an instrumental problem and the applied corrections made to explain the unexpectedly large amplitude of the red-filter light-curve observed with the BTr satellite in 2016 are fully confirmed by the 2018 results.

In this paper we introduce the MOBSTER collaboration and lay out its scientific goals. We present first results based on the analysis of nineteen previously known magnetic O, B and A stars observed in 2-minute cadence in sectors 1 and 2 of the Transiting Exoplanet Survey Satellite (TESS) mission. We derive precise rotational periods from the newly obtained light curves and compare them to previously published values. We also discuss the overall photometric phenomenology of the known magnetic massive and intermediate-mass stars and propose an observational strategy to augment this population by taking advantage of the high-quality observations produced by TESS.

We analyse time-series observations from the BRITE-Constellation of the well known $\beta$ Cephei type star $\theta$ Ophiuchi. Seven previously known frequencies were confirmed and nineteen new frequency peaks were detected. In particular, high-order g modes, typical for the SPB (Slowly Pulsating B-type star) pulsators, are uncovered. These low-frequency modes are also obtained from the 7-year SMEI light curve. If g modes are associated with the primary component of $\theta$ Oph, then our discovery allows, as in the case of other hybrid pulsators, to infer more comprehensive information on the internal structure. To this aim we perform in-depth seismic studies involving simultaneous fitting of mode frequencies, reproducing mode instability and adjusting the relative amplitude of the bolometric flux variations. To explain the mode instability in the observed frequency range a significant increase of the mean opacity in the vicinity of the $Z$-bump is needed. Moreover, constraints on mass, overshooting from the convective core and rotation are derived. If the low-frequency modes come from the speckle B5 companion then taking into account the effects of rotation is enough to explain the pulsational mode instability.

Binary systems where the axis of rotation (spin) of one or both components is tilted w.r.t. the axis of revolution are called misaligned binary systems. The angle of misalignment, obliquity, has been measured for a handful of stars and extrasolar planets to date. Here we present a mathematical framework for a complete and rigorous treatment of misalignment and introduce an extension to the public PHOEBE code that implements this framework. We discuss misalignment for the Roche geometry and introduce methods for computing stellar shapes, equilibrium (generalized Lagrange) points of the potential and minimal requirements for lobe existence. Efficient parametrization of misalignment is proposed in the plane-of-sky coordinates and implementation details in PHOEBE are given alongside the proof-of-concept toy model, comparison with a known misaligned binary DI Her, and comparison with a misaligned planetary system Kepler-13. We provide important mathematical details of the model in the Appendix. This paper accompanies the release of PHOEBE 2.1, which will be available soon from its website http://phoebe-project.org.

HD149404 is an evolved non-eclipsing O-star binary that has previously undergone a Roche lobe overflow interaction. Understanding some key properties of the system requires a determination of the orbital inclination and of the dimensions of the components. The BRITE-Heweliusz satellite was used to collect photometric data of HD149404. Additional photometry was retrieved from the SMEI archive. These data were analysed using a suite of period search tools. The orbital part of the lightcurve was modelled with the nightfall binary star code. The Gaia-DR2 parallax of HD149404 was used to provide additional constraints. The periodograms reveal a clear orbital modulation of the lightcurve with a peak-to-peak amplitude near 0.04 mag. The remaining non-orbital part of the variability is consistent with red noise. The lightcurve folded with the orbital period reveals ellipsoidal variations, but no eclipses. The minimum when the secondary star is in inferior conjunction is deeper than the other minimum due to mutual reflection effects between the stars. Combined with the Gaia-DR2 parallaxes, the photometric data indicate an orbital inclination in the range of 23{\deg} to 31{\deg} and a Roche lobe filling factor of the secondary larger than or equal to 0.96. The luminosity of the primary star is consistent with its present-day mass, whereas the more evolved secondary appears overluminous for its mass. We confirm that the primary's rotation period is about half the orbital period. Both features most probably stem from the past Roche lobe overflow episode.

In early-type Be stars, groups of nonradial pulsation (NRP) modes with numerically related frequencies may be instrumental for the release of excess angular momentum through mass-ejection events. Difference and sum/harmonic frequencies often form additional groups. The goal of this study is to find out whether a similar frequency pattern occurs in the cooler third-magnitude B7-8\,IIIe shell star $\nu$ Pup. Time-series analyses are performed of space photometry with BRITE-Constellation (2015, 2016/17, and 2017/18), SMEI (2003--011), and Hipparcos (1989-1993). Two IUE SWP and 27 optical echelle spectra spanning 20 years were retrieved from various archives. The optical spectra exhibit no anomalies or well-defined variabilities. A magnetic field was not detected. All three photometry satellites recorded variability near 0.656 c/d which is resolved into three features separated by ~0.0021 c/d. First harmonics form a second frequency group, also spaced by ~0.0021 c/d. The frequency spacing is very nearly but not exactly equidistant. Variability near 0.0021 c/d was not detected. The long-term frequency stability could be used to derive meaningful constraints on the properties of a putative companion star. The IUE spectra do not reveal the presence of a hot subluminous secondary. $\nu$\,Pup is another Be star exhibiting an NRP variability pattern with long-term constancy and underlining the importance of combination frequencies and frequency groups. The star is a good target for efforts to identify an effectively single Be star.

Photometric instabilities of $\beta$ Lyr were observed in 2016 by two red-filter BRITE satellites over more than 10 revolutions of the binary, with $\sim$100-minute sampling. Analysis of the time series shows that flares or fading events take place typically 3 to 5 times per binary orbit. The amplitudes of the disturbances (relative to the mean light curve, in units of the maximum out-of-eclipse light-flux, f.u.) are characterized by a Gaussian distribution with $\sigma=0.0130\pm0.0004$ f.u. Most of the disturbances appear to be random, with a tendency to remain for one or a few orbital revolutions, sometimes changing from brightening to fading or the reverse. Phases just preceding the center of the deeper eclipse showed the most scatter while phases around secondary eclipse were the quietest. This implies that the invisible companion is the most likely source of the instabilities. Wavelet transform analysis showed domination of the variability scales at phase intervals $0.05-0.3$ (0.65--4 d), with the shorter (longer) scales dominating in numbers (variability power) in this range. The series can be well described as a stochastic Gaussian process with the signal at short timescales showing a slightly stronger correlation than red noise. The signal de-correlation timescale $\tau=(0.068\pm0.018)$ in phase or $(0.88\pm0.23)$~d appears to follow the same dependence on the accretor mass as that observed for AGN and QSO masses 5--9 orders of magnitude larger than the $\beta$~Lyr torus-hidden component.

SMC AB 6 is the shortest-period (6.5d) Wolf-Rayet (WR) binary in the Small Magellanic Cloud, and is therefore crucial for the study of binary interaction and formation of WR stars at low metallicity. The WR component in AB 6 was previously found to be very luminous (logL=6.3[Lsun]) compared to its reported orbital mass (8Msun), placing it significantly above the Eddington limit. Through spectroscopy and orbital analysis of newly acquired optical data taken with UVES, we aim to understand the peculiar results reported for this system and explore its evolutionary history. Results: We find that AB 6 contains at least four stars. The 6.5d period WR binary comprises the WR primary (WN3:h, star A) and a rather rapidly rotating early O-type companion (O5.5 V, star B). Static N and He lines suggest the presence of an emission line star (O5.5 I(f), star C). Finally, narrow absorption lines portraying a long-term radial velocity variation show the existence of a fourth star (O7.5 V, star D). Star D appears to form a second 140d period binary together with a fifth stellar member, which is a B-type dwarf or a black hole. It is not clear that these additional components are bound to the WR binary. The WR star is found to be less luminous than previously thought (logL = 5.9[Lsun]) and, adopting 41Msun for star B, more massive (18Msun). Correspondingly, the WR star does not exceed the Eddington limit. We derive the initial masses of 60 and 40Msun for stars A and B and an age of 3.9 Myr for the system. The WR binary likely experienced nonconservative mass transfer in the past supported by the relatively rapid rotation of star B. Conclusion: Our study shows that AB 6 is a multiple -- probably quintuple -- system. This finding resolves the previously reported puzzle of the WR primary exceeding the Eddington limit and suggests that the WR star exchanged mass with its companion in the past.

$\eta$ Car is a massive, eccentric binary with a rich observational history. We obtained the first high-cadence, high-precision light curves with the BRITE-Constellation nanosatellites over 6 months in 2016 and 6 months in 2017. The light curve is contaminated by several sources including the Homunculus nebula and neighboring stars, including the eclipsing binary CPD$-$59$^\circ$2628. However, we found two coherent oscillations in the light curve. These may represent pulsations that are not yet understood but we postulate that they are related to tidally excited oscillations of $\eta$ Car's primary star, and would be similar to those detected in lower-mass eccentric binaries. In particular, one frequency was previously detected by van Genderen et al. and Sterken et al. through the time period of 1974 to 1995 through timing measurements of photometric maxima. Thus, this frequency seems to have been detected for nearly four decades, indicating that it has been stable in frequency over this time span. These pulsations could help provide the first direct constraints on the fundamental parameters of the primary star if confirmed and refined with future observations.

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.

Observations of 25 Ori much expand the picture derived of other early-type Be stars with BRITE and SMEI. Two instead of one difference frequencies rule the variability: (a) The lower one, 0.0129 c/d, is the frequency of events with full amplitudes of 100-200 mmag which may signal mass loss possibly driven by the higher one, 0.1777 c/d. (b) Much of the entire power spectrum is a tightly woven network of combination frequencies: (i) Below 0.25\,c/d, numerous frequencies are difference frequencies. (ii) Many frequencies above 2.5 c/d can be represented as sum frequencies and in a few cases as harmonics. (iii) Many frequencies between 1.1 and 1.75\,c/d can be portrayed as parents of combination frequencies. The number and fraction of combination frequencies increases steeply with decreasing amplitude and and accuracy of the frequency matching.

The BRITE Constellation of nanosatellites obtained mmag photometry of 28 Cygni for 11 months in 2014-2016. Observations with the Solar Mass Ejection Imager in 2003-2010 and 118 H$\alpha$ line profiles were added. For decades, 28 Cyg has exhibited four large-amplitude frequencies: two closely spaced frequencies of spectroscopically confirmed $g$ modes near 1.5 c/d, one slightly lower exophotospheric (Stefl) frequency, and at 0.05 c/d the difference frequency between the two g modes. This top-level framework is indistinguishable from eta Cen (Paper I), which is also very similar in spectral type, rotation rate, and viewing angle. The Stefl frequency is the only one that does not seem to be affected by the difference frequency. The amplitude of the latter undergoes large variations; around maximum the amount of near-circumstellar matter is increased, and the amplitude of the Stefl frequency grows by some factor. During such brightenings dozens of transient spikes appear in the frequency spectrum, concentrated in three groups. Only eleven frequencies were common to all years of BRITE observations. Be stars seem to be controlled by several coupled clocks, most of which are not very regular on timescales of weeks to months but function for decades. The combination of g modes to the low difference frequency and/or the atmospheric response to it appears significantly nonlinear. Like in eta Cen, the difference-frequency variability seems the main responsible for the modulation of the star-to-disc mass transfer in 28 Cyg. A hierarchical set of difference frequencies may reach the longest timescales known of the Be phenomenon.

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.

Stellar rotation affects the transport of chemical elements and angular momentum and is therefore a key process during stellar evolution, which is still not fully understood. This is especially true for massive stars, which are important for the chemical enrichment of the universe. It is therefore important to constrain their physical parameters and internal angular momentum distribution to calibrate stellar structure and evolution models. Stellar internal rotation can be probed through asteroseismic studies of rotationally split oscillations but such results are still quite rare, especially for stars more massive than the Sun. The SPB star HD201433 is known to be part of a single-lined spectroscopic triple system, with two low-mass companions orbiting with periods of about 3.3 and 154 d. Our results are based on photometric observations made by BRITE - Constellation and the SMEI on board the Coriolis satellite, high-resolution spectroscopy, and more than 96 years of radial velocity measurements. We identify a sequence of 9 rotationally split dipole modes in the photometric time series and establish that HD201433 is in principle a solid-body rotator with a very slow rotation period of 297+/-76 d. Tidal interaction with the inner companion has, however, significantly accelerated the spin of the surface layers by a factor of approximately one hundred. The angular momentum transfer onto the surface of HD201433 is also reflected by the statistically significant decrease of the orbital period of about 0.9 s during the last 96 years. Combining the asteroseismic inferences with the spectroscopic measurements and the orbital analysis of the inner binary system, we conclude that tidal interactions between the central SPB star and its inner companion have almost circularised the orbit but not yet aligned all spins of the system and have just begun to synchronise rotation.

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}$ .

The precision of photometric and spectroscopic observations has been systematically improved in the last decade, mostly thanks to space-borne photometric missions and ground-based spectrographs dedicated to finding exoplanets. The field of eclipsing binary stars strongly benefited from this development. Eclipsing binaries serve as critical tools for determining fundamental stellar properties (masses, radii, temperatures and luminosities), yet the models are not capable of reproducing observed data well either because of the missing physics or because of insufficient precision. This led to a predicament where radiative and dynamical effects, insofar buried in noise, started showing up routinely in the data, but were not accounted for in the models. PHOEBE (PHysics Of Eclipsing BinariEs; http://phoebe-project.org) is an open source modeling code for computing theoretical light and radial velocity curves that addresses both problems by incorporating missing physics and by increasing the computational fidelity. In particular, we discuss triangulation as a superior surface discretization algorithm, meshing of rotating single stars, light time travel effect, advanced phase computation, volume conservation in eccentric orbits, and improved computation of local intensity across the stellar surfaces that includes photon-weighted mode, enhanced limb darkening treatment, better reflection treatment and Doppler boosting. Here we present the concepts on which PHOEBE is built on and proofs of concept that demonstrate the increased model fidelity.

Magnetic B-type stars exhibit photometric variability due to diverse causes, and consequently on a variety of timescales. In this paper we describe interpretation of BRITE photometry and related ground-based observations of 4 magnetic B-type systems: $\epsilon$ Lupi, $\tau$ Sco, a Cen and $\epsilon$ CMa.

New high-cadence observations with BRITE covering many months confirm that coupled pairs of nonradial pulsation modes are widespread among early-type Be stars. With the difference frequency between the parental variations they may form a roughly sinusoidal variability or the amplitude may cyclicly vary. A first - amplified - beat pattern is also found. In all three cases the amplitudes of difference frequencies can exceed the amplitude sum of the base frequencies, and modulations of the star-to-circumstellar-disk mass-transfer rate may be associated with these slow variations. This suggests more strongly than any earlier observations that significant dissipation of pulsational energy in the atmosphere may be a cause of mass ejections from Be stars. A unifying interpretative concept is presented.

Extracted BRITE lightcurves must be carefully prepared and corrected for instrumental effects before a scientific analysis can be performed. Therefore, we have created a suite of Python routines to prepare and correct the lightcurves, which is publicly available. In this paper we describe the method and successive steps performed by these routines.

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.

Massive WR stars are evolved massive stars characterized by strong mass-loss. Hypothetically, they can form either as single stars or as mass donors in close binaries. About 40% of the known WR stars are confirmed binaries, raising the question as to the impact of binarity on the WR population. By performing a spectral analysis of all multiple WR systems in the SMC, we obtain the full set of stellar parameters for each individual component. Mass-luminosity relations are tested, and the importance of the binary evolution channel is assessed. The spectral analysis is performed with the PoWR model atmosphere code by superimposing model spectra that correspond to each component. Evolutionary channels are constrained using the BPASS evolution tool. Significant Hydrogen mass fractions (0.1 - 0.4) are detected in all WN components. A comparison with mass-luminosity relations and evolutionary tracks implies that the majority of the WR stars in our sample are not chemically homogeneous. The WR component in the binary AB 6 is found to be very luminous (Log L ~ 6.3 [Lsun]) given its orbital mass (~10 Msun), presumably because of observational contamination by a third component. Evolutionary paths derived for our objects suggest that Roche lobe overflow had occurred in most systems, affecting their evolution. However, the implied initial masses are large enough for the primaries to have entered the WR phase, regardless of binary interaction. Together with the results for the putatively single SMC WR stars, our study suggests that the binary evolution channel does not dominate the formation of WR stars at SMC metallicity.

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.

Empirical evidence for the involvement of nonradial pulsations (NRP's) in the mass loss from Be stars ranges from (i) a singular case (\object{$\mu$ Cen}) of repetitive mass ejections triggered by multi-mode beating to (ii) several photometric reports about enormous numbers of pulsation modes popping up during outbursts and on to (iii) effective single-mode pulsators. The BRITE Constellation of nanosatellites was used to obtain mmag photometry of the Be stars $\eta$ and \object{$\mu$ Cen}. In the low-inclination star \object{$\mu$ Cen}, light pollution by variable amounts of near-stellar matter prevented any new insights into the variability and other properties of the central star. In the equator-on star \object{$\eta$ Cen}, BRITE photometry and {\sc Heros} echelle spectroscopy from the 1990s reveal an intricate clockwork of star-disk interactions. The mass transfer is modulated with the frequency difference of two NRP modes and an amplitude three times as large as the amplitude sum of the two NRP modes. This process feeds a high-amplitude circumstellar activity running with the incoherent and slightly lower so-called \v{S}tefl frequency. The mass loss-modulation cycles are tightly coupled to variations in the value of the \v{S}tefl frequency and in its amplitude, albeit with strongly drifting phase differences. The observations are well described by the decomposition of the mass loss into a pulsation-related engine in the star and a viscosity-dominated engine in the circumstellar disk. Arguments are developed that large-scale gas-circulation flows occur at the interface. The propagation rates of these eddies manifest themselves as \v{S}tefl frequencies. Bursts in power spectra during mass-loss events can be understood as the noise inherent to these gas flows.

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 present an overview of four phase-constrained Chandra HETGS X-ray observations of Delta Ori A. Delta Ori A is actually a triple system which includes the nearest massive eclipsing spectroscopic binary, Delta Ori Aa, the only such object which 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 A 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 ksec and covering nearly the entire binary orbit. Companion papers discuss the X-ray variability seen in the Chandra spectra, present new space-based photometry and ground-based radial velocities simultaneous with the X-ray data to better constrain the system parameters, and model the effects of X-rays on the optical and UV spectrum. 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 \ion{Fe}{17} and \ion{Ne}{10} are inconsistent with model predictions, which may be an effect of resonance scattering

We present time-resolved and phase-resolved variability studies of an extensive X-ray high-resolution spectral dataset of the $\delta$ Orionis Aa binary system. The four observations, obtained with Chandra ACIS HETGS, have a total exposure time of ~479 ks and provide nearly complete binary phase coverage. Variability of the total X-ray flux in the range 5-25 $\AA$ is confirmed, with maximum amplitude of about +/-15% within a single ~125 ks observation. Periods of 4.76d and 2.04d are found in the total X-ray flux, as well as an apparent overall increase in flux level throughout the 9-day observational campaign. Using 40 ks contiguous spectra derived from the original observations, we investigate 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 phase=0.0 when the secondary $\delta$ Orionis Aa2 is at 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.

We report on both high-precision photometry from the 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.7d, and a more distant tertiary (O9 IV P > 400 yrs). This data was collected in concert with X-ray spectroscopy from the Chandra X-ray Observatory. Thanks to continuous coverage for 3 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 curve of Aa1. While we are unable to recover radial velocity 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 which show evidence for tidally induced pulsations.

Eclipsing systems of massive stars allow one to explore the properties of their components in great detail. We perform a multi-wavelength, non-LTE analysis of the three components of the massive multiple system $\delta$ Ori A, focusing on the fundamental stellar properties, stellar winds, and X-ray characteristics of the system. The primary's distance-independent parameters turn out to be characteristic for its spectral type (O9.5 II), but usage of the ${\rm \it Hipparcos}$ parallax yields surprisingly low values for the mass, radius, and luminosity. Consistent values follow only if $\delta$ Ori lies at about twice the ${\rm \it Hipparcos}$ distance, in the vicinity of the $\sigma$-Orionis cluster. The primary and tertiary dominate the spectrum and leave the secondary only marginally detectable. We estimate the V-band magnitude difference between primary and secondary to be $\Delta V \approx 2.\!\!^{\rm m}8$. The inferred parameters suggest the secondary is an early B-type dwarf ($\approx$ B1 V), while the tertiary is an early B-type subgiant ($\approx$ B0 IV). We find evidence for rapid turbulent velocities ($\sim 200$ km ${\rm s}^{-1}$) and wind inhomogeneities, partially optically thick, in the primary's wind. The bulk of the X-ray emission likely emerges from the primary's stellar wind ($\log L_{\text{X}} / L_{\text{Bol}} \approx -6.85$), initiating close to the stellar surface at $R_0 \sim 1.1\,R_*$. Accounting for clumping, the mass-loss rate of the primary is found to be $\log \dot{M} \approx -6.4\,[M_\odot\,{\rm yr}^{-1}]$, which agrees with hydrodynamic predictions, and provides a consistent picture along the X-ray, UV, optical and radio spectral domains.

We report on extended photometry of two pulsating sdB stars in close binaries. For both cases, we use rotational splitting of the pulsation frequencies to show that the sdB component rotates much too slowly to be in synchronous rotation. We use a theory of tidal interaction in binary stars to place limits on the mass ratios that are independent of estimates based on the radial velocity curves. The companions have masses below 0.26 M\odot. The pulsation spectra show the signature of high-overtone g-mode pulsation. One star, KIC 11179657, has a clear sequence of g-modes with equal period spacings as well as several periodicities that depart from that trend. KIC 02991403 shows a similar sequence, but has many more modes that do not fit the simple pattern.

We report on Kepler photometry of the hot sdB star B4 in the open cluster NGC 6791. We confirm that B4 is a reflection effect binary with an sdB component and a low-mass main sequence companion with a circular 0.3985 d orbit. The sdB star is a g-mode pulsator (a V1093 Her star) with periods ranging from 2384 s to 7643 s. Several of the pulsation modes show symmetric splitting by 0.62 microHz. Attributing this to rotational splitting, we conclude that the sdB component has a rotation period of approximately 9.63 d, indicating that tidal synchronization has not been achieved in this system. Comparison with theoretical synchronization time provides a discriminant between various theoretical models.

KPD 1930+2752 is a short-period pulsating subdwarf B (sdB) star. It is also an ellipsoidal variable with a known binary period just over two hours. The companion is most likely a white dwarf and the total mass of the system is close to the Chandresakhar limit. In this paper we report the results of Whole Earth Telescope (WET) photometric observations during 2003 and a smaller multisite campaign from 2002. From 355 hours of WET data, we detect 68 pulsation frequencies and suggest an additional 13 frequencies within a crowded and complex temporal spectrum between 3065 and 6343 $\mu$Hz (periods between 326 and 157 s). We examine pulsation properties including phase and amplitude stability in an attempt to understand the nature of the pulsation mechanism. We examine a stochastic mechanism by comparing amplitude variations with simulated stochastic data. We also use the binary nature of KPD 1930+2752 for identifying pulsation modes via multiplet structure and a tidally-induced pulsation geometry. Our results indicate a complicated pulsation structure that includes short-period ($\approx 16$ h) amplitude variability, rotationally split modes, tidally-induced modes, and some pulsations which are geometrically limited on the sdB star.