Pith. sign in

REVIEW 3 major objections 4 minor 5 cited by

In power-law viscoelastic media, active particles persist as superdiffusive t^(2-alpha_R) motion rather than ballistic t^2, and the persistence time is dramatically stretched, altering the standard link between propulsion and persistence.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · deepseek-v4-flash

2026-08-04 06:25 UTC pith:75EZZOLG

load-bearing objection Solid exact MSD for 2D fractional active particles, but the advertised t^{2−α_R} superdiffusive regime is a smooth crossover, not a power law. the 3 major comments →

arxiv 2512.20205 v1 pith:75EZZOLG submitted 2025-12-23 cond-mat.soft cond-mat.stat-mech

Active Brownian particles in power-law viscoelastic media

classification cond-mat.soft cond-mat.stat-mech MSC 82C3126A33
keywords active Brownian particlesviscoelasticityfractional Langevin equationpower-law memoryanomalous diffusionmean squared displacementpersistence timesuperdiffusion
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The paper studies an active Brownian particle in a viscoelastic medium whose memory decays as a power law, described by fractional Langevin equations with independent fractional orders for translation (alpha_T) and rotation (alpha_R). It derives an exact analytical expression for the mean squared displacement and shows that, when rotational memory is strong (alpha_R < 1), the persistent phase is no longer ballistic t^2 but superdiffusive t^(2-alpha_R). The persistence time and the crossover to normal diffusion are both stretched, with the effective diffusion time growing exponentially faster than the persistence time. This matters because many biological and soft materials, from gels to cytoplasm, exhibit power-law viscoelasticity, so the predictions are directly testable in microswimmer experiments.

Core claim

The central result is the analytic MSD for a fractional active Brownian particle: it combines a passive fractional term 2 D_T t^(alpha_T) with an active term expressed via incomplete gamma functions. In the intermediate persistent regime, the active contribution crosses over from ballistic t^2 to sub-ballistic superdiffusion t^(2-alpha_R), and the long-time diffusion coefficient becomes D_l = 2 v^2 Gamma(1+1/alpha_R)/(D_R/2)^(1/alpha_R). The competition between translational and rotational memory thus fundamentally modifies the standard active-matter relation between persistence and propulsion.

What carries the argument

The central object is the fractional Langevin equation with power-law memory kernel K_Q(t) ~ t^(-alpha_Q) for both translational and rotational degrees of freedom. The argument is carried by the Gaussian closure assumption: the angular displacement is taken to be Gaussian with variance <(phi-phi_0)^2> = D_R t^(alpha_R), giving the orientation autocorrelation <cos Delta phi> = exp(-D_R |t-t'|^(alpha_R)/2). This stretched-exponential memory, together with the power-law resolvent G_Q(t) = t^(alpha_Q-1)/(eta_Q Gamma(alpha_Q)), produces the incomplete-gamma-function structure of the MSD.

Load-bearing premise

The central MSD formula relies on the assumption that the angular displacement is Gaussian with zero mean, so that the orientation autocorrelation is exactly exp(-D_R |t-t'|^(alpha_R)/2); the paper itself notes that the rigorous Fokker-Planck construction for this non-Markovian process remains an open problem, and any non-Gaussianity would alter the persistence time and the superdiffusive phase.

What would settle it

Measure the MSD and the orientation autocorrelation of an active Brownian particle in a power-law viscoelastic medium (e.g., a colloidal gel or cytoplasm) across a range of alpha values. If the persistence phase exponent is not 2-alpha_R (with alpha_R from independent rheology), or if the orientation autocorrelation deviates from exp(-(t/tau)^alpha_R), the Gaussian closure and the predicted scaling would be falsified.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

Share X Bluesky LinkedIn Reddit HN

If this is right

  • If the central claim is correct, active particles in power-law viscoelastic media exhibit a persistent superdiffusive regime t^(2-alpha_R) instead of the standard ballistic t^2, which should be observable in microrheology experiments.
  • The persistence time tau_p = (2/D_R)^(1/alpha_R) Gamma(1+1/alpha_R) diverges as alpha_R -> 0, meaning orientation memory can become extremely long-lived in strongly viscoelastic environments.
  • The crossover to normal diffusion is delayed, with tau_l growing approximately as (4/alpha_R)^(1/alpha_R), so the enhanced-diffusion plateau may be pushed to timescales far beyond the persistence time.
  • When alpha_T and alpha_R are decoupled, the short-time subdiffusive regime and the persistent superdiffusive regime are governed by independent timescales, allowing a wider range of observable diffusion spectra.
  • The derived analytic MSD and numerical L1 scheme provide a toolkit for studying collective active behavior in viscoelastic media.

Where Pith is reading between the lines

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

  • The Gaussian closure is likely the main limitation: if the angular displacement is non-Gaussian, the orientation autocorrelation and the persistence time would change, so the predicted t^(2-alpha_R) exponent may be sensitive to that assumption; a rigorous non-Markovian Fokker-Planck treatment would settle this.
  • A direct experimental test would be to measure the MSD of a Brownian microswimmer in a well-characterized power-law gel and check whether the persistence phase scales as t^(2-alpha_R) with alpha_R extracted from independent rheology; the same alpha_R should also appear in the orientation autocorrelation.
  • The exponential divergence of tau_l relative to tau_p suggests that in many strongly viscoelastic biological environments, the apparent superdiffusion could persist over the entire experimental time window, meaning conventional 'enhanced diffusion' coefficients may be misleadingly large.
  • The framework could be extended to active particles with alignment interactions or to confined geometries, where the fractional memory may produce collective anomalous phases not captured by Markovian models.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 4 minor

Summary. The paper extends overdamped active Brownian motion to power-law viscoelastic media by replacing the ordinary friction with fractional Caputo derivatives of orders α_T and α_R and coupling the translational and rotational coordinates to independent fractional Gaussian noises whose amplitudes are fixed by fluctuation–dissipation relations. The authors derive the angular MSD, the orientation autocorrelation Eq. (11), and the total MSD Eq. (12), expressed as a passive subdiffusive term plus an active term involving incomplete gamma functions. They also construct an L1 discretization of the fractional Langevin equations and compare the numerics with the analytic formulas. The central interpretive claim is that the active MSD exhibits a distinct intermediate superdiffusive regime t^{2−α_R}, with an over-stretched persistence time relative to ordinary active Brownian motion.

Significance. If the claimed t^{2−α_R} regime existed as an extended power law, the paper would be a useful and nontrivial addition to the active-matter literature. The exact solution Eq. (12) and the explicit L1 numerical scheme are valuable and appear to be derived correctly; the model has no ad hoc free parameters beyond the fractional orders and standard transport coefficients. However, the headline regime claim is not supported by the exact solution, and the text's caveat about the Gaussian approximation is unnecessary because the angular process is exactly Gaussian under the stated linear fractional-Gaussian-noise model. The paper's lasting contribution is therefore the closed-form MSD and the observation that viscoelastic memory broadens the ballistic-to-diffusive crossover, rather than the existence of a new t^{2−α_R} asymptotic regime.

major comments (3)
  1. [Active diffusion section, Eq. (12), Appendix C] The claimed intermediate superdiffusive regime t^{2−α_R} is not contained in the exact solution. From the series in Appendix C, the active contribution is 2v^2 t^2 Σ_{k=0}∞ (-D_R t^{α_R}/2)^k / [k!(kα_R+1)(kα_R+2)]. The leading short-time term is v^2 t^2 and the long-time limit is linear in t; there is no term proportional to t^{2−α_R}. The local log-log slope is a smooth function that decreases monotonically from 2 to 1 and passes through the value 2−α_R at a single instant, not over a decade of time. For example, for α_R=0.3 the slope at z=D_R t^{α_R}/2 = 1, 2, 3 is approximately 1.81, 1.64, and 1.51, respectively, with no plateau near 1.7. Thus the statements in the text and Fig. 2 about a 'second transition from ballistic t^2 to sub-ballistic superdiffusive t^{2−α_R}' and a 'zoom into the t^{2−α_R} superdiffusive regime' are overinterpretations of a featureless crossover. The manuscr
  2. [Equation (11) and the paragraph introducing it] The text presents Eq. (11) as a Gaussian approximation and states that the corresponding Fokker–Planck construction is an open problem and generally non-Gaussian. For the model as defined, this is not an approximation. Equation (7) gives φ(t)−φ_0 = θ_R (ξ_R ∗ G_R)(t), where G_R(t)=t^{α_R−1}/(η_R Γ(α_R)) is a deterministic kernel and ξ_R is fractional Gaussian noise. A linear functional of a Gaussian process is Gaussian, so Δφ is exactly Gaussian with variance D_R t^{α_R}. Consequently, ⟨cos Δφ⟩ = exp(−D_R |Δt|^{α_R}/2) is exact within the model. The caveat should be removed or rephrased to refer to possible non-Gaussian generalizations of the noise; as written, the text unnecessarily weakens its own result and is internally inconsistent with the statement that the solution is a linear combination of Gaussian processes.
  3. [Numerical validation] The numerical comparison is presented as validation of the theoretical predictions, but it actually checks the L1 discretization against the analytic solution of the same equations. That is a legitimate check of the numerical scheme, but it is not an independent test of the physical model. In addition, the thermalization analysis in Fig. 1(a) shows that for small α_R the convergence to the overdamped limit occurs on a time scale O(10 t/τ_R), so the early-time numerical points for α_R=0.2–0.3 may not yet represent the stationary overdamped process. The manuscript should state this limitation explicitly when discussing the agreement between analytics and numerics.
minor comments (4)
  1. [Eq. (1) and Eq. (12)] The independence of the translational and rotational noises ξ_T and ξ_R is assumed implicitly but never stated. Since the total MSD in Eq. (12) is written as the sum of passive and active contributions, any cross-correlation between the two noises would add a term. Please state the independence assumption explicitly.
  2. [Appendix B] In Eq. (B2) the noise-amplitude replacement uses H_R for the general coordinate q; this should be H_Q or the corresponding subscript should be defined consistently. Also, the notation ζ_Q in the text after Eq. (2) is not used later.
  3. [Throughout] Several typos and misspellings should be corrected: 'directely' (p. 3), 'propuslion', 'persistance', and 'over-streched' (Conclusion). The phrase 'over-streched' should also be replaced by 'over-stretched'.
  4. [Fig. 1 and Fig. 2] The captions and text refer to 'analytical expressions' and 'numerical results' without always stating the parameter values used for D_T, D_R, η_Q, and the activity. Please provide the full parameter set in the figure captions or in the main text so that the comparisons are reproducible.

Circularity Check

0 steps flagged

No significant circularity: the MSD derivation follows from the stated fractional Langevin model, and the self-citations are peripheral.

full rationale

The central derivation is self-contained. The model inputs are the power-law memory kernels and fractional Gaussian noises in Eqs. (1)-(5), together with the fluctuation-dissipation relations in Eq. (6). The angular MSD is derived in Appendix B, the orientation autocorrelation follows from the Gaussian property of the linear rotational equation in Eq. (11), and the active MSD in Eq. (12) is obtained by direct integration in Appendix C. No parameter is fitted to the MSD to produce Eq. (12); the numerical solution in Eqs. (9)-(10) solves the same equations, so agreement is a consistency check rather than an independent empirical test, but it is not a fit masquerading as a prediction. The paper explicitly flags the Fokker-Planck construction for non-Markovian processes as an open problem, but this is a stated limitation of the formalism, not a circular step: the result is conditional on the model and the Gaussian property, not secretly imposed by it. The few self-citations ([43], [44], [57]) concern prethermal numerical effects and the open Fokker-Planck problem; they are not load-bearing for the MSD formula and no uniqueness theorem is imported. The skeptic's concern that the claimed t^{2-alpha_R} superdiffusive regime is not a true plateau in the series Eq. (C2) is a correctness/interpretation issue about the derived function, not a circularity, because the claim is an asserted reading of the derived series rather than an input to its derivation. Overall, no step reduces by construction to its own inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 5 axioms · 0 invented entities

The ledger is tidy: the paper introduces no new particles or forces. Its physics is a combination of established fractional Langevin machinery. The main new assumptions are the Gaussian closure for the orientation and the overdamped stationary-state limit. The H_R >= 1/2 constraint from the FDR is a modeling assumption rather than a validated prediction.

axioms (5)
  • domain assumption The system is overdamped and already at the stationary state at t = 0, so inertial terms are absent from Eq. (1).
    Necessary for the resolvent G_Q = t^{alpha-1}/eta Gamma(alpha) and the MSD D_Q t^alpha. The paper notes the stationary-state assumption in App. B.
  • ad hoc to paper The orientation is a Gaussian process with zero mean and variance sigma^2 = <(phi - phi0)^2> + |phi0|^2, giving Eq. (11).
    The text itself states this is a Gaussian approximation because the Fokker-Planck equation for the non-Markovian process is an open problem. This is the load-bearing approximation for the whole active MSD.
  • domain assumption The overdamped Caputo fractional Langevin equation (4) correctly describes the viscoelastic medium's response.
    The claim that the medium is a power-law relaxation material is imported from the GSE literature; the paper does not critique the GSE. The active force is inserted into the constitutive equation without a hydrodynamic derivation, so the level of modeling is asserted.
  • domain assumption Fluctuation-dissipation relation of the second kind, Eq. (6), linking noise correlation and memory kernel, holds for both translational and rotational components.
    This yields alpha_Q = 2 - 2H_Q and confines the Hurst exponents to [1/2,1), and the comparison with Ref. [59] at the end of the persistence section explicitly depends on it.
  • domain assumption The noise correlations for the angular and translational degrees of freedom are independent of each other, and mode-interference terms are neglected in the MSD.
    The MSD separates into translational and active contributions in Eq. (12); the cross-noise term is implicitly dropped. Whether the angular noise couples to the self-propulsion direction is an unstated modeling simplification.

pith-pipeline@v1.3.0-alltime-deepseek · 12110 in / 6192 out tokens · 51682 ms · 2026-08-04T06:25:52.496699+00:00 · methodology

0 comments
Cite this review

Pith. "Pith review of Active Brownian particles in power-law viscoelastic media." pith.science (2026). https://pith.science/paper/75EZZOLG

@misc{pith2026251220205,
  author       = {Pith},
  title        = {Pith review of: Active Brownian particles in power-law viscoelastic media},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/75EZZOLG}},
  note         = {Machine review of arXiv:2512.20205}
}
Share X Bluesky LinkedIn Reddit HN
read the original abstract

Many active particles are embedded in environments that exhibit viscoelastic properties. An important class of such media lacks a single characteristic relaxation timescale when subjected to a time-dependent stress. Rather, the stress response spans a broad continuum of timescales, a behavior naturally described by a scale-free, fractal-like power-law relaxation modulus. Using a generalization of the fractional Langevin equation, we investigate an active Brownian particle embedded in a power-law viscoelastic environment with translational and rotational dynamics governed by independent fractional orders. We solve the model analytically, develop a numerical scheme to validate the theoretical predictions, and provide tools that can be used in further studies. A rich variety of diffusion regimes emerges, which modify the intermediate-time behavior of the mean squared displacement. Notably, we find that the competition between translational and rotational contributions favors a superdiffusive persistence over the standard ballistic motion, and over-stretches its characteristic timescale, fundamentally altering the standard relation between persistence and propulsion in active matter.

Figures

Figures reproduced from arXiv: 2512.20205 by Cristiane Morais Smith, David Santiago Quevedo, Marjolein Dijkstra, Monica Conte.

Figure 1
Figure 1. Figure 1: FIG. 1. Dynamics of the angular component of an active [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Active diffusion regimes for [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Diffusion regimes for decoupled translational and ro [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗

discussion (0)

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

Forward citations

Cited by 5 Pith papers

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

  1. Atom-Photon Bound States in Fractal Photonic Lattices: Localization Length and Anomalous Diffusion

    quant-ph 2026-05 unverdicted novelty 7.0

    Atom-photon bound states in fractal photonic lattices exhibit localization length ξ ∼ Δ^{-1/d_w} governed by anomalous diffusion on the fractal.

  2. A particle-resolved rheological study of chirality transfer and odd transport

    cond-mat.stat-mech 2026-05 unverdicted novelty 6.0

    Collisions transfer chirality from an active bath to a passive tracer, yielding circular trajectories and odd transverse drift under force, rectified by nonlinear friction.

  3. Confined kinetics and heterogeneous diffusion driven by fractional Gaussian noise: A path integral approach

    cond-mat.stat-mech 2026-04 unverdicted novelty 6.0

    Path-integral analysis of multiplicative fractional Gaussian noise yields a Gaussian propagator in terms of the Lamperti transform and demonstrates that confinement induces an effective drift causing probability accum...

  4. Confined kinetics and heterogeneous diffusion driven by fractional Gaussian noise: A path integral approach

    cond-mat.stat-mech 2026-04 unverdicted novelty 6.0

    Path integral analysis of multiplicative fractional Gaussian noise yields a Gaussian propagator via Lamperti transform and reveals an effective drift causing probability buildup in low-noise regions under confinement.

  5. Passive memory reshapes active persistence

    cond-mat.soft 2026-05 unverdicted novelty 5.0

    Memory in the surrounding medium suppresses motility-induced phase separation in active particles when memory and persistence timescales match, but restores clustering for long memory via reduced short-time friction.

Reference graph

Works this paper leans on

61 extracted references · 3 linked inside Pith · cited by 4 Pith papers

  1. [1]

    te Vrugt, B

    M. te Vrugt, B. Liebchen, and M. E. Cates, What exactly is ’active matter’ ? (2025), arXiv:2507.21621

  2. [2]

    Bechinger, R

    C. Bechinger, R. Di Leonardo, H. L¨ owen, C. Reichhardt, G. Volpe, and G. Volpe, Active particles in complex and crowded environments, Rev. Mod. Phys.88, 045006 (2016)

  3. [3]

    Ramaswamy, The mechanics and statistics of active matter, Annu

    S. Ramaswamy, The mechanics and statistics of active matter, Annu. Rev. Condens. Matter Phys.1, 323 (2010)

  4. [4]

    Vicsek, A

    T. Vicsek, A. Czir´ ok, E. Ben-Jacob, I. Cohen, and O. Shochet, Novel type of phase transition in a system of self-driven particles, Phys. Rev. Lett.75, 1226 (1995)

  5. [5]

    Chat´ e, F

    H. Chat´ e, F. Ginelli, G. Gr´ egoire, F. Peruani, and F. Ray- naud, Modeling collective motion: variations on the vic- sek model, Eur. Phys. J. B64, 451 (2008)

  6. [6]

    Liu and M

    Z. Liu and M. Dijkstra, Collective dynamics of intelligent active Brownian particles with visual perception and ve- locity alignment in 3d: spheres, rods, and worms, Soft Matter21, 1529 (2025)

  7. [7]

    Budrikis, Collective motion strategies of sheep, Nat

    Z. Budrikis, Collective motion strategies of sheep, Nat. Rev. Phys.5, 82 (2023)

  8. [8]

    Turner, R

    L. Turner, R. Zhang, N. C. Darnton, and H. C. Berg, Vi- sualization of flagella during bacterial swarming, J. Bac- teriol.192, 3259 (2010)

  9. [9]

    Elgeti, R

    J. Elgeti, R. G. Winkler, and G. Gompper, Physics of mi- croswimmers—single particle motion and collective be- havior: a review, Rep. Prog. Phys.78, 056601 (2015)

  10. [10]

    Thakur and R

    S. Thakur and R. Kapral, Dynamics of self-propelled nanomotors in chemically active media, J. Chem. Phys. 135, 024509 (2011)

  11. [11]

    L¨ owen, Gigantic dynamical spreading and anoma- lous diffusion of jerky active particles (2025), arXiv:2507.08910

    H. L¨ owen, Gigantic dynamical spreading and anoma- lous diffusion of jerky active particles (2025), arXiv:2507.08910

  12. [12]

    Zhong, Anomalous collective diffusion in two- dimensional active particle systems, Phys

    W. Zhong, Anomalous collective diffusion in two- dimensional active particle systems, Phys. Rev. E112, 054112 (2025)

  13. [13]

    Stenhammar, D

    J. Stenhammar, D. Marenduzzo, R. J. Allen, and M. E. Cates, Phase behaviour of active Brownian particles: the role of dimensionality, Soft Matter10, 1489 (2014)

  14. [14]

    Surrey, F

    T. Surrey, F. N´ ed´ elec, S. Leibler, and E. Karsenti, Phys- ical properties determining self-organization of motors and microtubules, Science292, 1167 (2001)

  15. [15]

    Lakes,Viscoelastic Materials(Cambridge Univ

    R. Lakes,Viscoelastic Materials(Cambridge Univ. Press, Cambridge, 2009)

  16. [16]

    T. G. Mason, Estimating the viscoelastic moduli of com- plex fluids using the generalized Stokes–Einstein equa- tion, Rheol. Acta39, 371 (2000)

  17. [17]

    T. G. Mason, H. Gang, and D. A. Weitz, Diffusing-wave- spectroscopy measurements of viscoelasticity of complex fluids, J. Opt. Soc. Am. A14, 139 (1997)

  18. [18]

    Narinder, C

    N. Narinder, C. Bechinger, and J. R. Gomez-Solano, Memory-induced transition from a persistent random walk to circular motion for achiral microswimmers, Phys. Rev. Lett.121, 078003 (2018)

  19. [19]

    A. R. Sprenger, C. Bair, and H. L¨ owen, Active Brown- ian motion with memory delay induced by a viscoelastic medium, Phys. Rev. E105, 044610 (2022)

  20. [20]

    Bonfanti, J

    A. Bonfanti, J. L. Kaplan, G. Charras, and A. Kabla, Fractional viscoelastic models for power-law materials, Soft Matter16, 6002 (2020)

  21. [21]

    T. G. Mason and D. A. Weitz, Optical measurements of frequency-dependent linear viscoelastic moduli of com- plex fluids, Phys. Rev. Lett.74, 1250 (1995)

  22. [22]

    Sinkus, K

    R. Sinkus, K. Siegmann, T. Xydeas, M. Tanter, C. Claussen, and M. Fink, MR elastography of breast lesions: Understanding the solid/liquid duality can im- prove the specificity of contrast-enhanced mr mammog- raphy, Magn. Reson. Med.58, 1135 (2007)

  23. [23]

    S. Aney, P. Pandit, L. Ratke, and et al., On the origin of power-scaling exponents in silica aerogels, J. Sol-Gel Sci. Technol.114, 98 (2025)

  24. [24]

    S. Aime, L. Cipelletti, and L. Ramos, Power law vis- coelasticity of a fractal colloidal gel, J. Rheol.62, 1429 (2018)

  25. [25]

    T. S. K. Ng and G. H. McKinley, Power law gels at finite strains: The nonlinear rheology of gluten gels, J. Rheol. 52, 417 (2008)

  26. [26]

    Gittes, B

    F. Gittes, B. Schnurr, P. D. Olmsted, F. C. MacKintosh, and C. F. Schmidt, Microscopic viscoelasticity: Shear moduli of soft materials determined from thermal fluctu- ations, Phys. Rev. Lett.79, 3286 (1997)

  27. [27]

    Dichtl and E

    M. Dichtl and E. Sackmann, Colloidal probe study of short time local and long time reptational motion of semi- flexible macromolecules in entangled networks, New J. Phys.1, 18 (1999)

  28. [28]

    Caspi, R

    A. Caspi, R. Granek, and M. Elbaum, Enhanced diffu- sion in active intracellular transport, Phys. Rev. Lett.85, 5655 (2000)

  29. [29]

    Flenner, J

    E. Flenner, J. Das, M. C. Rheinst¨ adter, and I. Kosztin, Subdiffusion and lateral diffusion coefficient of lipid atoms and molecules in phospholipid bilayers, Phys. Rev. E79, 011907 (2009)

  30. [30]

    S. Joo, X. Durang, O.-c. Lee, and J.-H. Jeon, Anomalous diffusion of active Brownian particles cross-linked to a polymer network, Soft Matter16, 9188 (2020)

  31. [31]

    S. Mani, H. H. Winter, M. Silverstein, and et al., Power law relaxation in an interpenetrating polymer network, Colloid Polym. Sci.267, 1002 (1989)

  32. [32]

    Zhang, E

    Z. Zhang, E. Bouchbinder, and E. Lerner, The strain- stiffening critical exponents in polymer networks and their universality, J. Chem. Phys.163, 111102 (2025)

  33. [33]

    Maiti, Y

    A. Maiti, Y. Koyano, H. Kitahata, and K. K. Dey, Activity-induced diffusion recovery in crowded colloidal 7 suspensions, Phys. Rev. E109, 054607 (2024)

  34. [34]

    Ebata and K

    H. Ebata and K. e. a. Nishizawa, Single power-law rheology of crowded cytoplasm in living cells (2025), arXiv:2504.18922

  35. [35]

    Joo and J.-H

    S. Joo and J.-H. Jeon, Viscoelastic active diffusion gov- erned by nonequilibrium fractional Langevin equations: Underdamped dynamics and ergodicity breaking, Chaos, Solitons Fractals177, 114288 (2023)

  36. [36]

    Y. Kang, S. Seo, S. Kwon, and K. Kim, Fractional motion of an active particle in fractional generalized Langevin equations, Fractal Fract.9, 725 (2025)

  37. [37]

    Metzler and T

    R. Metzler and T. F. Nonnenmacher, Fractional relax- ation processes and fractional rheological models, Int. J. Plast.19, 941 (2003)

  38. [38]

    Hilfer,Applications of Fractional Calculus in Physics (World Scientific, Singapore, 2000)

    R. Hilfer,Applications of Fractional Calculus in Physics (World Scientific, Singapore, 2000)

  39. [39]

    Klatt, U

    D. Klatt, U. Hamhaber, P. Asbach, J. Braun, and I. Sack, Noninvasive assessment of the rheological behavior of human organs using multifrequency mr elastography: a study of brain and liver viscoelasticity, Phys. Med. Biol. 52, 7281 (2007)

  40. [40]

    F. C. Meral, T. J. Royston, and R. Magin, Fractional calculus in viscoelasticity: An experimental study, Com- mun. Nonlinear Sci. Numer. Simul.15, 939 (2010)

  41. [41]

    Metzler and J

    R. Metzler and J. Klafter, The random walk’s guide to anomalous diffusion: a fractional dynamics approach, Phys. Rep.339, 1 (2000)

  42. [42]

    Lutz, Fractional Langevin equation, Phys

    E. Lutz, Fractional Langevin equation, Phys. Rev. E64, 051106 (2001)

  43. [43]

    D. S. Quevedo, R. C. Verstraten, and C. Morais Smith, Emergent transient time crystal from a fractional Langevin equation with white and colored noise, Phys. Rev. A110, 052208 (2024)

  44. [44]

    R. C. Verstraten, R. F. Ozela, and C. M. Smith, Time glass: A fractional calculus approach, Phys. Rev. B103, L180301 (2021)

  45. [45]

    Caputo, Linear models of dissipation whose Q is al- most frequency independent—II, Geophys

    M. Caputo, Linear models of dissipation whose Q is al- most frequency independent—II, Geophys. J. Int13, 529 (1967)

  46. [46]

    Risken,The Fokker-Planck equation: methods of solu- tion and applications, Vol

    H. Risken,The Fokker-Planck equation: methods of solu- tion and applications, Vol. 2 (Springer. Heidelberg, Ger- many., 1989)

  47. [47]

    S. F. Kwok,Langevin and Fokker-Planck Equations and Their Generalizations: Descriptions and Solutions (World Scientific. New Jersey, USA, 2018)

  48. [48]

    P. Guo, C. Zeng, C. Li, and Y. Chen, Numerics for the fractional Langevin equation driven by the frac- tional Brownian motion, Fract. Calc. Appl. Anal.16, 123 (2013)

  49. [49]

    Murillo-Arcila, A

    M. Murillo-Arcila, A. Peris, and ´A. Vargas-Moreno, Dy- namics of the caputo fractional derivative, Fract. Calc. Appl. Anal. , 1 (2025)

  50. [50]

    Li and F

    C. Li and F. Zeng,Numerical methods for fractional cal- culus, Vol. 24 (CRC Press, 2015)

  51. [51]

    Flynn, fbm 0.3.0: Fractional Brownian motion real- izations

    C. Flynn, fbm 0.3.0: Fractional Brownian motion real- izations. (2019)

  52. [52]

    Wang, Long-time-correlation effects and biased anomalous diffusion, Phys

    K. Wang, Long-time-correlation effects and biased anomalous diffusion, Phys. Rev. A45, 833 (1992)

  53. [53]

    J. M. Porra, K.-G. Wang, and J. Masoliver, Generalized Langevin equations: Anomalous diffusion and probabil- ity distributions, Phys. Rev. E53, 5872 (1996)

  54. [54]

    Calvo and R

    I. Calvo and R. Sanchez, The path integral formulation of fractional Brownian motion for the general Hurst ex- ponent, J. Phys. A: Math. Theor.41, 282002 (2008)

  55. [55]

    H¨ anggi, Correlation functions and masterequations of generalized (non-markovian) langevin equations, Z

    P. H¨ anggi, Correlation functions and masterequations of generalized (non-markovian) langevin equations, Z. Phys. B: Condens. Matter31, 407 (1978)

  56. [56]

    Khan and A

    S. Khan and A. M. Reynolds, Derivation of a Fokker– Planck equation for generalized Langevin dynamics, Physica A350, 183 (2005)

  57. [57]

    D. S. Quevedo, F. S. Abril-Berm´ udez, and C. Morais Smith, Path integral and generalized Klein-Kramers equation for Langevin dynamics with fractional Gaussian noise, arXiv preprint (2026)

  58. [58]

    P. K. Ghosh, Y. Li, G. Marchegiani, and F. Marchesoni, Communication: Memory effects and active Brownian diffusion, J. Chem. Phys.143(2015)

  59. [59]

    J. R. Gomez-Solano and F. J. Sevilla, Active particles with fractional rotational Brownian motion, Jour. Stat. Mech: Theory Exp.2020, 063213 (2020)

  60. [60]

    Volpe, S

    G. Volpe, S. Gigan, and G. Volpe, Simulation of the ac- tive Brownian motion of a microswimmer, Am. J. Phys. 82, 659 (2014)

  61. [61]

    Mainardi and R

    F. Mainardi and R. Gorenflo, On Mittag-Leffler-type functions in fractional evolution processes, J. Comput. Appl. Math.118, 283 (2000). End Matter Appendix A. L1 discretization of the fractional active Langevin equations—On a uniform time gridt j =k∗h withhthe step size andk∈[0, . . . , n], the integral over the time interval in Eq. (3) can be split into ...