Your new experience awaits. Try the new design now and help us make it even better

BRIEF RESEARCH REPORT article

Front. Phys., 21 October 2025

Sec. Statistical and Computational Physics

Volume 13 - 2025 | https://doi.org/10.3389/fphy.2025.1667224

Equilibrium trajectories quantify second-order violations of the fluctuation–dissipation theorem without the need for a model

Juliana Caspers
Juliana Caspers1*Karthika Krishna KumarKarthika Krishna Kumar2Clemens BechingerClemens Bechinger2Matthias Krüger
Matthias Krüger1*
  • 1Institute for Theoretical Physics, Georg-August-Universität Göttingen, Göttingen, Germany
  • 2Fachbereich Physik, Universität Konstanz, Konstanz, Germany

Quantifying and characterizing fluctuations far away from equilibrium is a challenging task. We discuss and experimentally confirm a series expansion for a driven classical system, relating the different nonequilibrium cumulants of the observable conjugate to the driving protocol. This series is valid from micro- to macroscopic length scales, and it encompasses the fluctuation–dissipation theorem (FDT). We apply it in experiments of a Brownian probe particle confined and driven by an optical potential and suspended in a nonlinear and non-Markovian fluid. The expansion states that the form of the FDT remains valid away from equilibrium for Gaussian observables, up to the order presented. We show that this expansion agrees with that of a known fluctuation theorem up to an unresolved difference regarding moments versus cumulants.

Introduction

The fluctuation–dissipation theorem (FDT) [1, 2], connecting response and fluctuations of equilibrium systems, is of fundamental importance for condensed matter, fluids, plasmas, or electromagnetic fields [36]. One of its remarkable properties is its validity at any length scale, ranging from the nanoscale, for electric charges, to the macroscale, for macroscopic magnetization. It is, however, restricted to the linear regime, i.e., to situations close to equilibrium. Most previous research has been largely devoted to determining similar relations for nonequilibrium steady states [729] and for nonlinear responses [3051]. A typical observation in the found relations is the explicit appearance of microscopic details—sometimes referred to as frenetic components [52, 53] or information on the specific rule governing the time evolution [40]—often hampering a model-independent formulation and systematic changes in length scales such as coarse graining to macroscopic scales [46, 47, 49]. As a consequence, experimental tests and application of such relations have indeed been successful for systems with a small number of accessible Markovian degrees of freedom [51, 5457], for which the dynamics can be modeled.

In a different spirit, nonlinear fluctuation dissipation relations [3134, 58] and fluctuation theorems [33, 41, 5961] have been found, which can often be applied in the absence of a specific model. However, they have, to our knowledge, not been used to quantify the error of the FDT.

In this manuscript, we discuss and experimentally confirm a series expansion for a driven classical system, which relates the different nonequilibrium cumulants of the observable conjugate to the driving protocol, up to a certain order in driving velocity. This series (i) is valid from micro- to macroscopic length scales, (ii) is model-independent, and (iii) encompasses the FDT. We apply it in an experimental many-body system of a Brownian probe particle interacting with worm-like micelles and confined and driven by an optical potential. In these experiments, we demonstrate that the equilibrium third force cumulant quantifies the second-order deviation from the FDT under driving. Notably, our theoretical predictions demonstrate that the form of the FDT remains valid for purely Gaussian observables within the displayed order.

System and fluctuation series

Consider a classical system of stochastic degrees xt at time t (weakly) coupled to a heat bath at temperature T. The system’s potential energy U depends on xt and on a time-dependent deterministic protocol Xt, i.e., U(xt,Xt). The system is prepared in equilibrium at time t, with a protocol value X=Xt, for simplicity [62]. The time dependence of the protocol drives the system away from equilibrium.

The derivative of U in terms of Xt, Ft:=XtU(xt,Xt) is the observable conjugate to Xt. For example, if Xt is a position as in our experiments, Ft is (minus) the corresponding force. Ft can be micro- or macroscopic; for example, let X couple linearly to an observable A(x); U(xt,Xt)=U(xt,0)XtA(xt), i.e., Ft=A(xt). Thus, if A(x) is a macroscopic field, such as macroscopic magnetization, Ft is macroscopic. If A(x) is the position of a molecular particle, Ft is microscopic. The following remains valid if Xt enters U nonlinearly.

The statistical properties of Ft in this nonequilibrium situation are encoded in its cumulants and in its correlations with another state observable, Bt=B(xt,Xt), which we aim to study here. The well-known FDT connects the covariance in the unperturbed system and the first moment of B under weak driving [1, 2],

β2tdsẊsBt;Fseq=βBtBteq+O(Ẋ2),(1)

with β=1/kBT and Boltzmann constant kB. ; denotes the second cumulant, and similarly for higher orders below. Expectation values eq are evaluated using equilibrium trajectories with a protocol value fixed at Xt. Expectation values are measured in the driven system, i.e., under the given time-dependent protocol [62].

In Ref. [62], we derive identities connecting the nonequilibrium cumulants of Ft and Bt to different orders. These give rise to the following series involving the mentioned cumulants [62],

β2tdsẊsBt;Fs=βBtBteq+β32tdstdsẊsẊsBt;Fs;Fsβ46tdstdstdsẊsẊsẊs×Bt;Fs;Fs;Fs+O(Ẋ4),(2)

As presented in Ref. [62], this series expansion can also be obtained from a known fluctuation theorem [33, 41], albeit with an open question regarding cumulants versus moments.

Equation 2 is, as indicated, correct up to the fourth order in driving Ẋt, under the assumption of local detailed balance [63]. Expanding1 Equation 2 to the first order yields the FDT in Equation 1 so that it is included in Equation 2. To higher orders, first and second cumulants do not fulfill the FDT, and Equation 2 quantifies their difference in terms of third and fourth cumulants of F and B. Notably, the second to fourth lines of Equation 2 vanish for purely Gaussian distributed F and B so that first and second cumulants obey the FDT to the given nonequilibrium order. It is important to note that, in Equation 2, the protocol Ẋ appears as pre-factors and in the nonequilibrium cumulants themselves; i.e., the latter are evaluated under application of driving. As Equation 2 only requires measurement of F and B, we use the notion of model free.

Experimental setup

We exploit Equation 2 with experiments of Brownian particles interacting with micellar fluid. In particular, we use silica particles of diameter 1 µm suspended in a 5-mM equimolar solution of cetylpyridinium chloride monohydrate (CPyCl) and sodium salicylate (NaSal). At concentrations above the critical micellar concentration ( 4 mM), this fluid is known to form giant worm-like micelles, leading to a viscoelastic nonlinear behavior at ambient temperatures [64]; see SM. At 5 mM, we determine the relaxation time of the fluid from microrheological recoil experiments, where a particle is first driven with a constant external force which is then suddenly removed, to be (3±0.2)s [65]. A small amount of silica particles is added to the micellar solution, which is contained in a rectangular capillary with 100 µm height and kept at a temperature of 25 °C. This sample is placed on a custom-built optical tweezer setup that uses a Gaussian laser beam of wavelength 532 nm and a 100× oil immersion objective (NA = 1.45). The laser beam yields a potential U(xtXt), as shown in Figure 1a, centered at Xt, trapping one of the silica particles with coordinate xt. Ft=XtU(xtXt) is thus the force acting on the particle by the trapping potential (or vice versa). As the micellar degrees do not couple to X, they do not enter Ft explicitly, and knowing how they enter U is not required to apply Equation 2. Thus, the use of Equation 2 does not require the detection of the positions of micellar particles, and applying it to such a complex fluid demonstrates its strength. We consider BF, and made the potential asymmetric, to obtain a finite second-order response of βFt. This allows to test Equation 2 to the second order in our experiments, and it is achieved by a controlled lateral displacement of the vertically incident laser beam from the center of the objective lens (see SM). Notably, Equation 2 could also be tested in a purely viscous fluid, which also shows a second-order response due to the nonlinear potential. However, to demonstrate that Equation 2 is valid beyond simple systems, we have chosen the more challenging case of a micellar fluid.

Figure 1
Four-panel scientific chart illustrating different graphs. Panel (a) shows a potential energy profile with data points, including an inset with a schematic diagram of a particle confined in an asymmetric optical potential. Panel (b) depicts multiple lines showing β〈Fᵗ〉 versus time for different amplitudes ranging from 0.03 to 0.14 micrometers. Panel (c) displays variations in β²\(\tilde X\)ω(Fᵗ,\(\tilde F\)ᵗ). Panel (d) displays variations in third cumulant as a function of time. Panel (e) validates Equation 4 with curves representing first cumulant, second cumulant and sum of first cumulant and the third cumulant.

Figure 1. (a) Asymmetric optical potential U(x)=kBTlnP(x) felt by the probe particle (inset sketch), with P(x) the probability distribution with the trap at rest. (b) Mean force βFt, (c) force covariance β2X̂ωF̃t,Ft, and (d) equilibrium third cumulant β2(X̂ω)2F̃t;F̃t;Fteq/2, as functions of time, for driving frequency ω=8.4 rad/s and amplitudes X̂={0.03,0.06,0.08,0.09,0.14}µm as labeled. Tp=2πω. (e) Force covariance (solid line), mean force (dotted line), and the sum of mean force and third force cumulant (dashed line, Equation 4) for X̂ = 0.14 µm.

To apply the driving protocol, the sample cell is moved, whereas the optical trap remains stationary in our experiments. This is achieved using a piezo-driven stage, on which the sample is mounted and translated in an oscillating manner relative to the trap. In the fluid’s rest frame, this yields a periodic motion of the potential minimum Xt, i.e., the protocol,

Xt=X̂sin(ωt),(3)

with the amplitude X̂ and the frequency ω. Particle trajectories are recorded with a frame rate of 150 Hz using a video camera, and particle positions are determined using a custom MATLAB algorithm. To yield sufficient statistics, each protocol (X̂, ω) was measured over 1400s. We allowed the system to reach a steady state by recording trajectories only after at least five oscillation periods had passed. Thereafter, no further equilibration was visible in the data. Prior to each nonequilibrium protocol, we recorded particle trajectories for another 1000s with Xt at rest. These equilibrium data were used to check that the form of U does not vary between measurements and also to obtain the force cumulants under equilibrium conditions.

Data analysis

With the protocol of Equation 3, Equation 2 takes, expanded to the second order, the form

β2X̂ωF̃t;Ft=βFt+β3X̂2ω22F̃t;F̃t;Ft+O((X̂ω)3),(4)

where the tilde denotes the cosine transform, i.e., F̃ttdscos(ωs)Fs. We restrict the analysis to the lowest nontrivial, i.e., second, order and expand Equation 2 accordingly, also using Fteq=0. Notably, to extract the second-order contribution from the last term in Equation 4, the third cumulant in Equation 4 is replaced by its equilibrium version. This is because of the pre-factor X̂2ω2 and will be done in the following analysis.

The cumulants in Equation 4 depend on time t in a periodic manner, as shown in Figures 1b–d, for ω= 8.4 rad/s and driving2 amplitudes ranging from X̂= 0.03 µm (light) to X̂= 0.14 µm (dark). Figure 1b shows the mean force, which, as expected for a driven oscillator, is a periodic function with period Tp=2πω. For the smallest amplitude shown, the mean force is nearly harmonic with frequency ω, as expected from linear response. With growing amplitude, higher harmonics occur, as expected from nonlinear response. This asymmetric system shows the second-order response with expected frequencies of 2ω and 0ω.

Figure 1c shows the force covariance for the same parameters and color code. For small amplitude X̂, the curves in Figures 1b,c are equal within the experimental accuracy, as analyzed in detail below. For larger driving amplitude, the force covariance develops higher harmonics with signatures of second order. Very little is known about the properties of such nonequilibrium fluctuations, and quantifying these is difficult. It is notable that the curves in Figure 1c, for larger amplitudes, deviate from Figure 1b, a deviation that we claim to be quantified by Equation 4.

Figure 1d shows the third cumulant of force for the same parameters and color code. We have here restricted to the equilibrium cumulant as it appears in Equation 4, multiplied by (X̂ω)2. The curves in Figure 1d thus differ only because of the factor (X̂ω)2. They thus scale quadratically in driving velocity and only show frequencies of 2ω and 0ω.

Equation 4 states that in the shown range of amplitudes, the curves in Figure 1c are given by the sum of the curves in Figures 1b,d. For X̂=0.14 µm, the respective summed curve is shown as a dashed line together with the mean force and the force covariance in Figure 1e. The agreement is convincing and a confirmation of Equation 4.

To test this prediction systematically, we dissect the curves in Figures 1b–d into the contributions from harmonics with frequencies 0ω, ω, and 2ω, respectively, i.e., we expand the cumulant of order n into harmonics with the frequency mω,

βn(X̂ω)n1(F̃t;)n1Ft=m=0Am(n)sin(mωt+ϕm(n)),(5)

where the coefficients Am(n) depend on X̂ω. We set ϕ0(n)π/2 for consistency. Equation 4, projected on the harmonic of order m, yields relations between coefficients and phases for each m, which we can test.

Results

Figure 2 shows the coefficients Am(n) as a function of the driving amplitude X̂. The top panel gives the order m=1, which is considered to be linear in X̂ for the range shown, as expected from linear response. The graph shows the mean force (data points) and the force covariance (line). The latter is evaluated from equilibrium trajectories. The agreement in this panel, for the range of sufficiently small X̂, is expected from the FDT.

Figure 2
Three-panel graph showing data trends with error bars. Top panel: \( A_1\) versus \( \hat{X} \), blue squares (\( A_1^(1)\)) with line (\( A_1^(2)\)). Middle panel: \(A_0^{(1)}-A_0^{(2)}\), red circles, and \(-A_0^{(3)}\), black line with gray error band versus, \(\hat{X}\). Bottom panel: \( \Delta A_2 \), red circles, and \( -A_2^(3)\), gray error band, versus \( \hat{X} \). All panels show agreement between symbols and lines.

Figure 2. Coefficients Am corresponding to harmonics with the frequency mω, as a function of driving amplitude X̂, for ω=8.4 rad/s. Top panel: m=1, demonstrating FDT. Lower panels: difference in first and second cumulants (data points), and third cumulant (lines) for m=0 and m=2. The agreement confirms Equation 4. ΔA2[(A2(1))2+(A2(2))22A2(1)A2(2)cos(ϕ2(1)ϕ2(2))]1/2. Error bars and bands are obtained from partitioning trajectories into two pieces.

The center and lower panels in Figure 2 show the orders m=0 and m=2, respectively. In particular, these panels present the difference of first and second cumulants in Equation 4 (data points), i.e., βFtβ2X̂ωF̃t;Ft, together with the third cumulant (line), β3(X̂ω)2F̃t;F̃t;Fteq/2, evaluated from equilibrium measurements3. The latter is shown as a parabola with curvature obtained from the third force cumulant at equilibrium. The data points in this graph thus quantify the deviation from the FDT, with the line giving the prediction of Equation 4 for this deviation. The agreement is convincing for both m=0 and m=2, supporting the validity of Equation 4.

As data in the top panel of Figure 2 grow linearly and those in the center and lower panels grow quadratically with X̂, we fit a line and a parabola to obtain the, respective, slope and curvature for each m. The obtained values—divided by the respective power in ω—are shown in Figure 3 as a function of the frequency ω. We observe convincing agreement for the measured frequencies further supporting Equation 4. In tendency, the coefficients decrease with increasing ω. Notably, the statistical accuracy of the data points decreases with decreasing ω; with smaller ω, longer trajectories are required, especially in a fluid with pronounced memory, as the period of the cosine transform increases.

Figure 3
Three stacked graphs display data with curves and error bars. The top graph shows blue squares representing \( A_1^{(1)} \) and a black line for \( A_1^{(2)} \). The middle graph features red circles with error bars for \( A_0^{(1)} - A_0^{(2)} \) and a blackline with gray error band for \(-A_0^{(3)}\). The bottom graph has red circles with error bars for \(\Delta A_2\) and a black line with gray error band for \(-A_2^{(3)}\). The x-axis across all graphs marks frequency \(\omega\) in radians per second.

Figure 3. Coefficients Am, normalized as labeled, as a function of the frequency ω. The agreement in the top panel confirms the FDT, and the agreement in the lower panels confirms Equation 4. Each data point is obtained from averaging over driving amplitudes taking the respective scaling of Am with X̂ into account (see SM). Error bars are obtained from partitioning trajectories into two pieces (data points) and from the standard deviation between separate series of measurements (gray area).

Figure 4 provides the final test of Equation 4, namely, the phases ϕm(n) of Equation 5. These hardly depend on driving amplitude, and the shown data are averaged over the measured values of X̂. The top panel of Figure 4 shows the order m=1, i.e., the linear response, with convincing agreement. The phase angle for m=1 is small for the frequencies measured, indicating that the force Ft is almost in phase with the protocol Xt, as for an elastic material. The black curve, extracted from equilibrium data, shows that with smaller ω, the phase increases, presumably reaching π/2 in the limit of ω0. The slow increase with decreasing ω displays the slow nature of the investigated system.

Figure 4
Two graphs display phase angles versus angular frequency. The top graph shows a decreasing trend of \(\phi_1\) with blue squares as data points. The bottom graph shows \(\phi_2\) with a shaded area, red circles with error bars, indicating fluctuations around a steady line. Angular frequency is on a logarithmic scale for both graphs.

Figure 4. Phase angles ϕ1 and ϕ2 of the harmonics in Equation 5, as functions of ω. Top panel shows the phases of the linear response. Lower panel shows the phase of the second-order response, comparing the contributions of the terms in Equation 4. Δϕ2arctanA2(1)sinϕ2(1)A2(2)sinϕ2(2)A2(1)cosϕ2(1)A2(2)cosϕ2(2). Error bars are obtained from partitioning trajectories into two pieces. The gray error band is obtained as the standard deviation between separate series of measurements.

The lower panel shows the phase for m=2, i.e., to second order. Although the agreement between the line and data confirms Equation 4, the graph shows that the differences in phases of first and second cumulants are rather small. In other words, first and second cumulants deviate noticeable in amplitude, as shown in Figure 2, but not so much in phase.

Conclusion

We presented and tested a nonequilibrium fluctuation expansion for a driven classical system, emphasizing the validity on various length scales. Indeed, such relations are necessary, e.g., for a systematic coarse graining of nonequilibrium systems. The identity is confirmed for experiments of a Brownian particle interacting with a complex surrounding. Future work can explore other systems and aim to clarify the relation to the mentioned fluctuation theorems [33, 41]. It is also important to investigate the use of Equation 2 for treating systems far away from equilibrium.

This project was funded by the Deutsche Forschungsgemeinschaft (DFG), Grant No. SFB 1432 (Project ID 425217212)—Project C05.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions

JC: Writing – review and editing, Writing – original draft. KK: Writing – original draft, Writing – review and editing. CB: Writing – original draft, Writing – review and editing. MK: Writing – review and editing, Writing – original draft.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This project was funded by the Deutsche Forschungsgemeinschaft (DFG), Grant No. SFB 1432 (Project ID 425217212)—Project C05. The authors acknowledge support by the Open Access Publication Funds of the Göttingen University.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2025.1667224/full#supplementary-material

Footnotes

1The expansion in Equation 2 suggests the dimensionless expansion parameter βFtτtdsẋs with cumulant relaxation time τ.

2ω is determined from the power spectral density, thus carrying an error depending on the length of the measurement.

3As the phases ϕ2(1) and ϕ2(2) may differ, the coefficient ΔA2 of the difference of first and second cumulants is found via ΔA2A2(1)2+A2(2)22A2(1)A2(2)cosϕ2(1)ϕ2(2). For m = 0, i.e., the zero frequency contribution, there is no phase by definition, and the coefficients A0(n) can be compared directly.

References

1. Callen HB, Welton TA. Irreversibility and generalized noise. Phys Rev (1951) 83:34–40. doi:10.1103/physrev.83.34

CrossRef Full Text | Google Scholar

2. Kubo R. The fluctuation-dissipation theorem. Rep Prog Phys (1966) 29:255–84. doi:10.1088/0034-4885/29/1/306

CrossRef Full Text | Google Scholar

3. Landau LD, Lifshitz EM. Statistical physics. Pergamon (1980).

Google Scholar

4. Sitenko AG. Electromagnetic fluctuations in plasma. New York: Academic Press (1967).

Google Scholar

5. Rytov SM, Kravtsov IA, Tatarskii V. Principles of statistical radiophysics. Springer-Verlag (1987).

Google Scholar

6. Hansen JP, McDonald IR. Theory of simple liquids: with applications to soft matter. Elsevier Science (2013).

Google Scholar

7. Cugliandolo LF, Dean DS, Kurchan J. Fluctuation-dissipation theorems and entropy production in relaxational systems. Phys Rev Lett (1997) 79:2168–71. doi:10.1103/physrevlett.79.2168

CrossRef Full Text | Google Scholar

8. Ruelle D. General linear response formula in statistical mechanics, and the fluctuation-dissipation theorem far from equilibrium. Phys Lett A (1998) 245:220–4. doi:10.1016/s0375-9601(98)00419-8

CrossRef Full Text | Google Scholar

9. Crisanti A, Ritort F. Violation of the fluctuation–dissipation theorem in glassy systems: basic notions and the numerical evidence. J Phys A: Math Gen (2003) 36:R181–290. doi:10.1088/0305-4470/36/21/201

CrossRef Full Text | Google Scholar

10. Harada T, Sasa S-I. Equality connecting energy dissipation with a violation of the fluctuation-response relation. Phys Rev Lett (2005) 95:130602. doi:10.1103/physrevlett.95.130602

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Speck T, Seifert U. Restoring a fluctuation-dissipation theorem in a nonequilibrium steady state. EPL (2006) 74:391–6. doi:10.1209/epl/i2005-10549-4

CrossRef Full Text | Google Scholar

12. Deutsch JM, Narayan O. Energy dissipation and fluctuation response for particles in fluids. Phys Rev E (2006) 74:026112. doi:10.1103/physreve.74.026112

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Chetrite R, Falkovich G, Gawedzki K. Fluctuation relations in simple examples of non-equilibrium steady states. J Stat Mech (2008):P08005. doi:10.1088/1742-5468/2008/08/P08005

CrossRef Full Text | Google Scholar

14. Marconi UMB, Puglisi A, Rondoni L, Vulpiani A. Fluctuation–dissipation: response theory in statistical physics. Phys Rep (2008) 461:111–95. doi:10.1016/j.physrep.2008.02.002

CrossRef Full Text | Google Scholar

15. Saito K. Energy dissipation and fluctuation response in driven quantum langevin dynamics. Europhys Lett (2008) 83:50006. doi:10.1209/0295-5075/83/50006

CrossRef Full Text | Google Scholar

16. Baiesi M, Maes C, Wynants B. Fluctuations and response of nonequilibrium states. Phys Rev Lett (2009) 103:010602. doi:10.1103/physrevlett.103.010602

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Prost J, Joanny J-F, Parrondo JMR. Generalized fluctuation-dissipation theorem for steady-state systems. Phys Rev Lett (2009) 103:090601. doi:10.1103/physrevlett.103.090601

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Harada T. Macroscopic expression connecting the rate of energy dissipation with the violation of the fluctuation response relation. Phys Rev E (2009) 79:030106. doi:10.1103/physreve.79.030106

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Krüger M, Fuchs M. Fluctuation dissipation relations in stationary states of interacting brownian particles under shear. Phys Rev Lett (2009) 102:135701. doi:10.1103/physrevlett.102.135701

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Seifert U, Speck T. Fluctuation-dissipation theorem in nonequilibrium steady states. EPL (2010) 89:10007. doi:10.1209/0295-5075/89/10007

CrossRef Full Text | Google Scholar

21. Baiesi M, Maes C, Wynants B. The modified sutherland–einstein relation for diffusive non-equilibria. Proc R Soc A: Math Phys Eng Sci (2011) 467:2792–809. doi:10.1098/rspa.2011.0046

CrossRef Full Text | Google Scholar

22. Cugliandolo LF. The effective temperature. J Phys A: Math Theor (2011) 44:483001. doi:10.1088/1751-8113/44/48/483001

CrossRef Full Text | Google Scholar

23. Verley G, Chétrite R, Lacoste D. Modified fluctuation-dissipation theorem for general non-stationary states and application to the glauber–ising chain. J Stat Mech (2011):10025. doi:10.1088/1742-5468/2011/10/P10025

CrossRef Full Text | Google Scholar

24. Altaner B, Polettini M, Esposito M. Fluctuation-dissipation relations far from equilibrium. Phys Rev Lett (2016) 117:180601. doi:10.1103/physrevlett.117.180601

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Lippiello E, Baiesi M, Sarracino A. Nonequilibrium fluctuation-dissipation theorem and heat production. Phys Rev Lett (2014) 112:140602. doi:10.1103/physrevlett.112.140602

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Wu W, Wang J. Generalized fluctuation-dissipation theorem for non-equilibrium spatially extended systems. Front Phys (2020) 8:567523. doi:10.3389/fphy.2020.567523

CrossRef Full Text | Google Scholar

27. Caprini L. Generalized fluctuation–dissipation relations holding in non-equilibrium dynamics. J Stat Mech (2021):063202. doi:10.1088/1742-5468/abffd4

CrossRef Full Text | Google Scholar

28. Baldovin M, Caprini L, Puglisi A, Sarracino A, Vulpiani A. In nonequilibrium thermodynamics and fluctuation kinetics: modern trends and open questions. Springer International Publishing (2022).

Google Scholar

29. Johnsrud MK, Golestanian R, Generalized fluctuation dissipation relations for active field theories (2024), arXiv:2409.14977.

Google Scholar

30. Kubo R. Statistical-mechanical theory of irreversible processes. I. General theory and simple applications to magnetic and conduction problems. J Phys Soc Jpn (1957) 12:570–86. doi:10.1143/jpsj.12.570

CrossRef Full Text | Google Scholar

31. Bernard W, Callen HB. Irreversible thermodynamics of nonlinear processes and noise in driven systems. Rev Mod Phys (1959) 31:1017–44. doi:10.1103/revmodphys.31.1017

CrossRef Full Text | Google Scholar

32. Efremov GF, Eksp Z. A fluctuation dissipation theorem for nonlinear media. Theor Fiz (1968) 55:2322.

Google Scholar

33. Bochkov GN, Kuzovlev YE. Nonlinear fluctuation-dissipation relations and stochastic models in nonequilibrium thermodynamics. Phys A: Stat Mech Appl (1981) 106:443–79. doi:10.1016/0378-4371(81)90122-9

CrossRef Full Text | Google Scholar

34. Stratonovich RL. Nonlinear nonequilibrium thermodynamics I: Linear and nonlinear fluctuation-dissipation theorems. Springer (1992).

Google Scholar

35. Evans DJ, Morriss G. Statistical mechanics of nonequilibrium liquids. 2nd ed. Cambridge University Press (2008).

Google Scholar

36. Fuchs M, Cates ME. Integration through transients for Brownian particles under steady shear. J Phys Condens Matter (2005) 17:S1681–96. doi:10.1088/0953-8984/17/20/003

CrossRef Full Text | Google Scholar

37. Holsten T, Krüger M. Thermodynamic nonlinear response relation. Phys Rev E (2021) 103:032116. doi:10.1103/PhysRevE.103.032116

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Oppenheim I. Nonlinear response and the approach to equilibrium. Prog Theor Phys Supp (1989) 99:369–81. doi:10.1143/ptps.99.369

CrossRef Full Text | Google Scholar

39. Bouchaud J-P, Biroli G. Nonlinear susceptibility in glassy systems: a probe for cooperative dynamical length scales. Phys Rev B (2005) 72:064204. doi:10.1103/physrevb.72.064204

CrossRef Full Text | Google Scholar

40. Lippiello E, Corberi F, Sarracino A, Zannetti M. Nonlinear response and fluctuation-dissipation relations. Phys Rev E (2008) 78:041120. doi:10.1103/physreve.78.041120

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Andrieux D, Gaspard P. Quantum work relations and response theory. Phys Rev Lett (2008) 100:230404. doi:10.1103/physrevlett.100.230404

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Lucarini V, Colangeli M. Beyond the linear fluctuation-dissipation theorem: the role of causality. J Stat Mech (2012):p05013. doi:10.1088/1742-5468/2012/05/P05013

CrossRef Full Text | Google Scholar

43. Diezemann G. Nonlinear response theory for markov processes: simple models for glassy relaxation. Phys Rev E (2012) 85:051502. doi:10.1103/physreve.85.051502

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Wang E, Heinz U. Generalized fluctuation-dissipation theorem for nonlinear response functions. Phys Rev D (2002) 66:025008. doi:10.1103/physrevd.66.025008

CrossRef Full Text | Google Scholar

45. Andrieux D, Gaspard P. A fluctuation theorem for currents and non-linear response coefficients. J Stat Mech (2007):P02006. doi:10.1088/1742-5468/2007/02/P02006

CrossRef Full Text | Google Scholar

46. Colangeli M, Maes C, Wynants B. A meaningful expansion around detailed balance. J Phys A: Math Theor (2011) 44:095001. doi:10.1088/1751-8113/44/9/095001

CrossRef Full Text | Google Scholar

47. Basu U, Krüger M, Lazarescu A, Maes C. Frenetic aspects of second order response. Phys Chem Chem Phys (2015) 17:6653. doi:10.1039/C4CP04977B

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Basu U, Helden L, Krüger M. Extrapolation to nonequilibrium from coarse-grained response theory. Phys Rev Lett (2018) 120:180604. doi:10.1103/PhysRevLett.120.180604

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Maes C. Response theory: a trajectory-based approach. Front Phys (2020) 8:229. doi:10.3389/fphy.2020.00229

CrossRef Full Text | Google Scholar

50. Caspers J, Krüger M. Nonlinear langevin functionals for a driven probe. J Chem Phys (2024) 161:124109. doi:10.1063/5.0227674

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Helden L, Basu U, Krüger M, Bechinger C. Measurement of second-order response without perturbation. EPL (2016) 116:60003. doi:10.1209/0295-5075/116/60003

CrossRef Full Text | Google Scholar

52. Maes C. On the second fluctuation–dissipation theorem for nonequilibrium baths. J Stat Phys (2014) 154:705–22. doi:10.1007/s10955-013-0904-8

CrossRef Full Text | Google Scholar

53. Maes C. Frenesy: time-Symmetric dynamical activity in nonequilibria. Phys Rep (2020) 850:1–33. doi:10.1016/j.physrep.2020.01.002

CrossRef Full Text | Google Scholar

54. Gomez-Solano JR, Petrosyan A, Ciliberto S, Chetrite R, Gawedzki K. Experimental verification of a modified fluctuation-dissipation relation for a micron-sized particle in a nonequilibrium steady state. Phys Rev Lett (2009) 103:040601. doi:10.1103/physrevlett.103.040601

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Blickle V, Speck T, Lutz C, Seifert U, Bechinger C. Einstein relation generalized to nonequilibrium. Phys Rev Lett (2007) 98:210601. doi:10.1103/physrevlett.98.210601

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Mehl J, Blickle V, Seifert U, Bechinger C. Experimental accessibility of generalized fluctuation-dissipation relations for nonequilibrium steady states. Phys Rev E (2010) 82:032401. doi:10.1103/physreve.82.032401

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Gomez-Solano JR, Petrosyan A, Ciliberto S, Maes C. Fluctuations and response in a non-equilibrium micron-sized system. J Stat Mech (2011):P01008. doi:10.1088/1742-5468/2011/01/P01008

CrossRef Full Text | Google Scholar

58. Bochkov GN, Kuzovlev YE. General theory of thermal fluctuations in nonlinear systems. Zh Eksp Teor Fiz (1977) 72:238.

Google Scholar

59. Jarzynski C. Nonequilibrium equality for free energy differences. Phys Rev Lett (1997) 78:2690–3. doi:10.1103/physrevlett.78.2690

CrossRef Full Text | Google Scholar

60. Crooks GE. Entropy production fluctuation theorem and the nonequilibrium work relation for free energy differences. Phys Rev E (1999) 60:2721–6. doi:10.1103/physreve.60.2721

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Seifert U. Stochastic thermodynamics, fluctuation theorems and molecular machines. Rep Prog Phys (2012) 75:126001. doi:10.1088/0034-4885/75/12/126001

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Caspers J, Krüger M. Identities for nonlinear memory kernels. Phys Rev E (2025) 112:024124. doi:10.1103/pzs6-l3ws

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Maes C. Local detailed balance. Scipost Phys Lect Notes (2021) 32. doi:10.21468/scipostphyslectnotes.32

CrossRef Full Text | Google Scholar

64. Cates ME, Candau SJ. Statics and dynamics of worm-like surfactant micelles. J Phys Condens Matter (1990) 2:6869–92. doi:10.1088/0953-8984/2/33/001

CrossRef Full Text | Google Scholar

65. Gomez-Solano JR, Bechinger C. Transient dynamics of a colloidal particle driven through a viscoelastic fluid. New J Phys (2015) 17:103032. doi:10.1088/1367-2630/17/10/103032

CrossRef Full Text | Google Scholar

Keywords: fluctuation–dissipation theorem, nonequilibrium cumulants, Brownian probe particle, optical potential, nonlinear fluid, non-Markovian fluid, worm-like micelles, micellar fluid

Citation: Caspers J, Krishna Kumar K, Bechinger C and Krüger M (2025) Equilibrium trajectories quantify second-order violations of the fluctuation–dissipation theorem without the need for a model. Front. Phys. 13:1667224. doi: 10.3389/fphy.2025.1667224

Received: 16 July 2025; Accepted: 18 September 2025;
Published: 21 October 2025.

Edited by:

Andre P. Vieira, University of São Paulo, Brazil

Reviewed by:

Saravana Prakash Thirumuruganandham, SIT Health, Ecuador
Pedro Harunari, University of Luxembourg, Luxembourg

Copyright © 2025 Caspers, Krishna Kumar, Bechinger and Krüger. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Juliana Caspers, ai5jYXNwZXJzQHRoZW9yaWUucGh5c2lrLnVuaS1nb2V0dGluZ2VuLmRl; Matthias Krüger, bWF0dGhpYXMua3J1Z2VyQHVuaS1nb2V0dGluZ2VuLmRl

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.