To date, there is no explanation as to why disc-tracing CO first overtone (or `bandhead') emission is not a ubiquitous feature in low- to medium-resolution spectra of massive young stellar objects, but instead is only detected toward approximately 25 per cent of their spectra. In this paper, we investigate the hypothesis that only certain mass accretion rates result in detectable bandhead emission in the near infrared spectra of MYSOs. Using an analytic disc model combined with an LTE model of the CO emission, we find that high accretion rates ($\gtrsim 10^{-4}\,{\rm M}_{\odot}{\mathrm{yr}}^{-1}$) result in large dust sublimation radii, a larger contribution to the $K$-band continuum from hot dust at the dust sublimation radius, and therefore correspondingly lower CO emission with respect to the continuum. On the other hand, low accretion rates ($\lesssim10^{-6}\,{\rm M}_{\odot}{\mathrm{yr}}^{-1}$) result in smaller dust sublimation radii, a correspondingly smaller emitting area of CO, and thus also lower CO emission with respect to the continuum. In general, moderate accretion rates produce the most prominent, and therefore detectable, CO first overtone emission. We compare our findings to a recent near-infrared spectroscopic survey of MYSOs, finding results consistent with our hypothesis. We conclude that the detection rate of CO bandhead emission in the spectra of MYSOs could be the result of MYSOs exhibiting a range of mass accretion rates, perhaps due to the variable accretion suggested by recent multi-epoch observations of these objects.

We used Herschel Hi-GAL survey data to determine whether massive young stellar objects (MYSOs) are resolved at 70$\mu$m and to study their envelope density distribution. Our analysis of three relatively isolated sources in the l=30{\deg} and l=59{\deg} Galactic fields show that the objects are partially resolved at 70$\mu$m. The Herschel Hi-GAL survey data have a high scan velocity which makes unresolved and partially resolved sources appear elongated in the 70$\mu$m images. We analysed the two scan directions separately and examine the intensity profile perpendicular to the scan direction. Spherically symmetric radiative transfer models with a power law density distribution were used to study the circumstellar matter distribution. Single dish sub-mm data were also included to study how different spatial information affects the fitted density distribution. The density distribution which best fits both the 70$\mu$m intensity profile and SED has an average index of ~0.5. This index is shallower than expected and is probably due to the dust emission from bipolar outflow cavity walls not accounted for in the spherical models. We conclude that 2D axisymmetric models and Herschel images at low scan speeds are needed to better constrain the matter distribution around MYSOs.

Herbig Ae/Be stars lie in the mass range between low and high mass young stars, and therefore offer a unique opportunity to observe any changes in the formation processes that may occur across this boundary. This paper presents medium resolution VLT/X-Shooter spectra of six Herbig Ae/Be stars, drawn from a sample of 91 targets, and high resolution VLT/CRIRES spectra of five Herbig Ae/Be stars, chosen based on the presence of CO first overtone bandhead emission in their spectra. The X-Shooter survey reveals a low detection rate of CO first overtone emission (7 per cent), consisting of objects mainly of spectral type B. A positive correlation is found between the strength of the CO v=2-0 and Br {\gamma} emission lines, despite their intrinsic linewidths suggesting a separate kinematic origin. The high resolution CRIRES spectra are modelled, and are well fitted under the assumption that the emission originates from small scale Keplerian discs, interior to the dust sublimation radius, but outside the co-rotation radius of the central stars. In addition, our findings are in very good agreement for the one object where spatially resolved near-infrared interferometric studies have also been performed. These results suggest that the Herbig Ae/Be stars in question are in the process of gaining mass via disc accretion, and that modelling of high spectral resolution spectra is able to provide a reliable probe into the process of stellar accretion in young stars of intermediate to high masses.

Context. HD 85567 is an enigmatic object exhibiting the B[e] phenomenon, i.e. an infrared excess and forbidden emission lines in the optical. The object's evolutionary status is uncertain and there are conflicting claims that it is either a young stellar object or an evolved, interacting binary. Aims. To elucidate the reason for the B[e] behaviour of HD 85567, we have observed it with the VLTI and AMBER. Methods. Our observations were conducted in the K-band with moderate spectral resolution (R~1500, i.e. 200 km/s). The spectrum of HD 85567 exhibits Br gamma and CO overtone bandhead emission. The interferometric data obtained consist of spectrally dispersed visibilities, closure phases and differential phases across these spectral features and the K-band continuum. Results. The closure phase observations do not reveal evidence of asymmetry. The apparent size of HD 85567 in the K-band was determined by fitting the visibilities with a ring model. The best fitting radius, 0.8 +/- 0.3 AU, is relatively small making HD 85567 undersized in comparison to the size-luminosity relationship based on YSOs of low and intermediate luminosity. This has previously been found to be the case for luminous YSOs, and it has been proposed that this is due to the presence of an optically thick gaseous disc. We demonstrate that the differential phase observations over the CO bandhead emission are indeed consistent with the presence of a compact (~1 AU) gaseous disc interior to the dust sublimation radius. Conclusions. The observations reveal no sign of binarity. However, the data do indicate the presence of a gaseous disc interior to the dust sublimation radius. We conclude that the data are consistent with the hypothesis that HD 85567 is a YSO with an optically thick gaseous disc within a larger dust disc that is being photo-evaporated from the outer edge.

We present a detailed spectroscopic investigation of the young eruptive star PV Cep, to improve our understanding of its nature and characterise its circumstellar environment after its last outburst in 2004. The analysis of our medium-resolution spectroscopy in the near-IR (0.9-2.35 um), collected in 2012 at the Large Binocular Telescope with the IR spectrograph LUCIFER, allows us to infer the main stellar parameters (visual extinction, accretion luminosity, mass accretion and ejection rates), and model the inner disc, jet, and wind. The NIR spectrum displays several strong emission lines associated with accretion/ejection activity and circumstellar environment. Our analysis shows that the brightness of PV Cep is fading, as well as the mass accretion rate (2x10^-7 Msun/yr^-1 in 2012 vs ~5x10^-6 Msun/yr^-1 in 2004), which is more than one order of magnitude lower than in the outburst phase. Among the several emission lines, only the [FeII] intensity increased after the outburst. The observed [FeII] emission delineates blue- and red-shifted lobes, both with high- and low-velocity components, which trace an asymmetric jet and wind, respectively. The observed emission in the jet has a dynamical age of ~8 years, indicating that it was produced during the last outburst. The mass ejection rate in both lobes is 1.5x10^-7 Msun/yr^-1, approximately matching the high accretion rate observed during and immediately after the outburst . The observed jet/outflow asymmetries are consistent with an inhomogeneous medium. Our modelling of the CO emission hints at a small-scale gaseous disc ring, extending from ~0.2-0.4 AU to ~3 AU from the source, with an inner temperature of ~3000 K. Our HI lines modelling indicates that most of the observed emission comes from an expanding disc wind at Te=10000 K. The line profiles are strongly affected by scattering, disc screening, and outflow self-absorption.

We are investigating the circumstellar material for a sample of B[e] stars using high spectral resolution data taken in the optical and near-infrared regions with ESO/FEROS and ESO/CRIRES spectrographs, respectively. B[e] stars are surrounded by dense disks of still unknown origin. While optical emission lines from [O I] and [Ca II] reflect the disk conditions close to the star (few stellar radii), the near-infrared data, especially the CO band emission, mirror the characteristics in the molecular part of the disk farther away from the star (several AU). Based on our high resolution spectroscopic data, we seek to derive the density and temperature structure of the disks, as well as their kinematics. This will allow us to obtain a better understanding of their structure, formation history and evolution. Here we present our preliminary results.

Massive stars play an important role in many areas of astrophysics, but numerous details regarding their formation remain unclear. In this paper we present and analyse high resolution (R ~ 30,000) near-infrared 2.3 micron spectra of 20 massive young stellar objects from the RMS database, in the largest such study of CO first overtone bandhead emission to date. We fit the emission under the assumption it originates from a circumstellar disc in Keplerian rotation. We explore three approaches to modelling the physical conditions within the disc - a disc heated mainly via irradiation from the central star, a disc heated mainly via viscosity, and a disc in which the temperature and density are described analytically. We find that the models described by heating mechanisms are inappropriate because they do not provide good fits to the CO emission spectra. We therefore restrict our analysis to the analytic model, and obtain good fits to all objects that possess sufficiently strong CO emission, suggesting circumstellar discs are the source of this emission. On average, the temperature and density structure of the discs correspond to geometrically thin discs, spread across a wide range of inclinations. Essentially all the discs are located within the dust sublimation radius, providing strong evidence that the CO emission originates close to the central protostar, on astronomical unit scales. In addition, we show that the objects in our sample appear no different to the general population of MYSOs in the RMS database, based on their near- and mid-infrared colours. The combination of observations of a large sample of MYSOs with CO bandhead emission and our detailed modelling provide compelling evidence of the presence of small scale gaseous discs around such objects, supporting the scenario in which massive stars form via disc accretion.

HD 327083 is a sgB[e] star that forms a binary system with an orbital semi-major axis of ~1.7 AU. Our previous observations using the VLTI and AMBER in the medium resolution K-band mode spatially resolved the environment of HD 327083. The continuum visibilities obtained indicate the presence of a circumbinary disc. CO bandhead emission was also observed. However, due to the limited spectral resolution of the previous observations, the kinematic structure of the emitting material was not constrained. In this paper, we address this and probe the source of the CO emission with high spectral resolution and spatial precision. We have observed HD 327083 with high spectral resolution (25 & 6 km/s) using AMBER and CRIRES. The observations are compared to kinematical models to constrain the source of the emission. It is shown that the CO bandhead emission can be reproduced using a model of a Keplerian disc with an inclination and size consistent with our previous VLTI observations. The model is compared to AMBER differential phase measurements, which have a precision as high as 30-micro-arcseconds. A differential phase signal corresponding to 0.15 milli-arcseconds (~5 sigma) is seen over the bandhead emission, which is in excellent agreement with the model that fits the CRIRES observations. In comparison, a model of an equatorial outflow, as envisaged in the standard sgB[e] scenario, does not reproduce the observations well. The excellent agreement between the disc model and observations in the spatial and spectral domains is compelling evidence that the CO bandhead emission of HD 327083 originates in a circumbinary Keplerian disc. In contrast, the model of an equatorial outflow cannot reproduce the observations well. This suggests that the standard sgB[e] scenario is not applicable to HD 327083, which supports the hypothesis that the B[e] behaviour of HD 327083 is due to binarity (ABRIDGED).

We present a study of the kinematical properties of a small sample of nearby near-infrared bright massive and intermediate mass young stellar objects using emission lines sensitive to discs and winds. We show for the first time that the broad ($\sim500$kms$^{-1}$) symmetric line wings on the HI Brackett series lines are due to Stark broadening or electron scattering, rather than pure Doppler broadening due to high speed motion. The results are consistent with the presence of a very dense circumstellar environment. In addition, many of these lines show evidence for weak line self-absorption, suggestive of a wind or disc-wind origin for that part of the absorbing material. The weakness of the self-absorption suggests a large opening angle for such an outflow. We also study the fluorescent 1.688$\mu$m FeII line, which is sensitive to dense material. We fitted a Keplerian disc model to this line, and find reasonable fits in all bar one case, in agreement with previous finding for classical Be stars that fluorescent iron transitions are reasonable disc tracers. Overall the picture is one in which these stars still have accretion discs, with a very dense inner circumstellar environment which may be tracing either the inner regions of a disc, or of a stellar wind, and in which ionised outflow is also present. The similarity with lower mass stars is striking, suggesting that at least in this mass range they form in a similar fashion.

While the presence of discs around classical Be stars is well established, their origin is still uncertain. To understand what processes result in the creation of these discs and how angular momentum is transported within them, their physical properties must be constrained. This requires comparing high spatial and spectral resolution data with detailed radiative transfer modelling. We present a high spectral resolution, R~80,000, sub milli-arcsecond precision, spectroastrometric study of the circumstellar disc around the Be star beta CMi. The data are confronted with three-dimensional, NLTE radiative transfer calculations to directly constrain the properties of the disc. Furthermore, we compare the data to disc models featuring two velocity laws; Keperian, the prediction of the viscous disc model, and angular momentum conserving rotation. It is shown that the observations of beta CMi can only be reproduced using Keplerian rotation. The agreement between the model and the observed SED, polarisation and spectroastrometric signature of beta CMi confirms that the discs around Be stars are well modelled as viscous decretion discs.

Massive stars form whilst they are still embedded in dense envelopes. As a result, the roles of rotation, mass loss and accretion in massive star formation are not well understood. This study evaluates the source of the Q-band, lambda=19.5 microns, emission of massive young stellar objects (MYSOs). This allows us to determine the relative importance of rotation and outflow activity in shaping the circumstellar environments of MYSOs on 1000 AU scales. We obtained diffraction limited mid-infrared images of a sample of 20 MYSOs using the VLT/VISIR and Subaru/COMICS instruments. For these 8 m class telescopes and the sample selected, the diffraction limit, ~0.6", corresponds to approximately 1000 AU. We compare the images and the spectral energy distributions (SEDs) observed to a 2D, axis-symmetric dust radiative transfer model that reproduces VLTI/MIDI observations of the MYSO W33A. We vary the inclination, mass infall rate, and outflow opening angle to simultaneously recreate the behaviour of the sample of MYSOs in the spatial and spectral domains. The mid-IR emission of 70 percent of the MYSOs is spatially resolved. In the majority of cases, the spatial extent of their emission and their SEDs can be reproduced by the W33A model featuring an in-falling, rotating dusty envelope with outflow cavities. There is independent evidence that most of the sources which are not fit by the model are associated with ultracompact HII regions and are thus more evolved. We find that, in general, the diverse 20 micron morphology of MYSOs can be attributed to warm dust in the walls of outflow cavities seen at different inclinations. This implies that the warm dust in the outflow cavity walls dominates the Q-band emission of MYSOs. In turn, this emphasises that outflows are an ubiquitous feature of massive star formation.

HD 327083 is a luminous B type star which exhibits emission lines and an infrared excess and is therefore classified as a supergiant B[e] star. In addition, the star is the primary of a close binary system. It is not clear whether the B[e] behaviour of HD 327083 is related to its binarity or its evolutionary state. Here we address this issue by studying its circumstellar environment with high spatial resolution. To this end, we have observed HD 327083 with the VLTI and AMBER in the medium resolution K-band setting. 13CO bandhead emission is detected, confirming HD 327083 is a post-main sequence object. The observations spatially resolve the source of the NIR continuum and the Br-gamma and CO line emission. In addition, differential phase measurements allow us to probe the origin of the observed Br-gamma emission with sub-mas precision. Using geometrical models, we find that the visibilities and closure phases suggest that the close binary system is surrounded by a circum-binary disk. We also find that in the case of the binary HD 327083, the relative sizes of the continuum and Br-gamma emitting regions are different to those of a single supergiant B[e] star where the standard dual outflow scenario is thought to apply. These findings are consistent with the hypothesis that the mass loss of HD 327083 is related to its binary nature.

The majority of the pre-main-sequence Herbig Ae/Be stars reside in binary systems. As these systems are young, their properties may contain an imprint of the star formation process at intermediate masses. However, these systems are generally spatially unresolved, making it difficult to probe their circumstellar environment to search for manifestations of their formation process, such as accretion disks. Here we investigate the formation mechanism of Herbig Ae/Be (HAe/Be) binary systems by studying the relative orientation of their binary orbits and circumstellar disks. We present linear spectropolarimetric observations of HAe/Be stars over the Halpha line, which are used to determine the orientation of their circumstellar disks. In conjunction with data from the literature, we obtain a sample of 20 binaries with known disk position angles (PAs). We compare our disk PA data to a model to investigate whether HAe/Be binary systems and their disks are co-planar. Moreover, in the light of a recent suggestion that some HAe/Be star spectropolarimetric signatures may not necessarily be related to circumstellar disks, we re-assess the relationship between spectropolarimetric signatures and disk PAs. We do this by comparing spectropolarimetric and high spatial resolution observations of young stellar objects. We find that spectropolarimetric observations of pre-main-sequence stars do trace circumstellar disks. This finding is significant above the 3 sigma level. Using this result, our data are entirely consistent with the situation in which HAe/Be binary systems and circumstellar disks are co-planar, while random orientations can be rejected at the 2.2 sigma level. We conclude that HAe/Be binary systems and the disks around the primaries are likely to be co-planar, which is consistent with the notion that these systems form via monolithic collapse and subsequent disk fragmentation (abridged).

We apply mid-infrared spectro-interferometry to the massive young stellar object CRL2136. The observations were performed with the Very Large Telescope Interferometer and the MIDI instrument at a 42m baseline probing angular scales of 50 milli-arcseconds. We model the observed visibilities in parallel with diffraction-limited images at both 24.5micron and in the N-band (with resolutions of 0.6" and 0.3", respectively), as well as the spectral energy distribution. The arcsec-scale spatial information reveals the well-resolved emission from the dusty envelope. By simultaneously modelling the spatial and spectral data, we find that the bulk of the dust emission occurs at several dust sublimation radii (approximately 170 AU). This reproduces the high mid-infrared fluxes and at the same time the low visibilities observed in the MIDI data for wavelengths longward of 8.5micron. However, shortward of this wavelength the visibility data show a sharp up-turn indicative of compact emission. We discuss various potential sources of this emission. We exclude a dust disk being responsible for the observed spectral imprint on the visibilities. A cool supergiant star and an accretion disk are considered and both shown to be viable origins of the compact mid-infrared emission. We propose that CRL2136 is embedded in a dusty envelope, which truncates at several times the dust sublimation radius. A dust torus is manifest in the equatorial region. We find that the spectro-interferometric N-band signal can be reproduced by either a gaseous disk or a bloated central star. If the disk extends to the stellar surface, it accretes at a rate of 3.0 10^(-3) Msun/yr.

We present high resolution (R~50,000) spectroastrometry over the CO 1st overtone bandhead of a sample of seven intermediate/massive young stellar objects. These are primarily drawn from the red MSX source (RMS) survey, a systematic search for young massive stars which has returned a large, well selected sample of such objects. The mean luminosity of the sample is approximately 5 times 10^4 L_\odot, indicating the objects typically have a mass of ~15 solar masses. We fit the observed bandhead profiles with a model of a circumstellar disc, and find good agreement between the models and observations for all but one object. We compare the high angular precision (0.2-0.8 mas) spectroastrometric data to the spatial distribution of the emitting material in the best-fitting models. No spatial signatures of discs are detected, which is entirely consistent with the properties of the best-fitting models. Therefore, the observations suggest that the CO bandhead emission of massive young stellar objects originates in small-scale disks, in agreement with previous work. This provides further evidence that massive stars form via disc accretion, as suggested by recent simulations.

We present B and R band spectroastrometry of a sample of 45 Herbig Ae/Be stars in order to study their binary properties. All but one of the targets known to be binary systems with a separation of ~0.1-2.0 arcsec are detected by a distinctive spectroastrometric signature. Some objects in the sample exhibit spectroastrometric features that do not appear attributable to a binary system. We find that these may be due to light reflected from dusty halos or material entrained in winds. We present 8 new binary detections and 4 detections of an unknown component in previously discovered binary systems. The data confirm previous reports that Herbig Ae/Be stars have a high binary fraction, 74+/-6 per cent in the sample presented here. We use a spectroastrometric deconvolution technique to separate the spatially unresolved binary spectra into the individual constituent spectra. The separated spectra allow us to ascertain the spectral type of the individual binary components, which in turn allows the mass ratio of these systems to be determined. In addition, we appraise the method used and the effects of contaminant sources of flux. We find that the distribution of system mass ratios is inconsistent with random pairing from the Initial Mass Function, and that this appears robust despite a detection bias. Instead, the mass ratio distribution is broadly consistent with the scenario of binary formation via disk fragmentation.

The prototype of the Beta Cep class of pulsating stars, Beta Cep, rotates relatively slowly, and yet displays episodic Halpha emission. Such behaviour is typical of a classical Be star. For some time this posed a contradiction to our understanding of the Be phenomena as rapid rotation is thought to be a prerequisite for the emission phases of Be stars. Recent work has demonstrated that the Halpha emission is in fact due to a close companion (separation ~0.25") of the star. This resolves the apparent enigma if this close companion is a classical Be star. We aim to test the hypothesis that this close companion is a valid Be star by determining properties such as its spectral type and vsin i. We employed spectroastrometry to investigate the close binary system. Using the spectroastrometric signatures observed, we split the composite binary spectra into its constituent spectra in the B band and R band. The split spectra allow us to estimate spectral types of the binary components. We find that the primary of the close binary system has a spectral type of B2III and the secondary a spectral type of B5Ve. From the relationship between mass and spectral type, we determine the masses of the binary components to be M_p = 12.6+/-3.2 M_* and M_s = 4.4+/-0.7 M_* respectively. The spectroastrometric data allow some constraint on the orbit, and we suggest a moderate revision to the previously determined orbit. We confirm that the primary of the system is a slow rotator (vsin i =29 +43/-29 km/s), while the secondary rotates significantly faster, at a vsin i =230+/-45 km/s. We show that the close companion to the beta Cep primary is a valid classical Be star. It has a spectral type of B5Ve and is a relatively fast rotator.