The time-fractional kinetic equation for the non-equilibrium processes.

In this study, we consider the non-Markovian dynamics of the generic non-equilibrium kinetic process. We summarize the generalized master equation, the continuous and discrete forms of the time-fractional diffusion equation. Using path integral formulation, we generalized the solutions of the Markovian system to the non-Markovian for the non-equilibrium kinetic processes. Then, we obtain the time-fractional kinetic equation for the non-equilibrium systems in terms of free energy. Finally, we introduce a time-fractional equation to analyse time evolution of the open probability for the deformed voltage-gated ion-channel system as an example.

Introduction

A stochastic theory of the kinetics in univariant, non-equilibrium systems that undergoes phase transitions has been formulated in Refs.[1-5]. In these studies, the path integral formulation of the density distribution and conditional probability was given based on the Markovian master equation[6,7] and the kinetic transition probability[8-10].

Although the path integral approximation for the stochastic Markovian and non-Markovian processes are well known and well-studied topics in physics[11-16], to the best of our knowledge, the path integral formulation of the non-Markovian dynamics of the non-equilibrium kinetic processes has never been discussed so far. Therefore, in this study, we generalized the path integral solutions of the Markovian systems to the non-Markovian case for the non-equilibrium kinetic system according to the methods presented in Refs.[1-5].

Therefore, in this study, we consider a generic non-equilibrium kinetic system with non-Markovian dynamics. Following the formalism in Refs.[1-5] we construct the conditional probability for the non-Markovian dynamics by using path integral formulation. Then, we show that the non-Markovian processes in the non-equilibrium system kinetic systems lead to time-fractional kinetic equations. Additionally, we consider the simplest voltage-gated ion-channel system as an example to analyze the time-fractional dynamics of the open probability.

The outline of our paper is given as follow: First, we introduce the non-Markovian master equation. Then, in the following section, we obtain the time-fractional kinetic equation for the non-Markovian dynamics of the non-equilibrium kinetic process. Here, we also introduce the fractional kinetic equation of the deformed voltage-gated ion-channel system. Finally, we review the obtained results and present a brief discussion.

Generalized Master equation

First, we briefly give the information about the master equation for the non-Markovian dynamics which is called as generalized master equation (GME) in the literature[17-19]. The time evolution of finding probability P(x, t) can be presented on the one-dimensional space bywhere is the transfer probability kernel from the position x to . As seen from Eq. (1) that GME is very different from Markovian form. Here, we use the probability function of the generalized master equation kernel by reformulated the path integral formulation based on previous approximations[1-4].

The Fourier-Laplace transformation of Eq. (1) is written aswhere u is the Laplace variable, k is the wave number and denotes a Fourier convolution of the f and g functions. Dividing by u, after Laplace inversion and differentiation we obtain another representationof the Eq. (1). The new kernel of the master equation is given by , i.e. where can be presented with . In this representation, the kernels and are responsible for spatial correlations and memory in any stochastic process. These kernels are classified by the finite characteristic waiting time T and the finite jump length variance . It is known that, in the non-Markovian process, while the jump length variance is finite, the waiting time T diverges due to the spatial deformations, entropic restrictions or other memory effects.

For a non-Markovian processes is represented bywhere is Gamma function, is the macroscopic relaxation parameter and the exponent takes the value between . We note that Eq. (4) corresponds to the long-tailed waiting time probability distribution.

In the new situation, the new kernel is given by

The solution of the GME shows a strong dependence on its stochastic history. Therefore, the resulting equation is

Equation (6) includes the defining expression[20]where is the Riemann–Liouville fractional derivative[20]. Time-fractional master equation can be expressed in te form

Here we briefly summarize that, the non-Markovian dynamics in the continuum limit leads to the time-fractional differential equation Eq. (8). The discrete form of Eq. (8) is given by

Using Eq. (9) we will construct the path integral formulation of the conditional probability function P(x, t).

Time-fractional kinetic equation

It should be noted that the kernel of the spatial correlation for the kinetic processes can be defined in terms of Helmholtz free energywhere is the inverse temperature, is a constant and F is the free energy of the thermodynamic system. This definition of the kinetic transition probability in Eq. (10) suggested by Langer[8-10], which is an extension of a model proposed by Glauber[21] based on Zwanzig theory[22,23].

Afterwards, following previous theoretical schema[1-5] we construct the path integral definition of the probability function for the non-Markovian dynamics of the generic non-equilibrium kinetic. Thus, we write the master equation in Eq. (9) can be given in the formwhere is written aswhich causes the Kramers-Moyal expansion[24,25]where operates at the same time both and . The sums are over all possible values of the multi-indices m. It should also be noted that we set and in Eqs. (12) and  (13) for the convenience.

At this point, by using definition of the derivatives we can arrange the left side of Eq. (11)[1-5], and then we can write the probability function as

We define the Fourier transform of Eq. (14) as

This of course is represented bywhere is obtained from Eq. (13) by replacing with and x with [2-4]. On the other hand, the inverse Fourier transform of Eq. (16) is given by

Here, introducing Eqs. (16) into (17) and recognizing that for small the curly bracket in Eq. (16) is an exponential, in this case, we obtain

The kernel of Eq. (18) can be defined aswhere . This kernel represents the path integral formulation of the conditional probability for non-Markovian kinetics. We clearly see that the integral argument in Eq. (19) corresponds to Lagrangian of the system, which is given aswhere can be read as Hamiltonian. The path integral in Eq. (19) can be defined as the limit of the multiple integralwhen . We consider here that the transition between the small paths along the trajectory are independent of each other. Now, by using Eq. (10) we can write asor very small values, the integrations can be obtained as

Introducing Eqs. (23) into (21) and integrating over k, we get the time-fractional kernel

The path of extreme probability are those for which is extremized, namely, time-fractional action integral satisfy the conditionwhere is the time-fractional Lagrangian of the system. Euler-Lagrange equation in the action Eq. (25) can be solved due to k as

The solution of Eq. (26) gives the time-fractional kinetic equationwhere . As can be seen from Eq. (27) that the time evolution of the stochastic variable x is governed by integro-differential operator. The operator points out that the solution of the kinetic equation Eq. (27) possesses Mittag–Leffler function.

An example: the deformed ion-channel systems

In the simplest voltage-gated ion channel system, it is assumed that the channels are located on a two-dimensional membrane. At the equilibrium, the channels are open or closed depending on temperature and membrane voltage. Besides, the status of the channels depends on their intrinsic properties, the state of the channel is mainly determined by the membrane potential at a constant temperature. In these models, for simplicity, it is assumed that the channels are identical and distributed randomly on the membrane surface.

In this simple schema, the number of the channels is given by where denotes channels with energy and admits the open channels with energy , respectively. In the statistical framework, the Helmholtz energy of such channel model around equilibrium can be written as[26,27]where is the open probability under external potential V, corresponds to the number of configurations which is given by and z is the number of charges, and is the charge of electron. It is shown that the open probability can be found from first derivative of the Helmholtz free energy in Eq. (28) as to the variable x. Thus it can be written aswhere is the critical threshold voltage value. In the case of Markovian dynamics, it is assumed that the half of the channels on the membrane surface are open.

The open probability of the channel system in Eq. (29) is a well-known and well-studied topic in the literature. Furthermore, the kinetic behavior of the voltage gated ion channels was also examined within the framework of the Markovian formalism[26,27]. However, in the presence of the deformations such as a genetic mutation in the channel, we expect that ion channels have non-Markovian dynamics due to decoupling dynamics as stated in the previous section. As a result, channel conductivity is damaged and the cell may not be able to perform its previous tasks. If the channel kinetics is analyzed within the framework of the non-Markovian formalism, the time-fractional kinetic for the open probability of the ion channel system can be obtained as

As seen from this particular example, one can apply the fractional numerical integration method to Eq. (30) to obtain the time evolution of the open probability for the non-Markovian ion channel system far from steady-state. Hence, the time variation of the open probability and relaxation parameters of the non-Markovian ion channels can be obtained by using Eq. (30). The fractional derivative in Eq. (30) clearly indicates that the solution of the relaxation drastically deviates from the Markovian solution given in Refs.[26,27].

Conclusion

In this study, firstly, we briefly introduce the generalized master equation for non-Markovian dynamics. We also present the continuous and discrete forms of the time-fractional diffusion equation in the same section. In the subsequent section, we generalized the path integral solutions of the Markovian system to the non-Markovian for the non-equilibrium kinetic processes. Using path integral formulation we obtain the time-fractional kinetic equation for the non-equilibrium systems in terms of free energy. Here, we consider, as an example, a voltage-gated ion-channel system that behaves as a non-equilibrium system under cell voltage. We introduce the time-fractional kinetic equation for the open probability of the simple ion-channel system.

The non-Markovian dynamics and non-equilibrium behavior of the physical systems are very important two topics in physics. Indeed, many physical systems in nature may have one or both of these properties. Systems can be considered as systems which far from equilibrium that cannot reach equilibrium in very large time scales, and such systems can also be considered as systems that reach equilibrium in short time scales. On the other hand, we know that complex or disordered physical systems generally have non-Markovian dynamics. Therefore, it is very interesting to examine the time evolutions of kinetic systems with both non-equilibrium and Markovian dynamics.

The method presented in the study can be applied to the other physical systems to analyze the time-dependent evolution of the kinetic systems such as the interface dynamics of the phase transitions, kinetic flows, bacterial growth phenomena, other deformed ion channels, internet networks, and chemical kinetics[28-33].