Active dynamics and spatially coherent motion in chromosomes subject to enzymatic force dipoles.

Inspired by recent experiments on chromosomal dynamics, we introduce an exactly solvable model for the interaction between a flexible polymer and a set of motorlike enzymes. The enzymes can bind and unbind to specific sites of the polymer and produce a dipolar force on two neighboring monomers when bound. We study the resulting nonequilibrium dynamics of the polymer and find that the motion of the monomers has several properties that were observed experimentally for chromosomal loci: a subdiffusive mean-square displacement and the appearance of regions of correlated motion. We also determine the velocity autocorrelation of the monomers and find that the underlying stochastic process is not fractional Brownian motion. Finally, we show that the active forces swell the polymer by an amount that becomes constant for large polymers.

The influence of absorbing boundary conditions on the transition path time statistics.

We derive an analytical expression for the transition path time (TPT) distribution for a one-dimensional particle crossing a parabolic barrier. The solution is expressed in terms of the eigenfunctions and eigenvalues of the associated Fokker-Planck equation. The particle exhibits anomalous dynamics generated by a power-law memory kernel, which includes memoryless Markovian dynamics as a limiting case. Our result takes into account absorbing boundary conditions, extending existing results obtained for free boundaries. We show that TPT distributions obtained from numerical simulations are in excellent agreement with analytical results, while the typically employed free boundary conditions lead to a systematic overestimation of the barrier height. These findings may be useful in the analysis of experimental results on the transition path time. A web tool to perform this analysis is freely available.

Motion of kinesin in a viscoelastic medium.

Kinesin is a molecular motor that transports cargo along microtubules. The results of many in vitro experiments on kinesin-1 are described by kinetic models in which one transition corresponds to the forward motion and subsequent binding of the tethered motor head. We argue that in a viscoelastic medium like the cytosol of a cell this step is not Markov and has to be described by a nonexponential waiting time distribution. We introduce a semi-Markov kinetic model for kinesin that takes this effect into account. We calculate, for arbitrary waiting time distributions, the moment generating function of the number of steps made, and determine from this the average velocity and the diffusion constant of the motor. We illustrate our results for the case of a waiting time distribution that is Weibull. We find that for realistic parameter values, viscoelasticity decreases the velocity and the diffusion constant, but increases the randomness (or Fano factor).

The effect of active fluctuations on the dynamics of particles, motors and DNA-hairpins.

Inspired by recent experiments on the dynamics of particles and polymers in artificial cytoskeletons and in cells, we introduce a modified Langevin equation for a particle in an environment that is a viscoelastic medium and that is brought out of equilibrium by the action of active fluctuations caused by molecular motors. We show that within such a model, the motion of a free particle crosses over from superdiffusive to subdiffusive as observed for tracer particles in an in vitro cytoskeleton or in a cell. We investigate the dynamics of a particle confined by a harmonic potential as a simple model for the motion of the tethered head of kinesin-1. We find that the probability that the head is close to its binding site on the microtubule can be enhanced by a factor of two due to active forces. Finally, we study the dynamics of a particle in a double well potential as a model for the dynamics of DNA-hairpins. We show that the active forces effectively lower the potential barrier between the two minima and study the impact of this phenomenon on the zipping/unzipping rate.

Role of trapping and crowding as sources of negative differential mobility.

Increasing the crowding in an environment does not necessarily trigger negative differential mobility of strongly pushed particles. Moreover, the choice of the model, in particular the kind of microscopic jump rates, may be very relevant in determining the mobility. We support these points via simple examples and we therefore address recent claims saying that crowding in an environment is likely to promote negative differential mobility. Trapping of tagged particles enhanced by increasing the force remains the mechanism determining a drift velocity not monotonous in the driving force.

Dynamics of a polymer in an active and viscoelastic bath.

We study the dynamics of an ideal polymer chain in a viscoelastic medium and in the presence of active forces. The motion of the center of mass and of individual monomers is calculated. On time scales that are comparable to the persistence time of the active forces, monomers can move superdiffusively, while on larger time scales subdiffusive behavior occurs. The difference between this subdiffusion and that in the absence of active forces is quantified. We show that the polymer swells in response to active processes and determine how this swelling depends on the viscoelastic properties of the environment. Our results are compared to recent experiments on the motion of chromosomal loci in bacteria.

Transient behaviour of a polymer dragged through a viscoelastic medium.

We study the dynamics of a polymer that is pulled by a constant force through a viscoelastic medium. This is a model for a polymer being pulled through a cell by an external force, or for an active biopolymer moving due to a self-generated force. Using the Rouse model with a memory dependent drag force, we find that the center of mass of the polymer follows a subballistic motion. We determine the time evolution of the length and the shape of the polymer. Through an analysis of the velocity of the monomers, we investigate how the tension propagates through the polymer. We discuss how polymers can be used to probe the properties of a viscoelastic medium.

Exact current statistics of the asymmetric simple exclusion process with open boundaries.

Nonequilibrium systems are often characterized by the transport of some quantity at a macroscopic scale, such as, for instance, a current of particles through a wire. The asymmetric simple exclusion process (ASEP) is a paradigm for nonequilibrium transport that is amenable to exact analytical solution. In the present work, we determine the full statistics of the current in the finite size open ASEP for all values of the parameters. Our exact analytical results are checked against numerical calculations using density matrix renormalization group techniques.

Current fluctuations in the weakly asymmetric exclusion process with open boundaries.

The additivity principle allows a calculation of current fluctuations and associated density profiles in large diffusive systems. We apply this principle to the weakly asymmetric exclusion process (WASEP). We also calculate the cumulant generating function of the current and the density profiles associated with rare currents in finite systems using a numerical approach based on the density matrix renormalization group. Comparison of the two approaches allows us to verify the validity of the additivity principle and to get insight into the finite size scaling theory for current fluctuations in the WASEP. No evidence for a dynamical phase transition is found.

Geometry and topology of knotted ring polymers in an array of obstacles.

We study knotted polymers in equilibrium with an array of obstacles which models confinement in a gel or immersion in a melt. We find a crossover in both the geometrical and the topological behavior of the polymer. When the polymers' radius of gyration, RG, and that of the region containing the knot, R(G,k), are small compared to the distance b between the obstacles, the knot is weakly localized and RG scales as in a good solvent with an amplitude that depends on knot type. In an intermediate regime where R(G)>b>R(G,k), the geometry of the polymer becomes branched. When R(G,k) exceeds b, the knot delocalizes and becomes also branched. In this regime, RG is independent of knot type. We discuss the implications of this behavior for gel electrophoresis experiments on knotted DNA in weak fields.

The size of knots in polymers.

Circular DNA in viruses and bacteria is often knotted. While mathematically problematic, the determination of the knot size is crucial for the study of the physical and biological behaviour of long macromolecules. Here, we review work on the size distribution of these knots under equilibrium conditions. We discuss knot localization in good and poor solvents, or in polymers that are adsorbed on a surface. We also discuss recent evidence that knot size is a crucial quantity in relaxation processes of knotted polymers.

Density-matrix renormalization-group study of current and activity fluctuations near nonequilibrium phase transitions.

Cumulants of a fluctuating current can be obtained from a free-energy-like generating function, which for Markov processes equals the largest eigenvalue of a generalized generator. We determine this eigenvalue with the density-matrix renormalization group for stochastic systems. We calculate the variance of the current in the different phases, and at the phase transitions, of the totally asymmetric exclusion process. Our results can be described in the terms of a scaling ansatz that involves the dynamical exponent z . We also calculate the generating function of the dynamical activity (total number of configuration changes) near the absorbing-state transition of the contact process. Its scaling properties can be expressed in terms of known critical exponents.

Absorbing state phase transitions with quenched disorder.

Quenched disorder--in the sense of the Harris criterion--is generally a relevant perturbation at an absorbing state phase transition point. Here using a strong disorder renormalization group framework and effective numerical methods we study the properties of random fixed points for systems in the directed percolation universality class. For strong enough disorder the critical behavior is found to be controlled by a strong disorder fixed point, which is isomorph with the fixed point of random quantum Ising systems. In this fixed point dynamical correlations are logarithmically slow and the static critical exponents are conjecturedly exact for one-dimensional systems. The renormalization group scenario is confronted with numerical results on the random contact process in one and two dimensions and satisfactory agreement is found. For weaker disorder the numerical results indicate static critical exponents which vary with the strength of disorder, whereas the dynamical correlations are compatible with two possible scenarios. Either they follow a power-law decay with a varying dynamical exponent, like in random quantum systems, or the dynamical correlations are logarithmically slow even for a weak disorder. For models in the parity conserving universality class there is no strong disorder fixed point according to our renormalization group analysis.

Strong disorder fixed point in absorbing-state phase transitions.

The effect of quenched disorder on nonequilibrium phase transitions in the directed percolation universality class is studied by a strong disorder renormalization group approach and by density matrix renormalization group calculations. We show that for sufficiently strong disorder the critical behavior is controlled by a strong disorder fixed point and in one dimension the critical exponents are conjectured to be exact: beta=(3-sqrt[5])/2 and nu( perpendicular )=2. For disorder strengths outside the attractive region of this fixed point, disorder dependent critical exponents are detected. Existing numerical results in two dimensions can be interpreted within a similar scenario.