Fast Rotating Neutron Stars: Oscillations and Instabilities
- 1Theoretical Astrophysics, IAAT, University of Tübingen, Tübingen, Germany
- 2SISSA, INFN Sezione di Trieste, Trieste, Italy
- 3IFPU, Trieste, Italy
In this review article, we present the main results from our most recent research concerning the oscillations of fast rotating neutron stars. We derive a set of time evolution equations for the investigation of non-axisymmetric oscillations of rapidly rotating compact objects in full general relativity, taking into account the contribution of a dynamic spacetime. Using our code, which features high accuracy at comparably low computational expense, we are able to extract the frequencies of non-axisymmetric modes of compact objects with rotation rates up to the Kepler limit. We propose various universal relations combining bulk properties of isolated neutron stars as well as of binary systems before and after merger; these relations are independent of the true equation of state and may serve as a valuable tool for gravitational wave asteroseismology. We also present an introductory example using a Bayesian analysis.
Oscillations and instabilities of neutron stars were always considered among the promising sources for gravitational waves. The systematic study of non-axisymmetric neutron star oscillations began in the 1960s with the pioneering works of Thorne and collaborators (cf. Thorne and Campolattaro (1967); Thorne (1969) (and subsequent papers)), in which they laid out the equations governing the perturbations of compact stars in general relativity. The numerical solution of these equations has proven highly challenging and it has taken nearly two decades before Lindblom and Detweiler found an advantageous formulation of the eigenvalue problem that allowed them to determine the complex frequencies of the acoustic modes of sufficiently realistic stellar models (Lindblom and Detweiler (1983); Detweiler and Lindblom (1985)). These results did not conclude the investigation of these modes; in particular Chandrasekhar and Ferrari (1991a); Chandrasekhar and Ferrari (1991b) turned to the perturbations of relativistic stars and studied their oscillations as a scattering problem.
In the mid-1980s, inspired by a toy model Kokkotas and Schutz (1986) suggested that the dynamic spacetime of a neutron star exhibits its very own class of modes and christened them w-modes (Kokkotas and Schutz (1992)). The short damping times of these modes (comparable to the ones of black-holes) pose numerical challenges, but the application of the continued fraction method (Leaver (1985)) and numerical integration along anti-Stokes lines (Andersson et al. (1995)) improved considerably the accuracy in the calculation of both fluid and spacetime modes and led to the discovery of the new class of w-modes, the so-called wII-modes (Leins et al. (1993)). This advancement in computational accuracy allowed the study of many neutron star equations of state leading to the discovery of universal relations involving frequencies and damping times of oscillation modes. These universal relations provided ways for solving the inverse problem within the so-called gravitational wave asteroseismology, allowing to set constrains on the mass and/or radius of the neutron star, and eventually the nuclear equation of state, via the observation of oscillations (see, e.g., Andersson and Kokkotas (1996, 1998); Kokkotas et al. (2001); Benhar et al. (2004)).
With increasing computational power during the 1990s, the first attempts were made to increase the dimensionality of the hitherto (due to the restriction to spherical symmetry) purely radial problem and include time as a second dimension. The first successful time evolutions of perturbations of relativistic neutron stars were reported by Allen et al. (1998). Ruoff (2000, 2001) reformulated those evolution equations by means of the ADM-formalism (Arnowitt et al. (1962)) and introduced a non-uniform radial grid to ensure stable numerical evolution when using realistic equations of state.
The different oscillation patterns of a neutron star are characterized by their restoring force, e.g., p(pressure)-modes, g(gravity)-modes, i(Coriolis)-modes, s(shear)-modes or w(wave)-modes. The f-mode is the fundamental mode of the p-mode sequence and it is the oscillation mode most likely to be excited in violent processes such as neutron star formation by supernova core collapse (Torres-Forné et al. (2018; 2019)), the pre-merger interaction of neutron stars (see, e.g., Lai (1994); Kokkotas and Schafer (1995); Fuller and Lai (2011); Gold et al. (2012); Chirenti et al. (2017); Chaurasia et al. (2018); Suvorov and Kokkotas (2020); Kuan et al. (2021a,b)), the early-post-merger oscillations of the final object (Shibata and Taniguchi (2006); Bauswein and Janka (2012); Hotokezaka et al. (2013); Bernuzzi et al. (2014); Bauswein et al. (2014); Lehner et al. (2016); Bauswein et al. (2017); Rezzolla and Takami (2016); Takami et al. (2015); De Pietri et al. (2018); Breschi et al. (2019)). In the case that the merging neutron stars are of relatively small mass and the post-merger object is a fast spinning neutron star, the unstable f-mode oscillations can lead to its spin-down (Doneva et al. (2015)). The f-mode is associated with major density variations and thus can potentially be an emitter of copious amounts of gravitational radiation. The emission of gravitational waves is the primary reason for the mode’s rapid damping, at least for newly born neutron stars.
The efforts to associate the patterns of oscillations with the bulk parameters of the stars, e.g., their mass, radius or equation of state (henceforth EoS) was initiated in the mid-1990s and continued for almost two decades, advancing the field of gravitational wave asteroseismology (Andersson and Kokkotas (1996, 1998); Kokkotas et al. (2001); Benhar et al. (2004); Tsui and Leung (2005); Lau et al. (2010)). To date, very robust empirical relations have been derived for non-rotating neutron stars, connecting observables such as frequency, damping time, or moment of inertia I to the bulk properties; for example, relations of the form
In the era of gravitational wave astronomy, the various oscillation patterns (traced already in non-linear numerical simulations, e.g., Stergioulas et al. (2004); Kastaun (2006); Baiotti et al. (2008, 2009); Zink et al. (2010); Stergioulas et al. (2011); Bernuzzi et al. (2015)), if observed, can provide a wealth of information about the emitting sources, and their effects can leave their imprints both in the gravitational but also in the electromagnetic spectrum. Moreover, recent studies relate the f-mode frequencies to the Love numbers (Wen et al. (2019); Pratten et al. (2020); Andersson and Pnigouras (2021); Manoharan et al. (2021)) and even to the postmerger short gamma-ray bursts (Chirenti et al. (2019)).
Nearly all the previously mentioned studies were concerned with non-rotating stars. In nature, neutron stars will always rotate and their rotation rate may reach extreme values. The inclusion of rotation proves difficult since the extreme rotation rates that neutron stars may (and do) reach do not allow to neglect the star’s oblateness which removes the spherical symmetry from the system; this in turn makes the mathematical formulation much more involved. As a first approximation, rotation was treated perturbatively, too, which allowed considering even rotating stars as spherically symmetric (Hartle (1967)). In this so-called slow-rotation approximation, the perturbation equations gain considerably in complexity and have been written down in a gauge introduced by Regge and Wheeler (1957) first by Kojima (1992). Even though the problem remains one-dimensional, its solution is not straightforward as (among other technicalities) the outgoing-wave boundary condition at infinity is elusive. Notwithstanding, Andersson (1998) successfully applied this formalism and discovered that r-modes are prone to the so-called CFS-instability, named after their discoverers Chandrasekhar (1970) and Friedman and Schutz (1975, 1978), at any rotation rate. There has been continuing effort using the slow-rotation approximation concerning rotational modes (Lockitch et al. (2000; 2003)), and also employing different gauges (Ruoff et al. (2002); Vavoulidis et al. (2008); Passamonti et al. (2008)), but with the pressing need for frequencies of rapidly rotating neutron stars, the interest slowly faded.
Even though the slow-rotation approximation has proven fruitful in the understanding of neutron star physics, it is no longer applicable when considering neutron stars at arbitrary rotation rates, which is essential for nascent neutron stars or post-merger configurations in the immediate aftermath of a binary merger. Without the spherical symmetry of the problem, one has to account for at least two spatial dimensions which complicates the equations further and amplifies the computational expense; furthermore, it remains elusive how to formulate the outgoing-wave boundary condition at infinity for the spacetime perturbations in two spatial dimensions which essentially removes the possibility to formulate a corresponding eigenvalue problem. This issue can be circumvented by adopting the Cowling approximation (Cowling (1941)), in which the spacetime is considered static, also leading to a considerable simplification of the perturbation equations. Ignoring the impact of a dynamic spacetime (which is most severe for the quadrupolar f-mode), Yoshida and Eriguchi (1997, 1999) computed quadrupolar f-mode frequencies of rapidly rotating neutron stars and studied the associated CFS-instability in the late 1990s. Boutloukos and Nollert (2007) revived this approach and investigated the general properties of the spectrum of neutron stars regarding the acoustic and Coriolis-driven modes. As a further step toward a more general relativistic treatment, Yoshida (2012) revisited the problem in the conformal flatness approximation. However, with the mathematical difficulties of extending the eigenvalue formulation to include a dynamic spacetime, the focus shifted to study the oscillation spectra by evolving the perturbation equations in time.
Despite their complexity, the perturbation equations for rapidly rotating relativistic stars have been written down by Priou (1992), even though they were not approached numerically at that time. Gaertig and Kokkotas (2008, 2009) worked in the Cowling approximation and successfully extracted f- and g-mode frequencies of arbitrarily uniformly and fast rotating neutron stars and even of differentially rotating ones Krüger et al. (2010) by adding artificial viscosity (also known as Kreiss and Oliger (1973) dissipation) to their evolution equations in order to stabilise their time evolutions.
During the first decade of the new millennium, substantial advances were made in the time evolution of the unperturbed, non-linear Einstein equations, mostly driven by the aim to simulate compact binary mergers but also applicable to isolated neutron stars. These systems have hardly any symmetries that can be exploited to reduce the complexity of the problem, requiring to carry out the time evolutions on a three-dimensional grid. The upside of which is that essentially no constraints have to be placed on the rotational profile when simulating the dynamics of a neutron star. Such codes have been seen as a promising new approach to the calculation of mode frequencies of rapidly (and differentially) rotating neutron stars and already at the beginning of the decade, the frequencies of axisymmetric modes in the Cowling approximation (Font et al. (2000; 2001)) and those of (quasi-)radial modes in full general relativity (Font et al. (2002)) had been reported. The non-linear codes kept evolving and were used to generate mode frequencies of f-modes in the conformal flatness approximation (Dimmelmeier et al. (2006)) or those of inertial modes in the Cowling approximation (Kastaun (2008)). Baiotti et al. (2009). Not much later, the frequencies of non-axisymmetric modes in full general relativity of non-rotating polytropic neutron stars (Baiotti et al. (2009)) and soon those of rapidly rotating polytropic neutron stars (Zink et al. (2010)) were obtained from fully non-linear simulations. Even though successful, this approach to computing the frequencies of non-axisymmetric modes, however, has not been followed closely, which is also due to the computational expense associated with such numerical simulations and the accompanying limited accuracy.
In fact, from the point of view of gravitational wave detectability of oscillation modes, the most relevant scenarios are likely to involve rapidly rotating stars. Unfortunately, the aforementioned empirical relations cannot be trivially extended to rotating stars. Rotation splits the oscillation spectra in a similar fashion as the Zeeman splitting of the spectral lines due to the presence of magnetic fields. In rotating stars, the splitting leads to perturbations propagating in the direction of rotation (so-called co-rotating modes) and perturbations traveling in the opposite direction (counter-rotating modes). The oscillation frequency as observed by an observer at infinity will either increase or decrease depending on the propagation direction of the waves; for slow rotation there will be a shift of the form
Throughout this article, we employ units in which c = G = M⊚ = 1.
2 Perturbation Equations
2.1 Background Configuration
We are going to work with the Einstein equations along with the law for the conservation of energy-momentum,
where Gμν is the Einstein tensor and Tμν is the energy-momentum tensor.
We restrict ourselves to the study of the dynamics of small perturbations around an equilibrium configuration which allows us to linearise Eq. 1. We assume an axisymmetric, stationary background configuration for which the metric written in quasi-isotropic coordinates takes the form
Here, ν, ψ, μ, and ω are the four unknown metric potentials, depending only on r and θ.
We model the neutron star to be a perfect fluid without viscosity for which the corresponding energy-momentum tensor takes the form
where ϵ is the energy density, p is the pressure, and uμ the 4-velocity of the fluid. The only two non-vanishing components of the 4-velocity are linked via the star’s angular rotation rate, uφ = Ωut, and by means of the normalisation of the 4-velocity they are given by
After specifying an EoS, which may be a polytropic or a tabulated one, linking energy density and pressure to each other, we generate uniformly rotating equilibrium configurations using the
rns-code (Stergioulas and Friedman (1995); Nozawa et al. (1998); Stergioulas (1995)).
We considered sequences of neutron stars along which we keep either the central energy density constant or the baryon mass fixed (the latter are also known as evolutionary sequences) with rotation rates up to their respective mass-shedding limit. Our neutron star models were based on polytropic and realistic EoSs. We considered polytropic models with three different polytropic indices N = 0.6849, 0.7463, and 1.0. Furthermore, we employed piecewise-polytropic approximations, introduced by Read et al. (2009), for the four tabulated EoSs (APR4, H4, SLy, and WFF1) that we used. Our nonrotating configurations have gravitational masses M ∈ [1.17, 2.19]M⊚. Even though current astrophysical constraints play a role in our particular choice of EoSs, it is largely motivated by our desire to provide robust universal relations by covering a wide part of the parameter space.
2.2 Perturbation Equations
As usual in perturbative studies, we decompose the metric as
The metric perturbations are not unique but possess gauge freedom which can be utilised in different ways. Often, the gauge freedom is used to eliminate some of the spacetime perturbations, e.g., by using the well-known gauge by Regge and Wheeler (1957), and hence to reduce the number of perturbation equations. In our studies, however, we followed a different approach (which we will reason below) and opt for the Hilbert gauge, which is the gravitational equivalent to the well-known Lorenz gauge in electromagnetism, specified by
In the Hilbert gauge, the perturbed Einstein tensor takes the form.
Our choice of gauge, namely the Hilbert gauge, does not eliminate any of the metric perturbations; hence, we need 10 variables to describe the spacetime perturbations. As in previous studies of rotating neutron stars (Vavoulidis (2007); Krüger et al. (2010)), we use four variables for the fluid. In total, we need to evolve 14 variables in time. The evolution equations follow in a very straightforward manner from the perturbed Einstein equations
and the perturbed law for the conservation of energy-momentum
The evolution equations themselves are quite lengthy and not very enlightening and we refer the reader to Krüger and Kokkotas (2020b) for details in the derivation and their implementation.
3.1 Universal Relations for Single Neutron Stars
As shown in Krüger and Kokkotas (2020b), our code produces results in excellent agreement with previously published values (Zink et al. (2010); Chaurasia et al. (2018)) and our convergence tests demonstrate an accuracy of the obtained frequencies of 1–2%. In this article, we will provide some highlighted results in order to demonstrate the existence of asteroseismological relations of various types and we lay out the way that one can make use of these relations in analyzing gravitational wave signals.
More specifically, we will show different universal relations providing accurate estimates for the f-mode frequency given some bulk parameters of the star and vice versa. First, we observe a universal behavior of the f-mode frequency σi as observed in the inertial frame as a function of the star’s angular spinning frequency Ω along sequences of fixed central energy density models when we normalize both frequencies with the f-mode frequency σ0 of the corresponding non-rotating star. Figure 1 displays this behavior for more than 230 different neutron star models of each the co-rotating (i.e., stable) and counter-rotating (i.e., potentially unstable) branches of the f-mode for seven EoSs and various central energy densities (with corresponding central rest mass densities ρc ∈ [2.2, 7.3]ρ0, where ρ0 = 2.7 × 1014 g/cm3 is the nuclear saturation density); we model the universal behavior using the quadratic function
The results of a least squares fit are
FIGURE 1. Universal relations for the l =|m|=2 f-mode frequencies for sequences of constant central energy density as observed in the inertial frame. The graph shows the results from 21 such sequences (three sequences per EoS; three polytropic and four realistic EoSs). The potentially unstable f-mode branch displays a strikingly universal behavior; the largest deviations from a quadratic fit occur close to the mass-shedding limit of the sequences, mainly for the polytropic EoSs. Figure taken from Krüger and Kokkotas (2020a) with permission by APS.
We point out that our model predicts that the unstable branch of the quadrupolar f-mode becomes susceptible to the CFS instability once the angular rotation rate of the star exceeds
The stable branch of the f-mode can be fitted more accurately when switching to the comoving frame and considering sequences of constant baryon mass. The frequency σc observed in the comoving frame is related to the frequency observed in the inertial frame via σc = σi + mΩ/2π. We show our results for more than 120 different neutron star models using four realistic EoSs in Figure 2. We fit our results to the quadratic function
note that we use the Kepler velocity ΩK to normalize the star’s rotation rate in this formula. The results of a least squares fit are
FIGURE 2. Universal relations for the l =|m|=2 f-mode frequencies for sequences of constant baryon mass as observed in the comoving frame. The graph shows the results from 12 such sequences (three sequences per EoS). The stable f-mode branch displays universal behavior. Figure taken from Krüger and Kokkotas (2020a) with permission by APS.
In earlier studies for non-rotating models, fitting relations of the form
Another fitting relation which can easily be implemented in solving the inverse problem is incorporating the effective compactness
FIGURE 3. The scaled f-mode frequency of the potentially unstable branch in dependence of the effective compactness for different values of
The lines of constant
The fitting formula (12) has the advantage that it does not rely on specifically defined sequences of neutron stars, along which a particular property is held constant. For example, Eq. 11 depends on the f-mode frequency σ0 of the (in a very particular fashion) corresponding non-rotating configuration, which may not even exist in some cases (e.g., for supramassive neutron stars supported by rotation); the latter model, cf. Eq. 12, is satisfied with bulk properties of the star of which we want to know the oscillation frequency and vice versa. Another benefit of this formulation is that (as demonstrated by Doneva and Kokkotas (2015)) a similar formula can be derived for higher multipoles, i.e., l ≥ 3. Fitting formula (12) can be useful in imposing further constraints on the parameters of the postmerger objects since it combines the mass and spin of the resulting object with the f-mode frequency and, viaη, the moment of inertia I or the compactness M/R. Thus, the latter two can be further constrained by an observation of an f-mode signal, as mass and potentially spin can be extracted from the premerger and early postmerger analysis of the signal. The situation becomes more attractive if both co- and counter-rotating modes or other combination of modes are observed since only the mass of the postmerger object will be needed to constrain its parameters by using only the asteroseismological relations (see Gaertig and Kokkotas (2011); Doneva and Kokkotas (2015); Vretinaris et al. (2020)). This will be an independent yet complementary constraint in the estimation of the radius in addition to those based on the Love numbers (see Dudi et al. (2018); Abbott et al. (2018); De et al. (2018); Radice et al. (2018)).As a graphical illustration of the behavior of the f-mode frequency across the entire parameter space of stable equilibrium models, we present in Figure 4 the frequency of the counter-rotating branch as obtained from the time evolutions exemplary for EoS H4; the graph will be qualitatively similar for other EoSs. We constructed several hundred neutron star models across the entire M-Re plane along so-called evolutionary sequences; those are sequences of differently fast rotating neutron stars that share the same baryon mass. The two-dimensional plane of equilibrium models has four distinct boundaries: the static limit at the lower left along which the non-rotating models are located; second, the mass-shedding limit bounds the equatorial radius to the right; next, in the top left of the graph, the limit of stability with respect to quasi-radial perturbations connects the (non-rotating) maximum mass TOV model to (approximately) the heaviest model2; equilibrium configurations that are located to the left of this line will collapse to a black hole upon quasi-radial perturbation. Last, in line with current theory and observations (Suwa et al. (2018); Martinez et al. (2015)), we limit ourselves to neutron stars with masses M ≳ 1.17M⊚. The dots in Figure 4 depict the considered equilibrium models. The dash-dotted line slightly above M = 2.0 M⊚ is the evolutionary sequence with the baryon mass M0 = 2.3 M⊚ and it separates the supramassive neutron stars (above that line) which are merely supported by centrifugal force from those (below that line) which may spin down ending up at a stable non-rotating configuration. Other evolutionary sequences can be imagined by combining the dots with lines parallel to the dash-dotted one.
FIGURE 4. The frequency of the 2f2-mode for the entire EoS H4 is displayed color coded and some contour lines are shown. Each black dot indicates a neutron star model for which we have calculated its nonaxisymmetric mode frequencies. All neutron stars located above the (nearly horizontal) dash-dotted line are supramassive. The red (thick) contour line at 0.0 kHz separates the stable models from those that are susceptible to the CFS-instability. Figure taken from Krüger and Kokkotas (2020a) with permission by APS.
The graph shows that for this particular EoS each neutron star model may become CFS-unstable if it is sufficiently spun up. For heavier neutron stars, this happens considerably below the mass-shedding limit and sufficiently heavy supramassive models (approximately M ≳ 2.25M⊚) will inevitably be CFS-unstable (with respect to the quadrupolar f-mode); i.e. those stars will be stabilized merely by viscous mechanisms counteracting the CFS-instability.
3.2 Universal Relations Involving Long-Lived Remnants From BNS Mergers
Inspired by the works on universal relations for single neutron stars (Andersson and Kokkotas (1998); Benhar et al. (2004); Yagi and Yunes (2017)), the last 5 years have given rise to universal relations for BNS: they relate the pre-merger neutron stars to the early post-merger remnant, and have been developed using numerical relativity simulations (Bernuzzi et al. (2015); Rezzolla and Takami (2016); Kiuchi et al. (2020)).
These works have primarily focused on relating the tidal deformability of the pre-merger stars, which impact the dynamics of the pre-merger gravitational waves at leading order through the reduced tidal deformability or binary tidal deformability
Recently, Vretinaris et al. (2020) investigated empirical relations for BNS mergers based on the extensive CoRe data set of numerical relativity gravitational wave simulations (Dietrich et al. (2018)). Covering a wide range of mass ratios, they find an extensive set of universal relations involving the various peak frequencies of the post-merger gravitational wave signal, involving, e.g., the chirp mass and characteristic radius of a 1.6 M⊚ neutron star. In particular, they also find universal relations between the binary tidal deformability and the primary f-mode frequency of the post-merger signal (as in Kiuchi et al. (2020)), however, this time involving the chirp mass of the BNS. Other established universal relations connect the post-merger peak frequency to the tidal coupling constant (Bernuzzi et al. (2015); Breschi et al. (2019)).
The universal relation between the binary tidal deformability of the BNS and the stable, co-rotating f-mode frequency σs of the early, differentially rotating remnant proposed in Kiuchi et al. (2020), led us in Manoharan et al. (2021) to derive a similar relation for a potentially long-lived, uniformly rotating remnant: the relation takes the form
We also derive a relation for the potentially unstable, counter-rotating f-mode frequency of the long-lived remnant, presenting the possibility of predicting the onset of the earlier mentioned CFS-instability.
Combining these results with the universal relation, Eq. 11, for fast rotating neutron stars between the stable, co-rotating f-mode frequency and the effective compactness
that relates the pre-merger binary tidal deformability of the BNS with the effective compactness of the long-lived remnant. For rapidly rotating remnants, this relation achieves an average relative error of 2.4%.
Finally, by directly relating these quantities without going via the f-mode, we obtain a universal relation of the form
This relation achieves improved accuracy, reaching an average relative error of ∼ 1.5% for remnants with any rotation frequency. We show the quadratic fit in Figure 5 for the symmetric case of q = 1.
We finally also consider a direct relation between the binary tidal deformability and the compactness C = M/R of the long-lived remnant. Such a relation would allow the direct estimation of the remnant’s radius R using independent estimates of its gravitational mass. We propose a relation of the form
which, however, only achieves an accuracy an order of magnitude worse than for the effective compactness relation, reaching an average relative error of ∼ 8.8%. A graphic representation of this fit along with the data is shown in Figure 6.
The results presented in Manoharan et al. (2021) represent a first step towards finding universal relations between the pre-merger neutron stars and the potential long-lived remnant of a BNS merger using perturbative calculations. Our approach can be freely extended to, e.g., hot EoSs, phase transitions, as well as differential rotation for the remnant to, to cover, e.g., earlier parts of the post-merger phase.
4 Approaching the Inverse Problem
Having laid out the procedure to extract fluid oscillation frequencies from time evolution in a general relativistic framework and having developed various universal relations between bulk properties of neutron stars and their f-mode frequencies, we now turn to one of the possible applications. The field of gravitational wave asteroseismology attempts to invert the above outlined procedure and aims to constrain bulk properties of neutron stars such as mass and radius given one or more oscillations frequencies. The above derived universal relations are a valuable tool for this as they do not depend on the underlying and hitherto poorly understood nuclear equation of state. However, even without universal relations at hand, the inverse problem can be approached.
We will summarise an approach of the inverse problem using Bayesian methods as laid out by Völkel et al. (2021). Here, the fundamental idea is tested on the f-mode frequency, but it may easily be extended to other fluid modes or combined with their damping times. The use of the Bayesian framework is advantageous (compared to, e.g., an analytical inversion of the non-linear universal relation) as it directly allows to incorporate error bars of the measurement of the f-mode into the calculation. Besides the Bayesian approach and universal relations, there are in principle also semi-classical techniques, e.g., WKB theory, which can be used to address the inverse spectrum problem in simplified cases; see Völkel and Kokkotas (2019), for using axial w-modes of spherically symmetric neutron stars and Völkel and Kokkotas (2017) using a similar approach for ultra compact stars.
In a first step, we pick an EoS and a particular (rotating) neutron star model for which we then calculate two f-mode frequencies. We assume those to have a relative error of 3% (Gaussian), while our prior knowledge on M and R is uninformative. Two data points should in principle suffice to uniquely pinpoint one neutron star model (assuming a cold EoS). We show the result of the Bayesian analysis in Figure 7. The initially chosen model is depicted by the red cross at M = 1.8 M⊚ and R = 15 km. The blue shaded area in the big panel shows the correlations of mass and radius during sampling, whose posterior distributions are shown in the side panels. While the radius of the star can be reconstructed with fairly small error bars, the reconstructed mass is considerably less constrained. Nonetheless, the peak of the probability distribution nicely resembles the initially chosen model. If the f-mode frequencies were to be known more accurately, the corresponding error bars of reconstructed mass and radius would be smaller, too.
FIGURE 7. Here we compare the EoS method assuming the H4 EoS (blue) and the MPA1 EoS (orange). The diagonal panels show the sampled posterior distribution of M and R, while the off diagonal panel combines a scatter plot with logarithmic contour lines. The red cross and red line indicate the true H4 parameters that belong to the assumed 3% f-mode data. For both cases we assume that the mass is known by 10%. Figure taken from Völkel et al. (2021) with permission by APS.
This analysis required the assumption of a particular EoS and so far we used the same EoS to reconstruct the model that we also used to generate the underlying neutron star model. If we use a different EoS to perform the Bayesian analysis, we will obviously reconstruct a different model; this is shown in orange in Figure 7. The reconstructed star in this case is considerably smaller and lighter. Note that the apparent cut of the MPA1 EoS posteriors (in particular for the mass) is not an indication that this EoS has not been used for the injection, since a similar behavior can also be found for the H4 EoS, if the injected parameters are closer to the edge of the H4 EoS neutron star parameter space. One needs further information, such as the precise mass, the tidal deformability, or perhaps a damping time, in order to rule out this EoS.
Instead of making an educated guess for the “correct” EoS and reconstructing mass and radius as just described, we may also use the universal relation proposed in Eq. 12. We will not need to assume any EoS but the universal relation will provide us—given two f-mode frequencies and a good estimate of the star’s mass—with an estimate for the stars’s rotation rate Ω and its effective compactness η.
We again use the Bayesian method to invert the universal relation and we apply it to the same model as in the previous analysis; the H4 model with M = 1.8 M⊚ and R = 15 km. In Figure 8 we show the results of the analysis. On the x-axis, we show the star’s effective compactness and rotation rate, both normalised to the actual value corresponding to the equilibrium configuration. The two graphs then show the posterior distributions. We ran the same analysis twice, once assuming a 30% relative error on the prior mass (solid lines) and once assuming a 10% relative error (dashed lines). As expected, with a less informed prior on the mass, the posterior distribution for the desired quantities have a considerably larger variance.
FIGURE 8. Here we show the posterior distributions of the rotation rate Ω (top panel) and effective compactness η(bottom panel) normalized to the injected H4 values (Ω0, η0). Posteriors are obtained by using the EoS method with H4 EoS (blue), the MPA1 EoS (orange), as well as the universal relation method (green). Solid lines correspond to 30% relative error on the prior mass M and dashed lines to 10%. We indicate each mean of the posteriors as vertical lines. The f-mode relative error is assumed to be 3%. Figure taken from Völkel et al. (2021) with permission by APS.
The blue and orange curves correspond to the previously describe method where we assume a particular EoS. We have described this above and here we show the posteriors for η and 1Ω. The green curves show the posterior distribution employing the universal relation.
It is evident that the H4 EoS method (blue lines) and universal relation (UR) method (green lines) yield very similar results for the rotation rate Ω, while assuming the MPA1 EoS (orange lines) indicates a value that is larger than the correct one. Note that both observations hold independent of the specific prior knowledge of M assumed here (30 and 10%).
The situation for the effective compactness η is qualitatively different. First, the prior knowledge of M plays a big role for the UR method, but is less important for the EoS methods. For those we find that the correctly assumed H4 EoS is almost independent of uncertainties in M, while the posterior distribution obtained by the MPA1 EoS is shifted. Note that the rather different scaling behavior of the UR method is in agreement with the findings described above.
Finally, while the posteriors of Ω are very smooth, one observes small “bumps” for the H4 EoS, e.g., at η/η0 ≈ 1.02. We have verified that this likely is an artifact from the finite resolution and particular range of the used H4 f-mode data that is available to us. This directly sets the scale of how precise our currently implemented EoS data can be used to resolve the underlying parameters, which is of order percent level.
5 Summary and Outlook
We report the first extraction of frequencies of the l = |m| = 2 f-mode of general relativistic, rapidly rotating neutron stars without the commonly used slow-rotation or Cowling approximation to an extent that allows us to generalize the findings into universal relations. This concludes a long-standing open problem, building upon the effort from numerous studies throughout the past 5 decades.
For this, we have derived a set of time evolution equations governing the perturbations of the fluid of a rapidly rotating neutron star as well as its surrounding spacetime, derived in a perturbative framework in full general relativity. We have opted for the Hilbert gauge in order to arrive at a set of fully hyperbolic equations for the spacetime perturbations whose implementation does not pose major obstacles.
Convergence tests to study the accuracy of our code reveal that in a typical simulation, the obtained frequency usually deviates less than 1–2% from the limiting value or has an accuracy of about 10 Hz. As we model axisymmetric configurations only, we can reduce the problem to two spatial dimensions which drastically lowers the computational expense of our numerical time evolutions in comparison to non-linear codes that perform three-dimensional simulations. The evolution of the perturbations of a neutron star for 15 ms on a grid with 3,120 × 50 points, which facilitates a decent frequency resolution, requires only a few dozen of CPU hours; this enables us to study broad ranges of parameters and various EoSs. We expect to further reduce the computational cost even further by deriving a simplified set of perturbation equations.
As a result, we provide different universal relations for the frequencies of the l = |m| = 2 f-modes of uniformly rotating neutron stars at zero temperature which are independent of the EoS; the proposed formulae are calibrated to several hundred neutron star models that are constructed using both polytropic and realistic EoSs and are scattered across the entire parameter space of equilibrium solutions. Furthermore, it is possible to link the pre-merger binary tidal deformability to the effective compactness of the late post-merger remnant in a universal manner, too. Such universal relations will be an essential piece in the asteroseismological toolkit once the third-generation GW observatories will be able to pick up the ring-down and fluid ringing signal following the merger of a binary neutron star system; they allow to solve the inverse problem, leading to significantly tighter constraints for mass and radius of the postmerger object. For this task, it is elementary to have a smorgasbord of universal relations at hand, which allows to make a practical choice, depending on which observables are available or in which of the star’s properties one is interested. We will extend the present list of such universal relations in future articles, utilizing different combinations of bulk properties of the star; while we may obviously be (and already have been) inspired by previously published fitting formulae that were derived using different approximative frameworks, we need to be open-minded about models involving novel combinations of observables.
We also report the discovery of an accurate estimate for the onset of the CFS-instability when the f-mode frequency of the non-rotating member of the family is known and verified the corresponding critical value of T/|W|.
A natural extension of the present study will be a more comprehensive investigation of the spectrum of neutron stars (i.e. higher multipole f-modes, low p-modes and g-modes as well as w-modes) which may be excited in different astrophysical processes. Furthermore, we are going to extend our code to account for differentially rotating neutron stars and hot EoSs that are particularly relevant for nascent neutron stars or postmerger configurations in the immediate aftermath of a binary merger, both of which will have a considerable impact on the vibration frequencies or the onset of the CFS-instability (and via two further scaling parameters also on the universal relations) during a very short but dynamic interval of their lives.
Furthermore, a precise knowledge of the spectrum of compact objects is invaluable as also isolated neutron stars as well as those in inspiraling binary systems may possess high spin rates and various oscillation modes, which may be excited, e.g., via glitches or tidal coupling, may impact their electromagnetic emission or the dynamic of the whole system.
CK and KK contributed the majority of the manuscript. PM contributes the section on “Universal relations involving long-lived remnants from BNS mergers”. SV contributes the section on “Approaching the Inverse Problem”.
CK acknowledges financial support through DFG research Grant Nos. 413873357. SV acknowledges financial support provided under the European Union’s H2020 ERC Consolidator Grant “GRavity from Astrophysical to Microscopic Scales” grant agreement no. GRAMS-815673. A part of the computations were performed on Trillian, a Cray XE6m-200 supercomputer at UNH supported by the NSF MRI program under Grant Nos. PHY-1229408.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
2It is well-known that the uniformly rotating equilibrium models with the fastest rotation rate, the highest angular momentum, or the largest gravitational mass may be distinct and it depends on the EoS whether or not those are stable with respect to quasi-radial perturbations (see Cook et al. (1994)).
Abbott, B. P., Abbott, R., Abbott, T. D., Acernese, F., Ackley, K., Adams, C., et al. (2018). Erratum: Photothermal Engineering of Graphene Plasmons. Phys. Rev. Lett. 121, 161101. doi:10.1103/PhysRevLett.121.161101
Allen, G., Andersson, N., Kokkotas, K. D., and Schutz, B. F. (1998). Gravitational Waves from Pulsating Stars: Evolving the Perturbation Equations for a Relativistic star. Phys. Rev. D 58, 124012. doi:10.1103/physrevd.58.124012
Baiotti, L., Bernuzzi, S., Corvino, G., de Pietri, R., and Nagar, A. (2009). Gravitational-wave Extraction from Neutron-star Oscillations: Comparing Linear and Nonlinear Techniques. Phys. Rev. D 79, 024002. arXiv:0808.4002 [gr-qc]. doi:10.1103/physrevd.79.024002
Baiotti, L., Giacomazzo, B., and Rezzolla, L. (2008). Accurate Evolutions of Inspiralling Neutron-star Binaries: Assessment of the Truncation Error. Phys. Rev. D 78, 084033. doi:10.1103/physrevd.78.084033
Breschi, M., Bernuzzi, S., Zappa, F., Agathos, M., Perego, A., Radice, D., et al. (2019). Kilohertz Gravitational Waves from Binary Neutron star Remnants: Time-Domain Model and Constraints on Extreme Matter. Phys. Rev. D 100, 104029. arXiv:1908.11418 [gr-qc]. doi:10.1103/physrevd.100.104029
Chaurasia, S. V., Dietrich, T., Johnson-McDaniel, N. K., Ujevic, M., Tichy, W., and Brügmann, B. (2018). Gravitational Waves and Mass Ejecta from Binary Neutron star Mergers: Effect of Large Eccentricities. Phys. Rev. D 98, 104005. doi:10.1103/physrevd.98.104005
De Pietri, R., Feo, A., Font, J. A., Löffler, F., Maione, F., Pasquali, M., et al. (2018). Convective Excitation of Inertial Modes in Binary Neutron Star Mergers. Phys. Rev. Lett. 120, 221101. doi:10.1103/physrevlett.120.221101
Dietrich, T., Radice, D., Bernuzzi, S., Zappa, F., Perego, A., Brügmann, B., et al. (2018). CoRe Database of Binary Neutron star Merger Waveforms. Class. Quan. Grav. 35, 24LT01. doi:10.1088/1361-6382/aaebc0
Dimmelmeier, H., Stergioulas, N., and Font, J. A. (2006). Non-linear Axisymmetric Pulsations of Rotating Relativistic Stars in the Conformal Flatness Approximation. Monthly Notices R. Astronomical Soc. 368, 1609–1630. doi:10.1111/j.1365-2966.2006.10274.x
Doneva, D. D., Gaertig, E., Kokkotas, K. D., and Krüger, C. (2013). Gravitational Wave Asteroseismology of Fast Rotating Neutron Stars with Realistic Equations of State. Phys. Rev. D 88, 044052. doi:10.1103/physrevd.88.044052
Dudi, R., Pannarale, F., Dietrich, T., Hannam, M., Bernuzzi, S., Ohme, F., et al. (2018). Relevance of Tidal Effects and post-merger Dynamics for Binary Neutron star Parameter Estimation. Phys. Rev. D 98, 084061. doi:10.1103/physrevd.98.084061
Font, J. A., Dimmelmeier, H., Gupta, A., and Stergioulas, N. (2001). Axisymmetric Modes of Rotating Relativistic Stars in the Cowling Approximation. Monthly Notices R. Astronomical Soc. 325, 1463–1470. doi:10.1046/j.1365-8711.2001.04555.x
Font, J. A., Goodale, T., Iyer, S., Miller, M., Rezzolla, L., Seidel, E., et al. (2002). Three-dimensional General Relativistic Hydrodynamics. 2. Long Term Dynamics of Single Relativistic Stars. Phys. Rev. D 65, 084024. doi:10.1103/physrevd.65.084024
Font, J. A., Stergioulas, N., and Kokkotas, K. D. (2000). Non-linear Hydrodynamical Evolution of Rotating Relativistic Stars: Numerical Methods and Code Tests. MNRAS 313, 678–688. doi:10.1046/j.1365-8711.2000.03254.x
Glampedakis, K., and Gualtieri, L. (2018). “Gravitational Waves from Single Neutron Stars: An Advanced Detector Era Survey,” in Astrophysics and Space Science Library, Astrophysics and Space Science Library. Editors L. Rezzolla, P. Pizzochero, D. I. Jones, N. Rea, and I. Vidaña, 673–736. doi:10.1007/978-3-319-97616-7_12
Gold, R., Bernuzzi, S., Thierfelder, M., Brügmann, B., and Pretorius, F. (2012). Eccentric Binary Neutron star Mergers. Phys. Rev. D 86, 121501. arXiv:1109.5128 [gr-qc]. doi:10.1103/physrevd.86.121501
Hotokezaka, K., Kiuchi, K., Kyutoku, K., Muranushi, T., Sekiguchi, Y.-i., Shibata, M., et al. (2013). Remnant Massive Neutron Stars of Binary Neutron star Mergers: Evolution Process and Gravitational Waveform. Phys. Rev. D 88, 044026. arXiv:1307.5888 [astro-ph.HE]. doi:10.1103/physrevd.88.044026
Kastaun, W. (2006). High-resolution Shock Capturing Scheme for Ideal Hydrodynamics in General Relativity Optimized for Quasistationary Solutions. Phys. Rev. D 74, 124024. doi:10.1103/physrevd.74.124024
Kiuchi, K., Kawaguchi, K., Kyutoku, K., Sekiguchi, Y., and Shibata, M. (2020). Sub-radian-accuracy Gravitational Waves from Coalescing Binary Neutron Stars in Numerical Relativity. II. Systematic Study on the Equation of State, Binary Mass, and Mass Ratio. Phys. Rev. D 101, 084006. doi:10.1103/physrevd.101.084006
Kokkotas, K. D., Apostolatos, T. A., and Andersson, N. (2001). The Inverse Problem for Pulsating Neutron Stars: a 'fingerprint Analysis' for the Supranuclear Equation of State. Monthly Notices R. Astronomical Soc. 320, 307–315. doi:10.1046/j.1365-8711.2001.03945.x
Kuan, H.-J., Suvorov, A. G., and Kokkotas, K. D. (2021b). General-relativistic Treatment of Tidal G-Mode Resonances in Coalescing Binaries of Neutron Stars. II. As Triggers for Precursor Flares of Short Gamma-ray Bursts. arXiv e-prints, arXiv:2107.00533. arXiv:2107.00533 [astro-ph.HE].
Kuan, H.-J., Suvorov, A. G., and Kokkotas, K. D. (2021a). General-relativistic Treatment of Tidal G-Mode Resonances in Coalescing Binaries of Neutron Stars. I. Theoretical Framework and Crust Breaking. arXiv e-prints, arXiv:2106.16123. arXiv:2106.16123 [gr-qc].
Lehner, L., Liebling, S. L., Palenzuela, C., Caballero, O. L., O’Connor, E., Anderson, M., et al. (2016). Unequal Mass Binary Neutron star Mergers and Multimessenger Signals. Class. Quan. Grav. 33, 184002. arXiv:1603.00501 [gr-qc]. doi:10.1088/0264-9381/33/18/184002
Martinez, J. G., Stovall, K., Freire, P. C. C., Deneva, J. S., Jenet, F. A., McLaughlin, M. A., et al. (2015). Pulsar J0453+1559: a Double Neutron Star System with a Large Mass Asymmetry. ApJ 812, 143. doi:10.1088/0004-637x/812/2/143
Nozawa, T., Stergioulas, N., Gourgoulhon, E., and Eriguchi, Y. (1998). Construction of Highly Accurate Models of Rotating Neutron Stars - Comparison of Three Different Numerical Schemes. Astron. Astrophys. Suppl. Ser. 132, 431–454. doi:10.1051/aas:1998304
Priou, D. (1992). The Perturbations of a Fully General Relativistic and Rapidly Rotating Neutron star - I. Equations of Motion for the Solid Crust. Monthly Notices R. Astronomical Soc. 254, 435–452. doi:10.1093/mnras/254.3.435
Radice, D., Perego, A., Zappa, F., and Bernuzzi, S. (2018). GW170817: Joint Constraint on the Neutron Star Equation of State from Multimessenger Observations. ApJ 852, L29. doi:10.3847/2041-8213/aaa402
Read, J. S., Lackey, B. D., Owen, B. J., and Friedman, J. L. (2009). Constraints on a Phenomenologically Parametrized Neutron-star Equation of State. Phys. Rev. D 79, 124032. doi:10.1103/physrevd.79.124032
Ruoff, J., Stavridis, A., and Kokkotas, K. D. (2002). Evolution Equations for the Perturbations of Slowly Rotating Relativistic Stars. Monthly Notices R. Astronomical Soc. 332, 676–688. doi:10.1046/j.1365-8711.2002.05329.x
Shibata, M., and Taniguchi, K. (2006). Merger of Binary Neutron Stars to a Black Hole: Disk Mass, Short Gamma-ray Bursts, and Quasinormal Mode Ringing. Phys. Rev. D 73, 064027. doi:10.1103/physrevd.73.064027
Stergioulas, N., Apostolatos, T. A., and Font, J. A. (2004). Non-linear Pulsations in Differentially Rotating Neutron Stars: Mass-Shedding-Induced Damping and Splitting of the Fundamental Mode. MNRAS 352, 1089–1101. doi:10.1111/j.1365-2966.2004.07973.x
Stergioulas, N., Bauswein, A., Zagkouris, K., and Janka, H.-T. (2011). Gravitational Waves and Non-axisymmetric Oscillation Modes in Mergers of Compact Object Binaries. MNRAS 418, 427–436. doi:10.1111/j.1365-2966.2011.19493.x
Stergioulas, N. (1995). Center for Gravitation, Cosmology & Astrophysics. Available at: http://www.gravity.phys.uwm.edu/rns/.
Suvorov, A. G., and Kokkotas, K. D. (2020). Precursor Flares of Short Gamma-ray Bursts from Crust Yielding Due to Tidal Resonances in Coalescing Binaries of Rotating, Magnetized Neutron Stars. Phys. Rev. D 101, 083002. doi:10.1103/physrevd.101.083002
Torres-Forné, A., Cerdá-Durán, P., Passamonti, A., and Font, J. A. (2018). Towards Asteroseismology of Core-Collapse Supernovae with Gravitational-Wave Observations - I. Cowling Approximation. MNRAS 474, 5272–5286. doi:10.1093/mnras/stx3067
Torres-Forné, A., Cerdá-Durán, P., Passamonti, A., Obergaulinger, M., and Font, J. A. (2019). Towards Asteroseismology of Core-Collapse Supernovae with Gravitational Wave Observations - II. Inclusion of Space-Time Perturbations. MNRAS 482, 3967–3988. doi:10.1093/mnras/sty2854
Vretinaris, S., Stergioulas, N., and Bauswein, A. (2020). Empirical Relations for Gravitational-Wave Asteroseismology of Binary Neutron star Mergers. Phys. Rev. D 101, 084039. doi:10.1103/physrevd.101.084039
Wen, D.-H., Li, B.-A., Chen, H.-Y., and Zhang, N.-B. (2019). GW170817 Implications on the Frequency and Damping Time of f-mode Oscillations of Neutron Stars. prc 99, 045806. doi:10.1103/PhysRevC.99.045806
Yoshida, S. i., and Eriguchi, Y. (1997). Neutral Points of Oscillation Modes along Equilibrium Sequences of Rapidly Rotating Polytropes in General Relativity: Application of the Cowling Approximation. ApJ 490, 779–784. doi:10.1086/304918
Keywords: neutron star, equation of state, stars, rotation, stars, Bayesiah analysis, binary merger
Citation: Krüger CJ, Kokkotas KD, Manoharan P and Völkel SH (2021) Fast Rotating Neutron Stars: Oscillations and Instabilities. Front. Astron. Space Sci. 8:736918. doi: 10.3389/fspas.2021.736918
Received: 06 July 2021; Accepted: 15 September 2021;
Published: 28 September 2021.
Edited by:Brynmor Haskell, Nicolaus Copernicus Astronomical Center (PAS), Poland
Reviewed by:Ankan Sur, Nicolaus Copernicus Astronomical Center (PAS), Poland
Sebastiano Bernuzzi, Friedrich Schiller University Jena, Germany
Copyright © 2021 Krüger, Kokkotas, Manoharan and Völkel. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Christian J. Krüger, email@example.com,