The Young L Dwarf 2MASS J11193254-1137466 is a Planetary-Mass Binary

TL;DR

This paper reports the discovery of a planetary-mass binary brown dwarf system, 2MASS J11193254-1137466, which is the lowest-mass binary found to date and provides crucial data for testing models of young planetary-mass objects.

Contribution

The discovery of the lowest-mass binary brown dwarf system at 3.7 Jupiter masses, with detailed analysis of its properties and implications for planetary-mass object models.

Findings

2MASS J11193254-1137466 is a binary with 3.7 Jupiter masses each.

The system is likely a member of the TW Hydrae Association.

It exhibits a probable flux-reversal binary configuration.

Abstract

We have discovered that the extremely red, low-gravity L7 dwarf 2MASS J11193254-1137466 is a 0.14" (3.6 AU) binary using Keck laser guide star adaptive optics imaging. 2MASS J11193254-1137466 has previously been identified as a likely member of the TW Hydrae Association (TWA). Using our updated photometric distance and proper motion, a kinematic analysis based on the BANYAN II model gives an 82% probability of TWA membership. At TWA's 10$\pm$3 Myr age and using hot-start evolutionary models, 2MASS J11193254-1137466AB is a pair of $3.7^{+1.2}_{-0.9}$ $M_{\rm Jup}$ brown dwarfs, making it the lowest-mass binary discovered to date. We estimate an orbital period of $90^{+80}_{-50}$ years. One component is marginally brighter in $K$ band but fainter in $J$ band, making this a probable flux-reversal binary, the first discovered with such a young age. We also imaged the spectrally similar TWA L7 dwarf WISEA J114724.10-204021.3 with Keck and found no sign of binarity. Our evolutionary model-derived $T_{\rm eff}$ estimate for WISEA J114724.10-204021.3 is $\approx$230 K higher than for 2MASS J11193254-1137466AB, at odds with their spectral similarity. This discrepancy suggests that WISEA J114724.10-204021.3 may actually be a tight binary with masses and temperatures very similar to 2MASS J11193254-1137466AB, or further supporting the idea that near-infrared spectra of young ultracool dwarfs are shaped by factors other than temperature and gravity. 2MASS J11193254-1137466AB will be an essential benchmark for testing evolutionary and atmospheric models in the young planetary-mass regime.

Tables

Table 1. Table 1: Keck LGS AO Observations

Object	Date	Filter	$N$	$t_{\mathrm{int}}$	Airmass	FWHMaaFWHM and Strehl ratios were calculated from each image’s fitted PSF, except for the $Y$-band data for which we isolated the brighter object by rotating the image 180∘ about the fainter object and subtracting. The tabulated uncertainties are the rms of measurements from individual images.	Strehl ratioaaFWHM and Strehl ratios were calculated from each image’s fitted PSF, except for the $Y$-band data for which we isolated the brighter object by rotating the image 180∘ about the fainter object and subtracting. The tabulated uncertainties are the rms of measurements from individual images.	Separation	Position angle	$\Delta$mag
	(UT)			(sec)		(mas)		(mas)	(deg)
2MASS J11193254$-$1137466	2016 Nov 25	$K_{\mathrm{MKO}}$	3	60	1.26	$94\pm 3$	$0.148\pm 0.015$	$137.8\pm 1.7$	$239.2\pm 0.5$	$0.125\pm 0.010$
	2017 Mar 18	$K_{\mathrm{MKO}}$	6	60	1.17	$59\pm 2$	$0.41\pm 0.04$	$138.10\pm 0.10$	$239.08\pm 0.07$	$0.027\pm 0.010$
		$H_{\mathrm{MKO}}$	6	60	1.17	$51.3\pm 0.5$	$0.212\pm 0.011$	$138.08\pm 0.23$	$238.95\pm 0.07$	$0.008\pm 0.007$
		$J_{\mathrm{MKO}}$	6	120	1.17	$62\pm 9$	$0.062\pm 0.010$	$137.54\pm 0.27$	$238.72\pm 0.17$	$-0.097\pm 0.004$+
		$Y_{\mathrm{NIRC2}}$	5	120	1.18	$70\pm 27$	$0.036\pm 0.008$	$136.7\pm 3.5$	$238.4\pm 0.7$	$-0.094\pm 0.063$+
WISEA J114724.10$-$204021.3	2016 May 03	$K_{\mathrm{MKO}}$	4	60	1.33	$108\pm 11$	$0.088\pm 0.025$	$\cdots$	$\cdots$	$\cdots$

Table 2. Table 2: Properties of 2MASS J1119 − - 1137AB

Property	A component	B component	Ref.
	(northeast)	(southwest)
Observed
R.A.aaEpoch 54858.45 (MJD). (deg)	169.88521		1
Decl.aaEpoch 54858.45 (MJD). (deg)	$-11.62990$		1
${\mu }_{\alpha }\cos \delta$ (mas yr-1)	$-154.0\pm 4.0$		1
${\mu }_{\delta }$ (mas yr-1)	$-107.9\pm 1.8$		1
Radial velocity (km s-1)	$8.5\pm 3.3$		2
SpT	L7 vl-g		3,4,5
$Y$ (mag)	$19.045\pm 0.093$		2
$J$ (mag)	$17.330\pm 0.029$		2
$H$ (mag)	$15.884\pm 0.017$		2
$K$ (mag)	$14.658\pm 0.066$bbSynthetic photometry based on the SpeX prism spectrum (K15) and VHS $K_S$ photometry (K16).		3
$Y$ (mag)	$19.84\pm 0.10$	$19.75\pm 0.10$	3
$J$ (mag)	$18.13\pm 0.03$	$18.04\pm 0.03$	3
$H$ (mag)	$16.59\pm 0.02$	$16.60\pm 0.02$	3
$K$ (mag)	$15.40\pm 0.07$	$15.43\pm 0.07$	3
$Y-J$ (mag)	$1.71\pm 0.10$	$1.71\pm 0.10$	3
$J-H$ (mag)	$1.54\pm 0.03$	$1.43\pm 0.03$	3
$J-K$ (mag)	$2.73\pm 0.07$	$2.61\pm 0.07$	3
$\Delta \left(J-K\right)$ (mag)	$0.125\pm 0.011$		3
Estimated
$d_{\mathrm{phot}}$ (pc)	$26.4\pm 6.9$		3,8
$m-M$ (mag)	$2.11\pm 0.56$		3
$v_{\tan }$ (km s-1)	$23.6\pm 6.1$		3
Projected separation (AU)	$3.6\pm 0.9$		3
Semi-major axis (AU)	${3.9}_{-1.4}^{+1.9}$		3
Orbital Period (yr)	${90}_{-50}^{+80}$		3
log ($L_{\mathrm{bol}}/L_{\odot }$) (unresolved)	$-{4.44}_{-0.27}^{+0.21}$		3,6
log ($L_{\mathrm{bol}}/L_{\odot }$) (resolved)	$-{4.73}_{-0.27}^{+0.21}$	$-{4.74}_{-0.27}^{+0.21}$	3
Model-derived (Lyon/DUSTY) assuming TWA membership
Age (Myr)	$10\pm 3$ Myr		7
Mass ($M_{\mathrm{Jup}}$)	${3.7}_{-0.9}^{+1.2}$	${3.7}_{-0.9}^{+1.2}$	3
$T_{\mathrm{eff}}$ (K)	${1013}_{-109}^{+122}$	${1006}_{-109}^{+122}$	3
Model-derived (Lyon/DUSTY) assuming young field (vl-g)
Age (Myr)	$10-100$ Myr		8
Mass ($M_{\mathrm{Jup}}$)	${9.2}_{-1.9}^{+2.3}$	${9.0}_{-1.9}^{+2.4}$	3
$T_{\mathrm{eff}}$ (K)	${1065}_{-118}^{+133}$	${1059}_{-118}^{+133}$	3

Full text

brown dwarfs — binaries: close — stars: individual (2MASS J11193254$-$1137466, WISEA J114724.10$-$204021.3)

††slugcomment: Accepted by ApJ Letters, 2017 June 01

1 Introduction

Brown dwarfs with masses $\lesssim$15 $M_{\mathrm{Jup}}$ and ages $\lesssim$100 Myr lie at a nexus of astronomical interest. They represent the lowest-mass and youngest products of star formation, and as such offer rare empirical tests for evolutionary and atmospheric models. They are also the best field analogs to directly-imaged giant exoplanets, which are far more difficult to directly observe due to the glare of their host stars.

Brown dwarfs cool continuously as they age, and the resulting mass-age-luminosity degeneracy makes their physical properties challenging to infer without constraints on at least two of those three parameters. The atmospheres of young brown dwarfs exhibit clear spectral signatures of low gravity (Cruz et al., 2009; Allers & Liu, 2013), but the age calibration for these signatures lacks precision better than $\approx$100 Myr (e.g., Liu et al., 2016). Some of the lowest-mass objects have been identified as members of nearby young moving groups (e.g., Gagné et al., 2014), which provide much tighter age constraints and thus more precise mass estimates from evolutionary models than for ordinary field objects. Young binaries with small separations are even more useful as benchmarks, as their orbits can yield model-independent dynamical masses, providing exacting tests for models (e.g., Dupuy et al., 2009; Dupuy & Liu, 2017).

2MASS J11193254$-$1137466 (a.k.a. TWA 42; hereinafter 2MASS J1119$-$1137) was discovered by Kellogg et al. (2015, hereinafter K15) in a search for L and T dwarfs with unusual photometry. 2MASS J1119$-$1137 is an L7 dwarf with extremely red optical and near-IR colors along with spectral signatures of low gravity indicating youth (K15; Kellogg et al., 2016, hereinafter K16). K16 identified 2MASS J1119$-$1137 as a candidate member of the TW Hydrae Association (TWA; Webb et al., 1999), whose age implies a mass of only 4.3–7.6 $M_{\mathrm{Jup}}$ for this object. 2MASS J1119$-$1137 would be one of the two lowest-mass isolated members of TWA, comparable only to the L7 dwarf WISEA J114724.10$-$204021.3 (hereinafter WISEA J1147$-$2040; $6.6\pm 1.9$ $M_{\mathrm{Jup}}$; Schneider et al., 2016; Faherty et al., 2016), and among the lowest-mass free-floating brown dwarfs known.

We are conducting a high angular-resolution imaging survey of nearby brown dwarfs to identify binaries. In this Letter we show that 2MASS J1119$-$1137 is a nearly equal-flux binary with component masses in the planetary regime.

2 Observations

We observed 2MASS J1119$-$1137 on 2016 November 25 UT using the laser guide star adaptive optics (LGS AO) system at the Keck II Telescope (van Dam et al., 2006; Wizinowich et al., 2006). We used the facility infrared camera NIRC2 in its narrow field-of-view configuration, using the $R=13.5$ mag field star USNO-B1.0 0783-0249513 (Monet et al., 2003) located 35″ from 2MASS J1119$-$1137 for tip-tilt correction. Skies were mostly clear, with $K$-band seeing of $1.6^{\prime\prime}$ measured contemporaneously at UKIRT. We obtained 3 dithered images at $K$ band in which 2MASS J1119$-$1137 appeared to be an equal-flux binary. We observed 2MASS J1119$-$1137 again on 2017 March 18 UT in $YJHK$ bands using the same configuration, under clear skies with seeing $\approx 0.5^{\prime\prime}$ as measured by the differential image motion monitor (DIMM) at the Canada-France-Hawaii Telescope (CFHT). Details are in Table 1.

In addition, we observed WISEA J1147$-$2040 on 2016 May 18 UT with a similar configuration using the $R=17.4$ mag field star USNO-B1.0 0693-0264226 (Monet et al., 2003) located 66″ away for tip-tilt correction, under clear skies with seeing $\approx 0.5^{\prime\prime}$ from DIMM. The target appeared to be a single object at $0.11^{\prime\prime}$ resolution.

We reduced and analyzed our data using the methods described in, e.g., Liu et al. (2008) and Dupuy & Liu (2017). Briefly, we calibrated our images using flat fields and dark frames, performed sky subtraction, and registered and stacked images to form a final mosaic for each epoch and filter (Figure 1). We measured the relative astrometry and flux ratios of the binary by fitting a PSF model to the individual images, applying the NIRC2 pixel scale, orientation, and distortion correction from Service et al. (2016). For images in which the components were well separated (all but $Y$ band), we used the StarFinder software package (Diolaiti et al., 2000) to simultaneously solve for an empirical PSF and binary parameters. For $Y$ band, we used an analytical PSF of two elliptical three-component Gaussians. We applied additional corrections for differential aberration and atmospheric refraction. We used the rms of the measurements from individual images as the uncertainties on the separation, position angle (PA), and $\Delta$mag of the binary components (Table 1), adding the errors in plate scale (0.4%) and orientation (0.020 deg) from Service et al. (2016) in quadrature.

The NIRC2 $J$, $H$, and $K$ filters we used are from the Mauna Kea Observatories (MKO) photometric system (Simons & Tokunaga, 2002; Tokunaga et al., 2002), and the NIRC2 $Y$-band filter is described in Liu et al. (2012). The unresolved photometry reported in K16 is from the VISTA Hemisphere Survey (VHS; PI McMahon, Cambridge, UK), which uses MKO $J$ and $H$ filters but a non-MKO $K_{\rm S}$ filter. We used the IRTF/SpeX spectrum from K15 to calculate a synthetic $K_{\rm MKO}=14.658\pm 0.066$ mag for 2MASS J1119$-$1137, flux-calibrated with K16’s $K_{\rm S}$ magnitude. The VISTA and NIRC2 $Y$ filters are similar enough that no conversion was necessary.

We split the unresolved 2MASS J1119$-$1137 $YJH$ photometry from K16 and our synthetic $K$ magnitude into resolved photometry using our measured flux ratios (Table 3). $K$-band flux decreases monotonically with spectral type (e.g., Dupuy & Liu, 2012), and the northeast component of 2MASS J1119$-$1137 is slightly brighter in $K$, so we designate this object as the “A” component. We note a $0.098\pm 0.014$ mag difference between the $K$-band flux ratios measured at the two epochs. This may indicate variability in one or both components, or systematic errors unaccounted in our uncertainties. We use the 2017 March $K$-band flux ratio for our analysis, as photometry in all other bands was measured that night, and because of the better image quality.

We also identified a faint source in our $JHK$ images from March 2017, lying $3.79\pm 0.02^{\prime\prime}$ from 2MASS J1119$-$1137A at ${\rm PA}=76.47\pm 0.13$ deg. It was not well detected but appears to be a point source. Images from the DSS, 2MASS, SDSS, AllWISE, and Pan-STARRS1 surveys indicate no object at this location. We measure flux ratios of 4.8 mag in $J$, 5.0 mag in $H$, and 5.7 mag in $K$, relative to 2MASS J1119$-$1137A. At the same distance as 2MASS J1119$-$1137, this source’s $J=22.9$ mag would be consistent with known Y0 dwarfs, but its $(J-K)_{\rm MKO}=1.8$ mag color is $\gtrsim$3 mag too red (Leggett et al., 2017). It is almost certainly a background object. Its $(J-K)_{\rm MKO}$ color suggests an L dwarf, evolved star, or galaxy.

3 Results

\floattable

3.1 2MASS J1119$-$1137AB is comoving

Our $JHK$-band astrometry from March 2017 has a mean separation of $137.88\pm 0.34$ mas and PA $238.91\pm 0.20$ deg, with uncertainties estimated as in Section 2. The change in separation from November 2016 is $0.7\pm 1.5$ mas, consistent with no change. Using the proper motion and photometric distance of 2MASS J1119$-$1137 (Section 3.3), if 2MASS J1119$-$1137B were a stationary background object the separation would have decreased by $91.7\pm 9.8$ mas in March 2017, inconsistent by $9.4\sigma$ from our observation (Figure 2). In addition, images from DSS, 2MASS, SDSS, and Pan-STARRS1 showed no objects that could appear as a false close companion given the proper motion. We therefore conclude that 2MASS J1119$-$1137AB is a gravitationally bound binary.

3.2 Spectral type and gravity classification

We used the IRTF/SpeX prism spectrum for 2MASS J1119$-$1137 and the method of Allers & Liu (2013) to determine a spectral type of L7, concurring with previous work (K15, Faherty et al., 2016; Gagné et al., 2017). A higher-resolution $J$-band spectrum (K16) shows weakened K I absorption lines having equivalent widths consistent with vl-g classification for L7 dwarfs (Allers & Liu, 2013; Gagné et al., 2017). We therefore adopt L7 vl-g as the unresolved spectral type. The extremely red $(J-K)_{\rm MKO}$ colors of 2MASS J1119$-$1137AB are consistent with other low-gravity late-L dwarfs (Figure 3). Both the $(J-K)_{\rm MKO}$ colors and $K_{\rm MKO}$ magnitudes for the two components are similar, and the integrated-light spectrum shows no peculiarities that would suggest a blend of two objects with different spectral types. We conclude that both 2MASS J1119$-$1137A and 2MASS J1119$-$1137B have spectral types L7 vl-g.

3.3 Distance

We used the spectral-type-to-$M_{K_{\rm MKO}}$ relation for vl-g dwarfs from Liu et al. (2016) and our resolved photometry to estimate photometric distances of $26.3\pm 6.8$ pc and $26.6\pm 6.9$ pc for 2MASS J1119$-$1137A and 2MASS J1119$-$1137B, respectively. We adopt $26.4\pm 6.9$ pc as the distance to the system, giving the binary a projected separation of $3.6\pm 0.9$ AU. Our distance is consistent with K16’s kinematic distance of $28.9\pm 3.6$ pc assuming membership in TWA.

3.4 TWA Membership

K16 identified 2MASS J1119$-$1137 as a candidate TWA member, using the BANYAN II online tool (Malo et al., 2013; Gagné et al., 2014) to calculate an 88% membership probability with a contamination probability of 0.003% based on 2MASS J1119$-$1137’s sky position, proper motion, radial velocity, and youth. Faherty et al. (2016) analyzed moving group membership using four different tools and found $>$90% probabilities for TWA in three cases; LACEwING (Riedel et al., 2017) found a contrasting probability of 16%. Gagné et al. (2017) used the full BANYAN II analysis (including photometry) to calculate a 97% probability of membership. Liu et al. (2016) found that objects lacking parallaxes which have full BANYAN II probabilities $\gtrsim$80% tend to have memberships confirmed by subsequent parallax measurements for the well-established moving groups, including TWA.

We reassessed the moving group membership of 2MASS J1119$-$1137 using the BANYAN II online tool. For position and proper motion, we adopted the values from Best et al. (2017), which are calculated from 2MASS (Skrutskie et al., 2006) and Pan-STARRS1 $3\pi$ (Chambers et al., 2017; Magnier et al., 2017) astrometry and calibrated to the Gaia DR1 reference frame (Gaia Collaboration et al., 2016; Lindegren et al., 2016). Our $\mu_{\delta}$ differs by 35.5 mas yr*-1* ($2\sigma$) from K16’s value but is consistent with measurements by K15 and Gagné et al. (2017).

Using only our astrometry and the radial velocity from K16 as inputs to BANYAN II, we obtained a TWA membership probability of 73%. The decrease from K16’s 88% probability is due almost entirely to the difference in $\mu_{\delta}$. When we included our photometric distance, the probability of TWA membership increased to 82%.

Based on our BANYAN II results and the vl-g spectrum, 2MASS J1119$-$1137 is a very likely member of TWA. We note that a comparison in UVW and XYZ of 2MASS J1119$-$1137 to known members of moving groups using $\tilde{\chi}^{2}$ as a rubric (see Equation 1 of Liu et al., 2016) supports membership in TWA, $\beta$ Pictoris (0% probability from BANYAN II), and Tucana-Horologium (0%), and especially the AB Doradus (4%) moving group. A precise trigonometric distance, now underway by us at CFHT, is needed to firmly establish the membership.

3.5 Physical Properties

As a TWA member, 2MASS J1119$-$1137 would share the age of $10\pm 3$ Myr derived from the stellar members (Bell et al., 2015). We estimated the components’ masses and effective temperatures using bolometric luminosities and the Lyon/DUSTY hot-start evolutionary models (Chabrier et al., 2000). To obtain $L_{\rm bol}$, we adjusted the integrated-light value from Faherty et al. (2016) to our new photometric distance, and decomposed this into individual $L_{\rm bol}$ values using the binary’s $K_{\rm MKO}$ flux ratio. Using these $L_{\rm bol}$ and the assumed age, we then interpolated a mass and effective temperature from the DUSTY models for each component. We propagated the uncertainties on distance, flux ratio, and age into our calculations via Monte Carlo trials using normal distributions for each uncertainty, and we quote the resulting median and 68% confidence limits (Table 3). We estimate masses of $3.7^{+1.2}_{-0.9}$ $M_{\mathrm{Jup}}$ for both components. We did not use the more recent BHAC15 models (Baraffe et al., 2015) as they do not include masses below 0.01 $M_{\odot}$.

If 2MASS J1119$-$1137 is not a TWA member, its vl-g classification still constrains its age to $\approx$10–100 Myr (Liu et al., 2016). Our estimated masses and effective temperatures for this age range (uniformly distributed for error propagation) are shown in Table 3.

Using our masses, we considered the effect of 2MASS J1119$-$1137AB’s orbital motion on the radial velocity measured by K16. For an edge-on circular orbit with components at quadrature, we found a difference in radial velocities of $1.3^{+0.4}_{-0.3}$ km s*-1* for TWA masses or $2.0^{+0.6}_{-0.4}$ km s*-1* for field-vl-g masses. The $R\sim 6000$ (50 km s*-1*) spectrum of K16 would not resolve such orbital motion, so the radial velocity from K16 remains valid.

4 Discussion

Assuming 2MASS J1119$-$1137 is a member of TWA, its total mass is $7.4^{+2.5}_{-1.9}$ $M_{\mathrm{Jup}}$, making it the lowest-mass binary discovered to date. The individual masses of 2MASS J1119$-$1137A and 2MASS J1119$-$1137B also place them among the lowest-mass free-floating brown dwarfs, including WISEA J1147$-$2040 and the Y dwarfs (Dupuy & Kraus, 2013; Leggett et al., 2017). Even if 2MASS J1119$-$1137AB is actually a young field object, its total mass of $18.2^{+4.7}_{-3.8}$ $M_{\mathrm{Jup}}$ would still be among the lowest-mass binaries, surpassed only by the young binary DENIS-P J035726.9-441730 (14–15 $M_{\mathrm{Jup}}$; Bouy et al., 2003; Gagné et al., 2014) and possibly by the T9+Y0 binaries CFBDSIR J145829+101343 ($\approx$18–45 $M_{\mathrm{Jup}}$; Liu et al., 2011) and WISE J014656.66+423410.0 (Dupuy et al., 2015a). As an extremely low-mass young binary, 2MASS J1119$-$1137AB will be a crucial benchmark for tests of evolutionary and atmospheric models.

The isolation of 2MASS J1119$-$1137AB strongly suggests that it is a product of normal star-formation processes, which therefore must be capable of making binaries with $\lesssim$5 $M_{\mathrm{Jup}}$ components. 2MASS J1119$-$1137AB could be a fragment of a higher-order system that was ejected via dynamical interactions (Reipurth & Mikkola, 2015), although the lack of any confirmed member of TWA within 10∘ (projected separation $\approx 5$ pc) of 2MASS J1119$-$1137 makes this scenario unlikely. Formation of very low mass binaries in extended massive disks around Sun-like stars followed by ejection into the field has been proposed by, e.g., Stamatellos & Whitworth (2009), but disks of this type have not been observed.

Binary brown dwarfs can be monitored to map their orbits, which yield dynamical masses that stringently test evolutionary models. We estimated 2MASS J1119$-$1137AB’s orbital period using Kepler’s Third Law. We first used the projected separation and a conversion factor from Dupuy & Liu (2011, Table 6, assuming moderate discovery bias for very low-mass visual binaries) to estimate a semi-major axis of $3.9^{+1.9}_{-1.4}$ AU. Our model-derived masses assuming TWA membership give an orbital period of $90^{+80}_{-50}$ yr. A dynamical mass can be determined once $\approx$1/3 of an orbit has been observed (e.g., Dupuy & Liu, 2017), so 2MASS J1119$-$1137AB would yield a dynamical mass in $\approx$15–55 years. Assuming a field vl-g age, we estimate an orbital period of $60^{+50}_{-30}$ yr, yielding a mass in $\approx$10–35 years.

The integrated-light spectrum and photometry of 2MASS J1119$-$1137 are notably similar to those of the young L7 dwarf WISEA J1147$-$2040 (Gagné et al., 2017), implying similar temperatures and gravity. However, using the $L_{\rm bol}$ for WISEA J1147$-$2040 from Faherty et al. (2016) and the method from Section 3.5, we estimate a $T_{\rm eff}$ of $1242^{+73}_{-69}$ K, which is $\approx$230 K higher than our estimates for 2MASS J1119$-$1137A and 2MASS J1119$-$1137B. This discrepancy is particularly surprising given that both objects are very likely members of TWA and therefore should have the same age and composition. The simplest resolution is that WISEA J1147$-$2040 is also an equal-flux binary, unresolved in our images, with component temperatures very similar to those of 2MASS J1119$-$1137AB. Using the method from Section 3.3 we calculated $d_{\rm phot}=27.3\pm 6.9$ pc for a hypothetical equal-flux binary WISEA J1147$-$2040, leading to a maximum projected separation of $2.9\pm 0.8$ AU to remain unresolved in our images. Another possibility is that 2MASS J1119$-$1137AB is not a member of TWA and is older; our vl-g-age temperature estimates are only $\approx$50 K higher, but would then agree with the WISEA J1147$-$2040 estimate within uncertainties. A third intriguing possibility is that the low-resolution spectra of young red L dwarfs are driven at least partially by factors other than temperature and gravity (Allers & Liu, 2013; Liu et al., 2016), allowing coeval objects with differing masses and effective temperatures to have similar spectra.

Figure 3 compares the $J_{\rm MKO}$ vs. $(J-K)_{\rm MKO}$ position of 2MASS J1119$-$1137AB to the ultracool dwarf population, highlighting other low-gravity objects and substellar companions. 2MASS J1119$-$1137AB lies among other planetary-mass objects at the faint red end of the L dwarf sequence. 2MASS J1119$-$1137B is brighter in $J$ band and slightly fainter in $K$ band, making the system a probable flux-reversal binary. (The $K$ magnitudes for the two components are formally consistent within uncertainties, but 2MASS J1119$-$1137B is fainter in all nine of our individual $K$-band images.) The flux-reversal phenomenon is a hallmark of field-age L/T transition binaries (e.g., Gizis et al., 2003; Liu et al., 2006; Dupuy et al., 2015b), thought to occur when the cooler component reaches a temperature at which the clouds that suppress near-IR flux in L dwarfs begin to clear, reducing the $J$-band opacity relative to the warmer component (e.g., Burrows et al., 2006). The slightly bluer $(J-K)_{\rm MKO}$ color of 2MASS J1119$-$1137B implies that it is cooler than its primary. In field-age brown dwarfs this transition to bluer $J-K$ colors is typically seen at warmer temperatures ($\approx$1400 K; e.g., Dupuy & Liu, 2012). The potential flux-reversal of 2MASS J1119$-$1137AB suggests that it is beginning the transition at $T_{\rm eff}\approx 1000\pm 100$ K, an even lower temperature than the $\approx$1100$-1200$ K found for other low (but somewhat higher) mass L dwarfs (e.g. Metchev & Hillenbrand, 2006; Barman et al., 2011; Liu et al., 2013), implying a possible systematic correlation between mass and L/T transition temperature. Precise photometry and resolved spectroscopy of 2MASS J1119$-$1137AB with the Hubble and James Webb Space Telescopes will enable differential studies of the atmospheres of young planetary-mass objects and may yield insights into the L/T transition at young ages.

We thank the referee for an immediate and helpful review. The IRTF/SpeX spectrum for 2MASS J1119$-$1137 was retrieved from the SpeX Prism Library, maintained by Adam Burgasser at http://www.browndwarfs.org/spexprism. The Pan-STARRS1 Surveys have been made possible through contributions of the Institute for Astronomy, the University of Hawaii, the Pan-STARRS Project Office, the Max-Planck Society and its participating institutes, the Max Planck Institute for Astronomy, Heidelberg and the Max Planck Institute for Extraterrestrial Physics, Garching, The Johns Hopkins University, Durham University, the University of Edinburgh, Queen’s University Belfast, the Harvard-Smithsonian Center for Astrophysics, the Las Cumbres Observatory Global Telescope Network Incorporated, the National Central University of Taiwan, the Space Telescope Science Institute, the National Aeronautics and Space Administration under Grant No. NNX08AR22G issued through the Planetary Science Division of the NASA Science Mission Directorate, the National Science Foundation under Grant No. AST-1238877, the University of Maryland, Eotvos Lorand University (ELTE), and the Los Alamos National Laboratory. This work has made use of data from the European Space Agency (ESA) mission Gaia (http://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, http://www.cosmos.esa.int/web/gaia/dpac/consortium). WMJB and MCL received support from NSF grant AST-1518339. WMBJ, MCL, and EAM received support from NSF grant AST-1313455. Finally, the authors acknowledge the significance of the summit of Maunakea for the Native Hawaiian community. We are fortunate to have the opportunity to conduct observations from this mountain.

Keck, Pan-STARRS1, Gaia