The characterisation of the extinction curve in the near infrared (NIR) is fundamental to analyse the structure and stellar population of the Galactic centre (GC), whose analysis is hampered by the extreme interstellar extinction ($A_V\sim 30$ mag) that varies on arc-second scales. Recent studies indicate that the behaviour of the extinction curve might be more complex than previously assumed, pointing towards a variation of the extinction curve as a function of wavelength. We aim at analysing the variations of the extinction index, $\alpha$, with wavelength, line-of-sight, and absolute extinction, extending previous analysis to a larger area of the innermost regions of the Galaxy. We analysed the whole GALACTICNUCLEUS survey, a high-angular resolution ($\sim 0.2''$) $JHK_s$ NIR survey specially designed to observe the GC in unprecedented detail. It covers a region of $\sim 6000$\,pc$^2$, comprising fields in the nuclear stellar disc, the inner bulge, and the transition region between them. We applied two independent methods based on red clump (RC) stars to constrain the extinction curve and analysed its variation superseding previous studies. We used more than 165,000 RC stars and increased significantly the size of the regions analysed to confirm that the extinction curve varies with the wavelength. We estimated a difference $\Delta\alpha = 0.21\pm0.07$ between the obtained extinction indices, $\alpha_{JH}=2.44\pm0.05$ and $\alpha_{HK_s} = 2.23\pm0.05$. We also concluded that there is no significant variation of the extinction curve with wavelength, with the line-of-sight or the absolute extinction. Finally, we computed the ratios between extinctions, $A_J/A_H = 1.87\pm0.03$ and $A_{H}/A_{K_s} = 1.84\pm0.03$, consistent with all the regions of the GALACTICNUCLEUS catalogue.

[Abridged] Context: Radiation-driven mass loss plays a key role in the life-cycles of massive stars. However, basic predictions of such mass loss still suffer from significant quantitative uncertainties. Aims: We develop new radiation-driven, steady-state wind models for massive stars with hot surfaces, suitable for quantitative predictions of global parameters like mass-loss and wind-momentum rates. Methods: The simulations presented here are based on a self-consistent, iterative grid-solution to the spherically symmetric, steady-state equation of motion, using full NLTE radiative transfer solutions in the co-moving frame to derive the radiative acceleration. We do not rely on any distribution functions or parametrization for computation of the line force responsible for the wind driving. Results: In this first paper, we present models representing two prototypical O-stars in the Galaxy, one with a higher stellar mass M/Msun=59 and luminosity log10(L/Lsun)= 5.87 and one with a lower M/Msun= 27 and log10(L/Lsun)= 5.1. For these simulations, basic predictions for global mass-loss rates and velocity laws are given, and the influence from additional parameters like wind clumping and microturbulent speeds discussed. A key result is that our mass-loss rates are significantly lower than those predicted by the mass-loss recipes normally included in models of massive-star evolution. Conclusions: Our results support previous suggestions that Galactic O-star mass-loss rates may be overestimated in present-day stellar evolution models, and that new rates thus might be needed. Indeed, future papers in this series will incorporate our new models into such simulations of stellar evolution, extending the very first simulations presented here toward larger grids covering a range of metallicities, B supergiants across the bistability jump, and possibly also Wolf-Rayet stars.

The nuclear disc is a dense stellar structure at the centre of the Milky Way, with a radius of $\sim$150 pc. It has been a place of intense star formation in the past several tens of millions of years but its overall formation history has remained unknown up to now. Here we report the first detailed star formation history of this region. The bulk of its stars formed at least eight billion years ago. This initial activity was followed by a long period of quiescence that was ended by an outstanding event about 1 Gyr ago, during which roughly 5% of its mass formed in a time window $\sim$100 Myr, in what may arguably have been one of the most energetic events in the history of the Milky Way. Star formation continued subsequently on a lower level, creating a few percent of the stellar mass in the past $\sim$500 Myr, with an increased rate up to $\sim$30 Myr ago. Our findings contradict the previously accepted paradigm of quasi-continuous star formation at the centre of the Milky Way. The long quiescent phase agrees with the overall quiescent history of the Milky Way and suggests that our Galaxy's bar may not have existed until recently, or that gas transport through the bar was extremely inefficient during a long stretch of the Milky Way's life, and that the central black hole may have acquired most of its mass already in the early days of the Milky Way.

The star S0-2, orbiting the Galactic central massive black hole candidate Sgr A*, passed its pericenter in May 2018. This event is the first chance to detect the general relativistic (GR) effect of a massive black hole, free from non-gravitational physics. The observable GR evidence in the event is the difference between the GR redshift and the Newtonian redshift of photons coming from S0-2. Within the present observational precision, the 1st post-Newtonian (1PN) GR evidence is detectable. In this paper, we give a theoretical analysis of the time evolution of the 1PN GR evidence, under a presupposition that is different from used in previous papers. Our presupposition is that the GR/Newtonian redshift is always calculated with the parameter values (the mass of Sgr A*, the initial conditions of S0-2, and so on) determined by fitting the GR/Newtonian motion of S0-2 with the observational data. It is then revealed that the difference of the GR redshift and the Newtonian one shows two peaks before and after the pericenter passage. This double-peak-appearance is due to our presupposition, and reduces to a single peak if the same parameter values are used in both GR and Newtonian redshifts as considered in previous papers. In addition to this theoretical discussion, we report our observational data obtained with the Subaru telescope by 2018. The quality and the number of Subaru data in 2018 are not sufficient to confirm the detection of the double-peak-appearance.

Due to the extreme extinction towards the Galactic centre ($A_{V} \sim 30$ mag), its stellar population is mainly studied in the near-infrared (NIR) regime. Therefore, a proper analysis of the NIR extinction curve is necessary to fully characterise the stellar structure and population of the inner part of the galaxy. We studied the dependence of the extinction index ($\alpha_\lambda$) in the NIR on the line of sight, wavelength, and extinction. We used the GALACTICNUCLEUS imaging survey, a high angular resolution catalogue ($0.2''$) for the inner part of the Galaxy in $JHK_s$, and studied the spatial variation in the extinction index. We also applied two independent methods based on red clump stars to compute the extinction index between different bands and its variation with wavelength. We did not detect any significant line-of-sight or extinction variation in $\alpha$ within the studied region in the nuclear stellar disc. The extinction index between $JH$ and $HK_s$ differs by $0.19 \pm 0.05$. We obtained mean values for the extinction indices $\alpha_{JH} = 2.43\pm0.03$ and $\alpha_{HK_s} = 2.23\pm0.03$. The dependence of the extinction index on the wavelength could explain the differences obtained for $\alpha_\lambda$ in the literature since it was assumed constant for the NIR regime.

The high extinction and extreme source crowding of the central regions of the Milky Way are serious obstacles to the study of the structure and stellar population of the Galactic centre (GC). Existing surveys that cover the GC region (2MASS, UKIDSS, VVV, SIRIUS) do not have the necessary high angular resolution. Therefore, a high-angular-resolution survey in the near infrared is crucial to improve the state of the art. Here, we present the GALACTICNUCLEUS catalogue, a near infrared $JHK_s$ high-angular-resolution ($0.2''$) survey of the nuclear bulge of the Milky Way. We explain in detail the data reduction, data analysis, calibration, and uncertainty estimation of the GALACTICNUCLEUS survey. We assess the data quality comparing our results with previous surveys. We obtained accurate $JHK_s$ photometry for $\sim 3.3\times10^6$ stars in the GC detecting around 20 \% in $J$, 65 \% in $H,$ and 90 \% in $K_s$. The survey covers a total area of $\sim0.3$ square degrees, which corresponds to $\sim 6,000$ pc$^2$. The GALACTICNUCLEUS survey reaches 5\,$\sigma$ detections for $J \sim 22$ mag, $H \sim 21$ mag, and $K_s \sim 21$ mag. The uncertainties are below 0.05 mag at $J \sim 21$ mag, $H \sim 19$ mag, and $K_s \sim 18$ mag. The zero point systematic uncertainty is $\lesssim0.04$ mag in all three bands. We present colour-magnitude diagrams for the different regions covered by the survey.

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

Recent advances in laboratory spectroscopy lead to the claim of ionized Buckminsterfullerene (C60+) as the carrier of two diffuse interstellar bands (DIBs) in the near-infrared. However, irrefutable identification of interstellar C60+ requires a match between the wavelengths and the expected strengths of all absorption features detectable in the laboratory and in space. Here we present Hubble Space Telescope (HST) spectra of the region covering the C60+ 9348, 9365, 9428 and 9577 {\AA} absorption bands toward seven heavily-reddened stars. We focus in particular on searching for the weaker laboratory C60+ bands, the very presence of which has been a matter for recent debate. Using the novel STIS-scanning technique to obtain ultra-high signal-to-noise spectra without contamination from telluric absorption that afflicted previous ground-based observations, we obtained reliable detections of the (weak) 9365, 9428 {\AA} and (strong) 9577 {\AA} C60+ bands. The band wavelengths and strength ratios are sufficiently similar to those determined in the latest laboratory experiments that we consider this the first robust identification of the 9428 {\AA} band, and a conclusive confirmation of interstellar C60+.

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

We briefy describe the results of a K-band spectroscopic survey of over 500 highly reddened point-like objects on sightlines toward the Central Molecular Zone (CMZ) of the Galaxy. The goal was to find stars with featureless or nearly featureless spectra suitable for near- and mid-infrared absorption spectroscopy of the Galactic center's interstellar gas on sightlines spread across the CMZ. Until recently only a few such stars have been known outside of very localized sightlines in the vicinity of the Quintuplet and Central clusters. We have used Spitzer Space Telescope (GLIMPSE) and 2MASS photometry to select promising candidates, and over the last ten years have been acquiring low-resolution K-band spectra of them. As expected, the vast majority are cool and/or highly reddened red giants with complex photospheric spectra unsuitable for measuring faint interstellar lines. Approximately ten percent of them, whose observations are reported here, have featureless or nearly featureless spectra. Although not evenly distributed in Galactic longitude, these stars are scattered across the CMZ. Many of them are luminous stars that are deeply embedded in warm dust cocoons, and have K-band continua rising steeply to longer wavelengths. A significant fraction of them are hot stars of a variety of spectral types, including at least five newly discovered Wolf-Rayet stars. All of them should be suitable for spectroscopy of interstellar absorption lines at infrared wavelengths greater than 3 microns and many are also suitable at shorter wavelengths.

With both nebular- and stellar-derived abundances of $\lesssim$ 1/10 Zsun and low foreground extinction, Sextans A is a prime candidate to replace the Small Magellanic Cloud as reservoir of metal-poor massive stars and reference to study the metal-poor Universe. We report the discovery of two early-O type stars in Sextans A, the earliest O-stars with metallicity < 1/7 Zsun known to date, and two additional O9 stars. Colour-excess estimates towards individual targets, enabled by spectral typing, manifest that internal reddening is neither uniform nor negligible. The four targets define a new region of star formation far from the optically-brightest centre of the galaxy and from its conspicuous HII shells, but not devoid of neutral hydrogen. In fact, we detect a spatial correlation between OB-stars and HI in Sextans A and other dIrr's that leads us to propose that the neutral phase may be fundamental to star formation in low-density environments. According to the existing evidence at least two of the targets formed in isolation, thus suggestive of an stochastic sampling of the initial mass function that would enable low-mass galaxies like Sextans A to form very massive stars. The discovery of these four stars provide spatially-resolved, spectroscopic confirmation of recent findings suggesting that dwarf galaxies can sustain star formation despite the low density of the gas phase.

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.

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

As a young massive cluster in the Central Molecular Zone, the Arches cluster is a valuable probe of the stellar Initial Mass Function (IMF) in the extreme Galactic Center environment. We use multi-epoch Hubble Space Telescope observations to obtain high-precision proper motion and photometric measurements of the cluster, calculating cluster membership probabilities for stars down to 1.8 M$_{\odot}$ between cluster radii of 0.25 pc -- 3.0 pc. We achieve a cluster sample with just ~8% field contamination, a significant improvement over photometrically-selected samples which are severely compromised by the differential extinction across the field. Combining this sample with K-band spectroscopy of 5 cluster members, we forward model the Arches cluster to simultaneously constrain its IMF and other properties (such as age and total mass) while accounting for observational uncertainties, completeness, mass segregation, and stellar multiplicity. We find that the Arches IMF is best described by a 1-segment power law that is significantly top-heavy: $\alpha$ = 1.80 $\pm$ 0.05 (stat) $\pm$ 0.06 (sys), where dN/dm $\propto$ m$^{-\alpha}$, though we cannot discount a 2-segment power law model with a high-mass slope only slightly shallower than local star forming regions ($\alpha$ = 2.04$^{+0.14}_{-0.19}$ $\pm$ 0.04) but with a break at 5.8$^{+3.2}_{-1.2}$ $\pm$ 0.02 M$_{\odot}$. In either case, the Arches IMF is significantly different than the standard IMF. Comparing the Arches to other young massive clusters in the Milky Way, we find tentative evidence for a systematically top-heavy IMF at the Galactic Center.

How very massive stars form is still an open question in astrophysics. VFTS682 is among the most massive stars known, with an inferred initial mass of $\sim$150 $M_\odot$ . It is located in 30 Doradus at a projected distance of 29 pc from the central cluster R136. Its apparent isolation led to two hypotheses: either it formed in relative isolation or it was ejected dynamically from the cluster. We investigate the kinematics of VFTS682 as obtained by Gaia and Hubble Space Telescope astrometry. We derive a projected velocity relative to the cluster of $38 \pm 17 \mathrm{km \ s^{-1}}$ (1$\sigma$ confidence interval). Although the error bars are substantial, two independent measures suggest that VFTS682 is a runaway ejected from the central cluster. This hypothesis is further supported by a variety of circumstantial clues. The central cluster is known to harbor other stars more massive than 150 $M_\odot$ of similar spectral type and recent astrometric studies on VFTS16 and VFTS72 provide direct evidence that the cluster can eject some of its most massive members, in agreement with theoretical predictions. If future data confirm the runaway nature, this would make VFTS682 the most massive runaway star known to date.

The study of the inner region of the Milky Way's bulge is hampered by high interstellar extinction and extreme source crowding. Sensitive high angular resolution near-infrared imaging is needed to study stellar populations in such a complex environment. We use the 0.2$''$ angular resolution $JHK_s$ data from the GALACTICNUCLEUS survey to study the stellar population within two $8.0'\times 3.4'$ fields, about 0.6$^\circ$ and 0.4$^\circ$ to the Galactic north of the Milky Way's centre and to compare it with one in the immediate surroundings of Sagittarius A*. We also characterise the extinction curve of the two fields. The average interstellar extinction to the outer and the inner field is $A_{K_s} \sim 1.20 \pm 0.08$ mag and $\sim 1.48 \pm 0.10$ mag, respectively. We present $K_{s}$ luminosity functions that are complete down to at least 2 mag below the red clump (RC). We detect a feature in the luminosity functions that is fainter than the RC by $0.80\pm0.03$ and $0.79\pm0.02$ mag, respectively, in the $K_s$-band. It runs parallel to the reddening vector. We identify the feature as the red giant branch bump. Fitting $\alpha$-enhanced BaSTI luminosity functions to our data, we find that a single old stellar population of $\sim12.8 \pm 0.6$ Gyr and $Z = 0.040 \pm 0.003$ provides the best fit. We obtain that the stellar population in the innermost bulge is old, similar to the one at larger distances from the Galactic plane, and that its metallicity increases down to distances as short as about 60 pc from the centre of the Milky Way. Comparing it with previous known values at larger latitudes ($|b|>2^\circ$), our results favour a flattening of the gradient at $|b|<2^\circ$. As a secondary result we obtain that the extinction index in the studied regions agrees the value of $\alpha = 2.30\pm0.08$, derived in Nogueras-Lara et al. 2018 for the very Galactic centre.

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

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

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

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.

Measurements in the infrared wavelength domain allow us to assess directly the physical state and energy balance of cool matter in space, thus enabling the detailed study of the various processes that govern the formation and early evolution of stars and planetary systems in galaxies over cosmic time. Previous infrared missions, from IRAS to Herschel, have revealed a great deal about the obscured Universe, but sensitivity has been limited because up to now it has not been possible to fly a telescope that is both large and cold. SPICA is a mission concept aimed at taking the next step in mid- and far-infrared observational capability by combining a large and cold telescope with instruments employing state-of-the-art ultra-sensitive detectors. The mission concept foresees a 2.5-meter diameter telescope cooled to below 8 K. With cooling provided by mechanical coolers instead of depending on a limited cryogen supply, the mission lifetime can extend significantly beyond the required three years. SPICA offers instrumentation with spectral resolving powers ranging from R ~50 through 11000 in the 17-230 $\mu$m domain as well as R~28.000 spectroscopy between 12 and 18 $\mu$m. Additionally SPICA will provide efficient 30-37 $\mu$m broad band mapping, and polarimetric imaging in the 100-350 $\mu$m range. SPICA will provide unprecedented spectroscopic sensitivity of ~5 x $10^{-20}$ W/m$^2$ (5$\sigma$/1hr) - at least two orders of magnitude improvement over what has been attained to date. With this exceptional leap in performance, new domains in infrared astronomy will become accessible, allowing us, for example, to unravel definitively galaxy evolution and metal production over cosmic time, to study dust formation and evolution from very early epochs onwards, and to trace the formation history of planetary systems.

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

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

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

IR spectroscopy in the range 12-230 micron with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes that govern the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope (JWST) and the new generation of Extremely Large Telescopes (ELTs) at shorter wavelengths and the Atacama Large Millimeter Array (ALMA) at longer wavelengths. SPICA, with its 2.5-m telescope actively-cooled to below 8K, will obtain the first spectroscopic determination, in the mid-IR rest-frame, of both the star-formation rate and black hole accretion rate histories of galaxies, reaching lookback times of 12 Gyr, for large statistically significant samples. Densities, temperatures, radiation fields and gas-phase metallicities will be measured in dust-obscured galaxies and active galactic nuclei (AGN), sampling a large range in mass and luminosity, from faint local dwarf galaxies to luminous quasars in the distant Universe. AGN and starburst feedback and feeding mechanisms in distant galaxies will be uncovered through detailed measurements of molecular and atomic line profiles. SPICA's large-area deep spectrophotometric surveys will provide mid-IR spectra and continuum fluxes for unbiased samples of tens of thousands of galaxies, out to redshifts of z~6. Furthermore, SPICA spectroscopy will uncover the most luminous galaxies in the first few hundred million years of the Universe, through their characteristic dust and molecular hydrogen features.

The physical processes driving the chemical evolution of galaxies in the last $\sim 11\, \rm{Gyr}$ cannot be understood without directly probing the dust-obscured phase of star-forming galaxies and active galactic nuclei. This phase, hidden to optical tracers, represents the bulk of star formation and black hole accretion activity in galaxies at $1 < z < 3$. Spectroscopic observations with a cryogenic infrared (IR) observatory like SPICA will be sensitive enough to peer through the dust-obscured regions of galaxies and access the rest-frame mid- to far-IR range in galaxies at high-$z$. This wavelength range contains a unique suite of spectral lines and dust features that serve as proxies for the abundances of heavy elements and the dust composition, providing tracers with a feeble response to both extinction and temperature. In this work, we investigate how SPICA observations could be exploited to understand key aspects in the chemical evolution of galaxies: the assembly of nearby galaxies based on the spatial distribution of heavy element abundances, the global content of metals in galaxies reaching the knee of the luminosity function up to $z \sim 3$, and the dust composition of galaxies at high-$z$. Possible synergies with facilities available in the late 2020s are also discussed.

The Galactic Centre is of fundamental astrophysical interest, but existing near-infrared surveys fall short to cover it adequately. Here we introduce the GALACTICNUCLEUS survey, a JHKs imaging survey of the centre of the Milky Way with a 0.2" angular resolution. We present the observations of Field 1 of our survey, centred approximately on SgrA* with an approximate size of 7.95 ' x 3.43 '. We describe the observational set-up and data reduction pipeline and discuss the quality of the data. Finally, we present some preliminary analysis of the data. The data were acquired with the near-infrared camera HAWK-I at the ESO VLT. Short readout times in combination with the speckle holography algorithm allowed us to produce final images with a stable, Gaussian PSF of 0.2" FWHM. Astrometric calibration is achieved via the VVV survey and photometric calibration is based on the SIRIUS/IRSF survey. The quality of the data is assessed by comparison between observations of the same field with different detectors of HAWK-I and at different times. We reach 5 sigma detection limits of approximately J = 22, H = 21, and Ks = 20. The photometric uncertainties are less than 0.05 at J < 20, H < 17 and Ks < 16. We distinguish five stellar populations in the colour-magnitude diagrams; three of them appear to belong to foreground spiral arms, and the other two correspond to a high- and a low-extinction star groups at the Galactic Centre. We use our data to analyse the near-infrared extinction curve and conclude that it can be described very well by a power-law with an index of JHKs = 2.31 +- 0.03. We do not find any evidence that this index depends on the position along the line-of-sight, or on the absolute value of the extinction. We produce extinction maps that show the clumpiness of the ISM at the Galactic Centre. Finally, we estimate that the majority of the stars have solar or super-solar metallicity.

During the next closest approach of the orbiting star S2/S0-2 to the Galactic supermassive black hole (SMBH), it is estimated that RV uncertainties of ~ 10 km/s allow us to detect post-Newtonian effects throughout 2018. To evaluate an achievable uncertainty in RV and its stability, we have carried out near-infrared, high resolution (R ~ 20,000) spectroscopic monitoring observations of S2 using the Subaru telescope and the near-infrared spectrograph IRCS from 2014 to 2016. The Br-gamma absorption lines are used to determine the RVs of S2. The RVs we obtained are 497 km/s, 877 km/s, and 1108 km/s in 2014, 2015, and 2016, respectively. The statistical uncertainties are derived using the jackknife analysis. The wavelength calibrations in our three-year monitoring are stable: short-term (hours to days) uncertainties in RVs are < 0.5 km/s, and a long-term (three years) uncertainty is 1.2 km/s. The uncertainties from different smoothing parameter, and from the partial exclusion of the spectra, are found to be a few km/s. The final results using the Br-gamma line are 497 +- 17 (stat.) +- 3 (sys.) km/s in 2014, 877 +- 15 (stat.) +- 4 (sys.) km/s in 2015, and 1108 +- 12 (stat.) +- 4 (sys.) km/s in 2016. When we use two He I lines at 2.113\mum in addition to Br-gamma, the mean RVs are 513 km/s and 1114 km/s for 2014 and 2016, respectively. The standard errors of the mean are 16.2 km/s (2014) and 5.4 km/s (2016), confirming the reliability of our measurements. The difference between the RVs estimated by Newtonian mechanics and general relativity will reach about 200 km/s near the next pericenter passage in 2018. Therefore our RV uncertainties of 13 - 17 km/s with Subaru enable us to detect the general relativistic effects in the RV measurements with more than 10 sigma in 2018.

We report the detection of a number of emission lines in the 1.0--2.4~$\mu$m spectra of four of the five bright infrared dust-embedded stars at the center of the Galactic center's Quintuplet Cluster. Spectroscopy of the central stars of these objects is hampered not only by the large interstellar extinction that obscures all objects in the Galactic center, but also by the large amounts of warm circumstellar dust surrounding each of the five. The pinwheel morphologies of the dust observed previously around two of them are indicative of Wolf-Rayet colliding wind binaries; however, infrared spectra of each of the five have until now revealed only dust continua steeply rising to long wavelengths and absorption lines and bands from interstellar gas and dust. The emission lines detected, from ionized carbon and from helium, are broad and confirm that the objects are dusty late-type carbon Wolf-Rayet stars.

Due to recent advances in laboratory spectroscopy, the first optical detection of a very large molecule has been claimed in the diffuse interstellar medium (ISM): C60+ (ionized Buckminsterfullerene). Confirming the presence of this molecule would have significant implications regarding the carbon budget and chemical complexity of the ISM. Here we present results from a new method for ultra-high signal-to-noise (S/N) spectroscopy of background stars in the near infrared (at wavelengths 0.9-1 micron), using the Hubble Space Telescope Imaging Spectrograph (STIS) in a previously untested `STIS scan' mode. The use of HST provides the crucial benefit of eliminating the need for error-prone telluric correction methods in the part of the spectrum where the C60+ bands lie, and terrestrial water vapor contamination is severe. Our STIS spectrum of the heavily-reddened B0 star BD63\,1964 reaches an unprecedented S/N for this instrument ($\sim600-800$), allowing the detection of the diffuse interstellar band (DIB) at 9577 \AA\ attributed to C60+ as well as new DIBs in the near-IR. Unfortunately, the presence of overlapping stellar lines, and the unexpected weakness of the C60+ bands in this sightline, prevents conclusive detection of the weaker C60+ bands. A probable correlation between the 9577 \AA\ DIB strength and interstellar radiation field is identified, which suggests that more strongly-irradiated interstellar sightlines will provide the optimal targets for future C60+ searches.

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

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

FS CMa stars are low-luminosity objects showing the B[e] phenomenon whose evolutionary origin is yet to be unraveled. Various binary-related hypotheses have been recently proposed, two of them involving the spiral-in evolution of the binary orbit. The latter occurs more often in dense stellar environments like young massive clusters (YMCs). Hence, a systematic study of FS CMa stars in YMCs would be crucial to find out how these objects are created. In YMCs, two FS CMa stars have been confirmed and three candidates have been found through a search method based on narrow-band photometry at Paschen-$\alpha$ and the neighboring continuum. We apply this method to archival data from the Paschen-$\alpha$ survey of the Galactic Center region, yielding a new candidate in the Quintuplet cluster. Limitations of this method and other alternatives are briefly discussed.

Theoretically, rotation-induced chemical mixing in massive stars has far reaching evolutionary consequences, affecting the sequence of morphological phases, lifetimes, nucleosynthesis, and supernova characteristics. Using a sample of 72 presumably single O-type giants to supergiants observed in the context of the VLT-FLAMES Tarantula Survey (VFTS), we aim to investigate rotational mixing in evolved core-hydrogen burning stars initially more massive than $15\,M_\odot$ by analysing their surface nitrogen abundances. Using stellar and wind properties derived in a previous VFTS study, we constrained the nitrogen abundance by fitting the equivalent widths of relatively strong lines that are sensitive to changes in the abundance of this element. Given the quality of the data, we constrained the nitrogen abundance in 38 cases; for 34 stars only upper limits could be derived, which includes almost all stars rotating at $v_\mathrm{e}\sin i >200\,\mathrm{km s^{-1}}$. We analysed the nitrogen abundance as a function of projected rotation rate $v_\mathrm{e}\sin i$ and confronted it with predictions of rotational mixing. The upper limits on the nitrogen abundance of the rapidly rotating stars are not in apparent violation with theoretical expectations. However, we found a group of N-enhanced slowly-spinning stars that is not in accordance with predictions of rotational mixing in single stars. Among O-type stars with (rotation-corrected) gravities less than $\log\,g_c = 3.75$ this group constitutes 30$-$40 percent of the population. We found a correlation between nitrogen and helium abundance which is consistent with expectations, suggesting that, whatever the mechanism that brings N to the surface, it displays CNO-processed material.

The Dragonfish Nebula has been recently claimed to be powered by a superluminous but elusive OB association. Instead, systematic searches in near-infrared photometric surveys have found many other cluster candidates on this sky region. Among these, the first confirmed young massive cluster was Mercer 30, where Wolf-Rayet stars were found. We perform a new characterization of Mercer 30 with unprecedented accuracy, combining NICMOS/HST and VVV photometric data with multi-epoch ISAAC/VLT H- and K-band spectra. Stellar parameters for most of spectroscopically observed cluster members are found through precise non-LTE atmosphere modeling with the CMFGEN code. Our spectrophotometric study for this cluster yields a new, revised distance of d = (12.4 +- 1.7) kpc and a total of Q = 6.70 x 10^50 Lyman ionizing photons. A cluster age of (4.0 +- 0.8) Myr is found through isochrone fitting, and a total mass of (1.6 +- 0.6) x 10^4 Msol is estimated thanks to our extensive knowledge of the post-main-sequence population. As a consequence, membership of Mercer 30 to the Dragonfish star-forming complex is confirmed, allowing us to use this cluster as a probe for the whole complex, which turns out to be extremely large (400 pc across) and located at the outer edge of the Sagittarius-Carina spiral arm (11 kpc from the Galactic Center). The Dragonfish complex hosts 19 young clusters or cluster candidates (including Mercer 30 and a new candidate presented in this work) and an estimated minimum of 9 field Wolf-Rayet stars. The sum of all these contributions accounts for, at least, 73% of the Dragonfish Nebula ionization and leaves little or no room for the alleged superluminous OB association; alternative explanations are discussed.

Rich in HII regions, giant molecular clouds are natural laboratories to study massive stars and sequential star formation. The Galactic star forming complex W33 is located at l=~12.8deg and at a distance of 2.4 kpc, has a size of ~10 pc and a total mass of (~0.8 - ~8.0) X 10^5 Msun. The integrated radio and IR luminosity of W33 - when combined with the direct detection of methanol masers, the protostellar object W33A, and protocluster embedded within the radio source W33 main - mark the region out as a site of vigorous ongoing star formation. In order to assess the long term star formation history, we performed an infrared spectroscopic search for massive stars, detecting for the first time fourteen early-type stars, including one WN6 star and four O4-7 stars. The distribution of spectral types suggests that this population formed during the last ~2-4 Myr, while the absence of red supergiants precludes extensive star formation at ages 6-30 Myr. This activity appears distributed throughout the region and does not appear to have yielded the dense stellar clusters that characterize other star forming complexes such as Carina and G305. Instead, we anticipate that W33 will eventually evolve into a loose stellar aggregate, with Cyg OB2 serving as a useful, albeit richer and more massive, comparator. Given recent distance estimates, and despite a remarkably similar stellar population, the rich cluster Cl 1813-178 located on the north-west edge of W33 does not appear to be physically associated with W33.

FS CMa stars are low-luminosity objects showing the B[e] phenomenon whose evolutionary state remains a puzzle. These stars are surrounded by compact disks of warm dust of unknown origin. Hitherto, membership of FS CMa stars to coeval populations has never been confirmed. The discovery of low-luminosity line emitters in the young massive clusters Mercer 20 and Mercer 70 prompts us to investigate the nature of such objects. We intend to confirm membership to coeval populations in order to characterize these emission-line stars through the cluster properties. Based on ISAAC/VLT medium-resolution spectroscopy and NICMOS/HST photometry of massive cluster members, new characterizations of Mercer 20 and Mercer 70 are performed. Coevality of each cluster and membership of the newly-discovered B[e] objects are investigated using our observations as well as literature data of the surroundings. Infrared excess and narrow-band photometric properties of the B[e] stars are also studied. We confirm and classify 22 new cluster members, including Wolf-Rayet stars and blue hypergiants. Spectral types (O9-B1.5 V) and radial velocities of B[e] objects are compatible with the remaining cluster members, while emission features of Mg II, Fe II], and [Fe II] are identified in their spectra. The ages of these stars are 4.5 and 6 Myr, and they show mild infrared excesses. We confirm the presence of FS CMa stars in the coeval populations of Mercer 20 and Mercer 70. We discuss the nature and evolutionary state of FS CMa stars, discarding a post-AGB nature and introducing a new hypothesis about mergers. A new search method for FS CMa candidates in young massive clusters based on narrow-band Paschen-alpha photometry is proposed and tested in photometric data of other clusters, yielding three new candidates.

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

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

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

The evolution and fate of very massive stars (VMS) is tightly connected to their mass-loss properties. Their initial and final masses differ significantly as a result of mass loss. VMS have strong stellar winds and extremely high ionising fluxes, which are thought to be critical sources of both mechanical and radiative feedback in giant Hii regions. However, how VMS mass-loss properties change during stellar evolution is poorly understood. In the framework of the VLT-Flames Tarantula Survey (VFTS), we explore the mass-loss transition region from optically thin O to denser WNh star winds, thereby testing theoretical predictions. To this purpose we select 62 O, Of, Of/WN, and WNh stars, an unprecedented sample of stars with the highest masses and luminosities known. We perform a spectral analysis of optical VFTS as well as near-infrared VLT/SINFONI data using the non-LTE radiative transfer code CMFGEN to obtain stellar and wind parameters. For the first time, we observationally resolve the transition between optically thin O and optically thick WNh star winds. Our results suggest the existence of a kink between both mass-loss regimes, in agreement with recent MC simulations. For the optically thick regime, we confirm the steep dependence on the Eddington factor from previous theoretical and observational studies. The transition occurs on the MS near a luminosity of 10^6.1Lsun, or a mass of 80...90Msun. Above this limit, we find that - even when accounting for moderate wind clumping (with f = 0.1) - wind mass-loss rates are enhanced with respect to standard prescriptions currently adopted in stellar evolution calculations. We also show that this results in substantial helium surface enrichment. Based on our spectroscopic analyses, we are able to provide the most accurate ionising fluxes for VMS known to date, confirming the pivotal role of VMS in ionising and shaping their environments.

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.

We present the first quantitative UV spectroscopic analysis of resolved OB stars in IC1613. Because of its alleged very low metallicity (<~1/10 Zo, from HII regions), studies in this Local Group dwarf galaxy could become a significant step forward from the SMC towards the extremely metal-poor massive stars of the early Universe. We present HST-COS data covering the ~1150-1800{\AA} wavelength range with resolution R~2500. We find that the targets do exhibit wind features, and these are similar in strength to SMC stars. Wind terminal velocities were derived from the observed PCygni profiles with the SEI method. The vinf-Z relationship has been revisited. The terminal velocity of IC1613 O-stars is clearly lower than Milky Way counterparts, but there is no clear difference between IC1613 and SMC or LMC analogue stars. We find no clear segregation with host galaxy in the terminal velocities of B-supergiants, nor in the vinf/vesc ratio of the whole OB star sample in any of the studied galaxies. Finally, we present first evidence that the Fe-abundance of IC1613 OB stars is similar to the SMC, in agreement with previous results on red supergiants. With the confirmed ~1/10 solar oxygen abundances of B-supergiants, our results indicate that IC1613's [alpha/Fe] ratio is sub-solar.

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

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

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

We compile the first comprehensive census of hot luminous stars in the 30 Doradus (30 Dor) star forming region of the LMC. The census extends to a radius of 10arcmin (150pc) from the central cluster, R136. Stars were selected photometrically and combined with the latest spectral types. 1145 candidate hot luminous stars were identified of which >700 were considered genuine early type stars that contribute to feedback. We assess the spectroscopic completeness to be 85% in outer regions (>5pc) but fall to 35% in the vicinity of R136, giving a total of 500 hot luminous stars with spectroscopy. Stellar calibrations and models were used to obtain their physical parameters before integrated values were compared to global observations and the population synthesis code, Starburst99. The 31 W-R and Of/WN stars made large contributions to the total ionising and wind luminosities of ~40% and ~50%, respectively. Stars with Minit>100Msun also showed high contributions to the global feedback, ~25% in both cases. Such massive stars are not accounted for by the current Starburst99 code, which underestimated the ionising and wind luminosities of R136 by factors of ~2 and ~9, respectively. The census inferred a SFR of 0.073+/-0.04Msun/yr for 30 Dor, typically higher than results from popular SFR calibrations. However, it remained consistent with a far-UV luminosity tracer and a combined Halpha and mid-infrared tracer, but only after correcting for Halpha extinction. The global ionising output exceeded measurements from the associated gas and dust, suggesting ~6(+55/-6)% of ionising photons escape the region. When studying very luminous star forming regions, it is therefore essential to include the most massive stars to ensure a reliable energy budget. If 30 Dor is typical of other large star forming regions, estimates of the SFR will be underpredicted if this escape fraction is not accounted for.(abridged)

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

The development of infrared observational facilities has revealed a number of massive stars in obscured environments throughout the Milky Way and beyond. The determination of their stellar and wind properties from infrared diagnostics is thus required to take full advantage of the wealth of observations available in the near and mid infrared. However, the task is challenging. This session addressed some of the problems encountered and showed the limitations and successes of infrared studies of massive stars.

During the last decade, hundreds of young massive cluster candidates have been detected in the disk of the Milky Way. We investigate one of these candidates, Mercer 81, which was discovered through a systematic search for stellar overdensities, with follow-up NICMOS/HST infrared narrow-band photometry to find emission-line stars and confirm it as a massive cluster. Surprisingly, the brightest stars turned out to be a chance alignment of foreground stars, while a real massive cluster was found among some fainter stars in the field. From a first spectroscopic study of four emission-line stars (ISAAC/VLT), it follows that Mercer 81 is a very massive young cluster, placed at the far end of the Galactic bar. Additionally, in this work we present some unpublished spectra from a follow-up observation program which confirm that the cluster hosts several Nitrogen-rich Wolf-Rayet stars (WN) and blue supergiants.

Context. Results from the theory of radiatively driven winds are incorporated in stellar evolutionary and population synthesis models, and used in our interpretation of the observations of the deep Universe. Yet, the theory has been confirmed only until Small Magellanic Cloud metallicities. Analyses of O-stars at lower metallicities are needed to prove the theory. Aims. We have observed GHV-62024, an O6.5 IIIf star in the low-metallicity galaxy IC1613. According to a previous preliminary analysis this star could challenge the radiatively driven wind theory at low metallicities. Methods. Our observations were obtained with VIMOS at VLT, at R~2000 and were analysed using the latest version of the model atmosphere code FASTWIND, which includes N III Results. We obtain the stellar parameters and conclude that the star follows the average wind momentum-luminosity relationship (WLR) expected for its metallicity, but with a high value for the exponent of the wind velocity law, beta. We suggest that this high value may be reached because GHV-62024 could be a fast rotator seen at a low inclination angle. While the derived beta value does not change by adopting a lower wind terminal velocity, a wrong $V_\infty$ has a clear impact on the position of the star in the WLR diagram. The N and He abundances are very high, consistent with strong CNO mixing that could have been caused by the fast rotation, although we cannot discard a different origin. We find again the well-known mass-discrepancy. Conclusions. We conclude that the star follows the WLR expected for its metallicity. The results are consistent with GHV-62024 being a fast rotator seen close to pole-on, contaminated at the surface with CNO products and with a wind structure altered by the fast rotation without modifying the global WLR. We suggest that this could be a general property of fast rotators.

This is the third paper in a series aiming at the analysis of nitrogen abundances in O-type stars, to enable further constraints on the early evolution of massive stars. We provide first theoretical predictions for the NIV4058/NIII4640 emission line ratio in dependence of various parameters, and confront them with results from the analysis of a sample of early LMC/SMC O-stars. Stellar and wind parameters are determined by line profile fitting of H/He/N lines, exploiting the helium and nitrogen ionization balance. Corresponding synthetic spectra are calculated using the NLTE atmospheric code FASTWIND. Though there is a monotonic relationship between the emission line ratio and Teff, all other parameters being equal, theoretical predictions indicate additional dependencies, most notably, on the nitrogen abundance. These basic predictions are confirmed by results from atmospheric code CMFGEN. The effective temperatures for the earliest O-stars, inferred from the nitrogen ionization balance, are partly considerably hotter than indicated by previous studies. Consistent with earlier results, effective temperatures increase from supergiants to dwarfs for all spectral types in the LMC. The relation between observed NIV4058/NIII4640 emission line ratio and Teff, for a given luminosity class, turned out to be quite monotonic for our sample stars, and fairly consistent with our model predictions. The scatter within a spectral sub-type is mainly produced by abundance effects. Our findings suggest that the Walborn et al. (2002) classification scheme is able to provide a meaningful relation between spectral type and Teff, provided that it is possible to discriminate for the luminosity class. This might be difficult to achieve in low-Z environments such as the SMC, owing to rather low wind-strengths. According to our predictions, the major bias of the classification scheme is due to nitrogen content.

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]

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]

We present a near-infrared study of the candidate star cluster Mercer 81, located at the centre of the G338.4+0.1 HII region, and close to the TeV gamma-ray source HESS 1640-465. Using HST/NICMOS imaging and VLT/ISAAC spectroscopy we have detected a compact and highly extincted cluster of stars, though the bright stars in the centre of the field are in fact foreground objects. The cluster contains nine stars with strong Paschen-alpha emission, one of which we identify as a Wolf-Rayet (WR) star, as well as an A-type supergiant. The line-of-sight extinction is very large, $A_{V}\sim 45$, illustrating the challenges of locating young star clusters in the Galactic Plane. From a quantitative analysis of the WR star we argue for a cluster age of 3.7$^{+0.4}_{-0.5}$\,Myr, and, assuming that all emission-line stars are WRs, a cluster mass of $\ga 10^4$\msun. A kinematic analysis of the cluster's surrounding HII-region shows that the cluster is located in the Galactic disk at a distance of 11$\pm$2\,kpc. This places the cluster close to where the far end of the Bar intersects the Norma spiral arm. This cluster, as well as the nearby cluster [DBS2003]179, represent the first detections of active star cluster formation at this side of the Bar, in contrast to the near side which is well known to have recently undergone a $\sim 10^6$\msun\ starburst episode.

The spectrum of any star viewed through a sufficient quantity of diffuse interstellar material reveals a number of absorption features collectively called diffuse interstellar bands (DIBs). The first DIBs were reported 90 years ago, and currently well over 500 are known. None of them has been convincingly identified with any specific element or molecule, although recent studies suggest that the DIB carriers are polyatomic molecules containing carbon. Most of the DIBs currently known are at visible and very near-infrared wavelengths, with only two previously known at wavelengths beyond one micron (10,000 Angstroms), the longer of which is at 1.318 microns. Here we report the discovery of thirteen diffuse interstellar bands in the 1.5-1.8 micron interval on high extinction sightlines toward stars in the Galactic centre. We argue that they originate almost entirely in the Galactic Centre region, a considerably warmer and harsher environment than where DIBs have been observed previously. The relative strengths of these DIBs toward the Galactic Centre and the Cygnus OB2 diffuse cloud are consistent with their strengths scaling mainly with extinction by diffuse material.

This is the second paper in a series aiming at the analysis of nitrogen abundances in O-type stars, to enable further constraints on the early evolution of massive stars. We investigate the NIV lambda4058 emission line formation, provide nitrogen abundances for a substantial O-star sample in the LMC, and compare our findings with recent predictions from stellar evolutionary models. Stellar and wind parameters are determined by line profile fitting of hydrogen, helium and nitrogen lines, based on synthetic spectra calculated by FASTWIND. We derive nitrogen abundances for 20 O- and 5 B-stars, by analyzing all nitrogen lines present in the available optical spectra. The dominating process responsible for emission at NIV lambda4058 in O-stars is the strong depopulation of the lower level of the transition, which increases as a function of Mdot. Unlike the NIII triplet emission, resonance lines do not play a role for typical mass-loss rates and below. The bulk of our sample O-stars seems to be strongly nitrogen-enriched, and a clear correlation of nitrogen and helium enrichment is found. By comparing the nitrogen abundances as a function of vsini ('Hunter-plot') with tailored evolutionary calculations, we identify a considerable number of highly enriched objects at low rotation. Due to the low initial abundance, the detection of strong Nitrogen enrichment in the bulk of O-stars indicates that efficient mixing takes place already during the very early phases of stellar evolution of LMC O-stars.

This is the first paper in a series dealing with optical Nitrogen spectroscopy of O-type stars, aiming at the analysis of Nitrogen abundances. We implemented a new Nitrogen model atom into the NLTE atmosphere/spectrum synthesis code FASTWIND, and compare the resulting optical NIII lines at 4634/40/42 A with other predictions, mostly from Mihalas & Hummer (1973, ApJ 179, 827,`MH'), and from the alternative code CMFGEN. Using similar model atmospheres as MH (not blanketed and wind-free), we are able to reproduce their results, in particular the triplet emission lines. According to MH, these should be strongly related to dielectronic recombination (DR) and the drain by certain two-electron transitions. However, using realistic, fully line-blanketed atmospheres at solar abundances, the key role of DR controlling these emission features is superseded -- for O-star conditions -- by the strength of the stellar wind and metallicity. In the case of wind-free models, the resulting lower ionizing EUV-fluxes severely suppress the emission. As the mass-loss rate is increased, pumping through the NIII resonance line(s) in the presence of a near-photospheric velocity field results in a net optical triplet line emission. A comparison with results from CMFGEN is mostly satisfactory, except for the range 30 kK < Teff < 35 kK, where CMFGEN triggers the triplet emission at lower Teff than FASTWIND. This effect could be traced down to line overlap effects between the NIII and OIII resonance lines that so far cannot be simulated by FASTWIND. Since the efficiency of DR and `two electron drain' strongly depends on the degree of line-blanketing/-blocking, we predict the emission to become stronger in a metal-poor environment, though lower wind-strengths and Nitrogen abundances might counteract this effect. Weak winded stars should display less triplet emission than stars with `normal' winds.

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.

Context: The discovery of new massive star clusters and massive stellar populations in previously known clusters in our Galaxy by means of infrared studies has changed our view of the Milky Way from an inactive to an active star-forming machine. Within this scenario, we present a near-infrared spectrophotometric study of the stellar content of the compact H II region Sh2-152. Aims: We aim to determine the distance, extinction, age, and mass of Sh2-152, using for the first time near-infrared stellar classification for several sources in the region. Methods: Using our near-infrared (J, H, and Ks) photometry and the colour-magnitude diagram for the cluster field, we selected 13 bright stars, candidate members of the reddened cluster's main sequence, for H- and K- spectroscopy and spectral classification. This near-infrared information was complemented with an optical spectrum of the ionizing central star to confirm its spectral nature. Results: From the 13 spectroscopically observed stars, 5 were classified as B-type, 3 as G-type, 2 were young stellar objects (YSOs), and 3 remained unclassified (because of the poor data quality). The cluster's extinction varies from A_{Ks}=0.5 to 2.6 magnitudes (A_V=4.5 to 24 magnitudes) and the distance is estimated to be 3.21\pm0.21 kpc. The age of the cluster is younger than 9.4 Myr and the lower limit to the total mass of the cluster is (2.45\pm0.79)\cdot10^3 M_{Sun}. We compare the number of ionizing photons emitted from the OB-type stars with the Lyman continuum photons derived from the radio observations and conclude that both quantities are consistent for the central region of Sh2-152. In contrast, the main ionizing source of the lower region remains unidentified.

Context. Mass-loss, occurring through radiation driven supersonic winds, is a key issue throughout the evolution of massive stars. Two outstanding problems are currently challenging the theory of radiation-driven winds: wind clumping and the weak-wind problem. Aims. We seek to obtain accurate mass-loss rates of OB stars at different evolutionary stages to constrain the impact of both problems in our current understanding of massive star winds. Methods. We perform a multi-wavelength quantitative analysis of a sample of ten Galactic OB-stars by means of the atmospheric code CMFGEN, with special emphasis on the L-band window. A detailed investigation is carried out on the potential of Br\alpha\ and Pf\gamma\ as mass-loss and clumping diagnostics. Results. For objects with dense winds, Br\alpha\ samples the intermediate wind while Pf\gamma\ maps the inner one. In combination with other indicators (UV, H\alpha, Br\gamma) these lines enable us to constrain the wind clumping structure and to obtain "true" mass-loss rates. For objects with weak winds, Br\alpha\ emerges as a reliable diagnostic tool to constrain the mass-loss rates. The emission component at the line Doppler-core superimposed on the rather shallow Stark absorption wings reacts very sensitively to mass loss already at very low mass-loss values. On the other hand, the line wings display similar sensitivity to mass loss as H\alpha, the classical optical mass loss diagnostics. Conclusions. Our investigation reveals the great diagnostic potential of L-band spectroscopy to derive clumping properties and mass-loss rates of hot star winds. We are confident that Br\alpha\ will become the primary diagnostic tool to measure very low mass-loss rates with unprecedented accuracy

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

The VLT-FLAMES Tarantula Survey (VFTS) is an ESO Large Programme that has obtained multi-epoch optical spectroscopy of over 800 massive stars in the 30 Doradus region of the Large Magellanic Cloud (LMC). Here we introduce our scientific motivations and give an overview of the survey targets, including optical and near-infrared photometry and comprehensive details of the data reduction. One of the principal objectives was to detect massive binary systems via variations in their radial velocities, thus shaping the multi-epoch observing strategy. Spectral classifications are given for the massive emission-line stars observed by the survey, including the discovery of a new Wolf-Rayet star (VFTS 682, classified as WN5h), 2' to the northeast of R136. To illustrate the diversity of objects encompassed by the survey, we investigate the spectral properties of sixteen targets identified by Gruendl & Chu from Spitzer photometry as candidate young stellar objects or stars with notable mid-infrared excesses. Detailed spectral classification and quantitative analysis of the O- and B-type stars in the VFTS sample, paying particular attention to the effects of rotational mixing and binarity, will be presented in a series of future articles to address fundamental questions in both stellar and cluster evolution.

Young massive (M >10^4 Msun) stellar clusters are a good laboratory to study the evolution of massive stars. Only a dozen of such clusters are known in the Galaxy. Here we report about a new young massive stellar cluster in the Milky Way. Near-infrared medium-resolution spectroscopy with UIST on the UKIRT telescope and NIRSPEC on the Keck telescope, and X-ray observations with the Chandra and XMM satellites, of the Cl 1813-178 cluster confirm a large number of massive stars. We detected 1 red supergiant, 2 Wolf-Rayet stars, 1 candidate luminous blue variable, 2 OIf, and 19 OB stars. Among the latter, twelve are likely supergiants, four giants, and the faintest three dwarf stars. We detected post-main sequence stars with masses between 25 and 100 Msun. A population with age of 4-4.5 Myr and a mass of ~10000 Msun can reproduce such a mixture of massive evolved stars. This massive stellar cluster is the first detection of a cluster in the W33 complex. Six supernova remnants and several other candidate clusters are found in the direction of the same complex.

The Tarantula survey is an ESO Large Programme which has obtained multi-epochs spectroscopy of over 800 massive stars in the 30 Dor region in the Large Magelanic Cloud. Here we briefly describe the main drivers of the survey and the observational material derived.

Recent IR surveys of the Galactic plane have revealed a large number of candidate Luminous Blue Variables. In order to verify these classifications we have been undertaking a long term spectroscopic and photometric monitoring campaign supplemented with tailored non-LTE model atmosphere analysis. Here we present a brief overview of selected aspects of this program, highlighting the prospects for identification, classification and quantitative analysis of LBVs in the near-IR spectral window.

The theory of radiatively driven winds successfully explains the key points of the stellar winds of hot massive stars. However, there is an apparent break-down of this paradigm at L/Lsun<5.2: the stellar wind momentum is smaller than predicted for low luminosity early-type stars from metal poor environments, and there are also some Galactic examples. In this work we explore whether magnetic fields are playing a role.

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

We report on a successful, simultaneous observation and modelling of the millimeter (mm) to near-infrared (NIR) flare emission of the Sgr A* counterpart associated with the supermassive black hole at the Galactic centre (GC). We present a mm/sub-mm light curve of Sgr A* with one of the highest quality continuous time coverages and study and model the physical processes giving rise to the variable emission of Sgr A*.

The characterisation of the stellar population toward young high-mass star-forming regions allows to constrain fundamental cluster properties like distance and age. These are essential when using high-mass clusters as probes to conduct Galactic studies. NGC 7538 is a star-forming region with an embedded stellar population only unearthed in the near-infrared. We present the first near-infrared spectro-photometric study of the candidate high-mass stellar content in NGC 7538. We obtained H and K spectra of 21 sources with both the multi-object and long-slit modes of LIRIS at the WHT, and complement these data with sub-arcsecond JHKs photometry of the region using the imaging mode of the same instrument. We find a wide variety of objects within the studied stellar population of NGC 7538. Our results discriminate between a stellar population associated to the HII region, but not contained within its extent, and several pockets of more recent star formation. We report the detection of CO bandhead emission toward several sources as well as other features indicative of a young stellar nature. We infer a spectro-photometric distance of 2.7+-0.5 kpc, an age spread in the range 0.5-2.2 Myr and a total mass ~1.7x10^3 Msun for the older population.

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

The aims of this work are to provide accurate photometry in multiple near-infrared broadband filters, to determine the power-law index of the extinction-law toward the central parsec of the Galaxy, to provide measurements of the absolute extinction toward the Galactic center, and finally to measure the spatial variability of extinction on arcsecond scales.We use adaptive optics observations of the central parsec of the Milky Way. Absolute values for the extinction in the H, Ks, and L'-bands as well as of the power-law indices of the H to Ks and Ks to L' extinction-laws are measured based on the well-known properties of red clump stars. Extinction maps are derived based on H-Ks and Ks-L' colors. We present Ks-band photometry for ~7700 stars (H and L' photometry for a subset). From a number of recently published values we compute a mean distance of the Galactic center of R_0=8.03+-0.15 kpc, which has an uncertainty of just 2%. Based on this R_0 and on the RC method, we derive absolute mean extinction values toward the central parsec of the Galaxy of A_H=4.48+-0.13 mag, A_Ks=2.54+-0.12$ mag, and A_L'=1.27+-0.18 mag. We estimate values of the power-law indices of the extinction-law of alpha_{H-Ks}=2.21+-0.24 and alpha_{Ks-L'}=1.34+-0.29. A Ks-band extinction map for the Galactic center is computed based on this extinction law and on stellar H-Ks colors. Mean extinction values in a circular region with 0.5" radius centered on Sagittarius A* are A_{H, SgrA*}=4.35+-0.12, A_{Ks, SgrA*}=2.46+-0.03, and A_{L', SgrA*}=1.23+-0.08.

The massive post-Main Sequence star W243 in the galactic starburst cluster Westerlund 1 has undergone a spectral transformation from a B2Ia supergiant devoid of emission features in 1981 to an A-type supergiant with a rich emission-line spectrum by 2002/03. We used VLT/UVES and VLT/FLAMES to obtain high-resolution spectra on six epochs in 2003/04 (UVES) and ten epochs in 2008/09 (FLAMES). These spectra are used alongside other VLT/FLAMES and NTT/EMMI spectra to follow the evolution of W243 from 2002 to 2009. W243 displays a complex, time-varying spectrum with emission lines of Hydrogen, Helium and Lyman-pumped metals, forbidden lines of Nitrogen and Iron, and a large number of absorption lines from neutral and singly-ionized metals. Many lines are complex emission/absorption blends. The LBV has a current temperature of ~8500K (spectral type A3Ia+), and displays signs of photospheric pulsations and weak episodic mass loss. Nitrogen is highly overabundant, with Carbon and Oxygen depleted, indicative of surface CNO-processed material and considerable previous mass-loss, although current time-averaged mass-loss rates are low. The emission-line spectrum forms at large radii, when material lost by the LBV in a previous mass-loss event is ionized by an unseen hot companion. Monitoring of the near-infrared spectrum suggests that the star has not changed significantly since it finished evolving to the cool state, close to the Humphreys-Davidson limit, in early 2003. [ABRIDGED]

We investigated the nature of the massive star [OMN2000] LS1 and used these results to constrain the history of star formation within the host complex W51. A combination of near-IR spectroscopy and non-LTE model atmosphere analysis was used to derive the physical properties of [OMN2000] LS1 and a combination of theoretical evolutionary calculations and Monte Carlo simulations to apply limits on the star formation history of W51. The spectrum of [OMN2000] LS1 is consistent with that of a P Cygni supergiant, but with a temperature of 13.2-13.7kK and log(L/L_sun)<5.75, it is significantly cooler, less luminous, and massive than proposed by previous authors. The presence of such a star within W51 shows that star formation has been underway for at least 3Myr, while the formation of massive O stars is still on going. We find no evidence of internally triggered, sequential star formation within W51, and favour the suggestion that star formation has proceeded at multiple indepedent sites within the GMC. [ABRIDGED]

The Tarantula Survey is an ambitious ESO Large Programme that has obtained multi-epoch spectroscopy of over 1,000 massive stars in the 30 Doradus region of the Large Magellanic Cloud. Here we introduce the scientific motivations of the survey and give an overview of the observational sample. Ultimately, quantitative analysis of every star, paying particular attention to the effects of rotational mixing and binarity, will be used to address fundamental questions in both stellar and cluster evolution.

We review recent developments regarding radiation driven mass loss from OB-stars. We first summarize the fundamental theoretical predictions, and then compare these to observational results (including the VLT-FLAMES survey of massive stars). Especially we focus on the mass loss-metallicity dependence and on the so-called bi-stability jump. Subsequently we concentrate on two urgent problems, weak winds and wind clumping, that have been identified from various diagnostics and that challenge our present understanding of radiation driven winds. We discuss the problems of ``measuring'' mass-loss rates from weak winds and the potential of the near infrared, Br-alpha line as a tool to enable a more precise quantification, and comment on physical explanations for mass-loss rates that are much lower than predicted by the standard model. Wind clumping, conventionally interpreted as the consequence of a strong instability inherent to radiative line-driving, has severe implications for the interpretation of observational diagnostics, since derived mass-loss rates are usually overestimated when clumping is present but ignored in the analyses. Simplified techniques to account for clumping indicate overestimates by factors of 2 to 10, or even more. If actually true, these results would have a dramatic impact on the evolution of, and the feedback from, massive stars. We discuss ongoing attempts (including own work) to interprete the corresponding observations in terms of more sophisticated models. By allowing for porosity in density and velocity space, and for a non-void inter-clump medium, such models might require only moderate reductions of mass-loss rates.

- Context: We present the results of a near-infrared photometric and spectroscopic study of the star forming region G61.48+0.09. - Aims: The purpose of this study is to characterize the stellar content of the cluster and to determine its distance, extinction, age and mass. - Methods: The stellar population was studied by using color-magnitude diagrams to select twenty promising cluster members, for which follow up spectroscopy was done. The observed spectra allowed a spectral classification of the stars. - Results: Two stars have emission lines, twelve are G-type stars, and six are late-O or early-B stars. - Conclusions: The cluster's extinction varies from A_{K_S} = 0.9 to A_{K_S} = 2.6, (or A_{V}~8 to A_{V}~23). G61.48+0.09 is a star forming region located at 2.5+/-0.4 Kpc. The cluster is younger than 10 Myr and has a minimum stellar mass of 1500+/-500 Solar masses. However, the actual total mass of the cluster remains undetermined, as we cannot see its whole stellar content.

The location of the Scutum Red-Supergiant (RSG) clusters at the end of the Galactic Bar makes them an excellent probe of the Galaxy's secular evolution; while the clusters themselves are ideal testbeds in which to study the predictions of stellar evolutionary theory. To this end, we present a study of the RSGs' surface abundances using a combination of high-resolution H-band spectroscopy and spectral synthesis analysis. We provide abundance measurements for elements C, O, Si, Mg, Ti, and Fe. We find that the surface abundances of the stars studied are consistent with CNO burning and deep, rotationally enhanced mixing. The average a/Fe ratios of the clusters are solar, consistent with a thin-disk population. However, we find significantly sub-solar Fe/H ratios for each cluster, a result which strongly contradicts a simple extrapolation of the Galactic metallicity gradient to lower Galacto-centric distances. We suggest that a simple one-dimensional parameterization of the Galaxy's abundance patterns is insufficient at low Galactocentric distances, as large azimuthal variations may be present. Indeed, we show that the abundances of O, Si and Mg are consistent with independent measurements of objects in similar locations in the Galaxy. In combining our results with other data in the literature, we present evidence for large-scale (~kpc) azimuthal variations in abundances at Galacto-centric distances of 3-5kpc. While we cannot rule-out that this observed behaviour is due to systematic offsets between different measurement techniques, we do find evidence for similar behaviour in a study of the barred-spiral galaxy NGC4736 which uses homogeneous methodology. We suggest that these azimuthal abundance variations could result from the intense but patchy star formation driven by the potential of the central bar.

The Galactic Centre (GC) has experienced a high degree of recent star-forming activity, as evidenced by the large number of massive stars currently residing there. The relative abundances of chemical elements in the GC may provide insights into the origins of this activity. Here, we present high-resolution $H$-band spectra of two Red Supergiants in the GC (IRS~7 and VR~5-7), and in combination with spectral synthesis we derive abundances for Fe and C, as well as other $\alpha$-elements Ca, Si, Mg Ti and O. We find that the C-depletion in VR~5-7 is consistent with the predictions of evolutionary models of RSGs, while the heavy depletion of C and O in IRS~7's atmosphere is indicative of deep mixing, possibly due to fast initial rotation and/or enhanced mass-loss. Our results indicate that the {\it current} surface Fe/H content of each star is slightly above Solar. However, comparisons to evolutionary models indicate that the {\it initial} Fe/H ratio was likely closer to Solar, and has been driven higher by H-depletion at the stars' surface. Overall, we find $\alpha$/Fe ratios for both stars which are consistent with the thin Galactic disk. These results are consistent with other chemical studies of the GC, given the precision to which abundances can currently be determined. We argue that the GC abundances are consistent with a scenario in which the recent star-forming activity in the GC was fuelled by either material travelling down the Bar from the inner disk, or from the winds of stars in the inner Bulge -- with no need to invoke top-heavy stellar Initial Mass Functions to explain anomalous abundance ratios.

We report on a successful, simultaneous observation and modeling of the sub-millimeter to near-infrared flare emission of the Sgr A* counterpart associated with the super-massive black hole at the Galactic center. Our modeling is based on simultaneous observations that have been carried out on 03 June, 2008 using the NACO adaptive optics (AO) instrument at the ESO VLT and the LABOCA bolometer at the APEX telescope. Inspection and modeling of the light curves show that the sub-mm follows the NIR emission with a delay of 1.5+/-0.5 hours. We explain the flare emission delay by an adiabatic expansion of the source components.

Mass loss is a key process in the evolution of massive stars, and must be understood quantitatively to be successfully included in broader astrophysical applications. In this review, we discuss various aspects of radiation driven mass loss, both from the theoretical and the observational side. We focus on winds from OB-stars, with some excursions to the Luminous Blue Variables, Wolf- Rayet stars, A-supergiants and Central Stars of Planetary Nebulae. After reca- pitulating the 1-D, stationary standard model of line-driven wind, extensions accounting for rotation and magnetic fields are discussed. The relevance of the so-called bi-stability jump is outlined. We summarize diagnostical methods to infer wind properties from observations, and compare the results with theore- tical predictions, featuring the massloss-metallicity dependence. Subsequently, we concentrate on two urgent problems which challenge our present understanding of radiation driven winds: weak winds and wind- clumping. We discuss problems of measuring mass-loss rates from weak winds and the potential of NIR- spectroscopy. Wind-clumping has severe implications for the interpretation of observational diagnostics, as derived mass-loss rates can be overestimated by factors of 2 to 10 if clumping is ignored, and we describe ongoing attempts to allow for more uniform results. We point out that independent arguments from stellar evolution favor a moderate reduction of present- day mass-loss rates. We also consider larger scale wind structure, interpreted in terms of co-rotating interacting regions, and complete this review with a discussion of recent progress on the X-ray line emission from massive stars, highlighting as to how far the analysis of such X-ray line emission can give further clues regarding an adequate description of wind clumping. (Abridged abstract)

At the center of the Milky Way, with a distance of ~8 kpc, the compact source Sagittarius A* (SgrA*) can be associated with a super massive black hole of ~4x10^6 solar masses. SgrA* shows strong variability from the radio to the X-ray wavelength domains. Here we report on simultaneous NIR/sub-millimeter/X-ray observations from May 2007 that involved the NACO adaptive optics (AO) instrument at the European Southern Observatory's Very Large Telescope, the Australian Telescope Compact Array (ATCA), the US mm-array CARMA, the IRAM 30m mm-telescope, and other telescopes. We concentrate on the time series of mm/sub-mm data from CARMA, ATCA, and the MAMBO bolometer at the IRAM 30m telescope.

We report on recent near-infrared (NIR) and X-ray observations of Sagittarius A* (Sgr A*), the electromagnetic manifestation of the ~4x10^6 solar masses super-massive black hole (SMBH) at the Galactic Center. The goal of these coordinated multi-wavelength observations is to investigate the variable emission from Sgr A* in order to obtain a better understanding of the underlying physical processes in the accretion flow/outflow. The observations have been carried out using the NACO adaptive optics (AO) instrument at the European Southern Observatory's Very Large Telescope (July 2005, May 2007) and the ACIS-I instrument aboard the Chandra X-ray Observatory (July 2005). We report on a polarized NIR flare synchronous to a 8x1033 erg/s X-ray flare in July 2005, and a further flare in May 2007 that shows the highest sub-flare to flare contrast observed until now. The observations can be interpreted in the framework of a model involving a temporary disk with a short jet. In the disk component flux density variations can be explained due to hot spots on relativistic orbits around the central SMBH. The variations of the sub-structures of the May 2007 flare are interpreted as a variation of the hot spot structure due to differential rotation within the disk.

We present a measurement of metallicity in the Galactic center Quintuplet Cluster made using quantitative spectral analysis of two Luminous Blue Variables (LBVs). The analysis employs line-blanketed NLTE wind/atmosphere models fit to high-resolution near-infrared spectra containing lines of H, HeI, SiII, MgII, and FeII. We are able to break the H/He ratio vs. mass-loss rate degeneracy found in other LBVs and to obtain robust estimates of the He content of both objects. Our results indicate solar iron abundance and roughly twice solar abundance in the alpha-elements. These results are discussed within the framework of recent measurements of oxygen and carbon composition in the nearby Arches Cluster and iron abundances in red giants and supergiants within the central 30 pc of the Galaxy. The relatively large enrichment of alpha-elements with respect to iron is consistent with a history of more nucleosynthesis in high mass stars than the Galactic disk.

We present the 3.5m SPICA space telescope, a proposed Japanese-led JAXA-ESA mission scheduled for launch around 2017. The actively cooled (<5 K), single aperture telescope and monolithic mirror will operate from ~3.5 to ~210 um and will provide superb sensitivity in the mid- and far-IR spectral domain (better than JWST at lambda > 18 um). SPICA is one of the few space missions selected to go to the next stage of ESA's Cosmic Vision 2015-2025 selection process. In this White Paper we present the main specifications of the three instruments currently baselined for SPICA: a mid-infrared (MIR) coronagraph (~3.5 to ~27 um) with photometric and spectral capabilities (R~200), a MIR wide-field camera and high resolution spectrometer (R~30,000), and a far-infrared (FIR ~30 to ~210 um) imaging spectrometer - SAFARI - led by a European consortium. We discuss their capabilities in the context of MIR direct observations of exo-planets (EPs) and multiband photometry/high resolution spectroscopy observations of transiting exo-planets. We conclude that SPICA will be able to characterize the atmospheres of transiting exo-planets down to the super-Earth size previously detected by ground- or space-based observatories. It will also directly detect and characterize Jupiter/Neptune-size planets orbiting at larger separation from their parent star (>5-10 AU), by performing quantitative atmospheric spectroscopy and studying proto-planetary and debris disks. In addition, SPICA will be a scientific and technological precursor for future, more ambitious, IR space missions for exo-planet direct detection as it will, for example, quantify the prevalence exo-zodiacal clouds in planetary systems and test coronographic techniques, cryogenic systems and lightweight, high quality telescopes. (abridged)

We present results from K band slit scan observations of a ~20''x20'' region of the Galactic center (GC) in two separate epochs more than five years apart. The high resolution (R>=14,000) observations allow the most accurate radial velocity and acceleration measurements of the stars in the central parsec of the Galaxy. Detected stars can be divided into three groups based on the CO absorption band heads at ~2.2935 microns and the He I lines at ~2.0581 microns and ~2.112, 2.113 microns: cool, narrow-line hot and broad-line hot. The radial velocities of the cool, late-type stars have approximately a symmetrical distribution with its center at ~-7.8(+/-10.3) km/s and a standard deviation ~113.7(+/-10.3) km/s. Although our statistics are dominated by the brightest stars, we estimate a central black hole mass of 3.9(+/-1.1) million solar masses, consistent with current estimates from complete orbits of individual stars. Our surface density profile and the velocity dispersion of the late type stars support the existence of a low density region at the Galactic center suggested by earlier observations. Many hot, early-type stars show radial velocity changes higher than maximum values allowed by pure circular orbital motions around a central massive object, suggesting that the motions of these stars greatly deviate from circular orbital motions around the Galactic center. The correlation between the radial velocities of the early type He I stars and their declination offsets from Sagittarius A* suggests that a systematic rotation is present for the early-type population. No figure rotation around the Galactic center for the late type stars is supported by the new observations.

Stellar dynamics indicate the presence of a super massive 3-4x10^6 Msun solm black hole at the Galactic Center. It is associated with the variable radio, near-infrared, and X-ray counterpart Sagittarius A* (SgrA*). The goal is the investigation and understanding of the physical processes responsible for the variable emission from SgrA*. The observations have been carried out using the NACO adaptive optics (AO) instrument at the European Southern Observatory's Very Large Telescope (July 2005, May 2007) and the ACIS-I instrument aboard the Chandra X-ray Observatory (July 2005). We find that for the July 2005 flare the variable and polarized NIR emission of SgrA* occurred synchronous with a moderately bright flare event in the X-ray domain with an excess 2 - 8 keV luminosity of about 8x10^33erg/s. We find no time lag between the flare events in the two wavelength bands with a lower limit of less than 10 minutes. The May 2007 flare shows the highest sub-flare to flare contrast observed until now. It provides evidence for a variation in the profile of consecutive sub-flares. We confirm that highly variable and NIR polarized flare emission is non-thermal and that there exists a class of synchronous NIR/X-ray flares. We find that the flaring state can be explained via the synchrotron self-Compton (SSC) process involving up-scattered sub-millimeter photons from a compact source component. The observations can be interpreted in a model involving a temporary disk with a short jet. In the disk component the flux density variations can be explained due to spots on relativistic orbits around the central super massive black hole (SMBH). The profile variations for the May 2007 flare are interpreted as a variation of the spot structure due to differential rotation within the disk.

We present new high-resolution near-IR spectroscopy and OH maser observations to investigate the population of cool luminous stars of the young massive Galactic cluster RSGC1. Using the 2.293\micron CO-bandhead feature, we make high-precision radial velocity measurements of 16 of the 17 candidate Red Supergiants (RSGs) identified by Figer et al. We show that F16 and F17 are foreground stars, while we confirm that the rest are indeed physically-associated RSGs. We determine that Star F15, also associated with the cluster, is a Yellow Hypergiant based on its luminosity and spectroscopic similarity to $\rho$ Cas. Using the cluster's radial velocity, we have derived the kinematic distance to the cluster and revisited the stars' temperatures and luminosities. We find a larger spread of luminosities than in the discovery paper, consistent with a cluster age 30% older than previously thought (12$\pm$2Myr), and a total initial mass of $(3\pm1) \times 10^{4}$\msun. The spatial coincidence of the OH maser with F13, combined with similar radial velocities, is compelling evidence that the two are related. Combining our results with recent SiO and H$_2$O maser observations, we find that those stars with maser emission are the most luminous in the cluster. From this we suggest that the maser-active phase is associated with the end of the RSG stage, when the luminosity-mass ratios are at their highest.

We report on the unprecedented Red Supergiant (RSG) population of a massive young cluster, located at the base of the Scutum-Crux Galactic arm. We identify candidate cluster RSGs based on {\it 2MASS} photometry and medium resolution spectroscopy. With follow-up high-resolution spectroscopy, we use CO-bandhead equivalent width and high-precision radial velocity measurements to identify a core grouping of 26 physically-associated RSGs -- the largest such cluster known to-date. Using the stars' velocity dispersion, and their inferred luminosities in conjuction with evolutionary models, we argue that the cluster has an initial mass of $\sim$40,000\msun, and is therefore among the most massive in the galaxy. Further, the cluster is only a few hundred parsecs away from the cluster of 14 RSGs recently reported by Figer et al (2006). These two RSG clusters represent 20% of all known RSGs in the Galaxy, and now offer the unique opportunity to study the pre-supernova evolution of massive stars, and the Blue- to Red-Supergiant ratio at uniform metallicity. We use GLIMPSE, MIPSGAL and MAGPIS survey data to identify several objects in the field of the larger cluster which seem to be indicative of recent region-wide starburst activity at the point where the Scutum-Crux arm intercepts the Galactic bulge. Future abundance studies of these clusters will therefore permit the study of the chemical evolution and metallicity gradient of the Galaxy in the region where the disk meets the bulge.

[Abridged] We present a comprehensive study of the metallicity dependence of the mass-loss rates in stationary stellar winds of hot massive stars. Assuming a power-law dependence of mass loss on metallicity, Mdot \propto Z^{m}, and adopting a theoretical relation between the terminal velocity and metallicity, v_inf \propto Z^{0.13} (Leitherer et al.), we find m = 0.83 +/- 0.16 for non-clumped outflows from an analysis of the wind momentum luminosity relation (WLR) for stars more luminous than 10^{5.2} Lsun. Within the errors, this result agrees with the prediction of m = 0.69 +/- 0.10 from Vink et al. For the high luminosity stars we find the mass loss rates to be greater than the predictions, implying wind clumping factors in their line-forming regions of ~4. For lower luminosity stars, the winds are so weak that their strengths cannot be reliably derived from optical lines, and one must rely on analysis of UV lines. In the low-luminosity domain the Galactic WLR is found to be much steeper than expected from theory, leading to a discrepancy between UV mass-loss rates and the predictions by a factor 100 at luminosities of L ~ 10^{4.75} Lsun, the origin of which is unknown. We emphasize that even if the current mass-loss rates of hot luminous stars are overestimated as a result of wind clumping, the degree of clumping would likely to be independent of metallicity, so the scalings derived in this study are expected to remain correct.

[Abridged] We have studied the optical spectra of 28 O- and early B-type stars in the Large Magellanic Cloud, 22 of which are associated with the young star-forming region N11. Stellar parameters are determined using an automated fitting method, combining the stellar atmosphere code FASTWIND with the genetic-algorithm optimisation routine PIKAIA. Results for stars in the LH9 and LH10 associations of N11 are consistent with a sequential star formation scenario, in which activity in LH9 triggered the formation of LH10. Our sample contains four stars of spectral type O2, of which the hottest is found to be ~49-54 kK (cf. ~45-46 kK for O3 stars). The masses of helium-enriched dwarfs and giants are systematically lower than those implied by non-rotating evolutionary tracks. We interpret this as evidence for efficient rotationally-enhanced mixing, leading to the surfacing of primary helium and to an increase of the stellar luminosity. This result is consistent with findings for SMC stars by Mokiem et al. For bright giants and supergiants no such mass-discrepancy is found, implying that these stars follow tracks of modestly (or non-)rotating objects. Stellar mass-loss properties were found to be intermediate to those found in massive stars in the Galaxy and the SMC, and comparisons with theoretical predictions at LMC metallicity yielded good agreement over the luminosity range of our targets, i.e. 5.0 < log L/L(sun) < 6.1.

We have analyzed H and K_s-band images of the Arches cluster obtained using the NIRC2 instrument on Keck with the laser guide star adaptive optics (LGS AO) system. With the help of the LGS AO system, we were able to obtain the deepest ever photometry for this cluster and its neighborhood, and derive the background-subtracted present-day mass function (PDMF) down to 1.3 Msun for the 5 arcsec-9 arcsec annulus of the cluster. We find that the previously reported turnover at 6 Msun is simply due to a local bump in the mass function (MF), and that the MF continues to increase down to our 50 % completeness limit (1.3 Msun) with a power-law exponent of Gamma = -0.91 for the mass range of 1.3 < M/Msun < 50. Our numerical calculations for the evolution of the Arches cluster show that the Gamma values for our annulus increase by 0.1-0.2 during the lifetime of the cluster, and thus suggest that the Arches cluster initially had Gamma of -1.0 ~ -1.1, which is only slightly shallower than the Salpeter value.

Model atmosphere fits to high-resolution optical spectra of Wray 977 confirm the B hypergiant classification of the massive companion to the X-ray pulsar GX301-2. The models give a radius of 62 Rsun, an effective temperature of 18,100 K and a luminosity of 5 x 10^^5 Lsun. The deduced mass-loss rate and terminal velocity of the stellar wind are 10^^-5 Msun/yr and 305 km/s, respectively. The interstellar Na I D absorption indicates that Wray 977 is located behind the first intersection with the Sagittarius-Carina spiral arm (1-2.5 kpc) and probably belongs to the stellar population of the Norma spiral arm at a distance of 3-4 kpc. The spectra obtained with UVES on the Very Large Telescope (VLT) cover a full orbit of the system, including periastron passage, from which we derive the radial-velocity curve of the B hypergiant. The measured radial-velocity amplitude is 10 +/-3 km/s yielding a mass ratio q = 0.046 +/- 0.014. The absence of an X-ray eclipse results in a lower limit to the mass of Wray 977 of 39 Msun. An upper limit of 53 Msun is derived for the mass of Wray 977 adopting a maximum neutron star mass of 2.5 Msun. The ``spectroscopic'' mass of Wray 977 is 43 +/- 10 Msun, consistent with the range in mass derived from the binarity constraints. The mass of the neutron star is 1.85 +/- 0.6 Msun. Time series of spectral lines formed in the dense stellar wind indicate the presence of a gas stream trailing the neutron star in its orbit. The long-term behaviour of the H alpha equivalent width exhibits strong variations in wind strength; the sampling of the data is insufficient to conclude whether a relation exists between wind mass-loss rate and pulsar spin period.

Using adaptive optics at the Gemini North telescope we have obtained a K-band spectrum of the star near the center of the luminous Galactic center bowshock IRS8, as well as a spectrum of the bowshock itself. The stellar spectrum contains emission and absorption lines characteristic of an O5-O6 giant or supergiant. The wind from such a star is fully capable of producing the observed bowshock. However, both the early spectral type and the apparently young age of the star, if it is single, mark it as unique among hot stars within one parsec of the center.

We have studied the optical spectra of a sample of 31 O- and early B-type stars in the Small Magellanic Cloud, 21 of which are associated with the young massive cluster NGC 346. Stellar parameters are determined using an automated fitting method. Comparison with predictions of stellar evolution that account for stellar rotation does not result in a unique age, though most stars are best represented by an age of 1-3 Myr. The present day v_sini distribution of the 21 dwarf stars in our sample is consistent with an underlying rotational velocity (v_r) distribution that can be characterised by a mean velocity of about 160-190 km/s and an effective half width of 100-150 km/s. The v_r distribution must include a small percentage of slowly rotating stars. If predictions of the time evolution of the equatorial velocity for massive stars within the environment of the SMC are correct, the young age of the cluster implies that this underlying distribution is representative for the initial rotational velocity distribution. The location in the Hertzsprung-Russell diagram of the stars showing helium enrichment is in qualitative agreement with evolutionary tracks accounting for rotation, but not for those ignoring v_r. The mass loss rates of the SMC objects having luminosities of log L/L_sun > 5.4 are in excellent agreement with predictions. However, for lower luminosity stars the winds are too weak to determine M_dot accurately from the optical spectrum. Two of three spectroscopically classified Vz stars from our sample are located close to the theoretical zero age main sequence, as expected.

Recent calculations and analyses of O star spectra have revealed discrepancies between theory and observations, and between different theoretical calculations, for the strength of optical HeI singlet transitions.We investigate the source of these discrepancies. Using a non-LTE radiative transfer code we have undertaken detailed test calculations for a range of O star properties. Our principal test model has parameters similar to those of the O9V star, 10 Lac. We show that the discrepancies arise from uncertainties in the radiation field in the HeI resonance transition near 584Angs. The radiation field at 584Angs. is influenced by model assumptions, such as the treatment of line-blanketing and the adopted turbulent velocity, and by the FeIV atomic data. It isshown that two FeIV transitions near 584Angs can have a substantial influence on the strength of the HeI singlet transitions. Because of the difficulty of modeling the HeI singlet lines, particularly in stars with solar metalicity, the HeI triplet lines should be preferred in spectral analyses. These lines are much less sensitive to model assumptions.

We report the discovery of an extraordinarily massive young cluster of stars in the Galaxy, having an inferred total initial cluster mass comparable to the most massive young clusters in the Galaxy. Using {\it IRMOS}, {\it 2MASS}, and {\it Spitzer} observations, we conclude that there are 14 red supergiants in the cluster, compared with five, in what was previously thought to be the richest Galactic cluster of such stars. We infer spectral types from near-infrared spectra that reveal deep CO bandhead absorption that can only be fit by red supergiants. We identify a gap of $\Delta${\it K}$_s$$\sim$4 magnitudes between the stars and the bulk of the other stars in the region that can only be fit by models if the brightest stars in the cluster are red supergiants. We estimate a distance of 5.8~\kpc to the cluster by associating an OH maser with the envelope of one of the stars. We also identify a ``yellow'' supergiant of G6~I type in the cluster. Assuming a Salpeter IMF, we infer an initial cluster mass of 20,000 to 40,000~\Msun for cluster ages of 7-12~\Myr. Continuing with these assumptions, we find 80% of the intial mass and 99% of the number of stars remain at the present time. We associate the cluster with an x-ray source (detected by {\it ASCA} and {\it Einstein}), a recently discovered very high energy $\gamma$-ray source (detected by {\it INTEGRAL} and {\it HESS}), and several non-thermal radio sources, finding that these objects are likely related to recent supernovae in the cluster. In particular, we claim that the cluster has produced at least one recent supernova remnant with properties similar to the Crab nebula. It is not unlikely to find such a source in this cluster, given our estimated supernova rate of one per 40,000 to 80,000~{\it yr}.

We present the results of a spectroscopic analysis of the Trapezium cluster stars inside the Orion Nebula. The rotational velocities have been obtained using Fourier analysis method, finding agreement with values derived from the usual method, based on linewidth measurements. The rotational velocity derived for theta 1 Ori C through this method is consistent with the variability of some of its spectral features that have a period of 15.42 days. By means of the fit of H, HeI and HeII observed profiles with FASTWIND synthetic profiles, stellar parameters and wind characteristics have been derived. This methodology let us estimate the errors associated with these parameters. It is found that macroturbulence effects have to be included for a good fit to the HeI-II lines in the spectrum of theta 1 Ori C. By means of a very accurate study, oxygen abundances have been derived for the three B0.5V stars theta 1 Ori A, D and theta 2 Ori B. Final abundances are consistent with the nebular gas-phase results presented in Esteban et al. (2004) and are lower than those given by Cunha and Lambert (1994). Our results suggest a lower dust depletion factor of oxygen than previous estimations for the Orion nebula.

This paper presents L-band ISAAC and NAOS/CONICA (VLT) spectroscopic observations of the IRS~3-IRS~13 Galactic Center region. The ISAAC data allowed us to build the first spectroscopic data cube of the region in the L-band domain corrected for the foreground extinction. Maps of the water ice and hydrocarbon absorption line strength were then derived. These maps are important diagnostics of the interstellar and circumstellar medium (resp. ISM and CSM). They support our previous results that the absorption features are most probably occuring in the local Galactic center medium and can be associated with the individual sources. Moreover, turbulence seems to affect the studied region of the minispiral which appears like a mixture of a dense and diffuse medium. The Br alpha and Pf gamma emission line maps allowed us to find three sources with broad lines corresponding to a FWHM deconvolved line width of about 1100 km/s and moving towards us with a radial velocity of about -300km/s. These sources are most probably new Wolf-Rayet type stars located in projection to the north and west of IRS 3. Their derived radial velocities and proper motions show that only two of them might belong to the two rotating disks of young stars reported by Genzel et al. (2003) and Levin & Beloborodov (2003)...

We present the first automated fitting method for the quantitative spectroscopy of O- and early B-type stars with stellar winds. The method combines the non-LTE stellar atmosphere code FASTWIND from Puls et al. (2005) with the genetic algorithm based optimizing routine PIKAIA from Charbonneau (1995), allowing for a homogeneous analysis of upcoming large samples of early-type stars (e.g. Evans et al. 2005). In this first implementation we use continuum normalized optical hydrogen and helium lines to determine photospheric and wind parameters. We have assigned weights to these lines accounting for line blends with species not taken into account, lacking physics, and/or possible or potential problems in the model atmosphere code. We find the method to be robust, fast, and accurate. Using our method we analysed seven O-type stars in the young cluster Cyg OB2 and five other Galactic stars with high rotational velocities and/or low mass loss rates (including 10lac, zeta Oph, and tau Sco) that have been studied in detail with a previous version of FASTWIND. The fits are found to have a quality that is comparable or even better than produced by the classical ``by eye'' method. We define errorbars on the model parameters based on the maximum variations of these parameters in the models that cluster around the global optimum. Using this concept, for the investigated dataset we are able to recover mass-loss rates down to ~6d-8 Msun/yr to within an error of a factor of two. Comparison of our derived spectroscopic masses with those derived from stellar evolutionary models are in very good agreement, i.e. based on the limited sample that we have studied we do not find indications for a mass discrepancy.

We introduce a new survey of massive stars in the Galaxy and the Magellanic Clouds using the Fibre Large Array Multi-Element Spectrograph (FLAMES) instrument at the Very Large Telescope (VLT). Here we present observations of 269 Galactic stars with the FLAMES-Giraffe Spectrograph (R ~ 25,000), in fields centered on the open clusters NGC 3293, NGC 4755 and NGC 6611. These data are supplemented by a further 50 targets observed with the Fibre-Fed Extended Range Optical Spectrograph (FEROS, R = 48,000). Following a description of our scientific motivations and target selection criteria, the data reduction methods are described; of critical importance the FLAMES reduction pipeline is found to yield spectra that are in excellent agreement with less automated methods. Spectral classifications and radial velocity measurements are presented for each star, with particular attention paid to morphological peculiarities and evidence of binarity. These observations represent a significant increase in the known spectral content of NGC 3293 and NGC 4755, and will serve as standards against which our subsequent FLAMES observations in the Magellanic Clouds will be compared.

We report the discovery of additional hot and massive stars in the cluster surrounding the soft gamma repeater SGR 1806-20, based upon UKIRT and Keck near-infrared spectroscopy. Of the newly identified stars, three are Wolf-Rayet stars of types WC8, WN6, and WN7, and a fourth star is an OB supergiant. These three stars, along with four previously discovered, imply a cluster age of \~3.0-4.5 Myr, based on the presence of WC stars and the absence of red supergiants. Assuming coevality, this age suggests that the progenitor of SGR 1806-20 had an initial mass greater than ~50 Msun. This is consistent with the suggestion that SGRs are post-supernovae end states of massive progenitors, and may suggest that only massive stars evolve into magnetars that produce SGRs. It also suggests that very massive stars can evolve into neutron stars, not just black holes, as recently predicted by theory. The cluster age also provides constraints on the very high mass object, LBV 1806-20.

We present a quantitative spectral analysis of five very massive stars in the Arches cluster, located near the Galactic center, to determine stellar parameters, stellar wind properties and, most importantly, metallicity content. The analysis uses a new technique, presented here for the first time, and uses line-blanketed NLTE wind/atmosphere models fit to high-resolution near-infrared spectra of late-type nitrogen-rich Wolf-Rayet stars and OfI+ stars in the cluster. It relies on the fact that massive stars reach a maximum nitrogen abundance that is related to initial metallicity when they are in the WNL phase. We determine the present-day nitrogen abundance of the WNL stars in the Arches cluster to be 1.6% (mass fraction) and constrain the stellar metallicity in the cluster to be solar. This result is invariant to assumptions about the mass-luminosity relationship, the mass-loss rates, and rotation speeds. In addition, from this analysis, we find the age of the Arches cluster to be 2-2.5Myr, assuming coeval formation.

Despite much theoretical and observational progress, there is no known firm upper limit to the masses of stars. Our understanding of the interplay between the immense radiation pressure produced by massive stars in formation and the opacity of infalling material is subject to theoretical uncertainties, and many observational claims of ``the most massive star'' have failed the singularity test. LBV 1806-20 is a particularly luminous object, L~10^6 Lsun, for which some have claimed very high mass estimates (M_initial>200 Msun), based, in part, on its similarity to the Pistol Star. We present high-resolution near-infrared spectroscopy of LBV 1806-20, showing that it is possibly a binary system with components separated in velocity by ~70 kms. If correct, then this system is not the most massive star known, yet it is a massive binary system. We argue that a binary, or merged, system is more consistent with the ages of nearby stars in the LBV 1806-20 cluster. In addition, we find that the velocity of V_LSR=36 kms is consistent with a distance of 11.8 kpc, a luminosity of 10^6.3 Lsun, and a system mass of ~130 Msun.

We have analyzed six OB supergiants and one giant covering spectral types from O3 to B1 in the Galactic OB association Cyg OB2 by means of an updated version of FASTWIND (SAntolaya-Rey, Puls, Herrero, 1997, A&A 323, 348) that includes an approximate treatment of metal line blocking and blanketing. This large coverage in spectral type allows us to derive a new temperature scale for Galactic O supergiants that is lower than the one obtained by using pure H--He models, either plane-parallel and hydrostatic or spherical with mass-loss. The lower temperatures are thus a combined effect of line blanketing and the large mass-loss rates. In some cases, the newly derived effective temperature is reduced by up to 8000 K. Changes are larger for earlier stars with large mass--loss rates. As a consequence, luminosities are modified as well, which results in a lower number of emerging ionizing photons and reduces the mass discrepancy. Although there are still significant differences between spectroscopic and evolutionary masses, we do not find any obvious systematic pattern of those differences. We derive mass--loss rates and the corresponding wind momentum--luminosity relation for the analyzed stars. Although consistent with previous results by Puls et al. (1996, A&A 305, 171) for Galactic stars, our relation is better defined due to a reduction of errors related to stellar distances and points to a possible separation between extreme Of stars (Of+, Of*) and stars with more moderate morphologies. However this finding is only tentative, as the statistics are still scarce.

We have discovered an Ofpe/WN9 (WN11 following Smith et al.) star in the Sculptor spiral galaxy NGC 300, the first object of this class found outside the Local Group, during a recent spectroscopic survey of blue supergiant stars obtained at the ESO VLT. The light curve over a five-month period in late 1999 displays a variability at the 0.1 mag level. The intermediate resolution spectra (3800-7200 A) show a very close resemblance to the Galactic LBV AG Car during minimum. We have performed a detailed non-LTE analysis of the stellar spectrum, and have derived a chemical abundance pattern which includes H, He, C, N, O, Al, Si and Fe, in addition to the stellar and wind parameters. The derived stellar properties and the He and N surface enrichments are consistent with those of other Local Group WN11 stars in the literature, suggesting a similar quiescent or post-LBV evolutionary status.

We present and use new spectra and narrow-band images, along with previously published broad-band images, of stars in the Arches cluster to extract photometry, astrometry, equivalent width, and velocity information. The data are interpreted with a wind/atmosphere code to determine stellar temperatures, luminosities, mass-loss rates, and abundances. We have doubled the number of known emission-line stars, and we have also made the first spectroscopic identification of the main sequence for any population in the Galactic Center. We conclude that the most massive stars are bona-fide Wolf-Rayet (WR) stars and are some of the most massive stars known, having M_{initial} > 100 Msun, and prodigious winds, Mdot > 10^{-5} Msun yr^{-1}, that are enriched with helium and nitrogen; with these identifications, the Arches cluster contains about 5% of all known WR stars in the Galaxy. We find an upper limit to the velocity dispersion of 22 kms^{-1}, implying an upper limit to the cluster mass of 7(10^4) Msun within a radius of 0.23 pc; we also estimate the bulk heliocentric velocity of the cluster to be v_{cluster,odot} approximately +95 kms^{-1}.

In an attempt to determine the nature of the enigmatic cocoon stars in the Quintuplet Cluster, we have obtained mid-infrared imaging and spectrophotometry of the cluster, using the CAM and SWS instruments on ISO, using SpectroCam-10 on the Palomar 5m telescope, and NICMOS on HST. The spectra show smooth continua with various dust and ice absorption features. These features are all consistent with an interstellar origin, and there is no clear evidence for any circumstellar contribution to these features. We find no spectral line or feature that could elucidate the nature of these sources. Detailed modeling of the silicate absorption features shows that they are best reproduced by the mu Cep profile, which is typical of the interstellar medium, with tau(sil) \sim 2.9. The high spatial resolution mid-IR images show that three of the five cocoon stars have spatially extended and asymmetric envelopes, with diameters of \sim 20,000 AUs. A reddening law similar to that of Lutz (1999) but with silicate features based on the mu Cep profile and normalized to our value of tau(sil) is used to deredden the observed spectrophotometry. The dereddened energy distributions are characterised by temperatures of 750-925 K, somewhat cooler than determined from near IR data alone. Models of optically thin and geometrically thick dust shells, as used by Williams et al. (1987) for very dusty, late-type WC stars, reproduce the observed SEDs from 4 to 17 mic, and imply shell luminosities of log(L/L(sun)) \sim 4.5-4.9 for the brightest four components. An analysis of the various suggestions proposed to explain the nature of the cocoon stars reveals serious problems with all the hypotheses, and the nature of these sources remains an enigma.

We report optical and near-infrared spectra of the central star of the radio source G79.29+0.46, a candidate luminous blue variable. The spectra contain numerous narrow (FWHM < 100 km/s emission lines of which the low-lying hydrogen lines are the strongest, and resemble spectra of other LBVc's and B[e] supergiants. A few prominent infrared lines are unidentified. The terminal wind speed is determined from H alpha to be 110 km/s. The strength of H alpha implies the presence of a very dense wind. Extended emission from H alpha and [N II] was detected but appears to be associated with the Cygnus X region rather than the radio source. Both diffuse interstellar bands and interstellar absorption lines are present in the optical spectrum of the central star, suggesting that there are both diffuse and molecular cloud components to the extinction and implying a minimum distance of 1 kpc and minimum luminosity of \~10^5 L(sun) for the star. The new spectra and their analysis indicate a low excitation, dense, and slowly expanding wind and support the LBVc classification.

The formation and evolution of young low mass stars are characterized by important processes of mass loss and accretion ocurring in the innermost regions of their placentary circumstellar disks. Because of the large obscuration of these disks at optical and infrared wavelengths in the early protostellar stages (Class 0 Sources), they were previously detected only at radio wavelengths using interferometric techniques. We have detected with the Infrared Space Observatory (ISO) the mid-infrared emission associated with the Class 0 protostar VLA1 in the HH1-2 region located in the Orion nebula. The emission arises in the three wavelength windows at 5.3, 6.6 and 7.5 micras where the absorption due to ices and silicates has a local minimum that exposes the central parts of the youngest protostellar systems to mid-infrared investigations. The mid-infrared emission arises from a central source with 4 AU diameter at an averaged temperature of 700 K, deeply embedded in a dense region with a visual extinction of Av=80-100mag.

H and K band moderate resolution and 4 $\mu$m high resolution spectra have been obtained for FMM#362, a bright star in the Quintuplet Cluster near the Galactic Center. The spectral features in these bands closely match those of the Pistol Star, a luminous blue variable and one of the most luminous stars known. The new spectra and previously-obtained photometry imply a very high luminosity for FMM#362, L $\geq 10^6$ \Lsun, and a temperature of 10,000 - 13,000 K. Based on its luminosity, temperature, photometric variability, and similarities to the Pistol Star, we conclude that FMM#362 is a luminous blue variable.

We present initial results of a spectroscopic study of the Pistol and of the cocoon stars in the Quintuplet Cluster. From ISOCAM CVF 5--17 micron spectroscopy of the field of the Pistol Star, we have discovered a nearly spherical shell of hot dust surrounding this star, a probable LBV. This shell is most prominent at lambda > 12 micron, and its morphology clearly indicates that the shell is stellar ejecta. Emission line images show that most of the ionised material is along the northern border of this shell, and its morphology is very similar to that of the Pistol HII region (Yusef-Zadeh & Morris, 1987, AJ, 94, 1178). We thus confirm that the ionisation comes from very hot stars in the core of the Quintuplet Cluster. An SWS spectrum of the Pistol Nebula indicates a harder ionising radiation than could be provided by the Pistol Star, but which is consistent with ionisation from Wolf-Rayet stars in the Quintuplet Cluster. The CVF 5--17 micron spectra of the cocoon stars in the Quintuplet do not show any emission feature that could help elucidate their nature.