Energy transfer in a molecular motor in the Kramers regime

Nonequilibrium master equation for interacting Brownian particles in a deep-well periodic potential

Employing a creation and annihilation operator formulation, we derive an approximate many-body master equation describing discrete hopping from the more general continuous description of Brownian motion on a deep-well nonequilibrium periodic potential. The many-body master equation describes interactions of arbitrary strength and range arising from a "top-hat" two-body interaction potential. We show that this master equation reduces to the well-known asymmetric simple exclusion process and the zero range process in certain regimes. We also use the creation and annihilation operator formalism to derive results for the steady-state drift and the number fluctuations in special cases, including the unexplored limit of weak interparticle interactions.

Enhanced diffusion and the eigenvalue band structure of Brownian motion in tilted periodic potentials

We consider enhanced diffusion for Brownian motion on a tilted periodic potential. Expressing the effective diffusion in terms of the eigenvalue band structure, we establish a connection between band gaps in the eigenspectrum and enhanced diffusion. We explain this connection for a simple cosine potential with a linear force and then generalize to more complicated potentials including one-dimensional potentials with multiple frequency components and nonseparable multidimensional potentials. We find that potentials with multiple band gaps in the eigenspectrum can lead to multiple maxima or broadening of the force-diffusion curve. These features are likely to be observable in experiments.

Thermodynamic uncertainty relations and molecular-scale energy conversion

The thermodynamic uncertainty relation (TUR) is a universal constraint for nonequilibrium steady states that requires the entropy production rate to be greater than the relative magnitude of current fluctuations. It has potentially important implications for the thermodynamic efficiency of molecular-scale energy conversion in both biological and artificial systems. An alternative multidimensional thermodynamic uncertainty relation (MTUR) has also been proposed. In this paper we apply the TUR and the MTUR to a description of molecular-scale energy conversion that explicitly contains the degrees of freedom exchanging energy via a time-independent multidimensional periodic potential. The TUR and the MTUR are found to be universal lower bounds on the entropy generation rate and provide upper bounds on the thermodynamic efficiency. The TUR is found to provide only a weak bound while the MTUR provides a much tighter constraint by taking into account correlations between degrees of freedom. The MTUR is found to provide a tight bound in the near or far from equilibrium regimes but not in the intermediate force regime. Collectively, these results demonstrate that the MTUR is more appropriate than the TUR for energy conversion processes, but that both diverge from the actual entropy generation in certain regimes.

Reconstructing free-energy landscapes for cyclic molecular motors using full multidimensional or partial one-dimensional dynamic information

Diffusion on a free-energy landscape is a fundamental framework for describing molecular motors. In the landscape framework, energy conversion between different forms of energy, e.g., chemical and mechanical, is explicitly described using multidimensional nonseparable potential landscapes. We present a k-space method for reconstructing multidimensional free-energy landscapes from stochastic single-molecule trajectories. For a variety of two-dimensional model potential landscapes, we demonstrate the robustness of the method by reconstructing the landscapes using full dynamic information, i.e., simulated two-dimensional stochastic trajectories. We then consider the case where the stochastic trajectory is known only along one dimension. With this partial dynamic information, the reconstruction of the full two-dimensional landscape is severely limited in the majority of cases. However, we reconstruct effective one-dimensional landscapes for the two-dimensional model potentials. We discuss the interpretation of the one-dimensional landscapes and identify signatures of energy conversion. Finally, we consider the implications of these results for biological molecular motors experiments.

Analysing single-molecule trajectories to reconstruct free-energy landscapes of cyclic motor proteins

Stochastic trajectories measured in single-molecule experiments have provided key insights into the microscopic behaviour of cyclic motor proteins. However, the fundamental free-energy landscapes of motor proteins are currently only able to be determined by computationally intensive numerical methods that do not take advantage of available single-trajectory data. In this paper we present a robust method for analysing single-molecule trajectories of cyclic motor proteins to reconstruct their free-energy landscapes. We use simulated trajectories on model potential landscapes to show the reliable reconstruction of the potentials. We determine the accuracy of the reconstruction method for common precision limitations and show that the method converges logarithmically. These results are then used to determine the experimental precision required to reconstruct a potential with a desired accuracy. The key advantages of the method are that it is simple to implement, is free of numerical difficulties that plague existing methods and is easily generalizable to higher dimensions.

Tight-binding derivation of a discrete-continuous description of mechanochemical coupling in a molecular motor

We present a tight-binding derivation of a discrete-continuous description of mechanochemical coupling in a molecular motor. Our derivation is based on the continuous diffusion equation for overdamped Brownian motion on a time-independent tilted periodic potential in two dimensions. The free-energy potential is nonseparable to allow coupling between the chemical and mechanical degrees of freedom. We formally discretize the chemical coordinate by expanding in Wannier states that are localized along the chemical coordinate and parametrized along the mechanical coordinate. A discrete-continuous equation is derived that is valid for anisotropic systems with weak mechanochemical coupling and deep potential wells along the chemical coordinate. The discrete-continuous description is consistent with established theoretical models of molecular motors with discrete chemical states but is constrained by the underlying continuous two-dimensional potential. In particular, we derive analytic expressions for the effective potential along the mechanical coordinate and for the rate of thermal hopping between chemical states. We determine the thermodynamic efficiency of energy conversion and find that, for a molecular motor with one chemical state per cycle, the derived discrete-continuous equation can accurately describe mechanochemical coupling but cannot describe energy conversion.

Reconstructing free-energy landscapes for nonequilibrium periodic potentials

We present a method for reconstructing the free-energy landscape of overdamped Brownian motion on a tilted periodic potential. Our approach exploits the periodicity of the system by using the k-space form of the Smoluchowski equation and we employ an iterative approach to determine the nonequilibrium tilt. We reconstruct landscapes for a number of example potentials to show the applicability of the method to both deep and shallow wells and near-to- and far-from-equilibrium regimes. The method converges logarithmically with the number of Fourier terms in the potential.

Self-induced temperature gradients in Brownian dynamics

Brownian systems often surmount energy barriers by absorbing and emitting heat to and from their local environment. Usually, the temperature gradients created by this heat exchange are assumed to dissipate instantaneously. Here we relax this assumption to consider the case where Brownian dynamics on a time-independent potential can lead to self-induced temperature gradients. In the same way that externally imposed temperature gradients can cause directed motion, these self-induced gradients affect the dynamics of the Brownian system. The result is a coupling between the local environment and the Brownian subsystem. We explore the resulting dynamics and thermodynamics of these coupled systems and develop a robust method for numerical simulation. In particular, by focusing on one-dimensional situations, we show that self-induced temperature gradients reduce barrier-crossing rates. We also consider a heat engine and a heat pump based on temperature gradients induced by a Brownian system in a nonequilibrium potential.

Numerical study of the tight-binding approach to overdamped Brownian motion on a tilted periodic potential

We present a numerical study of the tight-binding approach to overdamped Brownian motion on a tilted periodic potential. In the tight-binding method the probability density is expanded on a basis of Wannier states to transform the Smoluchowski equation to a discrete master equation that can be interpreted in terms of thermal hopping between potential minima. We calculate the Wannier states and hopping rates for a variety of potentials, including tilted cosine and ratchet potentials. For deep potential minima the Wannier states are well localized and the hopping rates between nearest-neighbor states are qualitatively well described by Kramers' escape rate. The next-nearest-neighbor hopping rates are negative and must be negligible compared to the nearest-neighbor rates for the discrete master equation treatment to be valid. We find that the validity of the master equation extends beyond the quantitative applicability of Kramers' escape rate.

Local discretization method for overdamped Brownian motion on a potential with multiple deep wells

We present a general method for transforming the continuous diffusion equation describing overdamped Brownian motion on a time-independent potential with multiple deep wells to a discrete master equation. The method is based on an expansion in localized basis states of local metastable potentials that match the full potential in the region of each potential well. Unlike previous basis methods for discretizing Brownian motion on a potential, this approach is valid for periodic potentials with varying multiple deep wells per period and can also be applied to nonperiodic systems. We apply the method to a range of potentials and find that potential wells that are deep compared to five times the thermal energy can be associated with a discrete localized state while shallower wells are better incorporated into the local metastable potentials of neighboring deep potential wells.

Intrinsic irreversibility limits the efficiency of multidimensional molecular motors

We consider the efficiency limits of Brownian motors able to extract work from the temperature difference between reservoirs or from external thermodynamic forces. These systems can operate in a variety of modes, including as isothermal engines, heat engines, refrigerators, and heat pumps. We derive analytical results showing that certain classes of multidimensional Brownian motor, including the Smoluchowski-Feynman ratchet, are unable to attain perfect efficiency (Carnot efficiency for heat engines). This demonstrates the presence of intrinsic irreversibilities in their operating mechanism. We present numerical simulations showing that in some cases the loss process that limits efficiency is associated with vortices in the probability current.

Tight-binding approach to overdamped Brownian motion on a bichromatic periodic potential

We present a theoretical treatment of overdamped Brownian motion on a time-independent bichromatic periodic potential with spatially fast- and slow-changing components. In our approach, we generalize the Wannier basis commonly used in the analysis of periodic systems to define a basis of S states that are localized at local minima of the potential. We demonstrate that the S states are orthonormal and complete on the length scale of the periodicity of the fast-changing potential, and we use the S-state basis to transform the continuous Smoluchowski equation for the system to a discrete master equation describing hopping between local minima. We identify the parameter regime where the master equation description is valid and show that the interwell hopping rates are well approximated by Kramers' escape rate in the limit of deep potential minima. Finally, we use the master equation to explore the system dynamics and determine the drift and diffusion for the system.

Landauer's blowtorch effect as a thermodynamic cross process: Brownian cooling

The local heating of a selected region in a double-well potential alters the relative stability of the two wells and gives rise to an enhancement of population transfer to the cold well. We show that this Landauer's blowtorch effect may be considered in the spirit of a thermodynamic cross process linearly connecting the flux of particles and the thermodynamic force associated with the temperature difference and consequently ensuring the existence of a reverse cross effect. This reverse effect is realized by directing the thermalized particles in a double-well potential by application of an external bias from one well to the other, which suffers cooling.

Thermal fluctuation statistics in a molecular motor described by a multidimensional master equation

We present a theoretical investigation of thermal fluctuation statistics in a molecular motor. Energy transfer in the motor is described using a multidimensional discrete master equation with nearest-neighbor hopping. In this theory, energy transfer leads to statistical correlations between thermal fluctuations in different degrees of freedom. For long times, the energy transfer is a multivariate diffusion process with constant drift and diffusion. The fluctuations and drift align in the strong-coupling limit enabling a one-dimensional description along the coupled coordinate. We derive formal expressions for the probability distribution and simulate single trajectories of the system in the near- and far-from-equilibrium limits both for strong and weak coupling. Our results show that the hopping statistics provide an opportunity to distinguish different operating regimes.