# Primeval very low-mass stars and brown dwarfs – IV. New L subdwarfs, Gaia astrometry, population properties, and a blue brown dwarf binary

Z. H. Zhang, M. C. Galvez-Ortiz, D. J. Pinfield, A. J. Burgasser, N. Lodieu, H. R. A. Jones, E. L. Martín, B. Burningham, D. Homeier, F. Allard, et al.

<figure><img src="image_1.png" /></figure>

► To cite this version:

Z. H. Zhang, M. C. Galvez-Ortiz, D. J. Pinfield, A. J. Burgasser, N. Lodieu, et al. Primeval very low-mass stars and brown dwarfs – IV. New L subdwarfs, Gaia astrometry, population properties, and a blue brown dwarf binary. Monthly Notices of the Royal Astronomical Society, 2018, 480 (4), pp.5447-5474. 10.1093/mnras/sty2054. hal-02280770

HAL Id: hal-02280770

https://hal.science/hal-02280770v1

Submitted on 6 Sep 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

# Primeval very low-mass stars and brown dwarfs – IV. New L subdwarfs, Gaia astrometry, population properties, and a blue brown dwarf binary

Z. H. Zhang (张曾华) ID,<fnref n="1" />,<fnref n="2" />,<fnref n="3" />★† M. C. Galvez-Ortiz,<fnref n="4" /> D. J. Pinfield,<fnref n="5" /> A. J. Burgasser,<fnref n="6" />

N. Lodieu,<fnref n="2" />,<fnref n="3" /> H. R. A. Jones,<fnref n="5" /> E. L. Martín,<fnref n="4" /> B. Burningham,<fnref n="5" /> D. Homeier ID,<fnref n="7" />

F. Allard,<fnref n="8" /> M. R. Zapatero Osorio,<fnref n="4" /> L. C. Smith,<fnref n="5" />,<fnref n="9" /> R. L. Smart ID,<fnref n="10" />

B. López Martí,<fnref n="11" /> F. Marocco<fnref n="5" />,<fnref n="12" />‡ and R. Rebolo<fnref n="2" />,<fnref n="3" />

Accepted 2018 July 27. Received 2018 June 28; in original form 2018 March 28

## ABSTRACT

We present 27 new L subdwarfs and classify five of them as esdL and 22 as sdL. Our L subdwarf candidates were selected with the UKIRT Infrared Deep Sky Survey and Sloan Digital Sky Survey. Spectroscopic follow-up was carried out primarily with the OSIRIS spectrograph on the Gran Telescopio Canarias. Some of these new objects were followed up with the X-shooter instrument on the Very Large Telescope. We studied the photometric properties of the population of known L subdwarfs using colour-spectral type diagrams and colour-colour diagrams, by comparison with L dwarfs and main sequence stars, and identified new colour spaces for L subdwarf selection/study in current and future surveys. We further discussed the brown dwarf transition-zone and the observational stellar/substellar boundary. We found that about one-third of 66 known L subdwarfs are substellar objects, with two-thirds being very low-mass stars. We also present the Hertzsprung-Russell diagrams, spectral type-absolute magnitude corrections, and tangential velocities of 20 known L subdwarfs observed by the Gaia astrometry satellite. One of our L subdwarf candidates, ULAS J233227.03+123452.0, is a mildly metal-poor spectroscopic binary brown dwarf: a ~L6p dwarf and a ~T4p dwarf. This binary is likely a thick disc member according to its kinematics.

**Key words:** binaries: spectroscopic – brown dwarfs – stars: chemically peculiar – stars: Population II – stars: subdwarfs.

## 1 INTRODUCTION

Very low-mass stars (VLMS) and brown dwarfs (BDs) are classified using M, L, T, and Y types according to spectral morphology that is dominated by temperature-dependent chemistry and thermal

properties (Bessell 1991; Kirkpatrick, Henry & McCarthy 1991; Jones et al. 1994; Tsuji et al. 1996; Kirkpatrick et al. 1999a, 2012; Martín et al. 1999; Burgasser et al. 2002, 2003b; Cushing et al. 2011). The spectral types of field VLMS extend from M to early-type L (Dieterich et al. 2014; Dupuy & Liu 2017). Most massive BDs have spectral types between late-type M and mid-type L depending on their age (Zhang et al. 2017b). For 0.1–1 Gyr massive BDs are late-type M, evolving to later type as they cool (L type for

<footnote n="1">¹GEPI, Observatoire de Paris, Université PSL, CNRS, 5 Place Jules Janssen, F-92190 Meudon, France</footnote>

<footnote n="2">²Instituto de Astrofísica de Canarias, E-38205 La Laguna, Tenerife, Spain</footnote>

<footnote n="3">³Universidad de La Laguna, Dept. Astrofísica, E-38206 La Laguna, Tenerife, Spain</footnote>

<footnote n="4">⁴Centro de Astrobiología (CSIC-INTA), Ctra. Ajalvir km 4, E-28850 Torrejón de Ardoz, Madrid, Spain</footnote>

<footnote n="5">⁵Centre for Astrophysics Research, Science and Technology Research Institute, University of Hertfordshire, Hatfield AL10 9AB, UK</footnote>

<footnote n="6">⁶Center for Astrophysics and Space Science, University of California San Diego, La Jolla, CA 92093, USA</footnote>

<footnote n="7">⁷Zentrum für Astronomie der Universität Heidelberg, Landessternwarte, Königstuhl 12, D-69117 Heidelberg, Germany</footnote>

<footnote n="8">⁸Univ Lyon, ENS de Lyon, Univ Lyon 1, CNRS, Centre de Recherche Astrophysique de Lyon, UMR5574, F-69007 Lyon, France</footnote>

<footnote n="9">⁹Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK</footnote>

<footnote n="10">¹⁰Istituto Nazionale di Astrofisica, Osservatorio Astronomico di Torino, Strada Osservatorio 20, I-10025 Pino Torinese, Italy</footnote>

<footnote n="11">¹¹Saint Louis University - Madrid Campus, Avenida del Valle 34, E-28003 Madrid, Spain</footnote>

<footnote n="12">¹²Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Dr., Pasadena, CA 91109, USA</footnote>

<footnote>* E-mail: zenghuazhang@hotmail.com</footnote>

<footnote>† PSL Fellow.</footnote>

<footnote>‡ NASA Postdoctoral Program Fellow.</footnote>
~1–10 Gyr age, and T type for ages >10 Gyr; e.g. fig. 8 of Burrow
s et al. 2001).

The classification of VLMS and BDs with subsolar metallicity
is complicated by the population's wide metallicity range, which
leads to significant diversity in the strength of spectral features.
And available sample size also places some limitations on the scope
of classification (Gizis 1997; Kirkpatrick 2005; Burgasser, Cruz &
Kirkpatrick 2007; Lépine, Rich & Shara 2007; Jao et al. 2008;
Dhital et al. 2012; Kirkpatrick et al. 2014; Zhang et al. 2017a).
However, previous work has established a set of subdwarf classes
that are indicative of metallicity. Gizis (1997) refined M classifica-
tion by setting up three subclasses: dwarf (dM), subdwarf (sdM),
and extreme subdwarf (esdM). An additional intermediate subclass
d/sd for late-type M and L subdwarfs was also suggested by Bur-
gasser et al. (2007), aimed at intermediate metallicity to further en-
compasses the Galactic thick disc population. Lépine et al. (2007)
revised the M dwarf classification into four subclasses: dM, sdM,
esdM, and ultra-subdwarf (usdM).

Kirkpatrick (2005) proposed a three-parameter classification
strategy for late-type M, L, and T dwarfs to indicate their metal-
licity, temperature/clouds, and gravity. For example, a metal-poor
halo object 2MASS J05325346+8246465 (Burgasser et al. 2003a)
was classified esdL7 and an intermediate metal-poor object SDSS
J141624.12+134827.4 (SD1416; Bowler, Liu & Dupuy 2010; Burn-
ingham et al. 2010; Schmidt et al. 2010) was classified sdL7 in Kirk-
patrick et al. (2010). And more recently Zhang et al. (2017a) classi-
fies L subdwarfs using three subclasses: sdL, esdL, and usdL, based
on the relative strength of subsolar metallicity sensitive spectral
features across the optical and near-infrared (NIR). The metallicity
ranges of the usdL, esdL, and sdL subclasses are approximately
[Fe/H] ≲ −1.7: −1.7 ≲ [Fe/H] ≲ −1.0; and −1.0 ≲ [Fe/H] ≲
−0.3, respectively (see Zhang et al. 2017a).

The L subdwarf population is composed of metal-deficient low-
mass stars and high-mass BDs with effective temperature (T<sub>eff</sub>) in
the range of ~1300-2700 K. They have strong FeH absorption
bands, weak or absent VO and CO bands, and enhanced collision-
induced H<sub>2</sub> absorption (CIA H<sub>2</sub>; Bates 1952; Saumon et al. 2012), as
summarized in Zhang et al. (2017a). Sample identification has been
enabled by modern large-scale optical and NIR surveys: the Two
Micron All Sky Survey (2MASS; Skrutskie et al. 2006), the Sloan
Digital Sky Survey (SDSS; York et al. 2000), the UKIRT Infrared
Deep Sky Survey (UKIDSS; Lawrence et al. 2007), and the *Wide-
field Infrared Survey Explorer* (WISE; Wright et al. 2010), with a
total of 39 L subdwarfs previously reported in the literature (Lépine
et al. 2002; Burgasser et al. 2003a, 2004a; Burgasser 2004b; Scholz,
Lodieu & McCaughrean 2004a; Scholz et al. 2004b; Burgasser
& Kirkpatrick 2006; Gizis & Harvin 2006; Cushing et al. 2009;
Sivarami et al. 2009; Bowler et al. 2010; Burningham et al. 2010;
Kirkpatrick et al. 2010, 2014, 2016; Lodieu et al. 2010, 2012, 2017;
Schmidt et al. 2010; Luhman & Sheppard 2014; Schneider et al.
2016; Zhang et al. 2017a,b, 2018; Smith et al. 2018).

In general, subdwarfs with sdM–sdL subclass are kinematically
associated with the Galactic thick disc, while those with esdM–
esdL and usdM–usdL subclasses have kinematics of the Galactic
halo (Burgasser et al. 2008b; Cushing et al. 2009; Zhang et al.
2013, 2017a). By comparison to the population of known L dwarfs,
L subdwarfs in the solar neighbourhood are rare. This is partly
because they are predominantly thick disc and halo objects (cf.
Reddy, Lambert & Allende Prieto 2006 have shown that the frac-
tions of thick disc and halo stars in the solar neighbourhood are
7 and 0.6 per cent), but also because L subdwarfs occupy a nar-
rower mass range than L dwarfs (Zhang et al. 2017b). Lodieu et al.

(2017) report a UKIDSS/SDSS surface density (for late M and
early L dwarfs) of around 0.04 per deg², pointing to potential for
a much larger (hundreds strong) detectable population across sur-
veyed sky. Inroads into this population are important for the detailed
study of observed diversity and population make-up amongst the
L subdwarfs.

This is the fourth paper in a series titled *Primeval very low-
mass stars and brown dwarfs*. The first paper reported the discovery
of six new L subdwarfs, defined a new L subdwarf classification
scheme, and studied the atmospheric properties of ultra-cool sub-
dwarfs based on 22 late-type M and L subdwarfs (Zhang et al.
2017a, hereafter Paper I). In the second paper, we presented the
most metal-poor substellar object, and a procedure to distinguish
massive halo BDs from the least-massive stars. We also found that
mid-type L to early-type T subdwarfs of the Galactic halo are lo-
cated in a substellar subdwarf gap, known also as the halo BD
transition-zone. This zone covers a narrow mass range but spans a
wide T<sub>eff</sub> range due to unsteady nuclear fusion (Zhang et al. 2017b,
hereafter Paper II). In the third paper, we presented the discovery
of three new halo transitional BDs (Zhang et al. 2018, hereafter
Paper III). Here, we present the discovery of 27 new L subdwarfs
and their population properties, as well as a spectroscopic blue BD
binary. Observations are presented in Section 2. Classification is
carried out in Section 3. Section 4 assesses L subdwarf population
properties. Section 5 presents Gaia observations of L subdwarfs.
Section 6 sums up and presents our conclusions.

2 OBSERVATION

2.1 Candidate selection

The L subdwarf candidates in our programme were generally
selected from the UKIDSS Large Area Survey (LAS) and the
SDSS, following the selection criteria/procedure described in Paper
I. In addition two extra candidates without SDSS detections
were selected from a ULAS proper motion catalogue, covering
1500 deg² to an approximate 5σ depth of *J* = 19.6 (Smith et al.
2014). These two objects (ULAS J130710.22+151103.4 and ULAS
J144151.55+043738.5) both have relatively blue *J* − *K* colour and
proper motion higher than 0.4 arcsec per year.

Our programmatic candidate sample consisted of 64 objects of
which five had been confirmed as L subdwarfs by other studies. We
have spectroscopically followed up the remaining 59 candidates,
which has in total resulted in the confirmation/study of 34 L subd-
warfs (an overall 61 per cent confirmation rate). Eight of these were
published in Papers I and III, and the remaining 27 are reported
here. Table 1 presents optical and NIR photometry for all known L
subdwarfs, which includes the 34 studied by our programme plus 31
from elsewhere in the literature. The 24 programmatic candidates
that are not L subdwarfs consist of fourteen late M dwarfs, one L1
dwarf, six late type sdM subdwarfs, 1 esdM7 subdwarf, a probable
L+T unresolved binary (see Section 3.4), and one galaxy.

2.2 GTC spectroscopy

Candidates have been followed up with optical spectroscopy us-
ing the Optical System for Imaging and low-Resolution Integrated
Spectroscopy (OSIRIS; Cepa et al. 2000) instrument on the Gran
Telescopio Canarias (GTC) since 2012. Table 2 shows a summary
of our GTC spectroscopic observations. The spectra were mostly
obtained using the R500R grism, with three making use of the
R300R grism. These provide resolving power of around 500 and
<table><caption>Table 1. Optical and NIR photometry of L subdwarfs.</caption><thead><tr><th>Name</th><th>SpT</th><th>SDSS i</th><th>SDSS z</th><th>Y (MKO)</th><th>J (MKO)</th><th>H (MKO)</th><th>K (MKO)</th><th>Ref<sup>a</sup></th></tr></thead><tbody><tr><td>SSSSS J1010448.47+153501.9</td><td>usdL1.5</td><td>20.37 ± 0.05</td><td>19.29 ± 0.06</td><td>18.48 ± 0.05</td><td>17.93 ± 0.05</td><td>18.06 ± 0.11</td><td>18.08 ± 0.17</td><td>17.24</td></tr><tr><td>SSSPJM J10130734–1356204<sup>b</sup></td><td>usdL0</td><td>17.19 ± 0.01</td><td>16.09 ± 0.01</td><td>15.16 ± 0.01</td><td>14.56 ± 0.01</td><td>14.39 ± 0.01</td><td>14.30 ± 0.01</td><td>19.5</td></tr><tr><td>SDSS J125637.13–022452.4</td><td>usdL3</td><td>19.40 ± 0.02</td><td>17.71 ± 0.02</td><td>16.77 ± 0.01</td><td>16.08 ± 0.01</td><td>16.05 ± 0.01</td><td>16.08 ± 0.01</td><td>21.6</td></tr><tr><td>ULAS J135058.85+081506.8</td><td>usdL3</td><td>21.22 ± 0.08</td><td>19.47 ± 0.06</td><td>18.66 ± 0.05</td><td>17.93 ± 0.04</td><td>18.07 ± 0.10</td><td>17.95 ± 0.15</td><td>15.11</td></tr><tr><td>2MASS J16262034+3925190</td><td>usdL4</td><td>17.90 ± 0.01</td><td>16.16 ± 0.01</td><td>–</td><td>14.37 ± 0.11</td><td>14.52 ± 0.10</td><td>14.46 ± 0.09</td><td>3</td></tr><tr><td>WISEA J213409.15+713236.1</td><td>usdL0.5</td><td>–</td><td>–</td><td>–</td><td>16.12 ± 0.15</td><td>16.30 ± 0.28</td><td>≥16.64</td><td>12</td></tr><tr><td>ULAS J230711.01+014447.1</td><td>usdL4.5</td><td>22.51 ± 0.24</td><td>19.81 ± 0.09</td><td>18.99 ± 0.08</td><td>18.15 ± 0.06</td><td>18.34 ± 0.12</td><td>18.17 ± 0.18</td><td>25</td></tr><tr><td>WISEA J001450.17–083823.4</td><td>esdL0</td><td>17.43 ± 0.01</td><td>16.05 ± 0.01</td><td>–</td><td>14.42 ± 0.11</td><td>14.00 ± 0.11</td><td>13.74 ± 0.10</td><td>11.14</td></tr><tr><td>WISEA J020201.25–313645.2<sup>c</sup></td><td>esdL0.5</td><td>18.44 ± 0.01</td><td>16.99 ± 0.01</td><td>15.89 ± 0.01</td><td>15.15 ± 0.01</td><td>14.96 ± 0.01</td><td>14.80 ± 0.01</td><td>11</td></tr><tr><td>ULAS J020858.62+020657.0</td><td>esdL3</td><td>21.55 ± 0.08</td><td>19.83 ± 0.07</td><td>18.76 ± 0.05</td><td>18.00 ± 0.04</td><td>17.88 ± 0.13</td><td>17.62 ± 0.16</td><td>25</td></tr><tr><td>WISEA J030601.66–033059.0</td><td>esdL1</td><td>–</td><td>–</td><td>–</td><td>14.39 ± 0.11</td><td>14.11 ± 0.11</td><td>13.96 ± 0.11</td><td>11.14</td></tr><tr><td>ULAS J033351.10+001405.8</td><td>esdL0</td><td>19.24 ± 0.01</td><td>17.87 ± 0.02</td><td>16.81 ± 0.01</td><td>16.11 ± 0.01</td><td>15.77 ± 0.01</td><td>15.50 ± 0.02</td><td>16</td></tr><tr><td>WISEA J043535.82+211508.9</td><td>esdL1</td><td>18.59 ± 0.01</td><td>16.96 ± 0.01</td><td>–</td><td>14.96 ± 0.11</td><td>14.73 ± 0.12</td><td>14.57 ± 0.12</td><td>11.14</td></tr><tr><td>2MASS J05325346+824646.5</td><td>esdL7</td><td>20.36 ± 0.05</td><td>17.59 ± 0.02</td><td>–</td><td>15.08 ± 0.11</td><td>14.96 ± 0.10</td><td>14.92 ± 0.09</td><td>1</td></tr><tr><td>2MASS J06164006–6407194</td><td>esdL6</td><td>–</td><td>–</td><td>–</td><td>16.34 ± 0.16</td><td>16.34 ± 0.25</td><td>≥16.39</td><td>8</td></tr><tr><td>ULAS J111429.54+072809.5</td><td>esdL0</td><td>20.62 ± 0.05</td><td>19.26 ± 0.06</td><td>18.29 ± 0.03</td><td>17.59 ± 0.03</td><td>17.26 ± 0.04</td><td>17.12 ± 0.07</td><td>26</td></tr><tr><td>ULAS J124425.75+102439.3</td><td>esdL0.5</td><td>19.49 ± 0.02</td><td>18.01 ± 0.02</td><td>16.98 ± 0.01</td><td>16.26 ± 0.01</td><td>16.00 ± 0.01</td><td>15.77 ± 0.02</td><td>17</td></tr><tr><td>SDSS J124410.11+273625.8</td><td>esdL0.5</td><td>20.40 ± 0.05</td><td>19.12 ± 0.06</td><td>18.28 ± 0.04</td><td>17.58 ± 0.03</td><td>17.32 ± 0.05</td><td>17.12 ± 0.06</td><td>16</td></tr><tr><td>ULAS J135216.31+312327.0</td><td>esdL0.5</td><td>20.01 ± 0.04</td><td>18.66 ± 0.05</td><td>17.69 ± 0.02</td><td>16.93 ± 0.02</td><td>16.66 ± 0.03</td><td>16.41 ± 0.04</td><td>26</td></tr><tr><td>SDSS J141405.74–014202.7</td><td>esdL0</td><td>19.85 ± 0.03</td><td>18.45 ± 0.03</td><td>17.50 ± 0.03</td><td>16.81 ± 0.02</td><td>16.45 ± 0.03</td><td>16.14 ± 0.03</td><td>17</td></tr><tr><td>SSSPJM J144420.67–201922.2<sup>d</sup></td><td>esdL1</td><td>–</td><td>–</td><td>13.23 ± 0.01</td><td>12.45 ± 0.01</td><td>12.19 ± 0.11</td><td>11.93 ± 0.01</td><td>18.12</td></tr><tr><td>ULAS J145234.65+043738.4</td><td>esdL0.5</td><td>20.50 ± 0.06</td><td>18.92 ± 0.05</td><td>18.11 ± 0.03</td><td>17.28 ± 0.03</td><td>16.79 ± 0.04</td><td>16.62 ± 0.06</td><td>26</td></tr><tr><td>ULAS J151913.03–000030.0</td><td>esdL4</td><td>21.42 ± 0.09</td><td>19.26 ± 0.06</td><td>18.19 ± 0.03</td><td>17.21 ± 0.02</td><td>17.07 ± 0.03</td><td>16.97 ± 0.04</td><td>23</td></tr><tr><td>2MASS J16403197+1231068</td><td>esdL0</td><td>19.01 ± 0.02</td><td>17.56 ± 0.02</td><td>–</td><td>15.90 ± 0.13</td><td>15.65 ± 0.15</td><td>15.49 ± 0.18</td><td>2.9</td></tr><tr><td>WISEA J204027.30+695924.1</td><td>esdL0.5</td><td>–</td><td>–</td><td>–</td><td>13.68 ± 0.13</td><td>13.36 ± 0.12</td><td>13.09 ± 0.11</td><td>11.14</td></tr><tr><td>ULAS J223302.03+062030.8</td><td>esdL0.5</td><td>21.12 ± 0.06</td><td>19.68 ± 0.06</td><td>18.60 ± 0.04</td><td>17.90 ± 0.05</td><td>17.67 ± 0.08</td><td>17.44 ± 0.11</td><td>26</td></tr><tr><td>ULAS J231924.35+052524.5</td><td>esdL1</td><td>20.41 ± 0.05</td><td>18.99 ± 0.05</td><td>18.22 ± 0.03</td><td>17.33 ± 0.02</td><td>17.08 ± 0.04</td><td>16.93 ± 0.06</td><td>26</td></tr><tr><td>2MASS J00412179+3547133</td><td>sdL0.5</td><td>20.28 ± 0.03</td><td>18.33 ± 0.02</td><td>–</td><td>15.89 ± 0.13</td><td>15.78 ± 0.18</td><td>15.14 ± 0.15</td><td>2</td></tr><tr><td>WISEA J005757.65+201304.0</td><td>sdL7</td><td>21.19 ± 0.11</td><td>18.97 ± 0.06</td><td>–</td><td>16.26 ± 0.14</td><td>15.52 ± 0.14</td><td>14.99 ± 0.16</td><td>11.14</td></tr><tr><td>WISEA J011639.05–165420.5<sup>e</sup></td><td>sdL0</td><td>19.49 ± 0.02</td><td>17.66 ± 0.02</td><td>–</td><td>15.76 ± 0.12</td><td>15.38 ± 0.14</td><td>14.94 ± 0.16</td><td>20</td></tr><tr><td>ULAS J011824.89+034130.4</td><td>sdL0</td><td>21.93 ± 0.20</td><td>20.47 ± 0.21</td><td>19.10 ± 0.07</td><td>18.18 ± 0.05</td><td>17.73 ± 0.07</td><td>17.56 ± 0.11</td><td>26</td></tr><tr><td>WISEA J130132.66–104732.4</td><td>sdL0</td><td>19.33 ± 0.03</td><td>17.50 ± 0.02</td><td>–</td><td>15.58 ± 0.12</td><td>15.18 ± 0.14</td><td>14.80 ± 0.15</td><td>20</td></tr><tr><td>ULAS J021258.08+064115.9</td><td>sdL1</td><td>21.11 ± 0.08</td><td>19.38 ± 0.08</td><td>18.20 ± 0.03</td><td>17.43 ± 0.03</td><td>17.06 ± 0.03</td><td>16.78 ± 0.05</td><td>26</td></tr><tr><td>ULAS J021642.96+004005.7</td><td>sdL4</td><td>22.10 ± 0.16</td><td>19.99 ± 0.10</td><td>18.41 ± 0.05</td><td>17.30 ± 0.03</td><td>16.96 ± 0.04</td><td>16.51 ± 0.04</td><td>23</td></tr><tr><td>ULAS J023803.12+054526.1</td><td>sdL0</td><td>20.48 ± 0.04</td><td>18.57 ± 0.03</td><td>17.26 ± 0.02</td><td>16.43 ± 0.01</td><td>16.02 ± 0.02</td><td>15.59 ± 0.02</td><td>26</td></tr><tr><td>2MASS J06453153–6646120</td><td>sdL8</td><td>–</td><td>–</td><td>–</td><td>15.54 ± 0.13</td><td>14.77 ± 0.12</td><td>14.40 ± 0.13</td><td>10</td></tr><tr><td>ULAS J075335.23+200622.4</td><td>sdL0</td><td>19.76 ± 0.02</td><td>17.90 ± 0.02</td><td>16.68 ± 0.01</td><td>15.87 ± 0.01</td><td>15.48 ± 0.01</td><td>15.09 ± 0.01</td><td>26</td></tr><tr><td>ULAS J082206.61+044101.8</td><td>sdL0</td><td>20.33 ± 0.04</td><td>18.51 ± 0.03</td><td>17.12 ± 0.02</td><td>16.29 ± 0.01</td><td>15.96 ± 0.02</td><td>15.53 ± 0.02</td><td>26</td></tr><tr><td>WISEA J101329.72–724619.2<sup>f</sup></td><td>sdL2?</td><td>–</td><td>–</td><td>–</td><td>15.89 ± 0.01</td><td>15.49 ± 0.15</td><td>14.94 ± 0.02</td><td>12</td></tr><tr><td>2MASS J11582077+0435014</td><td>sdL7</td><td>21.08 ± 0.05</td><td>18.18 ± 0.02</td><td>16.61 ± 0.01</td><td>15.43 ± 0.01</td><td>14.88 ± 0.01</td><td>14.37 ± 0.01</td><td>10</td></tr><tr><td>ULAS J123142.99+015045.4</td><td>sdL0</td><td>21.38 ± 0.11</td><td>19.45 ± 0.09</td><td>18.46 ± 0.04</td><td>17.54 ± 0.03</td><td>17.21 ± 0.04</td><td>16.78 ± 0.05</td><td>26</td></tr><tr><td>ULAS J124104.75–000531.4</td><td>sdL0</td><td>22.41 ± 0.27</td><td>20.42 ± 0.17</td><td>19.12 ± 0.08</td><td>18.46 ± 0.10</td><td>18.08 ± 0.12</td><td>18.04 ± 0.20</td><td>26</td></tr><tr><td>ULAS J124947.04+095019.8</td><td>sdL1</td><td>20.40 ± 0.04</td><td>18.66 ± 0.04</td><td>17.62 ± 0.02</td><td>16.83 ± 0.02</td><td>16.40 ± 0.03</td><td>16.12 ± 0.04</td><td>23</td></tr><tr><td>ULAS J125226.62+092920.1</td><td>sdL0</td><td>20.78 ± 0.05</td><td>19.03 ± 0.05</td><td>17.69 ± 0.02</td><td>16.87 ± 0.02</td><td>16.47 ± 0.02</td><td>16.05 ± 0.03</td><td>26</td></tr><tr><td>VVV J12564163–6202039</td><td>sdL3<sup>g</sup></td><td>19.70 ± 0.08</td><td>–</td><td>17.10 ± 0.02</td><td>16.10 ± 0.01</td><td>15.89 ± 0.02</td><td>15.72 ± 0.03</td><td>22</td></tr><tr><td>ULAS J130710.22+151103.4</td><td>sdL8</td><td>–</td><td>–</td><td>19.31 ± 0.07</td><td>18.14 ± 0.04</td><td>17.53 ± 0.07</td><td>17.24 ± 0.09</td><td>26</td></tr><tr><td>ULAS J133348.27+273505.5</td><td>sdL1</td><td>20.52 ± 0.05</td><td>18.76 ± 0.04</td><td>17.47 ± 0.02</td><td>16.62 ± 0.01</td><td>16.27 ± 0.02</td><td>15.98 ± 0.02</td><td>23</td></tr><tr><td>ULAS J133836.97–022910.7</td><td>sdL7</td><td>22.53 ± 0.28</td><td>20.10 ± 0.15</td><td>18.56 ± 0.06</td><td>17.37 ± 0.03</td><td>16.81 ± 0.04</td><td>16.37 ± 0.05</td><td>23</td></tr><tr><td>ULAS J134206.86+053724.9</td><td>sdL0.5</td><td>21.95 ± 0.16</td><td>19.62 ± 0.08</td><td>18.38 ± 0.04</td><td>17.43 ± 0.03</td><td>17.01 ± 0.03</td><td>16.56 ± 0.04</td><td>26</td></tr><tr><td>ULAS J134423.98+280603.8</td><td>sdL4</td><td>22.67 ± 0.28</td><td>19.85 ± 0.10</td><td>18.40 ± 0.03</td><td>17.19 ± 0.02</td><td>16.72 ± 0.03</td><td>16.13 ± 0.03</td><td>26</td></tr><tr><td>ULAS J134749.79+333601.7</td><td>sdL0</td><td>19.87 ± 0.03</td><td>18.07 ± 0.02</td><td>16.66 ± 0.01</td><td>15.85 ± 0.01</td><td>15.46 ± 0.01</td><td>15.27 ± 0.02</td><td>23</td></tr><tr><td>ULAS J134852.93+101611.8</td><td>sdL0</td><td>20.76 ± 0.05</td><td>18.99 ± 0.04</td><td>17.54 ± 0.02</td><td>16.69 ± 0.01</td><td>16.32 ± 0.01</td><td>15.87 ± 0.02</td><td>26</td></tr><tr><td>ULAS J135359.58+011856.7</td><td>sdL0</td><td>21.49 ± 0.14</td><td>19.37 ± 0.08</td><td>18.26 ± 0.04</td><td>17.36 ± 0.03</td><td>16.89 ± 0.03</td><td>16.52 ± 0.04</td><td>26</td></tr><tr><td>WISEA J135501.90–825838.9<sup>f</sup></td><td>sdL5?</td><td>–</td><td>–</td><td>–</td><td>16.32 ± 0.01</td><td>15.37 ± 0.17</td><td>14.90 ± 0.01</td><td>12</td></tr><tr><td>ULAS J141203.85+121609.9</td><td>sdL5</td><td>21.29 ± 0.09</td><td>19.14 ± 0.05</td><td>17.54 ± 0.02</td><td>16.33 ± 0.01</td><td>15.85 ± 0.01</td><td>15.43 ± 0.02</td><td>26</td></tr><tr><td>SDSS J141624.12+134827.4</td><td>sdL7</td><td>18.38 ± 0.01</td><td>15.91 ± 0.01</td><td>14.26 ± 0.01</td><td>12.99 ± 0.01</td><td>12.47 ± 0.01</td><td>12.05 ± 0.01</td><td>7.10</td></tr><tr><td>ULAS J141832.35+025323.0</td><td>sdL0</td><td>20.11 ± 0.04</td><td>18.27 ± 0.03</td><td>16.86 ± 0.01</td><td>16.00 ± 0.01</td><td>15.61 ± 0.01</td><td>15.19 ± 0.01</td><td>26</td></tr><tr><td>ULAS J144151.55+043738.5</td><td>sdL4</td><td>–</td><td>–</td><td>18.45 ± 0.04</td><td>17.23 ± 0.03</td><td>16.95 ± 0.04</td><td>16.34 ± 0.04</td><td>26</td></tr><tr><td>ULAS J151649.84+083607.1</td><td>sdL5</td><td>22.65 ± 0.24</td><td>20.25 ± 0.10</td><td>18.74 ± 0.04</td><td>17.35 ± 0.02</td><td>16.71 ± 0.03</td><td>16.26 ± 0.03</td><td>26</td></tr><tr><td>ULAS J154638.34–011213.0</td><td>sdL3</td><td>22.10 ± 0.14</td><td>19.98 ± 0.11</td><td>18.56 ± 0.05</td><td>17.51 ± 0.04</td><td>17.21 ± 0.05</td><td>16.95 ± 0.08</td><td>26</td></tr><tr><td>2MASS J17561080+2815238</td><td>sdL1</td><td>–</td><td>–</td><td>–</td><td>14.66 ± 0.11</td><td>14.19 ± 0.11</td><td>13.79 ± 0.10</td><td>10</td></tr><tr><td>LSR J182611.3+301419.1</td><td>sdL0</td><td>–</td><td>–</td><td>–</td><td>11.61 ± 0.11</td><td>11.22 ± 0.10</td><td>10.78 ± 0.10</td><td>13</td></tr><tr><td>ULAS J223440.80+001002.6</td><td>sdL1</td><td>22.05 ± 0.14</td><td>20.16 ± 0.12</td><td>18.97 ± 0.07</td><td>17.63 ± 0.04</td><td>17.15 ± 0.05</td><td>16.90 ± 0.07</td><td>26</td></tr></tbody></table>
<table><caption>**Table 1** – continued</caption><thead><tr><th>Name</th><th>SpT</th><th>SDSS i</th><th>SDSS z</th><th>Y (MKO)</th><th>J (MKO)</th><th>H (MKO)</th><th>K (MKO)</th><th>Ref<sup>a</sup></th></tr></thead><tbody><tr><td>ULAS J225902.14+115602.1</td><td>sdL0</td><td>21.30 ± 0.16</td><td>19.23 ± 0.11</td><td>17.97 ± 0.02</td><td>17.05 ± 0.02</td><td>16.63 ± 0.03</td><td>16.21 ± 0.03</td><td>26</td></tr><tr><td>ULAS J230256.53+121310.2</td><td>sdL0</td><td>21.74 ± 0.15</td><td>19.46 ± 0.08</td><td>18.30 ± 0.04</td><td>17.49 ± 0.03</td><td>17.05 ± 0.06</td><td>16.73 ± 0.05</td><td>26</td></tr><tr><td>ULAS J230443.30+093423.9</td><td>sdL0</td><td>21.65 ± 0.16</td><td>19.27 ± 0.08</td><td>18.21 ± 0.03</td><td>17.23 ± 0.02</td><td>16.74 ± 0.03</td><td>16.27 ± 0.03</td><td>26</td></tr></tbody></table>

Notes. <sup>a</sup> Reference: 1. Burgasser et al. (2003a); 2. Burgasser et al. (2004a); 3. Burgasser (2004b); 4. Burgasser & Kirkpatrick (2006); 5. Burgasser et al. (2007); 6. Burgasser et al. (2009); 7. Bowler et al. (2010); Burningham et al. (2010); Schmidt et al. (2010); 8. Cushing et al. (2009); 9. Gizis & Harvin (2006); 10. Kirkpatrick et al. (2010); 11. Kirkpatrick et al. (2014); 12. Kirkpatrick et al. (2016); 13. Lépine et al. (2002); 14. Luhman & Sheppard (2014); 15. Lodieu et al. (2010); 16. Lodieu et al. (2012); 17. Lodieu et al. (2017); 18. Scholz et al. (2004a); 19. Scholz et al. (2004b); 20. Schneider et al. (2016); 21. Sivaraman et al. (2009); 22. Smith et al. (2018); 23. Paper I; 24. Paper II; 25. Paper III; 26. This work. <sup>b</sup>ATLAS i, z and VHS YJHK photometry. <sup>c</sup>ATLAS i, z and VIKING YJH photometry. <sup>d</sup>VHS YJHK photometry. <sup>e</sup>ATLAS i, z photometry. <sup>f</sup>VHS JK photometry. The other NIR photometry of objects with Y-band detection are from UKIDSS. J, H, and K photometry of the other objects without UKIDSS Y-band detection are converted from 2MASS with equations (1–3). <sup>g</sup>VVV J12564163–6202039 is an sdL3 subdwarf according to a new optical-NIR spectrum we obtained with the X-shooter.

300, respectively. The spectra were reduced using standard procedures with IRAF<fnref n="1" />. The HgAr, Ne, and Xe arcs were used on the wavelength calibration of our spectra. Spectral flux calibrations were achieved with standard stars (Oke 1974, 1990; Hamuy et al. 1994; Bohlin, Colina & Finley 1995). Spectral features caused by telluric absorptions in these spectra were not corrected. Five objects were observed twice to gain additional exposure time, and for the purposes of this paper we combined the spectra of these objects from both epochs.

## 2.3 VLT spectroscopy

We have also obtained optical-NIR spectroscopy of five of our candidates, a known sdL8 subdwarf and three radial velocity standard L dwarfs using the X-shooter spectrograph (Vernet et al. 2011) on the Very Large Telescope (VLT) since 2014. Table 3 shows a summary of our VLT spectroscopic observations. These X-shooter spectra were observed in an ABBA nodding mode with a 1.2 arcsec slit providing a resolving power of 6700 in the VIS arm and 4000 in the NIR arm. These spectra were reduced with ESO Reflex in a lazy mode (Freudling et al. 2013). The flat field used for the science data was also used for flat-fielding the flux standard star spectrum (by set flat strategy to false). Telluric corrections were achieved for both VIS and NIR arms using telluric standards that were observed right after or before our targets at a very close airmass. To extract telluric spectra for correction, we fitted and normalized spectra of telluric standards stars to their continuum and removed non-telluric features (e.g. hydrogen absorption lines) with IRAF SPLOT. These X-shooter spectra were smoothed by 101 and 51 pixels in the VIS and NIR arms, respectively, for display in figures, which reduced the resolving power to about 600 and increased their signal-to-noise ratio (SNR per pixel) by 10 and 7 times.

## 3 CLASSIFICATION

In Paper I, we classified L subdwarfs into sdL, esdL, and usdL subclasses. The subtypes of L subdwarfs are based on the comparison of their optical spectra to those of L dwarfs. The subclasses of L subdwarfs are based on their metallicity sensitive spectral features in the optical and NIR.

### 3.1 Optical classification

To assign spectral type/class to our candidates, we compared their optical spectra to those of established L subdwarf standards and identified the ones that provided the closest match to the overall spectral profile and to known spectral features sensitive to low metallicity. We also compared optical spectra of some candidates to those of L dwarf standards, when we could not find a close match to available L subdwarf standards.

We used TiO and VO absorption bands in the optical region to distinguish L subdwarfs from L dwarfs, taking into account the standard/model spectral comparisons shown in fig. 10 of Paper I. This figure shows that the TiO absorption bands at around 710, 775, and 850 nm are getting stronger from the dL to sdL subclasses (from [Fe/H] = 0.0 to -0.5), and then getting gradually weaker from the sdL to usdL subclasses (from [Fe/H] = -0.5 to -2.0). The strong TiO absorption bands are a signature of early-type sdL, which is not the case for early type L dwarfs. The VO absorption band at 788–810 nm gets weaker from early-type dL to sdL subclasses and disappears in the esdL subclass ([Fe/H] ≲ -1.0). Furthermore, a characteristic difference between early-type sdL and early-type esdL subclasses can be seen in the 770–810 nm range, where early-type esdL spectra follow a sloping straight line but early-type sdL spectra follow a dipping curve due to stronger TiO and weaker VO absorption.

Fig. 1 shows six of our new L subdwarfs that have stronger metal-poor spectral features with very weak or absent VO absorption bands at around 800 nm. Five of these have been classified as esdL, with one (UL0212+06) classified as sdL (but close to the sdL/esdL border). The VO absorption band strength of these subdwarfs is similar to the synthetic spectral models with [Fe/H] ≤ -1.0. By comparison, Fig. 2 shows 20 of our new L subdwarfs whose spectral features are indicative of slightly higher metallicity than those in Fig. 1. The spectra of these objects compared well with known sdL subdwarfs or showed typical sdL spectral features (strong TiO absorption bands).

#### 3.1.1 Six L subdwarfs with strong metal-poor features

ULAS J111429.54+072809.5 (UL1114+07) was classified esdL0, since it compares well with ULAS J033350.84+001406.1 (UL0333; Lodieu et al. 2012), an esdL0 classified in Paper I.

ULAS J135216.31+312327.0 (UL1352+31), ULAS J145234.65+043738.4 (UL1452+04), and ULAS J223302.03+062030.8 (UL2233+06) were all classified esdL05, since they compare well with ULAS J124425.90+102441.9 (UL1244; Lodieu et al. 2012), an esdL0.5 classified in Paper I. We also note that UL2233+06 compared best to UL1244 at

<footnote n="1">¹IRAF is distributed by the National Optical Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., under contract with the National Science Foundation.</footnote>
<table><caption>**Table 2.** Summary of the characteristics of the spectroscopic observations made with the OSIRIS on the GTC. R500R and R300R grisms cover a wavelength range of 480–1020nm. A 0.8 arcsec slit used for all observations.</caption><thead><tr><th>ULAS Name</th><th>Proposal ID</th><th>UT date</th><th>Seeing(arcsec)</th><th>Airmass</th><th>Grism</th><th>T<sub>int</sub> (s)</th><th>SpT</th><th>Std (SpT)<sup>Ref.</sup></th></tr></thead><tbody><tr><td>J111429.54+072809.5</td><td>GTC39-12B</td><td>2013-01-17</td><td>0.67</td><td>1.24</td><td>R500R</td><td>900 x 1</td><td>esdL0</td><td>GD 140 (DA2.2)<sup>3</sup></td></tr><tr><td>J135216.31+312327.0</td><td>GTC46-14A</td><td>2014-03-12</td><td>0.80</td><td>1.15</td><td>R500R</td><td>600 x 1</td><td>esdL0.5</td><td>GD 153 (DA1.2)<sup>1</sup></td></tr><tr><td>J145234.65+043738.4</td><td>GTC39-12B</td><td>2013-01-20</td><td>1.00</td><td>1.14</td><td>R500R</td><td>900 x 1</td><td>esdL0.5</td><td>G191-B2B (DA.8)<sup>3</sup></td></tr><tr><td>J223302.03+062030.8</td><td>GTC39-12B</td><td>2012-09-02</td><td>0.70</td><td>1.11</td><td>R500R</td><td>1800 x 1</td><td>esdL0.5</td><td>G191-B2B (DA.8)<sup>3</sup></td></tr><tr><td>J231924.35+052524.5</td><td>GTC39-12B</td><td>2012-09-02</td><td>0.60</td><td>1.13</td><td>R500R</td><td>900 x 1</td><td>esdL1</td><td>G191-B2B (DA.8)<sup>3</sup></td></tr><tr><td>J231924.35+052524.5</td><td>GTC46-14A</td><td>2014-07-06</td><td>0.90</td><td>1.40</td><td>R500R</td><td>900 x 1</td><td>esdL1</td><td>GD 248 (DA5)<sup>4</sup></td></tr><tr><td>J011840.73+084424.7</td><td>GTC63-13A</td><td>2013-08-01</td><td>0.80</td><td>1.09</td><td>R500R</td><td>1200 x 3</td><td>sdL0</td><td>Ross 640 (DZA5.5)<sup>3</sup></td></tr><tr><td>J021258.08+064115.9</td><td>GTC39-12B</td><td>2012-09-01</td><td>0.70</td><td>1.08</td><td>R500R</td><td>1800 x 1</td><td>sdL0.5</td><td>Ross 640 (DZA5.5)<sup>3</sup></td></tr><tr><td>J023803.12+054526.1</td><td>GTC80-15A</td><td>2015-08-23</td><td>0.70</td><td>1.18</td><td>R300R</td><td>900 x 1</td><td>sdL0</td><td>Ross 640 (DZA5.5)<sup>3</sup></td></tr><tr><td>J075335.23+200622.4</td><td>GTC39-12B</td><td>2012-12-06</td><td>0.75</td><td>1.15</td><td>R500R</td><td>900 x 1</td><td>sdL0</td><td>G191-B2B (DA.8)<sup>3</sup></td></tr><tr><td>J082206.61+044101.8</td><td>GTC39-12B</td><td>2012-12-06</td><td>0.75</td><td>1.11</td><td>R500R</td><td>900 x 1</td><td>sdL0.5</td><td>G191-B2B (DA.8)<sup>3</sup></td></tr><tr><td>J123142.99+015045.4</td><td>GTC80-15A</td><td>2015-03-16</td><td>0.90</td><td>1.52</td><td>R300R</td><td>900 x 1</td><td>sdL0</td><td>Hilt 600 (B1)<sup>2</sup></td></tr><tr><td>J124104.75-000531.4</td><td>GTC63-13A</td><td>2013-05-08</td><td>0.70</td><td>1.68</td><td>R500R</td><td>1200 x 3</td><td>sdL0</td><td>GD 190 (DB2)</td></tr><tr><td>J125226.62+092920.1</td><td>GTC80-15A</td><td>2015-03-16</td><td>0.70</td><td>1.21</td><td>R300R</td><td>900 x 1</td><td>sdL0.5</td><td>Hilt 600 (B1)<sup>2</sup></td></tr><tr><td>J130710.22+151103.4</td><td>GTC46-14A</td><td>2014-07-26</td><td>0.50</td><td>1.80</td><td>R500R</td><td>1200 x 3</td><td>sdL8</td><td>GD 153 (DA1.2)<sup>1</sup></td></tr><tr><td>J134206.86+053724.9</td><td>GTC46-14A</td><td>2014-07-25</td><td>0.80</td><td>1.56</td><td>R500R</td><td>1200 x 2</td><td>sdL1</td><td>GD 153 (DA1.2)<sup>1</sup></td></tr><tr><td>J134423.98+280603.8</td><td>GTC63-13A</td><td>2013-04-08</td><td>0.90</td><td>1.06</td><td>R500R</td><td>1200 x 3</td><td>sdL4</td><td>GD 190 (DB2)<sup>3</sup></td></tr><tr><td>J134852.93+101611.8</td><td>GTC63-13A</td><td>2013-03-19</td><td>0.90</td><td>1.37</td><td>R500R</td><td>900 x 1</td><td>sdL0</td><td>L 1363-3 (DQ6)<sup>3</sup></td></tr><tr><td>J135359.58+011856.7</td><td>GTC46-14A</td><td>2014-07-28</td><td>0.90</td><td>1.73</td><td>R500R</td><td>900 x 3</td><td>sdL0</td><td>GD 153 (DA1.2)<sup>1</sup></td></tr><tr><td>J141203.85+121609.9</td><td>GTC39-12B</td><td>2013-01-20</td><td>1.00</td><td>1.09</td><td>R500R</td><td>1800 x 1</td><td>sdL5</td><td>G191-B2B (DA.8)<sup>3</sup></td></tr><tr><td>J141832.35+025323.0</td><td>GTC39-12B</td><td>2013-01-17</td><td>0.67</td><td>1.12</td><td>R500R</td><td>900 x 1</td><td>sdL0</td><td>GD 140 (DA2.2)<sup>3</sup></td></tr><tr><td>J141832.35+025323.0</td><td>GTC46-14A</td><td>2014-03-12</td><td>0.80</td><td>1.19</td><td>R500R</td><td>600 x 1</td><td>sdL0</td><td>GD 153 (DA1.2)<sup>1</sup></td></tr><tr><td>J144151.55+043738.5</td><td>GTC46-14A</td><td>2014-07-19</td><td>0.60</td><td>1.34</td><td>R500R</td><td>1200 x 3</td><td>sdL4</td><td>GD 153 (DA1.2)<sup>1</sup></td></tr><tr><td>J151649.84+083607.1</td><td>GTC46-14A</td><td>2014-07-25</td><td>0.70</td><td>1.36</td><td>R500R</td><td>1200 x 3</td><td>sdL5</td><td>GD 153 (DA1.2)<sup>1</sup></td></tr><tr><td>J154638.34-011213.0</td><td>GTC63-13A</td><td>2013-05-07</td><td>0.90</td><td>1.17</td><td>R500R</td><td>1200 x 3</td><td>sdL3</td><td>L 1363-3 (DQ6)<sup>3</sup></td></tr><tr><td>J223440.80+001002.6</td><td>GTC63-13A</td><td>2013-07-17</td><td>0.80</td><td>1.15</td><td>R500R</td><td>1200 x 3</td><td>sdL1</td><td>L 1363-3 (DQ6)<sup>3</sup></td></tr><tr><td>J225902.14+115602.1</td><td>GTC63-13A</td><td>2013-07-18</td><td>0.95</td><td>1.18</td><td>R500R</td><td>900 x 2</td><td>sdL0</td><td>GD 140 (DA2.2)<sup>3</sup></td></tr><tr><td>J225902.14+115602.1</td><td>GTC46-14A</td><td>2014-07-20</td><td>0.60</td><td>1.06</td><td>R500R</td><td>1200 x 2</td><td>sdL0</td><td>GD 153 (DA1.2)<sup>1</sup></td></tr><tr><td>J230256.53+121310.2</td><td>GTC63-13A</td><td>2013-07-18</td><td>0.80</td><td>1.08</td><td>R500R</td><td>1350 x 2</td><td>sdL0</td><td>GD 190 (DB2)<sup>3</sup></td></tr><tr><td>J230443.30+093423.9</td><td>GTC63-13A</td><td>2013-07-13</td><td>0.70</td><td>1.08</td><td>R500R</td><td>1350 x 2</td><td>sdL0</td><td>LDS 749B (DB4)<sup>4</sup></td></tr><tr><td>J233227.03+123452.0</td><td>GTC46-14A</td><td>2014-08-24</td><td>0.80</td><td>1.07</td><td>R500R</td><td>1200 x 3</td><td>L6p+T4p</td><td>G158-100 (DC)<sup>4</sup></td></tr><tr><td>J135122.15+141914.9</td><td>GTC63-13A</td><td>2013-05-08</td><td>0.70</td><td>1.53</td><td>R500R</td><td>900 x 1</td><td>esdM7</td><td>GD 190 (DB2)<sup>3</sup></td></tr><tr><td>J002009.35+160451.2</td><td>GTC63-13A</td><td>2013-07-31</td><td>0.80</td><td>1.69</td><td>R500R</td><td>900 x 1</td><td>sdM9</td><td>GD 190 (DB2)<sup>3</sup></td></tr><tr><td>J010756.85+100811.3</td><td>GTC63-13A</td><td>2013-08-09</td><td>0.90</td><td>1.54</td><td>R500R</td><td>900 x 2</td><td>sdM7</td><td>GD 140 (DA2.2)<sup>3</sup></td></tr><tr><td>J020628.22+020255.6</td><td>GTC39-12B</td><td>2012-08-28</td><td>0.60</td><td>1.12</td><td>R500R</td><td>900 x 1</td><td>sdM7</td><td>Ross 640 (DZA5.5)<sup>3</sup></td></tr><tr><td>J024035.36+060629.3</td><td>GTC39-12B</td><td>2012-09-01</td><td>0.70</td><td>1.09</td><td>R500R</td><td>1800 x 1</td><td>sdM7</td><td>Ross 640 (DZA5.5)<sup>3</sup></td></tr><tr><td>J143154.18-004114.3</td><td>GTC46-14A</td><td>2014-07-13</td><td>0.65</td><td>1.23</td><td>R500R</td><td>900 x 1</td><td>sdM9</td><td>Ross 640 (DZA5.5)<sup>3</sup></td></tr><tr><td>J143517.18-014713.1</td><td>GTC63-13A</td><td>2013-05-08</td><td>0.70</td><td>1.67</td><td>R500R</td><td>900 x 1</td><td>sdM9</td><td>GD 190 (DB2)<sup>3</sup></td></tr><tr><td>J001747.60+130757.1</td><td>GTC63-13A</td><td>2013-07-30</td><td>0.70</td><td>1.15</td><td>R500R</td><td>1500 x 2</td><td>galaxy</td><td>—</td></tr><tr><td>J001837.37+020015.7</td><td>GTC63-13A</td><td>2013-07-30</td><td>0.90</td><td>1.13</td><td>R500R</td><td>1200 x 3</td><td>M9</td><td>L 1363-3 (DQ6)<sup>3</sup></td></tr><tr><td>J001931.33+063111.0</td><td>GTC39-12B</td><td>2012-08-26</td><td>0.65</td><td>1.23</td><td>R500R</td><td>900 x 3</td><td>M9</td><td>G191-B2B (DA.8)<sup>3</sup></td></tr><tr><td>J004602.85+091131.2</td><td>GTC63-13A</td><td>2013-08-01</td><td>0.80</td><td>1.48</td><td>R500R</td><td>1500 x 2</td><td>M9</td><td>Ross 640 (DZA5.5)<sup>3</sup></td></tr><tr><td>J004716.65+161242.4</td><td>GTC63-13A</td><td>2013-07-31</td><td>0.70</td><td>1.60</td><td>R500R</td><td>1350 x 2</td><td>M9</td><td>GD 190 (DB2)<sup>3</sup></td></tr><tr><td>J011711.98-005213.4</td><td>GTC63-13A</td><td>2013-07-29</td><td>0.90</td><td>1.18</td><td>R500R</td><td>900 x 1</td><td>M7</td><td>GD 190 (DB2)<sup>3</sup></td></tr><tr><td>J125938.50+301500.2</td><td>GTC63-13A</td><td>2013-05-08</td><td>0.70</td><td>1.48</td><td>R500R</td><td>900 x 1</td><td>M9</td><td>GD 190 (DB2)<sup>3</sup></td></tr><tr><td>J134505.85+342441.8</td><td>GTC63-13A</td><td>2013-03-17</td><td>0.80</td><td>1.16</td><td>R500R</td><td>900 x 1</td><td>M9</td><td>L 1363-3 (DQ6)<sup>3</sup></td></tr><tr><td>J134505.85+342441.8</td><td>GTC46-14A</td><td>2014-07-13</td><td>0.60</td><td>1.20</td><td>R500R</td><td>900 x 1</td><td>M9</td><td>Ross 640 (DZA5.5)<sup>3</sup></td></tr><tr><td>J205721.89+005628.7</td><td>GTC63-13A</td><td>2013-05-07</td><td>0.70</td><td>1.33</td><td>R500R</td><td>1350 x 2</td><td>M7</td><td>GD 140 (DA2.2)<sup>3</sup></td></tr><tr><td>J214816.13+012225.1</td><td>GTC63-13A</td><td>2013-07-16</td><td>0.80</td><td>1.13</td><td>R500R</td><td>1500 x 2</td><td>M7</td><td>GD 140 (DA2.2)<sup>3</sup></td></tr><tr><td>J223123.44+010025.1</td><td>GTC63-13A</td><td>2013-07-13</td><td>0.60</td><td>1.15</td><td>R500R</td><td>900 x 1</td><td>M6</td><td>LDS 749B (DB4)<sup>4</sup></td></tr><tr><td>J223623.17+034344.5</td><td>GTC63-13A</td><td>2013-07-05</td><td>0.80</td><td>1.39</td><td>R500R</td><td>1200 x 3</td><td>M9</td><td>GD 190 (DB2)<sup>3</sup></td></tr><tr><td>J224054.61+030902.0</td><td>GTC63-13A</td><td>2013-07-05</td><td>0.60</td><td>1.11</td><td>R500R</td><td>1200 x 1</td><td>M7</td><td>GD 140 (DA2.2)<sup>3</sup></td></tr><tr><td>J224749.77+053207.9</td><td>GTC46-14A</td><td>2014-07-19</td><td>0.70</td><td>1.12</td><td>R500R</td><td>1200 x 3</td><td>L1</td><td>GD 153 (DA1.2)<sup>1</sup></td></tr><tr><td>J231949.36+044559.5</td><td>GTC63-13A</td><td>2013-07-16</td><td>0.80</td><td>1.10</td><td>R500R</td><td>1200 x 1</td><td>M7</td><td>GD 140 (DA2.2)<sup>3</sup></td></tr><tr><td>J231949.36+044559.5</td><td>GTC46-14A</td><td>2014-07-20</td><td>0.60</td><td>1.21</td><td>R500R</td><td>1200 x 2</td><td>M7</td><td>GD 153 (DA1.2)<sup>1</sup></td></tr><tr><td>J233211.22+045554.2</td><td>GTC63-13A</td><td>2013-08-08</td><td>0.80</td><td>1.18</td><td>R500R</td><td>900 x 2</td><td>M6</td><td>Ross 640 (DZA5.5)<sup>3</sup></td></tr></tbody></table>

Note. References for flux calibration standard stars are 1. Bohlin et al. (1995); 2. Hamuy et al. (1994); 3. Oke (1974); 4. Oke (1990).

600–840 nm. And UL1352+31 has slightly more flux than UL1244 at 770–810 nm, which is indicative of a slightly lower metallicity. Since UL1244 has [Fe/H] ≈ −1.5 (Paper I), UL1352+31 is close to the esdL/usdL boundary.

ULAS J021258.08+064115.9 (UL0212+06) was classified sdL1, since it compares well to the sdL1 subdwarf SDSS J133348.24+273508.8 (SD1333; Paper I). It may have a flatter plateau over 0.738–0.757 µm than SD1333, though this could be
<table><caption>**Table 3.** Summary of the characteristics of the spectroscopic observations made with X-shooter. Wavelength ranges for the VIS and NIR arms are 530–1020 and 990–2480 nm. 1.2 arcsec slits are used for both VIS and NIR arms. DENIS-P J025503.3–470049, 2MASS J08354256–0819237, and DENIS-P J144137.3–4094559 were observed as RV standards.</caption><thead><tr><th>Name</th><th>SpT</th><th>UT date</th><th>Seeing (arcsec)</th><th>Airm</th><th>T<sub>int</sub> (VIS)</th><th>T<sub>int</sub> (NIR)</th><th>Telluric (SpT)</th><th>Airm</th></tr></thead><tbody><tr><td>ULAS J024035.3+060629.3</td><td>sdM7</td><td>2015-01-17</td><td>0.97</td><td>1.42</td><td>12 × 237 s</td><td>12 × 250 s</td><td>HD 22686 (A0 V)</td><td>1.27</td></tr><tr><td>ULAS J141203.85+121609.9</td><td>sdL5</td><td>2015-02-25</td><td>0.96</td><td>1.26</td><td>12 × 287 s</td><td>12 × 300 s</td><td>HIP 76069 (B9 V)</td><td>1.33</td></tr><tr><td>ULAS J233227.03+123452.0</td><td>L6p+T4p</td><td>2015-09-10</td><td>1.03</td><td>1.28</td><td>12 × 285 s</td><td>12 × 300 s</td><td>HIP 117927 (B9V)</td><td>1.20</td></tr><tr><td>ULAS J223302.03+062030.8</td><td>esdL0.5</td><td>2015-09-11</td><td>1.20</td><td>1.24</td><td>12 × 285 s</td><td>12 × 300 s</td><td>HIP 105164 (B7 III)</td><td>1.09</td></tr><tr><td>ULAS J130710.22+151103.4</td><td>sdL8</td><td>2018-04-23</td><td>1.03</td><td>1.34</td><td>12 × 285 s</td><td>12 × 300 s</td><td>HIP 61257 (B9.5 V)</td><td>1.15</td></tr><tr><td>2MASS J06453153–6646120</td><td>sdL8</td><td>2016-02-19</td><td>1.09</td><td>1.35</td><td>8 × 290 s</td><td>8 × 300 s</td><td>HD 77281 (A3 IV)</td><td>1.26</td></tr><tr><td>DENIS-P J025503.3–470049</td><td>L9</td><td>2015-08-12</td><td>1.57</td><td>1.44</td><td>8 × 140 s</td><td>8 × 150 s</td><td>HIP 105164 (B7 III)</td><td>1.36</td></tr><tr><td>2MASS J08354256–0819237</td><td>L5</td><td>2015-10-31</td><td>0.62</td><td>1.31</td><td>8 × 140 s</td><td>8 × 150 s</td><td>HIP 38734 (B9 V)</td><td>1.33</td></tr><tr><td>DENIS-P J144137.3–094559</td><td>L0.5</td><td>2016-03-19</td><td>0.66</td><td>1.11</td><td>8 × 140 s</td><td>8 × 150 s</td><td>HIP 76836 (B9.5 V)</td><td>1.04</td></tr></tbody></table>

due to noise or telluric effects. The flux of UL0212+06 beyond 880 nm is not well calibrated due to the lack of second-order correction for the OSIRIS spectrum. Although UL0212+06 (and SD1333) are classified as sdL they have weaker 800 nm VO absorption than any other sdL dwarfs, and thus lie close to the sdL/esdL border.

ULAS J231924.35+052524.5 (UL2319+05) has a similar spectral profile to SD1333, but has weaker TiO absorption bands and slightly more flux over 770–810 nm (with a sloping straight line morphology) indicative of a lower metallicity. We therefore classified this object as esdL1.

patrick et al. 1999a) at wavelengths below 860 nm. The metal-poor sensitivity of stronger TiO absorption bands below 860 nm become insignificant in the spectra of late-type sdL subclasses. However, metal hydride (CrH, FeH) absorption for late-type sdL subdwarfs is stronger than that in late-type L dwarfs. UL1307+15 has stronger 870 nm FeH absorption than 2M1632, indicating a lower metallicity. UL1307+15 compares well with the sdL8 subdwarf 2MASS J06453153–6646120 (2M0645; Kirkpatrick et al. 2010), particularly in the red optical (including the 870 nm FeH band), and it was thus classified sdL8.

### 3.1.2 Eighteen sdL0–8 subdwarfs

Thirteen of our objects in Fig. 2 compare well with the sdL0 subdwarf SDSS J134749.74+333601.7 (SD1347; Paper I) and were thus classified sdL0. These objects are ULAS J011840.73+084424.7, ULAS J023803.12+054526.1, ULAS J075335.23+200622.4, ULAS J082206.61+044101.8, ULAS J123142.99+015045.4, ULAS J124104.75 – 000531.4, ULAS J125226.62+092920.1, ULAS J134852.93+101611.8, ULAS J135359.58+011856.7, ULAS J141832.35+025323.0, ULAS J225902.14+115602.1, ULAS J230256.53+121310.2, and ULAS J230443.30+093423.9.

ULAS J134206.86+053724.9 (UL1342+05) was classified sdL0.5, since it compares well with 2MASS J00412179+3547133 (2M0041; Burgasser et al. 2004a), an sdL0.5 classified in Paper I.

ULAS J223440.80+001002.6 (UL2234+00) was classified sdL1, since it compares well to the sdL1 subdwarf 2MASS J17561080+2815238 (2M1756; Kirkpatrick et al. 2010), and shows stronger TiO absorption at around 710, 780, and 850 nm than the L1 dwarf standard 2MASS J14392836+1929149 (2M1439; Kirkpatrick et al. 1999a).

ULAS J134423.98+280603.8 (UL1344+28) and ULAS J144151.55+043738.5 (UL1441+04) compared well to the L4 standard 2MASS J1155009+230706 (2M1155; Kirkpatrick et al. 1999a). However, both have somewhat stronger TiO absorption at 850 nm (than this L4 dwarf) that is a signature of the sdL subclass. UL1344+28 and UL1441+04 compared well with the sdL4 subdwarf ULAS J021642.97+004005.6 (UL0216; Paper I) and were thus classified sdL4. We note that UL1344+28 has slightly weaker 850 nm TiO absorption than UL0216, and thus a slightly higher metallicity than UL0216. UL1441+04 was previously reported as a blue L1 dwarf by an independent search for L and T dwarfs (Marocco et al. 2015).

ULAS J130710.22+151103.4 (UL1307+15) compares well to the L8 dwarf standard 2MASS J16322911+1904407 (2M1632; Kirkpatrick et al. 1999a).

### 3.1.3 Three sdL subdwarfs with new subtypes

Three objects did not closely compare with any of the L subdwarf or L dwarf standards, but showed clear evidence of L subdwarf nature.

ULAS J154638.34–011213.0 (UL1546–01) compares reasonably to the optical spectral profile of the L3 dwarf standard 2MASS J03554191+2257016 (2M0355; Kirkpatrick et al. 1999a), but has stronger TiO absorption bands at around 710 and 850 nm. Therefore, we classified this object as sdL3.

ULAS J141203.85+121609.9 (UL1412+12) and ULAS J151649.84+083607.1 (UL1516+08) both compare reasonably to the L5 dwarf standard 2MASS J15074769–1627386 (2M1507; Reid et al. 2000), but have stronger TiO absorption at around 850 nm. We therefore classified these objects as sdL5. UL1516+08 has a slightly weaker 850 nm TiO band compared to UL1412+12, and likely has somewhat higher metallicity. These objects were both previously reported as blue L4 and L5 dwarfs, respectively, by Marocco et al. (2015).

We propose UL1412+12 as an sdL5 standard, since we have obtained good-quality optical and NIR spectra with the GTC and VLT (Section 3.2). Note that a previously classified sdL5 object, 2MASS J06164006–6407194 (2M0616; Cushing et al. 2009), was re-classified as esdL6 in Paper I. Also, ULAS J135058.85+081506.8 (Lodieu et al. 2010) was originally classified as sdL5 but was re-classified as sdL3.5–4 in Lodieu et al. (2017) and as usdL3 in Papers I and III. And a suspected sdL5 subdwarf, WISEA J135501.90–825838.9 (WI1355; Kirkpatrick et al. 2010), seems less secure as a spectroscopic standard, since its *J* − *H* and *J* − *K* colours are similar to those of L5 dwarfs.

From Table 1, we can see that we still missing some spectral subtypes of L subdwarfs. These missing subtypes including sdL2, sdL6, sdL9; esdL2, esdL5, esdL8–9; and usdL2 and usdL5–9. Most of them are late L subtypes as late L subdwarfs are fainter and more difficult to identify with current facilities.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Optical spectra of 6 new L subdwarfs (black) showing strong optical metal-poor features in comparison to the L0 dwarf (magenta) 2MASP J0345432+254023 (2M0435; Kirkpatrick et al. 1999a), and three known L subdwarfs (red): UL0333, UL1244 (Lodieu et al. 2012), and SD1333 (Paper I). Spectra are normalized near 800 nm. Spectra of UL0333, UL1244, and SD1333 were shifted by -150 to -300 km s⁻¹ to align their Na I lines with those of our objects (except for UL0212+06). Telluric absorption is indicated with shaded grey bands but not corrected. Lighter and thicker shaded bands indicate regions with weaker and stronger telluric effects.</figcaption></figure>

## 3.2 NIR constraints

L dwarfs/subdwarfs have rather diverse NIR spectral features because they are distributed into a wide range of metallicity and gravity. Objects that are more metal-poor have lower opacity, thus have more flux at shorter wavelength, and tend to be more compact (i.e. higher gravity). NIR spectral fluxes are suppressed gradually as metallicity decreases due to enhanced CIA H₂ (see fig. 9 in Paper I). The optical-NIR morphology of L subdwarfs is thus very characteristic. And within the NIR itself, it is relatively easier to distinguish sdL from dL subclass using NIR spectral features (which are more significant than optical features such as TiO absorptions for sdL subclasses). NIR spectra are also very useful for distinguishing usdL from esdL subclasses because the spectral variations caused by the metallicity differences are larger in the NIR than in the optical.

UL2233+06 was observed under our X-shooter follow-up programme. The spectrum of this object has SNR of ~8 at 820 nm and ~6 at 1300 nm, and was smoothed by 101 (VIS) and 51 (NIR) pixels for display. Fig. 3 shows the optical to NIR spectrum of UL2233+06 compared (via optical nomalization) to those of the sdL0.5 subdwarf 2MASS J00412179+3547133 (2M0041; Burgasser et al. 2004a; Paper I) and the L0 dwarf 2MASP J0345432+254023 (SD0345; Kirkpatrick et al. 1999a). UL2233+06 has much stronger NIR flux suppression than the sdL0.5 subdwarf, and its CO absorption is absent. This is entirely consistent with the esdL0.5 classification of UL2233+06, and the dL-sdL-esdL sequence clearly shows the gradually increasing level of NIR suppression.

Note that Fig. 3 does not show the true relative flux of these subclasses. Fig. 16 in Paper I shows that esdL and dL subclasses with close spectral type have similar H-band absolute magnitudes. Thus relative flux is best shown by normalising spectra in the H band, as in Fig. 4. This figure shows the same three objects from Fig. 3 but with the additional L0γ dwarf 2MASS J01415823-4633574 (2M0141; Kirkpatrick et al. 2006). Although these objects all have spectral subtype L0–L0.5 (classified in optical) they have very different physical parameters. UL2233+06 is older with lower metallicity, but relatively higher mass and $T_{\text{eff}}$. While at the other extreme the L0γ dwarf is young with higher metallicity, cooler $T_{\text{eff}}$, and much lower mass and gravity (Filippazzo et al. 2015; Faherty et al. 2016).

UL1412+12 was also observed by our X-shooter follow-up programme. The original spectrum has SNR of ~6 at 910 nm and ~30 at 1300 nm, and was smoothed by 101 (VIS) and 51 (NIR) pixels for display. Fig. 5 shows that UL1412+12 has largely suppressed NIR flux and a weaker CO absorption band compared to the L5 dwarf standard SDSS J083506.16+195304.4 (SD0835; Chiu et al. 2006; Kirkpatrick et al. 2010). This is as expected for the sdL5 classification (from the optical).

We obtained the optical to NIR spectra of UL1307+15 and 2M0645 under our X-shooter follow-up programme. The original spectrum of UL1307+15 has SNR of ~4 at 900 nm and ~7 at 1300 nm. The original spectrum of 2M0645 has SNR of ~20 at 900 nm and ~45 at 1300 nm. Fig. 6 shows that UL1307+15 has a similar spectral profile to 2M0645 from optical to H band. We note that UL1307+15 has more K-band flux and slightly narrower H-band spectral profile than 2M0645. UL1307+15 is worth of further investigation (e.g. on binarity).

To confirm the subdwarf status of the other 19 L subdwarfs (Fig. 2) without NIR spectra, we instead compared their optical to NIR photometric spectral energy distributions (SEDs) to those of known L dwarfs and subdwarfs with the same optical spectral subtypes (Fig. 7). These 19 L subdwarfs all show significantly suppressed NIR SEDs compared to those of the L dwarfs.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Optical spectra of 21 new L subdwarfs (black) compared to known L subdwarfs (red) and dwarfs (magenta). Spectra are normalized near 840 nm. Spectra of SD1347, 2M0041, and UL0216 are from Paper I. Spectra of 2M1439, 2M0355, 2M1155, and 2M1632 are from Kirkpatrick et al. (1999a). The spectrum of 2M1756 is from Kirkpatrick et al. (2010) and was corrected for telluric absorption. The spectrum of 2M1507 is from Reid et al. (2000). 2M0645 was observed by our X-shooter follow-up programme.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 3.** The X-shooter spectrum of UL2233+06 compared with optical to NIR spectra of 2M0041 (sdL0.5) and SD0345 (L0). NIR and optical spectra of 2M0041 are from Burgasser et al. (2004a) and Paper I, respectively. The spectrum of SD0345 is from Burgasser & McElwain (2006). Telluric absorption is indicated with grey shading and has been corrected in our VLT spectra. Lighter and thicker shaded bands indicate regions with weaker and stronger telluric effects.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 4.** NIR spectra of UL2233+06, 2M0041 and 2M0345 normalized in the H band. The young L0γ dwarf 2M0141 is also plotted for comparison.</figcaption></figure>

The 13 new L0 subdwarfs have similar SEDs to LSR J18261.1+301419.1 (LSR1826; Lépine et al. 2002), which was classified as sdL0 in Paper I. UL1342+05 shows a similar SED to the sdL0.5 2M0041. Although there is not sdL3 subdwarf to compare with UL1546−01, it can be seen that the UL1546−01 SED is more suppressed relative to the L3 dwarfs, and by a level similar to the suppression amongst the sdL0s (when compared to the L0 dwarfs). UL1344+28 and UL1441+04 have similar SEDs to the sdL4 Ul0216. UL1412+12 has a slightly stronger NIR suppression than UL1516+08 indicating a slightly lower metallicity. And this is also supported by its slightly stronger 850 nm TiO absorption band when compared to UL1516+08.

<figure><img src="image_8.png" /><figcaption>**Figure 5.** The optical–NIR spectrum of UL1412+12 compared to that of an L5 standard (SD0835) from Chiu et al. (2006).</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Optical–NIR spectrum of UL1307+15 compared to the L8 2M1632 and the sdL8 2M0645. The X-shooter spectra of UL1307+15 and 2M0645 were smoothed by 101 (VIS) and 51 (NIR) pixels. The spectrum of 2M1632 is from Burgasser (2007).</figcaption></figure>

## 3.3 Astrometry and radial velocity

One of our new L subdwarfs (UL0753+20) was observed by the Gaia (Gaia Collaboration et al. 2018) astrometric survey. Eleven of our new L subdwarfs (including UL0753+20) were observed by the UKIDSS second epoch survey, and are therefore included in the UKIDSS proper motion catalogue of Smith et al. (2014). Proper motions of the rest of our objects were measured from SDSS, UKIDSS, and PS1 epochs following the procedure described in Zhang et al. (2009). We estimated spectroscopic distances for these 27 new L subdwarfs based on the relationship between spectral type and H-band absolute magnitude, which is less sensitive to metallicity/subclass for L subdwarfs. The relationship was derived from 20 L subdwarfs observed by the Gaia mission (Gaia Collaboration et al. 2016, see Section 5). Table 4 shows the proper motions, distances and tangential velocities of these 27 new L subdwarfs. They are distributed at distances from ~50 to 250pc. Seventeen of them
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Optical–NIR photometric spectral energy distributions of 19 new L subdwarfs (black lines and blue diamonds) and an L1 dwarf (UL2247+05; magenta line and cyan hexagons), compared to those of L dwarfs (grey lines and green crosses) and known L subdwarfs (red lines and magenta open circles). The spectrum of LSR1826 is from Burgasser et al. (2004a). The spectrum of 2M1756 is from Kirkpatrick et al. (2010). The spectrum of UL0216 is from Paper I. Spex spectra of L0–L5 dwarfs used for comparison are of 2MASS J13313310+3407583 (L0), 2MASS J01340281+0508125 (L0; Kirkpatrick et al. 2010); 2MASSI J2107316–030733 (L0), SDSS J202820.32+005226.5 (L3), 2MASSI J1104012+195921 (L4; Burgasser et al. 2004a); 2MASSW J0228110+253738 (L0), 2MASS J12212770+0257198 (L0), 2MASSW J0208183+254253 (L1), SDSS J104842.84+011158.5 (L1; Burgasser et al. 2008a); 2MASS J20343769+0827009 (L1), 2MASS J08234818+2428577 (L3), 2MASSW J1146345+223053 (L3), 2MASS J22425317+2542573 (L3), 2MASS J13571237+1428398 (L4), 2MASS J13571237+1428398 (L4), 2MASS J17461199+5034036 (L5; Burgasser et al. 2010); SDSSp J053951.99–005902.0 [L5; M. Cushing (unpublished); Fan et al. 2000]; SD0835 (L5; Chiu et al. 2006).</figcaption></figure>
<table><caption>**Table 4.** Proper motions, distances, and tangential velocities of our 27 new L subdwarfs.</caption><thead><tr><th>Name</th><th>μ<sub>RA</sub><br>(mas yr<sup>-1</sup>)</th><th>μ<sub>Dec</sub><br>(mas yr<sup>-1</sup>)</th><th>Distance<br>(pc)</th><th>V<sub>tan</sub><br>(km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td>UL0118+03</td><td>14 ± 12</td><td>-80 ± 12</td><td>209<sup>+20</sup><sub>-18</sub></td><td>80<sup>+18</sup><sub>-18</sub></td></tr><tr><td>UL0212+06</td><td>2 ± 4</td><td>-429 ± 5</td><td>139<sup>+13</sup><sub>-12</sub></td><td>282<sup>+27</sup><sub>-25</sub></td></tr><tr><td>UL0238+05</td><td>186 ± 18</td><td>-145 ± 15</td><td>95<sup>+9</sup><sub>-8</sub></td><td>106<sup>+15</sup><sub>-14</sub></td></tr><tr><td>UL0753+20</td><td>-36 ± 1</td><td>-191 ± 1</td><td>74<sup>+7</sup><sub>-6</sub></td><td>68<sup>+72</sup><sub>-62</sub>pt</td></tr><tr><td>UL0822+04</td><td>35 ± 7</td><td>-154 ± 5</td><td>92<sup>+9</sup><sub>-8</sub></td><td>69<sup>+8</sup><sub>-7</sub></td></tr><tr><td>UL1114+07</td><td>-17 ± 9</td><td>-306 ± 6</td><td>211<sup>+20</sup><sub>-18</sub></td><td>306<sup>+31</sup><sub>-29</sub></td></tr><tr><td>UL1231+01</td><td>-225 ± 6</td><td>22 ± 4</td><td>164<sup>+16</sup><sub>-14</sub></td><td>176<sup>+18</sup><sub>-16</sub></td></tr><tr><td>UL1241-00</td><td>-61 ± 9</td><td>-45 ± 5</td><td>245<sup>+23</sup><sub>-21</sub></td><td>89<sup>+15</sup><sub>-14</sub></td></tr><tr><td>UL1252+09</td><td>-299 ± 7</td><td>6 ± 6</td><td>117<sup>+11</sup><sub>-10</sub></td><td>166<sup>+17</sup><sub>-15</sub></td></tr><tr><td>UL1307+15</td><td>-391 ± 14</td><td>-124 ± 13</td><td>86<sup>+8</sup><sub>-8</sub></td><td>168<sup>+8</sup><sub>-17</sub></td></tr><tr><td>UL1342+05</td><td>-106 ± 8</td><td>-212 ± 5</td><td>142<sup>+14</sup><sub>-12</sub></td><td>160<sup>+17</sup><sub>-15</sub></td></tr><tr><td>UL1344+28</td><td>-345 ± 11</td><td>120 ± 10</td><td>88<sup>+8</sup><sub>-8</sub></td><td>153<sup>+16</sup><sub>-15</sub></td></tr><tr><td>UL1348+10</td><td>-283 ± 11</td><td>-179 ± 9</td><td>109<sup>+10</sup><sub>-10</sub></td><td>174<sup>+18</sup><sub>-17</sub></td></tr><tr><td>UL1352+31</td><td>542 ± 8</td><td>-111 ± 29</td><td>148<sup>+14</sup><sub>-13</sub></td><td>388<sup>+42</sup><sub>-40</sub></td></tr><tr><td>UL1353+01</td><td>-93 ± 12</td><td>-37 ± 10</td><td>142<sup>+14</sup><sub>-12</sub></td><td>67<sup>+12</sup><sub>-12</sub></td></tr><tr><td>UL1412+12</td><td>-227 ± 10</td><td>-3 ± 10</td><td>54<sup>+5</sup><sub>-5</sub></td><td>58<sup>+7</sup><sub>-6</sub></td></tr><tr><td>UL1418+02</td><td>-3 ± 7</td><td>-277 ± 6</td><td>79<sup>+8</sup><sub>-7</sub></td><td>103<sup>+10</sup><sub>-10</sub></td></tr><tr><td>UL1441+04</td><td>-264 ± 8</td><td>-508 ± 8</td><td>98<sup>+9</sup><sub>-9</sub></td><td>266<sup>+26</sup><sub>-24</sub></td></tr><tr><td>UL1452+04</td><td>-46 ± 9</td><td>-291 ± 11</td><td>157<sup>+15</sup><sub>-14</sub></td><td>220<sup>+23</sup><sub>-22</sub></td></tr><tr><td>UL1516+08</td><td>-171 ± 8</td><td>44 ± 5</td><td>79<sup>+8</sup><sub>-7</sub></td><td>66<sup>+7</sup><sub>-7</sub></td></tr><tr><td>UL1546-01</td><td>-45 ± 7</td><td>-115 ± 6</td><td>122<sup>+12</sup><sub>-11</sub></td><td>71<sup>+9</sup><sub>-8</sub></td></tr><tr><td>UL2233+06</td><td>320 ± 6</td><td>-80 ± 5</td><td>236<sup>+22</sup><sub>-20</sub></td><td>370<sup>+36</sup><sub>-33</sub></td></tr><tr><td>UL2234+00</td><td>-34 ± 9</td><td>-96 ± 7</td><td>145<sup>+14</sup><sub>-13</sub></td><td>70<sup>+11</sup><sub>-10</sub></td></tr><tr><td>UL2259+11</td><td>34 ± 19</td><td>-115 ± 13</td><td>126<sup>+12</sup><sub>-11</sub></td><td>72<sup>+15</sup><sub>-15</sub></td></tr><tr><td>UL2302+12</td><td>-123 ± 7</td><td>-191 ± 6</td><td>153<sup>+15</sup><sub>-13</sub></td><td>165<sup>+17</sup><sub>-16</sub></td></tr><tr><td>UL2304+09</td><td>-21 ± 13</td><td>-46 ± 4</td><td>133<sup>+13</sup><sub>-12</sub></td><td>32<sup>+9</sup><sub>-9</sub></td></tr><tr><td>UL2319+05</td><td>509 ± 11</td><td>-401 ± 11</td><td>167<sup>+16</sup><sub>-14</sub></td><td>513<sup>+50</sup><sub>-46</sub></td></tr></tbody></table>

have tangential velocities higher than 100 km s⁻¹. Seven of them have tangential velocities above 200 km s⁻¹, which includes these six L subdwarfs with the strongest metal-poor features shown in Fig. 1.

We note that UL0212 (the sdL in Fig. 1 that is close to the sdL/esdL border) has tangential velocity of 329<sup>+71</sup><sub>-59</sub> km s⁻¹, which is similar to the space motions of the five esdL subdwarfs in Fig. 1. Thus, the most metal-poor objects of the sdL subclass (UL0212 and SD1333) both have typical halo kinematics.

We measured the radial velocities (RV) of UL2332+06, UL1412+12, UL1307+15, and 2M0645 based on their X-shooter spectra. To measure the RVs of L subdwarfs in our X-shooter follow-up programme, we also observed three RV standard L dwarfs [DENIS-P J144137.3-094559 (Martín et al. 1999), an L0.5 dwarf with an RV of -27.9 ± 1.2 km s⁻¹ (Bailer-Jones 2004); 2MASS J08354256-0819237 (Cruz et al. 2003), an L5 dwarf with an RV of 29.89 ± 0.06 km s⁻¹ (Blake, Charbonneau & White 2010); and DENIS-P J025503.3-470049 (Martín et al. 1999), an L9 dwarf with an RV of 17.5 ± 2.8 km s⁻¹ (Zapatero Osorio et al. 2007)]. All X-shooter spectra were smoothed by 21 and 11 pixels in the VIS and NIR, respectively, to increase their SNR. We first measured the RV differences between our objects and the RV standards with the closest spectral subtype using cross-correlation on their strong absorption lines (Na I and Cs I in the optical and K I in the

NIR). Measured RV differences were then corrected for barycentric effects and converted into final RV using the known RVs of the standards. The final RVs of UL2332+06, UL1412+12, UL1307+15, and 2M0645 are -164 ± 15, -57 ± 12, -60 ± 18, and -33 ± 10 km s⁻¹, respectively.

The esdL+usdL and the sdL subclasses are approximately kinematically associated with the halo population and thick disc populations, respectively. The ratio between the number of esdL+usdL and sdL subdwarfs in our UKIDSS-SDSS sample is much higher than that between halo and thick disc stars measured by Reddy et al. (2006). This is likely due to our selection bias, as esdL or usdL subdwarfs can more easily be picked out by their extreme colours. There are likely a lot more sdL subdwarfs observed by existing photometric surveys, but not yet confirmed by spectroscopy. This arises because sdL subdwarf identification is more problematic than esdL and usdL, due to contamination from scattered M and L dwarfs in colour-colour diagrams. There are 66 L subdwarfs known to-date, including 7 usdL, 20 esdL, and 39 sdL subdwarfs. Meanwhile, there are a few thousand photometric candidates (e.g. Skrzypek, Warren & Faherty 2016) and spectroscopically confirmed L dwarfs to-date (e.g. Best et al. 2018). The ratio between known esdL+usdL subdwarfs and L dwarfs is comparable to the fraction of halo population in the solar neighbourhood (e.g. 0.6 per cent; Reddy et al. 2006).

## 3.4 A blue BD binary

ULAS J233227.03+123452.0 (UL2332+12) was selected as one of our ultra-cool subdwarf candidates. This object has previously been identified by Skrzypek et al. (2016), who flagged it as a candidate binary from SED fitting. Their subsequent spectroscopic assessment classified it as a T1.5p, and following a more detailed NIR spectroscopic analysis they concluded it was more likely to be a metal-poor T0 dwarf. We re-assess the nature of this object based on our full optical-NIR spectroscopy.

Fig. 8 shows that the optical spectrum of UL2332+12 compares well with the L7 dwarf standard 2MASS J08503593+1057156 (2M0850; Kirkpatrick et al. 1999a) below 860 nm, but has relatively more flux at longer wavelength. In the red optical UL2332+12 compares better to the sdL7 subdwarf SD1416. However, UL2332+12 and SD1416 are quite different in the NIR.

UL2332+12 was observed by our X-shooter follow-up programme, and Fig. 9 shows the full optical-NIR spectrum. The original spectrum of UL2332+12 has SNR of ~20 at 900 nm and ~45 at 1300 nm, and was smoothed by 101 (VIS) and 51 (NIR) pixels for display. UL2332+12 has a similar NIR spectral profile to the sdL7 subdwarf SD1416, but has deeper H₂O absorption bands (a feature characteristic of later T dwarfs). UL2332+12 has a similar NIR spectrum to the T2 standard SDSS J075840.33+324723.4 (SD0758; Knapp et al. 2004); however, it has more flux at around 850 nm and no pronounced CH₄ absorption. Indeed, the full optical-NIR spectrum of UL2332+12 does not compare very well to any single L dwarf/subdwarf or T dwarf, but does have more in common with the spectrum of DENIS-P J225210.73-173013.4 (DE2252; Kendall et al. 2004; Reid et al. 2006), an L4+T3.5 spectroscopic binary. L dwarfs with unresolved T dwarf companions can have bluer NIR colours (e.g. fig. 14 of Zhang et al. 2010) and peculiar spectral features compared to normal L dwarfs (Burgasser et al. 2010; Bardalez Gagliuffi et al. 2014; Marocco et al. 2015; Manjavacas et al. 2016).

We compared the NIR spectrum of UL2332+12 to those of synthesized spectral binaries constructed using L and T spectra from the Spex library, as described in Burgasser et al. (2010). Fig. 10 shows
<figure><img src="image_11.png" /><figcaption>**Figure 8.** Optical spectra of UL2332+12 observed with GTC and VLT compared to those of SD1416 (from SDSS) and 2M0850 (from Kirkpatrick et al. 1999a). Spectra are normalized at 910 nm.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 9.** Optical–NIR spectrum of UL2332+12 compared to those of SD1416 (from Schmidt et al. 2010), SD0758 (from Burgasser et al. 2008c), and DE2252 (from Reid et al. 2006). Spectra are normalized at 1300 nm.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Figure 10.** The closest matching single object to UL2332+12 (top panel), and the best-fitting synthesized spectroscopic binary combination (bottom panel). The X-shooter spectrum of UL2332+12 (black) were smoothed by 401 (VIS) and 201 (NIR) pixels. The red spectrum in the bottom panel is synthesized from the L6 (blue) and T4 (magenta) dwarfs. Telluric absorption regions are excluded in our fits (indicated with grey bands). The spectrum of 2MASS J11582077+0435014 is from Kirkpatrick et al. (2010). The spectrum of SDSS J133148.92–011651.4 is from Burgasser et al. (2010). The spectrum of 2MASSI J2254188+312349 (2M2254) is from Burgasser et al. (2004a).</figcaption></figure>

best-fitting synthesized spectrum, which combines the L dwarf SDSS J133148.92–011651.4 (SD1331; Hawley et al. 2002) with the T4 dwarf (2MASSI J2254188+312349; Burgasser et al. 2002). SD1331 is optically classified L6 (Hawley et al. 2002) but also has an NIR classification of L8p (Knapp et al. 2004), a peculiar blue L dwarf. This synthesized combination is a much better match than the best-fitting single object (2MASS J11582077+0435014.sdL7; Kirkpatrick et al. 2010). We assessed our 78 best-fitting synthesized spectra to constrain spectral type estimates for primary (L6.1 ± 0.5) and the secondary (T4.0 ± 1.7) components. And we note that these best-fitting combinations all include SD1331 (or other blue L dwarfs) as the primary.

Fig. 11 compares the optical–NIR spectrum of SD1331 to that of other optically classified L6 dwarfs, demonstrating the ‘blue’ peculiarity of this object (see also Marocco et al. 2015). By direct comparison UL2332+12 has bluer *J* − *K* and redder *i* − *J* colours than SD1331, consistent with an additional unresolved mid-T companion. In colour-space SD1331 lies ~mid-way between the dL and sdL sequences (e.g. it is indicated by a magenta open pentagon in Fig. 17, a figure that will be fully discussed in Section 4.3), and is on the BD side of the stellar-substellar boundary. Objects like SD1331 are thus likely to be mildly metal-poor BDs, and our best-fitting explanation for UL2332+12 is that it is a mildly metal-poor L6+T4 blue BD binary. Note that there is not mildly metal-poor mid-T dwarf in the Spex library available for our synthesized spectroscopic fitting. However, this does not affect our conclusion, as the effect of slightly lower metallicity on the shape and flux ratio of *J* and *H* band spectra of T dwarfs is negligible (e.g. Pinfield
<figure><img src="image_14.png" /><figcaption>**Figure 11.** NIR spectrum of SD1331 compared to those of the L6 dwarfs from Burgasser et al. (2010, 2MASS J22120703+3430351; 2MASS J10365305–3441380; 2MASS J15150083+4847416) and Reid et al. (2006, 2MASS J10101480–0406499).</figcaption></figure>

et al. 2012). UL2332+12 has slightly less *K*-band flux than the synthesized spectrum. This indicates that the *T* type companion of UL2332+12 is slightly more metal poor than the *T*4 dwarf used in the synthesized spectrum.

The *H*-band absolute magnitudes (MKO) of L6 and *T*4 dwarfs are around 12.90 and 14.10 mag according to Dupuy & Liu (2012). The combined *H*-band absolute magnitude of an L6+*T*4 binary is thus around 12.59 mag, suggesting a distance of $58_{-10}^{+12}$ pc for UL2332+12 as a binary. A mildly metal-poor L6 dwarf could have a slightly brighter (~0.2 mag) *H*-band absolute magnitude compared to a normal L6 dwarf, which would lead to a slightly increased distance. However, such an increase (~5 pc) is much less than our estimated distance uncertainty and is thus not significant.

We measured the proper motion and RV of UL2332+12 following the procedure described in Section 3.3. UL2332+12 has proper motion $\mu_{RA} = 420 \pm 10$ mas yr$^{-1}$ and $\mu_{Dec} = -107 \pm 13$ mas yr$^{-1}$, which is measured from its UKIDSS and PS1 epochs with a baseline of 2.74 yr. The tangential velocity of UL2332+12 is thus $118_{-21}^{+25}$ km s$^{-1}$. The RV of UL2332+12 was measured from its X-shooter spectrum and an L9 type RV standard (DENIS-P J025503.3–470049), yielding an RV of $-35 \pm 8$ km s$^{-1}$. The space motion of UL2332+12 is thus $U = -85 \pm 18$, $V = -86 \pm 14$, and $W = -28 \pm 13$ km s$^{-1}$. We calculated the halo membership probabilities of UL2332+12 based on kinematics and stellar population fractions of the thin disc (0.93), thick disc (0.07), and halo (0.006) in the solar neighbourhood (Reddy et al. 2006). The thin disc and thick disc membership probabilities of UL2332+12 are 22 and 77 per cent, respectively.

A summary of the measured and best-fitting properties of UL2332+12 are presented in Table 5.

## 3.5 M subdwarfs and M–L dwarfs amongst our candidates

About one-third (22) of our candidate L subdwarfs were identified as M6–L1 dwarfs or M7–M9 subdwarfs through spectroscopic follow-up. Table 6 shows spectral types and photometry for these objects.

ULAS J224749.77+053207.9 (UL2247+05) has *i* − *J* and *J* − *K* colours similar to those of mid-type sdL. However, its optical spectrum compares well with the L1 dwarf 2M1349 (see Fig. 12).

<table><caption>**Table 5.** Properties of a spectral binary UL2332+12. Note the spectral type is from synthesized fit.</caption><thead><tr><th>Parameter</th><th>UL2332+12</th></tr></thead><tbody><tr><td>α (J2000)</td><td>23<sup>h</sup>32<sup>m</sup>27<sup>s</sup>03</td></tr><tr><td>δ (J2000)</td><td>+12°34′52.′′0</td></tr><tr><td>Epoch</td><td>2009-09-04</td></tr><tr><td>SDSS i</td><td>22.64 ± 0.27</td></tr><tr><td>SDSS z</td><td>19.91 ± 0.11</td></tr><tr><td>PS1 i</td><td>22.00 ± 0.17</td></tr><tr><td>PS1 z</td><td>20.23 ± 0.04</td></tr><tr><td>PS1 y</td><td>19.13 ± 0.03</td></tr><tr><td>UKIDSS Y</td><td>18.10 ± 0.04</td></tr><tr><td>UKIDSS J</td><td>16.90 ± 0.02</td></tr><tr><td>UKIDSS H</td><td>16.39 ± 0.03</td></tr><tr><td>UKIDSS K</td><td>15.88 ± 0.03</td></tr><tr><td>WISEW1</td><td>15.18 ± 0.04</td></tr><tr><td>WISEW2</td><td>14.77 ± 0.07</td></tr><tr><td>Spectral type</td><td>L6p+T4p</td></tr><tr><td>Distance (pc)</td><td>58<sup>+12</sup><sub>-10</sub></td></tr><tr><td>μ<sub>RA</sub>(mas yr<sup>-1</sup>)</td><td>420 ± 10</td></tr><tr><td>μ<sub>Dec</sub>(mas yr<sup>-1</sup>)</td><td>-107 ± 13</td></tr><tr><td>V<sub>tan</sub> (km s<sup>-1</sup>)</td><td>118<sup>+25</sup><sub>-21</sub></td></tr><tr><td>RV (km s<sup>-1</sup>)</td><td>-35 ± 8</td></tr><tr><td>U (km s<sup>-1</sup>)</td><td>-85 ± 18</td></tr><tr><td>V (km s<sup>-1</sup>)</td><td>-86 ± 14</td></tr><tr><td>W (km s<sup>-1</sup>)</td><td>-28 ± 13</td></tr></tbody></table>

Its *i*<sub>P1</sub> − *J* and *J* − *K* colours are consistent with early-type L dwarfs, and we note that its SDSS *i*-band magnitude has a much larger uncertainty than the PS1 *i*<sub>P1</sub>-band magnitude. Fig. 7 shows that UL2247+05 has a slightly bluer SEDs than other L1 dwarfs. UL2247+05 has a proper motion of $\mu_{RA} = -56 \pm 14$ mas yr$^{-1}$ and $\mu_{Dec} = -56 \pm 15$ mas yr$^{-1}$ (measured using the method described in Section 3.3). Its distance is $146_{-25}^{+30}$ pc based on spectral type and *H*-band magnitude, leading to a tangential velocity of $54_{-17}^{+18}$ km s$^{-1}$.

When considering the M subdwarfs, there are several classification schemes for the late-types. In the most recent sdM scheme (Lépine et al. 2007), the metallicity consistency of the latest sdM subclasses has not been fully tested across all subtypes, due to the lack of wide binaries containing both early- and late-type M subdwarfs. We have also noticed that the strength of the metallicity features are usually under estimated for late-type M subdwarfs compared to early-type M subdwarfs (section 4.1 of Paper I). As a result, late-type M stars with the same metallicity as early-type esdM subdwarfs can be classified as sdM, and some late-type M stars with the same metallicity as early-type sdM subdwarfs can be classified as M dwarfs. Zhang et al. (2013) noticed that ~18 per cent of SDSS M dwarfs classified on the scheme of Lépine et al. (2007) have thick disc or halo like kinematics ($V < -100$ km s$^{-1}$). Furthermore, different subclass names are used for objects with similar metallicity in different classification schemes. The late-type d/sdM subclass of Burgasser et al. (2007) would be expected to have a similar metallicity range as the early-type sdM subclass of Lépine et al. (2007).

We classified our M subdwarfs using their optical metal-poor spectral features, such as weak VO absorption near 800 nm. This represents an extension of our L subdwarf scheme (Paper I) into the late-type sdMs, and allows us to consistently carry out sdM/sdL classification of our programmatic sample. Fig. 13 shows the spectra of these 21 late-type M stars compared to stacked SDSS M6–M9 spectral templates (Bochanski et al. 2007). Seven of these objects show weaker VO absorption bands than the M
<table><caption>**Table 6.** Photometric properties of new M–L dwarfs/subdwarfs identified in this work.</caption><thead><tr><th>Name</th><th>SpT</th><th>SDSS <i>i</i></th><th>SDSS <i>z</i></th><th>UKIDSS <i>Y</i></th><th>UKIDSS <i>J</i></th><th>UKIDSS <i>H</i></th><th>UKIDSS <i>K</i></th></tr></thead><tbody><tr><td>ULAS J135122.15+141914.9</td><td>esdM7</td><td>19.24 ± 0.03</td><td>18.02 ± 0.03</td><td>17.01 ± 0.01</td><td>16.31 ± 0.01</td><td>15.92 ± 0.01</td><td>15.63 ± 0.02</td></tr><tr><td>ULAS J002009.35+160451.2</td><td>sdM9</td><td>20.81 ± 0.06</td><td>19.21 ± 0.05</td><td>18.07 ± 0.03</td><td>17.32 ± 0.02</td><td>16.90 ± 0.04</td><td>16.65 ± 0.04</td></tr><tr><td>ULAS J010756.85+100811.3</td><td>sdM7</td><td>21.29 ± 0.09</td><td>20.00 ± 0.11</td><td>18.80 ± 0.06</td><td>18.08 ± 0.07</td><td>17.76 ± 0.08</td><td>17.50 ± 0.10</td></tr><tr><td>ULAS J020628.22+020255.6</td><td>sdM7</td><td>20.79 ± 0.04</td><td>19.48 ± 0.05</td><td>18.45 ± 0.04</td><td>17.65 ± 0.03</td><td>17.28 ± 0.07</td><td>17.10 ± 0.10</td></tr><tr><td>ULAS J024035.36+060629.3</td><td>sdM7</td><td>21.07 ± 0.15</td><td>19.82 ± 0.08</td><td>18.66 ± 0.05</td><td>17.99 ± 0.04</td><td>17.65 ± 0.06</td><td>17.48 ± 0.12</td></tr><tr><td>ULAS J143154.18–004114.3</td><td>sdM9</td><td>20.86 ± 0.06</td><td>19.25 ± 0.05</td><td>17.92 ± 0.03</td><td>17.19 ± 0.03</td><td>16.84 ± 0.04</td><td>16.48 ± 0.05</td></tr><tr><td>ULAS J143517.18–014713.1</td><td>sdM8</td><td>19.79 ± 0.04</td><td>18.31 ± 0.03</td><td>17.20 ± 0.02</td><td>16.39 ± 0.01</td><td>15.97 ± 0.02</td><td>15.62 ± 0.02</td></tr><tr><td>ULAS J001837.37+020015.7</td><td>M9</td><td>21.98 ± 0.16</td><td>20.19 ± 0.11</td><td>19.10 ± 0.08</td><td>17.99 ± 0.06</td><td>17.58 ± 0.07</td><td>17.15 ± 0.08</td></tr><tr><td>ULAS J001931.33+063111.0</td><td>M9</td><td>21.79 ± 0.13</td><td>19.76 ± 0.08</td><td>18.47 ± 0.05</td><td>17.53 ± 0.04</td><td>17.07 ± 0.04</td><td>16.71 ± 0.05</td></tr><tr><td>ULAS J004602.85+091131.2</td><td>M9</td><td>22.32 ± 0.26</td><td>20.28 ± 0.16</td><td>19.07 ± 0.07</td><td>18.14 ± 0.05</td><td>17.66 ± 0.08</td><td>17.29 ± 0.10</td></tr><tr><td>ULAS J004716.65+161242.4</td><td>M9</td><td>21.86 ± 0.16</td><td>20.76 ± 0.28</td><td>19.22 ± 0.08</td><td>18.44 ± 0.10</td><td>17.99 ± 0.12</td><td>18.22 ± 0.26</td></tr><tr><td>ULAS J011711.98–005213.4</td><td>M7</td><td>20.67 ± 0.04</td><td>19.45 ± 0.06</td><td>18.21 ± 0.03</td><td>17.53 ± 0.03</td><td>17.13 ± 0.05</td><td>16.92 ± 0.07</td></tr><tr><td>ULAS J125938.50+301500.2</td><td>M9</td><td>19.62 ± 0.03</td><td>17.87 ± 0.02</td><td>16.41 ± 0.01</td><td>15.57 ± 0.01</td><td>15.18 ± 0.01</td><td>14.74 ± 0.01</td></tr><tr><td>ULAS J134505.85+342441.8</td><td>M9</td><td>20.64 ± 0.04</td><td>18.91 ± 0.04</td><td>17.61 ± 0.03</td><td>16.77 ± 0.02</td><td>16.45 ± 0.02</td><td>16.04 ± 0.03</td></tr><tr><td>ULAS J205721.89+005628.7</td><td>M7</td><td>21.41 ± 0.09</td><td>19.95 ± 0.08</td><td>18.90 ± 0.08</td><td>17.98 ± 0.06</td><td>17.65 ± 0.08</td><td>17.26 ± 0.08</td></tr><tr><td>ULAS J214816.13+012225.1</td><td>M7</td><td>21.74 ± 0.13</td><td>20.30 ± 0.13</td><td>19.34 ± 0.10</td><td>18.39 ± 0.08</td><td>17.97 ± 0.11</td><td>18.07 ± 0.20</td></tr><tr><td>ULAS J223123.44+010025.1</td><td>M6</td><td>20.84 ± 0.05</td><td>19.73 ± 0.07</td><td>18.47 ± 0.05</td><td>17.53 ± 0.04</td><td>17.46 ± 0.06</td><td>17.29 ± 0.10</td></tr><tr><td>ULAS J223623.17+034344.5</td><td>M9</td><td>22.22 ± 0.16</td><td>20.25 ± 0.11</td><td>19.03 ± 0.06</td><td>18.17 ± 0.05</td><td>17.91 ± 0.06</td><td>17.51 ± 0.10</td></tr><tr><td>ULAS J224054.61+030902.0</td><td>M7</td><td>21.04 ± 0.05</td><td>19.87 ± 0.08</td><td>18.85 ± 0.07</td><td>18.02 ± 0.06</td><td>17.59 ± 0.07</td><td>17.46 ± 0.10</td></tr><tr><td>ULAS J224749.77+053207.9</td><td>L1</td><td>22.44 ± 0.23</td><td>20.12 ± 0.12</td><td>18.60 ± 0.04</td><td>17.46 ± 0.03</td><td>16.90 ± 0.04</td><td>16.43 ± 0.04</td></tr><tr><td>ULAS J231949.36+044559.5</td><td>M7</td><td>21.34 ± 0.13</td><td>19.95 ± 0.13</td><td>18.94 ± 0.06</td><td>18.16 ± 0.06</td><td>17.70 ± 0.10</td><td>17.81 ± 0.18</td></tr><tr><td>ULAS J233211.22+045554.2</td><td>M6</td><td>21.28 ± 0.17</td><td>20.04 ± 0.14</td><td>18.95 ± 0.08</td><td>18.16 ± 0.07</td><td>17.98 ± 0.08</td><td>17.66 ± 0.10</td></tr></tbody></table>

<figure><img src="image_15.png" /><figcaption>**Figure 12.** The optical spectrum of UL2247+05 compared to the L1 2M1439 (Kirkpatrick et al. 1999a).</figcaption></figure>

dwarf templates. ULAS J135122.15+141914.9 has the strongest metal-poor features with the absence of 800 nm VO absorption, and was classified as esdM7. ULAS J002009.35+160451.2, ULAS J020628.22+020255.6, and ULAS J024035.36+060629.3 (UL0240+06) have a weaker 800 nm VO absorption bands than the M dwarf templates, and were classified as sdM9, sdM7, and sdM7, respectively. ULAS J010756.85+100811.3, ULAS J143154.18–004114.3, and ULAS J143517.18–014713.1 have weaker 800 nm VO absorption bands than the M dwarf templates, but slightly stronger than UL0240+06. We classified these objects as sdM subdwarfs, although they could be classified as dM or d/sdM using existing classification schemes. The other 14 objects were classified as late-type M dwarfs since they compare well with M dwarf templates.

Fig. 14 shows the X-shooter spectrum of UL0240+06. The original spectrum of UL0240+06 has SNR of ~2 at 910 nm and ~4 at

1300 nm, and was smoothed by 101 (VIS) and 51 (NIR) pixels for display. UL0240+06 shows suppressed NIR flux compared to the M7 dwarf VB 8 (Kirkpatrick et al. 1991). UL0240+06 has a similar spectral profile to the d/sdM7 object 2MASS J13593574+3031039 (2M1359; Burgasser et al. 2004a).

## 4 L SUBDWARF PHOTOMETRY AND STELLAR/SUBSTELLAR MIX

The label ‘subdwarf’ originally acknowledged the location of these objects below the stellar main sequence in the Hertzsprung–Russell diagram (HRD), due to their lower opacity and bluer colour resulting from subsolar metallicity. However, the three metallicity subclasses of L subdwarfs exhibit a range of different characteristics, sometimes showing quite substantial differences in spectral energy distribution from their L dwarf counterparts (see figs 14 and 15 of Paper I). In this section, we explore the optical to infrared colours of L subdwarfs across the different subtypes and subclasses, compare them with L dwarf and main sequence star populations, and assess the stellar/substellar mix amongst L subdwarfs.

### 4.1 A photometric sample of L subdwarfs

To study the photometric diversity of L subdwarfs, we joined our UKIDSS-SDSS sample with L subdwarfs collected from the literature to form a larger sample of 66 L subdwarfs. We gathered their photometry: SDSS *i*, *z*, UKIDSS (MKO) *Y*, *J*, *H*, *K*. We also used *i* and *z* photometry from the VLT Survey Telescope’s (VST) ATLAS survey (Shanks et al. 2015) for a few objects that are not covered by SDSS. We used NIR photometry from the Visible and Infrared Survey Telescope for Astronomy (VISTA) surveys for a few objects that were not observed by UKIDSS. We also gathered W1, W2, and W3 photometry from the all-sky WISE data release, and *i*<sub>P1</sub>, *z*<sub>P1</sub>, and *y*<sub>P1</sub> photometry from Pan-STARRS (PS1; Chambers et al. 2016). This compilation of optical to infrared photometry of 66 L subdwarfs is presented in Tables 1 and 7.

Some known L subdwarfs that were discovered by 2MASS and WISE were not observed by UKIDSS or VISTA’s Hemisphere Sur-
<figure><img src="image_16.png" /><figcaption>Figure 13. Optical spectra of newly discovered M subdwarfs (left panel) and dwarfs (middle and right panels) compared to SDSS M dwarf template spectra (Bochanski et al. 2007). The VO absorption band at around 800 nm is highlighted with a blue band in the left-hand panel. Spectra are normalized at 825 nm.</figcaption></figure>

vey (VHS) and Kilo-Degree Infrared Galaxy Survey (VIKING). To convert 2MASS photometry ($J_2$, $H_2$, and $K_2$) into the MKO system, we determined polynomial relationships between MKO and 2MASS photometric differences and spectral type (SpT) for known L0–T0 dwarfs. Ideally, these relations would be based on L subdwarf measurements, but data availability requires us to take the next best approach, and use L dwarfs:

$$J - J_2 = -0.0230 - 0.0023 \times \text{SpT}; (0.1062), \qquad (1)$$

$$H - H_2 = 0.0188 + 0.0029 \times \text{SpT}; (0.1027), \qquad (2)$$

$$K - K_2 = -0.0977 + 0.0068 \times \text{SpT}; (0.0938), \qquad (3)$$

where SpT is 10 for L0, 15 for L5, and 20 for T0. The root mean square (rms) of these fits are given after each equation.

Some L subdwarfs were not in the SDSS footprint but were observed by PS1. We thus studied the correlation between SDSS and PS1 photometry. We found that SDSS and PS1 photometry are similar in the *i* band (<0.04 for main sequence stars), but correlate with spectral type in the *z* band. This is presumably because the SDSS and PS1 *i*-band filters have similar wavelength coverage and transmission profiles, though their *z*-band filters are different. These differences between SDSS and PS1 *z*-band magnitudes (*z* −
<figure><img src="image_17.png" /><figcaption>Figure 14. The X-shooter spectrum of UL0240+06 compared to those of 2M1359 (d/sdM7; Burgasser et al. 2004a) and VB 8 (M7; Kirkpatrick et al. 1991; Burgasser et al. 2008a).</figcaption></figure>

$z_{P1}$) are around −0.05 for FGK dwarfs, −0.11 for M0–M3 dwarfs, and −0.43 for L dwarfs. We determined polynomial relationships between SDSS–PS1 photometric differences and spectral type for known L0–T0 dwarfs:

$$i - i_{P1} = -0.0565 + 0.0043 \times \text{SpT}; (0.0944), \quad (4)$$

$$z - z_{P1} = -0.3092 - 0.0086 \times \text{SpT}; (0.0886). \quad (5)$$

Fig. 15 shows our relationships between SDSS–PS1 and MKO–2MASS photometric differences and spectral types for L0–T0 dwarfs.

Four of 66 known L subdwarfs (2M0041, WISEA J005757.65+201304.0, UL1241–00, and WI1355) are excluded from our sample due to problematic photometry or ambiguous spectral type. We thus have 62 L subdwarfs in our photometric sample.

## 4.2 L subdwarf colour–spectral type relations

To identify the best colours for distinguishing and characterising L subdwarfs with different spectral types and subclasses, we assessed various optical to near- and mid-infrared colours. Fig. 16 shows optical to infrared colours of sdL, esdL, and usdL subdwarfs compared to third-order polynomial fits of spectral type–colour relations for L dwarfs. The coefficients of these polynomial fits and rms are presented in Table 8. L dwarfs used for these polynomial fits are collected from DwarfArchives.org, and cross-matched within the SDSS, UKIDSS, PS1, and WISE data bases. Note that a few objects near the WISE detection limit have photometric uncertainties of 0.6 mag.

The $J - H$, $J - K$, $J - W2$, $y_{P1} - H$, $z_{P1} - H$, $y_{P1} - K$, $Y - W1$, and $J - W1$ are the best metallicity indicators, as they can be used to separate sdL, esdL, and usdL subclasses from L dwarfs. The $i - J$, $z - J$, $i_{P1} - y_{P1}$, $z_{P1} - y_{P1}$, $i_{P1} - J$, $y_{P1} - J$, and $i_{P1} - H$ can be used to separate esdL and usdL from L dwarfs, but cannot distinguish sdL from L dwarfs well.

The esdL subclass has bluer $H - W1$ and $H - W2$ colours compared to L dwarfs, but the usdL subclass tends to have relatively redder $H - W1$ and $H - W2$ colours. The most metal-poor usdL subdwarf, SD0104+15 (Paper II), shows a redder $H - W2$ and much

redder $H - W1$ compared to L dwarfs. The esdL and usdL subclasses are redder in $K - W1$ than the L dwarfs. Note that SD0104+15 is much redder in $K - W1$ than the L dwarfs.

The $i - z$, $i_{P1} - z_{P1}$, $H - W2$, $W1 - W2$, and $W2 - W3$ colours are not good metallicity indicators. The spectral subtypes of L subdwarfs are based on red optical spectra ($i$ and $z$ bands) by comparison to L dwarfs (Paper I). Therefore, L subdwarfs have similar $i - z$ colours to L dwarfs with the same subtype. However, early-type esdL and usdL subdwarfs have a different optical spectral profile to early-type L dwarfs, and have bluer $i - z$ and $i_{P1} - z_{P1}$ colours compared to L dwarfs. L subdwarfs also appear to be slightly redder in $W2 - W3$ than the L dwarfs.

## 4.3 L subdwarfs in colour–colour plots

Colour–colour plots are often used to distinguish celestial objects of different types in large-scale photometric surveys. Particular objects occupy a variety of colour spaces, with varying level of contamination from other populations (e.g. fig. 2 of Zhang et al. 2013).

Fig. 17 shows the $i - J$ versus $J - K$ colour–colour plot for sdL, esdL, and usdL subdwarfs compared to L dwarfs and main-sequence stars. The figure also shows the BT-Dusty model predicted colours of L subdwarfs (Allard, Homeier & Freytag 2014). The different metallicity ranges of these four subclasses are well represented by their broad-band optical-NIR colours, and the dL, sdL, esdL, and usdL sequences are quite distinct from late M through the L types. The $i - J$ and $J - K$ colours are very useful for selecting L subdwarf candidates. The $i - J$ colour is a good indicator of $T_{\text{eff}}$ and can be used to separate L subdwarfs from main sequence stars, while the $J - K$ colour is a good indicator of metallicity and can be used to separate L subdwarfs from L dwarfs.

Fig. 18 (left) shows the $z - J$ versus $J - H$ colour–colour plot for L subdwarfs and dwarfs. The main-sequence population and BT-Dusty model grid with $\log g = 5.5$ are plotted for comparison. The dL, sdL, esdL, and usdL sequences are revealed in the $z - J$ versus $J - H$ plane, although not as well as in Fig. 17. Fig. 18 (right) shows the $J - W2$ versus $J - K$ colour–colour plot showing dL, dT, sdL, esdL, and usdL populations. The L subdwarf populations are generally separated from each other by their $J - W2$ and $J - K$ colours, although they are very close to the main sequence. Late-type L subdwarfs begin to merge with the T dwarfs on the right side of the diagram.

To further explore the potential of future L subdwarf searches and characterization, Fig. 19 compares our L subdwarf photometric sample to L dwarfs and main-sequence stars in twelve optical/NIR, NIR, and NIR/mid-IR colour–colour diagrams (constructed using the PS1, MKO, and WISE photometric systems). $J - K$ or $J - H$ colour are used for 11 of these plots, and are sensitive to metallicity. As for Figs 17 and 18, subdwarfs of different subclasses in these eleven plots are generally well separated. The longer base-line optical/NIR colours (e.g. $i_{P1} - J$ and $z_{P1} - J$) are relatively better at separating objects with different spectral type. The last three colour–colour plots in Fig. 19 use the $K - W1$ colour and are very good at isolating the most metal-poor L subdwarfs (usdLs), which have red $K - W1$ colour (also see Fig. 16).

## 4.4 The BD transition-zone

The luminosity (or $T_{\text{eff}}$) is a decreasing function of metallicity above the hydrogen-burning minimum mass (HBMM), but an increasing function of metallicity below it. A transition-zone between VLMS and degenerate BDs is revealed as a waterfall-like feature on the 10
<table><caption>Table 7. PS1 and WISE photometry of L subdwarfs.</caption><thead><tr><th>Name</th><th>PS1 name</th><th>i<sub>P1</sub></th><th>z<sub>P1</sub></th><th>γ<sub>P1</sub></th><th>W1</th><th>W2</th></tr></thead><tbody><tr><td>SSSPM J10130734-1356204</td><td>J101307404-135632.796</td><td>17.23 ± 0.01</td><td>16.25 ± 0.01</td><td>15.92 ± 0.01</td><td>13.80 ± 0.03</td><td>13.60 ± 0.03</td></tr><tr><td>SDSSS J010448.447+153501.9</td><td>J010448.655+153459.325</td><td>20.52 ± 0.02</td><td>19.49 ± 0.02</td><td>19.09 ± 0.03</td><td>16.61 ± 0.08</td><td>16.36 ± 0.25</td></tr><tr><td>SDSS J125637.13-022452.4</td><td>J125636.715-022456.107</td><td>19.45 ± 0.02</td><td>18.10 ± 0.02</td><td>17.51 ± 0.01</td><td>15.21 ± 0.04</td><td>15.01 ± 0.08</td></tr><tr><td>ULAS J135058.85+081506.8</td><td>J135058.782+081505.607</td><td>21.26 ± 0.17</td><td>19.76 ± 0.02</td><td>19.32 ± 0.04</td><td>17.45 ± 0.20</td><td>≥16.45</td></tr><tr><td>2MASS J16262034+3925190</td><td>J162618.820+392522.418</td><td>17.92 ± 0.01</td><td>16.34 ± 0.01</td><td>15.79 ± 0.01</td><td>13.48 ± 0.02</td><td>13.14 ± 0.03</td></tr><tr><td>WISEA J213409.15+713236.1</td><td>J213409.144+713236.852</td><td>19.09 ± 0.01</td><td>17.99 ± 0.02</td><td>17.58 ± 0.02</td><td>15.23 ± 0.03</td><td>14.95 ± 0.04</td></tr><tr><td>ULAS J230711.01+014447.1</td><td>J230711.008+014446.289</td><td>21.69 ± 0.16</td><td>20.16 ± 0.04</td><td>19.57 ± 0.08</td><td>17.47 ± 0.18</td><td>≥16.79</td></tr><tr><td>WISEA J001450.17-083823.4</td><td>J001450.175-083823.247</td><td>17.40 ± 0.06</td><td>16.26 ± 0.01</td><td>15.77 ± 0.01</td><td>13.43 ± 0.03</td><td>13.20 ± 0.03</td></tr><tr><td>WISEA J020201.25-313645.2</td><td>J020201.198-313649.894</td><td>18.48 ± 0.01</td><td>—</td><td>—</td><td>14.31 ± 0.03</td><td>13.96 ± 0.04</td></tr><tr><td>ULAS J020858.62+020657.0</td><td>J020858.628+020656.920</td><td>21.39 ± 0.03</td><td>19.99 ± 0.03</td><td>19.50 ± 0.06</td><td>17.17 ± 0.13</td><td>16.63 ± 0.29</td></tr><tr><td>ULAS J033351.10+001405.8</td><td>J033351.342+001405.465</td><td>19.22 ± 0.01</td><td>18.05 ± 0.01</td><td>17.54 ± 0.01</td><td>15.08 ± 0.04</td><td>14.77 ± 0.07</td></tr><tr><td>WISEA J030601.66-033059.0</td><td>J030601.672-033059.676</td><td>18.02 ± 0.02</td><td>16.57 ± 0.01</td><td>15.99 ± 0.01</td><td>13.43 ± 0.03</td><td>13.18 ± 0.03</td></tr><tr><td>WISEA J043535.82+211508.9</td><td>J043535.926+211507.487</td><td>18.60 ± 0.01</td><td>17.15 ± 0.01</td><td>16.56 ± 0.01</td><td>14.00 ± 0.03</td><td>13.67 ± 0.04</td></tr><tr><td>2MASS J05325346+8246465</td><td>J053309.378+824622.210</td><td>20.31 ± 0.06</td><td>18.09 ± 0.01</td><td>16.95 ± 0.01</td><td>13.82 ± 0.03</td><td>13.26 ± 0.03</td></tr><tr><td>2MASS J06164006-6407194</td><td>—</td><td>—</td><td>—</td><td>—</td><td>15.65 ± 0.03</td><td>15.18 ± 0.04</td></tr><tr><td>ULAS J111429.54+072809.5</td><td>J111429.525+072809.033</td><td>20.69 ± 0.02</td><td>19.55 ± 0.02</td><td>19.06 ± 0.03</td><td>16.84 ± 0.12</td><td>16.85 ± 0.39</td></tr><tr><td>SDSS J124410.11+273625.8</td><td>J124410.063+273624.064</td><td>20.40 ± 0.02</td><td>19.32 ± 0.01</td><td>18.94 ± 0.03</td><td>16.68 ± 0.09</td><td>16.94 ± 0.40</td></tr><tr><td>ULAS J124425.75+102439.3</td><td>J124425.617+102436.751</td><td>19.47 ± 0.03</td><td>18.21 ± 0.01</td><td>17.71 ± 0.01</td><td>15.45 ± 0.04</td><td>15.14 ± 0.09</td></tr><tr><td>ULAS J135216.31+312327.0</td><td>J135216.314+312326.963</td><td>19.97 ± 0.03</td><td>18.81 ± 0.01</td><td>18.34 ± 0.02</td><td>16.14 ± 0.06</td><td>15.89 ± 0.13</td></tr><tr><td>SDSS J141405.74-014202.7</td><td>J141405.632-014204.911</td><td>19.88 ± 0.01</td><td>18.73 ± 0.01</td><td>18.20 ± 0.02</td><td>—</td><td>—</td></tr><tr><td>SSSPM J144420.67-201922.2</td><td>J144418.155-201946.066</td><td>16.08 ± 0.01</td><td>14.69 ± 0.01</td><td>14.08 ± 0.02</td><td>11.47 ± 0.02</td><td>11.21 ± 0.02</td></tr><tr><td>ULAS J145234.65+043738.4</td><td>J145234.651+043737.276</td><td>20.76 ± 0.04</td><td>19.45 ± 0.01</td><td>18.95 ± 0.04</td><td>16.15 ± 0.05</td><td>16.23 ± 0.18</td></tr><tr><td>ULAS J151913.03-000030.0</td><td>J151913.024-000031.791</td><td>21.36 ± 0.07</td><td>19.63 ± 0.02</td><td>18.86 ± 0.02</td><td>16.25 ± 0.06</td><td>15.72 ± 0.14</td></tr><tr><td>2MASS J16403197+1231068</td><td>J164031.759+123104.908</td><td>19.15 ± 0.01</td><td>17.96 ± 0.02</td><td>17.39 ± 0.01</td><td>15.03 ± 0.04</td><td>14.85 ± 0.07</td></tr><tr><td>WISEA J204027.30+695924.1</td><td>J204027.336+695924.279</td><td>16.89 ± 0.01</td><td>15.65 ± 0.01</td><td>15.10 ± 0.01</td><td>12.68 ± 0.02</td><td>12.46 ± 0.02</td></tr><tr><td>ULAS J223302.03+062030.8</td><td>J223302.105+062030.621</td><td>21.15 ± 0.04</td><td>19.89 ± 0.03</td><td>19.38 ± 0.02</td><td>16.92 ± 0.12</td><td>≥16.45</td></tr><tr><td>ULAS J231924.35+052524.5</td><td>J231924.396+052524.161</td><td>20.78 ± 0.06</td><td>19.51 ± 0.03</td><td>18.89 ± 0.02</td><td>16.44 ± 0.08</td><td>16.12 ± 0.19</td></tr><tr><td>2MASS J00412179+3547133</td><td>J004121.680+354712.791</td><td>20.21 ± 0.04</td><td>18.70 ± 0.02</td><td>17.82 ± 0.01</td><td>14.74 ± 0.03</td><td>14.45 ± 0.05</td></tr><tr><td>WISEA J005757.65+201304.0</td><td>—</td><td>—</td><td>—</td><td>—</td><td>14.32 ± 0.03</td><td>13.87 ± 0.04</td></tr><tr><td>WISEA J011639.05-165420.5</td><td>—</td><td>—</td><td>—</td><td>—</td><td>14.42 ± 0.03</td><td>13.99 ± 0.04</td></tr><tr><td>ULAS J011824.89+034130.4</td><td>J011824.899+034130.329</td><td>21.20 ± 0.05</td><td>20.19 ± 0.04</td><td>19.55 ± 0.10</td><td>—</td><td>—</td></tr><tr><td>WISEA J013012.66-104732.4</td><td>—</td><td>—</td><td>—</td><td>—</td><td>14.42 ± 0.03</td><td>14.00 ± 0.04</td></tr><tr><td>ULAS J021258.08+064115.9</td><td>J021258.090+064113.734</td><td>21.00 ± 0.07</td><td>19.60 ± 0.01</td><td>19.03 ± 0.02</td><td>16.31 ± 0.06</td><td>15.98 ± 0.18</td></tr><tr><td>ULAS J021642.96+004005.7</td><td>J021642.932+004005.180</td><td>21.77 ± 0.12</td><td>20.30 ± 0.02</td><td>19.25 ± 0.02</td><td>15.96 ± 0.05</td><td>15.75 ± 0.13</td></tr><tr><td>SDSS J023803.12+054526.1</td><td>J023803.161+054525.696</td><td>20.46 ± 0.02</td><td>18.91 ± 0.01</td><td>18.08 ± 0.01</td><td>15.30 ± 0.04</td><td>15.13 ± 0.09</td></tr><tr><td>2MASS J06453153-6646120</td><td>—</td><td>—</td><td>—</td><td>—</td><td>13.76 ± 0.02</td><td>13.31 ± 0.02</td></tr><tr><td>ULAS J075335.23+200622.4</td><td>J075335.224+200621.615</td><td>19.71 ± 0.01</td><td>18.23 ± 0.01</td><td>17.46 ± 0.01</td><td>14.60 ± 0.03</td><td>14.27 ± 0.05</td></tr><tr><td>ULAS J082206.61+044101.8</td><td>J082206.621+044100.926</td><td>20.33 ± 0.03</td><td>18.80 ± 0.03</td><td>17.98 ± 0.02</td><td>15.10 ± 0.03</td><td>14.93 ± 0.10</td></tr><tr><td>WISEA J101329.72-724619.2</td><td>—</td><td>—</td><td>—</td><td>—</td><td>14.59 ± 0.03</td><td>14.17 ± 0.04</td></tr><tr><td>2MASS J115820.877+0435014</td><td>J115821.180+043451.812</td><td>20.73 ± 0.03</td><td>18.65 ± 0.01</td><td>17.57 ± 0.01</td><td>13.70 ± 0.03</td><td>13.36 ± 0.03</td></tr><tr><td>ULAS J123142.99+015045.4</td><td>J123142.963+015045.403</td><td>21.32 ± 0.04</td><td>19.93 ± 0.02</td><td>19.19 ± 0.02</td><td>16.15 ± 0.06</td><td>15.81 ± 0.18</td></tr><tr><td>ULAS J124104.75-000531.4</td><td>J124104.737-000531.675</td><td>—</td><td>21.00 ± 0.06</td><td>20.16 ± 0.06</td><td>—</td><td>—</td></tr><tr><td>ULAS J124947.04+095019.8</td><td>J124946.950+095018.252</td><td>20.41 ± 0.01</td><td>19.01 ± 0.01</td><td>18.36 ± 0.02</td><td>15.96 ± 0.05</td><td>15.71 ± 0.13</td></tr><tr><td>ULAS J130710.22+151103.4</td><td>J130710.063+151102.511</td><td>—</td><td>21.27 ± 0.15</td><td>—</td><td>16.25 ± 0.07</td><td>15.99 ± 0.18</td></tr><tr><td>ULAS J133348.27+273505.5</td><td>J133348.294+273503.885</td><td>20.45 ± 0.02</td><td>18.94 ± 0.01</td><td>18.27 ± 0.02</td><td>—</td><td>—</td></tr><tr><td>ULAS J133836.96+022910.7</td><td>J133836.963-022911.152</td><td>—</td><td>20.55 ± 0.05</td><td>19.54 ± 0.05</td><td>15.78 ± 0.05</td><td>15.62 ± 0.12</td></tr><tr><td>ULAS J134206.86+053724.9</td><td>J134206.849+053724.420</td><td>21.55 ± 0.04</td><td>20.05 ± 0.04</td><td>19.23 ± 0.04</td><td>16.26 ± 0.06</td><td>16.27 ± 0.21</td></tr><tr><td>ULAS J134423.98+280603.8</td><td>J134423.954+280603.887</td><td>22.08 ± 0.09</td><td>20.36 ± 0.04</td><td>19.30 ± 0.05</td><td>15.61 ± 0.04</td><td>15.44 ± 0.10</td></tr><tr><td>ULAS J134749.79+333601.7</td><td>J134749.789+333601.640</td><td>19.85 ± 0.01</td><td>18.35 ± 0.01</td><td>17.53 ± 0.01</td><td>14.69 ± 0.03</td><td>14.33 ± 0.05</td></tr><tr><td>ULAS J134852.93+101611.8</td><td>J134852.814+101610.787</td><td>20.69 ± 0.02</td><td>19.22 ± 0.02</td><td>18.46 ± 0.03</td><td>15.47 ± 0.04</td><td>15.14 ± 0.08</td></tr><tr><td>ULAS J125226.62+092920.1</td><td>J125226.519+092920.150</td><td>20.74 ± 0.03</td><td>19.25 ± 0.02</td><td>18.47 ± 0.01</td><td>15.80 ± 0.05</td><td>15.52 ± 0.14</td></tr><tr><td>ULAS J135359.58+011856.7</td><td>J135359.573+011856.654</td><td>21.49 ± 0.04</td><td>19.94 ± 0.03</td><td>19.08 ± 0.03</td><td>15.87 ± 0.05</td><td>15.77 ± 0.13</td></tr><tr><td>WISEA J135501.90-825838.9</td><td>—</td><td>—</td><td>—</td><td>—</td><td>14.12 ± 0.03</td><td>13.55 ± 0.03</td></tr><tr><td>ULAS J141203.85+121609.9</td><td>J141203.779+121609.818</td><td>21.35 ± 0.03</td><td>19.44 ± 0.02</td><td>18.39 ± 0.03</td><td>14.77 ± 0.03</td><td>14.49 ± 0.05</td></tr><tr><td>SDSS J141624.12+134827.4</td><td>J141624.149+134827.768</td><td>18.35 ± 0.01</td><td>16.30 ± 0.01</td><td>15.21 ± 0.01</td><td>11.36 ± 0.02</td><td>11.03 ± 0.02</td></tr><tr><td>ULAS J141832.35+025323.0</td><td>J141832.360+025322.309</td><td>20.10 ± 0.01</td><td>18.56 ± 0.01</td><td>17.69 ± 0.01</td><td>14.86 ± 0.03</td><td>14.43 ± 0.05</td></tr><tr><td>ULAS J144151.55+043738.5</td><td>J144151.466+043735.872</td><td>21.71 ± 0.14</td><td>20.21 ± 0.03</td><td>19.29 ± 0.02</td><td>15.92 ± 0.05</td><td>15.54 ± 0.10</td></tr><tr><td>ULAS J151649.84+083607.1</td><td>J151649.778+083607.405</td><td>—</td><td>20.65 ± 0.05</td><td>19.57 ± 0.05</td><td>15.67 ± 0.05</td><td>15.43 ± 0.11</td></tr><tr><td>ULAS J154638.34-011213.0</td><td>J154638.320-011213.746</td><td>21.93 ± 0.19</td><td>20.37 ± 0.07</td><td>19.32 ± 0.05</td><td>—</td><td>—</td></tr><tr><td>2MASS J17561080+2815238</td><td>J175610.126+281517.780</td><td>19.03 ± 0.01</td><td>17.38 ± 0.01</td><td>16.41 ± 0.01</td><td>13.38 ± 0.02</td><td>13.07 ± 0.03</td></tr><tr><td>LSR J182611.3+301419.1</td><td>J182608.778+301410.509</td><td>15.48 ± 0.01</td><td>13.97 ± 0.01</td><td>13.22 ± 0.01</td><td>10.40 ± 0.02</td><td>10.04 ± 0.02</td></tr><tr><td>ULAS J223440.80+001002.6</td><td>J223440.785+001002.087</td><td>21.83 ± 0.06</td><td>20.51 ± 0.03</td><td>19.64 ± 0.08</td><td>16.47 ± 0.08</td><td>16.59 ± 0.33</td></tr><tr><td>ULAS J225902.14+115602.1</td><td>J225902.146+115601.795</td><td>21.16 ± 0.02</td><td>19.65 ± 0.03</td><td>18.76 ± 0.01</td><td>15.96 ± 0.06</td><td>15.81 ± 0.15</td></tr><tr><td>ULAS J230256.53+121310.2</td><td>J230256.527+121309.986</td><td>21.26 ± 0.04</td><td>19.84 ± 0.02</td><td>19.08 ± 0.04</td><td>—</td><td>—</td></tr><tr><td>ULAS J230443.30+093423.9</td><td>J230443.305+093423.871</td><td>21.52 ± 0.07</td><td>19.96 ± 0.10</td><td>19.02 ± 0.02</td><td>15.82 ± 0.05</td><td>15.39 ± 0.10</td></tr></tbody></table>
<figure><img src="image_18.png" /><figcaption>**Figure 15.** Polynomial fits to SDSS–PS1 and MKO–2MASS magnitude differences as a function of spectral type, for L0–T0 dwarfs (red crosses). The green shaded area indicates the rms of these fits. L and T dwarfs observed by SDSS and PS1 have *i* or *z* band photometric precision better than 0.12 mag. Those observed by 2MASS have *J*-, *H*-, and *K*- band photometric precision better than 0.16 mag.</figcaption></figure>

Gyr mass–*T*<sub>eff</sub> isochrones (fig. 5 of Burrows et al. 2001). We used the intersection points between the 10 Gyr mass–*T*<sub>eff</sub> isochrones with different metallicities to define the steady HBMM and identified the corresponding *T*<sub>eff</sub> of steady HBMM at 10 Gyr as a function of metallicity. We also found that the transition-zone occupied a very large area in the *T*<sub>eff</sub> versus [Fe/H] parameter space (fig 9 of Zhang et al. 2017b).

Fig. 20 shows solar metallicity isochrones of VLMS and substel-
lar objects from 0.01 to 10 Gyr. We also over plotted VLMS and BD
binaries with known dynamical masses in/near the BD transition-
zone. VLMS on the right-hand side of the transition-zone have
steady hydrogen fusion and could maintain their luminosity/*T*<sub>eff</sub>
throughout their very long lifetimes. Meanwhile, degenerate BDs
on the left-hand side of the transition-zone have no hydrogen fusion
to provide energy, thus cool and dim over time. In other words,
the luminosity/*T*<sub>eff</sub> distinction between degenerate BD and VLMS
populations is extended over time and the transition-zone is the
stretched area between them (also see fig 5 of Zhang et al. 2018).

Transitional BDs in the transition-zone are fully convective and
have sporadical unsteady hydrogen fusion in their cores which par-
tially contributes to their luminosity (Paper II & III). As the initial
thermal energy of a transitional BD is slowly dissipating over time,
the unsteady hydrogen fusion slowly becomes the dominate energy
sources to maintain its luminisoty. Since the core temperature of
a transitional BD declines slowly over time. The efficiency of the
fusion also declines very slowly over time. The fuison is very sen-
sitive to the mass of transitional BDs. Therefore, field transitional
BDs at a certain age could span in a wide T<sub>eff</sub> range within a narrow

mass range. Note such a T<sub>eff</sub> range is even wider at older age or
lower metallicity. The long-lasting unsteady hydrogen fusion in the
cores of transitional BDs slowed down their cooling speed and also
made them different from both VLMS and degenerate BDs.

Most of field degenerate BDs would have evolved into the mid-
to late-type T dwarf domain by their ages of ~0.5–5 Gyr. The
rest of them with relatively younger age or higher mass would
be crossing the L dwarf domain. This is the main reason for that L
dwarfs have a much lower number density than T dwarfs in the solar
neighbourhood (Burgasser 2007; Kirkpatrick et al. 2012). L dwarfs
are composed of VLMS, transitional BDs, and degenerate BDs
which are difficult to distinguish by observation without knowing
their mass or age. Field transitional BDs are mostly L dwarfs, but
could have spectral types of late-type M and early-type T, depending
on their mass and age. The BD transition-zone, which has not
drawn much attention in the past, further blurred the observational
stellar/substellar boundary in addition to the mass/age degeneracy.
This is why we have not reached an agreement on the observational
stellar/substellar boundary among field population.

4.5 The observational stellar/substellar boundary

Transitional BDs in the field are difficult to identify by observation
because of the mass/age degeneracy and relatively small luminos-
ity/temperature distinction (e.g. GD 165B; Kirkpatrick et al. 1999b).
However, transitional BDs in the halo are distributed in a very wide
range of luminosity/temperature after over ~10 Gyr of cooling to
reduce their initial thermal energy. Nine known halo transitional
BDs identified by their atmospheric parameters are summarized in
Paper III. They are located in a very broad halo BD transition-zone
in the T<sub>eff</sub> versus [Fe/H] parameter space.

To identify an empirical observational boundary between metal-
poor stars and transitional BDs, we compared known metal-poor
transitional BDs to least-massive stars in the halo, thick disc (in
fig. 9 of Paper II), and the solar neighbourhood (Dieterich et al.
2014) by their relative locations on colour-colour plots presented
in Section 4.3. Fig. 17 shows that metal-poor transitional BDs can
be well separated from stars by their *i* − *J* and *J* − *K* colours that
are sensitive to subtype (i.e. *T*<sub>eff</sub>) and subclass (i.e. [Fe/H]), respec-
tively. Therefore, we draw an empirical stellar-substellar boundary
in Fig. 17 that is shown as a black dashed broken line.

The stellar–substellar boundary was drawn, where we can sep-
arate known transitional BDs from least-massive stars as much as
possible by their *i* − *J* and *J* − *K* colours. Then we lead this boundary
to the right side of the two least-massive field stars suggested by Di-
eterich et al. (2014). Eight out of these nine known halo transitional
BDs in Paper III lie below/redward of this stellar–substellar bound-
ary. However, we note one apparent inconsistency, since UL0208
appears just above the boundary (blue open circle in the middle
of Fig. 17). We believe this is likely due to the uncertain *K*-band
photometry of this object (giving a larger *J* − *K* error), and further
note that this object lies below the boundary in Fig. 18 as a result of
its NIR *J* − *H* colour. It can also be seen that there are a significant
number of esdL stars clustered just above/blueward of the boundary,
and considerably fewer esdL subdwarfs below/redward of it. This
dearth of objects in the colour space on the substellar side of the
boundary represents a substellar subdwarf gap spanning mid L to
early T types. This region (the transition-zone) covers a wide T<sub>eff</sub>
range but a narrow mass range, which results in relatively fewer
objects compared to halo VLMS.

We also draw stellar/substellar boundaries on other colour–colour
plots composed of T<sub>eff</sub> and [Fe/H] sensitive colours (Figs 18 and 19).
<figure><img src="image_19.png" /><figcaption>**Figure 16.** Optical to infrared colours of L subdwarfs. Red hexagons, blue circles, and black diamonds represent sdL, esdL, and usdL subclasses. Third-order polynomial fits (black lines) and their rms (green shaded areas) for colour–spectral type relations of L dwarfs are also plotted for comparison.</figcaption></figure>
<table><caption>**Table 8.** Coefficients of third-order polynomial fits of colours as a function of spectral types (SpT) and second-order fits of absolute magnitudes for L0–L9 dwarfs in Figs 16 and 22. The fits are defined as $Colour/M = \sum_{i=0}^{3} c_i \times (SpT)^i$, $SpT = 10$ for L0 and $SpT = 19$ for L9. The rms of polynomial fits are listed in the last column. The *i* and *z* photometry is assumed to be on the SDSS (AB) system, and the *J*, *H*, and *K* photometry is assumed to be on the MKO (Vega) system.</caption><thead><tr><th>Colour/M</th><th>c<sub>0</sub></th><th>c<sub>1</sub></th><th>c<sub>2</sub></th><th>c<sub>3</sub></th><th>rms</th></tr></thead><tbody><tr><td>i - z</td><td>9.049</td><td>-1.544</td><td>1.025 × 10<sup>-1</sup></td><td>-1.979 × 10<sup>-3</sup></td><td>0.174</td></tr><tr><td>i - J</td><td>9.434</td><td>-1.227</td><td>9.063 × 10<sup>-2</sup></td><td>-1.888 × 10<sup>-3</sup></td><td>0.203</td></tr><tr><td>z - J</td><td>3.583</td><td>-0.408</td><td>4.134 × 10<sup>-2</sup></td><td>-1.175 × 10<sup>-3</sup></td><td>0.157</td></tr><tr><td>Y - J</td><td>-2.951</td><td>0.757</td><td>-4.367 × 10<sup>-2</sup></td><td>8.076 × 10<sup>-4</sup></td><td>0.106</td></tr><tr><td>J - H</td><td>2.564</td><td>-0.528</td><td>4.487 × 10<sup>-2</sup></td><td>-1.134 × 10<sup>-3</sup></td><td>0.093</td></tr><tr><td>J - K</td><td>4.205</td><td>-0.852</td><td>7.452 × 10<sup>-2</sup></td><td>-1.940 × 10<sup>-3</sup></td><td>0.146</td></tr><tr><td>Y - W1</td><td>4.176</td><td>-0.820</td><td>9.090 × 10<sup>-2</sup></td><td>-2.592 × 10<sup>-3</sup></td><td>0.285</td></tr><tr><td>J - W1</td><td>5.271</td><td>-1.169</td><td>9.641 × 10<sup>-2</sup></td><td>-2.388 × 10<sup>-3</sup></td><td>0.145</td></tr><tr><td>H - W1</td><td>7.892</td><td>-1.722</td><td>1.437 × 10<sup>-1</sup></td><td>-3.589 × 10<sup>-3</sup></td><td>0.213</td></tr><tr><td>K - W1</td><td>3.985</td><td>-0.915</td><td>7.167 × 10<sup>-2</sup></td><td>-1.703 × 10<sup>-3</sup></td><td>0.098</td></tr><tr><td>J - W2</td><td>11.065</td><td>-2.324</td><td>1.837 × 10<sup>-1</sup></td><td>-4.418 × 10<sup>-3</sup></td><td>0.258</td></tr><tr><td>W1 - W2</td><td>2.045</td><td>-0.371</td><td>2.416 × 10<sup>-2</sup></td><td>-4.647 × 10<sup>-4</sup></td><td>0.092</td></tr><tr><td>H - W2</td><td>8.883</td><td>-1.879</td><td>1.447 × 10<sup>-1</sup></td><td>-3.419 × 10<sup>-3</sup></td><td>0.192</td></tr><tr><td>W2 - W3</td><td>12.055</td><td>-2.088</td><td>1.314 × 10<sup>-1</sup></td><td>-2.658 × 10<sup>-3</sup></td><td>0.903</td></tr><tr><td>i<sub>P1</sub> - z<sub>P1</sub></td><td>5.857</td><td>-0.836</td><td>4.562 × 10<sup>-2</sup></td><td>-5.525 × 10<sup>-4</sup></td><td>0.202</td></tr><tr><td>i<sub>P1</sub> - y<sub>P1</sub></td><td>5.333</td><td>-0.560</td><td>2.874 × 10<sup>-2</sup></td><td>-2.007 × 10<sup>-4</sup></td><td>0.169</td></tr><tr><td>z<sub>P1</sub> - y<sub>P1</sub></td><td>1.654</td><td>-0.192</td><td>1.600 × 10<sup>-2</sup></td><td>-3.962 × 10<sup>-4</sup></td><td>0.070</td></tr><tr><td>i<sub>P1</sub> - J</td><td>10.030</td><td>-1.310</td><td>9.259 × 10<sup>-2</sup></td><td>-1.875 × 10<sup>-3</sup></td><td>0.392</td></tr><tr><td>z<sub>P1</sub> - J</td><td>1.456</td><td>0.065</td><td>1.232 × 10<sup>-2</sup></td><td>-5.777 × 10<sup>-4</sup></td><td>0.177</td></tr><tr><td>y<sub>P1</sub> - J</td><td>0.481</td><td>0.116</td><td>5.786 × 10<sup>-3</sup></td><td>-3.881 × 10<sup>-4</sup></td><td>0.127</td></tr><tr><td>i<sub>P1</sub> - H</td><td>9.785</td><td>-1.220</td><td>9.273 × 10<sup>-2</sup></td><td>-1.951 × 10<sup>-3</sup></td><td>0.412</td></tr><tr><td>y<sub>P1</sub> - H</td><td>1.724</td><td>-0.138</td><td>3.165 × 10<sup>-2</sup></td><td>-1.087 × 10<sup>-3</sup></td><td>0.189</td></tr><tr><td>z<sub>P1</sub> - H</td><td>2.847</td><td>-0.221</td><td>4.053 × 10<sup>-2</sup></td><td>-1.332 × 10<sup>-3</sup></td><td>0.228</td></tr><tr><td>y<sub>P1</sub> - K</td><td>2.444</td><td>-0.283</td><td>4.933 × 10<sup>-2</sup></td><td>-1.623 × 10<sup>-3</sup></td><td>0.259</td></tr><tr><td>M<sub>G</sub> (Gaia)</td><td>11.264</td><td>0.546</td><td>-3.690 × 10<sup>-3</sup></td><td>-</td><td>0.428</td></tr><tr><td>M<sub>BP</sub> (Gaia)</td><td>12.819</td><td>0.587</td><td>-6.725 × 10<sup>-3</sup></td><td>-</td><td>0.944</td></tr><tr><td>M<sub>BP</sub> (Gaia)</td><td>9.218</td><td>0.622</td><td>-7.554 × 10<sup>-3</sup></td><td>-</td><td>0.421</td></tr><tr><td>M<sub>i</sub> (PS1)</td><td>13.490</td><td>0.106</td><td>1.459 × 10<sup>-2</sup></td><td>-</td><td>0.400</td></tr><tr><td>M<sub>z</sub> (PS1)</td><td>9.250</td><td>0.602</td><td>-7.531 × 10<sup>-3</sup></td><td>-</td><td>0.369</td></tr><tr><td>M<sub>y</sub> (PS1)</td><td>8.581</td><td>0.566</td><td>-6.708 × 10<sup>-3</sup></td><td>-</td><td>0.361</td></tr></tbody></table>

The corresponding kink points on these boundaries in different plots have the same value for each colour. So that the boundaries in different colour–colour plots are consistent. Metal-poor transitional BDs are well separated from stars across the sdL, esdL, and usdL populations in Fig. 18 and these first nine plots of Fig. 19. The last three plots in Fig. 19 use the *K* – *W1* colour could only separate transitional BDs of esdL/usdL subclasses and invalid for sdL subclass.

Using the stellar/substellar boundaries from these plots in Figs 17–19, we estimate that about 22 of these 62 L subdwarfs in our photometric sample are likely metal-poor transitional BDs based on their spectral subtype/subclass and multiband colours. These substellar subdwarfs generally have spectral types of L3 or later except the usdL1.5 SD0104+15 that has extremely low metallicity corresponding to a higher HBMM.

Note that sdL3/esdL3/usdL1.5 are likely the earliest spectral types for metal-poor transitional BDs that are in the substellar subdwarf gap between VLMS and degenerate BDs. Halo degenerate BDs would have temperature below ~1000K and have cooled into the T- and Y-type region (Paper II; Burrows et al. 2001). Field degenerate BDs could have much earlier spectral types depending on their ages. Degenerate BDs do not have unsteady hydrogen fusion, but the most massive ones above ~0.055 $M_\odot$ (Fig. 20) would be able to fuse lithium. Therefore, ultra-cool dwarfs with lithium absorption lines in their spectra would fall in the degenerate BD domain (e.g. Teide 1; Rebolo et al. 1996). However, lithium absorption lines are not expected in the spectrum of halo

degenerate BDs because of their low temperature (section 3.5 of Paper III).

The least-massive field stars above the HBMM have spectral types of around L2.5 (2MASSI J0523382–140302; Dieterich et al. 2014) and L3 (2MASSI J1017075+130839B; Dupuy & Liu 2017). The latest spectral types of stars in our L subdwarf sample are sdL1/esdL1/usdL0. Note we do not have sdL2/esdL2/usdL1 subdwarfs in our sample; thus, we do not know if they could be stars, as they would be very close to the stellar boundary if they exist. The subtypes of L subdwarfs are assigned by comparing their optical spectra with those of L dwarf standards. L subdwarfs are hotter, more luminous, and more massive than L dwarfs with the same spectral type. Meanwhile, the HBMM is higher at lower metallicity (Papers II and III). As a consequence of these two facts, the latest spectral types of stars is slightly shifting to earlier subtypes across the dL, sdL, esdL, and usdL subclasses with decreasing metallicity.

## 5 GAIA OBSERVATIONS OF L SUBDWARFS

The ESA’s Gaia (Gaia Collaboration et al. 2016) astrometric survey have measured precise parallaxes and proper motions (Lindegren et al. 2018) for ~1.332 billion stars in its second data release (DR2; Gaia Collaboration et al. 2018). The Gaia survey has three optical pass-bands ($G$, $G_{BP}$, and $G_{RP}$) that are not very sensitive for ultra-cool dwarfs. Ultra-cool subdwarfs is slightly easier to be detected by Gaia than dwarfs with the same spectral types. As ultra-cool
<figure><img src="image_20.png" /><figcaption>**Figure 17.** The *i* − *J* versus *J* − *K* colours of L subdwarfs. Green open circles, red hexagons, blue circles, and black diamonds represent dL, sdL, esdL, and usdL populations, respectively. Error bars are smaller than the symbol size for some objects. Grey dots are 5000 point sources selected from a 10 deg² area of ULAS-SDSS-PS1-WISE sky with 14 < *J* < 16, which represents main-sequence stars (mostly with spectral types of FGK and M0–M4). The broken dashed line indicates an empirical stellar–substellar boundary (Paper II). The *i*-band magnitudes of UL2307 and VLMS (magenta filled five-pointed stars) in the local field (Dieterich et al. 2014) have been converted from *i*<sub>P1</sub> with equation (4). We converted 2MASS magnitudes to MKO magnitudes for L subdwarfs not observed in UKIDSS. The blue open circle on the right indicates 2M0616 with an estimated *i* − *J* colour. The blue open circle in the middle indicates UL0208 with uncertain *K*-band photometry. The BT-Dusty model grid with log g = 5.5 is plotted for comparison (with T<sub>eff</sub> and [Fe/H] indicated). M and L subtypes are marked at corresponding locations by their average colours. UL2332+12 (magenta open square) and SD1331 (magenta open pentagon) are joined by a magenta dashed line (see Section 3.4).</figcaption></figure>

<figure><img src="image_21.png" /><figcaption>**Figure 18.** The *z* − *J* versus *J* − *H* (left) and *J* − *W2* versus *J* − *K* (right) colour–colour plots for the L subdwarfs. Symbols are as described in Fig. 17, with additional magenta open squares to represent T dwarfs.</figcaption></figure>
<figure><img src="image_22.png" /><figcaption>Figure 19. Optical to infrared colour–colour plots for L subdwarfs. Symbols are as described in Fig. 17.</figcaption></figure>
<figure><img src="image_23.png" /><figcaption>**Figure 20** Solar metallicity isochrones for VLMS and BDs at 0.01, 0.05, 0.1, 0.5, 1, 5, and 10 Gyr. The black solid, red dashed and dotted lines are isochrones from Burrows et al. (1997), Baraffe et al. (2003), and Baraffe et al. (2015), respectively. A blue shaded band indicates the BD transition-zone at solar metallicity, which has a small shift between different model predications. M (black diamonds), L (green circles), and T (magenta squares) dwarfs in binary systems (joined with cyan lines) with dynamical mass measurements (Dupuy & Liu 2017; Lazorenko & Sahlmann 2018) are also over plotted. Blue open circles and yellow hexagons indicate objects with and without lithium detection in their spectra, respectively.</figcaption></figure>

subdwarfs have brighter magnitude in the optical (see fig. 3 in Zhang et al. 2013).

Twenty of 66 known L subdwarfs in Table 1 were observed by Gaia. Table 9 shows the Gaia astrometry and photometry of these 20 L subdwarfs. As our UKIDSS-SDSS survey is relatively deep, only two of our L subdwarfs are in Gaia DR2; UL0753+20 is from this work and ULAS J134749.79+333601.7 was presented in Paper I. The rest are previous known L subdwarfs.

Fig. 21 shows the HRD of 20 L subdwarfs in comparison to field stars and L dwarfs in Gaia DR2. L subdwarfs have slightly bluer $G - G_{RP}$ colour than L dwarfs; however, L subdwarfs do not appear as ‘sub’ dwarfs in the HRD of $G$-band absolute magnitude ($M_G$) versus $G - G_{RP}$. Because L subdwarfs not only have bluer $G - G_{RP}$ but also brighter $G$-band absolute magnitude than L dwarfs therefore are shifted to upper left towards to M dwarfs. Early- to mid-type L subdwarfs could be as bright as M dwarfs at $G$ band. L subdwarfs appear as ‘sub’ dwarfs in the HRD of $z_{P1} - y_{P1}$, $z_{P1} - J$, and $J - K$ colours at which they are much bluer than L dwarfs (see Fig. 16).

To better understand the properties of L subdwarfs and estimate the distance of L subdwarfs that are not in Gaia DR2, we studied the correlation between spectral type and absolute magnitude of L0–7 subdwarfs. Fig. 22 shows the first-order polynomial fits of relationships between their spectral types and Gaia, PS1, and MKO absolute magnitudes. Table 10 shows the coefficients of these fittings (see Table 8 for L0–9 dwarfs). The $G_{BP}$ band magnitudes of known L subdwarfs/dwarfs are likely close to Gaia’s detection limit; therefore, these $M_{GBP}$ in Fig. 22 (b) might not be reliable. In general, the $G$ to $J$ absolute magnitudes get brighter from dL to sdL and esdL/usdL subclasses. The esdL/usdL subclasses have similar $M_H$ to L dwarfs but fainter $M_K$. These are similar to what were shown in fig. 3 of Zhang et al. (2013). However, the $M_H$ of sdL subclass seems fainter at L0 but brighter at L7 than dL, esdL/usdL subclasses. The $M_K$ of sdL subclass is also fainter at L0 but brighter at L7 than esdL/usdL subclasses.

Fig. 23 shows the tangential velocities of 20 L subdwarfs observed by Gaia. These field stars in Fig. 23 have an median $V_{tan}$ of ~ 36 km s⁻¹. Five of these seven sdL subdwarfs have $V_{tan} > 75$ km s⁻¹. All the esdL/usdL subdwarfs have $V_{tan} > 120$ km s⁻¹, and most of them are between 200 and 400 km s⁻¹. The esdL and usdL subclasses generally have halo kinematics, which is consisted to the esdM/usdM subclasses (on the classification of Lépine et al. 2007). The sdL subclass mostly have thick disc kinematics. There are a few relatively more metal-poor ([Fe/H] ≈ −1) sdL subdwarfs have halo kinematics: UL0212+06, SD1333 (Paper I), and VVV J12564163–6202039 (Smith et al. 2018).

## 6 CONCLUSIONS AND FUTURE DIRECTIONS

We present the discovery of 27 L subdwarfs including 5 esdL and 22 sdL subdwarfs. These new objects were classified according to the L subdwarf classification scheme presented in Paper I. Six of these L subdwarfs have spectral types between L3 and L8 and are likely substellar objects, while the other 21 L0–L1 subdwarfs are likely VLMS. We measured their proper motions and estimated their spectroscopic distances. We also measured the RV of three that have X-shooter spectroscopy. Our SDSS-UKIDSS programme has confirmed/classified 35 L subdwarfs in total (amongst a full known population of 66), including 11 probable BDs and 24 VLMS.

We also interpret one of our candidates (UL2332+12) as a mildly metal-poor unresolved binary consisting of a blue ~L6p primary and a ~T4p secondary. UL2332+12 has a high probability of thick disc membership by its kinematics. Metal-poor BD binaries are rare, but their properties and binary fraction may be very useful for our understanding of substellar formation in the early Galaxy (e.g. Stamatellos & Whitworth 2009; Bate 2014).

We have assessed optical to mid-infrared colours of the L subdwarf population, using colour–spectral type and 2-colour diagrams, comparing with both L dwarfs and main sequence stars. We found that L subdwarfs of different metallicity subclasses can be well separated from L dwarfs and main sequence stars using a range of optical/infrared colours. Colour spaces have been identified in which preferential selection can be made of the full range of L subdwarf subclasses, as well as separating stellar and substellar subdwarfs. This analysis shows that the photometric systems employed by the PS1, VST, and VISTA surveys provide strong future potential for expansion of the known L subdwarf population out to its metallicity extremes, for which (based on the current sample) around a third will be substellar objects with the remaining two thirds VLMS.

We also note that the PS1 $i_{P1}$, $z_{P1}$, and $y_{P1}$ filters have similar wavelength coverage and transmission profiles to those planned for the next decade’s Large Synoptic Survey Telescope (LSST; LSST Science Collaboration et al. 2017) and the Chinese Space Station Optical Survey (CSS-OS). The LSST will observe the southern hemisphere in six pass-bands to single-visit depths of 23.4, 22.2, and 21.6 mag and co-added depths of 26.4, 25.2, and 24.4 mag in the *i*, *z*, and *y* bands, respectively, from 2022. A large area in the northern sky missed by the LSST will be covered by the CSS-OS. The CSS-OS will observe ~ 17 500deg² (|*b*| > 15 deg and |*Dec.*| > 20 deg) of the sky in seven bands to depths of 25.9, 25.2, and 24.4 mag in *i*, *z*, and *y*, respectively. The LSST could provide relative parallax distances for objects well brighter than its single-visit depth
<table><caption>**Table 9.** Gaia DR2 astrometry and photometry of 20 L subdwarfs.</caption><thead><tr><th>Name</th><th>Gaia DR2 ID</th><th>π (mas)</th><th>μ<sub>RA</sub> (mas yr<sup>-1</sup>)</th><th>μ<sub>Dec</sub> (mas yr<sup>-1</sup>)</th><th>G</th><th>G<sub>BP</sub></th><th>G<sub>RP</sub></th></tr></thead><tbody><tr><td>WISEA J001450.17-083823.4</td><td>2429285424478485504</td><td>20.14 ± 0.28</td><td>1477.22 ± 0.55</td><td>-257.96 ± 0.28</td><td>18.13</td><td>20.54</td><td>16.61</td></tr><tr><td>WISEA J011639.05-165420.5</td><td>2358397882610264960</td><td>16.45 ± 0.68</td><td>602.21 ± 1.34</td><td>67.57 ± 1.16</td><td>20.10</td><td>21.05</td><td>18.45</td></tr><tr><td>WISEA J020201.25-313645.2</td><td>50195829494784822400</td><td>15.14 ± 0.34</td><td>-155.09 ± 0.54</td><td>-1152.44 ± 0.36</td><td>19.18</td><td>21.36</td><td>17.59</td></tr><tr><td>WISEA J030601.66-033059.0</td><td>5183457632811832960</td><td>24.61 ± 0.30</td><td>331.48 ± 0.51</td><td>-1290.79 ± 0.52</td><td>18.57</td><td>20.59</td><td>16.96</td></tr><tr><td>ULAS J033351.10+001405.8</td><td>3264554308968041728</td><td>8.62 ± 0.67</td><td>772.09 ± 0.93</td><td>-65.70 ± 1.03</td><td>19.95</td><td>21.23</td><td>18.36</td></tr><tr><td>WISEA J043535.82+211508.9</td><td>144711230753602048</td><td>18.30 ± 0.55</td><td>871.71 ± 0.95</td><td>-964.06 ± 0.61</td><td>19.16</td><td>21.57</td><td>17.57</td></tr><tr><td>2MASS J05325346+8246465</td><td>558122277038055808</td><td>40.24 ± 0.64</td><td>2038.34 ± 1.09</td><td>-1663.73 ± 1.23</td><td>20.19</td><td>21.80</td><td>18.44</td></tr><tr><td>ULAS J075335.23+200622.4</td><td>673344017522527488</td><td>12.12 ± 0.84</td><td>-36.37 ± 1.45</td><td>-190.54 ± 0.97</td><td>20.28</td><td>21.22</td><td>18.59</td></tr><tr><td>SSSPM J10130734-1356204</td><td>3753081236588280192</td><td>17.87 ± 0.22</td><td>68.07 ± 0.34</td><td>-1029.79 ± 0.30</td><td>17.97</td><td>20.37</td><td>16.55</td></tr><tr><td>ULAS J124425.75+102439.3</td><td>3927262643141395584</td><td>9.49 ± 0.80</td><td>-441.87 ± 1.79</td><td>-545.11 ± 0.78</td><td>20.18</td><td>21.13</td><td>18.57</td></tr><tr><td>SDSS J125637.13-022452.4</td><td>3685444645661181696</td><td>12.55 ± 0.72</td><td>-511.25 ± 1.33</td><td>-300.64 ± 0.83</td><td>20.05</td><td>21.59</td><td>18.40</td></tr><tr><td>VVV J12564163-6202039</td><td>5863122429178232704</td><td>14.94 ± 1.68</td><td>-1129.81 ± 3.25</td><td>21.29 ± 2.56</td><td>20.75</td><td>-</td><td>-</td></tr><tr><td>ULAS J134749.79+333601.7</td><td>1458522725665649536</td><td>13.45 ± 0.90</td><td>62.43 ± 1.11</td><td>-13.61 ± 1.14</td><td>20.38</td><td>21.42</td><td>18.63</td></tr><tr><td>SDSS J141624.12+134827.4</td><td>1227133699053734528</td><td>107.56 ± 0.30</td><td>85.69 ± 0.69</td><td>129.07 ± 0.47</td><td>18.33</td><td>22.23</td><td>16.63</td></tr><tr><td>2MASS J16262034+3925190</td><td>1332410734823626240</td><td>32.49 ± 0.23</td><td>-1374.66 ± 0.35</td><td>237.36 ± 0.44</td><td>18.38</td><td>21.33</td><td>16.78</td></tr><tr><td>2MASS J16403197+1231068</td><td>4460894909281150208</td><td>10.15 ± 0.64</td><td>-222.25 ± 1.14</td><td>-134.90 ± 0.91</td><td>19.87</td><td>21.04</td><td>18.31</td></tr><tr><td>2MASS J17561080+2815238</td><td>4584405146372926720</td><td>28.94 ± 0.37</td><td>-613.74 ± 0.57</td><td>-413.40 ± 0.67</td><td>19.40</td><td>21.34</td><td>17.71</td></tr><tr><td>LSR J182611.3+301419.1</td><td>4588438567346043776</td><td>90.00 ± 0.11</td><td>-2290.54 ± 0.17</td><td>-683.13 ± 0.18</td><td>15.95</td><td>19.60</td><td>14.35</td></tr><tr><td>WISEA J204027.30+695924.1</td><td>2271357312343219456</td><td>30.30 ± 0.13</td><td>1558.38 ± 0.27</td><td>1697.36 ± 0.25</td><td>17.55</td><td>20.87</td><td>16.00</td></tr><tr><td>WISEA J213409.15+713236.1</td><td>2272533033868183168</td><td>9.21 ± 0.43</td><td>466.05 ± 0.68</td><td>466.05 ± 0.68</td><td>19.87</td><td>21.41</td><td>18.33</td></tr></tbody></table>

<figure><img src="image_24.png" /><figcaption>**Figure 21.** The Hertzsprung-Russell diagram for L subdwarfs in comparison to field objects. Symbols are as described in Fig. 17. Transitional BDs are indicated with magenta circles. Grey dots are objects selected from Gaia DR2, PS1, and LAS with distance <100pc, 180° < RA < 220° and 0° < Dec. < 20°. The two grey sequences are white dwarfs (left) and main-sequence stars (right). Some field stars are scattered mostly because they are too bright in the PS1 and UKIDSS fields.</figcaption></figure>

using the Gaia's reference frame. Meanwhile, the CSS-OS have much deeper single-visit depth than the LSST, therefore, is better in detecting higher proper motion cool objects beyond LSST's single-visit depth.

Furthermore, ESA's Euclid (Laureijs et al. 2011) space survey telescope is scheduled to launch in 2021, and aiming to observe half of the sky in four pass-bands to depths of 25 mag in VIS band and 24 mag in Y, J, and H bands. Its slit-less spectroscopy will observe the 0.92–1.85 µm wavelength to a depth of H ≈ 19.5 mag and could be used to identify T subdwarfs by their Y/J index (Burgasser, Burrows & Kirkpatrick 2006; Mace et al. 2013a). The

NASA's Wide-Field InfraRed Survey Telescope (WFIRST; Spergel et al. 2015) is planning to observe 2000 deg² of the sky in its high-latitude survey to depths of Y = 26.7, J = 26.9, H = 26.7, and F184 = 26.2 from mid-2020s.

These future optical and NIR sky surveys will provide great opportunity to study large numbers of extremely metal-poor L subdwarfs and halo degenerate BDs (esdT/Y and usdT/Y types) in the near future. In particular, the *z* to *H* bands (covered by these facilities) will probe very large volumes for L subdwarfs, with the most extreme examples remaining reasonably bright in these bands.
<figure><img src="image_25.png" /><figcaption>**Figure 22.** Correlations between spectral types and absolute magnitudes. The blue, red, and green lines indicate our polynomial fits for esdL/usdL, sdL, and dL subclasses. Their rms are indicated by shaded areas. The polynomial fits for field dwarfs in J, H, and K bands are from Dupuy & Liu (2012). Other symbols are as described in Fig. 17.</figcaption></figure>

The *K*-band flux of ultra-cool objects are the most sensitive wavelength to metallicity disparity. The *J* − *K* colour is very useful in the identification and characterization of ultra-cool subdwarfs, particularly for T subdwarfs that have similar *J*- and *H*-band spectral profile to T dwarfs but stronger suppressed flux in *K* band (e.g. Burningham et al. 2014). However, the *K*-band filter is not included in current survey strategies of both *Euclid* and *WFIRST*. The *Euclid* and *WFIRST* surveys will gain a lot more impacts on the science of ultra-cool objects, if it could include a *K*-band filter and extend the red cut-off wavelength of their NIR spectroscopy from ~1.9 to 2.2 µm. The *WFIRST* has a larger aperture size than the *Euclid*, and would have a better capability in the *K*-band detection of T subdwarfs, which become very faint.

We further discussed the BD transition-zone and properties of transitional BDs following G113 and G114. Degenerate BDs have an essentially different evolution from VLMS but their observational distinction are blurred by transitional BDs. Firstly, because the existence of the BD transition-zone was not widely realised.

Secondly, least-massive stars, transitional and degenerate BDs are mixed in the L dwarf domain. Although, in the L dwarf domain, young or massive degenerate BDs are crossing, field transitional BDs are making a long stay or slowly crossing, meanwhile older least-massive stars are more permanent (at early-type L). The BD transition-zone worth in-depth studies by modelling and observation, and that would help us to better understand observations of substellar populations. For example, there is a lack of objects at the L/T transition (Burgasser et al. 2007). This is firstly because the rapid evolution of BD atmospheres at ~1200 K stretched the spectral subtype sampling. Secondly, the L/T transition is at the bottom of the BD transition-zone and next to the abundant degenerate BDs that crossed the BD transition-zone.

Twenty of these 66 known L subdwarfs were observed by the Gaia with precise astrometry. L subdwarfs do appear as ‘sub’ dwarfs on the HRD with some specific colours (e.g. $z_{P1} - y_{P1}$, $z_{P1} - J$, and $J - K$). Their absolute magnitudes are brighter in optical to J band and fainter in K band than L dwarfs. The esdL and usdL subclasses generally have halo kinematics and the sdL subclass
<table><caption>**Table 10.** Coefficients of first-order polynomial fits of absolute magnitude as a function of spectral types (SpT) for L0–L7 subdwarfs in Fig. 22. The fits are defined as $M = c_0 + c_1 \times \text{SpT}$, SpT = 10 for L0, and SpT = 17 for L7.</caption><thead><tr><th><i>M</i><sub>abs</sub></th><th><i>c</i><sub>0</sub></th><th><i>c</i><sub>1</sub></th><th>rms</th><th>Subclass</th></tr></thead><tbody><tr><td><i>G</i> (Gaia)</td><td>12.441</td><td>0.3501</td><td>0.267</td><td>sdL</td></tr><tr><td><i>G</i> (Gaia)</td><td>10.191</td><td>0.4488</td><td>0.316</td><td>esdL/usedL</td></tr><tr><td><i>G</i><sub>BP</sub> (Gaia)</td><td>10.745</td><td>0.6878</td><td>0.884</td><td>sdL</td></tr><tr><td><i>G</i><sub>BP</sub> (Gaia)</td><td>12.339</td><td>0.4307</td><td>0.703</td><td>esdL/usedL</td></tr><tr><td><i>G</i><sub>RP</sub> (Gaia)</td><td>10.689</td><td>0.3620</td><td>0.223</td><td>sdL</td></tr><tr><td><i>G</i><sub>RP</sub> (Gaia)</td><td>8.627</td><td>0.4512</td><td>0.259</td><td>esdL/usedL</td></tr><tr><td><i>M</i><sub>i</sub> (PS1)</td><td>11.094</td><td>0.4394</td><td>0.274</td><td>sdL</td></tr><tr><td><i>M</i><sub>i</sub> (PS1)</td><td>8.385</td><td>0.5546</td><td>0.388</td><td>esdL/usedL</td></tr><tr><td><i>M</i><sub>z</sub> (PS1)</td><td>10.345</td><td>0.3627</td><td>0.241</td><td>sdL</td></tr><tr><td><i>M</i><sub>z</sub> (PS1)</td><td>8.569</td><td>0.4197</td><td>0.291</td><td>esdL/usedL</td></tr><tr><td><i>M</i><sub>y</sub> (PS1)</td><td>9.888</td><td>0.3248</td><td>0.184</td><td>sdL</td></tr><tr><td><i>M</i><sub>y</sub> (PS1)</td><td>8.890</td><td>0.3424</td><td>0.201</td><td>esdL/usedL</td></tr><tr><td><i>M</i><sub>J</sub> (MKO)</td><td>9.279</td><td>0.2242</td><td>0.200</td><td>sdL</td></tr><tr><td><i>M</i><sub>J</sub> (MKO)</td><td>7.883</td><td>0.3016</td><td>0.110</td><td>esdL/usedL</td></tr><tr><td><i>M</i><sub>H</sub> (MKO)</td><td>8.972</td><td>0.2157</td><td>0.167</td><td>sdL</td></tr><tr><td><i>M</i><sub>H</sub> (MKO)</td><td>7.339</td><td>0.3299</td><td>0.106</td><td>esdL/usedL</td></tr><tr><td><i>M</i><sub>K</sub> (MKO)</td><td>8.647</td><td>0.2144</td><td>0.187</td><td>sdL</td></tr><tr><td><i>M</i><sub>K</sub> (MKO)</td><td>6.760</td><td>0.3685</td><td>0.124</td><td>esdL/usedL</td></tr></tbody></table>

<figure><img src="image_26.png" /><figcaption>**Figure 23.** Tangential velocities of L subdwarfs (symbols are as described in Fig. 17) compared to field stars (grey dots as in Fig. 21). Note, error bars are smaller than the symbol size for some L subdwarfs. Transitional BDs are indicated with magenta circles. SDSS spectroscopic confirmed esdM and usedM subdwarfs in Gaia DR2 are indicated as green open circles and squares, respectively. These esdM and usedM have distances of <200pc and were shortened by 0.6 times for better comparison in the plot.</figcaption></figure>

mostly have thick disc kinematics, which is consisted to M subdwarfs. Five of these 20 L subdwarf in Gaia DR2 are transitional BDs. Degenerate BDs or T subdwarfs are too faint for the Gaia survey. However, we could get the precise astrometry of a T subdwarf if it has a bright companion observed by the Gaia. For example, WISE J200520.38+542433.9 (Mace et al. 2013b) is an sdT8 subdwarf companion (separated by 188.5 arcsec) to Wolf 1130 that has a Gaia DR2 distance of 16.5587 ± 0.0094 pc. Note that the uncertainty of the Gaia distance of Wolf 1130 is about two-thirds of the projected separation from its cool companion.

More wide binary systems contain both metal-poor BDs and FGKM stars are expected to be discovered in the near future (e.g.

Marocco et al. 2017). Such systems can be used as benchmarks to characterize metal-poor BDs and test new ultra-cool atmosphere models and substellar evolutionary models. As we could have precise measurements of distances (from Gaia) and abundances (from high resolution spectroscopy) of these bright FGKM primaries, which could be applied to their metal-poor BD companions.

## ACKNOWLEDGEMENTS

Based on observations made with the GTC, installed in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias, in the island of La Palma. Based on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere under ESO programmes 094.C-0202, 095.C-0878, 096.C-0130, 096.C-0974, and 0101.C-0626.

This work is based in part on data obtained as part of the UKIDSS. The UKIDSS project is defined in Lawrence et al. (2007). UKIDSS uses the UKIRT Wide Field Camera (WFCAM; Casali et al. 2007). The photometric system is described in Hewett et al. (2006), and the calibration is described in Hodgkin et al. (2009). The pipeline processing and science archive are described in Irwin et al. (2004) and Hambly et al. (2008).

Funding for the SDSS has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Aeronautics and Space Administration, the National Science Foundation, the U.S. Department of Energy, the Japanese Monbukagakusho, and the Max Planck Society. The SDSS website is http://www.sdss.org/. The SDSS is managed by the Astrophysical Research Consortium (ARC) for the Participating Institutions. The Participating Institutions are The University of Chicago, Fermilab, the Institute for Advanced Study, the Japan Participation Group, The Johns Hopkins University, Los Alamos National Laboratory, the Max-Planck-Institute for Astronomy (MPIA), the Max-Planck-Institute for Astrophysics (MPA), New Mexico State University, University of Pittsburgh, Princeton University, the United States Naval Observatory, and the University of Washington.

The Pan-STARRS1 Surveys (PS1) and the PS1 public science archive have been made possible through contributions by 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, the 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 Grant No. AST-1238877, the University of Maryland, Eotvos Lorand University (ELTE), the Los Alamos National Laboratory, and the Gordon and Betty Moore Foundation.

This publication makes use of data products from the 2MASS (Skrutskie et al. 2006) and the VLT Survey Telescope ATLAS survey (Shanks et al. 2015). This publication has made use of data from the VISTA Hemisphere Survey, ESO Programme, 179.A-2010 (PI: McMahon), and the VIKING survey from VISTA at the ESO Paranal Observatory, programme ID 179.A-2004. Data processing has been contributed by the VISTA Data Flow System at CASU, Cambridge and WFAU, Edinburgh. 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 Insti-
tute of Technology, funded by the National Aeronautics and Space
Administration.

This work presents results from the European Space Agency
(ESA) space mission Gaia. Gaia data are being processed by the
Gaia Data Processing and Analysis Consortium (DPAC). Funding
for the DPAC is provided by national institutions, in particular
the institutions participating in the Gaia MultiLateral Agreement
(MLA). The Gaia mission website is https://www.cosmos.esa.int/g
aia. The Gaia archive website is https://archives.esac.esa.int/gaia.

This research has benefitted from the M, L, T, and Y dwarf
compendium housed at http://DwarfArchives.org. This research has
benefited from the SpeX Prism Spectral Libraries, maintained by
Adam Burgasser at http://www.brownwarfs.org/spexprim. This
publication makes use of VOSA, developed under the Spanish
Virtual Observatory project supported from the Spanish MICINN
through grant AyA2008-02156.

ZHZ is supported by the PSL fellowship. DJP and HRAJ ac-
knowledge support from the UK's Science and Technology Fa-
cilities Council, grant numbers ST/M001008/1, ST/N001818/1,
and ST/J001333/1. MCGO acknowledges the financial support
of the Spanish Virtual Observatory project supported from the
Spanish ministry of Education and Competitiveness (MINECO)
through grant AYA2014-55216. NL is funded by the project num-
ber AYA2015-69350-C3-2-P from Spanish Ministry of Economy
and Competitiveness (MINECO). EM is funded by the MINECO
under grant AYA2015-69350-C3-1-P. DH is supported by Sonder-
forschungsbereich SFB 881 'The Milky Way System' (subproject
A4) of the German Research Foundation (DFG). FA received fund-
ing from the French 'Programme National de Physique Stellaire'
(PNPS) and the 'Programme National de Planétologie' of CNRS
(INSU). The computations of atmosphere models were performed
in part on the Milky Way supercomputer, which is funded by the
Deutsche Forschungsgemeinschaft (DFG) through the Collaborative
Research Centre (SFB 881) 'The Milky Way System' (sub-
project Z2) and hosted at the University of Heidelberg Computing
Centre, and at the Pôle Scientifique de Modélisation Numérique
(PSMN) at the École Normale Supérieure (ENS) in Lyon, and at the
Gesellschaft für Wissenschaftliche Datenverarbeitung Göttingen in
collaboration with the Institut für Astrophysik Göttingen. The au-
thors thank the referee, J. Davy Kirkpatrick, for the useful and
constructive comments.

REFERENCES

Allard F., Homeier D., Freytag B., 2014, ASInC, 11, 33

Bailer-Jones C. A. L., 2004, A&A, 419, 703

Baraffe I., Chabrier G., Barman T. S., Allard F., Hauschildt P. H., 2003, A&A, 402, 701

Baraffe I., Homeier D., Allard F., Chabrier G., 2015, A&A, 577, A42

Bardalez Gagliuffi D. C. et al., 2014, ApJ, 794, 143

Bate M. R., 2014, MNRAS, 442, 285

Bates D. R., 1952, MNRAS, 112, 40

Bessell M. S., 1991, AJ, 101, 662

Best W. M. J. et al., 2018, ApJS, 234, 37

Blake C. H., Charbonneau D., White R. J., 2010, ApJ, 723, 684

Bochanski J. J., West A. A., Hawley S. L., Covey K. R., 2007, AJ, 133, 531

Bohlin R. C., Colina L., Finley D. S., 1995, AJ, 110, 1316

Bowler B. P., Liu M. C., Dupuy T. J., 2010, ApJ, 710, 45

Burgasser A. J., 2004b, ApJ, 614, L73

Burgasser A. J., 2007, ApJ, 659, 655

Burgasser A. J., Burrows A., Kirkpatrick J. D., 2006, ApJ, 639, 1095

Burgasser A. J., Cruz K. L., Cushing M., Gelino C. R., Looper D. L., Faherty J. K., Kirkpatrick J. D., Reid I. N., 2010, ApJ, 710, 1142

Burgasser A. J., Cruz K. L., Kirkpatrick J. D., 2007, ApJ, 657, 494

Burgasser A. J., Kirkpatrick J. D., 2006, ApJ, 645, 1485

Burgasser A. J., Kirkpatrick J. D., Liebert J., Burrows A., 2003b, ApJ, 594, 510

Burgasser A. J., Liu M. C., Ireland M. J., Cruz K. L., Dupuy T. J., 2008a, ApJ, 681, 579

Burgasser A. J., Liu M. C., Ireland M. J., Cruz K. L., Dupuy T. J., 2008c, ApJ, 681, 579

Burgasser A. J., McElwain M. W., 2006, AJ, 131, 1007

Burgasser A. J., McElwain M. W., Kirkpatrick J. D., Cruz K. L., Tinney C. G., Reid I. N., 2004a, AJ, 127, 2856

Burgasser A. J., Vrba F. J., Lépine S., Munn J. A., Luginbuhl C. B., Henden A. A., Guetter H. H., Canzian B. C., 2008b, ApJ, 672, 1159

Burgasser A. J., Witte S., Helling C., Sanderson R. E., Bochanski J. J., Hauschildt P. H., 2009, ApJ, 697, 148

Burgasser A. J. et al., 2002, ApJ, 564, 421

Burgasser A. J. et al., 2003a, ApJ, 592, 1186

Burningham B., Smith L., Cardoso C. V., Lucas P. W., Burgasser A. J., Jones H. R. A., Smart R. L., 2014, MNRAS, 440, 359

Burningham B. et al., 2010, MNRAS, 404, 1952

Burrows A., et al., 1997, ApJ, 491, 856

Burrows A., Hubbard W. B., Lunine J. I., Liebert J., 2001, Rev. Mod. Phys., 73, 719

Casali M. et al., 2007, A&A, 467, 777

Cepa J. et al., 2000, SPIE, 4008, 623

Chambers K. C. et al., 2016, preprint (arXiv:1612.05560)

Chiu K., Fan X., Leggett S. K., Golimowski D. A., Zheng W., Geballe T. R., Schneider D. P., Brinkmann J., 2006, AJ, 131, 2722

Cruz K. L., Reid I. N., Liebert J., Kirkpatrick J. D., Lowrance P. J., 2003, AJ, 126, 2421

Cushing M. C., Looper D., Burgasser A. J., Kirkpatrick J. D., Faherty J., Cruz K. L., Sweet A., Sanderson R. E., 2009, ApJ, 696, 986

Cushing M. C. et al., 2011, ApJ, 743, 50 1171

Dhital S., West A. A., Stassun K. G., Bochanski J. J., Massey A. P., Bastien F. A., 2012, AJ, 143, 67

Dieterich S. B., Henry T. J., Jao W.-C., Winters J. G., Hosey A. D., Riedel A. R., Subasavage J. P., 2014, AJ, 147, 94

Dupuy T. J., Liu M. C., 2012, ApJS, 201, 19

Dupuy T. J., Liu M. C., 2017, ApJS, 231, 15

Faherty J. K. et al., 2016, ApJS, 225, 10

Fan X. et al., 2000, AJ, 119, 928

Filippazzo J. C., Rice E. L., Faherty J., Cruz K. L., Van Gordon M. M., Looper D. L., 2015, ApJ, 810, 158

Freudling W., Romaniello M., Bramich D. M., Ballester P., Forchi V., García-Dabló C. E., Moehler S., Neeser M. J., 2013, A&A, 559, A96

Gaia Collaboration et al., 2018, A&A, 616, A1

Gaia Collaboration et al., 2016, A&A, 595, A1

Gizis J. E., 1997, AJ, 113, 806

Gizis J. E., Harvin J., 2006, AJ, 132, 2372

Hambly N. C. et al., 2008, MNRAS, 384, 637

Hamuy M., Suntzeff N. B., Heathcote S. R., Walker A. R., Gigoux P., Phillips M. M., 1994, PASP, 106, 566

Hawley S. L. et al., 2002, AJ, 123, 3409

Hewett P. C., Warren S. J., Leggett S. K., Hodgkin S. T., 2006, MNRAS, 367, 454

Hodgkin S. T., Irwin M. J., Hewett P. C., Warren S. J., 2009, MNRAS, 394, 675

Irwin M. J., 2004, in Iye M., Moorwood A. F., eds, Proc. SPIE Conf. Ser. Vol. 5493, Optical and IR Telescope Instrumentation and Detectors. SPIE, Bellingham, p. 411

Jao W.-C., Henry T. J., Beaulieu T. D., Subasavage J. P., 2008, AJ, 136, 840

Jones H. R. A., Longmore A. J., Jameson R. F., Mountain C. M., 1994, MNRAS, 267, 413
Kendall T. R., Delfosse X., Martín E. L., Forveille T., 2004, A&A, 416, L17

Kirkpatrick J. D., 2005, ARA&A, 43, 195

Kirkpatrick J. D., Allard F., Bida T., Zuckerman B., Becklin E. E., Chabrier G., Baraffe I., 1999b, ApJ, 519, 834

Kirkpatrick J. D., Barman T. S., Burgasser A. J., McGovern M. R., McLean I. S., Tinney C. G., Lowrance P. J., 2006, ApJ, 639, 1120

Kirkpatrick J. D., Henry T. J., McCarthy D. W., Jr, 1991, ApJS, 77, 417

Kirkpatrick J. D. et al., 1999a, ApJ, 519, 802

Kirkpatrick J. D. et al., 2010, ApJS, 190, 100

Kirkpatrick J. D. et al., 2012, ApJ, 753, 156

Kirkpatrick J. D. et al., 2014, ApJ, 783, 122

Kirkpatrick J. D. et al., 2016, ApJS, 224, 36

Knapp G. R. et al., 2004, AJ, 127, 3553

Laureijs R. et al., 2011, preprint (arXiv:1110.3193)

Lawrence A. et al., 2007, MNRAS, 379, 1599

Lazorenko P. F., Sahlmann J., 2018, preprint (arXiv:1808.07835)

Lindegren L. et al., 2018, preprint (arXiv:1804.09366)

Lodieu N., Espinoza Contreras M., Zapatero Osorio M. R., Solano E., Aberasturi M., Martín E. L., 2012, A&A, 542, A105

Lodieu N., Espinoza Contreras M., Zapatero Osorio M. R., Solano E., Aberasturi M., Martín E. L., Rodrigo C., 2017, A&A, 598, A92

Lodieu N., Zapatero Osorio M. R., Martín E. L., Solano E., Aberasturi M., 2010, ApJ, 708, L107

LSST Science Collaboration et al., 2017, preprint (arXiv:1708.04058)

Luhman K. L., Sheppard S. S., 2014, ApJ, 787, 126

Lépine S., Rich R. M., Neill J. D., Caulet A., Shara M. M., 2002, ApJ, 581, L47

Lépine S., Rich R. M., Shara M. M., 2007, ApJ, 669, 1235

Mace G. N. et al., 2013a, ApJS, 205, 6

Mace G. N. et al., 2013b, ApJ, 777, 36

Manjavacas E. et al., 2016, MNRAS, 455, 1341

Marocco F. et al., 2015, MNRAS, 449, 3651

Marocco F. et al., 2017, MNRAS, 470, 4885

Martín E. L., Delfosse X., Basri G., Goldman B., Forveille T., Zapatero Osorio M. R., 1999, AJ, 118, 2466

Oke J. B., 1974, ApJS, 27, 21

Oke J. B., 1990, AJ, 99, 1621

Pinfield D. J. et al., 2012, MNRAS, 422, 1922

Rebolo R., Martin E. L., Basri G., Marcy G. W., Zapatero-Osorio M. R., 1996, ApJ, 469, L53

Reddy B. E., Lambert D. L., Allende Prieto C., 2006, MNRAS, 367, 1329

Reid I. N., Kirkpatrick J. D., Gizis J. E., Dahn C. C., Monet D. G., Williams R. J., Liebert J., Burgasser A. J., 2000, AJ, 119, 369

Reid I. N., Lewitus E., Burgasser A. J., Cruz K. L., 2006, ApJ, 639, 1114

Saumon D., Marley M. S., Abel M., Frommhold L., Freedman R. S., 2012, ApJ, 750, 74

Schmidt S. J., West A. A., Burgasser A. J., Bochanski J. J., Hawley S. L., 2010, AJ, 139, 1045

Schneider A. C., Greco J., Cushing M. C., Kirkpatrick J. D., Mainzer A., Gelino C. R., Fajardo-Acosta S. B., Bauer J., 2016, ApJ, 817, 112

Scholz R.-D., Lehmann I., Matute I., Zinnecker H., 2004b, A&A, 425, 519

Scholz R.-D., Lodieu N., McCaughrean M. J., 2004a, A&A, 428, L25

Shanks T. et al., 2015, MNRAS, 451, 4238

Sivarani T., Lépine S., Kembhavi A. K., Gupchup J., 2009, ApJ, 694, L140

Skrutskie M. F. et al., 2006, AJ, 131, 1163

Skrzypek N., Warren S. J., Faherty J. K., 2016, A&A, 589, A49

Smith L., Lucas P. W., Burningham B., Jones H. R. A., Smart R. L., Andrei A. H., Catalán S., Pinfield D. J., 2014, MNRAS, 437, 3603

Smith L. C. et al., 2018, MNRAS, 474, 1826

Spergel D. et al., 2015, preprint (arXiv:1503.03757)

Stamatellos D., Whitworth A. P., 2009, MNRAS, 392, 413

Tsuji T., Ohnaka K., Aoki W., Nakajima T., 1996, A&A, 308, L29

Vernet J. et al., 2011, A&A, 536, A105

Wright E. L. et al., 2010, AJ, 140, 1868

York D. G. et al., 2000, AJ, 120, 1579

Zapatero Osorio M. R., Martín E. L., Béjar V. J. S., Bouy H., Deshpande R., Wainscoat R. J., 2007, ApJ, 666, 1205

Zhang Z. H., Homeier D., Pinfield D. J., Lodieu N., Jones H. R. A., Allard F., Pavlenko Y. V., 2017b, MNRAS, 468, 261 (Paper II)

Zhang Z. H. et al., 2009, A&A, 497, 619

Zhang Z. H. et al., 2010, MNRAS, 404, 1817

Zhang Z. H. et al., 2013, MNRAS, 434, 1005

Zhang Z. H. et al., 2017a, MNRAS, 464, 3040 (Paper I)

Zhang Z. H. et al., 2018, MNRAS, 479, 1383 (Paper III)

## SUPPORTING INFORMATION

Supplementary data are available at MNRAS online.

Please note: Oxford University Press is not responsible for the content or functionality of any supporting materials supplied by the authors. Any queries (other than missing material) should be directed to the corresponding author for the article.

This paper has been typeset from a *TEX/LATEX* file prepared by the author.