Issue |
A&A
Volume 497, Number 2, April II 2009
|
|
---|---|---|
Page(s) | 409 - 421 | |
Section | Interstellar and circumstellar matter | |
DOI | https://doi.org/10.1051/0004-6361/200811018 | |
Published online | 24 February 2009 |
Long-wavelength observations of debris discs around sun-like stars![[*]](/icons/foot_motif.png)
V. Roccatagliata1 - Th. Henning1 - S. Wolf2,1 - J. Rodmann3 - S. Corder4 - J. M. Carpenter4 - M. R. Meyer5 - D. Dowell6
1 - Max-Planck-Institut für Astronomie (MPIA), Königstuhl 17, 69117 Heidelberg, Germany
2 -
University of Kiel, Institute of Theoretical Physics and Astrophysics, Leibnizstrasse 15, 24098 Kiel, Germany
3 -
Research and Scientific Support Department, ESA/ESTEC, 2201 AZ Noordwijk, The Netherlands
4 -
California Institute of Technology, Department of Astronomy, MS 105-24, Pasadena, CA 91125, USA
5 -
Steward Observatory, The University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721, USA
6 -
Jet Propulsion Laboratory, California Institute of Technology, Mail Stop 169-506, Pasadena, CA 91109, USA
Received 23 September 2008 / Accepted 2 February 2009
Abstract
Context. Tracing the evolution of debris discs is essential for our understanding of the architecture of planetary system. Even though the evolution of their inner discs has been recently studied with the Spitzer Space Telescope at mid- to far-infrared wavelengths, the outer discs are best characterised by sensitive millimetre observations.
Aims. The goal of our study is to understand the evolution timescale of circumstellar debris discs, and the physical mechanisms responsible for such evolution around solar-type stars. In addition, we make a detailed characterisation of the detected debris discs.
Methods. Two deep surveys of circumstellar discs around solar-type stars at different ages were carried out at 350 m with the CSO and at 1.2 mm with the IRAM 30-m telescope. The dust disc masses were computed from the millimetre emission, where the discs are optically thin. Theoretically, the mass of the disc is expected to decrease with time. To test this hypothesis, we performed the generalised Kendall's tau correlation and three different two-sample tests. A characterisation of the detected debris discs has been obtained by computing the collision and Poynting-Robertson timescales and by modelling the spectral energy distribution.
Results. The Kendall's tau correlation yields a probability of 76
that the mass of debris discs and their age are correlated. Similarly, the three two-sample tests give a probability between 70 and 83
that younger and older debris systems belong to different parent populations in terms of dust mass. We detected submillimetre/millimetre emission from six debris discs, enabling a detailed SED modelling.
Conclusions. Our results on the correlation and evolution of dust mass as a function of age are conditioned by the sensitivity limit of our survey. Deeper millimetre observations are needed to confirm the evolution of debris material around solar-like stars. In the case of the detected discs, the comparison between collision and Poynting-Robertson timescales supports the hypothesis that these discs are dominated by collisions. All detected debris disc systems show the inner part evacuated from small micron-sized grains.
Key words: circumstellar matter - planetary systems: formation - stars: late-type - Kuiper Belt
1 Introduction
As part of the star formation process in molecular cloud cores,
circumstellar discs form by conserving the initial angular momentum.
These so-called primordial discs have typical masses
(10
)
and sizes comparable to what is
expected for the primitive solar nebula, and they provide the environment
and material from which planets are expected to form (Beckwith et al. 2000; Natta et al. 2007). The first stages of planet formation occur close to the disc mid-plane where sub-micron grains grow to millimetre sizes and settle (e.g. Beckwith
et al. 2000). Larger planetesimals may successively
form by collisions between grains and through gravitational instabilities
(e.g. Johansen et al. 2007; Henning 2008). The primordial gas
in circumstellar discs is thought to last for 1-10 Myr (e.g. Haisch et al. 2001; Lawson et al. 2004).
After the gas is dispersed from the discs, the dust grains are removed
from the optically thin disc by radiation pressure, Poynting-Robertson drag
or dust sublimation on timescales shorter than the stellar pre-main
sequence lifetime (Backman & Paresce 1993; Meyer et al. 2007).
Nevertheless, an infrared excess above the stellar photosphere has been
discovered around more than 300 pre- and main-sequence stars.
These stars are thought to host a remnant debris disc.
Backman & Paresce (1993) suggested that the infrared excess of these objects
comes from small dust particles that are the products of collisions of larger bodies. Thus, a debris disc should consist of a ``second generation'' of small dust particles, larger bodies, and possibly, planets.
The formation of large bodies is thought to be faster in the inner part of the disc (Kenyon & Bromley 2004a). As a consequence, new small dust particles will be produced by collisions first in the inner region of the disc and, later, further out (Kenyon & Bromley 2004a). During this phase, small dust particles are continuously produced in the disc. In this perspective, transitional objects are of particular interest. These may represent a transitional phase between the primordial optically thick disc phase and the debris optically thin disc phase. Their infrared excess suggests the presence of an inner hole in the disc where the small grain population, which emits in the near-IR, is missing. Some of them are still accreting material onto the central star, while others show intense collisional activity in their inner part, possibly induced by the formation of planets (Currie et al. 2007; Kenyon & Bromley 2004b).
While the dissipation of debris discs takes place by radiation pressure
and/or Poynting-Robertson drag and the total dust mass decreases with time,
at early stages the stellar wind drag can be an important dust removal
mechanism. Depending on which of these processes is the main driver of dust
removal, we expect a power-law dependence with a different index (Dominik
& Decin 2003; Wyatt et al. 2007a,b; Löhne et al. 2008)
and a different dust distribution if a planetesimal belt at a given
distance is creating dust of just one size (Wyatt 2005).
Spangler et al. (2001) measured the fractional
luminosity (
)
of pre- and main sequence stars with
ISO. They found that the fractional luminosity, which is proportional
to the dust present in the disc, decreases with time as t-1.76.
While the analysis of Spangler et al. (2001) included also
pre-main sequence stars, Habing et al. (2001) analysed the
ISO and IRAS data of nearby main-sequence stars alone: they concluded that
most of the stars reach the main-sequence surrounded by a disc which decays
in about 400 Myr. Decin et al. (2003) reviewed the different
studies on the time dependency of Vega-like systems done with ISO/IRAS. They
concluded that the observations only indicated a large spread in fractional
luminosity for stars at most ages, that the maximum excess was independent of
time and that the debris disc phenomenon was more common in younger stars. An
analysis of the near to far infrared excess, carried out with Spitzer, led to
the finding of a general decrease of the fractional luminosity with time as
t-1, starting later at longer wavelengths (Su et al. 2006).
This suggested a disc clearing more efficiently in the inner parts
(Su et al. 2006). Deviations from these trends have been
associated with a recent or still ongoing collisional cascade that
produces a small dust grain population which is rapidly removed by radiation
pressure or by the action of the Poynting-Robertson drag
(Song et al. 2005; Wyatt et al. 2007b).
A complementary picture of the evolution of the outer part of the debris
discs is provided by (sub)millimetre observations. These show a decline in
the dust mass as a function of age and significant evolution of the
circumstellar dust mass within the first 10 Myr and between 10 Myr and a
few Gyr (e.g. Liu et al. 2004; Carpenter et al. 2005, hereafter C05).
In this paper we present the results of two deep sub-millimetre and
millimetre surveys of circumstellar discs around solar-type stars with ages
between 3 Myr and 3 Gyr. The characteristics of our sample are presented in Sect. 2.1. The
observations carried out at 350 m and at 1.2 mm are reported in
Sects. 2.2 and 2.3. The results from our two surveys are presented
in Sect. 2.4. The discussion of the dust disc mass evolution as a function of the age of
the systems, compared with the history of our solar system, is presented in
Sect. 3. A detailed analysis of the debris disc properties by modelling their
spectral energy distribution (SED) is provided in Sect. 4 and finally
the summary is presented in Sect. 5.
2 Observations
2.1 Sample
Nearby sources (d<150 pc) were selected from The Formation and Evolution of Planetary Systems (FEPS) Spitzer Legacy Program (Meyer et al. 2006). The FEPS sample contains 314 stars with stellar masses between 0.5 and
2
,
ages spanning the range from 3 Myr to 3 Gyr and spectral
types between K7 and F5. In addition, 14 objects were initially included
in the FEPS sample because of an infrared excess previously detected with
IRAS/ISO. Since these sources can introduce a bias in the FEPS
sample with respect to the presence of a disc, they are not taken
into account in any statistical analysis in this paper.
The characteristics of the entire star sample have been extensively analysed by
the FEPS team (e.g. Meyer et al. 2006; Mamajek & Hillenbrand 2008).
The stellar ages are based on pre-main sequence tracks for stars younger than 10 Myr, X-ray activity and the strength of the CaII H and K emission for nearby solar-type stars
(50 pc), and the association of stars with clusters or star-forming regions of known age. A new calibration presented in Mamajek & Hillenbrand (2008) has been applied. Final ages will be presented by Hillenbrand et al. (in preparation; see also Mamajek & Hillenbrand 2008). The distances of the sources were determined by the FEPS team on the basis of Hipparcos parallaxes for nearby stars and kinematic distances for stars associated with young moving groups and associations.
Our total sample contains 141 sources of which:
- 16 targets were observed at 350
m with the CSO (Caltech Submillimetre Observatory) 10-m telescope (see Col. 6 Table 1);
- 40 targets were observed at 1.2 mm with the IRAM (Institut de
Radioastronomie Millimétrique) 30-m telescope (see Col. 7 Table 1);
- 121 targets are taken from the SEST (Swedish-ESO Submillimetre Telescope) and OVRO (Owens Valley Radio Observatory) surveys performed by Carpenter at. al (2005), of which 17 have been re-observed during our surveys.

Table 1:
The stellar properties, e.g. spectral type, temperature,
luminosity and age range are reported in the first columns (from Meyer et al. 2006, 2008, Hillenbrand et al. 2008). CSO 350 m and IRAM 1.2 mm measured fluxes and upper limits are reported in the sixth and seventh column. The detected sources are highlighted in bold. In the ninth column are the references to the stellar distances (listed in the eighth column).
2.2 CSO observations at 350
m
Observations of the sources listed in Table 1 were carried out at the CSO between 17 and 21 April 2005, using the Submillimetre High Angular Resolution Camera II (SHARC-II) with the 350 m filter. The full width at half maximum (FWHM) for the 10.4 m telescope beam profile is 8.5
at this wavelength. The median zenith opacity at 225 GHz was 0.04. The typical
on-source integration times varied from 40 min to 2 h depending on
source elevation and zenith opacity, although a deeper 3.5 h integration
was obtained for HD 107146 (see Corder et al. 2009). Pointing and instrumental flux calibrators were selected from the CSO SHARC-II list of calibrators. Data imaging and
reduction were carried out using the Comprehensive Reduction Utility for
SHARC-II. The errors associated with each detection are the quadratic sum of the
instrumental error (between 5.3 and 45 mJy) and the calibration errors
(between 12
and 22
of the flux). Instrumental errors are computed
based on the rms computed in a sky annulus, and propagated over the
photometric aperture. Calibration errors were determined based on
repeatability of the flux calibrators over the entire observing run.
2.3 IRAM observations at 1.2 mm
Millimetre continuum observations of 40 sources at 1.2 mm were
carried out between May and December 2005 (Projects num. 013-05, 141-05)
and between January and April 2008 (Project num. 185-07) at the IRAM
30-m telescope at Pico Veleta with MAMBO2, the 117 channel IRAM
bolometer. The FWHM of the 30-m telescope beam is 11
at 1.2 mm.
The median zenith opacity at 230 GHz was 0.5 during the first run between May and
November 2005, and 0.3 during the second run carried out in December 2005 and
during the last run in 2008. Pointing and calibrator sources were selected
from a standard list provided by IRAM.
The observations were carried out in on-off integration mode using a
standard calibration. An rms of 0.7 mJy was obtained in one hour (blocks of three times 20 min) of on-off integration. The weather conditions were stable during
most of the nights and the calibrators were observed within a few hours from the science targets.
The data were reduced with the data reduction package MOPSIC (Zylka 2006).
Observations of the same target at different dates were coadded to increase the
signal to noise. The uncertainties associated with each flux are the quadratic sum of the
statistical error and the calibration error. The statistical error is the rms of the flux
over the entire telescope field of view, while the calibrator error is based on the agreement
between the predicted calibrator fluxes and those measured over 48 h
before and after the observations. We obtained statistical errors between
0.8 and 2.8 mJy and calibration errors between 11
and 16
in flux.
The final uncertainties are dominated by the statistical error in the IRAM observations.
2.4 Detections and upper limits
The fluxes (in the case of detections) and the 3
upper limits (in
the case of non-detections) measured at 350
m and 1.2 mm are listed in
Table 1. The detection errors and the upper limits include only
the statistical error. The calibration uncertainty is reported
in the note to Table 1. The 3
upper limits were computed as 3 times the
statistical error described above. 31 sources have been observed with IRAM,
7 sources only with CSO and 9 sources with both telescopes.
We detected 11 sources in total: 5 sources with IRAM, 8 with the CSO. Two
of the sources were detected with both IRAM and CSO.
Five of them are primordial discs younger than
10 Myr, and six are debris discs.
All sources were unresolved with the exception of HD 107146 which had a
deconvolved size at 350
m of
with an
uncertainty of about
per axis (Corder et al. 2009).
3 Disc dust mass
The main goal of this project is to analyse the temporal evolution of circumstellar debris discs and the mechanism of their dissipation around solar-type stars.
The dust mass in debris discs is expected to decrease with time. Depending on which process (collisions or Poynting-Robertson drag) is the main driver of dust removal, we expect a different power-law dependence of the dust mass as a function of the age of the system (Dominik & Decin 2003).
In this analysis we include the dust masses computed from the millimetre observations of nearby solar-type stars obtained with IRAM presented here and with SEST and OVRO (from C05).
The disc is assumed to be isothermal and optically thin ()
at millimetre wavelengths.
The dust disc mass can then be derived from
the millimetre flux, using the following equation:
where






In order to be consistent with C05, we compute the dust masses assuming ,
k0= 2 cm2 g-1 at 1.3 mm (Beckwith et al. 1990) and
K.
The uncertainty in the dust mass is computed
using error propagation. In the following, we discuss the principal systematic errors which are
present in the final dust mass uncertainty.
The mass absorption coefficient is poorly constrained. We consider as
limit the commonly used value for debris discs of k0=1.7 cm2 g-1 at
800 m adopted by Zuckerman & Becklin (1993),
which corresponds to k0=1.0 cm2 g-1 at 1.3 mm.
The mass absorption coefficient alone introduces at least a factor 2
in the final uncertainty associated to the mass.
In order to be consistent with the previous analysis of C05, we adopted a
temperature of 40 K. This represents a compromise between the
temperature associated with cold dust in protoplanetary discs (
K,
e.g. Beckwith et al. 1990) and warm dust
in debris discs (
K; Zuckerman & Song 2004a).
Since the mass is proportional to
,
the temperature
alone introduces a factor of 2 in the uncertainty of the masses.
Another source of uncertainty in the final mass error is introduced by the
source's distance. Nevertheless, the parallax of most of the sources older than 20 Myr
(which will be in the following analysed), has been measured by Hipparcos, and the distance uncertainty is directly computed from the
parallax uncertainty. This contributes an uncertainty of at least 10
of the dust masses.
![]() |
Figure 1: Disc mass versus age of the system. Filled circles are the masses determined from the IRAM detections (this work) and from the detections presented by C05. Upper limits are shown as filled triangles (new IRAM observations), empty and filled triangles (objects measured by C05 and re-observed here with higher sensitivity) and down arrows (from C05). The sources observed in both studies are connected by dashed lines. The vertical line separates the sources younger than 20 Myr and evolved systems older than 20 Myr. |
Open with DEXTER |
The 3
upper limits to the dust mass have been computed when the
sources were not detected with a signal to noise ratio of at least three.
The mass of the newly observed sources younger than 20 Myr is also plotted
as reference.
Taking into account the uncertainty on distance, temperature, measured
flux and mass absorption coefficient, the final uncertainty is
55
of the dust mass.
3.1 Evolution of the dust mass
Mass detections and upper limits are presented in Fig. 1 as a
function of age. In total, we plotted 132 upper limits and 9 detections
computed from the millimetre measurements. 17 sources have been observed
in both surveys, C05 and with IRAM (connected by a dotted line in Fig. 1).
Due to the higher sensitivity of the IRAM bolometer compared to SEST
and OVRO, we detected 4 discs with dust masses down to
,
about one order of magnitude smaller than
those in C05.
Primordial discs are shown in Fig. 1 only as comparison with systems older than 20 Myr but they are not taken into account in any statistical analysis in this paper.
The other variable which can strongly affect our analysis is the age determination of the discs older than 20 Myr. Their ages derived using different techniques are summarised in Table 2. The final ages are computed as the average and the standard deviation is adopted as uncertainty (Table 2). A more detailed discussion of the ages of the whole FEPS sample will be presented by Hillenbrand et al. (in preparation; see also Mamajek & Hillenbrand 2008).
In the following subsection we present the statistical analysis aimed at constraining the physical characteristics of the dust in circumstellar discs at ages older than 20 Myr, in order to understand the timescale and the dissipation mechanism in debris discs.
Table 2:
Ages (in logarithmic scale) of the debris discs detected at the
millimetre wavelength, with IRAM (this work), SEST and OVRO (from C05).
They have been derived using different techniques: the strength of the
CaII H and K (
)
and Lithium (Li) emission, X-ray activity
(X-ray), and the rotational period-age (
)
relation
(L. Hillenbrand, private communication). The ages are averaged and the
standard deviation is adopted as error.
3.1.1 Method: statistical analysis
A quantitative analysis of the dust mass-age relation was performed for
systems older than 20 Myr. We
carried out a survival analysis (Feigelson & Nelson
1985) in order to include the information contained
in the dust mass upper limits derived from the millimetre non-detections.
We use the program ASURV Rev. 1.2 (Lavalley et al. 1992), which
implements the survival analysis methods presented in Feigelson & Nelson
(1985) and Isobe et al. (1986).
The generalised Kendall's tau correlation (Brown et al. 1974) was used to measure the degree of correlation
between age and disc mass among debris discs. Such a method includes the
analysis of both detections and upper limits.
The probability that the disc mass is correlated with the age is 76.
To determine whether there is a different dust mass distribution in debris
disc systems older than 20 Myr, we used two-sample tests. We applied
three different tests: the Gehan, logrank, and Peto-Prentice tests (e.g.
Feigelson & Nelson 1985).
In order to investigate the evolution of the debris dust we divided the
sample of objects older than 20 Myr in two sub-samples with almost the same
number of objects: the first with ages between 20 and 180 Myr, and
the second sub-sample with ages between
180 Myr and 5 Gyr.
The three two-sample tests give a probability between 70 and 83
that younger and older debris systems belong to different parent populations
in terms of dust mass.
The relatively low significance of our results comes from the low number of
detected sources, compared to the upper limits, and the narrow range in
age of the detected debris discs. This leads to a high probability to assign
a detected disc to a different sub-sample in age.
We could not perform any conclusive linear regression to the data points
since the detected objects span a very narrow range in age and masses, and in general
there are too many upper limits compared to detections.
3.2 Discussion
The detection limit of our survey does not allow us to draw any firm conclusion about the evolution of the debris dust around solar-like stars as a function of age. Even more sensitive millimetre observations are needed to find statistically reliable trends. Our sample is biased toward sources with infrared excess. However, with the current sensitivity limit, additional observations of sources without infrared excess, which have a lower probability to be detected in the millimetre, would only increase the number of upper limits. This would not change the result of our statistical analysis.
We recall that all the mass derivations are valid only if the mass absorption coefficient k, which is inversely proportional to the mass of the dust in the disc, remains constant with time.
While in the protoplanetary disc phase we would expect a strong change of the mass absorption coefficient with time, in the debris disc phase we expect an equilibrium phase which is characterised by an almost constant second generation of dust. The assumption of a roughly time-independent dust opacity seems to be reasonable as long as we consider discs in the debris phase.
4 SED analysis of debris discs detected at 350
m and/or 1.2 mm
The debris discs detected by our surveys offer the opportunity to characterise and compare the disc properties around solar type stars. This is the case for the following objects: HD 104860, HD 8907, HD 377, HD 107146, HD 61005 and HD 191089.
We compiled the spectral energy distributions from the infrared data
available in the literature and our new detections at sub-millimetre and/or millimetre
wavelengths. Synthetic photometric points at 13 and 33 m, derived from
the IRS spectrum (Hillenbrand et al. 2008), were also
included in the SEDs. These data points allow us to distinguish the
infrared excess of the debris disc from the stellar photosphere,
determined by a Kurucz model. The IRS spectrum is overplotted in
Figs. 2-7, but was not considered during
the SED modelling.
To further characterise the debris discs detected during our survey, we compute:
- the
index from their SED to derive the dust grain size present in the disc;
- the collisional and Poynting-Robertson drag timescales;
- the blowout grain size (
);
- best fit model parameters by modelling the SED.
index - In the millimetre wavelength range, where
,
the Planck function can be approximated with the Rayleigh-Jeans relation (
).
The disc is assumed to be optically thin with the optical depth
.
Under these conditions
.
At wavelengths
mm the mass absorption coefficient
is
,
where
is the dust opacity index.
Since the slope of the SED in the millimetre wavelength range is proportional to
,
the measure of the spectral index
enables us to directly
constrain the
index via the relation
and to derive the
index
.
The index
is 2 for interstellar medium grain sizes, between 0 and
2 for pebbles of the order of 1 mm in size, and 0 for even larger grains
(e.g. Beckwith et al. 2000).
The slope
of the SED in the sub-millimetre/millimetre wavelength
range is derived, using a linear regression that gives the errors in slope
of the best-fit. The values for the
index of the power law of the
dust absorption coefficient are always <1 (see Table 3),
which suggests the presence of dust with at least millimetre size particles
(e.g. Henning & Stognienko 1996; Rodmann et al. 2006).
Blowout grain size - Small grains are affected by the
interaction with the stellar radiation field. This causes a force
acting on the particles which is parametrised by the ratio of the
radiation force to the stellar gravity. We
compute the blowout grain size (
)
in the debris disc systems
using the equations presented in Hillenbrand et
al. (2008; based on Burns et al. 1979;
and Backman & Paresce 1993):
![]() |
(2) |
where,

Poynting-Robertson and collisional timescales -
Defining
(Wyatt et al. 1999-2005), where
and
are the Poynting-Robertson drag (
)
and collisional timescales respectively, it is possible to distinguish between discs in collisional
regime (
)
and Poynting-Robertson drag regime (
).
When
the disc is dense and the collisions occur faster than the P-R drag:
the dynamically bound dust remains at the same radial location as
the planetesimals, while unbound grains are blown out and their
surface density distribution falls off as r-1 (where r is
the distance from the star).
The Poynting-Robertson drag (
)
and collisional timescales
(
)
of the debris discs are computed, using the equations
presented in Hillenbrand et al. (2008; based on
Burns et al. 1979; and Backman & Paresce 1993):
![]() |
(3) |
where,


Under the assumption of circular orbits and completely destructive
collisions between grains of the same size, the collisional timescale
is computed using the equation:
![]() |
(4) |
where





The collisional timescale has been computed for grains of 10 m in size at the inner radius of the disc (computed by our model, see Table 4). According to Wyatt (2005), the ratio
is proportional to
,
where a is the grain size.
These values are summarised in Table 3.
Table 3:
Blowout grain size, P-R drag timescales of the debris discs detected in the CSO and/or IRAM surveys. The -index discs only for the stars
with the SED sampled in the submillimetre/millimetre range
(details in Sect. 4) is reported.
Table 4:
Model results, reduced
(
)
assuming a fixed outer radius (
)
of 150 AU and a maximum grain size (
)
of 3 mm.
SED modelling - The SEDs were modelled using the radiative transfer model of Wolf & Hillenbrand (2003).
The model assumes that the dust grains are compact spherical particles heated only by the direct stellar radiation. The interaction between stellar radiation and dust particles can be described by the following radiative processes: scattering (neglecting multiple scattering), absorption and re-emission of stellar radiation by dust grains. The disc is then checked to be optically thin to the stellar radiation and to the dust re-emission at all wavelengths. The contribution to the SED from a single dust grain results from the integration of the scattering and re-emission processes over all wavelengths. Finally the emergent spectrum is the sum of all the dust grain contributions.
The model parameters are the radiation emitted by the central star, the
total disc mass, the disc size, the radial density distribution, the
minimum and maximum grain size in the disc, the grain size distribution,
and the chemical composition of the grains.
The radial density distribution is described by a power-law,
where q=1 corresponds to a disc with a constant surface density
.
The grain size distribution is also described by a power-law
where the canonical value p=3.5 characterises a size distribution
initially produced by a collisional cascade. We always use this
grain size distribution in our models.
For each grain size a temperature distribution over the disc
size is computed.
The model parameter space has been extensively analysed by
Wolf & Hillenbrand (2003). They found that the
increase of the inner radius
causes a loss of the warm dust that
is mainly responsible for the near-infrared (NIR)/ mid-infrared (MIR) shape of
the SED; the excess emission over the stellar photosphere starts at longer
wavelengths. This is more pronounced when the minimum grain size
increases, because the main contribution to the NIR/MIR spectrum
at high temperatures comes from small grains. The presence of an inner gap
causes the loss of warm dust which is mainly responsible for the NIR/MIR
shape of the SED. Increasing the size of the gap, the MIR flux decreases
and the excess is shifted to longer wavelengths. Keeping the disc mass
constant, the flux in the millimetre region increases slightly, but
the net flux is smaller compared to a disc without an inner gap.
This is because the fraction of the stellar flux absorbed by a single
dust grain decreases with increasing radial distance from the
star.
We assumed astronomical silicate grains (optical data from Draine & Lee 1984; and Weingartner & Draine 2001) since they are expected to be the main component of the dust in protoplanetary discs (Pollack et al. 1994).
Since the outer radius cannot be constrained from the SED alone, it was assumed
to be 150 AU, consistent with the location of the most distant object of the
Kuiper belt in our Solar System (e.g. Gladman et al. 2002).
The choice of such a value is supported by a recent resolved millimetre map of one of the debris
discs analysed (Corder et al. 2009), where the
millimetre emission is extended up to 150 AU from the central star.
The maximum grains size is assumed to be 3 mm, since the
index measurements suggested the presence of at least
millimetre-sized particles in the discs; this is the grain size
which strongly contributes to the millimetre emission that we detect.
The Levenberg-Marquardt method was applied to find the best-fit parameters
of the model leaving the inner radius (
), minimum grain
size (
) and dust mass (
) as free parameters.
The reduced ,
associated with each best-fit model, was computed
for each fit normalising the standard statistic
to the number of
free parameters of the fit (Bevington 1969). These parameter uncertainties
are estimated from the covariance matrix (inverse of the
curvature function) by the Levenberg-Marquardt algorithm. The modelling results are listed in Table 4. The degeneracy of the SED models was checked using the
-confidence regions delineated by lines of constant
.
The model parameter space (
)
was explored around the best fit values. The 2-dimensional confidence regions in Figs. 2-7
are the projections of the 3-dimensional isosurfaces defined by
fixed
(e.g. Press et al. 1992).
The shape of the confidence regions confirm an intrinsic degeneracy
of our model probably between inner radius, minimum grain size and
dust mass.
![]() |
Figure 2:
SED of HD 104860 (upper plot a)) detections are
shown as empty diamonds while 3 |
Open with DEXTER |
![]() |
Figure 3: As Fig. 2 but for HD 8907. |
Open with DEXTER |
![]() |
Figure 4: As Fig. 2 but for HD 377. |
Open with DEXTER |
![]() |
Figure 5: As Fig. 2 but for HD 107146. |
Open with DEXTER |
![]() |
Figure 6: As Fig. 2 but for HD 61005. |
Open with DEXTER |
![]() |
Figure 7: As Fig. 2 but for HD191089. |
Open with DEXTER |
4.1 HD 104860
HD 104860 is an F8 star at a distance of 47 pc from the Sun.
It was classified as a Doradus variable candidate by
Handler (1999). Wichmann et al. (2003)
identified the object as a young zero age main sequence star using LiI
in absorption (Wichmann et al. 2003).
In Table 2 we summarise the ages derived from different
methods (L. Hillenbrand, private communication). The average age is
140 Myr,
the value that we adopted in our analysis.
The SED of HD 104860 is presented in Fig. 2. It includes the Spitzer/IRAC and Spitzer/MIPS photometry from 3.6 to 160 m, synthetic photometry
at 13 and 33
m obtained from the IRS low resolution spectrum
(Hillenbrand et al. (2008), our CSO and IRAM detections at
350
m and 1.2 mm, SCUBA detections at 450 and 850
m
(Najita & Williams 2005), and an OVRO
3.1 mm upper limit (Carpenter et al. 2005). The infrared-excess
emission detected by Spitzer suggests the presence of debris
material around this main sequence star (Carpenter et al. 2008).
Early modelling of the spectral energy distribution of HD 104860 was
presented by Najita & Williams (2005). Treating dust
grains as simple blackbody radiators, their sub-millimetre observations
at 450 and 850 m indicate a dust temperature of
K and a mass
of the disc of 0.16
.
Our new detection allows us to sample the SED at longer wavelengths and hence constrain the slope of the far-infrared to millimetre excess.
The best-fit model parameters are reported in Table 4 and
the computed best-fit model is plotted in Fig. 2. The
confidence regions of the fit (see lower plots in Fig. 2)
show some degree of degeneracy of the model.
The dust mass is constrained between 0.07
and 0.15
,
while the inner radius is located between 10 and 50 AU and the minimum
grain size is between 5 and 15
m.
The dust mass computed by the model (see Table 4) is of the
same order of magnitude as the masses computed using a simple black-body
model (Najita & Williams 2005) and Eq. (1).
The minimum grain size obtained is one order of magnitude higher than the
blowout grain size of the system.
The confidence regions delineated by constant
contours
suggest that a minimum grain size of the same order of magnitude of the
blowout size is a poor fit to the data.
4.2 HD 8907
HD 8907 is a F8 star, located at 34 pc from the Sun, which has already reached the main sequence. In our analysis we adopt an age of the system of 680 Myr, which is the average of the ages derived from different methods (Table 2, Hillenbrand, private communication).
The SED of HD 8907 is presented in Fig. 3. The ISO photometry at 60 and 90 m
are from Silverstone (2000), the SCUBA detections at 450 and 850
m from Najita & Williams (2005), IRAM detection at
1.2 mm from this work and the OVRO 2.6 and 3.1 mm upper limits from C05.
The IRAS detections (Helou & Walker 1988; empty
squares in Fig. 3), are systematically larger than the Spitzer photometry,
due to the bigger beam and they were not taken into account during the model fitting.
The infrared excess detected by ISO (Zuckerman & Song 2004a)
suggested the presence of debris dust. Under the assumption that the dust
particles radiate as a black body, a temperature of
60 K was inferred
for particles orbiting at
30 AU from the central star
(Zuckerman & Song 2004a). The SED infrared excess has
also been detected by Spitzer between 3.6 and 70
m by
Kim et al. (2005),
and they found an inner radius of
48 AU with a
temperature of
48 K.
Our new detection allows us to sample the SED at longer wavelengths
and hence constrain the slope of the far-infrared to millimetre excess.
The best-fit model parameters are reported in Table 4 and the
computed best-fit model is plotted in Fig. 3. The confidence
regions of the fit which show some degree of degeneracy of the model
(see lower plots in Fig. 3): the dust mass is constrained
between 0.03
and 0.1
,
while the inner radius is
located between 15 and 30 AU and the minimum grain size is between 5
and 10
m.
As in the case of HD 104860, the dust mass computed by the model of
HD 8907 (see Table 4) is comparable to the value computed
using Eq. (1). The minimum grain size is between 4 and 8 m.
This is one order of magnitude higher than the blowout grain size of
the system and is consistent with previous modelling
(Kim et al. 2005). The inner radius computed by our model
is poorly constrained since it has an uncertainty of 50%.
4.3 HD 377
HD 377 is a G2V star at a distance of 40 pc which has already reached the
main sequence. It is surrounded by a debris disc with a high
fractional luminosity (
),
(Moór et al. 2006).
The ages of HD 377 derived from different methods (L. Hillenbrand, private
communication) are summarised in Table 2. The average age is
110 Myr,
the value that we adopted in our analysis.
Pascucci et al. (2006) investigated the presence of gas in the disc of HD 377, but they did not detect any gas emission lines in the infrared nor in the millimetre wavelength range.
The SED was compiled from the literature including our new detection at 1.2 mm (Fig. 4). The OVRO 2.6 and 3.1 mm upper limits are from C05. The best-fit model parameters are reported in Table 4 and the computed best-fit model is plotted in Fig. 4. The confidence regions of the fit (see lower plots in Fig. 4) show a high degree of degeneracy of the model compared to the other debris discs modelled. This results from the poor constraints on the SED in the sub-millimetre region.
As for HD 104860 and HD 8907, the dust mass computed by the model (see
Table 4) is comparable to the value computed using Eq. (1).
The inner radius is the smallest in our sample, but the disc seems to lack
small dust grains with sizes up to 14 m. The minimum grain size
obtained is one order of magnitude higher than the blowout grain size
of the system.
4.4 HD 107146
HD 107146 is a G2 V star located 28.5 pc from the Sun. The star is a candidate periodic V-band photometric variable (Koen & Eyer 2002).
From the average between the ages derived with different methods
(Table 2), we obtain an age of 230 Myr.
HD 107146 is also the first debris disc around a solar-type star resolved
in scattered light by the Hubble Space Telescope
(Ardila et al. 2004). They resolved a ring-like disc with most of the material concentrated between 86 and 130 AU, using a coronograph of 1.8
in radius. The disc
was also marginally resolved at 450 and 850
m and
the observations suggested the presence of an inner hole with a radius larger than 31
AU (Williams et al. 2004). Such an inner hole would also explain the lack of an IRAS
25
m excess and a far-infrared excess at 60 and 100
m
(Silverstone 2000). C05 presented OVRO 3 mm images of this source.
HD 107146 has also been resolved during our survey at 350
m with CSO
and at 1.3 mm with CARMA (Corder et al. 2009). The
millimetre image shows a clumpy disc extended between
60 AU and
150 AU
with a symmetric peak of the emission at
97 AU from the central
star.
The SED was compiled from the literature, including our new detection at
350 m (Fig. 5) and the CARMA detection at 1.3 mm
(Corder et al. 2009).
HD 107146 is the only disc detected at 3.1 mm. The IRAS detections
(Helou & Walker 1988; empty squares in Fig. 5), are systematically
higher than the Spitzer photometry, due to the bigger beam and they
were not taken into account during the model fitting.
The best-fit model parameters are reported in Table 4
and the computed best-fit model is plotted in Fig. 5. The
degree of degeneracy of the model is shown in the confidence regions of
the fit (see lower plots in Fig. 5).
The dust mass is constrained between 0.08
and
0.15
,
while the inner radius is located between
7 and 15 AU and the minimum grain size is between 8 and 15
m.
It is interesting to notice that the inner radius computed by our model,
is more than 7 times smaller than the inner radius of the ring-like disc in
the HST images. This is because different grain populations are traced by
our model and the HST image: while the HST traces only the distribution of
the sub-micron size population (below
), our model includes
the re-emission from grains from micron to millimetre in size.
This explanation is also supported by the resolved images in the sub-millimetre
wavelength range of HD 107146 which suggest the presence of grains
at radii smaller compared to the ring-like structure seen in the HST image.
4.5 HD 61005
HD 61005 is a G3/G5V star with a Hipparcos distance of pc,
located in the local bubble, a region that is thought to be almost
completely free of diffuse dust clouds (Franco et al. 1990).
From an average of the chromospheric activity available from the
literature, and a new calibration of the chromospheric activity-age
relation (FEPS, private communication), we adopt an age of 135 Myr.
The SED was previously modelled by Hines et al. (2007):
under the assumption that dust particles of
10
m radiate as
a black body with a temperature of
50-70 K, they obtained a minimum
distance for the circumstellar material of
7 AU. This star was also
monitored during a radial velocity survey looking for planets around active
stars, but no planets were detected (Setiawan et al. 2008).
HST/NICMOS observations (Hines et al. 2007) reveal
dust-scattered starlight extending to distances of
240 AU from
the occulted star. The structure is strongly asymmetric about its major
axis, but is mirror-symmetric about its minor axis; morphologically, the
object resembles a wing-spread moth with the star as the head.
HST scattered light images reveal a swept shape of the disc with
an inner radius 10 AU and an outer radius up to 240 AU (Hines
et al. 2007). Such emission is thought to be associated
to the local interstellar medium that scatters the stellar light, so that the
movement of HD 61005 through the medium causes the swept shape of the
disc (Hines et al. 2007). These observations give
a complementary picture of the circumstellar material present in
the disc since different grain populations are traced in our model
and in the scattered light images.
The SED was compiled from the literature including our new detection at
350 m (Fig. 6) and the SEST upper limit at 1.2 mm from C05.
The best-fit model parameters are reported in
Table 4 and the computed best-fit model is plotted in
Fig. 6. The confidence regions of the fit (see lower plots in
Fig. 6) show that the best fit parameters represent only a local
minimum of the model parameters space. This is probably due to the fact that
the model assumptions are not valid in this case. As suggested from the HST
images, the object seems to be located in the local interstellar medium
material, which probably introduces an additional small grain population
which contributes to the scattered light in the SED. Our model assumes
that all the detected flux is stellar light scattered and/or absorbed and re-emitted
only by the dust particles in the disc.
4.6 HD 191089
HD 191089 is an F5V star with a Hipparcos distance of pc from the
Sun. This star was selected for the FEPS programme because it was known
to have a debris disc based on ISO and IRAS measurements and it was
not considered in our statistical analysis presented in Sect. 3.
An upper limit of the age (
3 Gyr) for the system was presented by
Nordström et al. (2004), using main-sequence
isochrones. We adopt an age of
300 Myr (Hillenbrand et al. 2008).
The SED was compiled from the literature (Fig. 7) including the Spitzer spectro-photometry (Carpenter et al. 2008), CSO detection at 350 m and IRAM upper limit
at 1.2 mm from this work and OVRO 3.1 mm upper limit from C05. The IRAS detections
(Helou & Walker 1988; empty squares in Fig. 7),
are systematically higher than the Spitzer photometry,
due to the bigger beam and they were not taken into
account during the model fitting. Since there is no
IRS spectrum of the source constraining the beginning
of the infrared excess close to the stellar photosphere,
the IRAS detection at 12
m was taken into account during the SED fitting.
The best-fit model parameters are reported in Table 4 and the
computed best-fit model is plotted in Fig. 7. The confidence
regions of the fit show some degree of degeneracy of the model (see lower
plots in Fig. 7).
We note that the model result is not significant at all. This is due to
the few data points which constrain the SED.
4.7 Modelling discussion
Comparing the collisional timescale to the Poynting-Robertson timescale, we conclude that all the detected debris disc are systems in the collision dominated regime.
The SED of the debris discs has been modelled leaving 3 parameters simultaneously free: inner radius, dust mass and minimum grain size.
We notice that the slope of the SEDs seems to change from the sub-millimetre to the millimetre wavelength range, but is well reproduced by our models within the errors. Including grains larger than 3 mm in the model would decrease the steepness of the slope in the sub-millimetre, but would not reproduce at the same time the sub-millimetre and millimetre fluxes.
A degeneracy between the inner radius and minimum grain size
can still exist:
in this case the -contours plots would represent only a local
minimum in the parameter space, instead of the absolute solution.
Only a resolved image of the entire grain population which delineates
the inner radius of the disc will allow us to break the degeneracy between
smallest grain size present in the disc and inner radius.
We now discuss how well constrained the fixed model parameters are
(i.e.
,
,
dust composition and size distribution) and in
particular if and how they can alter our conclusions about
,
and
.
We are not going to further
discuss the cases of HD 61005 and HD 191089 because our modelling was not
significant at all, nor the case of HD 107146 because the assumption on the outer
radius is supported by millimetre resolved images (Corder et al. 2009).
- The outer radius was assumed to be 150 AU. While in the case of
HD107146 such an assumption is supported by resolved millimetre images,
in the other cases, we need to investigate the robustness of our model
results. While the most distant objects of the Kuiper belt have been
observed up to 150 AU from the Sun (e.g. Gladman et al. 2002),
the edge of the Kuiper belt is commonly considered close to the 2:1 resonance
with Neptune at about 50 AU from the Sun (e.g. Trujillo & Brown 2001).
There is also no reason why all discs would be expected to have the same outer edge. For these
reasons we use five different outer radii (1000, 500, 300, 100
and 50 AU, or 60 AU when the fitting routine was not
able to find any reasonable inner radius smaller than 50 AU), leaving
free the dust mass (
), the inner radius (
) and the minimum grain size
. The results are summarised in Table 5.
Table 5: Model results leaving as free parameters using
We notice that, in the case of HD 104860 and HD 8907, changing the outer radius from 500 to 50 AU affects only the dust mass. In the case of HD 377 fixing the outer radius to 100 and 50 AU the model fitting is no longer satisfactory, while we find a better result increasing the outer radius up to 500 AU (see different,
and
for 6 different outer radii: 1000, 500, 300, 150, 100 and 50 AU.
in Table 5). The minimum grain size remains at least one order of magnitude bigger than the blowout size, despite the case of an outer radius of 1000 AU where it becomes of the same order of magnitude or smaller than the blowout size.
- The maximum grain size has been fixed to 3 mm; despite the fact that a
collisional cascade generates a distribution of grain size which does not
stop at the millimetre size, we include in our model the maximum particle
size which contributes significantly to the emission which we detect
during our observations. We remind the reader that the mass derived
by our model represents the mass of the grains smaller than 3 mm
present in the disc and that increasing the maximum grain size would increase
the disc mass.
- The dust has been assumed to be made up of astronomical silicate;
this is expected to be the main component in protoplanetary discs (Pollack
et al. 1994). Although a different composition of the dust
can be present in the disc, Spitzer infrared spectroscopic observations do not
allow us to constrain the material properties due to the lack of diagnostic features,
caused by large grain sizes (see the IRS spectra overplotted in the SEDs, Figs. 2-7).
A different dust composition would change the blow-out size as
well as the minimum grain size present in the disk computed by our model. A detailed
investigation of different dust compositions is beyond the goal of this paper.
- We finally recall that a more realistic grain size distribution can also affect our modelling results. Departures from the canonical size distribution we assumed, have been highlighted by collisional evolution studies of spatially resolved debris disks (e.g. Thébault & Augereau 2007).
5 Summary and conclusions
We have carried out two deep surveys at 350 m and 1.2 mm of
circumstellar discs around solar-type stars at ages between 3 Myr
and 3 Gyr. These new observations are an order of magnitude
more sensitive to dust emission than previous observations of C05.
The dust disc masses have been computed from the millimetre emission,
where the discs are assumed to be optically thin. A survival analysis
of the dust disc masses as a function of time has been carried out
of systems older than 20 Myr, including the dust mass upper limits.
The spectral energy distributions of the debris discs detected by
our sub-millimetre and millimetre surveys have been modelled. We
draw the following main conclusions from our work:
- 1.
- The Kendall's tau correlation yields a probability of 76
that the mass of debris discs and their age are correlated. Similarly, the three different two-sample tests gives a probability between 70 and 83
that younger and older debris systems belong to different parent populations in terms of dust mass. Our result on the relation between dust mass and age is limited by the sensitivity of our millimetre survey. Deeper millimetre observations are needed to confirm the evolution of debris material around solar-like stars.
- 2.
- The spectral energy distributions of the debris discs detected at
350
m and/or 1.2 mm were modelled. We found a degeneracy in the best fit parameters which will only be broken with high-resolution images which resolve the entire disc. Nevertheless, this approach allows us to identify debris discs with an inner region that has been evacuated from small micron-sized dust grains.
- 3.
- In the case of the detected debris discs, the comparison between collision and Poynting-Robertson timescales suggests that the debris discs are collision dominated.
Acknowledgements
All the authors wish to thank the entire FEPS team for their work. V.R. thanks A. Sicilia-Aguilar for helpful discussions on the IRAM data reduction, E. Feigelson for his suggestions about the survival analysis, L. Hillenbrand for some details about the stellar ages, S. J. Kim for her suggestions about the FEPS sample, A. Pasquali and A. Martinez-Sansigre for reading the manuscript. S.W. was supported at the MPIA by the German Research Foundation (DFG) through the Emmy Noether grant WO 857/ 2. J.R. wishes to thank Jesús Falcón-Barroso for a helpful discussion on how to project 3-d isosurfaces to lower dimensions. M.R.M. thanks support provided through the LAPLACE node of the NASA Astrobiology Institute. Research at the Caltech Submillimetre Observatory is supported by grant AST-0540882 from the National Science Foundation.
References
- Andrews, S. M., & Williams, J. P. 2005, ApJ, 631, 1134
- Andrews, S. M., & Williams, J. P. 2007, ApJ, 659, 705
- Ardila, D. R., Golimowski, D. A., Krist, J. E., et al. 2004, ApJ, 617, L147 In the text
- Backman, D. 2004, in Debris Disks and the Formation of Planets, ed. L. Caroff, et al. (San Francisco: ASP), ASP Conf. Ser., 324, 9 In the text
- Backman, D. E., & Paresce, F. 1993, Protostars and Planets III, 1253 In the text
- Beckwith, S. V. W., Sargent, A. I., Chini, R. S., & Güsten, R. 1990, AJ, 99, 924 In the text
- Beckwith, S. V. W., Henning, T., & Nakagawa, Y. 2000, Protostars and Planets IV, ed. V. Mannings, A. P. Boss, & S. S. Russell, 533 In the text
- Bevington, P. R. 1969, Data Reduction and Error Analysis for the Physical Sciences (New York: McGraw-Hill), 1 In the text
- Bottke, W. F., Durda, D. D., Nesvorný, D., et al. 2005, Icarus, 179, 63
- Brown, B. W. M., Hollander, M., & Korwar, R. M. 1974, Reliability and Biometry, ed. F. Proschan, & R. J. Serfling (Philadelphia: SIAM), 327 In the text
- Buckley, J., & James, I. 1979, Biometrika, 66, 429
- Burns, J. A., Lamy, P. L., & Soter, S. 1979, Icarus, 40, 1 In the text
- Calvet, N., D'Alessio, P., Watson, D. M., et al. 2005, ApJ, 630, L185
- Carpenter, J. M., Wolf, S., Schreyer, K., Launhardt, R., & Henning, T. 2005, AJ, 129, 1049 In the text
- Carpenter, J. M., Mamajek, E. E., Hillenbrand, L. A., & Meyer, M. R. 2006, ApJ, 651, L49
- Carpenter, J. M., Bouwman, J., Mamajek, E. E., et al. 2008, [arXiv:0810.1003] In the text
- Chen, C. H., Jura, M., Gordon, K. D., & Blaylock, M. 2005, ApJ, 623, 493
- Corder, S. A., Carpenter, J. M., Sargent, A. I., et al. 2009, ApJ, 690, L65 In the text
- Currie, T., Kenyon, S. J., Rieke, G., Balog, Z., & Bromley, B. C. 2007, ApJ, 663, L105 In the text
- Cutri, R. M., et al. 2003, The IRSA 2MASS All-Sky Point Source Catalog, NASA/IPAC Infrared Science Archive, http://irsa.ipac.caltech.edu/applications/Gator/
- de Zeeuw, P. T., Hoogerwerf, R., de Bruijne, J. H. J., Brown, A. G. A., & Blaauw, A. 1999, AJ, 117, 354
- Dominik, C., & Decin, G. 2003, ApJ, 598, 626 In the text
- Decin, G., Dominik, C., Waters, L. B. F. M., & Waelkens, C. 2003, ApJ, 598, 636 In the text
- Draine, B. T., & Lee, H. M. 1984, ApJ, 285, 89 In the text
- Feigelson, E. D., & Nelson, P. I. 1985, ApJ, 293, 192 In the text
- Franco, J., Tenorio-Tagle, G., & Bodenheimer, P. 1990, ApJ, 349, 126 In the text
- Ghez, A. M., Neugebauer, G., & Matthews, K. 1993, AJ, 106, 2005
- Gladman, B., Holman, M., Grav, T., et al. 2002, Icarus, 157, 269 In the text
- Habing, H. J., Dominik, C., Jourdain de Muizon, M., et al. 2001, A&A, 365, 545 In the text
- Haisch, K. E., Jr., Lada, E. A., & Lada, C. J. 2001, ApJ, 553, L153 In the text
- Handler, G. 1999, MNRAS, 309, L19 In the text
- Hartmann, W. K., Ryder, G., Dones, L., & Grinspoon, D. 2000, Origin of the earth and moon, ed. R. M. Canup, K. Righter and 69 collaborating authors (Tucson: University of Arizona Press), 493
- Helou, G., & Walker, D. W. 1988, Infrared astronomical satellite IRAS catalogues and atlases, 7, 1 In the text
- Henning Th. 2008, Phys. Scr. T, 130, 014019 In the text
- Henning, Th., & Stognienko, R. 1996, A&A, 311, 291 In the text
- Hillenbrand. L. A. 2005, [arXiv:astroph/0511083]
- Hillenbrand, L. A., Carpenter, J. M., Kim, J. S., et al. 2008, ApJ, 677, 630 In the text
- Hines, D. C., Schneider G., Hollenbach D., et al. S. 2007, ApJ, 671, 165 In the text
- Isobe, T., Feigelson, E. D., & Nelson, P. I. 1986, ApJ, 306, 490 In the text
- Johansen, A., Oishi, J. S., Low, M.-M. M., et al. 2007, Nature, 448, 1022 In the text
- Kenyon, S. J., & Bromley, B. C. 2001, AJ, 121, 538
- Kenyon, S. J., & Bromley, B. C. 2002, AJ, 123, 1757
- Kenyon, S. J., & Bromley, B. C. 2002, ApJ, 577, L35
- Kenyon, S. J., & Bromley, B. C. 2004, AJ, 127, 513 In the text
- Kenyon, S. J., & Bromley, B. C. 2004, ApJ, 602, L133 In the text
- Kenyon, S. J., & Hartmann, L. 1995, ApJS, 101, 117
- Kim, J. S., Hines, D.C., Backman, D. E., et al. 2005, ApJ, 632, 659 In the text
- Koen, C., & Eyer, L. 2002, VizieR Online Data Catalog, 733, 10045 In the text
- Lavalley, M., Isobe, T., & Feigelson, E. 1992, Astronomical Data Analysis Software and Systems I, ed. D. M. Worrall, C. Biemesderfer, & J. Barnes (San Francisco: ASP), 245 In the text
- Lawson, W. A., Lyo, A.-R., & Muzerolle, J. 2004, MNRAS, 351, L39 In the text
- Levison, H. F., Morbidelli, A., Gomes, R., & Backman, D. 2007, Protostars and Planets V, 669
- Liu, M. C., Matthews, B. C., Williams, J. P., & Kalas, P. G. 2004, ApJ, 608, 526 In the text
- Löhne, T., Krivov, A. V., & Rodmann, J. 2008, ApJ, 673, 1123 In the text
- Mamajek, E. E., Meyer, M. R., Hinz, P. M., et al. 2004, ApJ, 612, 496
- Mamajek, E. E., & Hillenbrand, L. A. 2008, ArXiv e-prints, 807, [arXiv:0807.1686]
- Meyer, M. R., Hillenbrand, L. A., Backman, D., et al. 2006, PASP, 118, 1690 In the text
- Meyer, M. R., Backman, D. E., Weinberger, A. J., & Wyatt, M. C. 2007, Protostars and Planets V, ed. B. Reipurth, D. Jewitt, & K. Keil (Tucson: University of Arizona Press), 951, 573 In the text
- Meyer, M. R., Carpenter, J. M., Mamajek, E. E., et al. 2008, ApJ, 673, L181 In the text
- Moór, A., Ábrahám, P., Derekas, A., et al. 2006, ApJ, 644, 525 In the text
- Najita, J., & Williams, J. P. 2005, ApJ, 635, 625
- Natta, A., Testi, L., Neri, R., Shepherd, D. S., & Wilner, D. J. 2004, A&A, 416, 179
- Natta, A., Testi, L., Calvet, N., et al. 2007, Protostars and Planets V, ed. B. Reipurth, D. Jewitt, & K. Keil (Tucson: University of Arizona Press), 767 In the text
- Nordström, B., Mayor, M., Andersen, J., et al. 2004, A&A, 418, 989 In the text
- Pan, X., Shao, M., & Kulkarni, S. R. 2004, Nature, 427, 326
- Pascucci, I., Gorti, U., Hollenbach, D., et al. 2006, ApJ, 651, 1177 In the text
- Perryman, M. A. C., Lindegren, L., Kovalevsky, J., et al. 1997, A&A, 323, L49
- Pinsonneault, M. H., Stauffer, J., Soderblom, D. R., King, J. R., & Hanson, R. B. 1998, ApJ, 504, 170
- Pollack, J. B., Hollenbach, D., Beckwith, S., et al. 1994, ApJ, 421, 615 In the text
- Prato, L. 2007, ApJ, 657, 338
- Preibisch, T., Guenther, E., Zinnecker, H., et al. 1998, A&A, 333, 619
- Preibisch, T., Brown, A. G. A., Bridges, T., Guenther, E., & Zinnecker, H. 2002, AJ, 124, 404
- Press, W. H., Flannery, B. P., Teukolsky, S. A., & Vetterling, W. T. 1992, Numerical Recipes in C, 2nd edition (Cambridge University Press), Sect. 15.6 In the text
- Rodmann, J., Henning, T., Chandler, C. J., Mundy, L. G., & Wilner, D. J. 2006, A&A, 446, 211 In the text
- Roques, F., Scholl, H., Sicardy, B., & Smith, B. A. 1994, Icarus, 108, 37 In the text
- Silverstone, M. D., Hines, D. C., Schneider, G., et al. 2006, ApJ, 650, 414
- Smith, B. A., Becklin, E. E., Schneider, G., et al. 1999, ApJ, 513, L127
- Setiawan, J., et al. 2008, Precision Spectroscopy in Astrophysics, 201 In the text
- Silverstone, M. D. 2000, Ph.D. Thesis In the text
- Silverstone, M. D., Meyer, M. R., Mamajek, E. E., et al. 2006, ApJ, 639, 1138
- Spangler, C., Sargent, A. I., Silverstone, M. D., Becklin, E. E., & Zuckerman, B. 2001, ApJ, 555, 932 In the text
- Benedict, G. F., Nelan, E., Soderblom, D. R., et al. 2004, 36, 735
- Song, I., Zuckerman, B., Weinberger, A. J., & Becklin, E. E. 2005, Nature, 436, 363 In the text
- Strom, S. E., Edwards, S., & Skrutskie, M. F. 1993, Protostars and Planets III, 837
- Rieke, G. H., Stransberry, J. A., Su, K. Y. L., et al. 2006, ApJ, 653, 675 In the text
- Thébault, P., & Augereau, J.-C. 2007, A&A, 472, 169
- Trujillo, C. A., & Brown, M. E. 2001, ApJ, 554, L95
- Weingartner, J. C., & Draine, B. T. 2001, ApJ, 548, 296 In the text
- Wichmann, R., Schmitt, J. H. M. M., & Hubrig, S. 2003, A&A, 399, 983 In the text
- Williams, J. P., & Andrews, S. M. 2006, ApJ, 653, 1480
- Williams, J. P., Najita, J., Liu, M. C., et al. 2004, ApJ, 604, 414 In the text
- Wolf, S., & Hillenbrand, L. A. 2003, ApJ, 596, 603 In the text
- Wright, J. T., Marcy, G. W., Butler, R. P., & Vogt, S. S. 2004, ApJS, 152, 261
- Wyatt, M. C. 2005, A&A, 433, 1007 In the text
- Wyatt, M. C., Dermott, S. F., Telesco, C. M., et al. 1999, ApJ, 527, 918 In the text
- Wyatt, M. C., Dent, W. R. F., & Greaves, J. S. 2003, MNRAS, 342, 876
- Wyatt, M. C., Smith, R., Greaves, J. S., et al. 2007a, ApJ, 658, 569 In the text
- Wyatt, M. C., Smith, R., Su, K. Y. L., et al. 2007b, ApJ, 663, 365 In the text
- Zacharias, N., Urban, S. E., Zacharias, M. I., et al. 2004, AJ, 127, 3043
- Zuckerman, B., & Becklin, E. E. 1993, ApJ, 414, 793
- Zuckerman, B., & Song, I. 2004a, ApJ, 603, 738 In the text
- Zuckerman, B., & Song, I. 2004b, ARA&A, 42, 685
- Zwahlen, N., North, P., Debernardi, Y., et al. 2004, A&A, 425, L45
Footnotes
- ... stars
- This work is based on observations made with the IRAM (Institut de Radioastonomie Millimétrique) 30-m telescope and the CSO (Caltech Submillimetre Observatory) 10-m telescope. IRAM is supported by INSU/CNRS (France), MPG (Germany), and IGN (Spain).
- ... reference
- Note however, that the assumption of optically thin emission at millimetre wavelength cannot be generalised to all primordial discs
All Tables
Table 1:
The stellar properties, e.g. spectral type, temperature,
luminosity and age range are reported in the first columns (from Meyer et al. 2006, 2008, Hillenbrand et al. 2008). CSO 350 m and IRAM 1.2 mm measured fluxes and upper limits are reported in the sixth and seventh column. The detected sources are highlighted in bold. In the ninth column are the references to the stellar distances (listed in the eighth column).
Table 2:
Ages (in logarithmic scale) of the debris discs detected at the
millimetre wavelength, with IRAM (this work), SEST and OVRO (from C05).
They have been derived using different techniques: the strength of the
CaII H and K (
)
and Lithium (Li) emission, X-ray activity
(X-ray), and the rotational period-age (
)
relation
(L. Hillenbrand, private communication). The ages are averaged and the
standard deviation is adopted as error.
Table 3:
Blowout grain size, P-R drag timescales of the debris discs detected in the CSO and/or IRAM surveys. The -index discs only for the stars
with the SED sampled in the submillimetre/millimetre range
(details in Sect. 4) is reported.
Table 4:
Model results, reduced
(
)
assuming a fixed outer radius (
)
of 150 AU and a maximum grain size (
)
of 3 mm.
Table 5:
Model results leaving as free parameters using
,
and
for 6 different outer radii: 1000, 500, 300, 150, 100 and 50 AU.
All Figures
![]() |
Figure 1: Disc mass versus age of the system. Filled circles are the masses determined from the IRAM detections (this work) and from the detections presented by C05. Upper limits are shown as filled triangles (new IRAM observations), empty and filled triangles (objects measured by C05 and re-observed here with higher sensitivity) and down arrows (from C05). The sources observed in both studies are connected by dashed lines. The vertical line separates the sources younger than 20 Myr and evolved systems older than 20 Myr. |
Open with DEXTER | |
In the text |
![]() |
Figure 2:
SED of HD 104860 (upper plot a)) detections are
shown as empty diamonds while 3 |
Open with DEXTER | |
In the text |
![]() |
Figure 3: As Fig. 2 but for HD 8907. |
Open with DEXTER | |
In the text |
![]() |
Figure 4: As Fig. 2 but for HD 377. |
Open with DEXTER | |
In the text |
![]() |
Figure 5: As Fig. 2 but for HD 107146. |
Open with DEXTER | |
In the text |
![]() |
Figure 6: As Fig. 2 but for HD 61005. |
Open with DEXTER | |
In the text |
![]() |
Figure 7: As Fig. 2 but for HD191089. |
Open with DEXTER | |
In the text |
Copyright ESO 2009
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.