Line-Depth Ratio in 0.97–1.32\,\text{mm}

Method to Estimate the Effective Temperatures of Late-Type Giants using Line-Depth Ratios in the Wavelength Range

Daisuke Taniguchi, Noriyuki Matsunaga, Naoto Kobayashi, Kei Fukue, Satoshi Hamano, Yuji Ikeda, Hideyo Kawakita, Sohei Kondo, Hiroaki Sameshima and Chikako Yasui
Department of Astronomy, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan
Laboratory of Infrared High-resolution spectroscopy (LiH), Koyama Astronomical Observatory, Kyoto Sangyo University, Motoyama,
Kamigamo, Kita-ku, Kyoto 603-8555, Japan
Institute of Astronomy, The University of Tokyo, 2-21-1 Osawa, Mitaka, Tokyo 181-0015, Japan
Kiso Observatory, The University of Tokyo, 10762-30 Mitake, Kiso-machi, Kiso-gun,Nagano 397-0101, Japan
Photocoding, 460-102 Iwakura-Nakamachi, Sakyo-ku, Kyoto 606-0025, Japan
Department of Astrophysics and Atmospheric Sciences, Faculty of Sciences, Kyoto Sangyo University, Motoyama, Kamigamo,
Kita-ku, Kyoto 603-8555, Japan
National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan
Accepted XXX. Received YYY; in original form ZZZ

The effective temperature, one of the most fundamental atmospheric parameters of a star, can be estimated using various methods, and here we focus on the method using line-depth ratios (LDRs). This method combines low- and high-excitation lines and makes use of relations between LDRs of these line pairs and the effective temperature. It has an advantage, for example, of being minimally affected by interstellar reddening, which changes stellar colours. We report 81 relations between LDRs and the effective temperature established with high-resolution, , spectra of nine G- to M-type giants in Y and J bands. Our analysis gives the first comprehensive set of LDR relations for this wavelength range. The combination of all these relations can be used to determine the effective temperatures of stars that have and to the precision of in the best cases.

techniques: spectroscopic – stars: late-type – stars: atmospheres – infrared: stars
pubyear: 2017pagerange: Method to Estimate the Effective Temperatures of Late-Type Giants using Line-Depth Ratios in the Wavelength Range Supporting Information

1 Introduction

Stellar atmosphere is characterized mainly by effective temperature (), surface gravity () and metallicity ([Fe/H]), of which has a particularly strong impact on the spectra and often is the first parameter to estimate during spectral analysis. Among many methods of measuring , we focus on the line-depth ratio (LDR) method in this paper. In cool stars (around the solar temperature and below), the depths of low-excitation lines of neutral atoms are sensitive to while those of high-excitation lines are relatively insensitive (Gray, 2008), and hence, their ratios are good temperature indicators (Gray & Johanson, 1991; Fukue et al., 2015, and references therein). This method has a few advantages: LDRs are not affected by interstellar reddening, and atmospheric parameters other than are expected to have no, or at least relatively weak, effects on the LDRs.

Sasselov & Lester (1990) used LDRs between \ionCi and \ionSii lines around to measure the temperatures and reddenings of Cepheids in a relative way. This pioneering work, however, had a limited impact on later applications due to the small numbers of the standard stars and line pairs. On the contrary, previous works on the LDR method mostly used visible spectral ranges(e.g. Kovtyukh et al., 2003, 2006; Kovtyukh, 2007). With the development of infrared facilities in recent years, it has become possible to make more comprehensive studies using high-resolution infrared spectra. Recently, Fukue et al. (2015) established the LDR method for H-band spectra, . In this wavelength range, however, the small number of low-excitation lines and severe blending between molecular lines and metal lines make it difficult to find a large number of useful LDR pairs. As a result, the precision of after taking an average of the values from available pairs is only for each star in the best case, while the LDR method in the visible range can achieve the precision of (Kovtyukh et al., 2006). Here we investigate Y- and J-band spectra to find the LDR– relations for the first time. This wavelength range has a number of both low- and high-excitation lines, e.g. Meléndez & Barbuy (1999), and thus one can expect many useful LDR pairs to achieve a high precision.

2 Observation

The ten targets listed in Table 1 were selected mainly from Gaia FGK benchmark stars (Blanco-Cuaresma et al., 2014; Jofré et al., 2014; Heiter et al., 2015) in addition to a few prototypes of red giants such as Arcturus. All of our targets are giants with the luminosity class of III or II, and have the spectral types between early G to early M (). Their metallicities are between and in [Fe/H] except Arcturus, .

Object HD Sp. Type [Fe/H] [dex] Obs. Date
ε Leo 84441 G1II 2014 Jan 23
ϰ Gem 62345 G8IIIa 2013 Dec 8
ε Vir 113226 G8III 2014 Jan 23
Pollux 62509 K0III 2013 Feb 28
μ Leo 85503 K2IIIb 2013 Feb 23
Alphard 81797 K3II-III 2013 Nov 30
Aldebaran 29139 K5III 2013 Feb 24
α Cet 18884 M1.5IIIa 2013 Nov 30
δ Oph 146051 M0.5III 2014 Jan 23
Arcturus 124897 K0III 2013 Feb 23

References: [1] Heiter et al. (2015) and [2] Prugniel et al. (2011).

Table 1: The target stars and their atmospheric parameters. Spectral types are from SIMBAD (Wenger et al., 2000). The temperatures of μ Leo and Aldebaran were estimated in both references (1) and (2), and we use the weighted means. The last column gives the dates of our WINERED observations.

We observed these ten stars between 2013 February 23 and 2014 January 23 using WINERED attached to a Nasmyth platform of the Araki Telescope at Koyama Astronomical Observatory of Kyoto Sangyo University in Japan (Ikeda et al., 2016). WINERED can collect spectra covering the wavelength range from to (z, Y and J bands) with the spectral resolution of with one integration. All of our targets are bright, , and the total integration time with each target within the slit was between and to achieve the signal-to-noise ratio (S/N) of or higher. We describe the resultant S/N in subsection 3.1. For every target star, we also observed a telluric standard star (an A0V star in most cases) for subtracting telluric absorption.

3 Data Analysis

3.1 Spectra Reduction

Basic data reduction was performed automatically using the pipeline software prepared by the WINERED team based on pyraf111pyraf is a prodict of the Space Telescope Science Institute, which is operated by AURA for NASA.. The pipeline process includes the standard analysis steps for echelle spectra: bad pixel masking, sky subtraction, flat fielding, scattered light subtraction, transformation of each two-dimensional echelle image into images with the space and wavelength axes orthogonal to each other for individual orders, spectrum extraction, wavelength calibration and continuum normalisation. In addition, it was necessary to correct for time-dependent wavelength shifts as follows. We used ThAr lamp data as the initial wavelength calibration in the pipeline software, but we found that some spectra had wavelength offsets probably caused by varying ambient temperature222In late 2016, we made an instrumental upgrade to minimize such offsets.. In this work, such offsets are not critical because we are not interested in radial velocities, but the adjustment of the wavelength scale of individual frames was necessary to avoid the artificial broadening of the line spread function in a combined spectrum of each target. The relative offsets between individual spectra were measured using almost isolated telluric absorption lines and were removed before combining the spectra. Then, telluric absorption lines in the combined spectrum of each target were removed using a spectrum of a telluric-standard star except the 53rd to 55th orders, to , in which almost no significant telluric lines are present. The wavelength scale of each telluric-corrected spectrum was converted to the rest frame of each star using its intrinsic absorption lines. Our wavelengths are in the standard air scale. Finally, we re-normalise the spectra using pyraf to adjust the continuum level to the unity. In the following discussions, we consider 12 orders listed in Table 2. We ignore other orders (42nd, 49th to 51st and 58th to 61st) which are around gaps between atmospheric windows. In the analysis to find and calibrate the LDR– relations, we ignore Arcturus because, according to Fukue et al. (2015), its low metallicity causes offsets from the trends between the LDR and the effective temperature. Table 2 shows the S/N of each order of each star calculated from the variance of the normalised counts in continuum regions, as an indicator of spectral quality. The spectral resolution was stable and quite constant for all of our spectra. The FWHM values measured with most stellar lines correspond to between and . Figure 1 shows a small range of the obtained spectra. This range has a line pair which shows a clear temperature dependency as we discuss below.

Object 57 56 55 54 53 52 48 47 46 45 44 43
ε Leo
ϰ Gem
ε Vir
μ Leo
α Cet
δ Oph
Table 2: S/N of the reduced spectra in each echelle order (52–57th in Y band and 43–48th in J band).
Figure 1: A part of the spectra in the 55th order are drawn in the order of decreasing effective temperature (from top) except Arcturus in the bottom which were not included in establishing the LDR– relations. The two lines in the graph show a clear dependency on the temperature (the high-excitation line on the left and the low-excitation line on the right).

3.2 Line Selection and Measurement of Line Depth

As the candidates of atomic lines to be used for LDRs, we consider the lines identified in our Arcturus spectrum. Here we outline its line identification, but details of the analysis and the entire line list will be published in a forthcoming paper (Ikeda et al. in preparation). Firstly, we produced a synthesized spectrum of Arcturus, whose parameters are assumed to be , , and (Smith et al., 2013) using the ATLAS9 atmospheric models and codes given by Kurucz (1993). As the input list of atomic lines for the synthetic spectrum, Vienna Atomic Line Database (VALD3; Ryabchikova et al., 2015) were used. Secondly, the observed and synthetic spectra were closely compared with each other to identify lines whose depth is larger than . Although the S/N of the Arcturus spectrum is higher than in all wavelengths, CN molecular lines often make the identification of shallow lines less certain. In such cases, to confirm the contribution of weak atomic lines around molecular lines, we further compared the observed spectrum with a synthetic spectrum but with only CN lines included. Thus we obtained approximately 800 atomic lines.

Among the lines identified in Arcturus, we here include relatively isolated neutral lines of Fe, Ti, Si, Cr, Ca, Ni, Mg, Na, Co, Al, Mn and K in our analysis for LDRs. Ion lines and lines of C, N and O were not included (Kovtyukh et al., 2006). Moreover, we excluded lines by the following criteria: deep lines which appear deeper than in Arcturus, and lines blended with molecular CN lines or other atomic lines which have different elements and/or excitation potentials. These criteria left 125 and 99 lines in Y and J bands, respectively. We refer to VALD3 for parameters of the lines such as excitation potential and oscillator strength.

For the 125 and 99 lines we selected, we used a quadratic function to fit three (or four) pixels near the bottom of each line rather than fitting the entire line profile, e.g., by the Gaussian or the Voigt function (Strassmeier & Schordan, 2000). We define as a line depth from the continuum level to the bottom of the fitted function, where indicates an ID number of the line and indicates an ID number of a star. Weak lines with are ignored in the following analysis.

3.3 Line Pair Selection

We calculated an LDR, , for any line pair, , of approximately 1500 pairs whose lines were both detected in more than four stars and have excitation potentials separated by more than . We treated each order independently, i.e. did not combine lines in different orders, and divided the depth of the low-excitation line, , by that of the high-excitation line, . For each line pair, , we plotted the effective temperatures against the common logarithms of LDRs, , and determined the regression line, , using the Weighted Total Least Squares method (see Markovsky & Van Huffel, 2007, for a review). The weight of a line pair, , of a star was calculated by , where and indicate the standard error of calculated by the S/N and that of literature for each star, respectively. The dispersion around each regression line is defined as . The effective temperature and its error of each star based on each relation are given as , of which the error was determined by variance–covariance matrix of the coefficients of the regression line and the error of . In the subsequent analysis, about 900 line pairs which have were not included. Moreover, 21 line pairs which have were excluded. A large fraction of the rejected pairs includes a weak iron line at () in Y band. This line may be blended by other line(s), which causes the unexpected dependency of LDR on the effective temperature, though the presence of such blended lines is not clear in available line lists. The rest of the line pairs with have \ionTii lines for both low- and high-excitation potentials 333We found relatively tight LDR– relations with the unexpected slopes for six pairs of \ionTii lines (for each pair, the first is the low-excitation line): , , , , and .. LDR– relations for the majority of \ionTii–\ionTii pairs are not tight enough to be selected, and their slopes can be negative, positive or not well determined. The reason of the opposite slopes for the small number of \ionTii–\ionTii pairs is unclear.

There is a large number of possible combinations of line pairs, and we selected a set of pairs as follows. One condition is that each line is used only once, i.e. not included in more than one line pairs. This ensures that statistical errors in values in different LDR– relation are independent of each other and makes it easy to calculate the statistical errors of combined effective temperatures ( in section 4). In fact, the final set of the relations we discuss in section 4 makes use of 162 lines, which is about of the selected lines. This condition for the line pair selection has no large impact on the statistical error, which can be better only by even if we use all the pre-selected pairs.

The basic idea of the process is to select line pairs which meet the following conditions as much as possible: (i) high precision in reproducing the effective temperatures of our sample stars, and (ii) small difference in wavelength between two lines of each line pair. With the latter condition fulfilled, the possibility for other instruments to detect both lines of a line pair in the same echelle order gets high and the error of the LDR introduced in the continuum normalisation is expected to be smaller.

First, we consider the undirected graph whose nodes correspond to absorption lines and edges connecting the nodes correspond to line pairs. Note that our analysis treats each order independently. The goal is to find the optimal matching, , that meets the above conditions for each order. We consider maximum matchings, the number of whose edges is as large as possible but no node is connected by more than one edges. In an ideal case, the size of a maximum matching, , corresponds to half the number of all the nodes in the original undirected graph, but this is not true in our case because many edges were rejected owing to, for example, large scatters around the corresponding LDR– relations. For a maximum matching, , of this undirected graph, each line pair of was named . For a given matching, each star’s effective temperature was redetermined as a weighted mean of and was named . The statistical error in was estimated by a weighted standard deviation of . Then, let be the difference in wavelength between two lines of a line pair, , and we consider the evaluation function as


represents the size of the error in redetermining the effective temperatures of the nine stars, without Arcturus included, for a given matching , while represents the wavelength difference of the line pairs and works as a penalty term. For a given value, there are different allowed combinations of line pairs which form different maximum matchings, and we select the one which gives the least as the optimal matching, . The coefficient will be determined, in the next section, by considering how and depend on it.

4 Result

We applied the procedure described above to our 224 nodes (i.e. lines) and 603 edges (i.e. line pairs). With changing the coefficient from to , we searched for the optimized matching, , which gives the smallest at each value and observed how affects the solutions. Figure 2 plots the values of and for with varying . By increasing , the weight of increases relative to in the evaluation function, and then the optimal matching and the values of and change. At around , is only slightly larger compared to the case with for most orders (the exception being the order 48) which optimizes the precision in the redetemined temperature by ignoring the difference in the wavelength, while improves to some extent by changing from to . The number of line pairs in the 48th order is only 4, and with the wavelength separations of the individual pairs get significantly small with only a small increase in . We selected , and Table 3 lists the resultant values of the evaluation functions and other parameters in our analysis.

Figure 2: The dependence of and on the parameter, , in Equation 1. The evaluation functions are defined in the text. Variations of (top row) and (bottom row) are presented for a group of echelle orders in each band (52nd to 57th in Y band on the left and 43rd to 48th in J band on the right). Each solid line represents the evaluation function and the vertical dashed line represents the value which we used for the final solutions. Legends show which colour corresponds to which echelle order. The colour version of this figure is available in the online journal.
Order 57 56 55 54 53 52 48 47 46 45 44 43 Combined
4 23 19 30 17 32 13 11 19 14 16 26 224
1 10 4 10 8 14 3 3 4 5 7 12 81
Table 3: Parameters of the selected line pairs in individual orders (52–57th in Y and 43–48th in J); the wavelength range from the WINERED web site (, the number of lines considered to make pairs, the number of line pairs which we included in the final solutions and the values of the evaluation functions. In addition, the last column gives the same parameters but for all the orders combined.

We obtained 47 LDR– relations in Y band and 34 in J band with . Table 4 and 5 list the line pairs and the relations, and Figure 3 shows some examples of the relations. We redetermined the effective temperatures using our final relations ( in Table 6) and compared them with in the literature in Figure 4. The difference, , has no clear dependency on or [Fe/H], suggesting that this LDR method is effective across the parameter range covered by the nine calibrating stars, i.e. and . In order to evaluate the power of the LDR method considering all the orders available, we consider the evaluation function of Equation 1 but including all the orders and the line pairs: i.e. we use instead of the temperatures for an individual order, , and consider the wavelength differences of all the line pairs. Thus calculated evaluation functions are and . These values are, for example, smaller than those444Fukue et al. (2015) used no evaluation function like ours, but we used their published data to calculate the values according to Equation 1. for the result by Fukue et al. (2015), and , who used nine line pairs in H band.

It is worthwhile to discuss for Arcturus, which was not included in the calibration. Fukue et al. (2015) found that Arcturus tends to show offsets from LDR– relations of other stars. They suggested that this is owing to the low metallicity and/or non-solar abundance ratios of Arcturus, , , by Jofré et al. (2015). The result in Table 6 appears to show that the temperature from our method () is consistent with the literature , but a closer look reveals the peculiarity of Arcturus’s . If we use only 22 line pairs composed of \ionTii (low-excitation) and \ionFei (high-excitation), we obtain the temperature of . If we use only 14 line pairs of \ionSii (high-excitation) and \ionFei (low-excitation), in contrast, the temperature is estimated at . These deviations can be understood as the effect of non-solar [Ti/Fe] and [Si/Fe] values. These offsets of the two major groups of line pairs cancelled out each other, and the average became close to the literature temperature , but giving a relatively large error, in the case of Arcturus. Such deviations are not observed in other objects in this study. Follow-up spectroscopic data for a larger sample of stars with known temperatures and abundances would be useful to further discuss the abundance dependences of the LDR– relations.

Low-excitation Line High-excitation Line LDR– Relation
ID Order [Å] El. E.P. [eV] [Å] El. E.P. [eV] [K] [K] flag
(1) 57
(2) 56 \ionTii \ionFei
(3) 56 \ionCri \ionFei 1
(4) 56 \ionTii \ionNai 1
(5) 56 \ionFei \ionFei
(6) 56 \ionTii \ionFei
(7) 56 \ionTii \ionFei
(8) 56 \ionTii \ionFei
(9) 56 \ionTii \ionFei
(10) 56 \ionTii \ionSii
(11) 56 \ionTii \ionFei 1
(12) 55 \ionFei \ionFei
(13) 55 \ionFei \ionNii
(14) 55 \ionTii \ionFei
(15) 55 \ionCai \ionFei
(16) 54 \ionFei \ionSii 1
(17) 54 \ionFei \ionSii 1
(18) 54 \ionNii \ionSii 1
(19) 54 \ionFei \ionFei 1
(20) 54 \ionTii \ionFei 1
(21) 54 \ionCai \ionSii 1
(22) 54 \ionFei \ionFei 1
(23) 54 \ionCoi \ionFei 1
(24) 54 \ionFei \ionFei 1
(25) 54 \ionTii \ionFei 1
(26) 53 \ionCri \ionNii 1
(27) 53 \ionTii \ionFei 1
(28) 53 \ionCri \ionFei 1
(29) 53 \ionTii \ionFei 1
(30) 53 \ionTii \ionSii 1
(31) 53 \ionFei \ionFei 1
(32) 53 \ionTii \ionCri 1
(33) 53 \ionFei \ionSii 1
(34) 52 \ionTii \ionFei
(35) 52 \ionFei \ionSii
(36) 52 \ionFei \ionSii
(37) 52 \ionNai \ionFei
(38) 52 \ionFei \ionFei
(39) 52 \ionTii \ionFei
(40) 52 \ionFei \ionSii
(41) 52 \ionTii \ionNii
(42) 52 \ionFei \ionSii
(43) 52 \ionTii \ionCai
(44) 52 \ionFei \ionSii
(45) 52 \ionNai \ionSii
(46) 52 \ionCri \ionMgi
(47) 52 \ionCri \ionSii
Table 4: List of low- and high-excitation lines and the LDR– relations in Y band. and represent the coefficients of , represents the number of stars used in fitting and represents the dispersion of the relation. If more than one lines are listed at the same wavelength in the original line catalogue, information of the line with the largest is adopted here. The flag in the last column indicates the line pairs that are used in the analysis without the telluric correction in section 5.
Low-excitation Line High-excitation Line LDR– Relation
ID Order [Å] El. E.P. [eV] [Å] El. E.P. [eV] [K] [K] flag
(48) 48 \ionFei \ionSii
(49) 48 \ionTii \ionFei
(50) 48 \ionTii \ionCai
(51) 47 \ionFei \ionSii
(52) 47 \ionCai \ionFei
(53) 47 \ionCri \ionMgi
(54) 46 \ionCai \ionSii
(55) 46 \ionTii \ionSii
(56) 46 \ionFei \ionSii
(57) 46 \ionFei \ionFei
(58) 45 \ionFei \ionSii 1
(59) 45 \ionTii \ionFei 1
(60) 45 \ionFei \ionMgi
(61) 45 \ionCri \ionMgi
(62) 45 \ionFei \ionSii
(63) 44 \ionTii \ionCoi
(64) 44 \ionTii \ionFei
(65) 44 \ionTii \ionNai
(66) 44 \ionTii \ionFei
(67) 44 \ionTii \ionFei
(68) 44 \ionTii \ionMgi
(69) 44 \ionTii \ionFei
(70) 43 \ionCri \ionFei
(71) 43 \ionCri \ionCai
(72) 43 \ionTii \ionFei
(73) 43 \ionTii \ionCai
(74) 43 \ionNii \ionFei
(75) 43 \ionCri \ionFei
(76) 43 \ionMni \ionSii
(77) 43 \ionTii \ionFei
(78) 43 \ionTii \ionSii
(79) 43 \ionFei \ionSii
(80) 43 \ionTii \ionAli
(81) 43 \ionCai \ionSii
Table 5: Same as Table 4 but for line pairs in J band.
Figure 3: Examples of the LDR– relations. The ID of the line pair is given in each panel. Arcturus indicated by a green open square in each panel was not used for the relation which was obtained with the other objects indicated by blue filled circles. Plots for all the 81 relations are available as online material — see Surpporting Information.
Figure 4: The differences between the redetermined effective temperatures by our LDR method and the literature temperatures are plotted against (left panel) and [Fe/H] (right panel). Symbols are same as in Figure 3. The size of the vertical error bars is dominated by the error in the literature temperature except for high-temperature stars (ε Leo and ϰ Gem) and Arcturus.

5 Discussion

It is useful to know what spectral quality is needed to measure the effective temperature based on our LDR method for planning future observations. Here, we discuss how the precision of depends on the S/N of a spectrum and on whether or not we can make the telluric correction.

First, considering a common S/N fixed for all orders, we calculate a feasible precision by , where is the temperature error for each star obtained in the same way described in section 3 except that the given S/N is used instead of the S/N obtained from our observed spectrum. It is straightforward to calculate the as a function of S/N, and we can estimate with which S/N the precision reaches two target values, and . As illustrated in Figure 5, in order to estimate the effective temperature of within the precision of , for example, one needs to obtain spectra with S/N higher than . Although required S/N values depend also on the metallicity and the resolution which change the line depths, Figure 5 can be used as a guideline for estimating the necessary observational times for solar-metal stars approximately using the instrument with the resolution . The statistical error shall get larger if a limited wavelength coverage allows one to use a smaller number of lines.

Figure 5: The value of S/N required for estimating with specific precision: (open diamond) and (cross). See more details in the text.

Next, we compare the precision which can be achieved with telluric-corrected or uncorrected spectra. We found that all the metal lines, used in this work, in the 53–55th orders and also a significant number of lines in other orders are free from blending with telluric lines for our objects. When the LDR method is applied to other stars, however, one needs to take radial velocities into account for considering the blending. For example, some of the lines in the 53rd to 55th orders could have been blended with telluric lines if our targets were offset by a particular redshift. For 23 line pairs marked in Table 4 and Table 5 (21 in Y band and two in J band), both lines of each pair have no telluric absorption line deeper than with the spectral range corresponding to the redshift within . We considered the absorption depths of a telluric standard spectrum of ο Aur taken on 2013 November 30. Table 6 compares the redetermined effective temperatures obtained with all line pairs in telluric-corrected spectra () and the values obtained with the above 23 line pairs in telluric-uncorrected spectra (). The error propagation was calculated in the same way as for as mentioned above and S/N values in Table 2 were used also for telluric-uncorrected spectra. We found that are consistent with . Even with significantly fewer line pairs, the precision of is still moderately high. Our LDR method can be used even without the telluric correction with slightly lower precision, which may give efficient access to stellar temperatures with less observational times and analytical works.

ε Leo (42) (14)
ϰ Gem (60) (19)
ε Vir (65) (20)
Pollux (73) (22)
μ Leo (80) (23)
Alphard (79) (23)
Aldebaran (78) (21)
α Cet (79) (21)
δ Oph (79) (21)
Arcturus (71) (18)
Table 6: The literature values of atmospheric parameters, the redetermined effective temperatures using all the available line pairs () and the ones based on uncorrected spectra (). The numbers in round brackets indicate the numbers of line pairs used for the LDR temperatures. The literature is same as Table 1.

The accuracy of the continuum normalisation is another important factor which may affect the accuracy of the LDR measurements. For example, Jofré et al. (2017) found that offsets between re-normalised continuum levels obtained with different methods are as large as in the worst case for optical lines of Arcturus (see their Fig. 4). If we add (or subtract) such a large offset to the depth of each line, the LDR temperatures could change by . However, the effect of the normalisation is expected to be smaller in the LDR method because the lines selected in this work tend to be in relatively uncrowded spectral parts for which the continuum is easily traced. In addition, the LDR– relations are based on the empirical calibration, and the line depths are not compared with a theoretical prediction in which the continuum is perfectly defined. Even if there were some offsets between the true and apparent continuum levels, the LDR values for stars with similar temperatures should be similar. In fact, the statistics of the LDR– relations we obtained indicates that the continuum normalisation has little impact on the precision of our LDRs. The scatters around individual LDR– relations are , and they can be explained by statistical errors according to the spectral quality, S/N, and the precision of the reference temperatures. Nonetheless, when LDRs are measured in spectra with a significantly different resolution, the continuum normalisation may cause systematic temperature offsets of the order of in the worst cases and the re-calibration with a homogeneous dataset is recommended.

6 Summary and Conclusion

We developed the method to estimate the effective temperature using the LDR in Y and J bands. This method enables us to estimate the effective temperatures of G, K and M giants with the precision of in the best cases using 81 line pairs. Roughly speaking, the S/N of is needed to estimate the effective temperature to the precision of for a solar-metal star with , while without the telluric correction the error would get two times larger. Although our set of line pairs is optimized for use with WINERED in terms of the resolution and the wavelength range of the individual orders, this set should also be useful for other instruments. Re-calibration of the LDR– relations would be useful for other instruments although LDRs are expected to be insensitive to the spectral resolution under no blending effects of neighbouring lines.


We acknowledge the useful comments from the referee, Ulrike Heiter. We are greatful to Masaomi Kinoshita, Kenshi Nakanishi, Tetsuya Nakaoka and Yoshiharu Shinnaka for observing a part of our targets. We also thank the staff of Koyama Astronomical Observatory for their support during our observation. This work has been supported by Grants-in-Aid, KAKENHI, from Japan Society for the Promotion of Science (JSPS; Nos. 16684001, 20340042, 21840052 and 26287028) and MEXT Supported Programs for the Strategic Research Foundation at Private Universities (Nos. S0801061 and S1411028). NK also acknowledges support through the Japan–India Scinece Cooperative Program between 2013 and 2018 under agreement between the JSPS and the Department of Science and Technology (DST) in India. KF is grateful to KAKENHI (16H07323) Grant-in-Aid for Research Activity Start-up.


  • Blanco-Cuaresma et al. (2014) Blanco-Cuaresma S., Soubiran C., Jofré P., Heiter U., 2014, A&A, 566, A98
  • Fukue et al. (2015) Fukue K., et al., 2015, ApJ, 812, 64
  • Gray (2008) Gray D. F., 2008, The Observation and Analysis of Stellar Photospheres
  • Gray & Johanson (1991) Gray D. F., Johanson H. L., 1991, PASP, 103, 439
  • Heiter et al. (2015) Heiter U., Jofré P., Gustafsson B., Korn A. J., Soubiran C., Thévenin F., 2015, A&A, 582, A49
  • Ikeda et al. (2016) Ikeda Y., et al., 2016, in Ground-based and Airborne Instrumentation for Astronomy VI. p. 99085Z, doi:10.1117/12.2230886
  • Jofré et al. (2014) Jofré P., et al., 2014, A&A, 564, A133
  • Jofré et al. (2015) Jofré P., et al., 2015, A&A, 582, A81
  • Jofré et al. (2017) Jofré P., et al., 2017, A&A, 601, A38
  • Kovtyukh (2007) Kovtyukh V. V., 2007, MNRAS, 378, 617
  • Kovtyukh et al. (2003) Kovtyukh V. V., Soubiran C., Belik S. I., Gorlova N. I., 2003, A&A, 411, 559
  • Kovtyukh et al. (2006) Kovtyukh V. V., Soubiran C., Bienaymé O., Mishenina T. V., Belik S. I., 2006, MNRAS, 371, 879
  • Kurucz (1993) Kurucz R., 1993, ATLAS9 Stellar Atmosphere Programs and 2 km/s grid. Kurucz CD-ROM No. 13.  Cambridge, Mass.: Smithsonian Astrophysical Observatory, 1993., 13
  • Markovsky & Van Huffel (2007) Markovsky I., Van Huffel S., 2007, Signal processing, 87, 2283
  • Meléndez & Barbuy (1999) Meléndez J., Barbuy B., 1999, ApJS, 124, 527
  • Prugniel et al. (2011) Prugniel P., Vauglin I., Koleva M., 2011, A&A, 531, A165
  • Ryabchikova et al. (2015) Ryabchikova T., Piskunov N., Kurucz R. L., Stempels H. C., Heiter U., Pakhomov Y., Barklem P. S., 2015, Phys. Scr., 90, 054005
  • Sasselov & Lester (1990) Sasselov D. D., Lester J. B., 1990, ApJ, 360, 227
  • Smith et al. (2013) Smith V. V., et al., 2013, ApJ, 765, 16
  • Strassmeier & Schordan (2000) Strassmeier K. G., Schordan P., 2000, Astronomische Nachrichten, 321, 277
  • Wenger et al. (2000) Wenger M., et al., 2000, A&AS, 143, 9

Supporting Information

Additional Supporting Information may be found in the online version of this article.
Figure 3. The LDR– plots for all the 81 relations.

Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description