High-resolution spectroscopic atlas of M subdwarfs. Effective temperature and metallicity 111Based on observations made with the ESO Very Large Telescope at the
Paranal Observatory under programme 087.D-0586A and 087.D-0586B 222Tables 1 and 2 along with reduced spectra used for the model comparison are available
at the CDS via anonymous ftp to
cdsarc.u-strasbg.fr (22.214.171.124) or via
Key Words.:Stars: low-mass – subdwarfs – atmospheres
Context:M subdwarfs are metal-poor and cool stars. They are important probes of the old galactic populations. However, they remain elusive because of their low luminosity. Observational and modelling efforts are required to fully understand their physics and to investigate the effects of metallicity in their cool atmospheres.
Aims:We performed a detailed study of a sample of subdwarfs to determine their stellar parameters and constrain the state-of-the art atmospheric models.
Methods:We present UVES/VLT high-resolution spectra of three late-K subdwarfs and 18 M subdwarfs. Our atlas covers the optical region from 6400 to the near-infrared at 8900 . We show spectral details of cool atmospheres at very high-resolution (R 40 000) and compare them with synthetic spectra computed from the recent BT-Settl atmosphere models.
Results:Our comparison shows that molecular features (TiO, VO, CaH), and atomic features (Fe I, Ti I, Na I, K I) are well-fitted by current models. We produce a relation of effective temperature versus spectral type over the entire subdwarf spectral sequence. The high-resolution of our spectra, unable us to perform a detailed comparison of the line profiles of individual elements such as Fe I, Ca II, and Ti I, and we determined accurate metallicities of these stars. These determinations in turn enable us to calibrate the relation between metallicity and molecular-band strength indices from low-resolution spectra.
Conclusions:The new generation of models is able to reproduce various spectral features of M subdwarfs. These high-resolution spectra allowed us to separate the atmospheric parameters (effective temperature, gravity, metallicity), which is not possible when using low-resolution spectroscopy or photometry.
M subdwarfs are metal-poor, low-luminosity dwarfs. Very few are known, and as a result the bottom-end of the main sequence is very poorly constrained for metal-poor stars. The search for M subdwarfs is hampered not only by the fact that metal-poor stars are rare and intrinsically faint, but also because late-type M subdwarfs do not show exceptionally red colours, as ultracool-M dwarfs and brown dwarfs do (Lépine et al., 2003a). Subdwarfs are typically very old (10 Gyr or older) and belong to the old Galactic populations: old disc, thick disc and spheroid, as shown by their spectroscopic features, kinematic properties and ages (Digby et al., 2003; Lépine et al., 2003a; Burgasser et al., 2003). Because the low-mass subdwarfs, with their extremely long nuclear-burning lifetimes, were presumably formed early in the Galaxy’s history, they are important tracers of the Galactic structure and chemical enrichment history. In addition, detailed studies of their complex spectral energy distributions give new insights on the role of metallicity in the opacity structure, chemistry, and evolution of cool atmospheres, and on fundamental questions on spectral classification and temperature-vs-luminosity scales. Gizis (1997) proposed a first classification of M subdwarfs (sdM) and extreme subdwarfs (esdM) based on TiO and CaH band strengths in low-resolution optical spectra. Lépine et al. (2007) have revised the adopted classification, and proposed a new classification for those most metal-poor, the ultra-subdwarfs (usdM). Jao et al. (2008) compared model grids with the optical spectra to characterize the M subdwarfs by three parameters: temperature, gravity, and metallicity, and thus gave an alternative classification scheme of subdwarfs.
Spectroscopic studies of sdM at high-resolution have proven to be a difficult task. In the low-temperature regime occupied by these stars, the optical spectrum is covered by a forest of molecular lines, hiding or blending most of the atomic lines used in spectral analysis. However, over the past decade, stellar atmosphere models of very low mass stars have made much progress by exploring metallicity effects (Allard et al., 2012, 2013). One of the most important recent improvement is the revision of the solar abundances (Asplund et al., 2009; Caffau et al., 2011).
Rapid progress in the investigation of cool atmospheres is expected thanks to the advent of eight-meter-class telescopes that allow high-resolution spectroscopy of these faint targets. In high-resolution spectra, access to weak lines allows us to determine the metallicity separately from the other main parameters gravity and effective temperature. The pressure changes equally affect all atmospheric parameters and hence the various absorption bands. The determination of gravity from the pressure-broadened wings is expected to be much more accurate than comparing colour ratios from photometry or low-resolution spectra. Thus it is necessary to achieve a very good fit in all important absorbers to determine atmospheric properties, because the chemical complexity of these atmospheres reacts sensitively to the main opacities. Descriptions of these stars therefore need validations by comparison with high-resolution spectroscopic observations.
Measuring metallicities for M dwarfs is also challenging. Over the entire M dwarf sequence, the effective temperature ranges from about 4000 K to 2400 K (see e.g. Rajpurohit et al., 2013). With decreasing temperature, the spectra show increasingly abundant diatomic and triatomic molecules. In particular, the TiO and HO bands have complex and extensive absorption structures, creating a pseudo-continuum that let pass only the strongest, often resonance, atomic lines. Woolf & Wallerstein (2005) and Woolf et al. (2009) obtained metallicities from high-resolution spectra by measuring equivalent widths of atomic lines in regions less dominated by molecular bands. Metallicities have also been obtained in binaries containing an M-type and a solar-type star from the much better understood spectra of the latter (Bonfils et al., 2005; Bean et al., 2006a, b). These results, together with the abundances obtained from high-resolution spectra, are combined to calibrate metallicities using photometry or molecular indices (Bonfils et al., 2005; Woolf et al., 2009; Casagrande et al., 2008).
Metal-poor stars are rare in the solar neighbourhood, and the current sample of local subdwarfs is very limited. High-resolution spectra of M dwarfs were shown by Tinney & Reid (1998) on the full optical range, but such observations are not available for the M subdwarfs on the whole temperature sequence. In this paper we present the first high-resolution optical atlas of stars that covers the whole sdM, esdM, and usdM sequence. It consists of 2 sdK, 11 sdM, 1 esdK, 5esdM, and 2 usdM observed with UVES at VLT. Using the most recent PHOENIX BT-Settl stellar atmosphere models, we performed a detailed comparison with our observed spectra. In this study we compare the models with the high-resolution spectra of sdM at subsolar metallicities and assign effective temperatures. We derive metallicities based on the best fit of synthetic spectra to the observed spectra and compare in detail the line profile of Fe I, Ca II, Ti I, and Na.
High-resolution spectroscopic observations are described in Section 2 and the model atmospheres used for comparison are presented in Section 3. The comparison between observed and synthetic spectra and the derived stellar parameters are given in Section 4. In Section 5, we present effective temperature versus colours and spectral-type relations, and the metallicity calibrations. The conclusion is given in Section 6.
2 High-resolution spectral atlas of M subdwarfs
High-resolution spectroscopy is a very important tool for understanding the physics of stars. Identifying of atomic and molecular absorption features in the spectra of M subdwarfs is important they can be used to constrain the main stellar parameters , , [Fe/H]. The determination of gravity from its effect on the pressure-dependent wings of the saturated atomic absorption lines (predominantly of the alkali elements) is expected to be much more accurate than comparing colour ratios from low-resolution spectra (Reiners et al., 2007).
We present a high-resolution spectral atlas of 21 very low mass objects that we selected to cover the entire M subdwarf spectral range, including 6 extreme and 2 ultra subdwarfs. The name, spectral type, and near-infrared photometry of these objects are given in Table 1. The photometry was taken from the Two Micron All Sky Survey (2MASS, Skrutskie et al., 2006). The spectral types and spectral indices were found in the literature.
|LHS 72||23 43 13.5||24 09 50||sdK4||9.61||9.04||8.82||–||–||–||a,–|
|LHS 73||23 43 16.6||24 11 14||sdK7||10.11||9.59||9.37||0.960||0.847||0.901||b,b|
|G 18-37||22 14 55.5||+05 42 37||esdK7||12.72||12.20||12.02||–||–||–||c,–|
|APMPM J2126-4454||21 26 23.9||44 53 34||sdM0||12.65||12.11||11.92||–||–||–||d,–|
|LHS 300||11 11 13.8||41 05 33||sdM0||10.48||10.01||9.80||1.031||0.886||0.943||e,e|
|LHS 401||15 39 39.1||55 09 10||sdM0.5||10.15||9.60||9.41||0.997||0.911||0.958||e,e|
|LHS 158||02 42 02.8||44 30 59||sdM1||10.43||9.94||9.73||0.732||0.639||0.829||e,e|
|LHS 320||12 02 33.7||+08 25 51||sdM2||10.74||10.18||9.99||0.620||0.512||0.751||f,f|
|LHS 406||15 43 18.4||20 15 31||sdM2||9.78||9.23||9.02||0.686||0.576||0.789||e,e|
|LHS 161||02 52 45.6||+01 55 50||esdM2||11.71||11.20||11.00||0.889||0.689||0.817||f,e|
|LP 771-87||03 07 34.0||17 36 38||usdM2||15.60||15.29||14.90||–||–||–||g,–|
|LHS 541||23 17 06.0||13 50 53||sdM3||13.03||12.56||12.41||0.730||0.467||0.664||h,h|
|LHS 272||09 43 46.3||17 47 07||sdM3||9.62||9.12||8.87||0.678||0.441||0.669||f,f|
|LP 707-15||01 09 54.1||10 12 13||esdM3||12.94||12.40||12.16||–||–||–||c,–|
|LSR 1755+1648||17 55 32.8||+16 48 59||sdM3.5||11.35||10.89||10.63||0.486||0.419||0.664||i,i|
|LHS 375||14 31 38.3||25 25 33||esdM4||12.15||11.67||11.51||0.829||0.372||0.547||f,f|
|LHS 1032||00 11 00.8||+04 20 25||usdM4.5||14.34||13.81||13.76||0.944||0.374||0.514||j,j|
|SSSPM J0500-5406||05 00 15.8||54 06 27||esdM6.5||14.44||14.12||13.97||0.755||0.220||0.331||h,h|
|LHS 377||14 39 00.3||+18 39 39||sdM7||13.19||12.73||12.48||0.232||0.205||0.396||f,f|
|APMPM J0559-2903||05 58 58.9||29 03 27||esdM7||14.89||14.45||14.46||0.600||0.210||0.320||l,l|
|SSSPM J1013-1356||10 13 07.3||13 56 20||sdM9.5||14.62||14.38||14.40||0.208||0.116||0.200||m,m|
(a) Rodgers & Eggen (1974); Bidelman (1985) – (b) Reylé et al. (2006) – (c) S. Lépine, private communication – (d) Scholz et al. (2002) – (e) Jao et al. (2008) – (f) Gizis (1997) – (g) Kirkpatrick et al. (2010) – (h) Dawson & De Robertis (2000) – (i) Lépine et al. (2003b) – (j) Lépine et al. (2007) – (k) Burgasser & Kirkpatrick (2006) – (l) Schweitzer et al. (1999) – (m) Scholz et al. (2004)
2.1 Observation and data reduction
The observations were carried out in visitor mode during April and September 2011 with the optical spectrometer UVES (Dekker et al., 2000) on the Very Large Telescope (VLT) at the European Southern Observatory (ESO) in Paranal, Chile. UVES was operated in dichroic mode using the red arm with non-standard setting centred at 830 nm. This setting covers the wavelength range 6400 -9000 , which contains various atomic lines such as Fe I, Ti I, K II, Na I and Ca II and is very useful for the spectral synthesis analysis. Between the two CCDs of the red arm, the spectra have a gap from 8200 to 8370 . The Na doublet at 8190 lies just blueward of the gap. The spectra were taken with a slit width of 1.0, yielding a nominal resolving power of . The signal-to-noise ratio varies over the wavelength region according to the object’s spectral energy distribution and detector efficiency. It reaches 30 to more than 100, depending on the magnitude of the targets in most of the spectral region so that the dense molecular and atomic absorption features are clearly discernible from noise. Data were reduced using the software called Reflex for UVES data, which runs standard ESO pipelines modules.
The spectra are shown in Fig. 1 for a representative sample along the sdM sequence, and in Fig. 2 for the esdM and usdM targets. The spectra were not corrected for terrestrial absorption by O and HO. The O features are quite regular and straightforward to distinguish from features in the stars themselves. This is not true for the HO absorption, which is complex and irregular (Tinney & Reid, 1998). We therefore show in the upper panels the spectrum of the reference star EG 21, where telluric absorptions appear. The main molecular and atomic features expected in the M subdwarfs are labelled in Fig.1 and Fig. 2.
2.2 Molecular features
The optical spectra are dominated by molecular absorption bands from metal oxide species such as titanium oxide (TiO), vanadium oxide (VO), and hydrides such as CaH and HO. They are the most significant opacity sources. However, due to the low metallicity of subdwarfs, they are TiO depleted. The primary effects are the strengthening of hydride bands and collision-induced absorption (CIA) by H, and the broadening of atomic lines (Allard et al., 1997). Unlike TiO, which produces distinctive band heads degraded on the red, VO produces more diffuse absorption. CaH hydride bands are significant opacity sources, but decrease in relative strength and become saturated with decreasing temperature.
2.3 Atomic lines
All the observed M subdwarfs show strong alkali lines in the observed wavelength range. They are massively pressure broadened, as expected from their high-gravity surface. Atomic features such as Ca II, K I, Rb I, Na I, Ti I, and Mg I are visible throughout the sequence, and their lines are prominent in almost all of the spectra. However, in regions where strong atmospheric absorption is present, it is difficult to measure the intensities of these lines.
The Na I lines at 8183 and 8194 are clearly visible in all the observed spectra and become broadened from hotter to cooler M subdwarfs. In our setting they appear just at the red end of the lower chip of the blue arm. The K I lines are very narrow for early-type subdwarfs but become very wide and smooth in late-type subdwarfs. The K I resonance lines at 7665 and 7698 govern the spectral shape of cool subdwarf spectra. The equivalent width of these K I lines are of several hundred . The ionized Ca II triplet lines at 8498 , 8542 , and 8662 are very strong in all the observed spectra. Their detailed study by Mallik (1997) shows that their strengths depend on stellar parameters like luminosity, temperature, and metallicity. They are ideal candidates to study their sensitivity to various stellar parameters in cool stars. Although the lower levels of the Ca II triplet lines are populated radiatively and are not collision controlled, they have been identified as very good luminosity probes by Mallik (1997). They are relatively free from blends and are little contaminated by telluric lines. In contrast, the Na I lines at 8183 and 8195 , which have also been used as luminosity probes in cool stars, have several atmospheric absorption lines in their vicinity (Alloin & Bica, 1989; Zhou, 1991).
3 Model atmospheres
We used the most recent BT-Settl models that were partially published in a review by Allard et al. (2012). These atmosphere models are computed with the PHOENIX multi-purpose atmosphere code version 15.5 (Hauschildt et al., 1997; Allard et al., 2001), solving the radiative transfer in 1D spherical symmetry, with the classical assumption of: hydrostatic equilibrium, convection using the mixing-length theory, chemical equilibrium, and a sampling treatment of the opacities. The models use a mixing length as derived by the radiation hydrodynamic simulations of Ludwig et al. (2002, 2006), and Freytag et al. (2012).
Compared with previous models by Allard et al. (2001), the current version of the BT-Settl model atmosphere uses the BT2 water-vapor line list computed by Barber et al. (2006), TiO, VO, CaH line lists by Plez (1998), MgH by Weck et al. (2003) and Skory et al. (2003), FeH and CrH by Dulick et al. (2003) and Chowdhury et al. (2006), CO by Tashkun et al. (2004), H CIA by Borysow et al. (2001); Abel et al. (2011), and CO by Goorvitch & Chackerian (1994a, b), to mention the most important line lists. We also included the H-He CIA from Borysow & Frommhold (1989) and Borysow et al. (1997), H-H from Gustafsson et al. (2003), and H-He (at least in the most recent models) from Gustafsson & Frommhold (2001). Detailed profiles for the alkali lines were also used (Allard et al., 2007). The reference solar elemental abundances used in this version of the BT-Settl models those measured by Caffau et al. (2011).
In general, the Unsold (1968) approximation is used, with the general exception of Na I, Si I, Ca I, and Fe I, where we instead used the respective correction factors found by Gustafsson et al. (2008) for the atomic damping constants with a correction factor to the widths of 2.5 for the non-hydrogenic atoms (Valenti & Piskunov, 1996). More accurate broadening data for neutral hydrogen collisions by Barklem et al. (2000) were included for several important atomic transitions such as the alkali, Ca I, and Ca II resonance lines. For molecular lines, we adopted average values (e. g. for water vapor lines) from the HITRAN database (Rothman et al., 2009), which are scaled to the local gas pressure and temperature
with a single temperature exponent of 0.5, to be compared with values ranging mainly from 0.3 to 0.6 for water transitions studied by Gamache et al. (1996). The HITRAN database gives widths for broadening in air, but Bailey & Kedziora-Chudczer (2012) found that these agree in general within 10 – 20% with those for broadening by a solar-composition hydrogen-helium mixture.
The BT-Settl grid extends from = 300 to 7000 K in steps of 100 K, = 2.5 to 5.5 in steps of 0.5, and [M/H]= -2.5 to 0.0 in steps of 0.5, which accounts for alpha-element enrichment. The alpha enhancement was taken as [alpha/Fe] = -0.4[Fe/H] for 0 [Fe/H] -1 and 0.4 for metallicities -1.0, which is also consistent with the choice of Gustafsson et al. (2008). These different prescriptions for alpha enhancement are rough estimates for the old thin disc and thick disc, respectively (see e.g. Neves et al., 2009, for a discussion on the abundance trends relative to Fe).
We linearly interpolated the grid at every 0.1 dex in and metallicity. For more details of the BT-Sett model atmosphere see Allard et al. (2012, 2013) andRajpurohit et al. (2012). The synthetic colours and spectra were distributed with a resolution of around via the PHOENIX web simulator444http://phoenix.ens-lyon.fr/simulator.
Fig. 3 shows BT-Settl synthetic spectra varying and [M/H]. Oxide bands that dominate in M dwarf spectra are weaker in the subdwarfs where the hydride bands dominate (Jao et al., 2008). They have complex and extensive band structures that leave no window for the true continuum and create a pseudo-continuum that only lets pass the strongest, often resonance atomic lines (Allard, 1990; Allard & Hauschildt, 1995). However, because of the lower metallicity of subdwarfs, the TiO bands are weaker, and the pseudo-continuum is brighter. This increases the contrast to the other opacities such as hydride bands and atomic lines, which feel the higher pressures of the deeper layers where they emerge from. We therefore see these molecular bands with more detaile and under more extreme gas pressure conditions than for M dwarfs.
Although the models were computed for different [M/H] values, in the following we use [Fe/H] everywhere. Indeed, we use only the iron lines for the metallicity determinations so that in fact we determine a [Fe/H] value.
4 Comparison with model atmospheres
We perform a comparison between observed and synthetic spectra computed from the BT Settl model to derive the physical parameters of our sample. Furthermore, the comparison with observed spectra is very crucial to reveal the inaccuracy or incompleteness of the opacities used in the model. Figures 4 and 5 show the comparison of the best-fit model for an sdM1 and a usdM4.5 star.
4.1 Molecular bands
The spectra agree well and reproduce the specific strengths of the TiO band heads at 6600 , 6700 , 7050 7680, and 8859 , and of the VO bands at 7000 , 7430 , and 7852 . The excellent match between models and observations over entire subdwarf sequence shows that the high-frequency pattern visible at this spectral resolution is the structure of the absorption band and not noise.
4.2 Atomic lines
The models predict the shape of the Na I doublet at 8194 and 8183 rather well but its strength is well fitted only poorly. In the sdM3.0 and later and in esdM, the observed lines are broader and shallower than those predicted by the models.
The qualitative behaviour of the K I doublet at 7665 and 7698 is well reproduced by the models, especially the strong pressure-broadening wings in the early sdM and esdM. In the sdM0, the cores of the observed K I lines are still visible as relatively narrow absorption minima embedded in wings extending a few tens to one hundred . This broader absorption component becomes saturated in sdM7 spectrum. The models also show this effect but do not reproduce it perfectly. This may be an indication (i) that the models do not yet produce correct densities of the neutral alkali metals in the uppermost part of the atmosphere, since the central parts of the lines form in the highest layers of the atmosphere, especially in the early sdM and esdM where alkali metals are not depleted strongly (Johnas et al., 2007), or (ii) that the depth of atomic lines can only be reproduced when accounting for a magnetic field and Zeeman broadening (Deen, 2013).
The models reproduce the strength and wings of Ca II triplet lines at 8498 , 8542 , and 8662 very well. We also have a good fit to the Ti I lines. They are located between 8400 to 9700 and belong to low-energy transitions that appear to be visible in such cool atmospheres (see also Reiners & Basri, 2006).
Absorption lines of Rb appear in spectra later than sdM7.0. They become stronger and wider towards lower temperature. Given the fact that they are embedded in pseudo-continua that are not always a good fit to the data, the behaviour of the two Rb lines at 7800 and 7948 is very well reproduced by the models.
4.3 Stellar parameter determination
The analysis using synthetic spectra requires the specification of several input parameters: effective temperature, surface gravity, and the overall metallicity with respect to the Sun. We first convolved the synthetic spectrum with a Gaussian kernel at the observed resolution and then interpolated the result with the observation. We performed a fitting between the grid of synthetic spectra and observed spectra in the wavelength range 6400 to 8900 , as shown by Önehag et al. (2012). We excluded the spectral region between 6860 to 6960 , 7550 to 7650 , and 8200 to 8430 because of atmospheric absorption. We let all the stellar parameters (, [Fe/H], ) vary.
The minimum value gives the best-fit parameters and their error bars were derived by allowing a 5 % deviation of the value from its minimum value. The acceptable parameters were finally inspected by comparing them with the observed spectra.
The TiO system at 6600 , 6700 , and 7100 is highly sensitive to (increasing in strength with decreasing ) and rather insensitive to variation in gravity. At a given , the band strengths change only slightly even for a large 0.5 change in gravity (in the log g =4.5-5.5 range expected for low-mass stars). At a given gravity, however, they vary significantly over a change of only 100 K in .
The surface gravity determination was checked on the width of gravity-sensitive atomic lines such as the K I and Na I doublets (see Fig. 6) as well as on the relative strength of metal-hydride bands such as CaH. The K I doublet at 7665 and 7699 and the Na I lines at 8183 and 8194 are particularly useful gravity discriminants for M dwarfs and subdwarfs. The overall line strength (central depth and equivalent width) increases with gravity as the decreasing ionization ratio due to higher electron pressure leaves more neutral alkali lines in the deeper atmosphere (Reiners, 2005). The width of the damping wings in addition increases due to the stronger pressure broadening, mainly by H, He, and H I collisions.
We also checked the metallicity determination on the spectral interval from 8440 to 8900 where molecular absorptions are lower and atomic
lines appear clearly (see Fig. 7). The synthetic spectrum represents the line profiles fairly well for elemental species such
as Ti, I Fe I, Ca II, and Mg I. The best-fit parameters (, , [Fe/H]) are given in Table 2.
|LHS 72||sdK4||3900 23||4.50.13||-1.4 0.27|
|LHS 73||sdK7||3800 39||4.50.13||-1.40.19|
|G 18-37||esdK7||3800 78||4.50.15||-1.30.44|
|APMPM J2126-4454||sdM0||3700 49||4.5 0.19||-1.30.23|
|LHS 300||sdM0||3800 39||4.50.17||-1.40.24|
|LHS 401||sdM0.5||3800 26||4.50.17||-1.40.28|
|LHS 158||sdM1||3600 48||4.50.17||-1.00.3|
|LHS 320||sdM2||3600 59||4.60.23||-0.60.31|
|LHS 406||sdM2||3600 40||4.70.24||-0.60.24|
|LHS 161||esdM2||3700 77||4.80.30||-1.20.36|
|LP 771-87||usdM2||3600 95||4.80.31||-1.40.52|
|LHS 541||sdM3||3500 76||5.10.31||-1.00.39|
|LHS 272||sdM3||3500 66||5.20.30||-0.70.37|
|LP 707-15||esdM3||3500 68||5.50.29||-0.50.36|
|LSR J1755+1648||sdM3.5||3400 52||5.10.31||-0.50.33|
|LHS 375||esdM4||3500 79||5.50.32||-1.10.31|
|LHS 1032||usdM4.5||3300 63||4.50.32||-1.70.25|
|LHS 377||sdM7||3100 32||5.30.25||-1.00.16|
|APMPM J0559-2903||esdM7||3200 68||5.40.34||-1.70.25|
|SSPM J1013-1356||sdM9.5||3000 0||5.50.05||-1.10.16|
The spectral type derived from Gizis (1997) is probably too late because of an observational problem that affects the extreme red slope of the spectrum (J. Gizis, private communication). This agrees with the warmer temperature we derived for this object.
The effective temperature versus near-infrared colours are shown in Fig. 8. The expected relations from evolution models (Baraffe et al., 1997, 1998) assuming an age of 10 Gyr and varying metallicities are also superimposed. The colours stand for different metallicities. The plot shows discrepancies between the models and our observations. It shows that the metallicities determined from the high-resolution spectral features are inconsistent with those we would infer from near-infrared colours. This inconsistency may be due to uncertainties in the CIA opacities or to outdated model interiors. It shows that broad-band colours are not sufficient to determine the parameters of M subdwarfs.
The relation of effective temperature versus spectral-type is shown in Fig. 9. The relation determined using UVES sample is compared with the scale of M dwarfs determined by Rajpurohit et al. (2013). The of subdwarfs is 200-300 K higher than of M dwarfs for the same spectral type except for hot temperatures. This is expected since the TiO bands are depleted with decreasing metallicity, and as a result, the pseudo-continuum is brighter and the flux is emitted from the hot deeper layer. A comparison with the earlier work from Gizis (1997) is also shown. Gizis (1997) determined the temperature by comparing the low-resolution optical spectra of a sample of sdM and esdM with the NextGen model atmosphere grid by Allard et al. (1997). The for subdwarfs agrees within 100 K. This difference is due to the incompleteness of the TiO and water vapor line lists used in the NextGen model atmospheres compared with the new BT-Settl models. Furthermore, this work allows us to extend the relation to the coolest M subdwarfs.
We compiled from the literature the spectral indices of TiO5, CaH2, and CaH3 computed from TiO and CaH band strengths on low-resolution spectra (see Table 1). The last column gives the references for spectral types and spectral indices. When possible, we took both the spectral type and the indices from the same paper. We note that the spectral indices may vary over several subtypes for some of the stars from one author to the other. Since Jao et al. (2008) did not give new spectral types in the scheme developed by Gizis (1997) that was extended by Lépine et al. (2007) (sd, esd, usd), which we used, we adopted the indices from Jao et al. (2008) and assigned sdM types only when no indices were available from the papers providing an sd, esd, usd classification.
Fig. 10 shows the CaH2+CaH3 versus TiO5. Lépine et al. (2003b) showed that such a diagram is useful to distinguish between the different object classes, sdM, esdM, and usdM. Our metallicity determinations are labelled in the diagram. It shows that the metallicity decreases from sdM to usdM as expected.
We also derived the parameter as defined by Lépine et al. (2007). This parameter is a combination of the TiO, CaH, and CaH spectral indices. Woolf et al. (2009) used a sample of 12 esdM and usdM with known metallicities to derive a correlation between and the metallicity, and showed that it can be used as a metallicity indicator. We plot the metallicity of our stars versus (Fig. 11) and superimpose the sample from Woolf et al. (2009) as well as their derived correlation (dashed line). We also define a similar correlation (solid line) using our sample,
Although the relation shows a strong dispersion in both samples, our metallicity determinations tend to be lower on average. This can be due to the different spectral indices adopted, or to the use of different atmosphere models. LHS 161, the only common star in both samples, illustrates the first assumption. Our metallicity determination (-1.2 dex) agrees with that of Woolf et al. (2009), (-1.3 dex) but the spectral indices adopted are very different: we used from the indices of Gizis (1997) whereas Woolf et al. (2009) used from the indices of Jao et al. (2008). Furthermore, Woolf et al. (2009) used NextGen models.
We presented a high-resolution optical spectral atlas for a sample of 18 M subdwarfs and 3 K subdwarfs, including 5 esdK,1 esdM and 2 usdM555The atlas will be made public through the Virtual Observatory tools.. We described various atomic and molecular features that appear in the spectra and their evolution with decreasing effective temperature and metallicity. We used the most recent BT-Settl model atmospheres, with revised solar abundances, to determine the scaled solar abundances of the subdwarfs. We compared them with the synthetic spectra produced by the BT-Settl atmosphere models and derived their fundamental stellar parameters. The accuracy of the atmospheric models involved in the metallicity determination can be inferred by looking at the fit to the individual atomic and molecular lines. These high-resolution spectra allowed us to separate the atmospheric parameters (effective temperature, gravity, metallicity), which is not possible when using broad-band photometry.
We determined the relation of effective temperature versus spectral type of M-subdwarfs and compared it with the previous study from Gizis (1997). Our relation agrees within 100 K and extends to the cooler spectral sequence. This work also contributes to calibrating the relation between metallicity and photometric colours and molecular band strengths. With calibration, it will be possible to estimate the metallicity of a large sample of subdwarfs, and obtain a meaningful statistical analysis.
The next step of this work is to investigate the effect of different alpha enhancements on the spectral energy distribution and the spectral-line strengths of M subdwarfs. This effect is expected to be as important as the revised oxygen abundances were to make the models fit solar metallicity M dwarfs. This is a limitation of the models because M subdwarfs are divided into two regimes, each with a rough prescription of alpha enhancements. More realistic abundance trends need to be applied to the models (Neves et al., 2009; Adibekyan et al., 2012). This atlas is appropriate for testing these prescriptions and their effect on the atmosphere spectra because they are precise enough to allow us to derive elemental abundances.
Acknowledgements.We acknowledge observing support from the ESO staff. We acknowledge financial support from "Programme National de Physique Stellaire" (PNPS) of CNRS/INSU, France. We thank the referee J. Gizis for his useful comments on the paper.
- Abel et al. (2011) Abel, M., Frommhold, L., Li, X., & Hunt, K. L. C. 2011, in 66th International Symposium On Molecular Spectroscopy
- Adibekyan et al. (2012) Adibekyan, V. Z., Sousa, S. G., Santos, N. C., et al. 2012, A&A, 545, A32
- Allard (1990) Allard, F. 1990, PhD thesis, PhD thesis. Ruprecht Karls Univ. Heidelberg, (1990)
- Allard & Hauschildt (1995) Allard, F. & Hauschildt, P. H. 1995, ApJ, 445, 433
- Allard et al. (1997) Allard, F., Hauschildt, P. H., Alexander, D. R., & Starrfield, S. 1997, ARA&A, 35, 137
- Allard et al. (2001) Allard, F., Hauschildt, P. H., Alexander, D. R., Tamanai, A., & Schweitzer, A. 2001, ApJ, 556, 357
- Allard et al. (2012) Allard, F., Homeier, D., & Freytag, B. 2012, Royal Society of London Philosophical Transactions Series A, 370, 2765
- Allard et al. (2013) Allard, F., Homeier, D., Freytag, B., et al. 2013, Memorie della Societa Astronomica Italiana Supplementi, 24, 128
- Allard et al. (2007) Allard, N. F., Kielkopf, J. F., & Allard, F. 2007, European Physical Journal D, 44, 507
- Alloin & Bica (1989) Alloin, D. & Bica, E. 1989, A&A, 217, 57
- Asplund et al. (2009) Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481
- Bailey & Kedziora-Chudczer (2012) Bailey, J. & Kedziora-Chudczer, L. 2012, MNRAS, 419, 1913
- Baraffe et al. (1997) Baraffe, I., Chabrier, G., Allard, F., & Hauschildt, P. H. 1997, A&A, 327, 1054
- Baraffe et al. (1998) Baraffe, I., Chabrier, G., Allard, F., & Hauschildt, P. H. 1998, A&A, 337, 403
- Barber et al. (2006) Barber, R. J., Tennyson, J., Harris, G. J., & Tolchenov, R. N. 2006, MNRAS, 368, 1087
- Barklem et al. (2000) Barklem, P. S., Piskunov, N., & O’Mara, B. J. 2000, A&A, 363, 1091
- Bean et al. (2006a) Bean, J. L., Sneden, C., Hauschildt, P. H., Johns-Krull, C. M., & Benedict, G. F. 2006a, ApJ, 652, 1604
- Bean et al. (2006b) Bean, J. L., Sneden, C., Hauschildt, P. H., Johns-Krull, C. M., & Benedict, G. F. 2006b, ApJ, 652, 1604
- Bidelman (1985) Bidelman, W. P. 1985, ApJS, 59, 197
- Bonfils et al. (2005) Bonfils, X., Forveille, T., Delfosse, X., et al. 2005, A&A, 443, L15
- Borysow & Frommhold (1989) Borysow, A. & Frommhold, L. 1989, ApJ, 341, 549
- Borysow et al. (2001) Borysow, A., Jørgensen, U. G., & Fu, Y. 2001, Journal of Quantitative Spectroscopy and Radiative Transfer, 68, 235
- Borysow et al. (1997) Borysow, A., Jorgensen, U. G., & Zheng, C. 1997, A&A, 324, 185
- Burgasser & Kirkpatrick (2006) Burgasser, A. J. & Kirkpatrick, J. D. 2006, ApJ, 645, 1485
- Burgasser et al. (2003) Burgasser, A. J., Kirkpatrick, J. D., Burrows, A., et al. 2003, ApJ, 592, 1186
- Caffau et al. (2011) Caffau, E., Ludwig, H.-G., Steffen, M., Freytag, B., & Bonifacio, P. 2011, Sol. Phys., 268, 255
- Casagrande et al. (2008) Casagrande, L., Flynn, C., & Bessell, M. 2008, MNRAS, 389, 585
- Chowdhury et al. (2006) Chowdhury, P. K., Merer, A. J., Rixon, S. J., Bernath, P. F., & Ram, R. S. 2006, Physical Chemistry Chemical Physics (Incorporating Faraday Transactions), 8, 822
- Dawson & De Robertis (2000) Dawson, P. C. & De Robertis, M. M. 2000, AJ, 120, 1532
- Deen (2013) Deen, C. P. 2013, AJ, 146, 51
- Dekker et al. (2000) Dekker, H., D’Odorico, S., Kaufer, A., Delabre, B., & Kotzlowski, H. 2000, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 4008, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, ed. M. Iye & A. F. Moorwood, 534–545
- Digby et al. (2003) Digby, A. P., Hambly, N. C., Cooke, J. A., Reid, I. N., & Cannon, R. D. 2003, MNRAS, 344, 583
- Dulick et al. (2003) Dulick, M., Bauschlicher, Jr., C. W., Burrows, A., et al. 2003, ApJ, 594, 651
- Freytag et al. (2012) Freytag, B., Steffen, M., Ludwig, H.-G., et al. 2012, J. Comp. Phys., 231, 919
- Gamache et al. (1996) Gamache, R. R., Lynch, R., & Brown, L. R. 1996, J. Quant. Spec. Radiat. Transf., 56, 471
- Gizis (1997) Gizis, J. E. 1997, AJ, 113, 806
- Goorvitch & Chackerian (1994a) Goorvitch, D. & Chackerian, Jr., C. 1994a, ApJS, 91, 483
- Goorvitch & Chackerian (1994b) Goorvitch, D. & Chackerian, Jr., C. 1994b, ApJS, 92, 311
- Gustafsson et al. (2008) Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951
- Gustafsson & Frommhold (2001) Gustafsson, M. & Frommhold, L. 2001, ApJ, 546, 1168
- Gustafsson et al. (2003) Gustafsson, M., Kristensen, L. E., Clénet, Y., et al. 2003, A&A, 411, 437
- Hauschildt et al. (1997) Hauschildt, P. H., Baron, E., & Allard, F. 1997, ApJ, 483, 390
- Jao et al. (2008) Jao, W.-C., Henry, T. J., Beaulieu, T. D., & Subasavage, J. P. 2008, AJ, 136, 840
- Johnas et al. (2007) Johnas, C. M. S., Hauschildt, P. H., Schweitzer, A., et al. 2007, A&A, 475, 1039
- Kirkpatrick et al. (2010) Kirkpatrick, J. D., Looper, D. L., Burgasser, A. J., et al. 2010, ApJS, 190, 100
- Lépine et al. (2003a) Lépine, S., Rich, R. M., & Shara, M. M. 2003a, AJ, 125, 1598
- Lépine et al. (2007) Lépine, S., Rich, R. M., & Shara, M. M. 2007, ApJ, 669, 1235
- Lépine et al. (2003b) Lépine, S., Shara, M. M., & Rich, R. M. 2003b, AJ, 126, 921
- Ludwig et al. (2002) Ludwig, H.-G., Allard, F., & Hauschildt, P. H. 2002, A&A, 395, 99
- Ludwig et al. (2006) Ludwig, H.-G., Allard, F., & Hauschildt, P. H. 2006, A&A, 459, 599
- Mallik (1997) Mallik, S. V. 1997, A&AS, 124, 359
- Neves et al. (2009) Neves, V., Santos, N. C., Sousa, S. G., Correia, A. C. M., & Israelian, G. 2009, A&A, 497, 563
- Önehag et al. (2012) Önehag, A., Heiter, U., Gustafsson, B., et al. 2012, A&A, 542, A33
- Plez (1998) Plez, B. 1998, A&A, 337, 495
- Rajpurohit et al. (2013) Rajpurohit, A. S., Reylé, C., Allard, F., et al. 2013, A&A, 556, A15
- Rajpurohit et al. (2012) Rajpurohit, A. S., Reylé, C., Schultheis, M., et al. 2012, A&A, 545, A85
- Reiners (2005) Reiners, A. 2005, Astronomische Nachrichten, 326, 930
- Reiners & Basri (2006) Reiners, A. & Basri, G. 2006, AJ, 131, 1806
- Reiners et al. (2007) Reiners, A., Homeier, D., Hauschildt, P. H., & Allard, F. 2007, A&A, 473, 245
- Reylé et al. (2006) Reylé, C., Scholz, R.-D., Schultheis, M., Robin, A. C., & Irwin, M. 2006, MNRAS, 373, 705
- Rodgers & Eggen (1974) Rodgers, A. W. & Eggen, O. J. 1974, PASP, 86, 742
- Rothman et al. (2009) Rothman, L. S., Gordon, I. E., Barbe, A., et al. 2009, J. Quant. Spec. Radiat. Transf., 110, 533
- Scholz et al. (2002) Scholz, R.-D., Ibata, R., Irwin, M., et al. 2002, MNRAS, 329, 109
- Scholz et al. (2004) Scholz, R.-D., Lehmann, I., Matute, I., & Zinnecker, H. 2004, A&A, 425, 519
- Schweitzer et al. (1999) Schweitzer, A., Scholz, R.-D., Stauffer, J., Irwin, M., & McCaughrean, M. J. 1999, A&A, 350, L62
- Skory et al. (2003) Skory, S., Weck, P. F., Stancil, P. C., & Kirby, K. 2003, ApJS, 148, 599
- Skrutskie et al. (2006) Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163
- Tashkun et al. (2004) Tashkun, S. A., Perevalov, V. I., Teffo, J.-L., et al. 2004, Proc. SPIE, 5311, 102
- Tinney & Reid (1998) Tinney, C. G. & Reid, I. N. 1998, MNRAS, 301, 1031
- Unsold (1968) Unsold, A. 1968, Physik der Sternatmospharen, MIT besonder Berucksichtigung der Sonne
- Valenti & Piskunov (1996) Valenti, J. A. & Piskunov, N. 1996, A&AS, 118, 595
- Weck et al. (2003) Weck, P. F., Schweitzer, A., Stancil, P. C., Hauschildt, P. H., & Kirby, K. 2003, ApJ, 584, 459
- Woolf et al. (2009) Woolf, V. M., Lépine, S., & Wallerstein, G. 2009, PASP, 121, 117
- Woolf & Wallerstein (2005) Woolf, V. M. & Wallerstein, G. 2005, MNRAS, 356, 963
- Zhou (1991) Zhou, X. 1991, A&A, 248, 367