Arxiv:Astro-Ph/0603219V2 7 Apr 2006

arXiv:astro-ph/0603219v2 7 Apr 2006 o 94 oia oln&Csel 96 eiuh & Sekiguchi 1996; Castelli “yellow” ( a & from to Meléndez Ram´ırez 2005) Bohlin Mod- & spread 2000; Fukugita Colina, large 1994; 1996). a lor Strobel show ( “blue” de still a (Cayrel (1994-2006) values values 0.685 shows ern to (1957-1982) literature 0.615 early from the which for g;adalae ihteSmnroPraet eAstronom de Permanente Seminario the with affiliated and ogy; oa“e”( (B-V) “red” a to h aionaIsiueo ehooy h nvriyo C a of partnership University NASA. the the scientific Technology, and of a Institute as California operated the is which Observatory, o n tr.Wt htmti tblt fjust of photometry; stability high-precision results photometric of 7%) capable a to is With IRAC up stars. 1%, (by K over inconsistent and yielded (Reach A Telescope 2005) for Space al. Spitzer Array the et Infrared the on of (IRAC) exam- calibration Camera For absolute recent systems. the photometric ple, space and terrestrial of inisEpcae fteUiesddNcoa ao eS de Mayor Nacional Universidad the Peru. of Marcos, Espaciales Ciencias y a,w ontytko cuaeytesolar particu- the accurately In know the yet not temperature. as do we to such lar, color of astrophysics, from calibration in transformation the for observables important which fundamental at- known, extremely well them their very makes tem- be age), (e.g., should mass, parameters parameters mospheric metallicity, fundamental luminosity, their in perature, Sun only the are to that present. chemicals is be of life probe may when signs will twins) produced that for missions (solar planets space Sun future Earth-like for the planets. targets habitable to best host the identical can Sun stars our Thus, like stars that proof planets? Earth-like for is: surveys planet lcrncades jorge,katie,[email protected] address: Electronic ≈ 2 1 wn ol eueu nteaslt u calibration flux absolute the in useful be could Twins ic oa wn r eysmlr(dal,identical) (ideally, similar very are twins solar Since terrestrial extrasolar for question important most The rpittpstuigL using Melendez typeset J. Preprint to: proofs Send .4 tazs&Vlag 94 omege l 2006) al. et Holmberg 1994; Valiauga & Straizys 0.64: loDprmn fAtooy aionaIsiueo Tech of Institute California Astronomy, of Department Also h aapeetdhri eeotie tteWM Keck W.M. the at obtained were herein presented data The ⊙ eerhSho fAtooy&Atohsc,M.SrmoOb Stromlo Mt. Astrophysics, & Astronomy of School Research hs tr hudb ie o roiyi xpae n EIsurve hos SETI the and habit to exoplanet their in percent) priority within few top planet headings: a given Subject terrestrial (within be a should similar hosting stars very the of these are that likelihood and Sco the twin 18 solar including second-closest and the 98618 is HD 98618 HD that show We hwdta 8Soi lotaslrti.I h rsn ok eus we work, present the In twin. solar by ago a decade almost a 10 is achieved = Sco was 18 milestone that important showed An date. to found color. ept h bevtoa ffr are u ntels e decade few last the in out carried effort observational the Despite ≈ ≈ .006:Sriy aiua19;Tay- 1994; Valiauga & Straizys 0.60-0.63: 5 .506:Ga 95 esl ta.1998) al. et Bessell 1995; Gray 0.65-0.68: ihSNKc IE pcr ocryotadffrnilaayi of analysis differential a out carry to spectra HIRES Keck S/N high ) hc erysassol esearched be should stars nearby which 1. A INTRODUCTION T E tl mltajv 6/22/04 v. emulateapj style X Sae 03.W r living are We 2003). (Seager tr:fnaetlprmtr tr:audne stars: – abundances stars: – parameters fundamental stars: pee u:fnaetlparameters fundamental Sun: – spheres D968 TRCOEYRSMLN U SUN OUR RESEMBLING CLOSELY STAR A 98618: HD og Mel Jorge endez ´ (B-V) edpof o .Melendez J. to: proofs Send 2 ai od-dn&Jos & Dodds-Eden Katie , alifornia color, mong ABSTRACT nol- an ´ ia evtr,Cte d,Wso re,AT21,Australia 2611, ACT Creek, Weston Rd., Cotter servatory, 05,crmshrcatvt,vsin v activity, Meléndez chromospheric & (Ram´ırez 2005), colors blue-optical-red-infrared very properties fundamental Sun. our with to both similar twins, as twin solar solar est a is Sco. 143436 candi- King, 18 HD twin recently, as that Very in four good suggested analog analyzed dwarfs. and solar (2005) G dates top Schuler 200 the & is about Boesgaard Sco of Tri- 18 sample & that Soubiran a (HD show Sco twin.” 18 solar (2004) that ever aud found closest found, (1997) “the solar Silva be is da to 146233) good yet & has “perfect Mello twin de a solar Porto perfect is a the analyzed Although of of stars none twin.” years that of 25 concluded number first and twins large the Strobel solar in de the for made Cayrel search the efforts Sun answering not. determine the or the reviewed in fairly anomalous (1996) is to of benefit Sun us the uniqueness allowing great whether 1998), the of Gustafsson of (e.g., be question will Sun. long-standing the star to similar host very importance stars ice greatest of “little the activity of a the is to monitor it related to Thus, ex- is 1976). of that (Eddy activity period age” a solar “Maunder the low 1976), for tremely Eddy case 1890; the (Maunder be may minimum” as changes, climate to errors by calibration. limited flux is absolute data its IRAC in of accuracy the however, prddHRSsetorp ihtehgetresolving (R highest power the the using with are spectrograph We HIRES Tele- upgraded I programs. Keck observing twin the other at solar during HIRES scope The using the observed 2005). being in is Fischer sample analysis & Valenti spectroscopic (e.g., previous literature from and allaxes - iue nodrt civ / fa least for at of observed S/N was a star achieve Each to order 2006. in minutes and 1-3 2005 in served chips. “red” and “green,” h addt oa wn eeslce ae ntheir on based selected were twins solar candidate The clos- the are Sco 18 and 98618 HD that show we Here our resembling closely stars of identification The related be to seem activity solar in variations Unusual ofr 6cniae lse nTbe1 aebe ob- been have 1) Table in (listed candidates 16 far, So ≈ 10 .Robles A. e ´ ,n efc oa wnhsbeen has twin solar perfect no s, 5 ys. udmna aaeeso both of parameters fundamental 2. ,covering ), ot eMlo&d iv,who Silva, da & Mello de Porto xrml ihrslto (R resolution high extremely e AAADANALYSIS AND DATA bezn.W ugs that suggest We zone. able tro u oa system, solar our of star t 6slrti candidates. twin solar 16 ciiy–sas atmo- stars: – activity 1 ≈ 0.4-0.8 µ i nte“blue,” the on m ipro par- Hipparcos , ≈ 300 2 Meléndez, Dodds-Eden, Robles per resolution element (or about 200 pixel−1) at 671 nm. during the same run and the data reduction and contin- For most stars, we fulfilled this requirement. During each uum normalization performed in the same way. On the run, a solar reflected spectrum was obtained by observing other hand, King et al. (2005) used a lunar spectrum either Ceres or Vesta. The data were reduced employ- obtained several years earlier, and their observations are ing the latest version of MAKEE (T. Barlow 2005, pri- of lower resolving power (R ≈ 45,000). The set of lines vate communication). Further data reductions (Doppler used is also different; while King et al. cover 647-675 nm, correction and continuum normalization) were performed our Fe I lines cover 522-675 nm. Note that part of the with IRAF. Sample spectra of the Sun, 18 Sco, HD 98618, region used by them has weak telluric contamination. and HD 143436 are shown in Figure 1. Orders falling on the “green” chip, which in our setup 3. DETAILED ANALYSIS OF 18 SCORPII AND HD 98618 covers ≈0.53-0.69 µm, were scrutinized for the best lines The top twin candidates, 18 Sco and HD 98618, were available for analysis. We selected lines of moderate analyzed in more detail to quantify their likeness to the strength (typically with line depths of 0.2-0.6, but weaker Sun by performing a differential spectroscopic equilib- and stronger lines were also included), unblended, rel- rium analysis with respect to the solar spectrum. The atively symmetric and free of telluric contamination. close proximity of our candidates to the Sun warrants Lines too close to the edges of the orders were discarded reliable relative chemical and physical parameters. because of the greater uncertainty in fitting the contin- The equivalent widths of 39 Fe I and 14 Fe II lines with uum in these regions. clean profiles and well-defined continuum regions were We employed two basic criteria to determine how close measured in the region 439-675 nm. The gf-values of a star is to the Sun. The first criterion assessed the simi- Fe I and Fe II lines were taken from the Oxford (e.g., larity between spectral lines in the sample star relative to Blackwell, Lynas-Gray & Smith 1995) and Hannover the Sun for different chemical elements. For each atomic groups (e.g., Bard & Kock 1994) and from Meléndez line, we determined the relative difference in line depth et al. (2006), respectively. Oxygen and lithium abun- δdr employing the 3 closest pixels to the line center. We dances were obtained from the O I triplet at 777 nm and critically selected 206 lines of 10 ions, with the largest the Li I profile at 670.7 nm. The interaction broadening I proportion (70) due to Fe . A median < δdr > was constants C6 were obtained from the collision broaden- determined for each ion and errors were estimated. A ing cross sections given by Barklem, Piskunov & O’Mara perfect solar twin should have < δdr > = 0 for lines of (2000) and Barklem & Aspelund-Johansson (2005). all the chemical elements, and although none was found, The LTE calculations were performed with the 2002 18 Sco and HD 98618 showed the < δdr > values closest version of MOOG (Sneden 1973), employing Kurucz to zero (i.e., solar) among the different chemical elements overshooting model atmospheres (Castelli et al. 1997). analyzed. The microturbulence vt was determined by requiring no The second criterion allowed us to identify how close a dependence of the iron abundance from Fe I lines (AF eI ) star is in temperature to the Sun. The principle behind against reduced equivalent width (Wλ/λ). this is the excitation equilibrium of Fe I lines, which are The zero point of the excitation and ionization very sensitive to effective temperature. For a star iden- equilibrium was defined by the analysis of the solar Fe I I tical to the Sun, the slope of a plot of δdr(Fe ) versus and Fe II lines: (dAF eI /dχexc)⊙ and (∆F eII−F eI )⊙ (= excitation potential χexc should be zero. In addition to AF eII - AF eI ). For the Sun we adopted the atmospheric I the Fe lines used in our first criterion, we selected four parameters (Teff , log g, [Fe/H] = 5777 K, 4.44, 0.0), lines around 525 nm from the “blue” chip, three of them which was also adopted as a first guess for the twin close to χexc = 0 eV and hence allowing a better deter- candidates. Then the parameters of the sample stars mination of the d(δdr)/d(χexc) slope. were iteratively changed in (∆Teff , ∆log g, ∆[Fe/H]) The results of the two criteria are shown in Table 1 and until we simultaneously achieved relative excitation and Figure 2, where it can be seen that 18 Sco is the best solar ionization equilibrium: twin out of our candidates, with HD 98618 also appearing very close to the Sun. Figure 2 also shows that, accord- (dAF eI /dχexc)∗ = (dAF eI /dχexc)⊙ ing to our criteria, HD 143436 is no match for 18 Sco, in contrast to the suggestion (King et al. 2005) that they (∆F eII−F eI )∗ = (∆F eII−F eI )⊙ may be equivalent. In addition, we estimated the relative depths of the Li feature at 6708 A.˚ HD 98618 and 18 Sco The stellar parameters so determined were used to ob- have similar Li profiles, which are both deeper than solar tain stellar masses and ages from Y2 isochrones (Demar- yet much closer to the Sun’s than the Li profile of HD que et al. 2004). Once the mass was estimated, a spec- 143436. Hence, also considering the Li abundance, HD troscopic luminosity was obtained from log (L/L⊙) = log 143436 is not a good solar twin. The spectroscopic log g (M/M⊙) - log(g/g⊙) + 4 log (Teff /Teff ⊙). (= 4.28) determined by King et al. (2005) is considerably We employed Hipparcos parallaxes to check both the lower than the trigonometric log g from Hipparcos par- surface gravities and luminosities obtained in the present allaxes (≈4.46, using the Teff and mass given by them), spectroscopic analysis. The trigonometric gravity (log which suggests errors in their analysis. Since the spec- gHip) and stellar luminosity (LHip) were calculated troscopic log g is too low, the spectroscopic Teff given by employing Hipparcos parallaxes, V, Mbol⊙ = 4.736 King et al. is probably also underestimated. (Bessell et al. 1998), and the bolometric correction It is important to note that our analysis is strictly differential and independent of model atmospheres, with BC = -1.6240 + 4.5066 θ -3.12936 θ2 (σ = 0.02 mag) both the sample and reflected solar spectrum (Vesta or Ceres) obtained employing the same HIRES setting and where θ = (5040 K)/Teff . BC was obtained from a HD 98618: A Star Similar to our Sun 3

fit to solar metallicity dwarfs from Alonso et al. (1995) 2%, Teff hotter by 0.7% and 1.1 %, metallicity higher and Blackwell & Lynas-Gray (1998); the relation above by 0.02 and 0.05 dex, and luminosity about 3% and 6 is valid from about 4200 to 7000 K. % higher, respectively. Both stars belong to the disk The spectroscopic and Hipparcos-based results are population, and, as the Sun, are orbiting the center of shown in Table 2, where one can see the excellent agree- our Galaxy at ≈8 kpc (Nordström et al. 2004). ment in both log g and luminosity, which lends confidence Both twin candidates seem to be rotating slightly faster to the results obtained in the present work. A compar- than the Sun and are about 10% younger than the Sun. ison with previous differential spectroscopic analyses is The higher Li abundance in 18 Sco and HD 98618 is per- shown in Table 3. All agree that the atmospheric pa- haps due to different amounts of pre-main-sequence de- rameters of 18 Sco and HD 98618 are very similar to the pletion (see, e.g., Ventura et al. 1998). Since 18 Sco and Sun, especially in Teff (1.1% or better). The differences HD 98618 are about 4 Gyr old, hypothetical terrestrial in log g and [Fe/H] are just within a few times 0.01 dex, planets may have had enough time to develop some kind almost the same as the typical uncertainties. Although of complex life, assuming the time-scale for formation of within the uncertainties all these works agree with each complex life is similar to Earth’s. Thus, we encourage other, the trigonometric gravities from Hipparcos paral- SETI programs to give top priority to these stars. laxes (Table 2) favor our results. Thanks to the long-term effort of the California and Rotation periods and stellar ages can be obtained from Carnegie Planet Search Project (see, e.g., Marcy et al. measurements of stellar activity. Table 2 shows the rel- 2005), we know that there are not any hot Jupiters ′ ative chromospheric activity indices (log RHK), chromo- around 18 Sco or HD 98618. According to Marcy et spheric ages and rotation periods estimated by Wright al. (2005) “among the nearest FGK main-sequence stars et al. (2004). A better rotation period for HD 146233 is (d < 40 pc) the yet undiscovered giant planets typically given by Frick et al. (2004), who analyzed time series of reside beyond 1 AU because giant planets within 1 AU chromospheric activity. have already been found.” Furthermore, Marcy et al. The rotation periods can be used to estimate a rota- estimate that undetected Jupiter-mass planets orbiting tional age (Barry 1988). The relation between age t and nearby FGK stars (<30 pc) may reside in orbits beyond rotation period P also depends on mass [log P = a log 3 AU or have masses less than MJup. Since the distance t + b log(M/M⊙) + c; Kawaler 1989), but since our at which an Earth-size planet would receive the same a stars have about 1M⊙, this relation simplifies to P ∝ t . amount of light from its star as Earth receives from the 1/2 Barry (1988) finds that a ≈ 1/e, and using a mean solar Sun is (L/L⊙) AU and the luminosities of the twin ◦ rotation period (equator to latitude 20 ) P⊙ ≈ 25.3 days candidates are very close to solar, the habitable zones (Howard 1984) and a solar main-sequence age of t⊙ = of 18 Sco and HD 98618 should be very close to 1 AU. 4.51 Gyr (Sackmann, Boothroyd & Kraemer 1993), we If further radial velocity observations show that the in- find t = (P /14.5)e, where P is in days and t is in Gyr. ner region around 18 Sco and HD 98618 is free of giant Using this relation, we estimated rotational ages (Table planets, then these stars have the potential to host terres- 2) adopting P from Frick et al. (2004) and Wright et al. trial planets around their habitable zones. We encourage (2004) for 18 Sco and HD 98618, respectively. long-term radial velocity monitoring of these stars, which 4. could reveal within a few years how far out gaseous giant DISCUSSION AND CONCLUSIONS planets, if any, reside within the exoplanetary systems As shown in Table 2, 18 Sco and HD 98618 are both these stars may host. very similar to the Sun, with the surface gravity being just 0.01 dex higher than solar, the mass above solar by

REFERENCES

Alonso, A., Arribas, S. & Mart´ınez-Roger, C. 1995, A&A, 297, 197 King, J. R., Boesgaard, A. M. & Schuler, S. C. 2005, AJ, 130, 2318 Bard, A. & Kock, M. 1994, A&A, 282, 1014 Luck, R. E. & Heiter, U. 2005, AJ, 129, 1063 Barry, D. C. 1988, ApJ, 334, 436 Marcy G. W. et al. 2005, ApJ, 619, 570 Barklem, P. S., Piskunov, N. & O’Mara, B.J. 2000, A&AS, 142, Maunder, E. W. 1890, MNRAS, 50, 251 467 Meléndez, J., Shchukina, N. G., Vasiljeva, I. E. & Ram´ırez, I. 2006, Barklem, P. S. & Aspelund-Johansson, J. 2005, A&A, 435, 373 ApJ, in press (astro-ph/0601256) Bessell, M. S., Castelli, F. & Plez, B. 1998, A&A, 333, 231 Nordström, B. et al. 2004, A&A, 418, 989 Blackwell, D. E.& Lynas-Gray, A. E. 1998, A&AS, 129, 505 Porto de Mello, G. F. & da Silva, L. 1997, ApJ, 482, L89 Blackwell, D. E., Lynas-Gray, A. E. & Smith, G. 1995, A&A, 296, Ram´ırez, I. & Meléndez, J. 2005, ApJ, 626, 465 217 Reach, W. T. et al. 2005, PASP, 117, 978 Castelli, F., Gratton, R. G. & Kurucz, R. L. 1997, A&A, 318, 841 (erratum 324, 432) Sackmann, I.-J., Boothroyd, A. I. & Kraemer, K. E. 1993, ApJ, Cayrel de Strobel, G. 1996, A&A Rev., 7, 243 418, 457 Colina, L., Bohlin, R. C., Castelli, F. 1996, AJ, 112, 307 Seager, S. 2003, Earth and Planetary Science Letters, 208, 113 Demarque, P., Woo, J.-H., Kim, Y.-C., Yi, S. K. 2004, ApJS, 155, Sekiguchi, M. & Fukugita, M. 2000, AJ, 120, 1072 667 Sneden, C. 1973, Ph.D. thesis, University of Texas Eddy, J. A. 1976, Science, 192, 1189 Soubiran, C. & Triaud, A. 2004, A&A, 418, 1089 Frick, P., Soon, W., Popova, E., Baliunas, S. 2004, NAR, 9, 599 Straizys, V. & Valiauga, G. 1994, Baltic Astronomy, 3, 282 Gray, D. F. 1995, PASP, 107, 120 Taylor, B. J. 1994, PASP, 106, 444 Gustafsson, B. 1998, Space Sci. Rev., 85, 419 Valenti, J. A. & Fischer, D. A. 2005, ApJS, 159, 141 Holmberg, J., Flynn, C., Portinari, L. 2006, MNRAS, in press Ventura, P., Zeppieri, A., Mazzitelli, I., D’Antona, F. 1998, A&A, (astro-ph/0511158) 331, 1011 Howard, R. 1984, ARA&A, 22, 131 Wright, J. T., Marcy, G. W., Butler, R. P., Vogt, S. S. 2004, ApJS, Kawaler, S. D. 1989, ApJ, 343, L65 152, 261 Note added in proof.- After this Letter was accepted for publication, M. C. Turnbull brought it to our attention 4 Meléndez, Dodds-Eden, Robles

TABLE 1 Slope d(δdr )/d(χexc) for Fe I, and Relative Differences in Line Depth < δdr > for Fe I, V I, and Ca I

Star Slope Fe I V I Ca I HD 9986 −0.0120 −0.042 −0.100 −0.019 HD 32963 −0.0075 −0.112 −0.301 −0.080 HD 33636 0.0048 0.220 0.445 0.116 HD 45184 −0.0063 0.010 0.098 0.004 HD 56124 0.0010 0.022 0.145 0.010 HD 71334 −0.0120 0.001 0.079 −0.027 HD 71148 −0.0047 0.037 0.145 0.021 HD 71881 0.0120 0.077 0.203 0.045 HD 98618 −0.0030 −0.003 −0.016 0.002 HD 112257 0.0046 −0.047 −0.292 −0.070 HD 143436 −0.0064 0.008 0.078 0.011 HD 146233 −0.0007 −0.001 0.024 −0.004 HD 147044 0.0040 0.090 0.270 0.055 HD 159222 −0.0160 −0.017 0.000 0.013 HD 186427 −0.0062 −0.079 −0.213 −0.053 HD 187123 −0.0120 −0.081 −0.153 −0.052

TABLE 2 Fundamental Parameters

Parameter (Star – Sun) 18 Sco HD 98618 − ∆ vt (km s 1) +0.08±0.15 +0.09±0.15 ∆ Teﬀ (K) +40 ± 30 +66 ± 30 ∆ log gspec (dex) +0.01 ± 0.04 +0.01± 0.04 ∆ log gHip (dex) +0.01 ± 0.02 +0.01± 0.03 ∆ log gadopted (dex) +0.01 ± 0.02 +0.01 ± 0.03 ∆ Lspec (L⊙) +0.02 ± 0.06 +0.04± 0.06 ∆ LHip (L⊙) +0.03 ± 0.03 +0.08± 0.07 ∆ Ladopted (L⊙) +0.03 ± 0.02 +0.06± 0.05 [Fe/H] (dex) +0.02 ± 0.03 +0.05± 0.03 [O/H] (dex) −0.03 ± 0.05 0.00 ± 0.04 [Li/H] (dex) +0.53 ± 0.09 +0.47 ± 0.09 ∆ mass (M⊙) +0.02 ± 0.03 +0.02 ± 0.03 ∆ ageisochro (Gyr) −0.8 ± 1.5 −1.1 ± 1.5 1 1 ∆ agechromos (Gyr) −0.3 +0.7 ∆ agerotation (Gyr) −1.1 −0.4 ∆ ageadopted (Gyr) −0.7 ± 0.4 −0.3 ± 0.9 ∆ rotation period (days) −2.52, -11 −11 ′ 1 1 ∆ logRHK (dex) 0.0 −0.05 ∆MV (mag) −0.04 ± 0.04 −0.09 ± 0.07 B-V 0.65 0.64 Distance(pc) 14.0 38.7

References. — (1) Wright et al. 2004; (2) Frick et al. 2004 that HD 98618 is included in the Catalog of Nearby Habitable Stellar Systems (M. C. Turnbull & J. C. Tarter, ApJS, 145, 181; ApJS, 149, 423 [2003]) for SETI searches with the Allen Telescope Array. It is ranked No. 106 in terms of its location on the H-R diagram, and considering its distance merits its ﬁnal ranking is No. 504. The present work suggests that HD 98618 should be given a higher ranking. HD 98618: A Star Similar to our Sun 5

TABLE 3 Comparison of Relative Spectroscopic Parameters ∆(Teﬀ /log g/[Fe/H])

18 Sco HD 98618 Ref. (K/dex/dex) (K/dex/dex) +40/+0.01/+0.02 +66/+0.01/+0.05 1 +12/+0.05/+0.05 2 +58/+0.03/+0.03 3 +14/−0.03/+0.03 +35/−0.02/+0.03 4

References. — (1) This work; (2) Porto de Mello & da Silva (1997); (3) Luck & Heiter (2005); (4) Valenti & Fischer (2005)

Fig. 1.— Left, sample spectra of the Sun and the candidate solar twins 18 Sco, HD 98618, and HD 143436; right, comparison of the solar-twin candidates (solid lines) with the Sun (dotted line) around 6090 A.˚ Since V I is very sensitive to Teﬀ , the weaker V I line in HD 143436 suggests that this star has a Teﬀ higher than solar, while both 18 Sco and HD 98618 are very similar to the Sun.

0.02 0.005

0.015

0.01 exc

χ 0.005 vs r d

δ 0 0 HD 146233 HD 98618 −0.005 HD146233

Slope of HD 186427 HD 143436 (smaller circle) −0.01 HD98618 −0.015

−0.005 −0.02 −0.25 −0.2 −0.15 −0.1 −0.05 0 0.05 0.1 0.15 0.2 −0.005 0 0.005 <δd > <δ d > r r HD45184

Fig. 2.— Slope d(δdr )/d(χexc) vs. the median relative diﬀerence in line depth < δdr > for Fe I lines (§2). A perfect solar twin should be at (0,0). The size of the circles is proportional to the Li I line depth at 6708 A,˚ with the smallest circles representing closer agreement to the solar Li I proﬁle. Left, the whole sample, with the inner rectangle representing the typical ±1σ error bars; right, a zoom to the inner rectangle of the left panel. The two best solar twin candidates, 18 Sco (HD 146233) and HD 98618, are labeled; HD 143436 (King et al. 2005) and HD 186427 (one of the most promising twin candidates in the early literature) are shown as well.