Emergence of Memory in Equilibrium versus Nonequilibrium Systems.

Experiments often probe observables that correspond to low-dimensional projections of high-dimensional dynamics. In such situations distinct microscopic configurations become lumped into the same observable state. It is well known that correlations between the observable and the hidden degrees of freedom give rise to memory effects. However, how and under which conditions these correlations emerge remain poorly understood. Here we shed light on two fundamentally different scenarios of the emergence of memory in minimal stationary systems, where observed and hidden degrees of freedom either evolve cooperatively or are coupled by a hidden nonequilibrium current. In the reversible setting the strongest memory manifests when the timescales of hidden and observed dynamics overlap, whereas, strikingly, in the driven setting maximal memory emerges under a clear timescale separation. Our results hint at the possibility of fundamental differences in the way memory emerges in equilibrium versus driven systems that may be utilized as a "diagnostic" of the underlying hidden transport mechanism.

Milestoning estimators of dissipation in systems observed at a coarse resolution.

Many nonequilibrium, active processes are observed at a coarse-grained level, where different microscopic configurations are projected onto the same observable state. Such "lumped" observables display memory, and in many cases, the irreversible character of the underlying microscopic dynamics becomes blurred, e.g., when the projection hides dissipative cycles. As a result, the observations appear less irreversible, and it is very challenging to infer the degree of broken time-reversal symmetry. Here we show, contrary to intuition, that by ignoring parts of the already coarse-grained state space we may-via a process called milestoning-improve entropy-production estimates. We present diverse examples where milestoning systematically renders observations "closer to underlying microscopic dynamics" and thereby improves thermodynamic inference from lumped data assuming a given range of memory, and we hypothesize that this effect is quite general. Moreover, whereas the correct general physical definition of time reversal in the presence of memory remains unknown, we here show by means of physically relevant examples that at least for semi-Markov processes of first and second order, waiting-time contributions arising from adopting a naive Markovian definition of time reversal generally must be discarded.

Controlling Uncertainty of Empirical First-Passage Times in the Small-Sample Regime.

We derive general bounds on the probability that the empirical first-passage time τ[over ¯]_{n}≡∑_{i=1}^{n}τ_{i}/n of a reversible ergodic Markov process inferred from a sample of n independent realizations deviates from the true mean first-passage time by more than any given amount in either direction. We construct nonasymptotic confidence intervals that hold in the elusive small-sample regime and thus fill the gap between asymptotic methods and the Bayesian approach that is known to be sensitive to prior belief and tends to underestimate uncertainty in the small-sample setting. We prove sharp bounds on extreme first-passage times that control uncertainty even in cases where the mean alone does not sufficiently characterize the statistics. Our concentration-of-measure-based results allow for model-free error control and reliable error estimation in kinetic inference, and are thus important for the analysis of experimental and simulation data in the presence of limited sampling.

Direct Route to Thermodynamic Uncertainty Relations and Their Saturation.

Thermodynamic uncertainty relations (TURs) bound the dissipation in nonequilibrium systems from below by fluctuations of an observed current. Contrasting the elaborate techniques employed in existing proofs, we here prove TURs directly from the Langevin equation. This establishes the TUR as an inherent property of overdamped stochastic equations of motion. In addition, we extend the transient TUR to currents and densities with explicit time dependence. By including current-density correlations we, moreover, derive a new sharpened TUR for transient dynamics. Our arguably simplest and most direct proof, together with the new generalizations, allows us to systematically determine conditions under which the different TURs saturate and thus allows for a more accurate thermodynamic inference. Finally, we outline the direct proof also for Markov jump dynamics.

Challenges in Inferring the Directionality of Active Molecular Processes from Single-Molecule Fluorescence Resonance Energy Transfer Trajectories.

We discuss some of the practical challenges that one faces in using stochastic thermodynamics to infer directionality of molecular machines from experimental single-molecule trajectories. Because of the limited spatiotemporal resolution of single-molecule experiments and because both forward and backward transitions between the same pairs of states cannot always be detected, differentiating between the forward and backward directions of, e.g., an ATP-consuming molecular machine that operates periodically, turns out to be a nontrivial task. Using a simple extension of a Markov-state model that is commonly employed to analyze single-molecule transition-path measurements, we illustrate how irreversibility can be hidden from such measurements but in some cases can be uncovered when non-Markov effects in low-dimensional single-molecule trajectories are considered.

Mathematical, Thermodynamical, and Experimental Necessity for Coarse Graining Empirical Densities and Currents in Continuous Space.

We present general results on fluctuations and spatial correlations of the coarse-grained empirical density and current of Markovian diffusion in equilibrium or nonequilibrium steady states on all timescales. We unravel a deep connection between current fluctuations and generalized time-reversal symmetry, providing new insight into time-averaged observables. We highlight the essential role of coarse graining in space from mathematical, thermodynamical, and experimental points of view. Spatial coarse graining is required to uncover salient features of currents that break detailed balance, and a thermodynamically "optimal" coarse graining ensures the most precise inference of dissipation. Defined without coarse graining, the fluctuations of empirical density and current are proven to diverge on all timescales in dimensions higher than one, which has far-reaching consequences for the central-limit regime in continuous space. We apply the results to examples of irreversible diffusion. Our findings provide new intuition about time-averaged observables and allow for a more efficient analysis of single-molecule experiments.

Non-ergodicity of a globular protein extending beyond its functional timescale.

Internal motions of folded proteins have been assumed to be ergodic, i.e., that the dynamics of a single protein molecule averaged over a very long time resembles that of an ensemble. Here, by performing single-molecule fluorescence resonance energy transfer (smFRET) experiments and molecular dynamics (MD) simulations of a multi-domain globular protein, cytoplasmic protein-tyrosine phosphatase (SHP2), we demonstrate that the functional inter-domain motion is observationally non-ergodic over the time spans 10-12 to 10-7 s and 10-1 to 102 s. The difference between observational non-ergodicity and simple non-convergence is discussed. In comparison, a single-strand DNA of similar size behaves ergodically with an energy landscape resembling a one-dimensional linear chain. The observed non-ergodicity results from the hierarchical connectivity of the high-dimensional energy landscape of the protein molecule. As the characteristic time for the protein to conduct its dephosphorylation function is ∼10 s, our findings suggest that, due to the non-ergodicity, individual, seemingly identical protein molecules can be dynamically and functionally different.

Erratum: Faster Uphill Relaxation in Thermodynamically Equidistant Temperature Quenches [Phys. Rev. Lett. 125, 110602 (2020)].

This corrects the article DOI: 10.1103/PhysRevLett.125.110602.

First-passage statistics of colloids on fractals: Theory and experimental realization.

In nature and technology, particle dynamics frequently occur in complex environments, for example in restricted geometries or crowded media. These dynamics have often been modeled invoking a fractal structure of the medium although the fractal structure was only indirectly inferred through the dynamics. Moreover, systematic studies have not yet been performed. Here, colloidal particles moving in a laser speckle pattern are used as a model system. In this case, the experimental observations can be reliably traced to the fractal structure of the underlying medium with an adjustable fractal dimension. First-passage time statistics reveal that the particles explore the speckle in a self-similar, fractal manner at least over four decades in time and on length scales up to 20 times the particle radius. The requirements for fractal diffusion to be applicable are laid out, and methods to extract the fractal dimension are established.

Thermodynamic Uncertainty Relation Bounds the Extent of Anomalous Diffusion.

In a finite system driven out of equilibrium by a constant external force the thermodynamic uncertainty relation (TUR) bounds the variance of the conjugate current variable by the thermodynamic cost of maintaining the nonequilibrium stationary state. Here we highlight a new facet of the TUR by showing that it also bounds the timescale on which a finite system can exhibit anomalous kinetics. In particular, we demonstrate that the TUR bounds subdiffusion in a single file confined to a ring as well as a dragged Gaussian polymer chain even when detailed balance is satisfied. Conversely, the TUR bounds the onset of superdiffusion in the active comb model. Remarkably, the fluctuations in a comb model evolving from a steady state behave anomalously as soon as detailed balance is broken. Our work establishes a link between stochastic thermodynamics and the field of anomalous dynamics that will fertilize further investigations of thermodynamic consistency of anomalous diffusion models.

Single-file diffusion in a bi-stable potential: Signatures of memory in the barrier-crossing of a tagged-particle.

We investigate memory effects in barrier-crossing in the overdamped setting. We focus on the scenario where the hidden degrees of freedom relax on exactly the same time scale as the observable. As a prototypical model, we analyze tagged-particle diffusion in a single file confined to a bi-stable potential. We identify the signatures of memory and explain their origin. The emerging memory is a result of the projection of collective many-body eigenmodes onto the motion of a tagged-particle. We are interested in the "confining" (all background particles in front of the tagged-particle) and "pushing" (all background particles behind the tagged-particle) scenarios for which we find non-trivial and qualitatively different relaxation behaviors. Notably and somewhat unexpectedly, at a fixed particle number, we find that the higher the barrier, the stronger the memory effects are. The fact that the external potential alters the memory is important more generally and should be taken into account in applications of generalized Langevin equations. Our results can readily be tested experimentally and may be relevant for understanding transport in biological ion-channels.

Faster Uphill Relaxation in Thermodynamically Equidistant Temperature Quenches.

We uncover an unforeseen asymmetry in relaxation: for a pair of thermodynamically equidistant temperature quenches, one from a lower and the other from a higher temperature, the relaxation at the ambient temperature is faster in the case of the former. We demonstrate this finding on hand of two exactly solvable many-body systems relevant in the context of single-molecule and tracer-particle dynamics. We prove that near stable minima and for all quadratic energy landscapes it is a general phenomenon that also exists in a class of non-Markovian observables probed in single-molecule and particle-tracking experiments. The asymmetry is a general feature of reversible overdamped diffusive systems with smooth single-well potentials and occurs in multiwell landscapes when quenches disturb predominantly intrawell equilibria. Our findings may be relevant for the optimization of stochastic heat engines.

Quantifying non-ergodicity of anomalous diffusion with higher order moments.

Anomalous diffusion is being discovered in a fast growing number of systems. The exact nature of this anomalous diffusion provides important information on the physical laws governing the studied system. One of the central properties analysed for finite particle motion time series is the intrinsic variability of the apparent diffusivity, typically quantified by the ergodicity breaking parameter EB. Here we demonstrate that frequently EB is insufficient to provide a meaningful measure for the observed variability of the data. Instead, important additional information is provided by the higher order moments entering by the skewness and kurtosis. We analyse these quantities for three popular anomalous diffusion models. In particular, we find that even for the Gaussian fractional Brownian motion a significant skewness in the results of physical measurements occurs and needs to be taken into account. Interestingly, the kurtosis and skewness may also provide sensitive estimates of the anomalous diffusion exponent underlying the data. We also derive a new result for the EB parameter of fractional Brownian motion valid for the whole range of the anomalous diffusion parameter. Our results are important for the analysis of anomalous diffusion but also provide new insights into the theory of anomalous stochastic processes.

Materials Nanoarchitecturing via Cation-Mediated Protein Assembly: Making Limpet Teeth without Mineral.

Teeth are designed to deliver high forces while withstanding the generated stresses. Aside from isolated mineral-free exception (e.g., marine polychaetes and squids), minerals are thought to be indispensable for tooth-hardening and durability. Here, the unmineralized teeth of the giant keyhole limpet (Megathura crenulata) are shown to attain a stiffness, which is twofold higher than any known organic biogenic structures. In these teeth, protein and chitin fibers establish a stiff compact outer shell enclosing a less compact core. The stiffness and its gradients emerge from a concerted interaction across multiple length-scales: packing of hydrophobic proteins and folding into secondary structures mediated by Ca2+ and Mg2+ together with a strong spatial control in the local fiber orientation. These results integrating nanoindentation, acoustic microscopy, and finite-element modeling for probing the tooth's mechanical properties, spatially resolved small- and wide-angle X-ray scattering for probing the material ordering on the micrometer scale, and energy-dispersive X-ray scattering combined with confocal Raman microscopy to study structural features on the molecular scale, reveal a nanocomposite structure hierarchically assembled to form a versatile damage-tolerant protein-based tooth, with a stiffness similar to mineralized mammalian bone, but without any mineral.

First passage time distribution in heterogeneity controlled kinetics: going beyond the mean first passage time.

The first passage is a generic concept for quantifying when a random quantity such as the position of a diffusing molecule or the value of a stock crosses a preset threshold (target) for the first time. The last decade saw an enlightening series of new results focusing mostly on the so-called mean and global first passage time (MFPT and GFPT, respectively) of such processes. Here we push the understanding of first passage processes one step further. For a simple heterogeneous system we derive rigorously the complete distribution of first passage times (FPTs). Our results demonstrate that the typical FPT significantly differs from the MFPT, which corresponds to the long time behaviour of the FPT distribution. Conversely, the short time behaviour is shown to correspond to trajectories connecting directly from the initial value to the target. Remarkably, we reveal a previously overlooked third characteristic time scale of the first passage dynamics mirroring brief excursion away from the target.

Signal focusing through active transport.

The accuracy of molecular signaling in biological cells and novel diagnostic devices is ultimately limited by the counting noise floor imposed by the thermal diffusion. Motivated by the fact that messenger RNA and vesicle-engulfed signaling molecules transiently bind to molecular motors and are actively transported in biological cells, we show here that the random active delivery of signaling particles to within a typical diffusion distance to the receptor generically reduces the correlation time of the counting noise. Considering a variety of signaling particle sizes from mRNA to vesicles and cell sizes from prokaryotic to eukaryotic cells, we show that the conditions for active focusing-faster and more precise signaling-are indeed compatible with observations in living cells. Our results improve the understanding of molecular cellular signaling and novel diagnostic devices.

Optimization and universality of Brownian search in a basic model of quenched heterogeneous media.

The kinetics of a variety of transport-controlled processes can be reduced to the problem of determining the mean time needed to arrive at a given location for the first time, the so-called mean first-passage time (MFPT) problem. The occurrence of occasional large jumps or intermittent patterns combining various types of motion are known to outperform the standard random walk with respect to the MFPT, by reducing oversampling of space. Here we show that a regular but spatially heterogeneous random walk can significantly and universally enhance the search in any spatial dimension. In a generic minimal model we consider a spherically symmetric system comprising two concentric regions with piecewise constant diffusivity. The MFPT is analyzed under the constraint of conserved average dynamics, that is, the spatially averaged diffusivity is kept constant. Our analytical calculations and extensive numerical simulations demonstrate the existence of an optimal heterogeneity minimizing the MFPT to the target. We prove that the MFPT for a random walk is completely dominated by what we term direct trajectories towards the target and reveal a remarkable universality of the spatially heterogeneous search with respect to target size and system dimensionality. In contrast to intermittent strategies, which are most profitable in low spatial dimensions, the spatially inhomogeneous search performs best in higher dimensions. Discussing our results alongside recent experiments on single-particle tracking in living cells, we argue that the observed spatial heterogeneity may be beneficial for cellular signaling processes.

First-passage statistics for aging diffusion in systems with annealed and quenched disorder.

Aging, the dependence of the dynamics of a physical process on the time ta since its original preparation, is observed in systems ranging from the motion of charge carriers in amorphous semiconductors over the blinking dynamics of quantum dots to the tracer dispersion in living biological cells. Here we study the effects of aging on one of the most fundamental properties of a stochastic process, the first-passage dynamics. We find that for an aging continuous time random walk process, the scaling exponent of the density of first-passage times changes twice as the aging progresses and reveals an intermediate scaling regime. The first-passage dynamics depends on ta differently for intermediate and strong aging. Similar crossovers are obtained for the first-passage dynamics for a confined and driven particle. Comparison to the motion of an aged particle in the quenched trap model with a bias shows excellent agreement with our analytical findings. Our results demonstrate how first-passage measurements can be used to unravel the age ta of a physical system.

Diffusion of finite-size particles in two-dimensional channels with random wall configurations.

Diffusion of chemicals or tracer molecules through complex systems containing irregularly shaped channels is important in many applications. Most theoretical studies based on the famed Fick-Jacobs equation focus on the idealised case of infinitely small particles and reflecting boundaries. In this study we use numerical simulations to consider the transport of finite-size particles through asymmetrical two-dimensional channels. Additionally, we examine transient binding of the molecules to the channel walls by applying sticky boundary conditions. We consider an ensemble of particles diffusing in independent channels, which are characterised by common structural parameters. We compare our results for the long-time effective diffusion coefficient with a recent theoretical formula obtained by Dagdug and Pineda [J. Chem. Phys., 2012, 137, 024107].

Soft collective fluctuations governing hydrophobic association.

The interaction between two associating hydrophobic particles has traditionally been explained in terms of the release of entropically frustrated hydration shell water molecules. However, this picture cannot account for the kinetics of hydrophobic association and is therefore not capable of providing a microscopic description of the hydrophobic interaction (HI). Here, Monte Carlo simulations of a pair of molecular-scale apolar solutes in aqueous solution reveal the critical role of collective fluctuations in the hydrogen bond (HB) network for the microscopic picture of the HI. The main contribution to the HI is the relaxation of solute-water translational correlations. The existence of a heat capacity maximum at the desolvation barrier is shown to arise from softening of non-HB water fluctuations and the relaxation of many-body correlations in the labile HB network. The microscopic event governing the kinetics of hydrophobic association has turned out to be a relatively large critical collective fluctuation in hydration water displacing a substantial fraction of HB clusters from the inner to the outer region of the first hydration shell.