Draft version October 7, 2020 Typeset using LATEX twocolumn style in AASTeX63 When Do Stalled Stars Resume Spinning Down? Advancing Gyrochronology with Ruprecht 147 Jason Lee Curtis,1,2 Marcel A. Agueros¨ ,1 Sean P. Matt,3 Kevin R. Covey,4 Stephanie T. Douglas,5 Ruth Angus,2,6,1 Steven H. Saar,5 Ann Marie Cody,7 Andrew Vanderburg,8,5 Nicholas M. Law,9 Adam L. Kraus,8 David W. Latham,5 Christoph Baranec,10 Reed Riddle,11 Carl Ziegler,12 Mikkel N. Lund,13 Guillermo Torres,5 Søren Meibom,5 Victor Silva Aguirre,13 and Jason T. Wright14 1Department of Astronomy, Columbia University, 550 West 120th Street, New York, NY 10027, USA 2Department of Astrophysics, American Museum of Natural History, Central Park West, New York, NY, USA 3University of Exeter, Department of Physics & Astronomy, Stocker Road, Exeter, EX4 4QL, UK 4Department of Physics & Astronomy, Western Washington University, Bellingham, WA 98225-9164, USA 5Center for Astrophysics | Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA 6Center for Computational Astrophysics, Flatiron Institute, 162 5th Avenue, Manhattan, NY, USA 7Bay Area Environmental Research Institute, 625 2nd Street, Ste. 209, Petaluma, CA 94952, USA 8Department of Astronomy, The University of Texas at Austin, Austin, TX 78712, USA 9Department of Physics and Astronomy, University of North Carolina, Chapel Hill, NC 27599, USA 10Institute for Astronomy, University of Hawai‘i at M¯anoa, 640 N. A‘oh¯ok¯uPl., Hilo, HI 96720-2700, USA 11Division of Physics, Mathematics, and Astronomy, California Institute of Technology, Pasadena, CA 91125, USA 12Dunlap Institute for Astronomy and Astrophysics, University of Toronto, 50 St. George Street, Toronto, Ontario M5S 3H4, Canada 13Stellar Astrophysics Centre, Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, DK-8000 Aarhus C, Denmark 14Penn State Extraterrestrial Intelligence Center / Center for Exoplanets and Habitable Worlds / Department of Astronomy & Astrophysics, 525 Davey Laboratory, The Pennsylvania State University, University Park, PA, 16802, USA (Received April 10, 2020; Accepted September 28, 2020) ABSTRACT Recent measurements of rotation periods (Prot) in the benchmark open clusters Praesepe (670 Myr), NGC 6811 (1 Gyr), and NGC 752 (1.4 Gyr) demonstrate that, after converging onto a tight sequence of slowly rotating stars in mass–period space, stars temporarily stop spinning down. These data also show that the duration of this epoch of stalled spin-down increases toward lower masses. To determine when stalled stars resume spinning down, we use data from the K2 mission and the Palomar Transient Factory to measure Prot for 58 dwarf members of the 2.7-Gyr-old cluster Ruprecht 147, 39 of which satisfy our criteria designed to remove short-period or near-equal-mass binaries. Combined with the Kepler Prot data for the approximately coeval cluster NGC 6819 (30 stars with M⋆ > 0.85 M⊙), our new measurements more than double the number of ≈2.5 Gyr benchmark rotators and extend this sample down to ≈0.55 M⊙. The slowly rotating sequence for this joint sample appears relatively flat arXiv:2010.02272v1 [astro-ph.SR] 5 Oct 2020 (22±2 days) compared to sequences for younger clusters. This sequence also intersects the Kepler intermediate period gap, demonstrating that this gap was not created by a lull in star formation. We calculate the time at which stars resume spinning down, and find that 0.55 M⊙ stars remain stalled for at least 1.3 Gyr. To accurately age-date low-mass stars in the field, gyrochronology formulae must be modified to account for this stalling timescale. Empirically tuning a core–envelope coupling model with open cluster data can account for most of the apparent stalling effect. However, alternative explanations, e.g., a temporary reduction in the magnetic braking torque, cannot yet be ruled out. Keywords: stellar ages – stellar rotation – open clusters – stellar evolution Corresponding author: Jason Lee Curtis [email protected] 2 Curtis et al. 10.0 1. INTRODUCTION Observations of Sun-like stars revealed that they spin Ruprecht 147 NGC 6819 down over time via magnetic braking. This led to the development of gyrochronology, a promising age-dating NGC 752 technique that uses rotation periods (Prot) as a clock (Barnes 2003).1 Recent measurements of P in the 1.0 NGC 6811 rot Hyades benchmark open clusters Praesepe (670 Myr; Douglas Praesepe et al. 2017), the Hyades (730 Myr; Douglas et al. 2019), Age (Gyr) M48 NGC 6811 (1.0 Gyr; Curtis et al. 2019a), and NGC 752 (1.4 Gyr; Ag¨ueros et al. 2018), however, have shown that M34 M35 the formula describing the process of stellar spin-down Pleiades Blanco 1 cannot be as simple as it once appeared. Instead of Prot 0.1 evolving continuously as a power law with a braking in- 10 100 1000 10000 dex that is constant in time and common to all stars n Distance (pc) (i.e., Prot ∝ t , with Skumanich 1972 originally propos- ing n =0.5), it is now clear that, after converging onto Figure 1. Age vs. distance for a selection of benchmark a slowly rotating sequence that illustrates the tight re- open clusters with Prot data. Red stars indicate clusters used in this study; gray stars represent other notable clusters lationship between mass and Prot in clusters older than with P data mentioned in this paper. Ruprecht 147 and ≈100 Myr, stars temporarily stop spinning down. Fur- rot NGC 6819 have similar ages, but Ruprecht 147 is much closer thermore, based on where the period sequences for these to Earth and its stars are ≈77× brighter. This allows us to benchmark clusters overlap, Ag¨ueros et al. (2018) con- measure Prot for Ruprecht 147 stars with much lower masses cluded that the duration of this epoch of stalled spin- than was possible with the Kepler survey of NGC 6819. down increases toward lower stellar masses. Most empirical gyrochronology relations (e.g., Barnes 2007, 2010) and angular momentum evolution mod- analyze photometric time series data from NASA’s K2 els (e.g., van Saders & Pinsonneault 2013; Matt et al. mission (Howell et al. 2014) and the Palomar Transient 2015; Gallet & Bouvier 2015; Lanzafame & Spada 2015), Factory (PTF; Law et al. 2009; Rau et al. 2009), and which predate these new observational findings, do not measure Prot for 58 dwarf members with masses as low account for the phenomenon of stalling, so ages inferred as M⋆ ≈ 0.55 M⊙. In Section 4, we identify 35 of these for low-mass stars with these models will likely be in- as benchmark rotators, including 23 likely single stars, correct. It is therefore imperative that we extend the discuss the impact of binarity on the color–period distri- benchmark sample of cluster rotators to lower masses bution, and present a rotation catalog for Ruprecht 147. and older ages. This will allow us to constrain the Our sample of rotators in Ruprecht 147 contains few > timescale for stalled spin-down and to eventually repair stars with masses ∼1 M⊙. To remedy this, we com- empirical gyrochronology models, so that accurate ages bine our sample with Prot measurements for the approx- can be inferred for low-mass stars in the field. These new imately coeval (2.5 Gyr) cluster NGC 6819, which was data are also needed to tune free parameters in the more surveyed during the primary Kepler mission (Meibom physically motivated class of gyrochronology models (for et al. 2015). Because stars in NGC 6819 were moni- tored 10-16× longer than those in Ruprecht 147, and a recent example using our new Prot data for NGC 6811 from Curtis et al. 2019a, see Spada & Lanzafame 2020). because the Kepler light curves contain negligible sys- tematics compared to our K2 Campaign 7 data, Meibom We present measurements of Prot for the 2.7-Gyr-old open cluster Ruprecht 147, the oldest nearby open clus- et al. (2015) were able to recover weaker amplitude sig- ter (see Figure 1). In Section 2, we discuss the prop- nals in the NGC 6819 data than we could detect in the erties of and membership criteria for this cluster, and Ruprecht 147 light curves. This enabled Meibom et al. then identify binaries among, and estimate stellar pa- (2015) to measure Prot for early G dwarfs in NGC 6819, rameters for, the targets of this study. In Section 3, we although the cluster’s larger distance modulus restricted their overall Prot sample to M⋆ > 0.85 M⊙. In this study, we contribute rotation periods for 20 benchmark 1 We define Sun-like broadly, to include all stars with radiative rotators with masses smaller than this lower limit for cores and convective envelopes; i.e., 0.4 . M⋆ . 1.3 M⊙, or spectral types early M to late F. These stars have solar-like dy- NGC 6819. namos and undergo magnetic braking. Advancing Gyrochronology with Ruprecht 147 3 Pairing the rotation data for NGC 6819 and We calculated the average interstellar reddening by 2 Ruprecht 147 allows us to describe stellar rotation at comparing the Gaia DR2 colors, (GBP − GRP), with 2.5-2.7 Gyr as a function of mass from late F down to spectroscopic temperatures (Teff) derived from high- M1 spectral types. Section 5 describes this procedure, resolution spectroscopy using SME for nearby, unred- which involves calculating the average interstellar red- dened field stars (Valenti & Fischer 2005) and the same dening toward NGC 6819 relative to Ruprecht 147, and six solar twin members of Ruprecht 147 used to mea- estimating their relative ages with gyrochronology using sure the cluster metallicity. We found a reddening value the portions of each sample that overlap significantly in of E(B − V )=0.099 ± 0.010, or equivalently an ex- mass (i.e., mid-to-late G dwarfs).