Search

We report the results of a spectroscopic survey of candidate T subdwarfs identified by the Backyard Worlds: Planet 9 program. Near-infrared spectra of 31 sources with red $J-W2$ colors and large $J$-band reduced proper motions show varying signatures of subsolar metallicity, including strong collision-induced H$_2$ absorption, obscured methane and water features, and weak K I absorption. These metallicity signatures are supported by spectral model fits and 3D velocities, indicating thick disk and halo population membership for several sources. We identify three new metal-poor T subdwarfs ([M/H] $\lesssim$ $-$0.5), CWISE J062316.19+071505.6, WISEA J152443.14$-$262001.8, and CWISE J211250.11-052925.2; and 19 new "mild" subdwarfs with modest metal deficiency ([M/H] $\lesssim$ $-$0.25). We also identify three metal-rich brown dwarfs with thick disk kinematics. We provide kinematic evidence that the extreme L subdwarf 2MASS J053253.46+824646.5 and the mild T subdwarf CWISE J113010.07+313944.7 may be part of the Thamnos population, while the T subdwarf CWISE J155349.96+693355.2 may be part of the Helmi stream. We define a metallicity classification system for T dwarfs that adds mild subdwarfs (d/sdT), subdwarfs (sdT), and extreme subdwarfs (esdT) to the existing dwarf sequence. We also define a metallicity spectral index that correlates with metallicities inferred from spectral model fits and iron abundances from stellar primaries of benchmark T dwarf companions. This expansion of the T dwarf classification system supports investigations of ancient, metal-poor brown dwarfs now being uncovered in deep imaging and spectroscopic surveys., Comment: 82 pages, 19 figures, accepted to ApJS

Stellar mass black holes are formed from the terminal collapse of massive stars if the ensuing neutrino shock is unable to eject the stellar envelope. Direct observations of black hole formation remain inconclusive. We report observations of M31-2014-DS1, a massive, hydrogen-depleted supergiant in the Andromeda galaxy identified via a mid-infrared brightening in 2014. Its total luminosity remained nearly constant for the subsequent thousand days, before fading dramatically over the next thousand days by $\gtrsim 10\times$ and $\gtrsim 10^4\times$ in total and visible light, respectively. Together with the lack of a detected optical outburst, the observations are explained by the fallback of the stellar envelope into a newly formed black hole, moderated by the injection of a $\sim 10^{48}$ erg shock. Unifying these observations with a candidate in NGC 6946, we present a concordant picture for the birth of stellar mass black holes from stripped massive stars., Comment: Submitted for review

We report the discovery of a high velocity, very low-mass star or brown dwarf whose kinematics suggest it is unbound to the Milky Way. CWISE J124909.08+362116.0 was identified by citizen scientists in the Backyard Worlds: Planet 9 program as a high proper motion ($\mu$ $=$ 0''9/yr) faint red source. Moderate resolution spectroscopy with Keck/NIRES reveals it to be a metal-poor early L subdwarf with a large radial velocity ($-$103$\pm$10 km/s), and its estimated distance of 125$\pm$8 pc yields a speed of 456$\pm$27 km/s in the Galactic rest frame, near the local escape velocity for the Milky Way. We explore several potential scenarios for the origin of this source, including ejection from the Galactic center $\gtrsim$3 Gyr in the past, survival as the mass donor companion to an exploded white dwarf. acceleration through a three-body interaction with a black hole binary in a globular cluster, and accretion from a Milky Way satellite system. CWISE J1249+3621 is the first hypervelocity very low mass star or brown dwarf to be found, and the nearest of all such systems. It may represent a broader population of very high velocity, low-mass objects that have undergone extreme accelerations., Comment: 17 pages, 3 figures, accepted for publication in Astrophysical Journal Letters

After decades of brown dwarf discovery and follow-up, we can now infer the functional form of the mass distribution within 20 parsecs, which serves as a constraint on star formation theory at the lowest masses. Unlike objects on the main sequence that have a clear luminosity-to-mass correlation, brown dwarfs lack a correlation between an observable parameter (luminosity, spectral type, or color) and mass. A measurement of the brown dwarf mass function must therefore be procured through proxy measurements and theoretical models. We utilize various assumed forms of the mass function, together with a variety of birthrate functions, low-mass cutoffs, and theoretical evolutionary models, to build predicted forms of the effective temperature distribution. We then determine the best fit of the observed effective temperature distribution to these predictions, which in turn reveals the most likely mass function. We find that a simple power law ($dN/dM \propto M^{-\alpha}$) with $\alpha \approx 0.5$ is optimal. Additionally, we conclude that the low-mass cutoff for star formation is $\lesssim0.005M_{\odot}$. We corroborate the findings of Burgasser (2004) which state that the birthrate has a far lesser impact than the mass function on the form of the temperature distribution, but we note that our alternate birthrates tend to favor slightly smaller values of $\alpha$ than the constant birthrate. Our code for simulating these distributions is publicly available. As another use case for this code, we present findings on the width and location of the subdwarf temperature gap by simulating distributions of very old (8-10 Gyr) brown dwarfs., Comment: 13 pages, 14 figures, 1 table, accepted to ApJ

We present the discovery of VHS J183135.58$-$551355.9 (hereafter VHS J1831$-$5513), an L/T transition dwarf identified as a result of its unusually red near-infrared colors ($J-K_{\rm S}=3.633\pm0.277$ mag; $J-W2=6.249\pm0.245$ mag) from the VISTA Hemisphere Survey and CatWISE2020 surveys. We obtain low resolution near-infrared spectroscopy of VHS J1831$-$5513 using Magellan/FIRE to confirm its extremely red nature and assess features sensitive to surface gravity (i.e., youth). Its near-infrared spectrum shows multiple CH$_{\rm 4}$ absorption features, indicating an exceptionally low effective temperature for its spectral type. Based on proper motion measurements from CatWISE2020 and a photometric distance derived from its $K_{\rm S}$-band magnitude, we find that VHS J1831$-$5513 is a likely ($\sim$85$\%$ probability) kinematic member of the $\beta$ Pictoris moving group. Future radial velocity and trigonometric parallax measurements will clarify such membership. Follow-up mid-infrared or higher resolution near-infrared spectroscopy of this object will allow for further investigation as to the cause(s) of its redness, such as youth, clouds, and viewing geometry., Comment: 14 pages, 4 figures; Accepted for publication in The Astronomical Journal

We present the discovery of 13 new widely separated T dwarf companions to M dwarf primaries, identified using WISE/NEOWISE data by the CatWISE and Backyard Worlds: Planet 9 projects. This sample represents a $\sim$60% increase in the number of known M+T systems, and allows us to probe the most extreme products of binary/planetary system formation, a discovery space made available by the CatWISE2020 catalog and the Backyard Worlds: Planet 9 effort. Highlights among the sample are WISEP J075108.79-763449.6, a previously known T9 thought to be old due to its SED, which we now find is part of a common-proper-motion pair with L 34-26 A, a well studied young M3 V star within 10 pc of the Sun; CWISE J054129.32-745021.5 B and 2MASS J05581644-4501559 B, two T8 dwarfs possibly associated with the very fast-rotating M4 V stars CWISE J054129.32-745021.5 A and 2MASS J05581644-4501559 A; and UCAC3 52-1038 B, which is among the widest late T companions to main sequence stars, with a projected separation of $\sim$7100 au. The new benchmarks presented here are prime $JWST$ targets, and can help us place strong constraints on formation and evolution theory of substellar objects as well as on atmospheric models for these cold exoplanet analogs., Comment: Accepted for publication in ApJ. 35 pages, 6 tables, 21 figures

We report the identification of 89 new systems containing ultracool dwarf companions to main sequence stars and white dwarfs, using the citizen science project Backyard Worlds: Planet 9 and cross-reference between Gaia and CatWISE2020. Thirty-two of these companions and thirty-three host stars were followed up with spectroscopic observations, with companion spectral types ranging from M7-T9 and host spectral types ranging from G2-M9. These systems exhibit diverse characteristics, from young to old ages, blue to very red spectral morphologies, potential membership to known young moving groups, and evidence of spectral binarity in 9 companions. Twenty of the host stars in our sample show evidence for higher order multiplicity, with an additional 11 host stars being resolved binaries themselves. We compare this sample's characteristics with those of the known stellar binary and exoplanet populations, and find our sample begins to fill in the gap between directly imaged exoplanets and stellary binaries on mass ratio-binding energy plots. With this study, we increase the population of ultracool dwarf companions to FGK stars by $\sim$42\%, and more than triple the known population of ultracool dwarf companions with separations larger than 1,000 au, providing excellent targets for future atmospheric retrievals., Comment: 61 pages, 11 figures, 11 tables. Accepted for publication in AJ

We present observations of the T8 dwarf 2MASS 0415-0935 with JWST's NIRSpec spectrograph using the G395H grating ($\sim$ 2.87 - 5.14 $\mu$m). We perform the first atmospheric retrieval analysis at the maximum spectral resolution of NIRSpec (R$\sim$2700) and combine the spectrum with previous observations to study the 0.9-20 $\mu$m spectral energy distribution. We obtain precise constraints on chemical abundances ($\sim$0.02 dex) for a number of species which complicate our understanding of disequilibrium chemistry, particularly for CO$_{2}$ and PH$_{3}$. Furthermore, we measure a $^{12}$CO/$^{13}$CO ratio of $\sim 97^{+9}_{-8}$, making 2MASS 0415-0935 the coldest ($\sim 760$ K) substellar object outside of our solar system with a measured $^{12}$CO/$^{13}$CO ratio. This work shows promise for similar observations with JWST to provide precise abundances of major chemical species as well as isotopologues, allowing for new tests of our understanding of the formation and atmospheres of substellar objects., Comment: 28 pages, 10 figures. Submitted to Nature Astronomy on Oct. 5th, 2023

Most tidal disruption events (TDEs) are currently found in time-domain optical and soft X-ray surveys, both of which are prone to significant obscuration. The infrared (IR), however, is a powerful probe of dust-enshrouded environments, and hence, we recently performed a systematic search of NEOWISE mid-IR data for nearby, obscured TDEs within roughly 200 Mpc. We identified 18 TDE candidates in galactic nuclei, using difference imaging to uncover nuclear variability amongst significant host galaxy emission. These candidates were selected based on the following IR light curve properties: (1) $L_\mathrm{W2}\gtrsim10^{42}$ erg s$^{-1}$ at peak, (2) fast rise, followed by a slow, monotonic decline, (3) no significant prior variability, and (4) no evidence for AGN activity in WISE colors. The majority of these sources showed no variable optical counterpart, suggesting that optical surveys indeed miss numerous obscured TDEs. Using narrow line ionization levels and variability arguments, we identified 6 sources as possible underlying AGN, yielding a total of 12 TDEs in our gold sample. This gold sample yields a lower limit on the IR-selected TDE rate of $2.0\pm0.3\times10^{-5}$ galaxy$^{-1}$ year$^{-1}$ ($1.3\pm0.2\times10^{-7}$ Mpc$^{-3}$ year$^{-1}$), which is comparable to optical and X-ray TDE rates. The IR-selected TDE host galaxies do not show a green valley overdensity nor a preference for quiescent, Balmer strong galaxies, which are both overrepresented in optical and X-ray TDE samples. This IR-selected sample represents a new population of dusty TDEs that have historically been missed by optical and X-ray surveys and helps alleviate tensions between observed and theoretical TDE rates and the so-called missing energy problem., Comment: 28 pages, 13 figures (+ 7 pages, 4 figures in appendix), accepted for publication in ApJ

We introduce the Wide-field Retrieval of Astrodata Program (WRAP), a tool created to aid astronomers in gathering photometric and astrometric data for point sources that may confuse simple cross-matching algorithms because of their faintness or motion. WRAP allows astronomers to correctly cross-identify objects with proper motion across multiple surveys by wedding the catalog data with its underlying images, thus providing visual confirmation of cross-associations in real time. Developed within the Backyard Worlds: Planet 9 citizen science project, WRAP aims to aid in the characterization of faint, high motion sources by this collaboration (and others)., Comment: 4 pages, 1 figure, Github and DOI link to program

A complete accounting of nearby objects -- from the highest-mass white dwarf progenitors down to low-mass brown dwarfs -- is now possible, thanks to an almost complete set of trigonometric parallax determinations from Gaia, ground-based surveys, and Spitzer follow-up. We create a census of objects within a Sun-centered sphere of 20-pc radius and check published literature to decompose each binary or higher-order system into its separate components. The result is a volume-limited census of $\sim$3,600 individual star formation products useful in measuring the initial mass function across the stellar ($<8 M_\odot$) and substellar ($\gtrsim 5 M_{Jup}$) regimes. Comparing our resulting initial mass function to previous measurements shows good agreement above 0.8$M_\odot$ and a divergence at lower masses. Our 20-pc space densities are best fit with a quadripartite power law, $\xi(M) = dN/dM \propto M^{-\alpha}$ with long-established values of $\alpha = 2.3$ at high masses ($0.55 < M < 8.00 M_\odot$) and $\alpha = 1.3$ at intermediate masses ($0.22 < M < 0.55 M_\odot$), but at lower masses we find $\alpha = 0.25$ for $0.05 < M <0.22 M_\odot$ and $\alpha = 0.6$ for $0.01 < M < 0.05 M_\odot$. This implies that the rate of production as a function of decreasing mass diminishes in the low-mass star/high-mass brown dwarf regime before increasing again in the low-mass brown dwarf regime. Correcting for completeness, we find a star to brown dwarf number ratio of, currently, 4:1, and an average mass per object of 0.41 $M_\odot$., Comment: Accepted for publication in The Astrophysical Journal Supplement Series. 123 pages with four ancillary files

Y dwarfs, the coolest known spectral class of brown dwarfs, overlap in mass and temperature with giant exoplanets, providing unique laboratories for studying low-temperature atmospheres. However, only a fraction of Y dwarf candidates have been spectroscopically confirmed. We present Keck/NIRES near-infrared spectroscopy of the nearby ($d \approx 6-8$ pc) brown dwarf CWISE J105512.11+544328.3. Although its near-infrared spectrum aligns best with the Y0 standard in the $J$-band, no standard matches well across the full $YJHK$ wavelength range. The CWISE J105512.11+544328.3 NH$_3$-$H$ = 0.427 $\pm$ 0.0012 and CH$_4$-$J$ = 0.0385 $\pm$ 0.0007 absorption indices and absolute Spitzer [4.5] magnitude of 15.18 $\pm$ 0.22 are also indicative of an early Y dwarf rather than a late T dwarf. CWISE J105512.11+544328.3 additionally exhibits the bluest Spitzer [3.6]$-$[4.5] color among all spectroscopically confirmed Y dwarfs. Despite this anomalously blue Spitzer color given its low luminosity, CWISE J105512.11+544328.3 does not show other clear kinematic or spectral indications of low metallicity. Atmospheric model comparisons yield a log(g) $\le$ 4.5 and $T_{\rm eff} \approx 500 \pm 150$ K for this source. We classify CWISE J105512.11+544328.3 as a Y0 (pec) dwarf, adding to the remarkable diversity of the Y-type population. JWST spectroscopy would be crucial to understanding the origin of this Y dwarf's unusual preference for low-gravity models and blue 3-5 $\mu$m color., Comment: accepted for publication in The Astrophysical Journal

We present three new brown dwarf spectral binary candidates: CWISE J072708.09$-$360729.2, CWISE J103604.84$-$514424.4, and CWISE J134446.62$-$732053.9, discovered by citizen scientists through the Backyard Worlds: Planet 9 project. Follow-up near-infrared spectroscopy shows that each of these objects is poorly fit by a single near-infrared standard. We constructed binary templates and found significantly better fits, with component types of L7+T4 for CWISE J072708.09$-$360729.2, L7+T4 for CWISE J103604.84$-$514424.4, and L7+T7 for CWISE J134446.62$-$732053.9. However, further investigation of available spectroscopic indices for evidence of binarity and large amplitude variability suggests that CWISE J072708.09$-$360729.2 may instead be a strong variability candidate. Our analysis offers tentative evidence and characterization of these peculiar brown dwarf sources, emphasizing their value as promising targets for future high-resolution imaging or photometric variability studies., Comment: Accepted to the Astronomical Journal

The Dark Energy Spectroscopic Instrument, consisting of 5020 robotic fiber positioners and associated systems on the Mayall telescope at Kitt Peak, Arizona, is carrying out a survey to measure the spectra of 40 million galaxies and quasars and produce the largest 3D map of the universe to date. The primary science goal is to use baryon acoustic oscillations to measure the expansion history of the universe and the time evolution of dark energy. A key function of the online control system is to position each fiber on a particular target in the focal plane with an accuracy of 11$\mu$m rms 2-D. This paper describes the set of software programs used to perform this function along with the methods used to validate their performance., Comment: 27 pages, 16 figures

Millions of quasar spectra will be collected by the Dark Energy Spectroscopic Instrument (DESI), leading to a four-fold increase in the number of known quasars. High accuracy quasar classification is essential to tighten constraints on cosmological parameters measured at the highest redshifts DESI observes ($z>2.0$). We present the spectral templates for identification and redshift estimation of quasars in the DESI Year 1 data release. The quasar templates are comprised of two quasar eigenspectra sets, trained on spectra from the Sloan Digital Sky Survey. The sets are specialized to reconstruct quasar spectral variation observed over separate yet overlapping redshift ranges and, together, are capable of identifying DESI quasars from $0.05 < z <7.0$. The new quasar templates show significant improvement over the previous DESI quasar templates regarding catastrophic failure rates, redshift precision and accuracy, quasar completeness, and the contamination fraction in the final quasar sample., Comment: to be published in Astronomical Journal; 21 pages, 8 figures

We detect the cross-correlation between 2.7 million DESI quasar targets across 14,700 deg$^2$ (180 quasars deg$^{-2}$) and Planck 2018 CMB lensing at $\sim$30$\sigma$. We use the cross-correlation on very large scales to constrain local primordial non-Gaussianity via the scale dependence of quasar bias. The DESI quasar targets lie at an effective redshift of 1.51 and are separated into four imaging regions of varying depth and image quality. We select quasar targets from Legacy Survey DR9 imaging, apply additional flux and photometric redshift cuts to improve the purity and reduce the fraction of unclassified redshifts, and use early DESI spectroscopy of 194,000 quasar targets to determine their redshift distribution and stellar contamination fraction (2.6%). Due to significant excess large-scale power in the quasar autocorrelation, we apply weights to mitigate contamination from imaging systematics such as depth, extinction, and stellar density. We use realistic contaminated mocks to determine the greatest number of systematic modes that we can fit, before we are biased by overfitting and spuriously remove real power. We find that linear regression with one to seven imaging templates removed per region accurately recovers the input cross-power, $f_{\textrm{NL}}$ and linear bias. As in previous analyses, our $f_{\textrm{NL}}$ constraint depends on the linear primordial non-Gaussianity bias parameter, $b_{\phi} = 2(b - p)\delta_c$ assuming universality of the halo mass function. We measure $f_{\textrm{NL}} = -26^{+45}_{-40}$ with $p=1.6$ $(f_{\textrm{NL}} = -18^{+29}_{-27}$ with $p=1.0$), and find that this result is robust under several systematics tests. Future spectroscopic quasar cross-correlations with Planck lensing lensing can tighten the $f_{\textrm{NL}}$ constraints by a factor of 2 if they can remove the excess power on large scales in the quasar auto power spectrum., Comment: 57 pages, 25 figures, updated to match accepted version

Using a sample of 361 nearby brown dwarfs, we have searched for 4.6$\mu$m variability indicative of large-scale rotational modulations or large-scale long-term changes on timescales of over 10 years. Our findings show no statistically significant variability in \textit{Spitzer} ch2 or \textit{WISE} W2 photometry. For \textit{Spitzer} the ch2 1$\sigma$ limits are $\sim$8 mmag for objects at 11.5 mag and $\sim$22 mmag for objects at 16 mag. This corresponds to no variability above 4.5$\%$ at 11.5 mag and 12.5$\%$ at 16 mag. We conclude that highly variable brown dwarfs, at least two previously published examples of which have been shown to have 4.6$\mu$m variability above 80 mmag, are very rare. While analyzing the data, we also developed a new technique for identifying brown dwarfs binary candidates in \textit{Spitzer} data. We find that known binaries have IRAC ch2 PRF (point response function) flux measurements that are consistently dimmer than aperture flux measurements. We have identified 59 objects that exhibit such PRF versus apertures flux differences and are thus excellent binary brown dwarf candidates., Comment: 18 pages, 11 figures, accepted for publication in AJ, and 1 figure set of 361 figures with the example being figure 8

The Galactic nova rate is intimately linked to our understanding of its chemical enrichment and progenitor channels of Type Ia supernovae. Yet past estimates have varied by more than an order of magnitude ($\approx10-300$ yr$^{-1}$) owing to limitations in both discovery methods as well as assumptions regarding the Galactic dust distribution and extragalactic stellar populations. Recent estimates utilizing synoptic near-infrared surveys have begun to provide a glimpse of a consensus ($\approx25-50$ yr$^{-1}$); however, a consistent estimate remains lacking. Here, we present the first all-sky search for Galactic novae using 8 years of data from the NEOWISE mid-infrared (MIR) survey. Operating at $3.4$ and $4.6$ $\mu$m where interstellar extinction is negligible, the 6-month cadence NEOWISE dataset offers unique sensitivity to discover slowly evolving novae across the entire Galaxy. Using a novel image subtraction pipeline together with systematic selection criteria, we identify a sample of 49 rapidly evolving MIR outbursts as candidate Galactic novae. While 27 of these sources are known novae, the remaining are previously missed nova candidates discovered in this work. The unknown events are spatially clustered along the densest and most heavily obscured regions of the Galaxy where previous novae are severely underrepresented. We use simulations of the NEOWISE survey strategy, the pipeline detection efficiency, and our criteria to derive a Galactic nova rate of $47.9^{+3.1}_{-8.3}$ yr$^{-1}$. The discovery of these exceptionally bright (yet overlooked) nova candidates confirm emerging suggestions that optical surveys have been highly incomplete in searches for Galactic novae, highlighting the potential for MIR searches in revealing the demographics of Galactic stellar outbursts., Comment: 13 pages, 4 figures, 1 table, 3 figures in Appendix. Submitted to MNRAS. Comments welcome

While the vast majority of Tidal Disruption Events (TDEs) has been identified by wide-field sky surveys in the optical and X-ray bands, recent studies indicate that a considerable fraction of TDEs may be dust obscured, and thus preferentially detected in the infrared (IR) wavebands. In this Letter, we present the discovery of a luminous mid-IR nuclear flare (termed WTP 14adbjsh) identified in a systematic transient search of archival images from the NEOWISE mid-IR survey. The source reached a peak luminosity of $L \simeq 10^{43} \text{erg s}^{-1}$ at 4.6 ${\mu}$m in 2015, before fading in the IR with a TDE-like $F \propto t^{-5/3}$ decline, radiating a total of more than $ 3\times 10^{51}$ erg in the last 7 years. The transient event took place in the nearby galaxy NGC 7392, at a distance of around 42 Mpc; yet, no optical or X-ray flare is detected. We interpret the transient as the nearest TDE candidate detected in the last decade, which was missed at other wavelengths due to dust obscuration, hinting at the existence of TDEs that have been historically overlooked. Unlike most previously detected TDEs, the transient was discovered in a star forming galaxy, corroborating earlier suggestions that dust obscuration suppresses significantly the detection of TDEs in these environments. Our results demonstrate that the study of IR-detected TDEs is critical in order to obtain a complete understanding of the physics of TDEs, and to conclude whether TDEs occur preferentially in a particular class of galaxies., Comment: 13 pages, 4 figures, accepted for publication by ApJL

A new class of white dwarfs, dubbed DAHe, that present Zeeman-split Balmer lines in emission has recently emerged. However, the physical origin of these emission lines remains unclear. We present here a sample of 21 newly identified DAHe systems and determine magnetic field strengths and (for a subset) periods which span the ranges of ~ 6.5 -- 147 MG and ~ 0.4 -- 36 h respectively. All but four of these systems were identified from the Dark Energy Spectroscopic Instrument (DESI) survey sample of more than 47000 white dwarf candidates observed during its first year of observations. We present detailed analysis of the new DAHe WDJ161634.36+541011.51 with a spin period of 95.3 min, which exhibits an anti-correlation between broadband flux and Balmer line strength that is typically observed for this class of systems. All DAHe systems cluster closely on the Gaia Hertzsprung-Russell diagram where they represent ~ 1 per cent of white dwarfs within that region. This grouping further solidifies their unexplained emergence at relatively late cooling times and we discuss this in context of current formation theories. Nine of the new DAHe systems are identifiable from SDSS spectra of white dwarfs that had been previously classified as featureless DC-type systems. We suggest high S/N, unbiased observations of DCs as a possible route for discovering additional DAHe systems., Comment: 19 pages, 10 Figures, accepted for publication in MNRAS

We present the discovery of CWISE J050626.96$+$073842.4 (CWISE J0506$+$0738), an L/T transition dwarf with extremely red near-infrared colors discovered through the Backyard Worlds: Planet 9 citizen science project. Photometry from UKIRT and CatWISE give a $(J-K)_{\rm MKO}$ color of 2.97$\pm$0.03 mag and a $J_{\rm MKO}-$W2 color of 4.93$\pm$0.02 mag, making CWISE J0506$+$0738 the reddest known free-floating L/T dwarf in both colors. We confirm the extremely red nature of CWISE J0506$+$0738 using Keck/NIRES near-infrared spectroscopy and establish that it is a low-gravity late-type L/T transition dwarf. The spectrum of CWISE J0506$+$0738 shows possible signatures of CH$_4$ absorption in its atmosphere, suggesting a colder effective temperature than other known, young, red L dwarfs. We assign a preliminary spectral type for this source of L8$\gamma$-T0$\gamma$. We tentatively find that CWISE J0506$+$0738 is variable at 3-5 $\mu$m based on multi-epoch WISE photometry. Proper motions derived from follow-up UKIRT observations combined with a radial velocity from our Keck/NIRES spectrum and a photometric distance estimate indicate a strong membership probability in the $\beta$ Pic moving group. A future parallax measurement will help to establish a more definitive moving group membership for this unusual object., Comment: Accepted to ApJ Letters

We introduce the DESI LOW-Z Secondary Target Survey, which combines the wide-area capabilities of the Dark Energy Spectroscopic Instrument (DESI) with an efficient, low-redshift target selection method. Our selection consists of a set of color and surface brightness cuts, combined with modern machine learning methods, to target low-redshift dwarf galaxies ($z$ < 0.03) between $19 < r < 21$ with high completeness. We employ a convolutional neural network (CNN) to select high-priority targets. The LOW-Z survey has already obtained over 22,000 redshifts of dwarf galaxies (M$_* < 10^9$ M$_\odot$), comparable to the number of dwarf galaxies discovered in SDSS-DR8 and GAMA. As a spare fiber survey, LOW-Z currently receives fiber allocation for just ~50% of its targets. However, we estimate that our selection is highly complete: for galaxies at $z < 0.03$ within our magnitude limits, we achieve better than 95% completeness with ~1% efficiency using catalog-level photometric cuts. We also demonstrate that our CNN selections $z<0.03$ galaxies from the photometric cuts subsample at least ten times more efficiently while maintaining high completeness. The full five-year DESI program will expand the LOW-Z sample, densely mapping the low-redshift Universe, providing an unprecedented sample of dwarf galaxies, and providing critical information about how to pursue effective and efficient low-redshift surveys., Comment: 24 pages, 14 figures, data to reproduce figures: https://zenodo.org/record/7422591

We present the discovery of CWISE J151044.74$-$524923.5, a wide low-mass companion to the nearby ($\sim$24.7 pc) system L 262-74, which was identified through the Backyard Worlds: Planet 9 citizen science project. We detail the properties of the system, and we assess that this companion is a mid-L dwarf, which will need to be verified spectroscopically. With an angular separation of 74\farcs3, we estimate a projected physical separation of $\sim$1837 au from the central system.

We describe the spectroscopic data processing pipeline of the Dark Energy Spectroscopic Instrument (DESI), which is conducting a redshift survey of about 40 million galaxies and quasars using a purpose-built instrument on the 4-m Mayall Telescope at Kitt Peak National Observatory. The main goal of DESI is to measure with unprecedented precision the expansion history of the Universe with the Baryon Acoustic Oscillation technique and the growth rate of structure with Redshift Space Distortions. Ten spectrographs with three cameras each disperse the light from 5000 fibers onto 30 CCDs, covering the near UV to near infrared (3600 to 9800 Angstrom) with a spectral resolution ranging from 2000 to 5000. The DESI data pipeline generates wavelength- and flux-calibrated spectra of all the targets, along with spectroscopic classifications and redshift measurements. Fully processed data from each night are typically available to the DESI collaboration the following morning. We give details about the pipeline's algorithms, and provide performance results on the stability of the optics, the quality of the sky background subtraction, and the precision and accuracy of the instrumental calibration. This pipeline has been used to process the DESI Survey Validation data set, and has exceeded the project's requirements for redshift performance, with high efficiency and a purity greater than 99 percent for all target classes., Comment: AJ, revised version, 55 pages, 55 figures, 4 tables

A key component of the Dark Energy Spectroscopic Instrument (DESI) survey validation (SV) is a detailed visual inspection (VI) of the optical spectroscopic data to quantify key survey metrics. In this paper we present results from VI of the quasar survey using deep coadded SV spectra. We show that the majority (~70%) of the main-survey targets are spectroscopically confirmed as quasars, with ~16% galaxies, ~6% stars, and ~8% low-quality spectra lacking reliable features. A non-negligible fraction of the quasars are misidentified by the standard spectroscopic pipeline but we show that the majority can be recovered using post-pipeline "afterburner" quasar-identification approaches. We combine these "afterburners" with our standard pipeline to create a modified pipeline to improve the overall quasar yield. At the depth of the main DESI survey both pipelines achieve a good-redshift purity (reliable redshifts measured within 3000 km/s) of ~99%; however, the modified pipeline recovers ~94% of the visually inspected quasars, as compared to ~86% from the standard pipeline. We demonstrate that both pipelines achieve an median redshift precision and accuracy of ~100 km/s and ~70 km/s, respectively. We constructed composite spectra to investigate why some quasars are missed by the standard spectroscopic pipeline and find that they are more host-galaxy dominated (i.e., distant analogs of "Seyfert galaxies") and/or dust reddened than the standard-pipeline quasars. We also show example spectra to demonstrate the overall diversity of the DESI quasar sample and provide strong-lensing candidates where two targets contribute to a single spectrum., Comment: Astronomical journal (in press). 26 pages, 15 figures, 9 tables, one of a suite of 8 papers detailing targeting for DESI. Figure data available from Zenodo (see paper for details)

The Dark Energy Spectroscopic Instrument (DESI) Survey has obtained a set of spectroscopic measurements of galaxies to validate the final survey design and target selections. To assist in these tasks, we visually inspect (VI) DESI spectra of approximately 2,500 bright galaxies, 3,500 luminous red galaxies (LRGs), and 10,000 emission line galaxies (ELGs), to obtain robust redshift identifications. We then utilize the VI redshift information to characterize the performance of the DESI operation. Based on the VI catalogs, our results show that the final survey design yields samples of bright galaxies, LRGs, and ELGs with purity greater than $99\%$. Moreover, we demonstrate that the precision of the redshift measurements is approximately 10 km/s for bright galaxies and ELGs and approximately 40 km/s for LRGs. The average redshift accuracy is within 10 km/s for the three types of galaxies. The VI process also helps improve the quality of the DESI data by identifying spurious spectral features introduced by the pipeline. Finally, we show examples of unexpected real astronomical objects, such as Ly$\alpha$ emitters and strong lensing candidates, identified by VI. These results demonstrate the importance and utility of visually inspecting data from incoming and upcoming surveys, especially during their early operation phases., Comment: 22 pages, 14 figures, 3 tables, one of a suite of 8 papers detailing targeting for DESI. ApJ accepted version with minor textual updates

The Dark Energy Spectroscopic Instrument (DESI) is carrying out a 5-year survey that aims to measure the redshifts of tens of millions of galaxies and quasars, including 8 million luminous red galaxies (LRGs) in the redshift range of $0.4

We describe the Milky Way Survey (MWS) that will be undertaken with the Dark Energy Spectroscopic Instrument (DESI) on the Mayall 4m telescope at the Kitt Peak National Observatory. Over the next 5 yr DESI MWS will observe approximately seven million stars at Galactic latitudes |b|>20 degrees, with an inclusive target selection scheme focused on the thick disk and stellar halo. MWS will also include several high-completeness samples of rare stellar types, including white dwarfs, low-mass stars within 100pc of the Sun, and horizontal branch stars. We summarize the potential of DESI to advance understanding of Galactic structure and stellar evolution. We introduce the final definitions of the main MWS target classes and estimate the number of stars in each class that will be observed. We describe our pipelines for deriving radial velocities, atmospheric parameters, and chemical abundances. We use ~500,000 spectra of unique stellar targets from the DESI Survey Validation program (SV) to demonstrate that our pipelines can measure radial velocities to ~1 km/s and [Fe/H] accurate to ~0.2 dex for typical stars in our main sample. We find the stellar parameter distributions from ~100 sq. deg of SV observations with >90% completeness on our main sample are in good agreement with expectations from mock catalogs and previous surveys., Comment: Accepted for publication in ApJ, 44 pages, 25 figures, 4 tables, one of a suite of 8 papers detailing targeting for DESI; v2 added links to data shown in figures, added citations to other DESI papers, corrected author list and minor typos; v3 fixed minor errors in Fig. 6 and clarified associated text; v4 updated to include minor changes in response to review

The Dark Energy Spectroscopic Instrument (DESI) will precisely constrain cosmic expansion and the growth of structure by collecting $\sim$40 million extra-galactic redshifts across $\sim$80\% of cosmic history and one third of the sky. The Emission Line Galaxy (ELG) sample, which will comprise about one-third of all DESI tracers, will be used to probe the Universe over the $0.6 < z < 1.6$ range, which includes the $1.1

The Dark Energy Spectroscopic Instrument (DESI) survey will measure large-scale structures using quasars as direct tracers of dark matter in the redshift range 0.92.1. We present several methods to select candidate quasars for DESI, using input photometric imaging in three optical bands (g, r, z) from the DESI Legacy Imaging Surveys and two infrared bands (W1, W2) from the Wide-field Infrared Explorer (WISE). These methods were extensively tested during the Survey Validation of DESI. In this paper, we report on the results obtained with the different methods and present the selection we optimized for the DESI main survey. The final quasar target selection is based on a Random Forest algorithm and selects quasars in the magnitude range 16.599% purity for a nominal effective exposure time of ~1000s. With a 310 per sq. deg. target density, the main selection allows DESI to select more than 200 QSOs per sq. deg. (including 60 quasars with z>2.1), exceeding the project requirements by 20%. The redshift distribution of the selected quasars is in excellent agreement with quasar luminosity function predictions., Comment: 21 pages, 21 figures, typos corrected, references added

The Dark Energy Spectroscopic Instrument (DESI) has embarked on an ambitious five-year survey to explore the nature of dark energy with spectroscopy of 40 million galaxies and quasars. DESI will determine precise redshifts and employ the Baryon Acoustic Oscillation method to measure distances from the nearby universe to z > 3.5, as well as measure the growth of structure and probe potential modifications to general relativity. In this paper we describe the significant instrumentation we developed for the DESI survey. The new instrumentation includes a wide-field, 3.2-deg diameter prime-focus corrector that focuses the light onto 5020 robotic fiber positioners on the 0.812 m diameter, aspheric focal surface. The positioners and their fibers are divided among ten wedge-shaped petals. Each petal is connected to one of ten spectrographs via a contiguous, high-efficiency, nearly 50 m fiber cable bundle. The ten spectrographs each use a pair of dichroics to split the light into three channels that together record the light from 360 - 980 nm with a resolution of 2000 to 5000. We describe the science requirements, technical requirements on the instrumentation, and management of the project. DESI was installed at the 4-m Mayall telescope at Kitt Peak, and we also describe the facility upgrades to prepare for DESI and the installation and functional verification process. DESI has achieved all of its performance goals, and the DESI survey began in May 2021. Some performance highlights include RMS positioner accuracy better than 0.1", SNR per \sqrt{\AA} > 0.5 for a z > 2 quasar with flux 0.28e-17 erg/s/cm^2/A at 380 nm in 4000s, and median SNR = 7 of the [OII] doublet at 8e-17 erg/s/cm^2 in a 1000s exposure for emission line galaxies at z = 1.4 - 1.6. We conclude with highlights from the on-sky validation and commissioning of the instrument, key successes, and lessons learned. (abridged), Comment: 78 pages, 32 figures, submitted to AJ

We present deep, full-sky maps built from Wide-field Infrared Survey Explorer (WISE) and NEOWISE exposures spanning the 2010 January - 2020 December time period. These coadds, which incorporate roughly 8 years of W1 (3.4 microns) and W2 (4.6 microns) imaging, are the deepest ever full-sky maps at wavelengths of 3-5 microns. Photometry based on these coadds will be a component of DESI Legacy Imaging Surveys DR10., Comment: published in RNAAS

While stars are often found in binary systems, brown dwarf binaries are much rarer. Brown dwarf--brown dwarf pairs are typically difficult to resolve because they often have very small separations. Using brown dwarfs discovered with data from the Wide-field Infrared Survey Explorer (WISE) via the Backyard Worlds: Planet 9 citizen science project, we inspected other, higher resolution, sky surveys for overlooked cold companions. During this process we discovered the brown dwarf binary system CWISE J0146$-$0508AB, which we find has a very small chance alignment probability based on the similar proper motions of the components of the system. Using follow-up near-infrared spectroscopy with Keck/NIRES, we determined component spectral types of L4 and L8 (blue), making CWISE J0146$-$0508AB one of only a few benchmark systems with a blue L dwarf. At an estimated distance of $\sim$40 pc, CWISE J0146$-$0508AB has a projected separation of $\sim$129 AU, making it the widest separation brown dwarf pair found to date. We find that such a wide separation for a brown dwarf binary may imply formation in a low-density star-forming region., Comment: Accepted for publication in ApJ Letters

In an effort to identify nearby and unusual cold objects in the solar neighborhood, we searched for previously unidentified moving objects using CatWISE2020 proper motion data combined with machine learning methods. We paired the motion candidates with their counterparts in 2MASS, UHS, and VHS. Then we searched for white dwarf, brown dwarf, and subdwarf outliers on the resulting color-color diagrams. This resulted in the discovery of 16 new dwarfs including two nearby M dwarfs (< 30 pc), a possible young L dwarf, a high motion early T dwarf and 3 later T dwarfs. This research represents a step forward in completing the census of the Sun's neighbors., Comment: 15 pages, 7 Figures, Submitted to The Astronomical Journal

Through the Backyard Worlds: Planet 9 citizen science project we discovered a late-type L dwarf co-moving with the young K0 star BD+60 1417 at a projected separation of 37" or 1662 AU. The secondary - CWISER J124332.12+600126.2 (W1243) - is detected in both the CatWISE2020 and 2MASS reject tables. The photometric distance and CatWISE proper motion both match that of the primary within ~1sigma and our estimates for chance alignment yield a zero probability. Follow-up near infrared spectroscopy reveals W1243 to be a very red 2MASS color(J-Ks=2.72), low-surface gravity source that we classify as L6 - L8gamma. Its spectral morphology strongly resembles that of confirmed late-type L dwarfs in 10 - 150 Myr moving groups as well as that of planetary mass companions. The position on near- and mid-infrared color-magnitude diagrams indicates the source is redder and fainter than the field sequence, a telltale sign of an object with thick clouds and a complex atmosphere. For the primary we obtained new optical spectroscopy and analyzed all available literature information for youth indicators. We conclude that the Li I abundance, its loci on color-magnitude and color-color diagrams, and the rotation rate revealed in multiple TESS sectors are all consistent with an age of 50 - 150 Myr. Using our re-evaluated age of the primary, the Gaia parallax along with the photometry and spectrum for W1243 we find a Teff=1303+/-31 K, logg=4.3+/-0.17 cm s-2, and a mass of 15+/-5 MJup. We find a physical separation of ~1662 AU and a mass ratio of ~0.01 for this system. Placing it in context with the diverse collection of binary stars, brown dwarf and planetary companions, the BD+60 1417 system falls in a sparsely sampled area where the formation pathway is difficult to assess., Comment: Accepted for publication in ApJ

We present the discovery of CWISE J052306.42$-$015355.4, which was found as a faint, significant proper motion object (0.52 $\pm$ 0.08 arcsec yr$^{-1}$) using machine learning tools on the unWISE re-processing on time series images from the Wide-field Infrared Survey Explorer. Using the CatWISE2020 W1 and W2 magnitudes along with a $J-$band detection from the VISTA Hemisphere Survey, the location of CWISE J052306.42$-$015355.4 on the W1$-$W2 vs. $J-$W2 diagram best matches that of other known, or suspected, extreme T subdwarfs. As there is currently very little knowledge concerning extreme T subdwarfs we estimate a rough distance of $\le$ 68 pc, which results in a tangential velocity of $\le$ 167 km s$^{-1}$, both of which are tentative. A measured parallax is greatly needed to test these values. We also estimate a metallicity of $-1.5 <$ [M/H] $< -0.5$ using theoretical predictions., Comment: 7 pages, 3 figures, accepted for publication in the AJ

Through the Backyard Worlds: Planet 9 citizen science project, we have identified a wide-separation ($\sim$10', $\sim$9900 au projected) substellar companion to the nearby ($\sim$17.5 pc), mid-M dwarf Ross 19. We have developed a new formalism for determining chance alignment probabilities based on the BANYAN $\Sigma$ tool, and find a 100% probability that this is a physically associated pair. Through a detailed examination of Ross 19A, we find that the system is metal-poor ([Fe/H]=$-$0.40$\pm$0.12) with an age of 7.2$^{+3.8}_{-3.6}$ Gyr. Combining new and existing photometry and astrometry, we find that Ross 19B is one of the coldest known wide-separation companions, with a spectral type on the T/Y boundary, an effective temperature of 500$^{+115}_{-100}$ K, and a mass in the range 15-40 $M_{\rm Jup}$. This new, extremely cold benchmark companion is a compelling target for detailed characterization with future spectroscopic observations using facilities such as the Hubble Space Telescope or James Webb Space Telescope., Comment: Accepted to ApJ

Continued follow-up of WISEA J153429.75-104303.3, announced in Meisner et al (2020), has proven it to have an unusual set of properties. New imaging data from Keck/MOSFIRE and HST/WFC3 show that this object is one of the few faint proper motion sources known with J-ch2 > 8 mag, indicating a very cold temperature consistent with the latest known Y dwarfs. Despite this, it has W1-W2 and ch1-ch2 colors ~1.6 mag bluer than a typical Y dwarf. A new trigonometric parallax measurement from a combination of WISE, Spitzer, and HST astrometry confirms a nearby distance of $16.3^{+1.4}_{-1.2}$ pc and a large transverse velocity of $207.4{\pm}15.9$ km/s. The absolute J, W2, and ch2 magnitudes are in line with the coldest known Y dwarfs, despite the highly discrepant W1-W2 and ch1-ch2 colors. We explore possible reasons for the unique traits of this object and conclude that it is most likely an old, metal-poor brown dwarf and possibly the first Y subdwarf. Given that the object has an HST F110W magnitude of 24.7 mag, broad-band spectroscopy and photometry from JWST are the best options for testing this hypothesis., Comment: 8 pages, 4 figures, accepted for publication in The Astrophysical Journal Letters

Schneider et al. (2020) presented the discovery of WISEA J041451.67-585456.7 and WISEA J181006.18-101000.5, which appear to be the first examples of extreme T-type subdwarfs (esdTs; metallicity <= -1 dex, T_eff <= 1400 K). Here we present new discoveries and follow-up of three T-type subdwarf candidates, with an eye toward expanding the sample of such objects with very low metallicity and extraordinarily high kinematics, properties that suggest membership in the Galactic halo. Keck/NIRES near-infrared spectroscopy of WISEA J155349.96+693355.2, a fast-moving object discovered by the Backyard Worlds: Planet 9 citizen science project, confirms that it is a mid-T subdwarf. With H_W2 = 22.3 mag, WISEA J155349.96+693355.2 has the largest W2 reduced proper motion among all spectroscopically confirmed L and T subdwarfs, suggesting that it may be kinematically extreme. Nevertheless, our modeling of the WISEA J155349.96+693355.2 near-infrared spectrum indicates that its metallicity is only mildly subsolar. In analyzing the J155349.96+693355.2 spectrum, we present a new grid of low-temperature, low-metallicity model atmosphere spectra. We also present the discoveries of two new esdT candidates, CWISE J073844.52-664334.6 and CWISE J221706.28-145437.6, based on their large motions and colors similar to those of the two known esdT objects. Finding more esdT examples is a critical step toward mapping out the spectral sequence and observational properties of this newly identified population., Comment: accepted for publication in ApJ; models available at https://doi.org/10.7910/DVN/SJRXUO

We present the discovery of CWISE J203546.35-493611.0, a peculiar M8 companion to the M4.5 star APMPM J2036-4936 discovered through the citizen science project Backyard Worlds: Planet 9. Given CWISE J203546.35-493611.0's proper motion ($\mu_{\alpha}$, $\mu_{\delta}$) = ($-$126$\pm$22, $-$478$\pm$23) and angular separation of 34.2$''$ from APMPM 2036-4936, we calculate a chance alignment probability of $1.15 \times 10^{-6}$. Both stars in this system appear to be underluminous, and the spectrum obtained for CWISE J203546.35-493611.0 shows a triangular H band. Further study of this system is warranted to understand these peculiarities., Comment: 3 pages, 1 figure