﻿ Closed Form Exact Solutions to the Higher Dimensional Fractional Schrodinger Equation via the Modified Simple Equation Method

Journal of Applied Mathematics and Physics
Vol.06 No.01(2018), Article ID:81681,13 pages
10.4236/jamp.2018.61009

Closed Form Exact Solutions to the Higher Dimensional Fractional Schrodinger Equation via the Modified Simple Equation Method

M. Nurul Islam1, M. Ali Akbar2

1Department of Mathematics, Islamic University, Kushtia, Bangladesh

2Department of Applied Mathematics, University of Rajshahi, Rajshahi, Bangladesh    Received: August 17, 2017; Accepted: January 8, 2018; Published: January 11, 2018

ABSTRACT

In this article, we investigate some exact wave solutions to the higher dimensional time-fractional Schrodinger equation, an important equation in quantum mechanics. The fractional Schrodinger equation further precisely describes the quantum state of a physical system changes in time. In order to determine the solutions a suitable transformation is considered to transmute the equations into a simpler ordinary differential equation (ODE) namely fractional complex transformation. We then use the modified simple equation (MSE) method to obtain new and further general exact wave solutions. The MSE method is more powerful and can be used in other works to establish completely new solutions for other kind of nonlinear fractional differential equations arising in mathematical physics. The affect of obtaining parameters for its definite values which are examined from the solutions of two dimensional and three dimensional time-fractional Schrodinger equations are discussed and therefore might be useful in different physical applications where the equations arise in this article.

Keywords:

Modified Simple Equation (MSE) Method, Fractional Differential Equation, Nonlinear Evolution Equations, Higher Dimensional Schrodinger Equation, Traveling Wave Transformation 1. Introduction

The differential equations with fractional order have recently become a valuable tool to the modeling of numerous tangible events and it has gained importance and popularity to the researchers. The modeling of a tangible incident plays an important role on the history of the previous time which can also be successfully achieved by using fractional calculus. The use of fractional differentiation for the mathematical modeling of real-world physical problems has been widespread in recent years, for example, the modeling of earthquake, the fluid dynamic traffic model with fractional derivatives, and measurement of viscoelastic material properties. Applications of fractional differential equations in other fields, like quantum mechanics, electricity, plasma physics, chemical kinematics, optical fibers and related area are also felt. The fractional calculus has dominated in every field of sciences and engineering. In quantum mechanics the fractional Schrodinger equation   is an equation that describes how the quantum state of a physical system changes in time. Thus, searching traveling wave solutions of fractional nonlinear evolution equations (NLEEs) plays a fundamental role. To know the internal mechanism of complex physical phenomena exact solutions of nonlinear fractional differential equations is very much important. As a result, recently some useful methods have been established and enhanced for obtaining exact solution to the fractional evolution equations such as, the extended direct algebraic function method   , the F-expansion method  , the Adomian decomposition method  , the homotopy perturbation method     , the tanh-function method  , the Sine-Cosine method  , the Jacobi elliptic method  , the finite difference method  , the variational iteration method   , the variational method  , the Fourier transform technique  , the modified decomposition method  , the Laplace transform technique  , the operational calculus method in  , the exp-function method   , the $\left({G}^{\prime }/G\right)$ -expansion method    , the modified simple equation method (MSE)  -  , the $\mathrm{exp}\left(-\phi \left(\eta \right)\right)$ -expansion method  , the sub equation method  , the multiple exp-function method   , the simplest equation method  , the direct algebraic function method     , the extended auxiliary equation method  etc.

The aim of this article is to examine the further general and new exact solution of higher dimensional time-fractional Schrodinger equation by making use of the modified simple equation method  -  and discuss effect of the included parameters to the obtained solutions. We also discuss that the attained solutions might be useful and realistic to analyze the fractional quantum system for the time fractional two and three dimensional Schrodinger equation. We also have studied the behavior of emerging parameters which affect the obtained solutions and also describe how the quantum state of a physical system changes in time. In earlier literature, the time fractional Schrodinger equation is investigated through the first integral method  , the F-expansion method  , the Fourier transformation method  , and the Laplace transformation method  .

The MSE method is a recently developed efficient, potential and rising method to investigate wave solutions to the nonlinear fractional equations. Its finding results are straightforward, efficient, systematic, and no need to use the symbolic computation software to manipulate the algebraic equations.

The rest of the article is processed as follows: In Section 2, we explain the Jumarie modified Riemann-Liouville derivative. In Section 3, we describe the outline of the MSE method. In Section 4, we investigate new and further general solutions to the time fractional equations mentioned above. In Section 5, we draw our conclusions.

2. Modified Riemann-Liouville Derivative

The Jumarie’s modified Riemann-Liouville derivative of order $\alpha$ is defined as follows  :

${D}_{x}^{\alpha }f\left(x\right)=\left\{\begin{array}{l}\frac{1}{\Gamma \left(-\alpha \right)}\frac{\text{d}}{\text{d}x}\underset{0}{\overset{x}{\int }}{\left(x-\xi \right)}^{-\alpha -1}\left(f\left(\xi \right)-f\left(0\right)\right)\text{d}\xi ,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\alpha <0,\\ \frac{1}{\Gamma \left(1-\alpha \right)}\frac{\text{d}}{\text{d}x}\underset{0}{\overset{x}{\int }}{\left(x-\xi \right)}^{-\alpha }\left(f\left(\xi \right)-f\left(0\right)\right)\text{d}\xi ,\text{\hspace{0.17em}}\text{\hspace{0.17em}}0<\alpha <1,\\ {\left[{f}^{\left(\alpha -n\right)}\left(x\right)\right]}^{\left(n\right)},\text{\hspace{0.17em}}\text{\hspace{0.17em}}n\le \alpha (2.1)

Some properties for the proposed modified Riemann-Liouville derivative are listed  as follows:

${D}_{x}^{\alpha }{x}^{\gamma }=\frac{\Gamma \left(\gamma +1\right)}{\Gamma \left(\gamma +1-\alpha \right)}{x}^{\gamma -\alpha },\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\gamma >0,$ (2.2)

${D}_{x}^{\alpha }\left(f\left(x\right)g\left(x\right)\right)=g\left(x\right){D}_{x}^{\alpha }f\left(x\right)+f\left(x\right){D}_{x}^{\alpha }g\left(x\right),$ (2.3)

${D}_{x}^{\alpha }f\left[g\left(x\right)\right]={f}_{g}\left[g\left(x\right)\right]{D}_{x}^{\alpha }g\left(x\right)={D}_{g}^{\alpha }f\left[g\left(x\right)\right]{\left({g}^{\prime }\left(x\right)\right)}^{\alpha },$ (2.4)

The above formulae play an important role in fractional calculus and also fractional differential equations.

3. Outline of the Method

Let us consider the nonlinear fractional evolution equation in the form:

$H\left(u,{D}_{t}^{\alpha }u,{D}_{x}^{\alpha }u,{D}_{t}^{2\alpha }u,{D}_{x}^{2\alpha }u,\cdots \right)=0,$ (3.1)

where $u=u\left(x,t\right)$ is wave function, $H$ is a polynomial in $u\left(x,t\right)$ and its partial derivatives, which consist of the highest order derivatives and nonlinear terms of the highest order, and the subscripts denote partial derivatives. To obtain the solution of (3.1) by using the MSE method  -  , we have to execute the subsequent steps:

Step 1: We assume, $u\left(x,t\right)=\varphi \left(\xi \right){\text{e}}^{i\eta }$ and the traveling wave variable,

$\xi =k\left(x-\frac{c{t}^{\alpha }}{\Gamma \left(1+\alpha \right)}\right)$ and $\eta =mx-\frac{w{t}^{\alpha }}{\Gamma \left(1+\alpha \right)}$ , (3.2)

permits us to transform the Equation (3.1) into the following ordinary differential equation (ODE):

$G\left(u,{u}^{\prime },{u}^{″},\cdots \right)=0$ , (3.3)

where $G$ is a polynomial in $u\left(\xi \right)$ and its derivatives, wherein ${u}^{\prime }\left(\xi \right)=\frac{\text{d}u}{\text{d}\xi }$ .

Step 2: We assume that the solution of Equation (3.3) can be revealed in the form:

$u\left(\xi \right)=\sum _{i=0}^{N}{a}_{i}{\left[\frac{{s}^{\prime }\left(\xi \right)}{s\left(\xi \right)}\right]}^{i},$ (3.4)

where ${a}_{i},\left(i=0,1,2,\cdots ,N\right)$ are unknown constants to be evaluated, such that ${a}_{N}\ne 0$ , and $s\left(\xi \right)$ is an unidentified function to be estimated. In Jacobi elliptic function method, $\left({G}^{\prime }/G\right)$ -expansion method, F-expansion method, Riccati equation method, extended tanh-function method etc., the solutions are pre-defined or the solutions are presented in terms of some well-known differential equations, but in the MSE method, $s\left(\xi \right)$ is neither pre-defined nor a solution of any pre-defined differential equation. This is the individuality and uniqueness of the MSE method. Therefore, some useful and realistic solutions might be obtained by this method.

Step 3: The positive integer N arises in Equation (3.4) can be determined by balancing the highest order nonlinear terms and the derivatives of highest order occur in Equation (3.4).

Step 4: Inserting (3.4) into (3.3) and simplifying for the function $s\left(\xi \right)$ , we

obtain a polynomial of $\left(\frac{1}{s\left(\xi \right)}\right)$ . From the resulted polynomial, we equate all

the coefficients of ${\left(s\left(\xi \right)\right)}^{-i},\left(i=0,1,2,\cdots ,N\right)$ to zero. This procedure yields a system of algebraic and differential equations which can be solved for determining ${a}_{i}\left(i=0,1,2,\cdots ,N\right)$ , $s\left(\xi \right)$ and the other necessary parameters. This completes the determination of the solutions to the Equation (3.1).

4. Applications of the Method

In this section, we will examine the new and further general useful solutions to the time fractional (2+1)-dimensional and (3+1)-dimensional Schrodinger equations.

4.1. The (2+1)-Dimensional Schrodinger Equation

In this section, we investigate some applicable close form traveling wave solutions to the time fractional (2+1)-dimensional Schrodinger equation by making use the MSE method. Let us consider the time fractional (2+1)-dimensional Schrodinger equation of the form:

$i{D}_{t}^{\alpha }u=\frac{1}{2}{u}_{xx}\left({u}_{xx}+{u}_{yy}\right)+pu+{|u|}^{2}u,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}0<\alpha \le 1,$ (4.1.1)

where $\alpha$ and $p$ are emerging parameters. The Schrodinger equation is a mathematical equation that describes the variation over time of a physical structure on the fractional quantum system, as for instance wave particle duality is noteworthy. We can use this equation as a mathematical formula for the study of quantum mechanical system (The equation is a mathematical formula for the study of quantum mechanical systems). By means of the traveling transformation (3.2), the Equation (4.1.1) is converted into the following nonlinear ODE:

$\left(2km-ck\right){\varphi }^{\prime }i+\left({k}^{2}{\varphi }^{″}+w\varphi -{m}^{2}\varphi -p\varphi -{\varphi }^{3}\right)=0$ , (4.1.2)

Equating real and imaginary part of Equation (4.1.2), we obtain

$\left(2km-ck\right){\varphi }^{\prime }=0$ , (4.1.3)

And

${k}^{2}{\varphi }^{″}+w\varphi -{m}^{2}\varphi -p\varphi -{\varphi }^{3}=0$ . (4.1.4)

As ${\varphi }^{\prime }\ne 0$ , from Equation (4.1.3), it can be easily obtained $c=2m$ .

Now, balancing the linear term of the highest order derivative term ${\varphi }^{″}$ and the nonlinear term of the highest order ${\varphi }^{3}$ occurring in (4.1.4), yields $N=1$ . Thus, the solution of Equation (4.1.4) is the form:

$\varphi \left(\xi \right)={a}_{0}+{a}_{1}\frac{{s}^{\prime }\left(\xi \right)}{s\left(\xi \right)}$ , (4.1.5)

where ${a}_{0}$ and ${a}_{1}$ are constants to be determined, such that ${a}_{1}\ne 0$ , and $s\left(\xi \right)$ is an unknown function to be evaluated. Now, it is simple to estimate the following:

${\varphi }^{\prime }\left(\xi \right)=-\frac{{a}_{1}{\left({s}^{\prime }\right)}^{2}}{{s}^{2}}+\frac{{a}_{1}{s}^{″}}{s}$ , (4.1.6)

${\varphi }^{″}\left(\xi \right)=\frac{2{a}_{1}{\left({s}^{\prime }\right)}^{3}}{{s}^{3}}-\frac{3{a}_{1}{s}^{\prime }{s}^{″}}{{s}^{2}}+\frac{{a}_{1}{s}^{‴}}{s}$ , (4.1.7)

Substituting the values of $\varphi \left(\xi \right),{\varphi }^{″}\left(\xi \right)$ from (4.1.5) and (4.1.7) into (4.1.4) and then equating the coefficients of ${s}^{0},{s}^{-1},{s}^{-2},{s}^{-3},{s}^{-4}$ to zero, we respectively obtain

$-{m}^{2}{a}_{0}-p{a}_{0}+w{a}_{0}+{a}_{0}^{3}=0$ (4.1.8)

$-{m}^{2}{a}_{1}{s}^{\prime }\left(\xi \right)-p{a}_{1}{s}^{\prime }\left(\xi \right)+w{a}_{1}{s}^{\prime }\left(\xi \right)-3{a}_{0}^{2}{a}_{1}{s}^{\prime }\left(\xi \right)+{k}^{2}{a}_{1}{s}^{‴}\left(\xi \right)=0$ (4.1.9)

$3{a}_{0}{a}_{1}^{2}{{s}^{\prime }}^{2}\left(\xi \right)+3{k}^{2}{a}_{1}{s}^{\prime }\left(\xi \right){s}^{″}\left(\xi \right)=0$ (4.1.10)

$2{k}^{2}{a}_{1}{{s}^{\prime }}^{3}\left(\xi \right)-{a}_{1}^{3}{{s}^{\prime }}^{3}\left(\xi \right)=0$ (4.1.11)

From (4.1.8) and (4.1.11), we attain

${a}_{0}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{a}_{0}=±\sqrt{-{m}^{2}-p+w}$ and ${a}_{1}=±\sqrt{2}k;\text{\hspace{0.17em}}{a}_{1}\ne 0$

Case 1: When ${a}_{0}=0$ , Equation (4.1.2) produces an absurd solution. Hence, the case is not accepted.

Case 2: When ${a}_{0}=\sqrt{-{m}^{2}-p+w}$ , ${a}_{1}=±\sqrt{2}k$ and $c=2m$ , then from Equations. (4.1.9) and (4.1.10), we obtain

$s\left(\xi \right)={c}_{1}-\frac{{c}_{2}{k}^{2}}{2N}{\text{e}}^{\mp 2M\xi }$ , where $M=\frac{\sqrt{-{m}^{2}-p+w}}{\sqrt{2}k},\text{\hspace{0.17em}}N={m}^{2}+p-w$ .

Substituting the value of ${a}_{0},{a}_{1}$ and $s\left(\xi \right)$ into the solution (4.1.5), it yields

$\varphi \left(\xi \right)=\frac{\sqrt{2}Mk\left({c}_{2}{k}^{2}{\text{e}}^{\mp 2M\xi }+2{c}_{1}N\right)}{\left(-{c}_{2}{k}^{2}{\text{e}}^{\mp 2M\xi }+2{c}_{1}N\right)}$ (4.1.12)

Converting the solution (4.1.12) from exponential to trigonometric function, we attain

$\varphi \left(\xi \right)=\frac{\sqrt{2}Mk\left(\left({c}_{2}{k}^{2}+2{c}_{1}N\right)\mathrm{cosh}\left(M\xi \right)±\left(-{c}_{2}{k}^{2}+2{c}_{1}N\right)\mathrm{sinh}\left(M\xi \right)\right)}{\left(-{c}_{2}{k}^{2}+2{c}_{1}N\right)\mathrm{cosh}\left(M\xi \right)±\left({c}_{2}{k}^{2}+2{c}_{1}N\right)\mathrm{sinh}\left(M\xi \right)}$ (4.1.13)

Since, ${c}_{1}$ and ${c}_{2}$ are arbitrary constant, one can choose their values arbitrarily. Therefore, if we choose, ${c}_{1}={k}^{2}$ and ${c}_{2}=±2N$ , from solution (4.1.13), we obtain

$\varphi \left(\xi \right)=±\sqrt{2}Mk\mathrm{coth}\left(M\xi \right)$ (4.1.14)

$\varphi \left(\xi \right)=±\sqrt{2}Mk\mathrm{tanh}\left(M\xi \right)$ (4.1.15)

Again when ${a}_{0}=-\sqrt{-{m}^{2}-p+w}$ , ${a}_{1}=±\sqrt{2}k$ and $c=2m$ , then from Equations (4.1.9) and (4.1.10), we obtain

$s\left(\xi \right)={c}_{1}-\frac{{c}_{2}{k}^{2}}{2N}{\text{e}}^{±2M\xi }$ , where $M=\frac{\sqrt{-{m}^{2}-p+w}}{\sqrt{2}k},\text{\hspace{0.17em}}N={m}^{2}+p-w$

Inserting the values of ${a}_{0},{a}_{1}$ and $s\left(\xi \right)$ into the solution (4.1.5), it yields

$\varphi \left(\xi \right)=-\frac{\sqrt{2}Mk\left({c}_{2}{k}^{2}{\text{e}}^{±2M\xi }+2{c}_{1}N\right)}{\left(-{c}_{2}{k}^{2}{\text{e}}^{±2M\xi }+2N{c}_{1}\right)}$ (4.1.16)

Transforming the solution (4.1.16) from exponential to trigonometric function, it becomes

$\varphi \left(\xi \right)=-\frac{\sqrt{2}Mk\left(\left({c}_{2}{k}^{2}+2{c}_{1}N\right)\mathrm{cosh}\left(M\xi \right)\mp \left(-{c}_{2}{k}^{2}+2{c}_{1}N\right)\mathrm{sinh}\left(M\xi \right)\right)}{\left(-{c}_{2}{k}^{2}+2{c}_{1}N\right)\mathrm{cosh}\left(M\xi \right)\mp \left({c}_{2}{k}^{2}+2{c}_{1}N\right)\mathrm{sinh}\left(M\xi \right)}$ (4.1.17)

Here ${c}_{1}$ and ${c}_{2}$ are arbitrary constants, so one can select their values arbitrarily. Thus, if we select, ${c}_{1}={k}^{2}$ and ${c}_{2}=±2N$ , from the solution (4.1.17), we obtain

$\varphi \left(\xi \right)=±\sqrt{2}Mk\mathrm{coth}\left(M\xi \right)$ (4.1.18)

And

$\varphi \left(\xi \right)=±\sqrt{2}Mk\mathrm{tanh}\left(M\xi \right)$ (4.1.19)

Therefore, combining the solutions (4.1.14), (4.1.15), (4.1.18) and (4.1.19), we obtain the required solutions for this case as follows:

$\varphi \left(\xi \right)=±\sqrt{2}Mk\mathrm{coth}\left(M\xi \right)$ (4.1.20)

And

$\varphi \left(\xi \right)=±\sqrt{2}Mk\mathrm{tanh}\left(M\xi \right)$ (4.1.21)

Now, making use of the fractional wave variable (3.2) into solution (4.1.20), we obtain

$\varphi \left(\xi \right)=±\sqrt{2}Mk\mathrm{coth}\left(Mk\left(x-\frac{c{t}^{\alpha }}{\Gamma \left(1+\alpha \right)}\right)\right)$ (4.1.22)

And the solution (4.1.21), becomes

$\varphi \left(\xi \right)=±\sqrt{2}Mk\mathrm{tanh}\left(Mk\left(x-\frac{c{t}^{\alpha }}{\Gamma \left(1+\alpha \right)}\right)\right)$ (4.1.23)

Substituting $M=\frac{\sqrt{-{m}^{2}-p+w}}{\sqrt{2}k}$ , solution (4.1.22) yields

$\varphi \left(\xi \right)=±\sqrt{-{m}^{2}-p+w}\mathrm{coth}\left(\frac{\sqrt{-{m}^{2}-p+w}}{\sqrt{2}}\left(x-\frac{c{t}^{\alpha }}{\Gamma \left(1+\alpha \right)}\right)\right)$ (4.1.24)

And the solution (4.1.23) yields

$\varphi \left(\xi \right)=±\sqrt{-{m}^{2}-p+w}\mathrm{tanh}\left(\frac{\sqrt{-{m}^{2}-p+w}}{\sqrt{2}}\left(x-\frac{c{t}^{\alpha }}{\Gamma \left(1+\alpha \right)}\right)\right)$ (4.1.25)

The solutions attained in (4.1.24) and (4.1.25) are new and further general than the existing solutions. If we choose alternative values of ${c}_{1}$ and ${c}_{2}$ , further closed form analytical solutions to the (2+1)-dimensional nonlinear time fractional Schrodinger equation can be extracted, but for simplicity and conciseness the remaining solutions have not been marked out.

4.2. The (3+1)-Dimensional Schrodinger Equation

In this section, we will use the MSE method to obtain new exact solution to the time-fractional two-dimensional Schrödinger equation. Consider the time- fractional two-dimensional Schrödinger equation is of the form:

$i{D}_{t}^{\alpha }u=\frac{1}{2}{u}_{xx}\left({u}_{xx}+{u}_{yy}+{u}_{zz}\right)+pu+{|u|}^{2}u,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}0<\alpha \le 1,$ (4.2.1)

where $\alpha$ and $p$ are emerging parameters. It arises as a description of the influence on the fractional quantum system. Using the traveling transformation (3.2), the Equation (4.2.1) becomes in the following nonlinear ODE:

$\left(6km-2ck\right){\varphi }^{\prime }i+\left(3{k}^{2}{\varphi }^{″}+2w\varphi -3{m}^{2}\varphi -2p\varphi -2{\varphi }^{3}\right)=0$ (4.2.2)

Separating real and imaginary part of Equation (4.2.2), we obtain

$\left(6km-2ck\right){\varphi }^{\prime }=0$ , (4.2.3)

and

$\left(3{k}^{2}{\varphi }^{″}+2w\varphi -3{m}^{2}\varphi -2p\varphi -2{\varphi }^{3}\right)=0$ . (4.2.4)

Balancing the highest order derivative term ${\varphi }^{″}$ and the nonlinear term of the highest order ${\varphi }^{3}$ occurring in (4.2.4) yields $N=1$ . Thus, the solution of Equation (4.2.4) takes formal form:

$\varphi \left(\xi \right)={a}_{0}+{a}_{1}\frac{{s}^{\prime }\left(\xi \right)}{s\left(\xi \right)}$ (4.2.5)

where ${a}_{0}$ and ${a}_{1}$ are constants to be evaluated such that ${a}_{1}\ne 0$ , and $s\left(\xi \right)$ is an unknown function to be determined. Now, it is simple to compute the followings:

${\varphi }^{\prime }\left(\xi \right)=-\frac{{a}_{1}{\left({s}^{\prime }\right)}^{2}}{{s}^{2}}+\frac{{a}_{1}{s}^{″}}{s}$ (4.2.6)

${\varphi }^{″}\left(\xi \right)=\frac{2{a}_{1}{\left({s}^{\prime }\right)}^{3}}{{s}^{3}}-\frac{3{a}_{1}{s}^{\prime }{s}^{″}}{{s}^{2}}+\frac{{a}_{1}{s}^{‴}}{s}$ (4.2.7)

Inserting the values of $\varphi \left(\xi \right),{\varphi }^{″}\left(\xi \right)$ from (4.2.5) and (4.2.7) into Equation (4.2.4) and then setting the coefficients of ${s}^{0},{s}^{-1},{s}^{-2},{s}^{-3},{s}^{-4}$ equal to zero, we respectively obtain

$-3{m}^{2}{a}_{0}-2p{a}_{0}+2w{a}_{0}-2{a}_{0}^{3}=0$ (4.2.8)

$-3{m}^{2}{a}_{1}{s}^{\prime }\left(\xi \right)-2p{a}_{1}{s}^{\prime }\left(\xi \right)+2w{a}_{1}{s}^{\prime }\left(\xi \right)-6{a}_{0}^{2}{a}_{1}{s}^{\prime }\left(\xi \right)+3{k}^{2}{a}_{1}{s}^{‴}\left(\xi \right)=0$ (4.2.9)

$-6{a}_{0}{a}_{1}^{2}{{s}^{\prime }}^{2}\left(\xi \right)-9{k}^{2}{a}_{1}{s}^{\prime }\left(\xi \right){s}^{″}\left(\xi \right)=0$ (4.2.10)

$6{k}^{2}{a}_{1}{{s}^{\prime }}^{3}\left(\xi \right)-2{a}_{1}^{3}{{s}^{\prime }}^{3}\left(\xi \right)=0$ (4.2.11)

As ${\varphi }^{\prime }\ne 0$ , from (4.2.3) it can be easily obtained $c=3m$ .

From Equations. (4.2.8) and Equation (4.2.11), we attain

${a}_{0}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{a}_{0}=±\frac{\sqrt{-3{m}^{2}-2p+2w}}{\sqrt{2}}$ and ${a}_{1}=±\sqrt{3}\text{ }k$ since ${a}_{1}\ne 0$ .

Case 1: When ${a}_{0}=0$ , Equation (4.2.2) provides an absurd solution. Hence, the case has not been accepted.

Case 2: When ${a}_{0}=\frac{\sqrt{-3{m}^{2}-2p+2w}}{\sqrt{2}}$ , ${a}_{1}=±\sqrt{3}k$ and $c=3m$ , from (4.2.9) and (4.2.10), we attain

$s\left(\xi \right)={c}_{1}-\frac{3{k}^{2}{c}_{2}}{N}{\text{e}}^{±2M\xi }$ , where $M=\frac{\sqrt{-3{m}^{2}-2p+2w}}{\sqrt{6}k}$ , $N=6{m}^{2}+4p-4w$

Therefore, substituting the values of ${a}_{0},{a}_{1}$ and $s\left(\xi \right)$ into the solution (4.2.5), we obtain

$\varphi \left(\xi \right)=\frac{\sqrt{3}Mk\left(3{c}_{2}{k}^{2}{\text{e}}^{\mp 2M\xi }+N{c}_{1}\right)}{-3{c}_{2}{k}^{2}{\text{e}}^{\mp 2M\xi }+N{c}_{1}}$ (4.2.12)

Converting the solution (4.2.12) from exponential to trigonometric function, we obtain

$\varphi \left(\xi \right)=-\frac{\sqrt{3}Mk\left(\left(3{c}_{2}{k}^{2}+{c}_{1}N\right)\mathrm{cosh}\left(M\xi \right)\mp \left(3{c}_{2}{k}^{2}-{c}_{1}N\right)\mathrm{sinh}\left(M\xi \right)\right)}{\left(3{c}_{2}{k}^{2}-{c}_{1}N\right)\mathrm{cosh}\left(M\xi \right)\mp \left(3{c}_{2}{k}^{2}+{c}_{1}N\right)\mathrm{sinh}\left(M\xi \right)}$ (4.2.13)

Here ${c}_{1}$ and ${c}_{2}$ are arbitrary constants. Since, ${c}_{1}$ and ${c}_{2}$ are arbitrary constants one might choose their values arbitrarily. Therefore, if we choose, and , from (4.2.13) we obtain (4.2.14)

and (4.2.15)

Again, when , and , from (4.2.9) and (4.2.10) we obtain , where Thus, from (4.2.5) we obtain (4.2.16)

Shifting the solution (4.2.16) from exponential to trigonometric function, we attain (4.2.17)

where and are integral constants. Since and are arbitrary constants, so one might pick their values randomly. Now, if we pick, and , from (4.2.17) we obtain (4.2.18)

and

(4.2.19)

Therefore, comparing the solutions (4.2.14), (4.2.15), (4.2.18) and (4.2.19), we obtain the next solutions:

(4.2.20)

and

(4.2.21)

Now, making use of the fractional wave variable (3.2) into solutions (4.2.20) and (4.2.21), we obtain

(4.2.22)

and

(4.2.23)

Putting the value, solutions (4.2.22) and (4.2.23) respectively become

(4.2.24)

and

(4.2.25)

The solutions attained (4.2.24) and (4.2.25) are new and more general than the existing solutions. If we choose alternative values of and, further closed form analytical solutions to the three dimensional nonlinear time-fractional Schrodinger equation can be extracted, but for simplicity and conciseness the residual solutions have not been marked out.

5. Conclusion

In this article, we have examined new and further general closed form solitons to time fractional two dimensional and three dimensional Schrodinger equations by means of the efficient technique known as modified simple equation (MSE) method. The solutions are attained in general form and definite values of the included parameters yield diverse known soliton solutions. The attained solutions might be useful to the influence on the fractional quantum system for the time fractional two dimensional and three dimensional Schrodinger equations. And we also have studied the behavior of emerging parameters which are affecting the physical system on the consider equations for the obtaining solutions. When the parameters take certain special values, the solitary waves are derived from the traveling waves. The established results show that the MSE method is more powerful, unified and can be used for many other fractional equations to get feasible solutions of the tangible incidents.

Cite this paper

Nurul Islam, M. and Ali Akbar, M. (2018) Closed Form Exact Solutions to the Higher Dimensional Fractional Schrodinger Equation via the Modified Simple Equation Method. Journal of Applied Mathematics and Physics, 6, 90-102. https://doi.org/10.4236/jamp.2018.61009

References

1. 1. Jiang, X., Qi, H. and Xu, M. (2011) Exact Solutions of Fractional Schrodinger-Like Equation with a Nonlocal Term. Journal of Mathematical Physics, 52, 042105. https://doi.org/10.1063/1.3576189

2. 2. Ablowitw, M.J. and Clarkson, P.A. (1991) Solutions, Nonlinear Evolution Equations and Inverse Scattering. Cambridge University Press, Cambridge. https://doi.org/10.1017/CBO9780511623998

3. 3. Seadawy, A.R. (2014) Stability Analysis for Zakharov-Kuznetsov Equation of Weakly Nonlinear Ion-Acoustic Wave in Plasma. Computers & Mathematics with Applications, 67, 172-180. https://doi.org/10.1016/j.camwa.2013.11.001

4. 4. Seadaway, A.R. (2016) Ion Acoustic Solitary Wave Solutions of Two-Dimensional Nonlinear Kadomtsev-Petviashvili-Burgers Equation in Quantum Plasma. Mathematical Methods in the Applied Sciences, 40, 1598-1607. https://doi.org/10.1002/mma.4081

5. 5. Zhao, Y.M. (2013) F-Expansion Method and Its Application for Finding New Exact Solutions to the Kudryashov-Sinelshchikov Equation. Journal of Applied Mathematics, Article ID 895760, 7 p.

6. 6. El-Sayad, A.M.A., El-Kallad, I.L. and Ziada, E.A.A. (2010) Adomian Solution of Multi-Dimensional Nonlinear Equations of Arbitrary Orders. International Journal of Applied Mathematics and Mechanics, 6, 38-52.

7. 7. Alomari, A.K. and Hashim, I. (2011) Analysis of Fully Developed Flow and Heat Transfer in a Vertical Channel with Prescribed Wall Heat Fluxes by the Homotopy Analysis Method. International Journal for Numerical Methods in Fluids, 67, 805-819. https://doi.org/10.1002/fld.2388

8. 8. Chen, Z. and Jiang, W. (2011) Piecewise Homotopy Perturbation Method for Solving Linear and Nonlinear Weakly Singular VIE of Second Order. Applied Mathematics and Computation, 217, 7790-7798. https://doi.org/10.1016/j.amc.2011.02.086

9. 9. Mehrabinezhad, M. and Saberi-Nadjafi, J. (2010) Application of He’s Homotopy Perturbation Method to Linear Programming Problem. International Journal of Computer Mathematics, 88, 341-347. https://doi.org/10.1080/00207160903443763

10. 10. Gupta, S., Kumar, D. and Singh, J. (2013) Application of He’s Homotopy Perturbation Method for Solving Nonlinear Wave-Like Equations with Variable Coefficients. International Journal of Advances in Applied Mathematics and Mechanics, 1, 65-79.

11. 11. Wazwaz, A.M. (2004) The Tanh-Function Method for Travelling Wave Solutions of Nonlinear Equations. Applied Mathematics and Computation, 154, 713-723. https://doi.org/10.1016/S0096-3003(03)00745-8

12. 12. Jawad, A.J.M. (2012) The Sine-Cosine Function Method for the Exact Solutions of Nonlinear Partial Differential Equations. International Journal of Recent Research and Applied Studies, 13, 186-191.

13. 13. Zheng, B. and Feng, Q. (2014) The Jacobi Elliptic Equation Method for Solving Fractional Partial Differential Equations. Abstract and Applied Analysis, 2014, Article ID 249071, 9 p.

14. 14. Dehghan, M. (2000) A Finite Difference Method for a Non-Local Boundary Value Problem for Two Dimensional Heat Equations. Applied Mathematics and Computation, 112, 133-142. https://doi.org/10.1016/S0096-3003(99)00055-7

15. 15. Wazwaz, A.M. (2007) The Variational Iteration Method for Analytic Treatment for Linear and Nonlinear ODEs. Applied Mathematics and Computation, 212, 120-134. https://doi.org/10.1016/j.amc.2009.02.003

16. 16. Wazwaz, A.M. (2007) The Variational Iteration Method: A Powerful Scheme for Handling Linear and Nonlinear Diffusion Equations. Computers & Mathematics with Applications, 54, 933-939. https://doi.org/10.1016/j.camwa.2006.12.039

17. 17. Helal, M.A. and Seadawy, A.R. (2009) Variational Method for the Derivative Nonlinear Schrodinger Equation with Computational Applications. Physica Scripta, 80, Article ID: 035004.

18. 18. Glushkov, E., Glushkov, N. and Chen, C.S. (2007) Semi Analytical Solution to Heat Transfer Problems using Fourier Transform Technique, Radial Basis Functions, and the Method of Fundamental Solutions. Numerical Heat Transfer, Part B: Fundamentals, 52, 1-28. https://doi.org/10.1080/10407790701443859

19. 19. Wazwaz, A.M. (2006) The Modified Decomposition Method and Pade Approximants for a Boundary Layer Equation in Unbounded Domain. Applied Mathematics and Computation, 177, 737-744. https://doi.org/10.1016/j.amc.2005.09.102

20. 20. Abassy, T.A., El-Tawil, M.A. and El-Zoheiry, H. (2007) Exact Solution of Some Nonlinear Partial Differential Equations using the Variational Iteration Method Linked with Laplace Transforms and the Pade Technique. Computers & Mathematics with Applications, 54, 940-954. https://doi.org/10.1016/j.camwa.2006.12.067

21. 21. Hilfer, R., Luchko, Y. and Tomovski, Z. (2009) Operational Method for the Solution of Fractional Differential Equations with Generalized Riemann-Liouville Fractional Derivatives. Fractional Calculus and Applied Analysis, 12, 299-318.

22. 22. Dai, C.Q. and Zhang, J.F. (2009) Application of He’s Exp-Function Method to the Stochastic mKdV Equation. International Journal of Nonlinear Sciences and Numerical Simulation, 10, 675-680. https://doi.org/10.1515/IJNSNS.2009.10.5.675

23. 23. Zhang, W.M. and Tian, L.X. (2009) Generalized Solitary Solution and Periodic Solution of the Combined KdV-mKdV Equation with Variable Coefficients using the Exp-Function Method. International Journal of Nonlinear Sciences and Numerical Simulation, 10, 711-715. https://doi.org/10.1515/IJNSNS.2009.10.6.711

24. 24. Akbar, M.A., Ali, N.H.M. and Zayed, E.M.E. (2012) A Generalized and Improved -Expansion Method for Nonlinear Evolution Equations. Mathematical Problems in Engineering, 2012, Article ID: 459879. https://doi.org/10.1155/2012/459879

25. 25. Zhang, S., Tong, J. and Wang, W. (2008) A Generalized -Expansion Method for the mKdV Equation with Variable Coefficients. Physics Letters A, 372, 2254-2257. https://doi.org/10.1016/j.physleta.2007.11.026

26. 26. Alam, M.N. and Akbar, M.A. (2014) Application of the New Approach of Generalized -Expansion Method to Find Exact Solutions of Nonlinear PDEs in Mathematical Physics. BIBECHANA, 10, 58-70.

27. 27. Bekir, A., Kaplan, M. and Guner, O. (2015) A Novel Modified Simple Equation Method and Its Application to Some Nonlinear Evolution Equation System. AIP Conference Proceedings, 1611, 30-36. https://doi.org/10.1063/1.4893799

28. 28. Zayed, E.M.E. (2013) The Enhanced Modified Simple Equation Method for Solving Nonlinear Evolution Equations with Variable Coefficients. AIP Conference Proceeding, 1558, 1999-2004.

29. 29. Akther, J. and Akbar, M.A. (2016) Solitary Wave Solution to Two Nonlinear Evolution Equations via the Modified Simple Equation Method. New Trends in Mathematical Sciences, 4, 12-26. https://doi.org/10.20852/ntmsci.2016422033

30. 30. Ayati, Z., Moradi, M. and Mirzazadeh, M. (2015) Application of Modified Simple Equation Method to Burgers, Huxley and Burgers-Huxeley Equations. Iranian Journal of Numerical Analysis and Optimization, 5, 59-73.

31. 31. Zayed, E.M.E. and Ai-Nowehy, A.G. (2016) The Modified Simple Equation Method, the Exp-Function Method, and the Method of Solution Ansatz for Solving the Long-Short Wave Resonance Equations. Zeitschrift für Naturforschung A, 71, 10.

32. 32. Khan, M.A., Akbar, M.A. and Belgacem, F.B.M. (2016) Solitary Wave Solutions for the Boussinesq and Fisher Equations by the Modified Simple Equation Method. Mathematics Letters, 2, 1-18.

33. 33. Zayed, E.M.E., Amer, Y.A. and Al-Nowehy, A.G. (2016) The Modified Simple Equation Method and the Multiple Ex-Function Method for Solving Nonlinear Fractional Sharma-Tasso-Olver Equation. Acta Mathematicae Applicatae Sinica, 32, 793-812. https://doi.org/10.1007/s10255-016-0590-9

34. 34. Khan, K. and Akbar, M.A. (2013) Exact and Solitary Wave Solutions for the Tzitzeica-Dodd Bullough and the Modified KdV-Zakharov-Kuznetsov Equations using the Modified Simple Equation Method. Ain Shams Engineering Journal, 4, 903-909. https://doi.org/10.1016/j.asej.2013.01.010

35. 35. Ali, A., Iqbal, M.A. and Mohyud-Din, S.T. (2016) Traveling Wave Solutions of Generalized Zakharov-Kuznetsov-Benjamin-Bona-Mahony and Simplified Modified Form of Camassa-Holm Equation Method by using Method. Egyptian Journal of Basic and Applied Sciences, 3, 134-140. https://doi.org/10.1016/j.ejbas.2016.01.001

36. 36. Yasar, E., Yildirim, Y. and Khalique, C.M. (2016) Lie Symmetry Analysis, Conservation Laws and Exact Solutions of the Seventh-Order Time Fractional Sawada-Kotera-Ito Equation. Results in Physics, 6, 322-328. https://doi.org/10.1016/j.rinp.2016.06.003

37. 37. Yildirim, Y., Yasar, E. and Adem, A.R. (2017) A Multiple Exp-Function Method for the Three Model Equations of Shallow Water Waves. Nonlinear Dynamics, 89, 2291-2297. https://doi.org/10.1007/s11071-017-3588-9

38. 38. Yildirim, Y. and Yasar, E. (2017) Multiple Exp-Function Method for Soliton Solutions of Nomlinear Evolution Equations. Chinese Physics B, 26, Article ID: 070201. https://doi.org/10.1088/1674-1056/26/7/070201

39. 39. Giresunlu, I.B., Ozkan, Y.S. and Yasar, E. (2017) On the Exact Solutions, Lie Symmetry Analysis, and Conservation Laws of Schamel-Korteweg-de Vries Equation. John Wiley Sons Ltd., Hoboken.

40. 40. Seadawy, A.R. (2015) Fractional Solitary Wave Solutions of the Nonlinear Higher-Order Extended KdV Equation in a Stratified Shear Flow: Part I. Computers & Mathematics with Applications, 70, 345-352. https://doi.org/10.1016/j.camwa.2015.04.015

41. 41. Seadawy, A.R. (2015) Three-Dimensional Nonlinear Modified Zakharov-Kuznetsov Equation of Ion-Acoustic Waves in a Magnetized Plasma. Computers & Mathematics with Applications, 71, 201-212.

42. 42. Seadawy, A.R. and El-Rashidy, K. (2013) Traveling Wave Solutions for Some Coupled Nonlinear Evolution Equations. Mathematical and Computer Modelling, 57, 1371-1379. https://doi.org/10.1016/j.mcm.2012.11.026

43. 43. Seadawy, A.R. (2016) Stability Analysis Solutions of Nonlinear Three-Dimensional Modified Kortweg-de Vries-Zakharov-Kuznetsov Equation in Magnetized Electron-Positron Plasma. Physica A, 455, 44-51. https://doi.org/10.1016/j.physa.2016.02.061

44. 44. Seadawy, A.R. (2017) Traveling Wave Solutions of a Weakly Nonlinear Two-Dimensional and High-Order Kadompsev-Petviashvili Dynamical Equation for Dispersive Shallow Water. The European Physical Journal Plus, 132, 29. https://doi.org/10.1140/epjp/i2017-11313-4

45. 45. Moosaei, H., Mirzazadeh, M. and Yildirim, A. (2011) Exact Solutions to the Perturbed Nonlinear Schrodinger’s Equation with Kerr Law Nonlinearity by using the First Integral Method. Nonlinear Analysis: Modelling and Control, 16, 332-339.

46. 46. Filiz, A., Ekici, M. and Sonmezoglu, A. (2014) F-Expansion Method and New Exact Solutions of the Schrodinger-KdV Equation. The Scientific World Journal, 2014, Article ID: 534063. https://doi.org/10.1155/2014/534063

47. 47. El-Tawil, A. and El-Zoheiry, H. (2007) Exact Solutions of Some Nonlinear Partial Differential Equations using the Variational Iteration Method Linked with Laplace Transforms and the Pade Technique. Computers & Mathematics with Applications, 54, 940-954. https://doi.org/10.1016/j.camwa.2006.12.067

48. 48. Jumarie, G. (2006) Modified Riemann-Liouville Derivative and Fractional Taylor Series of Non-Differentiable Functions Further Results. Computers & Mathematics with Applications, 51, 1367-1376. https://doi.org/10.1016/j.camwa.2006.02.001