Linearization Principle: The Geometric Origin of Nonlinear Fokker-Planck Equations

This paper presents a geometric derivation of the nonlinear Fokker-Planck equation by utilizing the $q$-logarithm as a natural coordinate system, thereby establishing a thermodynamic framework that links the dynamic index $q$ to a dual thermodynamic index $2-q$ to explain anomalous diffusion and power-law distributions without ad-hoc constraints.

The Gist

Imagine you are trying to understand how a crowd of people moves through a city. In a normal, calm city (standard physics), people walk in straight lines, bump into each other randomly, and spread out evenly over time. This is "normal diffusion."

But sometimes, in complex systems like stock markets, turbulent rivers, or even how viruses spread, things don't behave normally. People might suddenly sprint in one direction, get stuck in a corner, or form strange, heavy-tailed clusters. This is "anomalous diffusion."

For a long time, physicists have struggled to explain why this happens. They tried to force the math to fit the weird behavior by inventing complicated, "made-up" rules (like saying the wind changes direction based on how crowded the street is). These rules felt artificial and broke the fundamental laws of how particles usually interact with their environment.

Hiroki Suyari's paper proposes a much simpler, more elegant solution. He suggests that the weird behavior isn't because the rules of physics are broken, but because we are looking at the city through the wrong map.

Here is the breakdown of his idea using simple analogies:

1. The Wrong Map vs. The Right Map (The Growth Law)

Imagine you are trying to measure the growth of a plant.

• The Old Way: You measure the plant's height in inches every day. If the plant grows exponentially (doubling every day), your graph looks like a curve that shoots up to the sky. It's hard to predict.
• Suyari's Way: He suggests that for certain complex systems, the "natural" way to measure things isn't a straight line (inches), but a curved ruler.

He starts with a simple rule: "The rate of growth depends on the current size raised to a power ($q$)."

• If $q=1$, it's normal growth (a straight line on a log scale).
• If $q \neq 1$, it's "power-law" growth (the kind seen in complex systems).

Suyari realized that if you switch your coordinate system to a special "curved ruler" called the $q$-logarithm, that messy, curved growth suddenly looks like a straight line. It's like realizing that the Earth looks flat when you zoom in, but you need a globe (a curved map) to understand the whole world.

2. The "Linearization Principle"

This is the core of the paper. Suyari says: "Don't change the laws of physics; change the map."

In standard physics, the force pushing a particle (drift) is usually linear. In previous attempts to explain anomalous diffusion, scientists had to invent "nonlinear forces" (forces that change depending on how crowded the area is). This felt like cheating.

Suyari argues:

• If you use the $q$-logarithm as your natural map, the "drift" force remains linear and simple (just like in normal physics).
• The "nonlinearity" (the weird behavior) only appears because the map itself is curved.
• Analogy: Imagine driving a car on a curved road. If you look at a flat map, your path looks like a crazy, winding snake. But if you look at a map that is curved to match the road, your path is just a straight line. The car isn't doing anything weird; the map was just wrong.

3. The Magic Duality (The $q$ and $2-q$ Switch)

Here is the most fascinating part of the discovery. Suyari finds a beautiful symmetry, or duality, in nature:

• The Dynamic Index ($q$): This number describes how the particles move moment-to-moment (the local rules).
• The Thermodynamic Index ($2-q$): This number describes the overall stability and the final shape of the crowd (the global rules).

The Analogy:
Think of a dance party.

• $q$ is the rhythm the dancers are following right now. If the rhythm is fast and chaotic ($q > 1$), they move wildly.
• $2-q$ is the "energy cost" of the party. Even if the dancers are moving wildly, the overall energy of the room settles down in a specific pattern that follows a different rule.

Suyari proves that you don't need to invent "ghost particles" (called escort distributions in old theories) to make the math work. The system naturally balances itself using this $q$ vs. $2-q$ switch.

4. What Does This Actually Predict?

When you apply this new "curved map" to real-world problems, it predicts two famous things perfectly:

1. The Harmonic Oscillator (The Spring): If you trap a particle in a spring, it settles into a $q$-Gaussian distribution.

   • If $q=1$, it's a normal bell curve (Gaussian).
   • If $q > 1$, it has "fat tails" (extreme events are more likely).
   • If $q < 1$, it has "sharp edges" (the particle can't go beyond a certain point).
   • Why this matters: This explains why financial crashes or extreme weather events happen more often than standard bell curves predict.

2. The Free Particle (Anomalous Diffusion): If a particle is just wandering in open space:

   • It doesn't spread out at a normal speed.
   • It spreads at a speed determined by $q$.
   • This perfectly matches observations in turbulent fluids and biological systems.

The Big Takeaway

For decades, physicists tried to fix the "engine" of the car (the equations of motion) to explain why the car was driving weirdly. They added complex, ad-hoc parts that didn't make sense physically.

Suyari says: "The engine is fine. You just need to look at the road through a different lens."

By recognizing that the natural coordinate system for these complex systems is curved (the $q$-logarithm), he restores the simplicity of physics. The drift remains linear, the Einstein relation (the link between heat and movement) stays standard, and the weird, power-law behavior emerges naturally from the geometry of the space itself.

It's a reminder that sometimes, to understand the complexity of the universe, we don't need more complicated math—we just need to change our perspective.

Technical Summary

1. Problem Statement

The paper addresses a fundamental dilemma in statistical physics regarding the dynamical origins of anomalous diffusion and power-law distributions (specifically $q$-Gaussian distributions).

• The Conflict: While generalized entropies (like Tsallis entropy) successfully describe the static properties of these systems, deriving the corresponding dynamical equations (Nonlinear Fokker-Planck Equations, NFPEs) has been problematic.
• Existing Limitations: Previous phenomenological approaches to derive NFPEs often require:
  • State-dependent diffusion coefficients.
  • Nonlinear drift terms (forces that depend on the probability density $p$ itself), which violate the physical principle that particles should respond linearly to external potentials.
  • The use of escort distributions (auxiliary measures) to maintain mathematical consistency, which obscures the physical meaning of observables.

2. Methodology: The Linearization Principle

The author proposes a geometric derivation of the NFPE based on the Linearization Principle, rather than postulating a specific entropy functional.

• Fundamental Growth Law: The derivation starts with a nonlinear growth law for a state variable $y(x)$:
  $$\frac{dy}{dx} = y^q$$
  where $q$ is an interaction index ($q=1$ yields standard exponential growth; $q \neq 1$ yields power-law behavior).
• Natural Coordinates: The equation is linearized via the transformation $z = \ln_q y$, where $\ln_q$ is the $q$-logarithm:
  $$\ln_q(y) = \frac{y^{1-q} - 1}{1-q}$$
  The principle dictates that physical laws (diffusion and drift) must take linear forms within this natural coordinate system ($z$).
• Generalized Flux: Applying linear response theory in the $q$-logarithmic coordinate system:
  • The chemical potential is defined as $\mu_q = k_B T \ln_q p + V(x)$.
  • The flux $J$ is derived as $J = -\frac{p}{\gamma} \nabla \mu_q$.
  • Using the identity $\nabla(\ln_q p) = p^{-q} \nabla p$, the flux becomes:
    $$J = -D p^{1-q} \nabla p - \frac{1}{\gamma} p \nabla V$$
• Resulting Equation: Substituting the flux into the continuity equation yields the Nonlinear Fokker-Planck Equation (NFPE):
  $$\frac{\partial p}{\partial t} = \nabla \cdot \left( D p^{1-q} \nabla p + \frac{1}{\gamma} p \nabla V \right)$$
  Crucially, the drift term ($\nabla \cdot (p \nabla V)$) remains linear in $p$, preserving the standard response to external potentials.

3. Key Contributions

• Geometric Origin of Nonlinearity: The paper demonstrates that nonlinearity in the diffusion term arises naturally from the geometry of the state space (dictated by the growth law), not from ad-hoc phenomenological modifications.
• Preservation of the Einstein Relation: Unlike previous models, this framework preserves the standard Einstein relation ($D = k_B T / \gamma$) and the linearity of the drift term, avoiding fictitious nonlinear forces.
• Elimination of Escort Distributions: The theory provides a consistent thermodynamic framework without needing auxiliary escort distributions.
• Duality of Indices: The paper establishes a fundamental duality between the dynamic index $q$ (governing local dynamics) and the thermodynamic index $2-q$ (governing global stability).

4. Key Results

A. Stationary States and $q$-Gaussians

The stationary solution ($J=0$) of the derived NFPE yields the $q$-Gaussian distribution:
$$p_{st}(x) = \frac{1}{Z_q} \exp_q \left( -\beta_q V(x) \right)$$
where $\exp_q$ is the $q$-exponential function. For a harmonic potential $V(x) \propto x^2$, this recovers the power-law tails observed in complex systems.

B. The $q \leftrightarrow 2-q$ Duality and H-Theorem

• Free Energy Functional: The system minimizes a free energy functional $F[p]$ defined by a generalized entropy of index $2-q$:
  $$S_{2-q}[p] = \frac{1 - \int p^{2-q} dx}{(2-q) - 1}$$
• H-Theorem: The time derivative of the free energy is non-positive ($dF/dt \leq 0$), proving that the system monotonically approaches the stationary state.
• Significance: While the local dynamics are governed by $q$-arithmetic, the global thermodynamic stability is governed by the dual index $2-q$. This resolves the ambiguity in choosing entropic indices in non-extensive statistical mechanics.

C. Applications

1. Harmonic Oscillator: The theory predicts $q$-Gaussian stationary states.
   • $q > 1$: Heavy-tailed distributions (super-diffusive).
   • $q < 1$: Compact support distributions (sub-diffusive).
2. Free Particle ($V=0$): The equation reduces to the Porous Medium Equation (PME).
   • The mean square displacement follows an anomalous diffusion law: $\langle x^2(t) \rangle \propto t^{\frac{2}{3-q}}$.
   • This links the parameter $q$ directly to the macroscopic diffusion exponent.

5. Significance and Implications

• Theoretical Unification: The work bridges the gap between nonlinear dynamics and statistical mechanics by showing that the "nonlinear" behavior is a projection of linear dynamics onto a curved statistical manifold (Information Geometry).
• Physical Consistency: By maintaining linear drift terms, the model ensures that particles respond physically to external potentials, a feature often lost in previous phenomenological NFPEs.
• Quantum Connection: The author notes that this geometric framework serves as the foundation for understanding macroscopic states of 1D quantum fluids (referenced in companion paper [7]), linking classical anomalous diffusion to strongly correlated quantum systems.
• Future Directions: The framework provides a consistent basis for extending these analyses to multi-dimensional systems and complex non-equilibrium phenomena without relying on ad-hoc constraints.

In summary, Suyari's paper re-grounds the theory of anomalous diffusion in a rigorous geometric principle, demonstrating that the complex statistics of power-law systems emerge naturally from the geometry of the state space, thereby resolving long-standing inconsistencies in the dynamical formulation of non-extensive statistical mechanics.