{"id":"670e7357-cfff-467a-bece-c21cc20e1bf0","arxiv_id":"2507.13072","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"New orbital and shape models of Emma and Pulcova reveal a dichotomy in large binary asteroid systems, with Emma showing signs of a low-density outer shell and the two systems belonging to different formation groups.","lead":"This paper studies two large binary asteroids, (283) Emma and (762) Pulcova, by fitting their moons' orbits and modeling their shapes. It finds that Emma's interior may be layered and that large binary asteroid systems split into two populations with different formation histories.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The dichotomy's -0.98 correlation (p=10^-5) is computed from only 7 systems after splitting on eccentricity; the p-value is below the exact permutation minimum, so its statistical support is overstated and needs robustness testing.","rationale":"After reading the paper in good faith, the strongest quantitative evidence for the central claim (a dynamical dichotomy in large binary asteroids) is the -0.98 Pearson correlation of Section 7.1. The Emma internal structure result is presented as conditional ('Assuming the shape model to be accurate') and the authors transparently discuss shape-model ambiguities in Appendix A; the dichotomy, by contrast, is asserted with a very precise p-value and no caveats. That correlation is the load-bearing element for the paper's main conclusion of two distinct formation pathways. My concern is not that the correlation is necessarily false, but that its reported significance is invalid at n=7: the p=10^-5 cannot be reproduced by an exact test, and no robustness checks are provided. The sample is also subject to post-hoc selection on the dependent variable (eccentricity is used both to define the group and to measure the correlation). These are concrete, checkable statistical issues rather than differences of physical interpretation. I give credit for the new orbital solutions and the careful compilation of astrometric data, which are independently useful regardless of the dichotomy claim. If the correlation survives Spearman and leave-one-out testing, the central claim is strengthened; if not, the dichotomy is not yet established. The reader's weakest_assumption focused on Emma's shape model; I partially agree with that concern but consider the dichotomy statistics to be the more load-bearing issue for the paper's title claim.","tokens_in":33282,"tokens_out":15794,"duration_ms":173368,"concrete_test":"Recompute the eccentricity–b/a correlation for the 7 'atypical' systems in Figure 10 using (i) Spearman's rank correlation, (ii) leave-one-out deletion of each system, and (iii) an exact permutation test (7! = 5040 permutations) for Pearson's r. Report the resulting rho, the range of PCC across leave-one-out trials, and the exact two-sided p-value. Also report the sensitivity to the eccentricity threshold used to define the atypical group (e.g., move the cut from 0.03 to 0.02 or 0.05). If the Spearman rho is above -0.7, or the leave-one-out PCC changes by more than 0.05, or the exact p-value exceeds 0.05, the claimed 'strong correlation' is not robust and the dichotomy's statistical foundation is weakened.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 7.1 claims a Pearson correlation coefficient of -0.98 (p = 10^-5) between satellite eccentricity and primary elongation b/a for the 'atypical' (eccentric) binaries. This is the primary quantitative support for the paper's central dichotomy claim. The correlation is computed over only 7 systems (Figure 10), and those systems were assigned to the 'atypical' group using the same eccentricity variable that appears on the y-axis of the correlation. With n = 7, the smallest possible exact permutation p-value for a one-sided test is 1/7! ≈ 1.4×10^-4; the reported p = 10^-5 is therefore unattainable under any exact test and must be an asymptotic t-approximation that is invalid at this sample size. The paper does not report Spearman's rho, leave-one-out stability, or a permutation-based significance, and the b/a and e values are drawn from heterogeneous sources (Vernazza et al. 2021, DAMIT, Thomas et al. 1996, Dobrovolskis 1996) with no uncertainties propagated into the test. The dichotomy—the central claim of the paper—thus rests on a small-sample correlation whose strength and significance have not been robustly established.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper combines re-reduced archival adaptive-optics astrometry, new lightcurve shape modeling, and genoid orbital fits to study the binary asteroid systems (762) Pulcova and (283) Emma, and then places them in a broader population of large binary asteroids. For Pulcova the satellite orbit is found to be nearly circular and co-planar, leaving the primary's J2 essentially unconstrained; for Emma the orbit is eccentric and yields a dynamical J2 = 0.11 +/- 0.01, which the authors compare with a homogeneous-shape J2 of about 0.14 to argue for a non-homogeneous internal structure with a low-density outer shell. The paper also claims an overall dichotomy in large binary asteroids, with one group showing a strong correlation between primary elongation and satellite eccentricity, and proposes two distinct formation pathways.","tokens_in":33581,"tokens_out":7453,"duration_ms":88472,"significance":"If the Emma inference is sustained, it would be an important new datum on internal structure of a carbonaceous main-belt asteroid, and the proposed population dichotomy would be a useful organizing hypothesis for future observations and simulations. The paper has concrete strengths: it makes the full astrometric and photometric measurements available in tables, shows residual plots for the orbital fits, documents the shape-model construction, provides occultation predictions, and explicitly discusses many limitations of the data. The Emma orbital fit has impressively low residuals, and the Pulcova analysis honestly lays out the degeneracy of the circular co-planar case. However, the two headline conclusions rest on statistical and shape-model supports that need strengthening before the claims can be accepted as stated.","major_comments":[{"comment":"The central dichotomy claim rests on the Pearson correlation of r = -0.98 (p = 10^-5) between satellite eccentricity and primary elongation b/a for the seven 'atypical' systems. With n = 7, no exact permutation test can produce a one-sided p-value smaller than 1/7! ≈ 1.98 × 10^-4, so the reported p = 10^-5 cannot be a valid significance statement; it must come from an asymptotic approximation that is inappropriate at this sample size. The test also uses the same eccentricity variable to select the group as appears on the y-axis, the b/a values are drawn from heterogeneous sources without propagated uncertainties, and no Spearman correlation, permutation p-value, or leave-one-out analysis is reported. I request an exact permutation test, a rank-based correlation, a sensitivity analysis to the group-selection criterion, and a discussion of how the selection affects the inferred relationship.","section":"§7.1 (Fig. 10)"},{"comment":"Emma's non-homogeneous internal structure is a load-bearing conclusion, but the comparison between the dynamical J2 = 0.11 ± 0.01 and the shape-model J2 ≈ 0.14 is not robust to the shape and radius uncertainties documented in the paper itself. The authors state that their own new shape model has an oblateness incompatible with the orbit, that the adopted Viikinkoski et al. (2017) model is a poor fit to the AO images, and that concavities or cratering could mimic a low-density crust; the alternative radius gives J2 = 0.13, and the core/crust densities in Figure 7 are chosen in an inversion that reproduces the adopted J2. I request a systematic propagation of shape-model and radius uncertainty into the homogeneous-body J2, for example by considering all admissible shape models and occultation scalings, and a correspondingly conditional statement of the internal-structure result. As written, the abstract's claim of a 'significantly non-homogeneous internal structure' overstates what the current data can support.","section":"§4.2, §5, and Appendix A.2.1"},{"comment":"The two-population dichotomy is presented as an established result, but the classification into 'typical' and 'atypical' systems is made using the same orbital eccentricity that is then correlated with shape, and the family-size distinction in §7.2 is not statistically quantified: the 'large family sequence' in Figure 11 is identified by eye with very few points after excluding Hektor. The two-pathway formation conclusion is a plausible hypothesis, but the paper should demonstrate that the dichotomy is not an artifact of the chosen eccentricity threshold and should quantify the separation of the family-size distributions. I encourage the authors to include the full sample in a robustness test or to explicitly describe the selection effects that determine which binaries have measured b/a and eccentricity values.","section":"§7.1–§7.5"}],"minor_comments":[{"comment":"The semi-major axis for Emma's satellite is quoted as 588.3 ± 0.0 km; please report a nonzero uncertainty or explain how this value and its rounding were obtained.","section":"Table 2"},{"comment":"The main text uses several different radii for Emma (67 km, 71 km, and 74.5 km) without a single table stating which reference radius corresponds to the headline J2 = 0.11; please harmonize these values and the corresponding shape-model reference.","section":"§4.2 and Conclusions"},{"comment":"The figures showing the shape-eccentricity and family-size relations would be much more informative with uncertainty bars on b/a and family-size estimates; currently the strength of the visual correlations appears larger than the data quality warrants.","section":"Fig. 10 and Fig. 11"},{"comment":"Some entries in the astrometry table, such as ΔM = 100.0 or negative ΔM values, appear unphysical; please explain how these epochs are treated in the outlier rejection and whether they enter the orbital fits.","section":"Table C.2"},{"comment":"The phrase 'Pearson Correlation Coefficient test of the linearity' is imprecise: the Pearson coefficient measures linear association but is not a test of linearity; please rephrase.","section":"§7.1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript contains valuable observational data and a clearly presented orbital analysis, but the two central claims—Emma's non-homogeneous interior and the global dichotomy—are stronger than the current statistical and shape-model support. Both issues are fixable within the scope of the paper, so I recommend major revision rather than rejection. The p-value problem in §7.1 is a clear numerical issue that should be corrected with exact or rank-based tests, and the Emma shape-model uncertainty needs to be propagated into the J2 comparison before the internal-structure claim can be accepted."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"New orbital solutions are the solid part of this paper. The authors re-reduced all archival AO data, assembled 68 and 56 satellite positions over 20- and 10-year baselines, and produced well-documented genoid fits. For Pulcova they get a significantly revised orbit (period 4.14 vs 4.44 d, mass 1.865 vs 1.4e18 kg) and a new shape model. For Emma they get the first clear J2 detection, 0.11 +/- 0.01, with very good residuals. This is careful, reproducible observational work and the astrometry tables alone are useful.\n\nThe dichotomy is a new synthesis, but its statistical support has a real problem. The paper claims a Pearson correlation of -0.98 with p = 10^-5 between satellite eccentricity and primary elongation for the 7 'atypical' systems. With n = 7, the smallest possible exact permutation p-value is about 2e-4, so the reported p-value must come from an asymptotic approximation that is invalid at this sample size. This needs to be fixed with a permutation test or Spearman plus leave-one-out, and the heterogeneous b/a sources need uncertainties propagated. As it stands, the central population claim is overstated. The dichotomy itself is plausible, but the evidence is not as strong as the abstract suggests.\n\nThe Emma interior claim is the second soft spot. The dynamical J2 is well measured, but converting it to a non-homogeneous interior with a low-density crust requires the shape model to be right. The authors themselves say their new shape model had oblateness incompatible with the orbit, the Viikinkoski model is a poor fit to some AO images, and concavities or cratering could mimic a low-density shell. They do acknowledge this in the body text, but the abstract states 'significantly non-homogeneous internal structure' as a result. That should be a hypothesis, not a headline.\n\nThe paper deserves a serious referee. The orbital solutions and shape models are valuable and should be published, and the dichotomy is worth discussing even if the statistics need revision. The authors are honest about limitations, and the citation pattern is fine. I would send it to review, with clear requests for corrected statistics and softened interpretive language.\n\nFor the reading group: yes, it is a good case study in small-sample correlation and shape-model uncertainty.","headline":"New orbital solutions are the real product; the dichotomy's headline p-value is statistically impossible at n=7, and Emma's low-density crust is a shape-model-dependent hypothesis.","tokens_in":650,"tokens_out":1168,"would_cite":true,"duration_ms":36653,"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":"Large binary asteroid systems split into two dynamical populations, and the orbit of (283) Emma's satellite reveals a low-density outer shell on its primary.","keywords":["binary asteroids","asteroid satellites","adaptive optics","orbital dynamics","internal structure","asteroid families","quadrupole J2","collisional formation"],"falsifier":"A high-cadence stellar occultation of (283) Emma that constrains its oblateness well enough to show the shape-only $J_2$ is actually about 0.11, matching the orbit, would remove the need for the low-density shell.","tokens_in":33130,"feed_emoji":"🛰️","tokens_out":12893,"duration_ms":130024,"temperature":0.7,"pith_summary":"This paper tries to show that the population of large binary asteroid systems is not a single class, and that two similar-looking systems, (762) Pulcova and (283) Emma, trace back to different impact histories. By re-reducing twenty years of adaptive-optics images and fitting the moons' orbits, the authors find that Emma's moon precesses measurably, giving a gravitational quadrupole $J_2 = 0.11 \\pm 0.01$ that is significantly lower than the $J_2 \\approx 0.14$ expected if the primary were uniform. They read this as evidence that Emma has a low-density, void-rich outer shell around a denser core, consistent with a catastrophic impact and re-accumulation, and they place its bulk density at $0.9 \\pm 0.3\\ \\mathrm{g\\,cm^{-3}}$. Pulcova's moon, by contrast, sits on a nearly circular, co-planar orbit that leaves the internal structure ambiguous, and the authors assign Pulcova to the 'typical' group of circular binaries. Across the wider sample of large binaries, satellite eccentricity correlates strongly with primary elongation in one group, and the authors argue that two distinct formation pathways, sub-catastrophic and catastrophic impacts, produced the observed dichotomy.","feed_headline":"Asteroid binaries split into two formation families","feed_subtitle":"Emma's moon exposes a low-density shell; Pulcova's circular orbit belongs to the other group.","key_machinery":"The central comparison is between two estimates of the primary's gravity field: one computed from the shape model under the assumption of a homogeneous interior, via spherical-harmonic decomposition, and one measured from the satellite's orbit through the non-Keplerian signal. The key coefficient is $J_2$ (defined by $C_{20} = -J_2$ in the expansion of the gravitational potential), which measures the equatorial bulge of gravity; for a satellite in a moderately eccentric, nearly equatorial orbit, $J_2$ drives a nodal precession rate $\\omega_P = -3 D_p^2 J_2 \\omega \\cos\\Lambda \\,/\\, [8 (a (1-e^2))^2]$, so a measured precession directly yields $J_2$ once the semi-major axis, eccentricity, and inclination are known. The orbital fits are made with the genoid genetic algorithm, and the shape-derived multipoles are obtained by spherical-harmonic decomposition of the topographic shape model. The dichotomy diagnosis uses the axial ratio $b/a$ of the primary and the satellite eccentricity $e$ as the two population-defining observables, with the correlation between them separating the two groups.","core_discovery":"The paper's central claim is that large binaries split into two dynamical populations, and that Emma and Pulcova sit in different ones. The load-bearing new result is Emma: with 56 satellite positions spanning ten years, a purely Keplerian orbit cannot fit the data, and the genoid orbital search detects nodal precession at the level $J_2 = 0.11 \\pm 0.01$ for the adopted radius; an alternative shape model with a 133 km diameter gives $J_2 = 0.13$. Because the shape model of Emma, assuming a homogeneous interior, predicts $J_2 \\approx 0.14$, the dynamical value is roughly one-fifth lower, and the authors interpret the shortfall as a non-homogeneous interior: a dense core surrounded by a low-density outer shell that is at least 30% void and may be up to 80% void. Pulcova's orbital solution, from 68 positions over twenty years, is nearly circular and co-planar, so $J_2$ is effectively unconstrained (values up to about 0.13 fit the data); the authors report a mass of $1.865 \\pm 0.019 \\times 10^{18}\\ \\mathrm{kg}$ and a density of $1.4 \\pm 0.2\\ \\mathrm{g\\,cm^{-3}}$. At the population level, the paper compiles about thirteen large binaries and finds a strong anticorrelation between primary elongation $b/a$ and satellite eccentricity ($r = -0.98$, $p = 10^{-5}$) among the eccentric group, while the circular group occupies a narrow elongation range; the eccentric systems tend to have large families and slower rotation, the circular ones small or absent families and fast rotation. These two clusters are interpreted as end states of catastrophic versus sub-catastrophic impacts.","pith_inferences":["Editorial inference: If the dichotomy is real, a cheap way to classify newly discovered large binaries is to measure satellite eccentricity alone, then predict both primary elongation and the presence of a large family from the two-sequence relation.","Editorial inference: The paper's proposed factor-of-five mass criterion for distinguishing 'smashed target' from 'escaping ejecta' families could be applied to the whole asteroid-family catalog; if the dichotomy holds, the two groups should separate cleanly under that classification.","Editorial inference: The tidal-equilibrium condition derived in the paper (roughly $\\rho_p/\\rho_s \\approx 0.83$ for equal tidal quality factors) implies that the physical densities of primary and satellite differ systematically in eccentric systems, a prediction that future component-density measurements could check.","Editorial inference: The sample of thirteen systems is small, so the natural stress test is to add newly characterized large binaries; the dichotomy predicts that new members will fall near the existing two sequences rather than in the gap between them."],"forward_implications":["If Emma's two-layer interior is real, eccentric-satellite binaries with large families are rubble-pile re-accumulations with porous outer shells, so their bulk densities substantially underestimate the density of the solid material.","The $b/a$–$e$ correlation offers a cheap diagnostic: for eccentric binaries, measuring either primary elongation or satellite eccentricity predicts the other and flags which formation pathway a newly found system follows.","Pulcova-type circular binaries will remain stubborn for interior studies: without a detectable precession signal, orbital fits cannot distinguish homogeneous from layered interiors, so shape models must be improved by occultation or disk-resolved data before $J_2$ can be trusted.","The dichotomy predicts that undiscovered large binaries should cluster into two sequences in shape–eccentricity space rather than filling the plane, and that eccentric systems should typically harbor large asteroid families.","Stellar occultations are the paper's designated next test: the authors publish predictions for Emma and Pulcova through 2029, arguing that these are essential to break the shape-model degeneracy."],"supporting_citations":[{"why":"Supplies the adopted shape model of Emma, whose homogeneous-interior $J_2 \\approx 0.14$ is the baseline that the dynamical $J_2 = 0.11$ is compared against.","marker":"Viikinkoski et al. (2017)"},{"why":"Provides the DAMIT database entry from which the adopted Emma shape model is taken.","marker":"Ďurech et al. (2010)"},{"why":"Gives the earlier Emma orbital solution that this work extends with a ten-year baseline and the first reported $J_2$ detection.","marker":"Marchis et al. (2008b)"},{"why":"Provides the earlier Pulcova orbital solution (4.44-day period, 703 km semi-major axis) that the new twenty-year solution revises.","marker":"Marchis et al. (2008a)"},{"why":"Describes the genoid genetic algorithm used for all orbital fits in this study.","marker":"Vachier et al. (2012)"},{"why":"Supplies the Kalliope two-layer core-crust internal structure model that the paper adopts as the template for Emma's interior.","marker":"Ferrais et al. (2022)"},{"why":"Presents impact simulations showing elongated, fast-rotating remnants with abundant satellite material, used to support the 'typical' population's formation pathway.","marker":"Walsh et al. (2025)"},{"why":"Provides impact and re-accumulation simulations that distinguish sub-catastrophic from catastrophic satellite-forming collisions, the paper's two formation-pathway templates.","marker":"Durda et al. (2004)"},{"why":"Identifies the Emma family and its 841 members, used as evidence that Emma formed through a catastrophic impact.","marker":"Milani et al. (2019)"},{"why":"Provides the tidal-circularization mechanism invoked to explain the circular, co-planar orbits of the typical binary systems.","marker":"Nesvorný et al. (2020)"}],"fun_headline_variants":["Two formation paths for large asteroid moons","Emma and Pulcova show two asteroid binary classes","Large binaries divide into two dynamical groups","Asteroid moons reveal dual origins for binaries","Asteroid binary dichotomy points to two impact histories"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The Emma interior conclusion collapses if the adopted shape model misrepresents the primary's true oblateness, since the authors' own alternative shape model was orbit-incompatible and concavities or cratering could mimic a low-density crust.","fun_headline_variants_meta":{"raw":{"variants":["Two formation paths for large asteroid moons","Emma and Pulcova show two asteroid binary classes","Large binaries divide into two dynamical groups","Asteroid moons reveal dual origins for binaries","Asteroid binary dichotomy points to two impact histories"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001099,"raw_usage":{"total_tokens":4713,"prompt_tokens":1201,"completion_tokens":3512,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":817,"completion_tokens_details":{"reasoning_tokens":3443}},"tokens_in":817,"tokens_out":3512,"duration_ms":29816,"temperature":1.0,"reasoning_tokens":3443,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T16:31:00.377151+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A high-cadence stellar occultation of (283) Emma that constrains its oblateness well enough to show the shape-only $J_2$ is actually about 0.11, matching the orbit, would remove the need for the low-density shell.","supporting_citations":[{"cited_title":"2017, , 607, A117","cited_arxiv_id":null,"evidence_quote":"Supplies the adopted shape model of Emma, whose homogeneous-interior $J_2 \\approx 0.14$ is the baseline that the dynamical $J_2 = 0.11$ is compared against."},{"cited_title":"2012, , 543, A68","cited_arxiv_id":null,"evidence_quote":"Describes the genoid genetic algorithm used for all orbital fits in this study."},{"cited_title":"2022, , 662, A71","cited_arxiv_id":null,"evidence_quote":"Supplies the Kalliope two-layer core-crust internal structure model that the paper adopts as the template for Emma's interior."},{"cited_title":"Satellite formation around the largest asteroids","cited_arxiv_id":"2505.03325","evidence_quote":"Presents impact simulations showing elongated, fast-rotating remnants with abundant satellite material, used to support the 'typical' population's formation pathway."},{"cited_title":"D., Bottke , W","cited_arxiv_id":null,"evidence_quote":"Provides impact and re-accumulation simulations that distinguish sub-catastrophic from catastrophic satellite-forming collisions, the paper's two formation-pathway templates."},{"cited_title":"2019, , 622, A47","cited_arxiv_id":null,"evidence_quote":"Identifies the Emma family and its 841 members, used as evidence that Emma formed through a catastrophic impact."}],"review_version":1}