The migration and growth of protoplanets in protostellar discs

We investigate the gravitational interaction of a Jovian-mass protoplanet with a gaseous disc with aspect ratio and kinematic viscosity expected for the protoplanetary disc from which it formed. Different disc surface density distributions are investigated. We focus on the tidal interaction with the disc with the consequent gap formation and orbital migration of the protoplanet. Non-linear two-dimensional hydrodynamic simulations are employed using three independent numerical codes.
A principal result is that the direction of the orbital migration is always inwards and such that the protoplanet reaches the central star in a near-circular orbit after a characteristic viscous time‐scale of ∼10⁴ initial orbital periods. This is found to be independent of whether the protoplanet is allowed to accrete mass or not. Inward migration is helped by the disappearance of the inner disc, and therefore the positive torque it would exert, because of accretion on to the central star. Maximally accreting protoplanets reach about 4 Jovian masses on reaching the neighbourhood of the central star. Our results indicate that a realistic upper limit for the masses of closely orbiting giant planets is ∼5 Jupiter masses, if they originate in protoplanetary discs similar to the minimum-mass solar nebula. This is because of the reduced accretion rates obtained for planets of increasing mass.
Assuming that some process such as termination of the inner disc through a magnetospheric cavity stops the migration, the range of masses estimated for a number of close orbiting giant planets as well as their inward orbital migration can be accounted for by consideration of disc–protoplanet interactions during the late stages of giant planet formation.     

Dynamical behaviour of the primitive asteroid belt

In this paper we consider the dynamical evolution and orbital stability of objects in the asteroid belt. A simple physical model, including full gravitational perturbations from both giant planets, is used to compute the dynamical evolution of 1000 test particles simulating the primitive asteroids. The criterion of planet crossing (or close approach in the case of resonant objects) is used to reject particles from the simulation. 44 per cent of the particles survived for the whole time-span covered by the numerical integration (∼10⁹ yr).
The 4:1, 3:1 and to a lesser extent the 2:1 Kirkwood gaps are formed in ∼10⁷ yr of evolution, representing direct numerical evidence about their gravitational origin.
We found that the rms eccentricity and inclination of the sample experience a fast increase during the first 10⁶ yr. The final rms eccentricity is 0.11, ∼60 per cent smaller than the present rms eccentricity (0.17). Nevertheless, the gravitational action of the giant planets suffices to prevent the formation of large objects, allowing catastrophic collisions and the subsequent depletion of material from this zone of the Solar system. The excited eccentricity by Jupiter and Saturn may favour mutual encounters and the further increase of the relative velocities up to their present values.     

A Bayesian Kepler periodogram detects a second planet in HD 208487

An automatic Bayesian Kepler periodogram has been developed for identifying and characterizing multiple planetary orbits in precision radial velocity data. The periodogram is powered by a parallel tempering Markov chain Monte Carlo (MCMC) algorithm which is capable of efficiently exploring a multiplanet model parameter space. The periodogram employs an alternative method for converting the time of an observation to true anomaly that enables it to handle much larger data sets without a significant increase in computation time. Improvements in the periodogram and further tests using data from HD 208487 have resulted in the detection of a second planet with a period of 909⁸²₋₉₂ d, an eccentricity of 0.37^0.26_−0.20, a semimajor axis of 1.87^0.13_−0.14 au and an M sin  i = 0.45^+0.11_−0.13 M _J. The revised parameters of the first planet are period = 129.8 ± 0.4 d, eccentricity = 0.20 ± 0.09, semimajor axis = 0.51 ± 0.02 au and M sin  i = 0.41 ± 0.05  M _J. Particular attention is paid to several methods for calculating the model marginal likelihood which is used to compare the probabilities of models with different numbers of planets.     

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.     

Orbital migrations in planetesimal discs: N-body simulations and the resonant dynamical friction

We have performed N -body numerical simulations of the exchange of angular momentum between a massive planet and a 3D Keplerian disc of planetesimals. Our interest is directed at the study of the classical analytical expressions of the lineal theory of density waves, as representative of the dynamical friction in discs 'dominated by the planet' and the orbital migration of the planets with regard to this effect. By means of a numerical integration of the equations of motion, we have carried out a set of numerical experiments with a large number of particles  ( N ≥10 000)  , and planets with the mass of Jupiter, Saturn and one core mass of the giant planets in the Solar system  ( M _c=10 M_⊕)  . The torque, measured in a phase in which a 'steady forcing' is clearly measurable, yields inward migration in a minimum-mass solar disc  (Σ∼10 g cm^-2  ), with a characteristic drift time of ∼ a few 10⁶ yr. The planets predate the disc, but the orbital decay rate is not sufficient to allow accretion in a time-scale relevant to the formation of giant planets. We found reductions of the measured torque on the planet, with respect to the linear theory, by a factor of 0.38 for M _c, 0.04 for Saturn and 0.01 for Jupiter, due to the increase in the perturbation on the disc. The behaviour of planets whose mass is larger than M _c is similar to the one of type II migrators in gaseous discs. Our results suggest that, in a minimum mass, solar planetesimals disc, type I migrations occur for masses smaller than M _c, whereas for this mass value it could be a transition zone between the two types of migration.     

The total number of giant planets in debris discs with central clearings

Infrared spectra from the Spitzer Space Telescope ( SSC ) of many debris discs are well fit with a single blackbody temperature which suggest clearings within the disc. We assume that clearings are caused by orbital instability in multiple planet systems with similar configurations to our own. These planets remove dust-generating planetesimal belts as well as dust generated by the outer disc that is scattered or drifts into the clearing. From numerical integrations, we estimate a minimum planet spacing required for orbital instability (and so planetesimal and dust removal) as a function of system age and planet mass. We estimate that a 10⁸ yr old debris disc with a dust disc edge at a radius of 50 au hosted by an A star must contain approximately five Neptune mass planets between the clearing radius and the iceline in order to remove all primordial objects within it. We infer that known debris disc systems contain at least a fifth of a Jupiter mass in massive planets. The number of planets and spacing required is insensitive to the assumed planet mass. However, an order of magnitude higher total mass in planets could reside in these systems if the planets are more massive.     

Why are there so few hot Jupiters?

We use numerical simulations to model the migration of massive planets at small radii and compare the results with the known properties of 'hot Jupiters' (extrasolar planets with semimajor axes   a < 0.1  au). For planet masses   M _pl sin  i > 0.5 M _J  , the evidence for any 'pile-up' at small radii is weak (statistically insignificant), and although the mass function of hot Jupiters is deficient in high-mass planets as compared to a reference sample located further out, the small sample size precludes definitive conclusions. We suggest that these properties are consistent with disc migration followed by entry into a magnetospheric cavity close to the star. Entry into the cavity results in a slowing of migration, accompanied by a growth in orbital eccentricity. For planet masses in excess of 1 Jupiter mass we find eccentricity growth time-scales of a few ×10⁵ yr, suggesting that these planets may often be rapidly destroyed. Eccentricity growth appears to be faster for more massive planets which may explain changes in the planetary mass function at small radii and may also predict a pile-up of lower mass planets, the sample of which is still incomplete.     

A Bayesian periodogram finds evidence for three planets in HD 11964

A Bayesian multiplanet Kepler periodogram has been developed for the analysis of precision radial velocity data. The periodogram employs a parallel tempering Markov chain Monte Carlo algorithm. The HD 11964 data have been re-analysed using 1, 2, 3 and 4 planet models. Assuming that all the models are equally probable a priori, the three planet model is found to be ≥600 times more probable than the next most probable model which is a two planet model. The most probable model exhibits three periods of  38.02^+0.06_−0.05, 360⁺⁴₋₄ and 1924⁺⁴⁴₋₄₃ d  , and eccentricities of  0.22^+0.11_−0.22, 0.63^+0.34_−0.17 and 0.05^+0.03_−0.05  , respectively. Assuming the three signals (each one consistent with a Keplerian orbit) are caused by planets, the corresponding limits on planetary mass ( M sin  i ) and semimajor axis are     respectively. The small difference (1.3σ) between the 360-d period and one year suggests that it might be worth investigating the barycentric correction for the HD 11964 data.     

Excitation of oscillation modes by tides in close binaries: constraints on stellar and orbital parameters

The parameter space favourable for the resonant excitation of free oscillation modes by dynamic tides in close binary components is explored using qualitative considerations to estimate the order of magnitude of the tidal force and the frequency range covered by the tidally induced oscillations. The investigation is valid for slowly rotating stars with masses in the interval between 2 and  20 M_⊙  , and an evolutionary stage ranging from the beginning to the end of the main sequence. Oscillation modes with eigenfrequencies of the order of five times the inverse of the dynamical time-scale  τ_dyn  of the star, i.e. the lowest-order p -modes, the f -mode and the lowest-order g ⁺-modes, are found to be outside the favourable parameter space since their resonant excitation requires orbital eccentricities that are too high for the binary to stay detached when the components pass through the periastron of their relative orbit. Resonances between dynamic tides and g ⁺-modes with frequencies of the order of half of the inverse of the dynamical time-scale of the star on the other hand are found to be favourable for orbital periods up to  ∼200τ_dyn  , provided that the binary mass ratio q is larger than 1/3, and the orbital eccentricity e is larger than ∼0.25. This favourable range comes down to orbital periods of up to 5–12 d in the case of  2–20 M_⊙  zero-age main-sequence binary components, and orbital periods of up to 21–70 d in the case of terminal main-sequence binary components.     

The potential for Earth-mass planet formation around brown dwarfs

Recent observations point to the presence of structured dust grains in the discs surrounding young brown dwarfs, thus implying that the first stages of planet formation take place also in the substellar regime. Here, we investigate the potential for planet formation around brown dwarfs and very low-mass stars according to the sequential core accretion model of planet formation. We find that, for a brown dwarf mass 0.05 M_⊙, our models predict a maximum planetary mass of  ∼5   M_⊕  , orbiting with semimajor axis ∼ 1 au. However, we note that the predictions for the mass–semimajor axis distribution are strongly dependent upon the models chosen for the disc surface density profiles and the assumed distribution of disc masses. In particular, if brown dwarf disc masses are of the order of a few Jupiter masses, Earth-mass planets might be relatively frequent, while if typical disc masses are only a fraction of Jupiter mass, we predict that planet formation would be extremely rare in the substellar regime. As the observational constraints on disc profiles, mass dependencies and their distributions are poor in the brown dwarf regime, we advise caution in validating theoretical models only on stars similar to the Sun and emphasize the need for observational data on planetary systems around a wide range of stellar masses. We also find that, unlike the situation around solar-like stars, Type II migration is totally absent from the planet formation process around brown dwarfs, suggesting that any future observations of planets around brown dwarfs would provide a direct measure of the role of other types of migration.     

The orbit and mass of the Sagittarius dwarf galaxy

Possible orbital histories of the Sgr dwarf galaxy are explored. A special-purpose N -body code is used to construct the first models of the Milky Way–Sgr dwarf system in which both the Milky Way and the Sgr dwarf are represented by full N -body systems and followed for a Hubble time. These models are used to calibrate a semi-analytic model of the Sgr dwarf's orbit that enables us to explore a wider parameter space than is accessible to the N -body models. We conclude that the extant data on the Sgr dwarf are compatible with a wide range of orbital histories. At one extreme the Sgr dwarf initially possesses ∼10¹¹ M_⊙ and starts from a Galactocentric distance R _D(0)≳200 kpc. At the other extreme the Sgr dwarf starts with ∼10⁹ M_⊙ and R _D(0)∼60 kpc, similar to its present apocentric distance. In all cases the Sgr dwarf is initially dark matter dominated and the current velocity dispersion of the Sgr dwarf's dark matter is tightly constrained to be 21±2 km s⁻¹. This number is probably compatible with the smaller measured dispersion of the Sgr dwarf's stars because of (i) the dynamical difference between dark and luminous matter, and (ii) velocity anisotropy.     

Expected planets in globular clusters

We argue that all transient searches for planets in globular clusters have a very low detection probability. Planets of low-metallicity stars typically do not reside at small orbital separations. The dependence of planetary system properties on metallicity is clearly seen when the quantity   I _e ≡ M _p[ a (1 − e )]²  is considered;   M _p, a   and e are the planet mass, semimajor axis and eccentricity, respectively. In high-metallicity systems, there is a concentration of systems at high and low values of I _e , with a low-populated gap near   I _e ∼ 0.3 M _J au²  , where M _J is Jupiter's mass. In low-metallicity systems, the concentration is only at the higher range of I _e , with a tail to low values of I _e . Therefore, it is still possible that planets exist around main-sequence stars in globular clusters, although at small numbers because of the low metallicity, and at orbital periods of ≳10 d. We discuss the implications of our conclusions on the role that companions can play in the evolution of their parent stars in globular clusters, for example, influencing the distribution of horizontal branch stars on the Hertzsprung–Russell diagram of some globular clusters, and in forming low-mass white dwarfs.     

High-eccentricity planets from the Anglo-Australian Planet Search

We report Doppler measurements of the stars HD 187085 and HD 20782 which indicate two high eccentricity low-mass companions to the stars. We find HD 187085 has a Jupiter-mass companion with a ∼1000-d orbit. Our formal 'best-fitting' solution suggests an eccentricity of 0.47, however, it does not sample the periastron passage of the companion and we find that orbital solutions with eccentricities between 0.1 and 0.8 give only slightly poorer fits (based on rms and  χ²_ν  ) and are thus plausible. Observations made during periastron passage in 2007 June should allow for the reliable determination of the orbital eccentricity for the companion to HD 187085. Our data set for HD 20782 does sample periastron and so the orbit for its companion can be more reliably determined. We find the companion to HD 20782 has   M sin   i = 1.77 ± 0.22  M _Jup  , an orbital period of 595.86 ± 0.03 d and an orbit with an eccentricity of 0.92 ± 0.03. The detection of such high-eccentricity (and relatively low-velocity amplitude) exoplanets appears to be facilitated by the long-term precision of the Anglo-Australian Planet Search. Looking at exoplanet detections as a whole, we find that those with higher eccentricity seem to have relatively higher velocity amplitudes indicating higher mass planets and/or an observational bias against the detection of high-eccentricity systems.     

Dynamical relaxation and the orbits of low-mass extrasolar planets

We consider the evolution of a system containing a population of massive planets formed rapidly through a fragmentation process occurring on a scale on the order of 100 au and a lower mass planet that assembles in a disc on a much longer time-scale. During the formation phase, the inner planet is kept on a circular orbit owing to tidal interaction with the disc, while the outer planets undergo dynamical relaxation. Interaction with the massive planets left in the system after the inner planet forms may increase the eccentricity of the inner orbit to high values, producing systems similar to those observed.     

Do binaries in clusters form in the same way as in the field?

We examine the dynamical destruction of binary systems in star clusters of different densities. We find that at high densities  (10⁴– 10⁵ M_⊙ pc⁻³)  almost all binaries with separations  >10³  au are destroyed after a few crossing times. At low densities [     ], many binaries with separations  >10³  au are destroyed, and no binaries with separations  >10⁴  au survive after a few crossing times. Therefore, the binary separations in clusters can be used as a tracer of the dynamical age and past density of a cluster.
We argue that the central region of the Orion nebula cluster was ∼100 times denser in the past with a half-mass radius of only 0.1–0.2 pc as (i) it is expanding, (ii) it has very few binaries with separations  >10³  au and (iii) it is well mixed and therefore dynamically old.
We also examine the origin of the field binary population. Binaries with separations  <10²  au are not significantly modified in any cluster, therefore at these separations the field reflects the sum of all star formation. Binaries with separations in the range  10²– 10⁴  au are progressively more and more heavily affected by dynamical disruption in increasingly dense clusters. If most star formation is clustered, these binaries must be overproduced relative to the field. Finally, no binary with a separation  >10⁴  au can survive in any cluster and so must be produced by isolated star formation, but only if all isolated star formation produces extremely wide binaries.     

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.     

Extrasolar planets and the rotation and axisymmetric mass-loss of evolved stars

I examine the implications of the recently found extrasolar planets on the planet-induced axisymmetric mass-loss model for the formation of elliptical planetary nebulae (PNe). This model attributes the low departure from spherical mass-loss of upper asymptotic giant branch (AGB) stars to envelope rotation which results from deposition of orbital angular momentum of the planets. Since about half of all PNe are elliptical, i.e., have low equatorial to polar density contrast, it was predicted that about 50 per cent of all Sun-like stars have Jupiter-like planets around them, i.e., a mass about equal to that of Jupiter, M _J, or more massive. In the light of the new findings that only 5 per cent of Sun-like stars have such planets, and a newly proposed mechanism for axisymmetric mass-loss, the cool magnetic spots model, I revise this prediction. I predict that indeed ∼50 per cent of PN progenitors do have close planets around them, but the planets can have much lower masses, as low as ∼0.01 M _J, in order to spin-up the envelopes of AGB stars efficiently. To support this claim, I follow the angular momentum evolution of single stars with main-sequence mass in the range of 1.3–2.4 M_⊙ , as they evolve to the post-AGB phase. I find that single stars rotate much too slowly to possess any significant non-spherical mass-loss as they reach the upper AGB. It seems, therefore, that planets, in some cases even Earth-like planets, are sufficient to spin-up the envelope of these AGB stars for them to form elliptical PNe. The prediction that on average several such planets orbit each star, as in the Solar system, still holds.     

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.     

Substellar companions and isolated planetary-mass objects from protostellar disc fragmentation

Self-gravitating protostellar discs are unstable to fragmentation if the gas can cool on a time-scale that is short compared with the orbital period. We use a combination of hydrodynamic simulations and N -body orbit integrations to study the long-term evolution of a fragmenting disc with an initial mass ratio to the star of   M _disc/ M _*= 0.1  . For a disc that is initially unstable across a range of radii, a combination of collapse and subsequent accretion yields substellar objects with a spectrum of masses extending (for a Solar-mass star) up to  ≈0.01 M_⊙  . Subsequent gravitational evolution ejects most of the lower mass objects within a few million years, leaving a small number of very massive planets or brown dwarfs in eccentric orbits at moderately small radii. Based on these results, systems such as HD 168443 – in which the companions are close to or beyond the deuterium burning limit – appear to be the best candidates to have formed via gravitational instability. If massive substellar companions originate from disc fragmentation, while lower-mass planetary companions originate from core accretion, the metallicity distribution of stars which host massive substellar companions at radii of ∼1 au should differ from that of stars with lower mass planetary companions.     

Weak lensing signal in unified dark matter models

This paper is an extension of the work done by Pierens & Nelson in which they investigated the behaviour of a two-planet system embedded in a protoplanetary disc. They put a Jupiter mass gas giant on the internal orbit and a lower mass planet on the external one. We consider here a similar problem taking into account a gas giant with mass in the range 0.5 to  1 M _J  and a Super-Earth (i.e. a planet with mass  ≤10 M _⊕  ) as the outermost planet. By changing disc parameters and planet masses, we have succeeded in getting the convergent migration of the planets which allows for the possibility of their resonant locking. However, in the case in which the gas giant has the mass of Jupiter, before any mean-motion first-order commensurability could be achieved, the Super-Earth is caught in a trap when it is very close to the edge of the gap opened by the giant planet. This confirms the result obtained by Pierens & Nelson in their simulations. Additionally, we have found that, in a very thin disc, an apsidal resonance is observed in the system if the Super-Earth is captured in the trap. Moreover, the eccentricity of the small planet remains low, while the eccentricity of the gas giant increases slightly due to the imbalance between Lindblad and corotational resonances. We have also extended the work of Pierens & Nelson by studying analogous systems in which the gas giant is allowed to take sub-Jupiter masses. In this case, after conducting an extensive survey over all possible parameters, we have succeeded in getting the 1:2 mean-motion resonant configuration only in a disc with low aspect ratio and low surface density. However, the resonance is maintained just for a few thousand orbits. Thus, we conclude that for typical protoplanetary discs the mean-motion commensurabilities are rare if the Super-Earth is located on the external orbit relative to the gas giant.