Growth Equation of the General Fractional Calculus

TL;DR

This paper investigates the growth behavior of solutions to a generalized fractional differential equation involving a convolutional derivative, extending classical Mittag-Leffler function solutions and analyzing their asymptotic properties.

Contribution

It introduces a general convolutional derivative and derives the growth equation for solutions, extending the classical fractional calculus framework.

Findings

Derived the growth equation for solutions involving the general convolutional derivative.

Analyzed the asymptotic behavior of solutions as time approaches infinity.

Extended the classical Mittag-Leffler function solutions to a broader fractional calculus context.

Abstract

We consider the Cauchy problem $(\mathbb D_{(k)} u)(t)=\lambda u(t)$, $u(0)=1$, where $\mathbb D_{(k)}$ is the general convolutional derivative introduced in the paper (A. N. Kochubei, Integral Equations Oper. Theory {\bf 71} (2011), 583--600), $\lambda >0$. The solution is a generalization of the function $t\mapsto E_\alpha (\lambda t^\alpha)$ where $0<\alpha <1$, $E_\alpha$ is the Mittag-Leffler function. The asymptotics of this solution, as $t\to \infty$, is studied.

Equations

$$u^{\prime }(t)=zu(t),t>0;u(0)=1,$$

$$({\mathbb{D}}_{(k)}u)(t)=\frac{d}{dt}\underset{0}{\overset{t}{\int }}k(t-\tau )u(\tau )d\tau -k(t)u(0)$$

$$({\mathbb{D}}_{(k)}u)(t)=-\lambda u(t),t>0;u(0)=1,$$

$$({\mathbb{D}}_{(k)}w)(t,x)=\Delta w(t,x),t>0,x\in {\mathbb{R}}^n;w(0,x)=w_0(x),$$

$$k(t)=\frac{t^{-\alpha }}{\Gamma (1-\alpha )},t>0,$$

$$({\mathbb{D}}_{(k)}u)(t)=\lambda u(t),t>0;u(0)=1;$$

$$\mathcal{K}(p)=\underset{0}{\overset{\infty }{\int }}{e}^{-pt}k(t)dt.$$

$$\mathcal{K}(p)\to \infty ,\text{ as p→0};\mathcal{K}(p)\to 0,\text{ as p→∞};$$

$$p\mathcal{K}(p)\to 0,\text{ as p→0};p\mathcal{K}(p)\to \infty ,\text{ as p→∞};$$

$$\psi (z)=\frac{a}{z}+b+\underset{0}{\overset{\infty }{\int }}\frac{1}{z+t}\sigma (dt)$$

$$\underset{0}{\overset{\infty }{\int }}(1+t{)}^{-1}\sigma (dt)<\infty .$$

$$\mathcal{K}(p)=\underset{0}{\overset{\infty }{\int }}\frac{1}{z+t}\sigma (dt).$$

$$k(s)=\underset{0}{\overset{\infty }{\int }}{e}^{-ts}\sigma (dt),0<s<\infty .$$

$$(-1{)}^{n-1}{f}^{(n)}(z)\ge 0\text{for all n≥1,z>0}.$$

$$f(0+)=\underset{(0,\infty )\ni z\to 0}{\lim }f(z).$$

$$\sqrt{\frac{1+\cos \phi }{2}}\Phi (|p|)\le |\Phi (p)|\le \sqrt{\frac{2}{1+\cos \phi }}\Phi (|p|),\phi =\arg p.$$

$$g(s,p)=\mathcal{K}(p)e^{-s\Phi (p)},s>0,p>0.$$

$$g(s,p)=\underset{0}{\overset{\infty }{\int }}{e}^{-pt}G(s,t)dt.$$

$$\underset{0}{\overset{\infty }{\int }}G(s,t)ds=1\text{for each t}.$$

$$\Phi (p)\widetilde{u_{\lambda }}(p)-\lambda \widetilde{u_{\lambda }}(p)=\mathcal{K}(p),$$

$$\widetilde{u_{\lambda }}(p)=\frac{\mathcal{K}(p)}{\Phi (p)-\lambda },\text{if Φ(p)>λ.}$$

$$E(t,\lambda )=\underset{0}{\overset{\infty }{\int }}{e}^{\lambda s}G(s,t)ds,t>0.$$

$$\underset{0}{\overset{\infty }{\int }}{e}^{\lambda s}ds\underset{0}{\overset{\infty }{\int }}{e}^{-pt}G(s,t)dt=\underset{0}{\overset{\infty }{\int }}{e}^{\lambda s}g(s,\lambda )ds=\frac{\mathcal{K}(p)}{\Phi (p)-\lambda }$$

$$\Phi (p)=\underset{0}{\overset{\infty }{\int }}\frac{p}{p+t}\sigma (dt),$$

$$p_0(x+y)>p_0(x)+p_0(y)\text{for any x,y>0}.$$

$$\underset{t\in {\Sigma }_v}{\sup }\left|e^{-p_{0}t}{u}_{\lambda }(t)\right|<\infty ,p_0=p_0(\lambda ).$$

$$\underset{p\in p_0+{\Sigma }_{\rho +\frac{\pi }{2}}}{\sup }\left|(p-p_0)\widetilde{u_{\lambda }}(p)\right|<\infty .$$

$$\underset{1}{\overset{\infty }{\int }}\frac{ds}{s\Phi (s)}<\infty .$$

$$u_{\lambda }(t)=\frac{\lambda }{{\Phi }^{\prime }(p_0(\lambda ))p_0(\lambda )}{e}^{p_0(\lambda )t}+o(e^{p_0(\lambda )t}),t\to \infty .$$

$$\widetilde{u_{\lambda }}(p)=\frac{1}{p}\left(1+\frac{\lambda }{\Phi (p)-\lambda }\right).$$

Full text

Keywords generalized fractional derivatives, growth equation; Mittag-Leffler function

1 Introduction

In several models of dynamics of complex systems the time evolution for observed quantities has exponential asymptotics of two possible types. In the simplest cases these asymptotics are related with the solutions to the equations

[TABLE]

where we will consider positive and negative $z$ separately. For $z<0$ (the relaxation equation) the solution decays to zero. In particular models as, e.g., Glauber stochastic dynamics in the continuum, this corresponds to an exponential convergence to an equilibrium, see [16]. The case $z>0$ also may appear in applications. We can mention the contact model in the continuum where for the mortality below a critical value the density of population will grow exponentially fast [18, 17], as well as models of economic growth.

On the other hand, the observed behavior of specific physical and biological systems show an emergence of other time asymptotics which may be far from the exponential decay or growth. An attempt to obtain other relaxation characteristics is related with a use of generalized time derivatives in dynamical equations (see [14, 15]). In this way we may produce a wide spectrum of possible asymptotics to reflect a demand coming from applications [23].

The general fractional calculus introduced in [12] is based on a version of the fractional derivative, the differential-convolution operator

[TABLE]

where $k$ is a nonnegative locally integrable function satisfying additional assumptions, under which

(A) The Cauchy problem

[TABLE]

where $\lambda>0$, has a unique solution which is completely monotone;

(B) The Cauchy problem

[TABLE]

is solvable (under appropriate conditions for $w_{0}$) and possesses a fundamental solution, a kernel with the property of a probability density.

A class of functions $k$, for which (A) and (B) hold, was found in [12] and is described below. The simplest example is

[TABLE]

where $0<\alpha<1$, and for this case $\mathbb{D}_{(k)}$ is the Caputo-Djrbashian fractional derivative $\mathbb{D}^{(\alpha)}$. Another subclass is the one of distributed order derivatives; see [10] for the details.

Note that for the case where $k$ has the form (2), the solution of (1) is $u(t)=E_{\alpha}(-\lambda t^{\alpha})$ where $E_{\alpha}$ is the Mittag-Leffler function; see [9], Lemma 2.23 (page 98). This solution has a slow decay at infinity, due to the asymptotic property of the Mittag-Leffler function; see [7]. Note that using particular classes of fractional derivatives we observe several specific asymptotics for the solution of the equation (1) with $\lambda>0$. Some results in this direction were already obtained in [10, 11, 13]. A more detailed analysis of this problem will be performed in a forthcoming paper.

In this paper, we consider the Cauchy problem with the opposite sign in the right-hand side, that is

[TABLE]

as before, $\lambda>0$. In the case (2), we have $u(t)=E_{\alpha}(\lambda t^{\alpha})$ (see [9], Lemma 2.23 (page 98)), and, due to the well-known asymptotics of $E_{\alpha}$ [7], this is a function of exponential growth. The existence and uniqueness of an absolutely continuous solution of (3) follows from the results of [21] dealing with more general nonlinear equations. Here we study the asymptotic behavior of the solution of (3). Functions of this kind can be useful for fractional macroeconomic models with long dynamic memory; see [22] and references therein. Let us explain this in a little greater detail.

In modern macroeconomics, the most important are so-called growth models, which in the mathematical sense are reduced (for linear models) to the equation $u^{\prime}(t)=\lambda u(t)+f(t)$ with $\lambda>0$. In economics, an important role is played by processes with a distributed lag starting with Phillips’ works [19] (see also [1]) and long memory, starting with Granger’s work [8] (see also [3]).

If we assume the presence of effects of distributed lag (time delay) or fading memory in economic processes, then the fractional generalization of the linear classical growth models can be described by the fractional differential equation $D^{\alpha}u(t)=\lambda u(t)+f(t)$ with $\lambda>0,\alpha>0$. The fractional generalizations of well-known economics models have been first proposed for the Caputo-Djrbashian fractional derivative $D^{\alpha}$. Solving the problem in a more general case will allow us to describe accurately the conditions on the operator kernels (the memory functions), under which equations for models of economic growth with memory have solutions.

In general fractional calculus, which was proposed in [12] (see also [14]), the case $\lambda>0$ is not considered. The growth equation was considered in [10] for the special case of a distributed order derivative, where it was proved that a smooth solution exists and is monotone increasing.

In this article we propose correct mathematical statements for growth models with memory in more general case, for the general fractional derivative $\mathbb{D}_{(k)}$ with respect to the time variable. Their application can be useful for mathematical economics for the description of processes with long memory and distributed lag.

Note that the technique used below was developed initially in [13] for the use in the study of intermittency in fractional models of statistical mechanics.

2 Preliminaries

Our conditions regarding the function $k$ will be formulated in terms of its Laplace transform

[TABLE]

Denote $\Phi(p)=p\mathcal{K}(p)$.

We make the following assumptions leading to (A) and (B) (see [12]).

(*)

The Laplace transform (4) exists for all positive numbers $p$. The function $\mathcal{K}$ belongs to the Stieltjes class $\mathcal{S}$, and

[TABLE]

[TABLE]

Recall that the Stieltjes class consists of the functions $\psi$ admitting the integral representation

[TABLE]

where $a,b\geq 0$, $\sigma$ is a Borel measure on $[0,\infty)$, such that

[TABLE]

For a detailed exposition of the theory of Stieltjes functions including properties of the measure $\sigma$ see [20], especially Chapters 2 and 6.

In particular, for the Stieltjes function $\mathcal{K}$, the limit conditions (5), (6) imply the representation

[TABLE]

We can also write [12] that

[TABLE]

The function $\Phi$ belongs to the class $\mathcal{C}\mathcal{B}\mathcal{F}$ of complete Bernstein functions, a subclass of the class $\mathcal{B}\mathcal{F}$ of Bernstein functions. Recall that a function $f:\ (0,\infty)\to\mathbb{R}$ is called a Bernstein function, if $f\in C^{\infty}$, $f(z)\geq 0$ for all $z>0$, and

[TABLE]

so that the derivative of $f$ is completely monotone. A function $f$ belongs to $\mathcal{C}\mathcal{B}\mathcal{F}$, if it has an analytic continuation to the cut complex plane $\mathbb{C}\setminus(-\infty,0]$, such that $\operatorname{Im}z\cdot\operatorname{Im}f(z)\geq 0$, and there exists the real limit

[TABLE]

Both the classes $\mathcal{B}\mathcal{F}$ and $\mathcal{C}\mathcal{B}\mathcal{F}$ admit equivalent descriptions in terms of integral representations; see [20].

Below we will need the following inequality for complete Bernstein functions (Proposition 2.4 in [4]), valid, in particular, for the function $\Phi$. For any $p$ outside the negative real semi-axis, we have

[TABLE]

Solutions of the Cauchy problem (3) and a similar problem with the classical first order derivative are connected by the subordination identity (see [12]; for the case of the Caputo-Djrbashian derivative see [5]), an integral transformation with the kernel $G(s,t)$ constructed as follows.

Consider the function

[TABLE]

It is proved [12] that $g$ is a Laplace transform in the variable $t$ of the required kernel $G(s,t)$, that is

[TABLE]

$G$ is nonnegative, and

[TABLE]

3 Cauchy problem for the growth equation

Let us consider the Cauchy problem (3). If $u_{\lambda}(t)$ is its solution whose Laplace transform $\widetilde{u_{\lambda}}(p)$ exists for some $p$, then it follows from properties of the Laplace transform [6] that

[TABLE]

hence

[TABLE]

On the other hand, consider the function

[TABLE]

The existence of the integral in (12) for almost all $t>0$ is, by the Fubini-Tonelli theorem, a consequence of the absolute convergence of the repeated integral

[TABLE]

where $p>0$ is such that $\Phi(p)>\lambda$.

The above calculation shows that $E(t,\lambda)=u_{\lambda}(t)$, the solution of (3), and the identity (12) provides an integral representation of this solution.

A more detailed analysis of its properties is based on the analytic properties of the Stieltjes function $\mathcal{K}$ or, equivalently, of the complete Bernstein function $\Phi$; in particular, we use the representation

[TABLE]

which follows from (7). The measure $\sigma$ satisfies (8).

Since $\Phi$ is a Bernstein function, its derivative $\Phi^{\prime}$ is completely monotone. By our assumptions, $\Phi$ is not a constant function, so that $\Phi^{\prime}$ is not the identical zero. It follows from Bernstein’s description of completely monotone functions that $\Phi^{\prime}(p)\neq 0$ for any $p>0$ (see Remark 1.5 in [20]). Therefore $\Phi$ is strictly monotone, and for each $z>0$, there exists a unique $p_{0}=p_{0}(z)>0$, such that $\Phi(p_{0})=z$. The inequality $\Phi(p)>z$ is equivalent to the inequality $p>p_{0}(z)$. Since $\Phi$, as a complete Bernstein function, preserves the open upper and lower half-planes (in fact, this follows from (13), we have $\Phi(p)\neq z$ for any nonreal $p$.

It is proved in [13] that the function $p_{0}(z)$, $z>0$, is strictly superadditive, that is

[TABLE]

Proposition**.**

The solution $u_{\lambda}(t)$ of the Cauchy problem (3) admits a holomorphic continuation in the variable $t$ to a sector $\Sigma_{v}=\left\{re^{i\theta}:\ r>0,-v<\theta<v\right\}$, $0<v<\frac{\pi}{2}$, and

[TABLE]

Proof. It follows from (11) and (13) that the Laplace transform $\widetilde{u_{\lambda}}(p)$ is holomorphic in $p$ on any sector $p_{0}+\Sigma_{\rho+\frac{\pi}{2}}$, $0<\rho<\frac{\pi}{2}$. In addition,

[TABLE]

Now the assertion is implied by (15) and the duality theorem for holomorphic continuations of a function and its Laplace transform; see Theorem 2.6.1 in [2]. $\qquad\blacksquare$

Now we are ready to formulate and prove our main result.

Theorem**.**

Let the assumptions $(*)$ hold, and in addition,

[TABLE]

Then

[TABLE]

Proof. The representation (11) can be written as

[TABLE]

This implies the representation of $u_{\lambda}$ as $u_{\lambda}(t)=1+B_{\lambda}(t)$ where $B_{\lambda}$ has the Laplace transform

[TABLE]

for such $p$ that $\Phi(p)>\lambda$.

Using the inequality (9) we find that $|\Phi(p)|\geq\frac{1}{\sqrt{2}}\Phi(|p|)$ on any vertical line $\{p=\gamma+i\tau,\tau\in\mathbb{R}\}$ where $\gamma>p_{0}$. By our assumption (16), $\widetilde{B_{\lambda}}$ is absolutely integrable on such a line. In addition, it follows from (6) that $\widetilde{B_{\lambda}}(p)\to 0$, as $p\to\infty$ in the half-plane $\operatorname{Re}p>p_{0}$. These properties make it possible (see Theorem 28.2 in [6]) to write the inversion formula

[TABLE]

Denote

[TABLE]

where $0<r<p_{0}$. Then

[TABLE]

by virtue of (9), (16) and the Riemann-Lebesgue theorem.

On the other hand, we may write

[TABLE]

where the contour $\Lambda_{+}$ consists of the vertical rays $\{\operatorname{Re}p=r,\operatorname{Im}p\geq R\}$, $\{\operatorname{Re}p=\gamma,\operatorname{Im}p\geq R\}$, and the horizontal segment $\{r\leq\operatorname{Re}p\leq\gamma,\operatorname{Im}p=R\}$ ($R>0$), $\Lambda_{-}$ is a mirror reflection of $\Lambda_{+}$ with respect to the real axis, $\Lambda_{0}$ is the finite rectangle consisting of the vertical segments $\{\operatorname{Re}p=r,|\operatorname{Im}p|\leq R\}$, $\{\operatorname{Re}p=\gamma,|\operatorname{Im}p|\leq R\}$, and the horizontal segments $\{r\leq\operatorname{Re}p\leq\gamma,\operatorname{Im}p=\pm R\}$.

We have

[TABLE]

due to the Cauchy theorem, absolute integrability of the integrand on the vertical rays (see (16)) and the estimate

[TABLE]

where $\Pi_{h}=\{r\leq\operatorname{Re}p\leq\gamma,\operatorname{Im}p=h\}$, $h>R$. In a similar way, we prove that

[TABLE]

Due to the inequality $\Phi^{\prime}(p_{0})\neq 0$, there exists a complex neighborhood $W$ of the point $\lambda=\Phi(p_{0})$, in which $\Phi$ has a single-valued holomorphic inverse function $p=\psi(w)$, so that $\Phi(\psi(w))=w$ and $p_{0}=\psi(\lambda)$. In the above arguments, the numbers $r,\gamma,R$ were arbitrary. Now we choose $R$ and $\gamma-r$ so small that that the curvilinear rectangle $\Phi(\Lambda_{0})$ lies inside $W$. Making the change of variables $p=\psi(w)$ and using the Cauchy formula we find that

[TABLE]

Together with (18), this implies the required asymptotic relation (17). $\qquad\blacksquare$

Examples. 1) In the case (2) of the Caputo-Djrbashian fractional derivative of order $0<\alpha<1$, we have $u_{\lambda}(t)=E_{\alpha}(\lambda t^{\alpha})$, $\Phi(p)=p^{\alpha}$, and the condition (16) is satisfied. Here the asymptotics (17) coincides with the one given by the principal term of the asymptotic expansion of the Mittag-Leffler function. The above proof is different from the classical proof of the latter (see [7]).

2. Let us consider the case of a distributed order derivative with a weight function $\mu$, that is

[TABLE]

Suppose that $\mu\in C^{2}[0,1]$, $\mu(1)\neq 0$. In this case [10],

[TABLE]

and under the above assumptions,

[TABLE]

The condition (16) is satisfied, and our asymptotic result (17) is applicable in this case.

Acknowledgments

The authors are grateful to V. E. Tarasov for calling our attention to the problem studied in this paper, as well as for his advices regarding its text.

The work of the first-named author was funded in part under the budget program of Ukraine No. 6541230 “Support to the development of priority research trends” and under the research work ”Markov evolutions in real and p-adic spaces” of the Dragomanov National Pedagogical University of Ukraine.