TEMPERATURE STRUCTURE OF PROTOPLANETARY DISKS UNDERGOING LAYERED ACCRETION

Abstract

We calculate the temperature structures of protoplanetary disks (PPDs) around T Tauri stars heated by both incident starlight and viscous dissipation. We present a new algorithm for calculating the temperatures in disks in hydrostatic and radiative equilibrium, based on Rybicki's method for iteratively calculating the vertical temperature structure within an annulus. At each iteration, the method solves for the temperature at all locations simultaneously, and converges rapidly even at high (>>10{sup 4}) optical depth. The method retains the full frequency dependence of the radiation field. We use this algorithm to study for the first time disks evolving via the magnetorotational instability. Because PPD midplanes are weakly ionized, this instability operates preferentially in their surface layers, and disks will undergo layered accretion. We find that the midplane temperatures T{sub mid} are strongly affected by the column density {Sigma}{sub a} of the active layers, even for fixed mass accretion rate M-dot . Models assuming uniform accretion predict midplane temperatures in the terrestrial planet forming region several x 10{sup 2} K higher than our layered accretion models do. For M-dot < 10{sup -7} M{sub sun} yr{sup -1} and the column densities {Sigma}{sub a} < 10 g cm{sup -2} associated with layered accretion, diskmore » temperatures are indistinguishable from those of a passively heated disk. We find emergent spectra are insensitive to {Sigma}{sub a}, making it difficult to observationally identify disks undergoing layered versus uniform accretion.« less