Transit and Occultation Data Confirm a Decaying Orbit for WASP-12b

Research Article
Open access

Transit and Occultation Data Confirm a Decaying Orbit for WASP-12b

Zoutong Shen 1* , Hao Meng 2 , Yishan Wu 3 , Mei Zhang 4 , Yunhe Fei 5
  • 1 Carl Sagen Institute, Cornell University, Ithaca, NY    
  • 2 School of Physics, Huazhong University of Science and Technology, Wuhan, China    
  • 3 College of Sciences, Georgia Institute of Technology, Atlanta, GA 30332, USA    
  • 4 School of Physics, East China University of Science and Technology, Shanghai 200237, China    
  • 5 Nanjing Foreign Language School, Nanjing, China    
  • *corresponding author felixzshen@gmail.com
TNS Vol.107
ISSN (Print): 2753-8826
ISSN (Online): 2753-8818
ISBN (Print): 978-1-80590-087-0
ISBN (Online): 978-1-80590-088-7

Abstract

WASP-12b is the first reported hot Jupiter whose orbital period appears to be decreasing with time. This phenomenon can be interpreted using two models from earlier research: the decay of its orbit or the apsidal precession. This work chose the latest 118 transit and 6 occultation photometric observations conducted by TESS in Sectors 20, 43, 44, 45, 71, 72. They were combined with previous data to analyze the relationship between the orbit number, the mid-time of transit, and the mid-time of occultations (averaged over entire sectors to boost the signal-to-noise ratio), which favored the decay model with ∆BIC = 104.09. Our value of decay rate , and the reduced tidal quality factor when interpreting the decay as tidal orbital decay, are both aligned with earlier studies and suggest that the orbit of WASP-12b is in a state of rapid decay. Future observations and theoretical work are necessary to further explain the physical processes driving this decay and to predict the end of life of WASP-12b.

Keywords:

Exoplanets, Hot Jupiters, Transits and Occultation, Time Deviation

Shen,Z.;Meng,H.;Wu,Y.;Zhang,M.;Fei,Y. (2025). Transit and Occultation Data Confirm a Decaying Orbit for WASP-12b. Theoretical and Natural Science,107,111-124.
Export citation

1. Introduction

Almost 30 years ago, the first hot Jupiter was discovered. A hot Jupiter is a kind of planet with Jupiter-like size and an unusually short orbit radius[1] Different from the normal gas giant planet whose orbit is far from its host star since more solid material is available for the planet to form a huge core and accrete much gas as the temperature decrease, a hot Jupiter occurs closely to the host star, sometimes less than 0.1 AU. One theory explaining why the orbits of hot Jupiters are located so close to the host star is that the giant planet initially forms in the protoplanetary disk, migrates slowly through the spinning disk [2,3]. the center, and finally reaches close to the host star. When the disk disappears, a hot Jupiter eventually occurs. Another theory suggests that a giant planet, located far from the central star after the disk disappears and following a nearly circular orbit like Jupiter, is excited by other planets. As a result, the primitive circular orbit turns into an elliptical orbit and finally forms a circular orbit much closer to its host star due to tidal effects. In this process, the orbital energy is converted by friction into heat within the planet[4].

WASP-12b is a planet reported with decaying orbit [5, p. 20], [6, p. 201]. This hot Jupiter follows an almost circular orbit [7,8] with a 1.09-day orbital period [9] a huge size of 1.90 RJup, and a mass of \( 1.47{M_{Jup}} \) [9]. Recent research shows that this hot Jupiter is still undergoing orbit decay with a decay rate of \( -29±3{msyr^{-1}} \) [6]. The decrease speed is so fast that the equilibrium tidal model which shows orbital decay rate as \( \frac{P}{\dot{P}}=\frac{(m_{ba}^{2})}{(2{I_{star}})}{t_{sync}}≈1.2Gyr \) is too long comparing with the rate calculated before [10, p. 201], [11] Several assumptions occurred to explain this phenomenon.

One theory maintains that as WASP-12b gets closer to WASP-12, the host star will interact with its hot Jupiter by the increased tidal dissipation due to nonlinear wave-breaking of dynamical tide. During this process, the friction turns the orbital energy into heat and shrink the orbit. However, this theory requires the main star to be a sub-giant rather than a main-sequence star and Bailey and Goodman [10] suggests that WASP-12 is more likely to be a main-sequence star rather than a sub-giant.

Meanwhile, Millholland and Laughlin [12] claims that WASP-12b’s orbital period’s rapid decrease is the result of the planetary obliquity tides. Influenced by an unseen, small-mass planet whose orbit is outside that of the hot Jupiter’s, WASP-12b maintains a large obliquity due to the spin-orbit resonance with the smaller planet. As a result, the energy of resonance is converted into heat and stored in WASP-12b due to the effect of friction.

One of the interesting possibilities is the tidal orbital decay, where the angular momentum of the hot Jupiter is directly transferred to its host star due to the gravitational interaction between the planet and the star’s tidal bulge (equilibrium tide). The mechanisms that convert stellar tidal oscillations into heat within the star are unclear [13].

In this paper, we present new transits and occultations data of WASP-12b, as well as previous data (Section 2). Based on the data, we compare the orbital decay model and the apsidal precession model for the better fitted one (Section 3). Finally we discuss the possible tidal orbital decay and prospects for further observations (Section 4).

2. Data collection

2.1. Transit observation

We used light curve data of WASP-12b from the Transiting Exoplanet Survey Satellite (TESS). TESS has four identical 10-cm-diameter cameras. Each has a field of view of \( 24×24 {deg^{2}} \) , combining in a line to \( 2304 de{g^{2}} \) field of view, corresponding to a sector on celestial sphere. The cameras detect the flux of light on their CCDs, which form images covering \( 4096×16384 \) pixels. The CCDs produce continuous stream of images with an exposure time of 2 s, and the images will be binned into a cadence for several optional periods. TESS faces a sector for two spacecraft orbits, each lasting for 13.7 days. During every orbit, it spends about 2.2 days sending its data back to the earth, when observation is interrupted [14]

We firstly chose the data sets with 120s cadence between 2019 December and 2023 December, which includes 6 sectors. To fit the data point, we used a standard model [15] to the data from each transit.

Then, we linked all of the latest 132 transit data sets, used a single eclipse model with WASP-12 and ruled out the time scope where no light curve was collected. Meanwhile, we dismissed 14 sets of data not including enough data points(less than 80% of the data median number of all data), leaving 118 sets of valid transit data.

We fixed the orbital period on 1.091420078 days and let other system parameters vary freely to get the best fit. Finally, we got all the fitting curves of valid data sets, as well as their mid-transit time and uncertainty. We chose the mid-transit time T0 = 2456305.455609 as a reference [6], so that we calculated the epoch, namely the orbit number of all transit data. Table 1 consists of the fitted parameters for each sector, along with the combined parameters.

Table 1: The fitted parameters for lightcurve data in each sector and combined parameters

Parameters

S20

S43

S44

S45

S71

S72

All

Fitted Parameters

\( \frac{{R_{P}}}{{R_{star}}} \)

0.11960

±0.00048

0.12005

±0.00044

0.11975

±0.00043

0.11856

±0.00047

0.11871

±0.00048

0.11842

±0.00050

0.11924

±0.00019

\( T_{c}^{a} \)

855.54862

±0.00015

1486.26329

±0.00016

1512.67421

±0.00017

1538.42321

±0.00016

2247.22051

±0.00018

2273.67262

±0.00019

1580.10933

±0.00007

b

0.535

±0.028

0.575

±0.021

0.591

±0.021

0.562

±0.027

0.571

±0.022

0.530

±0.031

0.566

±0.010

\( \frac{a}{{R_{star}}} \)

3.002

±0.034

2.819

±0.038

2.820

±0.039

2.945

±0.040

2.801

±0.039

2.959

±0.046

2.894

±0.016

\( u_{1}^{b} \)

0.357

±0.026

0.375

±0.028

0.360

±0.028

0.373

±0.025

0.389

±0.027

0.360

±0.026

0.369

±0.011

Derived Parameters

\( i(deg) \)

79.87

±0.54

78.41

±0.47

78.07

±0.47

79.16

±0.53

78.23

±0.49

79.75

±0.60

78.89

±0.20

a \( BJ{D_{TDB}}-2457000 \)

b Limb darkening factor given by [15]

We combined the transit timings we obtained from TESS data and previous transit-timing data from Ivshina and Winn [16] to study the trend of mid-transit time along with the orbit number. A portion of both data are provided in Table 2.

Table 2: A portion of transit data

Type of event

Midtime(BJDTDB)

Uncertainty(Days)

Orbit number

Source

tran

2454515.524960

0.000430

-1640

previous work

tran

2455951.835340

0.000110

-324

previous work

tran

2456305.455360

0.000240

0

previous work

tran

2456663.441360

0.000190

328

previous work

tran

2457788.694640

0.000480

1359

previous work

tran

2458852.827461

0.000624

2334

this work

tran

2459503.312844

0.000773

2930

this work

tran

2460284.766783

0.000770

3646

this work

NOTE: The entire data of transit and occultation midtimes are available from the authors in a machine readable format.

2.2. Occultation observation

In comparison to the transits (primary eclipses), the individual secondary eclipse events are not detectable by TESS with a sufficiently high signal-to-noise ratio for reliable timing measurements. Therefore, instead of fitting the data for each individual occultation within each sector, we combined all full and interrupted occultation events within a TESS sector to calculate a single effective mid-occultation time.

We selected all valid data within one referenced transit duration from the predicted mid-occultation time, and retained only those data ranges that contain enough data points (more than 80% of the median number across all such data ranges). We median-normalized each of the secondary eclipse events and phase-folded them by centering them at half a phase later than our reference transit time. To further reduce noise, we binned the time with a bin size of 10, indicating a 10-minute cadence.

For each sector, we fitted a trapezoid dip with fixed ingress and egress time and full transit duration, as best fitted and derived by Wong et al [17]. We derived the out-of-transit base level by taking the median of the data outside the duration time of the occultation and within the selected folded time.

We fitted for the mid-occultation time, \( {T_{occ}} \) . The orbit number of \( {T_{occ}} \) for each sector was set to the secondary events closest to the mid-time of the time series. The best-fit timings for the six TESS sectors were presented in section 3.

3. Data analysis

We have gathered all the transit and occultation data, with the addition of the new data from section 2. For the sectors that are previously mentioned in other literature, we have decided to use our derived mid-time for methodological consistency. Most of this data was compiled by Yee et al [18] with the addition of occultations from CHEOPS from Wong et al [19]. We have added 118 new transits and 6 new occultations into the catalog from sectors 20, 43, 44, 45, 71, and 72. We want to emphasize that the data in the table is based on the BJDTDB system, calibrated to the same reference mid-time.

We applied three models to fit the transit and occultation data mentioned in Patra et al[6]. The first is linear model, assuming a circular orbit and constant orbital period:

\( {t_{tran}}(E)={t_{0}}+EP \)

\( {t_{occ}}(E)={t_{0}}+\frac{P}{2}+EP \)

Table 3: Occultation events for WASP-12b from this work

Midtime(BJDTDB) Orbit Numberσ(Flux)

2458854.465105

2335

0.004020

2459487.475208

2915

0.003770

2459513.682400

2939

0.003989

2459537.690823

2961

0.003879

2460246.032689

3610

0.003943

2460273.300323

3635

0.004118

where E is the epoch from a selected reference orbit, and t0 is the mid-transit time of this reference orbit [6]

The second is quadratic model, also assuming a circular orbit, yet the orbital period to be decaying constantly with time:

\( {t_{tran}}(N)={t_{0}}+{NP_{s}}-\frac{{eP_{a}}}{π}cosω \)

\( {t_{occ}}(N)={t_{0}}+\frac{{P_{a}}}{2}+{NP_{s}}+\frac{{eP_{a}}}{π}cosω \)

\( ω(N)={ω_{0}}+\frac{dω}{dN}N \)

\( {P_{s}}={P_{a}}(1-\frac{1}{(2π)}\frac{dω}{dN}) \)

where Ps is the sidereal period, Pa is the anomalistic period, and ω is the argument of pericenter.

We used the emcee code to conduct MCMC sampling of the posterior distribution in parameter space to stabilize the fitting result, following the approach outlined in Yee et al [20]. We employed broad uniform priors on all parameters and ran the MCMC with 100 walkers. The initial 40% of steps were discarded as burn-in. This process was run for over 10 autocorrelation times. We also examined the posteriors and calculating Geweke scores for each chain [21] to ensured convergence. The fit results are presented in Table 4, where we minimized χ2 to fit the model in all three cases.

As expected from previous literature [6,19,20], a linear model does not fit the data well; here, the linear model has a minimum χ2 of 1348.2, with 251 degree of freedom. Compared to the linear model, the χ2 in quadratic model and sinusoidal model has much smaller discrepancy, each with a χ2 of 330.27 and 423.499. While ∆χ2 = 93.23 itself is expected, it’s significantly larger than previous literature[6, p. 20] (∆χ2 = 5.5[6], ∆χ2 = 12.1[20]). This can be contributed to the longer observation period and further boost our confidence of the potential orbital decay of WASP-12b.

To show the superiority of the quadratic model which includes less parameters than sinusoidal model, we use Bayesian Information Criterion (BIC), defined as followed:

\( BIC≡{χ^{2}}+klogn \)

where k is the quantity of free parameters and n is the quantity of data points. We can see a ∆BIC in this case is 104.10. In order to better illustrate where the distinction comes from, we plotted out the residual plot for transits and occultations separately (see Figure 1 , using the linear model as baseline. It can be difficult to determine whether the decaying trend is solely part of the long-term precession when only examining transits. However, upon analyzing the occlusions, we observe a clear deviation from the predicted model. This strengthens our confidence in the quadratic model.

4. Discussion about tidal orbital decay

Regarding the tidal orbital decay we mentioned earlier, we used the simplified “constant phase lag” model [2] that considers the equilibrium tide. Under the assumptions that the orbit is circular (eccentricity = 0), the planet’s mass is constant, and the orbital angular momentum is significantly greater than stellar angular momentum, we applied equation (20) from Goldreich & Soter :

\( \frac{da}{dt}=-\frac{9}{2}{(\frac{G}{{M_{⋆}}})^{(\frac{1}{2})}}\frac{R_{⋆}^{5}}{{Q_{⋆}} \prime }\frac{{M_{p}}}{{a^{(\frac{11}{2})}}} \)

where a is the orbital semi-major axis, and Q⋆ is the “reduced” tidal quality factor described as \( Q_{⋆}^{ \prime }=\frac{3}{2}\frac{{M_{p}}}{{M_{⋆}}}{(\frac{a}{{R_{⋆}}})^{5}}\frac{n}{|\dot{n}|} \)

where Q is the tidal quality factor and k2 is the Love number. Applying Kepler’s third law, we could get the decay rate of period:

\( {t_{tran}}(E)={t_{0}}+{EP_{s}}-{eP_{a}}πcosω \)

\( {t_{occ}}(E)={t_{0}}+\frac{{P_{a}}}{2}+{EP_{s}}+\frac{{eP_{a}}}{π}cosω \)

\( ω(E)={ω_{0}}+\frac{dω}{dE}E \)

\( {P_{s}}={P_{a}}1-\frac{1}{2π}\frac{dω}{dE} \)

/word/media/image1.jpeg

Figure 1: Combined O-C (Observed minus Calculated) diagram for transit and occultation timing deviations of WASP-12b. The top panel shows transit timing deviations, while the bottom panel displays occultation timing deviations. Gray circles represent previous observations, and green squares indicate new observations from this work. The blue line represents the best-fit quadratic model, and the red line represents the best-fit sinusoidal model. The dashed black line at y=0 represents the linear ephemeris. Error bars represent 1σ uncertainties

Using the measured decay rate of period of 29.5ms/yr, we obtain the reduced tidal quality factor Q⋆ of WASP-12 to be:

\( {Q_{star}}≈1.66±0.06×{10^{5}} \) (7)

We could then calculate the decay rate of orbital energy and the decay rate of angular momentum

\( \frac{dE}{dt}=\frac{{(2π)^{(\frac{2}{3})}}{M_{p}}}{3}{(\frac{{GM_{⋆}}}{P})^{(\frac{2}{3})}}\frac{1}{P}\frac{dP}{dt}=-5.0×{10^{23}}W \)

\( \frac{dL}{dt}=\frac{{M_{p}}}{3{(2π)^{\frac{1}{3}}}}{(\frac{{GM_{⋆}}}{P})^{\frac{2}{3}}}\frac{dP}{dt}=-7.5×{10^{27}}kg{m^{2}}{s^{-2}} \)

where \( {M_{star}}=1.434×0.110{M_{sun}} \) , and \( {M_{planet}}=1.434×0.110{M_{sun}} \) [22]. These results are similar to Yee et al[13] , and are obvious while a planet’s orbit is decaying for tidal interation.

However, the magnitude of Q⋆ is smaller than expected one in earlier studies[23]. Millholland & Laughlin[12] indicates this might result from a likely-existing planet orbiting the same host star and perturbing the orbit of WASP-12b. Future observations on WASP-12 are necessary to verify this hypothesis.

5. Summary

Our analysis of WASP-12b’s orbit, incorporating the latest TESS data from Sectors 20, 43, 44, 45, 71, and 72, along with previously published observations, strongly supports the quadratic model. The key findings of our study are:

We analyzed 118 new transit observations and 6 new occultations, significantly expanding the available dataset for WASP-12b. The quadratic model provided the best fit to the data, with a ∆BIC of 104.09 compared to the linear model, strongly favoring orbital decay over apsidal precession. We calculated an orbital period decay rate of \( \dot{P}=-29.5\frac{ms}{yr} \) , consistent with previous studies and indicating a uniform decay of WASP-12b’s orbit. The reduced tidal quality factor for the host star was determined to be \( {Q^{ \prime }_{⋆}}=1.66±0.06×{10^{5}} \) aligning with earlier research and supporting the orbital decay hypothesis. These results reinforce the uniqueness of the WASP-12 system, as it remains one of the clearest examples of rapid orbital decay among hot Jupiters. The consistency of our findings with previous work strengthens the case for ongoing tidal interactions between WASP-12b and its host star.

Table 4: Timing model fit parameters

Parameter

Value (Uncertainty)

Linear Model

Period, P (days)

1.091419000(0)

Mid-Transit Time of Reference Orbit, t0

2456305.455462(27)

Ndof

251

\( χ_{min}^{2} \)

48.20224

BIC

1359.269022

Quadratic Model

Period of Reference Orbit, P (days)

1.091420000(0)

Mid-Transit Time of Reference Orbit, t0

2456305.455803(29)

Decay Rate, dP/dE (days/orbit)

−10.21(32) × 1010

Ndof

250

\( χ_{min}^{2} \)

330.47445

BIC

347.074620

Sinusoidal Model

Sidereal Period, Ps (days)

1.091584000(7)

Mid-Transit Time of Reference Orbit, t0

2456305.454717(116)

Eccentricity, e

0.00478(27)

Argument of Periastron, ω0 (rad)

2.29(6)

Precession Rate, dω/dE (rad/orbit)

0.000949(40)

Ndof

248

\( χ_{min}^{2} \)

423.49943

BIC

451.166375

NOTE—Uncertainties in parentheses are the 1-σ confidence intervals in the last two digits.

The rapid orbital decay of WASP-12b continues to challenge our understanding of star-planet interactions and tidal dissipation mechanisms. Future observations and theoretical work will be crucial to fully explain the physical processes driving this decay and to predict the ultimate fate of the WASP-12 system. This research not only advances our knowledge of WASP-12b but also contributes to our broader understanding of hot Jupiter evolution and the dynamics of close-in exoplanetary systems.

6. Statement of work

This research was a collaborative effort among multiple contributors. Zoutong Shen was the primary developer, responsible for all coding, data cleaning, analysis, and modeling, as well as partial data collection. Shen also drafted sections on data analysis and collection and created most of the figures and tables in the main text. Hao Meng outlined the paper structure, authored the Abstract, parts of the Data Collection and Analysis sections, equations, Discussion, and Acknowledgement, and organized all data collection materials. Yishan Wu contributed to the literature review and wrote portions of the Abstract, Introduction, Discussion, and equations. Mei Zhang assisted with the literature review and composed part of the Introduction. Yunhe Fei contributed to writing the data analysis and data collection sections. All authors reviewed and approved the final manuscript.

Acknowledgments

We appreciate Professor Joshua N. Winn, who provided necessary instruction about the theory on the transit and occultation of WASP-12b, and python codes to deal with the data. The authors are very grateful to teaching assistant Xinyue Huang, who provided inspiring guidance and discussion about the technical skills of the codes and paper. We also gratefully acknowledge all the anonymous researchers in the Science Processing Operations Center (SPOC) for providing clear and essential light curve data of WASP-12. We thank CIS, which organized and arranged the summer camp for the authors to discuss and finish this paper. The convenient place for this summer camp is supported by Shanghai Jiao Tong University.


References

[1]. M. Mayor and D. Queloz, “A Jupiter-mass companion to a solar-type star,” Nature, vol. 378, no. 6555, pp. 355–359, Nov. 1995, doi: 10.1038/378355a0.

[2]. P. Goldreich and S. Soter, “Q in the Solar System,” Icarus, vol. 5, pp. 375–389, Jan. 1966, doi: 10.1016/0019-1035(66)90051-0.

[3]. D. N. C. Lin and J. Papaloizou, “On the Tidal Interaction between Protoplanets and the Protoplanetary Disk. III. Orbital Migration of Protoplanets,” Astrophys. J., vol. 309, p. 846, Oct. 1986, doi: 10.1086/164653.

[4]. “Hot Jupiters: Origins, Structure, Atmospheres - Fortney - 2021 - Journal of Geophysical Research: Planets - Wiley Online Library.” Accessed: Sep. 04, 2024. [Online]. Available: https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2020JE006629

[5]. G. Maciejewski et al., “Departure from the constant-period ephemeris for the transiting exoplanet WASP-12 b,” Astron. Astrophys., vol. 588, p. L6, Apr. 2016, doi: 10.1051/0004-6361/201628312.

[6]. K. C. Patra, J. N. Winn, M. J. Holman, L. Yu, D. Deming, and F. Dai, “The Apparently Decaying Orbit of WASP-12b,” Astron. J., vol. 154, no. 1, p. 4, Jun. 2017, doi: 10.3847/1538-3881/aa6d75.

[7]. L. Hebb et al., “WASP-12b: The Hottest Transiting Extrasolar Planet Yet Discovered,” Astrophys. J., vol. 693, pp. 1920–1928, Mar. 2009, doi: 10.1088/0004-637X/693/2/1920.

[8]. N. Husnoo et al., “Orbital eccentricity of WASP-12 and WASP-14 from new radial velocity monitoring with SOPHIE,” Mon. Not. R. Astron. Soc., vol. 413, pp. 2500–2508, Jun. 2011, doi: 10.1111/j.1365-2966.2011.18322.x.

[9]. K. A. Collins, J. F. Kielkopf, and K. G. Stassun, “TRANSIT TIMING VARIATION MEASUREMENTS OF WASP-12b AND QATAR-1b: NO EVIDENCE OF ADDITIONAL PLANETS,” Astron. J., vol. 153, no. 2, p. 78, Jan. 2017, doi: 10.3847/1538-3881/153/2/78.

[10]. A. Bailey and J. Goodman, “Understanding WASP-12b,” Mon. Not. R. Astron. Soc., vol. 482, pp. 1872–1882, Jan. 2019, doi: 10.1093/mnras/sty2805.

[11]. J.-P. Zahn, “Stellar Tides,” in Tides in Astronomy and Astrophysics, J. Souchay, S. Mathis, and T. Tokieda, Eds., Berlin, Heidelberg: Springer, 2013, pp. 301–325. doi: 10.1007/978-3-642-32961-6_8.

[12]. S. Millholland and G. Laughlin, “Obliquity Tides May Drive WASP-12b’s Rapid Orbital Decay,” Astrophys. J., vol. 869, p. L15, Dec. 2018, doi: 10.3847/2041-8213/aaedb1.

[13]. I. Wong et al., “TESS Phase Curve of the Hot Jupiter WASP-19b,” Astron. J., vol. 159, p. 104, Mar. 2020, doi: 10.3847/1538-3881/ab6d6e.

[14]. G. R. Ricker et al., “Transiting Exoplanet Survey Satellite (TESS),” J. Astron. Telesc. Instrum. Syst., vol. 1, p. 014003, Jan. 2015, doi: 10.1117/1.JATIS.1.1.014003.

[15]. K. Mandel and E. Agol, “Analytic Light Curves for Planetary Transit Searches,” Astrophys. J., vol. 580, pp. L171–L175, Dec. 2002, doi: 10.1086/345520.

[16]. E. S. Ivshina and J. N. Winn, “TESS Transit Timing of Hundreds of Hot Jupiters,” Astrophys. J. Suppl. Ser., vol. 259, no. 2, p. 62, Apr. 2022, doi: 10.3847/1538-4365/ac545b.

[17]. I. Wong et al., “TESS Revisits WASP-12: Updated Orbital Decay Rate and Constraints on Atmospheric Variability,” Astron. J., vol. 163, no. 4, p. 175, Apr. 2022, doi: 10.3847/1538-3881/ac5680.

[18]. S. W. Yee et al., “The Orbit of WASP-12b Is Decaying,” Astrophys. J. Lett., vol. 888, no. 1, p. L5, Dec. 2019, doi: 10.3847/2041-8213/ab5c16.

[19]. I. Wong et al., “The Hubble PanCET Program: A Featureless Transmission Spectrum for WASP-29b and Evidence of Enhanced Atmospheric Metallicity on WASP-80b,” Astron. J., vol. 164, p. 30, Jul. 2022, doi: 10.3847/1538-3881/ac7234.

[20]. S. W. Yee et al., “The Orbit of WASP-12b is Decaying,” Astrophys. J. Lett., vol. 888, no. 1, p. L5, Jan. 2020, doi: 10.3847/2041-8213/ab5c16.

[21]. A. Gelman and D. B. Rubin, “Inference from Iterative Simulation Using Multiple Sequences,” Stat. Sci., vol. 7, pp. 457–472, Jan. 1992, doi: 10.1214/ss/1177011136.

[22]. A. Chakrabarty and S. Sengupta, “Precise Photometric Transit Follow-up Observations of Five Close-in Exoplanets: Update on Their Physical Properties,” Astron. J., vol. 158, no. 1, p. 39, Jul. 2019, doi: 10.3847/1538-3881/ab24dd.

[23]. G. I. Ogilvie, “Tidal Dissipation in Stars and Giant Planets,” Annu. Rev. Astron. Astrophys., vol. 52, no. Volume 52, 2014, pp. 171–210, Aug. 2014, doi: 10.1146/annurev-astro-081913-035941.


Cite this article

Shen,Z.;Meng,H.;Wu,Y.;Zhang,M.;Fei,Y. (2025). Transit and Occultation Data Confirm a Decaying Orbit for WASP-12b. Theoretical and Natural Science,107,111-124.

Data availability

The datasets used and/or analyzed during the current study will be available from the authors upon reasonable request.

Disclaimer/Publisher's Note

The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of EWA Publishing and/or the editor(s). EWA Publishing and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

About volume

Volume title: Proceedings of the 4th International Conference on Computing Innovation and Applied Physics

ISBN:978-1-80590-087-0(Print) / 978-1-80590-088-7(Online)
Editor:Ömer Burak İSTANBULLU, Marwan Omar, Anil Fernando
Conference website: https://2025.confciap.org/
Conference date: 17 January 2025
Series: Theoretical and Natural Science
Volume number: Vol.107
ISSN:2753-8818(Print) / 2753-8826(Online)

© 2024 by the author(s). Licensee EWA Publishing, Oxford, UK. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license. Authors who publish this series agree to the following terms:
1. Authors retain copyright and grant the series right of first publication with the work simultaneously licensed under a Creative Commons Attribution License that allows others to share the work with an acknowledgment of the work's authorship and initial publication in this series.
2. Authors are able to enter into separate, additional contractual arrangements for the non-exclusive distribution of the series's published version of the work (e.g., post it to an institutional repository or publish it in a book), with an acknowledgment of its initial publication in this series.
3. Authors are permitted and encouraged to post their work online (e.g., in institutional repositories or on their website) prior to and during the submission process, as it can lead to productive exchanges, as well as earlier and greater citation of published work (See Open access policy for details).

References

[1]. M. Mayor and D. Queloz, “A Jupiter-mass companion to a solar-type star,” Nature, vol. 378, no. 6555, pp. 355–359, Nov. 1995, doi: 10.1038/378355a0.

[2]. P. Goldreich and S. Soter, “Q in the Solar System,” Icarus, vol. 5, pp. 375–389, Jan. 1966, doi: 10.1016/0019-1035(66)90051-0.

[3]. D. N. C. Lin and J. Papaloizou, “On the Tidal Interaction between Protoplanets and the Protoplanetary Disk. III. Orbital Migration of Protoplanets,” Astrophys. J., vol. 309, p. 846, Oct. 1986, doi: 10.1086/164653.

[4]. “Hot Jupiters: Origins, Structure, Atmospheres - Fortney - 2021 - Journal of Geophysical Research: Planets - Wiley Online Library.” Accessed: Sep. 04, 2024. [Online]. Available: https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2020JE006629

[5]. G. Maciejewski et al., “Departure from the constant-period ephemeris for the transiting exoplanet WASP-12 b,” Astron. Astrophys., vol. 588, p. L6, Apr. 2016, doi: 10.1051/0004-6361/201628312.

[6]. K. C. Patra, J. N. Winn, M. J. Holman, L. Yu, D. Deming, and F. Dai, “The Apparently Decaying Orbit of WASP-12b,” Astron. J., vol. 154, no. 1, p. 4, Jun. 2017, doi: 10.3847/1538-3881/aa6d75.

[7]. L. Hebb et al., “WASP-12b: The Hottest Transiting Extrasolar Planet Yet Discovered,” Astrophys. J., vol. 693, pp. 1920–1928, Mar. 2009, doi: 10.1088/0004-637X/693/2/1920.

[8]. N. Husnoo et al., “Orbital eccentricity of WASP-12 and WASP-14 from new radial velocity monitoring with SOPHIE,” Mon. Not. R. Astron. Soc., vol. 413, pp. 2500–2508, Jun. 2011, doi: 10.1111/j.1365-2966.2011.18322.x.

[9]. K. A. Collins, J. F. Kielkopf, and K. G. Stassun, “TRANSIT TIMING VARIATION MEASUREMENTS OF WASP-12b AND QATAR-1b: NO EVIDENCE OF ADDITIONAL PLANETS,” Astron. J., vol. 153, no. 2, p. 78, Jan. 2017, doi: 10.3847/1538-3881/153/2/78.

[10]. A. Bailey and J. Goodman, “Understanding WASP-12b,” Mon. Not. R. Astron. Soc., vol. 482, pp. 1872–1882, Jan. 2019, doi: 10.1093/mnras/sty2805.

[11]. J.-P. Zahn, “Stellar Tides,” in Tides in Astronomy and Astrophysics, J. Souchay, S. Mathis, and T. Tokieda, Eds., Berlin, Heidelberg: Springer, 2013, pp. 301–325. doi: 10.1007/978-3-642-32961-6_8.

[12]. S. Millholland and G. Laughlin, “Obliquity Tides May Drive WASP-12b’s Rapid Orbital Decay,” Astrophys. J., vol. 869, p. L15, Dec. 2018, doi: 10.3847/2041-8213/aaedb1.

[13]. I. Wong et al., “TESS Phase Curve of the Hot Jupiter WASP-19b,” Astron. J., vol. 159, p. 104, Mar. 2020, doi: 10.3847/1538-3881/ab6d6e.

[14]. G. R. Ricker et al., “Transiting Exoplanet Survey Satellite (TESS),” J. Astron. Telesc. Instrum. Syst., vol. 1, p. 014003, Jan. 2015, doi: 10.1117/1.JATIS.1.1.014003.

[15]. K. Mandel and E. Agol, “Analytic Light Curves for Planetary Transit Searches,” Astrophys. J., vol. 580, pp. L171–L175, Dec. 2002, doi: 10.1086/345520.

[16]. E. S. Ivshina and J. N. Winn, “TESS Transit Timing of Hundreds of Hot Jupiters,” Astrophys. J. Suppl. Ser., vol. 259, no. 2, p. 62, Apr. 2022, doi: 10.3847/1538-4365/ac545b.

[17]. I. Wong et al., “TESS Revisits WASP-12: Updated Orbital Decay Rate and Constraints on Atmospheric Variability,” Astron. J., vol. 163, no. 4, p. 175, Apr. 2022, doi: 10.3847/1538-3881/ac5680.

[18]. S. W. Yee et al., “The Orbit of WASP-12b Is Decaying,” Astrophys. J. Lett., vol. 888, no. 1, p. L5, Dec. 2019, doi: 10.3847/2041-8213/ab5c16.

[19]. I. Wong et al., “The Hubble PanCET Program: A Featureless Transmission Spectrum for WASP-29b and Evidence of Enhanced Atmospheric Metallicity on WASP-80b,” Astron. J., vol. 164, p. 30, Jul. 2022, doi: 10.3847/1538-3881/ac7234.

[20]. S. W. Yee et al., “The Orbit of WASP-12b is Decaying,” Astrophys. J. Lett., vol. 888, no. 1, p. L5, Jan. 2020, doi: 10.3847/2041-8213/ab5c16.

[21]. A. Gelman and D. B. Rubin, “Inference from Iterative Simulation Using Multiple Sequences,” Stat. Sci., vol. 7, pp. 457–472, Jan. 1992, doi: 10.1214/ss/1177011136.

[22]. A. Chakrabarty and S. Sengupta, “Precise Photometric Transit Follow-up Observations of Five Close-in Exoplanets: Update on Their Physical Properties,” Astron. J., vol. 158, no. 1, p. 39, Jul. 2019, doi: 10.3847/1538-3881/ab24dd.

[23]. G. I. Ogilvie, “Tidal Dissipation in Stars and Giant Planets,” Annu. Rev. Astron. Astrophys., vol. 52, no. Volume 52, 2014, pp. 171–210, Aug. 2014, doi: 10.1146/annurev-astro-081913-035941.