We used 3.1 million spectroscopically labelled sources from the Sloan Digital Sky Survey (SDSS) to train an optimised random forest classifier using photometry from the SDSS and the Widefield Infrared Survey Explorer (WISE). We applied this machine learning model to 111 million previously unlabelled sources from the SDSS photometric catalogue which did not have existing spectroscopic observations. Our new catalogue contains 50.4 million galaxies, 2.1 million quasars, and 58.8 million stars. We provide individual classification probabilities for each source, with 6.7 million galaxies (13%), 0.33 million quasars (15%), and 41.3 million stars (70%) having classification probabilities greater than 0.99; and 35.1 million galaxies (70%), 0.72 million quasars (34%), and 54.7 million stars (93%) having classification probabilities greater than 0.9. Precision, Recall, and F1 score were determined as a function of selected features and magnitude error. We investigate the effect of class imbalance on our machine learning model and discuss the implications of transfer learning for populations of sources at fainter magnitudes than the training set. We used a non-linear dimension reduction technique, Uniform Manifold Approximation and Projection (UMAP), in unsupervised, semi-supervised, and fully-supervised schemes to visualise the separation of galaxies, quasars, and stars in a two-dimensional space. When applying this algorithm to the 111 million sources without spectra, it is in strong agreement with the class labels applied by our random forest model.

We present LOFAR observations at 150 MHz of the borderline FRI/FRII giant radio galaxy NGC 6251. This paper presents the most sensitive and highest-resolution images of NGC 6251 at these frequencies to date, revealing for the first time a low-surface-brightness extension to the northern lobe, and a possible backflow associated with the southern lobe. The integrated spectra of components of NGC 6251 are consistent with previous measurements at higher frequencies, similar to results from other LOFAR studies of nearby radio galaxies. We find the outer structures of NGC 6251 to be either at equipartition or slightly electron dominated, similar to those of FRII sources rather than FRIs; but this conclusion remains tentative because of uncertainties associated with the geometry and the extrapolation of X-ray measurements to determine the external pressure distribution on the scale of the outer lobes. We place lower limits on the ages of the extension of the northern lobe and the backflow of the southern lobe of $t \gtrsim 250$ Myr and $t \gtrsim 210$ Myr respectively. We present the first detection of polarisation at 150 MHz in NGC 6251. Taking advantage of the high Faraday resolution of LOFAR, we place an upper limit on the magnetic field in the group of $B < 0.2 (\Lambda_B / 10 {\rm kpc})^{-0.5} \mu$G for a coherence scale of $\Lambda_B < 60 {\rm kpc}$ and $B < 13 \mu$G for $\Lambda_B = 240$ kpc.

We present LOFAR data from 110--180~MHz of the merging galaxy cluster Abell 1682, alongside archival optical, radio and X-ray data. Our 6 arc-second resolution images at low frequencies reveal new structures associated with numerous radio galaxies in the cluster. At 20 arc-second resolution we see diffuse emission throughout the cluster over hundreds of kpc, indicating particle acceleration mechanisms are in play as a result of the cluster merger event and powerful active galactic nuclei. We show that a significant part of the cluster emission is from an old radio galaxy with very steep spectrum emission (having a spectral index of $\alpha < -2.5$). Furthermore we identify a new region of diffuse steep spectrum emission ($\alpha < -1.1$) as a candidate for a radio halo which is co-spatial with the centre of the cluster merger. We suggest its origin as a population of old and mildly relativistic electrons left over from radio galaxies throughout the cluster which have been re-accelerated to higher energies by shocks and turbulence induced by the cluster merger event. We also note the discovery of six new giant radio galaxies in the vicinity of Abell 1682.

The disc structure of the Milky Way is marked by a chemical dichotomy, with high-alpha and low-alpha abundance sequences, traditionally identified with the geometric thick and thin discs. This identification is aided by the old ages of the high-alpha stars, and lower average ages of the low-alpha ones. Recent large scale surveys such as APOGEE have provided a wealth of data on this chemical structure, including showing that an identification of chemical and geometric thick discs is not exact, but the origin of the chemical dichotomy has remained unclear. Here we demonstrate that a dichotomy arises naturally if the early gas-rich disc fragments, leading to some fraction of the star formation occurring in clumps of the type observed in high-redshift galaxies. These clumps have high star formation rate density. They, therefore, enrich rapidly, moving from the low-alpha to the high-alpha sequence, while more distributed star formation produces the low-alpha sequence. We demonstrate that this model produces a chemically-defined thick disc that has many of the properties of the Milky Way's thick disc. Because clump formation is common in high redshift galaxies, we predict that chemical bimodalities are common in massive galaxies.

The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168MHz survey of the entire northern sky for which observations are now 20% complete. We present our first full-quality public data release. For this data release 424 square degrees, or 2% of the eventual coverage, in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45$^\circ$00$'$00$''$ to 57$^\circ$00$'$00$''$) were mapped using a fully automated direction-dependent calibration and imaging pipeline that we developed. A total of 325,694 sources are detected with a signal of at least five times the noise, and the source density is a factor of $\sim 10$ higher than the most sensitive existing very wide-area radio-continuum surveys. The median sensitivity is S$_{\rm 144 MHz} = 71\,\mu$Jy beam$^{-1}$ and the point-source completeness is 90% at an integrated flux density of 0.45mJy. The resolution of the images is 6$''$ and the positional accuracy is within 0.2$''$. This data release consists of a catalogue containing location, flux, and shape estimates together with 58 mosaic images that cover the catalogued area. In this paper we provide an overview of the data release with a focus on the processing of the LOFAR data and the characteristics of the resulting images. In two accompanying papers we provide the radio source associations and deblending and, where possible, the optical identifications of the radio sources together with the photometric redshifts and properties of the host galaxies. These data release papers are published together with a further $\sim$20 articles that highlight the scientific potential of LoTSS.

The shape of low-frequency radio continuum spectra of normal galaxies is not well understood, the key question being the role of physical processes such as thermal absorption in shaping them. In this work we take advantage of the LOFAR Multifrequency Snapshot Sky Survey (MSSS) to investigate such spectra for a large sample of nearby star-forming galaxies. Using the measured 150MHz flux densities from the LOFAR MSSS survey and literature flux densities at various frequencies we have obtained integrated radio spectra for 106 galaxies. The spectra are explained through the use of a three-dimensional model of galaxy radio emission, and radiation transfer dependent on the galaxy viewing angle and absorption processes. Spectra of our galaxies are generally flatter at lower compared to higher frequencies but as there is no tendency for the highly inclined galaxies to have more flattened low-frequency spectra, we argue that the observed flattening is not due to thermal absorption, contradicting the suggestion of Israel & Mahoney (1990). According to our modelled radio maps for M51-like galaxies, the free-free absorption effects can be seen only below 30MHz and in the global spectra just below 20MHz, while in the spectra of starburst galaxies, like M82, the flattening due to absorption is instead visible up to higher frequencies of about 150MHz. Locally, within galactic disks, the absorption effects are distinctly visible in M51-like galaxies as spectral flattening around 100-200MHz in the face-on objects, and as turnovers in the edge-on ones, while in M82-like galaxies there are strong turnovers at frequencies above 700MHz, regardless of viewing angle. Our modelling suggests that the weak spectral flattening observed in the nearby galaxies studied here results principally from synchrotron spectral curvature due to cosmic ray energy losses and propagation effects.

We have observed and analysed the eclipses of the black widow pulsar J1810+1744 at low radio frequencies. Using LOw-Frequency ARray (LOFAR) and Westerbork Synthesis Radio Telescope observations between 2011--2015 we have measured variations in flux density, dispersion measure and scattering around eclipses. High-time-resolution, simultaneous beamformed and interferometric imaging LOFAR observations show concurrent disappearance of pulsations and total flux from the source during the eclipses, with a $3\sigma$ upper limit of 36 mJy ($<10\%$ of the pulsar's averaged out-of-eclipse flux density). The dispersion measure variations are highly asymmetric, suggesting a tail of material swept back due to orbital motion. The egress deviations are variable on timescales shorter than the 3.6 hr orbital period and are indicative of a clumpy medium. Additional pulse broadening detected during egress is typically $<20\%$ of the pulsar's spin period, showing no evidence of scattering the pulses beyond detectability in the beamformed data. The eclipses, lasting $\sim13\%$ of the orbit at 149 MHz, are shown to be frequency-dependent with total duration scaling as $\propto\nu^{-0.41\pm0.03}$. The results are discussed in the context of the physical parameters of the system, and an examination of eclipse mechanisms reveals cyclotron-synchrotron absorption as the most likely primary cause, although non-linear scattering mechanisms cannot be quantitatively ruled out. The inferred mass loss rate is a similar order-of-magnitude to the mean rate required to fully evaporate the companion in a Hubble time.

Feedback by radio-loud active galactic nuclei (AGN) in galaxy groups is not fully understood. Open questions include the duty cycle of the AGN, the spatial extent of the radio lobes, the effect they have on the intragroup medium, and the fate of the cosmic rays. We present the discovery of a 650 kpc-radio galaxy embedded in steep diffuse emission at $z = 0.18793 \pm 5 \times 10^{-5}$ located at the center of the galaxy group MaxBCG J199.31832+51.72503 using an observation from the LOFAR Two-meter Sky Survey (LoTSS) at the central frequency of 144 MHz. Subsequently, we performed a GMRT observation at the central frequency of 607 MHz to study the spectral properties of the source. The observations reveal a radio galaxy with a total radio power $P_{\rm tot, 1.4} \sim 2.1 \times 10^{24}$ W Hz$^{-1}$, exhibiting two asymmetrical jets and lobes. The derived spectral index map shows a steepening toward the inner regions and a steep-spectrum core region. We model the integrated radio spectrum, providing two possible interpretations: the radio source is evolved but still active or it is just at the end of its active phase. Finally, in the same field of view we have discovered Mpc-sized emission surrounding a close pair of AGN located at a redshift $z = 0.0587 \pm 2 \times 10^{-4}$ (SDSS J131544.56+521213.2 and SDSS J131543.99+521055.7) which could be a radio remnant source.

LOFAR observations at 144 MHz have revealed large-scale radio sources in the unrelaxed galaxy cluster Abell 1132. The cluster hosts diffuse radio emission on scales of $\sim$650 kpc near the cluster center and a head-tail (HT) radio galaxy, extending up to 1 Mpc, South of the cluster center. The central diffuse radio emission is not seen in NVSS, FIRST, WENSS, nor in C & D array VLA observations at 1.4 GHz, but is detected in our follow-up GMRT observations at 325 MHz. Using LOFAR and GMRT data, we determine the spectral index of the central diffuse emission to be $\alpha=-1.75\pm0.19$ ($S\propto\nu^{\alpha}$). We classify this emission as an ultra-steep spectrum radio halo and discuss the possible implications for the physical origin of radio halos. The HT radio galaxy shows narrow, collimated emission extending up to 1 Mpc and another 300 kpc of more diffuse, disturbed emission, giving a full projected linear size of 1.3 Mpc - classifying it as a giant radio galaxy (GRG) and making it the longest HT found to date. The head of the GRG coincides with an elliptical galaxy (SDSS J105851.01$+$564308.5) belonging to Abell 1132. In our LOFAR image, there appears to be a connection between the radio halo and the GRG. The turbulence that may have produced the halo may have also affected the tail of the GRG. In turn, the GRG may have provided seed electrons for the radio halo.

We present an extensively updated version of the purely ray-tracing 3D dust radiation transfer code DART-Ray. The new version includes five major upgrades : 1) a series of optimizations for the ray-angular density and the scattered radiation source function; 2) the implementation of several data and task parallelizations using hybrid MPI+OpenMP schemes; 3) the inclusion of dust self-heating; 4) the ability to produce surface brightness maps for observers within the models in HEALPix format; 5) the possibility to set the expected numerical accuracy already at the start of the calculation. We tested the updated code with benchmark models where the dust self-heating is not negligible. Furthermore, we performed a study of the extent of the source influence volumes, using galaxy models, which are critical in determining the efficiency of the DART-Ray algorithm. The new code is publicly available, documented for both users and developers, and accompanied by several programmes to create input grids for different model geometries and to import the results of N-body and SPH simulations. These programmes can be easily adapted to different input geometries, and for different dust models or stellar emission libraries.

Using a cosmological $N$-body simulation, we investigate the origin and distribution of stars in the intracluster light (ICL) of a Fornax-like cluster. In a dark matter only simulation we identify a halo which, at $z=0$, has $M_200 \simeq 4.1 \times 10^{13}M_{sun}$ and $r_{200} = 700kpc$, and replace infalling subhalos with models that include spheroid and disc components. As they fall into the cluster, the stars in some of these galaxies are stripped from their hosts, and form the ICL. We consider the separate contributions to the ICL from stars which originate in the haloes and the discs of the galaxies. We find that disc ICL stars are more centrally concentrated than halo ICL stars. The majority of the disc ICL stars are associated with one initially disc-dominated galaxy that falls to the centre of the cluster and is heavily disrupted, producing part of the cD galaxy. At radial distances greater than 200kpc, well beyond the stellar envelope of the cD galaxy, stars formerly from the stellar haloes of galaxies dominate the ICL. Therefore at large distances, the ICL population is dominated by older stars.

We present a new analysis of the widely used relation between cavity power and radio luminosity in clusters of galaxies with evidence for strong AGN feedback. We study the correlation at low radio frequencies using two new surveys - the First Alternative Data Release of the TIFR GMRT Sky Survey (TGSS ADR1) at 148 MHz and LOFAR's first all-sky survey, the Multifrequency Snapshot Sky Survey (MSSS) at 140 MHz. We find a scaling relation $P_{\rm cav} \propto L_{148}^{\beta}$, with a logarithmic slope of $\beta = 0.51 \pm 0.14$, which is in good agreement with previous results based on data at 327 MHz. The large scatter present in this correlation confirms the conclusion reached at higher frequencies that the total radio luminosity at a single frequency is a poor predictor of the total jet power. We show that including measurements at 148 MHz alone is insufficient to reliably compute the bolometric radio luminosity and reduce the scatter in the correlation. For a subset of four well-resolved sources, we examine the detected extended structures at low frequencies and compare with the morphology known from higher frequency images and Chandra X-ray maps. In Perseus we discuss details in the structures of the radio mini-halo, while in the 2A 0335+096 cluster we observe new diffuse emission associated with multiple X-ray cavities and likely originating from past activity. For A2199 and MS 0735.6+7421, we confirm that the observed low-frequency radio lobes are confined to the extents known from higher frequencies. This new low-frequency analysis highlights the fact that existing cavity power to radio luminosity relations are based on a relatively narrow range of AGN outburst ages. We discuss how the correlation could be extended using low frequency data from the LOFAR Two-metre Sky Survey (LoTSS) in combination with future, complementary deeper X-ray observations.

We report on the discovery in the LOFAR Multifrequency Snapshot Sky Survey (MSSS) of a giant radio galaxy (GRG) with a projected size of $2.56 \pm 0.07$ Mpc projected on the sky. It is associated with the galaxy triplet UGC 9555, within which one is identified as a broad-line galaxy in the Sloan Digital Sky Survey (SDSS) at a redshift of $0.05453 \pm 1 \times 10^{-5} $, and with a velocity dispersion of $215.86 \pm 6.34$ km/s. From archival radio observations we see that this galaxy hosts a compact flat-spectrum radio source, and we conclude that it is the active galactic nucleus (AGN) responsible for generating the radio lobes. The radio luminosity distribution of the jets, and the broad-line classification of the host AGN, indicate this GRG is orientated well out of the plane of the sky, making its physical size one of the largest known for any GRG. Analysis of the infrared data suggests that the host is a lenticular type galaxy with a large stellar mass ($\log~\mathrm{M}/\mathrm{M}_\odot = 11.56 \pm 0.12$), and a moderate star formation rate ($1.2 \pm 0.3~\mathrm{M}_\odot/\mathrm{year}$). Spatially smoothing the SDSS images shows the system around UGC 9555 to be significantly disturbed, with a prominent extension to the south-east. Overall, the evidence suggests this host galaxy has undergone one or more recent moderate merger events and is also experiencing tidal interactions with surrounding galaxies, which have caused the star formation and provided the supply of gas to trigger and fuel the Mpc-scale radio lobes.

The LOFAR Two-metre Sky Survey (LoTSS) is a deep 120-168 MHz imaging survey that will eventually cover the entire Northern sky. Each of the 3170 pointings will be observed for 8 hrs, which, at most declinations, is sufficient to produce ~5arcsec resolution images with a sensitivity of ~0.1mJy/beam and accomplish the main scientific aims of the survey which are to explore the formation and evolution of massive black holes, galaxies, clusters of galaxies and large-scale structure. Due to the compact core and long baselines of LOFAR, the images provide excellent sensitivity to both highly extended and compact emission. For legacy value, the data are archived at high spectral and time resolution to facilitate subarcsecond imaging and spectral line studies. In this paper we provide an overview of the LoTSS. We outline the survey strategy, the observational status, the current calibration techniques, a preliminary data release, and the anticipated scientific impact. The preliminary images that we have released were created using a fully-automated but direction-independent calibration strategy and are significantly more sensitive than those produced by any existing large-area low-frequency survey. In excess of 44,000 sources are detected in the images that have a resolution of 25arcsec, typical noise levels of less than 0.5 mJy/beam, and cover an area of over 350 square degrees in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45d00m00s to 57d00m00s).

Mergers of galaxy clusters are among the most energetic events in the Universe. These events have significant impact on the intra-cluster medium, depositing vast amounts of energy - often in the form of shocks - as well as heavily influencing the properties of the constituent galaxy population. Many clusters have been shown to host large-scale diffuse radio emission, known variously as radio haloes and relics. These sources arise as a result of electron (re-)acceleration in cluster-scale magnetic fields, although the processes by which this occurs are still poorly understood. We present new, deep radio observations of the high-redshift galaxy cluster MACS J0025.4$-$1222, taken with the GMRT at 325 MHz, as well as new analysis of all archival $Chandra$ X-ray observations. We aim to investigate the potential of diffuse radio emission and categorise the radio population of this cluster, which has only been covered previously by shallow radio surveys. We produce low-resolution maps of MACS J0025.4$-$1222 through a combination of uv-tapering and subtracting the compact source population. Radial surface brightness and mass profiles are derived from the $Chandra$ data. We also derive a 2D map of the ICM temperature. For the first time, two sources of diffuse radio emission are detected in MACS J0025.4$-$1222, on linear scales of several hundred kpc. Given the redshift of the cluster and the assumed cosmology, these sources appear to be consistent with established trends in power scaling relations for radio relics. The X-ray temperature map presents evidence of an asymmetric temperature profile and tentative identification of a temperature jump associated with one relic. We classify the pair of diffuse radio sources in this cluster as a pair of radio relics, given their consistency with scaling relations, location toward the cluster outskirts, and the available X-ray data.

Using $N$-body+SPH simulations of galaxies falling into a cluster, we study the evolution of their radial density profiles. When evolved in isolation, galaxies develop a type~II (down-bending) profile. In the cluster, the evolution of the profile depends on the minimum cluster-centric radius the galaxy reaches, which controls the degree of ram pressure stripping. If the galaxy falls to $\sim 50\%$ of the virial radius, then the profile remains type~II, but if the galaxy reaches down to $\sim 20\%$ of the virial radius the break weakens and the profile becomes more type~I like. The velocity dispersions are only slightly increased in the cluster simulations compared with the isolated galaxy; random motion therefore cannot be responsible for redistributing material sufficiently to cause the change in the profile type. Instead we find that the joint action of radial migration driven by tidally-induced spirals and the outside-in quenching of star formation due to ram pressure stripping alters the density profile. As a result, this model predicts a flattening of the age profiles amongst cluster lenticulars with type~I profiles, which can be observationally tested. %

Standard galaxy formation models predict that large-scale double-lobed radio sources, known as DRAGNs, will always be hosted by elliptical galaxies. In spite of this, in recent years a small number of spiral galaxies have also been found to host such sources. These so-called spiral DRAGNs are still extremely rare, with only $\sim 5$ cases being widely accepted. Here we report on the serendipitous discovery of a new spiral DRAGN in data from the Giant Metrewave Radio Telescope (GMRT) at 322 MHz. The host galaxy, MCG+07-47-10, is a face-on late-type Sbc galaxy with distinctive spiral arms and prominent bulge suggesting a high black hole mass. Using WISE infra-red and GALEX UV data we show that this galaxy has a star formation rate of 0.16-0.75 M$_{\odot}$yr$^{-1}$, and that the radio luminosity is dominated by star-formation. We demonstrate that this spiral DRAGN has similar environmental properties to others of this class, but has a comparatively low radio luminosity of $L_{\rm 1.4GHz}$ = 1.12$\times$10$^{22}$ W Hz$^{-1}$, two orders of magnitude smaller than other known spiral DRAGNs. We suggest that this may indicate the existence of a previously unknown low-luminosity population of spiral DRAGNS.

Recent analysis of the SDSS-III/APOGEE Data Release 12 stellar catalogue has revealed that the Milky Way's metallicity distribution function (MDF) changes shape as a function of radius, transitioning from being negatively skewed at small Galactocentric radii to positively skewed at large Galactocentric radii. Using a high resolution, N-body+SPH simulation, we show that the changing skewness arises from radial migration - metal-rich stars form in the inner disk and subsequently migrate to the metal-poorer outer disk. These migrated stars represent a large fraction (> 50%) of the stars in the outer disk; they populate the high metallicity tail of the MDFs and are, in general, more metal-rich than the surrounding outer disk gas. The simulation also reproduces another surprising APOGEE result: the spatially invariant high-[alpha/Fe] MDFs. This arises in the simulation from the migration of a population formed within a narrow range of radii (3.2+/-1.2 kpc) and time (8.8+/-0.6 Gyr ago), rather than from spatially extended star formation in a homogeneous medium at early times. These results point toward the crucial role radial migration has played in shaping our Milky Way.

We present the results of a four-month campaign searching for low-frequency radio transients near the North Celestial Pole with the Low-Frequency Array (LOFAR), as part of the Multifrequency Snapshot Sky Survey (MSSS). The data were recorded between 2011 December and 2012 April and comprised 2149 11-minute snapshots, each covering 175 deg^2. We have found one convincing candidate astrophysical transient, with a duration of a few minutes and a flux density at 60 MHz of 15-25 Jy. The transient does not repeat and has no obvious optical or high-energy counterpart, as a result of which its nature is unclear. The detection of this event implies a transient rate at 60 MHz of 3.9 (+14.7, -3.7) x 10^-4 day^-1 deg^-2, and a transient surface density of 1.5 x 10^-5 deg^-2, at a 7.9-Jy limiting flux density and ~10-minute time-scale. The campaign data were also searched for transients at a range of other time-scales, from 0.5 to 297 min, which allowed us to place a range of limits on transient rates at 60 MHz as a function of observation duration.

We examine metallicities, ages and orbital properties of halo stars in a Milky-Way like disk galaxy formed in the cosmological hydrodynamical MaGICC simulations. Halo stars were either accreted from satellites or they formed in situ in the disk or bulge of the galaxy and were then kicked up into the halo ("in situ/ kicked-up" stars). Regardless of where they formed both types show surprisingly similar orbital properties: the majority of stars of both types are on short-axis tubes with the same sense of rotation as the disk -- implying that a large fraction of satellites are accreted onto the halo with the same sense of angular momentum as the disk.

We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octaves (from 30 to 160 MHz). The broadband frequency coverage, together with the fast survey speed generated by LOFAR's multibeaming capabilities, make MSSS the first survey of the sort anticipated to be carried out with the forthcoming Square Kilometre Array (SKA). Two of the sixteen frequency bands included in the survey were chosen to exactly overlap the frequency coverage of large-area Very Large Array (VLA) and Giant Metrewave Radio Telescope (GMRT) surveys at 74 MHz and 151 MHz respectively. The survey performance is illustrated within the "MSSS Verification Field" (MVF), a region of 100 square degrees centered at J2000 (RA,Dec)=(15h,69deg). The MSSS results from the MVF are compared with previous radio survey catalogs. We assess the flux and astrometric uncertainties in the catalog, as well as the completeness and reliability considering our source finding strategy. We determine the 90% completeness levels within the MVF to be 100 mJy at 135 MHz with 108" resolution, and 550 mJy at 50 MHz with 166" resolution. Images and catalogs for the full survey, expected to contain 150,000-200,000 sources, will be released to a public web server. We outline the plans for the ongoing production of the final survey products, and the ultimate public release of images and source catalogs.

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.

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.

Stable, hydrogen-burning, M dwarf stars comprise about 75% of all stars in the Galaxy. They are extremely long-lived and because they are much smaller in mass than the Sun (between 0.5 and 0.08 MSun), their temperature and stellar luminosity are low and peaked in the red. We have re-examined what is known at present about the potential for a terrestrial planet forming within, or migrating into, the classic liquid-surface-water habitable zone close to an M dwarf star. Observations of protoplanetary disks suggest that planet-building materials are common around M dwarfs, but N-body simulations differ in their estimations of the likelihood of potentially-habitable, wet planets residing within their habitable zones, which are only ~ 1/5 to 1/50 of the width of that for a G star. Particularly in light of the claimed detection of the planets with masses as small as 5.5 and 7.5 MEarth orbiting M stars, there seems no reason to exclude the possibility of terrestrial planets. Tidally locked synchronous rotation within the narrow habitable zone doesn't necessarily lead to atmospheric collapse, and active stellar flaring may not be as much of an evolutionarily disadvantageous factor as has previously been supposed. We conclude that M dwarf stars may indeed be viable hosts for planets on which the origin and evolution of life can occur. A number of planetary processes such as cessation of geothermal activity, or thermal and non-thermal atmospheric loss processes may limit the duration of planetary habitability to periods far shorter than the extreme lifetime of the M dwarf star. Nevertheless, it makes sense to include M dwarf stars in programs that seek to find habitable worlds and evidence of life.

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.

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.