An Investigation of New Brown Dwarf Spectral Binary Candidates From the Backyard Worlds: Planet 9 Citizen Science Initiative

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.


INTRODUCTION
Brown dwarfs occupy a unique space between stars and planets, possessing masses below the threshold required for sustained hydrogen fusion in their cores (Kumar 1962(Kumar , 1963;;Hayashi & Nakano 1963).Their cool temperatures and intrinsic faintness make them difficult to detect, and nearby brown dwarfs are continuing to be discovered (e.g., Marocco et al. 2019;Best et al. 2020;Meisner et al. 2020a,b;Bardalez Gagliuffi et al. 2020;Kirkpatrick et al. 2021a,b;Schneider et al. 2021;Kota et al. 2022;Lodieu et al. 2022;Schapera et al. 2022).The Backyard Worlds: Planet 9 project (Kuchner et al. 2017) leverages the collaboration between citizen and professional scientists to locate nearby substellar objects by identifying and analyzing moving objects in images from the Wide-field Infrared Survey Explorer (Wright et al. 2010;Mainzer et al. 2014).Multi-epoch WISE images and motion measurements from WISE data (Cat-WISE 2020; Marocco et al. 2021) have allowed for numerous discoveries missed by previous surveys, many of which were based solely on infrared colors or lim-ited to nearby objects with very large proper motions.In addition to finding new nearby brown dwarfs, the Backyard Worlds project has been adept at discovering unusual substellar objects, such as old, low-metallicity subdwarfs (Schneider et al. 2020;Meisner et al. 2021;Brooks et al. 2022) and co-moving companions (Faherty et al. 2020;Rothermich et al. 2021;Jalowiczor et al. 2021;Faherty et al. 2021;Kiwy et al. 2021;Softich et al. 2022;Kiwy et al. 2022;Gramaize et al. 2022).In this study, we examine three spectrally peculiar brown dwarfs found in the Backyard Worlds program: CWISE J072708.09−360729.2 (W0727−3607), CWISE J103604.84−514424.4 (W1036−5144), and CWISE J134446.62−732053.9 (W1344−7320).These objects were discovered by citizen scientists Dan Caselden, Arttu Sainio, and Les Hamlet.Each were flagged as high-priority for follow-up observations because their photometry and estimated spectral types indicated that they may be nearby (d < 25 pc).Properties of these three sources are provided in Table 1. a All near-infrared photometry comes from the Vista Hemisphere Survey (VHS; McMahon et al. 2013), with the exception of the J-band magnitude of CWISE J134446.62−732053.9, which comes from 2MASS (Skrutskie et al. 2006).All proper motions listed are taken from the CatWISE 2020 catalog (Marocco et al. 2021).

ANALYSIS
To determine preliminary spectral types, we compared to the near-infrared L and T dwarf standards from Kirkpatrick et al. (2010) and Burgasser et al. (2006b), with the exception of the L7 standard, where we use 2MASSI J0825196+211552 as recommended in Cruz et al. (2018), as the original L7 near-infrared standard appears to have an unusually low surface gravity and age.The best fitting standards for each object are shown in the left column of Figure 1, and match reasonably well in the Jband region, but show significant discrepancies beyond ∼1.5 µm.To ensure that differences between standards and the observed spectra were not due to uncertainties incurred when stitching orders or some other observational systematic, we calculated synthetic photometric J − K S colors using the VISTA filter response curves.We find synthetic colors consistent with the photometric colors in Table 1 to within 1σ.For this reason, we constructed binary templates in an attempt to better fit the observed spectra.
To generate spectral binary templates, we acquired near-infrared spectra from the the SpeX Prism Library Analysis Toolkit (SPLAT; Burgasser & Splat Development Team 2017).We absolute flux calibrated each spectrum using UKIDSS Hemisphere Survey (UHS; Dye et al. 2018) or Vista Hemisphere Survey (VHS; McMahon et al. 2013) J-band photometry, and measured parallaxes (Tinney et al. 2003;Dupuy & Liu 2012;Manjavacas et al. 2013;Dahn et al. 2017;Best et al. 2020;Kirkpatrick et al. 2021a;Gaia Collaboration et al. 2023).Finally, we added the spectra together and normalized the result to the J-band peak between 1.27-1.29 µm.We compared the resulting templates to each observed a This object replaces the L2 near-infrared standard (Kelu-1), which is a resolved binary (Liu & Leggett 2005).b This object replaces the L5 near-infrared standard (2MASS J08350622+1953050), which has no measured parallax.c This object replaces the T0 near-infrared standard (2MASS J12074717+0244249), which is a suspected spectral binary (Burgasser et al. 2010).d This object replaces the T3 near-infrared standard (2MASS J12095613−1004008), which is a resolved binary (Liu et al. 2010) spectrum and found the best fits by calculating χ 2 ν values following Burgasser et al. (2010).

Spectral Index Calculations
Spectral indices are also used to identify potential objects of interest (e.g., young brown dwarfs, spectral binaries, photometric variables).We follow the methods and calculate the indices defined in Burgasser et al. (2002Burgasser et al. ( , 2006bBurgasser et al. ( , 2010) ) and Bardalez Gagliuffi et al. (2014) for our observed objects (Table 3) and a large comparison sample from the SPLAT library.Uncertainties were determined in a Monte Carlo fashion.We use binary index regions from Burgasser et al. (2010) 'and variability regions from Ashraf et al. (2022) in the following sections.

RESULTS AND DISCUSSION
The observed discrepancies between observed spectra and existing spectral standards can be attributed to various factors, such as binarity, variability, youth, and metallicity, or a combination of these influences, as discussed below.

Analysis of Binary Templates
The right column of Figure 1 compares the best-fitting binary templates to our sources.The discrepancies observed with the single fits beyond approximately 1.5 µm are greatly reduced with the binary templates, as verified by lower χ 2 ν values.A discussion of the results for each individual object follows.
CWISE J072708.09−360729.2:The best-fitting T4 single standard exhibits a significantly higher χ 2 ν = 503 compared to the L7+T4 binary template, which has χ 2 ν = 39.7.Following Burgasser et al. (2010), we evaluate our fits with a one-sided F -test, finding η SB = 12.7, well above the η SB >1.34 spectral binary confidence threshold given in that work.Although the binary template provides a superior fit overall, discrepancies persist, particularly at the H-band peak.We retain W0727−3607 as a spectral binary candidate, but discuss below that it is also a potential variable source (see Section 5).
CWISE J103604.84−514424.4:W1036−5144 shows a lower-than-expected peak in the H-band and more flux across the K-band compared to the T2 standard.These discrepancies are largely resolved in the L7+T4 binary fit, with only minor deviations in the J-band and in the longest wavelengths of the K-band.The binary fit yields a significantly lower χ 2 ν = 16.8, compared to the χ 2 ν = 107.4for the single fit, and yields η SB = 6.4,consistent with the spectral binary criteria specified in Burgasser et al. (2010).CWISE J134446.62−732053.9:The single standard fit to W1344−7320 has H-and K-bands that are underluminous, leading to a relatively poor χ 2 ν = 610.The L7+T7 binary template yields a significantly improved fit with χ 2 ν = 13.7.Although there is a slight discrepancy in Y -band, where the template is too bright compared to the observed spectrum, the remainder of the J-, H-, and K-band regions exhibit a well-matched morphology.Our χ 2 ν values correspond to η SB = 44.5, which again satisfies the spectral binary criteria specified in Burgasser et al. (2010).

Analysis of Binary Indices
The spectral indices shown to be indicative of binarity in Burgasser et al. (2010) are shown in Figure 2.That study designated sources that satisfied at least three index criteria as "strong" candidates, and two index criteria as "weak" candidates.
CWISE J072708.09−360729.2:W0727−3607 falls within two of five regions indicative of binarity, and are very close to the edges of these regions.This classifies the source as a weak candidate.Note that we exclude for consideration the H 2 O−J/H 2 O−H versus spectral type comparison as it is not applicable for spectral types later than T3.5.Along with the mediocre fit of the L7+T4 binary template, this result suggests a factor other than binarity is responsible for this object's unusual spectrum.
CWISE J103604.84−514424.4:W1036−5144 satisfies five of six regions, making it a strong binary candidate.Considering the well-fit binary template, binarity is a strong possibility for this source.
CWISE J134446.62−732053.9:W1344−7320 satisfies five of five binary regions (excluding H 2 O−J/H 2 O−H versus spectral type), making it a strong binary candidate.Again, the well-fit binary template of L7+T7 makes this object a likely spectral binary.

Distance Estimation
If these objects are indeed spectral binaries, then updated distance estimates can be obtained by scaling the best-fitting absolutely flux calibrated binary templates to each object's observed photometric magnitudes.Using J-band photometry from VHS for W0727−3607 and W1036−5144 and 2MASS (Skrutskie et al. 2006) for W1344−7320, we find estimated distances of 34±7 pc, 27±5 pc, and 20±4 pc, respectively.

Variability
There is evidence that some spectral binary candidates are instead single stars with inhomogeneous atmospheres as indicated via their photometric or spectroscopic variability (e.g., 2MASS J21392676+0220226;  Burgasser et al. (2010).
c Based on the index criteria defined in Ashraf et al. (2022).2002) Radigan et al. 2012;Khandrika et al. 2013).We applied the spectral index criteria introduced in Ashraf et al. (2022) to analyze the spectra of W0727−3607, W1036−5144, and W1344−7320 and evaluate their potential as strong variable sources.This method uses single low-resolution spectra to empirically identify spectral indices that may be indicative of variability based on known variable sources.The idea is that cloudy or patchy layers at different temperatures in a brown dwarf atmosphere that lead to photometric variability may have measurable effects in the emergent spectra of these objects.Note that we do not examine any indexindex criteria that use the H 2 O−K index, as the region of the numerator range (1.975−1.995µm) is especially noisy in our TripleSpec4.1 spectra.In Figure 3, there are nine total index-index correlation plots, as discussed below.Ashraf et al. (2022) only considered objects that satisfied all of their variability index criteria as "strong" candidates, and those that satisfied all but one of their outlined criteria as "weak" candidates.

Analysis of Variable Indices
CWISE J072708.09−360729.2:W0727−3607 falls within all of the designated regions indicating potential variability.Coupled with the weak spectral binary designation and less accurate match to binary templates, variability appear to be at least partly responsible for the unusual shape of this object's spectrum.We investigated potential variability using WISE single-exposure photometry following the methods in Schneider et al. (2023).We found no clear signs of variability in the single-exposure WISE data, with W1 and W2 photometric standard deviation values within 1σ of the median value for objects with similar magnitudes.However, we note that with magnitudes of W1 = 14.639±0.018mag and W2 = 13.862±0.15mag, W0727−3607 may be too faint to detect significant variability in WISE data, as no variables were identified in Schneider et al. (2023) with magnitudes as faint as these.
However, our investigation of variability indices suggests that CWISE J072708.09−360729.2 is a strong variability candidate.
The rarity of brown dwarf-brown dwarf binary systems, as consistently highlighted by statistical studies (Burgasser 2007;Radigan et al. 2013;Aberasturi et al. 2014;Opitz et al. 2016;Fontanive et al. 2018), underscores the importance of pinpointing more such systems.Moreover, the ability to identify spectral binaries offers the advantage that their orbital period are sufficiently short to be amenable to dynamical mass follow-up (Burgasser et al. 2012;Bardalez Gagliuffi et al. 2015;Burgasser et al. 2016;Sahlmann et al. 2020).Similarly, variable brown dwarfs provide opportunities to study their rotational and atmospheric dynamics, particularly through multi-wavelength studies (Buenzli et al. 2012;Apai et al. 2013;Radigan et al. 2013).Future work includes confirming the binary and/or variability nature of these objects through high resolution imaging and photometry, astrometric, and radial velocity monitoring.Increasing the sample of close binary and variable brown dwarfs has the potential to advance our understanding of substellar formation, evolution, and dynamics.Software: SPLAT (Burgasser & Splat Development Team 2017), Spextool (Vacca et al. 2003;Cushing et al. 2004) Facilities: SOAR, WISE

Figure 1 .
Figure1.Observed SOAR/TripleSpec4.1 spectra for W0727−3607 (top), W1036−5144 (middle), and W1344−7320 (bottom), with gray lines showing the full resolution and black lines showing the spectra smoothed to the resolution of the standards and binary templates (R≈150) according to their parallax-based absolute fluxes.The left-hand column depicts the best-fitting single standards (orange line) for the spectra of the three candidates while the right-hand column shows both the binary template (orange line) and the standard spectra that compose it (purple line for primary, blue line for secondary).Reduced χ 2 ν values are reported for each of the fits.and CWISE J134446.62−732053.9.By constructing binary templates and comparing them to the observed spectra, we found significantly better fits, revealing component types of T0+T7 for CWISE J072708.09−360729.2,L7+T4 for CWISE J103604.84−514424.4,and L7+T7 for CWISE J134446.62−732053.9.However, our investigation of variability indices suggests that CWISE J072708.09−360729.2 is a strong variability candidate.The rarity of brown dwarf-brown dwarf binary systems, as consistently highlighted by statistical studies(Burgasser 2007;Radigan et al. 2013;Aberasturi et al. 2014;Opitz et al. 2016;Fontanive et al. 2018), underscores the importance of pinpointing more such systems.Moreover, the ability to identify spectral binaries offers the advantage that their orbital period are sufficiently short to be amenable to dynamical mass follow-up(Burgasser et al. 2012;Bardalez Gagliuffi et al. 2015;Burgasser et al. 2016;Sahlmann et al. 2020).Similarly, variable brown dwarfs provide opportunities to study their rotational and atmospheric dynamics, particularly through multi-wavelength studies(Buenzli et al. 2012;Apai et al. 2013;Radigan et al. 2013).Future work includes confirming the binary and/or variability nature

Figure 3 .
Figure 3. Index-index plots highlighting the regions indicative of photometric variability defined in Ashraf et al. (2022) (green shaded regions).Background grey points are SPLAT objects with spectral types between L7 and T3, with our three observed targets indicated shown in the legend.

Table 1 .
Candidate Spectral Binary Properties a

Table 2 .
Properties of Objects Used to Create Binary Templates

Table 2 .
Spectral Index Values Based on the index criteria defined in Bardalez Gagliuffi et al. (2015)egions defined inBurgasser et al. (2010);Bardalez Gagliuffi et al. (2015)as indicative of spectral binarity (red shaded regions).Background grey points are SPLAT objects with spectral types between L3 and T8.The three targets of this study are shown in the legend.This publication makes use of data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration.This publication also makes use of data products from NEO-WISE, which is a project of the Jet Propulsion Laboratory/California Institute of Technology, funded by the Planetary Science Division of the National Aeronautics and Space Administration.