Entropy Production Rate of a One-Dimensional Alpha-Fractional Diffusion Process

In this paper, the one-dimensional α-fractional diffusion equation is revisited. This equation is a particular case of the timeand space-fractional diffusion equation with the quotient of the orders of the timeand space-fractional derivatives equal to one-half. First, some integral representations of its fundamental solution including the Mellin-Barnes integral representation are derived. Then a series representation and asymptotics of the fundamental solution are discussed. The fundamental solution is interpreted as a probability density function and its entropy in the Shannon sense is calculated. The entropy production rate of the stochastic process governed by the α-fractional diffusion equation is shown to be equal to one of the conventional diffusion equation.


Introduction
In the literature, many different kinds of the time-, space, and time-and space-fractional diffusion equations have been already introduced and analyzed.The fractional derivatives contained in the fractional diffusion equations are defined in the Riemann-Liouville, Weyl, Caputo, Riez, or Riesz-Feller sense to mention only some of the most used types of the fractional derivatives.From the mathematical viewpoint, all of these fractional diffusion equations can be seen as generalizations of the conventional diffusion equation and thus are worth to be investigated.On the other hand, it is not clear at all, what kinds of the fractional diffusion equations could/should be employed as mathematical models, say, for describing the phenomena of the anomalous diffusion (see, e.g., the recent survey paper [1] for about three hundred references to the relevant works).Usually, the anomalous diffusion processes are defined as those that do not longer follow the Gaussian statistics on the long time intervals.Especially, the linear time dependence of the mean squared displacement of the diffusing particles does not hold any more and has to be either replaced with a different (mainly power-law) dependence or the mean squared displacement does not exist at all.As a rule, the stochastic processes governed by different kinds of the fractional diffusion equations are non Gaussian and thus can be seen as potential candidates for the role of mathematical models for anomalous diffusion processes.
One more important characteristic of the diffusion processes is their entropy and the entropy production rate.The concept of entropy was first introduced in the macroscopic thermodynamics and then extended for description of some phenomena in statistical mechanics, information theory, ergodic theory of dynamical systems, etc. Historically, many definitions of entropy were proposed and applied in different knowledge areas.In this paper, we employ the statistical concept of entropy that goes back to Shannon and was introduced by him in the theory of communication and transmission of information (see [2]).The entropy of the processes governed by the time-and space-fractional diffusion equations has been discussed in [3][4][5][6], respectively.It is worth mentioning that according to [3,5] the entropy production rates for the time-and the space-fractional diffusion equations depend on the derivative order α of the time-or space-fractional derivative, respectively, and increase with increasing of α from 1 (diffusion) to 2 (wave propagation) that results in the so called entropy production paradox.In [7], entropy behavior of solutions to the one-dimensional neutral-fractional equation that contains fractional derivatives of the same order α, 1 ≤ α ≤ 2 both in space and in time, has been considered.It was shown in [7] that the entropy production rate of solutions to the neutral-fractional equation does not depend on the equation order α and is twice as much as the entropy production rate of solutions to the conventional diffusion equation.In this paper, we show that the entropy production rate of the fundamental solution to the α-fractional diffusion equation is exactly the same as in the case of the conventional diffusion equation.Thus the α-fractional diffusion equation combines the properties of the anomalous diffusion (the mean squared displacement of the diffusing particles does not exists) and of the conventional diffusion (the same entropy production rate) and could be considered to be a kind of a "natural fractionalization" of the diffusion equation.In this paper we restrict ourselves to the classical Shannon entropy, other kinds of generalized entropies (see e.g., [8][9][10] and the references therein) will be considered elsewhere.
From the mathematical viewpoint, the α-fractional diffusion equation is a particular case of the one-dimensional space-time fractional diffusion equation that has been considered in [11] in detail.The equation studied in [11] contains the Riesz-Feller derivative of order α ∈ (0, 2] and skewness θ and the Caputo fractional derivative of order β ∈ (0, 2].In particular, it was shown in [11] that the fundamental solution to the space-time fractional diffusion equation can be interpreted as a spatial probability density function evolving in time if The α-fractional diffusion equation we deal with in this paper corresponds to the case α = 2β, 0 ≤ β ≤ 1, θ = 0 in the space-time fractional diffusion equation considered in [11] and possesses some remarkable properties that do not hold true for solutions of the general equation.
The rest of the paper is organized as follows.In the 2nd section, the basic definitions, problem formulation, and some analytical results for the initial-value problems for the one-dimensional α-fractional equation are presented.Among other things, the Mellin-Barnes integral representation of the fundamental solution as well as its series representation and asymptotics are given.The last section is devoted to a probabilistic interpretation of the fundamental solution to the one-dimensional α-fractional diffusion equation.In particular, the Shannon entropy and the entropy production rate are calculated.The entropy production rate of the stochastic process governed by the α-fractional diffusion equation is shown to be independent on the equation order α and is exactly the same as the entropy production rate of the conventional diffusion process.

Problem Formulation
In this paper, we deal with the one-dimensional α-fractional diffusion equation in the form In the Equation (1), D α t is the Caputo time-fractional derivative of order α defined by I α , α ≥ 0 being the Riemann-Liouville fractional integral and Γ the Euler gamma function.For α = n, n ∈ IN, the Caputo fractional derivative coincides by definition with the derivative of order n.
For a sufficiently well-behaved function f , the Riesz fractional derivative −(−∆) α is defined as a pseudo-differential operator with the symbol −|κ| 2α (see e.g., [11][12][13][14]): F being the Fourier transform of a function f defined by the formula The Riesz fractional derivative (3) can be represented as a hypersingular integral under the condition 0 < α < 1 (see [14] for the case α = 1 2 and [12] for the general case) For α = 1 2 , the relation ( 4) can be interpreted in terms of the Hilbert transform where the integral is understood in the sense of the Cauchy principal value as first noted in [15] and then revisited and stated more precisely in [12,16].
Let us note that the Riesz fractional derivative is a symmetric operator with respect to the space variable x.Because of the relation −|κ| 2α = −(κ 2 ) α it can be formally interpreted as i.e., as a power of the self-adjoint and positive definite operator − d 2 dx 2 .For α = 1, Equation ( 1) is reduced to the one-dimensional diffusion equation.In what follows, we focus on the case 0 ≤ α < 1 because the case α = 1 (diffusion equation) is well studied in the literature.
In the rest of the paper, we consider the initial-value problem for the Equation (1).In doing so, we are mostly interested in behavior and properties of the fundamental solution (Green function) G α = G α (x, t) of the Equation (1), i.e., in its solution with the initial condition ϕ(x) = δ(x), δ being the Dirac delta function.

Fundamental Solution of the Alpha-Fractional Diffusion Equation
In this subsection, we follow the derivations presented in [11] for the more general case of the one-dimensional time-space fractional diffusion equation with some minor modifications.To determine the fundamental solution G α let us apply the Fourier transform to the Equation (1) and to the initial conditions (5) with ϕ(x) = δ(x).Using definition of the Riesz fractional derivative, for the Fourier transform Ĝα we get then the initial-value problem Ĝ(κ, 0) = 1 (6) for the fractional differential equation The unique solution of ( 6), ( 7) is given by the expression (see e.g., [17,18]) in terms of the Mittag-Leffler function E α that is defined as an convergent power series As follows from the well-known asymptotic formula the Fourier transform Ĝα belongs to L 1 (I R) with respect to κ under the condition α > 1 2 .In the further discussions we suppose that this conditions holds true.Then we can apply the inverse Fourier transform and get the representation that can be rewritten as the cos-Fourier transform: Now we are going to apply the technique of the Mellin integral transform to deduce the Mellin-Barnes representation of the fundamental solution.For the reader's convenience, some basic elements of the Mellin integral transform theory are presented below.
The Mellin integral transform of a function f is defined by the formula and the inverse Mellin integral transform by the formula The Mellin integral transform and its inversion exist under the following sufficient conditions (see e.g., [19] then the Mellin transform (12) of the function f exists and is analytic in the vertical strip If f is piecewise differentiable and f (r) r γ−1 ∈ L c (0, ∞) , then the Formula (13) holds true in all points where f is continuous.The integral in (13) must be understood in the sense of the Cauchy principal value.

Let us denote by
M ↔ the juxtaposition of a function f with its Mellin transform f * .With this notation the convolution theorem for the Mellin convolution reads as follows: Now let us return back to the integral representation (11) of the fundamental solution G α and consider the cases a) x = 0 and b) x = 0.
(a) For x = 0, the integral at the right-hand side of ( 11) can be interpreted as the Mellin integral transform of the Mittag-Leffler function at the point s = 1 2α .It converges under the condition α > 1 2 and its value is given by the formula (see e.g., [19]) that is in accordance with the above formula for the value of G α (0, t).
(b) In the case x = 0, we recognize that the integral at the right-hand side of ( 11) can be interpreted as the Mellin convolution of the functions κ at the point y = 1 |x| .Using the known Mellin integral transforms of the cos-function and the Mittag-Leffler function as well as some elementary properties of the Mellin integral transform (see [19,20]) we get the formulas: These formulas together with the convolution theorem and the inverse Mellin integral transform lead to the following Mellin-Barnes representation of the fundamental solution G α : The linear variables substitution s → 2s in the integral at the right-hand side of (15) leads to the representation and then to the formula by using the reflection formulas for the gamma function.
In the next section, we need one more Mellin-Barnes integral representation of the fundamental solution G α that is obtained from the Formula (16) using the linear variables substitution s → −s: −s ds, − min{1/2, α} < γ < 0.
(18) From Equation ( 18), a useful representation of the fundamental solution G α in terms of an auxiliary function L α defined by its Mellin-Barnes representation can be obtained.
Because the auxiliary function L α is defined in Equation ( 20) as an inverse Mellin transform, its Mellin transform is given by the formula Now we derive a series representation of the fundamental solution G α by employing the Formula ( 16) and the general theory of the Mellin-Barnes integrals ( [19]).To arrive at a series representation, the contour of integration in the integral at the right-hand side of ( 16) has to be transformed to the infinite loop L +∞ starting and ending at +∞ and encircling all poles s k = 1/2 + k, k = 0, 1, 2, . . . of the gamma function Γ 1  2 − s and all poles s k = α + αk, k = 0, 1, 2, . . . of the gamma function Γ 1 − s α .For the sake of simplicity let us restrict ourselves to the case of the simple poles, i.e., we suppose that the conditions 1/2 + k = α + αn are fulfilled for all k, n ∈ IN, i.e., that the derivative order α cannot be represented in the form α = 1/2+k 1+n , k, n ∈ IN.In particular, evidently it is the case if α is not a rational number or a rational number in the form α = p 2q+1 , p, q ∈ IN.Taking into account the known formula the Jordan lemma and the Cauchy residue theorem provide us with the desired series representation of G α : where Using the reflection formula for the gamma function and some elementary trigonometric formulas, the series for Σ 1 and Σ 2 can be represented in the following shorter form: In the case α = 1, we get and thus the well-known formula for the fundamental solution of the one-dimensional diffusion equation.
Let us now introduce another auxiliary variable, namely, z = 1 y = x 2 4t .Then the formulas from above can be rewritten in the following form: where It follows from the last formula that the asymptotic behavior of G α as z → 0 (e.g., as x → 0 with a fixed t > 0 or t → +∞ with a fixed x = 0) is of a power law type: We remember the readers that all derivations above are valid only under the condition α > 1/2 that we assumed to hold true.
To get the asymptotic behavior of the fundamental solution G α as z → +∞ (e.g., as |x| → +∞ with a fixed t > 0 or t → 0 with a fixed x = 0) we again employ the Mellin-Barnes representation (16).This time, the contour of integration in the integral at the right-hand side of ( 16) has to be transformed to the infinite loop L −∞ starting and ending at −∞ and encircling all poles s k = −αk, k = 0, 1, 2, . . . of the gamma function Γ s α .In doing so we first get an asymptotic series , z → +∞ and then the asymptotic formula

Entropy Production Rate of the Alpha-Fractional Diffusion Process
For α = 1, the fundamental solution G α is a Gaussian probability density function evolving in time: As has been shown in [11], the fundamental solution G α to the α-fractional diffusion Equation (1) can be interpreted as a probability density function evolving in time for all values of α between 0 and 1, too.Let us mention that it follows from the asymptotic Formula (25) that the second spatial moment of the probability density function G α does not exist and the mean squared displacement of the diffusing particles in the framework of the diffusion process that is governed by the α-fractional diffusion Equation ( 1) is not finite.Thus the α-fractional diffusion Equation (1) describes a kind of an anomalous diffusion.Still, we show in this section that the entropy production rate of a diffusion process that is governed by the α-fractional diffusion Equation ( 1) is exactly the same as the one of the conventional diffusion process.
Let us start with some definitions and examples.In the case of a one-dimensional continuous random variable with the probability density function p(x), x ∈ X ⊆ I R, we adopt the Shannon definition of the entropy: where the constant k can be set to be equal to one without loss of the generality.The Shannon entropy (26) is a special case of the more general definitions by Mathai, Tsallis or Rényi and these entropies will be considered elsewhere.
Let us mention that the Shannon entropy of a Gaussian random variable defined by the probability density function has the form S(N (µ; Thus the entropy increases with the width σ 2 of the probability density function N (µ; σ 2 ), i.e., the broader the distribution (uncertainty of the event), the larger the entropy, so that the Shannon entropy can be interpreted as a measure of uncertainty of an event that is governed by a probability density function p(x).
When a probability density function is time-dependent, the entropy (26) depends on time, too: For such time-dependent random processes, the entropy production rate R defined by is a very important characteristic that can be interpreted as a natural measure of the irreversibility of a process.Say, in the case of a diffusion process that is described by the one-dimensional diffusion equation with the diffusion coefficient taken to be equal to one and that is therefore governed by the Gaussian distribution N (0; 2t), the Formula (27) leads to the following result: i.e., the entropy production rate is strictly positive for t > 0 and the diffusion process can be classified as an irreversible process.Otherwise, a wave propagation described by the wave equation is a reversible process with the entropy production rate equal to zero for t > 0.
It is worth mentioning that the entropy production rates for the time-and the space-fractional diffusion equations that were calculated in [3] and [5], respectively, depend on the derivative order α and increase with increasing of α from 1 (diffusion) to 2 (wave propagation) that results in the so called entropy production paradox (see [3,5] for attempts of resolving this paradox).
To calculate the entropy of the process governed by the one-dimensional α-fractional diffusion Equation ( 1), let us employ the representation (19) of its fundamental solution in terms of the auxiliary function L α .Substituting ( 19) into (28) and after some elementary transformations, we get the following chain of equalities: where To determine the entropy production rate of G α , the constant A α is determined in explicit form.The integral that defines A α can be interpreted as the Mellin transform of the auxiliary function L α at the point s = 0.Because the Formula (21) for the Mellin transform of L α was derived under the condition − min{1/2, α} < (s) < 0, the Mellin transform of L α at the point s = 0 will be calculated as the limit of the right-hand side of (21) as s → 0. Thus we get the following chain of equalities: This formula along with the Formula (30) leads to the following expression for the entropy production rate R(t) of the α-fractional diffusion process described by the Equation (1): This formula shows that R(t) does not depend on the equation order α and is exactly the same as the entropy production rate of the conventional one-dimensional diffusion equation.

Conclusions
In this paper, a special case of the time-and space-fractional diffusion equation has been considered, namely, the case when the quotient of the orders of the time-and space-fractional derivatives is equal to one half as it is the case in the conventional diffusion equation.This special choice of the derivative orders leads to several important consequences.On the one hand, the fundamental solution to the α-fractional equation can be expressed via an auxiliary function of the argument x 2 4t like in the case of the conventional diffusion equation.Another important property of the α-fractional diffusion equation is that its entropy production rate is exactly the same as the one of the diffusion process.Thus the α-fractional diffusion equation can be considered to be a "natural fractionalization" of the diffusion equation.On the other hand, the mean squared displacement of the diffusing particles governed by the α-fractional diffusion equation is infinite and thus this equation describes an anomalous diffusion process.