Pith. sign in

REVIEW 5 major objections 5 minor 25 references

Aging Dynamics and Velocity Field Correlations in Three-Dimensional Uniformly Heated Granular Gases: A Molecular Dynamics Study

T0 review · 5 major / 5 minor · reviewed 2026-08-09 · deepseek-v4-flash

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

Pith's one-line read Velocity correlations in a uniformly heated granular gas age with waiting time.

desk verdict The aging claim rests on plotting against absolute time; without a fixed-lag analysis, the central conclusion is an artifact. read the letter →

arxiv 2501.19197 v1 pith:YOOY7XYF submitted 2025-01-31 physics.comp-ph cond-mat.stat-mech

Rameez Farooq Shah , Shikha Kumari , Syed Rashid Ahmad This is my paper · ORCID
classification physics.comp-phcond-mat.stat-mech
keywords aginggranulargasvelocityautocorrelationmoleculardynamicsevent-drivensimulationwhite-noisethermostatnonequilibriumsteadystateinelastichardspheres
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 uses event-driven molecular dynamics to simulate a three-dimensional gas of inelastic hard spheres kept out of equilibrium by a white-noise thermostat. The authors compute the two-time velocity autocorrelation function $C(\tau_w,\tau)$ and find that its normalized decay becomes slower as the waiting time $\tau_w$ grows, so the correlation depends on $\tau_w$ and $\tau$ separately rather than only on the elapsed time $\tau-\tau_w$. They take this explicit waiting-time dependence as evidence that the system exhibits aging. The result matters because it shows that a driven dissipative gas retains a growing memory of its earlier velocities, a signature usually associated with glasses and other non-equilibrium systems, and it offers a numerical benchmark for theories of granular hydrodynamics and aging.

What carries the argument

The central object is the two-time velocity autocorrelation function $C(\tau_w,\tau)=\frac{1}{N}\sum_{i}\mathbf{v}_i(\tau_w)\cdot\langle\mathbf{v}_i(\tau)\rangle$, evaluated after the system reaches its nonequilibrium steady state. The simulation uses event-driven molecular dynamics for $N=500{,}000$ identical inelastic hard spheres with periodic boundaries, and a white-noise thermostat implemented by adding a random velocity kick $\sqrt{A}\,\sqrt{\mathrm{d}t}\,\xi$ to each particle at every time step, followed by a shift to the center-of-mass frame. The diagnostic is the set of decay curves of $C(\tau_w,\tau)$ at fixed $\tau_w$: an equilibrium system would collapse them to a function of $\tau-\tau_w$, while here they separate by $\tau_w$, which the paper reads as aging. Haff's law for the cooling stage sets the temperature scale and frames the comparison with the freely cooling gas.

What would settle it

Measure the steady-state velocity distribution and the two-time autocorrelation using a true Gaussian white-noise thermostat and repeat the run with several significantly smaller time steps. If the normalized $C(\tau_w,\tau)$ curves collapse onto a single function of $\tau-\tau_w$, or if the decay rate changes systematically with $\mathrm{d}t$, the aging signature is not intrinsic to the heated granular gas. A collapse onto a master curve in $\tau/\tau_w$ would instead confirm genuine aging.

Watch

Extended reading notes

Core claim

The central discovery is that in a three-dimensional uniformly heated granular gas, the normalized velocity autocorrelation function decays exponentially at short times and then crosses over to a slower-than-exponential decay, with the decay rate decreasing as the waiting time $\tau_w$ increases. This two-time structure means the system's velocity memory cannot be described by a function of the time difference alone; the explicit dependence on $\tau_w$ is the paper's evidence of aging. The effect is most pronounced for the most inelastic case studied, $r=0.80$, where the decay is slowest, while weaker dissipation ($r=0.95$) shows faster decorrelation. The authors compare with the known behavior of freely cooling gases and with systems governed by a constant restitution coefficient, and report that the heated system's decay is consistently faster than in the constant-$r$ case while still displaying aging.

Load-bearing premise

The simulation heats the gas by adding uniform random velocity kicks at each discrete time step $\mathrm{d}t$ and treats this as equivalent to the Gaussian white-noise thermostat in the equations of motion, without specifying $\mathrm{d}t$ or showing that the steady state is the same; if this equivalence fails, the observed waiting-time dependence could be an artifact of the heating protocol rather than a genuine aging property of the granular gas.

Editorial extensions

If this is right

  • If the aging claim holds, two-time velocity correlations in a heated granular gas cannot be reduced to a function of elapsed time alone; the system's memory grows with its age.
  • Stronger dissipation (lower restitution coefficient $r$) produces a longer-lived velocity memory, so aging is more pronounced in more inelastic gases.
  • The early-time exponential decay followed by a slower-than-exponential tail means no single relaxation time describes the dynamics.
  • Because the decay with a variable, cooling-induced effective restitution is faster than with a fixed restitution coefficient, the heating protocol itself enters the aging rate and must be controlled in comparisons.

Reading between the lines

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

  • One could test whether the aging data obey a scaling collapse, for example $C(\tau_w,\tau)/C(\tau_w,\tau_w)$ as a function of $\tau/\tau_w$; the paper does not report such a collapse, but if it works it would connect this system to the aging phenomenology of glasses and spin glasses.
  • A direct measurement of the spatial velocity correlation length as a function of $\tau_w$ would test the paper's attribution of slow decay to growing velocity-field correlations; if the correlation length stays constant while the decay slows, the proposed mechanism would need revision.
  • Because the thermostat uses uniform rather than Gaussian kicks, repeating the simulation with a proper Gaussian white-noise integrator would separate heating-protocol effects from physics; the unspecified time step $\mathrm{d}t$ leaves convergence in step size as an open check.
  • If the aging is genuine, effective transport coefficients extracted from one-time averages in driven granular gases may depend on preparation time, which could matter for industrial and astrophysical applications of granular flow.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

5 major / 5 minor

Summary. The paper reports event-driven molecular dynamics simulations of 500,000 inelastic hard spheres in a cubic box with periodic boundary conditions and a white-noise thermostat. It defines a two-time velocity autocorrelation function C(τw, τ) and presents normalized curves for restitution coefficients r = 0.80, 0.85, 0.90, and 0.95 and waiting times τw = 10, 50, 100, and 200. The authors observe that the decay of the normalized autocorrelation becomes slower as τw increases and interpret this as evidence of aging. They also compare the decay with predictions for a constant restitution coefficient and conclude that the system retains stronger velocity memory at larger waiting times, especially for strong dissipation.

Significance. If the aging claim were established, the paper would provide useful numerical evidence for two-time velocity correlations in uniformly heated granular gases and would extend existing work on the homogeneous cooling state. The strengths of the manuscript are the large system size (N = 500,000), the average over 50 independent initial conditions, and the systematic variation of the restitution coefficient. However, the quantitative evidence for the central claim is currently missing: no fixed-lag collapse test is performed, no statistical uncertainties are shown, and the thermostat implementation is not clearly connected to the stated Gaussian white-noise model. The claim is therefore plausible but not demonstrated as presented; with additional analysis, the manuscript could become publishable.

major comments (5)
  1. [Section IV, Figs. 3–4] The central aging claim is not supported by the plotted data as presented. The normalized autocorrelation C̄(τw, τ)/C̄(τw, τw) is shown against the absolute time τ, so curves with different τw are horizontally shifted relative to one another. For any time-translation-invariant process, C(τw, τ) = C(0, τ − τw), and plotting against τ will produce different-looking curves whenever the decay is not purely exponential; the paper itself reports exponential early decay followed by slower-than-exponential decay. The observed widening separation is therefore exactly what a stationary process would produce, and it cannot be used as evidence of aging. To establish aging, the authors must show that at fixed lag Δ = τ − τw the correlation depends on τw, or equivalently that the curves do not collapse when plotted versus Δ. They should also verify that the system has reached a steady state before the earliest waiting time used; otherwise the τw dependence is transient relaxation toward the steady state rather than steady-state aging.
  2. [Section IV, Eq. (14)] The thermostat implementation is not equivalent to the stated Gaussian white-noise force of Eqs. (8)–(10) unless additional conditions are met. Equation (14) adds √A√dt ξ with ξ uniformly distributed on [−1/2, 1/2], whose variance is 1/12 rather than 1, and the time step dt is never specified. A fixed-timestep update is also at tension with the event-driven algorithm, which normally advances collision-to-collision. The authors should report dt, state how the uniform noise is mapped to the Gaussian noise (e.g., by matching second moments or by using a Gaussian random number generator), and justify that the resulting steady state is the same as that of the stated stochastic equations.
  3. [Section IV, paragraph after Eq. (15)] The text contains a direct contradiction about the restitution coefficient. It states that 'the system is not evolving with a constant coefficient of restitution; instead, the coefficient of restitution decreases as the system evolves due to the cooling of the granular gas,' but the simulations are described as one run for each fixed r = 0.80, 0.85, 0.90, and 0.95. If r is in fact time-dependent, the protocol must be specified; if r is fixed, the sentence should be removed. This ambiguity directly affects the comparison with constant-r results in Section V and must be resolved.
  4. [Section IV, Figs. 3–4] No statistical uncertainties are reported despite the claim of averages over 50 independent initial conditions. On the semilog and linear-log scales used, the differences between curves at successive waiting times can easily be within sampling error. Standard errors, confidence bands, or at least a statement of the run-to-run spread are needed to assess whether the observed τw dependence is significant.
  5. [Section III, Eqs. (12)–(13)] Equation (12) as written, C(τw, τ) = C(τw, τw) exp[(1 + ε)/(2d)(τ − τw)], grows with τ − τw because 1 + ε > 0, and therefore cannot represent the 'decay behavior' described in the text; the sign or normalization is likely wrong. Equation (13) has an unusual nested-exponential structure and is not used explicitly in the figures. These equations should be corrected and either derived or compared quantitatively with the simulation data.
minor comments (5)
  1. [Section III, Eq. (11)] The definition C(τw, τ) = (1/N) Σ_i v_i(τw)·⟨v_i(τ)⟩ is notationally incorrect: the ensemble average should act on the product v_i(τw)·v_i(τ), not on v_i(τ) alone, whose average vanishes by momentum conservation. Also, the bar symbol C̄ used in the figures is never defined.
  2. [Section IV, Figs. 3–4] The text says results for τw = 0 are plotted, but the figure legends show only τw = 10, 50, 100, and 200; please clarify what is actually shown.
  3. [Section II, Eq. (6) and simulation figures] The variable τ is introduced as a cumulative collision count in Eq. (6), but in Eq. (14) and in the figures it is treated as a continuous time; specify the conversion between real time and τ used in the simulations.
  4. [Abstract and Section III] The phrase 'dependence on both waiting time tw and correlation time t in an independent manner' is vague; the precise statement should be that C(τw, τ) is not a function of τ − τw alone.
  5. [Entire manuscript] There are numerous typographical errors (e.g., 'g as' in the abstract) and inconsistent notation between t and τ; the manuscript needs careful proofreading.

Circularity Check

0 steps flagged · score 0.0 of 10

No circular derivation: the aging conclusion is a definitional label applied to a measured correlation, not a fitted parameter or self-citation chain renamed as a prediction.

full rationale

The simulation and analysis chain is self-contained. Haff's law is derived in Eqs. (3)-(7) and then compared with the MD temperature data in Figs. 1-2; no parameter is fitted to the velocity-autocorrelation data and later renamed as a prediction. The aging conclusion is a definitional label: Section III states "The influence of τw on C(τw, τ) is termed aging," so the abstract's claim that the τw-dependence of C serves as "compelling evidence of the aging properties" is a semantic application of that definition, not a derivation that reduces to its own inputs. No Eq. (11)-(14) quantity is defined in terms of the aging conclusion, and no load-bearing self-citation supplies the central result. Concerns about whether the plotted τw-dependence could arise from plotting normalized C(τw,τ) against absolute τ in a stationary system are statistical/interpretive issues, not circular reasoning. Therefore no circularity is present.

Assumptions & free parameters 3 free parameters · 3 assumptions · 0 invented entities

The simulation depends on several unstated modeling choices: the amplitude of the thermostat noise A=0.001 (chosen ad hoc), the density n=0.02 (chosen ad hoc), and the unspecified time step dt for applying the noise. The equivalence between the uniform-kick thermostat in Eq. (14) and the Gaussian white-noise thermostat in Eq. (9) is assumed. The stationarity of the two-time correlations is inferred only from temperature stabilization.

free parameters (3)
  • Noise amplitude A = 0.001
    Chosen ad hoc; the paper does not justify this value or its relation to a target granular temperature.
  • Number density n = 0.02
    Chosen ad hoc; affects collision frequency and the resulting steady state.
  • Thermostat time step dt = not specified in text
    The thermostat is applied every dt, but dt is never stated, leaving the protocol underdetermined.
assumptions (3)
  • domain assumption Uniformly distributed random kicks at each dt are equivalent to Gaussian white noise in the continuum limit.
    Eq. (14) uses uniform random numbers, while Eqs. (9) describe Gaussian noise. The paper does not justify the equivalence for the finite dt used.
  • domain assumption The system reaches a non-equilibrium steady state by tau=1000.
    Fig. 2 shows temperature stabilization, but stationarity of two-time correlations is not demonstrated.
  • domain assumption The density field remains homogeneous throughout the simulation.
    Aging is interpreted via velocity field correlations without checking for clustering or other inhomogeneities.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Aging Dynamics and Velocity Field Correlations in Three-Dimensional Uniformly Heated Granular Gases: A Molecular Dynamics Study." pith.science (2026). https://pith.science/paper/YOOY7XYF

@misc{pith2026250119197,
  author       = {Pith},
  title        = {Pith review of: Aging Dynamics and Velocity Field Correlations in Three-Dimensional Uniformly Heated Granular Gases: A Molecular Dynamics Study},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/YOOY7XYF}},
  note         = {Machine review of arXiv:2501.19197}
}
abstract

We conduct a molecular dynamics simulation of an inelastic gas system utilizing an event-driven algorithm combined with a thermostat mechanism. Initially, the kinetic energy of the system experiences a decay before settling into a non-equilibrium steady state. To explore the aging characteristics, we analyze the velocity autocorrelation function, denoted as \( C(t_w, t) \). Our findings indicate that \( C(t_w, t) \) exhibits a dependence on both waiting time \( t_w \) and correlation time \( t \) in an independent manner. At the outset, \( C(t_w, t) \) demonstrates an exponential decay pattern. With increasing \( t_w \), a slower decay is observed, which can be attributed to the development of correlations in the velocity field. The explicit relationship of \( C(t_w, t) \) with respect to \( t_w \) serves as compelling evidence of the aging properties present in the system. These results deepen our comprehension of non-equilibrium statistical mechanics and the dynamics of dissipative systems. Our research has significant implications for a range of applications involving inelastic collisions, extending from granular materials to phenomena observed in astrophysics. The simulation methodology and the insights gained contribute to the wider field of complex systems, shedding light on the behavior of systems that are far from equilibrium, particularly those characterized by energy dissipation due to inelastic interactions.

Figures

Figures reproduced from arXiv: 2501.19197 by the authors.

Figure 1
Figure 1. FIG. 1: Time dependence of the granular temperature in [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2: Time dependence of the granular temperature in [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3: The normalized velocity autocorrelation function [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: FIG. 4: The normalized velocity autocorrelation function [PITH_FULL_IMAGE:figures/full_fig_p008_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

25 extracted references · 24 canonical work pages

  1. [1]

    H. M. Jaeger, S. R. Nagel, and R. P. Behringer, Rev. Mod. Phys. 68, 1259 (1996)

    1996

  2. [2]

    L. P. Kadanoff, Rev. Mod. Phys. 71, 435 (1999)

    work page 1999

  3. [3]

    I. S. Aranson and L. S. Tsimring, Rev. Mod. Phys. 78, 641 (2006)

    work page 2006

  4. [4]

    Duran, Sands, Powders and Grains (Springer-Verlag, Berlin, 2000)

    J. Duran, Sands, Powders and Grains (Springer-Verlag, Berlin, 2000)

    work page 2000

  5. [5]

    G. H. Ristow, Pattern Formation in Granular Materials (Springer-Verlag, Berlin, 2000)

    work page 2000

  6. [6]

    N. V. Brilliantov and T. P¨ oschel, Kinetic Theory of Granular Gases (Oxford University Press, Oxford, 2004)

    work page 2004

  7. [7]

    P. K. Haff, J. Fluid Mech. 134, 401 (1983)

    work page 1983

  8. [8]

    F. Melo, P. B. Umbanhowar, and H. L. Swinney, Phys. Rev. Lett. 75, 3838 (1995); P. B. Umbanhowar, F. Melo, and H. L. Swinney, Nature (London) 382, 793 (1996). 11

    work page 1995

Show all 25 references
  1. [9]

    Goldhirsch and G

    I. Goldhirsch and G. Zanetti, Phys. Rev. Lett. 70, 1619 (1993); I. Goldhirsch, M. L. Tan, and G. Zanetti, J. Sci. Comput. 8, 1 (1993)

    1993

  2. [10]

    McNamara and W

    S. McNamara and W. R. Young, Phys. Fluids A 4, 496 (1992); Phys. Rev. E 53, 5089 (1996)

    1992

  3. [11]

    J. J. Brey, F. Moreno, and J. W. Dufty, Phys. Rev. E 54, 445 (1996); J. J. Brey, F. Moreno, and M. J. Ruiz-Montero, Phys. Fluids 10, 2965 (1998); 10, 2976 (1998)

    1996

  4. [12]

    T. P. C. van Noije, M. H. Ernst, R. Brito, and J. A. G. Orza, Phys. Rev. Lett. 79, 411 (1997); T. P. C. van Noije, M. H. Ernst, and R. Brito, Phys. Rev. E 57, R4891 (1998)

    1997

  5. [13]

    Luding, M

    S. Luding, M. Huthmann, S. McNamara, and A. Zippelius, Phys. Rev. E 58, 3416 (1998); S. Luding and S. McNamara, Granular Matter 1, 113 (1998); S. Luding and H. J. Herrmann, Chaos 9, 673 (1999)

    1998

  6. [14]

    S. R. Ahmad and S. Puri, Europhys. Lett. 75, 56 (2006); S. R. Ahmad and S. Puri, Phys. Rev. E 75, 031302 (2007)

    2006

  7. [15]

    A. K. Dubey, A. Bodrova, S. Puri, and N. V. Brilliantov, Phys. Rev. E 87, 062202 (2013)

    2013

  8. [16]

    S. K. Das and S. Puri, Europhys. Lett. 61, 749 (2003); S. K. Das and S. Puri, Phys. Rev. E 68, 011302 (2003)

    2003

  9. [17]

    P. Das, S. Puri, and M. Schwartz, Granular Matter 20, 1-11 (2018)

    2018

  10. [18]

    T. P. C. van Noije and M. H. Ernst, Granular Matter 1, 57 (1998)

    1998

  11. [19]

    D. R. M. Williams and F. C. MacKintosh, Phys. Rev. E 54, R9-R12 (1996)

    1996

  12. [20]

    D. R. M. Williams, Physica A 233, 718 (1996)

    1996

  13. [21]

    Chapman and T

    S. Chapman and T. G. Cowling, The Mathematical Theory of Non-Uniform Gases (Cambridge University Press, New York, 1970)

    1970

  14. [22]

    Goldstein and M

    A. Goldstein and M. Shapiro, J. Fluid Mech. 282, 75 (1995)

    1995

  15. [23]

    N. V. Brilliantov and T. P¨ oschel, Europhys. Lett. 74, 424 (2006)

    2006

  16. [24]

    M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids (Oxford University Press, Oxford, 1987)

    1987

  17. [25]

    D. C. Rapaport, The Art of Molecular Dynamics Simulation , Second Edition (Cambridge University Press, Cambridge, 2005). 12

    2005

Pith tools

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