The study of oscillations of relativistic stars is motivated by the prospect of detecting such oscillations in electromagnetic or gravitational wave signals. In the same way that helioseismology is providing us with information about the interior of the Sun, the observational identification of oscillation frequencies of relativistic stars could constrain the high-density equation of state . The oscillations could be excited after a core collapse or during the final stages of a neutron star binary merger. Rapidly rotating relativistic stars can become unstable to the emission of gravitational waves.
When the oscillations of an equilibrium star are of small magnitude compared to its radius, it will suffice to approximate them as linear perturbations. Such perturbations can be described in two equivalent ways. In the Lagrangian approach, one studies the changes in a given fluid element as it oscillates about its equilibrium position. In the Eulerian approach, one studies the change in fluid variables at a fixed point in space. Both approaches have their strengths and weaknesses.
In the Newtonian limit, the Lagrangian approach has been used to develop variational principles [215, 118], but the Eulerian approach proved to be more suitable for numerical computations of mode frequencies and eigenfunctions [162, 217, 158, 160, 159]. Clement  used the Lagrangian approach to obtain axisymmetric normal modes of rotating stars, while nonaxisymmetric solutions were obtained in the Lagrangian approach by Imamura et al.  and in the Eulerian approach by Managan  and Ipser and Lindblom . While a lot has been learned from Newtonian studies, in the following we will focus on the relativistic treatment of oscillations of rotating stars.
© Max Planck Society and the author(s)