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.

The multiplicity properties of massive stars are one of the important outstanding issues in stellar evolution. Quantifying the binary statistics of all evolutionary phases is essential to paint a complete picture of how and when massive stars interact with their companions, and to determine the consequences of these interactions. We investigate the multiplicity of an almost complete census of red supergiant stars (RSGs) in NGC 330, a young massive cluster in the SMC. Using a combination of multi-epoch HARPS and MUSE spectroscopy, we estimate radial velocities and assess the kinematic and multiplicity properties of 15 RSGs in NGC 330. Radial velocities are estimated to better than +/-100 m/s for the HARPS data. The line-of-sight velocity dispersion for the cluster is estimated as 3.20 +0.69-0.52 km/s. When virial equilibrium is assumed, the dynamical mass of the cluster is log (M{dyn} /M{sun}) = 5.20+/-0.17, in good agreement with previous upper limits. We detect significant radial velocity variability in our multi-epoch observations and distinguish between variations caused by atmospheric activity and those caused by binarity. The binary fraction of NGC 330 RSGs is estimated by comparisons with simulated observations of systems with a range of input binary fractions. In this way, we account for observational biases and estimate the intrinsic binary fraction for RSGs in NGC 330 as f{RSG} = 0.3+/-0.1 for orbital periods in the range 2.3< log P [days] <4.3, with q>0.1. Using the distribution of the luminosities of the RSG population, we estimate the age of NGC 330 to be 45+/-5 Myr and estimate a red straggler fraction of 50%. We estimate the binary fraction of RSGs in NGC 330 and conclude that it appears to be lower than that of main-sequence massive stars, which is expected because interactions between an RSG and a companion are assumed to effectively strip the RSG envelope.

The orbits of the least chemically enriched stars open a window on the formation of our Galaxy when it was still in its infancy. The common picture is that these low-metallicity stars are distributed as an isotropic, pressure-supported component since these stars were either accreted from the early building blocks of the assembling Milky Way, or were later brought by the accretion of faint dwarf galaxies. Combining the metallicities and radial velocities from the Pristine and LAMOST surveys and Gaia DR2 parallaxes and proper motions for an unprecedented large and unbiased sample of very metal-poor stars at $[Fe/H]\leq-2.5$ we show that this picture is incomplete. This sample shows strong statistical evidence (at the $5.0\sigma$ level) of asymmetry in their kinematics, favouring prograde motion. Moreover, we find that $31\%$ of the stars that currently reside in the disk do not venture outside of the disk plane throughout their orbit. The discovery of this population implies that a significant fraction of stars with iron abundances $[Fe/H]\leq-2.5$ formed within or concurrently with the Milky Way disk and that the history of the disk was quiet enough to allow them to retain their disk-like orbital properties.

A majority of massive stars are part of binary systems, a large fraction of which will inevitably interact during their lives. Binary-interaction products (BiPs), i.e. stars affected by such interaction, are expected to be commonly present in stellar populations. BiPs are thus a crucial ingredient in the understanding of stellar evolution. We aim to identify and characterize a statistically significant sample of BiPs by studying clusters of 10-40 Myr, an age at which binary population models predict the abundance of BiPs to be highest. One example of such a cluster is NGC 330 in the Small Magellanic Cloud. Using MUSE WFM-AO observations of NGC 330, we resolve the dense cluster core for the first time and are able to extract spectra of its entire massive star population. We develop an automated spectral classification scheme based on the equivalent widths of spectral lines in the red part of the spectrum. We characterize the massive star content of the core of NGC 330 which contains more than 200 B stars, 2 O stars, 6 A-type supergiants and 11 red supergiants. We find a lower limit on the Be star fraction of $32 \pm 3\%$ in the whole sample. It increases to at least $46 \pm 10\%$ when only considering stars brighter than V=17 mag. We estimate an age of the cluster core between 35 and 40 Myr and a total cluster mass of $88^{+17}_{-18} \times 10^3 M_{\odot}$. We find that the population in the cluster core is different than the population in the outskirts: while the stellar content in the core appears to be older than the stars in the outskirts, the Be star fraction and the observed binary fraction are significantly higher. Furthermore, we detect several BiP candidates that will be subject of future studies.

The incidence of multiplicity in cool, luminous massive stars is relatively unknown compared to their hotter counterparts. Here we present radial velocity (RV) measurements and investigate the multiplicity properties of red supergiants (RSGs) in the 30~Doradus region of the Large Magellanic Cloud. We provide absolute RV measurements for our sample and estimate line-of-sight velocities for the Hodge 301 and SL 639 clusters, which agree well with those of hot stars in the same clusters. By combining results for the RSGs with those for nearby B-type stars, we estimate systemic velocities and velocity dispersions for the two clusters, obtaining estimates for their dynamical masses of $\log (M_{\rm dyn}/M_{\odot})=$3.8$\pm$0.3 for Hodge 301, and an upper limit of $\log (M_{\rm dyn}/M_{\odot})$<3.1$\pm$0.8 for SL 639, assuming Virial equilibrium. Analysis of the multi-epoch data reveals one RV-variable, potential binary candidate (VFTS744), which is likely a semi-regular variable asymptotic giant branch star. We estimate an upper limit on the observed binary fraction for our sample of 0.3, where we are sensitive to maximum periods for individual objects in the range of 1 to 10 000 days and mass-ratios above 0.3 depending on the data quality. From simulations of the RV measurements from binary systems given the current data we conclude that systems within the parameter range q>0.3, $\log$P[days]<3.5, would be detected by our variability criteria, at the 90% confidence level. The intrinsic binary fraction, accounting for observational biases, is estimated using simulations of binary systems with an empirically defined distribution of parameters where orbital periods are uniformly distributed in the 3.3<$\log$P[days]<4.3 range. A range of intrinsic binary fractions are considered; a binary fraction of 0.3 is found to best reproduce the observed data. [Abridged]

We estimate physical parameters for the late-type massive stars observed as part of the VLT-FLAMES Tarantula Survey (VFTS) in the 30 Doradus region of the Large Magellanic Cloud (LMC). The observational sample comprises 20 candidate red supergiants (RSGs) which are the reddest (($B-V$) $>$ 1 mag) and brightest ($V$ $<$ 16 mag) objects in the VFTS. We use optical and near-IR photometry to estimate their temperatures and luminosities, and introduce the luminosity-age diagram to estimate their ages. We derive physical parameters for our targets, including temperatures from a new calibration of $(J-K_{\rm s})_{0}$ colour for luminous cool stars in the LMC, luminosities from their $J$-band magnitudes (thence radii), and ages from comparisons with state-of-the-art evolutionary models. We show that interstellar extinction is a significant factor for our targets, highlighting the need to take it into account in analysis of the physical parameters of RSGs. We find that some of the candidate RSGs could be massive AGB stars. The apparent ages of the RSGs in the Hodge 301 and SL 639 clusters show a significant spread (12-24 Myr). We also apply our approach to the RSG population of the relatively nearby NGC 2100 cluster, finding a similarly large spread. We argue that the effects of mass-transfer in binaries may lead to more massive and luminous RSGs (which we call `red stragglers') than expected from single-star evolution, and that the true cluster ages correspond to the upper limit of the estimated RSG ages. In this way, the RSGs can serve as a new and potentially reliable age tracer in young star clusters. The corresponding analysis yields ages of 24$^{+5}_{-3}$ Myr for Hodge 301, 22$^{+6}_{-5}$ Myr for SL 639, and 23$^{+4}_{-2}$ Myr for NGC 2100.

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.

The Arches is one of the youngest, densest and most massive clusters in the Galaxy. As such it provides a unique insight into the lifecycle of the most massive stars known and the formation and survival of such stellar aggregates in the extreme conditions of the Galactic Centre. In a previous study we presented an initial stellar census for the Arches and in this work we expand upon this, providing new and revised classifications for ~30% of the 105 spectroscopically identified cluster members as well as distinguishing potential massive runaways. The results of this survey emphasise the homogeneity and apparent co-evality of the Arches and confirm the absence of H-free Wolf-Rayets of WC sub-type and predicted luminosities. The increased depth of our complete dataset also provides significantly better constraints on the main sequence population; with the identification of O9.5 V stars for the first time we now spectroscopically sample stars with initial masses ranging from ~16Msun to >120Msun. Indeed, following from our expanded stellar census we might expect >50 stars within the Arches to have been born with masses >60Msun, while all 105 spectroscopically confirmed cluster members are massive enough to leave relativistic remnants upon their demise. Moreover the well defined observational properties of the main sequence cohort will be critical to the construction of an extinction law appropriate for the Galactic Centre and consequently the quantitative analysis of the Arches population and subsequent determination of the cluster initial mass function.

The Quintuplet is one of the most massive galactic clusters known, but appears to host a diverse stellar population that is difficult to reconcile with an instantaneous formation event. We present HST photometry and VLT spectroscopy in order to improve observational constraints, finding the Quintuplet to be far more homogeneous than previously thought. O7-8 Ia and O9-B0 Ia supergiants form a smooth morphological sequence with a cohort of seven early-B hypergiants and six luminous blue variables and WN9-11h stars, which together comprise the richest population of such stars of any known stellar aggregate. No further H-free Wolf-Rayets were found, resulting in a 13:1 ratio for WC/WN stars. However a small population of late-O hypergiants and WN8-9ha stars was identified, while a subset of the supergiants are unexpectedly faint, suggesting they are both less massive and older than the greater cluster population. We find an impressive coincidence between cluster members preceding the H-free WR phase and the evolutionary predictions for a 60Msun star, suggesting an age of ~3.0-3.6Myr for the Quintuplet. Neither the late-O hypergiants nor the low luminosity supergiants are predicted; we suggest that the former either result from rapid rotators or are the products of binary driven mass-stripping, while the latter may be interlopers. The H-free WRs must have had initial masses >60Msun, but it is difficult to reconcile their observational properties with theoretical expectations. Since the WRs represent an evolutionary phase directly preceding core-collapse, they are are crucial to understanding both this process and the nature of the resultant relativistic remnants, providing unique constraints on the evolution and death of the most massive stars forming in the local, high metallicity Universe (Abridged).

We have carried out a spectroscopic variability survey of some of the most massive stars in the Arches cluster, using K-band observations obtained with SINFONI on the VLT. One target, F2, exhibits substantial changes in radial velocity; in combination with new KMOS and archival SINFONI spectra, its primary component is found to undergo radial velocity variation with a period of 10.483+/-0.002 d and an amplitude of ~350 km/s. A secondary radial velocity curve is also marginally detectable. We reanalyse archival NAOS-CONICA photometric survey data in combination with our radial velocity results to confirm this object as an eclipsing SB2 system, and the first binary identified in the Arches. We model it as consisting of an 82+/-12 M_sun WN8-9h primary and a 60+/-8 M_sun O5-6 Ia+ secondary, and as having a slightly eccentric orbit, implying an evolutionary stage prior to strong binary interaction. As one of four X-ray bright Arches sources previously proposed as colliding-wind massive binaries, it may be only the first of several binaries to be discovered in this cluster, presenting potential challenges to recent models for the Arches' age and composition. It also appears to be one of the most massive binaries detected to date; the primary's calculated initial mass of >~120 M_sun would arguably make this the most massive binary known in the Galaxy.

We present an abundance analysis of seven super-star clusters in the disk of M83. The near-infrared spectra of these clusters are dominated by Red Supergiants, and the spectral similarity in the J-band of such stars at uniform metallicity means that the integrated light from the clusters may be analysed using the same tools as those applied to single stars. Using data from VLT/KMOS we estimate metallicities for each cluster in the sample. We find that the abundance gradient in the inner regions of M83 is flat, with a central metallicity of [Z] = 0.21$\pm$0.11 relative to a Solar value of $Z_\odot$=0.014, which is in excellent agreement with the results from an analysis of luminous hot stars in the same regions. Compiling this latest study with our other recent work, we construct a mass-metallicity relation for nearby galaxies based entirely on the analysis of RSGs. We find excellent agreement with the other stellar-based technique, that of blue supergiants, as well as with temperature-sensitive (`auroral' or `direct') \hii-region studies. Of all the HII-region strong-line calibrations, those which are empirically calibrated to direct-method studies (N2 and O3N2) provide the most consistent results.

We present K-band Multi-Object Spectrograph (KMOS) observations of 18 Red Supergiant (RSG) stars in the Sculptor Group galaxy NGC 55. Radial velocities are calculated and are shown to be in good agreement with previous estimates, confirming the supergiant nature of the targets and providing the first spectroscopically confirmed RSGs in NGC 55. Stellar parameters are estimated for 14 targets using the $J$-band analysis technique, making use of state-of-the-art stellar model atmospheres. The metallicities estimated confirm the low-metallicity nature of NGC 55, in good agreement with previous studies. This study provides an independent estimate of the metallicity gradient of NGC 55, in excellent agreement with recent results published using hot massive stars. In addition, we calculate luminosities of our targets and compare their distribution of effective temperatures and luminosities to other RSGs, in different environments, estimated using the same technique.

We have obtained K-band Multi-Object Spectrograph (KMOS) near-IR spectroscopy for 14 red supergiant stars (RSGs) in the young massive star cluster NGC 2100 in the Large Magellanic Cloud (LMC). Stellar parameters including metallicity are estimated using the J-band analysis technique, which has been rigorously tested in the Local Universe. We find an average metallicity for NGC 2100 of [Z]=$-$0.43$\pm$0.10 dex, in good agreement with estimates from the literature for the LMC. Comparing our results in NGC 2100 with those for a Galactic cluster (at Solar-like metallicity) with a similar mass and age we find no significant difference in the location of RSGs in the Hertzsprung--Russell diagram. We combine the observed KMOS spectra to form a simulated integrated-light cluster spectrum and show that, by analysing this spectrum as a single RSG, the results are consistent with the average properties of the cluster. Radial velocities are estimated for the targets and the dynamical properties are estimated for the first time within this cluster. The data are consistent with a flat velocity dispersion profile, and with an upper limit of 3.9 km/s, at the 95% confidence level, for the velocity dispersion of the cluster. However, the intrinsic velocity dispersion is unresolved and could, therefore, be significantly smaller than the upper limit reported here. An upper limit on the dynamical mass of the cluster is derived as $M_{dyn}$ $\le$ $15.2\times10^{4}M_{\odot}$ assuming virial equilibrium.

We present a direct determination of the stellar metallicity in the close pair galaxy NGC~4038 (D = 20 Mpc) based on the quantitative analysis of moderate resolution KMOS/VLT spectra of three super star clusters (SSCs). The method adopted in our analysis has been developed and optimised to measure accurate metallicities from atomic lines in the $J$-band of single red supergiant (RSG) or RSG-dominated star clusters. Hence, our metallicity measurements are not affected by the biases and poorly understood systematics inherent to strong line HII methods which are routinely applied to massive data sets of galaxies. We find [Z]= +0.07 $\pm$ 0.03 and compare our measurements to HII strong line calibrations. Our abundances and literature data suggest the presence of a flat metallicity gradient, which can be explained as redistribution of metal-rich gas following the strong interaction.

We present a quantitative spectroscopic study of twenty-seven red supergiants in the Sculptor Galaxy NGC 300. J-band spectra were obtained using KMOS on the VLT and studied with state of the art synthetic spectra including NLTE corrections for the strongest diagnostic lines. We report a central metallicity of [Z]= -0.03 +/- 0.05 with a gradient of -0.083 +/- 0.014 [dex/kpc], in agreement with previous studies of blue supergiants and H II-region auroral line measurements. This result marks the first application of the J-band spectroscopic method to a population of individual red supergiant stars beyond the Local Group of galaxies and reveals the great potential of this technique.

Red Supergiants (RSGs) are cool (~4000K), highly luminous stars (L - 10^5 Lsun), and are among the brightest near-infrared (NIR) sources in star-forming galaxies. This makes them powerful probes of the properties of their host galaxies, such as kinematics and chemical abundances. We have developed a technique whereby metallicities of RSGs may be extracted from a narrow spectral window around 1{\mu}m from only moderate resolution data. The method is therefore extremely efficient, allowing stars at large distances to be studied, and so has tremendous potential for extragalactic abundance work. Here, we present an abundance study of the Large and Small Magellanic Clouds (LMC and SMC respectively) using samples of 9-10 RSGs in each. We find average abundances for the two galaxies of [Z]LMC = -0.37 +/- 0.14 and [Z]SMC = -0.53 +/- 0.16 (with respect to a Solar metallicity of Zsun=0.012). These values are consistent with other studies of young stars in these galaxies, and though our result for the SMC may appear high it is consistent with recent studies of hot stars which find 0.5-0.8dex below Solar. Our best-fit temperatures are on the whole consistent with those from fits to the optical-infrared spectral energy distributions, which is remarkable considering the narrow spectral range being studied. Combined with our recent study of RSGs in the Galactic cluster Per OB1, these results indicate that this technique performs well over a range of metallicities, paving the way for forthcoming studies of more distant galaxies beyond the Local Group.

We present near-IR spectroscopy of red supergiant (RSG) stars in NGC 6822, obtained with the new VLT-KMOS instrument. From comparisons with model spectra in the J-band we determine the metallicity of 11 RSGs, finding a mean value of [Z] = -0.52 $\pm$ 0.21 which agrees well with previous abundance studies of young stars and HII regions. We also find an indication for a low-significance abundance gradient within the central 1 kpc. We compare our results with those derived from older stellar populations and investigate the difference using a simple chemical evolution model. By comparing the physical properties determined for RSGs in NGC 6822 with those derived using the same technique in the Galaxy and the Magellanic Clouds, we show that there appears to be no significant temperature variation of RSGs with respect to metallicity, in contrast with recent evolutionary models.

We demonstrate how the metallicities of young super star clusters can be measured using novel spectroscopic techniques in the J-band. The near-infrared flux of super star clusters older than ~6 Myr is dominated by tens to hundreds of red supergiant stars. Our technique is designed to harness the integrated light of that population and produces accurate metallicities for new observations in galaxies above (M83) and below (NGC 6946) solar metallicity. In M83 we find [Z]= +0.28 +/- 0.14 dex using a moderate resolution (R~3500) J-band spectrum and in NGC 6496 we report [Z]= -0.32 +/- 0.20 dex from a low resolution spectrum of R~1800. Recently commissioned low resolution multiplexed spectrographs on the VLT (KMOS) and Keck (MOSFIRE) will allow accurate measurements of super star cluster metallicities across the disks of star-forming galaxies up to distances of 70 Mpc with single night observation campaigns using the method presented in this letter.

We report follow-up observations of 477 program Near-Earth Asteroids (NEAs) using nine telescopes of the EURONEAR network having apertures between 0.3 and 4.2 m. Adding these NEAs to our previous results we now count 739 program NEAs followed-up by the EURONEAR network since 2006. The targets were selected using EURONEAR planning tools focusing on high priority objects. Analyzing the resulting orbital improvements suggests astrometric follow-up is most important days to weeks after discovery, with recovery at a new opposition also valuable. Additionally we observed 40 survey fields spanning three nights covering 11 sq. degrees near opposition, using the Wide Field Camera on the 2.5m Isaac Newton Telescope (INT), resulting in 104 discovered main belt asteroids (MBAs) and another 626 unknown one-night objects. These fields, plus program NEA fields from the INT and from the wide field MOSAIC II camera on the Blanco 4m telescope, generated around 12,000 observations of 2,000 minor planets (mostly MBAs) observed in 34 square degrees. We identify Near Earth Object (NEO) candidates among the unknown (single night) objects using three selection criteria. Testing these criteria on the (known) program NEAs shows the best selection methods are our epsilon-miu model which checks solar elongation and sky motion and the MPC's NEO rating tool. Our new data show that on average 0.5 NEO candidates per square degree should be observable in a 2m-class survey (in agreement with past results), while an average of 2.7 NEO candidates per square degree should be observable in a 4m-class survey (although our Blanco statistics were affected by clouds). At opposition just over 100 MBAs (1.6 unknown to every 1 known) per square degree are detectable to R=22 in a 2m survey based on the INT data, while our two best ecliptic Blanco fields away from opposition lead to 135 MBAs (2 unknown to every 1 known) to R=23.