Nonlinear Outcome of Gravitational Instability in Disks with Realistic Cooling

Abstract

We consider the nonlinear outcome of gravitational instability in optically thick disks with a realistic cooling function. We use a numerical model that is local, razor-thin, and unmagnetized. External illumination is ignored. Cooling is calculated from a one-zone model using analytic fits to low temperature Rosseland mean opacities. The model has two parameters: the initial surface density Sigma0 and the rotation frequency Omega. We survey the parameter space and find: (1) The disk fragments when tc,eff Omega = 1, where tc,eff is an effective cooling time defined as the average internal energy of the model divided by the average cooling rate. This is consistent with earlier results that used a simplified cooling function. (2) The initial cooling time tc0 or a uniform disk with Q = 1 can differ by orders of magnitude from tc,eff in the nonlinear outcome. The difference is caused by sharp variations in the opacity with temperature. The condition tc0 Omega = 1 therefore does not necessarily indicate where fragmentation will occur. (3) The largest difference between tc,eff and tc0 is near the opacity gap, where dust is absent and hydrogen is largely molecular. (4) In the limit of strong illumination the disk is isothermal; we find that an isothermal version of our model fragments for Q < 1.4. Finally, we discuss some physical processes not included in our model, and find that most are likely to make disks more susceptible to fragmentation. We conclude that disks with tc,eff Omega < 1 do not exist.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…