Accessing Universal Relations of Binary Neutron Star Waveforms in Massive Scalar-Tensor Theory

We investigate how the quasiuniversal relations connecting tidal deformability with gravitational waveform characteristics and/or properties of individual neutron stars that were proposed in the literature within general relativity would be influenced in the massive Damour-Esposito-Farese-type scalar-tensor gravity. For this purpose, we systematically perform numerical relativity simulations of ∼120 binary neutron-star mergers with varying scalar coupling constants. Although only three neutron-star equations of state are adopted, a clear breach of universality can be observed in the datasets. In addition to presenting difficulties in constructing quasiuniversal relations in alternative gravity theories, we also briefly compare the impacts of non-general-relativity physics on the waveform features and those due to the first order or cross-over quantum chromodynamical phase transition.

Finding Structure in the Speed of Sound of Supranuclear Matter from Binary Love Relations

Analyses that connect observations of neutron stars with nuclear-matter properties can rely on equation-of-state insensitive relations. We show that the slope of the binary Love relations (between the tidal deformabilities of binary neutron stars) encodes the baryon density at which the speed of sound rapidly changes. Twin stars lead to relations that present a signature "hill," "drop," and "swoosh" due to the second (mass-radius) stable branch, requiring a new description of the binary Love relations. Together, these features can reveal new properties and phases of nuclear matter.

Compact Objects in Alternative Gravities

We address neutron stars and black holes in alternative gravities, after recalling their basic properties in General Relativity. Among the plethora of interesting alternative gravities we here focus on an interesting set of scalar-tensor theories. We discuss the phenomenon of spontaneous scalarization, that is, matter induced for neutron stars and curvature induced for black holes. Along with other relevant physical properties, we address the quasi-normal modes of these compact objects. In particular, we consider universal relations of neutron stars to largely reduce the dependence on the equation of state, and we briefly address the shadow of black holes.

Translating Neutron Star Observations to Nuclear Symmetry Energy via Deep Neural Networks

One of the most significant challenges involved in efforts to understand the equation of state of dense neutron-rich matter is the uncertain density dependence of the nuclear symmetry energy. In particular, the nuclear symmetry energy is still rather poorly constrained, especially at high densities. On the other hand, detailed knowledge of the equation of state is critical for our understanding of many important phenomena in the nuclear terrestrial laboratories and the cosmos. Because of its broad impact, pinning down the density dependence of the nuclear symmetry energy has been a long-standing goal of both nuclear physics and astrophysics. Recent observations of neutron stars, in both electromagnetic and gravitational-wave spectra, have already constrained significantly the nuclear symmetry energy at high densities. The next generation of telescopes and gravitational-wave observatories will provide an unprecedented wealth of detailed observations of neutron stars, which will improve further our knowledge of the density dependence of nuclear symmetry energy, and the underlying equation of state of dense neutron-rich matter. Training deep neural networks to learn a computationally efficient representation of the mapping between astrophysical observables of neutron stars, such as masses, radii, and tidal deformabilities, and the nuclear symmetry energy allows its density dependence to be determined reliably and accurately. In this work, we use a deep learning approach to determine the nuclear symmetry energy as a function of density directly from observational neutron star data. We show, for the first time, that artificial neural networks can precisely reconstruct the nuclear symmetry energy from a set of available neutron star observables, such as masses and radii as measured by, e.g., the NICER mission, or masses and tidal deformabilities as measured by the LIGO/VIRGO/KAGRA gravitational-wave detectors. These results demonstrate the potential of artificial neural networks to reconstruct the symmetry energy and the equation of state directly from neutron star observational data, and emphasize the importance of the deep learning approach in the era of multi-messenger astrophysics.

First Constraints on Nuclear Coupling of Axionlike Particles from the Binary Neutron Star Gravitational Wave Event GW170817

Light axion fields, if they exist, can be sourced by neutron stars due to their coupling to nuclear matter, and play a role in binary neutron star mergers. We report on a search for such axions by analyzing the gravitational waves from the binary neutron star inspiral GW170817. We find no evidence of axions in the sampled parameter space. The null result allows us to impose constraints on axions with masses below 10^{-11}  eV by excluding the ones with decay constants ranging from 1.6×10^{16} to 10^{18}  GeV at a 3σ confidence level. Our analysis provides the first constraints on axions from neutron star inspirals, and rules out a large region in parameter space that has not been probed by the existing experiments.

High-Order Multipole and Binary Love Number Universal Relations

Using a data set of approximately 2 million phenomenological equations of state consistent with observational constraints, we construct new equation-of-state-insensitive universal relations that exist between the multipolar tidal deformability parameters of neutron stars, Λl, for several high-order multipoles (l=5,6,7,8), and we consider finite-size effects of these high-order multipoles in waveform modeling. We also confirm the existence of a universal relation between the radius of the 1.4M⊙ NS, R1.4 and the reduced tidal parameter of the binary, Λ˜, and the chirp mass. We extend this relation to a large number of chirp masses and to the radii of isolated NSs of different mass M, RM. We find that there is an optimal value of M for every M such that the uncertainty in the estimate of RM is minimized when using the relation. We discuss the utility and implications of these relations for the upcoming LIGO O4 run and third-generation detectors.

Astrophysical and Theoretical Physics Implications from Multimessenger Neutron Star Observations

The Neutron Star Interior Composition Explorer (NICER) recently measured the mass and equatorial radius of the isolated neutron star PSR J0030+0451. We use these measurements to infer the moment of inertia, the quadrupole moment, and the surface eccentricity of an isolated neutron star for the first time, using relations between these quantities that are insensitive to the unknown equation of state of supranuclear matter. We also use these results to forecast the moment of inertia of neutron star A in the double pulsar binary J0737-3039, a quantity anticipated to be directly measured in the coming decade with radio observations. Combining this information with the measurement of the tidal Love number with LIGO/Virgo observations, we propose and implement the first theory-agnostic and equation-of-state-insensitive test of general relativity. Specializing these constraints to a particular modified theory, we find that consistency with general relativity places the most stringent constraint on gravitational parity violation to date, surpassing all other previously reported bounds by 7 orders of magnitude and opens the path for a future test of general relativity with multimessenger neutron star observations.

The I-Love-Q Relations for Superfluid Neutron Stars

The I-Love-Q relations are approximate equation-of-state independent relations that connect the moment of inertia, the spin-induced quadrupole moment, and the tidal deformability of neutron stars. In this paper, we study the I-Love-Q relations for superfluid neutron stars for a general relativistic two-fluid model: one fluid being the neutron superfluid and the other a conglomerate of all charged components. We study to what extent the two-fluid dynamics might affect the robustness of the I-Love-Q relations by using a simple two-component polytropic model and a relativistic mean field model with entrainment for the equation-of-state. Our results depend crucially on the spin ratio Ωn/Ωp between the angular velocities of the neutron superfluid and the normal component. We find that the I-Love-Q relations can still be satisfied to high accuracy for superfluid neutron stars as long as the two fluids are nearly co-rotating Ωn/Ωp≈1. However, the deviations from the I-Love-Q relations increase as the spin ratio deviates from unity. In particular, the deviation of the Q-Love relation can be as large as O(10%) if Ωn/Ωp differ from unity by a few tens of percent. As Ωn/Ωp≈1 is expected for realistic neutron stars, our results suggest that the two-fluid dynamics should not affect the accuracy of any gravitational waveform models for neutron star binaries that employ the relation to connect the spin-induced quadrupole moment and the tidal deformability.

A Gravitational-Wave Perspective on Neutron-Star Seismology

We provide a bird’s-eye view of neutron-star seismology, which aims to probe the extreme physics associated with these objects, in the context of gravitational-wave astronomy. Focussing on the fundamental mode of oscillation, which is an efficient gravitational-wave emitter, we consider the seismology aspects of a number of astrophysically relevant scenarios, ranging from transients (like pulsar glitches and magnetar flares), to the dynamics of tides in inspiralling compact binaries and the eventual merged object and instabilities acting in isolated, rapidly rotating, neutron stars. The aim is not to provide a thorough review, but rather to introduce (some of) the key ideas and highlight issues that need further attention.

Spinning Black Holes Fall in Love

The open question of whether a black hole can become tidally deformed by an external gravitational field has profound implications for fundamental physics, astrophysics, and gravitational-wave astronomy. Love tensors characterize the tidal deformability of compact objects such as astrophysical (Kerr) black holes under an external static tidal field. We prove that all Love tensors vanish identically for a Kerr black hole in the nonspinning limit or for an axisymmetric tidal perturbation. In contrast to this result, we show that Love tensors are generically nonzero for a spinning black hole. Specifically, to linear order in the Kerr black hole spin and the weak perturbing tidal field, we compute in closed form the Love tensors that couple the mass-type and current-type quadrupole moments to the electric-type and magnetic-type quadrupolar tidal fields. For a dimensionless spin ∼0.1, the nonvanishing quadrupolar Love tensors are ∼2×10^{-3}, thus showing that black holes are particularly "rigid" compact objects.

Combining Electromagnetic and Gravitational-Wave Constraints on Neutron-Star Masses and Radii

We perform a joint Bayesian inference of neutron-star mass and radius constraints based on GW170817, observations of quiescent low-mass x-ray binaries (QLMXBs), photospheric radius expansion x-ray bursting sources, and x-ray timing observations of J0030+0451. With this dataset, the form of the prior distribution still has an impact on the posterior mass-radius curves and equation of state (EOS), but this impact is smaller than recently obtained when considering QLMXBs alone. We analyze the consistency of the electromagnetic data by including an "intrinsic scattering" contribution to the uncertainties, and find only a slight broadening of the posteriors. This suggests that the gravitational-wave and electromagnetic observations of neutron-star structure are providing a consistent picture of the neutron-star mass-radius curve and the EOS.

Neutron Star Equation of State in Light of GW190814

The observation of gravitational waves from an asymmetric binary opens the possibility for heavy neutron stars, but these pose challenges to models of the neutron star equation of state. We construct heavy neutron stars by introducing nontrivial structure in the speed of sound sourced by deconfined QCD matter, which cannot be well recovered by spectral representations. Their moment of inertia, Love number, and quadrupole moment are very small, so a tenfold increase in sensitivity may be needed to test this possibility with gravitational waves, which is feasible with third generation detectors.

Green's functions for geophysics: a review

The Green's function (GF) method, which makes use of GFs, is an important and elegant tool for solving a given boundary-value problem for the differential equation from a real engineering or physical field. Under a concentrated source, the solution of a differential equation is called a GF, which is singular at the source location, yet is very fundamental and powerful. When looking at the GFs from different physical and/or engineering fields, i.e. assigning the involved functions to real physical/engineering quantities, the GFs can be scaled and applied to large-scale problems such as those involved in Earth sciences as well as to nano-scale problems associated with quantum nanostructures. GFs are ubiquitous and everywhere: they can describe heat, water pressure, fluid flow potential, electromagnetic (EM) and gravitational potentials, and the surface tension of soap film. In the undergraduate courses Mechanics of Solids and Structural Analysis, a GF is the simple influence line or singular function. Dropping a pebble in the pond, it is the circular ripple traveling on and on. It is the wave generated by a moving ship in the opening ocean or the atom vibrating on a nanoscale sheet induced by the atomic force microscopy. In Earth science, while various GFs have been derived, a comprehensive review is missing. Thus, this article provides a relatively complete review on GFs for geophysics. In section 1, the George Green's potential functions, GF definition, as well as related theorems and basic relations are briefly presented. In section 2, the boundary-value problems for elastic and viscoelastic materials are provided. Section 3 is on the GFs in full- and half-spaces (planes). The GFs of concentrated forces and dislocations in horizontally layered half-spaces (planes) are derived in section 4 in terms of both Cartesian and cylindrical systems of vector functions. The corresponding GFs in a self-gravitating and layered spherical Earth are presented in section 5 in terms of the spherical system of vector functions. The singularity and infinity associated with GFs in layered systems are analyzed in section 6 along with a brief review of various layer matrix methods. Various associated mathematical preliminaries are listed in appendix, along with the three sets of vector function systems. It should be further emphasized that, while this review is targeted at geophysics, most of the GFs and solution methods can be equally applied to other engineering and science fields. Actually, many GFs and solutions methods reviewed in this article are derived by engineers and scientists from allied fields besides geophysics. As such, the updated approaches of constructing and deriving the GFs reviewed here should be very beneficial to any reader.

Tests of General Relativity with GW170817

The recent discovery by Advanced LIGO and Advanced Virgo of a gravitational wave signal from a binary neutron star inspiral has enabled tests of general relativity (GR) with this new type of source. This source, for the first time, permits tests of strong-field dynamics of compact binaries in the presence of matter. In this Letter, we place constraints on the dipole radiation and possible deviations from GR in the post-Newtonian coefficients that govern the inspiral regime. Bounds on modified dispersion of gravitational waves are obtained; in combination with information from the observed electromagnetic counterpart we can also constrain effects due to large extra dimensions. Finally, the polarization content of the gravitational wave signal is studied. The results of all tests performed here show good agreement with GR.

GW170817: Measurements of Neutron Star Radii and Equation of State

On 17 August 2017, the LIGO and Virgo observatories made the first direct detection of gravitational waves from the coalescence of a neutron star binary system. The detection of this gravitational-wave signal, GW170817, offers a novel opportunity to directly probe the properties of matter at the extreme conditions found in the interior of these stars. The initial, minimal-assumption analysis of the LIGO and Virgo data placed constraints on the tidal effects of the coalescing bodies, which were then translated to constraints on neutron star radii. Here, we expand upon previous analyses by working under the hypothesis that both bodies were neutron stars that are described by the same equation of state and have spins within the range observed in Galactic binary neutron stars. Our analysis employs two methods: the use of equation-of-state-insensitive relations between various macroscopic properties of the neutron stars and the use of an efficient parametrization of the defining function p(ρ) of the equation of state itself. From the LIGO and Virgo data alone and the first method, we measure the two neutron star radii as R_{1}=10.8_{-1.7}^{+2.0}  km for the heavier star and R_{2}=10.7_{-1.5}^{+2.1}  km for the lighter star at the 90% credible level. If we additionally require that the equation of state supports neutron stars with masses larger than 1.97  M_{⊙} as required from electromagnetic observations and employ the equation-of-state parametrization, we further constrain R_{1}=11.9_{-1.4}^{+1.4}  km and R_{2}=11.9_{-1.4}^{+1.4}  km at the 90% credible level. Finally, we obtain constraints on p(ρ) at supranuclear densities, with pressure at twice nuclear saturation density measured at 3.5_{-1.7}^{+2.7}×10^{34}  dyn cm^{-2} at the 90% level.

Tidal Love numbers of neutron stars in f(R) gravity

The recent detection of gravitational waves from a neutron star merger was a significant step towards constraining the nuclear matter equation of state by using the tidal Love numbers (TLNs) of the merging neutron stars. Measuring or constraining the neutron star TLNs allows us in principle to exclude or constraint many equations of state. This approach, however, has the drawback that many modified theories of gravity could produce deviations from General Relativity similar to the deviations coming from the uncertainties in the equation of state. The first and the most natural step in resolving the mentioned problem is to quantify the effects on the TLNs from the modifications of General Relativity. With this motivation in mind, in the present paper we calculate the TLNs of (non-rotating) neutron stars in R 2 -gravity. More precisely, by solving numerically the perturbation equations, we calculate explicitly the polar and the axial l = 2 TLNs for three characteristic realistic equations of state and compare the results to General Relativity. Our results show that while the polar TLNs are slightly influenced by the R 2 modification of General Relativity, the axial TLNs can be several times larger (in terms of the absolute value) compared to the general relativistic case.

Rotating stars in relativity

Rotating relativistic stars have been studied extensively in recent years, both theoretically and observationally, because of the information they might yield about the equation of state of matter at extremely high densities and because they are considered to be promising sources of gravitational waves. The latest theoretical understanding of rotating stars in relativity is reviewed in this updated article. The sections on equilibrium properties and on nonaxisymmetric oscillations and instabilities in f-modes and r-modes have been updated. Several new sections have been added on equilibria in modified theories of gravity, approximate universal relationships, the one-arm spiral instability, on analytic solutions for the exterior spacetime, rotating stars in LMXBs, rotating strange stars, and on rotating stars in numerical relativity including both hydrodynamic and magnetohydrodynamic studies of these objects.

No-hair theorem for black holes in astrophysical environments

According to the no-hair theorem, static black holes are described by a Schwarzschild spacetime provided there are no other sources of the gravitational field. This requirement, however, is in astrophysical realistic scenarios often violated, e.g., if the black hole is part of a binary system or if it is surrounded by an accretion disk. In these cases, the black hole is distorted due to tidal forces. Nonetheless, the subsequent formulation of the no-hair theorem holds: The contribution of the distorted black hole to the multipole moments that describe the gravitational field close to infinity and, thus, all sources is that of a Schwarzschild black hole. It still has no hair. This implies that there is no multipole moment induced in the black hole and that its second Love numbers, which measure some aspects of the distortion, vanish as was already shown in approximations to general relativity. But here we prove this property for astrophysical relevant black holes in full general relativity.

The Confrontation between General Relativity and Experiment

The status of experimental tests of general relativity and of theoretical frameworks for analyzing them is reviewed and updated. Einstein's equivalence principle (EEP) is well supported by experiments such as the Eötvös experiment, tests of local Lorentz invariance and clock experiments. Ongoing tests of EEP and of the inverse square law are searching for new interactions arising from unification or quantum gravity. Tests of general relativity at the post-Newtonian level have reached high precision, including the light deflection, the Shapiro time delay, the perihelion advance of Mercury, the Nordtvedt effect in lunar motion, and frame-dragging. Gravitational wave damping has been detected in an amount that agrees with general relativity to better than half a percent using the Hulse-Taylor binary pulsar, and a growing family of other binary pulsar systems is yielding new tests, especially of strong-field effects. Current and future tests of relativity will center on strong gravity and gravitational waves.

Effectively universal behavior of rotating neutron stars in general relativity makes them even simpler than their Newtonian counterparts

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 λ (the I-Love-Q relations), independently of the equation of state of the compact object. In the present Letter 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.