Young accreting stars drive outflows that collimate into jets, which can be seen hundreds of au from their driving sources. Accretion and outflow activity cease with system age, and it is believed that magneto-centrifugally launched disk winds are critical agents in regulating accretion through the protoplanetary disk. Protostellar jets are well studied in classical T~Tauri stars ($M_\star\lesssim2\,M_\odot$), while few nearby ($d\lesssim150\,$pc) intermediate-mass stars ($M_\star=2-10\,M_\odot$), known as Herbig Ae/Be stars, have detected jets. We report VLT/MUSE observations of the Herbig~Ae/Be star HD~100546 and the discovery of a protostellar jet. The jet is similar in appearance to jets driven by low-mass stars and compares well with the jet of HD~163296, the only other known optical jet from a nearby Herbig~Ae/Be star. We derive a (one-sided) mass-loss rate in the jet of $\log \dot{M}_{jet} \sim -9.5$ (in $M_\odot$\,yr$^{-1}$) and a ratio of outflow to accretion of roughly $3\times10^{-3}$, which is lower than that of CTTS jets. The discovery of the HD~100546 jet is particularly interesting because the protoplanetary disk around HD~100546 shows a large radial gap, spiral structure, and might host a protoplanetary system. A bar-like structure previously seen in H$\alpha$ with VLT/SPHERE shares the jet position angle, likely represents the base of the jet, and suggests a jet-launching region within about 2\,au. We conclude that the evolution of the disk at radii beyond a few au does not affect the ability of the system to launch jets.

The intermediate-mass pre-main sequence Herbig Ae/Be stars are key to understanding the differences in formation mechanisms between low- and high-mass stars. The study of the general properties of these objects is hampered by the fact that few and mostly serendipitously discovered sources are known. Our goal is to identify new Herbig Ae/Be candidates to create a homogeneous and well defined catalogue of these objects. We have applied machine learning techniques to 4,150,983 sources with data from Gaia DR2, 2MASS, WISE, and IPHAS or VPHAS+. Several observables were chosen to identify new Herbig Ae/Be candidates based on our current knowledge of this class, which is characterised by infrared excesses, photometric variabilities, and H$\alpha$ emission lines. Classical techniques are not efficient for identifying new Herbig Ae/Be stars mainly because of their similarity with classical Be stars, with which they share many characteristics. By focusing on disentangling these two types of objects, our algorithm has also identified new classical Be stars. We have obtained a large catalogue of 8470 new pre-main sequence candidates and another catalogue of 693 new classical Be candidates with a completeness of $78.8\pm1.4\%$ and $85.5\pm1.2\%$, respectively. Of the catalogue of pre-main sequence candidates, at least 1361 sources are potentially new Herbig Ae/Be candidates according to their position in the Hertzsprung-Russell diagram. In this study we present the methodology used, evaluate the quality of the catalogues, and perform an analysis of their flaws and biases. For this assessment, we make use of observables that have not been accounted for by the algorithm and hence are selection-independent, such as coordinates and parallax based distances. The catalogue of new Herbig Ae/Be stars that we present here increases the number of known objects of the class by an order of magnitude.

Aims. To demonstrate that `INDICATE' is a powerful spatial analysis tool which when combined with kinematic data from Gaia DR2 can be used to robustly probe star formation history. Methods. We compared the dynamic & spatial distributions of young stellar objects (YSOs) at various evolutionary stages in NGC 2264 using Gaia DR2 proper motion data and INDICATE. Results. The dynamic & spatial behaviours of YSOs at different evolutionary stages are distinct. Dynamically, Class II YSOs predominately have non-random trajectories that are consistent with known substructures, whereas Class III YSOs have random trajectories with no clear expansion or contraction patterns. Spatially, there is a correlation between the evolutionary stage and source concentration: 69.4% of Class 0/I, 27.9% of Class II, and 7.7% of Class III objects are found to be clustered. The proportion of YSOs clustered with objects of the same class also follows this trend. Class 0/I objects are both found to be more tightly clustered with the general populous/objects of the same class than Class IIs and IIIs by a factor of 1.2/4.1 and 1.9/6.6, respectively. An exception to these findings is within 0.05deg of S Mon where Class III objects mimic the behaviours of Class II sources across the wider cluster region. Our results suggest (i) current YSOs distributions are a result of dynamical evolution, (ii) prolonged star formation has been occurring sequentially, and (iii) stellar feedback from S Mon is causing YSOs to appear as more evolved sources. Conclusions. Designed to provide a quantitative measure of clustering behaviours, INDICATE is a powerful tool with which to perform rigorous spatial analyses. Our findings are consistent with what is known about NGC 2264, effectively demonstrating that when combined with kinematic data from Gaia DR2 INDICATE can be used to robustly study the star formation history of a cluster.

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.

Aims: The inner regions of high-mass protostars are often invisible in the near-infrared. We aim to investigate the inner gaseous disc of IRAS11101-5829 through scattered light from the outflow cavity walls. Methods: We observed the environment of the high-mass young stellar object IRAS11101-5829 and the closest knots of its jet, HH135-136, with the VLT/SINFONI. We also retrieved archival data from the high-resolution long-slit spectrograph VLT/X-shooter. Results: We detect the first three bandheads of the $\upsilon=2-0$ CO vibrational emission for the first time in this object. It is coincident with continuum and Br$\gamma$ emission and extends up to $\sim10000$ au to the north-east and $\sim10 000$ au to the south-west. The line profiles have been modelled as a Keplerian rotating disc assuming a single ring in LTE. The model output gives a temperature of $\sim3000$ K, a CO column density of $\sim1\times10^{22}\mathrm{ cm^{-2}}$, and a projected Keplerian velocity $v_\mathrm{K}\sin i_\mathrm{disc} \sim 25\mathrm{ km s^{-1}}$, which is consistent with previous modelling in other high-mass protostars. In particular, the low value of $v_\mathrm{K}\sin i_\mathrm{disc}$ suggests that the disc is observed almost face-on, whereas the well-constrained geometry of the jet imposes that the disc must be close to edge-on. This apparent discrepancy is interpreted as the CO seen reflected in the mirror of the outflow cavity wall. Conclusions: From both jet geometry and disc modelling, we conclude that all the CO emission is seen through reflection by the cavity walls and not directly. This result implies that in the case of highly embedded objects, as for many high-mass protostars, line profile modelling alone might be deceptive and the observed emission could affect the derived physical and geometrical properties; in particular the inclination of the system can be incorrectly interpreted.

The rarity of young massive stars combined with the fact that they are often deeply embedded has limited the understanding of their formation. Ground based mid-infrared (IR) interferometry is one way of securing the spatial resolution required to study massive young stellar objects (MYSOs) and as the spatial-frequency coverage of such observations is often incomplete, direct-imaging can be supplementary to such a dataset. By consolidating these observations with modelling, the features of a massive protostellar environment can be constrained. This work simultaneously fits the aforementioned observations and a spectral energy distribution (SED) with a 2.5D radiative transfer model, providing an extensive view of the physical characteristics of the accreting regions of the MYSO G305.20+0.21. The high-resolution observations were obtained using the Very Large Telescope's MIDI and VISIR instruments, producing visibilities in the N-band and near-diffraction-limited imaging in the Q-band respectively. A model including a central protostar with a luminosity of $\sim$5$\times$10$^{4}$L$_{\odot}$ surrounded by a low-density bipolar cavity, a flared 1M$_{\odot}$ disk and an envelope provides a sufficient fit all three types of observation. The need to include a disk in the model implies that this MYSO follows a scaled-up version of the low-mass star formation process. The weak silicate absorption feature within the SED requires low-density envelope cavities to be successfully fit and is an atypical characteristic in comparison to previously studied MYSOs. Additionally, the inner radius of the disk must be three times the dust sublimation radius to satisfy the MIDI visibilities. The low density, low extinction environment implies the object is a more evolved MYSO and this combined with large inner radius of the disk suggests that it could be an example of a transitional disk around an MYSO.

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 introduce our new code MYOSOTIS (Make Your Own Synthetic ObservaTIonS) which is designed to produce synthetic observations from simulated clusters. The code can synthesise observations from both ground- and spaced-based observatories, for a range of different filters, observational conditions and angular/spectral resolution. In this paper, we highlight some of the features of MYOSOTIS, creating synthetic observations from young massive star clusters. Our model clusters are simulated using nbody6 code and have different total masses, half-mass radii, and binary fractions. The synthetic observations are made at the age of 2 Myr with Solar metallicity and under different extinction conditions. For each cluster, we create synthetic images of the Hubble Space Telescope (HST) in the visible (WFPC2/F555W) as well as Very Large Telescopes (VLT) in the nearIR (SPHERE/IRDIS/Ks). We show how MYOSOTIS can be used to look at mass function (MF) determinations. For this aim we re-estimate stellar masses using a photometric analysis on the synthetic images. The synthetic MF slopes are compared to their actual values. Our photometric analysis demonstrate that depending on the adopted filter, extinction, angular resolution and pixel sampling of the instruments, the power-law index of the underlying MFs can be shallower than the observed ones by at least 0.25 dex which is in agreement with the observed discrepancies reported in the literature, specially for young star clusters.

Aims. To present the new novel statistical clustering tool 'INDICATE' which assesses and quantifies the degree of spatial clustering of each object in a dataset, discuss its applications as a tracer of morphological stellar features in star forming regions, and to look for these features in the Carina Nebula (NGC 3372). Results. We successfully recover known stellar structure of the Carina Nebula, including the 5 young star clusters in this region. Four sub-clusters contain no, or very few, stars with a degree of association above random which suggests they may be fluctuations in the field rather than real clusters. In addition we find: (1) Stars in the NW and SE regions have significantly different clustering tendencies, which is reflective of differences in the apparent star formation activity in these regions. Further study is required to ascertain the physical origin of the difference; (2) The different clustering properties between these two regions are even more pronounced for OB stars; (3) There are no signatures of classical mass segregation present in the SE region - massive stars here are not spatially concentrated together above random; (4) Stellar concentrations are more frequent around massive stars than typical for the general population, particularly in the Tr14 cluster; (5) There is a relation between the concentration of OB stars and the concentration of (lower mass) stars around OB stars in the centrally concentrated Tr14 and Tr15, but no such relation exists in Tr16. We conclude this is due to the highly sub-structured nature of Tr16. Conclusions. INDICATE is a powerful new tool employing a novel approach to quantify the clustering tendencies of individual objects in a dataset within a user-defined parameter space. As such it can be used in a wide array of data analysis applications.

We present the first search for binary companions of Massive Young Stellar Objects (MYSOs) using AO-assisted $K$ band observations, with NaCo at the VLT. We have surveyed 32 MYSOs from the RMS catalogue, probing the widest companions, with a physical separation range of 400 - 46,000 au, within the predictions of models and observations for multiplicity of MYSOs. Statistical methods are employed to discern whether these companions are physical rather than visual binaries. We find 18 physical companions around 10 target objects, amounting to a multiplicity fraction of 31$\pm$8\% and a companion fraction of 53$\pm$9\%. For similar separation and mass ratio ranges, MYSOs seem to have more companions than T Tauri or O stars, respectively. This suggests that multiplicity increases with mass and decreases with evolutionary stage. We compute very rough estimates for the mass ratios from the $K$ band magnitudes, and these appear to be generally larger than 0.5. This is inconsistent with randomly sampling the IMF, as predicted by the binary capture formation theory. Finally, we find that MYSOs with binaries do not show any different characteristics to the average MYSO in terms of luminosity, distance, outflow or disc presence.

(Abridged) The detection of forming planets in disks around young stars remains elusive, and state-of-the-art observational techniques provide somewhat ambiguous results. It has been reported that the pre-transitional T Tauri star LkCa 15 could host three planets; candidate planet b is in the process of formation, as inferred from its H$\alpha$ emission. However, a more recent work casts doubts on the planetary nature of the previous detections. We have observed LkCa 15 with ISIS/WHT. The spectrograph's slit was oriented towards the last reported position of LkCa 15 b (parallel direction) and 90degr from that (perpendicular). The photocenter and full width half maximum (FWHM) of the Gaussians fitting the spatial distribution at H$\alpha$ and the adjacent continuum were measured. A well-known binary (GU CMa) was used as a calibrator to test the spectro-astrometric performance of ISIS/WHT, recovering consistent photocenter and FWHM signals. However, the photocenter shift predicted for LkCa 15 b is not detected, but the FWHM in H$\alpha$ is broader than in the continuum for both slit positions. Our simulations show that the photocenter and FWHM observations cannot be explained simultaneously by an accreting planet. In turn, both spectro-astrometric observations are naturally reproduced from a roughly symmetric Halpha emitting region centered on the star and extent comparable to the orbit originally attributed to the planet at several au. The extended H$\alpha$ emission around LkCa 15 could be related to a variable disk wind, but additional multi-epoch data and detailed modeling are necessary to understand its physical nature. Spectro-astrometry in H$\alpha$ is able to test the presence of accreting planets and can be used as a complementary technique to survey planet formation in circumstellar disks.

(abridged) The star formation rate (SFR) linearly correlates with the amount of dense gas mass (Mdg) involved in the formation of stars both for distant galaxies and clouds in our Galaxy. Similarly, the mass accretion rate (Macc) and the disk mass (Mdisk) of young, Class II stars are also linearly correlated. We plotted the corresponding observational data together, finding a statistically significant correlation that spans ~ 16 orders of magnitude. This probably represents one of the widest ranges of any empirical correlation known, encompassing galaxies that are several kpc in size, pc-size star-forming clouds within our Galaxy, down to young, pre-main sequence stars with au-size protoplanetary disks. We propose a bottom-up hypothesis suggesting that a relation between Macc and the total circumstellar mass surrounding Class 0/I sources (Mcs; disk+envelope) drives the correlation in protostellar-hosting clouds and cloud-hosting galaxies. This is consistent with the fact that the SFRs derived for clouds over a timescale of 2 Myr can be roughly recovered from the sum of instantaneous Macc values of the protostars embedded within them, implying that galactic SFRs averaged over ~ 10-100 Myr should be constant over this period too. Moreover, the sum of the Mcs values directly participating in the formation of the protostellar population in a cloud likely represents a non-negligible fraction of the Mdg within the cloud. If such fraction is ~ 1-35% of the Mdg associated with star-forming clouds and galaxies, then the global correlation for all scales has a near unity slope and an intercept consistent with the (proto-)stellar accretion timescale, Mcs/Macc. Therefore, an additional critical test of our hypothesis is that the Macc-Mdisk correlation for Class II stars should also be observed between Macc and Mcs for Class 0/I sources with similar slope and intercept.

We use Gaia Data Release 2 (DR2) to place 252 Herbig Ae/Be stars in the HR diagram and investigate their characteristics and properties. For all known Herbig Ae/Be stars with parallaxes in Gaia DR2, we collected their atmospheric parameters and photometric and extinction values from the literature. To these data we added near- and mid-infrared photometry, collected H$\alpha$ equivalent widths and line profiles, and their binarity status. In addition, we developed a photometric variability indicator from Gaia's DR2 information. We provide masses, ages, luminosities, distances, variabilities and infrared excesses homogeneously derived for the most complete sample of Herbig Ae/Be stars to date. We find that high mass stars have a much smaller infrared excess and have much lower optical variabilities compared to lower mass stars, with the break at around 7M$_{\odot}$. H$\alpha$ emission is generally correlated with infrared excess, with the correlation being stronger for infrared emission at wavelengths tracing the hot dust closest to the star. The variability indicator as developed by us shows that $\sim$25% of all Herbig Ae/Be stars are strongly variable. We observe that the strongly variable objects display doubly peaked H$\alpha$ line profiles, indicating an edge-on disk. The fraction of strongly variable Herbig Ae stars is close to that found for A-type UX Ori stars. It had been suggested that this variability is in most cases due to asymmetric dusty disk structures seen edge-on. The observation here is in strong support of this hypothesis. Finally, the difference in dust properties occurs at 7M$_{\odot}$, while various properties traced at UV/optical wavelengths differ at 3M$_{\odot}$. The latter has been linked to different accretion mechanisms at work whereas the differing infrared properties and variabilities are related to different (dust-)disk dispersal mechanisms.

We present a new radio-selected sample of PNe from the CORNISH survey. This is a radio continuum survey of the inner Galactic plane covering Galactic longitude, $10^\circ <l< 65^\circ$ and latitude, $|b| < 1^\circ$ with a resolution of 1.5$^{"}$ and sensitivity better than 0.4 mJy/beam. The radio regime, being unbiased by dust extinction, allows for a more complete sample selection, especially towards the Galactic mid-plane. Visual inspection of the CORNISH data, in combination with data from multi-wavelength surveys of the Galactic plane, allowed the CORNISH team to identify 169 candidate PNe. Here, we explore the use of multi-wavelength diagnostic plots and analysis to verify and classify the candidate PNe. We present the multi-wavelength properties of this new PNe sample. We find 90 new PNe, of which 12 are newly discovered and 78 are newly classified as PN. A further 47 previously suspected PNe are confirmed as such from the analysis presented here and 24 known PNe are detected. Eight sources are classified as possible PNe or other source types. Our sample includes a young sub-sample, with physical diameters < 0.12 pc, brightness temperatures (> 1000 K) and located closer than 7 kpc. Within this sample is a water-maser PN with a spectral index of $-0.55\pm 0.08$, which indicates non-thermal radio emission. Such a radio-selected sample, unaffected by extinction, will be particularly useful to compare with population synthesis models and should contribute to the understanding of the formation and evolution of PNe.

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.

HD 100546 is one of the few known pre-main-sequence stars that may host a planetary system in its disk. We analyze new VLT/SPHERE/ZIMPOL polarimetric images of HD 100546 with filters in H$\alpha$ and the adjacent continuum. We have probed the disk gap and the surface layers of the outer disk, covering a region < 500 mas (< 55 au at 109 pc) from the star, at an angular resolution of ~ 20 mas. Our data show an asymmetry: the SE and NW regions of the outer disk are more polarized than the SW and NE. This can be explained from a preferential scattering angle close to 90$^o$, consistent with previous polarization images. The outer disk extends from 13 $\pm$ 2 to 45 $\pm$ 9 au, with a position angle and inclination of 137 $\pm$ 5$^o$ and 44 $\pm$ 8$^o$. The comparison with previous estimates suggests that the disk inclination could increase with the stellocentric distance, although the different measurements are still consistent within the error bars. In addition, no direct signature of the innermost candidate companion is detected from polarimetry, confirming recent results based on intensity imagery. We set an upper limit to its mass accretion rate < 10$^{-8}$ M$_{\odot}$/yr for a sub-stellar mass of 15M$_{Jup}$. Finally, we report the first detection (> 3$\sigma$) of a ~ 20 au bar-like structure that crosses the gap through the central region of HD 100546. It is tentatively suggested that the bar could be dust dragged by infalling gas that radially flows from the outer disk. This could represent an exceptional case in which a small-scale radial inflow is observed in a single system. If this scenario is confirmed, it could explain the presence of atomic gas in the inner disk that would otherwise accrete on to the central star on a timescale of a few months/years, as previously indicated from spectro-interferometric data, and could be related with additional (undetected) planets.

We present medium-resolution (R~7000) near-infrared echelle spectroscopic data for 36 MYSOs drawn from the Red MSX Source (RMS) survey. This is the largest sample observed at this resolution at these wavelengths of MYSOs to date. The spectra are characterized mostly by emission from hydrogen recombination lines and accretion diagnostic lines. One MYSO shows photospheric HI absorption, a comparison with spectral standards indicates that the star is an A type star with a low surface gravity, implying that the MYSOs are probably swollen, as also suggested by evolutionary calculations. An investigation of the Br$\gamma$ line profiles shows that most are in pure emission, while 13$\pm$5 % display P Cygni profiles, indicative of outflow, while less than 8$\pm$4 % have inverse P Cygni profiles, indicative of infall. These values are comparable with investigations into the optically bright Herbig Be stars, but not with those of Herbig Ae and T Tauri stars, consistent with the notion that the more massive stars undergo accretion in a different fashion than lower mass objects which are undergoing magnetospheric accretion. Accretion luminosities and rates as derived from the Br$\gamma$ line luminosities agree with results for lower mass sources, providing tentative evidence for massive star formation theories based on scaling of low-mass scenarios. We present Br$\gamma$/Br12 line profile ratios exploiting the fact that optical depth effects can be traced as a function of Doppler shift across the lines. These show that the winds of MYSOs in this sample are nearly equally split between constant, accelerating, and decelerating velocity structures. There are no trends between the types of features we see and bolometric luminosities or near-infrared colours.

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.

Solar-mass stars form via circumstellar disk accretion (disk-mediated accretion). Recent findings indicate that this process is likely episodic in the form of accretion bursts, possibly caused by disk fragmentation. Although it cannot be ruled out that high-mass young stellar objects (HMYSOs; $M>$8 M$_\odot$, $L_{bol}>$5$\times$10$^3$ L$_\odot$) arise from the coalescence of their low-mass brethren, latest results suggest that they more likely form via disks. Accordingly, disk-mediated accretion bursts should occur. Here we report on the discovery of the first disk-mediated accretion burst from a $\sim$20 M$_\odot$ HMYSO. Our near-infrared images show the brightening of the central source and its outflow cavities. Near-infrared spectroscopy reveals emission lines typical of accretion bursts in low-mass protostars, but orders of magnitude more luminous. Moreover, the energy released and the inferred mass-accretion rate are also orders of magnitude larger. Our results identify disk accretion as the common mechanism of star formation across the entire stellar mass spectrum.

Yellow hypergiants are rare and represent a fast evolutionary stage of massive evolved stars. That evolutionary phase is characterised by a very intense mass loss, the understanding of which is still very limited. Here we report ALMA Compact Array observations of a 50$"$-mosaic toward the Fried Egg nebula, around one of the few Galactic yellow hypergiants IRAS 17163-3907. The emission from the $^{12}$CO J=2-1 line, H30$\alpha$ recombination line, and continuum is imaged at a resolution of $\sim$8$"$, revealing the morphology of the molecular environment around the star. The continuum emission is unresolved and peaks at the position of the star. The radio recombination line H30$\alpha$ shows unresolved emission at the star, with an approximately gaussian spectrum centered on a velocity of 21$\pm$3~km/s with a width of 57$\pm$6~km/s. In contrast, the CO 2-1 emission is complex and decomposes into several components beyond the contamination from interstellar gas in the line of sight. The CO spectrum toward the star is a broad plateau, centered at the systemic velocity of +18 km/s and with an expansion velocity of 100$\pm$10 km/s. Assuming isotropic and constant mass-loss, we estimate a mass-loss rate of 8$\pm$1.5 $\times10^{-5}$~M$_\odot$ yr$^{-1}$. At a radius of 25$"$ from the star, we detect CO emission associated with the dust ring previously imaged by {\it Herschel}. The kinematics of this ring, however, is not consistent with an expanding shell, but show a velocity gradient of $v_{sys} \pm$20 km/s. In addition, we find a puzzling bright feature radially connecting the star to the CO ring, at a velocity of +40 km/s relative to the star. This spur feature may trace a unidirectional ejection event from the star. Our ACA observations reveal the complex morphology around IRAS 17163 and illustrate the breakthroughs that ALMA will bring to the field of massive stellar evolution.

The Herbig Ae/Be stars (HAeBes) allow an exploration of the properties of Pre-Main Sequence(PMS) stars above the low-mass range ($<2{\rm\thinspace M_{\odot}}$) and those bordering the high-mass range ($>8{\rm\thinspace M_{\odot}}$). This paper is the second in a series exploring accretion in 91 HAeBes with Very Large Telescope/X-shooter spectra. Equivalent width measurements are carried out on 32 different lines, spanning the UV to NIR, in order to obtain their line luminosities. The line luminosities were compared to accretion luminosities, which were determined directly from measurements of an UV-excess. When detected, emission lines always demonstrate a correlation with the accretion luminosity, regardless of detection frequency. The average relationship between accretion luminosity and line luminosity is found to be ${\thinspace L_{\rm acc}}\propto{\thinspace L_{\rm line}}^{1.16 \pm 0.15}$. This is in agreement with the findings in Classical T Tauri stars, although the HAeBe relationship is generally steeper, particularly towards the Herbig Be mass range. Since all observed lines display a correlation with the accretion luminosity, all of them can be used as accretion tracers. This has increased the number of accretion diagnostic lines in HAeBes tenfold. However, questions still remain on the physical origin of each line, which may not be due to accretion.

Aims. The structural inhomogeneities and kinematics of massive star nebulae are tracers of their mass-loss history. We conduct a three-dimensional morpho-kinematic analysis of the ejecta of eta Car outside its famous Homunculus nebula. Methods. We carried out the first large-scale integral field unit observations of eta Car in the optical, covering a field of view of 1'x1' centered on the star. Observations with the Multi Unit Spectroscopic Explorer (MUSE) at the Very Large Telescope (VLT) reveal the detailed three-dimensional structure of eta Car's outer ejecta. Morpho-kinematic modeling of these ejecta is conducted with the code SHAPE. Results. The largest coherent structure in eta Car's outer ejecta can be described as a bent cylinder with roughly the same symmetry axis as the Homunculus nebula. This large outer shell is interacting with the surrounding medium, creating soft X-ray emission. We establish the shape and extent of the ghost shell in front of the southern Homunculus lobe and confirm that the NN condensation can best be modeled as a bowshock in the orbital/equatorial plane. Conclusions. The SHAPE modeling of the MUSE observations indicates that the kinematics of the outer ejecta measured with MUSE can be described by a spatially coherent structure, and this structure also correlates with the extended soft X-ray emission associated with the outer debris field. The ghost shell just outside the southern Homunculus lobe hints at a sequence of eruptions within the time frame of the Great Eruption from 1837-1858 or possibly a later shock/reverse shock velocity separation. Our 3D morpho-kinematic modeling and the MUSE observations constitute an invaluable dataset to be confronted with future radiation-hydrodynamics simulations. Such a comparison may shed light on the yet elusive physical mechanism responsible for eta Car-like eruptions.

This work presents CHARA/VEGA H$\alpha$ spectro-interferometry (R ~ 6000, and $\lambda$/2B ~ 1 mas) of HD 179218 and HD 141569, doubling the sample of Herbig Ae/Be (HAeBe) stars for which this type of observations is available so far. The observed H$\alpha$ emission is spatially unresolved, indicating that the size of the H$\alpha$ emitting region is smaller than ~ 0.21 and 0.12 au for HD 179218 and HD 141529 (~ 15 and 16 R*, respectively). This is smaller than for the two other HAeBes previously observed with the same instrumentation. Two different scenarios have been explored in order to explain the compact line emitting regions. A hot, several thousand K, blackbody disc is consistent with the observations of HD 179218 and HD 141569. Magnetospheric accretion (MA) is able to reproduce the bulk of the H$\alpha$ emission shown by HD 179218, confirming previous estimates from MA shock modelling with a mass accretion rate of 10^-8 Msun/yr, and an inclination to the line of sight between 30 and 50 degr. The H$\alpha$ profile of HD 141569 cannot be fitted from MA due to the high rotational velocity of this object. Putting the CHARA sample together, a variety of scenarios is required to explain the H$\alpha$ emission in HAeBe stars -compact or extended, discs, accretion, and winds-, in agreement with previous Br$\gamma$ spectro-interferometric observations.

The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.

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.

To probe the circumstellar environment of IRAS 13481-6124, a 20 M_sun high-mass young stellar object (HMYSO) with a parsec-scale jet and accretion disc, we investigate the origin of its Br\gamma-emission line through NIR interferometry. We present the first AMBER/VLTI observations of the Br\gamma-emitting region in an HMYSO at R~1500. Our AMBER/VLTI observations reveal a spatially and spectrally resolved Br\gamma-line in emission with a strong P Cygni profile, indicating outflowing matter with a terminal velocity of ~500 km/s. Visibilities, differential phases, and closure phases are detected in our observations within the spectral line and in the adjacent continuum. Both total visibilities (continuum plus line emitting region) and pure-line visibilities indicate that the Br\gamma-emitting region is more compact (2-4 mas in diameter or ~6-13 au at 3.2 kpc) than the continuum-emitting region (~5.4 mas or ~17 au). The absorption feature is also spatially resolved at the longest baselines (81 and 85 m) and has a visibility that is slightly smaller than the continuum-emitting region. The differential phases at the four longest baselines display an \u2018S\u2019-shaped structure across the line, peaking in the blue- and red-shifted high-velocity components. The calibrated photocentre shifts are aligned with the known jet axis, i.e they are probably tracing an ionised jet. The high-velocity components (v_r~100-500 km/s) are located far from the source, whereas the low-velocity components (0-100 km/s) are observed to be closer, indicating a strong acceleration of the gas flow in the inner 10 au. Finally, a non-zero closure phase along the continuum is detected. By comparing our observations with the synthetic images of the continuum around 2.16 um, we confirm that this feature originates from the asymmetric brightness distribution of the continuum owing to the inclination of the inner disc.

This chapter discusses the use and possibilities of optical and infrared interferometry to study star formation. The chapter starts with a brief overview of the star formation process and highlights the open questions from an observational point of view. These are found at the smallest scales, as this is, inevitably, where all the action such as accretion and outflows, occurs. We then use basic astrophysical concepts to assess which scales and conditions can be probed with existing interferometric set-ups for which we use the ESO/VLTI instrument suite as example. We will concentrate on the more massive stars observed at high resolution with continuum interferometry. Throughout, some of the most recent interferometric results are used as examples of the various processes discussed.

Not all stars exhibiting the optical spectral characteristics of B[e] stars share the same evolutionary stage. The Galactic B[e] star MWC 137 is a prime example of an object with uncertain classification, with previous work suggesting pre- and post-main sequence classification. Our goal is to settle this debate and provide reliable evolutionary classification. Integral field spectrograph observations with VLT MUSE of the cluster SH 2-266 are used to analyze the nature of MWC 137. A collimated outflow is discovered that is geometrically centered on MWC 137. The central position of MWC 137 in the cluster SH 2-266 within the larger nebula suggests strongly that it is a member of this cluster and that it is both at the origin of the nebula and the newly discovered jet. Comparison of the color-magnitude diagram of the brightest cluster stars with stellar evolutionary models results in a distance of about 5.2$\pm$1.4 kpc. We estimate that the cluster is at least 3 Myr old. The jet extends over 66" (1.7 pc) projected on the plane of the sky, shows several knots, and projected velocities of up to $\pm$450 km s$^{-1}$. From the Balmer emission line decrement of the diffuse intracluster nebulosity we determine E(B-V)=1.4 mag for the inner 1' cluster region. The spectral energy distribution of the brightest cluster stars yield a slightly lower extinction of E(B-V)~1.2 mag. The extinction towards MWC 137 is estimated to be E(B-V)~1.8 mag (A$_V$~5.6 mag). Our analysis of the optical and near-infrared color-magnitude and color-color diagrams suggests a post-main sequence stage of MWC 137. The existence of a jet in this object implies the presence of an accretion disk.

We present spatially and spectrally resolved Br-gamma emission around the planet-hosting, transitional Herbig Ae/Be star HD 100546. Aiming to gain insight into the physical origin of the line in possible relation to accretion processes, we carried out Br-gamma spectro-interferometry using AMBER/VLTI from three different baselines achieving spatial and spectral resolutions of 2-4 mas and 12000. The Br-gamma visibility is larger than that of the continuum for all baselines. Differential phases reveal a shift between the photocentre of the Br-gamma line -displaced 0.6 mas (0.06 au at 100 pc) NE from the star- and that of the K-band continuum emission -displaced 0.3 mas NE from the star. The photocentres of the redshifted and blueshifted components of the Br-gamma line are located NW and SE from the photocentre of the peak line emission, respectively. Moreover, the photocentre of the fastest velocity bins within the spectral line tends to be closer to that of the peak emission than the photocentre of the slowest velocity bins. Our results are consistent with a Br-gamma emitting region inside the dust inner rim (<0.25 au) and extending very close to the central star, with a Keplerian, disc-like structure rotating counter-clockwise, and most probably flared (25 deg). Even though the main contribution to the Br-gamma line does not come from gas magnetically channelled on to the star, accretion on to HD 100546 could be magnetospheric, implying a mass accretion rate of a few 10^-7 Msun/yr. This value indicates that the observed gas has to be replenished on time-scales of a few months to years, perhaps by planet-induced flows from the outer to the inner disc as has been reported for similar systems.

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.

Herbig Ae/Be stars span a key mass range that links low and high mass stars, and thus provide an ideal window from which to explore their formation. This paper presents VLT/X-Shooter spectra of 91 Herbig Ae/Be stars, HAeBes; the largest spectroscopic study of HAeBe accretion to date. A homogeneous approach to determining stellar parameters is undertaken for the majority of the sample. Measurements of the ultra-violet (UV) are modelled within the context of magnetospheric accretion, allowing a direct determination of mass accretion rates. Multiple correlations are observed across the sample between accretion and stellar properties: the youngest and often most massive stars are the strongest accretors, and there is an almost 1:1 relationship between the accretion luminosity and stellar luminosity. Despite these overall trends of increased accretion rates in HAeBes when compared to classical T Tauri stars, we also find noticeable differences in correlations when considering the Herbig Ae and Herbig Be subsets. This, combined with the difficulty in applying a magnetospheric accretion model to some of the Herbig Be stars, could suggest that another form of accretion may be occurring within the Herbig Be mass range.

Correlations between the accretion luminosity and emission line luminosities (L_acc and L_line) of pre-main sequence (PMS) stars have been published for many different spectral lines, which are used to estimate accretion rates. Despite the origin of those correlations is unknown, this could be attributed to direct or indirect physical relations between the emission line formation and the accretion mechanism. This work shows that all (near-UV/optical/near-IR) L_acc-L_line correlations are the result of the fact that the accretion luminosity and the stellar luminosity (L_star) are correlated, and are not necessarily related with the physical origin of the line. Synthetic and observational data are used to illustrate how the L_acc-L_line correlations depend on the L_acc-L_star relationship. We conclude that because PMS stars show the L_acc-L_star correlation immediately implies that L_acc also correlates with the luminosity of all emission lines, for which the L_acc-L_line correlations alone do not prove any physical connection with accretion but can only be used with practical purposes to roughly estimate accretion rates. When looking for correlations with possible physical meaning, we suggest that L_acc/L_star and L_line/L_star should be used instead of L_acc and L_line. Finally, the finding that L_acc has a steeper dependence on L_star for T-Tauri stars than for intermediate-mass Herbig Ae/Be stars is also discussed. That is explained from the magnetospheric accretion scenario and the different photospheric properties in the near-UV.

Through CO mm-line and optical spectroscopy, we investigate the properties of the Fried Egg nebula IRAS 17163-3907, which has recently been proposed to be one of the rare members of the yellow hypergiant class. The CO J=2-1 and J=3-2 emission arises from a region within 20" of the star and is clearly associated with the circumstellar material. The CO lines show a multi-component asymmetrical profile, and an unexpected velocity gradient is resolved in the east-west direction, suggesting a bipolar outflow. This is in contrast with the apparent symmetry of the dust envelope as observed in the infrared. The optical spectrum of IRAS 17163-3907 between 5100 and 9000 {\AA} was compared with that of the archetypal yellow hypergiant IRC+10420 and was found to be very similar. These results build on previous evidence that IRAS 17163-3907 is a yellow hypergiant.

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.

HD 142527 is a young pre-main sequence star with properties indicative of the presence of a giant planet or/and a low-mass stellar companion. We have analyzed an X-Shooter/Very Large Telescope spectrum to provide accurate stellar parameters and accretion rate. The analysis of the spectrum, together with constraints provided by the SED fitting, the distance to the star (140 +- 20 pc) and the use of evolutionary tracks and isochrones, lead to the following set of parameters T_eff = 6550 +- 100 K, log g = 3.75 +- 0.10, L_*/L_sun = 16.3 +- 4.5, M_*/M_sun = 2.0 +- 0.3 and an age of 5.0 +- 1.5 Myr. This stellar age provides further constrains to the mass of the possible companion estimated by Biller et al. (2012), being in-between 0.20 and 0.35 M_sun. Stellar accretion rates obtained from UV Balmer excess modelling, optical photospheric line veiling, and from the correlations with several emission lines spanning from the UV to the near-IR, are consistent to each other. The mean value from all previous tracers is 2 (+- 1) x 10^-7 M_sun yr^-1, which is within the upper limit gas flow rate from the outer to the inner disk recently provided by Cassasus et al. (2013). This suggests that almost all gas transferred between both components of the disk is not trapped by the possible planet(s) in-between but fall onto the central star, although it is discussed how the gap flow rate could be larger than previously suggested. In addition, we provide evidence showing that the stellar accretion rate of HD 142527 has increased by a factor ~ 7 on a timescale of 2-5 years.

An overview is presented of the recent advances in understanding the B[e] phenomenon among blue supergiant stars in light of high-angular resolution observations and with an emphasis on the results obtained by means of long baseline optical stellar interferometry. The focus of the review is on the circumstellar material and evolutionary phase of B[e] supergiants, but recent results on dust production in regular blue supergiants are also highlighted.

We initiated long-term optical interferometry monitoring of the diameters of unstable yellow hypergiants (YHG) with the goal of detecting both the long-term evolution of their radius and shorter term formation related to large mass-loss events. We observed HR5171 A with AMBER/VLTI. We also examined archival photometric data in the visual and near-IR spanning more than 60 years, as well as sparse spectroscopic data. HR5171A exhibits a complex appearance. Our AMBER data reveal a surprisingly large star for a YHG R*=1315+/-260Rsun\ (~6.1AU) at the distance of 3.6+/-0.5kpc. The source is surrounded by an extended nebulosity, and these data also show a large level of asymmetry in the brightness distribution of the system, which we attribute to a newly discovered companion star located in front of the primary star. The companion's signature is also detected in the visual photometry, which indicates an orbital period of Porb=1304+/-6d. Modeling the light curve with the NIGHTFALL program provides clear evidence that the system is a contact or possibly over-contact eclipsing binary. A total current system mass of 39^{+40}_{-22} solar mass and a high mass ratio q>10 is inferred for the system. The low-mass companion of HR5171 A is very close to the primary star that is embedded within its dense wind. Tight constraints on the inclination and vsini of the primary are lacking, which prevents us from determining its influence precisely on the mass-loss phenomenon, but the system is probably experiencing a wind Roche-Lobe overflow. Depending on the amount of angular momentum that can be transferred to the stellar envelope, HR5171 A may become a fast-rotating B[e]/Luminous Blue Variable (LBV)/Wolf-Rayet star. In any case, HR5171 A highlights the possible importance of binaries for interpreting the unstable YHGs and for massive star evolution in general.

Abridged: We have used the well-selected sample of ~1750 embedded, young, massive stars identified by the RMS survey to investigate the Galactic distribution of recent massive star formation. We describe the various methods used to assign distances extracted from the literature, and solve the distance ambiguities towards ~200 sources located within the Solar circle using archival HI data. These distances are used to calculate bolometric luminosities and estimate the survey completeness (~2x10^4 lsun). In total, we calculate the distance and luminosity of ~1650 sources, one third of which are above the survey's completeness threshold. Examination of the sample's longitude, latitude, radial velocities and mid-infrared images has identified ~120 small groups of sources, many of which are associated with well known star formation complexes, such as W43, W49 and W51. We compare the positional distribution of the sample with the expected locations of the spiral arms, assuming a model of the Galaxy consisting of four gaseous arms. The distribution of young massive stars in the Milky Way is spatially correlated with the spiral arms, with strong peaks in the source position and luminosity distributions at the arms' Galactocentric radii. The overall source and luminosity surface densities are both well correlated with the surface density of the molecular gas, which suggests that the massive star formation rate per unit molecular mass is approximately constant across the Galaxy. We estimate the total luminosity of the embedded massive star population to be ~0.76x10^8 lsun, 30% of which is associated with the ten most active star forming complexes. We measure the scale height as a function of Galactocentric distance and find that it increases only modestly from ~20-30 pc between 4 and 8 kpc, but much more rapidly at larger distances.

We present medium resolution ($R\approx5300$) $K$-band integral field spectroscopy of six MYSOs. The targets are selected from the RMS survey, and we used the NIFS on the Gemini North telescope. The data show various spectral line features including Br$\gamma$, CO, H$_2$, and \mbox{He\,{\sc i}}. The Br$\gamma$ line is detected in emission in all objects with $v_\mathrm{FWHM}\sim100$ -- 200 kms$^{-1}$. V645 Cyg shows a high-velocity P-Cygni profile between -800 kms$^{-1}$ and -300 kms$^{-1}$. We performed three-dimensional spectroastrometry to diagnose the circumstellar environment in the vicinity of the central stars using the Br$\gamma$ line. We measured the centroids of the velocity components with sub-mas precision. The centroids allow us to discriminate the blueshifted and redshifted components in a roughly east--west direction in both IRAS 18151--1208 and S106 in Br$\gamma$. This lies almost perpendicular to observed larger scale outflows. We conclude, given the widths of the lines and the orientation of the spectroastrometric signature, that our results trace a disc wind in both IRAS 18151--1208 and S106. The CO $\nu=2-0$ absorption lines at low $J$ transitions are detected in IRAS 18151--1208 and AFGL 2136. We analysed the velocity structure of the neutral gas discs. In IRAS 18151--1208, the absorption centroids of the blueshifted and redshifted components are separated in a direction of north-east to south-west, nearly perpendicular to that of the larger scale $H_2$ jet. The position-velocity relations of these objects can be reproduced with central masses of 30 $M_{\sun}$ for IRAS 18151--1208 and 20 $M_{\sun}$ for AFGL 2136. We also detect CO $\nu=2-0$ bandhead emission in IRAS 18151--1208, S106 and V645 Cyg. The results can be fitted reasonably with a Keplerian rotation model, with masses of 15, 20 and 20 $M_{\sun}$ respectively.

This work presents X-Shooter/VLT spectra of the prototypical, isolated Herbig Ae stars HD 31648 (MWC 480) and HD 163296 over five epochs separated by timescales ranging from days to months. Each spectrum spans over a wide wavelength range covering from 310 to 2475 nm. We have monitored the continuum excess in the Balmer region of the spectra and the luminosity of twelve ultraviolet, optical and near infrared spectral lines that are commonly used as accretion tracers for T Tauri stars. The observed strengths of the Balmer excesses have been reproduced from a magnetospheric accretion shock model, providing a mean mass accretion rate of 1.11 x 10^-7 and 4.50 x 10^-7 Msun yr^-1 for HD 31648 and HD 163296, respectively. Accretion rate variations are observed, being more pronounced for HD 31648 (up to 0.5 dex). However, from the comparison with previous results it is found that the accretion rate of HD 163296 has increased by more than 1 dex, on a timescale of ~ 15 years. Averaged accretion luminosities derived from the Balmer excess are consistent with the ones inferred from the empirical calibrations with the emission line luminosities, indicating that those can be extrapolated to HAe stars. In spite of that, the accretion rate variations do not generally coincide with those estimated from the line luminosities, suggesting that the empirical calibrations are not useful to accurately quantify accretion rate variability.

We present the Red MSX Source (RMS) Survey, the largest statistically selected catalog of young massive protostars and HII regions to date. We outline the construction of the catalog using mid and near infrared color selection, as well as the detailed follow up work at other wavelengths, and at higher spatial resolution in the infrared. We show that within the adopted selection bounds we are more than 90% complete for the massive protostellar population, with a positional accuracy of the exciting source of better than 2 arcseconds. We briefly summarize some of the results that can be obtained from studying the properties of the objects in the catalog as a whole, and find evidence that the most massive stars form: (i) preferentially nearer the Galactic centre than the anti-centre; (ii) in the most heavily reddened environments, suggestive of high accretion rates; and (iii) from the most massive cloud cores.

Disks are ubiquitous in stellar astronomy, and play a crucial role in the formation and evolution of stars. In this contribution we present an overview of the most recent results, with emphasis on high spatial and spectral resolution. We will start with a general discussion on direct versus indirect detection of disks, and then traverse the HR diagram starting with the pre-Main Sequence and ending with evolved stars.

We investigate the circumstellar dust properties of the oxygen-rich bipolar proto-planetary nebula IRAS 18276-1431 by means of two-dimensional radiative transfer simulations of the circumstellar dust shell. The model geometry is assumed to have a torus and an envelope. The parameters of the dust and the dust shell are constrained by comparing the SED and NIR intensity and polarisation data with the models. The polarisation in the envelope reaches 50 -- 60 % and is nearly constant in the H and K_S bands in the observations. This weak wavelength dependence of the polarisation can be reproduced with a grain size distribution function for the torus: 0.05 micron <= a with n(a)=a^{-(p=5.5)}exp(-a/{a_c=0.3 micron}). The power index p is significantly steeper than that for interstellar dust. Similar results have also been found in some other PPNs and suggest that mechanisms that grind down large particles may also have acted when the dust particles formed. The spectral opacity index beta is found to be 0.6+/-0.5 from the millimeter fluxes. This low value indicates the presence of large dust grains in the torus. We discuss two possible dust models for the torus. One has a size distribution function of 1.0 micron <= a <= a_max=5,000.0 micron with n(a)=a^{-(p=2.5)} and the other is 1.0 micron <= a <= a_max=10,000.0 micron with n(a)=a^{-(p=3.5)}. The former has beta of 0.633, but we are not able to find reasonable geometry parameters to fit the SED in the infrared. The latter has beta of 1.12, but reproduces the SED better over a wide wavelength range. With this dust model, the geometric parameters are estimated as follows: the inner and outer radii are 30 AU and 1000 AU and the torus mass is 3.0 M_sun. Assuming an expansion velocity of 15 kms^{-1}, the torus formation time and mass-loss rate are found to be \sim300 yrs and \sim10^{-2}M_sun yr^{-1} respectively.

Near-infrared H- and K-band spectra are presented for 247 objects, selected from the Red MSX Source (RMS) survey as potential young stellar objects (YSOs). 195 (~80%) of the targets are YSOs, of which 131 are massive YSOs (L_BOL > 5x10^3 L_solar), M > 8M_solar. This is the largest spectroscopic study of massive YSOs to date, providing a valuable resource for the study of massive star formation. In this paper we present our exploratory analysis of the data. The YSOs observed have a wide range of embeddedness (2.7 < A_V < 114), demonstrating that this study covers minimally obscured objects right through to very red, dusty sources. Almost all YSOs show some evidence for emission lines, though there is a wide variety of observed properties. The most commonly detected lines are Brgamma, H_2, fluorescent FeII, CO bandhead, [FeII] and HeI 2-1 2^1S-2^1P, in order of frequency of occurrence. In total, ~40% of the YSOs display either fluorescent FeII 1.6878um or CO bandhead emission (or both), indicative of a circumstellar disc; however, no correlation of the strength of these lines with bolometric luminosity was found. We also find that ~60% of the sources exhibit [FeII] or H_2 emission, indicating the presence of an outflow. Three quarters of all sources have Brgamma in emission. A good correlation with bolometric luminosity was observed for both the Brgamma and H_2 emission line strengths, covering 1 L_solar< L_BOL < 3.5x10^5 L_solar. This suggests that the emission mechanism for these lines is the same for low-, intermediate-, and high-mass YSOs, i.e. high-mass YSOs appear to resemble scaled-up versions of low-mass YSOs.

IRC +10420 is one of the few known massive stars in rapid transition from the Red Supergiant phase to the Wolf-Rayet or Luminous Blue Variable phase. The star has an ionised wind and using the Br gamma hydrogen recombination emission we assess the mass-loss on spatial scales of order 1 au. We present new VLT Interferometer AMBER data which are combined with all other AMBER data in the literature. The final dataset covers a position angle range of 180 degrees and baselines up to 110 meters. The spectrally dispersed visibilities, differential phases and line flux are conjointly analyzed and modelled. We also present AMBER/FINITO observations which cover a larger wavelength range and allow us to observe the Na I doublet at 2.2 micron. The data are complemented by X-Shooter data, which provide a higher spectral resolution view. The Brackett gamma line and the Na I doublet are both spatially resolved. After correcting the AMBER data for the fact that the lines are not spectrally resolved, we find that Br gamma traces a ring with a diameter of 4.18 milli-arcseconds. We consider a geometric model in which the Br gamma emission emerges from the top and bottom rings of an hour-glass shaped structure, viewed almost pole-on. It provides satisfactory fits to most visibilities and differential phases. The fact that we detect line emission from a neutral metal like Na I within the ionized region, a very unusual occurrence, suggests the presence of a dense pseudo-photosphere. The ionized wind can be reproduced with a polar wind, which could well have the shape of an hour-glass. The resolved Na I emission is found to occur on scales barely larger than the continuum. This fact and that many Yellow Hypergiants exhibit this comparatively rare emission hints at the presence of a "Yellow" or even "White Wall" in the Hertzsprung-Russell diagram, preventing them from visibly evolving to the blue.

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.

The CORNISH project is the highest resolution radio continuum survey of the Galactic plane to date. It is the 5 GHz radio continuum part of a series of multi-wavelength surveys that focus on the northern GLIMPSE region (10 deg < l < 65 deg), observed by the Spitzer satellite in the mid-infrared. Observations with the Very Large Array in B and BnA configurations have yielded a 1.5" resolution Stokes I map with a root-mean-squared noise level better than 0.4 mJy/beam. Here we describe the data-processing methods and data characteristics, and present a new, uniform catalogue of compact radio-emission. This includes an implementation of automatic deconvolution that provides much more reliable imaging than standard CLEANing. A rigorous investigation of the noise characteristics and reliability of source detection has been carried out. We show that the survey is optimised to detect emission on size scales up to 14" and for unresolved sources the catalogue is more than 90 percent complete at a flux density of 3.9 mJy. We have detected 3,062 sources above a 7-sigma detection limit and present their ensemble properties. The catalogue is highly reliable away from regions containing poorly-sampled extended emission, which comprise less than two percent of the survey area. Imaging problems have been mitigated by down-weighting the shortest spacings and potential artefacts flagged via a rigorous manual inspection with reference to the Spitzer infrared data. We present images of the most common source types found: regions, planetary nebulae and radio-galaxies. The CORNISH data and catalogue are available online at http://cornish.leeds.ac.uk

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.

We describe the motivation, design and implementation of the CORNISH survey, an arcsecond resolution radio continuum survey of the inner Galactic plane at 5 GHz using the Karl G. Jansky Very Large Array (VLA). It is a blind survey co-ordinated with the northern Spitzer GLIMPSE I region covering 10 deg < l < 65 deg and |b| < 1 deg at similar resolution. We discuss in detail the strategy that we employed to control the shape of the synthesised beam across this survey that covers a wide range of fairly low declinations. Two snapshots separated by 4 hours in hour angle kept the beam elongation to less that 1.5 over 75% of the survey area and less than 2 over 98% of the survey. The prime scientific motivation is to provide an unbiased survey for ultra-compact HII regions to study this key phase in massive star formation. A sensitivity around 2 mJy will allow the automatic distinction between radio loud and quiet mid-IR sources found in the Spitzer surveys. This survey has many legacy applications beyond star formation including evolved stars, active stars and binaries, and extragalactic sources. The CORNISH survey for compact ionized sources complements other Galactic plane surveys that target diffuse and non-thermal sources as well as atomic and molecular phases to build up a complete picture of the ISM in the Galaxy.

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 look for trends relating the mass accretion rate (Macc) and the stellar ages (t), spectral energy distributions (SEDs), and disk masses (Mdisk) for a sample of 38 HAeBe stars, comparing them to analogous correlations found for classical T Tauri stars. Our goal is to shed light on the timescale and physical processes that drive evolution of intermediate-mass pre-main sequence objects. Macc shows a dissipation timescale \tau = 1.3^{+1.0}_{-0.5} Myr from an exponential law fit, while a power law yields Macc(t) \propto t^{-\eta}, with \eta = 1.8^{+1.4}_{-0.7}. This result is based on our whole HAeBe sample (1-6 Msun), but the accretion rate decline most probably depends on smaller stellar mass bins. The near-IR excess is higher and starts at shorter wavelengths (J and H bands) for the strongest accretors. Active and passive disks are roughly divided by 2 x 10^{-7} Msun/yr. The mid-IR excess and the SED shape from the Meeus et al. classification are not correlated with Macc. We find Macc \propto Mdisk^{1.1 +- 0.3}. Most stars in our sample with signs of inner dust dissipation typically show accretion rates ten times lower and disk masses three times smaller than the remaining objects. The trends relating Macc with the near-IR excess and Mdisk extend those for T Tauri stars, and are consistent with viscous disk models. The differences in the inner gas dissipation timescale, and the relative position of the stars with signs of inner dust clearing in the Macc-Mdisk plane, could be suggesting a slightly faster evolution, and that a different process - such as photoevaporation - plays a more relevant role in dissipating disks in the HAeBe regime compared to T Tauri stars. Our conclusions must consider the mismatch between the disk mass estimates from mm fluxes and the disk mass estimates from accretion, which we also find in HAeBe stars.

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.

We present high-resolution HI data obtained using the Australia Telescope Compact Array to resolve the near/far distance ambiguities towards a sample of compact HII regions from the Red MSX Source (RMS) survey. The high resolution data are complemented with lower resolution archival HI data extracted from the Southern and VLA Galactic Plane surveys. We resolve the distance ambiguity for nearly all of the 105 sources where the continuum was strong enough to allow analysis of the HI absorption line structure. This represents another step in the determination of distances to the total RMS sample, which with over 1,000 massive young stellar objects and compact HII regions, is the largest and most complete sample of its kind. The full sample will allow the distribution of massive star formation in the Galaxy to be examined.

We performed a mid-infrared imaging survey of evolved stars in order to study the dust distribution in circumstellar envelopes around these objects and to better understand the mass-loss mechanism responsible for the formation of these envelopes. During this survey, we resolved for the first time the circumstellar environment of IRAS 17163-3907 (hereinafter IRAS17163), which is one of the brightest objects in the mid-infrared sky, but is surprisingly not well studied. Our aim is to determine the evolutionary status of IRAS 17163 and study its circumstellar environment in order to understand its mass-loss history. We obtained diffraction-limited images of IRAS 17163 in the mid-infrared using VISIR on the VLT. Optical spectra of the object allowed us to determine its spectral type and estimate its distance via the presence of diffuse interstellar bands. We show that IRAS 17163 is a Post-Red Supergiant, possibly belonging to the rare class of Yellow Hypergiants, and is very similar to the well studied object IRC +10420. Our mid-infrared images of IRAS 17163 are the first direct images of this bright mid-infrared source. These images clearly show the presence of a double dusty detached shell around the central star, due to successive ejections of material with a timescale of the order of 400 years and a total circumstellar mass larger than 4 M$_{\odot}$. This indicates that non quiescent mass-loss occurs during this phase of stellar evolution.

This work aims to derive accretion rates for a sample of 38 HAeBe stars. We apply magnetospheric accretion (MA) shock modelling to reproduce the observed Balmer excesses. We look for possible correlations with the strength of the Halpha, [OI]6300, and Brgamma emission lines. The median mass accretion rate is 2 x 10^-7 Msun yr^-1 in our sample. The model fails to reproduce the large Balmer excesses shown by the four hottest stars (T* > 12000 K). We derive Macc propto M*^5 and Lacc propto L*^1.2 for our sample, with scatter. Empirical calibrations relating the accretion and the Halpha, [OI]6300, and Brgamma luminosities are provided. The slopes in our expressions are slightly shallower than those for lower mass stars, but the difference is within the uncertainties, except for the [OI]6300 line. The Halpha 10% width is uncorrelated with Macc, unlike for the lower mass regime. The mean Halpha width shows higher values as the projected rotational velocities of HAe stars increase, which agrees with MA. The accretion rate variations in the sample are typically lower than 0.5 dex on timescales of days to months, Our data suggest that the changes in the Balmer excess are uncorrelated to the simultaneous changes of the line luminosities. The Balmer excesses and Halpha line widths of HAe stars can be interpreted within the context of MA, which is not the case for several HBes. The steep trend relating Macc and M* can be explained from the mass-age distribution characterizing HAeBe stars. The line luminosities used for low-mass objects are also valid to estimate typical accretion rates for the intermediate-mass regime under similar empirical expressions. However, we suggest that several of these calibrations are driven by the stellar luminosity.

The Red MSX Source (RMS) survey has identified a sample of ~1200 massive young stellar objects (MYSOs), compact and ultra compact HII regions from a sample of ~2000 MSX and 2MASS colour selected sources. We have used the 100 m Green Bank telescope to search for 22-24 GHz water maser and ammonia (1,1), (2,2) and (3,3) emission towards ~600 RMS sources located within the northern Galactic plane. We have identified 308 H2O masers which corresponds to an overall detection rate of ~50%. Abridged: We detect ammonia emission towards 479 of these massive young stars, which corresponds to ~80%. Ammonia is an excellent probe of high density gas allowing us to measure key parameters such as gas temperatures, opacities, and column densities, as well as providing an insight into the gas kinematics. The average kinetic temperature, FWHM line width and total NH3 column density for the sample are approximately 22 K, 2 km/s and 2x10^{15} cm^{-2}, respectively. We find that the NH3 (1,1) line width and kinetic temperature are correlated with luminosity and finding no underlying dependence of these parameters on the evolutionary phase of the embedded sources, we conclude that the observed trends in the derived parameters are more likely to be due to the energy output of the central source and/or the line width-clump mass relationship. The velocities of the peak H2O masers and the NH3 emission are in excellent agreement with each other, which would strongly suggest an association between the dense gas and the maser emission. Moreover, we find the bolometric luminosity of the embedded source and the isotropic luminosity of the H2O maser are also correlated. We conclude from the correlations of the cloud and water maser velocities and the bolometric and maser luminosity that there is a strong dynamical relationship between the embedded young massive star and the H2O maser.

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).

There is currently no accepted theoretical framework for the formation of the most massive stars, and the manner in which protostars continue to accrete and grow in mass beyond \sim10Msun is still a controversial topic. In this study we use several prescriptions of stellar accretion and a description of the Galactic gas distribution to simulate the luminosities and spatial distribution of massive protostellar population of the Galaxy. We then compare the observables of each simulation to the results of the Red MSX Source (RMS) survey, a recently compiled database of massive young stellar objects. We find that the observations are best matched by accretion rates which increase as the protostar grows in mass, such as those predicted by the turbulent core and competitive accretion (i.e. Bondi-Hoyle) models. These 'accelerating accretion' models provide very good qualitative and quantitative fits to the data, though we are unable to distinguish between these two models on our simulations alone. We rule out models with accretion rates which are constant with time, and those which are initially very high and which fall away with time, as these produce results which are quantitatively and/or qualitatively incompatible with the observations. To simultaneously match the low- and high-luminosity YSO distribution we require the inclusion of a 'swollen-star' pre-main-sequence phase, the length of which is well-described by the Kelvin-Helmholz timescale. Our results suggest that the lifetime of the YSO phase is \sim 10^5yrs, whereas the compact Hii-region phase lasts between \sim 2 - 4 \times 10^5yrs depending on the final mass of the star. Finally, the absolute numbers of YSOs are best matched by a globally averaged star-formation rate for the Galaxy of 1.5-2Msun/yr.

We analysed 337 multi-epoch optical spectra of 38 Herbig Ae/Be (HAeBe) stars to gain insights into the variability behaviour of the circumstellar (CS) atomic gas. Equivalent widths (EWs) and line fluxes of the Halpha, [OI]6300, HeI5876 and NaID lines were obtained for each spectrum; the Halpha line width at 10% of peak intensity (W10) and profile shapes were also measured and classified. The mean line strengths and relative variabilities were quantified for each star. Simultaneous optical photometry was used to estimate the line fluxes. We present a homogeneous spectroscopic database of HAeBe stars. The lines are variable in practically all stars and timescales, although 30 % of the objects show a constant EW in [OI]6300, which is also the only line that shows no variability on timescales of hours. The HeI5876 and NaID EW relative variabilities are typically the largest, followed by those in [OI]6300 and Halpha. The EW changes can be larger than one order of magnitude for the HeI5876 line, and up to a factor 4 for Halpha. The [OI]6300 and Halpha EW relative variabilities are correlated for most stars in the sample. The Halpha mean EW and W10 are uncorrelated, as are their relative variabilities. The Halpha profile changes in 70 % of the objects. The massive stars in the sample usually show more stable Halpha profiles with blueshifted self-absorptions and less variable 10% widths. Our data suggest multiple causes for the different line variations, but the [OI]6300 and Halpha variability must share a similar origin in many objects. The physical mechanism responsible for the Halpha line broadening does not depend on the amount of emission; unlike in lower-mass stars, physical properties based on the Halpha luminosity and W10 would significantly differ. Our results provide additional support to previous works that reported different physical mechanisms in Herbig Ae and Herbig Be stars.

We present a determination of the luminosity functions of massive young stellar objects (MYSOs) and compact (C)HII regions within the Milky Way Galaxy using the large, well-selected sample of these sources identified by the Red MSX Source (RMS) survey. The MYSO luminosity function decreases monotonically such that there are few with $L\gtrsim 10^{5}$Lsol, whilst the CHII regions are detected up to ~10$^{6}Lsol. The lifetimes of these phases are also calculated as a function of luminosity by comparison with the luminosity function for local main-sequence OB stars. These indicate that the MYSO phase has a duration ranging from 4x10$^{5}$ yrs for 10$^{4}$Lsol to ~7x10$^{4}$ yrs at 10$^{5}$Lsol, whilst the CHII region phase lasts of order 3x10$^{5}$ yrs or ~3-10% of the exciting star's main-sequence lifetime. MYSOs between 10$^{4} Lsol and ~10$^{5}$ Lsol are massive but do not display the radio continuum or near-IR \HI{} recombination line emission indicative of an HII region, consistent with being swollen due to high ongoing or recent accretion rates. Above ~10$^{5}$ Lsol the MYSO phase lifetime becomes comparable to the main-sequence Kelvin-Helmholtz timescale, at which point the central star can rapidly contract onto the main-sequence even if still accreting, and ionise a CHII region, thus explaining why few highly luminous MYSOs are observed.

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.

Context: The Red MSX Source (RMS) survey is returning a large sample of massive young stellar objects (MYSOs) and ultra-compact (UC) \HII{} regions using follow-up observations of colour-selected candidates from the MSX point source catalogue. Aims: To obtain the bolometric fluxes and, using kinematic distance information, the luminosities for young RMS sources with far-infrared fluxes. Methods: We use a model spectral energy distribution (SED) fitter to obtain the bolometric flux for our sources, given flux data from our work and the literature. The inputs to the model fitter were optimised by a series of investigations designed to reveal the effect varying these inputs had on the resulting bolometric flux. Kinematic distances derived from molecular line observations were then used to calculate the luminosity of each source. Results: Bolometric fluxes are obtained for 1173 young RMS sources, of which 1069 have uniquely constrained kinematic distances and good SED fits. A comparison of the bolometric fluxes obtained using SED fitting with trapezium rule integration and two component greybody fits was also undertaken, and showed that both produce considerable scatter compared to the method used here. Conclusions: The bolometric flux results allowed us to obtain the luminosity distributions of YSOs and UC\HII{} regions in the RMS sample, which we find to be different. We also find that there are few MYSOs with L $\geq$ 10$^{5}$\lsol{}, despite finding many MYSOs with 10$^{4}$\lsol{} $\geq$ L $\geq$ 10$^{5}$\lsol{}.

The Red MSX Source (RMS) survey has identified a large sample of massive young stellar objects (MYSOs) and ultra compact (UC) HII regions from a sample of ~2000 MSX and 2MASS colour selected sources. Using a recent catalogue of molecular clouds derived from the Boston University-Five College Radio Astronomy Observatory Galactic Ring Survey (GRS), and by applying a Galactic scaleheight cut off of 120 pc, we solve the distance ambiguity for RMS sources located within 18\degr < |l| > 54\degr. These two steps yield kinematic distances to 291 sources out of a possible 326 located within the GRS longitude range. Combining distances and integrated fluxes derived from spectral energy distributions, we estimate luminosities to these sources and find that > 90% are indicative of the presence of a massive star. We find the completeness limit of our sample is ~10^4 Lsun, which corresponds to a zero age main sequence (ZAMS) star with a mass of ~12 Msun. Selecting only these sources, we construct a complete sample of 196 sources. Comparing the properties of the sample of young massive stars with the general population, we find the RMS-clouds are generally larger, more massive, and more turbulent. We examine the distribution of this sub-sample with respect to the location of the spiral arms and the Galactic bar and find them to be spatially correlated. We identify three significant peaks in the source surface density at Galactocentric radii of approximately 4, 6 and 8 kpc, which correspond to the proposed positions of the Scutum, Sagittarius and Perseus spiral arms, respectively. Fitting a scale height to the data we obtain an average value of ~29+-0.5 pc, which agrees well with other reported values in the literature, however, we note a dependence of the scale height on galactocentric radius with it increases from 30 pc to 45 pc between 2.5 and 8.5 kpc.

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.

In order to understand the formation mechanism of the disks around Be stars it is imperative to have a good overview of both the differences and similarities between normal B stars and the Be stars. Here we investigate a previous report that there may be a large population of sub-arcsecond companions to Be stars. We present the first systematic, comparative imaging study of the binary properties of matched samples of B and Be stars observed using the same equipment. We obtained high angular resolution (0.07-0.1 arcsec) K band Adaptive Optics data of 40 B stars and 39 Be stars. The separations that can be probed range from 0.1 to 8 arcsec (corresponding to 20-1000 AU), and magnitude differences up to 10 magnitudes can in principle be covered. We detect 11 binaries out of 37 Be targets (corresponding to a binary fraction of 30 +/- 8%) and 10 binaries out of 36 B targets (29 +/- 8%). Further tests demonstrate that the B and Be binary systems are not only similar in frequency but also remarkably similar in terms of binary separations, flux differences and mass ratios. We find that any hypotheses invoking binary companions as responsible for the formation of a disk need the companions to be closer than 20 AU. Close companions are known to affect the circumstellar disks of Be stars, but as not all Be stars have been found to be close binaries, the data suggest that binarity can not be responsible for the Be phenomenon in all Be stars. Finally, the similarities of the binary parameters themselves also shed light on the Be formation mechanism. They strongly suggest that the initial conditions that gave rise to B and Be star formation must, to all intents and purposes, be similar. This in turn indicates that the Be phenomenon is not the result of a different star formation mechanism.

Context: The Red MSX Source (RMS) survey is a multi-wavelength campaign of follow-up observations of a colour-selected sample of candidate massive young stellar objects (MYSOs) in the galactic plane. This survey is returning the largest well-selected sample of MYSOs to date, while identifying other dust contaminant sources with similar mid-infrared colours including a large number of new ultra-compact (UC)HII regions. Aims:To measure the far-infrared (IR) flux, which lies near the peak of the spectral energy distribution (SED) of MYSOs and UCHII regions, so that, together with distance information, the luminosity of these sources can be obtained. Methods:Less than 50% of RMS sources are associated with IRAS point sources with detections at 60 micron and 100 micron, though the vast majority are visible in Spitzer MIPSGAL or IRAS Galaxy Atlas (IGA) images. However, standard aperture photometry is not appropriate for these data due to crowding of sources and strong spatially variable far-IR background emission in the galactic plane. A new technique using a 2-dimensional fit to the background in an annulus around each source is therefore used to obtain far-IR photometry for young RMS sources. Results:Far-IR fluxes are obtained for a total of 1113 RMS candidates identified as young sources. Of these 734 have flux measurements using IGA 60 micron and 100 micron images and 724 using MIPSGAL 70 micron images, with 345 having measurements in both data sets.

The circumstellar structure on 100 AU scales of the massive young stellar object W33A is probed using the VLTI and the MIDI instrument. N-band visibilities on 4 baselines are presented which are inconsistent with a spherically symmetric geometry. The visibility spectra and SED are simultaneously compared to 2D axi-symmetric dust radiative transfer models with a geometry including a rotationally flattened envelope and outflow cavities. We assume an O7.5 ZAMS star as the central source, consistent with the observed bolometric luminosity. The observations are also compared to models with and without (dusty and gaseous) accretion disks. A satisfactory model is constructed which reproduces the visibility spectra for each (u,v) point. It fits the silicate absorption, the mid-IR slope, the far-infrared peak, and the (sub)mm of the SED. It produces a 350 micron morphology consistent with observations. The 10 micron emission on 100 AU scales is dominated by the irradiated walls of the cavity sculpted by the outflow. The visibilities rule out the presence of dust disks with total (gas and dust) masses more than 0.01 Msun. However, optically thick accretion disks, interior to the dust sublimation radius, are allowed to accrete at rates equalling the envelope's mass infall rate (up to 10^(-3) Msun/yr) without substantially affecting the visibilities due to the extinction by the extremely massive envelope of W33A.

The Young Stellar Object (YSO) W33A is one of the best known examples of a massive star still in the process of forming. Here we present Gemini North ALTAIR/NIFS laser-guide star adaptive-optics assisted K-band integral-field spectroscopy of W33A and its inner reflection nebula. In our data we make the first detections of a rotationally-flattened outer envelope and fast bi-polar jet of a massive YSO at near-infrared wavelengths. The predominant spectral features observed are Br-gamma, H_2, and a combination of emission and absorption from CO gas. We perform a 3-D spectro-astrometric analysis of the line emission, the first study of its kind. We find that the object's Br-gamma emission reveals evidence for a fast bi-polar jet on sub-milliarcsecond scales, which is aligned with the larger-scale outflow. The hybrid CO features can be explained as a combination of hot CO emission arising in a disk close to the central star, while cold CO absorption originates in the cooler outer envelope. Kinematic analysis of these features reveals that both structures are rotating, and consistent with being aligned perpendicularly to both the ionised jet and the large-scale outflow. Assuming Keplerian rotation, we find that the circumstellar disk orbits a central mass of >10Msun, while the outer envelope encloses a mass of ~15Msun. Our results suggest a scenario of a central star accreting material from a circumstellar disk at the centre of a cool extended rotating torus, while driving a fast bi-polar wind. These results therefore provide strong supporting evidence for the hypothesis that the formation mechanism for high-mass stars is qualitatively similar to that of low-mass stars.

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.

Context: The Red MSX Source (RMS) survey has identified a large sample of candidate massive young stellar objects (MYSOs) and ultra compact (UC) HII regions from a sample of ~2000 MSX and 2MASS colour selected sources. Aims: To search for H2O masers towards a large sample of young high mass stars and to investigate the statistical correlation of H2O masers with the earliest stages of massive star formation. Methods: We have used the Mopra Radio telescope to make position-switched observations towards ~500 UCHII regions and MYSOs candidates identified from the RMS survey and located between 190\degr < l < 30\degr. These observations have a 4$\sigma$ sensitivity of ~1 Jy and a velocity resolution of ~0.4 km/s.} Results: We have detected 163 H2O masers, approximately 75% of which were previously unknown. Comparing the maser velocities with the velocities of the RMS sources, determined from 13CO observations, we have identified 135 RMS-H2O maser associations, which corresponds to a detection rate of ~27%. Taking into account the differences in sensitivity and source selection we find our detection rate is in general agreement with previously reported surveys. Conclusions: We find similar detection rates for UCHII regions and MYSOs candidates, suggesting that the conditions needed for maser activity are equally likely in these two stages of the star formation process. Looking at the detection rate as a function of distance from the Galactic centre we find it significantly enhanced within the solar circle, peaking at ~37% between 6-7 kpc, which is consistent with previous surveys of UC HII regions, possibly indicating the presence of a high proportion of more luminous YSOs and HII regions.

(abridged) As the favoured progenitors of long-duration gamma-ray bursts, massive stars may represent our best signposts of individual objects in the early Universe, but special conditions seem required to make these bursters, which might originate from the progenitor's rapid rotation and associated asymmetry. To obtain empirical constraints on the interplay between stellar rotation and wind asymmetry, we perform linear Halpha spectropolarimetry on a sample of 18 spectroscopically peculiar massive O stars, including OVz, Of?p, Oe, and Onfp stars, supplemented by an earlier sample of 20 O supergiants. Despite their rapid rotation (with vsin(i) up to 400 km/s) most O-type stars are found to be spherically symmetric, but with notable exceptions amongst specific object classes. We divide the peculiar O stars into four distinct categories: Groups III and IV include the Oe stars and Onfp stars, which are on the high-end tail of the O star rotation distribution and have in the past been claimed to be embedded in disks. Here we report the detection of a classical depolarization ``line effect'' in the Oe star HD 45314, but the overall incidence of line effects amongst Oe stars is significantly lower (1 out of 6) than amongst Be stars. The chance that the Oe and Be datasets are drawn from the same parent population is negligible (with 95% confidence). This implies there is as yet no evidence for a disk hypothesis in Oe stars, providing relevant constraints on the physical mechanism that is responsible for the Be phenomenon. Finally, we find that 3 out of 4 of the group IV Onfp stars show evidence for complex polarization effects, which are likely related to rapid rotation, and we speculate on the evolutionary links to B[e] stars.

(Abridged) Context: The Red MSX Source (RMS) survey is an ongoing multi-wavelength observational programme designed to return a large, well-selected sample of massive young stellar objects (MYSOs). We have identified $\sim$2000 MYSO candidates located throughout the Galaxy by comparing the colours of MSX and 2MASS point sources to those of known MYSOs. Aims: To identify the populations of UCHII regions and PNe within the sample and examine their Galactic distribution. Method: We have conducted high resolution radio continuum observations at 6 cm towards 659 MYSO candidates in the northern hemisphere ($10\degr< l < 250\degr$) using the VLA. In addition to these targeted observations we present archival data towards a further 315 RMS sources extracted from a previous VLA survey of the inner Galaxy. Results: We find radio emission towards 272 ($\sim$27% of the observed sample). Using results from other parts of our multi-wavelength survey we separate these RMS-radio associations into two distinct types of objects, classifying 51 as PNe and a further 208 as either compact or UC HII regions. Using this well selected sample of HII regions we estimate their Galactic scale height to be 0.6\degr. Conclusions: Using radio continuum and archival data we have identified 79 PNe and 391 HII regions within the northern RMS catalogue. We estimate the total fraction of contamination by PNe in the RMS sample is of order 10%. The sample of HII regions is probably the best representation to date of the Galactic population of HII regions as a whole.

We discuss VLTI AMBER and MIDI interferometry in addition to single-dish Subaru observations of massive young stellar objects. The observations probe linear size scales between 10 to 1000 AU for the average distance of our sources.

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.

Massive young stellar objects (MYSO) are surrounded by massive dusty envelopes. Our aim is to establish their density structure on scales of ~1000 AU, i.e. a factor 10 increase in angular resolution compared to similar studies performed in the (sub)mm. We have obtained diffraction-limited (0.6") 24.5 micron images of 14 well-known massive star formation regions with Subaru/COMICS. The images reveal the presence of discrete MYSO sources which are resolved on arcsecond scales. For many sources, radiative transfer models are capable of satisfactorily reproducing the observations. They are described by density powerlaw distributions (n(r) ~ r^(-p)) with p = 1.0 +/-0.25. Such distributions are shallower than those found on larger scales probed with single-dish (sub)mm studies. Other sources have density laws that are shallower/steeper than p = 1.0 and there is evidence that these MYSOs are viewed near edge-on or near face-on, respectively. The images also reveal a diffuse component tracing somewhat larger scale structures, particularly visible in the regions S140, AFGL 2136, IRAS 20126+4104, Mon R2, and Cep A. We thus find a flattening of the MYSO envelope density law going from ~10 000 AU down to scales of ~1000 AU. We propose that this may be evidence of rotational support of the envelope (abridged).

We present diffraction limited (0.6") 24.5micron Subaru/COMICS images of the red supergiant mu Cep. We report the detection of a circumstellar nebula, that was not detected at shorter wavelengths. It extends to a radius of at least 6" in the thermal infrared. On these angular scales, the nebula is roughly spherical, in contrast, it displays a pronounced asymmetric morphology closer in. We simultaneously model the azimuthally averaged intensity profile of the nebula and the observed spectral energy distribution using spherical dust radiative transfer models. The models indicate a constant mass-loss process over the past 1000 years, for mass-loss rates a few times 10^(-7) Msun/yr. This work supports the idea that at least part of the asymmetries in shells of evolved massive stars and supernovae may be due to the mass-loss process in the red supergiant phase.

We present interferometric and single-dish mid-infrared observations of a sample of massive young stellar objects (BN-type objects), using VLTI-MIDI (10 micron) and Subaru-COMICS (24.5 micron). We discuss the regions S140, Mon R2, M8E-IR, and W33A. The observations probe the inner regions of the dusty envelope at scales of 50 milli arcsecond and 0.6 arcsec (100-1000 AU), respectively. Simultaneous model fits to spectral energy distributions and spatial data are achieved using self-consistent spherical envelope modelling. We conclude that those MYSO envelopes that are best described by a spherical geometry, the commensurate density distribution is a powerlaw with index -1.0. Such a powerlaw is predicted if the envelope is supported by turbulence on the 100-1000AU scales probed with MIDI and COMICS, but the role of rotation at these spatial scales need testing.

The origin of the disks around Be stars is still not known. Further progress requires a proper parametrization of their structure, both spatially and kinematically. This is challenging as the disks are very small. Here we assess whether a novel method is capable of providing these data. We obtained spectro astrometry around the Pa beta line of two bright Be stars, alpha Col and zeta Tau, to search for disk signatures. The data, with a pixel to pixel precision of the centroid position of 0.3..0.4 milliarcsecond is the most accurate such data to date. Artefacts at the 0.85 mas level are present in the data, but these are readily identified as they were non-repeatable in our redundant datasets. This does illustrate the need of taking multiple data to avoid spurious detections. The data are compared with simple model simulations of the spectro astrometric signatures due to rotating disks around Be stars. The upper limits we find for the disk radii correspond to disk sizes of a few dozen stellar radii if they rotate Keplerian. This is very close to observationally measured and theoretically expected disk sizes, and this paper therefore demonstrates that spectro-astrometry, of which we present the first such attempt, has the potential to resolve the disks around Be stars.

Context: The Red MSX Source (RMS) survey is an ongoing multi-wavelength observational programme designed to return a large, high-resolution mid-infrared colour-selected sample of massive young stellar objects. Aims: A critical part of our follow-up programme is to conduct 13CO molecular line observations in order to determine kinematic distances to all of our MYSO candidates. These distances will allow us to identify and remove nearby low-mass YSOs and help in identifying evolved stars which are weak CO emitters. Method: We have used the 15 m James Clerk Maxwell Telescope (JCMT), the 13.7 m telescope of the Purple Mountain Observatory (PMO), the 20 m Onsala telescope and the 22m Mopra telescope to conduct molecular line observations towards 508 MYSOs candidates located in the 1st and 2nd Quadrants. Results: We detect 13CO emission towards 780 RMS sources which corresponds to approximately 84% of those observed (911). A total of 2595 emission components are detected above 3sigma level (typically T^*_{\rm{A}} > 0.3K), with multiple components being observed towards the majority of these sources -- 520 sources (~56%) -- with an average of ~4 molecular clouds detected along each line of sight. We have used archival CS (J=2-1) and maser velocities to resolve the component multiplicity towards 175 sources (~20%) and have derived a criterion which is used to identify the most likely component for a further 191 multiple component sources. Combined with the single component detections we have obtained unambiguous kinematic velocities for 638 of the 780 MYSOs candidates towards which CO is detected (~80% of the detections). Using the Galactic rotation curve we calculate kinematic distances for all detected components.

The yellow hypergiants are found in a stage between the massive Red Supergiants and the Wolf-Rayet stars. This review addresses current issues concerning the evolution of massive stars, concentrating on the transitional post-Red Supergiant phase. Few yellow hypergiants are known and even fewer show direct evidence for having evolved off the Red Supergiant branch. Indeed, only two such rare objects with clear evidence for having gone through of a previous mass losing phase are known, IRC +10420 and HD 179821. We will review their properties and present recent results employing near-infrared interferometry, integral field spectroscopy and polarimetry. Finally, their real-time evolution is discussed.

We present medium resolution spectropolarimetry and long term photo-polarimetry of two massive post-red supergiants, IRC +10420 and HD 179821. The data provide new information on their circumstellar material as well as their evolution. In IRC +10420, the polarization of the Halpha line is different to that of the continuum, which indicates that the electron-scattering region is not spherically symmetric. The observed long term changes in the polarimetry can be associated with an axi-symmetric structure, along the short axis of the extended reflection nebulosity. Long term photometry reveals that the star increased in temperature until the mid-nineties, after which the photospheric flux in the optical levelled off. As the photometric changes are mostly probed in the red, they do not trace high stellar temperatures sensitively. And so, it is not obvious whether the star has halted its increase in temperature or not. For HD 179821 we find no polarization effects across any absorption or emission lines, but observe very large polarization changes of order 5% over 15 years. Unexpectedly, during the same period, the optical photometry displayed modest variability at the 0.2 magnitude level. Several explanations for this puzzling fact are discussed. Most of which, involving asymmetries in the circumstellar material, seem to fail as there is no evidence for the presence of hot, dusty material close to the star. Alternatively, the variations can be explained by the presence of a non-radially pulsating photosphere. Changes in the photometry hint at an increase in temperature corresponding to a change through two spectral subclasses over the past ten years.

The paper investigates the milli-arcsecond scale structure of the present-day mass-loss of the post-Red Supergiant IRC+10420. We use three telescopes of the VLT Interferometer in combination with the AMBER near-infrared beam combiner to measure spectrally dispersed correlated fluxes in the K-band around the Br gamma transition. The resulting visibilities are compared to the predicted visibilities of emission structures with various simple models in order to infer the size of the observed emission region. The Br gamma line is resolved by VLTI+AMBER on all three baselines, with the maximum projected baseline extending 69 meter and a P.A. ranging between 10 and 30 degrees. A differential phase between line and continuum is detected on the longest baseline. The Br gamma emission region is found to have a diameter of 3.3 milli-arcseconds (FWHM), when compared to a Gaussian intensity distribution. A uniform disk and a ring-like intensity distribution do not fit the line visibilities. Comparing the AMBER equivalent width of Br gamma with measurements from various epochs, we find that the stellar photosphere contributes about 60% of the total continuum light at 2.2 micron. The remaining 40% continuum emission is found on scales larger than the 66mas AMBER field of view. Using simple arguments, and assuming optically thick line emission, we find that the line emitting region is elongated. We briefly discuss the possibilities whether such a structure is due to a bi-polar flow or a circumstellar disk. (Abridged).

Here we describe the Red MSX Source (RMS) survey which is the largest, systematic, galaxy-wide search for massive young stellar objects (MYSOs) yet undertaken. Mid-IR bright point sources from the MSX satellite survey have been followed-up with ground-based radio, millimetre, and infrared observations to identify the contaminating sources and characterise the MYSOs and UCHII regions. With the initial classification now complete the distribution of sources in the galaxy will be discussed, as well as some programmes being developed to exploit our sample.

We report on resolved interferometric observations with VLTI/MIDI of the massive young stellar object (MYSO) W33A. The MIDI observations deliver spectrally dispersed visibilities with values between 0.03 and 0.06, for a baseline of 45m over the wavelength range 8-13 micron. The visibilities indicate that W33A has a FWHM size of approximately 120AU (0.030'') at 8 micron which increases to 240AU at 13 micron, scales previously unexplored among MYSOs. This observed trend is consistent with the temperature falling off with distance. 1D dust radiative transfer models are simultaneously fit to the visibility spectrum, the strong silicate feature and the shape of the mid infrared spectral energy distribution (SED). For any powerlaw density distribution, we find that the sizes (as implied by the visibilities) and the stellar luminosity are incompatible. A reduction to a third of W33A's previously adopted luminosity is required to match the visibilities; such a reduction is consistent with new high resolution 70 micron data from Spitzer's MIPSGAL survey. We obtain best fits for models with shallow dust density distributions of r^(-0.5) and r^(-1.0) and for increased optical depth in the silicate feature produced by decreasing the ISM ratio of graphite to silicates and using optical grain properties by Ossenkopf et al. (1992).

Abridged abstract: The Red MSX Source (RMS) survey is an ongoing effort to return a large, well-selected sample of massive young stellar objects (MYSOs) within our Galaxy. A series of ground-based follow-up observations are being undertaken in order to remove contaminant objects from our list of 2000 candidates, and to begin characterising these MYSOs. As a part of these follow-up observations, high resolution (~1") mid-IR imaging aids the identification of contaminant objects which are resolved (UCHII regions, PN) as opposed to those which are unresolved (YSOs, evolved stars) as well as identifying YSOs near UCHII regions and other multiple sources. We present 10.4 micron imaging observations for 346 candidate MYSOs in the RMS survey in the Southern Hemisphere, primarily outside the region covered by the GLIMPSE Spitzer Legacy Survey. These were obtained using TIMMI2 on the ESO 3.6m telescope in La Silla, Chile. Our photometric accuracy is of order 0.05Jy, and our astrometric accuracy is 0.8", which is an improvement over the nominal 2" accuracy of the MSX PSC.

We present NAOMI/OASIS adaptive-optics assisted integral-field spectroscopy of the transitional massive hypergiant IRC +10420, an extreme mass-losing star apparently in the process of evolving from a Red Supergiant toward the Wolf-Rayet phase. To investigate the present-day mass-loss geometry of the star, we study the appearance of the line-emission from the inner wind as viewed when reflected off the surrounding nebula. We find that, contrary to previous work, there is strong evidence for wind axi-symmetry, based on the equivalent-width and velocity variations of H$\alpha$ and Fe {\sc ii} $\lambda$6516. We attribute this behaviour to the appearance of the complex line-profiles when viewed from different angles. We also speculate that the Ti {\sc ii} emission originates in the outer nebula in a region analogous to the Strontium Filament of $\eta$ Carinae, based on the morphology of the line-emission. Finally, we suggest that the present-day axisymmetric wind of IRC +10420, combined with its continued blueward evolution, is evidence that the star is evolving toward the B[e] supergiant phase.

We present the results of Monte Carlo mass-loss predictions for massive stars covering a wide range of stellar parameters. We critically test our predictions against a range of observed mass-loss rates -- in light of the recent discussions on wind clumping. We also present a model to compute the clumping-induced polarimetric variability of hot stars and we compare this with observations of Luminous Blue Variables, for which polarimetric variability is larger than for O and Wolf-Rayet stars. Luminous Blue Variables comprise an ideal testbed for studies of wind clumping and wind geometry, as well as for wind strength calculations, and we propose they may be direct supernova progenitors.

Abridged: The Red MSX Source (RMS) survey is an ongoing multi-wavelength observational programme designed to return a large, well-selected sample of massive young stellar objects (MYSOs). Here we present 13CO observations made towards 854 MYSOs candidates located in the 3rd and 4th quadrants. We detected 13CO emission towards a total of 751 of the 854 RMS sources observed (~88%). In total 2185 emission components are detected above 3$\sigma$ level. Multiple emission profiles are observed towards the majority of these sources - 455 sources (~60%) - with an average of ~4 molecular clouds along the line of sight. These multiple emission features make it difficult to assign a kinematic velocity to many of our sample. We have used archival CS (J=2-1) and maser velocities to resolved the component multiplicity towards 82 sources and have derived a criterion which is used to identify the most likely component for a further 202 multiple component sources. Combined with the single component detections we have obtained unambiguous kinematic velocities towards 580 sources (~80% of the detections). The 171 sources for which we have not been able to determine the kinematic velocity will require additional line data. Using the rotation curve of Brand and Blitz (1993) and their radial velocities we calculate kinematic distances for all components detected.

Clumping in the winds of massive stars may significantly reduce empirical mass-loss rates, and which in turn may have a large impact on our understanding of massive star evolution. Here, we investigate wind-clumping through the linear polarization induced by light scattering off the clumps. Through the use of an analytic wind clumping model, we predict the time evolution of the linear polarimetry over a large parameter space. We concentrate on the Luminous Blue Variables, which display the greatest amount of polarimetric variability and for which we recently conducted a spectropolarimetric survey. Our model results indicate that the observed level of polarimetric variability can be reproduced for two regimes of parameter space: one of a small number of massive, optically-thick clumps; and one of a very large number of low-mass clumps. Although a systematic time-resolved monitoring campaign is required to distinguish between the two scenarios, we currently favour the latter, given the short timescale of the observed polarization variability. As the polarization is predicted to scale linearly with mass-loss rate, we anticipate that all hot stars with very large mass-loss rates should display polarimetric variability. This is consistent with recent findings that intrinsic polarization is more common in stars with strong H$\alpha$ emission.

We present linear spectropolarimetric data for eight Herbig Be and four Herbig Ae stars at H alpha, H beta and H gamma. Changes in the linear polarisation are detected across all Balmer lines for a large fraction of the observed objects, confirming that the small-scale regions surrounding these objects are flattened (i.e. disk-like). Furthermore, all objects with detections show similar characteristics at the three spectral lines, despite differences in transition probability and optical depth going from H alpha to H gamma. A large fraction of early Herbig Be stars (B0-B3) observed show line depolarisation effects. However the early Herbig Ae stars (A0-A2), observed for comparison, show intrinsic line polarisation signatures. Our data suggest that the popular magnetic accretion scenario for T Tauri objects may be extended to Herbig Ae stars, but that it may not be extended to early Herbig Be stars, for which the available data are consistent with disc accretion.

The Red MSX Source (RMS) survey is a multi-wavelength programme of follow-up observations designed to distinguish between genuine massive young stellar objects (MYSOs) and other embedded or dusty objects, such as ultra compact (UC) HII regions, evolved stars and planetary nebulae (PNe), from a sample of ~2000 MYSOs candidates. These were identified by comparing the colours of objects from the MSX and 2MASS point source catalogues to those of known MYSOs, in order to develop colour selection criteria which have been used to produced the RMS sample of MYSOs candidates. Our ultimate aim is to produce a large unbiased sample of MYSOs (~500) with complementary multi-wavelength data with which to study their properties. Here we report the results of 826 radio continuum and 346 mid-IR imaging observations carried out using the Australia Telescope Compact Array and TIMMI2 on the ESO 3.6 m telescope respectively. These observations are aimed at identifying and removing contaminating sources.

Circumstellar disks play an important role in many stages of the evolution of stars. However, it is only possible to directly image circumstellar disks for a few of the nearest stars. For massive stars, the situation is even more difficult, as they are on average further away than the more numerous low mass stars. Here, we discuss spectropolarimetry which is a technique that can reveal the presence of disks on spatial scales of order stellar radii. It is therefore extremely useful in studying disks around young stars, while the method has prospects for quantitative modelling. In this review we will discuss the status on the study of (accretion) disks around intermediate- to high mass pre-main sequence stars and will concentrate on the linear spectropolarimetry aspect. We wil make a brief excursion to observations of massive evolved stars where spectropolarimetry has revealed the clumpy structure of stellar winds, and conclude with a brief outlook.

Be stars, which are characterised by intermittent emission in their hydrogen lines, are known to be fast rotators. This fast rotation is a requirement for the formation of a Keplerian disk, which in turn gives rise to the emission. However, the pulsating, magnetic B1IV star $\beta$ Cephei is a very slow rotator that still shows H$\alpha$ emission episodes like in other Be stars, contradicting current theories. We investigate the hypothesis that the H$\alpha$ emission stems from the spectroscopically unresolved companion of $\beta$ Cep. Spectra of the two unresolved components have been separated in the 6350-6850\AA range with spectro-astrometric techniques, using 11 longslit spectra obtained with ALFOSC at the Nordic Optical Telescope, La Palma. We find that the H$\alpha$ emission is not related to the primary in $\beta$ Cep, but is due to its 3.4 magnitudes fainter companion. This companion has been resolved by speckle techniques, but it remains unresolved by traditional spectroscopy. The emission extends from about $-$400 to +400 km s$^{-1}$. The companion star in its 90-year orbit is likely to be a classical Be star with a spectral type around B6-8. By identifying its Be-star companion as the origin of the H$\alpha$ emission behaviour, the enigma behind the Be status of the slow rotator $\beta$ Cep has been resolved.

We present new medium-resolution spectropolarimetric observations of HD 45677 in the B and R-bands. A change in polarisation is detected across Halpha, Hbeta and Hgamma confirming that the ionised region around the star is aspherical. (Q,U) points associated with these emission lines occur away from the continuum, defining a polarisation vector which points in the same direction for each of the lines at an average intrinsic polarisation angle of 164 +/- 3 degrees. These data were combined with past photometric and polarimetric data from the literature to investigate any variability. We find that HD 45677 is both photometrically and polarimetrically variable and that these changes are linked. We suggest that these variations may be caused by an aspherical blowout and by deriving a least-squares fit to the B-band polarimetric data in Q-U space, we show that the blowout occurs at an intrinsic polarisation angle of 175 +/- 1 degrees, along the same angle as the proposed geometry of the ionised region.

The Red MSX Source (RMS) survey is a multi-wavelength program of follow-up observations designed to distinguish between genuine massive young stellar objects (MYSOs) and other embedded or dusty objects, such as ultra compact (UC) HII regions, evolved stars and planetary nebulae (PNe). We have identified nearly 2000 massive YSOs candidates by comparing the colours of MSX and 2MASS point sources to those of known MYSOs. Unfortunately, there are several other types of embedded or dust enshrouded objects that have similar colours as MYSOs and contaminate our sample. Two sources of contamination are from UCHII regions and PNe, both of which can be identified from the radio emission emitted by their ionised nebulae. In order to identify UCHII regions and PNe that contaminate our sample we have conducted high resolution radio continuum observations at 3.6 and 6 cm of all southern MYSOs candidates ($235< l < 350$) using the Australia Telescope Compact Array (ATCA). These observations have a spatial resolution of $\sim$1--2\arcsec and typical image rms noise values of $\sim$0.3 mJy -- sensitive enough to detect a B0.5 star at the far side of the Galaxy. Of the 826 RMS sources observed we found 199 to be associated with radio emission, $\sim$25% of the sample. More significant than those sources associated with radio emission are the 627 RMS sources for which no radio emission was detected as these sources are still potential MYSOs. The Galactic distribution, morphologies and spectral indices of the radio sources associated with the RMS sources are consistent with these sources being UCHII regions. In addition to the 826 RMS sources observed we present observations of a further 190 fields. These observations were made towards MSX sources that passed cuts in earlier versions of the survey, but were later excluded.

The topic of wind-clumping has been the subject of much activity in recent years, due to the impact that it can have on derived mass-loss rates. Here we present an alternative method of investigating wind-clumping, that of polarimetry. We present simulations of the polarization produced by a clumpy wind, and argue that the observations may be reproduced just by statistical deviations from spherical symmetry when the outflow is only slightly fragmented. Here, the polarization scales with $\dot{M}$, which is consistent with observations of LBVs, WRs and O supergiants. Finally, we find clumping factors in the inner 2$R_{\star}$ of $\sim 2-3$, and speculate as to the clumping stratification of hot stars.

As part of the Red MSX Source Survey of Massive Young Stellar Objects (MYSOs) we have conducted multi-wavelength follow up observations of the well-known object V645 Cygni. We present our data on this object, whose near-infrared spectrum is exceptional and place these in context with previous observations. Our observations of V645 Cyg included near/mid infrared imaging observations, 13CO 2-1 line observations and high signal-to-noise velocity resolved near-infrared spectroscopy. The spectrum shows P-Cygni hydrogen Brackett emission, consistent with a high velocity stellar wind. A red-shifted emission component to a number of near-IR emission lines was also uncovered. This is associated with a similar component in the H alpha line. V645 Cyg is also found to have variable CO first overtone bandhead emission. The data clearly indicate that the outflow of V645 Cyg is variable. The unidentified feature in a previously published optical spectrum is identified with a receding outflow at 2000 km per second. The nature of this feature, which is found in hydrogen and helium atomic lines and CO molecular lines remains a puzzle.

We present the initial results of a 3-mm spectral line survey towards 83 methanol maser selected massive star-forming regions. Here we report observations of the J=5-4 and 6-5 rotational transitions of methyl cyanide (CH3CN) and the J=1-0 transition of HCO+and H13CO+. CH3CN emission is detected in 58 sources (70 %) of our sample). We estimate the temperature and column density for 37 of these using the rotational diagram method. The temperatures we derive range from 28-166 K, and are lower than previously reported temperatures, derived from higher J transitions. We find that CH3CN is brighter and more commonly detected towards ultra-compact HII (UCHII) regions than towards isolated maser sources. Detection of CH3CN towards isolated maser sources strongly suggests that these objects are internally heated and that CH3CN is excited prior to the UCHII phase of massive star-formation. HCO+ is detected towards 82 sources (99 % of our sample), many of which exhibit asymmetric line profiles compared to H13CO+. Skewed profiles are indicative of inward or outward motions, however, we find approximately equal numbers of red and blue-skewed profiles among all classes. Column densities are derived from an analysis of the HCO+ and H13CO+ line profiles. 80 sources have mid-infrared counterparts: 68 seen in emission and 12 seen in absorption as `dark clouds'. Seven of the twelve dark clouds exhibit asymmetric HCO+ profiles, six of which are skewed to the blue, indicating infalling motions. CH3CN is also common in dark clouds, where it has a 90 % detection rate.

We present high resolution spectro-astrometry of a sample of 28 Herbig Ae/Be and 3 F-type pre-main sequence stars. The spectro-astrometry is shown from both empirical and simulated data to be capable of detecting binary companions that are fainter by up to 6 magnitudes at separations larger than 0.1 arcsec. The nine targets that were previously known to be a binary are all detected. In addition, we report the discovery of 6 new binaries and present 5 further possible binaries. The resulting binary fraction of 68+/-11 per cent is the largest reported for any observed sample of Herbig Ae/Be stars, presumably because of the exquisite sensitivity of spectro-astrometry for detecting binary systems. The data hint that the binary frequency of the Herbig Be stars is larger than for the Herbig Ae stars. The appendix presents model simulations to assess the capabilities of spectro-astrometry and reinforces the empirical findings. Two objects, HD 87643 and Z CMa, display evidence for asymmetric outflows. Finally, the position angles of the binary systems have been compared with available orientations of the circumprimary disc and these appear to be co-planar. The alignment between the circumprimary discs and the binary systems strongly suggests that the formation of binaries with intermediate mass primaries is due to fragmentation as the alternative, stellar capture, does not naturally predict aligned discs. The aligment extends to the most massive B-type stars in our sample. This leads us to conclude that formation mechanisms that do result in massive stars, but predict random angles beween the binaries and the circumprimary disks, such as stellar collisions, are also ruled out for the same reason.

We investigate the use of HI data to resolve the near/far ambiguity in kinematic distances of massive young stellar object (MYSO) candidates. Kinematic distances were obtained from 13CO 1-0 spectral line observations with the Mopra Telescope towards 94 candidates selected from the Red MSX Source (RMS) survey in the fourth Galactic quadrant. HI data from the Southern Galactic Plane Survey (SPGS) was used in conjunction with the HI self-absorption technique to determine the near or far distance. We resolved the kinematic distance ambiguity to 70% of the sources. We can also simultaneously solve for any multiple line-of-sight component sources. We discuss the advantages and disadvantages of this technique in comparison with other methods, and also perform confidence checks on the reliability of using the HI self-absorption technique. We examined the projected location of these sources in both the Galactic plane and longitude-velocity diagrams to ascertain any recognisable spiral arm pattern. Although no obvious spiral pattern was found when compared to that proposed by Cordes & Lazio, far distance sources tended to lie on or near spiral arm loci. Near distance sources, however, had peculiar velocity uncertainties consistent with the separation between the spiral arms themselves. The longitude-velocity plot shows a more consistent picture, with tangent points of the spiral arm loci easily seen. We conclude that using the HI self-absorption technique to determine kinematic distance ambiguities is a quick and reliable method in most cases, with an 80% success rate in determining the correct designation of whether an object is at the near or far distance.

We present the first systematic spectropolarimetric study of Luminous Blue Variables (LBVs), and find that at least half those objects studied display evidence for intrinsic polarization -- a signature of significant inhomogeneity at the base of the wind. Furthermore, multi-epoch observations reveal that the polarization is variable in both strength and position angle. This evidence points away from a simple axi-symmetric wind structure \`{a} la the B[e] supergiants, and instead suggests a wind consisting of localised density enhancements, or `clumps'. We show with an analytical model that, in order to produce the observed variability, the clumps must be large, produced at or below the photosphere, and ejected on timescales of days. More details of LBV wind-clumping will be determined through further analysis of the model and a polarimetric monitoring campaign.

In order to study the inner parts of the circumstellar material around optically faint, infrared bright objects, we present the first medium-resolution spectropolarimetric data taken in the near-infrared. In this paper we discuss Pa beta line data of GL 490, a well-known embedded massive young stellar object, and of MWC 349A and MWC 342, two optically faint stars that are proposed to be in the pre-main sequence phase of evolution. As a check on the method, the classical Be star zeta Tau, known to display line polarization changes at optical wavelengths, was observed as well. Three of our targets show a ``line effect'' across Pa beta. For zeta Tau and MWC 349A this line effect is due to depolarisation by a circumstellar electron-scattering disk. In both cases, the position angle of the polarisation is consistent with that of the larger scale disks imaged at other wavelengths, validating infrared spectropolarimetry as a means to detect flattening on small scales. The tentative detection of a rotation in the polarization position angle at Pa beta in the embedded massive young stellar object GL 490 suggests the presence of a small scale rotating accretion disk with an inner hole - similar to those recently discovered at optical wavelengths in Herbig Ae and T Tauri stars.

Stars of all evolutionary phases have been found to have excess infrared emission due to the presence of circumstellar material. To identify such stars, we have positionally correlated the infrared MSX point source catalogue and the Tycho 2 optical catalogue. A near/mid infrared colour criteria has been developed to select infrared excess stars. The search yielded 1938 excess stars, over half (979) have never previously been detected by IRAS. The excess stars were found to be young objects such as Herbig Ae/Be and Be stars, and evolved objects such as OH/IR and carbon stars. A number of B type excess stars were also discovered whose infrared colours could not be readily explained by known catalogued objects.

We present the first systematic spectropolarimetric study of Luminous Blue Variables (LBVs) in the Galaxy and the Magellanic Clouds, in order to investigate the geometries of their winds. We find that at least half of our sample show changes in polarization across the strong H$\alpha$ emission line, indicating that the light from the stars is intrinsically polarized and therefore that asphericity already exists at the base of the wind. Multi-epoch spectropolarimetry on four targets reveals variability in their intrinsic polarization. Three of these, AG Car, HR Car and P Cyg, show a position angle (PA) of polarization which appears random with time. Such behaviour can be explained by the presence of strong wind-inhomogeneities, or `clumps' within the wind. Only one star, R 127, shows variability at a constant PA, and hence evidence for axi-symmetry as well as clumpiness. However, if viewed at low inclination, and at limited temporal sampling, such a wind would produce a seemingly random polarization of the type observed in the other three stars. Time-resolved spectropolarimetric monitoring of LBVs is therefore required to determine if LBV winds are axi-symmetric in general. The high fraction of LBVs ($>$ 50%) showing intrinsic polarization is to be compared with the lower $\sim$ 20-25 % for similar studies of their evolutionary neighbours, O supergiants and Wolf-Rayet stars. We anticipate that this higher incidence is due to the lower effective gravities of the LBVs, coupled with their variable temperatures within the bi-stability jump regime. This is also consistent with the higher incidence of wind asphericity that we find in LBVs with strong H$\alpha$ emission and recent (last $\sim$ 10 years) strong variability.

We present Halpha spectropolarimetry observations of a sample of 10 bright T Tauri stars, supplemented with new Herbig Ae/Be star data. A change in the linear polarisation across Halpha is detected in most of the T Tauri (9/10) and Herbig Ae (9/11) objects, which we interpret in terms of a compact source of line photons that is scattered off a rotating accretion disk. We find consistency between the position angle (PA) of the polarisation and those of imaged disk PAs from infrared and millimetre imaging and interferometry studies, probing much larger scales. For the Herbig Ae stars AB Aur, MWC 480 and CQ Tau, we find the polarisation PA to be perpendicular to the imaged disk, which is expected for single scattering. On the other hand, the polarisation PA aligns with the outer disk PA for the T Tauri stars DR Tau and SU Aur and FU Ori, conforming to the case of multiple scattering. This difference can be explained if the inner disks of Herbig Ae stars are optically thin, whilst those around our T Tauri stars and FU Ori are optically thick. Furthermore, we develop a novel technique that combines known inclination angles and our recent Monte Carlo models to constrain the inner rim sizes of SU Aur, GW Ori, AB Aur, and CQ Tau. Finally, we consider the connection of the inner disk structure with the orientation of the magnetic field in the foreground interstellar medium: for FU Ori and DR Tau, we infer an alignment of the stellar axis and the larger magnetic field direction.

Spectroastrometry is a technique which has the potential to resolve flux distributions on scales of milliarcseconds. In this study, we examine the application of spectroastrometry to binary point sources which are spatially unresolved due to the observational point spread function convolution. The technique uses measurements with sub-pixel accuracy of the position centroid of high signal-to-noise long-slit spectrum observations. With the objects in the binary contributing fractionally more or less at different wavelengths (particularly across spectral lines), the variation of the position centroid with wavelength provides some information on the spatial distribution of the flux. We examine the width of the flux distribution in the spatial direction, and present its relation to the ratio of the fluxes of the two components of the binary. Measurement of three observables (total flux, position centroid and flux distribution width) at each wavelength allows a unique separation of the total flux into its component parts even though the angular separation of the binary is smaller than the observations' point-spread function. This is because we have three relevant observables for three unknowns (the two fluxes, and the angular separation of the binary), which therefore generates a closed problem. This is a wholly different technique than conventional deconvolution methods, which produce information on angular sizes of the sampling scale. Spectroastrometry can produce information on smaller scales than conventional deconvolution, and is successful in separating fluxes in a binary object with a separation of less than one pixel. We present an analysis of the errors involved in making binary object spectroastrometric measurements and the separation method, and highlight necessary observing methodology.

In this paper we present multi-epoch observations of the pre-main sequence binary HK Ori. These data have been drawn from the EXPORT database and are complemented by high quality spectro-astrometric data of the system. The spectroscopic data appear to be very well represented by a combination of an A dwarf star spectrum superposed on a (sub-)giant G-type spectrum. The radial velocity of the system is consistent with previous determinations, and does not reveal binary motion, as expected for a wide binary. The spectral, photometric and polarimetric properties and variability of the system indicate that the active object in the system is a T Tauri star with UX Ori characteristics. The spectro-astrometry of HK Ori is sensitive down to milli-arcsecond scales and confirms the speckle interferometric results from Leinert et al. The spectro-astrometry allows with fair certainty the identification of the active star within the binary, which we suggest to be a G-type T Tauri star based on its spectral characteristics.

We present emission line polarimetry data and modelling relevant to the circumstellar geometry and kinematics around pre-main sequence stars. For a sample of both Herbig Ae/Be stars and T Tauri stars, we find that most show polarization changes across Halpha, implying that flattened structures are common on the smallest scales -- and over a range of stellar masses. We also present Monte Carlo calculations of spectral line profiles scattered in rotating accretion disks. We consider both the case of a central star that emits line photons uniformly, as well as via hot spots. Intriguingly, the switch between a uniform point source and a finite-sized star results in a marked difference in the position angle variation across the line. Our models demonstrate the diagnostic potential of line polarimetry in determining the disk inclination and the size of the inner hole -- a spatial scale no other technique currently accesses.

We have selected red MSX sources (RMS) that have the colours of massive young stellar objects (MYSOs). Our aim is to generate a large, systematically selected sample to address questions such as their luminosity function, lifetimes, clustering and triggering. Other objects such as UCHIIs, PN, PPN and AGB stars have similar IR colours and a large programme of ground-based follow-up observations is underway to identify and eliminate these from the sample of the red MSX sources. These include radio continuum observations, kinematic distances, ground-based mid-IR imaging, near-IR imaging and spectroscopy to distinguish. We report the progress of these campaigns on the 3000 candidates, with initial indications showing that a substantial fraction are indeed massive YSOs.

The aim of this article is to demonstrate the useful role that can be played by spectropolarimetric observations of young and evolved emission line stars that analyse the linearly polarized component in their spectra. At the time of writing, this demonstration has to be made on the basis of optical data since there is no common-user infrared facility, in operation, that offers the desired combination of spectral resolution and sensitivity. Here we focus on what can be learned from linear spectropolarimetry alone at reasonably high spectral resolution and at $10^3 < $S/N$ < 10^4$. And we remind that the near infrared (1--2 micron) has the potential to out-perform the optical as a domain to work in because of the greatly reduced interstellar obscuration at these wavelengths. This point has been reached at a time when theory, exploiting flexible Monte Carlo methods, is fast becoming a powerful tool. In short we have the complex phenomena, and the rise of the modelling capability to match -- good data are the missing link.

We present high resolution (lambda / Delta_lambda = 49000) echelle spectra of the intermediate mass, pre-main sequence stars BF Ori, SV Cep, WW Wul and XY Per. The spectra cover the range 3800-5900 angstroms and monitor the stars on time scales of months and days. All spectra show a large number of Balmer and metallic lines with variable blueshifted and redshifted absorption features superimposed to the photospheric stellar spectra. Synthetic Kurucz models are used to estimate rotational velocities, effective temperatures and gravities of the stars. The best photospheric models are subtracted from each observed spectrum to determine the variable absorption features due to the circumstellar gas; those features are characterized in terms of their velocity, v, dispersion velocity, Delta v, and residual absorption, R_max. The absorption components detected in each spectrum can be grouped by their similar radial velocities and are interpreted as the signature of the dynamical evolution of gaseous clumps with, in most cases, solar-like chemical composition. This infalling and outflowing gas has similar properties to the circumstellar gas observed in UX Ori, emphasizing the need for detailed theoretical models, probably in the framework of the magnetospheric accretion scenario, to understand the complex environment in Herbig Ae (HAe) stars. WW Vul is unusual because, in addition to infalling and outflowing gas with properties similar to those observed in the other stars, it shows also transient absorption features in metallic lines with no obvious counterparts in the hydrogen lines. This could, in principle, suggest the presence of CS gas clouds with enhanced metallicity around WW Vul. The existence of such a metal-rich gas component, however, needs to be confirmed by further observations and a more quantitative analysis.

We present Halpha spectropolarimetry of a large sample of pre-main sequence (PMS) stars of low and intermediate mass, and argue that the technique is a powerful tool in studying the circumstellar geometry around these objects. For the intermediate mass (2 -- 15 Msun) Herbig Ae/Be stars we find that 16 out of 23 show a line effect, which immediately implies that flattening is common among these objects. Furthermore, we find a significant difference in Halpha spectropolarimetry behaviour between the Herbig Be and Ae groups. For the Herbig Be stars, the concept of an electron scattering disc is shown to be a useful concept to explain the depolarizations seen in this spectral range. At lower masses, more complex Halpha polarimetry behaviour starts to appear. The concept of a compact source of Halpha emission that is formed close to the stellar surface, for instance by hot spots due to magnetospheric accretion, is postulated as a working hypothesis to qualitatively explain the Halpha spectropolarimetry behaviour around Herbig Ae and lower mass (M < 2 Msun) T Tauri stars. The striking resemblance in spectropolarimetric behaviour between the T Tauri star RY Tau and the Herbig Ae stars suggests a common origin of the polarized line photons, and hints that low and higher mass pre-main sequence stars may have more in common than had hitherto been suspected.

We present linear Halpha spectropolarimetry of the classical T Tauri star RY Tau. A change in the polarization percentage and position angle across Halpha is detected, which suggests that line photons are scattered in a rotating disc. We derive the position angle from the slope of the loop in the (Q,U) diagram and find it to be 146 +/- 3 degrees. This is perpendicular to the position angle of the disc of 48 +/- 5 degrees as deduced from submillimeter imaging by Koerner & Sargent (1995). This finding is consistent, as scattering off the imaged millimeter disc is expected to yield a polarization signature in a direction that is rotated by 90 degrees from this disc. The observed spectropolarimetric behaviour of RY Tau is reminiscent of that seen in a large group of Herbig Ae stars, suggesting a common circumstellar origin of the polarized photons.

We present the results of a high spectral resolution study of the circumstellar (CS) gas around the intermediate mass, pre-main sequence star UX Ori. The results are based on a set of 10 echelle spectra, monitoring the star on time scales of months, days and hours. A large number of transient blueshifted and redshifted absorption features are detected in the Balmer and in many metallic lines. A multigaussian fit is applied to determine for each transient absorption the velocity, v, dispersion velocity, Delta v, and the parameter R, which provides a measure of the absorption strength of the CS gas. The time evolution of those parameters is presented and discussed. A comparison of intensity ratios among the transient absorptions suggests a solar-like composition of the CS gas. This confirms previous results and excludes a very metal-rich environment as the cause of the transient features in UX Ori. The features can be grouped by their similar velocities into 24 groups, of which 17 are redshifted and 7 blueshifted. An analysis of the velocity of the groups allows us to identify them as signatures of the dynamical evolution of 7 clumps of gas, of which 4 represent accretion events and 3 outflow events. Most of the events decelerate at a rate of tenths of ms-2, while 2 events accelerate at approximately the same rate; one event is seen experiencing both an acceleration and a deceleration phase and lasts for a period of few days. This time scale seems to be the typical duration of outflowing and infalling events in UX Ori. The dispersion velocity and the relative aborption strength of the features do not show drastic changes during the lifetime of the events, which suggests they are gaseous blobs preserving their geometrical and physical identity.

We present Halpha spectropolarimetry observations of a sample of 23 Herbig Ae/Be stars. A change in the linear polarisation across Halpha is detected in a large fraction of the objects, which indicates that the regions around Herbig stars are flattened (disc-like) on small scales. A second outcome of our study is that the spectropolarimetric signatures for the Ae stars differ from those of the Herbig Be stars, with characteristics changing from depolarisation across Halpha in the Herbig Be stars, to line polarisations in the Ae group. The frequency of depolarisations detected in the Herbig Be stars (7/12) is particularly interesting as, by analogy to classical Be stars, it may be the best evidence to date that the higher mass Herbig stars are surrounded by flattened structures. For the Herbig Ae stars, 9 out of 11 show a line polarisation effect that can be understood in terms of a compact Halpha emission that is itself polarised by a rotating disc-like circumstellar medium. The spectropolarimetric difference between the Herbig Be and Ae stars may be the first indication that there is a transition in the Hertzsprung-Russell Diagram from magnetic accretion at spectral type A to disc accretion at spectral type B. Alternatively, the interior polarised line emission apparent in the Ae stars may be masked in the Herbig Be stars due to their higher levels of Halpha emission.

The combination of mid-infrared data from the MSX satellite mission and ground-based near-infrared photometry is used to characterise the properties of the mid-infrared population of the Galactic plane. The colours of the youngest sources still heavily embedded within their natal molecular clouds are in general different from evolved stars shrouded within their own dust shells. Our main motivation is to use MSX for an unbiased search for a large (~1000) sample of massive young stellar objects (MYSOs). A simple analysis shows that the MSX point source catalogue should contain most of the MYSOs within our Galaxy. We develop colour selection criteria using combined near- and mid-infrared data for MYSOs, which produces a list of 3071 objects, excluding the galactic centre region. The programme of follow-up observations already underway to separate the MYSOs from compact H II regions and other remaining objects is briefly described. We also show that these data can be used, just as IRAS data has been previously, to provide a separation between evolved stars with carbon rich and oxygen rich dust. These data may also be used to search for evidence of dust around normal main sequence stars, such as low mass pre-main sequence stars or the Vega-excess class of objects where debris disks are presumed to remain from the planet formation process. We discuss the accuracy and completeness of the MSX point source catalogue, and show that the errors present tend to be of a kind that is not significant for the main stellar populations we discuss in this paper.

H_alpha spectropolarimetry on Herbig Ae/Be stars shows that the innermost regions of intermediate mass (2 -- 15 M_sun) Pre-Main Sequence stars are flattened. This may be the best evidence to date that the higher mass Herbig Be stars are embedded in circumstellar discs. A second outcome of our study is that the spectropolarimetric signatures for the lower mass Herbig Ae stars differ from those of the higher mass Herbig Be stars. Depolarisations across H_alpha are observed in the Herbig Be group, whereas line polarisations are common amongst the Herbig Ae stars in our sample. These line polarisation effects can be understood in terms of a compact H_alpha source that is polarised by a rotating disc-like configuration. The difference we detect between the Herbig Be and Ae stars may be the first indication that there is a transition in the Hertzsprung-Russell Diagram from magnetic accretion at spectral type A to disc accretion at spectral type B. However, it is also possible that the compact polarised line component, present in the Herbig Ae stars, is masked in the Herbig Be stars due to their higher levels of H_alpha emission.

This paper presents optical UBVRI broadband photo-polarimetry of the EXPORT sample obtained at the 2.5m Nordic Optical Telescope. The database consists of multi-epoch photo-polarimetry of 68 pre-main-sequence and main-sequence stars. An investigation of the polarization variability indicates that 22 objects are variable at the 3sigma level in our data. All these objects are pre-main sequence stars, consisting of both T Tauri and Herbig Ae/Be objects while the main sequence, Vega type and post-T Tauri type objects are not variable. The polarization properties of the variable sources are mostly indicative of the UXOR-type behaviour; the objects show highest polarization when the brightness is at minimum. We add seven new objects to the class of UXOR variables (BH Cep, VX Cas, DK Tau, HK Ori, LkHa 234, KK Oph and RY Ori). The main reason for their discovery is the fact that our data-set is the largest in its kind, indicating that many more young UXOR-type pre-main sequence stars remain to be discovered. The set of Vega-like systems has been investigated for the presence of intrinsic polarization. As they lack variability, this was done using indirect methods, and apart from the known case of BD +31.643, the following stars were found to be strong candidates to exhibit polarization due to the presence of circumstellar disks: 51 Oph, BD +31.643C, HD 58647 and HD 233517.

We now understand how low mass stars evolve through cloud collapse and disc accretion, but whether higher mass stars are also disc accretors is as yet unknown. Spectropolarimetry observations can help in answering this basic question, as they probe the first few stellar radii around young stars.

We analyse a sample of 236 Cepheids from the Hipparcos catalog, using the method of ``reduced parallaxes'' in V, I, K and the reddening-free ``Wesenheit-index''. We compare our sample to those considered by Feast & Catchpole (1997) and Lanoix et al. (1999). We stress the importance of two systematic effects which influence the distance to the LMC: the slopes of the Galactic PL-relations and metallicity corrections. In order to assess the influence of these various effects, we present 27 distance moduli (DM) to the LMC, based on three different colours (V,I,K), three different slopes, and three different metallicity corrections. The DM to the LMC based on the parallax data can be summarised as follows. Based on the PL-relation in V and I, and the Wesenheit-index, the DM is 18.60 \pm 0.11 (\pm 0.08 slope)(^{+0.08}_{-0.15} metallicity), which is our current best estimate. Based on the PL-relation in K the DM is 18.52 \pm 0.18$ (\pm 0.03 slope) (\pm 0.06 metallicity)(^{+0.10}_{-0} sampling bias}). The random error is mostly due to the given accuracy of the Hipparcos parallaxes and the number of Cepheids in the respective samples. The terms between parentheses indicate the possible systematic uncertainties due to the slope of the Galactic PL-relations, the metallicity corrections, and in the K-band, due to the limited number of stars. Recent work by Sandage et al. (1999) indicates that the effect of metallicity towards shorter distances may be smaller in V and I than indicated here.

We present time-resolved spectroscopy of the Be star HD 76534, which was observed to have an Halpha outburst in 1995, when the line went from photospheric absorption to emission at a level of more than two times the continuum within 2.5 hours. To investigate the short-term behaviour of the spectrum of HD 76534 we have obtained 30 spectra within two hours real-time and searched for variations in the spectrum. Within the levels of statistical significance, no variability was found. Rather than periodic on short time scales, the Halpha behaviour seems to be commonly episodic on longer (>1 year) time scales, as an assessment of the existing data on the Halpha line and the Hipparcos photometry suggests. HD 76534 underwent only 1 photometric outburst in the 3 year span that the star was monitored by the Hipparcos satellite.

A parameter study of the spectral evolution of a typical post-AGB star, with particular emphasis on the evolution of the infrared colours, is presented. The models are based on the latest evolutionary tracks for hydrogen burning post-AGB stars. For such tracks the evolutionary rate is very dependent on the assumed mass loss rate as a function of time. We investigate this effect by modifying the mass loss prescription. The newly calculated evolutionary rates and density distributions are used to model the spectral evolution of a post-AGB star with the photo-ionization code Cloudy, including dust in the radiative transfer. Different assumptions for the dust properties and dust formation are considered. It is shown that by varying these parameters in a reasonable way, entirely different paths are followed in the IRAS colour-colour diagram. First of all, the effects of the evolution of the central star on the expanding dust shell can not be neglected. Also the dust properties and the definition of the end of the AGB phase have an important effect. The model tracks show that objects occupying the same location in the IRAS colour-colour diagram can have a different evolutionary past, and therefore the position in the IRAS colour-colour diagram alone can not a priori give a unique determination of the evolutionary status of an object. An alternative colour-colour diagram, the K-[12] vs. [12]-[25] diagram, is presented. The tracks in this diagram seem less affected by particulars of the grain emission. This diagram may be a valuable additional tool for studying post-AGB evolution.

We present the results of medium resolution (Delta v = 60 km/s) spectropolarimetric observations across Halpha of a sample of B[e] and Herbig Be objects. A change in linear polarization across Halpha is detected in a large fraction of the objects. Halpha in the spectra of HD 37806 and HD 50138 each consist of a double-peaked polarized line and a superposed unpolarized single emission peak, suggesting two distinct line-forming regions. Multiple observations of HD 45677 allow for the separation of electron and dust scattering effects for the first time: the data indicate that the dust-scattering region is clumpy. Two unexpected results are the non-detections of Halpha polarization changes in omega Ori, where depolarization has previously been detected, and in MWC 297, which exhibits source elongation at radio wavelengths. In omega Ori time variability is probably responsible such that this star's electron-scattering disk was much weakened at the time of observation. Two hypotheses are advanced that might explain the MWC 297 result. The general findings are that roughly half of the observed Herbig Be stars show polarization changes across Halpha, implying immediately that their ionized envelopes are not spherically symmetric. This pattern, if confirmed by observations of a larger sample, could indicate that the non-detection rate is simply a consequence of sampling randomly-oriented circumstellar disks able to scatter starlight within a few stellar radii. The stars classified as B[e] stars all show startling polarization changes across Halpha. The details in each case are different, but the widely accepted concept of dense Halpha emitting equatorial disks around these objects is supported.

We investigate the properties of K0V stars with Hipparcos parallaxes and spectral types taken from the Michigan Spectral Survey. The sample of 200 objects allows the empirical investigation of the magnitude selection (Malmquist) bias, which appears clearly present. By selecting those objects that are not affected by bias, we find a mean absolute magnitude of Mv~5.7, a downward revision from 5.9 mag. listed in Schmidt-Kaler (1982). Some objects have absolute magnitudes far brighter than Mv~5.7, and it is suggested that these objects (~20% of the total sample) are K0IV stars which may have been mis-classified as a K0V star. The presence of the Malmquist bias in even this high quality sample suggests that no sample can be expected to be bias-free.

Using the Planetary Camera on board the Hubble Space Telescope we have measured the projected separation of the binary components in the nucleus of the planetary nebula Abell 35 to be larger than 0.08 arcsecs but less than 0.14 arcsecs. The system was imaged in three filters centered at 2950, 3350 and 5785 Ang. The white dwarf primary star responsible for ionizing the nebula is half as bright as its companion in the 2950 Ang filter causing the source to be visibly elongated. The 3350 Ang setting, on the other hand, shows no elongation as a result of the more extreme flux ratio. The F300W data allows the determinination of the binary's projected separation. At the minimum distance of 160 parsec to the system, our result corresponds to 18+/-5 AU. This outcome is consistent with the wind accretion induced rapid rotation hypothesis, but cannot be reconciled with the binary having emerged from a common-envelope phase.

New high resolution spectroscopic and medium resolution spectropolarimetric data, complemented with optical broad and narrow band imaging, of the B[e] star HD 87643 are presented. The spectrum of HD 87643 exhibits the hybrid characteristics well known to be representative of the group of B[e] stars; a fast wind with an expansion velocity in excess of 1000 km/s is measured in the hydrogen and helium lines, while a slower component is traced by lower excitation lines and forbidden lines. Clues to the geometry of the rapidly expanding circumstellar shell are provided by the startling polarization changes across Halpha. Comparison with published schematic calculations indicates that the polarizing material is located in a slowly rotating, expanding disk structure. A hydrodynamical model is then presented whose results are consistent with the original two-wind concept for B[e] stars and exhibits kinematic properties that may well explain the observed spectral features in HD 87643. The model calculations use as input a B star undergoing mass loss, surrounded by an optically thick disk. The resulting configuration consists of a fast polar wind from the star and a slowly expanding disk wind. The model also predicts that the stellar wind at intermediate latitudes is slower and denser than in the polar region.

The theoretical prediction that trigonometric parallaxes suffer from a statistical effect, has become topical again now that the results of the Hipparcos satellite have become available. This statistical effect, the so-called Lutz-Kelker bias, causes measured parallaxes to be too large. This has the implication that inferred distances, and hence inferred luminosities are too small. Published analytic calculations of the Lutz-Kelker bias indicate that the inferred luminosity of an object is, on average, 30% too small when the error in the parallax is only 17.5%. Yet, this bias has never been determined empirically. In this paper we investigate whether there is such a bias by comparing the best Hipparcos parallaxes which ground-based measurements. We find that there is indeed a large bias affecting parallaxes, with an average and scatter comparable to predictions. We propose a simple method to correct for the LK bias, and apply it successfully to a sub-sample of our stars. We then analyze the sample of 26 `best' Cepheids used by Feast & Catchpole (1997) to derive the zero-point of the fundamental mode pulsators and leads to a distance modulus to the Large Magellanic Cloud - based on Cepheid parallaxes- of 18.56 +/- 0.08, consistent with previous estimates.

We have compared the B-V colour excess, E(B-V), obtained for a sample of five optically visible massive YSOs both from diffuse interstellar bands (DIBs) in their spectra and from their optical continuum slopes. Our targets are HD 200775, BD+40 4124, MWC 1080, MWC 297 and MWC 349A. First, E(B-V) towards each of the targets is derived by dereddening the observed continua to match those of B-type standard stars. A survey of DIBs in the spectra of the massive YSOs, and a control field star, then reveals that the DIBs are significantly weaker in the former than would be expected based on the total E(B-V) values. This result is strengthened by the finding that the DIBs in the control field star, HD 154445, have on average the strength expected from its continuum E(B-V). A rough estimate of the foreground reddening of intervening diffuse interstellar medium shows it to be smaller than the DIB-E(B-V), implying that at least part of the DIB carriers are formed within the parental molecular clouds in which the YSOs are embedded. The formation efficiency of the DIBs varies strongly however from cloud to cloud. The DIB-E(B-V) compares favourably with the total E(B-V) towards BD+40 4124, but is almost negligible in the line of sight towards MWC 297. Despite this general, but not unexpected, deficit we provide evidence that the DIB at 5849 Angstrom is a good tracer of total extinction in these lines of sight.

We have investigated the report of Tamura and Takeuti that the Halpha line of the F-type post-AGB star HD56126 is variable on time scales of minutes. To this end, HD56126 was observed on two occasions with the William Herschel Telescope. Seventeen, respectively thirty spectra were taken within time span of 1.5 hours in order to detect any short term variations. We find that the Halpha line profile changed strongly over the two month interval, but no evidence is found for short term variability. The variability Tamura and Takeuti claim to find is probably due to the low signal-to-noise in their spectra.