We consider the scattering of dark matter particles from superfluid liquid $^4$He, which has been proposed as a target for their direct detection. Focusing on dark matter masses below ~1 MeV, we demonstrate from sum-rule arguments that for momentum transfers of interest the dominant process is the creation of a single phonon, with the direct creation of two or more phonons much less likely. We show further that the anomalous dispersion of phonons in liquid $^4$He at low pressures [i.e., $d^2\omega (q)/dq^2 > 0$, where $q$ and $\omega(q)$ are the phonon momentum and energy] has the important consequence that a single phonon will decay over a relatively short distance into a shower of lower energy phonons centered on the direction of the original phonon. Thus the experimental challenge in this regime is to detect a shower of low energy phonons, not just a single phonon. Additional information from the distribution of phonons in such a shower could enhance the determination of the dark matter mass.

We present a much improved equation of state for neutron star matter, QHC19, with a smooth crossover from the hadronic regime at lower densities to the quark regime at higher densities. We now use the Togashi et al.~equation of state (Togashi:2017), a generalization of the Akmal-Pandharipande-Ravenhall equation of state of uniform nuclear matter, in the entire hadronic regime; the Togashi equation of state consistently describes non-uniform as well as uniform matter, and matter at beta equilibrium without the need for an interpolation between pure neutron and symmetric nuclear matter. We describe the quark matter regime at higher densities with the Nambu--Jona--Lasinio model, now identifying tight constraints on the phenomenological universal vector repulsion between quarks and the pairing interaction between quarks arising from the requirements of thermodynamic stability and causal propagation of sound. The resultant neutron star properties agree very well with the inferences of the LIGO/Virgo collaboration, from GW170817, of the pressure vs. baryon density, neutron star radii, and tidal deformabilities. The maximum neutron star mass allowed by QHC19 is 2.35 $M_\odot$, consistent with all neutron star mass determinations.

A moderately strong vector repulsion between quarks in dense quark matter is needed to explain how a quark core can support neutron stars heavier than two solar masses. We study this repulsion, parametrized by a four-fermion interaction with coupling g_V, in terms of non-perturbative gluon exchange in QCD in the Landau gauge. Matching the energy of quark matter, g_V n_q^2 (where n_q is the number density of quarks) with the quark exchange energy calculated in QCD with a gluon propagator parametrized by a finite gluon mass m_g and a frozen coupling alpha_s, at moderate quark densities, we find that gluon masses m_g in the range 200 - 600 MeV and alpha_s = 2 - 4 lead to a g_V consistent with neutron star phenomenology. Estimating the effects of quark masses and a color-flavor-locked (CFL) pairing gap, we find that g_V can be well approximated by a flavor-symmetric, decreasing function of density. We briefly discuss similar matchings for the isovector repulsion and for the pairing attraction.

As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations.

Neutron stars were first posited in the early thirties, and discovered as pulsars in late sixties; however, only recently are we beginning to understand the matter they contain. This talk describes the continuing development of a consistent picture of the liquid interiors of neutron stars, driven by four advances: observations of two heavy neutron stars with masses $\simeq$ 2.0 solar masses; inferences of masses and radii simultaneously for an increasing number of neutron stars in low mass X-ray binaries, and future determinations via the NICER observatory; the observation of the binary neutron star merger, GW170817, through gravitational waves as well as across the electromagnetic spectrum; and an emerging understanding in QCD of how nuclear matter can turn into deconfined quark matter in the interior. We describe the modern quark-hadron crossover equation of state, QHC18, and the corresponding neutron stars, which agree well with current observations.

We demonstrate that the observation of neutron stars with masses greater than one solar mass places severe demands on any exotic neutron decay mode that could explain the discrepancy between beam and bottle measurements of the neutron lifetime. If the neutron can decay to a stable, feebly-interacting dark fermion, the maximum possible mass of a neutron star is 0.7 solar masses, while all well-measured neutron star masses exceed one solar mass. The survival of $2 M_\odot$ neutron stars therefore indicates that any explanation beyond the Standard Model for the neutron lifetime puzzle requires dark matter to be part of a multi-particle dark sector with highly constrained interactions.

We review the equation of state of matter in neutron stars from the solid crust through the liquid nuclear matter interior to the quark regime at higher densities. We focus in detail on the question of how quark matter appears in neutron stars, and how it affects the equation of state. After discussing the crust and liquid nuclear matter in the core we briefly review aspects of microscopic quark physics relevant to neutron stars, and quark models of dense matter based on the Nambu--Jona-Lasinio framework, in which gluonic processes are replaced by effective quark interactions. We turn then to describing equations of state useful for interpretation of both electromagnetic and gravitational observations, reviewing the emerging picture of hadron-quark continuity in which hadronic matter turns relatively smoothly, with at most only a weak first order transition, into quark matter with increasing density. We review construction of unified equations of state that interpolate between the reasonably well understood nuclear matter regime at low densities and the quark matter regime at higher densities. The utility of such interpolations is driven by the present inability to calculate the dense matter equation of state in QCD from first principles. As we review, the parameters of effective quark models -- which have direct relevance to the more general structure of the QCD phase diagram of dense and hot matter -- are constrained by neutron star mass and radii measurements, in particular favoring large repulsive density-density and attractive diquark pairing interactions. We describe the structure of neutron stars constructed from the unified equations of states with crossover. Lastly we present the current equations of state -- called "QHC18" for quark-hadron crossover -- in a parametrized form practical for neutron star modeling.

We develop a unified description, via the Boltzmann equation, of damping of gravitational waves by matter, incorporating collisions. We identify two physically distinct damping mechanisms -- collisional and Landau damping. We first consider damping in flat spacetime, and then generalize the results to allow for cosmological expansion. In the first regime, maximal collisional damping of a gravitational wave, independent of the details of the collisions in the matter is, as we show, significant only when its wavelength is comparable to the size of the horizon. Thus damping by intergalactic or interstellar matter for all but primordial gravitational radiation can be neglected. Although collisions in matter lead to a shear viscosity, they also act to erase anisotropic stresses, thus suppressing the damping of gravitational waves. Damping of primordial gravitational waves remains possible. We generalize Weinberg's calculation of gravitational wave damping, now including collisions and particles of finite mass, and interpret the collisionless limit in terms of Landau damping. While Landau damping of gravitational waves cannot occur in flat spacetime, the expansion of the universe allows such damping by spreading the frequency of a gravitational wave of given wavevector.

We delineate the properties of QCD matter at baryon density n_B=1-10n_0 (n_0: nuclear saturation density), through the construction of neutron star equations of state that satisfy the neutron star mass-radius constraints as well as physical conditions on the speed of sound. The QCD matter is described in the 3-window modeling: at n_B < 2n_0 purely nuclear matter; at n_B > 5n_0 percolated quark matter; and at 2n_0 < n_B < 5n_0 matter intermediate between these two which are constructed by interpolation. Using a schematic quark model with effective interactions inspired from hadron and nuclear physics, we analyze the strength of interactions necessary to describe observed neutron star properties. Our finding is that the interactions should remain as strong as in the QCD vacuum, indicating that gluons at n_B =1-10 n_0 remain non-perturbative even after quark matter formation.

We present a comprehensive study of spectroscopic radius measurements of twelve neutron stars obtained during thermonuclear bursts or in quiescence. We incorporate, for the first time, a large number of systematic uncertainties in the measurement of the apparent angular sizes, Eddington fluxes, and distances, in the composition of the interstellar medium, and in the flux calibration of X-ray detectors. We also take into account the results of recent theoretical calculations of rotational effects on neutron star radii, of atmospheric effects on surface spectra, and of relativistic corrections to the Eddington critical flux. We employ Bayesian statistical frameworks to obtain neutron star radii from the spectroscopic measurements as well as to infer the equation of state from the radius measurements. Combining these with the results of experiments in the vicinity of nuclear saturation density and the observations of ~2 Msun neutron stars, we place strong and quantitative constraints on the properties of the equation of state between approximately 2-8 times the nuclear saturation density. We find that around M=1.5 Msun, the preferred equation of state predicts radii between 10.1 - 11.1 km. When interpreting the pressure constraints in the context of high density equations of state based on interacting nucleons, our results suggest a relatively weak contribution of the three-body interaction potential.

We construct an equation of state for massive neutron stars based on quantum chromodynamics phenomenology. Our primary purpose is to delineate the relevant ingredients of equations of state that simultaneously have the required stiffness and satisfy constraints from thermodynamics and causality. These ingredients are: (i) a repulsive density-density interaction, universal for all flavors; (ii) the color-magnetic interaction active from low to high densities; (iii) confining effects, which become increasingly important as the baryon density decreases; (iv) non-perturbative gluons, which are not very sensitive to changes of the quark density. We use the following "3-window" description: At baryon densities below about twice normal nuclear density, 2n_0, we use the Akmal-Pandharipande-Ravenhall (APR) equation of state, and at high densities, > (4-7)n_0, we use the three-flavor Nambu-Jona-Lasinio (NJL) model supplemented by vector and diquark interactions. In the transition density region, we smoothly interpolate the hadronic and quark equations of state in the chemical potential-pressure plane. Requiring that the equation of state approach APR at low densities, we find that the quark pressure in non-confining models can be larger than the hadronic pressure, unlike in conventional equations of state. We show that consistent equations of state of stiffness sufficient to allow massive neutron stars are reasonably tightly constrained, suggesting that gluon dynamics remains non-perturbative even at baryon densities ~10n_0.

We review the history of neutron star physics in the 1930s that is related to L. Landau. According to recollections of Rosenfeld (1974, Proc. 16th Solvay Conference on Physics, p. 174), Landau improvised the concept of neutron stars in a discussion with Bohr and Rosenfeld just after the news of the discovery of the neutron reached Copenhagen in February 1932. We present arguments that the discussion took place in March 1931, before the discovery of the neutron, and that they in fact discussed the paper written by Landau in Zurich in February 1931 but not published until February 1932 (Phys. Z. Sowjetunion, 1, 285). In his paper Landau mentioned the possible existence of dense stars which look like one giant nucleus; this can be regarded as an early theoretical prediction or anticipation of neutron stars, prior to the discovery of the neutron. The coincidence of the dates of the neutron's discovery and the paper's publication has led to an erroneous association of the paper with the discovery of the neutron. In passing, we outline the contribution of Landau to the theory of white dwarfs and to the hypothesis of stars with neutron cores.

We present the first astrophysical measurement of the pressure of cold matter above nuclear saturation density, based on recently determined masses and radii of three neutron stars. The pressure at higher densities are below the predictions of equations of state that account only for nucleonic degrees of freedom, and thus present a challenge to the microscopic theory of neutron star matter.

Recent observations of neutron star masses close to the maximum predicted by nucleonic equations of state begin to challenge our understanding of dense matter in neutron stars, and constrain the possible presence of quark matter in their deep interiors.

This short encyclopedia article, reviewing current information on neutron stars, is intended for a broad scientific audience.

We calculate the electrical conductivity in the early universe at temperatures below as well as above the electroweak vacuum scale, $T_c\simeq 100$GeV. Debye and dynamical screening of electric and magnetic interactions leads to a finite conductivity, $\sigma_{el}\sim T/\alpha\ln(1/\alpha)$, at temperatures well below $T_c$. At temperatures above, $W^\pm$ charge-exchange processes -- analogous to color exchange through gluons in QCD -- effectively stop left-handed charged leptons. However, right-handed leptons can carry current, resulting in $\sigma_{el}/T$ being only a factor $\sim \cos^4\theta_W$ smaller than at temperatures below $T_c$.

Observational identification of black holes as members of binary systems requires the knowledge of the upper limit on the gravitational mass of a neutron star. We use modern equations of state for neutron star matter, fitted to experimental nucleon-nucleon scattering data and the properties of light nuclei, to calculate, within the framework of Rhoades & Ruffini (1974), the minimum upper limit on a neutron star mass. Regarding the equation of state as valid up to twice nuclear matter saturation density, rho_{nm}, we obtain a secure upper bound on the neutron star mass equal to 2.9 solar masses. We also find that in order to reach the lowest possible upper bound of 2.2 solar masses, we need understand the physical properties of neutron matter up to a density of about 4 times rho_{nm}.