Very long-term relaxation of harmonic 1D self-gravitating systems

This paper numerically demonstrates that while non-degenerate one-dimensional self-gravitating systems relax on a timescale linear in particle number $N$, harmonic systems with degenerate orbital frequencies exhibit a significantly slower quadratic relaxation scaling ($N^2$), with partially degenerate systems transitioning between these regimes depending on the fraction of degenerate orbits.

The Gist

Here is an explanation of the paper "Very long-term relaxation of harmonic 1D self-gravitating systems," translated into everyday language with some creative analogies.

The Big Picture: A Cosmic Dance Floor

Imagine a giant, one-dimensional dance floor (a straight line) where thousands of people are dancing. These people are "particles" (like stars or atoms), and they are all holding hands with invisible rubber bands that pull them toward the center of the crowd. This is a self-gravitating system.

In the universe, these systems usually go through two phases:

1. The Violent Start: When the system is first formed, it's chaotic. Everyone is bumping into each other, and the shape changes rapidly until it settles into a somewhat stable pattern.
2. The Slow Drift: Once settled, the system looks frozen, but it's actually slowly changing over millions of years. Tiny, random nudges from other dancers cause everyone to slowly drift toward a perfect, calm equilibrium. This slow drift is called relaxation.

The Problem: The "Perfectly Tuned" Orchestra

For decades, physicists have had a rulebook (called the Landau and Balescu–Lenard theories) to predict how fast this "slow drift" happens. The rulebook says:

"The time it takes to relax is directly proportional to the number of dancers ($N$). If you double the crowd, it takes twice as long to settle down."

However, this rulebook has a loophole. It assumes that every dancer moves at a slightly different speed. If everyone moves at the exact same speed, the math breaks down.

The authors of this paper decided to test this loophole. They created a simulation of a Harmonic System—a special kind of dance floor where the "rubber bands" are tuned perfectly so that every single particle oscillates back and forth at the exact same frequency. It's like an orchestra where every instrument plays the exact same note at the exact same time.

The Discovery: The "Traffic Jam" Effect

The researchers ran massive computer simulations to see what happens in this perfectly tuned, "degenerate" system. Here is what they found:

1. The Standard Crowd (Non-Degenerate):
   If the dancers have different speeds (like a normal crowd), the relaxation time scales linearly ($N$). If you have 100 people, it takes $X$ time. If you have 1,000, it takes $10X$ time. This matches the old rulebook.

2. The Perfectly Tuned Crowd (Harmonic/Degenerate):
   When everyone moves at the exact same speed, the system gets stuck in a massive traffic jam.

   • The Result: The relaxation time doesn't just double; it quadruples when you double the crowd.
   • The Math: Instead of scaling with $N$, the time scales with $N^2$ (N squared).
   • The Analogy: Imagine a hallway where everyone walks at the exact same speed. If someone tries to change lanes or speed up, they can't because everyone else is perfectly synchronized with them. They are "locked" in place. To get out of this lock, they have to wait for a very rare, specific coincidence of movements. As the crowd gets bigger, these rare coincidences become exponentially harder to find.

The "Halfway" Crowd (Partially Degenerate)

The authors also tested a middle ground: a crowd where some people are perfectly synchronized, but others are moving at random speeds.

• Small Crowds: Even with a few random movers, the synchronized group dominates, and the system behaves like the "traffic jam" ($N^2$).
• Large Crowds: Once the crowd gets big enough, the random movers become numerous enough to break the lock. The system suddenly switches back to the normal behavior ($N$).
• The Lesson: The more synchronized (degenerate) the crowd is, the larger the group needs to be before it can "break free" and start relaxing normally.

Why Does This Matter? (The Real World)

You might ask, "Who cares about 1D lines of particles?" The authors explain that this is actually a key to understanding real galaxies:

• Dwarf Galaxies: Many small galaxies have "cores" in their centers that act like these harmonic systems. They are so dense and symmetric that the stars inside them move in a synchronized way.
• The "Stalling" Problem: Astronomers have been puzzled by why globular clusters (balls of stars) sometimes stop falling toward the center of a dwarf galaxy. They should crash into the center due to gravity, but they get stuck.
• The New Clue: This paper suggests that these cores are "thermodynamically blocked." Because the stars are so perfectly synchronized, they can't easily exchange energy or move around. They are stuck in a slow-motion traffic jam, which explains why they don't sink to the center as fast as we thought they would.

Summary in a Nutshell

• Normal Systems: Relaxation time grows linearly with size ($N$).
• Harmonic Systems (Perfectly Synchronized): Relaxation time grows quadratically with size ($N^2$). They are incredibly slow to change.
• The Takeaway: When a system is too perfectly ordered, it becomes "stuck." It takes a massive amount of time (and a huge number of particles) for the chaos of the universe to finally break the symmetry and let the system settle down.

The authors used a super-precise computer program (an "exact integrator") to prove this, showing that for these special systems, the old rulebooks were wrong, and the universe is even more patient (and stubborn) than we thought.

Technical Summary

Here is a detailed technical summary of the paper "Very long-term relaxation of harmonic 1D self-gravitating systems" by Tep, Fouvry, and Pichon.

1. Problem Statement

The paper addresses a fundamental gap in the kinetic theory of self-gravitating systems. Standard theories, specifically the Landau and Balescu–Lenard (BL) equations, describe the long-term relaxation of systems driven by two-body resonant interactions (Poisson shot noise). These theories predict that the relaxation time ($t_{rel}$) scales linearly with the number of particles ($N$), i.e., $t_{rel} \propto N t_{dyn}$.

However, these theories rely on the assumption that the system is non-degenerate, meaning the mapping from orbital actions to orbital frequencies has a non-zero Jacobian everywhere. In harmonic potentials, all particles share the exact same mean orbital frequency (dynamical degeneracy). In this regime:

• The resonance condition in the BL equation becomes ill-defined (the Dirac delta function $\delta(k\Omega - k'\Omega')$ becomes singular for all $k=k'$).
• The standard quasilinear kinetic theory breaks down, leaving the scaling of the relaxation time for harmonic systems theoretically unconstrained.

The authors aim to numerically determine the relaxation timescale of fully degenerate (harmonic) and partially degenerate 1D self-gravitating systems to understand how dynamical degeneracy alters long-term evolution.

2. Methodology

The authors employ a highly precise numerical approach to overcome the limitations of standard kinetic theory:

• System Model: They study 1D self-gravitating systems on an infinite line with equal-mass particles. The gravitational interaction kernel is $U(x, x') = G|x - x'|$.
• Integrator: They utilize an exact collision-driven N-body integrator.
  • In 1D, the force is constant between particle crossings (collisions).
  • The integrator solves the equations of motion exactly between collisions and handles particle crossings by swapping indices, velocities, and masses.
  • To manage the accumulation of round-off errors over the vast number of collisions required for long-term relaxation, they use 80-bit double-float precision (via the Julia DoubleFloats library), achieving relative energy/momentum conservation errors of $\sim 10^{-23}$.
• Potentials Investigated:
  • Plummer: Non-degenerate, infinite radial support.
  • Compact: Non-degenerate, compact radial support.
  • Harmonic: Fully degenerate (constant frequency), compact support.
  • Anharmonic: Partially degenerate, compact support. The fraction of non-degenerate orbits is controlled by a parameter $\epsilon$ ($0 \le \epsilon \le 1$).
• Relaxation Metric:
  • They define relaxation by measuring the distance between the instantaneous cumulative mass function $M(\le x, t)$ and the theoretical thermodynamic equilibrium $M_{th}(\le x)$.
  • The metric used is the Kolmogorov–Smirnov (KS) distance, $D_{KS}(t)$.
  • The relaxation time $t_{rel}$ is defined as the time when $D_{KS}(t)$ drops below a statistically significant threshold ($D_0 = 0.015$).
• Simulation Parameters: Simulations were run for particle numbers $21 \le N \le 141$, with ensemble averaging over$N_r$realizations such that$N \times N_r \approx 10^5$.

3. Key Contributions and Results

A. Relaxation Scaling in Fully Degenerate Systems (Harmonic)

The most significant finding is that harmonic systems relax on a timescale scaling quadratically with $N$:
$$t_{rel} \propto N^2 t_{dyn}$$
This contrasts sharply with the linear scaling ($t_{rel} \propto N t_{dyn}$) predicted by BL theory for non-degenerate systems. The authors note that this extremely slow relaxation cannot be explained by standard kinetic blocking (which usually yields $1/N^2$ fluxes in non-degenerate systems) because the harmonic system is "over-resonant" rather than "under-resonant." The phase mixing assumption, crucial for quasilinear expansions, fails completely here.

B. Behavior of Non-Degenerate Systems

• Plummer and Compact: Both systems exhibit the expected linear scaling $t_{rel} \propto N t_{dyn}$.
• Prefactor Difference: While the scaling is linear for both, the Compact system relaxes approximately 10 times slower than the Plummer system. This is attributed to "quasi-kinetic blocking": the Compact system has a narrower range of orbital frequencies, requiring higher-order resonances to satisfy the BL resonance condition, which are less efficient.

C. Transition in Partially Degenerate Systems (Anharmonic)

The study of Anharmonic systems (varying $\epsilon$) reveals a dynamical regime transition:

• Low $N$: For small particle numbers, even partially degenerate systems exhibit the quadratic scaling ($t_{rel} \propto N^2$) characteristic of the fully degenerate harmonic case.
• High $N$: As $N$ increases, the system transitions to the linear scaling ($t_{rel} \propto N$) typical of non-degenerate systems.
• Dependence on Degeneracy: The value of $N$ at which this transition occurs increases as the fraction of degenerate orbits increases (i.e., as $\epsilon \to 0$). This confirms that dynamical degeneracy delays the onset of standard kinetic relaxation.

4. Significance and Astrophysical Implications

• Theoretical Impact: The paper provides the first numerical evidence that dynamical degeneracy fundamentally alters the scaling of relaxation times in long-range interacting systems, invalidating the standard $O(N)$ prediction for harmonic cores. It suggests that the mechanism driving this slow relaxation is likely thermodynamic blocking (as proposed by Deme & Fouvry 2025), where adiabatic perturbations prevent heat and matter flows.
• Astrophysical Context:
  • Dwarf Galaxy Cores: The results offer new insights into the "core stalling" problem, where globular clusters fail to sink to the center of dwarf galaxies. The slow relaxation ($N^2$) in harmonic cores implies that dynamical friction is significantly less efficient than standard theory predicts, potentially explaining the persistence of off-center active galactic nuclei or stalled clusters.
  • Vertical Diffusion: The findings may refine models of vertical diffusion in galactic discs, where harmonic approximations are often used.
• Future Directions: The authors propose extending this work to 3D spherical clusters, investigating multi-mass systems, and developing renormalization group techniques or time-averaged kinetic theories to analytically describe the $N^2$ scaling observed.

Conclusion

This work demonstrates that the standard kinetic theory of self-gravitating systems fails for dynamically degenerate (harmonic) potentials. Through high-precision 1D simulations, the authors establish that harmonic systems relax on a timescale of $O(N^2 t_{dyn})$, a result driven by the breakdown of phase mixing and the ill-posed nature of the Balescu–Lenard equation in degenerate regimes. This has profound implications for understanding the long-term evolution of dense stellar systems and galactic cores.