CP Prediction from Residual $\mathbb Z_2^s$ and $\overline{\mathbb Z}_2^s$ Symmetries with JUNO First Data

The JUNO first data and the recent neutrino global fit results are implemented in the sum rule from the residual $\mathbb Z^s_2$ and $\overline{\mathbb Z}^s_2$ symmetries to make prediction of the leptonic Dirac CP phase $\delta_D$. Without involving model parameters, the probability distribution of $\delta_D$ can be readily obtained from the experimental measurements of the three mixing angles. We then confront the theoretical predictions with the global fit results for the CP phase as well as the T2K and NOvA joint analysis for their CP measurement to give the data preference of the two residual symmetries with Bayes factor for both normal and inverted orderings. We further extend our analysis to a two-dimensional probability distribution to fully explore the correlation between the CP phase $\delta_D$ and the atmospheric angle $\theta_a \equiv \theta_{23}$.