Asymptotic theory of gravity modes in rotating stars I. Ray dynamics

Context. The seismology of early-type stars is limited by our incomplete understanding of gravito-inertial modes. Aims. We develop a short-wavelength asymptotic analysis for gravito-inertial modes in rotating stars. Methods. The Wentzel-Kramers-Brillouin approximation was applied to the equations governing adiabatic small perturbations about a model of a uniformly rotating barotropic star. Results. A general eikonal equation, including the effect of the centrifugal deformation, is derived. The dynamics of axisymmetric gravito-inertial rays is solved numerically for polytropic stellar models of increasing rotation and analysed by describing the structure of the phase space. Three different types of phase-space structures are distinguished. The first type results from the continuous evolution of structures of the non-rotating integrable phase space. It is predominant in the low-frequency region of the phase space. The second type of structures are island chains associated with stable periodic rays. The third type of structures are large chaotic regions that can be related to the envelope minimum of the Brunt-Vaisala frequency. Conclusions. Gravito-inertial modes are expected to follow this classification, in which the frequency spectrum is a superposition of sub-spectra associated with these different types of phase-space structures. The detailed confrontation between the predictions of this ray-based asymptotic theory and numerically computed modes will be presented in a companion paper.