Tomoshige Miyaguchi

Distributional ergodicity in stored-energy-driven Lévy flights.

We study a class of random walk, the stored-energy-driven Lévy flight (SEDLF), whose jump length is determined by a stored energy during a trapped state. The SEDLF is a continuous-time random walk with jump lengths being coupled with the trapping times. It is analytically shown that the ensemble-averaged mean-square displacements exhibit subdiffusion as well as superdiffusion, depending on the coupling parameter. We find that time-averaged mean-square displacements increase linearly with time and the diffusion coefficients are intrinsically random, a manifestation of distributional ergodicity. The diffusion coefficient shows aging in subdiffusive regime, whereas it increases with the measurement time in superdiffusive regime.

Singular Behavior of Slow Dynamics of Single Excitable Cells

In various kinds of cultured cells, it has been reported that the membrane potential exhibits fluctuations with long-term correlations, although the underlying mechanism remains to be elucidated. A cardiac muscle cell culture serves as an excellent experimental system to investigate this phenomenon because timings of excitations can be determined over an extended time period in a noninvasive manner through visualization of contractions, although the properties of beat-timing fluctuations of cardiac muscle cells at the single-cell level remains to be fully clarified. In this article, we report on our investigation of spontaneous contractions of cultured rat cardiac muscle cells at the single-cell level. It was found that single cells exhibit several typical temporal patterns of contractions and spontaneous transitions among them. Detrended fluctuation analysis on the time series of interbeat intervals revealed the presence of 1/f(beta) noise at sufficiently large timescales. Furthermore, multifractality was also found in the time series of interbeat intervals. These experimental trends were successfully explained using a simple mathematical model, incorporating correlated noise into ionic currents. From these findings, it was established that singular fluctuations accompanying 1/f(beta) noise and multifractality are intrinsic properties of single cardiac muscle cells.

Fluctuation analysis of time-averaged mean-square displacement for the Langevin equation with time-dependent and fluctuating diffusivity.

The mean-square displacement (MSD) is widely utilized to study the dynamical properties of stochastic processes. The time-averaged MSD (TAMSD) provides some information on the dynamics which cannot be extracted from the ensemble-averaged MSD. In particular, the relative standard deviation (RSD) of the TAMSD can be utilized to study the long-time relaxation behavior. In this work, we consider a class of Langevin equations which are multiplicatively coupled to time-dependent and fluctuating diffusivities. Various interesting dynamics models such as entangled polymers and supercooled liquids can be interpreted as the Langevin equations with time-dependent and fluctuating diffusivities. We derive a general formula for the RSD of the TAMSD for the Langevin equation with the time-dependent and fluctuating diffusivity. We show that the RSD can be expressed in terms of the correlation function of the diffusivity. The RSD exhibits the crossover at the long time region. The crossover time is related to a weighted average relaxation time for the diffusivity. Thus the crossover time gives some information on the relaxation time of fluctuating diffusivity which cannot be extracted from the ensemble-averaged MSD. We discuss the universality and possible applications of the formula via some simple examples.

Ergodic properties of continuous-time random walks: Finite-size effects and ensemble dependences

The effects of spatial confinements and smooth cutoffs of the waiting time distribution in continuous-time random walks (CTRWs) are studied analytically. We also investigate dependences of ergodic properties on initial ensembles (i.e., distributions of the first waiting time). Here, we consider two ensembles: the equilibrium and a typical non-equilibrium ensembles. For both ensembles, it is shown that the time-averaged mean square displacement (TAMSD) exhibits a crossover from normal to anomalous diffusion due to the spacial confinement and this crossover does not vanish even in the long measurement time limit. Moreover, for the non-equilibrium ensemble, we show that the probability density function of the diffusion constant of TAMSD follows the transient Mittag-Leffler distribution, and that scatter in the TAMSD shows a clear transition from weak ergodicity breaking (an irreproducible regime) to ordinary ergodic behavior (a reproducible regime) as the measurement time increases. This convergence to ordinary ergodicity requires a long measurement time compared to common distributions such as the exponential distribution; in other words, the weak ergodicity breaking persists for a long time. In addition, it is shown that, besides the TAMSD, a class of observables also exhibits this slow convergence to ergodicity. We also point out that, even though the system with the equilibrium initial ensemble shows no aging, its behavior is quite similar to that for the non-equilibrium ensemble.

Anomalous diffusion in a quenched-trap model on fractal lattices.

Models with mixed origins of anomalous subdiffusion have been considered important for understanding transport in biological systems. Here one such mixed model, the quenched-trap model (QTM) on fractal lattices, is investigated. It is shown that both ensemble- and time-averaged mean-square displacements (MSDs) show subdiffusion with different scaling exponents, i.e., this system shows weak ergodicity breaking. Moreover, time-averaged MSD exhibits aging and converges to a random variable following the modified Mittag-Leffler distribution. It is also shown that the QTM on a fractal lattice cannot be reduced to the continuous-time random walks if the spectral dimension of the fractal lattice is less than 2.

Langevin equation with fluctuating diffusivity: A two-state model

Recently, anomalous subdiffusion, aging, and scatter of the diffusion coefficient have been reported in many single-particle-tracking experiments, though the origins of these behaviors are still elusive. Here, as a model to describe such phenomena, we investigate a Langevin equation with diffusivity fluctuating between a fast and a slow state. Namely, the diffusivity follows a dichotomous stochastic process. We assume that the sojourn time distributions of these two states are given by power laws. It is shown that, for a nonequilibrium ensemble, the ensemble-averaged mean-square displacement (MSD) shows transient subdiffusion. In contrast, the time-averaged MSD shows normal diffusion, but an effective diffusion coefficient transiently shows aging behavior. The propagator is non-Gaussian for short time and converges to a Gaussian distribution in a long-time limit; this convergence to Gaussian is extremely slow for some parameter values. For equilibrium ensembles, both ensemble-averaged and time-averaged MSDs show only normal diffusion and thus we cannot detect any traces of the fluctuating diffusivity with these MSDs. Therefore, as an alternative approach to characterizing the fluctuating diffusivity, the relative standard deviation (RSD) of the time-averaged MSD is utilized and it is shown that the RSD exhibits slow relaxation as a signature of the long-time correlation in the fluctuating diffusivity. Furthermore, it is shown that the RSD is related to a non-Gaussian parameter of the propagator. To obtain these theoretical results, we develop a two-state renewal theory as an analytical tool.

Phase Diagram in Stored-Energy-Driven Lévy Flight

Phase diagram based on the mean square displacement (MSD) and the distribution of diffusion coefficients of the time-averaged MSD for the stored-energy-driven Lévy flight (SEDLF) is presented. In the SEDLF, a random walker cannot move while storing energy, and it jumps by the stored energy. The SEDLF shows a whole spectrum of anomalous diffusions including subdiffusion and superdiffusion, depending on the coupling parameter between storing time (trapping time) and stored energy. This stochastic process can be investigated analytically with the aid of renewal theory. Here, we consider two different renewal processes, i.e., ordinary renewal process and equilibrium renewal process, when the mean trapping time does not diverge. We analytically show the phase diagram according to the coupling parameter and the power exponent in the trapping-time distribution. In particular, we find that distributional behavior of time-averaged MSD intrinsically appears in superdiffusive as well as normal diffusive regime even when the mean trapping time does not diverge.

Ultraslow convergence to ergodicity in transient subdiffusion.

We investigate continuous time random walks with truncated α-stable trapping times. We prove distributional ergodicity for a class of observables; namely, the time-averaged observables follow the probability density function called the Mittag-Leffler distribution. This distributional ergodic behavior persists for a long time, and thus the convergence to the ordinary ergodicity is considerably slower than in the case in which the trapping-time distribution is given by common distributions. We also find a crossover from the distributional ergodic behavior to the ordinary ergodic behavior.

Crossover time in relative fluctuations characterizes the longest relaxation time of entangled polymers

In entangled polymer systems, there are several characteristic time scales, such as the entanglement time and the disengagement time. In molecular simulations, the longest relaxation time (the disengagement time) can be determined by the mean square displacement (MSD) of a segment or by the shear relaxation modulus. Here, we propose the relative fluctuation analysis method, which is originally developed for characterizing large fluctuations, to determine the longest relaxation time from the center of mass trajectories of polymer chains (the time-averaged MSDs). Applying the method to simulation data of entangled polymers (by the slip-spring model and the simple reptation model), we provide a clear evidence that the longest relaxation time is estimated as the crossover time in the relative fluctuations.

Elucidating fluctuating diffusivity in center-of-mass motion of polymer models with time-averaged mean-square-displacement tensor

There have been increasing reports that the diffusion coefficient of macromolecules depends on time and fluctuates randomly. Here a method is developed to elucidate this fluctuating diffusivity from trajectory data. Time-averaged mean-square displacement (MSD), a common tool in single-particle-tracking (SPT) experiments, is generalized to a second-order tensor with which both magnitude and orientation fluctuations of the diffusivity can be clearly detected. This method is used to analyze the center-of-mass motion of four fundamental polymer models: the Rouse model, the Zimm model, a reptation model, and a rigid rodlike polymer. It is found that these models exhibit distinctly different types of magnitude and orientation fluctuations of diffusivity. This is an advantage of the present method over previous ones, such as the ergodicity-breaking parameter and a non-Gaussian parameter, because with either of these parameters it is difficult to distinguish the dynamics of the four polymer models. Also, the present method of a time-averaged MSD tensor could be used to analyze trajectory data obtained in SPT experiments.