Abstract
The thermal structure of a protoplanetary disc is regulated by the opacity that dust grains provide. However, previous works have often considered simplified prescriptions for the dust opacity in hydrodynamical disc simulations, for example, by considering only a single particle size. In the present work, we perform 2D hydrodynamical simulations of protoplanetary discs where the opacity is self-consistently calculated for the dust population, taking into account the particle size, composition, and abundance. We first compared simulations utilizing single grain sizes to two different multi-grain size distributions at different levels of turbulence strengths, parameterized through the α-viscosity, and different gas surface densities. Assuming a single dust size leads to inaccurate calculations of the thermal structure of discs, because the grain size dominating the opacity increases with orbital radius. Overall the two grain size distributions, one limited by fragmentation only and the other determined from a more complete fragmentation-coagulation equilibrium, give comparable results for the thermal structure. We find that both grain size distributions give less steep opacity gradients that result in less steep aspect ratio gradients, in comparison to discs with only micrometer-sized dust. Moreover, in the discs with a grain size distribution, the innermost (<5 AU) outward migration region is removed and planets embedded in such discs experience lower migration rates. We also investigated the dependency of the water iceline position on the alpha-viscosity (α), the initial gas surface density (Σg,0) at 1 AU and the dust-to-gas ratio (fDG) and find rice ∝ α0.61Σg,00.8fDG0.37 independently of the distribution used in the disc. The inclusion of the feedback loop between grain growth, opacities, and disc thermodynamics allows for more self-consistent simulations of accretion discs and planet formation.
Highlights
Protoplanetary discs surround young stars for the first few million years after their formation and they are the birthplaces of planetary systems
The structure of this paper is as follows: in Sect. 2 we describe the energy equations used in the hydrodynamical simulations, the new opacity module and the two grain size distributions that we included in the code
All of these suggested effects have been taken into account in the here presented work where we study the influence of α-viscosity, initial gas surface density, and total dust-to-gas ratio on the position of the iceline
Summary
Protoplanetary discs surround young stars for the first few million years after their formation and they are the birthplaces of planetary systems. It was shown by Tanaka et al (1996) that the ξ constant is independent of the specific parameters of the collisional outcome model, as long as it is self-similar, which in this case means that the outcome of impacts between dust grains depends on the masses of two colliding particles only through their ratio Such a description of a grain size distribution with only one power law is a simplification, since it only takes into account the coagulation/fragmentation equilibrium. Several works in the recent years aimed to couple the dust and gas components of protoplanetary discs in simulations and in most of the cases such models include a grain size distribution. We consider how grain dynamics and how grain size distributions affect the opacity and as a consequence the thermal structure of the disc in order to simulate the whole feedback loop.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.