Abstract
A quantum mechanical system that mimics the behavior of a classical harmonic oscillator in the quantum domain is called a simple harmonic quantum oscillator. The time-independent Schrödinger equation describes the quantum harmonic oscillator, and its eigenstates are quantized energy values that correspond to various energy levels. In this work, we first fractionalize the time-independent Schrödinger equation, and then we solve the generated problem with the use of the Adomian decomposition approach. It has been shown that fractional quantum harmonic oscillators can be handled effectively using the proposed approach, and their behavior can then be better understood. The effectiveness of the method is validated by a number of numerical comparisons.
Highlights
- Fractionalizing the time-independent Schrödinger equation using Caputo fractional differentiator.
- Solving the fractional-order Schrödinger equation with the use of the Adomian decomposition approach.
- The effectiveness of the Adomian decomposition approach is validated by a number of numerical comparisons.
1. Introduction
In recent years, fractional calculus – a branch of classical calculus that deals with instructions on non-integer integration and differentiation – has become a fascinating area of study. The notion of fractional operators, which emerged nearly concurrently with their classical equivalents, has garnered considerable attention owing to their multifarious applicability across multiple fields of mathematics and science. It has been demonstrated that fractional differential equations are effective instruments for understanding and simulating intricate engineering, chemical, and physical processes [1-5].
The Schrödinger equation is a key tool in the field of quantum mechanics for explaining how quantum systems behave. This equation, which was created by Erwin Schrödinger in 1926, uses time as an independent variable to describe quantum systems. Time-independent Schrödinger equations (TISE) and time-dependent Schrödinger equations (TDSE) are the two categories into which Schrödinger wave equations fall [5, 6].
Numerous authors have offered both analytical and numerical solutions for TISE and TDSE. The Schrödinger issue has been solved using modern analytical techniques such the Elzaki decomposition methodology, homotopy perturbation, and Adomian decomposition methods. On the other hand, techniques such as the Numerov algorithm for harmonic and linear potentials have been used to produce approximative numerical solutions.
In this work, we study the semi-analytical solution of the Adomian decomposition approach for time-independent Schrödinger equations. The ADM provides an accurate solution and is well-known for its effectiveness in solving Sturm-Liouville issues, ordinary and partial differential equations, nonlinear and stochastic problems, and more. When compared to other approaches, the ADM solves the TISE more quickly and accurately while still offering a series solution [7-10].
The structure of this article is as follows: In part 2, we establish the notations and present the essential concepts, including the fractional operators of the Riemann-Liouville integral and derivative, and the Caputo derivative, which will be utilized throughout this essay. Section 3, employs the ADM to solve the Schrödinger equation for the simple fractional harmonic oscillator, converting it into a Hermite differential equation. In Section 4, we depict some numerical comparisons between the ADM’s solution and the solution reported in some references. Finally, Section 5 summarizes the conclusion of this work.
2. Preliminaries and notion
In this section, we go over the basic ideas of fractional calculus in detail, which forms the basis of our investigation. We also explore the properties of certain operators that are important to our work [11-14].
Definition 2.1. The Riemann-Liouville fractional integral of of order , where , , , of a function is given as follows:
The following features of the Riemann-Liouville integral are worth noting:
Definition 2.2. Let , such that is positive integer and , the Riemann Liouville derivative of fractional of order is given as follows:
Definition 2.3. Suppose that , , . The Caputo fractional differential operator of order , is given as:
Remarks. The Caputo fractional derivative meets the following properties:
1. The power rule property, i.e.:
2. The constant property, i.e.:
3. The interpolation property, i.e.:
4. The linearity property, i.e.:
5. The non-commutation property, i.e.:
3. Theory
One typical method for modeling systems with viscoelasticity in damping terms is to incorporate fractional derivatives to introduce memory effects in the system. The mathematical framework offered by fractional calculus makes it possible to include memory effects and produce a more realistic depiction of the behavior of the system.
The damping term in the differential equation can be written as a fractional derivative of the displacement or velocity variable when modeling viscoelasticity damping with fractional derivatives. The memory effects in the damping system are captured by this fractional derivative, enabling a more accurate depiction of the system’s reaction to outside influences.
The core of this study is covered in this section, where we introduce our unique method for solving the Schrödinger equation for the basic fractional harmonic oscillator. In order to do this, we make use of the ADM, a reliable numerical technique that is well-known for its accuracy and efficiency in solving differential equations. Our objective is to transform the original Schrödinger equation into the Hermite differential equation, which is a more manageable form and has fractional order derivatives [15-17].
An important model for studying quantum systems with fractional characteristics is the simple fractional harmonic oscillator. Our goal is to achieve precise answers by utilizing the ADM, which will provide insight into the fascinating occurrences seen in the equation under consideration. The ADM is a perfect fit for this study because of its adaptability in handling nonlinearity and the complexities of fractional calculus [18-20].
In this context, the expression for the time-dependent Schrodinger's equation for a basic fractional harmonic oscillator is:
where , are constants, is the energy. It should be observed here that the above model corresponds to inertia terms in the second-order mechanistic counterpart model when alpha = 1. In tha same regard, we should also note that since we are dealing with a physical problem, we will not be considering the positive exponent solution that diverges at. It was noted, meanwhile, that Eq. (13) can be changed to the following form:
where . Now define and replace the dummy with . This reduces Eq. (14) to the form shown below:
where ,
The ADM is used to tackle this problem, by considering the following form:
Operating to both sides of Eq. (16) yields:
or:
Again, by taking to both sides of Eq. (17), we obtain:
The result is as follows:
where and are arbitrary constants. Following the ADM sense leads to the general solution for , which would be as:
where:
We will show some of the major calculations used above. The following steps show how we calculate up to . In this regard, we can find:
which means:
or:
The next term can be written as:
In other words, we have:
Once more, the next term can be written as:
This is equivalent to the following assertion:
Similarly, we can have:
which consequently gives:
If we continue in this manner, we got the desired solution reported in Eq. (19). To show the validity of our solution, we compare the ADM’s solution with the solution reported in [20] in Fig. 1 and Fig. 2 by taking 1 and = {0.5, 1.5}. Obviously, we can observe that the proposed solution coincides with the solution we compare, and this confirms the validity of the proposed solution.
Fig. 1ADM’s solution vs. the solution in [20] according to λ= 0.5
Fig. 2ADM’s solution vs. the solution in [20] according to λ= 1.5
Fig. 3Absolute error between the ADM’s solution and the solution in [20]
For more description, we plot in Fig. 3 and Fig. 4, the absolute errors gained from the performance comparison. These figures demonstrate the advantage of the suggested method of fractionalizing the model in question. However, the limitation of our proposed approach lies in the lack of capabilities to obtain a general form of that can generate all other components of the desired solution, and this will be left to the future for more consideration.
Fig. 4Absolute error between the ADM’s solution and the solution in [20] when λ= 0.5
4. Conclusions
In this study, a generalized form of the time-independent Schrödinger equation has been successfully proposed in its fractional simple harmonic quantum oscillator. The ADM has been used to obtain a semi-analytical solution for such an oscillator. This solution has proved its validity with an existing solution in [20] when a classical case of the generalized form is taken into consideration.
References
-
H. Goldstein, C. P. Poole, and J. L. Safko, Classical Mechanics. Boston: Addison Wesley, 2002.
-
T. Hamadneh et al., “General methods to synchronize fractional discrete reaction-diffusion systems applied to the glycolysis model,” Fractal and Fractional, Vol. 7, No. 11, p. 828, Nov. 2023, https://doi.org/10.3390/fractalfract7110828
-
G. B. Arfken and H. Weber, Mathematical Methods for Physicists. San Diego: Academic Press, 2001.
-
J. Biazar, E. Babolian, and R. Islam, “Solution of the system of ordinary differential equations by Adomian decomposition method,” Applied Mathematics and Computation, Vol. 147, No. 3, pp. 713–719, Jan. 2004, https://doi.org/10.1016/s0096-3003(02)00806-8
-
A.-M. Wazwaz, “Approximate solutions to boundary value problems of higher order by the modified decomposition method,” Computers and Mathematics with Applications, Vol. 40, No. 6-7, pp. 679–691, Sep. 2000, https://doi.org/10.1016/s0898-1221(00)00187-5
-
B. Zhang and J. Lu, “Exact solutions of homogeneous partial differential equation by a new Adomian decomposition method,” Procedia Environmental Sciences, Vol. 11, pp. 440–446, Jan. 2011, https://doi.org/10.1016/j.proenv.2011.12.070
-
S. Somali and G. Gokmen, “Adomian decomposition method for nonlinear Sturm-Liouville problems,” Surveys in Mathematics and its Applications, Vol. 2, pp. 11–20, 2007.
-
A. Cheniguel and A. Ayadi, “Solving heat equation by the Adomian decomposition method,” in Proceedings of the World Congress on Engineering, 2011.
-
M. Dehghan, “The use of Adomian decomposition method for solving the one-dimensional parabolic equation with non-local boundary specifications,” International Journal of Computer Mathematics, Vol. 81, No. 1, pp. 25–34, Jun. 2010, https://doi.org/10.1080/0020716031000112321
-
D. N. Khan Marwat and S. Asghar, “Solution of the heat equation with variable properties by two-step Adomian decomposition method,” Mathematical and Computer Modelling, Vol. 48, No. 1-2, pp. 83–90, Jul. 2008, https://doi.org/10.1016/j.mcm.2007.09.003
-
I. G. Rochdi Jebari and A. Boukricha, “Adomian decomposition method for solving nonlinear diffusion equation with convection term,” International Journal of Pure and Applied Sciences and Technology, Vol. 12, pp. 49–58, 2012.
-
F.-K. Yin, W.-Y. Han, and J.-Q. Song, “Modified Laplace decomposition method for Lane-Emden type differential equations,” International Journal of Applied Physics and Mathematics, Vol. 3, No. 2, pp. 98–102, Jan. 2013, https://doi.org/10.7763/ijapm.2013.v3.184
-
F. Mainardi, “Fractional Calculus,” Fractals and Fractional Calculus in Continuum Mechanics, pp. 291–348, Jan. 1997, https://doi.org/10.1007/978-3-7091-2664-6_7
-
Q. Wang, “Homotopy perturbation method for fractional KdV equation,” Applied Mathematics and Computation, Vol. 190, No. 2, pp. 1795–1802, Jul. 2007, https://doi.org/10.1016/j.amc.2007.02.065
-
R. L. Burden and J. D. Faires, Numerical Analysis. Pacific Grove: Brooks/Cole, 2011.
-
R. Courant, Methods of Mathematical Physics, Volume II: Partial Differential Equations. New York: Inter science Publishers Inc., 1966.
-
R. I. Nuruddeen, “Elzaki decomposition method and its applications in solving linear and nonlinear Schrodinger equations,” Sohag Journal of Mathematics, Vol. 4, No. 2, pp. 31–35, May 2017, https://doi.org/10.18576/sjm/040201
-
J. Biazar, R. Ansari, K. Hosseini, and P. Gholamin, “Solution of the linear and nonlinear Schrodinger equations using homotopy perturbation and Adomian decomposition methods,” International Mathematical Forum. Journal for Theory and Applications, Vol. 3, pp. 1891–1897, 2008.
-
A. Sadighi and D. D. Ganji, “Analytic treatment of linear and nonlinear Schrödinger equations: A study with homotopy-perturbation and Adomian decomposition methods,” Physics Letters A, Vol. 372, No. 4, pp. 465–469, Jan. 2008, https://doi.org/10.1016/j.physleta.2007.07.065
-
A. A. Anulo, “On analyzing numerical solution of time independent Schrodinger equation,” Global Scientific Journals, Vol. 7, No. 6, pp. 190–211, 2019.
-
A. K. Jaradat, A. A. Obeidat, M. A. Gharaibeh, and M. K. Hasan Qaseer, “Adomian decomposition approach to solve the simple harmonic quantum oscillator,” International Journal of Applied Engineering Research, Vol. 13, No. 2, pp. 1056–1059, 2018.
-
I. M. Batiha, A. Bataihah, A. A. Al-Nana, S. Alshorm, I. H. Jebril, and A. Zraiqat, “A numerical scheme for dealing with fractional initial value problem,” International Journal of Innovative Computing, Information and Control, Vol. 19, No. 3, pp. 763–774, Jan. 2023, https://doi.org/10.24507/ijicic.19.03.763
-
I. M. Batiha, S. Alshorm, A. Al-Husban, R. Saadeh, G. Gharib, and S. Momani, “The n-point composite fractional formula for approximating Riemann-Liouville integrator,” Symmetry, Vol. 15, No. 4, p. 938, Apr. 2023, https://doi.org/10.3390/sym15040938
-
I. M. Batiha, O. Y. Ababneh, A. A. Al-Nana, W. G. Alshanti, S. Alshorm, and S. Momani, “A numerical implementation of fractional-order PID controllers for autonomous vehicles,” Axioms, Vol. 12, No. 3, p. 306, Mar. 2023, https://doi.org/10.3390/axioms12030306
-
I. M. Batiha, S. Momani, S. Alshorm, and A. Ouannas, “Numerical solutions of stochastic differential equation using modified three-point fractional formula,” in 2023 International Conference on Fractional Differentiation and Its Applications (ICFDA), pp. 1–5, Mar. 2023, https://doi.org/10.1109/icfda58234.2023.10153192
-
I. M. Batiha et al., “Common fixed-point theorem in non-Archimedean menger PM-spaces using CLR property with application to functional equations,” IAENG International Journal of Applied Mathematics, Vol. 53, No. 1, 2023.
-
I. M. Batiha, Z. Chebana, T.-E. Oussaeif, A. Ouannas, S. Alshorm, and A. Zraiqat, “Solvability and dynamics of superlinear reaction diffusion problem with integral condition,” IAENG International Journal of Applied Mathematics, Vol. 53, No. 1, pp. 1–9, 2023.
-
I. M. Batiha, N. Alamarat, S. Alshorm, O. Y. Ababneh, and S. Momani, “Semi-analytical solution to a coupled linear incommensurate system of fractional differential equations,” Nonlinear Functional Analysis and Applications, Vol. 28, No. 2, pp. 449–471, Jun. 2023, https://doi.org/10.22771/nfaa.2023.28.02.09
-
I. M. Batiha, S. Alshorm, I. Jebril, A. Zraiqat, Zaid Momani, and S. Momani, “Modified 5-point fractional formula with Richardson extrapolation,” AIMS Mathematics, Vol. 8, No. 4, pp. 9520–9534, Jan. 2023, https://doi.org/10.3934/math.2023480
About this article
The authors have not disclosed any funding.
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.
Conceptualization, Iqbal M. Batiha; methodology, Iqbal H. Jebril and Abeer A. Al-Nana; validation, Shameseddin Alshorm; formal analysis, Iqbal M. Batiha.; investigation, Iqbal H. Jebril and Abeer A. Al-Nana; resources, Shameseddin Alshorm; data curation, Iqbal M. Batiha; writing-original draft, Abeer A. Al-Nana; visualization, Iqbal H. Jebril; supervision, Iqbal M. Batiha. All authors have read and agreed to the published version of the manuscript.
The authors declare that they have no conflict of interest.