Arxiv:2012.02211V1 [Astro-Ph.EP] 3 Dec 2020 Tected and Used to Constrain the Masses of Planets in the Sing Et Al
Total Page:16
File Type:pdf, Size:1020Kb
Draft version December 7, 2020 Typeset using LATEX twocolumn style in AASTeX63 Decaying Orbit of the Hot Jupiter WASP-12b: Confirmation with TESS Observations Jake D. Turner ,1 Andrew Ridden-Harper ,1 and Ray Jayawardhana 2 1Department of Astronomy and Carl Sagan Institute, Cornell University, Ithaca, New York 14853, USA 2Department of Astronomy, Cornell University, Ithaca, New York 14853, USA Submitted to AJ ABSTRACT Theory suggests that the orbits of some close-in giant planets should decay due to tidal interactions with their host stars. To date, WASP-12b is the only hot Jupiter reported to have a decaying orbit, at a rate of 29±2 msec year−1. We analyzed data from NASA's Transiting Exoplanet Survey Satellite (TESS) to verify that WASP-12b's orbit is indeed changing. We find that the TESS transit and occul- tation data are consistent with a decaying orbit with an updated period of 1.091420090±0.000000041 days and a decay rate of 32.53±1.62 msec year−1. We find an orbital decay timescale of τ = P/jP_ j = 2.90±0.14 Myr. If the observed decay results from tidal dissipation, the modified tidal quality factor 5 is Q'? = 1.39±0.15 ×10 , which falls at the lower end of values derived for binary star systems and hot Jupiters. Our result highlights the power of space-based photometry for investigating the orbital evolution of short-period exoplanets. Keywords: planets and satellites: gaseous planets { planet-star interactions 1. INTRODUCTION Penev et al. 2014; Barker 2020). As the planet's orbit de- Some exoplanetary systems exhibit temporal varia- cays over millions of years, its orbital period will change tions in their orbital parameters. Such systems provide by a small, yet potentially detectable amount. Obser- valuable insights into their physical properties and the vations of this decay are now possible because some processes that cause the variations. Light curves of tran- hot Jupiter systems have been monitored for decades. siting planets can be used to search for transit timing Such measurements enhance our understanding of the variations (TTVs; Agol et al. 2005; Agol & Fabrycky hot Jupiter population (e.g. Jackson et al. 2008; Hamer 2018) and, less frequently, for variations in transit dura- & Schlaufman 2019). tion (Agol & Fabrycky 2018) and the impact parameter Currently, WASP-12b is one of the few hot Jupiters (Herman et al. 2018; Szab´oet al. 2020). The presence of confirmed to have a varying period. It orbits a G0 star TTVs can indicate additional bodies in the system or an with a short period of one day (Hebb et al. 2009). Since unstable orbit resulting from tidal forces of the star (e.g., its discovery in 2009, WASP-12b's orbital parameters Miralda-Escud´e 2002; Mazeh et al. 2013). Also, TTVs (e.g., Campo et al. 2011; Maciejewski et al. 2016) and in tightly packed multi-planetary systems have been de- atmosphere (e.g., Croll et al. 2011; Stevenson et al. 2014; arXiv:2012.02211v1 [astro-ph.EP] 3 Dec 2020 tected and used to constrain the masses of planets in the Sing et al. 2016) have been studied extensively. Obser- system (e.g., Agol & Fabrycky 2018). vations with the Cosmic Origins Spectrograph (COS) on Theory suggests that short-period massive planets or- the Hubble Space Telescope reveal an early ingress in the biting stars with surface convective zones may exchange near-ultraviolet (Fossati et al. 2010; Haswell et al. 2012; energy with their host stars through tidal interactions, Nichols et al. 2015). An escaping atmosphere was sug- causing the host star to spin faster and the planet's orbit gested as the cause of the early ingress (e.g., Lai et al. to decay (e.g., Lin et al. 1996; Chambers 2009; Lai 2012; 2010; Bisikalo et al. 2013; Turner et al. 2016a). Ma- ciejewski et al.(2016) were the first to report evidence of a decreasing orbital period for WASP-12b. At the Corresponding author: Jake D. Turner time, the cause of the changing orbit was uncertain, and [email protected], [email protected] could be ascribed to orbital decay, apsidal precession, or 2 J.D. Turner, A. Ridden-Harper, R. Jayawardhana Figure 1. TESS light curve of WASP-12b in Sector 20. Top: Raw simple aperture photometry light curves. Bottom: Detrended Presearch Data Conditioning (PDC) time series. the Romer effect. Follow-up observations confirmed the 2. OBSERVATIONS AND DATA REDUCTION changing orbital period, with models slightly favoring TESS observed WASP-12b in Sector 20 (December orbital decay as the cause (Patra et al. 2017; Maciejew- 24, 2019 to January 20, 2020; Figure1). These observa- ski et al. 2018; Bailey & Goodman 2019). Recently, tions were processed by the Science Processing Opera- using new transit, occultation and radial velocity obser- tions Center (SPOC) pipeline (Jenkins et al. 2016). The vations, Yee et al.(2020) presented strong evidence that final product of the SPOC pipeline are light curves cor- WASP-12b's period variation is caused by orbital decay. rected for systematics that can be used to characterize The decay rate has been constrained to between 20-30 transiting planets. All of the SPOC data products are milliseconds/yr. publicly accessible from the Mikulski Archive for Space Motivated by these indications of WASP-12b's chang- Telescopes (MAST)1. The Presearch Data Condition- ing period, we use NASA's Transiting Exoplanet Survey ing (PDC) component of the SPOC pipeline corrects Satellite (TESS; Ricker et al. 2015) observations to ver- the light curves for pointing or focus related instrumen- ify its orbital decay, and derive updated orbital param- tal signatures, discontinuities resulting from radiation eters and planetary properties. TESS is well suited for events in the CCD detectors, outliers, and flux contam- our study because it provides high-precision time-series ination. We use the PDC light curve for the analysis in data, ideal for searching for TTVs (e.g. Hadden et al. 2019; von Essen et al. 2020; Ridden-Harper et al. 2020). 1 https://archive.stsci.edu/ Decaying Orbit of the hot Jupiter WASP-12b 3 Table 1. Physical properties of WASP-12b derived from the this paper, however a known issue2 results in the SPOC light curve modeling of the TESS data pipeline overestimating the uncertainties for some light curves, including that of WASP-12b. Therefore, it is rec- Parameter units value 1 σ uncertainty ommended that the uncertainties be estimated from the Rp/R∗ 0.1166445 0.0000097 scatter (Barclay, T., private communication), until the a/R∗ 3.1089 0.0014 data is reprocessed with an updated pipeline (data re- δocc 0.00048 0.00010 leased after DR 38 should be unaffected). To do this, we Inclination ◦ 84.955 0.037 used the standard deviation of the entire out-of-transit Duration mins 178.60 0.14 baseline as the uncertainty. Using smaller out-of-transit b 0.2734 0.0020 baselines surrounding only a few transits did not signif- R R 1.884 0.057 icantly change the standard deviation, indicating that b Jup Mb MJup 1.46 0.27 the noise is relatively constant. The raw light curve from −3 the SPOC and detrended PDC light curve are shown in ρb g cm 0.271 0.056 Figure1. log gb cgs 2.38 0.39 −3 ρ∗a g cm 0.477 0.026 3. DATA ANALYSIS Teq K 2551 56 3.1. Transit Modeling Θ 0.0260 0.0053 a au 0.02399 0.00072 We modeled the TESS transits of WASP-12b with the EXOplanet MOdeling Package (EXOMOP; Pearson et al. To obtain the final fitted parameters, the light curve of 3 2014; Turner et al. 2016c, 2017) to find a best-fit. WASP-12b was phase-folded at each derived mid-transit EXOMOP creates a model transit using the analytic equa- time and modeled with EXOMOP. The phase-folded light tions of Mandel & Agol(2002) and the data are mod- curve and model fit can be found in Figure2. We use the eled using a Differential Evolution Markov Chain Monte light curve model results combined with literature values Carlo (DE-MCMC; Eastman et al. 2013) analysis. To to calculate the planetary mass (Mb; Winn 2010), radius account for red noise in the light curve, EXOMOP incor- (Rb), density (ρb), equilibrium temperature (Teq), sur- porates the residual permutation, time-averaging, and face gravity (log gb; Southworth et al. 2007), orbital dis- wavelet methods. More detailed descriptions of EXOMOP tance (a), inclination (i), Safronov number (Θ; Safronov can be found in Pearson et al.(2014) and Turner et al. 1972; Southworth 2010), and stellar density (ρ∗a; Sea- (2016b). ger & Mall´en-Ornelas 2003). The planet properties and Each transit in the TESS data (lower panel of Figure transit ephemeris we derived for WASP-12b are shown 1) was modeled with EXOMOP independently. We used in Table1. All the planetary parameters are consistent 6 20 chains and 20 links for the DE-MCMC model and with their discovery values (Hebb et al. 2009) but their ensure chain convergence (Ford 2006) using the Gelman- precision is greatly improved. Rubin statistic (Gelman & Rubin 1992). The mid- 3.2. Occultation Modeling transit time (Tc), planet-to-star radius (Rp=R∗), incli- nation (i), and scaled semi-major axis (a=R∗) are set as The TESS occultation was also modeled with EXOMOP free parameters for every transit. The period (Porb) and to find the best-fit parameters. The light curve we fit linear and quadratic limb darkening coefficients are fixed was obtained by phase-folding all the data about the during the analysis.