Photodynamical analysis of the TOI-4504 system: Orbital precession brings planet d into transit revealing its deeply resonant state and similarity to GJ 876

Posterior samples Posterior samples from the photodynamical modeling.Filename: TOI-4504_Almenara2026_posterior.csvColumns: G: N-body gravitational constant (au^3 day^-2 M_Sun^-1) Tref: Reference time for the orbital elements (BJD_TDB) timestep: N-body time step (days) coordinate: Coordinate system m_s: Star mass (M_Sun) r_s: Star radius (au) m_p: Planet mass (M_Sun) r_p: Planet radius (au) a_p: Planet's semimajor axis (au) e_p: Planet's eccentricity i_p: Planet's orbital inclination (degrees) w_p: Planet's argument of pericentre (degrees) Node_p: Planet's longitude of the ascending node (degrees) M_p: Planet's mean anomaly (degrees) 'p' is the planet number, ranging from 1 (the innermost planet) to 3 (the outermost planet): 1: TOI-4504 b 2: TOI-4504 d 3: TOI-4504 c Transit forecast Transit timing predictions derived from the photodynamical modeling.Filename: TOI-4504_Almenara2026_forecast2030.csvColumns: 1. Planet number (see the description of the posterior file) 2. Epoch number of the transit 3. Transit time (BJD_TDB) 4. Lower 68.3% credible interval (days) 5. Upper 68.3% credible interval (days) Light curves Transit observations of TOI-4504. Filename: TOI-4504_[Mid-transit UT date]_[Telescope]_[Band].txtColumns: 1. Mid-Exposure Time (BJD_TDB) 2. Normalized Flux 3. Normalized Flux Error 4. Exposure Time (seconds) Additional columns : point spread function centroid shifts (pixels), full width at half maximum (pixels), peak flux value (ADU), and airmass