Advanced
Dust Around T Tauri Stars
Dust Around T Tauri Stars
Journal of Astronomy and Space Sciences. 2011. Dec, 28(4): 253-260
Copyright ©2011, The Korean Space Science Society
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : September 09, 2011
  • Accepted : November 11, 2011
  • Published : December 15, 2011
Download
PDF
e-PUB
PubReader
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Kyung-Won Suh
kwsuh@chungbuk.ac.kr
Young-Joo Kwon
Abstract
To reproduce the multiple broad peaks and the fine spectral features in the spectral energy distributions (SEDs) of T Tauri stars, we model dust around T Tauri stars using a radiative transfer model for multiple isothermal circumstellar dust shells. We calculate the radiative transfer model SEDs for multiple dust shells using the opacity functions for various dust grains at different temperatures. For six sample stars, we compare the model results with the observed SEDs including the Spitzer spectral data. We present model parameters for the best fit model SEDs that would be helpful to understand the overall structure of dust envelopes around classical T Tauri stars. We find that at least three separate dust components are required to reproduce the observed SEDs. For all the sample stars, an innermost hot (250-550 K) dust component of amorphous (silicate and carbon) and crystalline (corundum for all objects and forsterite for some objects) grains is needed. Crystalline forsterite grains can reproduce many fine spectral features of the sample stars. We find that crystal-line forsterite grains exist in cold regions (80-100 K) as well as in hot inner shells.
Keywords
1. INTRODUCTION
T Tauri stars are generally believed to be low-mass (0.1-2 M ) pre-main-sequence stars of spectral types F to M and surface effective temperatures of 3,000 K to 7,000 K. As a consequence of the star formation process they are surrounded by a gas and dust envelope and/or disk (Bertout 1989). Classical T Tauri stars (CTTS) are optically visible, as they possess a small fraction of their natal en-velope. In the spectral energy distribution (SED) classifi-cation, CTTS are Class II sources and their SEDs peak be-tween 1 and 10 ㎛ with moderate infrared (IR) excesses reaching beyond 100 ㎛.
The Spitzer space telescope (Spitzer) , launched in 2003, has the infrared spectrograph (IRS) with high sensitiv-ity and large spectral range (5-35 ㎛). The high resolu-tion IRS spectroscopic observations have revealed the detailed SEDs of many T Tauri stars in nearby molecular clouds where various amorphous and crystalline dust features are present (Olofsson et al. 2009). The infrared SEDs of the T Tauri stars show multiple broad peaks at 5-200 ㎛ (Sicilia-Aguilar et al. 2009).
In this paper, we model dust around six sample CTTS using the optical properties of various amorphous and crystalline dust grains at different temperatures. We use a radiative transfer model for multiple isothermal circum-stellar dust shells (Towers & Robinson 2009). Even though the model ignores the disk geometry, it is useful to inves-tigate the overall properties of the multiple dust com-ponents around the central star. We compare the model results with the observed SEDs of the stars, including the ground-based, infrared astronomical satellite ( IRAS ), Spitzer and AKARI data. Using the comparison, we deter-mine the model parameters for the best fit model SEDs that would help us to understand the overall structure of dust envelopes around CTTS.
2. SAMPLE STARS
In this paper, we choose six CTTS. For these stars, good quality observational data including the Spitzer IRS spec-tral data in a wide wavelength range are available. The sample stars are listed in Table 1 .
The high resolution Spitzer IRS spectroscopic observa-tions revealed the detailed SEDs of many T Tauri stars in nearby molecular clouds where various amorphous and crystalline dust features are present (Olofsson et al. 2009). The Spitzer has the infrared array camera (IRAC) designed for observations of faint sources and deep large-area sur-veys (Fazio et al. 2004). The IRAC has four channels that can obtain broadband images with high sensitivity at 3.6, 4.5, 5.8 and 8.0 ㎛ simultaneously.
The AKARI (Murakami et al. 2007) made an all sky sur-vey with the infrared camera (IRC) and far infrared sur-veyor (FIS). We can use the AKARI point source catalogue (PSC) data at two bands (9 and 18 ㎛) obtained by the IRC, and the bright source catalogue (BSC) data at four bands (65, 90, 140 and 160 ㎛) obtained by the FIS. For each object, we have cross-identified the AKARI source by finding the nearest one using the position information supplied the Set of Identifications, Measurements, and Bibliography for Astronomical Data (SIMBAD) database.
For each object, Table 1 lists the IRAS PSC number, the reference for the reduced Spitzer IRS data, the cross-iden-tified AKARI PSC number, the AKARI BSC number, and the cloud name with the distance.
To obtain the standard flux in W/m 2 for all the data, we use the zero-magnitude calibrating method. The zero-magnitude calibrating data are taken from the related ref-erences. The observed SEDs of the six stars are displayed in Fig. 1 . For IRAS and AKARI catalogue data, we use only good quality (q = 3) data.
3. DUST ENVELOPE MODEL CALCULATIONS
The infrared SEDs of the CTTS show multiple broad peaks at 5-200 ㎛ (Sicilia-Aguilar et al. 2009). A reason-able explanation for this would be that there are multiple
Sample of classical T Tauri stars.IRAS: infrared astronomical satellite, PSC: point source catalogue, IRS: infrared spectrograph, BSC: bright source catalogue.aFrom Makarov (2007),bfrom Mamajek et al. (1999),cfrom Bertout et al. (1999).
Lager Image
Sample of classical T Tauri stars. IRAS: infrared astronomical satellite, PSC: point source catalogue, IRS: infrared spectrograph, BSC: bright source catalogue. aFrom Makarov (2007), bfrom Mamajek et al. (1999), cfrom Bertout et al. (1999).
dust components radiating multiple regions of infrared wavelengths. It is difficult to reproduce the observed SEDs with a single component dust model (e.g., DUSTY code developed by Ivezić & Elitzur 1997). Many investiga-tors have tried to model the dust around T Tauri stars with a single component dust disk model (Miroshnichenko et al. 1999, Whitney et al. 2003), or using a simple Planck dust radiation law without considering absorption pro-cesses (Bouwman et al. 2008, Olofsson et al. 2010).
Suh (2011) modeled dust around Herbig Ae/Be stars using the radiative transfer model for multiple isothermal spherically symmetric circumstellar dust shells devel-oped by Towers & Robinson (2009), which assumes that each dust shell is in local thermodynamical equilibrium and that the temperature is constant. Dust grains in the innermost shell absorb the radiation from the central star and radiate at the equilibrium temperature. Dust grains in an outer shell absorb the radiation from the central star and the inner shell(s) and radiate at the equilibrium temperature. The scattering of light is ignored. These assumptions would be reasonable approximations for studying dust around T Tauri stars.
In this paper, we use the radiative transfer model for multiple isothermal circumstellar dust shells (Towers & Robinson 2009) to reproduce the observed multiple peaks and crystalline dust features. Though it is a rela-tively simple model, it can use the flexible parameters of dust properties to treat radiative processes through centrally heated multiple dust shells. The code would be useful to investigate the overall properties of complicated distribution of dust around a central star. We model dust envelopes around T Tauri stars using the optical proper-ties of various amorphous and crystalline dust grains at different temperatures.
Molecules in the atmospheres of CTTS make deep ab-sorption features in the near infrared region ( Fig. 1 ). In this paper, we ignore absorption processes by molecules and concentrate on dust around the central star.
For the central star, we assume simple blackbody ra-diation. For each object, we use the best fitting effective temperature and list in Table 2 . For all the sample stars, we assume that the luminosity is the same as the solar
Lager Image
Observations compared with model SEDs for the sample stars. SED: spectral energy distribution IRAS: infrared astronomical satellite PSC: point source catalogue IRS: infrared spectrograph.
luminosity ( L* = 1 L ). A change in the luminosity does not affect the shape of the output SED; it only affects the overall energy output.
- 3.1 Dust Opacity
The best fit model for the observed SED requires a proper combination of multiple isothermal components with different sets of diverse dust opacity functions. We have tried to use as many dust species as possible. We find that four dust species (amorphous silicate, amorphous carbon [AMC], crystalline corundum and crystalline for-sterite) are necessary to reproduce the SEDs of sample stars.
In this paper, we do not consider polycyclic aromatic hydrocarbon for the radiative transfer model calculations, because its thermal properties are not yet well known.
For amorphous silicate, we use the optical constants derived by Suh (1999) for cold silicate. For AMC, we use the optical constants derived by Suh (2000). For the two species, the extinction efficiency factors are calculated for
Lager Image
Dust opacity functions for the four species.
The model parameters of the multiple isothermal dust shells for the best fit model spectral energy distributions.1a: amorphous carbon (Suh 2000), s: amorphous silicate (Suh 1999), C: crystalline corundum (Tamanai et al. 2009), F: crystalline forsterite (Jäger et al. 1998).
Lager Image
The model parameters of the multiple isothermal dust shells for the best fit model spectral energy distributions. 1a: amorphous carbon (Suh 2000), s: amorphous silicate (Suh 1999), C: crystalline corundum (Tamanai et al. 2009), F: crystalline forsterite (Jäger et al. 1998).
spherical dust grains (Bohren & Huffman 1983) from the optical constants given in the references. The radii of the spherical dust grains are assumed to be 0.1 ㎛ uniformly.
For crystalline corundum, we use the extinction data of α-Al 2 O 3 (corundum-sample 1) obtained by Tamanai et al. (2009). For crystalline forsterite, we use the extinction data obtained by Jäger et al. (1998).
Dust opacity functions for the four species are dis-played in Fig. 2 . The crystalline forsterite grains produce conspicuous features at 10.0, 11.2, 16.3 19.5, 23.5, 27.5 and 33.5 ㎛. We had attempted to use other dust species (olivine, ice and oxides), but found they were not useful for this work.
- 3.2 Model SEDs
We perform various radiative transfer model calcu-lations in the wavelength range 0.01 to 36,000 ㎛. We choose 10 ㎛ as the fiducial wavelength that sets the scale of the dust optical depth (τ 10 ). We have computed the model SEDs for various optical depths of the multiple dust shells with different dust opacity.
For each object, we try to find the best fit model SED for the observations. Once we have a set of reason-able model parameters, we compare the model result with the observed SED and repeatedly revise the related parameter(s) until we get a satisfactory fit in the entire wavelength range. We may have to revise the parameters for the central star and the multiple dust shells repeatedly because of the correlated absorption processes.
The model parameters of the multiple isothermal dust shells for the best fit model SEDs are listed in Table 2 . For each object, the parameters of the central star (the black-body temperature) and those for multiple dust shells are listed. For each dust shell, the equilibrium dust tempera-ture, the dust optical depth (τ 10 ), the radius of the shell (r) in the unit of the radius of the central star (R * ) and the dust opacity function are listed. All of the above param-eters are input parameters. For each model, the code cal-culates the radii of the multiple dust shells and the model SED.
4. SED COMPARISON
The six panels of Fig. 1 show the best fit model SEDs compared with the observed SEDs for the six sample stars. For all of the sample stars, we use three or four separate dust shells with the four dust species ( Table 2 ). We could reproduce some of the observed fine spectral features with the theoretical models using the crystalline dust grains. Generally, a cold (24-50 K) outer dust shell with amorphous silicate and AMC grains reproduces the FIR region fairly well.
Haro 1-16 shows very conspicuous amorphous silicate features at 10 and 20 ㎛ because of abundant (68%) sili-cate grains in the innermost shell. The model can repro-duce the evident forsterite dust features at 23.5 and 27.5 ㎛ and the vague feature at 33.5 ㎛. Crystalline corun-dum grains in the hot (550 K) innermost shell improve the fit in the 10-20 ㎛ region.
J0843 shows weak crystalline forsterite dust features and unknown features. The model can reproduce the evi-dent forsterite dust features at 23.5, 27.5 and 33.5 ㎛.
RECX-5 looks to have the hottest (4,000 K) central star. The object shows the most conspicuous crystalline for-sterite dust features in a wide wavelength range. Of the six sample stars, only this object shows the evident 10.0, 11.2 and 16.3 ㎛ features from hot forsterite. The star also shows the conspicuous 23.5, 27.5 and 33.5 ㎛ features. All of those features are reproduced by the theoretical model. The content of forsterite in the warm (250 K) in-nermost shell is very high (20%).
RECX-9 shows many weak spectral features in a wide wavelength range. Some of those features are unknown. The model can reproduce the vague forsterite dust fea-tures at 23.5, 27.5 and 33.5 ㎛. Crystalline corundum grains in inner two shells improve the fit in the 10-30 ㎛ region.
VW Cha shows weak crystalline dust features. The crys-talline forsterite grains that exist only in the cold (90 K) shell can reproduce the features at 23.5, 27.5 and 33.5 ㎛.
WX Cha shows very weak crystalline forsterite features. The crystalline forsterite grains in the two cold (100 and 180 K) shells produce the vague features at 23.5, 27.5 and 33.5 ㎛.
For all the sample stars, crystalline corundum in the hot (250-550 K) innermost shell improves the fit in the 10-20 ㎛ region.
Even though we have tried to reproduce all the features and characteristics of the observed SEDs by the radiative transfer model, we could not reproduce some of them (e.g., some features in J0843 and RECX-9). There are three possible reasons for this. First, we may not have consid-ered some important dust species. Secondly, the dust opacity functions used for this work may need to be im-proved. Finally, the radiative transfer model could be too simple. A more sophisticated radiative transfer model, which can consider multiple components of dusty disk and shell with more dust species, would be able to repro-duce the SEDs better.
The size of dust grains only affects the overall energy output in the far infrared region. We find that changes of the dust size do not cause meaningful differences in the overall fitting for the sample stars.
5. DISCUSSION ON CRYSTALLINE DUST
For all the sample stars, the observed spectral features of crystalline silicate (forsterite) grains are reproduced by the model calculations. Because crystalline silicate grains show very sharp features, even a small content (about 2%) can be easily detectable.
Though a small content (2-3%) of crystalline corun-dum grains does not reproduce their own sharp spectral features directly, those in the hot (250-550 K) innermost shell can improve the fit in the 10-20 ㎛ region for all the sample stars. Similar effects were found for some Herbig Ae/Be stars (Suh 2011).
It is quite evident that crystalline silicate (forsterite) grains exist in cold (48-100 K) outer regions of many T Tauri stars (this paper, Juhasz et al. 2010). Though crys-talline silicates are abundant in many young stellar ob-jects (YSOs) and solar system comets, they are essentially missing from the interstellar medium (Juhász et al. 2010). It would be reasonable to assume that crystallization occurs in the low temperature envelopes (or disks) of YSOs.
Some known processes of crystallization (annealing and direct condensation from the gas phase) require high temperature (about 1,000 K) (Fabian et al. 2000). The vertical mixing by turbulence, in disks, can transfer crystalline grain components, which formed in the high temperatures and high ultraviolet radiation field of the disk atmosphere, into the disk midplane where the tem-perature is very low (Dullemond et al. 2006). Olofsson et al. (2009) argued that the abundant crystalline silicates found far from their presumed formation regions (hot regions) may suggest efficient outward radial transport mechanisms in the disks around T Tauri stars.
On the other hand, Carrez et al. (2002) and Kimura et al. (2008) suggested a mechanism of crystallization at low temperature by reporting that amorphous silicate grains were crystallized to forsterite by electron-beam irradia-tion. YSOs are known to undergo active and frequent flaring events in which electrons are accelerated. Therefore, the electron irradiation of dust in YSO environments could explain the origin of the crystalline silicate grains around T Tauri stars.
6. CONCLUSIONS
To reproduce the multiple broad peaks and fine spec-tral features in the SEDs of CTTS, we have modeled dust around CTTS using a radiative transfer model for mul-tiple isothermal circumstellar dust shells. By comparing the model results with the observed SEDs in a wide wave-length range for the six sample stars, we have presented the model parameters of the best fit model SEDs that will be helpful to understand the overall structure of dust en-velopes around T Tauri stars.
We have found that at least three separate dust compo-nents are required to reproduce the observed SEDs. For all the sample stars, an innermost hot (250-550 K) dust component of amorphous (silicate and carbon) and crys-talline (corundum for all objects and forsterite for some objects) grains is needed. We have found that crystalline forsterite grains can reproduce many fine spectral fea-tures of the sample stars.
A small content of crystalline corundum grains has been found to be present in all of the sample stars. Co-rundum appears to be one of the major dust components in YSOs.
For all the sample stars, the crystalline silicate (forst-erite) grains exist in cold (80-100 K) outer dust shells as well as in hot inner shells. Even though the reason for the existence of low temperature crystalline silicate is still un-certain, the electron irradiation of dust in YSO environ-ments or a process of outward radial transport could be possible scenarios.
Further investigations with sophisticated theoretical models using more various dust species could reveal use-ful information about the physical and chemical proper-ties of dust around T Tauri stars and the general environ-ments of star forming regions.
Acknowledgements
This work was supported by the research grant of the Chungbuk National University in 2010. This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France. This research is based on observations with AKARI, a JAXA project with the participation of ESA.
References
Bertout C 1989 T Tauri stars: wild as dust ARA&A 27 351 - 395
Bertout C , Robichon N , Arenou F 1999 Revisiting Hipparcos data for pre-main sequence stars A&A 352 574 - 586
Bohren CF , Huffman DR 1983 Absorption and scattering of light by small particles Wiley New York
Bouwman J , et al. , Henning Th , et al. , Hillenbrand LA , et al. , Meyer MR , et al. , Pas-cucci I , et al. 2008 The formation and evolution of planetary systems: grain growth and chemical processing of dust in T Tauri systems ApJ et al. 683 479 - 498
Bouwman J , et al. , Lawson WA , et al. , Dominik C , et al. , Feigelson ED , et al. , Hen-ning Th , et al. 2006 Binarity as a key factor in protoplanetary disk evolution: Spitzer disk census of the η chamaele-ontis cluster ApJ et al. 653 L57 - L60
Cambrésy L , et al. , Copet E , et al. , Epchtein N , et al. , de Batz B , et al. , Borsenberger J , et al. 1998 New young stellar object candidates in the Chamaeleon I molecular cloud discovered by DENIS A&A et al. 338 977 - 987
Carrez P , et al. , Demyk K , et al. , Leroux H , et al. , Cordier P , et al. , Jones AP , et al. 2002 Low-temperature crystallization of MgSiO3glasses under elec-tron irradiation: possible implications for silicate dust evolution in circumstellar environments M&PS et al. 37 1615 - 1622
Dullemond CP , Apai D , Walch S 2006 Crystalline silicates as a probe of disk formation history. ApJ 640 L67 - L70
Fabian D , Jäger C , Henning Th , Dorschner J , Mutschke H 2000 Steps toward interstellar silicate mineralogy. V. Thermal evolution of amorphous magnesium silicates and silica A&A 364 282 - 292
Fazio GG , et al. , Hora JL , et al. , Allen LE , et al. , Ashby MLN , et al. , Barmby P , et al. 2004 The infrared array camera (IRAC) for the Spitzer space telescope ApJS et al. 154 10 - 17
Gautier TN III , Rebull LM , Stapelfeldt KR , Mainzer A 2008 Spitzer-MIPS observations of the η Chamaeleontis young as-sociation ApJ 683 813 - 821
Ivezić A , Elitzur M 1997 Self-similarity and scaling behaviour of infrared emission from radiatively heated dust. I. Theo-ry. MNRAS 287 799 - 811
Jäger C , et al. , Molster FJ , et al. , Dorschner J , et al. , Henning Th , et al. , Mutschke H , et al. 1998 Steps toward interstellar silicate mineralogy. IV. The crystalline revolution. A&A et al. 339 904 - 916
Juhász A , et al. , Bouwman J , et al. , Henning Th , et al. , Acke B , et al. , van den Anck-er ME , et al. 2010 Dust Evolution in Protoplanetary Disks Around Herbig Ae/Be Stars- the Spitzer View. ApJ et al. 721 431 - 455
Kimura Y , Miyazaki Y , Kumamoto A , Saito M , Kaito C 2008 Char-acteristic low-temperature crystallization of amor-phous Mg-bearing silicate grains under electron ir-radiation. ApJ 680 L89 - L92
Lawson WA , Crause LA , Mamajek EE , Feigelson ED 2001 The η Chamaeleontis cluster: photometric study of the RO-SAT-detected weak-lined T Tauri stars. MNRAS 321 57 - 66
Lawson WA , Feigelson ED , Huenemoerder DP 1996 An improved HR diagram for Chamaeleon I pre-main-sequence stars. MNRAS 280 1071 - 1088
Luhman KL , et al. , Allen LE , et al. , Allen PR , et al. , Gutermuth RA , et al. , Hartmann L , et al. 2008 The disk population of the Chamaeleon I star-forming region ApJ et al. 675 1375 - 1406
Lyo A , Lawson WA , Bessell MS 2004 The spectroscopic character-istics of intermediate aged pre-main-sequence stars: the η Chamaeleontis cluster. MNRAS 355 363 - 373
Makarov VV 2007 Signatures of dynamical star formation in the ophiuchus association of pre-main-sequence stars. ApJ 670 1225 - 1233
Mamajek EE , Lawson WA , Feigelson ED 1999 The η Chamaeleon-tis cluster: a remarkable new nearby young open cluster. ApJ 516 L77 - L80
Megeath ST , Hartmann L , Luhman KL , Fazi GG 2005 Spitzer/IRAC photometry of the η Chameleontis association. ApJ 634 L113 - L116
Miroshnichenko A , Ivezić Z , Vinković D , Elitzur M 1999 Dust emission from Herbig Ae/Be stars: evidence for disks and envelopes. ApJ 520 L115 - L118
Mundt R , Bastian U 1980 UBV photometry of young emission-line objects. A&AS 39 245 - 250
Murakami H , et al. , Baba H , et al. , Barthel P , et al. , Clements DL , et al. , Cohen M , et al. 2007 The infrared astronomical mission AKARI. PASJ et al. 59 S369 - S376
Olofsson J , et al. , Augereau JC , et al. , van Dishoeck EF , et al. , Merin B , et al. , Lahuis F , et al. 2009 C2D Spitzer-IRS spectra of disks around T Tauri stars IV. Crystalline silicates. A&A et al. 507 327 - 345
Olofsson J , et al. , Augereau JC , et al. , van Dishoeck EF , et al. , Merin B , et al. , Grosso N , et al. 2010 C2D Spitzer-IRS spectra of disks around T Tauri stars. V. Spectral decomposition A&A et al. 520 A39 -
Shevchenko VS , Herbst W 1998 The search for rotational modula-tion of T Tauri stars in the ophiuchus dark cloud. AJ 116 1419 - 1431
Sicilia-Aguilar A , et al. , Bouwman J , et al. , Juhasz A , et al. , Henning Th , et al. , Roc-catagliata V , et al. 2009 The long-lived disks in the η Chamae-leontis cluster. ApJ et al. 701 1188 - 1203
Simon M , et al. , Ghez AM , et al. , Leinert Ch , et al. , Cassar L , et al. , Chen WP , et al. 1995 A lu-nar occultation and direct imaging survey of multiplicity in the Ophiuchus and Taurus starforming regions. ApJ et al. 443 625 - 637
Spangler C , Sargent AI , Silverstone MD , Becklin EE , Zuck-erma B 2001 Dusty debris around solar-type stars: temporal disk evolution. ApJ 555 932 - 944
Suh KW 2011 Dust around Herbig AE/Be stars JKAS 44 13 - 21
Suh KW 2000 Optical properties of the carbon dust grains in the en-velopes around AGB stars. MNRAS 315 740 - 750
Suh KW 1999 Optical properties of the silicate dust grains in the en-velopes around AGB stars. MNRAS 304 389 - 405
Tamanai A , et al. , Mutschke H , et al. , Blum J , et al. , Posch Th , et al. , Koike C , et al. 2009 Morphological effects on IR band profiles. Experimen-tal spectroscopic analysis with application to observed spectra of oxygen-rich AGB stars. A&A et al. 501 251 - 267
Towers IN , Robinson G 2009 A model for multiple isothermal cir-cumstellar dust shells. PhyS 80 015901 -
Whitney BA , Wood K , Bjorkman JE , Wolff MJ 2003 Two-dimen-sional radiative transfer in protostellar envelopes. I. Ef-fects of geometry on Class I sources. ApJ 591 1049 - 1063
Zacharias N , et al. , Monet DG , et al. , Levine SE , et al. , Urban SE , et al. , Gaume R , et al. The Naval Observatory Merged Astrometric Data-set (NOMAD). in American Astronomical Society 205th Meeting. San Diego CA Jan 2004 #48.15. et al. 9 - 13