A method for Lyapunov spectrum estimation using cloned dynamics and its application to the discontinuously-excited FitzHugh-Nagumo model

This work presents a new method to calculate the Lyapunov spectrum of dynamical systems based on the time evolution of initially small disturbed copies (clones) of the motion equations. In this approach, it is not necessary to construct the tangent space associated with the time evolution of linearized versions of motion equations, being the Lyapunov exponents directly estimated in terms of the rate of convergence or divergence of these disturbed clones with respect to the fiducial trajectory, there being periodic correction via the Gram-Schmidt Reorthonormalization procedure. The proposed method offers the possibility of partial estimation of the Lyapunov spectrum and can also be applied to nonsmooth dynamics, since the linearization procedure is no longer required. The idea is tested for representative continuous- and discrete-time dynamical systems and validated by means of comparison with the classical method to perform this calculation. To illustrate its applicability in the nonsmooth context, the largest Lyapunov exponent of the FitzHugh-Nagumo neuronal model under discontinuous periodic excitation is calculated taking the amplitude of stimulation as control parameter. This analysis reveals some complex behaviours for this simple neuronal model, which motivates relevant discussions about the possible role of chaos in the cognitive process.