Astrophysical charge-free black holes are known to satisfy no-hair relations through which all multipole moments can be specified in terms of just their mass and spin angular momentum. We here investigate the possible existence of no-hair-like relations among multipole moments for neutron stars and quark stars that are independent of their equation of state. We calculate the multipole moments of these stars up to hexadecapole order by constructing uniformly-rotating and unmagnetized stellar solutions to the Einstein equations. For slowly-rotating stars, we construct stellar solutions to quartic order in spin in a slow-rotation expansion, while for rapidly-rotating stars, we solve the Einstein equations numerically with the LORENE and RNS codes. We find that the multipole moments extracted from these numerical solutions are consistent with each other. We confirm that the current-dipole is related to the mass-quadrupole in an approximately equation of state independent fashion, which does not break for rapidly rotating neutron stars or quark stars. We further find that the current-octupole and the mass-hexadecapole moments are related to the mass-quadrupole in an approximately equation of state independent way to $\sim 10%$, worsening in the hexadecapole case. All of our findings are in good agreement with previous work that considered stellar solutions to leading-order in a weak-field expansion. The quartic in spin, slowly-rotating solutions found here allow us to estimate the systematic errors in the measurement of the neutron star's mass and radius with future X-ray observations, such as NICER and LOFT. We find that the effect of these quartic-in-spin terms on the quadrupole and hexadecapole moments and stellar eccentricity may dominate the error budget for very rapidly-rotating neutron stars. The new universal relations found here should help to reduce such systematic errors.

Black holes are said to have no hair because all of their multipole moments can be expressed in terms of just their mass, charge and spin angular momentum. The recent discovery of approximately equation-of-state-independent relations among certain multipole moments in neutron stars suggests that they are also approximately bald. We here explore the yet unknown origin for this universality. First, we investigate which region of the neutron star's interior and of the equation of state is most responsible for the universality. We find that the universal relation between the moment of inertia and the quadrupole moment is dominated by the star's outer-core, a shell of width (50-95)% of the total radius, which corresponds to the density range $(10^{14}$-$10^{15})$g/cm$^3$. Second, we study the impact on the universality of approximating stellar isodensity contours as self-similar ellipsoids. An analytical calculation in the Newtonian limit reveals that the shape of the ellipsoids does not affect the universal relations, but relaxing the self-similarity assumption can completely destroy it. Third, we investigate the eccentricity profiles of rotating relativistic stars and find that the ellipticity is roughly constant, with variations of roughly (20-30)% in the region that matters to the universal relations. Fourth, we repeat the above analysis for non-compact, regular stars and find that the ellipticity is not constant, with variations that easily exceed 100% and universality is lost. These findings suggest that universality arises as an emergent approximate symmetry: as one flows in the stellar-structure phase space from non-compact star region to the relativistic star region, the eccentricity variation inside stars decreases, leading to approximate self-similarity in their isodensity contours, which then leads to the universal behavior observed in their exterior multipole moments.

Recently it was shown that slowly rotating neutron stars exhibit an interesting correlation between their moment of inertia $I$, their quadrupole moment $Q$, and their tidal deformation Love number $\lambda$ (the I-Love-Q relations), independently of the equation of state of the compact object. In the present work a similar, more general, universality is shown to hold true for all rotating neutron stars within General Relativity; the first four multipole moments of the neutron star are related in a way independent of the nuclear matter equation of state we assume. By exploiting this relation, we can describe quite accurately the geometry around a neutron star with fewer parameters, even if we don't know precisely the equation of state. Furthermore, this universal behavior displayed by neutron stars, could promote them to a more promising class of candidates (next to black holes) for testing theories of gravity.

In this article we present some recent results on identifying correctly the relativistic multipole moments of numerically constructed spacetimes, and the consequences that this correction has on searching for appropriate analytic spacetimes that can approximate well the previously mentioned numerical spacetimes. We also present expressions that give the quadrupole and the spin octupole as functions of the spin parameter of a neutron star for various equations of state and in a range of masses for every equation of state used. These results are relevant for describing the exterior spacetime of rotating neutron stars that are made up of matter obeying realistic equations of state.

We have tested the appropriateness of two-soliton analytic metric to describe the exterior of all types of neutron stars, no matter what their equation of state or rotation rate is. The particular analytic solution of the vaccuum Einstein equations proved quite adjustable to mimic the metric functions of all numerically constructed neutron-star models that we used as a testbed. The neutron-star models covered a wide range of stiffness, with regard to the equation of state of their interior, and all rotation rates up to the maximum possible rotation rate allowed for each such star. Apart of the metric functions themselves, we have compared the radius of the innermost stable circular orbit $R_{\rm{ISCO}}$, the orbital frequency $\Omega\equiv\frac{d\phi}{dt}$ of circular geodesics, and their epicyclic frequencies $\Omega_{\rho}, \Omega_z$, as well as the change of the energy of circular orbits per logarithmic change of orbital frequency $\Delta\tilde{E}$. All these quantities, calculated by means of the two-soliton analytic metric, fitted with good accuracy the corresponding numerical ones as in previous analogous comparisons (although previous attempts were restricted to neutron star models with either high or low rotation rates). We believe that this particular analytic solution could be considered as an analytic faithful representation of the gravitation field of any rotating neutron star with such accuracy, that one could explore the interior structure of a neutron star by using this space-time to interpret observations of astrophysical processes that take place around it.

Identifying the relativistic multipole moments of a spacetime of an astrophysical object that has been constructed numerically is of major interest, both because the multipole moments are intimately related to the internal structure of the object, and because the construction of a suitable analytic metric that mimics a numerical metric should be based on the multipole moments of the latter one, in order to yield a reliable representation. In this note we show that there has been a widespread delusion in the way the multipole moments of a numerical metric are read from the asymptotic expansion of the metric functions. We show how one should read correctly the first few multipole moments (starting from the quadrupole mass-moment), and how these corrected moments improve the efficiency of describing the metric functions with analytic metrics that have already been used in the literature, as well as other consequences of using the correct moments.

We study small-amplitude, nonlinear pulsations of uniformly and differentially rotating neutron stars employing a two-dimensional evolution code for general-relativistic hydrodynamics. Using Fourier transforms at several points inside the star, both the eigenfrequencies and two-dimensional eigenfunctions of pulsations are extracted. The centrifugal forces and the degree of differential rotation have significant effects on the mode-eigenfunction. We find that near the mass-shedding limit, the pulsations are damped due to shocks forming at the surface of the star. This new damping mechanism may set a small saturation amplitude for modes that are unstable to the emission of gravitational waves. After correcting for the assumption of the Cowling approximation (used in our numerical code), we construct empirical relations that predict the range of gravitational-wave frequencies from quasi-periodic post-bounce oscillations in the core collapse of massive stars. We also find that the fundamental quasi-radial mode is split, at least in the Cowling approximation and mainly in differentially rotating stars, into two different sequences.

We study the problem of detecting, and infering astrophysical information from, gravitational waves from a pulsating neutron star. We show that the fluid f and p-modes, as well as the gravitational-wave w-modes may be detectable from sources in our own galaxy, and investigate how accurately the frequencies and damping rates of these modes can be infered from a noisy gravitational-wave data stream. Based on the conclusions of this discussion we propose a strategy for revealing the supranuclear equation of state using the neutron star fingerprints: the observed frequencies of an f and a p-mode. We also discuss how well the source can be located in the sky using observations with several detectors.

Gravitational-wave interferometers are expected to monitor the last three minutes of inspiral and final coalescence of neutron star and black hole binaries at distances approaching cosmological, where the event rate may be many per year. Because the binary's accumulated orbital phase can be measured to a fractional accuracy $\ll 10^{-3}$ and relativistic effects are large, the waveforms will be far more complex, carry more information, and be far harder to model theoretically than has been expected. Theorists must begin now to lay a foundation for extracting the waves' information.