The NICER Collaboration recently reported the measurement of the mass and radius of a pulsar PSR J0030+0451. We here use this new measurement to constrain one of the higher-order nuclear matter parameters $K_{\mathrm{sym,0}}$. We further combine the tidal measurement of the binary neutron star merger GW170817 by LIGO/Virgo to derive a joint 1-$\sigma$ constraint as $K_{\mathrm{sym,0}} = -102^{+71}_{-72}$ MeV. We believe this is the most reliable bound on the parameter to date under the assumption that there is no new physics above the saturation density which impacts neutron star observations.

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.

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.

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.

In this work, we study matter in the cores of proto-neutron stars, focusing on the impact of their composition on the stellar structure. We begin by examining the effects of finite temperature (through a fixed entropy per baryon) and lepton fraction on purely nucleonic matter by making use of the DSH model . We then turn our attention to a relativistic mean-field model containing exotic degrees of freedom, the Chiral Mean Field (CMF) model, again, under the conditions of finite temperature and trapped neutrinos. In the latter, since both hyperons and quarks are found in the cores of large-mass stars, their interplay and the possibility of mixtures of phases is taken into account and analyzed. Finally, we discuss how stellar rotation can affect our results.

The neutron star tidal deformability is a critical parameter which determines the pre-merger gravitational-wave signal in a neutron star merger. In this article, we show how neutron star tidal deformabilities behave in the presence of one or two sharp phase transition(s). We characterize how the tidal deformability changes when the properties of these phase transitions are modified in dense matter equation of state (EoS). Sharp phase transitions lead to the smallest possible tidal deformabilities and also induce discontinuities in the relation between tidal deformability and gravitational mass. These results are qualitatively unmodified by a modest softening of the phase transition. Finally, we test two universal relations involving the tidal deformability and show that their accuracy is limited by sharp phase transitions.

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.

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

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.

X-ray spectra of quiescent low-mass X-ray binaries containing neutron stars can be fit with atmosphere models to constrain the mass and the radius. Mass-radius constraints can be used to place limits on the equation of state of dense matter. We perform fits to the X-ray spectrum of a quiescent neutron star in the globular cluster M13, utilizing data from ROSAT, Chandra and XMM-Newton, and constrain the mass-radius relation. Assuming an atmosphere composed of hydrogen and a 1.4${\rm M}_{\odot}$ neutron star, we find the radius to be $R_{\rm NS}=12.2^{+1.5}_{-1.1}$ km, a significant improvement in precision over previous measurements. Incorporating an uncertainty on the distance to M13 relaxes the radius constraints slightly and we find $R_{\rm NS}=12.3^{+1.9}_{-1.7}$ km (for a 1.4${\rm M}_{\odot}$ neutron star with a hydrogen atmosphere), which is still an improvement in precision over previous measurements, some of which do not consider distance uncertainty. We also discuss how the composition of the atmosphere affects the derived radius, finding that a helium atmosphere implies a significantly larger radius.

We construct a new class of phenomenological equations of state for homogeneous matter for use in simulations of hot and dense matter in local thermodynamic equilibrium. We construct a functional form which respects experimental, observational and theoretical constraints on the nature of matter in various density and temperature regimes. Our equation of state matches (i) the virial coefficients expected from nucleon-nucleon scattering phase shifts, (ii) experimental measurements of nuclear masses and charge radii, (iii) observations of neutron star radii, (iv) theory results on the equation of state of neutron matter near the saturation density, and (v) theory results on the evolution of the EOS at finite temperatures near the saturation density. Our analytical model allows one to compute the variation in the thermodynamic quantities based on the uncertainties in the nature of the nucleon-nucleon interaction. Finally, we perform a correction to ensure the equation of state is causal at all densities, temperatures, and electron fractions.

Fitting models to data using Bayesian inference is quite common, but when each point in parameter space gives a curve, fitting the curve to a data set requires new nuisance parameters, which specify the metric embedding the one-dimensional curve into the higher-dimensional space occupied by the data. A generic formalism for curve fitting in the context of Bayesian inference is developed which shows how the aforementioned metric arises. The result is a natural generalization of previous works, and is compared to oft-used frequentist approaches and similar Bayesian techniques.

We present a quantitative analysis of superfluidity and superconductivity in dense matter from observations of isolated neutron stars in the context of the minimal cooling model. Our new approach produces the best fit neutron triplet superfluid critical temperature, the best fit proton singlet superconducting critical temperature, and their associated statistical uncertainties. We find that the neutron triplet critical temperature is likely $2.09^{+4.37}_{-1.41} \times 10^{8}$ K and that the proton singlet critical temperature is $7.59^{+2.48}_{-5.81} \times 10^{9}$ K. However, we also show that this result only holds if the Vela neutron star is not included in the data set. If Vela is included, the gaps increase significantly to attempt to reproduce Vela's lower temperature given its young age. Further including neutron stars believed to have carbon atmospheres increases the neutron critical temperature and decreases the proton critical temperature. Our method demonstrates that continued observations of isolated neutron stars can quantitatively constrain the nature of superfluidity in dense matter.

Observations of thermonuclear X-ray bursts from accreting neutron stars (NSs) in low-mass X-ray binary systems can be used to constrain NS masses and radii. Most previous work of this type has set these constraints using Planck function fits as a proxy: both the models and the data are fit with diluted blackbody functions to yield normalizations and temperatures which are then compared against each other. Here, for the first time, we fit atmosphere models of X-ray bursting NSs directly to the observed spectra. We present a hierarchical Bayesian fitting framework that uses state-of-the-art X-ray bursting NS atmosphere models with realistic opacities and relativistic exact Compton scattering kernels as a model for the surface emission. We test our approach against synthetic data, and find that for data that are well-described by our model we can obtain robust radius, mass, distance, and composition measurements. We then apply our technique to Rossi X-ray Timing Explorer observations of five hard-state X-ray bursts from 4U 1702-429. Our joint fit to all five bursts shows that the theoretical atmosphere models describe the data well but there are still some unmodeled features in the spectrum corresponding to a relative error of 1-5% of the energy flux. After marginalizing over this intrinsic scatter, we find that at 68% credibility the circumferential radius of the NS in 4U 1702-429 is R = 12.4+-0.4 km, the gravitational mass is M=1.9+-0.3 Msun, the distance is 5.1 < D/kpc < 6.2, and the hydrogen mass fraction is X < 0.09.

We analyze observations of eight quiescent low-mass X-ray binaries in globular clusters and combine them to determine the neutron star mass-radius curve and the equation of state of dense matter. We determine the effect that several uncertainties may have on our results, including uncertainties in the distance, the atmosphere composition, the neutron star maximum mass, the neutron star mass distribution, the possible presence of a hotspot on the neutron star surface, and the prior choice for the equation of state of dense matter. We find that the radius of a 1.4 solar mass neutron star is most likely from 10 to 14 km and that tighter constraints are only possible with stronger assumptions about the nature of the neutron stars, the systematics of the observations, or the nature of dense matter. Strong phase transitions are preferred over other models and interpretations of the data with a Bayes factor of 8 or more, and in this case, the radius is likely smaller than 12 km. However, radii larger than 12 km are preferred if the neutron stars have uneven temperature distributions.

Thermal states of neutron stars in soft X-ray transients (SXRTs) are thought to be determined by "deep crustal heating" in the accreted matter that drives the quiescent luminosity and cooling via emission of photons and neutrinos from the interior. In this study, we assume a global thermal steady-state of the transient system and calculate the heating curves (quiescent surface luminosity as a function of mean accretion rate) predicted from theoretical models, taking into account variations in the equations of state, superfluidity gaps, thickness of the light element layer and a phenomenological description of the direct Urca threshold. We further provide a statistical analysis on the uncertainties in these parameters, and compare the overall results with observations of several SXRTs, in particular the two sources containing the coldest (SAX J1808.4-3658) and the hottest (Aql X-1) neutron stars. Interpretation of the observational data indicates that the direct Urca process is required for the most massive stars and also suggests small superfluid gaps.

Two low mass neutron stars, J0737-3039B and the companion to J1756-2251, show strong evidence of being formed from the collapse of an ONeMg core in an electron capture supernova (ECSN) or in an ultra-stripped iron core collapse supernova (FeCCSN). Using three different systematically generated sets of equations of state we explore the relationship between the moment of inertia of J0737-3039A and the binding energy of the two low mass neutron stars. We find this relationship, a less strict variant of the recently discovered I-Love-Q relations, is nevertheless more robust than a previously explored correlation between the binding energy and the slope of the nuclear symmetry energy L. We find that, if either J0737-3039B or the J1756-2251 companion were formed in an ECSN, no more than 0.06 solar masses could have been lost from the progenitor core, more than four times the mass loss predicted by current supernova modeling. A measurement of the moment of inertia of J0737-3039A to within 10% accuracy from pulsar timing, possible within a decade, can discriminate between formation scenarios such as ECSN or ultra-stripped FeCCSN and, given current constraints on the predicted core mass loss, potentially rule them out. Using the I-Love-Q relations we find that an Advanced LIGO can potentially measure the neutron star tidal polarizability to equivalent accuracy in a neutron star-neutron star merger at a distance of 200 Mpc, thus obtaining similar constraints on the formation scenarios. Such information on the occurrence of ECSNe is important for population synthesis calculations, especially for estimating the rate of binary neutron star mergers and resulting electromagnetic and gravitational wave signals. Further progress needs to be made modeling the core collapse process that leads to low-mass neutron stars, particularly in making robust predictions for the mass loss from the progenitor core.

The bulk of the rare earth elements are believed to be synthesized in the rapid neutron capture process or $r$ process of nucleosynthesis. The solar $r$-process residuals show a small peak in the rare earths around $A\sim 160$, which is proposed to be formed dynamically during the end phase of the $r$ process by a pileup of material. This abundance feature is of particular importance as it is sensitive to both the nuclear physics inputs and the astrophysical conditions of the main $r$ process. We explore the formation of the rare earth peak from the perspective of an inverse problem, using Monte Carlo studies of nuclear masses to investigate the unknown nuclear properties required to best match rare earth abundance sector of the solar isotopic residuals. When nuclear masses are changed, we recalculate the relevant $\beta$-decay properties and neutron capture rates in the rare earth region. The feedback provided by this observational constraint allows for the reverse engineering of nuclear properties far from stability where no experimental information exists. We investigate a range of astrophysical conditions with this method and show how these lead to different predictions in the nuclear properties influential to the formation of the rare earth peak. We conclude that targeted experimental campaigns in this region will help to resolve the type of conditions responsible for the production of the rare earth nuclei, and will provide new insights into the longstanding problem of the astrophysical site(s) of the $r$ process.

Background: The nuclear symmetry energy is a fundamental ingredient in determining the equation of state (EOS) of neutron stars (NS). Recent terrestrial experiments constrain both its value and slope at nuclear saturation density, however, its value at higher densities is unknown. Assuming a Free Fermi-gas (FFG) model for the kinetic symmetry energy, the high-density extrapolation depends on a single parameter, the density dependence of the potential symmetry energy. The Correlated Fermi-gas (CFG) model improves on the FFG model by including the effects of short-range, correlated, high-momentum, nucleons in nuclear matter. Using the CFG model for the kinetic symmetry energy along with constraints from terrestrial measurements leads to a much softer density dependence for the potential symmetry energy. Purpose: Examine the ability of the FFG and CFG models to describe NS observables that are directly sensitive to the symmetry energy at high-density. Specifically, examine the ability of the CFG model, with its softer density dependence of the potential symmetry energy, to describe a two solar-mass NS. Methods: Using a Bayesian analysis of NS observables we compare the CFG and FFG models and examine the resulting parameters in the NS EOS and the density dependence of the potential symmetry energy. Results: Despite the large difference in the density dependence of the potential part of the symmetry energy, both models can describe the NS data and support a two solar-mass NS. The different density dependences has only a small effect on the NS EOS. Conclusions: While sensitive to the high-density values of the symmetry energy, NS observables alone are not enough to distinguish between the CFG and FFG models. This indicates that the NS EOS, obtain from Bayesian analysis of NS observables, is robust and is not sensitive to the exact nuclear model used for the kinetic symmetry energy.

This white paper informs the nuclear astrophysics community and funding agencies about the scientific directions and priorities of the field and provides input from this community for the 2015 Nuclear Science Long Range Plan. It summarizes the outcome of the nuclear astrophysics town meeting that was held on August 21-23, 2014 in College Station at the campus of Texas A&M University in preparation of the NSAC Nuclear Science Long Range Plan. It also reflects the outcome of an earlier town meeting of the nuclear astrophysics community organized by the Joint Institute for Nuclear Astrophysics (JINA) on October 9- 10, 2012 Detroit, Michigan, with the purpose of developing a vision for nuclear astrophysics in light of the recent NRC decadal surveys in nuclear physics (NP2010) and astronomy (ASTRO2010). The white paper is furthermore informed by the town meeting of the Association of Research at University Nuclear Accelerators (ARUNA) that took place at the University of Notre Dame on June 12-13, 2014. In summary we find that nuclear astrophysics is a modern and vibrant field addressing fundamental science questions at the intersection of nuclear physics and astrophysics. These questions relate to the origin of the elements, the nuclear engines that drive life and death of stars, and the properties of dense matter. A broad range of nuclear accelerator facilities, astronomical observatories, theory efforts, and computational capabilities are needed. With the developments outlined in this white paper, answers to long standing key questions are well within reach in the coming decade.

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.

The primary astrophysical source of the rare earth elements is the rapid neutron capture process ($r$ process). The rare earth peak that is seen in the solar $r$-process residuals has been proposed to originate as a pile-up of nuclei during the end of the $r$ process. We introduce a new method utilizing Monte Carlo studies of nuclear masses in the rare earth region, that includes self-consistently adjusting $\beta$-decay rates and neutron capture rates, to find the mass surfaces necessary for the formation of the rare earth peak. We demonstrate our method with two types of astrophysical scenarios, one corresponding conditions typical of core-collapse supernova winds and one corresponding to conditions typical of the ejection of the material from the tidal tails of neutron star mergers. In each type of astrophysical conditions, this method successfully locates a region of enhanced stability in the mass surface that is responsible for the rare earth peak. For each scenario, we find that the change in the mass surface has qualitatively different features, thus future measurements can shed light on the type of environment in which the $r$ process occurred.

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]

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.

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.

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.

There are several assumptions made in a standard $\chi^2$ analysis of data, including the frequent assumption that the likelihood function is well approximated by a multivariate Gaussian distribution. This article briefly reviews the standard approach and describes how Bayesian inference can be used to go beyond the assumption that the likelihood is Gaussian. Two separate types of analysis relevant to nuclear physics are used as test cases. The first is the determination of the equation of state of dense matter from neutron star mass and radius data. The second is the use of theoretical nuclear mass models to fit currently available data and predict the value of masses which have not yet been measured. For the problem of predicting nuclear masses, it is demonstrated that approximating the likelihood function with a Gaussian can produce biased predictions of unmeasured masses. Finally, the lessons learned from these fitting problems are used to propose a method for improving constraints on the nuclear symmetry energy.

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.

If the observed quasi-periodic oscillations in magnetar flares are partially confined to the crust, then the oscillation frequencies are unique probes of the nuclear physics of the neutron star crust. We study crustal oscillations in magnetars including corrections for a finite Alfven velocity. Our crust model uses a new nuclear mass formula that predicts nuclear masses with an accuracy very close to that of the finite range droplet model. This mass model for equilibrium nuclei also includes shell corrections and an updated neutron-drip line. We perturb our crust model to predict axial crust modes and assign them to observed giant flare quasi-periodic oscillation frequencies from the soft gamma-ray repeater SGR 1806-20. We find magnetar crusts that match observations for various magnetic field strengths, entrainment of the free neutron gas in the inner crust, and crust-core transition densities. We find that observations can be reconciled with smaller values of the symmetry energy slope parameter, L, if there is a significant amount of entrainment of the neutrons by the superfluid or if the crust-core transition density is large. We also find neutron star masses and radii which are in agreement with expectations from what is known about low-density matter from nuclear experiment. Matching observations with a field-free model we obtain the approximate values of M =1.35 Msun and R = 11.9 km. Matching observations using a model with the surface dipole field of SGR 1806-20 (B=2.4x10^15 G) we obtain the approximate values of M = 1.25 Msun and R = 12.4 km.

Neutron stars provide a fertile environment for exploring superfluidity under extreme conditions. It is not surprising that Cooper pairing occurs in dense matter since nucleon pairing is observed in nuclei as energy differences between even-even and odd-even nuclei. Since superfluids and superconductors in neutron stars profoundly affect neutrino emissivities and specific heats, their presence can be observed in the thermal evolution of neutron stars. An ever-growing number of cooling neutron stars, now amounting to 13 thermal sources, and several additional objects from which upper limits to temperatures can be ascertained, can now be used to discriminate among theoretical scenarios and even to dramatically restrict properties of nucleon pairing at high densities. In addition, observations of pulsars, including their spin-downs and glitch histories, additionally support the conjecture that superfluidity and superconductivity are ubiquitous within, and important to our understanding of, neutron stars.

We perform a systematic analysis of neutron star radius constraints from five quiescent low-mass X-ray binaries and examine how they depend on measurements of their distances and amounts of intervening absorbing material, as well as their assumed atmospheric compositions. We construct and calibrate to published results a semi-analytic model of the neutron star atmosphere which approximates these effects for the predicted masses and radii. Starting from mass and radius probability distributions established from hydrogen-atmosphere spectral fits of quiescent sources, we apply this model to compute alternate sets of probability distributions. We perform Bayesian analyses to estimate neutron star mass-radius curves and equation of state (EOS) parameters that best-fit each set of distributions, assuming the existence of a known low-density neutron star crustal EOS, a simple model for the high-density EOS, causality, and the observation that the neutron star maximum mass exceeds $2~M_\odot$. We compute the posterior probabilities for each set of distance measurements and assumptions about absorption and composition. We find that, within the context of our assumptions and our parameterized EOS models, some absorption models are disfavored. We find that neutron stars composed of hadrons are favored relative to those with exotic matter with strong phase transitions. In addition, models in which all five stars have hydrogen atmospheres are found to be weakly disfavored. Our most likely models predict neutron star radii that are consistent with current experimental results concerning the nature of the nucleon-nucleon interaction near the nuclear saturation density.

The nuclear symmetry energy is intimately connected with nuclear astrophysics. This contribution focuses on the estimation of the symmetry energy from experiment and how it is related to the structure of neutron stars. The most important connection is between the radii of neutron stars and the pressure of neutron star matter in the vicinity of the nuclear saturation density $n_s$. This pressure is essentially controlled by the nuclear symmetry energy parameters $S_v$ and $L$, the first two coefficients of a Taylor expansion of the symmetry energy around $n_s$. We discuss constraints on these parameters that can be found from nuclear experiments. We demonstrate that these constraints are largely model-independent by deriving them qualitatively from a simple nuclear model. We also summarize how recent theoretical studies of pure neutron matter can reinforce these constraints. To date, several different astrophysical measurements of neutron star radii have been attempted. Attention is focused on photospheric radius expansion bursts and on thermal emissions from quiescent low-mass X-ray binaries. While none of these observations can, at the present time, determine individual neutron star radii to better than 20% accuracy, the body of observations can be used with Bayesian techniques to effectively constrain them to higher precision. These techniques invert the structure equations and obtain estimates of the pressure-density relation of neutron star matter, not only near $n_s$, but up to the highest densities found in neutron star interiors. The estimates we derive for neutron star radii are in concordance with predictions from nuclear experiment and theory.

The temperature in the crust of an accreting neutron star, which comprises its outermost kilometer, is set by heating from nuclear reactions at large densities, neutrino cooling, and heat transport from the interior. The heated crust has been thought to affect observable phenomena at shallower depths, such as thermonuclear bursts in the accreted envelope. Here we report that cycles of electron capture and its inverse, $\beta^-$ decay, involving neutron-rich nuclei at a typical depth of about 150 m, cool the outer neutron star crust by emitting neutrinos while also thermally decoupling the surface layers from the deeper crust. This Urca mechanism has been studied in the context of white dwarfs and Type Ia supernovae, but hitherto was not considered in neutron stars, because previous models computed the crust reactions using a zero-temperature approximation and assumed that only a single nuclear species was present at any given depth. This thermal decoupling means that X-ray bursts and other surface phenomena are largely independent of the strength of deep crustal heating. The unexpectedly short recurrence times, of the order of years, observed for very energetic thermonuclear superbursts are therefore not an indicator of a hot crust, but may point instead to an unknown local heating mechanism near the neutron star surface.

Recent progress in quantum Monte Carlo with modern nucleon-nucleon interactions have enabled the successful description of properties of light nuclei and neutron-rich matter. As a demonstration, we show that the agreement between theoretical calculations of the charge form factor of 12C and the experimental data is excellent. Applying similar methods to isospin-asymmetric systems allows one to describe neutrons confined in an external potential and homogeneous neutron-rich matter. Of particular interest is the nuclear symmetry energy, the energy cost of creating an isospin asymmetry. Combining these advances with recent observations of neutron star masses and radii gives insight into the equation of state of neutron-rich matter near and above the saturation density. In particular, neutron star radius measurements constrain the derivative of the symmetry energy.

We review the calculation of the equation of state of pure neutron matter using quantum Monte Carlo (QMC) methods. QMC algorithms permit the study of many-body nuclear systems using realistic two- and three-body forces in a nonperturbative framework. We present the results for the equation of state of neutron matter, and focus on the role of three-neutron forces at supranuclear density. We discuss the correlation between the symmetry energy, the neutron star radius and the symmetry energy. We also combine QMC and theoretical models of the three-nucleon interactions, and recent neutron star observations to constrain the value of the symmetry energy and its density dependence.

White dwarfs, neutron stars and stellar mass black holes are key laboratories to study matter in most extreme conditions of gravity and magnetic field. The unprecedented effective area of Athena+ will allow us to advance our understanding of emission mechanisms and accretion physics over a wide range of mass accretion rates, starting from lower and sub-luminous quiescent X-ray binaries up to super-Eddington ultra-luminous sources. Athena+ will measure stellar black hole spins in a much higher number of binaries than achievable now, opening the possibility to study how spin varies with black hole history. The high throughput and energy resolution of the X-IFU will be instrumental in establishing how disc wind properties depend on accretion state, in determining wind launching mechanism and in quantifying the impact of the wind induced mass loss on binary evolution and environment. Triggers and high quality optical and radio data originating from large wide field contemporaneous instruments will provide essential complementary information on jet launching mechanisms and on the physics of rotation powered pulsars, for instance. In addition, Athena+ will furnish multiple, independent measurements of the neutron star mass/radius relation in a wide range of environments and conditions so as to constrain the debated equation of state.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

Many of the currently available equations of state for core-collapse supernova simulations give large neutron star radii and do not provide large enough neutron star masses, both of which are inconsistent with some recent neutron star observations. In addition, one of the critical uncertainties in the nucleon-nucleon interaction, the nuclear symmetry energy, is not fully explored by the currently available equations of state. In this article, we construct two new equations of state which match recent neutron star observations and provide more flexibility in studying the dependence on nuclear matter properties. The equations of state are also provided in tabular form, covering a wide range in density, temperature and asymmetry, suitable for astrophysical simulations. These new equations of state are implemented into our spherically symmetric core-collapse supernova model, which is based on general relativistic radiation hydrodynamics with three-flavor Boltzmann neutrino transport. The results are compared with commonly used equations of state in supernova simulations of 15 and 40 solar mass progenitors. We do not find any simple correlations between individual nuclear matter properties at saturation and the outcome of these simulations. However, the new equations of state lead to the most compact neutron stars among the relativistic mean-field models which we considered. The new models also obey the previously observed correlation between the time to black hole formation and the maximum mass of an s=4 neutron star.

The equation of state (EOS) of dense matter has been a long-sought goal of nuclear physics. Equations of state generate unique mass versus radius (M-R) relations for neutron stars, the ultra-dense remnants of stellar evolution. In this work, we determine the neutron star mass-radius relation and, based on recent observations of both transiently accreting and bursting sources, we show that the radius of a 1.4 solar mass neutron star lies between 10.4 and 12.9 km, independent of assumptions about the composition of the core. We show, for the first time, that these constraints remain valid upon removal from our sample of the most extreme transient sources or of the entire set of bursting sources; our constraints also apply even if deconfined quark matter exists in the neutron star core. Our results significantly constrain the dense matter EOS and are, furthermore, consistent with constraints from both heavy-ion collisions and theoretical studies of neutron matter. We predict a relatively weak dependence of the symmetry energy on the density and a value for the neutron skin thickness of lead which is less than 0.20 fm, results that are testable in forthcoming experiments.

A quasi-statistical equilibrium model is constructed to simulate the multicomponent composition of the crust of an accreting neutron star. The ashes of rp-process nucleosynthesis are driven by accretion through a series of electron captures, neutron emissions, and pycnonuclear fusions up to densities near the transition between the neutron star crust and core. A liquid droplet model which includes nuclear shell effects is used to provide nuclear masses far from stability. Reaction pathways are determined consistently with the nuclear mass model. The nuclear symmetry energy is an important uncertainty in the masses of the exotic nuclei in the inner crust and varying the symmetry energy changes the amount of deep crustal heating by as much as a factor of two.

The symmetry energy contribution to the nuclear Equation of State (EoS) impacts various phenomena in nuclear astrophysics, nuclear structure, and nuclear reactions. Its determination is a key objective of contemporary nuclear physics with consequences for the understanding of dense matter within neutron stars. We examine the results of laboratory experiments that have provided initial constraints on the nuclear symmetry energy and its density dependence at and somewhat below normal nuclear matter density. Some of these constraints have been derived from properties of nuclei. Others have been derived from the nuclear response to electroweak and hadronic probes. We also examine the most frequently used theoretical models that predict the symmetry energy and its slope. By comparing existing constraints on the symmetry pressure to theories, we demonstrate how the contribution of the three-body force, an essential ingredient in neutron matter models, can be determined.

Using a phenomenological form of the equation of state of neutron matter near the saturation density which has been previously demonstrated to be a good characterization of quantum Monte Carlo simulations, we show that currently available neutron star mass and radius measurements provide a significant constraint on the equation of state of neutron matter. At higher densities we model the equation of state using polytropes and a quark matter model, and we show that our results do not change strongly upon variation of the lower boundary density where these polytropes begin. Neutron star observations offer an important constraint on a coefficient which is directly connected to the strength of the three-body force in neutron matter, and thus some theoretical models of the three-body may be ruled out by currently available astrophysical data. In addition, we obtain an estimate of the symmetry energy of nuclear matter and its slope that can be directly compared to the experiment and other theoretical calculations.

The supernova remnant Cassiopeia A contains the youngest known neutron star which is also the first one for which real time cooling has ever been observed. In order to explain the rapid cooling of this neutron star, we first present the fundamental properties of neutron stars that control their thermal evolution with emphasis on the neutrino emission processes and neutron/proton superfluidity/superconductivity. Equipped with these results, we present a scenario in which the observed cooling of the neutron star in Cassiopeia A is interpreted as being due to the recent onset of neutron superfluidity in the core of the star. The manner in which the earlier occurrence of proton superconductivity determines the observed rapidity of this neutron star's cooling is highlighted. This is the first direct evidence that superfluidity and superconductivity occur at supranuclear densities within neutron stars.

We propose that the observed cooling of the neutron star in Cassiopeia A is due to enhanced neutrino emission from the recent onset of the breaking and formation of neutron Cooper pairs in the 3P2 channel. We find that the critical temperature for this superfluid transition is ~0.5x10^9 K. The observed rapidity of the cooling implies that protons were already in a superconducting state with a larger critical temperature. Our prediction that this cooling will continue for several decades at the present rate can be tested by continuous monitoring of this neutron star.

We determine an empirical dense matter equation of state from a heterogeneous dataset of six neutron stars: three type I X-ray bursters with photospheric radius expansion, studied by Ozel et al., and three transient low-mass X-ray binaries. We critically assess the mass and radius determinations from the X-ray burst sources and show explicitly how systematic uncertainties, such as the photospheric radius at touchdown, affect the most probable masses and radii. We introduce a parameterized equation of state and use a Markov Chain Monte Carlo algorithm within a Bayesian framework to determine nuclear parameters such as the incompressibility and the density dependence of the bulk symmetry energy. Using this framework we show, for the first time, that these parameters, predicted solely on the basis of astrophysical observations, all lie in ranges expected from nuclear systematics and laboratory experiments. We find significant constraints on the mass-radius relation for neutron stars, and hence on the pressure-density relation of dense matter. The predicted symmetry energy and the equation of state near the saturation density are soft, resulting in relatively small neutron star radii around 11-12 km for M=1.4 Msun. The predicted equation of state stiffens at higher densities, however, and our preferred model for X-ray bursts suggests that the neutron star maximum mass is relatively large, 1.9-2.2 Msun. Our results imply that several commonly used equations of state are inconsistent with observations.

We show that the fundamental seismic shear mode, observed as a quasi-periodic oscillation in giant flares emitted by highly-magnetized neutron stars, is particularly sensitive to the nuclear physics of the crust. The identification of an oscillation at ~ 30 Hz as the fundamental crustal shear mode requires a nuclear symmetry energy that depends very weakly on density near saturation. If the nuclear symmetry energy varies more strongly with density, then lower frequency oscillations, previously identified as torsional Alfven modes of the fluid core, could instead be associated with the crust. If this is the case, then future observations of giant flares should detect oscillations at around 18 Hz. An accurate measurement of the neutron skin thickness of lead will also constrain the frequencies predicted by the model.

Observationally inferred superburst ignition depths are shallower than models predict. We address this discrepancy by reexamining the superburst trigger mechanism. We first explore the hypothesis of Kuulkers et al. that exothermic electron captures trigger superbursts. We find that all electron capture reactions are thermally stable in accreting neutron star oceans and thus are not a viable trigger mechanism. Fusion reactions other than 12C + 12C are infeasible as well since the possible reactants either deplete at much shallower depths or have prohibitively large Coulomb barriers. Thus we confirm the proposal of Cumming & Bildsten and Strohmayer & Brown that 12C + 12C triggers superbursts. We then examine the 12C + 12C fusion rate. The reaction cross-section is experimentally unknown at astrophysically relevant energies, but resonances exist in the 12C + 12C system throughout the entire measured energy range. Thus it is likely, and in fact has been predicted, that a resonance exists near the Gamow peak energy ~ 1.5 MeV. For such a hypothetical 1.5 MeV resonance, we derive both a fiducial value and upper limit to the resonance strength and find that such a resonance could decrease the theoretically predicted superburst ignition depth by up to a factor of 4; in this case, observationally inferred superburst ignition depths would accord with model predictions for a range of plausible neutron star parameters. Said differently, such a resonance would decrease the temperature required for unstable 12C ignition at a column depth 10^12 g/cm^2 from 6 x 10^8 K to 5 x 10^8 K. Determining the existence of a strong resonance in the Gamow window requires measurements of the 12C + 12C cross-section down to a center-of-mass energy near 1.5 MeV, which is within reach of the proposed DUSEL facility.

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.

This white paper, directed to the Stars and Stellar Evolution panel, has three objectives: 1) to provide the Astro2010 Decadal Survey with a vista into the goals of the nuclear physics and nuclear astrophysics community; 2) to alert the astronomical community of joint opportunities for discoveries at the interface between nuclear physics and astronomy; and 3) to delineate efforts in nuclear physics and describe the observational and theoretical advances in astrophysics necessary to make progress towards answering the following questions in the Nuclear Science 2007 Long Range Plan: 1) What is the origin and distribution of the elements? 2) What are the nuclear reactions that power stars and stellar explosions? 3) What is the nature of dense matter? The scope of this white paper concerns the specific area of "low energy" nuclear astrophysics. We define this as the area of overlap between astrophysics and the study of nuclear structure and reactions. Of the questions listed above, two -- What is the origin of the elements? and What is the nature of dense matter? -- were specifically listed in the National Academies Study, "Connecting Quarks with the Cosmos".

Neutrinos emitted during the collapse, bounce and subsequent explosion provide information about supernova dynamics. The neutrino spectra are determined by weak interactions with nuclei and nucleons in the inner regions of the star, and thus the neutrino spectra are determined by the composition of matter. The composition of stellar matter at temperature ranging from $T=1-3$ MeV and densities ranging from $10^{-5}$ to 0.1 times the saturation density is explored. We examine the single-nucleus approximation commonly used in describing dense matter in supernova simulations and show that, while the approximation is accurate for predicting the energy and pressure at most densities, it fails to predict the composition accurately. We find that as the temperature and density increase, the single nucleus approximation systematically overpredicts the mass number of nuclei that are actually present and underestimates the contribution from lighter nuclei which are present in significant amounts.

We determine characteristic timescales for the viscous damping of r-mode oscillations in rapidly rotating compact stars that contain quark matter. We present results for the color-flavor-locked (CFL) phase of dense quark matter, in which the up, down and strange quarks are gapped, as well as the normal (ungapped) quark phase. While the ungapped quark phase supports a temperature window between 10^8 K and 5x10^9 K where the r-mode is damped even for rapid rotation, the r-mode in a rapidly rotating pure CFL star is not damped in the temperature range 10^10 K - 10^11 K. Rotating hybrid stars with quark matter cores display an instability window whose width is determined by the amount of quark matter present, and they can have large spin frequencies outside this window. Except at high temperatures T > 10^10 K, the presence of a quark phase allows for larger critical frequencies and smaller spin-periods compared to rotating neutron stars. If low-mass X-ray binaries contain a large amount of ungapped or CFL quark matter, then our estimates of the r-mode instability suggest that there should be a population of rapidly rotating binaries at frequencies greater than 1000 Hz which have not yet been observed.

{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 non-relativistic 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_F^4$ for density fluctuations and $v_F^2$ for spin fluctuations, where $v_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 which include tensor interactions where the neutron spin is not conserved, neutrino radiation from bremsstrahlung reactions occurs at order $v_F^0$. Consequently, even with the suppression factors due to superfluidity, this rate dominates near $T_C$. Present calculations of the pair-breaking emissivity are incomplete because they neglect the tensor component of the nucleon-nucleon interaction.

A fully self-consistent model of the neutron star inner crust based upon models of the nucleonic equation of state at zero temperature is constructed. The results nearly match those of previous calculations of the inner crust given the same input equation of state. The extent to which the uncertainties in the symmetry energy, the compressibility, and the equation of state of low-density neutron matter affect the composition of the crust are examined. The composition and pressure of the crust is sensitive to the description of low-density neutron matter and the nuclear symmetry energy, and the latter dependence is non-monotonic, giving larger nuclei for moderate symmetry energies and smaller nuclei for more extreme symmetry energies. Future nuclear experiments may help constrain the crust and future astrophysical observations may constrain the nuclear physics input.

Heavy-ion reactions provide a unique means to investigate the equation of state (EOS) of neutron-rich nuclear matter, especially the density dependence of the nuclear symmetry energy $E_{sym}(\rho)$. The latter plays an important role in understanding many key issues in both nuclear physics and astrophysics. Recent analyses of heavy-ion reactions have already put a stringent constraint on the $E_{sym}(\rho)$ around the saturation density. This subsequently allowed us to constrain significantly the radii and cooling mechanisms of neutron stars as well as the possible changing rate of the gravitational constant G.

Neutron star radii are primarily determined by the pressure of isospin asymmetric matter which is proportional to the slope of the nuclear symmetry energy. Available terrestrial laboratory data on the isospin diffusion in heavy-ion reactions at intermediate energies constrain the slope of the symmetry energy. Using this constraint, we show that the radius (radiation radius) of a 1.4 solar mass neutron star is between 11.5 (14.4) and 13.6 (16.3) km.

The existence of deconfined quark matter in the superdense interior of neutron stars is a key question that has drawn considerable attention over the past few decades. Quark matter can comprise an arbitrary fraction of the star, from 0 for a pure neutron star to 1 for a pure quark star, depending on the equation of state of matter at high density. From an astrophysical viewpoint, these two extreme cases are generally expected to manifest different observational signatures. An intermediate fraction implies a hybrid star, where the interior consists of mixed or homogeneous phases of quark and nuclear matter, depending on surface and Coulomb energy costs, as well as other finite size and screening effects. In this brief review article, we discuss what we can deduce about quark matter in neutron stars in light of recent exciting developments in neutron star observations. We state the theoretical ideas underlying the equation of state of dense quark matter, including color superconducting quark matter. We also highlight recent advances stemming from re-examination of an old paradigm for the surface structure of quark stars and discuss possible evolutionary scenarios from neutron stars to quark stars, with emphasis on astrophysical observations.

The symmetry energy of nucleonic matter is usually assumed to be quadratic in the isospin density. While this may be justified at sub-saturation densities, there is no need to enforce this restriction at super-saturation densities. The presence of a quartic term can strongly modify the critical density for the direct Urca process which leads to faster cooling of neutron stars. Neutron star cooling predictions which lie below the observational data can, for some equations of state, be repaired with a quartic term which effectively turns off the direct Urca process.

We construct strangelets, taking into account electrostatic effects, including Debye screening, and arbitrary surface tension sigma of the interface between vacuum and quark matter. We find that there is a critical surface tension sigma_crit below which large strangelets are unstable to fragmentation and below which quark star surfaces will fragment into a crystalline crust made of charged strangelets immersed in an electron gas. We derive a model-independent relationship between sigma_crit and two parameters that characterize any quark matter equation of state. For reasonable model equations of state, we find sigma_crit typically of order a few MeV/fm^2. If sigma <= sigma_crit, the size-distribution of strangelets in cosmic rays could feature a peak corresponding to the stable strangelets that we construct.

Nuclear reactions induced by stable and/or radioactive neutron-rich nuclei provide the opportunity to pin down the equation of state of neutron-rich matter, especially the density ($\rho$) dependence of its isospin-dependent part, i.e., the nuclear symmetry energy $E_{\rm sym}$. A conservative constraint, $32(\rho /\rho_{0})^{0.7} < E_{\rm sym}(\rho ) < 32(\rho /\rho _{0})^{1.1}$, around the nuclear matter saturation density $\rho_0$ has recently been obtained from the isospin diffusion data in intermediate energy heavy-ion collisions. We review this exciting result and discuss its consequences and implications on nuclear effective interactions, radii and cooling mechanisms of neutron stars.

The correlation between the thickness of the neutron skin in Pb-208, and the degree of isospin diffusion in heavy-ion collisions is examined. The same equation of state is used to compute the degree of isospin diffusion in an isospin-depedent transport model and the neutron skin thickness in the Hartree-Fock approximation. We find that skin thicknesses less than 0.15 fm are excluded by the isospin diffusion data.

We consider the effect of a six-fermion interaction of the {'}t Hooft form in the quark-quark channel on the ground state of matter at finite density. The coupling constant for this new term is varied within the limits suggested by naturalness. The flavor-mixing effects of the additional term destabilize the color-flavor-locked (CFL) and, to a lesser extent, the two-color superconducting (2SC) phases of quark matter, especially for positive values of the coupling. For some values of the coupling, the critical density for the CFL phase is nearly larger than the maximum density in the neutron stars. We comment on the implications for neutron star evolution.

We reexamine the surface composition of strange stars. Strange quark stars are hypothetical compact stars which could exist if strange quark matter was absolutely stable. It is widely accepted that they are characterized by an enormous density gradient ($~10^{26}$ g/cm$^4$) and large electric fields at surface. By investigating the possibility of realizing a heterogeneous crust, comprised of nuggets of strange quark matter embedded in an uniform electron background, we find that the strange star surface has a much reduced density gradient and negligible electric field. We comment on how our findings will impact various proposed observable signatures for strange stars.

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}$, the derivative of the nuclear symmetry energy at the same density and the radii of moderate mass neutron stars. Constraints on the symmetry properties of nuclear matter from the binding energies of nuclei are examined. The extent to which forthcoming neutron skin measurements will further delimit the symmetry properties is investigated. The impact of symmetry energy constraints for the mass and moment of inertia contained within neutron star crusts and the threshold density for the nucleon direct Urca process, all of which are potentially measurable, is explored. We also comment on the minimum neutron star radius, assuming that only nucleonic matter exists within the star.

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.

We investigate the consequences of enforcing local color neutrality on the color superconducting phases of quark matter by utilizing the Nambu-Jona-Lasinio model supplemented by diquark and the t'Hooft six-fermion interactions. In neutrino free matter at zero temperature, color neutrality guarantees that the number densities of u, d, and s quarks in the Color-Flavor-Locked (CFL) phase will be equal even with physical current quark masses. Electric charge neutrality follows as a consequence and without the presence of electrons. In contrast, electric charge neutrality in the less symmetric 2-flavor superconducting (2SC) phase with ud pairing requires more electrons than the normal quark phase. The free energy density cost of enforcing color and electric charge neutrality in the CFL phase is lower than that in the 2SC phase, which favors the formation of the CFL phase. With increasing temperature and neutrino content, an unlocking transition occurs from the CFL phase to the 2SC phase with the order of the transition depending on the temperature, the quark and lepton number chemical potentials. The astrophysical implications of this rich structure in the phase diagram, including estimates of the effects from Goldstone bosons in the CFL phase, are discussed.

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.

Neutrino fluxes from proto-neutron stars with and without quarks are studied. Observable differences become apparent after 10--20 s of evolution. Sufficiently massive stars containing negatively-charged, strongly interacting, particles collapse to black holes during the first minute of evolution. Since the neutrino flux vanishes when a black hole forms, this is the most obvious signal that quarks (or other types of strange matter) have appeared. The metastability timescales for stars with quarks are intermediate between those containing hyperons and kaon condensates.

Neutrino opacities important in the evolution of a proto-neutron star containing quark matter are studied. The results for pure quark matter are compared with limiting expressions previously derived, and are generalized to the temperatures, neutrino degeneracies and lepton contents encountered in a proto-neutron star's evolution. We find that the appearance of quarks in baryonic matter drastically reduces the neutrino opacity for a given entropy, the reduction being sensitive to the thermodynamic conditions in the mixed quark-hadron phase.

The main stages in the evolution of a neutron star, from its birth as a proto-neutron star, to its old age as a cold, catalyzed configuration, are described. A proto-neutron star is formed in the aftermath of a successful supernova explosion and its evolution is dominated by neutrino diffusion. Its neutrino signal is a valuable diagnostic of its internal structure and composition. During its transformation from a hot, lepton-rich to a cold, catalyzed remnant, the possibility exists that it can collapse into a black hole, which abruptly terminates neutrino emissions. The essential microphysics, reviewed herein, that controls its evolution are the equation of state of dense matter and its associated neutrino opacities. Several simulations of the proto-neutron star evolution, involving different assumptions about the composition of dense matter, are described. After its evolution into a nearly isothermal neutron star a hundred or so years after its birth, it may be observable through its thermal emission in X-rays during its life in the next million years. Its surface temperature will depend upon the rapidity of neutrino emission processes in its core, which depends on the composition of dense matter and whether or not its constituents exhibit superfluidity and superconductivity. Observations of thermal emission offer the best hope of a determination of the radius of a neutron star. The implications for the underlying dense matter equation of state of an accurate radius determination are explored.

Baryon and quark superfluidity in the cooling of neutron stars are investigated. Observations could constrain combinations of the neutron or Lambda-hyperon pairing gaps and the star's mass. However, in a hybrid star with a mixed phase of hadrons and quarks, quark gaps larger than a few tenths of an MeV render quark matter virtually invisible for cooling. If the quark gap is smaller, quark superfluidity could be important, but its effects will be nearly impossible to distinguish from those of other baryonic constituents.

The mixed phase of quarks and hadrons which might exist in the dense matter encountered in the varying conditions of temperature and trapped neutrino fraction in proto-neutron stars is studied. The extent that the mixed phase depends upon the thermodynamical parameters as well as on the stiffness of matter in the hadronic and quark phases is discussed. We show that hadronic equations of state that maximize the quark content of matter at a given {\it density} generally minimize the extent of the mixed phase region in a neutron star of a given mass, and that only in extreme cases could a pure quark star result. For both the Nambu Jona-Lasinio and MIT bag quark models, neutrino trapping inhibits the appearance of a mixed phase which leads to possible proto-neutron star metastability. The main difference between the two quark models is the small abundance of strange quarks in the former. We also demonstrate that $\partial T/\partial n<0$ along adiabats in the quark-hadron mixed phase, opposite to what is found for the kaon condensates-hadron mixed phase. This could lead to core temperatures which are significantly lower in stars containing quarks than in those not containing quarks.