{"id":"afa2313e-6a5f-45dd-8a58-4161afb56e9a","arxiv_id":"2506.09584","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Spatial extension of the Energy Transition Domain yields a database of millions of lunar ballistic captures, and a subset is transitioned into an ephemeris model as backup insertion trajectories for Lunar Trailblazer.","lead":"This paper builds a 3D database of ballistic capture trajectories around the Moon using the spatial Circular Restricted Three-Body Problem, then filters and transfers promising ones into a full Earth-Moon-Sun ephemeris model for NASA's Lunar Trailblazer mission backup insertion options. The value is a systematic, mission-tailored catalog of low-energy lunar capture paths that can be reached with modest trajectory corrections.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Printed Eq. (21) has a sign error before 2(1−μ)x; with that sign the Table 2 point cannot be on the ETD, so the derivation of all ETD initial conditions is suspect.","rationale":"The most load-bearing condition for the central claim is not the continuity of capture regions but the correctness of Eq. (21), because every database entry is generated from it. The printed equation appears to have a sign error that makes the stated example point not belong to the ETD, which would invalidate the entire IC generation and hence the database. I checked the derivation twice: expanding Eq. (19) and substituting into Eq. (4) gives a plus sign before 2(1−μ)x; the published minus sign makes |sinσ|>1 for the paper's own Table 2 point. This is an internal inconsistency, not a disagreement with prior consensus. The reader's continuity concern is genuine and worth testing, but it only threatens the comprehensiveness of the database; the sign error threatens every computed BC. I therefore recommend keeping the verdict CONDITIONAL, with the condition tightened to require correcting Eq. (21) and re-verifying that the database ICs satisfy ε2=0 and the stated C_J. The apparent rε=0.1 vs 0.3 discrepancy in Table 2 is a further reason to treat equation-level details as unverified.","tokens_in":28237,"tokens_out":23924,"duration_ms":225348,"concrete_test":"Independently re-derive Eq. (21) from Eqs. (19) and (4). Evaluate the Table 2 point with both the printed and the corrected sign; only the corrected sign satisfies |sinσ|≤1 for some cosζ>0 and reproduces the sphere-intersection geometry of Fig. 5. Then take one Appendix A initial condition, invert Eqs. (9), (12), (13), and (23), and check that the resulting (x,y,z,Γ,ζ) satisfies both ε2=0 and C_J. If the printed sign is confirmed, recompute a capture-set section such as C(Γ=0.52,z=0.2,ζ=0) with the corrected formula and compare with the published Fig. 9; a mismatch would show that the database is built from inconsistent ICs.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Section III.B, Eq. (21) is internally inconsistent with Eqs. (19) and (4). Substituting (19) into (4) and using x=x2+1−μ yields c=[2(1−μ)/r1 + 2(1−μ)x − (1−μ)^2 − C_J]/(2v2 cosζ), whereas the printed equation has −2(1−μ)x. The sign matters: for the ETD point in Table 2 (x2=−0.02, y2=−0.25, z=0.1, Γ=1), the printed formula gives |c/A|≈26/cosζ>1 for all admissible ζ, so Eq. (23) admits no η and the point is not on the ETD, contradicting Fig. 5 and Table 2. With the corrected sign, |c/A|≈0.72/cosζ≤1 within the allowed ζ range, so the point is on the ETD. Since Eq. (23) is the unique route from (x,y,z,Γ,ζ) to the initial conditions populating every capture set C(Γ,z,ζ), a wrong sign in Eq. (21) would miscompute every BC in the database. The continuity assumption in Algorithms 1–2 is a valid secondary completeness concern, but the sign error is more fundamental and should be resolved first.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper extends the Energy Transition Domain (ETD) method from the planar to the spatial Circular Restricted Three-Body Problem, defines spatial ballistic capture sets C(Γ,z,ζ), and constructs a database of roughly 200 million CR3BP capture initial conditions. It then introduces a mission-specific distance metric, filters candidates dynamically similar to the Lunar Trailblazer insertion phase, transitions them into an Earth–Moon–Sun ephemeris model (yielding about 20 million ephemeris BCs), and analyzes the resulting set, including sample trajectories with long capture durations, polar perilunes, and repeated close approaches. The central claims are that the ETD framework can be systematically extended to three dimensions to produce a comprehensive BC database, and that mission-relevant subsets can be transitioned into an ephemeris model to provide practical backup insertion options for Lunar Trailblazer.","tokens_in":28521,"tokens_out":10487,"duration_ms":113076,"significance":"If correct, the paper provides a large three-dimensional lunar ballistic-capture catalog and a repeatable pipeline from CR3BP initial conditions to ephemeris-based mission design. The two fmincon-validated transfers are a genuine strength because they give an independent check that the distance metric selects candidates reachable with tens of meters per second, and the Appendix A initial-condition tables support reproduction. The significance is currently tempered by a sign error in the printed core equations and by unverified algorithmic assumptions behind the completeness claim; once those are resolved, the framework would be a useful contribution to low-energy trajectory design.","major_comments":[{"comment":"The printed expression for c in Eq. (21) contains a sign error before the term 2(1−μ)x. Substituting Eq. (19) together with Eq. (9) into Eq. (4) gives c = [2(1−μ)/r1 + 2(1−μ)x − (1−μ)^2 − C_J]/(2 v2 cosζ), not the printed form with −2(1−μ)x. I verified the arithmetic on the Table 2 point (x2=−0.02, y2=−0.25, z=0.1, Γ=1): the printed formula yields |c/A| ≈ 26/cosζ > 1 for every admissible ζ, so Eq. (23) has no real solution and that point would not lie on the ETD, contradicting Fig. 5 and Table 2; with the corrected sign, |c/A| ≈ 0.72/cosζ, which falls below unity over part of the admissible ζ range. Because Eqs. (23), (9), (12), and (13) generate every initial condition that is propagated to build the capture sets C(Γ,z,ζ), this is a load-bearing error. Please correct the sign and provide a consistency check showing that the database entries were generated with the correct expression rather than with the printed one.","section":"Section III.B, Eqs. (21)-(24)"},{"comment":"The search-space initialization S ← C_{i,j−1,k} (Algorithm 1, line 9) and S ← C_{i,j,k−1} (Algorithm 2, line 13), followed by boundary expansion only while BCs are found on ∂S, implicitly assumes that the capture set varies continuously, and remains overlapping from one z- or ζ-slice to the next. The paper gives no test of this assumption. If a capture region becomes disconnected or jumps at a parameter step, the boundary expansion can stop before ever encountering it, and the claimed completeness of the database fails silently. Please add a validation, for example independent full-grid searches over selected (Γ,z,ζ) planes with a finer step, or a quantitative bound on how far capture regions can move between the adopted steps Δz=4×10^−3 and Δζ=1°.","section":"Section IV.B, Algorithms 1 and 2"},{"comment":"No resolution or convergence study is reported for the numerical parameters h=4×10^−4, d_O=2×10^−3, ΔΓ=0.02, and the doubled steps used in Algorithm 2. Since the expansion loops sample only boundary vertices (Algorithm 1, line 12), a capture feature smaller than the sampling spacing or located between vertices can be missed regardless of the offset. This is directly relevant to the paper's repeated claim of an 'extensive' and 'complete' database. Please report a convergence or sensitivity analysis, for example repeating a representative parameter slice at h/2 and h/4 and comparing the resulting capture sets and counts, or soften the completeness claim accordingly.","section":"Section IV.D, Table 3"}],"minor_comments":[{"comment":"In Eq. (15), the first two dotted components appear to carry the wrong subscript: the inverse transformation v2 = v + k × r2 should involve the synodic velocity components (ẋ, ẏ, ż), not the inertial components (ẋ2, ẏ2, ż2). Please clarify the notation.","section":"Section III.B, Eq. (15)"},{"comment":"The sentence stating that selecting all BCs below an inflated threshold 'can be guaranteed' to pick all candidates compatible with a required orbit is stronger than the approximations in Table 4 support; I recommend rephrasing this as a heuristic or providing a formal justification.","section":"Section V.C"},{"comment":"The roughly 2% of transitions that fail are attributed to ephemeris-model differences, but no spatial or temporal characterization of these failures is reported. Please clarify whether the failures are uniformly distributed or concentrated and whether they affect the final mission-specific subset.","section":"Section VI.C, Algorithm 3"},{"comment":"The fmincon validation of the two three-impulse transfers is a useful independent check, but the paper does not tabulate the maneuver epochs and magnitudes; adding them would allow independent verification of the reported Δv values.","section":"Section VII.B"},{"comment":"The percentages in Table 9 do not appear to sum to 100, and the categories in Table 8 mix 'total revolutions' with '2 revs ... 8 revs' columns; please clarify the categorization and rounding.","section":"Tables 8 and 9"},{"comment":"The initial-condition tables are a welcome reproducibility aid; the paper would benefit from also stating the integrator tolerances, force model settings, and Spice kernels used for the ephemeris propagations.","section":"Appendix A"}],"recommendation":"major_revision","confidential_remarks":"This is primarily an astrodynamics and mission-design paper rather than a numerical-analysis contribution; if the target journal is math.NA, the editor may wish to weigh scope suitability. I found no indication of duplicate publication, and the mission application appears novel relative to the authors' earlier planar ETD work."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the spatial extension of the ETD is a real step forward and the Lunar Trailblazer application is genuinely useful, but there is a sign error in Eq. (21) that needs to be sorted out before I would trust any number in the database.\n\nWhat is new: the planar ETD from [21] is extended to 3D via two-sphere intersection in velocity space, giving a clean one-parameter (ζ) family of initial conditions. The capture-set computation with polygonal growth in z and ζ is a reasonable way to handle a 5D parameter space, and the transition to ephemeris using rotopulsating frames is standard but well executed. The sample trajectories in Section VII are interesting, and the fmincon checks (66.1 and 65.8 m/s versus metric estimates 44.2 and 57.4 m/s) give some independent evidence that the distance metric is in the right ballpark. Credit where due: the geometric derivation in Section III.A is elegant, and the paper ships full initial conditions in Appendix A.\n\nSoft spots, in order. First and most important: Eq. (21) has a sign error. Substituting (19) into (4) gives c = [2(1−μ)/r1 + 2(1−μ)x − (1−μ)^2 − C_J]/(2 v2 cosζ), not the printed minus sign. With the printed sign, the Table 2 point (x2=−0.02, y2=−0.25, z=0.1, Γ=1) gives |c/A| ≈ 26/cosζ > 1, so Eq. (23) admits no η and that point is not on the ETD, contradicting Fig. 5 and Table 2. Since Eqs. (23)–(24) are the route from (x,y,z,Γ,ζ) to every initial condition in the capture sets, this has to be resolved: either the printed equation is a typo (and the code uses the correct sign), or the code has the bug and the database is built from initial conditions that do not satisfy the stated Jacobi constant. The paper as written cannot be correct as-is.\n\nSecond, the completeness claim in Algorithms 1–2 rests on an unstated continuity assumption: capture regions are assumed to vary continuously with z and ζ, so expanding the previous slice by a fixed offset captures everything. If a region appears or jumps at a parameter step, it is silently missed. That is a genuine but secondary concern; it is addressable with a finer grid or a conservative offset, and it does not invalidate the sample results.\n\nThird, the database itself is not released. The paper describes about 200 million BCs in CR3BP and ~20 million in ephemeris, but only summary statistics and a few sample ICs are given. For a 'database' paper, that is a significant omission; the authors should at least make the capture sets or a query tool available.\n\nIf the sign issue is a typo, the paper is a solid contribution for astrodynamics readers working on low-energy lunar missions. If it is not, the initial conditions are systematically wrong, and the conclusions about LTB backups are unsupported. My recommendation: send it to peer review, but require the authors to (a) fix or explain Eq. (21) with a clear derivation, (b) state and justify the continuity assumption, and (c) commit to releasing the database or a representative subset.","headline":"Good spatial extension and a practical mission study, but the sign error in Eq. (21) undermines the printed derivation of every ETD initial condition and must be resolved first.","tokens_in":29034,"tokens_out":5856,"would_cite":false,"duration_ms":51957,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["70F07","70F15"],"pacs":[],"model":"deepseek-v4-flash","headline":"Extending the Energy Transition Domain to three dimensions yields a database of millions of lunar ballistic captures and practical backup insertions for Lunar Trailblazer.","keywords":["ballistic capture","Energy Transition Domain","spatial circular restricted three-body problem","lunar trajectory design","ephemeris model transition","Lunar Trailblazer","low-energy insertion","capture database"],"falsifier":"On a slice where the three-body energy parameter is near the value where capture regions begin to split, run a dense uniform grid over the full Energy Transition Domain section at out-of-plane parameter values just past a step, and check whether every ballistic capture found lies inside the polygon produced by the offset-expansion algorithm; finding even one capture outside that boundary would refute the claimed comprehensiveness.","tokens_in":28027,"feed_emoji":"🌙","tokens_out":7650,"duration_ms":85774,"temperature":0.7,"pith_summary":"The paper claims that the Energy Transition Domain, a geometric filter for ballistic capture previously used in the planar circular restricted three-body problem, can be extended to full three-dimensional motion and turned into a large database of spatial ballistic captures. It applies this database to Lunar Trailblazer, selecting candidates whose Earth-escape orbits resemble the mission's nominal trajectory and converting them into a high-fidelity ephemeris model. If the method works, mission designers gain a ready catalog of lunar insertion backups reachable with tens of meters per second of correction, including multi-revolution, polar, and repeated-close-approach captures. The authors report roughly 20 million high-fidelity captures within their 60 m/s selection threshold, and they show sample three-impulse transfers from the nominal trajectory costing about 66 m/s.","feed_headline":"Millions of 3D lunar captures mapped for backup insertion","feed_subtitle":"A spatial energy-surface method finds low-cost, multi-revolution capture paths reachable from Lunar Trailblazer's nominal trajectory.","key_machinery":"The load-bearing object is the Energy Transition Domain: the set of positions in configuration space, at fixed Jacobi constant, where the two-body Kepler energy of the spacecraft with respect to the Moon is exactly zero, a necessary condition for temporary capture. In the spatial case, each such position gives a circle of admissible velocity vectors, obtained as the intersection of two velocity-space spheres, with the circle parameterized by the declination of the Moon-relative velocity. The database is built by tracking capture sets across the three-body energy parameter, out-of-plane position, and velocity declination using a polygonal-boundary expansion algorithm, with the planar symmetry of the circular restricted three-body problem used to halve the search. A mission-specific distance metric estimates the mono-impulsive velocity change needed to reshape the Earth-escape orbit of a candidate into the nominal Lunar Trailblazer orbit, and a rotopulsating-frame transformation converts selected initial conditions into the high-fidelity ephemeris model.","core_discovery":"The central discovery is that the spatial Energy Transition Domain has a simple geometric description in velocity space: for a fixed Jacobi constant, the condition of zero two-body energy with respect to the Moon makes the allowed velocity vectors the intersection of two spheres, a one-parameter family parameterized by the declination of the Moon-relative velocity. Sampling this family over grids in the three-body energy parameter, the out-of-plane position, and the velocity declination, then propagating under the circular restricted three-body problem, produces capture sets that can be expanded with polygonal offsets to form a spatial database of ballistic captures. A subset selected by a distance metric based on Earth-centered osculating orbital elements is then transformed into a Sun-Earth-Moon ephemeris model through a rotopulsating-frame transformation. The resulting mission-tailored set contains millions of captures for Lunar Trailblazer, including long-lived corridors, stable polar captures, and trajectories with two nearly identical polar perilune opportunities that could serve as backup insertion points.","pith_inferences":["Editorial extension: the same Energy Transition Domain to ephemeris pipeline could be used to generate contingency capture catalogs for future lunar orbiters, landers, or crewed vehicles, since the database query reduces to matching desired Earth-escape orbital elements.","Editorial extension: the claimed comprehensiveness of the database depends on capture regions varying continuously across parameter steps, so independent dense-grid spot checks at a few parameter slices would be needed to confirm that no disconnected capture regions were missed.","Editorial extension: the longest reported captures rely on the Moon's orbital eccentricity and solar perturbations in the ephemeris model, so their lifetimes may change in even higher-fidelity models that include lunar gravity harmonics or solar radiation pressure; testing that sensitivity is a natural next step.","Editorial extension: clustering or machine-learning analysis of the roughly 20 million entry database, which the paper suggests as future work, could reveal dynamical families such as distant-retrograde-orbit-like or butterfly-like captures and make the catalog easier to navigate for mission planners."],"forward_implications":["The reported database gives low-energy lunar missions a precomputed catalog of millions of spatial ballistic captures, so backup insertion options can be looked up without repeating the full search.","For Lunar Trailblazer specifically, the identified captures include corridors lasting 45 or more revolutions and trajectories with two closely matched polar perilunes, meaning a missed first insertion burn could be followed by another opportunity.","Because the selection is based on Earth-escape orbital elements, the same pipeline can be rerun for other lunar or planetary missions whose nominal low-energy trajectory is known.","The sample three-impulse transfers from the nominal Lunar Trailblazer trajectory to backup captures cost about 66 m/s, close to the distance metric's estimate of 44 to 57 m/s, indicating the metric can pre-filter candidates before full optimization.","Successive small braking maneuvers of 5 m/s at perilunes can stabilize a chaotic capture into a longer-lived lunar orbit, suggesting that insertion maneuvers can be distributed across multiple revolutions."],"supporting_citations":[{"why":"Supplies the planar Energy Transition Domain method, the ballistic capture definition, and the polygonal capture-set algorithm that this paper extends to spatial motion.","marker":"[21]"},{"why":"Provides the circular restricted three-body problem equations of motion and Jacobi constant formulation used throughout the analysis.","marker":"[23]"},{"why":"Establishes that ballistic capture requires energy above the L1 threshold, which defines the range of the three-body energy parameter in the search.","marker":"[24]"},{"why":"Provides the direction-cosine formulation for transforming circular restricted three-body problem states into an ephemeris model.","marker":"[25]"},{"why":"Supplies the rotopulsating-frame transformation and non-uniform time correspondence used to map Energy Transition Domain initial conditions into full-ephemeris states.","marker":"[26]"},{"why":"Provides the Spice toolkit and DE430 ephemerides used to obtain Earth, Moon, and Sun positions in the high-fidelity propagation model.","marker":"[22]"},{"why":"Gives the mono-impulsive velocity-change sub-cost formulas that define the distance metric for pre-filtering ballistic capture candidates.","marker":"[28]"}],"fun_headline_variants":["Spatial capture map yields millions of lunar backups","Millions of 3D captures chart backup paths to Moon","3D capture database expands lunar insertion backup options","Lunar Trailblazer gains millions of spatial capture backups"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The search assumes that every capture region at the next out-of-plane parameter step overlaps the region found at the previous step, because the search starts from the previous capture set and only expands by a fixed offset; a capture region that suddenly appears or jumps would be missed.","fun_headline_variants_meta":{"raw":{"variants":["Spatial capture map yields millions of lunar backups","Millions of 3D captures chart backup paths to Moon","3D capture database expands lunar insertion backup options","Lunar Trailblazer gains millions of spatial capture backups"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000852,"raw_usage":{"total_tokens":3707,"prompt_tokens":953,"completion_tokens":2754,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":569,"completion_tokens_details":{"reasoning_tokens":2689}},"tokens_in":569,"tokens_out":2754,"duration_ms":21825,"temperature":1.0,"reasoning_tokens":2689,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T04:45:05.128134+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"On a slice where the three-body energy parameter is near the value where capture regions begin to split, run a dense uniform grid over the full Energy Transition Domain section at out-of-plane parameter values just past a step, and check whether every ballistic capture found lies inside the polygon produced by the offset-expansion algorithm; finding even one capture outside that boundary would refute the claimed comprehensiveness.","supporting_citations":[{"cited_title":"Ballistic Capture Analysis using the Energy Transition Domain,","cited_arxiv_id":null,"evidence_quote":"Supplies the planar Energy Transition Domain method, the ballistic capture definition, and the polygonal capture-set algorithm that this paper extends to spatial motion."},{"cited_title":"H.,Anintroductionto themathematics andmethodsof astrodynamics, Aiaa, 1999, pp","cited_arxiv_id":null,"evidence_quote":"Provides the circular restricted three-body problem equations of motion and Jacobi constant formulation used throughout the analysis."},{"cited_title":"Low energy transit orbits in the restricted three-body problems,","cited_arxiv_id":null,"evidence_quote":"Establishes that ballistic capture requires energy above the L1 threshold, which defines the range of the three-body energy parameter in the search."},{"cited_title":"Trajectory refinement of three-body orbits in the real solar system model,","cited_arxiv_id":null,"evidence_quote":"Provides the direction-cosine formulation for transforming circular restricted three-body problem states into an ephemeris model."},{"cited_title":"Assessment of dynamical models for transitioning from the Circular Restricted Three-Body Problem to an ephemeris model with applications,","cited_arxiv_id":null,"evidence_quote":"Supplies the rotopulsating-frame transformation and non-uniform time correspondence used to map Energy Transition Domain initial conditions into full-ephemeris states."},{"cited_title":"Fundamentals of astrodynamics,","cited_arxiv_id":null,"evidence_quote":"Gives the mono-impulsive velocity-change sub-cost formulas that define the distance metric for pre-filtering ballistic capture candidates."}],"review_version":1}