On the evolution of vortices in massive protoplanetary discs

Abstract

It is expected that a pressure bump can be formed at the inner edge of a dead zone, and where vortices can develop through the Rossby Wave Instability (RWI). It has been suggested that self-gravity can significantly affect the evolution of such vortices. We present the results of 2D hydrodynamical simulations of the evolution of vortices forming at a pressure bump in self-gravitating discs with Toomre parameter in the range 4-30. We consider isothermal plus non-isothermal disc models that employ either the classical β prescription or a more realistic treatment for cooling. The main aim is to investigate whether the condensating effect of self-gravity can stabilize vortices in sufficiently massive discs. We confirm that in isothermal disc models with Q ≳ 15, vortex decay occurs due to the vortex self-gravitational torque. For discs with 3 ≲ Q ≲ 7, the vortex develops gravitational instabilities within its core and undergoes gravitational collapse, whereas more massive discs give rise to the formation of global eccentric modes. In non-isothermal discs with β cooling, the vortex maintains a turbulent core prior to undergoing gravitational collapse for β ≲ 0.1, whereas it decays if β ≥ 1. In models that incorporate both self-gravity and a better treatment for cooling, however, a stable vortex is formed with aspect ratio X ~ 3-4. Our results indicate that self-gravity significantly impacts the evolution of vortices forming in protoplanetary discs, although the thermodynamical structure of the vortex is equally important for determining its long-term dynamics.