skip to main content

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Friday, December 13 until 2:00 AM ET on Saturday, December 14 due to maintenance. We apologize for the inconvenience.


Title: The Hawaii Infrared Parallax Program. VI. The Fundamental Properties of 1000+ Ultracool Dwarfs and Planetary-mass Objects Using Optical to Mid-infrared Spectral Energy Distributions and Comparison to BT-Settl and ATMO 2020 Model Atmospheres
Abstract

We derive the bolometric luminosities (Lbol) of 865 field-age and 189 young ultracool dwarfs (spectral types M6–T9, including 40 new discoveries presented here) by directly integrating flux-calibrated optical to mid-infrared (MIR) spectral energy distributions (SEDs). The SEDs consist of low-resolution (R∼ 150) near-infrared (NIR; 0.8–2.5μm) spectra (including new spectra for 97 objects), optical photometry from the Pan-STARRS1 survey, and MIR photometry from the CatWISE2020 survey and Spitzer/IRAC. OurLbolcalculations benefit from recent advances in parallaxes from Gaia, Spitzer, and UKIRT, as well as new parallaxes for 19 objects from CFHT and Pan-STARRS1 presented here. Coupling ourLbolmeasurements with a new uniform age analysis for all objects, we estimate substellar masses, radii, surface gravities, and effective temperatures (Teff) using evolutionary models. We construct empirical relationships forLbolandTeffas functions of spectral type and absolute magnitude, determine bolometric corrections in optical and infrared bandpasses, and study the correlation between evolutionary model-derived surface gravities and NIR gravity classes. Our sample enables a detailed characterization ofBT-SettlandATMO2020 atmospheric model systematics as a function of spectral type and position in the NIR color–magnitude diagram. We find the greatest discrepancies between atmospheric and evolutionary model-derivedTeff(up to 800 K) and radii (up to 2.0RJup) at the M/L spectral type transition boundary. With 1054 objects, this work constitutes the largest sample to date of ultracool dwarfs with determinations of their fundamental parameters.

 
more » « less
PAR ID:
10478194
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
959
Issue:
1
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 63
Size(s):
Article No. 63
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Spectral energy distributions (SEDs) from X-ray to far-infrared (FIR) wavelengths are presented for a sample of 1246 X-ray-luminous active galactic nuclei (AGNs;L0.5–10 keV> 1043erg s−1), withzspec< 1.2, selected from Stripe 82X, COSMOS, and GOODS-N/S. The rest-frame SEDs show a wide spread (∼2.5 dex) in the relative strengths of broad continuum features at X-ray, ultraviolet (UV), mid-infrared (MIR), and FIR wavelengths. A linear correlation (log–log slope of 0.7 ± 0.04) is found betweenLMIRandLX. There is significant scatter in the relation between theLUVandLXowing to heavy obscuration; however, the most luminous and unobscured AGNs show a linear correlation (log–log slope of 0.8 ± 0.06) in the relation above this scatter. The relation betweenLFIRandLXis predominantly flat, but with decreasing dispersion atLX> 1044erg s−1. The ratio between the “galaxy-subtracted” bolometric luminosity and the intrinsicLXincreases from a factor of ∼10 to 70 from logLbol/(erg s−1) = 44.5 to 46.5. Characteristic SED shapes have been determined by grouping AGNs based on relative strengths of the UV and MIR emission. The averageL1μmis constant for the majority of these SED shapes, while AGNs with the strongest UV and MIR emission have elevatedL1μm, consistent with the AGN emission dominating their SEDs at optical and near-infrared wavelengths. A strong correlation is found between the SED shape and both theLXandLbol, such thatLbol/LX= 20.4 ± 1.8, independent of the SED shape. This is consistent with an evolutionary scenario of increasingLbolwith decreasing obscuration as the AGN blows away circumnuclear gas.

     
    more » « less
  2. Abstract

    We present a uniform forward-modeling analysis of 90 late-M and L dwarfs in nearby young (∼10–200 Myr) moving groups, the Pleiades, and the Hyades using low-resolution (R≈ 150) near-infrared (0.9–2.4μm) spectra and the BT-Settl model atmospheres. We derive the objects’ effective temperatures, surface gravities, radii, and masses by comparing our spectra to the models using a Bayesian framework with nested sampling and calculate the same parameters using evolutionary models. Assuming the evolutionary-based parameters are more robust, our spectroscopically inferred parameters from BT-Settl exhibit two types of systematic behavior for objects near the M-L spectral type boundary. Several objects are clustered aroundTeff≈ 1800 K andlogg5.5dex, implying impossibly large masses (150–1400MJup), while others are clustered aroundTeff≳ 3000 K andlogg3.0dex, implying unphysically low masses and unreasonably young ages. We find the fitted BT-Settl model spectra tend to overpredict the peakJ- andH-band flux for objects located near the M-L boundary, suggesting the dust content included in the model atmospheres is insufficient to match the observations. By adding an interstellar medium–like reddening law to the BT-Settl model spectra, we find the fits between models and observed spectra are greatly improved, with the largest reddening coefficients occurring at the M-L transition. This work delivers a systematic examination of the BT-Settl model atmospheres and constitutes the largest spectral analysis of benchmark late-M- and L-type brown dwarfs to date.

     
    more » « less
  3. Abstract

    Planetary-mass objects and brown dwarfs at the transition (Teff∼ 1300 K) from relatively red L dwarfs to bluer mid-T dwarfs show enhanced spectrophotometric variability. Multiepoch observations support atmospheric planetary-scale (Kelvin or Rossby) waves as the primary source of this variability; however, large spots associated with the precipitation of silicate and metal clouds have also been theorized and suggested by Doppler imaging. We applied both wave and spotted models to fit near-infrared (NIR), multiband (Y/J/H/K) photometry of SIMP J013656.5+093347 (hereafter SIMP0136) collected at the Canada–France–Hawaii Telescope using the Wide-field InfraRed Camera. SIMP0136 is a planetary-mass object (12.7 ± 1.0MJ) at the L/T transition (T2 ± 0.5) known to exhibit light-curve evolution over multiple rotational periods. We measure the maximum peak-to-peak variability of 6.17% ± 0.46%, 6.45% ± 0.33%, 6.51% ± 0.42%, and 4.33% ± 0.38% in theY,J,H, andKbands, respectively, and find evidence that wave models are preferred for all four NIR bands. Furthermore, we determine that the spot size necessary to reproduce the observed variations is larger than the Rossby deformation radius and Rhines scale, which is unphysical. Through the correlation between light curves produced by the waves and associated color variability, we find evidence of planetary-scale, wave-induced cloud modulation and breakup, similar to Jupiter’s atmosphere and supported by general circulation models. We also detect a 93.°8 ± 7.°4 (12.7σ) phase shift between theHKandJHcolor time series, providing evidence for complex vertical cloud structure in SIMP0136's atmosphere.

     
    more » « less
  4. Abstract

    Constraining L dwarf properties from their spectra is challenging. Near-infrared (NIR) spectra probe a limited range of pressures, while many species condense within their photospheres. Condensation creates two complexities: gas-phase species “rain out” (decreasing in abundances by many orders of magnitude) and clouds form. We designed tests using synthetic data to determine the best approach for retrieving L dwarf spectra, isolating the challenges in the absence of cloud opacity. We conducted atmospheric retrievals on synthetic cloud-free L dwarf spectra derived from the Sonora Bobcat models at SpeX resolution using a variety of thermal and chemical abundance profile parameterizations. For objects hotter than L5 (Teff∼ 1700 K), the limited pressure layers probed in the NIR are mostly convective; parameterized pressure–temperature (PT) profiles bias results and free, unsmoothed profiles should be used. Only when many layers both above and below the radiative-convective boundary are probed can parameterized profiles provide accurate results. Furthermore, a nonuniform abundance profile for FeH is needed to accurately retrieve bulk properties of early-to-mid L dwarfs. Nonuniform prescriptions for other gases in NIR retrievals may also be warranted near the L/T transition (CH4) and early Y dwarfs (Na and K). We demonstrate the utility of using realistic, self-consistent models to benchmark retrievals and suggest how they can be used in the future.

     
    more » « less
  5. ABSTRACT

    We present a sample of 19 583 ultracool dwarf candidates brighter than z ≤23 selected from the Dark Energy Survey DR2 coadd data matched to VHS DR6, VIKING DR5, and AllWISE covering ∼ 480 deg2. The ultracool candidates were first pre-selected based on their (i–z), (z–Y), and (Y–J) colours. They were further classified using a method that compares their optical, near-infrared, and mid-infrared colours against templates of M, L, and T dwarfs. 14 099 objects are presented as new L and T candidates and the remaining objects are from the literature, including 5342 candidates from our previous work. Using this new and deeper sample of ultracool dwarf candidates we also present: 20 new candidate members to nearby young moving groups and associations, variable candidate sources and four new wide binary systems composed of two ultracool dwarfs. Finally, we also show the spectra of 12 new ultracool dwarfs discovered by our group and presented here for the first time. These spectroscopically confirmed objects are a sanity check of our selection of ultracool dwarfs and photometric classification method.

     
    more » « less