Quark matter may appear due to a hadronic-quark transition in the core of a hybrid star. Quarkyonic matter is an approach in which both quarks and nucleons appear as quasi-particles in a crossover transition, and provides an explicit realization of early ideas concerning quark matter (e.g., the MIT bag model). This description has recently been employed by McLerran and Reddy to model chargeless (pure neutron) matter with an approach that has the virtue that the speed of sound rises quickly at a neutron-quark transition so as to satisfy observational constraints on the neutron star maximum mass ($\gtrsim2M_\odot$) and the radius of a $1.4M_\odot$ star ($R_{1.4}\lesssim 13.5$ km). Traditional models involving first-order transitions result in softer pressure-energy density relations that have difficulty satisfying these constraints except with very narrow choices of parameters. We propose a variation of quarkyonic matter involving protons and leptons whose energy can be explicitly minimized to achieve both chemical and beta equilibrium, which cannot be done in the chargeless formulation. Quarkyonic stellar models are able to satisfy observed mass and radius constraints with a wide range of model parameters, avoiding the obligatory fine-tuning of conventional hybrid star models, including requiring the transition density to be very close to the nuclear saturation density. Our formulation fits experimental and theoretical properties of the nuclear symmetry energy and pure neutron matter, and contains as few as three free parameters. This makes it an ideal tool for the study of high-density matter that is an efficient alternative to piecewise polytrope or spectral decomposition methods.

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

The possible detection of a compact object in the remnant of SN 1987A presents an unprecedented opportunity to follow its early evolution. The suspected detection stems from an excess of infrared emission from a dust blob near the compact object's predicted position. The infrared excess could be due to the decay of isotopes like 44Ti, accretion luminosity from a neutron star or black hole, magnetospheric emission or a wind originating from the spindown of a pulsar, or thermal emission from an embedded, cooling neutron star (NS 1987A). It is shown that the last possibility is the most plausible as the other explanations are disfavored by other observations and/or require fine-tuning of parameters. Not only are there indications the dust blob overlaps the predicted location of a kicked compact remnant, but its excess luminosity also matches the expected thermal power of a 30 year old neutron star. Furthermore, models of cooling neutron stars within the Minimal Cooling paradigm readily fit both NS 1987A and Cas A, the next-youngest known neutron star. If correct, a long heat transport timescale in the crust and a large effective stellar temperature are favored, implying relatively limited crustal n-1S0 superfluidity and an envelope with a thick layer of light elements, respectively. If the locations don't overlap, then pulsar spindown or accretion might be more likely, but the pulsar's period and magnetic field or the accretion rate must be rather finely tuned. In this case, NS 1987A may have enhanced cooling and/or a heavy-element envelope.

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

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

Neutron stars are not only of astrophysical interest, but are also of great interest to nuclear physicists, because their attributes can be used to determine the properties of the dense matter in their cores. One of the most informative approaches for determining the equation of state of this dense matter is to measure both a star's equatorial circumferential radius $R_e$ and its gravitational mass $M$. Here we report estimates of the mass and radius of the isolated 205.53 Hz millisecond pulsar PSR J0030+0451 obtained using a Bayesian inference approach to analyze its energy-dependent thermal X-ray waveform, which was observed using the Neutron Star Interior Composition Explorer (NICER). This approach is thought to be less subject to systematic errors than other approaches for estimating neutron star radii. We explored a variety of emission patterns on the stellar surface. Our best-fit model has three oval, uniform-temperature emitting spots and provides an excellent description of the pulse waveform observed using NICER. The radius and mass estimates given by this model are $R_e = 13.02^{+1.24}_{-1.06}$ km and $M = 1.44^{+0.15}_{-0.14}\ M_\odot$ (68%). The independent analysis reported in the companion paper by Riley et al. (2019) explores different emitting spot models, but finds spot shapes and locations and estimates of $R_e$ and $M$ that are consistent with those found in this work. We show that our measurements of $R_e$ and $M$ for PSR J0030$+$0451 improve the astrophysical constraints on the equation of state of cold, catalyzed matter above nuclear saturation density.

We report on Bayesian parameter estimation of the mass and equatorial radius of the millisecond pulsar PSR J0030$+$0451, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer (NICER) X-ray spectral-timing event data. We perform relativistic ray-tracing of thermal emission from hot regions of the pulsar's surface. We assume two distinct hot regions based on two clear pulsed components in the phase-folded pulse-profile data; we explore a number of forms (morphologies and topologies) for each hot region, inferring their parameters in addition to the stellar mass and radius. For the family of models considered, the evidence (prior predictive probability of the data) strongly favors a model that permits both hot regions to be located in the same rotational hemisphere. Models wherein both hot regions are assumed to be simply-connected circular single-temperature spots, in particular those where the spots are assumed to be reflection-symmetric with respect to the stellar origin, are strongly disfavored. For the inferred configuration, one hot region subtends an angular extent of only a few degrees (in spherical coordinates with origin at the stellar center) and we are insensitive to other structural details; the second hot region is far more azimuthally extended in the form of a narrow arc, thus requiring a larger number of parameters to describe. The inferred mass $M$ and equatorial radius $R_\mathrm{eq}$ are, respectively, $1.34_{-0.16}^{+0.15}$ M$_{\odot}$ and $12.71_{-1.19}^{+1.14}$ km, whilst the compactness $GM/R_\mathrm{eq}c^2 = 0.156_{-0.010}^{+0.008}$ is more tightly constrained; the credible interval bounds reported here are approximately the $16\%$ and $84\%$ quantiles in marginal posterior mass.

The LIGO/Virgo Consortium (LVC) released a preliminary announcement of a candidate gravitational wave signal, S190426c, that could have arisen from a black hole-neutron star merger. As the first such candidate system, it's properties such as masses and spin are of great interest. Although LVC policy prohibits disclosure of these properties in preliminary announcements, LVC does release the estimated probabilities that this system is in specific categories, such as binary neutron star, binary black hole and black hole-neutron star. LVC also releases information concerning relative signal strength, distance, and the probability that ejected mass or a remnant disc survived the merger. In the case of events with a finite probability of being in more than one category, such as is likely to occur with a black hole-neutron star merger, it is shown how to estimate the masses of the components and the spin of the black hole. This technique is applied to the source S190426c.

As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations.

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.

We use gravitational-wave observations of the binary neutron star merger GW170817 to explore the tidal deformabilities and radii of neutron stars. We perform Bayesian parameter estimation with the source location and distance informed by electromagnetic observations. We also assume that the two stars have the same equation of state; we demonstrate that for stars with masses comparable to the component masses of GW170817, this is effectively implemented by assuming that the stars' dimensionless tidal deformabilities are determined by the binary's mass ratio $q$ by $\Lambda_1/\Lambda_2 = q^6$. We investigate different choices of prior on the component masses of the neutron stars. We find that the tidal deformability and 90$\%$ credible interval is $\tilde{\Lambda}=222^{+420}_{-138}$ for a uniform component mass prior, $\tilde{\Lambda}=245^{+453}_{-151}$ for a component mass prior informed by radio observations of Galactic double neutron stars, and $\tilde{\Lambda}=233^{+448}_{-144}$ for a component mass prior informed by radio pulsars. We find a robust measurement of the common areal radius of the neutron stars across all mass priors of $8.9 \le \hat{R} \le 13.2$ km, with a mean value of $\langle \hat{R} \rangle = 10.8$ km. Our results are the first measurement of tidal deformability with a physical constraint on the star's equation of state and place the first lower bounds on the deformability and areal radii of neutron stars using gravitational waves.

Finite size effects in a neutron star merger are manifested, at leading order, through the tidal deformabilities (Lambdas) of the stars. If strong first-order phase transitions do not exist within neutron stars, both neutron stars are described by the same equation of state, and their Lambdas are highly correlated through their masses even if the equation of state is unknown. If, however, a strong phase transition exists between the central densities of the two stars, so that the more massive star has a phase transition and the least massive star does not, this correlation will be weakened. In all cases, a minimum Lambda for each neutron star mass is imposed by causality, and a less conservative limit is imposed by the unitary gas constraint, both of which we compute. In order to make the best use of gravitational wave data from mergers, it is important to include the correlations relating the Lambdas and the masses as well as lower limits to the Lambdas as a function of mass. Focusing on the case without strong phase transitions, and for mergers where the chirp mass M_chirp<1.4M_sun, which is the case for all observed double neutron star systems where a total mass has been accurately measured, we show that the dimensionless Lambdas satisfy Lambda_1/Lambda_2= q^6, where q=M_2/M_1 is the binary mass ratio; $M$ is mass of each star, respectively. Moreover, they are bounded by q^{n_-}>Lambda_1/Lambda_2> q^{n_{0+}+qn_{1+}}, where n_-<n_{0+}+qn_{1+}; the parameters depend only on M_chirp, which is accurately determined from the gravitational-wave signal. We also provide analytic expressions for the wider bounds that exist in the case of a strong phase transition. We argue that bounded ranges for Lambda_1/Lambda_2, tuned to M_chirp, together with lower bounds to Lambda(M), will be more useful in gravitational waveform modeling than other suggested approaches.

We propose the existence of a lower bound on the energy of pure neutron matter (PNM) on the basis of unitary-gas considerations. We discuss its justification from experimental studies of cold atoms as well as from theoretical studies of neutron matter. We demonstrate that this bound results in limits to the density-dependent symmetry energy, which is the difference between the energies of symmetric nuclear matter and PNM. In particular, this bound leads to a lower limit to the volume symmetry energy parameter $S_0$. In addition, for assumed values of $S_0$ above this minimum, this bound implies both upper and lower limits to the symmetry energy slope parameter $L$, which describes the lowest-order density dependence of the symmetry energy. A lower bound on the neutron-matter incompressibility is also obtained. These bounds are found to be consistent with both recent calculations of the energies of PNM and constraints from nuclear experiments. Our results are significant because several equations of state that are currently used in astrophysical simulations of supernovae and neutron star mergers, as well as in nuclear physics simulations of heavy-ion collisions, have symmetry energy parameters that violate these bounds. Furthermore, below the nuclear saturation density, the bound on neutron-matter energies leads to a lower limit to the density-dependent symmetry energy, which leads to upper limits to the nuclear surface symmetry parameter and the neutron-star crust-core boundary. We also obtain a lower limit to the neutron-skin thicknesses of neutron-rich nuclei. Above the nuclear saturation density, the bound on neutron-matter energies also leads to an upper limit to the symmetry energy, with implications for neutron-star cooling via the direct Urca process.

Recent developments in the theory of pure neutron matter and experiments concerning the symmetry energy of nuclear matter, coupled with recent measurements of high-mass neutron stars, now allow for relatively tight constraints on the equation of state of dense matter. We review how these constraints are formulated and describe the implications they have for neutron stars and core-collapse supernovae. We also examine thermal properties of dense matter, which are important for supernovae and neutron star mergers, but which cannot be nearly as well constrained at this time by experiment. In addition, we consider the role of the equation of state in medium-energy heavy-ion collisions.

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.

We explore the thermal properties of hot and dense matter using a model that reproduces the empirical properties of isospin symmetric and asymmetric bulk nuclear matter, optical model fits to nucleon-nucleus scattering data, heavy-ion flow data in the energy range 0.5-2 GeV/A, and the largest well-measured neutron star mass of 2 $\rm{M}_\odot$. Results of this model which incorporates finite range interactions through Yukawa type forces are contrasted with those of a zero-range Skyrme model that yields nearly identical zero-temperature properties at all densities for symmetric and asymmetric nucleonic matter and the maximum neutron star mass, but fails to account for heavy-ion flow data due to the lack of an appropriate momentum dependence in its mean field. Similarities and differences in the thermal state variables and the specific heats between the two models are highlighted. Checks of our exact numerical calculations are performed from formulas derived in the strongly degenerate and non-degenerate limits. Our studies of the thermal and adiabatic indices, and the speed of sound in hot and dense matter for conditions of relevance to core-collapse supernovae, the thermal evolution of neutron stars from their birth and mergers of compact binary stars reveal that substantial variations begin to occur at sub-saturation densities before asymptotic values are reached at supra-nuclear densities.

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.

We investigate the thermal properties of the potential model equation of state of Akmal, Pandharipande and Ravenhall. This equation of state approximates the microscopic model calculations of Akmal and Pandharipande, which feature a neutral pion condensate. We treat the bulk homogeneous phase for isospin asymmetries ranging from symmetric nuclear matter to pure neutron matter and for temperatures and densities relevant for simulations of core-collapse supernovae, proto-neutron stars, and neutron star mergers. Numerical results of the state variables are compared with those of a typical Skyrme energy density functional with similar properties at nuclear densities, but which differs substantially at supra-nuclear densities. Analytical formulas, which are applicable to non-relativistic potential models such as the equations of state we are considering, are derived for all state variables and their thermodynamic derivatives. A highlight of our work is its focus on thermal response functions in the degenerate and non-degenerate situations, which allow checks of the numerical calculations for arbitrary degeneracy. These functions are sensitive to the density dependent effective masses of neutrons and protons, which determine the thermal properties in all regimes of degeneracy. We develop the "thermal asymmetry free energy" and establish its relation to the more commonly used nuclear symmetry energy. We also explore the role of the pion condensate at supra-nuclear densities and temperatures. Tables of matter properties as functions of baryon density, composition (i.e., proton fraction) and temperature are being produced which are suitable for use in astrophysical simulations of supernovae and neutron stars.

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

Neutron stars are valuable laboratories for the study of dense matter. Recent observations have uncovered both massive and low-mass neutron stars and have also set constraints on neutron star radii. The largest mass measurements are powerfully influencing the high-density equation of state because of the existence of the neutron star maximum mass. The smallest mass measurements, and the distributions of masses, have implications for the progenitors and formation mechanisms of neutron stars. The ensemble of mass and radius observations can realistically restrict the properties of dense matter, and, in particular, the behavior of the nuclear symmetry energy near the nuclear saturation density. Simultaneously, various nuclear experiments are progressively restricting the ranges of parameters describing the symmetry properties of the nuclear equation of state. In addition, new theoretical studies of pure neutron matter are providing insights. These observational, experimental and theoretical constraints of dense matter, when combined, are now revealing a remarkable convergence.

One of the major uncertainties in the dense matter equation of state has been the nuclear symmetry energy. The density dependence of the symmetry energy is important in nuclear astrophysics, as it controls the neutronization of matter in core-collapse supernovae, the radii of neutron stars and the thicknesses of their crusts, the rate of cooling of neutron stars, and the properties of nuclei involved in r-process nucleosynthesis. We show that fits of nuclear masses to experimental masses, combined with other experimental information from neutron skins, heavy ion collisions, giant dipole resonances and dipole polarizabilities, lead to stringent constraints on parameters that describe the symmetry energy near the nuclear saturation density. These constraints are remarkably consistent with inferences from theoretical calculations of pure neutron matter, and, furthermore, with astrophysical observations of neutron stars. The concordance of experimental, theoretical and observational analyses suggests that the symmetry parameters S_v and L are in the range 29.0 - 32.7 MeV and 40.5 - 61.9 MeV, respectively, and that the neutron star radius, for a 1.4 M_sun star, is in the narrow window 10.7 km < R < 13.1 km (90% confidence). We can also set tight limits to the size of neutron star crusts and the fractional moment of inertia they contain, as well as the overall moment of inertia and quadrupole polarizability of 1.4 M_sun stars. Our results also have implications for the disk mass and ejected mass of compact mergers involving neutron stars.

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.

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.

The determination of neutron star masses is reviewed in light of a new measurement of 1.97 M$_\odot$ for PSR J1614-2230 and an estimate of 2.4 M$_\odot$ for the black widow pulsar. Using a simple analytic model related to the so-called maximally compact equation of state, model-independent upper limits to thermodynamic properties in neutron stars, such as energy density, pressure, baryon number density and chemical potential, are established which depend upon the neutron star maximum mass. Using the largest well-measured neutron star mass, 1.97 M$_\odot$, it is possible to show that the energy density can never exceed about 2 GeV, the pressure about 1.3 GeV, and the baryon chemical potential about 2.1 GeV. Further, if quark matter comprises a significant component of neutron star cores, these limits are reduced to 1.3 GeV, 0.9 GeV, and 1.5 GeV, respectively. We also find that the maximum binding energy of any neutron star is about 25% of the rest mass. Neutron matter properties and astrophysical constraints additionally imply an upper limit to the neutron star maximum mass of about 2.4 M$_\odot$. A measured mass of 2.4 M$\odot$ would be incompatible with hybrid star models containing {\it significant} proportions of exotica in the form of hyperons, Bose condensates or quark matter.

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

We 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 have redetermined the parallax and proper motion of the nearby isolated neutron star RX~J185635-3754. We used eight observations with the high resolution camera of the HST/ACS taken from 2002 through 2004. We performed the astrometric fitting using five independent methods, all of which yielded consistent results. The mean estimate of the distance is 123 (+11, -15) pc (1 sigma), in good agreement with our earlier published determination.

Gravitational waves from the final stages of inspiralling binary neutron stars are expected to be one of the most important sources for ground-based gravitational wave detectors. The masses of the components are determinable from the orbital and chirp frequencies during the early part of the evolution during which tidal effects provide small correction; however, during this phase the signal is relatively clean. The accumulated phase shift due to tidal corrections is characterized by a single quantity, the Love number, which is sensitive to the compactness parameter M/R and the star's internal structure, and its determination could constrain the star's radius. We show that the Love number of normal neutron stars are much different from those of self-bound strange quark matter stars and could therefore provide an important way to distinguish between these two classes of stars.

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.

The cores of neutron stars harbor the highest matter densities known to occur in nature, up to several times the densities in atomic nuclei. Similarly, magnetic field strengths can exceed the strongest fields generated in terrestrial laboratories by ten orders of magnitude. Hyperon-dominated matter, deconfined quark matter, superfluidity, even superconductivity are predicted in neutron stars. Similarly, quantum electrodynamics predicts that in strong magnetic fields the vacuum becomes birefringent. The properties of matter under such conditions is governed by Quantum Chromodynamics (QCD) and Quantum Electrodynamics (QED), and the close study of the properties of neutron stars offers the unique opportunity to test and explore the richness of QCD and QED in a regime that is utterly beyond the reach of terrestrial experiments. Experimentally, this is almost virgin territory.

The characteristic physical timescales near stellar-mass compact objects are measured in milliseconds. These timescales -- the free-fall time, the fastest stable orbital period, and stellar spin periods -- encode the fundamental physical properties of compact objects: mass, radius, and angular momentum. The characteristic temperature of matter in the vicinity of neutron stars is such that the principal electromagnetic window into their realms is the X-ray band. Because of these connections to the fundamental properties of neutron stars, X-ray timing studies remain today the most direct means of probing their structure and dynamics. While current X-ray observatories have revealed many relevant and fascinating phenomena, they lack the sensitivity to fully exploit them to uncover the fundamental properties of compact objects and their extreme physics. With this white paper, we summarize and highlight the science opportunities that will accompany an order-of-magnitude improvement in X-ray timing sensitivity, a goal attainable in the coming decade.

More than four decades after the discovery of pulsars, the composition of matter at their cores is still a mystery. This white paper summarizes how recent high-precision measurements of millisecond pulsar masses have introduced new experimental constraints on the properties of super-dense matter, and how continued timing of intriguing new objects, coupled with radio telescope surveys to discover more pulsars, might introduce significantly more stringent constraints.

We calculate Keplerian (mass shedding) configurations of rigidly rotating neutron stars and quark stars with crusts. We check the validity of empirical formula for Keplerian frequency, f_K, proposed by Lattimer & Prakash, f_K(M)=C (M/M_sun)^1/2 (R/10km)^-3/2, where M is the (gravitational) mass of Keplerian configuration, R is the (circumferential) radius of the non-rotating configuration of the same gravitational mass, and C = 1.04 kHz. Numerical calculations are performed using precise 2-D codes based on the multi-domain spectral methods. We use a representative set of equations of state (EOSs) of neutron stars and quark stars. We show that the empirical formula for f_K(M) holds within a few percent for neutron stars with realistic EOSs, provided 0.5 M_sun < M < 0.9 M_max,stat, where M_max,stat is the maximum allowable mass of non-rotating neutron stars for an EOS, and C=C_NS=1.08 kHz. Similar precision is obtained for quark stars with 0.5 M_sun < M < 0.9 M_max,stat. For maximal crust masses we obtain C_QS = 1.15 kHz, and the value of C_QS is not very sensitive to the crust mass. All our C's are significantly larger than the analytic value from the relativistic Roche model, C_Roche = 1.00 kHz. For 0.5 M_sun < M < 0.9 M_max,stat, the equatorial radius of Keplerian configuration of mass M, R_K(M), is, to a very good approximation, proportional to the radius of the non-rotating star of the same mass, R_K(M) = aR(M), with a_NS \approx a_QS \approx 1.44. The value of a_QS is very weakly dependent on the mass of the crust of the quark star. Both a's are smaller than the analytic value a_Roche = 1.5 from the relativistic Roche model.

We investigate how current and proposed observations of neutron stars can lead to an understanding of the state of their interiors and the key unknowns: the typical neutron star radius and the neutron star maximum mass. A theoretical analysis of neutron star structure, including general relativistic limits to mass, compactness, and spin rates is made. We consider observations made not only with photons, ranging from radio waves to X-rays, but also those involving neutrinos and gravity waves. We detail how precision determinations of structural properties would lead to significant restrictions on the poorly understood equation of state near and beyond the equilibrium density of nuclear matter.

Close binary systems of compact stars, due to the emission of gravitational radiation, may evolve into a phase in which the less massive star transfers mass to its companion. We describe mass transfer by using the model of Roche lobe overflow, in which mass is transferred through the first, or innermost, Lagrange point. Under conditions in which gravity is strong, the shapes of the equipotential surfaces and the Roche lobes are modified compared to the Newtonian case. We present calculations of the Roche lobe utilizing the second order post-Newtonian (2PN) approximation in the Arnowitt-Deser-Misner gauge. Heretofore, calculations of the Roche lobe geometry beyond the Newtonian case have not been available. Beginning from the general N-body Lagrangian derived by Damour and Schaffer, we develop the Lagrangian for a test particle in the vicinity of two massive compact objects. As an exact result for the transverse-traceless part of the Lagrangian is not available, we devise an approximation that is valid for regions close to the less massive star. We calculate the Roche lobe volumes, and provide a simple fitting formula for the effective Roche lobe radius analogous to that for the Newtonian case furnished by Eggleton. In contrast to the Newtonian case, in which the effective Roche radius depends only upon the mass ratio $q=m_1/m_2$, in the 2PN case the effective Roche lobe radius also depends on the ratio $z=2 (m_1+m_2)/a$ of the total mass and the orbital separation.

Binary mergers involving black holes and neutron stars have been proposed as major sources of gravitational waves, r--process nucleosynthesis, and gamma ray bursters. In addition, they represent an important, and possibly unique, observable that could distinguish between normal and self--bound neutron stars. These two families of stars have distinctly different mass--radius relationships resulting from their equations of state which can be revealed during their mergers if stable mass transfer ensues. We consider two cases of gravitational-radiation induced binary mergers: (i) a black hole and a normal neutron star, and (ii) a black hole and a self-bound strange quark matter star. We extend previous Newtonian analyses to incorporate the pseudo-general relativistic Paczy\'nski-Wiita potential or a potential correct to second--order post-Newtonian order in Arnowitt--Deser--Misner coordinates. These potentials are employed to study both the orbital evolution of the binary and the Roche lobe geometry that determines when and if stable mass transfer between the components is possible. The Roche lobe geometry with pseudo-general relativistic or post-Newtonian potentials has not heretofore been considered. Our analysis shows that differences in the evolution of normal neutron stars and strange quark matter stars are significant and could be detected in gravity waves. Both the amplitude and frequencies of the wave pattern are affected. In addition, details of the equation of state for either normal neutron stars or strange quark stars may be learned. A single merger could reveal one or two individual points of the mass-radius relation, and observations of several mergers could map a significant portion of this relation.

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.

We estimate that the moment of inertia of star A in the recently discovered double pulsar system PSR J0737-3039 may be determined after a few years of observation to something like 10% accuracy. This would enable accurate estimates of the radius of the star and the presure of matter in the vicinity of 1 to 2 times the nuclear saturation density, which would in turn provide strong constraints on the equation of state of neutron stars and the physics of their interiors.

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.

Neutron stars are some of the densest manifestations of massive objects in the universe. They are ideal astrophysical laboratories for testing theories of dense matter physics and provide connections among nuclear physics, particle physics and astrophysics. Neutron stars may exhibit conditions and phenomena not observed elsewhere, such as hyperon-dominated matter, deconfined quark matter, superfluidity and superconductivity with critical temperatures near ${10^{10}}$ kelvin, opaqueness to neutrinos, and magnetic fields in excess of $10^{13}$ Gauss. Here, we describe the formation, structure, internal composition and evolution of neutron stars. Observations that include studies of binary pulsars, thermal emission from isolated neutron stars, glitches from pulsars and quasi-periodic oscillations from accreting neutron stars provide information about neutron star masses, radii, temperatures, ages and internal compositions.

The mergers of black hole-neutron star binaries are calcuated using a pseudo-general relativistic potential that incorporates ${\mathcal O}(v^2/c^2)^3$ post-Newtonian corrections. Both normal matter neutron stars and self-bound strange quark matter stars are considered as black hole partners. As long as the neutron stars are not too massive relative to the black hole mass, orbital decay terminates in stable mass transfer rather than an actual merger. For a normal neutron star, mass transfer results in a widening of the orbit but the stable transfer ends before the minimum neutron star mass is reached. For a strange star, mass transfer does not result in an appreciable enlargement of the orbital separation, and the stable transfer continues until the strange star essentially disappears. These differences might be observable through their respective gravitational wave signatures.

We contrast the evolution of gravitational binary mergers for the two cases of a black hole and a normal neutron star, and a black hole and a self-bound strange quark matter star. In both cases inspiral continues until the Roche limit is reached, at which point it is expected that stable mass transfer to the balck hole ensues. Whereas a neutron star would then outspiral, the strange quark matter star barely moves out at all. Eventually, the strange quark star loses all its mass to the black hole, as compared to the neutron star whose outspiral continues until stable mass transfer terminates and inspiral resumes. These scenarios result in distinctly different gravitational wave signatures.

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.

New astrometric analysis of four WFPC2 images of the isolated neutron star RX J185635-3754 show that its distance is 117 +/- 12 pc, nearly double the originally published distance. At the revised distance, the star's age is 5 x 10^5 years, its space velocity is about 185 km/s, and its radiation radius inferred from thermal emission is approximately 15 km, in the range of many equations of state both with and without exotic matter. These measurements remove observational support for an extremely soft equation of state. The star's birthplace is still likely to be in the Upper Sco association, but a connection with zeta Oph is now unlikely.

We investigate the effects of very strong magnetic fields upon the equation of state of dense bayonic matter in which hyperons are present. In the presence of a magnetic field, the equation of state above nuclear density is significantly affected both by Landau quantization and magnetic moment interactions, but only for field strengths $B>5\times10^{18}$ G. The former tends to soften the EOS and increase proton and lepton abundances, while the latter produces an overall stiffening of the EOS. Each results in a supression of hyperons relative to the field-free case. The structure of a neutron star is, however, primarily determined by the magnetic field stress. We utilize existing general relativistic magneto-hydrostatic calculations to demonstrate that maximum average fields within a stable neutron are limited to values $B\le 1-3 \times10^{18}$ G. This is not large enough to significantly influence particle compositions or the matter pressure, unless fluctuations dominate the average field strengths in the interior or configurations with significantly larger field gradients are considered.

We discuss efforts to determine the mass, radius, and surface composition of the nearby compact object RX J185635-3754 from its multi-wavelength spectral energy distribution. We compute non-magnetized model atmospheres and emergent spectra for selected compositions and gravities, and discuss efforts to fit existing and new observational data from ROSAT, EUVE and the HST. The spectral energy distribution matches that expected from a heavy-element dominated atmosphere, but not from a uniform temperature blackbody. Non-magnetic light element atmospheres cannot be simultaneously reconciled with the optical and X-ray data. We extend previous studies, which were limited to one fixed neutron star mass and radius. For uniform temperature models dominated by heavy elements, the redshift z is constrained to be 0.3 < z < 0.4 and the best-fit mass and radius are M approx 0.9 solar masses and R approx 6 km (for a 61 pc distance). These values for M and R together are not permitted for any plausible equation of state, including that of a self-bound strange quark star. A simplified two-temperature model allows masses and radii up to about 50% larger, or a factor of 2 in the case of a black body. The observed luminosity is consistent with the thermal emission of an isolated neutron star no older than about 1 million years, the age inferred from available proper motion and parallax information

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.

Even the elusive neutrinos are trapped in matter, albeit transiently, in several astrophysical circumstances. Their interactions with the ambient matter not only reveal the properties of such exotic matter itself, but also shed light on the fundamental properties of the neutrinos themselves. The physical sites of interest include the early universe, supernovae, and newly-born neutron stars. Detection of neutrinos from these vastly different eras using the new generation of neutrino detectors holds great promise for enhancing our understanding of neutrino-matter interactions and astrophysical phenomena.

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.

We study static neutron stars with poloidal magnetic fields and a simple class of electric current distributions consistent with the requirement of stationarity. For this class of electric current distributions, we find that magnetic fields are too large for static configurations to exist when the magnetic force pushes a sufficient amount of mass off-center that the gravitational force points outward near the origin in the equatorial plane. (In our coordinates an outward gravitational force corresponds to $\partial\ln g_{tt}/\partial r>0$, where $t$ and $r$ are respectively time and radial coordinates and $g_{tt}$ is coefficient of $dt^2$ in the line element.) For the equations of state (EOSs) employed in previous work, we obtain configurations of higher mass than had been reported; we also present results with more recent EOSs. For all EOSs studied, we find that the maximum mass among these static configurations with magnetic fields is noticeably larger than the maximum mass attainable by uniform rotation, and that for fixed values of baryon number the maximum mass configurations are all characterized by an off-center density maximum.

We present simulations of the evolution of a proto-neutron star in which kaon-condensed matter might exist, including the effects of finite temperature and trapped neutrinos. The phase transition from pure nucleonic matter to the kaon condensate phase is described using Gibbs' rules for phase equilibrium, which permit the existence of a mixed phase. A general property of neutron stars containing kaon condensates, as well as other forms of strangeness, is that the maximum mass for cold, neutrino-free matter can be less than the maximum mass for matter containing trapped neutrinos or which has a finite entropy. A proto-neutron star formed with a baryon mass exceeding that of the maximum mass of cold, neutrino-free matter is therefore metastable, that is, it will collapse to a black hole at some time during the Kelvin-Helmholtz cooling stage. The effects of kaon condensation on metastable stars are dramatic. In these cases, the neutrino signal from a hypothetical galactic supernova (distance $\sim8.5$ kpc) will stop suddenly, generally at a level above the background in the SuperK and SNO detectors, which have low energy thresholds and backgrounds. This is in contrast to the case of a stable star, for which the signal exponentially decays, eventually disappearing into the background. We find the lifetimes of kaon-condensed metastable stars to be restricted to the range 40--70 s and weakly dependent on the proto-neutron star mass, in sharp contrast to the significantly larger mass dependence and range (1--100 s) of hyperon-rich metastable stars.

We study the equation of state (EOS) of kaon-condensed matter including the effects of temperature and trapped neutrinos. It is found that the order of the phase transition to a kaon-condensed phase, and whether or not Gibbs' rules for phase equilibrium can be satisfied in the case of a first order transition, depend sensitively on the choice of the kaon-nucleon interaction. The main effect of finite temperature, for any value of the lepton fraction, is to mute the effects of a first order transition, so that the thermodynamics becomes similar to that of a second order transition. Above a critical temperature, found to be at least 30--60 MeV depending upon the interaction, the first order transition disappears. The phase boundaries in baryon density versus lepton number and baryon density versus temperature planes are delineated. We find that the thermal effects on the maximum gravitational mass of neutron stars are as important as the effects of trapped neutrinos, in contrast to previously studied cases in which the matter contained only nucleons or in which hyperons and/or quark matter were considered. Kaon-condensed EOSs permit the existence of metastable neutron stars, because the maximum mass of an initially hot, lepton-rich protoneutron star is greater than that of a cold, deleptonized neutron star. The large thermal effects imply that a metastable protoneutron star's collapse to a black hole could occur much later than in previously studied cases that allow metastable configurations.

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.

The structure of neutron stars is considered from theoretical and observational perspectives. We demonstrate an important aspect of neutron star structure: the neutron star radius is primarily determined by the behavior of the pressure of matter in the vicinity of nuclear matter equilibrium density. In the event that extreme softening does not occur at these densities, the radius is virtually independent of the mass and is determined by the magnitude of the pressure. For equations of state with extreme softening, or those that are self-bound, the radius is more sensitive to the mass. Our results show that in the absence of extreme softening, a measurement of the radius of a neutron star more accurate than about 1 km will usefully constrain the equation of state. We also show that the pressure near nuclear matter density is primarily a function of the density dependence of the nuclear symmetry energy, while the nuclear incompressibility and skewness parameters play secondary roles. In addition, we show that the moment of inertia and the binding energy of neutron stars, for a large class of equations of state, are nearly universal functions of the star's compactness. These features can be understood by considering two analytic, yet realistic, solutions of Einstein's equations, due, respectively, to Buchdahl and Tolman. We deduce useful approximations for the fraction of the moment of inertia residing in the crust, which is a function of the stellar compactness and, in addition, the presssure at the core-crust interface.

The equation of state (EOS) of dense matter plays an important role in the supernova phenomenon, the structure of neutron stars, and in the mergers of compact objects (neutron stars and black holes). During the collapse phase of a supernova, the EOS at subnuclear densities controls the collapse rate, the amount of deleptonization and thus the size of the collapsing core and the bounce density. Properties of nuclear matter that are especially crucial are the symmetry energy and the nuclear specific heat. The nuclear incompressibility, and the supernuclear EOS, play supporting roles. In a similar way, although the maximum masses of neutron stars are entirely dependent upon the supernuclear EOS, other important structural aspects are more sensitive to the equation of state at nuclear densities. The radii, moments of inertia, and the relative binding energies of neutron stars are, in particular, sensitive to the behavior of the nuclear symmetry energy. The dependence of the radius of a neutron star on its mass is shown to critically influence the outcome of the compact merger of two neutron stars or a neutron star with a small mass black hole. This latter topic is especially relevant to this volume, since it stems from research prompted by the tutoring of David Schramm a quarter century ago.

We study the effects of very strong magnetic fields on the equation of state (EOS) in multicomponent, interacting matter by developing a covariant description for the inclusion of the anomalous magnetic moments of nucleons. For the description of neutron star matter, we employ a field-theoretical approach which permits the study of several models which differ in their behavior at high density. Effects of Landau quantization in ultra-strong magnetic fields ($B>10^{14}$ Gauss) lead to a reduction in the electron chemical potential and a substantial increase in the proton fraction. We find the generic result for $B>10^{18}$ Gauss that the softening of the EOS caused by Landau quantization is overwhelmed by stiffening due to the incorporation of the anomalous magnetic moments of the nucleons. In addition, the neutrons become completely spin polarized. The inclusion of ultra-strong magnetic fields leads to a dramatic increase in the proton fraction, with consequences for the direct Urca process and neutron star cooling. The magnetization of the matter never appears to become very large, as the value of $|H/B|$ never deviates from unity by more than a few percent. Our findings have implications for the structure of neutron stars in the presence of large frozen-in magnetic fields.

With the aim of constraining the structural properties of neutron stars and the equation of state of dense matter, we study sudden spin-ups, glitches, occurring in the Vela pulsar and in six other pulsars. We present evidence that glitches represent a self-regulating instability for which the star prepares over a waiting time. The angular momentum requirements of glitches in Vela indicate that at least 1.4% of the star's moment of inertia drives these events. If glitches originate in the liquid of the inner crust, Vela's `radiation radius' $R_\infty$ must exceed ~12 km for a mass of 1.4 solar masses. The isolated neutron star RX J18563-3754 is a promising candidate for a definitive radius measurement, and offers to further our understanding of dense matter and the origin of glitches.

With the aim of constraining the structural properties of neutron stars and the equation of state of dense matter, we study sudden spin-ups, glitches, occurring in the Vela pulsar and in six other pulsars. We present evidence that glitches represent a self-regulating instability for which the star prepares over a waiting time. The angular momentum requirements of glitches in Vela indicate that at least 1.4% of the star's moment of inertia drives these events. If glitches originate in the liquid of the inner crust, Vela's `radiation radius' must exceed ~12 km for a mass of 1.4 solar masses. Observational tests of whether other neutron stars obey this constraint will be possible in the near future.

An extensive study of the effects of correlations on both charged and neutral current weak interaction rates in dense matter is performed. Both strong and electromagnetic correlations are considered.The propagation of particle-hole interactions in the medium plays an important role in determining the neutrino mean free paths. The effects due to Pauli-Blocking and density, spin, and isospin correlations in the medium significantly reduce the neutrino cross sections. Due to the lack of experimental information at high density, these correlations are necessarily model dependent. For example, spin correlations in nonrelativistic models are found to lead to larger suppressions of neutrino cross sections compared to those of relativistic models. This is due to the tendency of the nonrelativistic models to develop spin instabilities. Notwithstanding the above caveats, and the differences between nonrelativistic and relativistic approaches such as the spin- and isospin-dependent interactions and the nucleon effective masses, suppressions of order 2--3, relative to the case in which correlations are ignored, are obtained. Neutrino interactions in dense matter are especially important for supernova and early neutron star evolution calculations. The effects of correlations for protoneutron star evolution are calculated. Large effects on the internal thermodynamic properties of protoneutron stars, such as the temperature, are found. These translate into significant early enhancements in the emitted neutrino energies and fluxes, especially after a few seconds. At late times, beyond about 10 seconds, the emitted neutrino fluxes decrease more rapidly compared to simulations without the effects of correlations, due to the more rapid onset of neutrino transparency in the protoneutron star

We study the thermal and chemical evolution during the Kelvin-Helmholtz phase of the birth of a neutron star, employing neutrino opacities that are consistently calculated with the underlying equation of state (EOS). Expressions for the diffusion coefficients appropriate for general relativistic neutrino transport in the equilibrium diffusion approximation are derived. The diffusion coefficients are evaluated using a field-theoretical finite temperature EOS that includes the possible presence of hyperons. The variation of the diffusion coefficients is studied as a function of EOS and compositional parameters. We present results from numerical simulations of protoneutron star cooling for internal stellar properties as well as emitted neutrino energies and luminosities. We discuss the influence of the initial stellar model, the total mass, the underlying EOS, and the addition of hyperons on the evolution of the protoneutron star and upon the expected signal in terrestrial detectors.

The early evolution of a protoneutron star depends both on the equation of state and neutrino interactions at high density. We identify the important sources of neutrino opacity and through model calculations show that in-medium effects play an important role in determining the neutrino mean free paths. The effects due to Pauli-blocking and many-body correlations due to strong interactions reduce the neutrino cross sections by large factors compared to the case in which these effects are ignored. We discuss these results in the context of neutrino transport in a protoneutron star.

We study the diffusive transport of neutrinos in a newly born neutron star to explore its sensitivity to dense matter properties. Energy and lepton number which are trapped during the catastrophic implosion diffuse out on the time scale of a few tens of seconds. Results for different dense matter models are presented.

We study the charged and neutral current weak interaction rates relevant for the determination of neutrino opacities in dense matter found in supernovae and neutron stars. We establish an efficient formalism for calculating differential cross sections and mean free paths for interacting, asymmetric nuclear matter at arbitrary degeneracy. The formalism is valid for both charged and neutral current reactions. Strong interaction corrections are incorporated through the in-medium single particle energies at the relevant density and temperature. The effects of strong interactions on the weak interaction rates are investigated using both potential and effective field-theoretical models of matter. We investigate the relative importance of charged and neutral currents for different astrophysical situations, and also examine the influence of strangeness-bearing hyperons. Our findings show that the mean free paths are significantly altered by the effects of strong interactions and the multi-component nature of dense matter. The opacities are then discussed in the context of the evolution of the core of a protoneutron star.

A protoneutron star is formed immediately after the gravitational collapse of the core of a massive star. At birth, the hot and high density matter in such a star contains a large number of neutrinos trapped during collapse. Trapped neutrinos generally inhibit the presence of exotic matter -- hyperons, a kaon condensate, or quarks. However, as the neutrinos diffuse out in about 10-15 s, the threshold for the appearance of strangeness is reduced; hence, the composition and the structure of the star can change significantly. The effect of exotic, negatively-charged, strangeness-bearing components is always to soften the equation of state, and the possibility exists that the star collapses to a black hole at this time. This could explain why no neutron star has yet been seen in the remnant of supernova SN1987A, even though one certainly existed when neutrinos were detected on Feb. 23, 1987. With new generation neutrino detectors it is feasible to test different theoretical scenarios observationally.

Immediately after they are born, neutron stars are characterized by an entropy per baryon of order unity and by the presence of trapped neutrinos. If the only hadrons in the star are nucleons, these effects slightly reduce the maximum mass relative to cold, catalyzed matter. However, if negatively charged particles in the form of hyperons, a kaon condensate, or quarks are also present, these effects result in an increase in the maximum mass of $\sim 0.2{\rm M}_{\odot}$ compared to that of a cold, neutrino-free star. This could lead to the delayed formation of a black hole; such a scenario is consistent with our present knowledge of SN1987A.

We approximate boson thermodynamic integrals as polynomials in two variables chosen to give the correct limiting expansions and to smoothly interpolate into other regimes. With 10 free parameters, an accuracy of better than 0.009\% is achieved for the pressure, internal energy density and the number density. We also revisit the fermion case, originally addressed by Eggleton, Faulkner and Flannery (1973), and substantially improve the accuracy of their fits.

We investigate the structure of neutron stars shortly after they are born, when the entropy per baryon is of order 1 or 2 and neutrinos are trapped on dynamical timescales. In all cases, the thermal effects for an entropy per baryon of order 2 or less are small when considering the maximum neutron star mass. Neutrino trapping, however, significantly changes the maximum mass due to the abundance of electrons. When matter is allowed to contain only nucleons and leptons, trapping decreases the maximum mass by an amount comparable to, but somewhat larger than, the increase due to finite entropy. When matter is allowed to contain strongly interacting negatively charged particles, in the form of strange baryons, a kaon condensate, or quarks, trapping instead results in an increase in the maximum mass of order $0.2M_\odot$, which adds to the effects of finite entropy. The presence of negatively-charged particles has two major implications. First, the value of the maximum mass will decrease during the early evolution of a neutron star as it loses trapped neutrinos, so that if a black hole forms, it either does so immediately after the bounce or it is delayed for a neutrino diffusion timescale of $\sim 10$ s. The latter case is most likely if the maximum mass of the hot star with trapped neutrinos is near $1.5M_\odot$. In the absence of negatively-charged hadrons, black hole formation would be due to accretion and therefore is likely to occur only immediately after bounce. Second, the appearance of hadronic negative charges results in a general softening of the equation of state that may be observable in the neutrino luminosities and average energies. Further, these additional negative charges decrease the electron fraction and may be observed in the relative excess of electron neutrinos compared to other neutrinos.

We investigate the possibility of kaon condensation in the dense interior of neutron stars through the s--wave interaction of kaons with nucleons. We include nucleon--nucleon interactions by using simple parametrizations of realistic forces, and include electrons and muons in $\beta$--equilibrium. The conditions under which kaon condensed cores undergo a transition to quark matter containing strange quarks are also established. The critical density for kaon condensation lies in the range (2.3--5.0)$~\rho_0$, where $\rho_0=0.16$~fm$^{-3}$ is the equilibrium density of nuclear matter. For too large a value of the strangeness content of the proton, matter with a kaon condensate is not sufficiently stiff to support the lower limit of $1.44 M_\odot$ for a neutron star. Kaon condensation dramatically increases the proton abundance of matter and even allows positrons to exist inside the core. We also consider the case when neutrinos are trapped in the matter. Neutrino trapping shifts both kaon condensation and the quark matter transition to higher densities than in the case of cold, catalyzed matter. A newly--formed neutron star is expected to have a rather low central density, the density rising only after mass accretion onto the star ends after a few seconds. Thus, it is likely that if kaon condensation and/or the quark--hadron phase transition occur, they do so only during or after the mass accretion and neutrino trapping stages. We suggest that neutrino observations from a galactic supernova may provide direct evidence for or against a condensate and/or a quark--hadron transition.