Planet Formation

1. Planet Formation Topic: Orbital dynamics and the restricted 3-body problem Lecture by: C.P. Dullemond

2. Objects of the Solar System Gas Giants Ice Giants Rocky Planets Earth Mercury Venus Mars Jupiter Saturn Uranus Neptune Planets 0.30- 0.47 AU e = 0.2 1.38- 1.67 AU e = 0.1 1 AU 0.72 AU 30 AU 9.05-10.1 AU e = 0.05 18-20 AU e = 0.05 4.95-5.46 AU e = 0.05 Ceres Pluto Dwarf- Planets + Haumea, Makemake, Eris 2.55- 2.99 AU e = 0.08 29.6- 48.9 AU e = 0.25 Asteroid Belt (2.0-3.4 AU) Kuiper Belt (30-50 AU) Asteroids & TNOs 1 AU 10 AU Logarithmic distance scale

3. Kepler orbits, eccentricity Kreisbahn (Mplanet << M*): Planet (Kepler Frequency) Sonne a (Kepler Velocity) („True anomaly“) y x

4. Kepler orbits, eccentricity Elliptic orbit(Mplanet << M*): Planet (Kepler Freq.) r Sonne Apoapsis (Aphelion) Periapsis (Perihelion) a e a Fokus a („Semi-major axis“) e („eccentricity“) („Meananomaly“) y („Eccentricanomaly“ , solvenumerically) x („True anomaly“)

5. Kepler orbits, eccentricity Elliptic orbit(Mplanet << M*): Planet (Kepler Freq.) r Sonne Apoapsis (Aphelion) Periapsis (Perihelion) a e a Fokus Total energy: y Note: Total energy depends only on a! x Angular momentum: Note: eccentricity thus leads to an angular momentum deficit.

6. Orbital elements Compared to the equatorial plane of the solar system (or exoplanetary system), and compared to a reference direction, you can uniquely orient the ellipse. Semi-major axis a, eccentricity e and true anomalyν, together with orientation (i, ω, Ω) are: orbital elements. From: Wikipedia: http://en.wikipedia.org/wiki/File:Orbit1.svg

7. Guiding-center and epicyclic motion If an orbit is almost circular, then we can describe this elliptic orbit as a circular orbit with epicyclic motion superposed on it. e=0.3 e=0.3

8. Guiding-center and epicyclic motion When using another orbital frequency as the reference frame than the orbital frequency of the particle, this epicyclic motion looks like: e=0.3 e=0.3

9. Guiding-center and epicyclic motion For small eccentricity the epicycle becomes an ellipse. As we will see (and use) later: these small deviations from Kepler can be described as velocity disturbances Δv. guiding center epicycle

10. Even kids are familiar with this... ;-)

11. What is the „restricted 3-body problem“? • Three bodies, one of which (M3) is a „test particle“: M3 <<< M2 < M1 • Bodies 1 and 2 only feel each other‘s gravity and thus perfectly follow a Kepler orbit. • Assume that body 1 and 2 are in perfect circular orbits • Put the center of coordinate system at center of mass • Body 3 feels bodies 1 and 2. • Resulting motions: • Some orbits are stable • Some orbits are unstable (body 3 gets ejected) • Some orbits are chaotic: Chaos theory! • Chaotic orbits are unpredictable on the long run.

12. Equations of motion for test particle Remember from chapter „Turbulence“ Section „Magnetorotational instability“ the equation of motion of a test particle in a rotating frame: (from chapter „Turbulence“) Now we do the same (though now we put x=0 at the center of mass of the entire system), we drop the fx and fy forces but now we include the forces of both the star and the planet. With the „effective potential“ given by: Coriolis forces Gravity and centrifugal forces Exercise: re-derive these equations.

13. L4 L3 L1 L2 L5 centrifugal kinetic energy Example: M2/M1=0.1 Effective potential, Lagrange points r1 r2 Effective potential in the co-rotating frame: Gravitational potential

14. L4 Effective potential in the co-rotating frame: L5 centrifugal kinetic energy Effective potential, Lagrange points r1 r2 L3 L1 L2 Gravitational potential Example: M2/M1=0.01

15. Jacobi‘s Integral The full 3-D set of equations is: now multiply by: and add them all up:

16. Jacobi‘s Integral This can be integrated once, to obtain: Traditionally the constant C is written as -½CJ: CJ is called Jacobi‘s constant or Jacobi‘s Integral of motion. For the restricted 3-body problem it is the only integral of motion, i.e. there exist no closed-form solutions.

17. Zero-velocity curves / surfaces >0 <0 Jacobi‘s constant is some kind of energy, sometimes called „Integral of relative energy“. It is the rotational equivalent of minus twice the total (potential + kinetic) energy of a test particle in a non- rotating system: Since the kinetic term <0, we know that a given particle on a given orbit (with a given constant CJ), can only reach points (x,y,z) where the effective potential obeys: Remember:

18. Zero-velocity curves / surfaces (allowed region in x,y,z) The boundaries of this region are called the „zero velocity curves“ (in 2-D) or „zero velocity surfaces“ (in 3-D). They are the potential lines. For CJ<0 no such restrictions exist (all points are allowed). More precisely: for CJ<min(-2Φeff) no such restriction exist. But for CJ>min(-2Φeff) there exists regions in (x,y,z) which are inaccessible for the particle. If CJ is sufficiently large, these inaccessible regions can even completely surround the star or planet system or both: = Not allowed region

19. Hill sphere (=Roche lobe) Not exactly a sphere, but approximately. It is the largest zero-velocity surface surrounding only the planet. It is the sphere of influence of the planet.

20. Meaning of Hill sphere • The Hill sphere plays a key role in planet formation: • To add mass to a planet, we must put the mass into the Hill sphere of the growing planet, because only then it that mass gravitationally bound to it. • Any circumplanetary disk or moons must be inside the planet‘s Hill radius • We will see later that the ratio of the Hill radius to the protoplanetary disk‘s thickness plays a key role in planet migration. • Any object that is larger than its own Hill sphere will be sheared apart by the tidal forces of the star. This will lead us to the definition of the „Roche density“ as the minimal density an object needs to remain gravitationally coherent and survive tidal forces. See next chapter.

21. Potential field lines as approximate orbits If then one can approximately write: This means that approximately the test particle moves along the potential field lines. To be more precise: the guiding center will do this; the test particle will epicycle around this guiding center. You can find orbits without epicycles, in which case the test particle indeed moves approximately along the field lines. It turns out to be also a fairly good approximation if the condition does not hold. This leads to a special set of orbits.

22. Kepler- and horseshoe orbits Kepler orbit around star only Kepler orbit around planet only Kepler orbit around both Horseshoe orbit

23. Trojans of Jupiter L4 and L5 are stable Lagrange points, while L1, L2 and L3 are not.