This study investigates the nonlinear Hamiltonian amplitude equation by using two analytical techniques, namely; the extended sinh-Gordon equation expansion method, and the extended rational sine-cosine/sinh-cosh methods. Some important wave solutions are successfully constructed, such as dark, bright, combined dark-bright, singular solitons, periodic and singular periodic wave solutions. The physical features of the acquired solutions are plotted to depict the clear dynamical behaviour of the reported results. All the acquired solutions have satisfied the original equation.
Citation: Abdullahi Yusuf, Tukur A. Sulaiman, Mustafa Inc, Sayed Abdel-Khalek, K. H. Mahmoud. M−truncated optical soliton and their characteristics to a nonlinear equation governing the certain instabilities of modulated wave trains[J]. AIMS Mathematics, 2021, 6(9): 9207-9221. doi: 10.3934/math.2021535
[1] | Yaning Li, Mengjun Wang . Well-posedness and blow-up results for a time-space fractional diffusion-wave equation. Electronic Research Archive, 2024, 32(5): 3522-3542. doi: 10.3934/era.2024162 |
[2] | Peng Gao, Pengyu Chen . Blowup and MLUH stability of time-space fractional reaction-diffusion equations. Electronic Research Archive, 2022, 30(9): 3351-3361. doi: 10.3934/era.2022170 |
[3] | Mustafa Aydin, Nazim I. Mahmudov, Hüseyin Aktuğlu, Erdem Baytunç, Mehmet S. Atamert . On a study of the representation of solutions of a Ψ-Caputo fractional differential equations with a single delay. Electronic Research Archive, 2022, 30(3): 1016-1034. doi: 10.3934/era.2022053 |
[4] | Mehmet Ali Özarslan, Ceren Ustaoğlu . Extended incomplete Riemann-Liouville fractional integral operators and related special functions. Electronic Research Archive, 2022, 30(5): 1723-1747. doi: 10.3934/era.2022087 |
[5] | Melih Cinar, Ismail Onder, Aydin Secer, Mustafa Bayram, Abdullahi Yusuf, Tukur Abdulkadir Sulaiman . A comparison of analytical solutions of nonlinear complex generalized Zakharov dynamical system for various definitions of the differential operator. Electronic Research Archive, 2022, 30(1): 335-361. doi: 10.3934/era.2022018 |
[6] | Bin Han . Some multivariate polynomials for doubled permutations. Electronic Research Archive, 2021, 29(2): 1925-1944. doi: 10.3934/era.2020098 |
[7] | Yong Zhou, Jia Wei He, Ahmed Alsaedi, Bashir Ahmad . The well-posedness for semilinear time fractional wave equations on RN. Electronic Research Archive, 2022, 30(8): 2981-3003. doi: 10.3934/era.2022151 |
[8] | María Guadalupe Morales, Zuzana Došlá, Francisco J. Mendoza . Riemann-Liouville derivative over the space of integrable distributions. Electronic Research Archive, 2020, 28(2): 567-587. doi: 10.3934/era.2020030 |
[9] | Yuanfei Li . A study on continuous dependence of layered composite materials in binary mixtures on basic data. Electronic Research Archive, 2024, 32(10): 5577-5591. doi: 10.3934/era.2024258 |
[10] | Yadan Shi, Yongqin Xie, Ke Li, Zhipiao Tang . Attractors for the nonclassical diffusion equations with the driving delay term in time-dependent spaces. Electronic Research Archive, 2024, 32(12): 6847-6868. doi: 10.3934/era.2024320 |
This study investigates the nonlinear Hamiltonian amplitude equation by using two analytical techniques, namely; the extended sinh-Gordon equation expansion method, and the extended rational sine-cosine/sinh-cosh methods. Some important wave solutions are successfully constructed, such as dark, bright, combined dark-bright, singular solitons, periodic and singular periodic wave solutions. The physical features of the acquired solutions are plotted to depict the clear dynamical behaviour of the reported results. All the acquired solutions have satisfied the original equation.
It is well-known that in the theory of linear operators there are three different points of view to represent semigroups [1,2,3]. One is the classical algebraic approach, i.e., a strongly continuous operator function T(⋅;A):R+→L(X) is called a C0-semigroup if T(0;A)=I and T(t;A)T(s;A)=T(t+s;A) for all t,s∈R+. The second approach involves the abstract Cauchy problem
{ddtu(t)=Au(t),t∈R+,u(0)=x, | (1.1) |
that is, A:D(A)⊆X→X generates a C0-semigroup T(t;A),t∈R+ if and only if there exists a unique mild solution u for every initial value x∈X. Furthermore, in this case, the closed-form representation of a mild solution is u(t)=T(t;A)x for all t∈R+ and x∈X. And the third approach is related to the Laplace integral transform, i.e., a strongly continuous operator function T(⋅;A):R+→L(X) is a C0-semigroup if and only if there exists an infinitesimal generator A:D(A)⊆X→X and ω∈R such that (ω,∞)∈ρ(A) and Re(λ)>ω,
R(λ;A)x=∞∫0e−λtT(t;A)xdt,x∈X. | (1.2) |
The same behaviour applies to strongly continuous cosine families [4,5,6]. In agreement with the purely algebraic approach, a strongly continuous operator function C(⋅:A):R→L(X) is called a cosine family if C(0;A)=I and C(t+s;A)+C(t−s;A)=2C(t;A)C(s;A) for all t,s∈R and the associated sine family S(⋅;A):R→L(X) is given by S(t;A)x=∫t0C(s;A)xds for each x∈X and t∈R. Based on the second approach, the following abstract initial value problem
{d2dt2u(t)=Au(t),t∈R,u(0)=x,ddtu(0)=y, | (1.3) |
has a unique mild solution u for each initial values x,y∈X if and only if A:D(A)⊆X→X is an infinitesimal generator of strongly continuous cosine and sine families C(t;A),t∈R and S(t;A),t∈R, respectively. Moreover, in this case, the closed-form representation of a mild solution is u(t)=C(t;A)x+S(t;A)y for all t∈R and x,y∈X. According to the last approach strongly continuous operator functions C(⋅;A):R+→L(X) and S(⋅;A):R+→L(X) are strongly continuous cosine and sine families if and only if there exists an infinitesimal generator A:D(A)⊆X→X and ω∈R such that (ω,∞)∈ρ(A) and Re(λ)>ω,
λR(λ2;A)x=∞∫0e−λtC(t;A)xdt,x∈X, | (1.4) |
R(λ2;A)x=∞∫0e−λtS(t;A)xdt,x∈X. | (1.5) |
Furthermore, we would like to draw the readers' attention to the research work of Vasil'ev, Krein and Piskarev [7,8] on operator theory, in particular on semigroup theory and cosine operator families in abstract spaces. These fundamental surveys on operator theory can serve as a fairly complete source of information on abstract differential equations and methods for solving them in recent years.
In recent years there has been a considerable increase in interest to fractional calculus and fractional differential equations, where the order of differentiation is permitted to be any real (or even complex) number, not just a natural number. This is largely due to the applications of fractional differential equations to problems in various fields of physics and engineering, namely electrical circuits [9], electrical power systems [10,11,12], vibration theory [13] and so on.
One of the emerging branches of fractional calculus and fractional-order differential equations are the fractional semigroups and fractional cosine (sine) families and their applications to the theory of fractional evolution equations, i.e., evolution equations in which the integer-order derivatives with respect to time are replaced by the derivatives of fractional-order. In [14], Bazhlekova has studied the following abstract Cauchy problem for a fractional evolution equation of Caputo type in a Banach space X and introduced the solution operator for (1.6) as follows:
{(Dα0+u)(t)=Au(t),t>0,u(0)=x,dkdtku(0)=0,k=1,2,…,m−1, | (1.6) |
where m−1<α<m, m∈N.
Definition 1.1. A family Tα(t;A), t∈R+ of bounded linear operators on X is called a solution operator for (1.6) if the following three conditions are satisfied:
(i) Tα(t;A) is strongly continuous for any t∈R+ and Tα(0;A)=I;
(ii) Tα(t;A)D(A)⊂D(A) and ATα(t;A)x=Tα(t;A)Ax for all x∈D(A) and t∈R+;
(iii) u(t)=Tα(t;A)x is a strong solution of (1.6) for every x∈D(A).
Note that the solution operator for an abstract Cauchy problem (1.6) was systematically studied in the thesis of Bazhlekova [14], and the results obtained generalise some facts of C0-semigroups (in case of 0<α<1) and cosine families (in case of 1<α<2). Furthermore, Lizama in [15] has reviewed the fundamental theory of solution operators associated to abstract linear evolution equations of fractional-order and provided their basic results concerning generation, analyticity and inversion problems. In [16], Jigen and Kexue have developed the notation of fractional semigroup to characterize solution operator for fractional abstract linear systems (1.6) of order 0<α<1 and discussed the well-posedness of the abstract Cauchy problem (1.6). It is shown that in [16] the problem (1.6) is well-posed if and only if its coefficient operator A generates a fractional semigroup. The same authors have also investigated similar results in [17] for Riemann-Liouville type fractional evolution equations of order α∈(0,1).
In terms of fractional anologue of strongly continuous cosine families, Kexue and Jigen in [18] and Mei et al. in [19] have developed an operator theory to study abstract Cauchy problems of order α∈(1,2) in terms of Caputo and generalized type fractional differentiation operators, respectively and proposed a generation theorem for exponentially bounded fractional solution operators-fractional strongly continuous sine families. Furthermore, in [20], Mei et al. have studied a new characteristic property of the Mittag-Leffler function Eα(atα), a∈R with α∈(1,2) and motivated by this property developed a new notion, named α-order cosine function. It is interesting to note that in [21] Mei et al. have introduced a novel one-parameter algebraic functional equation for fractional resolvent families, and using this functional equation, all fractional resolvent families except C0-semigroups can be obtained from this equation.
In addition, a subordination principle can be used to analyse fractional evolution equations, which is a useful tool to connect solution operators in the classical and fractional sense. Using a subordination formula, it is possible to construct new solutions from known solutions, e.g., solutions of fractional-order evolution equations from the solutions of the classical diffusion or wave equations (see [22,23,24,25]).
A strong inspiration for the study of perturbation theory for strongly continuous operator families, particularly, for cosine families comes from the fact that they have been proven to be useful tools for evolution equations in modeling many physical phenomena. The perturbed wave equation is a typical linear hyperbolic second-order partial differential equation which naturally arises when modeling phenomena of continuum mechanics such as sound, light, water or other kind of waves in acoustics, electro-magnetics, elasticity and fluid dynamics, etc. (see [26]).
The perturbation theory of linear operators in a Banach space has been explored to a considerable extent, most notably by Phillips [27], Travis & Webb [28] and Lutz [29]. In [27], Phillips intended to perturb the infinitesimal generator by adding to it a linear bounded operator and investigating some perturbation results for infinitesimal generators of strongly continuous C0-semigroups. In [28], Travis and Webb have proposed sufficient conditions for perturbed cosine operator families. In [29], Lutz has first studied the implication for the homogeneous initial value problem associated with the infinitesimal generator of the cosine function generated by bounded time-varying perturbations:
{d2dt2u(t)=(A+B(t))u(t),t∈R,u(0)=u0,u′(0)=u1, | (1.7) |
where B:R+→L(X) is a continuously differentiable function on R+.
In her dissertation, Bazhlekova showed that perturbation theory does not work in general for the solution operators of the perturbed abstract Cauchy problem with 0<α<1 (see counterexample in [14]). In the special case, a treatment of the perturbation problem from the point of view of the subjacent fractional abstract Cauchy problem in the case of 0<α<1 was given by El-Borai [30]. In the general case, however, perturbation properties are true for the fractional analogue of an abstract initial value problem (1.7) with 1<α<2 which has been established by Bazhlekova in [31]. She has proposed to uniquely determine a classical solution of the following homogeneous abstract Cauchy problem for the fractional evolution equation with Caputo derivative by time-dependent perturbations:
{(Dα0+u)(t)=(A+B(t))u(t),t>0,u(0)=u0,u′(0)=0, |
where B:R+→L(X) is a continuous function on R+.
Furthermore, the fractional analogue of the inhomogeneous abstract problem (1.7) was developed by Ahmadova et al. in [32] and the results presented in [32] extend those of [29,31] in several aspects. The authors in [32] have shown that the inhomogeneous abstract Cauchy problem with the infinitesimal generator A of fractional cosine families remains uniformly well-posed under bounded time-varying perturbations:
{(Dα0+u)(t)=(A+B(t))u(t)+h(t),t>0,u(0)=u0,u′(0)=u1, |
where B:R+→L(X) and h:R+→X are continuously differentiable functions on R+.
While there are several papers on operator theory for fractional evolution equations (for example, see [15,16,17,18,33]), there are few papers [30,31,32] dealing with perturbation theory for fractional abstract Cauchy problems. The lack of perturbation properties of the fractional strongly contionuous operator functions motivates us to develop new results on perturbation theory for evolution equations of fractional-order. The pioneering work on this subject in the classical sense was done by Phillips [27], Lutz [29] and Travis & Webb [28] and our development follows these approaches. Therefore, the main contributions of our work are described in detail below:
1) We propose some new perturbation results for fractional strongly continuous cosine families in a Banach space X;
2) We establish sufficient conditions such that A is the infinitesimal generator of a fractional strongly continuous cosine (sine) family in a Banach space X, and B is a bounded linear operator in X, then A+B is also the infinitesimal generator of a fractional strongly continuous cosine (sine) family in X;
3)Depending on commutativity condition of linear bounded operators, we obtain the elegant analytical representation formulas for perturbed Mittag-Leffler type functions generated by A+B where A,B∈L(X) are permutable and nonpermutable linear bounded operators;
4) We present an example in the context of a one-dimensional fractional perturbed wave equation to demonstrate the applicability of our theoretical results and we give some comparisons with the existing literature.
The paper contains significant innovations on the perturbation theory for fractional cosine families and is organized as follows. Section 2 is a preparatory section in which we recall the main definitions and results from functional analysis, fractional calculus, and operator theory of linear operators. In Section 3, we investigate some perturbation properties for fractional strongly continuous cosine families. We intend to perturb the infinitesimal generator by adding to it a linear bounded operator. Furthermore, we propose the elegant representation formulas for perturbed fractional uniformly continuous cosine families in a Banach space X. In Section 4, we present an example in the context of one-dimensional perturbed fractional wave equation to demonstrate the applicability of our theoretical results and we give some comparisons with existing ones. In Section 5, we give an overview of our main contributions and point out some open problems in the spirit of this research.
We start this section by briefly introducing the fundamentals of fractional calculus, the essential structure of an operator theory for strongly continuous cosine families and their fractional analogues. For the more salient details on these subjects, see the textbooks [2,34,35].
Let R+=[0,∞) and N denote the set of natural numbers with N0=N∪{0}. Let X be a Banach space equipped with the norm ‖⋅‖. We donote by L(X) the Banach algebra of all bounded linear operators on X and becomes a Banach space with regard to the norm ‖T‖=sup{‖Tx‖:‖x‖≤1}, for any T∈L(X). Let D(A) be the domain of A and ρ(A) be the resolvent set of A. The identity and zero operators on X are denoted by I∈L(X) and 0∈L(X), respectively.
We will use the following functional spaces [2] through the paper:
● C(R+,X) denotes the Banach space of continuous X-valued functions g:R+→X equipped with an infinity norm ‖g‖∞=supt∈R+‖g(t)‖;
● Cn(R+,X), n∈N denotes the Banach space of n-times continuously differentiable X-valued functions defined by
Cn(R+,X)={g∈Cn(R+,X):g(n)∈C(R+,X),n∈N} |
and equipped with an infinity norm ‖g‖∞=n∑k=0supt∈R+‖g(k)(t)‖. In addition, Cn(R+,X)⊂C(R+,X), n∈N;
● Lp(R+,X) denotes the Banach space of all (equivalence classes of) Bochner-measurable X-valued functions g:R+→X which are integrable and normed by ‖g‖Lp=(∫R+‖g(s)‖pds)1/p<∞ where 1≤p<∞;
● The Sobolev spaces can be defined in the following way [14] for n∈N and 1≤p<∞:
Wn,p(R+,X):={ϕ:∃h∈Lp(R+,X):ϕ(t)=n−1∑k=0cktkk!+tn−1(n−1)!∗h(t),t∈R+}.
Note that ck=ϕ(n)(0) and ϕ(n)(t)=h(t).
Let f:R+→R be an integrable scalar-valued function and let g:R+→X be a continuous X-valued function. We denote by
(f∗g)(t)=t∫0f(t−s)g(s)ds,t∈R+. |
the convolution operator of f and g. Furthermore, if T:R+→L(X) is a strongly continuous operator-valued map, we define
(T∗g)(t)=t∫0T(t−s)g(s)ds=t∫0T(s)g(t−s)ds,t∈R+. | (2.1) |
Definition 2.1. The Riemann-Liouville fractional integral of order α>0 for a function f∈L1(R+,X) is defined by
(Iαtf)(t)=(gα∗f)(t)=1Γ(α)t∫0(t−s)α−1f(s)ds,t>0. | (2.2) |
By setting (I0tf)(t)=f(t). For the sake of brevity we use the following notation for α≥0:
gα(t)={tα−1Γ(α),t>0,0,t≤0, |
where Γ:R+→R is the well-known Euler's gamma function. Note that g0(t)=0, since 1Γ(0)=0. These functions satisfy the semigroup property
(gα∗gβ)(t)=gα+β(t). | (2.3) |
Moreover, the Riemann-Liouville fractional integral operators {Iαt}α≥0 satisfy the semigroup property
IαtIβt=Iα+βt,α,β≥0. | (2.4) |
Definition 2.2. The Caputo fractional derivative of a function f∈Cn(R+,X) with fractional order n−1<α<n, n∈N is defined by
(Dα0+f)(t)=In−αtDntf(t),t>0. |
In particular, for α∈(1,2), the definition is given by for f∈C2(R+,X):
(Dα0+f)(t)=(g2−α∗f′′)(t)=1Γ(2−α)t∫0(t−s)1−αf′′(s)ds,t>0. |
Remark 2.1. Since f is an abstract function with values in X, the integrals which appear in Definition 2.1 and 2.2 are taken in Bochner's sense.
The Mittag-Leffler function is a natural generalization of the exponential function, first proposed as a single parameter function of one variable by using an infinite series [36].
Definition 2.3. [36] The classical Mittag-Leffler function is defined by
Eα(t)=∞∑k=0tkΓ(kα+1),α>0,t∈R. | (2.5) |
The two-parameter Mittag-Leffler function is given by
Eα,β(t)=∞∑k=0tkΓ(kα+β),α>0,β∈R,t∈R. | (2.6) |
It is important to note that
Eα,1(t)=Eα(t),E1(t)=exp(t),t∈R. |
Strongly continuous cosine and sine families of linear operators
Definition 2.4. See reference[37] 1) A one parameter family {C(t;A),t∈R}∈L(X) into itself is called a strongly continuous cosine family in the Banach space X if and only if
● C(0;A)=I;
● C(s+t;A)+C(s−t;A)=2C(s;A)C(t;A) for all s,t∈R;
● C(t;A)x is continuous in t on R for each fixed x∈X.
2) The corresponding strongly continuous sine family {S(t;A),t∈R}∈L(X) is defined by
S(t;A)x=t∫0C(s;A)xds,x∈X,t∈R. | (2.7) |
3) The linear operator A:D(A)⊆X→X defined by
Ax=2limt→0+C(t;A)x−xt2,x∈D(A)={x∈X:C(⋅;A)x∈C2(R,X)}, | (2.8) |
is called the infinitesimal generator of a strongly continuous cosine family {C(t;A),t∈R} and D(A) is a domain of A. It is known that the infinitesimal generator A is closed, densely-defined operator on X.
Theorem 2.1 (See reference [37]). Let {C(t;A):t∈R} be a strongly continuous cosine family in X satisfying ‖C(t;A)‖≤Meω|t| for all t∈R and let A be the infinitesimal generator of {C(t;A):t∈R}. Then for Re(λ)>ω, λ2 is in the resolvent set of A and
λR(λ2;A)x=∞∫0e−λtC(t;A)xdt,x∈X, | (2.9) |
R(λ2;A)x=∞∫0e−λtS(t;A)xdt,x∈X. | (2.10) |
Definition 2.5. See reference [33] Let α∈(1,2). A family Cα(⋅;A):R+→L(X) of all bounded linear operators on X is called a fractional strongly continuous cosine family if it satisfies the following hypotheses:
● Cα(t;A) is strongly continuous for all t∈R+ and Cα(0;A)=I;
● Cα(s;A)Cα(t;A)=Cα(t;A)Cα(s;A) for all s,t∈R+;
● The functional equation
Cα(s;A)IαtCα(t;A)−IαsCα(s;A)Cα(t;A)=IαtCα(t;A)−IαsCα(s;A)holdsforalls,t∈R+. | (2.11) |
The closed linear operator A is defined by
Ax:=Γ(α+1)limt→0+Cα(t;A)x−xtα,x∈D(A):={x∈X:Cα(⋅;A)∈C2(R+,X)}, | (2.12) |
where A is the infinitesimal generator of the fractional strongly continuous cosine family {Cα(t;A):t∈R+}.
Definition 2.6. See reference [33] A fractional strongly continuous cosine family {Cα(t;A):t∈R+} is said to be exponentially bounded if there exists constants M≥1, ω≥0 such that
‖Cα(t;A)‖≤Meωt,t∈R+. | (2.13) |
Definition 2.7. See reference [33] The fractional strongly continuous sine family Sα(⋅;A):R+→L(X) associated with Cα is defined by
Sα(t;A)x=t∫0Cα(s;A)xds,x∈X,t∈R+. | (2.14) |
This implies that
D1tSα(t;A)x=Cα(t;A)x,x∈X,t∈R+. | (2.15) |
Definition 2.8. See reference [33] The fractional strongly continuous Riemann-Liouville family Tα(⋅;A):R+→L(X) associated with Cα is defined by
Tα(t;A)x=Iα−1tCα(t;A)x=t∫0gα−1(t−s)Cα(s;A)xds,x∈X,t∈R+. | (2.16) |
Theorem 2.2 (See reference [33]). Let {Cα(t;A):t∈R+} be a fractional strongly continuous cosine family in X satisfying (2.13) and let A be the infinitesimal generator of {Cα(t;A):t∈R+}. Then for Re(λ)>ω, λα∈ρ(A), α∈(1,2), the following relations hold true:
λα−1R(λα;A)x=∞∫0e−λtCα(t;A)xdt,x∈X, | (2.17) |
λα−2R(λα;A)x=∞∫0e−λtSα(t;A)xdt,x∈X, | (2.18) |
R(λα;A)x=∞∫0e−λtTα(t;A)xdt,x∈X. | (2.19) |
In many concrete situations, the fractional evolution equation is given as a sum of several terms that have various physical meanings and various mathematical properties. While the mathematical analysis for each term may be straightforward, it is not entirely clear what happens after the summation. In the context of perturbed generators of fractional cosine families, we take this as a starting point. To study some perturbation results for fractional strongly continuous cosine families, we first prove the following lemma which establishes a connection between the resolvents of cosine (sine) families generated by A and A+B, respectively.
Note that the following lemma and corollary are true in more general case for closed linear operators and play an important role in the proof of Theorem 3.1.
Lemma 3.1. Let A be a closed linear operator on X to X and suppose B∈L(X) is such that ‖BR(λα;A)‖=θ<1 for some λα∈ρ(A), α∈(1,2). Then, A+B is a closed linear operator with domain D(A) and R(λα;A+B) exists and
R(λα;A+B)=∞∑n=0R(λα;A)[BR(λα;A)]n. | (3.1) |
Furthermore,
‖R(λα;A+B)−R(λα;A)‖≤‖R(λα;A)‖θ(1−θ)−1. | (3.2) |
Proof. Obviously, A+B is a closed linear operator with domain D(A). Since ‖BR(λα;A)‖=θ<1, we note that
R≡∞∑n=0R(λα;A)[BR(λα;A)]n=R(λα;A)[I−BR(λα;A)]−1, |
and hence that
[λαI−(A+B)]R=[I−BR(λα;A)][I−BR(λα;A)]−1=I. |
Moreover, the range of R is precisely D(A) since the range of [I−BR(λα;A)]−1 is X. Therefore, given x∈D(A) there exists a y such that x=Ry. Therefore, we attain that
R[λαI−(A+B)]x=R[λαI−(A+B)]Ry=Ry=x, |
so that R is both a left and a right inverse of [λαI−(A+B)]. The bound ‖R(λα;A+B)−R(λα;A)‖ comes directly from the expansion (3.1):
‖R(λα;A+B)−R(λα;A)‖=‖∞∑n=1R(λα;A)[BR(λα;A)]n‖≤‖R(λα;A)‖∞∑n=0‖BR(λα;A)‖n‖BR(λα;A)‖≤‖R(λα;A)‖(1−‖BR(λα;A)‖)−1‖BR(λα;A)‖=‖R(λα;A)‖θ(1−θ)−1. |
The proof is complete.
Corollary 3.1. Let A be a closed linear operator on X to X and suppose B∈L(X) is such that ‖BR(λα;A)‖=θ<1 for some λα∈ρ(A), α∈(1,2). Then, A+B is a closed linear operator with domain D(A) and λα−1R(λα;A+B), λα−2R(λα;A+B) exist and
λα−1R(λα;A+B)=∞∑n=0λα−1R(λα;A)[BR(λα;A)]n, | (3.3) |
λα−2R(λα;A+B)=∞∑n=0λα−2R(λα;A)[BR(λα;A)]n. | (3.4) |
Furthermore,
‖λα−1R(λα;A+B)−λα−1R(λα;A)‖≤‖λα−1R(λα;A)‖θ(1−θ)−1, | (3.5) |
‖λα−2R(λα;A+B)−λα−2R(λα;A)‖≤‖λα−2R(λα;A)‖θ(1−θ)−1. | (3.6) |
In the following theorem, we impose sufficient conditions such that A is the infinitesimal generator of a fractional strongly continuous cosine (sine) family in X, and B is a bounded linear operator in X, then A+B is also the infinitesimal generator of a fractional strongly continuous cosine (sine) family in X.
Theorem 3.1. Let A be the infinitesimal generator of the fractional strongly continuous cosine family {Cα(t;A):t∈R+} with (2.13), and let {Sα(t;A):t∈R+} denote the fractional strongly continuous sine family associated with Cα where α∈(1,2). Furthermore, suppose that B∈L(X). Then, fractional strongly continuous families of linear bounded operators Cα(⋅;A+B),Sα(⋅;A+B):R+→L(X) generated by A+B (defined on D(A)) can be represented by the series expansion:
Cα(t;A+B)x:=∞∑n=0Cα,n(t;A)x,x∈X, | (3.7) |
Sα(t;A+B)x:=∞∑n=0Sα,n(t;A)x,x∈X, | (3.8) |
where for t∈R+, x∈X and n∈N,
Cα,0(t;A)x:=Cα(t;A)x,Sα,0(t;A)x:=Sα(t;A)x,Cα,n(t;A)x:=t∫0Tα(t−s;A)BCα,n−1(s;A)xds, | (3.9) |
Sα,n(t;A)x:=t∫0Tα(t−s;A)BSα,n−1(s;A)xds. | (3.10) |
Proof. From (2.13), it follows that ‖Sα(t;A)‖≤Mteωt. This implies that Cα(t;A) and Sα(t;A) are strongly continuous on R+. Moreover, by using the formula (2.16), we obtain:
‖Tα(t;A)‖≤t∫0gα−1(t−s)‖Cα(s;A)‖ds≤Meωtt∫0(t−s)α−2Γ(α−1)ds=Meωtgα(t),t∈R+. | (3.11) |
Then, we suppose Cα,n(t;A) and Sα,n(t;A) are strongly continuous on R+ and it is true that for n∈N0:
‖Cα,n(t;A)‖≤Mn+1‖B‖neωtgnα+1(t), | (3.12) |
‖Sα,n(t;A)‖≤Mn+1‖B‖neωtgnα+2(t). | (3.13) |
Firstly, we start with the formula (3.12). Then, Tα(t−s;A)BCα,n(s;A) will be strongly continuous on [0,t] such that the integral defining Cα,n+1(t;A) exists in the strong topology. Moreover, in the case of (3.12), this is true by our remark above for n=0. By mathematical induction principle, we verify [32]:
‖Cα,n+1(t;A)‖≤t∫0‖Tα(t−s;A)‖‖B‖‖Cα,n(s;A)‖ds≤Mn+2‖B‖n+1t∫0gα(t−s)ew(t−s)eωsgnα+1(s)ds=Mn+2‖B‖n+1eωtg(n+1)α+1(t),t∈R+. | (3.14) |
The analogues procedure gives the bound for Sα,n(t;A), t∈R+, as follows:
‖Sα,n+1(t;A)‖≤t∫0‖Tα(t−s;A)‖‖B‖‖Sα,n(s;A)‖ds≤Mn+2‖B‖n+1t∫0gα(t−s)ew(t−s)eωsgnα+2(s)ds=Mn+2‖B‖n+1eωtg(n+1)α+2(t),t∈R+. | (3.15) |
Finally, for t1<t2, we have
‖Cα,n+1(t2;A)x−Cα,n+1(t1;A)x‖≤t1∫0‖[Tα(t2−s)−Tα(t1−s)]BCα,n(s;A)x‖ds+t2∫t1‖Tα(t2−s)‖‖B‖‖Cα,n(s;A)x‖ds,x∈X, | (3.16) |
‖Sα,n+1(t2;A)x−Sα,n+1(t1;A)x‖≤t1∫0‖[Tα(t2−s)−Tα(t1−s)]BSα,n(s;A)x‖ds+t2∫t1‖Tα(t2−s)‖‖B‖‖Sα,n(s;A)x‖ds,x∈X. | (3.17) |
As t1,t2→t0, the integrands in the first terms on the right of (3.16), (3.17) to zero boundedly and the integrands of the second terms of (3.16), (3.17) are bounded. Therefore, Cα,n+1(t;A)x, Sα,n+1(t;A)x are strongly continuous on R+ for each fixed x∈X. By induction, Cα,n(t;A)x, Sα,n(t;A)x are well-defined, strongly continuous, and satisfying (3.12) and (3.13), respectively, for n∈N0. Hence, the fractional-order families Cα(t;A+B), Sα(t;A+B) are uniformly convergent on compact subsets of R+ with respect to the operator norm topology. Moreover, the series are majorized by the following series expansions [32]:
||Cα(t;A+B)||≤MeωtEα(M||B||tα),t∈R+, | (3.18) |
||Sα(t;A+B)||≤MeωttEα,2(M||B||tα),t∈R+. | (3.19) |
Next, we claim that for Cα,n(t;A)x, n∈N for sufficiently large Re(λ) which holds the condition ‖BR(λα;A)‖<1 (it follows from the asymptotic expansion and continuity of the Mittag-Leffler function, see [14], Theorem 2.5, pp. 21):
∞∫0e−λtCα,n(t;A)xdt=λα−1R(λα;A)[BR(λα;A)]nx,x∈X. | (3.20) |
With the help of the mathematical induction principle, well-known Fubini's theorem for iterated integrals, commutative property of convolution operator (2.1) and relation (2.17), for n=1, we verify that
∞∫0e−λtCα,1(t;A)xdt=∞∫0e−λtt∫0Tα(t−s;A)BCα,0(s;A)xdsdt=∞∫0t∫0e−λtTα(t−s;A)BCα(s;A)xdsdt=∞∫0t∫0e−λtTα(s;A)BCα(t−s;A)xdsdt=∞∫0∞∫se−λtTα(s;A)BCα(t−s;A)xdtds=∞∫0∞∫0e−λ(s+t)Tα(s;A)BCα(t;A)xdtds=∞∫0e−λsTα(s;A)B∞∫0e−λtCα(t;A)xdtds=∞∫0e−λtTα(t;A)Bλα−1R(λα;A)xdt=λα−1R(λα;A)BR(λα;A)x,x∈X. |
In a recursive way, for n=2, we obtain that
∞∫0e−λtCα,2(t;A)xdt=∞∫0e−λtt∫0Tα(t−s;A)BCα,1(s;A)xdsdt=∞∫0t∫0e−λtTα(t−s;A)BCα,1(s;A)xdsdt=∞∫0t∫0e−λtTα(s;A)BCα,1(t−s;A)xdsdt=∞∫0∞∫se−λtTα(s;A)BCα,1(t−s;A)xdtds=∞∫0∞∫0e−λ(t+s)Tα(s;A)BCα,1(t;A)xdtds=∞∫0e−λsTα(s;A)B∞∫0e−λtCα,1(t;A)xdtds=∞∫0e−λtTα(t;A)Bλα−1R(λα;A)BR(λα;A)xdt=λα−1R(λα;A)[BR(λα;A)]2x,x∈X. |
Therefore, the relation (3.20) is true for arbitrary n∈N. Then, by the formulae (3.12), (3.18) and Corollary 3.1, we derive that for Cα(t;A+B)x, t∈R+, for chosen sufficiently large Re(λ) which holds the condition ‖BR(λα;A)‖<1,
λα−1R(λα;A+B)x=∞∫0e−λtCα(t;A+B)xdt=∞∑n=0λα−1R(λα;A)[BR(λα;A)]nx,x∈X. | (3.21) |
Therefore, for sufficiently large Re(λ), the Laplace transforms of both Cα(t;A+B) and ∞∑n=0Cα,n(t;A) are equal and hence by the uniqueness theorem, these two functions are equal.
In a similar way, by using the following formula (3.22), Corollary 3.1 and relations (2.18), (3.13), (3.19) for n∈N and sufficiently large Re(λ) which holds the condition ‖BR(λα;A)‖<1:
∞∫0e−λtSα,n(t;A)xdt=λα−2R(λα;A)[BR(λα;A)]nx,x∈X, | (3.22) |
we can derive the corresponding result for fractional strongly continuous sine family:
λα−2R(λα;A+B)x=∞∫0e−λtSα(t;A+B)xdt=∞∑n=0λα−2R(λα;A)[BR(λα;A)]nx,x∈X. | (3.23) |
The proof is complete.
Remark 3.1. Our results coincide with the classical ones whenever α=2 [29]. It should be note that in the case of α=2, Tα(t;A) and Sα(t;A) coincide with the strongly continuous sine function S(t;A) and correspondingly, Cα(t;A) coincide with the strongly continuous cosine function C(t;A). Furthermore, for α=2, two-parameter Mittag-Leffler type functions are converting to the hyperbolic cosine and sine functions, respectively:
E2,1(M||B||t2)=∞∑k=0Mk||B||kt2k(2k)!=cosh(√M||B||t),t∈R,tE2,2(M||B||t2)=∞∑k=0Mk||B||kt2k+1(2k+1)!=1√M||B||sinh(√M||B||t),t∈R. |
Theorem 3.2. Let A be the infinitesimal generator of the strongly continuous cosine family {C(t;A):t∈R} with ‖C(t;A)‖≤Meω|t| for all t∈R, and let {S(t;A):t∈R} denote the strongly continuous sine family associated with C. Furthermore, suppose that B∈L(X). Then, strongly continuous families of linear bounded operators C(⋅;A+B),S(⋅;A+B)∈L(X) generated by A+B (defined on D(A)) can be represented by the series expansion:
C(t;A+B)x:=∞∑n=0Cn(t,A)x,x∈X,S(t;A+B)x:=∞∑n=0Sn(t;A)x,x∈X, |
where for t∈R, x∈X and n∈N,
C0(t;A)x:=C(t;A)x,S0(t;A)x:=S(t;A)x,Cn(t;A)x:=t∫0S(t−s;A)BCn−1(s;A)xds, | (3.24) |
Sn(t;A)x:=t∫0S(t−s;A)BSn−1(s;A)xds. | (3.25) |
Moreover, Travis and Webb have proved anlogues results in [28] for the case of α=2. The authors in [28] have used instead of (3.24) the following recursion formula:
Cn(t;A)x:=t∫0C(t−s;A)BSn−1(s;A)xds,x∈X, | (3.26) |
which can be proven via integrating (formally) by parts formula in (3.24).
We end this chapter by adding some results on fractional analogues of uniformly continuous operator cosine and sine functions which are Mittag-Leffler type operator-valued functions generated by A+B where A,B are linear bounded operators in a Banach space X.
It is known that if Cα(⋅;A):R+→L(X) is fractional uniformly continuous cosine function then there is an A∈L(X) with
Cα(t;A)=Eα(Atα)=∞∑k=0AktkαΓ(kα+1),t∈R+. | (3.27) |
Moreover, the corresponding fractional uniformly continuous sine and Riemann-Liouville families Sα(⋅;A), Tα(⋅;A):R+→L(X) are defined by
Sα(t;A)=tEα,2(Atα)=∞∑k=0Aktkα+1Γ(kα+2),t∈R+, | (3.28) |
and
Tα(t;A)=tα−1Eα,α(Atα)=∞∑k=0Aktkα+α−1Γ(kα+α),t∈R+. | (3.29) |
In this case, we have A=Γ(α+1)limt→0+Cα(t;A)−Itα in the uniform operator topology and the domain D(A) coincides with the state space X i.e., D(A)=X.
For t∈R+, x∈X and n∈N, we define the following sequence of operator-valued functions via a recursive way:
Cα,0(t;A)x:=Cα(t;A)x=Eα(Atα)x,Cα,n(t;A)x:=t∫0(t−s)α−1Eα,α(A(t−s)α)BCα,n−1(s;A)xds, |
and
Sα,0(t;A)x:=Sα(t;A)x=tEα,2(Atα)x,Sα,n(t;A)x:=t∫0(t−s)α−1Eα,α(A(t−s)α)BSα,n−1(s;A)xds. |
Therefore, in this case, fractional uniformly continuous families of linear bounded operators Cα(⋅;A+B),Sα(⋅;A+B):R+→L(X) become the perturbation of a Mittag-Leffler type operator-valued functions defined on X:
Cα(t;A+B)x=Eα((A+B)tα)=∞∑n=0Cα,n(t;A)x,x∈X,Sα(t;A+B)x=tEα,2((A+B)tα)=∞∑n=0Sα,n(t;A)x,x∈X. |
It is interesting to note that depending of commutativity condition of linear bounded operators A,B∈L(X), we can derive the following elegant representation formulas for the perturbation of fractional uniformly continuous cosine and sine families. First, we obtain closed-form representations for perturbation of Mittag-Leffler type functions generated by A+B, where A,B∈L(X) are nonpermutable linear bounded operators in a Banach space X.
Theorem 3.3. Let A,B∈L(X) with non-zero commutator [A,B]:=AB−BA≠0. Then, the perturbation of Mittag-Leffler type operator-valued functions Cα(⋅;A+B),Sα(⋅;A+B):R+→L(X) can be represented by
Cα(t;A+B)x=Eα((A+B)tα)=∞∑k=0∞∑n=0QA,Bk,ntkα+nαΓ(kα+nα+1)x,x∈X, | (3.30) |
and
Sα(t;A+B)x=tEα,2((A+B)tα)x=∞∑k=0∞∑n=0QA,Bk,ntkα+nα+1Γ(kα+nα+2)x,x∈X, | (3.31) |
where QA,Bk,n∈L(X), k,n∈N0 is given by
QA,Bk,0:=Ak,QA,Bk,n:=k∑m=0Ak−mBQA,Bm,n−1,QA,B0,n:=Bn. |
Proof. By making use of the definition of fractional uniformly continuous cosine function Cα(⋅;A):R+→L(X) (3.27), we derive the basis case for n=0:
Cα,0(t;A)=Eα(Atα)=∞∑k=0AktkαΓ(kα+1)=∞∑k=0QA,Bk,0tkαΓ(kα+1),QA,Bk,0:=Ak,k∈N0. |
For n=1, by using the well-known Cauchy product formula and interchanging the order of summation and integration which is permissible in accordance with the uniform convergence of the series (3.27) and (3.29), we get:
Cα,1(t;A)=t∫0(t−s)α−1Eα,α(A(t−s)α)BCα(s;A)ds=∞∑k=0∞∑m=0AkBAmt∫0(t−s)kα+α−1Γ(kα+α)smαΓ(mα+1)ds=∞∑k=0∞∑m=0AkBAmtkα+mα+αΓ(kα+mα+α+1)=∞∑k=0k∑m=0Ak−mBAmtkα+αΓ(kα+α+1)=∞∑k=0k∑m=0Ak−mBQA,Bm,0tkα+αΓ(kα+α+1)=∞∑k=0QA,Bk,1tkα+αΓ(kα+α+1),QA,Bk,1:=k∑m=0Ak−mBQA,Bm,0,k∈N0. |
In a recursive way, for n=2, we derive that
Cα,2(t;A)=t∫0(t−s)α−1Eα,α(A(t−s)α)BCα,1(s;A)ds=∞∑k=0∞∑m=0AkBQA,Bm,1t∫0(t−s)kα+α−1Γ(kα+α)smα+αΓ(mα+α+1)ds=∞∑k=0∞∑m=0AkBQA,Bm,1tkα+mα+2αΓ(kα+mα+1α+2)=∞∑k=0k∑m=0Ak−mBQA,Bm,1tkα+2αΓ(kα+2α+1)=∞∑k=0QA,Bk,2tkα+2αΓ(kα+2α+1),QA,Bk,2:=k∑m=0Ak−mBQA,Bm,1,k∈N0. |
Eventually, for the n-th case, it yields that
Cα,n(t;A)=t∫0(t−s)α−1Eα,α(A(t−s)α)BCα,n−1(s;A)ds=∞∑k=0∞∑m=0AkBQA,Bm,n−1t∫0(t−s)kα+α−1Γ(kα+α)smα+(n−1)αΓ(mα+(n−1)α+1)ds=∞∑k=0∞∑m=0AkBQA,Bm,n−1tkα+mα+nαΓ(kα+mα+nα+1)=∞∑k=0k∑m=0Ak−mBQA,Bm,n−1tkα+nαΓ(kα+nα+1)=∞∑k=0QA,Bk,ntkα+nαΓ(kα+nα+1),QA,Bk,n:=k∑m=0Ak−mBQA,Bm,n−1,k∈N0. |
Therefore, for a perturbation of Mittag-Leffler operator-valued function Cα(⋅;A+B):R+→L(X), with the help of Cauchy product formula, we attain a desired result:
Cα(t;A+B)x=∞∑k=0QA,Bk,0tkαΓ(kα+1)x+∞∑k=0QA,Bk,1tkα+αΓ(kα+α+1)x+∞∑k=0QA,Bk,2tkα+2αΓ(kα+2α+1)x+…+∞∑k=0QA,Bk,ntkα+nαΓ(kα+nα+1)x+…=∞∑k=0∞∑n=0QA,Bk,ntkα+nαΓ(kα+nα+1)x=∞∑k=0k∑n=0QA,Bk−n,ntkαΓ(kα+1)x=∞∑k=0(A+B)ktkαΓ(kα+1)x=Eα((A+B)tα)x,x∈X. |
By making use of the definition of fractional uniformly continuous sine function Sα(⋅;A):R+→L(X) (3.28), we derive the basis case for n=0:
Sα,0(t;A)=tEα,2(Atα)=∞∑k=0Aktkα+1Γ(kα+2)=∞∑k=0QA,Bk,0tkα+1Γ(kα+2),QA,Bk,0:=Ak,k∈N0. |
For n=1, by making use of the above calculations, the well-known Cauchy product formula and interchanging the order of summation and integration which is permissible in accordance with the uniform convergence of the series (3.28) and (3.29), we get:
Sα,1(t;A)=t∫0(t−s)α−1Eα,α(A(t−s)α)BSα(s;A)ds=∞∑k=0∞∑m=0AkBAmt∫0(t−s)kα+α−1Γ(kα+α)smα+1Γ(mα+2)ds=∞∑k=0∞∑m=0AkBAmtkα+mα+α+1Γ(kα+mα+α+2)=∞∑k=0k∑m=0Ak−mBAmtkα+α+1Γ(kα+α+2)=∞∑k=0k∑m=0Ak−mBQA,Bm,0tkα+α+1Γ(kα+α+2)=∞∑k=0QA,Bk,1tkα+α+1Γ(kα+α+2),QA,Bk,1:=k∑m=0Ak−mBQA,Bm,0,k∈N0. |
In a recursive way, for n=2, we attain that
Sα,2(t;A)=t∫0(t−s)α−1Eα,α(A(t−s)α)BSα,1(s;A)ds=∞∑k=0∞∑m=0AkBQA,Bm,1t∫0(t−s)kα+α−1Γ(kα+α)smα+α+1Γ(mα+α+2)ds=∞∑k=0∞∑m=0AkBQA,Bm,1tkα+mα+2α+1Γ(kα+mα+2α+2)=∞∑k=0k∑m=0Ak−mBQA,Bm,1tkα+2α+1Γ(kα+2α+2)=∞∑k=0QA,Bk,2tkα+2α+1Γ(kα+2α+2),QA,Bk,2:=k∑m=0Ak−mBQA,Bm,1,k∈N0. |
Finally, for the n-th case, it follows that
Sα,n(t;A)=t∫0(t−s)α−1Eα,α(A(t−s)α)BSα,n−1(s;A)ds=∞∑k=0∞∑m=0AkBQA,Bm,n−1t∫0(t−s)kα+α−1Γ(kα+α)smα+(n−1)α+1Γ(mα+(n−1)α+2)ds=∞∑k=0∞∑m=0AkBQA,Bm,n−1tkα+mα+nα+1Γ(kα+mα+nα+2)=∞∑k=0k∑m=0Ak−mBQA,Bm,n−1tkα+nα+1Γ(kα+nα+2)=∞∑k=0QA,Bk,ntkα+nα+1Γ(kα+nα+2),QA,Bk,n:=k∑m=0Ak−mBQA,Bm,n−1,k∈N0. |
Therefore, for a perturbation of Mittag-Leffler type operator-valued function Sα(⋅;A+B):R+→L(X), by using the Cauchy product formula, we obtain a desired result:
Sα(t;A+B)x=∞∑k=0QA,Bk,0tkα+1Γ(kα+2)x+∞∑k=0QA,Bk,1tkα+α+1Γ(kα+α+2)x+∞∑k=0QA,Bk,2tkα+2α+1Γ(kα+2α+2)x+…+∞∑k=0QA,Bk,ntkα+nα+1Γ(kα+nα+2)x+…=∞∑k=0∞∑n=0QA,Bk,ntkα+nα+1Γ(kα+nα+2)x=∞∑k=0k∑n=0QA,Bk−n,ntkα+1Γ(kα+2)x=∞∑k=0(A+B)ktkα+1Γ(kα+2)x=tEα,2((A+B)tα)x,x∈X, |
where we have used that the following relation which is a generalization of binomial theorem for nonpermutable linear bounded operators A,B∈L(X)
k∑n=0QA,Bk−n,n=(A+B)k, |
where QA,Bk,n∈L(X), k,n∈N0 is given by
QA,Bk,0:=Ak,QA,Bk,n:=k∑m=0Ak−mBQA,Bm,n−1,QA,B0,n:=Bn. |
The proof is complete.
Remark 3.2. Further, it is interesting to note that the operator (matrix) construction QA,Bk,n, k,n∈N0 is satisfying the following crucial property which is a generalization of the well-known Pascal rule for nonpermutable linear bounded operators (matrices):
QA,Bk,n=AQA,Bk−1,n+BQA,Bk,n−1,k,n∈N0. |
To get more information regarding the properties of QA,Bk,n∈L(X), k,n∈N0, we refer to [38,39].
Secondly, we attain analytical representation formulas for perturbation of Mittag-Leffler type functions generated by A+B, where A,B∈L(X) are permutable linear bounded operators in a Banach space X.
Theorem 3.4. Let A,B∈L(X) with zero commutator [A,B]:=AB−BA=0. Then, the perturbation of Mittag-Leffler type operator-valued functions Cα(⋅;A+B),Sα(⋅;A+B):R+→L(X) can be represented by
Cα(t;A+B)x=Eα((A+B)tα)x=∞∑k=0∞∑n=0(k+nn)AkBntkα+nαΓ(kα+nα+1)x,x∈X, | (3.32) |
and
Sα(t;A+B)x=tEα,2((A+B)tα)x=∞∑k=0∞∑n=0(k+nn)AkBntkα+nα+1Γ(kα+nα+2)x,x∈X. | (3.33) |
Proof. By making use of the definition of fractional uniformly continuous cosine function Cα(⋅;A):R+→L(X) (3.27), we derive the basis case for n=0:
Cα,0(t;A)=Eα(Atα)=∞∑k=0AktkαΓ(kα+1)=∞∑k=0(k+00)AkB0tkαΓ(kα+1). |
For n=1, by using the condition AB=BA, well-known Cauchy product formula and interchanging the order of summation and integration which is permissible in accordance with the uniform convergence of the series (3.28) and (3.29), we get:
Cα,1(t;A)=t∫0(t−s)α−1Eα,α(A(t−s)α)BCα(s;A)ds=∞∑k=0∞∑m=0Ak+mBt∫0(t−s)kα+α−1Γ(kα+α)smαΓ(mα+1)ds=∞∑k=0k∑m=0AkBtkα+αΓ(kα+α+1)=∞∑k=0(k+1)AkBtkα+αΓ(kα+α+1)=∞∑k=0(k+11)AkBtkα+αΓ(kα+α+1). |
In a recursive way, for n=2, we derive that
Cα,2(t;A)=t∫0(t−s)α−1Eα,α(A(t−s)α)BCα,1(s;A)ds=∞∑k=0∞∑m=0AkB(m+11)AmBt∫0(t−s)kα+α−1Γ(kα+α)smα+αΓ(mα+α+1)ds=∞∑k=0∞∑m=0(m+11)Ak+mB2tkα+mα+2αΓ(kα+mα+2α+1)=∞∑k=0k∑m=0(m+11)AkB2tkα+2αΓ(kα+2α+1)=∞∑k=0(k+22)AkB2tkα+2αΓ(kα+2α+1). |
Eventually, for the n-th case, it yields that
Cα,n(t;A)=t∫0(t−s)α−1Eα,α(A(t−s)α)BCα,n−1(s;A)ds=∞∑k=0∞∑m=0AkB(m+n−1n−1)AmBn−1t∫0(t−s)kα+α−1Γ(kα+α)smα+(n−1)αΓ(mα+(n−1)α+1)ds=∞∑k=0∞∑m=0(m+n−1n−1)Ak+mBntkα+mα+nαΓ(kα+mα+nα+1)=∞∑k=0k∑m=0(m+n−1n−1)AkBntkα+nαΓ(kα+nα+1)=∞∑k=0(k+nn)AkBntkα+nαΓ(kα+nα+1). |
Therefore, for a perturbation of Mittag-Leffler operator-valued function Cα(⋅;A+B):R+→L(X), by using well-known Cauchy product formula and binomial theorem, we attain a desired result:
Cα(t;A+B)x=∞∑k=0(k+00)AkB0tkαΓ(kα+1)x+∞∑k=0(k+11)AkBtkα+αΓ(kα+α+1)x+∞∑k=0(k+22)AkB2tkα+2αΓ(kα+2α+1)x+…+∞∑k=0(k+nn)AkBntkα+nαΓ(kα+nα+1)x+…=∞∑k=0∞∑n=0(k+nn)AkBntkα+nαΓ(kα+nα+1)x=∞∑k=0k∑n=0(kn)Ak−nBntkαΓ(kα+1)x=∞∑k=0(A+B)ktkαΓ(kα+1)x=Eα((A+B)tα)x,x∈X. |
By making use of the definition of fractional uniformly continuous sine function Sα(⋅;A):R+→L(X) (3.28), we derive the basis case for n=0:
Sα,0(t;A)=tEα,2(Atα)=∞∑k=0Aktkα+1Γ(kα+2)=∞∑k=0(k+00)AkB0tkα+1Γ(kα+2). |
For n=1, by making use of the above calculations, the well-known Cauchy product formula and interchanging the order of summation and integration which is permissible in accordance with the uniform convergence of the series (3.28) and (3.29), we get:
Sα,1(t;A)=t∫0(t−s)α−1Eα,α(A(t−s)α)BSα(s;A)ds=∞∑k=0∞∑m=0AkB(m+00)AmB0t∫0(t−s)kα+α−1Γ(kα+α)smα+1Γ(mα+2)ds=∞∑k=0∞∑m=0(m+00)Ak+mBtkα+mα+α+1Γ(kα+mα+α+2)=∞∑k=0k∑m=0(m+00)AkBtkα+α+1Γ(kα+α+2)=∞∑k=0(k+11)AkBtkα+α+1Γ(kα+α+2). |
In a recursive way, for n=2, we derive that
Sα,2(t;A)=t∫0(t−s)α−1Eα,α(A(t−s)α)BSα,1(s;A)ds=∞∑k=0∞∑m=0AkB(m+11)AmBt∫0(t−s)kα+α−1Γ(kα+α)smα+α+1Γ(mα+α+2)ds=∞∑k=0∞∑m=0(m+11)Ak+mB2tkα+mα+2α+1Γ(kα+mα+2α+2)=∞∑k=0k∑m=0(m+11)AkB2tkα+2α+1Γ(kα+2α+2)=∞∑k=0(k+22)AkB2tkα+2α+1Γ(kα+2α+2). |
Eventually, for the n-th case, it yields that
Sα,n(t;A)=t∫0(t−s)α−1Eα,α(A(t−s)α)BSα,n−1(s;A)ds=∞∑k=0∞∑m=0AkB(m+n−1n−1)AmBn−1t∫0(t−s)kα+α−1Γ(kα+α)smα+(n−1)α+1Γ(mα+(n−1)α+2)ds=∞∑k=0∞∑m=0(m+n−1n−1)Ak+mBntkα+mα+nα+1Γ(kα+mα+nα+2)=∞∑k=0k∑m=0(m+n−1n−1)AkBntkα+nα+1Γ(kα=nα+2)=∞∑k=0(k+nn)AkBntkα+nα+1Γ(kα+nα+2). |
Therefore, for a perturbation of Mittag-Leffler type operator-valued function Sα(⋅;A+B):R+→L(X), by using the Cauchy product formula for double infinite series, we obtain a desired result:
Sα(t;A+B)x=∞∑k=0(k+00)AkB0tkα+1Γ(kα+2)x+∞∑k=0(k+11)AkBtkα+α+1Γ(kα+α+2)x+∞∑k=0(k+22)AkB2tkα+2α+1Γ(kα+2α+2)x+…+∞∑k=0(k+nn)AkBntkα+nα+1Γ(kα+nα+2)x+…=∞∑k=0∞∑n=0(k+nn)AkBntkα+nα+1Γ(kα+nα+2)x=∞∑k=0k∑n=0(kn)Ak−nBntkα+1Γ(kα+2)x=∞∑k=0(A+B)ktkα+1Γ(kα+2)x=tEα,2((A+B)tα)x,x∈X, |
where we have used the following well-known binomial theorem for permutable linear bounded operators A,B∈L(X)
k∑n=0(kn)Ak−nBn=(A+B)k. |
The proof is complete.
Remark 3.3. This case can also be obtained directly with the help the following property of the construction QA,Bk,n∈L(X), k,n∈N0 which is valid for permutable linear bounded operators (in particular, matrices) A,B∈L(X):
QA,Bk,n=(k+nn)AkBn. |
To get more information regarding this property of QA,Bk,n∈L(X), k,n∈N0, we refer to [38,39].
Let us take X:=L2(0,1) and W2,2(0,1):={ϕ:∃h∈L2(0,1):ϕ(t)=ct+t∗h(t),t∈(0,1)}. Note that c=ϕ′(0) and ϕ′′(t)=h(t) for any t∈(0,1). We consider the following abstract initial-boundary value problem with homogeneous Dirichlet boundary conditions for a one-dimensional perturbed fractional wave equation:
{∂αu(x,t)∂tα=a2∂2u(x,t)∂x2+bu(x,t),x∈(0,1),t>0,u(x,0)=f(x),∂u(x,0)∂t=g(x),x∈(0,1),u(0,t)=u(1,t)=0,t>0, | (4.1) |
where a,b∈R and ∂α∂tα is a fractional partial differentiation operator of order α∈(1,2) in Caputo's sense which is defined by
∂αu(x,t)∂tα=1Γ(2−α)t∫0(t−s)1−α∂2u(x,s)∂s2ds. |
We define the following linear operator A:D(A)⊆X→X as follows:
Au=a2∂2∂x2u,a∈R,u∈D(A), |
with the domain is given by
D(A):={u∈X:u∈W2,2(0,1),u(0)=u(1)=0}. |
Moreover, a linear bounded operator B:X→X is defined by Bu=bu for all b∈R and u∈X. It is obvious that A has eigenvalues of the form λn:=−a2n2π2 and the corresponding eigenvectors are given by un(x):=sin(nπx), x∈(0,1) for any n∈N.
Therefore, if f(x)=∞∑n=1cnsin(nπx) and g(x)=∞∑n=1dnsin(nπx), then by making use the formulas (3.32) and (3.33), the mild solution of abstract Cauchy problem (4.1) can be expressed with the help of perturbation of Mittag-Leffler type functions as below:
u(x,t)=∞∑n=1[Eα((b−a2n2π2)tα)cnsin(nπx)+tEα,2((b−a2n2π2)tα)dnsin(nπx)], | (4.2) |
and it follows that A+B is the infinitesimal generator of strongly continuous fractional cosine Cα(t;A+B) and sine Sα(t;A+B) families, u(x,t)=Cα(t;A+B)f(x)+Sα(t;A+B)g(x).
It is interesting to note that in the case of α=2, the mild solution of the following classical one-dimensional perturbed wave equation:
{∂2∂t2u(x,t)=a2∂2∂x2u(x,t)+bu(x,t),x∈(0,1),t>0,u(x,0)=f(x),ut(x,0)=g(x),x∈(0,1),u(0,t)=u(1,t)=0,t>0, | (4.3) |
can be expressed with the help of perturbed cosine and sine functions as follows:
u(x,t)=∞∑n=1[cos(√a2n2π2−b⋅t)cnsin(nπx)+1√a2n2π2−bsin(√a2n2π2−b⋅t)dnsin(nπx)]. | (4.4) |
Remark 4.1. Furthermore, in the case of b=0, by using a method of seperating variables the explicit representation of a mild solution to perturbed fractional and classical wave Eqs (4.1) and (4.3), respectively can be represented with the Fourier coefficients which is studied in [40]
u(x,t)=∞∑n=1[Eα(−a2n2π2tα)Bn(0)+tEα,2(−a2n2π2tα)B′n(0)]sin(nπx), |
and
u(x,t)=∞∑n=1[cos(anπt)Bn(0)+1anπsin(anπt)B′n(0)]sin(nπx). | (4.5) |
where Bn(0) and B′n(0) for n∈N are Fourier coefficients of f(x) and g(x), respectively, that is
Bn(0)=21∫0f(x)sin(nπx)dx,B′n(0)=21∫0g(x)sin(nπx)dx. |
The goal of this paper is to develop some perturbation results for infinitesimal generators of fractional strongly continuous cosine families. That is, we establish sufficient conditions such that A is the infinitesimal generator of a fractional strongly continuous cosine (sine) family in a Banach space X, and B is a bounded linear operator in X, then A+B is also the infinitesimal generator of a fractional strongly continuous cosine (sine) family in X. Meanwhile, our results coincide with the classical ones whenever α=2. Moreover, depending on commutativity condition of linear bounded operators, we derive the elegant analytical representation formulas for uniformly continuous perturbed fractional operator cosine and sine functions. At last, we present an illustrative example in the context of one-dimensional perturbed fractional wave equation to demonstrate the applicability of our theoretical results and we give some comparisons with the existing literature.
We conclude this paper with some open problems concerning with the perturbation of infinitesimal generators of the fractional strongly continuous cosine families and their applications.
● If A is the infinitesimal generator of a fractional strongly continuous cosine (sine) family and B is a closed (unbounded) linear operator in a Banach space X, is A+B the infinitesimal generator of a fractional strongly continuous cosine (sine) family? This problem has been considered in the classical case for α=2 in [7,28,41], but for α∈(1,2) it still remains an unsolved problem. One way of obtaining these desired results is to apply the principle of subordination, which serves to construct new results in the fractional sense on the basis of the classical results.
● The attained perturbation results can be applied to derive the explicit representation in closed-form for solutions of the following fractional functional evolution equation with α∈(1,2) in a Banach space X based on [42]:
{(Dα0+u)(t)=Au(t)+Bu(t−τ),t>0,τ>0,u(t)=φ(t),−τ≤t≤0, |
where A:D(A)⊆X→X is the infinitesimal generator of a fractional strongly continuous cosine (sine) family and B∈L(X). Moreover, φ:[−τ,0]→X is a continuous function which is describing the prehistory of the system. Note that these ideas are new even for the classical case with α=2.
Therefore, all these questions and more are open for consideration in our future research.
[1] | T. Dauxois, M. Peyrard, Physics of solitons, Cambridge: Cambridge University Press, 2006. |
[2] | Y. Kivshar, G. Agrawal, Optical solitons: from fibers to photonic crystals, London: Academic Press, 2003. |
[3] |
F. M. Mitschke, L. F. Mollenaurer, Experimental observation of interaction forces between solitons in optical fibers, Opt. Lett., 12 (1987), 355–357. doi: 10.1364/OL.12.000355
![]() |
[4] |
A. Antikainen, M. Erkintalo, J. Dudley, G. Genty, On the phase-dependent manifestation of optical rogue waves, Nonlinearity, 25 (2012), R73. doi: 10.1088/0951-7715/25/7/R73
![]() |
[5] |
M. Islam, C. Soccolich, J. Gordon, Ultrafast digital soliton logic gates, Opt. Quant. Electron., 24 (1992), S1215–S1235. doi: 10.1007/BF00624671
![]() |
[6] | J. D. Gibbon, A survey of the origins and physical importance of soliton equations, Philos. T. R. Soc. A, 315 (1985), 335–365. |
[7] | A. Yusuf, Symmetry analysis, invariant subspace and conservation laws of the equation for fluid flow in porous media, Int. J. Geom. Methods M., 17 (2020), 2050173. |
[8] |
T. A. Sulaiman, A. Yusuf, F. Tchier, M. Inc, FMO Tawfiq and F Bousbahi. Lie-Bäcklund symmetries, analytical solutions and conservation laws to the more general (2+1)-dimensional Boussinesq equation, Results Phys., 22 (2021), 103850. doi: 10.1016/j.rinp.2021.103850
![]() |
[9] |
A. I. Aliyu, Y. Li, M. Inc, A. Yusuf, B. Almohsen, Dynamics of solitons to the coupled sine-Gordon equation in nonlinear optics, Int. J. Mod. Phys. B, 35 (2021), 2150043. doi: 10.1142/S0217979221500430
![]() |
[10] |
H. I. Abdel-Gawad, M. Tantawy, M. Inc, A. Yusuf, Construction of rogue waves and conservation laws of the complex coupled Kadomtsev–Petviashvili equation, Int. J. Mod. Phys. B, 34 (2020), 2050115. doi: 10.1142/S0217979220501155
![]() |
[11] |
S. Singh, R. Sakthivel, M. Inc, A. Yusuf, K. Murugesan, Computing wave solutions and conservation laws of conformable time-fractional Gardner and Benjamin–Ono equations, Pramana-J. Phys., 95 (2021), 43. doi: 10.1007/s12043-020-02070-0
![]() |
[12] |
M. Arshad, A. R. Seadawy, D. Lu, Modulation stability and dispersive optical soliton solutions of higher order nonlinear Schrödinger equation and its applications in mono-mode optical fibers, Superlattice. Microst., 113 (2018), 419–429. doi: 10.1016/j.spmi.2017.11.022
![]() |
[13] | D. Lu, C. L. Liu, A sub-ODE method for generalized Gardner and BBM equation with nonlinear terms of any order, Appl. Math. Comput., 217 (2010), 1404–1407. |
[14] |
M. Arshad, A. R. Seadawy, D. Lu, Elliptic function and Solitary Wave Solutions of higherorder nonlinear Schrödinger dynamical equation with fourth-order dispersion and cubicquintic nonlinearity and its stability, Eur. Phys. J. Plus, 132 (2017), 371. doi: 10.1140/epjp/i2017-11655-9
![]() |
[15] |
D. J. Kaup, A. C. Newell, An exact solution for a derivative nonlinear Schrödinger equation, J. Math. Phys., 19 (1978), 798–801. doi: 10.1063/1.523737
![]() |
[16] |
M. Wadati, K. Sogo, Gauge transformations in soliton theory, J. Phys. Soc. Jpn., 52 (1983), 394–398. doi: 10.1143/JPSJ.52.394
![]() |
[17] |
S. Xu, J. He, L. Wang, The Darboux transformation of the derivative nonlinear Schrödinger equation, J. Phys. A: Math. Theor., 44 (2011), 305203. doi: 10.1088/1751-8113/44/30/305203
![]() |
[18] | D. Lu, B. Hong, New exact solutions for the (2 + 1)-dimensional generalized Broer-Kaup system, Appl. Math. Comput., 199 (2008), 572–580. |
[19] | W. Liu, Y. Zhang, J. He, Rogue wave on a periodic background for Kaup-Newell equation, Rom. Rep. Phys., 70 (2018), 106. |
[20] | G. P. Agrawal, Nonlinear fiber optics, 5 Eds., New York, 2013. |
[21] |
W. Maliet, Solitary wave solutions of nonlinear wave equation, Am. J. Phys., 60 (1992), 650–654. doi: 10.1119/1.17120
![]() |
[22] | M. J. Ablowitz, P. A. Clarkson, Solitons, nonlinear evolution equations and inverse scattering, New York: Cambridge University Press, 1991. |
[23] |
R. Hirota, Exact solutions of the Korteweg-de Vries equation for multiple collisions of solitons, Phys. Rev. Lett., 27 (1971), 1192–1194. doi: 10.1103/PhysRevLett.27.1192
![]() |
[24] |
H. Rezazadeh, M. Inc, D. Baleanu, New Solitary Wave Solutions for Variants of (3+1)-Dimensional Wazwaz-Benjamin-Bona-Mahony Equations, Front. Phys., 8 (2020), 332. doi: 10.3389/fphy.2020.00332
![]() |
[25] |
M. Inc, H. Rezazadeh, J. Vahidi, M. Eslami, M. A. Akinlar, M. N. Ali, et al., New solitary wave solutions for the conformable Klein-Gordon equation with quantic nonlinearity, AIMS Mathematics, 5 (2020), 6972–6984. doi: 10.3934/math.2020447
![]() |
[26] | S. Zhang, T. Xia, A generalized F-expansion method with symbolic computation exactly solving Broer-Kaup equations, Appl. Math. Comput., 189 (2007), 836–843. |
[27] |
Q. Zhao, L. H. Wu, Darboux transformation and explicit solutions to the generalized TD equation, Appl. Math. Lett., 67 (2017), 1–6. doi: 10.1016/j.aml.2016.11.012
![]() |
[28] |
D. B. Belobo, T. Dase, Solitary and Jacobi elliptic wave solutions of the generalized Benjamin-Bona-Mahony equation, Commun. Nonlinear Sci., 48 (2017), 270–277. doi: 10.1016/j.cnsns.2017.01.001
![]() |
[29] |
M. Arshad, D. Lu, M. Ur Rehman, I. Ahmed, A. M. Sultan, Optical solitary wave and elliptic function solutions of Fokas-Lenells equation in presence of perturbation terms and its modulation instability, Phys. Scripta, 94 (2019), 105202. doi: 10.1088/1402-4896/ab1791
![]() |
[30] |
W. P. Zhong, R. H. Xie, M. Belic, N. Petrovic, G. Chen, L. Yi, Exact spatial soliton solutions of the two-dimensional generalized nonlinear Schrödinger equation with distributed coefficients, Phys. Rev. A, 78 (2008), 023821. doi: 10.1103/PhysRevA.78.023821
![]() |
[31] |
W. P. Zhong, M. R. Belic, Y. Lu, T. Huang, Traveling and solitary wave solutions to the one-dimensional Gross-Pitaevskii equation, Phys. Rev. E, 81 (2010), 016605. doi: 10.1103/PhysRevE.81.016605
![]() |
[32] |
W. Zhong, W. P. Zhong, M. R. Belic, G. Cai, Embedded solitons in the (2 + 1)-dimensional sine-Gordon equation, Nonlinear Dynam., 100 (2020), 1519–1526. doi: 10.1007/s11071-020-05561-2
![]() |
[33] |
J. Manaan, M. F. Aghdaei, M. Khalilian, R. S. Jeddi, Application of the generalized G-expansion method for nonlinear PDEs to obtaining soliton wave solution, Optik, 135 (2017), 395–406. doi: 10.1016/j.ijleo.2017.01.078
![]() |
[34] |
M. Arshad, D. Lu, J. Wang, (N +1)-dimensional fractional reduced differential transform method for fractional order partial differential equations, Commun. Nonlinear Sci., 48 (2017), 509–519. doi: 10.1016/j.cnsns.2017.01.018
![]() |
[35] |
S. A. Khuri, A. Sayfy, Generalizing the variational iteration method for BVPs: proper setting of the correction functional, Appl. Math. Lett., 68 (2017), 68–75. doi: 10.1016/j.aml.2016.11.018
![]() |
[36] |
N. Savaissou, G. Betchewe, H. Rezazadeh, A. Bekir, S. Y. Doka, Exact optical solitons to the perturbed nonlinear Schrödinger equation with dual-power law of nonlinearity, Opt. Quant. Electron., 52 (2020), 1–16. doi: 10.1007/s11082-019-2116-1
![]() |
[37] |
N. Nasreen, D. Lu, M. Arshad, Optical solitons of nonlinear Schrödinger equation with second order spatiotemporal dispersion and its modulation instability, Optik, 161 (2018), 221–229. doi: 10.1016/j.ijleo.2018.02.043
![]() |
[38] |
Z. Pinar, H. Rezazadeh, M. Eslami, Generalized logistic equation method for Kerr law and dual power law Schrödinger equations, Opt. Quant. Electron., 52 (2020), 1–16. doi: 10.1007/s11082-019-2116-1
![]() |
[39] |
J. Y. Yang, W. X. Ma, Conservation laws of a perturbed Kaup-Newell equation, Mod. Phys. Lett. B, 30 (2016), 1650381. doi: 10.1142/S0217984916503814
![]() |
[40] |
M. A. Akbar, L. Akinyemi, S. W. Shao, A. Jhangeer, H. Rezazadeh, M. M. A. Khater, et al., Soliton solutions to the Boussinesq equation through sine-Gordon method and Kudryashov method, Results Phys., 25 (2021), 104228. doi: 10.1016/j.rinp.2021.104228
![]() |
[41] |
L. Akinyemi, M. Şenol, H. Rezazadeh, H. Ahmad, H. Wang, Abundant optical soliton solutions for an integrable (2+1)-dimensional nonlinear conformable Schrodinger system, Results Phys., 25 (2021), 104177. doi: 10.1016/j.rinp.2021.104177
![]() |
[42] |
J. C. He, Z. Y. Cong, Comment on revision of Kaup-Newell's works on IST for DNLS equation, Commun. Theor. Phys., 50 (2008), 1369–1374. doi: 10.1088/0253-6102/50/6/26
![]() |
[43] |
M. O. Al-Amr, New applications of reduced differential transform method, Alex. Eng. J., 53 (2014), 243–247. doi: 10.1016/j.aej.2014.06.005
![]() |
[44] |
S. Kumar, K. Singh, R. K. Gupta, Coupled Higgs field equation and Hamiltonian amplitude equation: Lie classical approach and (G′/G)-expansion method, Pramana-J. Phys., 79 (2012), 41–60. doi: 10.1007/s12043-012-0284-7
![]() |
[45] |
C. Y. Yang, W. Y. Li, W. T. Yu, M. L. Liu, Y. J. Zhang, G. L. Ma, et al., Amplification, reshaping, fission and annihilation of optical solitons in dispersion decreasing fiber, Nonlinear Dynam., 92 (2018), 203–213. doi: 10.1007/s11071-018-4049-9
![]() |
[46] |
M. Wadati, H, Segur, M. J. Ablowitz, A new Hamiltonian amplitude equation governing modulated wave instabilities, J. Phys. Soc. Jpn., 61 (1992), 1187. doi: 10.1143/JPSJ.61.1187
![]() |
[47] |
X. L. Yang, J. S. Tang, Travelling wave solutions for Konopelchenko-Dubrovsky equation using an extended sinh-Gordon equation expansion method, Commun. Theor. Phys., 50 (2008), 1047. doi: 10.1088/0253-6102/50/5/06
![]() |
[48] |
N. Mahaka, G. Akram, Extension of rational sine-cosine and rational sinh-cosh techniques to extract solutions for the perturbed NLSE with Kerr law nonlinearity, Eur. Phys. J. Plus, 134 (2019), 159. doi: 10.1140/epjp/i2019-12545-x
![]() |
[49] | A. C. Scott, Encyclopedia of nonlinear science, Routledge, Taylor and Francis Group, New York, 2005. |
[50] | P. Rosenau, What is a Compacton?, Notices of the AMS, 52 (2005), 738–739. |
[51] | E. W. Weisstein, Concise Encyclopedia of Mathematics, New York: CRC Press, 2002. |
1. | Arzu Ahmadova, Ismail Huseynov, Nazim I. Mahmudov, Perturbation theory for fractional evolution equations in a Banach space, 2022, 105, 0037-1912, 583, 10.1007/s00233-022-10322-1 | |
2. | Qiang Li, Jina Zhao, Extremal solutions for fractional evolution equations of order 1<γ<2, 2023, 8, 2473-6988, 25487, 10.3934/math.20231301 | |
3. | Javad A. Asadzade, Nazim I. Mahmudov, Delayed analogue of three-parameter pseudo-Mittag-Leffler functions and their applications to Hilfer pseudo-fractional time retarded differential equations, 2024, 65, 0022-2488, 10.1063/5.0159480 |