Violation of the fluctuation–dissipation theorem in glassy systems: basic notions and the numerical evidence

The Thermodynamics of Mind

To not only survive, but also thrive, the brain must efficiently orchestrate distributed computation across space and time. This requires hierarchical organisation facilitating fast information transfer and processing at the lowest possible metabolic cost. Quantifying brain hierarchy is difficult but can be estimated from the asymmetry of information flow. Thermodynamics has successfully characterised hierarchy in many other complex systems. Here, we propose the 'Thermodynamics of Mind' framework as a natural way to quantify hierarchical brain orchestration and its underlying mechanisms. This has already provided novel insights into the orchestration of hierarchy in brain states including movie watching, where the hierarchy of the brain is flatter than during rest. Overall, this framework holds great promise for revealing the orchestration of cognition.

Fluctuation-response relations for integrate-and-fire models with an absolute refractory period

We study the problem of relating the spontaneous fluctuations of a stochastic integrate-and-fire (IF) model to the response of the instantaneous firing rate to time-dependent stimulation if the IF model is endowed with a non-vanishing refractory period and a finite (stereotypical) spike shape. This seemingly harmless addition to the model is shown to complicate the analysis put forward by Lindner Phys. Rev. Lett. (2022), i.e., the incorporation of the reset into the model equation, the Rice-like averaging of the stochastic differential equation, and the application of the Furutsu-Novikov theorem. We derive a still exact (although more complicated) fluctuation-response relation (FRR) for an IF model with refractory state and a white Gaussian background noise. We also briefly discuss an approximation for the case of a colored Gaussian noise and conclude with a summary and outlook on open problems.

Violations of the fluctuation-dissipation theorem reveal distinct nonequilibrium dynamics of brain states

The brain is a nonequilibrium system whose dynamics change in different brain states, such as wakefulness and deep sleep. Thermodynamics provides the tools for revealing these nonequilibrium dynamics. We used violations of the fluctuation-dissipation theorem to describe the hierarchy of nonequilibrium dynamics associated with different brain states. Together with a whole-brain model fitted to empirical human neuroimaging data, and deriving the appropriate analytical expressions, we were able to capture the deviation from equilibrium in different brain states that arises from asymmetric interactions and hierarchical organization.

Entropy production and collective excitations of crystals out of equilibrium: The concept of entropons

We study the collective vibrational excitations of crystals under out-of-equilibrium steady conditions that give rise to entropy production. Their excitation spectrum comprises equilibriumlike phonons of thermal origin and additional collective excitations called entropons because each of them represents a mode of spectral entropy production. Entropons coexist with phonons and dominate them when the system is far from equilibrium while they are negligible in near-equilibrium regimes. The concept of entropons has been recently introduced and verified in a special case of crystals formed by self-propelled particles. Here we show that entropons exist in a broader class of active crystals that are intrinsically out of equilibrium and characterized by the lack of detailed balance. After a general derivation, several explicit examples are discussed, including crystals consisting of particles with alignment interactions and frictional contact forces.

Nonintegral form of the reciprocal relation associated with violation of the fluctuation response relation

We extend Onsager's reciprocal relation to systems in a nonequilibrium steady state. While Onsager's reciprocal relation concerns the kinetic (Onsager) coefficient, the extended reciprocal relation concerns violation of the fluctuation response relation (FRR) for mechanical and thermal perturbations. This extended relation holds at each frequency when the extent of the FRR violation is expressed in a frequency domain. This nonintegral form distinguishes the extended relation from previous relations expressed by integration over a frequency. To obtain this relation, we consider one-particle one-dimensional systems described by an overdamped Langevin equation with a force driving the system away from equilibrium. We assume a special property of the potential in the system. From this Langevin equation, we obtain the Fokker-Planck (FP) equation describing the time evolution of the distribution function of the particle. Using the FP equation, we calculate the responses of the particle velocity and heat current by applying time-dependent perturbations of the driving force and temperature. We express the extent of the FRR violation in terms of these responses with time correlation functions and expand them in powers of the FP operator. This reciprocal relation is valid far from equilibrium. One can also confirm this reciprocal relation through experiments with systems such as colloidal suspensions because the FRR violation can be experimentally observed.

Reorganization energy of electron transfer

The theory of electron transfer reactions establishes the conceptual foundation for redox solution chemistry, electrochemistry, and bioenergetics. Electron and proton transfer across the cellular membrane provide all energy of life gained through natural photosynthesis and mitochondrial respiration. Rates of biological charge transfer set kinetic bottlenecks for biological energy storage. The main system-specific parameter determining the activation barrier for a single electron-transfer hop is the reorganization energy of the medium. Both harvesting of light energy in natural and artificial photosynthesis and efficient electron transport in biological energy chains require reduction of the reorganization energy to allow fast transitions. This review article discusses mechanisms by which small values of the reorganization energy are achieved in protein electron transfer and how similar mechanisms can operate in other media, such as nonpolar and ionic liquids. One of the major mechanisms of reorganization energy reduction is through non-Gibbsian (nonergodic) sampling of the medium configurations on the reaction time. A number of alternative mechanisms, such as electrowetting of active sites of proteins, give rise to non-parabolic free energy surfaces of electron transfer. These mechanisms, and nonequilibrium population of donor-acceptor vibrations, lead to a universal phenomenology of separation between the Stokes shift and variance reorganization energies of electron transfer.

Nonlinear response theory for Markov processes. IV. The asymmetric double-well potential model revisited

The dielectric response of noninteracting dipoles is discussed in the framework of the classical model of stochastic reorientations in an asymmetric double-well potential (ADWP). In the nonlinear regime, this model exhibits some pecularities in the static response. We find that the saturation behavior of the symmetric double-well potential model does not follow the Langevin function and only in the linear regime are the standard results recovered. If a finite asymmetry is assumed, then the nonlinear susceptibilities are found to change the sign at a number of characteristic temperatures that depend on the magnitude of the asymmetry, as has been observed earlier for the third-order and fifth-order responses. If the kinetics of the barrier crossing in the ADWP model is described as a two-state model, then we can give analytical expressions for the values of the characteristic temperatures. The results for the response obtained from a (numerical) solution of the Fokker-Planck equation for the Brownian motion in a model ADWP behaves very similarly to the two-state model for high barriers. For small barriers no clear-cut timescale separation between the barrier crossing process and the intrawell relaxation exists and the model exhibits a number of timescales. In this case, the frequency-dependent linear susceptibility at low temperatures is dominated by the fast intrawell transitions and at higher temperatures by the barrier crossing kinetics. We find that for nonlinear susceptibilities the latter process appears to be more important and the intrawell transitions play only a role at the lowest temperatures.

From generalized Langevin stochastic dynamics to anomalous diffusion

Scaling methods are fundamental in all branches of physics. In stochastic process, we usually try to describe the long time behavior of a given time correlation function. In this work we investigate a scaling method for anomalous diffusion in systems with memory that produces good results for long and intermediate times. We will initially present a generalization of the diffusion exponent. Then, we present an asymptotic method to obtain an analytical expression for the diffusion coefficient by introducing a time scale factor λ(t). We found an exact expression for the function λ(t), which allows us to describe the diffusive process. For large times, λ(t) becomes a universal parameter determined by the diffusion exponent. In turn, the analytical results are then compared to the numerical results, with a good matching. Then, we'll show the practical effects of scaling. An important first result is that λ(t) quickly converges to a constant. Another very important point was the classification of new forms of diffusion due to the generalized exponent. In previous works, we verified the existence of ergodic ballistic diffusion with diffusion exponent α=2^{-}. Here, we verify the existence of the nonergodic ballistic diffusion type with the obtainment of the diffusion coefficient α=2. Finally, we show that the scaling works. This method is general and can be applied to various types of stochastic problems.

Conformational dynamics modulating electron transfer

Diffusional dynamics of the donor-acceptor distance are responsible for the appearance of a new time scale of diffusion over the distance of electronic tunneling in electron-transfer reactions. The distance dynamics compete with the medium polarization dynamics in the dynamics-controlled electron-transfer kinetics. The pre-exponential factor of the electron-transfer rate constant switches, at the crossover distance, between a distance-independent, dynamics-controlled plateau and exponential distance decay. The crossover between two regimes is controlled by an effective relaxation time slowed down by a factor exponentially depending on the variance of the donor-acceptor displacement. Flexible donor-acceptor complexes must show a greater tendency for dynamics-controlled electron transfer. Energy chains based on electron transport are best designed by placing the redox cofactors near the crossover distance.

Translation-rotation decoupling of tracers reflects medium-range crystalline order in two-dimensional colloid glasses

The dynamic heterogeneity and the translation-rotation decoupling are the dynamic signatures of glasses and supercooled liquids. Whether and how the dynamic heterogeneity would relate to the local structure of glasses has been a puzzle for decades. In this work we perform molecular dynamics simulations for tracers in both two-dimensional polydisperse colloids (2DPC) and two-dimensional binary colloids (2DBC). In 2DPC glasses, hexatic local structures develop at low enough temperatures and grow quickly along with the dynamic correlation length of the 2DPC, which is well known as the medium-range crystalline order (MRCO). In 2DBC glasses, on the other hand, any explicit local structure has not been reported to grow significantly with the dynamic correlation length at low temperatures. We introduce two different types of tracers into colloidal systems: A diamond tracer that resembles the MRCO of 2DPC glasses and a square tracer that is dissimilar to any local structure of glasses. The translation-rotation decoupling of the diamond tracer in 2DPC glasses is much more significant than that of the square tracer in the same 2DPC glasses. On the other hand, such a tracer shape-dependence of the decoupling is not observed in 2DBC glasses where the local hexatic structure does not develop significantly. We introduce a shape-dependency parameter of the decoupling and find that the shape-dependency parameter grows along with the dynamic correlation length in 2DPC glasses but not in 2DBC glasses. This illustrates that the dynamic heterogeneity and the translation-rotation decoupling of tracers could reveal the local structure that develops in glasses.

Thermodynamics of fluctuations based on time-and-space averages

We develop nonequilibrium theory by using averages in time and space as a generalized way to upscale thermodynamics in nonergodic systems. The approach offers a classical perspective on the energy dynamics in fluctuating systems. The rate of entropy production is shown to be explicitly scale dependent when considered in this context. We show that while any stationary process can be represented as having zero entropy production, second law constraints due to the Clausius theorem are preserved due to the fact that heat and work are related based on conservation of energy. As a demonstration, we consider the energy dynamics for the Carnot cycle and for Maxwell's demon. We then consider nonstationary processes, applying time-and-space averages to characterize nonergodic effects in heterogeneous systems where energy barriers such as compositional gradients are present. We show that the derived theory can be used to understand the origins of anomalous diffusion phenomena in systems where Fick's law applies at small length scales, but not at large length scales. We further characterize fluctuations in capillary-dominated systems, which are nonstationary due to the irreversibility of cooperative events.

Ewald sum corrections in simulations of ion and dipole solvation and electron transfer

Periodic boundary conditions and Ewald sums used in standard simulation protocols require finite-size corrections when the total charge of the simulated system is nonzero. Corrections for ion solvation were introduced by Hummer, Pratt, and García, [J. Chem. Phys. 107, 9275 (1997)]. The latter approach is extended here to derive finite-size correction for the Stokes-shift and reorganization energy applied to electron-transfer reactions. The same correction term, scaling inversely with the box size, adds to the reorganization energy from the energy-gap variance but is subtracted from the reorganization energy calculated from the Stokes shift. Finite-size corrections thus widen the gap between these two quantities, which were recently found to diverge for protein electron transfer. Corrections to the free energy of dipole solvation and the variance of the electric field scale as m²/L³ with the solute dipole m and the box size L.

Fluctuation-dissipation relations far from equilibrium: a case study

Fluctuation-dissipation relations or "theorems" (FDTs) are fundamental for statistical physics and can be rigorously derived for equilibrium systems. Their applicability to non-equilibrium systems is, however, debated. Here, we simulate an active microrheology experiment, in which a spherical colloid is pulled with a constant external force through a fluid, creating near-equilibrium and far-from-equilibrium systems. We characterize the structural and dynamical properties of these systems, and reconstruct an effective generalized Langevin equation (GLE) for the colloid dynamics. Specifically, we test the validity of two FDTs: The first FDT relates the non-equilibrium response of a system to equilibrium correlation functions, and the second FDT relates the memory friction kernel in the GLE to the stochastic force. We find that the validity of the first FDT depends strongly on the strength of the external driving: it is fulfilled close to equilibrium and breaks down far from it. In contrast, we observe that the second FDT is always fulfilled. We provide a mathematical argument why this generally holds for memory kernels reconstructed from a deterministic Volterra equation for correlation functions, even for non-stationary non-equilibrium systems. Motivated by the Mori-Zwanzig formalism, we therefore suggest to impose an orthogonality constraint on the stochastic force, which is in fact equivalent to the validity of this Volterra equation. Such GLEs automatically satisfy the second FDT and are unique, which is desirable when using GLEs for coarse-grained modeling.

Unifying fluctuation-dissipation temperatures of slow-evolving nonequilibrium systems from the perspective of inherent structures

For nonequilibrium systems, how to define temperature is one of the key and difficult issues to solve. Although effective temperatures have been proposed and studied to this end, it still remains elusive what they actually are. Here, we focus on the fluctuation-dissipation temperatures and report that such effective temperatures of slow-evolving systems represent characteristic temperatures of their equilibrium counterparts. By calculating the fluctuation-dissipation relation of inherent structures, we obtain a temperature-like quantity T IS For monocomponent crystal-formers, T IS agrees well with the crystallization temperature T c, while it matches with the onset temperature T on for glass-formers. It also agrees with effective temperatures of typical nonequilibrium systems, such as aging glasses, quasi-static shear flows, and quasi-static self-propelled flows. From the unique perspective of inherent structures, our study reveals the nature of effective temperatures and the underlying connections between nonequilibrium and equilibrium systems and confirms the equivalence between T on and T c.

The Einstein effective temperature can predict the tagged active particle density

We derive a distribution function for the position of a tagged active particle in a slowly varying in space external potential, in a system of interacting active particles. The tagged particle distribution has the form of the Boltzmann distribution but with an effective temperature that replaces the temperature of the heat bath. We show that the effective temperature that enters the tagged particle distribution is the same as the effective temperature defined through the Einstein relation, i.e., it is equal to the ratio of the self-diffusion and tagged particle mobility coefficients. This result shows that this effective temperature, which is defined through a fluctuation-dissipation ratio, is relevant beyond the linear response regime. We verify our theoretical findings through computer simulations. Our theory fails when an additional large length scale appears in our active system. In the system we simulated, this length scale is associated with long-wavelength density fluctuations that emerge upon approaching motility-induced phase separation.

Effects of dynein inhibitor on the number of motor proteins transporting synaptic cargos

Synaptic cargo transport by kinesin and dynein in hippocampal neurons was investigated by noninvasively measuring the transport force based on nonequilibrium statistical mechanics. Although direct physical measurements such as force measurement using optical tweezers are difficult in an intracellular environment, the noninvasive estimations enabled enumerating force-producing units (FPUs) carrying a cargo comprising the motor proteins generating force. The number of FPUs served as a barometer for stable and long-distance transport by multiple motors, which was then used to quantify the extent of damage to axonal transport by dynarrestin, a dynein inhibitor. We found that dynarrestin decreased the FPU for retrograde transport more than for anterograde transport. This result indicates the applicability of the noninvasive force measurements. In the future, these measurements may be used to quantify damage to axonal transport resulting from neuronal diseases, including Alzheimer's, Parkinson's, and Huntington's diseases.

Interconnections between networks acting like an external field in a first-order percolation transition

Many interdependent, real-world infrastructures involve interconnections between different communities or cities. Here we show how the effects of such interconnections can be described as an external field for interdependent networks experiencing a first-order percolation transition. We find that the critical exponents γ and δ, related to the external field, can also be defined for first-order transitions but that they have different values than those found for second-order transitions. Surprisingly, we find that both sets of different exponents (for first and second order) can even be found within a single model of interdependent networks, depending on the dependency coupling strength. Nevertheless, in both cases both sets satisfy Widom's identity, δ-1=γ/β, which further supports the validity of their definitions. Furthermore, we find that both Erdős-Rényi and scale-free networks have the same values of the exponents in the first-order regime, implying that these models are in the same universality class. In addition, we find that in k-core percolation the values of the critical exponents related to the field are the same as for interdependent networks, suggesting that these systems also belong to the same universality class.

Electron transfer in nonpolar media

Electron transfer in nonpolar media violates the temperature scaling predicted by the fluctuation–dissipation theorem.

On the fluctuation-dissipation relation in non-equilibrium and non-Hamiltonian systems

We review generalized fluctuation-dissipation relations, which are valid under general conditions even in "nonstandard systems," e.g., out of equilibrium and/or without a Hamiltonian structure. The response functions can be expressed in terms of suitable correlation functions computed in the unperturbed dynamics. In these relations, typically, one has nontrivial contributions due to the form of the stationary probability distribution; such terms take into account the interaction among the relevant degrees of freedom in the system. We illustrate the general formalism with some examples in nonstandard cases, including driven granular media, systems with a multiscale structure, active matter, and systems showing anomalous diffusion.

Configurational entropy of glass-forming liquids

The configurational entropy is one of the most important thermodynamic quantities characterizing supercooled liquids approaching the glass transition. Despite decades of experimental, theoretical, and computational investigation, a widely accepted definition of the configurational entropy is missing, its quantitative characterization remains fraught with difficulties, misconceptions, and paradoxes, and its physical relevance is vividly debated. Motivated by recent computational progress, we offer a pedagogical perspective on the configurational entropy in glass-forming liquids. We first explain why the configurational entropy has become a key quantity to describe glassy materials, from early empirical observations to modern theoretical treatments. We explain why practical measurements necessarily require approximations that make its physical interpretation delicate. We then demonstrate that computer simulations have become an invaluable tool to obtain precise, nonambiguous, and experimentally relevant measurements of the configurational entropy. We describe a panel of available computational tools, offering for each method a critical discussion. This perspective should be useful to both experimentalists and theoreticians interested in glassy materials and complex systems.

Investigation of multiple-dynein transport of melanosomes by non-invasive force measurement using fluctuation unit χ

Pigment organelles known as melanosomes disperse or aggregate in a melanophore in response to hormones. These movements are mediated by the microtubule motors kinesin-2 and cytoplasmic dynein. However, the force generation mechanism of dynein, unlike that of kinesin, is not well understood. In this study, to address this issue, we investigated the dynein-mediated aggregation of melanosomes in zebrafish melanophores. We applied the fluctuation theorem of non-equilibrium statistical mechanics to estimate forces acting on melanosomes during transport by dynein, given that the energy of a system is related to its fluctuation. Our results demonstrate that multiple force-producing units cooperatively transport a single melanosome. Since the force is generated by dynein, this suggests that multiple dyneins carry a single melanosome. Cooperative transport has been reported for other organelles; thus, multiple-motor transport may be a universal mechanism for moving organelles within the cell.

Application of the fluctuation theorem to motor proteins: from F1-ATPase to axonal cargo transport by kinesin and dynein

The fluctuation theorem is a representative theorem in non-equilibrium statistical physics actively studied in the 1990s. Relating to entropy production in non-equilibrium states, the theorem has been used to estimate the driving power of motor proteins from fluctuation in their motion. In this review, usage of the fluctuation theorem in experiments on motor proteins is illustrated for biologists, especially those who study mechanobiology, in which force measurement is a central issue. We first introduce the application of the fluctuation theorem in measuring the rotary torque of the rotary motor protein F1-ATPase. Next, as an extension of this application, a recent trial estimating the force generated during cargo transport in vivo by the microtubule motors kinesin and dynein is introduced. Elucidation of the physical mechanism of such transport is important, especially for neurons, in which deficits in cargo transport are deeply related to neuronal diseases. Finally, perspectives on the fluctuation theorem as a new technique in the field of neuroscience are discussed.

Two-tag correlations and nonequilibrium fluctuation-response relation in ageing single-file diffusion

Spatiotemporally correlated motions of interacting Brownian particles, confined in a narrow channel of infinite length, are studied in terms of statistical quantities involving two particles. A theoretical framework that allows analytical calculation of two-tag correlations is presented on the basis of the Dean-Kawasaki equation describing density fluctuations in colloidal systems. In the equilibrium case, the time-dependent Einstein relation holds between the two-tag displacement correlation and the response function corresponding to it, which is a manifestation of the fluctuation-dissipation theorem for the correlation of density fluctuations. While the standard procedure of closure approximation for nonlinear density fluctuations is known to be obstructed by inconsistency with the fluctuation-dissipation theorem, this difficulty is naturally avoided by switching from the standard Fourier representation of the density field to the label-based Fourier representation of the vacancy field. In the case of ageing dynamics started from equidistant lattice configuration, the time-dependent Einstein relation is violated, as the two-tag correlation depends on the waiting time for equilibration while the response function is not sensitive to it. Within linear approximation, however, there is a simple relation between the density (or vacancy) fluctuations and the corresponding response function, which is valid even if the system is out of equilibrium. This non-equilibrium fluctuation-response relation can be extended to the case of nonlinear fluctuations by means of closure approximation for the vacancy field.

Activated dynamics: An intermediate model between the random energy model and the p-spin model

To study the activated dynamics of mean-field glasses, which takes place on times of order exp(N), where N is the system size, we introduce a new model, the correlated random energy model (CREM), that allows for a smooth interpolation between the REM and the p-spin models. We study numerically and analytically the CREM in the intermediate regime between REM and p-spin. We fully characterize its energy landscape, which is like a golf course but, at variance with the REM, has metabasins (or holes) containing several configurations. We find that an effective description for the dynamics, in terms of traps, emerges, provided that one identifies metabasins in the CREM with configurations in the trap model.

Isomorph theory of physical aging

This paper derives and discusses the configuration-space Langevin equation describing a physically aging R-simple system and the corresponding Smoluchowski equation. Externally controlled thermodynamic variables like temperature, density, and pressure enter the description via the single parameter T_s/T, in which T is the bath temperature and T_s is the "systemic" temperature defined at any time t as the thermodynamic equilibrium temperature of the state point with density ρ(t) and potential energy U(t). In equilibrium, T_s ≅ T with fluctuations that vanish in the thermodynamic limit. In contrast to Tool's fictive temperature and other effective temperatures in glass science, the systemic temperature is defined for any configuration with a well-defined density, even if it is not close to equilibrium. Density and systemic temperature define an aging phase diagram, in which the aging system traces out a curve. Predictions are discussed for aging following various density-temperature and pressure-temperature jumps from one equilibrium state to another, as well as for a few other scenarios. The proposed theory implies that R-simple glass-forming liquids are characterized by the dynamic Prigogine-Defay ratio being equal to unity.

Derivation of the spin-glass order parameter from stochastic thermodynamics

A fluctuation relation is derived to extract the order parameter function q(x) in weakly ergodic systems. The relation is based on measuring and classifying entropy production fluctuations according to the value of the overlap q between configurations. For a fixed value of q, entropy production fluctuations are Gaussian distributed allowing us to derive the quasi-FDT so characteristic of aging systems. The theory is validated by extracting the q(x) in various types of glassy models. It might be generally applicable to other nonequilibrium systems and experimental small systems.

Ergodicity breaking of iron displacement in heme proteins

We present a model of the dynamical transition of atomic displacements in proteins. Increased mean-square displacement at higher temperatures is caused by the softening of the force constant for atomic/molecular displacements by electrostatic and van der Waals forces from the protein-water thermal bath. Displacement softening passes through a nonergodic dynamical transition when the relaxation time of the force-force correlation function enters, with increasing temperature, the instrumental observation window. Two crossover temperatures are identified. The lower crossover, presently connected to the glass transition, is related to the dynamical unfreezing of rotations of water molecules within nanodomains polarized by charged surface residues of the protein. The higher crossover temperature, usually assigned to the dynamical transition, marks the onset of water translations. All crossovers are ergodicity breaking transitions depending on the corresponding observation windows. Allowing stretched exponential relaxation of the protein-water thermal bath significantly improves the theory-experiment agreement when applied to solid protein samples studied by Mössbauer spectroscopy.

Conditional 1/f^{α} noise: From single molecules to macroscopic measurement

We demonstrate that the measurement of 1/f^{α} noise at the single molecule or nano-object limit is remarkably distinct from the macroscopic measurement over a large sample. The single-particle measurements yield a conditional time-dependent spectrum. However, the number of units fluctuating on the time scale of the experiment is increasing in such a way that the macroscopic measurements appear perfectly stationary. The single-particle power spectrum is a conditional spectrum, in the sense that we must make a distinction between idler and nonidler units on the time scale of the experiment. We demonstrate our results based on stochastic and deterministic models, in particular the well-known approach of superimposed Lorentzians, the blinking quantum dot model, and deterministic dynamics generated by a nonlinear mapping. Our results show that the 1/f^{α} spectrum is inherently nonstationary even if the macroscopic measurement completely obscures the underlying time dependence of the phenomena.

Nonlinear response theory for Markov processes. II. Fifth-order response functions

The nonlinear response of stochastic models obeying a master equation is calculated up to fifth order in the external field, thus extending the third-order results obtained earlier [G. Diezemann, Phys. Rev. E 85, 051502 (2012)PLEEE81539-375510.1103/PhysRevE.85.051502]. For sinusoidal fields the 5ω component of the susceptibility is computed for the model of dipole reorientations in an asymmetric double well potential and for a trap model with a Gaussian density of states. For most realizations of the models a hump is found in the higher-order susceptibilities. In particular, for the asymmetric double well potential model there are two characteristic temperature regimes showing the occurrence of such a hump as compared to a single characteristic regime in the case of the third-order response. In the case of the trap model the results strongly depend on the variable coupled to the field. As for the third-order response, the low-frequency limit of the susceptibility plays a crucial role with respect to the occurrence of a hump. The findings are discussed in light of recent experimental results obtained for supercooled liquids. The differences found for the third-order and the fifth-order response indicate that nonlinear response functions might serve as a powerful tool to discriminate among the large number of existing models for glassy relaxation.

Long-term aging behaviors in a model soft colloidal system

Colloidal and molecular systems share similar behaviors near to the glass transition volume fraction or temperature. Here, aging behaviors after volume fraction up-jump (induced by performing temperature down-jumps) conditions for a PS-PNIPAM/AA soft colloidal system were investigated using light scattering (diffusing wave spectroscopy, DWS). Both aging responses and equilibrium dynamics were investigated. For the aging responses, long-term experiments (100 000 s) were performed, and both equilibrium and non-equilibrium behaviors of the system were obtained. In the equilibrium state, as effective volume fraction increases (or temperature decreases), the colloidal dispersion displays a transition from the liquid to a glassy state. The equilibrium α-relaxation dynamics strongly depend on both the effective volume fraction and the initial mass concentration for the studied colloidal systems. Compared with prior results from our lab [X. Di, X. Peng and G. B. McKenna, J. Chem. Phys., 2014, 140, 054903], the effective volume fractions investigated spanned a wider range, to deeper into the glassy domain. The results show that the α-relaxation time τ_α of the samples aged into equilibrium deviate from the classical Vogel-Fulcher-Tammann (VFT)-type expectations and the super-Arrhenius signature disappears above the glass transition volume fraction. The non-equilibrium aging response shows that the time for the structural evolution into equilibrium and the α-relaxation time are decoupled. The DWS investigation of the aging behavior after different volume fraction jumps reveals a different non-equilibrium or aging behavior for the considered colloidal systems compared with either molecular glasses or the macroscopic rheology of a similar colloidal dispersions.

Aging Wiener-Khinchin theorem and critical exponents of 1/f^{β} noise

The power spectrum of a stationary process may be calculated in terms of the autocorrelation function using the Wiener-Khinchin theorem. We here generalize the Wiener-Khinchin theorem for nonstationary processes and introduce a time-dependent power spectrum 〈S_{t_{m}}(ω)〉 where t_{m} is the measurement time. For processes with an aging autocorrelation function of the form 〈I(t)I(t+τ)〉=t^{Υ}ϕ_{EA}(τ/t), where ϕ_{EA}(x) is a nonanalytic function when x is small, we find aging 1/f^{β} noise. Aging 1/f^{β} noise is characterized by five critical exponents. We derive the relations between the scaled autocorrelation function and these exponents. We show that our definition of the time-dependent spectrum retains its interpretation as a density of Fourier modes and discuss the relation to the apparent infrared divergence of 1/f^{β} noise. We illustrate our results for blinking-quantum-dot models, single-file diffusion, and Brownian motion in a logarithmic potential.

Nonlinear Response of Inertial Tracers in Steady Laminar Flows: Differential and Absolute Negative Mobility

We study the mobility and the diffusion coefficient of an inertial tracer advected by a two-dimensional incompressible laminar flow, in the presence of thermal noise and under the action of an external force. We show, with extensive numerical simulations, that the force-velocity relation for the tracer, in the nonlinear regime, displays complex and rich behaviors, including negative differential and absolute mobility. These effects rely upon a subtle coupling between inertia and applied force that induces the tracer to persist in particular regions of phase space with a velocity opposite to the force. The relevance of this coupling is revisited in the framework of nonequilibrium response theory, applying a generalized Einstein relation to our system. The possibility of experimental observation of these results is also discussed.

Polarizability of the active site of cytochrome c reduces the activation barrier for electron transfer

Enzymes in biology’s energy chains operate with low energy input distributed through multiple electron transfer steps between protein active sites. The general challenge of biological design is how to lower the activation barrier without sacrificing a large negative reaction free energy. We show that this goal is achieved through a large polarizability of the active site. It is polarized by allowing a large number of excited states, which are populated quantum mechanically by electrostatic fluctuations of the protein and hydration water shells. This perspective is achieved by extensive mixed quantum mechanical/molecular dynamics simulations of the half reaction of reduction of cytochrome c. The barrier for electron transfer is consistently lowered by increasing the number of excited states included in the Hamiltonian of the active site diagonalized along the classical trajectory. We suggest that molecular polarizability, in addition to much studied electrostatics of permanent charges, is a key parameter to consider in order to understand how enzymes work.

Communication: Microsecond dynamics of the protein and water affect electron transfer in a bacterial bc(1) complex

Cross-membrane electron transport between cofactors localized in proteins of mitochondrial respiration and bacterial photosynthesis is the source of all biological energy. The statistics and dynamics of nuclear fluctuations in these protein/membrane/water heterogeneous systems are critical for their energetic efficiency. The results of 13 μs of atomistic molecular dynamics simulations of the membrane-bound bc1 bacterial complex are analyzed here. The reaction is affected by a broad spectrum of nuclear modes, with the slowest dynamics in the range of time-scales ∼0.1-1.6 μs contributing half of the reaction reorganization energy. Two reorganization energies are required to describe protein electron transfer due to dynamical arrest of protein conformations on the observation window. This mechanistic distinction allows significant lowering of activation barriers for reactions in proteins.

Protein electron transfer: is biology (thermo)dynamic?

Simple physical mechanisms are behind the flow of energy in all forms of life. Energy comes to living systems through electrons occupying high-energy states, either from food (respiratory chains) or from light (photosynthesis). This energy is transformed into the cross-membrane proton-motive force that eventually drives all biochemistry of the cell. Life's ability to transfer electrons over large distances with nearly zero loss of free energy is puzzling and has not been accomplished in synthetic systems. The focus of this review is on how this energetic efficiency is realized. General physical mechanisms and interactions that allow proteins to fold into compact water-soluble structures are also responsible for a rugged landscape of energy states and a broad distribution of relaxation times. Specific to a protein as a fluctuating thermal bath is the protein-water interface, which is heterogeneous both dynamically and structurally. The spectrum of interfacial fluctuations is a consequence of protein's elastic flexibility combined with a high density of surface charges polarizing water dipoles into surface nanodomains. Electrostatics is critical to the protein function and the relevant questions are: (i) What is the spectrum of interfacial electrostatic fluctuations? (ii) Does the interfacial biological water produce electrostatic signatures specific to proteins? (iii) How is protein-mediated chemistry affected by electrostatics? These questions connect the fluctuation spectrum to the dynamical control of chemical reactivity, i.e. the dependence of the activation free energy of the reaction on the dynamics of the bath. Ergodicity is often broken in protein-driven reactions and thermodynamic free energies become irrelevant. Continuous ergodicity breaking in a dense spectrum of relaxation times requires using dynamically restricted ensembles to calculate statistical averages. When applied to the calculation of the rates, this formalism leads to the nonergodic activated kinetics, which extends the transition-state theory to dynamically dispersive media. Releasing the grip of thermodynamics in kinetic calculations through nonergodicity provides the mechanism for an efficient optimization between reaction rates and the spectrum of relaxation times of the protein-water thermal bath. Bath dynamics, it appears, play as important role as the free energy in optimizing biology's performance.

Linear response and modified fluctuation-dissipation relation in random potential

In this work, a physical system described by the Hamiltonian H(ω)=H(0)+V(ω)(t) consisting of a solvable model H(0) and external random and time-dependent potential V(ω)(t) is investigated. Under the conditions in which, for each realization, the potential changes smoothly so that the evolution of the system follows the Schrödinger dynamics, and that the average external potential with respect to all realizations is constant in time, an adjusted equilibrium state can be defined as a reference state and the mean dynamics can be derived from taking the average of the equation with respect to the configuration parameter ω. It provides extra contributions from the deviations of the Hamiltonian and evolves the state along the time by the Heisenberg and Liouville-von Neumann equations. Consequently, the Kubo formula and the fluctuation-dissipation relation (FDR) are modified in the sense that the contribution from the information of randomness and memory effects from the time dependence is also present. The modified Kubo formula now has a contribution from two terms. The first term is an antisymmetric cross correlation between two observables measured by a probe as expected, and the latter term is an accumulation of the propagation of the effects from the randomness. When the considered system is in the adjusted equilibrium state at the time the measurement probe interacts, the latter contribution vanishes, and the standard FDR is recovered.

Aging Wiener-Khinchin Theorem

The Wiener-Khinchin theorem shows how the power spectrum of a stationary random signal I(t) is related to its correlation function ⟨I(t)I(t+τ)⟩. We consider nonstationary processes with the widely observed aging correlation function ⟨I(t)I(t+τ)⟩∼t(γ)ϕ(EA)(τ/t) and relate it to the sample spectrum. We formulate two aging Wiener-Khinchin theorems relating the power spectrum to the time- and ensemble-averaged correlation functions, discussing briefly the advantages of each. When the scaling function ϕ(EA)(x) exhibits a nonanalytical behavior in the vicinity of its small argument we obtain the aging 1/f-type of spectrum. We demonstrate our results with three examples: blinking quantum dots, single-file diffusion, and Brownian motion in a logarithmic potential, showing that our approach is valid for a wide range of physical mechanisms.

Photosynthetic diode: electron transport rectification by wetting the quinone cofactor

We report 11 μs of molecular dynamics simulations of the electron-transfer reaction between primary and secondary quinone cofactors in the bacterial reaction center. The main question addressed here is the mechanistic reason for unidirectional electron transfer between chemically identical cofactors. We find that electron is trapped at the secondary quinone by wetting of the protein pocket following electron transfer on the time-scale shorter than the backward transition. This mechanism provides effective rectification of the electron transport, making the reaction center a molecular diode operating by cyclic charge-induced electrowetting.

Observable-dependence of the effective temperature in off-equilibrium diatomic molecular liquids

We discuss the observable-dependence of the effective temperature Teff, defined via the fluctuation-dissipation relation, of an out-of-equilibrium system composed by homonuclear dumbbell molecules. Teff is calculated by evaluating the fluctuation and the response for two observables associated, respectively, to translational and to rotational degrees of freedom, following a sudden temperature quench. We repeat our calculations for different dumbbell elongations ζ. At high elongations (ζ > 0.4), we find the same Teff for the two observables. At low elongations (ζ ⩽ 0.4), only for very deep quenches Teff coincides. The observable-dependence of Teff for low elongations and shallow quenches stresses the importance of a strong coupling between orientational and translational variables for a consistent definition of the effective temperature in glassy systems.

Functional significance of complex fluctuations in brain activity: from resting state to cognitive neuroscience

Behavioral studies have shown that human cognition is characterized by properties such as temporal scale invariance, heavy-tailed non-Gaussian distributions, and long-range correlations at long time scales, suggesting models of how (non observable) components of cognition interact. On the other hand, results from functional neuroimaging studies show that complex scaling and intermittency may be generic spatio-temporal properties of the brain at rest. Somehow surprisingly, though, hardly ever have the neural correlates of cognition been studied at time scales comparable to those at which cognition shows scaling properties. Here, we analyze the meanings of scaling properties and the significance of their task-related modulations for cognitive neuroscience. It is proposed that cognitive processes can be framed in terms of complex generic properties of brain activity at rest and, ultimately, of functional equations, limiting distributions, symmetries, and possibly universality classes characterizing them.

Nonequilibrium fluctuation-dissipation theorem and heat production

We use a relationship between response and correlation function in nonequilibrium systems to establish a connection between the heat production and the deviations from the equilibrium fluctuation-dissipation theorem. This scheme extends the Harada-Sasa formulation [Phys. Rev. Lett. 95, 130602 (2005)], obtained for Langevin equations in steady states, as it also holds for transient regimes and for discrete jump processes involving small entropic changes. Moreover, a general formulation includes two times and the new concepts of two-time work, kinetic energy, and of a two-time heat exchange that can be related to a nonequilibrium “effective temperature.” Numerical simulations of a chain of anharmonic oscillators and of a model for a molecular motor driven by adenosine triphosphate hydrolysis illustrate these points.

Nonequilibrium Brownian motion beyond the effective temperature

The condition of thermal equilibrium simplifies the theoretical treatment of fluctuations as found in the celebrated Einstein’s relation between mobility and diffusivity for Brownian motion. Several recent theories relax the hypothesis of thermal equilibrium resulting in at least two main scenarios. With well separated timescales, as in aging glassy systems, equilibrium Fluctuation-Dissipation Theorem applies at each scale with its own “effective” temperature. With mixed timescales, as for example in active or granular fluids or in turbulence, temperature is no more well-defined, the dynamical nature of fluctuations fully emerges and a Generalized Fluctuation-Dissipation Theorem (GFDT) applies. Here, we study experimentally the mixed timescale regime by studying fluctuations and linear response in the Brownian motion of a rotating intruder immersed in a vibro-fluidized granular medium. Increasing the packing fraction, the system is moved from a dilute single-timescale regime toward a denser multiple-timescale stage. Einstein’s relation holds in the former and is violated in the latter. The violation cannot be explained in terms of effective temperatures, while the GFDT is able to impute it to the emergence of a strong coupling between the intruder and the surrounding fluid. Direct experimental measurements confirm the development of spatial correlations in the system when the density is increased.

Characterization of the low-temperature properties of a simplified protein model

Prompted by results that showed that a simple protein model, the frustrated Gō model, appears to exhibit a transition reminiscent of the protein dynamical transition, we examine the validity of this model to describe the low-temperature properties of proteins. First, we examine equilibrium fluctuations. We calculate its incoherent neutron-scattering structure factor and show that it can be well described by a theory using the one-phonon approximation. By performing an inherent structure analysis, we assess the transitions among energy states at low temperatures. Then, we examine nonequilibrium fluctuations after a sudden cooling of the protein. We investigate the violation of the fluctuation-dissipation theorem in order to analyze the protein glass transition. We find that the effective temperature of the quenched protein deviates from the temperature of the thermostat, however it relaxes towards the actual temperature with an Arrhenius behavior as the waiting time increases. These results of the equilibrium and nonequilibrium studies converge to the conclusion that the apparent dynamical transition of this coarse-grained model cannot be attributed to a glassy behavior.

Nonequilibrium steady state and induced currents of a mesoscopically glassy system: interplay of resistor-network theory and Sinai physics

We introduce an explicit solution for the nonequilibrium steady state (NESS) of a ring that is coupled to a thermal bath, and is driven by an external hot source with log-wide distribution of couplings. Having time scales that stretch over several decades is similar to glassy systems. Consequently there is a wide range of driving intensities where the NESS is like that of a random walker in a biased Brownian landscape. We investigate the resulting statistics of the induced current I. For a single ring we discuss how sign of I fluctuates as the intensity of the driving is increased, while for an ensemble of rings we highlight the fingerprints of Sinai physics on the distribution of the absolute value of I.

Time asymmetry of the Kramers equation with nonlinear friction: fluctuation-dissipation relation and ratchet effect

We show by numerical simulations that the presence of nonlinear velocity-dependent friction forces can induce a finite net drift in the stochastic motion of a particle in contact with an equilibrium thermal bath and in an asymmetric periodic spatial potential. In particular, we study the Kramers equation for a particle subjected to Coulomb friction, namely a constant force acting in the direction opposite to the particle's velocity. We characterize the nonequilibrium irreversible dynamics by studying the generalized fluctuation-dissipation relation for this ratchet model driven by Coulomb friction.

Fluctuation relation for weakly ergodic aging systems

A fluctuation relation for aging systems is introduced and verified by extensive numerical simulations. It is based on the hypothesis of partial equilibration over phase-space regions in a scenario of entropy-driven relaxation. The relation provides a simple alternative method, amenable of experimental implementation, to measure replica symmetry breaking parameters in aging systems. The connection with the effective temperatures obtained from the fluctuation-dissipation theorem is discussed.