The collisional evolution of debris disks

András Gáspár, George H. Rieke, Zoltán Balog

Research output: Contribution to journalArticle

56 Citations (Scopus)

Abstract

We explore the collisional decay of disk mass and infrared emission in debris disks. With models, we show that the rate of the decay varies throughout the evolution of the disks, increasing its rate up to a certain point, which is followed by a leveling off to a slower value. The total disk mass falls off ∝t -0.35 at its fastest point (where t is time) for our reference model, while the dust mass and its proxy - the infrared excess emission - fades significantly faster (∝t -0.8). These later level off to a decay rate of M tot (t) t -0.08 and M dust (t) or L ir (t) t -0.6. This is slower than the ∝t -1 decay given for all three system parameters by traditional analytic models. We also compile an extensive catalog of Spitzer and Herschel 24, 70, and 100 μm observations. Assuming a log-normal distribution of initial disk masses, we generate model population decay curves for the fraction of stars harboring debris disks detected at 24 μm. We also model the distribution of measured excesses at the far-IR wavelengths (70-100 μm) at certain age regimes. We show general agreement at 24 μm between the decay of our numerical collisional population synthesis model and observations up to a Gyr. We associate offsets above a Gyr to stochastic events in a few select systems. We cannot fit the decay in the far-infrared convincingly with grain strength properties appropriate for silicates, but those of water ice give fits more consistent with the observations (other relatively weak grain materials would presumably also be successful). The oldest disks have a higher incidence of large excesses than predicted by the model; again, a plausible explanation is very late phases of high dynamical activity around a small number of stars. Finally, we constrain the variables of our numerical model by comparing the evolutionary trends generated from the exploration of the full parameter space to observations. Amongst other results, we show that erosive collisions are dominant in setting the timescale of the evolution and that planetesimals on the order of 100 km in diameter are necessary in the cascades for our population synthesis models to reproduce the observations.

Original languageEnglish (US)
Article number25
JournalAstrophysical Journal
Volume768
Issue number1
DOIs
StatePublished - May 1 2013

Fingerprint

debris
decay
dust
planetesimal
stars
protoplanets
leveling
synthesis
normal density functions
decay rates
catalogs
silicate
collision
silicates
cascades
ice
incidence
wavelength
timescale
trends

Keywords

  • circumstellar matter
  • infrared: stars
  • methods: numerical
  • planetary systems

ASJC Scopus subject areas

  • Space and Planetary Science
  • Astronomy and Astrophysics

Cite this

The collisional evolution of debris disks. / Gáspár, András; Rieke, George H.; Balog, Zoltán.

In: Astrophysical Journal, Vol. 768, No. 1, 25, 01.05.2013.

Research output: Contribution to journalArticle

Gáspár, András ; Rieke, George H. ; Balog, Zoltán. / The collisional evolution of debris disks. In: Astrophysical Journal. 2013 ; Vol. 768, No. 1.
@article{f59a1830b1794253b01917fe100d0fd2,
title = "The collisional evolution of debris disks",
abstract = "We explore the collisional decay of disk mass and infrared emission in debris disks. With models, we show that the rate of the decay varies throughout the evolution of the disks, increasing its rate up to a certain point, which is followed by a leveling off to a slower value. The total disk mass falls off ∝t -0.35 at its fastest point (where t is time) for our reference model, while the dust mass and its proxy - the infrared excess emission - fades significantly faster (∝t -0.8). These later level off to a decay rate of M tot (t) t -0.08 and M dust (t) or L ir (t) t -0.6. This is slower than the ∝t -1 decay given for all three system parameters by traditional analytic models. We also compile an extensive catalog of Spitzer and Herschel 24, 70, and 100 μm observations. Assuming a log-normal distribution of initial disk masses, we generate model population decay curves for the fraction of stars harboring debris disks detected at 24 μm. We also model the distribution of measured excesses at the far-IR wavelengths (70-100 μm) at certain age regimes. We show general agreement at 24 μm between the decay of our numerical collisional population synthesis model and observations up to a Gyr. We associate offsets above a Gyr to stochastic events in a few select systems. We cannot fit the decay in the far-infrared convincingly with grain strength properties appropriate for silicates, but those of water ice give fits more consistent with the observations (other relatively weak grain materials would presumably also be successful). The oldest disks have a higher incidence of large excesses than predicted by the model; again, a plausible explanation is very late phases of high dynamical activity around a small number of stars. Finally, we constrain the variables of our numerical model by comparing the evolutionary trends generated from the exploration of the full parameter space to observations. Amongst other results, we show that erosive collisions are dominant in setting the timescale of the evolution and that planetesimals on the order of 100 km in diameter are necessary in the cascades for our population synthesis models to reproduce the observations.",
keywords = "circumstellar matter, infrared: stars, methods: numerical, planetary systems",
author = "Andr{\'a}s G{\'a}sp{\'a}r and Rieke, {George H.} and Zolt{\'a}n Balog",
year = "2013",
month = "5",
day = "1",
doi = "10.1088/0004-637X/768/1/25",
language = "English (US)",
volume = "768",
journal = "Astrophysical Journal",
issn = "0004-637X",
publisher = "IOP Publishing Ltd.",
number = "1",

}

TY - JOUR

T1 - The collisional evolution of debris disks

AU - Gáspár, András

AU - Rieke, George H.

AU - Balog, Zoltán

PY - 2013/5/1

Y1 - 2013/5/1

N2 - We explore the collisional decay of disk mass and infrared emission in debris disks. With models, we show that the rate of the decay varies throughout the evolution of the disks, increasing its rate up to a certain point, which is followed by a leveling off to a slower value. The total disk mass falls off ∝t -0.35 at its fastest point (where t is time) for our reference model, while the dust mass and its proxy - the infrared excess emission - fades significantly faster (∝t -0.8). These later level off to a decay rate of M tot (t) t -0.08 and M dust (t) or L ir (t) t -0.6. This is slower than the ∝t -1 decay given for all three system parameters by traditional analytic models. We also compile an extensive catalog of Spitzer and Herschel 24, 70, and 100 μm observations. Assuming a log-normal distribution of initial disk masses, we generate model population decay curves for the fraction of stars harboring debris disks detected at 24 μm. We also model the distribution of measured excesses at the far-IR wavelengths (70-100 μm) at certain age regimes. We show general agreement at 24 μm between the decay of our numerical collisional population synthesis model and observations up to a Gyr. We associate offsets above a Gyr to stochastic events in a few select systems. We cannot fit the decay in the far-infrared convincingly with grain strength properties appropriate for silicates, but those of water ice give fits more consistent with the observations (other relatively weak grain materials would presumably also be successful). The oldest disks have a higher incidence of large excesses than predicted by the model; again, a plausible explanation is very late phases of high dynamical activity around a small number of stars. Finally, we constrain the variables of our numerical model by comparing the evolutionary trends generated from the exploration of the full parameter space to observations. Amongst other results, we show that erosive collisions are dominant in setting the timescale of the evolution and that planetesimals on the order of 100 km in diameter are necessary in the cascades for our population synthesis models to reproduce the observations.

AB - We explore the collisional decay of disk mass and infrared emission in debris disks. With models, we show that the rate of the decay varies throughout the evolution of the disks, increasing its rate up to a certain point, which is followed by a leveling off to a slower value. The total disk mass falls off ∝t -0.35 at its fastest point (where t is time) for our reference model, while the dust mass and its proxy - the infrared excess emission - fades significantly faster (∝t -0.8). These later level off to a decay rate of M tot (t) t -0.08 and M dust (t) or L ir (t) t -0.6. This is slower than the ∝t -1 decay given for all three system parameters by traditional analytic models. We also compile an extensive catalog of Spitzer and Herschel 24, 70, and 100 μm observations. Assuming a log-normal distribution of initial disk masses, we generate model population decay curves for the fraction of stars harboring debris disks detected at 24 μm. We also model the distribution of measured excesses at the far-IR wavelengths (70-100 μm) at certain age regimes. We show general agreement at 24 μm between the decay of our numerical collisional population synthesis model and observations up to a Gyr. We associate offsets above a Gyr to stochastic events in a few select systems. We cannot fit the decay in the far-infrared convincingly with grain strength properties appropriate for silicates, but those of water ice give fits more consistent with the observations (other relatively weak grain materials would presumably also be successful). The oldest disks have a higher incidence of large excesses than predicted by the model; again, a plausible explanation is very late phases of high dynamical activity around a small number of stars. Finally, we constrain the variables of our numerical model by comparing the evolutionary trends generated from the exploration of the full parameter space to observations. Amongst other results, we show that erosive collisions are dominant in setting the timescale of the evolution and that planetesimals on the order of 100 km in diameter are necessary in the cascades for our population synthesis models to reproduce the observations.

KW - circumstellar matter

KW - infrared: stars

KW - methods: numerical

KW - planetary systems

UR - http://www.scopus.com/inward/record.url?scp=84876777211&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84876777211&partnerID=8YFLogxK

U2 - 10.1088/0004-637X/768/1/25

DO - 10.1088/0004-637X/768/1/25

M3 - Article

AN - SCOPUS:84876777211

VL - 768

JO - Astrophysical Journal

JF - Astrophysical Journal

SN - 0004-637X

IS - 1

M1 - 25

ER -