A&A 461, 1195-1208 (2007)
DOI: 10.1051/0004-6361:20066171
M. J. Fogg - R. P. Nelson
Astronomy Unit, Queen Mary, University of London, Mile End Road, London E1 4NS, UK
Received 3 August 2006 / Accepted 9 October 2006
Abstract
Context. There are numerous extrasolar giant planets which orbit close to their central stars. These "hot-Jupiters'' probably formed in the outer, cooler regions of their protoplanetary disks, and migrated inward to 0.1 AU. Since these giant planets must have migrated through their inner systems at an early time, it is uncertain whether they could have formed or retained terrestrial planets.
Aims. We present a series of calculations aimed at examining how an inner system of planetesimals/protoplanets, undergoing terrestrial planet formation, evolves under the influence of a giant planet undergoing inward type II migration through the region bounded between 5-0.1 AU.
Methods. We have previously simulated the effect of gas giant planet migration on an inner system protoplanet/planetesimal disk using a N-body code which included gas drag and a prescribed migration rate. We update our calculations here with an improved model that incorporates a viscously evolving gas disk, annular gap and inner-cavity formation due to the gravitational field of the giant planet, and self-consistent evolution of the giant's orbit.
Results. We find that 60% of the solids disk survives by being scattered by the giant planet into external orbits. Planetesimals are scattered outward almost as efficiently as protoplanets, resulting in the regeneration of a solids disk where dynamical friction is strong and terrestrial planet formation is able to resume. A simulation that was extended for a few Myr after the migration of the giant planet halted at 0.1 AU, resulted in an apparently stable planet of
2
forming in the habitable zone. Migration-induced mixing of volatile-rich material from beyond the "snowline'' into the inner disk regions means that terrestrial planets that form there are likely to be water-rich.
Conclusions. We predict that hot-Jupiter systems are likely to harbor water-abundant terrestrial planets in their habitable zones. These planets may be detected by future planet search missions.
Key words: planets and satellites: formation - methods: N-body simulations - astrobiology
About one quarter of the extrasolar planetary systems discovered to
date contain a so-called "hot-Jupiter'' - a gas giant planet
orbiting within 0.1 AU of the central star (e.g. Butler et al. 2006).
It is improbable that these planets formed within such a hot region
of the protoplanetary disk, and is most likely that they originated
further out beyond the nebula snowline and moved inward via type II
migration (e.g Lin & Papaloizou 1986; Nelson et al. 2000; Ward 1997; Bryden et al. 1999). Since both giant
planet formation and type II migration require the nebular gas to
still be present, these processes are constrained to occur within
the first few million years of the disk lifetime (Haisch et al. 2001),
well within the 30-100 Myr timescale thought to be
required to complete the growth of terrestrial planets
(e.g. Halliday 2004; Chambers 2001; Kleine et al. 2002). Thus, these migrating
giant planets must have traversed the inner system, including its
habitable zone, before any planet formation there was complete,
raising the question of what effect such a disturbance would have on
the growth of terrestrial-like planets.
Initially, it was thought reasonable that the inward migration of a giant planet would be so disruptive of the material it passed through as to clear the swept zone completely, precluding the formation of any inner system planets. However, this view was a conservative assumption, often made in the support of speculative astrobiological arguments (e.g. Ward & Brownlee 2000; Lineweaver et al. 2004; Lineweaver 2001). Agreement as to the outcome failed to materialize from the first modeling studies of the process, the conclusions of which varied from the occurrence of terrestrial planets in hot-Jupiter systems being highly unlikely (Armitage 2003), through possible but rare (Mandell & Sigurdsson 2003), to commonplace (Raymond et al. 2005a). None of these studies, however, actually simulated terrestrial planet formation simultaneously with giant planet migration. Their disagreement about the likely outcome can be traced to assumptions made about the timing of giant planet formation and migration.
The first study to model inner system planetary accretion in the presence of a migrating giant was that of Fogg & Nelson (2005) (hereafter referred to as Paper I). This work used N-body simulations to examine oligarchic and giant-impact growth (Kokubo & Ida 2000) in a protoplanet/planetesimal disk based on the minimum mass solar nebula model of Hayashi (1981), extending between 0.4-4 AU. Five scenarios were considered, corresponding to five different ages for the inner planet forming disk at the point when a giant planet was assumed to form at 5 AU and migrate in to 0.1 AU.
In all five of their scenarios, Fogg & Nelson (2005) found that the majority
of the disk solids survived the passage of the giant planet, either
by being shepherded inward of the giant, or by being scattered by
the giant into excited exterior orbits. This partition of solid
material was shown to vary with the level of dissipative forces
present, which decline with disk maturity, favoring shepherding at
early times and scattering at late times. Within the portion of the
disk compacted inside the increasingly restricted volume interior to
the giant, accretion was found to speed up greatly at late times
resulting typically in a
planet
forming inside 0.1 AU. The similarity of these objects to the
recently identified class of "hot-Neptune'' planets
(Bonfils et al. 2005; Udry et al. 2006; McArthur et al. 2004; Butler et al. 2004; Vogt et al. 2005; Santos et al. 2004b; Rivera et al. 2005) was
noted and discussed. The fate of the material scattered into
external orbits was not subjected to further calculation, but it was
noted that, although the remaining solids surface density was
reduced from pre-migration-episode values, ample material remained
to provide for the eventual accretion of a set of external
terrestrial planets, including within each system's habitable zone.
Fogg & Nelson (2005) therefore concluded that the assumption that
hot-Jupiter systems are devoid of inner system terrestrial planets
is probably incorrect, and that planet formation and retention both
interior and exterior to the hot-Jupiter is possible. We note that
similar results, relating to the formation of planets interior to a
migrating giant, have been reported by Zhou et al. (2005).
In this paper, we extend the model introduced in Paper I by improving the realism with which gas dynamics is simulated. A 1D evolving viscous gas disk model is linked to the N-body code that: 1) allows the gas to deplete over time via viscous accretion onto the central star; 2) allows an annular gap to form in the vicinity of the giant planet; 3) includes the creation of an inner cavity due to dissipation of propagating spiral waves excited by the giant planet; 4) self-consistently drives the giant inward. Compared to the unevolving gas disk assumed in Paper I, this new model reduces the strength of dissipation present in all scenarios, especially in regions close to the central star and the giant. We examine and re-interpret the fate of the disk solids under these changed circumstances and look at whether our hypothesis of hot-Neptune formation remains robust. We also examine the post-migration evolution of the outer scattered disk of solids, and find that terrestrial planets do form in the habitable zone. Another issue we examine is the extent to which volatile-rich matter, originating from beyond the snowline, is driven into the inner system and is mixed into the surviving planet-forming material. We find that substantial mixing occurs, such that any terrestrial planets that form are likely to be water-rich bodies hosting deep, global oceans.
The plan of the paper is as follows. In Sect. 2 we outline the additions to our model and the initial conditions of the simulations; in Sect. 3 the results are presented and discussed; in Sect. 4 we consider some caveats and future model improvements, and in Sect. 5 we offer our conclusions.
We model planetary accretion using the Mercury 6 hybrid-symplectic integrator (Chambers 1999), run as an N + N'simulation, where we have N protoplanets embedded in a disk of N' "super-planetesimals'' - particles that represent an idealized ensemble of a much larger number of real planetesimals. The protoplanets (and the giant when it is introduced) interact gravitationally and can accrete with all the other bodies in the simulation, whereas the super-planetesimal population is non-self-interacting. These latter objects however are subject to a drag force from their motion relative to the nebular gas. A detailed outline of these aspects of our model is given in Sects. 2.1 and 2.2 of Paper I and we continue here to discuss the additional features we have introduced.
In Paper I we assumed a steady state gas disk model with a constant
surface density profile r-1.5. The migration rate of
the giant was prescribed from a calculation of the local viscous
evolution timescale. More realistically, the quantity of nebular gas
should decline due to viscous evolution and accretion onto the
central star, progressively depleting the inner disk. The gas
density should also decrease in the vicinity of the giant due to the
generation and dissipation of density waves at Lindblad resonance
positions, clearing an annular gap in a zone where the planetary
tidal torques dominate the intrinsic viscous torques of the disk. A
consistent calculation of the type II migration rate should involve
the back-reaction to these tidal torques.
To account for these processes we model the gas disk by solving numerically the disk viscous diffusion equation (Pringle 1981), with modifications included to account for the tidal influence of an embedded giant planet. Such a method has been used previously in studies that attempt to explain the statistical distribution of exoplanetary orbits through type II migration and disk dispersal (Alibert et al. 2005; Armitage et al. 2002; Trilling et al. 1998). The simplest technique for including the effect of the planet is the impulse approximation of Lin & Papaloizou (1986), where wave dissipation is assumed to occur close to the planet. A more sophisticated treatment of the problem is the WKB approximation (Takeuchi et al. 1996) which involves summing the torque contributions from a series of Lindblad resonances in the disk. The former technique was adopted in the studies cited above as it requires considerably less computation and generates comparable results. We follow this approach here, but in order to include the effect of non-local dissipation of waves that travel far into the disk we also include the WKB approximation torques due to the waves launched at the two innermost and outermost Lindblad resonances.
We assume a MMSN-type protoplanetary disk around a
star (Hayashi 1981). The initial surface density
of solids is:
The initial surface density of gas is:
Given a nebular radial temperature profile of
,
the sound speed of a solar composition
gas in cgs units is:
We solve the diffusion equation for
in the form:
The exchange of angular momentum between the planet and disk leads to a radial migration of the planet at a rate:
The rate of specific angular momentum transfer to the disk in the impulse approximation is given by Lin & Papaloizou (1986) as:
The total torque density exerted on the disk in the WKB approximation via the damping of waves excited by the planet is (Takeuchi et al. 1996):
The torque density is calculated from the radial gradient of the
angular momentum flux Fm(r):
The parameter
is a torque cutoff function
(Artymowicz 1993) given by:
The evolution of the nebular gas is computed by solving
Eq. (7) with an explicit finite-difference technique on
a grid with a cell width
.
The resulting
type II migration forces on the giant planet are computed from
Eq. (8) by deriving an instantaneous time scale,
,
and inserting this into Eq. (6) of Paper I. Strong eccentricity
damping for the giant planet is assumed, with the damping time scale
being 1/50th of the radial migration timescale. As in Paper I, we
have neglected the effects of type I migration and associated
eccentricity damping (Tanaka & Ward 2004; Tanaka et al. 2002; Ward 1997).
The gas disk adopted here extends from an inner radius
to an outer radius of
with an initial
profile given by Eq. (2). Since we are
considering a
disk,
and
the initial disk mass is
.
Note that this initial disk mass is
greater than the
of gas assumed in Paper I;
however since that amount was kept fixed, the gas present in the
simulations presented here falls below this level after
140 000 years. The boundary conditions for the
computation were
,
representing gas accretion onto the central star and at
the radial velocity of the gas was set to zero.
To correctly couple the evolving gas disk algorithm with the N-body
code, synchronization of their respective time-steps is necessary.
In each simulation sub-run (see Sect. 2.3) the symplectic
N-body time-step
was fixed whereas the gas
disk time-step
is adaptive and taken to be:
At early times,
the solid component of a young protoplanetary nebula will exhibit a
radial pattern of chemical composition, controlled by the
temperature-dependent condensation sequence of a variety of metals,
rock minerals and ices. As the planetary system grows and evolves,
phenomena such as dynamical spreading, gas drag induced orbital
decay and resonant interactions can cause a radial mixing of
material. According to one school of thought
(e.g. Morbidelli et al. 2000) the original matter that condensed in the
Earth's orbit is thought to have been dry (for an alternative
opinion see Drake & Righter 2002) and the origin of the Earth's water, and its
D/H ratio, can be explained if 10% of the planet's mass is
composed of carbonaceous chondrite-type material, originating from
between 2.5-4 AU, and
10% of the water gained thereby is
retained at the end of accretion. N-body simulations of terrestrial
planet formation from disks that extend out close to the orbit of
Jupiter are supportive of this idea and all demonstrate substantial
mixing of water rich material into the inner disk
(e.g. Raymond et al. 2004,2005a; Chambers 2001; Raymond et al. 2005b; Morbidelli et al. 2000).
In Paper I it was noted that one consequence of the inward migration of a giant planet is the shepherding of planetesimals that are damped by gas drag (Tanaka & Ida 1999), and the trapping of bodies at first order mean motion resonances. Hence the outer, more volatile-rich, regions of the protoplanetary disk are actively mixed into its inner regions. However, the composition of planetesimals and protoplanets, and their accretion products, were not logged in our simulations.
We amend this here by labeling all particles with a composition based on its original position in the disk, and summing the composition of protoplanets as they grow. In a similar manner to most other studies, we assume a crude three-phase initial radial composition with rocky material originating at <2 AU, material characteristic of chondritic meteorites between 2-2.7 AU, and trans-snowline material at >2.7 AU. For convenience, these phases are referred to as "dry'', "damp'' and "wet'' respectively but, in contrast to other studies, we do not assign an actual water mass fraction to them. Instead we use our results to make more generalized observations and predictions.
Table 1: Data describing initial solids disk set-up.
![]() |
Figure 1: Evolution of the mass of the gas disk. The mass of gas (in Jupiter masses) remaining at the launch point for each of the five migration scenarios is indicated. |
Open with DEXTER |
In common with Paper I, we assume a nominal age for the
protoplanetary disk of 0.5 Myr, this being the t = 0 start date
for our simulations, and a mass of three times that of the MMSN. Our
reasons for choosing a more massive nebula stem from the fact that
core-accretion theories of giant planet formation require an
enhanced density of solid material in order to grow a critical core
mass before the nebular gas is lost (Thommes et al. 2003; Lissauer 1987) which
might be supported by the observation that hot-Jupiters are found
predominantly around stars more metal-rich than the Sun
(Santos et al. 2004a; Fischer & Valenti 2005). As runs proceed, the initial gas disk
described by Eq. (2) evolves according to the algorithm
outlined in Sect. 2.1, initially without the presence of
a giant planet (i.e.
in Eq. (7)).
Simultaneously, Mercury 6 evolves an initial disk of solid
material, extending radially from 0.4-4.0 AU. We use the same
two-component initial solids disk as was used in Paper I, generated
according to the profile of Eq. (1) and in line with the
oligarchic growth picture of Kokubo & Ida (2000), which we assume to be a
reasonable description of the state of the inner disk at 0.5 Myr.
Data for this initial disk model are shown in Table 1
which gives, for zones interior and exterior to the snowline, values
for the total mass of solid material
,
the
number and mass of protoplanets N and
,
and
the number and mass of super-planetesimals N' and
.
Note that as super-planetesimals act as
statistical tracers for a much more numerous population of real
planetesimals, their mass is much greater than that of a real
planetesimal, but for the purposes of calculating gas drag, each
super-planetesimal is treated as if it is a single 10 km radius
object of realistic mass. The parameter
,
at the
foot of Table 1, is the mass fraction of the solids disk
contained in protoplanets and we use this here as a rough measure of
the evolution of the disk, taking
to denote
the transition between oligarchic and giant impact growth regimes.
From t = 0, we run our new model for 0.1, 0.25, 0.5, 1.0 and
1.5 Myr in the absence of the giant, with
days and a simulation inner
edge of
.
(Note that r* does not
necessarily denote the physical radius of the central star, but it
is the distance interior to which particles are removed from the
simulation and their masses added to the central star.) The
evolution of the nebular mass over this time span and the particular
nebular mass at each of these five epochs are shown in
Fig. 1. The gas surface density profiles resulting at
these times are shown in Fig. 2. It is apparent that
gas drains onto the central star very rapidly at first, as the
density gradient relaxes to a shallower profile. Compared to the
unevolved profile, order of magnitude reductions in gas density
occur within the disk's inner regions. Data for the evolved solid
components are given in Table 2 and include the values
of
,
the mass of the largest protoplanet to have
evolved in each case. (In contrast to Paper I we do not run our
model to 3.0 Myr, as there remains too little gas at this time to
provide for our giant planet's envelope.) The advance of planetary
growth with time is indicated by the progressive increase of
and
and the reduction in
particle numbers as super-planetesimals are accreted and
protoplanets merge. When compared to the equivalent data from Paper I, we note that no planetesimals are lost beyond the inner edge of
the simulation because of much lower gas densities near the central
star.
![]() |
Figure 2:
Evolution of the gas surface density within the inner 10 AU
of our simulated disk. The upper solid line is the
![]() ![]() ![]() |
Open with DEXTER |
Table 2: Overall solids disk data: after 0.1-1.5 Myr of evolution.
![]() |
Figure 3: Semi-major axis evolution of the giant planet in each scenario from the launch time (the top row in Table 2) to the time at which it arrives at 0.1 AU. |
Open with DEXTER |
The five type II migration scenarios studied here are constructed
from the five evolved disks indicated in Fig. 1 and
summarized in Table 2. Scenarios I-III take place
whilst the solids disk remains in its oligarchic growth phase
(
)
whereas Scenarios IV-V have just
entered the final giant impact stage of growth (
). A giant planet of mass
is placed into each
simulation at 5.0 AU after removing
of gas from
between 3-7 AU. The giant then proceeds to clear out an annular
gap in the disk and migrates inward according to the method
described in Sect. 2.1. The runs are halted once the
giant reaches 0.1 AU (see Fig. 3). For Scenarios I-III, this takes
110 000 - 120 000 years, close to the
prediction of Eq. (6). The process takes longer to
complete in the cases of Scenarios IV-V (
150 000 and
220 000 years respectively) because by the time of the appearance
of the giant planet the gas disk is substantially depleted and is
less effective at driving migration. In order to better model
processes when the giant migrates down to small radial distances, we
contract the simulation inner edge down to a realistic T Tauri star
radius:
.
The
initial timestep chosen for the symplectic integrator was
days, but it was necessary to reduce this
at late times as material is driven into closer orbits. Hence each
scenario was divided into a number of sequential sub-runs with
being adjusted at each re-start so as to keep
the timestep close to one tenth the orbital period of the innermost
object. Since planetesimals in this new model suffer less orbital
decay due to gas drag, it was possible to conduct these runs with a
higher value of
than in Paper I. However, this
advantage was negated at late times as the adaptive
falls steeply as the giant planet moves within
1 AU. From t = 0, the scenarios presented here each required
4-8 weeks of 2.8 GHz CPU-time for completion.
![]() |
Figure 4: Scenario I at 20 000 years after the start of giant planet migration, showing the mass, inclination and eccentricity of objects. Small black dots represent super-planetesimals; white filled circles are rocky protoplanets; grey filled circles are icy protoplanets and the large highlighted grey filled circle is the giant. Objects plotted between the dotted lines in the upper panel have orbits that intersect the orbit of the giant. The location of the 2:1, 3:2 and 4:3 resonances with the giant are indicated. Gas surface density is read on the right hand axis of the lower panel, the upper grey curve being the unevolved profile at t = 0 and the lower black curve being the current profile. |
Open with DEXTER |
![]() |
Figure 5: Scenario I at 60 000 years after the start of giant planet migration. The giant has now moved inward to 2.72 AU. Increasing excitation of the orbits of protoplanets captured at resonances is apparent, as is the build-up of matter scattered into external orbits. |
Open with DEXTER |
We begin describing the results by focussing on the evolution during the migration of the giant planet, discussing one specific case in detail before examining differences between the various runs. We then go onto describe the evolution of Scenario V after giant planet migration has halted, focussing on the issue of terrestrial planet formation in the scattered disk, and the likely composition of planets that form there.
The results of all scenarios showed a number of behavioral features in common. As in Paper I, we illustrate these first by describing the results of Scenario I in detail. We then proceed to examine how the results differ between scenarios (dependence on disk maturity) and how the results differ from those of Paper I (dependence on an evolving gas disk).
The typical effects of a migrating giant planet on an inner solids disk observed from Paper I were as follows: 1) shepherding of planetesimals; 2) capture of objects at first order mean motion resonances; 3) acceleration of accretion interior to the giant with possible hot-Neptune formation; and 4) generation of a scattered exterior disk. To a greater or lesser extent, these features were also observed in our new simulations. Four snapshots of the evolution of Scenario I are illustrated in Figs. 4-7 showing the mass, inclination and eccentricity of objects, and the gas surface density vs. semi-major axis. The original provenance of the protoplanets (interior or exterior to the snowline) is denoted by the shading of its symbol as described in the caption to Fig. 4. In the case of a merger between rocky and icy protoplanets, this shading is determined by that of the most massive of the pair.
![]() |
Figure 6:
Scenario I at 100 000 years after the start of giant planet
migration. The giant has now moved inward to 0.70 AU. Five protoplanets
are currently crossing the orbit of the giant. The scattered disk has
grown and a >
![]() |
Open with DEXTER |
An early stage in the evolution of Scenario I, 20 000 years after
the introduction of the giant planet, is shown in
Fig. 4. The giant has opened a 0.75 AU gap in
the gas and has migrated inward to 4.37 AU, shepherding the outer
disk edge at the 4:3 resonance, now at 3.61 AU. Capture of objects
at the 3:2 and 2:1 resonances, at 3.33 and 2.75 AU respectively, is
apparent from eccentricity spikes visible in the upper panel and a
clustering of protoplanets in the lower panel. Even at this early
phase, before the giant has entered the original confines of the
interior disk, one protoplanet and a handful of super-planetesimals
have been scattered into external orbits.
The system midway through the run, 60 000 years after the introduction of the giant planet, is shown in Fig. 5. The giant has now migrated to 2.72 AU, putting the positions of the 3:2 and 2:1 resonances at 2.07 and 1.71 AU respectively. Strong excitation of protoplanetary orbits is now apparent at these locations, as is the build-up of scattered material in external orbits. The primary mechanism of this expulsion is evident from the behavior of material captured at resonances. Continuous resonant pumping results in orbits becoming eccentric enough to eventually intersect the orbit of the giant planet. A series of close encounters with the giant then follows, eventually resulting in a final encounter where the object is accreted or expelled into a non-intersecting exterior orbit.
![]() |
Figure 7: Scenario I at 114 000 years after the start of giant planet migration. The giant planet has migrated to 0.1 AU. Most interior mass has been lost after the most massive interior protoplanet impacts the giant. 63% of the original solids disk mass now resides in exterior orbits. |
Open with DEXTER |
An advanced stage of Scenario I, 100 000 years after the
introduction of the giant planet is shown in Fig. 6.
The giant planet is now at 0.70 AU and the 3:2 and 2:1 resonances
are at 0.54 and 0.44 AU respectively. A substantial scattered
external disk has now formed and sufficient gas remains in this
early scenario to rapidly damp the orbits of scattered
planetesimals. An impression of the scattering process in action is
given by the five protoplanets currently crossing the giant's orbit.
The interior disk is compacted to high surface densities, but now
that strong first order resonances with the giant are influential
throughout its width, and gas densities have fallen by a factor of
>10, there is a noticeable dynamical stirring of its entire
remaining population. Nevertheless, accretion has speeded up in this
shepherded zone with the growth of one protoplanet of
1.17
at 0.53 AU outrunning that of its
neighbors.
![]() |
Figure 8:
Surface density evolution ( left hand panel) and accretion rates
( right hand panel) for Scenario I. Growing surface density peaks at the 2:1
and 3:2 resonances sweep through the inner system ahead of the giant.
Accretion rates increase after ![]() |
Open with DEXTER |
Scenario I is terminated at 114 000 years after the start of
migration when the giant planet arrives at 0.1 AU and the system is
illustrated at this point in Fig. 7. Two thirds of the
original solids disk mass has survived the migration episode - the
great majority of this residing in the scattered exterior disk. Most
of the remainder has been accreted by the giant planet. Just
of the original disk mass remains interior to the giant and
none of this is in the form of large bodies. The rapidly accreting
interior protoplanet seen in Fig. 6 continued its
inward progress close to the 3:2 resonance position with its orbit
being well-damped by strong collisional damping and dynamical
friction from planetesimals and smaller protoplanets (see
Fig. 10 in Sect. 3.1.4). However,
2000 years before its demise, the protoplanet drifted outward and
became captured at the 4:3 resonance. At this location, both
accretion and dynamical friction were reduced allowing the
protoplanet's orbit to become progressively more eccentric. At
109 220 years the protoplanet, now weighing in at
2.41
,
collided with and was accreted by the
giant. The five other less massive interior planets visible in
Fig. 6 grew very little, remaining between
,
and in due course one of them impacted the
giant and four were scattered.
To emphasize the above description of Scenario I, the surface
density evolution of the disk and its accretion rate are shown in
Fig. 8. The left hand panel shows the disk surface
density profile (obtained by summing all protoplanets and
super-planetesimals in 0.1 AU width bins) at 20 000, 40 000,
60 000, 80 000 and 100 000 years after the start of migration;
the right hand panel plots the amount of mass accreted onto
protoplanets only (including protoplanet mergers) every 1000 years
for the duration of the run. In the
plot, two
surface density enhancements are clearly visible as spikes at the
3:2 and 2:1 resonances and are seen to grow whilst moving inward. At
80 000 years, these have almost merged into one: the shepherded
portion of the original disk having by now been squeezed into a
dense ring. By 100 000 years, most of this mass is now confined
within 0.6 AU and
here has risen to
which is off the vertical scale in the
figure. This amounts to an increase by a factor of
10 over
the previous, undisturbed, disk surface density, but is only about
half the increase seen in the equivalent Scenario presented in Paper
I. The effect of this disk compaction process is visible in the
accretion rate plot. Mass accretion rises significantly after
80 000 years due to both the high values of
and the fact that much of this mass now resides in a zone where
dynamical times are shorter. However, the large, terminal, accretion
rate spike described in Paper I is not reproduced here (compare
Fig. 8 with Fig. 6 of Paper I). This is because close
to the end of that previous simulation a 15.65
hot-Neptune was assembled in a dramatic phase of runaway accretion
interior to 0.1 AU. In the case presented here the compacted
interior disk is only half as dense and is much less well damped
(note the large difference between the upper and lower gas density
curves in Fig. 7) and, whilst a protoplanet does grow
to 2.41
in this region, as described above, it
does not survive and is accreted by the giant planet.
We describe the formation and fate of interior planets
in more detail in Sect. 3.1.4.
The reason for running five scenarios through a progressively more
mature inner disk is to see if the timing of migration has any
systematic effect on the results. This is possible as the partitioning
of the solids disk between inner and outer remnants is influenced by
the level of damping that particles are subject to, which declines
with age. In Paper I, where
is fixed, this occurs as a
side effect of accretion: as planetesimals are accreted by
protoplanets (
increases), fewer small particles
remain that are subject to gas drag and which can exert dynamical
friction. In these latest simulations, since we now have an evolving
gas disk, the strength of gas drag on susceptible particles also
declines with time and is particularly marked close to the central
star and giant planet.
Table 3: Fate of the disk mass at the end of each run.
Data describing the fate of the solids disk mass at the end of each scenario are shown in Table 3. Disk mass that is lost is either accreted by the central star, ejected from the system, or accreted by the giant planet; that which survives is partitioned between bodies orbiting interior or exterior to the final orbit of the giant planet at 0.1 AU.
In all scenarios, a negligible quantity of mass was ejected or lost to the central star. However, a significant fraction of the disk mass (14-33%) was accreted by the giant, especially towards the end of the migration. At these late times, planetesimals are shepherded into the partially evacuated inner regions of the gas disk where gas drag is less effective at damping orbital perturbations from the giant and growing protoplanets. Once a planetesimal strays into the annular gap in the gas containing the giant, gas drag vanishes and accretion or scattering by the giant follows. This increased excitation of the shepherded planetesimal population, and the thinning down of their number, renders them less effective at damping protoplanetary orbits via dynamical friction and collisions. Hence, at late times the orbits of the remaining interior protoplanets also tend to destabilize, with one of the same two fates in store. A trend can be seen in Table 3 for the giant planet to accrete less material with disk maturity. This occurs because as the disk ages the gas density and the solids mass fraction in small bodies both decline, resulting in less dynamical dissipation of both planetesimals and protoplanets. Less matter is shepherded in such mature disks so there is less of an interior remnant for the giant to accrete from at late times.
In all scenarios, a large majority of the disk solids are found to
survive the migration episode - over two thirds of the original
inventory. Table 3 shows that there is essentially no
trend with disk maturity in the partitioning of surviving mass
between interior and exterior remnants. Just a few percent of the
mass remains interior to the giant in all cases. When the giant
planet migrates through a disk undergoing oligarchic growth
(
;
Scenarios I-III)
65% of the
original disk mass survives by being scattered into the exterior
disk. This fraction increases for disks undergoing giant
impact-style growth (
;
Scenarios IV-V) to
>80%, not because more mass remains at <0.1 AU but because
less mass is accreted by the giant at late times.
The salient dynamical behaviors of solids disk particles such as shepherding, resonant capture, scattering by the giant planet and eventual partition into interior and exterior remnants are observed generally in the results of both the present model and those of Paper I. However, the introduction of an evolving gas disk causes the relative predominance of these outcomes to differ. This is because both the principal sources of dissipation, dynamical friction and gas drag, fall with time, whereas in the model from Paper I only the former declines. There are five systematic differences between the results of Paper I and those presented here.
These tendencies are illustrated in Fig. 9 where the
percent of the original solids disk mass surviving at the end of the
simulations presented here are compared with the
simulation results from Paper I.
Interior and exterior remnants are plotted as dashed and solid lines
respectively; gray lines are the results of Paper I and black lines
are those of the present work. It can be seen that partitioning
varies strongly with disk maturity in the case of a steady-state gas
model and weakly, if at all, when the nebula gas is allowed to
viscously evolve. Since gas drag does not decline in the former
case, the shepherding of planetesimals is more influential and more
mass remains in the interior remnant in early scenarios, which is
mostly contained in surviving hot-Neptune-type planets. Late
scenarios from Paper I behave more similarly to the ones presented
here as a greater fraction of the solids disk mass is contained
within protoplanets, which do not feel the gas, and their source of
dissipation, dynamical friction, is weaker due to a decline in
super-planetesimal numbers.
![]() |
Figure 9:
Comparison of the results of Paper I (grey lines and ![]() ![]() ![]() |
Open with DEXTER |
The increased tendency in the present model for planetesimals to
scatter into the exterior disk can also be appreciated from
Fig. 9 by comparing, for the two models, the ratio of
the protoplanet mass fraction in the final exterior disk to that of
the original disk,
.
If the scattering process does
not discriminate between protoplanets and planetesimals, and there
is no further accretion after scattering, this ratio would remain at
.
Subsequent accretion will
raise this ratio, but not by much over the timescale considered,
especially as accretion rates are reduced in the exterior disk due
to its dynamically excited state. Reading from the right hand
y-axis, Fig. 9 shows that
for all of the present
model results whereas much higher values are obtained from Paper I.
Protoplanets are preferentially scattered into the exterior disk
when there is strong and unvarying gas drag, whereas fractionation
is much less marked when gas drag declines over time and falls to
near zero in the vicinity of the giant planet.
![]() |
Figure 10:
Blow-up of the inner 1 AU of Scenario I, showing eccentricity
vs semi-major axis 105 000 years after
the start of giant planet migration. Protoplanetary masses are indicated in
![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
We have performed additional calculations, very similar to those
presented here, in which we switched off the effect of dissipating
spiral waves travelling into the disk (i.e. T(r)=0 in
Eq. (7)). This term helps create an inner cavity in the
gas disk as the planet migrates inward. These calculations produced
interior planets ranging in mass 2-4
during
giant planet migration. A variety of final outcomes were noted: 1)
some of them were accreted by the giant planet; 2) others were
scattered externally into stable orbits at
AU; and 3) a
few remained in interior orbits, typically close to 0.076 AU (the
3:2 resonance). These results, taken alongside the five scenarios
already presented in this paper and those from Paper I show a
distinct trend: a strongly dissipating gas disk interior to the
giant planet leads to the formation of fairly massive,
hot-Neptune-like planets which survive; a gas disk of lower
density leads to the formation of lower mass interior planets that
often do not survive.
We note at this point that our models currently neglect some
potentially important sources of dissipation due to general
uncertainty about how planetary formation can proceed in their
presence, such as type I migration (Tanaka & Ward 2004; Tanaka et al. 2002; Ward 1997) and
the circularization of orbits due to stellar tides. Associated with
the former is strong eccentricity and inclination damping of low
mass planets which may facilitate the survival of inner planets.
Simulations are currently underway to examine this possibility.
Tidal damping of orbits however is unlikely to have a significant
effect on our results because all shepherded protoplanets are
scattered or accreted by the giant at
AU. Whilst
tidal circularization times at these radii are uncertain, due to the
uncertainty of Q, they probably range from
108 - 1010 yr (Goldreich & Soter 1966), orders of magnitude greater than the millenia
it takes for our giant planets to traverse the final 0.5 AU of their
migration.
Table 4: Data describing the external surviving protoplanets at the end of giant planet migration.
In Paper I it was noted that the exterior solids disks generated by
the giant planet migration were dynamically excited,
depleted of planetesimals, and
spread over a greater radial extent than the original disk. Whilst
subsequent planet formation still seemed possible in such a disk, it
was predicted to occur over significantly extended timescales.
However, in
the Paper I scenarios only
of
planetesimals were scattered into the exterior disk, whereas in the
present model this quantity rises to
.
In addition, planetesimals are not scattered as
widely as protoplanets and their excited orbits damp rapidly when
remote from the influence of the giant planet. Whilst the decline in
gas drag with time lessens this trend in later scenarios, in all
cases the inner regions of the scattered disk remain well populated.
There is a greater similarity between the current models in the outcome of
scattering of the protoplanetary population and data for the
external protoplanets are shown in Table 4, giving their
number, mean and maximum masses and orbital inclinations, and their
mean, minimum and maximum semi major axes and eccentricities. As
expected, protoplanetary numbers fall and masses rise with disk
maturity, an effect largely due to prior accretion before the
appearance of the giant planet. Allowing for stochastic events, such
as giant impacts and strong scatterings, these data are much the
same in Paper I, as are the mean semi-major axes and eccentricities
(
). However the
minimum values of the semi-major axes and eccentricities
(
and
)
are lower than in Paper I. This
is because late-shepherded protoplanets tend to scatter rather than
assembling into hot-Neptunes and their resultant orbits damp more
rapidly as many more planetesimals are available to exert dynamical
friction.
![]() |
Figure 11: The total solids mass between 0.75-1.75 AU, before and after the giant planet migration, plotted for each scenario. |
Open with DEXTER |
What do the results of the five scenarios presented here have to say
about the probablity of forming terrestrial planets in the scattered
disks? Considering the Scenarios I-V overall, if the factor of
long-term importance is the mass distribution, rather than initial
dynamics, then such planets should form and await discovery in
hot-Jupiter systems. This point is made in Fig. 11
where the total solids mass with semi-major axes between 0.75-1.75 AU, before and after giant planet migration, is plotted for
each scenario. A clear trend is visible for less matter to be found
in this region with increasing disk maturity as it is more widely
scattered. Mass dispersal in late scenarios, however, is partially
offset by the pre-existing inward evolution of material in more
mature disks which enhances the mass present in inner regions. It is
also offset by the fact that fairly massive protoplanets are
scattered into the external disk in late Scenarios (IV and V - see
Table 4), so that terrestrial planet formation in the
scattered disk has already received a significant boost from
accretion prior to and during migration. In all cases, more than
of planet forming material remains in the
"maximum greenhouse'' habitable zone (Kasting et al. 1993) of the system
after the passage of the giant planet.
Extending our runs for the additional 50-100 Myr it would take to
form a completed external planetary system is beyond our current
computational capabilities. The principal difficulty is the presence
of the hot-Jupiter at 0.1 AU which limits
to
an excessively low value. We have, however, extended Scenario V for
2 Myr and find that (apparently stable) terrestrial planets
do form in the habitable zone of this system.
As the gas density is low when migration of the giant planet halts
in Scenario V, we assume rapid removal of the remnant gas via
photoevaporation and evolve the system in the absence of gas.
![]() |
Figure 12:
Eccentricities of bodies within 4 AU at the end of Scenario
V ( top panel) and after a further 2 Myr of gas-free accretion
( bottom panel). Protoplanets are shown as white circles and are
labeled with their mass in
![]() |
Open with DEXTER |
What will be the physical nature of such planets? The issue of disk
and protoplanet composition after migration is discussed in detail
in Sect. 3.3. Here we note, however, that the 2.15 and 2.20
planets observed in Fig. 12 are composed
of
and
of material originating from beyond the
snowline, respectively. Assuming that trans-snowline planetesimals
and protoplanets contain
water, a naive prediction is that
these planets will contain between 20-30% water by mass. The
current water inventory of the Earth is estimated to be about 0.05% by mass. Evidently terrestrial planets forming in the habitable
zones of hot-Jupiter systems are likely to host deep global oceans
- essentially being "water-worlds'' (Raymond et al. 2006; Léger et al. 2004; Kuchner 2003),
even if significant loss of volatiles occurs during high-impact
accretion.
This prediction about water content, however, depends on the giant planet forming out beyond the inner edge of the snowline, rather than at its inner edge. The closer to the inner edge the giant forms, the smaller the amount of volatile-rich material that will be shepherded inward, and the lower the degree of volatile enrichment experienced by terrestrial planets forming after migration of the giant.
![]() |
Figure 13: Composition of the original solids disk in Scenario I ( top panel) compared with the composition of the scattered disk generated through giant planet migration ( bottom panel). The key is explained in the text. |
Open with DEXTER |
A result found in all Scenarios I-V is homogenization of the solids disk composition. This mixing occurs as the giant planet drives material inward from the outer part of the disk whilst generating the external disk via random scattering. An example is given in Fig. 13 which shows, for Scenario I, the composition of the solids disk before migration (top panel) and that of the surviving external disk after migration (bottom panel). In each case the total mass of solid material is plotted in 0.5 AU width bins with the histograms labeled "dry'', "damp'' and "wet'' representing rocky material originating at <2 AU, material characteristic of chondritic meteorites between 2-2.7 AU, and trans-snowline material at >2.7 AU, respectively. It is evident that a large amount of material from beyond the snowline is shepherded into the inner regions before being left behind. Compositional mixing is similar in other scenarios, although a little less smooth in Scenarios IV & V which have had more time to accrete substantial protoplanets from local material before the appearance of the giant planet.
In contrast to other studies of water delivery to terrestrial planets (e.g. Raymond et al. 2005a,2004; Chambers 2001; Raymond et al. 2006,2005b; Morbidelli et al. 2000), we do not assign an actual water mass fraction to our compositional phases. This is because our simulations do not extend to the completion of planetary accretion and our MMSN-type original solids surface density profile, which includes a large step-increase at 2.7 AU (see Eq. (1)), is different from those adopted in the above-cited papers making a detailed comparison difficult. Typically, we find that at the end point for early-migration scenarios (I-III) originally dry protoplanets that have found their way into the scattered disk contain 0-25% of trans-snowline material, but are surrounded by a large quantity of volatile-rich planetesimals from which to accrete further. In later scenarios (IV-V), where the surviving protoplanets are more mature, this range increases to 5-70%, and whilst there is less remaining mass in small bodies to sweep up in these systems, the accretion of volatile-rich material is not yet complete.
In general we can predict that the final terrestrial planets that form in hot-Jupiter systems are likely to be much more volatile-rich than the Earth, those in the habitable zone ending up with >10% of their material originating from beyond the snowline. Accretion in the external disk, however, may pass through an early high velocity phase before completion with some protoplanets potentially losing much of their volatile inventory during giant impacts (Asphaug et al. 2006; Canup & Pierazzo 2006). Since water-rich planetesimals are abundant throughout our external disks, however, protoplanets that have lost their original volatile endowment in a catastrophic event should be able to re-acquire some water before planetary accretion is complete. Dry terrestrial planets appear improbable in hot-Jupiter systems whereas Earth-like worlds and planets with deep global oceans (Raymond et al. 2006; Léger et al. 2004; Kuchner 2003) may well be commonplace.
The models presented here and in Paper I have only explored a small region of the parameter space relevant to this problem, and have inevitably adopted assumptions that simplify or omit potentially important physical processes. We discuss some of these issues below and their possible implications for our results.
We have assumed a radius for our giant planet that is
representative of a fully contracted state, which may not be
realistic so soon after its formation.
Giant planets of approximately Jovian mass contract to radii
2-3
during their rapid gas accretion phase
(Papaloizou & Nelson 2005), before cooling and contracting toward the
Jovian radius over longer time scales.
Adoption of a larger radius would probably result in greater
accretion of solids by the giant, with gas-drag-sensitive
planetesimals being preferentially accreted.
In this paper we have presented the results of simulations that model terrestrial planet formation during and after the migration of a gas-giant planet to form a "hot-Jupiter''. This work is an extension of our previous work (Fogg & Nelson 2005), with improvements being made by modelling the viscous evolution of the gas disk, gap-formation and inner cavity formation in the gas disk due to the gravitational influence of the planet, and self-consistent type II migration. A popular belief has been that hot-Jupiter systems are unlikely to host terrestrial planets, as migration of the giant planet through the terrestrial planet zone was expected to sweep that region of planet forming material. We find, however, that the majority of this mass survives the migration episode as an interior or exterior disk remnant from which terrestrial planet formation can resume. This occurs via a combination of shepherding of the original solids disk ahead of the giant, and random scattering of the majority of this compacted material into orbits external to the giant. The net effect is not a disappearance of planetary building blocks from the inner system, but rather a stirring and mixing of material originally formed at different radial distances.
Now that more realistic gas dynamics have been included in our model, generating partial cavity formation close to the central star, gap formation in the vicinity of the giant planet's orbit, consistent type II migration rates, and a decline in the overall mass of gas present with time, we have found the following qualitative differences between the results presented here and in our previous work: