Anomalous relaxation and dielectric response

Fractional electron transfer kinetics and a quantum breaking of ergodicity

The dissipative curve-crossing problem provides a paradigm for electron-transfer (ET) processes in condensed media. It establishes the simplest conceptual test bed to study the influence of the medium's dynamics on ET kinetics both on the ensemble level, and on the level of single particles. Single electron description is particularly important for nanoscaled systems like proteins, or molecular wires. Especially insightful is this framework in the semiclassical limit, where the environment can be treated classically, and an exact analytical treatment becomes feasible. Slow medium's dynamics is capable of enslaving ET and bringing it on the ensemble level from a quantum regime of nonadiabatic tunneling to the classical adiabatic regime, where electrons follow the nuclei rearrangements. This classical adiabatic textbook picture contradicts, however, in a very spectacular fashion to the statistics of single electron transitions, even in the Debye, memoryless media, also named Ohmic in the parlance of the famed spin-boson model. On the single particle level, ET always remains quantum, and this was named a quantum breaking of ergodicity in the adiabatic ET regime. What happens in the case of subdiffusive, fractional, or sub-Ohmic medium's dynamics, which is featured by power-law decaying dynamical memory effects typical, e.g., for protein macromolecules, and other viscoelastic media? Such a memory is vividly manifested by anomalous Cole-Cole dielectric response in such media. We address this question based both on accurate numerics and analytical theory. The ensemble theory remarkably agrees with the numerical dynamics of electronic populations, revealing a power-law relaxation tail even in a profoundly nonadiabatic electron transfer regime. In other words, ET in such media should typically display fractional kinetics. However, a profound difference with the numerically accurate results occurs for the distribution of residence times in the electronic states, both on the ensemble level and the level of single trajectories. Ergodicity is broken dynamically even in a more spectacular way than in the memoryless case. Our results question the applicability of all the existing and widely accepted ensemble theories of electron transfer in fractional, sub-Ohmic environments, on the level of single molecules, and provide a real challenge to face, both for theorists and experimentalists.

Viscoelastic subdiffusion in a random Gaussian environment

Viscoelastic subdiffusion governed by a fractional Langevin equation is studied numerically in a random Gaussian environment modeled by stationary Gaussian potentials with decaying spatial correlations. This anomalous diffusion is archetypal for living cells, where cytoplasm is known to be viscoelastic and a spatial disorder also naturally emerges. We obtain some first important insights into it within a model one-dimensional study. Two basic types of potential correlations are studied: short-range exponentially decaying and algebraically slow decaying with an infinite correlation length, both for a moderate (several k_BT, in the units of thermal energy), and strong (5-10k_BT) disorder. For a moderate disorder, it is shown that on the ensemble level viscoelastic subdiffusion can easily overcome the medium's disorder. Asymptotically, it is not distinguishable from the disorder-free subdiffusion. However, a strong scatter in single-trajectory averages is nevertheless seen even for a moderate disorder. It features a weak ergodicity breaking, which occurs on a very long yet transient time scale. Furthermore, for a strong disorder, a very long transient regime of logarithmic, Sinai-type diffusion emerges. It can last longer and be faster in the absolute terms for weakly decaying correlations as compared with the short-range correlations. Residence time distributions in a finite spatial domain are of a generalized log-normal type and are reminiscent also of a stretched exponential distribution. They can be easily confused for power-law distributions in view of the observed weak ergodicity breaking. This suggests a revision of some experimental data and their interpretation.

Sensing Magnetic Fields with Magnetosensitive Ion Channels

[-15]Magnetic nanoparticles are met across many biological species ranging from magnetosensitive bacteria, fishes, bees, bats, rats, birds, to humans. They can be both of biogenetic origin and due to environmental contamination, being either in paramagnetic or ferromagnetic state. The energy of such naturally occurring single-domain magnetic nanoparticles can reach up to 10-20 room k B T in the magnetic field of the Earth, which naturally led to supposition that they can serve as sensory elements in various animals. This work explores within a stochastic modeling framework a fascinating hypothesis of magnetosensitive ion channels with magnetic nanoparticles serving as sensory elements, especially, how realistic it is given a highly dissipative viscoelastic interior of living cells and typical sizes of nanoparticles possibly involved.

External dc bias-field effects in the nonlinear ac stationary response of dipolar particles in a mean-field potential

External dc bias-field effects on the nonlinear dielectric relaxation and dynamic Kerr effect of a system of permanent dipoles in a uniaxial mean-field potential are studied via the rotational Brownian motion model postulated in terms of the infinite hierarchy of differential-recurrence equations for the statistical moments f_{n}(t)=〈P_{n}〉(t) (the expectation value of the Legendre polynomials P_{n}). By solving these equations, the nonlinear dielectric and Kerr-effect ac stationary responses are evaluated for arbitrary dc field strength via perturbation theory in the ac field. Simple analytic equations based on the large separation of the time scales of the fast intrawell and slow overbarrier (interwell) relaxation processes are also derived.

Modeling magnetosensitive ion channels in the viscoelastic environment of living cells

We propose and study a model of hypothetical magnetosensitive ionic channels which are long thought to be a possible candidate to explain the influence of weak magnetic fields on living organisms ranging from magnetotactic bacteria to fishes, birds, rats, bats, and other mammals including humans. The core of the model is provided by a short chain of magnetosomes serving as a sensor, which is coupled by elastic linkers to the gating elements of ion channels forming a small cluster in the cell membrane. The magnetic sensor is fixed by one end on cytoskeleton elements attached to the membrane and is exposed to viscoelastic cytosol. Its free end can reorient stochastically and subdiffusively in viscoelastic cytosol responding to external magnetic field changes and can open the gates of coupled ion channels. The sensor dynamics is generally bistable due to bistability of the gates which can be in two states with probabilities which depend on the sensor orientation. For realistic parameters, it is shown that this model channel can operate in the magnetic field of Earth for a small number (five to seven) of single-domain magnetosomes constituting the sensor rod, each of which has a typical size found in magnetotactic bacteria and other organisms or even just one sufficiently large nanoparticle of a characteristic size also found in nature. It is shown that, due to the viscoelasticity of the medium, the bistable gating dynamics generally exhibits power law and stretched exponential distributions of the residence times of the channels in their open and closed states. This provides a generic physical mechanism for the explanation of the origin of such anomalous kinetics for other ionic channels whose sensors move in a viscoelastic environment provided by either cytosol or biological membrane, in a quite general context, beyond the fascinating hypothesis of magnetosensitive ionic channels we explore.

Cage effect for the velocity correlation functions of a Brownian particle in viscoelastic shear flows

The long-time limit behavior of velocity correlation functions (VCFs) for an underdamped Brownian particle in an oscillatory viscoelastic shear flow is investigated using the generalized Langevin equation with a power-law memory kernel. The influence of a fluctuating environment is modeled by an additive external fractional Gaussian noise. The exact expressions of the correlation functions of the fluctuating components of velocity for the Brownian particle in the shear plane have been calculated. Also, the particle's angular momentum is found. It is shown that in a certain region of the system parameters an interplay of the shear flow, memory effects, and external noise can induce a bounded long-time behavior of the VCFs, even in the shear flow direction, where in the case of internal noise the velocity process is subdiffusive, i.e., unbounded in time. Moreover, we find resonant behavior of the VCFs and the angular momentum versus the shear oscillation frequency, implying that they can be efficiently excited by oscillatory shear. The role of the initial positional distribution of particles is also discussed.

Anomalous transport in the crowded world of biological cells

A ubiquitous observation in cell biology is that the diffusive motion of macromolecules and organelles is anomalous, and a description simply based on the conventional diffusion equation with diffusion constants measured in dilute solution fails. This is commonly attributed to macromolecular crowding in the interior of cells and in cellular membranes, summarizing their densely packed and heterogeneous structures. The most familiar phenomenon is a sublinear, power-law increase of the mean-square displacement (MSD) as a function of the lag time, but there are other manifestations like strongly reduced and time-dependent diffusion coefficients, persistent correlations in time, non-Gaussian distributions of spatial displacements, heterogeneous diffusion and a fraction of immobile particles. After a general introduction to the statistical description of slow, anomalous transport, we summarize some widely used theoretical models: Gaussian models like fractional Brownian motion and Langevin equations for visco-elastic media, the continuous-time random walk model, and the Lorentz model describing obstructed transport in a heterogeneous environment. Particular emphasis is put on the spatio-temporal properties of the transport in terms of two-point correlation functions, dynamic scaling behaviour, and how the models are distinguished by their propagators even if the MSDs are identical. Then, we review the theory underlying commonly applied experimental techniques in the presence of anomalous transport like single-particle tracking, fluorescence correlation spectroscopy (FCS) and fluorescence recovery after photobleaching (FRAP). We report on the large body of recent experimental evidence for anomalous transport in crowded biological media: in cyto- and nucleoplasm as well as in cellular membranes, complemented by in vitro experiments where a variety of model systems mimic physiological crowding conditions. Finally, computer simulations are discussed which play an important role in testing the theoretical models and corroborating the experimental findings. The review is completed by a synthesis of the theoretical and experimental progress identifying open questions for future investigation.

Fractional-time random walk subdiffusion and anomalous transport with finite mean residence times: faster, not slower

Continuous time random walk (CTRW) subdiffusion along with the associated fractional Fokker-Planck equation (FFPE) is traditionally based on the premise of random clock with divergent mean period. This work considers an alternative CTRW and FFPE description which is featured by finite mean residence times (MRTs) in any spatial domain of finite size. Transient subdiffusive transport can occur on a very large time scale τ(c) which can greatly exceed mean residence time in any trap, τ(c) >>(τ), and even not being related to it. Asymptotically, on a macroscale transport becomes normal for t >> τ(c). However, mesoscopic transport is anomalous. Differently from viscoelastic subdiffusion no long-range anticorrelations among position increments are required. Moreover, our study makes it obvious that the transient subdiffusion and transport are faster than one expects from their normal asymptotic limit on a macroscale. This observation has profound implications for anomalous mesoscopic transport processes in biological cells because the macroscopic viscosity of cytoplasm is finite.

Fractional Brownian motors and stochastic resonance

We study fluctuating tilt Brownian ratchets based on fractional subdiffusion in sticky viscoelastic media characterized by a power law memory kernel. Unlike the normal diffusion case, the rectification effect vanishes in the adiabatically slow modulation limit and optimizes in a driving frequency range. It is shown also that the anomalous rectification effect is maximal (stochastic resonance effect) at optimal temperature and can be of surprisingly good quality. Moreover, subdiffusive current can flow in the counterintuitive direction upon a change of temperature or driving frequency. The dependence of anomalous transport on load exhibits a remarkably simple universality.

Generalized Langevin equation with multiplicative noise: temporal behavior of the autocorrelation functions

The temporal behavior of the mean-square displacement and the velocity autocorrelation function of a particle subjected to a periodic force in a harmonic potential well is investigated for viscoelastic media using the generalized Langevin equation. The interaction with fluctuations of environmental parameters is modeled by a multiplicative white noise, by an internal Mittag-Leffler noise with a finite memory time, and by an additive external noise. It is shown that the presence of a multiplicative noise has a profound effect on the behavior of the autocorrelation functions. Particularly, for correlation functions the model predicts a crossover between two different asymptotic power-law regimes. Moreover, a dependence of the correlation function on the frequency of the external periodic forcing occurs that gives a simple criterion to discern the multiplicative noise in future experiments. It is established that additive external and internal noises cause qualitatively different dependences of the autocorrelation functions on the external forcing and also on the time lag. The influence of the memory time of the internal noise on the dynamics of the system is also discussed.

Memory-enhanced energetic stability for a fractional oscillator with fluctuating frequency

The long-time limit behavior of the variance and the correlation function for the output signal of a fractional oscillator with fluctuating eigenfrequency subjected to a periodic force is considered. The influence of a fluctuating environment is modeled by a multiplicative white noise and by an additive noise with a zero mean. The viscoelastic-type friction kernel with memory is assumed as a power-law function of time. The exact expressions of stochastic resonance (SR) characteristics such as variance and signal-to-noise ratio (SNR) have been calculated. It is shown that at intermediate values of the memory exponent the energetic stability of the oscillator is significantly enhanced in comparison with the cases of strong and low memory. A multiresonancelike behavior of the variance and SNR as functions of the memory exponent is observed and a connection between this effect and the memory-induced enhancement of energetic stability is established. The effect of memory-induced energetic stability encountered in case the harmonic potential is absent, is also discussed.

Viscoelastic subdiffusion: from anomalous to normal

We study viscoelastic subdiffusion in bistable and periodic potentials within the generalized Langevin equation approach. Our results justify the (ultra)slow fluctuating rate view of the corresponding bistable non-Markovian dynamics which displays bursting and anticorrelation of the residence times in two potential wells. The transition kinetics is asymptotically stretched exponential when the potential barrier V0 several times exceeds thermal energy k(B)T [V(0) approximately (2-10)k(B)T] and it cannot be described by the non-Markovian rate theory (NMRT). The well-known NMRT result approximates, however, ever better with the increasing barrier height, the most probable logarithm of the residence times. Moreover, the rate description is gradually restored when the barrier height exceeds a fuzzy borderline which depends on the power-law exponent of free subdiffusion alpha . Such a potential-free subdiffusion is ergodic. Surprisingly, in periodic potentials it is not sensitive to the barrier height in the long time asymptotic limit. However, the transient to this asymptotic regime is extremally slow and it does profoundly depend on the barrier height. The time scale of such subdiffusion can exceed the mean residence time in a potential well or in a finite spatial domain by many orders of magnitude. All these features are in sharp contrast with an alternative subdiffusion mechanism involving jumps among traps with the divergent mean residence time in these traps.

Fractional Fokker-Planck subdiffusion in alternating force fields

The fractional Fokker-Planck equation for subdiffusion in time-dependent force fields is derived from the underlying continuous time random walk. Its limitations are discussed and it is then applied to the study of subdiffusion under the influence of a time-periodic rectangular force. As a main result, we show that such a force does not affect the universal scaling relation between the anomalous current and diffusion when applied to the biased dynamics: in the long-time limit, subdiffusion current and anomalous diffusion are immune to the driving. This is in sharp contrast with the unbiased case when the subdiffusion coefficient can be strongly enhanced, i.e., a zero-frequency response to a periodic driving is present.

Fractional Langevin equation: overdamped, underdamped, and critical behaviors

The dynamical phase diagram of the fractional Langevin equation is investigated for a harmonically bound particle. It is shown that critical exponents mark dynamical transitions in the behavior of the system. Four different critical exponents are found. (i) alpha_{c}=0.402+/-0.002 marks a transition to a nonmonotonic underdamped phase, (ii) alpha_{R}=0.441... marks a transition to a resonance phase when an external oscillating field drives the system, and (iii) alpha_{chi_{1}}=0.527... and (iv) alpha_{chi_{2}}=0.707... mark transitions to a double-peak phase of the "loss" when such an oscillating field present. As a physical explanation we present a cage effect, where the medium induces an elastic type of friction. Phase diagrams describing over and underdamped regimes, with or without resonances, show behaviors different from normal.

Critical exponent of the fractional Langevin equation

We investigate the dynamical phase diagram of the fractional Langevin equation and show that critical exponents mark dynamical transitions in the behavior of the system. For a free and harmonically bound particle the critical exponent alpha(c)=0.402+/-0.002 marks a transition to a nonmonotonic underdamped phase. The critical exponent alpha(R)=0.441... marks a transition to a resonance phase, when an external oscillating field drives the system. Physically, we explain these behaviors using a cage effect, where the medium induces an elastic type of friction. Phase diagrams describing the underdamped, the overdamped and critical frequencies of the fractional oscillator, recently used to model single protein experiments, show behaviors vastly different from normal.