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.

We present a detailed observational and modeling study of the hot core VLA 3 in the high-mass star-forming region AFGL 2591, which is a target region of the NOrthern Extended Millimeter Array (NOEMA) large program CORE. Using NOEMA observations at 1.37 mm with an angular resolution of ~0."42 (1 400 au at 3.33 kpc), we derived the physical and chemical structure of the source. We modeled the observed molecular abundances with the chemical evolution code MUSCLE (MUlti Stage ChemicaL codE). Results. With the kinetic temperature tracers CH3CN and H2CO we observe a temperature distribution with a power-law index of q = 0.41+-0.08. Using the visibilities of the continuum emission we derive a density structure with a power-law index of p = 1.7+-0.1. The hot core spectra reveal high molecular abundances and a rich diversity in complex molecules. The majority of the molecules have an asymmetric spatial distribution around the forming protostar(s), which indicates a complex physical structure on scales < 1 400 au. Using MUSCLE, we are able to explain the observed molecular abundance of 10 out of 14 modeled species at an estimated hot core chemical age of ~21 100 years. In contrast to the observational analysis, our chemical modeling predicts a lower density power-law index of p < 1.4. Reasons for this discrepancy are discussed. Conclusions. Combining high spatial resolution observations with detailed chemical modeling allows us to derive a concise picture of the physical and chemical structure of the famous AFGL 2591 hot core. The next steps are to conduct a similar analysis for the whole CORE sample, and then use this analysis to constrain the chemical diversity in high-mass star formation to a much greater depth.

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.

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.

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.

The past two decades have seen extensive surveys of the far-infrared to submillimeter continuum emission in the plane of our Galaxy. We line out prospects for the coming decade for corresponding molecular and atomic line surveys which are needed to fully understand the formation of the dense structures that give birth to clusters and stars out of the diffuse interstellar medium. We propose to work towards Galaxy wide surveys in mid-J CO lines to trace shocks from colliding clouds, Galaxy-wide surveys for atomic Carbon lines in order to get a detailed understanding of the relation of atomic and molecular gas in clouds, and to perform extensive surveys of the structure of the dense parts of molecular clouds to understand the importance of filaments/fibers over the full range of Galactic environments and to study how dense cloud cores are formed from the filaments. This work will require a large (50m) Single Dish submillimeter telescope equipped with massively multipixel spectrometer arrays, such as envisaged by the AtLAST project.

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.

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

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

We present high angular resolution 0.2 arcsec continuum and molecular emission line Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of G17.64+0.16 in Band 6 (220GHz) taken as part of a campaign in search of circumstellar discs around (proto)-O-stars. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. We detect SiO (5-4) emission that is marginally resolved and elongated in a direction perpendicular to the large-scale outflow seen in the 13CO (2-1) line using the main ALMA array in conjunction with the Atacama Compact Array (ACA). Morphologically, the SiO appears to represent a disc-like structure. Using parametric models we show that the position-velocity profile of the SiO is consistent with the Keplerian rotation of a disc around an object between 10-30Mo in mass, only if there is also radial expansion from a separate structure. The radial motion component can be interpreted as a disc wind from the disc surface. Models with a central stellar object mass between 20 and 30Mo are the most consistent with the stellar luminosity (100000 Lo) and indicative of an O-type star. The H30a millimetre recombination line (231.9GHz) is also detected, but spatially unresolved, and is indicative of a very compact, hot, ionised region co-spatial with the dust continuum core. Accounting for all observables, we suggest that G17.64 is consistent with a O-type young stellar object in the final stages of protostellar assembly, driving a wind, but that has not yet developed into a compact HII region. The existance and detection of the disc in G17.64 is likely related to its isolated and possibly more evolved nature, traits which may underpin discs in similar sources.

We present sensitive ($\sigma<10{\rm \mu 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} $.

The fragmentation mode of high-mass molecular clumps and the properties of the central rotating structures surrounding the most luminous objects have yet to be comprehensively characterised. Using the IRAM NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30-m telescope, the CORE survey has obtained high-resolution observations of 20 well-known highly luminous star-forming regions in the 1.37 mm wavelength regime in both line and dust continuum emission. We present the spectral line setup of the CORE survey and a case study for W3(H2O). At ~0.35" (700 AU at 2 kpc) resolution, the W3(H2O) clump fragments into two cores (West and East), separated by ~2300 AU. Velocity shifts of a few km/s are observed in the dense-gas tracer, CH3CN, across both cores, consistent with rotation and perpendicular to the directions of two bipolar outflows, one emanating from each core. The kinematics of the rotating structure about W3(H2O) W shows signs of differential rotation of material, possibly in a disk-like object. The observed rotational signature around W3(H2O) E may be due to a disk-like object, an unresolved binary (or multiple) system, or a combination of both. We fit the emission of CH3CN (12-11) K = 4-6 and derive a gas temperature map with a median temperature of ~165 K across W3(H2O). We create a Toomre Q map to study the stability of the rotating structures against gravitational instability. The rotating structures appear to be Toomre unstable close to their outer boundaries, with a possibility of further fragmentation in the differentially-rotating core W3(H2O) W. Rapid cooling in the Toomre-unstable regions supports the fragmentation scenario. Combining millimeter dust continuum and spectral line data toward the famous high-mass star-forming region W3(H2O), we identify core fragmentation on large scales, and indications for possible disk fragmentation on smaller spatial scales.

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.

Aims: We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during high-mass star formation on spatial scales of individual cores. Methods: Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30m telescope, we have observed in the IRAM large program CORE the 1.37mm continuum and spectral line emission at high angular resolution (~0.4'') for a sample of 20 well-known high-mass star-forming regions with distances below 5.5kpc and luminosities larger than 10^4Lsun. Results: We present the overall survey scope, the selected sample, the observational setup and the main goals of CORE. Scientifically, we concentrate on the mm continuum emission on scales on the order of 1000AU. We detect strong mm continuum emission from all regions, mostly due to the emission from cold dust. The fragmentation properties of the sample are diverse. We see extremes where some regions are dominated by a single high-mass core whereas others fragment into as many as 20 cores. A minimum-spanning-tree analysis finds fragmentation at scales on the order of the thermal Jeans length or smaller suggesting that turbulent fragmentation is less important than thermal gravitational fragmentation. The diversity of highly fragmented versus singular regions can be explained by varying initial density structures and/or different initial magnetic field strengths. Conclusions: The smallest observed separations between cores are found around the angular resolution limit which indicates that further fragmentation likely takes place on even smaller spatial scales. The CORE project with its numerous spectral line detections will address a diverse set of important physical and chemical questions in the field of high-mass star formation.

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.

We present the results of an outflow and infall survey towards a distance limited sample of 31 massive star forming regions drawn from the RMS survey. The presence of young, active outflows is identified from SiO (8-7) emission and the infall dynamics are explored using HCO$^+$/H$^{13}$CO$^+$ (4-3) emission. We investigate if the infall and outflow parameters vary with source properties, exploring whether regions hosting potentially young active outflows show similarities or differences with regions harbouring more evolved, possibly momentum driven, "fossil" outflows. SiO emission is detected towards approximately 46% of the sources. When considering sources with and without an SiO detection (i.e. potentially active and fossil outflows respectively), only the $^{12}$CO outflow velocity shows a significant difference between samples, indicating SiO is more prevalent towards sources with higher outflow velocities. Furthermore, we find the SiO luminosity increases as a function of the Herschel 70$\mu$m to WISE 22$\mu$m flux ratio, suggesting the production of SiO is prevalent in younger, more embedded regions. Similarly, we find tentative evidence that sources with an SiO detection have a smaller bolometric luminosity-to-mass ratio, indicating SiO (8-7) emission is associated with potentially younger regions. We do not find a prevalence towards sources displaying signatures of infall in our sample. However, the higher energy HCO$^+$ transitions may not be the best suited tracer of infall at this spatial resolution in these regions.

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

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

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.

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.

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

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.

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-{\mu}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-{\sigma} 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.

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.

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.

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.

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.

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.

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

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

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.

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.

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.

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.

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

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.

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

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.

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.

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.

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

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

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

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.

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

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.

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

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

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

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

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

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

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

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

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

We have surveyed a ~0.9-square-degree area of the W3 giant molecular cloud and star-forming region in the 850-micron continuum, using the SCUBA bolometer array on the James Clerk Maxwell Telescope. A complete sample of 316 dense clumps was detected with a mass range from around 13 to 2500 Msun. Part of the W3 GMC is subject to an interaction with the HII region and fast stellar winds generated by the nearby W4 OB association. We find that the fraction of total gas mass in dense, 850-micron traced structures is significantly altered by this interaction, being around 5% to 13% in the undisturbed cloud but ~25 - 37% in the feedback-affected region. The mass distribution in the detected clump sample depends somewhat on assumptions of dust temperature and is not a simple, single power law but contains significant structure at intermediate masses. This structure is likely to be due to crowding of sources near or below the spatial resolution of the observations. There is little evidence of any difference between the index of the high-mass end of the clump mass function in the compressed region and in the unaffected cloud. The consequences of these results are discussed in terms of current models of triggered star formation.

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

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

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

We present high spatial resolution radio and near-infrared hydrogen recombination line observations of the southern massive star-forming region G333.6-0.2. The 3.4-cm continuum peak is found slightly offset from the infrared source. The H90alpha spectra show for the first time a double peak profile at some positions. The complex velocity structure may be accounted for by champagne outflows, which may also explain the offset between the radio and infrared sources. The 2.17-um Brackett gamma image and H90alpha map are combined to construct an extinction map which shows a trend probably set by the blister nature of the HII region. The total number of Lyman continuum photons in the central 50-arcsec is estimated to be equivalent to that emitted by up to 19 O7V stars.

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

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

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

The ROSAT-ESO Flux Limited X-ray (REFLEX) galaxy cluster survey and the 2dF Galaxy Redshift Survey (2dFGRS) respectively comprise the largest, homogeneous X-ray selected cluster catalogue and completed galaxy redshift survey. In this work we combine these two outstanding datasets in order to study the effect of the large-scale cluster environment, as traced by X-ray luminosity, on the properties of the cluster member galaxies. We measure the LX-sigma relation from the correlated dataset and find it to be consistent with recent results found in the literature. Using a sample of 19 clusters with LX>=0.36*10^44 erg s^-1 in the (0.1-2.4 keV) band, and 49 clusters with lower X-ray luminosity, we find that the fraction of early spectral type (eta<=-1.4), passively-evolving galaxies is significantly higher in the high-LX sample within R200. We extend the investigation to include composite bJ cluster luminosity functions, and find that the characteristic magnitude of the Schechter-function fit to the early-type luminosity function is fainter for the high-LX sample compared to the low-LX sample (Delta M*=0.58+/-0.14). This seems to be driven by a deficit of such galaxies with M_bJ ~ -21. In contrast, we find no significant differences between the luminosity functions of star-forming, late-type galaxies. We believe these results are consistent with a scenario in which the high-LX clusters are more dynamically evolved systems than the low-LX clusters.

We present a power spectrum analysis of the final 2dF Galaxy Redshift Survey, employing a direct Fourier method. The sample used comprises 221,414 galaxies with measured redshifts. We investigate in detail the modelling of the sample selection. A new angular mask is derived, based on revisions to the photometric calibration. The redshift selection function is determined by dividing the survey according to rest-frame colour, and deducing a self-consistent treatment of k-corrections and evolution for each population. The covariance matrix for the power-spectrum estimates is determined using two different approaches to the construction of mock surveys which are used to demonstrate that the input cosmological model can be correctly recovered. We are confident that the 2dFGRS power spectrum can be used to infer the matter content of the universe. On large scales, our estimated power spectrum shows evidence for the `baryon oscillations' that are predicted in CDM models. Fitting to a CDM model, assuming a primordial $n_{s}=1$ spectrum, $h=0.72$ and negligible neutrino mass, the preferred parameters are $\Omega_{M} h = 0.168 \pm 0.016$ and a baryon fraction $\Omega_{b} /\Omega_{M} = 0.185\pm0.046$ (1$\sigma$ errors). The value of $\Omega_{M} h$ is $1\sigma$ lower than the $0.20 \pm 0.03$ in our 2001 analysis of the partially complete 2dFGRS. This shift is largely due to the signal from the newly-sampled regions of space, rather than the refinements in the treatment of observational selection. This analysis therefore implies a density significantly below the standard $\Omega_{M} =0.3$: in combination with CMB data from WMAP, we infer $\Omega_{M} =0.231\pm 0.021$. (Abridged.)

We present a spectroscopic study of the shape of the dust-extinction law between 1.0 and 2.2um towards a set of nine ultracompact HII regions with Av > 15 mag. We find some evidence that the reddening curve may tend to flatten at higher extinctions, but just over half of the sample has extinction consistent with or close to the average for the interstellar medium. There is no evidence of extinction curves significantly steeper than the standard law, even where water ice is present. Comparing the results to the predictions of a simple extinction model, we suggest that a standard extinction law implies a robust upper limit to the grain-size distribution at around 0.1 - 0.3um. Flatter curves are most likely due to changes in this upper limit, although the effects of flattening due to unresolved clumpy extinction cannot be ruled out.

We use the 2dF Galaxy Redshift Survey to measure the dependence of the bJ-band galaxy luminosity function on large-scale environment, defined by density contrast in spheres of radius 8h-1Mpc, and on spectral type, determined from principal component analysis. We find that the galaxy populations at both extremes of density differ significantly from that at the mean density. The population in voids is dominated by late types and shows, relative to the mean, a deficit of galaxies that becomes increasingly pronounced at magnitudes brighter than M_bJ-5log10h <-18.5. In contrast, cluster regions have a relative excess of very bright early-type galaxies with M_bJ-5log10h < -21. Differences in the mid to faint-end population between environments are significant: at M_bJ-5log10h=-18 early and late-type cluster galaxies show comparable abundances, whereas in voids the late types dominate by almost an order of magnitude. We find that the luminosity functions measured in all density environments, from voids to clusters, can be approximated by Schechter functions with parameters that vary smoothly with local density, but in a fashion which differs strikingly for early and late-type galaxies. These observed variations, combined with our finding that the faint-end slope of the overall luminosity function depends at most weakly on density environment, may prove to be a significant challenge for models of galaxy formation.

It is well known that the clustering of galaxies depends on galaxy type.Such relative bias complicates the inference of cosmological parameters from galaxy redshift surveys, and is a challenge to theories of galaxy formation and evolution. In this paper we perform a joint counts-in-cells analysis on galaxies in the 2dF Galaxy Redshift Survey, classified by both colour and spectral type, eta, as early or late type galaxies. We fit three different models of relative bias to the joint probability distribution of the cell counts, assuming Poisson sampling of the galaxy density field. We investigate the nonlinearity and stochasticity of the relative bias, with cubical cells of side 10Mpc \leq L \leq 45Mpc (h=0.7). Exact linear bias is ruled out with high significance on all scales. Power law bias gives a better fit, but likelihood ratios prefer a bivariate lognormal distribution, with a non-zero `stochasticity' - i.e. scatter that may result from physical effects on galaxy formation other than those from the local density field. Using this model, we measure a correlation coefficient in log-density space (r_LN) of 0.958 for cells of length L=10Mpc, increasing to 0.970 by L=45Mpc. This corresponds to a stochasticity sigma_b/bhat of 0.44\pm0.02 and 0.27\pm0.05 respectively. For smaller cells, the Poisson sampled lognormal distribution presents an increasingly poor fit to the data, especially with regard to the fraction of completely empty cells. We compare these trends with the predictions of semianalytic galaxy formation models: these match the data well in terms of overall level of stochasticity, variation with scale, and fraction of empty cells.

We select a sample of low-redshift (z ~ 0.1) E+A galaxies from the 2dF Galaxy Redshift Survey (2dFGRS). The spectra of these objects are defined by strong hydrogen Balmer absorption lines (H-delta, H-gamma, H-beta) combined with a lack of [OII] 3727A emission, together implying a recently-truncated burst of star formation. The E+A spectrum is thus a signpost to galaxies in the process of evolution. We quantify the local environments, clustering properties and luminosity function of the E+A galaxies. We find that the environments are consistent with the ensemble of 2dFGRS galaxies: low-redshift E+A systems are located predominantly in the field, existing as isolated objects or in poor groups. However, the luminosity distribution of galaxies selected using three Balmer absorption lines H-delta-gamma-beta appears more typical of ellipticals. Indeed, morphologically these galaxies are preferentially spheroidal (E/S0) systems. In a small but significant number we find evidence for recent major mergers, such as tidal tails. We infer that major mergers are one important formation mechanism for E+A galaxies, as suggested by previous studies. At low redshift the merger probability is high in the field and low in clusters, thus these recently-formed spheroidal systems do not follow the usual morphology-density relation for ellipticals. Regarding the selection of E+A galaxies: we find that basing the Balmer-line criterion solely on H-delta absorption leads to a significant sub-population of disk systems with detectable H-alpha emission. In these objects the [OII] emission is presumably either obscured by dust or present with a low signal-to-noise ratio, whilst the (H-gamma, H-beta) absorption features are subject to emission-filling.

We measure moments of the galaxy count probability distribution function in the two-degree field galaxy redshift survey (2dFGRS). The survey is divided into volume limited subsamples in order to examine the dependence of the higher order clustering on galaxy luminosity. We demonstrate the hierarchical scaling of the averaged p-point galaxy correlation functions, xibar_p, up to p=6. The hierarchical amplitudes, S_p = xibar_p/xibar_2^{p-1}, are approximately independent of the cell radius used to smooth the galaxy distribution on small to medium scales. On larger scales we find the higher order moments can be strongly affected by the presence of rare, massive superstructures in the galaxy distribution. The skewness S_3 has a weak dependence on luminosity, approximated by a linear dependence on log luminosity. We discuss the implications of our results for simple models of linear and non-linear bias that relate the galaxy distribution to the underlying mass.

We measure the redshift space reduced void probability function (VPF) for 2dFGRS volume limited galaxy samples covering the absolute magnitude range M_bJ-5logh=-18 to -22. Theoretically, the VPF connects the distribution of voids to the moments of galaxy clustering of all orders, and can be used to discriminate clustering models in the weakly non-linear regime. The reduced VPF measured from the 2dFGRS is in excellent agreement with the paradigm of hierarchical scaling of the galaxy clustering moments. The accuracy of our measurement is such that we can rule out, at a very high significance, popular models for galaxy clustering, including the lognormal distribution. We demonstrate that the negative binomial model gives a very good approximation to the 2dFGRS data over a wide range of scales, out to at least 20h-1Mpc. Conversely, the reduced VPF for dark matter in a LambdaCDM universe does appear to be lognormal on small scales but deviates significantly beyond \approx 4h-1Mpc. We find little dependence of the 2dFGRS reduced VPF on galaxy luminosity. Our results hold independently in both the north and south Galactic pole survey regions.

We measure the clustering of galaxy groups in the 2dFGRS Percolation-Inferred Galaxy Group (2PIGG) catalogue. The 2PIGG sample has 29,000 groups with at least two members. The clustering amplitude of the full 2PIGG catalogue is weaker than that of 2dFGRS galaxies, in agreement with theoretical predictions. We have subdivided the 2PIGG catalogue into samples that span a factor of 25 in median total luminosity. Our correlation function measurements span an unprecedented range of clustering strengths, connecting the regimes probed by groups fainter than L* galaxies and rich clusters. There is a steady increase in clustering strength with group luminosity; the most luminous groups are ten times more strongly clustered than the full 2PIGG catalogue. We demonstrate that the 2PIGG results are in very good agreement with the clustering of groups expected in the LCDM model.

We present the result of a decomposition of the 2dFGRS galaxy overdensity field into an orthonormal basis of spherical harmonics and spherical Bessel functions. Galaxies are expected to directly follow the bulk motion of the density field on large scales, so the absolute amplitude of the observed large-scale redshift-space distortions caused by this motion is expected to be independent of galaxy properties. By splitting the overdensity field into radial and angular components, we linearly model the observed distortion and obtain the cosmological constraint Omega_m^{0.6} sigma_8=0.46+/-0.06. The amplitude of the linear redshift-space distortions relative to the galaxy overdensity field is dependent on galaxy properties and, for L_* galaxies at redshift z=0, we measure beta(L_*,0)=0.58+/-0.08, and the amplitude of the overdensity fluctuations b(L_*,0) sigma_8=0.79+/-0.03, marginalising over the power spectrum shape parameters. Assuming a fixed power spectrum shape consistent with the full Fourier analysis produces very similar parameter constraints.

We use the two-degree field Galaxy Redshift Survey (2dFGRS) to test the hierarchical scaling hypothesis: namely, that the $p$-point galaxy correlation functions can be written in terms of the two point correlation function or variance. This scaling is expected if an initially Gaussian distribution of density fluctuations evolves under the action of gravitational instability. We measure the volume averaged $p$-point correlation functions using a counts in cells technique applied to a volume limited sample of 44,931 $L_*$ galaxies. We demonstrate that $L_{*}$ galaxies display hierarchical clustering up to order $p=6$ in redshift space. The variance measured for $L_{*}$ galaxies is in excellent agreement with the predictions from a $\Lambda$-cold dark matter N-body simulation. This applies to all cell radii considered, $0.3<(R/h^{-1}{\rm Mpc})<30$. However, the higher order correlation functions of $L_*$ galaxies have a significantly smaller amplitude than is predicted for the dark matter for $R<10h^{-1}$Mpc. This disagreement implies that a non-linear bias exists between the dark matter and $L_*$ galaxies on these scales. We also show that the presence of two rare, massive superclusters in the 2dFGRS has an impact on the higher-order clustering moments measured on large scales.

Complementary one-, two-, and three-dimensional tests for detecting the presence of substructure in clusters of galaxies are applied to recently obtained data from the 2dF Galaxy Redshift Survey. The sample of 25 clusters used in this study includes 16 clusters not previously investigated for substructure. Substructure is detected at or greater than the 99% CL level in at least one test for 21 of the 25 clusters studied here. From the results, it appears that low richness clusters commonly contain subclusters participating in mergers. About half of the clusters have two or more components within 0.5 h^{-1} Mpc of the cluster centroid, and at least three clusters (Abell 1139, Abell 1663, and Abell S333) exhibit velocity-position characteristics consistent with the presence of possible cluster rotation, shear, or infall dynamics. The geometry of certain features is consistent with influence by the host supercluster environments. In general, our results support the hypothesis that low richness clusters relax to structureless equilibrium states on very long dynamical time scales (if at all).

We present an analysis of the relative bias between early- and late-type galaxies in the Two-degree Field Galaxy Redshift Survey (2dFGRS). Our analysis examines the joint counts in cells between early- and late-type galaxies, using approximately cubical cells with sides ranging from 7h^{-1}Mpc to 42h^{-1}Mpc. We measure the variance of the counts in cells using the method of Efstathiou et al. (1990), which we find requires a correction for a finite volume effect. We fit lognormal models to the one-point density distribution and develop methods of dealing with biases in the recovered variances resulting from this technique. We directly fit deterministic models for the joint density distribution function, f(delta_E,delta_L), to the joint counts in cells using a maximum likelihood technique. Our results are consistent with a scale invariant relative bias factor on all scales studied. Linear bias is ruled out on scales less than l=28h^{-1}Mpc. A power-law bias model is a significantly better fit to the data on all but the largest scales studied; the relative goodness of fit of this model as compared to that of the linear bias model suggests that any nonlinearity is negligible for l>~40h^{-1}Mpc, consistent with the expectation from theory that the bias should become linear on large scales. (abridged)

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

We derive the fraction of blue galaxies in a sample of clusters at z < 0.11 and the general field at the same redshift. The value of the blue fraction is observed to depend on the luminosity limit adopted, cluster-centric radius and, more generally, local galaxy density, but it does not depend on cluster properties. Changes in the blue fraction are due to variations in the relative proportions of red and blue galaxies but the star formation rate for these two galaxy groups remains unchanged. Our results are most consistent with a model where the star formation rate declines rapidly and the blue galaxies tend to be dwarfs and do not favour mechanisms where the Butcher-Oemler effect is caused by processes specific to the cluster environment.

The construction of a catalogue of galaxy groups from the 2-degree Field Galaxy Redshift Survey (2dFGRS) is described. Groups are identified by means of a friends-of-friends percolation algorithm which has been thoroughly tested on mock versions of the 2dFGRS generated from cosmological N-body simulations. The tests suggest that the algorithm groups all galaxies that it should be grouping, with an additional 40% of interlopers. About 55% of the ~190000 galaxies considered are placed into groups containing at least two members of which ~29000 are found. Of these, ~7000 contain at least four galaxies, and these groups have a median redshift of 0.11 and a median velocity dispersion of 260km/s. This 2dFGRS Percolation-Inferred Galaxy Group (2PIGG) catalogue represents the largest available homogeneous sample of galaxy groups. It is publicly available on the WWW.

The 2dFGRS Percolation-Inferred Galaxy Group (2PIGG) catalogue of ~29000 objects is used to study the luminous content of galaxy systems of various sizes. Mock galaxy catalogues constructed from cosmological simulations are used to gauge the accuracy with which intrinsic group properties can be recovered. A Schechter function is found to be a reasonable fit to the galaxy luminosity functions in groups of different mass in the real data. The characteristic luminosity L* is larger for more massive groups. However, the mock data show that the shape of the recovered luminosity function is expected to differ from the true shape, and this must be allowed for when interpreting the data. The variation of halo mass-to-light ratio with group size is studied in both these wavebands. A robust trend of increasing M/L with increasing group luminosity is found in the 2PIGG data. From groups with L_bj=10^{10}Lsol to those 100 times more luminous, the typical bj-band M/L increases by a factor of 5, whereas the rf-band M/L grows by a factor of 3.5. These trends agree well with the simulations, which also predict a minimum M/L on a scale corresponding to the Local Group. Our data indicate that if such a minimum exists, then it must occur at L<~10^{10}Lsol, below the range accurately probed by the 2PIGG catalogue. According to the mock data, the bj M/Ls of the largest groups are expected to be approximately 1.1 times the global value. Assuming that this correction applies to the real data yields an estimate of Omega_m=0.26+/-0.03 (statistical).

The clustering properties of local, S_{1.4 GHz} > 1 mJy, radio sources are investigated for a sample of 820 objects drawn from the joint use of the FIRST and 2dF Galaxy Redshift surveys. To this aim, we present 271 new bj < 19.45 spectroscopic counterparts of FIRST radio sources to be added to those already introduced in Magliocchetti et al. (2002). The two-point correlation function for the local radio population is found to be entirely consistent with estimates obtained for the whole sample of 2dFGRS galaxies. We estimate the parameters of the real-space correlation function xi(r)=(r/r_0)^{-\gamma}, r_0=6.7^{+0.9}_{-1.1} Mpc and \gamma=1.6\pm 0.1, where h=0.7 is assumed. Different results are instead obtained if we only consider sources that present signatures of AGN activity in their spectra. These objects are shown to be very strongly correlated, with r_0=10.9^{+1.0}_{-1.2} Mpc and \gamma=2\pm 0.1, a steeper slope than has been claimed in other recent works. No difference is found in the clustering properties of radio-AGNs of different radio luminosity. These results show that AGN-fuelled sources reside in dark matter halos more massive than \sim 10^{13.4} M_{\sun}},higher the corresponding figure for radio-quiet QSOs. This value can be converted into a minimum black hole mass associated with radio-loud, AGN-fuelled objects of M_{BH}^{min}\sim 10^9 M_{\sun}. The above results then suggest -at least for relatively faint radio objects -the existence of a threshold black hole mass associated with the onset of significant radio activity such as that of radio-loud AGNs; however, once the activity is triggered, there appears to be no evidence for a connection between black hole mass and level of radio output. (abridged)

We analyse the observed correlation between galaxy environment and H-alpha emission line strength, using volume-limited samples and group catalogues of 24968 galaxies drawn from the 2dF Galaxy Redshift Survey (Mb<-19.5) and the Sloan Digital Sky Survey (Mr<-20.6). We characterise the environment by 1) Sigma_5, the surface number density of galaxies determined by the projected distance to the 5th nearest neighbour; and 2) rho1.1 and rho5.5, three-dimensional density estimates obtained by convolving the galaxy distribution with Gaussian kernels of dispersion 1.1 Mpc and 5.5 Mpc, respectively. We find that star-forming and quiescent galaxies form two distinct populations, as characterised by their H-alpha equivalent width, EW(Ha). The relative numbers of star-forming and quiescent galaxies varies strongly and continuously with local density. However, the distribution of EW(Ha) amongst the star-forming population is independent of environment. The fraction of star-forming galaxies shows strong sensitivity to the density on large scales, rho5.5, which is likely independent of the trend with local density, rho1.1. We use two differently-selected group catalogues to demonstrate that the correlation with galaxy density is approximately independent of group velocity dispersion, for sigma=200-1000 km/s. Even in the lowest density environments, no more than ~70 per cent of galaxies show significant H-alpha emission. Based on these results, we conclude that the present-day correlation between star formation rate and environment is a result of short-timescale mechanisms that take place preferentially at high redshift, such as starbursts induced by galaxy-galaxy interactions.

We reconstruct the underlying density field of the 2 degree Field Galaxy Redshift Survey (2dFGRS) for the redshift range 0.035<z<0.200 using the Wiener Filtering method. The Wiener Filter suppresses shot noise and accounts for selection and incompleteness effects. The method relies on prior knowledge of the 2dF power spectrum of fluctuations and the combination of matter density and bias parameters however the results are only slightly affected by changes to these parameters. We present maps of the density field in two different resolutions: 5 Mpc/h and 10 Mpc/h. We identify all major superclusters and voids in the survey. In particular, we find two large superclusters and two large local voids. A version of this paper with full set of colour maps can be found at http://www.ast.cam.ac.uk/~pirin.

We present new spectropolarimetry of a sample of nearby Compton-thin Seyfert 2 galaxies (ie those with N_H<10^23 cm^-2). We show that the detection rate of scattered broad Halpha in this sample is considerably higher than in Seyfert 2 galaxies as a whole. Our results also show that in this low obscuration set it is possible to find scattered broad Halpha even when the global properties of the galaxy are largely dominated by the host galaxy and not the active galactic nucleus. These results argue against the existence of a population of `pure' Seyfert 2 galaxies.

We present a detailed analysis of the two-point correlation function, from the 2dF Galaxy Redshift Survey (2dFGRS). We estimate the redshift-space correlation function, xi(s), from which we measure the redshift-space clustering length, s_0=6.82+/-0.28 Mpc/h. We also estimate the projected correlation function, Xi(sigma), and the real-space correlation function, xi(r), which can be fit by a power-law, with r_0=5.05+/-0.26Mpc/h, gamma_r=1.67+/-0.03. For r>20Mpc/h, xi drops below a power-law as is expected in the popular LCDM model. The ratio of amplitudes of the real and redshift-space correlation functions on scales of 8-30Mpc/h gives an estimate of the redshift-space distortion parameter beta. The quadrupole moment of xi on scales 30-40Mpc/h provides another estimate of beta. We also estimate the distribution function of pairwise peculiar velocities, f(v), including rigorously the effect of infall velocities, and find that it is well fit by an exponential. The accuracy of our xi measurement is sufficient to constrain a model, which simultaneously fits the shape and amplitude of xi(r) and the two redshift-space distortion effects parameterized by beta and velocity dispersion, a. We find beta=0.49+/-0.09 and a=506+/-52km/s, though the best fit values are strongly correlated. We measure the variation of the peculiar velocity dispersion with projected separation, a(sigma), and find that the shape is consistent with models and simulations. Using the constraints on bias from recent estimates, and taking account of redshift evolution, we conclude that beta(L=L*,z=0)=0.47+/-0.08, and that the present day matter density of the Universe is 0.3, consistent with other 2dFGRS estimates and independent analyses.

The 2dF Galaxy Redshift Survey (2dFGRS) has obtained spectra for 245591 sources, mainly galaxies, brighter than a nominal extinction-corrected magnitude limit of b_J=19.45. Reliable redshifts were measured for 221414 galaxies. The galaxies are selected from the extended APM Galaxy Survey and cover an area of approximately 1500 square degrees in three regions: an NGP strip, an SGP strip and random fields scattered around the SGP strip. This paper describes the 2dFGRS final data release of 30 June 2003 and complements Colless et al. (2001), which described the survey and the initial 100k data release. The 2dFGRS database and full documentation are available on the WWW at http://www.mso.anu.edu.au/2dFGRS/

We have calculated the two-point correlation functions in redshift space, xi(sigma,pi), for galaxies of different spectral types in the 2dF Galaxy Redshift Survey. Using these correlation functions we are able to estimate values of the linear redshift-space distortion parameter, beta = Omega_m^0.6/b, the pairwise velocity dispersion, a, and the real-space correlation function, xi(r), for galaxies with both relatively low star-formation rates (for which the present rate of star formation is less than 10% of its past averaged value) and galaxies with higher current star-formation activity. At small separations, the real-space clustering of passive galaxies is very much stronger than that of the more actively star-forming galaxies; the correlation-function slopes are respectively 1.93 and 1.50, and the relative bias between the two classes is a declining function of radius. On scales larger than 10 h^-1 Mpc there is evidence that the relative bias tends to a constant, b(passive)/b(active) ~ 1. This result is consistent with the similar degrees of redshift-space distortions seen in the correlation functions of the two classes -- the contours of xi(sigma,pi) require beta(active)=0.49+/-0.13, and beta(passive)=0.48+/-0.14. The pairwise velocity dispersion is highly correlated with beta. However, despite this a significant difference is seen between the two classes. Over the range 8-20 h^-1 Mpc, the pairwise velocity dispersion has mean values 416+/-76 km/s and 612+/-92 km/s for the active and passive galaxy samples respectively. This is consistent with the expectation from morphological segregation, in which passively evolving galaxies preferentially inhabit the cores of high-mass virialised regions.

We have determined the composite luminosity function (LF) for galaxies in 60 clusters from the 2dF Galaxy Redshift Survey. The LF spans the range $-22.5<M_{b_{\rm J}}<-15$, and is well-fitted by a Schechter function with ${M_{b_{\rm J}}}^{*}=-20.07\pm0.07$ and $\alpha=-1.28\pm0.03$ ($H_0$=100 km s$^{-1}$ Mpc$^{-1}$, $\Omega_M$=0.3, $\Omega_\Lambda$=0.7). It differs significantly from the field LF of \cite{mad02}, having a characteristic magnitude that is approximately 0.3 mag brighter and a faint-end slope that is approximately 0.1 steeper. There is no evidence for variations in the LF across a wide range of cluster properties. However the LF of early-type galaxies in clusters is both brighter and steeper than its field counterpart. The differences between the field and cluster LFs for the various spectral types can be qualitatively explained by the suppression of star formation in the dense cluster environment, together with mergers to produce the brightest early-type galaxies.

We present observations and models of the behaviour of the HI and HeI lines between 1.6 and 2.2um in a small sample of compact HII regions. As in our previous papers on planetary nebulae, we find that the `pure' 1.7007um 4^3D-3^3P and 2.16475um 7^(3,1)G-4^(3,1)F HeI recombination lines behave approximately as expected as the effective temperature of the central exciting star(s) increases. However, the 2.058um 2^1P-2^1S HeI line does not behave as the model predicts, or as seen in planetary nebulae. Both models and planetary nebulae showed a decrease in the HeI 2^1P-2^1S/HI Br gamma ratio above an effective temperature of 40000K. The compact HII regions do not show any such decrease. The problem with this line ratio is probably due to the fact that the photoionisation model does not account correctly for the high densities seen in these HII regions, and that we are therefore seeing more collisional excitation of the 2^1P level than the model predicts. It may also reflect some deeper problem in the assumed model stellar atmospheres. In any event, although the normal HeI recombination lines can be used to place constraints on the temperature of the hottest star present, the HeI 2^1P-2^1S/HI Br gamma ratio should not be used for this purpose in either Galactic HII regions or in starburst galaxies, and conclusions from previous work using this ratio should be regarded with extreme caution. We also show that the combination of the near infrared `pure' recombination line ratios with mid-infrared forbidden line data provides a good discriminant of the form of the far ultraviolet spectral energy distribution of the exciting star(s). From this we conclude that CoStar models are a poor match to the available data for our sources, though the more recent WM-basic models are a better fit.

We constrain flat cosmological models with a joint likelihood analysis of a new compilation of data from the cosmic microwave background (CMB) and from the 2dF Galaxy Redshift Survey (2dFGRS). Fitting the CMB alone yields a known degeneracy between the Hubble constant h and the matter density Omega_m, which arises mainly from preserving the location of the peaks in the angular power spectrum. This `horizon-angle degeneracy' is considered in some detail and shown to follow a simple relation Omega_m h^{3.4} = constant. Adding the 2dFGRS power spectrum constrains Omega_m h and breaks the degeneracy. If tensor anisotropies are assumed to be negligible, we obtain values for the Hubble constant h=0.665 +/- 0.047, the matter density Omega_m=0.313 +/- 0.055, and the physical CDM and baryon densities Omega_c h^2 = 0.115 +/- 0.009, Omega_b h^2 = 0.022 +/- 0.002 (standard rms errors). Including a possible tensor component causes very little change to these figures; we set a upper limit to the tensor-to-scalar ratio of r<0.7 at 95% confidence. We then show how these data can be used to constrain the equation of state of the vacuum, and find w<-0.52 at 95% confidence. The preferred cosmological model is thus very well specified, and we discuss the precision with which future CMB data can be predicted, given the model assumptions. The 2dFGRS power-spectrum data and covariance matrix, and the CMB data compilation used here, are available from http://www.roe.ac.uk/~wjp/

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

We constrain f_nu = Omega_nu / Omega_m, the fractional contribution of neutrinos to the total mass density in the Universe, by comparing the power spectrum of fluctuations derived from the 2dF Galaxy Redshift Survey with power spectra for models with four components: baryons, cold dark matter, massive neutrinos and a cosmological constant. Adding constraints from independent cosmological probes we find f_nu < 0.13 (at 95% confidence) for a prior of 0.1< Omega_m <0.5, and assuming the scalar spectral index n=1. This translates to an upper limit on the total neutrino mass and m_nu,tot < 1.8 eV for "concordance" values of Omega_m and the Hubble constant. Very similar results are obtained with a prior on Omega_m from Type Ia supernovae surveys, and with marginalization over n.

The 2dF Galaxy Redshift Survey is the first to observe more than 100,000 redshifts. This allows precise measurements of many of the key statistics of galaxy clustering, in particular redshift-space distortions and the large-scale power spectrum. This paper presents the current 2dFGRS results in these areas. Redshift-space distortions are detected with a high degree of significance, confirming the detailed Kaiser distortion from large-scale infall velocities, and measuring the distortion parameter beta equiv Omega_m^{0.6}/b = 0.43 +- 0.07. The power spectrum is measured to < 10% accuracy for k > 0.02 h Mpc^{-1}, and is well fitted by a CDM model with Omega_m h = 0.20 +- 0.03 and a baryon fraction of 0.15 +- 0.07. A joint analysis with CMB data requires Omega_m = 0.29 +- 0.05, assuming scalar fluctuations, but no priors on other parameters. Two methods are used to determine the large-scale bias parameter: an internal bispectrum analysis yields b = 1.04 +- 0.11, in very good agreement with the b = 1.10 +- 0.08 obtained from a joint 2dFGRS+CMB analysis, again assuming scalar fluctuations. These figures refer to galaxies of approximate luminosity 2L^*; luminosity dependence of clustering is detected at high significance, and is well described by b/b^* = 0.85 + 0.15(L/L^*).

We compare the amplitudes of fluctuations probed by the 2dF Galaxy Redshift Survey and by the latest measurements of the Cosmic Microwave Background anisotropies. By combining the 2dFGRS and CMB data we find the linear-theory rms mass fluctuations in 8 Mpc/h spheres to be sigma_8 = 0.73 +-0.05 (after marginalization over the matter density parameter Omega_m and three other free parameters). This normalization is lower than the COBE normalization and previous estimates from cluster abundance, but it is in agreement with some revised cluster abundance determinations. We also estimate the scale-independent bias parameter of present-epoch L_s = 1.9L_* APM-selected galaxies to be b(L_s,z=0) = 1.10 +- 0.08 on comoving scales of 0.02 < k < 0.15 h/Mpc. If luminosity segregation operates on these scales, L_* galaxies would be almost un-biased, b(L_*,z=0) = 0.96. These results are derived by assuming a flat Lambda-CDM Universe, and by marginalizing over other free parameters and fixing the spectral index n=1 and the optical depth due to reionization tau=0. We also study the best fit pair (Omega_m,b), and the robustness of the results to varying n and tau. Various modelling corrections can each change the resulting b by 5-15 per cent. The results are compared with other independent measurements from the 2dFGRS itself, and from the SDSS, cluster abundance and cosmic shear.

We compute the bispectrum of the 2dF Galaxy Redshift Survey (2dFGRS) and use it to measure the bias parameter of the galaxies. This parameter quantifies the strength of clustering of the galaxies relative to the mass in the Universe. By analysing 80 million triangle configurations in the wavenumber range 0.1 < k < 0.5 h/Mpc (i.e. on scales roughly between 5 and 30 Mpc/h) we find that the linear bias parameter is consistent with unity: b_1=1.04 pm 0.11, and the quadratic (nonlinear) bias is consistent with zero: b_2=-0.054 pm 0.08. Thus, at least on large scales, optically-selected galaxies do indeed trace the underlying mass distribution. The bias parameter can be combined with the 2dFGRS measurement of the redshift distortion parameter beta = Omega_m^{0.6}/b_1, to yield Omega_m = 0.27 pm 0.06 for the matter density of the Universe, a result which is determined entirely from this survey, independently of other datasets. Our measurement of the matter density of the Universe should be interpreted as Omega_m at the effective redshift of the survey (z=0.17).

We use more than 110500 galaxies from the 2dF galaxy redshift survey (2dFGRS) to estimate the b_J-band galaxy luminosity function at redshift z=0, taking account of evolution, the distribution of magnitude measurement errors and small corrections for incompletenessin the galaxy catalogue. Throughout the interval -16.5>M- 5log h>-22, the luminosity function is accurately described by a Schechter function with M* -5log h =-19.66+/-0.07, alpha=-1.21+/-0.03 and phistar=(1.61+/-0.08) 10^{-2} h^3/Mpc^3, giving an integrated luminosity density of rho_L=(1.82+/-0.17) 10^8 h L_sol/Mpc^3 (assuming an Omega_0=0.3, Lambda_0=0.7 cosmology). The quoted errors have contributions from the accuracy of the photometric zeropoint, large scale structure in the galaxy distribution and, importantly, from the uncertainty in the appropriate evolutionary corrections. Our luminosity function is in excellent agreement with, but has much smaller statistical errors than an estimate from the Sloan Digital Sky Survey (SDSS) data when the SDSS data are accurately translated to the b_J-band and the luminosity functions are normalized in the same way. We use the luminosity function, along with maps describing the redshift completeness of the current 2dFGRS catalogue, and its weak dependence on apparent magnitude, to define a complete description of the 2dFGRS selection function. Details and tests of the calibration of the 2dFGRS photometric parent catalogue are also presented.

We have measured the equivalent width of the H-alpha emission line for 11006 galaxies brighter than M_b=-19 (LCDM) at 0.05<z<0.1 in the 2dF Galaxy Redshift Survey (2dF), in the fields of seventeen known galaxy clusters. The limited redshift range ensures that our results are insensitive to aperture bias, and to residuals from night sky emission lines. We use these measurements to trace mustar, the star formation rate normalized to Lstar, as a function of distance from the cluster centre, and local projected galaxy density. We find that the distribution of mustar steadily skews toward larger values with increasing distance from the cluster centre, converging to the field distribution at distances greater than ~3 times the virial radius. A correlation between star formation rate and local projected density is also found, which is independent of cluster velocity dispersion and disappears at projected densities below ~1 galaxy (brighter than M_b=-19) per Mpc^2. This characteristic scale corresponds approximately to the mean density at the cluster virial radius. The same correlation holds for galaxies more than two virial radii from the cluster centre. We conclude that environmental influences on galaxy properties are not restricted to cluster cores, but are effective in all groups where the density exceeds this critical value. The present day abundance of such systems, and the strong evolution of this abundance, makes it likely that hierarchical growth of structure plays a significant role in decreasing the global average star formation rate. Finally, the low star formation rates well beyond the virialised cluster rule out severe physical processes, such as ram pressure stripping of disk gas, as being completely responsible for the variations in galaxy properties with environment.

We calculate the optical bj luminosity function of the 2dF Galaxy Redshift Survey (2dFGRS) for different subsets defined by their spectral properties. These spectrally selected subsets are defined using a new parameter, eta, which is a linear combination of the first two projections derived from a Principal Component Analysis. This parameter eta identifies the average emission and absorption line strength in the galaxy rest-frame spectrum and hence is a useful indicator of the present star formation. We use a total of 75,000 galaxies in our calculations, chosen from a sample of high signal-to-noise ratio, low redshift galaxies observed before January 2001. We find that there is a systematic steepening of the faint end slope (alpha) as one moves from passive (alpha = -0.54) to active (alpha = -1.50) star-forming galaxies, and that there is also a corresponding faintening of the rest-frame characteristic magnitude M* - 5\log_10(h) (from -19.6 to -19.2). We also show that the Schechter function provides a poor fit to the quiescent (Type 1) LF for very faint galaxies (M - 5log_10(h) fainter than -16.0), perhaps suggesting the presence of a significant dwarf population. The luminosity functions presented here give a precise confirmation of the trends seen previously in a much smaller preliminary 2dFGRS sample, and in other surveys. We also present a new procedure for determining self-consistent K-corrections and investigate possible fibre-aperture biases.

We present the bivariate brightness distribution (BBD) for the 2dF Galaxy Redshift Survey (2dFGRS) based on a preliminary subsample of 45,000 galaxies. The BBD is an extension of the galaxy luminosity function incorporating surface brightness information. It allows the measurement of the local luminosity density, j_B, and the galaxy luminosity and surface brightness distributions while accounting for surface brightness selection biases. The recovered 2dFGRS BBD shows a strong surface brightness-luminosity relation (M_B~2.4\mu_e). The luminosity-density is dominated by normal galaxies and the luminosity-density peak lies away from the selection boundaries implying that the 2dFGRS is complete and that luminous low surface brightness galaxies are rare. The final value we derive for the local luminosity-density, inclusive of surface brightness corrections, is: j_B=2.49+/-0.20x10^8 h L_solar Mpc^-3. Representative Schechter function parameters are: M*=-19.75+/-0.05, phi*=2.02+/-0.02x10^-2 and alpha=-1.09+/-0.03. Extending the conventional methodology to incorporate surface brightness selection effects has resulted in an increase in the luminosity-density of 37%.

We use redshift determinations and spectral analysis of galaxies in the 2dF Galaxy Redshift Survey to study the properties of local radio sources with $S\ge 1$ mJy. 557 objects drawn from the FIRST survey, corresponding to 2.3 per cent of the total radio sample, are found in the 2dFGRS cataloge within the area $9^h 48^m \simlt {\rm RA}({\rm 2000}) \simlt 14^h 32^m$ and $-2.77^\circ \simlt {\rm dec}({\rm 2000}) \simlt 2.25^\circ$, down to a magnitude limit $b_J= 19.45$. The excellent quality of 2dF spectra allows us to divide these sources into classes, according to their optical spectra. Radio galaxies (i.e. absorption systems, mainly belonging to the class of FR I sources, whose spectra may or may not show emission lines due to AGN activity) make up 63 per cent of the spectroscopic sample; they are characterized by relatively high radio-to-optical ratios, red colours and luminosities $10^{21}\simlt P_{1.4{\rm GHz}} W Hz^${-1}$ sr$^{-1}$ \simlt 10^{24}$. Actively star-forming galaxies instead contribute with 30 per cent to the total sample. These objects -- principally found at low redshifts ($z \simlt 0.1$) -- in general show low radio-to-optical ratios, blue colours and low radio luminosities. We also found 18 Seyfert 2 galaxies and 4 Seyfert 1's. Analysis of the local radio luminosity function shows that radio galaxies are well described by models that assume pure luminosity evolution, at least down to radio powers $P_{1.4{\rm GHz}}\simlt 10^{20.5} W Hz$^{-1}$ sr$^{-1}$ Late-type galaxies, whose relative contribution to the RLF is found to be lower than was predicted by previous work, present a luminosity function which is compatible with the IRAS galaxy LF, and therefore plausibly constitute the radio counterpart of the dusty spirals and starbursts that dominate the counts at 60 $\mu$m.

We have carried out a study of known clusters within the 2dF Galaxy Redshift Survey (2dFGRS) observed areas and have identified 431 Abell, 173 APM and 343 EDCC clusters. Precise redshifts, velocity dispersions and new centroids have been measured for the majority of these objects, and this information has been used to study the completeness of these catalogues, the level of contamination from foreground and background structures along the cluster's line of sight, the space density of the clusters as a function of redshift, and their velocity dispersion distributions. We find that the Abell and EDCC catalogues are contaminated at the level of about 10%, whereas the APM catalogue suffers only 5% contamination. If we use the original catalog centroids, the level of contamination rises to approximately 15% for the Abell and EDCC catalogues, showing that the presence of foreground and background groups may alter the richness of clusters in these catalogues. There is a deficiency of clusters at $z \sim 0.05$ that may correspond to a large underdensity in the Southern hemisphere. From the cumulative distribution of velocity dispersions for these clusters, we derive an upper limit to the space density of $\sigma > 1000 \kms$ clusters of $3.6 \times 10^{-6} \hdens$. This result is used to constrain models for structure formation; our data favour low-density cosmologies, subject to the usual assumptions concerning the shape and normalization of the power spectrum.

We have cross-matched the 1.4 GHz NRAO VLA Sky Survey (NVSS) with the first 210 fields observed in the 2dF Galaxy Redshift Survey (2dFGRS), covering an effective area of 325 square degrees (about 20% of the final 2dFGRS area). This yields a set of optical spectra of 912 candidate NVSS counterparts, of which we identify 757 as genuine radio IDs - the largest and most homogeneous set of radio-source spectra ever obtained. The 2dFGRS radio sources span the redshift range z=0.005 to 0.438, and are a mixture of active galaxies (60%) and star-forming galaxies (40%). About 25% of the 2dFGRS radio sources are spatially resolved by NVSS, and the sample includes three giant radio galaxies with projected linear size greater than 1 Mpc. The high quality of the 2dF spectra means we can usually distinguish unambiguously between AGN and star-forming galaxies. We have made a new determination of the local radio luminosity function at 1.4 GHz for both active and star-forming galaxies, and derive a local star-formation density of 0.022+/-0.004 solar masses per year per cubic Mpc. (Ho=50 km/s/Mpc).

We investigate the dependence of galaxy clustering on luminosity and spectral type using the 2dF Galaxy Redshift Survey (2dFGRS). Spectral types are assigned using the principal component analysis of Madgwick et al. We divide the sample into two broad spectral classes: galaxies with strong emission lines (`late-types'), and more quiescent galaxies (`early-types'). We measure the clustering in real space, free from any distortion of the clustering pattern due to peculiar velocities, for a series of volume-limited samples. The projected correlation functions of both spectral types are well described by a power law for transverse separations in the range 2 < (sigma/Mpc/h) < 15, with a marginally steeper slope for early-types than late-types. Both early and late types have approximately the same dependence of clustering strength on luminosity, with the clustering amplitude increasing by a factor of ~2.5 between L* and 4 L*. At all luminosities, however, the correlation function amplitude for the early-types is ~50% higher than that of the late-types. These results support the view that luminosity, and not type, is the dominant factor in determining how the clustering strength of the whole galaxy population varies with luminosity.

We present the first results on the history of star formation in the Universe based on the `cosmic spectrum', in particular, the volume-averaged, luminosity-weighted, stellar absorption line spectrum of present day galaxies from the 2dFGRS. This method is novel in that unlike previous studies it is not an estimator based on total luminosity density. The cosmic spectrum is fitted with models of population synthesis, tracing the history of star formation prior to the epoch of the observed galaxies, using a method we have developed which decouples continuum and spectral-line variations and is robust against spectrophotometric uncertainties. The cosmic spectrum can only be fitted with models incorporating chemical evolution and indicates there was a peak of star-formation rate in the past of at least three times the current value and that the increase back to z=1, assuming it scales as (1+z)^beta, has a strong upper limit of beta<5. We find in the general case there is some model degeneracy between star formation at low and high redshift. However, if we incorporate previous work on star formation at z<1 we can put strong upper limits on the SFR at z>1: e.g., if beta>2 then the SFR for 1<z<5 scales as (1+z)^alpha with alpha<2. This is equivalent to stating that no more than 80% of stars in the Universe formed at z>1. Our results are consistent with the best-fit results from compilations of cosmic SFR estimates based on UV luminosity density, which give 1.8<beta<2.9 and -1.0<alpha<0.7, and are also consistent with estimates of Omega_stars based on the K-band luminosity density.

We combine the 2MASS extended source catalogue and the 2dFGRS to produce an IR selected galaxy catalogue with 17,173 measured redshifts. We use this extensive dataset to estimate the J and K-band galaxy luminosity functions. The LFs are fairly well fit by Schechter functions with J: M*-5log h= -22.36+/-0.02, alpha= -0.93+/-0.04, Phi=0.0104+/-0.0016 h^3/Mpc^3 and K: M*-5log h= -23.44+/-0.03, alpha=-0.96+/-0.05, Phi=0.0108+/-0.0016 h^3/Mpc^3 (2MASS Kron magnitudes). These parameters assume a cosmological model with Omega=0.3 and Lambda=0.7. With datasets of this size, systematic rather than random errors are the dominant source of uncertainty in the determination of the LF. We carry out a careful investigation of possible systematic effects in our data. The surface brightness distribution of the sample shows no evidence that significant numbers of low surface brightness or compact galaxies are missed by the survey. We estimate the present-day distributions of B-K and J-K colours as a function of absolute magnitude and use models of the galaxy stellar populations, constrained by the observed optical and infrared colours, to infer the galaxy stellar mass function. Integrated over all galaxy masses, this yields a total mass fraction in stars (in units of the critical mass density) of Omega_*.h= (1.6+/-0.24)/10^3 for a Kennicutt IMF and Omega_*.h= (2.9+/-0.43)/10^3 for a Salpeter IMF. These values agree with those inferred from observational estimates of the star formation history of the universe provided that dust extinction corrections are modest.

The 2dF Galaxy Redshift Survey (2dFGRS) is designed to measure redshifts for approximately 250000 galaxies. This paper describes the survey design, the spectroscopic observations, the redshift measurements and the survey database. The 2dFGRS uses the 2dF multi-fibre spectrograph on the Anglo-Australian Telescope, which is capable of observing 400 objects simultaneously over a 2-degree diameter field. The source catalogue for the survey is a revised and extended version of the APM galaxy catalogue, and the targets are galaxies with extinction-corrected magnitudes brighter than b_J=19.45. The main survey regions are two declination strips, one in the southern Galactic hemisphere spanning 80deg x 15deg around the SGP, and the other in the northern Galactic hemisphere spanning 75deg x 10deg along the celestial equator; in addition, there are 99 fields spread over the southern Galactic cap. The survey covers 2000 sq.deg and has a median depth of z=0.11. Adaptive tiling is used to give a highly uniform sampling rate of 93% over the whole survey region. Redshifts are measured from spectra covering 3600A-8000A at a two-pixel resolution of 9.0A and a median S/N of 13 per pixel. All redshift identifications are visually checked and assigned a quality parameter Q in the range 1-5; Q>=3 redshifts are 98.4% reliable and have an rms uncertainty of 85 km/s. The overall redshift completeness for Q>=3 redshifts is 91.8%, but this varies with magnitude from 99% for the brightest galaxies to 90% for objects at the survey limit. The 2dFGRS database is available on the WWW at http://www.mso.anu.edu.au/2dFGRS

We investigate the dependence of the strength of galaxy clustering on intrinsic luminosity using the Anglo-Australian two degree field galaxy redshift survey (2dFGRS). The 2dFGRS is over an order of magnitude larger than previous redshift surveys used to address this issue. We measure the projected two-point correlation function of galaxies in a series of volume-limited samples. The projected correlation function is free from any distortion of the clustering pattern induced by peculiar motions and is well described by a power-law in pair separation over the range 0.1 < r /h Mpc < 10. The clustering of L* galaxies in real space is well fit by a correlation length r0 = 4.9 +/- 0.3 /h Mpc and power-law slope gamma = 1.71 +/- 0.06. The clustering amplitude increases slowly with absolute magnitude for galaxies fainter than M*, but rises more strongly at higher luminosities. At low luminosities, our results agree with measurements from the SSRS2 by Benoist et al. However, we find a weaker dependence of clustering strength on luminosity at the highest luminosities. The correlation function amplitude increases by a factor of 4.0 between $M_{b_{J}} -5\log_{10}h = -18$ and -22.5, and the most luminous galaxies are 3.0 times more strongly clustered than L* galaxies. The power-law slope of the correlation function shows remarkably little variation for samples spanning a factor of 20 in luminosity. Our measurements are in very good agreement with the predictions of the hierarchical galaxy formation models of Benson et al.

A number of studies have shown that the visibility of scattered broad emission lines in Seyfert 2 galaxies is strongly dependent on the IRAS f_60/f_25 flux ratio, where those Seyfert 2s with `warm' IRAS colours show polarised broad line emission. It is now clear that this effect is due to the increasing dominance of the galactic rather than the AGN emission at 60um in less luminous `cool' Seyfert 2s. However, we present evidence that the 25um emission is a good measure of the AGN luminosity for most Seyfert 2s. Using this result, we show that the visibility of scattered broad line emission has a dependence on the AGN luminosity. The observations can be interpreted self-consistently if the scale height of the scattering zone varies with central source luminosity whilst the scale height of the obscuring torus is approximately constant.

We perform a joint likelihood analysis of the power spectra of the 2dF Galaxy Redshift Survey (2dFGRS) and the cosmic microwave background (CMB) anisotropies under the assumptions that the initial fluctuations were adiabatic, Gaussian and well described by power laws with scalar and tensor indices of n_s and n_t. On its own, the 2dFGRS sets tight limits on the parameter combination Omega_m h, but relatively weak limits on the fraction of the cosmic matter density in baryons Omega_b/Omega_m. The CMB anisotropy data alone set poor constraints on the cosmological constant and Hubble constant because of a `geometrical degeneracy' among parameters. Furthermore, if tensor modes are allowed, the CMB data allow a wide range of values for the physical densities in baryons and cold dark matter. Combining the CMB and 2dFGRS data sets helps to break both the geometrical and tensor mode degeneracies. The values of the parameters derived here are consistent with the predictions of the simplest models of inflation, with the baryon density derived from primordial nucleosynthesis and with direct measurements of the Hubble parameter. In particular, we find strong evidence for a positive cosmological constant with a pm 2sigma range of 0.65 < Omega_Lambda < 0.85, completely independently of constraints on Omega_\Lambda derived from Type Ia supernovae.

This paper continues our study of the behaviour of near infrared helium recombination lines in planetary nebula. We find that the 1.7007um 4^3D-3^3P HeI line is a good measure of the HeI recombination rate, since it varies smoothly with the effective temperature of the central star. We were unable to reproduce the observed data using detailed photoionisation models at both low and high effective temperatures, but plausible explanations for the difference exist for both. We therefore conclude that this line could be used as an indicator of the effective temperature in obscured nebula. We also characterised the nature of the molecular hydrogen emission present in a smaller subset of our sample. The results are consistent with previous data indicating that ultraviolet excitation rather than shocks is the main cause of the molecular hydrogen emission in planetary nebulae.

We report the results of a spectropolarimetric survey of a complete far infrared selected sample of Seyfert 2 galaxies. We have found polarized broad Halpha emission in one new source, NGC5995. In the sample as a whole, there is a clear tendency for galaxies in which we have detected broad Halpha in polarized light to have warm mid--far infrared colours (F(60um)/F(25um)<4), in agreement with our previous results. However, a comparison of the optical, radio and hard x-ray properties of these systems leads us to conclude that this is a secondary consequence of the true mechanism governing our ability to see scattered light from the broad line region. We find a strong trend for galaxies showing such emission to lie above a critical value of the relative luminosity of the active core to the host galaxy (as measured from the [OIII] 5007A equivalent width) which varies as a function of the obscuring column density as measured from hard x-ray observations. The warmth of the infrared colours is then largely due to a combination of the luminosity of the active core, the obscuring column and the relative importance of the host galaxy in powering the far infrared emission, and not solely orientation as we inferred in our previous paper. Our data may also provide an explanation as to why the most highly polarized galaxies, which appear to have tori that are largely edge-on, are also the most luminous and have the most easily detectable scattered broad Halpha.

The 2dF Galaxy Redshift Survey has now measured in excess of 160000 galaxy redshifts. This paper presents the power spectrum of the galaxy distribution, calculated using a direct FFT-based technique. We argue that, within the k-space region 0.02<k<0.15 h Mpc^-1, the shape of this spectrum should be close to that of the linear density perturbations convolved with the window function of the survey. This window function and its convolving effect on the power spectrum estimate are analyzed in detail. By convolving model spectra, we are able to fit the power-spectrum data and provide a measure of the matter content of the universe. Our results show that models containing baryon oscillations are mildly preferred over featureless power spectra. Analysis of the data yields 68% confidence limits on the total matter density times the Hubble parameter \Omega_m h = 0.20 +/- 0.03, and the baryon fraction \Omega_b/\Omega_m = 0.15 +/- 0.07, assuming scale-invariant primordial fluctuations.

The large-scale structure in the distribution of galaxies is thought to arise from the gravitational instability of small fluctuations in the initial density field of the universe. A key test of this hypothesis is that superclusters of galaxies in the process of formation should generate systematic infall of other galaxies. This would be evident in the pattern of recessional velocities, causing an anisotropy in the inferred spatial clustering of galaxies. Here we report a precise measurement of this clustering, using the redshifts of more than 141,000 galaxies from the two-degree-field galaxy redshift survey. We determine the parameter beta = Omega^{0.6}/b = 0.43 +- 0.07, where Omega is the total mass-density parameter and b is a measure of the `bias' of the luminous galaxies in the survey. Combined with the anisotropy of the cosmic microwave background, our results favour a low-density universe with Omega approximately 0.3.

We report on the discovery of a binary pulsar, PSR J1740-3052, during the Parkes multibeam survey. Timing observations of the 570-ms pulsar at Jodrell Bank and Parkes show that it is young, with a characteristic age of 350 kyr, and is in a 231-day, highly eccentric orbit with a companion whose mass exceeds 11 M_sun. An accurate position for the pulsar was obtained using the Australia Telescope Compact Array. Near-infrared 2.2-um observations make with the telescopes at the Siding Spring observatory reveal a late-type star coincident with the pulsar position. However, we do not believe that this star is the pulsar's companion, because a typical star of this spectral type and required mass would extend beyond the pulsar's orbit. Furthermore, the measured advance of periastron of the pulsar suggests a more compact companion, for example, a main-sequence star with radius only a few times that of the sun. Such a companion is also more consistent with the small dispersion measure variations seen near periastron. Although we cannot conclusively rule out a black-hole companion, we believe the companion is probably an early B star, making the system similar to the binary PSR J0045-7319.

We announce here the public availability of the Edinburgh/Durham Southern Galaxy Catalogue (EDSGC, http://www.edsgc.org). This objective galaxy catalogue was constructed using the COSMOS micro-densitometer at the Royal Observatory Edinburgh and constitutes one of the largest digitized galaxy surveys currently in existence. The EDSGC contains a total of 1,495,877 galaxies (each with 27 image parameters) covering an contiguous area of 1182 sq deg centered on the South Galactic Pole. The data consists of photographic bj magnitudes calibrated via CCD sequences which provide a plate-to-plate accuracy of bj~0.1. Extensive external checks have demonstrated that the global EDSGC photometry is free of large-scale systematic gradients and is therefore, ideal for studying the distribution of galaxies on large angular scales. Independent spectroscopy of EDSGC galaxies has shown that the accuracy of the star-galaxy separation is consistent with earlier visual checks and that only 12% of EDSGC galaxies are potentially mis-classified. This paper is intended to provide a summary of the essential details of the catalogue's design and construction as well as provide a brief summary of the main scientific achievements using the EDSGC over the last ten years.

We have observed a large sample of compact planetary nebulae in the near-infrared to determine how the 2^1P-2^1S HeI line at 2.058um varies as a function of stellar effective temperature, Teff. The ratio of this line with HI Br g at 2.166um has often been used as a measure of the highest Teff present in a stellar cluster, and hence on whether there is a cut-off in the stellar initial mass function at high masses. However, recent photoionisation modelling has revealed that the behaviour of this line is more complex than previously anticipated. Our work shows that in most aspects the photoionisation models are correct. In particular, we confirm the weakening of the 2^1P-2^1S as Teff increases beyond 40000K. However, in many cases the model underpredicts the observed ratio when we consider the detailed physical conditions in the individual planetary nebulae. Furthermore, there is evidence that there is still significant 2^1P-2^1S HeI line emission even in the planetary nebulae with very hot (Teff>100000K) central stars. It is clear from our work that this ratio cannot be considered as a reliable measure of effective temperature on its own.

We present near-infrared broad-band polarization images of the nuclear regions of the Circinus galaxy in the J, H and K bands. For the first time the south-eastern reflection cone is detected in polarized light, which is obscured at optical wavelengths behind the galactic disk. This biconical structure is clearly observed in J and H band polarized flux whilst in the K band a more compact structure is detected. Total flux J-K and H-K colour maps reveal a complex colour gradient toward the south-east direction (where the Circinus galactic disk is nearer to us). We find enhanced extinction in an arc shaped structure, at about 200pc from the nucleus, probably part of the star-formation ring. We model the polarized flux images with the scattering and torus model of Young etal, with the same basic input parameters as used by Alexander etal in the spectropolarimetry modelling of Circinus. The best fit to the polarized flux is achieved with a torus radius of ~16pc, and a visual extinction Av, through the torus, to the near-infrared emission regions of >66 mags.

We describe the 2dF Galaxy Redshift Survey (2dFGRS), and the current status of the observations. In this exploratory paper, we apply a Principal Component Analysis to a preliminary sample of 5869 galaxy spectra and use the two most significant components to split the sample into five spectral classes. These classes are defined by considering visual classifications of a subset of the 2dF spectra, and also by comparing to high quality spectra of local galaxies. We calculate a luminosity function for each of the different classes and find that later-type galaxies have a fainter characteristic magnitude, and a steeper faint-end slope. For the whole sample we find M*=-19.7 (for Omega=1, H_0=100 km/sec/Mpc), alpha=-1.3, phi*=0.017. For class 1 (`early-type') we find M*=-19.6, alpha=-0.7, while for class 5 (`late-type') we find M*=-19.0, alpha=-1.7. The derived 2dF luminosity functions agree well with other recent luminosity function estimates.

We present a new observation of the compact HII region, G29.96-0.02, that allows us to compare the velocity structure in the ionised gas and surrounding molecular gas directly. This allows us to remove most of the remaining ambiguity about the nature of this source. In particular, the comparison of the velocity structure present in the 4S-3P HeI lines with that found in the 1-0 S(1) of molecular hydrogern convincingly rules out a bow shock as being important to the kinematics of this source. Our new observation therefore agrees with our previous conclusion, drawn from a velocity resolved HI Br gamma map, that most of the velocity structure in G29.96-0.02 can largely be explained as a result of a champagne flow model. We also find that the best simple model must invoke a powerful stellar wind to evacuate the `head' of the cometary HII region of ionised gas. However, residual differences between model and data tend to indicate that no single simple model can adequately explain all the observed features.

We present the results of a series of observations of the near- and mid-infrared polarisation properties of the Seyfert 2 galaxy NGC1068. Our data agree well with previously published results in showing the need for a separate polarisation mechanism in the near infrared apart from scattering. We find that the predictions of a simple model in which this component arises through absorptive dichroism due to aligned grains within the extended warm (~400K) dust fits the data reasonably if the obscured background source is itself due to dust emission (at temperature >1000K). By considering the change of polarisation with wavelength we show that the extinction to this hot dust region is in the range AV=20-40. Consideration of the observed data then leads us to the conclusion that if viewed face-on, NGC1068 would have a strong near-infrared excess similar to Seyfert 1 galaxies. Comparison with other independent measures of the extinction to the active nucleus itself lead us to the conclusion that the hot dust must provide screening equivalent to at least AV=40, and possibly much higher. We speculate that this component alone may be the `classical' torus discussed in terms of the unified model, and the more extensive mid-infrared emission may arise from circumnuclear molecular cloud material, and dust in the ionisation cones.

The radio jet axis of NGC 1068 is characterised by energetic activity from x-ray to radio wavelengths. Detailed kinematic and polarization studies have shown that this activity is confined to bipolar cones centered on the AGN which intersect the plane of the disk. Thus, molecular clouds at 1 kpc distance along this axis are an important probe of the nuclear ionizing luminosity and spectrum. Extended MIR emission coincident with the clouds is reasonably understood by dust heated to high temperatures by the nuclear radiation field. This model predicts that the nuclear spectrum is quasar-like (power law + blue excess) with a luminosity 2-5 times higher than inferred by Pier et al. Consequently, there is little or no polyaromatic hydrocarbon (PAH) emission associated with the radio-axis molecular clouds. We review this model in the light of new observations. A multi-waveband collage is included to illustrate the possible orientations of the double cones to our line of sight and the galaxian plane.

We have re-examined the nature of the cluster galaxy luminosity function using the data from the Edinburgh-Durham Southern Galaxy Catalogue and the Edinburgh-Milano Redshift Survey. We derive a best fit luminosity function over the range M(bj)=-18 to -21, for a composite sample of 22 of the richer clusters that has M*=-20.16+/-0.02 and alpha=-1.22+/-0.04. The dominant error in these values results from the choice of background subtraction method. From extensive simulations we can show that when the LF is fitted over this narrow range, it is difficult to discriminate against bright values of M*in the single cluster fits, but that faint values provide a strong test of the universality of the luminosity function. We find that all the individual cluster data are well fit by a Schechter function with alpha fixed at -1.25, and that <=10% of these have fitted values of M* that disagree from the average at the 99% confidence level. We further show that fitting only a single parameter Schechter function to composite subsets of the data can give erroneous results for the derived M*. By considering two parameter fits, the results of Monte-Carlo simulations and direct two-sample chi-squared tests we conclude that there is only weak evidence for differences between the data when broken down into subsets based on physical properties (Bautz-Morgan class, richness, velocity dispersion): from our simulations, only the evidence for a difference between subsets based on velocity dispersion may in fact be significant. We find no evidence at all that a Schechter function is not a good model for the intrinsic cluster luminosity function over this absolute magnitude range. Models that invoke strong evolution of galaxy luminosity of all galaxies within clusters are inconsistent with our results.

We present complete, low resolution IJHK spectroscopy of the ultracompact HII region, G45.12+0.13. From the observed HI line strengths, we derive a near infrared extinction law that is slightly steeper than the average. After correction with this extinction law, we find good agreement between the observed line ratios of HeI, Fe+, Fe++, S+ and S++ and the available atomic data. Our data show that the density within the core of G45.12+0.13 must be at least 10^4/cm^3. This is consistent with the known radio structure of the HII region and in considerable disagreement with previous work using mid and far infrared lines. There must also be considerable opacity in the HeI 2 3P-2 3S transition, and we show how the observed strengths of the other HeI lines are consistent with this. From modelling the photoionisation structure, we find good agreement with most of the observed data if the hottest star present has Teff < 42000K. Consideration of the helium ionisation state places a lower limit on this value so that we can also constrain Teff > 38000K. Discrepancies still exist between some of the observed and model line ratios, but the most obvious tend to be the mid-IR observations.

We have mapped the velocity structure of the cometary compact HII region G29.96-0.02 using long-slit echelle spectra of the HI Br gamma line. This technique detects line emission over a much wider area at the necessary spatial resolution compared to radio recombination line observations. Significant structure in both the velocity centroids and the line widths is seen over the entire nebula. Large line widths are seen ahead of the bow and in the tail which may be due to turbulent motions in shocked and interface regions respectively. We construct analytic models of the density and velocity structure in order to attempt to distinguish between the bow shock and champagne flow models which have been put forward to explain the cometary morphology of many compact HII regions. The bow shock model is unable to explain the large velocity gradient that we see right across the tail of the cometary region which can only be explained by the streaming motions towards low density regions in the champagne model. However, our approximation to the champagne model is also not able to fit all of the features of the data. More realistic versions of this model which include the effects of stellar winds and density gradients may be able to provide a better match to these data.

Supernovae remnants are known to be luminous sources of infrared [FeII] emission. By studying how the luminosity scales with age, environment and other relevant factors, we can construct an [FeII] luminosity function for supernovae remnants. This will enable us to predict supernovae rates in starburst galaxies that are too distant for individual remnants to be resolved. First, however, we require accurate luminosities for a sample of remnants of varying ages, and in varying physical environments. As part of this project we have carried out an initial study of a small sample of evolved (ages greater than a few thousand years) remnants in M33. From these data we tentatively conclude that there is evidence for the peak luminosity in the [FeII] lines of these sources to arise in a narrow range of ages. In other respects, the M33 remnants are similar to their galactic and Magellanic Cloud counterparts in the observed peak luminosity. From this, and internal evidence as to the environment present in these regions, we conclude that the luminosity of evolved remnants is only marginally dependent on density and metallicity.

We present the redshifts of the galaxies and galaxy clusters which form the Edinburgh--Milano (EM) cluster redshift survey. A total of 777 galaxy redshifts have been measured in 94 clusters extracted from the digitised Edinburgh Durham Cluster Catalogue. We also present the radial velocity dispersions for 37 clusters. Observational and data reduction techniques are discussed in detail, together with the strategy adopted to determine the mean redshift of a cluster and to identify and discard plausible phantom clusters. Some 10% of our clusters show heavy contamination, indicating that projection is a serious problem for optically selected rich clusters. The median velocity dispersion estimated for a sub--sample of richness R \geq 1 clusters is 742 \pm63 km/sec From a simple comparison with Omega=1 Cold Dark Matter models of structure formation, these results favour a biasing parameter b=1.6-2.0 and are inconsistent with a bias outside of the range b=1.3-2.5.

We present here the results of a statistical search for cluster alignments using the Edinburgh/Milano cluster redshift survey. This survey is a unique cluster database which has been objectively constructed to help minimise the systematic biases associated with previous optical cluster catalogues. We find some evidence for cluster alignments out to spatial separations of $<10\mpc$, however, it is not statistically significant. On larger scales, we find no evidence, statistically significant or not, for cluster alignments. These results are in most disagreement with the recent observations of West and Plionis; both of whom see significant cluster alignments out to $\simeq30\mpc$ and beyond in the Abell \& Lick catalogues of clusters. Our findings are consistent with other searches for cluster alignments that do not involve these catalogues.