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.

This work presents a spectroscopic study of 163 Herbig Ae/Be stars. Amongst these, we present new data for 30 objects. Stellar parameters such as temperature, reddening, mass, luminosity and age are homogeneously determined. Mass accretion rates are determined from $\rm H\alpha$ emission line measurements. Our data is complemented with the X-Shooter sample from previous studies and we update results using Gaia DR2 parallaxes giving a total of 78 objects with homogeneously determined stellar parameters and mass accretion rates. In addition, mass accretion rates of an additional 85 HAeBes are determined. We confirm previous findings that the mass accretion rate increases as a function of stellar mass, and the existence of a different slope for lower and higher mass stars respectively. The mass where the slope changes is determined to be $3.98^{+1.37}_{-0.94}\,\rm M_{\odot}$. We discuss this break in the context of different modes of disk accretion for low- and high mass stars. Because of their similarities with T Tauri stars, we identify the accretion mechanism for the late-type Herbig stars with the Magnetospheric Accretion. The possibilities for the earlier-type stars are still open, we suggest the Boundary Layer accretion model may be a viable alternative. Finally, we investigated the mass accretion - age relationship. Even using the superior Gaia based data, it proved hard to select a large enough sub-sample to remove the mass dependency in this relationship. Yet, it would appear that the mass accretion does decline with age as expected from basic theoretical considerations.

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.

We present Halpha linear spectropolarimetry of a large sample of Herbig Ae/Be stars. Together with newly obtained data for 17 objects, the sample contains 56 objects, the largest such sample to date. A change in linear polarization across the Halpha line is detected in 42 (75 %) objects, which confirms the previous finding that the circumstellar environment around these stars on small spatial scales has an asymmetric structure, which is typically identified with a disk. A second outcome of this research is that we confirm that Herbig Ae stars are similar to T Tauri stars in displaying a line polarization effect, while depolarization is more common among Herbig Be stars. This finding had been suggested previously to indicate that Herbig Ae stars form in the same manner than T Tauri stars through magnetospheric accretion. It appears that the transition between these two differing polarization line effects occurs around the B7-B8 spectral type. This would in turn not only suggest that Herbig Ae stars accrete in a similar fashion as lower mass stars, but also that this accretion mechanism switches to a different type of accretion for Herbig Be stars. We report that the magnitude of the line effect caused by electron scattering close to the stars does not exceed 2%. Only a very weak correlation is found between the magnitude of the line effect and the spectral type or the strength of the Halpha line. This indicates that the detection of a line effect only relies on the geometry of the line-forming region and the geometry of the scattering electrons.

We present the results of spectropolarimetric observations of 12 Herbig Ae/Be objects. Our data have the largest spectropolarimetric wavelength coverage, 4560 {\AA} to 9480 {\AA}, published to date. A change in linear polarisation across the H{\alpha} line, is detected in all objects. Such a line effect reveals the fact that stellar photons are scattered off free electrons that are not distributed in a spherically symmetric volume, suggesting the presence of small disks around these accreting objects. Thanks to the large wavelength coverage, we can report that H{\alpha} is the spectral line in the optical wavelength range that is most sensitive to revealing deviations from spherical symmetry, and the one most likely to show a line effect across the polarisation in such cases. Few other spectral lines display changes in polarisation across the line. In addition, H{\alpha} is the only line which shows an effect across its absorption component in some sources. We present a scenario explaining this finding and demonstrate that the detection of the line effect strongly relies on the number of photons scattered into our line of sight. We highlight the special case of R Mon, which is the only object in our sample to show many lines with a polarisation effect, which is much stronger than in all other objects. Given that the object and its nebulosity is spatially resolved, we argue that this is due to scattering of the stellar and emission spectrum off circumstellar dust.

The number of well-studied early-type pre-main-sequence objects is very limited, hampering the study of massive star formation from an observational point of view. Here, we present the results of VLT/FORS2 spectropolarimetric and VLT/X-shooter spectroscopic observations of two recently recognised candidate Herbig Be stars, PDS 27 and PDS 37. Through analysis of spectral lines and photometry, we find that these two objects are hot, 17500 $\pm$ 3500~K, have large radii, 17.0 $\pm$ 4.0 and 25.8 $\pm$ 5.0~${\rm R}_\odot$, and are very massive, 15.3 (+5.4, -4.4) and 21.1 (+11.0, -5.3)~M$_{\odot}$ for PDS 27 and PDS 37, respectively. This suggests that these two objects are very young in their evolution and may become O-type stars. Their youth is supported by their high accretion rates of the order of $10^{-3}$--$10^{-4.5}~{\rm M}_\odot$/yr. A change in linear polarisation across the absorption component of H$\alpha$ is detected in both objects. This change indicates that the circumstellar environment close to the star, at scales of several stellar radii, has a flattened structure, which we identify as an inner accretion disk. Strong variability is seen in both objects in many lines as further indication of an active circumstellar environment.