Early formation of moons around large trans-Neptunian objects via giant impacts

Hidenori Genda; Ryuki Hyodo; Sota Arakawa

arxiv: 1906.10833 · v1 · pith:WMTAAJGXnew · submitted 2019-06-26 · 🌌 astro-ph.EP · physics.geo-ph

Early formation of moons around large trans-Neptunian objects via giant impacts

Sota Arakawa , Ryuki Hyodo , Hidenori Genda This is my paper

Pith reviewed 2026-05-25 15:30 UTC · model grok-4.3

classification 🌌 astro-ph.EP physics.geo-ph

keywords giant impactstrans-Neptunian objectsTNO satelliteshydrodynamic simulationstidal evolutionmoon formationsolar system formation

0 comments

The pith

Giant impacts formed all moons of large trans-Neptunian objects early on, with the bodies starting molten and later becoming rigid.

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

The paper runs hydrodynamic simulations of giant impacts between large trans-Neptunian objects to test whether such collisions can produce the observed satellite systems. The resulting secondary-to-primary mass ratios cover a broad range that overlaps with real TNO satellites. Matching the current spin and orbital periods plus low eccentricities requires the systems to tidally evolve first as fluid-like bodies and later as rigid ones. This leads to the conclusion that every large TNO satellite formed via giant impact before Neptune migrated outward, at a time when the objects were fully or partially molten.

Core claim

Through systematic hydrodynamic simulations of giant impacts, the secondary-to-primary mass ratio in the resulting satellite systems varies over a wide range that overlaps with observed values. The distributions of spin and orbital periods together with small eccentricities are reproducible only when the bodies undergo tidal evolution that begins while fluid-like and ends while rigid. All satellites of large TNOs therefore originated from giant impacts in the early solar system, prior to Neptune's outward migration, while the parent bodies were fully or partially molten.

What carries the argument

Hydrodynamic simulations of giant impacts that generate variable secondary-to-primary mass ratios, followed by a two-stage tidal evolution model from fluid-like to rigid states.

If this is right

Observed mass ratios among TNO satellites arise directly from the outcomes of giant impacts.
Matching periods and eccentricities requires an initial fluid-like phase of tidal evolution followed by a rigid phase.
All large TNO satellites formed via giant impacts before Neptune's outward migration.
The parent bodies were fully or partially molten during the era of these impacts.

Where Pith is reading between the lines

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

The scenario ties satellite formation timing to the dynamical rearrangement of the outer solar system.
Similar impact-driven formation could apply to other distant bodies that experienced comparable collisions.
Compositional measurements of TNO satellites might show chemical signatures of an early molten state.

Load-bearing premise

The current spin, orbital periods, and eccentricities of these satellite systems can be explained only by tidal evolution that starts fluid-like and ends rigid.

What would settle it

A large TNO satellite system whose mass ratio or spin/orbital properties lie outside the range produced by any combination of the impact simulations and fluid-then-rigid tidal evolution.

Figures

Figures reproduced from arXiv: 1906.10833 by Hidenori Genda, Ryuki Hyodo, Sota Arakawa.

**Figure 1.** Figure 1: Snapshots of a giant impact between two differentiated bodies. Both the target and the impactor were differentiated bodies with ice mass fractions, fice, of 0.5 (50 wt.% ice and 50 wt.% basalt). The target mass, Mtar, and the impactor mass, Mimp, were Mtar = 4 × 1021 kg and Mimp = 2 × 1021 kg, respectively. The impact velocity, vimp, was 1.05vesc, and the impact angle, θimp, was 75° (where 0° represents a … view at source ↗

**Figure 2.** Figure 2: Summary of the range of outcomes for the simulated giant impacts. (a) outcomes for giant impacts between differentiated planetary bodies with masses of 4 × 1021 kg and 2 × 1021 kg. Both the impactor and target were differentiated with a fice of 0.5. (b) outcomes for giant impacts between icy undifferentiated planetary bodies with masses of 4 × 1021 kg and 2 × 1021 kg. The size (large, medium, or small) of … view at source ↗

**Figure 3.** Figure 3: Initial distribution of periapsis distance qini and eccentricity eini before tidal evolution. The size (large, medium, small) of each circle represents γsp (10-1 ≦ γsp for large circles, 10-2 ≦ γsp < 10-1 for medium-sized circles, and 10-3 ≦ γsp < 10-2 for small circles). (a) the case for differentiated bodies. (b) the case for undifferentiated bodies. We then performed semi-analytical tidal evolution calc… view at source ↗

**Figure 4.** Figure 4: Final eccentricity after 4.5 Gyr-tidal evolution (efin) for different times taken for intact moons to become rigid bodies (tfluid). The size of each circle represents the γsp value, as in [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗

read the original abstract

Recent studies have revealed that all large (over 1000 km in diameter) trans-Neptunian objects (TNOs) form satellite systems. Although the largest Plutonian satellite, Charon, is thought to be an intact fragment of an impactor directly formed via a giant impact, whether giant impacts can explain the variations in secondary-to-primary mass ratios and spin/orbital periods among all large TNOs remains to be determined. Here we systematically perform hydrodynamic simulations to investigate satellite formation via giant impacts. We find that the simulated secondary-to-primary mass ratio varies over a wide range, which overlaps with observed mass ratios. We also reveal that the satellite systems' current distribution of spin/orbital periods and small eccentricity can be explained only when their spins and orbits tidally evolve: initially as fluid-like bodies, but finally as rigid bodies. These results suggest that all satellites of large TNOs were formed via giant impacts in the early stage of solar system formation, before the outward migration of Neptune, and that they were fully or partially molten during the giant impact era.

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.

Desk Editor's Note private letter to a colleague

The hydro runs show giant impacts can hit the observed TNO satellite mass ratios, but the 'only when fluid then rigid' tidal claim lacks tests against other sequences.

read the letter

The main point is that a grid of giant impact simulations produces secondary-to-primary mass ratios that overlap the observed range for large TNOs, giving a single formation channel that covers the whole population rather than just Pluto-Charon. The work also ties the systems to an early molten state before Neptune's migration. That overlap is the concrete new output from the runs. The systematic variation of mass ratio, velocity, and angle is the part that moves the field forward from case-by-case studies. It supplies a plausible mechanism without obvious circularity, since the mass ratios come from the hydro code rather than being tuned to the targets. The tidal-evolution sequence is presented as necessary to match periods and low eccentricities, which is a reasonable inference from the data but rests on the assumption that other rheology paths were checked. The stress-test note flags that the abstract does not show explicit comparisons to always-rigid or always-fluid cases, so the uniqueness is not yet demonstrated. If the full methods section includes those runs and resolution tests, the claim holds; otherwise it is only that one path works. The pre-migration timing is stated but follows from the impact epoch rather than new evidence in the simulations. This paper is for people working on outer solar system satellite formation and TNO collisional histories. A reader who needs quantitative impact outcomes will get usable numbers from the grid. It deserves a serious referee because the core simulation result is new and addresses an open question, even if the tidal section needs tightening on alternatives.

Referee Report

2 major / 1 minor

Summary. The manuscript reports results from hydrodynamic simulations of giant impacts between large trans-Neptunian objects. These simulations produce secondary-to-primary mass ratios spanning a wide range that overlaps observed values for TNO satellite systems. The authors further argue that the observed distribution of spin/orbital periods together with low eccentricities can be reproduced only via a specific two-stage tidal evolution in which the bodies first behave as fluid-like and later as rigid. This leads to the claim that all satellites of large TNOs formed via giant impacts early in solar-system history, prior to Neptune's outward migration, while the bodies were fully or partially molten.

Significance. If the central results hold, the work supplies a single formation channel that accounts for the satellite systems of every large TNO, tying their origin to the giant-impact epoch and implying a molten state at that time. It thereby links satellite properties directly to the dynamical architecture of the outer solar system before Neptune's migration.

major comments (2)

[Abstract] Abstract: the assertion that the observed periods and eccentricities 'can be explained only when' spins and orbits tidally evolve first as fluid-like bodies and finally as rigid bodies is load-bearing for the interpretation, yet the manuscript provides no explicit comparisons against alternative rheology sequences (always-rigid, always-fluid, or different transition epochs).
[Abstract] Abstract: the conclusion that 'all satellites of large TNOs were formed via giant impacts' rests on both the mass-ratio overlap and the uniqueness of the tidal path; the reported overlap is described only qualitatively ('varies over a wide range') without quantitative coverage metrics or sensitivity tests to the listed free parameters (impact velocity, angle, mass ratio, initial thermal/spin states).

minor comments (1)

The abstract refers to 'systematic' hydrodynamic simulations but supplies no reference to resolution tests, convergence criteria, or tabulated parameter ranges that would allow independent assessment of robustness.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive report. We address the two major comments below. Both identify legitimate gaps in the current presentation, and we will revise the manuscript to incorporate the requested comparisons and quantitative metrics.

read point-by-point responses

Referee: [Abstract] Abstract: the assertion that the observed periods and eccentricities 'can be explained only when' spins and orbits tidally evolve first as fluid-like bodies and finally as rigid bodies is load-bearing for the interpretation, yet the manuscript provides no explicit comparisons against alternative rheology sequences (always-rigid, always-fluid, or different transition epochs).

Authors: We agree that explicit comparisons to alternative rheology sequences are needed to substantiate the uniqueness claim. The main text contains tidal-evolution calculations for the fluid-to-rigid pathway, but does not show side-by-side results for always-rigid, always-fluid, or alternate transition epochs. In the revised manuscript we will add a dedicated subsection and accompanying figure that directly compares the final spin/orbital periods and eccentricities obtained under each rheology sequence, confirming that only the two-stage fluid-then-rigid evolution reproduces the observed distribution. revision: yes
Referee: [Abstract] Abstract: the conclusion that 'all satellites of large TNOs were formed via giant impacts' rests on both the mass-ratio overlap and the uniqueness of the tidal path; the reported overlap is described only qualitatively ('varies over a wide range') without quantitative coverage metrics or sensitivity tests to the listed free parameters (impact velocity, angle, mass ratio, initial thermal/spin states).

Authors: The referee is correct that the abstract and results section describe the mass-ratio overlap only qualitatively. Although the simulation suite varies the listed parameters, we do not report quantitative coverage fractions or systematic sensitivity tests. We will revise the results section to include (i) histograms or cumulative distribution functions of secondary-to-primary mass ratios across the full parameter grid, (ii) the fraction of runs that fall within the observed TNO range, and (iii) sensitivity plots showing how the overlap changes with impact velocity, angle, mass ratio, and initial thermal/spin state. The abstract will be updated to reflect these quantitative findings. revision: yes

Circularity Check

0 steps flagged

No significant circularity; results from independent hydrodynamic simulations

full rationale

The derivation relies on hydrodynamic impact simulations that generate a range of secondary-to-primary mass ratios overlapping observations, plus separate tidal evolution modeling to match periods and eccentricities. Neither reduces by construction to the target data via definition, fitting, or self-citation chain. The 'only when' phrasing for fluid-then-rigid evolution is a modeling outcome rather than a self-definitional or fitted-input prediction. No load-bearing self-citation or ansatz smuggling is exhibited in the text. This is the common honest finding of a self-contained simulation study.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on the validity of hydrodynamic impact modeling and on a two-stage tidal-evolution model whose parameters are not independently constrained in the abstract; multiple impact conditions (velocity, angle, mass ratio) are varied to produce the reported overlap with observations.

free parameters (2)

Impact velocity, angle, and mass ratio
Varied across simulation runs to generate the reported range of secondary-to-primary mass ratios
Initial thermal and spin states of colliding bodies
Chosen to set the transition from fluid-like to rigid behavior in the tidal-evolution phase

axioms (2)

domain assumption Hydrodynamic equations govern the outcome of giant impacts between icy bodies at the relevant scales
Foundation of all reported simulation results
domain assumption Tidal evolution proceeds first as fluid-like bodies and later as rigid bodies
Required to match the observed spin/orbital periods and eccentricities

pith-pipeline@v0.9.0 · 5728 in / 1599 out tokens · 35552 ms · 2026-05-25T15:30:53.967992+00:00 · methodology

discussion (0)

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/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We also reveal that the satellite systems’ current distribution of spin/orbital periods and small eccentricity can be explained only when their spins and orbits tidally evolve: initially as fluid-like bodies, but finally as rigid bodies.
IndisputableMonolith/Foundation/AbsoluteFloorClosure.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

systematically perform hydrodynamic simulations to investigate satellite formation via giant impacts

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

5 extracted references · 5 canonical work pages · 1 internal anchor

[1]

primary", and the second largest body the

Analysis of the masses of planetary bodies 13 To obtain the masses of planetary bodies from the collision outcome data, we used a friends-of-friends algorithm14 to identify clumps of SPH particles. If the distance between two particles was less than a critical value (lFOF) we defined these particles as belonging to the same clump. lFOF is given by lFOF = ...

work page
[2]

The tidal lag caused by friction leads to angular momentum exchange, which also leads to spin and orbital evolution

Orbital Evolution of Formed Satellites Both the planet and satellite raise tides on each other. The tidal lag caused by friction leads to angular momentum exchange, which also leads to spin and orbital evolution. The tidal evolution of Pluto-Charon has been investigated by several previous studies34,35. In this study, we used O(e6) tidal evolution equatio...

work page
[3]

Monaghan, J. J. 1992, ARA&A, 30, 543. 31. Genda, H., Fujita, T., Kobayashi, H., Tanaka, H., & Abe, Y. 2015, Icarus, 262, 58. 32. Tillotson, J. H. 1962, General Atomic Report, GA-3216. 33. Agnor, C., & Asphaug, E. 2004, ApJL, 613, L157. 34. Ward, W. R., & Canup, R. M. 2006, Science, 313, 1107. 35. Cheng, W. H., Lee, M. H., & Peale, S. J. 2014, Icarus, 233,...

work page 1992
[4]

+ 𝒪(𝑒U)1 − 3𝑛2 𝑀&𝑀- 2𝑅-𝑎5M 𝑘.,-𝑄- Q732𝑒

TIDAL HEATING ON INTACT MOONS When the primary rotates quickly and the secondary is in synchronous state, the orbit-averaged variation of the semimajor axis can be simplified as follows: 1𝑎d𝑎d𝑡 = 3𝑛2 𝑀-𝑀& 2𝑅&𝑎5M 𝑘.,&𝑄& /2 + 27𝑒. + 𝒪(𝑒U)1 − 3𝑛2 𝑀&𝑀- 2𝑅-𝑎5M 𝑘.,-𝑄- Q732𝑒. + 𝒪(𝑒U)W. Here we discuss whether satellites in solid-state can turn into fluid-like bo...

work page
[5]

CATASTROPHIC DISRUPTION OF PRIMORDIAL MOONS AND THE FORMATION OF HAUMEAN MOONS Haumea has a short spin period of 3.92 hours11, and it has two satellites named Hi’iaka and Namaka. The masses of Haumea, Hi’iaka, and Namaka are 4.0 × 1021 kg, 1.8 × 1019 kg, and 1.8 × 1018 kg, respectively12, and the secondary-to-primary mass ratio (γsp) of the Haumean system...

work page internal anchor Pith review Pith/arXiv arXiv 2007

[1] [1]

primary", and the second largest body the

Analysis of the masses of planetary bodies 13 To obtain the masses of planetary bodies from the collision outcome data, we used a friends-of-friends algorithm14 to identify clumps of SPH particles. If the distance between two particles was less than a critical value (lFOF) we defined these particles as belonging to the same clump. lFOF is given by lFOF = ...

work page

[2] [2]

The tidal lag caused by friction leads to angular momentum exchange, which also leads to spin and orbital evolution

Orbital Evolution of Formed Satellites Both the planet and satellite raise tides on each other. The tidal lag caused by friction leads to angular momentum exchange, which also leads to spin and orbital evolution. The tidal evolution of Pluto-Charon has been investigated by several previous studies34,35. In this study, we used O(e6) tidal evolution equatio...

work page

[3] [3]

Monaghan, J. J. 1992, ARA&A, 30, 543. 31. Genda, H., Fujita, T., Kobayashi, H., Tanaka, H., & Abe, Y. 2015, Icarus, 262, 58. 32. Tillotson, J. H. 1962, General Atomic Report, GA-3216. 33. Agnor, C., & Asphaug, E. 2004, ApJL, 613, L157. 34. Ward, W. R., & Canup, R. M. 2006, Science, 313, 1107. 35. Cheng, W. H., Lee, M. H., & Peale, S. J. 2014, Icarus, 233,...

work page 1992

[4] [4]

+ 𝒪(𝑒U)1 − 3𝑛2 𝑀&𝑀- 2𝑅-𝑎5M 𝑘.,-𝑄- Q732𝑒

TIDAL HEATING ON INTACT MOONS When the primary rotates quickly and the secondary is in synchronous state, the orbit-averaged variation of the semimajor axis can be simplified as follows: 1𝑎d𝑎d𝑡 = 3𝑛2 𝑀-𝑀& 2𝑅&𝑎5M 𝑘.,&𝑄& /2 + 27𝑒. + 𝒪(𝑒U)1 − 3𝑛2 𝑀&𝑀- 2𝑅-𝑎5M 𝑘.,-𝑄- Q732𝑒. + 𝒪(𝑒U)W. Here we discuss whether satellites in solid-state can turn into fluid-like bo...

work page

[5] [5]

CATASTROPHIC DISRUPTION OF PRIMORDIAL MOONS AND THE FORMATION OF HAUMEAN MOONS Haumea has a short spin period of 3.92 hours11, and it has two satellites named Hi’iaka and Namaka. The masses of Haumea, Hi’iaka, and Namaka are 4.0 × 1021 kg, 1.8 × 1019 kg, and 1.8 × 1018 kg, respectively12, and the secondary-to-primary mass ratio (γsp) of the Haumean system...

work page internal anchor Pith review Pith/arXiv arXiv 2007