Gas accretion by planetary cores
We present accretion rates obtained from three-dimensional self-gravitating radiation hydrodynamical models of giant planet growth. We investigate the dependence of accretion rates upon grain opacity and core/protoplanet mass. The accretion rates found for low mass cores are inline with the results of previous one-dimensional models that include radiative transfer.
address=School of Physics, University of Exeter, Stocker Road, Exeter, EX4 4QL
The core accretion model of planet growth has been devised to explain the formation of giant planets such as Jupiter and Saturn [Perri and Cameron (1974), Mizuno et al. (1978)]. Through the coagulation of solid grains, planetesimals are built which then, through a so far uncertain process, bind to form large km scale bodies. When such bodies collide they shatter, but more often than not are held together by their self-gravity [Safronov and Zvjagina (1969)], allowing them to grow through collisions to reach many Earth masses. These large solid masses exert a considerable gravitational tug on the solar nebula in which they are immersed, causing them to slowly gather a gaseous atmosphere. If the solid cores are sufficiently massive (¿ 10 ) they may trigger runaway gas accretion. This occurs when the accreted gaseous atmosphere obtains a mass similar to the solid core, causing it to collapse from its hydrostatic state. The collapse leads to a rapid in fall of material from the surrounding nebula to replace the condensing gas. It is believed that such a process allows the accretion of hundreds of Earth masses of gas.
Three-dimensional hydrodynamical calculations examining the interaction of a protoplanet with a disc have usually assumed a locally isothermal equation of state, and/or do not resolve the gas flow far inside the planet’s Hill radius (Bryden et al., 1999; Kley, 1999; Lubow et al., 1999; Nelson et al., 2000; Masset, 2002; Bate et al., 2003; D’Angelo et al., 2003a, 2005; Paardekooper and Mellema, 2006; D’Angelo et al., 2006). Recent hydrodynamical calculations have also considered different equations of state and/or radiation transport (D’Angelo et al., 2003a; Morohoshi and Tanaka, 2003; Klahr and Kley, 2006; Paardekooper and Mellema, 2006, 2008). Fouchet & Mayer Fouchet and Mayer (2008) include both radiative transfer and self-gravity, but once again the resolution of gas flow near to the planetary core is limited. In this work we present results from three-dimensional self-gravitating radiative transfer hydrodynamical calculations that are able to resolve gas flow down to realistic core radii. We vary the core/protoplanet masses and the grain opacity of the nebula to examine the effects that these properties have upon the gas accretion rates of the growing planets.
2 Model setup
We model 10 - 333 solid cores/protoplanets embedded in a circumstellar disc with a surface density of 75 orbiting a 1 star at 5.2 AU. The initial surface density has a profile, whilst the initial temperature profile of the disc is , yielding a constant disc scaleheight with radius of . We vary the grain opacity of the disc from the interstellar grain opacity (IGO) down to 0.1% of this value to examine the effects this has upon the disc cooling.
To achieve high resolution near a protoplanet at reasonable computational cost we model only a small section of the protoplanetary disc centred on the planetary core. Our section size is , and radians; for details see Ayliffe and Bate (2009).
3 Results & Discussion
The accretion rates for the different mass cores/protoplanets that we modelled under differing opacities are shown in the left panel of figure 1. The trend in accretion rate against core/protoplanet mass is as would be expected, increasing as the central bodies mass increases. There is a reduction in accretion rates for the highest mass protoplanets caused by the evacuation of a large gap in the encompassing circumstellar disc by these protoplanets. Such a gap reduces the material available in the feeding zone of the protoplanet. The accretion rates also show an increasing trend with reduced grain opacities for the low mass cores. Reduced grain opacities allow the infalling material to cool more readily via radiation, and so the thermal pressure supporting it against capture is reduced. The only exception to this is at the lowest grain opacity of 0.1% IGO for the lowest mass core where the radiation and gas fields interact so little that the two fields decouple, preventing effective gas cooling via radiation. For the high mass protoplanets the grain opacity has very little impact on the accretion rates as the process is gravitationally dominated.
The right hand panel of figure 1 compares the accretion on to a 10 solid core in our models with the accretion onto 5 and 15 cores in the one-dimensional models of Papaloizou & Nelson Papaloizou and Nelson (2005). It can be seen that the results show good agreement, suggesting that for at least the phase of accretion preceding runaway growth, the transition from one-dimensional models to three-dimensional models does not have a significant impact on the accretion rates.
Figure 2 illustrates the differences between radiation hydrodynamical and isothermal calculations. Features in the vicinity of the planet that are well established under isothermal conditions, such as spiral shocks and circumplanetary discs, are weakened in the radiative transfer calculations. The thermal pressure acts to smear out these features. It can also be seen in figure 2 that reducing the grain opacity leads to a model that is more similar to the isothermal case. The increased cooling seen at lower opacities reduces the thermal pressure, thus mitigating its impact on the shocks and circumplanetary discs.
- F. Perri, and A. G. W. Cameron, Icarus 22, 416–425 (1974).
- H. Mizuno, K. Nakazawa, and C. Hayashi, Progress of Theoretical Physics 60, 699–710 (1978).
- V. S. Safronov, and E. V. Zvjagina, Icarus 10, 109–+ (1969).
- G. Bryden, X. Chen, D. N. C. Lin, R. P. Nelson, and J. C. B. Papaloizou, ApJ 514, 344–367 (1999).
- W. Kley, MNRAS 303, 696–710 (1999).
- S. H. Lubow, M. Seibert, and P. Artymowicz, ApJ 526, 1001–1012 (1999).
- R. P. Nelson, J. C. B. Papaloizou, F. Masset, and W. Kley, MNRAS 318, 18–36 (2000).
- F. S. Masset, A&A 387, 605–623 (2002).
- M. R. Bate, S. H. Lubow, G. I. Ogilvie, and K. A. Miller, MNRAS 341, 213–229 (2003).
- G. D’Angelo, T. Henning, and W. Kley, ApJ 599, 548–576 (2003a).
- G. D’Angelo, M. R. Bate, and S. H. Lubow, MNRAS 358, 316–332 (2005).
- S.-J. Paardekooper, and G. Mellema, A&A 459, L17–L20 (2006).
- G. D’Angelo, S. H. Lubow, and M. R. Bate, ApJ 652, 1698–1714 (2006).
- K. Morohoshi, and H. Tanaka, MNRAS 346, 915–923 (2003).
- H. Klahr, and W. Kley, A&A 445, 747–758 (2006).
- S.-J. Paardekooper, and G. Mellema, A&A 478, 245–266 (2008).
- L. Fouchet, and L. Mayer, ArXiv e-prints 806 (2008).
- B. A. Ayliffe, and M. R. Bate, MNRAS 393, 49–64 (2009).
- J. C. B. Papaloizou, and R. P. Nelson, A&A 433, 247–265 (2005).
- G. D’Angelo, W. Kley, and T. Henning, ApJ 586, 540–561 (2003b).