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.

The structure of cometary dust is a tracer of growth processes in the formation of planetesimals. Instrumentation on board the Rosetta mission to comet 67P/Churyumov- Gerasimenko captured dust particles and analysed them in situ. However, these deposits are a product of a collision within the instrument. We conducted laboratory experiments with cometary dust analogues, simulating the collection process by Rosetta instruments (specifically COSIMA, MIDAS). In Paper I we reported that velocity is a key driver in determining the appearance of deposits. Here in Paper II we use materials with different monomer sizes, and study the effect of tensile strength on the appearance of deposits. We find that mass transfer efficiency increases from $\sim$1 up to $\sim$10% with increasing monomer diameter from 0.3 $\mu$m to 1.5 $\mu$m (i.e. tensile strength decreasing from $\sim$12 to $\sim$3 kPa), and velocities increasing from 0.5 to 6 m/s. Also, the relative abundance of small fragments after impact is higher for material with higher tensile strength. The degeneracy between the effects of velocity and material strength may be lifted by performing a closer study of the deposits. This experimental method makes it possible to estimate the mass transfer efficiency in the COSIMA instrument. Extrapolating these results implies that more than half of the dust collected during the Rosetta mission has not been imaged. We analysed two COSIMA targets containing deposits from single collisions. The collision that occurred closest to perihelion passage led to more small fragments on the target.

Cometary dust provides a unique window on dust growth mechanisms during the onset of planet formation. Measurements by the Rosetta spacecraft show that the dust in the coma of comet 67P/Churyumov-Gerasimenko has a granular structure at size scales from sub-um up to several hundreds of um, indicating hierarchical growth took place across these size scales. However, these dust particles may have been modified during their collection by the spacecraft instruments. Here we present the results of laboratory experiments that simulate the impact of dust on the collection surfaces of COSIMA and MIDAS, instruments onboard the Rosetta spacecraft. We map the size and structure of the footprints left by the dust particles as a function of their initial size (up to several hundred um) and velocity (up to 6 m/s). We find that in most collisions, only part of the dust particle is left on the target; velocity is the main driver of the appearance of these deposits. A boundary between sticking/bouncing and fragmentation as an outcome of the particle-target collision is found at v ~ 2 m/s. For velocities below this value, particles either stick and leave a single deposit on the target plate, or bounce, leaving a shallow footprint of monomers. At velocities > 2 m/s and sizes > 80 um, particles fragment upon collision, transferring up to 50 per cent of their mass in a rubble-pile-like deposit on the target plate. The amount of mass transferred increases with the impact velocity. The morphologies of the deposits are qualitatively similar to those found by the COSIMA instrument.

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.

HD 50138 is a B[e] star surrounded by a large amount of circumstellar gas and dust. Its spectrum shows characteristics which may indicate either a pre- or a post-main-sequence system. Mapping the kinematics of the gas in the inner few au of the system contributes to a better understanding of its physical nature. We present the first high spatial and spectral resolution interferometric observations of the Br-gamma line of HD~50138, obtained with VLTI/AMBER. The line emission originates from a region more compact (up to 3 au) than the continuum-emitting region. Blue- and red-shifted emission originates from the two different hemispheres of an elongated structure perpendicular to the polarization angle. The velocity of the emitting medium decreases radially. An overall offset along the NW direction between the line- and continuum-emitting regions is observed. We compare the data with a geometric model of a thin Keplerian disk and a spherical halo on top of a Gaussian continuum. Most of the data are well reproduced by this model, except for the variability, the global offset and the visibility at the systemic velocity. The evolutionary state of the system is discussed; most diagnostics are ambiguous and may point either to a post-main-sequence or a pre-main-sequence nature.

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

At low redshift, a handful of gamma-ray bursts (GRBs) have been discovered with peak luminosities ($L_{\rm iso} < 10^{48.5}~\rm{erg\,s}^{-1}$) substantially lower than the average of the more distant ones ($L_{\rm iso} > 10^{49.5}~\rm{erg\,s}^{-1}$). The properties of several low-luminosity (low-$L$) GRBs indicate that they can be due to shock break-out, as opposed to the emission from ultrarelativistic jets. Owing to this, it is highly debated how both populations are connected, and whether there is a continuum between them. The burst at redshift $z=0.283$ from 2012 April 22 is one of the very few examples of intermediate-$L$ GRBs with a $\gamma$-ray luminosity of $L\sim10^{48.9}~\rm{erg\,s}^{-1}$ that have been detected up to now. Together with the robust detection of its accompanying supernova SN 2012bz, it has the potential to answer important questions on the origin of low- and high-$L$ GRBs and the GRB-SN connection. We carried out a spectroscopy campaign using medium- and low-resolution spectrographs at 6--10-m class telescopes, covering the time span of 37.3 days, and a multi-wavelength imaging campaign from radio to X-ray energies over a duration of $\sim270$ days. Furthermore, we used a tuneable filter centred at H$\alpha$ to map star formation in the host galaxy and the surrounding galaxies. We used these data to extract and model the properties of different radiation components and incorporate spectral-energy-distribution fitting techniques to extract the properties of the host galaxy. Modelling the light curve and spectral energy distribution from the radio to the X-rays revealed the blast-wave to expand with an initial Lorentz factor of $\Gamma_0\sim60$, low for a high-$L$ GRB, and that the afterglow had an exceptional low peak luminosity-density of $\lesssim2\times10^{30}~\rm{erg\,s}^{-1}\,\rm{Hz}^{-1}$ in the sub-mm. [Abridged]

Herbig Ae/Be stars are intermediate-mass pre-main sequence stars surrounded by circumstellar dust disks. Some are observed to produce jets, whose appearance as a sequence of shock fronts (knots) suggests a past episodic outflow variability. This "jet fossil record" can be used to reconstruct the outflow history. We present the first optical to near-infrared (NIR) VLT/X-shooter spectra of the jet from the Herbig Ae star HD 163296. We determine physical conditions in the knots, as well as their kinematic "launch epochs". Knots are formed simultaneously on either side of the disk, with a regular interval of ~16 yr. The velocity dispersion versus jet velocity and the energy input are comparable in both lobes. However, the mass loss rate, velocity, and shock conditions are asymmetric. We find Mjet/Macc ~ 0.01-0.1, consistent with magneto-centrifugal jet launching models. No evidence for dust is found in the high-velocity jet, suggesting it is launched within the sublimation radius (<0.5 au). The jet inclination measured from proper motions and radial velocities confirms it is perpendicular to the disk. A tentative relation is found between the structure of the jet and the photometric variability of the source. Episodes of NIR brightening were previously detected and attributed to a dusty disk wind. We report for the first time significant optical fadings lasting from a few days up to a year, coinciding with the NIR brightenings. These are likely caused by dust lifted high above the disk plane; this supports the disk wind scenario. The disk wind is launched at a larger radius than the high-velocity atomic jet, although their outflow variability may have a common origin. No significant relation between outflow and accretion variability could be established. Our findings confirm that this source undergoes periodic ejection events, which may be coupled with dust ejections above the disk plane.

We present the results of a quantitative spectroscopic analysis of the oxygen-sequence Wolf- Rayet star DR1 in the low-metallicity galaxy IC 1613. Our models suggest that the strong oxygen emission lines are the result of the high temperature of this WO3 star and do not necessarily reflect a more advanced evolutionary stage than WC stars.

Context. Oxygen sequence Wolf-Rayet (WO) stars are thought to represent the final evolutionary stage of the most massive stars. The characteristic strong O vi emission possibly originates from an enhanced oxygen abundance in the stellar wind. Alternatively, the O vi emission can be caused by the high temperature of these stars, in which case the WO stars are the high-temperature extension of the more common carbon sequence Wolf-Rayet (WC) stars. Aims. By constraining the physical properties and evolutionary status of DR1, a WO star in the low-metallicity Local Group dwarf galaxy IC 1613 and one of only two objects of its class known in a SMC-like metallicity environment, we aim to investigate the nature of WO stars and their evolutionary connection with WC stars. Methods. We use the non-LTE atmosphere code cmfgen to model the observed spectrum of DR1 and to derive its stellar and wind parameters. We compare our values with other studies of WC and WO stars, as well as with the predictions of evolutionary models. We also model the surrounding nebula using the photo-ionization code cloudy. Results. The oxygen and carbon abundances that we obtain are comparable to values found for WC stars. The temperature and luminosity are, however, higher than those of WC stars. DR1 is embedded in the hottest known H ii region in the Local Group. The nebular properties can be consistently reproduced by cloudy models adopting DR1 as central ionizing source. Conclusions. Comparison of the abundances and temperature of DR1 with core helium-burning models show that DR1 is currently well into the second half of helium burning. If the properties of DR1 are representative for the WO class, it would imply that WO stars are the high-temperature and high-luminosity extension of the WC stars, and do not necessarily represent a later evolutionary stage.

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.

AA Tau has been monitored for more than 20 years since 1987, exhibiting a nearly constant brightness level of V=12.5 mag. We report here that in 2011 it suddenly faded away, becoming 2 magnitudes fainter in the V-band, and has remained in this deep state since then. We report new optical and near-IR photometry and spectroscopy obtained during the fading event. The system appears much redder and fainter than it was in the bright state. Also, the 8.2d photometric period continuously observed for more than 20 years is not seen during most of the deep state. The analysis of the system's brightness and colors suggests that the visual extinction on the line of sight has increased by about 3-4 magnitudes in the deep state. At optical wavelengths, the system appears to be dominated by scattered light, probably originating from the upper surface layers of a highly inclined circumstellar disk. The profiles of the Balmer lines have significantly changed as well, with the disappearance of a central absorption component regularly observed in the bright state. We ascribe this change to the scattering of the system's spectrum by circumstellar dust. Noticeably, the mass accretion rate in the inner disk and onto the central star has not changed as the system faded away. We conclude that the deepening of the AA Tau system is due to a sudden increase of circumstellar dust extinction on the line of sight without concomitant change in the accretion rate. We suggest that the enhanced obscuration may be produced by a non-axisymmetric overdense region in the disk, located at a distance of 7.7 AU or more, that was recently brought on the line of sight by its keplerian motion around the central star.

R144 is a WN6h star in the 30 Doradus region. It is suspected to be a binary because of its high luminosity and its strong X-ray flux, but no periodicity could be established so far. Here, we present new Xshooter multi-epoch spectroscopy of R144 obtained at the ESO Very Large Telescope (VLT). We detect variability in position and/or shape of all the spectral lines. We measure radial velocity variations with an amplitude larger than 250 km/s in NIV and NV lines. Furthermore, the NIII and NV line Doppler shifts are anti-correlated and the NIV lines show a double-peaked profile on six of our seven epochs. We thus conclude that R144 is a double-lined spectroscopic binary. Possible orbital periods range from 2 to 6 months, although a period up to one year is allowed if the orbit is highly eccentric. We estimate the spectral types of the components to be WN5-6h and WN6-7h, respectively. The high luminosity of the system (log Lbol/Lsun ~ 6.8) suggests a present-day total mass content in the range of about 200 to 300 Msun, depending on the evolutionary stage of the components. This makes R144 the most massive binary identified so far, with a total mass content at birth possibly as large as 400 Msun. We briefly discuss the presence of such a massive object 60 pc away from the R136 cluster core in the context of star formation and stellar dynamics.

Jets around low- and intermediate-mass young stellar objects (YSOs) contain a fossil record of the recent accretion and outflow activity of their parent star-forming systems. We aim to understand whether the accretion/ejection process is similar across the entire stellar mass range of the parent YSOs. To this end we have obtained VLT/X-shooter spectra of HH 1042 and HH 1043, two newly discovered jets in the massive star-forming region RCW 36. HH 1042 is associated with the intermediate-mass YSO 08576nr292. Over 90 emission lines are detected in the spectra. High-velocity (up to 220 km/s) blue- and redshifted emission from a bipolar flow is observed in typical shock tracers. Low-velocity emission from the background cloud is detected in nebular tracers, including lines from high ionization species. We applied combined optical and infrared spectral diagnostic tools in order to derive the physical conditions (density, temperature, and ionization) in the jets. The measured mass outflow rates are Mjet ~ 10^-7 Msun/yr. We measure a high accretion rate for HH 1042 (Macc ~ 10^-6 Msun/yr) and Mjet/Macc ~ 0.1, comparable to low-mass sources and consistent with models for magneto-centrifugal jet launching. The knotted structure and velocity spread in both jets are interpreted as fossil signatures of a variable outflow rate. The mean velocities in both lobes of the jets are comparable, but the variations in Mjet and velocity in the two lobes are not symmetric, suggesting that the launching mechanism on either side of the accretion disk is not synchronized. For HH 1042, we have constructed an interpretative physical model with a stochastic or periodic outflow rate and a description of a ballistic flow as its constituents. The knotted structure and velocity spread can be reproduced qualitatively with the model, indicating that the outflow velocity varies on timescales on the order of 100 yr.

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 present the optical to near-infrared spectrum of MAXI J1659-152, during the onset of its 2010 X-ray outburst. The spectrum was obtained with X-shooter on the ESO - Very Large Telescope (VLT) early in the outburst simultaneous with high quality observations at both shorter and longer wavelengths. At the time of the observations, the source was in the low-hard state. The X-shooter spectrum includes many broad (~2000 km/s), double-peaked emission profiles of H, HeI, HeII, characteristic signatures of a low-mass X-ray binary during outburst. We detect no spectral signatures of the low-mass companion star. The strength of the diffuse interstellar bands results in a lower limit to the total interstellar extinction of Av ~ 0.4 mag. Using the neutral hydrogen column density obtained from the X-ray spectrum we estimate Av ~1 mag. The radial-velocity structure of the interstellar NaI D and CaII H & K lines results in a lower limit to the distance of ~ 4 +/- 1 kpc, consistent with previous estimates. With this distance and Av, the dereddened spectral energy distribution represents a flat disk spectrum. The two subsequent 10 minute X-shooter spectra show significant variability in the red wing of the emission-line profiles, indicating a global change in the density structure of the disk, though on a timescale much shorter than the typical viscous timescale of the disk.

The optical to near-infrared (300 - 2500 nm) spectrum of the candidate massive Young Stellar Object (YSO) B275, embedded in the star forming region M17, has been observed with X-shooter on the ESO Very Large Telescope. The spectrum includes both photospheric absorption lines and emission features (H and Ca triplet emission lines, 1st and 2nd overtone CO bandhead emission), as well as an infrared excess indicating the presence of a (flaring) circumstellar disk. The strongest emission lines are double-peaked with a peak separation ranging between 70 and 105 km s^-1 and provide information on the physical structure of the disk. The underlying photospheric spectrum is classified as B6-B7, significantly cooler than a previous estimate based on modeling of the spectral energy distribution. This discrepancy is solved by allowing for a larger stellar radius (i.e. a bloated star), and thus positioning the star above the main sequence. This constitutes the first firm spectral classification of an early-B pre-main-sequence (PMS) star. We discuss the position of B275 in the Hertzsprung-Russell diagram in terms of PMS evolution. Although the position is consistent with PMS tracks of heavily accreting protostars (Macc > 10^-5 Msun yr^-1), the fact that the photosphere of the object is detectable suggests that the current mass-accretion rate is not very high.

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.