Cold Debris Disk Research Articles

The MIRI Exoplanets Orbiting White dwarfs (MEOW) Survey: Mid-infrared Excess Reveals a Giant Planet Candidate around a Nearby White Dwarf

The MIRI Exoplanets Orbiting White dwarfs survey is a cycle 2 JWST program to search for exoplanets around dozens of nearby white dwarfs via infrared excess and direct imaging. In this Letter, we present the detection of mid-infrared excess at 18 and 21 μm toward the bright (V = 11.4) metal-polluted white dwarf WD 0310–688. The source of the IR excess is almost certainly within the system; the probability of background contamination is <0.1%. While the IR excess could be due to an unprecedentedly small and cold debris disk, it is best explained by a 3.0−1.9+5.5 M Jup cold (248 −61+84 K) giant planet orbiting the white dwarf within the forbidden zone (the region where planets are expected to be destroyed during the star’s red giant phase). We constrain the source of the IR excess to an orbital separation of 0.1–2 au, marking the first discovery of a white dwarf planet candidate within this range of separations. WD 0310–688 is a young remnant of an A- or late B-type star, and at just 10.4 pc, it is now the closest white dwarf with a known planet candidate. Future JWST observations could distinguish the two scenarios by either detecting or ruling out spectral features indicative of a planet atmosphere.

Planet populations inferred from debris discs

We know little about the outermost exoplanets in planetary systems because our detection methods are insensitive to moderate-mass planets on wide orbits. However, debris discs can probe the outer-planet population because dynamical modelling of observed discs can reveal properties of perturbing planets. We use four sculpting and stirring arguments to infer planet properties in 178 debris-disc systems from the ISPY, LEECH, and LIStEN planet-hunting surveys. Similar analyses are often conducted for individual discs, but we consider a large sample in a consistent manner. We aim to predict the population of wide-separation planets, gain insight into the formation and evolution histories of planetary systems, and determine the feasibility of detecting these planets in the near future. We show that a ‘typical’ cold debris disc likely requires a Neptune- to Saturn-mass planet at 10–100 au, with some needing Jupiter-mass perturbers. Our predicted planets are currently undetectable, but modest detection-limit improvements (e.g. from JWST) should reveal many such perturbers. We find that planets thought to be perturbing debris discs at late times are similar to those inferred to be forming in protoplanetary discs, so these could be the same population if newly formed planets do not migrate as far as currently thought. Alternatively, young planets could rapidly sculpt debris before migrating inwards, meaning that the responsible planets are more massive (and located farther inwards) than debris-disc studies assume. We combine self-stirring and size-distribution modelling to show that many debris discs cannot be self-stirred without having unreasonably high masses; planet- or companion-stirring may therefore be the dominant mechanism in many (perhaps all) debris discs. Finally, we provide catalogues of planet predictions and identify promising targets for future planet searches.

From Pebbles and Planetesimals to Planets and Dust: The Protoplanetary Disk–Debris Disk Connection

Abstract The similar orbital distances and detection rates of debris disks and the prominent rings observed in protoplanetary disks suggest a potential connection between these structures. We explore this connection with new calculations that follow the evolution of rings of pebbles and planetesimals as they grow into planets and generate dusty debris. Depending on the initial solid mass and planetesimal formation efficiency, the calculations predict diverse outcomes for the resulting planet masses and accompanying debris signature. When compared with debris disk incidence rates as a function of luminosity and time, the model results indicate that the known population of bright cold debris disks can be explained by rings of solids with the (high) initial masses inferred for protoplanetary disk rings and modest planetesimal formation efficiencies that are consistent with current theories of planetesimal formation. These results support the possibility that large protoplanetary disk rings evolve into the known cold debris disks. The inferred strong evolutionary connection between protoplanetary disks with large rings and mature stars with cold debris disks implies that the remaining majority population of low-mass stars with compact protoplanetary disks leaves behind only modest masses of residual solids at large radii and evolves primarily into mature stars without detectable debris beyond 30 au. The approach outlined here illustrates how combining observations with detailed evolutionary models of solids strongly constrains the global evolution of disk solids and underlying physical parameters such as the efficiency of planetesimal formation and the possible existence of invisible reservoirs of solids in protoplanetary disks.

Detection of H2 in the TWA 7 System: A Probable Circumstellar Origin

Abstract Using HST–COS far-UV spectra, we have discovered warm molecular hydrogen in the TWA 7 system. TWA 7, a ∼9 Myr old M2.5 star, has a cold debris disk and has previously shown no signs of accretion. Molecular hydrogen is expected to be extremely rare in a debris disk. While molecular hydrogen can be produced in starspots or the lower chromospheres of cool stars such as TWA 7, fluxes from progressions that get pumped by the wings of Lyα indicate that this molecular hydrogen could be circumstellar and thus that TWA 7 is accreting at very low levels and may retain a reservoir of gas in the near circumstellar environment.

Multi-wavelength, spatially resolved modelling of HD 48682’s debris disc

ABSTRACT Asteroids and comets (planetesimals) are created in gas- and dust-rich protoplanetary discs. The presence of these planetesimals around main-sequence stars is usually inferred from the detection of excess continuum emission at infrared wavelengths from dust grains produced by destructive processes within these discs. Modelling of the disc structure and dust grain properties for those discs is often hindered by the absence of any meaningful constraint on the location and spatial extent of the disc. Multi-wavelength, spatially resolved imaging is thus invaluable in refining the interpretation of these systems. Observations of HD 48682 at far-infrared (Spitzer, Herschel) and sub-millimetre (JCMT, SMA) wavelengths indicated the presence of an extended, cold debris disc with a blackbody temperature of 57.9 ± 0.7 K. Here, we combined these data to perform a comprehensive study of the disc architecture and its implications for the dust grain properties. The deconvolved images revealed a cold debris belt, verified by combining a 3D radiative transfer dust continuum model with image analysis to replicate the structure using a single, axisymmetric annulus. A Markov chain Monte Carlo analysis calculated the maximum likelihood of HD48682’s disc radius ($R_{\rm disc} = 89^{+17}_{-20}~$ au), fractional width ($\Delta R_{\rm disc} = 0.41^{+0.27}_{-0.20}$), position angle ($\theta = 66{_{.}^{\circ}} 3^{+4.5}_{-4.9}$), and inclination ($\phi = 112{_{.}^{\circ}} 5^{+4.2}_{-4.2}$). HD 48682 has been revealed to host a collisionally active, broad disc whose emission is dominated by small dust grains, smin ∼ 0.6 μm, and a size distribution exponent of 3.60 ± 0.02.

Twisted debris: how differential secular perturbations shape debris disks

Context. Resolved images suggest that asymmetric structures are a common feature of cold debris disks. While planets close to these disks are rarely detected, their hidden presence and gravitational perturbations provide plausible explanations for some of these features. Aims. To put constraints on the properties of yet undetected planetary companions, we aim to predict what features such a planet imprints in debris disks undergoing continuous collisional evolution. Methods. We discuss the basic equations, analytic approximations and timescales governing collisions, radiation pressure and secular perturbations. In addition, we combine our numerical model of the collisional evolution of the size and spatial distributions in debris disks with the gravitational perturbation by a single planet. Results. We find that the distributions of orbital elements in the disks are strongly dependent on grain sizes. Secular precession is differential with respect to involved semi-major axes and grain sizes. This leads to observable differences between the big grains tracing the parent belt and the small grains in the trailing halo. Observations at different wavelengths can be used to constrain the properties of a possible planet.

Constraints on HD 113337 fundamental parameters and planetary system

Context. HD 113337 is a main-sequence F6V field star more massive than the Sun. This star hosts one confirmed giant planet and possibly a second candidate, detected by radial velocities (RVs). The star also hosts a cold debris disc detected through the presence of an infrared excess, making it an interesting system to explore. Aims. We aim to bring new constraints on the star’s fundamental parameters, debris disc properties, and planetary companion(s) by combining complementary techniques. Methods. We used the VEGA interferometer on the CHARA array to measure the angular diameter of HD 113337. We derived its linear radius using the parallax from the Gaia Second Data Release. We computed the bolometric flux to derive its effective temperature and luminosity, and we estimated its mass and age using evolutionary tracks. Then, we used Herschel images to partially resolve the outer debris disc and estimate its extension and inclination. Next, we acquired high-contrast images of HD 113337 with the LBTI to probe the ~10–80 au separation range. Finally, we combined the deduced contrast maps with previous RVs of the star using the MESS2 software to bring upper mass limits on possible companions at all separations up to 80 au. We took advantage of the constraints on the age and inclination brought by fundamental parameter analysis and disc imaging, respectively, for this analysis. Results. We derive a limb-darkened angular diameter of 0.386 ± 0.009 mas that converts into a linear radius of 1.50 ± 0.04 R⊙ for HD 113337. The fundamental parameter analysis leads to an effective temperature of 6774 ± 125 K and to two possible age solutions: one young within 14–21 Myr and one old within 0.8–1.7 Gyr. We partially resolve the known outer debris disc and model its emission. Our best solution corresponds to a radius of 85 ± 20 au, an extension of 30 ± 20 au, and an inclination within 10–30° for the outer disc. The combination of imaging contrast limits, published RV, and age and inclination solutions allows us to derive a first possible estimation of the true masses of the planetary companions: ~7−2+4 MJup for HD 113337 b (confirmed companion) and ~16−3+10 MJup for HD 113337 c (candidate companion). We also constrain possible additional companions at larger separations.

Giant Planets Can Act as Stabilizing Agents on Debris Disks

Abstract We have explored the evolution of a cold debris disk under the gravitational influence of dwarf-planet-sized objects (DPs), both in the presence and absence of an interior giant planet. Through detailed long-term numerical simulations, we demonstrate that when the giant planet is not present, DPs can stir the eccentricities and inclinations of disk particles, in linear proportion to the total mass of the DPs; on the other hand, when the giant planet is included in the simulations, the stirring is approximately proportional to the mass squared. This creates two regimes: below a disk mass threshold (defined by the total mass of DPs), the giant planet acts as a stabilizing agent of the orbits of cometary nuclei, diminishing the effect of the scatterers; above the threshold, the giant contributes to the dispersion of the particles.

Exozodiacal clouds: hot and warm dust around main sequence stars

AbstractA warm/hot dust component (at temperature 300 K) has been detected around 20% of A, F, G, K stars. This component is called ‘exozodiacal dust’ as it presents similarities with the zodiacal dust detected in our solar system, even though its physical properties and spatial distribution can be significantly different. Understanding the origin and evolution of this dust is of crucial importance, not only because its presence could hamper future detections of Earth-like planets in their habitable zones, but also because it can provide invaluable information about the inner regions of planetary systems. In this review, we present a detailed overview of the observational techniques used in the detection and characterisation of exozodiacal dust clouds (‘exozodis’) and the results they have yielded so far, in particular regarding the incidence rate of exozodis as a function of crucial parameters such as stellar type and age, or the presence of an outer cold debris disc. We also present the important constr...

Inner mean-motion resonances with eccentric planets: a possible origin for exozodiacal dust clouds

High levels of dust have been detected in the immediate vicinity of many stars, both young and old. A promising scenario to explain the presence of this short-lived dust is that these analogues to the Zodiacal cloud (or exozodis) are refilled in situ through cometary activity and sublimation. As the reservoir of comets is not expected to be replenished, the presence of these exozodis in old systems has yet to be adequately explained. It was recently suggested that mean-motion resonances (MMR) with exterior planets on moderately eccentric ($\mathrm{e_p}\gtrsim 0.1$) orbits could scatter planetesimals on to cometary orbits with delays of the order of several 100 Myr. Theoretically, this mechanism is also expected to sustain continuous production of active comets once it has started, potentially over Gyr-timescales. We aim here to investigate the ability of this mechanism to generate scattering on to cometary orbits compatible with the production of an exozodi on long timescales. We combine analytical predictions and complementary numerical N-body simulations to study its characteristics. We show, using order of magnitude estimates, that via this mechanism, low mass discs comparable to the Kuiper Belt could sustain comet scattering at rates compatible with the presence of the exozodis which are detected around Solar-type stars, and on Gyr timescales. We also find that the levels of dust detected around Vega could be sustained via our proposed mechanism if an eccentric Jupiter-like planet were present exterior to the system's cold debris disc.

How to design a planetary system for different scattering outcomes: giant impact sweet spot, maximizing exocomets, scattered discs

This paper considers the dynamics of scattering of planetesimals or planetary embryos by a planet on a circumstellar orbit. We classify six regions in the planet's mass versus semimajor axis parameter space according to the dominant outcome for scattered objects: ejected, accreted, remaining, escaping, Oort Cloud, depleted Oort Cloud. We use these outcomes to consider which planetary system architectures maximise the observability of specific signatures, given that signatures should be detected first around systems with optimal architectures (if such systems exist in nature). Giant impact debris is most readily detectable for 0.1-10Mearth planets at 1-5au, depending on detection method and spectral type. While A stars have putative giant impact debris at 4-6au consistent with this sweet spot, that of FGK stars is typically <<1au contrary to expectations; an absence of 1-3au giant impact debris could indicate a low frequency of terrestrial planets there. Three principles maximise cometary influx from exo-Kuiper belts: a chain of closely separated planets interior to the belt, none of which is a Jupiter-like ejector; planet masses not increasing strongly with distance (for a net inward torque on comets); ongoing replenishment of comets, possibly by embedded low-mass planets. A high Oort Cloud comet influx requires no ejector and architectures that maximise the Oort Cloud population. Cold debris disks are usually considered classical Kuiper belt analogues. Here we consider the possibility of detecting scattered disk analogues, which could be betrayed by a broad radial profile and lack of small grains, as well as spherical 100-1000au mini-Oort Clouds. Some implications for escaping planets around young stars, detached planets akin to Sedna, and the formation of super-Earths, are also discussed.

DETECTION AND CHARACTERIZATION OF EXTRASOLAR PLANETS THROUGH MEAN-MOTION RESONANCES. I. SIMULATIONS OF HYPOTHETICAL DEBRIS DISKS

ABSTRACT The gravitational influence of a planet on a nearby disk provides a powerful tool for detecting and studying extrasolar planetary systems. Here we demonstrate that gaps can be opened in dynamically cold debris disks at the mean-motion resonances of an orbiting planet. The gaps are opened away from the orbit of the planet itself, revealing that not all disk gaps need contain a planetary body. These gaps are large and deep enough to be detectable in resolved disk images for a wide range of reasonable disk-planet parameters, though we are not aware of any such gaps detected to date. The gap shape and size are diagnostic of the planet location, eccentricity and mass, and allow one to infer the existence of unseen planets, as well as many important parameters of both seen and unseen planets in these systems. We present expressions to allow the planetary mass and semimajor axis to be calculated from observed gap width and location.

Herschel/PACS photometry of transiting-planet host stars with candidate warm debris disks

Dust in debris disks is produced by colliding or evaporating planetesimals, remnants of the planet formation process. Warm dust disks, known by their emission at < 24 micron, are rare (4% of FGK main sequence stars) and especially interesting because they trace material in the region likely to host terrestrial planets, where the dust has a very short dynamical lifetime. Statistical analyses of the source counts of excesses as found with the mid-IR Wide Field Infrared Survey Explorer (WISE) suggest that warm-dust candidates found for the Kepler transiting-planet host-star candidates can be explained by extragalactic or galactic background emission aligned by chance with the target stars. These statistical analyses do not exclude the possibility that a given WISE excess could be due to a transient dust population associated with the target. Here we report Herschel/PACS 100 and 160 micron follow-up observations of a sample of Kepler and non-Kepler transiting-planet candidates' host stars, with candidate WISE warm debris disks, aimed at detecting a possible cold debris disk in any of them. No clear detections were found in any one of the objects at either wavelength. Our upper limits confirm that most objects in the sample do not have a massive debris disk like that in beta Pic. We also show that the planet-hosting star WASP-33 does not have a debris disk comparable to the one around eta Crv. Although the data cannot be used to rule out rare warm disks around the Kepler planet-hosting candidates, the lack of detections and the characteristics of neighboring emission found at far-IR wavelengths support an earlier result suggesting that most of the WISE-selected IR excesses around Kepler candidate host stars are likely due to either chance alignment with background IR-bright galaxies and/or to interstellar emission.

Complex satellite systems: a general model of formation from rings

Summary of the AbstractWe present our model of formation of regular satellites by the spreading of a cold debris disk beyond the Roche radius. By numerical simulations and analytical calculations, we show that this process explains the peculiar properties of Saturn's small and mid-sized moons, in particular their mass - distance distribution. This process seems to also account for the structures of the satellite systems of Uranus, Neptune, and the Earth, suggesting that they used to have massive rings.

A dearth of small particles in debris disks

A prescription for the fragment size distribution resulting from dust grain collisions is essential when modelling a range of astrophysical systems, such as debris disks and planetary rings. While the slope of the fragment size distribution and the size of the largest fragment are well known, the behaviour of the distribution at the small size end is theoretically and experimentally poorly understood. This leads debris disk codes to generally assume a limit equal to, or below, the radiation blow-out size. We use energy conservation to analytically derive a lower boundary of the fragment size distribution for a range of collider mass ratios. Focussing on collisions between equal-sized bodies, we apply the method to debris disks. For a given collider mass, the size of the smallest fragments is found to depend on collision velocity, material parameters, and the size of the largest fragment. We provide a physically motivated recipe for the calculation of the smallest fragment, which can be easily implemented in codes for modelling collisional systems. For plausible parameters, our results are consistent with the observed predominance of grains much larger than the blow-out size in Fomalhaut's main belt and in the Herschel cold debris disks.

Vega's hot dust from icy planetesimals scattered inwards by an outward-migrating planetary system

Abstract Vega has been shown to host multiple dust populations, including both hot exozodiacal dust at sub-au radii and a cold debris disc extending beyond 100 au. We use dynamical simulations to show how Vega's hot dust can be created by long-range gravitational scattering of planetesimals from its cold outer regions. Planetesimals are scattered progressively inwards by a system of 5–7 planets from 30 to 60 au to very close-in. In successful simulations, the outermost planets are typically Neptune mass. The back-reaction of planetesimal scattering causes these planets to migrate outwards and continually interact with fresh planetesimals, replenishing the source of scattered bodies. The most favourable cases for producing Vega's exozodi have negative radial mass gradients, with sub-Saturn- to Jupiter-mass inner planets at 5–10 au and outer planets of 2.5 − 20 M⊕ . The mechanism fails if a Jupiter-sized planet exists beyond ∼15 au because the planet preferentially ejects planetesimals before they can reach the inner system. Direct-imaging planet searches can therefore directly test this mechanism.

Extrasolar planets and brown dwarfs around A-F type stars

In the frame of the search for extrasolar planets and brown dwarfs around early-type main-sequence stars, we present the detection of a giant planet around the young F-type star HD113337. We estimated the age of the system to be 150 +100/-50 Myr. Interestingly, an IR excess attributed to a cold debris disk was previously detected on this star. The SOPHIE spectrograph on the 1.93m telescope at Observatoire de Haute-Provence was used to obtain ~300 spectra over 6 years. We used our SAFIR tool, dedicated to the spectra analysis of A and F stars, to derive the radial velocity variations. The data reveal a 324.0 +1.7/-3.3 days period that we attribute to a giant planet with a minimum mass of 2.83 +- 0.24 Mjup in an eccentric orbit with e=0.46 +- 0.04. A long-term quadratic drift, that we assign to be probably of stellar origin, is superimposed to the Keplerian solution.

A TREND BETWEEN COLD DEBRIS DISK TEMPERATURE AND STELLAR TYPE: IMPLICATIONS FOR THE FORMATION AND EVOLUTION OF WIDE-ORBIT PLANETS

Cold debris disks trace the limits of planet formation or migration in the outer regions of planetary systems, and thus have the potential to answer many of the outstanding questions in wide-orbit planet formation and evolution. We characterized the infrared excess spectral energy distributions of 174 cold debris disks around 546 main-sequence stars observed by both Spitzer IRS and MIPS. We found a trend between the temperature of the inner edges of cold debris disks and the stellar type of the stars they orbit. This argues against the importance of strictly temperature-dependent processes (e.g. non-water ice lines) in setting the dimensions of cold debris disks. Also, we found no evidence that delayed stirring causes the trend. The trend may result from outward planet migration that traces the extent of the primordial protoplanetary disk, or it may result from planet formation that halts at an orbital radius limited by the efficiency of core accretion.

A Trend Between Cold Debris Disk Temperature and Stellar Type: Implications for the Formation and Evolution of Wide-Orbit Planets

AbstractCold debris disks have the potential to answer many outstanding questions in wide-orbit planet formation and evolution. We characterized the infrared excess SEDs of 174 cold debris disks with Spitzer IRS and MIPS. We found a trend between the temperature of the disks and the stellar type of the stars they orbit. This argues against the importance of strictly temperature-dependent processes (e.g. ice lines) in setting the dimensions of cold debris disks. We also found no evidence that delayed stirring causes the trend. The trend may result from outward planet migration that traces the extent of the primordial protoplanetary disk, or from planet formation that halts at an orbital radius limited by the efficiency of core accretion. For the full details of this work, see Ballering et al. (2013).

The nearby population of M-dwarfs with WISE: a search for warm circumstellar dust

Circumstellar debris disks are important because of their connection to planetary systems. An efficient way to identify these systems is through their infrared excess. Most studies so far concentrated on early-type or solar-type stars, but less effort has gone into investigating M dwarfs. We characterize the mid-infrared photometric behavior of M dwarfs and search for infrared excess in nearby M dwarfs taken from the volume-limited RECONS sample using data from the WISE satellite and the 2MASS catalog. Our sample consists of 85 sources encompassing 103 M dwarfs. We derive empirical infrared colors from these data and discuss their errors. Based on this, we check the stars for infrared excess and discuss the minimum excess we would be able to detect. Other than the M8.5 dwarf SCR 1845-6357 A, where the excess is produced by a known T6 companion, we detect no excesses in any of our sample stars. The limits we derive for the 22um excess are slightly higher than the usual detection limit of 10-15% for Spitzer studies, but including the [12]-[22] color in our analysis allows us to derive tight constraints on the fractional dust luminosity L_dust/L_star. We show that this result is consistent with M dwarf excesses in the mid-inrared being as frequent as excesses around earlier-type stars. The low detection rate could be an age effect. We also present a tentative excess detection at 22um around the known cold debris disk M dwarf AU Mic, which is not part of our statistical sample. There is still no clear detection of a mid-infrared excess around any old (>30 Myr) main-sequence M dwarf. It is unclear whether this is due to a different dust evolution around M dwarfs or whether this is an age effect combined with the diffculties involved in searching M dwarfs for infrared excesses. A significantly larger sample of well-studied M dwarfs is required to solve this question.