An angular-momentum preserving dissipative model for the point-mass N -body problem

This paper proposes an angular-momentum preserving dissipative model for the N-body problem that reduces to homographic equations for central configurations, analyzes the resulting two-body dynamics via Poincaré compactification, and demonstrates that the dissipation does not affect periapsis precession.

The Gist

Imagine the universe as a giant, cosmic dance floor. Usually, when we think of planets and stars dancing to the tune of gravity, we imagine a perfect, frictionless ballroom where the dancers never get tired, never slow down, and their energy is perfectly conserved. This is the "conservative" view of the universe.

But in reality, space isn't empty or frictionless. There are tides, gas clouds, and internal friction that act like a subtle drag on the dancers. This paper proposes a new, simplified way to model that drag.

Here is the breakdown of what the authors did, using everyday analogies:

1. The New Rule: "The Energy Siphon"

In most models of space friction (like a car driving through mud), the drag slows things down and messes up their spin. If you push a spinning top while dragging it, it might wobble and fall over.

The authors created a special "magic friction" for their model. Imagine a dancer who is being slowed down by a gentle hand on their shoulder. This hand removes their energy (they get tired and slow down), but it is very careful not to touch their spin.

• The Result: The dancer loses speed and falls inward, but they keep their total spin (angular momentum) exactly the same.
• Why it matters: In the real world, tidal forces (like the Moon pulling on Earth's oceans) do exactly this. They drain energy but keep the total spin of the Earth-Moon system constant.

2. The "Magic Number" (The Special Case)

The authors tested different ways this friction could work. They found that if the friction gets stronger in a very specific way as objects get closer (specifically, if the friction depends on distance in a "cubic" way, or $d=3$), something magical happens.

The Analogy: Imagine a complex dance routine involving 100 people (the N-body problem). Usually, predicting how 100 people move while holding hands is a nightmare. But with this specific "magic friction," the whole group suddenly starts moving like a single, simple pair of dancers.

• The Math: The complex equations for 100 bodies collapse into the same simple equations we use for just two bodies (like the Earth and the Sun).
• The Takeaway: This specific type of friction makes the universe behave in a surprisingly predictable, "homographic" way, where the whole system shrinks or expands uniformly while spinning.

3. The Two-Body Dance: The "Spiral In"

The authors then zoomed in on just two bodies (like a planet and a star) to see exactly how this friction changes their path. They used a mathematical trick called "Poincaré compactification," which is like taking a map of the entire infinite universe and folding it onto a single, finite disk so you can see the edges.

What they found:

• The Conservative Case (No Friction): A planet can orbit forever in an ellipse, or fly off into space forever. It's like a ball rolling on a perfect, endless hill.
• The Dissipative Case (With Friction):
  • The Trap: If a planet starts with the right amount of energy, it doesn't fly away. Instead, it gets "captured." It spirals inward, slowly losing energy, until it settles into a perfect, circular orbit.
  • The Escape: If a planet is moving too fast, it can still escape, but it has to fight harder against the friction to get out.
  • The "Ghost" Point: They found a strange, mathematical point at "infinity" (the edge of the map) that acts like a funnel. All the paths that are about to escape or get captured pass near this point first. It's like a cosmic air traffic control tower that directs all incoming and outgoing traffic.

4. The Big Picture: What Happens to the Solar System?

The authors ran simulations to see how this affects the shape of orbits over time.

• The Orbit Shrinks: The distance between the bodies (the semi-major axis) always gets smaller. The planet slowly drifts closer to the star.
• The Orbit Rounds Out: The orbit becomes more and more circular. The "wobbly" elliptical paths smooth out into perfect circles.
• The Spin Stays: Even though the planet is getting closer and the orbit is changing, the total spin of the system remains locked.

Real-World Application:
This model helps explain why some things in space are getting closer (like Mercury slowly spiraling into the Sun) while others might be drifting apart (like the Earth and Moon, though that involves more complex rotation effects not fully covered here). It suggests that over millions of years, friction turns chaotic, wobbly orbits into neat, circular, synchronized dances.

Summary

The paper introduces a smart friction model that drains energy but saves the spin. They discovered that with a specific type of friction, the chaotic dance of many stars simplifies into the elegant dance of two. They mapped out exactly how planets spiral into perfect circles, showing that while the universe is full of chaos, this specific type of friction acts as a gentle hand, guiding wandering bodies into stable, circular orbits.

Technical Summary

1. Problem Statement

The paper addresses the dynamical evolution of the $N$-body problem under the influence of dissipative forces, specifically those arising from tidal effects.

• Context: In conservative gravitational systems, orbits are either regular or chaotic, with bodies potentially escaping, colliding, or maintaining relative equilibria. In dissipative systems (common in nature due to tides, drag, etc.), energy is lost, often leading to spin-orbit coupling, circularization, or collisions.
• The Gap: Most existing dissipative models (e.g., Poynting-Robertson drag, atmospheric drag) introduce a damping term proportional to velocity ($-\alpha \dot{r}$). While realistic for certain scenarios, these forces typically remove both energy and angular momentum from the system.
• Objective: The authors propose a mathematical model for dissipation that removes energy but strictly preserves the total angular momentum of the system. This is crucial for modeling planetary-scale tidal interactions where angular momentum conservation is a fundamental constraint.

2. Methodology

The authors develop a theoretical framework combining analytical mechanics, dynamical systems theory, and averaging methods.

A. The Dissipative Force Model

The authors introduce a pairwise dissipative force $f_{diss}^{ij}$ between masses $m_i$ and $m_j$:
$$f_{diss}^{ij} = -\alpha G \frac{m_i m_j}{|r_{ij}|^d} (\dot{r}_{ij} \cdot \hat{r}_{ij}) \hat{r}_{ij}$$

• Key Features:
  • Central Force: The force acts along the line connecting the bodies ($\parallel \hat{r}_{ij}$), ensuring conservation of total angular momentum.
  • Velocity Dependence: It depends on the radial component of the relative velocity, making it dissipative.
  • Homogeneity: The force is homogeneous of degree $d$ in position and velocity. The parameter $d$ is treated as a free parameter, though $d=8$ corresponds to Mignard's tidal model.
  • Symmetry: The force breaks time-reversal symmetry (dissipation) but maintains rotational symmetry.

B. Analysis of Central Configurations (CCs)

The study focuses on Central Configurations, where the acceleration of each body is proportional to its position relative to the center of mass.

• Homographic Motion: The authors restrict the analysis to planar motions where the configuration evolves via scaling ($s(t)$) and rotation ($\theta(t)$).
• Special Case $d=3$: They demonstrate that when the radial dependence exponent is $d=3$, the dissipative force term becomes mathematically proportional to the gradient of the gravitational potential. This creates an additional symmetry, reducing the complex $N$-body dissipative equations to a set of homographic equations equivalent to the dissipative two-body (Kepler) problem.

C. Topological Analysis (Poincaré Compactification)

To understand the global behavior of the dissipative two-body system, the authors perform a Poincaré compactification of the phase space.

• They map the infinite phase space $(r, v)$ onto a compact disk (Poincaré disk) using coordinate transformations (charts $U_1, U_2, V_2$).
• This allows for the rigorous analysis of dynamics at infinity (scattering/escape orbits) and the identification of equilibrium points, including those introduced by the regularization process.

D. Averaging Method

The authors apply averaging theory over Keplerian motion (mean anomaly) to derive long-term secular evolution equations.

• They transform the system into Delaunay variables $(M, \Lambda, \varpi, e)$ and compute the averaged dissipative terms using Hansen coefficients.
• This yields differential equations for the evolution of the semi-major axis ($a$), eccentricity ($e$), and periapsis longitude ($\varpi$).

3. Key Contributions and Results

A. Equivalence to the Dissipative Kepler Problem

For the specific case of $d=3$, the authors prove that the homographic evolution of any planar central configuration in an $N$-body system is governed by the same equations as a dissipative two-body system. This significantly simplifies the study of complex $N$-body tidal evolution.

B. Phase Space Topology and Bifurcations

The Poincaré compactification reveals a rich topological structure dependent on the dissipation strength $\alpha$:

• Conservative Case ($\alpha=0$): The phase space contains a separatrix (parabolic orbit) dividing bounded (elliptic) and unbounded (hyperbolic) orbits.
• Weak Dissipation: A "weak-center-manifold" emerges from a degenerate equilibrium at infinity.
• Strong Dissipation: As $\alpha$ increases, a topological transition occurs. The separatrix reconnects with the center manifold.
  • Result: All incoming orbits (scattering trajectories) are eventually captured and spiral toward the circular orbit attractor. The region of "escaping" orbits shrinks and eventually disappears for high dissipation, meaning all bodies with sufficient time will circularize.

C. Secular Evolution (Averaged Equations)

The averaged equations provide explicit rates of change for orbital elements:

1. Angular Momentum ($C$): Remains strictly constant ($\dot{C} = 0$).
2. Eccentricity ($e$): Monotonically decreases ($\dot{e} < 0$). The system evolves toward a circular orbit ($e \to 0$).
3. Semi-major Axis ($a$): Decreases over time ($\dot{a} < 0$). The final semi-major axis is determined solely by the initial angular momentum: $a_{final} = C^2 / \gamma$.
4. Periapsis Precession: The averaged equation for the longitude of periapsis is $\dot{\varpi} = 0$.
   • Significance: The model predicts no secular precession of the periapsis due to this specific type of dissipation. This distinguishes it from other drag models that might induce precession.

D. Physical Implications

• Earth-Moon vs. Mercury-Sun: The model explains why the Earth-Moon distance is increasing (which requires angular momentum transfer from Earth's rotation, not just orbital dissipation) but fails to explain it via pure orbital dissipation (which always shrinks $a$). Conversely, it correctly models the shrinking distance of Mercury toward the Sun due to tidal dissipation.
• Transient Dynamics: The paper highlights that while final states (collisions or circularization) are inevitable, systems may spend billions of years in "transient" semi-equilibrium states, resembling stability.

4. Significance

• Theoretical Rigor: The paper provides a rare analytical model that decouples energy dissipation from angular momentum loss, a critical feature for accurate tidal modeling in celestial mechanics.
• Mathematical Insight: The discovery that $d=3$ yields a symmetry reducing $N$-body central configurations to the two-body problem is a significant mathematical simplification.
• Topological Clarity: The use of Poincaré compactification offers a complete global picture of the phase space, clarifying the conditions under which bodies are captured versus those that escape, and how dissipation strength alters these boundaries.
• Predictive Power: The result of zero periapsis precession offers a specific, testable prediction for systems dominated by this type of tidal dissipation, distinguishing it from models involving atmospheric drag or radiation pressure.

In summary, the paper establishes a robust mathematical framework for understanding how tidal dissipation drives celestial systems toward circular, relative equilibrium states while conserving angular momentum, revealing complex topological transitions in the process.