We have started an ambitious program to determine if the full diversity of extinction laws is real or if some of it is due to calibration or methodological issues. Here we start by analyzing the information on NIR extinction in a 2MASS stellar sample with good quality photometry and very red colours. We calculate the extinction at 1 $\mu$m, $A_1$, and the power-law exponent, $\alpha$ ($A_\lambda = A_1 \lambda^{-\alpha}$), for the 2MASS stars located in the extinction trajectory in the $H-K$ vs. $J-H$ plane expected for red giants with $A_1 > 5$ mag. We test the validity of the assumption about the nature of those stars, whether a single or multiple values of $\alpha$ are needed, and the spatial variations of the results. Most ($\sim$83%) of those stars can indeed be explained by high-extinction red giants and the rest is composed of extinguished AGB stars (mostly O-rich), blended sources, and smaller numbers of other objects, a contaminant fraction that can be reduced with the help of Gaia} DR2 data. Galactic red giants experience a NIR extinction with $\alpha\sim 2.27$ and an uncertainty of a few hundredths of a magnitude. There is no significant spread in $\alpha$ even though our sample is widely distributed and has a broad range of extinctions. Differences with previous results are ascribed to the treatment of non-linear photometric effects and/or the contaminant correction. Future research should concentrate in finding the correct functional form for the NIR extinction law. In the appendix we detail the treatment of non-linear photometric effects in the 2MASS bands.

We conducted an observational campaign towards one of the most massive and luminous colliding wind binaries in the Galaxy, HD~93129A, close to its periastron passage in 2018. During this time the source was predicted to be in its maximum of high-energy emission. Here we present our data analysis from the X-ray satellites \textit{Chandra} and \textit{NuSTAR} and the $\gamma$-ray satellite \textit{AGILE}. High-energy emission coincident with HD~93129A was detected in the X-ray band up to $\sim$18~keV, whereas in the $\gamma$-ray band only upper limits were obtained. We interpret the derived fluxes using a non-thermal radiative model for the wind-collision region. We establish a conservative upper limit for the fraction of the wind kinetic power that is converted into relativistic electron acceleration, $f_\mathrm{NT,e} < 0.02$. In addition, we set a lower limit for the magnetic field in the wind-collision region as $B_\mathrm{WCR} > 0.3$~G. We also argue a putative interpretation of the emission from which we estimate $f_\mathrm{NT,e} \approx 0.006$ and $B_\mathrm{WCR} \approx 0.5$~G. We conclude that multi-wavelength, dedicated observing campaigns during carefully selected epochs are a powerful tool for characterising the relativistic particle content and magnetic field intensity in colliding wind binaries.

CONTEXT. Many O-type stars have nearby companions whose presence hamper their characterization through spectroscopy. AIMS. We want to obtain spatially resolved spectroscopy of close O-type visual binaries to derive their spectral types. METHODS. We use the Space Telescope Imaging Spectrograph (STIS) of the Hubble Space Telescope (HST) to obtain long-slit blue-violet spectroscopy of eight Galactic O-type stars with nearby visual companions and use spatial-profile fitting to extract the separate spectra. We also use the ground-based Galactic O-Star Spectroscopic Survey (GOSSS) to study more distant visual components. RESULTS. We spatially resolve seven of the eight systems, present spectra for their components, and obtain their spectral types. Those seven multiple systems are iota Ori Aa,Ab,B, 15 Mon Aa,Ab,C, tau CMa Aa,Ab,B,C,D,E, HD 206 267 Aa,Ab,C,D, HD 193 443 A,B, HD 16 429 Aa,Ab, and IU Aur A,B. This is the first time that spatially resolved spectroscopy of the close visual binaries of those systems is obtained. We establish the applicability of the technique as a function of the separation and magnitude difference of the binary.

Binary systems composed of a red supergiant and a B star are useful probes of stellar evolution. We have searched the literature to create a catalog of 108 Galactic systems of such type, which is presented here.

LB-1 has recently been proposed to be a binary system at 4 kpc consisting of a B star of 8 Msol and a massive stellar black hole of 70 Msol. This finding challenges our current theories of massive star evolution and formation of BHs at solar metallicity. Our objective is to derive the effective temperature, surface gravity and chemical composition of the B-type component in order to determine its nature and evolutionary status and, indirectly, to constrain the mass of the BH. We use the non-LTE stellar atmosphere code FASTWIND to analyse new and archival high resolution data. We determine (Teff, logg) values of (14000$\pm500$ K, 3.50$\pm0.15$ dex) that, combined with the Gaia parallax, implies a spectroscopic mass, from logg, of $3.2^{+2.1}_{-1.9}$ Msol and an evolutionary mass, assuming single star evolution, of $5.2^{+0.3}_{-0.6}$ Msol. We determine an upper limit of 8 km/s for the projected rotational velocity and derive the surface abundances, finding the star to have a silicon abundance below solar, to be significantly enhanced in nitrogen and iron, and depleted in carbon and magnesium. Complementary evidence derived from a photometric extinction analysis and Gaia yields similar results for Teff and logg and a consistent distance around 2~kpc. We propose that the B star is a slightly evolved main sequence star of 3-5 Msol with surface abundances reminiscent of diffusion in late B/A chemically peculiar stars with low rotational velocities. There is also evidence for CN-processed material in its atmosphere. These conclusions rely critically on the distance inferred from the Gaia parallax. The goodness of fit of the Gaia astrometry also favours a high-inclination orbit. If the orbit is edge-on and the B star has a mass of 3-5 Msol, the mass of the dark companion would be 4-5 Msol, which would be easier to explain with our current stellar evolutionary models.

[ABRIDGED] CONTEXT: Gaia DR2 has opened a trove of astrometric and photometric data for Galactic clusters close to the Sun. Lucky imaging has been an operational technique to measure the positions of visual binary systems for a decade and a half, enough to apply its results to the calculation of orbits of some massive multiple systems within ~1 kpc of the Sun. AIMS: We are measuring distances to Galactic stellar groups containing O stars and I start with two of them: Collinder 419 in Cygnus and NGC 2264 in Monoceros. I also aim to derive new astrometric orbits for the Aa,Ab components for the main ionizing sources for both clusters: HD 193 322 and 15 Mon, respectively. METHODS: First, I present a method that uses Gaia DR2 photometry, positions, proper motions, and parallaxes to obtain the membership and distance of a stellar group and apply it to Collinder 419 and NGC 2264. Second, I present a new code that calculates astrometric orbits by searching the whole 7-parameter orbit space and apply it to HD 193 322 Aa,Ab and 15 Mon Aa,Ab using as input literature data from the Washington Double Star Catalog (WDS) and the AstraLux measurements recently presented by Ma\'iz Apell\'aniz et al. (2019). RESULTS: I obtain Gaia DR2 distances of 1006+37-34 pc for Collinder 419 and 719+-16 pc for NGC 2264, with the main contribution to the uncertainties coming from the spatial covariance of the parallaxes. The two NGC 2264 subclusters are at the same distance (within the uncertainties) and they show a significant relative proper motion. The distances are shown to be robust. HD 193 322 Aa,Ab follows an eccentric (e = 0.58+0.03-0.04) orbit with a period of 44+-1 a and the three stars it contains have a total mass of 76.1+9.9-7.4 M_Sol. The orbit of 15 Mon Aa,Ab is even more eccentric (e = 0.770+0.023-0.030), with a period of 108+-12 a and a total mass of 45.1+3.6-3.3 M_Sol for its two stars.

We present the photometric calibration of the twelve optical passbands observed by the Javalambre Photometric Local Universe Survey (J-PLUS). The proposed calibration method has four steps: (i) definition of a high-quality set of calibration stars using Gaia information and available 3D dust maps; (ii) anchoring of the J-PLUS gri passbands to the Pan-STARRS photometric solution, accounting for the variation of the calibration with the position of the sources on the CCD; (iii) homogenization of the photometry in the other nine J-PLUS filters using the dust de-reddened instrumental stellar locus in (X - r) versus (g - i) colours, where X is the filter to calibrate. The zero point variation along the CCD in these filters was estimated with the distance to the stellar locus. Finally, (iv) the absolute colour calibration was obtained with the white dwarf locus. We performed a joint Bayesian modelling of eleven J-PLUS colour-colour diagrams using the theoretical white dwarf locus as reference. This provides the needed offsets to transform instrumental magnitudes to calibrated magnitudes outside the atmosphere. The uncertainty of the J-PLUS photometric calibration, estimated from duplicated objects observed in adjacent pointings and accounting for the absolute colour and flux calibration errors, are ~19 mmag in u, J0378 and J0395, ~11 mmag in J0410 and J0430, and ~8 mmag in g, J0515, r, J0660, i, J0861, and z. We present an optimized calibration method for the large area multi-filter J-PLUS project, reaching 1-2% accuracy within an area of 1 022 square degrees without the need for long observing calibration campaigns or constant atmospheric monitoring. The proposed method will be adapted for the photometric calibration of J-PAS, that will observe several thousand square degrees with 56 narrow optical filters.

[ABRIDGED] AIMS. MONOS is collecting information and studying O-type spectroscopic binaries with delta > -20 deg. In this 1st paper we describe the sample and provide spectral classifications and additional information for objects with previous SB or EB orbits. In the future we will test the validity of previous solutions and calculate new SB orbits. METHODS.The spectra in this paper have 2 sources: GOSSS, which is obtaining blue-violet R~2500 spectroscopy of thousands of massive stars, and LiLiMaRlin, a library of libraries of high-resolution spectroscopy of massive stars obtained from 4 different surveys and additional data from our own observing programs and public archives. We also use lucky images from AstraLux. RESULTS. We present homogeneous spectral classifications for 92 O-type spectroscopic multiple systems and 10 optical companions. We discuss the visual multiplicity of each system with the support of AstraLux images and additional sources. For 11 O-type objects and for 6 B-type objects we present their first GOSSS spectral classifications. For 2 known EBs we detect SB2/SB1 lines for the first time, to which we add a third system already reported by us recently. For 2 previous SB1 systems we detect their SB2 nature for the first time and give their first separate spectral classifications, something we also do for a 3rd object just recently identified as a SB2. We also detect 9 new astrometric companions and provide updated information on several others. For sigma Ori AaAbB we provide spectral classifications for the 3 components with a single observation for the first time thanks to a lucky spectroscopy observation obtained close to the Aa,Ab periastron and for theta^1 Ori CaCb we add it to the class of Galactic Of?p stars, raising the number of its members to 6. Our sample of O-type spectroscopic binaries contains more triple/higher-order systems than double systems.

This paper describes the characterization of the GALANTE photometric system, a seven intermediate- and narrow-band filter system with a wavelength coverage from 3000 $\r{A}$ to 9000 $\r{A}$ . We describe the photometric system presenting the full sensitivity curve as a product of the filter sensitivity, CCD, telescope mirror, and atmospheric transmission curves, as well as some first- and second-order moments of this sensitivity function. The GALANTE photometric system is composed of four filters from the J-PLUS photometric system, a twelve broad-to-narrow filter system, and three exclusive filters, specifically designed to measure the physical parameters of stars such as effective temperature $T_{\rm eff}$, $\log(g)$, metallicity, colour excess $E(4405-5495)$, and extinction type $R_{5495}$. Two libraries, the Next Generation Spectral Library (NGSL) and the one presented in Ma\'iz Apell\'aniz & Weiler (2018), have been used to determine the transformation equations between the Sloan Digital Sky Survey ($\textit{SDSS}$) $\textit{ugriz}$ photometry and the GALANTE photometric system. We will use this transformation to calibrate the zero points of GALANTE images. To this end, a preliminary photometric calibration of GALANTE has been made based on two different $\textit{griz}$ libraries ($\textit{SDSS}$ DR12 and ATLAS All-Sky Stellar Reference Catalog, hereinafter $\textit{RefCat2}$). A comparison between both zero points is performed leading us to the choice of $\textit{RefCat2}$ as the base catalogue for this calibration, and applied to a field in the Cyg OB2 association.

We present a detailed spectroscopic study of Herschel 36 A (H36A), the main stellar component of the massive multiple system Herschel 36 in the Hourglass Nebula, based on high-resolution optical spectra obtained along an 11 years span. The three stellar components present in the spectrum of H36A are separated by means of a spectral disentangling technique. Individual spectral classifications are improved, and high precision orbital solutions for the inner and the outer orbits are calculated. H36A is confirmed to be a hierarchical triple system composed of a close massive binary (Ab1+Ab2, O9.5 V+B0.7 V) in wide orbit around a third O-type star (Aa, O7.5 Vz). The inner-pair orbit is characterized by a period of 1.54157 +/- 0.00006 days, and semi-amplitudes of 181.2 +/- 0.7 and 295.4 +/- 1.7 km/s. The outer orbit has a period of 492.81 +/- 0.69 days, and semi-amplitudes of 62.0 +/- 0.6 and 42.4 +/- 0.8 km/s. Inner and outer orbits are not coplanar, having a relative inclination of at least 20 degrees. Dynamical minimum masses of 20.6 +/- 0.8 Msun, 18.7 +/- 1.1 Msun, and 11.5 +/- 1.1 Msun are derived for the Aa, Ab1, and Ab2 components, respectively, in reasonable agreement with the theoretical calibrations.

In this paper, we present high-resolution HST/COS observations of the extreme magnetic O star NGC 1624-2. These represent the first ultraviolet spectra of this archetypal object. We examine the variability of its wind-sensitive resonance lines, comparing it to that of other known magnetic O stars. In particular, the observed variations in the profiles of the CIV and SiIV doublets between low state and high state are the largest observed in any magnetic O-type star, consistent with the expected properties of NGC 1624-2's magnetosphere. We also observe a redshifted absorption component in the low state, a feature not seen in most stars. We present preliminary modelling efforts based on the Analytic Dynamical Magnetosphere (ADM) formalism, demonstrating the necessity of using non-spherically symmetric models to determine wind/magnetospheric properties of magnetic O stars.

GALANTE is an optical photometric survey with seven intermediate/narrow filters that has been covering the Galactic Plane since 2016 using the Javalambre T80 and Cerro Tololo T80S telescopes. The P.I.s of the northern part (GALANTE NORTE) are Emilio J. Alfaro and Jes\'us Ma\'iz Apell\'aniz. and the P.I. of the southern part (GALANTE SUR) is Rodolfo H. Barb\'a. The detector has a continuous 1.4 degr x 1.4 degr field of view with a sampling of 0.55"/pixel and the seven filters are optimized to detect obscured early-type stars. The survey includes long, intermediate, short, and ultrashort exposure times to reach a dynamical range close to 20 magnitudes, something never achieved for such an optical project before. The characteristics of GALANTE allow for a new type of calibration scheme using external Gaia, Tycho-2, and 2MASS data that has already led to a reanalysis of the sensitivity of the Gaia G filter. We describe the project and present some early results. GALANTE will identify the majority of the early-type massive stars within several kpc of the Sun and measure their amount and type of extinction. It will also map the Halpha nebular emission, identify emission-line stars, and do other studies of low- and intermediate-mass stars.

We present results on two different Gaia-related topics. First, we describe our efforts to calibrate the three Gaia photometric passbands G, G_BP, and G_RP. We have built a new spectrophotometric HST/STIS library and used it to derive new sensitivity curves and zero points for the three bands, including recipes on how to correct some cases. Second, we present our results on Galactic runaway stars using Gaia DR1 proper motions: we detect 76 runaway stars, 17 (possibly 19) of them not previously identified as such.

LiLiMaRlin is a library of libraries of massive-star high-resolution optical spectra built by collecting data from [a] our spectroscopic surveys (OWN, IACOB. NoMaDS, and CAF\'E-BEANS) and programs and [b] searches in public archives. The current version has 18 077 spectra of 1665 stars obtained with seven different telescopes (HET 9.2 m, NOT 2.56 m, CAHA 2.2 m, MPG/ESO 2.2 m, OHP 1.93 m, Mercator 1.2 m, and Stella 1.2 m). All the spectra have been filtered to eliminate misidentifications and bad-quality ones, uniformly reprocessed, and placed on a common format. We present applications of this library of libraries to the analysis of spectroscopic binaries (OWN and MONOS, see poster by E. Trigueros P\'aez at this meeting) and the study of the interstellar medium (CollDIBs). We discuss our plans for the future.

CONTEXT: The second Gaia data release (DR2) took place on April 2018. DR2 included photometry for more than 1.3 10^9 sources in G, BP, and RP. Even though Gaia DR2 photometry is very precise, there are currently three alternative definitions of the sensitivity curves that show significative differences. AIMS: The aim of this paper is to improve the quality of the input calibration data to produce new compatible definitions of the three bands and to identify the reasons for the discrepancies between previous definitions. METHODS: We have searched the HST archive for STIS spectra with G430L+G750L data obtained with wide apertures and combined them with the CALSPEC library to produce a high quality SED library of 122 stars with a broad range of colors, including three very red stars. We have used it to compute new sensitivity curves for G, BP, and RP using a functional analytical formalism. RESULTS: The new curves are significantly better than the two previous attempts, REV and WEI. For G we confirm the existence of a systematic bias in magnitude and correct a color term present in REV. For BP we confirm the need to define two magnitude ranges with different sensitivity curves and measure the cut between them at G = 10.87 with a significant increase in precision. The new curves also fit the data better than either REV or WEI. For RP we obtain a sensitivity curve that better fits the STIS spectra and we find that the differences with previous attempts reside in a systematic effect between ground-based and HST spectral libraries. Additional evidence from color-color diagrams indicate that the new sensitivity curve is more accurate. Nevertheless, there is still room for improvement in the accuracy of the sensitivity curves because of the current dearth of good-quality red calibrators: adding more to the sample should be a priority before Gaia data release 3 takes place.

CONTEXT. 2MASS is the reference survey in the NIR part of the spectrum given its whole-sky coverage, large dynamic range, and proven calibration uniformity. However, previous studies disagree in the value of the zero points (ZPs) for its three bands JHK at the hundredth of a magnitude level. The disagreement should become more noticeable now that Gaia provides whole-sky optical photometry calibrated below that level. AIMS. We want to establish the value of the 2MASS ZPs based on NICMOS/HST spectrophotometry of the CALSPEC standard stars and test it with the help of Gaia DR2 parallaxes. METHODS. We have computed the synthetic JHK photometry for a sample of stars using the HST CALSPEC spectroscopic standards and compared it with their 2MASS magnitudes to evaluate the ZPs. We have tested our results by analysing a sample of FGK dwarfs with excellent 2MASS photometry and accurate Gaia DR2 parallaxes.} RESULTS. The Vega ZPs for 2MASS J, H, and K are found to be -0.025$\pm$0.005 mag, 0.004$\pm$0.005 mag, and -0.015$\pm$0.005 mag, respectively. The analysis of the FGK sample indicates that the new ZPs are more accurate than previous ones.

Farr and Mandel reanalyse our data, finding initial-mass-function slopes for high mass stars in 30 Doradus that agree with our results. However, their reanalysis appears to underpredict the observed number of massive stars. Their technique results in more precise slopes than in our work, strengthening our conclusion that there is an excess of massive stars above $30\,\mathrm{M}_\odot$ in 30 Doradus.

The 30 Doradus (30 Dor) nebula in the Large Magellanic Cloud (LMC) is the brightest HII region in the Local Group and a prototype starburst similar to those found in high redshift galaxies. It is thus a stepping stone to understand the complex formation processes of stars in starburst regions across the Universe. Here, we have studied the formation history of massive stars in 30 Dor using masses and ages derived for 452 mainly OB stars from the spectroscopic VLT-FLAMES Tarantula Survey (VFTS). We find that stars of all ages and masses are scattered throughout 30 Dor. This is remarkable because it implies that massive stars either moved large distances or formed independently over the whole field of view in relative isolation. We find that both channels contribute to the 30 Dor massive star population. Massive star formation rapidly accelerated about 8 Myr ago, first forming stars in the field before giving birth to the stellar populations in NGC 2060 and NGC 2070. The R136 star cluster in NGC 2070 formed last and, since then, about 1 Myr ago, star formation seems to be diminished with some continuing in the surroundings of R136. Massive stars within a projected distance of 8 pc of R136 are not coeval but show an age range of up to 6 Myr. Our mass distributions are well populated up to $200\,\mathrm{M}_\odot$. The inferred IMF is shallower than a Salpeter-like IMF and appears to be the same across 30 Dor. By comparing our sample of stars to stellar models in the Hertzsprung-Russell diagram, we find evidence for missing physics in the models above $\log L/\mathrm{L}_\odot=6$ that is likely connected to enhanced wind mass loss for stars approaching the Eddington limit. [abridged]

The IACOB and OWN surveys are two ambitious complementary observational projects which have made available a large multi-epoch spectroscopic database of optical high resolution spectra of Galactic massive O-type stars. As a first step in the study of the full sample of (more than 350) O stars surveyed by the IACOB/OWN projects, we have performed the quantitative spectroscopic analysis of a subsample of 128 stars included in the modern grid of O-type standards for spectral classification. We use semi-automatized tools to determine the set of spectroscopic parameters that can be obtained from the optical spectrum of O-type stars. We also benefit from the multi-epoch character of the surveys to perform a spectroscopic variability study of the sample, accounting for spectroscopic binarity and variability of the main wind diagnostic lines. We provide a general overview of the stellar and wind parameters of this reference sample, and updated recipes for the SpT\,--\,Teff/log g calibrations for Galactic O-type stars. We evaluate our semi-automatized analysis strategy with $\sim$40 stars from the literature, and find a good agreement. The agreement between the synthetic spectra associated with fastwind best fitting models and the observed spectra is good for most targets, but 46 stars present a particular behavior of the wind diagnostic lines that cannot be reproduced by our grid of spherically symmetric unclumped models. These are potential targets of interest for more detailed investigations of clumpy winds and/or the existence of additional circumstellar components. Last, our variability study has led to the detection of signatures of spectroscopic binarity in 27\% of the stars and small amplitude radial velocity variations in the photospheric lines of another 30\%. Additionally, 31\% of the investigated stars show variability in the wind diagnostic lines.

CONTEXT.The first Gaia Data Release (DR1) significantly improved the previously available proper motions for the majority of the Tycho-2 stars. AIMS. We want to detect runaway stars using Gaia DR1 proper motions and compare our results with previous searches. METHODS. Runaway O stars and BA supergiants are detected using a 2-D proper-motion method. The sample is selected using Simbad, spectra from our GOSSS project, literature spectral types, and photometry processed using CHORIZOS. RESULTS. We detect 76 runaway stars, 17 (possibly 19) of them with no prior identification as such, with an estimated detection rate of approximately one half of the real runaway fraction. An age effect appears to be present, with objects of spectral subtype B1 and later having travelled for longer distances than runaways of earlier subtypes. We also tentatively propose that the fraction of runaways is lower among BA supergiants that among O stars but further studies using future Gaia data releases are needed to confirm this. The frequency of fast rotators is high among runaway O stars, which indicates that a significant fraction of them (and possibly a majority) is produced in supernova explosions.

J-PLUS is an ongoing 12-band photometric optical survey, observing thousands of square degrees of the Northern hemisphere from the dedicated JAST/T80 telescope at the Observatorio Astrof\'isico de Javalambre. T80Cam is a 2 sq.deg field-of-view camera mounted on this 83cm-diameter telescope, and is equipped with a unique system of filters spanning the entire optical range. This filter system is a combination of broad, medium and narrow-band filters, optimally designed to extract the rest-frame spectral features (the 3700-4000\AA\ Balmer break region, H$\delta$, Ca H+K, the G-band, the Mgb and Ca triplets) that are key to both characterize stellar types and to deliver a low-resolution photo-spectrum for each pixel of the sky observed. With a typical depth of AB $\sim 21.25$ mag per band, this filter set thus allows for an indiscriminate and accurate characterization of the stellar population in our Galaxy, it provides an unprecedented 2D photo-spectral information for all resolved galaxies in the local universe, as well as accurate photo-z estimates ($\Delta\,z\sim 0.01-0.03$) for moderately bright (up to $r\sim 20$ mag) extragalactic sources. While some narrow band filters are designed for the study of particular emission features ([OII]/$\lambda$3727, H$\alpha$/$\lambda$6563) up to $z < 0.015$, they also provide well-defined windows for the analysis of other emission lines at higher redshifts. As a result, J-PLUS has the potential to contribute to a wide range of fields in Astrophysics, both in the nearby universe (Milky Way, 2D IFU-like studies, stellar populations of nearby and moderate redshift galaxies, clusters of galaxies) and at high redshifts (ELGs at $z\approx 0.77, 2.2$ and $4.4$, QSOs, etc). With this paper, we release $\sim 36$ sq.deg of J-PLUS data, containing about $1.5\times 10^5$ stars and $10^5$ galaxies at $r<21$ mag.

CONTEXT: Many massive stars have nearby companions whose presence hamper their characterization through spectroscopy. AIMS: We want to obtain spatially resolved spectroscopy of close massive visual binaries to derive their spectral types. METHODS: We obtain a large number of short long-slit spectroscopic exposures of five close binaries under good seeing conditions, select those with the best characteristics, extract the spectra using multiple-profile fitting, and combine the results to derive spatially separated spectra. RESULTS: We demonstrate the usefulness of Lucky Spectroscopy by presenting the spatially resolved spectra of the components of each system, in two cases with separations of only ~0.3". Those are delta Ori Aa+Ab (resolved in the optical for the first time) and sigma Ori AaAb+B (first time ever resolved). We also spatially resolve 15 Mon AaAb+B, zeta Ori AaAb+B (both previously resolved with GOSSS, the Galactic O-Star Spectroscopic Survey), and eta Ori AaAb+B, a system with two spectroscopic B+B binaries and a fifth visual component. The systems have in common that they are composed of an inner pair of slow rotators orbited by one or more fast rotators, a characteristic that could have consequences for the theories of massive star formation.

The Tarantula region in the Large Magellanic Cloud contains the richest population of spatially resolved massive O-type stars known so far. This unmatched sample offers an opportunity to test models describing their main-sequence evolution and mass-loss properties. Using ground-based optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS), we aim to determine stellar, photospheric and wind properties of 72 presumably single O-type giants, bright giants and supergiants and to confront them with predictions of stellar evolution and of line-driven mass-loss theories. We apply an automated method for quantitative spectroscopic analysis of O stars combining the non-LTE stellar atmosphere model {{\sc fastwind}} with the genetic fitting algorithm {{\sc pikaia}} to determine the following stellar properties: effective temperature, surface gravity, mass-loss rate, helium abundance, and projected rotational velocity. We present empirical effective temperature versus spectral subtype calibrations at LMC-metallicity for giants and supergiants. In the spectroscopic and classical Hertzsprung-Russell diagrams, our sample O stars are found to occupy the region predicted to be the core hydrogen-burning phase by Brott et al. (2011) and K\"{o}hler et al. (2015). Except for five stars, the helium abundance of our sample stars is in agreement with the initial LMC composition. The aforementioned five stars present moderate projected rotational velocities (i.e., $v_{\mathrm{e}}\,\sin\,i\,<\,200\,\mathrm{km\,s^{-1}}$) and hence do not agree with current predictions of rotational mixing in main-sequence stars. Adopting theoretical results for the wind velocity law, we find modified wind momenta for LMC stars that are $\sim$0.3 dex higher than earlier results. [Due to the limitation of characters, the abstract appearing here is slightly shorter than that in the PDF file.]

The 30 Doradus star-forming region in the Large Magellanic Cloud is a nearby analogue of large star-formation events in the distant Universe. We determine the recent formation history and the initial mass function (IMF) of massive stars in 30 Doradus based on spectroscopic observations of 247 stars more massive than 15 solar masses ($\mathrm{M}_\odot$). The main episode of massive star formation started about $8\,\mathrm{Myr}$ ago and the star-formation rate seems to have declined in the last $1\,\mathrm{Myr}$. The IMF is densely sampled up to $200\,\mathrm{M}_\odot$ and contains $32\pm12\%$ more stars above $30\,\mathrm{M}_\odot$ than predicted by a standard Salpeter IMF. In the mass range $15-200\,\mathrm{M}_\odot$, the IMF power-law exponent is $1.90^{+0.37}_{-0.26}$, shallower than the Salpeter value of 2.35.

[ABRIDGED] Context. O stars are excellent tracers of the intervening ISM because of their high luminosity, blue intrinsic SED, and relatively featureless spectra. We are currently conducting GOSSS, which is generating a large sample of O stars with accurate spectral types within several kpc of the Sun. Aims. To obtain a global picture of the properties of dust extinction in the solar neighborhood based on optical-NIR photometry of O stars with accurate spectral types. Methods. We have processed a photometric set with the CHORIZOS code to measure the amount and type of extinction towards 562 O-type stellar systems. We have tested three different families of extinction laws and analyzed our results with the help of additional archival data. Results. The Ma\'iz Apell\'aniz et al. (2014) family of extinction laws provides a better description of Galactic dust that either the Cardelli et al. (1989) or Fitzpatrick (1999) families, so it should be preferentially used. In many cases O stars and late-type stars experience similar amounts of extinction at similar distances but some O stars are located close to the molecular clouds left over from their births and have larger extinctions than the average for nearby late-type populations. In qualitative terms, O stars experience a more diverse extinction than late-type stars, as some are affected by the small-grain-size, low-R_5495 effect of molecular clouds and others by the large-grain-size, high-R_5495 effect of H II regions. Late-type stars experience a narrower range of grain sizes or R_5495, as their extinction is predominantly caused by the average, diffuse ISM. We propose that the reason for the existence of large-grain-size, high-R_5495 regions in the ISM in the form of H II regions and hot-gas bubbles is the selective destruction of small dust grains by EUV photons and possibly by thermal sputtering by atoms or ions.

Context. On September 2016 the first data from Gaia were released (DR1). The first release included photometry for over 10^9 sources in the very broad G system. Aims. To test the correspondence between G magnitudes in DR1 and the synthetic equivalents derived using spectral energy distributions from observed and model spectrophotometry. To correct the G passband curve and to measure the zero point in the Vega system. Methods. I have computed the synthetic G and Tycho-2 BV photometry for a sample of stars using the Next Generation Spectral Library (NGSL) and the Hubble Space Telescope (HST) CALSPEC spectroscopic standards. Results. I have found that the nominal G passband curve is too blue for the DR1 photometry, as shown by the presence of a color term in the comparison between observed and synthetic magnitudes. A correction to the passband applying a power law in lambda with an exponent of 0.783 eliminates the color term. The corrected passband has a Vega zero point of 0.070$\pm$0.004 magnitudes.

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

The VLT-FLAMES Tarantula Survey has observed hundreds of O-type stars in the 30 Doradus region of the Large Magellanic Cloud (LMC). We study the properties of 105 apparently single O-type dwarfs. To determine stellar and wind parameters, we used the IACOB-GBAT package, an automatic procedure based on a large grid of atmospheric models calculated with the FASTWIND code. In addition to classical techniques, we applied the Bayesian BONNSAI tool to estimate evolutionary masses. We provide a new calibration of effective temperature vs. spectral type for O-type dwarfs in the LMC, based on our homogeneous analysis of the largest sample of such objects to date and including all spectral subtypes. Good agreement with previous results is found, although the sampling at the earliest subtypes could be improved. Rotation rates and helium abundances are studied in an evolutionary context. We find that most of the rapid rotators (vsini higher than 300 km/s ) in our sample have masses below 25 MSun and intermediate rotation-corrected gravities (log gc between 3.9 and 4.1). Such rapid rotators are scarce at higher gravities (i.e. younger ages) and absent at lower gravities (larger ages). This is not expected from theoretical evolutionary models, and does not appear to be due to a selection bias in our sample. We compare the estimated evolutionary and spectroscopic masses, finding a trend that the former is higher for masses below 20 MSun. This can be explained as a consequence of limiting our sample to the O-type stars, and we see no compelling evidence for a systematic mass discrepancy. For most of the stars in the sample we were unable to estimate the wind-strength parameter (hence mass-loss rates) reliably, particularly for objects with luminosity lower than logL/LSun about 5.1. Ultraviolet spectroscopy is needed to undertake a detailed investigation of the wind properties of these dwarfs.

We present the first extensive spectroscopic study of the global population in star clusters Trumpler~16, Trumpler~14 and Collinder~232 in the Carina Nebula, using data from the Gaia-ESO Survey, down to solar-mass stars. In addition to the standard homogeneous Survey data reduction, a special processing was applied here because of the bright nebulosity surrounding Carina stars. We find about four hundred good candidate members ranging from OB types down to slightly sub-solar masses. About one-hundred heavily-reddened early-type Carina members found here were previously unrecognized or poorly classified, including two candidate O stars and several candidate Herbig Ae/Be stars. Their large brightness makes them useful tracers of the obscured Carina population. The spectroscopically-derived temperatures for nearly 300 low-mass members allows the inference of individual extinction values, and the study of the relative placement of stars along the line of sight. We find a complex spatial structure, with definite clustering of low-mass members around the most massive stars, and spatially-variable extinction. By combining the new data with existing X-ray data we obtain a more complete picture of the three-dimensional spatial structure of the Carina clusters, and of their connection to bright and dark nebulosity, and UV sources. The identification of tens of background giants enables us also to determine the total optical depth of the Carina nebula along many sightlines. We are also able to put constraints on the star-formation history of the region, with Trumpler~14 stars found to be systematically younger than stars in other sub-clusters. We find a large percentage of fast-rotating stars among Carina solar-mass members, which provide new constraints on the rotational evolution of pre-main-sequence stars in this mass range.

We have detected 13 new runaway-star candidates of spectral type O combining the TGAS (Tycho-Gaia Astrometric Solution) proper motions from Gaia Data Release 1 (DR1) and the sample from GOSSS (Galactic O-Star Spectroscopic Survey). We have also combined TGAS and Hipparcos proper motions to check that our technique recovers many of the previously known O-type runaways in the sample.

We present the science cases and technological discussions that came from the workshop entitled "Finding the UV-Visible Path Forward" held at NASA GSFC June 25-26, 2015. The material presented outlines the compelling science that can be enabled by a next generation space-based observatory dedicated for UV-visible science, the technologies that are available to include in that observatory design, and the range of possible alternative launch approaches that could also enable some of the science. The recommendations to the Cosmic Origins Program Analysis Group from the workshop attendees on possible future development directions are outlined.

We have used (a) HST ACS imaging and STIS spectroscopy, (b) ground-based PIONIER/VLT long-baseline interferometry, and (c) ground-based spectroscopy from different instruments to study the orbit of the extreme multiple system HD 93 129 Aa,Ab, which is composed of (at least) two very massive stars in a long-period orbit with e>0.92 that will pass through periastron in 2017/2018. In several ways, the system is an eta Car precursor. Around the time of periastron passage the two very strong winds will collide and generate an outburst of non-thermal hard X-ray emission without precedent in an O+O binary since astronomers have been able to observe above Earth's atmosphere. A coordinated multiwavelength monitoring in the next two years will enable a breakthrough understanding of the wind interactions in such extreme close encounters. Furthermore, we have found evidence that HD 93 129 Aa may be a binary system itself. In that case, we could witness a three-body interaction that may yield a runaway star or a stellar collision close to or shortly after the periastron passage. Either of those outcomes would be unprecedented, as they are predicted to be low-frequency events in the Milky Way.

We present the updates of the Galactic O-Star Catalog (GOSC) that we have undertaken in the last two years: new spectral types, more objects, additional information, and coordination with CDS. We also present updates for the Galactic O-Star Spectroscopic Survey (GOSSS). A new paper (GOSSS-III) has been published and ~1000 targets have been observed since 2014. Four new setups have been added to our lineup and for two of them we have already obtained over 100 spectra: with OSIRIS at the 10.4 m GTC we are observing northern dim stars and with FRODOspec at the 2.0 m Liverpool Telescope we are observing northern bright stars. Finally, we also make available new versions of MGB, the spectral classification tool associated with the project, and of the GOSSS grid of spectroscopic standards.

We present our progress on the study of extinction laws along three diferent lines. [a] We compare how well different families of extinction laws fit existing photometric data for Galactic sightlines and we find that the Ma\'iz Apell\'aniz et al. (2014) family provides better results than those of Cardelli et al. (1989) or Fitzpatrick (1999). [b] We describe the HST/STIS spectrophotometry in the 1700-10 200 Angstrom range that we are obtaining for several tens of sightlines in 30 Doradus with the purpose of deriving an improved wavelength-detailed family of extinction laws. [c] We present the study we are conducting on the behavior of the extinction law in the infrared by combining 2MASS and WISE photometry with Spitzer and ISO spectrophotometry.

This is the third installment of GOSSS, a massive spectroscopic survey of Galactic O stars, based on new homogeneous, high signal-to-noise ratio, R~2500 digital observations selected from the Galactic O-Star Catalog (GOSC). In this paper we present 142 additional stellar systems with O stars from both hemispheres, bringing the total of O-type systems published within the project to 590. Among the new objects there are 20 new O stars. We also identify 11 new double-lined spectroscopic binaries (SB2s), of which 6 are of O+O type and 5 of O+B type, and an additional new tripled-lined spectroscopic binary (SB3) of O+O+B type. We also revise some of the previous GOSSS classifications, present some egregious examples of stars erroneously classified as O-type in the past, introduce the use of luminosity class IV at spectral types O4-O5.5, and adapt the classification scheme to the work of Arias et al. (2016).

We study the O-type star HD 161853, which has been noted as a probable double-lined spectroscopic binary system. We secured high-resolution spectra of HD 161853 during the past nine years. We separated the two components in the system and measured their respective radial velocities for the first time. We confirm that HD 161853 is an $\sim$1 Ma old binary system consisting of an O8 V star ($M_{\rm A,RV} \geq 22$ M$_\odot$) and a B1--3 V star ($M_{\rm B,RV} \geq 7.2$ M$_\odot$) at about 1.3 kpc. From the radial velocity curve, we measure an orbital period $P$ = 2.66765$\pm$0.00001 d and an eccentricity $e$ = 0.121$\pm$0.007. Its $V$-band light curve is constant within 0.014 mag and does not display eclipses, from which we impose a maximum orbital inclination $i=54$ deg. HD 161853 is probably associated with an H II region and a poorly investigated very young open cluster. In addition, we detect a compact emission region at 50 arcsec to HD 161853 in 22$\mu$m-WISE and 24$\mu$m-Spitzer images, which may be identified as a dust wave piled up by the radiation pressure of the massive binary system.

Aims. We present the first orbital solution for the O-type supergiant star HD 74194, which is the optical counterpart of the supergiant fast X-ray transient IGR J08408-4503. Methods. We measured the radial velocities in the optical spectrum of HD 74194, and we determined the orbital solution for the first time. We also analysed the complex H{\alpha} profile. Results. HD 74194 is a binary system composed of an O-type supergiant and a compact object in a short-period ($P=9.5436\pm0.0002$ d) and high-eccentricity ($e=0.63\pm0.03$) orbit. The equivalent width of the H{\alpha} line is not modulated entirely with the orbital period, but seems to vary in a superorbital period ($P=285\pm10$ d) nearly 30 times longer than the orbital one.

The contact phase expected to precede the coalescence of two massive stars is poorly characterized due to the paucity of observational constraints. Here we report on the discovery of VFTS 352, an O-type binary in the 30 Doradus region, as the most massive and earliest spectral type overcontact system known to date. We derived the 3D geometry of the system, its orbital period $P_{\rm orb}=1.1241452(4)$ d, components' effective temperatures -- $T_1=42\,540\pm280$ K and $T_2=41\,120\pm290$ K -- and dynamical masses -- $M_1=28.63\pm0.30 M_{\odot}$ and $M_2 = 28.85\pm0.30 M_{\odot}$. Compared to single-star evolutionary models, the VFTS 352 components are too hot for their dynamical masses by about 2700 and 1100 K, respectively. These results can be explained naturally as a result of enhanced mixing, theoretically predicted to occur in very short-period tidally-locked systems. The VFTS 352 components are two of the best candidates identified so far to undergo this so-called chemically homogeneous evolution. The future of VFTS 352 is uncertain. If the two stars merge, a very rapidly rotating star will be produced. Instead, if the stars continue to evolve homogeneously and keep shrinking within their Roche Lobes, coalescence can be avoided. In this case, tides may counteract the spin down by winds such that the VFTS 352 components may, at the end of their life, fulfill the requirements for long gamma-ray burst (GRB) progenitors in the collapsar scenario. Independently of whether the VFTS 352 components become GRB progenitors, this scenario makes VFTS 352 interesting as a progenitor of a black hole binary, hence as a potential gravitational wave source through black hole-black hole merger.

Context: Nearly one century after their discovery, the carrier(s) of Diffuse Interstellar Bands is/are still unknown and there are few sightlines studied in detail for a large number of DIBs. Aims: We want to study the ISM sightlines towards LS III +46 11 and LS III +46 12, two early-O-type stellar systems, and LS III +46 11 B, a mid-B-type star. The three targets are located in the stellar cluster Berkeley 90 and have a high extinction. Methods: We use the multi-epoch high-S/N optical spectra presented in paper I (Ma\'iz Apell\'aniz et al. 2015), the extinction results derived there, and additional spectra. Results: We have measured equivalent widths, velocities, and FWHMs for a large number of absorption lines in the rich ISM spectrum in front of Berkeley 90. The absorbing ISM has at least two clouds at different velocities, one with a lower column density (thinner) in the K I lines located away from Berkeley 90 and another one with a higher column density (thicker) associated with the cluster. The first cloud has similar properties for both O-star sightlines but the second one is thicker for LS III +46 11. The comparison between species indicate that the cloud with a higher column density has a denser core, allowing us to classify the DIBs in a sigma-zeta scale, some of them for the first time. The LS III +46 12 sightline also has a high-velocity redshifted component.

I present results on the interstellar medium towards the O stars observed in five optical spectroscopic surveys: GOSSS, OWN, IACOB, NoMaDS, and CAF\'E-BEANS. I have measured both the amount [$E(4405-5495)$] and type [$R_{5495}$] of extinction towards several hundreds of Galactic O stars and verified that the Ma\'iz Apell\'aniz et al. (2014) family of extinction laws provides a significantly better fit to optical+NIR Galactic extinction than either the Cardelli et al. (1989) or the Fitzpatrick (1999) families. $R_{5495}$ values are concentrated between 3.0 and 3.5 but for low values of $E(4405-5495)$ there is a significant population with larger $R_{5495}$ associated with H II regions. I have also measured different DIBs and I have found that $W$(5797)/$W$(5780) is anticorrelated with $R_{5495}$, a sign that extreme $\zeta$ clouds are characterized not only by low ionization environments (as opposed to $\sigma$ clouds) but also by having a larger fraction of small dust grains. The equivalent width of the "Gaia DIB" (8621 \AA) is strongly correlated with $E(4405-5495)$, as expected, and its behavior appears to be more $\sigma$-like than $\zeta$-like. We have also started analyzing some individual sightlines in detail.

Extremely broad emission wings at H$\beta$ and H$\alpha$ have been found in VFTS data for five very luminous BA supergiants in or near 30 Doradus in the Large Magellanic Cloud. The profiles of both lines are extremely asymmetrical, which we have found to be caused by very broad diffuse interstellar bands (DIBs) in the longward wing of H$\beta$ and the shortward wing of H$\alpha$. These DIBs are well known to interstellar but not to many stellar specialists, so that the asymmetries may be mistaken for intrinsic features. The broad emission wings are generally ascribed to electron scattering, although we note difficulties for that interpretation in some objects. Such profiles are known in some Galactic hyper/supergiants and are also seen in both active and quiescent Luminous Blue Variables. No prior or current LBV activity is known in these 30 Dor stars, although a generic relationship to LBVs is not excluded; subject to further observational and theoretical investigation, it is possible that these very luminous supergiants are approaching the LBV stage for the first time. Their locations in the HRD and presumed evolutionary tracks are consistent with that possibility. The available evidence for spectroscopic variations of these objects is reviewed, while recent photometric monitoring does not reveal variability. A search for circumstellar nebulae has been conducted, with an indeterminate result for one of them.

We investigate the multiplicity properties of 408 B-type stars observed in the 30 Doradus region of the Large Magellanic Cloud with multi-epoch spectroscopy from the VLT-FLAMES Tarantula Survey (VFTS). We use a cross-correlation method to estimate relative radial velocities from the helium and metal absorption lines for each of our targets. Objects with significant radial-velocity variations (and with an amplitude larger than 16 km/s) are classified as spectroscopic binaries. We find an observed spectroscopic binary fraction (defined by periods of <10^3.5 d and mass ratios >0.1) for the B-type stars, f_B(obs) = 0.25 +/- 0.02, which appears constant across the field of view, except for the two older clusters (Hodge 301 and SL 639). These two clusters have significantly lower fractions of 0.08 +/- 0.08 and 0.10 +/- 0.09, respectively. Using synthetic populations and a model of our observed epochs and their potential biases, we constrain the intrinsic multiplicity properties of the dwarf and giant (i.e. relatively unevolved) B-type stars in 30 Dor. We obtain a present-day binary fraction f_B(true) = 0.58 +/- 0.11, with a flat period distribution. Within the uncertainties, the multiplicity properties of the B-type stars agree with those for the O stars in 30 Dor from the VFTS.

In this work we present high resolution spectroscopic data of the giant star-forming region of N11, obtained with the GIRAFFE instrument at the Very Large Telescope. By using this data set, we find that most of the H$\alpha$ emission lines profiles in this complex can be fitted by a single Gaussian, however, multiple emission line profiles can be observed in the central region of N11. By adding all the spectra, we derive the integrated H$\alpha$ profile of this complex, which displays a width ($\sigma$) of about 12 km s$^{-1}$ (corrected by instrumental and thermal width). We find that a single Gaussian fit on the integrated H$\alpha$ profile leaves remaining wings, which can be fitted by a secondary broad Gaussian component. In addition, we find high velocity features, which spatially correlate with soft diffuse X-ray emission.

Context: It appears that most (if not all) massive stars are born in multiple systems. At the same time, the most massive binaries are hard to find due to their low numbers throughout the Galaxy and the implied large distances and extinctions. AIMS: We want to study: [a] LS III +46 11, identified in this paper as a very massive binary; [b] another nearby massive system, LS III +46 12; and [c] the surrounding stellar cluster, Berkeley 90. Methods: Most of the data used in this paper are multi-epoch high-S/N optical spectra though we also use Lucky Imaging and archival photometry. The spectra are reduced with devoted pipelines and processed with our own software, such as a spectroscopic-orbit code, CHORIZOS, and MGB. Results: LS III +46 11 is identified as a new very-early-O-type spectroscopic binary [O3.5 If* + O3.5 If*] and LS III +46 12 as another early O-type system [O4.5 V((f))]. We measure a 97.2-day period for LS III +46 12 and derive minimum masses of 38.80$\pm$0.83 M_Sol and 35.60$\pm$0.77 M_Sol for its two stars. We measure the extinction to both stars, estimate the distance, search for optical companions, and study the surrounding cluster. In doing so, a variable extinction is found as well as discrepant results for the distance. We discuss possible explanations and suggest that LS III +46 12 may be a hidden binary system, where the companion is currently undetected.

The IACOB project is an ambitious long-term project which is contributing to step forward in our knowledge about the physical properties and evolution of Galactic massive stars. The project aims at building a large database of high-resolution, multi-epoch, spectra of Galactic OB stars, and the scientific exploitation of the database using state-of-the-art models and techniques. In this proceeding, we summarize the latest updates of the IACOB spectroscopic database and highlight some of the first scientific results from the IACOB project; we also announce the first data release and the first public version of the iacob-broad tool for the line-broadening characterization of OB-type spectra.

We report magnetic and spectroscopic observations and modeling of the Of?p star CPD -28 2561. Using more than 75 new spectra, we have measured the equivalent width variations and examined the dynamic spectra of photospheric and wind-sensitive spectral lines. A period search results in an unambiguous 73.41 d variability period. High resolution spectropolarimetric data analyzed using Least-Squares Deconvolution yield a Zeeman signature detected in the mean Stokes V profile corresponding to phase 0.5 of the spectral ephemeris. Interpreting the 73.41 d period as the stellar rotational period, we have phased the equivalent widths and inferred longitudinal field measurements. The phased magnetic data exhibit a weak sinusoidal variation, with maximum of about 565 G at phase 0.5, and a minimum of about -335 G at phase 0.0, with extrema approximately in phase with the (double-wave) Halpha equivalent width variation. Modeling of the Halpha equivalent width variation assuming a quasi-3D magnetospheric model produces a unique solution for the ambiguous couplet of inclination and magnetic obliquity angles: (i, beta) or (beta, i)=(35 deg,90 deg). Adopting either geometry, the longitudinal field variation yields a dipole polar intensity Bd=2.6\pm 0.9~kG, consistent with that obtained from direct modelling of the Stokes V profiles. We derive a wind magnetic confinement parameter eta*\simeq 100, leading to an Alfv\'en radius RA\simeq 3-5~R*, and a Kepler radius RK\simeq 20~R*. This supports a physical scenario in which the Halpha emission and other line variability have their origin in an oblique, co-rotating 'dynamical magnetosphere' structure resulting from a magnetically channeled wind. Nevertheless, the details of the formation of spectral lines and their variability within this framework remain generally poorly understood.

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.

The stellar population of the 30 Doradus star-forming region in the Large Magellanic Cloud contains a subset of apparently single, rapidly rotating O-type stars. The physical processes leading to the formation of this cohort are currently uncertain. One member of this group, the late O-type star VFTS 399, is found to be unexpectedly X-ray bright for its bolometric luminosity - in this study we aim to determine its physical nature and the cause of this behaviour. We find VFTS 399 to be an aperiodic photometric variable with an apparent near-IR excess. Its optical spectrum demonstrates complex emission profiles in the lower Balmer series and select HeI lines - taken together these suggest an OeBe classification. The highly variable X-ray luminosity is too great to be produced by a single star, while the hard, non-thermal nature suggests the presence of an accreting relativistic companion. Finally, the detection of periodic modulation of the X-ray lightcurve is most naturally explained under the assumption that the accretor is a neutron star. VFTS 399 appears to be the first high-mass X-ray binary identified within 30 Dor. Comparison of the current properties of VFTS 399 to binary-evolution models suggests a progenitor mass in excess of 25Msun for the putative neutron star, which may host a magnetic field comparable in strength to those of magnetars. VFTS 399 is now the second member of the cohort of rapidly rotating `single' O-type stars in 30 Dor to show evidence of binary interaction resulting in spin-up, suggesting that this may be a viable evolutionary pathway for the formation of a subset of this stellar population.

We present spectral classifications for 438 B-type stars observed as part of the VLT-FLAMES Tarantula Survey (VFTS) in the 30 Doradus region of the Large Magellanic Cloud. Radial velocities are provided for 307 apparently single stars, and for 99 targets with radial-velocity variations which are consistent with them being spectroscopic binaries. We investigate the spatial distribution of the radial velocities across the 30 Dor region, and use the results to identify candidate runaway stars. Excluding potential runaways and members of two older clusters in the survey region (SL 639 and Hodge 301), we determine a systemic velocity for 30 Dor of 271.6 +/- 12.2 km/s from 273 presumed single stars. Employing a 3-sigma criterion we identify nine candidate runaway stars (2.9% of the single stars with radial-velocity estimates). The projected rotational velocities of the candidate runaways appear to be significantly different to those of the full B-type sample, with a strong preference for either large (>345 km/s) or small (<65 km/s) rotational velocities. Of the candidate runaways, VFTS 358 (classified B0.5: V) has the largest differential radial velocity (-106.9 +/- 16.2 km/s), and a preliminary atmospheric analysis finds a significantly enriched nitrogen abundance of 12+log(N/H) > ~8.5. Combined with a large rotational velocity (vsini = 345 +/- 22 km/s), this is suggestive of past binary interaction for this star.

TLUSTY non-LTE model atmosphere calculations have been used to determine atmospheric parameters and nitrogen (N) abundances for 34 single and 18 binary B-type supergiants (BSGs). The effects of flux contribution from an unseen secondary were considered for the binary sample. We present the first systematic study of the incidence of binarity for a sample of BSGs across the theoretical terminal age main sequence (TAMS). To account for the distribution of effective temperatures of the BSGs it may be necessary to extend the TAMS to lower temperatures. This is consistent with the derived distribution of mass discrepancies, projected rotational velocities (vsini) and N abundances, provided that stars cooler than this temperature are post RSG objects. For the BSGs in the Tarantula and previous FLAMES surveys, most have small vsini. About 10% have larger vsini (>100 km/s) but surprisingly these show little or no N enhancement. All the cooler BSGs have low vsini of <70km/s and high N abundance estimates, implying that either bi-stability braking or evolution on a blue loop may be important. A lack of cool binaries, possibly reflects the small sample size. Single star evolutionary models, which include rotation, can account for the N enhancement in both the single and binary samples. The detailed distribution of N abundances in the single and binary samples may be different, possibly reflecting differences in their evolutionary history. The first comparative study of single and binary BSGs has revealed that the main sequence may be significantly wider than previously assumed, extending to Teff=20000K. Some marginal differences in single and binary atmospheric parameters and abundances have been identified, possibly implying non-standard evolution for some of the sample. This sample as a whole has implications for several aspects of our understanding of the evolution of BSGs. Full abstract in paper

CAF\'E-BEANS is an on-going survey running on the 2.2 m telescope at Calar Alto. For more than two years, CAF\'E-BEANS has been collecting high-resolution spectra of early-type stars with the aim of detecting and characterising spectroscopic binaries. The main goal of this project is a thorough characterisation of multiplicity in high-mass stars by detecting all spectroscopic and visual binaries in a large sample of Galactic O-type stars, and solving their orbits. Our final objective is eliminating all biases in the high-mass-star IMF created by undetected binaries.

We provide a complete characterization of the astrophysical properties of the $\sigma$ Ori Aa,Ab,B hierarchical triple system, and an improved set of orbital parameters for the highly eccentric $\sigma$ Ori Aa,Ab spectroscopic binary. We compiled a spectroscopic dataset comprising 90 high-resolution spectra covering a total time span of 1963 days. We applied the Lehman-Filh\'es method for a detailed orbital analysis of the radial velocity curves and performed a combined quantitative spectroscopic analysis of the {$\sigma$ Ori Aa,Ab,B} system by means of the stellar atmosphere code FASTWIND. We used our own plus other available information on photometry and distance to the system for measuring the radii, luminosities, and spectroscopic masses of the three components. We also inferred evolutionary masses and stellar ages using the Bayesian code BONNSAI. The orbital analysis of the new radial velocity curves led to a very accurate orbital solution of the $\sigma$ Ori Aa,Ab pair. We provided indirect arguments indicating that $\sigma$ Ori B is a fast rotating early-B dwarf. The FASTWIND+BONNSAI analysis showed that the Aa,Ab pair contains the hottest and most massive components of the triple system while $\sigma$ Ori B is a bit cooler and less massive. The derived stellar ages of the inner pair are intriguingly younger than the one widely accepted for the $\sigma$ Orionis cluster, at 3$\pm$1 Ma. The outcome of this study will be of key importance for a precise determination of the distance to the $\sigma$ Orionis cluster, the interpretation of the strong X-ray emission detected for $\sigma$ Ori Aa,Ab,B, and the investigation of the formation and evolution of multiple massive stellar systems and substellar objects.

In this poster we present three developments related to the Galactic O-Star Spectroscopic Survey (GOSSS). First, we are making public the first version of MGB, an IDL code that allows the user to compare oberved spectra to a grid of spectroscopic standards to measure spectral types, luminosity classes, rotation indexes, and spectral qualifiers. Second, we present the associated grid of standard stars for the spectral types O2 to O9.7, with several improvements over the original GOSSS grid of Sota et al. (2011). Third, we present a list of egregious classification errors in SIMBAD: stars that are or have been listed there as being of O type but that in reality are late-type stars.

We have conducted a 1.5 year-long variability study of the stars in the Cygnus OB2 association, the region in the northern hemisphere with the highest density of optically visible massive stars. The survey was conducted using four pointings in the Johnson $R$ and $I$ bands with a 35 cm Meade LX200-ACF telescope equipped with a 3.2 Mpixel SBIG ST10-XME CCD camera and includes 300+ epochs in each filter. A total of 1425 objects were observed with limiting magnitudes of 15 in $R$ and 14 in $I$. The photometry was calibrated using reference stars with existing $UBVJHK$ photometry. Bright stars have precisions better than 0.01 magnitudes, allowing us to detect 52 confirmed and 19 candidate variables, many of them massive stars without previous detections as variables. Variables are classified as eclipsing, pulsating, irregular/long period, and Be. We derive the phased light curves for the eclipsing binaries, with periods ranging from 1.3 to 8.5 days.

We present a multi-filter HST/WFPC2 UV-optical study of the Hourglass region in M8. We have extracted the stellar photometry of the sources in the area and obtained the separations and position angles of the Herschel 36 multiple system: for Herschel 36 D we detect a possible orbital motion between 1995 and 2009. We have combined our data with archival IUE spectroscopy and measured the Herschel 36 extinction law, obtaining a different result from that of Cardelli et al. (1989) due to the improvement in the quality of the optical-NIR data, in agreement with the results of Ma\'iz Apell\'aniz et al. (2014). A large fraction of the UV flux around Herschel 36 arises from the Hourglass and not directly from the star itself. In the UV the Hourglass appears to act as a reflection nebula located behind Herschel 36 along the line of sight. Finally, we also detect three new Herbig-Haro objects and the possible anisotropic expansion of the Hourglass Nebula.

We have recently derived a family of extinction laws for 30 Doradus that provides better fits to the optical photometry of obscured stars in the Galaxy and the LMC. Simultaneously, we are extending our Galactic O-Star Spectroscopic Survey (GOSSS) to fainter, more extinguished stars to obtain accurate spectral types for massive stars with more than 6 magnitudes of $V$-band extinction. I have combined both lines of research with 2MASS, WISE, and Spitzer photometry to obtain the 1-10 micron extinction law for O stars in the solar neighborhood. I present these results and compare them with the extinction laws in the same wavelength range derived from late-type stars and H II regions. I also discuss plans to extend the newly derived optical-IR extinction laws to the UV.

Context: Massive stars are extremely important for the evolution of the galaxies; there are large gaps in our understanding of their properties and formation, however, mainly because they evolve rapidly, are rare, and distant. It may well be that almost all massive stars are born as triples or higher multiples, but their large distances require very high angular resolution to directly detect the companions at milliarcsecond scales. Aims: Herschel36 is a young massive system located at 1.3 kpc. It has a combined smallest predicted mass of 45 M_sun. Multi-epoch spectroscopic data suggest the existence of at least three gravitationally bound components. Two of them, system Ab, are tightly bound in a spectroscopic binary, and the third one, component Aa, orbits in a wider orbit. Our aim was to image and obtain astrometric and photometric measurements of components Aa and Ab using, for the first time, long-baseline optical interferometry to further constrain its nature. Methods: We observed Herschel 36 with the near-infrared instrument AMBER attached to the ESO VLT Interferometer, which provides an angular resolution of approx. 2 mas. We used the code BSMEM to perform the interferometric image reconstruction. We fitted the interferometric observables using proprietary IDL routines and the code LitPro. Results: We imaged the Aa+Ab components of Herschel 36 in H and K filters. Component Ab is located at a projected distance of 1.81 mas, at a position angle of approx. 222 deg. east of north, the flux ratio between components Aa and Ab is close to one. The small measured angular separation indicates that system Ab and Ab may be approaching the periastron of their orbits. These results, only achievable with long-baseline near-infrared interferometry, constitute the first step toward a thorough understanding of this massive triple system.

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.

Rotation is a key parameter in the evolution of massive stars, affecting their evolution, chemical yields, ionizing photon budget, and final fate. We determined the projected rotational velocity, $v_e\sin i$, of $\sim$330 O-type objects, i.e. $\sim$210 spectroscopic single stars and $\sim$110 primaries in binary systems, in the Tarantula nebula or 30 Doradus (30\,Dor) region. The observations were taken using VLT/FLAMES and constitute the largest homogeneous dataset of multi-epoch spectroscopy of O-type stars currently available. The most distinctive feature of the $v_e\sin i$ distributions of the presumed-single stars and primaries in 30 Dor is a low-velocity peak at around 100\,$\rm{km s^{-1}}$. Stellar winds are not expected to have spun-down the bulk of the stars significantly since their arrival on the main sequence and therefore the peak in the single star sample is likely to represent the outcome of the formation process. Whereas the spin distribution of presumed-single stars shows a well developed tail of stars rotating more rapidly than 300\,$\rm{km s^{-1}}$, the sample of primaries does not feature such a high-velocity tail. The tail of the presumed-single star distribution is attributed for the most part -- and could potentially be completely due -- to spun-up binary products that appear as single stars or that have merged. This would be consistent with the lack of such post-interaction products in the binary sample, that is expected to be dominated by pre-interaction systems. The peak in this distribution is broader and is shifted toward somewhat higher spin rates compared to the distribution of presumed-single stars. Systems displaying large radial velocity variations, typical for short period systems, appear mostly responsible for these differences.

It has been suggested that the compact open cluster VdBH 222 is a young massive distant object. We set out to characterise VdBH 222 using a comprehensive set of multi-wavelength observations. We obtained multi-band optical (UBVR) and near-infrared (JHK) photometry of the cluster field, as well as multi-object and long-slit optical spectroscopy for a large sample of stars in the field. We applied classical photometric analysis, as well as more sophisticated methods using the CHORIZOS code, to determine the reddening to the cluster. We then plotted dereddened HR diagrams and determined cluster parameters via isochrone fitting. We have identified a large population of luminous supergiants confirmed as cluster members via radial velocity measurements. We find nine red supergiants (plus one other candidate) and two yellow supergiants. We also identify a large population of OB stars. Ten of them are bright enough to be blue supergiants. The cluster lies behind 7.5 mag of extinction for the preferred value of Rv of 2.9. Isochrone fitting allows for a narrow range of ages between 12 and 16 Ma. The cluster radial velocity is compatible with distances of 6 and 10 kpc. The shorter distance is inconsistent with the age range and Galactic structure. The longer distance implies an age of 12 Ma and a location not far from the position where some Galactic models place the far end of the Galactic Bar. VdBH 222 is a young massive cluster with a likely mass greater than 20000 Msolar . Its population of massive evolved stars is comparable to that of large associations, such as Per OB1. Its location in the inner Galaxy, presumably close to the end of the Galactic bar, adds to the increasing evidence for vigorous star formation in the inner regions of the Milky Way.

Context: The commonly used extinction laws of Cardelli et al. (1989) have limitations that, among other issues, hamper the determination of the effective temperatures of O and early B stars from optical+NIR photometry. Aims: We aim to develop a new family of extinction laws for 30 Doradus, check their general applicability within that region and elsewhere, and apply them to test the feasibility of using optical+NIR photometry to determine the effective temperature of OB stars. Methods: We use spectroscopy and NIR photometry from the VLT-FLAMES Tarantula Survey and optical photometry from HST/WFC3 of 30 Doradus and we analyze them with the software code CHORIZOS using different assumptions such as the family of extinction laws. Results: We derive a new family of optical+NIR extinction laws for 30 Doradus and confirm its applicability to extinguished Galactic O-type systems. We conclude that by using the new extinction laws it is possible to measure the effective temperatures of OB stars with moderate uncertainties and only a small bias, at least up to E(4405-5495) ~ 1.5 mag.

We present the second installment of GOSSS, a massive spectroscopic survey of Galactic O stars, based on new homogeneous, high signal-to-noise ratio, R ~ 2500 digital observations from both hemispheres selected from the Galactic O-Star Catalog (GOSC). In this paper we include bright stars and other objects drawn mostly from the first version of GOSC, all of them south of delta = -20 degrees, for a total number of 258 O stars. We also revise the northern sample of paper I to provide the full list of spectroscopically classified Galactic O stars complete to B = 8, bringing the total number of published GOSSS stars to 448. Extensive sequences of exceptional objects are given, including the early Of/WN, O Iafpe, Ofc, ON/OC, Onfp, Of?p, and Oe types, as well as double/triple-lined spectroscopic binaries. The new spectral subtype O9.2 is also discussed. The magnitude and spatial distributions of the observed sample are analyzed. We also present new results from OWN, a multi-epoch high-resolution spectroscopic survey coordinated with GOSSS that is assembling the largest sample of Galactic spectroscopic massive binaries ever attained. The OWN data combined with additional information on spectroscopic and visual binaries from the literature indicate that only a very small fraction (if any) of the stars with masses above 15-20 M_Sol are born as single systems. In the future we will publish the rest of the GOSSS survey, which is expected to include over 1000 Galactic O stars.

AIMS. We test the hypothesis of O Vz stars (characterized by having HeII4686 stronger in absorption than other He lines in their blue-violet spectra) being at a younger evolutionary stage than are normal O-type dwarfs. METHODS. We have performed a quantitative spectroscopic analysis of a sample of 38 O Vz and 46 O V stars, identified by the VLT-FLAMES Tarantula Survey in the 30 Doradus region of the Large Magellanic Cloud (LMC). We obtained the stellar and wind parameters of both samples using the FASTWIND stellar atmosphere code and the IACOB-GBAT grid-based automatic tool. In the framework of a differential study, we compared the physical and evolutionary properties of both samples, regarding Teff, logg, logQ and logL. We also investigated the predictions of the FASTWIND code about the O Vz phenomenon. RESULTS. We find a differential distribution of objects in terms of effective temperature, with O Vz stars dominant at intermediate values. The O Vz stars in 30 Doradus tend to be younger and less luminous, and they have weaker winds than the O V stars, but we also find examples with ages of 2-4 Myr and with luminosities and winds that are similar to those of normal O dwarfs. Moreover, the O Vz stars do not appear to have higher gravities than the O V stars. In addition to effective temperature and wind strength, our FASTWIND predictions indicate how important it is to take other stellar parameters (gravity and projected rotational velocity) into account for correctly interpreting the O Vz phenomenon. CONCLUSIONS. In general, the O Vz stars appear to be on or very close to the ZAMS, but there are some examples where the Vz classification does not necessarily imply extreme youth. In particular, the presence of O Vz stars in our sample at more evolved phases than expected is likely a consequence of modest O-star winds owing to the low-metallicity environment of the LMC.

The 30 Doradus (30\,Dor) region in the Large Magellanic Cloud (also known as the Tarantula Nebula) is the nearest massive starburst region, containing the richest sample of massive stars in the Local Group. It is the best possible laboratory to investigate aspects of the formation and evolution of massive stars. Here, we focus on rotation which is a key parameter in the evolution of these objects. We establish the projected rotational velocity, $v_{e}\sin i$, distribution of an unprecedented sample of 216 radial velocity constant ($\rm{\Delta RV\, \leq\, 20 \,km s^{-1}}$) O-type stars in 30\,Dor observed in the framework of the VLT-FLAMES Tarantula Survey (VFTS). The distribution of $v_{e}\sin i$ shows a two-component structure: a peak around 80 $\rm{km s^{-1}}$ and a high-velocity tail extending up to $\sim$600 $\rm{km s^{-1}}$. Around 75% of the sample has 0 $\leq\, v_{e}\sin i \leq$ 200 $\rm{km s^{-1}}$ with the other 25% distributed in the high-velocity tail. The presence of the low-velocity peak is consistent with that found in other studies of late-O and early-B stars. The high-velocity tail is compatible with expectations from binary interaction synthesis models and may be predominantly populated by post-binary interaction, spun-up, objects and mergers. This may have important implications for the nature of progenitors of long-duration gamma ray bursts.

Aims. Using ground based multi-object optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS), we aim to establish the (projected) rotational velocity distribution for a sample of 216 presumably single O-type stars in 30 Doradus (30 Dor). Methods. We measured projected rotational velocities, \vrot, by means of a Fourier transform method and a profile fitting method applied on a set of isolated spectral lines. We also used an iterative deconvolution procedure to infer the probability density, $\rm{P(\veq)}$, of the equatorial rotational velocity, \veq. Results. The distribution of \vrot\ shows a two-component structure: a peak around 80 \kms\ and a high-velocity tail extending up to $\sim$600 \kms. This structure is also present in the inferred distribution $\rm{P(\veq)}$ with around 80% of the sample having 0 $<$ \veq\, $\leq\, 300$ \kms\ and the other 20% distributed in the high-velocity region. Conclusions. Most of the stars in our sample rotate with a rate less than 20%\ of their break-up velocity. For the bulk of the sample, mass-loss in a stellar wind and/or envelope expansion is not efficient enough to significantly spin down these stars within the first few Myr of evolution. The presence of a sizeable population of fast rotators is compatible with recent population synthesis computations that investigate the influence of binary evolution on the rotation rate of massive stars. Despite the fact that we have excluded stars that show significant radial velocity variations, our sample may have remained contaminated by post-interaction binary products. The fact that the high-velocity tail may be preferentially (and perhaps even exclusively), populated by post-binary interaction products, has important implications for the evolutionary origin of systems that produce gamma-ray bursts.

The Galactic O-Star Spectroscopic Survey (GOSSS) is obtaining high quality R~2500 blue-violet spectroscopy of all Galactic stars ever classified as of O type with B < 12 and a significant fraction of those with B = 12-14. As of June 2013, we have obtained, processed, and classified 2653 spectra of 1593 stars, including all of the sample with B < 8 and most of the sample with B = 8-10, making GOSSS already the largest collection of high quality O-star optical spectra ever assembled by a factor of 3. We discuss the fraction of false positives (stars classified as O in previous works that do not belong to that class) and the implications of the observed magnitude distribution for the spatial distribution of massive stars and dust within a few kpc of the Sun. We also present new spectrograms for some of the interesting objects in the sample and show applications of GOSSS data to the study of the intervening ISM. Finally, we present the new version of the Galactic O-Star Catalog (GOSC), which incorporates the data in GOSSS-DR1, and we discuss our plans for MGB, an interactive spectral classification tool for OB stars.

We are using five different surveys to compile the largest sample of diffuse interstellar band (DIB) measurements ever collected. GOSSS is obtaining intermediate-resolution blue-violet spectroscopy of ~2500 OB stars, of which 60% have already been observed and processed. The other four surveys have already collected multi-epoch high-resolution optical spectroscopy of 700 OB stars with different telescopes, including the 9 m Hobby-Eberly Telescope in McDonald Observatory. Some of our stars are highly-extinguished targets for which no good-quality optical spectra have ever been published. For all of the targets in our sample we have obtained accurate spectral types, measured non-DIB ISM lines, and compiled information from the literature to calculate the extinction. Here we present the first results of the project, the properties of twenty DIBs in the 4100-5500 {\AA} range. We clearly detect a couple of previously elusive DIBs at 4170 {\AA} and 4591 {\AA}; the latter could have coronene and ovalene cations as carriers.

We present high-quality VLT-FLAMES optical spectroscopy of the nebular gas in the giant star-forming region 30 Doradus. In this paper, the first of a series, we introduce our observations and discuss the main kinematic features of 30 Dor, as revealed by the spectroscopy of the ionized gas in the region. The primary data set consists of regular grid of nebular observations, which we used to produce a spectroscopic datacube of 30 Dor, centered on the massive star cluster R136 and covering a field-of-view of 10'x10'. The main emission lines present in the datacube are from Halpha and [NII]6548,6584. The Halpha emission-line profile varies across the region from simple single-peaked emission to complex, multiple-component profiles, suggesting that different physical mechanisms are acting on the excited gas. To analyse the gas kinematics we fit Gaussian profiles to the observed Halpha features. Unexpectedly, the narrowest Halpha profile in our sample lies close to the supernova remnant 30 Dor B. We present maps of the velocity field and velocity dispersion across 30 Dor, finding five previously unclassified expanding structures. These maps highlight the kinematic richness of 30 Dor (e.g. supersonic motions), which will be analysed in future papers.

Building on results from the Magnetism in Massive Stars (MiMeS) project, this paper shows how a two-parameter classification of massive-star magnetospheres in terms of the magnetic wind confinement (which sets the Alfv\'en radius RA) and stellar rotation (which sets the Kepler co-rotation radius RK) provides a useful organisation of both observational signatures and theoretical predictions. We compile the first comprehensive study of inferred and observed values for relevant stellar and magnetic parameters of 64 confirmed magnetic OB stars with Teff > 16 kK. Using these parameters, we locate the stars in the magnetic confinement-rotation diagram, a log-log plot of RK vs. RA. This diagram can be subdivided into regimes of centrifugal magnetospheres (CM), with RA > RK, vs. dynamical magnetospheres (DM), with RK > RA. We show how key observational diagnostics, like the presence and characteristics of Halpha emission, depend on a star's position within the diagram, as well as other parameters, especially the expected wind mass-loss rates. In particular, we identify two distinct populations of magnetic stars with Halpha emission: namely, slowly rotating O-type stars with narrow emission consistent with a DM, and more rapidly rotating B-type stars with broader emission associated with a CM. For O-type stars, the high mass-loss rates are sufficient to accumulate enough material for line emission even within the relatively short free-fall timescale associated with a DM: this high mass-loss rate also leads to a rapid magnetic spindown of the stellar rotation. For the B-type stars, the longer confinement of a CM is required to accumulate sufficient emitting material from their relatively weak winds, which also lead to much longer spindown timescales. [abbreviated]

We are conducting a study on the imprint of the ISM on optical spectra based on two types of ongoing spectroscopic massive-star surveys: on the one hand, intermediate-resolution (R = 2500) green-blue spectra for ~3000 stars obtained with the Galactic O Star Spectroscopic Survey (GOSSS). On the other hand, high-resolution (R = 23 000 - 65 000) optical spectra for 600 stars obtained from three different surveys, OWN, IACOB, and NoMaDS. The R = 2500 data allows us to reach a larger sample with an average larger extinction while the R = 23 000 - 65 000 sample provides access to more diffuse interstellar bands (DIBs) and allows for the resolution in velocity of some ISM features. For each spectrum we are measuring the equivalent widths, FWHMs, and central wavelengths of 10-40 DIBs and interstellar lines (e.g. Ca II H+K, Na I D1+D2) and, in the case of GOSSS, the existence of an H II region around the star. We have also derived from auxiliary data or compiled from the literature values for the reddening, extinction law, H I column density, parallax, and H alpha emission. All of this constitutes the most complete collection ever of optical information on the ISM within 3 kpc of the Sun. We are analyzing the correlations between all of the collected quantities to discriminate between different possible origins of the DIBs.

This paper presents a first observational investigation of the faint Of?cp star NGC 1624-2, yielding important new constraints on its spectral and physical characteristics, rotation, magnetic field strength, X-ray emission and magnetospheric properties. Modeling the spectrum and spectral energy distribution, we conclude that NGC 1624-2 is a main sequence star of mass M {\simeq} 30 M{\odot}, and infer an effective temperature of 35 {\pm} 2 kK and log g = 4.0 {\pm} 0.2. Based on an extensive time series of optical spectral observations we report significant variability of a large number of spectral lines, and infer a unique period of 157.99 {\pm} 0.94 d which we interpret as the rotational period of the star. We report the detection of a very strong - 5.35 {\pm} 0.5 kG - longitudinal magnetic field <Bz>, coupled with probable Zeeman splitting of Stokes I profiles of metal lines confirming a surface field modulus <B> of 14 {\pm} 1 kG, consistent with a surface dipole of polar strength >~ 20 kG. This is the largest magnetic field ever detected in an O-type star, and the first report of Zeeman splitting of Stokes I profiles in such an object. We also report the detection of reversed Stokes V profiles associated with weak, high-excitation emission lines of O iii, which we propose may form in the close magnetosphere of the star. We analyze archival Chandra ACIS-I X-ray data, inferring a very hard spectrum with an X-ray efficiency log Lx/Lbol = -6.4, a factor of 4 larger than the canonical value for O-type stars and comparable to that of the young magnetic O-type star {\theta}1 Ori C and other Of?p stars. Finally, we examine the probable magnetospheric properties of the star, reporting in particular very strong magnetic confinement of the stellar wind, with {\eta}* {\simeq} 1.5 {\times} 10^4, and a very large Alfven radius, RAlf = 11.4 R*.

We describe observations in the nearby universe (<100 Mpc) with a 10-m or larger space-based telescope having imaging and spectral capabilities in the range 912-9000 \AA that would enable advances in the fields of massive stars, young populations, and star-forming galaxies, that are essential for achieving the Cosmic Origins Program objectives i) how are the chemical elements distributed in galaxies and dispersed in the circumgalactic and intergalactic medium; and ii) when did the first stars in the universe form, and how did they influence their environments. We stress the importance of observing hundreds of massive stars and their descendants individually, which will make it possible to separate the many competing factors that influence the observed properties of these systems (mass, composition, convection, mass-loss, rotation rate, binarity, magnetic fields, and cluster mass).

The last two decades have produced a proliferation of stellar atmosphere grids, evolutionary tracks, and isochrones which are available to the astronomical community from different internet services. However, it is not straightforward (at least for an inexperienced user) to manipulate those models to answer questions of the type: What is the spectral energy distribution of a 9000 K giant? What about its J-band magnitude for different metallicities? What can I tell about the mass of a star if I know that its unreddened B-V color is -0.05 and its luminosity in solar units is 10^5? The answers to those questions are indeed in the models but a series of transformations and combinations involving different variables and models are required to obtain them. To make the available knowledge more user friendly, I have combined a number of state-of-the-art sources to create a 3-D (effective temperature, luminosity, and metallicity) grid of stellar models for which I provide calibrated SEDs and magnitudes as well as auxiliary variables such as mass and age. Furthermore, I have generated a grid of extinguished magnitudes using the recent Ma\'iz Apell\'aniz et al. (2012) extinction laws and incorporated them into the Bayesian code CHORIZOS (Ma\'iz Apell\'aniz 2004).

I consider whether we can significantly improve the Cardelli et al. (1989) family of extinction laws using new data and techniques. There are six different aspects that need to be treated: The use of monochromatic quantities, the three different wavelength regimes (NIR, optical and UV), the sample, and the photometric calibration. Excluding the behavior in the NIR and UV, I discuss the other four aspects and propose a new family of extinction laws derived from VLT/FLAMES and HST/WFC3 data.

The Radiative Transfer Model (RTM) and the retrieval algorithm, incorporated in the SCIATRAN 2.2 software package developed at the Institute of Remote Sensing/Institute of Enviromental Physics of Bremen University (Germany), allows to simulate, among other things, radiance/irradiance spectra in the 2400-24 000 {\AA} range. In this work we present applications of RTM to two case studies. In the first case the RTM was used to simulate direct solar irradiance spectra, with different water vapor amounts, for the study of the water vapor content in the atmosphere above Sierra Nevada Observatory. Simulated spectra were compared with those measured with a spectrometer operating in the 8000-10 000 {\AA} range. In the second case the RTM was used to generate telluric model spectra to subtract the atmospheric contribution and correct high-resolution stellar spectra from atmospheric water vapor and oxygen lines. The results of both studies are discussed.

Detailed studies of resolved young massive star clusters are necessary to determine their dynamical state and evaluate the importance of gas expulsion and early cluster evolution. In an effort to gain insight into the dynamical state of the young massive cluster R136 and obtain the first measurement of its velocity dispersion, we analyse multi-epoch spectroscopic data of the inner regions of 30 Doradus in the Large Magellanic Cloud (LMC) obtained as part of the VLT-FLAMES Tarantula Survey. Following a quantitative assessment of the variability, we use the radial velocities of non-variable sources to place an upper limit of 6 km/s on the line-of-sight velocity dispersion of stars within a projected distance of 5 pc from the centre of the cluster. After accounting for the contributions of undetected binaries and measurement errors through Monte Carlo simulations, we conclude that the true velocity dispersion is likely between 4 and 5 km/s given a range of standard assumptions about the binary distribution. This result is consistent with what is expected if the cluster is in virial equilibrium, suggesting that gas expulsion has not altered its dynamics. We find that the velocity dispersion would be ~25 km/s if binaries were not identified and rejected, confirming the importance of the multi-epoch strategy and the risk of interpreting velocity dispersion measurements of unresolved extragalactic young massive clusters.

Although it has important ramifications for both the formation of star clusters and their subsequent dynamical evolution, rotation remains a largely unexplored characteristic of young star clusters (few Myr). Using multi-epoch spectroscopic data of the inner regions of 30 Doradus in the Large Magellanic Cloud (LMC) obtained as part of the VLT-FLAMES Tarantula Survey, we search for rotation of the young massive cluster R136. From the radial velocities of 36 apparently single O-type stars within a projected radius of 10 pc from the centre of the cluster, we find evidence, at the 95% confidence level, for rotation of the cluster as a whole. We use a maximum likelihood method to fit simple rotation curves to our data and find a typical rotational velocity of ~3 km/s. When compared to the low velocity dispersion of R136, our result suggests that star clusters may form with at least ~20% of the kinetic energy in rotation.

Whether massive stars can occasionally form in relative isolation or if they require a large cluster of lower-mass stars around them is a key test in the differentiation of star formation theories as well as how the initial mass function of stars is sampled. Previous attempts to find O-type stars that formed in isolation were hindered by the possibility that such stars are merely runaways from clusters, i.e., their current isolation does not reflect their birth conditions. We introduce a new method to find O-type stars that are not affected by such a degeneracy. Using the VLT-FLAMES Tarantula Survey and additional high resolution imaging we have identified stars that satisfy the following constraints: 1) they are O-type stars that are not detected to be part of a binary system based on RV time series analysis; 2) they are designated spectral type O7 or earlier ; 3) their velocities are within 1\sigma of the mean of OB-type stars in the 30 Doradus region, i.e. they are not runaways along our line-of-sight; 4) the projected surface density of stars does not increase within 3 pc towards the O-star (no evidence for clusters); 5) their sight lines are associated with gaseous and/or dusty filaments in the ISM, and 6) if a second candidate is found in the direction of the same filament with which the target is associated, both are required to have similar velocities. With these criteria, we have identified 15 stars in the 30 Doradus region, which are strong candidates for being high-mass stars that have formed in isolation. Additionally, we employed extensive MC stellar cluster simulations to confirm that our results rule out the presence of clusters around the candidates. Eleven of these are classified as Vz stars, possibly associated with the zero-age main sequence. We include a newly discovered W-R star as a candidate, although it does not meet all of the above criteria.

We present, for the first time, an optical spectroscopic data cube of the giant star-forming region 30 Doradus, obtained with the GIRAFFE on the VLT at Paranal Observatory. The main emission lines present in this data cube correspond to H{\alpha}, [NII] 6548 {\AA} and [NII] 6584 {\AA}. By using this data set, we found that H{\alpha} presents from simple to multiple profiles, which suggests that different physical mechanisms act in different ways on the excited gas in 30 Doradus. We found, at least, three unclassified large expanding structures. These structures correlate with peaks in the X-ray distribution. Given the excellent signal-to-noise ratio and the large spatial coverage of this data cube, we have studied in detail the kinematics of 30 Doradus, showing the importance of the small scale phenomena on the integrated properties of 30 Doradus.

Context: The mean UV extinction law for the Small Magellanic Cloud (SMC) is usually taken as a template for low-metallicity galaxies. However, its current derivation is based on only five stars, thus placing doubts on its universality. An increase in the number of targets with measured extinction laws in the SMC is necessary to determine its possible dependence on parameters such as metallicity and star-forming activity. Aims: To measure the UV extinction law for several stars in the quiescent molecular cloud SMC B1-1. Methods: We obtained HST/STIS slitless UV spectroscopy of a 25"x25" field of view and we combined it with ground-based NIR and visible photometry of the stars in the field. The results were processed using the Bayesian photometric package CHORIZOS to derive the visible-NIR extinction values for each star. The unextinguished Spectral Energy Distributions (SEDs) obtained in this way were then used to derive the UV extinction law for the four most extinguished stars. We also recalculated the visible-NIR extinction for the five SMC stars with preexisting UV extinction laws. Results: The UV extinction law for four SMC B1-1 stars within several pc of each other differs significantly from star to star. The 2175 {\AA} bump is moderately strong in one, weak in two, and absent in the fourth.

The IACOB spectroscopic survey of Galactic OB stars is an ambitious observational project aimed at compiling a large, homogeneous, high-resolution database of optical spectra of massive stars observable from the Northern hemisphere. The quantitative spectroscopic analysis of this database, complemented by the invaluable information provided by Gaia (mainly regarding photometry and distances), will represent a major step forward in our knowledge of the fundamental physical characteristics of Galactic massive stars. In addition, results from this analysis will be of interest for other scientific questions to be investigated using Gaia observations. In this contribution we outline the present status of the IACOB spectroscopic database and indicate briefly some of the synergy links between the IACOB and Gaia scientific projects.

There are several ongoing massive-star (mostly of spectral type O) surveys that are significantly increasing the quality and quantity of the spectroscopic information about these objects. Here we discuss and present results for two of them, GOSSS and NoMaDS. We also discuss recent and future developments on the Galactic O-Star Catalog and announce the upcoming availability of Marxist Ghost Buster, an IDL code that attacks spectral classification (hence the name) by using an interactive comparison with spectral libraries.

VFTS 682 is located in an active star-forming region, at a projected distance of 29 pc from the young massive cluster R136 in the Tarantula Nebula of the Large Magellanic Cloud. It was previously reported as a candidate young stellar object, and more recently spectroscopically revealed as a hydrogen-rich Wolf-Rayet (WN5h) star. Our aim is to obtain the stellar properties, such as its intrinsic luminosity, and to investigate the origin of VFTS 682. To this purpose, we model optical spectra from the VLT-FLAMES Tarantula Survey with the non-LTE stellar atmosphere code CMFGEN, as well as the spectral energy distribution from complementary optical and infrared photometry. We find the extinction properties to be highly peculiar (RV ~4.7), and obtain a surprisingly high luminosity log(L/Lsun) = 6.5 \pm 0.2, corresponding to a present-day mass of ~150Msun. The high effective temperature of 52.2 \pm 2.5kK might be explained by chemically homogeneous evolution - suggested to be the key process in the path towards long gamma-ray bursts. Lightcurves of the object show variability at the 10% level on a timescale of years. Such changes are unprecedented for classical Wolf-Rayet stars, and are more reminiscent of Luminous Blue Variables. Finally, we discuss two possibilities for the origin of VFTS 682: (i) the star either formed in situ, which would have profound implications for the formation mechanism of massive stars, or (ii) VFTS 682 is a slow runaway star that originated from the dense cluster R136, which would make it the most massive runaway known to date.

We report the discovery that R139 in 30 Doradus is a massive spectroscopic binary system. Multi-epoch optical spectroscopy of R139 was obtained as part of the VLT-FLAMES Tarantula Survey, revealing a double-lined system. The two components are of similar spectral types; the primary exhibits strong C III 4650 emission and is classified as an O6.5 Iafc supergiant, while the secondary is an O6 Iaf supergiant. The radial-velocity variations indicate a highly eccentric orbit with a period of 153.9 days. Photometry obtained with the Faulkes Telescope South shows no evidence for significant variability within an 18 month period. The orbital solution yields lower mass limits for the components of M1sin^3 i = 78 \pm 8 Msun and M2sin^3 i = 66 \pm 7 Msun. As R139 appears to be the most massive binary system known to contain two evolved Of supergiants, it will provide an excellent test for atmospheric and evolutionary models.

We show that star formation in the giant HII region NGC 5471 has been ongoing during the past 100 Myr. Using HST/WFPC2 F547M and F675W, ground based JHKs, and GALEX FUV and NUV images, we have conducted a photometric study of the star formation history in the massive giant extragalactic HII region NGC 5471 in M101. We perform a photometric study of the colour-magnitude diagram (CMD) of the resolved stars, and an integrated analysis of the main individual star forming clusters and of NGC 5471 as a whole. The integrated UV-optical-nIR photometry for the whole region provides two different reference ages, 8 Myr and 60 Myr, revealing a complex star formation history, clearly confirmed by the CMD resolved stellar photometry analysis. The spatial distribution of the stars shows that the star formation in NGC 5471 has proceeded along the whole region during, at least, the last 100 Myr. The current ionizing clusters are enclosed within a large bubble, which is likely to have been produced by the stars that formed in a major event < 20 Myr ago.

In this poster we show our preliminary analysis of DIBs (Diffuse Interstellar Bands) and other interstellar absorption lines with the purpose of understanding their origin and their relationship with extinction. We use the biggest Galactic O-star blue-violet spectroscopic sample ever (GOSSS, see contribution by Ma\'iz Apell\'aniz at this meeting). This sample allows a new insight on this topic because of the adequacy of O-star spectra, the sample number (700 and increasing, 400 used here), and their distribution in the MW disk. We confirm the high correlation coefficients between different DIBs and E(B - V), though the detailed behavior of each case shows small differences. We also detect a moderately low correlation coefficient between Ca II {\lambda}3934 (Ca K) and E(B - V) with a peculiar spatial distribution that we ascribe to the relationship between line saturation and velocity profiles for Ca II {\lambda}3934.

We present the first installment of a massive spectroscopic survey of Galactic O stars, based on new, high signal-to-noise ratio, R~2500 digital observations from both hemispheres selected from the Galactic O-Star Catalog of Ma\'iz Apell\'aniz et al. (2004) and Sota et al. (2008). The spectral classification system is rediscussed and a new atlas is presented, which supersedes previous versions. Extensive sequences of exceptional objects are given, including types Ofc, ON/OC, Onfp, Of?p, Oe, and double-lined spectroscopic binaries. The remaining normal spectra bring this first sample to 184 stars, which is close to complete to B=8 and north of delta = -20 and includes all of the northern objects in Ma\'iz Apell\'aniz et al. (2004) that are still classified as O stars. The systematic and random accuracies of these classifications are substantially higher than previously attainable, because of the quality, quantity, and homogeneity of the data and analysis procedures. These results will enhance subsequent investigations in Galactic astronomy and stellar astrophysics. In the future we will publish the rest of the survey, beginning with a second paper that will include most of the southern stars in Ma\'iz Apell\'aniz et al. (2004).

The Galactic O-Star Spectroscopic Survey (GOSSS) is an ambitious project that is observing all known Galactic O stars with B < 13 in the blue-violet part of the spectrum with R-2500. It is based on version 2 of the most complete catalog to date of Galactic O stars with accurate spectral types (v1, Ma\'iz Apell\'aniz et al. 2004 ;v2, Sota et al. 2008). Given the large amount of data that we are getting (more than 150 nights of observations at three different observatories in the last 4 years) we have developed an automatic spectroscopic reduction pipeline. This pipeline has been programmed in IDL and automates the process of data reduction. It can operate in two modes: automatic data reduction (quicklook) or semi-automatic data reduction (full). In "quicklook", we are able to get rectified and calibrated spectra of all stars of a full night just minutes after the observations. The pipeline automatically identifies the type of image and applies the standard reduction procedure (bias subtraction, flat field correction, application of bad pixel mask, ...). It also extracts all spectra of the stars in one image (including close visual binaries), aligns and merges all spectra of the same star (to increase the signal to noise ratio and to correct defects such as cosmic rays), calibrates in wavelength and rectifies the continuum. The same operations are performed in full mode, but allowing the user to adjust the parameters used in the process.

The Galactic O-Star Spectroscopic Survey (GOSSS) is a project that is observing all known Galactic O stars with B < 13 (~2000 objects) in the blue-violet part of the spectrum with R~2500. It also includes two companion surveys (a spectroscopic one at R~1500 and a high resolution imaging one). It is based on v2.0 of the Galactic O star catalog (v1, Ma\'iz Apell\'aniz et al. 2004; v2, Sota et al. 2008). We have completed the first part of the main project. Here we present results on the first 400 objects of the sample.

We present a massive spectroscopic survey of Galactic O stars, GOSSS, based on new, high signal-to-noise ratio, R~2500 blue-violet digital observations from both hemispheres. The sample size and selection criteria; the relationship between GOSSS, the Galactic O-Star Catalog (GOSC), and three sister surveys (OWN, IACOB, and Lucky Imaging); the current status; and our plans for the future are discussed. We also show some of our first results, which include the new Ofc category, two new examples of Of?p stars, a new atlas for O stars, and the introduction of the O9.7 type for luminosity classes III to V. Finally, our scientific objectives are discussed.

CONTEXT. Massive stars have high-multiplicity fractions, and many of them have still undetected components, thus hampering the study of their properties. AIMS. I study a sample of massive stars with high angular resolution to better characterize their multiplicity. METHODS. I observed 138 fields that include at least one massive star with AstraLux, a lucky imaging camera at the 2.2 m Calar Alto telescope. I also used observations of 3 of those fields with ACS/HRC on HST to obtain complementary information and to calibrate the AstraLux data. The results were compared with existing information from the Washington Double Star Catalog, Tycho-2, 2MASS, and other literature results. RESULTS. I discover 16 new optical companions of massive stars, the majority of which are likely to be physically bound to their primaries. I also improve the accuracy for the separation and magnitude difference of many previously known systems. In a few cases the orbital motion is detected when comparing the new data with existing ones and constraints on the orbits are provided. CONCLUSIONS. The analysis indicate that the majority of the AstraLux detections are bound pairs. For a range of separations of 0.1"-14" and magnitude differences lower than 8, I find that the multiplicity fraction for massive stars is close to 50%. When objects outside those ranges are included, the multiplicity fraction should be considerably higher.

I analyze the stochastic effects introduced by the sampling of the stellar initial mass function (SIMF) in the derivation of the individual masses and the cluster mass function (CMF) from broad-band visible-NIR unresolved photometry. The classical method of using unweighted UBV photometry to simultaneously establish ages and extinctions of stellar clusters is found to be unreliable for clusters older than approx. 30 Ma, even for relatively large cluster masses. On the other hand, augmenting the filter set to include longer-wavelength filters and using weights for each filter increases the range of masses and ages that can be accurately measured with unresolved photometry. Nevertheless, a relatively large range of masses and ages is found to be dominated by SIMF sampling effects that render the observed masses useless, even when using UBVRIJHK photometry.

It is currently common to use spatially unresolved multi-filter broad-band photometry to determine the masses of individual stellar clusters (and hence the cluster mass function, CMF). I analyze the stochastic effects introduced by the sampling of the stellar initial mass function (SIMF) in the derivation of the individual masses and the CMF and I establish that such effects are the largest contributor to the observational uncertainties. An analytical solution, valid in the limit where uncertainties are small, is provided to establish the range of cluster masses over which the CMF slope can be obtained with a given accuracy. The validity of the analytical solution is extended to higher mass uncertainties using Monte Carlo simulations and the Gamma approximation. The value of the Poisson mass is calculated for a large range of ages and a variety of filters for solar-metallicity clusters measured with single-filter photometry. A method that uses the code CHORIZOS is presented to simultaneously derive masses, ages, and extinctions. The classical method of using unweighted UBV photometry to simultaneously establish ages and extinctions of stellar clusters is found to be unreliable for clusters older than approx. 30 Ma, even for relatively large cluster masses. On the other hand, augmenting the filter set to include longer-wavelength filters and using weights for each filter increases the range of masses and ages that can be accurately measured with unresolved photometry. Nevertheless, a relatively large range of masses and ages is found to be dominated by SIMF sampling effects that render the observed masses useless, even when using UBVRIJHK photometry. A revision of some literature results affected by these effects is presented and possible solutions for future observations and analyses are suggested.

We use the new reduction of the Hipparcos data (van Leeuwen 2007) and a self-consistent distance determination technique for Lutz-Kelker limited samples to obtain distances to the massive stars in the solar vicinity. The distance uncertainties for the nearest massive stars have been substantially reduced with respect to those derived from the old Hipparcos reduction. In two cases (gamma2 Vel and theta2 Ori A) we have been able to verify that the new values are in good agreement with recent determinations with alternative methods. We also derive new values for the vertical displacement of the Sun with respect to the Galactic Plane and for the scale height of the thin disk from the spatial distribution of massive stars around us.

I discuss two of the possible sources of biases in the determination of the IMF: binning and the existence of unresolved components. The first source is important for clusters with a small number of stars detected in a given mass bin while the second one is relevant for all clusters located beyond the immediate solar neighborhood. For both cases I will present results of numerical simulations and I will discuss strategies to correct for their effects. I also present a brief description of a third unrelated bias source.

When calculating IMFs for young clusters, one has to take into account that (a) most massive stars are born in multiple systems (b) most IMFs are derived from data that cannot resolve such systems, and (c) multiple chance superpositions between members are expected to happen if the cluster is too distant. In this article I use numerical experiments to model the consequences of those phenomena on the observed color-magnitude diagrams and the IMFs derived from them. Real multiple systems affect the observed or apparent massive-star MF slope little but can create a significant population of apparently ultramassive stars. Chance superpositions produce only small biases when the number of superimposed stars is low but, once a certain number threshold is reached, they can affect both the observed slope and the apparent stellar upper mass limit. I apply those experiments to two well known massive young clusters in the Local Group, NGC 3603 and R136. In both cases I show that the observed population of stars with masses above 120 solar masses can be explained by the effects of unresolved objects, mostly real multiple systems for NGC 3603 and a combination of real and chance-alignment multiple systems for R136. Therefore, the case for the reality of a stellar upper mass limit at solar or near-solar metallicities is strengthened, with a possible value even lower than 150 solar masses. An IMF slope somewhat flatter than Salpeter or Kroupa with gamma between -1.6 and -2.0 is derived for the central region of NGC 3603, with a significant contribution to the uncertainty arising from the imprecise knowledge of the distance to the cluster. The IMF at the very center of R136 cannot be measured with the currently available data but the situation could change with new HST observations. (abridged)

The Galactic O star catalog (GOS) is an ambitious project to provide as much information regarding these types of objects as possible. The first version of the catalog (GOS v1) included data for 378 stars with precise spectral classification. It was intended to be complete up to V < 8, but also included many stars fainter than that limit. In this new version, we include a second list with more than 700 stars that have sometimes been classified as O stars. The catalog includes cross-references to other catalogs, spectral classification (various references for each star on the second list), coordinates (obtained in most cases from Tycho-2 data, but some from 2MASS and USNO B1), astrometric distances for some of the nearer stars, Optical and NIR photometry (Tycho-2, Johnson, Stromgren, Cousin and 2MASS), group membership, runaway character, and multiplicity information.

I have recently combined HST/STIS spectrophotometry with existing photometric data to analyze the calibration of three standard optical photometry systems: Tycho-2 B_TV_T, Stromgren uvby, and Johnson UBV. In this contribution I summarize those results, present new ones for 2MASS JHK_s, and combine them with recent literature results to generate a uniform set of zero points for six photometric systems, the above mentioned plus Cousins RI and SDSS ugriz. With the exception of the latter system, the zero points use the new Vega spectrum presented at this meeting by Ralph Bohlin. I also discuss the implementation of these results in CHORIZOS, a Bayesian photometric code that compares multi-filter observational data with spectral energy distributions to solve the inverse problem of finding the models which are compatible with the observations.

It has been traditionally stated that it is not possible to use optical/NIR photometry to measure the temperatures of O stars. In this contribution we describe the steps required to overcome the hurdles that have prevented this from happening in the past and we present our preliminary results for the low-extinction case.

We are using HST GO programs 10205, 10602, and 10898 to test the stellar upper mass limit in the solar vicinity by attempting to detect optical close companions, thus lowering the calculated evolutionary masses. We have observed with ACS/HRC all the known (as of early 2005) Galactic O2/3/3.5 stars. We also have observations with HST/FGS and ground-based spectroscopy from LCO and CASLEO. Here we discuss our results for Pismis 24 and HD 93129A.

Is there a stellar upper mass limit? Recent statistical work seems to indicate that there is and that it is in the vicinity of 150 solar masses. In this paper we use HST and ground-based data to investigate the brightest members of the cluster Pismis 24, one of which (Pismis 24-1) was previously inferred to have a mass greater than 200 solar masses, in apparent disagreement with that limit. We determine that Pismis 24-1 is composed of at least three objects, the resolved Pismis 24-1SW and the unresolved spectroscopic binary Pismis 24-1NE. The evolutionary zero-age masses of Pismis 24-1SW, the unresolved system Pismis 24-1NE, and the nearby star Pismis 24-17 are all approximately 100 solar masses, very large but under the stellar upper mass limit.

I use HST spectrophotometry to analyze the calibration of 3 optical photometry systems: Tycho-2 B_T V_T, Str\"omgren u v b y, and Johnson U B V. For Tycho-2, I revisit the analysis of Ma\'iz Apell\'aniz (2005) to include the new recalibration of grating/aperture corrections, vignetting, and charge transfer inefficiency effects produced by the STIS group and to consider the consequences of both random and systematic uncertainties. The new results reaffirm the good quality of both the Tycho-2 photometry and the HST spectrophotometry, but yield a slightly different value for ZP_B_T-V_T. For the Str\"omgren v, b, and y and the Johnson B and V filters I find that the published sensitivity curves are consistent with the available photometry and spectrophotometry and I derive new values for the associated ZP_b-y, ZP_m_1, and ZP_B-V. The situation is different for the Str\"omgren u and the Johnson U filters. There I find that the published sensitivity curves yield results that are inconsistent with the available photometry and spectrophotometry, likely caused by an incorrect treatment of atmospheric effects in the short-wavelength end. I reanalyze the data to produce new average sensitivity curves for those two filters and new values for ZP_c_1 and ZP_U-B. The new computation of synthetic U-B and B-V colors uses a single B sensitivity curve, which eliminates the previous unphysical existence of two different definitions. Finally, I find that if one uses values from the literature where uncertainties are not given, reasonable estimates for these are 1-2% for Str\"omgren b-y, m_1, and c_1 and 2-3% for Johnson B-V and U-B. The use of the results in this article should lead to a significant reduction of systematic errors when comparing synthetic photometry models with real colors and indices.

I study the importance of the accurate calibration of photometric systems in order to produce meaningful comparisons between the observed colors + magnitudes and model SEDs. Possible sources of errors are discussed and two examples are analyzed. I show that well-calibrated Tycho-2 photometry is stable and precise enough for such comparisons. On the contrary, the available calibrations for Johnson UBV photometry yield relative large systematic errors, which has prompted me to develop a new, more precise calibration. The advantages of multicolor photometry over the standard single-color + magnitude diagrams for the derivation of physical properties of stars (elimination of degeneracies, inclusion of multiple parameters, avoidance of linearizing approximations, possibility of a more precise treatment of errors) are discussed through the use of CHORIZOS, a code developed specifically for this purpose.

We discuss possible sources of biases in the determination of the initial mass function (IMF) introduced by the binning of the data, the uncertainty in the determinations of masses, and the existence of unresolved multiple systems. Those three effects tend to produce IMFs that are flatter than the real one. We analyze the importance of each effect and suggest techniques that minimize or eliminate the biases. We also report the detection of the first astrometric binary system composed of two very-early O-type stars, HD 93129 A.

We detect and quantify significant numerical biases in the determination of the slope of power laws with Salpeter (or similar) indices from uniformly-binned data using chi-square minimization. The biases are caused by the correlation between the number of stars per bin and the assigned weights and are especially important when the number of stars per bin is small. This result implies the existence of systematic errors in the values of IMFs calculated in this way. We propose as an alternative using variable-size bins and dividing the stars evenly among them. Such variable-size bins yield very small biases that are only weakly dependent on the number of stars per bin. Furthermore, we show that they allow for the calculation of reliable IMFs with only a small total number of stars. Therefore, they are a preferred alternative to the standard uniform-size binning.

We present a method designed to correct for Lutz-Kelker effects in distance-limited samples. The method allows for the calculation of distances to individual objects and, at the same time, provides a fit to a parameterized, self-consistent spatial distribution of the population. An example using Hipparcos data is presented and the relevance to Gaia is also discussed.