Using Lagrangian descriptors to reveal the phase space structure of dynamical systems described by fractional differential equations: Application to the Duffing oscillator

TL;DR

This paper demonstrates how Lagrangian descriptors can be used to analyze the phase space of fractional-order dynamical systems, specifically applying it to the Duffing oscillator with different fractional derivatives and backward-time integration methods.

Contribution

It introduces a novel application of Lagrangian descriptors to fractional differential equations and compares backward-time integration techniques for revealing phase space structures.

Findings

Nonlocal implicit inverse method effectively reveals phase space structure.

Differences in backward-time integration methods impact Lagrangian descriptor results.

The study provides insights into the physical significance of these differences.

Abstract

We showcase the utility of the Lagrangian descriptors method in qualitatively understanding the underlying dynamical behavior of dynamical systems governed by fractional-order differential equations. In particular, we use the Lagrangian descriptors method to study the phase space structure of the unforced and undamped Duffing oscillator when its time evolution is governed by fractional-order differential equations. In our study, we implement two types of fractional derivatives, namely the standard Gr\"unwald-Letnikov method, which is a finite difference approximation of the Riemann-Liouville fractional derivative, and a Gr\"unwald-Letnikov method with a correction term that approximates the Caputo fractional derivative. While there is no issue with forward-time integrations needed for the evaluation of Lagrangian descriptors, we discuss in detail ways to perform the non-trivial task of backward-time integrations and implement two methods for this purpose: a `nonlocal implicit inverse' technique and a `time-reverse inverse' approach. We analyze the differences in the Lagrangian descriptors results due to the two backward-time integration approaches, discuss the physical significance of these differences, and eventually argue that the nonlocal implicit inverse implementation of the Gr\"unwald-Letnikov fractional derivative manages to reveal the phase space structure of fractional-order dynamical systems correctly.