Explanation of the ATLAS Z-peaked excess in the NMSSM

Recently the ATLAS collaboration reported a $3\sigma$ excess in the leptonic-$Z+jets+E_{T}^{miss}$ channel. This may be interpreted in the Next-to-Minimal Supersymmetric Standard Model (NMSSM) by gluino pair production with the decay chain $\tilde{g} \to q \bar{q} \tilde{\chi}_2^0 \to q \bar{q} Z \tilde{\chi}_1^0$, where $\tilde{\chi}_1^0$ and $\tilde{\chi}_2^0$ denote the lightest and the next-to-lightest neutralinos with singlino and bino as their dominant components respectively. After exploring the relevant parameter space of the NMSSM by considering the constraints from the ATLAS searches for $jets + E_{T}^{miss}$ signals, we conclude that the NMSSM is able to explain the excess at $1 \sigma$ level with the number of the signal events reaching its measured central value in optimal cases, and the best explanation comes from a compressed spectrum such as $m_{\tilde{g}} \simeq 650 {\rm GeV}$, $m_{\tilde{\chi}_2^0} \simeq 565 {\rm GeV}$ and $m_{\tilde{\chi}_1^0} \simeq 465 {\rm GeV}$. We also check the consistency of the ATLAS results with the null result of the CMS on-$Z$ search. We find that under the CMS limits at $95\%$ C.L., the event number of the ATLAS on-$Z$ signal can still reach 11 in our scenario, which is about $1.2 \sigma$ away from the measured central value.