arXiv:2101.02242v1 [astro-ph.EP] 6 Jan 2021 Swaelmen ob bandwl eoeteluc.I hswr,we between work, analysis this bench-marking In a launch. of results the the necessary before present is well this obtained and be of properties to knowledge star accurate hosting meaningful an planet strategy A requires the list. selection targets target the the input of is the choice of Ariel definition of the goal for of scientific achievement the population the for diversified of factor key and A in the exoplanets. large observe launched transiting a will be Ariel of to mission, atmospheres 4-year expected its is During it 2028. and mission science Abstract date Accepted: / date Received: Sanna Nicoletta ail Danielski Camilla naBrucalassi a metho Anna targets: spectroscopic Ariel different between for analysis parameters comparison stellar of Determination ta ulig em vne awl apslDdo l Didcot l Campus (UK) Harwell 0QR Avenue, OX11 Fermi Data Building, Exochemistry Atlas Space for Centre UCL Portugal Danielski Porto, C. 4150-762 Estrelas, das Rua CAUP, Espa¸co, Universi Porto, Ciˆenciasdo do Astrof´ısica e Delgado-Mena de E. Instituto Adibekyan, V. Sousa, S. Italy Catone, Porzio Monte 00044 33, Frascati di Via Roma di Astronomico INAF-Osservatorio Biazzo K. Italy (Firenze), Fiorentino Sesto 50019 1, Sansone G. via Stu Universit`a degli Astronomia, Firenze, di e Fisica di Dipartimento Casali G. Italy Firenze, 50125 5, Fermi Enrico Swaelmen, Largo Arcetri, der di Van Astronomico M. INAF-Osservatorio Sanna Casali, N. G. Rainer, Magrini, M. L. Tsantaki, M. [email protected] E-mail: +39-055-220039 Fax: 243 Firenze +39-055-2752 50125 Tel.: 5, Fermi Enrico Largo Arcetri, Astronomico INAF-Osservatorio Brucalassi A. wl eisre yteeditor) the by inserted be No. (will manuscript Astronomy Experimental re a enslce stenx S M4 ESA next the as selected been has Ariel · oiaRainer Monica · · ai Tsantaki Maria ai Biazzo Katia · adnAdibekyan Vardan · · id Casali Giada ar Magrini Laura dade di · · rthnso h nenlsrcueo h exoplanet the of structure the internal have we the density of bulk hints the From first obtained. planet is the and measurements, measured mass be Doppler can spectroscopic radius from planet tran- the the From light-curve, exoplanets. sit of structure understand- internal on the focused ing been have studies these first, At ob- detected the jects. of history evolutionary and formation, of era an isation into research exoplanet ushered searches velocity has radial and [32], of transit space success ESPRESSO ground-based, recent The and SPIRoU[6]. [29] and CARMENES[36] HARPS [24], [66], NGTS TRAPPIST and [35], WASP [11] like along surveys, [37], CHEOPS ground-based PLATO with as few [38], such surveys, next TESS space the upcoming [51,33], the in Gaia discovered by be years to expected are Earth-size to planets down orbits. Jupiter-size and of an sizes thousands masses, showing Moreover, of 3000 detected terms than in been diversity more have incredible (with stars) planets their 4000 transiting than more far, So Introduction 1 Spectroscopy IR and Optical Keywords properties. for planet factor the key a of is accuracy turn the in and character- stars for exoplanet-host crucial ising is precise parameters and stellar of accurate derivation Ariel Homogeneous, complete the homogeneously sample. of to reference parameters used stellar be the will determine that con- to method aim We a sample. solidate reference Ariel targets the of number to selected belonging deter- a for to parameters used stellar mine techniques spectroscopic different three ahe a der Van Mathieu ls Delgado-Mena Elisa · egoSousa Sergio tde ihtega oivsiaetenature, the investigate to goal the with studies xpae atmospheres Exoplanet · · · pc missions Space character- ds. · 2 Anna Brucalassi et al. and the gas/ice/rock ratios. However, to have a reli- allowing to measure planetary atmospheric signals of able estimate of the density, an accurate knowledge of 10-100 ppm relative to the star. Such small signals re- both radius (with a precision of up to 5%) and plane- quire an exact knowledge of the host star spectrum, tary mass (up to 10%) is necessary [65,8]. Additionally, at least at the same level of the planetary signal, to the absolute value of planetary radius and mass relies map any stellar intrinsic variation (i.e. due to magnetic on the precise determination of the radius and mass activity and convective turbulence) in order to avoid of the exoplanet-host star. The derivation of these last misleading results with planetary features [39]. two values is in turn strongly connected to the effec- Another important point is that information on the tive temperature (Teff ), surface gravity (log g), and the host star composition is critical to separate the sig- metallicity of the star. Thus, the planetary properties natures left on the planet by its formation, evolution are critically dependent on their stellar host properties and migration processes, from those due to the spe- [57,50,2]. cific chemistry of the host star [63]. Indeed, recent stud- Furthermore, transiting planets provide us one of ies suggest that planetary O/H, C/H, C/O ratios and the best ways of characterising their atmospheres. In- metallicity with respect to the stellar values could pro- transit spectroscopy as well as secondary transit stud- vide stronger constraints on the planet formation region ies [12,54,44,20,62,68,34] using space observatories, as and their migration mechanisms (see [27] and references Hubble and Spitzer, and some ground-based observa- therein for a recent review), but similar considerations tories, have yielded the detection of some important apply to other elements (e.g., N, S, Ti, Al, [63]). molecules present in the planetary atmospheres for a Finally, an increasing number of studies have pointed limited number of targets, or have identified the pres- towards the existence of correlations between the prop- ence of clouds, probing the thermal structure and pro- erties of the host stars and the characteristics and fre- viding some constraints on the planet properties. How- quency of their planetary systems. In this respect, the ever, the data available is still too sparse to provide a correlation between the stellar metallicity and the fre- consistent interpretation and the achieved results point quency of giant planets [41,64,49], the connection be- out the main limitations of the existing facilities: very tween radius vs metallicity [13,43], eccentricity vs metal- narrow wavelength coverage, observations usually not licity [3,67], the role of the abundances of other ele- simultaneous for a wider spectral range with the intro- ments [4,16,17] in the host stars are only few examples duction of systematic noise, insufficient time allocated of different results that take a clear shape as the new to exoplanet science, and more in general the lack of planet discoveries increase, shading light on many de- a dedicated space-based exoplanet spectroscopy mis- tails still missing concerning planet formation and evo- sion. Thus, our current knowledge of exoplanetary at- lution. mospheric and thermal characteristics is still very lim- Such works rely upon homogeneously and precisely de- ited. rived stellar parameters. Therefore, homogeneous deriva- Ariel (Atmospheric Remote-sensing Infrared Exo- tion of stellar parameters using high-quality data is planet Large-survey) has been selected as the next ESA- crucial for characterising exoplanet-host stars [42,5,47], M4 science mission [55] and it is expected to be launched and in turn, is fundamental to improve the accuracy of in 2028. During its 4-year mission, Ariel will observe the planet properties. the atmospheres of a statistically representative sam- It is important to underline how a well-defined target ple (∼1000) of transiting gaseous (Jupiters, Saturns, selection strategy and the definition of the input target Neptunes) and rocky (super-Earths and sub-Neptunes) list with a statistically significant dataset in the range of planets using transit spectroscopy in the 1.10-7.8µm relevant stellar/planet parameters have a fundamental spectral range and three narrow-bands photometry in role for maximising the scientific yield and the achieve- the optical. The wavelength range proposed covers all ment of scientific goals of Ariel. A meaningful choice the expected major atmospheric gases from H2O, CO2, of the targets requires an accurate study of the stel- CH4, NH3, HCN, H2S up to the more exotic metallic lar properties that need to be derived in advance and compounds, such as TiO, VO, and condensed species. continuously updated as the mission approaches launch Ariel is designed as a dedicated survey mission for tran- and the target list evolves with the new exoplanet dis- sit, eclipse and phase-curved spectroscopy, providing a coveries. homogeneous dataset, with a consistent pipeline and In this context, we have started a benchmarking anal- an well-defined target selection strategy, maximising ysis between three different spectroscopic techniques the scientific yield. Focusing on transit, eclipse spec- used to determine stellar parameters for selected tar- troscopy, the methods are based on the differential anal- gets belonging to the Ariel Reference Sample [19]. Our ysis of the star and planet spectra in and out of transit, goal is to consolidate a method that will be applied to Title Suppressed Due to Excessive Length 3 homogeneously determine the stellar parameters for the complete Ariel Reference Sample. More generally, we refer also to the works of [26] and [45] where compar- 50 isons between the results from different spectroscopic 40 techniques are discussed.
For a global approach to characterize the stars of 30 Ariel Reference Sample see also [15], where an overview on the methods used to determine stellar fundamental 20 Star Sample Star parameters, elemental abundances, activity indices, and stellar ages for the Ariel Reference Sample is given and 10 in particular, results for the homogeneous estimation 0 of elemental abundances of Al, Mg, Si, C, N, and the 4 6 8 10 12 14 16 activity indices S and log(R’ HK) are presented. SweetCat Vmag In the following sections we describe the star sam- ple analysed, we give an overview of the methodology applied to derive the stellar parameters and then we 6 present a comparative analysis of the results. 8
10 2 Star Sample
12
On the bases of Ariel capabilities, a list of targets (Ariel SweetCatVmag Reference Sample) to be observed during the primary 14 mission life was prepared for the Phase A [69,19]. The sample includes ∼1000 potential targets with stellar 0 200 400 600 800 1000 types FGKM (typically brighter than K = 11 mag) and S/N planetary parameters in a range of size between Jupiter Fig. 1 Top panel: Distribution of the star sample analyzed down to Earth-like planets, temperature in the range as a function of V mag listed in the SWEET-Cat catalog. between 500K-2500K and bulk density between 0.10- Bottom panel: V mag vs S/N for our star sample. 10.0g/cm3. We started our analysis cross-matching the Ariel Table 1 Spectrograph information: resolving power, spectral ranges and number of stars considered. Reference Sample with the stars available in the SWEET- Cat (Stars With ExoplanETs Catalogue 1) and we se- Instrument Resolving Power Spectral Range N lected all the common targets with parameters source (λ/∆λ) (A˚) of Stars ”flag=1”, which are the stars analysed homogeneously ESPADONS 80000 3700-10500 2 FEROS 48000 3600-9200 38 (see [42]). Thus, our starting sub-sample includes 155 FIES 67000 3700-7300 7 FGK stars in a range of 5 3 Derivation of Stellar Parameters rors due to the uncertainties in the stellar parameters. The convergence criteria are not fixed a priori but are 3.1 SWEET-Cat based on the quality of the spectra. The code is de- scribed in [28]. For our work, first, we have measured SWEET-Cat is a catalog of stellar parameters taken in the EWs with DOOp (DAOSPEC Option Optimiser, general from the literature for planet hosting stars listed [14]) an automatic tool developed within the Gaia-ESO 2 within the Extrasolar Planets Encyclopaedia . This survey. The code is based on DAOSPEC code [53] and catalogue is continuously updated when new planets are uses Gaussian fit to measure EWs. We adopted the line announced and new stellar parameters derived. In par- list of the Gaia-ESO survey [23], which includes atomic ticular, for stars with spectra acquired with high reso- parameters (log g and damping coefficient) for a large lution spectrographs (mainly HARPS, UVES, FEROS) number of lines in the spectral range 4200-6800 A.˚ We and high signal-to-noise ratios (mostly S/N>100), the adopted the MARCS model atmospheres (plane paral- stellar parameters are obtained in a homogeneous way lel and spherical) [22]. (flagging the targets with 1, as mentioned above). The method used to derive stellar parameters is described in, e.g., [41] and [49]. Briefly, local thermodynamic equi- 3.3 FASMA analysis librium (LTE) condition is assumed, a grid of Kurucz plane-parallel model atmospheres ([25], ATLAS9), and The analysis with FASMA is based on the spectral syn- the ARES code [48] for measuring line equivalent widths thesis technique wrapped around the radiative transfer (EWs) are considered. Stellar parameters (effective tem- code, MOOG (version 2019; [46]). perature Teff , surface gravity log g, microturbulence ve- FASMA creates synthetic spectra on-the-fly and deliv- locity Vmicro, and iron abundance [Fe/H]) are derived ers the best-fit parameters (effective temperature, sur- using the MOOG code (version 2002; [46]) and the line face gravity, iron abundance, and projected rotational list by [49]. Excitation and ionization equilibria of the velocity) after a non-linear least-squares fit (Levenberg- Fe I and Fe II weak lines were imposed to derive stellar Marquardt algorithm). The line list is mainly comprised parameters. The errors on the atmospheric parameters of iron lines initially obtained from VALD3 [40]. The re- were derived as in [41] and references therein. gions of the spectral synthesis are defined within small The stellar masses listed in the SWEET-Cat cata- intervals around these iron lines (±2A).˚ The atomic log were derived with the calibration presented in [56], data are calibrated to match the spectra of the Sun using as input the spectroscopic parameters. According and Arcturus but given the strongest constrains on the to [42], a correction was applied for the cases in which solar parameters, we gave higher weights to the Sun. the calibration gives values between 0.7 and 1.3M⊙.The The damping parameters are based on the ABO the- errors for these mass values are computed as in [42] by ory [7] when available, or in any other case, we use the means of a Monte Carlo analysis. For each case 10 000 Blackwell approximation. The model atmospheres are random values of effective temperature, surface grav- interpolated from the ATLAS grid [25] in LTE. The ity, and stellar metallicity were drawn from a Gaussian minimization is a two step procedure where initially we distribution. From the resulting mass distribution, the obtain the best-fit values using fixed solar macro- and central value for the mass and 1-sigma uncertainty were microturbulence. Then, we refine our results in a sec- derived. ond step, with updated macro- and micro-turbulence based on empirical relations. Macroturbulence velocity is set by the calibration of [18] and microturbulence is 3.2 FAMA analysis set based on calibrations for either dwarfs [60] or giants [1]. The uncertainties are derived from the covariance The aim of FAMA is to allow the computation of the at- matrix constructed by the nonlinear least-squares fit. mospheric parameters and abundances of a large num- The methodology is described in detail in [59,58] where ber of stars using measurements of equivalent widths it is tested for both giant and dwarf samples, including (EWs) as automatic and as independent of any sub- stars with larger rotation. jective approach as possible. It is based on the simul- taneous search for three equilibria: excitation equilib- rium, ionization balance, and the relationship between 4 Comparison among the methods log n(FeI) and EW/λ. FAMA also evaluates the statis- tical errors on individual element abundances and er- We compare our results seeking possible outliers and 2 http://exoplanets.eu/ trends among the parameters obtained by the three Title Suppressed Due to Excessive Length 5 Table 2 Targets. Spectra source. Stellar Parameters (Teff , log g and [Fe/H]) and errors obtained by FAMA and FASMA for each star of the analysed sample. FAMA FASMA Star Spectra Teff eTeff log g elog g [Fe/H] e[Fe/H] Teff eTeff log g elog g [Fe/H] e[Fe/H] (K) (K) (dex) (dex) (dex) (dex) (K) (K) (dex) (dex) (dex) (dex) 55Cnc UVES 5358 150 4.44 0.20 0.49 0.01 5295 20 4.455 0.03 0.359 0.01 CoRoT-10 HARPS 4908 150 4.31 0.20 0.26 0.01 4972 10 4.439 0.03 0.262 0.01 CoRoT-19 HARPS 6331 150 4.07 0.20 0.15 0.06 6313 40 4.283 0.06 0.148 0.02 methods. In this preliminary work, we explore the re- SWEET-cat are increasing with log g. The differences sults of the three analysis, comparing with external are particularly high for higher surface gravity (log g benchmarks, as the surface gravities derived from light- >4.6dex). In the log g range of ∼ 4 − 4.6 dex, the trend curves in literature and from Gaia parallaxes and pho- is almost negligible. In addition, FAMA, on average, tometry, to identify the ranges of best performances. gives lower gravities, of about ∼−0.18 dex. 4.1 SWEET-Cat vs FAMA/FASMA: Teff 4.3 SWEET-Cat vs FAMA/FASMA: [Fe/H] Fig. 2 (top panels) shows the difference in effective tem- perature (Teff ) between SWEET-Cat and FAMA (left Fig. 2 (bottom panels) show respectively the compar- panel) or SWEET-Cat and FASMA (right panel) as a ison of [Fe/H] derived using the FAMA (left panels) function of the SWEET-Cat Teff . The color scale rep- method and FASMA (right panels) vs the results given resents log g for SWEET-Cat. in SWEET-Cat. The colour scale represents Teff for The bulk of the Teff results for the SWEET-Cat SWEET-Cat. and FAMA is in good agreement (within 1-σ), although Considering the typical errors on metallicity (∼0.10 there are some outliers: 28 stars show a difference in Teff dex), the results are quite satisfactory: there are almost ∼ over 1σ level but still within 417K (3σ level). Among no trend in the differences vs [Fe/H] and the means are these outliers, we noticed that 7 stars present large vsini close to zero (slightly negative for FAMA). (>10km/s). The Teff spread seems more evident when This is one of the most important results of the anal- considering stars with the lower and higher Teff in the ysis and we need to remark it. Despite the differences considered range. In addition, there is an offset between in the derived stellar parameters, the three methods the two methods, with the mean value (red dashed cen- converge to very similar final metallicities [Fe/H]. The tral line) for the difference of ∼70 K. three methods rely, indeed, on different line lists, with Focusing now on the FASMA analysis, also in this case, different set of atomic data, which make them to con- the bulk of the Teff results for the SWEET-Cat and verge on slightly different sets of Teff and log g val- FASMA is in good agreement (within 1σ), although ues, which should satisfy the conditions of excitation outliers are also present. For stars with the lower and equilibrium and ionization balance. Teff and log g tend higher Teff in the range of interest, the dispersion be- to vary along a local minimum, in which the value of tween the two methods increases. An offset between [Fe/H] is usually more constant. the two methods is present with a mean value (red In Fig. 3, we show the differences in metallicity of dashed line) in the difference of ∼34K. We note here the FAMA and FASMA results with respect to SWEET- that some of the outliers are not in common with the Cat, as a function of the differences in the derived sur- outliers found using the FAMA method. face gravities. The symbol are colour-codes by the dif- ferences in the derived Teff , again, with respect to the 4.2 SWEET-Cat vs FAMA/FASMA: log g SWEET-Cat ones. Removing the systematic offsets (+0.05 dex from Fig. 2 (central panels) show respectively the comparison FAMA and -0.05 dex for FASMA with respect to the of log g derived using the FAMA (left panel) or FASMA SWEET-Cat scale), the differences in [Fe/H] are negli- (right panel) methods vs the results given in SWEET- gible for variation in the surface gravity within ±0.3 dex Cat. The colour scale represents Teff for SWEET-Cat. with respect to the SWEET-Cat ones. In the compar- There are trends in both the FAMA and FASMA re- ison between FAMA and SWEET-cat [Fe/H] even for sults: the differences between their log g and those of differences in log g∼ -0.3 dex, the agreement in metallic- 6 Anna Brucalassi et al. 5.0 5.0 4 4