First we devised a method for calculating the inner disk radius. The traditional approach was to use the dust sublimation temperature Tsub as a unique well-defined boundary condition, based on the assumption that the temperature drops with the distance from the star. However, it has been shown numerically (Kama et al. 2009; Vinković 2012) and theoretically (Vinković 2006) that the temperature of the big grains that populate the surface of the inner disk actually increases with distance within the optically thin dusty zone (see sketch in Fig. 1) before it starts to monotonically decrease within the optically thick dust. This temperature inversion effect results in a local temperature maximum within the dust cloud (see Fig. 3 in Kama et al. 2009), where dust can overheat. A precise calculation of the dust distribution that includes this effect requires information on the spatial dust distribution within the optically thin zone. Since the density structure of the optically thin dust zone cannot be geometrically constrained, we are left with a highly ambiguous concept of the inner disk rim.
Vinković (2006) proposed using two inner disk radii that distinguish between two coexisting – optically thin and thick – disk zones around a star of luminosity L∗: (A.1)where Rin is the disk radius in the midplane, Lacc is the accretion luminosity described below and Ψ = 2 and Ψ ~ 1.2 are used for the optically thick and thin inner radii, respectively. Here we considered only optically thick disks with Tsub as the dust sublimation temperature of the inner disk rim surface. The sublimation temperature depends on the gas density ρgas. Based on the fit to the data from Pollack et al. (1994), we used: (A.2)This is sufficient for the range of gas densities used in our model. A more detailed description of Tsub is given by Kama et al. (2009).
The surface of the inner disk is directly exposed to stellar heating, which creates a vertical temperature structure very close to isothermal in the narrow spatial region of our interest. The high-resolution radiative transfer calculation by Vinković (2012) showed that this is a good approximation for the inner rim of the disk. This approximation can only overestimate the dusty disk height because in reality the temperature slightly drops with the height above midplane. The radial temperature change is approximated with (A.3)The gas is in vertical hydrostatic equilibrium when the gas density structure is described as
where ρ0 = ρgas(Rin,0) is the midplane gas density at the inner disk rim, Hp is the disk scale height, cs is the isothermal sound speed, ΩK is the Keplerian angular velocity at the midplane, M∗ is the stellar mass, μgas = 2.33 is the mean molecular weight in units of the proton mass mp, and k is Boltzmann’s constant.
Takeuchi & Lin (2002) showed that the vertical component of the gas velocity is much smaller than the horizontal component vr,gas, given by (A.8)Accretion of the gas onto the star contributes with its luminosity to the total energy that heats the inner rim. To accommodate for this effect we used the accretion luminosity (Calvet & Gullbring 1998; Muzerolle et al. 2004) in Eq. (A.1) (A.9)where M∗ is the stellar mass, ξ ~ 0.8 is the correction factor for the magnetospheric truncation of the disk, G is the gravitational constant, and R∗ is the stellar radius.
The accretion rate Ṁacc can be calculated from the vertical profile of the gas density (Eq. (A.4)) and the horizontal velocity (Eq. (A.8)) (A.10)which exists as inflow only if (A.11)The accretion rate depends on the distance from the star, which drives the disk evolution, except for q = −0.5 when the radial velocity is independent of r and the disk is in steady-state. We focused our attention on a very small part of the disk around Rin, where Ṁacc varies slowly. Hence, we define the accretion rate as Ṁacc(Rin). Equation (A.10) also shows that Ṁacc ∝ ρ0. We used this to derive ρ0 from Ṁacc, which is in turn derived from Lacc ∝ Ṁacc in Eq. (A.9).
The local viscous dissipation of energy in the accretion disks contributes to the local disk temperature, and we checked for this effect. The net temperature at Rin is the combination (Ciesla & Cuzzi 2006) of the stellar irradiation (Eq. (A.1), with Tsub replaced by Tirr) and viscous heating (A.12)where σ is the Stefan-Boltzmann constant. If , we used Eq. (A.1) to derive Rin and if , we used Eq. (A.12). These conditions have to be checked iteratively when calculating Rin.
The long-term stability of the inner disk rim depends on the radial migration of the dust particles. We followed the prescription by Takeuchi & Lin (2002) where dust particles are driven by gas drag and their radial velocity depends on the distance from the midplane. The source of this vertical velocity variation is the gas pressure gradient, which causes two effects: the gas rotation is different from Keplerian rotation above the midplane, and the gas drag force depends on gas density. Without the gas drag the dust particles would follow Keplerian orbits. We considered particles smaller than the mean free path of gas molecules. This condition is described with Epstein’s gas drag law that gives the gas drag stopping time:
We considered turbulent disks where the gas turbulent motion stirs up dust particles to higher altitudes to prevent dust sedimentation. The equilibrium distribution of dust particles is reached when the sedimentation due to gas drag and the diffusion due to turbulence are balanced. Using the α-prescription for the viscous effect of turbulence (Shakura & Sunyaev 1973; Pringle 1981), the dust density distribution is (Takeuchi & Lin 2002) (A.15)where ρd(r) is the midplane dust density profile. The Schmidt number Sc is a measure of the coupling between the particles and the gas. Sc ~ 1 for the small particles used in our paper and increases to infinity for very big particles. This means that very big particles are not influenced by the gas turbulence and they accumulate in the midplane. The expression of Sc is a matter of debate, and we adopted the expression suggested by Youdin & Lithwick (2007) for diffusion driven by anisotropic turbulence (A.16)Takeuchi & Lin (2003) showed that the stellar radiation pressure force slowly erodes the dusty disk surface, but they did not investigate the impact of this force on the surface of inner disk rim. We used their methodology to include the radiation pressure force in the dust velocity profiles. We ignored the complicated radiation pressure force from the near-infrared photons from the disk interior (Vinković 2009), but this force can only enhance the disk erosion and reduce the height of the optically thick disk. The stellar radiation pressure force compares to the stellar gravity by a factor of (Vinković 2009) (A.17)where ⟨ Qext ⟩ is the dust extinction coefficient averaged over the stellar spectrum. We used a stellar black-body spectrum to calculate ⟨ Qext ⟩ from olivine dust grains. However, the result is close
to ⟨ Qext ⟩ ~ 2 because the stellar radiation peaks at wavelengths shorter than the size of big grains that populate the inner disk surface. The radial component of the dust velocity with included radiation pressure was derived by Takeuchi & Lin (2003) as
Our goal is to find the relative height zmax/Rin that separates the inflow below this height from the outflow vr,dust(Rin,zmax) > 0above zmax. We varied the parameters listed in Table 1 and iteratively searched for ρ0 consistent with Ṁacc, which in turn has to be consistent with Lacc. The pseudocode describing the computational procedure is the following:
© ESO, 2014