Issue 
A&A
Volume 575, March 2015



Article Number  A26  
Number of page(s)  18  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201424310  
Published online  16 February 2015 
Benchmark stars for Gaia Fundamental properties of the Population II star HD 140283 from interferometric, spectroscopic, and photometric data^{⋆,}^{⋆⋆,}^{⋆⋆⋆}
^{1}
Institut d’Astrophysique Spatiale, Université Paris XI, UMR 8617,
CNRS,
Bâtiment 121,
91405
Orsay Cedex,
France
email:
orlagh.creevey@ias.upsud.fr
^{2}
Laboratoire Lagrange, Université de Nice SophiaAntipolis, UMR
7293, CNRS, Observatoire de la Côte d’Azur, Nice, France
^{3}
Institutionen för fysik och astronomi, Uppsala universitet,
Box 516,
751 20
Uppsala,
Sweden
^{4}
MaxPlanckInstitute for Astronomy (MPIA),
Königstuhl 17, 69117
Heidelberg,
Germany
^{5}
Lowell Observatory, 1400 W. Mars Hill Road, Flagstaff, AZ
86001,
USA
^{6}
Department of Astronomy, Yale University,
New Haven, CT
06511,
USA
^{7}
LESIA, Observatoire de Paris, CNRS UMR 8109, UPMC, Université
Paris Diderot, PSL, 5 place Jules
Janssen, 92195
Meudon,
France
^{8}
UMIFCA, CNRS/INSU, France (UMI 3386), and Dept. de Astronomía,
Universidad de Chile, Santiago, Chile
^{9}
Université Grenoble Alpes, IPAG, 38000
Grenoble,
France
^{10}
CNRS, IPAG, 38000
Grenoble,
France
^{11}
Georgia State University, PO Box 3969, Atlanta
GA
303023969,
USA
^{12}
CHARA Array, Mount Wilson Observatory,
91023
Mount Wilson
CA,
USA
Received: 30 May 2014
Accepted: 20 November 2014
Metalpoor halo stars are important astrophysical laboratories that allow us to unravel details about many aspects of astrophysics, including the chemical conditions at the formation of our Galaxy, understanding the processes of diffusion in stellar interiors, and determining precise effective temperatures and calibration of coloureffective temperature relations. To address any of these issues the fundamental properties of the stars must first be determined. HD 140283 is the closest and brightest metalpoor Population II halo star (distance = 58 pc and V = 7.21), an ideal target that allows us to approach these questions, and one of a list of 34 benchmark stars defined for Gaia astrophysical parameter calibration. In the framework of characterizing these benchmark stars, we determined the fundamental properties of HD 140283 (radius, mass, age, and effective temperature) by obtaining new interferometric and spectroscopic measurements and combining them with photometry from the literature. The interferometric measurements were obtained using the visible interferometer VEGA on the CHARA array and we determined a 1D limbdarkened angular diameter of θ_{1D} = 0.353 ± 0.013 milliarcsec. Using photometry from the literature we derived the bolometric flux in two ways: a zero reddening solution (A_{V} = 0.0 mag) of F_{bol} of 3.890 ± 0.066 × 10^{8} erg s^{1} cm^{2}, and a maximum of A_{V} = 0.1 mag solution of 4.220 ± 0.067 × 10^{8} erg s^{1} cm^{2}. The interferometric T_{eff} is thus between 5534 ± 103 K and 5647 ± 105 K and its radius is R = 2.21 ± 0.08R_{⊙}. Spectroscopic measurements of HD 140283 were obtained using HARPS, NARVAL, and UVES and a 1D LTE analysis of Hα line wings yielded T_{eff}_{spec} = 5626 ± 75 K. Using finetuned stellar models including diffusion of elements we then determined the mass M and age t of HD 140283. Once the metallicity has been fixed, the age of the star depends on M, initial helium abundance Y_{i}, andmixinglength parameter α, only two of which are independent. We derive simple equations to estimate one from the other two. We need to adjust α to much lower values than the solar one (~2) in order to fit the observations, and if A_{V} = 0.0 mag then 0.5 ≤ α ≤ 1. We give an equation to estimate t from M, Y_{i} (α), and A_{V}. Establishing a reference α = 1.00 and adopting Y_{i} = 0.245 we derive a mass and age of HD 140283: M = 0.780 ± 0.010M_{⊙} and t = 13.7 ± 0.7 Gyr (A_{V} = 0.0 mag), or M = 0.805 ± 0.010M_{⊙} and t = 12.2 ± 0.6 Gyr (A_{V} = 0.1 mag). Our stellar models yield an initial (interior) metalhydrogen mass fraction of [ Z/X ] _{i} = −1.70 and log g = 3.65 ± 0.03. Theoretical advances allowing us to impose the mixinglength parameter would greatly improve the redundancy between M,Y_{i}, and age, while from an observational point of view, accurate determinations of extinction along with asteroseismic observations would provide critical information allowing us to overcome the current limitations in our results.
Key words: stars: fundamental parameters / stars: Population II / techniques: interferometric / techniques: spectroscopic / surveys / stars: individual: HD 140283
Based on NARVAL and HARPS data obtained within the Gaia DPAC (Data Processing and Analysis Consortium) and coordinated by the GBOG (GroundBased Observations for Gaia) working group, and on data retrieved from the ESOADP database.
Full Table 12 is only available at the CDS via anonymous ftp to cdsarc.ustrasbg.fr (130.79.128.5) or via http://cdsarc.ustrasbg.fr/vizbin/qcat?J/A+A/575/A26
© ESO, 2015
1. Introduction
The determination of accurate and precise stellar properties (radius R, effective temperature T_{eff}, age, log g, ...) of metalpoor halo stars is a primary requirement for addressing many astrophysical questions from stellar to galactic physics. For example, improving our knowledge of stellar interior processes (diffusion and Li depletion e.g. Bonifacio & Molaro 1998; Lebreton 2000; Charbonnel & Primas 2005; Korn et al. 2007; Meléndez et al. 2010) and evolution models has important consequences for determining the precise ages of stars and clusters (e.g. Grundahl et al. 2000), and thus for determining the age, and formation history of our Galaxy, e.g. Yamada et al. (2013), Haywood et al. (2013).
The effective temperature T_{eff} is a particularly important fundamental quantity to determine. This has consequences for deriving masses and ages of stars through a HertzsprungRussell (HR) diagram analysis, and determinations of absolute abundances require accurate T_{eff} preferably a priori along with accurate knowledge of surface gravities log g. However, determining T_{eff} is difficult, and calibrating this scale has been the subject of many recent important studies, e.g. Casagrande et al. (2010). While much progress has been made, the temperature scales still remain to be fully validated at the lowmetallicity and highmetallicity regimes. This is particularly important in the context of largescale photometric surveys where many low metallicity stars are observed and provide important tests of stellar structure conditions quite different to those of our Sun, and hence of Galactic structure and evolution.
With recent developments from asteroseismology, in particular for cool FGK stars using the spacebased telescopes CoRoT and Kepler, determinations of precise fundamental properties for large samples of stars is possible, with direct consequences for galactic astrophysics. Mean densities can be determined with precisions of the order of ~1% and in combination with other data, masses and ages reach precisions of 5% and <10% (Metcalfe et al. 2009; Silva Aguirre et al. 2013; Lebreton & Goupil 2014; Metcalfe et al. 2014). Asteroseismic scaling relations which predict masses, radii and surface gravity from simple relations using global seismic quantities (Brown & Gilliland 1994; Kjeldsen & Bedding 1995) have a huge potential for stellar population studies and galactoarchaeology. However, while these relations have been validated in some regimes such as main sequence stars (Metcalfe et al. 2014; Huber et al. 2012; Creevey et al. 2013), they have yet to be validated in the metalpoor regime (e.g. Epstein et al. 2014). Independent determinations of radii and masses of these stars provide valuable tests of the widelyapplicable relations.
Combining high precision distances that Gaia^{1} will yield with robust determinations of predicted angular diameters from temperaturecolour relations (e.g. Kervella et al. 2004; Boyajian et al. 2014) yields one of the most constraining observables for stellar models: the radius. These various arguments clearly demonstrate the need for very thorough studies of the most fundamental properties of nearby bright stars.
Observation log for HD 140283 using the VEGA interferometer on the CHARA array in threetelescope mode (telescopes E1, W1, and W2).
The Gaia mission (e.g. Perryman 2005) was successfully launched at the end of 2013. In addition to distances and kinematic information, it will also deliver stellar properties for one billion objects (BailerJones et al. 2013). In preparation for this mission a set of priority bright stars has been defined^{2}. Extensive observations and analyses have been made on these benchmark stars, with the aim of using them to define (and refine) the stellar models that will be used for characterizing the one billion objects Gaia will observe. One of these benchmark stars is HD 140283 and a thorough analysis of its fundamental properties is timely.
HD 140283 (BD10 4149, HIP 76976, α = 15h43m03s, δ = −10°56′01′′, [ l,b ] = [ 356.31°,33.61° ]) is a metalpoor Galactic Halo (Population II) star. It is bright (V = 7.2 mag) and nearby (58 pc) and thus a benchmark for stellar and galactic astrophysics. HD 140283 has been the subject of numerous studies, especially in more recent years where more sophisticated stellar atmosphere models have been used to determine accurate metallicities and α abundances (e.g. Charbonnel & Primas 2005) and to study neutroncapture elements to understand better heavy element nucleosynthesis (Siqueira Mello et al. 2012; Gallagher et al. 2010; Collet et al. 2009). Concerning the fundamental properties of this star, Bond et al. (2013) published a refined parallax π = 17.15 + −0.14 mas, an uncertainty one fifth of that determined by the Hipparcos satellite (van Leeuwen 2007, 17.16 ± 0.68 mas), and assuming zero reddening determined an age of 14.5 Gyr. More recently, VandenBerg et al. (2014) determined an age for HD 140283 of 14.27 Gyr. These ages are slightly larger than the adopted age of the Universe (13.77 Gyr, Bennett et al. 2013) but within their quoted error bars. The new better precision in the parallax from Bond et al. (2013) reduces the derived uncertainty in the radius (this work) by a factor of 30% and leaves the radius uncertainty dominated almost entirely by its angular diameter.
The objective of this work is to determine the radius, mass, age, effective temperature, luminosity, and surface gravity of HD 140283. We present the very first interferometric measurements of this object obtained with the visible interferometer VEGA (Mourard et al. 2011) mounted on the CHARA array (ten Brummelaar et al. 2005) in California, USA (Sect. 2). Our observations also show the capabilities of this instrument to operate at a magnitude of V = 7.2 (without adaptive optics) and at very high angular resolution (~0.35 milliarcsec, mas), one of the smallest angular diameters to be measured. We combine these interferometric data with multiband photometry (Sect. 3), to derive the fundamental properties of HD 140283 (Sect. 4). We then analyse highresolution spectra (Sect. 5) to determine its T_{eff} while imposing log g from our analysis. In Sect. 7 we use the stellar evolution code CESAM to interpret our observations along with literature data (Sect. 6) to infer the stellar model properties (mass, age, initial abundance). We discuss the limitations in our observations, models, and analysis (Sect. 8) and we conclude by listing the next important observational and theoretical steps for overcoming these limitations.
2. CHARA/VEGA interferometric observations and the angular diameter of HD 140283
Interferometric observations of HD 140283 were performed on four nights during 2012 and 2014 using the VEGA instrument on the CHARA Array and the instrument CLIMB (Sturmann et al. 2010) for 3T group delay tracking. The telescopes E1, W1, and W2 were used for all of the observations, providing baselines of roughly 100 m, 221 m, and 313 m. Table 1 summarizes these observations. In this table, the sequence refers to the way the observations were made where a typical calibrated point consists of observations of a calibrator star, the target star, and again a calibrator star. We obtained a total of five calibrated points. The extracted V^{2} (squared visibility) measurements for the target star alone are the instrumental V^{2} and these need to be calibrated by stars with known diameters. We used a total of three different calibrator stars whose predicted uniform disk angular diameters in the R band are given in the caption of the table. These were estimated using surface brightness relations for V and V − K as provided by the SearchCal tool of JMMC (Bonneau et al. 2006). SearchCal also provides the limbdarkenedtouniformdisk converted angular diameters, and the latter are used for calibrating the data. The average seeing during the full observation sequence is given in the following column by the Fried parameter r0. This parameter gives the typical length scale over which the atmosphere can be considered uniform. Perfect conditions can have r0 of the order of 20 cm, while poor conditions have r0 of approximately 5 cm. Our observations were conducted in poor to average seeing conditions (low altitude and nonoptimal observing season), resulting in larger errors on the data and the loss of some points (see below).
2.1. Extraction of squared visibilities
We used the standard V2 reduction procedure of the VEGA instrument as described in Mourard et al. (2011) to process the data. Depending on the quality of the data, the observations were processed in a number of bands of different width. The wavelength coverage of the VEGA R/I band is ~45 nm. For the observations from 2012 we processed the data in four spectral bands of 10 nm each, centred on 705 nm, 715 m, 725 nm, and 735 nm. Because of the lower signaltonoise ratio (S/N), for the 2014 observations we processed the data in one band of 20 nm centred on 710 nm^{3}. Table 1 summarizes this information under the column heading Band processed/width. The number of potential V^{2} points depends on the number of bands processed, the number of calibrated points, and the number of telescopes pairs (in this case three) for each observation sequence. In some cases (variable seeing, poor S/N on target or calibrator) the visibility calibration fails. Thus, in the final column of Table 1 we give the number of V^{2} points that were successfully extracted from each observation sequence, resulting in a total of 19 usable V^{2} points. Some of the data were processed independently and the resulting V^{2} varied only in the third decimal place with no consequence for the angular diameter derivation. The calibrated visibilities have an associated intrinsic statistical error σ_{STAT} and an external error σ_{EXT} coming from the uncertainty on the angular diameter of the calibrators, and the total error is given as the quadratic sum of the two components (Mourard et al. 2012^{4}). By observing on different nights and by using different calibrators we reduced the possibility of systematic errors in our analysis. Table 2 lists the squared visibilities V^{2}, along with the errors separated into the statistical σ_{STAT} and external errors σ_{EXT}, the projected baselines^{5} Bp, and the effective wavelengths λ of the 19 V^{2} points.
2.2. From visibilities to an angular diameter
To determine the uniform disk and 1D angular diameters, we performed nonlinear leastsquares fits of the squaredvisibility data to visibility functions using the LevenbergMarquardt algorithm (Press et al. 1992). Table 3 lists the uniform disk θ_{UD}, 1D limbdarkened θ_{1D}, and 3D limbdarkened θ_{3D} angular diameters and uncertainties, all explained in the following paragraphs.
2.2.1. Uniform disk angular diameter
Fitting the data to a uniform disk angular diameter ignoring the external errors yields θ_{UD} = 0.338 ± 0.011 mas with a . Including the external errors as gives θ_{UD} = 0.340 ± 0.012 mas with a .
2.2.2. 1D limbdarkened angular diameter
The 1D limbdarkened disk function is given by (1)where μ is the wavelength dependent limbdarkening coefficient, J_{n} is the Bessel function of order n, and x = πBθ_{1D}λ^{1}. The value of μ was obtained by interpolating the linear limbdarkening coefficients from Claret et al. (2012) (“u” in their work) for the stellar parameters T_{eff} = 5500 and 5750 K, log g = 3.75 and [M/H] = –2.10 (see Sect. 7). We adopted the coefficients halfway between the R and I tabulated values, which correspond to the effective wavelengths of the measurements. The resulting values vary between 0.4795 and 0.4870 for the range of T_{eff} found in this work. Using the values of log g = 3.50 or [M/H] = –2.0 yields limbdarkening coefficients that vary in the third decimal place and these result in a change in θ_{1D} in the fourth decimal place.
Angular diameter determination for HD 140283.
The value of μ can in principle also be fitted, however, the brightness contrast is larger at the border of the projected disk of the star that corresponds to visibilities close to and just after the first zero point. Our data do not extend to the first zero point nor is the coverage at lower visibility sufficient to constrain μ. Just as in the case of many other studies we adopt modeldependent values.
A fit to a 1D limbdarkened disk angular diameter adopting the errors σ = σ_{STAT} yields θ_{1D} = 0.353 ± 0.012 mas with . The data and the corresponding visibility curve are shown in the top panel of Fig. 1, while the residuals scaled by the individual errors on the measurements are shown in the bottom panel.
At such small angular diameters with calibrators of similar size the external error could have an important effect on our derived diameter. To show explicitly the effect of this error we performed Monte Carlotype simulations. We simulated 10 000 sets of visibility points and using the original statistical errors as the uncertainties in the data we fitted these data to 1D limbdarkened visibility functions. A simulated set of visibility points was obtained by adding a random number r drawn from a Gaussian distribution with a mean of zero and a standard deviation of 1, scaled by the external errors as follows: for each point i in each set. For the 10 000 simulated sets we obtained 10 000 fitted angular diameters θ_{1D}.
Fig. 1 Top: squared visibility measurements and the fitted visibility curve corresponding to a 1D diameter θ_{1D} = 0.353 mas. Bottom: data residuals scaled by the errors on the measurements. The symbols represent the night the data were taken: ⋄ = 2012 April 18, °= 2012 May 21, □= 2014 May 03, and ⋆ = 2014 May 04 (Sect. 2.2.2). 

Open with DEXTER 
Fig. 2 A Gaussian curve representing the fitted angular diameter and the 68% probability region (± σ_{θ,STAT}) due to the statistical errors (light grey shaded region). The vertical line illustrates the centre of the Gaussian corresponding to the fitted diameter θ. The inner distribution shows fitted angular diameters from 10 000 simulations by randomly adding the external error to the measurements. The dark grey shaded region illustrates the 68% probability region due to the external errors (± σ_{θ,EXT}, Sect. 2.2.2). 

Open with DEXTER 
The distribution of angular diameters from the simulations is shown in Fig. 2 by the inner grey distribution, where the binsize of 0.0007 mas is determined by the FreedmanDiaconis rule^{6}. The darkshaded region corresponds to the 68% probability region (± σ_{θ,EXT}, where the subscript “θ” refers to the uncertainty on the derived diameter as opposed to the measurement errors), and represents a total of 0.008 mas or σ_{θ,EXT} = 0.004 (see Table 3). We also overplot a Gaussian function representing the fitted angular diameter θ_{1D} and the 68% probability region corresponding to the statistical uncertainty ± σ_{θ,STAT} (light grey shaded region). In this case the external errors do not contribute significantly to the total error, but if the statistical errors were much smaller then it would be of interest to reduce this external error. The uncertainty on the final adopted 1D limbdarkened diameter is obtained by adding the statistical and external error in quadrature and this yields θ_{1D} = 0.353 ± 0.013 mas (see Table 6), in agreement with the uncertainty obtained when fitting the data using the total measurement errors (). Any systematic affecting the determination of the diameter can only be evaluated by comparing this result with an independently determined angular diameter.
2.2.3. 3D limbdarkened angular diameter
We also calculated the limbdarkened diameter considering realistic 3D simulations of surface convection. Generally θ_{UD}<θ_{3D}<θ_{1D} (Bigot et al. 2006, 2011; Chiavassa et al. 2009, 2012), which is due to a smoother hydrodynamical temperature gradient in the surface layers compared with the hydrostatic gradient, although for nonevolved stars this difference may not be significant even in the visible wavelengths which are more sensitive to this effect. Casagrande et al. (2014) recently pointed out a very small discrepancy between their absolute T_{eff} scale from the infrared flux method (IRFM) and measured angular diameters and attributed it to the fact that 1D interferometric diameters are usually quoted and not 3D ones. In most cases the precision in the measurements is not sufficient to distinguish between the two; however, their work did show an almost perfect agreement between the T_{eff} derived from a 3D angular diameter of the metalpoor giant HD 122563 from Creevey et al. (2012) and their IRFM value (a difference of 8 K in T_{eff}). For this reason we performed a 3D analysis. The method follows that of Bigot et al. (2011) and we refer the readers to this paper for details. Briefly, the 3D limbdarkened profiles are obtained by computing the full 3D monochromatic line transfer within the wavelength range of the VEGA red camera at different inclination angles. We then averaged the monochromatic limbdarkened profiles within each bandwidth (10−20 nm) around the central wavelengths. The parameters that define the 3D model of granulation were 5650 K (temporal average), log g = 3.65 and [M/H] = −2.0. By comparing the 3D limbdarkened intensity profile directly with the 1D profile using a limbdarkening law and the adopted 1D coefficients, we find an average difference of <1% between the two profiles for over 99% of the linear diameter of the star. A fit to the data yields θ_{3D} = 0.352 ± 0.013 mas very close to the 1D profile. The difference between the 3D and 1D is negligible for HD 140283 since its limb darkening is very weak because of the lack of metals in its atmosphere. For the rest of the work, we adopt the 1D value to derive the stellar parameters (see Table 6).
Fig. 3 Photometric measurements converted to absolute flux for HD 140283, with an example of a fit to a model spectrum (Sect. 3). 

Open with DEXTER 
3. Bolometric flux
Knowledge of the bolometric flux F_{bol} is required to calculate the T_{eff} using the angular diameter or the luminosity L using the parallax. It can be estimated by using bolometric corrections with the observed V magnitude of the star, or by using empirical formulae. A more reliable way to determine it is by fitting the spectral energy distribution (SED), as long as there is sufficient wavelength cover. We adopt the SED approach to determine F_{bol} and in this section we describe this analysis and compare it to literature values.
One of the biggest difficulties in determining the bolometric flux is the unknown reddening to the star. For stars that are close enough we can often make the assumption that interstellar extinction is nonexistent, or A_{V} = 0.0 mag. However, we also know that the distribution of gas and dust is not homogenous in the Galaxy and even stars that are close may suffer some degree of reddening. In fact, Bond et al. (2013) discussed this point for HD 140283 and remarked that by using HD 140283 as a standard candle and assuming a small degree of reddening yielded a distance to the globular cluster M 92 in much better agreement than without reddening.
HD 140283 is not very distant, and while some authors do indeed confirm zero reddening (Casagrande et al. 2010, 2011; Ramírez et al. 2010), several authors have also reported nonzero reddening to this star. Fuhrmann (1998) cited several sources of extinction measurements: E_{b − y} = 0.056, 0.020, 0.043 mag from McMillan et al. (1976); Schuster & Nissen (1989); Hauck & Mermilliod (1990), and E_{B − V} = 0.02 mag from Ryan (1989). Using E_{B − V} ~ 1.35 × E_{b − y} and A_{V} = 3.1 × E_{B − V} we obtain a range of A_{V} between 0 and 0.24 mag. Additionally, Arenou et al. (1992) estimates A_{V} = 0.13 mag for HD 140283’s galactic coordinates and distance, while Lallement et al. (2014) obtain A_{V} = 0.0087 mag (Puspitarini, priv. comm.). The adoption of one particular result for A_{V} then becomes arbitrary. In this work we determine A_{V} in the direction of the star along with F_{bol} using a SED fitting method, described below. However, we also fix A_{V} = 0.0 mag and discuss our results considering both possible scenarios.
3.1. Spectral energy distribution fitting methods
We used two different SED approaches to determine F_{bol} for HD 140283, by using a compilation of 324 literature photometry measurements of our target. The data were converted to flux measurements using the proper zeropoints, e.g. Fukugita et al. (1995), and the flux measurements are shown in Fig. 3 along with a model fit^{7}. The first fitting approach (Approach 1) is the method employed in von Braun et al. (2012) initially described in van Belle et al. (2008), and the second approach (Approach 2) is an independent fitting method described in Appendix A. The two methods were developed independently of each other by different authors.
Both approaches are based on fitting literature photometric data to libraries of stellar spectra. Approach 1 was initially developed to use the Pickles library of empirical spectra (Pickles 1998b), and here this is referred to as method 1A. In this work we also replaced the Pickles library in Approach 1 (method 1A) by two other libraries: (1B) the BASEL library of semiempirical spectra (Lejeune et al. 1997) and (1C) the PHOENIX library of synthetic spectra (see Hauschildt et al. 1999 and references therein), made available by Husser et al. (2013). Approach 2 implemented only the BASEL (2B) and the PHOENIX (2C) libraries. The main difference between the two approaches is that the former evaluates values based on fixed points along the spectral library, while the latter performs interpolation among the stellar parameters. Below we describe the specifics of each method.

1A
Method 1A performs a χ^{2} minimization fit of the literature photometry measurements of our target to the wellknown empirical spectral templates published in Pickles (1998b). The Pickles (1998a) spectral templates span 0.2−2.5 μm and provide blackbody interpolation across wavelength ranges without data. Beyond 2.5 μm, extrapolation is done along a blackbody curve of the input temperature of the spectral template. Numerical integration of the scaled template yields the bolometric flux of the star. The Pickles library has the advantage that it is purely empirical and so the spectra resemble the SED of true stars. However, the most metalpoor spectral template in this library has [Fe/H] = –0.6. The lack of an adequate (metalpoor) spectral template is indicated by the relatively high values found. It is this high value that motivated the implementation of the semiempirical and pure synthetic spectral libraries. Table 4 lists the best F_{bol}, A_{V}, and values corresponding to fits to the data imposing A_{V} = 0.0 mag and then also allowing A_{V} to be fitted.
Table 4Bolometric flux and interstellar reddening from fitting literature photometryconvertedtoflux measurements to the Pickles (A), BASEL (B), and PHOENIX (C) libraries of stellar spectra using two independent approaches (1 and 2).

1B
The BASEL library spans a longer wavelength region of 9–160 000 nm, hence extrapolation is not necessary. We calculated values based on a set of spectral templates with the following parameters: T_{eff} ranging from 5250–6000 K in steps of 250 K; log g = 3.0, 3.5, 4.0; and [Fe/H] = −2.0, −2.5, and −3.0. The best results are shown in Table 4 where we also indicate the T_{eff} of the best template spectrum. In both reddening scenarios, these correspond to log g = 3.5 and [Fe/H] = −2.5. We note that the nonzero reddening result corresponds to a spectral template at the edge of our parameter space and thus it should be treated with caution.

1C
Method 1C implemented the mediumresolution PHOENIX spectra, which cover a wavelength region of 300−2500 nm. Extrapolation beyond 2500 nm provides the missing flux, just as for method 1A. Using the highresolution spectra (see method 2C) we estimated the flux contribution between 50 and 300 nm as this accounts for ~4% of the total flux, and thus cannot be neglected. We inspected synthetic spectra ranging from 5400−6300 K in steps of 100 K, using log g = 3.50, [Fe/H] = –2.0 and –3.0 and α = + 0.4 (αenhanced elements).

2B
Method 2B implements the full BASEL library of semiempirical spectra and so no restrictions on the parameter ranges are imposed. Because of degeneracies between parameters, log g and [Fe/H] were fixed at 3.65 and –2.5, respectively. Varying log g and [Fe/H] by 0.5 each results in insignificant changes in the results (~0.001 of units given, equivalent to <0.05σ). To determine A_{V} we fixed it at a range of input values while fitting only T_{eff} and the scaling factor θ_{s} (see Appendix A) and chose the value that returned the best . Fixing T_{eff} at different values while fitting A_{V} yielded the same results. We tested that all of our results were insensitive to the initial parameters.

2C
This final method implemented the highresolution αenhanced (+0.4) PHOENIX library of spectra which span a wavelength range of 50−5500 nm. The parameter space was restricted to a T_{eff} between 5000 K and 6500 K, [Fe/H] = –3.0 and –2.0, and log g = 3.5−4.0, and in the final iteration log g and [Fe/H] were fixed to 3.65 and –2.5, respectively.
3.2. Bolometric flux of HD 140283
The results from each fitting method are given in Table 4. For the zero reddening results, there is excellent agreement between approaches 1 and 2, using both B and C, with differences between them of under 0.007 × 10^{8} erg s^{1} cm^{2} (<1σ). The disagreement between method 1A and the other methods is expected because of the known low metallicity of our target and the high confirms this. From these results the selection of one particular result between approaches 1 or 2 becomes arbitrary. The difference between the results B and C are of the order of 0.060 × 10^{8} erg s^{1} cm^{2} for both approaches 1 and 2, indicating simply that a systematic error arises from the choice of libraries of stellar spectra. This 0.06 difference is also consistent with the scatter among results from the literature for a zero reddening solution (see below). Because we choose just one solution, we add the 0.06 difference in quadrature to the derived uncertainty.
Considering F_{bol} with nonzero reddening, we find that methods 1A, 1B, and 2B are in agreement to within the errors. While all of the F_{bol} do differ there is consistency among them. Method C results are systematically lower than method B, just as for the zero reddening results, and the difference between method 1C and 2C stem directly from the fitted A_{V}.
There is no clear evidence to have a preferred value of A_{V}. For this work we therefore adopt A_{V} = 0.0 mag along with one of the largest nonzero values. We may then consider that our results fall somewhere between the two extremes. We choose to adopt the results from approach 2 because it interpolates in the template spectra to the best T_{eff} for any combination of log g and [Fe/H]. We then arbitrarily choose to use the results from method C (PHOENIX). The adopted F_{bol} is given in the lower part of Table 4 where the extra source of error arising from the choice of libraries has been accounted for, and these values are also reproduced in the reference Table 6 of observed stellar parameters derived in this paper.
Several authors have also provided estimates of F_{bol}, either from direct determinations or from applications of their empirically derived formula. In Table 5 we summarize these determinations, which assume zero reddening. Calculating a standard deviation of the results yields a scatter of 0.053 (10^{8} erg s^{1} cm^{2}). This value is within the representative 0.06 value that we add in quadrature to the uncertainty.
Determinations of F_{bol} by other authors assuming zero reddening.
4. T_{eff}, R, g, L, and bolometric corrections
Table 6 lists θ_{1D} (hereafter θ) and F_{bol} derived in this work. Combining θ, F_{bol}, and the StephanBoltzmann equation allows us to solve for T_{eff}, (2)where σ_{SB} is the StephanBoltzmann constant. The luminosity L is calculated from (3)where d is the distance and we adopt the value of the parallax from Bond et al. (2013). The radius is calculated from (4)or equivalently it is obtained from the StephanBoltzmann equation if in this equation T_{eff} is obtained from θ. The surface gravity g can also be estimated by imposing a value of mass M and using the standard equation g = GM/R^{2} where G is the gravitational constant. Here we adopt M = 0.80M_{⊙} with a conservative error of 0.10 M_{⊙}, where the uncertainty is large enough to safely assume no model dependence. Table 6 provides the reference list of these derived stellar properties.
The bolometric magnitude is calculated directly from the luminosity where we adopt a solar bolometric magnitude M_{bol, ⊙} = 4.74 mag. Considering the V magnitude (see Table 11), our derived A_{V}, along with the parallax we can then calculate the absolute V magnitude M_{V}. The resulting bolometric corrections in the V band are –0.179 mag and –0.168 mag for the zero and nonzero reddening solutions, respectively. These are all listed in Table 6.
For the T_{eff} values derived in this work we find tabulated bolometric corrections (BCs) from Flower (1996) using the corrected coefficients from Torres (2010) of –0.13 and –0.11 mag, respectively (zero and nonzero reddening values). These should be compared with –0.19 and –0.18 in order to be consistent with the adopted bolometric magnitude of the Sun (Flower assumes M_{bol, ⊙} = 4.73, see Torres 2010; we assume M_{bol, ⊙} = 4.74).
Converting 2MASS K_{s} (see Table 11) to JohnsonGlass K using K = K_{s} + 0.032^{8}, we derive an observed (V − K)_{0} = 1.590 and 1.491 mag. The bolometric corrections tabulated for K by Houdashelt et al. (2000) where we adopt log g = 3.50 and [Fe/H] = –2.50 along with our derived T_{eff} are BC_{K} = 1.486 and 1.424 mag. The bolometric correction in V is then given by BC_{V} = BC_{K} – (V − K)_{0} which yields –0.104 and –0.066 mag, where they adopt M_{bol, ⊙} = 4.71. For both Flower and Houdashelt, the BCs from the zero reddening case are in better agreement with our values, although the differences are significant.
Masana et al. (2006) derives a bolometric correction for HD 140283 of BC_{V} = − 0.202 ± 0.053 assuming M_{bol, ⊙} = 4.74. This value is consistent with ours within the errors.
Stellar properties of HD 140283 derived in this paper by adopting π = 17.15 ± 0.14 mas from Bond et al. (2013).
5. 1D LTE determination of T_{eff}
With the primary aim of testing the agreement between the interferometric and a spectroscopicallyderived T_{eff} for metalpoor stars, we performed a 1D LTE spectroscopic analysis of HD 140283 using Hα line profiles.
5.1. Observations and method
Spectroscopic observations of HD 140283 were obtained using three different instruments (see Table 7) as part of a spectral library produced by BlancoCuaresma et al. (2014). The HARPS spectrograph is mounted on the ESO 3.6 m Telescope (Mayor et al. 2003), and the spectra were reduced by the HARPS Data Reduction Software (version 3.1). The NARVAL spectrograph is located at the 2 m Telescope Bernard Lyot (Pic du Midi, Aurière 2003). The data from NARVAL were reduced with the LibreESpRIT pipeline (Donati et al. 1997). The UVES spectrograph is hosted by Unit Telescope 2 of ESO’s VLT (Dekker et al. 2000). Two UVES spectra were retrieved for HD 140283, one from the Advanced Data Products collection of the ESO Science Archive Facility (reduced by the standard UVES pipeline version 3.2, Ballester et al. 2000), and one from the UVES Paranal Observatory Project (UVES POP) library (Bagnulo et al. 2003, processed with data reduction tools specifically developed for that project). All spectra have been convolved to a spectral resolution of R = λ/ Δλ = 70 000, from a higher original resolution.
The spectrum calculations and the fitting procedure were done with the Spectroscopy Made Easy (SME) package^{9} (version 360, 2013 June 04, Valenti & Piskunov 1996; Valenti & Fischer 2005). This tool performs an automatic parameter optimization using a LevenbergMarquardt chisquare minimization algorithm. Synthetic spectra are computed by a builtin spectrum synthesis code for a set of global model parameters and spectral line data. A subset of the global parameters is varied to find the parameter set which gives the best agreement between observations and calculations. The required model atmospheres are interpolated in a grid of MARCS models (Gustafsson et al. 2008).
An important step in spectrum analysis with SME is the definition of a line mask, which specifies the pixels in the observed spectrum that should be used to calculate the chisquare. We defined three different line masks, covering different windows in the wings of Hα which are free of telluric and stellar lines (judged from the observed spectra). Line mask 1 is similar to that used by Cayrel et al. (2011) – four narrow windows on either side of the line centre, reaching rather far out into the wings. Line mask 2 is similar to that used for the SME analysis of the UVES spectra obtained in the GaiaESO survey (Bergemann et al., in prep.) – two somewhat broader windows on either side, closer to the line centre than for Line mask 1. Line mask 3 is similar to that used by Ruchti et al. (2013), covering as much of the line wings as possible between 1 and 10 Å from the line centre on both sides. The line masks are visualized in Fig. 4, together with the NARVAL spectrum.
In an independent step from the parameter optimization, SME normalizes the spectrum to the local continuum by fitting a straight line through selected points in the observed spectrum defined by a continuum mask. We defined two sets of continuum masks. The narrow continuum mask has two windows about 22 and 29 Å from the line centre on the blue and red side, respectively. The wide continuum mask has two windows about 32 and 34 Å from the line centre on the blue and red side, respectively. The continuum mask windows are between 0.5 and 0.9 Å wide.
Information on observed spectra used for the fit: instrument/archive name, date of observation, mean S/N.
Fig. 4 Top: observed NARVAL spectrum (black) of the Hα profile of HD 140283 and fitted spectrum (red) using the three narrow continuum masks (shaded regions) for log g = 3.65 and [Fe/H] = –2.46. Bottom: residuals between the observed and fitted spectrum for each line mask, with the black dots indicating the line masks. See Sect. 5 for details. 

Open with DEXTER 
Atomic data for spectral lines contributing to flux in Line mask 3 for HD 140283.
The atomic data were taken from the line list compiled for the GaiaESO public spectroscopic survey (Gilmore et al. 2012). In the spectral interval from 6531 Å to 6598 Å, the list contains 70 atomic lines in addition to Hα, but only three of them have a possible contribution to the flux in the Line mask 3 wavelength regions for HD 140283. The atomic data for these lines are given in Table 8. The gfvalue used for the Hα line is 0.710 (Baker 2008, see Wiese & Fuhr 2009). The broadening of the Hα line by collisions with neutral hydrogen was calculated using the theory of Barklem et al. (2000), which is an improvement on the theory of Ali & Griem (1966). The more recent calculations by Allard et al. (2008) extend the description of the selfbroadening for Hα even further, but the resulting line profile is very similar to that using the Barklem et al. (2000) theory. The differences are smaller than the uncertainties in the observations.
5.2. Results
For the present study we adopted the observed log g derived in this work (log g = 3.65 ± 0.06, see Sect. 4). For [Fe/H] we adopted the mean and standard deviation of the reported values from the literature between 1990 and 2011, as listed in the PASTEL catalogue (Soubiran et al. 2010), of [Fe/H] = −2.46 ± 0.14. Line broadening from microturbulence v_{micro}, macroturbulence, v_{macro}, and vsini (rotation) were all fixed, v_{micro} = 1.3 km s^{1}, v_{macro} = 1.5 km s^{1}, and vsini = 5.0 km s^{1}, and changing them had no impact on the final result.
The initial value for the free parameter T_{eff} was chosen to be 6000 K (using an initial value of 5000 K resulted in the same solution for T_{eff}, to within 10 K). The results using log g = 3.65 and [Fe/H] = –2.46, rounded to the nearest 10 K, are given in Table 9 for each mask and each observation set along with its corresponding . The mean and scatter among these values is T_{eff} = 5642 ± 63 K.
Bestfitting T_{eff} values rounded to 10 K and reduced chisquare () for the different observations and masks for HD 140283.
As can be seen from Table 9, the best and comparable values are obtained using the NARVAL and UVES archive data, although the differences between their T_{eff} are of the order of 100 K. The observations that show the least sensitivity to the different masks are also the UVES archive data with a scatter of only 28 K, while those with the most sensitiviy are the NARVAL data (50 K).
Mean T_{eff} values fixing log g and [Fe/H] at their mean and ± 1σ values.
To illustrate the effect of changing log g and [Fe/H], in Table 10 we show the mean T_{eff} of the 24 results (averaging over the four sets of observations and six masks each) for nine combinations of log g and [Fe/H] considering the adopted errors in these parameters. The final column shows the typical scatter s found for each value given in that row. The T_{eff} are almost identical for [Fe/H] = –2.46 and –2.32, while the results for [Fe/H] = –2.60 decrease by ~30−40 K. To get a more realistic estimate of the error due to using different values of log g and [Fe/H], we performed Monte Carlo simulations and determined 24 values of T_{eff} (four observations, six masks) by adopting a randomly chosen log g and [Fe/H], where these were drawn from a normal distribution shifted and scaled by the means and standard deviations as given above. The simulations were repeated 100 times. The mean and standard deviation of the 2400 fitted T_{eff} is 5626 ± 75 K^{10}. This is the value that we adopt, and it is listed in the reference Table 6.
The interferometric T_{eff} falls between a mean value of 5534 K and 5647 K ± 105 K depending on the amount of reddening between the top of our Earth’s atmosphere and the star. We can then conclude that the 1D LTE determination of T_{eff} from Hα line wing fitting yields results that are consistent with the interferometric values and are thus valid for stars at this lower metallicity regime within the typical uncertainties for spectroscopic observations.
Adopting our spectroscopic T_{eff}, along with F_{bol}, and π, we determine R using the StefanBoltzmann equation, R = 2.14 ± 0.04R_{⊙} and 2.23 ± 0.04R_{⊙}, without and with reddening, respectively, or the angular diameter of 0.341 mas and 0.356 mas, respectively, ignoring the distance. Again assuming M = 0.80 ± 0.10M_{⊙}, we obtain an evolutionmodelindependent determination of log g, 3.68 ± 0.06 and 3.64 ± 0.06, respectively, where the error is dominated by the imposed error on mass. These spectroscopically derived values are also given in Table 6.
6. Previously observed stellar parameters
Using the General Catalogue of Photometric Data (Hauck et al. 1995), we found several sources of V magnitudes for HD 140283 more recent than 1980. These are 7.194, 7.212, 7.21, 7.22, 7.22, and 7.20 mag (Griersmith 1980; Cousins 1984; Norris et al. 1985; Mermilliod 1986; Carney & Latham 1987; Upgren et al. 1992). Adopting the mean and standard deviation of these we obtain V = 7.21 ± 0.01 mag. Using all of the available photometric measurements that date back as far as 1955 we obtain a mean V = 7.214 ± 0.018 mag with a mean error of 0.051 mag (these were the data used in Sect. 3). The infrared magnitudes from the 2MASS catalogue (Skrutskie et al. 2006) are H = 5.696 ± 0.036 mag and K_{s} = 5.588 ± 0.017 mag. The V and K_{s} magnitudes are given in Table 11.
We can apply the surfacebrightnesscolour relations calibrated by Kervella et al. (2004) and Boyajian et al. (2014) to HD 140283, to check the validity of these relations in the low metallicity regime. Combining the B magnitude from Morel & Magnenat (1978)B = 7.71 with the V and the infrared magntiudes H and K_{s}, the Kervella et al. (2004) relations yield consistent 1D limbdarkened angular diameters between the different colour indices of θ_{1D} = 0.321 mas, while considering a reddening correction corresponding to A_{V} = 0.10 mag. Considering no reddening correction at all results in a marginal change of the predicted angular diameter (θ_{1D} = 0.325 mas). Using the more newly calibrated relations from Boyajian et al. (2014) we also obtain θ_{1D} = 0.325 mas, where such relations yield a scatter of the order of 5%, which in this case corresponds to an uncertainty of approximately 0.016 mas. Our measured value is θ_{1D} = 0.353 ± 0.013, a difference of just over 2σ. This difference could perhaps be explained by the surface brightness relations being calibrated with more metalrich stars.
Jofré et al. (2014) recently presented a homogenous analysis of 34 FGK benchmark stars to derive their metallicities. They use several methods to determine [Fe/H] including the SME method presented here. For HD 140283 they adopt T_{eff} = 5720 K and log g = 3.67 and then derive an NLTE corrected [Fe/H] by combining individual line abundances of neutral lines. They also study the sensitivity of [Fe/H] to T_{eff} and find that a difference of 120 K results in a [Fe/H] difference of 0.04 dex. Applying this correction to the spectroscopic T_{eff} found here (~–90 K) yields a NLTE [Fe/H] = −2.39 and an LTE [Fe/H] = −2.46, in line with the value used in Sect. 5, but with lower uncertainties of 0.07 dex. We note that in the initial analysis presented in their paper, individual methods based on a global fit to the spectrum yielded slightly lower LTE [Fe/H] between −2.44 and −2.57. These last results have also been found by Gallagher et al. (2010) and Siqueira Mello et al. (2012): [Fe/H] = −2.59 and −2.60 adopting T_{eff} = 5750 K and log g = 3.70.
Thévenin (1998) derive [Fe/H] and abundances for the α elements using a T_{eff} = 5600 K and log g = 3.2 dex; [Fe/H] = −2.50, [Mg/H] = −2.05, [Si/H] = −2.10, [Ca/H] = −2.10, and [Ti/H] = −2.35, giving an average/median value of [α/Fe] = +0.35/+0.40. Adopting the NLTE [Fe/H] of −2.39 (by adapting the Jofré et al. 2014 value) along with [α/Fe] from Thévenin (1998) yields a global metallicity [M/H] = −2.12 / − 2.08 adopting the coefficients provided by Salaris et al. (1993) or [M/H] = −2.14 / − 2.10 adopting those provided by Grevesse & Noels (1993) as given in Yi et al. (2001). To account for the uncertainty on the α elements along with the 0.14 dex on [Fe/H] as discussed in Sect. 5, we adopt a conservative 0.20 dex as an error bar on [M/H]. Finally, Bonifacio & Molaro (1998) derive a Li abundance of A(Li) = 2.146, and Charbonnel & Primas (2005) derive A(Li) = 2.26−2.28. These parameters are all summarized in Table 11.
7. Interpretation of new data using stellar evolutionary models
We used the CESAM stellar evolution and structure code (CESAM2k; Morel 1997; Morel & Lebreton 2008) to derive the modeldependent properties of mass M, age, t, initial helium abundance mass fraction Y_{i}, and initial chemical composition Z_{i}/X_{i}, where Z and X denote heavy metals and hydrogen with Z_{i} + Y_{i} + X_{i} = 1. The initial abundances are at zeroage main sequence (ZAMS).
7.1. CESAM2K physics
The input physics of the models consist of the equation of state by Eggleton et al. (1973) with Coulomb corrections, the OPAL opacities (Rogers & Iglesias 1992) supplemented with Alexander & Ferguson (1994) molecular opacities. The p−p chain and CNOcycle nuclear reactions were calculated using the NACRE rates (Angulo 1999). We adopted the solar abundances of Grevesse & Noels (1993) (Z/X_{⊙} = 0.0243) and used the MARCS model atmospheres (Gustafsson et al. 2003) for a metallicity of [M/H] = –2.0 for the boundary conditions and reconstruction of the atmosphere. Convection in the outer envelope is treated by using the mixinglength theory described by Eggleton (1972), where l = αH_{p} is the mixinglength that tends to 0 as the radiative or convective borders are reached, H_{p} is the pressure scale height, and α is an adjustable parameter. Calibrating this value for the solar parameters we obtain α ~ 2.0. Microscopic diffusion was taken into account and follows the treatment described by Burgers (1969), and extra mixing is included by employing a parameter, Re_{ν}, as prescribed by Morel & Thévenin (2002) which slows down the depletion of helium and heavy elements during evolution. Nordlander et al. (2012) showed that the abundance difference between turnoff stars and red giants in the globular cluster NGC 6397 ([Fe/H] = −2.0) is typically 0.2 dex. We chose a value of the extramixing parameter Re_{ν} = 5, to produce a similar abundance difference. Figure 5 illustrates the evolution from ZAMS of the surface metallicity with T_{eff} by employing Re_{ν} = 5. As can be seen the initial surface metallicity is quite different to the observed one, which is illustrated by the arrows. For the metallicity of the star we adopt a [M/H] (≡[ Z/X ] _{s}) = −2.10 ± 0.20 as argued in Sect. 6. Using Z/X_{⊙} = 0.0243, without considering diffusion of chemical elements in the models, we obtain Z_{i}/X_{i} = 0.00022 + 0.00013 − 0.00008. From Fig. 5 it can be seen that the observed metallicity is lower than the initial one, and so the first constraint we have on the initial metallicity is Z_{i}/X_{i}> 0.00022.
Fig. 5 The evolution of the observed (surface) metal abundance, shown as a function of T_{eff}, for a 0.80 M_{⊙} model with Y_{i} = 0.245 including the effects of diffusion and extra mixing for two initial metalhydrogen ratio mass fractions (Z_{i}/X_{i} = 0.0005, 0.0007). At T_{eff} = 5600 K, we highlight the values of the luminosity for both of the models. See Sect. 7.1 for details. 

Open with DEXTER 
Each stellar evolution model is defined by a set of input model parameters – mass M, initial helium content Y_{i}, initial metal to hydrogen ratio Z_{i}/X_{i}, age t, the mixinglength parameter α, and the extramixing coefficient Re_{ν} – and these result in model observables, such as a model T_{eff} and a model L. By varying the input parameters we aimed to find models that fitted the derived L, T_{eff}, and [M/H] as given in Tables 6 and 11.
7.2. Approximate parameters of HD 140283 using L
The luminosity L provides a very important constraint on the evolution stage of the star. Because it is a measure of the energy production in the star, we can ignore the mixinglength parameter α which is relevant only to the outer convective envelope and hence R and T_{eff}. In Fig. 6 we show the evolution of luminosity as a function of age, for models of different mass (denoted with a number) and changing Y_{i} and Z_{i}/X_{i}, with the reference model of 0.80 M_{⊙} having Y_{i} = 0.245 and Z_{i}/X_{i} = 0.0005. We also highlight the error bars associated with L when we adopt the reddened and unreddened bolometric flux (reddened is a higher value). We have limited the age range to a minimum expected value, i.e. the minimal age of the old galactic clusters ~10 Gyr.
Using only the luminosity, Fig. 6 illustrates the following: 1) the mass of the star is between 0.78 and 0.84 M_{⊙} if we adopt a minimum Y_{i} = 0.245, consistent with a primordial value predicted by Big Bang nucleosynthesis. 2) Increasing Z_{i}/X_{i} by 0.0002 (typical error, see next section) changes the age by about 0.3 Gyr. 3) The effect of increasing Y_{i} to an upper limit of 0.260 (=Y_{i} + 3ΔY, where ΔY = 0.005 is indicated in Fig. 6) leads to a corresponding decrease in the age of the reference model by ~1.5 Gyr if everything else is left fixed, and imposing an approximate age of the Universe would decrease the lower limit in the mass to 0.75 M_{⊙}. Figure 6 depicts the approximate correlations between M and Y_{i}, given just the luminosity constraint. In this figure, we have denoted an upper age limit of ~14 Gyr (approximately the age of the Universe) by the dotted horizontal bars for each Y_{i} which imposes a lower limit in mass. We then depict the upper limit in mass when we assume that Halo stars should have begun forming at approximately >12 Gyr (continuous lines), and when we relax this constraint to the age of the youngest globular clusters of >10 Gyr (dashed lines). Using this assumption along with the assumption that the initial helium abundance is low (Y_{i}< 0.26), we obtain a modeldependent mass of HD 140283 of between 0.75 and 0.84 M_{⊙}. This then implies log g = 3.65 ± 0.03 by imposing R from this work.
Fig. 6 Left: evolutionary tracks for models of different masses (labelled by their solar masses) for Y_{i} = 0.245 and Z_{i}/X_{i} = 0.0005 plotted as a function of age. The dashed lines illustrate changing the M = 0.80M_{⊙} model in Y_{i} (ΔY) and Z_{i}/X_{i} (ΔZ). The unreddened and reddened luminosity constraints are indicated by the error bars. Right: approximate correlations between M and the initial helium abundance Y_{i} assuming an age >12 Gyr (continuous) and then extending it to >10 Gyr (dashed), given the luminosity constraint. See Sect. 7.2 for details. 

Open with DEXTER 
7.3. Refining the parameters using the angular diameter constraint
7.3.1. Constraining the mixinglength parameter α
The radius and the T_{eff} of a star are extremely sensitive to the adopted mixinglength parameter α which parametrizes convection in the convective envelope. It is an adjustable parameter and it needs to be tuned to fit the observational constraints. However, because of the correlations among different parameters and sets of input physics to models, in many cases it is left as a fixed parameter, and its value is set to the value obtained by calibration of the solar parameters with the same physics. This is also the case for many isochrones that are available publicly. While adopting α as the solarcalibrated value may be valid for stars that are similar to the Sun, for metalpoor evolved stars (different T_{eff}, log g, and [Fe/H]) this is certainly not the case as has been shown by Casagrande et al. (2007), for example, who suggest a downward revision for metalpoor stars, or Kervella et al. (2008) for 61 Cyg, or more recently by Creevey et al. (2012) for two metalpoor stars, and from asteroseismology for a few tens of stars, e.g. Metcalfe et al. (2012), Mathur et al. (2012), Bonaca et al. (2012). Fixing α with a solar calibrated value for HD 140283 yields no models consistent with observations. However, as we know that the age of HD 140283 is high (>10 Gyr) and its mass is relatively welldetermined as argued in the above paragraph (0.75 <M< 0.84M_{⊙}), along with its metallicity, α is one of the free parameters that we can in fact constrain, as well as further constraining its age and mass.
In order to reproduce the R and T_{eff} of HD 140283, it was necessary to decrease the value of α to about 1.00 (solar value ~2.0). Figure 7 illustrates how the value of α is constrained by the stellar evolution models. We show the HR diagram depicting the position of the error box of the reddened solution. Assuming Y_{i} = 0.245 the mass is between 0.78 and 0.84 M_{⊙} (see Fig. 6, right panel). We plot as continuous lines a 0.78 (blue), 0.80 (black), and 0.82 (green) M_{⊙} model with a mixinglength parameter of α = 1.0. Increasing α to 1.30 results in the 0.78 M_{⊙} evolution track falling inside the observed constraints (blue dashed). Decreasing α to 0.70 for the 0.82 M_{⊙} also results in the track falling within the error box (green dashed). Directly from the models we can establish a plausible range of α for a given mass for Y_{i} = 0.245:

0.78 M_{⊙} implies 1.10 <α< 1.50,

0.80 M_{⊙} implies 0.85 <α< 1.25, and

0.82 M_{⊙} implies 0.5 <α< 0.9.
The red dashed line above the error box shows a 0.84 M_{⊙} model with a value of α = 0.5. We were unable to make exolutionary tracks with values of α< 0.5, thus establishing a modeldependent lower limit on its value, and restricting the mass of the model to M ≤ 0.83M_{⊙}.
We can summarize M as a function of α by using the following relation M = 0.885 − 0.080α (± 0.01) M_{⊙} or alternatively as α = 11.013 − 12.449M (± 0.2) where the errors are the 1σ uncorrelated errors. With a mass between 0.78 and 0.83 M_{⊙} and α between 0.5 and 1.50, we set our reference model to the mean mass M = 0.805 ± 0.01M_{⊙} which implies α = 1.00 ± 0.20. This is the solution adopted for Y_{i} = 0.245 (Table 12), and for comparison with other solutions we adopt α = 1.00 as the reference value.
Fig. 7 Correlations between M and α for the reddened solution (error box) imposing Y_{i} = 0.245, illustrating the plausible ranges for these parameters imposed by stellar evolution. Similar colours have the same mass and the continuous lines have α = 1.0. See Sect. 7.3.1 for details. 

Open with DEXTER 
7.3.2. Mass, initial helium, and age of the nonreddened star
Figure 8, left panel shows the HR diagram error box and stellar models that satisfy or are close to the unreddened solution. The most central model corresponds to a mass of M = 0.78M_{⊙}, an initial helium abundance of Y_{i} = 0.245, initial metallicity Z/X_{i} = 0.0005, and α = 1.00. We also highlight models of other masses (blue) with the same Y_{i}, Z/X_{i}, and α that pass through or close to the error box. Models considering different α and Y_{i} (green dashed lines) are also shown. We apply the approach from the previous section to establish a relationship between M and α, and we also take Y_{i} into account. Given two of M,α,Y_{i} the third can be derived as
where these equations must satisfy 0.77 ≤ M ≤ 0.82, 0.5 ≤ α ≤ 1.5 and 0.245 ≤ Y_{i} ≤ 0.260. Using the reference α = 1.0 we give two solutions in Table 12 corresponding to Y_{i} = 0.245 (model 1, M1), and Y_{i} = 0.255 (model 2, M2), resulting in masses of 0.78 and 0.77 M_{⊙}.
To determine the age and other properties corresponding to these models we consider the properties of all of the models that pass through the 1σ error box. These properties are given in the lower part of Table 12 along with their uncertainties. For M1 we obtain an age of 13.68 Gyr with [M/H] = –2.03.
Fig. 8 HR diagrams illustrating the error box corresponding to the observations derived in this work for HD 140283. The left and right panels show the results when we adopt A_{V} = 0.00 and 0.10 mag. The dotted box in the left panel shows the error box when the distance is increased by 1σ, also illustrated by the arrow. In the right panel (nonzero reddening solution) we also show the A_{V} = 0.00 error box (dotted box), located in the lower right part of the figure. The central models correspond to Y_{i} = 0.245,α = 1.00, and Z/X_{i} = 0.0005 (with Re_{ν} = 5). The optimal masses are M = 0.780 and M = 0.805M_{⊙}, and are shown as the blue continuous lines passing near the centre of the error boxes. We also illustrate other models by changing the mass (blue) and the other parameters (green). The red dotted line corresponds to the radius derived in this paper. See Sect. 7.3 for details. 

Open with DEXTER 
7.3.3. Mass, initial helium and age of the reddened star
Figure 8 right panel illustrates the optimal models that satisfy the L, T_{eff} and [M/H] constraints (central continuous error box) for the reddened solution. The dotted error box shows the position of the unreddened constraints. Following the approach above we establish equations for deriving one of M,α,Y_{i} as a function of the other two. where the equations must satisfy 0.78 ≤ M ≤ 0.83, 0.5 ≤ α ≤ 1.5 and 0.245 ≤ Y_{i} ≤ 0.260. In Sect. 7.3.1 we found an optimal model of M = 0.805 ± 0.010M_{⊙} for Y_{i} = 0.245, and adopting the reference α = 1.00, we obtain M = 0.795M_{⊙} for Y_{i} = 0.255. The ages of these models are 12.17 and 11.95 Gyr respectively with metallities of [M/H] = –2.06 and –2.07. These properties are summarized in the lower part of Table 12 under the column headings M3 and M4.
Fig. 9 Relationship between the age and the mass of the models that satisfies the observational constraints. The circles/diamonds correspond to the Y_{i} = 0.245 / 0.255 solution, and the red denotes the A_{V} = 0.10 mag solution. The crosses denote the values corresponding to the centres of the error boxes in the HR diagram (Fig. 8) and the solutions given in Table 12. See Sect. 7.3.4 for details. 

Open with DEXTER 
7.3.4. The global solution
The determination of precise stellar properties that cannot be observed is extremely difficult, particularly because of the strong correlations between M,α, and Y_{i}, and consequently the age. The luminosity constraint along with the stellar models and the assumptions on the allowed age range of the star refined the mass to 0.75 <M< 0.84 and consequently log g to 3.65 ± 0.03. The angular diameter then further constrained α to within a range of 0.5−1.5 and M to 0.77 <M< 0.83, but most importantly it provided a tight correlation between them, along with Y_{i}. Our models have a surface metallicity of [M/H] ~− 2.05, in good agreement with the observed one (–2.10).
Determining the age of the star depends on the adopted values of of M,α,Y_{i}, where only two are independent and the third is calculated from Eqs. (5)–(10). The full range of mass values (by varying Y_{i} and α) has the largest impact on the determination of age, followed by Y_{i}. The value of α only indirectly influences the age through the required change of M and Y_{i}. To attempt to quantify this effect, in Fig. 9 we show the ages of the models that fit the observational constraints for a range of parameters. We varied the mass (xaxis) and Y_{i} (circles: Y_{i} = 0.245, diamonds: Y_{i} = 0.255) within the allowed correlations as shown in Fig. 6. For each combination we adjusted α to within its allowed range until the evolution track passed through the centre of the error box. This was done for the reddened solution (shown in red) and the unreddened solution (black). The figure shows the age of the models for these M,Y_{i} combinations. The crosses show the solutions in Table 12 (the central models in Fig. 8).
Allowing for all plausible combinations of M,α,Y_{i} results in a possible age of HD 140283 varying between 10.5 Gyr and 14.0 Gyr, with a lower limit of just under 11 Gyr for the zero reddened solution. The required values of α are systematically lower for the zero reddened solution with a typical value of 0.5 ≤ α ≤ 1.00 while the reddened solution requires 0.7 ≤ α ≤ 1.3.
We can parametrize the age t as a function of M,Y_{i}, and A_{V} using the following equation (11)where 0 ≤ A_{V} ≤ 0.10, 0.245 ≤ Y_{i} ≤ 0.260, 0.76 ≤ M ≤ 0.83, and the combination of M and Y_{i} necessarily implies α. The value of α can be approximated by (12)and is restricted to values between 0.5 and 1.5 (lower values for higher masses and vice versa). The full range of ages and α for all combinations of M,Y_{i},A_{V} are given in the online version.
Here it is clear that in order to determine the age of the star with the best precision we need to restrict the values of one of the parameters M,α,Y_{i} (see Sect. 8). Resolving the reddening problem would also relieve this degeneracy.
Assuming that the mixinglength parameter does not vary much within the atmospheric parameters of our models (log g = 3.65, [M/H] ~− 2.05, T_{eff}~5530–5660 K), and by requiring that the different solutions adopt a similar value, we can judge the impact of reddening and Y_{i} on the mass and age determination. Table 12 shows four solutions that adopt α = 1.00 for two values of Y_{i} and the two extreme values of A_{V}. The effect of increasing Y_{i} by 0.01 is a decrease in the mass and the age by 0.01 M_{⊙} and ~0.2 Gyr. The biggest difference comes from the adopted reddening, giving a 1.5 Gyr difference between the solutions with A_{V} = 0.0 and 0.1 mag. In principle we may expect to be able to impose a value of α in the near future either from the calibration of this parameter with large samples of targets, e.g. from asteroseismology (Bonaca et al. 2012; Metcalfe et al. 2014), or from advances in 3D numerical techniques, e.g. Trampedach & Stein (2011). Imposing this parameter externally would reduce the total range of mass, Y_{i}, and possibly extinction, thus delivering a more reliable determination of its age.
The uncertainties that we give in Table 12 for M, Z_{i}/X_{i}, and α are the 1σ uncorrelated uncertainties. These are obtained by fixing the parameters and varying the one of interest until we reach the edge of the error box. The uncertainties in the other parameters are obtained by varying the abovementioned parameters up to 1σ at the same time and requiring that they remain within the error box, e.g. for solution M3 the mass took values of between (0.815−0.790)/2 M_{⊙} (rounded), while α varied between 1.2 and 0.8, and Z_{i}/X_{i} between 0.0003 and 0.0007.
8. Discussion
8.1. Observational results
In this work we determined three important quantities for HD 140283: the angular diameter θ, the bolometric flux F_{bol}, and a 1D LTE T_{eff} derived from Hα line wings. Combining the first two observational results with a wellmeasured parallax from Bond et al. (2013) yields the stellar parameters R, L, T_{eff}, and log g. In this work we found that the biggest contributor to our systematic error in L and T_{eff} is the existence or not of interstellar reddening. The adoption of zero or nonzero reddening leads to a change in L and T_{eff} by ~0.3 L_{⊙} and 100 K, respectively. The determination of the spectroscopic T_{eff} led to a value closer to the reddened solution (5626 K). If we used this spectroscopic value to deduce F_{bol} and thus reddening (using the angular diameter), we would obtain F_{bol} = 4.16 × 10^{8} erg s^{1} cm^{2}, giving a flux excess of ~0.27 × 10^{8} and corresponding to A_{V}~ 0.084 mag or E(B − V) = 0.027 mag by adopting an extinction law of R(V) = 3.1. This value would be in agreement with the value proposed by Bond et al. (2013).
Fig. 10 Spectroscopic determinations of T_{eff} and [Fe/H] from the literature (coloured symbols) and our spectroscopic determination using Hα line wings (black circle). The shaded area illustrates the interferometric determination of T_{eff} from the nonzero reddening solution, while the dashed lines illustrate the boundaries of the interferometric T_{eff} from the zero reddening solution. 

Open with DEXTER 
8.2. 1D LTE spectroscopic T_{eff}
In this work we determined a 1D LTE spectroscopic T_{eff} using Hα line wings while imposing log g = 3.65 ± 0.06 and [Fe/H] = −2.46 ± 0.14. The log g was imposed from this work using the interferometric R and an estimate of M = 0.80M_{⊙}, later confirmed with models to be correct. The question we posed was whether such an analysis is capable of determining T_{eff} of metalpoor stars where generally missing physics in atmospheric models can lead to a large variety of T_{eff} determinations with correlations between the parameters. In Fig. 10 we show determinations of the spectroscopic parameters from various authors as provided by the PASTEL catalogue (Soubiran et al. 2010). We show the [Fe/H] determinations against T_{eff} with a colour and symbol code corresponding to different log g. The blue filled circles are those with 3.6 < log g< 3.7 corresponding to the constraints we can impose from this work (in many cases their log g was constrained prior to the analysis). The black circle marks our determination of T_{eff}. The shaded background illustrates the ± 1σ determination of the interferometric T_{eff} assuming reddening, and the dashed lines illustrate the boundaries with an assumption of zero reddening. If we can show that we need to consider reddening then we can conclude that the spectroscopic analysis is indeed capable of reproducing the T_{eff} of such metalpoor stars, if we can impose a reliable constraint such as log g.
Recently Ruchti et al. (2013) studied systematic biases in T_{eff} determinations, and they found that the T_{eff} determined from Hα were hotter than those derived from the angular diameter for metalpoor stars by 40–50 K (including HD 140283). Our results are in agreement with this statement if we assume some mild absorption to the star, e.g. A_{V} ~ 0.05 mag.
For HD 140283, Ruchti et al. (2013) used surface brightness relations to estimate the interferometric T_{eff} for this star (5720 K), and they obtained an HαT_{eff} = 5775 K, 150 K hotter than ours. There are a few explanations for this difference: 1) the log g and [Fe/H] values adopted affect the derived T_{eff} (see Table 10); however, they did not specify these values in their work; 2) the observations give different results – our analysis with NARVAL data yields higher T_{eff} than the other three sets; 3) the theoretical profiles were calculated using different model atmospheres and code; and 4) selfbroadening was computed using a different (older) theory. We believe that our analysis using different line masks and four sets of observations should yield a more reliable result. The only way to unveil the origin of the systematic difference is to use each other’s observations and compare the results.
8.3. Initial metallicity and the extramixing parameter in stellar models
The mean surface metallicity of our most reliable models is [ Z/X ] _{s} = − 2.05, with an initial value of [ Z/X ] _{i} = − 1.70 or Z_{i}/X_{i} = 0.0005 (the surface metallicity is mostly set by the initial value with some slight variation due to α). This value was obtained by adopting a diffusion parameter in the stellar models, one that we determined by comparing the difference in surface metallicity abundance of stars at their turnoff stage and the base of the giant branch with observed results from Nordlander et al. (2012) of the globular cluster NGC 6397. However, as it still remains an adjustable parameter in the stellar models, it could be incorrect. Adopting a much smaller value results in a maximum change in surface metallicity during main sequence evolution of up to 1.0 dex (see e.g. Fig. 5 which shows a maximum change of 0.55 dex). This would require an initial value of the order of –1.40 dex. This scenario is improbable and would result in a surface metallicity difference between turnoff and giant stars in disagreement with Nordlander et al. (2012). Adopting a larger value for the parameter would reduce the maximum change in surface metallicity during evolution, and converging on a surface value of –2.10 dex would require the initial value to be decreased, for example to –1.90 dex or Z_{i}/X_{i} = 0.0002. For a star with an age almost the age of the Galaxy (adopting zero reddening) a lower initial metallicity could be more likely. However, a star 1 billion years younger could be consistent with a higher initial metallicity. The impact of this change of physics and metallicity on the resulting stellar properties is currently limited by our observational errors. At the moment we cannot explore this further and we are required to impose the external constraints.
8.4. Improving the precision and accuracy of the age
It is clear that in order to determine the most accurate and precise age for HD 140283, it is necessary to solve the extinction problem, and to reduce the span of masses that pass through the error box. The second question can be addressed by obtaining a more precise determination of the angular diameter. From Fig. 8 one can immediately see that by reducing the error bar in T_{eff} by a factor of two, fewer models would satisfy the constraints. As an example, the right panel shows that models of masses between 0.790 and 0.805 (a total span of 0.015 M_{⊙}) would be the only models that pass through the error box of half of its size. The parallax and F_{bol} already contribute very little to the error bar for R and T_{eff}, and so the only option is to obtain a precision in θ of the order of 0.007 mas. By cutting both the external and the statistical errors in half, this precision could be achieved. The current precision is just under 4%, which is excellent considering that we are already working at the limits of angular resolution using the CHARA array with VEGA and VEGA’s sensitivity limits. However, the telescopes of the CHARA array will be equipped with adaptative optic (AO) systems in the next one or two years. This will improve both the sensitivity and measurement precision of VEGA. The gain in sensitivity would allow us to observe fainter (so smaller) calibrators, hence reduce the external error affecting the calibrated visibilities. New observations of HD140283 with CHARA/VEGA or CHARA/FRIEND (future instrument) equipped with AO would allow us to achieve the necessary precision.
Fig. 11 Relative difference in the sound speed profile between models M3 and M1 (continuous) and M3 and M4 (dashed). 

Open with DEXTER 
8.5. Improving the mass and age from asteroseismology
Another way to determine the mass of the star is through the detection and interpretation of stellar oscillation frequencies. Acoustic oscillation frequencies (sound waves) are sensitive to the sound speed profile of the star and thus its density structure and mean density. Because R has been measured, indepth asteroseismic analysis would allow a very high precision modeldependent determination of the mass. In Fig. 11 we show the relative differences between the sound speed profiles (from 0.02 to 0.8 the stellar radius) for the M3 and M1 models, (M3−M1)/M3. This represents the difference between the reddened and unreddened solution for Y_{i} = 0.245 (continuous). At the same time, we also show the difference between the reddened solutions with different initial helium abundance, M3 and M4 (dashed). The largest differences for both cases are found very close to the centre of the star and near the outer convection zone. These differences cause changes in the frequencies of the oscillation modes that can be detected, e.g. a radial highorder mode (l = 0,n = 24) has a predicted frequency of 1240.2 μHz for M3 and this same mode for M1 is 1251.0 μHz, although for the M4 model (different Y_{i}) it is 1240.7 μHz, not distinguishable from the M3 solution at current asteroseismic modelling limitations (e.g. Kjeldsen et al. 2008b). Precision in asteroseismic frequencies are well below 1.0 μHz, even for a tenday groundbased observational campaign. Stellar oscillations have been detected in stars of this mass. At this stage of evolution, the socalled gmodes are very sensitive to the core conditions and hence should be able to distinguish even more clearly between both models. However, we need to verify whether there are sufficient metallic lines in the spectrum to accurately detect radial velocity shifts (the oscillations).
8.5.1. Detection of oscillations in metalpoor stars
The star ν Ind has a similar mass to HD 140283, but has a higher metallicity ([M/H] = –1.3). This star was the subject of a radial velocity asteroseismic campaign using telescopes in Australia and Chile and solarlike oscillations were detected (Bedding et al. 2006; Carrier et al. 2007). The oscillations were analysed and the amplitudes of the frequencies were in the range 50−170 cm s^{1} (the solar amplitudes are of the order of 23 cm s^{1}, see also Kjeldsen et al. (2008a) for a more recent analysis of the amplitudes of both of these stars). A recent asteroseismic observational campaign of a different star using the SOPHIE spectrograph at the Observatoire de Haute Provence in February 2013 allowed us to reach a precision of 20−25 cm s^{1} in the power spectrum for a total of 8 good nights out of 12 assigned. The detections of such oscillation frequencies are clearly very feasible using different instruments or telescopes. Oscillations are detectable in metalpoor stars, but how metalpoor is observationally possible? HD 140283 would be an ideal target for continued observations, because the nondetection of oscillations would impose an upper limit on the amplitudes of the modes in metalpoor stars, while the detection of oscillation frequencies would provide valuable constraints on the mass and age of this star, and stellar interior models.
9. Conclusions
In this work we determined stellar properties of HD 140283 by measuring the angular diameter of the star and performing bolometric flux fitting. The measured diameter is θ_{1D} = 0.353 ± 0.013 mas and coupling this with a parallax of π = 17.15 ± 0.14 mas from Bond et al. (2013) we obtain a linear radius of the star R = 2.21 ± 0.08R_{⊙} and log g = 3.65 ± 0.06. The bolometric flux fitting resulted in two solutions, one where interstellar extinction A_{V} = 0.0 mag, and the other with a maximum nonzero value of 0.1 mag. Adopting these two values we derived L (4.12 or 4.47 L_{⊙}), T_{eff} (5534 or 5647 K), bolometric corrections, and the absolute magnitude M_{V} (Table 6). We also determined a spectroscopic T_{eff} using a 1D LTE analysis fitting Hα line wings and we found a value of 5626 K, a value more compatible with the interferometric T_{eff} assuming a small amount of reddening. If F_{bol} is indeed a higher value because of interstellar reddening, this is a very important result as it shows that a 1D LTE analysis of Hα wings yields accurate T_{eff} at this lowmetallicity regime. Our measured angular diameter is larger by about 2σ than that predicted by the surface brightness relations of Kervella et al. (2004), suggesting that indeed more determinations of θ are needed for such metalpoor and/or evolved objects to more accurately calibrate this scale, especially important with the advent of Gaia and with the availability of precise distance measurements for a billion stars in the next five years.
By interpreting the observational results using stellar models we derived a strict relationship between mass M, mixinglength parameter α, and initial helium abundance Y_{i} (Eqs. (5)–(10)), where only two are independent. We further determined the age t of the star for all possible combinations of these parameters; Eq. (11) parametrizes t as a function of the adopted M, Y_{i} (implying α), and reddening A_{V}. Fixing α = 1.0 we studied the impact of changing Y_{i} and A_{V} on the final solution. For A_{V} = 0.10 mag we obtain M = 0.805 ± 0.010M_{⊙}, corresponding to an age of 12.17 ± 0.69 Gyr. For A_{V} = 0.00 mag we obtain a slightly lower mass of M = 0.780 ± 0.010M_{⊙}, corresponding to an age of 13.68 ± 0.61 Gyr (fixing Y_{i} = 0.245). The impact of Y_{i} is much smaller, causing a decrease in M and t of 0.01 M_{⊙} and 0.2 Gyr. Table 12 summarizes some of our modelling results, while the online version gives the full range of parameter solutions.
By adopting zero or very little reddening Bond et al. (2013) and VandenBerg et al. (2014) derive an age of ~14.3 Gyr in a completely independent manner. These results were obtained by adopting an oxygen abundance which was derived with a higher T_{eff} than the one determined in this work. In their work, the mixinglength parameter was not explored in their models and this would shift the evolution tracks towards lower T_{eff}, thus providing a match to the observational data at a lower age. Casagrande et al. (2011) in an independent manner derive an age of 13.8 Gyr and 13.9 Gyr assuming no reddening and using the BaSTI (Cordier et al. 2007) and PADOVA (Bertelli et al. 2008) isochrones. These determinations of the age of HD 140283 give an indication of the systematic error we can expect on stellar ages by using different approaches and assumptions on both the data and models.
For all of our stellar models we needed to adapt the mixinglength parameter to a value much lower than the Sun, i.e. from 2.0 to 1.0 using the Eggleton (1972) treatment of mixinglength. If reddening is nonexistent, then α takes a lower value (0.5 ≤ α ≤ 1.0) on average. This result reinforces previously published results of the nonuniversality of the mixinglength parameter, e.g. Miglio & Montalbán (2005), one that is generally fixed in published stellar evolution tracks/isochrones. We note that our somewhat lower T_{eff} (compared to some recent determinations of 5700 K and higher) along with the observed Li abundances provides an important reference point for studying Li versus T_{eff} for metalpoor stars (see e.g. Charbonnel & Primas 2005, their Fig. 12) and consequently Li surface depletion scenarios.
The next crucial steps to better understand this star and to decorrelate the stellar properties of mass, age, initial helium abundance, and to look at the treatment of diffusion in such stars, is to reduce the uncertainty on the angular diameter by at least a factor of two, and to perform asteroseismic observations to better recover the internal profile of the star and hence its mass and age. Predictions of the mixinglength parameter from external constraints would also lead to a large improvement in determining mass, initial helium, and age. All of these are realistic with current numerical capabilities and current or nearfuture instrumentation. One of the biggest improvements to the age and mass determinations, however, relies on the exact knowledge of extinction between the Earth and this star.
The projected baseline Bp depends on a number of factors including the physical distance between the telescope pair (baseline), the angles of the baseline compared to the northeast direction and the altitude of the object over the horizon, see e.g. Nardetto, N., Ph.D. Thesis, 2005, p. 57; http://tel.archivesouvertes.fr
Acknowledgments
This paper is dedicated to the wonderful astronomer Olivier Chesneau. We thank Isabelle TallonBosc for detailed discussions regarding uncertainties in interferometric data. We thank the VEGA observing group at the Observatoire de la Côte d’Azur (OCA) in Nice, France (remote observing) and at CHARA, California, USA. The work by the SAM (http://www.astro.uu.se/~ulrike/GaiaSAM/) group in defining a set of benchmark stars is greatly appreciated. The CHARA Array is funded by the National Science Foundation through NSF grant AST0908253 and by Georgia State University through the College of Arts and Sciences. T.S.B. acknowledges support provided through NASA grant ADAP120172. U.H. acknowledges support from the Swedish National Space Board (Rymdstyrelsen). P.K. acknowledges the support of PHASE, the high angular resolution partnership between ONERA, Observatoire de Paris, CNRS and the Université Denis Diderot, Paris 7. We acknowledge financial support from the “Programme National de Physique Stellaire” (PNPS) of CNRS/INSU, France. During part of this research OLC was a Henri Poincaré Fellow at OCA, which is funded by the Conseil Général des AlpesMaritimes and OCA.
References
 Alexander, D. R., & Ferguson, J. W. 1994, ApJ, 437, 879 [NASA ADS] [CrossRef] [Google Scholar]
 Ali, A. W., & Griem, H. R. 1966, Phys. Rev., 144, 366 [NASA ADS] [CrossRef] [Google Scholar]
 Allard, N. F., Kielkopf, J. F., Cayrel, R., & van’t VeerMenneret, C. 2008, A&A, 480, 581 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Alonso, A., Arribas, S., & MartinezRoger, C. 1996, A&AS, 117, 227 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Angulo, C. 1999, AIP Conf. Ser., 495, 365 [NASA ADS] [CrossRef] [Google Scholar]
 Arenou, F., Grenon, M., & Gomez, A. 1992, A&A, 258, 104 [NASA ADS] [Google Scholar]
 Aurière, M. 2003, in EAS Pub. Ser. 9, eds. J. Arnaud, & N. Meunier, 105 [Google Scholar]
 Bagnulo, S., Jehin, E., Ledoux, C., et al. 2003, The Messenger, 114, 10 [NASA ADS] [Google Scholar]
 BailerJones, C. A. L., Andrae, R., Arcay, B., et al. 2013, A&A, 559, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Baker, J. 2008, Natl. Inst. Stand. Technol. Technical Note, 1612 [Google Scholar]
 Ballester, P., Modigliani, A., Boitquin, O., et al. 2000, The Messenger, 101, 31 [NASA ADS] [Google Scholar]
 Barklem, P. S., Piskunov, N., & O’Mara, B. J. 2000, A&A, 355, L5 [NASA ADS] [Google Scholar]
 Bedding, T. R., Butler, R. P., Carrier, F., et al. 2006, ApJ, 647, 558 [NASA ADS] [CrossRef] [Google Scholar]
 Bennett, C. L., Larson, D., Weiland, J. L., et al. 2013, ApJS, 208, 20 [NASA ADS] [CrossRef] [Google Scholar]
 Bertelli, G., Girardi, L., Marigo, P., & Nasi, E. 2008, A&A, 484, 815 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bigot, L., Kervella, P., Thévenin, F., & Ségransan, D. 2006, A&A, 446, 635 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bigot, L., Mourard, D., Berio, P., et al. 2011, A&A, 534, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 BlancoCuaresma, S., Soubiran, C., Jofré, P., & Heiter, U. 2014, A&A, 566, A98 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bonaca, A., Tanner, J. D., Basu, S., et al. 2012, ApJ, 755, L12 [NASA ADS] [CrossRef] [Google Scholar]
 Bond, H. E., Nelan, E. P., VandenBerg, D. A., Schaefer, G., & Harmer, D. 2013, in Am. Astron. Soc. Meet. Abstr., 221, 443.08 [Google Scholar]
 Bonifacio, P., & Molaro, P. 1998, ApJ, 500, L175 [NASA ADS] [CrossRef] [Google Scholar]
 Bonneau, D., Clausse, J.M., Delfosse, X., et al. 2006, A&A, 456, 789 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Boyajian, T. S., van Belle, G., & von Braun, K. 2014, AJ, 147, 47 [NASA ADS] [CrossRef] [Google Scholar]
 Brown, T. M., & Gilliland, R. L. 1994, ARA&A, 32, 37 [NASA ADS] [CrossRef] [Google Scholar]
 Burgers, J. M. 1969, Flow Equations for Composite Gases (Academic Press Inc) [Google Scholar]
 Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 [NASA ADS] [CrossRef] [Google Scholar]
 Carney, B. W., & Latham, D. W. 1987, AJ, 93, 116 [NASA ADS] [CrossRef] [Google Scholar]
 Carrier, F., Kjeldsen, H., Bedding, T. R., et al. 2007, A&A, 470, 1059 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Casagrande, L., Flynn, C., Portinari, L., Girardi, L., & Jimenez, R. 2007, MNRAS, 382, 1516 [NASA ADS] [CrossRef] [Google Scholar]
 Casagrande, L., Ramírez, I., Meléndez, J., Bessell, M., & Asplund, M. 2010, A&A, 512, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Casagrande, L., Schönrich, R., Asplund, M., et al. 2011, A&A, 530, A138 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Casagrande, L., Portinari, L., Glass, I. S., et al. 2014, MNRAS [Google Scholar]
 Cayrel, R., van’t VeerMenneret, C., Allard, N. F., & Stehlé, C. 2011, A&A, 531, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Charbonnel, C., & Primas, F. 2005, A&A, 442, 961 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chiavassa, A., Plez, B., Josselin, E., & Freytag, B. 2009, A&A, 506, 1351 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [Google Scholar]
 Chiavassa, A., Bigot, L., Kervella, P., et al. 2012, A&A, 540, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Claret, A., Hauschildt, P. H., & Witte, S. 2012, A&A, 546, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Collet, R., Nordlund, Å., Asplund, M., Hayek, W., & Trampedach, R. 2009, Mem. Soc. Astron. It., 80, 719 [NASA ADS] [Google Scholar]
 Cordier, D., Pietrinferni, A., Cassisi, S., & Salaris, M. 2007, AJ, 133, 468 [NASA ADS] [CrossRef] [Google Scholar]
 Cousins, A. W. J. 1984, South African Astronomical Observatory Circular, 8, 59 [NASA ADS] [Google Scholar]
 Creevey, O. L., Thévenin, F., Boyajian, T. S., et al. 2012, A&A, 545, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Creevey, O. L., Thévenin, F., Basu, S., et al. 2013, MNRAS, 431, 2419 [NASA ADS] [CrossRef] [Google Scholar]
 Dekker, H., D’Odorico, S., Kaufer, A., Delabre, B., & Kotzlowski, H. 2000, in SPIE Conf. Ser. 4008, eds. M. Iye, & A. F. Moorwood, 534 [Google Scholar]
 Donati, J., Semel, M., Carter, B. D., Rees, D. E., & Collier Cameron, A. 1997, MNRAS, 291, 658 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Drozdowski, R., Ignaciuk, M., Kwela, J., & Heldt, J. 1997, Z. Phys. D, 41, 125 [NASA ADS] [CrossRef] [Google Scholar]
 Eggleton, P. P. 1972, MNRAS, 156, 361 [NASA ADS] [Google Scholar]
 Eggleton, P. P., Faulkner, J., & Flannery, B. P. 1973, A&A, 23, 325 [NASA ADS] [Google Scholar]
 Epstein, C. R., Elsworth, Y. P., Johnson, J. A., et al. 2014, ApJ, 785, L28 [NASA ADS] [CrossRef] [Google Scholar]
 Flower, P. J. 1996, ApJ, 469, 355 [NASA ADS] [CrossRef] [Google Scholar]
 Fuhrmann, K. 1998, A&A, 330, 626 [NASA ADS] [Google Scholar]
 Fukugita, M., Shimasaku, K., & Ichikawa, T. 1995, PASP, 107, 945 [NASA ADS] [CrossRef] [Google Scholar]
 Gallagher, A. J., Ryan, S. G., García Pérez, A. E., & Aoki, W. 2010, A&A, 523, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gilmore, G., Randich, S., Asplund, M., et al. 2012, The Messenger, 147, 25 [NASA ADS] [Google Scholar]
 González Hernández, J. I., & Bonifacio, P. 2009, A&A, 497, 497 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gray, D. F. 2005, The Observation and Analysis of Stellar Photospheres (Cambridge University Press) [Google Scholar]
 Grevesse, N., & Noels, A. 1993, in Origin and Evolution of the Elements, eds. N. Prantzos, E. VangioniFlam, & M. Casse, 15 [Google Scholar]
 Griersmith, D. 1980, AJ, 85, 789 [NASA ADS] [CrossRef] [Google Scholar]
 Grundahl, F., VandenBerg, D. A., Bell, R. A., Andersen, M. I., & Stetson, P. B. 2000, AJ, 120, 1884 [NASA ADS] [CrossRef] [Google Scholar]
 Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2003, in Stellar Atmosphere Modeling, eds. I. Hubeny, D. Mihalas, & K. Werner, ASP Conf. Ser., 288, 331 [Google Scholar]
 Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hauck, B., & Mermilliod, M. 1990, A&AS, 86, 107 [NASA ADS] [Google Scholar]
 Hauck, B., Nitschelm, C., Mermilliod, M., & Mermilliod, J. C. 1995, VizieR Online Data Catalog: II/167 [Google Scholar]
 Hauschildt, P. H., Allard, F., Ferguson, J., Baron, E., & Alexander, D. R. 1999, ApJ, 525, 871 [NASA ADS] [CrossRef] [Google Scholar]
 Haywood, M., Di Matteo, P., Lehnert, M. D., Katz, D., & Gómez, A. 2013, A&A, 560, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Houdashelt, M. L., Bell, R. A., & Sweigart, A. V. 2000, AJ, 119, 1448 [NASA ADS] [CrossRef] [Google Scholar]
 Huber, D., Ireland, M. J., Bedding, T. R., et al. 2012, ApJ, 760, 32 [NASA ADS] [CrossRef] [Google Scholar]
 Husser, T.O., Wendevon Berg, S., Dreizler, S., et al. 2013, A&A, 553, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Jofré, P., Heiter, U., Soubiran, C., et al. 2014, A&A, 564, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kervella, P., Thévenin, F., Di Folco, E., & Ségransan, D. 2004, A&A, 426, 297 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kervella, P., Mérand, A., Pichon, B., et al. 2008, A&A, 488, 667 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kjeldsen, H., & Bedding, T. R. 1995, A&A, 293, 87 [NASA ADS] [Google Scholar]
 Kjeldsen, H., Bedding, T. R., Arentoft, T., et al. 2008a, ApJ, 682, 1370 [NASA ADS] [CrossRef] [Google Scholar]
 Kjeldsen, H., Bedding, T. R., & ChristensenDalsgaard, J. 2008b, ApJ, 683, L175 [NASA ADS] [CrossRef] [Google Scholar]
 Korn, A. J., Grundahl, F., Richard, O., et al. 2007, ApJ, 671, 402 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
 Kurucz, R. L. 2010, Robert L. Kurucz online database of observed and predicted atomic transitions, http://kurucz.harvard.edu/atoms/2200/gfemq2200.pos, http://kurucz.harvard.edu/atoms/2201/gfemq2201.pos [Google Scholar]
 Lallement, R., Vergely, J.L., Valette, B., et al. 2014, A&A, 561, A91 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lebreton, Y. 2000, ARA&A, 38, 35 [NASA ADS] [CrossRef] [Google Scholar]
 Lebreton, Y., & Goupil, M. J. 2014, A&A, 569, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lejeune, T., Cuisinier, F., & Buser, R. 1997, VizieR Online Data Catalog: J/A&AS/125/229 [Google Scholar]
 Masana, E., Jordi, C., & Ribas, I. 2006, A&A, 450, 735 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mathur, S., Metcalfe, T. S., Woitaszek, M., et al. 2012, ApJ, 749, 152 [NASA ADS] [CrossRef] [Google Scholar]
 May, M., Richter, J., & Wichelmann, J. 1974, A&AS, 18, 405 [NASA ADS] [Google Scholar]
 Mayor, M., Pepe, F., Queloz, D., et al. 2003, The Messenger, 114, 20 [NASA ADS] [Google Scholar]
 McMillan, R. S., Breger, M., Ferland, G. J., & Loumos, G. L. 1976, PASP, 88, 495 [NASA ADS] [CrossRef] [Google Scholar]
 Meléndez, J., Ramírez, I., Casagrande, L., et al. 2010, Ap&SS, 328, 193 [NASA ADS] [CrossRef] [Google Scholar]
 Mermilliod, J.C. 1986, Catalogue of Eggen’s UBV data [Google Scholar]
 Metcalfe, T. S., Creevey, O. L., & ChristensenDalsgaard, J. 2009, ApJ, 699, 373 [NASA ADS] [CrossRef] [Google Scholar]
 Metcalfe, T. S., Chaplin, W. J., Appourchaux, T., et al. 2012, ApJ, 748, L10 [NASA ADS] [CrossRef] [Google Scholar]
 Metcalfe, T. S., Creevey, O. L., Doğan, G., et al. 2014, ApJS, 214, 27 [NASA ADS] [CrossRef] [Google Scholar]
 Miglio, A., & Montalbán, J. 2005, A&A, 441, 615 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Morel, M., & Magnenat, P. 1978, A&AS, 34, 477 [NASA ADS] [Google Scholar]
 Morel, P. 1997, A&AS, 124, 597 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
 Morel, P., & Lebreton, Y. 2008, Ap&SS, 316, 61 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Morel, P., & Thévenin, F. 2002, A&A, 390, 611 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mourard, D., Bério, P., Perraut, K., et al. 2011, A&A, 531, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mourard, D., Challouf, M., Ligi, R., et al. 2012, in SPIE Conf. Ser., 8445 [Google Scholar]
 Nordlander, T., Korn, A. J., Richard, O., & Lind, K. 2012, ApJ, 753, 48 [NASA ADS] [CrossRef] [Google Scholar]
 Norris, J., Bessell, M. S., & Pickles, A. J. 1985, ApJS, 58, 463 [NASA ADS] [CrossRef] [Google Scholar]
 O’Donnell, J. E. 1994, ApJ, 422, 158 [NASA ADS] [CrossRef] [Google Scholar]
 Perryman, M. A. C. 2005, in Astrometry in the Age of the Next Generation of Large Telescopes, eds. P. K. Seidelmann, & A. K. B. Monet, ASP Conf. Ser., 338, 3 [Google Scholar]
 Pickles, A. J. 1998a, VizieR Online Data Catalog: J/PASP/110/863 [Google Scholar]
 Pickles, A. J. 1998b, PASP, 110, 863 [NASA ADS] [CrossRef] [Google Scholar]
 Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical recipes in FORTRAN, The art of scientific computing (Cambridge University Press) [Google Scholar]
 Ramírez, I., Collet, R., Lambert, D. L., Allen de Prieto, C., & Asplund, M. 2010, ApJ, 725, L223 [NASA ADS] [CrossRef] [Google Scholar]
 Rogers, F. J., & Iglesias, C. A. 1992, ApJS, 79, 507 [NASA ADS] [CrossRef] [Google Scholar]
 Ruchti, G. R., Bergemann, M., Serenelli, A., Casagrande, L., & Lind, K. 2013, MNRAS, 429, 126 [NASA ADS] [CrossRef] [Google Scholar]
 Ryan, S. G. 1989, AJ, 98, 1693 [NASA ADS] [CrossRef] [Google Scholar]
 Salaris, M., Chieffi, A., & Straniero, O. 1993, ApJ, 414, 580 [NASA ADS] [CrossRef] [Google Scholar]
 Schuster, W. J., & Nissen, P. E. 1989, A&A, 221, 65 [NASA ADS] [Google Scholar]
 Silva Aguirre, V., Basu, S., Brandão, I. M., et al. 2013, ApJ, 769, 141 [NASA ADS] [CrossRef] [Google Scholar]
 Siqueira Mello, C., Barbuy, B., Spite, M., & Spite, F. 2012, A&A, 548, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 [NASA ADS] [CrossRef] [Google Scholar]
 Soubiran, C., Le Campion, J.F., Cayrel de Strobel, G., & Caillo, A. 2010, A&A, 515, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Sturmann, J., ten Brummelaar, T., Sturmann, L., & McAlister, H. A. 2010, in SPIE Conf. Ser., 7734 [Google Scholar]
 ten Brummelaar, T. A., McAlister, H. A., Ridgway, S. T., et al. 2005, ApJ, 628, 453 [NASA ADS] [CrossRef] [Google Scholar]
 Thévenin, F. 1998, VizieR Online Data Catalog:III/193 [Google Scholar]
 Torres, G. 2010, AJ, 140, 1158 [NASA ADS] [CrossRef] [Google Scholar]
 Trampedach, R., & Stein, R. F. 2011, ApJ, 731, 78 [NASA ADS] [CrossRef] [Google Scholar]
 Upgren, A. R., Zhao, H., & Lee, J. T. 1992, AJ, 104, 1203 [NASA ADS] [CrossRef] [Google Scholar]
 Valenti, J. A., & Fischer, D. A. 2005, ApJS, 159, 141 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Valenti, J. A., & Piskunov, N. 1996, A&AS, 118, 595 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 van Belle, G. T., van Belle, G., CreechEakman, M. J., et al. 2008, ApJS, 176, 276 [NASA ADS] [CrossRef] [Google Scholar]
 van Leeuwen, F. 2007, A&A, 474, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 VandenBerg, D. A., Bond, H. E., Nelan, E. P., et al. 2014, ApJ, 792, 110 [NASA ADS] [CrossRef] [Google Scholar]
 von Braun, K., Boyajian, T. S., Kane, S. R., et al. 2012, ApJ, 753, 171 [NASA ADS] [CrossRef] [Google Scholar]
 Wahlgren, G. M. 2010, in EAS Pub. Ser. 43, eds. R. Monier, B. Smalley, G. Wahlgren, & P. Stee, 91 [Google Scholar]
 Wiese, W. L., & Fuhr, J. R. 2009, J. Phys. Chem. Ref. Data, 38, 565 [NASA ADS] [CrossRef] [Google Scholar]
 Yamada, S., Suda, T., Komiya, Y., Aoki, W., & Fujimoto, M. Y. 2013, MNRAS, 436, 1362 [NASA ADS] [CrossRef] [Google Scholar]
 Yi, S., Demarque, P., Kim, Y.C., et al. 2001, ApJS, 136, 417 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Bolometric flux fitting approach 2 (method 2B and 2C)
We describe in detail the second bolometric fluxfitting approach that was implemented in this work. The method incorporates a nonlinear leastsquares minimization algorithm (LevenbergMarqwardt) to find the best scaled interpolated spectrum that fits a set of observed flux points or photometric magnitudes (in this work we just use the flux data). The fitting method requires on input a set of initial parameters – T_{eff}, [Fe/H], log g, E(B − V), and a scaling factor θ_{s} – which define the characteristics of the spectra, the reddening to apply, and ratio of the stellar radius to the star’s distance. These parameters are optimized within the algorithm and the bolometric flux is calculated by numerical integration of the scaled optimal spectrum. The method is designed to allow implementation of any spectral library. In this work we use the BASEL (Lejeune et al. 1997) and PHOENIX libraries (Hauschildt et al. 1999; Husser et al. 2013).
Because of the strong degeneracy between the stellar parameters, in practice we fix three of the values T_{eff}, log g, [Fe/H], A_{V}, while fitting the fourth along with θ_{s}, so only two stellar parameters are fitted at any one time. Changing log g and [Fe/H] by 0.5 dex, which is much larger than the typical errors on these values, has a negligible influence on the final calculated bolometric flux F_{bol}. Generally these parameters are the ones that are fixed. Changing A_{V}, however, does have an important effect as explained below.
Within the fitting process, each spectrum is linearly interpolated between adjacent parameter points to the desired parameter values as proposed by the minimization algorithm. We have tested for different interpolation schemes and because of the relatively dense grid spacing (e.g. 100 K or 250 K in T_{eff}, 0.5 in log g, 0.5 in [Fe/H]) the choice and order of interpolation is irrelevant. The scale factor θ_{s} is also proposed at the same time by the minimization algorithm. A value is then calculated between the observed data and the newly proposed scaled spectrum. If the value improves then the minimization algorithm continues in the same direction. When the fails to improve after four iterations the fitting stops, but convergence is usually reached within only three or four iterations. The final results are not sensitive to the initial parameters; for example, initial T_{eff} values differing by 1000 K yielded consistent results.
To account for reddening we apply reddening laws to the spectra using the Interactive Data Language (IDL) routine ccm_unred and the input E(B − V), along with a choice of extinction laws of R_{V}, which we fix at 3.1 in this work. This routine is based on the Cardelli et al. (1989) coefficients and updated with O’Donnell (1994) in the UV. The values of T_{eff} and A_{V} are highly degenerate parameters if no external constraints are imposed, and so to determine A_{V} a fitting process is begun allowing only T_{eff} and θ_{s} to be fitted, but imposing a range of values of E(B − V). The fit yielding the smallest defines the adopted E(B − V) and thus A_{V} (using R_{V}). The same results are obtained by fixing T_{eff} at different values while fitting E(B − V). If the observed data give a best match to a reddened spectrum, the intrinsic (unreddened) flux of the star is thus larger, indicating a more luminous star.
All Tables
Observation log for HD 140283 using the VEGA interferometer on the CHARA array in threetelescope mode (telescopes E1, W1, and W2).
Bolometric flux and interstellar reddening from fitting literature photometryconvertedtoflux measurements to the Pickles (A), BASEL (B), and PHOENIX (C) libraries of stellar spectra using two independent approaches (1 and 2).
Stellar properties of HD 140283 derived in this paper by adopting π = 17.15 ± 0.14 mas from Bond et al. (2013).
Information on observed spectra used for the fit: instrument/archive name, date of observation, mean S/N.
Atomic data for spectral lines contributing to flux in Line mask 3 for HD 140283.
Bestfitting T_{eff} values rounded to 10 K and reduced chisquare () for the different observations and masks for HD 140283.
All Figures
Fig. 1 Top: squared visibility measurements and the fitted visibility curve corresponding to a 1D diameter θ_{1D} = 0.353 mas. Bottom: data residuals scaled by the errors on the measurements. The symbols represent the night the data were taken: ⋄ = 2012 April 18, °= 2012 May 21, □= 2014 May 03, and ⋆ = 2014 May 04 (Sect. 2.2.2). 

Open with DEXTER  
In the text 
Fig. 2 A Gaussian curve representing the fitted angular diameter and the 68% probability region (± σ_{θ,STAT}) due to the statistical errors (light grey shaded region). The vertical line illustrates the centre of the Gaussian corresponding to the fitted diameter θ. The inner distribution shows fitted angular diameters from 10 000 simulations by randomly adding the external error to the measurements. The dark grey shaded region illustrates the 68% probability region due to the external errors (± σ_{θ,EXT}, Sect. 2.2.2). 

Open with DEXTER  
In the text 
Fig. 3 Photometric measurements converted to absolute flux for HD 140283, with an example of a fit to a model spectrum (Sect. 3). 

Open with DEXTER  
In the text 
Fig. 4 Top: observed NARVAL spectrum (black) of the Hα profile of HD 140283 and fitted spectrum (red) using the three narrow continuum masks (shaded regions) for log g = 3.65 and [Fe/H] = –2.46. Bottom: residuals between the observed and fitted spectrum for each line mask, with the black dots indicating the line masks. See Sect. 5 for details. 

Open with DEXTER  
In the text 
Fig. 5 The evolution of the observed (surface) metal abundance, shown as a function of T_{eff}, for a 0.80 M_{⊙} model with Y_{i} = 0.245 including the effects of diffusion and extra mixing for two initial metalhydrogen ratio mass fractions (Z_{i}/X_{i} = 0.0005, 0.0007). At T_{eff} = 5600 K, we highlight the values of the luminosity for both of the models. See Sect. 7.1 for details. 

Open with DEXTER  
In the text 
Fig. 6 Left: evolutionary tracks for models of different masses (labelled by their solar masses) for Y_{i} = 0.245 and Z_{i}/X_{i} = 0.0005 plotted as a function of age. The dashed lines illustrate changing the M = 0.80M_{⊙} model in Y_{i} (ΔY) and Z_{i}/X_{i} (ΔZ). The unreddened and reddened luminosity constraints are indicated by the error bars. Right: approximate correlations between M and the initial helium abundance Y_{i} assuming an age >12 Gyr (continuous) and then extending it to >10 Gyr (dashed), given the luminosity constraint. See Sect. 7.2 for details. 

Open with DEXTER  
In the text 
Fig. 7 Correlations between M and α for the reddened solution (error box) imposing Y_{i} = 0.245, illustrating the plausible ranges for these parameters imposed by stellar evolution. Similar colours have the same mass and the continuous lines have α = 1.0. See Sect. 7.3.1 for details. 

Open with DEXTER  
In the text 
Fig. 8 HR diagrams illustrating the error box corresponding to the observations derived in this work for HD 140283. The left and right panels show the results when we adopt A_{V} = 0.00 and 0.10 mag. The dotted box in the left panel shows the error box when the distance is increased by 1σ, also illustrated by the arrow. In the right panel (nonzero reddening solution) we also show the A_{V} = 0.00 error box (dotted box), located in the lower right part of the figure. The central models correspond to Y_{i} = 0.245,α = 1.00, and Z/X_{i} = 0.0005 (with Re_{ν} = 5). The optimal masses are M = 0.780 and M = 0.805M_{⊙}, and are shown as the blue continuous lines passing near the centre of the error boxes. We also illustrate other models by changing the mass (blue) and the other parameters (green). The red dotted line corresponds to the radius derived in this paper. See Sect. 7.3 for details. 

Open with DEXTER  
In the text 
Fig. 9 Relationship between the age and the mass of the models that satisfies the observational constraints. The circles/diamonds correspond to the Y_{i} = 0.245 / 0.255 solution, and the red denotes the A_{V} = 0.10 mag solution. The crosses denote the values corresponding to the centres of the error boxes in the HR diagram (Fig. 8) and the solutions given in Table 12. See Sect. 7.3.4 for details. 

Open with DEXTER  
In the text 
Fig. 10 Spectroscopic determinations of T_{eff} and [Fe/H] from the literature (coloured symbols) and our spectroscopic determination using Hα line wings (black circle). The shaded area illustrates the interferometric determination of T_{eff} from the nonzero reddening solution, while the dashed lines illustrate the boundaries of the interferometric T_{eff} from the zero reddening solution. 

Open with DEXTER  
In the text 
Fig. 11 Relative difference in the sound speed profile between models M3 and M1 (continuous) and M3 and M4 (dashed). 

Open with DEXTER  
In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.