We explore constraints on the equation of state of neutron-rich matter based on microscopic calculations up to nuclear densities and observations of neutron stars. In a previous work arXiv:1303.4662 we showed that predictions based on modern nuclear interactions derived within chiral effective field theory and the observation of 2-solar-mass neutron stars result in a robust uncertainty range for neutron star radii and the equation of state over a wide range of densities. In this work we extend this study, employing both the piecewise polytrope extension from arXiv:1303.4662 as well as the speed of sound model of arXiv:1812.08188, and show that moment of inertia measurements of neutron stars can significantly improve the constraints on the equation of state and neutron star radii.

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 consider the elastic constants of phases with nonspherical nuclei, so-called pasta phases, predicted to occur in the inner crust of a neutron star. First, we treat perfectly ordered phases and give numerical estimates for lasagna and spaghetti when the pasta elements are spatially uniform: the results are in order-of-magnitude agreement with the numerical simulations of Caplan, Schneider, and Horowitz, Phys. Rev. Lett. 121, 132701 (2018). We then turn to pasta phases without long-range order and calculate upper (Voigt) and lower (Reuss) bounds on the effective shear modulus and find that the lower bound is zero, but the upper bound is nonzero. To obtain better estimates, we then apply the self-consistent formalism and find that this predicts that the shear modulus of the phases without long-range order is zero if the pasta elements are spatially uniform. In numerical simulations, the pasta elements are found to be modulated spatially and we show that this modulation is crucial to obtaining a nonzero elastic moduli for pasta phases without long-range order. In the self-consistent formalism we find that, for lasagna, the effective shear modulus is linear in the elastic constants that do not vanish when the pasta elements are spatially uniform while, for spaghetti, it varies as the square root of these elastic constants. We also consider the behavior of the elastic constant associated with a homologous strain (hydrostatic compression) of the structure of the pasta phases without long-range order.

A number of properties of dense matter can be understood semiquantitatively in terms of simple physical arguments. We begin with the outer parts of neutron stars, and consider the density at which pressure ionization occurs, the density at which electrons become relativistic, the density at which neutrons drip out of nuclei, and the size of the equilibrium nucleus in dense matter. Subsequently, we treat the so-called "pasta" phases expected to occur at densities just below the density at which the transition from the crust to the liquid core of a neutron star occurs. We then consider aspects of superfluidity in dense matter. Estimates of pairing gaps in homogeneous nuclear matter are given, and the effect of the dense medium on the interaction between nucleons is described. Finally, we turn to superfluidity in the crust of neutron stars and especially the neutron superfluid density, an important quantity in the theory of sudden speedups of the rotation rate of some pulsars.

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.

Calculations of the effects of band structure on the neutron superfluid density in the crust of neutron stars made under the assumption that the effects of pairing are small [N. Chamel, Phys. Rev. C 85, 035801 (2012)] lead to moments of inertia of superfluid neutrons so small that the crust alone is insufficient to account for the magnitude of neutron star glitches. Inspired by earlier work on ultracold atomic gases in an optical lattice, we investigate fermions with attractive interactions in a periodic lattice in the mean-field approximation. The effects of band structure are suppressed when the pairing gap is of order or greater than the strength of the lattice potential. By applying the results to the inner crust of neutron stars, we conclude that the reduction of the neutron superfluid density is considerably less than previously estimated and, consequently, it is premature to rule out models of glitches based on neutron superfluidity in the crust.

The interaction between nuclei in the inner crust of neutron stars consists of two contributions, the so-called "direct" interaction and an "induced" one due to density changes in the neutron fluid. For large nuclear separations $r$ the contributions from nuclear forces to each of these terms are shown to be nonzero. In the static limit they are equal in magnitude but have opposite signs and they cancel exactly. We analyze earlier results on effective interactions in the light of this finding. We consider the properties of long-wavelength collective modes and, in particular, calculate the degree of mixing between the lattice phonons and the phonons in the neutron superfluid. Using microscopic theory, we calculate the net non-Coulombic contribution to the nucleus--nucleus interaction and show that, for large $r$, the leading term is due to exchange of two phonons and varies as $1/r^7$: it is an analog of the Casimir--Polder interaction between neutral atoms.

In the two decades since the appearance of the book "Bose-Einstein Condensation" in 1995, there have been a number of developments in our understanding of dense matter. After a brief overview of neutron star structure and the Bose-Einstein condensed phases that have been proposed, we describe selected topics, including neutron and proton pairing gaps, the physics of the inner crust of neutron stars, where a neutron fluid penetrates a lattice of nuclei, meson condensates, and pairing in dense quark matter. Especial emphasis is placed on basic physical effects and on connections to the physics of cold atomic gases.

We survey the current status of understanding of pairing and superfluidity of neutrons and protons in neutron stars from a theoretical perspective, with emphasis on basic physical properties. During the past two decades, the blossoming of the field of ultracold atomic gases and the development of quantum Monte Carlo methods for solving the many-body problem have been two important sources of inspiration, and we shall describe how these have given insight into neutron pairing gaps. The equilibrium properties and collective oscillations of the inner crust of neutron stars, where neutrons paired in a $^1$S$_0$ state coexist with a lattice of neutron-rich nuclei, are also described. While pairing gaps are well understood at densities less than one tenth of the nuclear saturation density, significant uncertainties exist at higher densities due to the complicated nature of nucleon-nucleon interactions, the difficulty of solving the many-body problem under these conditions, and the increasing importance of many-nucleon interactions. We also touch more briefly on the subject of pairing of neutrons in other angular momentum states, specifically the $^3$P$_2$ state, as well as pairing of protons.

Elastic properties of the solid regions of neutron star crusts and white dwarfs play an important role in theories of stellar oscillations. Matter in compact stars is presumably polycrystalline and, since the elastic properties of single crystals of such matter are very anisotropic, it is necessary to relate elastic properties of the polycrystal to those of a single crystal. We calculate the effective shear modulus of polycrystalline matter with randomly oriented crystallites using a self-consistent theory that has been very successful in applications to terrestrial materials and show that previous calculations overestimate the shear modulus by approximately 28%.

At low energies nucleon-nucleon interactions are resonant and therefore supernova matter at subnuclear densities has many similarities to atomic gases with interactions dominated by a Feshbach resonance. We calculate the rates of neutrino processes involving nucleon-nucleon collisions and show that these are enhanced in mixtures of neutrons and protons at subnuclear densities due to the large scattering lengths. As a result, the rate for neutrino pair bremsstrahlung and absorption is significantly larger below 10^{13} g cm^{-3} compared to rates used in supernova simulations.

In the standard picture of the crust of a neutron star, matter there is simple: a body-centered-cubic (bcc) lattice of nuclei immersed in an essentially uniform electron gas. We show that at densities above that for neutron drip ($\sim4\times10^11$) g cm$^{-3}$ or roughly one thousandth of nuclear matter density, the interstitial neutrons give rise to an attractive interaction between nuclei that renders the lattice unstable. We argue that the likely equilibrium structure is similar to that in displacive ferroelectric materials such as BaTiO$_3$. As a consequence, properties of matter in the inner crust are expected to be much richer than previously appreciated and we mention consequences for observable neutron star properties.

Microscopic calculations of neutron matter based on nuclear interactions derived from chiral effective field theory, combined with the recent observation of a 1.97 +- 0.04 M_sun neutron star, constrain the equation of state of neutron-rich matter at sub- and supranuclear densities. We discuss in detail the allowed equations of state and the impact of our results on the structure of neutron stars, the crust-core transition density, and the nuclear symmetry energy. In particular, we show that the predicted range for neutron star radii is robust. For use in astrophysical simulations, we provide detailed numerical tables for a representative set of equations of state consistent with these constraints.

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 investigate neutrino processes for conditions reached in simulations of core-collapse supernovae. Where neutrino-matter interactions play an important role, matter is partially degenerate, and we extend earlier work that addressed the degenerate regime. We derive expressions for the spin structure factor in neutron matter, which is a key quantity required for evaluating rates of neutrino processes. We show that, for essentially all conditions encountered in the post-bounce phase of core-collapse supernovae, it is a very good approximation to calculate the spin relaxation rates in the nondegenerate limit. We calculate spin relaxation rates based on chiral effective field theory interactions and find that they are typically a factor of two smaller than those obtained using the standard one-pion-exchange interaction alone.

We show that microscopic calculations based on chiral effective field theory interactions constrain the properties of neutron-rich matter below nuclear densities to a much higher degree than is reflected in commonly used equations of state. Combined with observed neutron star masses, our results lead to a radius R = 9.7 - 13.9 km for a 1.4 M_{solar} star, where the theoretical range is due, in about equal amounts, to uncertainties in many-body forces and to the extrapolation to high densities.

A simple description of superfluid hydrodynamics in the inner crust of a neutron star is given. Particular attention is paid to the effect of the lattice of nuclei on the properties of the superfluid neutrons, and the effects of entrainment, the fact that some fraction of the neutrons are locked to the motion of the protons in nuclei.

We calculate neutrino processes involving two nucleons at subnuclear densities using chiral effective field theory. Shorter-range noncentral forces reduce the neutrino rates significantly compared with the one-pion exchange approximation currently used in supernova simulations. For densities rho < 10^{14} g cm^{-3}, we find that neutrino rates are well constrained by nuclear interactions and nucleon-nucleon scattering data. As an application, we calculate the mean-square energy transfer in scattering of a neutrino from nucleons and find that collision processes and spin-dependent mean-field effects dominate over the energy transfer due to nucleon recoil.

We consider the spin response of a normal Fermi liquid with noncentral interactions under conditions intermediate between the collisionless and hydrodynamic regimes. This problem is of importance for calculations of neutrino properties in dense matter. By expressing the deviation of the quasiparticle distribution function from equilibrium in terms of eigenfunctions of the transport equation under the combined influence of collisions and an external field, we derive a closed expression for the spin-density--spin-density response function and compare its predictions with that of a relaxation time approximation. Our results indicate that the relaxation time approximation is reliable for neutrino properties under astrophysically relevant conditions.

We present a unified approach to neutrino processes in nucleon matter based on Landau's theory of Fermi liquids that includes one- and two-quasiparticle-quasihole pair states as well as mean-field effects. We show how rates of neutrino processes involving two nucleons may be calculated in terms of the collision integral in the Landau transport equation for quasiparticles. Using a relaxation time approximation, we solve the transport equation for density and spin-density fluctuations and derive a general form for the response functions. We apply our approach to neutral-current processes in neutron matter, where the spin response function is crucial for calculations of neutrino elastic and inelastic scattering, neutrino-pair bremsstrahlung and absorption from strongly-interacting nucleons. We calculate the relaxation rates using modern nuclear interactions and including many-body contributions, and find that rates of neutrino processes are reduced compared with estimates based on the one-pion exchange interaction, which is used in current simulations of core-collapse supernovae.

We analyse the long-wavelength response of a normal Fermi liquid using Landau theory. We consider contributions from intermediate states containing one additional quasiparticle-quasihole pair as well as those from states containing two or more additional quasiparticle-quasihole pairs. For the response of an operator corresponding to a conserved quantity, we show that the behavior of matrix elements to states with more than one additional quasiparticle-quasihole pair at low excitation energies $\omega$ varies as $1/\omega$. It is shown how rates of processes involving transitions to two quasiparticle-quasihole states may be calculated in terms of the collision integral in the Landau transport equation for quasiparticles.

We use Landau's theory of a normal Fermi liquid to derive expressions for the static response of a system with a general tensor interaction that conserves the total spin and the total angular momentum of the quasiparticle-quasihole pair. The magnetic susceptibility is calculated in detail, with the inclusion of the center of mass tensor and cross vector terms in addition to the exchange tensor one. We also introduce a new parametrization of the tensor Landau parameters which significantly reduces the importance of high angular harmonic contributions. For nuclear matter and neutron matter we find that the two most important effects of the tensor interaction are to give a contribution from multipair states and to renormalize the magnetic moments. Response to a weak probe may be calculated using similar methods, replacing the magnetic moments with the matrix elements of the weak charges.

Observation of cooling neutron stars can potentially provide information about the states of matter at supernuclear densities. We review physical properties important for cooling such as neutrino emission processes and superfluidity in the stellar interior, surface envelopes of light elements due to accretion of matter and strong surface magnetic fields. The neutrino processes include the modified Urca process, and the direct Urca process for nucleons and exotic states of matter such as a pion condensate, kaon condensate, or quark matter. The dependence of theoretical cooling curves on physical input and observations of thermal radiation from isolated neutron stars are described. The comparison of observation and theory leads to a unified interpretation in terms of three characteristic types of neutron stars: high-mass stars which cool primarily by some version of the direct Urca process; low-mass stars, which cool via slower processes; and medium-mass stars, which have an intermediate behavior. The related problem of thermal states of transiently accreting neutron stars with deep crustal burning of accreted matter is discussed in connection with observations of soft X-ray transients.

We analyse the effect of non-central forces on the magnetic susceptibility of degenerate Fermi systems. These include the presence of contributions from transitions to states containing more than one quasiparticle-quasihole pair, which cannot be calculated within the framework of Landau Fermi-liquid theory, and renormalization of the quasiparticle magnetic moment, as well as explicit non-central contributions to the quasiparticle interaction. Consequently, the relationship between the Landau parameters and the magnetic susceptibility for Fermi systems with non-central forces is considerably more complicated than for systems with central forces. We use sum-rule arguments to place a lower bound on the contribution to the static susceptibility coming from transitions to multipair states.

The study of tensor correlations and pions in dense matter in astronomical contexts has recently acquired renewed importance for a number of reasons. One is the discovery by Akmal, and Pandharipande (Phys. Rev. C56 2261 (1997)) of a striking increase in tensor correlations in cold matter at densities not far above that of isospin-symmetric nuclear matter at saturation. A second is that, according to a calculation by Mayle and Wilson, softening of the equation of state by the appearance of thermal pions can play an important role. A third is that the strength of tensor correlations is a crucial ingredient in calculating rates of neutrino emission and scattering, which are needed in calculations of stellar collapse and of the cooling of neutron stars. We describe some recent results for the response of nuclear matter when tensor interactions are taken into account, and indicate some unresolved issues.

A review of properties of matter in the interior of neutron stars is given. Particular attention is paid to recent many-body theory calculations of the properties of dense matter. Among topics discussed are the strong increase of tensor correlations at relatively low densities, the "relativistic boost term" in the interaction, and the sensitivity of properties of neutron star models to three-body forces.

Neutrino-pair bremsstrahlung by relativistic degenerate electrons in a neutron-star crust at densities (10^9 - 1.5x10^{14}) g/cm^3 is analyzed. The processes taken into account are neutrino emission due to Coulomb scattering of electrons by atomic nuclei in a Coulomb liquid, and electron-phonon scattering and Bragg diffraction (the static-lattice contribution) in a Coulomb crystal. The static-lattice contribution is calculated including the electron band-structure effects for cubic Coulomb crystals of different types and also for the liquid crystal phases composed of rod- and plate-like nuclei in the neutron-star mantle (at 10^{14} - 1.5x10^{14} g/cm^3). The phonon contribution is evaluated with proper treatment of the multi-phonon processes which removes a jump in the neutrino bremsstrahlung emissivity at the melting point obtained in previous works. Below 10^{13} g/cm^3, the results are rather insensitive to the nuclear form factor, but results for the solid state near the melting point are affected significantly by the Debye-Waller factor and multi-phonon processes. At higher densities, the nuclear form factor becomes more significant. A comparison of the various neutrino generation mechanisms in neutron star crusts shows that electron bremsstrahlung is among the most important ones.

Recent calculations indicate that in the outer parts of neutron stars nuclei are rod-like or slab-like, rather than roughly spherical. We consider the elastic properties of these phases, and argue that they behave as liquid crystals, rather than rigid solids. We estimate elastic constants and discuss implications of our results for neutron star behavior.

We calculate the rate of energy emission by bremsstrahlung of neutrino pairs by electrons moving in the crystalline lattice of ions in dense matter in the crust of a neutron star. Since the periodic potential in the solid gives rise to electronic band gaps which can be as large as about 1 MeV, it is necessary in estimating the bremsstrahlung rate at low temperatures to take into account band structure in detail. We find that, in the densest parts of the inner crust of a neutron star, neutrino emission at temperatures of about $2\times 10^9$K or less is much suppressed compared with earlier estimates that treated the electron-lattice interaction perturbatively, and conclude that neutrino pair bremsstrahlung by electrons in the crusts of neutron stars is much less important for neutron star thermal evolution than was previously thought.

A brief review is given of neutrino emission processes in dense matter, with particular emphasis on recent developments. These include direct Urca processes for nucleons and hyperons, which can give rise to rapid energy loss from the stellar core without exotic matter, and the effect of band structure on neutrino bremsstrahlung from electrons in the crust, which results in much lower energy losses by this process than had previously been estimated.

We consider the rate at which energy is emitted by neutrinos from the dense interior of neutron stars containing a Bose condensate of pions or kaons. The rates obtained are larger, by a factor of 2, than those found earlier, and are consistent with those found for the direct Urca processes.

We demonstrate that band-structure effects suppress bremsstrahlung of neutrino pairs by electrons in the crusts of neutron stars at temperatures of the order of $5\times 10^9 \,{\rm K}$ and below. Taking this into account, together with the fact that recent work indicates that the masses of neutron star crusts are considerably smaller than previously estimated, we find neutrino pair bremsstrahlung to be much less important for the thermal evolution of neutron stars than earlier calculations suggested.