We explore the chaotic dynamics of the mass deformed ABJM model. To do so, we first perform a dimensional reduction of this model from 2 + 1-to 0 + 1dimensions, considering that the fields are spatially uniform. Working in the 't Hooft limit and tracing over ansatz configurations involving fuzzy two spheres, which are described in terms of the GRVV matrices with collective time dependence, we obtain a family of reduced effective Lagrangians and demonstrate that they have chaotic dynamics by computing the associated Lyapunov spectrum. In particular, we analyze in detail, how the largest Lyapunov exponent, λ L , changes as a function of E/N 2 . Depending on the structure of the effective potentials, we find either, where γ N (k, µ) is a constant determined in terms of the Chern-Simons coupling k, the mass µ, and the matrix level N . Using our results, we investigate the temperature dependence of the largest Lyapunov exponents and give upper bounds on the temperature above which λ L values comply with the MSS bound, λ L ≤ 2πT , and below which it will eventually be violated.

Starting from an 𝑆𝑈 (𝑁) matrix quantum mechanics model with massive deformation terms and by introducing an ansatz configuration involving fuzzy four-and two-spheres with collective time dependence, we obtain a family of effective Hamiltonians, 𝐻 𝑛 , (𝑁 = 1 6 (𝑛 + 1) (𝑛 + 2) (𝑛 + 3)) and examine their emerging chaotic dynamics. Through numerical work, we model the variation of the largest Lyapunov exponents as a function of the energy and find that they vary either as ∝ (𝐸 − (𝐸 𝑛 ) 𝐹 ) 1/4 or ∝ 𝐸 1/4 , where (𝐸 𝑛 ) 𝐹 stand for the energies of the unstable fixed points of the phase space. We use our results to put upper bounds on the temperature above which the Lyapunov exponents comply with the Maldacena-Shenker-Stanford (MSS) bound, 2𝜋𝑇, and below which it will eventually be violated.