arxiv: 2605.13703 · v1 · submitted 2026-05-13 · ❄️ cond-mat.mtrl-sci · cond-mat.stat-mech

Recognition: unknown

Nonlinear dynamic elastic moduli from equilibrium stress fluctuations

F. E. Garbuzov , Y. M. Beltukov

Authors on Pith no claims yet

Pith reviewed 2026-05-14 17:46 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci cond-mat.stat-mech

keywords fluctuation formulasnonlinear dynamic moduliviscoelastic responsetransient-time correlationsstress tensormolecular dynamicsanharmonic elasticityBorn terms

0 comments

The pith

Equilibrium stress fluctuations determine the nonlinear dynamic elastic moduli.

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

The paper derives transient-time correlation function expressions for the nonlinear time-dependent elastic moduli that control anharmonic viscoelastic response under finite strains. These expressions are obtained from equilibrium molecular dynamics by using the DOLLS/SLLOD equations of motion for irrotational flow and involve only time correlations of the stress tensor together with Born-kinetic terms. The same formulas recover the established quasi-static moduli and the linear dynamic moduli in the appropriate limits. A reader would care because the approach removes the need to impose explicit nonequilibrium deformation protocols when computing these material properties from simulation.

Core claim

Starting from the DOLLS/SLLOD equations of motion for irrotational flow, the authors obtain transient-time correlation function expressions for both linear and nonlinear dynamic moduli in terms of equilibrium time correlations of the stress tensor and Born-kinetic terms; these expressions recover the known quasi-static and linear dynamic results in the appropriate limits.

What carries the argument

Transient-time correlation functions of the stress tensor and Born-kinetic terms obtained from the DOLLS/SLLOD equations of motion for irrotational flow.

Load-bearing premise

The DOLLS/SLLOD equations of motion for irrotational motion correctly capture the nonlinear response under finite time-dependent strains.

What would settle it

Compute the nonlinear dynamic moduli from the derived equilibrium correlation formulas and compare them with the same moduli obtained from direct nonequilibrium simulations that apply the corresponding finite time-dependent strains; systematic mismatch would falsify the formulas.

read the original abstract

Fluctuation formulas for elastic and viscoelastic moduli allow their computation from equilibrium molecular dynamics simulations, avoiding explicit nonequilibrium deformation protocols. While such expressions are well established for the quasi-static moduli, and also the linear dynamic moduli, no fluctuation formula exists for the nonlinear time-dependent moduli that govern anharmonic viscoelastic response under finite time-dependent strains. In this work we derive transient-time correlation function expressions for both the linear and the nonlinear dynamic moduli, starting from the DOLLS/SLLOD equations of motion for irrotational motion. The resulting formulas involve equilibrium time correlations of the stress tensor and Born-kinetic terms, and they recover the known quasi-static and linear dynamic results in the appropriate limits.

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.

Desk Editor's Note private letter to a colleague

Derives new TTCF expressions for nonlinear dynamic moduli from equilibrium fluctuations using DOLLS/SLLOD, recovering known limits without added parameters.

read the letter

The main point is that this paper derives transient-time correlation function expressions for nonlinear dynamic elastic moduli. These describe anharmonic viscoelastic response under finite time-dependent strains, and the abstract states no such formulas existed before. The derivation starts from the DOLLS/SLLOD equations for irrotational motion and writes the moduli in terms of equilibrium correlations of the stress tensor and Born-kinetic terms. That is the actual advance. It does well by recovering the established quasi-static and linear dynamic results in the proper limits. This check is useful and shows internal consistency without introducing fitting parameters or circular steps. The approach stays within standard methods for equilibrium molecular dynamics, which avoids running explicit nonequilibrium deformations. The soft spots are modest but real. The work is presented as a derivation, so there are no numerical tests or example calculations shown to check convergence, accuracy at finite strains, or practical implementation details. The central assumption that the chosen equations of motion suffice for nonlinear response is common in the subfield, but it would help to see at least one concrete validation. No load-bearing flaws appear in the stated logic. This paper is for computational materials scientists who already use fluctuation methods for linear viscoelasticity and want to extend them to nonlinear cases in solids or soft matter. Readers who care about equilibrium routes to anharmonic properties will get direct value. The thinking is clear and the claim is specific, so the paper deserves a serious referee to examine the full equations and suggest validation steps.

Referee Report

0 major / 3 minor

Summary. The manuscript derives transient-time correlation function (TTCF) expressions for both the linear and nonlinear dynamic elastic moduli from equilibrium stress fluctuations. Starting from the DOLLS/SLLOD equations of motion under irrotational flow, the formulas are expressed in terms of equilibrium time correlations involving the stress tensor and Born-kinetic terms; the derivation is shown to recover the established quasi-static and linear dynamic limits in the appropriate regimes.

Significance. If the central derivation is correct, the work fills a notable gap by providing the first fluctuation-based route to nonlinear, time-dependent moduli governing anharmonic viscoelasticity under finite strains. This would enable efficient equilibrium MD computations of finite-strain dynamic response without explicit nonequilibrium deformation protocols, with clear utility for computational materials science.

minor comments (3)

The abstract and introduction refer to 'Born-kinetic terms' without an explicit definition or reference to their standard form; add a short paragraph in §2 or §3 that recalls their expression from the stress tensor derivative.
In the recovery of the linear limit (presumably around Eq. (22) or equivalent), the vanishing of the nonlinear contributions should be shown algebraically rather than stated, to make the reduction fully transparent.
The manuscript would benefit from a brief discussion of the irrotational-flow restriction and its implications for applicability to general shear or rotational flows, even if only as a limitations paragraph.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive summary and significance assessment of the derivation of TTCF expressions for linear and nonlinear dynamic elastic moduli. The recommendation for minor revision is noted; absent any specific major comments, we interpret this as a request for minor clarifications or improvements that we will address in the revised manuscript.

Circularity Check

0 steps flagged

Derivation from established DOLLS/SLLOD equations shows no circularity

full rationale

The paper derives TTCF expressions for linear and nonlinear dynamic moduli directly from the standard DOLLS/SLLOD equations of motion under irrotational strain. No self-definitional steps appear (no quantity defined in terms of the target result), no parameters are fitted to data and then relabeled as predictions, and no load-bearing self-citations or uniqueness theorems from the same authors are invoked. The formulas are shown to recover known quasi-static and linear limits, confirming the derivation adds independent content rather than reducing to its inputs by construction. This matches the default expectation of a non-circular derivation in the subfield.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the applicability of the DOLLS/SLLOD framework for irrotational motion and the assumption that nonlinear moduli are expressible via equilibrium stress correlations; no free parameters or new entities are mentioned.

axioms (1)

domain assumption DOLLS/SLLOD equations of motion for irrotational motion correctly describe the strained system
Explicitly stated as the starting point for the derivation in the abstract.

pith-pipeline@v0.9.0 · 5414 in / 1130 out tokens · 27563 ms · 2026-05-14T17:46:29.519075+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

23 extracted references

[1]

Strain fluctuations and elastic constants,

M. Parrinello and A. Rahman, “Strain fluctuations and elastic constants,”J. Chem. Phys., vol. 76, no. 5, pp. 2662–2666, 1982

1982
[2]

Molecular dynamics calculation of elastic constants for a crystalline system in equilibrium,

J. R. Ray, M. C. Moody, and A. Rahman, “Molecular dynamics calculation of elastic constants for a crystalline system in equilibrium,”Phys. Rev. B, vol. 32, no. 2, pp. 733– 735, 1985

1985
[3]

Generalized expressions for the calculation of elastic constants by computer simulation,

J. F. Lutsko, “Generalized expressions for the calculation of elastic constants by computer simulation,”J. Appl. Phys., vol. 65, no. 8, pp. 2991–2997, 1989

1989
[4]

Isothermal elastic constants for argon. Theory and Monte Carlo calculations,

D. Squire, A. Holt, and W. Hoover, “Isothermal elastic constants for argon. Theory and Monte Carlo calculations,”Physica, vol. 42, no. 3, pp. 388–397, 1969

1969
[5]

Statistical mechanics of elastic moduli,

F. Bavaud, Ph. Choquard, and J. R. Fontaine, “Statistical mechanics of elastic moduli,” J Stat Phys, vol. 42, no. 3-4, pp. 621–646, 1986

1986
[6]

Athermal nonlinear elastic constants of amorphous solids,

S. Karmakar, E. Lerner, and I. Procaccia, “Athermal nonlinear elastic constants of amorphous solids,”Phys. Rev. E, vol. 82, no. 2, p. 026105, 2010

2010
[7]

Viscoelastic properties of crystals,

S. R. Williams and D. J. Evans, “Viscoelastic properties of crystals,”J. Chem. Phys., vol. 131, no. 2, p. 024115, 2009

2009
[8]

Lennard-Jones triple-point bulk and shear viscosities. Green-Kubo theory, Hamiltonian 10 mechanics, and nonequilibrium molecular dynamics,

W. G. Hoover, D. J. Evans, R. B. Hickman, A. J. C. Ladd, W. T. Ashurst, and B. Moran, “Lennard-Jones triple-point bulk and shear viscosities. Green-Kubo theory, Hamiltonian 10 mechanics, and nonequilibrium molecular dynamics,”Phys. Rev. A, vol. 22, no. 4, pp. 1690–1697, 1980

1980
[9]

Nonlinear-response theory for steady planar Couette flow,

D. J. Evans and G. P. Morriss, “Nonlinear-response theory for steady planar Couette flow,”Phys. Rev. A, vol. 30, no. 3, pp. 1528–1530, 1984

1984
[10]

Non-Newtonian molecular dynamics,

D. J. Evans and G. Morriss, “Non-Newtonian molecular dynamics,”Computer Physics Reports, vol. 1, no. 6, pp. 297–343, 1984

1984
[11]

D. J. Evans and G. P. Morriss,Statistical mechanics of nonequilibrium liquids. Cambridge: Cambridge University Press, 2nd ed., 2008

2008
[12]

On the number dependence of viscosity in three dimensional fluids,

D. J. Evans, G. Morriss, and L. Hood, “On the number dependence of viscosity in three dimensional fluids,”Molecular Physics, vol. 68, no. 3, pp. 637–646, 1989

1989
[13]

Self-diffusion in equilibrium and sheared liquid mixtures by molecular dynamics,

W. Sandberg and D. Heyes, “Self-diffusion in equilibrium and sheared liquid mixtures by molecular dynamics,”Molecular Physics, vol. 85, no. 3, pp. 635–649, 1995

1995
[14]

Decane under shear: A molecular dynamics study using reversible NVT-SLLOD and NPT-SLLOD algorithms,

C. J. Mundy, J. I. Siepmann, and M. L. Klein, “Decane under shear: A molecular dynamics study using reversible NVT-SLLOD and NPT-SLLOD algorithms,”The Journal of Chemical Physics, vol. 103, no. 23, pp. 10192–10200, 1995

1995
[15]

Molecular dynamics simulation of entangled polymers in shear flow,

T. Aoyagi and M. Doi, “Molecular dynamics simulation of entangled polymers in shear flow,”Computational and Theoretical Polymer Science, vol. 10, no. 3-4, pp. 317–321, 2000

2000
[16]

Unraveling the pressure- viscosity behavior and shear thinning in glycerol using atomic scale molecular dynamics simulations,

V. Fadaei Naeini, M. Bj¨ orling, J. Larsson, and R. Larsson, “Unraveling the pressure- viscosity behavior and shear thinning in glycerol using atomic scale molecular dynamics simulations,”Journal of Molecular Liquids, vol. 390, p. 122990, 2023

2023
[17]

A proper approach for nonequilibrium molecular dynamics simulations of planar elongational flow,

C. Baig, B. J. Edwards, D. J. Keffer, and H. D. Cochran, “A proper approach for nonequilibrium molecular dynamics simulations of planar elongational flow,”J. Chem. Phys., vol. 122, no. 11, p. 114103, 2005

2005
[18]

A simple, direct derivation and proof of the validity of the SLLOD equations of motion for generalized homogeneous flows,

P. J. Daivis and B. D. Todd, “A simple, direct derivation and proof of the validity of the SLLOD equations of motion for generalized homogeneous flows,”The Journal of Chemical Physics, vol. 124, no. 19, p. 194103, 2006

2006
[19]

The Statistical Mechanical Theory of Transport Processes. IV. The Equations of Hydrodynamics,

J. H. Irving and J. G. Kirkwood, “The Statistical Mechanical Theory of Transport Processes. IV. The Equations of Hydrodynamics,”J. Chem. Phys., vol. 18, no. 6, pp. 817–829, 1950

1950
[20]

Crystal Structure and Pair Potentials: A Molecular- Dynamics Study,

M. Parrinello and A. Rahman, “Crystal Structure and Pair Potentials: A Molecular- Dynamics Study,”Phys. Rev. Lett., vol. 45, no. 14, pp. 1196–1199, 1980

1980
[21]

Generalization of nonlinear Murnaghan elastic model for viscoelastic materials,

F. Garbuzov and Y. Beltukov, “Generalization of nonlinear Murnaghan elastic model for viscoelastic materials,”International Journal of Non-Linear Mechanics, vol. 159, p. 104598, 2024. 11

2024
[22]

Small Finite Deformations of Viscoelastic Solids,

A. C. Pipkin, “Small Finite Deformations of Viscoelastic Solids,”Rev. Mod. Phys., vol. 36, no. 4, pp. 1034–1041, 1964

1964
[23]

W. N. Findley, J. S. Lai, and K. Onaran,Creep and relaxation of nonlinear viscoelastic materials, with an introduction to linear viscoelasticity. Amsterdam: Elsevier, 1976. 12

1976