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.

The NICER collaboration recently published a joint estimate of the mass and the radius of PSR J0030+0451, derived via X-ray pulse-profile modeling. Raaijmakers et al. (2019) explored the implications of this measurement for the dense matter equation of state (EOS) using two parameterizations of the high-density EOS: a piecewise-polytropic model, and a model based on the speed of sound in neutron stars. In this work we obtain further constraints on the EOS following this approach, but we also include information about the tidal deformability of neutron stars from the gravitational wave signal of the compact binary merger GW170817. We compare the constraints on the EOS to those set by the recent measurement of a 2.14 solar mass pulsar, included as a likelihood function approximated by a Gaussian, and find a small increase in information gain. To show the flexibility of our method, we also explore the possibility that GW170817 was a neutron star-black hole merger, which yields weaker constraints on the EOS.

We investigate how a phase transition from neutron-star matter to spin-polarized neutron matter affects the equation of state and mass-radius relation of neutron stars. While general extension schemes for the equation of state allow for high pressures inside neutron stars, we find that a phase transition to spin-polarized neutron matter excludes extreme regimes. Hence, such a transition limits the maximum mass of neutron stars to lie below 2.6$-$$2.9 \, M_{\odot}$, depending on the microscopic nuclear forces used, while significantly larger masses could be reached without these constraints. These limits are in good agreement with recent constraints extracted from the neutron-star merger GW170817 and its electromagnetic counterpart. Assuming the description in terms of spin-polarized neutron matter to be valid in the center of neutron stars, we find that stars with a large spin-polarized domain in their core are ruled out by GW170817.

We investigate the impact of different properties of the nuclear equation of state in core-collapse supernovae, with a focus on the proto-neutron-star contraction and its impact on the shock evolution. To this end, we introduce a range of equations of state that vary the nucleon effective mass, incompressibility, symmetry energy, and nuclear saturation point. This allows us to point to the different effects in changing these properties from the Lattimer and Swesty to the Shen et al. equations of state, the two most commonly used equations of state in simulations. In particular, we trace the contraction behavior to the effective mass, which determines the thermal nucleonic contributions to the equation of state. Larger effective masses lead to lower pressures at nuclear densities and a lower thermal index. This results in a more rapid contraction of the proto-neutron star and consequently higher neutrino energies, which aids the shock evolution to a faster explosion.

Both the mass and radius of the millisecond pulsar PSR J0030+0451 have been inferred via pulse-profile modeling of X-ray data obtained by NASA's NICER mission. In this Letter we study the implications of the mass-radius inference reported for this source by Riley et al. (2019) for the dense matter equation of state (EOS), in the context of prior information from nuclear physics at low densities. Using a Bayesian framework we infer central densities and EOS properties for two choices of high-density extensions: a piecewise-polytropic model and a model based on assumptions of the speed of sound in dense matter. Around nuclear saturation density these extensions are matched to an EOS uncertainty band obtained from calculations based on chiral effective field theory interactions, which provide a realistic description of atomic nuclei as well as empirical nuclear matter properties within uncertainties. We further constrain EOS expectations with input from the current highest measured pulsar mass; together, these constraints offer a narrow Bayesian prior informed by theory as well as laboratory and astrophysical measurements. The NICER mass-radius likelihood function derived by Riley et al. (2019) using pulse-profile modeling is consistent with the highest-density region of this prior. The present relatively large uncertainties on mass and radius for PSR J0030+0451 offer, however, only a weak posterior information gain over the prior. We explore the sensitivity to the inferred geometry of the heated regions that give rise to the pulsed emission, and find a small increase in posterior gain for an alternative (but less preferred) model. Lastly, we investigate the hypothetical scenario of increasing the NICER exposure time for PSR J0030+0451.

We exploit the many-body self-consistent Green's function method to analyze finite-temperature properties of infinite nuclear matter and to explore the behavior of the thermal index used to simulate thermal effects in equations of state for astrophysical applications. We show how the thermal index is both density and temperature dependent, unlike often considered, and we provide an error estimate based on our ${\it ab~initio}$ calculations. The inclusion of many-body forces is found to be critical for the density dependence of the thermal index. We also compare our results to a parametrization in terms of the density dependence of the nucleon effective mass. Our study questions the validity of predictions made for the gravitational-wave signal from neutron-star merger simulations with a constant thermal index.

Understanding the dense matter equation of state at extreme conditions is an important open problem. Astrophysical observations of neutron stars promise to solve this, with NICER poised to make precision measurements of mass and radius for several stars using the waveform modelling technique. What has been less clear, however, is how these mass-radius measurements might translate into equation of state constraints and what are the associated equation of state sensitivities. We use Bayesian inference to explore and contrast the constraints that would result from different choices for the equation of state parametrization; comparing the well-established piecewise polytropic parametrization to one based on physically motivated assumptions for the speed of sound in dense matter. We also compare the constraints resulting from Bayesian inference to those from simple compatibility cuts. We find that the choice of equation of state parametrization and particularly its prior assumptions can have a significant effect on the inferred global mass-radius relation and the equation of state constraints. Our results point to important sensitivities when inferring neutron star and dense matter properties. This applies also to inferences from gravitational wave observations.

We present first results on the scalar WIMP-pion coupling from 1 t$\times$yr of exposure with the XENON1T experiment. This interaction is generated when the WIMP couples to a virtual pion exchanged between the nucleons in a nucleus. In contrast to most non-relativistic operators, these pion-exchange currents can be coherently enhanced by the total number of nucleons, and therefore may dominate in scenarios where spin-independent WIMP-nucleon interactions are suppressed. Moreover, for natural values of the couplings, they dominate over the spin-dependent channel due to their coherence in the nucleus. Using the signal model of this new WIMP-pion channel, no significant excess is found, leading to an upper limit cross section of $6.4\times10^{-46}$ cm$^2$ (90 % confidence level) at 30 GeV/c$^2$ WIMP mass.

We present an efficient Monte Carlo framework for perturbative calculations of infinite nuclear matter based on chiral two-, three-, and four-nucleon interactions. The method enables the incorporation of all many-body contributions in a straightforward and transparent way, and makes it possible to extract systematic uncertainty estimates by performing order-by-order calculations in the chiral expansion as well as the many-body expansion. The versatility of this new framework is demonstrated by applying it to chiral low-momentum interactions, exhibiting a very good many-body convergence up to fourth order. Following these benchmarks, we explore new chiral interactions up to next-to-next-to-next-to-leading order (N$^3$LO). Remarkably, simultaneous fits to the triton and to saturation properties can be achieved, while all three-nucleon low-energy couplings remain natural. The theoretical uncertainties of nuclear matter are significantly reduced when going from next-to-next-to-leading order to N$^3$LO.

In this White Paper we present the potential of the Enhanced X-ray Timing and Polarimetry (eXTP) mission for determining the nature of dense matter; neutron star cores host an extreme density regime which cannot be replicated in a terrestrial laboratory. The tightest statistical constraints on the dense matter equation of state will come from pulse profile modelling of accretion-powered pulsars, burst oscillation sources, and rotation-powered pulsars. Additional constraints will derive from spin measurements, burst spectra, and properties of the accretion flows in the vicinity of the neutron star. Under development by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

The standard interpretation of direct-detection limits on dark matter involves particular assumptions of the underlying WIMP-nucleus interaction, such as, in the simplest case, the choice of a Helm form factor that phenomenologically describes an isoscalar spin-independent interaction. In general, the interaction of dark matter with the target nuclei may well proceed via different mechanisms, which would lead to a different shape of the corresponding nuclear structure factors as a function of the momentum transfer $q$. We study to what extent different WIMP-nucleus responses can be differentiated based on the $q$-dependence of their structure factors (or "form factors"). We assume an overall strength of the interaction consistent with present spin-independent limits and consider an exposure corresponding to XENON1T-like, XENONnT-like, and DARWIN-like direct detection experiments. We find that, as long as the interaction strength does not lie too much below current limits, the DARWIN settings allow a conclusive discrimination of many different response functions based on their $q$-dependence, with immediate consequences for elucidating the nature of dark matter.

Searches for invisible Higgs decays at the Large Hadron Collider constrain dark matter Higgs-portal models, where dark matter interacts with the Standard Model fields via the Higgs boson. While these searches complement dark matter direct-detection experiments, a comparison of the two limits depends on the coupling of the Higgs boson to the nucleons forming the direct-detection nuclear target, typically parameterized in a single quantity $f_N$. We evaluate $f_N$ using recent phenomenological and lattice-QCD calculations, and include for the first time the coupling of the Higgs boson to two nucleons via pion-exchange currents. We observe a partial cancellation for Higgs-portal models that makes the two-nucleon contribution anomalously small. Our results, summarized as $f_N=0.308(18)$, show that the uncertainty of the Higgs-nucleon coupling has been vastly overestimated in the past. The improved limits highlight that state-of-the-art nuclear physics input is key to fully exploiting experimental searches.

We calculate the properties of isospin-asymmetric nuclear matter based on chiral nucleon-nucleon (NN) and three-nucleon (3N) interactions. To this end, we develop an improved normal-ordering framework that allows to include general 3N interactions starting from a plane-wave partial-wave-decomposed form. We present results for the energy per particle for general isospin asymmetries based on a set of different Hamiltonians, study their saturation properties, the incompressibility, symmetry energy, and also provide an analytic parametrization for the energy per particle as a function of density and isospin asymmetry.

Neutron matter is an ideal laboratory for nuclear interactions derived from chiral effective field theory since all contributions are predicted up to next-to-next-to-next-to-leading order (N$^3$LO) in the chiral expansion. By making use of recent advances in the partial-wave decomposition of three- nucleon (3N) forces, we include for the first time N$^3$LO 3N interactions in many-body perturbation theory (MBPT) up to third order and in self-consistent Green's function theory (SCGF). Using these two complementary many-body frameworks we provide improved predictions for the equation of state of neutron matter at zero temperature and also analyze systematically the many-body convergence for different chiral EFT interactions. Furthermore, we present an extension of the normal-ordering framework to finite temperatures. These developments open the way to improved calculations of neutron-rich matter including estimates of theoretical uncertainties for astrophysical applications.

We generalize our virial approach to study the neutral current neutrino response of nuclear matter at low densities. In the long-wavelength limit, the virial expansion makes model-independent predictions for neutrino-nucleon scattering rates and the density S_V and spin S_A responses. We find S_A is significantly reduced from one even at low densities. We provide a simple fit S_A^f(n,T,Y_p) of the axial response as a function of density n, temperature T and proton fraction Y_p. This fit reproduces our model independent virial results at low densities and reproduces the Burrows and Sawyer random phase approximation (RPA) results at high densities. Our fit can be incorporated into supernova simulations in a straight forward manner. Preliminary one dimensional supernova simulations suggest that the reduction in the axial response may enhance neutrino heating rates in the gain region during the accretion phase of a core-collapse supernovae.

We present solutions of the BCS gap equation in the channels ${}^1S_0$ and ${}^3P_2-{}^3F_2$ in neutron matter based on nuclear interactions derived within chiral effective field theory (EFT). Our studies are based on a representative set of nonlocal nucleon-nucleon (NN) plus three-nucleon (3N) interactions up to next-to-next-to-next-to-leading order (N$^3$LO) as well as local and semilocal chiral NN interactions up to N$^2$LO and N$^4$LO, respectively. In particular, we investigate for the first time the impact of subleading 3N forces at N$^3$LO on pairing gaps and also derive uncertainty estimates by taking into account results for pairing gaps at different orders in the chiral expansion. Finally, we discuss different methods for obtaining self-consistent solutions of the gap equation. Besides the widely-used quasi-linear method by Khodel et al. we demonstrate that the modified Broyden method is well applicable and exhibits a robust convergence behavior. In contrast to Khodel's method it is based on a direct iteration of the gap equation without imposing an auxiliary potential and is straightforward to implement.

Neutrino-pair production and annihilation through nucleon-nucleon bremsstrahlung is included in current supernova simulations by rates that are based on the one-pion-exchange approximation. Here we explore the consequences of bremsstrahlung rates based on a modern nuclear interactions for proto-neutron star cooling and the corresponding neutrino emission. We find that despite a reduction of the bremsstrahlung emission by a factor of 2-5 in the neutrinospheric region, models with the improved treatment exhibit only $\lesssim$5% changes of the neutrino luminosities and an increase of $\lesssim$0.7 MeV of the average energies of the radiated neutrino spectra, with the largest effects for the antineutrinos of all flavors and at late times. Overall, the proto-neutron star cooling evolution is slowed down modestly by $\lesssim$0.5-1 s.

We propose a formalism for the analysis of direct-detection dark-matter searches that covers all coherent responses for scalar and vector interactions and incorporates QCD constraints imposed by chiral symmetry, including all one- and two-body WIMP-nucleon interactions up to third order in chiral effective field theory. One of the free parameters in the WIMP-nucleus cross section corresponds to standard spin-independent searches, but in general different combinations of new-physics couplings are probed. We identify the interference with the isovector counterpart of the standard spin-independent response and two-body currents as the dominant corrections to the leading spin-independent structure factor, and discuss the general consequences for the interpretation of direct-detection experiments, including minimal extensions of the standard spin-independent analysis. Fits for all structure factors required for the scattering off xenon targets are provided based on state-of-the-art nuclear shell-model calculations.

One of the primary science goals of the next generation of hard X-ray timing instruments is to determine the equation of state of the matter at supranuclear densities inside neutron stars, by measuring the radius of neutron stars with different masses to accuracies of a few percent. Three main techniques can be used to achieve this goal. The first involves waveform modelling. The flux we observe from a hotspot on the neutron star surface offset from the rotational pole will be modulated by the star's rotation, giving rise to a pulsation. Information about mass and radius is encoded into the pulse profile via relativistic effects, and tight constraints on mass and radius can be obtained. The second technique involves characterising the spin distribution of accreting neutron stars. The most rapidly rotating stars provide a very clean constraint, since the mass-shedding limit is a function of mass and radius. However the overall spin distribution also provides a guide to the torque mechanisms in operation and the moment of inertia, both of which can depend sensitively on dense matter physics. The third technique is to search for quasi-periodic oscillations in X-ray flux associated with global seismic vibrations of magnetars (the most highly magnetized neutron stars), triggered by magnetic explosions. The vibrational frequencies depend on stellar parameters including the dense matter equation of state. We illustrate how these complementary X-ray timing techniques can be used to constrain the dense matter equation of state, and discuss the results that might be expected from a 10m$^2$ instrument. We also discuss how the results from such a facility would compare to other astronomical investigations of neutron star properties. [Modified for arXiv]

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 present a common chiral power-counting scheme for vector, axial-vector, scalar, and pseudoscalar WIMP-nucleon interactions, and derive all one- and two-body currents up to third order in the chiral expansion. Matching our amplitudes to non-relativistic effective field theory, we find that chiral symmetry predicts a hierarchy amongst the non-relativistic operators. Moreover, we identify interaction channels where two-body currents that so far have not been accounted for become relevant.

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.

We calculate neutrino absorption rates due to charged-current reactions $\nu_e+n \rightarrow e^- + p $ and $\bar{\nu}_e+p \rightarrow e^+ + n $ in the outer regions of a newly born neutron star called the neutrino-sphere. To improve on recent work which has shown that nuclear mean fields enhance the $\nu_e$ cross-section and suppress the $\bar{\nu}_e$ cross-section, we employ realistic nucleon-nucleon interactions that fit measured scattering phase shifts. Using these interactions we calculate the momentum-, density-, and temperature-dependent nucleon self-energies in the Hartree-Fock approximation. A potential derived from chiral effective field theory and a pseudo-potential constructed to reproduce nucleon-nucleon phase shifts at the mean-field level are used to study the equilibrium proton fraction and the charged-current rates are studied in detail. We compare our results to earlier calculations obtained using phenomenological mean-field models and to those obtained in the virial expansion valid at low density. We find that for typical ambient conditions in the neutrino-sphere, $T=5-10$ MeV and $\rho =10^{11}-10^{13}$ g/cm$^3$, the difference between the $\nu_{e}$ and $\bar{\nu}_e$ absorption rates is much larger than predicted earlier. Our results have important implications for heavy-element nucleosynthesis in supernovae and for supernova neutrino detection.

We study the properties of spin-polarized neutron matter at next-to-next-to-next-to-leading order in chiral effective field theory, including two-, three-, and four-neutron interactions. The energy of spin-polarized neutrons is remarkably close to a non-interacting system at least up to saturation density, where interaction effects provide less than 10% corrections. This shows that the physics of neutron matter is similar to a unitary gas well beyond the scattering-length regime. Implications for energy-density functionals and for a possible ferromagnetic transition in neutron stars are discussed. Our predictions can be tested with lattice QCD, and we present results for varying pion mass.

We study the structure factors for spin-independent WIMP scattering off xenon based on state-of-the-art large-scale shell-model calculations, which are shown to yield a good spectroscopic description of all experimentally relevant isotopes. Our results are based on the leading scalar one-body currents only. At this level and for the momentum transfers relevant to direct dark matter detection, the structure factors are in very good agreement with the phenomenological Helm form factors used to give experimental limits for WIMP-nucleon cross sections. In contrast to spin-dependent WIMP scattering, the spin-independent channel, at the one-body level, is less sensitive to nuclear structure details. In addition, we explicitly show that the structure factors for inelastic scattering are suppressed by ~ 10^{-4} compared to the coherent elastic scattering response. This implies that the detection of inelastic scattering will be able to discriminate clearly between spin-independent and spin-dependent scattering. Finally, we provide fits for all calculated structure factors.

We perform state-of-the-art large-scale shell-model calculations of the structure factors for elastic spin-dependent WIMP scattering off 129,131Xe, 127I, 73Ge, 19F, 23Na, 27Al, and 29Si. This comprehensive survey covers the non-zero-spin nuclei relevant to direct dark matter detection. We include a pedagogical presentation of the formalism necessary to describe elastic and inelastic WIMP-nucleus scattering. The valence spaces and nuclear interactions employed have been previously used in nuclear structure calculations for these mass regions and yield a good spectroscopic description of these isotopes. We use spin-dependent WIMP-nucleus currents based on chiral effective field theory (EFT) at the one-body level and including the leading long-range two-body currents due to pion exchange, which are predicted in chiral EFT. Results for all structure factors are provided with theoretical error bands due to the nuclear uncertainties of WIMP currents in nuclei.

We study the impact of the nucleon-$\alpha$ P-wave resonances on neutrino-pair bremsstrahlung. Because of the non-central spin-orbit interaction, these resonances lead to an enhanced contribution to the nucleon spin structure factor for temperatures $T \lesssim 4$ MeV. If the $\alpha$-particle fraction is significant and the temperature is in this range, this contribution is competitive with neutron-neutron bremsstrahlung. This may be relevant for neutrino production in core-collapse supernovae or other dense astrophysical environments. Similar enhancements are expected for resonant non-central nucleon-nucleus interactions.

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.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

We investigate asymmetric nuclear matter with two- and three-nucleon interactions based on chiral effective field theory, where three-body forces are fit only to light nuclei. Focusing on neutron-rich matter, we calculate the energy for different proton fractions and include estimates of the theoretical uncertainty. We use our ab-initio results to test the quadratic expansion around symmetric matter with the symmetry energy term, and confirm its validity for highly asymmetric systems. Our calculations are in remarkable agreement with an empirical parametrization for the energy density. These findings are very useful for astrophysical applications and for developing new equations of state.

Direct dark matter detection focuses on elastic scattering of dark matter particles off nuclei. In this study, we explore inelastic scattering where the nucleus is excited to a low-lying state of 10-100 keV, with subsequent prompt de-excitation. We calculate the inelastic structure factors for the odd-mass xenon isotopes based on state-of-the-art large-scale shell-model calculations with chiral effective field theory WIMP-nucleon currents. For these cases, we find that the inelastic channel is comparable to or can dominate the elastic channel for momentum transfers around 150 MeV. We calculate the inelastic recoil spectra in the standard halo model, compare these to the elastic case, and discuss the expected signatures in a xenon detector, along with implications for existing and future experiments. The combined information from elastic and inelastic scattering will allow to determine the dominant interaction channel within one experiment. In addition, the two channels probe different regions of the dark matter velocity distribution and can provide insight into the dark halo structure. The allowed recoil energy domain and the recoil energy at which the integrated inelastic rates start to dominate the elastic channel depend on the mass of the dark matter particle, thus providing a potential handle to constrain its mass.

Neutron matter presents a unique system in chiral effective field theory (EFT), because all many-body forces among neutrons are predicted to next-to-next-to-next-to-leading order (N3LO). We discuss perturbative and first Quantum Monte Carlo (QMC) calculations of neutron matter with chiral EFT interactions and their astrophysical impact for the equation of state and neutron stars.

We calculate the chiral condensate in neutron matter at zero temperature based on nuclear forces derived within chiral effective field theory. Two-, three- and four-nucleon interactions are included consistently to next-to-next-to-next-to-leading order (N3LO) of the chiral expansion. We find that the interaction contributions lead to a modest increase of the condensate, thus impeding the restoration of chiral symmetry in dense matter and making a chiral phase transition in neutron-rich matter unlikely for densities that are not significantly higher than nuclear saturation density.

The neutron-matter equation of state constrains the properties of many physical systems over a wide density range and can be studied systematically using chiral effective field theory (EFT). In chiral EFT, all many-body forces among neutrons are predicted to next-to-next-to-next-to-leading order (N3LO). We present details and additional results of the first complete N3LO calculation of the neutron-matter energy, which includes the subleading three-nucleon as well as the leading four-nucleon forces, and provides theoretical uncertainties. In addition, we discuss the impact of our results for astrophysics: for the supernova equation of state, the symmetry energy and its density derivative, and for the structure of neutron stars. Finally, we give a first estimate for the size of the N3LO many-body contributions to the energy of symmetric nuclear matter, which shows that their inclusion will be important in nuclear structure calculations.

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.

Neutron matter presents a unique system for chiral effective field theory (EFT), because all many-body forces among neutrons are predicted to next-to-next-to-next-to-leading order (N3LO). We present the first complete N3LO calculation of the neutron matter energy. This includes the subleading three-nucleon (3N) forces for the first time and all leading four-nucleon (4N) forces. We find relatively large contributions from N3LO 3N forces. Our results provide constraints for neutron-rich matter in astrophysics with controlled theoretical uncertainties.

We present the first study of three-nucleon (3N) forces for proton-rich nuclei along the N=8 and N=20 isotones. Our results for the ground-state energies and proton separation energies are in very good agreement with experiment where available, and with the empirical isobaric multiplet mass equation. We predict the spectra for all N=8 and N=20 isotones to the proton dripline, which agree well with experiment for 18Ne, 19Na, 20Mg and 42Ti. In all other cases, we provide first predictions based on nuclear forces. Our results are also very promising for studying isospin symmetry breaking in medium-mass nuclei based on chiral effective field theory.

Chiral effective field theory (EFT) provides a systematic expansion for the coupling of WIMPs to nucleons at the momentum transfers relevant to direct cold dark matter detection. We derive the currents for spin-dependent WIMP scattering off nuclei at the one-body level and include the leading long-range two-body currents, which are predicted in chiral EFT. As an application, we calculate the structure factor for spin-dependent WIMP scattering off 129,131Xe nuclei, using nuclear interactions that have been developed to study nuclear structure and double-beta decays in this region. We provide theoretical error bands due to the nuclear uncertainties of WIMP currents in nuclei.

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

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.

Neutron-star (NS) merger simulations are conducted for 38 representative microphysical descriptions of high-density matter in order to explore the equation-of-state dependence of the postmerger ring-down phase. The formation of a deformed, oscillating, differentially rotating very massive NS is the typical outcome of the coalescence of two stars with 1.35 $M_{\odot}$ for most candidate EoSs. The oscillations of this object imprint a pronounced peak in the gravitational-wave (GW) spectra, which is used to characterize the emission for a given model. The peak frequency of this postmerger GW signal correlates very well with the radii of nonrotating NSs, and thus allows to constrain the high-density EoS by a GW detection. In the case of 1.35-1.35 $M_{\odot}$ mergers the peak frequency scales particularly well with the radius of a NS with 1.6 $M_{\odot}$, where the maximum deviation from this correlation is only 60 meters for fully microphysical EoSs which are compatible with NS observations. Combined with the uncertainty in the determination of the peak frequency it appears likely that a GW detection can measure the radius of a 1.6 $M_{\odot}$ NS with an accuracy of about 100 to 200 meters. We also uncover relations of the peak frequency with the radii of nonrotating NSs with 1.35 $M_{\odot}$ or 1.8 $M_{\odot}$, with the radius or the central energy density of the maximum-mass Tolman-Oppenheimer-Volkoff configuration, and with the pressure or sound speed at a fiducial rest-mass density of about twice nuclear saturation density. Furthermore, it is found that a determination of the dominant postmerger GW frequency can provide an upper limit for the maximum mass of nonrotating NSs. The prospects for a detection of the postmerger GW signal and a determination of the dominant GW frequency are estimated to be in the range of 0.015 to 1.2 events per year with the upcoming Advanced LIGO detector.

Understanding and predicting the formation of shell structure from nuclear forces is a central challenge for nuclear physics. While the magic numbers N=2,8,20 are generally well understood, N=28 is the first standard magic number that is not reproduced in microscopic theories with two-nucleon forces. In this Letter, we show that three-nucleon forces give rise to repulsive interactions between two valence neutrons that are key to explain 48Ca as a magic nucleus, with a high 2+ excitation energy and a concentrated magnetic dipole transition strength. The repulsive three-nucleon mechanism improves the agreement with experimental binding energies.

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.

We calculate the properties of neutron matter and highlight the physics of chiral three-nucleon forces. For neutrons, only the long-range 2 pi-exchange interactions of the leading chiral three-nucleon forces contribute, and we derive density-dependent two-body interactions by summing the third particle over occupied states in the Fermi sea. Our results for the energy suggest that neutron matter is perturbative at nuclear densities. We study in detail the theoretical uncertainties of the neutron matter energy, provide constraints for the symmetry energy and its density dependence, and explore the impact of chiral three-nucleon forces on the S-wave superfluid pairing gap.

The limit of neutron-rich nuclei, the neutron drip-line, evolves regularly from light to medium-mass nuclei except for a striking anomaly in the oxygen isotopes. This anomaly is not reproduced in shell-model calculations derived from microscopic two-nucleon forces. Here, we present the first microscopic explanation of the oxygen anomaly based on three-nucleon forces that have been established in few-body systems. This leads to repulsive contributions to the interactions among excess neutrons that change the location of the neutron drip-line from $^{28}$O to the experimentally observed $^{24}$O. Since the mechanism is robust and general, our findings impact the prediction of the most neutron-rich nuclei and the synthesis of heavy elements in neutron-rich environments.

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.

Nuclear theory has entered an exciting era. This is due to advances on many fronts, including the development of effective field theory and the renormalization group for nuclear forces, advances in ab-initio methods for nuclear structure, an effort to develop a universal density functional based on microscopic interactions, and the application of large-scale computing resources. I discuss their impact, recent highlights and the frontiers in understanding and predicting nuclei and the structure of strongly-interacting matter based on chiral interactions.

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 study the composition of the outer layers of a protoneutron star and show that light nuclei are present in substantial amounts. The composition is dominated by nucleons, deuterons, tritons and alpha particles; 3He is present in smaller amounts. This composition can be studied in laboratory experiments with new neutron-rich radioactive beams that can reproduce similar densities and temperatures. After including the corresponding neutrino interactions, we demonstrate that light nuclei have a small impact on the average energy of the emitted electron neutrinos, but are significant for the average energy of antineutrinos. During the early post-explosion phase, the average energy of electron antineutrinos is slightly increased, while at later times during the protoneutron star cooling it is reduced by about 1 MeV. The consequences of these changes for nucleosynthesis in neutrino-driven supernova outflows are discussed.

We calculate the neutron matter equation of state at finite temperature based on low-momentum two- and three-nucleon interactions. The free energy is obtained from a loop expansion around the Hartree-Fock energy, including contributions from normal and anomalous diagrams. We focus on densities below saturation density with temperatures T <= 10 MeV and compare our results to the model-independent virial equation of state and to variational calculations. Good agreement with the virial equation of state is found at low density. We provide simple estimates for the theoretical error, important for extrapolations to astrophysical conditions.

We extend the virial equation of state to include 3H and 3He nuclei, and predict significant mass-three fractions near the neutrinosphere in supernovae. While alpha particles are often more abundant, we demonstrate that energy transfer cross-sections for muon and tau neutrinos at low densities are dominated by breakup of the loosely-bound 3H and 3He nuclei. The virial coefficients involving A=3 nuclei are calculated directly from the corresponding nucleon-3H and nucleon-3He scattering phase shifts. For the neutral-current inelastic cross-sections and the energy transfer cross sections, we perform ab-initio calculations based on microscopic two- and three-nucleon interactions and meson-exchange currents.

We study in detail the dependence of the 1S0 superfluid pairing gap on nuclear interactions and on charge-independence breaking at the BCS level. Starting from chiral effective-field theory and conventional nucleon-nucleon (NN) interactions, we use the renormalization group to generate low-momentum interactions V_{low k} with sharp and smooth regulators. The resulting BCS gaps are well constrained by the NN scattering phase shifts, and the cutoff dependence is very weak for sharp or sufficiently narrow smooth regulators with cutoffs Lambda > 1.6 fm^{-1}. It is therefore likely that the effect of three-nucleon interactions on 1S0 superfluidity is small at the BCS level. The charge dependence of nuclear interactions has a 10% effect on the pairing gap.

We discuss superfluidity in neutron matter, with particular attention to induced interactions and to universal properties accessible with cold atoms.

We calculate the neutron matter equation of state at finite temperature based on low-momentum nucleon-nucleon and three-nucleon interactions. Our results are compared to the model-independent virial equation of state and to variational calculations. We provide a simple estimate for the theoretical error, important for extrapolations to astrophysical conditions.

We generalize our virial approach to study spin-polarized neutron matter and the consistent neutrino response at low densities. In the long-wavelength limit, the virial expansion makes model-independent predictions for the density and spin response, based only on nucleon-nucleon scattering data. Our results for the neutrino response provide constraints for random-phase approximation or other model calculations, and we compare the virial vector and axial response to response functions used in supernova simulations. The virial expansion is suitable to describe matter near the supernova neutrinosphere, and this work extends the virial equation of state to predict neutrino interactions in neutron matter.

We present the virial equation of state of low-density nuclear matter composed of neutrons, protons and alpha particles. The virial equation of state is model-independent, and therefore sets a benchmark for all nuclear equations of state at low densities. We calculate the second virial coefficients for nucleon-nucleon, nucleon-alpha and alpha-alpha interactions directly from the relevant binding energies and scattering phase shifts. The virial approach systematically takes into account contributions from bound nuclei and the resonant continuum, and consequently provides a framework to include strong-interaction corrections to nuclear statistical equilibrium models. The virial coefficients are used to make model-independent predictions for a variety of properties of nuclear matter over a range of densities, temperatures and compositions. Our results provide constraints on the physics of the neutrinosphere in supernovae. The resulting alpha particle concentration differs from all equations of state currently used in supernova simulations. Finally, the virial equation of state greatly improves our conceptual understanding of low-density nuclear matter.

We present a model-independent description of low-density neutron matter based on the virial expansion. The virial equation of state provides a benchmark for all nuclear equations of state at densities and temperatures where the interparticle separation is large compared to the thermal wavelength. We calculate the second virial coefficient directly from the nucleon-nucleon scattering phase shifts. Our results for the pressure, energy, entropy and the free energy correctly include the physics of the large neutron-neutron scattering length. We find that, as in the universal regime, thermodynamic properties of neutron matter scale over a wide range of temperatures, but with a significantly reduced interaction coefficient compared to the unitary limit.

We discuss how the renormalization group can be used to derive effective nuclear interactions. Starting from the model-independent low-momentum interaction V_{low k}, we successively integrate out high-lying particle and hole states from momentum shells around the Fermi surface as proposed by Shankar. The renormalization group approach allows for a systematic calculation of induced interactions and yields similar contributions to the scattering amplitude as the two-body parquet equations. We review results for the 1S0 and 3P2 superfluid pairing gaps as well as the spin dependence of effective interactions in neutron matter. Implications for the cooling of neutron stars are discussed.

We revisit the emissivity from neutrino pair bremsstrahlung in neutron-neutron scattering, nn -> nn nu nubar, which was calculated from the one-pion exchange potential including correlation effects by Friman and Maxwell. Starting from the free-space low-momentum nucleon-nucleon interaction V_{low k}, we include tensor, spin-orbit and second-order medium-induced non-central contributions to the scattering amplitude in neutron matter. We find that the screening of the nucleon-nucleon interaction reduces the emissivity from neutrino bremsstrahlung for densities below nuclear matter density. We discuss the implications of medium modifications for the cooling of neutron stars via neutrino emission, taking into account recent results for the polarization effects on neutron superfluidity.

Renormalization group methods can be applied to the nuclear many-body problem using the approach proposed by Shankar. We start with the two-body low momentum interaction V_{low k} and use the RG flow from the particle-hole channels to calculate the full scattering amplitude in the vicinity of the Fermi surface. This is a new straightforward approach to the many-body problem which is applicable also to condensed matter systems without long-range interactions, such as liquid 3He. We derive the one-loop renormalization group equations for the quasiparticle interaction and the scattering amplitude at zero temperature. The RG presents an elegant method to maintain all momentum scales and preserve the antisymmetry of the scattering amplitude. As a first application we solve the RG equations for neutron matter. The resulting quasiparticle interaction includes effects due to the polarization of the medium, the so-called induced interaction of Babu and Brown. We present results for the Fermi liquid parameters, the equation of state of neutron matter and the 1S0 superfluid pairing gap.