A. Isella ^{1,2} - A. Natta ^{1}
1 - Osservatorio Astrofisico di Arcetri, INAF, Largo E. Fermi 5,
50125 Firenze, Italy
2 -
Dipartimento di Fisica, Universitá di Milano, via Celoria 16,
20133 Milano, Italy
Received 26 January 2005 / Accepted 29 March 2005
Abstract
This paper discusses the properties of the inner puffed-up rim that
forms in circumstellar disks when dust evaporates. We argue that the
rim shape is controlled by a fundamental property of circumstellar
disks, namely their very large vertical density gradient, through the
dependence of grain evaporation temperature on gas density. As a
result, the bright side of the rim is curved, rather
than vertical, as expected when a constant evaporation
temperature is assumed. We have computed a number of rim models that
take into account this effect in a self-consistent way. The results
show that the curved rim (as the vertical rim) emits most of its
radiation in the near and mid-IR, and provides a simple explanation for
the observed values of the near-IR excess (the "3 m bump'' of
Herbig Ae stars). Contrary to the vertical rim, for curved rims the
near-IR excess does not depend much on the inclination, being maximum
for face-on objects.
We then computed synthetic images of the curved rim seen
under different inclinations;
face-on rims are seen as bright, centrally symmetric rings
on the sky; increasing the inclination, the rim takes an elliptical
shape, with one side brighter than the other.
Key words: accretion, accretion disks - radiative transfer - circumstellar matter - planetary systems: protoplanetary disks - stars: pre-main sequence - infrared: stars
The structure of the inner regions of circumstellar disks associated with pre-main sequence stars is the subject of intense research. Interferometers working in the near infrared are providing the first direct information on the morphology of disks on scales of fractions of AU. They show that in the majority of cases the observed visibility curves are not well reproduced by flared disk models; rather, they are consistent with the emission of a ring of uniform brightness, of radius similar to the dust evaporation distance from the star (Millan-Gabet et al. 2001; Tuthill et al. 2001).
The interferometric results provide strong support for the idea that the inner disk structure deviates substantially from that of a flared disk because dust evaporation introduces a strong discontinuity in the opacity, which results in a "puffed-up'' rim at the dust destruction radius, where dust is exposed directly to the heating stellar radiation. The idea of a puffed-up inner rim was proposed by Natta et al. (2001) and developed further by Dullemond et al. (2001, hereafter DDN01) for Herbig Ae stars, to account for the shape of the near-infrared excess of these stars (the "3-m bump''). These authors pointed out that the rim also had the right properties to explain the early interferometric results of Millan-Gabet et al. (2001). Recent theoretical work by Muzerolle et al. (2004) has shown that the condition required to produce a puffed-up inner rim are indeed likely to exist in most Herbig Ae and T Tauri stars. The concept of such an inner rim has been widely used to interpret near-IR interferometric data for Herbig and T Tauri stars (Eisner et al. 2004; Muzerolle et al. 2003; Colavita et al. 2003; Eisner et al. 2003; Monnier & Millan-Gabet 2002; Millan-Gabet et al. 2001). Its effects on the disk structure and emission at larger radii have been discussed by Dullemond & Dominik (2004), who propose that the classification of Herbig Ae stars in two groups, based on the shape of the far-infrared excess (Meeus et al. 2001), can be interpreted as differences between objects where the outer disk emerges from the shadow of the inner rim and objects where this never happens.
In spite of its success in accounting for a variety of observations, the actual structure of the rim has not been much discussed. DDN01 adopted for their models a very crude approximation, namely that the illuminated side of the rim is "vertical'', and that its photospheric height is controlled by radial heat diffusion behind the rim. Such a model, taken at face value, has the obvious disadvantage that the rim emission vanishes for objects seen face-on, for which the projection on the line of sight of the rim surface is null, and for objects seen edge-on, where the rim obscures its own emission. This is clearly inconsistent with observations of the SED, which show that all the Herbig Ae stars with disks have similar near-IR excess, regardless of their inferred inclination (Natta et al. 2001; Dominik et al. 2003).
The vertical shape of the illuminated face of the rim is clearly not physical, as pointed out already by DDN01. Several effects are likely to "bend'' the rim: among them, one can expect that radiation pressure on dust grains or dynamical instability, due to self-shadowing effects, could modify the illuminated face of the rim (Dullemond 2000; DDN01). None of these suggestions, however, has been explored further.
In this paper, we will discuss in detail a different process, not mentioned so far, which depends exclusively on the basic physics of dust evaporation, i.e., on the dependence of the evaporation temperature on gas density. Circumstellar disks are characterized by a very large variation of the density in the vertical direction, so that the dust evaporation temperature varies by several hundred degrees in a few scale heights; moving vertically away from the disk midplane along the rim, dust will evaporate at lower and lower temperatures, i.e., further away from the central star. This very simple effect curves significantly the inner face of the rim, as we will describe in the following.
The paper is organized as follows. In Sect. 2 we describe the model we use to compute the rim shape and its observational properties. The results are presented in Sect. 3, where we discuss also how the rim depends on dust properties. A discussion of the results, in view of the existing observations, follows in Sect. 4, and a summary is given in Sect. 5.
Our model of the inner rim of passive-irradiated flaring disks joins the two different analytical methods to solve the structure of a circumstellar disk proposed respectively by Calvet et al. (1991, 1992, hereafter C92) and Chiang & Goldreich (1997, hereafter CG97). The temperature in the rim atmosphere is determined using the analytical solution of the problem of the radiation transfer as in C92, neglecting the heating term due to the mass accretion. The vertical structure of the rim is then computed in a way derived from CG97 and DDN01, adding a relation between the dust vaporization temperature and the gas density as proposed in Pollack et al. (1994). As a result we obtain a curved model for the inner rim whose features are described in the next paragraph.
Although the expressions for the dust temperature derive from a first order solution of the radiation transfer equation, we found (see Appendix) good agreement with the correct numerical result in most cases. To zero order, we can compute the rim structure avoiding proper radiation transfer calculations.
In the limit where the incident angle of the stellar radiation onto the disk surface is , the equations of the temperature for and (where is the optical depth for the emitted radiation) are formally equal to the optically thin and midplane temperatures introduced by CG97 (two-layer approximation). Therefore, while the two layer approximation is useful to study the structure and the emission features (e.g. silicate features at ) of the flaring part of the disk (as in DDN01), it must be abandoned in modeling the inner rim, since . Nevertheless we can adapt the CG97 treatment of the vertical structure of the disk to the inner rim using the appropriate expressions for the dust temperature.
In order to clarify this concept and to introduce the relation between the vaporization temperature of dust and the gas density, the basic equations are briefly summarized. We refer to the cited works for a physical discussion of the equations.
We suppose that the disk is heated only by the stellar radiation and
we call
the incident angle beetwen the radiation and the disk
surface. The incident beam is absorbed exponentially as it penetrates
the dusts and if
is the optical depth for the emitted
radiation, the dust temperature
is given by the relation (C92)
In the flaring part of the disk, for which ,
Eq. (1) can be rewritten in terms of the two
layer approximation, proposed in CG97, in which the interior of the
disk (with
)
is heated to the temperature
by
half of the stellar flux, while the other half of the stellar flux is
reemitted by the superficial layer heated at the optically thin
temperature :
Since the dust in the rim is by
definition close to the evaporation temperature, the value of
is fixed and depends only on the grain absorbing cross
section. Moreover, to first order, since the difference between
and
is never very large, at the inner radius ()
the ratio between these two temperature depends only on
and is given by:
Figure 1: Width of the region over which 99% of the stellar radiation is absorbed as function of . | |
Open with DEXTER |
With the assumption that the disk is in hydrostatic equilibrium in the
gravitational field of the central star and that is isothermal in the
vertical direction z, the gas density distribution is expressed by
the Gaussian relation
Note that the interior temperature of the disk
depends on
the incident angle
between the stellar radiation and the disk
surface, through Eq. (5). The quantity
accounts for the projection of a disk annulus on a
plane perpendicular to the incident radiation and is given by the
relation
Figure 2: Photospherical height H ( ) of the inner rim for different values of , as labelled. The curves have been computed for a star with K, , . Each curve ends where the rim becomes optically thin to the stellar radiation. | |
Open with DEXTER |
We can follow CG97 and DDN01 to obtain a self-consistent solution of
Eqs. (2), (5), (6), (8),
(9), (11), (13) to determine the
structure of the disk, as long as dust evaporation can be
neglected. When dust evaporation is important, DDN01 have developed an
approximate solution of the rim/disk structure under the assumption of
constant
.
The inner rim has a vertical surface toward the
star located at the dust evaporation distance
.
The vertical
photospherical height H depends on the dust density behind the rim;
since the rim is higher than the flaring disk, it casts a shadows over
the disk. In this shadowed region, assuming that there is no external
heating except the stellar radiation, the pressure height h depends
only on the radial heating diffusion. The exact determination of the
structure of the rim in this diffusive region would thus require to
solve the problem of radiation transport in two dimensions (see
Dullemond 2002). Since this goes well beyond our aims, we adopt
the approximated relation used in DDN01
We now introduce the relation between the evaporation temperature of
the dust and the gas density. The physical reasons for this effect
can be easily understood thinking of evaporation as the process
by which equilibrium
between the gas pressure and the surface tension of the dust grains is
reached: the higher the gas density, the higher will be the evaporation
temperature. We adopt for the dust in the disk the model proposed by
Pollack et al. (1994). In this model, the grains with the
higher evaporation temperature are the silicates, which will therefore
determine the location of the rim. Their evaporation temperature (see
Table 3 in Pollack et al.) varies with the gas density roughly as a
power law, of the kind:
To obtain a self-consistent determination of the structure of the curved rim, we implemented a numerical method able to solve Eq. (16) together with the set of Eqs. (2), (5), (6), (8), (9), (11), (13). The distance of dust evaporation in the midplane is computed for z=0 in the set of equations and is taken as the starting point of the radial grid on which the rim structure is computed. As a result of the calculations, we obtain the location in the (R,z) plane of the rim surfaces characterized by a constant value of the optical depth. The surface for is thus the evaporation surface of dust grains, for we obtain the surface relative to the photospheric height H, defined through the Eq. (13), while for we obtain the emitting surface of the rim.
Moving on the rim surface away from the star, the incident angle decreases and, when it approaches zero, the determination of the rim shape becomes very difficult. This is mainly due to the fact that the described solution for the radiation transfer neglects the heat diffusion between contiguous annulus of the rim. Therefore both the midplane temperature and the pressure height h of the rim goes unphysically to zero for , according to Eqs. (5) and (9). To avoid this unrealistic behaviour, we use the approximated relation discussed previously (see Eq. (15)) to determine the pressure height h in the region where the diffusion is the dominant heating source. The transition distance between the region of the rim heated by the star and those heated by the diffusion is determined imposing continuity of . The photospherical height H is then determined as for the vertical rim.
Using the model described in the previous section, we compute the structure of the inner rim for a disk heated by a star with temperature K, mass and luminosity . We take a disk surface density , and a dust-to-gas mass ratio dust/gas = 0.01. In our models, this value of corresponds to a midplane gas density of about at 0.5 AU from the star. Note that the results are not very sensitive to the exact value of , as long as the inner disk remains very optically thick.
The dust properties are those of the astronomical silicates of Weingartner & Draine (2001). In our models, we consider that all the grains have the same size and characterize their properties with the quantity , the ratio of the mean Planck opacity at the evaporation temperature to that at the stellar temperature (see Eq. (2)). The evaporation temperature of silicate grains vary from 1600 K, for gas density of , to 1000 K for (see Eq. (16)). For the Weingartner silicates and a vaporization temperature K, for grains of radius , and grows to values of about unity for grain radii > .
The shape of the rim is shown in Fig. 2, which plots the locus of (i.e., the photospheric height H) as function of Rfor different values of . The ending point of the rim ( , ) is when the rim becomes optically thin at the stellar radiation.
For , corresponding to silicate grains bigger than , the inner radius , the outer radius and the maximum photospheric height of the rim , all vary very little, with values AU, AU and AU. For smaller values of , the rim becomes steeper and the inner radius increases. For , corresponding to grains with radius , the rim has AU, AU and AU.
The shape of the rim can be roughly characterized by the ratio of its maximum height over width . This quantity, which is nominally infinity in a vertical rim, becomes in the curved rim models 0.6 for and is 0.8 for . In other words, as grains grow, the inner rim approaches the star but the bending of the surface varies very little.
For any given value of ,
the dust temperature along the rim
is not constant (as in the vertical rim) but decreases from values of
about 1400 K, typical of silicate evaporation temperatures at
density of
,
to 1200 K for density of
.
Figure 3 plots the effective temperature of the rim
(i.e.,
)
along its surface for different values of .
For
the effective
temperature is equal to the vaporization temperature; for
the effective temperature is generally lower since the
dust temperature decreases with the optical depth, according to
Eq. (1). This variation, however, is sufficiently small
that most of the rim emission occurs in the near-IR, as for the
vertical rim.
Figure 3: Effective temperature (i.e., ) along the rim surface for different values of , as labelled. Each curve is plotted as function of , where is the distance of the rim from the star on the disk midplane. Note that, in fact, not only R but also z changes along each curve, as shown in Fig. 2. For is almost independent of the grain opacity. | |
Open with DEXTER |
The fraction of stellar luminosity intercepted by the rim, given by the
ratio
,
varies from ,
for small values of
,
to ,
for
.
Assuming that the rim is in thermal equilibrium, the intercepted
radiation is equal to the total emitted flux
.
For each
,
we have computed the spectral energy distribution of the
rim emission for different inclination angles. As discussed in
Sect. 2, the rim emission is computed as that of a blackbody
at the local temperature along the rim
surface. We will
come back in the Appendix to this assumption, which, in any case,
gives a very good approximation to the global properties of the rim
emission (see also Muzerolle et al. 2003, C92). Most of the
emission, as expected from the range of temperatures, occurs in the
near-IR. We have computed the fraction of the stellar luminosity
re-emitted by the rim in the wavelength range 1.25-7.0 m for
different values of
as function of the inclination
(Fig. 4). This near-IR excess peaks at zero
inclination, where has values between
(smaller grains) and
(larger grains). As the inclination increases, the near-IR
excess decreases slowly, reaching values between 5% and 8%,
depending on .
For inclination higher then
the rim emission is self-absorbed. Note that for the large grains,
with
,
the near-IR excess becomes almost
independent of .
Figure 4: Near infrared emission of the curved rim (integrated between and normalized to ) for different inclination angles and different values of . Starting from the bottom, the curves refer to (long-dashed), 0.14 (short-dashed), 0.25 (dotted), (short-dash-dotted), 0.92 (long-dash-dotted). For the last two values of the emission is almost the same and the lines overlap. The continuum line plots the emission of the vertical rim, computed for . | |
Open with DEXTER |
The behaviour described in Fig. 4 is very different from what one obtains in the case of the vertical rim. As described in the previous section, neglecting the dependence of the evaporation temperature of grains from gas density results in the inner face of the rim to be vertical (as in DDN01). With and (equal to the evaporation temperature at the inner radius of the curved rim), the inner radius is the same as for the curved rim; the photospheric height, evaluated using the approximation described by Eq. (15), is also the same ( AU); the fraction of stellar luminosity intercepted by the vertical rim is . However, the value of observed for different inclination angles is very different, with a very strong (and opposite) dependence on i. (Fig. 4). In particular, vanishes for face-on disks, and is maximum (of the order of 20%) for very inclined systems, just before self-absorption sets in.
There are also differences between the curved and vertical
rim models concerning the shape of the predicted 3 m bump.
While the vertical rim has a constant temperature of 1400 K over all
its surface, in the curved rim the temperature varies from 1400 K on
the midplane to about 1200 K at the outer edge (see Fig. 3) and the SED is broader than a single-temperature black
body. In practice, however, this is only a minor effect.
Figure 5: Synthetic images of the curved rim for different values of and inclination angle ( is face-on); comparison with the vertical rim. The left and central panels show the images of the curved rim calculated for (silicate grains with radius ) and ( ). The right panels show the images for the verical rim calculated for . The surface brightness of the rim, plotted in colors, is computed for a wavelength of . The stellar parameters are as in Sect. 3, the distance is . The vertical rim for has zero brightness. Note that as discussed in Sect. 3 for the images of the curved rim remain the same. | |
Open with DEXTER |
Figure 5 shows synthetic images of the curved rim at m for grains with (radius of about 0.1 m) and (radius of about 1.3 m). As described in Sect. 3.1 (see Figs. 2 and 3), the inner radius of the rim is larger for smaller grains and the surface brightness is lower, due to the lower effective temperature of the emitting surface.
For the comparison with the vertical model of the rim, the left panels of Fig. 5 show the images of the vertical rim, calculated for . The largest difference is at low inclinations: for (face-on disks) the vertical rim vanishes, as the projected emitting surface along the line of sight is zero, while the curved rim has a centrally symmetric ring shape. For a distance of 144 pc, the inner radius of the ring, for , is about 4 milliarcsecond (mas); its brightness peaks practically at ( ) and decreases slowly (by about a factor of two) outward, until at it drops to very low values. At this distance, the width of the rim (roughly its FWHM) is about 0.8 mas.
For higher inclination the central symmetry is lost and the projected image is an ellipse with one edge brighter than the other. In general, however, the brightness distribution of the rounded rim is much more symmetric than that of the vertical rim.
The results shown so far have been computed assuming that all the grains have the same composition and size. This is in practice unrealistic, and one wonders how the results will change if grains with different properties are present. Different grains absorb differently the stellar radiation and reach very different equilibrium temperatures. The degree of complexity of the radiation transfer problem increases remarkably (see Wolfire & Cassinelli 1986, 1987) and there is no approximate solution such as Eq. (1). In fact, to the best of our knowledge even numerical solutions are not available to describe accurately the transition region where some grains are cooler than evaporation and survive while others do not.
We can however estimate the effect of a grain mixture in the following way. Let us consider for simplicity a classical MRN grain size distribution, characterized by a power law a^{-q} with a the radius of the dust grains, q=3.5, and a varying between a minimum value m and a maximum value . From what we have discussed in Sect. 2, if is smaller than (corresponding to for the adopted silicate grains), all the grains have a maximum temperature at and the higher the value of , the smaller the evaporation distance from the star. As soon as the largest grains can survive, the stellar radiation is rapidly absorbed and all the other grains will also survive. In this case, the shape of the inner rim is controlled by the largest grains in the distribution. The shapes in Fig. 2 and the emitted fluxes shown in Fig. 4 should not change significantly as long as one interprets as relative to the largest grains.
The situation is more complex if , since the grains with can survive near the star only in an optically thin regime. However, our single-grain models show that there is very little difference in the location and shape of the rim as soon as . This suggests that varying above the "transition'' value (i.e., the value of at which the dependence of T on changes from decreasing to increasing) will not change the rim properties any further, at least to zero order.
It is also likely that grains have not just different size but also different chemical composition. In this case, the rim location and properties are determined by the dust species with the highest evaporation temperature, as long as its contribution to the disk opacity is sufficient to make the disk optically thick. As discussed in Sect. 2, in the Pollack et al. (1994) dust model silicates have the highest evaporation temperature. This is why we have performed all our calculations for silicate grains. However, one should keep in mind that, in some dust models, graphite contributes most of the opacity at short wavelengths; since graphite has a much higher evaporation temperature ( K) than silicates, its presence would move the rim much closer to the central star. Although the details of the rim shape may change, its curving, which is caused only by the dependence of on the density, will not disappear.
The curvature of the inner side of the rim has important consequences on the observable near-IR excess, as shown in Fig. 4, which shows that the predicted near-IR excess ranges between about 10 and 20%, and that there is no significant dependence on the inclination of the disk with respect to the observer. These results compare very well with the existing observations. Figure 6 shows the observed values for a well-studied sample of Herbig Ae stars from Natta et al. (2001) and Dominik et al. (2003); of a total of 16 objects, only one (HD 142527) has , and one (HD 169142) <0.09.
The observations of the near-IR SEDs of Herbig Ae stars have also shown that there is no systematic variation of the near-IR excess with the inclination of the disk. The sample of Fig. 6 includes some rather face-on objects such as AB Aur ( , ; Fukagawa et al. 2004) and HD 163296 ( , ; Grady et al. 2000) and some UXORS variable, such as UX Ori, WW Vul and CQ Tau, which are generally considered to be close to edge-on, and have values of . This, which has been a puzzle for the vertical rim, finds a natural physical explanation in our models, which predict that the near-IR excess depends little on the inclination and that, in particular, does not vanish for face-on objects. Actually, our models predict that the largest excesses should be seen in face-on objects, and it would be interesting to explore in detail if this is indeed the case. However, for such a study to be significant, one would need a large sample of objects with known inclinations, which is at present not available.
Another interesting aspect of our results is the dependence of
on grain properties. If taken at face-value, one
should expect that only objects with relatively large grains (greater
than about 1 m) can have large values of
,
and
that low values of the near-IR excess can only occur in face-on
objects with small grains. These properties of the rim emission are
potentially important for a better understanding of grain properties
in the inner disk and should be pursued further in the future,
combining observations of disks with well known inclinations with
models that explore a broader range of grain properties than we
considered in this paper.
Figure 6: Histogram of the observed values of for a sample of 16 Herbig Ae stars from Natta et al. (2001) and Dominik et al. (2003). Note that the Dominik et al. values are computed between 2 and 7 m, while those from Natta et al. between 1.25 and 7 m. The uncertainties are in most cases rather large, due to the subtraction of the photospheric flux, which contributes significantly to the observed ones at short wavelengths, and to variability (see Natta et al. for some examples). | |
Open with DEXTER |
One side product of our models is
,
the distance from the star to the
rim on the disk midplane. For a face-on disk,
practically coincide with the position of the peak of the near-IR brightness.
For fixed values of the stellar parameters,
depends on the dust properties; its value
is shown in Fig. 7
for different values of .
Larger grains
can survive closer to the star than smaller grains and, because of the
inversion in the temperature gradient,
does not change much
with the grain size once this exceeds the value for which
,
so that for
,
AU.
scales roughly as
,
so that one can expect
AU for
.
These values, especially for small grain sizes, are larger
than the predictions of simple, optically thin calculations, because
Eq. (1) correctly includes the effect of the diffuse radiation
field. Note that if backward scattering cannot be neglected, the
temperature at
will be even higher, increasing further the
size of the inner hole (see Appendix).
Figure 7: Behaviour of the inner radius of the rim for different values of , evaluated for the star and disk parameters described in Sect. 3. Note that for (corresponding to silicate grains of m) the radius of the inner rim is almost constant around 0.5 AU. | |
Open with DEXTER |
Values of have been derived for a number of objects from near-IR interferometric observations (see, for example, Akeson et al. 2000, 2002; Millan-Gabet et al. 2001; Monnier & Millan-Gabet 2002; Colavita et al. 2003; Eisner et al. 2003, 2004, 2005). The results depend not only on the assumed disk model, but also on other quantities such as the stellar SED, the disk inclination etc., so that a direct comparison with our values of can be misleading. We note, however, that AU for intermediate-mass objects is roughly consistent with the observations, suggesting that in many objects grains have grown to sizes m (see also Monnier & Millan-Gabet 2002).
This conclusion, of course, needs to be taken with great caution. There are a number of effects that can change the model-predicted , in addition to different grain sizes. For example, on one hand smaller values of can be due to some low-density gas in the inner disk hole, able to absorb the UV continuum from the star (Monnier & Millan-Gabet 2002; Akeson et al. 2005). Also, the presence of graphite, which has an evaporation temperature of about 2000 K, can decrease . On the other hand, the presence of accretion at a significant rate can increase by heating grains in the inner disk to temperatures higher than those produced by the photospheric radiation alone (Muzerolle et al. 2004).
In practice, one needs to compare in detail the model predictions of the quantities measured with interferometers (visibility curves, their dependence on baseline and hour angle, phase measurements) for specific, well-known objects. We are currently performing such a study with the aim of assessing the impact of curved rim models on the interpretation of current (and future) interferometric observations. In particular, we will include a discussion of the effect of the asymmetries of the rim projected images when not face-on. The results will be presented in a forthcoming paper (Isella et al. in prep.).
To the best of our knowledge, the only image of the inner region of a circumstellar disk is that obtained by Tuthill et al. (2001), for LkHa101 using Keck in the H and K bands. The results are reminescent of our images, showing an elliptical ring with a side much brighter than the other. The similarity is very interesting but one should keep in mind that LkHa101 is an early B star with a luminosity between 500 and 50 000 , depending on its distance. Our model may not apply to such an object (e.g., Monnier et al. 2005).
In this paper, we discuss the properties of the inner puffed-up rim that forms in circumstellar disks when dust evaporates. The existence of a rim has been claimed starting from the work of Natta et al. (2001), both on theoretical and observational grounds. Here, we investigate the shape of the illuminated face of the rim. We argue that this shape is controlled by a fundamental property of circumstellar disks, namely their very large vertical density gradient, through the dependence of grain evaporation temperature on density. As a result, the bright side of the rim is naturally curved, rather than vertical, as expected when a constant evaporation temperature is assumed.
We have computed a number of rim models that take into account this effect in a self-consistent way. A number of approximations have been necessary to perform the calculations, and we discuss their validity. The basic result (i.e., the curved shape of the rim illuminated face) appears to be quite robust.
For a given star, the rim properties depend mostly on the properties of the grains, and very little on those of the disk itself, for example the exact value of the surface density. The distance of the rim from the star is determined by the evaporation temperature (at the density of the disk midplane) of the dust species that has the highest evaporation temperature, as long as its opacity is sufficient to make the disk very optically thick; in the model of Pollack et al. (1994) of the dust in accretion disks, silicates have the highest evaporation temperature. Therefore, we have assumed in our models dust made of astronomical silicates, and varied their size over a large range of values. We find that the rim properties do not depend on size as soon as m; the values of the rim radii observed with interferometers suggest that in many pre-main-sequence disk grains have grown to sizes of 1-few m at least.
The curved rim (as the vertical rim) emits most of its radiation in the near and mid-IR, and provides a simple explanation for the observed values of the near-IR excess (the "3 m bump'' of Herbig Ae stars). Unlike the vertical rim, for curved rims the near-IR excess does not depend much on the inclination, being maximum for face-on objects and only somewhat smaller for highly inclined ones. This is in agreement with the apparent similarity of the observed near-IR SED between objects seen face-on and close to edge-on.
We have computed synthetic images of the curved rim seen under different inclinations. Face-on rims are seen as bright, centrally symmetric rings on the sky; increasing the inclination, the rim takes an elliptical shape, with one side brighter than the other. However, the brightness distribution of curved rims remains at any inclination much more centrally symmetric than that of vertical ones. In a forthcoming paper (Isella et al. in preparation) we will discuss the application of the curved rim models to the interpretation of near-IR interferometric observations of disks.
Acknowledgements
We would like to thank Endrik Krügel for having allowed us to make use of his radiation transfer code, Leonardo Testi, Kees Dullemond, Carsten Dominik, Giuseppe Lodato and Giuseppe Bertin for useful discussions. The authors acknowledge partial support by MIUR COFIN grant 2003/027003-001.
In this appendix we discuss in detail some of the assumptions on
radiation transfer made in building our rim models. We use as
templates the results of a radiation transfer code developed by
E. Krügel, which is described in Habart et al. (2004). The code
considers a plane-parallel slab of dust, illuminated on one side by a
star with properties as in Sect. 3.
Figure A.1: Temperature as function of the optical depth to the stellar radiation for small grains ( ) and very large grains ( ). The dashed curves are the results of the radiation transfer code, the solid lines the temperature obtained from Eq. (1). |
We show in Fig. A.1 the comparison of the temperature derived from Eq. (1) with the results of the radiation transfer code for two cases, one with (small grains) and one with (very large grains). In both cases scattering is neglected. The temperature is plotted as function of the optical depth to the stellar radiation . One can see that Eq. (1) gives the correct values of T for and for . For intermediate values of , the results are in both cases accurate within 10%. This is more than adequate for the purposes of this paper.
Figure A.2: Comparison of the SED computed with the radiation transfer code (dotted lines) and as a black-body at from Eq. (1) (solid line) for and , as labelled. |
Figure A.3: Temperature profile for grains with , when scattering is included ( ; dashed line) or suppressed ( ; solid line). |