CFBDS J005910.90-011401.3: reaching the T-Y Brown Dwarf transition? Based in part on observations obtained with MegaPrime/MegaCam, a joint project of CFHT and CEA/DAPNIA, at the Canada-France-Hawaii Telescope (CFHT) which is operated by the National Research Council (NRC) of Canada, the Institut National des Sciences de l’Univers of the Centre National de la Recherche Scientifique (CNRS) of France, and the University of Hawaii. This work is based in part on data products produced at TERAPIX and the Canadian Astronomy Data Centre as part of the Canada-France-Hawaii Telescope Legacy Survey, a collaborative project of NRC and CNRS. Also based on observations obtained at the Gemini Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the Gemini partnership: the National Science Foundation (United States), the Science and Technology Facilities Council (United Kingdom), the National Research Council (Canada), CONICYT (Chile), the Australian Research Council (Australia), CNPq (Brazil) and SECYT (Argentina) under programs GN-2007A-Q-201 and GN-2007B-Q-3. Also based on observations made with ESO Telescopes at the La Silla Observatory under programmes 078.C-0629 and 078.A-0651. Some of the data presented herein were obtained at the W.M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California, and the National Aeronautics and Space Administration. The Observatory was made possible by the generous financial support of the W.M. Keck Foundation.

CFBDS J005910.90-011401.3: reaching the T-Y Brown Dwarf transition? thanks: Based in part on observations obtained with MegaPrime/MegaCam, a joint project of CFHT and CEA/DAPNIA, at the Canada-France-Hawaii Telescope (CFHT) which is operated by the National Research Council (NRC) of Canada, the Institut National des Sciences de l’Univers of the Centre National de la Recherche Scientifique (CNRS) of France, and the University of Hawaii. This work is based in part on data products produced at TERAPIX and the Canadian Astronomy Data Centre as part of the Canada-France-Hawaii Telescope Legacy Survey, a collaborative project of NRC and CNRS. Also based on observations obtained at the Gemini Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the Gemini partnership: the National Science Foundation (United States), the Science and Technology Facilities Council (United Kingdom), the National Research Council (Canada), CONICYT (Chile), the Australian Research Council (Australia), CNPq (Brazil) and SECYT (Argentina) under programs GN-2007A-Q-201 and GN-2007B-Q-3. Also based on observations made with ESO Telescopes at the La Silla Observatory under programmes 078.C-0629 and 078.A-0651. Some of the data presented herein were obtained at the W.M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California, and the National Aeronautics and Space Administration. The Observatory was made possible by the generous financial support of the W.M. Keck Foundation.

P. Delorme Laboratoire d’Astrophysique de Grenoble, Observatoire de Grenoble, Université Joseph Fourier, CNRS, UMR 571 Grenoble, France    X. Delfosse Laboratoire d’Astrophysique de Grenoble, Observatoire de Grenoble, Université Joseph Fourier, CNRS, UMR 571 Grenoble, France    L. Albert Canada-France-Hawaii Telescope Corporation, 65-1238 Mamaloha Highway, Kamuela, HI 96743, U.S.A.    E. Artigau Gemini Observatory, Southern Operations Center, Association of Universities for Research in Astronomy, Inc., Casilla 603, La Serena, Chile.    T. Forveille Laboratoire d’Astrophysique de Grenoble, Observatoire de Grenoble, Université Joseph Fourier, CNRS, UMR 571 Grenoble, France Institute for Astronomy, 2680 Woodlawn Drive, Honolulu, HI 96822-1839, USA.    C. Reylé Observatoire de Besançon, Institut UTINAM, University of Franche-Comté, CNRS-UMR 6213, BP 1615, 25010 Besançon Cedex, France    F. Allard Centre de Recherche Astrophysique de Lyon, UMR 5574: CNRS, Université de Lyon, Ecole Normale Supérieure de Lyon, 46 allée d’Italie, 69364 Lyon Cedex 07, France    D. Homeier Institut für Astrophysik, Georg-August-Universität, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany    A.C. Robin Observatoire de Besançon, Institut UTINAM, University of Franche-Comté, CNRS-UMR 6213, BP 1615, 25010 Besançon Cedex, France    C.J. Willott Physics Department, University of Ottawa, 150 Louis Pasteur, MacDonald Hall, Ottawa, ON K1N 6N5, Canada    Michael C. Liu Institute for Astronomy, 2680 Woodlawn Drive, Honolulu, HI 96822-1839, USA. Alfred P. Sloan Research Fellow    T.J. Dupuy Institute for Astronomy, 2680 Woodlawn Drive, Honolulu, HI 96822-1839, USA. Philippe.Delorme@obs.ujf-grenoble.fr
Key Words.:
Low mass stars: brown dwarfs – spectroscopy – Surveys
offprints: P. Delorme,

We report the discovery of CFBDS J005910.90-011401.3 (hereafter CFBDS0059), the coolest brown dwarf identified to date. We found CFBDS0059 using and images from the Canada-France-Hawaii Telescope (CFHT), and present optical and near-infrared photometry, Keck laser guide star adaptive optics imaging, and a complete near-infrared spectrum, from 1.0 to 2.2 m. A side to side comparison of the near-infrared spectra of CFBDS0059 and ULAS J003402.77-005206.7 (hereafter ULAS0034), previously the coolest known brown dwarf, indicates that CFBDS0059 is K cooler. We estimate a temperature of K and gravity of . Evolutionary models translate these parameters into an age of 1-5 Gyr and a mass of M. We estimate a photometric distance of 13pc, which puts CFBDS0059 within easy reach of accurate parallax measurements. Its large proper motion suggests membership in the older population of the thin disk. The spectra of both CFBDS0059 and ULAS J0034 shows probable absorption by a wide ammonia band on the blue side of the -band flux peak. If, as we expect, that feature deepens further for still lower effective temperatures, its appearance will become a natural breakpoint for the transition between the T spectral class and the new Y spectral type. CFBDS0059 and ULAS J0034 would then be the first Y0 dwarfs.

1 Introduction

Observed stellar and substellar atmospheres cover a continuum of physical conditions from the hottest stars ( 100 000 K) to the coolest known brown dwarfs (previously ULAS J003402.77-005206.7 (hereafter ULAS0034), T8 warren.2007). There remains however a sizeable temperature gap between the 600-700 K ULAS0034 and the 100 K giant planets of the Solar System. Many of the currently known extrasolar planets populate this temperature interval, characterized by complex atmospheric physics: matter and radiation in these cold, dense, and turbulent atmospheres couple into a very dynamical mix, where molecules and dust form and dissipate. Current atmosphere models are rather uncertain in this unexplored temperature range and they will significantly benefit from observational constraints. Two major physical transitions are expected to occur between  700 K and  400 K and strongly alter the emergent near-infrared spectra (burrows.03): NH becomes an abundant atmospheric constituent and its near-infrared bands become major spectral features, and water clouds form and deplete HO from the gas phase. The corresponding near-infrared spectral changes are likely to be sufficiently drastic that the creation of a new spectral type will be warranted (kirkpatrick.2000). kirkpatrick.1999; kirkpatrick.2000 reserved the “Y” letter for the name of that putative new spectral type.

To help fill this temperature gap, we conduct the Canada France Brown Dwarf Survey (CFBDS, (delorme.2008)), which uses MegaCam (Boulade.2003proc) and images to select very cool brown dwarfs (and high redshift quasars) on their very red colour. We present here our coolest brown dwarf discovery to date, CFBDS J005910.90-011401.3 (hereafter CFBDS0059), a T8 dwarf with evidence for near-infrared NH absorption. Section 2 describes its discovery, and presents our follow-up observations: , , , , and photometry and astrometry of CFBDS0059 and (as a reference) ULAS0034, laser guide star adaptive optics imaging and a near infrared spectrum of the new brown dwarf. Section 3 discusses the kinematics and the dynamical population membership of CFBDS0059. Section 4 compares the spectrum of CFBDS0059 with those of Gl 570D (T7.5), 2MASS J0415-09 (T8) and ULAS0034, and in the light of synthetic spectra uses that comparison to determine its effective temperature, gravity and metallicity. We also examine in that section the new spectral features which appear below 700K, in particular an NH band, and discuss new spectral indices for spectral classification beyond T8. Finally, section 5 summarizes our findings and sketches our near-future plans.

2 Cfbds

Field ultracold brown dwarfs are intrinsically very faint, and as a result they can only be identified in sensitive wide-field imaging surveys. They are most easily detected in the near infrared, and one could thus naively expect them to be most easily identified in that wavelength range. Brown dwarf spectra however very much differ from a black body, and their considerable structure from deep absorption lines and bands produces broad-band pure near-infrared colours that loop back to the blue. At modest signal to noise ratios, those colours are not very distinctive. Brown dwarfs are therefore more easily recognized by including at least one photometric band under 1 m. At those shorter wavelengths their spectra have extremely steep spectral slopes, and the resulting very red and colours easily stand out.

As discussed in detail in delorme.2008, the CFBDS survey brown dwarf identification is a two-step process:

  • we first select brown dwarfs candidates on their red colour in MegaCam images which cover several hundred square degrees.

  • -band pointed observations of these candidates then discriminate actual brown dwarfs from artefacts, and astrophysical contaminants.

The selection takes advantage of the wide field of the MegaCam camera 111http://www-dapnia.cea.fr/Sap/Phys/Sap/Activites/Projets/Megacam on the 3.6m CFHT telescope, and of the trove of deep observational material obtained with that instrument. We rely on existing images from the Very Wide component of the Canada France Hawaii Telescope Legacy Survey (CFHTLS-VW) and, for different fields, on existing images from the Red sequence Cluster Survey 2 (RCS2). We then match those with either new or new exposures to obtain and pairs. The two parent surveys also obtain and images, which for the RCS2 survey are contemporaneous with the exposure. We don’t use those as primary selection tools, but the contemporaneous exposures from the RCS2 survey provide a welcome check that an apparently red source was not, instead, a variable brighter at all wavelenghts at the epochs. Since all our fields have low galactic extinction, the only other astrophysical point sources with a similarly red are quasars at redshifts of z, which represent the other motivation of our program (e.g. willott.2007). We discriminate between quasars and brown dwarfs with -band photometry obtained on several 2 to 4m-class telescopes (delorme.2008, for details). These targeted follow-up observations also reject a number of remaining unflagged artefacts and they provide a refined spectral type estimate, thanks to the much higher signal-to-noise ratio which we typically achieve on than on .

3 Observations

We first identified CFBDS0059 as a brown dwarf candidate when comparing a 360 s RCS2 exposure from 2005 December 27 with a 500 seconds CFBDS exposure from 2006 August 31. CFBDS0059 is undetected in the  image to (), in spite of a strong detection (). The RCS2 survey obtains contemporaneous , and images, and we checked the and exposures for a counterpart. These images, which were obtained within 50 minutes of the observation, show no object at the position of CFBDS0059. This essentially excludes that the detection was due to a variable or slowly moving object with neutral colours. The () lower limit was thus secure, and made CFBDS0059 a very strong candidate for follow-up.

Figure 1: -band finder chart for CFBDS0059. The field of view is 2 by 2 arcminutes.

3.1 -band photometry and near-infrared spectroscopy

Our initial -band imaging consists of five 20-second dithered exposures with the SOFI near-infrared camera on the ESO NTT telescope at La Silla on 2006 November 12. As discussed below, the photometric system of that instrument is non-standard. We used a modified version of the jitter utility within the ESO Eclipse package (devillard.97) to subtract the background and coadd the five exposures. We extracted photometry from the resulting image using PSF fitting within Source Extractor (Bertin.1996 and Bertin & Delorme, in preparation) and obtain =. The resulting colour confirmed CFBDS0059 as a strong very late-T dwarf candidate and we triggered -band spectroscopic observations with NIRI (hodapp.03) at Gemini-North. Those were obtained in queue mode on 2007 July 30 and immediately confirmed the very cool nature of CFBDS0059. We then requested and -band spectroscopy with the same instrument, which was obtained on 2007 September 1.

All spectra were obtained through a -wide slit, which produces a resolving power of . The object is dithered along the slit. The -band spectrum is the sum of 16 second integrations, while the and band spectra each are the sum of 9 second integrations. Consecutive image pairs are pair subtracted, flat fielded using a median combined spectral flat and corrected for both spectral and spatial distortions. Spectra are extracted using a positive and a negative extraction box matched to the trace profile. A first wavelength calibration was obtained with argon lamp arc spectra taken at the end of the sequence, and the wavelength scale was then fine tuned to match the atmospheric OH-lines. Individual spectra extracted from image pairs were then normalized and median combined into final spectra. Per-pixel S/N of 25, 40 and 5 where achieved on the J, H and K-band peaks respectively. For all 3 wavelength settings the A0 star HIP10512 was observed immediately after the science observations to calibrate the instrumental spectral response and the telluric transmission.

3.2 Additional near-infrared and optical Photometry

Target Spectral band
CFBDS0059 -
ULAS0034 -
warren.2007 -
CFBDS0059 -
ULAS0034 -
warren.2007 -
Table 1: Wide band and CFHT photometry measured with Megacam, MKO-system , , , and measured with WIRCam and associated colors. For comparison we include the warren.2007, , , , and photometry of ULAS0034.

The filter of the SOFI camera on the NTT has a quite non-standard bandpass, for which the large colour corrections that result from the strong structrure in T dwarf spectra (e.g. stephens.04) have not been fully characterized. To compare the spectral energy distribution (SEDs) of CFBDS0059 and ULAS0034 (warren.2007) we therefore prefered to obtain additional near-infrared wide-band photometry with WIRCam (puget.2004) on CFHT, which uses standard Mauna Kea Observatory infrared filters (simons.02; tokunaga.02; tokunaga.05, MKO system). The observations (2007 August 1st and 5th) used dithering patterns of  arcsec amplitude for total (respectively individual) exposure times of 300 (60), 150 (30), 300 (15) and 720 (20) seconds for the , , and bands. The skies were photometric and the seeing varied between 0.8 and 1.0”.

Table 1 summarizes the magnitudes of CFBDS0059 and ULAS0034 in all available bands. The WIRCam photometry of ULAS0034 agrees with the warren.2007 measurements within better than 1 for band and within 1.5 for . The warren.2007 K-band measurement used a filter, while our WIRCam measurement uses the narrower and bluer filter. The 0.15 mag difference between these two observations is approximately consistent with the stephens.04 prediction for the effect of these different filter bandpasses at late-T spectral types. Similarly, the better short-wavelength quantum efficiency of the WIRCam detector can qualitatively explain the 0.2 magnitude (2) discrepancy between our photometry and the warren.2007 WFCam measurement. The near-IR colours of the two brown dwarfs are similar, except which is 0.5 magnitude bluer for CFBDS0059 than for ULAS0034. We will interpret the implications of this low flux when we examine the near-infrared spectrum.

3.3 Astrometry

CFBDS0059 and ULAS0034 are serendipitously just 6.3 degrees apart on the sky, and at similar photometric distances from Earth since they have similar spectral types and apparent magnitudes , and we initially entertained the idea that they might, perhaps, be part of a common moving group. The proper motion of CFBDS0059 however, measured between our 2005 Megacam and 2007 WIRCam images, and uncorrected for its parallactic motion, is ”/yr, ”/yr (Table 2). ULAS0034 on the other hand moves by ”/yr and ”/yr (warren.2007). The two proper motions are thus sufficiently different that the two brown dwarfs are clearly unrelated. We checked for main sequence common proper motion companions to CFBDS0059, which would have provided welcome age and metallicity constraints (e.g. scholz.2003) but did not find any match within a 10 arcminute radius.

RA(J2000) DEC(J2000)
hh:mm:ss dd:mm:ss
00:59:10.903 -01:14:01.13
Table 2: Astrometry of CFBDS0059. (Epoch:August 5th, 2007)

3.4 Keck Laser Guide Star Adaptive Optics Imaging

To search for binarity, we imaged CFBDS0059 on 16 January 2008 UT using the laser guide star adaptive optics (LGS AO) system (vandam.2006; wizinowich.2006) of the 10-meter Keck II Telescope on Mauna Kea, Hawaii. Conditions were photometric with better than average seeing. We used the facility IR camera NIRC2 with its narrow field-of-view camera, which produces an image scale of  mas/pixel (pravdo.2006) and a field of view. The LGS provided the wavefront reference source for AO correction, with the exception of tip-tilt motion. The LGS brightness was equivalent to a  mag star, as measured by the flux incident on the AO wavefront sensor. Tip-tilt aberrations and quasi-static changes in the image of the LGS as seen by the wavefront sensor were measured contemporaneously with a second, lower-bandwidth wavefront sensor monitoring the  mag field star USNO-B1.0 0887-0010532 (monet.2003), located 32″ away from CFBDS0059. The sodium laser beam was pointed at the center of the NIRC2 field-of-view for all observations.

We obtained a series of dithered images, offsetting the telescope by a few arcseconds, with a total integration time of 1080s. We used the filter, which has a central wavelength of 1.592 m and a width of 0.126 m. This filter is positioned near the -band flux peak emitted by late-T dwarfs. The images were reduced in a standard fashion. We constructed flat fields from the differences of images of the telescope dome interior with and without continuum lamp illumination. Then we created a master sky frame from the median average of the bias-subtracted, flat-fielded images and subtracted it from the individual images. Images were registered and stacked to form a final mosaic, with a full-width at half-maximum of 0.09″ and a Strehl ratio of 0.05. No companions were clearly detected in a region centered on CFBDS0059.

We determined upper limits from the direct imaging by first smoothing the final mosaic with an analytical representation of the PSF’s radial profile, modeled as the sum of multiple gaussians. We then measured the standard deviation in concentric annuli centered on the science target, normalized by the peak flux of the targets, and adopted 10 as the flux ratio limits for any companions. These limits were verified with implantation of fake companions into the image using translated and scaled versions of the science target.

Figure 2: Adaptive Optic detections limits as a function of the separation from CFBDS0059. Dashed lines show the equivalent mass and temperature for a given contrast. The first value assumes an age of 1 Gyr and the second an age of 5 Gyr.

Figure 2 presents the final upper limits on any companions. We employed the ”COND” models of baraffe.03 to convert the limits into companion masses, for assumed ages of 1 and 5 Gyr and a photometric distance estimate of 13 pc. We assumed any cooler companions would have similar colors to CFBDS0059.

4 Kinematics

Figure 3: Thin disk probability membership contours in proper motion space from the Besancon stellar population model. The contours are generated for synthetic stars with distances between 10 and 20 pc, belonging to the thin disk (small dots) and the thick disk (small stars, with the density of the latter increased by a factor of 10 for display purposes). Less than one halo star would appear on the plot. Based on their measured proper motion, the likelihood that CFBDS0059 (large open square) and ULAS0034 (large open triangle) belongs to the thin disk are 95% and 99%.

We estimate a spectrophotometric distance for CFBDS0059 by adopting M, based on an approximate T9/Y0 spectral type (discussed below) and on an extrapolation of the M versus spectral type relation of knapp.04 beyond T8 (2MASSJ, hereafter 2M0415). This extrapolation is consistent with the chabrier.00 models which predict between brown dwarfs of 750K and 625K (like 2M0415 and CFBDS0059). The resulting d =  pc has significant systematic uncertainties, because spectral typing beyond T8 is just being defined, and especially because the linear 1 magnitude/subtype decline seen at earlier subtypes may not continue beyond T8. The adaptive optics observations exclude any companion of similar luminosity beyond 1.2 AU, but CFBDS0059 could still of course be a tighter binary. Its small distance fortunately puts CFBDS0059 within easy reach of modern parallax measurements.

We use the Besancon stellar population model (robin.2003) to generate synthetic stars between 10 pc and 20 pc for the thin (dots) and thick (star symbol) disk populations at the galactic position of CFBDS0059. Fig 3 shows their proper motions together with those of CFBDS0059 and ULAS0034. The contour lines show the probability that an object with a given proper motion belongs to the thin disk rather than the thick disk (halo membership probabilities are negligible). CFBDS0059, at its probable distance, is well within the 95% probability thin disk membership region, and ULAS0034 is within the 99% probability region. In spite of its somewhat high proper motion for an object beyond 10 pc, CFBDS0059 therefore most likely belongs to the thin disk. The mean age of the simulated stars in the region of the proper motion diagram occupied by CFBDS0059 is 4 Gyr, suggesting that it is an older member of the thin disk. That age is consistent with the 1 to 5 Gyr range derived below from comparison to COND evolutionary models (baraffe.03). As any kinematic age for an individual star, this determination has large error bars, but it suggests that CFBDS0059 might be older than ULAS0034.

5 Spectral comparison and atmospheric parameters

Fig 4 and Fig 5 present our spectrum of CFBDS0059, together with those of ULAS0034 (warren.2007, T8), 2M0415 (burgasser.03, T8) and Gl 570D (burgasser.2000, T7.5), which successively were the coolest known brown dwarfs. Thanks to their earlier discovery, 2M0415 and Gl 570 have the best characterized atmospheric parameters (saumon.06; saumon.07), and they provide the most solid baseline for a differential study. The warren.2007 spectrum of ULAS0034, kindly communicated by S. Leggett, was obtained with GNIRS on Gemini South and its  = 500 resolution matches that of our NIRI spectrum of CFBDS0059. We downloaded the burgasser.03; burgasser.02 OSIRIS spectra of the two other brown dwarfs from the Ultracool Dwarf Catalog222http://www.iac.es/galeria/ege/catalogo_espectral/index.html, and degraded their original spectral resolution of to match that of the GNIRS and NIRI spectra. Stronger telluric absorption from the lower altitude telescopes explains the wider blanked regions in the corresponding spectra, but doesn’t measurably affect any comparison: as illustrated by CFBDS0059, late-T dwarfs have essentially negligible flux wherever telluric HO absorption matters. Because the OSIRIS spectra do not cover the band, we complement them by lower resolution spectra from geballe.01 and knapp.04 for .

Figure 4: 0.9 m - 2.3 m spectra of CFBDS0059 and the three other coolest brown dwarfs. The spectra are normalized to unit flux densities at their 1.27m.peak, and vertically offset for clarity. The main T-dwarf spectral features are labeled. The temperatures of Gl 570D and 2M0415 are from the careful spectroscopic analyses of saumon.06 and saumon.07. Those of CFBDS0059 and ULAS0034 are from our versus index (Fig. 6).
Figure 5: , and -bands spectra of the four cool brown dwarfs. The green vertical lines in the 1.15 m - 1.35 m panel mark the two potassium doublets. The grey band in the 1.45 m - 1.65 m panel indicates the approximate limits of the new absorption bands, which we discuss in the text.

5.1 Atmospheric parameters

Atmospheric parameters of ultracool dwarfs are ideally determined from a combination of near and mid-IR information (e.g. saumon.06; saumon.07), but low resolution near-infrared spectra alone provide a useful proxy when mid-IR photometry and spectra are not (yet) available (e.g. burgasser.2006; leggett.2007). burgasser.2006 used a grid of solar metallicity cool brown dwarfs to calibrate two spectral ratios, HO- and , which respectively measure the strength of HO absorption at 1.15 m and the flux ratio between the and peaks, to T and . warren.2007 however found that HO- essentially saturates below T = 750K, and therefore chose not to use this spectral index for spectral types later than T8. They demonstrate on the other hand that the combination of the index with the width of -band peak, parametrised by their index, becomes a good T and diagnostic at T900 K, and remains useful significantly below  750K. We adopt their method.

Table 3 lists our measurement of these two indices for CFBDS0059, and Fig 6 compares them with the warren.2007 measurements for Gl 570D, HD3651B, 2M0415 and ULAS0034. To derive T and from the indices we use model indices from solar-abundance BT-settl atmospheric models (warren.2007; allard.03, Allard et al. 2008 in prep). The model with NH at chemical equilibrium abundance clearly produces too much absorption in the blue side of the -band, confirming the finding of saumon.06; saumon.07 that non-equilibrium processes keep the NH partial pressure well below its equilibrium value. We then use models that keep the abundances of NH and N at a fixed value in all parts of the atmosphere where the reaction timescale exceeds the mixing timescale, which typically occurs at the 600-800K temperature level. These ”quenched” models agree much better with the observed band shape.

As a first order correction for the remaining imperfections of the theoretical spectra, the model indices are shifted into agreement with the measurements of 2M0415 at the [T=750K; =5.00 and [M/H]=0] determined for that brown dwarf by saumon.07. The T=800K and =5.35 resulting from this calibration for Gl 570D are consistent with the T=800-820K and =5.1-5.25 derived by saumon.06 from a complete spectral analysis. For HD3651B, T=820-890K and =5.1-5.3 resulting from this calibration are roughly consistent with the T=780-840K and =5.1-5.5 derived by liu.07.

CFBDS0059 and ULAS0034 have very similar indices, but the new brown dwarf has a significatively smaller index. Visual comparison of the two spectra (Fig. 4) confirms that CFBDS0059 does have a weaker -band peak than any of the 3 comparison cool brown dwarfs. As widely discussed in the recent literature (e.g. liu.07, Fig. 3 in burgasser.2006, or Fig. 3 in leggett.2007), for a fixed metallicity a weaker -band peak is evidence for either a lower temperature or a higher gravity. The index lifts this degeneracy: it indicates, again assuming identical chemical compositions for the two brown dwarfs, that CFBDS0059 is cooler by 5015 K and has a 0.150.1 higher than ULAS0034.

As also discussed by warren.2007, the above uncertainties only reflect the random errors in the spectral indices. They are appropriate when comparing two very similar objects, like CFBDS0059 and ULAS0034, since systematic errors then cancel out. They must on the other hand be increased to compute absolute effective temperatures and gravity: one then needs to account for the uncertainties on the 2M0415 parameters which anchor the Fig 6 grid (TK and ; saumon.07), and for the uncertainties in the atmospheric models which may distort the grid between its anchor point [T=750K; =5.00 and [M/H]=0] and the 600K region of interest here. We conservatively adopt T=62050K and =4.750.3.

This 2-parameter analysis obviously cannot determine all three main atmospheric parameters (T,, and metallicity). As discussed by warren.2007, it actually determines the temperature with no ambiguity but leaves a combination of and undetermined, and they demonstrated that in the versus plot metallicity is degenerate with surface gravity, with . CFBDS0059 is thus definitely cooler than ULAS0034, but from versus diagram it can have either higher surface gravity or lower metallicity. This degeneracy affects the full JHK-band spectrum, where any metallicity vs gravity difference is at most very subtle. It is however lifted by the shape of the -band peak (Figs. 3 of burgasser.2006 or leggett.2007), since lower metallicity shifts the -band flux density peak of submetallic brown dwarfs significantly blueward. Fig 5 shows no such shift, and the two objects therefore have similar metallicities.

Fig. 7 overlays the observed CFBDS0059 spectrum with the synthetic spectrum for the closest point of the solar metallicity atmospheric model grid. Except on the red side of the -band, model and observations agree well, boosting our confidence in the derived atmospheric parameters. The main remaining predictive shortcoming of the models is their overestimated absorption on the red side of the -band peak. The principal opacity source in this region is the methane band centred at , for which comprehensive theoretical predictions are available, but only for transitions from the vibrational ground state (as will be discussed in detail in Homeier et al., in preparation). To make up for the missing absorption from higher bands, which constitutes a significant fraction of the opacity at brown dwarf temperatures, a constant empirical correction factor was used. This correction must in turn lead to some overestimate of the CH absorption as we reach the lower end of the T dwarf temperature range. Another possible source of errors are uncertainties in the models’ temperature profile. The BT-Settl models self-consistently describe gravitational settling, convective turbulence, advection, condensation, coalescence and coagulation of condensates to predict the formation and vertical extent of cloud layers (allard.03; helling.08). In late T dwarfs these clouds are predicted to reside deep in the optically thick part of the atmosphere. Their opacity is thus not directly visible in the spectrum, but it may still impact the thermal structure, and thus the relative abundance especially of temperature-sensitive species like CH.

Another (less serious) disagreement between the models and the observed spectra occurs in the band. The models overestimate the flux on the blue side of the -peak, and they imperfectly reproduce the general shape of the peak. As discussed below, the opacities in that band are dominated by pressure-broadened wings of the 0.77  K I line on the blue side and CH on the red side.

Figure 6: versus indices. The grid represents indices measured on solar metallicity BT-settl model spectra, shifted into agreement of the [T=750K; =5.00 and [M/H]=0] with the measured indices of 2M0415. The error bars represent the 1 uncertainties of the measured spectral indices.

For baraffe.03 evolutionary models, the T570-670K and 4.45-5.05 determined above translate into an age of 1-5 Gyr and a mass of 15M (for 1 Gyr) to 30M (for 5 Gyr). The kinematics of CFBDS0059 suggests that it belongs to an older population, and therefore slightly favour a higher mass and older age.

Figure 7: Overlay of the CFBDS0059 spectrum with the solar metallicity [T=600-650K; =4.75] BT-settl synthetic spectrum. The two spectra are scaled to agree at their 1.27m flux peak. The “quenched NH3” models are chemical equilibrium models which enforce a constant abundance of ammonia in the cooler regions of the atmosphere.

5.2 Individual spectral features

Fig. 5 zooms on the , and -band peaks of the four cool brown dwarf spectra. The published OSIRIS spectra of Gl 570D and 2M0415 do not cover the band, which instead is plotted from the lower resolution spectra of geballe.01 and knapp.04. For easier comparison, the CFBDS0059 and ULAS0034 -band spectra are smoothed to that resolution.

Direct comparison of the four spectra can be used to shed light on incipient new features and atmospheric chemistry. Features which are seen in both CFBDS0059 and ULAS0034 are likely to be real even when their significance is modest in each object, and those which are absent or weaker in the two hotter brown dwarfs, can reasonably be assigned to low temperature molecules. Conversely, features which disapear in the two cooler objects trace higher temperature species.

As discussed above, the -band spectra of CFBDS0059 and ULAS0034 do not differ much. Given the strong sensitivity of that band to [M/H] that implies that the two objects have similar chemical compositions. The shape of the -band peaks of these two coolest brown dwarfs on the other hand differ from that of Gl 570D and 2M0415, with the CFBDS0059 and ULAS0034 peaks extending further into the blue. The dominant absorbers in the blue wings of the -band peak is the pressure-broadened wing of the 0.77 m K I line (e.g. burgasser.2006), which must weaken as K I depletes from the gas phase under T=700K. As anticipated by leggett.2007, the slope of the blue side of the -band peak therefore shows good potential as an effective temperature diagnostics beyond spectral type T8.

The strength of the J-band K I doublet is a good gravity estimator in ultracool dwarfs (e.g. knapp.04), because an increased pressure at a fixed temperature favors KCl over K (lodders.99) and consequently weakens atomic potassium features. At T 750-800K the the -band K I doublets remain weakly visible and useful as a gravity proxy (Fig. 7 of knapp.04). At T 700K on the other hand, the K I doublets have completely vanished at the resolution of the current spectra (Fig. 5), even at the probably lower gravity of ULAS0034. Potassium is thus mostly converted to KCl (or perhaps other compounds) in the relevant photospheric layers.

The strongest new feature is wide absorption on the blue side of the -band, at . It is conspicuous in CFBDS0059 and well detected in ULAS0034, and with hindsight is weakly visible in the 2M0415 spectrum (Fig. 5). It is however clearly stronger at T 700K. To visually emphasize this broad feature, we bin the spectra to R=100 and overlay the four H-band spectra (Fig. 8, left panel). Absorption sets in at and becomes deeper for . These wavelengths overlap with strong HO and NH bands. Either molecule could a priori be responsible for the absorption.

Figure 8: Left: -band spectrum of the four cool brown dwarfs binned to R=100. The spectra are normalized at m. The integration intervals for the numerator and denominator of our proposed NH index are marked. Right: BT-Settl synthetic spectra for [K; ] with and without near-infrared NH opacity; the NH abundance is at its chemical equilibrium value.

5.3 Near infrared ammonia signatures

Those molecules unfortunately have imperfect opacity data, and the NH laboratory line lists in particular are incomplete below m. Computed ammonia opacities are therefore strictly lower limits. leggett.2007 compare synthetic spectra computed with and without NH opacity, using the irwin.1999 line list for m, and find that ammonia absorption in cold brown dwarfs strongly depletes the blue wing of the band (their Fig.10). Similarly, Fig. 9 of saumon.2000 plots synthetic -band spectra with and without NH opacity, and find differences in two wavelength ranges: the NH-rich model is significantly more absorbed for m and it has weaker but significant absorption in the range.

Fig 8 right panel plots two BT-Settl models for [K; ], without any near-infrared NH opacity, and with NH opacity for that molecule at its chemical equilibrium abundance. As discussed above the BT-Settl models do not reproduce the observed -peak shape very well, and a quantitative comparison is thus difficult. The comparison of the two models nonetheless confirms the saumon.2000 conclusion that ammonia produces strong absorption below m and weaker residual out to 1.595m. These model predictions qualitatively match the behaviour seen in Fig 8, left panel.

To emphasize the changes in brown dwarfs spectra when their effective temperature decreases from 800 to 600 K, we plot in Fig 9 the ratio of the spectra of CFBDS0059 and Gl 570D. The signal to noise ratio of the resulting -band spectrum is too low for detailed analysis, and we therefore focus on the , and flux peaks. To avoid confusion from changes in the temperature-sensitive methane bands, we also mostly ignore the parts of the spectrum affected by CH absorption bands, hatched in dark and light grey for respectively stronger and weaker bands. Fig. 10 shows the equivalent plot for ULAS0034, which is very similar.

The -band spectrum ratio prominently shows the new absorption band, which outside the CH band closely matches the 300 K NH transmission spectrum of irwin.1999. Both spectra are strongly absorbed between 1.49 and 1.52 m and rebound from 1.52 to 1.57m. Water absorption, by contrast, is a poor match to the features of spectrum ratio. The strongest water absorption (as computed from the HITRAN molecular database for a 600 K temperature) occurs below 1.49m, at significantly bluer wavelengths than the CFBDS0059 absorption feature.

Some weaker but still significant bands of the irwin.1999 laboratory ammonia spectrum occur in the band. Those again match features of the CFBDS0059/Gl 570D flux ratio, but that agreement is much less conclusive: water and ammonia absorptions overlap on the red side -band peak, and CH absorption affects the blue side of the peak. A 1.25-1.27 m feature is seen on both flux ratios and on the ammonia transmission, and could be due to ammonia since it is clear of any strong HO absorption band. The slight wavelength shift between the laboratory and observed features however leaves that association uncertain. Detailed synthetic spectra based on fully reliable opacities will be needed to decide whether NH absorption matters in the band at the effective temperature of CFBDS0059. The main pattern in the -band is a blue slope, which reflects the weaker pressure-broadened K I wing in the cooler brown dwarf. The weak 1.03m NH band is not seen.

Figure 9: Flux ratio between CFBDS0059 and Gl 570D (black), together with the laboratory room temperature transmission spectrum of NH (irwin.1999) (red, top panel) and the 600 K HO transmission spectrum computed from the HITRAN molecular database (red, bottom panel). The grey bands mark the parts of the spectrum affected by strong (dark grey) or moderate (light grey) CH absorption.
Figure 10: Flux ratio between ULAS0034 and Gl 570D (black). The overlays repeat those of Fig 9.

Ammonia is easily detected in mid-infrared SPITZER spectra for all spectral types cooler than T2 (Roellig.04; Cushing.06), though significantly weaker than initially expected because mixing from lower atmospheric levels reduces its abundance in the high atmosphere below the local equilibrium value (saumon.06). Weak near-infrared absorption by ammonia has been tentatively detected by saumon.2000 in the T7p dwarfs Gl 229B, but CFBDS0059 and ULAS0034 provide the first incontrovertible evidence for a strong near-infrared NH band in brown dwarf spectra.

This conclusion contrasts with warren.2007 finding possible but inconclusive evidence for ammonia in ULAS0034. The main difference between the two analyses is that warren.2007 focused on a higher resolution search, at a necessarily lower signal to noise ratio, for individual NH lines between 1.5 and 1.58m. We instead looked for the global signature of the absorption band, which only becomes obvious when looking at the full -band spectrum.

5.4 Spectral type

Table 3 lists for CFBDS0059 the spectral indices used by the spectral classification scheme of burgasser.06b, which refines the previous schemes of geballe.02 and burgasser.02. These indices would imply a T8 classification, identical to that of 2M0415. As discussed above however, the near-infrared spectrum of CFBDS0059 demonstrates that it is over 100 K cooler than 2M0415 and shows clearly different spectral features. Based on the new indices we present later, CFBDS0059 should be assigned a later spectral type. The almost identical burgasser.06b indices of the two brown dwarfs instead reflect those indices measuring HO and CH absorption bands which saturate and lose their effective temperature sensitivity at the T8 spectral type of 2M0415. Beyond T8 the burgasser.06b classification scheme therefore needs to be extended, with new spectral indicators that do not saturate until significantly later spectral types.

Spectral numerator denominator value spectral value spectral
indice wavelength wavelength type type
m CFBDS0059 ULAS0034
HO- 1.14-1.165 1.26-1.285 0.0290.005 T8 0.0120.007 T8
CH- 1.315-1.34 1.26-1.285 0.1650.005 T8 0.0140.009 T8
HO- 1.48-1.52 1.56-1.60 0.1190.008 T8 0.1330.010 T8
CH- 1.635-1.675 1.56-1.60 0.0840.002 T8 0.0960.006 T8
CH- 2.215-2.255 2.08-2.12 0.1280.037 T7 0.0910.015 T7
2.06-2.10 1.25-1.29 0.1010.002 0.1280.003
W 1.18-1.23 1.26-1.285 0.2570.004 0.2620.006
Table 3: Measured spectral classification indices for CFBDS0059 and ULAS 0034. The first six indices form the base of the burgasser.06b classification scheme, and the table includes the corresponding spectral classification on that scale. The index is a recent addition proposed by warren.2007.

Fully defining this extension is beyond the scope of the present paper, since two known objects beyond T8 are not enough to explore spectral variability, but one can nonetheless start exploring. Since the main new feature is NH absorption in the blue wings of the -band peak, we define a new NH- index as

(1)

Its numerator and denominator range are plotted in Fig 8. The numerator integrates the flux within the main NH band and its denominator measures the bulk of the -band peak (we note that in cooler objects, the denominator could be affected by some NH absorption; its integration boundaries might thus need to be refined after such objects have been discovered). We compute this index for Gl 570D, HD3651B, 2M0415, ULAS0034 and CFBDS0059 (Table 4), and find that it strongly decreases from Gl570D to ULAS0034 and CFBDS0059 (which have very similar NH-).

Gl570D HD3651B 2M0415 ULAS0034 CFBDS0059
NH- 0.6720.008 0.660.005 0.6250.003 0.5160.008 0.5260.005
Spec. Type T7.5 T7.5 T8 T9 or Y0 T9 or Y0
Table 4: NH- indices for the five brown dwarfs discussed in this paper. The spectral classification either follows burgasser.06b or is derived from the NH-H and W indices, as discussed in the text

Over the limited effective temperature range spanned by Gl570D, HD3651B, 2M0415, ULAS0034 and CFBDS0059, and as far as one can infer from just 5 examples, the NH- and W indices correlate strongly (Fig. 11). The numerator of W is centered at wavelengths where both ammonia (Fig. 10 of leggett.2007) and CH have significant opacity, and future modeling work should be able to establish whether the two indices probe the same molecule or not. Since the near-infrared spectra of ULAS0034 and CFBDS0059 differ significantly more from that of the T8 2M0415 than the latter differs from the T7.5 Gl570D (as quantitatively demonstrated by Fig. 11), it is natural to assign a full spectral subtype to the interval between the two coolest brown dwarfs and 2M0415. By that reasoning, and if ULAS0034 and CFBDS0059 are considered as T dwarfs, their spectral type should be T9, or perhaps slightly later.

Figure 11: NH-H index versus index. The error bars represent the 1 uncertainties of the measured spectral indices.

The T spectral class however is quite unlikely to remain the last spectral type, since for sufficiently low effective temperatures atmospheric models predict major changes in visible and near-infrared brown dwarf spectra: NH bands are predicted to appear in, and eventually to dominate, the near-infrared spectrum, the strong pressure-broadened optical lines of Na I and K I are predicted to disappear as those atomic species get incorporated into molecules and solids, and water clouds are predicted to form and to largely deplete water from the gas phase (burrows.03; kirkpatrick.05). Since spectral classification, for mostly practical reasons, is traditionally based on optical and near-infrared spectra, such a major transition will justify the introduction of a new spectral type, for which the Y letter has long been reserved (kirkpatrick.1999; kirkpatrick.2000). If the m NH band keeps deepening as the effective temperature decreases further, and eventually becomes a major spectral feature, its appearance at 650 K will become a natural transition between the T and Y spectral classes. ULAS0034 and CFBDS0059 would then be the first Y dwarfs, and the prototypes for Y0 brown dwarfs, rather than T9. That decision will to some extent remain a matter of convention, but it must in any case wait until larger numbers of similarly cool brown dwarfs can document spectral trends in finer detail, and preferably over a wider effective temperature range.

6 Summary and conclusions

We have reported the discovery of CFBDS0059, a very cool brown dwarf, discovered in the CFBDS survey (delorme.2008). Its effective temperature is 5015 K cooler than that of ULAS0034, most likely making it the coolest brown dwarf known at the present time. High spatial resolution imaging establishes than CFBDS0059 has no similarly bright companion beyond 0.09”, and no companion with a contrast under 3.5 magnitude beyond 0.3” (respectively 1.2 and 3.9 AU at the 13 pc photometric distance). Its kinematics suggest, with significant error bars, a 4 Gyr age at which CFBDS0059 would be a 30M brown dwarf. The atmospheric parameters of CFBDS0059 however are compatible with any age from 5 Gyr down to 1 Gyr, for which its mass would be 15M. A trigonometric parallax measurement together with mid-infrared photometry and spectroscopy with SPITZER will significantly refine its physical parameters, as demonstrated by saumon.07 for slightly warmer brown dwarfs.

We assign absorption in the blue wing of the -band peaks of both ULAS 0034 and CFBDS 0059 to an NH band. If that assignment is confirmed, and if, as we expect, the band deepens at still lower effective temperatures, its development would naturally define the scale of the proposed Y spectral class. ULAS 0034 and CFBDS 0059 would then become the prototypes of the Y0 sub-class.

The CFBDS survey has to date identified two brown dwarfs later than T8, CFBDS0059 and ULAS0034 (which we identified independently of warren.2007, delorme.2008) in the analysis of approximately 40% of its final 1000 square degree coverage. We therefore expect to find another few similarly cool objects, and hopefully one significantly cooler one.

CFBDS0059 and ULAS0034 provide a peek into the atmospheric physics for conditions that start approching those in giant planets, and the future discoveries that can be expected from CFBDS, ULAS, and Pan-STARRS will further close the remaining gap. They also bring into a sharper light the remaining imperfections of the atmopsheric models, and emphasize in particular the importance of more complete opacity data. Our analysis relies on a room temperature absorption NH spectrum, but higher excitation bands than can be excited at 300 K must matter in T = 600 K brown dwarfs. The eventual identification of ammonia absorption in the band will also need complete opacity information for HO and CH and full spectral synthesis, since the bands of the three molecules overlap in that spectral range.

The spectral indices that define the T dwarf spectral class saturate below 700 K, and new ones will be needed at lower effective temperatures. We introduce one here, NH-, which measures the likely NH absorption in the band. Together with the index of warren.2007 and the slope of the blue side of the -band peak (leggett.2007), it will hopefully define a good effective temperature sequence. Metallicity and gravity diagnostics are less immediately apparent, but will need to be identified as well.

Acknowledgements.
We are grateful to our referee, Sandy Leggett, for her very detailed report and numerous suggestions which significantly improved this paper. We would like to thank the observers and queue coordinators who carried our service observations at CFHT (programs 05BC05,06AC20,07BD97) and Gemini-North (GN-2007A-Q-201, GN-2007B-Q-3). We also thanks the NTT and Keck Observatory support astronomers for their help during the observations which led to these results. We thank S. Leggett and S Warren for providing their spectrum of ULAS0034 in a convenient numerical format, and Sandy Leggett for communicating spectra of 2M0415 and Gl 570D and the room temperature absorption spectrum of NH. We would also like to thank David Ehrenreich for providing us the computed 600 K water absorption. Financial support from the ”Programme National de Physique Stellaire” (PNPS) of CNRS/INSU, France, is gratefully acknowledged. MCL acknowledges support for this work from NSF grant AST-0507833 and an Alfred P. Sloan Research Fellowship.

References

Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
   
Add comment
Cancel
Loading ...
176861
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description