The inwards scattering of planetesimals towards white dwarfs is expected to be a stochastic process with variability on human time-scales. The planetesimals tidally disrupt at the Roche radius, producing dusty debris detectable as excess infrared emission. When sufficiently close to the white dwarf, this debris sublimates and accretes on to the white dwarf and pollutes its atmosphere. Studying this infrared emission around polluted white dwarfs can reveal how this planetary material arrives in their atmospheres. We report a near-infrared monitoring campaign of 34 white dwarfs with infrared excesses with the aim to search for variability in the dust emission. Time series photometry of these white dwarfs from the United Kingdom Infrared Telescope (Wide Field Camera) in the J, H and K bands were obtained over baselines of up to three years. We find no statistically significant variation in the dust emission in all three near-infrared bands. Specifically, we can rule out variability at ~1.3% for the 13 white dwarfs brighter than 16th mag in K band, and at ~10% for the 32 white dwarfs brighter than 18th mag over time-scales of three years. Although to date two white dwarfs, SDSS J095904.69-020047.6 and WD 1226+110, have shown K band variability, in our sample we see no evidence of new K band variability at these levels. One interpretation is that the tidal disruption events which lead to large variabilities are rare, occur on short time-scales, and after a few years the white dwarfs return to being stable in the near-infrared.

Metals observed in the atmospheres of white dwarfs suggest that many have recently accreted planetary bodies. In some cases, the compositions observed suggest the accretion of material dominantly from the core (or the mantle) of a differentiated planetary body. Collisions between differentiated exoplanetesimals produce such fragments. In this work, we take advantage of the large numbers of white dwarfs where at least one siderophile (core-loving) and one lithophile (rock-loving) species have been detected to assess how commonly exoplanetesimals differentiate. We utilise N-body simulations that track the fate of core and mantle material during the collisional evolution of planetary systems to show that most remnants of differentiated planetesimals retain core fractions similar to their parents, whilst some are extremely core-rich or mantle-rich. Comparison with the white dwarf data for calcium and iron indicates that the data are consistent with a model in which $66^{+4}_{-6}\%$ have accreted the remnants of differentiated planetesimals, whilst $31^{+5}_{-5}\%$ have Ca/Fe abundances altered by the effects of heating (although the former can be as high as $100\%$, if heating is ignored). These conclusions assume pollution by a single body and that collisional evolution retains similar features across diverse planetary systems. These results imply that both collisions and differentiation are key processes in exoplanetary systems. We highlight the need for a larger sample of polluted white dwarfs with precisely determined metal abundances to better understand the process of differentiation in exoplanetary systems.

Practically all known planet hosts will evolve into white dwarfs, and large parts of their planetary systems will survive this transition - the same is true for the solar system beyond the orbit of Mars. Spectroscopy of white dwarfs accreting planetary debris provides the most accurate insight into the bulk composition of exo-planets. Ground-based spectroscopic surveys of ~260, 000 white dwarfs detected with Gaia will identify >1000 evolved planetary systems, and high-throughput high-resolution space-based ultraviolet spectroscopy is essential to measure in detail their abundances. So far, evidence for two planetesimals orbiting closely around white dwarfs has been obtained, and their study provides important constraints on the composition and internal structure of these bodies. Major photometric and spectroscopic efforts will be necessary to assemble a sample of such close-in planetesimals that is sufficiently large to establish their properties as a population, and to deduce the architectures of the outer planetary systems from where they originated. Mid-infrared spectroscopy of the dusty disks will provide detailed mineralogical information of the debris, which, in combination with the elemental abundances measured from the white dwarf spectroscopy, will enable detailed physical modelling of the chemical, thermodynamic, and physical history of the accreted material. Flexible multi-epoch infrared observations are essential to determine the physical nature, and origin of the variability observed in many of the dusty disks. Finally, the direct detection of the outer reservoirs feeding material to the white dwarfs will require sensitive mid- and far-infrared capabilities.

We hypothesise that differences in the temperatures at which the rocky material condensed out of the nebula gas can lead to differences in the composition of key rocky species (e.g., Fe, Mg, Si, Ca, Al, Na) and thus planet bulk density. Such differences in the observed bulk density of planets may occur as a function of radial location and time of planet formation. In this work we show that the predicted differences are on the cusp of being detectable with current instrumentation. In fact, for HD 219134, the 10 % lower bulk density of planet b compared to planet c could be explained by enhancements in Ca, Al rich minerals. However, we also show that the 11 % uncertainties on the individual bulk densities are not sufficiently accurate to exclude the absence of a density difference as well as differences in volatile layers. Besides HD 219134 b, we demonstrate that 55 Cnc e and WASP-47 e are similar candidates of a new Super-Earth class that have no core and are rich in Ca and Al minerals which are among the first solids that condense from a cooling proto-planetary disc. Planets of this class have densities 10-20% lower than Earth-like compositions and may have very different interior dynamics, outgassing histories and magnetic fields compared to the majority of Super-Earths.

The infrared dust emission from the white dwarf GD 56 is found to rise and fall by 20% peak-to-peak over 11.2 yr, and is consistent with ongoing dust production and depletion. It is hypothesized that the dust is produced via collisions associated with an evolving dust disk, temporarily increasing the emitting surface of warm debris, and is subsequently destroyed or assimilated within a few years. The variations are consistent with debris that does not change temperature, indicating that dust is produced and depleted within a fixed range of orbital radii. Gas produced in collisions may rapidly re-condense onto grains, or may accrete onto the white dwarf surface on viscous timescales that are considerably longer than Poynting-Robertson drag for micron-sized dust. This potential delay in mass accretion rate change is consistent with multi-epoch spectra of the unchanging Ca II and Mg II absorption features in GD 56 over 15 yr, although the sampling is sparse. Overall these results indicate that collisions are likely to be the source of dust and gas, either inferred or observed, orbiting most or all polluted white dwarfs.

The most heavily polluted white dwarfs often show excess infrared radiation from circumstellar dust disks, which are modeled as a result of tidal disruption of extrasolar minor planets. Interaction of dust, gas, and disintegrating objects can all contribute to the dynamical evolution of these dust disks. Here, we report on two infrared variable dusty white dwarfs, SDSS J1228+1040 and G29-38. For SDSS J1228+1040, compared to the first measurements in 2007, the IRAC [3.6] and [4.5] fluxes decreased by 20% by 2014 to a level also seen in the recent 2018 observations. For G29-38, the infrared flux of the 10 $\mu$m silicate emission feature became 10% stronger between 2004 and 2007, We explore several scenarios that could account for these changes, including tidal disruption events, perturbation from a companion, and runaway accretion. No satisfactory causes are found for the flux drop in SDSS J1228+1040 due to the limited time coverage. Continuous tidal disruption of small planetesimals could increase the mass of small grains and concurrently change the strength of the 10 $\mu$m feature of G29-38. Dust disks around white dwarfs are actively evolving and we speculate that there could be different mechanisms responsible for the temporal changes of these disks.

White dwarfs that have accreted rocky planetary bodies provide unique insights regarding the bulk composition of exoplanetary material. The analysis presented here uses observed white dwarf atmospheric abundances to constrain both where in the planetary system the pollutant bodies originated, and the geological and collisional history of the pollutant bodies. At least one, but possibly up to nine, of the 17 systems analysed have accreted a body dominated by either core-like or mantle-like material. The approximately even spread in the core mass fraction of the pollutants and the lack of crust-rich pollutants in the 17 systems studied here suggest that the pollutants are often the fragments produced by the collision of larger differentiated bodies. The compositions of many pollutants exhibit trends related to elemental volatility, which we link to the temperatures and, thus, the locations at which these bodies formed. Our analysis found that the abundances observed in 11 of the 17 systems considered are consistent with the compositions of nearby stars in combination with a trend related to elemental volatility. The even spread and large range in the predicted formation location of the pollutants suggests that pollutants arrive in white dwarf atmospheres with a roughly equal efficiency from a wide range of radial locations. Ratios of elements with different condensation temperatures such as Ca/Mg, Na/Mg, and O/Mg distinguish between different formation temperatures, whilst pairs of ratios of siderophilic and lithophilic elements such as Fe/Mg, Ni/Mg and Al/Mg, Ca/Mg distinguish between temperature dependent trends and geological trends.

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.

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

At least 25% of white dwarfs show atmospheric pollution by metals, sometimes accompanied by detectable circumstellar dust/gas discs or (in the case of WD 1145+017) transiting disintegrating asteroids. Delivery of planetesimals to the white dwarf by orbiting planets is a leading candidate to explain these phenomena. Here, we study systems of planets and planetesimals undergoing planet-planet scattering triggered by the star's post-main sequence mass loss, and test whether this can maintain high rates of delivery over the several Gyr that they are observed. We find that low-mass planets (Earth to Neptune mass) are efficient deliverers of material and can maintain the delivery for Gyr. Unstable low-mass planetary systems reproduce the observed delayed onset of significant accretion, as well as the slow decay in accretion rates at late times. Higher-mass planets are less efficient, and the delivery only lasts a relatively brief time before the planetesimal populations are cleared. The orbital inclinations of bodies as they cross the white dwarf's Roche limit are roughly isotropic, implying that significant collisional interactions of asteroids, debris streams and discs can be expected. If planet-planet scattering is indeed responsible for the pollution of white dwarfs, many such objects, and their main-sequence progenitors, can be expected to host (currently undetectable) super-Earth planets on orbits of several au and beyond.

A warm/hot dust component (at temperature $>$ 300K) has been detected around $\sim$ 20% of stars. This component is called "exozodiacal dust" as it presents similarities with the zodiacal dust detected in our Solar System, even though its physical properties and spatial distribution can be significantly different. Understanding the origin and evolution of this dust is of crucial importance, not only because its presence could hamper future detections of Earth-like planets in their habitable zones, but also because it can provide invaluable information about the inner regions of planetary systems. In this review, we present a detailed overview of the observational techniques used in the detection and characterisation of exozodiacal dust clouds ("exozodis") and the results they have yielded so far, in particular regarding the incidence rate of exozodis as a function of crucial parameters such as stellar type and age, or the presence of an outer cold debris disc. We also present the important constraints that have been obtained, on dust size distribution and spatial location, by using state-of-the-art radiation transfer models on some of these systems. Finally, we investigate the crucial issue of how to explain the presence of exozodiacal dust around so many stars (regardless of their ages) despite the fact that such dust so close to its host star should disappear rapidly due to the coupled effect of collisions and stellar radiation pressure. Several potential mechanisms have been proposed to solve this paradox and are reviewed in detail in this paper. The review finishes by presenting the future of this growing field.

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.

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.

Context: Extended circumstellar emission has been detected within a few 100 milli-arcsec around > 10% of nearby main sequence stars using near-infrared interferometry. Follow-up observations using other techniques, should they yield similar results or non-detections, can provide strong constraints on the origin of the emission. They can also reveal the variability of the phenomenon. Aims: We aim to demonstrate the persistence of the phenomenon over time scales of a few years and to search for variability of our previously detected excesses. Methods: Using VLTI/PIONIER in H band we have carried out multi-epoch observations of the stars for which a near-infrared excess was previously detected with the same observing technique and instrument. The detection rates and distribution of the excesses from our original survey and the follow-up observations are compared statistically. A search for variability of the excesses in our time series is carried out based on the level of the broadband excesses. Results: In 12 of 16 follow-up observations, an excess is re-detected with a significance of > 2 sigma, and in 7 of 16 follow-up observations significant excess (> 3 sigma) has been re-detected. We statistically demonstrate with very high confidence that the phenomenon persists for the majority of the systems. We also present the first detection of potential variability in two sources. Conclusions: We conclude that the phenomenon responsible for the excesses persists over time scales of a few years for the majority of the systems. However, we also find that variability intrinsic to a target can cause it to have no significant excess at the time of a specific observation.

The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.

The gap between two component debris disks is often taken to be carved by intervening planets scattering away the remnant planetesimals. We employ N-body simulations to determine how the time needed to clear the gap depends on the location of the gap and the mass of the planets. We invert this relation, and provide an equation for the minimum planet mass, and another for the expected number of such planets, that must be present to produce an observed gap for a star of a given age. We show how this can be combined with upper limits on the planetary system from direct imaging non-detections (such as with GPI or SPHERE) to produce approximate knowledge of the planetary system.

Comets have been invoked in numerous studies as a potentially important source of dust and gas around stars, but none has studied the thermo-physical evolution, out-gassing rate, and dust ejection of these objects in such stellar systems. We investigate the thermo-physical evolution of comets in exo-planetary systems in order to provide valuable theoretical data required to interpret observations of gas and dust. We use a quasi 3D model of cometary nucleus to study the thermo-physical evolution of comets evolving around a single star from 0.1 to 50 AU, whose homogeneous luminosity varies from 0.1 to 70 solar luminosities. This paper provides mass ejection, lifetimes, and the rate of dust and water gas mass productions for comets as a function of the distance to the star and stellar luminosity. Results show significant physical changes to comets at high stellar luminosities. The models are presented in such a manner that they can be readily applied to any planetary system. By considering the examples of the Solar System, Vega and HD 69830, we show that dust grains released from sublimating comets have the potential to create the observed (exo)zodiacal emission. We show that observations can be reproduced by 1 to 2 massive comets or by a large number of comets whose orbits approach close to the star. Our conclusions depend on the stellar luminosity and the uncertain lifetime of the dust grains. We find, as in previous studies, that exozodiacal dust disks can only survive if replenished by a population of typically sized comets renewed from a large and cold reservoir of cometary bodies beyond the water ice line. These comets could reach the inner regions of the planetary system following scattering by a (giant) planet.

Metal pollution in white dwarf atmospheres is likely to be a signature of remnant planetary systems. Most explanations for this pollution predict a sharp decrease in the number of polluted systems with white dwarf cooling age. Observations do not confirm this trend, and metal pollution in old (1-5 Gyr) white dwarfs is difficult to explain. We propose an alternative, time-independent mechanism to produce the white dwarf pollution. The orbit of a wide binary companion can be perturbed by Galactic tides, approaching close to the primary star for the first time after billions of years of evolution on the white dwarf branch. We show that such a close approach perturbs a planetary system orbiting the white dwarf, scattering planetesimals onto star-grazing orbits, in a manner that could pollute the white dwarf's atmosphere. Our estimates find that this mechanism is likely to contribute to metal pollution, alongside other mechanisms, in up to a few percent of an observed sample of white dwarfs with wide binary companions, independent of white dwarf age. This age independence is the key difference between this wide binary mechanism and others mechanisms suggested in the literature to explain white dwarf pollution. Current observational samples are not large enough to assess whether this mechanism makes a significant contribution to the population of polluted white dwarfs, for which better constraints on the wide binary population are required, such as those that will be obtained in the near future with Gaia.

Exozodiacal dust is warm or hot dust found in the inner regions of planetary systems orbiting main sequence stars, in or around their habitable zones. The dust can be the most luminous component of extrasolar planetary systems, but predominantly emits in the near- to mid-infrared where it is outshone by the host star. Interferometry provides a unique method of separating this dusty emission from the stellar emission. The visitor instrument PIONIER at the Very Large Telescope Interferometer (VLTI) has been used to search for hot exozodiacal dust around a large sample of nearby main sequence stars. The results of this survey are summarised: 9 out of 85 stars show excess exozodiacal emission over the stellar photospheric emission.

Several lines of evidence indicate a non-chondritic composition for Bulk Earth. If Earth formed from the accretion of chondritic material, its non-chondritic composition, in particular the super-chondritic 142Nd/144Nd and low Mg/Fe ratios, might be explained by the collisional erosion of differentiated planetesimals during its formation. In this work we use an N-body code, that includes a state-of-the-art collision model, to follow the formation of protoplanets, similar to proto-Earth, from differentiated planetesimals (> 100 km) up to isolation mass (> 0.16 M_Earth). Collisions between differentiated bodies have the potential to change the core-mantle ratio of the accreted protoplanets. We show that sufficient mantle material can be stripped from the colliding bodies during runaway and oligarchic growth, such that the final protoplanets could have Mg/Fe and Si/Fe ratios similar to that of bulk Earth, but only if Earth is an extreme case and the core is assumed to contain 10% silicon by mass. This may indicate an important role for collisional differentiation during the giant impact phase if Earth formed from chondritic material.

Context. Detecting and characterizing circumstellar dust is a way to study the architecture and evolution of planetary systems. Cold dust in debris disks only traces the outer regions. Warm and hot exozodiacal dust needs to be studied in order to trace regions close to the habitable zone. Aims. We aim to determine the prevalence and to constrain the properties of hot exozodiacal dust around nearby main-sequence stars. Methods. We search a magnitude limited (H < 5) sample of 92 stars for bright exozodiacal dust using our VLTI visitor instrument PIONIER in the H-band. We derive statistics of the detection rate with respect to parameters such as the stellar spectral type and age or the presence of a debris disk in the outer regions of the systems. We derive more robust statistics by combining our sample with the results from our CHARA/FLUOR survey in the K-band. In addition, our spectrally dispersed data allows us to put constraints on the emission mechanism and the dust properties in the detected systems. Results. We find an over-all detection rate of bright exozodiacal dust in the H-band of 11% (9 out of 85 targets) and three tentative detections. The detection rate decreases from early type to late type stars and increases with the age of the host star. We do not confirm the tentative correlation between the presence of cold and hot dust found in our earlier analysis of the FLUOR sample alone. Our spectrally dispersed data suggest that either the dust is extremely hot or the emission is dominated by the scattered light in most cases. The implications of our results for the target selection of future terrestrial planet finding missions using direct imaging are discussed.

25%-50% of all white dwarfs (WDs) host observable and dynamically active remnant planetary systems based on the presence of close-in circumstellar dust and gas and photospheric metal pollution. Currently-accepted theoretical explanations for the origin of this matter include asteroids that survive the star's giant branch evolution at au-scale distances and are subsequently perturbed onto WD-grazing orbits following stellar mass loss. In this work we investigate the tidal disruption of these highly-eccentric (e > 0.98) asteroids as they approach and tidally disrupt around the WD. We analytically compute the disruption timescale and compare the result with fully self-consistent numerical simulations of rubble piles by using the N-body code PKDGRAV. We find that this timescale is highly dependent on the orbit's pericentre and largely independent of its semimajor axis. We establish that spherical asteroids readily break up and form highly eccentric collisionless rings, which do not accrete onto the WD without additional forces such as radiation or sublimation. This finding highlights the critical importance of such forces in the physics of WD planetary systems.

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.

Vega has been shown to host multiple dust populations, including both hot exo-zodiacal dust at sub-AU radii and a cold debris disk extending beyond 100 AU. We use dynamical simulations to show how Vega's hot dust can be created by long-range gravitational scattering of planetesimals from its cold outer regions. Planetesimals are scattered progressively inward by a system of 5-7 planets from 30-60 AU to very close-in. In successful simulations the outermost planets are typically Neptune-mass. The back-reaction of planetesimal scattering causes these planets to migrate outward and continually interact with fresh planetesimals, replenishing the source of scattered bodies. The most favorable cases for producing Vega's exo-zodi have negative radial mass gradients, with sub-Saturn- to Jupiter-mass inner planets at 5-10 AU and outer planets of 2.5 to 20 Earth masses. The mechanism fails if a Jupiter-sized planet exists beyond ~15 AU because the planet preferentially ejects planetesimals before they can reach the inner system. Direct-imaging planet searches can therefore directly test this mechanism.

High levels of exozodiacal dust have been observed in the inner regions of a large fraction of main sequence stars. Given the short lifetime of the observed small dust grains, these 'exozodis' are difficult to explain, especially for old (>100 Myr) stars. The exozodiacal dust may be observed as excess emission in the mid-infrared, or using interferometry. We hypothesise that exozodi are sustained by planetesimals scattered by planets inwards from an outer planetesimal belt, where collision timescales are long. In this work, we use N-body simulations to show that the outwards migration of a planet into a belt, driven by the scattering of planetesimals, can increase, or sustain, the rate at which planetesimals are scattered from the outer belt to the exozodi region. We hypothesise that this increase is sufficient to sustain the observed exozodi on Gyr timescales. No correlation between observations of an outer belt and an exozodi is required for this scenario to work, as the outer belt may be too faint to detect. If planetesimal driven migration does explain the observed exozodi, this work suggests that the presence of an exozodi indicates the presence of outer planets and a planetesimal belt.

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.

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

Imaging of debris disks has found evidence for both eccentric and offset disks. One hypothesis is that these provide evidence for massive perturbers that sculpt the observed structures. One such disk was recently observed in the far-IR by the Herschel Space Observatory around $\zeta^2$ Ret. In contrast with previously reported systems, the disk is significantly eccentric, and the system is Gyr-old. We aim to investigate the long-term evolution of eccentric structures in debris disks caused by a perturber on an eccentric orbit. Both analytical predictions and numerical N-body simulations are used to investigate the observable structures that could be produced by eccentric perturbers. The long-term evolution of the disk geometry is examined, with particular application to the $\zeta^2$ Ret system. In addition, synthetic images of the disk are produced for comparison with Herschel observations. We show that an eccentric companion can produce both the observed offsets and eccentric disks. Such effects are not immediate and we characterise the timescale required for the disk to develop to an eccentric state. For the case of $\zeta^2$ Ret, we place limits on the mass and orbit of the companion required to produce the observations. Synthetic images show that the pattern observed around $\zeta^2$ Ret can be produced by an eccentric disk seen close to edge-on, and allow us to bring additional constraints on the disk parameters of our model (disk flux, extent). We determine that eccentric planets or stellar companions can induce long-lived eccentric structures in debris disks. Observations of such eccentric structures provide potential evidence of the presence of such a companion in a planetary system. We consider the example of $\zeta^2$ Ret, whose observed eccentric disk can be explained by a distant companion at tens of AU, on an eccentric orbit ($e_p\gtrsim 0.3$).

The nearby star Fomalhaut harbours a cold, moderately eccentric dust belt with a sharp inner edge near 133 au. A low-mass, common proper motion companion (Fom b), was discovered near the inner edge and was identified as a planet candidate that could account for the belt morphology. However, the most recent orbit determination based on four epochs of astrometry over eight years reveals a highly eccentric orbit that appears to cross the belt in the sky plane projection. We perform here a full orbital determination based on the available astrometric data to independently validate the orbit estimates previously presented. Adopting our values for the orbital elements and their associated uncertainties, we then study the dynamical interaction between the planet and the dust ring, to check whether the proposed disk sculpting scenario by Fom b is plausible. We used a dedicated MCMC code to derive the statistical distributions of the orbital elements of Fom b. Then we used symplectic N-body integration to investigate the dynamics of the dust belt, as perturbed by a single planet. Different attempts were made assuming different masses for Fom b. We also performed a semi-analytical study to explain our results. Our results are in good agreement with others regarding the orbit of Fom b. We find that the orbit is highly eccentric, is close to apsidally aligned with the belt, and has a moderate mutual inclination relative to the belt plane of. If coplanar, this orbit crosses the disk. Our dynamical study then reveals that the observed planet could sculpt a transient belt configuration with a similar eccentricity to what is observed, but it would not be simultaneously apsidally aligned with the planet. This transient configuration only occurs a short time after the planet is placed on such an orbit (assuming an initially circular disk), a time that is inversely proportional to the planet's mass, and that is in any case much less than the 440 Myr age of the star. We constrain how long the observed dust belt could have survived with Fom b on its current orbit, as a function of its possible mass. This analysis leads us to conclude that Fom b is likely to have low mass, that it is unlikely to be responsible for the sculpting of the belt, and that it supports the hypothesis of a more massive, less eccentric planet companion Fom c.

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.

[Abridged] Debris disks are extrasolar analogs to the solar system planetesimal belts. The star Fomalhaut harbors a cold debris belt at 140 AU as well as evidence of a warm dust component, which is suspected of being a bright analog to the solar system's zodiacal dust. Interferometric observations obtained with the VLTI and the KIN have identified near- and mid-infrared excesses attributed to hot and warm exozodiacal dust in the inner few AU of the star. We performed parametric modeling of the exozodiacal disk using the GRaTeR radiative transfer code to reproduce the interferometric data, complemented by mid- to far-infrared measurements. A detailed treatment of sublimation temperatures was introduced to explore the hot population at the sublimation rim. We then used an analytical approach to successively testing several source mechanisms. A good fit to the data is found by two distinct dust populations: (1) very small, hence unbound, hot dust grains confined in a narrow region at the sublimation rim of carbonaceous material; (2) bound grains at 2 AU that are protected from sublimation and have a higher mass despite their fainter flux level. We propose that the hot dust is produced by the release of small carbon grains following the disruption of aggregates that originate from the warm component. A mechanism, such as gas braking, is required to further confine the small grains for a long enough time. In situ dust production could hardly be ensured for the age of the star, so the observed amount of dust must be triggered by intense dynamical activity. Fomalhaut may be representative of exozodis that are currently being surveyed worldwide. We propose a framework for reconciling the hot exozodi phenomenon with theoretical constraints: the hot component of Fomalhaut is likely the tip of the iceberg since it could originate from a warm counterpart residing near the ice line.

Excess emission, associated with warm, dust belts, commonly known as exozodis, has been observed around a third of nearby stars. The high levels of dust required to explain the observations are not generally consistent with steady-state evolution. A common suggestion is that the dust results from the aftermath of a dynamical instability, an event akin to the Solar System's Late Heavy Bombardment. In this work, we use a database of N-body simulations to investigate the aftermath of dynamical instabilities between giant planets in systems with outer planetesimal belts. We find that, whilst there is a significant increase in the mass of material scattered into the inner regions of the planetary system following an instability, this is a short-lived effect. Using the maximum lifetime of this material, we determine that even if every star has a planetary system that goes unstable, there is a very low probability that we observe more than a maximum of 1% of sun-like stars in the aftermath of an instability, and that the fraction of planetary systems currently in the aftermath of an instability is more likely to be limited to <0.06. This probability increases marginally for younger or higher mass stars. We conclude that the production of warm dust in the aftermath of dynamical instabilities is too short-lived to be the dominant source of the abundantly observed exozodiacal dust.

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.

Vega and Fomalhaut, are similar in terms of mass, ages, and global debris disk properties; therefore, they are often referred as "debris disk twins". We present Spitzer 10-35 um spectroscopic data centered at both stars, and identify warm, unresolved excess emission in the close vicinity of Vega for the first time. The properties of the warm excess in Vega are further characterized with ancillary photometry in the mid infrared and resolved images in the far-infrared and submillimeter wavelengths. The Vega warm excess shares many similar properties with the one found around Fomalhaut. The emission shortward of ~30 um from both warm components is well described as a blackbody emission of ~170 K. Interestingly, two other systems, eps Eri and HR 8799, also show such an unresolved warm dust using the same approach. These warm components may be analogous to the solar system's zodiacal dust cloud, but of far greater. The dust temperature and tentative detections in the submillimeter suggest the warm excess arises from dust associated with a planetesimal ring located near the water-frost line and presumably created by processes occurring at similar locations in other debris systems as well. We also review the properties of the 2 um hot excess around Vega and Fomalhaut, showing that the dust responsible for the hot excess is not spatially associated with the dust we detected in the warm belt. We suggest it may arise from hot nano grains trapped in the magnetic field of the star. Finally, the separation between the warm and cold belt is rather large with an orbital ratio >~10 in all four systems. In light of the current upper limits on the masses of planetary objects and the large gap, we discuss the possible implications for their underlying planetary architecture, and suggest that multiple, low-mass planets likely reside between the two belts in Vega and Fomalhaut.

High levels of exozodiacal dust are observed around a growing number of main sequence stars. The origin of such dust is not clear, given that it has a short lifetime against both collisions and radiative forces. Even a collisional cascade with km-sized parent bodies, as suggested to explain outer debris discs, cannot survive sufficiently long. In this work we investigate whether the observed exozodiacal dust could originate from an outer planetesimal belt. We investigate the scattering processes in stable planetary systems in order to determine whether sufficient material could be scattered inwards in order to retain the exozodiacal dust at its currently observed levels. We use N-body simulations to investigate the efficiency of this scattering and its dependence on the architecture of the planetary system. The results of these simulations can be used to assess the ability of hypothetical chains of planets to produce exozodi in observed systems. We find that for older (>100Myr) stars with exozodiacal dust, a massive, large radii (>20AU) outer belt and a chain of tightly packed, low-mass planets would be required in order to retain the dust at its currently observed levels. This brings into question how many, if any, real systems possess such a contrived architecture and are therefore capable of scattering at sufficiently high rates to retain exozodi dust on long timescales.

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.

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.

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.

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.

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.