On the nature of the galactic early-B hypergiants

Aims. 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) of spectral type B1-4 Ia+; Cyg OB2 #12, ζ1 Sco, HD 190603 and BP Cru. These are combined with an analysis of their long-term spectroscopic and photometric behaviour in order to determine their evolutionary status. Methods. Quantitative analysis of UV–radio photometric and spectroscopic datasets was undertaken with a non-LTE model atmosphere code in order to derive physical parameters for comparison with apparently closely related objects, such as B supergiants (BSGs) and luminous blue variables (LBVs), and theoretical evolutionary predictions. Results. The long-term photospheric and spectroscopic datasets compiled for the early-B HGs revealed that they are remarkably stable over long periods ( ≥ 40 yrs), with the possible exception of ζ1 Sco prior to the 20th century; in contrast to the typical excursions that characterise LBVs. Quantitative analysis of ζ1 Sco, HD 190603 and BP Cru yielded physical properties intermediate between BSGs and LBVs; we therefore suggest that BHGs are the immediate descendants and progenitors (respectively) of such stars, for initial masses in the range ~30−60 M⊙. Comparison of the properties of ζ1 Sco with the stellar population of its host cluster/association NGC 6231/Sco OB1 provides further support for such an evolutionary scenario. In contrast, while the wind properties of Cyg OB2 #12 are consistent with this hypothesis, the combination of extreme luminosity and spectroscopic mass (~110 M⊙) 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, instead traversing an O4-6Ia → O4-6Ia+ → WN7-9ha pathway, the properties of Cyg OB2 #12 are therefore difficult to understand under current evolutionary paradigms. Finally, we note that as with AG Car in its cool phase, despite exceeding the HD limit, the properties of Cyg OB2 #12 imply that it lies below the Eddington limit – thus we conclude that the HD limit does not define a region of the HR diagram inherently inimical to the presence of massive stars.

The young open cluster Berkeley 55

We present UBV photometry of the highly reddened and poorly studied open cluster Berkeley 55, revealing an important population of B-type stars and several evolved stars of high luminosity. Intermediate-resolution far-red spectra of several candidate members confirm the presence of one F-type supergiant and six late supergiants or bright giants. The brightest blue stars are mid-B giants. Spectroscopic and photometric analyses indicate an age 50 ± 10 Myr. The cluster is located at a distance d ≈ 4 kpc, consistent with other tracers of the Perseus Arm in this direction. Berkeley 55 is thus a moderately young open cluster with a sizable population of candidate red (super)giant members, which can provide valuable information about the evolution of intermediate-mass stars.

A massive association around the obscured open cluster RSGC3

Context. Four clusters of red supergiants have been discovered in a region of the Milky Way close to base of the Scutum-Crux Arm and the tip of the Long Bar. Population synthesis models indicate that they must be very massive to harbour so many supergiants. If the clusters are physically connected, this Scutum Complex would be the largest and most massive star-forming region ever identified in the Milky Way. Aims. The spatial extent of one of these clusters, RSGC3, has not been investigated. In this paper we explore the possibility that a population of red supergiants could be located in its vicinity. Methods. We utilised 2MASS JHKS photometry to identify candidate obscured luminous red stars in the vicinity of RSGC3. We observed a sample of candidates with the TWIN spectrograph on the 3.5-m telescope at Calar Alto, obtaining intermediate-resolution spectroscopy in the 8000−9000 Å range. We re-evaluated a number of classification criteria proposed in the literature for this spectral range and found that we could use our spectra to derive spectral types and luminosity classes. Results. We measured the radial velocity of five members of RSGC3, finding velocities similar to the average for members of Stephenson 2. Among the candidates observed outside the cluster, our spectra revealed eight M-type supergiants at distances <18′ from the centre of RSGC3, distributed in two clumps. The southern clump is most likely another cluster of red supergiants, with reddening and age identical to RSGC3. From 2MASS photometry, we identified four likely supergiant members of the cluster in addition to the five spectroscopically observed. The northern clump may be a small cluster with similar parameters. Photometric analysis of the area around RSGC3 suggests the presence of a large (>30) population of red supergiants with similar colours. Conclusions. Our data suggest that the massive cluster RSGC3 is surrounded by an extended association, which may be very massive ( ≳ 105 M⊙). We also show that supergiants in the Scutum Complex may be characterised via a combination of 2MASS photometry and intermediate-to-high-resolution spectroscopy in the Z band.

A VLT/FLAMES survey for massive binaries in Westerlund 1. III. The WC9d binary W239 and implications for massive stellar evolution

Context. 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 & orbital separation – are still poorly constrained. Aims. 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. Methods. To accomplish this we have utilised our spectroscopic data in conjunction with multi-year optical and near-IR photometric observations in order to search for binary signatures. Comparison of these data to synthetic non-LTE model atmosphere spectra were used to derive the fundamental properties of the WC9 primary. Results. We found W239 to have an orbital period of only ~5.05 days, making it one of the most compact WC binaries yet identified. Analysis of the long term near-IR lightcurve reveals a significant flare between 2004-6. We interpret this as evidence for a third massive stellar component in the system in a long period (>6 yr), eccentric orbit, with dust production occuring at periastron leading to the flare. The presence of a near-IR excess characteristic of hot (~1300 K) dust at every epoch is consistent with the expectation that the subset of persistent dust forming WC stars are short (<1 yr) period binaries, although confirmation will require further observations. Non-LTE model atmosphere analysis of the spectrum reveals the physical properties of the WC9 component to be fully consistent with other Galactic examples. Conclusions. 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 – likely via case A mass transfer or a contact configuration – 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.

Circumstellar emission in Be/X-ray binaries of the Magellanic Clouds and the Milky Way

Aims. We study the optical and near-infrared colour excesses produced by circumstellar emission in a sample of Be/X-ray binaries. Our main goals are exploring whether previously published relations, valid for isolated Be stars, are applicable to Be/X-ray binaries and computing the distance to these systems after correcting for the effects of the circumstellar contamination. Methods. Simultaneous UBVRI photometry and spectra in the 3500−7000 Å spectral range were obtained for 11 optical counterparts to Be/X-ray binaries in the LMC, 5 in the SMC and 12 in the Milky Way. As a measure of the amount of circumstellar emission we used the Hα equivalent width corrected for photospheric absorption. Results. We find a linear relationship between the strength of the Hα emission line and the component of E(B − V) originating from the circumstellar disk. This relationship is valid for stars with emission lines weaker than EW ≈ −15   Å. Beyond this point, the circumstellar contribution to E(B − V) saturates at a value ≈0.17   mag. A similar relationship is found for the (V − I) near infrared colour excess, albeit with a steeper slope and saturation level. The circumstellar excess in (B − V) is found to be about five times higher for Be/X-ray binaries than for isolated Be stars with the same equivalent width EW(Hα), implying significant differences in the physical properties of their circumstellar envelopes. The distance to Be/X-ray binaries (with non-shell Be star companions) can only be correctly estimated by taking into account the excess emission in the V band produced by free-free and free-bound transitions in the circumstellar envelope. We provide a simple method to determine the distances that includes this effect.

On the nature of candidate luminous blue variables in M 33

Context. Luminous blue variables (LBVs) are a class of highly unstable stars that have been proposed to play a critical role in massive stellar evolution as well as being the progenitors of some of the most luminous supernovae known. However the physical processes underlying their characteristic instabilities are currently unknown. Aims. In order to provide observational constraints on this behaviour we have initiated a pilot study of the population of (candidate) LBVs in the Local Group galaxy M 33. Methods. To accomplish this we have obtained new spectra of 18 examples within M 33. These provide a baseline of ≥ 4 yr with respect to previous observations, which is well suited to identifying LBV outbursts. We also employed existing multi-epoch optical and mid-IR surveys of M 33 to further constrain the variability of the sample and search for the presence of dusty ejecta. Results. Combining the datasets reveals that spectroscopic and photometric variability appears common, although in the majority of cases further observations will be needed to distinguish between an origin for this behavour in short lived stochastic wind structure and low level photospheric pulsations or coherent long term LBV excursions. Of the known LBVs we report a hitherto unidentified excursion of M 33 Var C between 2001-5, while the transition of the WNLh star B517 to a cooler B supergiant phase between 1993−2010 implies an LBV classification. Proof-of-concept quantitative model atmosphere analysis is provided for Romano’s star; the resultant stellar parameters being consistent with the finding that the LBV excursions of this star are accompanied by changes in bolometric luminosity. The combination of temperature and luminosity of two stars, the BHG [HS80] 110A and the cool hypergiant B324, appear 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 η Car are identified. The combined dataset suggests that the criteria employed to identify candidate LBVs results in a heterogeneous sample, also containing 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 progenitor of SN 2008S and the 2008 NGC 300 transient (albeit suffering less intrinsic extinction). Conclusions. The results of such a multiwavelength observational approach, employing multiplexing spectrographs and supplemented with quantitative model atmosphere analysis, appears to show considerable promise in both identifying and characterising the physical properties of LBVs as well as other short lived phases of massive stellar evolution.

NGC 7419 as a template for red supergiant clusters

Context. The open cluster NGC 7419 is known to contain five red supergiants and a very high number of Be stars. However, there are conflicting reports about its age and distance that prevent a useful comparison with other clusters. Aims. We intend to obtain more accurate parameters for NGC 7419, using techniques different from those of previous authors, so that it may be used as a calibrator for more obscured clusters. Methods. We obtained Strömgren photometry of the open cluster NGC 7419, as well as classification spectroscopy of ~20 stars in the area. We then applied standard analysis and classification techniques. Results. We find a distance of 4 ± 0.4 kpc and an age of 14 ± 2 Myr for NGC 7419. The main-sequence turn-off is found at spectral type B1, in excellent agreement. We identify 179 B-type members, implying that there are more than 1200 M⊙ in B stars at present. Extrapolating this to lower masses indicates an initial cluster mass of between 7000 and 10 000 M⊙, depending on the shape of the initial mass function. We find a very high fraction (≈40%) of Be stars around the turn-off, but very few Be stars at lower masses. We also report for the first time a strong variability in the emission characteristics of Be stars. We verified that the parameters of the red supergiant members can be used to obtain accurate cluster parameters. Conclusions. NGC 7419 is sufficiently massive to serve as a testbed for theoretical predictions and as a template to compare more obscured clusters. The distribution of stars above the main-sequence turn-off is difficult to accommodate with current evolutionary tracks. Though the presence of five red supergiants is marginally consistent with theoretical expectations, the high number of Be stars and very low number of luminous evolved B stars hint at some unknown physical factor that is not considered in current synthesis models.

The supergiant B[e] star LHA 115-S 18 - binary and/or luminous blue variable?

Context. 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. Aims. In order to investigate both issues we have analysed the long term behaviour of the canonical sgB[e] star LHA 115-S 18. Methods. We employed the OGLE II-IV lightcurve to search for (a-)periodic variability and supplemented these data with new and historic spectroscopy. Results. In contrast to historical expectations for sgB[e] stars, S18 is both photometrically and spectroscopically highly variable. The lightcurve is characterised by rapid aperiodic ` aring' 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 He ii emission, X-ray emission implies the presence of an unseen binary companion. Conclusions. 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 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.

Massive open star clusters using the VVV survey III. A young massive cluster at the far edge of the Galactic bar

Context. Young massive clusters are key to map the Milky Way’s structure, and near-infrared large area sky surveys have contributed strongly to the discovery of new obscured massive stellar clusters. Aims. We present the third article in a series of papers focused on young and massive clusters discovered in the VVV survey. This article is dedicated to the physical characterization of VVV CL086, using part of its OB-stellar population. Methods. We physically characterized the cluster using JHKS near-infrared photometry from ESO public survey VVV images, using the VVV-SkZ pipeline, and near-infrared K-band spectroscopy, following the methodology presented in the first article of the series. Results. Individual distances for two observed stars indicate that the cluster is located at the far edge of the Galactic bar. These stars, which are probable cluster members from the statistically field-star decontaminated CMD, have spectral types between O9 and B0 V. According to our analysis, this young cluster (1.0 Myr < age < 5.0 Myr) is located at a distance of 11+5-6 kpc, and we estimate a lower limit for the cluster total mass of (2.8+1.6-1.4) · 103 M⊙. It is likely that the cluster contains even earlier and more massive stars.

VdBH 222: a starburst cluster in the inner Milky Way

Context. It has been suggested that the compact open cluster VdBH 222 is a young massive distant object. Aims. We set out to characterise VdBH 222 using a comprehensive set of multi-wavelength observations. Methods. We obtained multi-band optical (UBVR) and near-infrared (JHKS) 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. Results. 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 = 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 ≈ 12 Ma and a location not far from the position where some Galactic models place the far end of the Galactic bar. Conclusions. VdBH 222 is a young massive cluster with a likely mass >20 000 M⊙. 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.

MY Camelopardalis, a very massive merger progenitor

Context. The early-type binary MY Cam belongs to the young open cluster Alicante 1, embedded in Cam OB3. Aims. MY Cam consists of two early-O type main-sequence stars and shows a photometric modulation suggesting an orbital period slightly above one day. We intend to confirm this orbital period and derive orbital and stellar parameters. Methods. Timing analysis of a very exhaustive (4607 points) light curve indicates a period of 1.1754514 ± 0.0000015 d. High-resolution spectra and the cross-correlation technique implemented in the todcor program were used to derive radial velocities and obtain the corresponding radial velocity curves for MY Cam. Modelling with the stellar atmosphere code fastwind was used to obtain stellar parameters and create templates for cross-correlation. Stellar and orbital parameters were derived using the Wilson-Devinney code, such that a complete solution to the binary system could be described. Results. The determined masses of the primary and secondary stars in MY Cam are 37.7 ± 1.6 and 31.6 ± 1.4M⊙, respectively. The corresponding temperatures, derived from the model atmosphere fit, are 42 000 and 39 000 K, with the more massive component being hotter. Both stars are overfilling their Roche lobes, sharing a common envelope. Conclusions. MY Cam contains the most massive dwarf O-type stars found so far in an eclipsing binary. Both components are still on the main sequence, and probably not far from the zero-age main sequence. The system is a likely merger progenitor, owing to its very short period.

Orbital and Physical Properties of the σ Ori Aa, Ab, B Triple System

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

An XMM-Newton view of FeKα in high-mass X-ray binaries

We present a comprehensive analysis of the whole sample of available XMM-Newton observations of high-mass X-ray binaries (HMXBs) until August 2013, focusing on the FeKα emission line. This line is key to better understanding the physical properties of the material surrounding the X-ray source within a few stellar radii (the circumstellar medium). We collected observations from 46 HMXBs and detected FeKα in 21 of them. We used the standard classification of HMXBs to divide the sample into different groups. We find that (1) different classes of HMXBs display different qualitative behaviours in the FeKα spectral region. This is visible especially in SGXBs (showing ubiquitous Fe fluorescence but not recombination Fe lines) and in γ Cass analogues (showing both fluorescent and recombination Fe lines). (2) FeKα is centred at a mean value of 6.42 keV. Considering the instrumental and fits uncertainties, this value is compatible with ionization states that are lower than Fe xviii. (3) The flux of the continuum is well correlated with the flux of the line, as expected. Eclipse observations show that the Fe fluorescence emission comes from an extended region surrounding the X-ray source. (4) We observe an inverse correlation between the X-ray luminosity and the equivalent width of FeKα (EW). This phenomenon is known as the X-ray Baldwin effect. (5) FeKα is narrow (σline< 0.15 keV), reflecting that the reprocessing material does not move at high speeds. We attempt to explain the broadness of the line in terms of three possible broadening phenomena: line blending, Compton scattering, and Doppler shifts (with velocities of the reprocessing material V ~ 1000 km s-1). (6) The equivalent hydrogen column (NH) directly correlates to the EW of FeKα, displaying clear similarities to numerical simulations. It highlights the strong link between the absorbing and the fluorescent matter. (7) The observed NH in supergiant X-ray binaries (SGXBs) is in general higher than in supergiant fast X-ray transients (SFXTs). We suggest two possible explanations: different orbital configurations or a different interaction compact object – wind. (8) Finally, we analysed the sources IGR J16320-4751 and 4U 1700-37 in more detail, covering several orbital phases. The observed variation in NH between phases is compatible with the absorption produced by the wind of their optical companions. The results clearly point to a very important contribution of the donor’s wind in the FeKα emission and the absorption when the donor is a supergiant massive star.

The little-studied cluster Berkeley 90 I. LS III +46 11: a very massive O3.5 If* + O3.5 If* binary

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

Open cluster Dolidze 25: Stellar parameters and the metallicity in the Galactic anticentre

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