A publishing partnership

2MASS J06164006−6407194: THE FIRST OUTER HALO L SUBDWARF

, , , , , , , and

Published 2009 April 20 © 2009. The American Astronomical Society. All rights reserved.
, , Citation Michael C. Cushing et al 2009 ApJ 696 986 DOI 10.1088/0004-637X/696/1/986

0004-637X/696/1/986

ABSTRACT

We present the serendipitous discovery of an L subdwarf in the Two Micron All Sky Survey (2MASS) J06164006−6407194, in a search of the 2MASS for T dwarfs. Its spectrum exhibits features indicative of both a cool and metal poor atmosphere including a heavily pressure-broadened K i resonant doublet, Cs i and Rb i lines, molecular bands of CaH, TiO, CrH, FeH, and H2O, and enhanced collision induced absorption of H2. We assign 2MASS J0616−6407 a spectral type of sdL5 based on a comparison of its red optical spectrum to that of near solar-metallicity L dwarfs. Its high proper motion (μ = 1farcs405 ± 0.008 yr−1), large radial velocity (Vrad = 454 ± 15 km s−1), estimated u, v, w velocities (94, −573, 125) km s−1 and Galactic orbit with an apogalacticon at ∼29 kpc are indicative of membership in the outer halo making 2MASS J0616−6407 the first ultracool member of this population.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

Ultracool subdwarfs are metal-poor stars and brown dwarfs that have spectral types later than ∼sdM7/esdM7 (Burgasser et al. 2005) and Teff ≲ 3000 K (Leggett et al. 2000). Since cool subdwarfs are typically members of the thick disk or halo population, they exhibit large space motions relative to the Sun (〈V〉 = −202 km s−1; Gizis 1997) and thus are often identified as high proper-motion stars (e.g., Giclas et al. 1971, 1978; Luyten 1979; Pokorny et al. 2004; Lépine & Shara 2005). Due to their intrinsic faintness and paucity in the solar neighborhood, only ∼40 ultracool subdwarfs are currently known (Burgasser et al. 2007; Lépine & Scholz 2008). Such small numbers stand in stark contrast to the over one thousand near solar-metallicity ultracool dwarfs (spectral types later than M7 and Teff ≲ 2400 K) that have been discovered in the Two Micron All Sky Survey (2MASS; Skrutskie et al. 2006), the Sloan Digital Sky Survey (SDSS; York et al. 2000), and the Deep Near-Infrared Southern Sky Survey (Epchtein et al. 1997), many of which populate the L and T spectral classes6 (Kirkpatrick 2005).

Of the roughly 40 ultracool subdwarfs known, only three, 2MASS J0532346 + 8246465 (hereafter 2MASS J0532+8246; Burgasser et al. 2003), 2MASS J16262034 + 3925190 (hereafter 2MASS J1626+3925; Burgasser 2004), and SDSS J125637–022452 (hereafter SDSS J1256–0224; Sivarani et al. 2009) are classified as L subdwarfs.7 The spectra of all three L subdwarfs exhibit Cs i and Rb i lines, pressure-broadened K i lines, CaH, FeH, CrH and H2O bands, and enhanced collision induced H2 absorption (CIA, e.g., Linsky 1969; Borysow et al. 1997), all features indicative of very cool (Kirkpatrick et al. 1999; Leggett et al. 2001; Reid et al. 2001) and metal-depleted atmospheres (Bessell 1982; Saumon et al. 1994). In addition, they have large proper motions (0farcs6 to 2farcs6 yr−1) and radial velocities (−130 to −260 km s−1) consistent with halo membership. Burgasser et al. (2008) recently measured the parallax of 2MASS J0532+8246 and confirmed, based on its Galactic orbit, that it is indeed a member of the (inner) halo. Although no spectral classification scheme exists for L subdwarfs, Burgasser et al. (2007, 2009) have tentatively assigned spectral types of sdL3.5, sdL4, and sdL7 to SDSS J1256−0224, 2MASS J1626+3925, and 2MASS J0532+8246, respectively, based on a comparison with the red optical spectra of the near solar metallicity L dwarf standards (Kirkpatrick et al. 1999).

In this paper, we present the discovery of a fourth L subdwarf, 2MASS J06164006−6407194 (hereafter 2MASS J0616−6407) found in a search of the 2MASS database for T dwarfs. The discovery of 2MASS J0616−6407, along with follow-up near-infrared imaging and spectroscopy, and red-optical spectroscopy is presented in Section 2. In Section 3, we discuss the spectral classification and kinematics of 2MASS J0616−6407.

2. OBSERVATIONS

2.1. Source Identification

2MASS J0616−6407 was identified as a T dwarf candidate by Looper et al. (2007) in the course of a search for T dwarfs in the 2MASS Database. It was selected based on its very blue near-infrared color (JKs < 0.02) and lack of any optical counterpart within 5'' in the USNO-A2.0 catalog or the Digital Sky Survey (DSS) I and II images. A finder chart for 2MASS J0616−6407 is given in Figure 1 and pertinent properties are listed in Table 1. Below we describe follow-up near-infrared imaging and spectroscopy, and red-optical spectroscopy.

Figure 1.

Figure 1. J-band image of the field centered around 2MASS J0616−6407 obtained with the CPAPIR on the 1.5 m CTIO on 2008 January 17 (UT). The field is 4' on a side and north is up and east is to the left.

Standard image High-resolution image

Table 1. Properties of 2MASS J0616−6407

Parameter Value
Observed
α(J2000.0)a 6h16m40fs06
δ(J2000.0)a −64°07'19farcs4
2MASS J 16.403 ± 0.113 mag
2MASS H 16.275 ± 0.228 mag
2MASS Ks > 16.381b mag
μα (J2000.0) 1405 ± 8 mas yr−1
μδ (J2000.0) −51 ± 18 mas yr−1
Vrad 454 ± 15 km s−1
Estimated
d 57 ± 9 pc
Vtan  379 ± 60 km s−1
uc 94 ± 10 km s−1
v −573 ± 31 km s−1
w 125 ± 53 km s−1

Notes. a2MASS coordinates at epoch 1998.95 (UT). bThe read flag is 0 and the quality flag is U indicating that no flux was detected at the position and that the value is an upper limit. cWe follow the convention that the u is positive toward the Galactic center (l = 0, b = 0).

Download table as:  ASCIITypeset image

2.2. Near-Infrared Spectroscopy

2MASS J0616−6407 was observed with the Ohio State InfraRed Imager/Spectrometer (OSIRIS) mounted at the 4.1 m SOAR Telescope on 2007 March 12 (UT). We used the 1farcs0 wide slit which provides a resolving power of R ≡ λ/Δλ∼ 1400 across the 1 to 2.2 μm wavelength range in three spectral orders. A series of six 180 s integrations were obtained along the 24'' long slit. The A0 V star HD 53191 was also observed to correct for telluric absorption and to flux-calibrate the final spectrum. The data were reduced using a modified version (version 3.4) of the Spextool data reduction package (Cushing et al. 2004). Wavelength calibration was accomplished using the OH airglow lines in the science frames. The spectra from the three orders were then corrected for telluric absorption using the observations of the A0 V standard stars and the technique described in Vacca et al. (2003). The three spectral orders were then merged into a single spectrum covering the 1.2–2.3 μm wavelength range. The final spectrum was then rebinned to one pixel per resolution element in order to increase the signal-to-noise ratio (S/N) of the spectrum.

The spectrum of 2MASS J0616−6407 is shown in Figure 2 (middle). The S/N of the spectrum is low, ranging from ∼8 in the J and H bands to ∼4 in the K band, but is nevertheless adequate to confirm that 2MASS J0616−6407 is not a T dwarf (the JKs color of 2MASS J0616−6407 corresponds to a spectral type of ∼T4) since the spectrum lacks the prominent CH4 bands found in the spectra of T dwarfs (Burgasser et al. 2002; Geballe et al. 2002). Also shown in Figure 2 are the spectra of the two other L subdwarfs with published spectra, 2MASS J0532+8246 (Burgasser et al. 2003) and 2MASS J1626+3925 (Burgasser 2004). All three spectra exhibit blue continua due to enhanced collision induced absorption (CIA) of H2 centered at 2.4 μm and broad absorption bands of H2O centered at 1.4 and 1.9 μm indicating that 2MASS J0616−6407 is an L subdwarf.

Figure 2.

Figure 2. Near-infrared spectra of 2MASS J1626+3925  (sdL4; Burgasser et al. 2007), 2MASS J0616−6407, and 2MASS J0532+8246  (sdL7; Burgasser et al. 2003). The spectra are normalized to unity at 1.3 μm and offset by constants (dotted lines). Prominent atomic and molecular absorption features are indicated. The S/N of the 2MASS J0616−6407 spectrum is not high enough to detect the 1.25 μm K i doublet.

Standard image High-resolution image

2.3. Red Optical Spectroscopy

2.3.1. Magellan/LDSS-3

Red optical spectra of 2MASS J0616−6407 were obtained on 2006 May 10 (UT) with LDSS-3 on Magellan. The VPH-red grism (660 lines mm−1) with a 0farcs75 wide (4 pixels) longslit mask was used, with the slit aligned to the parallactic angle. This configuration provides 6050 Å–10500 Å spectra across the entire chip with an average resolving power of R ≈ 1800 and dispersion along the chip of ∼1.2 Å pixel−1. The OG590 longpass filter was used to eliminate second-order light shortward of 6000 Å. Two slow-read exposures of 1800 s each were obtained at an airmass of 1.49–1.60. We also observed the G2 V star HD 44611 immediately after and at a roughly similar airmass (1.77) for telluric absorption correction. The flux standard LTT 7987 (aka GJ 2147; Hamuy et al. 1994) was observed on the previous night using an identical slit and grism combination. All spectral observations were accompanied by HeNeAr arc lamp and flat-field quartz lamp exposures for dispersion and pixel response calibration.

LDSS-3 data were reduced in the IRAF8 environment (Tody 1986). Raw images were first corrected for amplifier bias voltage, stitched together, and subtracted by a median-combined set of slow-read bias frames taken during the afternoon. These processed images were then divided by a median-combined, bias-subtracted and normalized set of flat field frames. The LTT 7987 and HD 44611 spectra were optimally extracted first using the APALL task with background subtraction. The spectrum of 2MASS J0616−6407 was then extracted using the G star dispersion trace as a template. Dispersion solutions were determined from arc lamp spectra extracted using the same dispersion trace; solutions were accurate to ∼0.08 pixels, or ∼0.1 Å. Flux calibration was determined using the tasks STANDARD and SENSFUNC with observations of LTT 7987, adequate over the spectral range 6000 Å–10000 Å. Corrections to telluric O2 (6855 Å–6955 Å B-band, 7580 Å–7740 Å A-band) and H2O (7160 Å–7340 Å, 8125 Å–8350 Å, 9270 Å–9680 Å) absorption bands were determined by linearly interpolating over these features in the G dwarf spectrum, dividing by the uncorrected spectrum, and multiplying the result with the spectrum of 2MASS J0616−6407. The two spectra of 2MASS J0616−6407 were then coadded to improve S/N. We measured the radial velocity of 2MASS J0616−6407 using the Rb i (7948 Å), Na i (8183, 8195 Å), and Cs i (8521 Å) lines. The observed line positions determined by Gaussian fits were compared against the rest wavelengths (Ralchenko et al. 2008) to derive a heliocentric radial velocity (after correcting for the Earth's orbital velocity) of 480 ± 30 km s−1, where the error in the radial velocities arises primarily from the scatter in the line positions. We defer further discussion of the spectrum to the following section.

2.3.2. Gemini/GMOS Spectroscopy

2MASS J0616−6407 was also observed with the Gemini Multi-Object Spectrograph (Hook et al. 2004) mounted on the Gemini South Telescope on 2007 September 13 (UT). We used the 0farcs75 wide slit and R400 grating which yields a spectrum from 5900 to 10100 Å  with a resolving power of R∼ 1300. The OG515 filter was used to suppress stray light from shorter wavelengths. A total of five 1800 s exposures were acquired at two different central wavelengths and at two positions along the slit. Standard calibration frames including flat-field and CuAr lamp images as well as observations of the G2 V star HD 60402 for telluric correction were also obtained. The flux standard EG 21 (Hamuy et al. 1994) was observed on 2007 September 14 (UT) with the same instrument setup.

The data were first reduced using the IRAF Gemini GMOS package. The raw frames were first bias-subtracted (with a bias created using GBIAS) and flat-fielded (with a normalized flat field created with GSFLAT). The OH sky lines were then subtracted and then spectra were optimally extracted using the GSEXTRACT routine. Wavelength calibration was achieved using the CuAr lamp images. The spectra were then flux calibrated using the observations of EG 21 and the GSCALIBRATE routine. Telluric correction was achieved by linearly interpolating over regions of O2 and H2O absorption in the spectrum of HD 60402, dividing by the raw spectrum of HD 60402, and then multiplying the results into the flux-calibrated spectrum of 2MASS J0616−6407.

Unfortunately, the resulting spectrum exhibits strong noise spikes at the wavelengths of nearly all the OH sky lines. Fringing in the detector makes the subtraction of the OH sky lines using low-order polynomials difficult. We therefore re-extracted the raw spectrum of 2MASS J0616−6407 so that we could first pair subtract the images with 2MASS J0616−6407 at different positions along the slit. Although this did not completely eliminate the strong residuals, it did improve the spectrum dramatically. Finally we measured a radial velocity of 445 ±18 km s−1 using the same technique described in the previous section. The radial velocity measurements from the LDSS-3 and GMOS spectra agree within the errors and a weighted average of the two values gives 454 ± 15 km s−1.

The LDSS-3 and GMOS spectra are shown in Figure 3. The GMOS spectrum has a S/N of 10 to 20 while the LDSS-3 has a lower S/N. The spectra exhibit band heads of CaH (7035 Å), TiO (7053, 8432 Å), CrH (8611 Å), FeH (8692, 9896 Å), and atomic lines of K i (7665, 7699 Å), Rb i (7800, 7948 Å), Na i (8183, 8195 Å), and Cs i (8521 Å). The heavily pressure-broadened K i doublet, strong CrH and FeH band heads, and Cs i and Rb i lines are hallmark spectral features of L dwarfs (Kirkpatrick et al. 1999) while the presence of the CaH and TiO band heads are consistent with known late-type M and L subdwarf spectra (Burgasser et al. 2007). Overall the spectra agree well except longward of ∼9200 Å where the LDSS-3 spectrum is systematically higher than the GMOS spectrum. Burgasser et al. (2007) found a similar offset between LDSS-3 and GMOS spectra of ultracool M subdwarfs and ascribed it to an unknown flux calibration error and/or the fact that the 9200 Å telluric H2O band was not corrected in their GMOS data. Our GMOS spectrum of 2MASS J0616−6407 has been corrected for telluric absorption so we can eliminate this as the cause of the mismatch. Since the cause remains unknown, we restrict further analysis to λ< 9200 Å.

Figure 3.

Figure 3. Red optical spectra of 2MASS J0616−6407 obtained with LDSS-3/Magellan (black) and GMOS/Gemini (red). The spectra are normalized to unity at 8100 Å.   Prominent atomic and molecular absorption features are indicated.

Standard image High-resolution image

2.4. Near-Infrared Imaging

An image of 2MASS J0616−6407 was taken on 2008 Jan 17 (UT) using the 1.5 m telescope at CTIO with the Caméra PAnoramique Proche InfraRouge (CPAPIR; Artigau et al. 2004) wide field-IR imager in the J band. The object was observed close to the meridian at an airmass of 1.2 and with seeing conditions ∼1farcs5. We used a nine-point dither pattern with 10'' offsets. Each exposure was ∼60 s amounting to a total integration time of ∼9 minutes. Science frames were first sky-subtracted using a sky frame created by median-combining all of the science data taken on a given night and then flat-fielded using a normalized dome flat. Individual frames were shifted and stacked to form the final combined image. The reduced science frame was astrometrically calibrated using the 2MASS Point Source catalog. We used a basic, six parameter, least-squares, linear transformation to obtain an astrometric solution for the CPAPIR image. The field-of-view of CPAPIR is ∼30'' so there were over 1600 2MASS reference stars to calibrate CPAPIR to 2MASS. We limited our reference stars to the 450 stars with 12 < J < 15 as objects in this intermediate magnitude range transformed with the smallest residuals from epoch to epoch. The solution reference stars were required to transform with total absolute residuals against 2MASS of < 0.2 pixels (see Faherty et al. 2009 for more details). We measure a proper motion of μα = +1404 ± 8 mas yr−1 and μα = −51 ± 18 mas yr−1 over the 9.1 yr baseline between the CPAPIR image and the discovery 2MASS image. The errors were computed by combining the residuals of the astrometric solution, which accounted for 0.16 pixels, with the plate scale of CPAPIR (∼1'' pixel−1) under the given baseline. The positional uncertainty for 2MASS J0616−6407 was also calculated by comparing the residuals of transforming the x, y position for our target over consecutive dithered images. However this uncertainty was negligible compared to the contribution from the astrometric solution.

3. ANALYSIS

3.1. Spectral Classification

Determining the spectral type of 2MASS J0616−6407 is not straightforward because an L subdwarf spectral classification scheme has yet to be rigorously defined in either the red optical or the near-infrared due to the few examples known. The depth of the 1.4 H2O band and the shape of the H band (see Figure 2) suggest that 2MASS J0616−6407 has a near-infrared spectral type intermediate between that of 2MASS J1626+3925 and 2MASS J0532+8246, perhaps sdL4 to sdL6. In order to derive a spectral type based on the red optical spectrum, we followed Burgasser et al. (2007) and compared the GMOS spectrum of 2MASS J0616−6407 to the spectra of field L dwarf spectral standards (Kirkpatrick et al. 1999). The best match is the L5 dwarf DENIS-P J1228.2−1547 resulting in an optical spectral type of sdL5 for 2MASS J0616−6407. Note that a similar comparison in the near-infrared is impossible because L dwarf spectral standards have yet to be defined at these wavelengths (see however J. Kirkpatrick et al. 2009, in preparation). Figure 4 shows the red optical spectra of both 2MASS J0616−6407 and DENIS -P J1228.2−1547. The strengths of the CrH and FeH band heads, width of the K i doublet, and overall 8000 Å-9000 Å slope of the two spectra agree well. Also shown in Figure 4 are the spectra of 2MASS J0532+8246 and 2MASS J1626+3925 along with their respective L dwarf standards. The strength of the features in the spectrum of 2MASS J0616−6407 are intermediate between that of 2MASS J1626+3925 and 2MASS J0532+8246, in good agreement with its derived spectral type of sdL5.

Figure 4.

Figure 4. Red optical spectra of 2MASS J1626+3925  (sdL4; Burgasser et al. 2007), 2MASS J0616−6407, and 2MASS J0532+8246  (sdL7; Burgasser et al. 2003) compared to those of 2MASS 1155009+230706 (L4; Kirkpatrick et al. 1999), DENIS-P J1228.2−1547 (L5; Kirkpatrick et al. 1999), and DENIS-P J0205.4−1159 (L7; Kirkpatrick et al. 1999), respectively. The spectra are normalized to unity and offset by constants (dotted lines). Prominent atomic and molecular absorption features are indicated.

Standard image High-resolution image

3.2. Kinematics

Although the high proper motion and large radial velocity of 2MASS J0616−6407 suggest that 2MASS J0616−6407 is a high-velocity star and thus possibly a member of the halo population, a distance is required in order to compute space velocities. We estimated a spectrophotometric distance to 2MASS J0616−6407 in the following way. First we performed a weighted linear fit of the 2MASS MJ, MH, and MKs values with respect to spectral type of the ultracool subdwarfs with known parallaxes (Monet et al. 1992; Dahn et al. 2008b; Schilbach et al. 2009; Burgasser et al. 2008), LHS 377 (sdM7), LSR 2036 + 5059 (sdM7.5), LSR J1425 + 7102 (sdM8), SSSPM 1013 − 1356 (sdM9.5), 2MASS J1626+3925 (sdL4) and 2MASS J0532+8246 (sdL7) and find,

Equation (1)

Equation (2)

Equation (3)

where SpType = 7 for M7, SpType = 10 for L0, etc., and Cov is the covariance matrix of the fit. For those dwarfs with two parallax measurements, we used the weighted average of the two values to compute their absolute magnitudes. Figure 5 shows the absolute magnitudes of the subdwarfs along with those of single field M and L dwarfs from Dahn et al. (2002) and references therein. In comparison to the field dwarfs, the L subdwarfs appear overluminous in the J band and underluminous in the K band. If we assume that the effective temperature scale of the L dwarfs and subdwarfs are similar (see however Burgasser & Kirkpatrick 2006; Burgasser et al. 2009), we can ascribe this behavior to a change in metallicity. At a fixed Teff and g, L dwarfs become brighter in the J band and fainter in the K band with decreasing metallicity (Burrows et al. 2006) due primarily to the increasing importance of CIA H2 absorption in the K band.

Figure 5.

Figure 5. 2MASS MJ, MH, and MKs values for LHS 377 (sdM7), LSR 2036+5059 (sdM7.5), LSR J1425+7102 (sdM8), SSSPM 1013−1356 (sdM9.5), 2MASS J1626+3925 (sdL4), and 2MASS J0532+8246 (sdL7) are shown in red. The linear fit to these values and the ±1 σ range are shown as solid and dotted lines. Shown in black are a sample of field dwarfs from Dahn et al. (2002) and references therein that have not been identified as binaries in high-resolution imaging. All spectral types are based on optical spectroscopy.

Standard image High-resolution image

We estimate a distance of 57 ± 9 pc for 2MASS J0616−6407 using the absolute J and H magnitudes derived using Equations (1) and (2) for a spectral type of sdL5 (2MASS J0616−6407 has only an upper limit in the Ks band). The value is the weighted average of the two spectrophotometric distance estimates and the error includes the covariance terms and an error of ±0.5 a subtype in the spectral type. The estimated distance and proper motion of 2MASS J0616−6407 imply a tangential velocity Vtan  of 379 ± 60 km s−1 which is roughly an order of magnitude greater than the median value of field L dwarfs (Vrba et al. 2004; Faherty et al. 2009) and confirms that 2MASS J0616−6407 is a member of the halo population. Indeed, Faherty et al. (2009) found no M, L, and T dwarfs with Vtan  > 170 km s−1 in their sample of 634 L and T dwarfs and 456 M dwarfs. We have also computed (u, v, w)9 velocities of 2MASS J0616−6407 with respect to the Local Standard of Rest (LSR) following Johnson & Soderblom (1987) updated with the J2000.0 galactic coordinate transformations of Murray (1989). Following the convention that the u is positive toward the Galactic center (l = 0, b = 0), we find (u, v, w)=(94, −573, 125)±(10, 31, 53) km s−1 after correcting for a solar velocity with respect to the LSR of (10.00, 5.25, 7.17) km s−1 (Dehnen & Binney 1998). The LSR moves at ΘLSR=220 km s−1 (Θ and v point in the direction of Galactic rotation, l = 90, b = 0; Kerr & Lynden-Bell 1986) in the rest frame of the Galaxy so 2MASS J0616−6407 has a highly retrograde orbit with Θ = −353 km s−1. The stars in the inner halo of the Galaxy have slightly prograde orbits with Θ=0 to 50 km s−1 while the stars in the outer halo have retrograde orbits with Θ = −40 to −70 km s−1 (Carollo et al. 2007). The rotation of 2MASS J0616−6407 around the center of the Galaxy clearly indicates that 2MASS J0616−6407 is a member of the outer halo.

Based on its (u, v, w) velocities we also computed the Galactic orbit of 2MASS J0616−6407 following Burgasser et al. (2009). We examined two simplified versions of the Galactic potential model described in Dehnen & Binney (1998) composed of spherically-symmetric halo and bulge mass distributions and an axisymmetric, thin exponential disk. The two Galactic mass models were parameterized according to Binney & Tremaine (2008, Table 2.3), which fit the measured rotation curve of the Galaxy but bracket the allowable ratios of disk to halo mass at the Solar radius. The orbit of 2MASS J0616−6407 was integrated using a second-order leapfrog method (kick-drift-kick) with a constant timestep of 1 kyr over ±1 Gyr centered on the present epoch. Energy was conserved to better than 1 part in 10−4 over the full length of the simulation, with the error dominated by the resolution of the grid on which the disk force and potential were interpolated. The z-component of angular momentum was conserved to 1 part in 10−13.

Figure 6 shows the orbit of 2MASS J0616-6407 based on the halo-dominant mass model, plotted in the rest frame of the Galaxy. The z-coordinate is defined to be positive toward the Galactic center to align with our definition of u, and the Sun is located at (−8.5, 0, 0.027) kpc (Kerr & Lynden-Bell 1986; Chen et al. 2001). Prograde motion in this reference frame is counterclockwise, so the orbit of 2MASS J0616−6407 is clearly retrograde, as well as wide and eccentric, with R extending from 8 to 36 kpc (e ∼ 0.7 from max/min R) over the simulation period. The orbit of 2MASS J0616−6407 also exhibits a fairly broad range of inclinations, with z spanning ±10 kpc (max |i| ∼ 15 degrees). For the disk-dominant mass model, the orbit is even wider (extending out to 40 kpc) but with similar eccentricity and inclination range. The breadth of this orbit and its retrograde motion support the conclusion that 2MASS J0616−6407 is a member of the outer halo, as this population dominates at R > 14–20 kpc (Carollo et al. 2007). We note that our conclusions remain unchanged if 2MASS J0616−6407 is an equal magnitude binary.

Figure 6.

Figure 6. Galactic orbit of 2MASS J0616−6407 over ±1 Gyr, based on the halo-dominant mass model of Binney & Tremaine (Model I; 2008). The upper panels show the orbit in (x, y) and (r, z) coordinates. The x-axis is along the line connecting the Galactic center and the Sun while the z-axis denotes height above the Galactic plane. The lower panels show the time evolution of r and z, with past motion denoted by dashed lines and future motion denoted by solid lines. The position of the Sun is denoted with a solid circle.

Standard image High-resolution image

4. DISCUSSION

The Galactic orbit of 2MASS J0616−6407 is distinctly different from that of other late-type subdwarfs (e.g., LSR J0822 + 1700, 2MASS J0532+8246, LSR J1425 + 7102, SDSS J1256−0224) which all have perigalacticons of ≲1 kpc and apogalacticons near the Sun (Lépine et al. 2004; Burgasser et al. 2008, 2009; Dahn et al. 2008a; Burgasser et al. 2009). As noted by Burgasser et al. (2009), these other subdwarfs have kinematic properties which are consistent with membership in the inner halo of the Galaxy. 2MASS J0616−6407 is therefore the first late-type subdwarf that is a member of the outer halo and therefore may represent a new class of L subdwarfs. In particular, the inner halo comprises stars with a peak metallicity of [Fe/H] ∼−1.6 while the outer halo has a peak near −2.2 (Carollo et al. 2007). This suggests that 2MASS J0616−6407 may be even more metal-poor than the other late-type subdwarfs known. However, the techniques needed to derive accurate metallicities of late-type subdwarfs are in their infancy (e.g., Burgasser et al. 2007, 2009) and consequently an accurate measurement of the metallicity of 2MASS J0616−6407 must await future work. Finally we note that the discovery of additional outer halo ultracool subdwarfs will be difficult given their intrinsic faintness and the large amount of time that they spend away from the solar neighborhood.

We thank Bill Vacca for fruitful discussions about the reduction of optical spectra and John Bochanski and Andrew West for clarifications on the Galactic orbit geometry. This publication makes use of data from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center, and funded by the National Aeronautics and Space Administration and the National Science Foundation, the SIMBAD database, operated at CDS, Strasbourg, France, NASA's Astrophysics Data System Bibliographic Services, the M, L, and T dwarf compendium housed at DwarfArchives.org and maintained by Chris Gelino, Davy Kirkpatrick, and Adam Burgasser, and the NASA/ IPAC Infrared Science Archive, which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. OSIRIS is a collaborative project between the Ohio State University and Cerro Tololo Inter-American Observatory (CTIO) and was developed through NSF grants AST 90-16112 and AST 92-18449. CTIO is part of the National Optical Astronomy Observatory (NOAO), based in La Serena, Chile. NOAO is operated by the Association of Universities for Research in Astronomy (AURA), Inc. under cooperative agreement with the National Science Foundation.

Facilities: Gemini (GMOS), SOAR (OSIRIS), LDSS-3 (Magellan), CTIO (CPAPIR).

Footnotes

  • A compendium of known L and T dwarfs is given at http://DwarfArchives.org.

  • LSR 1610−0040 was originally classified as an L subdwarf by Lépine et al. (2003) but observations by Cushing & Vacca (2006), Reiners & Basri (2006), and Dahn et al. (2008b) indicate that it is probably a peculiar M dwarf/subdwarf binary.

  • IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation.

  • We follow the convention that lower case letters refer to velocities with respect to the LSR while upper case letters refer to velocities with respect to the Sun (Mihalas & Binney 1981, page 608).

Please wait… references are loading.
10.1088/0004-637X/696/1/986