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.Comment: 11 Pages Accepted to MNRA

Hamiltonian model of capture into mean motion resonance

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

Steady-state evolution of debris disks around A stars

In this paper a simple analytical model for the steady-state evolution of debris disks due to collisions is confronted with Spitzer observations of main sequence A stars. All stars are assumed to have planetesimal belts with a distribution of initial masses and radii. In the model disk mass is constant until the largest planetesimals reach collisional equilibrium whereupon the mass falls off oc 1/t. We find that the detection statistics and trends seen at both 24 and 70um can be fitted well by the model. While there is no need to invoke stochastic evolution or delayed stirring to explain the statistics, a moderate rate of stochastic events is not ruled out. Potentially anomalous systems are identified by a high dust luminosity compared with the maximum permissible in the model (HD3003, HD38678, HD115892, HD172555). Their planetesimals may have unusual properties (high strength or low eccentricity) or this dust could be transient. While transient phenomena are also favored for a few systems in the literature, the overall success of our model, which assumes planetesimals in all belts have the same strength, eccentricity and maximum size, suggests a large degree of uniformity in the outcome of planet formation. The distribution of planetesimal belt radii, once corrected for detection bias, follows N(r) oc r^{-0.8+-0.3} for 3-120AU. Since the inner edge is often attributed to an unseen planet, this provides a unique constraint on the planetary systems of A stars. It is also shown that P-R drag may sculpt the inner edges of A star disks close to the Spitzer detection threshold (HD2262, HD19356, HD106591, HD115892). This model can be readily applied to the interpretation of future surveys, and predictions are made for the upcoming SCUBA-2 survey, including that >17% of A stars should be detectable at 850um.Comment: Accepted by Ap

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.Comment: Accepted for publication in MNRAS; 24 pages, 19 figure

Predicting the frequencies of diverse exo-planetary systems

Extrasolar planetary systems range from hot Jupiters out to icy comet belts more distant than Pluto. We explain this diversity in a model where the mass of solids in the primordial circumstellar disk dictates the outcome. The star retains measures of the initial heavy-element (metal) abundance that can be used to map solid masses onto outcomes, and the frequencies of all classes are correctly predicted. The differing dependences on metallicity for forming massive planets and low-mass cometary bodies are also explained. By extrapolation, around two-thirds of stars have enough solids to form Earth-like planets, and a high rate is supported by the first detections of low-mass exo-planets.Comment: 5 pages, 2 figures; accepted by MNRA

ALMA and Herschel Observations of the Prototype Dusty and Polluted White Dwarf G29-38

ALMA Cycle 0 and Herschel PACS observations are reported for the prototype, nearest, and brightest example of a dusty and polluted white dwarf, G29-38. These long wavelength programs attempted to detect an outlying, parent population of bodies at 1-100 AU, from which originates the disrupted planetesimal debris that is observed within 0.01 AU and which exhibits L_IR/L = 0.039. No associated emission sources were detected in any of the data down to L_IR/L ~ 1e-4, generally ruling out cold dust masses greater than 1e24 - 1e25 g for reasonable grain sizes and properties in orbital regions corresponding to evolved versions of both asteroid and Kuiper belt analogs. Overall, these null detections are consistent with models of long-term collisional evolution in planetesimal disks, and the source regions for the disrupted parent bodies at stars like G29-38 may only be salient in exceptional circumstances, such as a recent instability. A larger sample of polluted white dwarfs, targeted with the full ALMA array, has the potential to unambiguously identify the parent source(s) of their planetary debris.Comment: 8 pages, 5 figures and 1 table. Accepted to MNRA

On the observability of resonant structures in planetesimal disks due to planetary migration

We present a thorough study of the impact of a migrating planet on a planetesimal disk, by exploring a broad range of masses and eccentricities for the planet. We discuss the sensitivity of the structures generated in debris disks to the basic planet parameters. We perform many N-body numerical simulations, using the symplectic integrator SWIFT, taking into account the gravitational influence of the star and the planet on massless test particles. A constant migration rate is assumed for the planet. The effect of planetary migration on the trapping of particles in mean motion resonances is found to be very sensitive to the initial eccentricity of the planet and of the planetesimals. A planetary eccentricity as low as 0.05 is enough to smear out all the resonant structures, except for the most massive planets. The planetesimals also initially have to be on orbits with a mean eccentricity of less than than 0.1 in order to keep the resonant clumps visible. This numerical work extends previous analytical studies and provides a collection of disk images that may help in interpreting the observations of structures in debris disks. Overall, it shows that stringent conditions must be fulfilled to obtain observable resonant structures in debris disks. Theoretical models of the origin of planetary migration will therefore have to explain how planetary systems remain in a suitable configuration to reproduce the observed structures.Comment: 16 pages, 13 figures. Accepted for publication in A&

Transience of hot dust around sun-like stars

There is currently debate over whether the dust content of planetary systems is stochastically regenerated or originates in planetesimal belts evolving in steady state. In this paper a simple model for the steady state evolution of debris disks due to collisions is developed and confronted with the properties of the emerging population of 7 sun-like stars that have hot dust <10AU. The model shows there is a maximum possible disk mass at a given age, since more massive primordial disks process their mass faster. The corresponding maximum dust luminosity is f_max=0.00016r^(7/3)/t_age. The majority (4/7) of the hot disks exceed this limit by >1000 and so cannot be the products of massive asteroid belts, rather the following systems must be undergoing transient events characterized by an unusually high dust content near the star: eta Corvi, HD69830, HD72905 and BD+20307. It is also shown that the hot dust cannot originate in a recent collision in an asteroid belt, since there is also a maximum rate at which collisions of sufficient magnitude to reproduce a given dust luminosity can occur. Further it is shown that the planetesimal belt feeding the dust in these systems must be located further from the star than the dust, typically at >2AU. Other notable properties of the 4 hot dust systems are: two also have a planetesimal belt at >10AU (eta Corvi and HD72905); one has 3 Neptune mass planets at <1AU (HD69830); all exhibit strong silicate features in the mid-IR. We consider the most likely origin for the dust in these systems to be a dynamical instability which scattered planetesimals inwards from a more distant planetesimal belt in an event akin to the Late Heavy Bombardment in our own system, the dust being released from such planetesimals in collisions and possibly also sublimation.Comment: 16 pages, accepted by ApJ, removed HD128400 as hot dust candidat

Research returns redux: a meta-analysis of the returns to agricultural R&D

A total of 294 studies of returns to agricultural R&D (including extension) were compiled and these studies provide 1,858 separate estimates of rates of return. This includes some extreme values, which are implausible. When the highest and lowest 2.5 percent of the rates of return were set aside, the estimated annual rates of return averaged 73 percent overall–88 percent for research only, 45 percent for research and extension, and 79 percent for extension only. But these averages reveal little meaningful information from a large and diverse body of literature, which provides rate-of-return estimates that are often not directly comparable. The purpose of this study was to go behind the averages, and try to account for the sources of differences, in a meta-analysis of the studies of returns to agricultural R&D. The results conform with the theory and prior beliefs in many ways. Several features of the methods used by research evaluators matter, in particular assumptions about lag lengths and the nature of the research-induced supply shift.Rate of return., Agricultural research.,