TSWJ The Scientific World Journal 1537-744X 2356-6140 Hindawi Publishing Corporation 10.1155/2014/964643 964643 Research Article Modified Fractional Variational Iteration Method for Solving the Generalized Time-Space Fractional Schrödinger Equation http://orcid.org/0000-0002-4764-3594 Hong Baojian 1, 2 Lu Dianchen 1 Inc Mustafa 1 Faculty of Science Jiangsu University Zhenjiang, Jiangsu 212013 China ujs.edu.cn 2 Department of Basic Causes Nanjing Institute of Technology Nanjing 211167 China jit.edu.cn 2014 392014 2014 08 05 2014 24 07 2014 27 07 2014 4 9 2014 2014 Copyright © 2014 Baojian Hong and Dianchen Lu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Based on He’s variational iteration method idea, we modified the fractional variational iteration method and applied it to construct some approximate solutions of the generalized time-space fractional Schrödinger equation (GFNLS). The fractional derivatives are described in the sense of Caputo. With the help of symbolic computation, some approximate solutions and their iterative structure of the GFNLS are investigated. Furthermore, the approximate iterative series and numerical results show that the modified fractional variational iteration method is powerful, reliable, and effective when compared with some classic traditional methods such as homotopy analysis method, homotopy perturbation method, adomian decomposition method, and variational iteration method in searching for approximate solutions of the Schrödinger equations.

1. Introduction

In the past decades, due to the numerous applications of fractional differential equations (FDEs) in the areas of nonlinear science , many important phenomena can be described successfully using the FDEs models such as materials and processes , engineering and physics , dielectric polarization , and quantitative finance . Searching for solutions of these FDEs plays an important and significant role in all aspects of this subject. But because of the complexity of nonlinear terms and fractional derivative, it is very difficult for us to obtain the exact analytic solutions of most FDEs, so approximate and numerical methods must be considered. A great deal of efforts have been proposed for these problems, like the homotopy analysis method (HAM) , the homotopy perturbation method (HPM) , the adomian decomposition method (ADM) , the generalized differential transform method , and so forth .

The variational iteration method (VIM) established in 1999 by He in  is thoroughly used by many researchers to construct the approximate solutions of a wide variety of scientific and engineering models [12, 13]. After some modifications, the fractional variational iteration method (FVIM) was applied to fractional differential equations by He and many authors . The motivation of this paper is to construct some analytical approximate solutions for the GFNLS powerfully. Firstly, we give some modifications for the FVIM and extend the application of the FVIM. Secondly, we use the modified fractional variational iteration method (MFVIM) to the GFNLS and compare the efficiency of MFVIM with some other traditional perturbation methods. The results show that MFVIM gives rapid and standard convergence to the exact solution if such a solution exists.

We give some basic definitions and properties of the fractional calculus theory which are used further in this paper; we define the following fractional integral and derivatives [20, 21].

Definition 1.

A real function f(x) is said to be in the space Cμ, where μR, x>0, if there exists a real number p(>μ) such that f(x)=xpf1(x), where f1(x)C[0,) and it is said to be in the space Cμm if and only if f(m)Cμ, mN.

Definition 2.

The Riemann-Liouville fractional integral operator of order α>0 for a function f(x)Cμ, μ-1, is defined as follows: (1)Jxαf(x)=1Γ(α)0x(x-ξ)α-1f(ξ)dξ=1Γ(α+1)0xf(ξ)(dξ)α,α>0,x>0,Jx0f(x)=f(x). Also one has the following properties: (2)JαJβf(x)=Jα+βf(x),JαJβf(x)=JβJαf(x),Jαxγ=Γ(γ+1)Γ(α+γ+1)xα+γ.

Definition 3.

For α>0, x>0, and f(x)C-1n, the Caputo fractional derivative operator of order α on the whole space is defined as follows: (3)Dαf(x)=Jn-αDnf(x)={1Γ(n-α)0x(x-ξ)n-α-1f(n)(ξ)dξ,FFFn-1<α<n,nN,d(n)f(x)dxn,α=n.

Also one has the following properties: (4)DαC=0,(Cis  a  constant),Dαxγ={Γ(γ+1)Γ(γ-α+1)xγ-α,γ>α-1,0,γα-1,JαDαf(x)=f(x)-k=0n-1f(k)(0+)xkk!,n-1<α<n,DαJαf(x)=f(x).

2. Analysis of the MFVIM and the FGNLS

Consider the following generalized time and space fractional nonlinear Schrödinger equation with variable coefficients [22, 23]: (5)iαutα+a2βux2β+v(x)u+γ|u|2u=0,t>0,0<α,β1,u(x,0)=f(x), where u=u(x,t), αu/tα=Dtαu, 2βu/x2β=Dxβ(Dxβu), v(x) is the trapping potential, and a, γ are the slowly increasing dispersion coefficient and nonlinear coefficient, respectively. If we select α=β=1, v(x)=0, this equation turns to the famous nonlinear Schrödinger equations in optical fiber .

According to the FVIM , we can build a correction functional for (5) as follows: (6)un+1=un+1Γ(1+α)×0tλ(τ,x)[iαuntα+a2βu~nx2βFF,+v(x)u~n+γu~n|u~n|22βu~nx2β](dτ)α, with the initial condition u0=u(x,0)=f(x), where λ(t,x) is a general Lagrange’s multiplier which can be identified optimally with the variational theory. The function u~n is a restricted variation which means δu~n=0. Therefore, we first determine the Lagrange multiplier λ(t,x) that will be identified optimally via integration by parts . The successive approximations un+1, n0, of the solution u(x,t) will be readily obtained through λ(t,x) and any selective function u0. The initial values are usually used for choosing the zeroth approximation u0. With λ(t,x) determined, then several approximations uk, k=1,2, follow immediately. Consequently, the exact solution may be procured by using u=limnun. The convergence of FVIM has been proved in . In this paper, notice that (5) is a complex differential equation with complex modulus term |u|2, as we all know, a complex function u(ξ) can be written as c(ξ)eiθ(ξ), where c(ξ) and θ(ξ) are real functions, noticed that |u(ξ)|2=|c(ξ)|2, we can give some modification for the iteration formulation (6), assume that limn|u~n|2=|u|2=|u0|2, we get the MFVIM for (5). This modification should enhance rapidly the efficiency of our iteration.

In what follows, in order to illustrate the strength of this method, we will apply the MFVIM to some models about (5).

3. Approximate Solutions for the FGNLS Example 4.

We first consider the time-fractional NLS equation : (7)iαutα+a2ux2+γu|u|2=0,t>0,0<α1,u(x,0)=f(x). The correction functional for (7) reads (8)un+1=un+1Γ(1+α)×0tλ(τ,x)[iαuntα+a2u~nx2+γu~n|u~n|2](dτ)α. Making the above correction functional stationary, (9)δun+1=δun+λ(t,x)iδun-1Γ(1+α)0t[iDτ(α)λ(τ,x)δun](dτ)α. After getting the coefficients of δun to zero we can determine the Lagrange multiplier (10)λ=i. We produce the iteration formulation as follows: (11)un+1=un+iΓ(1+α)0t[iαunτα+a2unx2+γun|u0|2](dτ)α. As stated before, we can select u0=u(x,0)=Asechx; using the iteration (11) and the mathematica software, we obtain the following successive approximations: (12)u0=Asechx,u1=Asechx+aAi(sechx-2sech3x)1Γ(1+α)tα+A3iγsech3x1Γ(1+α)tα=Asechx[1+aitαΓ(1+α)]+(-2a+A2γ)Aisech3x1Γ(1+α)tα=Asechx[1+aitαΓ(1+α)],(A2=2aγ),u2=Asechx+aAi(sechx-2sech3x)×[1Γ(1+α)tα+ai1Γ(1+2α)t2α]+γA3isech3x[1Γ(1+α)tα+ai1Γ(1+2α)t2α]=Asechx[1+aitαΓ(1+α)+1Γ(1+2α)(aitα)2],u3=Asechx[1+aitαΓ(1+α)+(aitα)2Γ(1+2α)+(aitα)3Γ(1+3α)],un=Asechxk=0n1Γ(1+kα)(aitα)k. The exact solution of (7) is (13)u=Asechxlimnk=0n1Γ(1+kα)(aitα)k=±2aγsechxEα(aitα), where Eα(aitα) is the Mittag-Leffler function. If we let α=1 in (13), the exact solution of the regular NLS equation (7) can be obtained as follows: (14)u|α=1=±2aγsechxeiat.

Example 5.

We now consider the time-space fractional NLS equation [30, 31]: (15)iαutα+a2βux2β+2au|u|2=0,t>0,0<α,β1,u(x,0)=eix. With the similar process, we get the iteration formulation as follows: (16)un+1=un+iΓ(1+α)×0t[iαunτα+a2βunx2β+2aun|u0|2](dτ)α. Using the iteration (16) and the mathematica software, we obtain the following successive approximations: (17)u0=eix,u1=eix[1+(2+eiπβ)iatαΓ(1+α)]=eix(1+c1iatαΓ(1+α)),c1=2+eiπβ,u2=eix{1+(2+eiπβ)iatαΓ(1+α)+[(2+eiπβ)eiπβ+2(2+eiπβ)]i2a2t2αΓ(1+2α)(2+eiπβ)iatαΓ(1+α)}=eix[1+c1iatαΓ(1+α)+c12i2a2t2αΓ(1+2α)],u3=eix{1+(2+eiπβ)iatαΓ(1+α)+(2c1+c1eiπβ)i2a2t2αΓ(1+2α)+(2c12+c12eiπβ)i3a3t3αΓ(1+3α)(2+eiπβ)iatαΓ(1+α)}=eix[1+c1iatαΓ(1+α)+c12i2a2t2αΓ(1+2α)+c13i3a3t3αΓ(1+3α)],un=eixk=0n1Γ(1+kα)[i(2+eiπβ)atα]k. The exact solution of (15) is (18)u=eixlimnk=0n1Γ(1+kα)[i(2+eiπβ)atα]k=eixEα[i(2+eiπβ)atα]. If we let α=1 and let β=1 in (18), the exact solution of the regular NLS equation (15) can be obtained as follows: (19)u|(α=1,β=1)=ei(x+at).

Remark 6.

The solution (18) is more standard than the result (3.18) in . If one selects a=1 or a=1/2, the solution (19) is the same as the result (49) in , the result (3.21) in , and the result (29) in , but one can find that this iteration is much more standard and powerful than the HAM, the ADM, and the VIM mentioned in .

Example 7.

Consider the following time-space fractional NLS equation : (20)iαutα+122βux2β-ucos2x-|u|2u=0,t>0,0<α1,u(x,0)=sinx. With the similar process, we get the iteration formulation as follows: (21)un+1=un+iΓ(1+α)×0t[iαunτα+122βunx2β-uncos2x-un|u0|2](dτ)α. If we select u0=u(x,0)=sinx, using the iteration (21) and the mathematica software, we obtain the following successive approximations: (22)u0=sinx,u1=sinx+[12sin(x+πβ)-sinx]itαΓ(1+α),u2=sinx+[12sin(x+πβ)-sinx]itαΓ(1+α)+[14sin(x+2πβ)-sin(x+πβ)+sinx]×i2t2αΓ(1+2α),u3=u2+[18sin(x+3πβ)-34sin(x+2πβ)+32sin(x+πβ)-sinx]i3t3αΓ(1+3α),un=k=0nck(x)(itα)kΓ(1+kα), where ck(x)=ck,0sinx+ck,1sin(x+πβ)+ck,2sin(x+2πβ)++ck,ksin(x+kπβ), ck,0=(-1)k, ck,1=(1/2)ck-1,1-ck-1,0,,ck,k-1=(1/2)ck-1,k-2-ck-1,k-1, and ck,k=(1/2)ck-1,k-1, k2, c0,0=1; c1,0=-1, c1,1=(1/2); c2,0=1, c2,1=-1, c2,2=(1/4); .

The exact solution of (20) is (23)u=limnk=0nck(x)(itα)kΓ(1+kα)=sin[xβΓ(1+β)]Exp[-32itαΓ(1+α)]. If we let α=1 and let β=1 in (23), the exact solution of the regular NLS equation (20) can be obtained as follows. (24)u|(α=1,β=1)=sinxe-(3/2)it. If we select u0=u(x,0)=cosx, with the same process, we can also obtain the following exact solution of (20): (25)u=icos[xβΓ(1+β)]Exp[-12itαΓ(1+α)].

Remark 8.

If one selects β=1, the solution (23) is more standard than the result (5.10) in . The solutions (23) and (25) are new exact solutions for (20) to our knowledge.

Comparisons between the real part of some numerical results and the exact solution (23) are summarized in Tables 1 and 2, and the simulations for u4, uabs, and u are plotted in Figures 1 and 2, which shows that the MFVIM produced a rapidly convergent series.

Comparison between the real part of u4 and u when α=β=1.

x t Approximate solution u4appr Exact solution Absolute error
1 0.4 0.6945501509 0.6944959727 0.0000541782
5 0.4 −0.7914960963 −0.7914343559 0.0000617404
1 0.3 0.7577097797 0.7577001100 9.66968 × 10 - 6
15 0.3 0.5855572741 0.5855498014 7.47272 × 10 - 6
12 0.2 −0.5126082301 −0.5126076876 6.16597 × 10 - 6
3 0.2 0.13481723570 0.13481709305 1.42655 × 10 - 7
2 0.1 0.8990870113 0.8990869969 1.43796 × 10 - 8
0.2 0.1 0.1964384915 0.1964384884 3.14175 × 10 - 9

Comparison between the real part of u4 and u when α=0.7,β=0.9.

x t Approximate solution u4appr Exact solution Absolute error
1 0.4 0.5356787528 0.5565185584 0.02083980558
2 0.4 0.5092408025 0.6018550029 0.0926142003
1 0.3 0.6273816439 0.6535277868 0.02614614289
2 0.3 0.6242882934 0.7067670289 0.08247873543
1 0.2 0.7153235950 0.7417557951 0.02643220012
0.2 0.2 0.1978265247 0.2080382267 0.01021170194
0.2 0.1 0.1989280524 0.2288404399 0.02991238755
12 0.1 0.8741657132 0.8902825407 0.01611682753

Comparison between the real part of u4 and the exact solution u.

Plots of the absolute error uabs when α=0.7 and β=0.9.

4. Summary

In this paper, the MFVIM is used for finding approximate and exact solutions of the GFNLS equation with Caputo derivative. The obtained results indicate that the MFVIM is effective, convenient, and powerful method for solving nonlinear fractional complex differential equations when comparing it with some other traditional asymptotic decomposition methods such as HAM, VIM, and ADM. We believe that these methods should play an important role for finding exact and approximate solutions in the mathematical physics.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This paper is supported by the National Nature Science Foundation of China (no. 61070231); the Outstanding Personal Program in Six Fields of Jiangsu (2009188); and the General Program of Innovation Foundation of NanJing Institute of Technology (Grant no. CKJB201218).

Dalir M. Bashour M. Applications of fractional calculus Applied Mathematical Sciences 2010 4 1021 1032 MR2607355 2-s2.0-77953332285 Podlubny I. Fractional Differential Equations 1999 New York, NY, USA Academic Press MR1658022 Miller K. S. Ross B. An Introduction to the Fractional Calculus and Fractional Differential Equations 1993 New York, NY, USA John Wiley & Sons Sun H. H. Abdelwahab A. A. Onaral B. Linear approximation of transfer function with a pole of fractional order IEEE Transactions on Automatic Control 1984 29 5 441 444 2-s2.0-0021424253 Laskin N. Fractional market dynamics Physica A 2000 287 3-4 482 492 10.1016/S0378-4371(00)00387-3 MR1802416 2-s2.0-0034517069 Panda S. Bhowmik A. Das R. Repaka R. Martha S. C. Application of homotopy analysis method and inverse solution of a rectangular wet fin Energy Conversion and Management 2014 80 305 318 Gepreel K. A. The homotopy perturbation method applied to the nonlinear fractional Kolmogorov-Petrovskii-Piskunov equations Applied Mathematics Letters 2011 24 8 1428 1434 10.1016/j.aml.2011.03.025 MR2793646 2-s2.0-79955470464 Song L. N. Wang W. G. A new improved Adomian decomposition method and its application to fractional differential equations Applied Mathematical Modelling 2013 37 3 1590 1598 10.1016/j.apm.2012.03.016 MR3002246 2-s2.0-84870243484 Odibat Z. Momani S. Erturk V. S. Generalized differential transform method: application to differential equations of fractional order Applied Mathematics and Computation 2008 197 2 467 477 10.1016/j.amc.2007.07.068 MR2400669 ZBL1141.65092 2-s2.0-39449122795 Inc M. Akgül A. Numerical solution of seventh-order boundary value problems by a novel method Abstract and Applied Analysis 2014 2014 9 745287 10.1155/2014/745287 MR3191062 He J. H. Variational iteration method—a kind of non-linear analytical technique: some examples International Journal of Non-Linear Mechanics 1999 34 4 699 708 2-s2.0-0000092673 Noor M. A. Mohyud-Din S. T. Variational iteration method for solving higher-order nonlinear boundary value problems using He's polynomials International Journal of Nonlinear Sciences and Numerical Simulation 2008 9 141 156 Momani S. Abuasad S. Application of He's variational iteration method to Helmholtz equation Chaos, Solitons & Fractals 2006 27 5 1119 1123 10.1016/j.chaos.2005.04.113 MR2164845 2-s2.0-24944474278 He J. H. Some applications of nonlinear fractional differential equations and their approximations Bulletin of Science, Technology & Society 1999 15 2 86 90 He J. A short remark on fractional variational iteration method Physics Letters A 2011 375 38 3362 3364 10.1016/j.physleta.2011.07.033 MR2826245 2-s2.0-84860388823 Merdan M. A numeric-analytic method for time-fractional Swift-Hohenberg (S-H) equation with modified Riemann-Liouville derivative Applied Mathematical Modelling. Simulation and Computation for Engineering and Environmental Systems 2013 37 6 4224 4231 10.1016/j.apm.2012.09.003 MR3020566 ZBL06198513 2-s2.0-84872605462 Wu G. C. A fractional variational iteration method for solving fractional nonlinear differential equations Computers & Mathematics with Applications 2011 61 8 2186 2190 10.1016/j.camwa.2010.09.010 MR2785581 2-s2.0-79953735813 Wu G. C. Baleanu D. Discrete fractional logistic map and its chaos Nonlinear Dynamics 2014 75 1-2 283 287 10.1007/s11071-013-1065-7 MR3144852 Wu G. C. Baleanu D. Variational iteration method for the Burgers' flow with fractional derivatives—new Lagrange multipliers Applied Mathematical Modelling 2013 37 9 6183 6190 10.1016/j.apm.2012.12.018 MR3039000 2-s2.0-84875378030 Inc M. The approximate and exact solutions of the space- and time-fractional Burgers equations with initial conditions by variational iteration method Journal of Mathematical Analysis and Applications 2008 345 1 476 484 10.1016/j.jmaa.2008.04.007 MR2422665 ZBL1146.35304 2-s2.0-43949121726 Sakar M. G. Erdogan F. Yıldırım A. Variational iteration method for the time-fractional Fornberg-Whitham equation Computers & Mathematics with Applications 2012 63 9 1382 1388 10.1016/j.camwa.2012.01.031 MR2912063 2-s2.0-84859756253 Ganjiani M. Solution of nonlinear fractional differential equations using homotopy analysis method Applied Mathematical Modelling. Simulation and Computation for Engineering and Environmental Systems 2010 34 6 1634 1641 10.1016/j.apm.2009.09.011 MR2592597 ZBL1193.65147 2-s2.0-75449103955 Rida S. Z. El-Sherbiny H. M. Arafa A. A. M. On the solution of the fractional nonlinear Schrödinger equation Physics Letters A 2008 372 5 553 558 10.1016/j.physleta.2007.06.071 MR2378723 2-s2.0-38049021511 Hao R. Li L. Li Z. Xue W. Zhou G. A new approach to exact soliton solutions and soliton interaction for the nonlinear Schrödinger equation with variable coefficients Optics Communications 2004 236 1–3 79 86 10.1016/j.optcom.2004.03.005 2-s2.0-2342454377 Chen Y. Li B. An extended subequation rational expansion method with symbolic computation and solutions of the nonlinear Schrödinger equation model Nonlinear Analysis: Hybrid Systems 2008 2 2 242 255 10.1016/j.nahs.2006.04.008 MR2400948 2-s2.0-40949140530 Li B. Chen Y. On exact solutions of the nonlinear Schrödinger equations in optical fiber Chaos, Solitons and Fractals 2004 21 1 241 247 10.1016/j.chaos.2003.10.029 MR2033700 2-s2.0-0942299377 Almeida R. Torres D. F. M. Calculus of variations with fractional derivatives and fractional integrals Applied Mathematics Letters 2009 22 12 1816 1820 10.1016/j.aml.2009.07.002 MR2558546 2-s2.0-70349490327 Yang S. Xiao A. Su H. Convergence of the variational iteration method for solving multi-order fractional differential equations Computers & Mathematics with Applications 2010 60 10 2871 2879 10.1016/j.camwa.2010.09.044 MR2734326 2-s2.0-78049255335 Ganjiani M. Solution of nonlinear fractional differential equations using homotopy analysis method Applied Mathematical Modelling 2010 34 6 1634 1641 10.1016/j.apm.2009.09.011 MR2592597 2-s2.0-75449103955 Herzallah M. A. E. Gepreel K. A. Approximate solution to the time-space fractional cubic nonlinear Schrodinger equation Applied Mathematical Modelling 2012 36 11 5678 5685 10.1016/j.apm.2012.01.012 MR2956775 2-s2.0-84863627476 Wazwaz A. A study on linear and nonlinear Schrodinger equations by the variational iteration method Chaos, Solitons and Fractals 2008 37 4 1136 1142 10.1016/j.chaos.2006.10.009 ZBL1148.35353 2-s2.0-40649100751