Polynomial mixing of the critical Ising model on sparse Erdos-Renyi graphs

We consider the stochastic Ising model on sparse Erdos-Renyi graphs $G(n,d/n)$ with $d>1$ at the critical temperature $\beta_c=\tanh^{-1}(d^{-1})$ and prove that with high probability, the mixing time is at most polynomial in $n$. Our approach combines the recent stochastic localization framework of Chen and Eldan, which yields spectral gap bounds in the well-behaved bulk of the graph, together with classical results on the relaxation time of Glauber dynamics on trees to handle regions where we cannot apply the Chen-Eldan method directly because of atypically large local neighborhoods.