Photoevaporation by stellar ionizing radiation is believed to play an important role in the dispersal of disks around young stars. The mass loss model for dust-free disks developed by Hollenbach et al. is currently regarded as a conventional one and
has been used in a wide variety of studies. However, the rate in this model was derived by the crude so-called 1+1D approximation of ionizing radiation transfer, which assumes that diffuse radiation propagates in a direction vertical to the disk. In this study, we revisit the photoevaporation of dust-free disks by solving the 2D axisymmetric radiative transfer for steady-state disks. Unlike that solved by the conventional model, we determine that direct stellar radiation is more important than the diffuse field at the disk surface. The radial density distribution at the ionization boundary is represented by the single power-law with an index -3/2 in contrast to the conventional double power-law. For this distribution, the photoevaporation rate from the entire disk can be written as a function of the ionizing photon emissivity, Phi_EUV, from the central star and the disk outer radius, r_d, as follows: Mdot_PE = 5.4 x 10^-5 x (Phi_EUV/10^49 sec^-1)^1/2 x (r_d/1000 AU)^1/2 Msun/yr. This new rate depends on the outer disk radius rather than on the gravitational radius as in the conventional model, caused by the enhanced contribution to the mass loss from the outer disk annuli. In addition, we discuss its applications to present-day as well as primordial star formation.
The temperature in the optically thick interior of protoplanetary discs is essential for the interpretation of millimeter observations of the discs, for the vertical structure of the discs, for models of the disc evolution and the planet formation, a
nd for the chemistry in the discs. Since large icy grains have a large albedo even in the infrared, the effect of scattering of the diffuse radiation in the discs on the interior temperature should be examined. We have performed a series of numerical radiation transfer simulations including isotropic scattering by grains with various typical sizes for the diffuse radiation as well as for the incident stellar radiation. We also have developed an analytic model including isotropic scattering to understand the physics concealed in the numerical results. With the analytic model, we have shown that the standard two-layer approach is valid only for grey opacity (i.e. grain size $ga10$ micron) even without scattering. A three-layer interpretation is required for grain size $la10$ micron. When the grain size is 0.1--10 micron, the numerical simulations show that isotropic scattering reduces the temperature of the disc interior. This reduction is nicely explained by the analytic three-layer model as a result of the energy loss by scatterings of the incident stellar radiation and of the warm diffuse radiation in the disc atmosphere. For grain size $ga10$ micron (i.e. grey scattering), the numerical simulations show that isotropic scattering does not affect the interior temperature. This is nicely explained by the analytic two-layer model; the energy loss by scattering in the disc atmosphere is exactly offset by the green-house effect due to scattering of the cold diffuse radiation in the interior.
