Evolution of migrating protoplanets heated by pebble accretion

We study the interactions in a protoplanetary system consisting of a gas disk, a pebble disk and embedded low-mass protoplanets. The hydrodynamic simulations are performed using a new code based on 2D FARGO (Masset 2000) which we call FARGO_THORIN (http://sirrah.troja.mff.cuni.cz/~chrenko/). The code treats the hydrodynamics of gas and pebbles within a two-fluid approximation, accounts for the heating and cooling processes in the gaseous component (including heating due to pebble accretion) and propagates the planets in 3D using a high-order integration scheme (IAS15; Rein & Spiegel 2015). Our aim is to investigate how pebble accretion alters the orbital evolution of protoplanets undergoing Type-I migration.First, we demonstrate that pebble accretion can heat the protoplanets so that their luminosity induces the heating torque (Benítez-Llambay et al. 2015) and the hot-trail effect (Chrenko et al. 2017; Eklund & Masset 2017). The heating torque is always positive and alters the migration rates and directions profoundly, thus changing the position of planet traps and deserts. The hot-trail effect, on the other hand, pumps the eccentricity of initially circular orbits up to e ~ h. After becoming eccentric, the protoplanets exhibit reduced probability of resonant locking during the migration and moreover, their close encounters become more frequent and provide more opportunities for scattering or merger events. The mergers can be massive enough to become giant planet cores. We discuss the importance of the excited eccentricities and violent orbital evolution for the extrasolar planet population synthesis. Finally, we present an extended model with flux-mean opacities caused by a coupled disk of coagulating dust grains with a realistic size distribution. The aim of this model is to constrain possible pathways of migrating planets towards the inner rim of the protoplanetary disk.