Arxiv:2107.01213V1 [Astro-Ph.EP] 2 Jul 2021
Total Page:16
File Type:pdf, Size:1020Kb
Draft version August 17, 2021 Typeset using LATEX twocolumn style in AASTeX62 Hα and Ca ii Infrared Triplet Variations During a Transit of the 23 Myr Planet V1298 Tau c Adina D. Feinstein,1, ∗ Benjamin T. Montet,2, 3 Marshall C. Johnson,4 Jacob L. Bean,1 Trevor J. David,5, 6 Michael A. Gully-Santiago,7 John H. Livingston,8 and Rodrigo Luger5 1Department of Astronomy and Astrophysics, University of Chicago, 5640 S. Ellis Ave, Chicago, IL 60637, USA 2School of Physics, University of New South Wales, Sydney, NSW 2052, Australia 3UNSW Data Science Hub, University of New South Wales, Sydney, NSW 2052, Australia 4Las Cumbres Observatory, 6740 Cortona Drive, Suite 102, Goleta, CA 93117, USA 5Center for Computational Astrophysics, Flatiron Institute, 162 Fifth Ave, New York, NY 10010, USA 6Department of Astrophysics, American Museum of Natural History, New York, NY 10024, USA 7Department of Astronomy, The University of Texas at Austin, 2515 Speedway Boulevard, Austin, TX 78712, USA 8Department of Astronomy, University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan ABSTRACT Young transiting exoplanets (< 100 Myr) provide crucial insight into atmospheric evolution via photoevaporation. However, transmission spectroscopy measurements to determine atmospheric com- position and mass loss are challenging due to the activity and prominent stellar disk inhomogeneities present on young stars. We observed a full transit of V1298 Tau c, a 23 Myr, 5.59 R⊕ planet orbiting a young K0-K1.5 solar analogue with GRACES on Gemini-North. We were able to measure the Doppler tomographic signal of V1298 Tau c using the Ca ii infrared triplet (IRT) and find a projected obliquity of λ = 5◦ ± 15◦. The tomographic signal is only seen in the chromospherically driven core of the Ca ii IRT, which may be the result of star-planet interactions. Additionally, we find that excess ab- sorption of the Hα line decreases smoothly during the transit. While this could be a tentative detection of hot gas escaping the planet, we find this variation is consistent with similar timescale observations of other young stars that lack transiting planets over similar timescales. We show this variation can also be explained by the presence of starspots with surrounding facular regions. More observations both in- and out-of the transits of V1298 Tau c are required to determine the nature of the Ca ii IRT and Hα line variations. Keywords: Stellar activity, exoplanet astronomy, exoplanet atmospheres, transmission spectroscopy Features: Affiliated Python scripts/Jupyter notebooks - 1. INTRODUCTION tion, and measuring atmospheric mass-loss rate due to Young planets provide valuable insights into the early photoevaporation. stages of planet formation and evolution. Through tran- sit missions like K2 (Howell et al. 2014) and TESS (Tran- 1.1. Obliquity Measurements siting Exoplanet Survey Satellite; Ricker et al. 2014) the population of young transiting planets (here defined Close-in transiting exoplanets have brought into ques- arXiv:2107.01213v2 [astro-ph.EP] 16 Aug 2021 as < 100 Myr) has grown to approximately a dozen tion the mechanisms in which they get there. Measuring (David et al. 2016; Mann et al. 2016; Benatti et al. 2019; the obliquity, λ, or the relative angle between the spin Newton et al. 2019; Mann et al. 2020; Plavchan et al. of the star and orbit of the planet (Lai 2014), can yield 2020; Rizzuto et al. 2020). These planets are key to insights into different migratory paths. A low-obliquity understanding planetary migration, atmospheric evolu- hints at a smooth disk migration history (Goldreich & Tremaine 1979; Ford 2014), while binary companions or nearby stars in stellar birth clusters can torque planets Corresponding author: Adina D. Feinstein; to a high obliquity over millions of years (Fabrycky & 7 @afeinstein20; @afeinstein20; [email protected] Tremaine 2007). Such measurements have only just begun for young ∗ NSF Graduate Research Fellow systems. The recently discovered 20 Myr planet 2 Feinstein et al. AU Mic b (Plavchan et al. 2020) has been the subject of tions, Shkolnik & Barman(2014) found that early M significant follow-up. All independent studies have mea- stars remain at high, saturated levels of NUV and FUV sured a low obliquity (Addison et al. 2020; Hirano et al. flux for the first ∼650 Myr. The same relation was found 2020; Martioli et al. 2020; Palle et al. 2020), suggesting of K stars (Richey-Yowell et al. 2019).The activity life- the planet formed beyond the ice-line within the pro- time for low-mass stars is longer than that of solar-type toplanetary disk and experienced a smooth migration stars. However, the higher levels of XUV flux makes inward. photoevaporation more efficient for planets around low- Similar studies of spin-orbit alignment have been com- mass stars (Rogers & Owen 2021), thus the evolution of pleted for DS Tuc Ab, a planet in a slightly older (35- planetary radii is faster (Owen & Wu 2013). 45 Myr) system, in a known binary (Newton et al. 2019; An increase in magnetic activity leads to elevated Montet et al. 2020; Zhou et al. 2020). The results from UV/X-Ray radiation in young stars, manifests itself in both studies concluded the planet has a low obliquity larger spot coverage, and results in increased photomet- and that the migration was not significantly torqued ric variability (Feigelson & Montmerle 1999). Several through Kozai-Lidov oscillations by DS Tuc B. This is detailed studies of T Tauri stars have estimated that additionally supported by the fact that the age of the spots cover 29-41% of the less spotted hemisphere and system is younger than the timescale required for Kozai- 61-67% of the more spotted hemisphere (Grankin 1999), Lidov interactions (Montet et al. 2020). V1298 Tau has with the extreme case of LkCa 4 with of ∼ 80% spot no known companions and thus any spin-orbit misalign- coverage (Gully-Santiago et al. 2017). Larger studies ment would hint towards potential torquing of the disk of young stars have also seen a lack of phase depen- due to neighboring stars in the birth cluster. dence for flares, another proxy for magnetic activity, in- dicating similar levels of high spot coverage across both hemispheres (Feinstein et al. 2020b). For detailed spec- 1.2. Photoevaporation & Young Stellar Activity troscopic studies of young planets (e.g. studying their spin-orbit alignment or atmospheric characterization), Photoevaporation of young close-in planets is one of it is crucial to understand the underlying starspot and the leading explanations for the \radius gap" of transit- active region coverage. ing planets (Owen & Wu 2017; Fulton et al. 2017; Van Eylen et al. 2018; Bean et al. 2021). This gap is between a bimodal distribution of small planets, centered at 1.3 and 2.4 R⊕ with a distinct lack of planets in between. However, early studies did not account for the ages of 1.3. V1298 Tau the systems; Berger et al.(2020) presented the first anal- V1298 Tau is a 23 Myr pre-main sequence K star ysis of the planet-radius distribution while accounting (M?/M = 1.101; R?/R = 1.345) hosting four known for such ages. They found tentative evidence of an age transiting planets (David et al. 2019a,b) observed in K2 dependence of the planet radius distribution, demon- Campaign 4 (Howell et al. 2014). V1298 Tau is the first strating that the number of super-Earths (R/R⊕ < 1:8) known young transiting multi-planet system. V1298 Tau increases with age, assuming a fixed radius gap loca- was first identified as a young star via strong X-ray emis- tion. This could be attributed to a core-powered mass- sion observed with ROSAT and strong Li i 670.8 nm loss mechanism (Gupta & Schlichting 2019), which may absorption (Wichmann et al. 1996, 2000). V1298 Tau continue for the first 1 Gyr (Rogers & Owen 2021). belongs to a young association in the foreground of the However, the general contributions of atmospheric re- Taurus-Auriga star-forming complex (Oh et al. 2017) moval from photoevaporation and core-powered mass- with an estimated age of 20-30 Myr based on comparison loss is currently unknown. While photoevaporation rate to empirical and theoretical isochrones (Luhman 2018; is a function of flux, these young transiting planets have David et al. 2019b). V1298 Tau c (R ∼ 5:59R⊕) and d short orbital periods and allow us to probe the most (R ∼ 6:41R⊕) are predicted to evolve to sub-Neptunes dramatic environments for atmospheric removal. (< 4R⊕) or super-Earths (< 2R⊕), depending on the Atmospheric photoevaporation is driven by the high XUV evolution of the host star, after 5 Gyr (David et al. energy irradiation from young stars (Lammer et al. 2003; 2019a; Poppenhaeger et al. 2021). The light curve for Owen 2019). The effects are particularly dramatic in the V1298 Tau also shows strong starspot variability, with a first 100 Myr where stellar XUV emission is higher and peak-to-trough amplitude of ∼ 6% (David et al. 2019b). planets are still contracting, which makes their atmo- Therefore, any in-transit spectra need to be carefully spheres more vulnerable to mass loss (Preibisch et al. searched for the signatures of contaminating ARs and 2005; Feigelson et al. 2007). Using GALEX observa- spots. V1298 Tau Stellar Activity 3 V1298 Tau b is seen to be well aligned through both erage of 400-1000 nm and a resolving power of R ∼ a clear Rossiter-McLaughlin signal and Doppler tomo- 81,000. To obtain an optimal signal-to-noise of ap- graphic analyses (Johnson et al. in prep). Additionally, proximately 300 at 600 nm, we observed V1298 Tau at David et al. (in prep) reports detections of excess Hα 360 second exposures.