Modeling and dynamics near irregular elongated asteroids

E. Mart\'inez; J.L.Zapata; J.Vidarte

REVIEW 2 major objections 4 minor 11 references

Modeling and dynamics near irregular elongated asteroids

T0 review · 2 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Linear-density segment yields a closed-form gravity model for elongated asteroids

desk verdict The linear-density segment potential has a sign error in c3 that invalidates the A≠0 dynamics; the homogeneous limit is fine. read the letter →

arxiv 2411.14240 v1 pith:UMM3RNBF submitted 2024-11-21 math.DS

E. Martínez , J.Vidarte , J.L.Zapata This is my paper · ORCID

classification math.DS MSC 70F1537N0570F16

keywords straight-segmentmodellineardensityclosed-formgravitationalpotentialHamiltoniandynamicscircularorbitsquasi-periodicPoincarésectionselongatedasteroids

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

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.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Colleague,

The 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].

This 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.

What 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.

Who 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.

Referee Report

2 major / 4 minor

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.

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 (2)

[Eqs. (3)-(4), Section 2] 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.
[Theorem 4.2, Section 4] 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.

minor comments (4)

[Theorem 4.2, Section 4] 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.
[Proposition 2.1(iii), Section 2] 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.
[Figures 4-7, Section 5] 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.
[Remark 2, Eq. (19)] 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.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the closed-form potential is derived from the defining line integral, and the orbit theorems are independent existence results, not fitted predictions.

full rationale

The paper's derivation chain is not circular. The closed-form potential (Eq. 4) is obtained by direct quadrature of the defining line integral (Eq. 2) after the stated substitution u = 2Ls - L - cbar, with coefficients c1...c4 explicitly defined in Eq. (3). No parameter is fitted to data and no quantity is renamed as a prediction. The A=0 limit is checked against the known homogeneous-segment potential from independent works [3, 1, 8], which is legitimate external support. The main orbit results are existence theorems: Theorem 4.1 uses the implicit function theorem around the known A=0 circular orbit of [8], and Theorem 4.2 solves the algebraic equations F1=F2=0 for a given s*. The Poincare sections in Section 5 are simulations of the Hamiltonian (7)-(8), not inputs used to define the potential; thus they cannot make the derivation circular. The numerical assertion that d+(s;A) lies in (0,2) is unsupported in the text, but an unproved numerical bound is a rigor gap, not circular reasoning. Likewise, the endpoint-labeling and sign concern about Eq. (4) raised in the skeptic summary would be a correctness or algebraic-error issue, not a circularity issue, because the paper does not assume the truth of Eq. (4) when deriving it. There are no load-bearing self-citations: the cited prior works [8, 9] are by Riaguas, Elipe, and Lara, not by the present authors. Accordingly, the circularity score is 0.

Assumptions & free parameters 2 free parameters · 4 assumptions · 0 invented entities

No empirical fitting is involved; A and c are model inputs. The analysis leans on standard potential theory and the implicit function theorem, but the circular-orbit theorem relies on an unproved numerical range for d+, and the potential derivation requires an endpoint-labeling consistency that the paper does not satisfy.

free parameters (2)

A (dimensionless density-slope parameter) = 0 <= A <= 1/3
Chosen model input controlling the asymmetry of the linear density; not fitted to any asteroid data but central to the one-parameter Hamiltonian family.
c = P_theta (angular momentum integral) = positive constant, e.g. c = 0.1, 0.35, 0.9 used in figures
Indexes the families of circular and reduced-periodic orbits; an arbitrary conserved constant of the reduced system, not an observable-constrained parameter.

assumptions (4)

standard math The implicit function theorem applies at (s0, 0, A=0) with nonzero Jacobian determinant (Theorem 4.1 proof).
Used to assert a unique circular orbit for sufficiently small A; a standard theorem, but the existence of the local branch does not by itself extend to all A.
domain assumption The density remains positive along the segment, giving -M/(2L^2) < alpha < M/(2L^2) and after scaling 0 <= A <= 1/3.
Physical positivity constraint; assumed throughout and used to bound s+3Ad > 0 in Proposition 3.2.
ad hoc to paper The root d+(s;A) stays in the interval (0,2) for every s in (2, infinity) and every A in (0,1/3).
Theorem 4.2 restricts to d* = d+ and needs this range for r(s,d) to be real and the orbit to be physical. The paper supports it only by limits and numerical plots, not by a proof.
ad hoc to paper The line integral (3) can be evaluated with an internally consistent endpoint labeling for r1, r2, and the parameter s.
The substitution u = 2Ls - L - bar c and the denominator coefficients c3,c4 use opposite conventions for which endpoint sits at s=0; this consistency is load-bearing for Eq. (4).

how reviews work

0 comments

Cite this review

Pith. "Pith review of Modeling and dynamics near irregular elongated asteroids." pith.science (2026). https://pith.science/paper/UMM3RNBF

@misc{pith2026241114240,
  author       = {Pith},
  title        = {Pith review of: Modeling and dynamics near irregular elongated asteroids},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/UMM3RNBF}},
  note         = {Machine review of arXiv:2411.14240}
}

read the original abstract

We investigate the qualitative characteristics of a test particle attracted to an irregular elongated body, modeled as a non-homogeneous straight segment with a variable linear density. By deriving the potential function in closed form, we formulate the Hamiltonian equations of motion for this system. Our analysis reveals a family of periodic circular orbits parameterized by angular momentum. Additionally, we utilize the axial symmetry resulting from rotations around the segment's axis to consider the corresponding reduced system. This approach identifies several reduced-periodic orbits by analyzing appropriate Poincar\'e sections. These periodic orbits are then reconstructed into quasi-periodic orbits within the full dynamical system.

Figures

Figures reproduced from arXiv: 2411.14240 by the authors.

**Figure 1.** The fixed segment in different reference frames. Left: The reference frame {O; x, y, z} with the segment centered on the origin and lying along the x-axis. Right: The new frame {O; u, v, w} results from a translation in the x-axis locating the center of mass at the new origin. To simplify the denominator, we express PQ⃗ = r1+s r12, where r12 = r2−r1 = 2L with s ∈ (0, 1), which leads us to |PQ| ⃗ = p r 2 1 + s 2 r 2 … view at source ↗

**Figure 2.** Left and right figures shows the surfaces d+(s; A) and d−(s; A), for s ∈ (2, 100) and the parameter A ∈ [0, 1/3). Green and blue planes are d = 2 and d = −2 respectively. In the domain of these figures we observe the brown surface, d+(s) ∈ (0, 2) and the yellow plane, d−(s) ∈ (−∞, −2). (17) d±(s; A) = −2 ± r 3As ln s+2 s−2 − 6A 2 + 4 − 36A2 3A ln s+2 s−2 . Note that, since A ∈ [0, 1/3) and s ∈ (2, +∞) the … view at source ↗

**Figure 3.** This figure shows the plot of the periodic orbit, The yellow straight line is a representation of the segment. In contrast with the constant density case A = 0, the periodic orbit does not lie on the ηζ-plane. For A = 0.25, the initial conditions of this orbit are r = 4.8926, θ = 0, x = −0.5042, Pr = 0, Pθ = 3.1023, Px = 0. Theorem 4.2 (Circular orbit). For each fixed s = s ∗ , we consider (s ∗ , d∗ , c∗ ) as they w… view at source ↗

Figures from the paper (6 more)

**Figure 4.** Figure 4: Poincar´e sections for A = 0: From left to right, c = 0.1, 0.35, 0.9. These sections reproduce the analysis from [9] and are constructed using the same units as in that study. Specifically, the length unit is taken as 2L. 5. Quasi-periodic orbits around a non-homogeneo…

**Figure 5.** Figure 5: Poincar´e sections for A ̸= 0: First and second rows with A = 1/8, 1/4, respectively. In both cases, from left to right, c = 0.35, 0.7, 1, 1.8. These sections are constructed using the unit scaling defined in (6). 8 [PITH_FULL_IMAGE:figures/full_fig_p008_5.png]

**Figure 6.** Figure 6: The relative equilibria Ei for i = 1, . . . , 5 correspond with examples of quasiperiodic orbits for A = 1/8. From left to right, we consider c = 0.7 and c = 1 respectively. Reduced-periodic orbits will be detected through suitable Poincar´e sections, allowing us to a…

**Figure 7.** Figure 7: The relative equilibria Ei for i = 6, . . . , 10 correspond with examples of quasi-periodic orbits for A = 1/4. From left to right, we consider c = 0.7 and c = 1 respectively. Here, we compute Poincar´e sections in the plane (r, Pr) for x = 0 and Px > 0 [PITH_FULL_IMA…

**Figure 8.** Figure 8: The first row corresponds to the equilibria E1 = (1.10845, 0, −0.0398045, 1.34386) and the first figure displays 2D projection in the rx-plane, and the following figures are two different views of the 3D orbit in configuration space, illustrating the orbit over several…

**Figure 9.** Figure 9: The first row corresponds to the equilibria E9 = (1.68132, 0, −0.46653, 0.900399). The first figure displays 2D projection in the rxplane, and the following figures are two different views of the 3D orbit in configuration space, illustrating the orbit over several rev…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

11 extracted references · 11 canonical work pages

[1]

Duboshin, On one particular case of the problem of the translational-rotational motion of two bodies, Soviet Astronomy, 3 (1959), p

G. Duboshin, On one particular case of the problem of the translational-rotational motion of two bodies, Soviet Astronomy, 3 (1959), p. 154

work page 1959
[2]

Geissler, J.-M

P. Geissler, J.-M. Petit, D. D. Durda, R. Greenberg, W. Bottke, M. Nolan, and J. Moore , Erosion and ejecta reaccretion on 243 ida and its moon , Icarus, 120 (1996), pp. 140–157

work page 1996
[3]

O. D. Kellogg , Foundations of Potential Theory , Springer, 1967

work page 1967
[4]

MacMillan, The Theory of the Potential , Dover Publications, Inc., New York, 1958

work page 1958
[5]

Najid and E

N.-E. Najid and E. H. Elourabi , Equilibria and stability around a straight rotating segment with a parabolic profile of mass density , The Open Astronomy Journal, 5 (2012), pp. 19–25. 11

work page 2012
[6]

Najid, E

N.-E. Najid, E. H. Elourabi, and M. Zegoumou , Potential generated by a massive inhomogeneous straight segment, Research in Astronomy and Astrophysics, 11 (2011), p. 345

work page 2011
[7]

Prieto-Llanos and M

T. Prieto-Llanos and M. A. Gomez-Tierno , Stationkeeping at libration points of natural elongated bodies, Journal of Guidance, Control, and Dynamics, 17 (1994), pp. 787–794

work page 1994
[8]

Riaguas, Thesis doctoral, Universidad de Zaragoza, (1999)

A. Riaguas, Thesis doctoral, Universidad de Zaragoza, (1999)

work page 1999

Show all 11 references

[9]

Riaguas, A

A. Riaguas, A. Elipe, and M. Lara , Periodic orbits around a massive straight segment , Celestial Mechanics and Dynamical Astronomy, 73 (1999), pp. 169–178

1999
[10]

R. A. Werner and D. J. Scheeres , Exterior gravitation of a polyhedron derived and compared with harmonic and mascon gravitation representations of asteroid 4769 castalia , Celestial Mechanics and Dynamical Astronomy, 65 (1996), pp. 313–344

1996
[11]

X. Zeng, Y. Zhang, Y. Yu, and X. Liu , The dipole segment model for axisymmetrical elongated asteroids, The Astronomical Journal, 155 (2018), p. 85. 1Departamento de Matem´atica, F acultad de Ciencias, Universidad del B´ıo-B´ıo, Casilla 5-C, Concepci´on, Chile 2Departamento de...

2018

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.

[1] [1]

Duboshin, On one particular case of the problem of the translational-rotational motion of two bodies, Soviet Astronomy, 3 (1959), p

G. Duboshin, On one particular case of the problem of the translational-rotational motion of two bodies, Soviet Astronomy, 3 (1959), p. 154

work page 1959

[2] [2]

Geissler, J.-M

P. Geissler, J.-M. Petit, D. D. Durda, R. Greenberg, W. Bottke, M. Nolan, and J. Moore , Erosion and ejecta reaccretion on 243 ida and its moon , Icarus, 120 (1996), pp. 140–157

work page 1996

[3] [3]

O. D. Kellogg , Foundations of Potential Theory , Springer, 1967

work page 1967

[4] [4]

MacMillan, The Theory of the Potential , Dover Publications, Inc., New York, 1958

work page 1958

[5] [5]

Najid and E

N.-E. Najid and E. H. Elourabi , Equilibria and stability around a straight rotating segment with a parabolic profile of mass density , The Open Astronomy Journal, 5 (2012), pp. 19–25. 11

work page 2012

[6] [6]

Najid, E

N.-E. Najid, E. H. Elourabi, and M. Zegoumou , Potential generated by a massive inhomogeneous straight segment, Research in Astronomy and Astrophysics, 11 (2011), p. 345

work page 2011

[7] [7]

Prieto-Llanos and M

T. Prieto-Llanos and M. A. Gomez-Tierno , Stationkeeping at libration points of natural elongated bodies, Journal of Guidance, Control, and Dynamics, 17 (1994), pp. 787–794

work page 1994

[8] [8]

Riaguas, Thesis doctoral, Universidad de Zaragoza, (1999)

A. Riaguas, Thesis doctoral, Universidad de Zaragoza, (1999)

work page 1999