key: cord-0846054-lsoi34df authors: Al-Raeei, Marwan; El-Daher, Moustafa Sayem title: Numerical simulation of the space dependent fractional Schrödinger equation for London dispersion potential type date: 2020-07-18 journal: Heliyon DOI: 10.1016/j.heliyon.2020.e04495 sha: bcaf7a65a28dbda9f0c18682eba8ccae13077467 doc_id: 846054 cord_uid: lsoi34df In this study, we apply the definition of one of the fractional derivatives definitions of increasing values of the variable, which is the fractional derivative of Riemann-Liouville, and the numerical-integral methods to find numerical solutions of the fractional Schrödinger equation with the time-independent form for Van Der Walls potential type. We use the dimensionless formalism of the fractional Schrödinger equation in the space-dependent form in case of London dispersion potential in the stationary state. The solutions are found for multiple values of the space-dependent fractional Schrödinger equation parameter with a certain value of the energy. We find that the numerical solutions are physically acceptable for some values of the space dependent fractional parameter of the fractional Schrödinger equation but are not physically acceptable for others for a specific case. The numerical solutions can be applied for the systems that obey London dispersion potential type, which is resulted from the polarization of the instantaneous multi-poles of two moieties, such as soft materials systems and fluids of the inert gases. The linear fractional Schr€ odinger equation or the fractional Schr€ odinger equation is a type of the linear Schr€ odinger equations used in the fractional quantum mechanics and this equation use one of the physical potentials to give the probability. The fractional Schr€ odinger equation in the general space representation form with a space fractional parameter α is given in the linear formalism as follows [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19] : Where h is Planck constant, Ψ(r,t) is the wave function of the system in space-representation, j is the imaginary unit and b H α is Hamiltonian operator in the fractional type which is defined by the following formula: K α is a coefficient, U(r) is the interaction potential of the system and ðÀℏ 2 ΔÞ α=2 is the space fractional operator which is given by: The Φ(Π,t) is the wave function of the system in the momentumrepresentation which relates to the wave function in spacerepresentation by Fourier transforms. In one dimension, the fractional Schr€ odinger equation in the space-dependent form is written as: There is another type of linear fractional Schr€ odinger equation called the time-dependent fractional Schr€ odinger equation and the two types can be derived from the general fractional Schr€ odinger equation which depending on two fractional parameters but in this work we focus on the space-dependent type in the stationary state. In most of literatures, the term fractional Schr€ odinger equation refers to the space-dependent fractional Schr€ odinger equation type so in this work we use this convention. The analytical or exact solutions of some applications of the fractional Schr€ odinger equation in the space form are found in multiple physical cases, the first case is free particle in one and three dimensions [7] , the second case is infinite potential well and barrier penetration [14] , the third case is δ function [7, 8] , the fourth case is linear potential field and Coulomb's potential in one dimension [7, 8, 14] and the fifth case is the energy level of stationary states for one dimension fractional oscillator by semi classical approximation [7] , in addition to a number of other mathematical solutions. In the following, we find an algorithm to solve the space-dependent fractional Schr€ odinger equation in one dimension in case of potential of London dispersion forces which is important in case of some macromolecular such as polymers and proteins, the simple relation that gives potential of London dispersion forces between two atoms or molecules is given as follows [20, 21, 22, 23, 24, 25, 26] : Where I is the first ionization energy, λ is the dipole polarizablitiy and C ij is a constant can be calculated from the first ionization energy and dipole The fractional Schr€ odinger equation in the fractional quantum mechanics is very similar to the traditional Schr€ odinger equation in the traditional quantum mechanics. It is well known that the traditional Schr€ odinger equation is considered the most important equation in quantum mechanics where we can find a lot of properties for a certain system described by a specific potential by solving the traditional Schr€ odinger equation either by an exact solution or by a numerical solution. The same thing can be seen from the fractional Schrodinger equation in the fractional quantum mechanics because if we find a solution of the fractional Schrodinger equation, we can determine a lot of properties of the system. In this study, we find a numerical solutions of the fractional Schrodinger equation in the stationary state in case of London dispersion potential type which is considered an important potential for describing the weak forces between particles which is resulted from the fluctuating polarization of neighbour particles in a specific system which has an interaction between the instantaneous multi-poles particles such as polymers, proteins and inert gases. This study is arranged as follows: in the second section, we illustrate the algorithm which we use for finding the numerical solutions of the fractional Schr€ odinger for London dispersion potential type. In the third section, we insert computational results of the study and the discussion of these results and in the last section, we illustrate a conclusion of the study. The fractional Schr€ odinger equation in the space-dependent form for London dispersion potential can be written as follows: As we see, the Hamiltonian of this case does not explicitly depend on the time and then we write the special solution of the last equation as follows: Ψ ðr; tÞ ¼ e Àj Et ℏ ψðrÞ (7) This leads to: b H α ψðrÞ ¼ E α ψðrÞ (8) Where E α and ψ(r) are the eigenvalues and eigenvectors of the Hamiltonian operator. The previous equation in the dimensionless formalism is written as follows: D α r * ψðr * Þ þ λ * r *À6 ψðr * Þ þ E * ψðr * Þ ¼ 0 Where r* is the dimensionless distance, D α r * is the fractional derivative operator of the α order with respect to the dimensionless distance, ψ(r*) is the wave function in the stationary formalism as a function of the dimensionless distance, λ* is the reduced dipole polarizablitiy and E* is the eigenvalues of the energy in the dimensionless formalism. For purpose of finding the numerical solutions of the previous equation, we apply the fractional derivative of Riemann-Liouville [27, 28, 29] and the Runge-Kutta methods. We use Wolfram Mathematica programing language to code the method and we apply it for multiple values of the space-dependent fractional parameter. We apply the method for the range ]0,10000a 0 ] of the distance. We insert and plot some values of the calculated results in the following section. The full results are ready when requested. Figure 1 shows the algorithm used for this method. We apply the previous method to compute the wave functions of the stationary state, however, if we need to find the total wave function, we apply equation. 7 to find the full wave function. The method must be used only with 1<α < 2 because we only use Riemann-Liouville definition of the fractional derivatives for the last range of the space dependent fractional parameter. So, the method can be applied for all values of the fractional parameter of the space-dependent fractional Schrodinger equation except α equals to the limited value. The wave functions must be physical so it must be continuous, single-valued and limited, i. e. when the distance increases to infinity, the wave function should be decreased. The method can be applied in case of systems that obeys London dispersion forces which have induced multi-poles between moieties such as soft materials system like proteins and some other polymers and Figure 13 . The probability versus distance for α ¼ 1.90. macromolecules. In this work, we applied the numerical method in case of Polyethylene [24, 30] for multiple values of the space dependent fractional parameter in case of the ground state. The values of the wave function were calculated for values of the space dependent fractional parameter of the fractional Schr€ odinger equation from α ¼ 1.10 to α ¼ 1.90. Some of the results for α ¼ 1.10 and α ¼ 1.20 were inserted in Table 1 and for α ¼ 1.30 and α ¼ 1.40 in Table 2 where a is Bohr radius. However, for the full values of the wave functions for the previous space dependent fractional parameters of the fractional Schr€ odinger equation, see supplementary martials file of this article. As we notice from Table 1 , the wave functions values increase with distance for the two values of the space dependent fractional parameter of the fractional Schr€ odinger equation α ¼ 1.10 and α ¼ 1.20, and the same thing we can see from all values less than α ¼ 1.20 which mean they are not acceptable for physical reasons. While, for the other values of the space dependent parameters of the fractional Schr€ odinger equation, i. e. for the values more than α ¼ 1.20 as it can be seen from Table 2 for two values, the values of the wave functions decrease when the distance increase. Also, we can see from the two tables and from the tables illustrated in the supplementary martials file that the most values of the wave functions for the previous values of the space dependent fractional parameters are complex. Besides, we see from Table 2 that the values of the wave function for α ¼ 1.40 go to zero for the greatest values of the distance. We plotted the amplitude of wave function versus the distance for the following four values of the space dependent fractional parameter of the fractional Schr€ odinger equation α ¼ 1.10, α ¼ 1.20, α ¼ 1.21, α ¼ 1.40, α ¼ 1.50, α ¼ 1.60 and α ¼ 1.90, the results were depicted in Figure 2 , We see from the previous figures and tables, all results of the wave functions for the space-dependent fractional parameter of the fractional Schr€ odinger equation greater than 1.2 are limited so they do not give unlimited probability while the results for the space-dependent fractional parameter smaller than or equals to 1.2 increase with the distance. We used a numerical algorithm to solve the space-dependent fractional Schr€ odinger equation in the stationary state for Van Der Walls potential type based on Riemann-Liouville definition of the fractional derivatives in addition to the Runge-Kutta quadrature method. We used the dimensionless form of the space-dependent fractional Schr€ odinger equation in case of London dispersion potential because it is more useful for reducing the consumed time in the computational calculations. We did the computational calculations in Wolfram Mathematica software. We found the wave functions values (Table 1, Table 2 Figure 12 and Figure 13 ) in case of system obeys London dispersion potential in the ground state, for multiple values of the space dependent fractional parameter except for the limited value of the space dependent fractional parameter. We found that the values of the space-dependent fractional parameter which give limited wave functions are greater than 1.2 while the values those are smaller than or equals 1.2 give increasing wave functions so these values are not accepted for physical reasons. We believe that the numerical method we used in this study can be applied for multiple systems interact via London dispersion potential type such as soft martials system and fluid of the inert gases, all we need to apply the method: the parameters of the systems. So, if we find the numerical solutions of the fractional Schr€ odinger equation of London dispersion potential type, i. e. finding the probability, we can find a lot of physical properties of the considered system by applying the integration over the required properties using the probability functions which were found by the numerical solutions. Marwan Al-Raeei: Conceived and designed the experiments; Performed the experiments; Analyzed and interpreted the data; Contributed reagents, materials, analysis tools or data; Wrote the paper. Moustafa Sayem El-Daher: Contributed reagents, materials, analysis tools or data; Wrote the paper. This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors. An algorithm for fractional Schr€ odinger equation in case of Morse potential Fractionals and Quantum Mechanics A numerical method for fractional Schr€ odinger equation of Lennard-Jones potential Fractional quantum mechanics and L evy path integrals On: new optical soliton solutions for nonlinear complex fractional Schr€ odinger equation via new auxiliary equation method and novel (G 0 / G)-expansion method Fractional quantum mechanics in polariton condensates with velocity-dependent mass Fractional Quantum Mechanics Some physical applications of fractional Schr€ odinger equation Fractional Klein-Gordon equation composed of Jumarie fractional dervatives and its interpretation by a smoothness parameter A novel approach for solving fractional Fisher equation using differential transform method A study of fractional Schrodinger equation composed of Jumarie fractional dervative An Improved Collocation Method for Multi-Dimensional Space-Time Variable-Order Fractional Schr€ odinger Equations The second hyperpolarizability of systems described by the space-fractional Schr€ odinger equation Some solutions to the space fractional Schr€ odinger equation using momentum representation method The forecasting of COVID-19 with mortality using SIRD epidemic model for the United States, Russia, China, and the Syrian Arab Republic Structural and quantum properties of van der Waals cluster near the unitary regime Nanoparticles in a nanochannel: van der Waals interaction and diffusion Higher-dimensional fractional time independent Schrodinger equation via fractional dervative with generalised pseudoharmonic potential Significant decrease in thermal conductivity of multi-walled carbon nanotube induced by interwall van der Waals interactions Physics for the Prearatory Year of Medical Colleges Analytical formula of heat capacity in soft matter materials using Lennard-Jones potential Analytical static structure factor for a two-component system interacting via van der Waals potential Using Methods of Statistical Mechanics in the Study of Soft Condensed Matter Materials and Complex Structures Temperature dependence of the specific volume of Lennard-Jones potential and applying in case of polymers and other materials Thermodynamic properties of fluids with Lennard-Jones-Gauss potential from computer simulation and the coupling parameter series expansion Electronic and optical properties of vanderWaals vertical heterostructures based on two-dimensional transition metal dichalcogenides: firstprinciplescalculations Basic Theory of Fractional Differential Equations Fractional Calculus and Fractional Processes with Applications to Financial Economics Single electron states in polyethylene The authors declare no conflict of interest. Supplementary content related to this article has been published online at https://doi.org/10.1016/j.heliyon.2020.e04495.