Massive stars like company. However, low-mass companions have remained extremely difficult to detect at angular separations ($\rho$) smaller than 1" (~1000-3000 au considering typical distance to nearby massive stars) given the large brightness contrast between the companion and the central star. Constraints on the low-mass end of the companions mass-function for massive stars are however needed, for example to help distinguishing between various scenarios for the formation of massive stars. To obtain statistically significant constraint on the presence of low-mass companions beyond the typical detection limit of current surveys ($\Delta \mathrm{mag} \lesssim 5$ at $\rho \lesssim 1$"), we initiated a survey of O and Wolf-Rayet stars in the Carina region using the SPHERE coronagraphic instrument on the VLT. In this first paper, we aim to introduce the survey, to present the methodology and to demonstrate the capability of SPHERE for massive stars using the multiple system QZ~Car. High-contrast imaging techniques, such as angular- and spectral-differential imaging techniques as well as PSF-fitting, were applied to detect and measure the relative flux of companions in each spectral channel of the instrument. We detected 19 sources around the QZ~Car system with detection limits of 9~mag at $\rho > 200$~mas for IFS and as faint as 13~mag at $\rho$ > 1".8 for IRDIS (corresponding to sub-solar masses for potential companions). All but two are reported here for the first time. Based on this proof of concept, we showed that VLT/SPHERE allows us to reach the sub-solar mass regime of the companion mass function. This paves the way for this type of observation with a large sample of massive stars to provide novel constraints on the multiplicity of massive stars in a region of the parameter space that has remained inaccessible so far.

Our multi-view metric learning framework enables robust characterization of star categories by directly learning to discriminate in a multi-faceted feature space, thus, eliminating the need to combine feature representations prior to fitting the machine learning model. We also demonstrate how to extend standard multi-view learning, which employs multiple vectorized views, to the matrix-variate case which allows very novel variable star signature representations. The performance of our proposed methods is evaluated on the UCR Starlight and LINEAR datasets. Both the vector and matrix-variate versions of our multi-view learning framework perform favorably --- demonstrating the ability to discriminate variable star categories.

We present the construction of a novel time-domain signature extraction methodology and the development of a supporting supervised pattern detection algorithm. We focus on the targeted identification of eclipsing binaries that demonstrate a feature known as the O'Connell effect. Our proposed methodology maps stellar variable observations to a new representation known as distribution fields (DFs). Given this novel representation, we develop a metric learning technique directly on the DF space that is capable of specifically identifying our stars of interest. The metric is tuned on a set of labeled eclipsing binary data from the Kepler survey, targeting particular systems exhibiting the O'Connell effect. The result is a conservative selection of 124 potential targets of interest out of the Villanova Eclipsing Binary Catalog. Our framework demonstrates favorable performance on Kepler eclipsing binary data, taking a crucial step in preparing the way for large-scale data volumes from next-generation telescopes such as LSST and SKA.

The Next Generation Very Large Array will have excellent sensitivity for detecting the thermal emission from massive stars and from red giants. This will allow direct imaging of the winds for a large number of hot massive stars. It will also allow using the radio emission for the massive stars as a way to detect stars to allow high resolution measurements can be made, even with large extinction. A few examples of the utility of the high resolution measurements are given: dynamics of globular clusters with red giants, detection of intermediate mass stripped stars in binaries, and measurement of masses of stars in massive binaries.

Massive stars can be found in wide (hundreds to thousands AU) binaries with other massive stars. We use $N$-body simulations to show that any bound cluster should always have approximately one massive wide binary: one will probably form if none are present initially; and probably only one will survive if more than one are present initially. Therefore any region that contains many massive wide binaries must have been composed of many individual subregions. Observations of Cyg OB2 show that the massive wide binary fraction is at least a half (38/74) which suggests that Cyg OB2 had at least 30 distinct massive star formation sites. This is further evidence that Cyg OB2 has always been a large, low-density association. That Cyg OB2 has a normal high-mass IMF for its total mass suggests that however massive stars form they 'randomly sample' the IMF (as the massive stars did not 'know' about each other).

We present VLT/MUSE observations of NGC 2070, the dominant ionizing nebula of 30 Doradus in the LMC, plus HST/STIS spectroscopy of its central star cluster R136. Integral Field Spectroscopy (MUSE) and pseudo IFS (STIS) together provides a complete census of all massive stars within the central 30x30 parsec^2 of the Tarantula. We discuss the integrated far-UV spectrum of R136, of particular interest for UV studies of young extragalactic star clusters. Strong HeII 1640 emission at very early ages (1-2 Myr) from very massive stars cannot be reproduced by current population synthesis models, even those incorporating binary evolution and very massive stars. A nebular analysis of the integrated MUSE dataset implies an age of ~4.5 Myr for NGC 2070. Wolf-Rayet features provide alternative age diagnostics, with the primary contribution to the integrated Wolf-Rayet bumps arising from R140 rather than the more numerous H-rich WN stars in R136. Caution should be used when interpreting spatially extended observations of extragalactic star-forming regions.

We introduce a HST/STIS stellar census of R136a, the central ionizing star cluster of 30 Doradus. We present low resolution far-ultraviolet STIS/MAMA spectroscopy of R136 using 17 contiguous 52x0.2 arcsec slits which together provide complete coverage of the central 0.85 parsec (3.4 arcsec). We provide spectral types of 90% of the 57 sources brighter than m_F555W = 16.0 mag within a radius of 0.5 parsec of R136a1, plus 8 additional nearby sources including R136b (O4\,If/WN8). We measure wind velocities for 52 early-type stars from CIV 1548-51, including 16 O2-3 stars. For the first time we spectroscopically classify all Weigelt & Baier members of R136a, which comprise three WN5 stars (a1-a3), two O supergiants (a5-a6) and three early O dwarfs (a4, a7, a8). A complete Hertzsprung-Russell diagram for the most massive O stars in R136 is provided, from which we obtain a cluster age of 1.5+0.3_-0.7 Myr. In addition, we discuss the integrated ultraviolet spectrum of R136, and highlight the central role played by the most luminous stars in producing the prominent HeII 1640 emission line. This emission is totally dominated by very massive stars with initial masses above ~100 Msun. The presence of strong HeII 1640 emission in the integrated light of very young star clusters (e.g A1 in NGC 3125) favours an initial mass function extending well beyond a conventional upper limit of 100 Msun. We include montages of ultraviolet spectroscopy for LMC O stars in the Appendix. Future studies in this series will focus on optical STIS/CCD medium resolution observations.

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 present results of a high angular resolution survey of massive OB stars in the Cygnus OB2 association that we conducted with the Fine Guidance Sensor 1R (FGS1r) on the Hubble Space Telescope. FGS1r is able to resolve binary systems with a magnitude difference delta-V < 4 down to separations as small as 0.01 arcsec. The sample includes 58 of the brighter members of Cyg OB2, one of the closest examples of an environment containing a large number of very young and massive stars. We resolved binary companions for 12 targets and confirmed the triple nature of one other target, and we offer evidence of marginally resolved companions for two additional stars. We confirm the binary nature of 11 of these systems from complementary adaptive optics imaging observations. The overall binary frequency in our study is 22% to 26% corresponding to orbital periods ranging from 20 - 20,000 years. When combined with the known short-period spectroscopic binaries, the results supports the hypothesis that the binary fraction among massive stars is > 60%. One of the new discoveries is a companion to the hypergiant star MT 304 = Cyg OB2-12, and future measurements of orbital motion should provide mass estimates for this very luminous star.

We present a catalog of 56 candidate intermediate mass eclipsing binary systems extracted from the 3rd data release of the All Sky Automated Survey. We gather pertinent observational data and derive orbital properties, including ephemerides, for these systems as a prelude to anticipated spectroscopic observations. We find that 37 of the 56, or ~66% of the systems are not identified in the Simbad Astronomical Database as known binaries. As a specific example, we show spectroscopic data obtained for the system HI Mon (B0 V + B0.5 V) observed at key orbital phases based on the computed ephemeris, and we present a combined spectroscopic and photometric solution for the system and give stellar parameters for each component.

We present an examination of high resolution, ultraviolet spectroscopy from Hubble Space Telescope of the photospheric spectrum of the O-supergiant in the massive X-ray binary HD 226868 = Cyg X-1. We analyzed this and ground-based optical spectra to determine the effective temperature and gravity of the O9.7 Iab supergiant. Using non-local thermodynamic equilibrium (non-LTE), line blanketed, plane parallel models from the TLUSTY grid, we obtain T_eff = 28.0 +/- 2.5kK and log g > 3.00 +/- 0.25, both lower than in previous studies. The optical spectrum is best fit with models that have enriched He and N abundances. We fit the model spectral energy distribution for this temperature and gravity to the UV, optical, and IR fluxes to determine the angular size of and extinction towards the binary. The angular size then yields relations for the stellar radius and luminosity as a function of distance. By assuming that the supergiant rotates synchronously with the orbit, we can use the radius - distance relation to find mass estimates for both the supergiant and black hole as a function of the distance and the ratio of stellar to Roche radius. Fits of the orbital light curve yield an additional constraint that limits the solutions in the mass plane. Our results indicate masses of 23^{+8}_{-6} M_sun for the supergiant and 11^{+5}_{-3} M_sun for the black hole.

We present results from Hubble Space Telescope UV spectroscopy of the massive X-ray binary system, HD226868 = Cyg X-1. The spectra were obtained at both orbital conjunction phases in two separate runs in 2002 and 2003 when the system was in the X-ray high/soft state. The stellar wind lines suffer large reductions in strength when the black hole is in the foreground due to the X-ray ionization of the wind ions. A comparison of HST and archival IUE spectra shows that similar photoionization effects occur in both the X-ray states. We constructed model UV wind line profiles assuming that X-ray ionization occurs everywhere in the wind except the zone where the supergiant blocks the X-ray flux. The good match between the observed and model profiles indicates that the wind ionization extends to near to the hemisphere of the supergiant facing the X-ray source. The H-alpha emission strength is generally lower in the high/soft state compared to the low/hard state, but the He II 4686 emission is relatively constant between states. The results suggest that mass transfer in Cyg X-1 is dominated by a focused wind flow that peaks along the axis joining the stars and that the stellar wind contribution is shut down by X-ray photoionization effects. The strong stellar wind from the shadowed side of the supergiant will stall when Coriolis deflection brings the gas into the region of X-ray illumination. This stalled gas component may be overtaken by the orbital motion of the black hole and act to inhibit accretion from the focused wind. The variations in the strength of the shadow wind component may then lead to accretion rate changes that ultimately determine the X-ray state.

We use the first data release from the SINGG H-alpha survey of HI-selected galaxies to study the quantitative behavior of the diffuse, warm ionized medium (WIM) across the range of properties represented by these 109 galaxies. The mean fraction f_WIM of diffuse ionized gas in this sample is 0.59+/- 0.19, slightly higher than found in previous samples. Since lower surface-brightness galaxies tend to have higher f_WIM, we believe that most of this difference is due to selection effects favoring large, optically-bright, nearby galaxies with high star-formation rates. As found in previous studies, there is no appreciable correlation with Hubble type or total star-formation rate. However, we find that starburst galaxies, defined here by an H-alpha surface brightness > 2.5x 10^39 erg s^-1 kpc^-2 within the H-alpha half-light radius, do show much lower fractions of diffuse H-alpha emission. The cause apparently is not dominated by a lower fraction of field OB stars. However, it is qualitatively consistent with an expected escape of ionizing radiation above a threshold star-formation rate, predicted from our model in which the ISM is shredded by pressure-driven supernova feedback. The HI gas fractions in the starburst galaxies are also lower, suggesting that the starbursts are consuming and ionizing all the gas, and thus promoting regions of density-bounded ionization. If true, these effects imply that some amount of Lyman continuum radiation is escaping from most starburst galaxies, and that WIM properties and outflows from mechanical feedback are likely to be pressure-driven. However, in view of previous studies showing that the escape fraction of ionizing radiation is generally low, it is likely that other factors also drive the low fractions of diffuse ionized gas in starbursts.

We present the results of an H-alpha monitoring campaign on the BeXRB and microquasar system LS I +61 303. We use radial velocity measurements of HeI lines in our spectra to re-evaluate the orbital elements and to better establish the time of periastron. We list equivalent widths and other parameters for the H-alpha emission line and discuss the orbital phase related variations observed. We call attention to a dramatic episode of emission weakening that occurred in less than a day that probably resulted from exposure to a transient source of ionizing radiation. We argue that the increase in H-alpha and X-ray emission following periastron probably results from the creation of an extended density wave in the disk created by tidal forces. We also discuss estimates of the size of the disk from the H-alpha equivalent width measurements, and we suggest that the disk radius from the average equivalent width corresponds to a resonant truncation radius of the disk while the maximum equivalent width corresponds to a radius limited by the separation of the stars at periastron. We note that a nearby faint companion is probably an unrelated foreground object.