{"id":"2eb9f186-5c85-494d-ac55-dde2ebcedb0d","arxiv_id":"2608.05457","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"The debiased mean plane of the Plutinos sits at inclination 1.9 degrees and node 46.8 degrees in the J2000 frame, separate from the local Laplace plane and consistent with the invariable plane, while the g+ and g- Kozai subgroups differ from each other.","lead":"Astronomers measured the average orbital plane of five groups of Kuiper belt objects locked in resonances with Neptune, using a debiased statistical method. The Plutinos, the largest resonant group, have a mean plane that differs from the standard Laplace plane prediction but is consistent with the solar system's invariable plane.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The abstract's high-confidence g+/g- mean-plane difference is unsupported: the body calls it tentative and the debiasing assumption is explicitly questionable for vZLK librators.","rationale":"The reader's verdict is CONDITIONAL, and my read does not move it. The paper's central result that the all-Plutino mean plane is far from the local Laplace plane is supported by two independent debiasing methods (SCT25 and VM17) that agree closely, and by a large synthetic sample showing the same qualitative pattern; that part of the claim is robust. The load-bearing weakness is the g+/g- mean-plane difference, which the abstract presents as high-confidence even though the body explicitly labels it tentative and warns that the SCT25 separability assumption is questionable for vZLK librators. The non-overlap of 99.7% credible contours is not a formal statistical test, and the summary admits the result depends on an assumption that may fail. A concrete test using the existing synthetic sample with a realistic selection function would settle whether the SCT25 method recovers the true subgroup mean planes when the separability assumption is violated. Because the body already contains the necessary caveats, the correct action is to revise the abstract and summary to match those caveats, rather than to reject the work. The condition is that the g+/g- claim be downgraded to tentative and that the formal statistical comparison be supplied or the claim withdrawn.","tokens_in":19252,"tokens_out":4715,"duration_ms":44553,"concrete_test":"Use the synthetic Lawler et al. (2025) Plutino population described in Appendix A to test the SCT25 estimator under the specific conditions that the concern targets: apply a realistic survey selection function (for example, an ecliptic-latitude and magnitude bias model matching current Plutino catalogs) to create a biased subsample of roughly 446 objects, then run the SCT25 method separately on the g+ and g- subsets. Compare the recovered mean planes to the true vMF mean planes computed in Table 4. If either recovery error is a substantial fraction of the claimed 12.27 degree separation (say, greater than 2 degrees), the observed g+/g- difference cannot be taken at face value. Repeating over many bootstrap realizations would provide the distribution of the recovery bias.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central novelty includes the claim that the g+ and g- vZLK-librating Plutino subgroups have different mean planes 'at high statistical confidence.' The evidence in Section 3 is that each subgroup's 95% credible region lies outside the other's 99.7% contour, but the text immediately states that this is only a visual observation and not a statistical comparison of means, so no P-value is attached. The Summary repeats the 99.7% language but adds that the result is tentative because the SCT25 method's key separability assumption (inclination independent of energy and scalar angular momentum) may fail for vZLK librators, where eccentricity and inclination exchange. The abstract omits these caveats and elevates a tentative, assumption-dependent result to 'high statistical confidence.' This matters because the 12.27 degree separation and the 144 degree node offset are the only evidence that Kozai libration center correlates with mean plane orientation. If the separability assumption fails, the debiased g+ and g- mean planes are biased and the separation could be partly or wholly a selection artifact.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper measures the debiased mean orbital planes of five resonant Kuiper belt populations (the 3:2 Plutinos and the 5:3, 7:4, 2:1, and 5:2 resonances) using the SCT25 maximum-likelihood vMF method, with VM17 and vMF comparisons for the Plutinos. The headline results are that the Plutino mean plane (i = 1.88°, Omega = 46.83°, theta95% = 2.04°) is strongly separated from the local Laplace plane but is found consistent with the invariable plane under the SCT25 method, and that the g+ and g− vZLK-librating Plutino subgroups appear to have mutually inclined mean planes. The paper also reports vMF concentration parameters for the relative inclination distributions and compares the observed Plutino results with a synthetic Plutino population.","tokens_in":19463,"tokens_out":10779,"duration_ms":99669,"significance":"If the Laplace-plane separation is robust, this is an important empirical constraint: resonant populations need not align with the linear secular forced plane, and the Plutinos' proximity to the nu18 secular resonance is a plausible cause. The paper is strong in applying three independent estimators (vMF, SCT25, VM17) to the main sample, in providing code and data, and in explicitly discussing the assumptions and limitations of the SCT25 method. The robust core is the rejection of the Laplace plane, which is consistent across SCT25 and VM17. However, two of the abstract's most prominent statements are less secure: the invariable-plane consistency is method-dependent (VM17 gives P = 0.024), and the high-confidence g+/g− difference is explicitly described in the body as tentative and rests on an assumption the authors themselves flag as questionable for vZLK librators. These issues are fixable by rephrasing and additional analysis, but they are load-bearing for the paper's central claims.","major_comments":[{"comment":"The abstract's claim that the Plutino mean plane is “indistinguishable from the invariable plane” is not method-independent and is contradicted by one of the paper's own estimators. Table 1 gives P = 0.17 for the SCT25 estimate, which is a failure to reject at the 95% level rather than positive evidence for the null; VM17 gives P = 0.024, a 95% rejection; and the vMF estimate gives P = 0.0001. Since Section 5 states that SCT25 and VM17 yield “substantially similar results,” the disagreement on this headline point needs to be resolved or reported. The abstract and Summary item 1 should say “consistent with the invariable plane under the SCT25 method” and should quote the VM17 P-value, or otherwise qualify the word “indistinguishable.”","section":"§2, Table 1; §5, Summary item 1"},{"comment":"The abstract's “at high statistical confidence” for the g+ and g− mean-plane difference overstates the evidence in Section 3. The 12.27-degree separation and 144-degree node offset are supported only by the statement that each subgroup's mean plane lies outside the other's 99.7% credibility contour; the text immediately notes that “this visual observation is not a statistical comparison of means.” The same paragraph warns that the SCT25 separability assumption may fail for vZLK librators because of eccentricity-inclination exchange. Appendix A strengthens the caution: in the synthetic Plutino population the analogous g+/g− separation is significant at the 95% level but not at the 99.7% level. The Summary's “appear different at the 99.7% level” and the abstract's “high statistical confidence” should be replaced with “tentatively different” or “suggestive,” and the caveat about the SCT25 assumption should accompany the claim in every location, including the abstract.","section":"§3, Table 2; §5, Summary item 2"},{"comment":"The SCT25 P-values and credible regions in Tables 1–3 are conditional on the fitted concentration parameter gamma_est; Section 2 acknowledges this and states that marginalizing over gamma was computationally prohibitive. Because gamma_est is estimated from the same sample, the 95% ellipses in Figures 1, 2, and 4 and the P-values in Tables 1–3 understate parametric uncertainty. This is especially relevant for the g± subgroups with n = 71 and 66 and gamma_est = 22.0 and 17.1, where the width of the credible region is the primary basis for the g+/g− claim. The authors should at least provide a sensitivity check on the g+/g− separation as gamma is varied (for example, a profile or approximate marginalization) before using the 99.7% wording in the abstract.","section":"§2, Eqs. (1)–(4)"}],"minor_comments":[{"comment":"The abstract uses “95% confidence” for the two-degree uncertainty, while Section 2 explicitly adopts “credibility” for the SCT25 conditional posterior; the terminology should be aligned.","section":"Abstract; §2"},{"comment":"The text reports “ecliptic mean inclinations of 5.74 degrees for the g+ subset and 7.17 degrees for the g− subset,” but Table 2 lists the ecliptic inclinations as 6.93 degrees and 7.13 degrees; the quoted numbers are actually the separations from the invariable plane in the same table. Correct the wording or the table.","section":"§3, text near Table 2"},{"comment":"The three credible ellipses in Figure 2 are not individually labeled; the caption and legend should clearly identify which ellipse corresponds to the all-Plutinos, g+, and g− estimates.","section":"§3, Figure 2"},{"comment":"The sentence “It is also widely separated from the Laplace plane” has an unclear antecedent; rephrase to name the mean plane of all Plutinos explicitly.","section":"§3, paragraph after Figure 3"}],"recommendation":"major_revision","confidential_remarks":"I agree with the stress-test assessment: the g+/g− concern lands, and the VM17 discrepancy is real. The central Laplace-plane separation is a solid empirical result, and the paper is not circular in its method usage. The abstract currently overstates two of the paper's own results; if the authors restrict the abstract to the robust Laplace-plane separation and clearly label the g± result as tentative, I would view the paper as nearly ready."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The headline result is real: the Plutinos' mean plane sits about 4.8 degrees from the local Laplace plane, and both debiasing methods (SCT25 and VM17) agree on that separation. That is a solid measurement, and the first empirical look at mean planes for the 5:3, 7:4, 2:1, and 5:2 resonances, plus separate planes for the g+ and g- vZLK librators, is genuinely new. The paper ships code and data on GitHub, the methods are described carefully, and the authors are unusually upfront about the SCT25 method's limitations (conditional rather than marginal posterior, separability assumption). That is good practice.\n\nBut the abstract oversells two things. First, 'at high statistical confidence' for the g+/g- plane difference: in the body that is a visual 99.7% contour overlap, explicitly 'not a statistical comparison of means,' and the Summary itself calls it tentative because the SCT25 separability assumption is questionable for vZLK librators. There is no P-value attached. That is not 'high statistical confidence.' Second, 'indistinguishable from the invariable plane' for the Plutinos: that is only true for SCT25 (P=0.17). Their own VM17 estimate gives P=0.024, i.e., distinct at 95%. The body reports both, but the abstract resolves the tension in favor of SCT25 without saying so. A reader relying on the abstract gets a stronger conclusion than the paper's own statistics support.\n\nThe comparison with the Lawler et al. synthetic Plutino population is a nice check: the same ~142 degree node offset between g+ and g- appears there, though with much smaller mutual inclination. That suggests the pattern may be real rather than a pure selection artifact, but the observed separation of 12.27 degrees is much larger, and the synthetic sample also has a different inclination distribution, so I would not hang too much on it yet.\n\nOverall, this is a careful, reproducible measurement paper for Kuiper belt dynamics. It deserves a serious referee. The central Laplace-plane-vs-Plutino result will hold up. The abstract needs to be rewritten to match the hedging in the body, and the g+/g- claim should be labeled tentative, not high-confidence, until a proper statistical comparison of the two means (ideally with a marginal posterior over gamma) is done. I would engage with it.","headline":"Solid resonant-KBO mean-plane measurements, but the abstract overstates both the g+/g- difference and the invariable-plane consistency.","tokens_in":20023,"tokens_out":3378,"would_cite":true,"duration_ms":31594,"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":"The debiased mean plane of the observed Plutinos is indistinguishable from the Solar System's invariable plane and clearly separate from the local Laplace plane; the two Kozai-librating subgroups have distinct mean planes.","keywords":["Plutinos","mean orbital plane","Laplace plane","invariable plane","Kuiper belt","mean motion resonance","von Zeipel-Lidov-Kozai resonance","debiased statistics"],"falsifier":"Run a simulated Plutino population with a known forced plane through the real survey selection function and apply the SCT25 estimator: if the recovered mean plane does not fall within its own 95% credible region of the injected plane, the bias-cancellation argument fails. With LSST's expected 11.7-fold increase in Plutino detections the 95% region should shrink by roughly a factor of 3.4, so a future measurement that moves the mean plane away from the invariable plane would also contradict the current claim.","tokens_in":19021,"feed_emoji":"🪐","tokens_out":16520,"duration_ms":130310,"temperature":0.7,"pith_summary":"This paper tries to establish where the orbital planes of resonant Kuiper belt objects actually point, and whether they follow the predictions of linear secular theory. Applying a bias-insensitive maximum-likelihood method to the observed Plutinos, the population in Neptune's 3:2 mean motion resonance, it finds a debiased mean plane at $1.88^\\circ$ inclination with ascending node $46.83^\\circ$ in the J2000 frame, a 95% credible region about $2^\\circ$ wide, statistically indistinguishable from the Solar System's invariable plane and clearly separate from the local Laplace plane. It also reports that the two subgroups whose perihelion arguments librate around $+90^\\circ$ and $-90^\\circ$ have mean planes separated by $12.27^\\circ$, with the authors cautioning that the debiasing assumption is least safe for these Kozai-librating groups. For four other Neptunian resonances (5:3, 7:4, 2:1, and 5:2), the measured mean planes are consistent with the local Laplace plane, Neptune's plane, and the invariable plane. If these measurements hold, resonant populations become usable as independent probes of unmodeled gravitational perturbations, since no analytical theory currently predicts the forced plane inside a mean motion resonance.","feed_headline":"Plutinos' mean plane matches the invariable plane","feed_subtitle":"Bias-corrected Plutino plane aligns with the invariable plane, not the predicted Laplace plane.","key_machinery":"The load-bearing tool is the SCT25 maximum-likelihood estimator of a debiased mean plane. It assumes the source population's orbital poles follow a von Mises-Fisher distribution (a rotationally symmetric bell-shaped distribution on the unit sphere) about a single mean pole, and that inclination is separable from orbital energy and scalar angular momentum. With those assumptions, the log-likelihood of the observed orbit-normal unit vectors $\\hat{\\mathbf{J}}_i$ is $\\log L = \\gamma\\,\\hat{\\mathbf{m}}\\cdot \\sum_i \\hat{\\mathbf{J}}_i - \\sum_i \\log I_0(\\gamma \\sin\\theta_i)$, where $\\hat{\\mathbf{m}}$ is the mean-plane normal, $\\theta_i$ the angle between the object's position vector and $\\hat{\\mathbf{m}}$, and $I_0$ the modified Bessel function of order zero. The survey's pointing biases cancel between numerator and normalization, leaving only the relative-inclination term; maximizing the likelihood gives the mean plane, and a grid search at fixed $\\gamma$ gives the 95% credible region. The physical mechanism under test is the Laplace plane from linear Laplace-Lagrange secular theory, which works for non-resonant populations but has no established analog inside mean motion resonances.","core_discovery":"On its own terms, the paper's central discovery is that the mean orbital plane of the observed Plutinos, corrected for survey bias with the SCT25 method, lies at $i=1.88^\\circ$, $\\Omega=46.83^\\circ$ (J2000) with a 95% credible region of radius $2.04^\\circ$. That plane is statistically consistent with the Solar System invariable plane ($P=0.17$) but distinct from Neptune's orbital plane ($P=0.018$) and from the local Laplace plane given by linear Laplace-Lagrange theory ($P\\approx0$). The paper further claims that the two von Zeipel-Lidov-Kozai (vZLK)-librating subgroups, whose arguments of perihelion librate about $+90^\\circ$ and $-90^\\circ$, have mutually different mean planes, separated by $12.27^\\circ$ and lying outside each other's 99.7% credible contours, though this result is explicitly flagged as tentative because the method's separability assumption may fail for these subgroups. For the 5:3, 7:4, 2:1, and 5:2 resonances, the measured mean planes are not statistically distinguishable from the local Laplace plane, Neptune's plane, or the invariable plane.","pith_inferences":["Beyond the paper's stated results, the near-180-degree difference in the $g+$ and $g-$ subgroups' nodal longitudes suggests that the two Kozai libration centers couple to opposite orbital-plane orientations; that coupling is absent from existing circular, co-planar secular models and would be a new dynamical prediction to test.","The paper's simulated Plutino comparison shows a qualitatively similar but much smaller $g+/g-$ plane separation (about 1 degree), so the observed 12-degree separation, if it holds, would imply that real Kozai-librating Plutinos are dynamically special rather than representative of a uniformly filled resonance.","A testable consequence of the paper's reasoning is that secular-resonance proximity, not resonance occupancy, drives the deviation from the Laplace plane: binning the 5:2 resonant population by semimajor axis or measuring other resonances near $\\nu_{18}$ would check this.","The credible regions are conditional on a fixed concentration parameter, so marginalizing over the concentration parameter would likely widen them; the reported subgroup separation should be read as an optimistic bound until the method is extended."],"forward_implications":["The Plutino forced plane is not the Laplace plane: the debiased mean plane of the 3:2 resonance population is consistent with the invariable plane, consistent with the idea that the resonance and the nearby $\\nu_{18}$ secular resonance warp the local forced plane.","Kozai libration state is tied to plane orientation: the $g+$ and $g-$ subgroups have mean planes separated by 12.27 degrees, with their ascending nodes separated by 144 degrees measured on the whole population's mean plane.","The four other measured resonant groups (5:3, 7:4, 2:1, 5:2) show no detectable anomaly: their mean planes are consistent with the local Laplace plane, Neptune's plane, and the invariable plane at roughly the 95% level.","The intrinsic inclination distributions of resonant KBOs are wider than the apparent ones: the debiased von Mises fits imply Rayleigh widths of about 11 degrees for all Plutinos and 12 to 14 degrees for the two Kozai subgroups.","A larger sample will settle which reference plane is correct: LSST is expected to detect about 11.7 times as many Plutinos, shrinking the credible region by roughly a factor of 3.4 and making the subgroup separation testable."],"supporting_citations":[{"why":"Supplies the SCT25 maximum-likelihood method for bias-insensitive mean-plane estimation and the false-alarm probability used for credible regions.","marker":"Siraj et al. (2025)"},{"why":"Provides the comparison VM17 method and the previous measurement of the warped Kuiper belt mean plane that the new estimates are checked against.","marker":"Volk and Malhotra (2017)"},{"why":"Supplies the dynamical classification catalog of 453 Plutinos and other resonant TNOs that defines the observed sample.","marker":"Volk and Van Laerhoven (2024)"},{"why":"Identifies the doubly resonant Plutinos and their +90 and -90 degree argument-of-perihelion libration centers, defining the g+ and g- subgroups.","marker":"Malhotra and Ito (2025)"},{"why":"Provides the synthetic 4 Gyr Plutino population used as a model comparison for the observed mean planes.","marker":"Lawler et al. (2025)"},{"why":"Supplies the machine-learning-classified Kuiper belt objects and the non-resonant mean-plane measurements used for context.","marker":"Matheson and Malhotra (2023)"},{"why":"Establishes the warped Laplace plane of the Kuiper belt from linear secular theory, the main reference plane the Plutino result is contrasted with.","marker":"Chiang and Choi (2008)"},{"why":"Supplies the linear Laplace-Lagrange secular theory and the (q,p) notation used to define the Laplace plane and report mean planes.","marker":"Murray and Dermott (1999)"},{"why":"Provides the earlier von Mises-Fisher fit to Plutino orbital poles that the debiased method revisits.","marker":"Matheson et al. (2023)"}],"fun_headline_variants":["Plutinos tilt toward invariable plane, not Laplace","Kuiper belt's Plutinos follow invariable plane","Plutino plane matches invariable, not Laplace","Resonant KBOs' plane matches invariable, not predicted"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire debiased mean-plane measurement rests on the assumption that the source population's orbital orientations are symmetrically clustered around a single pole and are statistically independent of the orbits' sizes, shapes, and how the survey pointed; for the Kozai-librating subgroups, which trade eccentricity for inclination, the paper itself warns this independence is questionable.","fun_headline_variants_meta":{"raw":{"variants":["Plutinos tilt toward invariable plane, not Laplace","Kuiper belt's Plutinos follow invariable plane","Plutino plane matches invariable, not Laplace","Resonant KBOs' plane matches invariable, not predicted"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000536,"raw_usage":{"total_tokens":2669,"prompt_tokens":1136,"completion_tokens":1533,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":752,"completion_tokens_details":{"reasoning_tokens":1464}},"tokens_in":752,"tokens_out":1533,"duration_ms":9932,"temperature":1.0,"reasoning_tokens":1464,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T12:43:07.690433+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a simulated Plutino population with a known forced plane through the real survey selection function and apply the SCT25 estimator: if the recovered mean plane does not fall within its own 95% credible region of the injected plane, the bias-cancellation argument fails. With LSST's expected 11.7-fold increase in Plutino detections the 95% region should shrink by roughly a factor of 3.4, so a future measurement that moves the mean plane away from the invariable plane would also contradict the current claim.","supporting_citations":[{"cited_title":"Lido: Exploring the stable plutino parameter space","cited_arxiv_id":null,"evidence_quote":"Provides the synthetic 4 Gyr Plutino population used as a model comparison for the observed mean planes."}],"review_version":1}