pith. sign in

arxiv: 1907.02389 · v1 · pith:3SHGGRX2new · submitted 2019-07-04 · ❄️ cond-mat.soft · cond-mat.stat-mech

Dynamical collective memory in fluidized granular materials

A. Plati , A. Baldassarri , A. Gnoli , G. Gradenigo , A. Puglisi This is my paper

Pith reviewed 2026-05-25 09:11 UTC · model grok-4.3

classification ❄️ cond-mat.soft cond-mat.stat-mech
keywords granular materialscollective rotationsuperdiffusiondynamical heterogeneityvibrated granular mediaprobe diffusioncollective memory
0
0 comments X

The pith

A persistent collective rotational mode in dense vibrated granular media causes the superdiffusive motion of an embedded probe particle.

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

The paper sets up simulations of a vibrated granular medium that match experimental probe diffusion data across densities. It finds that at high density and low granular temperature a macroscopic fraction of the grains begins to rotate slowly as a collective unit, with direction switches occurring only after very long times. This rotational mode directly produces the probe's superdiffusion at large times. The same motion shows dynamical heterogeneity during the intermediate cage regime and then a sharp drop in fluctuations once superdiffusion appears. A reader would care because the finding supplies a concrete mechanism linking collective particle motion to anomalous transport in a driven, dissipative system.

Core claim

The central claim is that at high density and low granular temperature a persistent collective rotational mode emerges: a macroscopic fraction of the medium slowly rotates and randomly switches direction after very long times. This rotational mode of the host medium is the origin of the probe's superdiffusion. Collective motion is accompanied by a kind of dynamical heterogeneity at intermediate times in the cage stage, followed by a strong reduction of fluctuations at late times when superdiffusion sets in.

What carries the argument

persistent collective rotational mode of the granular medium

If this is right

  • The probe's superdiffusion is produced by the host medium's rotational mode rather than by single-particle dynamics.
  • Dynamical heterogeneity appears specifically during the intermediate cage stage.
  • Fluctuations in the medium decrease sharply once the superdiffusive regime begins.
  • The rotational mode appears only above a threshold density and below a threshold granular temperature.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Similar slow collective rotation could be searched for in other driven dissipative systems such as vibrated colloids or active suspensions at comparable packing fractions.
  • The reduction of fluctuations at late times suggests a possible route to long-term memory or ordering in granular flows that could be tested by varying container geometry.
  • If the mode persists in three-dimensional realizations, it would offer a simple explanation for anomalous diffusion reported in some industrial granular mixers.

Load-bearing premise

The numerical simulation quantitatively reproduces the experimental observations of the probe.

What would settle it

A simulation or experiment that shows probe superdiffusion while the granular medium exhibits no macroscopic rotational motion would falsify the claim.

Figures

Figures reproduced from arXiv: 1907.02389 by A. Baldassarri, A. Gnoli, A. Plati, A. Puglisi, G. Gradenigo.

Figure 1
Figure 1. Figure 1: A: setup of the experiment and of the simulation. B [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Dynamics of the probe. A: msd in simulations with [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Comparison probe vs. collective rotation of the [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Collective variable Θc(t) at two different values of N and Γ = 39.8, compared with its “slow component”, see text for definition. A) Msd; B) Power spectra of Ωc. whole simulation. We also run longer simulations, not shown here, confirming that sudden turns occur also at N = 2600, but with tcoll 103 seconds. The interesting connection between spatial rearrangements and changes of rotation speed or direction… view at source ↗
Figure 6
Figure 6. Figure 6: a) Two grains are going to collide with given linear and angular velocities. b) The contact between the grains [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Distributions of displacements for the collective position Θ [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Effect of reducing the diameter of the particles by a factor 2 and increasing [PITH_FULL_IMAGE:figures/full_fig_p009_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Left: evolution with time of the percentage of crystallized particles (see text for the definition). Right: mean-squared [PITH_FULL_IMAGE:figures/full_fig_p010_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Left: effect of inelasticity in the interaction among particles at constant granular temperature. The collective [PITH_FULL_IMAGE:figures/full_fig_p010_10.png] view at source ↗
read the original abstract

Recent experiments with rotational diffusion of a probe in a vibrated granular media revealed a rich scenario, ranging from the dilute gas to the dense liquid with cage effects and an unexpected superdiffusive behavior at large times. Here we setup a simulation that reproduces quantitatively the experimental observations and allows us to investigate the properties of the host granular medium, a task not feasible in the experiment. We discover a persistent collective rotational mode which emerges at high density and low granular temperature: a macroscopic fraction of the medium slowly rotates, randomly switching direction after very long times. Such a rotational mode of the host medium is the origin of probe's superdiffusion. Collective motion is accompanied by a kind of dynamical heterogeneity at intermediate times (in the cage stage) followed by a strong reduction of fluctuations at late times, when superdiffusion sets in.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper sets up a numerical simulation of vibrated granular media that is asserted to quantitatively reproduce experimental observations of a probe particle's rotational diffusion across regimes from dilute gas to dense liquid (including cage effects and long-time superdiffusion). Using the simulation, the authors identify a persistent collective rotational mode—a macroscopic fraction of the medium rotating coherently with rare, long-time direction switches—at high density and low granular temperature, and attribute the probe's superdiffusion to this mode. Collective motion is also linked to dynamical heterogeneity at intermediate times followed by reduced fluctuations at late times.

Significance. If the quantitative match between simulation and experiment holds across all regimes (especially long-time superdiffusion), the work supplies a concrete mechanistic origin for the observed superdiffusion via an emergent collective rotational mode that is inaccessible in experiment. The simulation approach is credited for enabling extraction of host-medium properties and for direct comparison to independent experimental data rather than internal fitting.

major comments (2)
  1. [Abstract and results on validation] Abstract and results section on simulation-experiment comparison: the central claim that the rotational mode is the origin of superdiffusion requires the simulation to quantitatively reproduce experimental probe statistics (MSD and related observables) in the late-time superdiffusive regime. The manuscript asserts this reproduction but does not specify the compared observables, quantitative tolerances, parameter-tuning protocol, or statistical measures of agreement; without these, it is impossible to assess whether the match is robust or whether the identified rotation could be a simulation artifact.
  2. [Discussion of collective rotational mode] Section discussing the rotational mode and its causal link: while the mode is identified, the manuscript must demonstrate that it drives the probe motion (e.g., via velocity cross-correlations between probe and medium or by showing that suppressing the mode eliminates superdiffusion). The current argument rests on temporal coincidence rather than an explicit causal test.
minor comments (2)
  1. [Methods] Notation for granular temperature and density parameters should be defined explicitly on first use and kept consistent with experimental definitions.
  2. [Figures] Figure captions for MSD and rotation plots should include the number of independent runs and error estimation method.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments. We address the two major points below and will revise the manuscript to strengthen the presentation of the simulation validation and the evidence linking the rotational mode to superdiffusion.

read point-by-point responses

  1. Referee: [Abstract and results on validation] Abstract and results section on simulation-experiment comparison: the central claim that the rotational mode is the origin of superdiffusion requires the simulation to quantitatively reproduce experimental probe statistics (MSD and related observables) in the late-time superdiffusive regime. The manuscript asserts this reproduction but does not specify the compared observables, quantitative tolerances, parameter-tuning protocol, or statistical measures of agreement; without these, it is impossible to assess whether the match is robust or whether the identified rotation could be a simulation artifact.

    Authors: We agree that the manuscript would benefit from explicit details on the quantitative comparison. The original text states that the simulation reproduces the experimental observations but does not list the precise observables, tolerances, or fitting protocol. In the revised manuscript we will add a dedicated paragraph (or subsection) specifying: the observables compared (full MSD curves in all regimes, including the long-time superdiffusive exponent; rotational diffusion coefficients; cage-effect signatures), the parameter-tuning protocol (matching packing fraction and granular temperature to the experimental values via vibration amplitude and frequency), and quantitative agreement measures (relative deviation in the superdiffusive slope <15% across independent runs, with error bands from multiple realizations). This will allow readers to judge the robustness of the match. revision: yes

  2. Referee: [Discussion of collective rotational mode] Section discussing the rotational mode and its causal link: while the mode is identified, the manuscript must demonstrate that it drives the probe motion (e.g., via velocity cross-correlations between probe and medium or by showing that suppressing the mode eliminates superdiffusion). The current argument rests on temporal coincidence rather than an explicit causal test.

    Authors: The current manuscript identifies the mode through direct averaging of particle velocities and notes its temporal coincidence with the onset of superdiffusion. We acknowledge that an explicit causal demonstration is desirable. In the revision we will add velocity cross-correlation functions between the probe and the surrounding particles, which show sustained positive correlations precisely during the coherent-rotation intervals. However, performing additional simulations in which the rotational mode is deliberately suppressed (e.g., via altered boundaries or constraints) lies outside the scope of the present computational study; such tests would require a new set of runs and are not feasible within the current work. The added correlation analysis, together with the existing evidence, will provide a stronger mechanistic link. revision: partial

Circularity Check

0 steps flagged

No circularity: simulation validated against independent experiment; rotational mode extracted from validated model.

full rationale

The paper's central step is setting up a numerical simulation whose output is asserted to match experimental probe statistics quantitatively, after which the simulation is used to identify a collective rotational mode inaccessible to direct experiment. No derivation reduces a claimed prediction to a fitted parameter by construction, no self-citation chain is invoked to justify uniqueness, and no ansatz is smuggled via prior work. The match to external data supplies independent content, so the derivation chain remains self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Only the abstract is available; no explicit free parameters, axioms, or invented entities can be extracted.

pith-pipeline@v0.9.0 · 5680 in / 1052 out tokens · 25215 ms · 2026-05-25T09:11:47.740304+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

48 extracted references · 48 canonical work pages

  1. [1]

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

    work page 1996

  2. [2]

    Andreotti, Y

    B. Andreotti, Y. Forterre, and O. Pouliquen, Granular Media (Cambridge University Press, 20013)

    work page

  3. [3]

    Brilliantov, P

    N. Brilliantov, P. Krapivsky, A. Bodrova, F. Spahn, H. Hayakawa, V. Stadnichuk, and J. Schmidt, Proc. Nat. Acad. Sci. USA 112, 9536 (2015)

    work page 2015

  4. [4]

    P¨ oschel and N

    T. P¨ oschel and N. V. Brilliantov, eds.,Granular Gas Dy- namics (Springer, Berlin, Germany, 2003) pp. 293–316

    work page 2003

  5. [5]

    A. R. Abate and D. J. Durian, Phys. Rev. Lett. 101, 245701 (2008)

    work page 2008

  6. [6]

    C. Ness, R. Mari, and M. E. Cates, Sci. Adv. 4, eaar3296 (2018)

    work page 2018

  7. [7]

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

    work page 2004

  8. [8]

    Puglisi, Transport and Fluctuations in Granular Flu- ids (Springer-Verlag, 2015)

    A. Puglisi, Transport and Fluctuations in Granular Flu- ids (Springer-Verlag, 2015)

    work page 2015

  9. [9]

    Forterre and O

    Y. Forterre and O. Pouliquen, Annu. Rev. Fluid Mech. 40, 1 (2008)

    work page 2008

  10. [10]

    Puglisi, A

    A. Puglisi, A. Gnoli, G. Gradenigo, A. Sarracino, and D. Villamaina, J. Chem. Phys. 014704, 136 (2012)

    work page 2012

  11. [11]

    B. Kou, Y. Cao, J. Li, C. Xia, Z. Li, H. Dong, A. Zhang, J. Zhang, W. Kob, and Y. Wang, Nature 551, 360 (2017)

    work page 2017

  12. [12]

    Puglisi, A

    A. Puglisi, A. Baldassarri, and A. Vulpiani, J. Stat. Mech. 2007, P08016 (2007)

    work page 2007

  13. [13]

    Sarracino, D

    A. Sarracino, D. Villamaina, G. Gradenigo, and A. Puglisi, Europhys. Lett. 92, 34001 (2010)

    work page 2010

  14. [14]

    Gnoli, A

    A. Gnoli, A. Puglisi, A. Sarracino, and A. Vulpiani, Plos One 9, e93720 (2014)

    work page 2014

  15. [15]

    Feitosa and N

    K. Feitosa and N. Menon, Physical review letters 88, 198301 (2002)

    work page 2002

  16. [16]

    Scalliet, A

    C. Scalliet, A. Gnoli, A. Puglisi, and A. Vulpiani, Phys. Rev. Lett. 114, 198001 (2015)

    work page 2015

  17. [17]

    Lasanta and A

    A. Lasanta and A. Puglisi, J. Chem. Phys. 143, 064511 (2015)

    work page 2015

  18. [18]

    Briand, M

    G. Briand, M. Schindler, and O. Dauchot, Phys. Rev. Lett. 120, 208001 (2018)

    work page 2018

  19. [19]

    B. Kou, Y. Cao, J. Li, C. Xia, Z. Li, H. Dong, A. Zhang, J. Zhang, W. Kob, and Y. Wang, Phys. Rev. Lett. 121, 018002 (2018)

    work page 2018

  20. [20]

    Plimpton, J

    S. Plimpton, J. Comp. Phys 117, 1 (1995)

    work page 1995

  21. [21]

    [20, 25, 35–47]

    See Supplemental Material at [URL will be inserted by publisher] for details on analytical and numerical results, which includes Refs. [20, 25, 35–47]

    work page

  22. [22]

    Sarracino, D

    A. Sarracino, D. Villamaina, G. Costantini, and A. Puglisi, J. Stat. Mech. , P04013 (2010)

    work page 2010

  23. [23]

    Gnoli, A

    A. Gnoli, A. Puglisi, and H. Touchette, Europhys. Lett. 102, 14002 (2013)

    work page 2013

  24. [24]

    Our analysis revealed that the behavior of the granu- lar medium is largely independent from the presence or absence of the blade, however a blade-free simulation seemed us cleaner with respect to assessing the minimal ingredients of the observed dynamics

    work page

  25. [25]

    Lechenault, R

    F. Lechenault, R. Candelier, O. Dauchot, J.-P. Bouchaud, and G. Biroli, Soft Matt. 6, 3059 (2010)

    work page 2010

  26. [26]

    M. D. Ediger, Annu. Rev. Phys. Chem. 51, 99 (2000)

    work page 2000

  27. [27]

    Berthier, G

    L. Berthier, G. Biroli, J.-P. Bouchaud, L. Cipelletti, and W. van Saarloos, Dynamical heterogeneities in glasses, colloids, and granular media , Vol. 150 (OUP Oxford, 2011)

    work page 2011

  28. [28]

    Dauchot, G

    O. Dauchot, G. Marty, and G. Biroli, Phys. Rev. Lett. 95, 265701 (2005)

    work page 2005

  29. [29]

    A. S. Keys, A. R. Abate, S. C. Glotzer, and D. J. Durian, Nature Physics 3, 260 (2007)

    work page 2007

  30. [30]

    Candelier, O

    R. Candelier, O. Dauchot, and G. Biroli, Phys. Rev. Lett. 102, 088001 (2009)

    work page 2009

  31. [31]

    D. V. T. G. G. Gradenigo, A. Sarracino and A. Puglisi, J. Stat. Mech. L12002 (2010)

    work page 2010

  32. [32]

    Cavagna, I

    A. Cavagna, I. Giardina, and T. S. Grigera, Phys. Rep. 728, 1 (2018)

    work page 2018

  33. [33]

    Manacorda and A

    A. Manacorda and A. Puglisi, Phys. Rev. Lett. 119, 208003 (2017)

    work page 2017

  34. [34]

    Manacorda, Lattice Models for Fluctuating Hydrody- namics in Granular and Active Matter (Springer, 2018)

    A. Manacorda, Lattice Models for Fluctuating Hydrody- namics in Granular and Active Matter (Springer, 2018)

    work page 2018

  35. [35]

    H. P. Zhang and H. A. Makse, Phys. Rev. E 72, 1 (2005)

    work page 2005

  36. [36]

    L. E. Silbert, D. Erta, G. S. Grest, T. C. Halsey, D. Levine, and S. J. Plimpton, Phys. Rev. E64, 1 (2001)

    work page 2001

  37. [37]

    N. V. Brilliantov, F. Spahn, J. M. Hertzsch, and T. P¨ oschel, Phys. Rev. E53, 5382 (1996)

    work page 1996

  38. [38]

    V. L. Popov, Contact Mechanics and Friction (Springer- Verlag, Berlin, 2010)

    work page 2010

  39. [39]

    Di Renzo and F

    A. Di Renzo and F. P. Di Maio, Chem. Eng. Sci. 59, 525 (2004)

    work page 2004

  40. [40]

    T. S. T. P¨ oschel, Computational Granular Dynamics (Springer, Berlin, 2005)

    work page 2005

  41. [41]

    Rackl and K

    M. Rackl and K. J. Hanley, Powd. Tech. 307, 73 (2017)

    work page 2017

  42. [42]

    Lechner and C

    W. Lechner and C. Dellago, J. Chem. Phys. 129, 114707 (2008)

    work page 2008

  43. [43]

    Russo and H

    J. Russo and H. Tanaka, Sci. Reps. 2, 505 (2012)

    work page 2012

  44. [44]

    Zaccarelli, C

    E. Zaccarelli, C. Valeriani, E. Sanz, W. C. K. Poon, M. E. Cates, and P. N. Pusey, Phys. Rev. Lett. 103, 135704 (2009)

    work page 2009

  45. [45]

    Pusey, E

    P. Pusey, E. Zaccarelli, C. Valeriani, E. Sanz, W. C. Poon, and M. E. Cates, Philos. Trans. Royal Soc. A 367, 4993 (2009)

    work page 2009

  46. [46]

    Hydrodynamic long-time tails,

    W. Brenig, “Hydrodynamic long-time tails,” in Statistical Theory of Heat (Springer, 1989) p. 149

    work page 1989

  47. [47]

    Fiege, T

    A. Fiege, T. Aspelmeier, and A. Zippelius, Phys. Rev. Lett. 102, 098001 (2009). 6 Figure 6. a) Two grains are going to collide with given linear and angular velocities. b) The contact between the grains starts with a normal and tangential relative velocity ⃗ gN 12 and⃗ gT

    work page 2009

  48. [48]

    connected

    c) An example of surface deformation and definition of the instantaneous normal compression ξ(t). d) Microscopic interpretation of the model: superficial asperities define a normal and a tangential surfaces of contact that enable respectively the transmission of impulse and angular momentum [37]. SM: SUPPLEMENT AL MA TERIAL INTERACTION MODEL FOR THE GRAIN CO...

    work page 2000