Structure of Passive Circumstellar Disks: Beyond the Two-Temperature Approximation

Roman R. Rafikov,Fabio De Colle

DOI: https://doi.org/10.1086/504833

2006-01-31

Abstract:Structure and spectral energy distributions (SEDs) of externally irradiated circumstellar disks are often computed on the basis of the two-temperature model of Chiang & Goldreich. We refine these calculations by using a more realistic temperature profile which is continuous at all optical depths and thus goes beyond the two-temperature model. It is based on the approximate solution of the radiation transfer in the disk obtained from the frequency-integrated moment equations in the Eddington approximation. We come up with a simple procedure (``constant g_z approximation'') for treating the vertical structure of the disk in regions where its optical depth to stellar radiation is high. This allows us to obtain expressions for the vertical profiles of density and pressure at every point in the disk and to determine the shape of its surface. Armed with these analytical results we calculate the full radial structure of the disk and demonstrate that it favorably agrees with the results of direct numerical calculations. We also describe a simple and efficient way of the SED calculation based on our adopted temperature profile. Resulting spectra provide very good match (especially at short wavelengths) to the results of more detailed (but also more time-consuming) SED calculations solving the full frequency- and angle-dependent radiation transfer within the disk.

Astrophysics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to improve the computational methods for the structure and spectral energy distribution (SEDs) of circumstellar disks around stars, going beyond the traditional two - temperature model. Specifically:

1. **Improving the temperature profile**: The traditional method uses a two - temperature model to describe the temperature distribution in the disk, namely the outer "super - hot" layer and the isothermal mid - plane region. This method has a discontinuous temperature jump at an optical depth of 1, which is not in line with the actual situation. This paper proposes a more realistic temperature profile that is continuous at all optical depths, based on an approximate solution of radiative transfer.

2. **Handling the vertical structure**: In order to handle the vertical structure of the disk, especially in the high - optical - depth region, the author introduces the "constant gz approximation" method. This method allows obtaining the vertical profiles of density and pressure at any point in the disk and determining the shape of the disk surface.

3. **Calculating the radial structure**: Using the above - mentioned analysis results, the author calculates the complete radial structure of the disk and shows its good agreement with the results of direct numerical calculations.

4. **Calculating the spectral energy distribution (SED)**: Based on the improved temperature profile, the author provides a simple and efficient way to calculate the SED, especially in the short - wavelength range, which is very consistent with the results of detailed radiative transfer calculations.

### Formula summary

- **Optical depth**:
\[
\tau_\star(z) \approx \alpha^{-1} \kappa_P(T_\star) \Sigma(z)
\]
where \(\rho(z)\) is the local gas density, \(\Sigma(z)=\int_z^\infty \rho(z') dz'\) is the disk surface density above a given height \(z\), \(\kappa_P(T)\) is the Planck - mean extinction coefficient, and \(\alpha\) is the expansion angle.

- **Temperature distribution**:
\[
T^4(z)=\phi T_\star^4\left(\frac{R_\star}{a}\right)^2\left[\frac{\alpha}{2\psi_{ex}} \psi_{sh}+\frac{q(T)}{4} e^{-\tau_\star(z)}\right]
\]
where \(q(T)=\frac{\kappa_P(T_\star)}{\kappa_P(T)}\) is the ratio of the extinction coefficients at temperatures \(T\) and \(T_\star\), \(\phi\) is the integral of the stellar surface area visible from the disk surface, and \(\tau_\star(z)\) and \(\alpha\) are given by equations (2) and (3) respectively.

- **Exposed - layer temperature**:
\[
T\approx T_\star\left(\frac{R_\star}{2a}\right)^{1/2}\left[\phi q(T) e^{-\tau_\star(z)}\right]^{1/4}
\]

- **Shielded - layer temperature**:
\[
T_{sh}(a)=T_\star\left(\frac{\phi \alpha}{2 \psi_{ex} \psi_{sh}}\right)^{1/4}\left(\frac{R_\star}{a}\right)^{1/2}
\]

Through these improvements, the author is able to describe the structure and spectral characteristics of the disk more accurately, especially the spectral energy distribution in the short - wavelength range.