Extend KEH framework to differential rotation
Develop a fully general relativistic computational framework within the Komatsu–Eriguchi–Hachisu (KEH) method to model axisymmetric neutron stars with differential rotation by solving the angular momentum relation j(Ω) = A^2(Ω_c − Ω) (Eq. 39) self-consistently, rather than imposing the rigid-rotation limit Ω = Ω_c, to obtain equilibrium configurations for finite A.
References
To correctly describe differential rotation in a fully general relativistic framework, Eq.~39 should be solved through a self-consistent iterative procedure. In the present study, however, we restrict ourselves to the limiting case of large $A$, which corresponds to rigid-body rotation, i.e., $\Omega = \Omega_c$. We will focus on exploring the relationship between rotational effects and properties of the nuclear equation of state, leaving the extension to differential rotation for future work.