Using deep neural networks to compute the mass of forming planets
Abstract
Context. Computing the mass of planetary envelopes and the critical mass beyond which planets accrete gas in a runaway fashion is important for studying planet formation, in particular, for planets up to the Neptunemass range. This computation in principle requires solving a set of differential equations, the internal structure equations, for some boundary conditions (pressure, temperature in the protoplanetary disc where a planet forms, core mass, and the rate of accretion of solids by the planet). Solving these equations in turn proves to be timeconsuming and sometimes numerically unstable.
Aims: The aim is to provide a way to approximate the result of integrating the internal structure equations for a variety of boundary conditions.
Methods: We computed a set of internal planetary structures for a very large number (millions) of boundary conditions, considering two opacities: that of the interstellar medium, and a reduced opacity. This database was then used to train deep neural networks (DNN) in order to predict the critical core mass and the mass of planetary envelopes as a function of the boundary conditions.
Results: We show that our neural networks provide a very good approximation (at the percent level) of the result obtained by solving interior structure equations, but the required computer time is much shorter. The difference with the real solution is much smaller than the difference that is obtained with some analytical formulas that are available in the literature, which only provide the correct order of magnitude at best. We compare the results of the DNN with other popular machinelearning methods (random forest, gradient boost, support vector regression) and show that the DNN outperforms these methods by a factor of at least two.
Conclusions: We show that some analytical formulas that can be found in various papers can severely overestimate the mass of planets and therefore predict the formation of planets in the Jupitermass regime instead of the Neptunemass regime. The python tools that we provide allow computing the critical mass and the mass of planetary envelopes in a variety of cases, without the requirement of solving the internal structure equations. These tools can easily replace previous analytical formulas and provide far more accurate results.
A copy of the Python code is available at the CDS via anonymous ftp to http://cdsarc.ustrasbg.fr (ftp://130.79.128.5) or via http://cdsarc.ustrasbg.fr/vizbin/qcat?J/A+A/626/A21
 Publication:

Astronomy and Astrophysics
 Pub Date:
 June 2019
 DOI:
 10.1051/00046361/201834942
 arXiv:
 arXiv:1903.00320
 Bibcode:
 2019A&A...626A..21A
 Keywords:

 methods: numerical;
 planets and satellites: formation;
 Astrophysics  Earth and Planetary Astrophysics
 EPrint:
 accepted in A&