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

Escalating observations of exo-minor planets and their destroyed remnants both passing through the solar system and within white dwarf planetary systems motivate an understanding of the orbital history and fate of exo-Kuiper belts and planetesimal discs. Here we explore how the structure of a 40-1000 au annulus of planetesimals orbiting inside of a solar system analogue that is itself initially embedded within a stellar cluster environment varies as the star evolves through all of its stellar phases. We attempt this computationally challenging link in four parts: (1) by performing stellar cluster simulations lasting 100 Myr, (2) by making assumptions about the subsequent quiescent 11 Gyr main-sequence evolution, (3) by performing simulations throughout the giant branch phases of evolution, and (4) by making assumptions about the belt's evolution during the white dwarf phase. Throughout these stages, we estimate the planetesimals' gravitational responses to analogues of the four solar system giant planets, as well as to collisional grinding, Galactic tides, stellar flybys, and stellar radiation. We find that the imprint of stellar cluster dynamics on the architecture of $\gtrsim 100$ km-sized exo-Kuiper belt planetesimals is retained throughout all phases of stellar evolution unless violent gravitational instabilities are triggered either (1) amongst the giant planets, or (2) due to a close ($\ll 10^3$ au) stellar flyby. In the absence of these instabilities, these minor planets simply double their semimajor axis while retaining their primordial post-cluster eccentricity and inclination distributions, with implications for the free-floating planetesimal population and metal-polluted white dwarfs.

Resolved debris disk features (e.g., warps, offsets, edges and gaps, azimuthal asymmetries, radially thickened rings, scale heights) contain valuable information about the underlying planetary systems, such as the posited planet's mass, semi-major axis, and other orbital parameters. Most existing models assume a single planet is sculpting the disk feature, but recent observations of mature planetary systems (e.g., by radial velocity surveys or \textit{Kepler}) have revealed that many planets reside in multi-planet systems. Here we investigate if/how planet properties inferred from single-planet models are compromised when multiple planets reside in the system. For each disk feature, we build a two-planet model that includes a planet b with fixed parameters and a planet c with a full range of possible parameters. We investigate these two-planet systems and summarize the configurations for which assuming a single planet (i.e., planet b) leads to significantly flawed inferences of that planet's properties. We find that although disk features are usually primarily dominated by a single planet, when using single-planet models we are at risk of misinterpreting planet properties by orders of magnitude in extreme cases. Specifically, we are at high risk of misinterpreting planet properties from disk warps; at moderate risk from disk edges and gaps, radially thickened rings, and scale height features; and at low risk from host star-disk center offsets and azimuthal asymmetries. We summarize situations where we can infer the need to use a multi-planet model instead of a single-planet one from disk morphology dissimilarities.

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

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

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

Today, Pluto and Eris are the largest and most massive Trans-Neptunian Objects respectively. They are believed to be the last remnants of a population of planetesimals that has been reduced by >99% since the time of its formation. This reduction implies a primordial population of hundreds or thousands of Pluto-mass objects, and a mass-number distribution that could have extended to hundreds of Lunas, dozens of Mars, and several Earths. Such lost protoplanets would have left signatures in the dynamics of the present-day Trans-Neptunian Populations, and we statistically limit their primordial number by considering the survival of ultra-wide binary TNOs, the Cold Classical Kuiper belt, and the resonant population. We find that if the primordial mass-number distribution extended to masses greater than Pluto (~1e-3 Earth masses), it must have turned downwards to be no more top-heavy than roughly equal mass per log size, a significant deviation from the distribution observed between 1e-5 and 1e-3 Earth masses. We compare these limits to the predicted mass-number distribution of various planetesimal and proto-planet growth models. The limits derived here provide a test for future models of planetesimal formation.

Giant, wide-separation planets often lie in the gap between multiple, distinct rings of circumstellar debris: this is the case for the HR\,8799 and HD\,95086 systems, and even the solar system where the Asteroid and Kuiper belts enclose the four gas and ice giants. In the case that a debris disk, inferred from an infrared excess in the SED, is best modelled as two distinct temperatures, we infer the presence of two spatially separated rings of debris. Giant planets may well exist between these two belts of debris, and indeed could be responsible for the formation of the gap between these belts. We observe 24 such two-belt systems using the VLT/SPHERE high contrast imager, and interpret our results under the assumption that the gap is indeed formed by one or more giant planets. A theoretical minimum mass for each planet can then be calculated, based on the predicted dynamical timescales to clear debris. The typical dynamical lower limit is $\sim$0.2$M_J$ in this work, and in some cases exceeds 1$M_J$. Direct imaging data, meanwhile, is typically sensitive to planets down to $\sim$3.6$M_J$ at 1'', and 1.7$M_J$ in the best case. Together, these two limits tightly constrain the possible planetary systems present around each target, many of which will be detectable with the next generation of high-contrast imagers.

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

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

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

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

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

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

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.

Transiting circumbinary planets discovered by Kepler provide unique insight into binary star and planet formation. Several features of this new found population, for example the apparent pile-up of planets near the innermost stable orbit, may distinguish between formation theories. In this work, we determine how planet-planet scattering shapes planetary systems around binaries as compared to single stars. In particular, we look for signatures that arise due to differences in dynamical evolution in binary systems. We carry out a parameter study of N-body scattering simulations for four distinct planet populations around both binary and single stars. While binarity has little influence on the final system multiplicity or orbital distribution, the presence of a binary dramatically effects the means by which planets are lost from the system. Most circumbinary planets are lost due to ejections rather than planet-planet or planet-star collisions. The most massive planet in the system tends to control the evolution. Systems similar to the only observed multi-planet circumbinary system, Kepler-47, can arise from much more tightly packed, unstable systems. Only extreme initial conditions introduce differences in the final planet populations. Thus, we suggest that any intrinsic differences in the populations are imprinted by formation.

Large Kuiper Belt Objects are conventionally thought to have formed out of a massive planetesimal belt that is a few thousand times its current mass. Such a picture, however, is incompatible with multiple lines of evidence. Here, we present a new model for the conglomeration of Cold Classical Kuiper belt objects, out of a solid belt only a few times its current mass, or a few percent of the solid density in a Minimum Mass Solar Nebula. This is made possible by depositing most of the primordial mass in grains of size centimetre or smaller. These grains collide frequently and maintain a dynamically cold belt out of which large bodies grow efficiently: an order-unity fraction of the solid mass can be converted into large bodies, in contrast to the ~0.1% efficiency in conventional models. Such a light belt may represent the true outer edge of the Solar system, and it may have effectively halted the outward migration of Neptune. In addition to the high efficiency, our model can also produce a mass spectrum that peaks at an intermediate size, similar to the observed Cold Classicals, if one includes the effect of cratering collisions. In particular, the observed power-law break observed at ~30 km for Cold Classicals, one that has been interpreted as a result of collisional erosion, may be primordial in origin.

We study the efficiency of forming large bodies, starting from a sea of equal-sized planetesimals. This is likely one of the earlier steps of planet formation and relevant for the formation of the asteroid belt, the Kuiper belt and extra-solar debris disks. Here we consider the case that the seed planetesimals do not collide frequently enough for dynamical collisional to be important (the collisionless limit), using a newly constructed conglomeration code, and by carefully comparing numerical results with analytical scalings. In the absence of collisional cooling, as large bodies grow by accreting small bodies, the velocity dispersion of the small bodies ($u$) is increasingly excited. Growth passes from the well-known run-away stage (when $u$ is higher than the big bodies' hill velocity) to the newly discovered trans-hill stage (when $u$ and big bodies both grow, but $u$ remains at the big bodies' hill velocity). We find, concurring with the analytical understandings developed in Lithwick (2014), as well as previous numerical studies, that a size spectrum $dn/dR \propto R^{-4}$ results, and that the formation efficiency, defined as mass fraction in bodies much larger than the initial size, is $\sim {\mathrm{ a \,\, few}}\times R_\odot/a$, or $\sim 10^{-3}$ at the distance of the Kuiper belt. We argue that this extreme inefficiency invalidates the conventional conglomeration model for the formation of both our Kuiper belt and extra-solar debris disks. New theories, possibly involving direct gravitational collapse, or strong collisional cooling of small planetesimals, are required.

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

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

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

The origin of trace hydrogen in white dwarfs (WDs) with He-dominated atmospheres is a long-standing problem, one that cannot satisfactorily be explained by the historically-favoured hypothesis of accretion from the interstellar medium. Here we explore the possibility that the gradual accretion of exo-Oort cloud comets, which are a rich source of H, contributes to the apparent increase of trace H with WD cooling age. We determine how often remnant exo-Oort clouds, freshly excited from post-main-sequence stellar mass loss, dynamically inject comets inside the WD's Roche radius. We improve upon previous studies by considering a representative range of single WD masses (0.52-1.00 Solar masses) and incorporating different cloud architectures, giant branch stellar mass loss, stellar flybys, Galactic tides and a realistic escape ellipsoid in self-consistent numerical simulations that integrate beyond 8 Gyr ages of WD cooling. We find that about 10^{-5} of the material in an exo-Oort cloud is typically amassed onto the WD, and that the H deposits accumulate even as the cloud dissipates. This accumulation may account for the relatively large amount of trace H, 10^{22}-10^{25} g, that is determined frequently among WDs with cooling ages >= 1 Gyr. Our results also reaffirm the notion that exo-Oort cloud comets are not the primary agents of the metal budgets observed in polluted WD atmospheres.

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

The recent discovery of planets orbiting main sequence binaries will provide crucial constraints for theories of binary and planet formation. The formation pathway for these planets is complicated by uncertainties in the formation mechanism of the host stars. In this paper, we compare the dynamical states of single and binary star planetary systems. Specifically, we pose two questions: (1) What does it mean for a circumbinary system to be dynamically packed? (2) How many systems are required to differentiate between a population of packed or sparse planets? We determine when circumbinary systems become dynamically unstable as a function of the separation between the host-stars and the inner planet, and the first and second planets. We show that these represent unique stability constraints compared to single-star systems. We find that although the existing Kepler data is insufficient to distinguish between a population of packed or sparse circumbinary systems, a more thorough study of circumbinary TTVs combined with an order of magnitude increase in the number of systems may prove conclusive. Future space missions such as TESS provide the best opportunity for increasing the sample size.

Observations of dusty debris disks can be used to test theories of planetesimal coagulation. Planetesimals of sizes up to a couple thousand kms are embedded in these disks and their mutual collisions generate the small dust grains that are observed. The dust luminosities, when combined with information on the dust spatial extent and the system age, can be used to infer initial masses in the planetesimal belts. Carrying out such a procedure for a sample of debris disks around Sun-like stars, we reach the following two conclusions. First, if we assume that colliding planetesimals satisfy a primordial size spectrum of the form dn/ds ~ s^{-q}, observed disks strongly favor a value of q between 3.5 and 4, while both current theoretical expectations and statistics of Kuiper belt objects favor a somewhat larger value. Second, number densities of planetesimals are two to three orders of magnitude higher in detected disks than in the Kuiper belt, for comparably-sized objects. This is a surprise for the coagulation models. It would require a similar increase in the solid surface density of the primordial disk over that of the Minimum Mass Solar Nebula, which is unreasonable. Both of our conclusions are driven by the need to explain the presence of bright debris disks at a few Gyrs of age.