Issue |
A&A
Volume 512, March-April 2010
|
|
---|---|---|
Article Number | A11 | |
Number of page(s) | 5 | |
Section | Interstellar and circumstellar matter | |
DOI | https://doi.org/10.1051/0004-6361/200912743 | |
Published online | 23 March 2010 |
Full two-dimensional radiative transfer modelling of the transitional disk LkCa 15
(Research Note)
G. D. Mulders1,3 - C. Dominik1,4 - M. Min2
1 - Astronomical Institute ``Anton Pannekoek'', University of Amsterdam, PO Box 94249, 1090 GE Amsterdam, The Netherlands
2 -
Astronomical institute Utrecht, University of Utrecht, PO Box 80000, 3508 TA Utrecht, The Netherlands
3 -
SRON Netherlands Institute for Space Research, PO Box 800, 9700 AV, Groningen, The Netherlands
4 -
Afdeling Sterrenkunde, Radboud Universiteit Nijmegen, Postbus 9010, 6500 GL Nijmegen, The Netherlands
Received 22 June 2009 / Accepted 5 January 2010
Abstract
Context. With the legacy of Spitzer and current advances in
(sub)mm astronomy, a considerable number of so-called ``transitional''
disks has been identified which are believed to contain gaps or have
developped large inner holes, some filled with dust. This may indicate
that complex geometries may be a key feature in disk evolution that has
to be understood and modelled correctly. The disk around LkCa
15 is such a disk, with a large gap ranging from 5-46
AU, as identified by Espaillat et al. (2007, ApJ, 670, L135) using 1+1D
radiative transfer modelling. To fit the spectral energy distribution
(SED), they propose two possible scenarios for the inner
(<5 AU) disk - optically thick or optically thin - and one
scenario for the outer disk.
Aims. We use the gapped disk of LkCa 15 as a case in point to
illustrate the importance of 2D radiative transfer in transitional
disks by showing how the vertical dust distribution in dust-filled
inner holes determines not only the radial optical depth but also the
outer disk geometry.
Methods. We use MCMax, a 2D radiative transfer code with a
self-consistent vertical density and temperature structure, to model
the SED of LkCa 15.
Results. We identify two possible geometries for the inner and
outer disk that are both different from those in Espaillat et al.
(2007, ApJ, 670, L135). An inner disk in hydrostatic equilibrium
reprocesses enough starlight to fit the near infrared flux, but also
casts a shadow on the inner rim of the outer disk. This requires the
outer disk scale height to be high enough to rise out of the shadow. An
optically thin inner disk does not cast such a shadow, and the SED can
be fitted with a smaller outer disk scale height. For the dust in the
inner regions to become optically thin however, the scale height would
have to be so much higher than its hydrostatic equilibrium value that
it effectively becomes a dust shell. It is currently unclear if a
physical mechanism exists which could provide for such a configuration.
Conclusions. We find that the radial optical depth of dust
within the inner hole of LkCa 15 is controlled by its vertical
distribution. If it turns optically thick, the outer disk scale height
must be increased to raise the outer disk out of the inner disk's
shadow.
Key words: stars: pre-main sequence - stars: individual: LkCa 15 - protoplanetary disks - radiative transfer - circumstellar matter
1 Introduction
Protoplanetary disks are the main sites of planet formation. Within them, small dust grains grow from micron to millimeter sizes and settle to the midplane where they are believed to eventually form planetary systems. Sufficiently massive planets locally deplete the disk of gas and dust and may open a gap or enlarge the inner hole. Recently, considerable numbers of so-called transitional disks have been identified on the basis of their low near infrared excess (Najita et al. 2007), indicating that their inner regions may be (partially) depleted of dust.LkCa
15 is such a transitional disk, but with a more complex
geometry: it has an inner disk, a gap and an outer disk. The outer
disk, seen at an inclination of 51,
extends from 46 to 800 AU (Piétu et al. 2006). The Spitzer data show that the inner hole is not devoid of dust between 0.1 and 5 AU (Espaillat et al. 2007,
hereafter Esp07). On the basis of 1+1D RT modelling,Esp07 identify two
possible scenarios for the inner regions: one where the dust is
optically thick and one where it is optically thin. Additional
observations show that this dust is most likely present in the form of
a small optically thick inner disk at the dust evaporating radius (Espaillat et al. 2008, hereafter Esp08). The gap inbetween is mostly empty.
In this paper we will model both geometries proposed by Esp07 using a full 2D radiative transfer code and focus on the vertical structure of the inner and outer disk. Our modelling approach differs from that of Esp07 in one important point: We combine the several disk components in one model before we perform the radiative transfer and create the spectral energy distribution (SED), whereas they create the SED after performing the radiative transfer on each seperate component. Although their approach is valid if the seperate components do not influence each other, we will show that at least for LkCa 15 this is not the case.
2 Disk model
The disk model we use here is MCMax (Min et al. 2009),
a 2D Monte Carlo radiative transfer code with a selfconsistent vertical
structure. The Monte Carlo radiative transfer simulation is based on
that of Bjorkman & Wood (2001)
and calculates the dust temperature for a given density structure. It
also calculates the vertical structure by directly integrating the
equation of hydrostatic equilibrium using the Monte Carlo temperature.
Because the vertical structure calculation requires the temperature in
the midplane to be accurately determined, MCMax treats those regions
with an analytical diffusion approximation (see Min et al. (2009)
for details). By iterating between the Monte Carlo simulation and the
hydrostatic equilibrium calculation a self-consistent solution is
found. This typically takes four or five iterations. MCMax also
includes a modified random-walk approximation that speeds up the
calculations in regions with extreme optical depths by making multiple
interaction steps in a single computation (Min et al. 2009, see also Fleck & Canfield 1984).
The typical run time is less than three minutes on a single core CPU
for a model with a radial optical depth in the midplane of
using six iterations with
photon packages each.
Because models with well-mixed dust and gas in vertical hydrostatic
equilibrium tend to overpredict the disk's surface height and far
infrared flux, dust settling and grain growth are necessary ingredients
to fit SEDs (e.g. D'Alessio et al. 2006).
For the present study we will describe them with two parameters that
lower the disk surface: The mass fraction of large grains
and a sedimentation parameter
.
The sedimentation parameter
describes the settling of small dust particles and is implemented by reducing the dust scale height with respect to the selfconsistent gas scale height:
.
A similar global sedimentation parameter has also been used by Ratzka et al. (2007).
Both parameters have a different effect on the surface shape and SED
and are not degenerate. The sedimentation parameter is global and
reduces the scale height in both the inner and outer disk and therefore
scales the SED accordingly. The fraction of large, mm-sized grains only
affects the surface height in the outer disk: increasing
reduces the flaring angle of the outer disk and mostly affects the far infrared flux.
We will use a bimodial size distribution here, consisting of small (0.1
)
and big (2 mm) grains. For the dust opacity of the small particles we use that of 0.1
grains, consisting of 80% silicate (Dorschner et al. 1995) and 20% amorphous carbon (Preibisch et al. 1993). We do not model them as compact spheres, but as irregular shaped particles (Min et al. 2005).
The opacity of the big particles is modelled using 2 mm sized irregular
shaped silicate particles. The total opacity is determined by the
fitting parameter
,
the fraction of mass in big particles. To settle the big grains to the midplane, we introduce a stratification parameter
that reduces the scale height of large grains with respect to that of small grains:
.
Table 1: Stellar and best-fit parameters for the two different geometries.
For the stellar photosphere we use a main-sequence Kurucz model with
the stellar mass, luminosity and effective temperature from Kenyon & Hartmann (1996) (Table 1). This photosphere is fitted to the de-reddened photometry between 0.5
and 2
by scaling the distance.
3 Results
For each of the geometries (Fig. 1) we follow a different approach. In the optically thick case, we model the inner disk first to constrain the size of the shadow cast on the outer disk and proceed by modelling the outer disk. Note that although we change the fit parameters for one component at a time, we always perform the radiative transfer including all components. In the optically thin case, the inner disk produces no shadow and we can model the outer disk first. We then proceed by including an inner disk so we can explore the effect of the vertical dust distibution on the radial optical depth.
3.1 Optically thick inner disk
![]() |
Figure 1: Sketches (not to scale) of the two possible geometries for LkCa 15 from this paper ( top panel) and - for comparison - that from Esp07 ( bottom panel). Depicted are: a) A model with a 1 AU optically thick inner disk and an outer disk that rises out of its shadow. b) A model with a 5 AU optically thin spherical dust shell and a flatter outer disk with a fully illuminated inner rim. c) The model proposed by Esp07, a combination of the inner disk from a) with the fully illuminated outer disk of b). |
Open with DEXTER |
![]() |
Figure 2:
Spectral Energy Distribution of LkCa 15 (diamonds) with 1 |
Open with DEXTER |
![]() |
Figure 3:
Radial |
Open with DEXTER |
In the optically thick case, the dust mass in the inner disk is hard to
constrain from the SED. Therefore we use the same surface density
powerlaw (p=1) as for the outer disk, which yields
,
assuming a dust-to-gas ratio of 1:100. We fix the inner radius at the
dust evaporating radius at 0.1 AU and fit an outer radius for the inner
disk of 1.0 AU, yielding a dust mass for the inner disk of 4
.
Our best fit (Fig. 2)
does not require any sedimentation, and does not constrain
stratification or big particle fraction, for which we will use the fit
values from the outer disk.
This choice of inner disk parameters is not unique, but the
size of the shadow cast on the outer disk is. The height at the outer
edge of the inner disk (
,
Fig. 3)
sets the total amount of energy reprocessed by the inner disk, which is
constrained by the SED. The size of the shadow is therefore
approximately equal for all solutions and independent of the exact
inner disk geometry.
With the size of the shadow cast by the inner disk constrained,
we proceed by fitting the outer disk. We find a disk mass of 0.3
using p=1. This disk mass is much higher than previously reported, for instance by Andrews & Williams (2005).
This can be attributed to the millimeter opacity of our big particles,
which is a factor of five lower at 1.3 mm. Corrected for this, the
disk mass for LkCa
15 does not differ with from that of Andrews & Williams (2005).
We lower the outer disk surface by using a sedimentation parameter of
and a big particle fraction of 94%. We also settle the big grains to the midplane (
)
to fit the relatively low flux at 350
.
A model without stratification (
)
overpredicts this flux a little, but is still consistent with the large
errors at that wavelength. The inner rim of the outer disk is almost
completely shadowed by the inner disk, and only the disk surface is
illuminated (Fig. 3).
The reduction in surface height compared to a well-mixed model without
big particles is four, comparable to that of other T-Tauri stars (D'Alessio et al. 2006).
3.2 Optically thin inner disk
To investigate the effects of the optical depth of the inner disk on
the outer disk structure, we first construct a disk model without an
inner disk. In this model, the inner rim of the outer disk is fully
illuminated, and its contribution to the SED from 25
to 100
will be much larger compared to a shadowed outer disk. To fit the far infrared flux, we need stronger sedimentation (
)
and more big particles (
= 99%), so the outer disk becomes twice as flat (Fig. 3,
dotted line). We find a total reduction in the surface height of a
factor of eight, significantly higher than in other T-Tauri stars,
which was also reported by Esp07.
To fit the complete SED, we follow Esp07 and add an inner disk
between 0.1 and 5 AU. If we assume this disk is in vertical
hydrostatic equilibrium, we fit a dust mass in small grains of 2
as constrained by the flux in the near and mid infrared. This inner disk however is not optically thin to starlight. Its radial
surface reaches the same height (
)
as the outer disk, casting a shadow over the entire outer disk which
results in a far infrared flux drop by a factor of ten (dotted line in
the lower panel of Fig. 2).
One way to minimize the effect of this shadow is to make the inner disk extremely flat (
)
(dashed line in the lower panel of Fig. 2).
Although this model fits the far infrared flux reasonably well, the
inner disk does not reprocess enough of the stellar radiation to
reproduce the near infrared excess. Another way to avoid the shadow is
to increase the disk scale height much beyond its hydrostatic
equilibrium value. By spreading the same amount of dust over a larger
solid angle, its radial optical depth decreases. For the material to
become optically thin, the scale height has to be so high that we will
model it as a spherical shell around the star (solid line in the lower
panel of Fig. 2), but the exact geometry is hard to constrain as long as the material remains optically thin. We need only 1
of dust in this shell because it does not shield itself and is
therefore hotter. It is also in our line of sight and reddens the
starlight by an additional Av=0.25.
4 Discussion and conclusions
We have studied the radiative transfer effects in the gapped disk of LkCa 15 using MCMax, a self-consistent 2D radiative transfer code. We show that only an optically thin dust shell or an optically thick inner disk can explain the SED in the near and mid infrared. The outer disk geometry differs significantly depending on the inner disk properties. With an optically thick inner disk, the outer disk must rise out of the shadow cast by the inner disk. With an optically thin dust shell, there is no shadow and the outer disk becomes twice as flat.
These models are similar to but not the same as those proposed in Esp07. The latter do not take into account the shadow cast by an optically thick inner disk on the outer disk. They find the same z/r for the inner disk rim (0.01 AU / 0.12 AU = 0.08) as for the outer disk rim (4 AU / 46 AU = 0.08). This means that the outer disk wall cannot be illuminated, as it lies in the shadow of the inner disk (see Fig. 2, dotted line). The model with an optically thin inner disk can be reproduced, with the exception that to become optically thin, its scale height needs to be so high that it effectively becomes a dust shell.
Esp08 find the near infrared excess to be a blackbody and conclude that the inner disk has to be optically thick. On the other hand, Vinkovic et al. (2006) concluded that compact dust halos can explain the near infrared bumps in Herbig Ae/Be stars and that their existence agrees well with interferometric data.
An open question remains as to why the inner disk would have so high a scale height that it becomes an optically thin dust shell. Such a high scale height for small dust grains could be due to a halo of unknown origin, maybe residual infall, but such an infall would not lead to a dust cloud confined to <5 AU. Another possibility could be violent scattering and disruption of planetesimals that could create a cloud of dust with a scale height determined by dynamic processes rather than the hydrostatic equilibrium of gas. However, there is currently no additional evidence supporting this scenario for LkCa 15. Furthermore an optically thick inner disk leads to a model with outer disk settling requirements similar to what is seen in other T Tauri stars. Therefore we select the model with the optically thick inner disk as the preferred one.
![]() |
Figure 4: 20 micron images of our best-fit disk models of LkCa 15. The left panel shows a model with an optically thick inner disk of 1 AU that casts a shadow on the rim of the outer disk. The right panel shows a 5 AU optically thin dust shell and a fully illuminated outer rim. Note that if the rim of the outer disk is shadowed, much more radiation emerges from the disk surface. |
Open with DEXTER |
To distinguish between geometries with a dust disk or a dust shell - or
anything inbetween - a near-infrared interferometer can be used to
directly probe the inner regions. But because a shadow from the inner
disk connects the inner and outer disk geometry, imaging of the outer
disk can also reveal the nature of the inner disk (Fig. 4).
If the inner disk is optically thick, radiation from the outer disk
will be more extended because the contribution from the shadowed inner
rim of the outer disk is small, while the disk surface has a larger
flaring angle and contributes to the 20 m
flux as well. For an optically thin inner dust shell, the emission will
be more centrally peaked because the outer disk is flatter and its rim
fully illuminated.
We can draw one general conclusion about transitional objects: the radial optical depth - and therefore the vertical structure - of dust within an inner hole determines the illumination of the outer disk and ergo its geometry. This is especially important in the light of observations with the Spitzer space telescope, where large numbers of disks with dust-filled inner holes have been identified (Najita et al. 2007).
AcknowledgementsThis research project is financially supported by a joint grant from the Netherlands Research School for Astronomy (NOVA) and the Netherlands Institute for Space Research (SRON). M. Min acknowledges financial support from the Netherlands Organisation for Scientific Research (NWO) through a Veni grant.
References
- Andrews, S. M., & Williams, J. P. 2005, ApJ, 631, 1134 [NASA ADS] [CrossRef] [Google Scholar]
- Bjorkman, J. E., & Wood, K. 2001, ApJ, 554, 615 [NASA ADS] [CrossRef] [Google Scholar]
- D'Alessio, P., Calvet, N., Hartmann, L., Franco-Hernández, R., & Servín, H. 2006, ApJ, 638, 314 [NASA ADS] [CrossRef] [Google Scholar]
- Dorschner, J., Begemann, B., Henning, T., Jaeger, C., & Mutschke, H. 1995, A&A, 300, 503 [NASA ADS] [Google Scholar]
- Espaillat, C., Calvet, N., D'Alessio, P., Hernández, J., Qi, C., Hartmann, L., Furlan, E., & Watson, D. M. 2007, ApJ, 670, L135 [NASA ADS] [CrossRef] [Google Scholar]
- Espaillat, C., Calvet, N., Luhman, K. L., Muzerolle, J., & D'Alessio, P. 2008, ApJ, 682, L125 [NASA ADS] [CrossRef] [Google Scholar]
- Fleck, Jr., J. A., & Canfield, E. H. 1984, J. Computational Phys., 54, 508 [NASA ADS] [CrossRef] [Google Scholar]
- Kenyon, S. J., & Hartmann, L. 1996, VizieR Online Data Catalog, 210, 10117 [NASA ADS] [Google Scholar]
- Min, M., Hovenier, J. W., & de Koter, A. 2005, A&A, 432, 909 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Min, M., Dullemond, C. P., Dominik, C., de Koter, A., & Hovenier, J. W. 2009, A&A, 497, 155 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Najita, J. R., Strom, S. E., & Muzerolle, J. 2007, MNRAS, 378, 369 [NASA ADS] [CrossRef] [Google Scholar]
- Piétu, V., Dutrey, A., Guilloteau, S., Chapillon, E., & Pety, J. 2006, A&A, 460, L43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Preibisch, T., Ossenkopf, V., Yorke, H. W., & Henning, T. 1993, A&A, 279, 577 [NASA ADS] [Google Scholar]
- Ratzka, T., Leinert, C., Henning, T., Bouwman, J., Dullemond, C. P., & Jaffe, W. 2007, A&A, 471, 173 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Robitaille, T. P., Whitney, B. A., Indebetouw, R., & Wood, K. 2007, ApJS, 169, 328 [NASA ADS] [CrossRef] [Google Scholar]
- Simon, M., Dutrey, A., & Guilloteau, S. 2000, ApJ, 545, 1034 [NASA ADS] [CrossRef] [Google Scholar]
- Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 [NASA ADS] [CrossRef] [Google Scholar]
- Vinkovic, D., Ivezic, Z., Jurkic, T., & Elitzur, M. 2006, ApJ, 636, 348 [NASA ADS] [CrossRef] [Google Scholar]
- Weaver, W. B., & Jones, G. 1992, ApJS, 78, 239 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Table 1: Stellar and best-fit parameters for the two different geometries.
All Figures
![]() |
Figure 1: Sketches (not to scale) of the two possible geometries for LkCa 15 from this paper ( top panel) and - for comparison - that from Esp07 ( bottom panel). Depicted are: a) A model with a 1 AU optically thick inner disk and an outer disk that rises out of its shadow. b) A model with a 5 AU optically thin spherical dust shell and a flatter outer disk with a fully illuminated inner rim. c) The model proposed by Esp07, a combination of the inner disk from a) with the fully illuminated outer disk of b). |
Open with DEXTER | |
In the text |
![]() |
Figure 2:
Spectral Energy Distribution of LkCa 15 (diamonds) with 1 |
Open with DEXTER | |
In the text |
![]() |
Figure 3:
Radial |
Open with DEXTER | |
In the text |
![]() |
Figure 4: 20 micron images of our best-fit disk models of LkCa 15. The left panel shows a model with an optically thick inner disk of 1 AU that casts a shadow on the rim of the outer disk. The right panel shows a 5 AU optically thin dust shell and a fully illuminated outer rim. Note that if the rim of the outer disk is shadowed, much more radiation emerges from the disk surface. |
Open with DEXTER | |
In the text |
Copyright ESO 2010
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.