# An Efficient Method for Time-Fractional Coupled Schrödinger System

Authors: Hossein Aminikhah; A. Refahi Sheikhani; Hadi Rezazadeh
Journal: International Journal of Partial Differential Equations
 (2014)

---

## Abstract

We present a new technique to obtain the solution of time-fractional coupled Schrödinger system. The fractional derivatives are considered in Caputo sense. The proposed scheme is based on Laplace transform and new homotopy perturbation method. To illustrate the power and reliability of the method some examples are provided. The results obtained by the proposed method show that the approach is very efficient and simple and can be applied to other partial differential equations.


---

## Body

## 1. Introduction

The intuitive idea of fractional order calculus is as old as integer order calculus. It can be observed from a letter that was written by Leibniz toĹHôpital. The fractional order calculus is a generalization of the integer order calculus to a real or complex number. Fractional differential equations are used in many branches of sciences, mathematics, physics, chemistry, and engineering. Applications of fractional calculus and fractional-order differential equations include dielectric relaxation phenomena in polymeric materials [1], transport of passive tracers carried by fluid flow in a porous medium in groundwater hydrology [2], transport dynamics in systems governed by anomalous diffusion [3, 4], and long-time memory in financial time series [5] and so on [6, 7]. In particular, recently, much attention has been paid to the distributed-order differential equations and their applications in engineering fields that both integer-order systems and fractional-order systems are special cases of distributed-order systems. The reader may refer to [8–10].Several schemes have been developed for the numerical solution of differential equations. The homotopy perturbation method was proposed by He [11] in 1999. This method has been used by many mathematicians and engineers to solve various functional equations. Homotopy method was further developed and improved by He and applied to nonlinear oscillators with discontinuities [12], nonlinear wave equations [13], and boundary value problems [14]. It can be said that He’s homotopy perturbation method is a universal one and is able to solve various kinds of nonlinear functional equations. For example, it was applied to nonlinear Schrödinger equations [15], to nonlinear equations arising in heat transfer [16], and to other equations [17–20]. In this method, the solution is considered to be an infinite series which usually converges rapidly to exact solutions. In this paper we introduce a new form of homotopy perturbation and Laplace transform methods by extending the idea of [21].We extend the homotopy perturbation and Laplace transform method to solve the time-fractional coupled Schrödinger system. The nonlinear time-fractional coupled Schrödinger partial differential system is as [22]
(1)iDtμu+iux+uxx+u+v+ν1f(|u|2,|v|2)u=0,iDtμv-ivx+vxx+u-v+ν2g(|u|2,|v|2)v=0,
where u,v are unknown functions, ν1,ν2 are real constants, and 0<μ≤1 is a parameter describing the order of the fractional Caputo derivative. f and g are arbitrary (smooth) nonlinear real functions. Nonlinear Schrödinger system is one of the canonical nonlinear equations in physics, arising in various fields such as nonlinear optics, plasma physics, and surface waves [23].This paper is organized as follows. In Section2, we recall some basic definitions and results dealing with the fractional calculus and Laplace transform which are later used in this paper. In Section 3 the homotopy perturbation method is described. The basic idea behind the new method is illustrated in Section 4. Finally, in Section 5, the application of homotopy perturbation and Laplace transform method for solving time-fractional coupled Schrödinger systems are presented.


## 2. Preliminaries and Notations

Some basic definitions and properties of the fractional calculus theory are used in this paper.Definition 1.
A real functionf(t), t>0, is said to be in the space Cμ, μ∈R, if there exists a real number p>μ such that f(t)=tpf1(t), where f1(t)∈C[0,∞) and it is said to be in the space Cμn if and only if f(n)∈Cμ, n∈N. Clearly Cμ∈Cβ if β>μ.Definition 2.
The left-sided Riemann-Liouville fractional integral operator of orderμ>0, of a function f∈Cμ, μ≥-1, is defined as follows:
(2)Jμf(t)={1Γ(μ)∫0t‍(t-τ)μ-1f(τ)dτ,μ>0,f(t),μ=0,
where Γ(·) is the well-known Gamma function.Some of the most important properties of operatorJμ, for f(t)∈Cμ, μ,β≥0, and γ>-1 are as follows:
(3)(1)JμJβf(t)=J(μ+β)f(t),(2)JμJβf(t)=JβJμf(t),(3)Jμtγ=Γ(γ+1)Γ(μ+γ+1)tμ+γ.Definition 3.
Amongst a variety of definitions for fractional order derivatives, Caputo fractional derivative has been used [24, 25] as it is suitable for describing various phenomena, since the initial values of the function and its integer order derivatives have to be specified, so Caputo fractional derivative of function f(t) is defined as
(4)Dμf(t)=Jn-μf(t)=1Γ(n-μ)∫0t‍f(n)(τ)(t-τ)μ-n+1dτ,
where n-1<μ≤n, n∈N, t≥0, and f(t)∈C-1n.In this paper, we have considered time-fractional coupled Schrödinger system, where the unknown functionu(x,t) is assumed to be a causal function of fractional derivatives which are taken in Caputo sense as follows.Definition 4.
The Caputo time-fractional derivative operator of orderμ>0 is defined as
(5)Dtμu(x,t)=Jtn-μu(x,t)=1Γ(n-μ)∫0t‍(t-τ)n-μ-1∂nu(x,τ)∂τndτ,n-1<μ≤n.Definition 5.
The Laplace transform of a functionu(x,t), t>0 is defined as
(6)L[u(x,t)]=∫0∞‍e-stu(x,t)dt,
where s can be either real or complex. The Laplace transform L[u(x,t)] of the Caputo derivative is defined as
(7)L{Dtμu(x,t)}=sμL{u(x,t)}-∑k=0n-1‍u(k)(x,0+)sμ-1-k,n-1<μ≤n.Lemma 6.
Ifn-1<μ≤n, n∈N, and k≥0, then we have
(8)L-1{1sμL{tkμΓ(kμ+1)}}=tμ(k+1)Γ(μ(k+1)+1),
where L-1 is the inverse Laplace transform.The Mittag-Leffler function plays a very important role in the fractional differential equations and in fact it was introduced by Mittag-Leffler in 1903 [26]. The Mittag-Leffler function Eμ(z) with μ>0 is defined by the following series representation:
(9)Eμ(z)=∑n=0∞‍znΓ(μn+1),
where z∈C. For μ=1, (9) becomes
(10)E1(z)=∑n=0∞‍znΓ(n+1)=∑n=0∞‍znn!=ez.
The key result that indicates why Mittag-Leffler functions are so important in fractional calculus is the following theorem. It essentially states that the eigenfunctions of Caputo differential operators may be expressed in terms of Mittag-Leffler functions.Theorem 7 (see [27]).
Ifμ>0 and λ∈R, then we have
(11)Dtμ(Eμ(λzμ))=λEμ(λzμ).


## 3. The Homotopy Perturbation Method

For the convenience of the reader, we will first present a brief account of homotopy perturbation method. Let us consider the following differential equation:(12)A(u)-f(r)=0,r∈Ω,
with boundary conditions
(13)B(u,∂u∂n)=0,r∈Γ,
where A is a general differential operator, B is a boundary operator, f(r) is a known analytic function, and Γ is the boundary of the domain Ω.The operatorA can be generally divided into two parts L and N, where L is linear, while N is nonlinear. Therefore, (12) can be written as follows:
(14)L(u)+N(u)-f(r)=0.
By using homotopy technique, one can construct a homotopy y(r,p):Ω×[0,1]→R which satisfies
(15)H(y,p)=(1-p)[L(y)+L(u0)]+p[A(y)-f(r)]=0,p∈[0,1],r∈Ω,
which is equivalent to
(16)H(y,p)=L(y)-L(u0)+pL(u0)+p[N(y)-f(r)]=0,
where p∈[0,1] is an embedding parameter and u0 is an initial guess approximation of (12) which satisfies the boundary conditions. Clearly, we have
(17)H(y,0)=L(y)-L(u0)=0,H(y,1)=A(y)-f(r)=0.
Thus, the changing process of p from 0 to 1 is just that of y(r,p) from u0(r) to y(r). In topology this is called deformation and L(y)-L(u0) and A(y)-f(r) are called homotopic. If, the embedding parameter p, (0≤p≤1) is considered as a small parameter, applying the classical perturbation technique, we can naturally assume that the solution of (15) and (16) can be given as a power series in p; that is,
(18)y=y0+py1+p2y2+⋯.
According to homotopy perturbation method, the approximation solution of (12) can be expressed as a series of the power of p; that is,
(19)u=lim⁡p→1y=y0+y1+y2+⋯.
The convergence of series (19) has been proved by He in his paper [11]. It is worth noting that the major advantage of homotopy perturbation method is that the perturbation equation can be freely constructed in many ways by homotopy in topology and the initial approximation can also be freely selected. Moreover, the construction of the homotopy for the perturbed problem plays a very important role for obtaining desired accuracy [28].


## 4. Basic Ideas of the Homotopy Perturbation and Laplace Transform Method

To illustrate the basic ideas of this method, we consider the general form of a system of nonlinear fractional partial differential equations:(20)Dtμui(x,t)-Ai(ui(x,t))-fi(x,t)=0,μ∈(0,1],i=1,2,…,n,
with initial conditions
(21)ui(x,0)=ui,i=1,2,…,n,
where Ai are operators and fi(x,t) are known analytical functions. The operators Ai can be divided into two parts, Li and Ni, where Li are the linear operators and Ni are nonlinear operators. Therefore, (20) can be rewritten as
(22)Dtμui(x,t)-Li(ui(x,t))-Ni(ui(x,t))-fi(x,t)=0,i=1,…,n.
By the new homotopy perturbation method [29], we construct the following homotopies:
(23)H(Vi(x,t),p)=(1-p)[DtμVi(x,t)-ui,0(x,t)]+p[DμVi(x,t)-Li(Vi(x,t))-Ni(Vi(x,t))-fi(x,t)]=0,i=1,…,n,
or equivalently
(24)H(Vi(x,t),p)=DtμVi(x,t)-ui,0(x,t)+p[ui,0(x,t)-Li(Vi(x,t))-Ni(Vi(x,t))-fi(x,t)]=0,i=1,…,n,
where p∈[0,1] is an embedding parameter and ui,0(x,t) are initial approximations for the solution of (20). Clearly, we have from (23) and (24)
(25)H(Vi(x,t),0)=DtμVi(x,t)-ui,0(x,t)=0,i=1,…,n,H(Vi(x,t),1)=DtμVi(x,t)-Li(Vi(x,t))-Ni(Vi(x,t))-fi(x,t)=0,i=1,…,n.
By applying Laplace transform on both sides of (24), we have
(26)L{DtμVi(x,t)-ui,0(x,t)+p[ui,0(x,t)-Li(Vi(x,t))-Ni(Vi(x,t))-fi(x,t)]}=0,i=1,…,n.
Using (7), we derive
(27)sμL{Vi(x,t)}-sμ-1Vi(x,0)=L{ui,0(x,t)-p[ui,0(x,t)-Li(Vi(x,t))-Ni(Vi(x,t))-fi(x,t)]},i=1,…,n,
or
(28)L{Vi(x,t)}=1sVi(x,0)+1sμ{L{ui,0(x,t)-p[ui,0(x,t)-Li(Vi(x,t))-Ni(Vi(x,t))-fi(x,t)]}},i=1,…,n.
By applying inverse Laplace transform on both sides of (28), we have
(29)Vi(x,t)=Vi(x,0)+L-1{1sμL{ui,0(x,t)-p[ui,0(x,t)-Li(Vi(x,t))-Ni(Vi(x,t))-fi(x,t)]}1sμ},i=1,…,n.
According to the homotopy perturbation method, we can first use the embedding parameter p as a small parameter and assume that the solution of (29) can be written as a power series in p as follows:
(30)Vi(x,t)=∑j=0∞‍pjVi,j(x,t),i=1,…,n,
where Vi,j(x,t), i=0,…,n, j=1,… are functions which should be determined. Suppose that the initial approximation of the solutions of (20) is in the following form:
(31)ui,0(x,t)=∑j=0∞‍ai,j(x)Γ(jμ+1)tjμ,i=1,…,n,
where ai,j(x), for j=1,2,…, i=1,…,n, are functions which must be computed. Substituting (30) and (31) into (29) and equating terms with identical powers of p we obtain the following set of equations:
(32)p0:Vi,0(x,t)=Vi(x,0)+L-1{1sμL{∑j=0∞‍ai,j(x)Γ(jμ+1)tjμ}},p1:Vi,1(x,t)=L-1{∑n=0∞‍1sμL×{-∑n=0∞‍ai,j(x)Γ(jμ+1)tjμ+[Li(Vi,0(x,t))+Ni(Vi,0(x,t))+fi(x,t)]∑n=0∞‍}},p2:Vi,2(x,t)=L-1{1sμL{Li(Vi,0(x,t),Vi,1(x,t))+Ni(Vi,0(x,t),Vi,1(x,t))}1sμ},p3:Vi,3(x,t)=L-1{1sμL{Li(Vi,0(x,t),Vi1(x,t),Vi,2(x,t))+Ni(Vi,0(x,t),Vi,1(x,t),Vi,2(x,t))}1sμ},⋮pk:Vi,k(x,t)=L-1{1sμL{Li(Vi,0(x,t),Vi,1(x,t),…,Vi,k-1(x,t))+Ni(Vi,0(x,t),Vi,1(x,t),…,Vi,k-1(x,t))}1sμ},⋮Now if we solve these equations in such a way thatVi,1(x,t)=0, then (32) yield
(33)Vi,2(x,t)=Vi,3(x,t)=⋯=0,i=1,…,n.
Therefore the exact solution is obtained by
(34)ui(x,t)=Vi(x,t)=Vi(x,0)+L-1{1sμL{∑j=0∞‍ai,j(x)Γ(jμ+1)tjμ}}=Vi(x,0)+∑j=0∞‍ai,j(x)Γ(μ(j+1))tμ(j+1),i=1,…,n.


## 5. Example

To illustrate the power and reliability of the method for the time-fractional coupled Schrödinger system some examples are provided. The results reveal that the method is very effective and simple.Example 8.
Consider the following linear time-fractional coupled Schrödinger system:(35)iDtμu+iux+uxx+2u+v=f1(x,t),iDtμv-ivx+vxx+u-3v=f2(x,t),
subject to the following initial conditions:
(36)u(x,0)=v(x,0)=0,
where f1(x,t) and f2(x,t) are of the form
(37)f1(x,t)=i22μ+1tμe2iπxΓ(μ+3/2)π(2μ+1)+(-2πe2iπx-4π2e2iπx+3e2iπx)t2μ,f2(x,t)=i22μ+1tμe2iπxΓ(μ+3/2)π(2μ+1)+(2πe2iπx-4π2e2iπx-2e2iπx)t2μ,
with the exact solutions
(38)u(x,t)=v(x,t)=e2iπxt2μ,
where 0<μ≤1.
To solve (35) by the homotopy perturbation and Laplace transform method, we construct the following homotopy:
(39)iDtμU=u0(x,t)-p[u0(x,t)+iUx+Uxx+2U+V-f1(x,t)],iDtμV=v0(x,t)-p[v0(x,t)-iVx+Vxx+U-3V-f2(x,t)].
Applying the Laplace transform on both sides of (39), we have
(40)sμL{U}-sμ-1U(x,0)=-iL{u0(x,0)-p[u0(x,t)+iUx+Uxx+2U+V-f1(x,t)]},sμL{V}-sμ-1V(x,0)=-iL{v0(x,0)-p[v0(x,t)-iVx+Vxx+U-3V-f2(x,t)]},
or
(41)L{U}=1sU(x,0)-isμ{L{u0(x,0)-p[u0(x,t)+iUx+Uxx+2U+V-f1(x,t)]}},L{V}=1sV(x,0)-isμ{L{v0(x,0)-p[v0(x,t)-iVx+Vxx+U-3V-f2(x,t)]}}.
The inverse Laplace transform of (41) and the initial conditions U(x,0)=V(x,0)=0 lead us to
(42)U=-iL-1{1sμ{L{u0(x,0)-p[u0(x,t)+iUx+Uxx+2U+V-f1(x,t)]}}1sμ},V=-iL-1{1sμ{L{v0(x,0)-p[v0(x,t)-iVx+Vxx+U-3V-f2(x,t)]}}1sμ}.
Suppose that the solution is expanded as (30); substituting (30) into (42), collecting the same powers of p, and equating each coefficient of p to zero yield
(43)p0:{U0(x,t)=-iL-1{1sμ{L{u0(x,t)}}},V0(x,t)=-iL-1{1sμ{L{v0(x,t)}}},p1:{U1(x,t)=iL-1{1sμ×{L{u0(x,t)+iU0x+U0xx+2U0+V0-f1(x,t)}}1sμ},V1(x,t)=iL-1{1sμ×{L{v0(x,t)-iV0x+V0xx+U0-3V0-f2(x,t)V0xx}V0xx}1sμ},p2:{U2(x,t)=iL-1{1sμ{L{iU1x+U1xx+2U1+V1}}},V2(x,t)=iL-1{1sμ{L{-iV1x+V1xx+U1-3V1}}},⋮pj:{Uj(x,t)=iL-1{1sμ{L{iUj-1x+Uj-1xx+2Uj-1+Vj-1}}},Vj(x,t)=iL-1{1sμ{L{-iVj-1x+Vj-1xx+Uj-1-3Vj-1}}},⋮
Assumeu0(x,t)=∑n=0∞‍((an(x)/Γ(nμ+1))tnμ), v0(x,t)=∑n=0∞‍((bn(x)/Γ(nμ+1))tnμ). Solving the above equations, for U1(x,t), V1(x,t), leads to the result
(44)U1(x,t)=(ia0(x)Γ(μ+1))tμ+((e2iπxΓ(2μ+1)+d2dx2a0(x)+b0(x)+ia1(x)+iddxa0(x)+2a0(x)d2dx2)×(Γ(2μ+1))-1i(-id2dx2a2(x)+a3(x)+ddxa2(x))t2μ+(i(-id2dx2a1(x)+ddxa1(x)+a2(x)-2ia1(x)-ib1(x)+e2iπxΓ(2μ+1)(-3+4π2+2π)d2dx2)×(Γ(3μ+1))-1i(-id2dx2a2(x)+a3(x)+ddxa2(x))t3μ+(i(-id2dx2a2(x)+a3(x)+ddxa2(x)-2ia2(x)-ib2(x)d2dx2)×(Γ(4μ+1))-1i(-id2dx2a2(x)+a3(x)+ddxa2(x))t4μ+⋯,V1(x,t)=(ib0(x)Γ(μ+1))tμ+((e2iπxΓ(2μ+1)+d2dx2b0(x)+ib1(x)+a0(x)-iddxb0(x)-3b0(x)d2dx2)×(Γ(2μ+1))-1(e2iπxΓ(2μ+1)+d2dx2b0(x)+ib1(x)+a0(x))t2μ+(i(d2dx2-ddxb1(x)+b2(x)-ia1(x)+3ib1(x)-id2dx2b1(x)+2e2iπxΓ(2μ+1)(1+2π2-π))×(Γ(3μ+1))-1i(d2dx2-ddxb1(x)+b2(x)-ia1(x)+3ib1(x))t3μ+(i(d2dx2-ddxb2(x)+b3(x)-ia2(x)+3ib2(x)-id2dx2b2(x))×(Γ(4μ+1))-1)t4μ+⋯.
By the vanishing of U1(x,t), V1(x,t) the coefficients an(x), bn(x)(n=0,1,2,…) are determined to be
(45)a0(x)=0,a1(x)=ie2iπxΓ(2μ+1),a2(x)=a3(x)=⋯=0,b0(x)=0,b1(x)=ie2iπxΓ(2μ+1),b2(x)=b3(x)=⋯=0.
Therefore, the solutions of (35) are
(46)u(x,t)=U0(x,t)=-i(a0(x)tμΓ(μ+1)+a1(x)t2μΓ(2μ+1)+a2(x)t3μΓ(3μ+1)+a3(x)t4μΓ(4μ+1)+⋯)=e2iπxt2μ,v(x,t)=V0(x,t)=-i(b0(x)tμΓ(μ+1)+b1(x)t2μΓ(2μ+1)+b2(x)t3μΓ(3μ+1)+b3(x)t4μΓ(4μ+1)+⋯)=e2iπxt2μ,
which are the exact solutions. Now, if we put μ=1 in (46), we obtain u(x,t)=e2iπxt2, v(x,t)=e2iπxt2 which is the exact solution of the given coupled Schrödinger system (35).Example 9.
Consider the following nonlinear time-fractional coupled Schrödinger system:(47)iDtμu+iux+uxx+u+v+2(|u|2+|v|2)u=f1(x,t),iDtμv-ivx+vxx+u-v+4(|u|2+|v|2)v=f2(x,t),
subject to the following initial conditions:
(48)u(x,0)=v(x,0)=0,
where f1(x,t) and f2(x,t) have the following form:
(49)f1(x,t)=(i22μ+1tμΓ(μ+3/2)π(2μ+1)+4t6μ)eix,f2(x,t)=(i22μ+1tμΓ(μ+3/2)π(2μ+1)+8t6μ)eix,
with the exact solutions
(50)u(x,t)=v(x,t)=eixt2μ,
where 0<μ≤1.
To solve (47) by the homotopy perturbation and Laplace transform method, we construct the following homotopy:
(51)iDtμU=u0(x,t)-p[(|U|2+|V|2)u0(x,t)+iUx+Uxx+U+V+2(|U|2+|V|2)U-f1(x,t)],iDtμV=v0(x,t)-p[(|U|2+|V|2)v0(x,t)-iVx+Vxx+U-V+4(|U|2+|V|2)V-f2(x,t)].
Applying Laplace transform on both sides of (51), we have
(52)sμL{U}-sμ-1U(x,0)=-iL{(|U|2+|V|2)u0(x,0)-p[(|U|2+|V|2)u0(x,t)+iUx+Uxx+U+V+2(|U|2+|V|2)U-f1(x,t)]},sμL{V}-sμ-1V(x,0)=-iL{(|U|2+|V|2)v0(x,0)-p[(|U|2+|V|2)v0(x,t)-iVx+Vxx+U-V+4(|U|2+|V|2)V-f2(x,t)]},
or
(53)L{U}=1sU(x,0)-isμ{(|U|2+|V|2)L{(|U|2+|V|2)u0(x,0)-p[(|U|2+|V|2)u0(x,t)+iUx+Uxx+U+V+2(|U|2+|V|2)U-f1(x,t)]}},L{V}=1sV(x,0)-isμ{(|U|2+|V|2)L{(|U|2+|V|2)v0(x,0)-p[(|U|2+|V|2)v0(x,t)-iVx+Vxx+U-V+4(|U|2+|V|2)V-f2(x,t)]}}.
The inverse Laplace transform of (53) and the initial conditions U(x,0)=V(x,0)=0 lead us to
(54)U=-iL-1{1sμ{(|U|2+|V|2)L{(|U|2+|V|2)u0(x,0)-p[(|U|2+|V|2)u0(x,t)+iUx+Uxx+U+V+2(|U|2+|V|2)U-f1(x,t)]}}1sμ},V=-iL-1{1sμ{(|U|2+|V|2)L{(|U|2+|V|2)v0(x,0)-p[(|U|2+|V|2)v0(x,t)-iVx+Vxx+U-V+4(|U|2+|V|2)V-f2(x,t)]}}1sμ}.
Suppose that the solution is expanded as (30); substituting (30) into (54), collecting the same powers of p, and equating each coefficient of p to zero yield
(55)p0:{U0(x,t)=-iL-1{1sμ{L{u0(x,t)}}},V0(x,t)=-iL-1{1sμ{L{v0(x,t)}}},p1:{U1(x,t)=iL-1{1sμ{L{u0(x,t)+iU0x+U0xx+U0+V0+2(U0U¯0+V0V¯0)U0-f1(x,t)}}1sμ},V1(x,t)=iL-1{1sμ{L{v0(x,t)-iV0x+V0xx+U0-V0+4(U0U¯0+V0V¯0)V0-f2(x,t)}}1sμ},p2:{U2(x,t)=iL-1{1sμ{L{iU1x+U1xx+2U1+V1+2(U1U¯1+V1V¯1)U1}}1sμ},V2(x,t)=iL-1{1sμ{L{-iV1x+V1xx+U1-3V1+4(U1U¯1+V1V¯1)V1}}1sμ},⋮pj:{Uj(x,t)=iL-1{1sμ{L×{iUj-1x+Uj-1xx+2Uj-1+Vj-1+2(Uj-1U¯j-1+Vj-1V¯j-1)Uj-1}}1sμ},Vj(x,t)=iL-1{1sμ{L{-iVj-1x+Vj-1xx+Uj-1-3Vj-1+4(Uj-1U¯j-1+Vj-1V¯j-1)Vj-1}}1sμ},⋮
Assumeu0(x,t)=∑n=0∞‍((an(x)/Γ(nμ+1))tnμ), v0(x,t)=∑n=0∞‍((bn(x)/Γ(nμ+1))tnμ). Solving the above equations, for U1(x,t), V1(x,t), leads to the result
(56)U1(x,t)=(ia0(x)Γ(μ+1))tμ+(((ddxa0(x))eixμ4μΓ(μ)Γ(μ+12)+ia1(x)π+a0(x)π+πd2dx2a0(x)+πb0(x)+i(ddxa0(x))π)×(Γ(2μ+1)π)-1((ddxa0(x))eixμ4μΓ(μ)Γ(μ+12)+ia1(x)π+a0(x)π)t2μ+(i(a2(x)+ddxa1(x)-id2dx2a1(x)-ia1(x)-ib1(x)d2dx2)×(Γ(3μ+1))-1)t3μ+i((a3(x)+ddxa2(x)-id2dx2a2(x)-ia2(x)-ib2(x)d2dx2)×(Γ(4μ+1))-1-iΓ(μ+13)Γ(μ+23)×(a¯0(x)a0(x)+b0(x)b¯0(x))31/2+3μa0(x)×(Γ(4μ+1)(Γ(μ+1))2π)-1(a3(x)+ddxa2(x)-id2dx2a2(x)-ia2(x))t4μ+⋯,V1(x,t)=(ib0(x)Γ(μ+1))tμ+((eixμ4μΓ(μ)Γ(μ+12)+a0(x)π-πb0(x)+πd2dx2b0(x)+ib1(x)π-i(ddxb0(x))π(μ+12))×(Γ(2μ+1)π)-1(eixμ4μΓ(μ)Γ(μ+12)+a0(x)π-πb0(x))t2μ+(i(d2dx2b2(x)-ia1(x)+ib1(x)-ddxb1(x)-id2dx2b1(x))×(Γ(3μ+1))-1)t3μ+i((-ddxb2(x)-ia2(x)-id2dx2b2(x)+ib2(x)+b3(x)d2dx2)×(Γ(4μ+1))-1-2iΓ(μ+13)Γ(μ+23)×(a0(x)a¯0(x)+b0(x)b¯0(x))31/2+3μb0(x)×(Γ(4μ+1)(Γ(μ+1))2π)-1(-ddxb2(x)-ia2(x)-id2dx2b2(x)+ib2(x))t4μ+⋯.
By the vanishing of U1(x,t), V1(x,t) the coefficients an(x), bn(x)(n=0,1,2,…) are determined to be
(57)a0(x)=0,a1(x)=ieixΓ(2μ+1),a2(x)=a3(x)=⋯=0,b0(x)=0,b1(x)=ieixΓ(2μ+1),b2(x)=b3(x)=⋯=0.
Therefore, the exact solutions of the system of (47) can be expressed as
(58)u(x,t)=U0(x,t)=-i(a0(x)tμΓ(μ+1)+a1(x)t2μΓ(2μ+1)+a2(x)t3μΓ(3μ+1)+a3(x)t4μΓ(4μ+1)+⋯)=eixt2μ,v(x,t)=V0(x,t)=-i(b0(x)tμΓ(μ+1)+b1(x)t2μΓ(2μ+1)+b2(x)t3μΓ(3μ+1)+b3(x)t4μΓ(4μ+1)+⋯)=eixt2μ.
Now, if we put μ=1 in (58), we obtain u(x,t)=eixt2, v(x,t)=eixt2 which is the exact solution of the given coupled Schrödinger system (47).Example 10.
Consider the following nonlinear time-fractional coupled Schrödinger system:(59)iDtμu+iux+uxx+u+v+2(|u|2+|v|2)u=f1(x,t),iDtμv-ivx+vxx+u-v-3(|u|2-|v|2)v=f2(x,t),
subject to the following initial conditions:
(60)u(x,0)=eix,v(x,0)=e-ix,
where f1(x,t) and f2(x,t) have the following form:
(61)f1(x,t)=(ieix+4eixEμ2(tμ)-2isinx)Eμ(tμ),f2(x,t)=(ie-ix-2e-ix+2isinx)Eμ(tμ),
with the exact solutions
(62)u(x,t)=eixEμ(tμ),v(x,t)=e-ixEμ(tμ),
where 0<μ≤1.
To solve (59) by the LTNHPM, we construct the following homotopy:
(63)iDtμU=u0(x,t)-p[(|U|2+|V|2)u0(x,t)+iUx+Uxx+U+V+2(|U|2+|V|2)U-f1(x,t)],iDtμV=v0(x,t)-p[(|U|2+|V|2)v0(x,t)-iVx+Vxx+U-V-3(|U|2-|V|2)V-f2(x,t)].
Applying Laplace transform on both sides of (63), we have
(64)sμL{U}-sμ-1U(x,0)=-iL{(|U|2+|V|2)u0(x,0)-p[(|U|2+|V|2)u0(x,t)+iUx+Uxx+U+V+2(|U|2+|V|2)U-f1(x,t)]},sμL{V}-sμ-1V(x,0)=-iL{(|U|2+|V|2)v0(x,0)-p[(|U|2+|V|2)v0(x,t)-iVx+Vxx+U-V-3(|U|2-|V|2)V-f2(x,t)]},
or
(65)L{U}=1sU(x,0)-isμ{(|U|2+|V|2)L{(|U|2+|V|2)u0(x,0)-p[(|U|2+|V|2)u0(x,t)+iUx+Uxx+U+V+2(|U|2+|V|2)U-f1(x,t)]}},L{V}=1sV(x,0)-isμ{(|U|2+|V|2)L{(|U|2+|V|2)v0(x,0)-p[(|U|2+|V|2)v0(x,t)-iVx+Vxx+U-V-3(|U|2-|V|2)V-f2(x,t)]}}.
The inverse Laplace transform (65) and the initial conditions U(x,0)=eix, V(x,0)=e-ix, lead us to
(66)U=eix-iL-1{1sμ{(|U|2+|V|2)L{(|U|2+|V|2)u0(x,0)-p[(|U|2+|V|2)u0(x,t)+iUx+Uxx+U+V+2(|U|2+|V|2)U-f1(x,t)]}}1sμ},V=e-ix-iL-1{1sμ{(|U|2+|V|2)L{(|U|2+|V|2)v0(x,0)-p[(|U|2+|V|2)v0(x,t)-iVx+Vxx+U-V-3(|U|2-|V|2)V-f2(x,t)]}}1sμ}.
Suppose that the solution is expanded as (30); substituting (30) into (66), collecting the same powers of p, and equating each coefficient of p to zero yield
(67)p0:{U0(x,t)=eix+L-1{-isμ{L{u0(x,t)}}},V0(x,t)=e-ix+L-1{-isμ{L{v0(x,t)}}},p1:{U1(x,t)=iL-1{1sμ{L{u0(x,t)+iU0x+U0xx+U0+V0+2(U0U¯0+V0V¯0)U0-f1(x,t)}}1sμ},V1(x,t)=iL-1{1sμ{L{v0(x,t)-iV0x+V0xx+U0-V0-3(U0U¯0-V0V¯0)V0-f2(x,t)}}1sμ},p2:{U2(x,t)=iL-1{1sμ{L{iU1x+U1xx+2U1+V1+2(U1U¯1+V1V¯1)U1}}1sμ},V2(x,t)=iL-1{1sμ{L{-iV1x+V1xx+U1-3V1-3(U1U¯1-V1V¯1)V1}}1sμ},⋮pj:{Uj(x,t)=iL-1{1sμ{L{iUj-1x+Uj-1xx+2Uj-1+Vj-1+2(Uj-1U¯j-1+Vj-1V¯j-1)Uj-1}}1sμ},Vj(x,t)=iL-1{1sμ{L{-iVj-1x+Vj-1xx+Uj-1-3Vj-1-3(Uj-1U¯j-1-Vj-1V¯j-1)Vj-1}}1sμ},⋮
Assumeu0(x,t)=∑n=0∞‍((an(x)/Γ(nμ+1))tnμ), v0(x,t)=∑n=0∞‍((bn(x)/Γ(nμ+1))tnμ). Solving the above equations, for U1(x,t), V1(x,t), we obtain the result
(68)U1(x,t)=(i(a0(x)-ieix)Γ(μ+1))tμ+(i(a1(x)-7ia0(x)+2ib¯0(x)-ib0(x)-id2dx2a0(x)+ddxa0(x)-e-ix-(11+i)eix+2ie2ix(-b0(x)+a¯0(x))d2dx2)×(Γ(2μ+1))-1d2dx2)t2μ+⋯,V1(x,t)=(i(b0(x)-ie-ix)Γ(μ+1))tμ+(i(b1(x)-2ib0(x)-ddxb0(x)-3ia¯0(x)-ia0(x)-id2dx2b0(x)-eix+3ie-2ix(b¯0(x)+a0(x))+(3-i)e-ixddx)×(Γ(2μ+1))-1ddx)t2μ+⋯.
By the vanishing of U1(x,t), V1(x,t) the coefficients an(x), bn(x)(n=0,1,2,…) are determined to be
(69)a0(x)=a1(x)=a2(x)=a3(x)=⋯=ieix,b0(x)=b1(x)=b2(x)=b3(x)=⋯=ie-ix.
Therefore, the solutions of (59) are
(70)u(x,t)=U0(x,t)=eix-i(a0(x)tμΓ(μ+1)+a1(x)t2μΓ(2μ+1)+a2(x)t3μΓ(3μ+1)+a3(x)t4μΓ(4μ+1)+⋯)=eixEμ(tμ),v(x,t)=V0(x,t)=e-ix-i(b0(x)tμΓ(μ+1)+b1(x)t2μΓ(2μ+1)+b2(x)t3μΓ(3μ+1)+b3(x)t4μΓ(4μ+1)+⋯)=e-ixEμ(tμ).
Now, if we put μ=1 in (70), we obtain u(x,t)=eix+t, v(x,t)=e-ix+t which is the exact solution of the given coupled Schrödinger system (59).


## 6. Conclusion

In this paper, we have introduced a combination of Laplace transform and homotopy perturbation methods for solving fractional Schrödinger equations which we called homotopy perturbation and Laplace transform method. In this scheme, the solution considered to be a Taylor series which converges rapidly to the exact solution of the nonlinear equation. As shown in the three examples of this paper, a clear conclusion can be drawn from the results that the homotopy perturbation and Laplace transform method provide an efficient method to handle nonlinear partial differential equations of fractional order. The computations associated with the examples were performed using Maple 13.



---

*Source: 137470-2014-07-15.xml*