Search

We add an ensemble of nuclei to the equation of state for homogeneous nucleonic matter to generate a new set of models suitable for astrophysical simulations of core-collapse supernovae and neutron star mergers. We implement empirical constraints from (i) nuclear mass measurements, (ii) proton-proton scattering phase shifts, and (iii) neutron star observations. Our model is also guided by microscopic many-body theory calculations based on realistic nuclear forces, including the zero-temperature neutron matter equation of state from quantum Monte Carlo simulations and thermal contributions to the free energy from finite-temperature many-body perturbation theory. We ensure that the parameters of our model can be varied while preserving thermodynamic consistency and the connection to experimental or observational data, thus providing a probability distribution of the astrophysical hot and dense matter equation of state. We compare our results with those obtained from other available equations of state. While our probability distributions indeed represent a large number of possible equations of state, we cannot yet claim to have fully explored all of the uncertainties, especially with regard to the structure of nuclei in the hot and dense medium., Comment: 13 pages, 9 figures

The gravitational wave emission from the merging binary neutron star system GW170817 arrived full of tidal information which can be used to probe the fundamental ultra-dense nuclear physics residing in these stars. In previous work, we used two-dimensional correlations between nuclear matter parameters and tidal deformabilities of neutron stars applying specifically to GW170817 to derive constraints on the former. Here, we extend this analysis by finding similar correlations for varying chirp masses, the dominant determining factor in the frequency evolution of the inspiral, such that one can apply the same method to future detections. We estimate how accurately one can measure nuclear parameters with future gravitational wave interferometers and show how such measurements can be improved by combining multiple events. We find that bounds on the nuclear parameters with future observations can improve from the current one with GW170817 only by $\sim 30\%$ due to the existence of systematic errors caused mainly by the remaining uncertainty in the equation of state near and just above the nuclear saturation density. We show that such systematic errors can be reduced by considering multidimensional correlations among nuclear parameters and tidal deformabilities with various neutron star masses., Comment: 16 pages, 17 figures; new citations, hybrid star discussion paragraph added; Used conditional probability distributions on tidal deformability and nuclear parameters, rather than the two-dimensional PDFs

Recent observations of neutron-star properties, in particular the recent detection of gravitational waves emitted from binary neutron stars, GW 170817, open the way to put strong constraints on nuclear interactions. In this paper, we review the state of the art in calculating the equation of state of strongly interacting matter from first principle calculations starting from microscopic interactions among nucleons. We then review selected properties of neutron stars that can be directly compared with present and future observations., Comment: 28 pages, 9 figures, Special Issue Article for Journal of Physics G

X-ray observations of transiently accreting neutron stars during quiescence provide information about the structure of neutron star crusts and the properties of dense matter. Interpretation of the observational data requires an understanding of the nuclear reactions that heat and cool the crust during accretion, and define its nonequilibrium composition. We identify here in detail the typical nuclear reaction sequences down to a depth in the inner crust where the mass density is 2E12 g/cm^3 using a full nuclear reaction network for a range of initial compositions. The reaction sequences differ substantially from previous work. We find a robust reduction of crust impurity at the transition to the inner crust regardless of initial composition, though shell effects can delay the formation of a pure crust somewhat to densities beyond 2E12 g/cm^3. This naturally explains the small inner crust impurity inferred from observations of a broad range of systems. The exception are initial compositions with A >= 102 nuclei, where the inner crust remains impure with an impurity parameter of Qimp~20 due to the N = 82 shell closure. In agreement with previous work we find that nuclear heating is relatively robust and independent of initial composition, while cooling via nuclear Urca cycles in the outer crust depends strongly on initial composition. This work forms a basis for future studies of the sensitivity of crust models to nuclear physics and provides profiles of composition for realistic crust models., 25 Pages, accepted for publication in Ap. J

The tidal measurement of gravitational waves from the binary neutron star merger event GW170817 allows us to probe nuclear physics that suffers less from astrophysical systematics compared to neutron star radius measurements with electromagnetic wave observations. A recent work found strong correlation among neutron-star tidal deformabilities and certain combinations of nuclear parameters associated with the equation of state. These relations were then used to derive bounds on such parameters from GW170817 assuming that the relations and neutron star masses are known exactly. Here, we expand on this important work by taking into account a few new considerations: (1) a broader class of equations of state; (2) correlations with the mass-weighted tidal deformability that was directly measured with GW170817; (3) how the relations depend on the binary mass ratio; (4) the uncertainty from equation of state variation in the correlation relations; (5) adopting the updated tidal deformability measurement from GW170817. Upon these new considerations, we find GW170817 bounds on nuclear parameters (the incompressibility $K_0$, its slope $M_0$ and the curvature of symmetry energy $K_{\mathrm{sym},0}$ at nuclear saturation density) to be 81 MeV $\leq K_0 \leq$ 362 MeV, 1556 MeV $\leq M_0 \leq$ 4971 MeV, and -254 MeV $\leq K_\mathrm{sym,0} \leq$ 27 MeV, which are more conservative than previously found with systematic errors more properly taken into account., Comment: 15 pages, 14 figures; Updated the section on computing the EoS scattering uncertainty to include a more comprehensive error analysis; Fixed issue with incorrect EoSs plotted in Fig. 3; Fixed a small mistake on nuclear matter bounds; Used conditional probability distributions rather than 2D ones to compute nuclear parameter bounds

Using a model for the equation of state and composition of dense matter and the magnitude of singlet proton superconductivity and triplet neutron superfluidity, we perform the first simultaneous fit of neutron star masses and radii determined from observations of quiescent low-mass x-ray binaries and luminosities and ages determined from observations of isolated neutron stars. We find that the Vela pulsar strongly determines the values inferred for the superfluid/superconducting gaps and the neutron star radius. We find, regardless of whether or not the Vela pulsar is included in the analysis, that the threshold density for the direct Urca process lies between the central density of 1.7 and 2 solar mass neutron stars. We also find that two solar mass stars are unlikely to cool principally by the direct Urca process because of the suppression by neutron triplet superfluidity., Comment: 8 pages, 4 figures

The cooling phase of thermonuclear (type-I) X-ray bursts can be used to constrain the neutron star (NS) compactness by comparing the observed cooling tracks of bursts to accurate theoretical atmosphere model calculations. By applying the so-called cooling tail method, where the information from the whole cooling track is used, we constrain the mass, radius, and distance for three different NSs in low-mass X-ray binaries 4U 1702-429, 4U 1724-307, and SAX J1810.8-260. Care is taken to only use the hard state bursts where it is thought that only the NS surface alone is emitting. We then utilize a Markov chain Monte Carlo algorithm within a Bayesian framework to obtain a parameterized equation of state (EoS) of cold dense matter from our initial mass and radius constraints. This allows us to set limits on various nuclear parameters and to constrain an empirical pressure-density relation for the dense matter. Our predicted EoS results in NS radius between 10.5-12.8 km (95% confidence limits) for a mass of 1.4 $M_{\odot}$. Due to systematic errors and uncertainty in the composition these results should be interpreted as lower limits for the radius., Comment: 23 pages, 12 figures, revised version, accepted for publication in A&A

We investigate constraints on neutron star structure arising from the assumptions that neutron stars have crusts, that recent calculations of pure neutron matter limit the equation of state of neutron star matter near the nuclear saturation density, that the high-density equation of state is limited by causality and the largest high-accuracy neutron star mass measurement, and that general relativity is the correct theory of gravity. We explore the role of prior assumptions by considering two classes of equation of state models. In a first, the intermediate- and high-density behavior of the equation of state is parameterized by piecewise polytropes. In the second class, the high-density behavior of the equation of state is parameterized by piecewise continuous line segments. The smallest density at which high-density matter appears is varied in order to allow for strong phase transitions above the nuclear saturation density. We critically examine correlations among the pressure of matter, radii, maximum masses, the binding energy, the moment of inertia, and the tidal deformability, paying special attention to the sensitivity of these correlations to prior assumptions about the equation of state. It is possible to constrain the radii of $1.4~\mathrm{M}_{\odot}$ neutron stars to a be larger than 10 km, even without consideration of additional astrophysical observations, for example, those from photospheric radius expansion bursts or quiescent low-mass X-ray binaries. We are able to improve the accuracy of known correlations between the moment of inertia and compactness as well as the binding energy and compactness. We also demonstrate the existence of a correlation between the neutron star binding energy and the moment of inertia., Comment: Eur. Phys. J. A (in press)

Neutron star (binary neutron star and neutron star - black hole) mergers are believed to produce short-duration gamma-ray bursts. They are also believed to be the dominant source of gravitational waves to be detected by the advanced LIGO and the dominant source of the heavy r-process elements in the universe. Whether or not these mergers produce short-duration GRBs depends sensitively on the fate of the core of the remnant (whether, and how quickly, it forms a black hole). In this paper, we combine the results of merger calculations and equation of state studies to determine the fate of the cores of neutron star mergers. Using population studies, we can determine the distribution of these fates to compare to observations. We find that black hole cores form quickly only for equations of state that predict maximum non-rotating neutron star masses below 2.3-2.4 solar masses. If quick black hole formation is essential in producing gamma-ray bursts, LIGO observed rates compared to GRB rates could be used to constrain the equation of state for dense nuclear matter., Comment: 12 pages, submitted to ApJ

We perform a systematic assessment of models for the equation of state (EOS) of dense matter in the context of recent neutron star mass and radius measurements to obtain a broad picture of the structure of neutron stars. We demonstrate that currently available neutron star mass and radius measurements provide strong constraints on moments of inertia, tidal deformabilities, and crust thicknesses. A measurement of the moment of inertia of PSR J0737-3039A with 10% error, without any other information from observations, will constrain the EOS over a range of densities to within 50%$-$60%. We find tidal deformabilities between 0.6 and $6\times 10^{36}$ g cm$^{2}$ s$^{2}$ (to 95% confidence) for $M=1.4~\mathrm{M}_{\odot}$, and any measurement which constrains this range will provide an important constraint on dense matter. The crustal fraction of the moment of inertia can be as large as 10% for $M=1.4~\mathrm{M}_{\odot}$ permitting crusts to have a large enough moment of inertia reservoir to explain glitches in the Vela pulsar even with a large amount of superfluid entrainment. Finally, due to the uncertainty in the equation of state, there is at least a 40% variation in the thickness of the crust for a fixed mass and radius, which implies that future simulations of the cooling of a neutron star crust which has been heated by accretion will need to take this variation into account., Comment: 7 pages, 2 figures

The minimal cooling paradigm for neutron star cooling assumes that enhanced cooling due to neutrino emission from any direct Urca process, due either to nucleons or to exotica such as hyperons, Bose condensates, or deconfined quarks, does not occur. This scenario was developed to replace and extend the so-called standard cooling scenario to include neutrino emission from the Cooper pair breaking and formation processes that occur near the critical temperature for superfluid/superconductor pairing. Recently, it has been found that Cooper-pair neutrino emission from the vector channel is suppressed by a large factor compared to the original estimates that violated vector current conservation. We show that Cooper-pair neutrino emission remains, nevertheless, an efficient cooling mechanism through the axial channel. As a result, the elimination of neutrino emission from Cooper-paired nucleons through the vector channel has only minor effects on the long-term cooling of neutron stars within the minimal cooling paradigm. We further quantify precisely the effect of the size of the neutron 3P2 gap and demonstrate that consistency between observations and the minimal cooling paradigm requires that the critical temperature T_c for this gap covers a range of values between T_c^min < 0.2 x 10^9 K up to T_c^max > 0.5 \times 10^9 K in the core of the star. In addition, it is required that young neutron stars have heterogenous envelope compositions: some must have light-element compositions and others must have heavy-element compositions. Unless these two conditions are fulfilled, about half of the observed young cooling neutron stars are inconsistent with the minimal cooling paradigm and provide evidence for the existence of enhanced cooling., 13 pages, 9 figures, submitted to ApJ

Constraints on the EoS for symmetric matter (equal neutron and proton numbers) at supra-saturation densities have been extracted from energetic collisions of heavy ions. Collisions of neutron-deficient and neutron-rich heavy ions now provide initial constraints on the EoS of neutron-rich matter at sub-saturation densities. Comparisons are made to other available constraints., Comment: 7 pages, 4 figures, Erice 2008 proceeding

We calculate the neutrino emissivity of superfluid neutron matter in the inner crust of neutron stars. We find that neutrino emission due to fluctuations resulting from the formation of Cooper pairs at finite temperature is highly suppressed in nonrelativistic systems. This suppression of the pair-breaking emissivity in a simplified model of neutron matter with interactions that conserve spin is of the order of v{sub F}{sup 4} for density fluctuations and v{sub F}{sup 2} for spin fluctuations, where v{sub F} is the Fermi velocity of neutrons. The larger suppression of density fluctuations arises because the dipole moment of the density distribution of a single component system does not vary in time. For this reason, we find that the axial current response (spin fluctuations) dominates. In more realistic models of neutron matter that include tensor interactions where the neutron spin is not conserved, neutrino radiation from bremsstrahlung reactions occurs at order v{sub F}{sup 0}. Consequently, even with the suppression factors due to superfluidity, this rate dominates near T{sub C}. Present calculations of the pair-breaking emissivity are incomplete because they neglect the tensor component of the nucleon-nucleon interaction.

The pervasive role of the nuclear symmetry energy in establishing some nuclear static and dynamical properties, and in governing some attributes of neutron star properties is highlighted., Comment: 8 pages, to appear in the proceedings of the International Symposium on Exotic States of Nuclear Matter (EXOCT 2007), Catania, Italy, 11-15 Jun 2007

The roles of isospin asymmetry in nuclei and neutron stars are investigated using a range of potential and field-theoretical models of nucleonic matter. The parameters of these models are fixed by fitting the properties of homogeneous bulk matter and closed-shell nuclei. We discuss and unravel the causes of correlations among the neutron skin thickness in heavy nuclei, the pressure of beta-equilibrated matter at a density of 0.1 fm$^{-3}$, and the radii of moderate mass neutron stars. The influence of symmetry energy on observables in heavy-ion collisions is summarized., Comment: 8 pages, 5 figures; Proceedings for the 21st Winter Workshop on Nuclear Dynamics, Breckenridge, Colorado, February 5-12, 2005; To appear in Heavy Ion Physics

A new classification of neutron star cooling scenarios, involving either ``minimal'' cooling or ``enhanced'' cooling is proposed. The minimal cooling scenario replaces and extends the so-called standard cooling scenario to include neutrino emission from the Cooper pair breaking and formation process. This emission dominates that due to the modified Urca process for temperatures close to the critical temperature for superfluid pairing. Minimal cooling is distinguished from enhanced cooling by the absence of neutrino emission from any direct Urca process, due either to nucleons or to exotica. Within the minimal cooling scenario, theoretical cooling models can be considered to be a four parameter family involving the equation of state of dense matter, superfluid properties of dense matter, the composition of the neutron star envelope, and the mass of the neutron star. Consequences of minimal cooling are explored through extensive variations of these parameters. Results are compared with the inferred properties of thermally-emitting neutron stars in order to ascertain if enhanced cooling occurs in any of them. All stars for which thermal emissions have been clearly detected are at least marginally consistent with the lack of enhanced cooling. The two pulsars PSR 0833-45 (Vela) and PSR 1706-44 would require enhanced cooling in case their ages and/or temperatures are on the lower side of their estimated values whereas the four stars PSR 0656+14, PSR 1055-52, Geminga, and RX J0720.4-3125 may require some source of internal heating in case their age and/or luminosity are on the upper side of their estimated values. The new upper limits on the thermal luminosity of PSR J0205+6449 and RX J0007.0+7302 are indicative of the occurrence of some enhanced neutrino emission beyond the minimal scenario., Comment: Version to appear in ApJ Supplements. Minor modifications in text and discussion of updated data with new figures

Recent Chandra observations have raised expectations that the objects RX J185635-3754 and 3C58 are either bare strange quark stars or stars with extended quark cores. However, these observations can also be interpreted in terms of normal neutron stars. Essential requirements for either explanation are that they simultaneously account for the observed (i) spectral features (i.e., thermal or nonthermal, possible lines), (ii) bounds on the inferred mass (M) and radius (R), and (iii) the cooling curves (effective temperatures T and luminosities L vs. age). Compact stars offer the promise of delineating QCD at finite baryon density (Lattimer and Prakash, 2001) in a fashion complementary to relativistic heavy ion collider experiments, which largely address the finite temperature, but baryon poor regime., Comment: 4 pages, 1 figure, uses espcrc1.sty (included)