Macroscopic fluctuation theory of interacting Brownian particles
Pith reviewed 2026-05-21 16:44 UTC · model grok-4.3
The pith
Macroscopic fluctuation theory combined with equilibrium collective diffusion gives exact large-scale dynamics for interacting Brownian particles.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
By combining it with standard results of equilibrium statistical mechanics for the collective diffusion coefficient, the MFT gives access to the exact large-scale dynamical properties of the system, both in- and out-of-equilibrium. In particular, we obtain exact results for dynamical correlations between the density and the current of particles.
What carries the argument
Macroscopic fluctuation theory framework incorporating the collective diffusion coefficient from equilibrium statistical mechanics to capture interaction effects at large scales.
If this is right
- Exact results for dynamical correlations between the density and the current of particles in and out of equilibrium.
- Precise description of these correlations for one-dimensional models such as the Calogero and Riesz gases and nearest-neighbor interaction systems.
- Characterization of tracer diffusion with the single-file constraint for arbitrary pairwise interactions.
- Quantitative results for higher-dimensional systems using virial expansion methods.
Where Pith is reading between the lines
- The framework may enable similar exact results for other dynamical quantities in interacting particle systems where equilibrium coefficients are known.
- It suggests a general way to bridge equilibrium statistical mechanics with non-equilibrium fluctuation theories for transport calculations.
- Potential applications include modeling of fluctuations in colloidal suspensions or other many-body systems with pairwise forces.
Load-bearing premise
The collective diffusion coefficient obtained from equilibrium statistical mechanics can be inserted directly into the macroscopic fluctuation theory framework for arbitrary pairwise interactions without additional out-of-equilibrium corrections or validity restrictions.
What would settle it
Comparison of the predicted dynamical correlations with results from direct numerical simulations of the particle trajectories in a specific model, such as the Riesz gas or tethered particles.
Figures
read the original abstract
We apply the macroscopic fluctuation theory (MFT) to study the large-scale dynamical properties of Brownian particles with arbitrary pairwise interaction. By combining it with standard results of equilibrium statistical mechanics for the collective diffusion coefficient, the MFT gives access to the exact large-scale dynamical properties of the system, both in- and out-of-equilibrium. In particular, we obtain exact results for dynamical correlations between the density and the current of particles. For one-dimensional systems, this allows us to obtain a precise description of these correlations for emblematic models, such as the Calogero and Riesz gases, and for systems with nearest-neighbor interactions such as the Rouse chain of hardcore particles or the recently introduced model of tethered particles. Tracer diffusion with the single-file constraint (but for arbitrary pairwise interaction) is also studied. For higher-dimensional systems, we quantitatively characterize these dynamical correlations by relying on standard methods such as the virial expansion.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper applies macroscopic fluctuation theory (MFT) to Brownian particles with arbitrary pairwise interactions. It combines MFT with the collective diffusion coefficient from equilibrium statistical mechanics to claim exact access to large-scale dynamical properties in and out of equilibrium, particularly dynamical correlations between density and current. Applications include 1D models like Calogero and Riesz gases, Rouse chain, tethered particles, single-file tracer diffusion, and higher-D via virial expansion.
Significance. If the results hold, this provides a valuable framework for obtaining exact dynamical correlations in interacting particle systems without additional parameters, extending MFT to general interactions and enabling precise descriptions for important models in statistical mechanics. The approach could be useful for both theoretical understanding and quantitative characterization of non-equilibrium dynamics.
major comments (1)
- The central claim relies on direct insertion of the equilibrium collective diffusion coefficient D(ρ) into the MFT equations to obtain exact out-of-equilibrium density-current correlations for arbitrary pairwise potentials. No explicit derivation is supplied demonstrating that higher-order gradient corrections or interaction-induced memory effects vanish in the hydrodynamic limit (e.g., for Calogero 1/r² or Riesz 1/r^α potentials). This justification is load-bearing for the exactness assertion both in and out of equilibrium.
minor comments (2)
- The abstract and introduction could more clearly delineate which results follow from the MFT-equilibrium combination versus standard hydrodynamic limits.
- For the 1D applications, additional references to prior exact solutions for Calogero and Riesz gases would strengthen the context.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for their positive evaluation of its significance. We address the single major comment below.
read point-by-point responses
-
Referee: The central claim relies on direct insertion of the equilibrium collective diffusion coefficient D(ρ) into the MFT equations to obtain exact out-of-equilibrium density-current correlations for arbitrary pairwise potentials. No explicit derivation is supplied demonstrating that higher-order gradient corrections or interaction-induced memory effects vanish in the hydrodynamic limit (e.g., for Calogero 1/r² or Riesz 1/r^α potentials). This justification is load-bearing for the exactness assertion both in and out of equilibrium.
Authors: We thank the referee for highlighting this point. Macroscopic fluctuation theory is constructed precisely in the hydrodynamic scaling limit, where fields are coarse-grained over lengths much larger than any microscopic scale (particle size or interaction range). In this regime, higher-order gradient corrections are systematically suppressed by powers of the inverse macroscopic length and do not contribute to the leading large-scale correlations. The underlying microscopic dynamics are overdamped Langevin equations, which are Markovian; consequently, no interaction-induced memory kernels survive at the hydrodynamic level. The equilibrium collective diffusion coefficient D(ρ) is obtained from the static structure factor (or equivalently from the second derivative of the free-energy functional) and already encodes the full effect of arbitrary pairwise interactions on the equilibrium fluctuations. Inserting this D(ρ) into the MFT action therefore yields the exact leading-order dynamical correlations without additional parameters. For the long-range cases (Calogero 1/r² and Riesz 1/r^α), the same hydrodynamic argument applies once D(ρ) is computed from the known equilibrium thermodynamics of these models. We will add a concise paragraph in Section II (or a short appendix) that recalls the standard hydrodynamic derivation of MFT for interacting Brownian particles and explicitly states why higher-order and memory corrections are negligible. revision: yes
Circularity Check
No significant circularity; derivation combines independent equilibrium inputs with MFT framework
full rationale
The paper applies the macroscopic fluctuation theory (MFT) to Brownian particles with arbitrary pairwise interactions and inserts the collective diffusion coefficient D(ρ) obtained from standard equilibrium statistical mechanics. This step is presented as yielding exact large-scale dynamical properties and density-current correlations both in and out of equilibrium, including for models such as the Calogero gas and Riesz gas. No quoted equation or derivation step reduces a prediction to a fitted parameter or self-citation by construction; the equilibrium results are invoked as external, independent inputs rather than derived within the paper. The framework for one- and higher-dimensional cases relies on this combination without evidence that out-of-equilibrium results are tautologically forced by the inputs. The derivation chain remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Macroscopic fluctuation theory applies to systems of Brownian particles with arbitrary pairwise interactions.
- domain assumption Equilibrium statistical mechanics supplies the exact collective diffusion coefficient that can be combined with MFT.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
D(ρ)=μ₀P'(ρ) ... fluctuation-dissipation relation 2kBT D(ρ)/σ(ρ)=f''(ρ) ... σ(ρ)=2μ₀kBTρ
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
MFT action S[ρ,H] ... Euler-Lagrange equations (13) for optimal q,p
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Forward citations
Cited by 1 Pith paper
-
Universal tracer statistics in single-file transport
Tracer position statistics in single-file hard-rod transport are universal across diffusive and ballistic dynamics for one-time distributions.
Reference graph
Works this paper leans on
-
[1]
V0(x) = 0, to illustrate the method
Independent particles We first consider the simplest and well-known case, which is the one of independent particles, i.e. V0(x) = 0, to illustrate the method. In this case the partition function (17) becomes ZN,V (β) = 1 N! V ℓ0 N ,(28) where the volume V is the length of the system. Using Stirling’s formula, we obtain that the free energy den- sity (15) ...
-
[2]
For this model, the grand potential has been computed explic- itly [ 59]
The Calogero potential We now consider particles interacting pairwise via the Calogero potential V0(x) = g x2 ,(32) which has been widely studied [ 14, 15, 55–58]. For this model, the grand potential has been computed explic- itly [ 59]. We can thus use this result to compute the diffusion coefficient. Explicitly, we introduce the grand canonical partitio...
-
[3]
Arbitrary nearest-neighbor interactions We now consider a variation of the model (1) in which a particle interacts only with its two nearest neighbors (one on each side), so that the evolution of the positions is now given by dxi dt =−µ 0V ′ 0(xi −x i−1)−µ 0V ′ 0(xi −x i+1) + p 2D0 ηi . (45) This model can be seen as an approximation of the origi- nal mod...
-
[4]
The gas of hard rods We first consider, as a direct illustration of the re- sult (61), the case of the gas of hard rods of length ℓ. This corresponds to a nearest-neighbor interaction (45), with the potential V0(x) = +∞for|x|< ℓ 0 for|x|> ℓ .(64) Applying the formalism of Section III B 3, the diffusion coefficient (61) is expressed in terms of ˆV(β, s) = ...
-
[5]
The gas of sticky hard rods Another classical model of nearest-neighbour interact- ing particles in one dimension is given by the sticky hard rods, corresponding to the interaction potential V0 defined by [60, 64] e−βV0(x) = Θ(x−ℓ) +γδ(x−ℓ),(68) where ℓ is the size of the particles and γ a parameter which controls the adhesiveness of the particles. This p...
-
[6]
The Rouse chain with hardcore repulsion We now consider the important model of the Rouse chain, here with additional and physically relevant hard- core repulsion. The standard Rouse model (without hard- core interaction) has been introduced to model a poly- mer, in which the monomers are attached by harmonic springs [11]. We derive here the collective dif...
-
[7]
The model of tethered particles with hardcore repulsion We consider now the recently introduced model of teth- ered random walkers [67], corresponding to hardcore parti- cles of size ℓ, attached to the next particle with a cable of length ∆ −ℓ > 0, so that the centers of two neighbouring particles cannot be separated by a distance larger than ∆. This corr...
-
[8]
A comment on the effect of the size of the particles Let us consider the case of particles of length ℓ > 0 interacting via nearest-neighbour interaction, with an arbitrary interaction potential V0(x) = +∞for|x|< ℓ , ˜V0(x−ℓ) for|x|> ℓ . (84) 10 0.0 0.2 0.4 0.6 0.8 1.0 0 20 40 60 80 100 D( ) Rouse V(x) = 1 2 x2 at T=1 simulations = 0 = 1 = 2 (a) 0.0 0.2 0....
-
[9]
Virial expansion The virial expansion gives the equation of state P (ρ) as a power series in the densityρ[51], P(ρ) kBT =ρ+B 2(T)ρ 2 +B 3(T)ρ 3 +O(ρ 4),(91) where B2 and B3 are the virial coefficients, which are explicitly given by B2(T) =− 1 2V Z V ddx1 Z V ddx2 e−βV0(x1−x2) −1 ≃ − 1 2 Z ddx1 e−βV0(x1) −1 ,(92) and B3(T) =− 1 3V Z V ddx1 Z V ddx2 Z V ddx...
-
[10]
Virial expansion for the Riesz gas The Riesz gas is a well-studied model of interacting par- ticles (see the review [ 58]), corresponding to the pairwise potential [69] V(x) = g ||x||s .(96) Except for the specific case s = 2 in d = 1, corresponding to the Calogero gas discussed above, the equation of state for the Riesz gas is not known. We can neverthel...
-
[11]
(100) In the case of the Calogero gas (i.e
Γ −d s +O(ρ 2) # . (100) In the case of the Calogero gas (i.e. d = 1 and s = 2), this reduces to the leading order of (44), as expected. We have also checked numerically that the third virial coeffi- cient (93) gives the correct next term of (44), although we have not been able to perform the integrals analytically
-
[12]
Pressure from the pair correlation function Another possibility to determine the pressure is to use the relation with the pair correlation functiong(r), P(ρ) =ρk BT−ρ 2 πd/2 dΓ( d 2) Z ∞ 0 rdV ′ 0(r)g(r)dr ,(101) where d is again the dimension. In this expression, g(r) depends on the density ρ, although it is not written explicitly. Since various approxim...
-
[13]
However, for t̸ = t′ the density presents long-range spatial correlations
Note that for t = t′ we recover the local equilibrium correlations computed at t = t′ = 0 (147), as expected. However, for t̸ = t′ the density presents long-range spatial correlations. The result (159) is for now written in terms of the first order ρ1 of the actual macroscopic density ρ. Using the relation (149) this becomes, ⟨ρ(x, t)ρ(x′, t′)⟩c ≃ Λ→∞ Λ−d...
-
[14]
Spohn,Large scale dynamics of interacting particles (Springer Berlin, Heidelberg, 1991)
H. Spohn,Large scale dynamics of interacting particles (Springer Berlin, Heidelberg, 1991)
work page 1991
-
[15]
M. R. Evans and T. Hanney, J. Phys. A38, R195 (2005)
work page 2005
- [16]
- [17]
-
[18]
T. Chou, K. Mallick, and R. K. Zia, Rep. Prog. Phys. 74, 116601 (2011)
work page 2011
-
[19]
L. Bertini, A. De Sole, D. Gabrielli, G. Jona-Lasinio, and C. Landim, Rev. Mod. Phys.87, 593 (2015)
work page 2015
- [20]
- [21]
-
[22]
K. Mallick, Physica A418, 17 (2015), proceedings of the 13th International Summer School on Fundamental Problems in Statistical Physics
work page 2015
-
[23]
K. Mallick, H. Moriya, and T. Sasamoto, Phys. Rev. Lett. 129, 040601 (2022)
work page 2022
- [24]
-
[25]
J.-F. Dufrˆ eche, O. Bernard, and P. Turq, J. Mol. Liq. 118, 189 (2005), contributions to the 28th International Conference on Solution Chemistry
work page 2005
-
[26]
Q. Liu, S. Huang, and Z. Suo, Phys. Rev. Lett.114, 224301 (2015)
work page 2015
- [27]
- [28]
- [29]
-
[30]
D. S. Dean, J. Phys. A29, L613 (1996)
work page 1996
- [31]
-
[32]
V. D´ emery, O. B´ enichou, and H. Jacquin, New J. Phys. 16, 053032 (2014)
work page 2014
-
[33]
D. S. Dean and R. Podgornik, Phys. Rev. E89, 032117 (2014)
work page 2014
- [34]
- [35]
- [36]
-
[37]
M. Kr¨ uger, A. Solon, V. D´ emery, C. M. Rohwer, and D. S. Dean, J. Chem. Phys.148(2018)
work page 2018
-
[38]
S. Mahdisoltani and R. Golestanian, Phys. Rev. Lett.126, 158002 (2021)
work page 2021
-
[39]
D. Venturelli, P. Illien, A. Grabsch, and O. B´ enichou, arXiv:2411.09326 (2025), 10.48550/arXiv.2411.09326
-
[40]
P. L. Muzzeddu, E. Kalz, A. Gambassi, A. Sharma, and R. Metzler, New J. Phys.27, 033025 (2025)
work page 2025
- [41]
- [42]
- [43]
-
[44]
P. L. Krapivsky, K. Mallick, and T. Sadhu, Phys. Rev. Lett.113, 078101 (2014)
work page 2014
-
[45]
P. L. Krapivsky, K. Mallick, and T. Sadhu, J. Stat. Mech. 2015, P09007 (2015)
work page 2015
-
[46]
P. L. Krapivsky, K. Mallick, and T. Sadhu, J. Stat. Phys. 160, 885 (2015)
work page 2015
- [47]
- [48]
-
[49]
A. Grabsch, A. Poncet, P. Rizkallah, P. Illien, and O. B´ enichou, Sci. Adv.8, eabm5043 (2022)
work page 2022
-
[50]
E. Bettelheim, N. R. Smith, and B. Meerson, Phys. Rev. Lett.128, 130602 (2022)
work page 2022
- [51]
-
[52]
P. Rizkallah, A. Grabsch, P. Illien, and O. B´ enichou, J. Stat. Mech.2023, 013202 (2023)
work page 2023
-
[53]
A. Grabsch, P. Rizkallah, A. Poncet, P. Illien, and O. B´ enichou, Phys. Rev. E107, 044131 (2023)
work page 2023
- [54]
- [55]
- [56]
-
[57]
S. Saha and T. Sadhu, arXiv:2501.03164 (2025), 10.48550/arXiv.2501.03164
-
[58]
A. Grabsch, D. Venturelli, and O. B´ enichou, Phys. Rev. Lett.135, 137102 (2025)
work page 2025
- [59]
-
[60]
C. Van den Broeck, F. Lostak, and H. N. W. Lekkerkerker, J. Chem. Phys.74, 2006 (1981)
work page 2006
- [61]
- [62]
-
[63]
B. U. Felderhof, J. Chem. Phys.131, 064504 (2009)
work page 2009
-
[64]
T. L. Hill,An introduction to statistical thermodynamics (Courier Corporation, 1986)
work page 1986
- [65]
- [66]
-
[67]
Note that the same argument also holds for an under- damped dynamics
- [68]
- [69]
-
[70]
A. P. Polychronakos, J. Phys. A39, 12793 (2006)
work page 2006
- [71]
-
[72]
Classical and quantum partition functions of the Calogero—Moser—Sutherland model,
P. Choquard, “Classical and quantum partition functions of the Calogero—Moser—Sutherland model,” inCalogero— Moser— Sutherland Models, edited by J. F. van Diejen and L. Vinet (Springer New York, New York, NY, 2000) pp. 117–125
work page 2000
- [73]
- [74]
-
[75]
Q.-H. Wei, C. Bechinger, and P. Leiderer, Science287, 625 (2000)
work page 2000
-
[76]
B. Lin, M. Meron, B. Cui, S. A. Rice, and H. Diamant, Phys. Rev. Lett.94, 216001 (2005)
work page 2005
-
[77]
R. J. Baxter, J. Chem. Phys.49, 2770 (1968)
work page 1968
-
[78]
S. Schweers, A. P. Antonov, A. Ryabov, and P. Maass, Phys. Rev. E107, L042102 (2023)
work page 2023
-
[79]
S. B. Yuste and A. Santos, J. Stat. Phys.72, 703 (1993)
work page 1993
-
[80]
S. Bravo Yuste, A. Baumgaertner, and E. Abad, J. Chem. Phys.162, 224114 (2025)
work page 2025
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.