Nonlinear outcome of gravitational instability in disks with realistic cooling

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 Sigma(0) and the rotation frequency Omega. We survey the parameter space and find the following. (1) The disk fragments when [[tau(c)]] Omega similar to 1, where [[tau(c)]] 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 tau(c0) for a uniform disk with Q = 1 can differ by orders of magnitude from [[tau(c)]] in the nonlinear outcome. The difference is caused by sharp variations in the opacity with temperature. The condition tau(c0)Omegasimilar to 1 therefore does not necessarily indicate where fragmentation will occur. (3) The largest difference between [[tau(c)]] and tau(c0) 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. nd that an isothermal version of our model fragments for Q less than or similar to 1.4. Finally, we discuss some physical processes not included in our model and. nd that most are likely to make disks more susceptible to fragmentation. We conclude that disks with [[tau(c)]] Omega less than or similar to 1 do not exist.