Stability maps illustrating survival and ejection probabilities for exomoons orbiting a close-in giant planet, factoring in the planet’s J₂ quadrupole moment and relativistic corrections over 1000 P_p. The colour scale indicates the selenity, which shows the moon’s likelihood of remaining bound to the planet under external perturbations. The left panel shows stability in the a_s − e_s plane, with the Roche limit and Hill radius as key markers. The right panel displays stability in the a_s − M_s plane, reflecting stability variations based on satellite mass and planet distance.