{"id":"b7e6c849-72ec-4e36-8a4b-8837649aed6a","arxiv_id":"2505.11637","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Stellar spin-down can knock outer companions to hot Jupiters out of the transiting plane, biasing transit surveys against them, while inner companions remain detectable.","lead":"This paper models how the spin-down of a host star can tilt nearby outer companions of hot Jupiters out of the transiting plane, while inner companions stay put. It maps which orbital configurations hide these outer planets and shows that higher stellar obliquity makes them even harder to detect.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The decoupling bias claim requires planets to be in place before the star's J2 crossed the secular resonance; the paper states this in §5.1 but never tests whether the observed hot Jupiter systems satisfy it, so the headline observational conclusion is conditional on early assembly.","rationale":"Read in good faith: the paper makes a narrow dynamical claim—inner companions remain in the hot Jupiter's transiting plane, outer companions can be decoupled by the decaying stellar quadrupole, and obliquity amplifies this. The analytical criterion (Eq. 6) and the N-body results are internally consistent, and the qualitative asymmetry is supported by prior work the authors cite. The weakest link is not the numerics but the applicability of the mechanism to the observed sample. Section 5.1 explicitly states that the system must have assembled early enough to cross the J2 resonance, yet the paper never estimates stellar ages or J2(t) for the six hosts, and the integrations all start at J2 = 10^-3, which corresponds to a very young, rapidly rotating star. If hot Jupiters predominantly arrive late via high-eccentricity migration, as the paper's own introduction suggests (Zink & Howard 2023), the J2 resonance would already be in the past and the claimed transit-selection bias would not operate. The numerical test with a lower initial J2 would directly show whether the bias depends on early assembly; a complementary gyrochronology calculation would show whether real systems satisfy the condition. This does not invalidate the paper, but it means the headline observational conclusion is conditional on formation pathway, matching the reader's conditional verdict. The J2 ramp's lack of a physical stellar evolution model and the absence of uncertainty quantification are secondary concerns; they affect confidence but not the logical structure. No fatal internal inconsistency was found.","tokens_in":16859,"tokens_out":32131,"duration_ms":335972,"concrete_test":"Re-run the Section 4.1 case-study suites with all parameters unchanged except that the J2 ramp starts at 10^-4 instead of 10^-3 (i.e., the system assembles after the star has already spun down through the relevant resonance), keeping the same logarithmic range, duration, and obliquities; if the outer-planet transit-probability curves in Figure 5 no longer decline with increasing lambda, the decoupling bias is confined to systems assembled before J2 fell below about 10^-3, and the headline claim must be restricted to early-assembled (in situ/disk-migration) hot Jupiters.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 5.1 states the mechanism requires \"the system must have assembled early enough that the planets were in roughly their observed orbital locations when the star's J2 attained the value that would have caused the onset of a secular resonance.\" The headline claim—that missing outer companions to hot Jupiters are a transit-selection bias driven by stellar spin-down and obliquity—therefore holds only if the current architecture was in place before J2 dropped below the resonance value. The paper does not quantify this for the six observed systems, nor does it connect the J2 = 10^-3 initial condition used in the integrations to a physical stellar age. Because the paper itself cites Zink & Howard (2023) arguing that most hot Jupiters arrive late via high-eccentricity tidal migration, a substantial fraction of systems may have assembled after J2 was already below resonance; for those systems the modeled decoupling would not occur and the observed dearth of outer companions would not be explained by this mechanism. This is an explicit scope limitation in §5.1, but it is not tested, making it the weakest load-bearing premise of the observational conclusion.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies whether stellar spin-down, through a decaying stellar quadrupole moment J2, can remove adjacent planetary companions of hot Jupiters from the transiting plane. Using second-order secular theory and N-body integrations with REBOUND/REBOUNDx, the authors simulate three-planet configurations based on WASP-47 and then inject a test outer companion into the six known hot-Jupiter systems with inner companions. They report that inner companions almost always remain transiting with the hot Jupiter, while outer companions can be driven to mutually non-transiting configurations, with higher stellar obliquity increasing this effect. The central observational claim is that the observed dearth of outer companions to hot Jupiters is partly a transit-selection bias produced by stellar evolution.","tokens_in":17057,"tokens_out":3709,"duration_ms":40217,"significance":"If the conclusions are robust, the paper provides a concrete dynamical explanation for the striking architectural asymmetry in the small observed sample: all known hot-Jupiter systems with adjacent companions except WASP-47 show inner companions only. The modeling is a genuine forward calculation: observed masses, radii, and periods are inputs, and no parameter is adjusted to reproduce the observed absence of outer companions. The authors appropriately extend earlier work by Spalding & Batygin (2016) and Brefka & Becker (2021) to realistic three-planet hot-Jupiter systems and to the current observed census, and the use of standard, publicly available integration tools is a strength. The result, if it survives the robustness tests requested below, would be of interest to exoplanet demographics and to planning searches for non-transiting companions. However, the quantitative predictive content is currently limited by the ad hoc J2 decay schedule and by an explicitly stated but untested early-assembly condition.","major_comments":[{"comment":"The stellar quadrupole is implemented as a fixed logarithmic ramp from J2 = 10^-3 to 10^-8 over 10^5 years, with no tie to a physical stellar spin-down model, to stellar age, or to the specific stellar parameters in Table 1. Because the secular resonance is encountered during this ramp, the final inclinations depend on the crossing rate and on the initial J2 value; a different decay law, or an initial J2 set by Eq. (4) with realistic k2 and rotation evolution, could shift the boundaries between the decoupling, stable-transiting, and oscillating regimes in Figure 3. The authors should add sensitivity tests over initial J2 and decay rate, or map the adopted schedule onto stellar ages and spin-down tracks, and show that the headline conclusion for outer companions is not an artifact of this choice.","section":"§3.1, Figure 3"},{"comment":"The paper correctly states that the decoupling mechanism requires the planets to be in roughly their observed orbital locations while the stellar J2 is still large enough for the system to cross the secular resonance. This is a load-bearing condition for the observational interpretation, yet the manuscript never quantifies whether the six observed systems, or the simulated WASP-47 analogues, plausibly satisfy it. Since the authors themselves cite Zink & Howard (2023) for the view that most hot Jupiters arrive late via high-eccentricity tidal migration, a substantial fraction of hot-Jupiter systems may have assembled after J2 had already decayed below the resonance value; for those systems the modeled decoupling would not operate. Please either add a quantitative assessment of the assembly-time requirement (for instance, using stellar spin-down ages and migration timescales) or explicitly restrict the observational claim to systems that assembled early.","section":"§5.1"},{"comment":"The numerical validation currently consists only of varying the integration length. The results also depend on choices that are not varied: zero initial eccentricities, a test outer companion mass fixed at 10% of the hot Jupiter mass with zero initial inclination, and, in the WASP-47 suite, a fixed stellar obliquity of 10 degrees. Given that Eq. (6) depends sensitively on the companion mass and semimajor-axis ratio, the authors should demonstrate robustness of the transit-probability statistics to these choices, or at minimum state which conclusions are qualitative and which are quantitative.","section":"§3.1, §4.1"}],"minor_comments":[{"comment":"The definition of the sampling distribution is unclear and appears to have a units error: after defining ζHJ as a dimensionless scaling factor for the hot Jupiter semimajor axis, the text states ζHJ ∼ Uniform(aUSP,0,1.4), which mixes a semimajor-axis value with a dimensionless bound. Please restate the sampling procedure unambiguously.","section":"§3.1"},{"comment":"The sentence 'The term b(1)3/2 is a convergent Laplace coefficient. The m represents the mass of each planet...' is duplicated verbatim; one copy should be removed.","section":"§2, after Eq. (3)"},{"comment":"Equation (6) has typesetting issues (the resonance condition appears as 'J2 res≈' and the brackets are unbalanced), and the notation mixes m1, m2, and mj in the numerator. Please typeset the equation cleanly and verify that the masses are defined consistently with the two-planet criterion of Spalding & Batygin (2016).","section":"Eq. (6)"},{"comment":"The phrase 'the mean planet of the planets' should read 'the mean plane of the planets.'","section":"§4.1"},{"comment":"The caption and axis labels are confusing: the diagonal line labeled aHJ = aUSP is not explained as a forbidden boundary, and the color scale for the transit probability is not defined in the caption. Please clarify.","section":"Figure 3"},{"comment":"The table header contains the typo 'T ransit' for 'Transit', and the footnote for WASP-47 c should make explicit that m sin i is a lower limit on the mass everywhere it is used.","section":"Table 1"}],"recommendation":"major_revision","confidential_remarks":"The core mechanism is not new, and the paper's main contribution is the three-planet generalization and application to the observed sample. That contribution is worth publishing, but only after the J2 schedule and the early-assembly condition are addressed, since both directly control whether the observational conclusion holds."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know about 2505.11637. First, the central dynamical asymmetry — inner companions stay in the transiting plane, outer companions can be decoupled by the J2 secular resonance as the star spins down — is not new. Spalding & Batygin (2016) derived the resonance criterion, and Brefka & Becker (2021) and Chen et al. (2022) showed the decoupling behavior. The paper cites these correctly and doesn't overclaim the mechanism itself. Second, the actual contribution is a clean numerical parameter sweep: 498 WASP-47-like three-planet realizations, plus case studies injecting a test outer companion into the six known inner-companion systems under five stellar obliquities. That is a legitimate, useful extension of the existing program.\n\nWhat the paper does well: the numerical setup is standard Rebound with GR and J2, the transit probability maps in Figure 3 are informative, and the system-by-system application in Figure 5 gives concrete predictions — higher obliquity and larger a_outer/a_HJ lower the chance of seeing the outer companion in transit. The paper is also honest about where it doesn't apply: the warm Jupiter caveat in Section 5.2 and the explicit assembly-timing condition in Section 5.1. That honesty is appreciated.\n\nThe biggest soft spot is exactly that Section 5.1 condition. The mechanism only works if the planets were already in roughly their current orbits when J2 crossed the resonance value. The paper states this but never tests it against the observed systems. Given Zink & Howard (2023) argue most hot Jupiters arrive late via tidal migration, a substantial fraction of these systems may have assembled after J2 was already too low. The paper's headline observational conclusion — that the dearth of outer companions is partly a spin-down bias — is therefore conditional on early assembly, and the paper doesn't quantify what fraction of hot Jupiters could satisfy that. That's a real limitation, but the paper is explicit about it; it's a scope issue, not a hidden flaw.\n\nSecond-order issues: the logarithmic J2 decay from 1e-3 to 1e-8 over 1e5 years is imposed rather than tied to a stellar spin-down model. They validate against integration length only, not against the physical schedule. The transit probability maps have no error bars, but for a parameter study that's minor. No code or data is provided, which makes reproduction harder.\n\nWho this is for: anyone working on hot Jupiter companion demographics or transit detection biases. This paper doesn't overturn anything, but it sharpens the question and provides system-specific visibility maps that future RV or phase-curve searches can use. It deserves a serious referee. The ideas are not novel enough for a high-impact venue, but it's publishable in a specialized journal after addressing the assembly-timing caveat more quantitatively and ideally tying J2 to physical spin-down. I'd send it to review.","headline":"A solid, honest application of an already-known secular mechanism to hot Jupiter companion detectability; the central asymmetry is not new, but the system-by-system transit maps are useful, and the paper is upfront about its main limitation.","tokens_in":17623,"tokens_out":3179,"would_cite":true,"duration_ms":26762,"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":"Stellar spin-down hides outer planets around hot Jupiters","keywords":["hot Jupiters","planetary companions","secular resonance","stellar obliquity","stellar spin-down","transit probability","J2 quadrupole moment","exoplanet dynamics"],"falsifier":"Compare the transit-detection rate of outer companions in hot Jupiter systems split by measured stellar obliquity: if high-obliquity hosts do not show a deficit of transiting outer companions relative to low-obliquity hosts, the proposed bias is not operating at observable strength. A complementary check is that RV-detected outer companions should preferentially appear around high-obliquity stars if the bias is real.","tokens_in":16554,"feed_emoji":"🪐","tokens_out":5970,"duration_ms":56234,"temperature":0.7,"pith_summary":"The paper argues that the apparent dearth of outer planetary companions to hot Jupiters is in part an observational bias caused by stellar evolution. As a host star spins down, its quadrupole moment $J_2$ decreases, and for a hot Jupiter with an exterior companion the pair can pass through a secular resonance that tilts the outer planet out of the transiting plane while leaving the hot Jupiter and any inner companion transiting. Higher stellar obliquity amplifies this effect. If correct, transit surveys will systematically miss outer companions in high-obliquity systems, so the known sample of hot Jupiter multi-planet systems is not representative of what is actually there.","feed_headline":"Stellar spin-down hides outer planets around hot Jupiters","feed_subtitle":"As host stars age, their shrinking quadrupole moment can tilt exterior companions out of the transiting plane.","key_machinery":"The central mechanism is the time-varying stellar quadrupole moment $J_2 = (k_2/3)(\\Omega_{\\rm star}/\\Omega_{\\rm breakup})^2$, which decays as the star spins down. In the secular theory the $J_2$ term enters the diagonal of the planet-planet interaction matrix $B$, and as $J_2$ falls the eigenfrequencies of the hot Jupiter-outer companion pair can cross, exciting large mutual inclination oscillations. The crossing criterion, derived from Laplace-Lagrange theory, predicts that only the outer companion pair is vulnerable; the paper tests this prediction with N-body integrations that logarithmically decrease $J_2$ from $10^{-3}$ to $10^{-8}$ and compute final transit probabilities.","core_discovery":"On its own terms, the paper establishes that in a three-planet architecture with a hot Jupiter flanked by an inner and an outer nearby companion, the transit probability of the inner planet is nearly immune to the combined effect of stellar obliquity and the decaying stellar quadrupole moment, while the transit probability of the outer planet is not. The reason is that the hot Jupiter-outer companion pair encounters a secular resonance as $J_2$ decays, whereas the inner pair does not. Across 498 WASP-47-like N-body realizations and case studies of the six known hot Jupiter-plus-inner-companion systems, the outer companion frequently attains a non-mutually-transiting configuration, with the frequency increasing with stellar obliquity. The paper reads this as a dynamical explanation for why WASP-47 is the only known system with an outer adjacent companion, and as a caution that missing outer companions do not imply missing planets.","pith_inferences":["If the mechanism operates, radial-velocity or phase-curve searches for non-transiting outer companions should find more companions around high-obliquity hot Jupiter hosts than around low-obliquity hosts; this is a testable prediction the paper does not make.","The same decoupling physics may apply to other compact multi-planet systems with a high-mass inner planet and a low-mass outer planet, not only hot Jupiters, because the resonance criterion is mass-ratio dependent.","A completeness model folding in this obliquity-dependent bias could reconcile the apparent rarity of outer companions with formation models that predict them."],"forward_implications":["Inner companions to hot Jupiters in stable configurations are expected to remain mutually transiting; their absence in a system is not attributable to this spin-down mechanism.","Outer companions can be present yet non-transiting, so occurrence rates of outer companions derived from transit surveys are lower limits.","Higher stellar obliquity should correlate with lower detection rates of outer companions in transit.","When an outer companion is observed in transit, the system likely either assembled late or has low stellar obliquity, offering a timing constraint on hot Jupiter migration.","Warm Jupiter systems, with larger separations, are not expected to show this bias because $J_2$ precession falls steeply with distance."],"supporting_citations":[{"why":"Supplies the analytic criterion for when a decaying stellar $J_2$ excites a secular resonance between a hot Jupiter and an outer companion.","marker":"Spalding & Batygin 2016"},{"why":"Demonstrates the decoupling behavior of planets from the transiting plane under stellar spin-down, which the paper extends to hot Jupiter systems.","marker":"Brefka & Becker 2021"},{"why":"Provides the discovery and orbital parameters of WASP-47, the exemplar system with both inner and outer adjacent companions used throughout the simulations.","marker":"Becker et al. 2015"},{"why":"Shows similar inclination-decoupling dynamics, used as a comparison for the regimes identified in the numerical results.","marker":"Chen et al. 2022"},{"why":"Gives the physical expression for $J_2$ in terms of stellar rotation, which the paper adopts to model spin-down.","marker":"Ward et al. 1976"},{"why":"Provides the Laplace-Lagrange secular theory and the $B$-matrix formalism on which the analytic model is built.","marker":"Murray & Dermott 1999"},{"why":"Provides the numerical implementation of general relativity and $J_2$ effects used to evolve the simulated systems.","marker":"Tamayo et al. 2019"},{"why":"Provides the N-body integrator used for all numerical simulations in the paper.","marker":"Rein & Liu 2012"}],"fun_headline_variants":["Outer companions to hot Jupiters get tilted away by stellar aging","Stellar spin-down tilts outer companions out of transit","Aging stars hide outer planets around hot Jupiters","Why we see few outer companions to hot Jupiters","Hot Jupiter outer neighbors vanish from transit view"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the planets were already near their currently observed orbits when the host star's $J_2$ crossed the value that triggers the resonance; if hot Jupiters arrived late via migration, the decoupling would not have affected those systems.","fun_headline_variants_meta":{"raw":{"variants":["Outer companions to hot Jupiters get tilted away by stellar aging","Stellar spin-down tilts outer companions out of transit","Aging stars hide outer planets around hot Jupiters","Why we see few outer companions to hot Jupiters","Hot Jupiter outer neighbors vanish from transit view"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000827,"raw_usage":{"total_tokens":3616,"prompt_tokens":945,"completion_tokens":2671,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":561,"completion_tokens_details":{"reasoning_tokens":2590}},"tokens_in":561,"tokens_out":2671,"duration_ms":19466,"temperature":1.0,"reasoning_tokens":2590,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T20:51:24.986130+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compare the transit-detection rate of outer companions in hot Jupiter systems split by measured stellar obliquity: if high-obliquity hosts do not show a deficit of transiting outer companions relative to low-obliquity hosts, the proposed bias is not operating at observable strength. A complementary check is that RV-detected outer companions should preferentially appear around high-obliquity stars if the bias is real.","supporting_citations":[],"review_version":1}