Gas detection around main sequence stars is becoming more common with around 20 systems showing the presence of CO. However, more detections are needed, especially around later spectral type stars to better understand the origin of this gas and refine our models. To do so, we carried out a survey of 10 stars with predicted high likelihoods of secondary CO detection using ALMA in band 6. We looked for continuum emission of mm-dust as well as gas emission (CO and CN transitions). The continuum emission was detected in 9/10 systems for which we derived the discs' dust masses and geometrical properties, providing the first mm-wave detection of the disc around HD 106906, the first mm-wave radius for HD 114082, 117214, HD 15745, HD 191089 and the first radius at all for HD 121191. A crucial finding of our paper is that we detect CO for the first time around the young 10-16 Myr old G1V star HD 129590, similar to our early Sun. The gas seems colocated with its planetesimal belt and its total mass is likely between $2-10 \times 10^{-5}$ M$_\oplus$. This first gas detection around a G-type main-sequence star raises questions as to whether gas may have been released in the Solar System as well in its youth, which could potentially have affected planet formation. We also detected CO gas around HD 121191 at a higher S/N than previously and find that the CO lies much closer-in than the planetesimals in the system, which could be evidence for the previously suspected CO viscous spreading owing to shielding preventing its photodissociation. Finally, we make estimates for the CO content in planetesimals and the HCN/CO outgassing rate (from CN upper limits), which we find are below the level seen in Solar System comets in some systems.

$\pi$ Men hosts a transiting super Earth ($P\approx6.27$ d, $m\approx4.82$ $M_{\oplus}$, $R\approx2.04$ $R_{\oplus}$) discovered by TESS and a cold Jupiter ($P\approx2093$ d, $m \sin I\approx10.02$ $M_{\rm{Jup}}$, $e\approx0.64$) discovered from radial velocity. We use Gaia DR2 and Hipparcos astrometry to derive the star's velocity caused by the orbiting planets and constrain the cold Jupiter's sky-projected inclination ($I_b=41-65^{\circ}$). From this we derive the mutual inclination ($\Delta I$) between the two planets, and find that $49^{\circ}< \Delta I < 131^{\circ}$ (1$\sigma$), and $28^{\circ} < \Delta I < 152^{\circ}$ (2$\sigma$). We examine the dynamics of the system using $N$-body simulations, and find that potentially large oscillations in the super Earth's eccentricity and inclination are suppressed by General Relativistic precession. However, nodal precession of the inner orbit around the invariable plane causes the super Earth to only transit between 7-22 per cent of the time, and to usually be observed as misaligned with the stellar spin axis. We repeat our analysis for HAT-P-11, finding a large $\Delta I$ between its close-in Neptune and cold Jupiter and similar dynamics. $\pi$ Men and HAT-P-11 are prime examples of systems where dynamically hot outer planets excite their inner planets, with the effects of increasing planet eccentricities, planet-star misalignments, and potentially reducing the transit multiplicity. Formation of such systems likely involves scattering between multiple giant planets or misaligned protoplanetary discs. Future imaging of the faint debris disc in $\pi$ Men and precise constraints on the stellar spin orientation would provide strong tests for these formation scenarios.

We present results of using individual galaxies' effective redshift probability density information as a method of identifying potential catastrophic outliers in empirical photometric redshift estimation. In the course we develop a method of modification of the redshift distribution of training sets to improve both the baseline accuracy of high redshift (z>1.5) estimation as well as catastrophic outlier mitigation. We demonstrate these using two real test data sets and one simulated test data set spanning a wide redshift range (0<z<4). We present these results in the context of "photometric redshift realism" and aim to show that the methods and results presented here can inform a 'prescription' that can be applied as a realistic photometric redshift estimation scenario for a hypothetical large-scale survey. We find that with appropriate optimization, we can identify a large percentage (>30%) of catastrophic outlier galaxies while simultaneously incorrectly flagging only a small percentage (<7% and in many cases <3%) of non-outlier galaxies as catastrophic outliers. We find also that our training set redshift distribution modification results in a significant decrease (>10%) in the percentage of outlier galaxies greater than z=1.5 with only a small increase (<3%) in the percentage of outlier galaxies less than z=1.5 compared to the unmodified training set. In addition, we find that this modification can in some cases decrease the percentage of incorrectly identified non-outlier galaxies by almost 20%, while in other cases cause only a small (<1%) increase in this metric.

The SEDs of some nearby stars show mid-infrared excesses from warm habitable zone dust, known as exozodiacal dust. This dust may originate in collisions in a planetesimal belt before being dragged inwards. This paper presents an analytical model for the size distribution of particles at different radial locations in such a scenario, considering evolution due to destructive collisions and Poynting-Robertson (P-R) drag. Results from more accurate but computationally expensive numerical simulations of this process are used to validate the model and fit its free parameters. The model predicts 11 $\mu$m excesses ($R_{11}$) for discs with a range of dust masses and planetesimal belt radii using realistic grain properties. We show that P-R drag should produce exozodiacal dust levels detectable with the Large Binocular Telescope Interferometer (LBTI) ($R_{11} > 0.1\%$) in systems with known outer belts; non-detection may indicate dust depletion, e.g. by an intervening planet. We also find that LBTI could detect exozodiacal dust dragged in from a belt too faint to detect at far-infrared wavelengths, with fractional luminosity $f\sim 10^{-7}$ and radius $\sim 10-80$ au. Application to systems observed with LBTI shows that P-R drag can likely explain most (5/9) of the exozodiacal dust detections in systems with known outer belts; two systems ($\beta$ Uma and $\eta$ Corvi) with bright exozodi may be due to exocomets. We suggest that the three systems with exozodiacal dust detections but no known belt may have cold planetesimal belts too faint to be detectable in the far-infrared. Even systems without outer belt detections could have exozodiacal dust levels $R_{11} > 0.04\%$ which are problematic for exo-Earth imaging.

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.

The Large Binocular Telescope Interferometer (LBTI) enables nulling interferometric observations across the N band (8 to 13 um) to suppress a star's bright light and probe for faint circumstellar emission. We present and statistically analyze the results from the LBTI/HOSTS (Hunt for Observable Signatures of Terrestrial Systems) survey for exozodiacal dust. By comparing our measurements to model predictions based on the Solar zodiacal dust in the N band, we estimate a 1 sigma median sensitivity of 23 zodis for early type stars and 48 zodis for Sun-like stars, where 1 zodi is the surface density of habitable zone (HZ) dust in the Solar system. Of the 38 stars observed, 10 show significant excess. A clear correlation of our detections with the presence of cold dust in the systems was found, but none with the stellar spectral type or age. The majority of Sun-like stars have relatively low HZ dust levels (best-fit median: 3 zodis, 1 sigma upper limit: 9 zodis, 95% confidence: 27 zodis based on our N band measurements), while ~20% are significantly more dusty. The Solar system's HZ dust content is consistent with being typical. Our median HZ dust level would not be a major limitation to the direct imaging search for Earth-like exoplanets, but more precise constraints are still required, in particular to evaluate the impact of exozodiacal dust for the spectroscopic characterization of imaged exo-Earth candidates.

The presence of CO gas around 10-50 Myr old A stars with debris discs has sparked debate on whether the gas is primordial or secondary. Since secondary gas released from planetesimals is poor in H$_2$, it was thought that CO would quickly photodissociate never reaching the high levels observed around the majority of A stars with bright debris discs. Kral et al. 2019 showed that neutral carbon produced by CO photodissociation can effectively shield CO and potentially explain the high CO masses around 9 A stars with bright debris discs. Here we present a new model that simulates the gas viscous evolution, accounting for carbon shielding and how the gas release rate decreases with time as the planetesimal disc loses mass. We find that the present gas mass in a system is highly dependant on its evolutionary path. Since gas is lost on long timescales, it can retain a memory of the initial disc mass. Moreover, we find that gas levels can be out of equilibrium and quickly evolving from a shielded onto an unshielded state. With this model, we build the first population synthesis of gas around A stars, which we use to constrain the disc viscosity. We find a good match with a high viscosity ($\alpha\sim0.1$), indicating that gas is lost on timescales $\sim1-10$ Myr. Moreover, our model also shows that high CO masses are not expected around FGK stars since their planetesimal discs are born with lower masses, explaining why shielded discs are only found around A stars. Finally, we hypothesise that the observed carbon cavities could be due to radiation pressure or accreting planets.

Low mass stars might offer today the best opportunities to detect and characterise planetary systems, especially those harbouring close-in low mass temperate planets. Among those stars, TRAPPIST-1 is exceptional since it has seven Earth-sized planets, of which three could sustain liquid water on their surfaces. Here we present new and deep ALMA observations of TRAPPIST-1 to look for an exo-Kuiper belt which can provide clues about the formation and architecture of this system. Our observations at 0.88 mm did not detect dust emission, but can place an upper limit of 23$\mu$Jy if the belt is smaller than 4 au, and 0.15 mJy if resolved and 100 au in radius. These limits correspond to low dust masses of $\sim10^{-5}-10^{-2}$ $M_\oplus$, which are expected after 8 Gyr of collisional evolution unless the system was born with a $>20$ $M_\oplus$ belt of 100 km-sized planetesimals beyond 40 au or suffered a dynamical instability. This $20$ $M_\oplus$ mass upper limit is comparable to the combined mass in TRAPPIST-1 planets, thus it is possible that most of the available solid mass in this system was used to form the known planets. A similar analysis of the ALMA data on Proxima Cen leads us to conclude that a belt born with a mass $\gtrsim1$ $M_\oplus$ in 100 km-sized planetesimals could explain its putative outer belt at 30 au. We recommend that future characterisations of debris discs around low mass stars should focus on nearby and young systems if possible.

The Habitable Exoplanet Observatory, or HabEx, has been designed to be the Great Observatory of the 2030s. For the first time in human history, technologies have matured sufficiently to enable an affordable space-based telescope mission capable of discovering and characterizing Earthlike planets orbiting nearby bright sunlike stars in order to search for signs of habitability and biosignatures. Such a mission can also be equipped with instrumentation that will enable broad and exciting general astrophysics and planetary science not possible from current or planned facilities. HabEx is a space telescope with unique imaging and multi-object spectroscopic capabilities at wavelengths ranging from ultraviolet (UV) to near-IR. These capabilities allow for a broad suite of compelling science that cuts across the entire NASA astrophysics portfolio. HabEx has three primary science goals: (1) Seek out nearby worlds and explore their habitability; (2) Map out nearby planetary systems and understand the diversity of the worlds they contain; (3) Enable new explorations of astrophysical systems from our own solar system to external galaxies by extending our reach in the UV through near-IR. This Great Observatory science will be selected through a competed GO program, and will account for about 50% of the HabEx primary mission. The preferred HabEx architecture is a 4m, monolithic, off-axis telescope that is diffraction-limited at 0.4 microns and is in an L2 orbit. HabEx employs two starlight suppression systems: a coronagraph and a starshade, each with their own dedicated instrument.

Dippers are a common class of young variable star exhibiting day-long dimmings with depths of up to several tens of percent. A standard explanation is that dippers host nearly edge-on (70 deg) protoplanetary discs that allow close-in (< 1 au) dust lifted slightly out of the midplane to partially occult the star. The identification of a face-on dipper disc and growing evidence of inner disc misalignments brings this scenario into question. Thus we uniformly (re)derive the inclinations of 24 dipper discs resolved with (sub-)mm interferometry from ALMA. We find that dipper disc inclinations are consistent with an isotropic distribution over 0-75 deg, above which the occurrence rate declines (likely an observational selection effect due to optically thick disc midplanes blocking their host stars). These findings indicate that the dipper phenomenon is unrelated to the outer (>10 au) disc resolved by ALMA and that inner disc misalignments may be common during the protoplanetary phase. More than one mechanism may contribute to the dipper phenomenon, including accretion-driven warps and "broken" discs caused by inclined (sub-)stellar or planetary companions.

The accretion of planetesimals onto white dwarf atmospheres allows determination of the composition of this polluting material. This composition is usually inferred from observed pollution levels by assuming it originated from a single body. This paper instead uses a stochastic model wherein polluting planetesimals are chosen randomly from a mass distribution, finding that the single body assumption is invalid in >20% of cases. Planetesimal compositions are modelled assuming parent bodies that differentiated into core, mantle and crust components. Atmospheric levels of Ca, Mg and Fe in the model are compared to a sample of 230 DZ white dwarfs for which such pollution is measured. A good fit is obtained when each planetesimal has its core, mantle and crust fractions chosen independently from logit-normal distributions which lead to average mass fractions of $f_\text{Cru}=0.15$, $f_\text{Man}=0.49$ and $f_\text{Cor}=0.36$. However, achieving this fit requires a factor 4 depletion of Mg relative to stellar material. This depletion is unlikely to originate in planetesimal formation processes, but might occur from heating while the star is on the giant branch. Alternatively the accreted material has stellar abundance, and either the inferred low Mg abundance was caused by an incorrect assumption that Mg sinks slower than Ca and Fe, or there are unmodelled biases in the observed sample. Finally, the model makes predictions for the timescale on which the observed pollutant composition varies, which should be the longer of the sinking and disc timescales, implying variability on decadal timescales for DA white dwarfs.

This paper considers how planetesimal impacts affect planetary atmospheres. Atmosphere evolution depends on the ratio of gain from volatiles to loss from atmosphere stripping f_v; for constant bombardment, atmospheres with f_v<1 are destroyed in finite time, but grow linearly with time for f_v>1. An impact outcome prescription is used to characterise how f_v depends on planetesimal impact velocities, size distribution and composition. Planets that are low mass and/or close to the star have atmospheres that deplete in impacts, while high mass and/or distant planets grow secondary atmospheres. Dividing these outcomes is an fv=1 impact shoreline analogous to Zahnle & Catling's cosmic shoreline. The impact shoreline's location depends on assumed impacting planetesimal properties, so conclusions for the atmospheric evolution of a planet like Earth with f_v~1 are only as strong as those assumptions. Application to the exoplanet population shows the gap in the planet radius distribution at ~1.5R_earth is coincident with the impact shoreline, which has a similar dependence on orbital period and stellar mass to the observed gap. Given sufficient bombardment, planets below the gap would be expected to lose their atmospheres, while those above could have atmospheres enhanced in volatiles. The level of atmosphere alteration depends on the total bombardment a planet experiences, and so on the system's (usually unknown) other planets and planetesimals, though massive distant planets would have low accretion efficiency. Habitable zone planets around lower luminosity stars are more susceptible to atmosphere stripping, disfavouring M stars as hosts of life-bearing planets if Earth-like bombardment is conducive to the development of life.

Extrasolar debris disks are the dust disks found around nearby main sequence stars arising from the break-up of asteroids and comets orbiting the stars. Far-IR surveys (e.g., with Herschel) showed that ~20% of stars host detectable dust levels. While dust temperatures suggest a location at 10s of au comparable with our Kuiper belt, orders of magnitude more dust is required implying a planetesimal population more comparable with the primordial Kuiper belt. High resolution imaging (e.g., with ALMA) has mapped the nearest and brightest disks, providing evidence for structures shaped by an underlying planetary system. Some of these are analogous to structures in our own Kuiper belt (e.g., the hot and cold classical, resonant, scattered disk and cometary populations), while others have no Solar System counterpart. CO gas is seen in some debris disks, and inferred to originate in the destruction of planetesimals with a similar volatile-rich composition to Solar System comets. This chapter reviews our understanding of extrasolar Kuiper belts and of how our own Kuiper belt compares with those of neighbouring stars.

The Wide Field Infrared Survey Telescope (WFIRST) Coronagraph Instrument (CGI) will be the first high-performance stellar coronagraph using active wavefront control for deep starlight suppression in space, providing unprecedented levels of contrast, spatial resolution, and sensitivity for astronomical observations in the optical. One science case enabled by the CGI will be taking images and(R~50)spectra of faint interplanetary dust structures present in the habitable zone of nearby sunlike stars (~10 pc) and within the snow-line of more distant ones(~20pc), down to dust density levels commensurate with that of the solar system zodiacal cloud. Reaching contrast levels below~10-7 for the first time, CGI will cross an important threshold in debris disks physics, accessing disks with low enough optical depths that their structure is dominated by transport phenomena than collisions. Hence, CGI results will be crucial for determining how exozodiacal dust grains are produced and transported in low-density disks around mature stars. Additionally, CGI will be able to measure the brightness level and constrain the degree of asymmetry of exozodiacal clouds around individual nearby sunlike stars in the optical, at the ~10x solar zodiacal emission level. This information will be extremely valuable for optimizing the observational strategy of possible future exo-Earth direct imaging missions, especially those planning to operate at optical wavelengths, such as Habitable Exoplanet Observatory (HabEx) and the Large Ultraviolet/Optical/Infrared Surveyor (LUVOIR).

Exoplanet science is one of the most thriving fields of modern astrophysics. A major goal is the atmospheric characterization of dozens of small, terrestrial exoplanets in order to search for signatures in their atmospheres that indicate biological activity, assess their ability to provide conditions for life as we know it, and investigate their expected atmospheric diversity. None of the currently adopted projects or missions, from ground or in space, can address these goals. In this White Paper we argue that a large space-based mission designed to detect and investigate thermal emission spectra of terrestrial exoplanets in the MIR wavelength range provides unique scientific potential to address these goals and surpasses the capabilities of other approaches. While NASA might be focusing on large missions that aim to detect terrestrial planets in reflected light, ESA has the opportunity to take leadership and spearhead the development of a large MIR exoplanet mission within the scope of the "Voyage 2050" long-term plan establishing Europe at the forefront of exoplanet science for decades to come. Given the ambitious science goals of such a mission, additional international partners might be interested in participating and contributing to a roadmap that, in the long run, leads to a successful implementation. A new, dedicated development program funded by ESA to help reduce development and implementation cost and further push some of the required key technologies would be a first important step in this direction. Ultimately, a large MIR exoplanet imaging mission will be needed to help answer one of mankind's most fundamental questions: "How unique is our Earth?"

We compile a sample of 341 binary and multiple star systems with the aim of searching for and characterising Kuiper belt-like debris discs. The sample is assembled by combining several smaller samples studied in previously published work with targets from two unpublished Herschel surveys. We find that 38 systems show excess emission at 70 or 100 $\mu$m suggestive of a debris disc. While nine of the discs appear to be unstable to perturbations from their host binary based on a simple analysis of their inferred radii, we argue that the evidence for genuine instability is not strong, primarily because of uncertainty in the true disc radii, uncertainty in the boundaries of the unstable regions, and orbital projection effects. The binary separation distributions of the disc-bearing and disc-free systems are different at a confidence level of $99.4\%$, indicating that binary separation strongly influences the presence of detectable levels of debris. No discs are detected for separations between $\sim$25 and 135 au; this is likely a result of binaries whose separations are comparable with typical disc radii clearing out their primordial circumstellar or circumbinary material via dynamical perturbations. The disc detection rate is $19^{+5}_{-3}\%$ for binaries wider than 135 au, similar to published results for single stars. Only $8^{+2}_{-1}\%$ of systems with separations below 25 au host a detectable disc, which may suggest that planetesimal formation is inhibited in binaries closer than a few tens of au, similar to the conclusions of studies of known planet-hosting binaries.

Debris disks are extrasolar analogs to our own Kuiper Belt and they are detected around at least 17% of nearby Sun-like stars. The morphology and dynamics of a disk encode information about its history, as well as that of any exoplanets within the system. We used ALMA to obtain 1.3 mm observations of the debris disk around the nearby F5V star HD 170773. We image the face-on ring and determine its fundamental parameters by forward-modeling the interferometric visibilities through a Markov Chain Monte Carlo approach. Using a symmetric Gaussian surface density profile, we find a 71 $\pm$ 4 au wide belt with a radius of 193$^{+2}_{-3}$ au, a relatively large radius compared to most other millimeter-resolved belts around late A / early F type stars. This makes HD 170773 part of a group of four disks around A and F stars with radii larger than expected from the recently reported planetesimal belt radius - stellar luminosity relation. Two of these systems are known to host directly imaged giant planets, which may point to a connection between large belts and the presence of long-period giant planets. We also set upper limits on the presence of CO and CN gas in the system, which imply that the exocomets that constitute this belt have CO and HCN ice mass fractions of <77% and <3%, respectively, consistent with Solar System comets and other exocometary belts.

Debris disks are exoplanetary systems containing planets, minor bodies (such as asteroids and comets) and debris dust. Unseen planets are presumed to perturb the minor bodies into crossing orbits, generating small dust grains that are detected via remote sensing. Debris disks have been discovered around main sequence stars of a variety of ages (from 10 Myr to several Gyr) and stellar spectral types (from early A-type to M-type stars). As a result, they serve as excellent laboratories for understanding whether the architecture and the evolution of our Solar System is common or rare. This white paper addresses two outstanding questions in debris disk science: (1) Are debris disk minor bodies similar to asteroids and comets in our Solar System? (2) Do planets separate circumstellar material into distinct reservoirs and/or mix material during planet migration? We anticipate that SOFIA/HIRMES, JWST, and WFIRST/CGI will greatly improve our understanding of debris disk composition, enabling the astronomical community to answer these questions. However, we note that despite their observational power, these facilities will not provide large numbers of detections or detailed characterization of cold ices and silicates in the Trans Neptunian zone. Origins Space Telescope is needed to revolutionize our understanding of the bulk composition and mixing in exoplanetary systems.

The vertical distribution of dust in debris disks is sensitive to the number and size of large planetesimals dynamically stirring the disk, and is therefore well-suited for constraining the prevalence of otherwise unobservable Uranus and Neptune analogs. Information regarding stirring bodies has previously been inferred from infrared and optical observations of debris disk vertical structure, but theoretical works predict that the small particles traced by short-wavelength observations will be `puffed up' by radiation pressure, yielding only upper limits. The large grains that dominate the disk emission at millimeter wavelengths are much less sensitive to the effects of stellar radiation or stellar winds, and therefore trace the underlying mass distribution more directly. Here we present ALMA 1.3 mm dust continuum observations of the debris disk around the nearby M star AU Mic. The 3 au spatial resolution of the observations, combined with the favorable edge-on geometry of the system, allows us to measure the vertical thickness of the disk. We report a scale height-to-radius aspect ratio of $h = 0.031_{-0.004}^{+0.005}$ between radii of $\sim 23$ au and $\sim 41$ au. Comparing this aspect ratio to a theoretical model of size-dependent velocity distributions in the collisional cascade, we find that the perturbing bodies embedded in the local disk must be larger than about 400 km, and the largest perturbing body must be smaller than roughly 1.8 M$_\odot$. These measurements rule out the presence of a gas giant or Neptune analog near the $\sim 40$ au outer edge of the debris ring, but are suggestive of large planetesimals or an Earth-sized planet stirring the dust distribution.

Evidence for exocomets, icy bodies in extrasolar planetary systems, has rapidly increased over the past decade. Volatiles are detected through the gas that exocomets release as they collide and grind down within their natal belts, or as they sublimate once scattered inwards to the regions closest to their host star. Most detections are in young, 10 to a few 100 Myr-old systems that are undergoing the final stages of terrestrial planet formation. This opens the exciting possibility to study exocomets at the epoch of volatile delivery to the inner regions of planetary systems. Detection of molecular and atomic gas in exocometary belts allows us to estimate molecular ice abundances and overall elemental abundances, enabling comparison with the Solar Nebula and Solar System comets. At the same time, observing star-grazing exocomets transiting in front of their star (for planetary systems viewed edge-on) and exozodiacal dust in the systems' innermost regions gives unique dynamical insights into the inward scattering process producing delivery to inner rocky planets. The rapid advances of this budding subfield of exoplanetary science will continue in the short term with the upcoming JWST, WFIRST and PLATO missions. In the longer term, the priority should be to explore the full composition of exocomets, including species crucial for delivery and later prebiotic synthesis. Doing so around an increasingly large population of exoplanetary systems is equally important, to enable comparative studies of young exocomets at the epoch of volatile delivery. We identify the proposed LUVOIR and Origins flagship missions as the most promising for a large-scale exploration of exocometary gas, a crucial component of the chemical heritage of young exo-Earths.

Many white dwarf stars show signs of having accreted smaller bodies, implying that they may host planetary systems. A small number of these systems contain gaseous debris discs, visible through emission lines. We report a stable 123.4min periodic variation in the strength and shape of the CaII emission line profiles originating from the debris disc around the white dwarf SDSSJ122859.93+104032.9. We interpret this short-period signal as the signature of a solid body held together by its internal strength.

Thousands of exoplanets have been found with many widely different from the ones in our own system. Despite the success, systems with planets in wide orbits analogous to those of Jupiter and Saturn, in the critical first several hundred million years of evolution, are virtually unexplored. Where are the low-mass planets that are hidden from our exoplanet detection techniques? Is our Solar System's planetary architecture unique? High-fidelity debris disk images offer an effective method to answer these questions. We can use them to study the formation and evolution of low-mass planets from youth to the age of the Solar System, providing snapshots of the complex processes and valuable insights into the formation and migration history of giant planets at wide orbits. This white paper focuses on resolving debris structures in thermal emission that is applicable to a large unbiased sample. We summarize the properties of the known debris disks and assess the feasibility of resolving them within our current and future infrared and millimeter facilities by adopting uniform criteria. JWST and the 9-m Origins Space Telescope are the most promising missions in the coming decades to resolve almost half of the known disks at high fidelity. Resolved debris structures at multiple wavelengths and at all stages of evolution would reveal the properties of unseen planet populations, enabling a unique demographic study of overall planet formation and evolution.

If the Solar system had a history of planet migration, the signature of that migration may be imprinted on the populations of asteroids and comets that were scattered in the planets' wake. Here, we consider the dynamical and collisional evolution of inner Solar system asteroids which join the Oort cloud. We compare the Oort cloud asteroid populations produced by migration scenarios based on the `Nice' and `Grand Tack' scenarios, as well as a null hypothesis where the planets have not migrated, to the detection of one such object, C/2014 S3 (PANSTARRS). Our simulations find that the discovery of C/2014 S3 (PANSTARRS) only has a greater than one percent chance of occurring if the Oort cloud asteroids evolved on to Oort cloud orbits when the Solar system was not more than about one million years old, as this early transfer to the Oort cloud is necessary to keep the amount of collisional evolution low. We argue this only occurs when a giant (greater than thirty Earth masses) planet orbits at 1 ~ 2 au, and thus our results strongly favour a `Grand Tack'-like migration having occurred early in the Solar system's history.

In this paper we present new ALMA observations towards the proto-planet hosting transitional disc of Herbig Ae/Be star HD 100546. This includes resolved 1.3 mm continuum, $^{13}$CO and the first detection of C$^{18}$O in this disc, which displays azimuthal asymmetry in regions spatially coincident with structures previously identified in HST images related to spiral arms. The lower limit on the mass of the dust disc is calculated to be 9.6x10$^{-4}$M$_\odot$. A firm lower-limit on the total gas mass calculated from optically thin, mid-plane tracing C$^{18}$O (2-1) emission is 0.018M$_\odot$ assuming ISM abundances. These mass estimates provide an estimate of gas-to-dust ratio in the disc of 19, the ratio will increase if C$^{18}$O is relatively under-abundant in the disc compared to CO and H2. Through deprojection and azimuthal averaging of the image plane we detect 1.3 mm continuum emission out to 290+/-10 au,$^{13}$CO to 390+/-10 au and C$^{18}$O to 300+/-10au. We measure a radially increasing millimetre spectral index between wavelengths of 867$\mu$m and 1.3 mm, which shows that grain sizes increase towards the star, with solid particles growing to cm scales in the inner disc.

The inclination distribution of the Kuiper belt provides unique constraints on its origin and dynamical evolution, motivating vertically resolved observations of extrasolar planetesimal belts. We present ALMA observations of millimeter emission in the near edge-on planetesimal belt around $\beta$ Pictoris, finding that the vertical distribution is significantly better described by the sum of two Gaussians compared to a single Gaussian. This indicates that, as for the Kuiper belt, the inclination distribution of $\beta$ Pic's belt is better described by the sum of dynamically hot and cold populations rather than a single component. The hot and cold populations have RMS inclinations of 8.9$^{+0.7}_{-0.5}$ and 1.1$^{+0.5}_{-0.5}$ degrees. We also report that an axisymmetric belt model provides a good fit to new and archival ALMA visibilities, and confirm that the midplane is misaligned with respect to $\beta$ Pic b's orbital plane. However, we find no significant evidence for either the inner disk tilt observed in scattered light and CO emission or the South-West/North-East (SW/NE) asymmetry previously reported for millimeter emission. Finally, we consider the origin of the belt's inclination distribution. Secular perturbations from $\beta$ Pic b are unlikely to provide sufficient dynamical heating to explain the hot population throughout the belt's radial extent, and viscous stirring from large bodies within the belt alone cannot reproduce the two populations observed. This argues for an alternative or additional scenario, such as planetesimals being born with high inclinations, or the presence of a `$\beta$ Pic c' planet, potentially migrating outwards near the belt's inner edge.

Nearly all young stars are initially surrounded by `protoplanetary' discs of gas and dust, and in the case of single stars at least 30\% of these discs go on to form planets. The process of protoplanetary disc formation can result in initial misalignments, where the disc orbital plane is different to the stellar equator in single star systems, or to the binary orbital plane in systems with two stars. A quirk of the dynamics means that initially misaligned `circumbinary' discs -- those that surround two stars -- are predicted to evolve to one of two possible stable configurations, one where the disc and binary orbital planes are coplanar, and one where they are perpendicular (a `polar' configuration). Prior work has found coplanar circumbinary discs, but no polar examples were known until now. Here we report the first discovery of a protoplanetary circumbinary disc in the polar configuration, supporting the predictions that such discs should exist. The disc shows some characteristics that are similar to discs around single stars, and that are attributed to dust growth. Thus, the first stages of planet formation appear able to proceed in polar circumbinary discs.

In the last few years, multiwavelength observations have revealed the ubiquity of gaps/rings in circumstellar discs. Here we report the first ALMA observations of HD 92945 at 0.86 mm, that reveal a gap at about 73$\pm$3 au within a broad disc of planetesimals that extends from 50 to 140 au. We find that the gap is $20^{+10}_{-8}$ au wide. If cleared by a planet in situ, this planet must be less massive than 0.6 $M_\mathrm{Jup}$, or even lower if the gap was cleared by a planet that formed early in the protoplanetary disc and prevented planetesimal formation at that radius. By comparing opposite sides of the disc we also find that the disc could be asymmetric. Motivated by the asymmetry and the fact that planets might be more frequent closer to the star in exoplanetary systems, we show that the gap and asymmetry could be produced by two planets interior to the disc through secular resonances. These planets excite the eccentricity of bodies at specific disc locations, opening radial gaps in the planetesimal distribution. New observations are necessary to confirm if the disc is truly asymmetric, thus favouring the secular resonance model, or if the apparent asymmetry is due to a background galaxy, favouring the in-situ planet scenario. Finally, we also report the non-detection of CO and HCN gas confirming that no primordial gas is present. The CO and HCN non-detections are consistent with the destruction of volatile-rich Solar System-like comets.

The quest for other habitable worlds and the search for life among them are major goals of modern astronomy. One way to make progress towards these goals is to obtain high-quality spectra of a large number of exoplanets over a broad range of wavelengths. While concepts currently investigated in the United States are focused on visible/NIR wavelengths, where the planets are probed in reflected light, a compelling alternative to characterize planetary atmospheres is the mid-infrared waveband (5-20~$\mu$m). Indeed, mid-infrared observations provide key information on the presence of an atmosphere, the surface conditions (e.g., temperature, pressure, habitability), and the atmospheric composition in important species such as H$_2$O, CO$_2$, O$_3$, CH$_4$, and N$_2$O. This information is essential to investigate the potential habitability of exoplanets and to make progress towards the search for life in the universe. Obtaining high-quality mid-infrared spectra of exoplanets from the ground is however extremely challenging due to the overwhelming brightness and turbulence of Earth's atmosphere. In this paper, we present a concept of space-based mid-infrared interferometer that can tackle this observing challenge and discuss the main technological developments required to launch such a sophisticated instrument.

The properties of dust in the interstellar medium (ISM) nearest the Sun are poorly understood because the low column densities of dust toward nearby stars induce little photometric reddening, rendering the grains largely undetectable. Stellar polarimetry offers one pathway to deducing the properties of this diffuse material. Here we present multi-wavelength aperture polarimetry measurements of seven bright stars chosen to probe interstellar polarization near the edge of the Local Hot Bubble (LHB) - an amorphous region of relatively low density interstellar gas and dust extending ~70-150 pc from the Sun. The measurements were taken using the HIgh Precision Polarimetric Instrument (HIPPI) on the 3.9-m Anglo-Australian Telescope. HIPPI is an aperture stellar polarimeter with a demonstrated sensitivity of 4.3 parts-per-million (ppm). Of the stars observed two are polarized to a much greater degree than the others; they have a wavelength of maximum polarization ($\lambda_{max}$) of ~550 $\pm$ 20 nm - similar to that of stars beyond the LHB - and we conclude that they are in the wall of the LHB. The remaining five stars have polarizations of ~70 to 160 ppm, of these four have a much bluer $\lambda_{max}$, ~350 $\pm$ 50 nm. Bluer values of $\lambda_{max}$ may indicate grains shocked during the evolution of the Loop I Superbubble. The remaining star, HD 4150 is not well fit by a Serkowski curve, and may be intrinsically polarized.

Despite being $>10$Myr, there are $\sim$10 debris discs with as much CO gas as in protoplanetary discs. Such discs have been assumed to be "hybrid", i.e., with secondary dust but primordial gas. Here we show that both the dust and gas in such systems could instead be secondary, with the high CO content caused by accumulation of neutral carbon (C$^0$) that shields CO from photodissociating; i.e., these could be "shielded secondary discs". New ALMA observations are presented of HD131835 that detect $\sim 3 \times 10^{-3}$ M$_\oplus$ of C$^0$, the majority 40-200au from the star, in sufficient quantity to shield the previously detected CO. A simple semi-analytic model for the evolution of CO, C and O originating in a volatile-rich planetesimal belt shows how CO shielding becomes important when the viscous evolution is slow (low $\alpha$ parameter) and/or the CO production rate is high. Shielding by C$^0$ may also cause the CO content to reach levels at which CO self-shields, and the gas disc may become massive enough to affect the dust evolution. Application to the HD 131835 observations shows these can be explained if $\alpha \sim 10^{-3}$; an inner cavity in C$^0$ and CO may also mean the system has yet to reach steady state. Application to other debris discs with high CO content finds general agreement for $\alpha=10^{-3}$ to $0.1$. The shielded secondary nature of these gas discs can be tested by searching for C$^0$, as well as CN, N$_2$ and CH$^{+}$, which are also expected to be shielded by C$^0$.

This paper discusses an algorithm for detecting single transits in photometric time-series data. Specifically, we aim to identify asymmetric transits with ingress that is more rapid than egress, as expected for cometary bodies with a significant tail. The algorithm is automated, so can be applied to large samples and only a relatively small number of events need to be manually vetted. We applied this algorithm to all long cadence light curves from the Kepler mission, finding 16 candidate transits with significant asymmetry, 11 of which were found to be artefacts or symmetric transits after manual inspection. Of the 5 remaining events, four are the 0.1% depth events previously identified for KIC 3542116 and 11084727. We identify HD 182952 (KIC 8027456) as a third system showing a potential comet transit. All three stars showing these events have H-R diagram locations consistent with $\sim$100Myr-old open cluster stars, as might be expected given that cometary source regions deplete with age, and giving credence to the comet hypothesis. If these events are part of the same population of events as seen for KIC 8462852, the small increase in detections at 0.1% depth compared to 10% depth suggests that future work should consider whether the distribution is naturally flat, or if comets with symmetric transits in this depth range remain undiscovered. Future searches relying on asymmetry should be more successful if they focus on larger samples and young stars, rather than digging further into the noise.

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.

WD 1145+017 is currently the only white dwarf known to exhibit periodic transits of planetary debris as well as absorption lines from circumstellar gas. We present the first simultaneous fast optical spectrophotometry and broad-band photometry of the system, obtained with the Gran Telescopio Canarias (GTC) and the Liverpool Telescope (LT), respectively. The observations spanned $5.5$ h, somewhat longer than the $4.5$-h orbital period of the debris. Dividing the GTC spectrophotometry into five wavelength bands reveals no significant colour differences, confirming grey transits in the optical. We argue that absorption by an optically thick structure is a plausible alternative explanation for the achromatic nature of the transits that can allow the presence of small-sized ($\sim\mu$m) particles. The longest ($87$ min) and deepest ($50$ per cent attenuation) transit recorded in our data exhibits a complex structure around minimum light that can be well modelled by multiple overlapping dust clouds. The strongest circumstellar absorption line, Fe II $\lambda$5169, significantly weakens during this transit, with its equivalent width reducing from a mean out-of-transit value of $2$ \AA\ to $1$ \AA\ in-transit, supporting spatial correlation between the circumstellar gas and dust. Finally, we made use of the Gaia Data Release 2 and archival photometry to determine the white dwarf parameters. Adopting a helium-dominated atmosphere containing traces of hydrogen and metals, and a reddening $E(B-V)=0.01$ we find $T_\mathrm{eff}=15\,020 \pm 520$ K, $\log g=8.07\pm0.07$, corresponding to $M_\mathrm{WD}=0.63\pm0.05\ \mbox{$\mathrm{M}_{\odot}$}$ and a cooling age of $224\pm30$ Myr.

[abridged] The presence of large amounts of dust in the habitable zones of nearby stars is a significant obstacle for future exo-Earth imaging missions. We executed an N band nulling interferometric survey to determine the typical amount of such exozodiacal dust around a sample of nearby main sequence stars. The majority of our data have been analyzed and we present here an update of our ongoing work. We find seven new N band excesses in addition to the high confidence confirmation of three that were previously known. We find the first detections around Sun-like stars and around stars without previously known circumstellar dust. Our overall detection rate is 23%. The inferred occurrence rate is comparable for early type and Sun-like stars, but decreases from 71% [+11%/-20%] for stars with previously detected mid- to far-infrared excess to 11% [+9%/-4%] for stars without such excess, confirming earlier results at high confidence. For completed observations on individual stars, our sensitivity is five to ten times better than previous results. Assuming a lognormal luminosity function of the dust, we find upper limits on the median dust level around all stars without previously known mid to far infrared excess of 11.5 zodis at 95% confidence level. The corresponding upper limit for Sun-like stars is 16 zodis. An LBTI vetted target list of Sun-like stars for exo-Earth imaging would have a corresponding limit of 7.5 zodis. We provide important new insights into the occurrence rate and typical levels of habitable zone dust around main sequence stars. Exploiting the full range of capabilities of the LBTI provides a critical opportunity for the detailed characterization of a sample of exozodiacal dust disks to understand the origin, distribution, and properties of the dust.

Cold outer debris belts orbit a significant fraction of stars, many of which are planet-hosts. Radiative forces from the star lead to dust particles leaving the outer belts and spiralling inwards under Poynting-Robertson drag. We present an empirical model fitted to N-body simulations that allows the fate of these dust particles when they encounter a planet to be rapidly calculated. High mass planets eject most particles, whilst dust passes low mass planets relatively unperturbed. Close-in, high mass planets (hot Jupiters) are best at accreting dust. The model predicts the accretion rate of dust onto planets interior to debris belts, with mass accretions rates of up to hundreds of kilograms per second predicted for hot Jupiters interior to outer debris belts, when collisional evolution is also taken into account. The model can be used to infer the presence and likely masses of as yet undetected planets in systems with outer belts. The non-detection of warm dust with the Large Binocular Telescope Interferometer (LBTI) around Vega could be explained by the presence of a single Saturn mass planet, or a chain of lower mass planets. Similarly, the detection of warm dust in such systems implies the absence of planets above a quantifiable level, which can be lower than similar limits from direct imaging. The level of dust detected with LBTI around beta Leo can be used to rule out the presence of planets more massive than a few Saturn masses outside of ~5au.

The TRAPPIST-1 system is unique in that it has a chain of seven terrestrial Earth-like planets located close to or in its habitable zone. In this paper, we study the effect of potential cometary impacts on the TRAPPIST-1 planets and how they would affect the primordial atmospheres of these planets. We consider both atmospheric mass loss and volatile delivery with a view to assessing whether any sort of life has a chance to develop. We ran N-body simulations to investigate the orbital evolution of potential impacting comets, to determine which planets are more likely to be impacted and the distributions of impact velocities. We consider three scenarios that could potentially throw comets into the inner region (i.e within 0.1au where the seven planets are located) from an (as yet undetected) outer belt similar to the Kuiper belt or an Oort cloud: Planet scattering, the Kozai-Lidov mechanism and Galactic tides. For the different scenarios, we quantify, for each planet, how much atmospheric mass is lost and what mass of volatiles can be delivered over the age of the system depending on the mass scattered out of the outer belt. We find that the resulting high velocity impacts can easily destroy the primordial atmospheres of all seven planets, even if the mass scattered from the outer belt is as low as that of the Kuiper belt. However, we find that the atmospheres of the outermost planets f, g and h can also easily be replenished with cometary volatiles (e.g. $\sim$ an Earth ocean mass of water could be delivered). These scenarios would thus imply that the atmospheres of these outermost planets could be more massive than those of the innermost planets, and have volatiles-enriched composition.

While detecting low mass exoplanets at tens of au is beyond current instrumentation, debris discs provide a unique opportunity to study the outer regions of planetary systems. Here we report new ALMA observations of the 80-200 Myr old Solar analogue HD 107146 that reveal the radial structure of its exo-Kuiper belt at wavelengths of 1.1 and 0.86 mm. We find that the planetesimal disc is broad, extending from 40 to 140 au, and it is characterised by a circular gap extending from 60 to 100 au in which the continuum emission drops by about 50%. We also report the non-detection of the CO J=3-2 emission line, confirming that there is not enough gas to affect the dust distribution. To date, HD 107146 is the only gas-poor system showing multiple rings in the distribution of millimeter sized particles. These rings suggest a similar distribution of the planetesimals producing small dust grains that could be explained invoking the presence of one or more perturbing planets. Because the disk appears axisymmetric, such planets should be on circular orbits. By comparing N-body simulations with the observed visibilities we find that to explain the radial extent and depth of the gap, it would be required the presence of multiple low mass planets or a single planet that migrated through the disc. Interior to HD 107146's exo-Kuiper belt we find extended emission with a peak at ~20 au and consistent with the inner warm belt that was previously predicted based on 22$\mu$m excess as in many other systems. This warm belt is the first to be imaged, although unexpectedly suggesting that it is asymmetric. This could be due to a large belt eccentricity or due to clumpy structure produced by resonant trapping with an additional inner planet.

The star KIC 8462852 (Boyajian's Star) displays both fast dips of up to 20% on time scales of days, plus long-term secular fading by up to 19% on time scales from a year to a century. We report on CCD photometry of KIC 8462852 from 2015.75 to 2018.18, with 19,176 images making for 1,866 nightly magnitudes in BVRI. Our light curves show a continuing secular decline (by 0.023 +- 0.003 mags in the B-band) with three superposed dips with duration 120-180 days. This demonstrates that there is a continuum of dip durations from a day to a century, so that the secular fading is seen to be by the same physical mechanism as the short-duration Kepler dips. The BVRI light curves all have the same shape, with the slopes and amplitudes for VRI being systematically smaller than in the B-band by factors of 0.77 +- 0.05, 0.50 +- 0.05, and 0.31 +- 0.05. We rule out any hypothesis involving occultation of the primary star by any star, planet, solid body, or optically thick cloud. But these ratios are the same as that expected for ordinary extinction by dust clouds. This chromatic extinction implies dust particle sizes going down to ~0.1 micron, suggesting that this dust will be rapidly blown away by stellar radiation pressure, so the dust clouds must have formed within months. The modern infrared observations were taken at a time when there was at least 12.4% +- 1.3% dust coverage (as part of the secular dimming), and this is consistent with dimming originating in circumstellar dust.

HD141569 is a pre-main sequence star with a disc uniquely placed between protoplanetary and debris discs, similar to the older `hybrid' type discs. This work aims to place the mass and spatial structure of the disc midplane in the context of the debris, hybrid and protoplanetary discs. We observed HD~141569 with ALMA in 1.3~mm continuum and $^{13}$CO (2-1). This is the first detection and image of the optically thin gas emission from the midplane of this disc. In continuum emission, we detect a combination of an unresolved central peak and a ring of millimetre emission at 220$\pm$10~au, slightly interior to one of the rings discovered in scattered light. The minimum dust mass of the ring is 0.13$\pm$0.02~M$_{\oplus}$ while the unresolved millimetre peak at the stellar location is predominantly thermal emission due to a minimum of 1.2$\pm$0.2~M$_{\oplus}$ of dust. $^{13}$CO is distributed asymmetrically around the stellar position with a peak at 1.1" distance and a P.A. of -33$^\circ$. The gas is detected as far as 220$\pm$10~au, a radial separation the same as that of the mm ring. Assuming optically thin emission and standard ISM abundances, we use our $^{13}$CO data to derive the gas mass in the disc of (6.0$\pm$0.9) $\times 10^{-4}~$M$_\odot$. Comparison to published $^{12}$CO data shows that $^{12}$CO is optically thick, explaining why estimates based on $^{12}$CO underestimated the gas mass.

This paper shows how the dips and secular dimming in the KIC8462852 light curve can originate in circumstellar material distributed around a single elliptical orbit (e.g., exocomets). The expected thermal emission and wavelength dependent dimming is derived for different orbital parameters and geometries, including dust that is optically thick to stellar radiation, and for a size distribution of dust with realistic optical properties. We first consider dust distributed evenly around the orbit, then show how to derive its uneven distribution from the optical light curve and to predict light curves at different wavelengths. The fractional luminosity of an even distribution is approximately the level of dimming times stellar radius divided by distance from the star at transit. Non-detection of dust thermal emission for KIC8462852 thus provides a lower limit on the transit distance to complement the 0.6au upper limit imposed by 0.4day dips. Unless the dust distribution is optically thick, the putative 16% century-long secular dimming must have disappeared before the WISE 12micron measurement in 2010, and subsequent 4.5micron observations require transits at >0.05au. However, self-absorption of thermal emission removes these constraints for opaque dust distributions. The passage of dust clumps through pericentre is predicted to cause infrared brightening lasting 10s of days and dimming during transit, such that total flux received decreases at wavelengths <5micron, but increases to potentially detectable levels at longer wavelengths. We suggest that lower dimming levels than seen for KIC8462852 are more common in the Galactic population and may be detected in future transit surveys.

Exocomets scattered by planets have been invoked to explain observations in multiple contexts, including the frequently found near- and mid-infrared excess around nearby stars arising from exozodiacal dust. Here we investigate how the process of inward scattering of comets originating in an outer belt, is affected by the architecture of a planetary system, to determine whether this could lead to observable exozodi levels or deliver volatiles to inner planets. Using N-body simulations, we model systems with different planet mass and orbital spacing distributions in the 1-50 AU region. We find that tightly packed ($\Delta a_\mathrm{p}<20 R_\mathrm{Hill,m}$) low mass planets are the most efficient at delivering material to exozodi regions (5-7% of scattered exocomets end up within 0.5 AU at some point), although the exozodi levels do not vary by more than a factor of ~7 for the architectures studied here. We suggest that emission from scattered dusty material in between the planets could provide a potential test for this delivery mechanism. We show that the surface density of scattered material can vary by two orders of magnitude (being highest for systems of low mass planets with medium spacing), whilst the exozodi delivery rate stays roughly constant, and that future instruments such as JWST could detect it. In fact for $\eta$ Corvi, the current Herschel upper limit rules out the scattering scenario by a chain of $\lesssim$30 M$_\oplus$ planets. Finally, we show that exocomets could be efficient at delivering cometary material to inner planets (0.1-1% of scattered comets are accreted per inner planet). Overall, the best systems at delivering comets to inner planets are the ones that have low mass outer planets and medium spacing ($\sim20 R_\mathrm{Hill,m}$).

Several tens of white dwarfs are known to host circumstellar discs of dusty debris, thought to arise from the tidal disruption of rocky bodies originating in the star's remnant planetary system. This paper investigates the evolution of such discs if they are very massive, as may be the case if their progenitor was a terrestrial planet, moon, or dwarf planet. Assuming the discs are physically thin and flat, like Saturn's rings, their evolution is governed by Poynting-Robertson drag or viscous spreading, where the disc's effective viscosity is due to self-gravity wakes. For discs with masses >10^26 g, located in the outer parts of the tidal disruption zone, viscous spreading dominates the evolution, and mass is transported both in- and outwards. When outwards-spreading material flows beyond the Roche limit, it coagulates into new (minor) planets in a process analogous to the ongoing formation of moonlets at the outer edge of Saturn's rings. The newly formed bodies migrate outwards by exchanging angular momentum with the disc and coalesce into larger objects through mutual collisions. Eventually, the disc's Roche-limit overflow recycles tens of percent of the original disc mass; most ends up in a single large body near 2:1 mean-motion resonance with the disc's outer edge. Hence, the recycling of a tidally disrupted super-Earth, for example, could yield an Earth-mass planet on a ~10-h orbit, located in the habitable zone for 2-to-10-Gyr-old white dwarfs. The recycling process also creates a population of smaller bodies just outside the Roche limit, which may explain the minor planets recently postulated to orbit WD 1145+017.

Debris disks are the dust disks found around ~20% of nearby main sequence stars in far-IR surveys. They can be considered as descendants of protoplanetary disks or components of planetary systems, providing valuable information on circumstellar disk evolution and the outcome of planet formation. The debris disk population can be explained by the steady collisional erosion of planetesimal belts; population models constrain where (10-100au) and in what quantity (>1Mearth) planetesimals (>10km in size) typically form in protoplanetary disks. Gas is now seen long into the debris disk phase. Some of this is secondary implying planetesimals have a Solar System comet-like composition, but some systems may retain primordial gas. Ongoing planet formation processes are invoked for some debris disks, such as the continued growth of dwarf planets in an unstirred disk, or the growth of terrestrial planets through giant impacts. Planets imprint structure on debris disks in many ways; images of gaps, clumps, warps, eccentricities and other disk asymmetries, are readily explained by planets at >>5au. Hot dust in the region planets are commonly found (<5au) is seen for a growing number of stars. This dust usually originates in an outer belt (e.g., from exocomets), although an asteroid belt or recent collision is sometimes inferred.

Resolved observations of millimetre-sized dust, tracing larger planetesimals, have pinpointed the location of 26 Edgeworth-Kuiper belt analogs. We report that a belt's distance $R$ to its host star correlates with the star's luminosity $L_{\star}$, following $R\propto L^{0.19}_{\star}$ with a low intrinsic scatter of $\sim$17%. Remarkably, our Edgeworth-Kuiper belt in the Solar System and the two CO snow lines imaged in protoplanetary disks lie close to this $R$-$L_{\star}$ relation, suggestive of an intrinsic relationship between protoplanetary disk structures and belt locations. To test the effect of bias on the relation, we use a Monte Carlo approach and simulate uncorrelated model populations of belts. We find that observational bias could produce the slope and intercept of the $R$-$L_{\star}$ relation, but is unable to reproduce its low scatter. We then repeat the simulation taking into account the collisional evolution of belts, following the steady state model that fits the belt population as observed through infrared excesses. This significantly improves the fit by lowering the scatter of the simulated $R$-$L_{\star}$ relation; however, this scatter remains only marginally consistent with the one observed. The inability of observational bias and collisional evolution alone to reproduce the tight relationship between belt radius and stellar luminosity could indicate that planetesimal belts form at preferential locations within protoplanetary disks. The similar trend for CO snow line locations would then indicate that the formation of planetesimals and/or planets in the outer regions of planetary systems is linked to the volatility of their building blocks, as postulated by planet formation models.

The HOSTS (Hunt for Observable Signatures of Terrestrial Systems) survey searches for dust near the habitable zones (HZs) around nearby, bright main sequence stars. We use nulling interferometry in N band to suppress the bright stellar light and to probe for low levels of HZ dust around the 30 stars observed so far. Our overall detection rate is 18%, including four new detections, among which are the first three around Sun-like stars and the first two around stars without any previously known circumstellar dust. The inferred occurrence rates are comparable for early type and Sun-like stars, but decrease from 60 (+16/-21)% for stars with previously detected cold dust to 8 (+10/-3)% for stars without such excess, confirming earlier results at higher sensitivity. For completed observations on individual stars, our sensitivity is five to ten times better than previous results. Assuming a lognormal excess luminosity function, we put upper limits on the median HZ dust level of 13 zodis (95% confidence) for a sample of stars without cold dust and of 26 zodis when focussing on Sun-like stars without cold dust. However, our data suggest that a more complex luminosity function may be more appropriate. For stars without detectable LBTI excess, our upper limits are almost reduced by a factor of two, demonstrating the strength of LBTI target vetting for future exo-Earth imaging missions. Our statistics are so far limited and extending the survey is critical to inform the design of future exo-Earth imaging surveys.

We present the results of a Herschel survey of 21 late-type stars that host planets discovered by the radial velocity technique. The aims were to discover new disks in these systems and to search for any correlation between planet presence and disk properties. In addition to the known disk around GJ 581, we report the discovery of two new disks, in the GJ 433 and GJ 649 systems. Our sample therefore yields a disk detection rate of 14%, higher than the detection rate of 1.2% among our control sample of DEBRIS M-type stars with 98% confidence. Further analysis however shows that the disk sensitivity in the control sample is about a factor of two lower in fractional luminosity than for our survey, lowering the significance of any correlation between planet presence and disk brightness below 98%. In terms of their specific architectures, the disk around GJ 433 lies at a radius somewhere between 1 and 30au. The disk around GJ 649 lies somewhere between 6 and 30au, but is marginally resolved and appears more consistent with an edge-on inclination. In both cases the disks probably lie well beyond where the known planets reside (0.06-1.1au), but the lack of radial velocity sensitivity at larger separations allows for unseen Saturn-mass planets to orbit out to $\sim$5au, and more massive planets beyond 5au. The layout of these M-type systems appears similar to Sun-like star + disk systems with low-mass planets.

This paper presents a study of circumstellar debris around Sun-like stars using data from the Herschel DEBRIS Key Programme. DEBRIS is an unbiased survey comprising the nearest ~90 stars of each spectral type A-M. Analysis of the 275 F-K stars shows that excess emission from a debris disc was detected around 47 stars, giving a detection rate of 17.1+2.6-2.3 per cent, with lower rates for later spectral types. For each target a blackbody spectrum was fitted to the dust emission to determine its fractional luminosity and temperature. The derived under- lying distribution of fractional luminosity versus blackbody radius in the population showed that most detected discs are concentrated at f ~ 10^-5 and at temperatures corresponding to blackbody radii 7-40 AU, which scales to ~40 AU for realistic dust properties (similar to the current Kuiper belt). Two outlying populations are also evident; five stars have exceptionally bright emission ( f > 5x10^-5), and one has unusually hot dust < 4 AU. The excess emission distributions at all wavelengths were fitted with a steady-state evolution model, showing these are compatible with all stars being born with a narrow belt that then undergoes collisional grinding. However, the model cannot explain the hot dust systems - likely originating in transient events - and bright emission systems - arising potentially from atypically massive discs or recent stirring. The emission from the present-day Kuiper belt is predicted to be close to the median of the population, suggesting that half of stars have either depleted their Kuiper belts (similar to the Solar System), or had a lower planetesimal formation efficiency.

HR8799 is a benchmark system for direct imaging studies. It hosts two debris belts, which lie internally and externally to four giant planets. This paper considers how the four known planets and a possible fifth planet, interact with the external population of debris through N-body simulations. We find that when only the known planets are included, the inner edge of the outer belt predicted by our simulations is much closer to the outermost planet than recent ALMA observations suggest. We subsequently include a fifth planet in our simulations with a range of masses and semi-major axes, which is external to the outermost known planet. We find that a fifth planet with a mass and semi-major axis of 0.1$\mathrm{M_J}$ and 138au predicts an outer belt that agrees well with ALMA observations, whilst remaining stable for the lifetime of HR8799 and lying below current direct imaging detection thresholds. We also consider whether inward scattering of material from the outer belt can input a significant amount of mass into the inner belt. We find that for the current age of HR8799, only $\sim$1\% of the mass loss rate of the inner disk can be replenished by inward scattering. However we find that the higher rate of inward scattering during the first $\sim$10Myr of HR8799 would be expected to cause warm dust emission at a level similar to that currently observed, which may provide an explanation for such bright emission in other systems at $\sim10$Myr ages.

The young A0V star HR 4796A is host to a bright and narrow ring of dust, thought to originate in collisions between planetesimals within a belt analogous to the Solar System's Edgeworth-Kuiper belt. Here we present high spatial resolution 880$\mu$m continuum images from the Atacama Large Millimeter Array. The 80au radius dust ring is resolved radially with a characteristic width of 10au, consistent with the narrow profile seen in scattered light. Our modelling consistently finds that the disk is also vertically resolved with a similar extent. However, this extent is less than the beam size, and a disk that is dynamically very cold (i.e. vertically thin) provides a better theoretical explanation for the narrow scattered light profile, so we remain cautious about this conclusion. We do not detect $^{12}$CO J=3-2 emission, concluding that unless the disk is dynamically cold the CO+CO$_2$ ice content of the planetesimals is of order a few percent or less. We consider the range of semi-major axes and masses of an interior planet supposed to cause the ring's eccentricity, finding that such a planet should be more massive than Neptune and orbit beyond 40au. Independent of our ALMA observations, we note a conflict between mid-IR pericenter-glow and scattered light imaging interpretations, concluding that models where the spatial dust density and grain size vary around the ring should be explored.

The exocometary origin of CO gas has been confirmed in several extrasolar Kuiper belts, with CO ice abundances consistent with Solar System comets. We here present a molecular survey of the $\beta$ Pictoris belt with the Submillimeter Array (SMA), reporting upper limits for CN, HCN, HCO$^+$, N$_2$H$^+$ and H$_2$CO, as well as for H$_2$S, CH$_3$OH, SiO and DCN from archival ALMA data. Non-detections can be attributed to rapid molecular photodissociation due to the A-star's strong UV flux. CN is the longest-lasting and most easily detectable molecule after CO in this environment. We update our NLTE excitation model to include UV fluorescence, finding it plays a key role in CO and CN excitation, and use it to turn the SMA CN/CO flux ratio constraint into an upper limit of $<2.5$% on the HCN/(CO+CO$_2$) ratio of outgassing rates. This value is consistent with, but at the low end of, the broad range observed in Solar System comets. If sublimation dominates outgassing, then this low value may be caused by decreased outgassing for the less volatile molecule HCN compared to CO. If instead UV photodesorption or collisional vaporization of unbound grains dominates outgassing, then this low ratio of rates would imply a low ice abundance ratio, which would in turn indicate a variation in cometary cyanide abundances across planetary systems. To conclude, we make predictions for future molecular surveys and show that CN and HCN should be readily detectable with ALMA around $\beta$ Pictoris for Solar-System-like exocometary compositions.

We present a photometric detection of the first brightness dips of the unique variable star KIC 8462852 since the end of the Kepler space mission in 2013 May. Our regular photometric surveillance started in October 2015, and a sequence of dipping began in 2017 May continuing on through the end of 2017, when the star was no longer visible from Earth. We distinguish four main 1-2.5% dips, named "Elsie," "Celeste," "Skara Brae," and "Angkor", which persist on timescales from several days to weeks. Our main results so far are: (i) there are no apparent changes of the stellar spectrum or polarization during the dips; (ii) the multiband photometry of the dips shows differential reddening favoring non-grey extinction. Therefore, our data are inconsistent with dip models that invoke optically thick material, but rather they are in-line with predictions for an occulter consisting primarily of ordinary dust, where much of the material must be optically thin with a size scale <<1um, and may also be consistent with models invoking variations intrinsic to the stellar photosphere. Notably, our data do not place constraints on the color of the longer-term "secular" dimming, which may be caused by independent processes, or probe different regimes of a single process.

White dwarf WD 1145+017 is orbited by several clouds of dust, possibly emanating from actively disintegrating bodies. These dust clouds reveal themselves through deep, broad, and evolving transits in the star's light curve. Here, we report two epochs of multi-wavelength photometric observations of WD 1145+017, including several filters in the optical, K$_\mathrm{s}$ and 4.5 $\mu$m bands in 2016 and 2017. The observed transit depths are different at these wavelengths. However, after correcting for excess dust emission at K$_\mathrm{s}$ and 4.5 $\mu$m, we find the transit depths for the white dwarf itself are the same at all wavelengths, at least to within the observational uncertainties of $\sim$5%-10%. From this surprising result, and under the assumption of low optical depth dust clouds, we conclude that there is a deficit of small particles (with radii $s \lesssim$ 1.5 $\mu$m) in the transiting material. We propose a model wherein only large particles can survive the high equilibrium temperature environment corresponding to 4.5 hr orbital periods around WD 1145+017, while small particles sublimate rapidly. In addition, we evaluate dust models that are permitted by our measurements of infrared emission.

According to the current paradigm of circumstellar disk evolution, gas-rich primordial disks evolve into gas-poor debris disks compose of second-generation dust. To explore the transition between these phases, we searched for $^{12}$CO, $^{13}$CO, and C$^{18}$O emission in seven dust-rich debris disks around young A-type stars, using ALMA in Band 6. We discovered molecular gas in three debris disks. In all these disks, the $^{12}$CO line was optically thick, highlighting the importance of less abundant molecules in reliable mass estimates. Supplementing our target list by literature data, we compiled a volume-limited sample of dust-rich debris disks around young A-type stars within 150 pc. We obtained a CO detection rate of 11/16 above a $^{12}$CO J=2$-$1 line luminosity threshold of $\sim 1.4 \times 10 ^4$ Jykms$^{-1}$pc$^2$ in the sample. This high incidence implies that the presence of CO gas in bright debris disks around young A-type stars is likely more the rule than the exception. Interestingly, dust-rich debris disks around young FG-type stars exhibit, with the same detectability threshold as for A-type stars, significantly lower gas incidence. While the transition from protoplanetary to debris phase is associated with a drop of dust content, our results exhibit a large spread in the CO mass in our debris sample, with peak values comparable to those in protoplanetary Herbig Ae disks. In the particularly CO-rich debris systems the gas may have primordial origin, characteristic of a hybrid disk.

The so-called "dipper" stars host circumstellar disks and have optical and infrared light curves that exhibit quasi-periodic or aperiodic dimming events consistent with extinction by transiting dusty structures orbiting in the inner disk. Most of the proposed mechanisms explaining the dips---i.e., occulting disk warps, vortices, and forming planetesimals---assume nearly edge-on viewing geometries. However, our analysis of the three known dippers with publicly available resolved sub-mm data reveals disks with a range of inclinations, most notably the face-on transition disk J1604-2130 (EPIC 204638512). This suggests that nearly edge-on viewing geometries are not a defining characteristic of the dippers and that additional models should be explored. If confirmed by further observations of more dippers, this would point to inner disk processes that regularly produce dusty structures far above the outer disk midplane in regions relevant to planet formation.

To test alternative hypotheses for the behavior of KIC 8462852, we obtained measurements of the star over a wide wavelength range from the UV to the mid-infrared from October 2015 through December 2016, using Swift, Spitzer and at AstroLAB IRIS. The star faded in a manner similar to the long-term fading seen in Kepler data about 1400 days previously. The dimming rate for the entire period reported is 22.1 +\- 9.7 milli-mag/yr in the Swift wavebands, with amounts of 21.0 +\- 4.5 mmag in the groundbased B measurements, 14.0 +\- 4.5 mmag in V, and 13.0 +\- 4.5 in R, and a rate of 5.0 +\- 1.2 mmag/yr averaged over the two warm Spitzer bands. Although the dimming is small, it is seen at >= 3 sigma by three different observatories operating from the UV to the IR. The presence of long-term secular dimming means that previous SED models of the star based on photometric measurements taken years apart may not be accurate. We find that stellar models with T_{eff} = 7000 - 7100 K and A_V ~ 0.73 best fit the Swift data from UV to optical. These models also show no excess in the near-simultaneous Spitzer photometry at 3.6 and 4.5 microns, although a longer wavelength excess from a substantial debris disk is still possible (e.g., as around Fomalhaut). The wavelength dependence of the fading favors a relatively neutral color (i.e., R_V >= 5, but not flat across all the bands) compared with the extinction law for the general ISM (R_V = 3.1), suggesting that the dimming arises from circumstellar material.

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

Recent ALMA observations present mounting evidence for the presence of exocometary gas released within Kuiper belt analogues around nearby main sequence stars. This represents a unique opportunity to study their ice reservoir at the younger ages when volatile delivery to planets is most likely to occur. We here present the detection of CO J=2-1 emission co-located with dust emission from the cometary belt in the 440 Myr-old Fomalhaut system. Through spectro-spatial filtering, we achieve a 5.4$\sigma$ detection and determine that the ring's sky-projected rotation axis matches that of the star. The CO mass derived ($0.65-42 \times10^{-7}$ M$_{\oplus}$) is the lowest of any circumstellar disk detected to date, and must be of exocometary origin. Using a steady state model, we estimate the CO+CO$_2$ mass fraction of exocomets around Fomalhaut to be between 4.6-76%, consistent with Solar System comets and the two other belts known to host exocometary gas. This is the first indication of a similarity in cometary compositions across planetary systems that may be linked to their formation scenario and is consistent with direct ISM inheritance. In addition, we find tentative evidence that $(49\pm 27)$% of the detected flux originates from a region near the eccentric belt's pericentre. If confirmed, the latter may be explained through a recent impact event or CO pericentre glow due to exocometary release within a steady state collisional cascade. In the latter scenario, we show how the azimuthal dependence of the CO release rate leads to asymmetries in gas observations of eccentric exocometary belts.

We present ALMA mosaic observations at 1.3 mm (223 GHz) of the Fomalhaut system with a sensitivity of 14 $\mu$Jy/beam. These observations provide the first millimeter map of the continuum dust emission from the complete outer debris disk with uniform sensitivity, enabling the first conclusive detection of apocenter glow. We adopt a MCMC modeling approach that accounts for the eccentric orbital parameters of a collection of particles within the disk. The outer belt is radially confined with an inner edge of $136.3\pm0.9$ AU and width of $13.5\pm1.8$ AU. We determine a best-fit eccentricity of $0.12\pm0.01$. Assuming a size distribution power law index of $q=3.46\pm 0.09$, we constrain the dust absorptivity power law index $\beta$ to be $0.9<\beta<1.5$. The geometry of the disk is robustly constrained with inclination $65.\!\!^\circ6\pm0.\!\!^\circ3$, position angle $337.\!\!^\circ9\pm0.\!\!^\circ3$, and argument of periastron $22.\!\!^\circ5\pm4.\!\!^\circ3$. Our observations do not confirm any of the azimuthal features found in previous imaging studies of the disk with HST, SCUBA, and ALMA. However, we cannot rule out structures $\leq10$ AU in size or which only affect smaller grains. The central star is clearly detected with a flux density of $0.75\pm0.02$ mJy, significantly lower than predicted by current photospheric models. We discuss the implications of these observations for the directly imaged Fomalhaut b and the inner dust belt detected at infrared wavelengths.

A decade of surveys has hinted at a possible higher occurrence rate of debris discs in systems hosting low mass planets. This could be due to common favourable forming conditions for rocky planets close in and planetesimals at large radii. In this paper we present the first resolved millimetre study of the debris disc in the 4.6 Gyr old multiplanet system 61 Vir, combining ALMA and JCMT data at 0.86 mm. We fit the data using a parametric disc model, finding that the disc of planetesimals extends from 30 AU to at least 150 AU, with a surface density distribution of millimetre sized grains with a power law slope of 0.1$^{+1.1}_{-0.8}$. We also present a numerical collisional model that can predict the evolution of the surface density of millimetre grains for a given primordial disc, finding that it does not necessarily have the same radial profile as the total mass surface density (as previous studies suggested for the optical depth), with the former being flatter. Finally, we find that if the planetesimal disc was stirred at 150 AU by an additional unseen planet, that planet should be more massive than 10 M$_{\oplus}$ and lie between 10-20 AU. Lower planet masses and semi-major axes down to 4 AU are possible for eccentricities $\gg$ 0.1.

This paper uses observations of dusty debris discs, including a growing number of gas detections in these systems, to test our understanding of the origin and evolution of this gaseous component. It is assumed that all debris discs with icy planetesimals create second generation CO, C and O gas at some level, and the aim of this paper is to predict that level and assess its observability. We present a new semi-analytical equivalent of the numerical model of Kral et al. (2016) allowing application to large numbers of systems.That model assumes CO is produced from volatile-rich solid bodies at a rate that can be predicted from the debris disc's fractional luminosity. CO photodissociates rapidly into C and O that then evolve by viscous spreading. This model provides a good qualitative explanation of all current observations, with a few exceptional systems that likely have primordial gas. The radial location of the debris and stellar luminosity explain some non-detections, e.g. close-in debris (like HD 172555) is too warm to retain CO, while high stellar luminosities (like $\eta$ Tel) result in short CO lifetimes. We list the most promising targets for gas detections, predicting $>15$ CO detections and $>30$ CI detections with ALMA, and tens of CII and OI detections with future far-IR missions. We find that CO, CI, CII and OI gas should be modelled in non-LTE for most stars, and that CO, CI and OI lines will be optically thick for the most gas-rich systems. Finally, we find that radiation pressure, which can blow out CI around early-type stars, can be suppressed by self-shielding.

This paper considers whether the population of known transiting exoplanets provides evidence for additional outer planets on inclined orbits, due to the perturbing effect of such planets on the orbits of inner planets. As such, we develop a semi-analytical method for calculating the probability that two mutually inclined planets are observed to transit. We subsequently derive a simplified analytical form to describe how the mutual inclination between two planets evolves due to secular interactions with a wide orbit inclined planet and use this to determine the mean probability that the two inner planets are observed to transit. From application to Kepler-48 and HD-106315 we constrain the inclinations of the outer planets in these systems (known from RV). We also apply this work to the so called Kepler Dichotomy, which describes the excess of single transiting systems observed by Kepler. We find 3 different ways of explaining this dichotomy: some systems could be inherently single, some multi-planet systems could have inherently large mutual inclinations, while some multi-planet systems could cyclically attain large mutual inclinations through interaction with an inclined outer planet. We show how the different mechanisms can be combined to fit the observed populations of Kepler systems with one and two transiting planets. We also show how the distribution of mutual inclinations of transiting two planet systems constrains the fraction of two planet systems that have perturbing outer planets, since such systems should be preferentially discovered by Kepler when the inner planets are coplanar due to an increased transit probability.

This prospective chapter gives our view on the evolution of the study of circumstellar discs within the next 20 years from both observational and theoretical sides. We first present the expected improvements in our knowledge of protoplanetary discs as for their masses, sizes, chemistry, the presence of planets as well as the evolutionary processes shaping these discs. We then explore the older debris disc stage and explain what will be learnt concerning their birth, the intrinsic links between these discs and planets, the hot dust and the gas detected around main sequence stars as well as discs around white dwarfs.

Infrared excesses around metal polluted white dwarfs have been associated with the accretion of dusty, planetary material. This work analyses the available infrared data for an unbiased sample of white dwarfs and demonstrates that no more than 3.3% can have a wide, flat, opaque dust disc, extending to the Roche radius, with a temperature at the disc inner edge of $T_{in} = 1,400$K, the standard model for the observed excesses. This is in stark contrast to the incidence of pollution of about 30%. We present four potential reasons for the absence of an infrared excess in polluted white dwarfs, depending on their stellar properties and inferred accretion rates: i) their dust discs are opaque, but narrow, thus evading detection if more than 85% of polluted white dwarfs have dust discs narrower than $\delta r< 0.04r$, ii) their dust discs have been fully consumed, which only works for the oldest white dwarfs with sinking timescales longer than hundreds of years, iii) their dust is optically thin, which can supply low accretion rates of $<10^7 $gs$^{-1}$ if dominated by PR-drag, and higher accretion rates, if inwards transport of material is enhanced, for example due to the presence of gas, iv) their accretion is supplied by a pure gas disc, which could result from the sublimation of optically thin dust for T* > 20, 000K. Future observations sensitive to faint infrared excesses or the presence of gas, can test the scenarios presented here, thereby better constraining the nature of the material fuelling accretion in polluted white dwarfs.

RZ Psc is a young Sun-like star, long associated with the UXor class of variable stars, which is partially or wholly dimmed by dust clumps several times each year. The system has a bright and variable infrared excess, which has been interpreted as evidence that the dimming events are the passage of asteroidal fragments in front of the host star. Here, we present a decade of optical photometry of RZ Psc and take a critical look at the asteroid belt interpretation. We show that the distribution of light curve gradients is non-uniform for deep events, which we interpret as possible evidence for an asteroidal fragment-like clump structure. However, the clumps are very likely seen above a high optical depth mid-plane, so the disk's bulk clumpiness is not revealed. While circumstantial evidence suggests an asteroid belt is more plausible than a gas-rich transition disk, the evolutionary status remains uncertain. We suggest that the rarity of Sun-like stars showing disk-related variability may arise because i) any accretion streams are transparent, and/or ii) turbulence above the inner rim is normally shadowed by a flared outer disk.

We present the first scattered-light images of the debris disk around 49 ceti, a ~40 Myr A1 main sequence star at 59 pc, famous for hosting two massive dust belts as well as large quantities of atomic and molecular gas. The outer disk is revealed in reprocessed archival Hubble Space Telescope NICMOS F110W images, as well as new coronagraphic H band images from the Very Large Telescope SPHERE instrument. The disk extends from 1.1" (65 AU) to 4.6" (250 AU), and is seen at an inclination of 73degr, which refines previous measurements at lower angular resolution. We also report no companion detection larger than 3 M_Jup at projected separations beyond 20 AU from the star (0.34"). Comparison between the F110W and H-band images is consistent with a grey color of 49 ceti's dust, indicating grains larger than >2microns. Our photometric measurements indicate a scattering efficiency / infrared excess ratio of 0.2-0.4, relatively low compared to other characterized debris disks. We find that 49 ceti presents morphological and scattering properties very similar to the gas-rich HD 131835 system. From our constraint on the disk inclination we find that the atomic gas previously detected in absorption must extend to the inner disk, and that the latter must be depleted of CO gas. Building on previous studies, we propose a schematic view of the system describing the dust and gas structure around 49 ceti and hypothetic scenarios for the gas nature and origin.

Uncovering the ingredients and the architecture of planetary systems is a very active field of research that has fuelled many new theories on giant planet formation, migration, composition, and interaction with the circumstellar environment. We aim at discovering and studying new such systems, to further expand our knowledge of how low-mass companions form and evolve. We obtained high-contrast H-band images of the circumstellar environment of the F5V star HD206893, known to host a debris disc never detected in scattered light. These observations are part of the SPHERE High Angular Resolution Debris Disc Survey (SHARDDS) using the InfraRed Dual-band Imager and Spectrograph (IRDIS) installed on VLT/SPHERE. We report the detection of a source with a contrast of 3.6x10^{-5} in the H-band, orbiting at a projected separation of 270 milliarcsecond or 10 au, corresponding to a mass in the range 24 to 73 Mjup for an age of the system in the range 0.2 to 2 Gyr. The detection was confirmed ten months later with VLT/NaCo, ruling out a background object with no proper motion. A faint extended emission compatible with the disc scattered light signal is also observed. The detection of a low-mass companion inside a massive debris disc makes this system an analog of other young planetary systems such as beta Pictoris, HR8799 or HD95086 and requires now further characterisation of both components to understand their interactions.

We present the first resolved image of the debris disk around the 16+/-8 Myr old star, HD 114082. The observation was made in the H-band using the SPHERE instrument. The star is at a distance of 92+/-6 pc in the Lower Centaurus Crux association. Using a Markov Chain Monte Carlo analysis, we determined that the debris is likely in the form of a dust ring with an inner edge of 27.7+2.8/-3.5 au, position angle -74+0.5/-1.5 deg, and an inclination with respect to the line of sight of 6.7+3.8/-0.4 deg. The disk imaged in scattered light has a surface density declining with radius like ~r^(-4), steeper than expected for grain blowout by radiation pressure. We find only marginal evidence (2 sigma) of eccentricity, and rule out planets more massive than 1.0 Mjup orbiting within 1 au of the ring's inner edge, since such a planet would have disrupted the disk. The disk has roughly the same fractional disk luminosity (Ldisk/L*=3.3x10^(-3)) as HR4796A and Beta Pictoris, however it was not detected by previous instrument facilities most likely because of its small angular size (radius~0.4"), low albedo (~0.2) and low scattering efficiency far from the star due to high scattering anisotropy. With the arrival of extreme adaptive optics systems like SPHERE and GPI, the morphology of smaller, fainter and more distant debris disks are being revealed, providing clues to planet-disk interactions in young protoplanetary systems.

While most of the known debris discs present cold dust at tens of AU, a few young systems exhibit hot dust analogous to the Zodiacal dust. $\eta$ Corvi is particularly interesting as it is old and it has both, with its hot dust significantly exceeding the maximum luminosity of an in-situ collisional cascade. Previous work suggested that this system could be undergoing an event similar to the Late Heavy Bombardment (LHB) soon after or during a dynamical instability. Here we present ALMA observations of $\eta$ Corvi with a resolution of 1."2 (~22au) to study its outer belt. The continuum emission is consistent with an axisymmetric belt, with a mean radius of 152au and radial FWHM of 46au, which is too narrow compared to models of inward scattering of an LHB-like scenario. Instead, the hot dust could be explained as material passed inwards in a rather stable planetary configuration. We also report a 4sigma detection of CO at ~ 20au. CO could be released in situ from icy planetesimals being passed in when crossing the H$_2$O or CO$_2$ ice lines. Finally, we place constraints on hidden planets in the disc. If a planet is sculpting the disc's inner edge, this should be orbiting at 75-100au, with a mass of 3-30 M$_\oplus$ and an eccentricity < 0.08. Such a planet would be able to clear its chaotic zone on a timescale shorter than the age of the system and scatter material inwards from the outer belt to the inner regions, thus feeding the hot dust.

This paper considers the dynamics of scattering of planetesimals or planetary embryos by a planet on a circumstellar orbit. We classify six regions in the planet's mass versus semimajor axis parameter space according to the dominant outcome for scattered objects: ejected, accreted, remaining, escaping, Oort Cloud, depleted Oort Cloud. We use these outcomes to consider which planetary system architectures maximise the observability of specific signatures, given that signatures should be detected first around systems with optimal architectures (if such systems exist in nature). Giant impact debris is most readily detectable for 0.1-10Mearth planets at 1-5au, depending on detection method and spectral type. While A stars have putative giant impact debris at 4-6au consistent with this sweet spot, that of FGK stars is typically <<1au contrary to expectations; an absence of 1-3au giant impact debris could indicate a low frequency of terrestrial planets there. Three principles maximise cometary influx from exo-Kuiper belts: a chain of closely separated planets interior to the belt, none of which is a Jupiter-like ejector; planet masses not increasing strongly with distance (for a net inward torque on comets); ongoing replenishment of comets, possibly by embedded low-mass planets. A high Oort Cloud comet influx requires no ejector and architectures that maximise the Oort Cloud population. Cold debris disks are usually considered classical Kuiper belt analogues. Here we consider the possibility of detecting scattered disk analogues, which could be betrayed by a broad radial profile and lack of small grains, as well as spherical 100-1000au mini-Oort Clouds. Some implications for escaping planets around young stars, detached planets akin to Sedna, and the formation of super-Earths, are also discussed.

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.

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.

This paper presents a self-consistent model for the evolution of gas produced in the debris disc of $\beta$ Pictoris. Our model proposes that atomic carbon and oxygen are created from the photodissociation of CO, which is itself released from volatile-rich bodies in the debris disc due to grain-grain collisions or photodesorption. While the CO lasts less than one orbit, the atomic gas evolves by viscous spreading resulting in an accretion disc inside the parent belt and a decretion disc outside. The temperature, ionisation fraction and population levels of carbon and oxygen are followed with the photodissociation region model Cloudy, which is coupled to a dynamical viscous $\alpha$ model. We present new gas observations of $\beta$ Pic, of C I observed with APEX and O I observed with Herschel, and show that these along with published C II and CO observations can all be explained with this new model. Our model requires a viscosity $\alpha$ > 0.1, similar to that found in sufficiently ionised discs of other astronomical objects; we propose that the magnetorotational instability is at play in this highly ionised and dilute medium. This new model can be tested from its predictions for high resolution ALMA observations of C I. We also constrain the water content of the planetesimals in $\beta$ Pic. The scenario proposed here might be at play in all debris discs and this model could be used more generally on all discs with C, O or CO detections.

An increasing number of observations have shown that gaseous debris discs are not an exception. However, until now we only knew of cases around A stars. Here we present the first detection of 12CO (2-1) disc emission around an F star, HD 181327, obtained with ALMA observations at 1.3 mm. The continuum and CO emission are resolved into an axisymmetric disc with ring-like morphology. Using a Markov chain Monte Carlo method coupled with radiative transfer calculations we study the dust and CO mass distribution. We find the dust is distributed in a ring with a radius of 86.0 +- 0.4 AU and a radial width of 23.2 +- 1.0 AU. At this frequency the ring radius is smaller than in the optical, revealing grain size segregation expected due to radiation pressure. We also report on the detection of low level continuum emission beyond the main ring out to ~200 AU. We model the CO emission in the non-LTE regime and we find that the CO is co-located with the dust, with a total CO gas mass ranging between 1.2x10^-6 Mearth and 2.9x10^-6 Mearth, depending on the gas kinetic temperature and collisional partners densities. The CO densities and location suggest a secondary origin, i.e. released from icy planetesimals in the ring. We derive a CO+CO2 cometary composition that is consistent with Solar system comets. Due to the low gas densities it is unlikely that the gas is shaping the dust distribution.

Debris discs are typically revealed through excess emission at infrared wavelengths. Most discs exhibit excess at mid- and far-infrared wavelengths, analogous to the solar system's Asteroid and Edgeworth-Kuiper belts. Recently, stars with strong (1 per cent) excess at near-infrared wavelengths were identified through interferometric measurements. Using the HIgh Precision Polarimetric Instrument (HIPPI), we examined a sub-sample of these hot dust stars (and appropriate controls) at parts-per-million sensitivity in SDSS g' (green) and r' (red) filters for evidence of scattered light. No detection of strongly polarized emission from the hot dust stars is seen. We therefore rule out scattered light from a normal debris disk as the origin of this emission. A wavelength dependent contribution from multiple dust components for hot dust stars is inferred from the dispersion (difference in polarization angle in red and green) of southern stars. Contributions of 17 ppm (green) and 30 ppm (red) are calculated, with strict 3 sigma upper limits of 76 and 68 ppm, respectively. This suggests weak hot dust excesses consistent with thermal emission, although we cannot rule out contrived scenarios, e.g. dust in a spherical shell or face on discs. We also report on the nature of the local interstellar medium, obtained as a byproduct of the control measurements. Highlights include the first measurements of the polarimetric colour of the local interstellar medium and discovery of a southern sky region with a polarization per distance thrice the previous maximum. The data suggest the wavelength of maximum polarization is bluer than typical.

Giant impacts refer to collisions between two objects each of which is massive enough to be considered at least a planetary embryo. The putative collision suffered by the proto-Earth that created the Moon is a prime example, though most Solar System bodies bear signatures of such collisions. Current planet formation models predict that an epoch of giant impacts may be inevitable, and observations of debris around other stars are providing mounting evidence that giant impacts feature in the evolution of many planetary systems. This chapter reviews giant impacts, focussing on what we can learn about planet formation by studying debris around other stars. Giant impact debris evolves through mutual collisions and dynamical interactions with planets. General aspects of this evolution are outlined, noting the importance of the collision-point geometry. The detectability of the debris is discussed using the example of the Moon-forming impact. Such debris could be detectable around another star up to 10Myr post-impact, but model uncertainties could reduce detectability to a few 100yr window. Nevertheless the 3% of young stars with debris at levels expected during terrestrial planet formation provide valuable constraints on formation models; implications for super-Earth formation are also discussed. Variability recently observed in some bright disks promises to illuminate the evolution during the earliest phases when vapour condensates may be optically thick and acutely affected by the collision-point geometry. The outer reaches of planetary systems may also exhibit signatures of giant impacts, such as the clumpy debris structures seen around some stars.

The Large Binocular Telescope Interferometer (LBTI) is a versatile instrument designed for high-angular resolution and high-contrast infrared imaging (1.5-13 microns). In this paper, we focus on the mid-infrared (8-13 microns) nulling mode and present its theory of operation, data reduction, and on-sky performance as of the end of the commissioning phase in March 2015. With an interferometric baseline of 14.4 meters, the LBTI nuller is specifically tuned to resolve the habitable zone of nearby main-sequence stars, where warm exozodiacal dust emission peaks. Measuring the exozodi luminosity function of nearby main-sequence stars is a key milestone to prepare for future exoEarth direct imaging instruments. Thanks to recent progress in wavefront control and phase stabilization, as well as in data reduction techniques, the LBTI demonstrated in February 2015 a calibrated null accuracy of 0.05% over a three-hour long observing sequence on the bright nearby A3V star beta Leo. This is equivalent to an exozodiacal disk density of 15 to 30 zodi for a Sun-like star located at 10pc, depending on the adopted disk model. This result sets a new record for high-contrast mid-infrared interferometric imaging and opens a new window on the study of planetary systems.

Over the duration of the Kepler mission, KIC8462852 was observed to undergo irregularly shaped, aperiodic dips in flux of up to $\sim 20$\%. The dipping activity can last for between 5 and 80 days. We characterize the object with high-resolution spectroscopy, spectral energy distribution fitting, radial velocity measurements, high-resolution imaging, and Fourier analyses of the Kepler light curve. We determine that KIC8462852 is a typical main-sequence F3 V star that exhibits no significant IR excess, and has no very close interacting companions. In this paper, we describe various scenarios to explain the dipping events observed in the Kepler light curve. We confirm that the dipping signals in the data are not caused by any instrumental or data processing artifact, and thus are astrophysical in origin. We construct scenario-independent constraints on the size and location of a body in the system that is needed to reproduce the observations. We deliberate over several assorted stellar and circumstellar astrophysical scenarios, most of which have problems explaining the data in hand. By considering the observational constraints on dust clumps in orbit around a normal main-sequence star, we conclude that the scenario most consistent with the data in hand is the passage of a family of exocomet or planetesimal fragments, all of which are associated with a single previous break-up event, possibly caused by tidal disruption or thermal processing. The minimum total mass associated with these fragments likely exceeds $10^{-6}$~\mearth, corresponding to an original rocky body of $>100$~km in diameter. We discuss the necessity of future observations to help interpret the system.

This paper considers secular interactions within multi-planet systems. In particular we consider dynamical evolution of known planetary systems resulting from an additional hypothetical planet on an eccentric orbit. We start with an analytical study of a general two-planet system, showing that a planet on an elliptical orbit transfers all of its eccentricity to an initially circular planet if the two planets have comparable orbital angular momenta. Application to the single Super-Earth system HD38858 shows that an additional hypothetical planet below current radial velocity (RV) constraints with {\textit{Msini}}=3-10M$_\oplus$, semi-major axis 1-10au and eccentricity 0.2-0.8 is unlikely to be present from the eccentricity that would be excited in the known planet (albeit cyclically). However, additional planets in proximity to the known planet could stabilise the system against secular perturbations from outer planets. Moreover these additional planets can have an {\textit{Msini}} below RV sensitivity and still affect their neighbours. For example, application to the two Super-Earth system 61Vir shows that an additional hypothetical planet cannot excite high eccentricities in the known planets, unless its mass and orbit lie in a restricted area of parameter space. Inner planets in HD38858 below RV sensitivity would also modify conclusions above about excluded parameter space. This suggests that it may be possible to infer the presence of additional stabilising planets in systems with an eccentric outer planet and an inner planet on an otherwise suspiciously circular orbit. This reinforces the point that the full complement of planets in a system is needed to assess its dynamical state.

We present ten young ($\lesssim$10 Myr) late-K and M dwarf stars observed in K2 Campaign 2 that host protoplanetary disks and exhibit quasi-periodic or aperiodic dimming events. Their optical light curves show $\sim$10-20 dips in flux over the 80-day observing campaign with durations of $\sim$0.5-2 days and depths of up to $\sim$40%. These stars are all members of the $\rho$ Ophiuchus ($\sim$1 Myr) or Upper Scorpius ($\sim$10 Myr) star-forming regions. To investigate the nature of these "dippers" we obtained: optical and near-infrared spectra to determine stellar properties and identify accretion signatures; adaptive optics imaging to search for close companions that could cause optical variations and/or influence disk evolution; and millimeter-wavelength observations to constrain disk dust and gas masses. The spectra reveal Li I absorption and H$\alpha$ emission consistent with stellar youth (<50 Myr), but also accretion rates spanning those of classical and weak-line T Tauri stars. Infrared excesses are consistent with protoplanetary disks extending to within $\sim$10 stellar radii in most cases; however, the sub-mm observations imply disk masses that are an order of magnitude below those of typical protoplanetary disks. We find a positive correlation between dip depth and WISE-2 excess, which we interpret as evidence that the dipper phenomenon is related to occulting structures in the inner disk, although this is difficult to reconcile with the weakly accreting aperiodic dippers. We consider three mechanisms to explain the dipper phenomenon: inner disk warps near the co-rotation radius related to accretion; vortices at the inner disk edge produced by the Rossby Wave Instability; and clumps of circumstellar material related to planetesimal formation.

Lambda Boo stars are predominately A-type stars with solar abundant C, N, O, and S, but up to 2 dex underabundances of refractory elements. The stars' unusual surface abundances could be due to a selective accretion of volatile gas over dust. It has been proposed that there is a correlation between the Lambda Boo phenomenon and IR-excesses which are the result of a debris disk or interstellar medium (ISM) interaction providing the accreting material. We observe 70 or 100 and 160 $\mu$m excess emission around 9 confirmed Lambda Boo stars with the Herschel Space Observatory, to differentiate whether the dust emission is from a debris disk or an ISM bow wave. We find that 3/9 stars observed host well resolved debris disks. While the remaining 6/9 are not resolved, they are inconsistent with an ISM bow wave based on the dust emission being more compact for its temperature and predicted bow wave models produce hotter emission than what is observed. We find the incidence of bright IR-excesses around Lambda Boo stars is higher than normal A-stars. To explain this given our observations, we explore Poynting-Robertson (PR) drag as a mechanism of accretion from a debris disk but find it insufficient. As an alternative, we propose the correlation is due to higher dynamical activity in the disks currently underway. Large impacts of planetesimals or a higher influx of comets could provide enough volatile gas for accretion. Further study on the transport of circumstellar material in relation to the abundance anomalies are required to explain the phenomenon through external accretion.

This paper discusses compelling science cases for a future long-baseline interferometer operating at millimeter and centimeter wavelengths, like the proposed Next Generation Vary Large Array (ngVLA). We report on the activities of the Cradle of Life science working group, which focused on the formation of low- and high-mass stars, the formation of planets and evolution of protoplanetary disks, the physical and compositional study of Solar System bodies, and the possible detection of radio signals from extraterrestrial civilizations. We propose 19 scientific projects based on the current specification of the ngVLA. Five of them are highlighted as possible Key Science Projects: (1) Resolving the density structure and dynamics of the youngest HII regions and high-mass protostellar jets, (2) Unveiling binary/multiple protostars at higher resolution, (3) Mapping planet formation regions in nearby disks on scales down to 1 AU, (4) Studying the formation of complex molecules, and (5) Deep atmospheric mapping of giant planets in the Solar System. For each of these projects, we discuss the scientific importance and feasibility. The results presented here should be considered as the beginning of a more in-depth analysis of the science enabled by such a facility, and are by no means complete or exhaustive.

We investigate the general interaction between an eccentric planet and a coplanar debris disc of the same mass, using analytical theory and n-body simulations. Such an interaction could result from a planet-planet scattering or merging event. We show that when the planet mass is comparable to that of the disc, the former is often circularised with little change to its semimajor axis. The secular effect of such a planet can cause debris to apsidally anti-align with the planet's orbit (the opposite of what may be naively expected), leading to the counter-intuitive result that a low-mass planet may clear a larger region of debris than a higher-mass body would. The interaction generally results in a double-ringed debris disc, which is comparable to those observed in HD 107146 and HD 92945. As an example we apply our results to HD 107146, and show that the disc's morphology and surface brightness profile can be well-reproduced if the disc is interacting with an eccentric planet of comparable mass (~10-100 Earth masses). This hypothetical planet had a pre-interaction semimajor axis of 30 or 40 au (similar to its present-day value) and an eccentricity of 0.4 or 0.5 (which would since have reduced to ~0.1). Thus the planet (if it exists) presently resides near the inner edge of the disc, rather than between the two debris peaks as may otherwise be expected. Finally we show that disc self-gravity can be important in this mass regime and, whilst it would not affect these results significantly, it should be considered when probing the interaction between a debris disc and a planet.

The spatially resolved AU Mic debris disc is among the most famous and best-studied debris discs. We aim at a comprehensive understanding of the dust production and the dynamics of the disc objects with in depth collisional modelling including stellar radiative and corpuscular forces. Our models are compared to a suite of observational data for thermal and scattered light emission, ranging from the ALMA radial surface brightness profile at 1.3mm to polarisation measurements in the visible. Most of the data can be reproduced with a planetesimal belt having an outer edge at around 40au and subsequent inward transport of dust by stellar winds. A low dynamical excitation of the planetesimals with eccentricities up to 0.03 is preferred. The radial width of the planetesimal belt cannot be constrained tightly. Belts that are 5au and 17au wide, as well as a broad 44au-wide belt are consistent with observations. All models show surface density profiles increasing with distance from the star as inferred from observations. The best model is achieved by assuming a stellar mass loss rate that exceeds the solar one by a factor of 50. While the SED and the shape of the ALMA profile are well reproduced, the models deviate from the scattered light data more strongly. The observations show a bluer disc colour and a lower degree of polarisation for projected distances <40au than predicted by the models. The problem may be mitigated by irregularly-shaped dust grains which have scattering properties different from the Mie spheres used. From tests with a handful of selected dust materials, we derive a preference for mixtures of silicate, carbon, and ice of moderate porosity. We address the origin of the unresolved central excess emission detected by ALMA and show that it cannot stem from an additional inner belt alone. Instead, it should derive, at least partly, from the chromosphere of the central star.

A pathway to the formation of planetesimals, and eventually giant planets, may occur in concentrations of dust grains trapped in pressure maxima. Dramatic crescent-shaped dust concentrations have been seen in recent radio images at sub-mm wavelengths. These disk asymmetries could represent the initial phases of planet formation in the dust trap scenario, provided that grain sizes are spatially segregated. A testable prediction of azimuthal dust trapping is that progressively larger grains should be more sharply confined and furthermore the trapped grains should follow a distribution that is markedly different from the gas. However, gas tracers such as CO and the infrared emission from small grains are both very optically thick where the submm continuum originates, so observations have been unable to test the trapping predictions or to identify compact concentrations of larger grains required for planet formation by core-accretion. Here we report multifrequency observations of HD142527, from 34GHz to 700GHz, that reveal a compact concentration of ~cm-sized grains, with a few Earth masses, embedded in a large-scale crescent of ~mm-sized particles. The emission peaks at wavelengths shorter than ~1mm are optically thick and trace the temperature structure resulting from shadows cast by the inner regions. Given this temperature structure, we infer that the largest dust grains are concentrated in the 34 GHz clump. We conclude that dust trapping is efficient for approximately cm-sized grains and leads to enhanced concentrations, while the smaller grains largely reflect the gas distribution.

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.

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.

Imaging a star's companion at multiple epochs over a short orbital arc provides only four of the six coordinates required for a unique orbital solution. Probability distributions of possible solutions are commonly generated by Monte Carlo (MCMC) analysis, but these are biased by priors and may not probe the full parameter space. We suggest alternative methods to characterise possible orbits, which compliment the MCMC technique. Firstly the allowed ranges of orbital elements are prior-independent, and we provide means to calculate these ranges without numerical analyses. Hence several interesting constraints (including whether a companion even can be bound, its minimum possible semi-major axis and its minimum eccentricity) may be quickly computed using our relations as soon as orbital motion is detected. We also suggest an alternative to posterior probability distributions as a means to present possible orbital elements, namely contour plots of elements as functions of line of sight coordinates. These plots are prior-independent, readily show degeneracies between elements and allow readers to extract orbital solutions themselves. This approach is particularly useful when there are other constraints on the geometry, for example if a companion's orbit is assumed to be aligned with a disc. As examples we apply our methods to several imaged sub-stellar companions including Fomalhaut b, and for the latter object we show how different origin hypotheses affect its possible orbital solutions. We also examine visual companions of A- and G-type main sequence stars in the Washington Double Star Catalogue, and show that $\gtrsim50$ per cent must be unbound.

We report on the first nulling interferometric observations with the Large Binocular Telescope Interferometer (LBTI), resolving the N' band (9.81 - 12.41 um) emission around the nearby main-sequence star eta Crv (F2V, 1-2 Gyr). The measured source null depth amounts to 4.40% +/- 0.35% over a field-of-view of 140 mas in radius (~2.6\,AU at the distance of eta Corvi) and shows no significant variation over 35{\deg} of sky rotation. This relatively low null is unexpected given the total disk to star flux ratio measured by Spitzer/IRS (~23% across the N' band), suggesting that a significant fraction of the dust lies within the central nulled response of the LBTI (79 mas or 1.4 AU). Modeling of the warm disk shows that it cannot resemble a scaled version of the Solar zodiacal cloud, unless it is almost perpendicular to the outer disk imaged by Herschel. It is more likely that the inner and outer disks are coplanar and the warm dust is located at a distance of 0.5-1.0 AU, significantly closer than previously predicted by models of the IRS spectrum (~3 AU). The predicted disk sizes can be reconciled if the warm disk is not centrosymmetric, or if the dust particles are dominated by very small grains. Both possibilities hint that a recent collision has produced much of the dust. Finally, we discuss the implications for the presence of dust at the distance where the insolation is the same as Earth's (2.3 AU).

We present deep HST/STIS coronagraphic images of the Beta Pic debris disk obtained at two epochs separated by 15 years. The new images and the re-reduction of the 1997 data provide the most sensitive and detailed views of the disk at optical wavelengths as well as the yet smallest inner working angle optical coronagraphic image of the disk. Our observations characterize the large-scale and inner-disk asymmetries and we identify multiple breaks in the disk radial surface brightness profile. We study in detail the radial and vertical disk structure and show that the disk is warped. We explore the disk at the location of the Beta Pic b super-jupiter and find that the disk surface brightness slope is continuous between 0.5 and 2.0 arcsec, arguing for no change at the separations where Beta Pic b orbits. The two epoch images constrain the disk surface brightness evolution on orbital and radiation pressure blow-out timescales. We place an upper limit of 3% on the disk surface brightness change between 3-5 arcsec, including the locations of the disk warp, and the CO and dust clumps. We discuss the new observations in the context of high-resolution multi-wavelength images and divide the disk asymmetries in two groups: axisymmetric and non-axisymmetric. The axisymmetric structures (warp, large-scale butterfly, etc.) are consistent with disk structure models that include interactions of a planetesimal belt and a non-coplanar giant planet. The non-axisymmetric features, however, require a different explanation.

Dust grains migrating under Poynting-Robertson drag may be trapped in mean-motion resonances with planets. Such resonantly trapped grains are observed in the solar system. In extrasolar systems, the exozodiacal light produced by dust grains is expected to be a major obstacle to future missions attempting to directly image terrestrial planets. The patterns made by resonantly trapped dust, however, can be used to infer the presence of planets, and the properties of those planets, if the capture and evolution of the grains can be modelled. This has been done with N-body methods, but such methods are computationally expensive, limiting their usefulness when considering large, slowly evolving grains, and for extrasolar systems with unknown planets and parent bodies, where the possible parameter space for investigation is large. In this work, we present a semi-analytic method for calculating the capture and evolution of dust grains in resonance, which can be orders of magnitude faster than N-body methods. We calibrate the model against N-body simulations, finding excellent agreement for Earth to Neptune mass planets, for a variety of grain sizes, initial eccentricities, and initial semimajor axes. We then apply the model to observations of dust resonantly trapped by the Earth. We find that resonantly trapped, asteroidally produced grains naturally produce the `trailing blob' structure in the zodiacal cloud, while to match the intensity of the blob, most of the cloud must be composed of cometary grains, which owing to their high eccentricity are not captured, but produce a smooth disk.

The Hunt for Observable Signatures of Terrestrial planetary Systems (HOSTS) on the Large Binocular Telescope Interferometer will survey nearby stars for faint emission arising from ~300 K dust (exozodiacal dust), and aims to determine the exozodiacal dust luminosity function. HOSTS results will enable planning for future space telescopes aimed at direct spectroscopy of habitable zone terrestrial planets, as well as greater understanding of the evolution of exozodiacal disks and planetary systems. We lay out here the considerations that lead to the final HOSTS target list. Our target selection strategy maximizes the ability of the survey to constrain the exozodi luminosity function by selecting a combination of stars selected for suitability as targets of future missions and as sensitive exozodi probes. With a survey of approximately 50 stars, we show that HOSTS can enable an understanding of the statistical distribution of warm dust around various types of stars and is robust to the effects of varying levels of survey sensitivity induced by weather conditions.

The protoplanetary disks of Herbig Ae stars eventually dissipate leaving a tenuous debris disk comprised of planetesimals and dust, as well as possibly gas and planets. This paper uses the properties of 10-20Myr A star debris disks to consider the protoplanetary to debris disk transition. The physical distinction between these two classes is argued to rest on the presence of primordial gas in sufficient quantities to dominate the motion of small dust grains (not the secondary nature of the dust or its level of stirring). This motivates an observational classification based on the dust spectrum, empirically defined so that A star debris disks require fractional excesses <3 at 12um and <2000 at 70um. We also propose a hypothesis to test, that the main sequence planet/planetesimal structures are already in place (but obscured) during the protoplanetary disk phase. This may be only weakly true if planetary architectures change until frozen during disk dispersal, or completely false if planets and planetesimals form during disk dispersal. Five steps in the transition are discussed: (i) carving an inner hole to form a transition disk; (ii) depletion of mm-sized dust in outer disk, noting the importance of determining whether this mass ends up in planetesimals or is collisionally depleted; (iii) final clearing of inner regions, noting that many mechanisms replenish moderate hot dust levels at later phases, and likely also operate in protoplanetary disks; (iv) disappearence of gas, noting recent discoveries of primordial and secondary gas in debris disks that highlight our ignorance and its impending enlightenment by ALMA; (v) formation of ring-like planetesimal structures, noting these are shaped by interactions with planets, and that the location of planetesimals in protoplanetary disks may be unrelated to the dust concentrations therein that are set by gas interactions.

In recent years, gas has been observed in an increasing number of debris discs, though its nature remains to be determined. Here, we analyse CO molecular excitation in optically thin debris discs, and search ALMA Cycle-0 data for CO J=3-2 emission in the Fomalhaut ring. No significant line emission is observed; we set a 3-$\sigma$ upper limit on the integrated line flux of 0.16 Jy km s$^{-1}$. We show a significant dependency of the CO excitation on the density of collisional partners $n$, on the gas kinetic temperature $T_k$ and on the ambient radiation field $J$, suggesting that assumptions widely used for protoplanetary discs (e.g. LTE) do not necessarily apply to their low density debris counterparts. When applied to the Fomalhaut ring, we consider a primordial origin scenario where H$_2$ dominates collisional excitation of CO, and a secondary origin scenario dominated by e$^-$ and H$_2$O. In either scenario, we obtain a strict upper limit on the CO mass of 4.9 $\times$ 10$^{-4}$ M$_{\oplus}$. This arises in the non-LTE regime, where the excitation of the molecule is determined solely by the well-known radiation field. In the secondary scenario, assuming any CO present to be in steady state allows us to set an upper limit of $\sim$55% on the CO/H$_2$O ice ratio in the parent planetesimals. This could drop to $\sim$3% if LTE applies, covering the range observed in Solar System comets (0.4-30%). Finally, in light of our analysis, we present prospects for CO detection and characterisation in debris discs with ALMA.

Habitable zone dust levels are a key unknown that must be understood to ensure the success of future space missions to image Earth analogues around nearby stars. Current detection limits are several orders of magnitude above the level of the Solar System's Zodiacal cloud, so characterisation of the brightness distribution of exo-zodi down to much fainter levels is needed. To this end, the large Binocular Telescope Interferometer (LBTI) will detect thermal emission from habitable zone exo-zodi a few times brighter than Solar System levels. Here we present a modelling framework for interpreting LBTI observations, which yields dust levels from detections and upper limits that are then converted into predictions and upper limits for the scattered light surface brightness. We apply this model to the HOSTS survey sample of nearby stars; assuming a null depth uncertainty of 10$^{-4}$ the LBTI will be sensitive to dust a few times above the Solar System level around Sun-like stars, and to even lower dust levels for more massive stars.

The Oort cloud is usually thought of as a collection of icy comets inhabiting the outer reaches of the Solar system, but this picture is incomplete. We use simulations of the formation of the Oort cloud to show that ~4% of the small bodies in the Oort cloud should have formed within 2.5 au of the Sun, and hence be ice-free rock-iron bodies. If we assume these Oort cloud asteroids have the same size distribution as their cometary counterparts, the Large Synoptic Survey Telescope should find roughly a dozen Oort cloud asteroids during ten years of operations. Measurement of the asteroid fraction within the Oort cloud can serve as an excellent test of the Solar system's formation and dynamical history. Oort cloud asteroids could be of particular concern as impact hazards as their high mass density, high impact velocity, and low visibility make them both hard to detect and hard to divert or destroy. However, they should be a rare class of object, and we estimate globally catastrophic collisions should only occur about once per billion years.

We present a study of debris disks whose spectra are well modelled by dust emission at two different temperatures. These disks are typically assumed to be a sign of multiple belts, which in only a few cases have been confirmed via high resolution observations. We first compile a sample of two-temperature disks to derive their properties, summarised by the ratios of the warm and cool component temperatures and fractional luminosities. The ratio of warm to cool temperatures is constant in the range 2-4, and the temperatures of both warm and cool components increases with stellar mass. We then explore whether this emission can arise from dust in a single narrow belt, with the range of temperatures arising from the size variation of grain temperatures. This model can produce two-temperature spectra for Sun-like stars, but is not supported where it can be tested by observed disk sizes and far-IR/mm spectral slopes. Therefore, while some two-temperature disks arise from single belts, it is probable that most have multiple spatial components. These disks are plausibly similar to the outer Solar System's configuration of Asteroid and Edgeworth-Kuiper belts separated by giant planets. Alternatively, the inner component could arise from inward scattering of material from the outer belt, again due to intervening planets. In either case, we suggest that the ratio of warm/cool component temperatures is indicative of the scale of outer planetary systems, which typically span a factor of about ten in radius.

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.

We investigate the interaction between an eccentric planet and a less massive external debris disc. This scenario could occur after planet-planet scattering or merging events. We characterise the evolution over a wide range of initial conditions, using a suite of n-body integrations combined with theory. Planets near the disc mid-plane remove the inner debris region, and surviving particles form an eccentric disc apsidally aligned with the planet. The inner disc edge is elliptical and lies just beyond the planet's orbit. Moderately inclined planets ($i_{\rm plt} \gtrsim 20^\circ$ for $e_{\rm plt} = 0.8$) may instead sculpt debris into a bell-shaped structure enveloping the planet's orbit. Finally some highly inclined planets ($i_{\rm plt} \sim 90^\circ$) can maintain a disc orthogonal to the planet's plane. In all cases disc particles undergo rapid evolution, whilst the overall structures evolve more slowly. The shapes of these structures and their density profiles are characterised. The width of the chaotic zone around the planet's orbit is derived in the coplanar case using eccentric Hill radius arguments. This zone is cleared within approximately ten secular or diffusion times (whichever is longer), and debris assumes its final shape within a few secular times. We quantify the planet's migration and show it will almost always be small in this mass regime. Our results may be used to characterise unseen eccentric planets using observed debris features.

Fomalhaut is a triple system, with all components widely separated (~1E5 au). Such widely separated binaries are thought to form during cluster dissolution, but that process is unlikely to form such a triple system. We explore an alternative scenario, where A and C form as a tighter binary from a single molecular cloud core (with semimajor axis ~1E4 au), and B is captured during cluster dispersal. We use N-body simulations augmented with the Galactic tidal forces to show that such a system naturally evolves into a Fomalhaut-like system in about half of cases, on a timescale compatible with the age of Fomalhaut. From initial non-interacting orbits, Galactic tides drive cycles in B's eccentricity that lead to a close encounter with C. After several close encounters, typically lasting tens of millions of years, one of the stars is ejected. The Fomalhaut-like case with both components at large separations is almost invariably a precursor to the ejection of one component, most commonly Fomalhaut C. By including circumstellar debris in a subset of the simulations, we also show that such an evolution usually does not disrupt the coherently eccentric debris disk around Fomalhaut A, and in some cases can even produce such a disk. We also find that the final eccentricity of the disk around A and the disk around C are correlated, which may indicate that the dynamics of the three stars stirred C's disk, explaining its unusual brightness.

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)

We consider the observational signatures of giant impacts between planetary embryos. While the debris released in the impact remains in a clump for only a single orbit, there is a much longer lasting asymmetry caused by the fact that all debris must pass through the collision-point. The resulting asymmetry is stationary, it does not orbit the star. The debris is concentrated in a clump at the collision-point, with a more diffuse structure on the opposite side. The asymmetry lasts for typically around 1000 orbital periods of the progenitor, which can be several Myr at distances of ~50 AU. We describe how the appearance of the asymmetric disc depends on the mass and eccentricity of the progenitor, as well as viewing orientation. The wavelength of observation, which determines the grain sizes probed, is also important. Notably, the increased collision rate of the debris at the collision-point makes this the dominant production site for any secondary dust and gas created. For dust small enough to be removed by radiation pressure, and gas with a short lifetime, this causes their distribution to resemble a jet emanating from the (stationary) collision-point. We suggest that the asymmetries seen at large separations in some debris discs, like Beta Pictoris, could be the result of giant impacts. If so this would indicate that planetary embryos are present and continuing to grow at several tens of AU at ages of up to tens of Myr.

Main sequence stars, like the Sun, are often found to be orbited by circumstellar material that can be categorized into two groups, planets and debris. The latter is made up of asteroids and comets, as well as the dust and gas derived from them, which makes debris disks observable in thermal emission or scattered light. These disks may persist over Gyrs through steady-state evolution and/or may also experience sporadic stirring and major collisional breakups, rendering them atypically bright for brief periods of time. Most interestingly, they provide direct evidence that the physical processes (whatever they may be) that act to build large oligarchs from micron-sized dust grains in protoplanetary disks have been successful in a given system, at least to the extent of building up a significant planetesimal population comparable to that seen in the Solar System's asteroid and Kuiper belts. Such systems are prime candidates to host even larger planetary bodies as well. The recent growth in interest in debris disks has been driven by observational work that has provided statistics, resolved images, detection of gas in debris disks, and discoveries of new classes of objects. The interpretation of this vast and expanding dataset has necessitated significant advances in debris disk theory, notably in the physics of dust produced in collisional cascades and in the interaction of debris with planets. Application of this theory has led to the realization that such observations provide a powerful diagnostic that can be used not only to refine our understanding of debris disk physics, but also to challenge our understanding of how planetary systems form and evolve.

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.

This paper explores how the stochastic accretion of planetesimals onto white dwarfs would be manifested in observations of their atmospheric pollution. Archival observations of pollution levels for unbiased samples of DA and non-DA white dwarfs are used to derive the distribution of accretion rates, confirming that rates become systematically lower as sinking time is decreased, with no discernable dependence on cooling age. The accretion rates expected from planetesimals that are all the same mass (ie, a mono-mass distribution) are explored both analytically and using a Monte Carlo model, quantifying how measured accretion rates inevitably depend on sinking time, since different sinking times probe different times since the last accretion event. However, that dependence is so dramatic that a mono-mass distribution can be excluded. Consideration of accretion from a broad distribution of planetesimal masses uncovers an important conceptual difference: accretion is continuous (rather than stochastic) for planetesimals below a certain mass, and the accretion of such planetesimals determines the rate typically inferred from observations; smaller planetesimals dominate the rates for shorter sinking times. A reasonable fit to the observationally inferred accretion rate distributions is found with model parameters consistent with a collisionally evolved mass distribution up to Pluto-mass, and an underlying accretion rate distribution consistent with that expected from descendants of debris discs of main sequence A stars. With these parameters, while both DA and non-DA white dwarfs accrete from the same broad planetesimal distribution, this model predicts that the pollution seen in DAs is dominated by the continuous accretion of <35km objects, and that in non-DAs by >35km objects (though the dominant size varies between stars by around an order of magnitude from this value)...

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.

HD 166191 has been identified by several studies as hosting a rare and extremely bright warm debris disc with an additional outer cool disc component. However, an alternative interpretation is that the star hosts a disc that is currently in transition between a full gas disc and a largely gas-free debris disc. With the help of new optical to mid-IR spectra and Herschel imaging, we argue that the latter interpretation is supported in several ways: i) we show that HD 166191 is co-moving with the ~4 Myr-old Herbig Ae star HD 163296, suggesting that the two have the same age, ii) the disc spectrum of HD 166191 is well matched by a standard radiative transfer model of a gaseous protoplanetary disc with an inner hole, and iii) the HD 166191 mid-IR silicate feature is more consistent with similarly primordial objects. We note some potential issues with the debris disc interpretation that should be considered for such extreme objects, whose lifetime at the current brightness is mush shorter than the stellar age, or in the case of the outer component requires a mass comparable to the solid component of the Solar nebula. These aspects individually and collectively argue that HD 166191 is a 4-5 Myr old star that hosts a gaseous transition disc. Though it does not argue in favour of either scenario, we find strong evidence for 3-5 um disc variability. We place HD 166191 in context with discs at different evolutionary stages, showing that it is a potentially important object for understanding the protoplanetary to debris disc transition.

Debris discs are commonly detected orbiting main-sequence stars, yet little is known regarding their fate as the star evolves to become a giant. Recent observations of radial velocity detected planets orbiting giant stars highlight this population and its importance for probing, for example, the population of planetary systems orbiting intermediate mass stars. Our Herschel survey observed a subset of the Johnson et al program subgiants, finding that 4/36 exhibit excess emission thought to indicate debris, of which 3/19 are planet-hosting stars and 1/17 are stars with no current planet detections. Given the small numbers involved, there is no evidence that the disc detection rate around stars with planets is different to that around stars without planets. Our detections provide a clear indication that large quantities of dusty material can survive the stars' main-sequence lifetime and be detected on the subgiant branch, with important implications for the evolution of planetary systems and observations of polluted or dusty white dwarfs. Our detection rates also provide an important constraint that can be included in models of debris disc evolution.

We present Herschel far-infrared and submillimeter maps of the debris disk associated with the HR 8799 planetary system. We resolve the outer disk emission at 70, 100, 160 and 250 um and detect the disk at 350 and 500 um. A smooth model explains the observed disk emission well. We observe no obvious clumps or asymmetries associated with the trapping of planetesimals that is a potential consequence of planetary migration in the system. We estimate that the disk eccentricity must be <0.1. As in previous work by Su et al. (2009), we find a disk with three components: a warm inner component and two outer components, a planetesimal belt extending from 100 - 310 AU, with some flexibility (+/- 10 AU) on the inner edge, and the external halo which extends to ~2000 AU. We measure the disk inclination to be 26 +/- 3 deg from face-on at a position angle of 64 deg E of N, establishing that the disk is coplanar with the star and planets. The SED of the disk is well fit by blackbody grains whose semi-major axes lie within the planetesimal belt, suggesting an absence of small grains. The wavelength at which the spectrum steepens from blackbody, 47 +/- 30 um, however, is short compared to other A star debris disks, suggesting that there are atypically small grains likely populating the halo. The PACS longer wavelength data yield a lower disk color temperature than do MIPS data (24 and 70 um), implying two distinct halo dust grain populations.

Increasing numbers of sub-stellar companions are now being discovered via direct imaging. Orbital elements for some of these objects have been derived using star--companion astrometry, and several of these appear to have eccentricities significantly greater than zero. We show that stellar motion caused by an undetected inner body may result in the companion elements derived in such a way being incorrect, which could lead to an overestimation of the eccentricity. The magnitude of this effect is quantified in several regimes and we derive the maximum eccentricity error a third body could introduce in a general form, which may be easily applied to any imaged system. Criteria for identifying systems potentially susceptible to this scenario are presented, and we find that around half of the planets/companion brown dwarfs currently imaged could be liable to these errors when their orbital elements are derived. In particular, this effect could be relevant for systems within 100 pc with companions at >50 AU, if they also harbour an unseen $\sim$ 10 Jupiter mass object at $\sim$ 10 AU. We use the Fomalhaut system as an example and show that a 10% error could be induced on the planet's eccentricity by an observationally allowed inner mass, which is similar in size to the current error from astrometry.

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.

Our improving understanding of the life cycle of planetary systems prompts investigations of the role of the Galactic environment before, during and after Asymptotic Giant Branch (AGB) stellar evolution. Here, we investigate the interplay between stellar mass loss, Galactic tidal perturbations, and stellar flybys for evolving stars which host one planet, smaller body or stellar binary companion and reside in the Milky Way's bulge or disc. We find that the potential evolutionary pathways from a main sequence (MS) to a white dwarf (WD) planetary system are a strong function of Galactocentric distance only with respect to the prevalence of stellar flybys. Planetary ejection and collision with the parent star should be more common towards the bulge. At a given location anywhere in the Galaxy, if the mass loss is adiabatic, then the secondary is likely to avoid close flybys during AGB evolution, and cannot eventually escape the resulting WD because of Galactic tides alone. Partly because AGB mass loss will shrink a planetary system's Hill ellipsoid axes by about 20 to 40 per cent, Oort clouds orbiting WDs are likely to be more depleted and dynamically excited than on the MS.

The Debris Disk Explorer (DDX) is a proposed balloon-borne investigation of debris disks around nearby stars. Debris disks are analogs of the Asteroid Belt (mainly rocky) and Kuiper Belt (mainly icy) in our Solar System. DDX will measure the size, shape, brightness, and color of tens of disks. These measurements will enable us to place the Solar System in context. By imaging debris disks around nearby stars, DDX will reveal the presence of perturbing planets via their influence on disk structure, and explore the physics and history of debris disks by characterizing the size and composition of disk dust. The DDX instrument is a 0.75-m diameter off-axis telescope and a coronagraph carried by a stratospheric balloon. DDX will take high-resolution, multi-wavelength images of the debris disks around tens of nearby stars. Two flights are planned; an overnight test flight within the United States followed by a month-long science flight launched from New Zealand. The long flight will fully explore the set of known debris disks accessible only to DDX. It will achieve a raw contrast of 10^-7, with a processed contrast of 10^-8. A technology benefit of DDX is that operation in the near-space environment will raise the Technology Readiness Level of internal coronagraphs, deformable mirrors, and wavefront sensing and control, all potentially needed for a future space-based telescope for high-contrast exoplanet imaging.

Recent results suggest that the two planets in the HD 82943 system are inclined to the sky plane by 20 +/- 4deg. Here, we show that the debris disk in this system is inclined by 27 +/- 4deg, thus adding strength to the derived planet inclinations and suggesting that the planets and debris disk are consistent with being aligned at a level similar to the Solar System. Further, the stellar equator is inferred to be inclined by 28 +/- 4deg, suggesting that the entire star - planet - disk system is aligned, the first time such alignment has been tested for radial velocity discovered planets on ~AU wide orbits. We show that the planet-disk alignment is primordial, and not the result of planetary secular perturbations to the disk inclination. In addition, we note three other systems with planets at >10AU discovered by direct imaging that already have good evidence of alignment, and suggest that empirical evidence of system-wide star - planet - disk alignment is therefore emerging, with the exception of systems that host hot Jupiters. While this alignment needs to be tested in a larger number of systems, and is perhaps unsurprising, it is a reminder that the system should be considered as a whole when considering the orientation of planetary orbits.

We present the first characterisation of the 12um warm dust ("exo-Zodi") luminosity function around Sun-like stars, focussing on the dustiest systems that can be identified by WISE. We detect six new warm dust candidates, five of which have unknown ages. We show that the dustiest old (>Gyr) systems like BD+20 307 are 1 in 10,000 occurrences. Bright warm dust is more common around young (<120Myr) systems, with a ~1% occurrence rate. We show that a two component in situ model where all stars have initially massive warm disks and in which warm debris is also generated at some random time along the stars' main-sequence lifetime, perhaps due to a collision, can explain the observations. However, if all stars only have initially massive warm disks these would not be visible at Gyr ages, and random collisions on the main-sequence are too infrequent to explain the high disk occurrence rate for young stars. That is, neither component can explain the observations on their own. Despite these conclusions, we cannot rule out an alternative model in which comets are scattered from outer regions because the distribution of systems with the appropriate dynamics is unknown. Our in situ model predicts the fraction of stars with exo-Zodi bright enough to cause problems for future exo-Earth imaging attempts is at least ~10%, and is higher for populations of stars younger than a few Gyr. This prediction of ~10% applies to old stars because bright systems like BD+20 307 imply a population of fainter systems that were once bright but are now decaying through fainter levels. Our prediction should be strongly tested by the Large Binocular Telescope Interferometer, providing valuable input for more detailed evolution models. A detection fraction lower than our prediction could indicate that the hot dust in systems like BD+20 307 has a cometary origin due to the quirks of the planetary dynamics.

We present Herschel spatially resolved images of the debris disc orbiting the subgiant Kappa CrB. Not only are these the first resolved images of a debris disc orbiting a subgiant, but Kappa CrB is a rare example of an intermediate mass star where a detailed study of the structure of the planetary system can be made, including both planets and planetesimal belt(s). The only way to discover planets around such stars using the radial velocity technique is to observe 'retired' A stars, which are cooler and slower rotators compared to their main-sequence counterparts. A planetary companion has already been detected orbiting the subgiant Kappa CrB, with revised parameters of m sin i = 2.1MJ and apl = 2.8AU (Johnson et al. 2008a). We present additional Keck I HIRES radial velocity measurements that provide evidence for a second planetary companion, alongside Keck II AO imaging that places an upper limit on the mass of this companion. Modelling of our Herschel images shows that the dust is broadly distributed, but cannot distinguish between a single wide belt (from 20 to 220AU) or two narrow dust belts (at around 40 and 165AU). Given the existence of a second planetary companion beyond approximately 3AU it is possible that the absence of dust within approximately 20AU is caused by dynamical depletion, although the observations are not inconsistent with depletion of these regions by collisional erosion, which occurs at higher rates closer to the star.

Exoplanets have been observed at many stages of their host star's life, including the main sequence (MS), subgiant and red giant branch stages. Also, polluted white dwarfs (WDs) likely represent dynamically active systems at late times. Here, we perform 3-body simulations which include realistic post-MS stellar mass loss and span the entire lifetime of exosystems with two massive planets, from the endpoint of formation to several Gyr into the WD phase of the host star. We find that both MS and WD systems experience ejections and star-planet collisions (Lagrange instability) even if the planet-planet separation well-exceeds the analytical orbit-crossing (Hill instability) boundary. Consequently, MS-stable planets do not need to be closely-packed to experience instability during the WD phase. This instability may pollute the WD directly through collisions, or, more likely, indirectly through increased scattering of smaller bodies such as asteroids or comets. Our simulations show that this instability occurs predominately between tens of Myr to a few Gyrs of WD cooling.

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.

Spectral modeling of the large infrared excess in the Spitzer IRS spectra of HD 172555 suggests that there is more than 10^19 kg of sub-micron dust in the system. Using physical arguments and constraints from observations, we rule out the possibility of the infrared excess being created by a magma ocean planet or a circumplanetary disk or torus. We show that the infrared excess is consistent with a circumstellar debris disk or torus, located at approximately 6 AU, that was created by a planetary scale hypervelocity impact. We find that radiation pressure should remove submicron dust from the debris disk in less than one year. However, the system's mid-infrared photometric flux, dominated by submicron grains, has been stable within 4 percent over the last 27 years, from IRAS (1983) to WISE (2010). Our new spectral modeling work and calculations of the radiation pressure on fine dust in HD 172555 provide a self-consistent explanation for this apparent contradiction. We also explore the unconfirmed claim that 10^47 molecules of SiO vapor are needed to explain an emission feature at 8 um in the Spitzer IRS spectrum of HD 172555. We find that unless there are 10^48 atoms or 0.05 Earth masses of atomic Si and O vapor in the system, SiO vapor should be destroyed by photo-dissociation in less than 0.2 years. We argue that a second plausible explanation for the 8 um feature can be emission from solid SiO, which naturally occurs in submicron silicate "smokes" created by quickly condensing vaporized silicate.

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.

We describe a search for infra-red excess emission from dusty circumstellar material around 180,000 stars observed by the Kepler and WISE missions. This study is motivated by i) the potential to find bright warm disks around planet host stars, ii) a need to characterise the distribution of rare warm disks, and iii) the possible identification of candidates for discovering transiting dust concentrations. We find about 8,000 stars that have excess emission, mostly at 12um. The positions of these stars correlate with the 100um background level so most of the flux measurements associated with these excesses are spurious. We identify 271 stars with plausible excesses by making a 5MJy/sr cut in the IRAS 100um emission. The number counts of these excesses, at both 12 and 22um, have the same distribution as extra-Galactic number counts. Thus, although some excesses may be circumstellar, most can be explained as chance alignments with background galaxies. The one exception is a 22um excess associated with a relatively nearby A-type star that we were able to confirm because the disk occurrence rate is independent of stellar distance. Despite our low detection rate, these results place valuable upper limits on the distribution of large mid-infrared excesses; e.g. fewer than 1:1000 stars have 12um excesses (F_ obs/F_star) larger than a factor of five. In contrast to previous studies, we find no evidence for disks around 1790 stars with candidate planets (we attribute one significant 12um excess to a background galaxy), and no evidence that the disk distribution around planet hosts is different to the bulk population. Higher resolution imaging of stars with excesses is the best way to rule out galaxy confusion and identify more reliable disk candidates among Kepler stars. A similar survey to ours that focusses on nearby stars would be well suited to finding the distribution of rare warm disks.

We study the evolution of debris created in the giant impacts expected during the final stages of terrestrial planet formation. The starting point is the debris created in a simulation of the Moon-forming impact. The dynamical evolution is followed for 10 Myr including the effects of Earth, Venus, Mars and Jupiter. The spatial distribution evolves from a clump in the first few months to an asymmetric ring for the first 10 kyr and finally becoming an axisymmetric ring by about 1 Myr after the impact. By 10 Myr after the impact 20% of the particles have been accreted onto Earth and 17% onto Venus, with 8% ejected by Jupiter and other bodies playing minor roles. However, the fate of the debris also depends strongly on how fast it is collisionally depleted, which depends on the poorly constrained size distribution of the impact debris. Assuming that the debris is made up of 30% by mass mm-cm-sized vapour condensates and 70% boulders up to 500 km, we find that the condensates deplete rapidly on ~1000 yr timescales, whereas the boulders deplete predominantly dynamically. By considering the luminosity of dust produced in collisions within the boulder-debris distribution we find that the Moon-forming impact would have been readily detectable around other stars in Spitzer 24 micron surveys for around 25 Myr after the impact, with levels of emission comparable to many known hot dust systems. The vapour condensates meanwhile produce a short-lived, optically thick, spike of emission. We use these surveys to make an estimate of the fraction of stars that form terrestrial planets, F_TPF. Since current terrestrial planet formation models invoke multiple giant impacts, the low fraction of 10-100 Myr stars found to have warm (~150 K) dust implies that F_TPF ~<10%.

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.

A re-evaluation of time-averaged accretion rates at DBZ-type white dwarfs points to historical, time-averaged rates significantly higher than the currently observed episodes at their DAZ counterparts. The difference between the ongoing, instantaneous accretion rates witnessed at DAZ white dwarfs, which often exceed 1e8 g/s, and those inferred over the past 1e5-1e6 yr for the DBZ stars can be a few orders of magnitude, and therefore must result from high-rate episodes of tens to hundreds of years so they remain undetected to date. This paper explores the likelihood that such brief, intense accretion episodes of gas-phase material can account for existing data. For reasonable assumptions about the circumstellar gas, accretion rates approaching or exceeding 1e15 g/s are possible, similar to rates observed in quiescent cataclysmic variables, and potentially detectable with future x-ray missions or wide-field monitoring facilities. Gaseous debris that is prone to such rapid accretion may be abundant immediately following a tidal disruption event via collisions and sublimation, or if additional bodies impinge upon an extant disk. Particulate disk matter accretes at or near the Poynting-Robertson drag rate for long periods between gas-producing events, consistent with rates inferred for dusty DAZ white dwarfs. In this picture, warm DAZ stars without infrared excesses have rates consistent with accretion from particulate disks that remain undetected. This overall picture has implications for quasi-steady state models of accretion and the derived chemical composition of asteroidal debris in DBZ white dwarfs.

SPICES (Spectro-Polarimetric Imaging and Characterization of Exoplanetary Systems) is a five-year M-class mission proposed to ESA Cosmic Vision. Its purpose is to image and characterize long-period extrasolar planets and circumstellar disks in the visible (450 - 900 nm) at a spectral resolution of about 40 using both spectroscopy and polarimetry. By 2020/22, present and near-term instruments will have found several tens of planets that SPICES will be able to observe and study in detail. Equipped with a 1.5 m telescope, SPICES can preferentially access exoplanets located at several AUs (0.5-10 AU) from nearby stars ($<$25 pc) with masses ranging from a few Jupiter masses to Super Earths ($\sim$2 Earth radii, $\sim$10 M$_{\oplus}$) as well as circumstellar disks as faint as a few times the zodiacal light in the Solar System.

We present new MIDI interferometric and VISIR spectroscopic observations of HD113766 and HD172555. Additionally we present VISIR 11um and 18um imaging observations of HD113766. These sources represent the youngest (16Myr and 12Myr old respectively) debris disc hosts with emission on <<10AU scales. We find that the disc of HD113766 is partially resolved on baselines of 42-102m, with variations in resolution with baseline length consistent with a Gaussian model for the disc with FWHM of 1.2-1.6AU (9-12mas). This is consistent with the VISIR observations which place an upper limit of 0."14 (17AU) on the emission, with no evidence for extended emission at larger distances. For HD172555 the MIDI observations are consistent with complete resolution of the disc emission on all baselines of lengths 56-93m, putting the dust at a distance of >1AU (>35mas). When combined with limits from TReCS imaging the dust at ~10um is constrained to lie somewhere in the region 1-8AU. Observations at ~18um reveal extended disc emission which could originate from the outer edge of a broad disc, the inner parts of which are also detected but not resolved at 10um, or from a spatially distinct component. These observations provide the most accurate direct measurements of the location of dust at 1-8AU that might originate from the collisions expected during terrestrial planet formation. These observations provide valuable constraints for models of the composition of discs at this epoch and provide a foundation for future studies to examine in more detail the morphology of debris discs.

We present models for the formation of terrestrial planets, and the collisional evolution of debris disks, in planetary systems that contain multiple unstable gas giants. We previously showed that the dynamics of the giant planets introduces a correlation between the presence of terrestrial planets and debris disks. Here we present new simulations that show that this connection is qualitatively robust to changes in: the mass distribution of the giant planets, the width and mass distribution of the outer planetesimal disk, and the presence of gas in the disk. We discuss how variations in these parameters affect the evolution. Systems with equal-mass giant planets undergo the most violent instabilities, and these destroy both terrestrial planets and the outer planetesimal disks that produce debris disks. In contrast, systems with low-mass giant planets efficiently produce both terrestrial planets and debris disks. A large fraction of systems with low-mass outermost giant planets have stable gaps between these planets that are frequently populated by planetesimals. Planetesimal belts between outer giant planets may affect debris disk SEDs. If Earth-mass seeds are present in outer planetesimal disks, the disks radially spread to colder temperatures. We argue that this may explain the very low frequency of > 1 Gyr-old solar-type stars with observed 24 micron excesses. Among the (limited) set of configurations explored, the best candidates for hosting terrestrial planets at ~1 AU are stars older than 0.1-1 Gyr with bright debris disks at 70 micron but with no currently-known giant planets. These systems combine evidence for rocky building blocks, with giant planet properties least likely to undergo destructive dynamical evolution. We predict an anti-correlation between debris disks and eccentric giant planets, and a positive correlation between debris disks and terrestrial planets.

The Sun will eventually lose about half of its current mass nonlinearly over several phases of post-main sequence evolution. This mass loss will cause any surviving orbiting body to increase its semimajor axis and perhaps vary its eccentricity. Here, we use a range of Solar models spanning plausible evolutionary sequences and assume isotropic mass loss to assess the possibility of escape from the Solar System. We find that the critical semimajor axis in the Solar System within which an orbiting body is guaranteed to remain bound to the dying Sun due to perturbations from stellar mass loss alone is approximately 1,000 AU - 10,000 AU. The fate of objects near or beyond this critical semimajor axis, such as the Oort Cloud, outer scattered disc and specific bodies such as Sedna, will significantly depend on their locations along their orbits when the Sun turns off of the main sequence. These results are applicable to any exoplanetary system containing a single star with a mass, metallicity and age which are approximately equal to the Sun's, and suggest that few extrasolar Oort Clouds could survive post-main sequence evolution intact.

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 scattering of small bodies by planets is an important dynamical process in planetary systems. We present an analytical model to describe this process using the simplifying assumption that each particle's dynamics is dominated by a single planet at a time. As such the scattering process can be considered as a series of three body problems during each of which the Tisserand parameter with respect to the relevant planet is conserved. This constrains the orbital parameter space into which a particle can be scattered. Such arguments have previously been applied to the process by which comets are scattered to the inner Solar System from the Kuiper belt. Our analysis generalises this for an arbitrary planetary system. For particles scattered from an outer belt directly along a chain of planets, based on the initial value of the Tisserand parameter, we find that it is possible to (i) determine which planets can eject the particles from the system, (ii) define a minimum stellar distance to which particles can be scattered, and (iii) constrain range of particle inclinations (and hence the disc height) at different distances. Applying this to the Solar System, we determine that the planets are close to optimally separated for scattering particles between them. Concerning warm dust found around stars that also have Kuiper belt analogues, we show that, if there is to be a dynamical link between the outer and inner regions, then certain architectures for the intervening planetary system are incapable of producing the observations. Furthermore we show that for certain planetary systems, comets can be scattered from an outer belt, or with fewer constraints, from an Oort cloud analogue, onto star-grazing orbits, in support of a planetary origin to the metal pollution and dustiness of some nearby white dwarfs.

We have analyzed Spitzer and NASA/IRTF 2 - 35 \mum spectra of the warm, ~350 K circumstellar dust around the nearby MS star {\eta} Corvi (F2V, 1.4 \pm 0.3 Gyr). The spectra show clear evidence for warm, water- and carbon-rich dust at ~3 AU from the central star, in the system's Terrestrial Habitability Zone. Spectral features due to ultra-primitive cometary material were found, in addition to features due to impact produced silica and high temperature carbonaceous phases. At least 9 x 10^18 kg of 0.1 - 100 \mum warm dust is present in a collisional equilibrium distribution with dn/da ~ a^-3.5, the equivalent of a 130 km radius KBO of 1.0 g/cm^3 density and similar to recent estimates of the mass delivered to the Earth at 0.6 - 0.8 Gyr during the Late Heavy Bombardment. We conclude that the parent body was a Kuiper-Belt body or bodies which captured a large amount of early primitive material in the first Myrs of the system's lifetime and preserved it in deep freeze at ~150 AU. At ~1.4 Gyr they were prompted by dynamical stirring of their parent Kuiper Belt into spiraling into the inner system, eventually colliding at 5-10 km/sec with a rocky planetary body of mass \leq MEarth at ~3 AU, delivering large amounts of water (>0.1% of MEarth's Oceans) and carbon-rich material. The Spitzer spectrum also closely matches spectra reported for the Ureilite meteorites of the Sudan Almahata Sitta fall in 2008, suggesting that one of the Ureilite parent bodies was a KBO.

The orbit of a planet is surrounded by a chaotic zone wherein nearby particles' orbits are chaotic and unstable. Wisdom (1980) showed that the chaos is driven by the overlap of mean motion resonances which occurs within a distance (da/a)_chaos = 1.3 mu^2/7 of the planet's orbit. However, the width of mean motion resonances grows with the particles' eccentricity, which will increase the width of the chaotic zone at higher eccentricities. Here we investigate the width of the chaotic zone using the iterated encounter map and N-body integrations. We find that the classical prescription works well for particles on low-eccentricity orbits. However, above a critical eccentricity, dependent upon the mass of the planet, the width of the chaotic zone increases with eccentricity. An extension of Wisdom's analytical arguments then shows that, above the critical eccentricity, the chaotic zone width is given by (da/a)_chaos = 1.8 e^1/5 mu^1/5, which agrees well with the encounter map results. The critical eccentricity is given by e_cr = 0.21 mu^3/7. This extended chaotic zone results in a larger cleared region when a planet sculpts the inner edge of a debris disc comprised of eccentric planetesimals. Hence, the planet mass estimated from the classical chaotic zone may be erroneous. We apply this result to the HR8799 system, showing that the masses of HR8799b inferred from the truncation of the disc may vary by up to 50% depending on the disc particles' eccentricities. With a disc edge at 90AU, the necessary mass of planet b to cause the truncation would be 8-10 Jovian masses if the disc particles have low eccentricities (<0.02), but only 4-8 Jovian masses if the disc particles have higher eccentricities. Our result also has implications for the ability of a planet to feed material into an inner system, a process which may explain metal pollution in White Dwarf atmospheres.

The main-sequence solar-type star HD69830 has an unusually large amount of dusty debris orbiting close to three planets found via the radial velocity technique. In order to explore the dynamical interaction between the dust and planets, we have performed multi-epoch photometry and spectroscopy of the system over several orbits of the outer dust. We find no evidence for changes in either the dust amount or its composition, with upper limits of 5-7% (1 $\sigma$ per spectral element) on the variability of the {\it dust spectrum} over 1 year, 3.3% (1 $\sigma$) on the broad-band disk emission over 4 years, and 33% (1 $\sigma$) on the broad-band disk emission over 24 years. Detailed modeling of the spectrum of the emitting dust indicates that the dust is located outside of the orbits of the three planets and has a composition similar to main-belt, C-type asteroids asteroids in our solar system. Additionally, we find no evidence for a wide variety of gas species associated with the dust. Our new higher SNR spectra do not confirm our previously claimed detection of H$_2$O ice leading to a firm conclusion that the debris can be associated with the break-up of one or more C-type asteroids formed in the dry, inner regions of the protoplanetary disk of the HD69830 system. The modeling of the spectral energy distribution and high spatial resolution observations in the mid-infrared are consistent with a $\sim$ 1 AU location for the emitting material.

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.

We describe a search for dust created in collisions between the Saturnian irregular satellites using archival \emph{Spitzer} MIPS observations. Although we detected a degree scale Saturn-centric excess that might be attributed to an irregular satellite dust cloud, we attribute it to the far-field wings of the PSF due to nearby Saturn. The Spitzer PSF is poorly characterised at such radial distances, and we expect PSF characterisation to be the main issue for future observations that aim to detect such dust. The observations place an upper limit on the level of dust in the outer reaches of the Saturnian system, and constrain how the size distribution extrapolates from the smallest known (few km) size irregulars down to micron-size dust. Because the size distribution is indicative of the strength properties of irregulars, we show how our derived upper limit implies irregular satellite strengths more akin to comets than asteroids. This conclusion is consistent with their presumed capture from the outer regions of the Solar System.

Mounting discoveries of extrasolar planets orbiting post-main sequence stars motivate studies aimed at understanding the fate of these planets. In the traditional "adiabatic" approximation, a secondary's eccentricity remains constant during stellar mass loss. Here, we remove this approximation, investigate the full two-body point-mass problem with isotropic mass loss, and illustrate the resulting dynamical evolution. The magnitude and duration of a star's mass loss combined with a secondary's initial orbital characteristics might provoke ejection, modest eccentricity pumping, or even circularisation of the orbit. We conclude that Oort clouds and wide-separation planets may be dynamically ejected from 1-7 Solar-mass parent stars during AGB evolution. The vast majority of planetary material which survives a supernova from a 7-20 Solar-mass progenitor will be dynamically ejected from the system, placing limits on the existence of first-generation pulsar planets. Planets around >20 Solar-mass black hole progenitors may easily survive or readily be ejected depending on the core collapse and superwind models applied. Material ejected during stellar evolution might contribute significantly to the free-floating planetary population.

The eccentric orbits of the known extrasolar giant planets provide evidence that most planet-forming environments undergo violent dynamical instabilities. Here, we numerically simulate the impact of giant planet instabilities on planetary systems as a whole. We find that populations of inner rocky and outer icy bodies are both shaped by the giant planet dynamics and are naturally correlated. Strong instabilities -- those with very eccentric surviving giant planets -- completely clear out their inner and outer regions. In contrast, systems with stable or low-mass giant planets form terrestrial planets in their inner regions and outer icy bodies produce dust that is observable as debris disks at mid-infrared wavelengths. Fifteen to twenty percent of old stars are observed to have bright debris disks (at wavelengths of ~70 microns) and we predict that these signpost dynamically calm environments that should contain terrestrial planets.

Circumstantial evidence suggests that most known extra-solar planetary systems are survivors of violent dynamical instabilities. Here we explore how giant planet instabilities affect the formation and survival of terrestrial planets. We simulate planetary system evolution around Sun-like stars from initial conditions that comprise: an inner disk of planetesimals and planetary embryos, three giant planets at Jupiter-Saturn distances, and a massive outer planetesimal disk. We then calculate dust production rates and debris disk SEDs assuming that each planetesimal particle represents an ensemble of smaller bodies in collisional equilibrium. We predict a strong correlation between the presence of terrestrial planets and debris disks, mediated by the giant planets. Strong giant planet instabilities destroy all rocky material - including fully-formed terrestrial planets if the instabilities occur late - along with the icy planetesimals. Stable or weakly unstable systems allow terrestrial planets to accrete and significant dust to be produced in their outer regions. Stars older than ~100 Myr with bright cold dust emission (at ~70 microns) signpost the dynamically calm environments conducive to efficient terrestrial accretion. We predict that while the typical eccentricities of terrestrial planets are small, there should exist a novel class of terrestrial planet system whose single planet undergoes large amplitude oscillations in eccentricity and inclination. By scaling to the observed semimajor axis distribution of giant exoplanets, we estimate that terrestrial exoplanets in the same systems should be a few times more abundant at 0.5 AU than giant or terrestrial exoplanets at 1 AU. Finally, we discuss the Solar System, which appears to be unusual in combining a rich terrestrial planet system with a low dust content.

We present a new scheme for determining the shape of the size distribution, and its evolution, for collisional cascades of planetesimals undergoing destructive collisions and loss processes like Poynting-Robertson drag. The scheme treats the steady state portion of the cascade by equating mass loss and gain in each size bin; the smallest particles are expected to reach steady state on their collision timescale, while larger particles retain their primordial distribution. For collision-dominated disks, steady state means that mass loss rates in logarithmic size bins are independent of size. This prescription reproduces the expected two phase size distribution, with ripples above the blow-out size, and above the transition to gravity-dominated planetesimal strength. The scheme also reproduces the expected evolution of disk mass, and of dust mass, but is computationally much faster than evolving distributions forward in time. For low-mass disks, P-R drag causes a turnover at small sizes to a size distribution that is set by the redistribution function (the mass distribution of fragments produced in collisions). Thus information about the redistribution function may be recovered by measuring the size distribution of particles undergoing loss by P-R drag, such as that traced by particles accreted onto Earth. Although cross-sectional area drops with 1/age^2 in the PR-dominated regime, dust mass falls as 1/age^2.8, underlining the importance of understanding which particle sizes contribute to an observation when considering how disk detectability evolves. Other loss processes are readily incorporated; we also discuss generalised power law loss rates, dynamical depletion, realistic radiation forces and stellar wind drag.

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.

A quarter of DA white dwarfs are metal polluted, yet elements heavier than helium sink down through the stellar atmosphere on timescales of days. Hence, these white dwarfs must be currently accreting material containing heavy elements. Here, we consider whether the scattering of comets or asteroids from an outer planetary system, following stellar mass loss on the asymptotic giant branch, can reproduce these observations. We use N-body simulations to investigate the effects of stellar mass loss on a simple system consisting of a planetesimal belt whose inner edge is truncated by a planet. Our simulations find that, starting with a planetesimal belt population fitted to the observed main sequence evolution, sufficient mass is scattered into the inner planetary system to explain the inferred heavy element accretion rates. This assumes that some fraction of the mass scattered into the inner planetary system ends up on star-grazing orbits, is tidally disrupted and accreted onto the white dwarf. The simulations also reproduce the observed decrease in accretion rate with cooling age and predict accretion rates in old (>1Gyr) white dwarfs, in line with observations. The efficiency we assumed for material scattered into the inner planetary system to end up on star-grazing orbits is based on a Solar-like planetary system, since the simulations show that a single planet is not sufficient. Although the correct level of accretion is reproduced, the simulations predict a higher fraction of accreting white dwarfs than observed. This could indicate that evolved planetary systems are less efficient at scattering bodies onto star-grazing orbits or that dynamical instabilities post-stellar mass loss cause rapid planetesimal belt depletion for a significant fraction of systems.

Mean motion resonances are a common feature of both our own Solar System and of extrasolar planetary systems. Bodies can be trapped in resonance when their orbital semi-major axes change, for instance when they migrate through a protoplanetary disc. We use a Hamiltonian model to thoroughly investigate the capture behaviour for first and second order resonances. Using this method, all resonances of the same order can be described by one equation, with applications to specific resonances by appropriate scaling. We focus on the limit where one body is a massless test particle and the other a massive planet. We quantify how the the probability of capture into a resonance depends on the relative migration rate of the planet and particle, and the particle's eccentricity. Resonant capture fails for high migration rates, and has decreasing probability for higher eccentricities. More massive planets can capture particles at higher eccentricities and migration rates. We also calculate libration amplitudes and the offset of the libration centres for captured particles, and the change in eccentricity if capture does not occur. Libration amplitudes are higher for larger initial eccentricity. The model allows for a complete description of a particle's behaviour as it successively encounters several resonances. We discuss implications for several scenarios: (i) Planet migration through gas discs trapping other planets or planetesimals in resonances. (ii) Planet migration through a debris disc. (iii) Dust migration through PR drag. The Hamiltonian model will allow quick interpretation of the resonant properties of extrasolar planets and Kuiper Belt Objects, and will allow synthetic images of debris disc structures to be quickly generated, which will be useful for predicting and interpreting disc images made with ALMA, Darwin/TPF or similar missions. [Abridged]

The population of debris discs on the main sequence is well constrained, however very little is known about debris discs around evolved stars. In this work we provide a theoretical framework that considers the effects of stellar evolution on debris discs; firstly considering the evolution of an individual disc from the main sequence through to the white dwarf phase, then extending this to the known population of debris discs around main sequence A stars. It is found that discs around evolved stars are harder to detect than on the main sequence. In the context of our models discs should be detectable with Herschel or Alma on the giant branch, subject to the uncertain effect of sublimation on the discs. The best chances are for hot young white dwarfs, fitting nicely with the observations e.g the helix nebula (Su et al. 2007) and 9 systems presented by Chu & Bilikova.

Mean motion resonances are a common feature of both our own Solar System and of extrasolar planetary systems. Bodies can be trapped in resonance when their orbital semi-major axes change, for instance when they migrate through a protoplanetary disc. We use a Hamiltonian model to thoroughly investigate the capture behaviour for first and second order resonances. Using this method, all resonances of the same order can be described by one equation, with applications to specific resonances by appropriate scaling. We focus on the limit where one body is a massless test particle and the other a massive planet. We quantify how the the probability of capture into a resonance depends on the relative migration rate of the planet and particle, and the particle's eccentricity. Resonant capture fails for high migration rates, and has decreasing probability for higher eccentricities, although for certain migration rates, capture probability peaks at a finite eccentricity. We also calculate libration amplitudes and the offset of the libration centres for captured particles, and the change in eccentricity if capture does not occur. Libration amplitudes are higher for larger initial eccentricity. The model allows for a complete description of a particle's behaviour as it successively encounters several resonances. The model is applicable to many scenarios, including (i) Planet migration through gas discs trapping other planets or planetesimals in resonances; (ii) Planet migration through a debris disc; (iii) Dust migration through PR drag. Full details can be found in \cite{2010submitted}. (Abridged)

Since the 1980's it has been becoming increasingly clear that the Solar System's irregular satellites are collisionally evolved. We derive a general model for the collisional evolution of an irregular satellite swarm and apply it to the Solar System and extrasolar planets. Our model reproduces the Solar System's complement of observed irregulars well, and suggests that the competition between grain-grain collisions and Poynting-Robertson (PR) drag helps set the fate of the dust. Because swarm collision rates decrease over time the main dust sink can change with time, and may help unravel the accretion history of synchronously rotating regular satellites that show brightness asymmetries. Some level of dust must be present on AU scales around the Solar System's giant planets, which we predict may be at detectable levels. We also predict whether dust produced by extrasolar circumplanetary swarms can be detected. The coronagraphic instruments on JWST will have the ability to detect the dust generated by these swarms, which are most detectable around planets that orbit at tens of AU from the youngest stars. Because the collisional decay of swarms is relatively insensitive to planet mass, swarms can be much brighter than their host planets and allow discovery of Neptune-mass planets that would otherwise remain invisible. This dust may have already been detected. The observations of the planet Fomalhaut b can be explained as scattered light from dust produced by the collisional decay of an irregular satellite swarm around a 10 Earth-mass planet. Such a swarm comprises about 5 Lunar masses worth of irregular satellites. Finally, we consider what happens if Fomalhaut b passes through Fomalhaut's main debris ring, which allows the circumplanetary swarm to be replenished through collisions with ring planetesimals. (abridged)

We have obtained new resolved images of the well-studied HR 4796A dust ring at 18 and 25 microns with the 8-meter Gemini telescopes. These images confirm the previously observed spatial extent seen in mid-IR, near-IR, and optical images of the source. We detect brightness and temperature asymmetries such that dust on the NE side is both brighter and warmer than dust in the SW. We show that models of so-called pericenter glow account for these asymmetries, thus both confirming and extending our previous analyses. In this scenario, the center of the dust ring is offset from the star due to gravitational perturbations of a body with an eccentric orbit that has induced a forced eccentricity on the dust particle orbits. Models with 2-micron silicate dust particles and a forced eccentricity of 0.06 simultaneously fit the observations at both wavelengths. We also show that parameters used to characterize the thermal-emission properties of the disk can also account for the disk asymmetry observed in shorter-wavelength scattered-light images.

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.

Two thirds of the F star members of the 12 Myr old Beta Pictoris Moving Group (BPMG) show significant excess emission in the mid-infrared, several million years after the expected dispersal of the protoplanetary disc. Theoretical models of planet formation suggest that this peak in the mid-infrared emission could be due to the formation of Pluto-sized bodies in the disc, which ignite the collisional cascade and enhance the production of small dust. Here we present resolved mid-infrared imaging of the disc of HD191089 (F5V in the BPMG) and consider its implications for the state of planet formation in this system. HD191089 was observed at 18.3 microns using T-ReCS on Gemini South and the images were compared to models of the disc to constrain the radial distribution of the dust. The emission observed at 18.3 microns is shown to be significantly extended beyond the PSF at a position angle of 80 degrees. This is the first time dust emission has been resolved around HD191089. Modelling indicates that the emission arises from a dust belt from 28-90 AU, inclined at 35 degrees from edge on with very little emission from the inner 28AU of the disc, indicating the presence of an inner cavity. The steep slope of the inner edge is more consistent with truncation by a planet than with ongoing stirring. A tentative brightness asymmetry F(W)/F(E)=0.80+/-0.12 (1.8 sigma) between the two sides of the disc could be evidence for perturbations from a massive body on an eccentric orbit in the system.

While the population of main sequence debris discs is well constrained, little is known about debris discs around evolved stars. This paper provides a theoretical framework considering the effects of stellar evolution on debris discs, particularly the production and loss of dust within them. Here we repeat a steady state model fit to disc evolution statistics for main sequence A stars, this time using realistic grain optical properties, then evolve that population to consider its detectability at later epochs. Our model predicts that debris discs around giant stars are harder to detect than on the main sequence because radiation pressure is more effective at removing small dust around higher luminosity stars. Just 12% of first ascent giants within 100pc are predicted to have discs detectable with Herschel at 160um. However this is subject to the uncertain effect of sublimation on the disc, which we propose can thus be constrained with such observations. Our model also finds that the rapid decline in stellar luminosity results in only very young white dwarfs having luminous discs. As such systems are on average at larger distances they are hard to detect, but we predict that the stellar parameters most likely to yield a disc detection are a white dwarf at 200pc with cooling age of 0.1Myr, in line with observations of the Helix Nebula. Our model does not predict close-in (<0.01AU) dust, as observed for some white dwarfs, however we find that stellar wind drag leaves significant mass (~10^{-2}Msolar), in bodies up to ~10m in diameter, inside the disc at the end of the AGB phase which may replenish these discs.

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.

(Abridged) Studies of debris discs have shown that most systems are analogous to the EKB. In this study we aim to determine how many IRAS 25um excesses towards A stars are real, and investigate where the dust lies. We observe with TIMMI2, VISIR, Michelle and TReCS a sample of A and B-type main sequence stars reported as having mid-IR excess. We constrain the location of the debris through combined modelling of the emission spectrum and a modelling technique designed to constrain the radial extent of emission in mid-IR imaging. We independently confirm the presence of warm dust around 3 of the candidates: HD3003, HD80950 and eta Tel. For the binary HD3003 a stability analysis indicates the dust is either circumstellar and lying at ~4 AU with the binary orbiting at >14AU, or the dust lies in an unstable location; there is some evidence for temporal evolution of its excess emission on a ~20 year timescale. For 7 of the targets we present quantitative limits on the location of dust around the star. We demonstrate that the disc around HD71155 must have spatially distinct components at 2 and 60AU. We model the limits of current instrumentation and show that most of the known A star debris discs which could be readily resolved at 18um on 8m instruments have been resolved. Limits from unresolved imaging can help distinguish between competing models of the disc emission, but resolved imaging is key to the determination of the disc location. Modelling of the detection limits for extended emission can be useful for targeting future observational campaigns. MIRI on the JWST will be able to resolve most of the known A star debris disc population. METIS on the E-ELT will provide the opportunity to explore the hot disc population more thoroughly by detecting extended emission where calibration accuracy limits disc detection through photometry, reaching levels below 1 zodi for stars at <10pc.

This paper considers the evidence that debris disks are self-stirred by the formation of Plutos. A model for the dust produced during self-stirring is applied to statistics for A stars. As there is no significant difference between excesses of A-stars <50Myr old, we focus on reproducing the broad trends, the "rise and fall" of the fraction of stars with excesses. Using a population model, we find that the statistics and trends can be reproduced with a self-stirring model of planetesimal belts with radii distributed between 15-120AU. Disks must have this 15AU minimum radius to show a peak in disk fraction, rather than a monotonic decline. Populations of extended disks with fixed inner and/or outer radii fail to fit the statistics, due mainly to the slow 70um evolution as stirring moves further out in the disk. This conclusion, that debris disks are narrow belts, is independent of the significance of 24um trends for young A-stars. We show that the statistics can also be reproduced with a model in which disks are stirred by secular perturbations from a nearby eccentric planet. Detailed imaging is therefore the best way to characterise the stirring mechanism. From a more detailed look at beta Pictoris Moving Group and TW Hydrae Association A-stars we find that the disk around beta Pictoris is likely the result of secular stirring by the proposed planet at ~10AU; the structure of the HR 4796A disk also points to sculpting by a planet. The two other stars with disks, HR 7012 and eta Tel, possess transient hot dust, though the outer eta Tel disk is consistent with a self-stirred origin. Planet formation provides a natural explanation for the belt-like nature of debris disks, with inner regions cleared by planets that may also stir the disk, and the outer edges set by where planetesimals can form. [abridged]

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 habitability of planets is strongly affected by impacts from comets and asteroids. Indications from the ages of Moon rocks suggest that the inner Solar System experienced an increased rate of impacts roughly 3.8 Gya known as the Late Heavy Bombardment (LHB). Here we develop a model of how the Solar System would have appeared to a distant observer during its history based on the Nice model of Gomes et al. (2005). We compare our results with observed debris discs. We show that the Solar System would have been amongst the brightest of these systems before the LHB. Comparison with the statistics of debris disc evolution shows that such heavy bombardment events must be rare occurring around less than 12% of Sun-like stars.

Models for the steady state collisional evolution of low eccentricity planetesimal belts identify debris disks with hot dust at 1AU, like eta Corvi and HD69830, as anomalous since collisional processing should have removed most of the planetesimal mass over their >1 Gyr lifetimes. This paper looks at the effect of large planetesimal eccentricities (e>>0.3) on their collisional lifetime and the amount of mass that can remain at late times M_{late}. For an axisymmetric planetesimal disk with common pericentres and eccentricities e, we find that M_{late} \propto e^{-5/3}(1+e)^{4/3}(1-e)^{-3}. For a scattered disk-like population (i.e., common pericentres), in the absence of dynamical evolution, the mass evolution at late times would be as if only planetesimals with the largest eccentricity were present. Despite the increased remaining mass, higher eccentricities do not increase the hot emission from the collisional cascade until e>0.99, partly because most collisions occur near pericentre thus increasing the dust blow-out diameter. However, at high eccentricities (e>0.97) the blow-out population extending out from pericentre may be detectable above the collisional cascade; higher eccentricities also increase the probability of witnessing a recent collision. All of the imaging and spectroscopic constraints for eta Corvi can be explained with a single planetesimal population with pericentre at 0.75AU, apocentre at 150AU, and mass 5M_\oplus; however, the origin of such a high eccentricity population remains challenging. The mid-IR excess to HD69830 can be explained by the ongoing destruction of a debris belt produced in a recent collision in an eccentric planetesimal belt, but the lack of far-IR emission requires small bound grains to be absent from the parent planetesimal belt, possibly due to sublimation.

We consider the formation and migration of protoplanetary embryos in disks around the stars in tight binary systems (separations ~ 20 AU. In such systems, the initial stages of runaway embryo formation are expected to only take place within some critical disk radius a_{crit}, due to the perturbing effect of the binary companions (Thebault et al. 2009). We perform n-body simulations of the evolution of such a population of inner-disk embryos surrounded by an outer-disk of smaller planetesimals. Taking Alpha Centauri-B as our fiducial reference example in which a_{crit} ~ 0.7 AU, and using a Minimum Mass Nebular Model with $\Sigma \propto a^{-3/2}$, we find that within 10^6 yrs (10^7 yrs), systems will on average contain embryos which have migrated out to 0.9 AU (1.2 AU), with the average outer-most body having a mass of 0.2 M_{earth} 0.4 M_{earth}. Changes to increase the surface density of solids or to use a flatter profile both produce increased embryo migration and growth. At a given time, the relative change in semi-major axis of the outer-most embryo in these simulations is found to be essentially independent of a_{crit}, and we note that little further embryo migration takes place beyond 10^7 years. We conclude that the suppression of runaway growth outside a_{crit} does not mean that the habitable zones in such tight binary systems will be devoid of detectable, terrestrial mass planets, even if a_{crit} lies significantly interior to the inner edge of the habitable zone.

Detectable debris discs are thought to require dynamical excitation (`stirring'), so that planetesimal collisions release large quantities of dust. We investigate the effects of the secular perturbations of a planet, which may lie at a significant distance from the planetesimal disc, to see if these perturbations can stir the disc, and if so over what time-scale. The secular perturbations cause orbits at different semi-major axes to precess at different rates, and after some time t_cross initially non-intersecting orbits begin to cross. We show that t_cross is proportional to a_disc^(9/2)/(m_pl e_pl a_pl^3), where m_pl, e_pl and a_pl are the mass, eccentricity, and semi-major axis of the planet, and a_disc is the semi-major axis of the disc. This time-scale can be faster than that for the growth of planetesimals to Pluto's size within the outer disc. We also calculate the magnitude of the relative velocities induced amongst planetesimals and infer that a planet's perturbations can typically cause destructive collisions out to 100's of AU. Recently formed planets can thus have a significant impact on planet formation in the outer disc which may be curtailed by the formation of giant planets much closer to the star. The presence of an observed debris disc does not require the presence of Pluto-sized objects within it, since it can also have been stirred by a planet not in the disc. For the star epsilon Eridani, we find that the known RV planet can excite the planetesimal belt at 60 AU sufficiently to cause destructive collisions of bodies up to 100 km in size, on a time-scale of 40 Myr.

Although 70 % of the stars in the Galaxy are M-dwarfs, thermal emission searches for cold debris disks have been conducted mostly for A-type and solar-type stars. We report on new lambda=1.2 mm continuum observations of thirty M-dwarfs, using the MAMBO-2 bolometer array camera at the IRAM 30m telescope. For a statistical analysis, we combine these data with our prior SCUBA and MAMBO-2 observations of 20 other M-dwarfs. Our total sample divides in M-dwarfs in moving groups, with relatively young ages, and in nearby M-dwarfs with unknown ages. Only one cold debris disk (GJ842.2) was detected significantly. We compare the implied disk abundance constraints with those found in two comparable submillimeter surveys of 10 to 190 Myr old A- and FGK-type stars. For the 19 youngest (ages less than 200 Myr) M-dwarfs in our sample, we derive a cold disk fraction of 5.3^{+10.5}_{-5.0} %, compared to 15 +/-11.5 % for FGK-stars and 22^{+33}_{-20} % for A-stars. Hence, for this age group, there is an apparent trend of fewer cold disks for later stellar types, i.e., lower star masses. Although its statistical significance is marginal, this trend is strengthened by the deeper sensitivity of observations in the M-dwarf sample. We derive a cold disk fraction of < 10 % for the older (likely a few Gyr) M-dwarfs in our sample. Finally, although inconclusively related to a debris disk, we present the complex millimeter structure found around the position of the M1.5 dwarf GJ526 in our sample.

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.

The Nice model of Gomes et al. (2005) suggests that the migration of the giant planets caused a planetesimal clearing event which led to the Late Heavy Bombardment (LHB) at 880 Myr. Here we investigate the IR emission from the Kuiper belt during the history of the Solar System as described by the Nice model. We describe a method for easily converting the results of n-body planetesimal simulations into observational properties (assuming black-body grains and a single size distribution) and further modify this method to improve its realism (using realistic grain properties and a three-phase size distribution). We compare our results with observed debris discs and evaluate the plausibility of detecting an LHB-like process in extrasolar systems. Recent surveys have shown that 4% of stars exhibit 24 um excess and 16% exhibit 70 um excess. We show that the Solar System would have been amongst the brightest of these systems before the LHB at both 24 and 70 um. We find a significant increase in 24 um emission during the LHB, which rapidly drops off and becomes undetectable within 30 Myr, whereas the 70 um emission remains detectable until 360 Myr after the LHB. Comparison with the statistics of debris disc evolution shows that such depletion events must be rare occurring around less than 12% of Sun-like stars and with this level of incidence we would expect approximately 1 of the 413 Sun-like, field stars so far detected to have a 24 um excess to be currently going through an LHB. We also find that collisional processes are important in the Solar System before the LHB and that parameters for weak Kuiper belt objects are inconsistent with the Nice model interpretation of the LHB.

Most of the known debris discs exhibit cool dust in regions analogous to the Edgeworth-Kuiper Belt. However, a rare subset show hot excess from within a few AU, which is often inferred to be transient. We examine 2 such sources to place limits on their location to help distinguish between different interpretations for their origin. We use MIDI on the VLTI to observe the debris discs around eta Corvi and HD69830 using baseline lengths from 44-130m. New VISIR observations of HD69830 at 18.7um are also presented. These observations are compared with disc models to place limits on disc size. The visibility functions measured with MIDI for both sources show significant variation with wavelength across 8-13um in a manner consistent with the disc flux being well resolved, notably with a dip at 10-11.5um due to the silicate emission feature. The average ratio of visibilities measured between 10-11.5um and 8-9um is 0.934+/-0.015 for HD69830 and 0.880+/-0.013 for eta Corvi over the 4 baselines for each source, a departure of 4 and 9sigma from that expected if the discs were unresolved. HD69830 is unresolved by VISIR at 18.7um. The combined limits from MIDI and 8m imaging constrain the warm dust to lie within 0.05-2.4AU for HD69830 and 0.16-2.98AU for eta Corvi. These results represent the first resolution in the mid-IR of dust around main sequence stars. The constraints placed on the location of the dust are consistent with radii predicted by SED modelling. Tentative evidence for a common position angle for the dust at 1.7AU with that at 150AU around eta Corvi, which might be expected if the hot dust is fed from the outer disc, demonstrates the potential of this technique for constraining the origin of the dust and more generally for the study of dust in the terrestrial regions of main sequence stars.

The fine dust detected by IR emission around the nearby Beta Pic analogue star HD172555 is very peculiar. The dust mineralogy is composed primarily of highly refractory, non-equilibrium materials, with approximately three-quarters of the Si atoms in silica (SiO2) species. Tektite and obsidian lab thermal emission spectra (non-equilibrium glassy silicas found in impact and magmatic systems) are required to fit the data. The best-fit model size distribution for the observed fine dust is dn/da = a-3.95 +/- 0.10. This steep a size distribution, with abundant micron-sized particles, argues for a fresh source of material within the last 0.1 Myr. The location of the dust with respect to the star is at 5.8 +/- 0.6 AU (equivalent to 1.9 +/- 0.2 AU from the Sun), within the terrestrial planet formation region but at the outer edge of any possible terrestrial habitability zone. The mass of fine dust is 4 x 10^19 - 2 x 10^20 kg, equivalent to a 150 - 200 km radius asteroid. Significant emission features centered at 4 and 8 um due to fluorescing SiO gas are also found. Roughly 10^22 kg of SiO gas, formed by vaporizing silicate rock, is also present in the system, and a separate population of very large, cool grains, massing 10^21 - 10^22 kg and equivalent to the largest sized asteroid currently found in the Solar System's main asteroid belt, dominates the solid circumstellar material by mass. The makeup of the observed dust and gas, and the noted lack of a dense circumstellar gas disk, strong primary x-ray activity, or an extended disk of Beta-meteroids argues that the source of the observed circumstellar materials is a giant hypervelocity (> 10 km sec^-1) impact between large rocky planetesimals, similar to the ones which formed the Moon and which stripped the surface crustal material off of Mercury's surface.

HD 69830 exhibits radial velocity variations attributed to three planets as well as infrared emission attributed to a warm debris disk. Previous studies have developed models for the planet migration and mass growth (Alibert et al. 2006) and the replenishment of warm grains (Wyatt et al. 2007). We perform n-body integrations in order to explore the implications of these models for: 1) the excitation of planetary eccentricity, 2) the accretion and clearing of a putative planetesimal disk, 3) the distribution of planetesimal orbits following migration, and 4) the implications for the origin of the IR emission. We find that: i) It is not possible to explain the observed planetary eccentricities (e~0.1) purely as the result of planetary perturbations during migration unless the planetary system is nearly face-on. ii) The rate of accretion of planetesimals onto planets in our n-body simulations is significantly different to that assumed in the semi-analytic models, suggesting that one cannot successfully treat planetesimal accretion in the simplified manner of Alibert et al. (2006). iii) Eccentricity damping of planetesimals does not act as an insurmountable obstacle to the existence of an excited eccentric disk: All simulations result in at least 25 Earth-masses of material remaining bound in the region ~1-9 AU, even after all three planets have migrated through the region. iv) Gas drag works to size-sort the planetesimals, with the largest bodies preferentially occupying the highest eccentricity and longest-lived orbits. Further work will be required to understand whether these eccentricity distributions are high enough to explain the level of dust emission observed despite mass loss via steady state collisional evolution. [abridged]

60% of the A star members of the 12 Myr old beta Pictoris moving group (BPMG) show significant excess emission in the mid-infrared, several million years after the proto-planetary disk is thought to disperse. Theoretical models suggest this peak may coincide with the formation of Pluto-sized planetesimals in the disk, stirring smaller bodies into collisional destruction. Here we present resolved mid-infrared imaging of the disk of eta Tel (A0V in the BPMG) and consider its implications for the state of planet formation in this system. eta Tel was observed at 11.7 and 18.3um using T-ReCS on Gemini South. The resulting images were compared to simple disk models to constrain the radial distribution of the emitting material. The emission observed at 18.3um is shown to be significantly extended beyond the PSF along a position angle 8 degrees. This is the first time dust emission has been resolved around eta Tel. Modelling indicates that the extension arises from an edge-on disk of radius 0.5 arcsec (~24 AU). Combining the spatial constraints from the imaging with those from the spectral energy distribution shows that >50% of the 18um emission comes from an unresolved dust component at ~4 AU. The radial structure of the eta Tel debris disk is reminiscent of the Solar System, suggesting that this is a young Solar System analogue. For an age of 12Myr, both the radius and dust level of the extended cooler component are consistent with self-stirring models for a protoplanetary disk of 0.7 times minimum mass solar nebula. The origin of the hot dust component may arise in an asteroid belt undergoing collisional destruction, or in massive collisions in ongoing terrestrial planet formation.

We present the 3.5m SPICA space telescope, a proposed Japanese-led JAXA-ESA mission scheduled for launch around 2017. The actively cooled (<5 K), single aperture telescope and monolithic mirror will operate from ~3.5 to ~210 um and will provide superb sensitivity in the mid- and far-IR spectral domain (better than JWST at lambda > 18 um). SPICA is one of the few space missions selected to go to the next stage of ESA's Cosmic Vision 2015-2025 selection process. In this White Paper we present the main specifications of the three instruments currently baselined for SPICA: a mid-infrared (MIR) coronagraph (~3.5 to ~27 um) with photometric and spectral capabilities (R~200), a MIR wide-field camera and high resolution spectrometer (R~30,000), and a far-infrared (FIR ~30 to ~210 um) imaging spectrometer - SAFARI - led by a European consortium. We discuss their capabilities in the context of MIR direct observations of exo-planets (EPs) and multiband photometry/high resolution spectroscopy observations of transiting exo-planets. We conclude that SPICA will be able to characterize the atmospheres of transiting exo-planets down to the super-Earth size previously detected by ground- or space-based observatories. It will also directly detect and characterize Jupiter/Neptune-size planets orbiting at larger separation from their parent star (>5-10 AU), by performing quantitative atmospheric spectroscopy and studying proto-planetary and debris disks. In addition, SPICA will be a scientific and technological precursor for future, more ambitious, IR space missions for exo-planet direct detection as it will, for example, quantify the prevalence exo-zodiacal clouds in planetary systems and test coronographic techniques, cryogenic systems and lightweight, high quality telescopes. (abridged)

The number of stars that are known to have debris disks is greater than that of stars known to harbour planets. These disks are detected because dust is created in the destruction of planetesimals in the disks much in the same way that dust is produced in the asteroid belt and Kuiper belt in the solar system. For the nearest stars the structure of their debris disks can be directly imaged, showing a wide variety of both axisymmetric and asymmetric structures. A successful interpretation of these images requires a knowledge of the dynamics of small bodies in planetary systems, since this allows the observed dust distribution to be deconvolved to provide information on the distribution of larger objects, such as planetesimals and planets. This chapter reviews the structures seen in debris disks, and describes a disk dynamical theory which can be used to interpret those observations. In this way much of the observed structures, both axisymmetric and asymmetric, can be explained by a model in which the dust is produced in a planetesimal belt which is perturbed by a nearby, as yet unseen, planet. While the planet predictions still require confirmation, it is clear that debris disks have the potential to provide unique information about the structure of extrasolar planetary systems, since they can tell us about planets analogous to Neptune and even the Earth. Significant failings of the model at present are its inability to predict the quantity of small grains in a system, and to explain the origin of the transient dust seen in some systems. Given the complexity of planetary system dynamics and how that is readily reflected in the structure of a debris disk, it seems inevitable that the study of debris disks will play a vital role in our understanding of extrasolar planetary systems.

(ABRIDGED) Studies of debris disks have shown that most systems are analogous to the Edgeworth-Kuiper Belt. However a rare subset of sun-like stars possess dust which lies in the terrestrial planet region. In this study we aim to determine how many sources with apparent mid-IR excess are truly hosts of warm dust, and investigate where the dust must lie. We observed using mid-IR imaging with TIMMI2, VISIR and MICHELLE a sample of FGK main sequence stars reported to have hot dust. A new modelling approach was developed to determine the constraints that can be set on the radial extent of excess emission. We confirm the presence of warm dust around 3 of the candidates (eta Corvi, HD145263 and HD202406), and present constraints on the emitting dust regions. Of 2 alternative models for the eta Corvi excess emission, we find that a model with 1 hot dust component at <3 AU (combined with the known submm dust population) fits the data better at the 2.6sigma level than an alternative model with 2 populations of dust in the mid-IR. We identify several systems which have a companion (HD65277 and HD79873) or background object (HD53246, HD123356 and HD128400) responsible for their mid-infrared excess, and for 3 other systems we were able to rule out a point-like source near the star at the level of excess observed in lower resolution observations (HD12039, HD69830 and HD191089). Hot dust sources are either young and possibly primordial or transitional, or have relatively small radius steady-state planetesimal belts, or they are old and luminous with transient emission. High resolution imaging can be used to constrain the location of the disk and help to discriminate between different models of disk emission. For some small disks, interferometry is needed to resolve the disk location.

We present an analysis of the gas-poor circumstellar material in the HD 113766 binary system (F3/F5, 10 - 16 Myr), recently observed by the Spitzer Space Telescope. For our study we have used the infrared mineralogical model derived from observations of the Deep Impact experiment. We find the dust dominated by warm, fine (~1 um) particles, abundant in Mg-rich olivine, crystalline pyroxenes, amorphous silicates, Fe-rich sulfides, amorphous carbon, and colder water-ice. The warm dust material mix is akin to an inner main belt asteroid of S-type composition. The ~440 K effective temperature of the warm dust implies that the bulk of the observed material is in a narrow belt ~1.8 AU from the 4.4 L_solar central source, in the terrestrial planet-forming region and habitable zone of the system (equivalent to 0.9 AU in the solar system). The icy dust lies in 2 belts, located at 4-9 AU and at 30 - 80 AU. The lower bound of warm dust mass in 0.1 - 20 um, dn/da ~ a^-3.5 particles is very large, at least 3 x 10^20 kg, equivalent to a 320 km radius asteroid of 2.5 g cm^-3 density. Assuming 10m largest particles present, the lower bound of warm dust mass is at least 0.5 M_Mars The dust around HD 113766A originates from catastrophic disruption of terrestrial planet embryo(s) and subsequent grinding of the fragments, or from collisions in a young, extremely dense asteroid belt undergoing aggregation. The persistence of the strong IR excess over the last two decades argues for a mechanism to provide replenishment of the circumstellar material on yearly timescales.

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

There is evidence for the existence of massive planets at orbital radii of several hundred AU from their parent stars where the timescale for planet formation by core accretion is longer than the disc lifetime. These planets could have formed close to their star and then migrated outwards. We consider how the transfer of angular momentum by viscous disc interactions from a massive inner planet could cause significant outward migration of a smaller outer planet. We find that it is in principle possible for planets to migrate to large radii. We note, however, a number of effects which may render the process somewhat problematic.

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.

We present photometric detections of dust emission at 850 and 450 micron around the pre-main sequence M1 dwarf TWA 7 using the SCUBA camera on the James Clerk Maxwell Telescope. These data confirm the presence of a cold dust disk around TWA 7, a member of the TW Hydrae Association. Based on the 850 micron flux, we estimate the mass of the disk to be 18 lunar masses (0.2 Earth masses) assuming a mass opacity of 1.7 cm^2/g with a temperature of 45 K. This makes the TWA 7 disk (d=55 pc) an order of magnitude more massive than the disk reported around AU Microscopii (GL 803), the closest (9.9 pc) debris disk detected around an M dwarf. This is consistent with TWA 7 being slightly younger than AU Mic. We find that the mid-IR and submillimeter data require the disk to be comprised of dust at a range of temperatures. A model in which the dust is at a single radius from the star, with a range of temperatures according to grain size, is as effective at fitting the emission spectrum as a model in which the dust is of uniform size, but has a range of temperatures according to distance. We discuss this disk in the context of known disks in the TW Hydrae Association and around low-mass stars; a comparison of masses of disks in the TWA reveals no trend in mass or evolutionary state (gas-rich vs. debris) as a function of spectral type.

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.

The dust disks observed around mature stars are evidence that plantesimals are present in these systems on spatial scales that are similar to that of the asteroids and the KBOs in the Solar System. These dust disks (a.k.a. ``debris disks'') present a wide range of sizes, morphologies and properties. It is inferred that their dust mass declines with time as the dust-producing planetesimals get depleted, and that this decline can be punctuated by large spikes that are produced as a result of individual collisional events. The lack of solid state features indicate that, generally, the dust in these disks have sizes larger than approximately 10 microns, but exceptionally, strong silicate features in some disks suggest the presence of large quantities of small grains, thought to be the result of recent collisions. Spatially resolved observations of debris disks show a diversity of structural features, such as inner cavities, warps, offsets, brightness asymmetries, spirals, rings and clumps. There is growing evidence that, in some cases, these structures are the result of the dynamical perturbations of a massive planet. Our Solar System also harbors a debris disk and some of its properties resemble those of extra-solar debris disks. From the cratering record, we can infer that its dust mass has decayed with time, and that there was at least one major ``spike'' in the past during the Late Heavy Bombardment. This offers a unique opportunity to use extra-solar debris disks to shed some light in how the Solar System might have looked in the past. Similarly, our knowledge of the Solar System is influencing our understanding of the types of processes which might be at play in the extra-solar debris disks.

We have used the infrared mineralogical model derived from the Spitzer IRS observations of the Deep Impact experiment to study the nature of the dust in the debris found around the K0V star HD69830 (Beichman et al. 2005). Using a robust approach to determine the bulk average mineralogical composition of the dust, we show it to be substantially different from that found for comet 9P/Tempel 1 (Lisse et al. 2006) or C/Hale-Bopp 1995 O1 and comet-dominated YSO HD100546 (Lisse et al. 2007). Lacking in carbonaceous and ferrous materials but including small icy grains, the composition of the HD 69830 dust resembles that of a disrupted P or D-type asteroid. The amount of mass responsible for the observed emission is the equivalent of a 30 km radius, 2500 kg m-3 sphere, while the radiative temperature of the dust implies that the bulk of the observed material is at ~1.0 AU from the central source, coincident with the 2:1 and 5:2 mean motion resonances of the outermost of 3 Neptune-sized planets detected by Lovis et al. (2006). In our solar system, P and D-type asteroids are both large and numerous in the outer main belt and near Jupiter (e.g. the Hildas and Trojans) and have undergone major disruptive events to produce debris disk-like structures (c.f. the Karin and Veritas families 5-8 Myrs ago). The short-lived nature of the small and icy dust implies that the disruption occurred within the last year, or that replenishment due to ongoing collisional fragmentation is occurring.

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.

Debris disks are believed to be related to planetesimals left over around stars after planet formation has ceased. The frequency of debris disks around M-dwarfs which account for 70% of the stars in the Galaxy is unknown while constrains have already been found for A- to K-type stars. We have searched for cold debris disks around 32 field M-dwarfs by conducting observations at lambda = 850 microns with the SCUBA bolometerarray camera at the JCMT and at lambda = 1.2mm with the MAMBO array at the IRAM 30-m telescopes. This is the first survey of a large sample of M-dwarfs conducted to provide statistical constraints on debris disks around this type of stars. We have detected a new debris disk around the M0.5 dwarf GJ842.2 at lambda = 850 microns, providing evidence for cold dust at large distance from this star (~ 300AU). By combining the results of our survey with the ones of Liu et al. (2004), we estimate for the first time the detection rate of cold debris disks around field M-dwarfs with ages between 20 and 200 Myr. This detection rate is 13^{+6}_{-8} % and is consistent with the detection rate of cold debris disks (9 - 23 %) around A- to K-type main sequence stars of the same age. This is an indication that cold disks may be equally prevalent across stellar spectral types.

Over the past 10 years abundant evidence has emerged that many (if not all) stars are born with circumstellar disks. Understanding the evolution of post-accretion disks can provide strong constraints on theories of planet formation and evolution. In this review, we focus on developments in understanding: a) the evolution of the gas and dust content of circumstellar disks based on observational surveys, highlighting new results from the Spitzer Space Telescope; b) the physical properties of specific systems as a means to interpret the survey results; c) theoretical models used to explain the observations; d) an evolutionary model of our own solar system for comparison to the observations of debris disks around other stars; and e) how these new results impact our assessment of whether systems like our own are common or rare compared to the ensemble of normal stars in the disk of the Milky Way.

This paper considers the distribution of dust which originates in the break-up of planetesimals that are trapped in resonance with a planet. There are three grain populations with different spatial distributions: (I) large grains have the clumpy resonant distribution of the planetesimals; (II) moderate sized grains are no longer in resonance and have an axisymmetric distribution; (III) small grains are blown out of the system by radiation pressure and have a distribution which falls off oc 1/r, however these grains can be further divided into subclasses: (IIIa) grains produced from pop I that exhibit trailing spiral structure emanating from the clumps; and (IIIb) grains produced from pop II that have an axisymmetric distribution. Since observations in different wavebands are sensitive to different sized dust grains, multi-wavelength imaging can be used to test models for the origin of debris disk structure. For example, a disk with no blow-out grains would appear clumpy in the sub-mm, but smooth at mid- to far-IR wavelengths. The wavelength of the transition is indicative of the mass of the perturbing planet. The size distribution of Vega's disk is modeled in the light of the recent Spitzer observations. The origin of the large quantities of pop III grains seen by Spitzer must be in the destruction of the grains seen in the sub-mm, and so at high resolution and sensitivity the mid- to far-IR structure of Vega's disk is predicted to include spiral structure emanating from the sub-mm clumps.

This paper discusses the impact of introducing a planet on an eccentric orbit into a planetesimal disk. That planet's secular perturbations cause the orbits of the planetesimals to evolve in such a way that at any one time planetesimals at the same distance from the star have common pericentres and eccentricities. This causes the surface density distribution of an extended planetesimal disk to exhibit two spirals, one exterior the other interior to the planet's orbit. These two spirals unwind in different directions and their structure is described by two parameters: the time since the planet was introduced and the planet's eccentricity. At late times the spirals become tightly wound and the offset centre of symmetry of the pericentre glow approximation is recovered. Comparison with spiral structure seen in the HD141569 disk shows that its spiral at 325 AU is similar to that caused by introducing a planet 5 Myr ago with a mass 0.2-2M_Jup orbiting at 235-250 AU with an eccentricity of 0.05-0.2; likewise a Saturn mass planet at 150 AU would cause structure like that seen at 200 AU. More definitive statements about any planets orbiting HD141569 from this model could be made once the effect of the binary companion on the disk is known, and once the disk's structure has been better characterised down to 100 AU, including the location of the star within the disk. The relatively young age of this system (~5 Myr) means that if giant planets really do exist at hundreds of AU from HD141569, this provides a unique opportunity to set constraints on the mechanism by which those planets came to be at such large distances.

We present high spatial resolution mid- and far-infrared images of the Vega debris disk obtained with the Multiband Imaging Photometer for Spitzer (MIPS). The disk is well resolved and its angular size is much larger than found previously. The radius of the disk is at least 43" (330 AU), 70"(543 AU), and 105" (815 AU) in extent at 24, 70 and 160 um, respectively. The disk images are circular, smooth and without clumpiness at all three wavelengths. The radial surface brightness profiles imply an inner boundary at a radius of 11"+/-2" (86 AU). Assuming an amalgam of amorphous silicate and carbonaceous grains, the disk can be modeled as an axially symmetric and geometrically thin disk, viewed face-on, with the surface particle number density following an r^-1 power law. The disk radiometric properties are consistent with a range of models using grains of sizes ~1 to ~50 um. We find that a ring, containing grains larger than 180 um and at radii of 86-200 AU from the star, can reproduce the observed 850 um flux, while its emission does not violate the observed MIPS profiles. This ring could be associated with a population of larger asteroidal bodies analogous to our own Kuiper Belt. Cascades of collisions starting with encounters amongthese large bodies in the ring produce the small debris that is blown outward by radiation pressure to much larger distances where we detect its thermal emission. The dust production rate is >~10^15 g/s based on the MIPS results. This rate would require a very massive asteroidal reservoir for the dust to be produced in a steady state throughout Vega's life. Instead, we suggest that the disk we imaged is ephemeral and that we are witnessing the aftermath of a large and relatively recent collisional event, and subsequent collisional cascade.

This paper considers a simple model in which dust produced in a planetesimal belt migrates in toward the star due to P-R drag suffering destructive collisions with other dust grains on the way. Assuming the dust is all of the same size, the resulting surface density distribution can be derived analytically and depends only on the parameter eta0=5000tau*sqrt(Mstar/r)/beta; this parameter can be determined observationally with the hypothesis that beta=0.5. For massive belts in which eta0>>1 dust is confined to the planetesimal belt, while the surface density of more tenuous belts, in which eta0<<1, is constant with distance from the star. The emission spectrum of dust from planetesimal belts at different distances from different mass stars shows that the dust belts which have been detected to date should have eta0>>1; dust belts with eta0<<1 are hard to detect as they are much fainter than the stellar photosphere. This is confirmed for a sample of 37 debris disk candidates for which eta0 was determined to be >10. This means that these disks are so massive that mutual collisions prevent dust from reaching the inner regions of these systems and P-R drag can be ignored when studying their dynamics. Models for the formation of structure in debris disks by the trapping of particles into planetary resonances by P-R drag should be reconsidered. However, since collisions do not halt 100% of the dust, this means that in the absence of planetary companions debris disk systems should be populated by small quantities of hot dust which may be detectable in the mid-IR. Even in disks with eta0<<1 the temperature of dust emission is shown to be a reliable tracer of the planetesimal distribution meaning that inner holes in the dust distribution imply a lack of colliding planetesimals in the inner regions.

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.

We present new sub-mm observations and modelling of Vega excess stars, using realistic dust grain models. For resolved disks, we find that different objects require very different dust grain properties in order to simultaneously fit the image data and SED. Fomalhaut and Vega require solid dust grains, whilst HR4796 and HD141569 can only be fitted using porous grains. The older stars tend to have less porous grains than younger stars, which may indicate that collisions have compacted the dust grains. Eps Eri appears to be deficient in small dust grains compared to our best fitting model. This may be due to factors which affect the size distribution of grains close to the radiation pressure blowout limit. Alternatively, this discrepancy may be due to some external influence on the disk (e.g. a planet). When the model is applied to unresolved targets, an estimate of the disk size can be made. However, the large diversity in dust composition for the resolved disks means that we cannot make a reliable assumption as to the composition of the grains in an unresolved disk, and there is corresponding uncertainty in the disk size. In addition, the poor fit for Eps Eri shows that the model cannot always account for the SED even if the disk size is known. These two factors mean that it may not be possible to determine a disk's size without actually resolving it.

This paper describes a model which can explain the observed clumpy structures of debris disks. Clumps arise because after a planetary system forms its planets migrate due to angular momentum exchange with the remaining planetesimals. Outward migration of the outermost planet traps planetesimals outside its orbit into its resonances and resonant forces cause azimuthal structure in their distribution. The model is based on numerical simulations of planets of different masses, Mpl, migrating at different rates, dapl/dt, through a dynamically cold (e<0.01) planetesimal disk initially at a semimajor axis a. Trapping probabilities and the resulting azimuthal structures are presented for a planet's 2:1, 5:3, 3:2, and 4:3 resonances. Seven possible dynamical structures are identified from migrations defined by mu=Mpl/Mstar and theta=dapl/dt*sqrt(a/Mstar). Application of this model to the 850um image of Vega's disk shows its two clumps of unequal brightness can be explained by the migration of a Neptune-mass planet from 40 to 65AU over 56Myr; tight constraints are set on possible ranges of these parameters. The clumps are caused by planetesimals in the 3:2 and 2:1 resonances; the asymmetry arises because of the overabundance of planetesimals in the 2:1(u) over the 2:1(l) resonance. The similarity of this migration to that proposed for our own Neptune hints that Vega's planetary system may be much more akin to the solar system than previously thought. Predictions are made which would substantiate this model, such as the orbital motion of the clumpy pattern, the location of the planet, and the presence of lower level clumps.

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.

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.

This paper presents a model for the outcome of collisions between planetesimals in a debris disk and assesses the impact of collisional processes on the structure and size distribution of the disk. The model is presented by its application to Fomalhaut's collisionally replenished dust disk; a recent 450 micron image of this disk shows a clump embedded within it with a flux ~5 per cent of the total. The following conclusions are drawn: (i) SED modelling is consistent with Fomalhaut's disk having a collisional cascade size distribution extending from bodies 0.2 m in diameter down to 7 micron-sized dust. (ii) Collisional lifetime arguments imply that the cascade starts with planetesimals 1.5-4 km in diameter. Any larger bodies must be predominantly primordial. (iii) Constraints on the timescale for the ignition of the cascade are consistent with these primordial planetesimals having a distribution that extends up to 1000km, resulting in a disk mass of 5-10 times the minimum mass solar nebula. (iv) The debris disk is expected to be intrinsically clumpy, since planetesimal collisions result in dust clumps. The intrinsic clumpiness of Fomalhaut's disk is below current detection limits, but could be detectable by future observatories such as the ALMA, and could provide the only way of determining the primordial planetesimal population. (v) The observed clump could have originated in a collision between two runaway planetesimals, both larger than 1400 km diameter. It is unlikely that we should witness such an event unless both the formation of these runaways and the ignition of the collisional cascade occurred within the last ~10 Myr. (vi) Another explanation for Fomalhaut's clump is that ~5 per cent of the planetesimals in the ring are trapped in 1:2 resonance with a planet orbiting at 80 AU.

Recent advances in astronomical instrumentation mean that we are now able to image the thermal emission from the disks of dust around main sequence stars that may be the fossil remnants of planetary formation. These observations imply that the distribution of dust in the debris disks is neither smooth nor symmetrical; e.g., mid-infrared images of the disk of dust around the young A0V star HR 4796A show two lobes of emission, one of which may be ~5% brighter than the other. The observed structure of the debris disk in the solar system, i.e., the zodiacal cloud, also contains asymmetries: it has an offset center of symmetry, it is warped, and there is an asymmetric ring of dust co-orbiting with the Earth. Since the zodiacal cloud's asymmetries have been shown to be signatures of the gravitational perturbations of the solar system's planets, it is hoped that it may be possible to indirectly detect extrasolar planetary systems by their signatures in debris disk observations. This dissertation uses the physical processes that affect the evolution of debris material in the solar system to create a generalized model for the evolution of circumstellar debris material. It then shows how planetary perturbations affect that evolution, thereby causing the signatures of planets seen in the structure of the zodiacal cloud. This model can be used to provide a quantitative interpretation of debris disk observations, and the necessary modeling techniques are demonstrated by their application to observations of the HR 4796 disk. As well as determining the large scale structure of the HR 4796 disk, the modeling shows how a small body (>10 Earth mass) in the HR 4796 system that is on an orbit with an eccentricity larger than 0.02 could be the cause of the observed brightness asymmetry.

We present new 10.8 and 18.2 micron images of HR 4796A, a young A0V star that was recently discovered to have a spectacular, nearly edge-on, circumstellar disk prominent at ~20 microns (Jayawardhana et al. 1998; Koerner et al. 1998). These new images, obtained with OSCIR at Keck II, show that the disk's size at 10 microns is comparable to its size at 18 microns. Therefore, the 18 micron-emitting dust may also emit some, or all, of the 10 micron radiation. Using these multi-wavelength images, we determine a "characteristic" diameter of 2-3 microns for the mid-infrared-emitting dust particles if they are spherical and composed of astronomical silicates. Particles this small are expected to be blown out of the system by radiation pressure in a few hundred years, and therefore these particles are unlikely to be primordial. Dynamical modeling of the disk (Wyatt et al. 2000) indicates that the disk surface density is relatively sharply peaked near 70 AU, which agrees with the mean annular radius deduced by Schneider et al. (1999) from their NICMOS images. We present evidence (~1.8 sigma significance) for a brightness asymmetry that may result from the presence of the hole and the gravitational perturbation of the disk particle orbits by the low-mass stellar companion or a planet. This "pericenter glow," which must still be confirmed, results from a very small (a few AU) shift of the disk's center of symmetry relative to the central star HR 4796A; one side of the inner boundary of the annulus is shifted towards HR 4796A, thereby becoming warmer and more infrared-emitting. The possible detection of pericenter glow implies that the detection of even complex dynamical effects of planets on disks is within reach.

Recent images of the disks of dust around the young stars HR 4796A and Fomalhaut show, in each case, a double-lobed feature that may be asymmetric (one lobe may be brighter than the other). A symmetric double-lobed structure is that expected from a disk of dust with a central hole that is observed nearly edge-on (i.e., close to the plane of the disk). This paper shows how the gravitational influence of a second body in the system with an eccentric orbit would cause a brightness asymmetry in such a disk by imposing a "forced eccentricity" on the orbits of the constituent dust particles, thus shifting the center of symmetry of the disk away from the star and causing the dust near the forced pericenter of the perturbed disk to glow. Dynamic modeling of the HR 4796 disk shows that its 5% brightness asymmetry could be the result of a forced eccentricity as small as 0.02 imposed on the disk by either the binary companion HR 4796B, or by an unseen planet close to the inner edge of the disk. Since it is likely that a forced eccentricity of 0.01 or higher would be imposed on a disk in a system in which there are planets, but no binary companion, the corresponding asymmetry in the disk's structure could serve as a sensitive indicator of these planets that might otherwise remain undetected.