Tidal Deformabilities and Radii of Neutron Stars from the Observation of GW170817

Minimum Neutron Star Mass in Neutrino-Driven Supernova Explosions

Supernova theory has struggled to explain the lightest known neutron star candidate with an accurate mass determination, the 1.174M_{⊙} companion in the eccentric compact binary system J0453+1559. To improve the theoretical lower limit for neutron star birth masses, we perform 3D supernova simulations for five stellar models close to the minimum mass for iron core collapse. We obtain a record-low neutron star mass of 1.192M_{⊙} and a substantial kick of ∼100  km s^{-1}. Given residual uncertainties in stellar evolution, a neutron star origin for the 1.174M_{⊙} object remains plausible.

Self-Consistent Picture of the Mass Ejection from a One Second Long Binary Neutron Star Merger Leaving a Short-Lived Remnant in a General-Relativistic Neutrino-Radiation Magnetohydrodynamic Simulation

We perform a general-relativistic neutrino-radiation magnetohydrodynamic simulation of a one second-long binary neutron star merger on the Japanese supercomputer Fugaku using about 85 million CPU hours with 20 736 CPUs. We consider an asymmetric binary neutron star merger with masses of 1.2M_{⊙} and 1.5M_{⊙} and a "soft" equation of state SFHo. It results in a short-lived remnant with the lifetime of ≈0.017  s, and subsequent massive torus formation with the mass of ≈0.05M_{⊙} after the remnant collapses to a black hole. For the first time, we find that after the dynamical mass ejection, which drives the fast tail and mildly relativistic components, the postmerger mass ejection from the massive torus takes place due to the magnetorotational instability-driven turbulent viscosity in a single simulation and the two ejecta components are seen in the distributions of the electron fraction and velocity with distinct features.

Degeneracy in the Inference of Phase Transitions in the Neutron Star Equation of State from Gravitational Wave Data

Gravitational wave (GW) detections of binary neutron star inspirals will be crucial for constraining the dense matter equation of state (EOS). We demonstrate a new degeneracy in the mapping from tidal deformability data to the EOS, which occurs for models with strong phase transitions. We find that there exists a new family of EOS with phase transitions that set in at different densities and that predict neutron star radii that differ by up to ∼500  m but that produce nearly identical tidal deformabilities for all neutron star masses. Next-generation GW detectors and advances in nuclear theory may be needed to resolve this degeneracy.

Probing the Incompressibility of Nuclear Matter at Ultrahigh Density through the Prompt Collapse of Asymmetric Neutron Star Binaries

Using 250 neutron star merger simulations with microphysics, we explore for the first time the role of nuclear incompressibility in the prompt collapse threshold for binaries with different mass ratios. We demonstrate that observations of prompt collapse thresholds, either from binaries with two different mass ratios or with one mass ratio but combined with the knowledge of the maximum neutron star mass or compactness, will constrain the incompressibility at the maximum neutron star density K_{max} to within tens of percent. This otherwise inaccessible measure of K_{max} can potentially reveal the presence of hyperons or quarks inside neutron stars.

Temperature Effects on Core g-Modes of Neutron Stars

Neutron stars provide a unique physical laboratory in which to study the properties of matter at high density and temperature. We study a diagnostic of the composition of high-density matter, namely, g-mode oscillations, which are driven by buoyancy forces. These oscillations can be excited by tidal forces and couple to gravitational waves. We extend prior results for the g-mode spectrum of cold neutron star matter to high temperatures that are expected to be achieved in neutron star mergers using a parameterization for finite-temperature effects on equations of state recently proposed by Raithel, Özel and Psaltis. We find that the g-modes of canonical mass neutron stars (≈1.4M⊙) are suppressed at high temperatures, and core g-modes are supported only in the most massive (≥2M⊙) of hot neutron stars.

Ranking Love Numbers for the Neutron Star Equation of State: The Need for Third-Generation Detectors

Gravitational-wave measurements of the tidal deformability in neutron-star binary coalescences can be used to infer the still unknown equation of state (EOS) of dense matter above the nuclear saturation density. By employing a Bayesian-ranking test, we quantify the ability of current and future gravitational-wave observations to discriminate among families of nuclear-physics based EOS which differ in particle content and ab initio microscopic calculations. While the constraining power of GW170817 is limited, we show that even twenty coalescences detected by LIGO-Virgo at design sensitivity are not enough to discriminate between EOS with similar softness but distinct microphysics. However, just a single detection with a third-generation detector such as the Einstein Telescope or Cosmic Explorer will rule out several families of EOS with very strong statistical significance and can discriminate among models which feature similar softness, hence, constraining the properties of nuclear matter to unprecedented levels.

The Equation of State of Neutron-Rich Matter at Fourth Order of Chiral Effective Field Theory and the Radius of a Medium-Mass Neutron Star

We report neutron star predictions based on our most recent equations of state. These are derived from chiral effective field theory, which allows for a systematic development of nuclear forces, order by order. We utilize high-quality two-nucleon interactions and include all three-nucleon forces up to fourth order in the chiral expansion. Our ab initio predictions are restricted to the domain of applicability of chiral effective field theory. However, stellar matter in the interior of neutron stars can be up to several times denser than normal nuclear matter at saturation, and its composition is essentially unknown. Following established practices, we extend our microscopic predictions to higher densities matching piecewise polytropes. The radius of the average-size neutron star, about 1.4 solar masses, is sensitive to the pressure at normal densities, and thus it is suitable to constrain ab initio theories of the equation of state. For this reason, we focus on the radius of medium-mass stars. We compare our results with other theoretical predictions and recent constraints.

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.

Nuclear Physics and Astrophysics Constraints on the High Density Matter Equation of State

(1) This review has been written in memory of Steven Moszkowski who unexpectedly passed away in December 2020. It has been inspired by our many years of discussions. Steven’s enthusiasm, drive and determination to understand atomic nuclei in simple terms of basic laws of physics was infectious. He sought the fundamental origin of nuclear forces in free space, and their saturation and modification in nuclear medium. His untimely departure left our job unfinished but his legacy lives on. (2) Focusing on the nuclear force acting in nuclear matter of astrophysical interest and its equation of state (EoS), we take several typical snapshots of evolution of the theory of nuclear forces. We start from original ideas in the 1930s moving through to its overwhelming diversity today. The development is supported by modern observational and terrestrial data and their inference in the multimessenger era, as well as by novel mathematical techniques and computer power. (3) We find that, despite the admirable effort both in theory and measurement, we are facing multiple models dependent on a large number of variable correlated parameters which cannot be constrained by data, which are not yet accurate, nor sensitive enough, to identify the theory closest to reality. The role of microphysics in the theories is severely limited or neglected, mostly deemed to be too difficult to tackle. (4) Taking the EoS of high-density matter as an example, we propose to develop models, based, as much as currently possible, on the microphysics of the nuclear force, with a minimal set of parameters, chosen under clear physical guidance. Still somewhat phenomenological, such models could pave the way to realistic predictions, not tracing the measurement, but leading it.

Progress in Constraining Nuclear Symmetry Energy Using Neutron Star Observables Since GW170817

The density dependence of nuclear symmetry energy is among the most uncertain parts of the Equation of State (EOS) of dense neutron-rich nuclear matter. It is currently poorly known especially at suprasaturation densities partially because of our poor knowledge about isovector nuclear interactions at short distances. Because of its broad impacts on many interesting issues, pinning down the density dependence of nuclear symmetry energy has been a longstanding and shared goal of both astrophysics and nuclear physics. New observational data of neutron stars including their masses, radii, and tidal deformations since GW170817 have helped improve our knowledge about nuclear symmetry energy, especially at high densities. Based on various model analyses of these new data by many people in the nuclear astrophysics community, while our brief review might be incomplete and biased unintentionally, we learned in particular the following: (1) The slope parameter L of nuclear symmetry energy at saturation density ρ0 of nuclear matter from 24 new analyses of neutron star observables was about L≈57.7±19 MeV at a 68% confidence level, consistent with its fiducial value from surveys of over 50 earlier analyses of both terrestrial and astrophysical data within error bars. (2) The curvature Ksym of nuclear symmetry energy at ρ0 from 16 new analyses of neutron star observables was about Ksym≈−107±88 MeV at a 68% confidence level, in very good agreement with the systematics of earlier analyses. (3) The magnitude of nuclear symmetry energy at 2ρ0, i.e., Esym(2ρ0)≈51±13 MeV at a 68% confidence level, was extracted from nine new analyses of neutron star observables, consistent with the results from earlier analyses of heavy-ion reactions and the latest predictions of the state-of-the-art nuclear many-body theories. (4) While the available data from canonical neutron stars did not provide tight constraints on nuclear symmetry energy at densities above about 2ρ0, the lower radius boundary R2.01=12.2 km from NICER’s very recent observation of PSR J0740+6620 of mass 2.08±0.07M⊙ and radius R=12.2–16.3 km at a 68% confidence level set a tight lower limit for nuclear symmetry energy at densities above 2ρ0. (5) Bayesian inferences of nuclear symmetry energy using models encapsulating a first-order hadron–quark phase transition from observables of canonical neutron stars indicated that the phase transition shifted appreciably both L and Ksym to higher values, but with larger uncertainties compared to analyses assuming no such phase transition. (6) The high-density behavior of nuclear symmetry energy significantly affected the minimum frequency necessary to rotationally support GW190814’s secondary component of mass (2.50–2.67) M⊙ as the fastest and most massive pulsar discovered so far. Overall, thanks to the hard work of many people in the astrophysics and nuclear physics community, new data of neutron star observations since the discovery of GW170817 have significantly enriched our knowledge about the symmetry energy of dense neutron-rich nuclear matter.

Constraining Dense Matter Physics Using f-Mode Oscillations in Neutron Stars

In this paper, an investigation of the role of nuclear saturation parameters on f-mode oscillations in neutron stars is performed within the Cowling approximation. It is found that the uncertainty in the effective nucleon mass plays a dominant role in controlling the f-mode frequencies. The effect of the uncertainties in saturation parameters on previously-proposed empirical relations of the frequencies with astrophysical observables relevant for asteroseismology are also investigated. These results can serve as an important tool for constraining the nuclear parameter space and understand the behaviour of dense nuclear matter from the future detection of f-modes.

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.

Multimessenger Binary Mergers Containing Neutron Stars: Gravitational Waves, Jets, and γ-Ray Bursts

Neutron stars (NSs) are extraordinary not only because they are the densest form of matter in the visible Universe but also because they can generate magnetic fields ten orders of magnitude larger than those currently constructed on earth. The combination of extreme gravity with the enormous electromagnetic (EM) fields gives rise to spectacular phenomena like those observed on August 2017 with the merger of a binary neutron star system, an event that generated a gravitational wave (GW) signal, a short γ -ray burst (sGRB), and a kilonova. This event serves as the highlight so far of the era of multimessenger astronomy. In this review, we present the current state of our theoretical understanding of compact binary mergers containing NSs as gleaned from the latest general relativistic magnetohydrodynamic simulations. Such mergers can lead to events like the one on August 2017, GW170817, and its EM counterparts, GRB 170817 and AT 2017gfo. In addition to exploring the GW emission from binary black hole-neutron star and neutron star-neutron star mergers, we also focus on their counterpart EM signals. In particular, we are interested in identifying the conditions under which a relativistic jet can be launched following these mergers. Such a jet is an essential feature of most sGRB models and provides the main conduit of energy from the central object to the outer radiation regions. Jet properties, including their lifetimes and Poynting luminosities, the effects of the initial magnetic field geometries and spins of the coalescing NSs, as well as their governing equation of state, are discussed. Lastly, we present our current understanding of how the Blandford-Znajek mechanism arises from merger remnants as the trigger for launching jets, if, when and how a horizon is necessary for this mechanism, and the possibility that it can turn on in magnetized neutron ergostars, which contain ergoregions, but no horizons.

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.

Probing the Nuclear Equation of State from the Existence of a ∼2.6 M⊙ Neutron Star: The GW190814 Puzzle

On 14 August 2019, the LIGO/Virgo collaboration observed a compact object with mass ∼2.59−0.09+0.08M⊙, as a component of a system where the main companion was a black hole with mass ∼23M⊙. A scientific debate initiated concerning the identification of the low mass component, as it falls into the neutron star–black hole mass gap. The understanding of the nature of GW190814 event will offer rich information concerning open issues, the speed of sound and the possible phase transition into other degrees of freedom. In the present work, we made an effort to probe the nuclear equation of state along with the GW190814 event. Firstly, we examine possible constraints on the nuclear equation of state inferred from the consideration that the low mass companion is a slow or rapidly rotating neutron star. In this case, the role of the upper bounds on the speed of sound is revealed, in connection with the dense nuclear matter properties. Secondly, we systematically study the tidal deformability of a possible high mass candidate existing as an individual star or as a component one in a binary neutron star system. As the tidal deformability and radius are quantities very sensitive on the neutron star equation of state, they are excellent counters on dense matter properties. We conjecture that similar isolated neutron stars or systems may exist in the universe and their possible future observation will shed light on the maximum neutron star mass problem.

Isospin Symmetry Breaking Effects on the Mass-Radius Relation of a Neutron Star

Isospin symmetry breaking effects on the mass-radius relation of a cold, non-accreting neutron star are studied on the basis of two Skyrme Energy Density Functionals (EDFs). One functional contains isospin symmetry breaking terms other than those typically included in Skyrme EDFs while its counterpart is of standard form. Both functionals are based on the same fitting protocol except for the observables and pseudo-observables sensitive to the isospin symmetry breaking channel. The quality of those functionals is similar in the description of terrestrial observables but choosing either of them has a non-negligible effect on the mass-radius relation and tidal deformability of a neutron star. Further investigations are needed to clarify the effects of isospin symmetry breaking on these and other observables of neutron stars that are, and will become, available.

Inference of the Neutron Star Equation of State from Cosmological Distances

Finite-size effects on the gravitational wave signal from a neutron star merger typically manifest at high frequencies where detector sensitivity decreases. Proposed sensitivity improvements can give us access both to stronger signals and to a myriad of weak signals from cosmological distances. The latter will outnumber the former and the relevant part of the signal will be redshifted towards the detector's most sensitive band. We study the redshift dependence of information about neutron star matter and find that single-scale properties, such as the star radius or the postmerger frequency, are better measured from the distant weak sources from z∼1.

Neutron Star Properties: Quantifying the Effect of the Crust–Core Matching Procedure

The impact of the equation of state (EoS) crust-core matching procedure on neutron star (NS) properties is analyzed within a meta-modeling approach. Using a Taylor expansion to parametrize the core equation of state (EoS) and the SLy4 crust EoS, we create two distinct EoS datasets employing two matching procedures. Each EoS describes cold NS matter in a β equilibrium that is thermodynamically stable and causal. It is shown that the crust-core matching procedure affects not only the crust-core transition but also the nuclear matter parameter space of the core EoS, and thus the most probable nuclear matter properties. An uncertainty of as much as 5% (8%) on the determination of low mass NS radii (tidal deformability) is attributed to the complete matching procedure, including the effect on core EoS. By restricting the analysis, imposing that the same set of core EoS is retained in both matching procedures, the uncertainty on the NS radius drops to 3.5% and below 1.5% for 1.9M⊙. Moreover, under these conditions, the crust-core matching procedure has a strong impact on the Love number k2, of almost 20% for 1.0M⊙ stars and 7% for 1.9M⊙ stars, but it shows a very small impact on the tidal deformability Λ, below 1%.

Equation of State Constraints from the Threshold Binary Mass for Prompt Collapse of Neutron Star Mergers

Using hydrodynamical simulations for a large set of high-density matter equations of state (EOSs), we systematically determine the threshold mass M_{thres} for prompt black-hole formation in equal-mass and asymmetric neutron star (NS) mergers. We devise the so far most direct, general, and accurate method to determine the unknown maximum mass of nonrotating NSs from merger observations revealing M_{thres}. Considering hybrid EOSs with hadron-quark phase transition, we identify a new, observable signature of quark matter in NS mergers. Furthermore, our findings have direct applications in gravitational wave searches, kilonova interpretations, and multimessenger constraints on NS properties.

A lower bound on the maximum mass if the secondary in GW190814 was once a rapidly spinning neutron star

The recent detection of GW190814 featured the merger of a binary with a primary having a mass of $\sim 23\, \mathrm{ M}_{\odot }$ and a secondary with a mass of $\sim 2.6\, \mathrm{ M}_{\odot }$. While the primary was most likely a black hole, the secondary could be interpreted as either the lightest black hole or the most massive neutron star ever observed, but also as the indication of a novel class of exotic compact objects. We here argue that although the secondary in GW190814 is most likely a black hole at merger, it needs not be an ab-initio black hole nor an exotic object. Rather, based on our current understanding of the nuclear-matter equation of state, it can be a rapidly rotating neutron star that collapsed to a rotating black hole at some point before merger. Using universal relations connecting the masses and spins of uniformly rotating neutron stars, we estimate the spin, $0.49_{-0.05}^{+0.08} \lesssim \chi \lesssim 0.68_{-0.05}^{+0.11}$, of the secondary – a quantity not constrained so far by the detection – and a novel strict lower bound on the maximum mass, $M_{_{\mathrm{TOV}}}\gt 2.08^{+0.04}_{-0.04}\, \, \mathrm{ M}_{\odot }$ and an optimal bound of $M_{_{\mathrm{TOV}}}\gt 2.15^{+0.04}_{-0.04}\, \, \mathrm{ M}_{\odot }$, of non-rotating neutron stars, consistent with recent observations of a very massive pulsar. The new lower bound also remains valid even in the less likely scenario in which the secondary neutron star never collapsed to a black hole.

Fast Rotating Relativistic Stars: Spectra and Stability without Approximation

We study oscillations and instabilities of relativistic stars using perturbation theory in general relativity and take into account the contribution of a dynamic spacetime. We present the oscillation spectrum as well as the critical values for the onset of the secular CFS instability of neutron stars, and propose universal relations for gravitational wave asteroseismology, which may help constrain the neutron star radius and/or the nuclear equation of state. The results are relevant for all stages during a neutron star's life but especially to nascent or remnant objects following a binary merger.

The Equation of State of Nuclear Matter: From Finite Nuclei to Neutron Stars

Background. We investigate possible correlations between neutron star observables and properties of atomic nuclei. In particular, we explore how the tidal deformability of a 1.4 solar mass neutron star, M1.4, and the neutron-skin thickness of 48Ca and 208Pb are related to the stellar radius and the stiffness of the symmetry energy. Methods. We examine a large set of nuclear equations of state based on phenomenological models (Skyrme, NLWM, DDM) and ab initio theoretical methods (BBG, Dirac–Brueckner, Variational, Quantum Monte Carlo). Results: We find strong correlations between tidal deformability and NS radius, whereas a weaker correlation does exist with the stiffness of the symmetry energy. Regarding the neutron-skin thickness, weak correlations appear both with the stiffness of the symmetry energy, and the radius of a M1.4. Our results show that whereas the considered EoS are compatible with the largest masses observed up to now, only five microscopic models and four Skyrme forces are simultaneously compatible with the present constraints on L and the PREX experimental data on the 208Pb neutron-skin thickness. We find that all the NLWM and DDM models and the majority of the Skyrme forces are excluded by these two experimental constraints, and that the analysis of the data collected by the NICER mission excludes most of the NLWM considered. Conclusion. The tidal deformability of a M1.4 and the neutron-skin thickness of atomic nuclei show some degree of correlation with nuclear and astrophysical observables, which however depends on the ensemble of adopted EoS.

Binary Neutron Star Merger Simulations with a Calibrated Turbulence Model

Magnetohydrodynamic (MHD) turbulence in neutron star (NS) merger remnants can impact their evolution and multi-messenger signatures, complicating the interpretation of present and future observations. Due to the high Reynolds numbers and the large computational costs of numerical relativity simulations, resolving all the relevant scales of the turbulence will be impossible for the foreseeable future. Here, we adopt a method to include subgrid-scale turbulence in moderate resolution simulations by extending the large-eddy simulation (LES) method to general relativity (GR). We calibrate our subgrid turbulence model with results from very-high-resolution GRMHD simulations, and we use it to perform NS merger simulations and study the impact of turbulence. We find that turbulence has a quantitative, but not qualitative, impact on the evolution of NS merger remnants, on their gravitational wave signatures, and on the outflows generated in binary NS mergers. Our approach provides a viable path to quantify uncertainties due to turbulence in NS mergers.

Kilonovae

The coalescence of double neutron star (NS-NS) and black hole (BH)-NS binaries are prime sources of gravitational waves (GW) for Advanced LIGO/Virgo and future ground-based detectors. Neutron-rich matter released from such events undergoes rapid neutron capture (r-process) nucleosynthesis as it decompresses into space, enriching our universe with rare heavy elements like gold and platinum. Radioactive decay of these unstable nuclei powers a rapidly evolving, approximately isotropic thermal transient known as a "kilonova", which probes the physical conditions during the merger and its aftermath. Here I review the history and physics of kilonovae, leading to the current paradigm of day-timescale emission at optical wavelengths from lanthanide-free components of the ejecta, followed by week-long emission with a spectral peak in the near-infrared (NIR). These theoretical predictions, as compiled in the original version of this review, were largely confirmed by the transient optical/NIR counterpart discovered to the first NS-NS merger, GW170817, discovered by LIGO/Virgo. Using a simple light curve model to illustrate the essential physical processes and their application to GW170817, I then introduce important variations about the standard picture which may be observable in future mergers. These include ∼ hour-long UV precursor emission, powered by the decay of free neutrons in the outermost ejecta layers or shock-heating of the ejecta by a delayed ultra-relativistic outflow; and enhancement of the luminosity from a long-lived central engine, such as an accreting BH or millisecond magnetar. Joint GW and kilonova observations of GW170817 and future events provide a new avenue to constrain the astrophysical origin of the r-process elements and the equation of state of dense nuclear matter.