Context. The fate of a massive star during the latest stages of its evolution is highly dependent on its mass-loss rate/geometry and therefore knowing the geometry of the circumstellar material close to the star and its surroundings is crucial. Aims. We aim to study the nature (i.e. geometry, rates) of mass-loss episodes. In this context, yellow hypergiants are great targets. Methods. We analyse a large set of optical/near-infrared data, in spectroscopic and photometric (X-shooter/VLT), spectropolarimetric (ISIS/WHT), and interferometric GRAVITY-AMBER/VLTI) modes, toward the yellow hypergiant IRAS 17163-3907. We present the first model-independent reconstructed images of IRAS 17163-3907 at these wavelengths at milli-arcsecond scales. Lastly, we apply a 2D radiative transfer model to fit the dereddened photometry and the radial profiles of published VISIR images at 8.59 {\mu}m, 11.85 {\mu}m and 12.81 {\mu}m simultaneously, adopting the revised Gaia distance (DR2). Results. The interferometric observables around 2 {\mu}m show that the Br{\gamma} emission is more extended and asymmetric than the Na i and the continuum emission. In addition to the two known shells surrounding IRAS 17163-3907 we report on the existence of a third hot inner shell with a maximum dynamical age of only 30 yr. Conclusions. The interpretation of the presence of Na i emission at closer distances to the star compared to Br{\gamma} has been a challenge in various studies. We argue that the presence of a pseudophotosphere is not needed, but it is rather an optical depth effect. The three observed distinct mass-loss episodes are characterised by different mass-loss rates and can inform the theories on mass-loss mechanisms, which is a topic still under debate. We discuss these in the context of photospheric pulsations and wind bi-stability mechanisms.

Binarity and multiplicity appear to be a common outcome in star formation. In particular, the binary fraction of massive (OB-type) stars can be very high. In many cases, the further stellar evolution of these stars is affected by binary interactions at some stage during their lifetime. The origin of this high binarity and the binary parameters are poorly understood because observational constraints are scarce, which is predominantly due to a dearth of known young massive binary systems. We aim to identify and describe massive young binary systems in order to fill in the gaps of our knowledge of primordial binarity of massive stars, which is crucial for our understanding of massive star formation. We observed the two massive young stellar objects (MYSOs) PDS 27 and PDS 37 at the highest spatial resolution provided by VLTI/PIONIER in the H-band (1.3 mas). We applied geometrical models to fit the observed squared visibilities and closure phases. In addition, we performed a radial velocity analysis using published VLT/FORS2 spectropolarimetric and VLT/X-shooter spectroscopic observations. Our findings suggest binary companions for both objects at 12 mas (30 au) for PDS 27 and at 22-28 mas (42-54 au) for PDS 37. This means that they are among the closest MYSO binaries resolved to date. Our data spatially resolve PDS 27 and PDS 37 for the first time, revealing two of the closest and most massive ($>$8 M$_\odot$) YSO binary candidates to date. PDS 27 and PDS 37 are rare but great laboratories to quantitatively inform and test the theories on formation of such systems.

It is not well known what drives the chemistry of a protostellar envelope, in particular the role of the stellar mass and the outflows on its chemical enrichment. We study the chemical structure of NGC 1333 IRAS 4A in order to (i) investigate the influence of the outflows on the chemistry, (ii) constrain the age of our object, (iii) compare it with a typical high-mass protostellar envelope. In our analysis we use JCMT line mapping and HIFI pointed spectra. To study the influence of the outflow on the degree of deuteration, we compare JCMT maps of HCO+ and DCO+ with non-LTE (RADEX) models in a region that spatially covers the outflow activity of IRAS 4A. To study the envelope chemistry, we derive empirical molecular abundance profiles for the observed species using the radiative transfer code (RATRAN) and adopting a 1D dust density/temperature profile from the literature. We compare our best-fit observed abundance profiles with the predictions from the time dependent gas grain chemical code (ALCHEMIC). The CO, HCN, HNC and CN abundance require an enhanced UV field which points towards an outflow cavity. The abundances (wrt H2) are 1 to 2 orders of magnitude lower than those observed in the high mass protostellar envelope (AFGL 2591), while they are found to be similar within factors of a few with respect to CO. Differences in UV radiation may be responsible for such chemical differentiation, but temperature differences seem a more plausible explanation. The CH3OH modeled abundance profile points towards an age of > 4x10^4 yrs for IRAS 4A. The spatial distribution of H2D+ differs from that of other deuterated species, indicating an origin from a foreground colder layer (<20 K). The observed abundances can be explained by passive heating towards the high mass protostellar envelope, while the presence of UV cavity channels become more important toward the low mass protostellar envelope.

Protostellar evolution, following the formation of the protostar is becoming reasonably well characterized, but the evolution from a prestellar core to a protostar is not well known, although the first hydrostatic core (FHSC) must be a pivotal step. NGC 1333 IRAS 4C is a potentially very young object, that we directly compare with the nearby Class 0 IRAS 4A and IRAS 4B. Observational constraints are provided by spectral imaging from the JCMT Spectral Legacy Survey (330-373 GHz) and continuum and line observations from CARMA. We present integrated intensity and velocity maps of several species, including CO, H2CO and CH3OH. The velocity of an observed outflow, the degree of CO depletion, the deuteration of DCO+/HCO+ and gas kinetic temperatures are observational signatures that we present. We report differences between the three sources in four aspects: a) the kinetic temperature is much lower towards IRAS 4C, b) the line profiles of the detected species show strong outflow activity towards IRAS 4A and IRAS 4B but not towards IRAS 4C, c) the HCN/HNC is < 1 towards IRAS 4C, which confirms the cold nature of the source, d) the degree of CO depletion and the deuteration are the lowest towards the warmest of the sources, IRAS 4B. IRAS 4C seems to be in a different evolutionary state than the IRAS 4A and IRAS 4B sources. We can probably exclude the FHSC stage due to the relatively low Lsmm/Lbol (~6\%) and we investigate the earliest accretion phase of Class 0 stage and the transition between Class 0 to Class I. Our results do not show a consistent scenario for either case, with the major issue being the absence of outflow activity and the cold nature of IRAS 4C. The number of FHSC candidates in Perseus is ~10 times higher than current models predict, which suggests that the lifespan of these objects is greater than 1000 years, possibly due to an accretion rate lower than 4x10^-5 msun/yr.

The gas-solid budget of carbon in protoplanetary disks is related to the composition of the cores and atmospheres of the planets forming in them. The key gas-phase carbon carriers CO, C$^{0}$ and C$^{+}$ can now be observed in disks. The gas-phase carbon abundance in disks has not yet been well characterized, we aim to obtain new constraints on the [C]/[H] ratio in a sample of disks, and to get an overview of the strength of [CI] and warm CO emission. We carried out a survey of the CO$\,6$--$5$ and [CI]$\,1$--$0$ and $2$--$1$ lines towards $37$ disks with APEX, and supplemented it with [CII] data from the literature. The data are interpreted using a grid of models produced with the DALI code. We also investigate how well the gas-phase carbon abundance can be determined in light of parameter uncertainties. The CO$\,6$--$5$ line is detected in $13$ out of $33$ sources, the [CI]$\,1$--$0$ in $6$ out of $12$, and the [CI]$\,2$--$1$ in $1$ out of $33$. With deep integrations, the first unambiguous detections of [CI]~$1$--$0$ in disks are obtained, in TW~Hya and HD~100546. Gas-phase carbon abundance reductions of a factor $5$--$10$ or more can be identified robustly based on CO and [CI] detections. The atomic carbon detection in TW~Hya confirms a factor $100$ reduction of [C]/[H]$_{\rm gas}$ in that disk, while the data are consistent with an ISM-like carbon abundance for HD~100546. In addition, BP~Tau, T~Cha, HD~139614, HD~141569, and HD~100453 are either carbon-depleted or gas-poor disks. The low [CI]~$2$--$1$ detection rates in the survey mostly reflect insufficient sensitivity to detect T~Tauri disks. The Herbig~Ae/Be disks with CO and [CII] upper limits below the models are debris disk like systems. A roughly order of magnitude increase in sensitivity compared to our survey is required to obtain useful constraints on the gas-phase [C]/[H] ratio in most of the targeted systems.

We present a study of the properties of the O-type, massive eclipsing binary 2MASS J13130841-6239275 located in the outskirts of the Danks 2 cluster in the G305 star-forming complex, using near-infrared spectroscopy from VLT/ISAAC. We derive the masses and radii to be 24.5$\pm$0.9 M$_{\odot}$ and 9.2$\pm$0.1 R$_{\odot}$ for the primary and 21.7$\pm$0.8 M$_{\odot}$ and 8.7$\pm$0.1 R$_{\odot}$ for the secondary component. In addition, we evaluate the sensitivity of our parameters to the choice of the spectral features used to determine the radial velocities. Both components appear to be main-sequence O6.5$-$O7 type stars at an age of $\sim$5 Myr, which is in agreement with the age of the cluster. A high visual extinction of A$_{5495}$=11.9$\pm$0.1 mag is reported, which is likely attributed to the cold molecular gas contaminating the north-east region of the cluster. By fitting the spectral energy distribution of the system to the available $BVI_{c}JHK_{s}$ photometry, we determine a distance to the system of 3.52$\pm$0.08 kpc with a precision of 2$\%$, which is the most well-determined distance to the Danks 2 cluster and the host complex reported in the literature.

Intermediate mass protostarsprovide a bridge between theories of low- and high-mass star formation. Emerging molecular outflows can be used to determine the influence of fragmentation and multiplicity on protostellar evolution through the correlation of outflow forces of intermediate mass protostars with the luminosity. The aim of this paper is to derive outflow forces from outflows of six intermediate mass protostellar regions and validate the apparent correlation between total luminosity and outflow force seen in earlier work, as well as remove uncertainties caused by different methodology. By comparing CO 6--5 observations obtained with APEX with non-LTE radiative transfer model predictions, optical depths, temperatures, densities of the gas of the molecular outflows are derived. Outflow forces, dynamical timescales and kinetic luminosities are subsequently calculated. Outflow parameters, including the forces, were derived for all sources. Temperatures in excess of 50 K were found for all flows, in line with recent low-mass results. However, comparison with other studies could not corroborate conclusions from earlier work on intermediate mass protostars which hypothesized that fragmentation enhances outflow forces in clustered intermediate mass star formation. Any enhancement in comparison with the classical relation between outflow force and luminosity can be attributed the use of a higher excitation line and improvement in methods; They are in line with results from low-mass protostars using similar techniques. The role of fragmentation on outflows is an important ingredient to understand clustered star formation and the link between low and high-mass star formation. However, detailed information on spatial scales of a few 100 AU, covering all individual members is needed to make the necessary progress.

We try to understand the gas heating and cooling in the S 140 star forming region by spatially and spectrally resolving the distribution of the main cooling lines with GREAT/SOFIA. We mapped the fine structure lines of [OI] (63 {\mu}m) and [CII] (158 {\mu}m) and the rotational transitions of CO 13-12 and 16-15 with GREAT/SOFIA and analyzed the spatial and velocity structure to assign the emission to individual heating sources. We measure the optical depth of the [CII] line and perform radiative transfer computations for all observed transitions. By comparing the line intensities with the far-infrared continuum we can assess the total cooling budget and measure the gas heating efficiency. The main emission of fine structure lines in S 140 stems from a 8.3'' region close to the infrared source IRS 2 that is not prominent at any other wavelength. It can be explained by a photon-dominated region (PDR) structure around the embedded cluster if we assume that the [OI] line intensity is reduced by a factor seven due to self-absorption. The external cloud interface forms a second PDR at an inclination of 80-85 degrees illuminated by an UV field of 60 times the standard interstellar radiation field. The main radiation source in the cloud, IRS 1, is not prominent at all in the fine structure lines. We measure line-to-continuum cooling ratios below 10^(-4), i.e. values lower than in any other Galactic source, rather matching the far-IR line deficit seen in ULIRGs. In particular the low intensity of the [CII] line can only be modeled by an extreme excitation gradient in the gas around IRS 1. We found no explanation why IRS 1 shows no associated fine-structure line peak, while IRS 2 does. The inner part of S 140 mimics the far-IR line deficit in ULIRGs thereby providing a template that may lead to a future model.

In dense parts of interstellar clouds (> 10^5 cm^-3), dust & gas are expected to be in thermal equilibrium, being coupled via collisions. However, previous studies have shown that the temperatures of the dust & gas may remain decoupled even at higher densities. We study in detail the temperatures of dust & gas in the photon-dominated region S 140, especially around the deeply embedded infrared sources IRS 1-3 and at the ionization front. We derive the dust temperature and column density by combining Herschel PACS continuum observations with SOFIA observations at 37 $\mu$m and SCUBA at 450 $\mu$m. We model these observations using greybody fits and the DUSTY radiative transfer code. For the gas part we use RADEX to model the CO 1-0, CO 2-1, 13CO 1-0 and C18O 1-0 emission lines mapped with the IRAM-30m over a 4' field. Around IRS 1-3, we use HIFI observations of single-points and cuts in CO 9-8, 13CO 10-9 and C18O 9-8 to constrain the amount of warm gas, using the best fitting dust model derived with DUSTY as input to the non-local radiative transfer model RATRAN. We find that the gas temperature around the infrared sources varies between 35 and 55K and that the gas is systematically warmer than the dust by ~5-15K despite the high gas density. In addition we observe an increase of the gas temperature from 30-35K in the surrounding up to 40-45K towards the ionization front, most likely due to the UV radiation from the external star. Furthermore, detailed models of the temperature structure close to IRS 1 show that the gas is warmer and/or denser than what we model. Finally, modelling of the dust emission from the sub-mm peak SMM 1 constrains its luminosity to a few ~10^2 Lo. We conclude that the gas heating in the S 140 region is very efficient even at high densities, most likely due to the deep UV penetration from the embedded sources in a clumpy medium and/or oblique shocks.

Only a small number of high mass stars (> 30 Mo) have fundamental parameters (i.e. masses and radii) measured with high enough accuracy from eclipsing binaries to constrain formation and evolutionary models of massive stars. This work aims to increase this limited sample, by studying the 4 massive eclipsing binary candidates discovered by Bonanos in the young massive cluster Westerlund 1. We present new follow-up echelle spectroscopy of these binaries and models of their light and radial velocity curves. We obtain fundamental parameters (i.e. masses, radii) for the 8 component stars, finding masses that span a range of 10-40 Mo, and contributing accurate fundamental parameters for 1 additional very massive star, the 33 Mo component of W13. WR77o is found to have a ~40 Mo companion, which provides a second dynamical constraint on the mass of the progenitor of the magnetar known in the cluster. We also use W13 to estimate the first, direct, eclipsing binary distance to Westerlund 1 and therefore the magnetar, and find it to be at 3.7 +/- 0.6 kpc. Our results confirm previous evidence for a high mass for the progenitor of the magnetar. In addition, the availability of eclipsing binaries with accurate parameters opens the way for direct, independent, high precision eclipsing binary distance measurements to Westerlund 1.

We present fundamental parameters of 4 massive eclipsing binaries in the young massive cluster Westerlund 1. The goal is to measure accurate masses and radii of their component stars, which provide much needed constraints for evolutionary models of massive stars. Accurate parameters can further be used to determine a dynamical lower limit for the magnetar progenitor and to obtain an independent distance to the cluster. Our results confirm and extend the evidence for a high mass for the progenitor of the magnetar.

Westerlund 1 is one of the most massive young clusters known in the Local Group, with an age of 3-5 Myr. It contains an assortment of rare evolved massive stars, such as blue, yellow and red supergiants, Wolf-Rayet stars, a luminous blue variable, and a magnetar, as well as 4 massive eclipsing binary systems (Wddeb, Wd13, Wd36, WR77o, see Bonanos 2007). The eclipsing binaries present a rare opportunity to constrain evolutionary models of massive stars, the distance to the cluster and furthermore, to determine a dynamical lower limit for the mass of a magnetar progenitor. Wddeb, being a detached system, is of great interest as it allows determination of the masses of 2 of the most massive unevolved stars in the cluster. We have analyzed spectra of all 4 eclipsing binaries, taken in 2007-2008 with the 6.5 meter Magellan telescope at Las Campanas Observatory, Chile, and present fundamental parameters (masses, radii) for their component stars.