Bateman–feshbach Tikochinsky and Caldirola–kanai Oscillators with New Fractional Differentiation

In this work, the study of the fractional behavior of the Bateman–Feshbach–Tikochinsky and Caldirola–Kanai oscillators by using different fractional derivatives is presented. We obtained the Euler–Lagrange and the Hamiltonian formalisms in order to represent the dynamic models based on the Liouville–Caputo, Caputo–Fabrizio–Caputo and the new fractional derivative based on the Mittag–Leffler kernel with arbitrary order α. Simulation results are presented in order to show the fractional behavior of the oscillators, and the classical behavior is recovered when α is equal to 1.


Introduction
Several phenomenological models of dissipative systems have been proposed, such as the Bateman-Feshbach-Tikochinsky (BFT) or Caldirola-Kanai (CK) oscillators, the first model consists of a damped and an amplified oscillator, and this one-dimensional system exhibits an exponentially increasing mass with a Lagrangian given by Bateman [1][2][3][4][5].Both quantum damped oscillators have been studied as a model to understand dissipation in quantum theory [6].Bateman suggested the time-dependent Hamiltonian [2] and Caldirola the time dependent Hamiltonian to describe damped oscillations [4].The Caldirola-Kanai oscillator is an open system whose parameters such as mass and frequency are all time dependent, while the Bateman-Feshbach-Tikochinsky oscillator is a closed system whose total energy is conserved and the dissipated energy from the damped oscillator is transferred to amplified one [7,8].The fractional Hamiltonians are non-local and they are associated with dissipative systems [8].There are few definitions of operators with fractional order, the Liouville-Caputo fractional derivative involving a kernel with singularity, and this definition is based on the power law and present singularity at the origin [9].Recently, in order to solve the problem of singularity at the origin, Caputo and Fabrizio used the exponential decay law to construct a derivative with no singularity; however, the used kernel was local [10][11][12][13][14][15][16][17][18].Thus, Atangana and Baleanu used the generalized Mittag-Leffler function to construct a derivative with no-singular and non-local kernel [19][20][21][22].In this paper, we obtain alternative representations of the BFT and CK oscillators by using the Liouville-Caputo, Caputo-Fabrizio-Caputo and the new fractional derivative based in Mittag-Leffler kernel with arbitrary order α.Numerical solutions are based in a Crank-Nicholson scheme.

Fractional Operators
The Adams method is a multi-step method, and this method uses the information of all the previous values, y i , y i−1 , y i−m+1 , in order to calculate y i+1 .This is the difference between the Adams method and the single-step methods, such as the Heun, Taylor and Runge-Kutta numerical schemes, which use only the last value to calculate the next one.There are two types of Adams methods, the Adams-Bashforth and the Adams-Moulton.The combination of these methods results in the predictor-corrector Adams-Bashforth-Moulton Method [23][24][25][26].
The generalization of this method for any order of derivative is called the fractional Adams-Bashforth Method [23] C where α > 0 and C 0 D α t is the Liouville-Caputo operator Equation ( 1) satisfies the following Volterra integral equation The fractional Adams method to solve (1) has been studied firstly by Diethelm, Ford and Freed [24], and this solution scheme is derived as follows: a j,w+1 g(t j , f j ) + a w+1,w+1 g(t w+1 , f P k+1 ) .
The fractional operator proposed by Caputo and Fabrizio in Liouville-Caputo sense (CFC) is expressed as follows [10]: where B(α) = B(0) = B(1) = 1 (is a normalization function).In this sense, the Laplace transform is given by The fractional derivative based in Mittag-Leffler kernel (Atangana-Baleanu fractional operator in Liouville-Caputo sense, ABC) is given as where E α is a Mittag-Leffler function [19].The fractional integral is defined as The Laplace transform of (7) produces 3. Applications

Bateman-Feshbach-Tikochinsky Oscillator
The classical Lagrangian of the BFT oscillator is given by where q 1 is the damped harmonic oscillator coordinate and q 2 corresponds to the time-reversed counterpart, and the parameters m, ρ, K are time independent.The fractional Lagrangian ( 10) is given by and the Lagrange model of fractional order is Now, we can get the generalized momentum as follows: where L F is the Lagrangian of fractional order and i = 1, 2.
The two generalized momentums are given by Applying the Legendre transformation, we obtain the Hamiltonian of fractional order Using the Equation ( 15), we have We define ω = K − ρ 2 4m and the Hamiltonian takes the form The fractional Hamilton model of the BFT oscillator is given by Now, we consider the fractional operators of Liouville-Caputo, Caputo-Fabrizio-Caputo and the fractional derivative based in the Mittag-Leffler kernel.

•
First case.In the Liouville-Caputo sense, we have The numerical approximation of ( 19) is obtained using the algorithm (4).
• Second case.In the Caputo-Fabrizio-Caputo sense, where • Third case.For the fractional derivative based on the Mittag-Leffler kernel, we used the numerical approximation scheme developed in [20] AB where and the system ( 18) is represented by

Numerical Simulations
Figures 1-3 shows the position q 1 = x 2 (t), q 2 = x 1 (t), a D α t x 1 (t) = x 3 (t) and a D α t x 2 (t) = x 4 (t) for systems (19), ( 20) and ( 23), respectively.For the simulation, the following values were considered: m = 5, ρ = 2, K = 0.1 and different values of α, the total simulation time considered is 5 s, and the computational step 1 × 10 −3 .The initial conditions x 1 (0) = 1, x 2 (0) = 0.1, x 3 (0) = 1 and x 4 (0) = 0.5 were considered.The results show that by keeping the parameters constant and by varying α, we obtain different results.The reported results illustrate that the fractional approach is more suitable to describe the complex dynamics of the investigated model.

Caldirola-Kanai Oscillator
We consider a harmonic CK oscillator whose mass depends on time such that m(t) = m exp(sin βγt), in this case, the Lagrangian is given by where m depends explicitly on time, and β and γ are variable parameter and damping factors.
The fractional Lagrangian ( 24) is given by and The generalized momentum is where L F is the Lagrangian of fractional order of ( 24) with i = 1, q 1 = q and p 1 = p.The Hamiltonian of fractional order is obtained using the Legendre transformation where The fractional Hamilton model of the CK oscillator is given by Now, we consider the fractional operators of Liouville-Caputo, Caputo-Fabrizio-Caputo and the fractional derivative based on the Mittag-Leffler kernel.

•
First case.In the Liouville-Caputo sense, we have The numerical approximation of (32) is obtained using algorithm (4).
• Second case.In the Caputo-Fabrizio-Caputo sense, the Adams-Moulton rule for system (31) is given by where • Third case.For the fractional derivative based on the Mittag-Leffler kernel, we have

Numerical Simulations
Figures 4-6 depicted the numerical evaluation of (32)-(34) in Liouville-Caputo, Caputo-Fabrizio-Caputo and the fractional derivative based on the Mittag-Leffler kernel, respectively, considering different values of ω(t) and fractional order γ, for all cases a = 0 and b = 1, and the total simulation time considered is one second and computational step 1 × 10 −5 .It is clear from the figures that the behaviors of the fractional equations strongly depend on the order α of the fractional derivatives, in addition to the form of the function w(t).

Conclusions
Alternative representations of the Bateman-Feshbach-Tikochinsky and Caldirola-Kanai oscillators were studied using fractional operators of Liouville-Caputo type.We derive new solutions of these models using an iterative scheme and via a Crank-Nicholson scheme.The Liouville-Caputo fractional derivative involves a kernel with singularity, and this definition is based on the power law and present singularity at the origin.Recently, Caputo and Fabrizio solved the problem of singularity at the origin and used the exponential decay law to construct a derivative with no singularity; however, the used kernel is local.This derivative therefore has an advantage over the Liouville-Caputo derivative because the full effect of the memory can be portrayed.Atangana and Baleanu suggested two fractional derivatives based on the generalized Mittag-Leffler function.These derivatives with fractional order in Liouville-Caputo and Riemann-Liouville sense have non-singular and non-local kernel and preserve the benefits of the Riemann-Liouville, Liouville-Caputo and Caputo-Fabrizio operators.
Using these fractional operators, the results show that, by keeping the parameters constant and by varying α, we obtain different behaviors.The reported results illustrate that the fractional approach is more suitable to describe the complex dynamics of the investigated models.Finally, we observe novel behaviors that cannot be obtained with standard models and using local derivatives.