pith. sign in

arxiv: 2602.04545 · v2 · submitted 2026-02-04 · ⚛️ physics.plasm-ph · cond-mat.stat-mech

Nonlinear Dynamical Friction from the Doppler-Shifted Equilibrium Memory Kernel

N. R. Sree Harsha , Zhenyuan Yu , Chuang Ren , Virginia Billings , Michael Huang This is my paper

Pith reviewed 2026-05-16 07:22 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph cond-mat.stat-mech
keywords dynamical frictionmemory kernelgeneralized Langevin equationnon-equilibrium steady stateDoppler shiftfluctuation-dissipation theoremplasma dragparticle-in-cell simulation
0
0 comments X

The pith

An equilibrium memory kernel, when Doppler-shifted for particle velocity, fully describes nonlinear dynamical friction in non-equilibrium plasma states.

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

The paper shows that a memory kernel extracted from equilibrium force fluctuations via the fluctuation-dissipation theorem is enough to model test-particle drag once a Doppler shift is applied for the particle's motion relative to the plasma. This kernel is inserted into the generalized Langevin equation, which then reproduces the full non-equilibrium steady-state dynamics without extra driving terms. A reader would care because the approach replaces expensive non-equilibrium simulations with cheaper equilibrium runs plus an analytic shift, while automatically recovering non-Markovian effects such as mass renormalization and oscillatory relaxation. In the Markovian limit the same kernel reduces exactly to the classic Chandrasekhar stopping-power formula. Direct particle-in-cell simulations confirm the predicted kernel shape and the resulting drag force.

Core claim

The equilibrium-derived memory kernel, obtained from the stochastic force autocorrelation in a thermal state and Doppler-shifted by the test-particle velocity, is sufficient inside the generalized Langevin equation to generate the correct nonlinear friction force and all associated non-Markovian phenomena in a non-equilibrium steady state.

What carries the argument

The Doppler-shifted equilibrium memory kernel: the equilibrium force autocorrelation function shifted in time by the relative velocity between test particle and plasma background.

If this is right

  • The Chandrasekhar stopping-power formula emerges directly as the Markovian limit of the same kernel.
  • Non-Markovian signatures such as effective mass renormalization and oscillatory relaxation appear automatically from the kernel moments.
  • First-principles equilibrium simulations alone yield quantitative predictions for non-equilibrium friction coefficients.
  • The framework supplies a computationally cheaper route to complex plasma drag problems than full non-equilibrium runs.

Where Pith is reading between the lines

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

  • The same Doppler-shift construction may apply to other velocity-dependent transport problems where equilibrium fluctuations are easier to sample than the driven state.
  • If the kernel shape is known analytically or from short equilibrium runs, the method could reduce the cost of repeated drag calculations across a range of test-particle speeds.
  • The approach highlights that the dominant non-equilibrium correction for a moving test particle is kinematic (the Doppler shift) rather than a change in the underlying fluctuation spectrum.

Load-bearing premise

The equilibrium memory kernel, once Doppler-shifted, already contains every correction needed for the non-equilibrium drag without requiring state-dependent or driving-specific additions.

What would settle it

A particle-in-cell run at a chosen velocity in which the measured drag force or the extracted memory-kernel oscillations deviate measurably from the values computed from the Doppler-shifted equilibrium kernel.

Figures

Figures reproduced from arXiv: 2602.04545 by Chuang Ren, Michael Huang, N. R. Sree Harsha, Virginia Billings, Zhenyuan Yu.

Figure 1
Figure 1. Figure 1: Comparison of the test particle velocity decay for [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Comparison of the test particle velocity decay for [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
read the original abstract

We present a statistical mechanics framework for modeling equilibrium friction coefficients using the Generalized Langevin Equation (GLE). We show that the kernel, obtained via the Fluctuation-Dissipation Theorem (FDT) from the stochastic force autocorrelation measured in a thermal equilibrium state, is sufficient to model the dynamics of the system in a Non-Equilibrium Steady State (NESS). This approach provides a computationally efficient path to modeling complex equilibrium friction problems. We apply this framework to the canonical problem of test particle drag in a uniform plasma. The GLE formalism is shown to naturally capture non-Markovian phenomena through the moments of the kernel, including an effective mass renormalization and oscillatory relaxation. We demonstrate that the standard Chandrasekhar stopping power formula arises naturally as the Markovian limit of this equilibrium memory kernel. These theoretical predictions are quantitatively validated by direct Particle-in-Cell simulations, which confirm the predicted oscillatory structure of the memory kernel. This work thus establishes a practical method for predicting equilibrium friction properties from first-principles equilibrium simulations.

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 claims that a memory kernel derived via the Fluctuation-Dissipation Theorem from the stochastic force autocorrelation in thermal equilibrium, when Doppler-shifted for test-particle velocity, is sufficient to model nonlinear dynamical friction in a non-equilibrium steady state (NESS) using the Generalized Langevin Equation. Applied to test-particle drag in uniform plasma, the framework captures non-Markovian effects including effective mass renormalization and oscillatory relaxation; the Chandrasekhar stopping power emerges as the Markovian limit, and the oscillatory kernel structure is validated by direct Particle-in-Cell simulations.

Significance. If the central claim holds, the work supplies a computationally efficient route to equilibrium friction coefficients from first-principles equilibrium data alone, extending to NESS dynamics without full non-equilibrium runs. It recovers a known Markov limit and supplies an independent PIC check of kernel structure, offering a practical bridge between equilibrium statistical mechanics and plasma drag problems.

major comments (2)
  1. [Abstract and derivation of NESS dynamics] Abstract and derivation of NESS dynamics: the central claim that the equilibrium-derived, Doppler-shifted kernel K(t) fully determines NESS friction without additional state-dependent or driving-induced terms is not demonstrated by direct comparison of the predicted non-Markovian drag trajectory (including effective mass renormalization) against driven simulations; the rigid frequency shift may miss v-dependent wake and resonant-particle alterations to the spectrum.
  2. [PIC validation section] PIC validation section: the abstract and validation paragraphs report confirmation of oscillatory kernel structure but supply no quantitative error bars, exclusion criteria, or direct comparison of the full drag force time series, leaving the quantitative accuracy of the NESS prediction unassessed.
minor comments (2)
  1. [Abstract] Abstract lacks any numerical measure of agreement between the GLE prediction and PIC data.
  2. [Introduction/Methods] Notation for the memory kernel and its Doppler shift should be defined explicitly with an equation number at first use.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for providing constructive feedback. We address each of the major comments below and outline the revisions we plan to make.

read point-by-point responses

  1. Referee: [Abstract and derivation of NESS dynamics] Abstract and derivation of NESS dynamics: the central claim that the equilibrium-derived, Doppler-shifted kernel K(t) fully determines NESS friction without additional state-dependent or driving-induced terms is not demonstrated by direct comparison of the predicted non-Markovian drag trajectory (including effective mass renormalization) against driven simulations; the rigid frequency shift may miss v-dependent wake and resonant-particle alterations to the spectrum.

    Authors: The derivation presented in the manuscript demonstrates that the Doppler-shifted equilibrium kernel, when inserted into the GLE, yields the NESS dynamics without requiring additional terms, as the memory kernel encodes the necessary correlations. However, we agree that an explicit numerical validation comparing the full predicted drag trajectory (including mass renormalization effects) to driven simulations would strengthen the presentation. We will include such a comparison in the revised manuscript. Concerning the rigid frequency shift, this approximation follows from transforming the equilibrium force autocorrelation to the test-particle frame; v-dependent alterations to the wake are partially captured by the velocity-dependent shift, though we acknowledge that a fully self-consistent non-equilibrium spectrum would require additional modeling beyond the current framework. revision: yes

  2. Referee: [PIC validation section] PIC validation section: the abstract and validation paragraphs report confirmation of oscillatory kernel structure but supply no quantitative error bars, exclusion criteria, or direct comparison of the full drag force time series, leaving the quantitative accuracy of the NESS prediction unassessed.

    Authors: We thank the referee for pointing this out. The current validation focuses on the oscillatory structure of the kernel extracted from equilibrium simulations. In the revised manuscript, we will add quantitative error bars to the kernel comparisons, specify the criteria used for identifying the oscillatory features, and include a direct comparison of the GLE-predicted drag force time series against the driven PIC simulations to assess the quantitative accuracy of the NESS predictions. revision: yes

Circularity Check

0 steps flagged

Equilibrium FDT kernel applied to NESS via Doppler shift; independent PIC validation prevents circularity

full rationale

The derivation extracts the memory kernel K(t) directly from the equilibrium force autocorrelation via the standard Fluctuation-Dissipation Theorem, then applies a rigid Doppler shift for test-particle velocity to obtain the NESS friction. The Markovian limit is shown to recover the known Chandrasekhar stopping power by direct reduction of the GLE moments, and the oscillatory kernel structure is confirmed by separate Particle-in-Cell simulations that are not used to fit or tune the target drag result. No load-bearing self-citation, self-definitional step, or fitted-input-renamed-as-prediction is present; the central claim that the equilibrium kernel suffices for NESS is tested against an external benchmark rather than being true by construction. This yields only a minor self-citation risk at most and keeps the overall circularity low.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The framework rests on standard statistical mechanics results without introducing new free parameters or postulated entities.

axioms (2)
  • domain assumption Fluctuation-Dissipation Theorem applies to the stochastic force autocorrelation in thermal equilibrium
    Invoked to obtain the memory kernel from equilibrium measurements.
  • domain assumption Generalized Langevin Equation governs the test-particle dynamics in both equilibrium and NESS
    Basis for extending the equilibrium kernel to non-equilibrium drag.

pith-pipeline@v0.9.0 · 5488 in / 1268 out tokens · 28320 ms · 2026-05-16T07:22:16.557694+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

34 extracted references · 34 canonical work pages

  1. [1]

    using the Force Autocorrelation Function (FACF) averaged over the ensemble: γ(t) = Q2 M kBT ⟨E(t)·E(0)⟩ ensemble.(22) This procedure yielded the non-Markovian friction ker- nel solely from equilibrium noise. In the second set of simulations, we introduced active test particles fully coupled to the Maxwell solver, allow- ing them to excite wakes and experi...

    work page

  2. [2]

    Onsager,Reciprocal Relations in Irreversible Pro- cesses

    L. Onsager,Reciprocal Relations in Irreversible Pro- cesses. I, Phys. Rev.37, 405 (1931)

    work page 1931

  3. [3]

    M. S. Green,Markoff Random Processes and the Sta- tistical Mechanics of Time-Dependent Phenomena. II. Irreversible Processes in Fluids, J. Chem. Phys.22, 398 7 (1954)

    work page 1954

  4. [4]

    H. B. Callen and T. A. Welton,Irreversibility and Gen- eralized Noise, Phys. Rev.83, 34 (1951)

    work page 1951

  5. [5]

    Kubo,The fluctuation-dissipation theorem, Rep

    R. Kubo,The fluctuation-dissipation theorem, Rep. Prog. Phys.29, 255 (1966)

    work page 1966

  6. [6]

    D. J. Evans and G. P. Morriss,Statistical Mechanics of Nonequilibrium Liquids(Cambridge University Press, Cambridge, 2008)

    work page 2008

  7. [7]

    Zwanzig,Memory Effects in Irreversible Thermody- namics, Phys

    R. Zwanzig,Memory Effects in Irreversible Thermody- namics, Phys. Rev.124, 983 (1961)

    work page 1961

  8. [8]

    Peter and J

    Th. Peter and J. Meyer-ter-Vehn,Energy loss of heavy ions in dense plasma, Phys. Rev. A43, 1998 (1991)

    work page 1998

  9. [9]

    Zwicknagel, C

    G. Zwicknagel, C. Toepffer, and P.-G. Reinhard,Stop- ping of heavy ions in plasmas at strong coupling, Phys. Rep.309, 117 (1999)

    work page 1999

  10. [10]

    Mori,Transport, Collective Motion, and Brownian Motion, Prog

    H. Mori,Transport, Collective Motion, and Brownian Motion, Prog. Theor. Phys.33, 423 (1965)

    work page 1965

  11. [11]

    Grabert,Projection Operator Techniques in Nonequilibrium Statistical Mechanics(Springer, Berlin, 1982)

    H. Grabert,Projection Operator Techniques in Nonequilibrium Statistical Mechanics(Springer, Berlin, 1982)

    work page 1982

  12. [12]

    N. G. van Kampen,Stochastic Processes in Physics and Chemistry(North-Holland, Amsterdam, 1992)

    work page 1992

  13. [13]

    C. W. Gardiner,Handbook of Stochastic Methods (Springer, Berlin, 1985)

    work page 1985

  14. [14]

    S. R. de Groot and P. Mazur,Non-Equilibrium Ther- modynamics(Dover, New York, 1984)

    work page 1984

  15. [15]

    Chandrasekhar,Dynamical Friction

    S. Chandrasekhar,Dynamical Friction. I. General Con- siderations, Astrophys. J.97, 255 (1943)

    work page 1943

  16. [16]

    Spitzer,Physics of Fully Ionized Gases(Interscience, New York, 1962)

    L. Spitzer,Physics of Fully Ionized Gases(Interscience, New York, 1962)

    work page 1962

  17. [17]

    E. M. Lifshitz and L. P. Pitaevskii,Physical Kinet- ics (Landau and Lifshitz Course of Theoretical Physics, Vol. 10)(Pergamon Press, Oxford, 1981)

    work page 1981

  18. [18]

    Ichimaru,Basic Principles of Plasma Physics(Ben- jamin/Cummings, Reading, MA, 1973)

    S. Ichimaru,Basic Principles of Plasma Physics(Ben- jamin/Cummings, Reading, MA, 1973)

    work page 1973

  19. [19]

    Balescu,Irreversible Processes in Ionized Gases, Phys

    R. Balescu,Irreversible Processes in Ionized Gases, Phys. Fluids3, 52 (1960)

    work page 1960

  20. [20]

    Lenard,On Bogoliubov’s Kinetic Equation for a Spa- tially Homogeneous Plasma, Ann

    A. Lenard,On Bogoliubov’s Kinetic Equation for a Spa- tially Homogeneous Plasma, Ann. Phys. (N.Y.)3, 242 (1960)

    work page 1960

  21. [21]

    Y. L. Klimontovich,Kinetic Theory of Nonideal Plas- mas(Academic Press, New York, 1975)

    work page 1975

  22. [22]

    T. S. Hahm and H. H. Kaang,Non-Markovian turbulent transport, Phys. Plasmas29, 102302 (2022)

    work page 2022

  23. [23]

    Daligault,Calculations of the electron-ion energy re- laxation rate in dense plasmas from equilibrium molec- ular dynamics simulations, Phys

    J. Daligault,Calculations of the electron-ion energy re- laxation rate in dense plasmas from equilibrium molec- ular dynamics simulations, Phys. Rev. E100, 043201 (2019)

    work page 2019

  24. [24]

    Hansen and I

    J.-P. Hansen and I. R. McDonald,Theory of Simple Liquids(Academic Press, London, 2013)

    work page 2013

  25. [25]

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

    work page 1987

  26. [26]

    M. S. Murillo,Viscosity of strongly coupled plasma, Phys. Rev. E62, 4115 (2000)

    work page 2000

  27. [27]

    C. E. Starrett, J. Daligault, and D. Saumon,Uncertain- ties in stopping power of dense plasmas, Phys. Rev. E 91, 013104 (2015)

    work page 2015

  28. [28]

    A. V. Ivlev et al.,Non-Newtonian viscosity of complex plasmas, Phys. Rev. Lett.98, 145003 (2007)

    work page 2007

  29. [29]

    Ichimaru,Statistical Plasma Physics, Vol

    S. Ichimaru,Statistical Plasma Physics, Vol. I: Basic Principles(Westview Press, Boulder, 2004)

    work page 2004

  30. [30]

    Kubo,Statistical-Mechanical Theory of Irreversible Processes

    R. Kubo,Statistical-Mechanical Theory of Irreversible Processes. I, J. Phys. Soc. Jpn.12, 570 (1957)

    work page 1957

  31. [31]

    Risken,The Fokker-Planck Equation(Springer- Verlag, Berlin, 1989)

    H. Risken,The Fokker-Planck Equation(Springer- Verlag, Berlin, 1989)

    work page 1989

  32. [32]

    T. D. Arberet al., Plasma Phys. Control. Fusion57, 113001 (2015)

    work page 2015

  33. [33]

    C. K. Birdsall and A. B. Langdon,Plasma Physics via Computer Simulation(Taylor & Francis, New York, 2004)

    work page 2004

  34. [34]

    S. X. Hu, V. N. Goncharov, P. B. Radha, S. P. Re- gan, and E. M. Campbell,A review on charged-particle transport modeling for laser direct-drive fusion, Phys. Plasmas31, 040501 (2024)

    work page 2024