pith. sign in

arxiv: 2603.11787 · v2 · pith:PWAYIGZUnew · submitted 2026-03-12 · 🌌 astro-ph.EP

Equilibrium figure of Haumea and possible detection by stellar occultation

C. Staelen , N. Rambaux , F. Chambat , J. C. Castillo-Rogez This is my paper

Pith reviewed 2026-05-15 12:22 UTC · model grok-4.3

classification 🌌 astro-ph.EP
keywords Haumeastellar occultationhydrostatic equilibriumequilibrium figuredwarf planetinterior structurecritical rotation
0
0 comments X

The pith

Haumea's observed shape matches hydrostatic models near critical rotation that pinch and deviate up to 110 km from an ellipsoid.

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

The paper checks whether the shape of dwarf planet Haumea measured by stellar occultation fits a body in hydrostatic equilibrium with a differentiated interior. Three families of models are tested, each built from a rocky core, an intermediate layer, and a volatile-rich outer shell that may be porous. The BALEINES code computes the equilibrium shapes of the layer boundaries and shows that the configurations closest to the observed outline sit near critical rotation. This produces a pinched figure instead of a clean ellipsoid. The result matters because it shows the occultation data can be explained by hydrostatic balance alone and identifies a specific future observation that can confirm or rule out the pinched geometry.

Core claim

Using the code BALEINES, which solves for the equilibrium figures of the boundaries between layers, we show that the hydrostatic models closest to the shape derived by stellar occultation approach a state of critical rotation, which translates into a pinched shape with large deviations from an ellipsoid (up to 110 km).

What carries the argument

BALEINES code that solves equilibrium figures of the layer boundaries for three-layer interior models of Haumea.

If this is right

  • Existing stellar occultation data and light curves cannot distinguish an ellipsoid from the pinched shape.
  • The pinched figure should be detectable if northern or southern limb chords are obtained during the May 4 2026 stellar occultation.
  • Haumea is close to the critical rotation rate for its mass and density distribution.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Confirmation would imply Haumea's spin is near the breakup limit set by its internal density profile.
  • The same layered modeling approach can be used to interpret shapes of other fast-spinning trans-Neptunian objects.
  • A detected pinch would favor the assumed interior layering over non-hydrostatic explanations such as a recent collision.

Load-bearing premise

Haumea is in hydrostatic equilibrium with one of the three assumed interior structures consisting of a rocky core, an intermediate layer of partial differentiation or organics, and a volatile-rich outer shell that may contain porosity.

What would settle it

Chords recorded on the northern or southern limbs of Haumea's shadow during the May 4 2026 stellar occultation that show no pinched shape or deviations smaller than about 50 km from an ellipsoid.

Figures

Figures reproduced from arXiv: 2603.11787 by C. Staelen, F. Chambat, J. C. Castillo-Rogez, N. Rambaux.

Figure 1
Figure 1. Figure 1: Possible hydrostatic models obtained in this work (pink region) compared with the Jacobi sequence (red line). The models of table 1 are plotted with black crosses; “D19” corresponds to the hydrostatic model preferred by Dunham et al. (2019). The best-fit ellipsoid obtained by Ortiz et al. (2017) is reported as well and plotted with 1, 2 or 3σ uncertainty regions. with these shapes are reported in Appendix … view at source ↗
Figure 3
Figure 3. Figure 3: Projections of the six hydrostatic representative solutions in the sky plane at the epoch of the 2017 occultation. u and v are the celestial East and North, respectively. ing a relative difference of roughly 5 %. In the formation sce￾nario hypothesized by Noviello et al. (2022), Haumea’s rotation was slowed through hydration of the rocky core by hydrother￾mal circulation that generated a redistribution of … view at source ↗
Figure 2
Figure 2. Figure 2: Hydrostatic shape models for Haumea in the (xOz) (left) and (yOz) (right) planes, corresponding to the configurations reported in ta￾ble 1. The dashed lines correspond to the ellipsoids with the same axis lengths as the layers’ boundary. 4.2. Haumea’s shape as a fossil figure As shown by configurations A, B and D, some models are fully compatible with the 2017 occultation and the observed rota￾tion period … view at source ↗
Figure 4
Figure 4. Figure 4: Projections of the hydrostatic shape of configuration C in the sky plane at the time of the 2017 occultation (left) and the future 2026 occultation (right). The projection of the ellipsoidal shape derived by Ortiz et al. (2017) is shown in white dotted line for comparison. The facet colours code the relative irradiance they receive, red being the highest and purple-blue the lowest. occultation, this would … view at source ↗
read the original abstract

The equilibrium figure of dwarf planet Haumea is studied to determine if the observed shape is compatible with a differentiated hydrostatic body. Three groups of interior models of Haumea are assumed, all with a rocky core and a volatile-rich outer shell that may contain some porosity. A third layer located between the core and the outer shell has a density suggesting partial differentiation or the presence of a large fraction of organic matter. Using the code BALEINES, which solves for the equilibrium figures of the boundaries between layers, we show that the hydrostatic models closest to the shape derived by stellar occultation approach a state of critical rotation, which translates into a pinched shape with large deviations from an ellipsoid (up to 110 km). The previous stellar occultation and light curves cannot distinguish between the ellipsoid and the pinched shape, but we predict this figure could be observable on the next stellar occultation of Haumea on May 4, 2026, if some chords are obtained in the northern or southern limbs of the shadow.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript examines whether the shape of dwarf planet Haumea derived from stellar occultation is consistent with hydrostatic equilibrium. It assumes three groups of three-layer interior models (rocky core, intermediate layer of uncertain composition, and volatile-rich outer shell that may include porosity) and uses the BALEINES code to compute equilibrium figures of the layer boundaries. The central result is that the models closest to the observed shape approach critical rotation, producing a pinched non-ellipsoidal figure with deviations up to 110 km from an ellipsoid; the authors predict this feature could be detected during the May 2026 occultation if northern or southern limb chords are obtained.

Significance. If the modeling holds, the work strengthens the case for differentiated hydrostatic equilibrium in Haumea and supplies a concrete, falsifiable prediction for the 2026 occultation. The use of a dedicated multi-layer equilibrium code (BALEINES) to generate non-ellipsoidal figures is a technical strength that goes beyond simple ellipsoidal approximations common in the literature.

major comments (2)
  1. [§3] §3 (Interior models and BALEINES implementation): The manuscript states that BALEINES solves for equilibrium boundaries but provides no explicit equations, convergence criteria, or numerical scheme. Without these details it is impossible to verify that the reported approach to critical rotation and the 110 km deviation are not artifacts of the discretization or boundary conditions.
  2. [Table 2] Table 2 (or equivalent parameter table): The ranges and specific values chosen for the three model groups (core density, intermediate-layer density and thickness, outer-shell porosity) are not tabulated with uncertainties or sensitivity tests. The claim that only models near critical rotation match the occultation shape therefore rests on an incompletely documented parameter search.
minor comments (2)
  1. [Abstract] The abstract and introduction use “pinched shape” without a quantitative definition (e.g., the maximum radial deviation or the flattening parameter at which the equatorial profile ceases to be convex).
  2. [Results] No error budget is given for the 110 km deviation figure; it is unclear whether this is the maximum, rms, or peak-to-peak value relative to the best-fit ellipsoid.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript and the recommendation for major revision. We address the two major comments below and will incorporate the requested improvements to enhance reproducibility and documentation.

read point-by-point responses

  1. Referee: §3 (Interior models and BALEINES implementation): The manuscript states that BALEINES solves for equilibrium boundaries but provides no explicit equations, convergence criteria, or numerical scheme. Without these details it is impossible to verify that the reported approach to critical rotation and the 110 km deviation are not artifacts of the discretization or boundary conditions.

    Authors: We agree that additional technical details on BALEINES are required. In the revised manuscript we will expand §3 with a new subsection that presents the governing equations (generalized Clairaut equation for multi-layer bodies under hydrostatic equilibrium), the numerical scheme (spectral expansion in spherical harmonics up to degree 20 with finite-difference radial discretization), and the convergence criteria (potential residual < 10^{-5} and shape change < 1 km between iterations). We will also note that BALEINES has been benchmarked against analytic ellipsoidal solutions for uniform-density bodies, confirming that the reported 110 km deviation arises from the approach to critical rotation rather than discretization artifacts. revision: yes

  2. Referee: Table 2 (or equivalent parameter table): The ranges and specific values chosen for the three model groups (core density, intermediate-layer density and thickness, outer-shell porosity) are not tabulated with uncertainties or sensitivity tests. The claim that only models near critical rotation match the occultation shape therefore rests on an incompletely documented parameter search.

    Authors: We acknowledge that the parameter exploration needs fuller documentation. We will revise Table 2 to list the full ranges explored for each parameter (core density 2.5–3.5 g cm^{-3}, intermediate-layer density 1.8–2.4 g cm^{-3} and thickness 50–150 km, outer-shell porosity 0–30 %), together with the specific values of the best-matching models. We will also add a sensitivity-analysis paragraph and an accompanying figure that quantifies how the equatorial pinch deviation varies with rotation rate, showing that deviations remain >50 km unless the spin period is within ~5 % of the critical value, thereby reinforcing that only near-critical models reproduce the occultation-derived shape. revision: yes

Circularity Check

0 steps flagged

No significant circularity in the hydrostatic modeling chain

full rationale

The derivation assumes three-layer interior density profiles as inputs, invokes the external BALEINES solver to compute layer-boundary equilibrium figures under hydrostatic equilibrium and rotation, then compares the resulting shapes to the independently measured occultation-derived figure of Haumea. Compatible models are shown to lie near critical rotation and exhibit a pinched non-ellipsoidal deviation; the 2026 occultation is offered as an independent future test. No equation reduces the output shape to a re-expression of the input densities or occultation data by construction, no self-citation supplies a load-bearing uniqueness theorem, and the central result remains falsifiable against new observations without circular reduction.

Axiom & Free-Parameter Ledger

3 free parameters · 2 axioms · 0 invented entities

The central claim rests on assumed layered densities and hydrostatic equilibrium without independent verification of those structures for Haumea.

free parameters (3)
  • Layer densities = not specified
    Densities for rocky core, intermediate layer, and volatile-rich shell are selected to match mass and observed shape.
  • Layer thicknesses or radii = not specified
    Boundaries between the three layers are adjusted within the models.
  • Porosity in outer shell = not specified
    Porosity fraction is included as a variable in the volatile-rich layer.
axioms (2)
  • domain assumption Haumea is in hydrostatic equilibrium
    Required for solving equilibrium figures of layer boundaries with the BALEINES code.
  • domain assumption Interior structure consists of three layers with specified density ranges
    Rocky core, intermediate layer indicating partial differentiation or organics, and volatile-rich outer shell possibly porous.

pith-pipeline@v0.9.0 · 5488 in / 1512 out tokens · 45161 ms · 2026-05-15T12:22:06.288804+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

  • IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    Using the code BALEINES, which solves for the equilibrium figures of the boundaries between layers, we show that the hydrostatic models closest to the shape derived by stellar occultation approach a state of critical rotation, which translates into a pinched shape with large deviations from an ellipsoid (up to 110 km).

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

23 extracted references · 23 canonical work pages

  1. [1]

    2024, Nature Communications, 15, 6202

    Barnouin, O., Ballouz, R.-L., Marchi, S., et al. 2024, Nature Communications, 15, 6202

    work page 2024

  2. [2]

    Barr, A. C. & Schwamb, M. E. 2016, MNRAS, 460, 1542 Brozovi´c, M. & Jacobson, R. A. 2024, AJ, 167, 256

    work page 2016

  3. [3]

    2023, Journal of Geophys- ical Research: Planets, 128, e2022JE007432

    Castillo-Rogez, J., Weiss, B., Beddingfield, C., et al. 2023, Journal of Geophys- ical Research: Planets, 128, e2022JE007432

    work page 2023

  4. [4]

    2023, A&A, 675, L2

    Chambat, F., Vermersch, J., & Rambaux, N. 2023, A&A, 675, L2

    work page 2023

  5. [5]

    T., Desch, S

    Dunham, E. T., Desch, S. J., & Probst, L. 2019, ApJ, 877, 41

    work page 2019

  6. [6]

    J., Grundy, W

    Holler, B. J., Grundy, W. M., Buie, M. W., & Noll, K. S. 2021, Icarus, 355, 114130

    work page 2021

  7. [7]

    L., Gaskell, R

    Jorda, L., Lamy, P. L., Gaskell, R. W., et al. 2012, Icarus, 221, 1089

    work page 2012

  8. [8]

    Kondratyev, B. P. 2016, Ap&SS, 361, 169

    work page 2016

  9. [9]

    Kondratyev, B. P. 2024, Astronomy Reports, 68, 1449

    work page 2024

  10. [10]

    2008, AJ, 135, 1749

    Lacerda, P., Jewitt, D., & Peixinho, N. 2008, AJ, 135, 1749

    work page 2008

  11. [11]

    1982, Deformations of an Elastic Earth (New York: Academic Press)

    Lanzano, P. 1982, Deformations of an Elastic Earth (New York: Academic Press)

    work page 1982

  12. [12]

    C., Brown, M

    Lockwood, A. C., Brown, M. E., & Stansberry, J. 2014, Earth Moon and Planets, 111, 127

    work page 2014

  13. [13]

    & Meynet, G

    Maeder, A. & Meynet, G. 2012, Rev. Mod. Phys., 84, 25

    work page 2012

  14. [14]

    B., Glein, C

    McKinnon, W. B., Glein, C. R., Bertrand, T., & Rhoden, A. R. 2020, arXiv preprint arXiv:2011.14030

    work page arXiv 2020

  15. [15]

    Murray, C. D. & Dermott, S. F. 1999, Solar System Dynamics (Cambridge: Cam- bridge University Press)

    work page 1999

  16. [16]

    L., Desch, S

    Noviello, J. L., Desch, S. J., Neveu, M., Proudfoot, B. C. N., & Sonnett, S. 2022, The Planetary Science Journal, 3, 225

    work page 2022

  17. [17]

    L., Santos-Sanz, P., Sicardy, B., et al

    Ortiz, J. L., Santos-Sanz, P., Sicardy, B., et al. 2017, Nature, 550, 219

    work page 2017

  18. [18]

    L., Thirouin, A., Campo Bagatin, A., et al

    Ortiz, J. L., Thirouin, A., Campo Bagatin, A., et al. 2012, MNRAS, 419, 2315

    work page 2012

  19. [19]

    Proudfoot, B. C. N., Ragozzine, D. A., Giforos, W., et al. 2024, The Planetary Science Journal, 5, 69

    work page 2024

  20. [20]

    L., Ortiz, J

    Rizos, J. L., Ortiz, J. L., Gutierrez, P. J., Navajas, I. M., & Lara, L. M. 2025, A&A, 703, A178

    work page 2025

  21. [21]

    & Huré, J.-M

    Staelen, C. & Huré, J.-M. 2026, A&A, accepted for publication

    work page 2026

  22. [22]

    2025, Philos

    Vermersch, J., Rambaux, N., & Sicardy, B. 2025, Philos. Trans. A Math. Phys. Eng. Sci., 383, 20240186

    work page 2025

  23. [23]

    Zharkov, V . N. & Trubitsyn, V . P. 1978, Physics of planetary interiors (Tucson: Pachart Publishing House) Article number, page 4 of 5 C. Staelen et al.: Equilibrium figure of Haumea and possible detection by stellar occultation Appendix A: Light curves The Hapke model was used for the reflectance of Haumea, with the same parameters as Lockwood et al. (2...

    work page 1978