@article{KharchenkoGoychuk2013,
  author    = {Kharchenko, Vasyl O. and Goychuk, Igor},
  title     = {Subdiffusive rocking ratchets in viscoelastic media transport optimization and thermodynamic efficiency in overdamped regime},
  series = {Physical review : E, Statistical, nonlinear and soft matter physics},
  volume    = {87},
  journal   = {Physical review : E, Statistical, nonlinear and soft matter physics},
  number    = {5},
  publisher = {American Physical Society},
  address   = {College Park},
  issn      = {1539-3755},
  doi       = {10.1103/PhysRevE.87.052119},
  pages     = {9},
  year      = {2013},
  abstract  = {We study subdiffusive overdamped Brownian ratchets periodically rocked by an external zero-mean force in viscoelastic media within the framework of a non-Markovian generalized Langevin equation approach and associated multidimensional Markovian embedding dynamics. Viscoelastic deformations of the medium caused by the transport particle are modeled by a set of auxiliary Brownian quasiparticles elastically coupled to the transport particle and characterized by a hierarchy of relaxation times which obey a fractal scaling. The most slowly relaxing deformations which cannot immediately follow to the moving particle imprint long-range memory about its previous positions and cause subdiffusion and anomalous transport on a sufficiently long time scale. This anomalous behavior is combined with normal diffusion and transport on an initial time scale of overdamped motion. Anomalously slow directed transport in a periodic ratchet potential with broken space inversion symmetry emerges due to a violation of the thermal detailed balance by a zero-mean periodic driving and is optimized with frequency of driving, its amplitude, and temperature. Such optimized anomalous transport can be low dispersive and characterized by a large generalized Peclet number. Moreover, we show that overdamped subdiffusive ratchets can sustain a substantial load and do useful work. The corresponding thermodynamic efficiency decays algebraically in time since the useful work done against a load scales sublinearly with time following to the transport particle position, but the energy pumped by an external force scales with time linearly. Nevertheless, it can be transiently appreciably high and compare well with the thermodynamical efficiency of the normal diffusion overdamped ratchets on sufficiently long temporal and spatial scales.},
  language  = {en}
}
@article{GoychukKharchenkoMetzler2014,
  author    = {Goychuk, Igor and Kharchenko, Vasyl O. and Metzler, Ralf},
  title     = {How Molecular Motors Work in the Crowded Environment of Living Cells: Coexistence and Efficiency of Normal and Anomalous Transport},
  series = {PLoS one},
  volume    = {9},
  journal   = {PLoS one},
  number    = {3},
  publisher = {PLoS},
  address   = {San Fransisco},
  issn      = {1932-6203},
  doi       = {10.1371/journal.pone.0091700},
  pages     = {7},
  year      = {2014},
  abstract  = {Recent experiments reveal both passive subdiffusion of various nanoparticles and anomalous active transport of such particles by molecular motors in the molecularly crowded environment of living biological cells. Passive and active microrheology reveals that the origin of this anomalous dynamics is due to the viscoelasticity of the intracellular fluid. How do molecular motors perform in such a highly viscous, dissipative environment? Can we explain the observed co-existence of the anomalous transport of relatively large particles of 100 to 500 nm in size by kinesin motors with the normal transport of smaller particles by the same molecular motors? What is the efficiency of molecular motors in the anomalous transport regime? Here we answer these seemingly conflicting questions and consistently explain experimental findings in a generalization of the well-known continuous diffusion model for molecular motors with two conformational states in which viscoelastic effects are included.},
  language  = {en}
}
@article{GoychukKharchenkoMetzler2014,
  author    = {Goychuk, Igor and Kharchenko, Vasyl O. and Metzler, Ralf},
  title     = {Molecular motors pulling cargos in the viscoelastic cytosol: how power strokes beat subdiffusion},
  series = {Physical chemistry, chemical physics : a journal of European Chemical Societies},
  volume    = {16},
  journal   = {Physical chemistry, chemical physics : a journal of European Chemical Societies},
  number    = {31},
  publisher = {Royal Society of Chemistry},
  address   = {Cambridge},
  issn      = {1463-9076},
  doi       = {10.1039/c4cp01234h},
  pages     = {16524 -- 16535},
  year      = {2014},
  abstract  = {The discovery of anomalous diffusion of larger biopolymers and submicron tracers such as endogenous granules, organelles, or virus capsids in living cells, attributed to the viscoelastic nature of the cytoplasm, provokes the question whether this complex environment equally impacts the active intracellular transport of submicron cargos by molecular motors such as kinesins: does the passive anomalous diffusion of free cargo always imply its anomalously slow active transport by motors, the mean transport distance along microtubule growing sublinearly rather than linearly in time? Here we analyze this question within the widely used two-state Brownian ratchet model of kinesin motors based on the continuous-state diffusion along microtubules driven by a flashing binding potential, where the cargo particle is elastically attached to the motor. Depending on the cargo size, the loading force, the amplitude of the binding potential, the turnover frequency of the molecular motor enzyme, and the linker stiffness we demonstrate that the motor transport may turn out either normal or anomalous, as indeed measured experimentally. We show how a highly efficient normal active transport mediated by motors may emerge despite the passive anomalous diffusion of the cargo, and study the intricate effects of the elastic linker. Under different, well specified conditions the microtubule-based motor transport becomes anomalously slow and thus significantly less efficient.},
  language  = {en}
}
@article{GoychukKharchenko2014,
  author    = {Goychuk, Igor and Kharchenko, Vasyl O.},
  title     = {Anomalous features of diffusion in corrugated potentials with spatial correlations: Faster than normal, and other surprises},
  series = {Physical review letters},
  volume    = {113},
  journal   = {Physical review letters},
  number    = {10},
  publisher = {American Physical Society},
  address   = {College Park},
  issn      = {0031-9007},
  doi       = {10.1103/PhysRevLett.113.100601},
  pages     = {5},
  year      = {2014},
  abstract  = {Normal diffusion in corrugated potentials with spatially uncorrelated Gaussian energy disorder famously explains the origin of non-Arrhenius exp[-sigma(2)/(k(B)T(2))] temperature dependence in disordered systems. Here we show that unbiased diffusion remains asymptotically normal also in the presence of spatial correlations decaying to zero. However, because of a temporal lack of self-averaging, transient subdiffusion emerges on the mesoscale, and it can readily reach macroscale even for moderately strong disorder fluctuations of sigma similar to 4 - 5k(B)T. Because of its nonergodic origin, such subdiffusion exhibits a large scatter in single-trajectory averages. However, at odds with intuition, it occurs essentially faster than one expects from the normal diffusion in the absence of correlations. We apply these results to diffusion of regulatory proteins on DNA molecules and predict that such diffusion should be anomalous, but much faster than earlier expected on a typical length of genes for a realistic energy disorder of several room k(B)T, or merely 0.05-0.075 eV.},
  language  = {en}
}
@misc{GoychukKharchenko2013,
  author    = {Goychuk, Igor and Kharchenko, Vasyl O.},
  title     = {Rocking subdiffusive ratchets},
  series = {Mathematical Modelling of Natural Phenomena},
  journal   = {Mathematical Modelling of Natural Phenomena},
  number    = {622},
  issn      = {1866-8372},
  doi       = {10.1051/mmnp/20138210},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-416138},
  pages     = {15},
  year      = {2013},
  abstract  = {We study origin, parameter optimization, and thermodynamic efficiency of isothermal rocking ratchets based on fractional subdiffusion within a generalized non-Markovian Langevin equation approach. A corresponding multi-dimensional Markovian embedding dynamics is realized using a set of auxiliary Brownian particles elastically coupled to the central Brownian particle (see video on the journal web site). We show that anomalous subdiffusive transport emerges due to an interplay of nonlinear response and viscoelastic effects for fractional Brownian motion in periodic potentials with broken space-inversion symmetry and driven by a time-periodic field. The anomalous transport becomes optimal for a subthreshold driving when the driving period matches a characteristic time scale of interwell transitions. It can also be optimized by varying temperature, amplitude of periodic potential and driving strength. The useful work done against a load shows a parabolic dependence on the load strength. It grows sublinearly with time and the corresponding thermodynamic efficiency decays algebraically in time because the energy supplied by the driving field scales with time linearly. However, it compares well with the efficiency of normal diffusion rocking ratchets on an appreciably long time scale.},
  language  = {en}
}
@article{GoychukGoychuk2015,
  author    = {Goychuk, Igor and Goychuk, Andriy},
  title     = {Stochastic Wilson-Cowan models of neuronal network dynamics with memory and delay},
  series = {New journal of physics : the open-access journal for physics},
  volume    = {17},
  journal   = {New journal of physics : the open-access journal for physics},
  publisher = {IOP Publ. Ltd.},
  address   = {Bristol},
  issn      = {1367-2630},
  doi       = {10.1088/1367-2630/17/4/045029},
  pages     = {17},
  year      = {2015},
  abstract  = {We consider a simple Markovian class of the stochastic Wilson-Cowan type models of neuronal network dynamics, which incorporates stochastic delay caused by the existence of a refractory period of neurons. From the point of view of the dynamics of the individual elements, we are dealing with a network of non-Markovian stochastic two-state oscillators with memory, which are coupled globally in a mean-field fashion. This interrelation of a higher-dimensional Markovian and lower-dimensional non-Markovian dynamics is discussed in its relevance to the general problem of the network dynamics of complex elements possessing memory. The simplest model of this class is provided by a three-state Markovian neuron with one refractory state, which causes firing delay with an exponentially decaying memory within the two-state reduced model. This basic model is used to study critical avalanche dynamics (the noise sustained criticality) in a balanced feedforward network consisting of the excitatory and inhibitory neurons. Such avalanches emerge due to the network size dependent noise (mesoscopic noise). Numerical simulations reveal an intermediate power law in the distribution of avalanche sizes with the critical exponent around -1.16. We show that this power law is robust upon a variation of the refractory time over several orders of magnitude. However, the avalanche time distribution is biexponential. It does not reflect any genuine power law dependence.},
  language  = {en}
}
@article{GoychukGoychuk2015,
  author    = {Goychuk, Igor and Goychuk, Andriy},
  title     = {Stochastic Wilson},
  series = {New journal of physics},
  volume    = {17},
  journal   = {New journal of physics},
  number    = {4},
  publisher = {Deutsche Physikalische Gesellschaft, Institute of Physics},
  address   = {Bad Honnef, London},
  issn      = {1367-2630},
  doi       = {10.1088/1367-2630/17/4/045029},
  year      = {2015},
  abstract  = {We consider a simple Markovian class of the stochastic Wilson-Cowan type models of neuronal network dynamics, which incorporates stochastic delay caused by the existence of a refractory period of neurons. From the point of view of the dynamics of the individual elements, we are dealing with a network of non-Markovian stochastic two-state oscillators with memory, which are coupled globally in a mean-field fashion. This interrelation of a higher-dimensional Markovian and lower-dimensional non-Markovian dynamics is discussed in its relevance to the general problem of the network dynamics of complex elements possessing memory. The simplest model of this class is provided by a three-state Markovian neuron with one refractory state, which causes firing delay with an exponentially decaying memory within the two-state reduced model. This basic model is used to study critical avalanche dynamics (the noise sustained criticality) in a balanced feedforward network consisting of the excitatory and inhibitory neurons. Such avalanches emerge due to the network size dependent noise (mesoscopic noise). Numerical simulations reveal an intermediate power law in the distribution of avalanche sizes with the critical exponent around -1.16. We show that this power law is robust upon a variation of the refractory time over several orders of magnitude. However, the avalanche time distribution is biexponential. It does not reflect any genuine power law dependence.},
  language  = {en}
}
@misc{GoychukGoychuk2015,
  author    = {Goychuk, Igor and Goychuk, Andriy},
  title     = {Stochastic Wilson},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-77814},
  year      = {2015},
  abstract  = {We consider a simple Markovian class of the stochastic Wilson-Cowan type models of neuronal network dynamics, which incorporates stochastic delay caused by the existence of a refractory period of neurons. From the point of view of the dynamics of the individual elements, we are dealing with a network of non-Markovian stochastic two-state oscillators with memory, which are coupled globally in a mean-field fashion. This interrelation of a higher-dimensional Markovian and lower-dimensional non-Markovian dynamics is discussed in its relevance to the general problem of the network dynamics of complex elements possessing memory. The simplest model of this class is provided by a three-state Markovian neuron with one refractory state, which causes firing delay with an exponentially decaying memory within the two-state reduced model. This basic model is used to study critical avalanche dynamics (the noise sustained criticality) in a balanced feedforward network consisting of the excitatory and inhibitory neurons. Such avalanches emerge due to the network size dependent noise (mesoscopic noise). Numerical simulations reveal an intermediate power law in the distribution of avalanche sizes with the critical exponent around -1.16. We show that this power law is robust upon a variation of the refractory time over several orders of magnitude. However, the avalanche time distribution is biexponential. It does not reflect any genuine power law dependence.},
  language  = {en}
}
@article{Goychuk2018,
  author    = {Goychuk, Igor},
  title     = {Viscoelastic subdiffusion in a random Gaussian environment},
  series = {Physical chemistry, chemical physics : a journal of European Chemical Societies},
  volume    = {20},
  journal   = {Physical chemistry, chemical physics : a journal of European Chemical Societies},
  number    = {37},
  publisher = {Royal Society of Chemistry},
  address   = {Cambridge},
  issn      = {1463-9076},
  doi       = {10.1039/c8cp05238g},
  pages     = {24140 -- 24155},
  year      = {2018},
  abstract  = {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 kBT, in the units of thermal energy), and strong (5-10kBT) 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.},
  language  = {en}
}
@article{Goychuk2017,
  author    = {Goychuk, Igor},
  title     = {Fractional Bhatnagar-Gross-Krook kinetic equation},
  series = {The European physical journal : B, Condensed matter and complex systems},
  volume    = {90},
  journal   = {The European physical journal : B, Condensed matter and complex systems},
  publisher = {Springer},
  address   = {New York},
  issn      = {1434-6028},
  doi       = {10.1140/epjb/e2017-80297-x},
  pages     = {13},
  year      = {2017},
  abstract  = {The linear Boltzmann equation approach is generalized to describe fractional superdiffusive transport of the Levy walk type in external force fields. The time distribution between scattering events is assumed to have a finite mean value and infinite variance. It is completely characterized by the two scattering rates, one fractional and a normal one, which defines also the mean scattering rate. We formulate a general fractional linear Boltzmann equation approach and exemplify it with a particularly simple case of the Bohm and Gross scattering integral leading to a fractional generalization of the Bhatnagar, Gross and Krook kinetic equation. Here, at each scattering event the particle velocity is completely randomized and takes a value from equilibrium Maxwell distribution at a given fixed temperature. We show that the retardation effects are indispensable even in the limit of infinite mean scattering rate and argue that this novel fractional kinetic equation provides a viable alternative to the fractional Kramers-Fokker-Planck (KFP) equation by Barkai and Silbey and its generalization by Friedrich et al. based on the picture of divergent mean time between scattering events. The case of divergent mean time is also discussed at length and compared with the earlier results obtained within the fractional KFP. Also a phenomenological fractional BGK equation without retardation effects is proposed in the limit of infinite scattering rates. It cannot be, however, rigorously derived from a scattering model, being rather clever postulated. It this respect, this retardationless equation is similar to the fractional KFP by Barkai and Silbey. However, it corresponds to the opposite, much more physical limit and, therefore, also presents a viable alternative.},
  language  = {en}
}
@misc{Goychuk2019,
  author    = {Goychuk, Igor},
  title     = {Comment on "Anomalous Escape Governed by Thermal 1/f Noise" Reply (R. K. Singh)},
  series = {Physical review letters},
  volume    = {123},
  journal   = {Physical review letters},
  number    = {23},
  publisher = {American Physical Society},
  address   = {College Park},
  issn      = {0031-9007},
  doi       = {10.1103/PhysRevLett.123.238902},
  pages     = {1},
  year      = {2019},
  language  = {en}
}
@article{Goychuk2018,
  author    = {Goychuk, Igor},
  title     = {Sensing magnetic fields with magnetosensitive ion channels},
  series = {Sensors},
  volume    = {18},
  journal   = {Sensors},
  number    = {3},
  publisher = {MDPI},
  address   = {Basel},
  issn      = {1424-8220},
  doi       = {10.3390/s18030728},
  pages     = {19},
  year      = {2018},
  abstract  = {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.},
  language  = {en}
}
@misc{Goychuk2016,
  author    = {Goychuk, Igor},
  title     = {Quantum ergodicity breaking in semi-classical electron transfer dynamics},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-102264},
  pages     = {11},
  year      = {2016},
  abstract  = {Can the statistical properties of single-electron transfer events be correctly predicted within a common equilibrium ensemble description? This fundamental in nanoworld question of ergodic behavior is scrutinized within a very basic semi-classical curve-crossing problem. It is shown that in the limit of non-adiabatic electron transfer (weak tunneling) well-described by the Marcus-Levich-Dogonadze(MLD) rate the answer is yes. However, in the limit of the so-called solvent-controlled adiabatic electron transfer, a profound breaking of ergodicity occurs. Namely, a common description based on the ensemble reduced density matrix with an initial equilibrium distribution of the reaction coordinate is not able to reproduce the statistics of single-trajectory events in this seemingly classical regime. For sufficiently large activation barriers, the ensemble survival probability in a state remains nearly exponential with the inverse rate given by the sum of the adiabatic curve crossing (Kramers) time and the inverse MLD rate. In contrast, near to the adiabatic regime, the single-electron survival probability is clearly non-exponential, even though it possesses an exponential tail which agrees well with the ensemble description. Initially, it is well described by a Mittag-Leffler distribution with a fractional rate. Paradoxically, the mean transfer time in this classical on the ensemble level regime is well described by the inverse of the nonadiabatic quantum tunneling rate on a single particle level. An analytical theory is developed which perfectly agrees with stochastic simulations and explains our findings.},
  language  = {en}
}
@article{Goychuk2017,
  author    = {Goychuk, Igor},
  title     = {Quantum ergodicity breaking in semi-classical electron transfer dynamics},
  series = {Physical chemistry, chemical physics : a journal of European Chemical Societies},
  volume    = {19},
  journal   = {Physical chemistry, chemical physics : a journal of European Chemical Societies},
  publisher = {Royal Society of Chemistry},
  address   = {Cambridge},
  issn      = {1463-9076},
  doi       = {10.1039/c6cp07206b},
  pages     = {3056 -- 3066},
  year      = {2017},
  abstract  = {Can the statistical properties of single-electron transfer events be correctly predicted within a common equilibrium ensemble description? This fundamental in nanoworld question of ergodic behavior is scrutinized within a very basic semi-classical curve-crossing problem. It is shown that in the limit of non-adiabatic electron transfer (weak tunneling) well-described by the Marcus-Levich-Dogonadze (MLD) rate the answer is yes. However, in the limit of the so-called solvent-controlled adiabatic electron transfer, a profound breaking of ergodicity occurs. Namely, a common description based on the ensemble reduced density matrix with an initial equilibrium distribution of the reaction coordinate is not able to reproduce the statistics of single-trajectory events in this seemingly classical regime. For sufficiently large activation barriers, the ensemble survival probability in a state remains nearly exponential with the inverse rate given by the sum of the adiabatic curve crossing (Kramers) time and the inverse MLD rate. In contrast, near to the adiabatic regime, the single-electron survival probability is clearly non-exponential, even though it possesses an exponential tail which agrees well with the ensemble description. Initially, it is well described by a Mittag-Leffler distribution with a fractional rate. Paradoxically, the mean transfer time in this classical on the ensemble level regime is well described by the inverse of the nonadiabatic quantum tunneling rate on a single particle level. An analytical theory is developed which perfectly agrees with stochastic simulations and explains our findings.},
  language  = {en}
}
@article{Goychuk2019,
  author    = {Goychuk, Igor},
  title     = {Fractional electron transfer kinetics and a quantum breaking of ergodicity},
  series = {Physical review : E, Statistical, nonlinear and soft matter physics},
  volume    = {99},
  journal   = {Physical review : E, Statistical, nonlinear and soft matter physics},
  number    = {5},
  publisher = {American Physical Society},
  address   = {College Park},
  issn      = {2470-0045},
  doi       = {10.1103/PhysRevE.99.052136},
  pages     = {21},
  year      = {2019},
  abstract  = {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.},
  language  = {en}
}
@article{Goychuk2019,
  author    = {Goychuk, Igor},
  title     = {Fractional Hydrodynamic Memory and Superdiffusion in Tilted Washboard Potentials},
  series = {Physical review letters},
  volume    = {123},
  journal   = {Physical review letters},
  number    = {18},
  publisher = {American Physical Society},
  address   = {College Park},
  issn      = {0031-9007},
  doi       = {10.1103/PhysRevLett.123.180603},
  pages     = {6},
  year      = {2019},
  abstract  = {Diffusion in tilted washboard potentials can paradoxically exceed free normal diffusion. The effect becomes much stronger in the underdamped case due to inertial effects. What happens upon inclusion of usually neglected fractional hydrodynamics memory effects (Basset-Boussinesq frictional force), which result in a heavy algebraic tail of the velocity autocorrelation function of the potential-free diffusion making it transiently superdiffusive? Will a giant enhancement of diffusion become even stronger, and the transient superdiffusion last even longer? These are the questions that we answer in this Letter based on an accurate numerical investigation. We show that a resonancelike enhancement of normal diffusion becomes indeed much stronger and sharper. Moreover, a long-lasting transient regime of superdiffusion, including Richardson-like diffusion, <delta x(2)(t)>proportional to t(3) and ballistic supertransport, <delta x(t)>proportional to t(2), is revealed.},
  language  = {en}
}
@article{Goychuk2019,
  author    = {Goychuk, Igor},
  title     = {Perfect anomalous transport of subdiffusive cargos by molecular motors in viscoelastic cytosol},
  series = {Biosystems : journal of biological and information processing sciences},
  volume    = {177},
  journal   = {Biosystems : journal of biological and information processing sciences},
  publisher = {Elsevier},
  address   = {Oxford},
  issn      = {0303-2647},
  doi       = {10.1016/j.biosystems.2018.11.004},
  pages     = {56 -- 65},
  year      = {2019},
  abstract  = {Multiple experiments show that various submicron particles such as magnetosomes, RNA messengers, viruses, and even much smaller nanoparticles such as globular proteins diffuse anomalously slow in viscoelastic cytosol of living cells. Hence, their sufficiently fast directional transport by molecular motors such as kinesins is crucial for the cell operation. It has been shown recently that the traditional flashing Brownian ratchet models of molecular motors are capable to describe both normal and anomalous transport of such subdiffusing cargos by molecular motors with a very high efficiency. This work elucidates further an important role of mechanochemical coupling in such an anomalous transport. It shows a natural emergence of a perfect subdiffusive ratchet regime due to allosteric effects, where the random rotations of a "catalytic wheel" at the heart of the motor operation become perfectly synchronized with the random stepping of a heavily loaded motor, so that only one ATP molecule is consumed on average at each motor step along microtubule. However, the number of rotations made by the catalytic engine and the traveling distance both scale sublinearly in time. Nevertheless, this anomalous transport can be very fast in absolute terms.},
  language  = {en}
}
@article{Goychuk2014,
  author    = {Goychuk, Igor},
  title     = {Life and death of stationary linear response in anomalous continuous time random walk dynamics},
  series = {Communications in theoretical physics : a series journal of the Chinese Physical Society (A)},
  volume    = {62},
  journal   = {Communications in theoretical physics : a series journal of the Chinese Physical Society (A)},
  number    = {4},
  publisher = {IOP Publ. Ltd.},
  address   = {Bristol},
  issn      = {0253-6102},
  pages     = {497 -- 504},
  year      = {2014},
  abstract  = {Linear theory of stationary response in systems at thermal equilibrium requires to find equilibrium correlation function of unperturbed responding system. Studies of the response of the systems exhibiting anomalously slow dynamics are often based on the continuous time random walk description (CTRW) with divergent mean waiting times. The bulk of the literature on anomalous response contains linear response functions like one by Cole-Cole calculated from such a CTRW theory and applied to systems at thermal equilibrium. Here we show within a fairly simple and general model that for the systems with divergent mean waiting times the stationary response at thermal equilibrium is absent, in accordance with some recent studies. The absence of such stationary response (or dying to zero non-stationary response in aging experiments) would confirm CTRW with divergent mean waiting times as underlying physical relaxation mechanism, but reject it otherwise. We show that the absence of stationary response is closely related to the breaking of ergodicity of the corresponding dynamical variable. As an important new result, we derive a generalized Cole-Cole response within ergodic CTRW dynamics with finite waiting time. Moreover, we provide a physically reasonable explanation of the origin and wide presence of 1/f noise in condensed matter for ergodic dynamics close to normal, rather than strongly deviating.},
  language  = {en}
}
@article{Goychuk2015,
  author    = {Goychuk, Igor},
  title     = {Modeling magnetosensitive ion channels in the viscoelastic environment of living cells},
  series = {Physical review : E, Statistical, nonlinear and soft matter physics},
  volume    = {92},
  journal   = {Physical review : E, Statistical, nonlinear and soft matter physics},
  number    = {4},
  publisher = {American Physical Society},
  address   = {College Park},
  issn      = {1539-3755},
  doi       = {10.1103/PhysRevE.92.042711},
  pages     = {18},
  year      = {2015},
  abstract  = {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.},
  language  = {en}
}
@misc{Goychuk2016,
  author    = {Goychuk, Igor},
  title     = {Molecular machines operating on the nanoscale: from classical to quantum},
  series = {Beilstein journal of nanotechnology},
  volume    = {7},
  journal   = {Beilstein journal of nanotechnology},
  publisher = {Beilstein-Institut zur F{\~A}\Prderung der Chemischen Wissenschaften},
  address   = {Frankfurt, Main},
  issn      = {2190-4286},
  doi       = {10.3762/bjnano.7.31},
  pages     = {328 -- 350},
  year      = {2016},
  abstract  = {The main physical features and operating principles of isothermal nanomachines in the microworld, common to both classical and quantum machines, are reviewed. Special attention is paid to the dual, constructive role of dissipation and thermal fluctuations, the fluctuation-dissipation theorem, heat losses and free energy transduction, thermodynamic efficiency, and thermodynamic efficiency at maximum power. Several basic models are considered and discussed to highlight generic physical features. This work examines some common fallacies that continue to plague the literature. In particular, the erroneous beliefs that one should minimize friction and lower the temperature for high performance of Brownian machines, and that the thermodynamic efficiency at maximum power cannot exceed one-half are discussed. The emerging topic of anomalous molecular motors operating subdiffusively but very efficiently in the viscoelastic environment of living cells is also discussed.},
  language  = {en}
}
@article{Goychuk2015,
  author    = {Goychuk, Igor},
  title     = {Anomalous transport of subdiffusing cargos by single kinesin motors: the role of mechano-chemical coupling and anharmonicity of tether},
  series = {Physical biology : a journal for the fundamental understanding of biological systems},
  volume    = {12},
  journal   = {Physical biology : a journal for the fundamental understanding of biological systems},
  number    = {1},
  publisher = {IOP Publ. Ltd.},
  address   = {Bristol},
  issn      = {1478-3967},
  doi       = {10.1088/1478-3975/12/1/016013},
  pages     = {14},
  year      = {2015},
  abstract  = {Here we generalize our previous model of molecular motors trafficking subdiffusing cargos in viscoelastic cytosol by (i) including mechano-chemical coupling between cyclic conformational fluctuations of the motor protein driven by the reaction of ATP hydrolysis and its translational motion within the simplest two-state model of hand-over-hand motion of kinesin, and also (ii) by taking into account the anharmonicity of the tether between the motor and the cargo (its maximally possible extension length). It is shown that the major earlier results such as occurrence of normal versus anomalous transport depending on the amplitude of binding potential, cargo size and the motor turnover frequency not only survive in this more realistic model, but the results also look very similar for the correspondingly adjusted parameters. However, this more realistic model displays a substantially larger thermodynamic efficiency due to a bidirectional mechano-chemical coupling. For realistic parameters, the maximal thermodynamic efficiency can transiently be about 50\% as observed for kinesins, and even larger, surprisingly also in a novel strongly anomalous (sub) transport regime, where the motor enzymatic turnovers become also anomalously slow and cannot be characterized by a turnover rate. Here anomalously slow dynamics of the cargo enforces anomalously slow cyclic kinetics of the motor protein.},
  language  = {en}
}