The inner regions of the discs of high-mass young stellar objects (HMYSOs) are still poorly known due to the small angular scales and the high visual extinction involved. We deploy near-infrared (NIR) spectro-interferometry to probe the inner gaseous disc in HMYSOs and investigate the origin and physical characteristics of the CO bandhead emission (2.3-2.4 $\mu$m). We present the first GRAVITY/VLTI observations at high spectral (R=4000) and spatial (mas) resolution of the CO overtone transitions in NGC 2024 IRS2. The continuum emission is resolved in all baselines and is slightly asymmetric, displaying small closure phases ($\leq$8$^{\circ}$). Our best ellipsoid model provides a disc inclination of 34$^{\circ}$$\pm$1$^{\circ}$, a disc major axis position angle of 166$^{\circ}$$\pm$1$^{\circ}$, and a disc diameter of 3.99$\pm$0.09 mas (or 1.69$\pm$0.04 au, at a distance of 423 pc). The small closure phase signals in the continuum are modelled with a skewed rim, originating from a pure inclination effect. For the first time, our observations spatially and spectrally resolve the first four CO bandheads. Changes in visibility, as well as differential and closure phases across the bandheads are detected. Both the size and geometry of the CO-emitting region are determined by fitting a bidimensional Gaussian to the continuum-compensated CO bandhead visibilities. The CO-emitting region has a diameter of 2.74$\pm^{0.08}_{0.07}$ mas (1.16$\pm$0.03 au), and is located in the inner gaseous disc, well within the dusty rim, with inclination and $PA$ matching the dusty disc geometry, which indicates that both dusty and gaseous discs are coplanar. Physical and dynamical gas conditions are inferred by modelling the CO spectrum. Finally, we derive a direct measurement of the stellar mass of $M_*\sim$14.7$^{+2}_{-3.6}$ M$_{\odot}$ by combining our interferometric and spectral modelling results.

Using MUSE AO data, we probe the ISM in the local spiral galaxy NGC 7793 at a spatial resolution of $\sim$ 10 pc. We identify HII regions and compile a catalogue of supernova remnants (SNRs), planetary nebulae (PNe) and Wolf Rayet stars (WR). We compute electron densities and temperatures from the [SII]6716/6731 and [SIII]6312/9069 line ratios. We study the properties of the ionised gas through BPT diagrams combined with gas velocity dispersion. We spectroscopically confirm 2 previously detected WR stars and a SNR and report the discovery of 7 WR, 1 SNR, and 2 PNe. The diffuse ionized gas (DIG) fraction is between $\sim$ 27 and 42% depending on the method used to define the HII region boundaries. In agreement with previous studies, we find that the DIG exhibits enhanced [SII]/H$\alpha$ and [NII]/H$\alpha$ ratios and a median temperature that is $\sim$ 3000 K higher than in HII regions. We also observe an apparent inverse correlation between temperature and H$\alpha$ surface brightness. Overall, the observed [SII]6716/6731 ratio is consistent within 1$\sigma$ with $n_e$ < 30 cm$^{-3}$, with an almost identical distribution for the DIG and HII regions. The velocity dispersion of the ionised gas indicates that the DIG has a higher degree of turbulence than the HII regions. Comparison with photoionisation and shock models reveals that the diffuse component can only partially be explained via shocks and that it is most likely consistent with photons leaking from density bounded HII regions or with radiation from evolved field stars. [abridged]

The process of massive star ($M\geq8~M_\odot$) formation is still poorly understood. Observations of massive young stellar objects (MYSOs) are challenging due to their rarity, short formation timescale, large distances, and high circumstellar extinction. Here, we present the results of a spectroscopic analysis of a population of MYSOs in the Large Magellanic Cloud (LMC). We took advantage of the spectral resolution and wavelength coverage of X-shooter (300-2500 nm), mounted on the European Southern Observatory Very Large Telescope, to detect characteristic spectral features in a dozen MYSO candidates near 30 Doradus, the largest starburst region in the Local Group hosting the most massive stars known. The X-shooter spectra are strongly contaminated by nebular emission. We used a scaling method to subtract the nebular contamination from our objects. We detect H$\alpha,\beta$, [O I] 630.0 nm, Ca II infrared triplet, [Fe II] 1643.5 nm, fluorescent Fe II 1687.8 nm, H$_2$ 2121.8 nm, Br$\gamma$, and CO bandhead emission in the spectra of multiple candidates. This leads to the spectroscopic confirmation of 10 candidates as bona fide MYSOs. We compare our observations with photometric observations from the literature and find all MYSOs to have a strong near-infrared excess. We compute lower limits to the brightness and luminosity of the MYSO candidates, confirming the near-infrared excess and the massive nature of the objects. No clear correlation is seen between the Br$\gamma$ luminosity and metallicity. Combining our sample with other LMC samples results in a combined detection rate of disk features such as fluorescent Fe II and CO bandheads which is consistent with the Galactic rate (40\%). Most of our MYSOs show outflow features.

Context: The identification and characterisation of populations of young massive stars in (giant) HII regions provides important constraints on i) the formation process of massive stars and their early feedback on the environment, and ii) the initial conditions for population synthesis models predicting the evolution of ensembles of stars. Aims: We identify and characterise the stellar populations of the following young giant HII regions: M8, G333.6-0.2, and NGC6357. Methods: We acquired H- and K-band spectra of around 200 stars using The K-band KMOS on the ESO Very Large Telescope. The targets for M8 and NGC6357 were selected from the MYStIX project, which combines X-ray observations with near-infrared and mid-infrared data. For G333.6-0.2, the sample selection is based on the near-infrared colours combined with X-ray data. We introduce an automatic spectral classification method in order to obtain temperatures and luminosities for the observed stars. We analyse the stellar populations using their photometric, astrometric, and spectroscopic properties and compared the position of the stars in the Hertzprung-Russell diagram with stellar evolution models to constrain their ages and mass ranges. Results: We confirm the presence of candidate ionising sources in the three regions and report new ones, including the first spectroscopically identified O stars in G333.6-0.2. In M8 and NGC6357, two populations are identified: i) OB main-sequence stars ($M > 5~\rm{M_{\odot}}$) and ii) pre-main sequence stars ($M\approx0.5-5~\rm{M_{\odot}}$). The ages of the clusters are $\sim$1-3~Myr, $< 3$~Myr, and $\sim$0.5-3~Myr for M8, G333.6-0.2, and NGC6357, respectively. We show that MYStIX selected targets have $>$ 90\% probability of being members of the HII region, whereas a selection based on near infrared (NIR) colours leads to a membership probability of only $\sim$70\%.

We have used the capability of the MUSE instrument to explore the impact of stellar feedback at large scales in Haro 11, a galaxy under extreme starburst condition and one of the first galaxies where Lyman continuum (LyC) has been detected. Using Ha, [OIII] and [OI] emission lines from deep MUSE observations, we have constructed a sequence of velocity-dependent maps of the Ha emission, the state of the ionised gas and a tracer of fast shocks. These allowed us to investigate the ionisation structure of the galaxy in 50 kms^2 bins over a velocity range of -400 to 350 kms. The ionised gas in Haro 11 is assembled by a rich arrangement of structures, such as superbubbles, filaments, arcs and galactic ionised channels, whose appearances change drastically with velocity. The central star forming knots and the star forming dusty arm are the main engines that power the strong mechanical feedback in this galaxy, although with different impact on the ionisation structure. Haro 11 appears to leak LyC radiation in many directions. We found evidence of a kpc-scale fragmented superbubble, that may have cleared galactic-scale channels in the ISM. Additionally, the southwestern hemisphere is highly ionised in all velocities, hinting at a density bound scenario. A compact kpc-scale structure of lowly ionised gas coincides with the diffuse Lya emission and the presence of fast shocks. Finally, we find evidence that a significant fraction of the ionised gas mass may escape the gravitational potential of the galaxy.

X-ray radiation from accreting compact objects is an important part of stellar feedback. The metal-poor galaxy ESO 338-4 has experienced vigorous starburst during the last 40 Myr and contains some of the most massive super star clusters in the nearby Universe. Given its starburst age and its star-formation rate, ESO 338-4 is one of the most efficient nearby manufactures of neutron stars and black holes, hence providing an excellent laboratory for feedback studies. We compared X-ray images and spectra obtained by XMM-Newton and Chandra telescopes with integral field spectroscopic VLT MUSE observations in the optical to constrain the nature of strong X-ray emitters. X-ray observations uncover three ultraluminous X-ray sources (ULXs) in ESO 338-4. The brightest among them, ESO 338~X-1, has X-ray luminosity in excess of 10^{40} erg/s. We speculate that ESO 338-4 is powered by accretion on an intermediate-mass (~300Msun) black hole. We show that X-ray radiation from ULXs and hot superbubbles strongly contributes to HeII ionization and general stellar feedback in this template starburst galaxy.

The interplay between the formation of stars, stellar feedback and cloud properties strongly influences the star formation history of giant molecular clouds. The formation of massive stars leads to a variety of stellar clusters, ranging from low stellar density OB associations to dense, gravitationally bound starburst clusters. We aimed at identifying the massive stellar content and reconstructing the star formation history of the W51 giant molecular cloud. We performed near-infrared imaging and K-band spectroscopy of the massive stars in W51. We analyzed the stellar populations using colour-magnitude and colour-colour diagrams and compared the properties of the spectroscopically identified stars with stellar evolution models. We derived the ages of the different sub-clusters in W51 and, based on our spectroscopy derived an age for W51 of 3 Myrs or less. The age of the P Cygni star LS1 and the presence of two still forming proto-clusters suggests that the star formation history of W51 is more complex than a single burst. We did not find evidence for triggered star formation and we concluded that the star formation in W51 is multi seeded. We finally concluded that W51 is a OB association where different sub-clusters form over a time span of at least 3-5 Myrs.

Stellar feedback strongly affects the interstellar medium (ISM) of galaxies. Stellar feedback in the first galaxies likely plays a major role in enabling the escape of LyC photons, which contribute to the re-ionization of the Universe. Nearby starburst galaxies serve as local analogues allowing for a spatially resolved assessment of the feedback processes in these galaxies. We characterize the feedback effects from the star clusters in the local high-redshift analogue ESO 338-IG04 on the ISM and compare the results with the properties of the most massive clusters. We use high quality VLT/MUSE optical integral field data to derive the physical properties of the ISM such as ionization, density, shocks, and perform new fitting of the spectral energy distributions of the brightest clusters in ESO 338-IG04 from HST imaging. ESO 338-IG04 has a large ionized halo which we detect to a distance of 9 kpc. We identify 4 Wolf-Rayet (WR) clusters based on the blue and red WR bump. We follow previously identified ionization cones and find that the ionization of the halo increases with distance. Analysis of the galaxy kinematics shows two complex outflows driven by the numerous young clusters in the galaxy. We find a ring of shocked emission traced by an enhanced [OI]/H$\alpha$ ratio surrounding the starburst and at the end of the outflow. Finally we detect nitrogen enriched gas associated with the outflow, likely caused by the WR stars in the massive star clusters. Photo-ionization dominates the central starburst and sets the ionization structure of the entire halo, resulting in a density bounded halo, facilitating the escape of LyC photons. Outside the central starburst, shocks triggered by an expanding super bubble become important. The shocks at the end of the outflow suggest interaction between the hot outflowing material and the more quiescent halo gas.

We report on the discovery of ionised gas filaments in the circum-galactic halo of the extremely metal-poor compact starburst SBS 0335-052E in a 1.5h integration with the MUSE integral-field spectrograph. We detect these features in H${\alpha}$ and [OIII] emission down to surface-brightness levels of $5 \times 10^{-19}$erg s$^{-1}$cm$^{-2}$arcsec$^{-2}$. The filaments have projected diameters of 2.1 kpc and extend more than 9 kpc to the north and north-west from the main stellar body. We also detect extended nebular HeII $\lambda$4686 emission that brightens towards the north-west at the rim of a star-burst driven super-shell, suggestive of a locally enhanced UV radiation field due to shocks. We also present a velocity field of the ionised gas. The filaments appear to connect seamlessly in velocity space to the kinematical disturbances caused by the shell. Similar to high-$z$ star-forming galaxies, the ionised gas in this galaxy is dispersion dominated. We argue that the filaments were created via feedback from the starburst and that these ionised structures in the halo may act as escape channels for Lyman continuum radiation in this gas-rich system.

The formation process of massive stars is still poorly understood. Massive young stellar objects (mYSOs) are deeply embedded in their parental clouds, they are rare and thus typically distant, and their reddened spectra usually preclude the determination of their photospheric parameters. M17 is one of the best studied HII regions in the sky, is relatively nearby, and hosts a young stellar population. With X-shooter on the ESO Very Large Telescope we have obtained optical to near-infrared spectra of candidate mYSOs, identified by Hanson et al. (1997), and a few OB stars in this region. The large wavelength coverage enables a detailed spectroscopic analysis of their photospheres and circumstellar disks. We confirm the pre-main sequence (PMS) nature of six of the stars and characterise the O stars. The PMS stars have radii consistent with being contracting towards the main sequence and are surrounded by a remnant accretion disk. The observed infrared excess and the (double-peaked) emission lines provide the opportunity to measure structured velocity profiles in the disks. We compare the observed properties of this unique sample of young massive stars with evolutionary tracks of massive protostars by Hosokawa & Omukai (2009), and propose that these mYSOs near the western edge of the HII region are on their way to become main-sequence stars ($\sim 6 - 20$ $M_{\odot}$) after having undergone high mass-accretion rates (${\dot{M}_{\rm acc}} \sim 10^{-4} - 10^{-3}$ $M_{\odot}$ $\rm yr^{-1}$). Their spin distribution upon arrival at the zero age main sequence (ZAMS) is consistent with that observed for young B stars, assuming conservation of angular momentum and homologous contraction.

The formation of massive stars remains poorly understood and little is known about their birth multiplicity properties. Here, we investigate the strikingly low radial-velocity dispersion measured for a sample of 11 massive pre- and near-main-sequence stars (sigma_rv = 5.6 +/- 0.2 km/s) in the young massive star forming region M17 to obtain first constraints on the multiplicity properties of young massive stellar objects. Methods: We compute the RV dispersion of synthetic populations of massive stars for various multiplicity properties and we compare the simulated sigma_rv distributions to the observed value. We specifically investigate two scenarios: a low binary fraction and a dearth of short-period binary systems. Results: Simulated populations with low binary fractions (f_bin = 0.12_{-0.09}^{+0.16}) or with truncated period distributions (P_cutoff > 9 months) are able to reproduce the low sigma_rv observed within their 68%-confidence intervals. Parent populations with f_bin > 0.42 or P_cutoff < 47 d can however be rejected at the 5%-significance level. Both constraints are contrast with the high binary fraction and plethora of short-period systems found in few Myr-old, OB-type populations. To explain the difference, the first scenario requires a variation of the outcome of the massive star formation process. In the the second scenario, compact binaries must form later on, and the cut-off period may be related to physical length-scales representative of the bloated pre-main-sequence stellar radii or of their accretion disks. Conclusions: If the obtained constraints are representative of the overall properties of massive young stellar objects, our results may provide support to a formation process in which binaries are initially formed at larger separations, then harden or migrate to produce the typical (untruncated) power-law period distribution observed in few Myr-old OB binaries.

Massive stars form on different scales ranging from large, dispersed OB associations to compact, dense starburst clusters. The complex structure of regions of massive star formation, and the involved short timescales provide a challenge for our understanding of their birth and early evolution. As one of the most massive and luminous star-forming region in our Galaxy, W49 is the ideal place to study the formation of the most massive stars. By classifying the massive young stars deeply embedded into the molecular cloud of W49, we aim to investigate and trace the star formation history of this region. We analyse near-infrared $K$-band spectroscopic observations of W49 from LBT/LUCI combined with $JHK$ images obtained with NTT/SOFI and LBT/LUCI. Based on $JHK$-band photometry and K-band spectroscopy the massive stars are placed in a Hertzsprung Russell diagram. By comparison with evolutionary models, their age and hence the star formation history of W49 can be investigated. Fourteen O type stars as well as two young stellar objects (YSOs) are identified by our spectroscopic survey. Eleven O-stars are main sequence stars with subtypes ranging from O3 to O9.5, with masses ranging from $\sim 20$ $M_{\odot}$ to $\sim 120$ $M_{\odot}$. Three of the O-stars show strong wind features, and are considered to be Of-type supergiants with masses beyond 100 $M_{\odot}$ . The two YSOs show CO emission, indicative for the presence of circumstellar disks in the central region of the massive cluster. The age of the cluster is estimated as $\sim1.5$ Myr, with star formation still ongoing in different parts of the region. The ionising photons from the central massive stars have not yet cleared the molecular cocoon surrounding the cluster. W49 is comparable to extragalactic star-forming regions and provides us with an unique possibility to study a starburst in detail.

We present integral field spectroscopic observations with the Potsdam Multi Aperture Spectrophotometer of all 14 galaxies in the $z\sim 0.1$ Lyman Alpha Reference Sample (LARS). We produce 2D line of sight velocity maps and velocity dispersion maps from the Balmer $\alpha$ (H$\alpha$) emission in our data cubes. These maps trace the spectral and spatial properties of the LARS galaxies' intrinsic Ly$\alpha$ radiation field. We show our kinematic maps spatially registered onto the Hubble Space Telescope H$\alpha$ and Lyman $\alpha$ (Ly$\alpha$) images. Only for individual galaxies a causal connection between spatially resolved H$\alpha$ kinematics and Ly$\alpha$ photometry can be conjectured. However, no general trend can be established for the whole sample. Furthermore, we compute non-parametric global kinematical statistics -- intrinsic velocity dispersion $\sigma_0$, shearing velocity $v_\mathrm{shear}$, and the $v_\mathrm{shear}/\sigma_0$ ratio -- from our kinematic maps. In general LARS galaxies are characterised by high intrinsic velocity dispersions (54\,km\,s$^{-1}$ median) and low shearing velocities (65\,km\,s$^{-1}$ median). $v_\mathrm{shear}/\sigma_0$ values range from 0.5 to 3.2 with an average of 1.5. Noteworthy, five galaxies of the sample are dispersion dominated systems with $v_\mathrm{shear}/\sigma_0 <1$ and are thus kinematically similar to turbulent star forming galaxies seen at high redshift. When linking our kinematical statistics to the global LARS Ly$\alpha$ properties, we find that dispersion dominated systems show higher Ly$\alpha$ equivalent widths and higher Ly$\alpha$ escape fractions than systems with $v_\mathrm{shear}/\sigma_0 > 1$. Our result indicates that turbulence in actively star-forming systems is causally connected to interstellar medium conditions that favour an escape of Ly$\alpha$ radiation.

The leakage of Lyman continuum photons from star forming galaxies is an elusive parameter. When observed, it provides a wealth of information on star formation in galaxies and the geometry of the interstellar medium, and puts constraints on the role of star forming galaxies in the reionization of the universe. H-alpha-selected galaxies at z~2 trace the highest star formation population at the peak of cosmic star formation history, providing a base for directly measuring Lyman continuum escape. Here we present this method, and highlight its benefits as well as caveats. We also use the method on 10 H-alpha emitters in the Chandra Deep Field South at z=2.2, also imaged with the Hubble Space Telescope in the ultraviolet. We find no individual Lyman continuum detections, and our stack puts a 5 sigma upper limit on the average absolute escape fraction of <24%, consistent with similar studies. With future planned observations, the sample sizes would rapidly increase and the method presented here would provide very robust constraints on the escape fraction.

The stellar content of young massive star clusters emit large amounts of Lyman continuum photons and inject momentum into the inter stellar medium (ISM) by the strong stellar winds of the most massive stars in the cluster. When the most massive stars explode as supernovae, large amounts of mechanical energy are injected in the ISM. A detailed study of the ISM around these massive cluster provides insights on the effect of cluster feedback. We present high quality integral field spectroscopy taken with VLT/MUSE of two starburst galaxies: ESO 338-IG04 and Haro 11. Both galaxies contain a significant number of super star clusters. The MUSE data provide us with an unprecedented view of the state and kinematics of the ionized gas in the galaxy allowing us to study the effect of stellar feedback on small and large spatial scales. We present our recent results on studying the ISM state of these two galaxies. The data of both galaxies show that the mechanical and ionization feedback of the super star clusters in the galaxy modify the state and kinematics of the ISM substancially by creating highly ionized bubbles around the cluster, making the central part of the galaxy highly ionized. This shows that the HII regions around the individual clusters are density bounded, allowing the ionizing photons to escape and ionize the ISM further out.

The Ly$\alpha$ line is an important diagnostic for star formation at high redshift, but interpreting its flux and line profile is difficult because of the resonance nature of Ly$\alpha$. Trends between the escape of Ly$\alpha$ photons and dust and properties of the interstellar medium (ISM) have been found, but detailed comparisons between Ly$\alpha$ emission and the properties of the gas in local high-redshift analogs are vital for understanding the relation between Ly$\alpha$ emission and galaxy properties. For the first time, we can directly infer the properties of the ionized gas at the same location and similar spatial scales of the extended Ly$\alpha$ halo around ESO 338-IG04. We obtained VLT/MUSE integral field spectra. We used ionization parameter mapping of the [SII]/[OIII] line ratio and the kinematics of H$\alpha$ to study the ionization state and kinematics of the ISM of ESO338-IG04. The velocity map reveals two outflows. The entire central area of the galaxy is highly ionized by photons leaking from the HII regions around the youngest star clusters. Three highly ionized cones have been identified, of which one is associated with an outflow. We propose a scenario where the outflows are created by mechanical feedback of the older clusters, while the highly ionized gas is caused by the hard ionizing photons emitted by the youngest clusters. A comparison with the Ly$\alpha$ map shows that the (approximately bipolar) asymmetries observed in the Ly$\alpha$ emission are consistent with the base of the outflows detected in H$\alpha$. No clear correlation with the ionization cones is found. The mechanical and ionization feedback of star clusters significantly changes the state of the ISM by creating ionized cones and outflows. The comparison with Ly$\alpha$ suggests that especially the outflows could facilitate the escape of Ly$\alpha$ photons [Abridged].

Very massive stars (M>100 M$_{\odot}$) are very rare objects, but have a strong influence on their environment. The formation of this kind of objects is of prime importance in star formation, but observationally still poorly constrained. We report on the identification of a very massive star in the central cluster of the star-forming region W49. We investigate near-infrared K-band spectroscopic observations of W49 from VLT/ISAAC together with JHK images obtained with NTT/SOFI and LBT/LUCI. We derive a spectral type of W49nr1, the brightest star in the dense core of the central cluster of W49. On the basis of its K-band spectrum, W49nr1 is classified as an O2-3.5If* star with a K-band absolute magnitude of -6.27$\pm$0.10 mag. The effective temperature and bolometric correction are estimated from stars of similar spectral type. After comparison to the Geneva evolutionary models, we find an initial mass between 100 M$_{\odot}$ and 180 M$_{\odot}$. Varying the extinction law results in a larger initial mass range of 90 - 250 M$_{\odot}$.

Embedded clusters like W3 Main are complex and dynamically evolving systems that represent an important phase of the star formation process. We aim at the characterization of the entire stellar content of W3 Main in a statistical sense to identify possible differences in evolutionary phase of the stellar populations and find clues about the formation mechanism of this massive embedded cluster. Methods. Deep JHKs imaging is used to derive the disk fraction, Ks-band luminosity functions and mass functions for several subregions in W3 Main. A two dimensional completeness analysis using artificial star experiments is applied as a crucial ingredient to assess realistic completeness limits for our photometry. We find an overall disk fraction of 7.7 $\pm$ 2.3%, radially varying from 9.4 $\pm$ 3.0 % in the central 1 pc to 5.6 $\pm$ 2.2 % in the outer parts of W3 Main. The mass functions derived for three subregions are consistent with a Kroupa and Chabrier mass function. The mass function of IRSN3 is complete down to 0.14 Msun and shows a break at M $\sim$ 0.5 Msun. We interpret the higher disk fraction in the center as evidence for a younger age of the cluster center. We find that the evolutionary sequence observed in the low-mass stellar population is consistent with the observed age spread among the massive stars. An analysis of the mass function variations does not show evidence for mass segregation. W3 Main is currently still actively forming stars, showing that the ionizing feedback of OB stars is confined to small areas ($\sim$ 0.5 pc). The FUV feedback might be influencing large regions of the cluster as suggested by the low overall disk fraction.

Massive stars play a dominant role in the process of clustered star formation, with their feedback into the molecular cloud through ionizing radiation, stellar winds and outflows. The formation process of massive stars is poorly constrained because of their scarcity, the short formation timescale and obscuration. By obtaining a census of the newly formed stellar population, the star formation history of the young cluster and the role of the massive stars within it can be unraveled. We aim to reconstruct the formation history of the young stellar population of the massive star-forming region RCW 36. We study several dozens of individual objects, both photometrically and spectroscopically, look for signs of multiple generations of young stars and investigate the role of the massive stars in this process. We obtain a census of the physical parameters and evolutionary status of the young stellar population. Using a combination of near-infrared photometry and spectroscopy we estimate ages and masses of individual objects. We identify the population of embedded young stellar objects (YSO) by their infrared colors and emission line spectra. RCW 36 harbors a stellar population of massive and intermediate-mass stars located around the center of the cluster. Class 0/I and II sources are found throughout the cluster. The central population has a median age of 1.1 +/- 0.6 Myr. Of the stars which could be classified, the most massive ones are situated in the center of the cluster. The central cluster is surrounded by filamentary cloud structures; within these, some embedded and accreting YSOs are found. Our age determination is consistent with the filamentary structures having been shaped by the ionizing radiation and stellar winds of the central massive stars. The formation of a new generation of stars is ongoing, as demonstrated by the presence of embedded protostellar clumps, and two exposed jets.

The spatial variations in PAH band intensities are normally attributed to the physical conditions of the emitting PAHs, however in recent years it has been suggested that such variations are caused mainly by extinction. To resolve this question, we have obtained near-infrared (NIR), mid-infrared (MIR) and radio observations of the compact H ii region IRAS 12063-6259. We use these data to construct multiple independent extinction maps and also to measure the main PAH features (6.2, 7.7, 8.6 and 11.2 {\mu}m) in the MIR. Three extinction maps are derived: the first using the NIR hydrogen lines and case B recombination theory; the second combining the NIR data with radio data; and the third making use of the Spitzer/IRS MIR observations to measure the 9.8 {\mu}m silicate absorption feature using the Spoon method and PAHFIT (as the depth of this feature can be related to overall extinction). The silicate absorption over the bright, southern component of IRAS 12063-6259 is almost absent while the other methods find significant extinction. While such breakdowns of the relationship between the NIR extinction and the 9.8 {\mu}m absorption have been observed in molecular clouds, they have never been observed for HII regions. We then compare the PAH intensity variations in the Spitzer/IRS data after dereddening to those found in the original data. It was found that in most cases, the PAH band intensity variations persist even after dereddening, implying that extinction is not the main cause of the PAH band intensity variations.

We aim to take a census of molecular hydrogen emission line objects (MHOs) in the Rho Ophiuchi molecular cloud and to make the first systematic proper motion measurements of these objects in this region. Deep H2 near-infrared imaging is performed to search for molecular hydrogen emission line objects. Multi-epoch data are used to derive the proper motions of the features of these objects, and the lengths and opening angles of the molecular hydrogen outflows. Our imaging covers an area of about 0.11 deg2 toward the L1688 core in the Rho Ophiuchi molecular cloud. In total, six new MHOs are discovered, 32 previously known MHOs are detected, and the proper motions for 86 features of the MHOs are measured. The proper motions lie in the range of 14 to 247 mas/yr, corresponding to transversal velocities of 8 to 140 km s-1 with a median velocity of about 35 km s-1. Based on morphology and proper motion measurements, 27 MHOs are ascribed to 21 driving sources. The molecular hydrogen outflows have a median length of about 0.04 pc and random orientations. We find no obvious correlation between H2 jet length, jet opening angle, and the evolutionary stage of the driving sources as defined by their spectral indices. We find that the fraction of protostars (23%) that drive molecular hydrogen outflows is similar to the one for Class II sources (15%). For most molecular hydrogen outflows, no obvious velocity variation along the outflow has been found. In Ophiuchus the frequency of occurrence of molecular hydrogen outflows has no strong dependency on the evolutionary stage of the driving source during the evolution from the protostellar stage to the classical T Tauri stage.

We present near-infrared JHKs imaging as well as K-band multi-object spectroscopy of the massive stellar content of W3 Main using LUCI at the LBT. We confirm 13 OB stars by their absorption line spectra in W3 Main and spectral types between O5V and B4V have been found. Three massive Young Stellar Objects are identified by their emission line spectra and near-infrared excess. From our spectrophotometric analysis of the massive stars and the nature of their surrounding HII regions we derive the evolutionary sequence of W3 Main and we find an age spread of 2-3 Myr.

We observed three high-mass star-forming regions in the W3 high-mass star formation complex with the Submillimeter Array and IRAM 30 m telescope. These regions, i.e. W3 SMS1 (W3 IRS5), SMS2 (W3 IRS4) and SMS3, are in different evolutionary stages and are located within the same large-scale environment, which allows us to study rotation and outflows as well as chemical properties in an evolutionary sense. While we find multiple mm continuum sources toward all regions, these three sub-regions exhibit different dynamical and chemical properties, which indicates that they are in different evolutionary stages. Even within each subregion, massive cores of different ages are found, e.g. in SMS2, sub-sources from the most evolved UCHII region to potential starless cores exist within 30 000 AU of each other. Outflows and rotational structures are found in SMS1 and SMS2. Evidence for interactions between the molecular cloud and the HII regions is found in the 13CO channel maps, which may indicate triggered star formation.

Recent studies of massive-star forming regions indicate that they can contain multiple generations of young stars. These observations suggest that star formation in these regions is sequential and/or triggered by a previous generation of (massive) stars. Here we present new observations of the star forming region RCW 36 in the Vela Molecular Ridge, hosting a young cluster of massive stars embedded in a molecular cloud complex. In the periphery of the cluster several young stellar objects (YSOs) are detected which produce bipolar jets (HH 1042 and HH 1043) demonstrating that these objects are still actively accreting. The VLT/X-shooter spectrum of the jet structure of HH 1042 provides detailed information on the physical conditions and kinematical properties of the jet plasma. From this information the YSO's accretion history can be derived. Combining the photometric and spectroscopic observations of RCW 36 gives insight into the formation process of individual stars and the star formation history of this young massive-star forming region.

We have started a campaign to identify massive star clusters inside bright molecular bubbles towards the Galactic Center. The CN15/16/17 molecular complex is the first example of our study. The region is characterized by the presence of two young clusters, DB10 and DB11, visible in the NIR, an ultra-compact HII region identified in the radio, several young stellar objects visible in the MIR, a bright diffuse nebulosity at 8\mu m coming from PAHs and sub-mm continuum emission revealing the presence of cold dust. Given its position on the sky (l=0.58, b=-0.85) and its kinematic distance of ~7.5 kpc, the region was thought to be a very massive site of star formation in proximity of the CMZ. The cluster DB11 was estimated to be as massive as 10^4 M_sun. However the region's properties were known only through photometry and its kinematic distance was very uncertain given its location at the tangential point. We aimed at better characterizing the region and assess whether it could be a site of massive star formation located close to the Galactic Center. We have obtained NTT/SofI JHKs photometry and long slit K band spectroscopy of the brightest members. We have additionally collected data in the radio, sub-mm and mid infrared, resulting in a quite different picture of the region. We have confirmed the presence of massive early B type stars and have derived a spectro-photometric distance of ~1.2 kpc, much smaller than the kinematic distance. Adopting this distance we obtain clusters masses of M(DB10) ~ 170 M_sun and M(DB11) ~ 275 M_sun. This is consistent with the absence of any O star, confirmed by the excitation/ionization status of the nebula. No HeI diffuse emission is detected in our spectroscopic observations at 2.113\mu m, which would be expected if the region was hosting more massive stars. Radio continuum measurements are also consistent with the region hosting at most early B stars.

We present near-infrared multi-object spectroscopy and JHKs imaging of the massive stellar content of the Galactic star-forming region W3 Main, obtained with LUCI at the Large Binocular Telescope. We confirm 15 OB stars in W3 Main and derive spectral types between O5V and B4V from their absorption line spectra. Three massive Young Stellar Objects are identified by their emission line spectra and near-infrared excess. The color-color diagram of the detected sources allows a detailed investigation of the slope of the near-infrared extinction law towards W3 Main. Analysis of the Hertzsprung Russell diagram suggests that the Nishiyama extinction law fits the stellar population of W3 Main best (E(J-H)/E(H-Ks) = 1.76 and R_(Ks) = 1.44). From our spectrophotometric analysis of the massive stars and the nature of their surrounding HII regions we derive the evolutionary sequence of W3 Main and we find evidence of an age spread of at least 2-3 Myr. While the most massive star (IRS2) is already evolved, indications for high-mass pre-main-sequence evolution is found for another star (IRS N1), deeply embedded in an ultra compact HII region, in line with the different evolutionary phases observed in the corresponding HII regions. We derive a stellar mass of W3 Main of (4 +- 1) 10^3 Msun, by extrapolating from the number of OB stars using a Kroupa IMF and correcting for our spectroscopic incompleteness. We have detected the photospheres of OB stars from the more evolved diffuse HII region to the much younger UCHII regions, suggesting that these stars have finished their formation and cleared away their circumstellar disks very fast. Only in the hyper-compact HII region (IRS5), the early type stars seem to be still surrounded by circumstellar material.

We present near-infrared JHKs imaging as well as K-band multi-object spectroscopy of the massive stellar content of W3 Main using LUCI at the LBT. We confirm 13 OB stars by their absorption line spectra in W3 Main and spectral types between O5V and B4V have been found. Three massive Young Stellar Objects are identified by their emission line spectra and near-infrared excess. From our spectrophotometric analysis of the massive stars and the nature of their surrounding HII regions we derive the evolutionary sequence of W3 Main and we find evidence of an age spread of at least 2-3 Myr. While the most massive star (IRS2) is already evolved, indications for high-mass pre--main-sequence evolution is found for another star (IRS N1), deeply embedded in an ultra compact HII region, in line with the different evolutionary phases observed in the corresponding HII regions. We have detected the photospheres of OB stars from the more evolved diffuse HII region to the much younger UCHII regions, suggesting that the OB stars have finished their formation and cleared away their possible circumstellar disks very fast. Only in the hyper-compact HII region (IRS5), the early type stars are still surrounded by circumstellar material.

Context. The formation and early evolution of high- and intermediate-mass stars towards the main sequence involves the interplay of stars in a clustered and highly complex environment. To obtain a full census of this interaction, the Formation and Early evolution of Massive Stars (FEMS) collaboration studies a well-selected sample of 10 high-mass star-forming regions. Aims. In this study we examine the stellar content of the high-mass star-forming region centered on IRAS 06084-0611 in the Monoceros R2 cloud. Methods. Using the near-infrared H- and K-band spectra from the VLT/SINFONI instrument on the ESO Very Large Telescope (VLT)and photometric near-infrared NTT/SOFI, 2MASS and Spitzer/IRAC data, we were able to determine the spectral types for the most luminous stars in the cluster. Results. Two very young and reddened massive stars have been detected by SINFONI: a massive Young Stellar Object (YSO) con- sistent with an early-B spectral type and a Herbig Be star. Furthermore, stars of spectral type G and K are detected while still in the Pre-Main Sequence (PMS) phase. We derive additional properties such as temperatures, extinctions, radii and masses. We present a Hertzsprung-Russell diagram and find most objects having intermediate masses between \sim1.5-2.5 M\odot. For these stars we derive a median cluster age of \sim4 Myr. Conclusions. Using Spitzer/IRAC data we confirm earlier studies that the younger class 0/I objects are centrally located while the class II objects are spread out over a larger area, with rough scale size radii of \sim0.5 pc and \sim1.25 pc respectively. Moreover, the presence of a massive YSO, an ultracompact H ii region and highly reddened objects in the center of the cluster suggest a much younger age of < 1 Myr. A possible scenario for this observation would be sequential star formation along the line of sight; from a cluster of intermediate-mass to high-mass stars.

Context: High-mass star formation has been a very active field over the last decade, however, most studies targeted regions of luminosities between 10^4 and 10^5 L_sun. Methods: We selected the W31 star-forming complex with a total luminosity of ~6x10^6 L_sun for a multi-wavelength spectral line and continuum study covering wavelengths from the near- and mid-infrared via (sub)mm wavelength observations to radio data in the cm regime. Results: While the overall structure of the multi-wavelength continuum data resembles each other well, there are several intriguing differences. The 24mum emission stemming largely from small dust grains follows tightly the spatial structure of the cm emission tracing the ionized free-free emission. Hence warm dust resides in regions that are spatially associated with the ionized hot gas (~10^4 K) of the HII regions. Furthermore, we find several evolutionary stages within the same complexes, ranging from infrared-observable clusters, via deeply embedded regions associated with active star formation traced by 24\,$\mu$m and cm emission, to at least one high-mass gas clump devoid of any such signature. The 13CO(2-1) and C18O(2-1) spectral line observations reveal a large kinematic breadth in the entire region with a total velocity range of approximately 90 km/s. While the average virial mass ratio for W31 is close to unity, the line width analysis indicates large-scale evolutionary differences between the southern and northern sub-regions (G10.2-0.3 and G10.3-0.1) of the whole W31 complex. The clump mass function - tracing cluster scales and not scales of individual stars - derived from the 875mum continuum data has a slope of 1.5+-0.3, consistent with previous cloud mass functions.

We present observations of the embedded massive young stellar object (YSO) candidate 08576nr292, obtained with X-shooter and SINFONI on the ESO Very Large Telescope (VLT). The flux-calibrated, medium-resolution X-shooter spectrum (300 - 2500 nm) includes over 300 emission lines, but no (photospheric) absorption lines and is consistent with a reddened disk spectrum. Among the emission lines are three hydrogen series and helium lines, both permitted and forbidden metal lines, and CO first-overtone emission. A representative sample of lines with different morphologies is presented. The H alpha and Ca II triplet lines are very strong, with profiles indicative of outflow and - possibly - infall, usually observed in accreting stars. These lines include a blue-shifted absorption component at ~ -125 km/s. The He I and metal-line profiles are double-peaked, with a likely origin in a circumstellar disk. The forbidden lines, associated with outflow, have a single blue-shifted emission component centered at -125 km/s, coinciding with the absorption components in H alpha and Ca II. SINFONI H- and K-band integral-field spectroscopy of the cluster environment demonstrates that the [Fe II] emission is produced by a jet originating at the location of 08576nr292. Because the spectral type of the central object cannot be determined, its mass remains uncertain. We argue that 08576nr292 is an intermediate-mass YSO with a high accretion rate Mdot_acc ~ 10^-5 - 10^-6 Msun/yr. These observations demonstrate the potential of X-shooter and SINFONI to study in great detail an accretion disk-jet system, rarely seen around the more massive YSOs.

We used the near-IR imager/spectrograph LUCIFER mounted on the Large Binocular Telescope (LBT) to image, with sub-arcsec seeing, the local dwarf starburst NGC 1569 in the JHK bands and HeI 1.08 micron, [FeII] 1.64 micron and Brgamma narrow-band filters. We obtained high-quality spatial maps of HeI, [FeII] and Brgamma emission across the galaxy, and used them together with HST/ACS images of NGC 1569 in the Halpha filter to derive the two-dimensional spatial map of the dust extinction and surface star formation rate density. We show that dust extinction is rather patchy and, on average, higher in the North-West (NW) portion of the galaxy [E_g(B-V) = 0.71 mag] than in the South-East [E_g(B-V) = 0.57 mag]. Similarly, the surface density of star formation rate peaks in the NW region of NGC 1569, reaching a value of about 4 x 10^-6 M_sun yr^-1 pc^-2. The total star formation rate as estimated from the integrated, dereddened Halpha luminosity is about 0.4 M_sun yr^-1, and the total supernova rate from the integrated, dereddened [FeII] luminosity is about 0.005 yr^-1 (assuming a distance of 3.36 Mpc). The azimuthally averaged [FeII]/Brgamma flux ratio is larger at the edges of the central, gas-deficient cavities (encompassing the super star clusters A and B) and in the galaxy outskirts. If we interpret this line ratio as the ratio between the average past star formation (as traced by supernovae) and on-going activity (represented by OB stars able to ionize the interstellar medium), it would then indicate that star formation has been quenched within the central cavities and lately triggered in a ring around them. The number of ionizing hydrogen and helium photons as computed from the integrated, dereddened Halpha and HeI luminosities suggests that the latest burst of star formation occurred about 4 Myr ago and produced new stars with a total mass of ~1.8 x 10^6 M_sun. [Abridged]

To understand evolutionary and environmental effects during the formation of high-mass stars, we observed three regions of massive star formation at different evolutionary stages that reside in the same natal molecular cloud. Methods. The three regions S255IR, S255N and S255S were observed at 1.3 mm with the Submillimeter Array (SMA) and followup short spacing information was obtained with the IRAM 30m telescope. Near infrared (NIR) H + K-band spectra and continuum observations were taken for S255IR with VLT-SINFONI to study the different stellar populations in this region. The combination of millimeter (mm) and near infrared data allow us to characterize different stellar populations within the young forming cluster in detail. While we find multiple mm continuum sources toward all regions, their outflow, disk and chemical properties vary considerably. The most evolved source S255IR exhibits a collimated bipolar outflow visible in CO and H2 emission, the outflows from the youngest region S255S are still small and rather confined in the regions of the mm continuum peaks. Also the chemistry toward S255IR is most evolved exhibiting strong emission from complex molecules, while much fewer molecular lines are detected in S255N, and in S255S we detect only CO isotopologues and SO lines. Also, rotational structures are found toward S255N and S255IR. Furthermore, a comparison of the NIR SINFONI and mm data from S255IR clearly reveal two different (proto) stellar populations with an estimated age difference of approximately 1 Myr. A multi-wavelength spectroscopy and mapping study reveals different evolutionary phases of the star formation regions. We propose the triggered outside-in collapse star formation scenario for the bigger picture and the fragmentation scenario for S255IR.

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.

We present VLT/SINFONI integral field spectroscopy of RCW 34 along with Spitzer/IRAC photometry of the surroundings. RCW 34 consists of three different regions. A large bubble has been detected on the IRAC images in which a cluster of intermediate- and low-mass class II objects is found. At the northern edge of this bubble, an HII region is located, ionized by 3 OB stars. Intermediate mass stars (2 - 3 Msun) are detected of G- and K- spectral type. These stars are still in the pre-main sequence (PMS) phase. North of the HII region, a photon-dominated region is present, marking the edge of a dense molecular cloud traced by H2 emission. Several class 0/I objects are associated with this cloud, indicating that star formation is still taking place. The distance to RCW 34 is revised to 2.5 +- 0.2 kpc and an age estimate of 2 - 1 Myrs is derived from the properties of the PMS stars inside the HII region. The most likely scenario for the formation of the three regions is that star formation propagates from South to North. First the bubble is formed, produced by intermediate- and low-mass stars only, after that, the HII region is formed from a dense core at the edge of the molecular cloud, resulting in the expansion as a champagne flow. More recently, star formation occurred in the rest of the molecular cloud. Two different formation scenarios are possible: (a) The bubble with the cluster of low- and intermediate mass stars triggered the formation of the O star at the edge of the molecular cloud which in turn induces the current star-formation in the molecular cloud. (b) An external triggering is responsible for the star-formation propagating from South to North. [abridged]

IRAS19410+2336 is a young massive star forming region with an intense outflow activity. We present here spatially resolved NIR spectroscopy which allows us to verify whether the H2 emission detected in this object originates from thermal emission in shock fronts or from fluorescence excitation by non-ionizing UV photons. Moreover, NIR spectroscopy also offers the possibility of studying the characteristics of the putative driving source(s) of the H2 emission by the detection of photospheric and circumstellar spectral features, and of the environmental conditions (e.g. extinction). We obtained long-slit, intermediate-resolution, NIR spectra of IRAS19410+2336 using LIRIS. As a complement, we also obtained J, H and K_s images with the Las Campanas 2.5m Du Pont Telescope, and archival mid-infrared (MIR) Spitzer images at 3.6, 4.5, 5.8 and 8.0 um. We confirm the shocked nature of the H2 emission, with an excitation temperature of about 2000 K. We have also identified objects with very different properties and evolutionary stages in IRAS19410+2336. The most massive source at millimeter wavelengths, mm1, with a mass of a few tens of solar masses, has a bright NIR (and MIR) counterpart. This suggests that emission is leaking at these wavelengths. The second most massive millimeter source, mm2, is only detected at lambda > 6 um, suggesting that it could be a high-mass protostar still in its main accretion phase. The NIR spectra of some neighboring sources show CO first-overtone bandhead emission which is associated with neutral material located in the inner regions of the circumstellar environment of YSOs.

Near-infrared imaging surveys of high-mass star-forming regions reveal an amazingly complex interplay between star formation and the environment (Churchwell et al. 2006; Alvarez et al. 2004). By means of near-IR spectroscopy the embedded massive young stars can be characterized and placed in the context of their birth site. However, so far spectroscopic surveys have been hopelessly incomplete, hampering any systematic study of these very young massive stars. New integral field instrumentation available at ESO has opened the possibility to take a huge step forward by obtaining a full spectral inventory of the youngest massive stellar populations in star-forming regions currently accessible. Simultaneously, the analysis of the extended emission allows the characterization of the environmental conditions. The Formation and Early Evolution of Massive Stars (FEMS) collaboration aims at setting up a large observing campaign to obtain a full census of the stellar content, ionized material, outflows and PDR's over a sample of regions that covers a large parameter space. Complementary radio, mm and infrared observations will be used for the characterization of the deeply embedded population. For the first eight regions we have obtained 40 hours of SINFONI observations. In this contribution, we present the first results on three regions that illustrate the potential of this strategy.

There is considerable controversy surrounding the nature of M1-78, a compact nebula located beyond the Perseus arm. It was first classified as a planetary nebula and is nowadays generally considered to be a compact HII region. To investigate the nature M1-78 further, we present a detailed spectroscopic study of M1-78 in the optical and near-infrared. M1-78 is a high-density nebula with substantial physical differences between its two main morphological zones: a bright arc to the SW and a blob of emission in the NE. Specifically, the blob in the NE has a higher electron temperature and visual extinction than the SW arc. The most important result, however, is the confirmation of a nitrogen enrichment in M1-78. This enrichment is stronger at the location of the NE blob and is correlated with a defficiency in the O abundance and a (dubious) He enrichment. Such an abundance pattern is typical of ejecta nebulae around evolved massive stars such as Wolf-Rayet and Luminous Blue Variable stars. The spatial variations in the physical conditions and chemical abundances and the presence of more than one possible ionizing source indicates, however, that M1-78 is better described as a combination of a compact HII region + ejecta. Finally, we detect H2 emission that extends over a large (~30 arcsec) area around the ionized nebula. Analysis of the near-infrared H2 lines indicates that the excitation mechanism is UV fluorescence.

Near-infrared surveys of high-mass star-forming regions start to shed light onto their stellar content. A particular class of objects found in these regions, the so-called massive Young Stellar Objects (YSOs) are surrounded by dense circumstellar material. Several near- and mid-infrared diagnostic tools are used to infer the physical characteristics and geometry of this circumstellar matter. Near-infrared hydrogen emission lines provide evidence for a disk-wind. The profiles of the first overtone of the CO band-heads, originating in the inner 10 AU from the central star, are well fitted assuming a keplerian rotating disk. The mid-infrared spectral energy distribution requires the presence of a more extended envelope containing dust at a temperature of about 200 K. CRIRES observations of CO fundamental absorption lines confirm the presence of a cold envelope. We discuss the evolutionary status of these objects.

We present the first multi-epoch radial velocity study of embedded young massive stars using near-infrared spectra obtained with ISAAC mounted at the ESO Very Large Telescope, with the aim to detect massive binaries. Our 16 targets are located in high-mass star-forming regions and many of them are associated to known ultracompact HI I regions, whose young age ensures that dynamic evolution of the clusters did not influence the intrinsic binarity rate. We identify two stars with about 90 km/s velocity differences between two epochs proving the presence of close massive binaries. The fact that 2 out of the 16 observed stars are binary systems suggests that at least 20% of the young massive stars are formed in close multiple systems, but may also be consistent with most if not all young massive stars being binaries. In addition, we show that the radial velocity dispersion of the full sample is about 35 km/s, significantly larger than our estimated uncertainty (25 km/s). This finding is consistent with similar measurements of the young massive cluster 30 Dor which might have a high intrinsic binary rate. Furthermore, we argue that virial cluster masses derived from the radial velocity dispersion of young massive stars may intrinsically overestimate the cluster mass due to the presence of binaries.

We have obtained high resolution (R = 10,000) K-band spectra of candidate young massive stars deeply embedded in high-mass star-forming regions. These objects were selected from a near-infrared survey of 44 regions of high-mass star-formation (Kaper et al, 2006). In these clusters, 38 OB stars are identified whose K-band spectra are dominated by photospheric emission. In almost all those stars, the K-band spectra are indistinguishable from field stars. However, in some stars the profile of the Bracket gamma line is different (less deep, or absent) from those of the O field stars. One of the explanations of these profiles might be an enhanced mass-loss.

High-quality K-band spectra of strongly reddened point sources, deeply embedded in (ultra-) compact HII regions, have revealed a population of 20 young massive stars showing no photospheric absorption lines, but sometimes strong Br-gamma emission. The Br-gamma equivalent widths occupy a wide range (from about 1 to over 100 A); the line widths of 100-200 km/s indicate a circumstellar rather than a nebular origin. The K-band spectra exhibit one or more features commonly associated with massive young stellar objects (YSOs) surrounded by circumstellar material: a very red colour (J-K) > 2, CO bandhead emission, hydrogen emission lines (sometimes doubly peaked), and FeII and/or MgII emission lines. The massive YSO distribution in the CMD suggests that the majority of the objects are of similar spectral type as the Herbig Be stars, but some of them are young O stars. The CO emission must come from a relatively dense (~10^{10} cm^{-3}) and hot (T~ 2000-5000 K) region, sufficiently shielded from the intense UV radiation field of the young massive star. The hydrogen emission is produced in an ionised medium exposed to UV radiation. The best geometrical solution is a dense and neutral circumstellar disk causing the CO bandhead emission, and an ionised upper layer where the hydrogen lines are produced. We present arguments that the circumstellar disk is more likely a remnant of the accretion process than the result of rapid rotation and mass loss such as in Be/B[e] stars.

High-quality K-band spectra of point sources, deeply embedded in massive star-forming regions, have revealed a population of 20 young massive stars showing no photospheric absorption lines, but only emission lines. The K-band spectra exhibit one or more features commonly associated with massive Young Stellar Objects surrounded by circumstellar material: a very red color (J-K) = 2, CO bandhead emission, hydrogen emission lines (sometimes doubly peaked), and FeII and/or MgII emission lines. The CO emission comes from a relatively dense (~10^10 cm^(-3)) and hot (T ~ 2000-5000 K) region, sufficiently shielded from the intense UV radiation field of the young massive star. Modeling of the CO-first overtone emission shows that the CO gas is located within 5 AU of the star. The hydrogen emission is produced in an ionized medium exposed to UV radiation. The best geometrical configuration is a dense and neutral circumstellar disk causing the CO bandhead emission, and an ionized upper layer where the hydrogen lines are produced. We argue that the circumstellar disk is likely a remnant of the accretion via a circumstellar disk.

We have obtained high resolution (R = 10,000) K-band spectra of candidate young massive stars deeply embedded in (ultra-) compact HII regions (UCHIIs). In these clusters, three types of objects are identified. The first type (38 objects) consists of ``naked'' OB stars whose K-band spectra are dominated by photospheric emission. We classify the K-band spectra of the OB-type cluster members using near-infrared classification criteria. The spectral classification provides an important constraint on the distance to the embedded cluster. The ionising power of the population thus derived is compared to the information obtained from the infrared and radio flux of these sources. In most cases these two different determinations of the ionising flux are consistent, from which we conclude that we have identified the ionising star(s) in about 50% of the embedded clusters. The second type (7 objects) are point sources associated with UCHII radio emission, that exhibit nebular emission lines in the near-infrared. Six of the objects in this group produce HeI emission indicative of an embedded O-type star. These objects are more embedded than the OB stars and probably do not dominate the infrared flux as measured by IRAS. The third type (20 objects) is characterised by broad (100--200 km/s) Brgamma emission and no photospheric absorption profiles. Bik et al (2005) show that these objects are massive YSO candidates surrounded by dense circumstellar disks.

We report on the observation of hot water vapor (steam) in the inner AU of a young massive star located in the star-forming region IRAS 08576-4334. The water lines are detected in a medium resolution (R=10,000) K-band spectrum taken by the infrared spectrometer ISAAC mounted on the VLT-ANTU. The water vapor is at a mean temperature of 1565+/-510 K, cooler than the hot CO found in the same object, which is at 1660K and the column density is N(H2O)=(2.5 +/- 0.4)x1E18 cm-2. The profile of both H2O and CO lines is best reproduced by the emission from a Keplerian disk. To interpret the data, we also investigate the formation of molecules and especially CO and water vapor in the inner hot and dense part of disks around young high mass stars using a pseudo time-dependent gas-phase chemical model. Molecules are rapidly photodissociated but this destruction is compensated by an efficient formation due to fast neutral-neutral reactions. The ability of CO molecules to self-shield significantly enhances its abundance. Water molecules are sufficiently abundant to be detectable. The observed H2O/CO ratio is reproduced by gas at 1600K and an enhanced UV field over gas density ratio I_UV/nH=1E(-4)- 1E(-6). The simulations support the presence of CO and H2O molecules in the inner disks around young massive stars despite the strong UV radiation and show that the OH radical plays an essential role in hot gas chemistry.

We report on the observation of CO bandhead emission around 51 Oph (dv=2). A high resolving power (R~10,000) spectrum was obtained with the infrared spectrometer ISAAC mounted on VLT-ANTU. Modeling of the profile suggests that the hot (Tgas=2000-4000K) and dense (n_H>1e10 cm-3) molecular material as probed by the CO bandhead is located in the inner AU of a Keplerian disk viewed almost edge-on. Combined with the observation of cooler gas (Tgas=500-900K) by ISO-SWS and the lack of cold material, our data suggest that the disk around 51 Oph is essentially warm and small. We demonstrate the presence of a dust-free inner disk that extents from the inner truncation radius until the dust sublimation radius. The disk around 51 Oph may be in a rare transition state toward a small debris disk object.

We present observations of CO overtone bandhead emission toward four massive Young Stellar Objects (spectral type O6--B5). The high signal-to-noise ratio K-band spectra were obtained with VLT-ISAAC at a resolution of 30 km/s, sufficient to resolve the bandheads, but not the individual J-lines. We are able to explain the shape of the lines by assuming a simple isothermal keplerian disk model seen at different inclinations. The gas temperature ranges from 1500 to 4500 K and the CO column density is between 0.1 and 4 10^21 cm^-2. The emission probably arises within the first few astronomical units of the disk, consistent with the high gas temperature. Our results indicate that molecules can survive close to a hot star and suggest that dense (n_H > 10^10 cm^-3) inner disks may be relatively common at an advanced stage of high-mass star formation.

We re-examine the nature of NGC2024-IRS2 in light of the recent discovery of the late O-type star, IRS2b, located 5 arcsec from IRS2. Using L-band spectroscopy, we set a lower limit of Av = 27.0 mag on the visual extinction towards IRS2. Arguments based on the nature of the circumstellar material, favor an Av of 31.5 mag. IRS2 is associated with the UCHII region G206.543-16.347 and the infrared source IRAS 05393-0156. We show that much of the mid-infrared emission towards IRS2, as well as the far infrared emission peaking at ~ 100 micron, do not originate in the direct surroundings of IRS2, but instead from an extended molecular cloud. Using new K-, L- and L'-band spectroscopy and a comprehensive set of infrared and radio continuum measurements from the literature, we apply diagnostics based on the radio slope, the strength of the infrared hydrogen recombination lines, and the presence of CO band-heads to constrain the nature and spatial distribution of the circumstellar material of IRS2. Using simple gaseous and/or dust models of prescribed geometry, we find strong indications that the infrared flux originating in the circumstellar material of IRS2 is dominated by emission from a dense gaseous disk with a radius of about 0.6 AU. At radio wavelengths the flux density distribution is best described by a stellar wind recombining at a radius of about 100 AU. Although NGC2024/IRS2 shares many similarities with BN-like objects, we do not find evidence for the presence of a dust shell surrounding this object. Therefore, IRS2 is likely more evolved.

A high quality, medium-resolution K-band spectrum has been obtained of the ultracompact HII region G29.96-0.02 with the Very Large Telescope (VLT). The slit was positioned along the symmetry axis of the cometary shaped nebula. Besides the spectrum of the embedded ionizing O star, the long-slit observation revealed the rich emission-line spectrum produced by the ionized nebula with sub-arcsec spatial resolution. The nebular spectrum includes Br$\gamma$, several helium emission lines and a molecular hydrogen line. A detailed analysis is presented of the variation in strength, velocity and width of the nebular emission lines along the slit. The results are consistent with previous observations, but the much better spatial resolution allows a critical evaluation of models explaining the cometary shape of the nebula. Our observations support neither the wind bow shock model nor the champagne flow model. The measured line ratios of the nebular hydrogen and helium lines are compared to predictions from case B recombination-line theory. The results indicate an electron temperature between 6400 and 7500 K, in good agreement with other determinations and the Galactocentric distance of 4.6 kpc. The He$^{+}$/H$^{+}$ ratio is practically constant over the slit; we argue that He is singly ionized throughout the nebula. We review the various observational constraints on the effective temperature of the ionizing star and show that these are in agreement with its K-band spectral type of O5-O6V.

We propose the late-O, early-B star IRS2b as the ionizing source of the Flame Nebula (NGC 2024). It has been clear that such a hot, massive star must be present in this heavily obscured region, and now it has been identified. New near-infrared photometry shows that IRS2b is the most luminous and hottest star in the young star cluster embedded in the center of NGC 2024. The near-infrared observations (5' x 5') cover ~90 % of the HII region detected in radio continuum radiation, making the probability very low that the ionizing star is not present in the field. A K-band spectrum of IRS2b obtained with ISAAC on the Very Large Telescope indicates that the spectral type of IRS2b is in the range O8V - B2V. Additional arguments indicate that its spectral type is likely closer to O8 than to B2. The corresponding amount of ionizing radiation is consistent with published radio continuum and recombination line observations.

We study the cluster population in a region of 3.2x3.2 kpc^2 in the inner spiral arms of the intergacting galaxy M51, at a distance of about 1 to 3 kpc from the nucleus, based on HST--WFPC2 images taken through five broadband and two narrowband filters. We found 877 cluster candidates and we derived their ages, initial masses and extinctions by comparing their energy distribution with the Starburst99 cluster models. We describe the 3 and 2-dimensional least-square energy fitting method that was used (3DEF, 2DEF). The lack of [OIII] emission in even the youngest clusters with strong H-alpha emission, indicates the absence of the most massive stars and suggests a mass upper limit of about 25 to 30 solar masses. The mass versus age distribution of the clusters shows a drastic decrease in the number of clusters with age, which indicates that cluster disruption is occurring on a timescale of about 10 Myr for low mass clusters. The cluster initial mass function for clusters younger than 10 Myr has an exponent of alpha = 2.0 (+- 0.05) We derived the cluster formation history from clusters with an initial mass larger than 10^4 solar masses. There is no evidence for a peak in the cluster formation rate within a factor two at about 200 to 400 Myr ago, i.e. at the time of the interaction with the companion galaxy NGC 5194.