Source author record

Mark C. Wyatt

Mark C. Wyatt appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

45works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

45 published item(s)

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) III: The vertical structure of debris disks

Debris disks -- collisionally sustained belts of dust and sometimes gas around main sequence stars -- are remnants of planet formation processes and are found in systems ${\gtrsim}10$ Myr old. Millimeter-wavelength observations are particularly important, as the grains probed by these observations are not strongly affected by radiation pressure and stellar winds, allowing them to probe the dynamics of large bodies producing dust. The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) is analyzing high-resolution observations of 24 debris disks to enable the characterization of debris disk substructures across a large sample for the first time. For the most highly inclined disks, it is possible to recover the vertical structure of the disk. We aim to model and analyze the most highly inclined systems in the ARKS sample in order to uniformly extract the vertical dust distributions for a sample of well-resolved debris disks. We employed both parametric and nonparametric methods to constrain the vertical dust distributions for the most highly inclined ARKS targets. We find a broad range of aspect ratios, revealing a wide diversity in vertical structure, with a range of best-fit parametric values of $0.0026 \leq h_{\rm HWHM} \leq 0.193$ and a median best-fit value of $h_{\rm HWHM}=0.021$. The results obtained by nonparametric modeling are generally consistent with the parametric modeling results. We find that five of the 13 disks are consistent with having total disk masses less than that of Neptune (17 $M_{\oplus}$), assuming stirring by internal processes (self-stirring and collisional and frictional damping). Furthermore, most systems show a significant preference for a Lorentzian vertical profile rather than a Gaussian.

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) IX: Gas-driven origin for the continuum arc in the debris disc of HD 121617

Debris discs were long considered to be largely gas-free environments governed by collisional fragmentation, gravitational stirring, and radiative forces. Recent CO detections show that gas is present, but its abundance and origin remain uncertain. The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) revealed a narrow gas and dust ring in the disc HD 121617 with an asymmetric arc 40% brighter than the rest of the ring. We aim to constrain the total gas mass in HD 121617 assuming the dust arc is produced by hydrodynamical gas-dust interactions. We used the Dusty FARGO-ADSG code, modelling dust as Lagrangian particles, including radiation pressure and dust feedback, and varying the total gas mass. Simulations were compared to observations using radiative transfer. An unstable gas ring creates a size-dependent radial and azimuthal dust trap whose efficiency depends on gas mass. Two models, with 50 and 5 Earth masses of gas, reproduce both the ALMA band 7 arc and the outward offset of the VLT/SPHERE scattered-light ring via gas drag and radiation pressure. We infer a conservative gas-mass range of 2.5 to 250 Earth masses. If the ALMA asymmetry is caused by gas drag, the required gas mass compared with the observed CO implies substantial H2, consistent with primordial gas. HD 121617 would then be a hybrid disc between protoplanetary and debris stages. Since a planet could also create an arc, future observations are needed to distinguish these scenarios.

preprint2023arXiv

Has the dust clump in the debris disk of Beta Pictoris moved?

The edge-on debris disk of the nearby young star Beta Pictoris shows an unusual brightness asymmetry in the form of a clump. The clump has been detected in both the mid-IR and CO and its origin has so far remained uncertain. Here we present new mid-IR observations of Beta Pic to track any motion of the dust clump. Together with previous observations, the data span a period of 12 years. We measured any projected displacement of the dust clump over the 12-yr period to be $0.2^{+1.3}_{-1.4}$ au away from the star based on the median and 1$σ$ uncertainty, and constrain this displacement to be <11 au at the 3$σ$ level. This implies that the observed motion is incompatible with Keplerian motion at the 2.8$σ$ level. It has been posited that a planet migrating outwards may trap planetesimals into a 2:1 resonance, resulting in the observed clump at pericentre of their orbits that trails the planet. The observed motion is also incompatible with such resonant motion at the 2.6$σ$ level. While Keplerian motion and resonant motion is still possible, the data suggest that the dust clump is more likely stationary. Such a stationary dust clump could originate from the collision or tidal disruption of a planet-sized body, or from secular perturbations due to a planet that create regions with enhanced densities in the disk.

preprint2022arXiv

ALMA's view of the M-dwarf GSC 07396-00759's edge-on debris disc: AU Mic's coeval twin

We present new ALMA Band 7 observations of the edge-on debris disc around the M1V star GSC 07396-00759. At ~20 Myr old and in the beta Pictoris Moving Group along with AU Mic, GSC 07396-00759 joins it in the handful of low mass M-dwarf discs to be resolved in the sub-mm. With previous VLT/SPHERE scattered light observations we present a multi-wavelength view of the dust distribution within the system under the effects of stellar wind forces. We find the mm dust grains to be well described by a Gaussian torus at 70 au with a FWHM of 48 au and we do not detect the presence of CO in the system. Our ALMA model radius is significantly smaller than the radius derived from polarimetric scattered light observations, implying complex behaviour in the scattering phase function. The brightness asymmetry in the disc observed in scattered light is not recovered in the ALMA observations, implying that the physical mechanism only affects smaller grain sizes. High resolution follow-up observations of the system would allow investigation into its unique dust features as well as provide a true coeval comparison for its smaller sibling AU Mic, singularly well observed amongst M-dwarfs systems.

preprint2022arXiv

An ALMA 1.3 millimeter Search for Debris Disks around Solar-type Stars in the Pleiades

Millimeter emission from debris disks around stars of different ages provides constraints on the collisional evolution of planetesimals. We present ALMA 1.3 millimeter observations of a sample of 76 Solar-type stars in the ~115 Myr old Pleiades star cluster. These ALMA observations complement previous infrared observations of this sample by providing sensitivity to emission from circumstellar dust at lower temperatures, corresponding to debris at radii comparable to the Kuiper Belt and beyond. The observations obtain a beam size of 1.5 arcsec (200 au) and a median rms noise of 54 mircoJy/beam, which corresponds to a fractional luminosity $L_{dust}/L_{star} \sim 10^{-4}$ for 40 K dust for a typical star in the sample. The ALMA images show no significant detections of the targeted stars. We interpret these limits in the context of a steady-state collisional cascade model for debris disk evolution that provides a good description of observations of the field population near the Sun but is not well-calibrated on younger populations.The ALMA non-detections of the Pleiades systems are compatible with the disk flux predictions of this model. We find no high fractional luminosity outliers from these ALMA data that could be associated with enhanced collisions resulting from activity not accounted for by steady-state evolution. However, we note that two systems (HII 1132 and HD 22680) show 24 micron excess much higher than the predictions of this model, perhaps due to unusually high dust production from dynamical events involving planets.

preprint2022arXiv

Polarization from Aligned Dust Grains in the $β$ Pic Debris Disk

We present 870 $μ$m ALMA polarization observations of thermal dust emission from the iconic, edge-on debris disk $β$ Pic. While the spatially resolved map does not exhibit detectable polarized dust emission, we detect polarization at the $\sim$3$σ$ level when averaging the emission across the entire disk. The corresponding polarization fraction is $P_\textrm{frac}$ = $0.51 \pm 0.19$%. The polarization position angle $χ$ is aligned with the minor axis of the disk, as expected from models of dust grains aligned via radiative alignment torques (RAT) with respect to a toroidal magnetic field ($B$-RAT) or with respect to the anisotropy in the radiation field ($k$-RAT). When averaging the polarized emission across the outer versus inner thirds of the disk, we find that the polarization arises primarily from the SW third. We perform synthetic observations assuming grain alignment via both $k$-RAT and $B$-RAT. Both models produce polarization fractions close to our observed value when the emission is averaged across the entire disk. When we average the models in the inner versus outer thirds of the disk, we find that $k$-RAT is the likely mechanism producing the polarized emission in $β$ Pic. A comparison of timescales relevant to grain alignment also yields the same conclusion. For dust grains with realistic aspect ratios (i.e., $s > 1.1$), our models imply low grain-alignment efficiencies.

preprint2022arXiv

Rave: A non-parametric method for recovering the surface brightness and height profiles of edge-on debris disks

Extrasolar analogues of the Solar System's Kuiper belt offer unique constraints on outer planetary system architecture. Radial features such as the sharpness of disk edges and substructures such as gaps may be indicative of embedded planets within a disk. Vertically, the height of a disk can constrain the mass of embedded bodies. Edge-on debris disks offer a unique opportunity to simultaneously access the radial and vertical distribution of material, however recovering either distribution in an unbiased way is challenging. In this study, we present a non-parametric method to recover the surface brightness profile (face-on surface brightness as a function of radius) and height profile (scale height as a function of radius) of azimuthally symmetric, edge-on debris disks. The method is primarily designed for observations at thermal emission wavelengths, but is also applicable to scattered light observations under the assumption of isotropic scattering. By removing assumptions on underlying functional forms, this algorithm provides more realistic constraints on disk structures. We also apply this technique to ALMA observations of the AU Mic debris disk and derive a surface brightness profile consistent with estimates from parametric approaches, but with a more realistic range of possible models that is independent of parametrisation assumptions. Our results are consistent with a uniform scale height of 0.8 au, but a scale height that increases linearly with radius is also possible.

preprint2020arXiv

Dust Populations in the Iconic Vega Planetary System Resolved by ALMA

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

preprint2020arXiv

Dust size and spatial distributions in debris discs: predictions for exozodiacal dust dragged in from an exo-Kuiper belt

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

preprint2020arXiv

Evidence for a high mutual inclination between the cold Jupiter and transiting super Earth orbiting $π$ Men

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

preprint2020arXiv

Solution to the debris disc mass problem: planetesimals are born small?

Debris belts on the periphery of planetary systems, encompassing the region occupied by planetary orbits, are massive analogues of the Solar system's Kuiper belt. They are detected by thermal emission of dust released in collisions amongst directly unobservable larger bodies that carry most of the debris disc mass. We estimate the total mass of the discs by extrapolating up the mass of emitting dust with the help of collisional cascade models. The resulting mass of bright debris discs appears to be unrealistically large, exceeding the mass of solids available in the systems at the preceding protoplanetary stage. We discuss this "mass problem" in detail and investigate possible solutions to it. These include uncertainties in the dust opacity and planetesimal strength, variation of the bulk density with size, steepening of the size distribution by damping processes, the role of the unknown "collisional age" of the discs, and dust production in recent giant impacts. While we cannot rule out the possibility that a combination of these might help, we argue that the easiest solution would be to assume that planetesimals in systems with bright debris discs were "born small", with sizes in the kilometre range, especially at large distances from the stars. This conclusion would necessitate revisions to the existing planetesimal formation models, and may have a range of implications for planet formation. We also discuss potential tests to constrain the largest planetesimal sizes and debris disc masses.

preprint2020arXiv

The HOSTS survey for exozodiacal dust: Observational results from the complete survey

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.

preprint2019arXiv

Extrasolar Kuiper Belts

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

preprint2016arXiv

Discovery of a low-mass companion inside the debris ring surrounding the F5V star HD206893

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

preprint2016arXiv

Insights into planet formation from debris disks: II. Giant impacts in extrasolar planetary systems

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

preprint2015arXiv

Constraining the orbits of sub-stellar companions imaged over short orbital arcs

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

preprint2015arXiv

Dynamical Constraints on Outer Planets in Super-Earth Systems

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

preprint2015arXiv

Kuiper belt structure around nearby super-Earth host stars

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

preprint2015arXiv

Target Selection for the LBTI Exozodi Key Science Program

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

preprint2015arXiv

The Inner Disk Structure, Disk-Planet Interactions, and Temporal Evolution in the Beta Pictoris System: A Two-Epoch HST/STIS Coronagraphic Study

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

preprint2014arXiv

Debris from giant impacts between planetary embryos at large orbital radii

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.

preprint2014arXiv

Do Two Temperature Debris Disks Have Multiple Belts?

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

preprint2014arXiv

Dynamical evolution of an eccentric planet and a less massive debris disc

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

preprint2014arXiv

Exo-zodi modelling for the Large Binocular Telescope Interferometer

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.

preprint2014arXiv

Five steps in the evolution from protoplanetary to debris disk

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

preprint2014arXiv

Observations, Modeling and Theory of Debris Disks

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

preprint2013arXiv

Herschel Observations of Debris Discs Orbiting Planet-hosting Subgiants

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.

preprint2013arXiv

Imaged sub-stellar companions: not as eccentric as they appear? The effect of an unseen inner mass on derived orbits

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

preprint2013arXiv

Simulations of two-planet systems through all phases of stellar evolution: implications for the instability boundary and white dwarf pollution

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.

preprint2013arXiv

Spatially Resolved Images of Dust Belt(s) Around the Planet-hosting Subgiant Kappa CrB

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.

preprint2013arXiv

The bright end of the exo-Zodi luminosity function: Disk evolution and implications for exo-Earth detectability

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

preprint2013arXiv

The Great Escape III: Placing post-main-sequence evolution of planetary and binary systems in a Galactic context

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

preprint2012arXiv

Debris disks as signposts of terrestrial planet formation. II Dependence of exoplanet architectures on giant planet and disk properties

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

preprint2012arXiv

Debris from terrestrial planet formation: the Moon-forming collision

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

preprint2012arXiv

Resolved Debris Discs Around A Stars in the Herschel DEBRIS Survey

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

preprint2012arXiv

The Solar System's Post-Main Sequence Escape Boundary

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

preprint2011arXiv

Debris disk size distributions: steady state collisional evolution with P-R drag and other loss processes

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

preprint2011arXiv

Debris disks as signposts of terrestrial planet formation

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

preprint2011arXiv

Dependence of a planet's chaotic zone on particle eccentricity: the shape of debris disc inner edges

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

preprint2011arXiv

Searching for Saturn's Dust Swarm: Limits on the size distribution of Irregular Satellites from km to micron sizes

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

preprint2011arXiv

The Great Escape: How Exoplanets and Smaller Bodies Desert Dying Stars

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.

preprint2010arXiv

A general model of resonance capture in planetary systems: First and second order resonances

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

preprint2010arXiv

Collisional Evolution of Irregular Satellite Swarms: Detectable Dust around Solar System and Extrasolar Planets

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

preprint2010arXiv

Resolved Imaging of the HD191089 Debris Disc

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

preprint2009arXiv

Outward Migration of Terrestrial Embryos in Binary Systems

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