{"id":"175ee915-913a-4e4f-a07a-de18bdb9e272","arxiv_id":"2412.05609","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":13,"one_line_summary":"TOI-4504 c shows the largest known transit timing variations, revealing a non-transiting Jovian companion in a 2:1 mean-motion resonance.","lead":"Astronomers found a Jupiter-sized planet whose transits arrive up to two days early or late, the largest timing wobble ever seen, and used it to uncover a hidden second giant planet. The pair appears locked in a 2:1 orbital resonance, a configuration that hints at how giant planets migrate inward.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 41.2 d RV signal for d crosses only the 10% FAP threshold, and Section 4.3 does not show the TTV data alone reject the alternative modes in Fig. A1; the d mass/resonance claims ride on that marginal peak and a 103 m/s jitter comparable to K_d.","rationale":"The paper is solid on the headline TTV detection: the transit times in Table 3 show a secure, roughly 2-day peak-to-node modulation, and the comparison with known TTV systems supports the record-amplitude claim. The stellar parameters, contamination checks, and stability integrations are all reasonable. The weakest point is the identification and characterization of TOI-4504 d. The reader flagged the marginal 41.2 d RV signal and jitter, and I agree that this is the critical assumption. I add that the paper does not demonstrate that the TTVs alone uniquely select the 40.56-day interior 2:1 solution over the multimodal alternatives admitted in Fig. A1, and the quoted ~930 d super-period is not transparently reproduced by Eq. (1) from the Table 4 osculating periods. A TTV-only fit with evidence comparison would settle whether the d mass and resonance geometry are independently constrained or rely on the borderline RV peak. Since the central discovery of a giant perturber is secure but the specific mass-period-resonance solution is not fully secured, the CONDITIONAL verdict should stand unchanged.","tokens_in":27192,"tokens_out":18301,"duration_ms":174985,"concrete_test":"Re-run the Section 4.3 global Nested Sampling search using the TESS TTVs only, dropping the FEROS RVs entirely, with the same broad priors as in Table A1. Compute the Bayesian evidence difference between the best interior 40.56-day 2:1 mode and the next-best modes visible in Fig. A1. If the TTV-only posterior does not recover the same 40.56 d / 1.42 MJ solution with Δln Z > 5 over the alternatives, the d mass and 2:1 resonance claims depend on the marginal 41.2 d RV signal and should be presented as provisional; if it does, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Sect. 4.2 the post-84 d GLS peak at 41.2 d is reported to cross only the 10% FAP level, and the final N-body RV residuals show no significant peaks. Table 4 lists K_d = 90.8 ± 2.5 m/s with an adopted RV jitter of 103.4 m/s, i.e. the jitter is comparable to the signal. The TTVs prove that a large perturbing companion exists, but the paper's own Fig. A1 shows the joint TTV+RV search is multimodal, with alternative period ratios capable of matching the data; the text asserts that the 40.56 d 2:1 solution is favored but reports no TTV-only fit and no Δln Z or ΔBIC separating it from the other modes. If the 41 d periodicity is partly or wholly activity or an unmodeled harmonic of the 84 d signal, the RV term can shift K_d and the derived mass, and the 'likely locked in 2:1 MMR' conclusion is only as strong as the uniqueness of the 40.56 d solution. This is the least secure link because the TTV amplitude mainly constrains a mass/Δ combination, not independently the period and eccentricity of d.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports a joint TTV and radial-velocity analysis of the TOI-4504 system, confirming a transiting warm Jupiter (TOI-4504 c, P ≈ 82.54 d, m ≈ 3.77 MJup) and an inner non-transiting Jovian planet (TOI-4504 d, P ≈ 40.56 d, m ≈ 1.42 MJup) claimed to be in a 2:1 mean-motion resonance. The TTVs of TOI-4504 c have a peak-to-node amplitude of about 2 days and a super-period near 930–950 d, making this the largest absolute TTV signal reported to date. The paper also characterizes an inner sub-Neptune (TOI-4504 b) with P ≈ 2.426 d and Rp ≈ 2.69 R⊕. The detection of the large-amplitude TTV is secure, but the specific period and mass of planet d rest on a multimodal joint fit and on an RV signal that only crosses the 10% false-alarm-probability threshold, with an adopted RV jitter comparable to the claimed RV semi-amplitude of d.","tokens_in":27587,"tokens_out":5189,"duration_ms":49606,"significance":"If the inferred architecture is correct, TOI-4504 provides a new benchmark for TTV amplitudes and a further example of a warm-Jupiter pair in a first-order mean-motion resonance, which is relevant to disk-migration theories. The paper's strengths include the secure detection of very large TTVs, the joint TTV+RV N-body modeling with publicly available tools, the explicit presentation of predicted future transit times, and a careful stellar characterization. However, the central claim about the period, mass, and resonant state of TOI-4504 d is not yet uniquely established, and the paper's own text acknowledges that the large unexplained RV scatter could not be fully modeled.","major_comments":[{"comment":"The global nested-sampling search for the period of TOI-4504 d is reported as multimodal, but the paper does not provide a quantitative comparison between the preferred 40.56 d, 2:1 solution and the other modes. The text says the 41-day solution leads to \"significantly better fits,\" but no Δln Z, ΔBIC, or equivalent evidence is given. Because the claimed mass and period of d, and hence the resonance interpretation, depend on the uniqueness of this solution, please report either the Bayes factor or BIC for the 40.56 d solution relative to the other modes in Fig. A1. In addition, a fit using only the TTV data (without RVs) should be shown to demonstrate whether the 40.56 d period is identified by the TTVs alone; otherwise the RV data, which are marginal, may be driving the period selection.","section":"§4.3 and Fig. A1"},{"comment":"The 41.2 d RV peak crosses only the 10% FAP threshold, and the fitted RV jitter is 103.4 m/s, which is comparable to K_d = 90.8 m/s. The paper itself states that the large RV scatter \"cannot be fully resolved\" and that the source of the jitter remains unclear. This makes the RV contribution to the d mass and period uncertain. Please quantify how the derived K_d, mass, and period of d change if the RVs are excluded from the fit or if the jitter is modeled more flexibly (e.g., with a Gaussian process or an activity-correlated component). If the mass is not robust to these choices, the text should temper the claims about the mass and period of planet d.","section":"§4.2 and Table 4"},{"comment":"The abstract and summary describe TOI-4504 c as having \"the largest TTV amplitude ever observed.\" The paper's own Fig. 10 shows that this statement holds for absolute peak-to-node amplitude, but not for the relative TTV amplitude (TTV amplitude divided by orbital period). Please clarify this distinction explicitly in the abstract and main text so that the record claim is not overstated.","section":"§1 and §5"}],"minor_comments":[{"comment":"The super-period of the TTVs is quoted inconsistently: the abstract says ~930 d, the Fig. 6 caption says 946.5 d, and §5 says ~2.9 years (about 1060 d). Please harmonize these values.","section":"Abstract, Fig. 6 caption, §5"},{"comment":"The text first states that after subtracting the 84 d signal \"no other significant signals were detected,\" but then refers to the 41.2 d peak as \"prominent\" despite it crossing only the 10% FAP threshold. This wording is contradictory; please specify that 10% FAP is not a significant detection threshold in the usual sense.","section":"§4.2"},{"comment":"The prior for the inclination of planet c is listed as U(89.7, 0.1), which appears to be a typo (likely N(89.7, 0.1), as in Table 4, or a misprinted uniform range). Please correct it.","section":"Table A1"},{"comment":"The derived eccentricity uncertainties (e.g., e_d = 0.0445 ± 0.0010) are quoted to three decimal places; given the 103 m/s jitter and the marginal RV detection of d, please verify that these uncertainties are not underestimated by the adopted white-noise jitter model.","section":"Table 4"},{"comment":"The eight non-detections of predicted transits of c are attributed to large TTVs. It would be helpful to explicitly show that these non-detections are consistent with the best-fit TTV model, for example by plotting the expected transit times against the observed windows; otherwise the non-detections could be interpreted as tension.","section":"§2.4"}],"recommendation":"major_revision","confidential_remarks":"The paper is timely and addresses a striking new TTV system, but the central claim regarding TOI-4504 d's period and mass is only as strong as the uniqueness of the 40.56 d solution. The authors should be encouraged to provide the requested model comparisons and TTV-only fits; without them, the resonance claim is not yet fully supported. The paper is otherwise well within the scope of the journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line up front: this is a genuine discovery paper. TOI-4504 c shows the largest TTV amplitude on record, the signal is huge relative to the measurement uncertainties, and there is definitely something out there pulling on it. The soft spot is the identity of that something: the 41-day period and 1.4 Jupiter mass for planet d come from a joint TTV+RV fit where the RV component is marginal and the parameter space is multimodal.\n\nWhat the paper does well: clean TESS TTV extraction across many sectors, sensible two-step modeling with Exo-Striker, careful stellar characterization, speckle imaging to rule out blends, and a stability analysis that draws from the MCMC posteriors rather than just the best fit. The context figure comparing TTV amplitudes with known systems is useful. The paper is also honest that the 103 m/s RV jitter is not fully explained and that more precise RVs are needed.\n\nThe weak spots, in order of how much they worry me. First, the 41.2 d RV signal only crosses the 10% false-alarm threshold, and the fitted jitter is comparable to K_d. The paper says the RVs 'conclusively confirm' two resonant planets, which is stronger than the evidence. The TTVs confirm a perturbing planet; they do not by themselves pin down its period at 40.56 d. Second, Figure A1 shows the global search is multimodal, but the paper gives no TTV-only fit and no Δln Z or ΔBIC separation between the 40.56 d mode and the alternatives. If an alternative mode is only a fraction of a log-likelihood away, the 2:1 resonance conclusion weakens. Third, the resonance libration angle is computed from the fitted parameters, so it is a property of the model, not an independent prediction. These are addressable with a bit more analysis, not fatal flaws.\n\nWho should read this: anyone working on TTV inversion, resonant giant planet formation, or warm Jupiter demographics. The record amplitude alone makes it a reference point for future TTV searches. I would send it to peer review. A good referee should ask for the TTV-only posterior, a quantitative model comparison between the competing modes, and a discussion of how activity-correlated noise could bias K_d. With those additions the paper will be solid.\n\nRecommendation: send out for review, require those revisions.","headline":"A genuine TTV record with a secure perturber, but the resonant architecture rests on a 41 d signal that the RVs alone barely support.","tokens_in":28252,"tokens_out":2269,"would_cite":true,"duration_ms":22373,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"TOI-4504 c shows the largest transit timing variations ever observed—about two days of wobble—and the wobble is caused by a hidden Jupiter-mass planet in a 2:1 resonance.","keywords":["transit timing variations","mean-motion resonance","warm Jupiter","non-transiting planet","exoplanet dynamics","TOI-4504","radial velocity"],"falsifier":"A decisive test is to obtain high-precision radial-velocity series around the predicted 40.56-day signal of TOI-4504 d: if the 41.2-day periodogram peak does not persist above the 1% false-alarm level once the ~103 m/s jitter is modeled together with stellar activity indicators, the planetary interpretation of the TTV signal is weakened; conversely, a clean detection at ~91 m/s semi-amplitude would confirm the hidden planet.","tokens_in":27002,"feed_emoji":"🪐","tokens_out":6953,"duration_ms":56438,"temperature":0.7,"pith_summary":"This paper argues that the transiting planet TOI-4504 c is the most extreme known case of transit timing variations, with its transits arriving up to roughly two days early or late in a regular ~930-day cycle. Using those timing variations together with radial-velocity measurements, the authors uncover a previously unknown non-transiting planet, TOI-4504 d, that is causing the wobble. The system contains three planets: a hot sub-Neptune and two warm Jupiter-mass gas giants locked in a 2:1 mean-motion resonance. The result matters because strongly resonant giant-planet pairs are direct evidence that these planets migrated smoothly inward through the protoplanetary disk.","feed_headline":"Record 2-day transit wobbles reveal a hidden giant planet","feed_subtitle":"The pair orbits in a 2:1 resonance, a key test for how giant planets migrate.","key_machinery":"The central object is the transit timing variation signal of TOI-4504 c: a roughly sinusoidal modulation of transit arrival times with ~930-day super-period and ~2-day peak-to-node amplitude, produced by the gravitational tug of the inner non-transiting planet d. The argument is carried by a Jacobi-coordinate N-body dynamical model that simultaneously fits space-based transit times and ground-based radial velocities, using a fast analytic transit-timing prescription for the TTV computation and a resonance-angle diagnostic ($\\theta_1 = \\lambda_c - 2\\lambda_d + \\omega_c$); the libration of $\\theta_1$ around $0^\\circ$ is the key signature that the two Jupiters are genuinely trapped in a 2:1 mean-motion resonance.","core_discovery":"The paper claims that the warm Jupiter TOI-4504 c, with osculating period $P = 82.54 \\pm 0.02$ days and dynamical mass $m = 3.77 \\pm 0.18\\,M_{\\rm J}$, shows transit timing variations with peak-to-node semi-amplitude of about 2 days and a super-period of about 930 days—the largest TTV amplitude ever recorded. The joint TTV plus radial-velocity dynamical model identifies the perturber as an interior non-transiting planet TOI-4504 d, with period $P = 40.56 \\pm 0.04$ days and dynamical mass $m = 1.42^{+0.07}_{-0.06}\\,M_{\\rm J}$, placing the pair in a first-order 2:1 mean-motion resonance. The paper also reports a transiting sub-Neptune TOI-4504 b with period $2.42614$ days and radius $2.69 \\pm 0.19\\,R_{\\oplus}$, whose mass is unconstrained by the data. Long-term N-body integrations of the giant pair show stability over 10 million years, with the first-order resonance angle $\\theta_1$ librating around $0^\\circ$, which the authors take as evidence of an active 2:1 resonance rather than a mere near-commensurability.","pith_inferences":["If the large libration amplitude (~65–73°) persists over many super-periods, the TOI-4504 pair could become a probe of resonant libration damping, which is expected to slowly shrink the libration amplitude over time.","The unmodeled ~100 m/s radial-velocity jitter is comparable in size to planet d's own signal, so a plausible extension is that additional planets or persistent stellar activity hide in the data; higher-precision velocities would either sharpen d's mass or reveal a fourth body.","Because the model implies planet d is more likely to transit than c, a dedicated transit search of the inner companion could detect it directly, and its TTV signal would be more than half again as large as the one observed for c.","Applying the same joint timing-and-velocity analysis to other 2:1 resonant warm-Jupiter pairs could turn the growing sample into a population-level test of disk-driven migration rates."],"forward_implications":["If the dynamical solution is correct, TOI-4504 becomes the benchmark for the largest observed transit-timing amplitude, roughly twice the previous record holder Kepler-30 b.","The pair of warm Jupiters in a 2:1 resonance adds to the small sample of resonant giant-planet systems whose growing occurrence supports smooth, disk-driven migration into resonance.","Because the perturber d is non-transiting, an observer whose line of sight made d transit would see TTVs more than 50% larger than those seen for c, making the resonance an even stronger timing signal from other viewing geometries.","The predicted transit times over the next decade give a direct observational target: future photometry near those windows can test the model's forecasts."],"supporting_citations":[{"why":"Establishes that planets in or near low-order mean-motion resonances exhibit the largest transit timing variations, which sets up the paper's expectation for a large TTV signal.","marker":"Agol et al. 2005"},{"why":"Provides the fast analytic transit-timing model used inside the N-body fit to compute the TTVs of TOI-4504 c.","marker":"Deck et al. 2014"},{"why":"Supplies the joint TTV-plus-radial-velocity N-body fitting scheme, applied here to a similar 2:1 resonant warm-giant system.","marker":"Trifonov et al. 2021"},{"why":"Supplies the Kepler transit-timing catalogue used to place TOI-4504 c's amplitude in context and to claim it is the largest known.","marker":"Holczer et al. 2016"},{"why":"Quantifies the previous record-holder Kepler-30 b's TTV semi-amplitude, which TOI-4504 c roughly doubles.","marker":"Panichi et al. 2018"},{"why":"Defines the resonance-angle libration criterion used to conclude that the pair is caught in an active 2:1 mean-motion resonance.","marker":"Lee 2004"}],"fun_headline_variants":["Record 2-day TTVs reveal a hidden giant planet","Largest transit wobble yet uncovers a resonant pair","2-day timing shifts expose an unseen giant","Hidden planet discovered via 2-day TTV super-period","Warm Jupiter's record wobble betrays a second giant"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The 41.2-day radial-velocity signal attributed to TOI-4504 d is only marginally significant, crossing just the 10% false-alarm threshold, and the fitted radial-velocity jitter of 103 m/s is as large as the 91 m/s semi-amplitude of planet d itself; the paper's masses and periods for d therefore depend on the assumption that this unexplained jitter does not bias the two-planet solution.","fun_headline_variants_meta":{"raw":{"variants":["Record 2-day TTVs reveal a hidden giant planet","Largest transit wobble yet uncovers a resonant pair","2-day timing shifts expose an unseen giant","Hidden planet discovered via 2-day TTV super-period","Warm Jupiter's record wobble betrays a second giant"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000693,"raw_usage":{"total_tokens":3302,"prompt_tokens":1280,"completion_tokens":2022,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":896,"completion_tokens_details":{"reasoning_tokens":1942}},"tokens_in":896,"tokens_out":2022,"duration_ms":15518,"temperature":1.0,"reasoning_tokens":1942,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T20:32:00.379616+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive test is to obtain high-precision radial-velocity series around the predicted 40.56-day signal of TOI-4504 d: if the 41.2-day periodogram peak does not persist above the 1% false-alarm level once the ~103 m/s jitter is modeled together with stellar activity indicators, the planetary interpretation of the TTV signal is weakened; conversely, a clean detection at ~91 m/s semi-amplitude would confirm the hidden planet.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the resonance-angle libration criterion used to conclude that the pair is caught in an active 2:1 mean-motion resonance."}],"review_version":1}