key: cord-0011820-nle4oisc
authors: Berezhkovskii, Alexander M.; Makarov, Dmitrii E.
title: Detailed balance for diffusion in a potential with trapping and forward–backward symmetry of trapping time distributions
date: 2020-06-14
journal: J Chem Phys
DOI: 10.1063/1.5142566
sha: 89f43b5e9d08ab3cedc590025b658a9ce25678c9
doc_id: 11820
cord_uid: nle4oisc

For particles diffusing in a potential, detailed balance guarantees the absence of net fluxes at equilibrium. Here, we show that the conventional detailed balance condition is a special case of a more general relation that works when the diffusion occurs in the presence of a distributed sink that eventually traps the particle. We use this relation to study the lifetime distribution of particles that start and are trapped at specified initial and final points. It turns out that when the sink strength at the initial point is nonzero, the initial and final points are interchangeable, i.e., the distribution is independent of which of the two points is initial and which is final. In other words, this conditional trapping time distribution possesses forward–backward symmetry.

Consider a particle diffusing along a one-dimensional coordinate x in a constraining potential U(x), U(x)| |x|→∞ → ∞. Its propagator (Green's function), denoted by G(x, t|x 0 ), is the probability density of finding the particle at point x at time t, conditional on that the particle started from x = x 0 at t = 0. The propagator satisfies the Smoluchowski equation,

subject to the initial condition G(x, 0|x 0 ) = δ(x − x 0 ). In Eq. (1), D is the particle diffusivity, and β = 1/(kBT), where kB and T are the Boltzmann constant and absolute temperature. This propagator describes the relaxation of the initial δ-distribution to the equilibrium one, 

Once G(x, t|x 0 ) is known, the time evolution of any initial distribution p(x 0 , 0) can be determined by computing the convolution of the two functions. The propagator satisfies the detailed balance condition,

which is a relationship between the propagators G(x, t|x 0 ) and G(x 0 , t|x) that guarantees that there are no net fluxes between arbitrary chosen points x and x 0 at equilibrium. Ultimately, detailed balance is a consequence of time reversibility of diffusion trajectories contributing to the propagator. [1] [2] [3] [4] When particle diffusion occurs in the presence of an arbitrary sink γ(x), γ(x) ≥ 0, Eq. (1) takes the form:

The Journal of Chemical Physics NOTE scitation.org/journal/jcp

Here, the propagator vanishes as t → ∞, since the particle is trapped by the sink,

Although there is no equilibrium in this case, the propagator nevertheless satisfies the detailed balance condition [Eq. (3)]. This is one of the main results of this work. We use this result to show that the trapping time distributions possess forward-backward symmetry. To explain what this symmetry is, consider an ensemble of particles starting from x = xi, where the sink strength is not equal to zero, γ(xi) > 0. Some of the particles are trapped by the sink γ(x) at a point x = x f . We use these ensemble members to find the distribution of their lifetime, denoted by px i →x f (t). We can repeat this procedure for an ensemble of particles starting from x = x f and find the lifetime distribution of the ensemble members trapped at the point x = xi, denoted by px f →x i (t). It turns out that these trapping time distributions are identical,

and hence possess forward-backward symmetry. This is another main result of this work. To avoid confusion, it is worth noting that the trapping time distributions differ from the conditional first passage time distributions discussed, e.g., in Ref. 5 . The equality in Eq. (6) is applicable for transitions between points, in which γ(x) > 0. As an example, consider channelfacilitated transport of diffusing particles between two compartments separated by a membrane. To a good approximation, particle dynamics in the channel can be described as one-dimensional diffusion along the channel axis in a potential of mean force. Particle escape from the channel is modeled by imposing radiation boundary conditions at the channel ends. 6,7 A particle entering the channel either returns to the same compartment from which it entered or passes through the channel and escapes on the opposite side of the membrane. Formally, one can replace the radiation boundaries at the channel ends by reflecting ones and place localized (δ-function) sinks near them with the sink strengths equal to the trapping rates entering the radiation boundary conditions. The particle is trapped by the sink located at the channel entrance in the first of the two above-mentioned scenarios and by the sink located at the channel exit in the second scenario. For particles passing through the channel, it has been shown that the distribution of their residence time in the channel is independent of the passage direction 8-10 and hence satisfies Eq. (6).

The passage time distribution remains finite even when the channel ends are perfectly absorbing boundaries and the passage probabilities vanish. In this limiting case, they describe distributions of the direct-transit time also referred to as the transition path time. 11 This time, which is the duration of a trajectory fragment during which the particle diffuses from its starting point to the end point without touching the starting point, has attracted a flurry of recent theoretical 12, 13 and experimental 14, 15 efforts. It is known that the direct-transit time distributions are also independent of the passage direction. 11, [16] [17] [18] [19] Thus, the identity of the trapping time distributions [Eq. (6)] is a generalization of known results and contains these results as special cases. The identities of the forward-backward transit time distributions discussed in Refs. 8-11 and 16-19 were proved assuming that the particle can be trapped only at the interval end points. Here, we show that this assumption can be relaxed, and the identities remain unchanged even if the particle can be trapped within the interval. As examples, we mention transport of molecular motors along microtubules and protein search for the target site on DNA, [20] [21] [22] where the detachment from the microtubule or a DNA may be modeled by introducing a sink.

Our interest in the trapping time distributions and related problems was stimulated by recent single-molecule and singlechannel experiments, in which such distributions can potentially be measured with high accuracy, 14,15,23-26 thereby providing new insights into the mechanisms of channel transport, biomolecular folding and binding, and other biophysical phenomena.

We now prove that the propagator G(x, t|x 0 ) of a particle diffusing in the presence of an arbitrary sink satisfies the detailed balance condition [Eq. (3)]. To this end, we write this propagator as

Substituting this into Eq. (4), one can check that g(x, t|x 0 ) satisfies a Schrödinger-like equation,

subject to the initial condition g(x, 0|x 0 ) = δ(x − x 0 ). Thus, g(x, t|x 0 ) is Green's function for Eq. (8), in which the Hamiltonian H is a Hermitian operator. The eigenfunction expansion (spectral representation) of g(x, t|x 0 ) is given by

where εn (εn < 0) and φn(x) are the eigenvalues and eigenfunctions ofĤ,Ĥ φn(x) = εnφn(x), n = 1, 2, . . . .

Substituting g(x, t|x 0 ) in Eq. (9) into Eq. (7), we obtain

Multiplying this by e −βU(x 0 ) , we arrive at

which is the detailed balance condition for nonequilibrium systems [Eq. (3)].

Note that, although we focus on the case of a constraining potential ( U(x)| |x|→∞ → ∞), the results in Eqs. (3) and (6) remain unchanged if we relax this assumption and allow the particle to escape to infinity. In this case, the eigenvalue spectrum of the Hamil-tonianĤ defined in Eq. (8) is continuous, and the summation in Eqs. (9) , (11) , and (12) should be replaced by integration. The rest of the above reasoning remains unchanged.

We illustrate our result with two examples. In the first, the sink strength is distributed uniformly along x, γ(x) = γ. Such a model may, for example, be used to describe the detachment of a molecular motor moving along its track, if it occurs with the same probability regardless of the motor position. Let G 0 (x, t|x 0 ) be Green's function in the absence of trapping, which satisfies Eq. (1). It is easy to verify that the solution of Eq. (4) in the case where the trapping rate is coordinateindependent is simply given by G(x, t|x 0 ) = e −γt G 0 (x, t|x 0 ), and since G 0 (x, t|x 0 ) satisfies detailed balance [Eq. (3)] so does Green's function G(x, t|x 0 ).

The second example represents the opposite extreme, with spatially localized trapping. Specifically, we consider the following coordinate dependence of the sink strength,

is Dirac's delta-function. Such a scenario arises, for example, when the dynamics of a particle inside a membrane channel (a < x < b) is considered, and sinks are introduced to describe the escape from the channel to the bulk on its either side. 6, 7 Green's function G(x, t|x 0 ) of such a system can be related to Green's function G 0 (x, t|x 0 ) in the absence of trapping using the Dyson equation-see, e.g., Ref. 27 , 

Applying this result to x = a,b and x 0 = a,b, it can be shown that

and a similar expression is obtained forĜ(a, s|b) by simply exchanging a and b. Unlike the case of the coordinate-independent sink strength, the temporal behavior of Green's functions in the case of two localized sinks is nontrivial and depends on both γa and γ b , as well as on the sink locations. Yet when taking a ratiô G(a, s|b)/Ĝ(b, s|a), we observe that it is identical to the ratiô G 0 (a, s|b)/Ĝ 0 (b, s|a) = exp{−β[U(a) − U(b)]}, which proves the detailed balance condition in the presence of trapping for this model.

To prove the forward-backward symmetry of the trapping time distributions, consider a particle starting from x = xi, where γ(xi) ≠ 0, at time t = 0. Let qx i →x f (t) be the flux trapped by the sink γ(x) at point x = x f at time t,

The integral of this flux over time from zero to infinity, denoted by Φx i →x f ,

is the spatial distribution of the location of the trapping point normalized to unity,

The time distribution for xi → x f trapping, denoted by px i →x f (t), is the flux qx i →x f (t) normalized to its integral over time Φx i →x f ,

Correspondingly, the time distribution for x f → xi trapping, denoted by px f →x i (t), is given by

Using the detailed balance condition [Eq. (3)], one can see that the trapping time distributions in Eqs. (19) and (20) are equal, as claimed in Eq. (6). In summary, detailed balance for diffusing particles [Eq. (3)] is conventionally discussed for systems with particle conservation. Here, it is shown that Eq. (3) is true even if the particles can be trapped by a sink γ(x), and hence, the conventional detailed balance condition is a special case where γ(x) = 0. We used the detailed balance condition [Eq. 6)]. Finally, we mention that the present study complements our recent work on forward-backward symmetry of transition path time distributions in nonequilibrium systems, which, in contrast to this work, is mainly focused on Markovian dynamics on networks. 28 

Stochastic Processes in Physics and Chemistry

The Fokker-Planck Equation

Handbook of Stochastic Methods for Physics

Statistical Mechanics

A Guide to First-Passage Processes

Single Molecule Science: Physical Principles and Models

Motor Proteins and Molecular Motors

Stochastic Processes in Cell Biology

Theoretical perspectives on biological machines

Chemical Dynamics in Condensed Phases