The Virial Equation of State of Low-Density Neutron 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.Comment: 7 pages, 6 figures, minor revisions, to appear in Phys. Lett.

The Neutrino Response of Low-Density Neutron Matter from the Virial Expansion

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.Comment: 8 pages, 5 figures, minor additions, to appear in Phys. Lett.

Neutron matter at finite temperature

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.Comment: 15 pages, 6 figure

Links Between Heavy Ion and Astrophysics

Heavy ion experiments provide important data to test astrophysical models. The high density equation of state can be probed in HI collisions and applied to the hot protoneutron star formed in core collapse supernovae. The Parity Radius Experiment (PREX) aims to accurately measure the neutron radius of $^{208}$Pb with parity violating electron scattering. This determines the pressure of neutron rich matter and the density dependence of the symmetry energy. Competition between nuclear attraction and coulomb repulsion can form exotic shapes called nuclear pasta in neutron star crusts and supernovae. This competition can be probed with multifragmentation HI reactions. We use large scale semiclassical simulations to study nonuniform neutron rich matter in supernovae. We find that the coulomb interactions in astrophysical systems suppress density fluctuations. As a result, there is no first order liquid vapor phase transition. Finally, the virial expansion for low density matter shows that the nuclear vapor phase is complex with significant concentrations of alpha particles and other light nuclei in addition to free nucleons.Comment: 8 pages, 6 figures. To be published in "Dynamics and Thermodynamics with Nucleon Degrees of Freedom", eds. P. Chomaz, F. Gulminelli, J. Natowitz, and S. Yennello, http://cyclotron.tamu.edu/wci3/wci_book.htm

A Second Relativistic Mean Field and Virial Equation of State for Astrophysical Simulations

We generate a second equation of state (EOS) of nuclear matter for a wide range of temperatures, densities, and proton fractions for use in supernovae, neutron star mergers, and black hole formation simulations. We employ full relativistic mean field (RMF) calculations for matter at intermediate density and high density, and the Virial expansion of a non-ideal gas for matter at low density. For this EOS we use the RMF effective interaction FSUGold, whereas our earlier EOS was based on the RMF effective interaction NL3. The FSUGold interaction has a lower pressure at high densities compared to the NL3 interaction. We calculate the resulting EOS at over 100,000 grid points in the temperature range $T$ = 0 to 80 MeV, the density range $n_B$ = 10$^{-8}$ to 1.6 fm$^{-3}$, and the proton fraction range $Y_p$ = 0 to 0.56. We then interpolate these data points using a suitable scheme to generate a thermodynamically consistent equation of state table on a finer grid. We discuss differences between this EOS, our NL3 based EOS, and previous EOSs by Lattimer-Swesty and H. Shen et al for the thermodynamic properties, composition, and neutron star structure. The original FSUGold interaction produces an EOS, that we call FSU1.7, that has a maximum neutron star mass of 1.7 solar masses. A modification in the high density EOS is introduced to increase the maximum neutron star mass to 2.1 solar masses and results in a slightly different EOS that we call FSU2.1. The EOS tables for FSU1.7 and FSU2.1 are available for download.Comment: updated version according to referee's comments. Phys. Rev. C in pres

Constraining mean-field models of the nuclear matter equation of state at low densities

An extension of the generalized relativistic mean-field (gRMF) model with density dependent couplings is introduced in order to describe thermodynamical properties and the composition of dense nuclear matter for astrophysical applications. Bound states of light nuclei and two-nucleon scattering correlations are considered as explicit degrees of freedom in the thermodynamical potential. They are represented by quasiparticles with medium-dependent properties. The model describes the correct low-density limit given by the virial equation of state (VEoS) and reproduces RMF results around nuclear saturation density where clusters are dissolved. A comparison between the fugacity expansions of the VEoS and the gRMF model provides consistency relations between the quasiparticles properties, the nucleon-nucleon scattering phase shifts and the meson-nucleon couplings of the gRMF model at zero density. Relativistic effects are found to be important at temperatures that are typical in astrophysical applications. Neutron matter and symmetric matter are studied in detail.Comment: 50 pages, 21 figure

Neutron rich matter, neutron stars, and their crusts

Neutron rich matter is at the heart of many fundamental questions in Nuclear Physics and Astrophysics. What are the high density phases of QCD? Where did the chemical elements come from? What is the structure of many compact and energetic objects in the heavens, and what determines their electromagnetic, neutrino, and gravitational-wave radiations? Moreover, neutron rich matter is being studied with an extraordinary variety of new tools such as Facility for Rare Isotope Beams (FRIB) and the Laser Interferometer Gravitational Wave Observatory (LIGO). We describe the Lead Radius Experiment (PREX) that is using parity violation to measure the neutron radius in 208Pb. This has important implications for neutron stars and their crusts. Using large scale molecular dynamics, we model the formation of solids in both white dwarfs and neutron stars. We find neutron star crust to be the strongest material known, some 10 billion times stronger than steel. It can support mountains on rotating neutron stars large enough to generate detectable gravitational waves. Finally, we describe a new equation of state for supernova and neutron star merger simulations based on the Virial expansion at low densities, and large scale relativistic mean field calculations.Comment: 10 pages, 2 figures, Plenary talk International Nuclear Physics Conference 2010, Vancouver, C

Multi-messenger observations of neutron rich matter

Neutron rich matter is central to many fundamental questions in nuclear physics and astrophysics. Moreover, this material is being studied with an extraordinary variety of new tools such as the Facility for Rare Isotope Beams (FRIB) and the Laser Interferometer Gravitational Wave Observatory (LIGO). We describe the Lead Radius Experiment (PREX) that uses parity violating electron scattering to measure the neutron radius in $^{208}$Pb. This has important implications for neutron stars and their crusts. We discuss X-ray observations of neutron star radii. These also have important implications for neutron rich matter. Gravitational waves (GW) open a new window on neutron rich matter. They come from sources such as neutron star mergers, rotating neutron star mountains, and collective r-mode oscillations. Using large scale molecular dynamics simulations, we find neutron star crust to be very strong. It can support mountains on rotating neutron stars large enough to generate detectable gravitational waves. Finally, neutrinos from core collapse supernovae (SN) provide another, qualitatively different probe of neutron rich matter. Neutrinos escape from the surface of last scattering known as the neutrino-sphere. This is a low density warm gas of neutron rich matter. Observations of neutrinos can probe nucleosyntheses in SN. Simulations of SN depend on the equation of state (EOS) of neutron rich matter. We discuss a new EOS based on virial and relativistic mean field calculations. We believe that combing astronomical observations using photos, GW, and neutrinos, with laboratory experiments on nuclei, heavy ion collisions, and radioactive beams will fundamentally advance our knowledge of compact objects in the heavens, the dense phases of QCD, the origin of the elements, and of neutron rich matter.Comment: 13 pages, 4 figures, Added discussion of dipole polarizability, pygmy resonances, and neutron skin