Showing 1–50 of 86 results for author: Lumsden, S L

The RMS survey: a census of massive YSO multiplicity in the $K$-band

Authors: Robert G. Shenton, Rebecca J. Houghton, René D. Oudmaijer, Simon P. Goodwin, Stuart L. Lumsden, Evgenia Koumpia, Maria Koutoulaki

Abstract: Close to 100 per cent of massive stars are thought to be in binary systems. The multiplicity of massive stars seems to be intrinsically linked to their formation and evolution, and Massive Young Stellar Objects are key in observing this early stage of star formation. We have surveyed three samples totalling hundreds of MYSOs ($>8M_\odot$) across the Galaxy from the RMS catalogue, using UKIDSS and… ▽ More Close to 100 per cent of massive stars are thought to be in binary systems. The multiplicity of massive stars seems to be intrinsically linked to their formation and evolution, and Massive Young Stellar Objects are key in observing this early stage of star formation. We have surveyed three samples totalling hundreds of MYSOs ($>8M_\odot$) across the Galaxy from the RMS catalogue, using UKIDSS and VVV point source data, and UKIRT $K-$band imaging to probe separations between 0.8-9 arcsec (approx 1000-100,000 au). We have used statistical methods to determine the binary statistics of the samples, and we find binary fractions of $64\pm 4$ per cent for the UKIDSS sample, $53\pm 4$ per cent for the VVV sample, and $49\pm 8$ per cent for the RMS imaging sample. Also we use the $J-$ and $K-$band magnitudes as a proxy for the companion mass, and a significant fraction of the detected systems have estimated mass ratios greater than 0.5, suggesting a deviation from the capture formation scenario which would be aligned with random IMF sampling. Finally, we find that YSOs located in the outer Galaxy have a higher binary fraction than those in the inner Galaxy. This is likely due to a lower stellar background density than observed towards the inner Galaxy, resulting in higher probabilities for visual binaries to be physical companions. It does indicate a binary fraction in the probed separation range of close to 100 per cent without the need to consider selection biases. △ Less

Submitted 28 November, 2023; originally announced November 2023.

Comments: 14 pages, 9 figures, accepted to MNRAS

Kinematics and stability of high-mass protostellar disk candidates at sub-arcsecond resolution -- Insights from the IRAM NOEMA large program CORE

Authors: Aida Ahmadi, H. Beuther, F. Bosco, C. Gieser, S. Suri, J. C. Mottram, R. Kuiper, Th. Henning, Á. Sánchez-Monge, H. Linz, R. E. Pudritz, D. Semenov, J. M. Winters, T. Möller, M. T. Beltrán, T. Csengeri, R. Galván-Madrid, K. G. Johnston, E. Keto, P. D. Klaassen, S. Leurini, S. N. Longmore, S. L. Lumsden, L. T. Maud, L. Moscadelli , et al. (6 additional authors not shown)

Abstract: The fragmentation mode of high-mass molecular clumps and the accretion processes that form the most massive stars ($M\gtrsim 8M_\odot$) are still not well understood. To this end, we have undertaken a large observational program (CORE) making use of interferometric observations from the Northern Extended Millimetre Array (NOEMA) for a sample of 20 luminous ($L>10^4L_\odot$) protostellar objects in… ▽ More The fragmentation mode of high-mass molecular clumps and the accretion processes that form the most massive stars ($M\gtrsim 8M_\odot$) are still not well understood. To this end, we have undertaken a large observational program (CORE) making use of interferometric observations from the Northern Extended Millimetre Array (NOEMA) for a sample of 20 luminous ($L>10^4L_\odot$) protostellar objects in the 1.37 mm wavelength regime in both continuum and line emission, reaching $\sim$0.4" resolution (800 au at 2 kpc). Using the dense gas tracer CH$_3$CN, we find velocity gradients across 13 cores perpendicular to the directions of bipolar molecular outflows, making them excellent disk candidates. Specific angular momentum ($j$) radial profiles are on average $\sim10^{-3}$ km /s pc and follow $j \propto r^{1.7}$, consistent with a poorly resolved rotating and infalling envelope/disk model. Fitting the velocity profiles with a Keplerian model, we find protostellar masses in the range of $\sim 10-25$ $M_\odot$. Modelling the level population of CH$_3$CN lines, we present temperature maps and find median gas temperatures in the range $70-210$ K. We create Toomre $Q$ maps to study the stability of the disks and find almost all (11 of 13) disk candidates to be prone to fragmentation due to gravitational instabilities at the scales probed by our observations. In particular, disks with masses greater than $\sim10-20\%$ of the mass of their host (proto)stars are Toomre unstable, and more luminous protostellar objects tend to have disks that are more massive and hence more prone to fragmentation. Our finings show that most disks around high-mass protostars are prone to disk fragmentation early in their formation due to their high disk to stellar mass ratio. This impacts the accretion evolution of high-mass protostars which will have significant implications for the formation of the most massive stars. △ Less

Submitted 3 May, 2023; v1 submitted 28 April, 2023; originally announced May 2023.

Comments: 27 pages, 12 figures, 6 appendices - accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 677, A171 (2023)

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

Authors: T. Irabor, M. G. Hoare, M. Burton, W. D. Cotton, P. Diamond, S. Dougherty, S. P. Ellingsen, R. Fender, G. A. Fuller, S. Garrington, P. F. Goldsmith, J. Green, A. G. Gunn, J. Jackson, S. Kurtz, S. L. Lumsden, J. Marti, I. McDonald, S. Molinari, T. J. Moore, M. Mutale, T. Muxlow, T. OBrien, R. D. Oudmaijer, R. Paladini , et al. (10 additional authors not shown)

Abstract: 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-So… ▽ More 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. △ Less

Submitted 5 January, 2023; originally announced January 2023.

First spatially resolved Na I and He I transitions towards an MYSO. Finding new tracers for the gaseous star/disc interface

Authors: Evgenia Koumpia, M. Koutoulaki, W. -J. de Wit, R. D. Oudmaijer, A. J. Frost, S. L. Lumsden, J. M. Pittard

Abstract: With steady observational advances, the formation of massive stars is being understood in more detail. Numerical models are converging on a scenario where accretion discs play a key role. Direct observational evidence of such discs at a few au scales is scarce, due to the rarity of such objects and the observational challenges, including the lack of adequate diagnostic lines in the near-IR. We pre… ▽ More With steady observational advances, the formation of massive stars is being understood in more detail. Numerical models are converging on a scenario where accretion discs play a key role. Direct observational evidence of such discs at a few au scales is scarce, due to the rarity of such objects and the observational challenges, including the lack of adequate diagnostic lines in the near-IR. We present the analysis of K-band spectro-interferometric observations toward the Massive Young Stellar Object IRAS 13481-6124, which is known to host an accreting dusty disc. Using GRAVITY on the VLTI, we trace the crucial au-scales of the warm inner interface between the star and the accretion dusty disc. We detect and spatially resolve the Na I doublet and He I transitions towards an object of this class for the first time. The new observations in combination with our geometric models allowed us to probe the smallest au-scales of accretion/ejection around an MYSO. We find that Na I originates in the disc at smaller radii than the dust disc and is more compact than any of the other spatially resolved diagnostics (Br$γ$, He I, and CO). Our findings suggest that Na I can be a new powerful diagnostic line in tracing the warm star/disc accreting interface of forming (massive) stars, while the similarities between He I and Br$γ$ point towards an accretion/ejection origin of He I △ Less

Submitted 23 November, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS Letters; 5 pages, 3 Figures

The Spatial Evolution of Young Massive Clusters III. Effect of the Gaia Filter on 2D Spatial Distribution Studies

Authors: Anne S. M. Buckner, Zeinab Khorrami, Marta González, Stuart L. Lumsden, Paul Clark, Estelle Moraux

Abstract: [Context.] Gaia is limited in the optical down to G~21 mag so it is essential to understand the biases introduced by a magnitude-limited sample on spatial distribution studies. [Aims.] We ascertain how sample incompleteness in Gaia observations of young clusters affects the local spatial analysis tool INDICATE and subsequently the perceived spatial properties of these clusters. [Methods.] We creat… ▽ More [Context.] Gaia is limited in the optical down to G~21 mag so it is essential to understand the biases introduced by a magnitude-limited sample on spatial distribution studies. [Aims.] We ascertain how sample incompleteness in Gaia observations of young clusters affects the local spatial analysis tool INDICATE and subsequently the perceived spatial properties of these clusters. [Methods.] We created a mock Gaia cluster catalogue from a synthetic dataset using the observation generating tool MYOSOTIS. The effect of cluster distance, uniform and variable extinction, binary fraction, population masking by the point spread function wings of high-mass members, and contrast sensitivity limits on the trends identified by INDICATE are explored. A comparison of the typical index values derived by INDICATE for members of the synthetic dataset and their corresponding mock Gaia catalogue observations is made to identify any significant changes. [Results.] We typically find only small variations in the pre- and post-observation index values of cluster populations, which can increase as a function of incompleteness percentage and binarity. No significant strengthening, or false signatures, of stellar concentrations are found but real signatures may be diluted. Conclusions drawn about the spatial behaviour of Gaia-observed cluster populations that are, and are not, associated with their natal nebulosity are reliable for most clusters but the perceived behaviours of individual members can change so INDICATE should be used as a measure of spatial behaviours between members as a function of their intrinsic properties (e.g. mass, age, object type), rather than to draw conclusions about any specific observed member. [Conclusions.] INDICATE is a robust spatial analysis tool to reliably study Gaia-observed young cluster populations within 1 kpc, up to a sample incompleteness of 83.3% and binarity of 50%. △ Less

Submitted 15 December, 2021; v1 submitted 23 November, 2021; originally announced November 2021.

Comments: Accepted for publication in A&A. 20 pages, 4 figures, 2 manuscript tables, appendix with 9 reader reference tables

Journal ref: A&A 659, A72 (2022)

Tying the geometrical traits of massive young stellar objects and their discs to a potential evolutionary sequence using infrared observations

Authors: A. J. Frost, R. D. Oudmaijer, S. L. Lumsden, W-J de Wit

Abstract: Young massive stars influence their surroundings from local to galactic scales, but the observational challenges associated with their distance and embedded nature has, until the recent decade, made high-resolution studies of these objects difficult. In particular, comparative analyses of massive young stellar object (MYSO) discs are currently lacking and our understanding of their evolution is li… ▽ More Young massive stars influence their surroundings from local to galactic scales, but the observational challenges associated with their distance and embedded nature has, until the recent decade, made high-resolution studies of these objects difficult. In particular, comparative analyses of massive young stellar object (MYSO) discs are currently lacking and our understanding of their evolution is limited. Here, we combine the results of two studies with the aim to attribute geometrical features to an evolutionary sequence for a sample of seven MYSOs. The time evolution is based on a near-IR spectral features, while the geometry is determined from a multi size-scale study of MYSOs. We find that MYSO discs with determined geometrical substructure turn out to be also spectroscopically more evolved. This implies that disc evolution and dispersal are occurring within MYSOs, similar to low-mass YSO disc evolution, despite their faster formation timescales. △ Less

Submitted 12 August, 2021; v1 submitted 11 August, 2021; originally announced August 2021.

A Galactic survey of radio jets from massive protostars

Authors: S. J. D. Purser, S. L. Lumsden, M. G. Hoare, S. Kurtz

Abstract: In conjunction with a previous southern-hemisphere work, we present the largest radio survey of jets from massive protostars to date with high-resolution, ($\sim 0.04^{\prime\prime}$) VLA observations towards two subsamples of massive star-forming regions of different evolutionary statuses: 48 infrared-bright, massive YSOs and 8 IRDCs containing 16 luminous (${\rm L_{bol}}>10^3\,{\rm L_\odot}$) co… ▽ More In conjunction with a previous southern-hemisphere work, we present the largest radio survey of jets from massive protostars to date with high-resolution, ($\sim 0.04^{\prime\prime}$) VLA observations towards two subsamples of massive star-forming regions of different evolutionary statuses: 48 infrared-bright, massive YSOs and 8 IRDCs containing 16 luminous (${\rm L_{bol}}>10^3\,{\rm L_\odot}$) cores. For $94\%$ of the MYSO sample we detect thermal radio ($α\geq -0.1$ whereby $S_ν\propto ν^α$) sources coincident with the protostar, of which $84\%$ (13 jets and 25 candidates) are jet-like. Radio luminosity is found to scale with ${\rm L_{bol}}$ similarly to the low-mass case supporting a common mechanism for jet production across all masses. Associated radio lobes tracing shocks are seen towards $52\%$ of jet-like objects and are preferentially detected towards jets of higher radio and bolometric luminosities, resulting from our sensitivity limitations. We find jet mass loss rate scales with bolometric luminosity as $\dot{m}_{\rm jet} \propto {\rm L_{bol}}^{0.9\pm0.2}$, thereby discarding radiative, line-driving mechanisms as the dominant jet-launching process. Calculated momenta show that the majority of jets are mechanically capable of driving the massive, molecular outflow phenomena since $p_{\rm jet}>p_{\rm outflow}$. Finally, from their physical extent we show that the radio emission can not originate from small, optically-thick \textsc{Hii} regions. Towards the IRDC cores, we observe increasing incidence rates/radio fluxes with age using the proxy of increasing luminosity-to-mass $\left( \frac{L}{M} \right)$ and decreasing infrared flux ratios $\left(\frac{S_{70 \rm μm}}{S_{24 \rm μm}}\right)$. Cores with $\frac{L}{M}<40\,{\rm L_{sol}}{\rm M_{sol}}^{-1}$ are not detected above ($5.8{\rm GHz}$) radio luminosities of $\sim1{\rm mJy\,kpc}^2$. △ Less

Submitted 16 March, 2021; originally announced March 2021.

Comments: To be published in MNRAS. 20 pages main text (4 tables, 14 figures), 90 pages appendices (6 tables, 66 figures)

Unveiling the traits of massive young stellar objects through a multi-scale survey

Authors: A. J. Frost, R. D. Oudmaijer, W. J. de Wit, S. L. Lumsden

Abstract: The rarity and deeply embedded nature of stars with masses larger than 8 solar masses has limited our understanding of their formation. Previous work has shown that complementing spectral energy distributions with interferometric and imaging data can probe the circumstellar environments of massive young stellar objects (MYSOs) well. However, complex studies of single objects often use different ap… ▽ More The rarity and deeply embedded nature of stars with masses larger than 8 solar masses has limited our understanding of their formation. Previous work has shown that complementing spectral energy distributions with interferometric and imaging data can probe the circumstellar environments of massive young stellar objects (MYSOs) well. However, complex studies of single objects often use different approaches in their analysis. Therefore the results of these studies cannot be directly compared. This work aims to obtain the physical characteristics of a sample of MYSOs at ~0.01" scales, at 0.1" scales, and as a whole, which enables us to compare the characteristics of the sources. We apply the same multi-scale method and analysis to a sample of MYSOs. High-resolution interferometric data, near-diffraction-limited imaging data, and a multi-wavelength spectral energy distribution are combined. By fitting simulated observables derived from 2.5D radiative transfer models of disk-outflow-envelope systems to our observations, the properties of the MYSOs are constrained. We find that the observables of all the MYSOs can be reproduced by models with disk-outflow-envelope geometries, analogous to the Class I geometry associated with low-mass protostars. The characteristics of the envelopes and the cavities within them are very similar across our sample. On the other hand, the disks seem to differ between the objects, in particular with regards to what we interpret as evidence of complex structures and inner holes. This is comparable to the morphologies observed for low-mass young stellar objects. A strong correlation is found between the luminosity of the central MYSO and the size of the transition disk-like inner hole for the MYSOs, implying that photoevaporation or the presence of binary companions may be the cause. △ Less

Submitted 9 February, 2021; originally announced February 2021.

Journal ref: A&A 648, A62 (2021)

A Multi-Epoch Study of Radio Continuum Emission from MassiveProtostars

Authors: W. O. Obonyo, S. L. Lumsden, M. G. Hoare, S. E. Kurtz, S. J. D Purser

Abstract: We report the results of the Jansky Very Large Array (JVLA) observation of five massive protostars at 6 and 22.2 GHz. The aim of the study was to compare their current fluxes and positions with previous observations to search for evidence of variability. Most of the observed sources present the morphologies of a thermal core, hosting the protostar and exhibiting no proper motion, and associated no… ▽ More We report the results of the Jansky Very Large Array (JVLA) observation of five massive protostars at 6 and 22.2 GHz. The aim of the study was to compare their current fluxes and positions with previous observations to search for evidence of variability. Most of the observed sources present the morphologies of a thermal core, hosting the protostar and exhibiting no proper motion, and associated non-thermal radio lobes that are characterised by proper motions and located away from the thermal core. Some of the protostars drive jets whose lobes have dissimilar displacement vectors, implying precession of the jets or the presence of multiple jet drivers. The jets of the protostars were found to have proper motions that lie in the range 170$ \leq v \leq $650 kms$^{-1}$, and precessions of periods 40$\leq p \leq$50 years and angles 2$\leq α\leq$10$^\circ$, assuming that their velocities $v=$500 kms$^{-1}$. The core of one of the sources, S255 NIRS 3, which was in outburst at the time of our observations, showed a significant change in flux compared to the other sources. Its spectral index decreased during the outburst, consistent with the model of an expanding gas bubble. Modelling the emission of the outburst as that of a new non-thermal lobe that is emerging from a thermal core whose emission enshrouds that of the lobe also has the potential to account for the increase in flux and a decrease in the spectral index of the source's outburst. △ Less

Submitted 23 December, 2020; originally announced December 2020.

Comments: 16 pages, 15 figures

S2D2: Small-scale Significant substructure DBSCAN Detection I. NESTs detection in 2D star-forming regions

Authors: Marta González, Isabelle Joncour, Anne S. M. Buckner, Zeinhab Khorrami, Estelle Moraux, Stuart L. Lumsden, Paul Clark, René D. Oudmaijer, José Manuel Blanco, Ignacio de la Calle, José María Herrera-Fernandez, Jesús J. Salgado, Luis Valero-Martín, Zoe Torres, Álvaro Hacar, Ana Ulla

Abstract: The spatial and dynamical structure of star-forming regions can help provide insights on stellar formation patterns. The amount of data from current and upcoming surveys calls for robust and objective procedures to detect structure, so the results can be statistically analysed and different regions compared. We provide the community with a tool able to detect the small scale significant structure,… ▽ More The spatial and dynamical structure of star-forming regions can help provide insights on stellar formation patterns. The amount of data from current and upcoming surveys calls for robust and objective procedures to detect structure, so the results can be statistically analysed and different regions compared. We provide the community with a tool able to detect the small scale significant structure, above random expectation, in star-forming regions, which could be the imprint of the stellar formation process. The tool makes use of the one point correlation function and of nearest neighbour statistics to determine the parameters for the DBSCAN algorithm. The procedure successfully detects significant small scale substructures in heterogeneous regions, fulfilling the goals it was designed for, and providing very reliable structures. The analysis of regions close to complete spatial randomness ($Q \in [0.7,0.87]$) shows that, even when some structure is present and recovered, it is hardly distinguishable from spurious detection in homogeneous regions due to projection effects. Interpretation should thus be done with care. For concentrated regions, we detect a main structure surrounded by smaller ones, corresponding to the core plus some Poisson fluctuations around it. We argue that these structures do not correspond to the small compact regions we are looking for. In some realistic cases, a more complete hierarchical, multi-scale analysis would be needed to capture the complexity of the region. We have developed implementations of our procedure, and a catalogue of the NESTs (Nested Elementary STructures) detected by it in four star-forming regions (Taurus, IC 348, Upper Scorpius, and Carina), which are publicly available to the community. Implementations of the 3D, and up to 6D versions of the procedure including proper movements are in progress, and will be provided as future work. △ Less

Submitted 20 November, 2020; originally announced November 2020.

Comments: 24 pages, 21 figures. Accepted for publication in A&A. Abstract abridged

Journal ref: A&A 647, A14 (2021)

Mass segregation and sequential star formation in NGC 2264 revealed by Herschel

Authors: T. Nony, J. -F. Robitaille, F. Motte, M. Gonzalez, I. Joncour, E. Moraux, A. Men'shchikov, P. Didelon, F. Louvet, A. S. M. Buckner, N. Schneider, S. L. Lumsden, S. Bontemps, Y. Pouteau, N. Cunningham, E. Fiorellino, R. Oudmaijer, P. André, B. Thomasson

Abstract: The mass segregation of stellar clusters could be primordial rather than dynamical. Despite the abundance of studies of mass segregation for stellar clusters, those for stellar progenitors are still scarce, so the question on the origin and evolution of mass segregation is still open. Our goal is to characterize the structure of the NGC 2264 molecular cloud and compare the populations of clumps an… ▽ More The mass segregation of stellar clusters could be primordial rather than dynamical. Despite the abundance of studies of mass segregation for stellar clusters, those for stellar progenitors are still scarce, so the question on the origin and evolution of mass segregation is still open. Our goal is to characterize the structure of the NGC 2264 molecular cloud and compare the populations of clumps and young stellar objects (YSOs) in this region whose rich YSO population has shown evidence of sequential star formation. We separated the Herschel column density map of NGC 2264 in three subregions and compared their cloud power spectra using a multiscale segmentation technique. We identified in the whole NGC 2264 cloud a population of 256 clumps with typical sizes of ~0.1 pc and masses ranging from 0.08 Msun to 53 Msun. Although clumps have been detected all over the cloud, the central subregion of NGC 2264 concentrates most of the massive, bound clumps. The local surface density and the mass segregation ratio indeed indicate a strong degree of mass segregation for the 15 most massive clumps, with a median $Σ_6$ three time that of the whole clumps population and $Λ_{MSR}$ about 8. We showed that this cluster of massive clumps is forming within a high-density cloud ridge, itself formed and probably still fed by the high concentration of gas observed on larger scales in the central subregion. The time sequence obtained from the combined study of the clump and YSO populations in NGC 2264 suggests that the star formation started in the northern subregion, that it is now actively developing at the center and will soon start in the southern subregion. Taken together, the cloud structure and the clump and YSO populations in NGC 2264 argue for a dynamical scenario of star formation. △ Less

Submitted 11 November, 2020; originally announced November 2020.

Comments: 26 pages, 19 figures. Accepted for publication in A&A, abstract abridged for arXiv

Journal ref: A&A 645, A94 (2021)

The Spatial Evolution of Young Massive Clusters II. Looking for Imprints of Star Formation in NGC 2264 with Gaia DR2

Authors: Anne S. M. Buckner, Zeinab Khorrami, Marta Gonzalez, Stuart L. Lumsden, Estelle Moraux, Rene D. Oudmaijer, Paul Clark, Isabelle Joncour, Jose Manuel Blanco, Ignacio de la Calle, Alvaro Hacar, Jose M. Herrera-Fernandez, Frederique Motte, Jesus Salgado, Luis Valero-Martin

Abstract: 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 beha… ▽ More 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. △ Less

Submitted 18 March, 2020; v1 submitted 28 February, 2020; originally announced February 2020.

Comments: 13 pages, 4 figures, 5 tables, Accepted for publication in A&A

Journal ref: A&A 636, A80 (2020)

Fragmentation, rotation and outflows in the high-mass star-forming region IRAS 23033+5951. A case study of the IRAM NOEMA large program CORE

Authors: F. Bosco, H. Beuther, A. Ahmadi, J. C. Mottram, R. Kuiper, H. Linz, L. Maud, J. M. Winters, T. Henning, S. Feng, T. Peters, D. Semenov, P. D. Klaassen, P. Schilke, J. S. Urquhart, M. T. Beltrán, S. L. Lumsden, S. Leurini, L. Moscadelli, R. Cesaroni, Á. Sánchez-Monge, A. Palau, R. Pudritz, F. Wyrowski, S. Longmore , et al. (1 additional authors not shown)

Abstract: The formation process of high-mass stars (>8M$_\odot$) is poorly constrained, particularly, the effects of clump fragmentation creating multiple systems and the mechanism of mass accretion onto the cores. We study the fragmentation of dense gas clumps, and trace the circumstellar rotation and outflows by analyzing observations of the high-mass (~500M$_\odot$) star-forming region IRAS 23033+5951. U… ▽ More The formation process of high-mass stars (>8M$_\odot$) is poorly constrained, particularly, the effects of clump fragmentation creating multiple systems and the mechanism of mass accretion onto the cores. We study the fragmentation of dense gas clumps, and trace the circumstellar rotation and outflows by analyzing observations of the high-mass (~500M$_\odot$) star-forming region IRAS 23033+5951. Using the Northern Extended Millimeter Array (NOEMA) in three configurations and the IRAM 30-m single-dish telescope at 220GHz, we probe the gas and dust emission at an angular resolution of ~0.45arcsec, corresponding to 1900au. In the mm continuum emission, we identify a protostellar cluster with at least four mm-sources, where three of them show a significantly higher peak intensity well above a signal-to-noise ratio of 100. Hierarchical fragmentation from large to small spatial scales is discussed. Two fragments are embedded in rotating structures and drive molecular outflows, traced by $^{13}$CO (2-1) emission. The velocity profiles across two of the cores are similar to Keplerian but are missing the highest velocity components close to the center of rotation, which is a common phenomena from observations like these, and other rotation scenarios are not excluded entirely. Position-velocity diagrams suggest protostellar masses of ~6 and 19M$_\sun$. Rotational temperatures from fitting CH$_3$CN ($12_K-11_K$) spectra are used for estimating the gas temperature and by that the disk stability against gravitational fragmentation, utilizing Toomre's $Q$ parameter. [We] identify only one candidate disk to be unstable against gravitational instability caused by axisymmetric perturbations. The dominant sources cover different evolutionary stages within the same maternal gas clump. The appearance of rotation and outflows of the cores are similar to those found in low-mass star-forming regions. △ Less

Submitted 9 July, 2019; originally announced July 2019.

Comments: 23 pages, 16 figures, 5 tables, accepted for publication in A&A

Journal ref: A&A 629, A10 (2019)

IRAS23385+6053: An embedded massive cluster in the making

Authors: R. Cesaroni, H. Beuther, A. Ahmadi, M. T. Beltran, T. Csengeri, R. Galvan-Madrid, C. Gieser, T. Henning, K. G. Johnston, P. D. Klaassen, R. Kuiper, S. Leurini, H. Linz, S. Longmore, S. L. Lumsden, L. T. Maud, L. Moscadelli, J. C. Mottram, A. Palau, T. Peters, R. E. Pudritz, A. Sanchez-Monge, P. Schilke, D. Semenov, S. Suri , et al. (4 additional authors not shown)

Abstract: This study is part of the project ``CORE'', an IRAM/NOEMA large program consisting of observations of the millimeter continuum and molecular line emission towards 20 selected high-mass star forming regions. We focus on IRAS23385+6053, which is believed to be the least evolved source of the CORE sample. The observations were performed at ~1.4 mm and employed three configurations of NOEMA and additi… ▽ More This study is part of the project ``CORE'', an IRAM/NOEMA large program consisting of observations of the millimeter continuum and molecular line emission towards 20 selected high-mass star forming regions. We focus on IRAS23385+6053, which is believed to be the least evolved source of the CORE sample. The observations were performed at ~1.4 mm and employed three configurations of NOEMA and additional single-dish maps, merged with the interferometric data to recover the extended emission. Our correlator setup covered a number of lines from well-known hot core tracers and a few outflow tracers. The angular (~0.45"$-$0.9") and spectral (0.5 km/s) resolutions were sufficient to resolve the clump in IRAS23385+6053 and investigate the existence of large-scale motions due to rotation, infall, or expansion. We find that the clump splits into six distinct cores when observed at sub-arcsecond resolution. These are identified through their 1.4 mm continuum and molecular line emission. We produce maps of the velocity, line width, and rotational temperature from the methanol and methyl cyanide lines, which allow us to investigate the cores and reveal a velocity and temperature gradient in the most massive core. We also find evidence of a bipolar outflow, possibly powered by a low-mass star. We present the tentative detection of a circumstellar self-gravitating disk lying in the most massive core and powering a large-scale outflow previously known in the literature. In our scenario, the star powering the flow is responsible for most of the luminosity of IRAS23385+6053 (~$3000~L_\odot$). The other cores, albeit with masses below the corresponding virial masses, appear to be accreting material from their molecular surroundings and are possibly collapsing or on the verge of collapse. We conclude that we are observing a sample of star-forming cores that is bound to turn into a cluster of massive stars. △ Less

Submitted 27 May, 2019; originally announced May 2019.

Journal ref: A&A 627, A68 (2019)

A search for non-thermal radio emission from jets of massive young stellar objects

Authors: W. O. Obonyo, S. L. Lumsden, M. G. Hoare, S. J. D. Purser, S. E. Kurtz, K. G. Johnston

Abstract: Massive young stellar objects (MYSOs) have recently been shown to drive jets whose particles can interact with either the magnetic fields of the jet or ambient medium to emit non-thermal radiation. We report a search for non-thermal radio emission from a sample of 15 MYSOs to establish the prevalence of the emission in the objects. We used their spectra across the L-, C- and Q-bands along with spe… ▽ More Massive young stellar objects (MYSOs) have recently been shown to drive jets whose particles can interact with either the magnetic fields of the jet or ambient medium to emit non-thermal radiation. We report a search for non-thermal radio emission from a sample of 15 MYSOs to establish the prevalence of the emission in the objects. We used their spectra across the L-, C- and Q-bands along with spectral index maps to characterise their emission. We find that about 50% of the sources show evidence for non-thermal emission with 40% showing clear non-thermal lobes, especially sources of higher bolometric luminosity. The common or IRAS names of the sources that manifest non-thermal lobes are; V645Cyg, IRAS 22134+5834, NGC 7538 IRS 9, IRAS 23262+640, AFGL 402d and AFGL 490. All the central cores of the sources are thermal with corresponding mass-loss rates that lie in the range 3X10^{-7} to 7X10^{-6} solar masses per year. Given the presence of non-thermal lobes in some of the sources and the evidence of non-thermal emission from some spectral index maps, it seems that magnetic fields play a significant role in the jets of massive protostars. Also noted is that some of the sources show evidence of binarity and variability. △ Less

Submitted 11 April, 2019; v1 submitted 9 April, 2019; originally announced April 2019.

Comments: Accepted for publication in MNRAS

MSC Class: Astronomy and Astrophysics

Journal ref: MNRAS, 2019

A multi-scale exploration of a massive young stellar object - a transition disk around G305.20+0.21?

Authors: A. J. Frost, R. D. Oudmaijer, W. J. de Wit, S. L. Lumsden

Abstract: 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 supplem… ▽ More 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. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 16 pages, 8 figures

Journal ref: A&A 625, A44 (2019)

The Spatial Evolution of Young Massive Clusters I. A New Tool to Quantitatively Trace Stellar Clustering

Authors: Anne S. M. Buckner, Zeinab Khorrami, Pouria Khalaj, Stuart L. Lumsden, Isabelle Joncour, Estelle Moraux, Paul Clark, Rene D. Oudmaijer, Jose Manuel Blanco, Ignacio de la Calle, Jose M. Herrera-Fernandez, Frederique Motte, Jesus J. Salgado, Luis Valero-Martin

Abstract: 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 Nebul… ▽ More 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. △ Less

Submitted 18 January, 2019; v1 submitted 7 January, 2019; originally announced January 2019.

Comments: Accepted for publication in A&A; 15 pages, 10 figures, 3 appendices

Journal ref: A&A 622, A184 (2019)

Investigating the temporal domain of massive ionized jets - I. A pilot study

Authors: S. J. D. Purser, S. L. Lumsden, M. G. Hoare, N. Cunningham

Abstract: We present sensitive ($σ<10{\rm μJy \, beam^{-1}}$), radio continuum observations using the Australian telescope compact array (ATCA) at frequencies of $6$ and $9$ GHz towards 4 MYSOs. From a previous, less sensitive work, these objects are known to harbour ionized jets associated with radio lobes, which result from shock processes. In comparison with that work, further emission components are det… ▽ More We present sensitive ($σ<10{\rm μJy \, beam^{-1}}$), radio continuum observations using the Australian telescope compact array (ATCA) at frequencies of $6$ and $9$ GHz towards 4 MYSOs. From a previous, less sensitive work, these objects are known to harbour ionized jets associated with radio lobes, which result from shock processes. In comparison with that work, further emission components are detected towards each MYSO. These include extended, direct, thermal emission from the ionized jet's stream, new radio lobes indicative of shocks close ($<10^5\, {\rm au}$) to the MYSO, 3 radio Herbig-Haro objects separated by up to $3.8 {\rm pc}$ from the jet's launching site and an IR-dark source coincident with CH$_3$OH maser emission. No significant, integrated flux variability is detected towards any jets or shocked lobes, and only one proper motion is observed ($1806\pm596\,{\rm km \, s^{-1}}$ parallel to the jet axis of G310.1420+00.7583A). Evidence for precession is detected in all four MYSOs with precession-periods and angles ranging from $66-15480 \, {\rm yr}$ and $6-36 ^\circ$ respectively. Should precession be the result of the influence from a binary companion, we infer orbital radii of $30-1800\, {\rm au} $. △ Less

Submitted 26 September, 2018; originally announced September 2018.

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 475, Issue 1, p.2-19, 2018

The Coordinated Radio and Infrared Survey for High-mass Star Formation. IV: A new radio selected sample of compact Galactic Planetary Nebulae

Authors: T. Irabor, M. G. Hoare, R. D. Oudmaijer, J. S. Urquhart, S. Kurtz, S. L. Lumsden, C. R. Purcell, A. A. Zijlstra, G. Umana

Abstract: 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 tow… ▽ More 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. △ Less

Submitted 1 August, 2018; v1 submitted 24 July, 2018; originally announced July 2018.

The Coordinated Radio and Infrared Survey for High-Mass Star Formation III. A catalogue of northern ultra-compact H II regions

Authors: I. E. Kalcheva, M. G. Hoare, J. S. Urquhart, S. Kurtz, S. L. Lumsden, C. R. Purcell, A. A. Zijlstra

Abstract: A catalogue of 239 ultra-compact HII regions (UCHIIs) found in the CORNISH survey at 5 GHz and 1.5" resolution in the region $10^{\circ} < l < 65^{\circ}, ~|b| < 1^{\circ}$ is presented. This is the largest complete and well-selected sample of UCHIIs to date and provides the opportunity to explore the global and individual properties of this key state in massive star formation at multiple waveleng… ▽ More A catalogue of 239 ultra-compact HII regions (UCHIIs) found in the CORNISH survey at 5 GHz and 1.5" resolution in the region $10^{\circ} < l < 65^{\circ}, ~|b| < 1^{\circ}$ is presented. This is the largest complete and well-selected sample of UCHIIs to date and provides the opportunity to explore the global and individual properties of this key state in massive star formation at multiple wavelengths. The nature of the candidates was validated, based on observational properties and calculated spectral indices, and the analysis is presented in this work. The physical sizes, luminosities and other physical properties were computed by utilising literature distances or calculating the distances whenever a value was not available. The near- and mid-infrared extended source fluxes were measured and the extinctions towards the UCHIIs were computed. The new results were combined with available data at longer wavelengths and the spectral energy distributions (SEDs) were reconstructed for 177 UCHIIs. The bolometric luminosities obtained from SED fitting are presented. By comparing the radio flux densities to previous observational epochs, we find about 5% of the sources appear to be time variable. This first high-resolution area survey of the Galactic plane shows that the total number of UCHIIs in the Galaxy is ~ 750 - a factor of 3-4 fewer than found in previous large area radio surveys. It will form the basis for future tests of models of massive star formation. △ Less

Submitted 25 March, 2018; originally announced March 2018.

Comments: 30 pages; accepted for publication in A&A on 21.03.2018

Journal ref: A&A 615, A103 (2018)

The protoplanetary system HD 100546 in H$α$ polarized light from SPHERE/ZIMPOL. A bar-like structure across the disk gap?

Authors: I. Mendigutía, R. D. Oudmaijer, A. Garufi, S. L. Lumsden, N. Huélamo, A. Cheetham, W. J. de Wit, B. Norris, F. A. Olguin, P. Tuthill

Abstract: 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$α$ 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 s… ▽ More 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$α$ 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$σ$) 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. △ Less

Submitted 31 October, 2017; originally announced November 2017.

Comments: 8 pages, 7 figures. Accepted in A&A

Journal ref: A&A 608, A104 (2017)

Medium resolution near-infrared spectroscopy of Massive Young Stellar Objects

Authors: R. Pomohaci, R. D. Oudmaijer, S. L. Lumsden, M. G. Hoare, I. Mendigutia

Abstract: 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… ▽ More 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$γ$ 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$γ$ 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$γ$/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. △ Less

Submitted 12 September, 2017; originally announced September 2017.

Comments: 30 pages, 21 figures (with appendices). Accepted for publication in Monthly Notices of the Royal Astronomical Society

Connecting low- and high-mass star formation: the intermediate-mass protostar IRAS 05373+2349 VLA 2

Authors: G. M. Brown, K. G. Johnston, M. G. Hoare, S. L. Lumsden

Abstract: Until recently, there have been few studies of the protostellar evolution of intermediate-mass (IM) stars, which may bridge the low- and high-mass regimes. This paper aims to investigate whether the properties of an IM protostar within the IRAS 05373+2349 embedded cluster are similar to that of low- and/or high-mass protostars. We carried out Very Large Array as well as Combined Array for Research… ▽ More Until recently, there have been few studies of the protostellar evolution of intermediate-mass (IM) stars, which may bridge the low- and high-mass regimes. This paper aims to investigate whether the properties of an IM protostar within the IRAS 05373+2349 embedded cluster are similar to that of low- and/or high-mass protostars. We carried out Very Large Array as well as Combined Array for Research in Millimeter Astronomy continuum and 12CO(J=1-0) observations, which uncover seven radio continuum sources (VLA 1-7). The spectral index of VLA 2, associated with the IM protostar is consistent with an ionised stellar wind or jet. The source VLA 3 is coincident with previously observed H2 emission line objects aligned in the north-south direction (P.A. -20 to -12 deg), which may be either an ionised jet emanating from VLA 2 or (shock-)ionised cavity walls in the large-scale outflow from VLA 2. The position angle between VLA 2 and 3 is slightly misaligned with the large-scale outflow we map at ~5-arcsec resolution in 12CO (P.A. ~30 deg), which in the case of a jet suggests precession. The emission from the mm core associated with VLA 2 is also detected; we estimate its mass to be 12-23 Msun, depending on the contribution from ionised gas. Furthermore, the large-scale outflow has properties intermediate between outflows from low- and high-mass young stars. Therefore, we conclude that the IM protostar within IRAS 05373+2349 is phenomenologically as well as quantitatively intermediate between the low- and high-mass domains. △ Less

Submitted 30 August, 2016; originally announced August 2016.

Comments: 11 pages, 7 figures and 4 tables. Accepted by MNRAS

Locally linear embedding: dimension reduction of massive protostellar spectra

Authors: J. L. Ward, S. L. Lumsden

Abstract: We present the results of the application of locally linear embedding (LLE) to reduce the dimensionality of dereddened and continuum subtracted near-infrared spectra using a combination of models and real spectra of massive protostars selected from the Red MSX Source survey database. A brief comparison is also made with two other dimension reduction techniques; Principal Component Analysis (PCA) a… ▽ More We present the results of the application of locally linear embedding (LLE) to reduce the dimensionality of dereddened and continuum subtracted near-infrared spectra using a combination of models and real spectra of massive protostars selected from the Red MSX Source survey database. A brief comparison is also made with two other dimension reduction techniques; Principal Component Analysis (PCA) and Isomap using the same set of spectra as well as a more advanced form of LLE, Hessian locally linear embedding. We find that whilst LLE certainly has its limitations, it significantly outperforms both PCA and Isomap in classification of spectra based on the presence/absence of emission lines and provides a valuable tool for classification and analysis of large spectral data sets. △ Less

Submitted 22 June, 2016; originally announced June 2016.

Comments: 8 pages, 7 figures. Accepted for publication in MNRAS 2016 June 21

A search for ionized jets towards massive young stellar objects

Authors: S. J. D. Purser, S. L. Lumsden, M. G. Hoare, J. S. Urquhart, N. Cunningham, C. R. Purcell, K. J. Brooks, G. Garay, A. E. Gúzman, M. A. Voronkov

Abstract: Radio continuum observations using the Australia telescope compact array at 5.5, 9.0, 17.0 and 22.8 GHz have detected free-free emission associated with 45 of 49 massive young stellar objects and HII regions. Of these, 26 sources are classified as ionized jets (12 of which are candidates), 2 as ambiguous jets or disc winds, 1 as a disc-wind, 14 as HII regions and 2 were unable to be categorised. C… ▽ More Radio continuum observations using the Australia telescope compact array at 5.5, 9.0, 17.0 and 22.8 GHz have detected free-free emission associated with 45 of 49 massive young stellar objects and HII regions. Of these, 26 sources are classified as ionized jets (12 of which are candidates), 2 as ambiguous jets or disc winds, 1 as a disc-wind, 14 as HII regions and 2 were unable to be categorised. Classification as ionized jets is based upon morphology, radio flux and spectral index, in conjunction with previous observational results at other wavelengths. Radio-luminosity and momentum are found to scale with bolometric luminosity in the same way as low-mass jets, indicating a common mechanism for jet production across all masses. In 13 of the jets, we see associated non-thermal/optically-thin lobes resulting from shocks either internal to the jet and/or at working surfaces. Ten jets display non-thermal (synchrotron emission) spectra in their lobes, with an average spectral index of -0.55 consistent with Fermi acceleration in shocks. This shows that magnetic fields are present, in agreement with models of jet formation incorporating magnetic fields. Since the production of collimated radio jets is associated with accretion processes, the results presented in this paper support the picture of disc-mediated accretion for the formation of massive stars with an upper-limit on the jet phase lasting approximately $6.5 \times 10^4 yr$. Typical mass loss rates in the jet are found to be $1.4 \times 10^{-5} M_\odot yr^{-1}$ with associated momentum rates of the order $(1-2) \times 10^{-2} M_\odot km s^{-1} yr^{-1}$. △ Less

Submitted 4 May, 2016; originally announced May 2016.

Submillimeter Array Observations of NGC 2264-C: Molecular Outflows and Driving Sources

Authors: Nichol Cunningham, Stuart L. Lumsden, Claudia J. Cyganowski, Luke T. Maud, Cormac Purcell

Abstract: We present 1.3mm Submillimeter Array (SMA) observations at $\sim$3$^{\prime\prime}$ resolution towards the brightest section of the intermediate/massive star forming cluster NGC 2264-C. The millimetre continuum emission reveals ten 1.3mm continuum peaks, of which four are new detections. The observed frequency range includes the known molecular jet/outflow tracer SiO (5-4), thus providing the firs… ▽ More We present 1.3mm Submillimeter Array (SMA) observations at $\sim$3$^{\prime\prime}$ resolution towards the brightest section of the intermediate/massive star forming cluster NGC 2264-C. The millimetre continuum emission reveals ten 1.3mm continuum peaks, of which four are new detections. The observed frequency range includes the known molecular jet/outflow tracer SiO (5-4), thus providing the first high resolution observations of SiO towards NGC 2264-C. We also detect molecular lines of twelve additional species towards this region, including CH$_3$CN, CH$_3$OH, SO, H$_2$CO, DCN, HC$_3$N, and $^{12}$CO. The SiO (5-4) emission reveals the presence of two collimated, high velocity (up to 30kms$^{-1}$ with respect to the systemic velocity) bi-polar outflows in NGC 2264-C. In addition, the outflows are traced by emission from $^{12}$CO, SO, H$_2$CO, and CH$_3$OH. We find an evolutionary spread between cores residing in the same parent cloud. The two unambiguous outflows are driven by the brightest mm continuum cores, which are IR-dark, molecular line weak, and likely the youngest cores in the region. Furthermore, towards the RMS source AFGL 989-IRS1, the IR-bright and most evolved source in NGC 2264-C, we observe no molecular outflow emission. A molecular line rich ridge feature, with no obvious directly associated continuum source, lies on the edge of a low density cavity and may be formed from a wind driven by AFGL 989-IRS1. In addition, 229GHz class I maser emission is detected towards this feature. △ Less

Submitted 18 February, 2016; originally announced February 2016.

Comments: 29 pages, 12 figures

The JCMT Plane Survey: early results from the l = 30 degree field

Authors: T. J. T. Moore, R. Plume, M. A. Thompson, H. Parsons, J. S. Urquhart, D. J. Eden, J. T. Dempsey, L. K. Morgan, H. S. Thomas, J. Buckle, C. M. Brunt, H. Butner, D. Carretero, A. Chrysostomou, H. M. deVilliers, M. Fich, M. G. Hoare, G. Manser, J. C. Mottram, C. Natario, F. Olguin, N. Peretto, D. Polychroni, R. O. Redman, A. J. Rigby , et al. (33 additional authors not shown)

Abstract: We present early results from the JCMT Plane Survey (JPS), which has surveyed the northern inner Galactic plane between longitudes l=7 and l=63 degrees in the 850-μm continuum with SCUBA-2, as part of the James Clerk Maxwell Telescope Legacy Survey programme. Data from the l=30 degree survey region, which contains the massive star-forming regions W43 and G29.96, are analysed after approximately 40… ▽ More We present early results from the JCMT Plane Survey (JPS), which has surveyed the northern inner Galactic plane between longitudes l=7 and l=63 degrees in the 850-μm continuum with SCUBA-2, as part of the James Clerk Maxwell Telescope Legacy Survey programme. Data from the l=30 degree survey region, which contains the massive star-forming regions W43 and G29.96, are analysed after approximately 40% of the observations had been completed. The pixel-to-pixel noise is found to be 19 mJy/beam, after a smooth over the beam area, and the projected equivalent noise levels in the final survey are expected to be around 10 mJy/beam. An initial extraction of compact sources was performed using the FellWalker method resulting in the detection of 1029 sources above a 5-σ surface-brightness threshold. The completeness limits in these data are estimated to be around 0.2 Jy/beam (peak flux density) and 0.8 Jy (integrated flux density) and are therefore probably already dominated by source confusion in this relatively crowded section of the survey. The flux densities of extracted compact sources are consistent with those of matching detections in the shallower ATLASGAL survey. We analyse the virial and evolutionary state of the detected clumps in the W43 star-forming complex and find that they appear younger than the Galactic-plane average. △ Less

Submitted 1 September, 2015; originally announced September 2015.

Comments: Accepted by MNRAS

A distance limited sample of massive molecular outflows

Authors: L. T. Maud, T. J. T. Moore, S. L. Lumsden, J. C. Mottram, J. S. Urquhart, M. G. Hoare

Abstract: We have observed 99 mid-infrared-bright, massive young stellar objects and compact HII regions drawn from the Red MSX source (RMS) survey in the J=3$-$2 transition of $^{12}$CO and $^{13}$CO, using the James Clerk Maxwell Telescope. 89 targets are within 6 kpc of the Sun, covering a representative range of luminosities and core masses. These constitute a relatively unbiased sample of bipolar molec… ▽ More We have observed 99 mid-infrared-bright, massive young stellar objects and compact HII regions drawn from the Red MSX source (RMS) survey in the J=3$-$2 transition of $^{12}$CO and $^{13}$CO, using the James Clerk Maxwell Telescope. 89 targets are within 6 kpc of the Sun, covering a representative range of luminosities and core masses. These constitute a relatively unbiased sample of bipolar molecular outflows associated with massive star formation. Of these, 59, 17 and 13 sources (66, 19 and 15 percent) are found to have outflows, show some evidence of outflow, and have no evidence of outflow, respectively. The time-dependent parameters of the high-velocity molecular flows are calculated using a spatially variable dynamic timescale. The canonical correlations between the outflow parameters and source luminosity are recovered and shown to scale with those of low-mass sources. For coeval star formation we find the scaling is consistent with all the protostars in an embedded cluster providing the outflow force, with massive stars up to $\sim$30 M$_{\odot}$ generating outflows. Taken at face value, the results support the model of a scaled-up version of the accretion-related outflow-generation mechanism associated with discs and jets in low-mass objects with time-averaged accretion rates of $\sim$10$^{-3}$ M$_{\odot}$ yr$^{-1}$ onto the cores. However, we also suggest an alternative model, in which the molecular outflow dynamics are dominated by the entrained mass and are unrelated to the details of the acceleration mechanism. We find no evidence that outflows contribute significantly to the turbulent kinetic energy of the surrounding dense cores. △ Less

Submitted 1 September, 2015; originally announced September 2015.

Journal ref: 2015MNRAS.453..645M

A distance limited sample of massive star forming cores from the RMS survey

Authors: L. T. Maud, S. L. Lumsden, T. J. T. Moore, J. C. Mottram, J. S. Urquhart, A. Cicchini

Abstract: We analyse C$^{18}$O ($J=3-$2) data from a sample of 99 infrared-bright massive young stellar objects (MYSOs) and compact HII regions that were identified as potential molecular-outflow sources in the Red MSX source (RMS) survey. We extract a distance limited (D $<$ 6 kpc) sample shown to be representative of star formation covering the transition between the source types. At the spatial resolutio… ▽ More We analyse C$^{18}$O ($J=3-$2) data from a sample of 99 infrared-bright massive young stellar objects (MYSOs) and compact HII regions that were identified as potential molecular-outflow sources in the Red MSX source (RMS) survey. We extract a distance limited (D $<$ 6 kpc) sample shown to be representative of star formation covering the transition between the source types. At the spatial resolution probed, Larson-like relationships are found for these cores, though the alternative explanation, that Larson's relations arise where surface-density-limited samples are considered, is also consistent with our data. There are no significant differences found between source properties for the MYSOs and HII regions, suggesting that the core properties are established prior to the formation of massive stars, which subsequently have little impact at the later evolutionary stages investigated. There is a strong correlation between dust-continuum and C$^{18}$O-gas masses, supporting the interpretation that both trace the same material in these IR-bright sources. A clear linear relationship is seen between the independently established core masses and luminosities. The position of MYSOs and compact HII regions in the mass-luminosity plane is consistent with the luminosity expected a cluster of protostars when using a $\sim$40 percent star-formation efficiency and indicates that they are at a similar evolutionary stage, near the end of the accretion phase. △ Less

Submitted 1 September, 2015; originally announced September 2015.

Journal ref: 2015MNRAS.452..637M

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

Authors: I. Mendigutía, R. D. Oudmaijer, E. Rigliaco, J. R. Fairlamb, N. Calvet, J. Muzerolle, N. Cunningham, S. L. Lumsden

Abstract: 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… ▽ More 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. △ Less

Submitted 9 July, 2015; originally announced July 2015.

Comments: 9 pages, 5 figures, 2 tables. Accepted in MNRAS

The RMS Survey: Ammonia mapping of the environment of massive young stellar objects

Authors: J. S. Urquhart, C. C. Figura, T. J. T. Moore, T. Csengeri, S. L. Lumsden, T. Pillai, M. A. Thompson, D. J. Eden, L. K. Morgan

Abstract: We present the results of ammonia observations towards 66 massive star forming regions identified by the Red MSX source survey. We have used the Green Bank Telescope and the K-band focal plane array to map the ammonia NH3 (1,1) and (2,2) inversion emission at a resolution of 30 arcsec in 8 arcmin regions towards the positions of embedded massive star formation. We have identified a total of 115 di… ▽ More We present the results of ammonia observations towards 66 massive star forming regions identified by the Red MSX source survey. We have used the Green Bank Telescope and the K-band focal plane array to map the ammonia NH3 (1,1) and (2,2) inversion emission at a resolution of 30 arcsec in 8 arcmin regions towards the positions of embedded massive star formation. We have identified a total of 115 distinct clumps, approximately two-thirds of which are associated with an embedded massive young stellar object or compact HII region, while the others are classified as quiescent. There is a strong spatial correlation between the peak NH3 emission and the presence of embedded objects. We derive the spatial distribution of the kinetic gas temperatures, line widths, and NH$_3$ column densities from these maps, and by combining these data with dust emission maps we estimate clump masses, H$_2$ column densities and ammonia abundances. The clumps have typical masses of ~1000 Msun and radii ~0.5 pc, line widths of ~2 km/s and kinetic temperatures of ~16-20 K. We find no significant difference between the sizes and masses of the star forming and quiescent subsamples; however, the distribution maps reveal the presence of temperature and line width gradients peaking towards the centre for the star forming clumps while the quiescent clumps show relatively uniform temperatures and line widths throughout. Virial analysis suggests that the vast majority of clumps are gravitationally bound and are likely to be in a state of global free fall in the absence of strong magnetic fields. The similarities between the properties of the two subsamples suggest that the quiescent clumps are also likely to form massive stars in the future, and therefore provide a excellent opportunity to study the initial conditions of massive pre-stellar and protostellar clumps. △ Less

Submitted 8 July, 2015; originally announced July 2015.

Comments: Accepted for publication by MNRAS. 25 pages, 21 figures and 6 tables. This paper also contains 100 pages of additional online images and tables that are not included here. A full resolution version of the paper including all of the online material can be found here http://rms.leeds.ac.uk/cgi-bin/public/RMS_PUBLICATIONS.cgi

ATLASGAL --- towards a complete sample of massive star forming clumps

Authors: J. S. Urquhart, T. J. T. Moore, T. Csengeri, F. Wyrowski, F. Schuller, M. G. Hoare, S. L. Lumsden, J. C. Mottram, M. A. Thompson, K. M. Menten, C. M. Walmsley, L. Bronfman, S. Pfalzner, C. König, M. Wienen

Abstract: By matching infrared-selected, massive young stellar objects (MYSOs) and compact HII regions in the RMS survey to massive clumps found in the submillimetre ATLASGAL survey, we have identified ~1000 embedded young massive stars between 280\degr < $\ell$ < 350\degr and 10degr < $\ell$ < 60\degr with |b|<1.5degr. Combined with an existing sample of radio-selected methanol masers and compact HII regio… ▽ More By matching infrared-selected, massive young stellar objects (MYSOs) and compact HII regions in the RMS survey to massive clumps found in the submillimetre ATLASGAL survey, we have identified ~1000 embedded young massive stars between 280\degr < $\ell$ < 350\degr and 10degr < $\ell$ < 60\degr with |b|<1.5degr. Combined with an existing sample of radio-selected methanol masers and compact HII regions, the result is a catalogue of ~1700 massive stars embedded within ~1300 clumps located across the inner Galaxy, containing three observationally distinct subsamples, methanol-maser, MYSO and HII-region associations, covering the most important tracers of massive star formation, thought to represent key stages of evolution. We find that massive star formation is strongly correlated with the regions of highest column density in spherical, centrally condensed clumps. We find no significant differences between the three samples in clump structure or the relative location of the embedded stars, which suggests that the structure of a clump is set before the onset of star formation, and changes little as the embedded object evolves towards the main sequence. There is a strong linear correlation between clump mass and bolometric luminosity, with the most massive stars forming in the most massive clumps. We find that the MYSO and HII-region subsamples are likely to cover a similar range of evolutionary stages and that the majority are near the end of their main accretion phase. We find few infrared-bright MYSOs associated with the most massive clumps, probably due to very short pre-main sequence lifetimes in the most luminous sources. △ Less

Submitted 19 June, 2014; originally announced June 2014.

Comments: Accepted for publication in MNRAS: 33 pages, 24 figures and 9 tables. Full versions of Fig. 3 and Tables 1, 3, 6 and 7 will only be available in electronic form via CDS. The resolution of Figs. 1, 5, 7 and 21 have been reduced but a higher resolution version of the paper is available upon request

The RMS Survey: Galactic distribution of massive star formation

Authors: J. S. Urquhart, C. C. Figura, T. J. T., Moore, M. G. Hoare, S. L. Lumsden, J. C. Mottram, M. A. Thompson, R. D. Oudmaijer

Abstract: 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 di… ▽ More 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. △ Less

Submitted 17 October, 2013; originally announced October 2013.

Comments: Accepted for publication in MNRAS. Paper consists of 17 pages, four table and 11 figures. The full versions of Tables 1,2 and 3 and Figures 1 and 4 will only be available from CDS

Near-infrared integral field spectroscopy of Massive Young Stellar Objects

Authors: K. Murakawa, S. L. Lumsden, R. D. Oudmaijer, B. Davies, H. E. Wheelright, M. G. Hoare, J. D. Ilee

Abstract: 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… ▽ More 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. △ Less

Submitted 24 September, 2013; originally announced September 2013.

Comments: 17 pages, 10 figures, accepted by MNRAS

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

Authors: S. L. Lumsden, M. G. Hoare, J. S. Urquhart, R. D. Oudmaijer, B. Davies, J. C. Mottram, H. D. B. Cooper, T. J. T. Moore

Abstract: 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… ▽ More 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. △ Less

Submitted 2 August, 2013; v1 submitted 1 August, 2013; originally announced August 2013.

Comments: Accepted for publication in ApJS. The full version of Table 1 will be made available electronically, but the most up-to-date version of our database will always be that accessible at http://rms.leeds.ac.uk

ATLASGAL --- properties of compact HII regions and their natal clumps

Authors: J. S. Urquhart, M. A. Thompson, T. J. T. Moore, C. R. Purcell, M. G. Hoare, F. Schuller, F. Wyrowski, T. Csengeri, K. M. Menten, S. L. Lumsden, S. Kurtz, C. M. Walmsley, L. Bronfman, L. K. Morgan, D. J. Eden, D. Russeil

Abstract: We present a complete sample of molecular clumps containing compact and ultra-compact (UC) HII regions between \ell=10\degr and 60\degr\ and $|b|<1\degr, identified by combining the the ATLASGAL submm and CORNISH radio continuum surveys with visual examination of archival infrared data. Our sample is complete to optically thin, compact and UCHII regions driven by a zero age main sequence star of s… ▽ More We present a complete sample of molecular clumps containing compact and ultra-compact (UC) HII regions between \ell=10\degr and 60\degr\ and $|b|<1\degr, identified by combining the the ATLASGAL submm and CORNISH radio continuum surveys with visual examination of archival infrared data. Our sample is complete to optically thin, compact and UCHII regions driven by a zero age main sequence star of spectral type B0 or earlier embedded within a 1,000 Msun clump. In total we identify 213 compact and UCHII regions, associated with 170 clumps. Unambiguous kinematic distances are derived for these clumps and used to estimate their masses and physical sizes, as well as the Lyman continuum fluxes and sizes of their embedded HII regions. We find a clear lower envelope for the surface density of molecular clumps hosting massive star formation of 0.05 g cm^{-2}, which is consistent with a similar sample of clumps associated with 6.7 GHz masers. The mass of the most massive embedded stars is closely correlated with the mass of their natal clump. Young B stars appear to be significantly more luminous in the ultraviolet than predicted by current stellar atmosphere models. The properties of clumps associated with compact and UCHII regions are very similar to those associated with 6.7 GHz methanol masers and we speculate that there is little evolution in the structure of the molecular clumps between these two phases. Finally, we identify a significant peak in the surface density of compact and UCHII regions associated with the W49A star-forming complex, noting that this complex is truly one of the most massive and intense regions of star formation in the Galaxy. △ Less

Submitted 15 July, 2013; originally announced July 2013.

Comments: Accepted for publication in MNRAS. Paper consists of 29 pages, 26 figures and 4 tables. The resolution of Figures 21 and 25 has been reduced to comply with file size restrictions. Full versions of Figures 11 and 12 and Tables 3 and 4 will only be available in the electronic version of the journal and through CDS

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

Authors: H. D. B. Cooper, S. L. Lumsden, R. D. Oudmaijer, M. G. Hoare, A. J. Clarke, J. S. Urquhart, J. C. Mottram, T. J. T. Moore, B. Davies

Abstract: 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… ▽ More 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. △ Less

Submitted 17 January, 2013; originally announced January 2013.

Comments: Accepted for publication in MNRAS

CO bandhead emission of massive young stellar objects: determining disc properties

Authors: J. D. Ilee, H. E. Wheelwright, R. D. Oudmaijer, W. J. de Wit, L. T. Maud, M. G. Hoare, S. L. Lumsden, T. J. T. Moore, J. S. Urquhart, J. C. Mottram

Abstract: 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 assumpt… ▽ More 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. △ Less

Submitted 3 December, 2012; originally announced December 2012.

Comments: 14 pages, 8 figures, accepted for publication in MNRAS

The Co-ordinated Radio and Infrared Survey for High-Mass Star Formation - II. Source Catalogue

Authors: C. R. Purcell, M. G. Hoare, W. D. Cotton, S. L. Lumsden, J. S. Urquhart, C. Chandler, E. B. Churchwell, P. Diamond, S. M. Dougherty, R. P. Fender, G. Fuller, S. T. Garrington, T. M. Gledhill, P. F. Goldsmith, L. Hindson, J. M. Jackson, S. E. Kurtz, J. Marti, T. J. T. Moore, L. G. Mundy, T. W. B. Muxlow, R. D. Oudmaijer, J. D. Pandian, J. M. Paredes, D. S. Shepherd , et al. (5 additional authors not shown)

Abstract: 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 Stoke… ▽ More 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 △ Less

Submitted 3 December, 2012; v1 submitted 29 November, 2012; originally announced November 2012.

Comments: 25 pages, 22 figures. To appear in ApJ. Sup

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

Authors: M. G. Hoare, C. R. Purcell, E. B. Churchwell, P. Diamond, W. D. Cotton, C. J. Chandler, S. Smethurst, S. E. Kurtz, L. G. Mundy, S. M. Dougherty, R. P. Fender, G. A. Fuller, J. M. Jackson, S. T. Garrington, T. R. Gledhill, P. F. Goldsmith, S. L. Lumsden, J. Martí, T. J. T. Moore, T. W. B. Muxlow, R. D. Oudmaijer, J. D. Pandian, J. M. Paredes, D. S. Shepherd, R. E. Spencer , et al. (4 additional authors not shown)

Abstract: 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 emp… ▽ More 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. △ Less

Submitted 16 August, 2012; originally announced August 2012.

Comments: 17 pages, 14 figures. Accepted for publication in PASP (16th-August-2012)

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

Authors: S. L. Lumsden, H. E. Wheelwright, M. G. Hoare, R. D. Oudmaijer, J. E. Drew

Abstract: 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 broad… ▽ More 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. △ Less

Submitted 10 May, 2012; originally announced May 2012.

Comments: 20 pages, 7 figures, accepted for publication in MNRAS

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

Authors: H. E. Wheelwright, W. J. de Wit, R. D. Oudmaijer, M. G. Hoare, S. L. Lumsden, T. Fujiyoshi, J. L. Close

Abstract: 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… ▽ More 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. △ Less

Submitted 20 February, 2012; originally announced February 2012.

Comments: Accepted for publication in A&A. The images in this version have been compressed. A high resolution version is available on request

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

Authors: J. S. Urquhart, M. G. Hoare, S. L. Lumsden, R. D. Oudmaijer, T. J. T. Moore, J. C. Mottram, H. D. B. Cooper, M. Mottram, H. C. Rogers

Abstract: 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… ▽ More 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. △ Less

Submitted 7 November, 2011; originally announced November 2011.

Comments: Accepted by MNRAS. This paper consists of 15 pages and contains 10 figures and 5 tables

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

Authors: J. S. Urquhart, L. K. Morgan, C. C. Figura, T. J. T. Moore, S. L. Lumsden, M. G. Hoare, R. D. Oudmaijer, J. C. Mottram, B. Davies, M. K. Dunham

Abstract: 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… ▽ More 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. △ Less

Submitted 20 July, 2011; originally announced July 2011.

Comments: 17 pages and 17 figures and 8 tables. Tables\,2 and 5 and full versions of Figs. 3 and 7 are only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.125.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/

The RMS Survey: Critical Tests of Accretion Models for the Formation of Massive Stars

Authors: Ben Davies, Melvin G. Hoare, Stuart L. Lumsden, Takashi Hosokawa, Rene D. Oudmaijer, James S. Urquhart, Joseph C. Mottram, Joseph Stead

Abstract: 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… ▽ More 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. △ Less

Submitted 19 May, 2011; originally announced May 2011.

Comments: 22 pages, 24 figures. Accepted for publication in MNRAS

The Luminosity Functions and Timescales of MYSOs and Compact HII regions

Authors: Joseph C. Mottram, Melvin G. Hoare, Ben Davies, Stuart L. Lumsden, Rene D. Oudmaijer, James S. Urquhart, Toby J. T. Moore, Heather D. B. Cooper, Joseph J. Stead

Abstract: 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… ▽ More 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. △ Less

Submitted 23 February, 2011; originally announced February 2011.

Comments: 16 pages in pre-print format, 4 figures, 1 table, accepted to ApJL

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

Authors: W. J. de Wit, M. G. Hoare, R. D. Oudmaijer, D. E. A Nuernberger, H. E. Wheelwright, S. L. Lumsden

Abstract: 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" an… ▽ More 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. △ Less

Submitted 1 December, 2010; originally announced December 2010.

Comments: Accepted for publication in A&A

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

Authors: Joseph C. Mottram, M. G. Hoare, J. S. Urquhart, S. L. Lumsden, R. D. Oudmaijer, T. P. Robitaille, T. J. T. Moore, B. Davies, J. Stead

Abstract: 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: W… ▽ More 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{}. △ Less

Submitted 10 September, 2010; v1 submitted 9 September, 2010; originally announced September 2010.

Comments: 12 pages, 12 figures, 3 tables, accepted to A&A. The full versions of tables 1 and 2 will be available via the CDS upon publication

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

Authors: J. S. Urquhart, T. J. T. Moore, M. G. Hoare, S. L. Lumsden, R. D. Oudmaijer, J. M. Rathborne, J. C. Mottram, B. Davies, J. J. Stead

Abstract: 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 1… ▽ More 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. △ Less

Submitted 10 September, 2010; v1 submitted 18 August, 2010; originally announced August 2010.

Comments: Accepted for publication by MNRAS. Paper consists of 15 pages including 12 figures and four tables. Full versions of Tables 2 and 3 will only be available online. The resolution of Figure 9 has been reduced - a full resolution version of the paper can be download from here: http://www.ast.leeds.ac.uk/cgi-bin/RMS/RMS_PUBLICATIONS.cgi

Probing discs around massive young stellar objects with CO first overtone emission

Authors: H. E. Wheelwright, R. D. Oudmaijer, W. J. de Wit, M. G. Hoare, S. L. Lumsden, J. S. Urquhart

Abstract: 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… ▽ More 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. △ Less

Submitted 19 July, 2010; originally announced July 2010.

Comments: Accepted for publication in MNRAS