On the solutions of fractional differential equations with modified Mittag-Leffler kernel and Dirac Delta function: Analytical results and numerical simulations

Abstract

In this paper we define, for the first time, the modified fractional derivative with Mittage-Leffler kernel of Riemann–Liouville (R-L) type of arbitrary order delta. We derive the infinite series representations for the modified derivatives of R-L and Caputo types and present a relationship between them. We also investigate the modified derivatives for the Dirac delta functions, and study related fractional differential equations. Explicit solutions were presented for linear fractional differential equations with constant coefficients via the Laplace transform. A fractional model with the modified derivative is considered and numerical simulations were presented.

Introduction

In recent years several fractional derivatives with nonsingular kernels were introduced, such as the Caputo-Fabrizio derivative which involves a kernel of exponential type [1] and the Atangana-Baleanu derivative which involves a Mittag-Leffler kernel [2]. These types of derivatives were utilized to model several dynamical systems and were extensively studied by several authors. It was noted that fractional differential equations (FDEs) with fractional derivatives involving nonsingular kernels admit a limitation in their solutions, see Lemmas 3.1, 3.2 and 4.1 in [3], Lemma 3.4 in [4], Proposition 2.1 in [5] and the extensive discussion in [6]. Mainly, extra necessary conditions were imposed to verify the solvability of such equations, which may effect the validity of such derivatives in modeling real life problems. To overcome this issue, recently Al-Refai and Baleanu in [7] have extended the Atangana-Baleanu derivative of Caputo type (ABC) to obtain the modified fractional derivative with Mittage-Leffler kernel of Caputo type which involves a singular kernel. They proved that in general and as a result they presented a non-zero solution to the Cauchy problem

(0.1)

such that Later on, the modified derivative of Caputo type was utilized to model several dynamical systems ([8–13]). A general hybrid coupled system of FDEs with the modified derivative was investigated in [14], which involves several dynamical systems in the literature as particular cases. Several numerical schemes were developed to tackle fractional order systems with the modified derivatives based on the Lagrange polynomials in [9], the Laplace Adomian decomposition method in [8], the Gaussian elimination combined with Taylor’s expansion in [10] and on operational matrix approach in [15]. The basic theory of systems of FDEs with the modified derivative was discussed in [16], analytical solutions in closed forms were obtained for constant coefficients systems and a numerical scheme based on the collocation method was developed for nonlinear systems. The idea of obtaining the modified derivative was extended to other types of fractional derivatives with nonsingular kernels, see [17, 18].

Definition 0.1. [7] Let the modified derivative with Mittag-Leffler kernel of Caputo type and order , is given by

(0.2)

where and denotes the convolution of two functions, and is a normalization function with

We stress on the fact that the kernel in the modified derivative possesses an integrable singularity at the origin. The integral operator corresponding to the modified derivative (0.2) is the same as the integral operator of the original ABC, and here we use the modified integral operator presented in [19]. Operators in the Caputo sense are widely used by researchers in both theoretical and application viewpoints, see [20, 21] The aim of the current study is to address the challenges in solving fractional differential equations that involve fractional derivatives with non-singular kernels. To the best of our knowledge, this is the first study on fractional differential equations with modified derivative of R-L type. We remark that the modified Atangana-Baleanu derivative of the RL-type is an extension of the Atangana-Baleanu derivative, designed to operate in a broader space. This allows for the solvability of related fractional differential equations without the need for additional, unnecessary conditions, as demonstrated in Sect 4.

In this paper, we present the modified fractional derivative with Mittage-Leffler kernels of R-L type and study related problems. Section 1 is devoted to the extension of the modified derivatives. In Sect 2, we present infinite series representations for the modified derivatives of R-L and Caputo types and present the relation between them. In Sect 3, we present the modified derivative of the Dirac delta function and solve related FDEs analytically. An application is presented in Sect 4, with numerical simulations. Finally, some concluding remarks and suggestions for future work are presented in Sect 5.

The following known formulas will be used throughout the text [22, 23].

(0.3)(0.4)(0.5)(0.6)(0.7)

where , denotes the Mittag-Leffler function of two parameters, and is the Laplace transform.

1 The modified fractional derivative with Mittage-Leffler kernel of RL-type

We start with the extension for Let the Atangana-Baleanu derivative of RL-type is given by

(1.1)

Let then integration by parts of the above equation and using (0.6) we have

Because and we arrive at

Definition 1.1. Let the modified fractional derivative with Mittage-Leffler kernel of RL-type and order is given by

(1.2)

where

We apply the above approach to obtain the modified higher order derivatives of order and Let we have

(1.3)

Integration by parts will lead to

Definition 1.2. For and the modified fractional derivative with Mittage-Leffler kernel of RL-type and order is given by

(1.4)

2 Infinite series representations

We derive infinite series representations of the modified derivatives in terms of the R-L fractional integral, and establish a relationship between them. We apply the following facts about the R-L fractional integral operator

(2.1)(2.2)

where denotes the integer derivative of order n, and is the R-L fractional integral of order

2.1 Modified fractional derivative with Mittage-Leffler kernel of RL-type

Using the results in Eqs (0.7) and (2.1) we arrive at the following infinite series representation of the modified fractional derivative with Mittage-Leffler kernel of RL-type

(2.3)

and for

(2.4)

2.2 Modified fractional derivative with Mittage-Leffler kernel of Caputo type

For the infinite series representation of the modified derivative of Caputo type was derived in [7]. For arbitrary we have

(2.5)

Using the results in (0.7) and (2.2) it holds that

Because we arrive at

(2.6)

Substituting the above result in Eq 2.5 will lead to

(2.7)

As a particular case of the above result and for and we have

(2.8)

which agrees with the result obtained in [7].

From the representation in Eqs 2.3 and (2.7) we arrive at the following relationship between the modified derivatives and .

Proposition 2.1. For and it holds that

(2.9)

As a particular case, and for we have

(2.10)

Because

(2.11)

we arrive at

(2.12)

Example 2.1. For the constant function u(t) = c₀, using the representation in Eq 2.3 we have

(2.13)

3 Dirac delta function

We recall that the Dirac delta function has the following form, namely

We remark here that there is no such real valued function with these properties, and the Dirac delta function is defined in the sense of distribution, or it can be captured by defining a special measure called the Dirac measure. The following properties hold true for the Dirac delta function

The above properties hold true for any interval where The R-L and Caputo derivatives of the Dirac delta were derived recently in [24, 25] and related FDEs were studied. The modified fractional derivative with Mittage-Leffler kernel of Caputo type for the Dirac delta function is given by

Because and

where H is the Heaviside function, then

(3.1)

Because then

(3.2)

Eq 3.1 leads to the fact that is a solution to the fractional initial value problem

Using the result in Eq 0.5 we arrive at

(3.3)

Proposition 3.1. For and , the solution of the FDE

(3.4)

is given by

(3.5)

where and provided that is well defined.

Proof: Because

applying the Laplace transform to (3.4) yields

Direct calculations will lead to

(3.6)

Applying the inverse Laplace transform we have

(3.7)

which completes the proof.

Proposition 3.2. For and , the solution of the FDE

(3.8)

is given by

(3.9)

where and provided that and is well defined.

Proof: Applying the Laplace transform to (3.8) will lead to

(3.10)

Thus,

(3.11)

Because the result follows by applying the inverse Laplace transform to (3.11).

Corollary 3.1. For and , the solution of the FDE

(3.12)

is given by

(3.13)

where

Proof: From Eq 2.12 we have − , and the results follow by substituting the above result in Eq 3.4.

Given that where H is the Heaviside function, then the solution of

(3.14)

is given by

(3.15)

where,

(3.16)

4 Numerical simulation

We consider the equation of the Resistor-Inductor circuit in the fractional case of the form

(4.1)

subject to the initial conditions . Here I(t) represents the current flowing through the circuit, − is the Dirac delta function input voltage, L = 1 is the inductance, and R = 2 is the resistance. The solution of the above system is given by Eq 3.13. The solution I(t) within the range and is depicted in Fig 1. Additionally, Fig 2 illustrates the solutions for in increments of 0.1. Figs 3 and 4 present the solutions with for in steps of 0.02, and for in steps of 0.1 respectively. The curves intersect at t = 0.571017 and t = 1.94714, which enriches the dynamics of the system by varying the order of the fractional derivative. All figures are plotted in the range of t between 0 and 3, and using the first 1000 terms of the infinite sum. We remark here that if we replace the modified derivative of Caputo type with the original ABC-derivative, then the problem with the initial condition admits no solutions.

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

Memory effects are introduced into the system by the parameter , which controls the order of the fractional derivative and affects the damping and transient responsiveness of the current I(t). While lower values () show slower decay and stronger memory effects, higher values () correlate to conventional integer-order dynamics with faster decay. The initial condition I(0) is also important; when , the current is driven only by the external input, but when , an interaction between the external stimulus and the initial stored energy is added, which enhances the dynamics of the system. Furthermore, as the dynamics converges across various values of , the curve intersections at t = 0.571017 and t = 1.94714 demonstrate universal behavior and indicate crucial points of balance within the circuit.

5 Conclusion and further work

We have introduced the modified derivative with Mittage-Leffler kernel of R-L type, and developed the basic theory of related fractional differential equations. Infinite series representations of the modified derivatives of R-L and Caputo types were derived and implemented to obtain a closed formula for the relationship between the two derivatives. The modified derivative of the Dirac delta function is presented and related resistor-inductor model is discussed. The solutions of the presented model enrich the dynamics of the system and indicate the efficiency of implementing the modified derivative in modeling real life problems. In one side, the problem admits solutions without the need of imposing extra conditions, and as the fractional derivative approaches 1, the solution coincides with the solution obtained by solving the associated differential equation with integer derivative of order 1. In future work, we aim to implement the modified derivative in modeling various dynamical systems and to develop suitable numerical schemes for integrating these systems.

References

1. 1. Caputo M, Fabrizio M. A new definition of fractional derivative without singular kernel. Progr Fract Differ Appl. 2015;1(2):73–85.

* View Article

* Google Scholar

2. 2. Atangana A, Baleanu D. New fractional derivatives with non-local and non- singular kernel: theory and application to heat transfer model. Thermal Sci. 2016;20(2):763–9.

* View Article

* Google Scholar

3. 3. Al-Refai M, Abdeljawad T. Analysis of the fractional diffusion equations with fractional derivative of non-singular kernel. Adv Diff Equ. 2017;:315.

4. 4. Al-Refai M. Fractional differential equations involving Caputo fractional derivative with Mittag-Leffler non-singular kernel: Comparison principles and applications. Electron J Differential Equations. 2018;36:1–10.

* View Article

* Google Scholar

5. 5. Al-Refai M. Fundamental results on systems of fractional differential equations involving Caputo-Fabrizio fractional derivative. Jordan J Math Stat. 2020;13(3):389–99.

* View Article

* Google Scholar

6. 6. Polyanin AD, Manzhirov AV. Handbook of integral equations. Second ed. Chapman and Hall/CRC. 2008.

7. 7. Al-Refai M, Baleanu D. On an extension of the operator with Mittag-Leffler kernel. Fractals. 2022;30(5):2240129.

* View Article

* Google Scholar

8. 8. Farman M, Jamil S, Riaz MB, Azeem M, Saleem M. Numerical and quantitative analysis of HIV/AIDS model with modified Atangana-Baleanu in Caputo sense derivative. Alex Eng J. 2023;66:31–42.

* View Article

* Google Scholar

9. 9. Khan H, Alzabut J, Alfwzan WF, Gulzar H. Nonlinear dynamics of a piecewise modified ABC fractional-order leukemia model with symmetric numerical simulations. Symmetry. 2023;15(1):Article 1338.

* View Article

* Google Scholar

10. 10. Taneja K, Deswal K, Kumar D, Baleanu D. Novel Numerical Approach for Time Fractional Equations with Nonlocal Condition. Numer Algor. 2023;95(3):1413–33.

* View Article

* Google Scholar

11. 11. Aldwoah KA, Almalahi MA, Hleili M, Alqarni FA, Aly ES, Shah K. Analytical study of a modified-ABC fractional order breast cancer model. J Appl Math Comput. 2024;70(1):3685–716.

* View Article

* Google Scholar

12. 12. Mehmet M, Arfan M, Sami A. Theoretical and numerical investigation of a modified ABC fractional operator for the spread of polio under the effect of vaccination. AIMS Biophys. 2024;11(1):97–120.

* View Article

* Google Scholar

13. 13. Zehra A, Jamil S, Farman M, Nisar KS. Modeling and analysis of Hepatitis B dynamics with vaccination and treatment with novel fractional derivative. PLoS One. 2024;19(7):e0307388. pmid:39024307

* View Article

* PubMed/NCBI

* Google Scholar

14. 14. Khan H, Alzabut J, Gómez-Aguilar JF, Alfwzan WF. A nonlinear perturbed coupled system with an application to chaos attractor. Results Phys. 2023;52:106891.

* View Article

* Google Scholar

15. 15. Partohaghighi M, Mortezaee M, Akgül A, Eldin SM. Numerical estimation of the fractional advection–dispersion equation under the modified Atanagana–Baleanu–Caputo derivative. Results Phys. 2023;49(106451).

* View Article

* Google Scholar

16. 16. Al-Refai M, Syam MI, Baleanu D. Analytical treatments to systems of fractional differential equations with modified Atangana-Baleanu derivative. Fractals. 2023;31(10):1–12.

* View Article

* Google Scholar

17. 17. Odibat Z. A new fractional derivative operator with a generalized exponential kernel. Nonlinear Dyn. 2024;112(17):15219–30.

* View Article

* Google Scholar

18. 18. Odibat Z, Al-Refai M, Baleanu D. On some properties of generalized cardinal sine kernel fractional operators: Advantages and applications of the extended operators. Chinese Journal of Physics. 2024;91:349–60.

* View Article

* Google Scholar

19. 19. Al-Refai M. Proper inverse operators of fractional derivatives with nonsingular kernels. Rend Circ Mat Palermo, II Ser. 2021;71(2):525–35.

* View Article

* Google Scholar

20. 20. Kamran A, Ahmadian A, Salahshour A, Salimi M. Robust numerical approximation of advection diffusion equations with nonsingular kernel derivative. Physica Scripta. 2021;96(12):ID 124015.

* View Article

* Google Scholar

21. 21. Liu X, Karman, Yao Y. Numerical approximation of Riccati fractional differential equation in the sense of Caputo-type fractional derivative. J Math. 2020.

22. 22. Kilbas AA, Srivastava HM, Trujillo JJ. Theory and applications of fractional differential equations. Amsterdam: Elsevier Science B.V.. 2006.

23. 23. Podlubny I. Fractional differential equations. San Diego: Academic Press. 1999.

24. 24. Camrud E. A novel approach to fractional calculus: Utilizing fractional integrals and derivatives of the Dirac delta function. Progr Fract Differ Appl. 2018;4(4):463–78.

* View Article

* Google Scholar

25. 25. Feng Z, Ye L, Zhang Y. On the fractional derivative of Dirac delta function and its application. Adv Math Phys. 2020.

Citation: Al-Refai M, Baleanu D, Alomari A (2025) On the solutions of fractional differential equations with modified Mittag-Leffler kernel and Dirac Delta function: Analytical results and numerical simulations. PLoS One 20(6): e0325897. https://doi.org/10.1371/journal.pone.0325897

About the Authors:

Mohammed Al-Refai

Roles: Conceptualization, Formal analysis, Methodology, Project administration, Writing – original draft, Writing – review & editing

E-mail: [email protected]

Affiliation: Department of Mathematics, Yarmouk University, Irbed, Jordan

ORICD: https://orcid.org/0000-0001-9399-6756

Dumitru Baleanu

Roles: Conceptualization, Funding acquisition, Investigation, Supervision

Affiliations: Department of Computer Science and Mathematics, Lebanese American University, Beirut, Lebanon, Institute of Space Sciences, Magurele-Bucharest, Romania

A.K. Alomari

Roles: Investigation, Software, Writing – review & editing

Affiliation: Department of Mathematics, Faculty of Science, Islamic University of Madinah, Madinah, Saudi Arabia

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

References

1. Caputo M, Fabrizio M. A new definition of fractional derivative without singular kernel. Progr Fract Differ Appl. 2015;1(2):73–85.

2. Atangana A, Baleanu D. New fractional derivatives with non-local and non- singular kernel: theory and application to heat transfer model. Thermal Sci. 2016;20(2):763–9.

3. Al-Refai M, Abdeljawad T. Analysis of the fractional diffusion equations with fractional derivative of non-singular kernel. Adv Diff Equ. 2017;:315.

4. Al-Refai M. Fractional differential equations involving Caputo fractional derivative with Mittag-Leffler non-singular kernel: Comparison principles and applications. Electron J Differential Equations. 2018;36:1–10.

5. Al-Refai M. Fundamental results on systems of fractional differential equations involving Caputo-Fabrizio fractional derivative. Jordan J Math Stat. 2020;13(3):389–99.

6. Polyanin AD, Manzhirov AV. Handbook of integral equations. Second ed. Chapman and Hall/CRC. 2008.

7. Al-Refai M, Baleanu D. On an extension of the operator with Mittag-Leffler kernel. Fractals. 2022;30(5):2240129.

8. Farman M, Jamil S, Riaz MB, Azeem M, Saleem M. Numerical and quantitative analysis of HIV/AIDS model with modified Atangana-Baleanu in Caputo sense derivative. Alex Eng J. 2023;66:31–42.

9. Khan H, Alzabut J, Alfwzan WF, Gulzar H. Nonlinear dynamics of a piecewise modified ABC fractional-order leukemia model with symmetric numerical simulations. Symmetry. 2023;15(1):Article 1338.

10. Taneja K, Deswal K, Kumar D, Baleanu D. Novel Numerical Approach for Time Fractional Equations with Nonlocal Condition. Numer Algor. 2023;95(3):1413–33.

11. Aldwoah KA, Almalahi MA, Hleili M, Alqarni FA, Aly ES, Shah K. Analytical study of a modified-ABC fractional order breast cancer model. J Appl Math Comput. 2024;70(1):3685–716.

12. Mehmet M, Arfan M, Sami A. Theoretical and numerical investigation of a modified ABC fractional operator for the spread of polio under the effect of vaccination. AIMS Biophys. 2024;11(1):97–120.

13. Zehra A, Jamil S, Farman M, Nisar KS. Modeling and analysis of Hepatitis B dynamics with vaccination and treatment with novel fractional derivative. PLoS One. 2024;19(7):e0307388. pmid:39024307

14. Khan H, Alzabut J, Gómez-Aguilar JF, Alfwzan WF. A nonlinear perturbed coupled system with an application to chaos attractor. Results Phys. 2023;52:106891.

15. Partohaghighi M, Mortezaee M, Akgül A, Eldin SM. Numerical estimation of the fractional advection–dispersion equation under the modified Atanagana–Baleanu–Caputo derivative. Results Phys. 2023;49(106451).

16. Al-Refai M, Syam MI, Baleanu D. Analytical treatments to systems of fractional differential equations with modified Atangana-Baleanu derivative. Fractals. 2023;31(10):1–12.

17. Odibat Z. A new fractional derivative operator with a generalized exponential kernel. Nonlinear Dyn. 2024;112(17):15219–30.

18. Odibat Z, Al-Refai M, Baleanu D. On some properties of generalized cardinal sine kernel fractional operators: Advantages and applications of the extended operators. Chinese Journal of Physics. 2024;91:349–60.

19. Al-Refai M. Proper inverse operators of fractional derivatives with nonsingular kernels. Rend Circ Mat Palermo, II Ser. 2021;71(2):525–35.

20. Kamran A, Ahmadian A, Salahshour A, Salimi M. Robust numerical approximation of advection diffusion equations with nonsingular kernel derivative. Physica Scripta. 2021;96(12):ID 124015.

21. Liu X, Karman, Yao Y. Numerical approximation of Riccati fractional differential equation in the sense of Caputo-type fractional derivative. J Math. 2020.

22. Kilbas AA, Srivastava HM, Trujillo JJ. Theory and applications of fractional differential equations. Amsterdam: Elsevier Science B.V.. 2006.

23. Podlubny I. Fractional differential equations. San Diego: Academic Press. 1999.

24. Camrud E. A novel approach to fractional calculus: Utilizing fractional integrals and derivatives of the Dirac delta function. Progr Fract Differ Appl. 2018;4(4):463–78.

25. Feng Z, Ye L, Zhang Y. On the fractional derivative of Dirac delta function and its application. Adv Math Phys. 2020.