Young stellar objects (YSOs) are surrounded by a gas and dust disk as a consequence of the star formation process (e.g., Hartmann 2009). The T Tauri stars are generally known to be low-mass (0.1-2 Msun) pre-main sequence (PMS) stars. The Classical T Tauri (CTT) stars are optically visible because they have relatively thinner dust envelopes (e.g., Bertout 1989). The Herbig Ae/Be (HAB) stars are generally believed to be intermediate-mass (2-10 Msun) PMS stars (e.g., Waters & Waelkens 1998).
To understand the physical and chemical properties of the circumstellar disk, high-resolution infrared (IR) spectral observations are required. The infrared space observatory (ISO) and Spitzer space telescope (Spitzer) have provided high resolution spectroscopic observational data revealing detailed spectral energy distributions (SEDs) of many YSOs, which show various amorphous and crystalline dust features (Meeus et al. 2001; Olofsson et al. 2009; Sicilia-Aguilar et al. 2009).
In this study, we present theoretical models for the circumstellar dust disks around eight sample YSOs (four CTT stars and four HAB stars) using the opacity functions of various amorphous and crystalline dust species. We use the radiative transfer model, CGPLUS, for the irradiated circumstellar dust disks developed by Dullemond et al. (2001). Even though this model uses a single dust component of the disk geometry ignoring dust shell components, it is useful for investigating the overall properties of the disk around the central star. We compare the model SEDs with the observed SEDs for the sample stars and obtain the theoretical model parameters that reproduce the observed SEDs.
In two previous studies (Suh 2011; Suh & Kwon 2011; hereafter SK2011), a radiative transfer model for multiple isothermal circumstellar dust shells was used to reproduce the observed SEDs for eight CTT stars and eight HAB stars. In this work, we choose four CTT stars and four HAB stars that can be more easily fitted by the model for a single dust disk, which is used in this work (see section 3). This better fit could be because the contribution from the dust shell components is relatively small compared to the contribution from the dust disk for the new eight sample stars.
For the eight sample stars listed in Table 1, good quality observational data, including the Spitzer or ISO spectral data in a wide wavelength range, are available. We use the Spitzer infrared spectrograph (IRS) data for the CTT stars (Olofsson et al. 2009; Sicilia-Aguilar et al. 2009). For the HAB stars, we use the ISO short wavelength spectrometer (SWS) (λ = 2.4-45.2 μm) and long wavelength spectrometer (LWS) (λ = 43-197 mm) data reduced by Meeus et al. (2001).
[Table 1.] Sample of the YSOs and the model parameters for the central star
Sample of the YSOs and the model parameters for the central star
For the sample stars, we use the data from the AKARI (Murakami et al. 2007) point source catalog (PSC) and Bright Source Catalog (BSC) data in six bands (9, 18, 65, 90, 140, and 160 μm). We use the PSC data from the two micron all sky survey (2MASS) (Skrutskie et al. 2006) in 1.25, 1.65, and 2.17 μm bands. We use the wide-field infrared survey explorer (WISE) (Wright et al. 2010) data in four bands (3.4, 4.6, 12, and 22 μm). We also use the ground-based infrared observational data obtained from various authors.
For each object, Table 1 lists the infrared astronomical satellite (IRAS) PSC, AKARI PSC, AKARI BSC numbers, and the model parameters for the central star (see section 3.3). In Figs. 1 and 2, we present the observed SEDs of the eight sample stars compared with the theoretical model SEDs (see section 3).
3. DUST DISK MODEL CALCULATIONS
The observed SEDs of YSOs were reproduced by using various models. Using a simple Planck dust radiation law, Bouwman et al. (2008) and Olofsson et al. (2010) tried to reproduce the observed SEDs of the YSOs. To overcome the limitation of the single component dust shell model (e.g., DUSTY developed by Ivezić & Elitzur 1997), SK2011 used a radiative transfer model for the multiple isothermal circumstellar dust shells to reproduce the observed SEDs of the YSOs. However, the dust envelopes of the YSOs are believed to have a disk geometry (e.g., Whitney et al. 2003).
3.1 The Radiative Transfer Model
In this work, we perform the radiative transfer model calculations using the CGPLUS model developed by Dullemond et al. (2001) for a passive irradiated circumstellar dust disk to reproduce the observed SEDs of the sample YSOs. The model is based on the flaring disk model given by Chiang & Goldreich (1997). However, the central regions of the disk are removed, and the inner rim of the disk is puffed up and is much hotter than the rest of the disk because it is directly exposed to the stellar flux. The circumstellar dust disk with an inner hole is irradiated and flared by the central star. The model can use various parameters of the disk geometry to consider radiative processes through the dust disk, but does not consider the scattering processes by dust or absorption processes by molecules. We model the dust disks around the YSOs using the opacity functions of various dust species (see section 3.2).
The dust density distribution uses a simple power law with a power law index (
The disk has six adjustable input parameters: the outer radius (
[Table 2.] Model parameters of the dust disk for the best fit model
Model parameters of the dust disk for the best fit model
For the central hot source, we assume that the object is a single star that emits simple blackbody radiation. The model parameters for the central star are the mass, luminosity, blackbody temperature, and distance (
Dust grains in the disk of the YSOs absorb and scatter the stellar radiation and re-emit the radiation at longer wavelengths. In this work, we use the seven dust species listed in Table 3. We do not consider the polycyclic aromatic hydrocarbons (PAH).
[Table 3.] Dust opacity functions used in this work
Dust opacity functions used in this work
For amorphous dust, we use two dust species: amorphous silicate and amorphous carbon (AMC). We use the optical constants derived by Suh (1999) for amorphous warm silicate. We use the optical constants derived by Suh (2000) for AMC. For crystalline olivine and water-ice, we use the optical constants obtained by Fabian et al. (2001) and those obtained by Bertie et al. (1969), respectively. For crystalline corundum, we use the optical constants of
For crystalline forsterite and enstatite, the MAC data obtained by Jäger et al. (1998) are used. Dust opacity functions for the seven species are displayed in Fig. 3. The crystalline dust grains show conspicuous emission features.
We perform model calculations (see section 3.1) in the wavelength range from 0.1 to 200 μm. We find a set of model parameters that produces a similar fit to the observed SED and revise the related parameter (
The best fit model SEDs are presented in Figs. 1 and 2. We list the model parameters of the central star and dust disk for the best fit models in Tables 1 and 2. Table 1 lists the model parameters of the central star: the mass (
Figs. 1 and 2 display the best fit theoretical model SEDs as well as the observed SEDs for the four CTT and four HAB sample stars. The theoretical model reproduces the observed SEDs, including the fine spectral features using crystalline dust as well as amorphous dust in the disk for most of the sample stars. The model SEDs, which do not use the crystalline dust, are also shown for comparison. When crystalline dust is not used,
Haro 1-16 appears to have very abundant (90%) silicate dust, and it shows prominent 10 and 20 μm features. The object seems to have a relatively hotter (3,500K) central star compared with the other CTT objects. The forsterite features at 23.5, 27.5, and 33.5 μm are reproduced by the model. Because of the limitation of the theoretical model, which uses a single disk component, the model fit in some wavelength ranges could not be improved.
For VW Cha, the crystalline forsterite features are reproduced by the best fit model, which requires the highest AMC content (55%). This object appears to have the smallest disk mass (3×10−4
For WX Cha, the crystalline features could be due to various crystalline dust species (forsterite, enstatite, and corundum). The model requires the largest disk mass (1×10−2
J0843 appears to be the nearest (D = 97 pc) object for which the model requires the lowest stellar luminosity (0.13
For AB Aur, the model uses only amorphous dust (silicate and AMC). Unlike the other objects, the crystalline dust species are not useful for reproducing the observed SED. The best fit model requires the largest stellar luminosity (47
HD 100546, which has the largest crystalline dust content (26%), shows emission features due to various crystalline dust species (forsterite, enstatite, water-ice, and corundum). We were not able to reproduce the wavelength range 3-8 μm, which is dominated by the PAH features. The object shows prominent crystalline forsterite features at 19.5, 23.5, 27.5, and 33.5 μm, but the best fit model produces much more prominent features. This could be because the forsterite grains exist in a separate cold outer dust component (about 100K) (Suh 2011), but the CGPLUS model can consider only a single disk component. This object requires the largest disk size (
For HD 144432, the model uses various crystalline dust species (forsterite, enstatite, and corundum) as well as silicate (80%) and AMC. The crystalline dust makes a minor improvement in the model fit.
For HD 150193, the best fit model uses crystalline olivine (10%) as well as silicate (85%) and AMC (5%). The crystalline olivine produces a slightly better fit in the 9-30 μm range.
We find that the theoretical model, CGPLUS, is useful for investigating the overall physical properties of the disk around the central star, even though it uses a single dust component of the disk geometry ignoring other dust disk or shell components (see Tables 1 and 2). The model parameters for the central star for the TTS stars are:
The model parameters for the dust disk for TTS stars are:
It is possible that the variety of the model parameters of the dust disk is not due to the different intrinsic property of the single central star, whether it is a CTT or a HAB star. Even though we have assumed that the central object is a single star for this work, it is believed that a major portion of the YSOs are binary and multiple systems, and these systems have significant dynamical effects on the structure of the dust disk (e.g., Bouwman et al. 2006; Hartmann 2009). The disk shapes for the sample YSOs could be significantly influenced by the binarity of the central object.
There could be three possible reasons for the unavoidable discrepancies between the theoretical and the observed SEDs. First, we may have not used some dust species that are present in the disks of the YSOs. Second, the dust opacity functions used in this work may need to be improved. Finally, the model for a single dust disk component used for this work could be too simple to reproduce all of the observed features.
Because emission features of crystalline dust grains are conspicuous, a small content (about 5-26%) can be easily observed (see Figs. 1 and 2). The crystalline grains can also improve the model fit in a wide spectral range for most objects. It is known that crystalline silicates are essentially missing from the interstellar medium, but they are abundant in many YSOs and solar system comets (e.g., Juhász et al. 2010; Suh 2014). Though the annealing crystallization mechanism suggested by Fabian et al. (2000) requires a high temperature (about 1,000K), other mechanisms that do not require a high temperature have also been suggested by various authors (SK2011). For most sample stars, we find that the dust temperature is low for a significant portion of the crystalline grains because the temperature of the outer regions of the disk is quite low (
SK2011 used a radiative transfer model for multiple isothermal circumstellar dust shells to reproduce the SEDs of the YSOs. Generally, the multiple dust shell model reproduced the observed SEDs of the YSOs better than the single dust disk model used for this work. Because the basic schemes of the two models (SK2011 and this work) are totally different, the dust compositions derived from the two models cannot be directly compared. Furthermore, SK2011 did not use the crystalline olivine dust that produces peaks at some wavelengths similar to those for crystalline corundum or forsterite (see Fig. 3). Therefore, even though the overall dust compositions appear fairly similar, the detailed dust compositions derived from the two models show some discrepancies.
The isothermal multiple dust shell model used by SK2011 would be more useful for identifying the dust composition because it reproduces the observed SEDs much better without providing detailed physical properties (density distributions) of the dust envelopes. In contrast, the single dust disk model (CGPLUS) used in this work makes reasonable estimations of the physical properties (the detailed density distribution and mass) of the dust disk. The model parameters derived in this work would be useful for preparing initial model parameters of a more complicated radiative transfer model, which can thoroughly consider multiple dust disk and shell components (e.g., Suh 2016).
To reproduce the SEDs of the YSOs, we have used a radiative transfer model for the irradiated circumstellar dust disk with various shapes and many dust species. For eight sample objects of the CTT and HAB stars, we have compared the theoretical model results with the observed SEDs described by the Spitzer and ISO spectral data. We have found that crystalline dust as well as amorphous dust in the disk can reproduce the observed SEDs fairly well for most of the sample stars.
We have presented the best fit model parameters of the dust disk for the eight sample stars after comparing the theoretical model SEDs with the observed SEDs. We have found that the model requires a higher mass, luminosity, and temperature for the central star for the HAB stars than for the CTT stars. Generally, the outer radius, the total mass, the thickness, and the rim size of the dust disk for the HAB stars are larger than those for the CTT stars. These disk model parameters could be significantly influenced by the binarity of the central object.
Because a major portion of the YSOs are believed to be binary and multiple systems, we may need to consider their dynamic effects in modeling the circumstellar dust envelopes around the YSOs. We expect that the model parameters of the dust disks obtained in this work would be useful for further investigations using a radiative transfer model that can consider multiple dust disk and shell components.