Brownian motion of supermassive black holes in galaxy cores

Carlo Nipoti; Luca Ciotti; Pierfrancesco Di Cintio

REVIEW 3 major objections 3 minor 1 cited by

Brownian motion of supermassive black holes in galaxy cores

T0 review · 3 major / 3 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The off-center position of M87's supermassive black hole can be explained by stochastic gravitational encounters with stars, with the black hole wandering roughly 6 parsecs from the galaxy center in 10 billion years.

desk verdict A plausible mechanism for SMBH wandering, but the headline '≈6 pc' depends on an unstated cutoff in the noise distribution. read the letter →

arxiv 1908.04283 v1 pith:BV4ALNDK submitted 2019-08-12 astro-ph.GA

Pierfrancesco Di Cintio , Luca Ciotti , Carlo Nipoti This is my paper · ORCID

classification astro-ph.GA

keywords stellardynamicsblackholephysicsn-bodysimulationsmethods:statisticalBrownianmotionLangevinequationHoltsmarkdistributionM87

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper argues that the observed off-center position of the supermassive black hole in M87 can be explained by purely dynamical 'Brownian' motion: random gravitational kicks from passing stars slowly push the black hole away from the galaxy's center. The authors build a semi-analytic Langevin model that combines dynamical friction with stochastic noise drawn from the Holtsmark distribution, and they validate it against direct N-body simulations of intermediate-mass black holes. For an M87-like galaxy they find the black hole wanders about 6 parsecs from the center over 10 billion years, matching the offset claimed for M87. If correct, the result suggests that some off-center supermassive black holes require no exotic formation mechanism, only ordinary stellar encounters.

What carries the argument

The load-bearing object is the Langevin equation for the black hole's acceleration: $\ddot{\mathbf{r}} = -\nabla \Phi_{\mathrm{tot}}(\mathbf{r}) - \eta(\mathbf{r},v)\,\mathbf{v} + \mathbf{F}(\mathbf{r})$, where $\Phi_{\mathrm{tot}}$ is the smooth stellar-plus-dark potential, $\eta$ is the dynamical-friction coefficient, and $\mathbf{F}$ is a stochastic acceleration per unit mass. The magnitude of $\mathbf{F}$ is sampled from the Holtsmark distribution, the classic fluctuating-force distribution for a homogeneous Poisson field of point masses, evaluated at the local stellar density. Because the Holtsmark distribution has infinite variance, the scheme imposes a cutoff on large forces to make it renormalizable, and the equations are integrated with a quasi-symplectic scheme. This combination lets the model reach mass ratios $M_{\mathrm{BH}}/m_* \sim 10^9$ that direct N-body simulations cannot, while still matching the radial displacement seen in N-body runs at mass ratio 100.

What would settle it

Run the same Langevin integration with the large-force cutoff on the Holtsmark distribution increased and decreased by a factor of 10; if the 10-Gyr displacement changes by more than a factor of a few, the predicted ~6 pc wander is not a robust consequence of the model.

Watch

Extended reading notes

Core claim

The central claim is that a supermassive black hole in a galactic core undergoes Brownian-like motion driven by gravitational encounters with stars, and that for an M87-like system this motion produces a displacement of roughly 6 pc after 10 Gyr. The paper demonstrates this by solving a Langevin equation with a Holtsmark-distributed stochastic force, and shows that the Holtsmark noise matches direct N-body simulations better than Gaussian noise for a Plummer cluster with mass ratio 100. For M87 parameters (MBH = 6e9 solar masses, total mass 3e12 solar masses, γ = 1.2, rc = 3 kpc), the model yields radii of order 6 pc, consistent with the off-centre displacement claimed for M87. The authors emphasize that this offset arises only from multiple dynamical collisions with stars.

Load-bearing premise

The prediction rests on assuming that the stochastic force on the black hole follows a Holtsmark distribution with a large-force cutoff, evaluated at the local stellar density, and that this description still holds when extrapolated from the tested mass ratio of 100 to the galactic ratio of about $10^9$.

Editorial extensions

If this is right

For an M87-like galaxy, the model predicts a black-hole displacement of about 6 pc after 10 Gyr of stellar encounters alone.
A stochastic-noise description of stellar encounters reproduces the radial displacement seen in direct N-body simulations better when the noise is drawn from a Holtsmark distribution than from a Gaussian.
Off-center supermassive black holes in massive elliptical galaxies do not require binary-black-hole recoil or asymmetric accretion to explain their displacements; ordinary two-body relaxation can do it.
The Langevin approach with Holtsmark noise allows exploration of supermassive-black-hole dynamics at mass ratios inaccessible to direct N-body simulations, up to about $10^9$.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

If the few-parsec wander is generic for massive ellipticals, the model implies a population of off-center supermassive black holes whose offset distribution could be compared with future high-angular-resolution samples; the paper does not compute this distribution.
The cutoff on the Holtsmark force is a free parameter; systematically varying it, calibrated against higher-resolution N-body runs, could turn the present order-of-magnitude estimate into a quantitative prediction.
The same Langevin machinery could be applied to intermediate-mass black holes in globular clusters, where direct N-body validation is possible, to predict observable offsets or ejection rates.
For cuspy galaxies with steeper central density slopes, the local stellar density near the center is higher, so one might expect larger stochastic kicks and faster wander; the paper only explores the case γ = 1.2.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Colleague,

This is a short proceedings paper, but it makes a concrete claim: a 6×10^9 solar-mass black hole in an M87-like galaxy wanders about 6 pc in 10 Gyr because of stellar encounters. The new bit is not Brownian motion itself — that is established — but applying the Langevin equation with Holtsmark noise to galactic SMBHs and getting a few-parsec displacement. The N-body comparison for an intermediate-mass black hole (mass ratio 100) is a genuine check, and the Holtsmark noise does reproduce the N-body wandering better than Gaussian noise, at least qualitatively.

The soft spot is real and it is the one the stress-test note flags. The Holtsmark distribution has infinite variance, the Mannella scheme needs a finite second moment, so the authors truncate at a large-F cutoff. They give no value, no algorithm, no physical prescription for that cutoff. The diffusion coefficient that sets the displacement scales like F_max^{1/2} with the tail they use, so the '≈6 pc' in Fig. 2 is not a unique prediction of the model as written. This is not a minor implementation detail; it is the parameter that controls the answer. The N-body calibration in Fig. 1 does not fix it, because the same cutoff is used in the stochastic runs and the match is visual. The extrapolation from mass ratio 100 to 10^9 is also uncontrolled. The M87 comparison is a consistency check against a contested offset, not a derivation.

None of this kills the underlying mechanism. The Langevin-plus-Holtsmark approach is reasonable, and a fuller paper with a stated cutoff (e.g. set by the finite size of stars or by the maximum single-star acceleration) and quantitative error bars could make the prediction meaningful. As it stands, the paper is honest about being preliminary and the limitation is stated, but the central number is not yet pinned down.

I would not cite the 6 pc value, but I might cite the method if it gets completed. Bring it to reading group only if you want to discuss how stochastic descriptions of encounters need regularization. I would send it to a referee if it were submitted as a full paper, because the mechanism deserves scrutiny and the N-body comparison is worth checking; but the referee should insist on a specified cutoff.

Referee Report

3 major / 3 minor

Summary. This paper develops a Langevin-equation model for the Brownian motion of a supermassive black hole (SMBH) in a galactic core, combining smooth stellar and dark matter potentials, Chandrasekhar dynamical friction, and a stochastic force drawn from the Holtsmark distribution for gravitational fluctuations. The model is tested against direct N-body simulations of a 100-to-1 mass-ratio black hole in a Plummer cluster, and then applied to an M87-like galaxy, yielding a predicted displacement of about 6 pc over 10 Gyr. The authors argue that this displacement is compatible with reported off-centering of M87's SMBH and that the stochastic approach is computationally advantageous.

Significance. If the predicted few-parsec wander were established, the paper would provide a physically transparent, purely stellar-dynamical mechanism to explain off-center SMBHs, with direct relevance to observations of M87 and other nearby galaxies. The modeling strategy is appealing, and the comparison with direct N-body simulations is a sensible way to test the stochastic approach. However, the central quantitative prediction currently depends on an unspecified truncation of the Holtsmark distribution, so the significance of the result as stated is not yet established. At this stage the paper is better viewed as a promising method and a preliminary consistency check rather than as a completed derivation.

major comments (3)

[Section 2, after Eq. (2.8)] The Holtsmark distribution has divergent standard deviation, and the manuscript states that a 'cut-off large F' is imposed, but no value or physical prescription for this cutoff is given. Because the noise amplitude in the Mannella scheme is set by the second moment of the truncated distribution, and because this second moment grows as F_max^{1/2} when the asymptotic form in Eq. (2.7) is used, the diffusion coefficient and hence the r(t) shown in Fig. 2 depend on the arbitrary F_max. The headline '≈ 6 pc' is therefore not a unique prediction of the model as written. The authors should specify how F_max is chosen (for example from a minimum impact parameter or local stellar density) and show the sensitivity of Fig. 2 to this choice.
[Section 3, Figure 1] The validation against N-body simulations is qualitative: the bottom panels show that the N-body distributions lie between the Gaussian and Holtsmark cases, and the text itself says the N-body behavior is 'somewhat intermediate.' No quantitative error metric or parameter inference is provided. More importantly, the N-body comparison cannot calibrate the cutoff because the same unspecified cutoff enters the stochastic runs that are being compared. A quantitative match, for example a goodness-of-fit statistic as a function of F_max, would be needed to claim that the stochastic model is validated.
[Section 3, application to M87] The extrapolation from the validated case MBH/m* = 100 (with N = 10^4 or 3 x 10^4) to a galactic SMBH with MBH/m* ~ 10^9 is not justified by any demonstrated scale invariance. The Holtsmark force is evaluated with the local stellar density, but the gravitational influence of the SMBH changes the stellar distribution and encounter statistics in its vicinity, and this effect may not be captured by the simple local-density prescription. The authors should provide a dimensional-analysis or numerical argument that the stochastic description remains valid over this seven-order-of-magnitude range in mass ratio.

minor comments (3)

[Figure 1 caption and Section 3] The caption of Figure 1 reports N = 10^4 particles for the N-body run, while the text of Section 3 says '3 x 10^4 stars'; these numbers should be reconciled.
[Section 3] In the sentence 'the case using the Holtsmak distribution better approaches the results', 'Holtsmak' is a typo for 'Holtsmark'.
[Section 3] The physical meaning of the adopted parameters should be stated explicitly: it is unclear whether Mgal = 3 x 10^12 M_sun is the total mass of the galaxy model or only the stellar mass, and how rc = 3 kpc is shared between the stellar and dark components.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the M87 few-parsec displacement is an integration output, not a fitted input.

full rationale

The derivation chain is self-contained in the following sense. Section 2 specifies a Langevin equation whose inputs are the dynamical friction coefficient (Eq. 2.3), the stellar and dark matter density profiles (Eqs. 2.1 and 2.2), and a stochastic force drawn from the Holtsmark distribution (Eq. 2.6). None of these inputs is defined in terms of the target M87 displacement. Section 3 validates the stochastic model against a direct N-body simulation at MBH/m*=100, which is an external benchmark rather than a calibration to the M87 offset. The galactic run uses M87-like parameters (Mgal=3e12 Msun, gamma=1.2, rc=3 kpc, dark-to-visible ratio ~6) taken from independent references (Wu & Tremaine 2006; Event Horizon Telescope Collaboration 2019), and the quoted r(10 Gyr) ~ 6 pc is presented as an output of the integration. The paper's own caveat that the Holtsmark distribution has infinite variance and must be truncated at a large-F cutoff (Section 2, immediately after Eq. 2.8) is a genuine robustness and underdetermination concern: the diffusion coefficient, and hence the wander amplitude, can depend on the unspecified cutoff. However, the text does not state that the cutoff was chosen to reproduce 6 pc or any observed offset, so this is a correctness risk rather than a circular reduction. The only self-citation, Alessandrini et al. (2014) for the spatial dependence of the dynamical friction coefficient, is a standard approximation and is not load-bearing for the headline result. No equation reduces the prediction to its inputs by construction, and no fitted parameter is renamed as a prediction.

Assumptions & free parameters 2 free parameters · 5 assumptions · 0 invented entities

The model's free parameters are the Holtsmark cutoff and Coulomb logarithm (both unspecified), which control the noise amplitude and friction strength. It relies on standard Chandrasekhar dynamical friction with a Maxwellian velocity distribution, a local homogeneous-Poisson noise model, a Markovian approximation, and the test-particle approximation for the SMBH. No invented entities are introduced.

free parameters (2)

Holtsmark noise cutoff F_max
The Holtsmark distribution has infinite variance; a large-F cutoff is introduced in Section 2 to make it normalizable for the Mannella integrator. Its value is not stated, and it directly controls the amplitude of the stochastic kicks and hence the resulting displacement.
Coulomb logarithm ln Lambda
Appears in the dynamical friction coefficient (Eq. 2.3); the value used is not specified. Typical choices are 10-15, but the displacement could depend on it.

assumptions (5)

domain assumption Dynamical friction can be described by the Chandrasekhar formula with a local Maxwellian velocity distribution (Eqs. 2.3, 2.5).
Assumes an isotropic, Maxwellian stellar velocity distribution with dispersion from Jeans equations; this neglects anisotropy and the SMBH's influence on the local distribution.
domain assumption The stochastic force follows the Holtsmark distribution for a homogeneous Poisson field, evaluated at the local stellar density (Eq. 2.6).
The galaxy has a strong density gradient, whereas the Holtsmark distribution assumes a uniform random distribution of field stars; using the local density is an approximation.
domain assumption The noise and friction are Markovian and delta-correlated in time (Eq. 2.8 and the delta-function correlation stated in Section 2).
Gravitational encounters have finite duration; the Markovian approximation is standard but may fail for high-mass ratios.
domain assumption The SMBH is a test particle; its gravitational potential does not modify the stellar distribution.
The SMBH's influence on the stellar cusp (e.g., core scouring and Brownian motion of the cusp) is ignored in the model.
domain assumption Spherical symmetry and Jeans equilibrium equations are used to compute the velocity dispersion.
The gamma-model is spherical; real galaxies like M87 are triaxial, which can affect the wandering dynamics.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Brownian motion of supermassive black holes in galaxy cores." pith.science (2026). https://pith.science/paper/BV4ALNDK

@misc{pith2026190804283,
  author       = {Pith},
  title        = {Pith review of: Brownian motion of supermassive black holes in galaxy cores},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/BV4ALNDK}},
  note         = {Machine review of arXiv:1908.04283}
}

abstract

We investigate the dynamics of supermassive black holes (SMBHs) in galactic cores by means of a semi-analytic model based on the Langevin equation, including dynamical friction and stochastic noise accounting for the gravitational interactions with stars. The model is validated against direct $N$-body simulations of intermediate-mass black holes in stellar clusters where a realistic number of particles is accessible. For the galactic case, we find that the SMBH experiences a Brownian-like motion with a typical displacement from the geometric center of the Galaxy of a few parsecs, for system parameters compatible with M87. \keywords{stellar dynamics, black hole physics, methods: n-body simulations, methods: statistical.

Figures

Figures reproduced from arXiv: 1908.04283 by the authors.

**Figure 1.** Top row: orbit projections in the (x, y) plane for a BH initially at rest at the centre of an isotropic Plummer model with N = 104 particles in a direct N−body integration (left), and stochastic simulation with Gaussian (middle) and Holtsmark (right) random force distribution. In all cases MBH/m∗ = 100. Bottom row: distribution of positions along x attained by the BH (left) and associated distributions of vx (right)… view at source ↗

**Figure 2.** Evolution of the distance from the galactic center of SMBH of MBH = 6 × 109M in a Mgal = 3 × 1012M γ−model with γ = 1.2 obtained with Holtsmark noise. starting at rest in a galaxy with stellar and dark matter density distributions (2.1), total mass Mgal = 3 × 1012M , γ = 1.2 (for both stellar and dark components), rc = 3 kpc, and a dark to visible matter ratio of ≈ 6, [parameters roughly corresponding to the case of… view at source ↗

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

On an empirical method to build near-equilibrium axisymmetric and triaxial galaxy models
astro-ph.GA 2026-07 conditional novelty 5.0 of 10

An improved adiabatic-squeezing method with an Ornstein-Uhlenbeck noise term lets modelers prescribe the final axial ratios of near-equilibrium triaxial N-body galaxy models.

Reference graph

Works this paper leans on

21 extracted references · 21 canonical work pages · cited by 1 Pith paper

[1]

Alessandrini E., Lanzoni B., Miocchi P., Ciotti L., Ferraro F. R. 2014, ApJ, 795, 169

work page 2014
[2]

2014a, MNRAS, 444, 3738

Arca-Sedda M., Capuzzo-Dolcetta R. 2014a, MNRAS, 444, 3738

work page
[3]

2014b, ApJ, 785, 51

Arca-Sedda M., Capuzzo-Dolcetta R. 2014b, ApJ, 785, 51

work page
[4]

2010, ApJ Letters, 717, L6

Batcheldor D., Robinson A., Axon D.J., Perlman E.S., Merritt, D. 2010, ApJ Letters, 717, L6

work page 2010
[5]

2008, Galactic Dynamics, Princeton University Press NJ 2nd ed

Binney J., Tremaine S. 2008, Galactic Dynamics, Princeton University Press NJ 2nd ed

work page 2008
[6]

2002, EPL, 57, 315

Bottaccio M., Amici A., Miocchi P., Capuzzo-Dolcetta R., Montuori M., Pietronero L. 2002, EPL, 57, 315

work page 2002
[7]

1943, ApJ, 97, 255

Chandrasekhar S. 1943, ApJ, 97, 255

work page 1943
[8]

1949, Reviews of Modern Physics, 21, 383

Chandrasekhar S. 1949, Reviews of Modern Physics, 21, 383

work page 1949

Show all 21 references

[9]

1942, ApJ, 95, 489

Chandrasekhar S., von Neumann J. 1942, ApJ, 95, 489

1942
[10]

1943, ApJ, 97, 1

Chandrasekhar S., von Neumann J. 1943, ApJ, 97, 1

1943
[11]

and Loeb A

Chatterjee P., Hernquist L. and Loeb A. 2002 ApJ, 572, 371

2002
[12]

The Event Horizon Telescope Collaboration 2019, ApJ Letters, 875, L1

2019
[13]

1999, EPL, 46, 127

Gabrielli A., Sylos Labini F., Pellegrini S. 1999, EPL, 46, 127

1999
[14]

Gebhardt K. et al. 2011 ApJ, 729, 13

2011
[15]

1919, Annalen der Physik, 363, 577

Holtsmark J. 1919, Annalen der Physik, 363, 577

1919
[16]

1986, Journal of Quantitative Spectroscopy and Radiative Transfer, 36, 1

Hummer D.G. 1986, Journal of Quantitative Spectroscopy and Radiative Transfer, 36, 1

1986
[17]

Kandrup H. E. 1980, Phys. Rep., 63, 1

1980
[18]

2004, Phys

Mannella R. 2004, Phys. Rev. E, 69, 041107

2004
[19]

2015, ApJ, 804, 52

Merritt D. 2015, ApJ, 804, 52

2015
[20]

Petrovskaya I. V. 1986, Soviet Astronomy Letters, 12, 237

1986
[21]

2006, ApJ, 643, 210

Wu, X., Tremaine, S. 2006, ApJ, 643, 210

2006

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.

[1] [1]

Alessandrini E., Lanzoni B., Miocchi P., Ciotti L., Ferraro F. R. 2014, ApJ, 795, 169

work page 2014

[2] [2]

2014a, MNRAS, 444, 3738

Arca-Sedda M., Capuzzo-Dolcetta R. 2014a, MNRAS, 444, 3738

work page

[3] [3]

2014b, ApJ, 785, 51

Arca-Sedda M., Capuzzo-Dolcetta R. 2014b, ApJ, 785, 51

work page

[4] [4]

2010, ApJ Letters, 717, L6

Batcheldor D., Robinson A., Axon D.J., Perlman E.S., Merritt, D. 2010, ApJ Letters, 717, L6

work page 2010

[5] [5]

2008, Galactic Dynamics, Princeton University Press NJ 2nd ed

Binney J., Tremaine S. 2008, Galactic Dynamics, Princeton University Press NJ 2nd ed

work page 2008

[6] [6]

2002, EPL, 57, 315

Bottaccio M., Amici A., Miocchi P., Capuzzo-Dolcetta R., Montuori M., Pietronero L. 2002, EPL, 57, 315

work page 2002

[7] [7]

1943, ApJ, 97, 255

Chandrasekhar S. 1943, ApJ, 97, 255

work page 1943

[8] [8]

1949, Reviews of Modern Physics, 21, 383

Chandrasekhar S. 1949, Reviews of Modern Physics, 21, 383

work page 1949