{"id":"74902cc7-e864-45fe-b0e6-2de24e2afe6f","arxiv_id":"2412.18661","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Outer giant planets can tilt tightly packed inner planets out of the transit plane, raising the observed gap complexity and potentially explaining a known dichotomy in exoplanet systems.","lead":"This paper uses orbital dynamics simulations to show that a giant planet orbiting further out can tilt the orbits of inner planets, making some of them disappear from transit view. This could explain why compact multi-planet systems with outer giants look like they have more irregular spacings in Kepler data.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Line-of-sight selection is the load-bearing assumption: the reported OG-driven increase in gap complexity is computed along a specially favorable viewing geometry, not the observer-averaged geometry of the Kepler sample.","rationale":"The reader's weakest-assumption analysis correctly identifies the favorable line-of-sight choice as the most load-bearing issue. The paper's own text in Section 2.3 and footnote 4 states that the chosen geometry is not the gap complexity a single observer would see, but rather the geometry from which a system is most likely to be classified as a multi-planet system. Since gap complexity is only defined for systems with at least three transiting planets, the observational comparison is inherently conditional on detection. The simulations condition on a maximum-transit LOS, which is precisely the subset of observer directions where the OG-induced inclination growth has the largest chance of still leaving at least three planets transiting. A realistic observer direction would more often place the system below the multi-transit threshold, removing those snapshots from the observed gap-complexity distribution. Thus the reported average increases in time-averaged gap complexity may not be the same quantity that He & Weiss (2023) measure. This is not an internal inconsistency: the Laplace-Lagrange machinery is standard, the code is public, and the N-body spot checks provide useful support. The issue is external validity of the statistical comparison. The proposed test would settle it by averaging over observer directions and applying the same N >= 3 transit selection as the observations. Since the reader already assigned a conditional verdict and this concern supports that verdict, no change in outcome is recommended.","tokens_in":12214,"tokens_out":3553,"duration_ms":39287,"concrete_test":"For the same 12,500 initial conditions, replace the fixed favorable LOS with a Monte Carlo sampling of observer directions uniform on the sphere (or isotropic, matching Kepler's random target orientations). At each time, compute C only when at least three planets transit with inclination below arctan(R*/a) for a specified stellar radius; otherwise treat the snapshot as unobserved. Compare the observer-averaged ensemble mean difference between systems with and without an OG against Figure 5. If the positive difference is not preserved under this selection, the secular mechanism fails as a quantitative explanation of the observed dichotomy.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that secular forcing from OGs 'could account for' the He & Weiss (2023) dichotomy requires that the simulated ensemble-mean change in gap complexity survives the same observational selection as the Kepler/KGPS sample. Instead, Section 2.3 (footnote 4) fixes the line of sight to be orthogonal to the line of nodes and in the mean inclination plane of the STIP, the geometry from which the system is most likely to be seen as a multi-planet system. The paper acknowledges this is not what a single observer sees. This matters because the gap complexity is undefined unless at least three planets transit; OG-induced inclination growth both raises C when a system remains multi-transiting and frequently drops the system below the N >= 3 threshold (e.g., Figure 3, grey regions). Under a realistic distribution of observer directions, many high-amplification systems would be excluded from the observed sample, and the conditional mean difference could shrink or reverse. The current Figure 5 averages only over the favorable LOS, so the +0.036 to +0.062 ensemble-mean increases are not yet shown to be the quantity that enters the observed OG versus no-OG comparison.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a dynamical explanation for the observed trend that systems of tightly packed inner planets (STIPs) with exterior giant companions (OGs) have higher observed gap complexity than those without. The authors use second-order Laplace–Lagrange secular theory to evolve the inclinations of idealized STIPs with and without an OG, then compute the time-averaged gap complexity along a line of sight chosen to maximize the chance of seeing at least three transiting planets. Across 12,500 simulation pairs, they find ensemble-mean increases in gap complexity of +0.036 to +0.062 when an OG is present (Figure 5), and argue that this secular mechanism can account for the He & Weiss (2023) dichotomy between OG and stellar-companion (SC) systems.","tokens_in":12557,"tokens_out":3994,"duration_ms":42040,"significance":"If the central claim holds, the paper would provide a plausible post-formation dynamical mechanism for a statistically significant observational dichotomy, connecting secular inclination forcing to a population-level observable. The study is not circular: no parameters are fitted to the He & Weiss dichotomy; the model inputs (period-ratio spacing, Rayleigh inclination distribution, companion mass and semi-major axis ranges) come from independent literature, and the result emerges generically from forward secular evolution. The paper also ships reproducible code (GitHub and Zenodo) and includes limited N-body validation in Section 2.4. However, the quantitative conclusion rests on a specially chosen viewing geometry and on an unspecified stellar radius, which currently limits the strength of the comparison to the observed Kepler/KGPS sample.","major_comments":[{"comment":"This is load-bearing because the abstract and conclusions directly compare the simulated time-averaged gap complexity to the observed C distributions.","section":"Section 2.3, footnote 4 and Figure 5"},{"comment":"A referee cannot reproduce the results without this parameter.","section":"Equation (18) and Section 2.3"},{"comment":"This issue is load-bearing because the paper's stated goal is to explain a statistical dichotomy, not merely to show that an increase is possible in some geometries.","section":"Section 4.1"}],"minor_comments":[{"comment":"The GitHub repository is mentioned, but the paper should be self-contained for this key check.","section":"Section 2.4"},{"comment":"This is a clarity issue, not a correctness issue.","section":"Section 3, Figure 5 caption"},{"comment":"This is a minor caveat for the idealized setup.","section":"Section 2.3, Equation (13)"},{"comment":"The table caption notes the mass approximation, but the small number of systems is not discussed.","section":"Table 1"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a timely and interesting question, and the secular framework is well chosen. The main concern is that the headline quantitative result is computed under a special viewing geometry that is not the one under which the Kepler sample is observed, and the paper does not yet demonstrate that the result survives an observer-averaged or detection-weighted treatment. The unspecified stellar radius is a smaller but still load-bearing omission. If the authors can add an observer-averaged calculation (or a convincing argument that the favorable LOS is representative), and specify R* with sensitivity checks, the paper would be suitable for publication. The qualitative comparison to He & Weiss is acceptable for a 'plausibility' claim, but the language in the abstract and conclusions slightly overstates the strength of the comparison."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper makes a plausible case that secular inclination forcing from an exterior giant can raise the observed gap complexity of a tightly packed inner system. That is genuinely new: nobody has computed the time evolution of the Gilbert-Fabrycky gap complexity under this specific secular driver, and the He & Weiss dichotomy has been floating without a post-formation mechanism. The Laplace-Lagrange machinery is standard, the code is public on GitHub and Zenodo, and the authors test their secular model against N-body integrations near the Hill stability boundary. They also do the sensible control comparison: stellar companions have a much smaller secular effective mass, which lines up with the observed SC null result. The paper is honest about its own limitations, especially in Section 4.2. I would not want to referee this from a standing start; it deserves referees who know the literature.\n\nThe soft spots are real but addressable. The line of sight is always chosen orthogonal to the line of nodes and in the mean inclination plane of the STIP, which the authors admit is not what any single observer sees. The ensemble-mean increases of +0.036 to +0.062 in Figure 5 are averaged only over intervals where at least three planets transit. That is a conditional average, and the conditioning is favorable: the same mechanism that raises gap complexity also frequently drops the multiplicity below three, and the paper does not show that the observed Kepler/KGPS sample is selected along such a favorable direction. A proper forward model of detection geometry, or even a sensitivity analysis over random observer directions, is needed before the quantitative claim to explain the dichotomy is secure. The stellar radius R* is also never specified, and the comparison to He & Weiss is qualitative rather than a real mock-sample fit. These are not fatal flaws; they are the difference between demonstrating a mechanism and demonstrating the explanation.\n\nI think the central dynamical insight holds: secular inclination amplification from an outer giant is a plausible contributor to the observed dichotomy, and the SC comparison is a nice consistency check. The paper deserves a serious referee. I would send it out, with the expectation that the authors add a quantitative observer-geometry treatment or explicitly weaken the claim to a proof-of-principle.","headline":"A credible secular mechanism for the OG gap-complexity dichotomy, but the quantitative case depends on a favorable line-of-sight choice that needs a proper observational selection model before the claim fully lands.","tokens_in":12988,"tokens_out":1526,"would_cite":true,"duration_ms":17652,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Outer giant planets can create irregular spacings in inner exoplanet systems after formation.","keywords":["exoplanet dynamics","gap complexity","secular perturbations","Laplace-Lagrange theory","peas in a pod","outer giant planets","transit geometry","Kepler/KGPS sample"],"falsifier":"Take the same 2,500-condition ensembles and compute the gap-complexity difference using lines of sight drawn from an isotropic distribution of observer orientations instead of always the favorable perpendicular geometry; if the ensemble-average difference between STIPs with and without outer giants is no longer positive, the proposed explanation of the observed dichotomy fails. A second, observational check is to compare the predicted fraction of STIP+OG systems that would appear as one- or two-transiting-planet systems with the multiplicity rates in the Kepler/KGPS sample; a large mismatch would rule out secular forcing as the dominant cause.","tokens_in":12013,"feed_emoji":"🪐","tokens_out":14115,"duration_ms":113518,"temperature":0.7,"pith_summary":"Kepler surveys show that systems of tightly packed inner planets (STIPs), the 'peas in a pod' systems, tend to have less regular orbital spacings when a giant planet orbits farther out, and more regular spacings when the outer companion is a star. This paper tests whether that observed dichotomy can be produced after the planets form, purely by the long-term secular gravitational tug of the outer giant on the inner system. Using Laplace-Lagrange secular theory, the authors evolve 2,500 initial condition sets in each of five ensembles (12,500 simulations) both with and without an outer giant, and compute the time-averaged gap complexity $\\langle \\tilde{C}\\rangle$ along a favorable observing direction. They find that outer giants raise the average gap complexity by between +0.036 and +0.062 depending on STIP multiplicity, with the biggest increases coming from massive, close-in giants. If correct, secular dynamics alone can explain the statistical gap-complexity dichotomy, with no need to invoke formation-stage effects.","feed_headline":"Stirred by outer giants, inner planets show irregular spacings","feed_subtitle":"Secular tugs from an outer giant knock inner planets out of transit, inflating the observed gap complexity.","key_machinery":"The load-bearing machinery is second-order Laplace-Lagrange secular theory for the inclination degrees of freedom, in which an $N\\times N$ matrix $\\mathbf{B}$ built from Laplace coefficients determines the eigenfrequencies and eigenmodes of each planet's inclination vector. The strength of the outer companion's forcing enters through the secular effective mass, $m_k \\alpha_{jk}\\bar{\\alpha}_{jk} b^{(1)}_{3/2}(\\alpha_{jk})$, which is about two orders of magnitude larger for a typical giant planet than for a stellar companion; this ratio is what lets the model reproduce both halves of the observed dichotomy. The observable is the gap complexity $C$ of Gilbert & Fabrycky, a convex complexity that combines the Shannon entropy and disequilibrium of the normalized log-period spacings $p^\\star_i$. For each of 12,500 simulations, the paper compares the time-averaged value $\\langle \\tilde{C}\\rangle$ with and without the outer giant, using the line of sight perpendicular to the line of nodes and lying in the mean inclination plane of the inner system. The secular approximation is checked against N-body integrations with a Wisdom-Holman integrator and found to match up to the Hill stability boundary.","core_discovery":"The paper's central claim is that secular perturbations from an exterior giant companion can account for the gap-complexity dichotomy that He & Weiss (2023) measured in the Kepler/KGPS sample. The giant adds a new mode to the Laplace-Lagrange inclination solution of the inner system; that mode amplifies the planets' mutual inclinations, so each planet spends more time tilted out of the transiting plane. Since gap complexity is computed only from the planets currently seen in transit, a missing planet changes the normalized log-spacing weights $p^\\star_i$ in the definition of $C$ and produces artificial gaps. Averaged over long integrations and over each 2,500-realization ensemble, the presence of an outer giant increases the time-averaged gap complexity by +0.036 ($N=4$, $I_{\\rm OG}=10^\\circ$) to +0.062 ($N=6$, $I_{\\rm OG}=10^\\circ$), while in some individual parameter regions the sign is reversed. The same framework explains why stellar companions do not show the effect: their secular effective mass is roughly two orders of magnitude smaller, so the induced inclination forcing is too weak.","pith_inferences":["The favorable-line-of-sight assumption implies a testable census prediction: if secular forcing is the cause, STIP+OG populations should show an elevated fraction of systems observed with fewer transiting planets than their true multiplicity, compared with STIP-only populations.","Repeating the calculation with isotropically distributed observer lines of sight could shrink or reverse the average +0.036 to +0.062 shift, so the strength of the explanation depends on how strongly Kepler-style detection selects systems with many transiting planets.","Because only second-order secular theory is used, eccentricity-inclination coupling effects such as Lidov-Kozai oscillations are omitted; including them could alter the gap-complexity evolution on long timescales, particularly for highly inclined outer giants."],"forward_implications":["If correct, the observed STIP+OG gap-complexity dichotomy can be explained without invoking a formation-stage mechanism; post-formation secular forcing alone raises the time-averaged gap complexity by the amount seen in the Kepler/KGPS sample.","The model predicts that gap-complexity enhancement should be strongest for massive, close-in outer giants and absent or weak for distant stellar companions, matching the sample positions of OG and SC systems in the secular-effective-mass plane.","Because individual simulations show both positive and negative changes, the population-level trend is not a deterministic statement about any one system; individual STIP+OG systems can have lower gap complexity than their no-giant counterparts.","Systems whose inner planets are frequently knocked out of transit will sometimes be observed as one- or two-planet systems, so gap-complexity samples are biased toward the sub-population that remains multi-transiting, an effect the paper identifies as a limitation of the metric."],"supporting_citations":[{"why":"Defines the observational dichotomy--STIPs with outer giants have higher gap complexity and STIPs with stellar companions do not--that this paper seeks to explain.","marker":"He & Weiss (2023)"},{"why":"Introduces the gap complexity metric $C$ and its normalization, used as the paper's output statistic.","marker":"Gilbert & Fabrycky (2020)"},{"why":"Supplies the second-order Laplace-Lagrange secular theory and the form of the inclination evolution equations.","marker":"Murray & Dermott (1999)"},{"why":"Establishes the 'peas in a pod' period-spacing relation (spacing parameter $\\mathcal{P}=1.03\\pm0.27$) used to build the model STIPs.","marker":"Weiss et al. (2018)"},{"why":"Provides the empirically fit Rayleigh(2.5 degrees) initial mutual-inclination distribution from which STIP inclinations are drawn.","marker":"Fabrycky et al. (2014)"},{"why":"celmech computes the inclination eigenvalues and eigenvectors of the secular matrix in the Python implementation.","marker":"Hadden & Tamayo (2022)"},{"why":"Provides the rebound N-body integration suite used to validate the secular approximation.","marker":"Rein & Liu (2012)"},{"why":"Supplies the Hill stability criterion that bounds the region where the secular and N-body results agree.","marker":"Gladman (1993)"},{"why":"Presents N-body simulations of the formation stage in which outer giants reduce inner-system gap complexity, the alternative explanation this paper argues against.","marker":"Kong et al. (2024)"}],"fun_headline_variants":["Outer giants create fake gaps in inner planetary systems","Giant neighbors tilt inner planets, hiding them from view","How outer giants scramble the spacing of inner planets","Secular tugs from outer giants inflate observed exoplanet gaps","Outer giant's gravity masks inner planets, boosting gap complexity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quantitative results assume the observer's line of sight is always the most favorable one, perpendicular to the line of nodes and in the mean inclination plane of the inner system, rather than a random or typical observing geometry.","fun_headline_variants_meta":{"raw":{"variants":["Outer giants create fake gaps in inner planetary systems","Giant neighbors tilt inner planets, hiding them from view","How outer giants scramble the spacing of inner planets","Secular tugs from outer giants inflate observed exoplanet gaps","Outer giant's gravity masks inner planets, boosting gap complexity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000235,"raw_usage":{"total_tokens":1494,"prompt_tokens":932,"completion_tokens":562,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":548,"completion_tokens_details":{"reasoning_tokens":480}},"tokens_in":548,"tokens_out":562,"duration_ms":4790,"temperature":1.0,"reasoning_tokens":480,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T04:35:45.006824+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the same 2,500-condition ensembles and compute the gap-complexity difference using lines of sight drawn from an isotropic distribution of observer orientations instead of always the favorable perpendicular geometry; if the ensemble-average difference between STIPs with and without outer giants is no longer positive, the proposed explanation of the observed dichotomy fails. A second, observational check is to compare the predicted fraction of STIP+OG systems that would appear as one- or two-transiting-planet systems with the multiplicity rates in the Kepler/KGPS sample; a large mismatch would rule out secular forcing as the dominant cause.","supporting_citations":[],"review_version":1}