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 Hα 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 M . 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.

We present Hα 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 Hα 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 Hα 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.

Context. The fate of a massive star during the latest stages of its evolution is highly dependent on its mass-loss rate and geometry and therefore knowing the geometry of the circumstellar material close to the star and its surroundings is crucial. Aims. We aim to provide insight into the nature (i.e. geometry, rates) of mass-loss episodes, and in particular, the connection between the observed asymmetries due to the mass lost in a fast wind or during a previous, prodigious mass-losing phase. In this context, yellow hypergiants offer a good opportunity to study mass-loss events. Methods. We analysed a large set of optical and near-infrared data in spectroscopic and photometric, spectropolarimetric, and interferometric (GRAVITY/VLTI) modes, towards the yellow hypergiant IRAS 17163−3907. We used X-shooter optical observations to determine the spectral type of this yellow hypergiant and we present the first model-independent, reconstructed images of IRAS 17163−3907 at these wavelengths tracing milli-arcsecond scales. Lastly, we applied a 2D radiative transfer model to fit the dereddened photometry and the radial profiles of published diffraction-limited VISIR images at 8.59 μm, 11.85 μm, and 12.81 μm simultaneously, adopting a revised distance determination using Gaia Data Release 2 measurements. Results. We constrain the spectral type of IRAS 17163−3907 to be slightly earlier than A6Ia (Teff ∼ 8500 K). The interferometric observables around the 2 μm window towards IRAS 17163−3907 show that the Brγ emission appears to be more extended and asymmetric than the Na I and the continuum emission. Interestingly, the spectrum of IRAS 17163−3907 around 2 μm shows Mg II emission that is not previously seen in other objects of its class. In addition, Brγ shows variability in a time interval of four months that is not seen towards Na I. Lastly, 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 2 μm continuum originates directly from the star and not from hot dust surrounding the stellar object. The observed spectroscopic variability of Brγ could be a result of variability in the mass-loss rate. The interpretation of the presence of Na I emission at closer distances to the star compared to Brγ has been a challenge in various studies. To address this, we examine several scenarios. We argue that the presence of a pseudo-photosphere, which was traditionally considered to be the prominent explanation, is not needed and that it is rather an optical depth effect. The three observed distinct mass-loss episodes are characterised by different mass-loss rates and can inform theories of mass-loss mechanisms, which is a topic still under debate both in theory and observations. We discuss these in the context of photospheric pulsations and wind bi-stability mechanisms.

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 ± 3500 K, have large radii, 17.0 ± 4.0 and 25.8 ± 5.0 R ⊙ , and are very massive, 15.3 (+5.4, -4.4) and 21.1 (+11.0, -5.3) M ⊙ 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 M ⊙ /yr. A change in linear polarisation across the absorption component of Hα 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 disc. Strong variability is seen in both objects in many lines as further indication of an active circumstellar environment.

Context. 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. Aims. 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. Methods. 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. Results. 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. Conclusions. Our data spatially resolve PDS 27 and PDS 37 for the first time, revealing two of the closest and most massive (>8 M ⊙ ) 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 the results of spectropolarimetric observations of 12 Herbig Ae/Be objects. Our data have the largest spectropolarimetric wavelength coverage, 4560Å to 9480Å, published to date. A change in linear polarisation across the Hα 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α 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α 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 fate of a massive star during the latest stages of its evolution is highly dependent on its mass-loss history and geometry, with the yellow hypergiants being key objects. We present near-IR interferometric observations of the famous yellow hypergiant IRC +10420 and blue spectra taken between 1994-2019. Our 2.2 μm GRAVITY/VLTI observations attain a spatial resolution of ∼5 stellar radii and spatially resolve the hot emission in the K-band tracing the gas via Na i doublet emission and the Brγ emission. Our geometric modelling reveals a compact neutral zone (Na i) which is slightly larger than the continuum but within an extended Brγ emitting region. Our study confirms an hour-glass geometry of the wind, but we find no signature of a companion at 7-800 au separations at the contrast limit of our observations (3.7 mag at 3σ) to explain this geometry. We report an evolution of the ejecta over 7 years, and constrain the opening angle of the hour-glass to be <10○. Lastly, we present the first blue optical spectra of IRC +10420 since 1994. The multi-epoch data indicate that the spectral type, and thus temperature, of the object has essentially remained constant during the intervening years. Therefore, the observed increase in temperature of 2000 K in less than two decades prior to 1994 is now halted. This suggests that this yellow hypergiant has ‘hit’ the White Wall in the HR-diagram preventing it from evolving blue-wards, and will likely undergo a major mass-loss event in the near future.