Source author record

Andrew W. Steiner

Andrew W. Steiner appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2021arXiv

Constraining nuclear matter parameters with GW170817

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.

preprint2021arXiv

Hot and Dense Matter Equation of State Probability Distributions for Astrophysical Simulations

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

preprint2020arXiv

Combining Electromagnetic and Gravitational-Wave Constraints on Neutron-Star Masses and Radii

We perform a joint Bayesian inference of neutron-star mass and radius constraints based on GW170817, observations of quiescent low-mass X-ray binaries (QLMXBs), photospheric radius expansion X-ray bursts (PREs), and X-ray timing observations of J0030+0451. With this data set, the form of the prior distribution still has an impact on the posterior mass-radius (MR) curves and equation of state (EOS), but this impact is smaller than recently obtained when considering QLMXBs alone. We analyze the consistency of the electromagnetic data by including an "intrinsic scattering" contribution to the uncertainties, and find only a slight broadening of the posteriors. This suggests that the gravitational-wave and electromagnetic observations of neutron-star structure are providing a consistent picture of the neutron-star mass-radius curve and the EOS.

preprint2020arXiv

Future Prospects for Constraining Nuclear Matter Parameters with Gravitational Waves

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.

preprint2020arXiv

Measuring Nuclear Matter Parameters with NICER and LIGO/Virgo

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-$σ$ 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.

preprint2016arXiv

Measuring the neutron star equation of state using X-ray timing

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]

preprint2016arXiv

White Paper on Nuclear Astrophysics

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.

preprint2015arXiv

Sound velocity bound and neutron stars

It has been conjectured that the velocity of sound in any medium is smaller than the velocity of light in vacuum divided by $\sqrt{3}$. Simple arguments support this bound in non-relativistic and/or weakly coupled theories. The bound has been demonstrated in several classes of strongly coupled theories with gravity duals and is saturated only in conformal theories. We point out that the existence of neutron stars with masses around two solar masses combined with the knowledge of the equation of state of hadronic matter at "low" densities is in strong tension with this bound.

preprint2015arXiv

The Fate of the Compact Remnant in Neutron Star Mergers

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.

preprint2015arXiv

Using Neutron Star Observations to Determine Crust Thicknesses, Moments of Inertia, and Tidal Deformabilities

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.

preprint2014arXiv

Constraints on the Symmetry Energy Using the Mass-Radius Relation of Neutron Stars

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.

preprint2014arXiv

Hypernuclei and the hyperon problem in neutron stars

The likely presence of $Λ$ baryons in dense hadronic matter tends to soften the equation of state to an extend that the observed heaviest neutron stars are difficult to explain. We analyze this "hyperon problem" with a phenomenological approach. First, we review what can be learned about the interaction of $Λ$ particle with dense matter from the observed hypernuclei and extend this phenomenological analysis to asymmetric matter. We add to this the current knowledge on non-strange dense matter, including its uncertainties, to conclude that the interaction between $Λ$s and dense matter has to become repulsive at densities below three times the nuclear saturation density.

preprint2014arXiv

Magnetar giant flare oscillations and the nuclear symmetry energy

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.

preprint2014arXiv

Moving Beyond Chi-Squared in Nuclei and Neutron Stars

There are several assumptions made in a standard $χ^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.

preprint2014arXiv

Neutron Star Masses and Radii from Quiescent Low-Mass X-ray Binaries

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.

preprint2014arXiv

Stellar Superfluids

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.

preprint2013arXiv

From nuclear structure to neutron stars

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.

preprint2012arXiv

Core-collapse supernova equations of state based on neutron star observations

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.

preprint2012arXiv

Deep Crustal Heating in a Multicomponent Accreted Neutron Star Crust

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.

preprint2012arXiv

The Neutron Star Mass-Radius Relation and the Equation of State of Dense Matter

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.

preprint2011arXiv

Rapid Cooling of the Neutron Star in Cassiopeia A Triggered by Neutron Superfluidity in Dense Matter

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.

preprint2011arXiv

Superfluid Neutrons in the Core of the Neutron Star in Cassiopeia A

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.

preprint2010arXiv

The Equation of State from Observed Masses and Radii of Neutron Stars

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.

preprint2008arXiv

Viscous damping of r-mode oscillations in compact stars with quark matter

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.

preprint2005arXiv

The Color-Superconducting {'}t Hooft Interaction

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.

preprint2004arXiv

Minimal Cooling of Neutron Stars: A New Paradigm

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.