🔭 astrophysics

A Stochastic Framework for the Spherical Jeans Equation Motivated by Scalar-Tensor Gravity

Motivated by scalar-tensor gravity, this paper develops a stochastic framework that models spatial fluctuations in the gravitational coupling as an additive noise term, transforming the spherical Jeans equation into a linear Itô stochastic differential equation to derive analytical expressions for the mean, variance, and covariance of radial and line-of-sight velocity dispersions in various dark matter halo models.

Original authors: Velasquez-toribio, A. M

Published 2026-06-11

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Velasquez-toribio, A. M

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to predict how fast stars are moving inside a giant, invisible cloud of dark matter (a "halo") that holds a galaxy together. For decades, astronomers have used a standard rulebook called the Jeans Equation to do this. Think of this rulebook like a perfect, rigid blueprint. It assumes that gravity is a fixed, unchanging force—like a constant weight on a scale—and that the galaxy is perfectly calm and still.

However, this new paper suggests that the blueprint might be too rigid. It proposes a "stochastic" (random) version of the rulebook, inspired by theories where gravity isn't just a fixed number, but a flexible field that can wiggle and fluctuate.

Here is the breakdown of the paper's ideas using simple analogies:

1. The Problem: Gravity Might Be "Wobbly"

In standard physics, gravity is like a heavy, solid anchor. But in some advanced theories (called Scalar-Tensor Gravity), gravity is more like a rubber band or a fluid. It can change strength depending on its surroundings.

The author suggests that inside a galaxy, there are tiny, unresolved fluctuations in this "gravity field." Imagine trying to walk across a floor that is mostly flat but has tiny, invisible bumps and dips that you can't see individually. These bumps represent the random fluctuations in gravity.

2. The Solution: A "Random Walk" for Stars

Instead of calculating a single, perfect speed for the stars, the author treats the stars' movement as a random walk along the radius of the galaxy.

The Old Way: You calculate one specific speed for every distance from the center. It's like drawing a single, smooth line on a graph.
The New Way: You calculate a cloud of possibilities. Because gravity is "wobbly," the speed of the stars isn't just one number; it's a range of likely numbers.

The paper uses a mathematical tool called a Stochastic Differential Equation. Think of this as a recipe that says: "Start with the standard speed, then add a little bit of random 'kick' at every step as you move inward."

3. The "Probability Band"

The most important result is that this randomness doesn't change the average speed. If you took a million different universes with slightly different gravity wobbles and averaged them, you would get the exact same result as the old, rigid rulebook.

However, the new method gives you a probability band (or a "fuzziness") around that average.

Analogy: Imagine aiming an arrow at a target. The old method says, "The arrow will hit the bullseye." The new method says, "The arrow will likely hit the bullseye, but it might land anywhere within a 2-inch circle around it."
This "circle" is the probability band. It tells astronomers: "We can't predict the exact speed of the stars at this distance, but we know it will fall within this specific range."

4. Testing the Theory with Different Galaxy Shapes

The author tested this new "wobbly gravity" idea on three common shapes of dark matter halos (the invisible clouds holding galaxies):

NFW: The standard, most common shape found in computer simulations.
Hernquist: A shape that is similar to NFW but has a finite edge (it doesn't go on forever).
Einasto: A shape that curves smoothly, changing its steepness as you go outward.

For all three shapes, the math worked out. The author showed that even with these different shapes, the "wobbly gravity" creates a predictable band of uncertainty around the star speeds.

5. Looking from the Outside (The "Line-of-Sight")

Astronomers can't see the 3D movement of stars inside a galaxy; they can only see the speed of stars moving toward or away from us (like looking at a car driving toward you on a highway).

The paper shows how to take the "wobbly" 3D speeds and project them onto this 2D view. The result is that the uncertainty band (the "fuzziness") also appears in the data we actually observe. This means that when astronomers measure star speeds, they should expect to see a spread of values, not just a single sharp line, and this spread could be caused by the natural fluctuations of gravity itself.

6. The "Pullback" Concept

The paper uses a clever way to visualize how these random speeds settle down. Imagine a river flowing from a wide mouth (the edge of the galaxy) toward a narrow source (the center).

If you drop a leaf at the edge, it might start at different spots.
As it flows inward, the current (the gravity) pushes all the leaves toward the same path, regardless of where they started.
The author calls this a "pullback attractor." It means that even if you are unsure about the conditions at the very edge of the galaxy, the random fluctuations inside will eventually guide the stars to a specific, predictable pattern of uncertainty as you get closer to the center.

Summary

This paper doesn't claim that gravity is broken or that our current models are wrong. Instead, it offers a more realistic way to handle the "noise" in the universe. It says: "Gravity isn't a perfectly smooth, fixed force; it has tiny, random ripples. If we account for those ripples, we don't get a single answer for how fast stars move; we get a range of likely answers, which matches the messy reality of the universe better."

The author notes that this specific math works best for calm, stable galaxies. For the wild, crashing edges of galaxy clusters where things are falling in rapidly, a different, more complex version of this math would be needed in the future.

Technical Summary: A Stochastic Framework for the Spherical Jeans Equation Motivated by Scalar-Tensor Gravity

Problem Statement
The standard spherical Jeans equation is a fundamental tool for describing the equilibrium dynamics of collisionless, self-gravitating systems, such as dark matter halos and elliptical galaxies. In its conventional application, the equation assumes a fixed, universal gravitational constant ( $G$ ). However, in scalar–tensor theories of gravity, the effective gravitational strength is not a rigid parameter but is controlled by a scalar field. If this scalar field exhibits spatial fluctuations, environmental dependence, or unresolved small-scale variations, the effective gravitational coupling ( $G_{\text{eff}}$ ) entering the dynamical equations may also fluctuate. The paper addresses the need for a controlled framework to quantify how these fluctuations in the gravitational sector propagate into the velocity-dispersion profiles of halos, moving beyond deterministic models that ignore such stochastic effects.

Methodology
The authors develop a stochastic framework by modeling unresolved spatial fluctuations of the scalar sector as an effective stochastic contribution to the gravitational coupling. The key methodological steps include:

Stochastic Coupling Formulation: The effective gravitational coupling is decomposed into a smooth, coarse-grained mean component and a stochastic fluctuation:
$G_{\text{eff}}(r, \omega) = \bar{G}_{\text{eff}}(r) + \Gamma_G(r)\xi(r, \omega)$
Here, $\Gamma_G(r)$ controls the radial amplitude of the fluctuation, and $\xi(r, \omega)$ is a zero-mean stochastic process. In the idealized white-noise limit, $\xi(r)$ is treated as the derivative of a Wiener process ( $dW_r/dr$ ), satisfying $\langle \xi(r)\xi(r') \rangle = \delta(r-r')$ .
Stochastic Differential Equation (SDE): Substituting this fluctuating coupling into the stationary spherical Jeans equation transforms it into a linear Itô stochastic differential equation for the radial velocity dispersion squared, $y(r) = \sigma_r^2(r)$ . The radial coordinate $r$ acts as an evolution parameter, defining a nonautonomous radial random flow rather than a time-evolution problem.
Fokker–Planck Analysis: The authors derive the associated Fokker–Planck equation to describe the probability density $P(y, r)$ of the velocity dispersion. Because the noise is additive (the diffusion coefficient is independent of $y$ ), the deterministic Jeans solution is recovered as the mean profile, while the stochastic sector generates a probability band around it.
Analytical Solutions for Halo Models: The framework is specialized to three standard halo density profiles: Navarro–Frenk–White (NFW), Hernquist, and Einasto. For each, the authors derive integral expressions for the mean, variance, and covariance of the radial velocity dispersion. These expressions rely on an integrating factor determined by the tracer density and anisotropy parameter.
Projection to Observables: The radial covariance structure is propagated to the projected line-of-sight velocity dispersion ( $\sigma_{\text{los}}^2$ ) using standard projection kernels, allowing for the calculation of uncertainty bands in observable kinematic quantities.

Key Contributions and Results

Derivation of Moment Equations: The paper provides closed integral expressions for the mean ( $\mu$ $μ$ ), variance ( $\Sigma^2$ $Σ^{2}$ ), and covariance ( $\text{Cov}$ $Cov$ ) of the radial velocity dispersion.
- The mean profile coincides exactly with the deterministic Jeans solution using the mean gravitational coupling.
- The variance is determined by the integral of the squared stochastic amplitude weighted by the integrating factor, explicitly showing how fluctuation amplitudes broaden the velocity-dispersion profile.
- The covariance between radii $r$ and $s$ is non-zero over the overlap of their integration intervals, reflecting the correlation induced by the stochastic forcing.
Gaussian Distribution: Due to the linearity of the equation and the Gaussian nature of the noise, the solution $y(r)$ is Gaussian distributed at any fixed radius (given a deterministic boundary condition).
Application to Specific Profiles:
- NFW: Requires a finite outer boundary due to logarithmic mass growth; kernels involve elementary functions.
- Hernquist: Features finite total mass and simpler rational kernels, serving as a useful benchmark.
- Einasto: Involves the lower incomplete gamma function for the enclosed mass but remains semi-analytical.
Pullback Attractor Interpretation: The authors interpret the stochastic flow as a "pullback attractor." Since the system is nonautonomous, the limiting random profile is selected by the inward propagation from the outer boundary. Different initial boundary conditions converge toward the same random profile, indicating a loss of memory of the imposed boundary data.
Numerical Validation: Monte Carlo realizations of the radial random flow are shown to agree with the analytical moment equations, validating the Fokker–Planck formulation.

Significance and Claims
The paper claims to provide a "semi-analytical framework" for assessing how effective gravitational fluctuations affect halo kinematic observables within the stationary Jeans regime. Its primary significance lies in:

Quantifying Uncertainty: It offers a method to generate probability bands for radial and projected velocity dispersions, rather than single deterministic curves, quantifying how scalar-field fluctuations broaden theoretical predictions.
Theoretical Consistency: It demonstrates that in an additive-noise model, the stochastic effects do not displace the mean profile but generate a probability distribution around the standard deterministic solution.
Scope Limitation: The authors explicitly restrict their scope to the stationary spherical Jeans regime (quasi-virialized regions). They note that extensions to cluster outskirts, splashback regions, or turnaround regions—where mean radial streaming velocities are significant—require a generalized Jeans equation and are left for future work. The framework is presented as a controlled way to propagate scalar-sector fluctuations into kinematic observables without invoking full non-equilibrium dynamics.