1
|
Kringelbach ML, Sanz Perl Y, Deco G. The Thermodynamics of Mind. Trends Cogn Sci 2024; 28:568-581. [PMID: 38677884 DOI: 10.1016/j.tics.2024.03.009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2024] [Revised: 03/16/2024] [Accepted: 03/18/2024] [Indexed: 04/29/2024]
Abstract
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.
Collapse
Affiliation(s)
- Morten L Kringelbach
- Centre for Eudaimonia and Human Flourishing, Linacre College, University of Oxford, Oxford, UK; Department of Psychiatry, University of Oxford, Oxford, UK; Center for Music in the Brain, Department of Clinical Medicine, Aarhus University, Aarhus, Denmark; International Centre for Flourishing, Universities of Oxford, Aarhus, and Pompeu Fabra, Oxford, UK.
| | - Yonatan Sanz Perl
- International Centre for Flourishing, Universities of Oxford, Aarhus, and Pompeu Fabra, Oxford, UK; Center for Brain and Cognition, Computational Neuroscience Group, Department of Information and Communication Technologies, Universitat Pompeu Fabra, Roc Boronat 138, Barcelona, Spain; Department of Physics, University of Buenos Aires, Buenos Aires, Argentina
| | - Gustavo Deco
- International Centre for Flourishing, Universities of Oxford, Aarhus, and Pompeu Fabra, Oxford, UK; Center for Brain and Cognition, Computational Neuroscience Group, Department of Information and Communication Technologies, Universitat Pompeu Fabra, Roc Boronat 138, Barcelona, Spain; Institució Catalana de la Recerca i Estudis Avançats (ICREA), Passeig Lluís Companys 23, Barcelona, Spain.
| |
Collapse
|
2
|
Puttkammer F, Lindner B. Fluctuation-response relations for integrate-and-fire models with an absolute refractory period. BIOLOGICAL CYBERNETICS 2024; 118:7-19. [PMID: 38261004 PMCID: PMC11068698 DOI: 10.1007/s00422-023-00982-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/09/2023] [Accepted: 12/11/2023] [Indexed: 01/24/2024]
Abstract
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.
Collapse
Affiliation(s)
- Friedrich Puttkammer
- Bernstein Center for Computational Neuroscience Berlin, Philippstr. 13, Haus 2, 10115, Berlin, Germany
- Physics Department of Humboldt University Berlin, Newtonstr. 15, 12489, Berlin, Germany
| | - Benjamin Lindner
- Bernstein Center for Computational Neuroscience Berlin, Philippstr. 13, Haus 2, 10115, Berlin, Germany.
- Physics Department of Humboldt University Berlin, Newtonstr. 15, 12489, Berlin, Germany.
| |
Collapse
|
3
|
Deco G, Lynn CW, Sanz Perl Y, Kringelbach ML. Violations of the fluctuation-dissipation theorem reveal distinct nonequilibrium dynamics of brain states. Phys Rev E 2023; 108:064410. [PMID: 38243472 DOI: 10.1103/physreve.108.064410] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2023] [Accepted: 10/16/2023] [Indexed: 01/21/2024]
Abstract
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.
Collapse
Affiliation(s)
- Gustavo Deco
- Center for Brain and Cognition, Computational Neuroscience Group, Universitat Pompeu Fabra, Roc Boronat 138, Barcelona 08018, Spain
- Institució Catalana de la Recerca i Estudis Avançats (ICREA), Passeig Lluís Companys 23, Barcelona 08010, Spain
| | - Christopher W Lynn
- Initiative for the Theoretical Sciences, Graduate Center, City University of New York, New York, New York 10016, USA and Joseph Henry Laboratories of Physics, Princeton University, Princeton, New Jersey 08544, USA
| | - Yonatan Sanz Perl
- Center for Brain and Cognition, Computational Neuroscience Group, Universitat Pompeu Fabra, Roc Boronat 138, Barcelona 08018, Spain
- Department of Physics, University of Buenos Aires, Buenos Aires 1428, Argentina and Paris Brain Institute (ICM), Paris 75013, France
| | - Morten L Kringelbach
- Centre for Eudaimonia and Human Flourishing, Linacre College, University of Oxford, Oxford OX3 9BX, United Kingdom; Department of Psychiatry, University of Oxford, Oxford OX3 7JX, United Kingdom; and Center for Music in the Brain, Department of Clinical Medicine, Aarhus University, Aarhus 8000, Denmark
| |
Collapse
|
4
|
Caprini L, Marini Bettolo Marconi U, Löwen H. Entropy production and collective excitations of crystals out of equilibrium: The concept of entropons. Phys Rev E 2023; 108:044603. [PMID: 37978682 DOI: 10.1103/physreve.108.044603] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2023] [Accepted: 09/21/2023] [Indexed: 11/19/2023]
Abstract
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.
Collapse
Affiliation(s)
- L Caprini
- Heinrich-Heine-Universität Düsseldorf, Institut für Theoretische Physik II: Weiche Materie, Universitätsstrasse, 40225 Düsseldorf, Germany
| | - U Marini Bettolo Marconi
- Physics Department, Scuola di Scienze e Tecnologie, Università di Camerino - via Madonna delle Carceri, 62032 Camerino, Italy
- Istituto Nazionale di Fisica Nucleare, Sezione di Perugia, Via A. Pascoli, 06123 Perugia, Italy
| | - H Löwen
- Heinrich-Heine-Universität Düsseldorf, Institut für Theoretische Physik II: Weiche Materie, Universitätsstrasse, 40225 Düsseldorf, Germany
| |
Collapse
|
5
|
Kasuga K, Yoshimori A. Nonintegral form of the reciprocal relation associated with violation of the fluctuation response relation. Phys Rev E 2023; 108:034109. [PMID: 37849163 DOI: 10.1103/physreve.108.034109] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2022] [Accepted: 08/15/2023] [Indexed: 10/19/2023]
Abstract
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.
Collapse
Affiliation(s)
- Kotaro Kasuga
- Department of Physics, Niigata University, Niigata 950-2181, Japan
| | - Akira Yoshimori
- Department of Physics, Niigata University, Niigata 950-2181, Japan
| |
Collapse
|
6
|
Abstract
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.
Collapse
Affiliation(s)
- Dmitry V Matyushov
- School of Molecular Sciences and Department of Physics, Arizona State University, PO Box 871504, Tempe, Arizona 85287-1504, USA.
| |
Collapse
|
7
|
Diezemann G. Nonlinear response theory for Markov processes. IV. The asymmetric double-well potential model revisited. Phys Rev E 2022; 106:064122. [PMID: 36671146 DOI: 10.1103/physreve.106.064122] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2022] [Accepted: 12/02/2022] [Indexed: 06/17/2023]
Abstract
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.
Collapse
Affiliation(s)
- Gregor Diezemann
- Department Chemie, Johannes Gutenberg Universität Mainz, 55128 Mainz, Germany
| |
Collapse
|
8
|
Ferreira RMS. From generalized Langevin stochastic dynamics to anomalous diffusion. Phys Rev E 2022; 106:054157. [PMID: 36559457 DOI: 10.1103/physreve.106.054157] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2022] [Accepted: 11/08/2022] [Indexed: 06/17/2023]
Abstract
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.
Collapse
Affiliation(s)
- Rogelma M S Ferreira
- Centro de Ciências Exatas e Tecnológicas, Universidade Federal do Recôncavo da Bahia, 44.380-000, Cruz das Almas, Bahia, Brasil
| |
Collapse
|
9
|
Matyushov DV. Conformational dynamics modulating electron transfer. J Chem Phys 2022; 157:095102. [DOI: 10.1063/5.0102707] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Dmitry V. Matyushov
- Departments of Physics and School of Molecular Sciences, Arizona State University, United States of America
| |
Collapse
|
10
|
|
11
|
Chun DJ, Oh Y, Sung BJ. Translation-rotation decoupling of tracers reflects medium-range crystalline order in two-dimensional colloid glasses. Phys Rev E 2021; 104:054615. [PMID: 34942845 DOI: 10.1103/physreve.104.054615] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2021] [Accepted: 11/17/2021] [Indexed: 11/07/2022]
Abstract
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.
Collapse
Affiliation(s)
- Dong Jae Chun
- Department of Chemistry, Sogang University, Seoul 04107, Republic of Korea
| | - Younghoon Oh
- Department of Chemistry, Sogang University, Seoul 04107, Republic of Korea
| | - Bong June Sung
- Department of Chemistry, Sogang University, Seoul 04107, Republic of Korea
| |
Collapse
|
12
|
McClure JE, Berg S, Armstrong RT. Thermodynamics of fluctuations based on time-and-space averages. Phys Rev E 2021; 104:035106. [PMID: 34654200 DOI: 10.1103/physreve.104.035106] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2020] [Accepted: 06/23/2021] [Indexed: 11/07/2022]
Abstract
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.
Collapse
Affiliation(s)
- James E McClure
- Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061, USA
| | - Steffen Berg
- Shell Global Solutions International B.V., Grasweg 31, 1031HW Amsterdam, The Netherlands
| | - Ryan T Armstrong
- School of Minerals and Energy Resources Engineering, University of New South Wales, Sydney 2052, Australia
| |
Collapse
|
13
|
Matyushov DV. Ewald sum corrections in simulations of ion and dipole solvation and electron transfer. J Chem Phys 2021; 155:114110. [PMID: 34551535 DOI: 10.1063/5.0061644] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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 m2/L3 with the solute dipole m and the box size L.
Collapse
Affiliation(s)
- Dmitry V Matyushov
- School of Molecular Sciences and Department of Physics, Arizona State University, P.O. Box 871504, Tempe, Arizona 85287-1504, USA
| |
Collapse
|
14
|
Jung G, Schmid F. Fluctuation-dissipation relations far from equilibrium: a case study. SOFT MATTER 2021; 17:6413-6425. [PMID: 34132298 PMCID: PMC8262459 DOI: 10.1039/d1sm00521a] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/07/2021] [Accepted: 06/06/2021] [Indexed: 06/12/2023]
Abstract
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.
Collapse
Affiliation(s)
- Gerhard Jung
- Institut für Theoretische Physik, Universität Innsbruck, Technikerstraße 21A, A-6020 Innsbruck, Austria.
| | - Friederike Schmid
- Institut für Physik, Johannes Gutenberg-Universität Mainz, 55099 Mainz, Germany.
| |
Collapse
|
15
|
Zhang J, Zheng W, Zhang S, Xu D, Nie Y, Jiang Z, Xu N. Unifying fluctuation-dissipation temperatures of slow-evolving nonequilibrium systems from the perspective of inherent structures. SCIENCE ADVANCES 2021; 7:eabg6766. [PMID: 34321210 PMCID: PMC8318365 DOI: 10.1126/sciadv.abg6766] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/22/2021] [Accepted: 06/11/2021] [Indexed: 06/13/2023]
Abstract
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.
Collapse
Affiliation(s)
- Jianhua Zhang
- Hefei National Laboratory for Physical Sciences at the Microscale, CAS Key Laboratory of Microscale Magnetic Resonance, Department of Physics, University of Science and Technology of China, Hefei 230026, P. R. China
| | - Wen Zheng
- Hefei National Laboratory for Physical Sciences at the Microscale, CAS Key Laboratory of Microscale Magnetic Resonance, Department of Physics, University of Science and Technology of China, Hefei 230026, P. R. China
| | - Shiyun Zhang
- Hefei National Laboratory for Physical Sciences at the Microscale, CAS Key Laboratory of Microscale Magnetic Resonance, Department of Physics, University of Science and Technology of China, Hefei 230026, P. R. China
| | - Ding Xu
- Hefei National Laboratory for Physical Sciences at the Microscale, CAS Key Laboratory of Microscale Magnetic Resonance, Department of Physics, University of Science and Technology of China, Hefei 230026, P. R. China
| | - Yunhuan Nie
- Hefei National Laboratory for Physical Sciences at the Microscale, CAS Key Laboratory of Microscale Magnetic Resonance, Department of Physics, University of Science and Technology of China, Hefei 230026, P. R. China
| | - Zhehua Jiang
- Hefei National Laboratory for Physical Sciences at the Microscale, CAS Key Laboratory of Microscale Magnetic Resonance, Department of Physics, University of Science and Technology of China, Hefei 230026, P. R. China
| | - Ning Xu
- Hefei National Laboratory for Physical Sciences at the Microscale, CAS Key Laboratory of Microscale Magnetic Resonance, Department of Physics, University of Science and Technology of China, Hefei 230026, P. R. China.
| |
Collapse
|
16
|
Shakerpoor A, Flenner E, Szamel G. The Einstein effective temperature can predict the tagged active particle density. J Chem Phys 2021; 154:184901. [PMID: 34241010 DOI: 10.1063/5.0049239] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Alireza Shakerpoor
- Department of Chemistry, Colorado State University, Fort Collins, Colorado 80523, USA
| | - Elijah Flenner
- Department of Chemistry, Colorado State University, Fort Collins, Colorado 80523, USA
| | - Grzegorz Szamel
- Department of Chemistry, Colorado State University, Fort Collins, Colorado 80523, USA
| |
Collapse
|
17
|
Hayashi K, Miyamoto MG, Niwa S. Effects of dynein inhibitor on the number of motor proteins transporting synaptic cargos. Biophys J 2021; 120:1605-1614. [PMID: 33617835 PMCID: PMC8204214 DOI: 10.1016/j.bpj.2021.02.018] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2020] [Revised: 02/02/2021] [Accepted: 02/08/2021] [Indexed: 11/27/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Kumiko Hayashi
- Department of Applied Physics, Graduate School of Engineering, Tohoku University, Sendai, Japan.
| | - Miki G Miyamoto
- Department of Applied Physics, Graduate School of Engineering, Tohoku University, Sendai, Japan
| | - Shinsuke Niwa
- Frontier Research Institute for Interdisciplinary Sciences and Graduate School of Life Science, Tohoku University, Sendai, Japan
| |
Collapse
|
18
|
Gross B, Sanhedrai H, Shekhtman L, Havlin S. Interconnections between networks acting like an external field in a first-order percolation transition. Phys Rev E 2020; 101:022316. [PMID: 32168699 DOI: 10.1103/physreve.101.022316] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2019] [Accepted: 02/05/2020] [Indexed: 11/07/2022]
Abstract
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.
Collapse
Affiliation(s)
- Bnaya Gross
- Department of Physics, Bar Ilan University, Ramat Gan, Israel
| | | | - Louis Shekhtman
- Network Science Institute, Northeastern University, Boston 02115, USA
| | - Shlomo Havlin
- Department of Physics, Bar Ilan University, Ramat Gan, Israel.,Institute of Innovative Research, Tokyo Institute of Technology, Midori-ku, Yokohama 226-8503, Japan
| |
Collapse
|
19
|
Abstract
Electron transfer in nonpolar media violates the temperature scaling predicted by the fluctuation–dissipation theorem.
Collapse
Affiliation(s)
- Dmitry V. Matyushov
- Department of Physics and School of Molecular Sciences
- Arizona State University
- Tempe
- USA
| |
Collapse
|
20
|
Sarracino A, Vulpiani A. On the fluctuation-dissipation relation in non-equilibrium and non-Hamiltonian systems. CHAOS (WOODBURY, N.Y.) 2019; 29:083132. [PMID: 31472486 DOI: 10.1063/1.5110262] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/15/2019] [Accepted: 08/08/2019] [Indexed: 06/10/2023]
Abstract
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.
Collapse
Affiliation(s)
- A Sarracino
- Dipartimento di Ingegneria, Università della Campania "L. Vanvitelli," via Roma 29, 81031 Aversa (CE), Italy
| | - A Vulpiani
- Dipartimento di Fisica, Università Sapienza-p.le A. Moro 2, 00185 Roma, Italy
| |
Collapse
|
21
|
Berthier L, Ozawa M, Scalliet C. Configurational entropy of glass-forming liquids. J Chem Phys 2019; 150:160902. [PMID: 31042883 DOI: 10.1063/1.5091961] [Citation(s) in RCA: 46] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Ludovic Berthier
- Laboratoire Charles Coulomb (L2C), Université de Montpellier, CNRS, Montpellier, France
| | - Misaki Ozawa
- Laboratoire Charles Coulomb (L2C), Université de Montpellier, CNRS, Montpellier, France
| | - Camille Scalliet
- Laboratoire Charles Coulomb (L2C), Université de Montpellier, CNRS, Montpellier, France
| |
Collapse
|
22
|
Hasegawa S, Sagawa T, Ikeda K, Okada Y, Hayashi K. Investigation of multiple-dynein transport of melanosomes by non-invasive force measurement using fluctuation unit χ. Sci Rep 2019; 9:5099. [PMID: 30911050 PMCID: PMC6433852 DOI: 10.1038/s41598-019-41458-w] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2018] [Accepted: 03/05/2019] [Indexed: 12/19/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Shin Hasegawa
- Department of Applied Physics, Graduate School of Engineering, Tohoku University, Sendai, Japan
| | - Takashi Sagawa
- Advanced ICT Research Institute, National Institute of Information and Communications Technology, Kobe, Japan
| | - Kazuho Ikeda
- Laboratory for Cell Dynamics Observation, Center for Biosystems Dynamics Research, RIKEN, Osaka, Japan
| | - Yasushi Okada
- Laboratory for Cell Dynamics Observation, Center for Biosystems Dynamics Research, RIKEN, Osaka, Japan.,Department of Physics and Universal Biology Institute, Graduate School of Science, The University of Tokyo, Tokyo, Japan.,Department of Physics, Universal Biology Institute, and the International Research Center for Neurointelligence (WPI-IRCN), The University of Tokyo, Tokyo, Japan
| | - Kumiko Hayashi
- Department of Applied Physics, Graduate School of Engineering, Tohoku University, Sendai, Japan.
| |
Collapse
|
23
|
Hayashi K. Application of the fluctuation theorem to motor proteins: from F 1-ATPase to axonal cargo transport by kinesin and dynein. Biophys Rev 2018; 10:1311-1321. [PMID: 30019204 PMCID: PMC6233339 DOI: 10.1007/s12551-018-0440-5] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2018] [Accepted: 07/02/2018] [Indexed: 12/27/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Kumiko Hayashi
- Department of Applied Physics, Graduate School of Engineering, Tohoku University, Sendai, Japan.
| |
Collapse
|
24
|
Ooshida T, Otsuki M. Two-tag correlations and nonequilibrium fluctuation-response relation in ageing single-file diffusion. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2018; 30:374001. [PMID: 30027890 DOI: 10.1088/1361-648x/aad4cc] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
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.
Collapse
Affiliation(s)
- Takeshi Ooshida
- Department of Mechanical and Physical Engineering, Tottori University, Tottori 680-8552, Japan
| | | |
Collapse
|
25
|
Baity-Jesi M, Achard-de Lustrac A, Biroli G. Activated dynamics: An intermediate model between the random energy model and the p-spin model. Phys Rev E 2018; 98:012133. [PMID: 30110833 DOI: 10.1103/physreve.98.012133] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2018] [Indexed: 11/07/2022]
Abstract
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.
Collapse
Affiliation(s)
- Marco Baity-Jesi
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | | | - Giulio Biroli
- Institut de Physique Théorique, Université Paris Saclay, CEA, CNRS, F-91191 Gif-sur-Yvette, France and Laboratoire de Physique Statistique, École Normale Supérieure, PSL Research University, 24 rue Lhomond, 75005 Paris, France
| |
Collapse
|
26
|
Abstract
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 Ts/T, in which T is the bath temperature and Ts 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, Ts ≅ 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.
Collapse
Affiliation(s)
- Jeppe C Dyre
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
27
|
Crisanti A, Picco M, Ritort F. Derivation of the spin-glass order parameter from stochastic thermodynamics. Phys Rev E 2018; 97:052103. [PMID: 29906861 DOI: 10.1103/physreve.97.052103] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2017] [Indexed: 06/08/2023]
Abstract
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.
Collapse
Affiliation(s)
- A Crisanti
- Dipartimento di Fisica, Università di Roma "La Sapienza", and Istituto dei Sistemi Complessi ISC-CNR, P.le Aldo Moro 2, 00185 Roma, Italy
| | - M Picco
- Sorbonne Université, CNRS, Université Pierre et Marie Curie, LPTHE, F-75005 Paris, France
| | - F Ritort
- Small Biosystems Lab, Condensed Matter Physics Department, Universitat de Barcelona, C/ Martí i Franquès, E-08028, Barcelona, Spain
- CIBER-BBN Center for Bioengineering, Biomaterials and Nanomedicine, Instituto de Salud Carlos III, Madrid, Spain
| |
Collapse
|
28
|
Seyedi S, Matyushov DV. Ergodicity breaking of iron displacement in heme proteins. SOFT MATTER 2017; 13:8188-8201. [PMID: 29082406 DOI: 10.1039/c7sm01561e] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
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.
Collapse
Affiliation(s)
- Salman Seyedi
- Department of Physics, Arizona State University, PO Box 871504, Tempe, Arizona 85287, USA
| | | |
Collapse
|
29
|
Leibovich N, Barkai E. Conditional 1/f^{α} noise: From single molecules to macroscopic measurement. Phys Rev E 2017; 96:032132. [PMID: 29346920 DOI: 10.1103/physreve.96.032132] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2017] [Indexed: 06/07/2023]
Abstract
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.
Collapse
Affiliation(s)
- N Leibovich
- Department of Physics, Institute of Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan 5290002, Israel
| | - E Barkai
- Department of Physics, Institute of Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan 5290002, Israel
| |
Collapse
|
30
|
Diezemann G. Nonlinear response theory for Markov processes. II. Fifth-order response functions. Phys Rev E 2017; 96:022150. [PMID: 28950644 DOI: 10.1103/physreve.96.022150] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2017] [Indexed: 06/07/2023]
Abstract
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.
Collapse
Affiliation(s)
- Gregor Diezemann
- Institut für Physikalische Chemie, Universität Mainz, Duesbergweg 10-14, 55128 Mainz, Germany
| |
Collapse
|
31
|
Li Q, Peng X, McKenna GB. Long-term aging behaviors in a model soft colloidal system. SOFT MATTER 2017; 13:1396-1404. [PMID: 28120996 DOI: 10.1039/c6sm02408d] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
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.
Collapse
Affiliation(s)
- Qi Li
- Department of Chemical Engineering, Texas Tech University, Lubbock, Texas 79409, USA.
| | - Xiaoguang Peng
- Department of Chemical Engineering, Texas Tech University, Lubbock, Texas 79409, USA.
| | - Gregory B McKenna
- Department of Chemical Engineering, Texas Tech University, Lubbock, Texas 79409, USA.
| |
Collapse
|
32
|
Leibovich N, Dechant A, Lutz E, Barkai E. Aging Wiener-Khinchin theorem and critical exponents of 1/f^{β} noise. Phys Rev E 2016; 94:052130. [PMID: 27967149 DOI: 10.1103/physreve.94.052130] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2016] [Indexed: 06/06/2023]
Abstract
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.
Collapse
Affiliation(s)
- N Leibovich
- Department of Physics, Institute of Nanotechnology and Advanced Materials, Bar Ilan University, Ramat-Gan 52900, Israel
| | - A Dechant
- Department of Physics, Institute of Nanotechnology and Advanced Materials, Bar Ilan University, Ramat-Gan 52900, Israel
- Department of Physics, Friedrich-Alexander-Universität Erlangen-Nürnberg, 91058 Erlangen, Germany
| | - E Lutz
- Department of Physics, Friedrich-Alexander-Universität Erlangen-Nürnberg, 91058 Erlangen, Germany
| | - E Barkai
- Department of Physics, Institute of Nanotechnology and Advanced Materials, Bar Ilan University, Ramat-Gan 52900, Israel
| |
Collapse
|
33
|
Sarracino A, Cecconi F, Puglisi A, Vulpiani A. Nonlinear Response of Inertial Tracers in Steady Laminar Flows: Differential and Absolute Negative Mobility. PHYSICAL REVIEW LETTERS 2016; 117:174501. [PMID: 27824440 DOI: 10.1103/physrevlett.117.174501] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/01/2016] [Indexed: 06/06/2023]
Abstract
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.
Collapse
Affiliation(s)
- A Sarracino
- CNR-ISC and Dipartimento di Fisica, Sapienza Università di Roma, piazzale A. Moro 2, 00185 Roma, Italy
| | - F Cecconi
- CNR-ISC and Dipartimento di Fisica, Sapienza Università di Roma, piazzale A. Moro 2, 00185 Roma, Italy
| | - A Puglisi
- CNR-ISC and Dipartimento di Fisica, Sapienza Università di Roma, piazzale A. Moro 2, 00185 Roma, Italy
| | - A Vulpiani
- Dipartimento di Fisica, Sapienza Università di Roma, piazzale A. Moro 2, 00185 Roma, Italy
| |
Collapse
|
34
|
Dinpajooh M, Martin DR, Matyushov DV. Polarizability of the active site of cytochrome c reduces the activation barrier for electron transfer. Sci Rep 2016; 6:28152. [PMID: 27306204 PMCID: PMC4910110 DOI: 10.1038/srep28152] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2016] [Accepted: 05/27/2016] [Indexed: 11/12/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Mohammadhasan Dinpajooh
- Department of Physics and School of Molecular Sciences, Arizona State University, PO Box 871504, Tempe, AZ 85287-1504, USA
| | - Daniel R Martin
- Department of Physics and School of Molecular Sciences, Arizona State University, PO Box 871504, Tempe, AZ 85287-1504, USA
| | - Dmitry V Matyushov
- Department of Physics and School of Molecular Sciences, Arizona State University, PO Box 871504, Tempe, AZ 85287-1504, USA
| |
Collapse
|
35
|
Martin DR, Matyushov DV. Communication: Microsecond dynamics of the protein and water affect electron transfer in a bacterial bc(1) complex. J Chem Phys 2016; 142:161101. [PMID: 25933744 DOI: 10.1063/1.4919222] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Daniel R Martin
- Department of Physics and Department of Chemistry and Biochemistry, Arizona State University, P.O. Box 871504, Tempe, Arizona 85287, USA
| | - Dmitry V Matyushov
- Department of Physics and Department of Chemistry and Biochemistry, Arizona State University, P.O. Box 871504, Tempe, Arizona 85287, USA
| |
Collapse
|
36
|
Matyushov DV. Protein electron transfer: is biology (thermo)dynamic? JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2015; 27:473001. [PMID: 26558324 DOI: 10.1088/0953-8984/27/47/473001] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
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.
Collapse
Affiliation(s)
- Dmitry V Matyushov
- Department of Physics and School of Molecular Sciences, Arizona State University, PO Box 871504, Tempe, AZ 85287-1504, USA
| |
Collapse
|
37
|
Sakuldee F, Suwanna S. Linear response and modified fluctuation-dissipation relation in random potential. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:052118. [PMID: 26651658 DOI: 10.1103/physreve.92.052118] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/06/2015] [Indexed: 06/05/2023]
Abstract
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.
Collapse
Affiliation(s)
- Fattah Sakuldee
- MU-NECTEC Collaborative Research Unit on Quantum Information, Department of Physics, Faculty of Science, Mahidol University, Bangkok 10400, Thailand
| | - Sujin Suwanna
- MU-NECTEC Collaborative Research Unit on Quantum Information, Department of Physics, Faculty of Science, Mahidol University, Bangkok 10400, Thailand
| |
Collapse
|
38
|
Leibovich N, Barkai E. Aging Wiener-Khinchin Theorem. PHYSICAL REVIEW LETTERS 2015; 115:080602. [PMID: 26340172 DOI: 10.1103/physrevlett.115.080602] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/06/2015] [Indexed: 06/05/2023]
Abstract
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.
Collapse
Affiliation(s)
- N Leibovich
- Department of Physics, Institute of Nanotechnology and Advanced Materials, Bar Ilan University, Ramat-Gan 52900, Israel
| | - E Barkai
- Department of Physics, Institute of Nanotechnology and Advanced Materials, Bar Ilan University, Ramat-Gan 52900, Israel
| |
Collapse
|
39
|
Medvedev GA, Caruthers JM. On Thermodynamic Consistency of a Stochastic Constitutive Model for Glassy Polymers. Ind Eng Chem Res 2015. [DOI: 10.1021/acs.iecr.5b01347] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Grigori A. Medvedev
- School of Chemical Engineering Purdue University West Lafayette, Indiana 47907, United States
| | - James M. Caruthers
- School of Chemical Engineering Purdue University West Lafayette, Indiana 47907, United States
| |
Collapse
|
40
|
Martin DR, Matyushov DV. Photosynthetic diode: electron transport rectification by wetting the quinone cofactor. Phys Chem Chem Phys 2015; 17:22523-8. [PMID: 26171665 DOI: 10.1039/c5cp03397g] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
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.
Collapse
Affiliation(s)
- Daniel R Martin
- Department of Physics, Arizona State University, PO Box 871504, Tempe, Arizona 85287, USA
| | | |
Collapse
|
41
|
Ninarello AS, Gnan N, Sciortino F. Observable-dependence of the effective temperature in off-equilibrium diatomic molecular liquids. J Chem Phys 2014; 141:194507. [PMID: 25416898 DOI: 10.1063/1.4901526] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
| | - Nicoletta Gnan
- CNR-ISC Uos "Sapienza," Piazzale A. Moro 2, 00186 Roma, Italy
| | - Francesco Sciortino
- Dipartimento di Fisica, "Sapienza" Università di Roma, Piazzale A. Moro 2, 00186 Roma, Italy
| |
Collapse
|
42
|
Papo D. Functional significance of complex fluctuations in brain activity: from resting state to cognitive neuroscience. Front Syst Neurosci 2014; 8:112. [PMID: 24966818 PMCID: PMC4052734 DOI: 10.3389/fnsys.2014.00112] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2014] [Accepted: 05/26/2014] [Indexed: 11/15/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- David Papo
- Computational Systems Biology Group, Center for Biomedical Technology, Universidad Politécnica de Madrid Madrid, Spain
| |
Collapse
|
43
|
Papo D. Measuring brain temperature without a thermometer. Front Physiol 2014; 5:124. [PMID: 24723893 PMCID: PMC3973909 DOI: 10.3389/fphys.2014.00124] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2014] [Accepted: 03/13/2014] [Indexed: 11/13/2022] Open
Affiliation(s)
- David Papo
- Computational Systems Biology Group, Center for Biomedical Technology, Universidad Politécnica de Madrid Madrid, Spain
| |
Collapse
|
44
|
Lippiello E, Baiesi M, Sarracino A. Nonequilibrium fluctuation-dissipation theorem and heat production. PHYSICAL REVIEW LETTERS 2014; 112:140602. [PMID: 24765939 DOI: 10.1103/physrevlett.112.140602] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/14/2013] [Indexed: 06/03/2023]
Abstract
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.
Collapse
Affiliation(s)
- E Lippiello
- Department of Mathematics and Physics, Second University of Naples, Via Vivaldi 43, 81100 Caserta, Italy
| | - M Baiesi
- Department of Physics and Astronomy, University of Padua, Via Marzolo 8, I-35131 Padova, Italy and INFN, Sezione di Padova, Via Marzolo 8, I-35131 Padova, Italy
| | - A Sarracino
- ISC-CNR and Dipartimento di Fisica, Università "Sapienza", piazzale A. Moro 2, 00185 Rome, Italy and Laboratoire de Physique Théorique de la Matière Condensée, CNRS UMR 7600, case courrier 121, Université Paris 6, 4 Place Jussieu, 75255 Paris Cedex, France
| |
Collapse
|
45
|
Gnoli A, Puglisi A, Sarracino A, Vulpiani A. Nonequilibrium Brownian motion beyond the effective temperature. PLoS One 2014; 9:e93720. [PMID: 24714671 PMCID: PMC3979703 DOI: 10.1371/journal.pone.0093720] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2014] [Accepted: 03/05/2014] [Indexed: 11/19/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Andrea Gnoli
- Istituto dei Sistemi Complessi - Consiglio Nazionale delle Ricerche, Rome, Italy
- Dipartimento di Fisica, Università “Sapienza”, Rome, Italy
| | - Andrea Puglisi
- Istituto dei Sistemi Complessi - Consiglio Nazionale delle Ricerche, Rome, Italy
- Dipartimento di Fisica, Università “Sapienza”, Rome, Italy
- * E-mail:
| | - Alessandro Sarracino
- Istituto dei Sistemi Complessi - Consiglio Nazionale delle Ricerche, Rome, Italy
- Dipartimento di Fisica, Università “Sapienza”, Rome, Italy
- Laboratoire de Physique Théorique de la Matière Condensée - Centre National de la Recherche Scientifique Unité mixte de recherche 7600, Université Paris 6, Paris, France
| | - Angelo Vulpiani
- Istituto dei Sistemi Complessi - Consiglio Nazionale delle Ricerche, Rome, Italy
- Dipartimento di Fisica, Università “Sapienza”, Rome, Italy
| |
Collapse
|
46
|
Hagmann JG, Nakagawa N, Peyrard M. Characterization of the low-temperature properties of a simplified protein model. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 89:012705. [PMID: 24580255 DOI: 10.1103/physreve.89.012705] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/22/2013] [Indexed: 06/03/2023]
Abstract
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.
Collapse
Affiliation(s)
- Johannes-Geert Hagmann
- Laboratoire de Physique, Ecole Normale Supérieure de Lyon, CNRS, 46 Allée d'Italie, 69364 Lyon, France
| | - Naoko Nakagawa
- Department of Mathematical Sciences, Ibaraki University, Mito, Ibaraki 310-8512, Japan
| | - Michel Peyrard
- Laboratoire de Physique, Ecole Normale Supérieure de Lyon, CNRS, 46 Allée d'Italie, 69364 Lyon, France
| |
Collapse
|
47
|
Hurowitz D, Rahav S, Cohen D. Nonequilibrium steady state and induced currents of a mesoscopically glassy system: interplay of resistor-network theory and Sinai physics. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 88:062141. [PMID: 24483419 DOI: 10.1103/physreve.88.062141] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/25/2013] [Indexed: 06/03/2023]
Abstract
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.
Collapse
Affiliation(s)
- Daniel Hurowitz
- Department of Physics, Ben-Gurion University of the Negev, Beer-Sheva, Israel
| | - Saar Rahav
- Schulich Faculty of Chemistry, Technion-Israel Institute of Technology, Haifa 32000, Israel
| | - Doron Cohen
- Department of Physics, Ben-Gurion University of the Negev, Beer-Sheva, Israel
| |
Collapse
|
48
|
Sarracino A. Time asymmetry of the Kramers equation with nonlinear friction: fluctuation-dissipation relation and ratchet effect. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 88:052124. [PMID: 24329231 DOI: 10.1103/physreve.88.052124] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/13/2013] [Indexed: 06/03/2023]
Abstract
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.
Collapse
Affiliation(s)
- A Sarracino
- CNR-ISC and Dipartimento di Fisica, Università Sapienza, p.le A. Moro 2, 00185 Roma, Italy
| |
Collapse
|
49
|
Klotsa D, Jack RL. Controlling crystal self-assembly using a real-time feedback scheme. J Chem Phys 2013; 138:094502. [DOI: 10.1063/1.4793527] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
50
|
Crisanti A, Picco M, Ritort F. Fluctuation relation for weakly ergodic aging systems. PHYSICAL REVIEW LETTERS 2013; 110:080601. [PMID: 23473124 DOI: 10.1103/physrevlett.110.080601] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/13/2012] [Indexed: 06/01/2023]
Abstract
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.
Collapse
Affiliation(s)
- A Crisanti
- Dipartimento di Fisica, Università di Roma La Sapienza, Piazzale Aldo Moro 2, 00185 Roma, Italy
| | | | | |
Collapse
|