Planetary migration to large radii

There is evidence for the existence of massive planets at orbital radii of several hundred au from their parent stars where the time-scale for planet formation by core accretion is longer than the disc lifetime. These planets could have formed close to their star and then migrated outwards. We consider how the transfer of angular momentum by viscous disc interactions from a massive inner planet could cause significant outward migration of a smaller outer planet. We find that it is in principle possible for planets to migrate to large radii. We note, however, a number of effects which may render the process somewhat problematic.     

Predictions for the frequency and orbital radii of massive extrasolar planets

We investigate the migration of massive extrasolar planets caused by gravitational interaction with a viscous protoplanetary disc. We show that a model in which planets form at 5 au at a constant rate, before migrating, leads to a predicted distribution of planets that is a steeply rising function of log( a ), where a is the orbital radius. Between 1 and 3 au, the expected number of planets per logarithmic interval in a roughly doubles. We demonstrate that, once selection effects are accounted for, this is consistent with current data, and then extrapolate the observed planet fraction to masses and radii that are inaccessible to current observations. In total, approximately 15 per cent of stars targeted by existing radial velocity searches are predicted to possess planets with masses  0.3< M _p sin( i )<10 M _J  and radii  0.1< a <5 au  . A third of these planets (around 5 per cent of the target stars) lie at the radii most amenable to detection via microlensing. A further  5–10  per cent of stars could have planets at radii of  5< a <8 au  that have migrated outwards. We discuss the probability of forming a system (akin to the Solar system) in which significant radial migration of the most massive planet does not occur. Approximately  10–15  per cent of systems with a surviving massive planet are estimated to fall into this class. Finally, we note that a smaller fraction of low-mass planets than high-mass planets is expected to survive without being consumed by the star. The initial mass function for planets is thus predicted to rise more steeply towards small masses than the observed mass function.     

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.     

Numerical simulations of type III planetary migration – III. Outward migration of massive planets

We present a numerical study of rapid, so-called type III migration for Jupiter-sized planets embedded in a protoplanetary disc. We limit ourselves to the case of outward migration, and study in detail its evolution and physics, concentrating on the structure of the corotation and circumplanetary regions, and processes for stopping migration. We also consider the dependence of the migration behaviour on several key parameters. We perform this study using global, two-dimensional hydrodynamical simulations with adaptive mesh refinement. We find that the outward-directed type III migration can be started if the initial conditions support the initial average non-dimensional migration rate bigger than one. Unlike the inward-directed migration, in the outward migration the migration rate increases due to the growing of the volume of the co-orbital region. We find the migration to be strongly dependent on the rate of the mass accumulation in the circumplanetary disc, leading to two possible regimes of migration, fast and slow. The structure of the co-orbital region and the stopping mechanism differs between these two regimes.     

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.     

Episodic accretion in magnetically layered protoplanetary discs

We study protoplanetary disc evolution assuming that angular momentum transport is driven by gravitational instability at large radii, and magnetohydrodynamic (MHD) turbulence in the hot inner regions. At radii of the order of 1 au such discs develop a magnetically layered structure, with accretion occurring in an ionized surface layer overlying quiescent gas that is too cool to sustain MHD turbulence. We show that layered discs are subject to a limit cycle instability, in which accretion on to the protostar occurs in ∼10⁴-yr bursts with ̇ ∼10⁻⁵ M_⊙ yr⁻¹, separated by quiescent intervals lasting ∼10⁵ yr where ̇ ≈10⁻⁸ M_⊙ yr⁻¹. Such bursts could lead to repeated episodes of strong mass outflow in young stellar objects. The transition to this episodic mode of accretion occurs at an early epoch ( t ≪1 Myr), and the model therefore predicts that many young pre-main-sequence stars should have low rates of accretion through the inner disc. At ages of a few Myr, the discs are up to an order of magnitude more massive than the minimum-mass solar nebula, with most of the mass locked up in the quiescent layer of the disc at r ∼1 au. The predicted rate of low-mass planetary migration is reduced at the outer edge of the layered disc, which could lead to an enhanced probability of giant planet formation at radii of 1–3 au.     

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  .     

The brown dwarf desert as a consequence of orbital migration

We show that the dearth of brown dwarfs in short-period orbits around Solar-mass stars – the brown dwarf desert – can be understood as a consequence of inward migration within an evolving protoplanetary disc. Brown dwarf secondaries forming at the same time as the primary star have masses which are comparable to the initial mass of the protoplanetary disc. Subsequent disc evolution leads to inward migration, and destruction of the brown dwarf, via merger with the star. This is in contrast with massive planets, which avoid this fate by forming at a later epoch when the disc is close to being dispersed. Within this model, a brown dwarf desert arises because the mass at the hydrogen-burning limit is coincidentally comparable to the initial disc mass for a Solar mass star. Brown dwarfs should be found in close binaries around very low mass stars, around other brown dwarfs, and around Solar-type stars during the earliest phases of star formation.     

On the migration of a system of protoplanets

The evolution of a system consisting of a protoplanetary disc with two embedded Jupiter-sized planets is studied numerically. The disc is assumed to be flat and non-self-gravitating; this is modelled by the planar (two-dimensional) Navier–Stokes equations. The mutual gravitational interaction of the planets and the star, and the gravitational torques of the disc acting on the planets and the central star are included. The planets have an initial mass of one Jupiter mass M _Jup each, and the radial distances from the star are one and two semimajor axes of Jupiter, respectively.
During the evolution a joint wide annular gap is created by the planets. Both planets increase their mass owing to accretion of gas from the disc: after about 2500 orbital periods of the inner planet it has reached a mass of 2.3  M _Jup, while the outer planet has reached a mass of 3.2  M _Jup. The net gravitational torques exerted by the disc on the planets result in an inward migration of the outer planet on time-scales comparable to the viscous evolution time of the disc. The semimajor axis of the inner planet remains constant as there is very little gas left in its vicinity to induce any migration. When the distance of close approach eventually becomes smaller than the mutual Hill radius, the eccentricities increase strongly and the system may become unstable.
If disc depletion occurs rapidly enough before the planets come too close to each other, a stable system similar to our own Solar system may remain. Otherwise the orbits may become unstable and produce systems like υ And.     

Metallicity, planetary formation and migration

Recent observations show a clear correlation between the probability of hosting a planet and the metallicity of the parent star. As radial velocity surveys are biased, however, towards detecting planets with short orbital periods, the probability–metallicity correlation could merely reflect a dependence of migration rates on metallicity. We investigated the possibility, but find no basis to suggest that the migration process is sensitive to the metallicity. The indication is, therefore, that a higher metallicity results in a higher probability for planet  formation .     

Do accretion discs regulate the rotation of young stars?

We present a photometric study of I -band variability in the young cluster IC 348. The main purpose of the study was to identify periodic stars. In all, we find 50 periodic stars, of which 32 were previously unknown. For the first time in IC 348, we discover periods in significant numbers of lower-mass stars  ( M < 0.25 M_⊙)  and classical T Tauri stars. This increased sensitivity to periodicities is a result of the enhanced depth and temporal density of our observations, compared with previous studies. The period distribution is at first glance similar to that seen in the Orion nebula cluster (ONC), with the higher-mass stars  ( M > 0.25 M_⊙)  showing a bi-modal period distribution concentrated around periods of 2 and 8 d, and the lower-mass stars showing a uni-modal distribution, heavily biased towards fast rotators. Closer inspection of the period distribution shows that the higher-mass stars show a significant dearth of fast rotators, compared to the ONC, whilst the low-mass stars are rotating significantly faster than those in Orion. We find no correlation between rotation period and K – L colour or Hα equivalent width.
We also present a discussion of our own IC 348 data in the context of previously published period distributions for the ONC, the Orion flanking fields and NGC 2264. We find that the previously claimed correlation between infrared excess and rotation period in the ONC might, in fact, result from a correlation between infrared excess and mass. We also find a marked difference in period distributions between NGC 2264 and IC 348, which presents a serious challenge to the disc-locking paradigm, given the similarity in ages and disc fractions between the two clusters.     

Evidence of a high-velocity ionized outflow in a second narrow-line quasar PG 0844+349

Following the discovery of X-ray absorption in a high-velocity outflow from the bright quasar PG 1211 + 143 we have searched for similar features in XMM–Newton archival data of a second (high accretion rate) quasar PG 0844+349. Evidence is found for several faint absorption lines in both the EPIC and RGS spectra, whose most likely identification with resonance transitions in H-like Fe, S and Ne implies an origin in highly ionized matter with an outflow velocity of order ∼0.2c. The line equivalent widths require a line-of-sight column density of   N _H∼ 4 × 10²³ cm⁻²  , at an ionization parameter of log  ξ∼ 3.7  . Assuming a radial outflow being driven by radiation pressure from the inner accretion disc, as suggested previously for PG 1211 + 143, the flow in PG 0844+349 is also likely to be optically thick, in this case within ∼25 Schwarzschild radii. Our analysis suggests that a high-velocity, highly ionized outflow is likely to be a significant component in the mass and energy budgets of active galactic nuclei accreting at or above the Eddington rate.