CO snow lines are stabilised by the vertical transport of volatiles

Abstract Volatile evolution in protoplanetary discs determines the compositional evolution of forming planets. Below their sublimation temperatures, volatiles freeze out from the vapour phase onto dust grains in the disc and transition to being dynamically-coupled to the dust component as opposed to the gas. The boundary between the ice and vapour phases is referred to as the snow line, when thought of as the mid-plane radius at which the phase transition occurs, or the snow surface, when viewed as a 2D (radial and vertical) structure in the disc. We investigate whether the CO snow line (and therefore snow surface) is thermally unstable and therefore liable to changes in its location during disc evolution using the disc evolution code cuDisc, to which we have added an ice-vapour chemistry solver. We find that the instability does lead to there being two steady-state stable equilibrium solutions for the snow surface when including the vertical structure. However, in dynamically-evolving simulations, the disc does not enter a limit-cycle — as seen in previous 1D models — due to the shape of the 2D snow surface and the vertical transport of volatiles. We therefore expect that dynamically evolution of snow lines due to instability is limited to transient, stochastic events rather than oscillatory behaviour with a regular period. However, we also expect the snow surface to evolve substantially during the disc lifetime solely due to changes in the thermal structure driven by evolution of the dust spatial structure and grain-size distribution — this we will explore in future models.