arxiv: v1 [physics.comp-ph] 2 Jan 2018

Transcription

1 arxiv: v1 [physics.comp-ph] 2 Jan 2018 Solution of the Nonlinear Schrödinger Equation with Defocusing Strength Nonlinearities Through the Laplace-Adomian Decomposition Method O. González-Gaxiola 1, Pedro Franco, R. Bernal-Jaquez Departamento de Matemticas Aplicadas y Sistemas, Universidad Autónoma Metropolitana-Cuajimalpa, Vasco de Quiroga 4871, Santa Fe, Cuajimalpa, Mexico City, Mexico 1 ogonzalez@correo.cua.uam.mx August 12, 2018 Abstract In this work we apply the Adomian decomposition method combined with the Laplace transform (LADM) in order to solve the 1-dimensional nonlinear Schrödinger equation whose nonlinear term presents a nonlinear defocusing strength that varies in the spatial direction. The suggested iterative scheme (LADM) finds the solution without any discretization, linearization or restrictive assumptions. Finally, three numerical examples are presented to demonstrate the reliability and accuracy of the method. 1 Introduction Obtaining the solutions of nonlinear partial differential equations that emerge from scientific and engineering problems constitutes one of the most important, challenging and essential tasks of mathematics that become crucial for the understanding of innumerable phenomena. Under this circumstances, great efforts have been done in order to develop different mathematical approaches that allow the solution of nonlinear equations in an analytical, semi-analytical or numerical way. Among these semi-analytical methods, the one developed by the mathematician George Adomian ( ), has proved to be very useful in many different contexts [4], [5] and [41]. The Adomian decomposition method (ADM) has the advantage that it converges to the exact solution in many important cases and applications. Besides, this method can be handled easily for a wide class of differential equations (ordinary and partial) both linear and nonlinear. In the other hand, the nonlinear Schrödinger equation (NLSE) is one of the most important models of mathematical-physics that describes non only the quantum world 1

2 but also describes many different phenomena in fields as diverse as nonlinear optics, plasma physics, wave propagation in oceans, protein dynamics, laser light propagation in materials with a refractive index sensitive to the incident wave amplitude among many others [23] and [36]. In the last decades, this equation has become fundamental to describe many aspects of the dynamics of the Bose-Einstein condensates such as vortices formation [33], interaction among condensates [24] as well as in the description of the interaction of atom with lasers [17], inusual behavior of solitons in non-homogeneous media and nonlinear excitations in material science [28]. In this paper we will use ADM combined with the Laplace transform in order to find approximate solutions to the nonlinear Schrödinger equation when we consider a function that acts as a defocusing strength of the nonlinear potential that varies in the spatial direction [22, 27, 11, 16]. We will see that the nonlinear term can be easily handled with the help of the so called, Adomian s polynomials. The ADM method combined with use of the Laplace transform, proposed originally in [29], plays an important role among semi-analytic methods and has been used to solve many different nonlinear problems as in [34] in which the method was used to solve a bio-mathematical problem, in [42] was employed to solve a problem from the mechanics of vibrations, in [30] was used to solve a nonlinear problem relevant in chemical engineering and in [38] was used to solve nonlinear integro-differential equations of the Volterra type. Our work is divided in several sections. In The Laplace Adomian Decomposition Method (LADM) section, we present, in a briefly and self-contained way, the LADM method and some references are presented to delve deeper into the subject and to study their mathematical foundations that are not part of the issues of this work. In Solution of the General Nonlinear Schrödinger Equation for Bright and Dark Solitons Through LADM section, we give a brief introduction to the model described by the nonlinear Schrödinger equation, whose nonlinearity term presents a defocusing strength function that varies in the spatial direction, and we will establish that LADM can be used to solve it. In Numerical Examples we will show by means of three examples, the quality and accuracy of our method, comparing the obtained results with the exact solutions that appear in the literature often without giving an explicit method to find them [12] and [25]. Finally, in the Conclusion and Discussion section, we summarize our findings and present our final conclusions. 2 Conclusion and Discussion In the present work we have shown that it is possible to obtain accurate approximate solutions of the nonlinear Schrödinger equation with a defocusing strength function in the nonlinear potential by means of the Laplace Adomian decomposition method (LADM). Our work has a twofold intension, in one hand we want to show that the LADM is a reliable method to solve hard nonlinear problems such as the ones treated here and we illustrate how to use the method through the solution of different examples and on the other hand we have chosen the nonlinear Schrödinger equation equation due to its importance in modeling a wide range of physical phenomena. The method presented in this work is quite competitive with other decomposition meth- 2

3 ods due to the application of the Laplace transform that avoids the integration of polynomials that in every step become of higher order, providing an elegant and easy to apply method that decreases considerably the amount of calculations. Nevertheless, we have to mention, as a part of the conclusions of this work, that the present method has some deficiencies, for instance, the method provides a series solution that has to be truncated in most of the cases (if we can not identify an analytic solution expressed as a series) and besides, the region in which the ADM solution converges to the exact solution (and hence the LADM converges as well) is limited (region in which, usually, convergence is very fast). If one needs the solution to be convergent in an ample region we need to include more terms in the ADM series solution. To show the accuracy and efficiency of this method, we have solved three examples, comparing our results with the exact solution of the equation that was obtained in [25] and [12]. Our results show that LADM produces highly accurate solutions for small values of time, for complicated nonlinear problems. In this manner we can conclude that the LADM is a notable non-sophisticated powerful tool that produces high quality approximate solutions for nonlinear partial differential equations using simple calculations and that attains convergence with only few terms. All the numerical work was accomplished with the Mathematica software package. References [1] Abbaoui, K., Cherruault, Y.: Convergence of Adomian s method applied to differential equations. Comput. Math. Appl. 28(5), (1994). doi: / (94) [2] Abbaoui, K., Cherruault, Y.: New ideas for proving convergence of decomposition methods. Comput. Math. Appl. 29(7), (1995). doi: / (95)00022-q [3] Abushammala, M., Khuri, S. A., Sayfy A.: A novel fixed point iteration method for the solution of third order boundary value problems. Applied Math. and Computation 271, (2015). doi: /j.amc [4] Adomian, G.: Nonlinear Stochastic Operator Equations. Academic Press, Orlando (1986) [5] Adomian, G.: Solving Frontier Problems of Physics: The Decomposition Method. Kluwer Academic Publishers, Boston, MA (1994) [6] Adomian, G., Rach, R.: Linear and nonlinear Schrödinger equations. Foundations of physics 21(8), (1991). doi: /BF [7] Ahmed, I., Chunlai, M., Zheng P.: Exact solution of the (2 + 1) dimensional hyperbolic nonlinear Schrödinger equation by Adomian decomposition method. Malaya J. Mat. 2(2), (2014) [8] Alomari, A. K., Noorani, M.S.M., Nazar, R.: Explicit series solutions of some linear and nonlinear Schrödinger equations via the homotopy analysis 3

4 method. Commun. Nonlinear Sci. Numer. Simul. 14(4), (2009). doi: /j.cnsns [9] Babolian, E., Javadi, Sh.: New method for calculating Adomian polynomials. Applied Math. and Computation 153, (2004). doi: /S (03) [10] Biazar, J., Ghazvini, H.: Exact solutions for non-linear Schrödinger equations by He s homotopy perturbation method. Phys. Lett. A. 366(1-2), (2007). doi: /j.physleta [11] Bloch, I., Dalibard, J., Zwerger, W.: Many-body physics with ultracold gases. Rev. Mod. Phys. 80(3), (2008). doi: /RevModPhys [12] Borovkova, O.V., Kartashov, Y. V., Torner, L., Malomed, B. A.: Bright solitons from defocusing nonlinearities. Phys. Rev. E 84, (R) (2011). doi: /PhysRevE [13] Bratsos, A., Ehrhardt, M., Famelis, I. Th.: A discrete Adomian decomposition method for discrete nonlinear Schrödinger equations. Appl. Math. Comput. 197(1), (2008). doi: /j.amc [14] Cherruault, Y.: Convergence of Adomian s method. Kybernetes 18(2), (1989) [15] Cherruault, Y., Adomian, G.: Decomposition methods: a new proof of convergence. Math. Comput. Modelling 18(12), (1993). doi: / (93)90233-O [16] Chin, C., Grimm, R., Julienne, P., Tiesinga, E.: Feshbach resonances in ultracold gases. Rev. Mod. Phys. 82(2), (2010). doi: /RevMod- Phys [17] Drummond P.D., Kheruntsyan, K. V.: Asymptotic solutions to the Gross- Pitaevskii gain equation: Growth of a Bose-Einstein condensate. Phys. Rev. A 63(1), (2000). doi: /PhysRevA [18] Duan, J. S.: Convenient analytic recurrence algorithms for the Adomian polynomials. Appl. Math. Comput. 217, (2011). doi: /j.amc [19] Duan, J. S.: New recurrence algorithms for the nonclassic Adomian polynomials. Appl. Math. Comput. 62, (2011). doi: /j.camwa [20] Edeki, S. O., Akinlabi, G. O., Adeosun, S. A.: Analytic and numerical solutions of time-fractional linear Schrödinger equation. Comm. in Math. and Appl. 7(1), 1-10 (2016) [21] Fariborzi, M. A., Naghshband, S.: On convergence of homotopy analysis method to solve the Schrödinger equation with a power law nonlinearity. Int. J. Industrial Math. 5(4), (2013) 4

5 [22] Fedichev, P.O., Kagan, Y., Shlyapnikov, G.V., Walraven, J. T. M.: Influence of nearly resonant light on the scattering length in low-temperature atomic gases. Phys. Rev. Lett. 77(14), (1996). doi: /PhysRevLett [23] Fordy, A.P.: Soliton Theory: A Survey of Results. Manchester University Press, Manchester (1990) [24] Hoston, W., You, L.: Interference of two condensates. Phys. Rev. A 53(6), 4254 (1996). doi: /PhysRevA [25] Hou Tsang, C., Malomed, B. A., Wing Chow, K.: Exact solutions for periodic and solitary matter waves in nonlinear lattices. Discrete and continuous dynamical systems-series S 4(5), (2011). doi: /dcdss [26] Hukriede, J., Runde, D., Kip, D.: Fabrication and application of holographic Bragg gratings in lithium niobate channel waveguides. J. Phys. D. 36, 1-16 (2003). [27] Inouye, S., Andrews, M.R., Stenger, Miesner, H. J., Stamper-Kurn, D. M., Ketterle, W: Observation of Feshbach resonances in a Bose-Einstein condensate. Nature 392, (1998) [28] Kartashov, Y. V., Malomed, B.A.,Torner, L.: Solitons in nonlinear lattices. Rev. Mod. Phys. 83(1), (2011). doi: /RevModPhys [29] Khuri, S. A.: A Laplace decomposition algorithm applied to a class of nonlinear differential equations. J. of Applied Mathematics 1(4), (2001). doi: /S X [30] Khuri, S. A.: A new approach to Bratu s problem. Appl. Math. Comput. 147, (2004). doi: /S (02) [31] Khuri, S. A., Sayfy A.: A Laplace variational iteration strategy for the solution of differential equations. Appl. Math. Letters 25(12), (2012). doi: /j.aml [32] Khuri, S. A.: A new approach to the cubic Schrödinger equation: An application of the decomposition technique. Appl. Math. Comput. 97, (1998) [33] Martikainen, J. P., Suominen, K. A., Santos, L., Schulte, T., Sanpera A.: Generation and evolution of vortex-antivortex pairs in Bose-Einstein condensates. Phys. Rev. A 64(6), (2001). doi: /PhysRevA [34] Ongun, M. Y.: The Laplace Adomian decomposition method for solving a model for HIV infection of CD4 + T cells. Mathematical and Computer Modelling 53, (2011). doi: /j.mcm [35] Sadighi, A., Ganji D. D.: Analytic treatment of linear and nonlinear Schrödinger equations: A study with homotopy-perturbation and Adomian decomposition methods. Phys. Letters A 372(4), (2008). doi: /j.physleta

6 [36] Vázquez, L., Streit, L., Pérez-García, V.M. (Eds.): Nonlinear Schrödinger and Klein-Gordon systems: Theory and Applications. World Scientific, Singapore (1996) [37] Wazwaz, A. M.: A new algorithm for calculating Adomian polynomials for nonlinear operators. Appl. Math. and Computation 111(1), (2000). doi: /s (99) [38] Wazwaz, A. M.: The combined Laplace transform-adomian decomposition method for handling nonlinear Volterra integro-differential equations. Appl. Math. Comput. 216(4), (2010). doi: /j.amc [39] Wazwaz, A. M., El-Sayed, S. M.: A new modification of the Adomian decomposition method for linear and nonlinear operators. Appl. Math. Comput. 122(3), (2001). doi: /S (00) [40] Wazwaz, A. M.: A reliable technique for solving linear and nonlinear Schrödinger equations by Adomian decomposition method. Bull. of the Intitute. of Math. Acad. Sinica 29(2), (2011) [41] Wazwaz, A.M.: Partial Differential Equations: Methods and Applications. Balkema Publishers, Lisse (2002) [42] Yusufoǧlu, E.: Numerical solution of Duffing equation by the Laplace decomposition algorithm. Appl. Math. Comput. 177, (2006). doi: /j.amc