In this paper, we investigate massive charged scalar perturbations in four-dimensional charged Lifshitz–AdS black holes with scalar hair within the framework of Einstein–Maxwell–Dilaton (EMD) gravity. Using the improved asymptotic iteration method (AIM), we compute the quasinormal modes (QNMs) and explore their dependence on key parameters, including the Lifshitz dynamical exponent z, the scalar field mass and charge, and the black hole charge, under various spatial curvature settings (k=0,±1). Our results reveal rich and sensitive behavior in both the real and imaginary parts of the QNMs. In particular, the decay rates can exhibit monotonic or non-monotonic dependence on the black hole charge, depending on the values of z, ms, and qs. These findings highlight the significant role of field and geometric parameters in governing the dynamical stability of Lifshitz black holes and offer insights into the perturbative properties of non-AdS holographic systems.
Zhang et al. (Sat,) studied this question.