Performance of chaos diagnostics based on Lagrangian descriptors. Application to the 4D standard map

Abstract: We investigate the ability of simple diagnostics based on Lagrangian descriptor (LD) computations of initially nearby orbits to detect chaos in conservative dynamical systems with phase space dimensionality higher than two. In particular, we consider the recently introduced methods of the difference ($D_L^n$) and the ratio ($R_L^n$) of the LDs of neighboring orbits, as well as a quantity ($S_L^n$) related to the finite-difference second spatial derivative of the LDs, and use them to determine the chaotic or regular nature of ensembles of orbits of a prototypical area-preserving map model, the 4-dimensional (4D) symplectic standard map. Using the distributions of the indices' values we determine appropriate thresholds to discriminate between regular and chaotic orbits, and compare the obtained characterization against that achieved by the Smaller Alignment Index (SALI) method of chaos detection, by recording the percentage agreement $P_A$ between the two classifications. We study the influence of various factors on the performance of these indices, and show that the increase of the final number of orbit iterations T and the order n of the indices (i.e. the dimensionality of the space where the considered nearby orbits lie), as well as the decrease of the distance $\sigma$ of neighboring orbits, increase the $P_A$ values along with the required computational effort. Balancing between these two factors we find appropriate T, n and $\sigma$ values, which allow the efficient use of the $D_L^n$, $R_L^n$ and $S_L^n$ indices as short time and computationally cheap chaos diagnostics achieving $P_A \gtrsim 90 \%$, with $D_L^n$ and $S_L^n$ having larger $P_A$ values than $R_L^n$. Our results show that the three LDs-based indices perform better for systems with large percentages of chaotic orbits.