Search

We aim to investigate the nature of time-variable X-ray sources detected in the {\it XMM-Newton} serendipitous survey. The X-ray light curves of objects in the {\it XMM-Newton} serendipitous survey were searched for variability and coincident serendipitous sources observed by {\it Chandra} were also investigated. Subsequent infrared spectroscopy of the counterparts to the X-ray objects that were identified using UKIDSS was carried out using {\it ISAAC} on the VLT. We found that the object 4XMM~J182531.5--144036 detected in the {\it XMM-Newton} serendipitous survey in April 2008 was also detected by {\it Chandra} as CXOU~J182531.4--144036 in July 2004. Both observations reveal a hard X-ray source displaying a coherent X-ray pulsation at a period of 781~s. The source position is coincident with a $K=14$ mag infrared object whose spectrum exhibits strong HeI and Br$\gamma$ emission lines and an infrared excess above that of early B-type dwarf or giant stars. We conclude that 4XMM~J182531.5--144036 is a Be/X-ray binary pulsar exhibiting persistent X-ray emission and is likely in a long period, low eccentricity orbit, similar to X Per., Comment: 6 pages, 8 figures, accepted for publication in Publications of the Astronomical Society of Australia

Serving as the progenitors of electromagnetic and gravitational wave transients, massive stars have received renewed interest in recent years. However, many aspects of their birth and evolution remain opaque, particularly in the context of binary interactions. The centre of our galaxy hosts a rich cohort of very massive stars, which appear to play a prominent role in the ecology of the region. In this paper we investigate the binary properties of the Arches cluster, which is thought to host a large number of very massive stars. A combination of multi-epoch near-IR spectroscopy and photometry was utilised to identify binaries. 13 from 36 cluster members meet our criteria to be classed as RV variable. Combining the spectroscopic data with archival radio and X-ray observations - to detect colliding wind systems - provides a lower limit to the binary fraction of ~43%; increasing to >50% for the O-type hypergiants and WNLha. Dynamical and evolutionary masses reveal the primaries to be uniformly massive (>50M$_{\odot}$). Where available, orbital analysis reveals a number of short period, highly eccentric binaries, which appear to be pre-interaction systems. Such systems are X-ray luminous, with 80% above an empirical bound of $(L_{\rm x}/L_{\rm bol})\sim10^{-7}$ and their orbital configurations suggest formation and evolution via a single star channel; however, we cannot exclude a binary formation channel for a subset. Qualitative comparison to surveys of lower mass OB-type stars confirms that the trend to an extreme binary fraction (>60%) extends to the most massive stars currently forming in the local Universe., Comment: Accepted in MNRAS; In Memoriam: This work is dedicated to the memory of our dearest friend and colleague Simon Clark

The galactic cluster Westerlund 1 contains a rich population of evolved, massive stars, and a high binary fraction has been inferred from previous multiwavelength observations. We use multi-epoch spectroscopy of a large sample of early-type stars to identify new binaries and binary candidates in the cluster. VLT/FLAMES was used to obtain spectra of ~100 OB stars over a 14-month baseline in 2008 and 2009, supplemented with follow-up observations in 2011 and 2013, and we identify 20 new OB I--III binaries, a WN9h binary, and a WC9d binary, greatly increasing the number of directly confirmed binary systems in Westerlund 1, while 12 O9--9.5 Iab--III stars are identified as candidate binaries. The 173.9 day SB1 W1030 represents the first longer-period system identified in the cluster, while the determination of a 53.95 day period for W44/L makes it the first Wolf-Rayet binary in Westerlund 1 with a confirmed orbital period greater than ten days. Our results suggest the binary fraction in the OB population is at least 40%, and may be significantly higher. (Abridged), Comment: Accepted for publication in Astronomy & Astrophysics, and dedicated to the memory of Dr Simon Clark

Since July 2014, the Gaia mission has been engaged in a high-spatial-resolution, time-resolved, precise, accurate astrometric, and photometric survey of the entire sky. Aims: We present the Gaia Science Alerts project, which has been in operation since 1 June 2016. We describe the system which has been developed to enable the discovery and publication of transient photometric events as seen by Gaia. Methods: We outline the data handling, timings, and performances, and we describe the transient detection algorithms and filtering procedures needed to manage the high false alarm rate. We identify two classes of events: (1) sources which are new to Gaia and (2) Gaia sources which have undergone a significant brightening or fading. Validation of the Gaia transit astrometry and photometry was performed, followed by testing of the source environment to minimise contamination from Solar System objects, bright stars, and fainter near-neighbours. Results: We show that the Gaia Science Alerts project suffers from very low contamination, that is there are very few false-positives. We find that the external completeness for supernovae, $C_E=0.46$, is dominated by the Gaia scanning law and the requirement of detections from both fields-of-view. Where we have two or more scans the internal completeness is $C_I=0.79$ at 3 arcsec or larger from the centres of galaxies, but it drops closer in, especially within 1 arcsec. Conclusions: The per-transit photometry for Gaia transients is precise to 1 per cent at $G=13$, and 3 per cent at $G=19$. The per-transit astrometry is accurate to 55 milliarcseconds when compared to Gaia DR2. The Gaia Science Alerts project is one of the most homogeneous and productive transient surveys in operation, and it is the only survey which covers the whole sky at high spatial resolution (subarcsecond), including the Galactic plane and bulge., Comment: 26 pages, 26 figures, accepted for publication in Astronomy & Astrophysics

R 144 is the visually brightest WR star in the Large Magellanic Cloud (LMC). R 144 was reported to be a binary, making it potentially the most massive binary thus observed. We perform a comprehensive spectral, photometric, orbital, and polarimetric analysis of R 144. R 144 is an eccentric (e=0.51) 74.2-d binary comprising two relatively evolved (age~2 Myr), H-rich WR stars. The hotter primary (WN5/6h, T=50 kK) and the cooler secondary (WN6/7h,T=45kK) have nearly equal masses. The combination of low rotation and H-depletion observed in the system is well reproduced by contemporary evolution models that include boosted mass-loss at the upper-mass end. The systemic velocity of R 144 and its relative isolation suggest that it was ejected as a runaway from the neighbouring R 136 cluster. The optical light-curve shows a clear orbital modulation that can be well explained as a combination of two processes: excess emission stemming from wind-wind collisions and double wind eclipses. Our light-curve model implies an orbital inclination of i=60.4+-1.5deg, resulting in accurately constrained dynamical masses of 74+-4 and 69+-4 Msun. Assuming that both binary components are core H-burning, these masses are difficult to reconcile with the derived luminosities (logL1,2 = 6.44, 6.39 [Lsun]), which correspond to evolutionary masses of the order of 110 and 100Msun, respectively. Taken at face value, our results imply that both stars have high classical Eddington factors of Gamma_e = 0.78+-0.1. If the stars are on the main sequence, their derived radii (~25Rsun) suggest that they are only slightly inflated, even at this high Eddington factor. Alternatively, the stars could be core-He burning, strongly inflated from the regular size of classical Wolf-Rayet stars (~1Rsun), a scenario that could help resolve the observed mass discrepancy., Comment: Accepted to A&A on the April 06 2021. 17 pages + 5 page appendix. Version 28.10.2021: corrected typo in Eq. (2); numerator and denominator were inversed on RHS (only typo; no impact on results)

Many galaxies host pronounced circumnuclear starbursts, fuelled by infalling gas. Such activity is expected to drive the secular evolution of the nucleus and generate super winds, while the intense radiation fields and extreme gas and cosmic ray densities present may act to modify the outcome of star formation with respect to more quiescent galactic regions. The centre of the Milky Way is the only example of this phenomenon where, by virtue of its proximity, individual stars may be resolved. Previous studies have revealed that it hosts a rich population of massive stars; these are located within three clusters, with an additional contingent dispersed throughout the Central Molecular Zone (CMZ). We employed VLT+KMOS to obtain homogeneous, high S/N spectroscopy of the later cohort for classification and quantitative analysis. Including previously identified examples, we found a total of 83 isolated massive stars within the Galactic Centre, which are biased towards objects supporting powerful stellar winds and/or extensive circumstellar envelopes. No further stellar clusters, or their tidally stripped remnants, were identified, although an apparent stellar overdensity was found to be coincident with the Sgr B1 star forming region. The cohort of isolated massive stars within the CMZ is comparable in size to that of the known clusters but, due to observational biases, is likely highly incomplete at this time. Combining both populations yields over 320 spectroscopically classified stars that are expected to undergo core collapse within the next 20Myr. Given that this is presumably an underestimate of the true number, the population of massive stars associated with the CMZ appears unprecedented amongst star formation complexes within the Milky Way, and one might anticipate that they play a substantial role in the energetics and evolution of the nuclear region., Comment: 28 pages, 14 figures, accepted for publication in Astronomy and Astrophysics

The role of massive stars is central to an understanding of galactic ecology. It is important to establish the details of how massive stars provide radiative, chemical, and mechanical feedback in galaxies. Central to these issues is an understanding of the evolution of massive stars, and the critical role of mass loss via strongly structured winds and stellar binarity. Ultimately, massive stellar clusters shape the structure and energetics of galaxies. We aim to conduct high-resolution, deep field mapping at 21cm of the core of the massive Cygnus OB2 association and to characterise the properties of the massive stars and colliding winds at this waveband. We used seven stations of the e-MERLIN radio facility, with its upgraded bandwidth and enhanced sensitivity to conduct a 21cm census of Cygnus OB2. Based on 42 hours of observations, seven overlapping pointings were employed over multiple epochs during 2014 resulting in 1 sigma sensitivities down to ~21microJy and a resolution of ~180mas. A total of 61 sources are detected at 21cm over a ~0.48deg x 0.48deg region centred on the heart of the Cyg OB2 association. Of these 61 sources, 33 are detected for the first time. We detect a number of previously identified sources including four massive stellar binary systems, two YSOs, and several known X-ray and radio sources. We also detect the LBV candidate (possible binary system) and blue hypergiant (BHG) star of Cyg OB2 #12. The 21cm observations secured in the COBRaS Legacy project provide data to constrain conditions in the outer wind regions of massive stars; determine the non-thermal properties of massive interacting binaries; examine evidence for transient sources, including those associated with young stellar objects; and provide unidentified sources that merit follow-up observations. The 21cm data are of lasting value and will serve in combination with other key surveys of Cyg OB2., Comment: 41 pages, 12 figures, accepted in A&A

Accurate stellar parameters of individual objects in binary systems are essential to constrain the effects of binarity on stellar evolution. These parameters serve as a prerequisite to probing existing and future theoretical evolutionary models. We aim to derive the atmospheric parameters of the 31 SB2s in the TMBM sample. This sample, composed of detached, semi-detached and contact systems with at least one of the components classified as an O star, is an excellent test-bed to study how binarity can impact our knowledge of the evolution of massive stars. 32 epochs of FLAMES/GIRAFFE spectra are analysed using spectral disentangling to construct the individual spectra of 62 components. We apply the CMFGEN atmosphere code to determine their stellar parameters and their He, C and N surface abundances. From these properties, we show that the effects of tides on chemical mixing are limited. Components on longer-period orbits show higher nitrogen enrichment at their surface than those on shorter-period orbits, in contrast to expectations of rotational or tidal mixing, implying that other mechanisms play a role in this process. Components filling their Roche lobe are mass donors. They exhibit higher nitrogen content at their surface and rotate more slowly than their companions. By accreting new material, their companions spin faster and are rejuvenated. Their locations in the N-vsini diagram tend to show that binary products are good candidates to populate the two groups of stars (slowly rotating, nitrogen-enriched and rapidly rotating non-enriched) that cannot be reproduced through single-star population synthesis. This sample is the largest sample of binaries to be studied in such a homogeneous way. The study of these objects gives us strong observational constraints to test theoretical binary evolutionary tracks., Comment: accepted to A&A Core of the paper 15 pages, Appendix: 41 pages

A high fraction of massive stars are found to be binaries but only a few of them are reported as photometrically variable. By studying the populations of SB2 in the 30 Doradus region, we found a subset of them that have photometry from the OGLE project and that display variations in their light curves related to orbital motions. The goal of this study is to determine the dynamical masses and radii of the 26 binary components to investigate the mass-discrepancy problem and to provide an empirical mass-luminosity relation for the LMC. We use the PHOEBE programme to perform a systematic analysis of the OGLE V and I light curves obtained for 13 binary systems in 30 Dor. We adopt Teff, and orbital parameters derived previously to obtain the inclinations of the systems and the parameters of the individual components. Three systems display eclipses in their light curves, while the others only display ellipsoidal variations. We classify two systems as over-contact, five as semi-detached, and four as detached. The two remaining systems have uncertain configurations due to large uncertainties on their inclinations. The fact that systems display ellipsoidal variations has a significant impact on the inclination errors. From the dynamical masses, luminosities, and radii, we provide LMC-based empirical mass-luminosity and mass-radius relations, and we compare them to other relations given for the Galaxy, the LMC, and the SMC. These relations differ for different mass ranges, but do not seem to depend on the metallicity regimes. We also compare the dynamical, spectroscopic, and evolutionary masses of the stars in our sample. While the dynamical and spectroscopic masses agree with each other, the evolutionary masses are systematically higher, at least for stars in semi-detached systems. This suggests that the mass discrepancy can be partly explained by past or ongoing interactions between the stars., Comment: Accepted to A&A The core of the paper 9 pages, 5 pages of Appendix

Massive stars and their stellar winds are important for a number of feedback processes. The mass lost in the stellar wind can help determine the end-point of the star as a NS or a BH. However, the impact of mass-loss on the post-Main Sequence evolutionary stage of massive stars is not well understood. Westerlund 1 is an ideal astrophysical laboratory in which to study massive stars and their winds in great detail over a large range of different evolutionary phases. Aims: We aim to study the radio emission from Westerlund 1, in order to measure radio fluxes from the population of massive stars, and determine mass-loss rates and spectral indices where possible. Methods: Observations were carried out in 2015 and 2016 with the Australia telescope compact array (ATCA) at 5.5 and 9 GHz using multiple configurations, with maximum baselines ranging from 750m to 6km. Results: 30 stars were detected in the radio from the fully concatenated dataset, 10 of which were WRs (predominantly late type WN stars), 5 YHGs, 4 RSGs, 1 LBV star, the sgB[e] star W9, and several O and B supergiants. New source detections in the radio were found for 5 WR stars, and 5 OB supergiants. These detections have led to evidence for 3 new OB supergiant binary candidates, inferred from derived spectral index limits. Conclusions: Spectral indices and index limits were determined for massive stars in Westerlund 1. For cluster members found to have partially optically thick emission, mass-loss rates were calculated. Under the approximation of a thermally emitting stellar wind and a steady mass-loss rate, clumping ratios were then estimated for 8 WRs. Diffuse radio emission was detected throughout the cluster. Detections of knots of radio emission with no known stellar counterparts indicate the highly clumped structure of this intra-cluster medium, likely shaped by a dense cluster wind., Comment: 22 pages, 10 figures, accepted for publication in A&A

The formation, properties, and evolution of massive stars remain subject to considerable uncertainty; impacting on fields as diverse as galactic feedback and the nature of the progenitors of both electromagnetic and gravitational wave transients. The clusters many such stars reside within provide a unique laboratory for addressing these issues, and in this work we provide a comprehensive stellar census of Westerlund 1 to underpin such efforts. 69 new members were identified via I-band spectroscopy, yielding a total cluster population of 166 stars with initial masses of ~25Msun to ~50Msun, with more massive stars already lost to supernova. The stellar population follows a smooth and continuous morphological sequence from late-O giant through to OB supergiant. Subsequently, the progression bifurcates, with one branch yielding mid- to late-B hypergiants and cool super-/hypergiants, and the other massive blue stragglers, prior to a diverse population of H-depleted Wolf-Rayets. A substantial population of O-type stars with anomalously broad Paschen series lines are seen, a property which we attribute to binarity. Binary interaction is clearly required to yield the uniquely rich cohort of hypergiants, which includes both mass-stripped primaries and rejuvenated secondaries/stellar mergers. As a consequence future observations of Wd1 and similar stellar aggregates hold out the prospect of characterising both single- and binary- evolutionary channels for massive stars and determining their relative contributions. This in turn will permit the physical properties of such objects at the point of core-collapse to be predicted; of direct relevance for understanding the formation of relativistic remnants such as the magnetars associated with Wd1 and other young massive clusters (Abridged)., Comment: 21 pages, 9 figures, accepted for publication in Astronomy and Astrophysics

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], Comment: accepted in A&A

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., Comment: 13 pages, 8 figures, accepted for publication in Astronomy and Astrophysics journal

Despite the first detection of X-rays from massive stars occurring four decades ago, the physical dependence of the emission mechanism(s) on the underlying stellar and binary properties of the emitters remains uncertain. The young massive cluster Westerlund 1 provides an ideal testbed for understanding this phenomenon, with over 50 cluster members detected in historical X-ray observations. In the decade since these data were obtained, significant new multi-epoch observations of the cluster have been made, allowing a fundamental reappraisal of the nature of both X-ray bright and dark stars. A total of 45 X-ray sources within Wd1 now have precise spectral classifications. These observations identify 16 candidate and confirmed massive binaries; by comparison 22 stars have X-ray properties that imply a contribution from a wind collision zone. X-ray emission appears confined to O9-B0.5 supergiants, Wolf-Rayets and a small group of highly luminous interacting/post-interaction binaries. Despite their presence in large numbers, no emission is seen from earlier, less evolved O stars or later, cooler B super-/hypergiants. We suppose that the lack of X-ray emission from O giants is due to their comparatively low bolometric luminosities if, as expected, they follow the canonical L_X/L_bol relation for hot stars. The transition away from X-ray emission for OB supergiants occurs at the location of the bistability jump; we speculate that below this limit, stellar wind velocities are insufficient for internal, X-ray emitting shocks to form. Our results are consistent with recent findings that massive binaries are not uniformly brighter than single stars of comparable luminosity and spectral type, although it is noteworthy that the brightest and hardest stellar X-ray sources within Wd1 are all either confirmed or candidate massive, interacting/post-interaction binaries (Abridged)., Comment: 10 pages, 2 figures, accepted for publication in Astronomy and Astrophysics

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., Comment: Accepted by A&A

The Galactic supergiant B[e] star CI Camelopardalis (CI Cam) was the first sgB[e] star detected during an X-ray outburst. The star brightened to $\sim$2 Crab in the X-ray regime within hours before decaying to a quiescent level in less than 2 weeks, clearly indicative of binarity. Since the outburst of CI Cam, a number of sgB[e] stars have been identified as X-ray overluminous for a single star (i.e. $L_X > 10^{-7}~L_{bol}$). This small population has recently expanded to include two Ultra Luminous X-ray sources (ULX), Holmberg II X-1 and NGC300 ULX-1/supernova imposter SN2010da. We revisit CI Cam to investigate its behaviour over several timescales and shed further light on the nature of the compact object in the system, its X-ray outburst in 1998 and the binary system parameters. We analyse archival XMM-Newton EPIC-pn spectra and light curves along with new data from Swift and NuSTAR. We also present high-resolution MERCATOR/HERMES optical spectra, including a spectrum taken 1.02 days after our NuSTAR observation. Despite being in quiescence, CI Cam is highly X-ray variable on timescales of days, both in terms of total integrated flux and spectral shape. We interpret these variations by invoking the presence of an accreting compact companion immersed in a dense, highly structured, aspherical circumstellar envelope. The differences in the accretion flux and circumstellar extinction represent either changes in this environment, triggered by variable mass loss from the star, or the local conditions to the accretor due to its orbital motion. We find no evidence for pulsations in the X-ray light curve. CI Cam has many similarities with SN2010da across mid-IR, optical and X-ray wavelengths suggesting that, subject to distance determination for CI Cam, if CI Cam was located in an external galaxy its 1998 outburst would have led to a classification as a supernova imposter., Comment: 13 pages, 8 Figures and 3 Tables. Accepted for Publication in A&A

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., Comment: 10 pages, 6 figures, accepted for publication in Astronomy and Astrophysics

Observational studies suggest that many OB stars are found within binary systems which may be expected to interact during their lifetimes. Significant mass transfer or merger of both components will modify evolutionary pathways, facilitating the production of exceptionally massive stars which will present as blue stragglers. Identification and characterisation of such objects is crucial if the outcomes of binary evolutionary channels are to be quantified. The massive cluster Westerlund 1 hosts a rich population of X-ray bright stars where the emission is thought to derive from the wind collision zones of massive binaries. Selected on this basis, we present the results of a multiwavelength analysis of the X-ray luminous O stars Wd1-27 and -30a. We find both to be early/mid-O hypergiants with luminosities, temperatures and masses significantly in excess of other early stars within Wd1, hence qualifying as massive blue stragglers. The nature of Wd1-27 remains uncertain but the detection of radial velocity changes and the X-ray properties of Wd1-30a suggest that it is a binary. Analysis of Gaia proper motion and parallactic data indicates that both are cluster members; we also provide a membership list for Wd1 based on this analysis. The presence of hypergiants of spectral types O to M within Wd1 cannot be understood via single-star evolution. We suppose that the early-B and mid-O hypergiants formed via binary-induced mass-stripping of the primary and mass-transfer to the secondary, respectively. This implies that for a subset of objects massive star-formation may be regarded as a two-stage process, with binary-driven mass-transfer or merger yielding stars with masses significantly in excess of their initial birth mass (Abridged)., Comment: 19 pages, 9 figures, accepted for publication in Astronomy and Astrophysics

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)., Comment: 27 pages, 12 figures, accepted for publication in Astronomy and Astrophysics

As part of a wider investigation of evolved massive stars in Galactic open clusters, we have spectroscopically identified three candidate classical Cepheids in the little-studied clusters Berkeley 51, Berkeley 55 and NGC 6603. Using new multi-epoch photometry, we confirm that Be 51 #162 and Be 55 #107 are bona fide Cepheids, with pulsation periods of 9.83+/-0.01 d and 5.850+/-0.005 d respectively, while NGC 6603 star W2249 does not show significant photometric variability. Using the period-luminosity relationship for Cepheid variables, we determine a distance to Be 51 of 5.3(+1.0,-0.8) kpc and an age of 44(+9,-8) Myr, placing it in a sparsely-attested region of the Perseus arm. For Be 55, we find a distance of 2.2+/-0.3 kpc and age of 63(+12,-11) Myr, locating the cluster in the Local arm. Taken together with our recent discovery of a long-period Cepheid in the starburst cluster VdBH222, these represent an important increase in the number of young, massive Cepheids known in Galactic open clusters. We also consider new Gaia (data release 2) parallaxes and proper motions for members of Be 51 and Be 55; the uncertainties on the parallaxes do not allow us to refine our distance estimates to these clusters, but the well-constrained proper motion measurements furnish further confirmation of cluster membership. However, future final Gaia parallaxes for such objects should provide valuable independent distance measurements, improving the calibration of the period-luminosity relationship, with implications for the distance ladder out to cosmological scales., Comment: 8 pages, 10 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society. Author details corrected

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., Comment: 15 pages, 18 figures, accepted for publication in Astronomy & Astrophysics. Paper II in a series; for Paper I see arXiv:1803.09567

Massive stars play an important role in both cluster and galactic evolution and the rate at which they lose mass is a key driver of both their own evolution and their interaction with the environment up to and including their SNe explosions. Young massive clusters provide an ideal opportunity to study a co-eval population of massive stars. We performed 3mm continuum observations with the Atacama Large Millimetre/submillimetre Array of the Galactic cluster Westerlund 1, to study the constituent massive stars and determine mass-loss rates for the diverse post-main sequence population. We detected emission from 50 stars in Westerlund 1, comprising all 21 Wolf-Rayets within the field of view, eight cool and 21 OB super-/hypergiants. Emission nebulae were associated with a number of the cool hypergiants while, unexpectedly, a number of hot stars also appear spatially resolved. We measured the mass-loss rates for a unique population of massive post-main sequence stars at every stage of evolution, confirming a significant increase as stars transition from OB supergiant to WR states. The range of spectral types exhibited provides a critical test of radiatively driven wind theory and the reality of the bi-stability jump. The extreme mass-loss rate inferred for the interacting binary Wd1-9 in comparison to other cluster members confirmed the key role binarity plays in massive stellar evolution. The presence of compact nebulae around a number of OB and WR stars is unexpected; by analogy to the cool super-/hypergiants we attribute this to confinement and sculpting of the stellar wind via interaction with the intra-cluster medium/wind. Given the morphology of core collapse SNe depend on the nature of the pre-explosion circumstellar environment, if this hypothesis is correct then the properties of the explosion depend not just on the progenitor, but also the environment in which it is located., Comment: 35 pages, 20 figures, accepted for publication in A&A

Located within the central region of the Galaxy, the Arches cluster appears to be one of the youngest, densest and most massive stellar aggregates within the Milky Way. As such it has the potential to be a uniquely instructive laboratory for the study of star formation in extreme environments and the physics of very massive stars. In order to determine the fundamental physical properties of both cluster and constituent stars, we provide and analyse new HST+VLT near-IR datasets. Stacking multiple epochs of spectroscopy results in the deepest view of the cluster ever obtained, allowing us to to identify candidate giant and main sequence stars for the first time. All cluster members are found to be WNLh or O stars, with the smooth and continuous progression in spectral morphologies from O super-/hypergiants through to the WNLh cohort implying a direct evolutionary connection. Importantly no H-free Wolf-Rayets are found, while no products of binary interaction/mass-transfer may be unambiguously identified, despite the presence of massive binaries within the Arches. We infer a main sequence turn-off around O4-5V, corresponding to ~30-38Msun, while the eclipsing binary F2 implies current masses of ~80Msun and ~60Msun for the WNLh and O hypergiant cohorts, respectively. A cluster age of ~2-3Myr is suggested by the location of the main-sequence turn-off. While the absence of H-free Wolf-Rayets argues against the prior occurrence of SNe, such an age does accommodate such events for exceptionally massive stars. Future progress requires quantitative analysis of cluster members combined with additional spectroscopic observations to better constrain the binary population; nevertheless it is already abundantly clear that the Arches offers an unprecedented insight into the formation, evolution and death of the most massive stars Nature allows to form in the local universe (Abridged)., Comment: 34 pages, 12 figures, accepted for publication in Astronomy and Astrophysics

We report new 5.5 GHz radio observations of the massive star cluster Westerlund 1, taken by the Australia Telescope Compact Array, detecting nine of the ten yellow hypergiants (YHGs) and red supergiants (RSGs) within the cluster. Eight of nine sources are spatially resolved. The nebulae associated with the YHGs Wd1-4a, -12a and -265 demonstrate a cometary morphology - the first time this phenomenon has been observed for such stars. This structure is also echoed in the ejecta of the RSGs Wd1-20 and -26; in each case the cometary tails are directed away from the cluster core. The nebular emission around the RSG Wd1-237 is less collimated than these systems but once again appears more prominent in the hemisphere facing the cluster. Considered as a whole, the nebular morphologies provide compelling evidence for sculpting via a physical agent associated with Westerlund 1, such as a cluster wind., Comment: 5 pages, 2 figures, Acccepted for publication in MNRAS

We present a detailed optical and X-ray study of the 2013 outburst of the Small Magellanic Cloud Be/X-ray binary SXP7.92, as well as an overview of the last 18 years of observations from OGLE, RXTE, Chandra and XMM-Newton. We revise the position of this source to RA(J2000)=00:57:58.4, Dec(J2000)=-72:22:29.5 with a $1\sigma$ uncertainty of 1.5 arcsec, correcting the previously reported position by Coe et al. (2009) by more than 20 arcminutes. We identify and spectrally classify the correct counterpart as a B1Ve star. The optical spectrum is distinguished by an uncharacteristically deep narrow Balmer series, with the H$\alpha$ line in particular having a distinctive shell profile, i.e. a deep absorption core embedded in an emission line. We interpret this as evidence that we are viewing the system edge on and are seeing self obscuration of the circumstellar disc. We derive an optical period for the system of 40.0$\pm0.3$ days, which we interpret as the orbital period, and present several mechanisms to describe the X-ray/Optical behaviour in the recent outburst, in particular the "flares" and "dips" seen in the optical light curve, including a transient accretion disc and an elongated precessing disc., Comment: 15 pages, 7 figures, 3 tables, accepted for publication in MNRAS

Mass-loss in massive stars plays a critical role in their evolution, although the precise mechanism(s) responsible - radiatively driven winds, impulsive ejection and/or binary interaction -remain uncertain. In this paper we present ALMA line and continuum observations of the supergiant B[e] star Wd1-9, a massive post-Main Sequence object located within the starburst cluster Westerlund 1. We find it to be one of the brightest stellar point sources in the sky at millimetre wavelengths, with (serendipitously identified) emission in the H41alpha radio recombination line. We attribute these properties to a low velocity (~100 km/s) ionised wind, with an extreme mass-loss rate 6.4x10^-5(d/5kpc)^1.5 Msol/yr. External to this is an extended aspherical ejection nebula indicative of a prior phase of significant mass-loss. Taken together, the millimetre properties of Wd1-9 show a remarkable similarity to those of the highly luminous stellar source MWC349A.We conclude that these objects are interacting binaries evolving away from the main sequence and undergoing rapid case-A mass transfer. As such they - and by extension the wider class of supergiant B[e] stars - may provide a unique window into the physics of a process that shapes the life-cycle of ~70% of massive stars found in binary systems., Comment: 6 pages, 3 figures, MNRAS accepted

Galactic starburst clusters play a twin role in astrophysics, serving as laboratories for the study of stellar physics and also delineating the structure and recent star formation history of the Milky Way. In order to exploit these opportunities we have undertaken a multi-epoch spectroscopic survey of the red supergiant dominated young massive clusters thought to be present at both near and far ends of the Galactic Bar. Significant spectroscopic variability suggestive of radial pulsations was found for the yellow supergiant VdBH 222 #505. Follow-up photometric investigations revealed modulation with a period of ~23.325d; both timescale and pulsational profile are consistent with a Cepheid classification. As a consequence #505 may be recognised as one of the longest period Galactic cluster Cepheids identified to date and hence of considerable use in constraining the bright end of the period/luminosity relation at solar metallicities. In conjunction with extant photometry we infer a distance of ~6kpc for VdBH222 and an age of ~20Myr. This results in a moderate reduction in both integrated cluster mass (~2x10^4Msun) and the initial stellar masses of the evolved cluster members (~10Msun). As such, VdBH222 becomes an excellent test-bed for studying the properties of some of the lowest mass stars observed to undergo type-II supernovae. Moreover, the distance is in tension with a location of VdBH 222 at the far end of the Galactic Bar. Instead a birthsite in the near 3kpc arm is suggested; providing compelling evidence of extensive recent star formation in a region of the inner Milky Way which has hitherto been thought to be devoid of such activity., Comment: 5 pages, 3 figures, accepted for publication in Astronomy and Astrophysics

We investigate the characteristics of two newly discovered short-period, double-lined, massive binary systems, VFTS 450 (O9.7$\;$II--Ib$\,$+$\,$O7::) and VFTS 652 (B1$\;$Ib$\,+\,$O9:$\;$III:). We perform model-atmosphere analyses to characterise the photospheric properties of both members of each binary (denoting the `primary' as the spectroscopically more conspicuous component). Radial velocities and optical photometry are used to estimate the binary-system parameters. We estimate $T_{\rm eff}=27$ kK, $\log{(g)}=2.9$ (cgs) for the VFTS 450 primary spectrum (34kK, 3.6: for the secondary spectrum); and $T_{\rm eff} = 22$kK, $\log{(g)}=2.8$ for the VFTS 652 primary spectrum (35kK, 3.7: for the secondary spectrum). Both primaries show surface nitrogen enrichments (of more than 1 dex for VFTS 652), and probable moderate oxygen depletions relative to reference LMC abundances. We determine orbital periods of 6.89d and 8.59d for VFTS 450 and VFTS 652, respectively, and argue that the primaries must be close to filling their Roche lobes. Supposing this to be the case, we estimate component masses in the range $\sim$20--50M$_\odot$. The secondary spectra are associated with the more massive components, suggesting that both systems are high-mass analogues of classical Algol systems, undergoing case-A mass transfer. Difficulties in reconciling the spectroscopic analyses with the light-curves and with evolutionary considerations suggest that the secondary spectra are contaminated by (or arise in) accretion disks.

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

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

We have analyzed a ~130 ks XMM-Newton observation of the dynamically confirmed black hole + Wolf-Rayet (BH+WR) X-ray binary (XB) IC10 X-1, covering ~1 orbital cycle. This system experiences periodic intensity dips every ~35 hours. We find that energy-independent evolution is rejected at a >5 sigma level. The spectral and timing evolution of IC10 X-1 are best explained by a compact disk blackbody and an extended Comptonized component, where the thermal component is completely absorbed and the Comptonized component is partially covered during the dip. We consider three possibilities for the absorber: cold material in the outer accretion disk, as is well documented for Galactic neutron star (NS) XBs at high inclination; a stream of stellar wind that is enhanced by traveling through the L1 point; and a spherical wind. We estimated the corona radius (r_ADC) for IC10 X-1 from the dip ingress to be ~1 E+6 km, assuming absorption from the outer disk, and found it to be consistent with the relation between r_m ADC and 1--30 keV luminosity observed in Galactic NS XBs that spans 2 orders of magnitude. For the other two scenarios, the corona would be larger. Prior BH mass (M_BH) estimates range over 23--38 M_Sun, depending on the inclination and WR mass. For disk absorption, the inclination, i, is likely to be ~60--80 degrees, with M_BH ~24--41 M_Sun. Alternatively, the L1-enhanced wind requires i ~80 degrees, suggesting ~24--33 M_Sun. For a spherical absorber, i ~ 40 degrees, and M_BH ~50--65 M_Sun., Comment: Accepted for publication in ApJ; 10 pages, 4 figures

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

The first soft gamma-ray repeater was discovered over three decades ago, and subsequently identified as a magnetar. However there is currently no consenus on the formation channel(s) of these objects. The presence of a magnetar in the starburst cluster Westerlund 1 implies a >40Msun progenitor, favouring formation in a binary that was disrupted at supernova. To test this hypothesis we searched for the putative pre-SN companion via a radial velocity survey to identify high-velocity runaways, with subsequent atmospheric analysis of the resultant candidate, Wd1-5. Wd1-5 appears overluminous for its spectroscopic mass and we find evidence of He- and N-enrichment, O-depletion, and critically C-enrichment, a combination of properties that is difficult to explain for a single star. We infer a pre-SN history for Wd1-5 which supposes an initial close binary comprising two stars of comparable (~41Msun+35Msun) masses. Efficient mass transfer from the initially more massive component leads to the mass-gainer evolving more rapidly, initiating luminous blue variable/common envelope evolution. Reverse, wind-driven mass transfer during its subsequent WC Wolf-Rayet phase leads to the carbon pollution of Wd1-5, before a type Ibc supernova disrupts the binary system. Under the assumption of a physical association between Wd1-5 and J1647-45, the secondary is identified as the magnetar progenitor; its common envelope evolutionary phase prevents spin-down of its core prior to SN and the seed magnetic field for the magnetar forms either in this phase or during the earlier episode of mass transfer in which it was spun-up. Our results suggest that binarity is a key ingredient in the formation of at least a subset of magnetars by preventing spin-down via core-coupling and potentially generating a seed magnetic field., Comment: 17 pages, 10 figures, accepted for publication in Astronomy and Astrophysics

An overview of a spectroscopic survey for massive stars in the direction of the Galactic giant molecular complex G23.3-0.3 is presented (Messineo et al. 2010, and 2014 A&A submitted). This region is interesting because it is rich in HII regions and supernova remnants (SNRs). A number of 38 early-type stars, a new luminous blue variable, and a dozen of red supergiants were detected. We identified the likely progenitors of the SNRs W41, G22.7-00.2, and G22.7583-0.4917., Comment: Massive Young Star Clusters Near and Far. From the Milky Way to Re-ionization - Fifth in the series of Guillermo Haro Conferences, Puebla, Mexico, December 2-6, 2013

We report the discovery of the B[e] star VFTS 822 in the 30 Doradus star-forming region of the Large Magellanic Cloud, classified by optical spectroscopy from the VLT-FLAMES Tarantula Survey and complementary infrared photometry. VFTS 822 is a relatively low-luminosity (log $L$ = 4.04 $\pm$ 0.25 $L_{\odot}$) B8[e] star. In this Letter, we evaluate the evolutionary status of VFTS 822 and discuss its candidacy as a Herbig B[e] star. If the object is indeed in the pre-main sequence phase, it would present an exciting opportunity to measure mass accretion rates at low metallicity spectroscopically, to understand the effect of metallicity on accretion rates., Comment: 5 pages, 4 Figures

The location of the supergiant (sg)B[e] star Wd1-9 within the coeval starburst cluster Westerlund 1 allows for its placement in a detailed post-Main Sequence evolutionary scheme and hence we have utilised a comprehensive multiwavelength dataset to determine its physical properties and relation to other sgB[e] stars and the evolved stellar population of Wd1. Wd1-9 is found to exhibit the rich optical emission line spectrum that is characteristic of sgB[e] stars, while dust mass, composition and disc geometry determined from mid-IR data resemble those of other well known examples such as R126. Extreme historical and ongoing mass loss is inferred from radio observations, while the X-ray properties of Wd1-9 imply the presence of high temperature plasma within the system and are directly comparable to a number of confirmed short-period colliding wind binaries within Wd1. The most complete explanation for the observational properties of Wd1-9 is that it is a massive interacting binary currently undergoing, or recently exited from, rapid Roche-lobe overflow, supporting the hypothesis that binarity mediates the formation of (a subset of) sgB[e] stars. The mass loss rate of Wd1-9 is consistent with such an assertion, while viable progenitor and descendent systems are present within Wd1 and other examples of sgB[e] binaries have previously been found. Moreover, the rarity of sgB[e] stars - only two examples are identified from a census of ~68 young massive Galactic clusters and associations containing ~600 post-Main Sequence stars - is explicable given the rapidity (~10^4yr) expected for this phase of massive binary evolution., Comment: 16 pages, 7 figures, accepted for publication in Astronomy and Astrophysics

Yellow hypergiants represent a short-lived evolutionary episode experienced by massive stars as they transit to and from a red supergiant phase. As such, their properties provide a critical test of stellar evolutionary theory. The galactic yellow hypergiant IRC +10420 is a cornerstone system for understanding this phase, since it is the strongest post-RSG candidate known, has demonstrated real-time evolution across the Hertzsprung-Russell diagram and been subject to extensive mass loss. Here we report on the discovery of a twin of IRC +10420 - IRAS 18357-0604. Optical and near-IR spectroscopy are used to investigate the physical properties of IRAS 18357-0604, while IR photometry probes the nature of its circumstellar environment. Pronounced spectral similarities between IRAS 18357-0604 and IRC +10420 suggest comparable temperatures and wind geometries, while photometric data reveals a similarly dusty circumstellar environment, although historical mass loss appears to have been heavier in IRC +10420. The systemic velocity implies a distance compatible with the red supergiant dominated complex at the base of the Scutum Crux arm; the resultant luminosity determination is consistent with a physical association but suggests a lower initial mass than inferred for IRC +10420 (<20Msun versus ~40Msun). Evolutionary predictions for the physical properties of supernova progenitors derived from ~18-20Msun stars - or ~12-15Msun stars that have experienced enhanced mass loss as red supergiants - compare favourably with those of IRAS 18357-0604, which in turn appears to be similar to the the progenitor of SN2011dh; it may therefore provide an important insight into the nature of the apparently H-depleted yellow hypergiant progenitors of some type IIb SNe., Comment: 7 pages, 5 figures, accepted for publication in Astronomy and Astrophysics

The mechanism by which supergiant (sg)B[e] stars support cool, dense dusty discs/tori and their physical relationship with other evolved, massive stars such as luminous blue variables is uncertain. In order to investigate both issues we have analysed the long term behaviour of the canonical sgB[e] star LHA 115-S 18. We employed the OGLE II-IV lightcurve to search for (a-)periodic variability and supplemented these data with new and historic spectroscopy. In contrast to historical expectations for sgB[e] stars, S18 is both photometrically and spectroscopically highly variable. The lightcurve is characterised by rapid aperiodic `flaring' throughout the 16 years of observations. Changes in the high excitation emission line component of the spectrum imply evolution in the stellar temperature - as expected for luminous blue variables - although somewhat surprisingly, spectroscopic and photometric variability appears not to be correlated. Characterised by emission in low excitation metallic species, the cool circumstellar torus appears largely unaffected by this behaviour. Finally, in conjunction with intense, highly variable HeII emission, X-ray emission implies the presence of an unseen binary companion. S18 provides observational support for the putative physical association of (a subset of) sgB[e] stars and luminous blue variables. Given the nature of the circumstellar environment of S18 and that luminous blue variables have been suggested as SN progenitors, it is tempting to draw a parallel to the progenitors of SN1987A and the unusual transient SN2009ip. Moreover the likely binary nature of S18 strengthens the possibility that the dusty discs/tori that characterise sgB[e] stars are the result of binary-driven mass-loss; consequently such stars may provide a window on the short lived phase of mass-transfer in massive compact binaries., Comment: 15 pages, 10 figures; accepted for publication in Astronomy and Astrophysics

We present a detailed X-ray study of the 2003 XMM-Newton observation of the High Mass X-ray Binary XTE J0421+560/CI Cam. The continuum of the X-ray spectrum is well described by a flat power law (Gamma=1.0+/-0.2) with a large intrinsic absorbing column (N_H=(4.4+/-0.5)\times10^{23} cm^{-2}). We have decomposed the broad iron line into 3 separate components: Fe\textsc{i}-K$\alpha$, Fe\textsc{i}-K$\beta$ and Fe\textsc{xxiv-xxv}K$\alpha$. It is unclear how both neutral and almost fully ionised iron can exist simultaneously, however we suggest this could be evidence that the compact object is embedded in the circumstellar material. This doesn't appear to be consistent with the X-ray flux and spectrum of the source, which has remained essentially unchanged since the initial outburst. The iron abundance implied by the ratio of the neutral Fe-K$\alpha$ and Fe-K$\beta$ is compatible with solar. We search for lags in the neutral Fe-K$\alpha$ with respect to the continuum and find marginal evidence for a lag at ~10 ks. We interpret this as the light crossing time of the torus which would suggest that the neutral iron is located at a radius of 10 AU. This result depends on several assumptions including the distance to the system, the inclination, the mass of system and the orbital period none of which are known with any great certainty. Better constraints on these system parameters and further observations of this system are required to confirm this result. We discuss the nature of this system in light of our results and place it in context with other binary B[e] stars., Comment: 9 pages, 8 figures, 3 tables accepted for publication in MNRAS

Aims. We analyze the multiplicity properties of the massive O-type star population. With 360 O-type stars, this is the largest homogeneous sample of massive stars analyzed to date. Methods. We use multi-epoch spectroscopy and variability analysis to identify spectroscopic binaries. We also use a Monte-Carlo method to correct for observational biases. Results. We observe a spectroscopic binary fraction of 0.35\pm0.03, which corresponds to the fraction of objects displaying statistically significant radial velocity variations with an amplitude of at least 20km/s. We compute the intrinsic binary fraction to be 0.51\pm0.04. We adopt power-laws to describe the intrinsic period and mass-ratio distributions: f_P ~ (log P)^\pi\ (with 0.15 < log P < 3.5) and f_q ~ q^\kappa\ with 0.1 < q=M_2/M_1 < 1.0. The power-law indexes that best reproduce the observed quantities are \pi = -0.45 +/- 0.30 and \kappa = -1.0\pm0.4. The obtained period distribution thus favours shorter period systems compared to an Oepik law. The mass ratio distribution is slightly skewed towards low mass ratio systems but remains incompatible with a random sampling of a classical mass function. The binary fraction seems mostly uniform across the field of view and independent of the spectral types and luminosity classes. The binary fraction in the outer region of the field of view (r > 7.8', i.e. approx117 pc) and among the O9.7 I/II objects are however significantly lower than expected from statistical fluctuations. Conclusions. Using simple evolutionary considerations, we estimate that over 50% of the current O star population in 30 Dor will exchange mass with its companion within a binary system. This shows that binary interaction is greatly affecting the evolution and fate of massive stars, and must be taken into account to correctly interpret unresolved populations of massive stars., Comment: Accepted for publication in A&A. 23 pages, 20 figures

We present a Herschel far-infrared study towards the rich massive star- forming complex G305, utilising PACS 70, 160 {\mu}m and SPIRE 250, 350, and 500 {\mu}m observations from the Hi-GAL survey of the Galactic plane. The focus of this study is to identify the embedded massive star-forming population within G305, by combining far-infrared data with radio continuum, H2O maser, methanol maser, MIPS, and Red MSX Source survey data available from previous studies. By applying a frequentist technique we are able to identify a sample of the most likely associations within our multi-wavelength dataset, that can then be identified from the derived properties obtained from fitted spectral energy distributions (SEDs). By SED modelling using both a simple modified blackbody and fitting to a comprehensive grid of model SEDs, some 16 candidate associations are identified as embedded massive star-forming regions. We derive a two-selection colour criterion from this sample of log(F70/F500)\geq 1 and log(F160/F350)\geq 1.6 to identify an additional 31 embedded massive star candidates with no associated star-formation tracers. Using this result we can build a picture of the present day star-formation of the complex, and by extrapolating an initial mass function, suggest a current population of \approx 2 \times 10^4 young stellar objects (YSOs) present, corresponding to a star formation rate (SFR) of 0.01-0.02 M\odot yr^-1. Comparing this resolved star formation rate, to extragalactic star formation rate tracers (based on the Kennicutt-Schmidt relation), we find the star formation activity is underestimated by a factor of \geq 2 in comparison to the SFR derived from the YSO population., Comment: Accepted by MNRAS, 16 pages, 8 figures, 3 tables

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], Comment: 23 pages, 14 figures, accepted for publication in Astronomy and Astrophysics

Despite their importance to a number of astrophysical fields, the lifecycles of very massive stars are still poorly defined. In order to address this shortcoming, we present a detailed quantitative study of the physical properties of four early-B hypergiants (BHGs); Cyg OB2 #12, zeta Sco, HD190603 and BP Cru. These are combined with an analysis of their long-term spectroscopic and photometric behaviour in order to determine their evolutionary status. The long-term datasets revealed that they are remarkably stable over long periods (>40yr), with the possible exception of zeta Sco prior to the 20th century, in contrast to the typical excursions that characterise luminous blue variables (LBVs). Zeta Sco, HD190603 and BP Cru possess physical properties intermediate between B supergiants and LBVs; we therefore suggest that BHGs are the immediate descendants and progenitors (respectively) of such stars (for initial masses in the range ~30-60Msun). In contrast, while the wind properties of Cyg OB2 #12 are consistent with this hypothesis, the combination of extreme luminosity and spectroscopic mass (~110Msun) and comparatively low temperature means it cannot be accommodated in such a scheme. Likewise, despite its co-location with several LBVs above the Humphreys-Davidson (HD) limit, the lack of long term variability and its unevolved chemistry apparently excludes such an identification. Since such massive stars are not expected to evolve to such cool temperatures, the properties of Cyg OB2 #12 are difficult to understand under current evolutionary paradigms. [ABRIDGED], Comment: 36 pages, 19 figures (of which 17 pages are online supplemental material). Accepted for publication in Astronomy and Astrophysics

We present an analysis of a peculiar supergiant B-type star (VFTS698/Melnick 2/Parker 1797) in the 30 Doradus region of the Large Magellanic Cloud which exhibits characteristics similar to the broad class of B[e] stars. We analyse optical spectra from the VLT-FLAMES survey, together with archival optical and infrared photometry and X-ray imaging to characterise the system. We find radial velocity variations of around 400 km/s in the high excitation Si IV, N III and He II spectra, and photometric variability of ~0.6 mag with a period of 12.7 days. In addition, we detect long-term photometric variations of ~0.25 mag, which may be due to a longer-term variability with a period of ~400 days. We conclude that VFTS698 is likely an interacting binary comprising an early B-type star secondary orbiting a veiled, more massive companion. Spectral evidence suggests a mid-to-late B-type primary, but this may originate from an optically-thick accretion disc directly surrounding the primary., Comment: 18 pages, 14 figures and 8 tables. Table 8 to be published online

We present wide-area radio continuum 5.5 and 8.8 GHz (5.5 and 3.4 cm) Australia Telescope Compact Array observations of the complex and rich massive star-forming region G305. The aim of this study is to perform an un-targeted survey of the region in search of the compact radio emission associated with ultra-compact (UC) HII regions. Observations presented here encompass the entire complex and have a maximum resolution of ~1.5x1.4" and sensitivity of ~0.07 mJy beam^-1. By applying a data reduction method that emphasises small-scale structure, we are able to detect 71 compact radio sources distributed throughout the observed field. To explore the nature of these compact radio sources we compare to mid-infrared data and in this way identify 56 background sources, eight stellar radio sources, a single bright-rimmed cloud and six candidate UCHII regions. The physical properties of these candidate UCHII regions are determined and reveal five candidates have peak properties consistent with known UCHII regions with source radii ranging from 0.04-0.1 pc, emission measures from 2.56--10.3x10^-6 pc cm^-6 and electron densities of 0.34--1.03x10^4 cm^-3. We comment on these sites of recent massive star formation within G305 and by comparing to other star formation tracers (masers, NH3, YSOs) build a picture of the star formation history of the region. Using these results we estimate a lower limit to the star formation rate for the region of ~0.003 Msun yr^-1., Comment: 13 pages, 4 figures

NGC346-013 is a peculiar double-lined eclipsing binary in the Small Magellanic Cloud discovered by the VLT-FLAMES survey of massive stars. Spectra obtained with VLT-FLAMES are used to construct a radial velocity curve and photometry obtained with the Faulkes Telescope South is then used to derive orbital parameters, while spectra of the secondary are compared with synthetic spectra from TLUSTY model atmospheres. The orbital period is found to be 4.20381(12) days, with masses of 19.1+/-1.0 and 11.9+/-0.6 Msun. The primary is a rapidly rotating late-O dwarf while the secondary, an early-B giant, displays near-synchronous rotation and has filled its Roche lobe, implying that it was originally the more massive component with recent mass transfer `spinning up' the primary to near-critical rotation. Comparison with synthetic spectra finds temperatures of 34.5kK and 24.5kK for the primary and secondary respectively, with the nitrogen abundance of the secondary enhanced compared to baseline values for the SMC, consistent with the predictions of models of interacting binaries. NGC346-013 likely evolved via non-conservative mass transfer in a system with initial masses ~22+15Msun, with the well-constrained orbital solution and atmospheric parameters making it an excellent candidate for tailored modelling with binary evolution codes. This system will form a cornerstone in constraining the physics of thermal timescale mass transfer, and the associated mass transfer efficiency, in massive close binary systems., Comment: 9 pages, 5 figures, accepted for publication in Astronomy and Astrophysics

The G305 HII complex (G305.4+0.1) is one of the most massive star forming structures yet identified within the Galaxy. It is host to many massive stars at all stages of formation and evolution, from embedded molecular cores to post main-sequence stars. Here, we present a detailed near-infrared analysis of the two central star clusters Danks 1 and Danks 2, using HST+NICMOS imaging and VLT+ISAAC spectroscopy. We find that the spectro-photometric distance to the clusters is consistent with the kinematic distance to the G305 complex, an average of all measurements giving a distance of 3.8\pm0.6kpc. From analysis of the stellar populations and the pre-main-sequence stars we find that Danks 2 is the elder of the two clusters, with an age of 3^{+3}_{-1}Myr. Danks 1 is clearly younger with an age of 1.5^{+1.5}_{-0.5}Myr, and is dominated by three very luminous H-rich Wolf-Rayet stars which may have masses \geq100\msun. The two clusters have mass functions consistent with the Salpeter slope, and total cluster masses of 8000\pm1500\msun\ and 3000\pm800\msun\ for Danks 1 and Danks 2 respectively. Danks 1 is significantly the more compact cluster of the two, and is one of the densest clusters in the Galaxy with $\log (\rho/M_{\odot}{\rm pc}^{-3}) = 5.5^{+0.5}_{-0.4}$. In addition to the clusters, there is a population of apparently isolated Wolf-Rayet stars within the molecular cloud's cavity. Our results suggest that the star-forming history of G305 began with the formation of Danks 2, and subsequently Danks 1, with the origin of the diffuse evolved population currently uncertain. Together, the massive stars at the centre of the G305 region appear to be clearing away what is left of the natal cloud, triggering a further generation of star formation at the cloud's periphery., Comment: Accepted for publication in MNRAS. 18 pages, 9 figures, 5 tables

We report near-infrared observations of the supergiant donor to the eclipsing high mass X-ray binary pulsar IGR J18027-2016. We aim to determine its spectral type and measure its radial velocity curve and hence determine the stellar masses of the components. ESO/VLT observations of the donor utilising the NIR spectrograph ISAAC were obtained in the H and K bands. The multi-epoch H band spectra were cross-correlated with RV templates in order to determine a radial solution for the system. The spectral type of the donor was confirmed as B0-1 I. The radial velocity curve constructed has a semi-amplitude of $23.8 \pm 3.1$ km s$^{-1}$. Combined with other measured system parameters, a dynamically determined neutron star mass of 1.4 $\pm$ 0.2 - 1.6 $\pm$ 0.3 M$_{\odot}$ is found. The mass range of the B0-B1 I donor was 18.6 $\pm$ 0.8 - 21.8 $\pm$ 2.4 ~M$_{\odot}$. These lower and upper limits were obtained under the assumption that the system is viewed edge-on (i = 90$^\circ$ with $\beta$ = 0.89) for the lower limit and the donor fills its Roche lobe ($\beta = 1$ with i = 73.1$^\circ$) for the upper limit respectively., Comment: 5 pages, 4 figures, Accepted for publication in A&A

There is growing evidence that a treatment of binarity amongst OB stars is essential for a full theory of stellar evolution. However the binary properties of massive stars - frequency, mass ratio and orbital separation - are still poorly constrained. In order to address this shortcoming we have undertaken a multiepoch spectroscopic study of the stellar population of the young massive cluster Westerlund 1. In this paper we present an investigation into the nature of the dusty Wolf-Rayet star and candidate binary W239. We determined an orbital period of only 5.05 days - making it one of the most compact WC binaries yet identified - while we interpret a significant near-IR flare between 2004-6 as evidence for a third massive stellar component in the system in a wide eccentric orbit. The presence of a near-IR excess characteristic of hot dust at every epoch is consistent with the expectation that the subset of persistent dust forming WC stars are short period binaries. In conjunction with previous results, the simultaneous presence of both short period Wolf-Rayet binaries and cool hypergiants within Wd 1 provides compelling evidence for a bifurcation in the post-Main Sequence evolution of massive stars due to binarity. Short period O+OB binaries will evolve directly to the Wolf-Rayet phase, either due to an episode of binary mediated mass loss or via chemically homogenous evolution. Conversely, long period binaries and single stars will instead undergo a red loop across the HR diagram via a cool hypergiant phase. Future analysis of the full spectroscopic dataset for Wd 1 will constrain the proportion of massive stars experiencing each pathway; hence quantifying the importance of binarity in massive stellar evolution up to and beyond supernova and the resultant production of relativistic remnants and X-ray binaries. [ABRIDGED], Comment: 11 pages, 6 figures, accepted for publication in A&A