The apparent lack of massive O-type stars near the zero-age main sequence (at ages < 2 Myr) is a topic widely discussed. Different explanations for this elusive detection have been proposed, but no firm conclusions have been reached yet. We reassess this empirical result benefiting from the high-quality spectroscopic observations of >400 Galactic O-type stars gathered by the IACOB and OWN surveys. We used temperatures and gravities from a iacob-gbat/fastwind spectroscopic analysis to locate our sample in the Kiel and spectroscopic HR diagrams. We evaluated the completeness of our sample of stars, observational biases using information from the Galactic O star catalog (GOSC), systematics of our methodology, and compare with other recent studies using smaller samples of Galactic O-type stars. We base our discussion on the spectroscopic HR diagram to avoid the use of uncertain distances. We performed a detailed study of the young cluster Trumpler-14 as an example of how Gaia cluster distances can help to construct the associated classical HR diagram. The apparent lack of massive O-type stars near the ZAMS with masses between 30 and 70 Msol persist even when spectroscopic results from a large, non-biased sample of stars are used. We do not find correlation between the dearth of stars and observational biases, limitations of our methodology, or the use of spectroscopic HR diagram instead of the classical one. Investigating the efficiency of mass accretion during the formation process we conclude that an adjustment of the accretion rate towards lower values could reconcile the hotter boundary of detected O-type stars and the theoretical birthline. Last, we discuss that the presence of a small sample of O2-O3.5 stars found closer to the ZAMS might be explained taking into account non-standard star evolution (e.g. binary interaction, mergers, or homogeneous evolution).

Lack of high-precision long-term continuous photometric data for large samples of stars has prevented the large-scale exploration of pulsational variability in the OB star regime. As a result, the candidates for in-depth asteroseismic modelling remained limited to a few tens of dwarfs. The TESS nominal space mission has surveyed the southern sky, yielding continuous data of at least 27 d for hundreds of OB stars. We aim to couple TESS data in the southern sky with spectroscopy to study the variability over mass and evolution. We focus mainly on the presence of coherent pulsation modes that may or may not be present in the theoretical instability domains and unravel all frequency behaviour in the amplitude spectra of the TESS data. We compose a sample of 98 OB-type stars observed by TESS in Sectors 1-13 and with available high-resolution spectroscopy gathered by the IACOB and OWN surveys. We present the short-cadence 2-min light curves of dozens of OB-type stars, that have one or more spectra in the IACOB or OWN database. Based on these light curves and their Lomb-Scargle periodograms we perform variability classification and frequency analysis, and place the stars in the spectroscopic Hertzsprung-Russell diagram to interpret the variability in an evolutionary context. We deduce diverse origins of the variability found in all of the 98 OB stars in the TESS data. Among these we find several new variable stars, including three hybrid pulsators, three eclipsing binaries, high frequency modes in a Be star, and potential heat-driven pulsations in two Oe stars. We identify stars for which future asteroseismic modelling is possible, provided mode identification is achieved. By comparing the position of the variables to theoretical instability strips we discuss the current shortcomings in non-adiabatic pulsation theory, and the distribution of pulsators in the upper Hertzsprung-Russell diagram.

The recent gravitational wave measurements have demonstrated the existence of stellar mass black hole binaries. It is essential for our understanding of massive star evolution to identify the contribution of binary evolution to the formation of double black holes. A promising way to progress is investigating the progenitors of double black hole systems and comparing predictions with local massive star samples such as the population in 30 Doradus in the Large Magellanic Cloud (LMC). Methods. To this purpose, we analyse a large grid of detailed binary evolution models at LMC metallicity with initial primary masses between 10 and 40 Msun, and identify which model systems potentially evolve into a binary consisting of a black hole and a massive main sequence star. We then derive the observable properties of such systems, as well as peculiarities of the OB star component. We find that about 3% of the LMC late O and early B stars in binaries are expected to possess a black hole companion, when assuming stars with a final helium core mass above 6.6 M to form black holes. While the vast majority of them may be X-ray quiet, our models suggest that these may be identified in spectroscopic binaries, either by large amplitude radial velocity variations ( > 50 km s ) and simultaneous nitrogen surface enrichment, or through a moderate radial velocity ( > 10 km/s ) and simultaneously rapid rotation of the OB star. The predicted mass ratios are such that main sequence companions could be excluded in most cases. A comparison to the observed OB+WR binaries in the LMC, Be/X-ray binaries, and known massive BH binaries supports our conclusion. We expect spectroscopic observations to be able to test key assumptions in our models, with important implications for massive star evolution in general, and for the formation of double-black hole mergers in particular.

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.

Star clusters are the building blocks of galaxies. They are composed of stars of nearly equal age and chemical composition, allowing us to use them as chronometers and as testbeds for gauging stellar evolution. It has become clear recently that massive stars are formed preferentially in close binaries, in which mass transfer will drastically change the evolution of the stars. This is expected to leave a significant imprint in the distribution of cluster stars in the Hertzsprung-Russell diagram. Our results, based on a dense model grid of more than 50,000 detailed binary-evolution calculations, indeed show several distinct, coeval main-sequence (MS) components, most notably an extended MS turnoff region, and a group of near-critical rotating stars that is spread over a large luminosity range on the red side of the classical MS. We comprehensively demonstrate the time evolution of the features in an animation, and we derive analytic expressions to describe these features. We find quantitative agreement with results based on recent photometric and spectroscopic observations. We conclude that while other factors may also be at play, binary evolution has a major impact on the MS morphology of young star clusters.

We report the discovery of a square axisymmetric circumstellar nebula around the emission-line star HD 93795 in archival Spitzer Space Telescope 24 micron data. We classify HD 93795 as an B9 Ia star using optical spectra obtained with the Southern African Large Telescope (SALT). A spectral analysis carried out with the stellar atmosphere code FASTWIND indicates that HD 93795 only recently left the main sequence and is evolving redward for the first time. We discuss possible scenarios for the origin of the nebula and suggest that HD 93795 was originally a binary system and that the nebula was formed because of merger of the binary components. We also discuss a discrepancy between distance estimates for HD 93795 based on the Gaia data and the possible membership of this star of the Car OB1 association, and conclude that HD 93795 could be at the same distance as Car OB1.

High-performance instruments at large ground-based telescopes have made integral field spectroscopy (IFS) a powerful tool for the study of extended objects such as galaxies, nebulae, or even larger survey fields on the sky. Here we discuss the capabilities of IFS for the study of resolved stellar populations, using the new method of PSF-fitting crowded field IFS, analogous to the well-established technique of crowded field photometry with image sensors. We review early pioneering work with first generation integral field spectrographs, the breakthrough achieved with the MUSE instrument at the ESO Very Large Telescope, the remarkable progress accomplished with MUSE in the study of globular clusters, and first results on nearby galaxies. We discuss the synergy of integral field spectrographs at 8-10m class telescopes with future facilities such as the Extremely Large Telescope (ELT).

It was recently claimed that the magnetic O-type star HD 54879 exhibits important radial velocity variability indicative of its presence in a spectroscopic binary. More remarkably, it was furthermore reported that the star underwent a short, sudden variation in spectral type and magnetic field. In this Letter we examine new Narval and ESPaDOnS data of this star in addition to the previously-published FORS2 data and conclude that both the reported velocity variations and the sudden spectral and magnetic changes are spurious.

We present the concept of BlueMUSE, a blue-optimised, medium spectral resolution, panoramic integral field spectrograph based on the MUSE concept and proposed for the Very Large Telescope. With an optimised transmission down to 350 nm, a larger FoV (1.4 x 1.4 arcmin$^2$) and a higher spectral resolution compared to MUSE, BlueMUSE will open up a new range of galactic and extragalactic science cases allowed by its specific capabilities, beyond those possible with MUSE. For example a survey of massive stars in our galaxy and the Local Group will increase the known population of massive stars by a factor $>$100, to answer key questions about their evolution. Deep field observations with BlueMUSE will also significantly increase samples of Lyman-alpha emitters, spanning the era of Cosmic Noon. This will revolutionise the study of the distant Universe: allowing the intergalactic medium to be detected unambiguously in emission, enabling the study of the exchange of baryons between galaxies and their surroundings. By 2030, at a time when the focus of most of the new large facilities (ELT, JWST) will be on the infra-red, BlueMUSE will be a unique facility, outperforming any ELT instrument in the Blue/UV. It will have a strong synergy with ELT, JWST as well as ALMA, SKA, Euclid and Athena.

NGC 2345 is a young open cluster hosting seven blue and red supergiants, low metallicity and a high fraction of Be stars which makes it a privileged laboratory to study stellar evolution. We aim to improve the determination of the cluster parameters and study the Be phenomenon. Our objective is also to characterise its seven evolved stars by deriving their atmospheric parameters and chemical abundances. We performed a complete analysis combining for the first time $ubvy$ photometry with spectroscopy as well as $Gaia$ Data Release 2. We obtained spectra with classification purposes for 76 stars and high-resolution spectroscopy for an in-depth analysis of the blue and red evolved stars. We identify a new red supergiant and 145 B-type likely members within a radius of 18.7$\pm$1.2 arcmin, which implies an initial mass, $M_{\textrm{cl}}\approx$5200 M$_{\odot}$. We find a distance of 2.5$\pm$0.2 kpc for NGC 2345, placing it at $R_{\textrm{GC}}$=10.2$\pm$0.2 kpc. Isochrone fitting supports an age of 56$\pm$13 Ma, implying masses around 6.5 M$_{\odot}$ for the supergiants. A high fraction of Be stars ($\approx$10$\%$) is found. From the spectral analysis we estimate for the cluster an average $v_{\textrm{rad}}$=$+58.6\pm0.5$ kms$^{-1}$ and a low metallicity, [Fe/H]=$-$0.28$\pm$0.07. We also have determined chemical abundances for Li, O, Na, Mg, Si, Ca, Ti, Ni, Rb, Y, and Ba for the evolved stars. The chemical composition of the cluster is consistent with that of the Galactic thin disc. One of the K supergiants, S50, is a Li-rich star, presenting an A(Li)$\approx$2.1. An overabundance of Ba is found, supporting the enhanced $s$-process. NGC 2345 has a low metallicity for its Galactocentric distance, comparable to typical LMC stars. It is massive enough to serve as a testbed for theoretical evolutionary models for massive intermediate-mass stars.

Cosmic History has witnessed the lives and deaths of multiple generations of massive stars, all of them invigorating their host galaxies with ionizing photons, kinetic energy, fresh material and stellar-mass black holes. Ubiquitous engines as they are, Astrophysics needs a good understanding of their formation, evolution, properties and yields throughout the history of the Universe, and with decreasing metal content mimicking the environment at the earliest epochs. Ultimately, a physical model that could be extrapolated to zero metallicity would enable tackling long-standing questions such as "What did the First, very massive stars of the Universe look like?" or "What was their role in the re-ionization of the Universe?". Yet, most our knowledge of metal-poor massive stars is drawn from one single point in metallicity. Massive stars in the Small Magellanic Cloud (SMC, $\sim 1/5 Z_{\odot}$) currently serve as templates for low-metallicity objects in the early Universe, even though significant differences with respect to massive stars with poorer metal content have been reported. This White Paper summarizes the current knowledge on extremely (sub-SMC) metal poor massive stars, highlighting the most outstanding open questions and the need to supersede the SMC as standard. A new paradigm can be built from nearby extremely metal-poor galaxies that make a new metallicity ladder, but massive stars in these galaxies are out of reach to current observational facilities. Such task would require an L-size mission, consisting of a 10m-class space telescope operating in the optical and the ultraviolet ranges. Alternatively, we propose that ESA unites efforts with NASA to make the LUVOIR mission concept a reality, thus continuing the successful partnership that made Hubble Space Telescope one of the greatest observatories of all time.

Almost all massive stars explode as supernovae and form a black hole or neutron star. The remnant mass and the impact of the chemical yield on subsequent star formation and galactic evolution strongly depend on the internal physics of the progenitor star, which is currently not well understood. The theoretical uncertainties of stellar interiors accumulate with stellar age, which is particularly pertinent for the blue supergiant phase. Stellar oscillations represent a unique method of probing stellar interiors, yet inference for blue supergiants is hampered by a dearth of observed pulsation modes. Here we report the detection of diverse variability in blue supergiants using the K2 and TESS space missions. The discovery of pulsation modes or an entire spectrum of low-frequency gravity waves in these stars allow us to map the evolution of hot massive stars towards the ends of their lives. Future asteroseismic modelling will provide constraints on ages, core masses, interior mixing, rotation and angular momentum transport. The discovery of variability in blue supergiants is a step towards a data-driven empirical calibration of theoretical evolution models for the most massive stars in the Universe.

Multiple generations of massive stars have lived and died during Cosmic History, invigorating host galaxies with ionizing photons, kinetic energy, fresh material and stellar-size black holes. At present, massive stars in the Small Magellanic Cloud (SMC) serve as templates for low-metallicity objects in the early Universe. However, recent results have highlighted important differences in the evolution, death and feedback of massive stars with poorer metal content that better matches the extremely low metallicity of previous Cosmic epochs. This paper proposes to supersede the SMC standard with a new metallicity ladder built from very metal-poor galaxies, and provides a brief overview of the technological facilities needed to this aim.

We report that LAMOST-HVS1 is a massive hyper-runaway subgiant star with mass of 8.3 Msun and super-Solar metallicity, ejected from the inner stellar disk of the Milky Way $\sim$ 33 Myr ago with the intrinsic ejection velocity of $568^{+19}_{-17}$ km/s (corrected for the streaming motion of the disk), based on the proper motion data from Gaia Data Release 2 (DR2) and high-resolution spectroscopy. The extremely large ejection velocity indicates that this star was not ejected by the supernova explosion of the binary companion. Rather, it was probably ejected by a 3- or 4-body dynamical interaction with more massive objects in a high-density environment. Such a high-density environment may be attained at the core region of a young massive cluster with mass of $\gtrsim 10^4$ Msun. The ejection agent that took part in the ejection of LAMOST-HVS1 may well be an intermediate mass black hole ($\gtrsim$ 100 Msun), a very massive star ($\gtrsim$ 100 Msun), or multiple ordinary massive stars ($\gtrsim$ 30 Msun). Based on the flight time and the ejection location of LAMOST-HVS1, we argue that its ejection agent or its natal star cluster is currently located near the Norma spiral arm. The natal star cluster of LAMOST-HVS1 may be an undiscovered young massive cluster near the Norma spiral arm.

We present the first stellar spectroscopy in the low-luminosity (M_V ~-9.3 mag), dwarf galaxy Leo P. Its significantly low oxygen abundance (3% solar) and relative proximity (~1.6 Mpc) make it a unique galaxy to investigate the properties of massive stars with near-primordial compositions akin to those in the early Universe. From our VLT-MUSE spectroscopy we find the first direct evidence for an O-type star in the prominent HII region, providing an important test case to investigate the potential environmental dependence of the upper end of the initial mass function in the dwarf galaxy regime. We classify 14 further sources as massive stars (and 17 more as candidate massive stars), most likely B-type objects. From comparisons with published evolutionary models we argue that the absolute visual magnitudes of massive stars in very metal-poor systems such as Leo P and I Zw 18 may be fainter by ~0.5 mag compared to Galactic stars. We also present spectroscopy of two carbon stars identified previously as candidate asymptotic-giant-branch stars. Two of three further candidate asymptotic-giant-branch stars display CaII absorption, confirming them as cool, evolved stars; we also recover CaII absorption in the stacked data of the next brightest 16 stars in the upper red giant branch. These discoveries will provide targets for future observations to investigate the physical properties of these objects and to calibrate evolutionary models of luminous stars at such low metallicity. The MUSE data also reveal two 100pc-scale ring structures in Halpha emission, with the HII region located on the northern edge of the southern ring. Lastly, we report serendipitous observations of 20 galaxies, with redshifts ranging from z=0.39, to a close pair of star-forming galaxies at z=2.5.

We report the discovery of two mid-infrared nebulae in the northern hemisphere with the Wide-field Infrared Survey Explorer and the results of optical spectroscopy of their central stars, BD+60 2668 (composed of two components, separated from each other by \approx 3 arcsec) and ALS 19653, with the Calar Alto 3.5-m telescope and the Southern African Large Telescope (SALT), respectively. We classify the components of BD+60 2668 as stars of spectral types B0.5 II and B1.5 III. ALS 19653 is indicated in the SIMBAD data base as a planetary nebula, while our observations show that it is a massive B0.5 Ib star, possibly in a binary system. Using the stellar atmosphere code FASTWIND, we derived fundamental parameters of the three stars as well as their surface element abundances, implying that all of them are either on the main sequence or only recently left it. This provides further evidence that massive stars can produce circumstellar nebulae while they are still relatively unevolved. We also report the detection of optical counterparts to the mid-infrared nebulae and a second, more extended optical nebula around ALS 19653, and present the results of SALT spectroscopy of both nebulae associated with this star. The possible origin of the nebulae is discussed.

We report the results of our study of the high-velocity (\approx160 km/s) runaway O star CPD-64 2731 and its associated horseshoe-shaped nebula discovered with the Wide-field Infrared Survey Explorer. Spectroscopic observations with the Southern African Large Telescope and spectral analysis indicate that CPD-64 2731 is a fast-rotating main-sequence O5.5 star with enhanced surface nitrogen abundance. We derive a projected rotational velocity of \approx300 km/s which is extremely high for this spectral type. Its kinematic age of \approx6 Myr, assuming it was born near the Galactic plane, exceeds its age derived from single star models by a factor of two. These properties suggest that CPD-64 2731 is a rejuvenated and spun-up binary product. The geometry of the nebula and the almost central location of the star within it argue against a pure bow shock interpretation for the nebula. Instead, we suggest that the binary interaction happened recently, thereby creating the nebula, with a cavity blown by the current fast stellar wind. This inference is supported by our results of 2D numerical hydrodynamic modelling.

Gravitational 3-body interaction among binary stars and the supermassive black hole (SMBH) at the center of the Milky Way occasionally ejects a hypervelocity star (HVS) with a velocity of ~1000 km/s. Due to the ejection location, such a HVS initially has negligible azimuthal angular momentum Lz ~ 0 kpc km/s. Even if the halo is mildly triaxial, Lz of a recently ejected nearby HVS remains negligible, since its flight time from the Galactic Center is too short to accumulate noticeable torque. However, if we make a wrong assumption about the Solar position and velocity, such a HVS would apparently have noticeable non-zero azimuthal angular momentum, due to the wrong reflex motion of the Sun. Conversely, with precise astrometric data for a nearby HVS, we can measure the Solar position and velocity by assuming that the HVS has zero azimuthal angular momentum. Based on this idea, here we propose a method to estimate the Galactocentric distance of the Sun R0 and the Galactocentric Solar azimuthal velocity Vsun by using a HVS. We demonstrate with mock data for a nearby HVS candidate that the Gaia astrometric data, along with the currently available constraint on Vsun/R0 from the proper motion measurement of Sgr A*, can constrain R0 and Vsun with uncertainties of ~0.27 kpc and ~7.8 km/s (or fractional uncertainties of 3 percent), respectively. Our method will be a promising tool to constrain (R0, Vsun), given that Gaia is expected to discover many nearby HVSs in the near future.

We use GAIA DR2 proper motions of the RIOTS4 field OB stars in the Small Magellanic Cloud (SMC) to study the kinematics of runaway stars. The data reveal that the SMC Wing has a systemic peculiar motion relative to the SMC Bar of (v_RA, v_Dec) = (62 +/-7, -18+/-5) km/s and relative radial velocity +4.5 +/- 5.0 km/s. This unambiguously demonstrates that these two regions are kinematically distinct: the Wing is moving away from the Bar, and towards the Large Magellanic Cloud with a 3-D velocity of 64 +/- 10 km/s. This is consistent with models for a recent, direct collision between the Clouds. We present transverse velocity distributions for our field OB stars, confirming that unbound runaways comprise on the order of half our sample, possibly more. Using eclipsing binaries and double-lined spectroscopic binaries as tracers of dynamically ejected runaways, and high-mass X-ray binaries (HMXBs) as tracers of runaways accelerated by supernova kicks, we find significant contributions from both populations. The data suggest that HMXBs have lower velocity dispersion relative to dynamically ejected binaries, consistent with the former corresponding to less energetic supernova kicks that failed to unbind the components. Evidence suggests that our fast runaways are dominated by dynamical, rather than supernova, ejections.

We present a comprehensive stellar atmosphere analysis of 329 O- and B-type stars in the Small Magellanic Cloud (SMC) from the RIOTS4 survey. Using spectroscopically derived effective temperature (Teff) and surface gravities, we find that classical Be stars appear misplaced to low Teff and high luminosity in the spectroscopic Hertzsprung-Russell diagram (sHRD). Together with the most luminous stars in our sample, the stellar masses derived from the sHRD for these objects are systematically larger than those obtained from the conventional HRD. This suggests that the well-known, spectroscopic mass-discrepancy problem may be linked to the fact that both groups of stars have outer envelopes that are nearly gravitationally unbound. The non-emission-line stars in our sample mainly appear on the main-sequence, allowing a first estimate of the terminal-age main-sequence (TAMS) in the SMC, which matches the predicted TAMS between 12 and 40$\,$M$_{\odot}$ at SMC metallicity. We further find a large underabundance of stars above $\sim 25\,$M$_{\odot}$ near the ZAMS, reminiscent of such earlier findings in the Milky Way and LMC.

We focus on early-B type stars with helium overabundance, for which the presence of a magnetic field has not previously been reported. The measurements were carried out using high-spectral-resolution spectropolarimetric observations obtained with the High Accuracy Radial velocity Planet Searcher (HARPS) in polarimetric mode, installed at the ESO La Silla 3.6m telescope. For five He-rich stars, the longitudinal magnetic field was detected for the first time. For one target, HD58260, the presence of a longitudinal magnetic field of the order of 1.8 kG has already been reported in the literature, but the magnetic field has remained constant over tens of years. Our measurement carried out using the polarimetric spectra obtained in 2015 March indicates a slight decrease of the longitudinal magnetic field strength compared to measurements reported in previous works. A search for periodic modulation in available photometric data allowed us to confidently establish a period of 2.64119$\pm$0.00420 d in archival ASAS3 data for CPD -27$^{\circ}$ 1791. No period could be determined for the other five stars. The obtained results support the scenario that all He-rich stars are detectably magnetic and form an extension of the Ap star phenomenon to higher temperatures.

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.

NGC 3105 is a young open cluster hosting blue, yellow and red supergiants. This rare combination makes it an excellent laboratory to constrain evolutionary models of high-mass stars. It is poorly studied and fundamental parameters such as its age or distance are not well defined. We intend to characterize in an accurate way the cluster as well as its evolved stars, for which we derive for the first time atmospheric parameters and chemical abundances. We identify 126 B-type likely members within a radius of 2.7$\pm$0.6 arcmin, which implies an initial mass, $M_{cl}\approx$4100 M$_{\odot}$. We find a distance of 7.2$\pm$0.7 kpc for NGC 3105, placing it at $R_{GC}$=10.0$\pm$1.2 kpc. Isochrone fitting supports an age of 28$\pm$6 Ma, implying masses around 9.5 M$_{\odot}$ for the supergiants. A high fraction of Be stars ($\approx$25 %) is found at the top of the main sequence down to spectral type b3. From the spectral analysis we estimate for the cluster a $v_{rad}$=+46.9$\pm$0.9 km s$^{-1}$ and a low metallicity, [Fe/H]=-0.29$\pm$0.22. We also have determined, for the first time, chemical abundances for Li, O, Na, Mg, Si, Ca, Ti, Ni, Rb, Y, and Ba for the evolved stars. The chemical composition of the cluster is consistent with that of the Galactic thin disc. An overabundance of Ba is found, supporting the enhanced $s$-process. NGC 3105 has a low metallicity for its Galactocentric distance, comparable to typical LMC stars. It is a valuable spiral tracer in a very distant region of the Carina-Sagittarius spiral arm, a poorly known part of the Galaxy. As one of the few Galactic clusters containing blue, yellow and red supergiants, it is massive enough to serve as a testbed for theoretical evolutionary models close to the boundary between intermediate and high-mass stars.

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.]

We present an 80-d long uninterrupted high-cadence K2 light curve of the B1Iab supergiant rho Leo (HD 91316), deduced with the method of halo photometry. This light curve reveals a dominant frequency of $f_{\rmrot}=0.0373$d$^{-1}$ and its harmonics. This dominant frequency corresponds with a rotation period of 26.8d and is subject to amplitude and phase modulation. The K2 photometry additionally reveals multiperiodic low-frequency variability ($<1.5 $d$^{-1}$) and is in full agreement with low-cadence high-resolution spectroscopy assembled during 1800 days. The spectroscopy reveals rotational modulation by a dynamic aspherical wind with an amplitude of about 20km s$^{-1}$ in the H$\alpha$ line, as well as photospheric velocity variations of a few km s$^{-1}$ at frequencies in the range 0.2 to 0.6 d$^{-1}$ in the SiIII 4567\AA\ line. Given the large macroturbulence needed to explain the spectral line broadening of the star, we interpret the detected photospheric velocity as due to travelling super-inertial low-degree large-scale gravity waves with dominant tangential amplitudes and discuss why $\rho$~Leo is an excellent target to study how the observed photospheric variability propagates into the wind.

We introduce VLT-MUSE observations of the central 2$'\times2'$ (30$\times$30 pc) of the Tarantula Nebula in the Large Magellanic Cloud. The observations provide an unprecedented spectroscopic census of the massive stars and ionised gas in the vicinity of R136, the young, dense star cluster located in NGC 2070, at the heart of the richest star-forming region in the Local Group. Spectrophotometry and radial-velocity estimates of the nebular gas (superimposed on the stellar spectra) are provided for 2255 point sources extracted from the MUSE datacubes, and we present estimates of stellar radial velocities for 270 early-type stars (finding an average systemic velocity of 271$\pm$41 km/s). We present an extinction map constructed from the nebular Balmer lines, with electron densities and temperatures estimated from intensity ratios of the [SII], [NII], and [SIII] lines. The interstellar medium, as traced by H$\alpha$ and [NII] $\lambda$6583, provides new insights in regions where stars are probably forming. The gas kinematics are complex, but with a clear bi-modal, blue- and red-shifted distribution compared to the systemic velocity of the gas centred on R136. Interesting point-like sources are also seen in the eastern cavity, western shell, and around R136; these might be related to phenomena such as runaway stars, jets, formation of new stars, or the interaction of the gas with the population of Wolf--Rayet stars. Closer inspection of the core reveals red-shifted material surrounding the strongest X-ray sources, although we are unable to investigate the kinematics in detail as the stars are spatially unresolved in the MUSE data. Further papers in this series will discuss the detailed stellar content of NGC 2070 and its integrated stellar and nebular properties.

Automatic classification methods applied to sky surveys have revolutionized the astronomical target selection process. Most surveys generate a vast amount of time series, or \quotes{lightcurves}, that represent the brightness variability of stellar objects in time. Unfortunately, lightcurves' observations take several years to be completed, producing truncated time series that generally remain without the application of automatic classifiers until they are finished. This happens because state of the art methods rely on a variety of statistical descriptors or features that present an increasing degree of dispersion when the number of observations decreases, which reduces their precision. In this paper we propose a novel method that increases the performance of automatic classifiers of variable stars by incorporating the deviations that scarcity of observations produces. Our method uses Gaussian Process Regression to form a probabilistic model of each lightcurve's observations. Then, based on this model, bootstrapped samples of the time series features are generated. Finally a bagging approach is used to improve the overall performance of the classification. We perform tests on the MACHO and OGLE catalogs, results show that our method classifies effectively some variability classes using a small fraction of the original observations. For example, we found that RR Lyrae stars can be classified with around 80\% of accuracy just by observing the first 5\% of the whole lightcurves' observations in MACHO and OGLE catalogs. We believe these results prove that, when studying lightcurves, it is important to consider the features' error and how the measurement process impacts it.

This paper presents a detailed analysis of the light and radial velocity curves of the semi-detached eclipsing binary system OGLE-LMC-ECL-09937. The system is composed of a hot, massive and luminous primary star of a late-O spectral type, and a more evolved, but less massive and luminous secondary, implying an Algol-type system that underwent a mass transfer episode. We derive masses of 21.04 +/- 0.34 M_Sun and 7.61 +/- 0.09 M_Sun and radii of 9.93 +/- 0.06 R_Sun and 9.18 +/- 0.04 R_Sun, for the primary and the secondary component, respectively, which make it the most massive known Algol-type system with masses and radii of the components measured with <2% accuracy. Consequently, the parameters of OGLE-LMC-ECL-09937 provide an important contribution to the sparsely populated high-mass end of the stellar mass distribution, and an interesting object for stellar evolution studies, being a possible progenitor of a binary system composed of two neutron stars.

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.

We provide an overview of Science Verification MUSE observations of NGC 2070, the central region of the Tarantula Nebula in the Large Magellanic Cloud. Integral-field spectroscopy of the central 2' x 2' region provides the first complete spectroscopic census of its massive star content, nebular conditions and kinematics. The star-formation surface density of NGC 2070 is reminiscent of the intense star-forming knots of high-redshift galaxies, with nebular conditions similar to low-redshift Green Pea galaxies, some of which are Lyman continuum leakers. Uniquely, MUSE permits the star-formation history of NGC 2070 to be studied from both spatially-resolved and integrated-light spectroscopy.

We present the analysis of near infrared (NIR), adaptive optics (AO) Subaru and archived HST imaging data of a region near the northern middle lobe (NML) of the Centaurus A (Cen A) jet, at a distance of $\sim15$ kpc north-east (NE) from the center of NGC5128. Low-pass filtering of the NIR images reveals strong -- $>3\sigma$ above the background mean -- signal at the expected position of the brightest star in the equivalent HST field. Statistical analysis of the NIR background noise suggests that the probability to observe $>3\sigma$ signal at the same position, in three independent measurements due to stochastic background fluctuations alone is negligible ($\leq10^{-7}\%$) and, therefore, that this signal should reflect the detection of the NIR counterparts of the brightest HST star. An extensive photometric analysis of this star yields $V-I$, visual-NIR, and NIR colors expected from a yellow supergiant (YSG) with an estimated age $\sim10^{+4}_{-3}$ Myr. Furthermore, the second and third brighter HST stars are, likely, also supergiants in Cen A, with estimated ages $\sim16^{+6}_{-3}$ Myr and $\sim25^{+15}_{-9}$ Myr, respectively. The ages of these three supergiants are in good agreement with the ages of the young massive stars that were previously found in the vicinity and are thought to have formed during the later phases of the jet-HI cloud interaction that appears to drive the star formation (SF) in the region for the past $\sim100$ Myr.

When a massive star in a binary system explodes as a supernova its companion star may be polluted with heavy elements from the supernova ejecta. Such a pollution had been detected in a handful of post-supernova binaries (Gonzalez Hernandez et al. 2011), but none of them is associated with a supernova remnant. We report the discovery of a solar-type star in a close, eccentric binary system with a neutron star within the young Galactic supernova remnant RCW 86. Our discovery implies that the supernova progenitor was a moving star, which exploded near the edge of its wind bubble and lost most of its initial mass due to common-envelope evolution shortly before core collapse. We find that the solar-type star is strongly polluted with calcium and other elements, which places the explosion within the class of calcium-rich supernovae -- faint and fast transients (Filippenko et al. 2003; Kasliwal et al. 2012}, whose origin is strongly debated (Kawabata et al. 2010; Waldman et al. 2011), and provides the first observational evidence that supernovae of this type can arise from core-collapse explosions (Kawabata et al. 2010; Moriya et al. 2010).

NGC 6067 is a young open cluster hosting the largest population of evolved stars among known Milky Way clusters in the 50-150 Ma age range. It thus represents the best laboratory in our Galaxy to constrain the evolutionary tracks of 5-7 M$_{\odot}$ stars. We have used high-resolution spectra of a large sample of bright cluster members (45), combined with archival photometry, to obtain accurate parameters for the cluster as well as stellar atmospheric parameters. We derive a distance of 1.78$\pm$0.12 kpc, an age of 90$\pm$20 Ma and a tidal radius of 14.8$^{6.8}_{3.2}$ arcmin. We estimate an initial mass above 5700 M$_{\odot}$, for a present-day evolved population of two Cepheids, two A supergiants and 12 red giants with masses $\approx$6 M$_{\odot}$. We also determine chemical abundances of Li, O, Na, Mg, Si, Ca, Ti, Ni, Rb, Y and Ba for the red clump stars. We find a supersolar metallicity, [Fe/H]=+0.19$\pm$0.05, and a homogeneus chemical composition, consistent with the Galactic metallicity gradient. The presence of a Li-rich red giant, star 276 with A(Li)=2.41, is also detected. An over-abundance of Ba is found, supporting the enhanced $s$-process. The ratio of yellow to red giants is much smaller than one, in agreement with models with moderate overshooting, but the properties of the cluster Cepheids do not seem consistent with current Padova models for supersolar metallicity.

The double distance and metallicity frontier marked by the SMC has been finally broken with the aid of powerful multi-object spectrographs installed at 8-10m class telescopes. VLT, GTC and Keck have enabled studies of massive stars in dwarf irregular galaxies of the Local Group with poorer metal content than the SMC. The community is working to test the predictions of evolutionary models in the low-metallicity regime, set the new standard for the metal-poor high-redshift Universe, and test the extrapolation of the physics of massive stars to environments of decreasing metallicity. In this paper, we review current knowledge on this topic.

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

HD 164492C is a spectroscopic triple stellar system that has been recently detected to possess a strong magnetic field. We have obtained high-resolution spectroscopic and spectropolarimetric observations over a timespan of two years and derived physical, chemical, and magnetic properties for this object. The system is formed by an eccentric close spectroscopic binary (Ca1-Ca2) with a period of 12.5 days, and a massive tertiary Cb. We calculated the orbital parameters of the close pair, reconstructed the spectra of the three components, and determined atmospheric parameters and chemical abundances by spectral synthesis. From spectropolarimetric observations, multi-epoch measurements of the longitudinal magnetic fields were obtained. The magnetic field is strongly variable on timescales of a few days, with a most probable period in the range of 1.4-1.6 days. Star Cb with Teff~25000 K is the apparently fastest rotator and the most massive star of this triple system and has anomalous chemical abundances with a marked overabundance of helium, 0.35+/-0.04 by number. We identified this star as being responsible for the observed magnetic field, although the presence of magnetic fields in the components of the Ca pair cannot be excluded. Star Ca1 with a temperature of about 24000 K presents a normal chemical pattern, while the least massive star Ca2 is a mid-B type star (Teff~15000 K) with an overabundance of silicon. The obtained stellar parameters of the system components suggest a distance of 1.5 kpc and an age of 10-15 Myr.

We report the results of spectroscopic observations and numerical modelling of the H II region IRAS 18153-1651. Our study was motivated by the discovery of an optical arc and two main-sequence stars of spectral type B1 and B3 near the centre of IRAS 18153-1651. We interpret the arc as the edge of the wind bubble (blown by the B1 star), whose brightness is enhanced by the interaction with a photoevaporation flow from a nearby molecular cloud. This interpretation implies that we deal with a unique case of a young massive star (the most massive member of a recently formed low-mass star cluster) caught just tens of thousands of years after its stellar wind has begun to blow a bubble into the surrounding dense medium. Our two-dimensional, radiation-hydrodynamics simulations of the wind bubble and the H II region around the B1 star provide a reasonable match to observations, both in terms of morphology and absolute brightness of the optical and mid-infrared emission, and verify the young age of IRAS 18153-1651. Taken together our results strongly suggest that we have revealed the first example of a wind bubble blown by a main-sequence B star.

The origin and evolution of magnetism in OB stars is far from being well understood. With approximately 70 magnetic OB stars known, any new object with unusual characteristics may turn out to be a key piece of the puzzle. We report the detection of an exceptionally strong magnetic field in the He-strong B2IV star CPD-62 2124. Spectropolarimetric FORS2 and HARPSpol observations were analysed by two independent teams and procedures, concluding on a strong longitudinal magnetic field of approximately 5.2 kG. The quantitative characterisation of the stellar atmosphere yields an effective temperature of 23650$\pm$250 K, a surface gravity of 3.95$\pm$0.10 dex and a surface helium fraction of 0.35$\pm$0.02 by number. The metal composition is in agreement with the cosmic abundance standard, except for Mg, Si and S, which are slightly non-solar. The strong and broad ($\sim$300 km/s) disc-like emission displayed by the H$\alpha$ line suggests a centrifugal magnetosphere supported by the strong magnetic field. Our results imply that CPD-62 2124 is an early B-type star hosting one of the strongest magnetic fields discovered to date, and one of the most evolved He-strong stars known, with a fractional main-sequence lifetime of approximately 0.6.

The "B fields in OB stars" (BOB) collaboration is based on an ESO Large Programme, to study the occurrence rate, properties, and ultimately the origin of magnetic fields in massive stars. In the framework of this programme, we carried out low-resolution spectropolarimetric observations of a large sample of massive stars using FORS2 installed at the ESO VLT 8-m telescope. We determined the magnetic field values with two completely independent reduction and analysis pipelines. Our in-depth study of the magnetic field measurements shows that differences between our two pipelines are usually well within 3sigma errors. From the 32 observations of 28 OB stars, we were able to monitor the magnetic fields in CPD-57 3509 and HD164492C, confirm the magnetic field in HD54879, and detect a magnetic field in CPD-62 2124. We obtain a magnetic field detection rate of 6+-3% for the full sample of 69 OB stars observed with FORS2 within the BOB programme. For the pre-selected objects with a v sin i below 60 km/s, we obtain a magnetic field detection rate of 5+-5%. We also discuss X-ray properties and multiplicity of the objects in our FORS2 sample with respect to the magnetic field detections.

In an era of large spectroscopic surveys of stars and big data, sophisticated statistical methods become more and more important in order to infer fundamental stellar parameters such as mass and age. Bayesian techniques are powerful methods because they can match all available observables simultaneously to stellar models while taking prior knowledge properly into account. However, in most cases it is assumed that observables are uncorrelated which is generally not the case. Here, we include correlations in the Bayesian code BONNSAI by incorporating the covariance matrix in the likelihood function. We derive a parametrisation of the covariance matrix that, in addition to classical uncertainties, only requires the specification of a correlation parameter that describes how observables co-vary. Our correlation parameter depends purely on the method with which observables have been determined and can be analytically derived in some cases. This approach therefore has the advantage that correlations can be accounted for even if information for them are not available in specific cases but are known in general. Because the new likelihood model is a better approximation of the data, the reliability and robustness of the inferred parameters are improved. We find that neglecting correlations biases the most likely values of inferred stellar parameters and affects the precision with which these parameters can be determined. For example, we apply our technique to massive OB stars, but emphasise that it is valid for any type of stars. For effective temperatures and surface gravities determined from atmosphere modelling, we find that masses can be underestimated on average by $0.5\sigma$ and mass uncertainties overestimated by a factor of about 2 when neglecting correlations. At the same time, the age precisions are underestimated over a wide range of stellar parameters. [abridged]

We aim to provide new empirical clues about macroturbulent spectral line broadening in O- and B-type stars to evaluate its physical origin. We use high-resolution spectra of ~430 stars with spectral types in the range O4-B9 (all luminosity classes). We characterize the line-broadening of adequate diagnostic metal lines using a combined FT and GOF technique. We perform a quantitative spectroscopic analysis of the whole sample using automatic tools coupled with a huge grid of FASTWIND models. We also incorporate quantitative information about line asymmetries to our observational description of the characteristics of the line-profiles, and present a comparison of the shape and type of line-profile variability found in a small sample of O stars and B supergiants with still undefined pulsational properties and B main sequence stars with variable line-profiles. We present a homogeneous and statistically significant overview of the (single snapshot) line-broadening properties of stars in the whole O and B star domain. We find empirical evidence of the existence of various types of non-rotational broadening agents acting in the realm of massive stars. Even though all of them could be quoted and quantified as a macroturbulent broadening from a practical point of view, their physical origin can be different. Contrarily to the early- to late-B dwarfs/giants, which present a mixture of cases in terms of line-profile shape and variability, the whole O-type and B supergiant domain (or, roughly speaking, stars with M_ZAMS > 15 M_sol) is fully dominated by stars with a remarkable non-rotational broadening component and very similar profiles (including type of variability). We provide some examples illustrating how this observational dataset can be used to evaluate scenarios aimed at explaining the existence of sources of non-rotational broadening in massive stars.

Macroturbulence, introduced as a fudge to reproduce the width and shape of stellar absorption lines, reflects gas motions in stellar atmospheres. While in cool stars, it is thought to be caused by convection zones immediately beneath the stellar surface, the origin of macroturbulence in hot stars is still under discussion. Recent works established a correlation between the turbulent-to-total pressure ratio inside the envelope of stellar models and the macroturbulent velocities observed in corresponding Galactic stars. To probe this connection further, we evaluated the turbulent pressure that arises in the envelope convective zones of stellar models in the mass range 1-125 Msun based on the mixing-length theory and computed for metallicities of the Large and Small Magellanic Cloud. We find that the turbulent pressure contributions in models with these metallicities located in the hot high-luminosity part of the Hertzsprung-Russel (HR) diagram is lower than in similar models with solar metallicity, whereas the turbulent pressure in low-metallicity models populating the cool part of the HR-diagram is not reduced. Based on our models, we find that the currently available observations of hot massive stars in the Magellanic Clouds appear to support a connection between macroturbulence and the turbulent pressure in stellar envelopes. Multidimensional simulations of sub-surface convection zones and a larger number of high-quality observations are necessary to test this idea more rigorously.

Low resolution (4.5 to 5 Angstroem) spectra of 58 blue supergiant stars distributed over the disk of the Magellanic spiral galaxy NGC 55 in the Sculptor group are analyzed by means of non-LTE techniques to determine stellar temperatures, gravities and metallicities (from iron peak and alpha-elements). A metallicity gradient of -0.22 +/- 0.06$ dex/R_25 is detected. The central metallicity on a logarithmic scale relative to the Sun is [Z] = -0.37 +\- 0.03. A chemical evolution model using the observed distribution of stellar and interstellar medium gas mass column densities reproduces the observed metallicity distribution well and reveals a recent history of strong galactic mass accretion and wind outflows with accretion and mass-loss rates of the order of the star formation rate. There is an indication of spatial inhomogeneity in metallicity. In addition, the relatively high central metallicity of the disk confirms that two extra-planar metal poor HII regions detected in previous work 1.13 to 2.22 kpc above the galactic plane are ionized by massive stars formed in-situ outside the disk. For a sub-sample of supergiants, for which Hubble Space Telescope photometry is available, the flux-weighted gravity--luminosity relationship is used to determine a distance modulus of 26.85 +\- 0.10 mag.

A significant fraction of massive main-sequence stars show strong, large-scale magnetic fields. The origin of these fields, their lifetimes, and their role in shaping the characteristics and evolution of massive stars are currently not well understood. We compile a catalogue of 389 massive main-sequence stars, 61 of which are magnetic, and derive their fundamental parameters and ages. The two samples contain stars brighter than magnitude 9 in the V band and range in mass between 5 and 100 Msun. We find that the fractional main-sequence age distribution of all considered stars follows what is expected for a magnitude limited sample, while that of magnetic stars shows a clear decrease towards the end of the main sequence. This dearth of old magnetic stars is independent of the choice of adopted stellar evolution tracks, and appears to become more prominent when considering only the most massive stars. We show that the decreasing trend in the distribution is significantly stronger than expected from magnetic flux conservation. We also find that binary rejuvenation and magnetic suppression of core convection are unlikely to be responsible for the observed lack of older magnetic massive stars, and conclude that its most probable cause is the decay of the magnetic field, over a time span longer than the stellar lifetime for the lowest considered masses, and shorter for the highest masses. We then investigate the spin-down ages of the slowly rotating magnetic massive stars and find them to exceed the stellar ages by far in many cases. The high fraction of very slowly rotating magnetic stars thus provides an independent argument for a decay of the magnetic fields.

Since its launch, the X-ray and gamma-ray observatory INTEGRAL satellite has revealed a new class of high-mass X-ray binaries (HMXB) displaying fast flares and hosting supergiant companion stars. Optical and infrared (OIR) observations in a multi-wavelength context are essential to understand the nature and evolution of these newly discovered celestial objects. The goal of this multiwavelength study (from ultraviolet to infrared) is to characterise the properties of IGR J16465-4507, to confirm its HMXB nature and that it hosts a supergiant star. We analysed all OIR, photometric and spectroscopic observations taken on this source, carried out at ESO facilities. Using spectroscopic data, we constrained the spectral type of the companion star between B0.5 and B1 Ib, settling the debate on the true nature of this source. We measured a high rotation velocity of v = 320 +/- 8 km/s from fitting absorption and emission lines in a stellar spectral model. We then built a spectral energy distribution from photometric observations to evaluate the origin of the different components radiating at each energy range. We finally show that, having accurately determined the spectral type of the early-B supergiant in IGR J16465-4507, we firmly support its classification as an intermediate supergiant fast X-ray transient (SFXT).

Approximately 10% of massive OBA main-sequence (MS) and pre-MS stars harbour strong, large-scale magnetic fields. At the same time there is a dearth of magnetic stars in close binaries. A process generating strong magnetic fields only in some stars must be responsible with the merging of pre-MS and MS stars being suggested as one such channel. Stars emerging from the coalescence of two MS stars are rejuvenated, appearing younger than they are. They can therefore be identified by comparison with reference clocks. Here we predict the rejuvenation of MS merger products over a wide range of masses and binary configurations calibrated to smoothed-particle-hydrodynamical merger models. We find that the rejuvenation is of the order of the nuclear timescale and is strongest in the lowest-mass mergers and the most evolved binary progenitors with the largest mass ratios. These predictions allow us to put constraints on the binary progenitors of merger products. We show that the magnetic stars HR 2949 and $\tau$ Sco are younger than the potential binary companion HR 2948 and the Upper-Sco association, respectively, making them promising merger candidates. We find that the age discrepancies and the potential binary progenitors of both are consistent with them being rejuvenated merger products, implying that their magnetic fields may originate from this channel. Searching for age discrepancies in magnetic stars is therefore a powerful way to explore which fraction of magnetic stars may have obtained their strong magnetic field in MS mergers and to improve our understanding of magnetism in massive stars and their remnants.

We report the detection of a magnetic field in the helium-strong star CPD-57 3509 (B2 IV), a member of the Galactic open cluster NGC3293, and characterise the star's atmospheric and fundamental parameters. Spectropolarimetric observations with FORS2 and HARPSpol are analysed using two independent approaches to quantify the magnetic field strength. A high-S/N FLAMES/GIRAFFE spectrum is analysed using a hybrid non-LTE model atmosphere technique. Comparison with stellar evolution models constrains the fundamental parameters of the star. We obtain a firm detection of a surface averaged longitudinal magnetic field with a maximum amplitude of about 1 kG. Assuming a dipolar configuration of the magnetic field, this implies a dipolar field strength larger than 3.3 kG. Moreover, the large amplitude and fast variation (within about 1 day) of the longitudinal magnetic field implies that CPD-57 3509 is spinning very fast despite its apparently slow projected rotational velocity. The star should be able to support a centrifugal magnetosphere, yet the spectrum shows no sign of magnetically confined material; in particular, emission in H{\alpha} is not observed. Apparently, the wind is either not strong enough for enough material to accumulate in the magnetosphere to become observable or, alternatively, some leakage process leads to loss of material from the magnetosphere. The quantitative spectroscopic analysis of the star yields an effective temperature and a logarithmic surface gravity of 23750+-250 K and 4.05+-0.10, respectively, and a surface helium fraction of 0.28+-0.02 by number. The surface abundances of C, N, O, Ne, S, and Ar are compatible with the cosmic abundance standard, whereas Mg, Al, Si, and Fe are depleted by about a factor of 2. This abundance pattern can be understood as the consequence of a fractionated stellar wind. CPD-57 3509 is one of the most evolved He-strong stars known.

About 10$\%$ of the massive main sequence stars have recently been found to host a strong, large scale magnetic field. Both, the origin and the evolutionary consequences of these fields are largely unknown. We argue that these fields may be sufficiently strong in the deep interior of the stars to suppress convection near the outer edge of their convective core. We performed parametrised stellar evolution calculations and assumed a reduced size of the convective core for stars in the mass range 16 M$_{\odot}$ to 28 M$_{\odot}$ from the zero age main sequence until core carbon depletion. We find that such models avoid the coolest part of the main sequence band, which is usually filled by evolutionary models that include convective core overshooting. Furthermore, our `magnetic' models populate the blue supergiant region during core helium burning, i.e., the post-main sequence gap left by ordinary single star models, and some of them end their life in a position near that of the progenitor of Supernova 1987A in the HR diagram. Further effects include a strongly reduced luminosity during the red supergiant stage, and downward shift of the limiting initial mass for white dwarf and neutron star formation.

The young open cluster Dolidze 25, in the direction of the Galactic Anticentre, has been attributed a very low metallicity, with typical abundances between $-0.5$ and $-0.7$ dex below solar. We intend to derive accurate cluster parameters and accurate stellar abundances for some of its members. We have obtained a large sample of intermediate- and high-resolution spectra for stars in and around Dolidze 25. We used the FASTWIND code to generate stellar atmosphere models to fit the observed spectra. We derive stellar parameters for a large number of OB stars in the area, and abundances of oxygen and silicon for a number of stars with spectral types around B0. We measure low abundances in stars of Dolidze 25. For the three stars with spectral types around B0, we find $0.3$ dex (Si) and $0.5$ dex (O) below the values typical in the solar neighbourhood. These values, even though not as low as those given previously, confirm Dolidze 25 and the surrounding H II region Sh2-284 as the most metal-poor star-forming environment known in the Milky Way. We derive a distance $4.5\pm0.3\:$kpc to the cluster ($r_{\textrm{G}}\approx12.3\:$kpc). The cluster cannot be older than $\sim3\:$Myr, and likely is not much younger. One star in its immediate vicinity, sharing the same distance, has Si and O abundances at most $0.15\:$dex below solar. The low abundances measured in Dolidze 25 are compatible with currently accepted values for the slope of the Galactic metallicity gradient, if we take into account that variations of at least $\pm0.15\:$dex are observed at a given radius. The area traditionally identified as Dolidze 25 is only a small part of a much larger star-forming region that comprises the whole dust shell associated with Sh2-284 and very likely several other smaller H II regions in its vicinity.

In this paper, we present the FATS (Feature Analysis for Time Series) library. FATS is a Python library which facilitates and standardizes feature extraction for time series data. In particular, we focus on one application: feature extraction for astronomical light curve data, although the library is generalizable for other uses. We detail the methods and features implemented for light curve analysis, and present examples for its usage.

The massive red supergiant (RSG) W26 in Westerlund 1 is one of a growing number of RSGs shown to have winds that are ionized from the outside in. The fate of this dense wind material is important for models of second generation star formation in massive star clusters. Mackey et al. (2014) showed that external photoionization can stall the wind of RSGs and accumulate mass in a dense static shell. We use 1D R-HD simulations of an externally photoionized wind to predict the Halpha and [NII] emission arising from photoionized winds both with and without a dense shell. We analyse spectra of the Halpha and [NII] emission in the environment around W26 and compare them with predicted synthetic emission. Simulations of slow winds that are decelerated into a dense shell show strongly limb-brightened line emission, with line radial velocities that are independent of the wind speed. Faster winds (>22 km/s) do not form a dense shell, have less limb-brightening, and the line radial velocity is a good tracer of the wind speed. The brightness of the [NII] and Halpha lines as a function of distance from W26 agrees reasonably well with observations when only the line flux is considered. The radial velocity disagrees, however: the brightest observed emission is blueshifted by ~25 km/s relative to the radial velocity of the star, whereas a spherically symmetric wind has the brightest emission at zero radial velocity. Our results show that the bright nebula surrounding W26 must be asymmetric; we suggest it is confined by external ram pressure from the wind of the nearby supergiant W9. We obtain a lower limit on the nitrogen abundance within the nebula of 2.35 times solar. The line ratio strongly favours photoionization over shock ionization, and so even if the observed nebula is pressure confined there should still be an ionization front and a photoionization-confined shell closer to the star.

We present a study of the properties of the O-type, massive eclipsing binary 2MASS J13130841-6239275 located in the outskirts of the Danks 2 cluster in the G305 star-forming complex, using near-infrared spectroscopy from VLT/ISAAC. We derive the masses and radii to be 24.5$\pm$0.9 M$_{\odot}$ and 9.2$\pm$0.1 R$_{\odot}$ for the primary and 21.7$\pm$0.8 M$_{\odot}$ and 8.7$\pm$0.1 R$_{\odot}$ for the secondary component. In addition, we evaluate the sensitivity of our parameters to the choice of the spectral features used to determine the radial velocities. Both components appear to be main-sequence O6.5$-$O7 type stars at an age of $\sim$5 Myr, which is in agreement with the age of the cluster. A high visual extinction of A$_{5495}$=11.9$\pm$0.1 mag is reported, which is likely attributed to the cold molecular gas contaminating the north-east region of the cluster. By fitting the spectral energy distribution of the system to the available $BVI_{c}JHK_{s}$ photometry, we determine a distance to the system of 3.52$\pm$0.08 kpc with a precision of 2$\%$, which is the most well-determined distance to the Danks 2 cluster and the host complex reported in the literature.

Within the context of the collaboration "B fields in OB stars (BOB)", we used the FORS2 low-resolution spectropolarimeter to search for a magnetic field in 50 massive stars, including two reference magnetic massive stars. Because of the many controversies of magnetic field detections obtained with the FORS instruments, we derived the magnetic field values with two completely independent reduction and analysis pipelines. We compare and discuss the results obtained from the two pipelines. We obtained a general good agreement, indicating that most of the discrepancies on magnetic field detections reported in the literature are caused by the interpretation of the significance of the results (i.e., 3-4 sigma detections considered as genuine, or not), instead of by significant differences in the derived magnetic field values. By combining our results with past FORS1 measurements of HD46328, we improve the estimate of the stellar rotation period, obtaining P = 2.17950+/-0.00009 days. For HD125823, our FORS2 measurements do not fit the available magnetic field model, based on magnetic field values obtained 30 years ago. We repeatedly detect a magnetic field for the O9.7V star HD54879, the HD164492C massive binary, and the He-rich star CPD -57 3509. We obtain a magnetic field detection rate of 6+/-4%, while by considering only the apparently slow rotators we derive a detection rate of 8+/-5%, both comparable with what was previously reported by other similar surveys. We are left with the intriguing result that, although the large majority of magnetic massive stars is rotating slowly, our detection rate is not a strong function of the stellar rotational velocity.

The major mass fraction of the envelope of hot luminous stars is radiatively stable. However, the partial ionisation of hydrogen, helium and iron gives rise to extended sub-surface convection zones in all of them. In this work, we investigate the effect of the pressure induced by the turbulent motion in these zones based on the mixing length theory, and search for observable consequences. We find that the turbulent pressure fraction can amount up to ~5% in OB supergiants, and to ~30% in cooler supergiants. The resulting structural changes are, however, not significantly affecting the evolutionary tracks compared to previous calculations. Instead, a comparison of macroturbulent velocities derived from high quality spectra of OB stars with the turbulent pressure fraction obtained in corresponding stellar models reveals a strong correlation of these two quantities. We discuss a possible physical connection, and conclude that turbulent pressure fluctuations may drive high-order oscillations, which - as conjectured earlier - manifest themselves as macroturbulence in the photospheres of hot luminous stars.

The number of magnetic stars detected among massive stars is small; nevertheless, the role played by the magnetic field in stellar evolution cannot be disregarded. Links between line profile variability, enhancements/depletions of surface chemical abundances, and magnetic fields have been identified for low-mass B-stars, but for the O-type domain this is almost unexplored. Based on FORS2 and HARPS spectropolarimetric data, we present the first detection of a magnetic field in HD54879, a single slowly rotating O9.7 V star. Using two independent and different techniques we obtained the firm detection of a surface average longitudinal magnetic field with a maximum amplitude of about 600 G, in modulus. A quantitative spectroscopic analysis of the star with the stellar atmosphere code FASTWIND results in an effective temperature and a surface gravity of 33000$\pm1000$ K and 4.0$\pm0.1$ dex. The abundances of carbon, nitrogen, oxygen, silicon, and magnesium are found to be slightly lower than solar, but compatible within the errors. We investigate line-profile variability in HD54879 by complementing our spectra with spectroscopic data from other recent OB-star surveys. The photospheric lines remain constant in shape between 2009 and 2014, although H$\alpha$ shows a variable emission. The H$\alpha$ emission is too strong for a standard O9.7 V and is probably linked to the magnetic field and the presence of circumstellar material. Its normal chemical composition and the absence of photospheric line profile variations make HD54879 the most strongly magnetic, non-variable single O-star detected to date.

Massive B-type stars with strong magnetic fields and fast rotation are very rare and provide a mystery for theories of both star formation and magnetic field evolution. Only two such stars, called sigma Ori E analogs, were previously known. Recently, a team involved in APOGEE, one of the Sloan Digital Sky Survey III programs, announced the discovery of two additional rigidly rotating magnetosphere stars, HD23478 and HD345439. The presence of magnetic fields in these newly discovered sigma Ori E analogs was not investigated in the past. In the framework of our ESO Large Programme, and one normal ESO programme, we carried out low-resolution FORS2 spectropolarimetric observations of HD23478 and HD345439. From the measurements using hydrogen lines, we discover a rather strong longitudinal magnetic field of the order of up to 1.5kG in HD23478, and up to 1.3kG using the entire spectrum. The analysis of HD345439 using four subsequent spectropolarimetric subexposures does not reveal the presence of a magnetic field at a significance level of 3sigma. On the other hand, the inspection of individual subexposures indicates that HD345439 may host a strong magnetic field, rapidly varying over 88 minutes. A hint at the fast rotation of HD345439 is also given by the behaviour of several metallic and He I lines in the low-resolution FORS2 spectra, showing profile variations already on such a short time scale.

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.

Within the context of the "B fields in OB stars (BOB)" collaboration, we used the HARPSpol spectropolarimeter to observe the early B-type stars beta CMa (HD44743; B1 II/III) and epsilon CMa (HD52089; B1.5 II). For both stars, we consistently detected the signature of a weak (<30 G in absolute value) longitudinal magnetic field. We determined the physical parameters of both stars and characterise their X-ray spectrum. For beta CMa, our mode identification analysis led to determining a rotation period of 13.6+/-1.2 days and of an inclination angle of the rotation axis of 57.6+/-1.7 degrees, with respect to the line of sight. On the basis of these measurements and assuming a dipolar field geometry, we derived a best fitting obliquity of ~22 degrees and a dipolar magnetic field strength (Bd) of ~100 G (60<Bd<230 G within 1 sigma), below what is typically found for other magnetic massive stars. For epsilon CMa we could only determine a lower limit on the dipolar magnetic field strength of 13 G. For this star, we determine that the rotation period ranges between 1.3 and 24 days. Both stars are expected to have a dynamical magnetosphere. We also conclude that both stars are most likely core hydrogen burning and that they have spent more than 2/3 of their main sequence lifetime. A histogram of the distribution of the dipolar magnetic field strength for the magnetic massive stars known to date does not show the magnetic field "desert" observed instead for intermediate-mass stars. The biases involved in the detection of (weak) magnetic fields in massive stars with the currently available instrumentation and techniques imply that weak fields might be more common than currently observed. Our results show that, if present, even relatively weak magnetic fields are detectable in massive stars and that more observational effort is probably still needed to properly access the magnetic field incidence.

The distribution of stars in the Hertzsprung-Russell diagram narrates their evolutionary history and directly assesses their properties. Placing stars in this diagram however requires the knowledge of their distances and interstellar extinctions, which are often poorly known for Galactic stars. The spectroscopic Hertzsprung-Russell diagram (sHRD) tells similar evolutionary tales, but is independent of distance and extinction measurements. Based on spectroscopically derived effective temperatures and gravities of almost 600 stars, we derive for the first time the observational distribution of Galactic massive stars in the sHRD. While biases and statistical limitations in the data prevent detailed quantitative conclusions at this time, we see several clear qualitative trends. By comparing the observational sHRD with different state-of-the-art stellar evolutionary predictions, we conclude that convective core overshooting may be mass-dependent and, at high mass ($\geq 15\,M_\odot$), stronger than previously thought. Furthermore, we find evidence for an empirical upper limit in the sHRD for stars with $T_{\rm{eff}}$ between 10000 and 32000 K and, a strikingly large number of objects below this line. This over-density may be due to inflation expanding envelopes in massive main-sequence stars near the Eddington limit.

We are now in an era of large spectroscopic surveys of OB-type stars. Quantitative spectroscopic analysis of these modern datasets is enabling us to review the physical properties of blue massive stars with robust samples, not only revisiting mean properties and general trends, but also incorporating information about the effects of second-order parameters. We investigate the spectral type -- effective temperature (SpT - Teff) calibration for O-type dwarfs, and its claimed dependence on metallicity, using statistically-meaningful samples of stars extracted from the IACOB and VFTS surveys. We perform a homogeneous differential spectroscopic analysis of 33 Galactic and 53 LMC O~dwarfs (spanning spectral types of O4 -- O9.7) using the IACOB-GBAT package, a $\chi^2$-fitting algorithm based on a large pre-computed grid of FASTWIND models, and standard techniques for the hydrogen/helium analysis of O-type stars. We compare the estimated effective temperatures and gravities as a function of (internally consistent) spectral classifications. While the general trend is that the temperature of a star increases with earlier spectral types and decreasing metallicity, we show that the large range of gravities found for O-type dwarfs -- spaning up to 0.45-0.50 dex in some spectral bins -- plays a critical role on the dependence of the effective temperature calibrations as a function of spectral type and metallicity. This result warns us about the use of SpT - Teff calibrations for O-dwarfs which ignore the effects of gravity, and highlights the risks of employing calibrations based on small samples. The effects of this scatter in gravities (evolutionary status) for O-type dwarfs should be included in future recipes which employ SpT - Teff calibrations.

The B fields in OB stars (BOB) survey is an ESO large programme collecting spectropolarimetric observations for a large number of early-type stars in order to study the occurrence rate, properties, and ultimately the origin of magnetic fields in massive stars. As of July 2014, a total of 98 objects were observed over 20 nights with FORS2 and HARPSpol. Our preliminary results indicate that the fraction of magnetic OB stars with an organised, detectable field is low. This conclusion, now independently reached by two different surveys, has profound implications for any theoretical model attempting to explain the field formation in these objects. We discuss in this contribution some important issues addressed by our observations (e.g., the lower bound of the field strength) and the discovery of some remarkable objects.

BD+6073 is the optical counterpart of the X-ray source IGR J00370+6122, a probable accretion-powered X-ray pulsar. The X-ray light curve of this binary system shows clear periodicity at 15.7 d, which has been interpreted as repeated outbursts around the periastron of an eccentric orbit. We obtained high-resolution spectra of BD+6073 at different epochs. We used the FASTWind code to generate a stellar atmosphere model to fit the observed spectrum and obtain physical magnitudes. The synthetic spectrum was used as a template for cross-correlation with the observed spectra to measure radial velocities. The radial velocity curve provided an orbital solution for the system. We have also analysed the RXTE/ASM and Swift/BAT light curves to confirm the stability of the periodicity. BD +6073 is a BN0.7 Ib low-luminosity supergiant located at an approximate distance of 3.1 kpc, in the CasOB4 association. We derive Teff=24000 K and log gc=3.0, and chemical abundances consistent with a moderately high level of evolution. The spectroscopic and evolutionary masses are consistent at the 1 sigma level with a mass of 15 solar masses. The recurrence time of the X-ray flares is the orbital period of the system. The NS is in a high eccentricity (e=0.56) orbit, and the X-ray emission is strongly peaked around orbital phase 0.2, though the observations are consistent with some level of X-ray activity happening at all orbital phases. The X-ray behaviour of IGR J00370+6122 is reminiscent of intermediate SFXTs, though its peak luminosity is rather low. The orbit is somewhat wider than those of classical persistent supergiant X-ray binaries, which, combined with the low luminosity of the mass donor, explains the low X-ray luminosity. IGR J00370+6122 will likely evolve towards a persistent supergiant system, highlighting the evolutionary connection between different classes of wind-accreting X-ray sources.

Recent magnetic field surveys in O- and B-type stars revealed that about 10% of the core-hydrogen-burning massive stars host large-scale magnetic fields. The physical origin of these fields is highly debated. To identify and model the physical processes responsible for the generation of magnetic fields in massive stars, it is important to establish whether magnetic massive stars are found in very young star-forming regions or whether they are formed in close interacting binary systems. In the framework of our ESO Large Program, we carried out low-resolution spectropolarimetric observations with FORS2 in 2013 April of the three most massive central stars in the Trifid nebula, HD164492A, HD164492C, and HD164492D. These observations indicated a strong longitudinal magnetic field of about 500-600G in the poorly studied component HD164492C. To confirm this detection, we used HARPS in spectropolarimetric mode on two consecutive nights in 2013 June. Our HARPS observations confirmed the longitudinal magnetic field in HD164492C. Furthermore, the HARPS observations revealed that HD164492C cannot be considered as a single star as it possesses one or two companions. The spectral appearance indicates that the primary is most likely of spectral type B1-B1.5V. Since in both observing nights most spectral lines appear blended, it is currently unclear which components are magnetic. Long-term monitoring using high-resolution spectropolarimetry is necessary to separate the contribution of each component to the magnetic signal. Given the location of the system HD164492C in one of the youngest star formation regions, this system can be considered as a Rosetta Stone for our understanding of the origin of magnetic fields in massive stars.

Star clusters are known as superb tools for understanding stellar evolution. In a quest for understanding the physical origin of magnetism and chemical peculiarity in about 7% of the massive main-sequence stars, we analysed two of the ten brightest members of the ~10 Myr old Galactic open cluster NGC 2264, the early B-dwarfs HD47887 and HD47777. We find accurate rotation periods of 1.95 and 2.64 days, respectively, from MOST photometry. We obtained ESPaDOnS spectropolarimetric observations, through which we determined stellar parameters, detailed chemical surface abundances, projected rotational velocities, and the inclination angles of the rotation axis. Because we found only small (<5 km/s) radial velocity variations, most likely caused by spots, we can rule out that HD47887 and HD47777 are close binaries. Finally, using the least-squares deconvolution technique, we found that both stars possess a large-scale magnetic field with an average longitudinal field strength of about 400 G. From a simultaneous fit of the stellar parameters we determine the evolutionary masses of HD47887 and HD47777 to be 9.4+/-0.7 M0 and 7.6+/-0.5 M0. Interestingly, HD47777 shows a remarkable helium underabundance, typical of helium-weak chemically peculiar stars, while the abundances of HD47887 are normal, which might imply that diffusion is operating in the lower mass star but not in the slightly more massive one. Furthermore, we argue that the rather slow rotation, as well as the lack of nitrogen enrichment in both stars, can be consistent with both the fossil and the binary hypothesis for the origin of the magnetic field. However, the presence of two magnetic and apparently single stars near the top of the cluster mass-function may speak in favour of the latter.

The hard X-ray source IGR J11215-5952 is a peculiar transient, displaying very short X-ray outbursts every 165 days. We obtained high-resolution spectra of the optical counterpart, HD 306414, at different epochs, spanning a total of three months, before and around the 2007 February outburst with the combined aims of deriving its astrophysical parameters and searching for orbital modulation. We fit model atmospheres generated with the fastwind code to the spectrum. We also cross-correlated each individual spectrum to the best-fit model to derive radial velocities. From its spectral features, we classify HD 306414 as B0.5 Ia. From the model fit, we find Teff = 24 700 K and log g = 2.7, in good agreement with the morphological classification. Using the interstellar lines in its spectrum, we estimate a distance to HD 306414 d > 7 kpc. Assuming this distance, we derive R* = 40 Rsol and Mspect = 30 Msol (consistent, within errors, with Mevol = 38 Msol). Radial velocity changes are not dominated by the orbital motion, and we find an upper limit on the semi-amplitude for the optical component Kopt < 11 +- 6 km/s. Large variations in the depth and shape of photospheric lines suggest the presence of strong pulsations, which may be the main cause of the radial velocity changes. Very significant variations, uncorrelated with those of the photospheric lines are seen in the shape and position of the Halpha emission feature around the time of the X-ray outburst, but large excursions are also observed at other times. HD 306414 is a normal B0.5 Ia supergiant. Its radial velocity curve is dominated by an effect that is different from binary motion, and is most likely stellar pulsations. The data available suggest that the X-ray outbursts are caused by the close passage of the neutron star in a very eccentric orbit, perhaps leading to localised mass outflow. (abridged).

We report the results of optical spectroscopy of a candidate evolved massive star in the Cygnus X region, TYC 3159-6-1, revealed via detection of its curious circumstellar nebula in archival data of the Spitzer Space Telecope. We classify TYC 3159-6-1 as an O9.5-O9.7 Ib star and derive its fundamental parameters by using the stellar atmosphere code FASTWIND. The He and CNO abundances in the photosphere of TYC 3159-6-1 are consistent with the solar abundances, suggesting that the star only recently evolved off the main sequence. Proper motion and radial velocity measurements for TYC 3159-6-1 show that it is a runaway star. We propose that Dolidze 7 is its parent cluster. We discuss the origin of the nebula around TYC 3159-6-1 and suggest that it might be produced in several successive episodes of enhanced mass-loss rate (outbursts) caused by rotation of the star near the critical, \Omega-limit.

{We aim to detect and interpret photometric and spectroscopic variability of the bright CoRoT B-type supergiant target HD\,46769 ($V=5.79$). We also attempt to detect a magnetic field in the target.} {We analyse a 23-day oversampled CoRoT light curve after detrending, as well as spectroscopic follow-up data, by using standard Fourier analysis and Phase Dispersion Minimization methods. We determine the fundamental parameters of the star, as well as its abundances from the most prominent spectral lines. We perform a Monte Carlo analysis of spectropolarimetric data to obtain an upper limit of the polar magnetic field, assumping a dipole field.} {In the CoRoT data, we detect a dominant period of 4.84\,d with an amplitude of 87\,ppm, and some of its (sub-)multiples. Given the shape of the phase-folded light curve and the absence of binary motion, we interpret the dominant variability in terms of rotational modulation, with a rotation period of 9.69\,d. Subtraction of the rotational modulation signal does not reveal any sign of pulsations. Our results are consistent with the absence of variability in the Hipparcos light curve. The spectroscopy leads to a projected rotational velocity of 72$\pm 2$\,km\,s$^{-1}$ and does not reveal periodic variability nor the need to invoke macroturbulent line broadening. No signature of a magnetic field is detected in our data. A field stronger than $\sim 500$\,G at the poles can be excluded, unless the possible non-detected field were more complex than dipolar.} {The absence of pulsations and of macroturbulence of this evolved B-type supergiant is placed into context of instability computations and of observed variability of evolved B-type stars.}

Recent studies have claimed the existence of very massive stars (VMS) up to 300 solar masses in the local Universe. As this finding may represent a paradigm shift for the canonical stellar upper-mass limit of 150 Msun, it is timely to discuss the status of the data, as well as the far-reaching implications of such objects. We held a Joint Discussion at the General Assembly in Beijing to discuss (i) the determination of the current masses of the most massive stars, (ii) the formation of VMS, (iii) their mass loss, and (iv) their evolution and final fate. The prime aim was to reach broad consensus between observers and theorists on how to identify and quantify the dominant physical processes.

We present the Spitzer Extragalactic Representative Volume Survey (SERVS), an 18 square degrees medium-deep survey at 3.6 and 4.5 microns with the post-cryogenic Spitzer Space Telescope to ~2 microJy (AB=23.1) depth of five highly observed astronomical fields (ELAIS-N1, ELAIS-S1, Lockman Hole, Chandra Deep Field South and XMM-LSS). SERVS is designed to enable the study of galaxy evolution as a function of environment from z~5 to the present day, and is the first extragalactic survey both large enough and deep enough to put rare objects such as luminous quasars and galaxy clusters at z>1 into their cosmological context. SERVS is designed to overlap with several key surveys at optical, near- through far-infrared, submillimeter and radio wavelengths to provide an unprecedented view of the formation and evolution of massive galaxies. In this paper, we discuss the SERVS survey design, the data processing flow from image reduction and mosaicing to catalogs, as well as coverage of ancillary data from other surveys in the SERVS fields. We also highlight a variety of early science results from the survey.

The quantitative study of the physical properties and chemical abundances of large samples of massive blue stars at different metallicities is a powerful tool to understand the nature and evolution of these objects. Their analysis beyond the Milky Way is challenging, nonetheless it is doable and the best way to investigate their behavior in different environments. Fulfilling this task in an objective way requires the implementation of automatic analysis techniques that can perform the analyses systematically, minimizing at the same time any possible bias. As part of the ARAUCARIA project we carry out the first quantitative spectroscopic analysis of a sample of 12 B-type supergiants in the galaxy NGC55 at 1.94 Mpc away. By applying the methodology developed in this work, we derive their stellar parameters, chemical abundances and provide a characterization of the present-day metallicity of their host galaxy. Based on the characteristics of the stellar atmosphere/line formation code FASTWIND, we designed and created a grid of models for the analysis of massive blue supergiant stars. Along with this new grid, we implemented a spectral analysis algorithm. Both tools were specially developed to perform fully consistent quantitative spectroscopic analyses of low spectral resolution of B-type supergiants in a fast and objective way. We present the main characteristics of our FASTWIND model grid and perform a number of tests to investigate the reliability of our methodology. The automatic tool is applied afterward to a sample of 12 B-type supergiant stars in NGC55, deriving the stellar parameters and abundances. The results indicate that our stars are part of a young population evolving towards a red supergiant phase. The derived chemical composition hints to an average metallicity similar to the one of the Large Magellanic Cloud, with no indication of a spatial trend across the galaxy.

Luminous blue variables (LBVs) are expected to play an important role in massive stellar evolution as well as being the progenitors of some of the most luminous supernovae known. In this paper we provide a multiwavelength study of the population of (candidate) LBVs identified within M33. New spectra provide an observational baseline of >4yr with respect to published data, which is well suited to identifying LBV outbursts. Multi-epoch optical and mid-IR surveys of M33 further constrain the variability of the sample and permit a search for dusty circumstellar ejecta. Spectroscopic and photometric variability appears common amongst the sample, although in many cases further observations will be required to determine its physical origin. Nevertheless, we report a new outburst of M33 Var C, while the transition of the WNLh star B517 to a cooler B supergiant phase between 1993-2010 confirms an LBV classification. Proof-of-concept quantitative analysis is provided for Romano's star; the results being consistent with the finding that its bolometric luminosity varies during its LBV excursions. The combination of the temperature and luminosity of two stars, the B hypergiant [HS80] 110A and the cool hypergiant B324, appears to be in violation of the empirical Humphreys-Davidson limit. Mid-IR observations demonstrate that a number of candidates appear associated with hot circumstellar dust, although no objects as extreme as Eta Carinae are identified. The combined multiwavelength dataset suggests that the population of LBVs studied is contaminated by stars demonstrating the B[e] phenomenon. Of these, a subset of optically faint, low luminosity stars associated with hot dust are of particular interest since they appear similar to the likely progenitors of SN 2008S and the 2008 NGC300 transient, albeit suffering less intrinsic extinction. [ABRIDGED]

We present a brief progress report of our project aimed at the investigation of the so-called macroturbulent broadening affecting the line-profiles of O and early B-type stars, and speculate about the possible use of this spectroscopic feature as a single snap-shot approach to investigate pulsations in massive stars.

We present the IACOB grid-based automatic tool for the quantitative spectroscopic analysis of O-stars. The tool consists of an extensive grid of FASTWIND models, and a variety of programs implemented in IDL to handle the observations, perform the automatic analysis, and visualize the results. The tool provides a fast and objective way to determine the stellar parameters and the associated uncertainties of large samples of O-type stars within a reasonable computational time.

This paper presents the first distance measurement to the massive, semi-detached, eclipsing binary LMC-SC1-105, located in the LH 81 association of the Large Magellanic Cloud (LMC). Previously determined parameters of the system are combined with new near-infrared photometry and a new temperature analysis to constrain the reddening toward the system, and determine a distance of 50.6+/-1.6 kpc (corresponding to a distance modulus of 18.52+/-0.07 mag), in agreement with previous eclipsing binary measurements. This is the sixth distance measurement to an eclipsing binary in the LMC, although the first to an O-type system. We thus demonstrate the suitability of O-type eclipsing binaries (EBs) as distance indicators. We suggest using bright, early-type EBs to measure distances along different sight lines, as an independent way to map the depth of the LMC and resolve the controversy about its three-dimensional structure.

The important role of metallicity on massive star evolution and the combination of multi-object spectrographs and 10m class telescopes, have lead to numerous systematic studies of massive stars in Local Group galaxies. While color based quests of blue massive stars are relatively successful, they must be confirmed with spectroscopy and usually lead to lists dominated by B-type modest-mass stars. We have developed a friends of friends code to find OB associations in Local Group galaxies, presented in Garcia et al. (2009). One of the key points of the method is the photometric criterion to choose candidate OB stars from the reddening-free Q parameter, that could be easily extended to include from GALEX to near-IR photometry. While not a new idea, one of our code's strong advantages is the automatic determination of evolutionary masses for the members, enabling a quick and more insightful choice of candidates for spectroscopy, and the identification of potential advanced evolutionary stages. We present our work on the very metal-poor irregular IC 1613 (Garcia et al. 2010). The association properties are not only a powerful aid towards finding the most interesting candidate massive stars, but also reveal the galaxy's structure and recent star formation history.

We present the IACOB spectroscopic database, an homogeneous set of high quality, high resolution spectra of Galactic O- and B-type stars obtained with the FIES spectrograph attached to the Nordic Optical Telescope. We also present some results from ongoing projects using the IACOB database.

As part of a long term observational project, we are investigating the macroturbulent broadening in O and B supergiants (Sgs) and its possible connection with spectroscopic variability phenomena and stellar oscillations. We present the first results of our project, namely firm observational evidence for a strong correlation between the extra broadening and photospheric line-profile variations in a sample of 13 Sgs with spectral types ranging from O9.5 to B8.

In this poster we present the results of our analyses of three early massive stars in IC 1613, whose spectra have been observed with VIMOS and analyzed with CMFGEN and FASTWIND. One of the targets resulted a possible LBV and the other two are Of stars with unexpectedly strong winds. The Of stars seem to be strongly contaminated by CNO products. Our preliminary results may represent a challenge for the theory of stellar atmospheres, but they still have to be confirmed by the analysis of more objects and a more complete coverage of the parameter space.

We present the IACOB spectroscopic database, the largest homogeneous database of high-resolution, high signal-to-noise ratio spectra of Northern Galactic OB-type stars compiled up to date. The spectra were obtained with the FIES spectrograph attached to the Nordic Optical Telescope. We briefly summarize the main characeristics and present status of the IACOB, first scientific results, and some future plans for its extension and scientific exploitation.

The spectra of O and B supergiants are known to be affected by a significant form of extra line broadening (usually referred to as macroturbulence) in addition to that produced by stellar rotation. Recent analyses of high resolution spectra have shown that the interpretation of this line broadening as a consequence of large scale turbulent motions would imply highly supersonic velocity fields in photospheric regions, making this scenario quite improbable. Stellar oscillations have been proposed as a likely alternative explanation. As part of a long term observational project, we are investigating the macroturbulent broadening in O and B supergiants and its possible connection with spectroscopic variability phenomena and stellar oscillations. In this letter, we present the first encouraging results of our project, namely firm observational evidence for a strong correlation between the extra broadening and photospheric line-profile variations in a sample of 13 supergiants with spectral types ranging from O9.5 to B8.

We set out to determine the ratio, q(IR), of rest-frame 8-1000um flux, S(IR), to monochromatic radio flux, S(1.4GHz), for galaxies selected at far-IR and radio wavelengths, to search for signs that the ratio evolves with redshift, luminosity or dust temperature, and to identify any far-IR-bright outliers - useful laboratories for exploring why the far-IR/radio correlation is generally so tight when the prevailing theory suggests variations are almost inevitable. We use flux-limited 250-um and 1.4-GHz samples, obtained in GOODS-N using Herschel (HerMES; PEP) and the VLA. We determine bolometric IR output using ten bands spanning 24-1250um, exploiting data from PACS and SPIRE, as well as Spitzer, SCUBA, AzTEC and MAMBO. We also explore the properties of an L(IR)-matched sample, designed to reveal evolution of q(IR) with z, spanning log L(IR) = 11-12 L(sun) and z=0-2, by stacking into the radio and far-IR images. For 1.4-GHz-selected galaxies, we see tentative evidence of a break in the flux ratio, q(IR), at L(1.4GHz) ~ 10^22.7 W/Hz, where AGN are starting to dominate the radio power density, and of weaker correlations with z and T(d). From our 250-um-selected sample we identify a small number of far-IR-bright outliers, and see trends of q(IR) with L(1.4GHz), L(IR), T(d) and z, noting that some of these are inter-related. For our L(IR)-matched sample, there is no evidence that q(IR) changes significantly as we move back into the epoch of galaxy formation: we find q(IR) goes as (1+z)^gamma, where gamma = -0.04 +/- 0.03 at z=0-2; however, discounting the least reliable data at z<0.5 we find gamma = -0.26 +/- 0.07, modest evolution which may be related to the radio background seen by ARCADE2, perhaps driven by <10uJy radio activity amongst ordinary star-forming galaxies at z>1.

The spectrum of O and B Supergiants is known to be affected by an important extra line-broadening (usually called macroturbulence) that adds to stellar rotation. Recent analysis of high resolution spectra has shown that the interpretation of this line-broadening as a consequence of large-scale turbulent motions would imply highly super-sonic velocity fields, making this scenario quite improbable. Stellar oscillations have been proposed as a likely alternative explanation. We present first encouraging results of an observational project aimed at investigating the $macroturbulent$ broadening in O and B Supergiants, and its possible connection with spectroscopic variability phenomena and stellar oscillations: a) all the studied B Supergiants show line profile variations, quantified by means of the first (< v >) and third velocity (< v^3 >) moments of the lines, b) there is a strong correlation between the peak-to-peak amplitudes of the < v > and < v^3 > variability and the size of the extra-broadening.

[abridged] Context: Very few examples of luminous blue variable (LBV) stars or LBV candidates (LBVc) are known, particularly at metallicities below the SMC. The LBV phase is crucial for the evolution of massive stars, and its behavior with metallicity is poorly known. V39 in IC 1613 is a well-known photometric variable, with B-band changes larger than 1mag. over its period. The star, previously proposed to be a projection of a Galactic W Virginis and an IC 1613 red supergiant, shows features that render it a possible LBVc. Method: We investigate mid-resolution blue and red VLT-VIMOS spectra of V39, covering a time span of 40 days, and perform a quantitative analysis of the combined spectrum using the model atmosphere code CMFGEN. Results: We identify strong Balmer and FeII P-Cygni profiles, and a hybrid spectrum resembling a B-A supergiant in the blue and a G-star in the red. No significant Vrad variations are detected, and the spectral changes are small over the photometric period. Our analysis places V39 in the low-luminosity part of the LBV and LBVc region, but it is also consistent with a sgB[e] star. Conclusions: The radial velocity indicates that V39 belongs to IC 1613. The lack of Vrad changes and spectroscopic variations excludes binary scenarios. The features observed are not consistent with a W Virginis star, and this possibility is also discarded. We propose that the star is a B-A LBVc or sgB[e] star surrounded by a thick disk precessing around it. If confirmed, V39 would be the lowest metallicity resolved LBV candidate known to date. Alternatively, it could represent a new transient phase of massive star evolution, an LBV impostor.

Several works have recently shown that there is an important extra line-broadening (usually called $macroturbulence$) affecting the spectra of O and B Supergiants that adds to stellar rotation. So far, the only (very recent) physical explanation for the appearance of $macroturbulence$ relates to oscillations. This is a plausible explanation, but no direct evidence confirming its validity has been presented yet. We recently started an observational project to obtain constraints on the time-scales of variability associated to this extra line-broadening and its possible origin. Our observational strategy consists of the study of a well selected group of O and B stars, for which we obtain time series of high-quality spectra. We present some preliminary results from our first campaign with FIES@NOT2.5m.

Context: Determining the parameters of massive stars is crucial to understand many processes in galaxies and the Universe, since these objects are important sources of ionization, chemical enrichment and momentum. 10m class telescopes enable us to perform detailed quantitative spectroscopic analyses of massive stars in other galaxies, sampling areas of different metallicity. Relating the stars to their environment is crucial to understand the physical processes ruling their formation and evolution. Aims: In preparation for the GTC, our goal is to build a catalogue of massive star candidates in the metal-poor irregular galaxy IC1613 with high astrometric accuracy, apt for the current generation of multi-object spectrographs. A census of OB associations in this galaxy is also needed, to provide important additional information about age and environment of the candidate OB stars. Methods: From INT-WFC observations, we have built an astrometric and photometric catalogue of stars in IC1613. Candidate blue massive stars are preselected from their colors. A friends-of-friends algorithm is developed to find their clustering in the galaxy. While a common physical origin for all the members of the associations cannot be ensured, this is a necessary first step to place candidate OB stars in a population context. Results: We have produced a deep catalogue of targets in IC1613 that covers a large field of view. To achieve high astrometric accuracy a new astrometric procedure is developed for the INT-WFC data. We have also built a catalogue of OB associations in IC1613. We have found that they concentrate in the central regions, specially in the HII bubbles. The study of extinction confirms that it is patchy, with local values of color-excess above the foreground value.

This is the first paper of a series devoted to studying the population of blue massive stars in NGC 55, a galaxy of the Sculptor group at a distance of about 2 Mpc. We have obtained optical (3300-6210A) low-resolution spectra of blue massive stars with VLT-FORS2, which we have classified with the aid of Milky Way and Magellanic Cloud standard stars. We have developed the first census of massive blue stars in NGC 55. A study of stellar radial velocities shows agreement with existing HI rotational velocity curve work and reveals the presence of one object with peculiar velocity. A qualitative study of the stellar metallicity suggests that its global distribution over NGC 55 is close to that of the LMC, as derived from previous studies. We present a catalogue with 164 classifications of blue massive stars in NGC 55. This catalogue is a first and necessary step for the subsequent quantitative study of blue massive stars in NGC 55 with state-of-the-art models.