Transit Timing of the White Dwarf-Cold Jupiter System WD 1856+534

Abstract: We present new transit timing measurements for the white dwarf-cold Jupiter system WD 1856+534, extending the baseline of observations from 311 epochs to 1498 epochs. The planet is unlikely to have survived the host star's red-giant phase at its present location and is likely too small for common-envelope evolution to take place. As such, a plausible explanation for the short semimajor axis is that the exoplanet started out on a much larger orbit and then spiraled inward through high-eccentricity tidal migration (HETM). A past transit-timing analysis found tentative evidence for orbital growth, which could have been interpreted as a residual effect of HETM, but we find the data are consistent with a constant-period model after adding 18 new transit measurements. We use the estimated period derivative $\dot{P} = 0.04\pm0.43$ ms yr$^{-1}$ to place a lower limit on the planetary tidal quality factor of $Q_p' \gtrsim 3.1 \times 10^6$, consistent with that of Jupiter in our own Solar System. We also test for the presence of companion planets in the system, which could have excited WD 1856 b onto an eccentric orbit via the Kozai-Lidov process, and ultimately rule out the presence of an additional planet with a mass greater than $4.1\,M_J$ and a period shorter than 1500 days. We find no evidence for nonzero eccentricity, with an upper limit of $e\lesssim10^{-2}$. If the planet indeed reached its current orbit through HETM, the low present-day eccentricity indicates that the migration process has now ceased, and any further orbital evolution will be governed solely by weak planetary tides.

• The paper finds no significant orbital decay for WD 1856+534b, ruling out ongoing major tidal evolution.
• Transit timing analysis using extended observations and MCMC fitting provided precise constraints on the orbital period and an eccentricity of less than 0.01.
• No evidence for companion planets was detected, limiting potential perturbers involved in the planet's high-eccentricity tidal migration.

Transit Timing Constraints for the White Dwarf–Cold Jupiter System WD 1856+534

Introduction

This work presents a comprehensive analysis of transit timing variations (TTVs) for WD 1856+534b, a Jupiter-sized planet on a short-period orbit around a degenerate white dwarf. The system is of utmost interest to the exoplanet community because WD 1856+534b occupies a `forbidden zone’: it orbits at only 0.0204 AU—well within the distance where it should have been engulfed during the host’s red-giant phase, according to stellar evolution models. The formation and present configuration of this system thus provide empirical constraints on high-eccentricity tidal migration (HETM), common-envelope (CE) evolution, and the viability of planetary companionship and migration mechanisms.

Observational Methods and Data Reduction

New transit light curves for WD 1856+534b were obtained with the Lick Nickel 1-m telescope, adding 20 transit epochs to the observational baseline, which now comprises 68 transits over 1498 orbital cycles. Data reduction was performed with AstroImageJ, and light curve modeling utilized the Mandel & Agol formalism with iterative treatment of limb-darkening parameters due to the grazing geometry and sampling limitations of individual transits. The stacking of light curves enabled improved constraints on geometric and limb-darkening coefficients.

Figure 1: Stacked light curve of WD 1856+534b from the extended dataset, demonstrating deep ($\sim$56%) grazing transits typical for white dwarf–planet systems.

Transit Timing Modeling and Orbital Change Constraints

Transit times were fit with both constant-period and quadratic (period derivative) models using MCMC sampling. The extension of the baseline allowed for highly sensitive detection of secular changes in the orbital period:

• Reference period: $P = 1.407939211 \pm 5\times10^{-9}$ days.
• Period derivative: $\dot{P} = 0.04 \pm 0.43$ ms yr${}^{-1}$; consistent with zero over the baseline.

No statistically significant evidence for orbital growth or decay was found, robustly ruling out ongoing major tidal evolution at the currently probed sensitivity.

Figure 2: New transit light curves and fitted models from Nickel telescope data demonstrate the stability and reproducibility of observed WD 1856+534b transits.

Tidal Quality Factor Limits and Dynamical Implications

The period derivative constraint was translated into lower limits for the tidal quality factors:

• Stellar tidal quality factor: $Q_*' \gtrsim 0.0034$
• Planetary tidal quality factor: $Q_p' \gtrsim 3.1 \times 10^6$

Given the extremely compact nature of the white dwarf, it is expected that planetary tides dominate any residual orbital evolution, consistent with Jupiter's $Q_p'$ in the Solar System.

Spin-orbit synchronization considerations, using planetary rotation rates and $Q_p'$ estimates, indicate that unless $Q_p'$ is anomalously high ($\gtrsim 10^9$), WD 1856+534b should already be tidally locked, and further tidal evolution is effectively halted.

Constraints on Companion Planets

TTVs were analyzed using grid searches and 7D MCMC sampling incorporating the TTVFast code. Phase space was explored over a range of possible companion planet masses ($1\,M_\oplus$ to $13\,M_J$) and periods (50–1500 days):

• No significant likelihood improvement for any companion; the best grid-search suggestion ($0.94\,M_J$, 99.5 days) produced a $\chi^2$ reduction insufficient to justify a companion.
• BIC penalizes complexity, with the constant-period model preferred over all alternative configurations with companion planets of $>4.1 \, M_J$ and periods $< 1500$ days.

  Figure 3: Posterior distributions from TTVFast-based MCMC analysis; higher probability density is observed at low companion mass and large period, where sensitivity is sharply reduced.

  

  Figure 4: Likelihood ratio log map for hypothetical companion planets; extensive regions of phase space (purple) are ruled out by transit data, while no viable region leads to a statistically justified companion.

  

  Figure 5: Zoomed-in likelihood ratio map confirms exclusion of companions down to $\sim1\,M_J$ at short periods, with no compelling evidence for closer planetary perturbers.

Limits on Eccentricity and Orbital Evolution

The non-detection of coherent TTVs constrains the orbital eccentricity:

• $e \lesssim 0.01$, as higher eccentricity would produce precession-induced TTVs well above the detectability threshold ($\gtrsim5$ s) over the baseline.
• MCMC fits incorporating apsidal precession confirm $e=0$ is strongly preferred.

Residual energy dissipation via tides ($dE/dt$) with current upper limits cannot power significant luminosity in the planet, and the observed temperature is not inconsistent with equilibrium irradiation.

Figure 6: O–C (observed minus calculated) timing residuals for the constant-period and optimal companion models; absence of significant structure excludes nonzero eccentricity and companion-induced periodic TTVs.

Implications for Formation Scenarios

CE evolution scenarios require an unfeasibly high energy transfer to eject the envelope, especially given the $\sim$5–6 $M_J$ measured mass for WD 1856b. This renders CE models untenable unless extreme assumptions about efficiency ($\alpha_{\rm CE}$) are made. HETM remains the most plausible scenario, requiring a source of strong eccentricity excitation—most likely a now-absent planetary, stellar, or flyby companion.

The current absence of detectable companions and negligible eccentricity implies that migration via HETM has ceased, leaving the planet on a nearly circular orbit governed only by weak planetary tides.

Conclusions

Extended transit timing observations of WD 1856+534b have decisively constrained the system’s present-day dynamical architecture:

• There is no evidence for orbital decay or growth, with period changes $<0.9$ ms yr${}^{-1}$ ($2\sigma$).
• No companion planet $>4.1\,M_J$ with period $<1500$ days exists, limiting classes of possible past perturbers required for HETM.
• Planetary tidal quality factors are consistent with those of Solar System gas giants.
• Eccentricity is $<0.01$, ruling out active dynamical migration and tidal heating.

The theoretical and observational constraints presented here reinforce the complexity of post-main-sequence planetary system evolution, and the WD 1856+534 system remains a touchstone for future theoretical modeling, further IR and high-cadence photometric monitoring, and the search for more white dwarf–planetary systems.