Quasi-periodic oscillations inferred during rare magnetar giant flare tails were initially interpreted as torsional oscillations of the neutron star (NS) crust, and have been more recently described as global core+crust perturbations. Similar frequencies are also present in high signal-to-noise magnetar short bursts. In magnetars, disturbances of the field are strongly coupled to the NS crust regardless of the triggering mechanism of short bursts. For low-altitude magnetospheric magnetar models of fast radio bursts (FRBs) associated with magnetar short bursts, such as the low-twist model, crustal oscillations may be associated with additional radio bursts in the encompassing short burst event (as recently suggested for SGR 1935+2154). Given the large extragalactic volume probed by wide-field radio transient facilities, this offers the prospect of studying NS crusts leveraging samples far more numerous than galactic high-energy magnetar bursts by studying statistics of sub-burst structure or clustered trains of FRBs. We sketch what may be learned and issues that may arise in such a program. Lower $l$-number eigenmodes (corresponding to FRB time intervals of $\sim5-50$ ms) are likely less susceptible than high-$l$ modes to confusion by systematic effects associated with the NS crust physics, magnetic field, and damping. They may be more promising in their utility, and also may corroborate models where FRBs arise from mature magnetars. Future observational characterization of such signals can also determine whether they can be employed as cosmological "standard oscillators" to constrain redshift, or can be used to constrain the mass of FRB-producing magnetars when reliable redshifts are available.

The increasing richness of data related to cold dense matter, from laboratory experiments to neutron-star observations, requires a framework for constraining the properties of such matter that makes use of all relevant information. Here, we present a rigorous but practical Bayesian approach that can include diverse evidence, such as nuclear data and the inferred masses, radii, tidal deformabilities, moments of inertia, and gravitational binding energies of neutron stars. We emphasize that the full posterior probability distributions of measurements should be used rather than, as is common, imposing a cut on the maximum mass or other quantities. Our method can be used with any parameterization of the equation of state (EOS). We use both a spectral parameterization and a piecewise polytropic parameterization with variable transition densities to illustrate the implications of current measurements and show how future measurements in many domains could improve our understanding of cold catalyzed matter. We find that different types of measurements will play distinct roles in constraining the EOS in different density ranges. For example, better symmetry energy measurements will have a major influence on our understanding of matter somewhat below nuclear saturation density but little influence above that density. In contrast, precise radius measurements or multiple tidal deformability measurements of the quality of those from GW170817 or better will improve our knowledge of the EOS over a broader density range.

In this work we show the results for numerical solutions of the relativistic Grad-Shafranov equation for a typical neutron star with 1.4 solar masses. We have studied the internal magnetic field considering both the poloidal and toroidal components, as well as the behavior of the field lines parametrized by the ratio between these components of the field.

Neutron star mergers can form a hypermassive neutron star (HMNS) remnant, which may be the engine of a short gamma ray burst (SGRB) before it collapses to a black hole, possibly several hundred milliseconds after the merger. During the lifetime of a HMNS, numerical relativity simulations indicate that it will undergo strong oscillations and emit GWs with frequencies of a few kilohertz, which are unfortunately too high for detection to be probable with Advanced LIGO. Here we discuss the current and future prospects for detecting these oscillations as modulation of the SGRB. The understanding of the physical mechanism responsible for the HMNS oscillations will provide information on the equation of state of the hot HMNS, and the observation of these frequencies in the SGRB data would give us insight into the emission mechanism of the SGRB.

The quasiperiodic oscillations found in the three giant flares of soft gamma-ray repeaters observed to date have been interpreted as crustal oscillations caused by a starquake following a dramatic rearrangement of the stellar magnetic field. Motivated by these observations, we study the influence of the magnetic field geometry in the frequencies of the torsional oscillations of magnetized neutron stars. We use realistic tabulated equations of state for the core and crust of the stars and model their magnetic field as a dipole plus a toroidal component, using the relativistic Grad-Shafranov equation. The frequencies of the torsional modes are obtained by the numerical solution of the eigenvalue problem posed by the linear perturbation equations in the Cowling approximation. Our results show how the asteroseismology of these stars becomes complicated by the degeneracy in the frequencies due to the large relevant parameter space. However, we are able to propose a testable scenario in which the rearrangement of the magnetic field causes an evolution in the frequencies. Finally, we show that there is a magnetic field configuration that maximizes the energy in the perturbation at linear order, which could be related to the trigger of the giant flare.

The discovery of quasi-periodic brightness oscillations (QPOs) in the X-ray emission accompanying the giant flares of the soft gamma-ray repeaters SGR 1806-20 and SGR 1900+14 has led to intense speculation about their nature and what they might reveal about the interiors of neutron stars. Here we take a fresh look at the giant flare data for SGR 1806-20, and in particular we analyze short segments of the post-peak emission using a Bayesian procedure that has not previously been applied to these data. We find at best weak evidence that any QPO persists for more than $\sim 1$ second; instead, almost all the data are consistent with a picture in which there are numerous independently-excited modes that decay within a few tenths of a second. This has interesting implications for the rapidity of decay of the QPO modes, which could occur by the previously-suggested mechanism of coupling to the MHD continuum. The strongest QPOs favor certain rotational phases, which might suggest special regions of the crust or of the magnetosphere. We also find several previously unreported QPOs in these data, which may help in tracking down their origin.

Gravitational wave and electromagnetic observations can provide new insights into the nature of matter at supra-nuclear densities inside neutron stars. Improvements in electromagnetic and gravitational wave sensing instruments continue to enhance the accuracy with which they can measure the masses, radii, and tidal deformability of neutron stars. These better measurements place tighter constraints on the equation of state of cold matter above nuclear density. In this article, we discuss a complementary approach to get insights into the structure of neutron stars by providing a model prediction for non-linear fundamental eigenmodes (f-modes) and their decay over time, which are thought to be induced by time-dependent tides in neutron star binaries. Building on pioneering studies that relate the properties of f-modes to the structure of neutron stars, we systematically study this link in the non-perturbative regime using models that utilize numerical relativity. Using a suite of fully relativistic numerical relativity simulations of oscillating TOV stars, we establish blueprints for the numerical accuracy needed to accurately compute the frequency and damping times of f-mode oscillations, which we expect to be a good guide for the requirements in the binary case. We show that the resulting f-mode frequencies match established results from linear perturbation theory, but the damping times within numerical errors depart from linear predictions. This work lays the foundation for upcoming studies aimed at a comparison of theoretical models of f-mode signatures in gravitational waves, and their uncertainties with actual gravitational wave data, searching for neutron star binaries on highly eccentric orbits, and probing neutron star structure at high densities.

After a long wait, gravitational wave astronomy has finally begun. Binary black hole mergers are being detected by LIGO and Virgo, and theorists are starting to receive a wealth of data to be analized. At this point we can at long last begin to test alternative theories of gravity and different models of compact objects. One powerful tool to do this is the perturbative analysis of background spacetimes. The objective of this brief review is to introduce the notion and analysis of black hole quasinormal modes, starting from the linear perturbation theory and including a brief discussion of numerical methods and astrophysical implications in the gravitational wave signals recently detected. With these basic ingredients, more sophisticated analyses and applications are possible.

The puzzling existence of a number of neutron stars that appear to be in the r-mode instability window supports further investigations of the r-mode frequencies, damping times and saturation amplitudes, especially now in light of the much anticipated and exciting results from LIGO and NICER. It has been suggested by Mahmoodifar and Strohmayer that coherent frequencies found in the RXTE data of the accreting millisecond pulsar XTE J1751-305 and during the 2001 superburst of 4U 1636-536 (one of the accreting sources to be observed by NICER) could, in fact, be r-modes. Based on these observations, we expand here the results of a previous work on relativistic, rotational and differential rotation corrections to the r-mode frequency in order to provide more accurate estimates of the $M-R$ relation for these neutron stars. Finally, we compare our results with predictions from a few realistic tabulated equations of state, providing further constraints that favor more compact models. We find that, if the observed frequencies indeed correspond to r-modes, then the masses of XTE J1751-305 and 4U 1636-536 should lie in the approximate ranges $1.48 - 1.56 M_{\odot}$ and $1.60 - 1.68 M_{\odot}$, respectively.

As gravitational wave instrumentation becomes more sensitive, it is interesting to speculate about subtle effects that could be analyzed using upcoming generations of detectors. One such effect that has great potential for revealing the properties of very dense matter is fluid oscillations of neutron stars. These have been found in numerical simulations of the hypermassive remnants of double neutron star mergers and of highly eccentric neutron star orbits. Here we focus on the latter and sketch out some ideas for the production, gravitational-wave detection, and analysis of neutron star oscillations. These events will be rare (perhaps up to several tens per year could be detected using third-generation detectors such as the Einstein Telescope or the Cosmic Explorer), but they would have unique diagnostic power for the analysis of cold, catalyzed, dense matter. Furthermore, these systems are unusual in that analysis of the tidally excited f-modes of the stars could yield simultaneous measurements of their masses, moments of inertia, and tidal Love numbers, using the frequency, damping time, and amplitude of the modes. They would thus present a nearly unique opportunity to test observationally the I-Love-Q relation. The analysis of such events will require significant further work in nuclear physics and general relativistic nonlinear mode coupling, and thus we discuss further directions that will need to be pursued. For example, we note that for nearly-grazing encounters, numerical simulations show that the energy delivered to the f-modes may be up to two orders of magnitude greater than predicted in the linear theory.

The interferometric LIGO detectors have recently measured the first direct gravitational-wave signal from what has been interpreted as the inspiral, merger and ringdown of a binary system of black holes. The signal-to-noise ratio of the measured signal is large enough to leave little doubt that it does refer to the inspiral of two massive and ultracompact objects, whose merger yields a rotating black hole. Yet, the quality of the data is such that some room is left for alternative interpretations that do not involve black holes, but other objects that, within classical general relativity, can be equally massive and compact, namely, gravastars. We here consider the hypothesis that the merging objects were indeed gravastars and explore whether the merged object could therefore be not a black hole but a rotating gravastar. After comparing the real and imaginary parts of the ringdown signal of GW150914 with the corresponding quantities for a variety of gravastars, and notwithstanding the very limited knowledge of the perturbative response of rotating gravastars, we conclude it is not possible to model the measured ringdown of GW150914 as due to a rotating gravastar.

We study the $f$-mode frequencies and damping times of nonrotating neutron stars (NS) in general relativity (GR) by solving the linearized perturbation equations, with the aim to establish "universal" relations that depend only weakly on the equations of state (EOS). Using a more comprehensive set of EOSs, we re-examine some proposed linearizations that describe the $f$-mode parameters in terms of mass and radius of the neutron star (NS), and we test a more recent proposal for expressing the $f$-mode parameters as quadratic functions of the effective compactness. Our extensive results for each equation of state considered allow us to study the accuracy of each proposal. In particular, we find that the damping time deviates quite considerably from the proposed linearization. We introduce a new universal relation for the product of the $f$-mode frequency and damping time as a function of the (ordinary) compactness, which proved to be more accurate. The relations using the effective compactness on the other hand also fit our data accurately. Our results show that the maximum oscillation frequency depends strongly on the EOS, such that the measurement of a high oscillation frequency would rule out several EOSs. Lastly, we compare the exact mode frequencies to those obtained in the Cowling approximation, and also to results obtained with a nonlinear evolution code, validating the implementations of the different approaches.

Gravastars have been recently proposed as potential alternatives to explain the astrophysical phenomenology traditionally associated to black holes, raising the question of whether the two objects can be distinguished at all. Leaving aside the debate about the processes that would lead to the formation of a gravastar and the astronomical evidence in their support, we here address two basic questions: Is a gravastar stable against generic perturbations? If stable, can an observer distinguish it from a black hole of the same mass? To answer these questions we construct a general class of gravastars and determine the conditions they must satisfy in order to exist as equilibrium solutions of the Einstein equations. For such models we perform a systematic stability analysis against axial-perturbations, computing the real and imaginary parts of the eigenfrequencies. Overall, we find that gravastars are stable to axial perturbations, but also that their quasi-normal modes differ from those of a black hole of the same mass and thus can be used to discern, beyond dispute, a gravastar from a black hole.