Recent numerical works have revealed the instability of many-body localized (MBL) phase in disordered quantum many-body systems with finite system sizes and over finite timescales. This instability originates from the rare Griffiths regions in the thermodynamic limit, which rapidly thermalize and induce an avalanche mechanism that destabilizes surrounding MBL regions. Here, we consider the \(\mathbb {Z}_2\) -preserving interacting Ising Majorana chain, whose phase diagram is enriched by the presence of an intermediate ergodic phase separating two distinct MBL phases characterized by different long-range orders. We investigate the dynamic characteristics of the model when coupled to an infinite bath under perturbation, and through the scaling behavior of the slowest thermalization rate, we find how the critical disorder strength in finite-size systems is affected by the avalanche mechanism. We also implement the embedded inclusion model and use the time evolution of mutual information between each spin and the artificial Griffiths region to probe the diffusion of the thermal bubble. We observe that in finite-sized systems, the critical disorder strength gradually drifts away from the center. Our results show that both the MBL paramagnetic phase and MBL spin-glass phase are unstable at finite sizes.
