The new NIKA2 camera at the IRAM 30m radiotelescope was used to observe three known debris disks in order to constrain the SED of their dust emission in the millimeter wavelength domain. We have found that the spectral index between the two NIKA2 bands (1mm and 2mm) is consistent with the Rayleigh-Jeans regime (lambda^{-2}), unlike the steeper spectra (lambda^{-3}) measured in the submillimeter-wavelength domain for two of the three disks $-$ around the stars Vega and HD107146. We provide a succesful proof of concept to model this spectral inversion in using two populations of dust grains, those smaller and those larger than a grain radius a0 of 0.5mm. This is obtained in breaking the slope of the size distribution and the functional form of the absorption coefficient of the standard model at a0. The third disk - around the star HR8799 - does not exhibit this spectral inversion but is also the youngest.

Debris discs are second generation dusty discs formed by collisions of planetesimals. Many debris discs have been found and resolved around hot and solar-type stars. However, only a handful have been discovered around M-stars, and the reasons for their paucity remain unclear. Here we check whether the sensitivity and wavelength coverage of present-day telescopes are simply unfavourable for detection of these discs or if they are truly rare. We approach this question by looking at the Herschel/DEBRIS survey that has searched for debris discs including M-type stars. Assuming that these cool-star discs are "similar" to those of the hotter stars in some sense (i.e., in terms of dust location, temperature, fractional luminosity, or mass), we check whether this survey should have found them. With our procedure we can reproduce the $2.1^{+4.5}_{-1.7}$% detection rate of M-star debris discs of the DEBRIS survey, which implies that these discs can indeed be similar to discs around hotter stars and just avoid detection. We then apply this procedure to IRAM NIKA-2 and ALMA bands 3, 6 and 7 to predict possible detection rates and give recommendations for future observations. We do not favour observing with IRAM, since it leads to detection rates lower than for the DEBRIS survey, with 0.6%-4.5% for a 15 min observation. ALMA observations, with detection rates 0.9%-7.4%, do not offer a significant improvement either, and so we conclude that more sensitive far-infrared and single dish sub-millimetre telescopes are necessary to discover the missing population of M-star debris discs.

We present here new observations of the eccentric debris ring surrounding the Gyr-old solar-type star HD 202628: at millimeter wavelengths with ALMA, at far-infrared wavelengths with \textit{Herschel}, and in scattered light with \textit{HST}. The ring inner edge is found to be consistent between ALMA and \textit{HST} data. As radiation pressure affects small grains seen in scattered-light, the ring appears broader at optical than at millimeter wavelengths. The best fit to the ring seen with ALMA has inner and outer edges at $143.1 \pm 1.7$ AU and $165.5 \pm 1.4$, respectively, and an inclination of $57.4^\circ \pm 0.4$ from face-on. The offset of the ring centre of symmetry from the star allows us to quantify its eccentricity to be $e=0.09_{-0.01}^{+0.02}$. This eccentric feature is also detected in low resolution \textit{Herschel}/PACS observations, under the form of a pericenter-glow. Combining the infrared and millimeter photometry, we retrieve a disk grain size distribution index of $\sim -3.4$, and therefore exclude in-situ formation of the inferred belt-shaping perturber, for which we provide new dynamical constraints. Finally, ALMA images show four point-like sources that exceed 100$\,\mu$Jy, one of them being just interior to the ring. Although the presence of a background object cannot be excluded, we cannot exclude either that this source is circumplanetary material surrounding the belt-shaper, in which case degeneracies between its mass and orbital parameters could be lifted, allowing us to fully characterize such a distant planet in this mass and age regime for the very first time.

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 can be seen as the left-overs of giant planet formation and the possible nurseries of rocky planets. While M-type stars out-number more massive stars we know very little about the time evolution of their circumstellar disks at ages older than $\sim 10$\,Myr. Sub-millimeter observations are best to provide first order estimates of the available mass reservoir and thus better constrain the evolution of such disks. Here, we present ALMA Cycle\,3 Band\,7 observations of the debris disk around the M2 star TWA\,7, which had been postulated to harbor two spatially separated dust belts, based on unresolved far-infrared and sub-millimeter data. We show that most of the emission at wavelengths longer than $\sim 300$\,$\mu$m is in fact arising from a contaminant source, most likely a sub-mm galaxy, located at about 6.6" East of TWA\,7 (in 2016). Fortunately, the high resolution of our ALMA data allows us to disentangle the contaminant emission from that of the disc and report a significant detection of the disk in the sub-millimeter for the first time with a flux density of 2.1$\pm$0.4 mJy at 870 $\mu$m. With this detection, we show that the SED can be reproduced with a single dust belt.

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.

HR~8799 is a young F0-type star with four directly imaged giant planets and two debris belts, one located exterior and another one interior to the region occupied by the planetary orbits. Having an architecture similar to that of our Solar System, but also revealing dissimilarities such as high masses of planets, a huge radial extent and a high mass of the outer debris belt, HR~8799 is considered to be a benchmark to test formation and evolution models of planetary systems. Here we focus on the outer debris ring and its relation to the planets. We demonstrate that the models of the outer disc, proposed previously to reproduce Herschel observations, are inconsistent with the ALMA data, and vice versa. In an attempt to find a physically motivated model that would agree with both observational sets, we perform collisional simulations. We show that a narrow planetesimal belt and a radiation pressure induced dust halo cannot account for the observed radial brightness profiles. A single, wide planetesimal disc does not reproduce the data either. Instead, we propose a two-population model, comprising a Kuiper-Belt-like structure of a low-eccentricity planetesimal population ("the classical Kuiper Belt") and a high-eccentricity population of comets ("scattered disc"). We argue that such a structure of the exo-Kuiper belt of HR~8799 could be explained with planet migration scenarios analogous to those proposed for the Kuiper Belt of the Solar System.

One of the defining properties of debris discs compared to protoplanetary discs used to be their lack of gas, yet small amounts of gas have been found around an increasing number of debris discs in recent years. These debris discs found to have gas tend to be both young and bright. In this paper we conduct a deep search for CO gas in the system HD 95086 - a 17 Myr old, known planet host that also has a debris disc with a high fractional luminosity of $1.5\times10^{-3}$. Using the Atacama Large Millimeter/submillimeter Array (ALMA) we search for CO emission lines in bands 3, 6 and 7. By implementing a spectro-spatial filtering technique, we find tentative evidence for CO $J$=2-1 emission in the disc located at a velocity, 8.5$\pm$0.2 km s$^{-1}$, consistent with the radial velocity of the star. The tentative detection suggests that the gas on the East side of the disc is moving towards us. In the same region where continuum emission is detected, we find an integrated line flux of 9.5$\pm$3.6 mJy km s$^{-1}$, corresponding to a CO mass of (1.4-13)$\times10^{-6}$ M$_\oplus$. Our analysis confirms that the level of gas present in the disc is inconsistent with the presence of primordial gas in the system and is consistent with second generation production through the collisional cascade.

The presence of a debris disc around the Gyr-old solar-type star $\zeta^2\,\mathrm{Reticuli}$ was suggested by the $\mathit{Spitzer}$ infrared excess detection. Follow-up observations with $\mathit{Herschel}$/PACS revealed a double-lobed feature, that displayed asymmetries both in brightness and position. Therefore, the disc was thought to be edge-on and significantly eccentric. Here we present ALMA/ACA observations in Band 6 and 7 which unambiguously reveal that these lobes show no common proper motion with $\zeta^2\,\mathrm{Reticuli}$. In these observations, no flux has been detected around $\zeta^2\,\mathrm{Reticuli}$ that exceeds the $3\sigma$ levels. We conclude that surface brightness upper limits of a debris disc around $\zeta^2\,\mathrm{Reticuli}$ are $5.7\,\mathrm{\mu Jy/arcsec^2}$ at 1.3 mm, and $26\,\mathrm{\mu Jy/arcsec^2}$ at 870 microns. Our results overall demonstrate the capability of the ALMA/ACA to follow-up $\mathit{Herschel}$ observations of debris discs and clarify the effects of background confusion.

When a protoplanetary disc loses gas, it leaves behind planets and one or more planetesimal belts. The belts get dynamically excited, either by planets ('planet stirring') or by embedded big planetesimals ('self-stirring'). Collisions between planetesimals become destructive and start to produce dust, creating an observable debris disc. Following Kenyon & Bromley (2008), it is often assumed that self-stirring starts to operate as soon as the first ~1000 km-sized embedded 'Plutos' have formed. However, state-of-the-art pebble concentration models robustly predict planetesimals between a few km and ~200 km in size to form in protoplanetary discs rapidly, before then slowly growing into Pluto-sized bodies. We show that the timescale, on which these planetesimals excite the disc sufficiently for fragmentation, is shorter than the formation timescale of Plutos. Using an analytic model based on the Ida & Makino (1993) theory, backed up by numerical integrations, we find the excitation timescale to be $T_{\rm excite} \approx 100 x_{\rm m}^{-1} M_*^{-3/2} a^3$ Myr, where $x_{\rm m}$ is the total mass of a protoplanetary disc progenitor in the units of the Minimum-Mass Solar Nebula, $a$ its radius in the units of 100 AU, and $M_*$ is the stellar mass in solar masses. These results are applied to a set of 23 debris discs that have been well resolved with ALMA or SMA. We find that the majority of these discs are consistent with being self-stirred. However, three large discs around young early-type stars do require planets as stirrers. These are 49 Cet, HD 95086, and HR 8799, of which the latter two are already known to have planets.

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.

Planets and minor bodies such as asteroids, Kuiper-belt objects and comets are integral components of a planetary system. Interactions among them leave clues about the formation process of a planetary system. The signature of such interactions is most prominent through observations of its debris disk at millimeter wavelengths where emission is dominated by the population of large grains that stay close to their parent bodies. Here we present ALMA 1.3 mm observations of HD 95086, a young early-type star that hosts a directly imaged giant planet b and a massive debris disk with both asteroid- and Kuiper-belt analogs. The location of the Kuiper-belt analog is resolved for the first time. The system can be depicted as a broad ($\Delta R/R \sim$0.84), inclined (30\arcdeg$\pm$3\arcdeg) ring with millimeter emission peaked at 200$\pm$6 au from the star. The 1.3 mm disk emission is consistent with a broad disk with sharp boundaries from 106$\pm$6 to 320$\pm$20 au with a surface density distribution described by a power law with an index of --0.5$\pm$0.2. Our deep ALMA map also reveals a bright source located near the edge of the ring, whose brightness at 1.3 mm and potential spectral energy distribution are consistent with it being a luminous star-forming galaxy at high redshift. We set constraints on the orbital properties of planet b assuming co-planarity with the observed disk.

We present the first scattered light image of the debris disk around HD 129590, a ~1.3 M$_\odot$ G1V member of the Scorpius Centaurus association with age ~10-16 Myr. The debris disk is imaged with the high contrast imaging instrument SPHERE at the Very Large Telescope, and is revealed by both the IRDIS and IFS subsytems, operating in the H and YJ bands respectively. The disk has a high infrared luminosity of $L_{\textrm{IR}}/L_{\textrm{star}}$~5$\times$10$^{-3}$, and has been resolved in other studies using ALMA. We detect a nearly edge on ring, with evidence of an inner clearing. We fit the debris disk using a model characterized by a single bright ring, with radius ~60-70 AU, in broad agreement with previous analysis of the target SED. The disk is vertically thin, and has an inclination angle of ~75$^\circ$. Along with other previously imaged edge-on disks in the Sco-Cen association such as HD 110058, HD 115600, and HD 111520, this disk image will allow of the structure and morphology of very young debris disks, shortly after the epoch of planet formation has ceased.

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

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

We present the first 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.

High levels of dust have been detected in the immediate vicinity of many stars, both young and old. A promising scenario to explain the presence of this short-lived dust is that these analogues to the Zodiacal cloud (or exozodis) are refilled in situ through cometary activity and sublimation. As the reservoir of comets is not expected to be replenished, the presence of these exozodis in old systems has yet to be adequately explained. It was recently suggested that mean-motion resonances (MMR) with exterior planets on moderately eccentric ($\mathrm{e_p}\gtrsim 0.1$) orbits could scatter planetesimals on to cometary orbits with delays of the order of several 100 Myr. Theoretically, this mechanism is also expected to sustain continuous production of active comets once it has started, potentially over Gyr-timescales. We aim here to investigate the ability of this mechanism to generate scattering on to cometary orbits compatible with the production of an exozodi on long timescales. We combine analytical predictions and complementary numerical N-body simulations to study its characteristics. We show, using order of magnitude estimates, that via this mechanism, low mass discs comparable to the Kuiper Belt could sustain comet scattering at rates compatible with the presence of the exozodis which are detected around Solar-type stars, and on Gyr timescales. We also find that the levels of dust detected around Vega could be sustained via our proposed mechanism if an eccentric Jupiter-like planet were present exterior to the system's cold debris disc.

We present 1.3 mm observations of the Sun-like star $\tau$ Ceti with the Atacama Large Millimeter/submillimeter Array (ALMA) that probe angular scales of $\sim1$'' (4 AU). This first interferometric image of the $\tau$ Ceti system, which hosts both a debris disk and possible multiplanet system, shows emission from a nearly face-on belt of cold dust with a position angle of $90^\circ$ surrounding an unresolved central source at the stellar position. To characterize this emission structure, we fit parametric models to the millimeter visibilities. The resulting best-fit model yields an inner belt edge of $6.2^{+9.8}_{-4.6}$ AU, consistent with inferences from lower resolution, far-infrared Herschel observations. While the limited data at sufficiently short baselines preclude us from placing stronger constraints on the belt properties and its relation to the proposed five planet system, the observations do provide a strong lower limit on the fractional width of the belt, $\Delta R/R > 0.75$ with $99\%$ confidence. This fractional width is more similar to broad disks such as HD 107146 than narrow belts such as the Kuiper Belt and Fomalhaut. The unresolved central source has a higher flux density than the predicted flux of the stellar photosphere at 1.3 mm. Given previous measurements of an excess by a factor of $\sim2$ at 8.7 mm, this emission is likely due to a hot stellar chromosphere.

Many debris discs reveal a two-component structure, with an outer Kuiper-belt analogue and a warm inner component whose origin is still a matter of debate. One possibility is that warm emission stems from an "asteroid belt" closer in to the star. We consider a scenario in which a set of giant planets is formed in an initially extended planetesimal disc. These planets carve a broad gap around their orbits, splitting up the disc into the outer and the inner belts. After the gas dispersal, both belts undergo collisional evolution in a steady-state regime. This scenario is explored with detailed collisional simulations involving realistic physics to describe a long-term collisional depletion of the two-component disc. We find that the inner disc may be able to retain larger amounts of material at older ages than thought before on the basis of simplified analytic models. We show that the proposed scenario is consistent with a suite of thermal emission and scattered light observational data for a bright two-temperature debris disc around a nearby solar-type star q$^1$ Eridani. This implies a Solar System-like architecture of the system, with an outer massive "Kuiper belt", an inner "asteroid belt", and a few Neptune- to Jupiter-mass planets in between.

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

The star HR 8799 hosts one of the largest known debris discs and at least four giant planets. Previous observations have found evidence for a warm belt within the orbits of the planets, a cold planetesimal belt beyond their orbits and a halo of small grains. With the infrared data, it is hard to distinguish the planetesimal belt emission from that of the grains in the halo. With this in mind, the system has been observed with ALMA in band 6 (1.34 mm) using a compact array format. These observations allow the inner edge of the planetesimal belt to be resolved for the first time. A radial distribution of dust grains is fitted to the data using an MCMC method. The disc is best fit by a broad ring between $145^{+12}_{-12}$ AU and $429^{+37}_{-32}$ AU at an inclination of $40^{+5}_{-6}${\deg} and a position angle of $51^{+8}_{-8}${\deg}. A disc edge at ~145 AU is too far out to be explained simply by interactions with planet b, requiring either a more complicated dynamical history or an extra planet beyond the orbit of planet b.

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

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.

Cool debris discs are a relic of the planetesimal formation process around their host star, analogous to the solar system's Edgeworth-Kuiper belt. As such, they can be used as a proxy to probe the origin and formation of planetary systems like our own. The Herschel Open Time Key Programmes "DUst around NEarby Stars" (DUNES) and "Disc Emission via a Bias-free Reconnaissance in the Infrared/Submillimetre" (DEBRIS) observed many nearby, sun-like stars at far-infrared wavelengths seeking to detect and characterize the emission from their circumstellar dust. Excess emission attributable to the presence of dust was identified from around $\sim$ 20% of stars. Herschel's high angular resolution ($\sim$ 7" FWHM at 100 $\mu$m) provided the capacity for resolving debris belts around nearby stars with radial extents comparable to the solar system (50 to 100 au). As part of the DUNES and DEBRIS surveys, we obtained observations of three debris disc stars, HIP 22263 (HD 30495), HIP 62207 (HD 110897), and HIP 72848 (HD 131511), at far-infrared wavelengths with the Herschel PACS instrument. Combining these new images and photometry with ancilliary data from the literature, we undertook simultaneous multi-wavelength modelling of the discs' radial profiles and spectral energy distributions using three different methodologies: single annulus, modified black body, and a radiative transfer code. We present the first far-infrared spatially resolved images of these discs and new single-component debris disc models. We characterize the capacity of the models to reproduce the disc parameters based on marginally resolved emission through analysis of two sets of simulated systems (based on the HIP 22263 and HIP 62207 data) with the noise levels typical of the Herschel images. We find that the input parameter values are recovered well at noise levels attained in the observations presented here.

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

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

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

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.

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.

We present 1.3 millimeter ALMA Cycle 0 observations of the edge-on debris disk around the nearby, ~10 Myr-old, M-type star AU Mic. These observations obtain 0.6 arcsec (6 AU) resolution and reveal two distinct emission components: (1) the previously known dust belt that extends to a radius of 40 AU, and (2) a newly recognized central peak that remains unresolved. The cold dust belt of mass about 1 lunar mass is resolved in the radial direction with a rising emission profile that peaks sharply at the location of the outer edge of the "birth ring" of planetesimals hypothesized to explain the midplane scattered light gradients. No significant asymmetries are discerned in the structure or position of this dust belt. The central peak identified in the ALMA image is ~6 times brighter than the stellar photosphere, which indicates an additional emission process in the inner regions of the system. Emission from a stellar corona or activity may contribute, but the observations show no signs of temporal variations characteristic of radio-wave flares. We suggest that this central component may be dominated by dust emission from an inner planetesimal belt of mass about 0.01 lunar mass, consistent with a lack of emission shortward of 25 microns and a location <3 AU from the star. Future millimeter observations can test this assertion, as an inner dust belt should be readily separated from the central star at higher angular resolution.

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.

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.

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.

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

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.

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.

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]