This white paper describes the science case for Very Long Baseline Interferometry (VLBI) and provides suggestions towards upgrade paths for the European VLBI Network (EVN). The EVN is a distributed long-baseline radio interferometric array, that operates at the very forefront of astronomical research. Recent results, together with the new science possibilities outlined in this vision document, demonstrate the EVN's potential to generate new and exciting results that will transform our view of the cosmos. Together with e-MERLIN, the EVN provides a range of baseline lengths that permit unique studies of faint radio sources to be made over a wide range of spatial scales. The science cases are reviewed in six chapters that cover the following broad areas: cosmology, galaxy formation and evolution, innermost regions of active galactic nuclei, explosive phenomena and transients, stars and stellar masers in the Milky Way, celestial reference frames and space applications. The document concludes with identifying the synergies with other radio, as well as multi-band/multi-messenger instruments, and provide the recommendations for future improvements. The appendices briefly describe other radio VLBI arrays, the technological framework for EVN developments, and a selection of spectral lines of astrophysical interest below 100 GHz. The document includes a glossary for non-specialists, and a list of acronyms at the end.

We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

The Vega planetary system hosts the archetype of extrasolar Kuiper belts, and is rich in dust from the sub-au region out to 100's of au, suggesting intense dynamical activity. We present ALMA mm observations that detect and resolve the outer dust belt from the star for the first time. The interferometric visibilities show that the belt can be fit by a Gaussian model or by power-law models with a steep inner edge (at 60-80 au). The belt is very broad, extending out to at least 150-200 au. We strongly detect the star and set a stringent upper limit to warm dust emission previously detected in the infrared. We discuss three scenarios that could explain the architecture of Vega's planetary system, including the new {ALMA} constraints: no outer planets, a chain of low-mass planets, and a single giant planet. The planet-less scenario is only feasible if the outer belt was born with the observed sharp inner edge. If instead the inner edge is currently being truncated by a planet, then the planet must be $\gtrsim$6 M$_{\oplus}$ and at $\lesssim71$ au to have cleared its chaotic zone within the system age. In the planet chain scenario, outward planet migration and inward scattering of planetesimals could produce the hot and warm dust observed in the inner regions of the system. In the single giant planet scenario, an asteroid belt could be responsible for the warm dust, and mean motion resonances with the planet could put asteroids on star-grazing orbits, producing the hot dust.

We present the highest resolution single-dish submillimetre observations of the detached shell source U Antliae to date. The observations were obtained at $450~\micron$ and $850~\micron$ with SCUBA-2 instrument on the James Clerk Maxwell Telescope as part of the Nearby Evolved Stars Survey. The emission at $850~\micron$ peaks at $40\arcsec$ with hints of a second peak seen at $\sim 20\arcsec$. The emission can be traced out to a radius of $56\arcsec$ at a $3\sigma$ level. The outer peak observed at $850~\micron$ aligns well with the peak observed at Herschel/PACS wavelengths. With the help of spectral energy distribution fitting and radiative transfer calculations of multiple-shell models for the circumstellar envelope, we explore the various shell structures and the variation of grain sizes along the in the circumstellar envelope. We determine a total shell dust mass of $(2.0 \pm 0.3) \times 10^{-5}$ M$_{\odot}$ and established that the thermal pulse which gave rise to the detached shell occurred 3500 $\pm$ 500 years ago.

We present the POL-2 850 $\mu$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $\mu$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $\mu$G and a criticality criterion $\lambda_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component.

We present the 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.

Consistent modeling of protoplanetary disks requires the simultaneous solution of both continuum and line radiative transfer, heating/cooling balance between dust and gas and, of course, chemistry. Such models depend on panchromatic observations that can provide a complete description of the physical and chemical properties and energy balance of protoplanetary systems. Along these lines we present a homogeneous, panchromatic collection of data on a sample of 85 T Tauri and Herbig Ae objects for which data cover a range from X-rays to centimeter wavelengths. Datasets consist of photometric measurements, spectra, along with results from the data analysis such as line fluxes from atomic and molecular transitions. Additional properties resulting from modeling of the sources such as disc mass and shape parameters. dust size and PAH properties are also provided for completeness. Targets were selected based on their properties data availability. Data from more than 50 different telescopes and facilities were retrieved and combined in homogeneous datasets directly from public data archives or after being extracted from more than 100 published articles. X-ray data for a subset of 56 sources represent an exception as they were reduced from scratch and are presented here for the first time. Compiled datasets along with a subset of continuum and emission-line models are stored in a dedicated database and distributed through a publicly accessible online system. All datasets contain metadata descriptors that allow to backtrack them to their original resources. The graphical user interface of the online system allows the user to visually inspect individual objects but also compare between datasets and models. It also offers to the user the possibility to download any of the stored data and metadata for further processing.

We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.

The European FP7 project DIANA has performed a coherent analysis of a large set of observations from protoplanetary disks by means of thermo-chemical disk models. The collected data include extinction-corrected stellar UV and X-ray input spectra (as seen by the disk), photometric fluxes, low and high resolution spectra, interferometric data, emission line fluxes, line velocity profiles and line maps. We define and apply a standardized modelling procedure to simultaneously fit all these data by state-of-the-art modelling codes (ProDiMo, MCFOST, MCMax) which solve the continuum and line radiative transfer, disk chemistry, and the heating & cooling balance for both the gas and the dust. We allow for up to two radial disk zones to obtain our best-fitting models that have about 20 free parameters. This approach is novel and unique in its completeness and level of consistency. In this paper, we present the results from pure SED fitting for 27 objects and from the all inclusive DIANA-standard models for 14 objects. We fit most infrared to millimeter emission line fluxes within a factor better than 3, simultaneously with SED, PAH features and radial brightness profiles extracted from images at various wavelengths. Our analysis shows a number of Herbig Ae and T Tauri stars with very cold and massive outer disks which are situated at least partly in the shadow of a tall and gas-rich inner disk. The disk masses derived are often in excess to previously published values, since these disks are partially optically thick even at millimeter wavelength and so cold that they emit less than in the Rayleigh-Jeans limit. Some line observations cannot be reproduced by the models, probably caused by foreground cloud absorption or object variability. Our data collection, the fitted physical disk parameters as well as the full model output are available at an online database (http://www.univie.ac.at/diana).

(No abstract for this journals: article commences: ) Fairness is a key issue in the careers of astronomers. I examine here the anecdota l suggestion that "you're more likely to get time if you're on the TAC", using public and published data for a large international telescope facility...

Debris disks, comprised of planetsimal belts and the dust and gas produced by their mutual collisions, are the longest-lived phase of circumstellar disks. Typically much fainter in emission than protoplanetary disks, debris disks can be found in associations in which some members still host protoplanetary disks rich in gas as well as around main sequence stars of all ages. There are even classes of debris disks around AGB stars and white dwarfs. Typically, these disks have been studied through dust emission though an increasing number of young disks are now found to host some gas which may be remnants of the protoplanetary disks or second-generation gas. The ngVLA will have a particular niche in the study of the dust population of these disks, since the size distribution of the dust can be derived from the spectral index of the spectral energy distribution from the far-infrared to the centimetre. With the longest lever arms provided by ngVLA, the size distribution can be characterized, testing models of collisional evolution. The ngVLA could also be capable of measuring the quantities of HI and OH emission associated with the disks. HI and OH are by-products of the dissociation of water, which is expected to be a primary component of outgassed material from cometary collisions. A probe of the water content of extrasolar systems has important implications for the delivery of water to terrestrial, potentially Earth-like planets.

Magnetic field is one of the key agents that play a crucial role in shaping molecular clouds and regulating star formation, yet the complete information on the magnetic field is not well constrained due to the limitations in observations. We study the magnetic field in the massive infrared dark cloud G035.39-00.33 from dust continuum polarization observations at 850 $\micron$ with SCUBA-2/POL-2 at JCMT. The magnetic field tends to be perpendicular to the densest part of the main filament (F$_{M}$), whereas it has a less defined relative orientation in the rest of the structure, where it tends to be parallel to some diffuse regions. A mean plane-of-the-sky magnetic field strength of $\sim$50 $\mu$G for F$_{M}$ is obtained using Davis-Chandrasekhar-Fermi method. Based on $^{13}$CO (1-0) line observations, we suggest a formation scenario of F$_{M}$ due to large-scale ($\sim$10 pc) cloud-cloud collision. Using additional NH$_3$ line data, we estimate that F$_{M}$ will be gravitationally unstable if it is only supported by thermal pressure and turbulence. The northern part of F$_{M}$, however, can be stabilized by a modest additional support from the local magnetic field. The middle and southern parts of F$_{M}$ are likely unstable even if the magnetic field support is taken into account. We claim that the clumps in F$_{M}$ may be supported by turbulence and magnetic fields against gravitational collapse. Finally, we identified for the first time a massive ($\sim$200 M$_{\sun}$), collapsing starless clump candidate, "c8", in G035.39-00.33. The magnetic field surrounding "c8" is likely pinched, hinting at an accretion flow along the filament.

Several interstellar environments produce 'anomalous microwave emission', with brightness-peaks at tens-of-gigahertz frequencies. The emission's origins are uncertain - rapidly-spinning nano-particles could emit electric-dipole radiation, but polycyclic aromatic hydrocarbons proposed as the carrier are now found not to correlate with Galactic signals. The difficulty is to identify co-spatial sources over long lines of sight. Here we identify anomalous microwave emission in three proto-planetary discs. These are the only known systems that host hydrogenated nano-diamonds, in contrast to very common detection of polycyclic aromatic hydrocarbons. Spectroscopy locates the nano-diamonds close to the host-stars, at physically-constrained temperatures. Developing disc models, we reproduce the emission with diamonds 0.75-1.1 nanometres in radius, holding less than or equal to 1-2 per cent of the carbon budget. The microwave-emission:stellar-luminosity ratios are approximately constant, allowing nano-diamonds to be ubiquitous but emitting below detection thresholds in many star-systems. This can unify the findings with similar-sized diamonds found within solar system meteorites. As nano-diamond spectral absorption is seen in interstellar sightlines, these particles are also a candidate for generating galaxy-scale anomalous microwave emission.

We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $\mu$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $\lambda$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical.

We present high spatial resolution observations of the continuum emission from the young multiple star system UZ Tau at frequencies from 6 to 340 GHz. To quantify the spatial variation of dust emission in the UZ Tau E circumbinary disk, the observed interferometric visibilities are modeled with a simple parametric prescription for the radial surface brightnesses at each frequency. We find evidence that the spectrum steepens with radius in the disk, manifested as a positive correlation between the observing frequency and the radius that encircles a fixed fraction of the emission ($R_{eff} \propto \nu^{0.34 \pm 0.08}$). The origins of this size--frequency relation are explored in the context of a theoretical framework for the growth and migration of disk solids. While that framework can reproduce a similar size--frequency relation, it predicts a steeper spectrum than is observed. Moreover, it comes closest to matching the data only on timescales much shorter ($\le 1$ Myr) than the putative UZ Tau age (~2-3 Myr). These discrepancies are the direct consequences of the rapid radial drift rates predicted by models of dust evolution in a smooth gas disk. One way to mitigate that efficiency problem is to invoke small-scale gas pressure modulations that locally concentrate drifting solids. If such particle traps reach high continuum optical depths at 30-340 GHz with a ~30-60% filling fraction in the inner disk ($r \lesssim20$ au), they can also explain the observed spatial gradient in the UZ Tau E disk spectrum.

We present 850 $\mu$m imaging polarimetry data of the $\rho$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $\rho$ Oph A core; some of them can be part of a coherent magnetic field structure in the $\rho$ Oph region. The results are consistent with previous observations of the brightest regions of $\rho$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components.

We present 850 and 450 micron observations of the dense regions within the Auriga-California molecular cloud using SCUBA-2 as part of the JCMT Gould Belt Legacy Survey to identify candidate protostellar objects, measure the masses of their circumstellar material (disk and envelope), and compare the star formation to that in the Orion A molecular cloud. We identify 59 candidate protostars based on the presence of compact submillimeter emission, complementing these observations with existing Herschel/SPIRE maps. Of our candidate protostars, 24 are associated with young stellar objects (YSOs) in the Spitzer and Herschel/PACS catalogs of 166 and 60 YSOs, respectively (177 unique), confirming their protostellar nature. The remaining 35 candidate protostars are in regions, particularly around LkHalpha 101, where the background cloud emission is too bright to verify or rule out the presence of the compact 70 micron emission that is expected for a protostellar source. We keep these candidate protostars in our sample but note that they may indeed be prestellar in nature. Our observations are sensitive to the high end of the mass distribution in Auriga-Cal. We find that the disparity between the richness of infrared star forming objects in Orion A and the sparsity in Auriga-Cal extends to the submillimeter, suggesting that the relative star formation rates have not varied over the Class II lifetime and that Auriga-Cal will maintain a lower star formation efficiency.

We present Atacama Large Millimeter/sub-millimeter Array (ALMA) Cycle 2 observations of the 1.3 mm dust continuum emission of the protoplanetary disc surrounding the T Tauri star Elias 24 with an angular resolution of $\sim 0.2"$ ($\sim 28$ au). The dust continuum emission map reveals a dark ring at a radial distance of $0.47"$ ($\sim 65$ au) from the central star, surrounded by a bright ring at $0.58"$ ($\sim 81$ au). In the outer disc, the radial intensity profile shows two inflection points at $0.71"$ and $0.87"$ ($\sim 99$ and $121$ au respectively). We perform global three-dimensional smoothed particle hydrodynamic gas/dust simulations of discs hosting a migrating and accreting planet. Combining the dust density maps of small and large grains with three dimensional radiative transfer calculations, we produce synthetic ALMA observations of a variety of disc models in order to reproduce the gap- and ring-like features observed in Elias 24. We find that the dust emission across the disc is consistent with the presence of an embedded planet with a mass of $\sim 0.7\, \mathrm{M_{\mathrm{J}}}$ at an orbital radius of $\sim$ 60 au. Our model suggests that the two inflection points in the radial intensity profile are due to the inward radial motion of large dust grains from the outer disc. The surface brightness map of our disc model provides a reasonable match to the gap- and ring-like structures observed in Elias 24, with an average discrepancy of $\sim$ 5% of the observed fluxes around the gap region.

Saturn's moon Enceladus has vents emerging from a sub-surface ocean, offering unique probes into the liquid environment. These vents drain into the larger neutral torus in orbit around Saturn. We present a methanol (CH3OH) detection observed with IRAM 30-m from 2008 along the line-of-sight through Saturn's E-ring. Additionally, we also present supporting observations from the Herschel public archive of water (ortho-H2O; 1669.9 GHz) from 2012 at a similar elongation and line-of-sight. The CH3OH 5(1,1)-4(1,1) transition was detected at 5.9 sigma confidence. The line has 0.43 km/s width and is offset by +8.1 km/s in the moon's reference frame. Radiative transfer models allow for gas cloud dimensions from 1750 km up to the telescope beam diameter ~73000 km. Taking into account the CH3OH lifetime against solar photodissociation and the redshifted line velocity, there are two possible explanations for the CH3OH emission: methanol is primarily a secondary product of chemical interactions within the neutral torus that (1) spreads outward throughout the E-ring or (2) originates from a compact, confined gas cloud lagging Enceladus by several km/s. We find either scenario to be consistent with significant redshifted H2O emission (4 sigma) measured from the Herschel public archive. The measured CH3OH:H2O abundance (> 0.5 per cent) significantly exceeds the observed abundance in the direct vicinity of the vents (~0.01 per cent), suggesting CH3OH is likely chemically processed within the gas cloud with methane (CH4) as its parent species.

We present the first ALMA observations of the closest known extrasolar debris disc. This disc orbits the star $\epsilon$ Eridani, a K-type star just 3.2pc away. Due to the proximity of the star, the entire disc cannot fit within the ALMA field of view. Therefore, the observations have been centred 18" North of the star, providing us with a clear detection of the northern arc of the ring, at a wavelength of 1.3mm. The observed disc emission is found to be narrow with a width of just 11-13AU. The fractional disc width we find is comparable to that of the Solar System's Kuiper Belt and makes this one of the narrowest debris discs known. If the inner and outer edges are due to resonances with a planet then this planet likely has a semi-major axis of 48AU. We find tentative evidence for clumps in the ring, although there is a strong chance that at least one is a background galaxy. We confirm, at much higher significance, the previous detection of an unresolved emission at the star that is above the level of the photosphere and attribute this excess to stellar chromospheric emission.

Debris discs are evidence of the ongoing destructive collisions between planetesimals, and their presence around stars also suggests that planets exist in these systems. In this paper, we present submillimetre images of the thermal emission from debris discs that formed the SCUBA-2 Observations of Nearby Stars (SONS) survey, one of seven legacy surveys undertaken on the James Clerk Maxwell telescope between 2012 and 2015. The overall results of the survey are presented in the form of 850 microns (and 450 microns, where possible) images and fluxes for the observed fields. Excess thermal emission, over that expected from the stellar photosphere, is detected around 49 stars out of the 100 observed fields. The discs are characterised in terms of their flux density, size (radial distribution of the dust) and derived dust properties from their spectral energy distributions. The results show discs over a range of sizes, typically 1-10 times the diameter of the Edgeworth-Kuiper Belt in our Solar System. The mass of a disc, for particles up to a few millimetres in size, is uniquely obtainable with submillimetre observations and this quantity is presented as a function of the host stars' age, showing a tentative decline in mass with age. Having doubled the number of imaged discs at submillimetre wavelengths from ground-based, single dish telescope observations, one of the key legacy products from the SONS survey is to provide a comprehensive target list to observe at high angular resolution using submillimetre/millimetre interferometers (e.g., ALMA, SMA).

We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments.

The presence of dusty debris around main sequence stars denotes the existence of planetary systems. Such debris disks are often identified by the presence of excess continuum emission at infrared and (sub-)millimetre wavelengths, with measurements at longer wavelengths tracing larger and cooler dust grains. The exponent of the slope of the disk emission at sub-millimetre wavelengths, `q', defines the size distribution of dust grains in the disk. This size distribution is a function of the rigid strength of the dust producing parent planetesimals. As part of the survey `PLAnetesimals around TYpical Pre-main seqUence Stars' (PLATYPUS) we observed six debris disks at 9-mm using the Australian Telescope Compact Array. We obtain marginal (~3-\sigma) detections of three targets: HD 105, HD 61005, and HD 131835. Upper limits for the three remaining disks, HD20807, HD109573, and HD109085, provide further constraint of the (sub-)millimetre slope of their spectral energy distributions. The values of q (or their limits) derived from our observations are all smaller than the oft-assumed steady state collisional cascade model (q = 3.5), but lie well within the theoretically expected range for debris disks q ~ 3 to 4. The measured q values for our targets are all < 3.3, consistent with both collisional modelling results and theoretical predictions for parent planetesimal bodies being `rubble piles' held together loosely by their self-gravity.

We present 450 and 850 micron submillimetre continuum observations of the IC5146 star-forming region taken as part of the JCMT Gould Belt Survey. We investigate the location of bright submillimetre (clumped) emission with the larger-scale molecular cloud through comparison with extinction maps, and find that these denser structures correlate with higher cloud column density. Ninety-six individual submillimetre clumps are identified using FellWalker and their physical properties are examined. These clumps are found to be relatively massive, ranging from 0.5to 116 MSun with a mean mass of 8 MSun and a median mass of 3.7 MSun. A stability analysis for the clumps suggest that the majority are (thermally) Jeans stable, with M/M_J < 1. We further compare the locations of known protostars with the observed submillimetre emission, finding that younger protostars, i.e., Class 0 and I sources, are strongly correlated with submillimetre peaks and that the clumps with protostars are among the most Jeans unstable. Finally, we contrast the evolutionary conditions in the two major star-forming regions within IC5146: the young cluster associated with the Cocoon Nebula and the more distributed star formation associated with the Northern Streamer filaments. The Cocoon Nebula appears to have converted a higher fraction of its mass into dense clumps and protostars, the clumps are more likely to be Jeans unstable, and a larger fraction of these remaining clumps contain embedded protostars. The Northern Streamer, however, has a larger number of clumps in total and a larger fraction of the known protostars are still embedded within these clumps.

The recent progress in instrumentation and telescope development has brought us different ways to observe protoplanetary disks, including interferometers, space missions, adaptive optics, polarimetry, and time- and spectrally-resolved data. While the new facilities have changed the way we can tackle the existing open problems in disk structure and evolution, there is a substantial lack of interconnection between different observing techniques and their user communities. Here, we explore the complementarity of some of the state-of-the-art observing techniques, and how they can be brought together in a collective effort to understand how disks evolve and disperse at the time of planet formation. This paper was born at the "Protoplanetary Discussions" meeting in Edinburgh, 2016. Its goal is to clarify where multi-wavelength observations of disks converge in unveiling disk structure and evolution, and where they diverge and challenge our current understanding. We discuss caveats that should be considered when linking results from different observations, or when drawing conclusions based on limited datasets (in terms of wavelength or sample). We focus on disk properties that are currently being revolutionized by multi-wavelength observations. Specifically: the inner disk radius, holes and gaps and their link to large-scale disk structures, the disk mass, and the accretion rate. We discuss how the links between them, as well as the apparent contradictions, can help us to disentangle the disk physics and to learn about disk evolution.

Gravitational forces are expected to excite spiral density waves in protoplanetary disks, disks of gas and dust orbiting young stars. However, previous observations that showed spiral structure were not able to probe disk midplanes, where most of the mass is concentrated and where planet formation takes place. Using the Atacama Large Millimeter/submillimeter Array we detected a pair of trailing symmetric spiral arms in the protoplanetary disk surrounding the young star Elias 2-27. The arms extend to the disk outer regions and can be traced down to the midplane. These millimeter-wave observations also reveal an emission gap closer to the star than the spiral arms. We argue that the observed spirals trace shocks of spiral density waves in the midplane of this young disk.

We present observations of the Cepheus Flare obtained as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Legacy Survey (GBLS) with the SCUBA-2 instrument. We produce a catalogue of sources found by SCUBA-2, and separate these into starless cores and protostars. We determine masses and densities for each of our sources, using source temperatures determined by the Herschel Gould Belt Survey. We compare the properties of starless cores in four different molecular clouds: L1147/58, L1172/74, L1251 and L1228. We find that the core mass functions for each region typically show shallower-than-Salpeter behaviour. We find that L1147/58 and L1228 have a high ratio of starless cores to Class II protostars, while L1251 and L1174 have a low ratio, consistent with the latter regions being more active sites of current star formation, while the former are forming stars less actively. We determine that, if modelled as thermally-supported Bonnor-Ebert spheres, most of our cores have stable configurations accessible to them. We estimate the external pressures on our cores using archival $^{13}$CO velocity dispersion measurements and find that our cores are typically pressure-confined, rather than gravitationally bound. We perform a virial analysis on our cores, and find that they typically cannot be supported against collapse by internal thermal energy alone, due primarily to the measured external pressures. This suggests that the dominant mode of internal support in starless cores in the Cepheus Flare is either non-thermal motions or internal magnetic fields.

Recent ALMA observations unveiled the structure of CO gas in the 23 Myr-old $\beta$ Pictoris planetary system, a component that has been discovered in many similarly young debris disks. We here present ALMA CO J=2-1 observations, at an improved spectro-spatial resolution and sensitivity compared to previous CO J=3-2 observations. We find that 1) the CO clump is radially broad, favouring the resonant migration over the giant impact scenario for its dynamical origin, 2) the CO disk is vertically tilted compared to the main dust disk, at an angle consistent with the scattered light warp. We then use position-velocity diagrams to trace Keplerian radii in the orbital plane of the disk. Assuming a perfectly edge-on geometry, this shows a CO scale height increasing with radius as $R^{0.75}$, and an electron density (derived from CO line ratios through NLTE analysis) in agreement with thermodynamical models. Furthermore, we show how observations of optically thin line ratios can solve the primordial versus secondary origin dichotomy in gas-bearing debris disks. As shown for $\beta$ Pictoris, subthermal (NLTE) CO excitation is symptomatic of H$_2$ densities that are insufficient to shield CO from photodissociation over the system's lifetime. This means that replenishment from exocometary volatiles must be taking place, proving the secondary origin of the disk. In this scenario, assuming steady state production/destruction of CO gas, we derive the CO+CO$_2$ ice abundance by mass in $\beta$ Pic's exocomets to be at most $\sim$6%, consistent with comets in our own Solar System and in the coeval HD181327 system.

LkCa 15 is an extensively studied star in the Taurus region known for its pre-transitional disk with a large inner cavity in dust continuum and normal gas accretion rate. The most popular hypothesis to explain the LkCa 15 data invokes one or more planets to carve out the inner cavity, while gas continues to flow across the gap from the outer disk onto the central star. We present spatially unresolved HCO+ J=4-3 observations of the LkCa 15 disk from the JCMT and model the data with the ProDiMo code. We find that: (1) HCO+ line-wings are clearly detected, certifying the presence of gas in the cavity within <50 AU of the star. (2) Reproducing the observed line-wing flux requires both a significant suppression of cavity dust (by a factor >10^4 compared to the ISM) and a substantial increase in the gas scale-height within the cavity (H_0/R_0 ~ 0.6). An ISM dust-to-gas ratio (d:g=10^-2) yields too little line-wing flux regardless of the scale-height or cavity gas geometry, while a smaller scale-height also under predicts the flux even with a reduced d:g. (3) The cavity gas mass is consistent with the surface density profile of the outer disk extended inwards to the sublimation radius (corresponding to mass M_d ~ 0.03 M_sun), and masses lower by a factor >10 appear to be ruled out.

Brown dwarf disks are excellent laboratories to test our understanding of disk physics in an extreme parameter regime. In this paper we investigate a sample of 29 well-characterized brown dwarfs and very low mass stars, for which Herschel far-infrared fluxes as well as (sub)-mm fluxes are available. We have measured new Herschel PACS fluxes for 11 objects and complement these with (sub)-mm data and Herschel fluxes from the literature. We analyze their spectral energy distributions in comparison with results from radiative transfer modeling. Fluxes in the far-infrared are strongly affected by the shape and temperature of the disk (and hence stellar luminosity), whereas the (sub)-mm fluxes mostly depend on disk mass. Nevertheless, there is a clear correlation between far-infrared and (sub)-mm fluxes. We argue that the link results from the combination of the stellar mass-luminosity relation and a scaling between disk mass and stellar mass. We find strong evidence of dust settling to the disk midplane. The spectral slopes between near- and far-infrared are mostly between $-0.5$ and $-1.2$ in our sample, comparable to more massive T Tauri stars, which may imply that the disk shapes are similar as well, though highly-flared disks are rare among brown dwarfs. We find that dust temperatures in the range of 7-15 K, calculated with $T\approx25\,(L/L_\odot)^{0.25}$ K, are appropriate for deriving disk masses from (sub)-mm fluxes for these low luminosity objects. About half of our sample hosts disks with at least one Jupiter mass, confirming that many brown dwarfs harbour sufficient material for the formation of Earth-mass planets in their midst.

Discs of dusty debris around main-sequence star indicate fragmentation of orbiting planetesimals, and for a few A-type stars, a gas component is also seen that may come from collisionally-released volatiles. Here we find the sixth example of a CO-hosting disc, around the 30Myr old A0-star HD 32297. Two more of these CO-hosting stars, HD 21997 and 49 Cet, have also been imaged in dust with SCUBA-2 within the SONS project. A census of 27 A-type debris hosts within 125 pc now shows 7/16 detections of carbon-bearing gas within the 5-50 Myr epoch, with no detections in 11 older systems. Such a prolonged period of high fragmentation rates corresponds quite well to the epoch when most of the Earth was assembled from planetesimal collisions. Recent models propose that collisional products can be spatially asymmetric if they originate at one location in the disc, with CO particularly exhibiting this behaviour as it can photodissociate in less than an orbital period. Of the six CO-hosting systems, only beta Pic is in clear support of this hypothesis. However, radiative transfer modelling with the ProDiMo code shows that the CO is also hard to explain in a proto-planetary disc context.

We present the JCMT Gould Belt Survey's first look results of the southern extent of the Orion A Molecular Cloud ($\delta \leq -5\mathrm{:}31\mathrm{:}27.5$). Employing a two-step structure identification process, we construct individual catalogues for large-scale regions of significant emission labelled as islands and smaller-scale subregions called fragments using the 850 $\mu$m continuum maps obtained using SCUBA-2. We calculate object masses, sizes, column densities, and concentrations. We discuss fragmentation in terms of a Jeans instability analysis and highlight interesting structures as candidates for follow-up studies. Furthermore, we associate the detected emission with young stellar objects (YSOs) identified by Spitzer and Herschel. We find that although the population of active star-forming regions contains a wide variety of sizes and morphologies, there is a strong positive correlation between the concentration of an emission region and its calculated Jeans instability. There are, however, a number of highly unstable subregions in dense areas of the map that show no evidence of star formation. We find that only $\sim$72\% of the YSOs defined as Class 0+I and flat-spectrum protostars coincide with dense 850 $\mu$m emission structures (column densities $>3.7\times10^{21}\mathrm{\:cm}^{-2}$). The remaining 28\% of these objects, which are expected to be embedded in dust and gas, may be misclassified. Finally, we suggest that there is an evolution in the velocity dispersion of young stellar objects such that sources which are more evolved are associated with higher velocities.

Epsilon Eridani is a nearby, young Sun-like star that hosts a ring of cool debris analogous to the solar system's Edgeworth-Kuiper belt. Early observations at (sub-)mm wavelengths gave tentative evidence of the presence of inhomogeneities in the ring, which have been ascribed to the effect of a putative low eccentricity planet, orbiting close to the ring. The existence of these structures have been recently challenged by high resolution interferometric millimeter observations. Here we present the deepest single-dish image of Epsilon Eridani at millimeter wavelengths, obtained with the Large Millimeter Telescope Alfonso Serrano (LMT). The main goal of these LMT observations is to confirm (or refute) the presence of non-axisymmetric structure in the disk. The dusty ring is detected for the first time along its full projected elliptical shape. The radial extent of the ring is not spatially resolved and shows no evidence, to within the uncertainties, of dust density enhancements. Additional features of the 1.1 mm map are: (i) the presence of significant flux in the gap between the ring and the star, probably providing the first exo-solar evidence of Poynting-Robertson drag, (ii) an unambiguous detection of emission at the stellar position with a flux significantly above that expected from Epsilon Eridani's photosphere, and (iii) the identification of numerous unresolved sources which could correspond to background dusty star-forming galaxies.

The dust emissivity spectral index, $\beta$, is a critical parameter for deriving the mass and temperature of star-forming structures, and consequently their gravitational stability. The $\beta$ value is dependent on various dust grain properties, such as size, porosity, and surface composition, and is expected to vary as dust grains evolve. Here we present $\beta$, dust temperature, and optical depth maps of the star-forming clumps in the Perseus Molecular Cloud determined from fitting SEDs to combined Herschel and JCMT observations in the 160 $\mu$m, 250 $\mu$m, 350 $\mu$m, 500 $\mu$m, and 850 $\mu$m bands. Most of the derived $\beta$, and dust temperature values fall within the ranges of 1.0 - 2.7 and 8 - 20 K, respectively. In Perseus, we find the $\beta$ distribution differs significantly from clump to clump, indicative of grain growth. Furthermore, we also see significant, localized $\beta$ variations within individual clumps and find low $\beta$ regions correlate with local temperature peaks, hinting at the possible origins of low $\beta$ grains. Throughout Perseus, we also see indications of heating from B stars and embedded protostars, as well evidence of outflows shaping the local landscape.

We present SCUBA-2 450{\mu}m and 850{\mu}m observations of the W40 complex in the Serpens-Aquila region as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) of nearby star-forming regions. We investigate radiative heating by constructing temperature maps from the ratio of SCUBA-2 fluxes using a fixed dust opacity spectral index, {\beta} = 1.8, and a beam convolution kernel to achieve a common 14.8" resolution. We identify 82 clumps ranging between 10 and 36K with a mean temperature of 20{\pm}3K. Clump temperature is strongly correlated with proximity to the external OB association and there is no evidence that the embedded protostars significantly heat the dust. We identify 31 clumps that have cores with densities greater than 105cm{^{-3}}. Thirteen of these cores contain embedded Class 0/I protostars. Many cores are associated with bright-rimmed clouds seen in Herschel 70 {\mu}m images. From JCMT HARP observations of the 12CO 3-2 line, we find contamination of the 850{\mu}m band of up to 20 per cent. We investigate the free-free contribution to SCUBA-2 bands from large-scale and ultracompact H ii regions using archival VLA data and find the contribution is limited to individual stars, accounting for 9 per cent of flux per beam at 450 {\mu}m or 12 per cent at 850 {\mu}m in these cases. We conclude that radiative heating has potentially influenced the formation of stars in the Dust Arc sub-region, favouring Jeans stable clouds in the warm east and fragmentation in the cool west.

Debris disks, the tenuous rocky and icy remnants of planet formation, are believed to be evidence for planetary systems around other stars. The JCMT/SCUBA-2 debris disk legacy survey 'SCUBA-2 Observations of Nearby Stars' (SONS) observed 100 nearby stars, amongst them HD~76582, for evidence of such material. Here we present imaging observations by JCMT/SCUBA-2 and \textit{Herschel}/PACS at sub-millimetre and far-infrared wavelengths, respectively. We simultaneously model the ensemble of photometric and imaging data, spanning optical to sub-millimetre wavelengths, in a self-consistent manner. At far-infrared wavelengths, we find extended emission from the circumstellar disk providing a strong constraint on the dust spatial location in the outer system, although the angular resolution is too poor to constrain the interior of the system. In the sub-millimetre, photometry at 450 and 850~$\mu$m reveal a steep fall-off that we interpret as a disk dominated by moderately-sized dust grains ($a_{\rm min}~=~36~\mu$m), perhaps indicative of a non-steady-state collisional cascade within the disk. A disk architecture of three distinct annuli, comprising an unresolved component at $\sim$ 20 au and outer components at 80 and 270 au, along with a very steep particle size distribution ($\gamma~=~5$), is proposed to match the observations.

The young star beta Pictoris is well known for its dusty debris disk, produced through the grinding down by collisions of planetesimals, kilometre-sized bodies in orbit around the star. In addition to dust, small amounts of gas are also known to orbit the star, likely the result from vaporisation of violently colliding dust grains. The disk is seen edge on and from previous absorption spectroscopy we know that the gas is very rich in carbon relative to other elements. The oxygen content has been more difficult to assess, however, with early estimates finding very little oxygen in the gas at a C/O ratio 20x higher than the cosmic value. A C/O ratio that high is difficult to explain and would have far-reaching consequences for planet formation. Here we report on observations by the far-infrared space telescope Herschel, using PACS, of emission lines from ionised carbon and neutral oxygen. The detected emission from C+ is consistent with that previously reported being observed by the HIFI instrument on Herschel, while the emission from O is hard to explain without assuming a higher-density region in the disk, perhaps in the shape of a clump or a dense torus, required to sufficiently excite the O atoms. A possible scenario is that the C/O gas is produced by the same process responsible for the CO clump recently observed by ALMA in the disk, and that the re-distribution of the gas takes longer than previously assumed. A more detailed estimate of the C/O ratio and the mass of O will have to await better constraints on the C/O gas spatial distribution.

To characterize the mechanisms of planet formation it is crucial to investigate the properties and evolution of protoplanetary disks around young stars, where the initial conditions for the growth of planets are set. Our goal is to study grain growth in the disk of the young, intermediate mass star HD163296 where dust processing has already been observed, and to look for evidence of growth by ice condensation across the CO snowline, already identified in this disk with ALMA. Under the hypothesis of optically thin emission we compare images at different wavelengths from ALMA and VLA to measure the opacity spectral index across the disk and thus the maximum grain size. We also use a Bayesian tool based on a two-layer disk model to fit the observations and constrain the dust surface density. The measurements of the opacity spectral index indicate the presence of large grains and pebbles ($\geq$1 cm) in the inner regions of the disk (inside $\sim$50 AU) and smaller grains, consistent with ISM sizes, in the outer disk (beyond 150 AU). Re-analysing ALMA Band 7 Science Verification data we find (radially) unresolved excess continuum emission centered near the location of the CO snowline at $\sim$90 AU. Our analysis suggests a grain size distribution consistent with an enhanced production of large grains at the CO snowline and consequent transport to the inner regions. Our results combined with the excess in infrared scattered light found by Garufi et al. (2014) suggests the presence of a structure at 90~AU involving the whole vertical extent of the disk. This could be evidence for small scale processing of dust at the CO snowline.

Theoretical models of grain growth predict dust properties to change as a function of protoplanetary disk radius, mass, age and other physical conditions. We lay down the methodology for a multi-wavelength analysis of (sub-)mm and cm continuum interferometric observations to constrain self-consistently the disk structure and the radial variation of the dust properties. The computational architecture is massively parallel and highly modular. The analysis is based on the simultaneous fit in the uv-plane of observations at several wavelengths with a model for the disk thermal emission and for the dust opacity. The observed flux density at the different wavelengths is fitted by posing constraints on the disk structure and on the radial variation of the grain size distribution. We apply the analysis to observations of three protoplanetary disks (AS 209, FT Tau, DR Tau) for which a combination of spatially resolved observations in the range ~0.88mm to ~10mm is available (from SMA, CARMA, and VLA), finding evidence of a decreasing maximum dust grain size (a_max) with radius. We derive large a_max values up to 1 cm in the inner disk between 15 and 30 AU and smaller grains with a_max~1 mm in the outer disk (R > 80AU). In this paper we develop a multi-wavelength analysis that will allow this missing quantity to be constrained for statistically relevant samples of disks and to investigate possible correlations with disk or stellar parameters.

We present a first look at the SCUBA-2 observations of three sub-regions of the Orion B molecular cloud: LDN 1622, NGC 2023/2024, and NGC 2068/2071, from the JCMT Gould Belt Legacy Survey. We identify 29, 564, and 322 dense cores in L1622, NGC 2023/2024, and NGC 2068/2071 respectively, using the SCUBA-2 850 micron map, and present their basic properties, including their peak fluxes, total fluxes, and sizes, and an estimate of the corresponding 450 micron peak fluxes and total fluxes, using the FellWalker source extraction algorithm. Assuming a constant temperature of 20 K, the starless dense cores have a mass function similar to that found in previous dense core analyses, with a Salpeter-like slope at the high-mass end. The majority of cores appear stable to gravitational collapse when considering only thermal pressure; indeed, most of the cores which have masses above the thermal Jeans mass are already associated with at least one protostar. At higher cloud column densities, above 1-2 x 10^23 cm^-2, most of the mass is found within dense cores, while at lower cloud column densities, below 1 x 10^23 cm^-2, this fraction drops to 10% or lower. Overall, the fraction of dense cores associated with a protostar is quite small (<8%), but becomes larger for the densest and most centrally concentrated cores. NGC 2023 / 2024 and NGC 2068/2071 appear to be on the path to forming a significant number of stars in the future, while L1622 has little additional mass in dense cores to form many new stars.

We propose a set of standard assumptions for the modelling of Class II and III protoplanetary disks, which includes detailed continuum radiative transfer, thermo-chemical modelling of gas and ice, and line radiative transfer from optical to cm wavelengths. We propose new standard dust opacities for disk models, we present a simplified treatment of PAHs sufficient to reproduce the PAH emission features, and we suggest using a simple treatment of dust settling. We roughly adjust parameters to obtain a model that predicts typical Class II T Tauri star continuum and line observations. We systematically study the impact of each model parameter (disk mass, disk extension and shape, dust settling, dust size and opacity, gas/dust ratio, etc.) on all continuum and line observables, in particular on the SED, mm-slope, continuum visibilities, and emission lines including [OI] 63um, high-J CO lines, (sub-)mm CO isotopologue lines, and CO fundamental ro-vibrational lines. We find that evolved dust properties (large grains) often needed to fit the SED, have important consequences for disk chemistry and heating/cooling balance, leading to stronger emission lines in general. Strong dust settling and missing disk flaring have similar effects on continuum observations, but opposite effects on far-IR gas emission lines. PAH molecules can shield the gas from stellar UV radiation because of their strong absorption and negligible scattering opacities. The observable millimetre-slope of the SED can become significantly more gentle in the case of cold disk midplanes, which we find regularly in our T Tauri models. We propose to use line observations of robust chemical tracers of the gas, such as O, CO, and H2, as additional constraints to determine some key properties of the disks, such as disk shape and mass, opacities, and the dust/gas ratio, by simultaneously fitting continuum and line observations.

We present new results from the Disks@EVLA program for two young stars: CY Tau and DoAr 25. We trace continuum emission arising from their circusmtellar disks from spatially resolved observations, down to tens of AU scales, at {\lambda} = 0.9, 2.8, 8.0, and 9.8 mm for DoAr25 and at {\lambda} = 1.3, 2.8, and 7.1 mm for CY Tau. Additionally, we constrain the amount of emission whose origin is different from thermal dust emission from 5 cm observations. Directly from interferometric data, we find that observations at 7 mm and 1 cm trace emission from a compact disk while millimeter-wave observations trace an extended disk structure. From a physical disk model, where we characterize the disk structure of CY Tau and DoAr 25 at wavelengths shorter than 5 cm, we find that (1) dust continuum emission is optically thin at the observed wavelengths and over the spatial scales studied, (2) a constant value of the dust opacity is not warranted by our observations, and (3) a high-significance radial gradient of the dust opacity spectral index, {\beta}, is consistent with the observed dust emission in both disks, with low-{\beta} in the inner disk and high-{\beta} in the outer disk. Assuming that changes in dust properties arise solely due to changes in the maximum particle size (amax), we constrain radial variations of amax in both disks, from cm-sized particles in the inner disk (R < 40 AU) to millimeter sizes in the outer disk (R > 80 AU). These observational constraints agree with theoretical predictions of the radial-drift barrier, however, fragmentation of dust grains could explain our amax(R) constraints if these disks have lower turbulence and/or if dust can survive high-velocity collisions.

We present far-infrared and submillimeter maps from the Herschel Space Observatory and the James Clerk Maxwell Telescope of the debris disk host star AU Microscopii. Disk emission is detected at 70, 160, 250, 350, 450, 500 and 850 micron. The disk is resolved at 70, 160 and 450 micron. In addition to the planetesimal belt, we detect thermal emission from AU Mic's halo for the first time. In contrast to the scattered light images, no asymmetries are evident in the disk. The fractional luminosity of the disk is $3.9 \times 10^{-4}$ and its mm-grain dust mass is 0.01 MEarth (+/- 20%). We create a simple spatial model that reconciles the disk SED as a blackbody of 53 +/- 2 K (a composite of 39 and 50 K components) and the presence of small (non-blackbody) grains which populate the extended halo. The best fit model is consistent with the "birth ring" model explored in earlier works, i.e., an edge-on dust belt extending from 8.8-40 AU, but with an additional halo component with an $r^{-1.5}$ surface density profile extending to the limits of sensitivity (140 AU). We confirm that AU Mic does not exert enough radiation force to blow out grains. For stellar mass loss rates of 10-100x solar, compact (zero porosity) grains can only be removed if they are very small, consistently with previous work, if the porosity is 0.9, then grains approaching 0.1 micron can be removed via corpuscular forces (i.e., the stellar wind).

We present new observations of the Kuiper belt analogues around HD 38858 and HD 20794, hosts of super-Earth mass planets within 1 au. As two of the four nearby G-type stars (with HD 69830 and 61 Vir) that form the basis of a possible correlation between low-mass planets and debris disc brightness, these systems are of particular interest. The disc around HD 38858 is well resolved with Herschel and we constrain the disc geometry and radial structure. We also present a probable JCMT sub-mm continuum detection of the disc and a CO J=2-1 upper limit. The disc around HD 20794 is much fainter and appears marginally resolved with Herschel, and is constrained to be less extended than the discs around 61 Vir and HD 38858. We also set limits on the radial location of hot dust recently detected around HD 20794 with near-IR interferometry. We present HARPS upper limits on unseen planets in these four systems, ruling out additional super-Earths within a few au, and Saturn-mass planets within 10 au. We consider the disc structure in the three systems with Kuiper belt analogues (HD 69830 has only a warm dust detection), concluding that 61 Vir and HD 38858 have greater radial disc extent than HD 20794. We speculate that the greater width is related to the greater minimum planet masses (10-20 $M_\oplus$ vs. 3-5 $M_\oplus$), arising from an eccentric planetesimal population analogous to the Solar System's scattered disc. We discuss alternative scenarios and possible means to distinguish among them.

Circumstellar dust disks have been observed around many nearby stars. However, many stars are part of binary or multiple stellar systems. A natural question arises regarding the presence and properties of such disks in systems with more than one star. To address this, we consider a sample of 449 systems (spectral types A-M) observed with the Herschel Space Observatory as part of the DEBRIS program. We have examined the stellar multiplicity of this sample by gathering information from the literature and performing an adaptive optics imaging survey at Lick Observatory. Five new companions were revealed with our program. In total, we identify 188 (42%) binary or multiple star systems. The multiplicity of the sample is examined with regards to the detection of circumstellar disks for stars of spectral types AFGK. In general, disks are less commonly detected around binaries than single stars, though the disk frequency is comparable among A stars regardless of multiplicity. However, this sample reveals the period distribution of disk-bearing binaries is consistent with that of non-disk binaries and with comparison field samples. We find that the properties of disks in binary systems are not statistically different from those around single stars. Although the frequency of disk-bearing FGK binaries may be lower than in single star systems, the processes behind disk formation and the characteristics of these disks are comparable among both populations.

We present 850$\mu$m and 450$\mu$m data from the JCMT Gould Belt Survey obtained with SCUBA-2 and characterise the dust attributes of Class I, Class II and Class III disk sources in L1495. We detect 23% of the sample at both wavelengths, with the detection rate decreasing through the Classes from I--III. The median disk mask is 1.6$\times 10^{-3}$M$_{\odot}$, and only 7% of Class II sources have disk masses larger than 20 Jupiter masses. We detect a higher proportion of disks towards sources with stellar hosts of spectral type K than spectral type M. Class II disks with single stellar hosts of spectral type K have higher masses than those of spectral type M, supporting the hypothesis that higher mass stars have more massive disks. Variations in disk masses calculated at the two wavelengths suggests there may be differences in dust opacity and/or dust temperature between disks with hosts of spectral types K to those with spectral type M.

The study of the planet-debris disk connection can shed light on the formation and evolution of planetary systems, and may help predict the presence of planets around stars with certain disk characteristics. In preliminary analyses of the Herschel DEBRIS and DUNES surveys, Wyatt et al. (2012) and Marshall et al. (2014) identified a tentative correlation between debris and low-mass planets. Here we use the cleanest possible sample out these surveys to assess the presence of such a correlation, discarding stars without known ages, with ages < 1 Gyr and with binary companions <100 AU, to rule out possible correlations due to effects other than planet presence. In our sample of 204 FGK stars, we do not find evidence that debris disks are more common or more dusty around stars harboring high-mass or low-mass planets compared to a control sample without identified planets, nor that debris disks are more or less common (or more or less dusty) around stars harboring multiple planets compared to single-planet systems. Diverse dynamical histories may account for the lack of correlations. The data show the correlation between the presence of high-mass planets and stellar metallicity, but no correlation between the presence of low-mass planets or debris and stellar metallicity. Comparing the observed cumulative distribution of fractional luminosity to those expected from a Gaussian distribution, we find that a distribution centered on the Solar system's value fits well the data, while one centered at 10 times this value can be rejected. This is of interest in the context of future terrestrial planet characterization because it indicates that there are good prospects for finding a large number of debris disk systems (i.e. with evidence of harboring the building blocks of planets) with exozodiacal emission low enough to be appropriate targets for an ATLAST-type mission to search for biosignatures.

In this paper we present the first observations of the Ophiuchus molecular cloud performed as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) with the SCUBA-2 instrument. We demonstrate methods for combining these data with previous HARP CO, Herschel, and IRAM N$_{2}$H$^{+}$ observations in order to accurately quantify the properties of the SCUBA-2 sources in Ophiuchus. We produce a catalogue of all of the sources found by SCUBA-2. We separate these into protostars and starless cores. We list all of the starless cores and perform a full virial analysis, including external pressure. This is the first time that external pressure has been included in this level of detail. We find that the majority of our cores are either bound or virialised. Gravitational energy and external pressure are on average of a similar order of magnitude, but with some variation from region to region. We find that cores in the Oph A region are gravitationally bound prestellar cores, while cores in the Oph C and E regions are pressure-confined. We determine that N$_{2}$H$^{+}$ is a good tracer of the bound material of prestellar cores, although we find some evidence for N$_{2}$H$^{+}$ freeze-out at the very highest core densities. We find that non-thermal linewidths decrease substantially between the gas traced by C$^{18}$O and that traced by N$_{2}$H$^{+}$, indicating the dissipation of turbulence at higher densities. We find that the critical Bonnor-Ebert stability criterion is not a good indicator of the boundedness of our cores. We detect the pre-brown dwarf candidate Oph B-11 and find a flux density and mass consistent with previous work. We discuss regional variations in the nature of the cores and find further support for our previous hypothesis of a global evolutionary gradient across the cloud from southwest to northeast, indicating sequential star formation across the region.

Earth will become uninhabitable within 2-3 Gyr as a result of the moving boundaries of the habitable zone caused by the increasing luminosity of the Sun. Predictions about the future of habitable conditions on Earth include a decline in species diversity and habitat extent, ocean loss and changes in the magnitudes of geochemical cycles. However, testing these predictions on the present-day Earth is difficult. The discovery of a planet that is a near analogue to the far future Earth could provide a means to test these predictions. Such a planet would need to have an Earth-like biosphere history, requiring it to have been in its system's habitable zone (HZ) for Gyr-long periods during the system's past, and to be approaching the inner-edge of the HZ at present. Here we assess the possibility of finding this very specific type of exoplanet and discuss the benefits of analysing older Earths in terms of improving our understanding of long-term geological and bio-geological processes. As an illustrative example, G stars within 10 parsecs are assessed as potential old-Earth-analogue hosts. Surface temperature estimates for hypothetical inner-HZ Earth analogues are used to determine whether any such planets in these systems would be at the right stage in their late-habitable lifetimes to exhibit detectable biosignatures. Predictions from planet formation studies and biosphere evolution models suggest that only 0.36% of G stars in the solar neighbourhood could host an old-Earth-analogue. However, if the development of an Earth-like biosphere is assumed to be rare, requiring a sequence of low-probability events to occur, then such planets are unlikely to be found in the solar neighbourhood - although 1000s could be present in the galaxy as a whole.

We present SCUBA-2 450micron and 850micron observations of the Serpens MWC 297 region, part of the JCMT Gould Belt Survey of nearby star-forming regions. Simulations suggest that radiative feedback influences the star-formation process and we investigate observational evidence for this by constructing temperature maps. Maps are derived from the ratio of SCUBA-2 fluxes and a two component model of the JCMT beam for a fixed dust opacity spectral index of beta = 1.8. Within 40 of the B1.5Ve Herbig star MWC 297, the submillimetre fluxes are contaminated by free-free emission with a spectral index of 1.03+-0.02, consistent with an ultra-compact HII region and polar winds/jets. Contamination accounts for 73+-5 per cent and 82+-4 per cent of peak flux at 450micron and 850micron respectively. The residual thermal disk of the star is almost undetectable at these wavelengths. Young Stellar Objects are confirmed where SCUBA-2 850micron clumps identified by the fellwalker algorithm coincide with Spitzer Gould Belt Survey detections. We identify 23 objects and use Tbol to classify nine YSOs with masses 0.09 to 5.1 Msun. We find two Class 0, one Class 0/I, three Class I and three Class II sources. The mean temperature is 15+-2K for the nine YSOs and 32+-4K for the 14 starless clumps. We observe a starless clump with an abnormally high mean temperature of 46+-2K and conclude that it is radiatively heated by the star MWC 297. Jeans stability provides evidence that radiative heating by the star MWC 297 may be suppressing clump collapse.

NGC1333 is a 1-2 Myr old cluster of stars in the Perseus molecular cloud. We used 850mu data from the Gould Belt Survey with SCUBA-2 on the JCMT to measure or place limits on disc masses for 82 Class II sources in this cluster. Eight disc-candidates were detected; one is estimated to have mass of about 9 Jupiter masses in dust plus gas, while the others host only 2-4 Jupiter masses of circumstellar material. None of these discs exceeds the threshold for the 'Minimum Mass Solar Nebula' (MMSN). This reinforces previous claims that only a small fraction of Class II sources at an age of 1-2 Myr has discs exceeding the MMSN threshold and thus can form a planetary system like our own. However, other regions with similarly low fractions of MMSN discs (IC348, UpSco, SigmaOri) are thought to be older than NGC1333. Compared with coeval regions, the exceptionally low fraction of massive discs in NGC1333 cannot easily be explained by the effects of UV radiation or stellar encounters. Our results indicate that additional environmental factors significantly affect disc evolution and the outcome of planet formation by core accretion.

Molecular gas constitutes the dominant mass component of protoplanetary discs. To date, these sources have not been studied comprehensively at the longest far-infrared and shortest submillimetre wavelengths. This paper presents Herschel SPIRE FTS spectroscopic observations toward 18 protoplanetary discs, covering the entire 450-1540 GHz (666-195 $\mu$m) range at R~400-1300. The spectra reveal clear detections of the dust continuum and, in six targets, a significant amount of spectral line emission primarily attributable to $^{12}$CO rotational lines. Other targets exhibit little to no detectable spectral lines. Low signal-to-noise detections also include signatures from $^{13}$CO, [CI] and HCN. For completeness, we present upper limits of non-detected lines in all targets, including low-energy transitions of H2O and CH$^+$ molecules. The ten $^{12}$CO lines that fall within the SPIRE FTS bands trace energy levels of ~50-500 K. Combined with lower and higher energy lines from the literature, we compare the CO rotational line energy distribution with detailed physical-chemical models, for sources where these are available and published. Our 13CO line detections in the disc around Herbig Be star HD 100546 exceed, by factors of ~10-30, the values predicted by a model that matches a wealth of other observational constraints, including the SPIRE $^{12}$CO ladder. To explain the observed $^{12}$CO/$^{13}$CO ratio, it may be necessary to consider the combined effects of optical depth and isotope selective (photo)chemical processes. Considering the full sample of 18 objects, we find that the strongest line emission is observed in discs around Herbig Ae/Be stars, although not all show line emission. In addition, two of the six T Tauri objects exhibit detectable $^{12}$CO lines in the SPIRE range.

$\tau$ Ceti is a nearby, mature G-type star very similar to our Sun, with a massive Kuiper Belt analogue (Greaves et al. 2004) and possible multiplanet system (Tuomi et al. 2013) that has been compared to our Solar System. We present Herschel Space Observatory images of the debris disk, finding the disk is resolved at 70 and 160 microns, and marginally resolved at 250 microns. The Herschel images and infrared photometry from the literature are best modelled using a wide dust annulus with an inner edge between 1-10 AU and an outer edge at ~55 AU, inclined from face-on by 35$\pm$10 degrees, and with no significant azimuthal structure. We model the proposed tightly-packed planetary system of five super-Earths and find that the innermost dynamically stable disk orbits are consistent with the inner edge found by the observations. The photometric modelling, however, cannot rule out a disk inner edge as close to the star as 1 AU, though larger distances produce a better fit to the data. Dynamical modelling shows that the 5 planet system is stable with the addition of a Neptune or smaller mass planet on an orbit outside 5 AU, where the Tuomi et al. analysis would not have detected a planet of this mass.

ALMA Cycle 0 and Herschel PACS observations are reported for the prototype, nearest, and brightest example of a dusty and polluted white dwarf, G29-38. These long wavelength programs attempted to detect an outlying, parent population of bodies at 1-100 AU, from which originates the disrupted planetesimal debris that is observed within 0.01 AU and which exhibits L_IR/L = 0.039. No associated emission sources were detected in any of the data down to L_IR/L ~ 1e-4, generally ruling out cold dust masses greater than 1e24 - 1e25 g for reasonable grain sizes and properties in orbital regions corresponding to evolved versions of both asteroid and Kuiper belt analogs. Overall, these null detections are consistent with models of long-term collisional evolution in planetesimal disks, and the source regions for the disrupted parent bodies at stars like G29-38 may only be salient in exceptional circumstances, such as a recent instability. A larger sample of polluted white dwarfs, targeted with the full ALMA array, has the potential to unambiguously identify the parent source(s) of their planetary debris.

In light of recent substantial updates to spectral type estimations and newly established intrinsic colours, effective temperatures, and bolometric corrections for pre-main sequence (PMS) stars, we re-address the theory of accretion-disc regulated stellar angular momentum (AM) evolution. We report on the compilation of a consistent sample of fully convective stars within two of the most well-studied and youngest, nearby regions of star formation: the Orion Nebula Cluster (ONC) and Taurus-Auriga. We calculate the average specific stellar AM ($j_{\star}$) assuming solid body rotation, using surface rotation periods gathered from the literature and new estimates of stellar radii and ages. We use published Spitzer IRAC fluxes to classify our stars as Class II or Class III and compare their $j_{\star}$ evolution. Our results suggest that disc dispersal is a rapid process that occurs at a variety of ages. We find a consistent $j_{\star}$ reduction rate between the Class II and Class III PMS stars which we interpret as indicating a period of accretion disc-regulated AM evolution followed by near-constant AM evolution once the disc has dissipated. Furthermore, assuming our observed spread in stellar ages is real, we find the removal rate of $j_{\star}$ during the Class II phase is more rapid than expected by contraction at constant stellar rotation rate. A much more efficient process of AM removal must exist, most likely in the form of an accretion-driven stellar wind or other outflow from the star-disc interaction region or extended disc surface.

We present sub-arcsecond resolution observations of continuum emission associated with the GG Tau quadruple star system at wavelengths of 1.3, 2.8, 7.3, and 50 mm. These data confirm that the GG Tau A binary is encircled by a circumbinary ring at a radius of 235 AU with a FWHM width of ~60 AU. We find no clear evidence for a radial gradient in the spectral shape of the ring, suggesting that the particle size distribution is spatially homogeneous on angular scales of ~0.1". A central point source, likely associated with the primary component (GG Tau Aa), exhibits a composite spectrum from dust and free-free emission. Faint emission at 7.3 mm is observed toward the low-mass star GG Tau Ba, although its origin remains uncertain. Using these measurements of the resolved, multifrequency emission structure of the GG Tau A system, models of the far-infrared to radio spectrum are developed to place constraints on the grain size distribution and dust mass in the circumbinary ring. The non-negligible curvature present in the ring spectrum implies a maximum particle size of 1-10 mm, although we are unable to place strong constraints on the distribution shape. The corresponding dust mass is 30-300 M_earth, at a temperature of 20-30 K. We discuss how this significant concentration of relatively large particles in a narrow ring at a large radius might be produced in a local region of higher gas pressures (i.e., a particle "trap") located near the inner edge of the circumbinary disk.

Many stars are surrounded by disks of dusty debris formed in the collisions of asteroids, comets and dwarf planets. But is gas also released in such events? Observations at submm wavelengths of the archetypal debris disk around $\beta$ Pictoris show that 0.3% of a Moon mass of carbon monoxide orbits in its debris belt. The gas distribution is highly asymmetric, with 30% found in a single clump 85AU from the star, in a plane closely aligned with the orbit of the inner planet, $\beta$ Pic b. This gas clump delineates a region of enhanced collisions, either from a mean motion resonance with an unseen giant planet, or from the remnants of a collision of Mars-mass planets.

The $\textit{Herschel}$ DEBRIS, DUNES and GT programmes observed 37 exoplanet host stars within 25 pc at 70, 100 and 160 $\mu$m with the sensitivity to detect far-infrared excess emission at flux density levels only an order of magnitude greater than that of the Solar system's Edgeworth-Kuiper belt. Here we present an analysis of that sample, using it to more accurately determine the (possible) level of dust emission from these exoplanet host stars and thereafter determine the links between the various components of these exoplanetary systems through statistical analysis. We have fitted the flux densities measured from recent \textit{Herschel} observations with a simple two parameter ($T_{d}$, $L_{\rm IR}/L_{\star}$) black body model (or to the 3-$\sigma$ upper limits at 100 $\mu$m). From this uniform approach we calculate the fractional luminosity, radial extent, dust temperature and disc mass. We then plotted the calculated dust luminosity or upper limits against the stellar properties, e.g. effective temperature, metallicity, age, and identified correlations between these parameters. A total of eleven debris discs are identified around the 37 stars in the sample. An incidence of ten cool debris discs around the Sun-like exoplanet host stars (29 $\pm$ 9 %) is consistent with the detection rate found by DUNES (20.2 $\pm$ 2.0 %). For the debris disc systems, the dust temperatures range from 20 to 80 K, and fractional luminosities ($L_{\rm IR}/L_{\star}$) between 2.4 $\times$10$^{-6}$ and 4.1 $\times$10$^{-4}$. In the case of non-detections, we calculated typical 3-$\sigma$ upper limits to the dust fractional luminosities of a few $\times10^{-6}$. We recover the previously identified correlation between stellar metallicity and hot Jupiter planets in our data set. We find a correlation between the increased presence of dust, lower planet masses and lower stellar metallicities. (abridged)

For over a decade, the structure of the inner cavity in the transition disk of TW Hydrae has been a subject of debate. Modeling the disk with data obtained at different wavelengths has led to a variety of proposed disk structures. Rather than being inconsistent, the individual models might point to the different faces of physical processes going on in disks, such as dust growth and planet formation. Our aim is to investigate the structure of the transition disk again and to find to what extent we can reconcile apparent model differences. A large set of high-angular-resolution data was collected from near-infrared to centimeter wavelengths. We investigated the existing disk models and established a new self-consistent radiative-transfer model. A genetic fitting algorithm was used to automatize the parameter fitting. Simple disk models with a vertical inner rim and a radially homogeneous dust composition from small to large grains cannot reproduce the combined data set. Two modifications are applied to this simple disk model: (1) the inner rim is smoothed by exponentially decreasing the surface density in the inner ~3 AU, and (2) the largest grains (>100 um) are concentrated towards the inner disk region. Both properties can be linked to fundamental processes that determine the evolution of protoplanetary disks: the shaping by a possible companion and the different regimes of dust-grain growth, respectively. The full interferometric data set from near-infrared to centimeter wavelengths requires a revision of existing models for the TW Hya disk. We present a new model that incorporates the characteristic structures of previous models but deviates in two key aspects: it does not have a sharp edge at 4 AU, and the surface density of large grains differs from that of smaller grains. This is the first successful radiative-transfer-based model for a full set of interferometric data.

We present new observations of the XZ Tau system made at high angular resolution (55 milliarcsec) with the Karl G. Jansky Very Large Array (VLA) at a wavelength of 7 mm. Observations of XZ Tau made with the VLA in 2004 appeared to show a triple system, with XZ Tau A resolved into two sources, XZ Tau A and XZ Tau C. The angular separation of XZ Tau A and C (0.09 arcsec) suggested a projected orbital separation of around 13 AU with a possible orbital period of around 40 yr. Our observations were obtained approximately 8 yr later, a fifth of this putative orbital period, and should therefore allow us to constrain the orbit of XZ Tau C, and evaluate the possibility that a recent periastron passage of C coincided with the launch of extended optical outflows from XZ Tau A. Despite improved sensitivity and resolution compared with previous observations, we find no evidence of XZ Tau C in our data. Components A and B are detected with a signal-to-noise ratio greater than ten; their orbits are consistent with previous studies, although the emission from XZ Tau A appears to be weaker. Three possible interpretations are offered: either XZ Tau C is transiting XZ Tau A, which is broadly consistent with the periastron passage hypothesis, or the emission seen in 2004 was that of a transient, or XZ Tau C does not exist. A fourth interpretation, that XZ Tau C was ejected from the system, is dismissed due to the lack of angular momentum redistribution in the orbits of XZ Tau A and XZ Tau B that would result from such an event. Our observations are insufficient to distinguish between the remaining possibilities, at least not until we obtain further VLA observations at a sufficiently later time. A further non-detection would allow us to reject the transit hypothesis, and the periastron passage of XZ Tau C as agent of XZ Tau A's outflows.

We present far-infrared and sub-millimeter images of the eta Crv debris disk system obtained with Herschel and SCUBA-2, as well as Hubble Space Telescope visible and near-infrared coronagraphic images. In the 70 micron Herschel image, we clearly separate the thermal emission from the warm and cold belts in the system, find no evidence for a putative dust population located between them, and precisely determine the geometry of the outer belt. We also find marginal evidence for azimuthal asymmetries and a global offset of the outer debris ring relative to the central star. Finally, we place stringent upper limits on the scattered light surface brightness of the outer ring. Using radiative transfer modeling, we find that it is impossible to account for all observed properties of the system under the assumption that both rings contain dust populations with the same properties. While the outer belt is in reasonable agreement with the expectations of steady-state collisional cascade models, albeit with a minimum grain size that is four times larger than the blow-out size, the inner belt appears to contain copious amounts of small dust grains, possibly below the blow-out size. This suggests that the inner belt cannot result from a simple transport of grains from the outer belt and rather supports a more violent phenomenon as its origin. We also find that the emission from the inner belt has not declined over three decades, a much longer timescale than its dynamical timescale, which indicates that the belt is efficiently replenished.

We mapped five massive star forming regions with the SCUBA-2 camera on the James Clerk Maxwell Telescope (JCMT). Temperature and column density maps are obtained from the SCUBA-2 450 and 850 $\mu$m images. Most of the dense clumps we find have central temperatures below 20 K with some as cold as 8 K, suggesting that they have no internal heating due to the presence of embedded protostars. This is surprising, because at the high densities inferred from these images and at these low temperatures such clumps should be unstable, collapsing to form stars and generating internal heating. The column densities at the clump centres exceed 10$^{23}$ cm$^{-2}$, and the derived peak visual extinction values are from 25-500 mag for $\beta$ = 1.5-2.5, indicating highly opaque centres. The observed cloud gas masses range from $\sim$ 10 to 10$^{3}$ M$_{\odot}$. The outer regions of the clumps follow an $r^{-2.36\pm0.35}$ density distribution and this power-law structure is observed outside of typically 10$^{4}$ AU. All these findings suggest that these clumps are high-mass starless clumps and most likely contain high-mass starless cores.

Fomalhaut is one of the most interesting and well studied nearby stars, hosting at least one planet, a spectacular debris ring, and two distant low-mass stellar companions (TW PsA and LP 876-10, a.k.a. Fomalhaut B & C). We observed both companions with Herschel, and while no disc was detected around the secondary, TW PsA, we have discovered the second debris disc in the Fomalhaut system, around LP 876-10. This detection is only the second case of two debris discs seen in a multiple system, both of which are relatively wide ($\gtrsim$3000 AU for HD 223352/40 and 158 kAU [0.77 pc] for Fomalhaut/LP 876-10). The disc is cool (24K) and relatively bright, with a fractional luminosity $L_{\rm disc}/L_\star = 1.2 \times 10^{-4}$, and represents the rare observation of a debris disc around an M dwarf. Further work should attempt to find if the presence of two discs in the Fomalhaut system is coincidental, perhaps simply due to the relatively young system age of 440 Myr, or if the stellar components have dynamically interacted and the system is even more complex than it currently appears.

Many nearby main-sequence stars have been searched for debris using the far-infrared Herschel satellite, within the DEBRIS, DUNES and Guaranteed-Time Key Projects. We discuss here 11 stars of spectral types A to M where the stellar inclination is known and can be compared to that of the spatially-resolved dust belts. The discs are found to be well aligned with the stellar equators, as in the case of the Sun's Kuiper belt, and unlike many close-in planets seen in transit surveys. The ensemble of stars here can be fitted with a star-disc tilt of ~<10 degrees. These results suggest that proposed mechanisms for tilting the star or disc in fact operate rarely. A few systems also host imaged planets, whose orbits at tens of AU are aligned with the debris discs, contrary to what might be expected in models where external perturbers induce tilts.

The biosignatures of life on Earth do not remain static, but change considerably over the planet's habitable lifetime. Earth's future biosphere, much like that of the early Earth, will consist of predominantly unicellular microorganisms due to the increased hostility of environmental conditions caused by the Sun as it enters the late stage of its main sequence evolution. Building on previous work, the productivity of the biosphere is evaluated during different stages of biosphere decline between 1 Gyr and 2.8 Gyr from present. A simple atmosphere-biosphere interaction model is used to estimate the atmospheric biomarker gas abundances at each stage and to assess the likelihood of remotely detecting the presence of life in low-productivity, microbial biospheres, putting an upper limit on the lifetime of Earth's remotely detectable biosignatures. Other potential biosignatures such as leaf reflectance and cloud cover are discussed.

We present very high spatial resolution deep radio continuum observations at 5 GHz (6 cm) made with e-MERLIN of the young stars DG Tau A and B. Assuming it is launched very close (~=1 au) from the star, our results suggest that the DG Tau A outflow initially starts as a poorly focused wind and undergoes significant collimation further along the jet (~=50 au). We derive jet parameters for DG Tau A and find an initial jet opening angle of 86 degrees within 2 au of the source, a mass-loss rate of 1.5x10^-8 solar masses/yr for the ionised component of the jet, and the total ejection/accretion ratio to range from 0.06-0.3. These results are in line with predictions from MHD jet-launching theories.

We present SCUBA-2 850um observations for 7 very low mass stars (VLMS) and brown dwarfs (BDs): 3 in Taurus, 4 in the TWA, and all classical T Tauri (cTT) analogs. We detect 2 of the 3 Taurus disks, but none of the TWA ones. Our 3sigma limits correspond to a dust mass of 1.2 MEarth in Taurus and a mere 0.2 MEarth in the TWA (3--10x deeper than previous work). We combine our data with other sub-mm/mm surveys of Taurus, rho Oph and the TWA to investigate trends in disk mass and grain growth during the cTT phase. We find : (1) The minimum disk outer radius required to explain the upper envelope of sub-mm/mm fluxes is 100 AU for intermediate-mass stars, solar-types and VLMS, and 20 AU for BDs. (2) While the upper envelope of disk masses increases with Mstar from BDs to VLMS to solar-types, no increase is seen from solar-type to intermediate-mass stars. We propose this is due to enhanced photoevaporation around intermediate masses. (3) Many disks around Taurus and rho Oph intermediate-mass and solar-type stars evince an opacity index beta of 0--1, indicating large grains. Of the only four VLMS/BDs in these regions with multi-wavelength data, three are consistent with large grains, though optically thick disks are not ruled out. (4) For the TWA VLMS (TWA 30A, B), combining our fluxes with Mdot and ages suggests substantial grain growth by 10 Myr. The degree of grain growth in the TWA BDs (2M1207A, SSPM1102) remains largely unknown. (5) A Bayesian analysis shows that mean(log[Mdisk/Mstar]) = -2.4, roughly constant all the way from intermediate-mass stars to VLMS/BDs, and (6) the disk mass in close solar-type Taurus binaries is significantly lower than in singles (by a factor of 10), while that in wide solar-type Taurus binaries is closer to that in singles (lower by a factor of 3). (7) We discuss the implications for planet formation, and for the dependence of Mdot on Mstar.

We present Herschel observations from the Herschel Gould Belt Survey and SCUBA-2 science verification observations from the JCMT Gould Belt Survey of the B1 clump in the Perseus molecular cloud. We determined the dust emissivity index using four different techniques to combine the Herschel PACS+SPIRE data at 160 - 500 microns with the SCUBA-2 data at 450 microns and 850 microns. Of our four techniques, we found the most robust method was to filter-out the large-scale emission in the Herschel bands to match the spatial scales recovered by the SCUBA-2 reduction pipeline. Using this method, we find beta ~ 2 towards the filament region and moderately dense material and lower beta values (beta > 1.6) towards the dense protostellar cores, possibly due to dust grain growth. We find that beta and temperature are more robust with the inclusion of the SCUBA-2 data, improving estimates from Herschel data alone by factors of ~ 2 for beta and by ~ 40% for temperature. Furthermore, we find core mass differences of < 30% compared to Herschel-only estimates with an adopted beta = 2, highlighting the necessity of long wavelength submillimeter data for deriving accurate masses of prestellar and protostellar cores.

We have carried out a survey of the NGC 2068 region in the Orion B molecular cloud using HARP on the JCMT, in the 13CO and C18O (J = 3-2) and H13CO+ (J = 4-3) lines. We used 13CO to map the outflows in the region, and matched them with previously defined SCUBA cores. We decomposed the C18O and H13CO+ into Gaussian clumps, finding 26 and 8 clumps respectively. The average deconvolved radii of these clumps is 6200 +/- 2000 AU and 3600 +/- 900 AU for C18O and H13CO+ respectively. We have also calculated virial and gas masses for these clumps, and hence determined how bound they are. We find that the C18O clumps are more bound than the H13CO+ clumps (average gas mass to virial mass ratio of 4.9 compared to 1.4). We measure clump internal velocity dispersions of 0.28 +/- 0.02 kms-1 and 0.27 +/- 0.04 kms-1 for C18O and H13CO+ respectively, although the H13CO+ values are heavily weighted by a majority of the clumps being protostellar, and hence having intrinsically greater linewidths. We suggest that the starless clumps correspond to local turbulence minima, and we find that our clumps are consistent with formation by gravoturbulent fragmentation. We also calculate inter-clump velocity dispersions of 0.39 +/- 0.05 kms-1 and 0.28 +/- 0.08 kms-1 for C18O and H13CO+ respectively. The velocity dispersions (both internal and external) for our clumps match results from numerical simulations of decaying turbulence in a molecular cloud. However, there is still insufficient evidence to conclusively determine the type of turbulence and timescale of star formation, due to the small size of our sample.

We present observations of the debris disk around gamma Doradus, an F1V star, from the Herschel Key Programme DEBRIS (Disc Emission via Bias-free Reconnaissance in the Infrared/Submillimetre). The disk is well-resolved at 70, 100 and 160 micron, resolved along its major axis at 250 micron, detected but not resolved at 350 micron, and confused with a background source at 500 micron. It is one of our best resolved targets and we find it to have a radially broad dust distribution. The modelling of the resolved images cannot distinguish between two configurations: an arrangement of a warm inner ring at several AU (best-fit 4 AU) and a cool outer belt extending from ~55 to 400 AU or an arrangement of two cool, narrow rings at ~70 AU and ~190 AU. This suggests that any configuration between these two is also possible. Both models have a total fractional luminosity of ~10^{-5} and are consistent with the disk being aligned with the stellar equator. The inner edge of either possible configuration suggests that the most likely region to find planets in this system would be within ~55 AU of the star. A transient event is not needed to explain the warm dust's fractional luminosity.

Debris disks have been found primarily around intermediate and solar mass stars (spectral types A-K) but rarely around low mass M-type stars. We have spatially resolved a debris disk around the remarkable M3-type star GJ581 hosting multiple planets using deep PACS images at 70, 100 and 160 microns as part of the DEBRIS Program on the Herschel Space Observatory. This is the second spatially resolved debris disk found around an M-type star, after the one surrounding the young star AU Mic (12 Myr). However, GJ 581 is much older (2-8 Gyr), and is X-ray quiet in the ROSAT data. We fit an axisymmetric model of the disk to the three PACS images and found that the best fit model is for a disk extending radially from 25+/-12 AU to more than 60 AU. Such a cold disk is reminiscent of the Kuiper Belt but it surrounds a low mass star (0.3 M_sol) and its fractional dust luminosity L_dust/L_* of \sim 10^-4 is much higher. The inclination limits of the disk found in our analysis make the masses of the planets small enough to ensure the long-term stability of the system according to some dynamical simulations. The disk is collisionally dominated down to submicron-sized grains and the dust cannot be expelled from the system by radiation or wind pressures because of the low luminosity and low X-ray luminosity of GJ581. We suggest that the correlation between low-mass planets and debris disks recently found for G-type stars also applies to M-type stars.Finally, the known planets, of low masses and orbiting within 0.3 AU from the star, cannot dynamically perturb the disk over the age of the star, suggesting that an additional planet exists at larger distance that is stirring the disk to replenish the dust.

Some planetary systems harbour debris disks containing planetesimals such as asteroids and comets. Collisions between such bodies produce small dust particles, the spectral features of which reveal their composition and, hence, that of their parent bodies. A measurement of the composition of olivine crystals has been done for the protoplanetary disk HD100546 and for olivine crystals in the warm inner parts of planetary systems. The latter compares well with the iron-rich olivine in asteroids (x<0.29). In the cold outskirts of the Beta Pictoris system, an analogue to the young Solar System, olivine crystals were detected but their composition remained undetermined, leaving unknown how the composition of the bulk of Solar System cometary olivine grains compares with that of extrasolar comets. Here we report the detection of the 69-micrometre-wavelength band of olivine crystals in the spectrum of Beta Pictoris. Because the disk is optically thin, we can associate the crystals with an extrasolar proto-Kuiper belt a distance of 15-45 astronomical units from the star (one astronomical unit is the Sun-Earth distance), determine their magnesium-rich composition (x=0.01+/-0.001) and show that they make up 3.6+/-1.0 per cent of the total dust mass. These values are strikingly similar to those for the dust emitted by the most primitive comets in the Solar System, even though Beta Pictoris is more massive and more luminous and has a different planetary system architecture.

We use data from the Disc Emission via a Bias-free Reconnaissance in the Infrared/Submillimetre (DEBRIS) survey, taken at 100 um with the Photoconductor Array Camera and Spectrometer instrument on board the Herschel Space Observatory, to make a cosmic variance independent measurement of the extragalactic number counts. These data consist of 323 small-area mapping observations performed uniformly across the sky, and thus represent a sparse sampling of the astronomical sky with an effective coverage of ~2.5 deg^2. We find our cosmic variance independent analysis to be consistent with previous count measurements made using relatively small area surveys. Furthermore, we find no statistically significant cosmic variance on any scale within the errors of our data. Finally, we interpret these results to estimate the probability of galaxy source confusion in the study of debris discs.

The future biosphere on Earth (as with its past) will be made up predominantly of unicellular microorganisms. Unicellular life was probably present for at least 2.5 Gyr before multicellular life appeared and will likely be the only form of life capable of surviving on the planet in the far future, when the ageing Sun causes environmental conditions to become more hostile to more complex forms of life. Therefore, it is statistically more likely that habitable Earth-like exoplanets we discover will be at a stage in their habitable lifetime more conducive to supporting unicellular, rather than multicellular life. The end stage of habitability on Earth is the focus of this work. A simple, latitude-based climate model incorporating eccentricity and obliquity variations is used as a guide to the temperature evolution of the Earth over the next 3 Gyr. This allows inferences to be made about potential refuges for life, particularly in mountains and cold-trap (ice) caves and what forms of life could live in these environments. Results suggest that in high latitude regions, unicellular life could persist for up to 2.8 Gyr from present. This begins to answer the question of how the habitability of Earth will evolve at local scales alongside the Sun's main sequence evolution and, by extension, how the habitability of Earth-like planets would evolve over time with their own host stars.

We present dust continuum observations of the protoplanetary disk surrounding the pre-main sequence star AS 209, spanning more than an order of magnitude in wavelength from 0.88 to 9.8 mm. The disk was observed with sub-arcsecond angular resolution (0.2"-0.5") to investigate radial variations in its dust properties. At longer wavelengths, the disk emission structure is notably more compact, providing model-independent evidence for changes in the grain properties across the disk. We find that physical models which reproduce the disk emission require a radial dependence of the dust opacity \kappa_{\nu}. Assuming that the observed wavelength-dependent structure can be attributed to radial variations in the dust opacity spectral index (\beta), we find that \beta(R) increases from \beta<0.5 at \sim20 AU to \beta>1.5 for R>80 AU, inconsistent with a constant value of \beta\ across the disk (at the 10\sigma\ level). Furthermore, if radial variations of \kappa_{\nu} are caused by particle growth, we find that the maximum size of the particle-size distribution (a_{max}) increases from sub-millimeter-sized grains in the outer disk (R>70 AU) to millimeter and centimeter-sized grains in the inner disk regions (R< 70 AU). We compare our observational constraint on a_{max}(R) with predictions from physical models of dust evolution in proto-planetary disks. For the dust composition and particle-size distribution investigated here, our observational constraints on a_{max}(R) are consistent with models where the maximum grain size is limited by radial drift.

The majority of debris discs discovered so far have only been detected through infrared excess emission above stellar photospheres. While disc properties can be inferred from unresolved photometry alone under various assumptions for the physical properties of dust grains, there is a degeneracy between disc radius and dust temperature that depends on the grain size distribution and optical properties. By resolving the disc we can measure the actual location of the dust. The launch of Herschel, with an angular resolution superior to previous far-infrared telescopes, allows us to spatially resolve more discs and locate the dust directly. Here we present the nine resolved discs around A stars between 20 and 40 pc observed by the DEBRIS survey. We use these data to investigate the disc radii by fitting narrow ring models to images at 70, 100 and 160 {\mu}m and by fitting blackbodies to full spectral energy distributions. We do this with the aim of finding an improved way of estimating disc radii for unresolved systems. The ratio between the resolved and blackbody radii varies between 1 and 2.5. This ratio is inversely correlated with luminosity and any remaining discrepancies are most likely explained by differences to the minimum size of grain in the size distribution or differences in composition. We find that three of the systems are well fit by a narrow ring, two systems are borderline cases and the other four likely require wider or multiple rings to fully explain the observations, reflecting the diversity of planetary systems.

We present resolved Herschel images of circumbinary debris disks in the alpha CrB (HD139006) and beta Tri (HD13161) systems. We find that both disks are consistent with being aligned with the binary orbital planes. Though secular perturbations from the binary can align the disk, in both cases the alignment time at the distances at which the disk is resolved is greater than the stellar age, so we conclude that the coplanarity was primordial. Neither disk can be modelled as a narrow ring, requiring extended radial distributions. To satisfy both the Herschel and mid-IR images of the alpha CrB disk, we construct a model that extends from 1-300AU, whose radial profile is broadly consistent with a picture where planetesimal collisions are excited by secular perturbations from the binary. However, this model is also consistent with stirring by other mechanisms, such as the formation of Pluto-sized objects. The beta Tri disk model extends from 50-400AU. A model with depleted (rather than empty) inner regions also reproduces the observations and is consistent with binary and other stirring mechanisms. As part of the modelling process, we find that the Herschel PACS beam varies by as much as 10% at 70um and a few % at 100um. The 70um variation can therefore hinder image interpretation, particularly for poorly resolved objects. The number of systems in which circumbinary debris disk orientations have been compared with the binary plane is now four. More systems are needed, but a picture in which disks around very close binaries (alpha CrB, beta Tri, and HD 98800, with periods of a few weeks to a year) are aligned, and disks around wider binaries (99 Her, with a 50 yr period) are misaligned, may be emerging. This picture is qualitatively consistent with the expectation that the protoplanetary disks from which the debris emerged are more likely to be aligned if their binaries have shorter periods.

It is not known whether brown dwarfs (stellar-like objects with masses less than the hydrogen-burning limit, 0.075 Msun) are formed in the same way as solar-type stars or by some other process. Here we report the clear-cut identification of a self-gravitating condensation of gas and dust with a mass in the brown-dwarf regime, made through millimeter interferometric observations. The level of thermal millimeter continuum emission detected from this object indicates a mass ~ 0.02-0.03 Msun, while the small radius < 460 AU and narrow spectral lines imply a dynamical mass of 0.015-0.02 Msun. The identification of such a pre-brown dwarf core supports models according to which brown dwarfs are formed in the same manner as hydrogen-burning stars.

This paper describes Herschel observations of the nearby (8.5pc) G5V multi-exoplanet host star 61 Vir at 70-500micron carried out by the DEBRIS survey. These reveal emission that is extended out to >15arcsec with a morphology that can be fitted by a nearly edge-on (77deg inclination) radially broad (from 30AU to >100AU) debris disk of fractional luminosity 2.7x10^-5, with two unrelated sources nearby that are more prominent at longer wavelengths. Chance alignment with a background object seen at 1.4GHz provides potential for confusion, but the star's 1.4"/yr proper motion allows Spitzer 70micron images to confirm that what we are interpreting as disk emission really is circumstellar. Although the exact shape of the disk's inner edge is not well constrained, the region inside 30AU must be significantly depleted in planetesimals. This is readily explained if there are additional planets in the 0.5-30AU region, but is also consistent with collisional erosion. We also find tentative evidence that the presence of detectable debris around nearby stars correlates with the presence of the lowest mass planets that are detectable in current radial velocity surveys. Out of an unbiased sample of the nearest 60 G stars, 11 are known to have planets, of which 6 (including 61 Vir) have planets that are all less massive than Saturn, and 4 of these have evidence for debris. The debris toward one of these planet-hosts (HD20794) is reported here for the first time. This fraction (4/6) is higher than that expected for nearby field stars (15%), and implies that systems that form low-mass planets are also able to retain bright debris disks. We suggest that this correlation could arise because such planetary systems are dynamically stable and include regions that are populated with planetesimals in the formation process where the planetesimals can remain unperturbed over Gyr timescales.

Fomalhaut is a young, nearby star that is suspected to harbor an infant planetary system, interspersed with one or more belts of dusty debris. We present far-infrared images obtained with the Herschel Space Observatory with an angular resolution between 5.7 and 36.7 arcsec at wavelengths between 70 and 500 micrometer. The images show the main debris belt in great detail. Even at high spatial resolution, the belt appears smooth. The region in between the belt and the central star is not devoid of material; thermal emission is observed here as well. Also at the location of the star, excess emission is detected. We use a dynamical model together with radiative-transfer tools to derive the parameters of the debris disk. We include detailed models of the interaction of the dust grains with radiation, for both the radiation pressure and the temperature determination. Comparing these models to the spatially resolved temperature information contained in the images allows us to place strong constraints on the presence of grains that will be blown out of the system by radiation pressure. We use this to derive the dynamical parameters of the system. The appearance of the belt points towards a remarkably active system in which dust grains are produced at a very high rate by a collisional cascade in a narrow region filled with dynamically excited planetesimals. Dust particles with sizes below the blow-out size are abundantly present. The equivalent of 2000 one-km-sized comets are destroyed every day, out of a cometary reservoir amounting to 110 Earth masses. From comparison of their scattering and thermal properties, we find evidence that the dust grains are fluffy aggregates, which indicates a cometary origin. The excess emission at the location of the star may be produced by hot dust with a range of temperatures, but may also be due to gaseous free-free emission from a stellar wind.

A list of 705 new Northern Celestial Hemisphere common proper motion pairs as derived from the Second U.S. Naval Observatory CCD Astrograph Catalog is presented, along with details of their separation, motion, brightness and colour. Objects of particular interest are noted, and the use of the catalogue for such work gauged.

The exoplanet host 51 Pegasi has a widely separated red dwarf companion. Two other, far more distant, stars are also co-moving with this star, showing that they are at least of common Galactic orbit due to a common origin.

Spectroscopically identified white dwarfs from the SDSS data release 1 are matched against SDSS data release 1 proper motion survey in the search for common proper motion companions

Aims. We present photometric data of debris disks around HIP 103389 (HD 199260), HIP 107350 (HN Peg, HD206860), and HIP 114948 (HD 219482), obtained in the context of our Herschel Open Time Key Program DUNES (DUst around NEarby Stars). Methods. We used Herschel/PACS to detect the thermal emission of the three debris disks with a 3 sigma sensitivity of a few mJy at 100 um and 160 um. In addition, we obtained Herschel/PACS photometric data at 70 um for HIP 103389. Two different approaches are applied to reduce the Herschel data to investigate the impact of data reduction on the photometry. We fit analytical models to the available spectral energy distribution (SED) data. Results. The SEDs of the three disks potentially exhibit an unusually steep decrease at wavelengths > 70 um. We investigate the significance of the peculiar shape of these SEDs and the impact on models of the disks provided it is real. Our modeling reveals that such a steep decrease of the SEDs in the long wavelength regime is inconsistent with a power-law exponent of the grain size distribution -3.5 expected from a standard equilibrium collisional cascade. In contrast, a very distinct range of grain sizes is implied to dominate the thermal emission of such disks. However, we demonstrate that the understanding of the data of faint sources obtained with Herschel is still incomplete and that the significance of our results depends on the version of the data reduction pipeline used. Conclusions. A new mechanism to produce the dust in the presented debris disks, deviations from the conditions required for a standard equilibrium collisional cascade (grain size exponent of -3.5), and/or significantly different dust properties would be necessary to explain the potentially steep SED shape of the three debris disks presented. (abridged)

MACHO project photometry shows the candidate Symbiotic Star V5485 Sagittarii's red giant component has a photometric orbital period of 270 days as evidenced by ellipsoidal variation of the lightcurve.

Optical data from the Catalina RealTime Transient Survey (CRTS) reveal an emergent eruptive in the embedded Infrared Star Cluster [C2001b] 11, within the Monoceros R2 Molecular Cloud.

The Gould Belt Legacy Survey will map star-forming regions within 500 pc, using HARP (Heterodyne Array Receiver Programme), SCUBA-2 (Submillimetre Common-User Bolometer Array 2) and POL-2 (Polarimeter 2) on the James Clerk Maxwell Telescope (JCMT). This paper describes HARP observations of the J = 3-2 transitions of 13CO and C18O towards Orion A. The 1500-resolution observations cover 5 pc of the Orion filament, including OMC1 (inc. BN-KL and Orion Bar), OMC 2/3 and OMC 4, and allow a comparative study of the molecular gas properties throughout the star-forming cloud. The filament shows a velocity gradient of ~1 km/s /pc between OMC 1, 2 and 3, and high velocity emission is detected in both isotopologues. The Orion Nebula and Bar have the largest masses and line widths, and dominate the mass and energetics of the high velocity material. Compact, spatially resolved emission from CH3CN, 13CH3OH, SO, HCOOCH3, C2H5OH, CH3CHO and CH3OCHO is detected towards the Orion Hot Core. The cloud is warm, with a median excitation temperature of ~24 K; the Orion Bar has the highest excitation temperature gas, at >80 K. The C18O excitation temperature correlates well with the dust temperature (to within 40%). The C18O emission is optically thin, and the 13CO emission is marginally optically thick; despite its high mass, OMC 1 shows the lowest opacities. A virial analysis indicates that Orion A is too massive for thermal or turbulent support, but is consistent with a model of a filamentary cloud that is threaded by helical magnetic fields. The variation of physical conditions across the cloud is reflected in the physical characteristics of the dust cores....continued

We present resolved Herschel images of a circumbinary debris disk in the 99 Herculis system. The primary is a late F-type star. The binary orbit is well characterised and we conclude that the disk is misaligned with the binary plane. Two different models can explain the observed structure. The first model is a ring of polar orbits that move in a plane perpendicular to the binary pericenter direction. We favour this interpretation because it includes the effect of secular perturbations and the disk can survive for Gyr timescales. The second model is a misaligned ring. Because there is an ambiguity in the orientation of the ring, which could be reflected in the sky plane, this ring either has near-polar orbits similar to the first model, or has a 30 degree misalignment. The misaligned ring, interpreted as the result of a recent collision, is shown to be implausible from constraints on the collisional and dynamical evolution. Because disk+star systems with separations similar to 99 Herculis should form coplanar, possible formation scenarios involve either a close stellar encounter or binary exchange in the presence of circumstellar and/or circumbinary disks. Discovery and characterisation of systems like 99 Herculis will help understand processes that result in planetary system misalignment around both single and multiple stars.

The potential for hosting photosynthetic life on Earth-like planets within binary/multiple stellar systems was evaluated by modelling the levels of photosynthetically active radiation (PAR) such planets receive. Combinations of M and G stars in: (i) close-binary systems; (ii) wide-binary systems and (iii) three-star systems were investigated and a range of stable radiation environments found to be possible. These environmental conditions allow for the possibility of familiar, but also more exotic forms of photosynthetic life, such as infrared photosynthesisers and organisms specialised for specific spectral niches.

In this paper we present a model of the Beta Leo debris disc, with an emphasis on modelling the resolved PACS images obtained as part of the Herschel key programme DEBRIS. We also present new SPIRE images of the disc at 250 microns, as well as new constraints on the disc from SCUBA-2, mid-IR and scattered light imaging. Combining all available observational constraints, we find three possible models for the Beta Leo (HD102647) debris disc: (i) A 2 component model, comprised of a hot component at 2 AU and a cold component from 15-70 AU. (ii) A 3 component model with hot dust at 2 AU, warm dust at 9 AU, and a cold component from 30-70 AU, is equally valid since the cold emission is not resolved within 30 AU. (iii) A somewhat less likely possibility is that the system consists of a single very eccentric planetesimal population, with pericentres at 2 AU and apocentres at 65 AU. Thus, despite the wealth of observational constraints significant ambiguities remain; deep mid-IR and scattered light imaging of the dust distribution within 30 AU seems the most promising method to resolve the degeneracy. We discuss the implications for the possible planetary system architecture; e.g., the 2 component model suggests planets may exist at 2-15 AU, while the 3 component model suggests planets between 2-30 AU with a stable region containing the dust belt at 9 AU, and there should be no planets between 2-65 AU in the eccentric planetesimal model. We suggest that the hot dust may originate in the disintegration of comets scattered in from the cold disc, and examine all A stars known to harbour both hot and cold dust to consider the possibility that the ratio of hot and cold dust luminosities is indicative of the intervening planetary system architecture.

Pluto's icy surface has changed colour and its atmosphere has swelled since its last closest approach to the Sun in 1989. The thin atmosphere is produced by evaporating ices, and so can also change rapidly, and in particular carbon monoxide should be present as an active thermostat. Here we report the discovery of gaseous CO via the 1.3mm wavelength J=2-1 rotational transition, and find that the line-centre signal is more than twice as bright as a tentative result obtained by Bockelee-Morvan et al. in 2000. Greater surface-ice evaporation over the last decade could explain this, or increased pressure could have caused the atmosphere to expand. The gas must be cold, with a narrow line-width consistent with temperatures around 50 K, as predicted for the very high atmosphere, and the line brightness implies that CO molecules extend up to approximately 3 Pluto radii above the surface. The upper atmosphere must have changed markedly over only a decade since the prior search, and more alterations could occur by the arrival of the New Horizons mission in 2015.

We present high angular resolution observations, taken with the Very Large Array (VLA) and Multiple Element Radio Linked Interferometer Network (MERLIN) radio telescopes, at 7mm and 4.4cm respectively, of the prototype Class 0 protostar VLA1623. At 7mm we detect two sources (VLA1623A & B) coincident with the two previously detected components at the centre of this system. The separation between the two is 1.2arcsec, or ~170AU at an assumed distance of 139pc. The upper limit to the size of the source coincident with each component of VLA1623 is ~0.7arcsec, in agreement with previous findings. This corresponds to a diameter of ~100AU at an assumed distance of 139pc. Both components show the same general trend in their broadband continuum spectra, of a steeper dust continuum spectrum shortward of 7mm and a flatter spectrum longward of this. We estimate an upper limit to the VLA1623A disc mass of <0.13Msol and an upper limit to its radius of ~50AU. The longer wavelength data have a spectral index of \alpha~0.6+/-0.3. This is too steep to be explained by optically thin free-free emission. It is most likely due to optically thick free-free emission. Alternatively, we speculate that it might be due to the formation of larger grains or planetesimals in the circumstellar disc. We estimate the mass of VLA1623B to be <0.15M$sol. We can place a lower limit to its size of ~30x7 AU, and an upper limit to its diameter of ~100AU. The longer wavelength data of VLA1623B also have a spectral index of \alpha~0.6+/-0.3. The nature of VLA1623B remains a matter of debate. It could be a binary companion to the protostar, or a knot in the radio jet from VLA1623A.

We present an analysis of debris disk data around Solar-type stars (spectral types F0-K5) using the steady-state analytical model of Wyatt et al. (2007). Models are fitted to published data from the FEPS (Meyer et al. 2006) project and various GTO programs obtained with MIPS on the Spitzer Space Telescope at 24 micron and 70 micron, and compared to a previously published analysis of debris disks around A stars using the same evolutionary model. We find that the model reproduces most features found in the data sets, noting that the model disk parameters for solar-type stars are different to those of A stars. Although this could mean that disks around Solar-type stars have different properties from their counterparts around earlier-type stars, it is also possible that the properties of disks around stars of different spectral types appear more different than they are because the blackbody disk radius underestimates the true disk radius by a factor $X_r$ which varies with spectral type. We use results from realistic grain modelling to quantify this effect for solar-type stars and for A stars. Our results imply that planetesimals around solar-type stars are on average larger than around A stars by a factor of a few but that the mass of the disks are lower for disks around FGK stars, as expected. We also suggest that discrepancies between the evolutionary timescales of 24 micron statistics predicted by our model and that observed in previous surveys could be explained by the presence of two-component disks in the samples of those surveys, or by transient events being responsible for the 24 micron emission of cold disks beyond a few Myr. Further study of the prevalence of two component disks, and of constraints on $X_r$, and increasing the size of the sample of detected disks, are important for making progress on interpreting the evolution of disks around solar-type stars.

The dangers inherent in the utilisation of easily available and readily analysed datasets in tandem with "blackbox" software applications without supporting taxonomic understanding and data nature familiarity is exemplified via examination of a misrepresentation of the results from a DEBIL analysis of OGLE II Galactic Bulge candidate variables.

We have carried out a pilot study for the SCUBA-2 'All-Sky' Survey, SASSy, a wide and shallow mapping project at 850 microns, designed to find rare objects, both Galactic and extragalactic. Two distinct sets of exploratory observations were undertaken, and used to test the SASSy approach and data reduction pipeline. The first was a 0.5 by 0.5 degrees map around the nearby galaxy NGC 2559. The galaxy was easily detected at 156 mJy, but no other convincing sources are present in the map. Comparison with other galaxies with similar wavelength coverage indicates that NGC 2559 has relatively warm dust. The second observations cover 1 square degree around the W5-E HII region. As well as diffuse structure in the map, a filtering approach was able to extract 27 compact sources with signal-to-noise greater than 6. By matching with data at other wavelengths we can see that the SCUBA-2 data can be used to discriminate the colder cores. Together these observations show that the SASSy project will be able to meet its original goals of detecting new bright sources which will be ideal for follow-up observations with other facilities.

We present a study of the impact of different model groups in the detection of circumstellar debris disks. Almost all previous studies in this field have used Kurucz model spectra to predict the stellar contribution to the flux at the wavelength of observation thus determining the existence of a disk excess. Only recently have other model groups or families like Marcs and NextGen-Phoenix become available to the same extent. This study aims to determine whether the predicted stellar flux of a disk target can change with the choice of model family - can a disk excess be present in the use of one model family whilst being absent from another? A simple comparison of Kurucz model spectra with Mrcs and NextGen model spectra of identical stellar parameters was conducted and differences were present at near-infrared wavelengths. Model spectra often do not extend in wavelength to that of observation and therefore extrapolation of the spectrum is required. In extrapolation of model spectra to the Spitzer MIPS passbands, prediction of the stellar contribution differed by 5 % at 70mum for F, G and early K spectral types with differences increasing to 15% for early M dwarfs. Analysis of the Spitzer MIPS 24mum observations of 37 F, G and K solar-like stars in the Pleiades cluster was conducted. In using Kurucz model spectra, 7 disk excesses were detected while only 3 and 4 excesses were detected in using Marcs and NextGen-Phoenix model spectra respectively.

The mass of solids in a young circumstellar disc may be the key factor in its efficiency in building planetesimals and planetary cores, and dust observed around young T Tauri and Herbig Ae stars can be used as a proxy for this initial solid content. The dust-mass distributions are taken from recent millimetre-wavelength data and fitted using survival analysis to take into account upper limits, and threshold disc-masses for building planets and belts of comets are estimated. Amongst A-stars, 20% gas giant and 55% debris disc systems are predicted, in good agreement with observations. For M-stars, the predicted and observed planet-frequencies agree at ~2-3%, and this low incidence is explained by a lack of massive discs. However, debris is predicted around approx. 14% of M-stars, while only ~2% such systems have so far been found. This suggests that deeper searches such as with Herschel and SCUBA-2 may find a cold disc population previously missed around these low-luminosity stars. Also, an estimate of the efficiency of building millimetre-detected dust into planetary cores suggests that about a third of M-stars could host an Earth-mass planet -- but as the dust is spread over large disc areas, such planets may orbit far from the star.

The thick disc contains stars formed within the first Gyr of Galactic history, and little is known about their planetary systems. The Spitzer MIPS instrument was used to search 11 of the closest of these old low-metal stars for circumstellar debris, as a signpost that bodies at least as large as planetesimals were formed. A total of 22 thick disc stars has now been observed, after including archival data, but dust is not found in any of the systems. The data rule out a high incidence of debris among star systems from early in the Galaxy's formation. However, some stars of this very old population do host giant planets, at possibly more than the general incidence among low-metal Sun-like stars. As the Solar System contains gas giants but little cometary dust, the thick disc could host analogue systems that formed many Gyr before the Sun.

It has recently been noted that many discs around T Tauri stars appear to comprise only a few Jupiter-masses of gas and dust. Using millimetre surveys of discs within six local star-formation regions, we confirm this result, and find that only a few percent of young stars have enough circumstellar material to build gas giant planets, in standard core accretion models. Since the frequency of observed exo-planets is greater than this, there is a `missing mass' problem. As alternatives to simply adjusting the conversion of dust-flux to disc mass, we investigate three other classes of solution. Migration of planets could hypothetically sweep up the disc mass reservoir more efficiently, but trends in multi-planet systems do not support such a model, and theoretical models suggest that the gas accretion timescale is too short for migration to sweep the disc. Enhanced inner-disc mass reservoirs are possible, agreeing with predictions of disc evolution through self-gravity, but not adding to millimetre dust-flux as the inner disc is optically thick. Finally, the incidence of massive discs is shown to be higher at the {\it proto}stellar stages, Classes 0 and I, where discs substantial enough to form planets via core accretion are abundant enough to match the frequency of exo-planets. Gravitational instability may also operate in the Class 0 epoch, where half the objects have potentially unstable discs of $\ga$30 % of the stellar mass. However, recent calculations indicate that forming gas giants inside 50 AU by instability is unlikely, even in such massive discs. Overall, the results presented suggest that the canonically 'proto-planetary' discs of Class II T Tauri stars {\bf have globally low masses in dust observable at millimetre wavelengths, and conversion to larger bodies (anywhere from small rocks up to planetary cores) must already have occurred.}

The Gould Belt Legacy Survey (GBS) on the JCMT has observed a region of 260 square arcminutes in 12CO J=3--2 emission, and a 190 square arcminute subset of this in 13CO and C18O towards the Serpens molecular cloud. We examine the global velocity structure of the non-outflowing gas, and calculate excitation temperatures and opacities. The large scale mass and energetics of the region are evaluated, with special consideration for high velocity gas. We find the cloud to have a mass of 203 solar masses, and to be gravitationally bound, and that the kinetic energy of the outflowing gas is approximately seventy percent of the turbulent kinetic energy of the cloud. We identify compact outflows towards some of the submillimetre Class 0/I sources in the region

We present results from the earliest observations of DEBRIS, a Herschel Key Programme to conduct a volume- and flux-limited survey for debris discs in A-type through M-type stars. PACS images (from chop/nod or scan-mode observations) at 100 and 160 micron are presented toward two A-type stars and one F-type star: beta Leo, beta UMa and eta Corvi. All three stars are known disc hosts. Herschel spatially resolves the dust emission around all three stars (marginally, in the case of beta UMa), providing new information about discs as close as 11 pc with sizes comparable to that of the Solar System. We have combined these data with existing flux density measurements of the discs to refine the SEDs and derive estimates of the fractional luminosities, temperatures and radii of the discs.

We obtained Herschel PACS and SPIRE images of the thermal emission of the debris disk around the A5V star {\beta} Pic. The disk is well resolved in the PACS filters at 70, 100, and 160 {\mu}m. The surface brightness profiles between 70 and 160 {\mu}m show no significant asymmetries along the disk, and are compatible with 90% of the emission between 70 and 160 {\mu}m originating in a region closer than 200 AU to the star. Although only marginally resolving the debris disk, the maps obtained in the SPIRE 250 - 500 {\mu}m filters provide full-disk photometry, completing the SED over a few octaves in wavelength that had been previously inaccessible. The small far-infrared spectral index ({\beta} = 0.34) indicates that the grain size distribution in the inner disk (<200AU) is inconsistent with a local collisional equilibrium. The size distribution is either modified by non-equilibrium effects, or exhibits a wavy pattern, caused by an under-abundance of impactors which have been removed by radiation pressure.

We present five band imaging of the Vega debris disc obtained using the Herschel Space Observatory. These data span a wavelength range of 70-500 um with full-width half-maximum angular resolutions of 5.6-36.9". The disc is well resolved in all bands, with the ring structure visible at 70 and 160 um. Radial profiles of the disc surface brightness are produced, and a disc radius of 11" (~ 85 AU) is determined. The disc is seen to have a smooth structure thoughout the entire wavelength range, suggesting that the disc is in a steady state, rather than being an ephemeral structure caused by the recent collision of two large planetesimals.

As part of a JCMT Legacy Survey of star formation in the Gould Belt, we present early science results for Taurus. CO J=3-2 maps have been secured along the north-west ridge and bowl, collectively known as L 1495, along with deep 13CO and C18O J=3-2 maps in two sub-regions. With these data we search for molecular outflows, and use the distribution of flows, HH objects and shocked H2 line emission features, together with the population of young stars, protostellar cores and starless condensations to map star formation across this extensive region. In total 21 outflows are identified. It is clear that the bowl is more evolved than the ridge, harbouring a greater population of T Tauri stars and a more diffuse, more turbulent ambient medium. By comparison, the ridge contains a much younger, less widely distributed population of protostars which, in turn, is associated with a greater number of molecular outflows. We estimate the ratio of the numbers of prestellar to protostellar cores in L 1495 to be ~ 1.3-2.3, and of gravitationally unbound starless cores to (gravitationally bound) prestellar cores to be ~ 1. If we take previous estimates of the protostellar lifetime of ~ 5 x 10^5 yrs, this indicates a prestellar lifetime of 9(+/-3) x 10^5 yrs. From the number of outflows we also crudely estimate the star formation efficiency in L 1495, finding it to be compatible with a canonical value of 10-15 %. We note that molecular outflow-driving sources have redder near-IR colours than their HH jet-driving counterparts. We also find that the smaller, denser cores are associated with the more massive outflows, as one might expect if mass build-up in the flow increases with the collapse and contraction of the protostellar envelope.

Numerous nearby FGK dwarfs possess discs of debris generated by collisions among comets. Here we fit the levels of dusty excess observed by Spitzer at 70$\umu$m and show that they form a rather smooth distribution. Taking into account the transition of the dust removal process from collisional to Poynting-Robertson drag, all the stars may be empirically fitted by a single population with many low-excess members. Within this ensemble, the Kuiper Belt is inferred to be such a low-dust example, among the last 10% of stars, with a small cometary population. Analogue systems hosting gas giant planets and a modest comet belt should occur for only a few per cent of Sun-like stars, and so terrestrial planets with a comparable cometary impact rate to the Earth may be uncommon. The nearest such analogue system presently known is HD154345 at 18pc, but accounting for survey completeness, a closer example should lie at around 10pc.

Debris discs - analogous to the Asteroid and Kuiper-Edgeworth belts in the Solar system - have so far mostly been identified and studied in thermal emission shortward of 100 um. The Herschel space observatory and the SCUBA-2 camera on the James Clerk Maxwell Telescope will allow efficient photometric surveying at 70 to 850 um, which allow for the detection of cooler discs not yet discovered, and the measurement of disc masses and temperatures when combined with shorter wavelength photometry. The SCUBA-2 Unbiased Nearby Stars (SUNS) survey and the DEBRIS Herschel Open Time Key Project are complimentary legacy surveys observing samples of ~500 nearby stellar systems. To maximise the legacy value of these surveys, great care has gone into the target selection process. This paper describes the target selection process and presents the target lists of these two surveys.

The Gould Belt Legacy Survey will survey nearby star-forming regions (within 500 pc), using HARP (Heterodyne Array Receiver Programme), SCUBA-2 (Submillimetre Common- User Bolometer Array 2) and POL-2 (Polarimeter 2) on the James Clerk Maxwell Telescope (JCMT). This paper describes the initial data obtained using HARP to observe 12CO, 13CO and C18O J = 3 - 2 towards two regions in Orion B, NGC 2024 and NGC 2071. We describe the physical characteristics of the two clouds, calculating temperatures and opacities utilizing all three isotopologues. We find good agreement between temperatures calculated from CO and from dust emission in the dense, energetic regions. We determine the mass and energetics of the clouds, and of the high-velocity material seen in 12CO emission, and compare the relative energetics of the high- and low-velocity material in the two clouds. We present a CLUMPFIND analysis of the 13CO condensations. The slope of the condensation mass functions, at the high-mass ends, is similar to the slope of the initial mass function.

An unbiased search for debris discs around nearby Sun-like stars is reported. Thirteen G-dwarfs at 12-15 parsecs distance were searched at 850 $\umu$m wavelength, and a disc is confirmed around HD 30495. The estimated dust mass is 0.008 M$_{\oplus}$ with a net limit $\la 0.0025$ M$_{\oplus}$ for the average disc of the other stars. The results suggest there is not a large missed population of substantial cold discs around Sun-like stars -- HD 30495 is a bright rather than unusually cool disc, and may belong to a few hundred Myr-old population of greater dust luminosity. The far-infared and millimetre survey data for Sun-like stars are well fitted by either steady state or stirred models, provided that typical comet belts are comparable in size to that in the Solar System.

This White Paper, submitted to the National Academy of Sciences' Astro2010 Decadal Review Committee, focuses on 2 central themes in the study of young brown dwarfs -- their formation mechanism and disk characteristics -- which are of direct relevance to fundamental questions of stellar and planetary origins and properties.

We report the first attempt to observe the secondary eclipse of a transiting extra-solar planet at radio wavelengths. We observed HD 189733 b with the Robert C. Byrd Green Bank Telescope of the NRAO over about 5.5 hours before, during and after secondary eclipse, at frequencies of 307 - 347 MHz. In this frequency range, we determine the 3-sigma upper limit to the flux density to be 81 mJy. The data are consistent with no eclipse or a marginal reduction in flux at the time of secondary eclipse in all subsets of our bandwidth; the strongest signal is an apparent eclipse at the 2-sigma level in the 335.2 - 339.3 MHz region. Our observed upper limit is close to theoretical predictions of the flux density of cyclotron-maser radiation from the planet.

We have imaged the disc of the young star HL Tau using the VLA at 1.3 cm, with 0.08" resolution (as small as the orbit of Jupiter). The disc is around half the stellar mass, assuming a canonical gas-mass conversion from the measured mass in large dust grains. A simulation shows that such discs are gravitationally unstable, and can fragment at radii of a few tens of AU to form planets. The VLA image shows a compact feature in the disc at 65 AU radius (confirming the `nebulosity' of Welch et al. 2004), which is interpreted as a localised surface density enhancement representing a candidate proto-planet in its earliest accretion phase. If correct, this is the first image of a low-mass companion object seen together with the parent disc material out of which it is forming. The object has an inferred gas plus dust mass of approximately 14 M(Jupiter), similar to the mass of a proto-planet formed in the simulation. The disc instability may have been enhanced by a stellar flyby: the proper motion of the nearby star XZ Tau shows it could have recently passed the HL Tau disc as close as ~600 AU.

The discovery of extra-solar planets is one of the greatest achievements of modern astronomy. The detection of planets with a wide range of masses demonstrates that extra-solar planets of low mass exist. In this paper we describe a mission, called Darwin, whose primary goal is the search for, and characterization of, terrestrial extrasolar planets and the search for life. Accomplishing the mission objectives will require collaborative science across disciplines including astrophysics, planetary sciences, chemistry and microbiology. Darwin is designed to detect and perform spectroscopic analysis of rocky planets similar to the Earth at mid-infrared wavelengths (6 - 20 micron), where an advantageous contrast ratio between star and planet occurs. The baseline mission lasts 5 years and consists of approximately 200 individual target stars. Among these, 25 to 50 planetary systems can be studied spectroscopically, searching for gases such as CO2, H2O, CH4 and O3. Many of the key technologies required for the construction of Darwin have already been demonstrated and the remainder are estimated to be mature in the near future. Darwin is a mission that will ignite intense interest in both the research community and the wider public.

The probability of a star hosting a planet that is detectable in radial velocity surveys increases Ppl(Z) oc 10^2Z, where Z is metallicity. Core accretion models reproduce this trend, since the protoplanetary disk of a high metallicity star has a high density of solids and so forms cores which accrete gas before the primordial gas disk dissipates. This paper considers the origin of the form of Ppl(Z). We introduce a simple model in which detectable planets form when the mass of solids in the protoplanetary disk, Ms, exceeds a critical value. In this model the form of Ppl(Z) is a direct reflection of the distribution of protoplanetary disk masses, Mg, and the observed Ppl(Z) is reproduced if P(Mg>Mg') oc 1/Mg'^2. We argue that a protoplanetary disk's sub-mm dust mass is a pristine indicator of the mass available for planet-building and find the observed sub-mm disk mass distribution is consistent with the observed Ppl(Z) if Ms>0.5M_J is required to form detectable planets. Any planet formation model which imposes a critical solid mass for planet formation would reproduce the observed Ppl(Z), and core accretion models are empirically consistent with a threshold criterion. We identify 7 protoplanetary disks which, by rigid application of this criterion, would be expected to form detectable planets. A testable prediction is that Ppl(Z) should flatten both for Z>0.5dex and as more distant and lower mass planets are discovered. Further, combining this model with one in which the evolution of a star's debris disk is also influenced by the solid mass in its protoplanetary disk, results in the prediction that debris disks detected around stars with planets should be more infrared luminous than those around stars without planets in tentative agreement with recent observations.

We present the scientific motivation and observing plan for an upcoming detection survey for debris disks using the James Clerk Maxwell Telescope. The SCUBA-2 Unbiased Nearby Stars (SUNS) Survey will observe 500 nearby main sequence and sub-giant stars (100 of each of the A, F, G, K and M spectral classes) to the 850 micron extragalactic confusion limit to search for evidence of submillimeter excess, an indication of circumstellar material. The survey distance boundaries are 8.6, 16.5, 22, 25 and 45 pc for M, K, G, F and A stars, respectively, and all targets lie between the declinations of -40 deg to 80 deg. In this survey, no star will be rejected based on its inherent properties: binarity, presence of planetary companions, spectral type or age. This will be the first unbiased survey for debris disks since IRAS. We expect to detect ~125 debris disks, including ~50 cold disks not detectable in current shorter wavelength surveys. A substantial amount of complementary data will be required to constrain the temperatures and masses of discovered disks. High resolution studies will likely be required to resolve many of the disks. Therefore, these systems will be the focus of future observational studies using a variety of observatories to characterize their physical properties. For non-detected systems, this survey will set constraints (upper limits) on the amount of circumstellar dust, of typically 200 times the Kuiper Belt mass, but as low as 10 times the Kuiper Belt mass for the nearest stars in the sample (approximately 2 pc).

This paper describes a James Clerk Maxwell Telescope (JCMT) legacy survey that has been awarded roughly 500 hrs of observing time to be carried out from 2007 to 2009. In this survey we will map with SCUBA-2 (Submillimetre Common User Bolometer Array 2) almost all of the well-known low-mass and intermediate-mass star-forming regions within 0.5 kpc that are accessible from the JCMT. Most of these locations are associated with the Gould Belt. From these observations we will produce a flux-limited snapshot of star formation near the Sun, providing a legacy of images, as well as point-source and extended-source catalogues, over almost 700 square degrees of sky. The resulting images will yield the first catalogue of prestellar and protostellar sources selected by submillimetre continuum emission, and should increase the number of known sources by more than an order of magnitude. We will also obtain CO maps with the array receiver HARP (Heterodyne Array Receiver Programme), in three CO isotopologues, of a large typical sample of prestellar and protostellar sources. We will then map the brightest hundred sources with the SCUBA-2 polarimeter (POL-2), producing the first statistically significant set of polarization maps in the submillimetre. The images and source catalogues will be a powerful reference set for astronomers, providing a detailed legacy archive for future telescopes, including ALMA, Herschel and JWST.

Disks in the 6 Myr old cluster eta Chamaeleontis were searched for emission from hot H2. Around the M3 star ECHAJ0843.3-7905 we detect circumstellar gas orbiting at ~2 AU. If the gas is UV-excited, the ro-vibrational line traces a hot gas layer supported by a disk of mass ~0.03Msolar, similar to the minimum mass solar nebula. Such a gas reservoir at 6 Myr would promote the formation and inwards migration of gas giant planets.

We report on multi-frequency linear polarization monitoring of 15 active galactic nuclei containing highly relativistic jets with apparent speeds from $\sim$4 $c$ to $>40c$. The measurements were obtained at optical, 1 mm, and 3 mm wavelengths, and at 7 mm with the Very Long Baseline Array. The data show a wide range in degree of linear polarization among the sources, from $<$1% to $>$30%, and interday polarization variability in individual sources. The polarization properties suggest separation of the sample into three groups with low, intermediate, and high variability of polarization in the core at 7 mm : LVP, IVP, and HVP, respectively. The groups are partially associated with the common classification of active galactic nuclei as radio galaxies and quasars with low optical polarization (LVP), BL Lacertae objects (IVP), and highly optically polarized quasars (HVP). Our study investigates correlations between total flux, fractional polarization, and polarization position angle at the different wavelengths. We interpret the polarization properties of the sources in the sample through models in which weak shocks compress turbulent plasma in the jet. The differences in the orientation of sources with respect to the observer, jet kinematics, and abundance of thermal matter external to the jet near the core can account for the diversity in the polarization properties. The results provide strong evidence that the optical polarized emission originates in shocks, most likely situated between the 3 mm and 7 mm VLBI cores. They also support the idea that the 1 mm core lies at the edge of the transition zone between electromagnetically dominated and turbulent hydrodynamical sections of the jet.

The sub-millimetre wavelength regime is perhaps the most poorly explored over large areas of the sky, despite the considerable effort that has been expended in making deep maps over small regions. As a consequence the properties of the sub-millimetre sky as a whole, and of rare bright objects in particular, remains largely unknown. Here we describe a forthcoming survey (the SCUBA-2 ``All-Sky'' Survey, or SASSy) designed to address this issue by making a large-area map of approximately one-fifth of the sky visible from the JCMT (4800 square degrees) down to a 1 sigma noise level of 30 mJy/beam. This map forms the pilot for a much larger survey, which will potentially map the remaining sky visible from the JCMT, with the region also visible to ALMA as a priority. SASSy has been awarded 500 hours for the 4800 square degree pilot phase and will commence after the commissioning of SCUBA-2, expected in early 2008.

Extrasolar planetary systems range from hot Jupiters out to icy comet belts more distant than Pluto. We explain this diversity in a model where the mass of solids in the primordial circumstellar disk dictates the outcome. The star retains measures of the initial heavy-element (metal) abundance that can be used to map solid masses onto outcomes, and the frequencies of all classes are correctly predicted. The differing dependences on metallicity for forming massive planets and low-mass cometary bodies are also explained. By extrapolation, around two-thirds of stars have enough solids to form Earth-like planets, and a high rate is supported by the first detections of low-mass exo-planets.

In this paper a simple analytical model for the steady-state evolution of debris disks due to collisions is confronted with Spitzer observations of main sequence A stars. All stars are assumed to have planetesimal belts with a distribution of initial masses and radii. In the model disk mass is constant until the largest planetesimals reach collisional equilibrium whereupon the mass falls off oc 1/t. We find that the detection statistics and trends seen at both 24 and 70um can be fitted well by the model. While there is no need to invoke stochastic evolution or delayed stirring to explain the statistics, a moderate rate of stochastic events is not ruled out. Potentially anomalous systems are identified by a high dust luminosity compared with the maximum permissible in the model (HD3003, HD38678, HD115892, HD172555). Their planetesimals may have unusual properties (high strength or low eccentricity) or this dust could be transient. While transient phenomena are also favored for a few systems in the literature, the overall success of our model, which assumes planetesimals in all belts have the same strength, eccentricity and maximum size, suggests a large degree of uniformity in the outcome of planet formation. The distribution of planetesimal belt radii, once corrected for detection bias, follows N(r) oc r^{-0.8+-0.3} for 3-120AU. Since the inner edge is often attributed to an unseen planet, this provides a unique constraint on the planetary systems of A stars. It is also shown that P-R drag may sculpt the inner edges of A star disks close to the Spitzer detection threshold (HD2262, HD19356, HD106591, HD115892). This model can be readily applied to the interpretation of future surveys, and predictions are made for the upcoming SCUBA-2 survey, including that >17% of A stars should be detectable at 850um.

Magnetic fields are an important but largely unknown ingredient of planetary nebulae. They have been detected in oxygen-rich AGB and post-AGB stars, and may play a role in the shaping of their nebulae. Here we present SCUBA sub-millimeter polarimetric observations of four bipolar planetary nebulae and post-AGB stars, including two oxygen-rich and two carbon-rich nebulae, to determine the geometry of the magnetic field by dust alignment. Three of the four sources (NGC 7027, NGC 6537 and NGC 6302) present a well-defined toroidal magnetic field oriented along their equatorial torus or disk. NGC 6302 may also show field lines along the bipolar outflow. CRL 2688 shows a complex field structure, where part of the field aligns with the torus, whilst an other part approximately aligns with the polar outflow. It also presents marked asymmetries in its magnetic structure. NGC 7027 shows evidence for a disorganized field in the south-west corner, where the SCUBA shows an indication for an outflow. The findings show a clear correlation between field orientation and nebular structure.

There is currently debate over whether the dust content of planetary systems is stochastically regenerated or originates in planetesimal belts evolving in steady state. In this paper a simple model for the steady state evolution of debris disks due to collisions is developed and confronted with the properties of the emerging population of 7 sun-like stars that have hot dust <10AU. The model shows there is a maximum possible disk mass at a given age, since more massive primordial disks process their mass faster. The corresponding maximum dust luminosity is f_max=0.00016r^(7/3)/t_age. The majority (4/7) of the hot disks exceed this limit by >1000 and so cannot be the products of massive asteroid belts, rather the following systems must be undergoing transient events characterized by an unusually high dust content near the star: eta Corvi, HD69830, HD72905 and BD+20307. It is also shown that the hot dust cannot originate in a recent collision in an asteroid belt, since there is also a maximum rate at which collisions of sufficient magnitude to reproduce a given dust luminosity can occur. Further it is shown that the planetesimal belt feeding the dust in these systems must be located further from the star than the dust, typically at >2AU. Other notable properties of the 4 hot dust systems are: two also have a planetesimal belt at >10AU (eta Corvi and HD72905); one has 3 Neptune mass planets at <1AU (HD69830); all exhibit strong silicate features in the mid-IR. We consider the most likely origin for the dust in these systems to be a dynamical instability which scattered planetesimals inwards from a more distant planetesimal belt in an event akin to the Late Heavy Bombardment in our own system, the dust being released from such planetesimals in collisions and possibly also sublimation.

The evidence for a rotation of the epsilon Eridani debris disc is examined. Data at 850 micron wavelength were previously obtained using the Submillimetre Common User Bolometer Array (SCUBA) over periods in 1997-1998 and 2000-2002. By chi-square fitting after shift and rotation operations, images from these two epochs were compared to recover proper motion and orbital motion of the disc. The same procedures were then performed on simulated images to estimate the accuracy of the results. Minima in the chi-square plots indicate a motion of the disc of approximately 0.6'' per year in the direction of the star's proper motion. This underestimates the true value of 1'' per year, implying that some of the structure in the disc region is not associated with epsilon Eridani, originating instead from background galaxies. From the chi-square fitting for orbital motion, a counterclockwise rotation rate of ~2.75 degrees per year is deduced. Comparisons with simulated data in which the disc is not rotating show that noise and background galaxies result in approximately Gaussian fluctuations with a standard deviation +/-1.5 degrees per year. Thus counterclockwise rotation of disc features is supported at approximately a 2-sigma level, after a 4-year time difference. This rate is faster than the Keplerian rate of 0.65 degrees per year for features at ~65 AU from the star, suggesting their motion is tracking a planet inside the dust ring. Future observations with SCUBA-2 can rule out no rotation of the epsilon Eridani dust clumps with ~4-sigma confidence. Assuming a rate of about 2.75 degrees per year, the rotation of the features after a 10-year period could be shown to be >1 degree per year at the 3-sigma level.

We describe results from a survey for J=3-2 12CO emission from visible stars with an infrared excess. The line is clearly detected in 21 objects, with molecular gas (>10^-3 Jupiter masses) common in targets with infrared excesses >0.01 (>56% of objects). Such high excesses indicate the presence of a disc of opening angle >12 degrees; within this, the optically thick disc prevents CO photodissociation. Two or three stars with associated CO have an excess <0.01, implying a disc opening angle <1 degree. Most line profiles are double-peaked or relatively broad. Model fits, assuming a Keplerian disc, indicate outer radii, R_out, of ~20-300 au. As many as 5 discs have outer radii smaller than the Solar System (50 au), and a further 4 have gas at radii <20 au. R_out is independent of the stellar spectral type (from K through to B9), but is correlated with total dust mass. R_out appears to decrease with time: discs around stars of age 3-7 Myr have a mean radius of ~210 au, whereas discs of age 7-20 Myr are a factor of 3 smaller. The only bona fide debris disc with detected CO is HD9672; this has a double peaked line profile and is the most compact gas disc observed, with a modelled radius 17 au). A fit to HD141569 suggests the gas lies in two rings of radii 90 and 250 au, similar to the scattered light structure. In both AB Aur and HD163296 the sizes of the molecular and dust scattering discs are also similar, suggesting that the gas and small dust grains are co-located.

We present sub-millimeter and mid-infrared images of the circumstellar disk around the nearby F2V star eta Corvi. The disk is resolved at 850um with a size of ~100AU. At 450um the emission is found to be extended at all position angles, with significant elongation along a position angle of 130+-10deg; at the highest resolution (9.3") this emission is resolved into two peaks which are to within the uncertainties offset symmetrically from the star at 100AU projected separation. Modeling the appearance of emission from a narrow ring in the sub-mm images shows the observed structure cannot be caused by an edge-on or face-on axisymmetric ring; the observations are consistent with a ring of radius 150+-20AU seen at 45+-25deg inclination. More face-on orientations are possible if the dust distribution includes two clumps similar to Vega; we show how such a clumpy structure could arise from the migration over 25Myr of a Neptune mass planet from 80-105AU. The inner 100AU of the system appears relatively empty of sub-mm emitting dust, indicating that this region may have been cleared by the formation of planets, but the disk emission spectrum shows that IRAS detected an additional hot component with a characteristic temperature of 370+-60K (implying a distance of 1-2AU). At 11.9um we found the emission to be unresolved with no background sources which could be contaminating the fluxes measured by IRAS. The age of this star is estimated to be ~1Gyr. It is very unusual for such an old main sequence star to exhibit significant mid-IR emission. The proximity of this source makes it a perfect candidate for further study from optical to mm wavelengths to determine the distribution of its dust.

The results of a submillimetre wavelength spectral line survey between 455.1 -- 507.4 GHz of the Orion-KL hot cloud core are reported. A total of 254 lines were detected to a main beam brightness temperature sensitivity ~ 1 - 3 K. The detected lines are identified as being associated with 30 different molecular species or their isotopomeric variants. The strongest line detected was the J = 4-3 transition of the CO molecule. Apart from abundant diatomic rotors such as CO and CS, the spectrum is dominated by SO, SO_2 and CH_3OH and large organic molecules such as (CH_3)_2O, CH$_3CN, C_2H_3CN, C_2H_5CN and HCOOCH_3 which make up ~ 72% of the total number of lines; unidentified lines ~ 13%; and other lines the remaining ~ 15% of the total. Rotational temperatures and column densities derived using standard rotation diagram analysis techniques were found to range from 70 - 300 K, and 10^14 - 10^17 cm^2 respectively. Accepted by Astronomy & Astrophysics 22 May 2003

We have observed 22 Lindroos stars with SCUBA on the JCMT to search for evidence of dust discs. Stars in this sample are the less massive companions of B-type primaries and have ages of 10-170Myr. Dust was detected around three of these stars (HD74067, HD112412 and HD99803B). The emission around HD74067 is centrally peaked and is approximately symmetrically distributed out to ~70". This emission either arises from a two component disc, one circumstellar and the other circumbinary with dust masses of 0.3 and >27Mearth respectively, or an unrelated background object. The other two detections we attribute to circumsecondary discs with masses of 0.04 and 0.3Mearth; a circumprimary disc is also present around HD112413 with a similar mass to that around the companion HD112412. Cross-correlation of our sample with the IRAS catalogs only showed evidence for dust emission at 25um and 60um toward one star (HD1438); none of the sub-mm detections were evident in the far-IR data implying that these discs are cold (<40K assuming beta=1). Our sub-mm detections are some of the first of dust discs surrounding evolved stars that were not detected by IRAS or ISO and imply that 9-14% of stars could harbour previously undetected dust discs that await discovery in unbiased sub-mm surveys. If these discs are protoplanetary remnants, rather than secondary debris discs, dust lifetime arguments show that they must be devoid of small <0.1mm grains. Thus it may be possible to determine the origin of these discs from their spectral energy distributions. The low dust masses for this sample support the picture that protoplanetary dust discs are depleted to the levels of the brightest debris discs (~1Mearth) within 10Myr.

A polarimeter has been built for use with the Submillimetre Common-User Bolometer Array (SCUBA), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. SCUBA is the first of a new generation of highly sensitive submillimetre cameras, and the UK/Japan Polarimeter adds a polarimetric imaging/photometry capability in the wavelength range 350 to 2000 microns. Early science results range from measuring the synchrotron polarization of the black hole candidate Sgr A* to mapping magnetic fields inferred from polarized dust emission in Galactic star-forming clouds. We describe the instrument design, performance, observing techniques and data reduction processes, along with an assessment of the current and future scientific capability.

We present new, high-sensitivity sub-millimeter observations towards 55 Cancri, a nearby G8 star with one, or possibly two, known planetary companion(s). Our 850 $\mu$m map, obtained with the SCUBA instrument on the James Clerk Maxwell Telescope, shows three peaks of emission at the 2.5 mJy level in the vicinity of the star's position. However, the observed peaks are 25$\arcsec$--40$\arcsec$ away from the star and a deep $R$-band optical image reveals faint point sources that coincide with two of the sub-millimeter peaks. Thus, we do not find evidence for dust emission spatially associated with 55 Cancri. The excess 60 $\mu$m emission detected with ISO may originate from one or more of the 850 $\mu$m peaks that we attribute to background sources. Our new results, together with the HST/NICMOS coronographic images in the near-infrared, place stringent limits on the amount of dust in this planetary system, and argue against the existence of a detectable circumstellar dust disk around 55 Cnc.

Imaging polarimetry of the 0.85mm continuum emission in the NGC 7538 region, obtained with the SCUBA Polarimeter, is presented. The polarization map is interpreted in terms of thermal radiation by magnetically aligned dust grains. Two prominent cores associated with IRS 1 and IRS 11, IRS 1(SMM) and IRS 11(SMM), are found in the surface brightness map. Although these cores look similar in surface brightness, their polarization shows striking differences. In IRS 11(SMM), the polarization vectors are extremely well ordered, and the degrees of polarization are quite high with an average of ~3.9 %. In IRS 1(SMM), on the other hand, the directions of polarization vectors are locally disturbed, and the degrees of polarization are much lower than those of IRS 11(SMM). These differences suggest that small scale fluctuations of the magnetic field are more prominent in IRS 1(SMM). This can be interpreted in terms of the difference in evolutionary stage of the cores. Inside IRS 1(SMM), which seems to be at a later evolutionary stage than IRS 11(SMM), substructures such as subclumps or a cluster of infrared sources have already formed. Small scale fluctuations in the magnetic field could have developed during the formation of these substructures. The distribution of magnetic field directions derived from our polarization map agrees well with those of molecular outflows associated with IRS 1(SMM) and IRS 11(SMM). Comparisons of energy densities between the magnetic field and the outflows show that the magnetic field probably plays an important role in guiding the directions of the outflows.

Molecular outflows from young protostars are widely believed to be collimated by magnetic fields, but there has been little observational evidence to support this hypothesis. Using the new technique of millimetre-wavelength spectro-polarimetry, we demonstrate the existence of a magnetic field in the NGC2024-FIR5 outflow lobe. The 1.3mm J=2-1 transition of carbon monoxide (CO) is polarized at a level of approximately 1%, in a direction within 10-15 degrees of the outflow axis. This agrees with theoretical models where the magnetic field channels the outflowing gas, and shows that the process can be effective as far as 0.1pc from the protostar.

We present first results from a submillimetre continuum survey of the Galactic Centre `Central Molecular Zone' (CMZ), made with SCUBA on the James Clerk Maxwell Telescope. SCUBA's scan-map mode has allowed us to make extremely wide-field maps of thermal dust emission with unprecedented speed and sensitivity. We also discuss some issues related to the elimination of artefacts in scan-map data. Our simultaneous 850/450 micron maps have a total size of approximately 2.8 x 0.5 degrees (400 x 75 pc) elongated along the galactic plane. They cover the Sgr A region-including Sgr A*, the circumnuclear disc, and the +20 km/s and +50 km/s clouds; the area around the Pistol; Sgr B2-the brightest feature on the map; and at their Galactic Western and Eastern edges the Sgr C and Sgr D regions. There are many striking features such as filaments and shell-like structures, as well as point sources such as Sgr A* itself. The total mass in the Central Molecular Zone is greater than that revealed in previous optically-thin molecular line maps by a factor of ~3, and new details are revealed on scales down to 0.33 pc across this 400 pc wide region.

We present the first published maps of magnetic fields in pre-stellar cores, to test theoretical ideas about the way in which the magnetic field geometry affects the star formation process. The observations are JCMT-SCUBA maps of 850 micron thermal emission from dust. Linear polarizations at typically ten or more independent positions in each of three objects, L1544, L183 and L43 were measured, and the geometries of the magnetic fields in the plane of the sky were mapped from the polarization directions. The observed polarizations in all three objects appear smooth and fairly uniform. In L1544 and L183 the mean magnetic fields are at an angle of around 30 degrees to the minor axes of the cores. The L43 B-field appears to have been influenced in its southern half, such that it is parallel to the wall of a cavity produced by a CO outflow from a nearby T Tauri star, whilst in the northern half the field appears less disturbed and has an angle of 44 degrees to the core minor axis. We briefly compare our results with published models of magnetized cloud cores and conclude that no current model can explain these observations simultaneously with previous ISOCAM data.

We report the detection of linear polarization from SgrA* at 750, 850, 1350 and 2000 microns which confirms the contribution of synchrotron radiation. From the lack of polarization at longer wavelengths it appears to arise in the millimetre/sub-millimetre excess. There are large position angle changes between the millimetre and sub-millimetre results and these are discussed in terms of a polarized dust contribution in the sub-millimetre and various synchrotron models. In the model which best explains the data the synchrotron radiation from the excess is self-absorbed in the millimetre region and becomes optically thin in the sub-millimetre. This implies that the excess arises in an extremely compact source to the order of 2 Schwarzschild radii.

The presence of debris disks around $\sim$ 1-Gyr-old main sequence stars suggests that an appreciable amount of dust may persist even in mature planetary systems. Here we report the detection of dust emission from 55 Cancri, a star with one, or possibly two, planetary companions detected through radial velocity measurements. Our observations at 850$\mu$m and 450$\mu$m imply a dust mass of 0.0008-0.005 Earth masses, somewhat higher than that in the the Kuiper Belt of our solar system. The estimated temperature of the dust grains and a simple model fit both indicate a central disk hole of at least 10 AU in radius. Thus, the region where the planets are detected is likely to be significantly depleted of dust. Our results suggest that far-infrared and sub-millimeter observations are powerful tools for probing the outer regions of extrasolar planetary systems.

Linearly polarized, non-masing, rotational lines have been detected for the first time in the interstellar medium. This effect occurs in molecular clouds with a magnetic field, and traces the field direction, offering an alternative technique to dust emission polarimetry. The line polarization mechanism is similar to that in masers, but with lower degrees of polarization (~1%). We have detected the effect towards the Galactic Centre and its surrounding `2 pc ring', and in the molecular clouds S140 and DR21 (tentatively), in the lines CO J=2-1 and J=3-2, and 13CO J=2-1. The deduced magnetic field directions agree well with previous dust polarimetry results, confirming that the line polarization is a real effect. This new technique will be useful in sources that are too faint for dust polarimetry, and can also be used to investigate 3-dimensional morphology of magnetic fields, where the velocity structure of the clouds is known.

Dust emission around the nearby star epsilon Eridani has been imaged using a new submillimetre camera (SCUBA at the JCMT). At 850 microns wavelength a ring of dust is seen, peaking at 60 AU from the star and with much lower emission inside 30 AU. The mass of the ring is at least 0.01 Earth masses in dust, while an upper limit of 0.4 Earth masses in molecular gas is imposed by CO observations. The total mass is comparable to the estimated amount of material, 0.04-0.3 Earth masses, in comets orbiting the Solar System. The most probable origin of the the ring structure is that it is a young analogue to the Kuiper Belt in our Solar System, and that the central region has been partially cleared by the formation of grains into planetesimals. Dust clearing around epsilon Eri is seen within the radius of Neptune's orbit, and the peak emission at 35-75 AU lies within the estimated Kuiper Belt zone of 30-100 AU radius. epsilon Eri is a main-sequence star of type K2V (0.8 Solar masses) with an estimated age of 0.5-1.0 Gyr, so this interpretation is consistent with the early history of the Solar System where heavy bombardment occurred up to approximately 0.6 Gyr. An unexpected discovery is substructure within the ring, and these asymmetries could be due to perturbations by planets.