Thermalization in massive deformations of Yang-Mills matrix models

We numerically study the classical evolution of a Yang-Mills matrix model with two distinct mass deformation terms, which can be contemplated as a massive deformation of the bosonic part of the BFSS model. Through numerical analysis, it is shown that when the simulations are started from a certain set of initial conditions, thermalization occurs. Besides, an estimation method is proposed to determine the approximate thermalization time. Using this method, we demonstrate that thermalization time vary logarithmically with increasing matrix size when the mass terms differ. Introducing a matrix configuration, we also obtain reduced actions and subsequently analyze how the thermalization time change as a function of the energy.