arXiv:1804.09082v3 [astro-ph.SR] 31 May 2019 upr o h oeaceinmcaimo lntformation. planet of mechanism accretion core the for support httepoaiiyo omto fagatpae sapoiaeyao a disk approximately protoplanetary is the planet in metals giant of a amount of total formation of probability the that in lntocrec nrae ihbt tla asadmetallic and mass stellar both with o increases sample updated occurrence an planet with giant together stars, a evolved previous of in sample overestimated our been for not have Stars”) A “Retired (the oainlvlcte rvd diinlcnrainta h mass the that confirmation additional Sim provide literature. velocities the rotational in claims to contrary distributions, velocity space o h rttm.We ednn orcin r osdrdi th in considered are -0.12 corrections a reddening and When s time. this first of the and velocities for subgiants space the 245 calculate of also We sample co parameters. an these physical on address analysis to spectroscopic order a In ve perform occurrence. the planet regarding and concerns mass recent stellar to th challenging led more has far which is counterparts, giants and subgiants for masses accurate lnt eed o nyo tla ealct [eH) u lotema the also but ([Fe/H]), metallicity stellar on only not depends planets yee sn L 2019 using 4, Typeset June version Draft 3 EIE TR EIIE:A PAE IN LNTOCCURREN PLANET GIANT UPDATED AN REVISITED: STARS A RETIRED 2 1 [email protected] avr-mtsna etrfrAtohsc,6 adnSt Garden 60 Astrophysics, for Center Harvard-Smithsonian eateto srnm n srpyis nvriyo Ch of University Astrophysics, and Astronomy of Department bevtoi ainl u eea o´ rsio 7 20 77, Jos´e Cristino, General Rua Observat´orio Nacional, xpae uvy feovdsashv rvddicesn evide increasing provided have stars evolved of surveys Exoplanet M A T ⊙ unGhezzi, Luan E X ffe sapidt h eut,temse ftesbinsaecon are subgiants the of masses the results, the to applied is offset manuscript UCINO TLA EALCT N MASS AND METALLICITY STELLAR OF FUNCTION A 1 tl nAASTeX61 in style ejmnT Montet, T. Benjamin ABSTRACT 2-0,Sa rs´va,Rod aer,R,Brazil; RJ, S˜ao Janeiro, Crist´ov˜ao,921-400, de Rio M et abig,M 23 USA 02138 MA Cambridge, reet, 2, ⋆ cg,54 .ElsAe hcg,I 03,USA 60637, IL Chicago, Ave, Ellis S. 5640 icago, 10 ∗ n onAhrJohnson Asher John and [ e/H F ] hscreainpoie additional provides correlation This . so ugat with subgiants of es aiyo h orlto between correlation the of racity crsw s IE pcr to spectra HIRES use we ncerns GMdaf,w ofimthat confirm we dwarfs, FGKM f ml nahmgnosmanner homogeneous a in ample aye.Uigteenwresults new these Using nalyses. ni sfrtermain-sequence their for is it an c httefraino giant of formation the that nce lry u esrmn ftheir of measurement our ilarly, acltoso tla masses stellar of calculations e t pt . M 2.0 to up ity s( ss eienwamshrcand atmospheric new derive M et-n ucino the of function ne-to-one ⋆ .Hwvr measuring However, ). itn ihtheir with sistent 3 ERT AS RATE CE M ⊙ eshow We . ⋆ ≥ 1 . M 6 ⊙ 2 Ghezzi, Montet & Johnson
Keywords: stars: fundamental parameters — stars: atmospheres — stars: evolution — planets: formation
1. INTRODUCTION
The ever increasing sample of exoplanets has facilitated many robust studies of how planet formation and evolution are affected by the physical properties of their host stars. One such result is the correlation between stellar metallicity and the occurrence rate of giant planets around FGK dwarf and subgiant stars (e.g., Gonzalez 1997; Fischer & Valenti 2005; Ghezzi et al. 2010a; Wang & Fischer 2015). Similar analyses for the stellar masses, however, are much more scarce and controversial.
The main limitation appears on the massive end (&1.2 M⊙) of the planet search samples, since the detection or confirmation of exoplanets around main-sequence A and early F stars using radial velocities (RVs) is more difficult due to their higher effective temperatures and larger rotational velocities. In order to overcome this issue, radial velocity surveys have focused on intermediate-mass stars. As stars evolve off of the main sequence towards the red giant branch (RGB), their effective temperatures decrease and their rotation velocities slow down, resulting in spectra from which precise Doppler shifts can be measured. To date, more than 100 exoplanets have been detected around more than 1000 evolved stars (Jofr´eet al. 2015; Niedzielski et al. 2015). In one of these efforts, Johnson et al.
(2010a) analyzed a sample of 246 subgiants from the SPOCS IV catalog, as well as 948 FGKM dwarfs and subgiants from the Keck M Dwarf Survey (Butler et al. 2006) and the original SPOCS catalog (Valenti & Fischer 2005). Using the complete sample of 1194 stars, with masses ranging from 0.2 M⊙to 2.0 M⊙, Johnson et al. (2010a) determined that the occurrence rate of giant planets increases approximately linearly with stellar mass, from ∼3% for M dwarfs, to ∼8% for FGK dwarfs,
∗ Sagan Fellow Retired A Stars Revisited 3
to ∼14% for stars more massive than 1.5 M⊙, the so-called ”Retired A stars.” This correlation appears to be supported, at least for the FGKM dwarfs, by results from transit surveys (Fressin et al. 2013; Gaidos et al. 2013). Johnson et al. (2010a) also confirmed that the planet-metallicity correlation holds for both M dwarfs and the Retired A Stars. This suggests that planet formation is significantly affected by the stellar mass and chemical composition—two stellar properties that provide crucial links to the physical properties of stars’ protoplanetary disks. These findings were questioned by Lloyd (2011), who pointed out that some discrepancies between the rotational velocities distributions of evolved planet hosts and field stars could be traced back to erroneous mass determinations. Moreover, Lloyd found that the numbers of massive subgiants (M⋆&
1.5 M⊙) observed by Johnson et al. (2010a) seems inconsistent with the predictions from Galactic stellar population models. More massive stars are formed less frequently according to the initial mass function (IMF) and also have a relatively rapid evolution across the Hertzprung Gap. Based on these arguments, Lloyd concludes that many of the Retired A Stars should in fact lie in the mass range 1.0–1.2 M⊙. The errors on the original masses could stem from systematic uncertainties on spectroscopically determined atmospheric parameters or ambiguities in the predictions from stellar evolution models. Johnson et al. (2013) addressed the second concern by generating Galactic stellar population models with TRILEGAL (Girardi et al. 2005). Using a magnitude-limited sample (instead of the volume- limited sample used by Lloyd 2011), they produced a simulated distribution that allows for a larger number of massive stars than Lloyd found—an effect similar to the Malmquist bias—and is thus consistent with the masses determined by Johnson et al. (2010a). Lloyd (2013), on the other hand, argues that the mass distribution is not sensitive to the Malmquist bias, but instead to the usage of different Galactic models or input parameters within a given model. Using the Besan¸con model
(Robin et al. 2003), Lloyd (2013) obtains a mass distribution with less massive subgiants that can not be reconciled with the observed data from Johnson et al. (2010a). The fraction of massive subgiants within the Retired A Star sample was also questioned by Schlaufman & Winn (2013). Based on an analysis that is independent of stellar evolution mod- 4 Ghezzi, Montet & Johnson els, they showed that the space velocity dispersion of the subgiants with planets were larger than those of a sample of main-sequence A5–F0 stars, but consistent with the that obtained for F5–G5 dwarfs. They thus conclude that the subgiants with planets are most likely the evolved counterparts of less massive main-sequence solar-type stars, instead of “Retired A Stars”. Their analysis also revealed that no more than 40% of the planet-hosting subgiants could have been A5–F0 stars more massive than ∼1.3 M⊙ while on the main sequence. In another study, Sousa et al. (2015b) pointed out that, while the masses they determined for planet-host stars were in general good agreement with values from the literature, there were some notable exceptions. In particular, differences of up to 100% were found for a few evolved stars. Note, however, that some of these inconsistencies are explained by erroneous mass determinations in the original papers, as explained by (Takeda & Tajitsu 2015, see Section 4.2).
More recent studies of the stellar mass dependence of planet occurrence have been performed by Reffert et al. (2015) and Jones et al. (2016). In the former study, the analysis of 373 G and K giants reveals that giant planet occurrence increases with stellar mass in the range 1.0 – 1.9 M⊙, consistent with the findings of Johnson et al. (2010a). Moreover, there seems to be a maximum for
+0.1 the occurrence rate as a function of stellar mass at 1.9−0.5 M⊙ and a rapid decrease for more massive stars. The latter study uses a sample of 166 giant stars and finds a consistent result: occurrence rate increases with stellar mass, with the maximum at 2.1 M⊙. In the previous works mentioned, masses for the evolved stars were estimated by comparing ob- served properties (effective temperature, metallicity, and luminosity or surface gravity) with grids of stellar evolution models (e.g., Bressan et al. 2012). Although the accuracy of this method has been thoroughly tested for main-sequence stars (e.g., Torres et al. 2010), results for sugbiants and giants are still debated, as is clear from the above discussion. In an additional effort to test if stellar evolution models are also reliable for evolved stars, Ghezzi & Johnson (2015) showed that PARSEC evolutionary tracks (Bressan et al. 2012) coupled with the code PARAM (da Silva et al. 2006) can recover model-independent masses from eclipsing binaries and asteroseismology within ∼4% for the interval ∼0.7 – 4.5 M⊙. Therefore, the method itself does not present any issues that would signif- Retired A Stars Revisited 5 icantly overestimate the masses of subgiants and giants. Nevertheless, incorrect input parameters could still lead to erroneous mass determinations. To avoid model-dependent issues, most recent efforts to solve the Retired A Stars controversy focused on the determination of empirical masses. Johnson et al. (2014) performed multiple mass measurements for the bright giant HD 185351 using spectroscopic, interferometric and asteroseismic data. Although discrepancies on the 2.6σ level were observed, all determined values confirm that
HD 185351 is more massive than 1.5 M⊙. Campante et al. (2017) determined an asteroseismic mass for the planet-hosting subgiant HD 212771 that is consistent with recent spectroscopic estimates and also with its classification as a Retired A Star. Stassun et al. (2017) determined virtually model independent masses (from empirical radii and density or spectroscopic surface gravity) for 358 planet- hosting stars, achieving a precision better than 15% for 134 of them. Within this most accurate sample, 30 stars lie in the HR diagram region typically occupied by the Retired A Stars and ∼80% of them have masses consistent with this classification. North et al. (2017) and Stello et al. (2017) investigated the masses of evolved stars using aster- oseismology. The former study show that five of the seven analyzed objects (1.0 – 1.7 M⊙) have spectroscopic masses slightly larger than the corresponding asteroseismic value. However, the au- thors claim the offset is not significant and highly dependent on the adopted literature mass. In a similar comparison, the latter study reveals a 15-20% offset for six (out of seven) stars with masses larger than 1.6 M⊙. Both studies highlight a large scatter among literature masses and trace the differences back to the input parameters and their uncertainties used in the determination. In this work, we revisit the Retired A Stars sample with the goal of checking if their original masses were in fact overestimated, thereby an artificial correlation between this physical parameter and the occurrence rate of giant planets. We also take the opportunity to investigate the concerns related to the kinematics of these stars. The paper is organized as follows. The sample is presented in Section 2 and new spectroscopic and kinematical analyzes are described in Section 3. Our results are validated in Section 4, including comparisons with the previous ones from Johnson et al. (2010a), other literature sources and other methods. In Section 5, we discuss the consistency between masses 6 Ghezzi, Montet & Johnson and space and rotational velocities and obtain an updated relation for the occurrence rate of giant planets as a function of stellar metallicity and mass. Finally, our concluding remarks are presented in Section 6.
2. SAMPLE AND DATA
Our sample of subgiants consists of 245 stars that have been continuously monitored by the Lick and Keck subgiant planet surveys since 2004 and 2007, respectively (see Table 1). The details regarding the sample selection are described in Johnson et al. (2006, 2010b). In summary, targets were selected from the Hipparcos catalog (ESA 1997; van Leeuwen 2007)) according to the following criteria: 0.5
< MV < 3.5 and 0.55 < B − V < 1.10 and V . 8.5. Stars lying less than 1 mag above the main sequence defined by Wright (2005) or in the clump (B − V > 0.8 and MV < 2.0) region are excluded. High-resolution spectra for stars in both the Lick and Keck surveys were obtained with the HIRES (High Resolution Echelle Spectrometer; Vogt et al. 1994) spectrograph on the Keck I 10-m telescope (Mauna Kea, Hawaii). Deckers B1, B2, B3, B5, C2 and E2 were used with a binning of 3x1. This instrumental setup produces resolutions R ≃ 50,000 – 100,0001 and an almost complete spectral coverage from ∼3600 A˚ to ∼7990 A,˚ except for an inter-detector gap from ∼6420 A˚ to ∼6543 A˚ and some inter-order spacings redward of 6600 A.˚ We use the “template” spectra taken without the iodine cell. They were reduced with the Keck pipeline following standard procedures. We measured the signal-to-noise (S/N) values using 112 apparent continuum regions between 5220 A˚ and 6860 A,˚ carefully selected using the spectra of the Sun (reflected off Vesta) and HD 185351 as references. Our targets have S/N & 100 and the typical value is ∼220 (see Table 1).
3. ANALYSIS AND RESULTS
3.1. The Line List
We compiled the initial line list for Fe I and Fe II from multiple sources: Sousa et al. (2008), Ghezzi et al. (2010a), Schuler et al. (2011), Tsantaki et al. (2013), Sousa et al. (2014), Liu et al.
1 https://www2.keck.hawaii.edu/inst/hires/slitres.html Retired A Stars Revisited 7
(2014) and Bedell et al. (2014). Only spectral features with λ > 5000 A˚ are selected because of line crowding for lower wavelengths. Lines in the intervals 6270 A˚ ≤ λ ≤ 6330 A˚ and λ> 6865 A˚ are removed in order to avoid possible contamination by telluric lines. Using the HIRES solar spectrum as a reference, we selected only those Fe I lines that were relatively isolated and unblended. We were more flexible for the Fe II case due to the more limited number of lines available for this species. The HIRES solar spectrum was also used as a reference to exclude lines located in the inter-detector or inter-order gaps.
We retrieved all atomic parameters (wavelength λ, excitation potential χ, log gf and van der Waals damping factor) from the Vienna Atomic Line Database (VALD; Ryabchikova et al. 2015 and ref- erences therein) on January 2016. We measured equivalent widths (EWs) for all lines on the Solar Flux Atlas (Kurucz et al. 1984) with the updated version of ARES (Sousa et al. 2015a). We adopted the following parameters: smoothder = 4, space = 3.0, rejt = 0.999, lineresol = 0.1 and miniline = 5. A few lines presented clearly wrong EWs and these were replaced by manual measurements done using the task splot in IRAF2. Using the above list, a Kurucz ATLAS9 ODFNEW model atmosphere (Castelli & Kurucz 2004)
3 for the Sun (effective temperature Teff = 5777 K, surface gravity log g= 4.44, metallicity [Fe/H] = 0.00 and microturbulence ξ = 1.00 km s−1) and the driver abfind of the July 2014 version of the Local Thermodynamic Equilibrium (LTE) line analysis code MOOG4 (Sneden 1973), we derived individual Fe abundances for all lines. Option 1 was used for the treatment of the damping in MOOG, i.e. tabulated van der Waals damping factors (Barklem et al. 2000; Barklem & Aspelund-Johansson 2005) are used when available within the code and values from VALD (see Table 2) are adopted otherwise. We removed the Fe I lines that returned absolute abundances A(Fe)5 lower than 7.20 or higher than 7.80. A similar cut was not applied to Fe II due to the more restricted number of lines for this species. We also excluded lines with reduced equivalent widths log(RW) ≡ log(EW/λ)
2 Image Reduction and Analysis Facility (IRAF) is distributed by the National Optical Astronomy Observatory (NOAO), which is operated by the Association of Universities for Research in Astronomy, Inc. (AURA) under coop- erative agreement with the National Science Foundation (NSF). 3 [Fe/H]= log(NF e/NH )⋆ - log(NF e/NH )⊙ 4 http://www.as.utexas.edu/ chris/moog.html 5 A(Fe) = log(NF e/NH )+12 8 Ghezzi, Montet & Johnson lower than -6.0 and larger than -4.8 in order to avoid too weak or saturated lines, respectively. Finally, we derived solar log gf values by imposing that all remaining lines returned A(Fe) = 7.50 (Asplund et al. 2009) after running MOOG with the same model atmosphere as above. They were used in the analysis of our sample (see Section 3.2) in order to improve the precision of the results, since laboratory and theoretical log gf values for many iron lines still suffer from large uncertainties (e.g., Mashonkina et al. 2011). The final line list contains 158 Fe I and 18 Fe II lines (see Table 2). Using the the updated version of ARES (Sousa et al. 2015a), we automatically measured equivalent widths for the lines in this final list for all 245 target stars. The parameters were the same as those used for the Solar Atlas, except for rejt, which was chosen according to the S/N value of each individual spectrum.
3.2. Atmospheric Parameters
We determined a homogeneous set of atmospheric parameters (Teff , log g, [Fe/H] and ξ) for our target stars using the standard spectroscopic method, which is based on the excitation and ionization equilibria of Fe I and Fe II. As for the Solar Atlas, the analysis was performed in LTE using the July
2014 version of MOOG with option 1 for the damping parameter. We used the ATLAS9 ODFNEW grid from Castelli & Kurucz (2004) to obtain interpolated model atmospheres. We obtained the final parameters for each star through an automated iterative process that has to simultaneously meet four criteria: zero slopes in the linear fits between A(Fe I) and χ (excitation equilibrium) and between A(Fe I) and log(RW); same average values of A(Fe I) and A(Fe II) (ioniza- tion equilibrium); and same value for the metallicity in the input model atmosphere and the output result from MOOG. The convergence in each criterion is achieved by successive changes in Teff , ξ, log g and [Fe/H], respectively. Note that [Fe/H] = hA(Fe I)⋆i - A(Fe)⊙, where we adopted A(Fe)⊙ = 7.50 (Asplund et al. 2009). The iterative process also included two rounds of 2σ clipping in order to remove lines that returned abundances too discrepant from the average values. We calculated uncertainties in the atmospheric parameters in the following way. The microturbu- lence was changed until the linear fit between A(Fe I) and log(RW) had a slope equal to the error of the zero slope from the converged solution. The difference between the this new ξ and the best Retired A Stars Revisited 9 value is adopted as the uncertainty in this parameter. A similar procedure was done for the effective temperature, but the linear fit between A(Fe I) and χ was used instead. Moreover, the contribution from the error in ξ is taken into account by varying this parameter by its uncertainty determined in the previous step and checking what is the slope produced in the linear fit between A(Fe I) and
χ. Then, Teff is varied until we recover a zero slope and the difference between this temperature and best one is taken as the contribution from the error in ξ. The uncertainty in the surface gravity is obtained by varying this parameter until the difference between A(Fe I) and A(Fe II) is equal to the standard deviation of the mean for the latter abundance, taken from the converged solution. The contribution of Teff to the error on log g was also considered (in a procedure similar to the one used above to estimate contribution of ξ to the uncertainty on Teff ). Finally, the uncertainty in [Fe/H] takes into account the standard deviation of mean for A(Fe I) and the variations caused by the errors in Teff , ξ and log g, all added in quadrature. As a first test for our method, we analyzed the HIRES spectrum of sunlight reflected off Vesta (S/N = 267) as if it was a regular target star. The only difference is that only one round of 2σ clipping was performed. Otherwise, many good lines would be excluded due to the low value of σ (as expected, since the log gf values were adjusted using the Sun as a reference). The determined atmospheric parameters are: Teff = 5778 ± 20 K, log g = 4.45 ± 0.07, [Fe/H] = 0.01 ± 0.01 and ξ = 1.009 ± 0.030 km s−1. They are all in excellent agreement with the canonical solar values. The final atmospheric parameters for our target stars are shown in Figure 1. The individual values as well as the associated uncertainties can be seen in Table 3. We should note that the relatively low errors presented on this table are the internal uncertainties of the spectroscopic analysis. More realistic external uncertainties for the standard spectroscopic method employed here (used in the discussion of Section 4.5.1) should be higher, with typical values of up to ∼100 K for Teff , ∼0.2 dex for log g, ∼0.2 dex for [Fe/H] and ∼0.2 km s−1 for ξ (e.g., Ghezzi et al. 2010a; Hinkel et al. 2016). For completeness, we show in Table 4 the sensitivities of the atmospheric parameters to these larger uncertainties for HD 185351 (first spectrum), taken as a representative star in our sample. 10 Ghezzi, Montet & Johnson