Search

Precise measurements of neutron star observables (such as mass and radius) allow one to constrain the equations of state for supranuclear matter and develop a stronger understanding of nuclear physics. The Neutron star Interior Composition ExploreR (NICER) tracks X-ray hotspots on rotating NSs and is able to infer precise information about the compactness of the star. Gravitational waves carry information about the tidal deformability (related to the tidal Love number) of neutron stars, which has been measured by the LIGO/Virgo/KAGRA collaboration. These two observables enjoy an approximately universal property between each other that is insensitive to the equations of state (the "universal Love-C relation"). In this paper, we focus on deriving two analytic expressions for the Love-C relations that are ready-to-use and improve upon previous analytic expressions. The first model is inspired by a Newtonian polytrope, whose perturbation to the gravitational potential can be found analytically. We extend this Newtonian model to the relativistic regime by providing a quadratic fit to the gravitational potential perturbation against stellar compactness. The second model makes use of the Tolman VII model and adopts a spectral expansion with Chebyshev polynomials, which converges faster than the Taylor expansions used in previous work. We find that the first model provides a more accurate description of the Love-C relation for realistic neutron stars than the second model, and it provides the best expression among all other analytic relations studied here in terms of describing the averaged numerical Love-C relation. These new models are not only useful in practice, but they also show the power and importance of analytic modeling of neutron stars., Comment: 11 pages, 9 figures

Tidal effects on compact objects provide profound theoretical insights into the structure of the field equations, and are wonderful probes of the equation of state of matter, the nature of black holes and of the underlying theory of gravity. The natural framework for understanding tides is a perturbative scheme. Here, we point out ambiguities in determining tidal response functions within such a framework, which may lead to bias in constraining physical parameters with gravitational-wave observations if the computed quantities are not properly linked to observables. We propose a Canonical Tidal Response Function (CTRF) definition to compare values of tides in theories beyond vacuum General Relativity in a unified manner. As an example, we provide black hole tidal response functions, including both conservative and dissipative pieces, in various theories of gravity. Tidal dissipation Love numbers for black holes are derived here for the first time in most of the non-Einsteinian theories considered in this paper., Comment: 24 pages

The response of astrophysical compact objects to external tidal fields carries valuable information on the nature of these objects, on the equation of state of matter, and on the underlying gravitational theory. In this work, we highlight subtleties in describing relativistic dynamical tidal responses that arise from ambiguities in the decomposition of a perturbed metric into external tidal and induced response pieces. Observables are unambiguous. However in practice, differences arising from implicit assumptions in the definition of tidal deformabilities may lead to a bias in constraining nuclear physics or gravitational theories, if not properly tied to observational data. We propose calibration of a tidal response function for any compact objects in vacuum General Relativity. Within this framework, the dynamical tidal Love numbers of a Schwarzschild black hole in both even and odd sectors vanish at any multipole order. The calibration allows one to define dynamical tidal deformabilities of relativistic stars, such as neutron stars, as the difference from the BH values (zero) under the unified definition in a simple manner. As a straightforward extension of our framework, we compute the next-to-leading dissipative tidal response of Schwarzschild black holes for the first time., Comment: 22 pages

Gravitational waves have been shown to provide new constraints on gravitational theories beyond general relativity (GR), especially in the strong field regime. Gravitational wave signals from galactic double white dwarfs, expected to be detected by the Laser Interferometer Space Antenna (LISA), also have the potential to place stringent bounds on certain theories that give rise to relatively large deviations from GR in less compact binaries, such as through scalar radiation. Nevertheless, the orbital evolution of close double white dwarf systems is also affected by various astrophysical effects, such as stellar rotation, tidal interactions, and magnetic interactions, which add complexity to the gravity tests. In this work, we employ the parametrized post-Einsteinian model to capture the leading beyond-GR effect on the signal and estimate the measurement uncertainties using the Fisher information matrix. We then study the systematic error caused by ignoring each astrophysical effect mentioned above on the parameter estimation. Our numerical results show that, to place bounds on the non-GR effects comparable to existing bounds from pulsar observations, tight priors on the mass of the binary and long observation time are required. At this level of sensitivity, we found that systematic errors from the astrophysical effects dominate statistical errors. The most significant effects investigated here are torques from tidal synchronization and magnetic unipolar induction for sufficiently large magnetic fields ($>10^7$G). Meanwhile, even the weaker astrophysical effects from quadrupolar deformations are of a similar order of magnitude as the statistical uncertainty, and hence cannot be ignored in the waveform model. We conclude that the astrophysical effects must be carefully accounted for in the parameter estimation to test gravity with galactic double white dwarfs detected by LISA., Comment: 14 pages, 8 figures, submitted to PRD

A massive black hole can develop a dark-matter overdensity, and the dark matter changes the evolution of a stellar-mass compact object inspiraling around the massive black hole through the dense dark-matter environment. Specifically, dynamical friction speeds up the inspiral of the compact object and causes feedback on the dark-matter distribution. These intermediate mass-ratio inspirals with dark matter are a source of gravitational waves (GWs), and the waves can dephase significantly from an equivalent system in vacuum. Prior work has shown that this dephasing needs to be modeled to detect the GWs from these systems with LISA (the Laser Interferometer Space Antenna); it also showed that the density and distribution of dark matter can be inferred from a GW measurement. In this paper, we study whether the parametrized post-Einsteinian (ppE) framework can be used to infer the presence of dark matter in these systems. We confirm that if vacuum waveform templates are used to model the GWs from an inspiral in a dark-matter halo, then the resulting parameter estimation is biased. We then apply the ppE framework to determine whether it can reduce the parameter-estimation biases, and we find that adding one ppE phase term to a waveform template eliminates the parameter-estimation biases (statistical errors become larger than the systematic ones), but the effective post-Newtonian order in the ppE framework must be specified without uncertainties. When the post-Newtonian order has uncertainty, we find that the systematic errors on the ppE and the binary's parameters exceed the statistical errors. Thus, the simplest ppE framework would not give unbiased results for these systems, and a further extension of it, or dedicated parameter estimation with gravitational waveforms that include dark-matter effects would be needed., Comment: 16 pages, 8 figures

Quasiperiodic oscillations (QPOs) have been recently discovered in the short gamma-ray bursts (GRBs) 910711 and 931101B. Their frequencies are consistent with those of the radial and quadrupolar oscillations of binary neutron star merger remnants, as obtained in numerical relativity simulations. These simulations reveal quasiuniversal relations between the remnant oscillation frequencies and the tidal coupling constant of the binaries. Under the assumption that the observed QPOs are due to these postmerger oscillations, we use the frequency-tide relations in a Bayesian framework to infer the source redshift, as well as the chirp mass and the binary tidal deformability of the binary neutron star progenitors for GRBs 910711 and 931101B. We further use this inference to estimate bounds on the mass-radius relation for neutron stars. By combining the estimates from the two GRBs, we find a 68\% credible range $R_{1.4}=12.48^{+0.41}_{-0.41}$~km for the radius of a neutron star with mass $M=1.4$~M$_\odot$, which is one of the tightest bounds to date.

This review is focused on tests of Einstein's theory of general relativity with gravitational waves that are detectable by ground-based interferometers and pulsar-timing experiments. Einstein's theory has been greatly constrained in the quasi-linear, quasi-stationary regime, where gravity is weak and velocities are small. Gravitational waves are allowing us to probe a complimentary, yet previously unexplored regime: the non-linear and dynamical \emph{extreme gravity regime}. Such a regime is, for example, applicable to compact binaries coalescing, where characteristic velocities can reach fifty percent the speed of light and gravitational fields are large and dynamical. This review begins with the theoretical basis and the predicted gravitational-wave observables of modified gravity theories. The review continues with a brief description of the detectors, including both gravitational-wave interferometers and pulsar-timing arrays, leading to a discussion of the data analysis formalism that is applicable for such tests. The review then discusses gravitational-wave tests using compact binary systems, and ends with a description of the first gravitational wave observations by advanced LIGO, the stochastic gravitational wave background observations by pulsar timing arrays, and the tests that can be performed with them.., Comment: Update to Living Rev.Rel. 16 (2013) 9, now with twice the content and three times the number of references! arXiv admin note: substantial text overlap with arXiv:1304.3473

Neutron stars (NSs) provide a unique laboratory to study matter under extreme densities. Recent observations from gravitational and electromagnetic waves have enabled constraints on NS properties, such as tidal deformability (related to the tidal Love number) and stellar compactness. Although each of these two NS observables depends strongly on the stellar internal structure, the relation between them (called the Love-C relation) is known to be equation-of-state insensitive. In this study, we investigate the effects of a possible crystalline phase in the core of hybrid stars (HSs) on the mass-radius and Love-C relations, where HSs are a subclass of NS models with a quark matter core and a nuclear matter envelope with a sharp phase transition in between. We find that both the maximum mass and the corresponding radius increase as one increases the stiffness of the quark matter core controlled by the speed of sound, while the density discontinuity at the nuclear-quark matter transition effectively softens the equations of state. Deviations of the Love-C relation for elastic HSs from that of fluid NSs become more pronounced with a larger shear modulus, lower transition pressure, and larger density gap and can be as large as 60%. These findings suggest a potential method for testing the existence of distinct phases within HSs, though deviations are not large enough to be detected with current measurements of the tidal deformability and compactness., Comment: 29 pages, 12 figures, Accepted for publication in General Relativity and Gravitation

Non-radial oscillation modes of a neutron star possess valuable information about its internal structure and nuclear physics. Starting from the quadrupolar order, such modes under general relativity are known as quasi-normal modes since they dissipate energy through gravitational radiation and their frequencies are complex. The stability of these modes is governed by the sign of the imaginary part of the frequency, which determines whether the mode would decay or grow over time. In this Letter, we develop a fully consistent framework in general relativity to study quasi-normal modes of neutron stars with anisotropic pressure, whose motivation includes strong internal magnetic fields and non-vanishing shear or viscosity. We employ parametrized models for the anisotropy and solve the perturbed Einstein field equations numerically. We find that, unlike the case for isotropic neutron stars, the imaginary parts of some of the pressure ($p$-)modes flip signs as the degree of anisotropy deviates from zero, depicting a transition from stable modes to unstable modes. This finding indicates that some anisotropic neutron star models are unstable, potentially restricting the form of sustained anisotropy., Comment: 8+6 pages, 2+1 figures, re-submitted to PRD

General relativity (GR) has proven to be a highly successful theory of gravity since its inception. The theory has thrivingly passed numerous experimental tests, predominantly in weak gravity, low relative speeds, and linear regimes, but also in the strong-field and very low-speed regimes with binary pulsars. Observable gravitational waves (GWs) originate from regions of spacetime where gravity is extremely strong, making them a unique tool for testing GR, in previously inaccessible regions of large curvature, relativistic speeds, and strong gravity. Since their first detection, GWs have been extensively used to test GR, but no deviations have been found so far. Given GR's tremendous success in explaining current astronomical observations and laboratory experiments, accepting any deviation from it requires a very high level of statistical confidence and consistency of the deviation across GW sources. In this paper, we compile a comprehensive list of potential causes that can lead to a false identification of a GR violation in standard tests of GR on data from current and future ground-based GW detectors. These causes include detector noise, signal overlaps, gaps in the data, detector calibration, source model inaccuracy, missing physics in the source and in the underlying environment model, source misidentification, and mismodeling of the astrophysical population. We also provide a rough estimate of when each of these causes will become important for tests of GR for different detector sensitivities. We argue that each of these causes should be thoroughly investigated, quantified, and ruled out before claiming a GR violation in GW observations., Comment: Review article; 1 figure; 1 table; comments welcome

The response of black holes to small perturbations is known to be partially described by a superposition of quasinormal modes. Despite their importance to enable strong-field tests of gravity, little to nothing is known about what overtones and quasinormal-mode amplitudes are like for black holes in extensions to general relativity. We take a first step in this direction and study what is arguably the simplest model that allows first-principle calculations to be made: a nonrotating black hole in an effective-field-theory extension of general relativity with cubic-in-curvature terms. Using a phase-amplitude scheme that uses analytical continuation and the Pr\"ufer transformation, we compute, for the first time, the quasinormal overtone frequencies (in this theory) and quasinormal-mode excitation factors (in any theory beyond general relativity). We find that the overtone quasinormal frequencies and their excitation factors are more sensitive than the fundamental mode to the lengthscale $l$ introduced by the higher-derivative terms in the effective field theory. We argue that a description of all overtones cannot be made within the regime of validity of the effective field theory, and we conjecture that this is a general feature of any extension to general relativity that introduces a new lengthscale. We also find that a parametrization of the modifications to the general-relativistic quasinormal frequencies in terms of the ratio between $l$ and the black hole's mass is somewhat inadequate, and we propose a better alternative. As an application, we perform a preliminary study of the implications of the breakdown, in the effective field theory, of the equivalence between the quasinormal mode spectra associated to metric perturbations of polar and axial parity of the Schwarzschild black hole in general relativity. We also present a simple justification for the loss of isospectrality., Comment: 21 pages, 9 figures, 1 table. (v2) Matches version published in Phys. Rev. D

Certain relations among neutron-star observables that are insensitive to the equation of state are known to exist. Such universal relations have been shown to be valid for cold and stationary neutron stars. Here, we study these relations in more dynamic scenarios: protoneutron stars and hypermassive neutron stars. First, we study protoneutron stars. We use an effective equation of state, extracted from three-dimensional core-collapse supernova simulations, to obtain the structure of spherically symmetric protoneutron stars. We then consider nonradial oscillations to compute their $f$-mode frequency ($f$), as well as slow rotation and small tidal deformation, to compute their moment of inertia ($I$), spin-induced quadrupole moment ($Q$), and Love number. We find that well-established universal relations for cold neutron stars involving these observables ($I$-Love-$Q$ and $f$-Love relations) are approximately valid for protoneutron stars, with a deviation below $\approx$ 10$\%$ for a postbounce time above $\approx$ 0.5 s, considering eight different supernova progenitors and the SFHo equation of state. Next, we study hypermassive neutron stars. We obtain a new universal relation between the $f$-mode frequency and the compactness of cold and nonrotating neutron stars, using bulk quantities. We show that this relation has an equation-of-state-variation of $\approx$ $3\%$, considering a set of ten equations of state. Using results from binary neutron star merger simulations, we study the evolution of hypermassive neutron stars on the $f$-$C$ plane, considering two different mass ratios and the SFHo equation of state. We find that the relation between the peak frequency of the gravitational-wave signal and the compactness from these hypermassive neutron stars deviates from the universal $f$-$C$ relation by 70 $-$ 80$\%$, when the peak frequency is taken directly as a proxy for the $f$-mode.

The Laser Interferometer Space Antenna (LISA) is the first scientific endeavour to detect and study gravitational waves from space. LISA will survey the sky for Gravitational Waves in the 0.1 mHz to 1 Hz frequency band which will enable the study of a vast number of objects ranging from Galactic binaries and stellar mass black holes in the Milky Way, to distant massive black-hole mergers and the expansion of the Universe. This definition study report, or Red Book, presents a summary of the very large body of work that has been undertaken on the LISA mission over the LISA definition phase., Comment: 155 pages, with executive summary and table of contents