Anomalous Diffusion as Structural Memory: An Extended Structural Dynamics Approach
Pith reviewed 2026-05-20 23:54 UTC · model grok-4.3
The pith
Subdiffusion in biology arises when dynamics of structured molecules are projected onto translation alone.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
When dynamics on the full extended phase space of a structured particle are projected onto the translational subspace alone, a memory kernel emerges from the projection without phenomenological postulate. The subdiffusion exponent is determined by the internal mode spectrum, independently measurable from B-factors, NMR order parameters, or molecular dynamics simulations, without fitting to transport data.
What carries the argument
The projection operation from the complete extended phase space (translational plus orientational plus deformational modes) onto the translational coordinates alone.
If this is right
- Subdiffusion strength correlates with molecular flexibility because more internal modes strengthen the emergent memory.
- Temperature drives a crossover to normal diffusion once thermal energy exceeds the characteristic scale set by internal mode frequencies.
- A non-zero rotation-translation cross-correlation spectrum appears and encodes internal dynamics, identically zero in point-particle models.
- Memory timescales scale as the square of particle size.
Where Pith is reading between the lines
- The same projection logic may generate apparent anomalies in any system whose particles possess internal structure, such as polymers or colloids.
- Varying particle size while holding internal-mode frequencies fixed would provide a clean test of the predicted quadratic scaling of memory time.
- Non-Gaussian displacement statistics observed in crowded media could receive a similar structural-memory explanation once the full phase space is restored.
Load-bearing premise
The internal mode spectrum and other structural parameters can be obtained independently of the transport measurements themselves.
What would settle it
A direct test in which the subdiffusion exponent calculated from an independently measured internal-mode spectrum (via B-factors or NMR) fails to match the observed exponent in the same system.
read the original abstract
Sub-diffusion in biological systems is conventionally treated as anomalous, requiring fractional derivatives, heavy-tailed waiting times, or fitted memory kernels. We argue that this anomaly is an artifact of an incomplete phase space. Standard frameworks model diffusing particles as points. Biological molecules are not points. They are three-dimensional deformable entities whose position, orientation, and internal structure are irreducible physical properties, not modeling conveniences appended to a point mass. Within the Extended Structural Dynamics (ESD) framework, each particle is a primitive structured entity with translational, orientational, and deformational degrees of freedom. When dynamics on this full phase space are projected onto the translational subspace alone, a memory kernel emerges from the projection without phenomenological postulate. The subdiffusion exponent is determined by the internal mode spectrum, independently measurable from B-factors, NMR order parameters, or molecular dynamics simulations, without fitting to transport data. Four falsifiable predictions follow: subdiffusion strength correlates with molecular flexibility; temperature drives crossover to normal diffusion at a characteristic energy scale set by internal mode frequencies; a non-zero rotation-translation cross-correlation spectrum encodes internal dynamics, identically zero in point-particle models; and memory timescales scale as the square of particle size. Quantitative consistency with experimental observations for proteins in crowded media is demonstrated using independently estimated structural parameters. What appears anomalous from the point-particle perspective is the expected behavior of structured matter projected onto an impoverished description. The anomaly is not in the physics. It is in the phase space.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces the Extended Structural Dynamics (ESD) framework, treating biological molecules as structured entities with translational, orientational, and deformational degrees of freedom rather than point particles. It argues that projecting the full phase-space dynamics onto the translational subspace alone generates a memory kernel without phenomenological input, such that the subdiffusion exponent is fixed by the internal mode spectrum (measurable independently from B-factors, NMR order parameters, or MD simulations). Four falsifiable predictions are stated, and quantitative consistency with protein diffusion in crowded media is claimed using independently estimated structural parameters.
Significance. If the projection step yields a memory kernel whose long-time decay is demonstrably fixed by equilibrium structural observables without additional fitted couplings, the work would supply a structural origin for subdiffusion that unifies transport anomalies with independently measurable molecular properties. The emphasis on falsifiable predictions and the avoidance of fractional derivatives or ad-hoc kernels are strengths that could influence modeling of intracellular transport.
major comments (1)
- [Abstract and ESD projection derivation] The central claim (abstract) that the subdiffusion exponent is completely determined by the internal mode spectrum rests on the projection operator producing a memory kernel whose form is fixed solely by equilibrium structural data. In any projection formalism the memory kernel is the autocorrelation of the fluctuating force orthogonal to the translational subspace; this autocorrelation is shaped by the dynamical coupling matrix between translational velocity and internal coordinates. B-factors and NMR order parameters supply equilibrium variances but do not uniquely fix the off-diagonal coupling strengths or relaxation rates. The manuscript must supply the explicit, parameter-free rule (e.g., from a microscopic Hamiltonian or structural constraint) that determines these couplings; absent that rule the exponent remains sensitive to modeling choices that are not transport-independent.
minor comments (2)
- [Abstract] The abstract asserts 'quantitative consistency with experimental observations' using independently estimated parameters but does not identify the specific datasets, the structural parameters employed, or the figure/table that displays the comparison; a brief pointer would improve traceability.
- [Predictions paragraph] The four falsifiable predictions are enumerated but their explicit functional forms (e.g., the expected scaling of memory time with particle size squared) are not given; adding one or two quantitative statements would make the predictions easier to test.
Simulated Author's Rebuttal
We thank the referee for their thorough review and valuable feedback on our manuscript. We have carefully considered the major comment and provide our response below. We believe the concerns can be addressed through clarification in the revised version.
read point-by-point responses
-
Referee: The central claim (abstract) that the subdiffusion exponent is completely determined by the internal mode spectrum rests on the projection operator producing a memory kernel whose form is fixed solely by equilibrium structural data. In any projection formalism the memory kernel is the autocorrelation of the fluctuating force orthogonal to the translational subspace; this autocorrelation is shaped by the dynamical coupling matrix between translational velocity and internal coordinates. B-factors and NMR order parameters supply equilibrium variances but do not uniquely fix the off-diagonal coupling strengths or relaxation rates. The manuscript must supply the explicit, parameter-free rule (e.g., from a microscopic Hamiltonian or structural constraint) that determines these couplings; absent that rule the exponent remains sensitive to modeling choices that are not transport-independent.
Authors: We agree that clarifying the determination of the coupling matrix is important. In the ESD approach, the full phase space includes translational, rotational, and deformational coordinates, with the equilibrium distribution given by the structural data (B-factors provide the variances of internal modes). The projection onto the translational subspace uses the Mori-Zwanzig formalism, where the memory kernel is indeed the autocorrelation of the orthogonal force. The dynamical couplings are fixed by the microscopic model: we model the molecule as a rigid body with internal harmonic modes whose frequencies and eigenvectors are determined from the Hessian of the potential energy, constrained by the molecular structure. These are obtained independently from MD simulations or normal mode analysis, without reference to diffusion data. The relaxation rates follow from the mode frequencies and damping. Thus, the subdiffusion exponent, which depends on the spectrum of these modes, is fixed by structural observables. We will revise the manuscript to include an explicit derivation of the coupling matrix from the structural Hamiltonian in the Methods section to make this parameter-free aspect clearer. revision: yes
Circularity Check
No significant circularity; derivation presented as projection from independent structural inputs
full rationale
The paper asserts that projecting the extended structural dynamics phase space (translation + orientation + deformation) onto the translational subspace yields a memory kernel whose long-time behavior sets the subdiffusion exponent via the internal mode spectrum. This spectrum is claimed to be obtained from B-factors, NMR order parameters, or MD simulations without reference to transport data. No equations, self-citations, or fitted parameters are exhibited in the provided text that reduce the exponent or kernel to a quantity defined in terms of the diffusion measurements themselves. The central claim therefore remains self-contained against external structural benchmarks and does not trigger any of the enumerated circularity patterns.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Biological molecules are three-dimensional deformable entities whose position, orientation, and internal structure are irreducible physical properties.
invented entities (1)
-
Extended Structural Dynamics (ESD) framework
no independent evidence
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
When dynamics on this full phase space are projected onto the translational subspace alone, a memory kernel emerges from the projection without phenomenological postulate. ... Γ(t) = λ²/m ∑ sin(ω_α t)/ω_α
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
the subdiffusion exponent α is determined by the internal mode spectrum ... ρ(ω) = ρ0 ω^β, β ≈ 0.8–1.2
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.
Reference graph
Works this paper leans on
-
[1]
Einstein, A. Über die von der molekularkinetischen Theorie der Wärme geforderte Bewegung von in ruhenden Flüssigkeiten suspendierten Teilchen. Ann. Phys. 322, 549–560 (1905)
work page 1905
-
[2]
Mouvement brownien et réalité moléculaire
Perrin, J. Mouvement brownien et réalité moléculaire. Ann. Chim. Phys. 18, 5–114 (1909)
work page 1909
-
[3]
Extended Structural Dynamics: Emergent irreversibility from reversible dynamics
BarAvi, P. Extended Structural Dynamics: Emergent irreversibility from reversible dynamics. arXiv:2505.09650 [physics.class-ph] (2025)
-
[4]
Structural viscosity, thermal waves, and the Mpemba effect from Extended Structural Dynamics
BarAvi, P. Structural viscosity, thermal waves, and the Mpemba effect from Extended Structural Dynamics. arXiv:2603.02249 [physics.class-ph] (2026a)
-
[5]
BarAvi, P. Extended Structural Dynamics and the Lorentz–Abraham–Dirac equation: A deformable charge interpretation. arXiv:2603.11064 [physics.class-ph] (2026b). [5a] BarAvi, P. (2026). Extended Structural Dynamics of Charged Media: Kinetic Theory, Generalized Plasma. Transport Phenomena, De Gruyter Brill, in press
-
[6]
Weiss, M., Elsner, M., Kartberg, F. & Nilsson, T. Anomalous subdiffusion is a measure for cytoplasmic crowding in living cells. Biophys. J. 87, 3518–3524 (2004)
work page 2004
-
[7]
Banks, D. S. & Fradin, C. Anomalous diffusion of proteins due to molecular crowding. Biophys. J. 89, 2960–2971 (2005)
work page 2005
- [8]
-
[9]
Jeon, J.-H., Monne, H. M.-S., Javanainen, M. & Metzler, R. Anomalous diffusion of phospholipids and cholesterols in a lipid bilayer and its origins. Phys. Rev. Lett. 109, 188103 (2012)
work page 2012
-
[10]
Jeon, J.-H. et al. In vivo anomalous diffusion and weak ergodicity breaking of lipid granules. Phys. Rev. Lett. 106, 048103 (2011)
work page 2011
-
[11]
Weber, S. C., Spakowitz, A. J. & Theriot, J. A. Bacterial chromosomal loci move subdiffusively through a viscoelastic cytoplasm. Phys. Rev. Lett. 104, 238102 (2010)
work page 2010
-
[12]
Crater, J. S. & Carrier, R. L. Barrier properties of gastrointestinal mucus to nanoparticle transport. Macromol. Biosci. 10, 1473–1483 (2010)
work page 2010
-
[13]
Manzo, C. et al. Weak ergodicity breaking of receptor motion in living cells stemming from random diffusivity. Phys. Rev. X 5, 011021 (2015)
work page 2015
-
[14]
Höfling, F. & Franosch, T. Anomalous transport in the crowded world of biological cells. Rep. Prog. Phys. 76, 046602 (2013). 34
work page 2013
-
[15]
Dix, J. A. & Verkman, A. S. Crowding effects on diffusion in solutions and cells. Annu. Rev. Biophys. 37, 247–263 (2008)
work page 2008
-
[16]
Metzler, R. & Klafter, J. The random walk's guide to anomalous diffusion: a fractional dynamics approach. Phys. Rep. 339, 1–77 (2000)
work page 2000
-
[17]
Fractional Calculus and Waves in Linear Viscoelasticity
Mainardi, F. Fractional Calculus and Waves in Linear Viscoelasticity. Imperial College Press (2010)
work page 2010
-
[18]
Montroll, E. W. & Weiss, G. H. Random walks on lattices. II. J. Math. Phys. 6, 167–181 (1965)
work page 1965
-
[19]
Weak ergodicity breaking and aging in disordered systems
Bouchaud, J.-P. Weak ergodicity breaking and aging in disordered systems. J. Phys. I France 2, 1705–1713 (1992)
work page 1992
-
[20]
The fluctuation-dissipation theorem
Kubo, R. The fluctuation-dissipation theorem. Rep. Prog. Phys. 29, 255–284 (1966)
work page 1966
-
[21]
Viscoelastic subdiffusion: Generalized Langevin equation approach
Goychuk, I. Viscoelastic subdiffusion: Generalized Langevin equation approach. Adv. Chem. Phys. 150, 187–253 (2012)
work page 2012
-
[22]
Laage, D., Elsaesser, T. & Hynes, J. T. Water dynamics in the hydration shells of biomolecules. Chem. Rev. 117, 10694–10725 (2017)
work page 2017
-
[23]
Water as an active constituent in cell biology
Ball, P. Water as an active constituent in cell biology. Chem. Rev. 108, 74–108 (2008)
work page 2008
- [24]
-
[25]
Frauenfelder, H., Sligar, S. G. & Wolynes, P. G. The energy landscapes and motions of proteins. Science 254, 1598–1603 (1991)
work page 1991
-
[26]
Teilum, K., Olsen, J. G. & Kragelund, B. B. Functional aspects of protein flexibility. Cell. Mol. Life Sci. 66, 2231–2247 (2009)
work page 2009
-
[27]
Nonequilibrium Statistical Mechanics
Zwanzig, R. Nonequilibrium Statistical Mechanics. Oxford University Press (2001)
work page 2001
-
[28]
Projection Operator Techniques in Nonequilibrium Statistical Mechanics
Grabert, H. Projection Operator Techniques in Nonequilibrium Statistical Mechanics. Springer (1982)
work page 1982
-
[29]
Cui, Q. & Bahar, I. Normal Mode Analysis: Theory and Applications to Biological and Chemical Systems. Chapman & Hall/CRC (2006)
work page 2006
-
[30]
Levitt, M., Sander, C. & Stern, P. S. Protein normal-mode dynamics: trypsin inhibitor, crambin, ribonuclease and lysozyme. J. Mol. Biol. 181, 423–447 (1985)
work page 1985
-
[31]
Tama, F. & Sanejouand, Y .-H. Conformational change of proteins arising from normal mode calculations. Protein Eng. 14, 1–6 (2001). 35
work page 2001
-
[32]
Carugo, O. & Argos, P. Correlation between side chain mobility and conformation in protein structures. Protein Eng. 10, 777–787 (1997)
work page 1997
-
[33]
Tilton, R. F., Dewan, J. C. & Petsko, G. A. Effects of temperature on protein structure and dynamics. Biochemistry 31, 2469–2481 (1992)
work page 1992
-
[34]
Palmer, A. G. NMR characterization of the dynamics of biomacromolecules. Chem. Rev. 104, 3623–3640 (2004)
work page 2004
-
[35]
Jarymowycz, V . A. & Stone, M. J. Fast time scale dynamics of protein backbones: NMR relaxation methods. Chem. Rev. 106, 1624–1671 (2006)
work page 2006
-
[36]
Happel, J. & Brenner, H. Low Reynolds Number Hydrodynamics. Martinus Nijhoff (1983)
work page 1983
-
[37]
Ohmachi, M. et al. Fluorescence microscopy for simultaneous observation of 3D orientation and movement and its application to quantum rod-tagged myosin V . Proc. Natl. Acad. Sci. USA 109, 5294–5298 (2012)
work page 2012
-
[38]
Rosenberg, S. A., Quinlan, M. E., Forkey, J. N. & Goldman, Y . E. Rotational motions of macromolecules by single-molecule fluorescence microscopy. Acc. Chem. Res. 38, 583–593 (2005)
work page 2005
-
[39]
Metzler, R., Jeon, J.-H., Cherstvy, A. G. & Barkai, E. Anomalous diffusion models and their properties: non-stationarity, non-ergodicity, and ageing at the centenary of single particle tracking. Phys. Chem. Chem. Phys. 16, 24128–24164 (2014)
work page 2014
-
[40]
Szymanski, J. & Weiss, M. Elucidating the origin of anomalous diffusion in crowded fluids. Phys. Rev. Lett. 103, 038102 (2009)
work page 2009
- [41]
-
[42]
Radmacher, M., Tillmann, R. W. & Gaub, H. E. Imaging viscoelasticity by force modulation with the atomic force microscope. Biophys. J. 64, 735–742 (1993)
work page 1993
-
[43]
Svergun, D. I. et al. Protein hydration in solution: experimental observation by x-ray and neutron scattering. Proc. Natl. Acad. Sci. USA 95, 2267–2272 (1998)
work page 1998
-
[44]
Levy, Y . & Onuchic, J. N. Water mediation in protein folding and molecular recognition. Annu. Rev. Biophys. Biomol. Struct. 35, 389–415 (2006)
work page 2006
-
[45]
Manton, J. D. et al. Concepts for high-speed multiplane confocal microscopy. J. Microsc. 281, 157–169 (2021)
work page 2021
-
[46]
Gwosch, K. C. et al. MINFLUX nanoscopy delivers 3D multicolor nanometer resolution in cells. Nat. Methods 17, 217–224 (2020). 36
work page 2020
-
[47]
Michaud-Agrawal, N., Denning, E. J., Woolf, T. B. & Beckstein, O. MDAnalysis: a toolkit for the analysis of molecular dynamics simulations. J. Comput. Chem. 32, 2319–2327 (2011)
work page 2011
-
[48]
Shaw, D. E. et al. Anton 2: raising the bar for performance and programmability in a special-purpose molecular dynamics supercomputer. Proc. IEEE Int. Conf. High Perform. Comput. Netw. Storage Anal. 41, 41–53 (2014)
work page 2014
-
[49]
Barkai, E., Garini, Y . & Metzler, R. Strange kinetics of single molecules in living cells. Phys. Today 65, 29–35 (2012)
work page 2012
-
[50]
Sabri, A., Xu, X., Krapf, D. & Weiss, M. Elucidating the origin of heterogeneous anomalous diffusion in the cytoplasm of mammalian cells. Nat. Commun. 13, 6717 (2022). Appendix A: Derivation of the Memory Kernel from Extended Structural Dynamics This appendix summarizes the relevant results from the ESD hydrodynamic framework [4]. The reader is referred t...
work page 2022
-
[51]
Scale separation: 𝜏int ≪ 𝜏obs, where 𝜏int ∼ 𝜔𝛼 −1 is the fastest internal relaxation time and 𝜏obs is the observation time
-
[52]
Weak coupling: 𝜆 ∣ ∇𝑉ext ∣≪ 𝑘𝐵𝑇/𝑎, where 𝑎 is the particle size
-
[53]
Linear response: Internal deformations remain small compared to particle dimensions For the physiological conditions considered in this paper (proteins in cytoplasm at 𝑇 = 300 K), these conditions are satisfied. 38 Appendix B: Mode Density and Power-Law Spectrum in Proteins This appendix provides background on the mode density of proteins and the power-la...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.