{"id":"e3b4112d-7134-4816-a769-7b267dc1cc77","arxiv_id":"2412.07700","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Hilda asteroids in the planar CRTBP are shown to lie on 2D invariant tori around stable periodic orbits, with frequencies near 0.5, and the authors propose using these frequencies for classification.","lead":"The authors model Hilda asteroids as test particles in the Sun-Jupiter restricted three-body problem and find that they move on quasi-periodic tori around a stable family of periodic orbits. They argue that membership in the Hilda group is better diagnosed by these tori's frequencies than by traditional two-body orbital elements.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Classification claim lacks a control: frequency-based Hilda membership is tested only on asteroids preselected by the orbital elements it claims to outperform.","rationale":"The reader's weakest assumption is that the planar projection suffices despite inclinations up to 20 degrees; that is a genuine and acknowledged limitation. My stress-test identifies a more fundamental evidentiary gap: the paper's headline classification claim is a comparative statement, but the experiment contains no comparison population. The sample is selected by the very orbital-element criterion the paper proposes to replace, so the observed confinement on tori can only show that the frequency method is consistent with the existing classification on that preselected set—not that it is better. This is not a numerical-error or internal-consistency objection: the frequency-analysis accuracy checks in Table I are impressive, the ERTBP continuation is a thoughtful robustness test, and the T-PSP evidence for invariant islands is credible. However, the claimed superiority of frequency-based membership requires a false-positive/false-negative benchmark, and none is presented. The planar projection concern could also be settled by extending the benchmark to spatial models for high-inclination objects, but the absence of a control undermines the central claim even before that extension. Therefore I keep the reader's CONDITIONAL verdict: the core dynamical picture is likely correct, but the abstract's classification claim should be downgraded or explicitly labeled as a proposal until a control experiment and a concrete decision rule are provided. Releasing the code and data would make the proposed test straightforward.","tokens_in":20525,"tokens_out":4663,"duration_ms":49372,"concrete_test":"Build a benchmark set containing (i) the existing Hilda candidates and (ii) a control sample of ~100 asteroids from the JPL Small Body Database with semi-major axis, eccentricity, and inclination inside or near the Section I.C box but not in the 3:2 resonance with Jupiter. Apply the exact Section II change of coordinates to the planar CRTBP, integrate each object to 2^20 adimensional time units, and run the same frequency analysis used in Section III and Appendix A. Plot the resulting (ω1,ω2) pairs against the family curves in Figures 9–10 and compute a confusion matrix for membership in the Hilda island, also recording how many candidates are discarded by the collision criterion. If the control asteroids overlap the Hilda frequency island, the claim fails; if they scatter outside while the Hildas remain inside, the claim is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The dynamical results may be sound, but the abstract's central claim—that frequencies in the planar CRTBP are a better way to decide Hilda membership than two-body orbital elements—is not supported by the experiment presented. Section I.C selects the Hilda sample using exactly the two-body orbital-element box (3.7–4.2 AU, e<0.3, i<20°) that the paper argues is inferior. The six representative asteroids and the 40–50 asteroids per T-PSP are all drawn from that preselected set. Thus the T-PSP and frequency analysis establish that most orbital-element-defined Hildas lie on or near two-dimensional tori; they do not establish that the planar frequencies separate Hildas from non-Hildas. No control population is tested: no near-boundary non-Hildas, no asteroids with similar a,e,i but outside the 3:2 resonance, and no quantified false-positive or false-negative rates. No operational decision boundary in (ω1,ω2) is given. Appendix A also discards asteroids that collide during integration, saying they are considered not quasi-periodic, but the number or fraction of discarded candidates is never reported. The planar-projection limitation noted in the Conclusions ('Some asteroids exhibit significant inclinations, making it essential to incorporate a three-dimensional model') is real, but it is secondary: even in the planar model, the comparative classification claim is untested.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper analyzes the Hilda asteroid group in the planar Circular Restricted Three-Body Problem (CRTBP) and the planar Elliptic RTBP (ERTBP). It numerically computes a family of stable periodic orbits in the CRTBP and, using Poincaré sections and frequency analysis, shows that six representative Hilda asteroids are surrounded by islands of two-dimensional quasi-periodic invariant tori. In the ERTBP, the same asteroids appear to lie on three-dimensional tori, with the Jupiter eccentricity adding a frequency equal to 1. The authors propose that membership in the Hilda class should be decided by the two dominant frequencies in the planar CRTBP rather than by two-body orbital elements.","tokens_in":20759,"tokens_out":6757,"duration_ms":58192,"significance":"The paper's dynamical machinery is careful and the confinement result for the preselected Hilda sample is credible: frequency combinations are checked at the 10^-14 level (Table I), the Taylor integration is run with local threshold 10^-16, and the CRTBP/ERTBP comparison is internally consistent. If limited to the statement that the orbital-element-defined Hilda sample lies in islands of quasi-periodic motion, the contribution is solid and relevant to the dynamical classification of asteroids. However, the abstract's stronger claim that the planar CRTBP frequencies are a better membership criterion than two-body elements is not supported by the experiments: the sample is preselected by the two-body element box, there is no control population, and no decision boundary is given. The paper itself notes in the Conclusions that a three-dimensional model is essential for high-inclination members, which further limits the planar frequency criterion.","major_comments":[{"comment":"The comparative classification claim is not tested. The Hilda sample used throughout the paper is selected in Section I.C using the two-body orbital-element box (3.7–4.2 AU, eccentricity below 0.3, inclination below 20°) taken from Zellner et al., and the six representative asteroids and the T-PSP groups of 40–50 asteroids per Jacobi-constant range in Section III.B are drawn from this preselected set. The paper therefore shows that most orbital-element-defined Hildas lie on or near two-dimensional tori, but it does not show that planar CRTBP frequencies separate Hildas from non-Hildas. No control population is integrated (no near-boundary objects, no asteroids with similar a, e, i outside the 3:2 resonance), no false-positive/false-negative rates are reported, and no operational decision boundary in (ω1, ω2) is provided. The abstract's claim that frequencies are 'much better' than two-body elements requires a comparative test on a mixed sample.","section":"Section I.C, Section III.B, Abstract"},{"comment":"The quasi-periodicity acceptance criterion is under-specified. Appendix A states that when the frequency-analysis output 'match[es] the input data' the motion is accepted as quasi-periodic, and that asteroids whose trajectories collide with the Sun or Jupiter are discarded as not quasi-periodic, but it does not report the number or fraction of discarded candidates in the T-PSP groups or in the full sample. Without this count, the conclusion that the Hilda group is 'confined' in an island of two-dimensional quasi-periodic solutions cannot be quantitatively assessed. Please give the tolerance used for the acceptance test and the fraction of discarded asteroids per group.","section":"Appendix A, Section III.B"},{"comment":"The planar projection is applied to asteroids with inclinations up to 20°, and the Conclusions explicitly state that 'Some asteroids exhibit significant inclinations, making it essential to incorporate a three-dimensional model into our analysis.' Since the proposed membership criterion is based on planar CRTBP frequencies, the effect of out-of-plane motion on the two dominant frequencies and on the island boundaries is unquantified; high-inclination members could be misclassified. Provide a quantitative estimate, for example by comparing the planar frequencies with the projection of a full three-dimensional integration for several high-inclination Hildas.","section":"Section II, Conclusions"}],"minor_comments":[{"comment":"In the discussion of the invariant relation (7), the statement that 'the integral term can be taken as zero' is correct only at the initial section f=0, not at f=2π; please clarify that the evaluation is made at f=0.","section":"Section IV.B"},{"comment":"The text and caption use 'medium value' where 'median value' is meant (the figure uses Q1, Q2, Q3 and the interquartile range).","section":"Section III (Figure 3)"},{"comment":"The epoch description is confusing: the coordinates are downloaded at MJD 59800 and must be propagated backwards to the chosen date MJD 58914; the phrase 'up to the date of interest' suggests forward propagation.","section":"Section I.C"},{"comment":"There are typos in captions: 'Analougous' in Figure 10 and 'collumn' in Figures 15 and 16.","section":"Figures 10, 15, 16"}],"recommendation":"major_revision","confidential_remarks":"The dynamical core is sound, but the abstract overstates the classification result. The missing control population is the key gap; it can be addressed by testing the frequency criterion on a mixed sample and reporting a decision boundary. The paper's own Conclusions are more cautious than the abstract, which suggests the overclaim is fixable. The paper is within scope for a dynamical-systems/celestial-mechanics journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe real content here is the tori picture, not the classification claim. The paper does a solid job showing that a sample of Hilda asteroids—selected by the usual orbital element box—lies on two-dimensional invariant tori of the planar CRTBP, and that the same picture survives Jupiter's eccentricity in the ERTBP. The frequency analysis is careful: Table I's integer combination errors around 1e-14, the Poincare sections are clean, and the continuation from CRTBP to ERTBP is a nice piece of work. I buy the dynamical result: Hildas are confined in islands of quasi-periodic motion around a stable periodic orbit family.\n\nThe soft spot is the abstract's last sentence: 'to decide if a given asteroid belongs to the Hilda class, it is much better to look at its frequencies...' That claim is not supported by anything in the paper. The sample is preselected by the same two-body orbital element box (3.7-4.2 AU, e<0.3, i<20°) the paper argues is inferior. No control population is tested: no near-boundary non-Hildas, no asteroids with similar a,e,i outside the 3:2 resonance, no false-positive/false-negative rates, no operational frequency boundary. So the paper demonstrates necessary conditions for a frequency-based classification, not the classification itself. This is not a fatal flaw in the dynamics; it's an overreach in the framing.\n\nTwo smaller issues. First, Appendix A discards asteroids that collide with Sun or Jupiter, but never reports how many were discarded—for a claim about 'the Hilda group,' the fraction matters. Second, the planar projection is acknowledged in the Conclusions as needing a 3D model, but the classification claim rests on planar frequencies; for inclinations up to 20°, that's a real unquantified effect. The authors note this, but it remains untested. Also, no code or data is released, only 'available upon reasonable request,' which makes the classification claim harder to verify independently.\n\nWho is this for? Celestial mechanicians and asteroid dynamics people. The tori computations are a genuine step for the Hilda group, and the method is portable to other resonant populations. The classification claim needs a control experiment: run the same frequency analysis on non-Hilda near-boundary asteroids, report the frequency spread of both groups, and give a decision boundary. That's a feasible revision, not a redo. I'd send this to review; a good referee will ask for exactly that. I wouldn't cite the classification claim as it stands, but I'd cite the tori result.","headline":"Solid tori computations for Hilda asteroids, but the frequency-based classification claim is tested only on the orbital-element-selected sample it claims to beat.","tokens_in":21294,"tokens_out":2777,"would_cite":true,"duration_ms":25429,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["37N05","70F07","70F15"],"pacs":[],"model":"deepseek-v4-flash","headline":"Hilda asteroids are best defined by their two in-plane frequencies in the Sun-Jupiter circular restricted three-body problem, not by a box of two-body orbital elements.","keywords":["Hilda asteroids","restricted three-body problem","quasi-periodic invariant tori","Poincaré section","frequency analysis","3:2 mean-motion resonance","Jacobi constant","asteroid classification"],"falsifier":"Take the full set of current Hilda candidates, propagate each in a three-dimensional model that includes Jupiter's eccentricity, and compare each object's two in-plane frequencies with the island boundaries found here: if any object whose planar frequencies sit inside the island escapes over a few million years, or if any object outside the island with in-island planar frequencies persists indefinitely, the criterion fails. A cheaper targeted check is to take a high-inclination member, with inclination near 20 degrees, and see whether its vertical motion moves its fundamental frequencies by more than the island's observed frequency width.","tokens_in":20317,"feed_emoji":"🪐","tokens_out":10352,"duration_ms":87692,"temperature":0.7,"pith_summary":"This paper tries to show that the Hilda asteroids, traditionally defined by a box of two-body orbital elements (semi-major axis between 3.7 and 4.2 AU, eccentricity below 0.3, inclination below 20 degrees), are better understood as a single dynamical object in the planar Sun-Jupiter circular restricted three-body problem. It finds a family of stable periodic orbits surrounded by islands of two-dimensional quasi-periodic motion, and shows that observed Hilda asteroids sit inside these islands, each tracing a closed curve in a suitable Poincaré section. The two dominant frequencies of that quasi-periodic motion lie in a narrow band near 0.5, and the same picture survives in the planar elliptic problem where Jupiter's eccentricity is included. If this is right, Hilda membership is a robust dynamical property of the Sun-Jupiter system, measurable by two frequencies, and the standard orbital-element classification can misassign boundary asteroids.","feed_headline":"Two frequencies, not orbital elements, tell a Hilda asteroid","feed_subtitle":"In a Sun-Jupiter three-body model the group sits in a stable two-frequency island, and Jupiter's eccentricity barely moves it.","key_machinery":"The load-bearing object is the family of stable periodic orbits around the Sun in the planar CRTBP, parameterized by Jacobi constant $C\\in[2.98,3.06]$, which acts as the skeleton of the Hilda group. Around each such orbit lies a Cantor family of two-dimensional invariant tori, a nearly continuous stack of tori with exponentially small gaps at resonances; their intersection with the section $\\Sigma$ is an invariant curve $\\varphi(\\theta)$ satisfying $P_C(\\varphi(\\theta))=\\varphi(\\theta+\\rho)$, computed as a truncated Fourier series by a Newton method. In the ERTBP, the central tool is the stroboscopic map taken at true anomaly $f=2\\pi k$, whose invariant curves represent two-dimensional tori of the flow, with rotation number fixed by the CRTBP period through $\\rho=2\\pi\\omega/\\omega_e=4\\pi^2/T$, and whose stability is decided by the generalized eigenvalue problem for $D_xP(\\varphi(\\theta))\\psi(\\theta)=\\lambda\\Gamma_\\rho\\psi(\\theta)$. Frequency analysis, using a Hanning-windowed discrete Fourier transform followed by collocation, extracts the two dominant frequencies and verifies that all other frequencies are integer combinations of them; the Thick-Poincaré Section Plot then lets the whole group of asteroids be compared by narrow ranges of $C$ without computing one section plot per asteroid.","core_discovery":"The paper's central claim is that the Hilda group occupies one coherent dynamical region rather than a box of orbital elements: a family of stable periodic orbits around the Sun in the planar Sun-Jupiter CRTBP, each surrounded by a Cantor family of two-dimensional invariant tori. For a fixed Jacobi constant $C$, the Poincaré section $\\Sigma=\\{y=0,\\dot y<0\\}$ shows each asteroid's trajectory as a single closed invariant curve, with neighboring curves forming concentric islands separated by a chaotic sea. The two fundamental frequencies of these tori are close to $0.5$, consistent with the underlying $4\\pi$ periodic orbit, and they vary only in a narrow range along the family; resonances appear as chains of islands and as horizontal segments or jumps in the frequency plot. Repeating the analysis in the planar ERTBP, where Jupiter's eccentricity acts as a $2\\pi$-periodic perturbation, turns the periodic orbits into two-dimensional tori and the two-dimensional tori into three-dimensional tori, with essentially the same two main frequencies. The authors conclude that membership in the Hilda class is better decided by these frequencies in the planar CRTBP than by two-body orbital elements.","pith_inferences":["A practical classifier could be built from the island's frequency boundaries: integrate a candidate in the planar CRTBP, extract its two dominant frequencies, and accept it as Hilda only if the pair falls inside the measured island; this would assign boundary objects like (164903) or (210340) unambiguously.","The frequency table's integer-combination checks, with residuals near $10^{-15}$ in the CRTBP and $10^{-14}$ in the ERTBP, suggest the quasi-periodic description is numerically sharp; the same tool could look for slow diffusion by testing whether the two frequencies drift over integrations longer than the $2^{20}$ time units used here.","If the planar criterion is adopted, high-inclination Hildas are the stress test: the vertical degree of freedom adds a third frequency, and the paper's plan to extend to three dimensions implies the planar island boundaries may move once inclinations are included.","The same machinery should transfer to other resonant asteroid groups, such as the 2:1 resonance or the Trojans, since the method only requires a stable periodic family and the surrounding tori."],"forward_implications":["Hilda membership becomes a dynamical statement about the Sun-Jupiter system: an asteroid belongs if its planar CRTBP trajectory lies on a closed curve in the Poincaré section at its Jacobi constant, inside the stable island.","The orbital-element box (3.7–4.2 AU, eccentricity below 0.3, inclination below 20 degrees) can misclassify boundary objects, and the frequency criterion removes the epoch dependence of those cuts.","Jupiter's eccentricity does not destroy the island: the ERTBP retains the same families as higher-dimensional tori with the same two main frequencies plus the forcing frequency 1, so the CRTBP result is robust to the main perturbation.","The two dominant frequencies of a Hilda asteroid sit close to approximately (0.505, 0.456) in the planar CRTBP, near the $4\\pi$ periodic orbit, giving an objective numerical signature for membership.","Resonances inside the island show up as chains of islands and as flats or jumps in the frequency plot, so the frequency analysis can separate asteroids trapped on resonant chains from those on regular quasi-periodic motion."],"supporting_citations":[{"why":"Provides the earlier long-period dynamical study of Hilda-type motion that this paper builds on and extends.","marker":"[1]"},{"why":"Supplies the equations of motion and the invariant relation for the elliptic restricted three-body problem used in the comparison.","marker":"[6]"},{"why":"Establishes the existence of a Cantor family of invariant tori around the stable periodic orbits that form the Hilda skeleton.","marker":"[5]"},{"why":"Gives the orbital-element box (3.7–4.2 AU, eccentricity below 0.3, inclination below 20 degrees) that the paper argues should be replaced by a frequency criterion.","marker":"[7]"},{"why":"Shows that lower-dimensional invariant tori persist under quasi-periodic perturbations, justifying the ERTBP family as an effectively continuous family.","marker":"[12]"},{"why":"Introduces frequency analysis, the method used to extract the two dominant frequencies of the quasi-periodic motions.","marker":"[20]"},{"why":"Provides the collocation Fourier method used to compute frequencies and amplitudes accurately in the numerical frequency analysis.","marker":"[18]"},{"why":"Gives the analytical error estimates for that collocation Fourier method, supporting the high-accuracy frequency checks.","marker":"[19]"},{"why":"Supplies the high-order Taylor integration method used to propagate asteroid orbits for $2^{20}$ time units.","marker":"[24]"}],"fun_headline_variants":["Hilda asteroids identified by two frequencies, not orbital elements","Two frequencies reveal Hilda asteroids in Sun-Jupiter model","Hilda class membership: frequencies beat orbital elements","How to tell a Hilda asteroid: look at two frequencies","Sun-Jupiter model: Hilda asteroids share a frequency signature"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire classification hangs on using only the flat, in-plane motion: asteroids with orbits tilted up to 20 degrees are projected onto the plane, and the paper does not quantify how much that projection shifts the frequencies or the size of the island.","fun_headline_variants_meta":{"raw":{"variants":["Hilda asteroids identified by two frequencies, not orbital elements","Two frequencies reveal Hilda asteroids in Sun-Jupiter model","Hilda class membership: frequencies beat orbital elements","How to tell a Hilda asteroid: look at two frequencies","Sun-Jupiter model: Hilda asteroids share a frequency signature"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000581,"raw_usage":{"total_tokens":2732,"prompt_tokens":939,"completion_tokens":1793,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":555,"completion_tokens_details":{"reasoning_tokens":1712}},"tokens_in":555,"tokens_out":1793,"duration_ms":10442,"temperature":1.0,"reasoning_tokens":1712,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T18:35:36.220874+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the full set of current Hilda candidates, propagate each in a three-dimensional model that includes Jupiter's eccentricity, and compare each object's two in-plane frequencies with the island boundaries found here: if any object whose planar frequencies sit inside the island escapes over a few million years, or if any object outside the island with in-island planar frequencies persists indefinitely, the criterion fails. A cheaper targeted check is to take a high-inclination member, with inclination near 20 degrees, and see whether its vertical motion moves its fundamental frequencies by more than the island's observed frequency width.","supporting_citations":[{"cited_title":"Schubart ,\\ title title Long-period effects in the motion of Hilda-type planets , \\ 10.1086/110605 journal journal Astron","cited_arxiv_id":null,"evidence_quote":"Provides the earlier long-period dynamical study of Hilda-type motion that this paper builds on and extends."},{"cited_title":"Szebehely ,\\ @noop title Theory of Orbits \\ ( publisher Academic Press ,\\ year 1967 ) NoStop","cited_arxiv_id":null,"evidence_quote":"Supplies the equations of motion and the invariant relation for the elliptic restricted three-body problem used in the comparison."},{"cited_title":"Jorba \\ and\\ author J","cited_arxiv_id":null,"evidence_quote":"Establishes the existence of a Cantor family of invariant tori around the stable periodic orbits that form the Hilda skeleton."},{"cited_title":"Zellner , author A","cited_arxiv_id":null,"evidence_quote":"Gives the orbital-element box (3.7–4.2 AU, eccentricity below 0.3, inclination below 20 degrees) that the paper argues should be replaced by a frequency criterion."},{"cited_title":"Jorba \\ and\\ author J","cited_arxiv_id":null,"evidence_quote":"Shows that lower-dimensional invariant tori persist under quasi-periodic perturbations, justifying the ERTBP family as an effectively continuous family."},{"cited_title":"G \\'o mez , author J","cited_arxiv_id":null,"evidence_quote":"Provides the collocation Fourier method used to compute frequencies and amplitudes accurately in the numerical frequency analysis."},{"cited_title":"G \\'o mez , author J","cited_arxiv_id":null,"evidence_quote":"Gives the analytical error estimates for that collocation Fourier method, supporting the high-accuracy frequency checks."}],"review_version":1}