{"id":"3b2e6b49-a255-4dfb-9e7d-9f17ff184126","arxiv_id":"2411.11615","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Using linearized energy-optimal control, the paper bounds the set of forced periodic trajectories near an L2 halo orbit and ranks deviation directions by cost.","lead":"This paper maps the nearby looping paths a low-thrust spacecraft can fly around a known Earth-Moon halo orbit and still return to its starting point each orbit, along with each path's fuel cost. It finds a practical answer to a mission design question: dipping closer to the lunar surface at perilune is one of the more expensive deviations.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Validation does not cover the J* amplitude or any direction besides the 4th eigenvector, so Table 1's semiaxes and the perilune-cost ranking rest on unvalidated extrapolation.","rationale":"Reader's weakest_assumption correctly identified the one-direction validation. I am sharpening it: even in the tested direction, the validation stops well short of J*. The linear algebra is internally consistent; the issue is an empirical gap at the operating point of all quantitative claims. A secondary concern is that the perilune conclusion uses eigenvector cost as a proxy for the constrained optimization over initial states, but the validation test would also expose whether per-direction rankings are trustworthy. Since the central derivation is sound and the gap is remediable, the appropriate verdict remains CONDITIONAL.","tokens_in":7120,"tokens_out":15847,"duration_ms":166927,"concrete_test":"Run the same Newton-Raphson validation as in Figure 1 for each of the five nonzero eigenvectors of E* (eigenvectors 2-6 in Table 1), with initial deviation magnitude set to the corresponding semiaxis a_i = sqrt(2J*/γ_i) at J* = 3.51e-4 DU²/TU³. Record relative error |J_computed - J_linear|/J_computed. If any direction exceeds 1% error, or if the ordering of J_computed across eigenvectors differs from the γ_i ordering, then the hyperellipsoid extents and the perilune-cost conclusion are not validated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central object, the hyperellipsoid (1/2)δx0^T E*δx0 ≤ J* (Eqs 20-22), is exact for the linearized system, but its use as a bound for the nonlinear CR3BP requires the linear approximation to be accurate up to J*. The only validation (Figure 1) is along the 4th eigenvector and over |dx|≤1e-3 DU, for which the largest computed cost is ≈1e-4 DU²/TU³. Yet Table 1 reports semiaxes for J*=3.51e-4 DU²/TU³: the 4th-eigenvector extent is 4.6e-3 DU and the 2nd-eigenvector extent is 1.98e-2 DU, roughly 5-20× the validated deviation range. Thus the paper's statement that J* 'falls in the range of costs with <1% error' is not supported by the shown data, and the extents/orderings in Table 1 are extrapolations. Nonlinear sensitivities are generically direction-dependent, so the cost ranking (e.g., 5th eigenvector = second-most-expensive, used for the perilune conclusion) could change in unvalidated directions.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies forced periodic trajectories in the circular restricted three-body problem (CR3BP) under energy-optimal low-thrust control. Using a linearized optimal-control formulation with a quadratic control cost, the authors derive that the set of initial-state deviations δx0 that can be returned to after one reference period with cost at most J* is a hyperellipsoid (1/2)δx0ᵀE*δx0 ≤ J*, with semi-axes given by the eigenvectors and eigenvalues of E* (Eqs. 20–22). They validate the linear estimate against a Newton-Raphson solver for deviations along one eigenvector (Figure 1), then apply the method to an L2 halo orbit in the Earth-Moon system, computing the semi-axis extents for a representative spacecraft (Table 1) and concluding that decreasing perilune distance is relatively expensive because it corresponds to the second-most-expensive eigenvector direction.","tokens_in":7420,"tokens_out":7642,"duration_ms":71661,"significance":"If the linearized ellipsoidal bound remains accurate at the amplitudes and directions used in Table 1, the paper offers a computationally cheap, systematic way to characterize low-thrust forced periodic trajectories and to rank the cost of state deviations in the CR3BP. The derivation is self-contained and standard: it follows from optimal-control theory, the state transition matrix, and an eigenvalue decomposition, with no fitted parameters other than the energy budget J*. The main strength is the transparent analytical framework, which could be useful for mission design. However, the quantitative conclusions depend on validation that currently covers only a small portion of the relevant state space, so the significance is conditional on the validation being extended.","major_comments":[{"comment":"The validation is performed only for deviations along the 4th eigenvector and for |dx| ≤ 1e-3 DU, corresponding to costs up to approximately 1e-4 DU²/TU³. Table 1, however, lists semi-axis extents for J* = 3.51e-4 DU²/TU³, with the 4th-eigenvector extent at 4.6e-3 DU and the 2nd-eigenvector extent at 1.98e-2 DU. The statement that J* \"falls in the range of costs with <1% error\" is therefore not supported by the displayed data. Please extend the validation to deviations at the semi-axis extents along each eigenvector, or at least along the eigenvectors used for the conclusions, and report the maximum error in J.","section":"Validation (Figure 1)"},{"comment":"The conclusion that reducing perilune distance is relatively expensive relies on both the time history of the 5th eigenvector and on the eigenvalue ordering that places it as the second-most-expensive direction. Nonlinear corrections to the cost are generically direction-dependent, and the error at the semi-axis amplitudes could differ between the 4th and 5th eigenvectors. Please verify the cost ranking at J* by computing nonlinear optimal costs for initial deviations along each of the five finite eigenvectors at their semi-axis extents.","section":"Investigation and Analysis (Table 1, Figure 3)"},{"comment":"The computation of E* requires inverting the 6x6 block Φx_λ. For a periodic reference orbit this block is related to the monodromy matrix and may be ill-conditioned due to the unit eigenvalue associated with the along-track direction. The paper does not report a condition number or any regularization. Because the near-zero eigenvalue of E* (first row of Table 1) indicates near-singularity, please include a brief conditioning analysis and discuss how the numerical eigenvalues and eigenvectors of E* are affected.","section":"Methods, Eq. (15)"}],"minor_comments":[{"comment":"Reference [7] contains a typo: \"Reahcable\" should be \"Reachable\"; several references also have extra spacing in \"V ol\" and \"T able\" that should be corrected.","section":"References"},{"comment":"The caption states that deviations are made in the direction of the \"4th eigenvector in Table 1\", but Table 1 lists eigenvectors by row order; this is ambiguous because rows are ordered by extent, not by eigenvector index. Please specify the eigenvector by its components or by a consistent label.","section":"Figure 1 caption"},{"comment":"The phrase \"maximum specific power\" is imprecise: J* has units of DU²/TU³, which are units of specific energy, not power. Consider using \"specific energy per orbit\" or \"control effort\" for clarity.","section":"After Eq. (24)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript reads as a conference-style paper (AAS 24-206) and would benefit from a more detailed validation section and a discussion of conditioning to meet the standard of a journal article. None of the requested changes are outside the scope of a revision; the central derivation is sound but the numerical claims need additional support."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper does one clean thing—it specializes the energy-limited reachable-set formalism from Kulik et al. to forced periodic orbits by imposing δxf = δx0, giving E* = [I I]E[I I]^T and a per-eigenvector cost ranking. That is a real, modest contribution. The derivation is standard and internally consistent, and the Newton-Raphson validation along the 4th eigenvector shows the linear estimate is excellent in that direction.\n\nNew and good: the forced-periodic constraint turns a 12-dimensional boundary-condition reachable set into a 6-dimensional hyperellipsoid in initial-state space; Table 1 is a useful cost-direction catalog for a concrete L2 halo reference; the observation that decreasing perilune distance is expensive is physically interesting and worth checking. The paper is honest about relying on prior work and does not oversell novelty.\n\nSoft spots: the validation is narrower than the conclusions. Figure 1 only tests deviations along the 4th eigenvector, up to |dx| = 1e-3 DU, where the largest computed cost is about 1e-4 DU²/TU³. But Table 1 uses J* = 3.51e-4 DU²/TU³, with semiaxes of 4.6e-3 and 1.98e-2 DU—roughly 5 to 20 times the validated range. So the sentence saying J* 'falls in the range of costs with <1% error' is not supported by the data shown. The cost ordering across eigenvectors is also an extrapolation; nonlinear sensitivities are direction-dependent, and the perilune conclusion rests on the linearized 5th-eigenvector ranking. That may well be right, but it is not validated here. Minor: no code or data included, and conditioning of the Φxλ inversion is not discussed.\n\nBottom line: this is a competent conference-level contribution that deserves referee time in astrodynamics, but the claims about accuracy and the cost ranking should be softened or backed by validation in more directions and at the J* amplitudes. I would not use Table 1 for design yet.","headline":"A clean specialization of reachable-set theory to forced periodic orbits, with a useful cost table—but validation covers only one direction and small amplitudes, so the headline numbers are extrapolations.","tokens_in":7901,"tokens_out":2555,"would_cite":false,"duration_ms":23352,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["70F07","70F15","70Q05"],"pacs":[],"model":"deepseek-v4-flash","headline":"Energy-optimal forced periodic trajectories near the L2 halo orbit form a hyperellipsoid in initial-state space, and lowering perilune distance is among the most expensive deviations.","keywords":["circular restricted three-body problem","low-thrust propulsion","forced periodic trajectories","energy-optimal control","reachable set","halo orbit","Earth-Moon L2","linear analysis"],"falsifier":"Take the same L2 halo reference and solve the full nonlinear energy-optimal two-point boundary value problem for initial deviations along the fifth eigenvector at magnitudes spanning the semiaxis extent; if the true cost deviates from $\\frac{1}{2}\\delta x_0^T E^* \\delta x_0$ by more than a few percent in that direction, the ellipsoid's boundary and the perilune-cost ranking would need to be re-drawn.","tokens_in":6979,"feed_emoji":"🌙","tokens_out":7819,"duration_ms":73320,"temperature":0.7,"pith_summary":"This paper aims to establish that, for low-thrust spacecraft near a periodic orbit in the circular restricted three-body problem, the set of initial states that can be forced back to the same state after one period with bounded energy cost is exactly a six-dimensional ellipsoid computed from linearized optimal control. Working near a natural L2 halo orbit in the Earth-Moon system, the authors show that the linear energy cost of such a forced periodic trajectory is a quadratic form in the initial deviation, and that the reachable set is therefore the hyperellipsoid described by that quadratic form. The eigenvectors of the defining matrix rank the directions of deviation by cost. Under a representative 50 mN, 1000 kg spacecraft thrust constraint, the analysis finds that decreasing perilune distance is relatively expensive, which matters for lunar-observation mission design.","feed_headline":"One ellipsoid prices every halo-orbit deviation","feed_subtitle":"Energy-optimal low-thrust orbits near Earth-Moon L2 form a 6-D ellipsoid; dipping toward perilune is pricey.","key_machinery":"The central object is $E^* = \\begin{bmatrix} I_6 & I_6 \\end{bmatrix} E \\begin{bmatrix} I_6 & I_6 \\end{bmatrix}^T$, where $E$ is built from the augmented state transition matrix $\\Phi(t_f,t_0)$ and the cost kernel $\\int (\\Phi^{\\lambda_v}_y)^T \\Phi^{\\lambda_v}_y \\, dt$ after eliminating the initial costates through the linearized boundary-value constraint. The eigenvalue decomposition of this symmetric positive semidefinite matrix converts the constrained optimal-control problem into ellipsoidal geometry: each eigenvector names a direction of initial-state deviation and each eigenvalue sets that direction's cost. The largest semiaxis, corresponding to the smallest eigenvalue, gives an almost-free along-track direction, while the fifth eigenvector, which lowers perilune distance, has the second-largest eigenvalue and therefore the second-highest cost.","core_discovery":"Under the linearized energy-optimal control model, the paper's central claim is that forced periodic trajectories near a natural periodic reference orbit are characterized by the symmetric matrix $E^*$, so that the cost to begin and end at the same initial deviation $\\delta x_0$ is $J = \\frac{1}{2}\\delta x_0^T E^* \\delta x_0$. The set of states with $J \\le J^*$ is therefore the hyperellipsoid $\\{\\delta x_0 : \\frac{1}{2}\\delta x_0^T E^* \\delta x_0 \\le J^*\\}$ with semiaxes $a_i = \\sqrt{2J^*/\\gamma_i}\\, w_i$ from the eigenpairs $(\\gamma_i, w_i)$ of $E^*$. For the chosen L2 halo orbit in the Earth-Moon system and a representative low-thrust spacecraft, the paper computes the semiaxis extents and shows numerically, for deviations along the fourth eigenvector, that the linear cost estimate tracks the full nonlinear energy-optimal cost to better than $0.1\\%$ in the displayed range. From the eigenvector geometry it concludes that reducing perilune distance is the second-most expensive direction, hence relatively expensive.","pith_inferences":["The validation is confined to one eigenvector direction; extending the nonlinear Newton-Raphson check to the fifth eigenvector and to off-axis combinations would test whether the ellipsoid's extremal extents are as large as predicted.","The theoretically infinite semiaxis along the first eigenvector suggests a phase-shift direction that linear theory treats as cost-free; in practice any real phase change will carry a small nonzero cost that this model misses.","The same $E^*$ construction should transfer to other reference orbits, such as L1 halos, near-rectilinear halo orbits, or L4/L5 families, where the cost ranking of directions may differ; the perilune conclusion is specific to this L2 halo reference.","The cost ordering across eigenvectors could change if the quadratic control cost were replaced with a mass-optimal or minimum-time objective, so the ranking is tied to the energy-optimal formulation."],"forward_implications":["For any thrust-limited spacecraft, the semiaxes of the ellipsoid give an immediate feasibility check for a candidate forced periodic trajectory without solving a two-point boundary value problem.","The eigenvector cost ranking supplies a design rule: cheap deviations lie along the large-semiaxis directions, while small-semiaxis directions, including the perilune-lowering fifth eigenvector, should be reserved for high-priority mission objectives.","Because the costates are not periodic when the position and velocity are periodic, the optimal thrust profile repeats only in the state space, not in the control space, so the control must be re-planned each period.","The same $E^*$ matrix gives a full six-dimensional reachable set, so mission designers can map position and velocity bounds at apolune, perilune, and intermediate phases.","The linear analysis makes the energy-constrained reachable set explicit, turning a family of nonlinear optimal-control problems into a single eigenvalue computation."],"supporting_citations":[{"why":"Supplies the energy-limited reachable set formalism that the paper specializes to forced periodic trajectories.","marker":"[7]"},{"why":"Provides the optimal-control costate equations and two-point boundary value formulation used to derive the cost quadratics.","marker":"[8]"},{"why":"Prior computation of reachable sets for energy-limited low-thrust spacecraft relative motion, extended here to periodic constraints.","marker":"[9]"},{"why":"Prior analysis of low-thrust reachable sets whose ellipsoidal geometry the paper adopts.","marker":"[10]"},{"why":"Provides the state transition matrix propagation approach used to precompute reference-trajectory STMs.","marker":"[11]"}],"fun_headline_variants":["Ellipsoid prices all halo orbit changes","Perilune dip is costly: one ellipsoid shows why","Energy-optimal orbits: cost is a hyperellipsoid","Six-D ellipsoid sizes every orbit deviation","Orbit tweak costs? One ellipsoid answers"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The linearized cost estimate is checked against the true nonlinear optimal cost in only one direction in state space, but the paper's ellipsoid sizes and its conclusion that lowering perilune is expensive assume that the same accuracy holds in every direction and at all deviations up to the energy limit.","fun_headline_variants_meta":{"raw":{"variants":["Ellipsoid prices all halo orbit changes","Perilune dip is costly: one ellipsoid shows why","Energy-optimal orbits: cost is a hyperellipsoid","Six-D ellipsoid sizes every orbit deviation","Orbit tweak costs? One ellipsoid answers"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000431,"raw_usage":{"total_tokens":2174,"prompt_tokens":891,"completion_tokens":1283,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":507,"completion_tokens_details":{"reasoning_tokens":1204}},"tokens_in":507,"tokens_out":1283,"duration_ms":12033,"temperature":1.0,"reasoning_tokens":1204,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T18:19:31.909934+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the same L2 halo reference and solve the full nonlinear energy-optimal two-point boundary value problem for initial deviations along the fifth eigenvector at magnitudes spanning the semiaxis extent; if the true cost deviates from $\\frac{1}{2}\\delta x_0^T E^* \\delta x_0$ by more than a few percent in that direction, the ellipsoid's boundary and the perilune-cost ranking would need to be re-drawn.","supporting_citations":[{"cited_title":"Kulik, M","cited_arxiv_id":null,"evidence_quote":"Provides the optimal-control costate equations and two-point boundary value formulation used to derive the cost quadratics."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Prior computation of reachable sets for energy-limited low-thrust spacecraft relative motion, extended here to periodic constraints."},{"cited_title":"Lee and I","cited_arxiv_id":null,"evidence_quote":"Prior analysis of low-thrust reachable sets whose ellipsoidal geometry the paper adopts."},{"cited_title":"Sun and J","cited_arxiv_id":null,"evidence_quote":"Provides the state transition matrix propagation approach used to precompute reference-trajectory STMs."}],"review_version":1}