Metallicity, debris discs and planets

We investigate the populations of main-sequence stars within 25 pc that have debris discs and/or giant planets detected by Doppler shift. The metallicity distribution of the debris sample is a very close match to that of stars in general, but differs with >99 per cent confidence from the giant planet sample, which favours stars of above average metallicity. This result is not due to differences in age of the two samples. The formation of debris-generating planetesimals at tens of au thus appears independent of the metal fraction of the primordial disc, in contrast to the growth and migration history of giant planets within a few au. The data generally fit a core accumulation model, with outer planetesimals forming eventually even from a disc low in solids, while inner planets require fast core growth for gas to still be present to make an atmosphere.     

Migration of giant planets in planetesimal discs

Planets orbiting a planetesimal circumstellar disc can migrate inward from their initial positions because of dynamical friction between planets and planetesimals. The migration rate depends on the disc mass and on its time evolution. Planets that are embedded in long-lived planetesimal discs, having total mass of 10⁻⁴– 0.01 M_⊙ , can migrate inward a large distance and can survive only if the inner disc is truncated or as a result of tidal interaction with the star. In this case the semimajor axis, a , of the planetary orbit is less than 0.1 au. Orbits with larger a are obtained for smaller values of the disc mass or for a rapid evolution (depletion) of the disc. This model may explain not only several of the orbital features of the giant planets that have been discovered in recent years orbiting nearby stars, but also the metallicity enhancement found in several stars associated with short-period planets.     

A search for protoplanetary discs around millisecond pulsars

Since the discovery of companions to B1257+12, it has been known that planets can exist around pulsars. Such planets may be formed in discs analogous to those around young stars, so we have searched for dust grain emission towards a sample of nine nearby millisecond pulsars. No emission is detected down to typical 2 σ limits of 5 mJy, at a wavelength of 850 μm. Using a model in which grains are heated by the pulsar spin-down luminosity, these dust flux limits correspond to disc masses of typically 10 Earth masses. The low dust limits show that nearby pulsar planets must already exist, rather than be in the process of forming, but only B1257+12 is known to have such planets. Planetary systems appear to occur around only a few per cent of pulsars and main-sequence stars, and are thus a rare phenomenon irrespective of circumstellar environment.     

Debris disc stirring by secular perturbations from giant planets

Detectable debris discs are thought to require dynamical excitation ('stirring'), so that planetesimal collisions release large quantities of dust. We investigate the effects of the secular perturbations of a planet, which may lie at a significant distance from the planetesimal disc, to see if these perturbations can stir the disc, and if so over what time-scale. The secular perturbations cause orbits at different semimajor axes to precess at different rates, and after some time   t _cross  initially non-intersecting orbits begin to cross. We show that   t _cross∝ a ^9/2_disc/( m _pl e _pl a ³_pl)  , where   m _pl, e _pl  and   a _pl  are the mass, eccentricity and semimajor axis of the planet, and   a _disc  is the semimajor axis of the disc. This time-scale can be faster than that for the growth of planetesimals to Pluto's size within the outer disc. We also calculate the magnitude of the relative velocities induced among planetesimals and infer that a planet's perturbations can typically cause destructive collisions out to 100 s of au. Recently formed planets can thus have a significant impact on planet formation in the outer disc which may be curtailed by the formation of giant planets much closer to the star. The presence of an observed debris disc does not require the presence of Pluto-sized objects within it, since it can also have been stirred by a planet not in the disc. For the star ε Eridani, we find that the known radial velocity planet can excite the planetesimal belt at 60 au sufficiently to cause destructive collisions of bodies up to 100 km in size, on a time-scale of 40 Myr.     

Debris discs around nearby solar analogues

An unbiased search for debris discs around nearby Sun-like stars is reported. 13 G-dwarfs at 12–15 parsec distance were searched at 850 μm wavelength, and a disc is confirmed around HD 30495. The estimated dust mass is  0.008 M_⊕  with a net limit  ≲0.0025 M_⊕  for the average disc of the other stars. The results suggest there is not a large missed population of substantial cold discs around Sun-like stars – HD 30495 is a bright rather than unusually cool disc, and may belong to a few hundred Myr old population of greater dust luminosity. The far-infrared and millimetre survey data for Sun-like stars are well fitted by either steady state or stirred models, provided that typical comet belts are comparable in size to that in the Solar system.     

Time-dependent models of the structure and stability of self-gravitating protoplanetary discs

Angular momentum transport within young massive protoplanetary discs may be dominated by self-gravity at radii where the disc is too weakly ionized to allow the development of the magneto-rotational instability. We use time-dependent one-dimensional disc models, based on a local cooling time calculation of the efficiency of transport, to study the radial structure and stability (against fragmentation) of protoplanetary discs in which self-gravity is the sole transport mechanism. We find that self-gravitating discs rapidly attain a quasi-steady state in which the surface density in the inner disc is high and the strength of turbulence very low (  α∼ 10⁻³  or less inside 5 au). Temperatures high enough to form crystalline silicates may extend out to several astronomical units at early times within these discs. None of our discs spontaneously develop regions that would be unambiguously unstable to fragmentation into substellar objects, though the outer regions (beyond 20 au) of the most massive discs are close enough to the threshold that fragmentation cannot be ruled out. We discuss how the mass accretion rates through such discs may vary with disc mass and with mass of the central star, and note that a determination of the     relation for very young systems may allow a test of the model.     

On the evolution of giant protoplanets forming in circumbinary discs

We present the results of hydrodynamic simulations of Jovian mass protoplanets that form in circumbinary discs. The simulations follow the orbital evolution of the binary plus protoplanet system acting under their mutual gravitational forces, and forces exerted by the viscous circumbinary disc. The evolution involves the clearing of the inner circumbinary disc initially, so that the binary plus protoplanet system orbits within a low density cavity. Continued interaction between disc and protoplanet causes inward migration of the planet towards the inner binary. Subsequent evolution can take three distinct paths: (i) the protoplanet enters the 4 : 1 mean motion resonance with the binary, but is gravitationally scattered through a close encounter with the secondary star; (ii) the protoplanet enters the 4 : 1 mean motion resonance, the resonance breaks, and the planet remains in a stable orbit just outside the resonance; (iii) when the binary has initial eccentricity   e _bin≥ 0.2  , the disc becomes eccentric, leading to a stalling of the planet migration, and the formation of a stable circumbinary planet.
These results have implications for a number of issues in the study of extrasolar planets. The ejection of protoplanets in close binary systems provides a source of 'free-floating planets', which have been discovered recently. The formation of a large, tidally truncated cavity may provide an observational signature of circumbinary planets during formation. The existence of protoplanets orbiting stably just outside a mean motion resonance (4 : 1) in the simulations indicate that such sites may harbour planets in binary star systems, and these could potentially be observed. Finally, the formation of stable circumbinary planets in eccentric binary systems indicates that circumbinary planets may not be uncommon.     

Three-dimensional calculations of high- and low-mass planets embedded in protoplanetary discs

We analyse the non-linear, three-dimensional response of a gaseous, viscous protoplanetary disc to the presence of a planet of mass ranging from 1 Earth mass (1 M_⊕) to 1 Jupiter mass (1 M_J) by using the zeus hydrodynamics code. We determine the gas flow pattern, and the accretion and migration rates of the planet. The planet is assumed to be in a fixed circular orbit about the central star. It is also assumed to be able to accrete gas without expansion on the scale of its Roche radius. Only planets with masses   M _p≳ 0.1 M_J  produce significant perturbations in the surface density of the disc. The flow within the Roche lobe of the planet is fully three-dimensional. Gas streams generally enter the Roche lobe close to the disc mid-plane, but produce much weaker shocks than the streams in two-dimensional models. The streams supply material to a circumplanetary disc that rotates in the same sense as the orbit of the planet. Much of the mass supply to the circumplanetary disc comes from non-coplanar flow. The accretion rate peaks with a planet mass of approximately 0.1 M_J and is highly efficient, occurring at the local viscous rate. The migration time-scales for planets of mass less than 0.1 M_J, based on torques from disc material outside the Roche lobes of the planets, are in excellent agreement with the linear theory of type I (non-gap) migration for three-dimensional discs. The transition from type I to type II (gap) migration is smooth, with changes in migration times of about a factor of 2. Starting with a core which can undergo runaway growth, a planet can gain up to a few M_J with little migration. Planets with final masses of the order of 10 M_J would undergo large migration, which makes formation and survival difficult.     

Maximum-likelihood method for estimating the mass and period distributions of extrasolar planets

We investigate the distribution of mass M and orbital period P of extrasolar planets, taking account of selection effects caused by the limited velocity precision and duration of existing surveys. We fit the data on 72 planets to a power-law distribution of the form  d n = CM ^−α P ^−β(d M / M )(d P / P )  , and find  α= 0.11 ± 0.10  ,  β=−0.27 ± 0.06  for   M ≲ 10 M _J  , where   M _J  is the mass of Jupiter. The correlation coefficient between these two exponents is −0.31, indicating that uncertainties in the two distributions are coupled. We estimate that 4 per cent of solar-type stars have companions in the range  1 M _J < M < 10 M _J  ,  2_d < P < 10 yr  .     

Dynamical relaxation and massive extrasolar planets

Following the suggestion of Black that some massive extrasolar planets may be associated with the tail of the distribution of stellar companions, we investigate a scenario in which 5 N 100 planetary mass objects are assumed to form rapidly through a fragmentation process occuring in a disc or protostellar envelope on a scale of 100 au. These are assumed to have formed rapidly enough through gravitational instability or fragmentation that their orbits can undergo dynamical relaxation on a time-scale of ∼100 orbits.
Under a wide range of initial conditions and assumptions, the relaxation process ends with either (i) one potential 'hot Jupiter' plus up to two 'external' companions, i.e. planets orbiting near the outer edge of the initial distribution; (ii) one or two 'external' planets or even none at all; (iii) one planet on an orbit with a semi-major axis of 10 to 100 times smaller than the outer boundary radius of the inital distribution together with an 'external' companion. Most of the other objects are ejected and could contribute to a population of free-floating planets. Apart from the potential 'hot Jupiters', all the bound objects are on orbits with high eccentricity, and also with a range of inclination with respect to the stellar equatorial plane. We found that, apart from the close orbiters, the probability of ending up with a planet orbiting at a given distance from the central star increases with the distance. This is because of the tendency of the relaxation process to lead to collisions with the central star. The scenario we envision here does not impose any upper limit on the mass of the planets. We discuss the application of these results to some of the more massive extrasolar planets.