## Signs of a faint disc population at polluted white dwarfs

C. Bergfors 1 /star , J. Farihi 1 † , P. Dufour 2 and M. Rocchetto 1

D´ epartement de Physique, Universit´ e de Montr´al, Montr ´al, QC H3C 3J7, Canada e e

1 University College London, Department of Physics and Astronomy, Gower Street, London WC1E 6BT, UK 2

Accepted 2014 July 31. Received 2014 July 29; in original form 2014 June 5

## ABSTRACT

Observations of atmospheric metals and dust discs around white dwarfs provide important clues to the fate of terrestrial planetary systems around intermediate mass stars. We present Spitzer Infrared Array Camera observations of 15 metal polluted white dwarfs to investigate the occurrence and physical properties of circumstellar dust created by the disruption of planetary bodies. We find subtle infrared excess emission consistent with warm dust around KUV15519+1730 and HS2132+0941, and weaker excess around the DZ white dwarf G24558, which, if real, makes it the coolest white dwarf known to exhibit a 3 6 . µ m excess and the first DZ star with a bright disc. All together our data corroborate a picture where (1) discs at metal-enriched white dwarfs are commonplace and most escape detection in the infrared (possibly as narrow rings), (2) the discs are long lived, having lifetimes on the order of 10 6 yr or longer, and (3) the frequency of bright, infrared detectable discs decreases with age, on a timescale of roughly 500 Myr, suggesting large planetesimal disruptions decline on this same timescale.

Key words: circumstellar matter - planetary systems - white dwarfs

## 1 INTRODUCTION

Planetary systems around evolved stars and white dwarfs provide keys to understanding the history of planet formation in our Galaxy, and their fate - including that of our Solar System. Although no planetary mass objects have yet been discovered around single white dwarfs, several planets are known to orbit evolved stars and stellar remnants, e.g. pulsars (Wolszczan &amp; Frail 1992), binary white dwarf + neutron star systems (e.g. Sigurdsson et al. 2003), giants and subgiants (e.g. Frink et al. 2002; Johnson et al. 2007; Sato et al. 2008; Bowler et al. 2010). These discoveries, together with the fact that more than 95% of all stars are destined to become white dwarfs, imply that at least some planetary systems survive post-main sequence evolution. Based on the prevalence of white dwarfs that have metal polluted atmospheres deriving from the accretion of planetary material, at least ≈ 30 per cent still have rocky minor (or major) planets (Sion et al. 2009; Zuckerman et al. 2010; Koester, G¨ ansicke &amp; Farihi 2014).

mospheres consisting of pure hydrogen (DA type) or helium (nonDA) . Owing to their high surface gravity, the time-scale for any 1 metals to sink out of the atmosphere is short compared to the evolutionary time-scale for white dwarfs cooler than T eff /lessorsimilar 25 000 K, for which radiative forces are negligible (Koester 2009). The metal diffusion time-scale varies from only a few days for hot and young DAwhite dwarfs to a few 10 6 yr for helium and old hydrogen white dwarfs. Any heavy element enrichment in the atmospheres of cool white dwarfs therefore implies relatively recent accretion.

The presence and abundance pattern of metals in the atmospheres of white dwarfs, combined with often-observed infrared excesses, imply accretion of circumstellar material created by the disruption of minor planets within the stellar Roche limit. Since the first observation of infrared excess from a circumstellar disc around a (single) white dwarf (G29-38; Zuckerman &amp; Becklin 1987), clues to the origin of the circumstellar dust have accumulated, primarily from spectroscopic observations. Most white dwarfs have thin at-

- /star E-mail: c.bergfors@ucl.ac.uk

† Ernest Rutherford Fellow

The composition of the accreted material provides vital clues to its source. While early discoveries of polluted white dwarfs (e.g. Lacombe et al. 1983; Zuckerman et al. 2003) could not exclusively determine the origin of the accreted material - interstellar clouds, comets or asteroids have been suggested - the evidence has since become compelling that the accretion of planetary debris is the cause of the photospheric pollution and the only source that can adequately explain different observational properties. For instance, if accretion from the interstellar medium were responsible for the pollution it must have occurred relatively recently for DA stars. This theory is strongly disfavoured due to the lack of nearby dense interstellar clouds, and the relative absence of hydrogen in helium white dwarfs (see e.g. Farihi et al. 2010a).

The strongest support for accretion of disrupted planetesimals

1 White dwarfs are classified as D (degenerate) followed by identified optical or ultraviolet spectral features; in our sample we have types A (hydrogen) and B (helium), with Z denoting metal lines. DZ stars have helium-rich atmospheres but are too cool for He I absorption to be visible in their spectra (McCook&amp; Sion 1999).

comes from the inferred composition of the circumstellar material, which has been shown to be silicate rich and carbon poor - in contrast to the interstellar medium which is rich in volatiles such as carbon. The first infrared spectra with Spitzer of the dusty white dwarfs G29-38 and GD 362 (Reach et al. 2005a; Jura et al. 2007a) showed strong silicate emission at 9 -11 µ m consistent in shape with zodiacal dust emission, and no signs of polycyclic aromatic hydrocarbon (PAH) emission which is abundant in the interstellar medium at mid-infrared wavelengths. These features in the infrared spectra have since been confirmed for all white dwarfs with circumstellar discs that have been observed with the Spitzer Infrared Spectrograph (Jura, Farihi &amp; Zuckerman 2009), and additional silicate emission features were also found at 18 -20 µ m for G29-38 (Reach et al. 2009).

The metal abundances of the stars mirror the circumstellar material. Analysis of the heavily polluted atmosphere of GD 362 revealed 15 heavy elements, the relative abundances of which indicated that the debris had a composition that was similar to that of the Earth-Moon system (Zuckerman et al. 2007). In total 19 distinct metals have been detected in the atmospheres among several polluted white dwarfs from high resolution spectroscopic observations (e.g. Zuckerman et al. 2007; Melis et al. 2011; Dufour et al. 2012; Kawka &amp; Vennes 2012; Xu et al. 2014).

To date, 30 previously published metal polluted white dwarfs show excess infrared emission, and in some cases metallic gas emission (G¨ ansicke et al. 2006, 2008; G¨ ansicke, Marsh &amp; Southworth 2007; Farihi et al. 2012a; Melis et al. 2012), indicative of circumstellar debris. These discs are in general similar -the observed infrared excesses imply emission from hot circumstellar dust and a lack of cool dust outside of the stellar Roche limit ( r ≈ 1 R /circledot ). The now standard model for formation of these circumstellar discs involves a (as of yet undetected) planet perturbing a minor planetary body into a highly eccentric orbit (Debes &amp; Sigurdsson 2002; Bonsor, Mustill &amp; Wyatt 2011; Debes, Walsh &amp; Stark 2012) which is tidally disrupted when it comes within the Roche limit. The planetary material is viscously spread into a vertically thin, optically thick flat disc of temperature T ∼ 1000 K with an inner edge near a radius where silicate grains would rapidly sublimate (Jura 2003; Rafikov &amp; Garmilla 2012), from where it is accreted on to the white dwarf, causing the observed metal pollution. The composition of the accreted material therefore provides a way of measuring the bulk compositions of the disrupted minor or even major planets .

## 2 SCIENTIFIC MOTIVATION

To date, more than 200 white dwarfs have been searched for infrared excess, however dust emission from circumstellar discs has only been detected in those that have metal polluted atmospheres (von Hippel et al. 2007; Mullally et al. 2007; Farihi et al. 2012a). Previous studies have shown that while a significant fraction of polluted white dwarfs show infrared excess emission from circumstellar dust, the majority have no observable discs (Debes, Sigurdsson &amp; Hansen 2007; Farihi, Jura &amp; Zuckerman 2009). One possibility for the lack of observed disc emission around metal polluted white dwarfs is that the material is mainly gaseous from high speed collisions in discs with low particle densities (Jura, Farihi &amp; Zuckerman 2007b). Alternatively, for helium white dwarfs which have long diffusion time, it is possible that the metals detected in the photosphere have been accreted from discs that have since dissipated.

Jura (2008) proposed a model in which massive dust discs, and hence those with observable infrared emission, originate from the disruption of large asteroids, while for the majority of metal polluted white dwarfs continuous accretion of small asteroids creates a mainly gaseous disc with low surface dust density and hence modest infrared excess. In fact, previous studies have detected relatively subtle excesses around white dwarfs, which can be modelled as narrow dust rings (Farihi, Zuckerman &amp; Becklin 2008; Farihi et al. 2010b).

Observationally, the infrared bright disc frequency appears to be correlated with both the accretion rate (averaged over time for helium white dwarfs) and the white dwarf cooling age such that younger white dwarfs with high mass accretion rate are more prone to show infrared excess (Farihi et al. 2009; Xu &amp; Jura 2012). This correlation is however somewhat biased since it is observationally more challenging to detect metal pollution in the warmest white dwarfs compared to cooler ones, resulting in a smaller number of relatively young metal polluted stars observed with Spitzer . Discriminating between these two groups of polluted white dwarfs those with infrared detectable dust and those without - could possibly indicate whether the pollutant is one single large asteroid or a collection of several smaller ones (Jura 2008), with implications for the interpretation of the heavy element abundance patters in a spectral analysis.

We observed 15 metal polluted white dwarfs with Spitzer IRACin search for infrared excess indicative of warm circumstellar dust. The aim of this study is to extend previous surveys by observing additional metal polluted white dwarfs, and especially investigate bright disc frequency among cooler stars. In Section 3 we present the Spitzer IRAC and complementary ground based observations and data reduction. The results and analysis are presented in Section 4, including the detection of infrared excess consistent with warm circumstellar dust for two, possibly three, of our targets: KUV 15519+1730, G245-58 and HS 2132+0941. We discuss the properties of the newly discovered infrared excesses in the context of disc frequency and lifetime in Section 5, and summarise the results of the survey in Section 6.

## 3 OBSERVATIONS AND DATA REDUCTION

## 3.1 Target selection and Spitzer IRAC observations

Our 15 targets were selected in 2010 as the bulk of metal enhanced nearby and relatively bright white dwarfs that had not been observed with the Spitzer Infrared Array Camera (IRAC). The majority of the targets are moderately polluted with [Ca/H(e)] /lessorsimilar -7 5 . and mainly with effective temperatures T eff /lessorequalslant 10 000 K. These targets belong to a region of temperature and Ca abundance phase space that had not been well explored in previous studies. Three targets, KUV 15519+1730, HS 2132+0941 and HE 2230-1230, are warm ( T eff /greaterorsimilar 13 000 K) and highly metal polluted ([Ca/H(e)] /greaterorsimilar -7 5 . ) white dwarfs, thus likely to have observable excess infrared emission from circumstellar dust.

Warm Spitzer (Werner et al. 2004) infrared observations were carried out with IRAC (Fazio et al. 2004a) in Cycle 7 within programs 70037 and 70116. 30 medium size dithers of 30 s integrations were obtained in the cycling pattern for each target in IRAC channels 1 ( 3 6 . µ m) and 2 ( 4 5 . µ m). When combined to a 900 s integration mosaic the many dithers efficiently reduce artefacts such as

Table 1. List of IRAC imaging targets. Coordinates are measured from the Cycle 7 IRAC images (equinox J2000.0).

| WD       | Name           | R. A. (h ms)   | Dec ( ◦ ′ ′′ )   | SpT   | T eff ( K )   | [Ca/H(e)]   | Ref.   |
|----------|----------------|----------------|------------------|-------|---------------|-------------|--------|
| 0122-227 | KUV01223-2245  | 01 24 44.8     | - 22 29 07       | DZ    | 10 000        | - 9 . 4     | 1      |
| 0246+734 | G245-58        | 02 51 51.8     | +73 41 34        | DZ    | 7500          | - 9 . 7     | 2      |
| 0543+579 | GD 290         | 05 47 56.5     | +57 59 22        | DAZ   | 8140          | - 10 . 3    | 3      |
| 0625+100 | G105-B2B       | 06 27 37.6     | +10 02 13        | DZ    | 8800          | - 9 . 5     | 4      |
| 0816-310 | SCR0818 - 3110 | 08 18 40.5     | +31 10 29        | DZ    | 6460          | - 9 . 2     | 5      |
| 0840-136 | LP786-1        | 08 42 48.3     | - 13 47 13       | DZ    | 4870          | - 10 . 9    | 5      |
| 1009-184 | WT1759         | 10 12 01.5     | - 18 43 33       | DZ    | 6040          | - 10 . 6    | 5      |
| 1055-039 | G163-28        | 10 57 47.8     | - 04 13 32       | DZ    | 6500          | - 10 . 7    | 4      |
| 1338-311 | CE354          | 13 41 25.7     | - 31 24 50       | DZ    | 8210          | - 10 . 0    | 6      |
| 1425+540 | G200-40        | 14 27 35.8     | +53 48 30        | DBAZ  | 14 490        | - 9 . 3     | 7, 8   |
| 1551+175 | KUV15519+1730  | 15 54 09.1     | +17 21 24        | DBAZ  | 15 550        | - 6 . 5     | 8      |
| 1821-131 | GJ2135         | 18 24 04.4     | - 13 08 49       | DAZ   | 7030          | - 10 . 7    | 3      |
| 2132+096 | HS2132+0941    | 21 34 50.8     | +09 55 20        | DAZ   | 13 200        | - 7 . 7     | 9      |
| 2157-574 | WT853          | 22 00 45.1     | - 57 11 24       | DAZ   | 7220          | - 8 . 0     | 2, 10  |
| 2230-125 | HE2230 - 1230  | 22 33 38.8     | - 12 15 28       | DAZ   | 20 300        | - 6 . 3     | 9      |

References: (1) Friedrich et al. (2000); (2) this paper; (3) Zuckerman et al. (2003); (4) Dupuis, Fontaine &amp; Wesemael (1993); (5) Giammichele, Bergeron &amp; Dufour (2012); (6) Dufour et al. (2007); (7) Kenyon et al. (1988); (8) Bergeron et al. (2011); (9) Koester et al. (2005); (10) Subasavage et al. (2007)

cosmic rays and bad pixels. Table 1 lists the observed targets, their coordinates at the epoch of observation, spectral types, published temperatures and calcium abundances.

Table 2. IRAC Ch1 and Ch2 photometry.

## 3.2 Spitzer IRAC photometry

The IRAC basic calibrated data (BCD) reduction pipeline produces flux-calibrated data frames from which photometry with a few per cent accuracy can be performed either directly on the individual images or after mosaic combination. We combined the additional artefact mitigated, corrected BCD (CBCD) images to a create a coadded mosaic image using the MOPEX post-BCD software 2 with a pixel size 0.6 arcsec pixel -1 .

We perform photometry following Farihi et al. (2008). The main effects that contribute to increase photometric uncertainties arise from the array location of the target, and the pixel phase effect where the position of the point-spread function (PSF) centroid on the peak pixel affects the photometric result (Reach et al. 2005b; Hora et al. 2008). Both the array-location and pixel-phase effects are expected to average out for well dithered data such as ours. Colour corrections are assumed to be negligible.

Aperture and PSF-fitting photometry was performed on the mosaic combined images for each target and filter using the IRAF APPHOT and DAOPHOT packages. To avoid photometric confusion, we perform photometry either using a small aperture of radius r = 4 pixel and 24-40 pixel sky annulus, with the corrections recommended in the IRAC Instrument Handbook, or PSF-fitting on the mosaic combined images. We calculate 1 σ errors of the aperture photometry by measuring the standard deviation in the sky, and adopt the errors reported by DAOPHOT for the PSF-fitting photometry, adding in quadrature a 5 per cent absolute calibration er-

2 Details on the reduction pipeline and additional software can be found in the IRAC Instrument Handbook and related documentation at the Spitzer Science Center (SSC) website at http://irsa.ipac.caltech.edu/data/SPITZER/docs/irac/

| WD         | F 3 . 6 µm ( µ Jy)   | F 4 . 5 µm ( µ Jy)   | IR excess   |
|------------|----------------------|----------------------|-------------|
| 0122-227   | 52 ± 3               | 38 ± 2               | N           |
| 0246+734 a | 65 ± 5               | 47 ± 5               | Y           |
| 0543+579   | 172 ± 9              | 106 ± 6              | N           |
| 0625+100   | 82 ± 8               | 56 ± 7               | N           |
| 0816-310   | 382 ± 19             | 248 ± 13             | N           |
| 0840-136   | 550 ± 28             | 366 ± 18             | N           |
| 1009-184   | 474 ± 24             | 291 ± 15             | N           |
| 1055-039   | 129 ± 7              | 88 ± 5               | N           |
| 1338-311   | 48 ± 5               | 32 ± 4               | N           |
| 1425+540   | 163 ± 8              | 102 ± 5              | N           |
| 1551+175   | 29 ± 2               | 26 ± 2               | Y           |
| 1821-131 b | -                    | -                    | -           |
| 2132+096 a | 98 ± 5               | 79 ± 5               | Y           |
| 2157-574   | 252 ± 13             | 152 ± 8              | N           |
| 2230-125   | 53 ± 3               | 34 ± 2               | N           |

a Weighted average of fluxes extracted with DAOPHOT and APEX (see Section 3.4).

b Accurate fluxes could not be extracted for this target due to severe confusion.

ror to each measurement (see Farihi et al. 2008). Additional photometry was performed for a few targets using the Astronomical Point-source Extractor (APEX) in MOPEX to deblend some targets or confirm an infrared excess (see Section 3.4). The Spitzer IRAC photometry is listed in Table 2.

## 3.3 Ground based near-infrared photometry

Additional ground-based observations were obtained for some Spitzer IRAC targets in the near-infrared JHK s bands with the Long-slit Intermediate Resolution Infrared Spectrograph (LIRIS;

Figure 1. Spitzer IRAC 0 . ′′ 6 pixel -1 images in 3 6 . µ m(Ch 1) and 4 5 . µ m(Ch 2) of G245-58. The images are shown on identical square root scales. The target star appears slightly elongated at 4 5 . µ m, most likely due to a background object.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0229v2-figure-0c8df228-4b3b-43bf-8569-c18dd5d8f606.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0229v2-figure-c05a6293-dfb7-4087-9660-bdc5bb7d350a.jpg)

Manchado et al. 1998) at the 4.2 m William Herschel Telescope (WHT)in 2011 March and October, and with Son of ISAAC (SOFI; Moorwood, Cuby &amp; Lidman 1998) at the 3.6 m New Technology Telescope (NTT) in 2011 August and 2012 January. The science targets were observed in a continuous dither pattern with individual exposure times of 30, 15, and 10 s at J , H , and K s respectively for a total integration time of 270 s in each filter. Each night three standard star fields from the ARNICA catalogue (Hunt et al. 1998) were observed for flux calibration.

is thus likely associated with circumstellar dust around the white dwarf.

Standard data reduction was performed. The median sky subtracted frames were shifted and average combined into a single image on which aperture photometry with IRAF APPHOT or, for targets with nearby sources, DAOPHOT, was performed using standard aperture radii between 3 8 . and 4 3 . arcsec and sky annuli in the range 5 0 . -8 6 . arcsec. Smaller apertures were used for the faintest sources ( m&gt; 15 5 . mag), with corrections derived from bright stars in the same image. The flux calibration was good to 5 per cent or better on each night.

## 3.4 Photometric extraction of blended sources

G245-58 appears slightly elongated at 4 5 . µ m and therefore blending with a background object was suspected (Fig. 1). The additional source is brighter at 4 5 . µ m, where it appears relatively extended and diffuse, and hence it is likely a background galaxy.

In order to avoid contamination from the neighbouring source we performed PSF-fitting photometry on the science target using both DAOPHOT and APEX, in addition to aperture photometry, with the results listed in Table 3. A comparison between the aperture and PSF-fitting photometry showed that while the aperture photometry gives slightly higher flux densities, the values are consistent with those from the PSF-fitting photometry within errors. We also performed aperture photometry on the residual images generated by DAOPHOT after subtraction of the white dwarf, and found a residual flux density of 2 -4 per cent of the extracted target flux in both channels. This is well within the error bars and identical to that found for isolated targets analysed in the same manner. The potential contamination by the background source is therefore much smaller than the total flux excess measured compared to a white dwarf atmospheric model (see Section 4.3) and the infrared excess

As an additional test we compared the position relative to other stars in the field in the LIRIS JHK s and the 3 6 . and 4 5 . µ m IRAC images, since a shift of the PSF peak could indicate contamination of the photometry by the background source. We find no significant shift and conclude that the flux we extracted and measured in the IRAC images is that of the target. The PRF-fitting photometry with APEX gives marginally different values, and we adopt a weighted average from the two methods.

Based on our previous experience, approximately 1 out of 20 white dwarfs observed with IRAC needs some photometric deblending, but are typically recoverable (e.g. Farihi et al. 2010b). Our adopted flux measurement for G245-58 at 4 5 . µ m is 47 µ Jy, corresponding to 16.5 magnitudes (Vega). In order to assess the probability that the flux we measure is that of a background source and not the science target, we use IRAC 4 5 . µ m source counts from Fazio et al. (2004b). For the worst case scenario in which the galaxy contributes most of the flux measured at this wavelength, there are less than 2 × 10 4 galaxies deg -2 within a 1 mag interval centred at 17.0 mag. With a photometric aperture radius of 2 4 . arcsec, we find a probability of less than 1 per cent that the measured flux is due to a sufficiently bright background galaxy.

Two of our targets, WD 0625+100 and WD 1338-311, suffer from some blending with nearby stars. In both cases PSF-fitting photometry successfully de-convolved both point sources. No infrared excess is observed for these targets in IRAC, but flux contamination is the likely cause of excess seen in the Wide-field Infrared Survey Explorer ( WISE ; Wright et al. 2010) 3 4 . µ m flux of WD0625+100 (see Fig. 2). We also perform photometry with APEX and DAOPHOT for HS2132+0941 to confirm the infrared excess which is only seen at 4 5 . µ m. The photometry using both methods and the adopted values are listed in Table 3.

## 4 RESULTS AND ANALYSIS

## 4.1 Spectral Energy Distribution

Models of pure hydrogen and helium atmosphere white dwarfs (Koester 2010) for log g = 8 0 . were scaled to fit available photometric data in the optical from the Sloan Digital Sky Survey (SDSS;

Figure 2. Spectral energy distributions for target stars. Photometric data are from GALEX , SDSS, Subasavage et al. (2007, 2008), 2MASS, Table 4, WISE and IRAC (see Tables 2 and 4). The WISE 4 6 . µ mphotometry (light grey) have been slightly offset from the central wavelength in the plot for clarity.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0229v2-figure-c8c21aaa-4d1f-4dc6-b6fd-322e9b8b732f.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0229v2-figure-305e756a-883a-44cc-9120-dc0e2346099a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0229v2-figure-ea521837-b843-430a-b18a-a0e182baefaf.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0229v2-figure-ed793389-10f7-4f5d-b58e-7242252ab486.jpg)

Table 3. IRAF and APEX photometry of G245-58 and HS 2132+0941.

Table 4. Ground-based JHK s photometry. All errors are 5 per cent

| Wavelength ( µ m )   | APEX ( µ Jy )   | DAOPHOT ( µ Jy )   | Adopted flux ( µ Jy )   |
|----------------------|-----------------|--------------------|-------------------------|
| G245-58              | G245-58         | G245-58            | G245-58                 |
| 3 . 6                | 65 ± 4          | 64 ± 5             | 65 ± 5                  |
| 4 . 5                | 45 ± 3          | 50 ± 5             | 47 ± 5                  |
| HS2132+0941          | HS2132+0941     | HS2132+0941        | HS2132+0941             |
| 3 . 6                | 99 ± 5          | 96 ± 6             | 98 ± 5                  |
| 4 . 5                | 81 ± 4          | 77 ± 5             | 79 ± 5                  |

Stoughton et al. 2002), Subasavage et al. (2007, 2008), and in the near-infrared from Table 4 or the Two Micron All Sky Survey (2MASS; Skrutskie et al. 2006). We also plot in Fig. 2-4 ultraviolet photometry from the Galaxy Evolution Explorer ( GALEX ; Martin et al. 2005) if available, with 10 per cent error bars. The farultraviolet (FUV) and near-ultraviolet (NUV) data are however not

| WD       |   J (mag) |   H (mag]) |   K s (mag) | Instrument /Telescope   |
|----------|-----------|------------|-------------|-------------------------|
| 0122-227 |     16.83 |      16.82 |       16.83 | SOFI/NTT                |
| 0246+734 |     16.93 |      16.85 |       16.84 | LIRIS/WHT               |
| 0625+100 |     16.35 |      16.32 |       16.3  | LIRIS/WHT               |
| 1055-039 |     15.98 |      15.92 |       15.84 | SOFI/NTT                |
| 1338-311 |     17.03 |      16.99 |       17.06 | SOFI/NTT                |
| 1551+175 |     17.52 |      17.64 |       17.56 | SOFI/NTT                |
| 1821-131 |     14.51 |      14.28 |       14.21 | SOFI/NTT                |
| 2132+096 |     16.15 |      16.17 |       16.17 | SOFI/NTT                |

considered for the model fitting as the fluxes may be affected by metal blanketing and extinction (Koester et al. 2011). When fitting the models we assumed error bars of 5 per cent for the SDSS photometry, errors as listed in catalogue for 2MASS data (minimum 5 per cent), and 3 per cent error bars for Johnson V and Cousins R and I photometric data from Subasavage et al. (2007, 2008). Figs 2

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0229v2-figure-91b7dfad-0375-4ea8-8808-c5b0b4ea94d5.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0229v2-figure-d6c9b723-3ab8-4830-8c7b-27fe9e2b9ede.jpg)

Figure 2. - continued

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0229v2-figure-34363b8e-d3e2-44df-8e8d-91b433ff2ce7.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0229v2-figure-82a8e521-b972-4bbf-a391-467f23fd1a84.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0229v2-figure-5d0bbb8d-7861-460e-9457-a3e44ea709d4.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0229v2-figure-e3a4c43f-9225-49d1-9101-837ea5cdf47b.jpg)

and 3 show the spectral energy distributions of the targets without infrared excess. Photometry from WISE at 3 4 . and 4 6 . µ m is plotted for comparison in these figures if available; it is not however considered for model fitting. WISE data are insufficient in spatial resolution and sensitivity for these observational goals (Melis et al. 2011; Rocchetto et al., in preparation).

Two of the 15 targets show clear infrared excess compared to the white dwarf atmospheric model at 3 6 . and 4 5 . µ m: KUV15519+1730 shows a 4 σ excess at both 3 6 . and 4 5 . µ m, and a 3 σ excess at 2 2 . µ m, compared to the best model fit to the optical and near-infrared data. We detect a less significant ( 2 5 . σ ≈ 20 -25 per cent) excess in both wavelength bands for G245-58 in what is likely the most subtle disc emission detected at at white dwarf. HS2132+0941 exhibits a 4 σ excess at 4 5 . µ m compared to the expected photospheric flux, but no significant excess at shorter wavelengths. The spectral energy distributions of these targets are shown in Fig. 4.

We model the infrared emission of these three stars using opaque and geometrically flat dust rings (Jura 2003), which well describe the thermal dust emission observed at metal-rich white dwarfs. For each star the stellar radius to distance ratio R /d ∗ is calculated using the cooling models of the Montr´ eal group 3 for the adopted stellar atmospheric model temperature and assuming log g = 8 0 . (Holberg &amp; Bergeron 2006; Kowalski &amp; Saumon 2006; Bergeron et al. 2011; Tremblay, Bergeron &amp; Gianninas 2011), and three free parameters are varied to achieve a best fit to the near-infrared and IRAC data: the inner and outer disc temperature, T in and T out , and the inclination i of the disc. Corresponding disc radii are calculated following Chiang &amp; Goldreich (1997).

While the inner edge temperature can be relatively well constrained by the emission at shorter wavelengths ( 2 2 . and 3 6 . µ m), there is significant degeneracy in these models and they should be considered descriptive rather than definitive. For subtle excesses as found for these three targets, the thermal dust emission can be equally well fitted with relatively wide rings at high inclination, or as narrow rings at modest inclination (Girven et al. 2012). We therefore present two sets of models for each target with excess emission in Table 5 and Fig. 4. These model disc parameters do not represent the actual disc properties, but give representative, physically realistic disc temperatures consistent with the data. Importantly, all realistic models place these discs completely within the Roche limit, and thus consistent with disrupted planetesimals. The targets and their infrared excesses are discussed more in detail in the following subsections.

## 4.2 KUV15519+1730

This helium white dwarf shows traces of hydrogen in its spectrum together with Ca H &amp; K lines and is therefore classified as spectral type DBAZ (Wegner &amp; Swanson 1990; Bergeron et al. 2011). The excess in each of the 3 6 . and 4 5 . µ m wavelength bands is derived compared to a pure helium atmosphere white dwarf model of T eff = 15500 K. The star also shows a ≈ 3 σ excess in the K s band data obtained with SOFI at the NTT. The white dwarf is isolated and there is no indication of photometric contamination. Fig. 4 shows the disc spectral energy distribution for the two sets of models using the model input parameters listed in Table 5.

3 http://www.astro.umontreal.ca/ ∼ bergeron/CoolingModels/

Figure 3. (See Fig. 2 caption). The photometry for WD 0840-136 is best fitted with a T = 5750 K black body (Section 4.5).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0229v2-figure-1c6e22c1-c896-4b21-b558-47ff8c6ed6d9.jpg)

## 4.3 G245-58

Using a pure helium atmosphere model, the optical and nearinfrared photometry of this source is best fitted with T eff = 8250 K, and is shown in Fig. 4. Compared to this model, there is a 2 5 . σ infrared excess in both IRAC channels. We adopt the weighted average flux obtained with DAOPHOT and APEX for the 3 6 . and 4 5 . µ m IRAC photometry (see Section 3.4 and Table 3). The 2 5 . σ excess is derived using conservative absolute calibration errors of 5 per cent. The likelihood that this is a real excess is strengthened in two ways. First, Reach et al. (2005b) report an absolute IRAC calibration error of less than 3 per cent, which, if adopted here, would increase the infrared excess to &gt; σ 3 . Second, the excess is observed in both IRAC channels while the background source is more prominent at 4 5 . µ m.

Taken at face value, the best fitting disc model to the weak excess at G245-58 confines the circumstellar dust into a markedly narrow ring - even for the high inclination model, the ring is smaller than four white dwarf radii. If the excess proves to be real, this is the first DZ star discovered with infrared excess, and the coolest with excess 3 6 . µ m emission.

We also fit metal-rich models of Dufour et al. (2007) to the short wavelength photometry of this source and derive a slightly cooler effective temperature of T eff = 7500 K and [Ca/He]= -9 7 . (see Table 1). These values are used in the following discussion and in Fig. 5 and 6.

## 4.4 HS2132+0941

We fit a T eff = 13000 K atmospheric model to optical and near infrared JHK s photometry for this hydrogen white dwarf. The star shows a 4 σ infrared excess at 4 5 . µ m but no significant excess at shorter wavelengths (see Fig. 4). This is uncommon - thermal disc emission is in general notable at white dwarfs at least in the 3 6 . µ m wavelength band and sometimes also in K s . It is however not unique. Observations of the cool white dwarf G166-58 with Spitzer IRAC revealed excess emission at 5 8 . and 7 9 . µ m only (Farihi et al. 2008), corresponding to dust no warmer than 400 K. For PG1225-079, excess emission was observed only at 7.9 µ m (Farihi et al. 2010b). Since we only see excess flux in one wavelength band for HS 2132+0941, the outer radius of the disc can not

Figure 4. (See Fig. 2 caption.) Disc model fits to the stars with infrared excess. The dashed red and blue lines show the disc emission model for highly inclined and narrow rings respectively (see Section 4.1 and Table 5), while the dash-dotted lines give the total emission from stellar atmosphere and disc combined.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0229v2-figure-1b662702-720f-4672-9696-ec455196508e.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0229v2-figure-0e524fe4-cb22-454b-9a42-0b66d1038319.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0229v2-figure-ba3fd9e7-e1c2-4409-83cf-2eae4acbc652.jpg)

be well constrained. Until further observations are able to better constrain the nature of the excess, the best fitting disc model parameters are listed in Table 5.

## 4.5 WD0840-136

This cool DZ star has an estimated temperature of T eff ≈ 4900 K (Giammichele et al. 2012), while our coolest He atmosphere model has a temperature of T eff = 6000 K. We therefore fit a blackbody function to the photometric measurements in JHK s from 2MASS and in VRI from Subasavage et al. (2007), finding decent agreement for T = 5750 K which is shown in Fig 3. This illustrates the well known behaviour of He atmospheric models with metals (and thus free electrons), which generally yield colours similar to a metal-free (or blackbody) model of higher temperature.

## 5 DISCUSSION

## 5.1 Disc fraction as a function of cooling age

Metal polluted white dwarfs are common, yet only a fraction have observable disc emission. While the presence of metals in the atmospheres of helium white dwarfs may be explained as remnants of accretion from depleted discs due to the long diffusion times, the lack of discs around polluted hydrogen white dwarfs provides more of a puzzle.

Previous studies have shown that the fraction of polluted white dwarfs with discs observable as infrared excess is related to both metal accretion rate and effective temperature (i.e. cooling age) of the stars (e.g. von Hippel et al. 2007; Farihi et al. 2009; Xu &amp; Jura 2012). Fig. 5 depicts the time averaged mass accretion (Farihi et al. 2012b; Girven et al. 2012) as a function of effective temperature for all polluted white dwarfs observed by Spitzer IRAC. Adding our observations to already published data yield 104 stars, and corroborates trends noted in previous studies: infrared excess appears to be positively correlated with higher temperatures and younger ages, while there is a lack of detected excesses around cooler and older white dwarfs. However, because warmer white dwarfs have more opaque atmospheres than cooler ones, only the strongest metal absorption is detectable in their atmospheres from the ground, introducing a bias in the trend of accretion rate with temperature (von Hippel et al. 2007; Zuckerman et al. 2010). Koester et al. (2014) obtained spectra of 85 DA white dwarfs with the Cosmic Origins Spectrograph on the Hubble Space Telescope and, when combined with ground-based observations, found no trend in accretion rate over cooling ages ≈ 40 Myr to ≈ 2 Gyr, or T eff = 5000 -23 000 K.

Fig. 6 depicts the fraction of polluted white dwarfs that have been observed with Spitzer IRAC for which disc emission was detected to those without. This distribution is somewhat biased since only the strongest pollution is detected for the warmest white dwarfs. High metal pollution is likely correlated with the most massive and therefore easily detected discs, thereby increasing the disc fraction for warm, young white dwarfs. Even so, from Fig. 5 it is clear that there is a trend of decreasing disc detections with decreasing temperature, even if considering only the most heavily polluted white dwarfs with [Ca/H(e)] &gt; -7 5 . , consistent with previous findings (Farihi et al. 2009; Xu &amp; Jura 2012). The fraction of metal enhanced white dwarfs with observable discs drops sharply around T eff ≈ 10 000 K, or around 500 Myr, suggesting the disruption of large asteroids subsides on this timescale.

## 5.2 A few discs around old white dwarfs

Jura (2008) provided a model in which the lack of infrared excess emission among metal polluted white dwarfs is plausibly explained

Table 5. Representative disc parameters obtained by fitting a model with high inclination or narrow radial extent.

| WD                        | T eff ( K )               | R ∗ /d ( 10 - 12 )        | T in ( K )                | T out ( K )               | r in ( R WD )             | r out ( R WD )            | i (deg)                   |
|---------------------------|---------------------------|---------------------------|---------------------------|---------------------------|---------------------------|---------------------------|---------------------------|
| Model 1: High inclination | Model 1: High inclination | Model 1: High inclination | Model 1: High inclination | Model 1: High inclination | Model 1: High inclination | Model 1: High inclination | Model 1: High inclination |
| 0246+734                  | 8250                      | 3.73                      | 1200                      | 900                       | 7.8                       | 11.4                      | 80                        |
| 1551+175                  | 15500                     | 1.80                      | 1400                      | 900                       | 14.7                      | 26.5                      | 85                        |
| 2132+096                  | 13000                     | 3.84                      | 800                       | 400                       | 24.6                      | 61.9                      | 85                        |
| Model 2: Narrow rings     | Model 2: Narrow rings     | Model 2: Narrow rings     | Model 2: Narrow rings     | Model 2: Narrow rings     | Model 2: Narrow rings     | Model 2: Narrow rings     | Model 2: Narrow rings     |
| 0246+734                  | 8250                      | 3.73                      | 1000                      | 970                       | 9.9                       | 10.4                      | 60                        |
| 1551+175                  | 15500                     | 1.80                      | 1300                      | 1200                      | 16.2                      | 18.1                      | 60                        |
| 2132+096                  | 13000                     | 3.84                      | 600                       | 520                       | 36.0                      | 43.6                      | 60                        |

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0229v2-figure-40ed5cde-6850-4de8-a6ba-227cec1253b6.jpg)

Figure 5. Mass accretion rate as a function of effective temperature for white dwarfs with known Ca abundance that have been observed with Spitzer IRAC. Hydrogen white dwarfs with infrared excesses are plotted as red stars, and as black circles otherwise. Similarly, helium white dwarfs with infrared excess are plotted as blue stars, and as grey circles otherwise.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0229v2-figure-f00cbd3e-f458-4e9c-8eb4-250e32ab4f2d.jpg)

by accretion from tenuous or gaseous discs created as small asteroids enter a pre-existing disc at a non-zero inclination. In this model, erosion of the debris solids by sputtering creates a gaseous disc of vaporised asteroidal material, providing metals to be accreted by the host star but no observable infrared excess. Assuming a power-law mass distribution proportional to M -2 , similar to that of the solar system asteroid belt, implies that the majority of asteroids perturbed and destroyed within the tidal radius of the star are small, and massive discs may be expected to appear primarily for young systems in which larger asteroids (more massive than the pre-existing disc) are still relatively prevalent. In this picture, the observed trend of low detectable disc fraction for white dwarfs cooler than T eff /lessorsimilar 10 000 K, supports a scenario where large planetesimal disruptions become rare.

Before this study, only one white dwarf significantly cooler than 10 000 K was known to have an infrared excess: the hydrogen white dwarf G166-58 (Farihi et al. 2008). The excess emission of this star is unusual in being detectable beginning at 6 µ m. The possible discovery in this study of a weak infrared excess from a narrow dust ring around the cool DZ white dwarf G245-58, together with the subtle excess of G166-58, strengthens the idea that detectable discs become rare with age, but that there is a population of white

Figure 6. Fraction of polluted white dwarfs observed to have infrared excess with Spitzer IRAC. Light grey bins represent the total number of observed white dwarfs, dark grey those with observed disc emission. While the fraction of polluted white dwarfs with discs is higher for warmer and younger white dwarfs, the discoveries in this study add to the number of cooler white dwarfs with faint excesses, suggesting there is a population of white dwarfs which still retain dust at old age.

dwarfs that still retain bright dust discs at ages of ≈ 1 Gyr. The long metal diffusion time-scales for helium white dwarfs at this age imply that discs can prevail over time-scales of the order 10 6 years or longer, consistent with the findings of Girven et al. (2012) and the model described in Bochkarev &amp; Rafikov (2011).

The fractional luminosity of discs around polluted white dwarfs has the potential to evolve with cooling age. Farihi (2011) first noted that the highest L IR /L ∗ occur only for relatively cool, dusty white dwarfs. This trend is consistent with an increasing area available to emitting solids interior to the fixed stellar Roche radius with decreasing stellar luminosity and temperature; the radius at which grains will rapidly sublimate moves inwards as a white dwarf ages (Rocchetto et al., in preparation). The relatively faint excesses around both G166-58 and G245-58 are remarkable in two ways. First, two detections among the coolest 42 stars (the three lowest T eff bins in Fig. 6) clearly indicate a low frequency of infrared excess, especially considering these stars have nearly maximal areas available to emitting dust. Second, these two excesses are among

the most subtle L IR /L ∗ detected, and indicate that typically only tenuous rings exist at these ages.

## 5.3 Subtle discs/narrow rings are common

Prior studies suggested that narrow rings may be present around metal polluted stars without observed infrared excess (Farihi et al. 2010b), and the three new discs discovered here corroborate this view. In fact, only Spitzer observations were capable of confirming dust in these particular systems; strong infrared excesses would have been previously detected in the K -band from the ground (Kilic et al. 2006), or with WISE (Hoard et al. 2013). Owing to the degeneracy between the spatial extent and inclination for faint discs such as these, we provided two representative disc models: one narrow at moderate inclination and the other more extended but at high inclination. However, it is geometrically unlikely that all three of these discs have high inclinations, and therefore the subtle excesses imply that the disc material is confined to narrow rings around each star. This is also likely for previously published infrared excesses with low L IR /L ∗ .

The disc model for G245-58 essentially represents the limit of what is confidently detectable using Spitzer IRAC. While 0.5 stellar radii ( 1 R WD ≈ 8800 km) is the most narrow ring inferred for any polluted white dwarf, it is still comparable to many planetary rings (Esposito 1993); e.g. the main ring of Jupiter is 6500 km wide. Dust rings as narrow as a few tenths of an Earth radius are suggested by the collective infrared properties of polluted white dwarfs, but are undetectable with current instrumentation. Our data support the existence of a population of white dwarfs with very faint dust discs, most of which are currently beyond detection limits and implying that dust discs are more common than what is observed. Future missions like the James Webb Space Telescope ( JWST ) and the Space Infrared telescope for Cosmology and Astrophysics ( SPICA ) with infrared spectroscopic capabilities should detect many more faint discs.

## 6 SUMMARYAND CONCLUSIONS

We observed 15 metal polluted white dwarfs with Spitzer IRAC to search for infrared excess emission indicative of circumstellar dust. Such excesses were found for two or possibly three white dwarfs in our target sample: KUV 15519+1730, HS 2132+0941 and G24558. Notably, all the dust emissions are subtle in L IR /L ∗ ; such subtle excesses not only imply narrow rings around these three stars, but also for faint discs discovered in previous studies. All together these data suggest an undetected population of weakly emitting discs at polluted white dwarfs, possibly resembling the planetary rings of the outer solar system.

The detection of discs at 1 Gyr white dwarfs supports previous findings that discs can persist for timescales of Myr or longer. For the infrared bright discs so far detected, the frequency is highest for warm and young white dwarfs, and clearly declines on a timescale of 500 Myr, possibly implying that the reservoir of large asteroids further out in the system is depleted on this timescale.

## ACKNOWLEDGMENTS

JF gratefully acknowledges the support of the STFC via an Ernest Rutherford Fellowship. CB and MR acknowledge support from the STFC. We are grateful to D. Koester for providing the white dwarf atmospheric models used. Balmer/Lyman lines in the models were calculated with the modified Stark broadening profiles of Tremblay &amp; Bergeron (2009), kindly made available by the authors. The authors thank the anonymous reviewer for suggestions that improved the manuscript. This work is based on observations made with the Spitzer Space Telescope, which is operated by JPL/Caltech under a contract with NASA. This work has made use of data products from the SDSS, which is managed by the Astrophysical Research Consortium for the Participating Institutions, and from 2MASS, which is a joint project of the University of Massachusetts and IPAC/Caltech funded by NASA and NSF. This research has made use of the VizieR catalogue access tool and the SIMBAD database operated at CDS, Strasbourg, France, and the SAO/NASA Astrophysics Data System.

## REFERENCES

| Bergeron, P. et al., 2011, ApJ, 737, 28 Bochkarev, K. V., Rafikov, R. R., 2011, ApJ, 741,                                                                                                 |
|-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| 36 Bonsor, A., Mustill, A. J., Wyatt, M. C., 2011, MNRAS, 414, 930 Bowler, B. P. et al., 2010, ApJ, 709, 396 Chiang, E. I., Goldreich, P., 1997, ApJ, 490, 368                            |
| Debes, J. H., Sigurdsson, S., 2002, ApJ, 572, 556                                                                                                                                         |
| Debes, J. H., Sigurdsson, S., Hansen, B., 2007, AJ, 134, 1662                                                                                                                             |
| Debes, J. H., Walsh, K. J., Stark, C., 2012, ApJ, 747, 148                                                                                                                                |
| Dufour, P. et al., 2007, ApJ, 663, 1291                                                                                                                                                   |
| Dufour, P., Kilic, M., Fontaine, G., Bergeron, P., Melis, C., Bochanski, J., 2012, ApJ, 749, 6                                                                                            |
| Dupuis, J., Fontaine, G., Wesemael, F., 1993, ApJS, 87, 345                                                                                                                               |
| Esposito, L. W., 1993, AREPS, 21, 487                                                                                                                                                     |
| Farihi, J., Zuckerman, B., Becklin, E. E., 2008, ApJ, 674, 431                                                                                                                            |
| Farihi, J., Jura, M., Zuckerman, B., 2009, ApJ, 694, 805 Farihi, J., Barstow, M. A., Redfield, S., Dufour, P., Hambly,                                                                    |
| N. C., 2010a, MNRAS, 404, 2123                                                                                                                                                            |
| Farihi, J., Jura, M., Lee, J.-E., Zuckerman, B., 2010b, ApJ, 714, 1386                                                                                                                    |
| Farihi, J., 2011, in White Dwarf Atmospheres and Circumstellar Environments, ed. D. W. Hoard (Berlin: Wiley-VHC), pp 117- 171                                                             |
| Farihi, J., G¨nsicke, a B. T., Steele, P. R., Girven, J., Burl eigh, M. R., Breedt, E., Koester, D., 2012a, MNRAS, 421, 1635                                                              |
| Farihi, J., G¨nsicke, a B. T., Wyatt, M. C., Girven, J., Pringe, l J. E., King, A. R., 2012b, MNRAS, 424, 464 Fazio, G. G. et al., 2004a, ApJS, 154, 10                                   |
| Friedrich, S., Koester, D., Christlieb, N., Reimers, D., L., 2000, A&A, 363, 1040                                                                                                         |
| Fazio, G. G. et al., 2004b, ApJS, 154, 39                                                                                                                                                 |
| Wisotzki,                                                                                                                                                                                 |
| Frink, S., Mitchell, D. S., Quirrenbach, A., Fischer, D. A., Marcy, G. W., Butler, R. P., 2002, ApJ, 576, 478                                                                             |
| G¨nsicke, a B. T., Marsh, T. R., Southworth, J., Rebassa- Mansergas, A., 2006, Science, 314, 1908 G¨nsicke, a B. T., Marsh, T. R., Southworth, J., 2007, MNRAS,                           |
| 380, L35 a B. T., Koester, D., Marsh, T. R., Rebassa-Manser gas, A., Southworth, J., 2008, MNRAS, 391, L103                                                                               |
| G¨nsicke,                                                                                                                                                                                 |
| Giammichele, N., Bergeron, P., Dufour, P., 2012, ApJS, 199, 29 Girven, J., Brinkworth, C. S., Farihi, J., G¨nsicke, a B. T., Hoard, D. W., Marsh, T. R., Koester, D., 2012, ApJ, 749, 154 |
| Hoard, D. W., Debes, J. H., Wachter, S., Leisawitz, D. T.,                                                                                                                                |
| M., 2013, ApJ, 770, 21                                                                                                                                                                    |
| Cohen, Holberg, J. B., Bergeron, P., 2006, AJ, 132, 1221                                                                                                                                  |

| Hora, J. L., et al., 2008, PASP, 120, 1233 Hunt, L. K., Mannucci, F., Testi, L., Migliorini, S., Stanga, R. M., Baffa, C., Lisi, F., Vanzi, L., 1998, AJ, 115, 2594 J. A., et al., 2007, ApJ, 665, 785 M., 2003, ApJL, 584, L91                                                                                                |
|--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Johnson, Jura, Jura, M., 2008, AJ, 135, 1785                                                                                                                                                                                                                                                                                   |
| Jura, M., Farihi, J., Zuckerman, B., 2007b, ApJ, 663, 1285 Jura, M., Farihi, J., Zuckerman, B., Becklin, E. E., 2007a, AJ, 1927                                                                                                                                                                                                |
| 133, Jura, M., Farihi, J., Zuckerman, B., 2009, AJ, 137, 3191                                                                                                                                                                                                                                                                  |
| Kawka, A., Vennes, S., 2012, A&A, 538, 13 Kenyon, S. J., Shipman, H. L., Sion, E. M., Aannestad, P. A., 1988, ApJL, 328, L65                                                                                                                                                                                                   |
| Kilic, M., von Hippel, T., Leggett, S. K., Winget, D. E., 2006, ApJ,                                                                                                                                                                                                                                                           |
| 646, 474 Koester, D., 2009, A&A, 498, 517 Koester, D., 2010, in Prada Moroni, P. G., Bono, G., School of As- trophysics 'F. Lucchin' X Cycle, First Course 2008, Vol. 81 of                                                                                                                                                    |
| A., et al., in Fowler A. M. ed., Infrared Astronomical Instrumentation, Vol. 3354 of SPIE Conference Series, LIRIS: a long-slit intermediate-resolution infrared spectrograph for the WHT. pp 448-455 Martin, D. C., et al., 2005, ApJL, 619, L1 McCook, G. P., Sion, E. M., 1999, ApJS, 121,1                                 |
| 272, 660                                                                                                                                                                                                                                                                                                                       |
| Manchado,                                                                                                                                                                                                                                                                                                                      |
| Melis, C., Farihi, J., Dufour, P., Zuckerman, B., Burgasser, A. J.,                                                                                                                                                                                                                                                            |
| 91, 9 Mullally, F., Kilic, M., Reach, W. T., Kuchner, M. J., von Hippel, T., Burrows, A., Winget, D. E., 2007, ApJS, 171, 206 R. R., Garmilla, J. A., 2012, ApJ, 760, 123 W. T., Kuchner, M. J., von Hippel, T., Burrows, A., lally, F., Kilic, M., Winget, D. E., 2005a, ApJL, 635, L161 W. T., et al., 2005b, PASP, 117, 978 |
| Rafikov, Reach, W. T., Lisse, C., von Hippel, T., Mullally, F., 2009,                                                                                                                                                                                                                                                          |
| Reach, Reach, 693, 697 B., et al., 2008, PASJ, 60, 1317 S., Richer, H. B., Hansen, B. M., Stairs, I.                                                                                                                                                                                                                           |
| Sato, Sigurdsson, Thorsett, S. E., 2003, Science, 301, 193 Sion, E. M., Holberg, J. B., Oswalt, T. D., McCook, G.                                                                                                                                                                                                              |
| Wasatonic, R., 2009, AJ, 138, 1681                                                                                                                                                                                                                                                                                             |
| Skrutskie, M. F., et al., 2006, AJ, 131, 1163 Stoughton, C., et al., 2002, AJ, 123, 485 Subasavage, J. P., Henry, T. J., Bergeron, P., Dufour, P., N. C., Beaulieu, T. D., 2007, AJ, 134, 252                                                                                                                                  |
| Subasavage, J. P., Henry, T. J., Bergeron, P., Dufour, P.,                                                                                                                                                                                                                                                                     |
| N. C., 2008, AJ, 136, 899 Tremblay, P.-E., Bergeron, P., 2009, ApJ, 696, 1755 Tremblay, P.-E., Bergeron, P., Gianninas, A., 2011, ApJ, 730,                                                                                                                                                                                    |
| 128                                                                                                                                                                                                                                                                                                                            |
| Mul-                                                                                                                                                                                                                                                                                                                           |
| ApJ, H., P.,                                                                                                                                                                                                                                                                                                                   |
| Hambly, Hambly,                                                                                                                                                                                                                                                                                                                |

| Wegner, G., Swanson, S. R., 1990, AJ, 100, 1274                                     |
|-------------------------------------------------------------------------------------|
| Werner, M. W., et al., 2004, ApJS, 154, 1                                           |
| Wolszczan, A., Frail, D. A., 1992, Nature, 355, 145                                 |
| Wright, E. L., et al., 2010, AJ, 140, 1868                                          |
| Xu, S., Jura, M., 2012, ApJ, 745, 88                                                |
| Xu, S., Jura, M., Koester, D., Klein, B., Zuckerman, B., 2014, ApJ, 783, 79         |
| Zuckerman, B., Becklin, E. E., 1987, Nature, 330, 138                               |
| Zuckerman, B., Koester, D., Reid, I. N., H¨nsch, u M., 2003, ApJ, 596, 477          |
| Zuckerman, B., Koester, D., Melis, C., Hansen, B. M., Jura, M., 2007, ApJ, 671, 872 |
| Zuckerman, B., Melis, C., Klein, B., Koester, D., Jura, M., 2010, ApJ, 722, 725     |
Carolina Bergfors Jay Farihi Patrick Dufour Marco Rocchetto
Signs of a faint disc population at polluted white dwarfs
Observations of atmospheric metals and dust discs around white dwarfs provide
important clues to the fate of terrestrial planetary systems around
intermediate mass stars. We present Spitzer IRAC observations of 15 metal
polluted white dwarfs to investigate the occurrence and physical properties of
circumstellar dust created by the disruption of planetary bodies. We find
subtle infrared excess emission consistent with warm dust around KUV 15519+1730
and HS 2132+0941, and weaker excess around the DZ white dwarf G245-58, which,
if real, makes it the coolest white dwarf known to exhibit a 3.6 micron excess
and the first DZ star with a bright disc. All together our data corroborate a
picture where 1) discs at metal-enriched white dwarfs are commonplace and most
escape detection in the infrared (possibly as narrow rings), 2) the discs are
long lived, having lifetimes on the order of 10^6 yr or longer, and 3) the
frequency of bright, infrared detectable discs decreases with age, on a
timescale of roughly 500 Myr, suggesting large planetesimal disruptions decline
on this same timescale.