{"id":"3f73ffea-d8f7-457d-80a4-ed8a30aee352","arxiv_id":"2411.14240","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":2,"one_line_summary":"A closed-form potential for a straight-segment asteroid model with linear density is proposed and used to find periodic orbits, but the potential formula has a sign inconsistency for asymmetric densities.","lead":"This paper derives a closed-form gravitational potential for a test particle near an elongated asteroid modeled as a straight segment with density varying linearly along its length, then uses it to find periodic and quasi-periodic orbits. The potential formula appears to contain an endpoint-labeling sign error for asymmetric densities, so the orbit results may not describe the intended physical model.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The closed-form potential Eq. (4) does not follow from Eq. (3) under the endpoint labels stated in Eq. (1): the parametrization puts s=0 at the left endpoint while c3,c4 encode the right endpoint, so the Hamiltonian and orbit theorems describe a different potential.","rationale":"The paper's central contribution is the closed-form potential (4) for a straight segment with linearly varying density and the Hamiltonian dynamics derived from it. The load-bearing condition is that Eq. (4) is the actual value of the integral in Eq. (3) under the definitions in Eq. (1). That condition fails: the parametrization and the denominator coefficients label the endpoints inconsistently. The c1=0 limit matches the known homogeneous potential, which validates only the logarithmic part and hides the problem in the novel asymmetric term. Numerical quadrature confirms the mismatch for the explicit example above. Because Eqs. (7)-(8), Theorems 4.1-4.2, and the Poincare sections are derived from Eq. (4), all downstream dynamical results inherit this defect. A secondary gap is the unproved numerical assertion that d+(s;A) remains in (0,2), but that matters only after the potential is corrected. The reader's weakest_assumption identifies the same endpoint-labeling issue, so I agree with it. The paper includes no machine-checked proofs or reproducible code that would independently corroborate the formulas. The verdict should remain REJECT.","tokens_in":9339,"tokens_out":11817,"duration_ms":97916,"concrete_test":"Independently evaluate the line integral in Eq. (3) using u=2Ls-L-cbar and the endpoint distances exactly as defined in Eq. (1), by symbolic integration or high-precision quadrature, and compare it with Eq. (4). A minimal spot check: set L=1, alpha=0.2, beta=1, P=(0.5,1,0); the true integral is approximately -1.706, whereas Eq. (4) with the Eq. (1) labels gives approximately -1.645. If this mismatch reproduces, the potential, Hamiltonian, and all orbit computations in Sections 3-5 must be redone before the central claims can be assessed.","verdict_should_be":"REJECT","load_bearing_attack":"Eq. (3) is obtained from Eq. (2) by the change u=2Ls-L-cbar, so s=0 corresponds to the left endpoint u=-L-cbar and s=1 to the right endpoint u=L-cbar. Under Eq. (1), however, r1 is the distance to the right endpoint E1=(L,0,0): the expression (-L+u+cbar)^2 vanishes at u=L-cbar. Likewise, r2 is the distance to the left endpoint. Therefore the denominator at s=0 must be r2, not r1, so the correct coefficients are c4=r2^2/(4L^2) and c3=(-4L^2-r2^2+r1^2)/(4L^2), the negative of the c3 printed in Eq. (3). Direct quadrature for L=1, alpha=0.2, beta=1, P=(0.5,1,0) gives the true line integral as approximately -1.706, while Eq. (4) with the Eq. (1) definitions of r1,r2 gives approximately -1.645. This is not a small numerical error; it is a structural sign/labeling failure of the asymmetric term. The homogeneous limit c1=0 matches the known segment potential and masks the error. Since the Hamiltonian (7)-(8), Theorems 4.1-4.2, and the Poincare sections in Section 5 are all built on Eq. (4), they concern a gravitational field different from the stated linear-density segment.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper models an elongated asteroid as a non-homogeneous straight segment with linearly varying density, derives a claimed closed-form gravitational potential, and uses it to formulate a Hamiltonian system. It proves the existence of a branch of circular orbits for small density-asymmetry parameter A via the implicit function theorem and reports a larger family of circular orbits for all A, together with Poincaré sections that show quasi-periodic orbits.","tokens_in":9668,"tokens_out":11670,"duration_ms":96561,"significance":"The model is potentially useful because it offers an explicit, fast-to-evaluate potential for asymmetric elongated bodies, complementing the constant-density segment and the dipole model. The paper correctly recovers the known constant-density limit and presents a compact formulation in terms of the variables s and d. However, the central potential formula is currently derived with inconsistent endpoint labeling, so the dynamical conclusions as written do not apply to the stated physical model; furthermore, the global existence claim for circular orbits rests on unproved numerical evidence.","major_comments":[{"comment":"The closed-form potential (4) does not follow from the integral (3) under the endpoint definitions in (1). The change of variables u=2Ls-L-cbar puts s=0 at the left endpoint u=-L-cbar and s=1 at the right endpoint u=L-cbar, while Eq. (1) defines r1 and r2 as the distances to the right and left endpoints, respectively. Therefore the denominator at s=0 must be r2, not r1. The printed coefficients c4=r1^2/(4L^2) and c3=(-4L^2-r1^2+r2^2)/(4L^2) are those appropriate to starting at the right endpoint. Direct numerical quadrature for L=1, alpha=0.2, beta=1, P=(0.5,1,0) gives the line integral (3) as approximately -1.706, while Eq. (4) with the definitions in (1) gives approximately -1.645. Since the Hamiltonian (7)-(8), the equations of motion (10)-(13), and the results of Sections 4 and 5 are all built on Eq. (4), they describe a different gravitational field from the stated linear-density segment.","section":"Eqs. (3)-(4), Section 2"},{"comment":"The claim that a unique circular orbit exists for every s* is not rigorously established. The paper states 'we have numerical evidence' that d+(s;A) lies in (0,2) and d-(s;A) lies in (-infinity,-2) for all s in (2,infinity) and A in [0,1/3), but no proof is supplied. The proof of Theorem 4.2 begins with 'Assuming the hypothesis d=d* in (0,2)', which is precisely the assumption that needs to be proven. The two limits at s=2 and infinity do not by themselves guarantee the bound on the whole interval. Without this bound, the existence and uniqueness of the circular orbit family is conditional on numerical observation rather than a theorem.","section":"Theorem 4.2, Section 4"}],"minor_comments":[{"comment":"The statement announces a 'unique circular orbit', but the proof obtains two possibilities, c* and -c*. Clarify whether uniqueness is up to the sign of c.","section":"Theorem 4.2, Section 4"},{"comment":"The strict inequalities -M/(2L^2) < alpha < M/(2L^2) are stated, but if zero density at one endpoint is physically admissible, the closed interval is the correct condition; please clarify.","section":"Proposition 2.1(iii), Section 2"},{"comment":"The axes and section conditions are not always clear in the Poincaré sections; please add axis labels and explicitly state the section plane and direction in each caption.","section":"Figures 4-7, Section 5"},{"comment":"The sentence stating that d'(0) is 'strictly decreasing in x' is ambiguous about the direction of monotonicity; rephrase to indicate that d'(0) is positive and decreases to 0 as x increases.","section":"Remark 2, Eq. (19)"}],"recommendation":"major_revision","confidential_remarks":"The endpoint-labeling error in Eq. (4) is a genuine, load-bearing mistake, but it appears correctable by relabeling endpoints and recomputing the derived formulas. The authors should also supply a rigorous proof of the bound on d+(s;A) if they wish to keep Theorem 4.2 in its current form. I recommend major revision rather than rejection because the overall approach is promising and the specific errors are identifiable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe linear-density segment model is a good idea, but the central formula is wrong. The potential in Eq. (4) does not follow from the integral in Eq. (3) with the endpoint labels as defined. The fault is a sign error in c3. With s=0 at the left endpoint (E2), the denominator squared at s=0 must be r2^2, not r1^2, so c4=r2^2/(4L^2) and c3=(r1^2-r2^2-4L^2)/(4L^2). The printed c3 has the opposite sign on (r2^2-r1^2). For a concrete check, take L=1, α=0.2, β=1, P=(0.5,1,0): the true line integral is about -1.706, while Eq. (4) gives about -1.645. The homogeneous limit c1=0 masks the error because the offending term vanishes, which is why the A=0 results still match [8,9].\n\nThis is not a minor typo. The Hamiltonian (7)-(8), Theorems 4.1-4.2, and the Poincaré sections all describe a different potential than the stated linear-density segment. Since the whole point of the paper is the asymmetric A≠0 case, the load-bearing result is invalid.\n\nWhat is good: the idea is a natural and potentially useful extension of the constant-density segment, and the paper is honest in flagging the numerical range for d+(s;A) as evidence rather than proof. That secondary gap, plus the approximate initial conditions for the quasi-periodic orbits, would be fixable with more work. But the sign error needs to be corrected and all computations redone.\n\nWho this is for: the asteroid-dynamics and celestial mechanics crowd would value a corrected version; the closed-form potential would be a fast alternative to polyhedral and mascon models. As it stands, I would not cite it. Still, it deserves peer review: the error is subtle and the framework is otherwise coherent, so a referee can give clear, actionable revision guidance.","headline":"The linear-density segment potential has a sign error in c3 that invalidates the A≠0 dynamics; the homogeneous limit is fine.","tokens_in":10204,"tokens_out":15627,"would_cite":false,"duration_ms":111023,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["70F15","37N05","70F16"],"pacs":[],"model":"deepseek-v4-flash","headline":"Linear-density segment yields a closed-form gravity model for elongated asteroids","keywords":["straight-segment model","linear density","closed-form gravitational potential","Hamiltonian dynamics","circular orbits","quasi-periodic orbits","Poincaré sections","elongated asteroids"],"falsifier":"Evaluate the integral in Eq. (3) numerically at a generic off-axis point using $r_1$ and $r_2$ exactly as defined in Eq. (1), and compare it with the closed form in Eq. (4); a disagreement would show the formula as written is not the potential. Independently, scan $d_+(s;A)$ over $s\\in(2,\\infty)$ and $A\\in[0,1/3)$; if the branch leaves $(0,2)$, the uniqueness argument for circular orbits breaks down.","tokens_in":9052,"feed_emoji":"☄️","tokens_out":7525,"duration_ms":66418,"temperature":0.7,"pith_summary":"This paper tries to establish that the gravitational field of an irregular elongated asteroid can be modeled by a non-homogeneous straight segment whose linear density varies along its length, and that this model has a closed-form potential. From that potential the authors build a Hamiltonian for a test particle and reduce it, using axial symmetry, to two degrees of freedom plus a conserved angular momentum. They prove the existence of circular orbits parameterized by angular momentum and display quasi-periodic orbits obtained from Poincaré sections. If the derivation holds, the model would give a fast, analytic alternative to polyhedral or mascon representations for elongated bodies whose mass distribution is asymmetric.","feed_headline":"Linear density gives closed-form asteroid gravity","feed_subtitle":"One asymmetry parameter reproduces circular and quasi-periodic orbits near elongated bodies.","key_machinery":"The load-bearing machinery is the change of variables $s=R_1+R_2$ and $d=R_1-R_2$, where $R_1$ and $R_2$ are the scaled distances from the particle to the two endpoints. In these variables the potential becomes $U(Q;A)=3Ad-\\frac{1}{4}(3Ads+4)\\ln\\frac{s+2}{s-2}$, the equations of motion split cleanly, and the axial symmetry of the segment makes the polar angle cyclic so $P_\\theta=c$ is a conserved parameter. Critical points of the reduced $(r,x)$ flow are roots of two algebraic equations $F_1=0$ and $F_2=0$ in $(s,d)$; Theorem 4.1 applies the implicit function theorem at the known $A=0$ circular orbit $s_0$, and Theorem 4.2 uses the numerically observed branch $d_+(s;A)\\in(0,2)$ to extend the result. Poincaré sections at $x=0$, $P_x>0$ then convert the search for quasi-periodic orbits into a two-dimensional map study.","core_discovery":"The central discovery is that a segment with density $\\sigma(x)=\\alpha x+\\beta$ produces a potential $V$ that depends only on the two endpoint distances $r_1$ and $r_2$, with the closed form in Eq. (4); after a symplectic rescaling the dynamics reduces to a one-parameter Hamiltonian $H(Q,P;A)$ with $0\\le A\\le 1/3$. The paper claims that for each fixed $s^*=R_1+R_2$ there is a circular orbit (with prograde and retrograde copies) whose position shifts along the segment as $A$ grows, in contrast to the constant-density case where the orbit lies in the perpendicular plane. Theorems 4.1 and 4.2 establish this circular family, and the Poincaré sections of Section 5 identify reduced-periodic orbits that lift to quasi-periodic orbits of the full system. The authors locate the model as the linear-density counterpart to existing constant-density and quadratic-density segment models, and as one of only two closed-form options that can represent asymmetric mass distributions.","pith_inferences":["A consistently labeled version of the closed-form potential would likely extend to piecewise-linear densities on a chain of segments, giving a multi-segment closed-form model for more complex asteroid shapes.","The one-parameter family makes an inverse problem natural: fit $A$, $M$, and $L$ to observed orbital data around an elongated asteroid, then test whether residuals are compatible with this model.","The Poincaré sections suggest that low-angular-momentum orbits are largely chaotic, so practical station-keeping around an asymmetric elongated body would be safest at high $P_\\theta$.","Comparing the linear-density segment's equipotentials with a polyhedral model of a specific asteroid would quantify the error introduced by the linear-density assumption; the paper leaves this quantitative calibration to future work."],"forward_implications":["Asymmetric elongated asteroids can be modeled with one extra parameter $A$ without leaving the realm of closed-form potentials, making the model cheap to evaluate for orbit computations.","For each angular momentum $c$ the reduced system has a circular solution; lifting it gives bounded, roughly circular orbits around bodies whose mass distribution is lopsided.","Reduced-periodic orbits found in the Poincaré sections become quasi-periodic orbits in three dimensions, so the model predicts bounded trajectories near the segment over many revolutions.","The $A=0$ limit reproduces the known constant-density segment results, so the linear-density model is a direct generalization that can be checked against existing calculations.","The parameter $A$ controls the shift of the circular orbit off the perpendicular symmetry plane, giving a measurable signature of mass asymmetry for mission design or remote sensing."],"supporting_citations":[{"why":"Supplies the constant-density closed-form potential and the $A=0$ circular orbit used as the base point for the implicit-function-theorem argument.","marker":"[8]"},{"why":"Provides the Poincaré-section methodology and the $A=0$ quasi-periodic orbit results that Section 5 reproduces and generalizes.","marker":"[9]"},{"why":"Gives the classical closed-form potential for a homogeneous straight segment, the $c_1=0$ limit of Eq. (4).","marker":"[1]"},{"why":"Is the standard potential-theory reference underpinning the line-integral derivation of the potential.","marker":"[3]"},{"why":"Is the classical potential-theory text used for the straight-segment potential.","marker":"[4]"},{"why":"Models the quadratic-density segment whose symmetry limitation motivates the linear-density model.","marker":"[6]"},{"why":"Is the dipole-segment model for axisymmetric elongated asteroids, the alternative the paper compares for asymmetric mass distributions.","marker":"[11]"}],"fun_headline_variants":["Linear density segment gives closed-form asteroid gravity","Closed-form potential for asymmetric asteroid segment","Asteroid orbits from linear density: exact solution","Linear density model reveals asteroid orbit families","Asymmetric segment yields exact asteroid dynamics"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument stands on the assumption that the integration variable's endpoints are labeled consistently between the density substitution and the denominator coefficients of Eq. (3), and on the numerically observed fact that the branch $d_+(s;A)$ stays between $0$ and $2$; if either is wrong, the closed-form potential or the uniqueness of the circular orbit fails.","fun_headline_variants_meta":{"raw":{"variants":["Linear density segment gives closed-form asteroid gravity","Closed-form potential for asymmetric asteroid segment","Asteroid orbits from linear density: exact solution","Linear density model reveals asteroid orbit families","Asymmetric segment yields exact asteroid dynamics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000183,"raw_usage":{"total_tokens":1265,"prompt_tokens":846,"completion_tokens":419,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":462,"completion_tokens_details":{"reasoning_tokens":354}},"tokens_in":462,"tokens_out":419,"duration_ms":4347,"temperature":1.0,"reasoning_tokens":354,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:27:21.318343+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the integral in Eq. (3) numerically at a generic off-axis point using $r_1$ and $r_2$ exactly as defined in Eq. (1), and compare it with the closed form in Eq. (4); a disagreement would show the formula as written is not the potential. Independently, scan $d_+(s;A)$ over $s\\in(2,\\infty)$ and $A\\in[0,1/3)$; if the branch leaves $(0,2)$, the uniqueness argument for circular orbits breaks down.","supporting_citations":[{"cited_title":"Riaguas, Thesis doctoral, Universidad de Zaragoza, (1999)","cited_arxiv_id":null,"evidence_quote":"Supplies the constant-density closed-form potential and the $A=0$ circular orbit used as the base point for the implicit-function-theorem argument."},{"cited_title":"Riaguas, A","cited_arxiv_id":null,"evidence_quote":"Provides the Poincaré-section methodology and the $A=0$ quasi-periodic orbit results that Section 5 reproduces and generalizes."},{"cited_title":"Duboshin, On one particular case of the problem of the translational-rotational motion of two bodies, Soviet Astronomy, 3 (1959), p","cited_arxiv_id":null,"evidence_quote":"Gives the classical closed-form potential for a homogeneous straight segment, the $c_1=0$ limit of Eq. (4)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Is the standard potential-theory reference underpinning the line-integral derivation of the potential."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Is the classical potential-theory text used for the straight-segment potential."},{"cited_title":"Najid, E","cited_arxiv_id":null,"evidence_quote":"Models the quadratic-density segment whose symmetry limitation motivates the linear-density model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Is the dipole-segment model for axisymmetric elongated asteroids, the alternative the paper compares for asymmetric mass distributions."}],"review_version":1}