Arxiv:2006.04811V2 [Astro-Ph.SR] 16 Jun 2020

Draft version June 18, 2020 Typeset using LATEX twocolumn style in AASTeX61 BAFFLES: BAYESIAN AGES FOR FIELD LOWER-MASS STARS S Adam Stanford-Moore,1 Eric L. Nielsen,1, 2 Robert J. De Rosa,1, 3 Bruce Macintosh,1 and Ian Czekala4, ∗ 1Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, Stanford, CA 94305, USA 2Department of Astronomy, New Mexico State University, P.O. Box 30001, MSC 4500, Las Cruces, NM 88003, USA 3European Southern Observatory, Alonso de Córdova 3107, Vitacura, Santiago, Chile 4Department of Astronomy, 501 Campbell Hall, University of California, Berkeley, CA 94720-3411, USA ABSTRACT Age is a fundamental parameter of stars, yet in many cases ages of individual stars are presented without robust estimates of the uncertainty. We have developed a Bayesian framework, BAFFLES, to produce the age posterior for a 0 star from its calcium emission strength (log(RHK )) or lithium abundance (Li EW) and B − V color. We empirically determine the likelihood functions for calcium and lithium as functions of age from literature measurements of stars in benchmark clusters with well-determined ages. We use a uniform prior on age which reflects a uniform star formation rate. The age posteriors we derive for several test cases are consistent with literature ages found from other methods. BAFFLES represents a robust method to determine the age posterior probability distribution for any field star with 0 0:45 ≤ B − V ≤ 0:9 and a measurement of RHK and/or 0:35 ≤ B − V ≤ 1:9 and measured Li EW. We compile colors, 0 RHK , and Li EW from over 2630 nearby field stars from the literature and present the derived BAFFLES age posterior for each star. Keywords: Stellar activity (1580), Stellar ages (1581), Field stars (2103), Bayesian statistics (1900) arXiv:2006.04811v2 [astro-ph.SR] 16 Jun 2020 Corresponding author: S Adam Stanford-Moore [email protected] ∗ NASA Hubble Fellowship Program Sagan Fellow 2 1. INTRODUCTION ized particles in its stellar wind, which carry away angu- Age, along with mass and metallicity, is a fundamental lar momentum (Kraft 1967; Weber & Davis 1967), and parameter of stars. Accurate stellar ages are needed in as a result the magnetic field strength and thus calcium a wide variety of astronomical studies, including galac- emission strength generally decrease with age. 0 tic evolution, globular clusters, open clusters, star form- The index RHK is a measure of the flux in the nar- ing regions, stellar multiples, brown dwarf companions, row emission line in the core of the Calcium II H and ˚ ˚ and planetary systems. For direct imaging exoplanet K absorption lines at ∼3968 A and ∼3934 A respec- 0 surveys, such as the Gemini Planet Imager Exoplanet tively (Noyes et al. 1984; Wright et al. 2004). RHK is Survey (Macintosh et al. 2018; Nielsen et al. 2019), stel- derived from an intermediate index, the S index, which lar age is important at all stages of the survey. First, represents the ratio of the narrow emission flux to the while selecting target stars, younger stars are preferred, background continuum flux. S provides a relative com- since their planets will be inherently brighter and eas- parison of emission strength, yet includes both chromo- ier to detect. Second, the mass for an imaged planet is spheric and photospheric contributions and is dependent derived from the age of the host star using evolutionary on B − V (as well as age). Therefore, to remove the de- models that link mass, age, and luminosity (e.g. Allard pendencies on B − V the S index is transformed by two 2014; Baraffe et al. 2015), and the dominant measure- empirically determined polynomials in B − V , resulting 0 ment uncertainty in deriving mass is from age (Bowler in RHK (Noyes et al. 1984; Wright et al. 2004), where the 2016). Third, age is a requirement for measuring the polynomials have been calibrated over a B − V range of occurrence rates of planets. Translating sensitivity in 0.45 to 0.90, corresponding to an approximate spectral apparent brightness to mass sensitivity requires the age type range of F6 to K2. In addition to the long-term of each observed star. Thus completeness to planets as decline in activity over time, the S value for a single a function of mass, a key ingredient for occurrence rate, star also varies by ∼10% over that star's activity cycle relies heavily on precise ages for the entire sample (e.g. (Wright et al. 2004). Bowler 2016; Nielsen et al. 2013, 2019). Yet unlike the ages of stars in co-eval groups, like open 1.1.2. Lithium Equivalent Width clusters or moving groups, the ages of field stars are diffi- The strength of the lithium absorption line traces the cult to determine robustly. For stellar clusters with well- amount of lithium present in the photosphere of a star. determined membership lists, the main sequence turnoff When stars initially form, their primordial lithium abun- is used to robustly determine the age (e.g. Goudfrooij dances are similar to the abundance from Big Bang nu- et al. 2014; Cummings & Kalirai 2018). The lithium de- cleosynthesis, with number densities of ∼10−9 that of pletion boundary is applicable to both clusters and more hydrogen (Sestito & Randich 2005). Over time stars sparse moving groups, with the reddest objects in an as- deplete their primordial lithium by nuclear burning in sociation with detectable lithium absorption setting the the core and convective mixing, so that measurements overall age (e.g. Burke et al. 2004; Soderblom 2010). of remaining surface lithium correlate with stellar age For isolated field stars, however, a less robust set of ob- (Soderblom 2010; Skumanich 1971). For stars cooler servables that track age are available, including spectro- than ∼7000 K, lithium abundance can be measured scopic indicators (e.g. Skumanich 1971; Wright et al. based on the equivalent width (EW) of the absorption of 2004), gyrochronology (e.g. Kraft 1967; Barnes 2009), the lithium doublet at 6708 A(˚ Soderblom 2010); hotter and asteroseismology (e.g. Cunha et al. 2007). Here we stars (OBA spectral types) have ionized their lithium, present a Bayesian method to determine age through and have negligible 6708 A˚ absorption even with no two spectral indicators: calcium emission strength and lithium burning. the depth of the lithium absorption line. Lithium's two isotopes, 6Li and the more abundant 7Li, burn at temperatures of 2.2 million K and 2.6 1.1. Empirical Age Indicators million K, respectively. Since stellar surface tempera- 1.1.1. Calcium Emission Strength tures are much lower (∼2500 K for low-mass M stars Calcium emission strength, as given by the index and ∼46000 K for high-mass O-stars), in order to burn, 0 RHK , is connected to the strength of a star's magnetic lithium must be brought into hotter layers via convec- field through the stellar dynamo (Noyes et al. 1984). tion (Soderblom et al. 1990). As a result, the rate of The rotation of the star and convection within induces a lithium depletion largely depends on the depth of the magnetic field whose strength is proportional to the rate convection zone, allowing lower-mass stars { which while of rotation (Noyes et al. 1984; Skumanich 1971). Over having lower surface temperatures have much deeper time the star's rotation inevitably slows as it ejects ion- convective layers { to deplete lithium faster than higher- 3 mass stars (Soderblom et al. 1990). In addition to con- to each, add 0.1 dex of Gaussian noise to each value, vection it is thought that slow mixing induced by rota- and use the polynomial to convert back to age. The tion and angular momentum loss may affect lithium de- returned age distribution has a significant spike at .1 pletion (Sestito & Randich 2005), so that lithium abun- Gyr, which becomes more prominent as the amplitude dance is a function of age, spectral type, and the initial of the Gaussian noise is increased. Soderblom et al. rotation rate and rotational evolution of the star. (1991) attempted to correct for this effect by adjusting the polynomial fit at large ages, by constraining it with 0 the nearby star sample and assuming that sample had 1.2. Functional fits to RHK and Li EW evolution a uniform star formation rate. Mamajek & Hillenbrand Previous studies have taken advantage of the corre- (2008) advocates for an activity/rotation/age relation lation between R0 and age to create empirical fits of HK instead, which flattens out the age distribution of the mean cluster log(R0 ) vs. log cluster age (Soderblom HK volume-limited sample out to 6 Gyr. Here, we present et al. 1991; Donahue 1993; Lachaume et al. 1999; Ma- an explicit prior uniform in age when creating age pos- majek & Hillenbrand 2008). However while these poly- terior probability density functions to address this issue. nomial fits allow one to find an expression for age as a Furthermore, the median age estimates for separate function of R0 , the polynomial makes no direct pre- HK methods (e.g. R0 and Li EW) are difficult to rig- diction of uncertainty in the age derivation. Soderblom HK orously combine without precise uncertainty estimates. et al.(1991) found the standard deviation of their stel- Previous works have, for example, simply averaged the lar data around their power-law fit to be ∼0.2 dex and ages obtained from R0 and lithium e.g. Nielsen & concluded that ages predicted from their fit would be HK Close 2010. accurate to ∼50%. A similar approach has been done Brandt et al.(2014b) develops a Bayesian method to for lithium as well (finding an average fit to clusters, combine the age PDF of a star's likely moving group and assigning a single age to a star based on its loca- with its posterior PDF from indicators of chromospheric tion relative to the cluster fits), e.g.

Load more