Time-dependent, long-term hydrodynamic simulations of the inner protoplanetary disk

Abstract

Context. The spin evolution of young protostars, surrounded by an accretion disk, still poses problems for observations and theoretical models. In recent studies, the importance of the magnetic star-disk interaction for stellar spin evolution has been elaborated. The accretion disk in these studies, however, is only represented by a simplified model and important features are not considered. Aims. A more realistic representation of the accretion disk is indispensable for a better understanding of the star-disk interaction and the stellar spin evolution. The aim of this study is to investigate the influence of a hydrodynamic disk evolution on the stellar rotational period and vice versa during the accretion phase. Methods. We combined the implicit hydrodynamic TAPIR disk code with a stellar spin evolution model. The influence of stellar magnetic fields on the disk dynamics, the radial position of the inner disk radius, as well as the influence of stellar rotation on the disk were calculated self-consistently. Results. Within a defined parameter space, we can reproduce the majority of fast and slow rotating stars observed in young stellar clusters. Additionally, the back reaction of different stellar spin evolutionary tracks on the disk can be analyzed. Disks around fast rotating stars are located closer to the star. Consequently, the disk midplane temperature in the innermost disk region increases significantly compared to slow rotating stars. We can show the effects of stellar rotation on episodic accretion outbursts. The higher temperatures of disks around fast rotating stars result in more outbursts and a longer outbursting period over the disk lifetime. Conclusions. The combination of a long-term hydrodynamic disk and a stellar spin evolution model allows the inclusion of previously unconsidered effects such as the back-reaction of stellar rotation on the long-term disk evolution and the occurrence of accretion outbursts. However, a wider parameter range has to be studied to further investigate these effects. Additionally, a possible interaction between our model and a more realistic stellar evolution code (e.g., the MESA code) can improve our understanding of the stellar spin evolution and its effects on the pre-main sequence star.