Source author record

R. D. Oudmaijer

R. D. Oudmaijer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

49works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

49 published item(s)

preprint2023arXiv

The Co-Ordinated Radio and Infrared Survey for High-Mass Star Formation. V. The CORNISH-South Survey and Catalogue

We present the first high spatial resolution radio continuum survey of the southern Galactic plane. The CORNISH project has mapped the region defined by $295^{\circ} < l < 350^{\circ}$; $|b| < 1^{\circ}$ at 5.5-GHz, with a resolution of 2.5$^{''}$ (FWHM). As with the CORNISH-North survey, this is designed to primarily provide matching radio data to the Spitzer GLIMPSE survey region. The CORNISH-South survey achieved a root mean square noise level of $\sim$ 0.11 mJy beam$^{-1}$, using the 6A configuration of the Australia Telescope Compact Array (ATCA). In this paper, we discuss the observations, data processing and measurements of the source properties. Above a 7$σ$ detection limit, 4701 sources were detected, and their ensemble properties show similar distributions with their northern counterparts. The catalogue is highly reliable and is complete to 90 per cent at a flux density level of 1.1 mJy. We developed a new way of measuring the integrated flux densities and angular sizes of non-Gaussian sources. The catalogue primarily provides positions, flux density measurements and angular sizes. All sources with IR counterparts at 8$μm$ have been visually classified, utilizing additional imaging data from optical, near-IR, mid-IR, far-IR and sub-millimetre galactic plane surveys. This has resulted in the detection of 524 H II regions of which 255 are ultra-compact H II regions, 287 planetary nebulae, 79 radio stars and 6 massive young stellar objects. The rest of the sources are likely to be extra-galactic. These data are particularly important in the characterization and population studies of compact ionized sources such as UCHII regions and PNe towards the Galactic mid-plane.

preprint2022arXiv

Tracing a decade of activity towards a yellow hypergiant. The spectral and spatial morphology of IRC+10420 at au scales

The fate of a massive star during the latest stages of its evolution is highly dependent on its mass-loss history and geometry, with the yellow hypergiants being key objects to study those phases of evolution. We present near-IR interferometric observations of the famous yellow hypergiant IRC +10420 and blue spectra taken between 1994-2019. Our 2.2 $μ$m GRAVITY/VLTI observations attain a spatial resolution of $\sim$5 stellar radii and probe the hot emission in the K-band tracing the gas via Na i double emission and the Br$γ$ emission. The observed configurations spatially resolve the 2.2 $μ$m continuum as well as the Br$γ$ and the Na i emission lines. Our geometric modelling demonstrates the presence of a compact neutral zone (Na i) which is slightly larger than the continuum but within an extended Br$γ$ emitting region. Our geometric models of the Br$γ$ emission confirm an hour-glass geometry of the wind. To explain this peculiar geometry we investigate the presence of a companion at 7-800 au separations and find no signature at the contrast limit of our observations (3.7 mag at 3$σ$). We report an evolution of the ejecta over a time span of 7 years, which allows us to constrain the opening angle of the hour-glass geometry at $<$10$^\circ$. Lastly, we present the first blue optical spectra of IRC +10420 since 1994. The multi-epoch data indicate that the spectral type, and thus temperature, of the object has essentially remained constant during the intervening years. This confirms earlier conclusions that following an increase in temperature of 2000 K in less than two decades prior to 1994, the temperature increase has halted. This suggests that this yellow hypergiant has "hit" the White Wall in the HR-diagram preventing it from evolving blue-wards, and will likely undergo a major mass-loss event in the near future.

preprint2020arXiv

Catalogue of new Herbig Ae/Be and classical Be stars. A machine learning approach to Gaia DR2

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$α$ 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.

preprint2020arXiv

Discovery of a jet from the single HAe/Be star HD 100546

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$α$ 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.

preprint2020arXiv

Optical/near-infrared observations of the Fried Egg Nebula: Multiple shell ejections on a 100 yr timescale from a massive yellow hypergiant

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 μm, 11.85 μm and 12.81 μm simultaneously, adopting the revised Gaia distance (DR2). Results. The interferometric observables around 2 μm show that the Brγ 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γ 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.

preprint2016arXiv

Dissecting a SN impostor's circumstellar medium: MUSEing about the SHAPE of eta Car's outer ejecta

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.

preprint2016arXiv

Linear spectropolarimetry across the optical spectrum of Herbig Ae/Be stars

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

preprint2016arXiv

The compact H$α$ emitting regions of the Herbig Ae/Be stars HD 179218 and HD 141569 from CHARA spectro-interferometry

This work presents CHARA/VEGA H$α$ spectro-interferometry (R ~ 6000, and $λ$/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$α$ emission is spatially unresolved, indicating that the size of the H$α$ 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$α$ 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$α$ 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$α$ emission in HAeBe stars -compact or extended, discs, accretion, and winds-, in agreement with previous Br$γ$ spectro-interferometric observations.

preprint2016arXiv

Tracing jet emission at the base of a high-mass YSO. First AMBER/VLTI observations of the Brγemission in IRAS 13481-6124

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γ-emission line through NIR interferometry. We present the first AMBER/VLTI observations of the Brγ-emitting region in an HMYSO at R~1500. Our AMBER/VLTI observations reveal a spatially and spectrally resolved Brγ-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γ-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.

preprint2015arXiv

A spectroscopic survey of Herbig Ae/Be stars with X-Shooter I: Stellar parameters and accretion rates

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.

preprint2015arXiv

High-resolution Br-gamma spectro-interferometry of the transitional Herbig Ae/Be star HD 100546: a Keplerian gaseous disc inside the inner rim

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.

preprint2015arXiv

Investigating the nature of the Fried Egg nebula: CO mm-line and optical spectroscopy of IRAS 17163-3907

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

preprint2015arXiv

On the origin of the correlations between the accretion luminosity and emission line luminosities in pre-main sequence stars

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.

preprint2015arXiv

Spectroscopy and Linear Spectropolarimetry of the early Herbig Be Stars PDS 27 and PDS 37

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$α$ 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.

preprint2014arXiv

Stellar parameters and accretion rate of the transition disk star HD 142527 from X-Shooter

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.

preprint2013arXiv

Accretion variability of Herbig Ae/Be stars observed by X-Shooter. HD 31648 and HD 163296

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.

preprint2013arXiv

Investigation of dust properties of the proto-planetary nebula IRAS 18276-1431

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.

preprint2013arXiv

Near-infrared integral field spectroscopy of Massive Young Stellar Objects

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$γ$, CO, H$_2$, and \mbox{He\,{\sc i}}. The Br$γ$ 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$γ$ 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$γ$. 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 $ν=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 $ν=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.

preprint2013arXiv

The Red MSX Source Survey: the Massive Young Stellar Population of our Galaxy

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.

preprint2013arXiv

The RMS Survey: Galactic distribution of massive star formation

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.

preprint2013arXiv

The RMS Survey: Near-IR Spectroscopy of Massive Young Stellar Objects

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.

preprint2012arXiv

Accretion-related properties of Herbig Ae/Be stars. Comparison with T Tauris

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 τ= 1.3^{+1.0}_{-0.5} Myr from an exponential law fit, while a power law yields Macc(t) \propto t^{-η}, with η= 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.

preprint2012arXiv

AMBER and CRIRES observations of the binary sgB[e] star HD 327083: evidence of a gaseous disc traced by CO bandhead emission

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

preprint2012arXiv

Probing the envelopes of massive young stellar objects with diffraction limited mid-infrared imaging

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.

preprint2012arXiv

Probing the properties of Be star discs with spectroastrometry and NLTE radiative transfer modelling: beta CMi

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.

preprint2012arXiv

The Co-ordinated Radio and Infrared Survey for High Mass Star Formation (The CORNISH Survey) - I. Survey Design

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.

preprint2012arXiv

Tracers of Discs and Winds around Intermediate and High Mass Young Stellar Objects

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$μ$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.

preprint2012arXiv

Unveiling the structure and kinematics of B[e] stars' disks from FEROS and CRIRES spectra

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.

preprint2012arXiv

VLTI/AMBER observations of the binary B[e] supergiant HD 327083

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.

preprint2011arXiv

A double detached shell around a post-Red Supergiant:IRAS 17163-3907, the Fried Egg nebula

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.

preprint2011arXiv

Accretion rates and accretion tracers of Herbig Ae/Be stars

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.

preprint2011arXiv

On the alignment between the circumstellar disks and orbital planes of Herbig Ae/Be binary systems

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

preprint2011arXiv

Optical spectroscopic variability of Herbig Ae/Be stars

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.

preprint2011arXiv

The RMS Survey: Ammonia and water maser analysis of massive star forming regions

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.

preprint2011arXiv

The RMS Survey: Resolving kinematic distance ambiguities towards a sample of compact HII regions using HI absorption

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.

preprint2010arXiv

Mid-infrared interferometry towards the massive young stellar object CRL 2136: inside the dust rim

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.

preprint2010arXiv

The RMS Survey: Distribution and properties of a sample of massive young stars

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.

preprint2010arXiv

The RMS Survey: The Bolometric Fluxes and Luminosity Distributions of Young Massive Stars

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{}.

preprint2009arXiv

On the presence and absence of disks around O-type stars

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

preprint2009arXiv

The mass ratio and formation mechanisms of Herbig Ae/Be star binary systems

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.

preprint2009arXiv

The origin of mid-infrared emission in massive young stellar objects: multi-baseline VLTI observations of W33A

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.

preprint2009arXiv

The RMS Survey: 6 cm continuum VLA observations towards candidate massive YSOs in the northern hemisphere

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

preprint2009arXiv

The RMS Survey: Far-Infrared Photometry of Young Massive Stars

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.

preprint2009arXiv

The RMS Survey: H2O masers towards a sample of southern hemisphere massive YSO candidates and ultra compact HII regions

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$σ$ 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.

preprint2008arXiv

Young massive stars and their environment in the mid-infrared at high angular resolution

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.

preprint2006arXiv

Evidence for variable outflows in the Young Stellar Object V645 Cygni

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.

preprint2006arXiv

On the H$α$ emission from the $β$ Cephei system

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 $β$ Cephei is a very slow rotator that still shows H$α$ emission episodes like in other Be stars, contradicting current theories. We investigate the hypothesis that the H$α$ emission stems from the spectroscopically unresolved companion of $β$ Cep. Spectra of the two unresolved components have been separated in the 6350-6850Årange with spectro-astrometric techniques, using 11 longslit spectra obtained with ALFOSC at the Nordic Optical Telescope, La Palma. We find that the H$α$ emission is not related to the primary in $β$ 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$α$ emission behaviour, the enigma behind the Be status of the slow rotator $β$ Cep has been resolved.

preprint1999arXiv

The spectral evolution of post-AGB stars

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.

preprint1998arXiv

The absolute magnitude of K0V stars from HIPPARCOS parallaxes

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.