Though individual stellar parameters of compact stars usually demonstrate obvious dependence on the equation of state (EOS), EOS-insensitive universal formulas relating these parameters remarkably exist. In the present paper, we explore the interrelationship between two such formulas, namely the $f$-$I$ relation connecting the $f$-mode quadrupole oscillation frequency $\omega_2$ and the moment of inertia $I$, and the $I$-Love-$Q$ relations relating $I$, the quadrupole tidal deformability $\lambda_2$, and the quadrupole moment $Q$, which have been proposed by Lau, Leung, and Lin [Astrophys. J. {\bf 714}, 1234 (2010)] and Yagi and Yunes [Science {\bf 341}, 365 (2013)], respectively. A relativistic universal relation between $\omega_l$ and $\lambda_l$ with the same angular momentum $l=2,3,\ldots$, the so-called "diagonal $f$-Love relation" that holds for realistic compact stars and stiff polytropic stars, is unveiled here. An in-depth investigation in the Newtonian limit is further carried out to pinpoint its underlying physical mechanism and hence leads to a unified $f$-$I$-Love relation. We reach the conclusion that these EOS-insensitive formulas stem from a common physical origin --- compact stars can be considered as quasiincompressible when they react to slow time variations introduced by $f$-mode oscillations, tidal forces and rotations.

The recent discovery of the universal I-Love-Q relations connecting the moment of inertia, tidal deformability, and the spin-induced quadrupole moment of compact stars is intriguing and totally unexpected. In this paper, we provide numerical evidence showing that the universality can be attributed to the incompressible limit of the I-Love-Q relations. The fact that modern equations of state are stiff, with an effective adiabatic index larger than about two, above the nuclear density range is the key to establishing the universality for neutron stars and quark stars with typical compactness from about 0.1 to 0.3. On the other hand, the I-Love-Q relations of low-mass neutron stars near the minimum mass limit depend more sensitively on the underlying equation of state because these stars are composed mainly of softer matter at low densities. However, the I-Love-Q relations for low-mass quark stars can still be represented accurately by the incompressible limit. We also study the I-Love relation connecting the moment of inertia and tidal deformability analytically in Newtonian gravity and show why the I-Love-Q relation is weakly dependent on the underlying equation of state and can be attributed to its incompressible limit.

Due to our ignorance of the equation of state (EOS) beyond nuclear density, there is still no unique theoretical model for neutron stars (NSs). It is therefore surprising that universal EOS-independent relations connecting different physical quantities of neutron stars can exist. Lau et al. [ApJ, 714, 1234 (2010)] found that the frequency of the $f$-mode oscillation, the mass, and the moment of inertia are connected by universal relations. More recently, Yagi and Yunes [Science, 341, 365 (2013)] discovered the I-Love-Q universal relations among the mass, the moment of inertia, the Love number, and the quadrupole moment. In this paper, we study these universal relations in the Eddington-inspired Born-Infeld (EiBI) gravity. This theory differs from general relativity (GR) significantly only at high densities due to the nonlinear coupling between matter and gravity. It thus provides us an ideal case to test how robust the universal relations of NSs are with respect to the change of the gravity theory. Thanks to the apparent EOS formulation of EiBI gravity developed recently by Delsate and Steinhoff [Phys. Rev. Lett., 109, 021101 (2012)], we are able to study the universal relations in EiBI gravity using the same techniques as those in GR. We find that the universal relations in EiBI gravity are essentially the same as those in GR. Our work shows that, within the currently viable coupling constant, there exists at least one modified gravity theory that is indistinguishable from GR in view of the unexpected universal relations.

We study how generic phase transitions taking place in compact stars constructed in the framework of the Eddington-inspired Born-Infeld (EiBI) gravity can lead to anomalous behavior of these stars. For the case with first-order phase transitions, compact stars in EiBI gravity with a positive coupling parameter $\kappa $ exhibit a finite region with constant pressure, which is absent in general relativity. However, for the case with a negative $\kappa $, an equilibrium stellar configuration cannot be constructed. Hence, EiBI gravity seems to impose stricter constraints on the microphysics of stellar matter. Besides, in the presence of spatial discontinuities in the sound speed $c_s$ due to phase transitions, the Ricci scalar is spatially discontinuous and contains $\delta$-function singularities proportional to the jump in $c_s^2$ acquired in the associated phase transition.

We study the hydrostatic equilibrium structure of compact stars in the Eddington-inspired Born-Infeld gravity recently proposed by Banados and Ferreira [Phys. Rev. Lett. 105, 011101 (2010)]. We also develop a framework to study the radial perturbations and stability of compact stars in this theory. We find that the standard results of stellar stability still hold in this theory. The frequency square of the fundamental oscillation mode vanishes for the maximum-mass stellar configuration. The dependence of the oscillation mode frequencies on the coupling parameter \kappa of the theory is also investigated. We find that the fundamental mode is insensitive to the value of \kappa, while higher order modes depend more strongly on \kappa.