1 Introduction
In the early 1980s, optical polarimetric observations [13] revealed the existence of elongated and flattened circumstellar dust material around some Pre-Main-Sequence (PMS) stars such as the low-mass TTauri stars, similar to the Sun when it was about years old. A few years later, observations from the InfraRed Astronomical Satellite (IRAS) found significant infrared (IR) excesses around many TTauri stars [24], showing the existence of cold circumstellar dust around these objects. These exciting discoveries motivated several groups to model the Spectral Energy Distribution (SED) of these young stars and to propose an evolutionary scenario [1].
IR excess around TTauri stars was integrated as the signature of a dust disk orbiting the central stars, the material being the residue of the molecular cloud which formed the central star. With total masses () ranging from 0.001 to a 0.1 M⊙, most of these disks contain enough gas (H2) to allow, in theory, the formation of a giant planet. Hence, they are often called ‘protoplanetary disks’. During this phase, the dust emission is optically thick in the Near-IR (NIR) and the central young star still accretes material from its disk. Such disks are then also called ‘accretion disks’.
Proto-planetary disks orbiting around young stars of a few million years old are now routinely observed by large single-aperture telescopes and interferometers from the radio up to the optical domain. Large optical telescopes provide images of the disk surface where the dust grains scatter the impinging stellar light. When the dust disk is observed close to edge-on, these images reveal that the disks are flaring, as expected from the theory of disks in hydrostatic equilibrium [23]. In the meantime, spectro-imaging of the CO gas provided by millimetre/submillimetre (mm/submm) arrays have shown that protoplanetary disks are in Keplerian rotation around their central star and that disks can be as large as 500–800 AU, in radius. At the distance of the closest star forming regions ( pc), the angular radii of these disks are usually within the range 1– or 150–450 AU. However, being mostly heated by the central star, their spectral energy distribution is such that the apparent angular size drops very fast when the observing wavelength goes from the millimeter domain to the optical one. Typically the apparent radius of the thermal emission of a hot dust disk at 10 m would not exceed a few AU or while it would be around 1– at mm.
Fig. 1 is a sketch showing a ‘classical’ accretion disk orbiting a TTauri star of a few million years and located at the Taurus distance (150 pc). From the observations, one can then naturally distinguish three parts in disks. The very inner disk, located close to its inner radius, is extending up to radius of –3 AU and is only resolved by optical and NIR interferometry. The inner disk is extending up to a radius of AU and would somewhat correspond in our system to the inner solar system, located up to the first Kuiper Belt objects. The outer disk corresponds to the relatively cold (–50 K) large dust and gas disk usually traced by CO rotation line emission and observed by millimeter and sub-millimeter interferometers such as the IRAM array (Fig. 2), the SMA or OVRO.
More generally, these last years have revealed the potential of interferometric techniques as powerful tools to resolve and model such disks. This is true both for small hot inner disks and large colder outer disks. It appears that only aperture synthesis (interferometry) can nowadays technically achieve the angular resolution needed.
2 Disk properties: from large to smaller scales
The next section presents some major highlights resulting from the technical improvements of large telescopes and interferometers.
2.1 Gas disk observed with millimeter/submillimeter arrays
Current mm/submm interferometers provide images of the dust and gas outer disks at an angular resolution of or 150 AU at the distance of the nearest star forming regions. Molecular hydrogen is, by far, the major constituent of the disk (with a dust-to-gas ratio taken as ) but it remains difficult to observe essentially because it does not possess a dipolar rotational spectrum. Its quadrupolar rotational lines, in the MID-IR domain, mostly trace material at hot/warm temperatures in the disk surface close to the star and do not sample the bulk of the cold gas located along the mid-plane. After H2, carbon monoxide is the most abundant molecule (with abundance CO/H2 of order ). Its first rotation lines are observable with current mm arrays allowing to trace the outer gas disk properties. Since the density in disks is very high ( cm−3), the -0 and -1 CO lines are thermalized by collision with H2 in almost the whole disk. Hence, a simple model of a Keplerian disk assuming LTE conditions is sufficient to retrieve the CO disk properties [9].
CO maps not only reveal that disks are in Keplerian rotation around TTauri [18] and Herbig Ae stars [20], they also allow us to retrieve the physical conditions in outer disks. Comparing resolved CO maps to a disk model by performing a minimization of the disk parameters gives quantitative information on the density and temperature radial and vertical distributions, provided the analysis properly takes into account the transfer function of the interferometer [11,14].
The radial profiles of the temperature deduced from12CO images appear to be consistent with models of stellar heating in flared disks, the outer disk being colder than the inner one. Moreover, since the and -0 and -1 lines have different opacities, they sample different disk layers and allow us to probe the existence of a vertical temperature gradient. Dartois et al. [8] have shown that the ‘CO disk surface’, traced by , is significantly hotter than the mid-plane traced by the optically thinner -0 transition, in the DM Tau disk. This effect appears in the region of the disk where the dust is still optically thick in absorption to the stellar radiation while it is already optically thin to its own emission, around –200 AU in the DM Tau case. Beyond AU where the dust becomes optically thin to both processes, the temperature profile appears vertically isothermal. These observations were the first direct confirmation of models of the vertical thermal structure in disks (e.g., [5–7]). More recently, Piétu et al. [22] have extended this type of analysis to several objects, confirming the existence of vertical and radial temperature gradients. They also find that the turbulence appears to be subsonic in the outer part of disks [22].
It is also important to mention that CO in the gas phase has been observed in all TTauri disks at a temperature which appears to be below the CO freeze out point (17 K). This suggests that the vertical turbulence should play an important role by providing chemical mixing between the vertical molecular layers.
Searching for molecules other than CO remains sensitivity limited to the most abundant molecules found in molecular clouds. So far, there are only a few attempts to survey a large set of molecules in protoplanetary disks [10,17,27,11]. In addition to and C18O, only HCO+, H , DCO+, CS, HCN, HNC, DCN, CN, H2CO, N2H+and C2H have been detected. This domain will also benefit from ALMA, the large submm interferometer currently in construction in Chili.
2.2 MID interferometry: mineralogy of the inner disk
Continuum observations of the dust, both at mm and NIR wavelengths, imply that the dust in such objects has evolved compared to that found in the interstellar medium. In particular, the particle size has started to increase from sub-micron size up to a few centimeters [25]. One key problem in planetary formation is the understanding of the dust evolution in term of grain size, vertical and radial distribution, composition and chemical nature from the inner disk to the outer disk. For example, in the interstellar medium (ISM), the silicate grains are amorphous while in our Solar system, comets and meteorites present crystalline silicates. The new generation of optical, Near-IR and MID-IR interferometers begin to provide direct insights on the inner disks where planetary formation is thought to occur. The observations discussed below present a qualitative but major improvement in the knowledge of the dust properties in inner disks.
Van Boekel et al. [26] recently observed three disks orbiting young stars of M⊙ and a few million year old with MIDI, the 10 m interferometer operated on the Very Large Telescope Interferometer (VLTI) by ESO. They compare the spectra of the silicate bands at 10 m of the very inner disks (radius –4 AU) observed with the VLTI with the spectra of the disks outside –4 AU. The comparison reveals that in the inner disks the silicate grains appear more crystallized (1) than any other dust observed around young stellar objects and (2) than in the outer part of the disks. Fig. 3 shows the spectra. In the very inner disks, the spectra present the three peaks which are characteristics of crystalline silicates; the same peaks are also observed in the spectra of comets such as Halley or Hale-Bopp. The spectra in the outer part of the disks are more similar to those observed in the ISM, characteristics of amorphous silicates. The presence of crystalline silicates in inner disks of a few million years old suggests that silicates crystallize very early in the evolution of the disks, likely before the formation of terrestrial planets.
So far, these recent data bring the first observational constraints on the properties of dust grains at the scale of a few AU in disks. These kinds of observations will likely provide in the next years the first detailed analysis of the dust mineralogy in inner disks, leading to a better understanding of the dust evolution in the earliest phases of the planetary formation.
More generally, recent progress in observational astrophysics begins to provide the sensitivity and the angular resolution needed to study the regions of proto-planetary disks where planetary formation should occur.
3 Dissipation of dust and gas disks
We have seen in the previous section that the dust found around disks a few million years old has strongly evolved compared to that found in the ISM and in star forming regions. The IR excess is characteristic of the amount of dust found in inner disks; the observation of this excess is often used as a tracer of the inner disk evolution while the submm/mm flux can be reasonably considered as a tracer of the dust in the outer disk. Tracing the evolution of the gas (the bulk of the disk mass) is difficult because it is essentially based on the observation of the low J CO rotational lines. These CO lines also strongly emit in molecular clouds where most of these stars are, at least partly, embedded. As a consequence, the high opacity of the CO lines in the clouds does not allow astronomers to detect CO disks, except if the disks are located at the periphery. It is then difficult to get significant statistics on the gas component evolution. This problem will likely disappear in a few years as soon as the submm array ALMA provides enough sensitivity to systematically survey the rarer isotopes of CO.
3.1 The lifetime of disks
Observations of numerous TTauri stars located in star forming regions of various ages have been made by several groups. Fig. 4 is a compilation by [15] of the IR excess (traced by the H–K excess sampling the very inner disk) versus the age of the stars. There are many uncertainties associated to this kind of analysis, the most important being (1) the errors on the estimate of the IR excess and (2) the derivation of the stellar age (which is taken as the median apparent age of the individual objects belonging to the same cluster). However, there is a clear trend which can be seen on this plot: at age around 1 million years, 80–90% of the TTauri stars are surrounded by an innermost dust disk while the fraction of inner disks drops to 30–40% at 2–3 million years and beyond 5 million years, almost all disks have disappeared.
Comparison of the IR excess with the observations of the dust observed in the mm/submm domain [3] and CO data [12] suggest that the dissipation of the inner and outer disks are roughly simultaneous. However, since the statistics in the mm/submm domain are poor, it is difficult to derive quantitative numbers, yet.
Nevertheless, taking into account the various types of surveys and their results, it seems reasonable to conclude that the dissipation is globally a fast process since more than half of the inner disks disappear after 2–3 Myr.
3.2 Dissipation of disks
There are many physical reasons which should explain the apparent dispersal of disks:
- • the formation of a planetary system;
- • the formation of a binary or a multiple stellar system;
- • the tidal truncation by stellar encounters;
- • the photo-evaporation of the disk.
Large optical/NIR surveys of star forming clusters have shown that many stars form in binary and/or multiple systems (e.g., [19]). The formation of a binary or a multiple stellar system strongly affects the original distribution of material which is tidally disrupted. For example, observations of the TTauri binary system GG Tau [9] have shown that the gas and the dust are located inside a circumbinary disk orbiting the binary while very small individual disks exist around each stellar component, inside the Roche lobe. Hence, even if the innermost disks can exist, statistics provided by studies such as those presented above, may be partially biased by the presence in the samples of unresolved multiple stellar systems which are not yet known.
In some clusters, such as Ophiuchus or Orion, the stellar density can be as high as 2– pc−3 (eg [16] and external tidal truncations of disks by stellar encounters should be a relatively common phenomenon. For a matter of sensitivity and angular resolution, we do not have yet enough resolved images of disks in these clusters to estimate the importance of this effect. It is, however, reasonable to take it seriously into account, at least as a phenomenon which should be responsible for the existence of smaller disks in dense stellar clusters.
The photo-evaporation of disks is invoked as a possible phenomenon responsible for the disruption of inner disks [2]. UV photons coming from the young stars reach the disk surface where they can create a ionization front, heating the material behind. The ionized gas becomes gravitationally unbound and flows in a wind away from the disk surface. This phenomenon is observed for the proplyds orbiting TTauri stars located in the Orion Trapezium region [4], close to the OB association. The strong ionizing UV flux from the OB stars photo-evaporates the outer disks in their surroundings. Alexander and collaborators have shown that in some cases photo-evaporation should clear the inner disk very efficiently, creating a central hole, and disrupt the whole disk on a time-scale of about years after a disk lifetime of 2–3 Myr. This physical process has then to be taken into account because the time-scales (both lifetime and dispersal) are those which are observed for TTauri disks, as illustrated in the Section 3.1. However, there is yet no clear observational evidences of the importance of this process for the disk clearing. Dedicated observations, such as sensitive and resolved mm/submm images of inner disks, have to be performed in the next years in order to disentangle between the various scenarii which can be responsible for the existence of an inner hole in disks: photo-evaporation or planetary formation.
Moreover, whatever the frequencies of the various physical processes described in this section are, all of them may affect planet formation since they have a direct impact on the dust and gas distribution close to the stars with the same time-scales.
4 Towards planetary formation: the LkCa15 case
The extension of the IRAM array baselines performed in summer 2006 provides an angular resolution of which corresponds to AU at the Taurus distance; such a linear resolution is of the order of our Solar System size. Using this new configuration of the IRAM array, Piétu et al. [21] have observed the disks surrounding MWC480 and LkCa15, two young stars of a few million years old and masses of and 1 M⊙, respectively. Fig. 5 presents the images of the thermal emission of the dust at 2.8 and 1.3 mm around MWC480 and LkCa15. The emission from the dust disk of MWC480 is strongly centrally peaked, as expected from a standard density distribution, while the central part of the LkCa15 disk appears partly devoid of dust. The best modelling of the data confirms the visual impression in both cases. The dust disk of LkCa15 exhibits a central cavity of radius AU. Such a cavity is not detected at NIR wavelengths because, contrary to the mm/submm range, the high opacity of the dust in the NIR precludes the detection of cavities corresponding to low or moderate density contrast (–100). The most probable explanation for such a hole is the presence of a massive planet of 5–10 masses of Jupiter orbiting at AU from the star or a low-mass companion ( M⊙).
The LkCa15 disk has been observed at mm wavelengths with lower angular resolutions (ranging from to 200 to 70 AU) by several interferometers. None of these observations were suggesting such a wide cavity. This example illustrates how an increase in linear resolution (and sensitivity) can change the physics of an astrophysical object, as soon as the spatial scales needed to study the physical processes are reached.
5 Summary and open questions
In the last fifteen years, resolved observations, from the mm to the optical domain, of the circumstellar disks orbiting TTauri stars have opened the process leading to the formation of planetary systems to direct investigations. Today, we are far from understanding the whole process but the characterization of the physical and chemical properties of these proto-planetary disks represents a major step since planets are built from the disk material. Our current knowledge allows us to tell that:
- • many disks can have large outer radii, within the range 100–800 AU;
- • disks are usually in Keplerian rotation;
- • optical images reveal that disks are flaring;
- • vertical temperature gradients in proto-planetary disks appear to be compatible with heating by the central star, in agreement with disk models;
- • both mm and optical data show that grains in disks are more evolved than in the parent molecular clouds, with sizes up to a few centimeters;
- • NIR surveys suggest that inner disk typical lifetimes are 2–3 million years, beyond 5 million years, most proto-planetary disks may have disappeared.
We have seen through several examples that multi-wavelength approaches are required to get a global understanding of the physical properties of these disks. Optical to MID-IR interferometry has started to unveil the physics of the inner disk. The inner disk will be imaged by the next generation of instruments. With its capability of reaching AU-scale angular resolution at wavelengths where dust has a moderate opacity, ALMA is expected to be the premium instrument to probe these regions and allow astronomers to constrain the models of planet formation. ALMA will also have the potential of providing much higher angular resolution images of the emission of molecular lines allowing detailed modelling of the gas content and of the disk kinematics.
The forthcoming generation of astronomical facilities will likely change our views on planetary formation. In particular, they will bring new information on the following fundamental questions:
- • what is the impact of the parent cloud and its vicinity on the disk properties and its evolution?
- • what is the feedback of the young star on the gas and dust disk and its evolution (photo-evaporation - photo-chemistry at the disk surface) and how does it affect planetary formation?
- • what are the time-scales for grain growth and for vertical settlement of the dust on the disk mid-plane?
- • more generally, what is the mass distribution (gas + dust) radially and vertically in proto-planetary disks and how does it evolve with time?
- • what are the indirect/direct evidences of the formation of a planet in a proto-planetary disk?
Finally, the comparison of these proto-planetary disks which are expected to be representative of the first phases of a planetary system with our own Solar system remains fundamental, as a guideline for each other domain. In the coming years, the understanding of the dynamical and physical/chemical evolution of the dust and gas content in disks will naturally provide a better estimate of the chronology of the planetary formation and a better understanding of the formation of our own Solar System.
Acknowledgements
Anne Dutrey thanks her collaborators: S. Guilloteau, M. Simon, E. Dartois and V. Piétu for a fruitful long-term collaboration. Drs L. Hillenbrand and R. Van Boekel are acknowledged for providing figures.