We have obtained a very deep exposure (813 ks) of $\zeta\,$Puppis (O4 supergiant) with the Chandra/HETG Spectrometer. Here we report on analysis of the 1-9 \r{A} region, especially well suited for Chandra, which has a significant contribution from continuum emission between well separated emission lines from high-ionization species. These data allow us to study the hottest plasma present through the continuum shape and emission line strengths. Assuming a powerlaw emission measure distribution which has a high-temperature cut-off, we find that the emission is consistent with a thermal spectrum having a maximum temperature of 12 MK. This implies an effective wind shock velocity of $900\,\mathrm{km\,s^{-1}}$, well below the wind terminal speed of $2250\,\mathrm{km\,s^{-1}}$. For X-ray emission which forms close to the star, the speed and X-ray flux are larger than can be easily reconciled with strictly self-excited line-deshadowing-instability models, suggesting a need for a fraction of the wind to be accelerated extremely rapidly right from the base. This is not so much a dynamical instability as a nonlinear response to changing boundary conditions.

In circumstellar gas, the complex organic molecule methanol has been found almost exclusively around young stellar objects, and is thus regarded as a signpost of recent star formation. Here we report the first probable detection of methanol around an evolved high-mass star, in the complex circumstellar environment around the Luminous Blue Variable $\eta$ Carinae, while using ALMA to investigate molecular cloud conditions traced by CO (2-1) in an orbit phase of the massive binary preceding the 2020 periastron. Favoring methanol over a $^{13}$CS alternative, the emission originates from hot ($T_{\rm{gas}} \simeq$ 700 K) material, $\sim$2$''$ (0.02 pc) across, centered on the dust-obscured binary in contrast to the CO which traces inner layers of the extended massive equatorial torus, and is accompanied by prominent absorption in a cooler ($T_{\rm{gas}} \simeq$ 110 K) layer of gas. We also report detections of water in $Herschel$/HIFI observations at 557 GHz and 988 GHz. The methanol abundance is several to 50 times higher than observed towards several lower mass stars, while water abundances are similar to those observed in cool, dense molecular clouds. The very high methanol:water abundance ratio in the core of $\eta$ Carinae may suggest methanol formation processes similar to Fischer-Tropsch-type catalytic reactions on dust grains. These observations prove that complex molecule formation can occur in the chemically evolved environments around massive stars in the end stages of their evolution, given sufficient gas densities and shielding conditions as may occur in material around massive interacting companions and merger remnants.

Context. EE Cep is one of few eclipsing binary systems with a dark, dusty disk around an invisible object similar to {\epsilon} Aur. The system is characterized by grey and asymmetric eclipses every 5.6 yr, with a significant variation in their photometric depth, ranging from ~ 0 m .5 to ~ 2 m .0. Aims. The main aim of the observational campaign of the EE Cep eclipse in 2014 was to test the model of disk precession (Galan et al. 2012). We expected that this eclipse would be one of the deepest with a depth of ~ 2 m .0. Methods. We collected multicolor observations from almost 30 instruments located in Europe and North America. This photometric data covers 243 nights during and around the eclipse. We also analyse the low- and high-resolution spectra from several instruments. Results. The eclipse was shallow with a depth of 0 m .71 in V-band. The multicolor photometry illustrates small color changes during the eclipse with a total amplitude of order ~ +0 m . 15 in B-I color index. The linear ephemeris for this system is updated by including new times of minima, measured from the three most recent eclipses at epochs E = 9, 10 and 11. New spectroscopic observations were acquired, covering orbital phases around the eclipse, which were not observed in the past and increased the data sample, filling some gaps and giving a better insight into the evolution of the H {\alpha} and NaI spectral line profiles during the primary eclipse. Conclusions. The eclipse of EE Cep in 2014 was shallower than expected 0 m .71 instead of ~ 2 m . 0. This means that our model of disk precession needs revision.

Rapid rotation is a fundamental characteristic of classical Be stars and a crucial property allowing for the formation of their circumstellar disks. Past evolution in a mass and angular momentum transferring binary system offers a plausible solution to how Be stars attained their fast rotation. Although the subdwarf remnants of mass donors in such systems should exist in abundance, only a few have been confirmed due to tight observational constraints. An indirect method of detecting otherwise hidden companions is offered by their effect on the outer parts of Be star disks, which are expected to be disrupted or truncated. In the context of the IR and radio continuum excess radiation originating in the disk, the disk truncation can be revealed by a turndown in the spectral energy distribution due to reduced radio flux levels. In this work we search for signs of spectral turndown in a sample of 57 classical Be stars with radio data, which include new data for 23 stars and the longest wavelength detections so far ($\lambda \approx$ 10\,cm) for 2 stars. We confidently detect the turndown for all 26 stars with sufficient data coverage (20 of which are not known to have close binary companions). For the remaining 31 stars, data are inconclusive as to whether the turndown is present or not. The analysis suggests that many if not all Be stars have close companions influencing their outer disks. If confirmed to be subdwarf companions, the mass transfer spin-up scenario might explain the existence of the vast majority of classical Be stars.

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

Cosmic-ray acceleration has been a long-standing mystery and despite more than a century of study, we still do not have a complete census of acceleration mechanisms. The collision of strong stellar winds in massive binary systems creates powerful shocks, which have been expected to produce high-energy cosmic-rays through Fermi acceleration at the shock interface. The accelerated particles should collide with stellar photons or ambient material, producing non-thermal emission observable in X-rays and gamma-rays. The supermassive binary star eta Carinae drives the strongest colliding wind shock in the solar neighborhood. Observations with non-focusing high-energy observatories indicate a high energy source near eta Carinae, but have been unable to conclusively identify eta Carinae as the source because of their relatively poor angular resolution. Here we present the first direct focussing observations of the non-thermal source in the extremely hard X-ray band, which is found to be spatially coincident with the star within several arc-seconds. These observations show that the source of non-thermal X-rays varies with the orbital phase of the binary, and that the photon index of the emission is similar to that derived through analysis of the gamma-ray spectrum. This is conclusive evidence that the high-energy emission indeed originates from non-thermal particles accelerated at colliding wind shocks.

We present a detailed study of the 2017 eruption of the classical nova ASASSN-17pf (LMCN 2017-11a), which is located in the Large Magellanic Cloud, including data from AAVSO, ASAS-SN, SALT, SMARTS, SOAR, and the Neil Gehrels \textit{Swift} Observatory. The optical light-curve is characterized by multiple maxima (flares) on top of a slowly evolving light-curve (with a decline time, $t_2>$ 100 d). The maxima correlate with the appearance of new absorption line systems in the optical spectra characterized by increasing radial velocities. We suggest that this is evidence of multiple episodes of mass-ejection with increasing expansion velocities. The line profiles in the optical spectra indicate very low expansion velocities (FWHM $\sim$ 190 km s$^{-1}$), making this nova one of the slowest expanding ever observed, consistent with the slowly evolving light-curve. The evolution of the colors and spectral energy distribution show evidence of decreasing temperatures and increasing effective radii for the pseudo-photosphere during each maximum. The optical and infrared light-curves are consistent with dust formation 125 days post-discovery. We speculate that novae showing several optical maxima have multiple mass-ejection episodes leading to shocks that may drive $\gamma$-ray emission and dust formation.

The interacting binary Eta Carinae remains one of the most enigmatic massive stars in our Galaxy despite over four centuries of observations. In this work, its light curve from the ultraviolet to the near-infrared is analysed using spatially resolved HST observations and intense monitoring at the La Plata Observatory, combined with previously published photometry. We have developed a method to separate the central stellar object in the ground-based images using HST photometry and applying it to the more numerous ground-based data, which supports the hypothesis that the central source is brightening faster than the almost-constant Homunculus. After detrending from long-term brightening, the light curve shows periodic orbital modulation ($\Delta V$ $\sim$ 0.6 mag) attributed to the wind-wind collision cavity as it sweeps around the primary star and it shows variable projected area to our line-of-sight. Two quasi-periodic components with time scales of 2-3 and 8-10 yr and low amplitude, $\Delta V$ $<$ 0.2 mag, are superimposed on the brightening light-curve, being the only stellar component of variability found, which indicates minimal stellar instability. Moreover, the light curve analysis shows no evidence of `shell ejections' at periastron. We propose that the long-term brightening of the stellar core is due to the dissipation of a dusty clump in front of the central star, which works like a natural coronagraph. Thus, the central stars appear to be more stable than previously thought since the dominant variability originates from a changing circumstellar medium. We predict that the brightening phase, due mainly to dust dissipation, will be completed around 2032 $\pm$ 4 yr, when the star will be brighter than in the 1600's by up to $\Delta V$ $\sim$ 1 mag.

We present a census of all known and candidate luminous blue variables in the Local Group.

We report the discovery of two massive eccentric systems with BRITE data, $\tau$ Ori and $\tau$ Lib, showing heartbeat effects close to the periastron passage. $\tau$ Lib exhibits shallow eclipses that will soon vanish due to the apsidal motion in the system. In neither system, tidally excited oscillations were detected.

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

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 MOST space telescope observed four Cepheid variables so far, all of different subtypes. Here we summarize the results obtained and the possible ways to continue to study these stars.

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.

Massive binaries (MBs) play a crucial role in the Universe. Knowing the distributions of their orbital parameters (OPs) is important for a wide range of topics, from stellar feedback to binary evolution channels, from the distribution of supernova types to gravitational wave progenitors, yet, no direct measurements exist outside the Milky Way. The Tarantula Massive Binary Monitoring was designed to help fill this gap by obtaining multi-epoch radial velocity monitoring of 102 MBs in the 30 Dor. In this paper, we analyse 32 VLT/FLAMES observations of 93 O- and 7 B-type binaries. We performed a Fourier analysis and obtained orbital solutions for 82 systems: 51 single- and 31 double-lined spectroscopic binaries. Overall, the OPs and binary fraction are remarkably similar across the 30 Dor region and compared to existing Galactic samples (GSs). This indicates that within these domains environmental effects are of second order in shaping the properties of MBs. A small difference is found in the distribution of orbital periods (OrbPs), which is slightly flatter (in log space) in 30 Dor than in the Galaxy, although this may be compatible within error estimates and differences in the fitting methodology. Also, OrbPs in 30 Dor can be as short as 1.1 d; somewhat shorter than seen in GSs. Equal mass binaries q>0.95 in 30 Dor are all found outside NGC 2070 the very young and massive cluster at 30 Dor's core. One outstanding exception however is the fact that earliest spectral types (O2-O7) tend to have shorter OrbPs than latter (O9.2-O9.7). Our results point to a relative universality of the incidence rate of MBs and their OPs in the metallicity range from solar ($Z_{\odot}$) to about $0.5Z_{\odot}$. This provides the first direct constraints on MB properties in massive star-forming galaxies at the Universes peak of star formation at redshifts z~1 to 2, which are estimated to have $Z~0.5Z_{\odot}$.

{\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 late-type Be star $\beta$ CMi is remarkably stable compared to other Be stars that have been studied. This has led to a realistic model of the outflowing Be disk by Klement et al. These results showed that the disk is likely truncated at a finite radius from the star, which Klement et al.~suggest is evidence for an unseen binary companion in orbit. Here we report on an analysis of the Ritter Observatory spectroscopic archive of $\beta$ CMi to search for evidence of the elusive companion. We detect periodic Doppler shifts in the wings of the H$\alpha$ line with a period of 170 d and an amplitude of 2.25 km s$^{-1}$, consistent with a low-mass binary companion ($M\approx 0.42 M_\odot$). We then compared the small changes in the violet-to-red peak height changes ($V/R$) with the orbital motion. We find weak evidence that it does follow the orbital motion, as suggested by recent Be binary models by Panoglou et al. Our results, which are similar to those for several other Be stars, suggest that $\beta$ CMi may be a product of binary evolution where Roche lobe overflow has spun up the current Be star, likely leaving a hot subdwarf or white dwarf in orbit around the star. Unfortunately, no direct sign of this companion star is found in the very limited archive of {\it International Ultraviolet Explorer} spectra.

Space-based photometric measurements first revealed low-amplitude irregularities in the pulsations of Cepheid stars, but their origins and how commonly they occur remain uncertain. To investigate this phenomenon, we present MOST space telescope photometry of two Cepheids. V473 Lyrae is a second-overtone, strongly modulated Cepheid, while U Trianguli Australis is a Cepheid pulsating simultaneously in the fundamental mode and first overtone. The nearly continuous, high-precision photometry reveals alternations in the amplitudes of cycles in V473 Lyr, the first case of period doubling detected in a classical Cepheid. In U TrA, we tentatively identify one peak as the $f_X$ or 0.61-type mode often seen in conjunction with the first radial overtone in Cepheids, but given the short length of the data, we cannot rule out that it is a combination peak instead. Ground-based photometry and spectroscopy were obtained to follow two modulation cycles in V473 Lyr and to better specify its physical parameters. The simultaneous data yield the phase lag parameter (the phase difference between maxima in luminosity and radial velocity) of a second-overtone Cepheid for the first time. We find no evidence for a period change in U TrA or an energy exchange between the fundamental mode and the first overtone during the last 50 years, contrary to earlier indications. Period doubling in V473 Lyr provides a strong argument that mode interactions do occur in some Cepheids and we may hypothesise that it could be behind the amplitude modulation, as recently proposed for Blazhko RR Lyrae stars.

We present the first SB2 orbital solution and disentanglement of the massive Wolf-Rayet binary R145 (P = 159d) located in the Large Magellanic Cloud. The primary was claimed to have a stellar mass greater than 300Msun, making it a candidate for the most massive star known. While the primary is a known late type, H-rich Wolf-Rayet star (WN6h), the secondary could not be so far unambiguously detected. Using moderate resolution spectra, we are able to derive accurate radial velocities for both components. By performing simultaneous orbital and polarimetric analyses, we derive the complete set of orbital parameters, including the inclination. The spectra are disentangled and spectroscopically analyzed, and an analysis of the wind-wind collision zone is conducted. The disentangled spectra and our models are consistent with a WN6h type for the primary, and suggest that the secondary is an O3.5 If*/WN7 type star. We derive a high eccentricity of e = 0.78 and minimum masses of M1 sin^3 i ~ M2 sin^3 i ~ 13 +- 2 Msun, with q = M2 / M1 = 1.01 +- 0.07. An analysis of emission excess stemming from a wind-wind collision yields a similar inclination to that obtained from polarimetry (i = 39 +- 6deg). Our analysis thus implies M1 = 53^{+40}_{-20} and M2 = 54^{+40}_{-20} Msun, excluding M1 > 300Msun. A detailed comparison with evolution tracks calculated for single and binary stars, as well as the high eccentricity, suggest that the components of the system underwent quasi-homogeneous evolution and avoided mass-transfer. This scenario would suggest current masses of ~ 80 Msun and initial masses of Mi,1 ~ 105 and Mi,2 ~ 90Msun, consistent with the upper limits of our derived orbital masses, and would imply an age of ~2.2 Myr.

Context. The mass loss from massive stars is not understood well. Eta Car is a unique object for studying the massive stellar wind during the LBV phase. It is also an eccentric binary with a period of 5.54 yr. The nature of both stars is uncertain, although we know from X-ray studies that there is a wind-wind collision whose properties change with orbital phase. Methods. Observations of Eta Car were carried out with the ESO VLTI and the AMBER instrument between approximately five and seven months before the August 2014 periastron passage. Velocity-resolved aperture-synthesis images were reconstructed from the spectrally dispersed interferograms. Interferometric studies can provide information on the binary orbit, the primary wind, and the wind collision. Results. We present velocity-resolved aperture-synthesis images reconstructed in more than 100 different spectral channels distributed across the Br Gamma 2.166 micrometer emission line. The intensity distribution of the images strongly depends on wavelength. At wavelengths corresponding to radial velocities of approximately -140 to -376 km/s measured relative to line center, the intensity distribution has a fan-shaped structure. At the velocity of -277 km/s, the position angle of the symmetry axis of the fan is ~ 126 degree. The fan-shaped structure extends approximately 8.0 mas (~ 18.8 au) to the southeast and 5.8 mas (~ 13.6 au) to the northwest, measured along the symmetry axis at the 16% intensity contour. The shape of the intensity distributions suggests that the obtained images are the first direct images of the innermost wind-wind collision zone. Therefore, the observations provide velocity-dependent image structures that can be used to test three-dimensional hydrodynamical, radiative transfer models of the massive interacting winds of Eta Car.

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.

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.

Supernova (SN) 2015bh (or SNhunt275) was discovered in NGC 2770 on 2015 February with an absolute magnitude of Mr ~ -13.4 mag, and was initially classified as a SN impostor. Here we present the photometric and spectroscopic evolution of SN 2015bh from discovery to late phases (~ 1 yr after). In addition, we inspect archival images of the host galaxy up to ~ 21 yr before discovery, finding a burst ~ 1 yr before discovery, and further signatures of stellar instability until late 2014. Later on, the luminosity of the transient slowly increases, and a broad light curve peak is reached after about three months. We propose that the transient discovered in early 2015 could be a core-collapse SN explosion. The pre-SN luminosity variability history, the long-lasting rise and faintness first light curve peak suggests that the progenitor was a very massive, unstable and blue star, which exploded as a faint SN because of severe fallback of material. Later on, the object experiences a sudden brightening of 3 mag, which results from the interaction of the SN ejecta with circumstellar material formed through repeated past mass-loss events. Spectroscopic signatures of interaction are however visible at all epochs. A similar chain of events was previously proposed for the similar interacting SN 2009ip.

Eta Carinae, the closest, active, massive binary containing a highly unstable LBV, exhibits expanding, compressed wind shells, seen in emission, that are spatially and spectrally resolved by \hst/\stis. Starting in June 2009, these structures were mapped across its 5.54-year, highly elliptical, binary orbit to follow temporal changes in the light of [Fe~III] 4659\AA\ and [Fe~II] 4815\AA. The emissions trace portions of fossil wind-shells, that were formed by wind-wind interactions across each cycle. Over the high-ionization state, dense arcs, photo-ionized by far ultraviolet radiation from the hot secondary, are seen in [Fe~III]. Other arcs, ionized by mid-ultraviolet radiation from the primary star, are seen in [Fe II]. The [Fe III] structures tend to be interior to [Fe II] structures that trace extensive, less disturbed primary wind. During the brief periastron passage when the secondary plunges deep into the primary's extremely dense wind, on the far side of primary star, high-ionization [Fe III] structures fade and reappear in [Fe II]. Multiple fossil wind-structures were traced across the 5.7-year monitoring interval. The strong similarity of the expanding [Fe II] shells suggests that the wind and photo-ionization properties of the massive binary have not changed substantially from one orbit to the next over the past several orbital cycles. These observations trace structures that can be used to test three-dimensional hydrodynamical and radiative-transfer models of massive, interacting winds. They also provide a baseline for following future changes in Eta Carinae, especially of its winds and photoionization properties.

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.

We have monitored the massive binary star $\eta$ Carinae with the CTIO/SMARTS 1.5~m telescope and CHIRON spectrograph from the previous apastron passage of the system through the recent 2014.6 periastron passage. Our monitoring has resulted in a large, homogeneous data set with an unprecedented time-sampling, spectral resolving power, and signal-to-noise. This allowed us to investigate temporal variability previously unexplored in the system and discover a kinematic structure in the P Cygni absorption troughs of neutral helium wind lines. The features observed occurred prior to the periastron passage and are seen as we look through the trailing arm of the wind-wind collision shock cone. We show that the bulk of the variability is repeatable across the last five periastron passages, and that the absorption occurs in the inner 230 AU of the system. In addition, we found an additional, high-velocity absorption component super-imposed on the P Cygni absorption troughs that has been previously un-observed in these lines, but which bears resemblance to the observations of the He~I $\lambda$10830 \AA\ feature across previous cycles. Through a comparison of the current smoothed particle hydrodynamical simulations, we show that the observed variations are likely caused by instabilities in the wind-wind collision region in our line-of-sight, coupled with stochastic variability related to clumping in the winds.

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.

The evolved, massive highly eccentric binary system, eta Carinae, underwent a periastron passage in the summer of 2014. We obtained two coordinated X-ray observations with XMM-Newton and NuSTAR during the elevated X-ray flux state and just before the X-ray minimum flux state around this passage. These NuSTAR observations clearly detected X-ray emission associated with eta Car extending up to ~50 keV for the first time. The NuSTAR spectrum above 10 keV can be fit with the bremsstrahlung tail from a kT ~6 keV plasma. This temperature is Delta kT ~2 keV higher than those measured from the iron K emission line complex, if the shocked gas is in collisional ionization equilibrium. This result may suggest that the companion star's pre-shock wind velocity is underestimated. The NuSTAR observation near the X-ray minimum state showed a gradual decline in the X-ray emission by 40% at energies above 5 keV in a day, the largest rate of change of the X-ray flux yet observed in individual eta Car observations. The column density to the hardest emission component, NH ~1e24 cm-2, marked one of the highest values ever observed for eta Car, strongly suggesting increased obscuration of the wind-wind colliding X-ray emission by the thick primary stellar wind prior to superior conjunction. Neither observation detected the power-law component in the extremely hard band that INTEGRAL and Suzaku observed prior to 2011. If the non-detection by NuSTAR is caused by absorption, the power-law source must be small and located very near the WWC apex. Alternatively, it may be that the power-law source is not related to either eta Car or the GeV gamma-ray source.

{\eta} Carinae is an extremely massive binary system in which rapid spectrum variations occur near periastron. Most notably, near periastron the He II $\lambda 4686$ line increases rapidly in strength, drops to a minimum value, then increases briefly before fading away. To understand this behavior, we conducted an intense spectroscopic monitoring of the He II $\lambda 4686$ emission line across the 2014.6 periastron passage using ground- and space-based telescopes. Comparison with previous data confirmed the overall repeatability of EW(He II $\lambda 4686$), the line radial velocities, and the timing of the minimum, though the strongest peak was systematically larger in 2014 than in 2009 by 26%. The EW(He II $\lambda 4686$) variations, combined with other measurements, yield an orbital period $2022.7\pm0.3$ d. The observed variability of the EW(He II $\lambda 4686$) was reproduced by a model in which the line flux primarily arises at the apex of the wind-wind collision and scales inversely with the square of the stellar separation, if we account for the excess emission as the companion star plunges into the hot inner layers of the primary's atmosphere, and including absorption from the disturbed primary wind between the source and the observer. This model constrains the orbital inclination to $135^\circ$-$153^\circ$, and the longitude of periastron to $234^\circ$-$252^\circ$. It also suggests that periastron passage occurred on $T_0 = 2456874.4\pm1.3$ d. Our model also reproduced EW(He II $\lambda 4686$) variations from a polar view of the primary star as determined from the observed He II $\lambda 4686$ emission scattered off the Homunculus nebula.

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.

The SDSS III APOGEE survey recently identified two new $\sigma$ Ori E type candidates, HD 345439 and HD 23478, which are a rare subset of rapidly rotating massive stars whose large (kGauss) magnetic fields confine circumstellar material around these systems. Our analysis of multi-epoch photometric observations of HD 345439 from the KELT, SuperWASP, and ASAS surveys reveals the presence of a $\sim$0.7701 day period in each dataset, suggesting the system is amongst the faster known $\sigma$ Ori E analogs. We also see clear evidence that the strength of H-alpha, H I Brackett series lines, and He I lines also vary on a $\sim$0.7701 day period from our analysis of multi-epoch, multi-wavelength spectroscopic monitoring of the system from the APO 3.5m telescope. We trace the evolution of select emission line profiles in the system, and observe coherent line profile variability in both optical and infrared H I lines, as expected for rigidly rotating magnetosphere stars. We also analyze the evolution of the H I Br-11 line strength and line profile in multi-epoch observations of HD 23478 from the SDSS-III APOGEE instrument. The observed periodic behavior is consistent with that recently reported by Sikora and collaborators in optical spectra.

Eta Carinae is the nearest example of a supermassive, superluminous, unstable star. Mass loss from the system is critical in shaping its circumstellar medium and in determining its ultimate fate. Eta Car currently loses mass via a dense, slow stellar wind and possesses one of the largest mass loss rates known. It is prone to episodes of extreme mass ejection via eruptions from some as-yet unspecified cause; the best examples of this are the large-scale eruptions which occurred in 19th century. Eta Car is a colliding wind binary in which strong variations in X-ray emission and in other wavebands are driven by the violent collision of the wind of eta Car-A and the fast, less dense wind of an otherwise hidden companion star. X-ray variations are the simplest diagnostic we have to study the wind-wind collision and allow us to measure the state of the stellar mass loss from both stars. We present the X-ray lightcurve over the last 20 years from ROSAT observations and monitoring with the Rossi X-ray Timing Explorer and the X-ray Telescope on the Swift satellite. We compare and contrast the behavior of the X-ray emission from the system over that timespan, including surprising variations during the 2014 X-ray minimum.

We report on high-resolution spectroscopy of the 2009.0 spectroscopic event of $\eta$ Carinae collected via SMARTS observations using the CTIO 1.5 m telescope and echelle spectrograph. Our observations were made almost every night over a two-month interval around the photometric minimum of $\eta$ Car associated with the periastron passage of a hot companion. The photoionizing flux of the companion and heating related to colliding winds causes large changes in the wind properties of the massive primary star. Here we present an analysis of temporal variations in a sample of spectral lines that are clearly formed in the wind of the primary star. These lines are affected by a changing illumination of the flux of the secondary star during the periastron passage. We document the sudden onset of blue-shifted absorption that occurred in most of the lines near or slightly after periastron, and we argue that these absorption components are seen when we view the relatively undisturbed wind of the foreground primary star. We present time series measurements of the net equivalent width of the wind lines and of the radial velocities of the absorption trough minima and the emission peak midpoints. Most lines decrease in emission strength around periastron, and those high excitation lines formed close to the primary exhibit a red-ward velocity excursion. We show how these trends can be explained using an illuminated hemisphere model that is based on the idea that the emission originates primarily from the side of the primary facing the hot companion.

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.

The multiple star system of delta Orionis is one of the closest examples of a system containing a luminous O-type, bright giant star (component Aa1). It is often used as a spectral-type standard and has the highest observed X-ray flux of any hot-star binary. The main component Aa1 is orbited by two lower mass stars, faint Aa2 in a 5.7 day eclipsing binary, and Ab, an astrometric companion with an estimated period of 346 years. Generally the flux from all three stars is recorded in ground-based spectroscopy, and the spectral decomposition of the components has proved difficult. Here we present HST/STIS ultraviolet spectroscopy of delta Ori A that provides us with spatially separated spectra of Aa and Ab for the first time. We measured radial velocities for Aa1 and Ab in two observations made near the velocity extrema of Aa1. We show tentative evidence for the detection of the Aa2 component in cross-correlation functions of the observed and model spectra. We discuss the appearance of the UV spectra of Aa1 and Ab with reference to model spectra. Both stars have similar effective temperatures, but Ab is fainter and is a rapid rotator. The results will help in the interpretation of ground-based spectroscopy and in understanding the physical and evolutionary parameters of these massive stars.

A classical nova occurs when material accreting onto the surface of a white dwarf in a close binary system ignites in a thermonuclear runaway. Complex structures observed in the ejecta at late stages could result from interactions with the companion during the common envelope phase. Alternatively, the explosion could be intrinsically bipolar, resulting from a localized ignition on the surface of the white dwarf or as a consequence of rotational distortion. Studying the structure of novae during the earliest phases is challenging because of the high spatial resolution needed to measure their small sizes. Here we report near-infrared interferometric measurements of the angular size of Nova Delphini 2013, starting from one day after the explosion and continuing with extensive time coverage during the first 43 days. Changes in the apparent expansion rate can be explained by an explosion model consisting of an optically thick core surrounded by a diffuse envelope. The optical depth of the ejected material changes as it expands. We detect an ellipticity in the light distribution, suggesting a prolate or bipolar structure that develops as early as the second day. Combining the angular expansion rate with radial velocity measurements, we derive a geometric distance to the nova of 4.54 +/- 0.59 kpc from the Sun.

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 present the results of an all-sky survey made with the Fine Guidance Sensor on Hubble Space Telescope to search for angularly resolved binary systems among the massive stars. The sample of 224 stars is comprised mainly of Galactic O- and B-type stars and Luminous Blue Variables, plus a few luminous stars in the Large Magellanic Cloud. The FGS TRANS mode observations are sensitive to detection of companions with an angular separation between 0."01 and 1."0 and brighter than $\triangle m = 5$. The FGS observations resolved 52 binary and 6 triple star systems and detected partially resolved binaries in 7 additional targets (43 of these are new detections). These numbers yield a companion detection frequency of 29% for the FGS survey. We also gathered literature results on the numbers of close spectroscopic binaries and wider astrometric binaries among the sample, and we present estimates of the frequency of multiple systems and the companion frequency for subsets of stars residing in clusters and associations, field stars, and runaway stars. These results confirm the high multiplicity fraction, especially among massive stars in clusters and associations. We show that the period distribution is approximately flat in increments of log P. We identify a number of systems of potential interest for long term orbital determinations, and we note the importance of some of these companions for the interpretation of the radial velocities and light curves of close binaries that have third companions.

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.

We present the first high angular resolution observations in the nearinfrared H-band (1.6 microns) of the Luminous Blue Variable star P Cygni. We obtained six-telescope interferometric observations with the CHARA Array and the MIRC beam combiner. These show that the spatial flux distribution is larger than expected for the stellar photosphere. A two component model for the star (uniform disk) plus a halo (two-dimensional Gaussian) yields an excellent fit of the observations, and we suggest that the halo corresponds to flux emitted from the base of the stellar wind. This wind component contributes about 45% of the H-band flux and has an angular FWHM = 0.96 mas, compared to the predicted stellar diameter of 0.41 mas. We show several images reconstructed from the interferometric visibilities and closure phases, and they indicate a generally spherical geometry for the wind. We also obtained near-infrared spectrophotometry of P Cygni from which we derive the flux excess compared to a purely photospheric spectral energy distribution. The H-band flux excess matches that from the wind flux fraction derived from the two component fits to the interferometry. We find evidence of significant near-infrared flux variability over the period from 2006 to 2010 that appears similar to the variations in the H-alpha emission flux from the wind. Future interferometric observations may be capable of recording the spatial variations associated with temporal changes in the wind structure.

We report on a high angular resolution survey of circumstellar disks around 24 northern sky Be stars. The K-band continuum survey was made using the CHARA Array long baseline interferometer (baselines of 30 to 331 m). The interferometric visibilities were corrected for the flux contribution of stellar companions in those cases where the Be star is a member of a known binary or multiple system. For those targets with good uv coverage, we used a four parameter Gaussian elliptical disk model to fit the visibilities and to determine the axial ratio, position angle, K-band photospheric flux contribution, and angular diameter of the disk major axis. For the other targets with relatively limited uv coverage, we constrained the axial ratio, inclination angle, and or disk position angle where necessary in order to resolve the degeneracy between possible model solutions. We also made fits of the ultraviolet and infrared spectral energy distributions to estimate the stellar angular diameter and infrared flux excess of each target. The mean ratio of the disk diameter (measured in K-band emission) to stellar diameter (from SED modeling) is 4.4 among the 14 cases where we reliably resolved the disk emission, a value which is generally lower than the disk size ratio measured in the higher opacity Halpha emission line. We estimated the equatorial rotational velocity from the projected rotational velocity and disk inclination for 12 stars, and most of these stars rotate close to or at the critical rotational velocity.

\gamma Cas is the prototypical classical Be star and is best known for its variable hard X-ray emission. To elucidate the reasons for this emission, we mounted a multiwavelength campaign in 2010 centered around 4 XMM observations. The observational techniques included long baseline optical interferometry (LBOI), monitoring by an Automated Photometric Telescope and Halpha observations. Because gamma Cas is also known to be in a binary, we measured Halpha radial velocities and redetermined its period as 203.55+/-0.2 days and an eccentricity near zero. The LBOI observations suggest that the star's decretion disk was axisymmetric in 2010, has an inclination angle near 45^o, and a larger radius than previously reported. The Be star began an "outburst" at the beginning of our campaign, made visible by a disk brightening and reddening during our campaign. Our analyses of the new high resolution spectra disclosed many attributes found from spectra obtained in 2001 (Chandra) and 2004 (XMM). As well as a dominant hot 14 keV thermal component, these familiar ones included: (i) a fluorescent feature of Fe K stronger than observed at previous times, (ii) strong lines of N VII and Ne XI lines indicative of overabundances, and (iii) a subsolar Fe abundance from K-shell lines but a solar abundance from L-shell ions. We also found that 2 absorption columns are required to fit the continuum. While the first one maintained its historical average of 1X10^21 cm^-2, the second was very large and doubled to 7.4X10^23 cm^-2 during our X-ray observations. Although we found no clear relation between this column density and orbital phase, it correlates well with the disk brightening and reddening both in the 2010 and earlier observations. Thus, the inference from this study is that much (perhaps all?) of the X-ray emission from this source originates behind matter ejected by gamma Cas into our line of sight.

We have executed a survey of nearby, main sequence A, F, and G-type stars with the CHARA Array, successfully measuring the angular diameters of fortyfour stars with an average precision of ~ 1.5%. We present new measures of the bolometric flux, which in turn leads to an empirical determination of the effective temperature for the stars observed. In addition, these CHARA-determined temperatures, radii, and luminosities are fit to Yonsei-Yale model isochrones to constrain the masses and ages of the stars. These results are compared to indirect estimates of these quantities obtained by collecting photometry of the stars and applying them to model atmospheres and evolutionary isochrones. We find that for most cases, the models overestimate the effective temperature by ~ 1.5-4%, when compared to our directly measured values. The overestimated temperatures and underestimated radii in these works appear to cause an additional offset in the star's surface gravity measurements, which consequently yield higher masses and younger ages, in particular for stars with masses greater than ~1.3 Msun. Additionally, we compare our measurements to a large sample of eclipsing binary stars, and excellent agreement is seen within both data sets. Finally, we present temperature relations with respect to (B-V) and (V-K) color as well as spectral type showing that calibration of effective temperatures with errors ~ 1% is now possible from interferometric angular diameters of stars.

The hot star HDE 326823 is a candidate transition-phase object that is evolving into a nitrogen-enriched Wolf-Rayet star. It is also a known low-amplitude, photometric variable with a 6.123 d period. We present new, high and moderate resolution spectroscopy of HDE 326823, and we show that the absorption lines show coherent Doppler shifts with this period while the emission lines display little or no velocity variation. We interpret the absorption line shifts as the orbital motion of the apparently brighter star in a close, interacting binary. We argue that this star is losing mass to a mass gainer star hidden in a thick accretion torus and to a circumbinary disk that is the source of the emission lines. HDE 326823 probably belongs to a class of objects that produce short-period WR+O binaries.

We report on the preliminary results of a three-year monitoring campaign of Galactic and Magellanic LBVs with both moderate and high resolution spectroscopy. We have collected more than 500 moderate-resolution spectra of 20 stars in the red portion of the optical spectrum, which includes the wind-sensitive transitions of H$\alpha$ and He I 5876 and 6678 Angstroms. We summarize our detailed study of 15 years of high resolution H$\alpha$ spectroscopy of the prototypical luminous blue variable, P Cygni. We report on the discovery of discrete absorption components in P Cygni's H$\alpha$ profile, and we discuss their relationship to structure in the wind. These results are compared to our recent high resolution interferometric imaging of the $H$-band emitting region surrounding the star. We discuss recent observations of $\eta$ Car, highlighting its unusual recovery from the previous periastron passage. Our results on HDE 326823 indicate that the star has a short period binary orbit, with Roche Lobe overflow onto an unseen massive companion.

The CHARA Array is a six 1-m telescope optical and near infrared interferometer located at the Mount Wilson Observatory in southern California and operated by the Center for High Angular Resolution Astronomy of Georgia State University. The CHARA Array has been in regular scientific operation since 2005 and now has over 55 publications in the refereed literature, including two in Science and one in Nature. The Array now supports seven beam combiners ranging from 0.5 microns up to 2.3 microns and combing from 2 to 4 beams at a time. An upgrade to a full 6 beam combiner is now underway and fringes with all six telescopes were achieved soon after the meeting.We present some of the more recent results from the CHARA-Array.

P Cygni is a prototype of the Luminous Blue Variables (or S Doradus variables), and the star displays photometric and emission line variability on a timescale of years (known as the "short S Doradus phase" variations). Here we present new high resolution H-alpha spectroscopy of P Cyg that we combine with earlier spectra and concurrent V-band photometry to document the emission and continuum flux variations over a 24 y time span. We show that the emission and continuum fluxes vary in concert on timescales of 1.6 y and longer, but differ on shorter timescales. The H-alpha profile shape also varies on the photometric timescales, and we describe the observed co-variations of the emission peak and absorption trough properties. We argue that the episodes of photometric and emission brightening are caused by increases in the size of the emission region that are related to variations in wind mass loss rate and outflow speed. We find evidence of blueward accelerating, Discrete Absorption Components (DACs) in the absorption trough of the H-alpha profile, and these features have slower accelerations and longer durations than those observed in other lines. The DAC strengths also appear to vary on the photometric timescales, and we suggest that the propagation of the DAC-related wind structures is closely related to changes in the overall wind mass loss rate and velocity.

Deneb is often considered the prototypical A-type supergiant, and is one of the visually most luminous stars in the Galaxy. A-type supergiants are potential extragalactic distance indicators, but the variability of these stars needs to be better characterized before this technique can be considered reliable. We analyzed 339 high resolution echelle spectra of Deneb obtained over the five-year span of 1997 through 2001 as well as 370 Stromgren photometric measurements obtained during the same time frame. Our spectroscopic analysis included dynamical spectra of the H-alpha profile, H-alpha equivalent widths, and radial velocities measured from Si II 6347, 6371. Time-series analysis reveals no obvious cyclic behavior that proceeds through multiple observing seasons, although we found a suspected 40 day period in two, non-consecutive observing seasons. Some correlations are found between photometric and radial velocity data sets, and suggest radial pulsations at two epochs. No correlation is found between the variability of the H-alpha profiles and that of the radial velocities or the photometry. Lucy (1976) found evidence that Deneb was a long period single-lined spectroscopic binary star, but our data set shows no evidence for radial velocity variations caused by a binary companion.

We present interferometric observations of the Be star Zeta Tau obtained using the MIRC beam combiner at the CHARA Array. We resolved the disk during four epochs in 2007-2009. We fit the data with a geometric model to characterize the circumstellar disk as a skewed elliptical Gaussian and the central Be star as a uniform disk. The visibilities reveal a nearly edge-on disk with a FWHM major axis of ~ 1.8 mas in the H-band. The non-zero closure phases indicate an asymmetry within the disk. Interestingly, when combining our results with previously published interferometric observations of Zeta Tau, we find a correlation between the position angle of the disk and the spectroscopic V/R ratio, suggesting that the tilt of the disk is precessing. This work is part of a multi-year monitoring campaign to investigate the development and outward motion of asymmetric structures in the disks of Be stars.

We present spectrophotometric data from 0.4 to 4.2 microns for bright, northern sky, Be stars and several other types of massive stars. Our goal is to use these data with ongoing, high angular resolution, interferometric observations to model the density structure and sky orientation of the gas surrounding these stars. We also present a montage of the H-alpha and near-infrared emission lines that form in Be star disks. We find that a simplified measurement of the IR excess flux appears to be correlated with the strength of emission lines from high level transitions of hydrogen. This suggests that the near-IR continuum and upper level line fluxes both form in the inner part of the disk, close to the star.

We report on H-alpha spectroscopy of the 2009.0 spectroscopic event of eta Carinae collected via SMARTS observations using the CTIO 1.5 m telescope and echelle spectrograph. Our observations were made almost every night over a two month interval around the predicted minimum of eta Car. We observed a significant fading of the line emission that reached a minimum seven days after the X-ray minimum. About 17 d prior to the H-alpha flux minimum, the H-alpha profile exhibited the emergence of a broad, P Cygni type, absorption component (near a Doppler shift of -500 km/s) and a narrow absorption component (near -144 km/s and probably associated with intervening gas from the Little Homunculus Nebula). All these features were observed during the last event in 2003.5 and are probably related to the close periastron passage of the companion. We argue that these variations are consistent with qualitative expectations about changes in the primary star's stellar wind that result from the wind-wind collision with a massive binary companion and from atmospheric eclipses of the companion.

We present a radial velocity study of the rapidly rotating B-star Regulus that indicates the star is a single-lined spectroscopic binary. The orbital period (40.11 d) and probable semimajor axis (0.35 AU) are large enough that the system is not interacting at present. However, the mass function suggests that the secondary has a low mass (M_2 > 0.30 M_sun), and we argue that the companion may be a white dwarf. Such a star would be the remnant of a former mass donor that was the source of the large spin angular momentum of Regulus itself.