-
Strengthening nuclear symmetry energy constraints using multiple resonant shattering flares of neutron stars with realistic mass uncertainties
Authors:
Duncan Neill,
David Tsang,
William G. Newton
Abstract:
With current and planned gravitational-wave (GW) observing runs, coincident multimessenger timing of Resonant Shattering Flares (RSFs) and GWs may soon allow for neutron star (NS) asteroseismology to be used to constrain the nuclear symmetry energy, an important property of fundamental nuclear physics that influences the composition and equation of state of NSs. In this work we examine the effects…
▽ More
With current and planned gravitational-wave (GW) observing runs, coincident multimessenger timing of Resonant Shattering Flares (RSFs) and GWs may soon allow for neutron star (NS) asteroseismology to be used to constrain the nuclear symmetry energy, an important property of fundamental nuclear physics that influences the composition and equation of state of NSs. In this work we examine the effects of combining multiple RSF detections on these symmetry energy constraints, and consider how realistic uncertainties in the masses of the progenitor NSs may weaken them. We show that the detection of subsequent multimessenger events has the potential to substantially improve constraints beyond those obtained from the first, and that this improvement is insensitive to the mass of the NSs which produce the RSFs and its uncertainty. This sets these asteroseismic constraints apart from bulk NS properties such as radius, for which the NS mass is highly important, meaning that any multimessenger RSF and GW events can equally improve our knowledge of fundamental physics.
△ Less
Submitted 6 March, 2024;
originally announced March 2024.
-
Constraining Nuclear Symmetry Energy with Multi-messenger Resonant Shattering Flares
Authors:
Duncan Neill,
Rebecca Preston,
William G. Newton,
David Tsang
Abstract:
Much effort is devoted to measuring the nuclear symmetry energy through neutron star (NS) and nuclear observables. Since matter in the NS core may be non-hadronic, observables like radii and tidal deformability may not provide reliable constraints on properties of nucleonic matter. We demonstrate that coincident timing of a resonant shattering flare (RSF) and gravitational wave signal during binar…
▽ More
Much effort is devoted to measuring the nuclear symmetry energy through neutron star (NS) and nuclear observables. Since matter in the NS core may be non-hadronic, observables like radii and tidal deformability may not provide reliable constraints on properties of nucleonic matter. We demonstrate that coincident timing of a resonant shattering flare (RSF) and gravitational wave signal during binary NS inspiral probes the crust-core transition region and provides constraints on the symmetry energy comparable to terrestrial nuclear experiments. We show that nuclear masses, RSFs and measurements of NS radii and tidal deformabilities constrain different density ranges of the EOS, providing complementary probes.
△ Less
Submitted 1 August, 2022;
originally announced August 2022.
-
Horizons: Nuclear Astrophysics in the 2020s and Beyond
Authors:
H. Schatz,
A. D. Becerril Reyes,
A. Best,
E. F. Brown,
K. Chatziioannou,
K. A. Chipps,
C. M. Deibel,
R. Ezzeddine,
D. K. Galloway,
C. J. Hansen,
F. Herwig,
A. P. Ji,
M. Lugaro,
Z. Meisel,
D. Norman,
J. S. Read,
L. F. Roberts,
A. Spyrou,
I. Tews,
F. X. Timmes,
C. Travaglio,
N. Vassh,
C. Abia,
P. Adsley,
S. Agarwal
, et al. (140 additional authors not shown)
Abstract:
Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilit…
▽ More
Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilities across an ever growing number of disciplines and subfields that need to be integrated. We take a holistic view of the field discussing the unique challenges and opportunities in nuclear astrophysics in regards to science, diversity, education, and the interdisciplinarity and breadth of the field. Clearly nuclear astrophysics is a dynamic field with a bright future that is entering a new era of discovery opportunities.
△ Less
Submitted 16 May, 2022;
originally announced May 2022.
-
Ensembles of unified crust and core equations of state in a nuclear-multimessenger astrophysics environment
Authors:
William G. Newton,
Lauren Balliet,
Srdan Budimir,
Gabriel Crocombe,
Brianna Douglas,
Thomas Blake Head,
Zach Langford,
Luis Rivera,
Josh Sanford
Abstract:
We present an ensemble of unified neutron star crust and core equations of state, constructed using an extended Skyrme energy density functional through the crust and outer core, and appended by two piecewise polytropes at higher densities. The equations of state are parameterized by the first three coefficients in the density expansion of the symmetry energy $J,L$ and $K_{\rm sym}$, the moment of…
▽ More
We present an ensemble of unified neutron star crust and core equations of state, constructed using an extended Skyrme energy density functional through the crust and outer core, and appended by two piecewise polytropes at higher densities. The equations of state are parameterized by the first three coefficients in the density expansion of the symmetry energy $J,L$ and $K_{\rm sym}$, the moment of inertia of a 1.338 M$_{\odot}$ star $I_{1.338}$ and the maximum neutron star mass $M_{\rm max}$. We construct an ensemble with uniform priors on all five parameters, and then apply data filters to the ensemble to explore the effect of combining neutron skin data from PREX with astrophysical measurements of radii and tidal deformabilities from NICER and LIGO/VIRGO. Neutron skins are calculated directly using the EDFs. We demonstrate that both the nuclear data and astrophysical data play a role in constraining crust properties such as the mass, thickness and moment of inertia of the crust and the nuclear pasta layers therein, and that astrophysical data better constrains $K_{\rm sym}$ than PREX data.
△ Less
Submitted 9 March, 2022; v1 submitted 22 December, 2021;
originally announced December 2021.
-
Resonant Shattering Flares in Black Hole-Neutron Star and Binary Neutron Star Mergers
Authors:
Duncan Neill,
David Tsang,
Hendrik van Eerten,
Geoffrey Ryan,
William G. Newton
Abstract:
Resonant Shattering flares (RSFs) are bursts of gamma-rays expected to be triggered by tidal resonance of a neutron star (NS) during binary inspiral. They are strongly dependent on the magnetic field strength at the surface of the NS. By modelling these flares as being the result of multiple colliding relativistic shells launched during the resonance window, we find that the prompt non-thermal gam…
▽ More
Resonant Shattering flares (RSFs) are bursts of gamma-rays expected to be triggered by tidal resonance of a neutron star (NS) during binary inspiral. They are strongly dependent on the magnetic field strength at the surface of the NS. By modelling these flares as being the result of multiple colliding relativistic shells launched during the resonance window, we find that the prompt non-thermal gamma-ray emission may have luminosity up to a few $\times10^{48}\text{ erg/s}$, and that a broad-band afterglow could be produced. We compute the expected rates of detectable RSFs using the BPASS population synthesis code, with different assumptions about the evolution of surface magnetic field strengths before merger. We find the rate of detectable RSFs to be $\sim 0.0001-5$ per year for BHNS mergers and $\sim 0.0005-25$ per year for NSNS mergers, with the lower bound corresponding to surface-field decay consistent with magneto-thermal evolution in purely crustal fields, while the upper bounds are for systems which have longer-lived surface magnetic fields supported by flux frozen into the superconducting core. If some of the observed SGRB precursor flares are indeed RSFs, this suggests the presence of a longer-lived surface field for some fraction of the neutron star population, and that we could expect RSFs to be the most common detectable EM counterpart to GW detections of BHNS mergers. The non-detection of a RSF prior to GRB170817A provides an upper bound on the magnetic fields of the progenitor NSs of $B_{\rm surf}\sim 10^{13.5} \text{ G}$.
△ Less
Submitted 10 June, 2022; v1 submitted 5 November, 2021;
originally announced November 2021.
-
Glassy quantum nuclear pasta in neutron star crusts
Authors:
William G. Newton,
Mark Alexander Kaltenborn,
Sarah Cantu,
Shuxi Wang,
Amber Stinson,
Jirina Rikovska Stone
Abstract:
We investigate the nuclear pasta phases in neutron star crusts by conducting a large number of three-dimensional Hartree-Fock+BCS calculations at densities leading to the crust-core transition. We survey the shape parameter space of pasta at constant pressure. Spaghetti, waffles, lasagna, bi-continuous phases and cylindrical holes occupy local minima in the resulting Gibbs energy surfaces. The bi-…
▽ More
We investigate the nuclear pasta phases in neutron star crusts by conducting a large number of three-dimensional Hartree-Fock+BCS calculations at densities leading to the crust-core transition. We survey the shape parameter space of pasta at constant pressure. Spaghetti, waffles, lasagna, bi-continuous phases and cylindrical holes occupy local minima in the resulting Gibbs energy surfaces. The bi-continuous phase, in which both the neutron gas and nuclear matter extend continuously in all dimensions and therefore protons are delocalized, appears over a large range of depths. Our results support the idea that nuclear pasta is a glassy system. Multiple pasta configurations coexist in a given layer of the crust. At a characteristic temperature, of order $10^8$-$10^9$K, different phases become frozen into domains whose sizes we estimate to be 1-50 times the lattice spacing and over which the local density and electron fraction can vary. Above this temperature, there is very little long-range order and matter is an amorphous solid. Electron scattering off domain boundaries may contribute to the disorder resistivity of the pasta phases. Annealing of the domains may occur during cooling; repopulating of local minima during crustal heating might lead to temperature dependent transport properties in the deep layers of the crust. We identify 4 distinct regions: (1) nuclear pasta first appears as a local minima, but spherical nuclei are the ground state; (2) nuclear pasta become the absolute minimum, but spherical nuclei are still a local minimum (3) only nuclear pasta appears in local minima, and protons are still localized in at least one dimension (4) only pasta appears, and protons are delocalized. The whole pasta region can occupy up to 70% of the crust by mass and 40% by thickness, and the layer in which protons are delocalized could occupy 45% of the crust mass and 25% of its thickness.
△ Less
Submitted 23 April, 2021;
originally announced April 2021.
-
Resonant Shattering Flares as Multimessenger Probes of the Nuclear Symmetry Energy
Authors:
Duncan Neill,
William G. Newton,
David Tsang
Abstract:
The behaviour of the nuclear symmetry energy near saturation density is important for our understanding of dense nuclear matter. This density dependence can be parameterised by the nuclear symmetry energy and its derivatives evaluated at nuclear saturation density. In this work we show that the core-crust interface mode of a neutron star is sensitive to these parameters, through the (density-weigh…
▽ More
The behaviour of the nuclear symmetry energy near saturation density is important for our understanding of dense nuclear matter. This density dependence can be parameterised by the nuclear symmetry energy and its derivatives evaluated at nuclear saturation density. In this work we show that the core-crust interface mode of a neutron star is sensitive to these parameters, through the (density-weighted) shear-speed within the crust, which is in turn dependent on the symmetry energy profile of dense matter. We calculate the frequency at which the neutron star quadrupole ($\ell = 2$) crust-core interface mode must be driven by the tidal field of its binary partner to trigger a Resonant Shattering Flare (RSF). We demonstrate that coincident multimessenger timing of an RSF and gravitational wave chirp from a neutron star merger would enable us to place constraints on the symmetry energy parameters that are competitive with those from current nuclear experiments.
△ Less
Submitted 11 March, 2021; v1 submitted 18 December, 2020;
originally announced December 2020.
-
Superconducting phases in a two-component microscale model of neutron star cores
Authors:
Toby S. Wood,
Vanessa Graber,
William G. Newton
Abstract:
We identify the possible ground states for a mixture of two superfluid condensates (one neutral, the other electrically charged) using a phenomenological Ginzburg-Landau model. While this framework is applicable to any interacting condensed-matter mixture of a charged and a neutral component, we focus on nuclear matter in neutron star cores, where proton and neutron condensates are coupled via non…
▽ More
We identify the possible ground states for a mixture of two superfluid condensates (one neutral, the other electrically charged) using a phenomenological Ginzburg-Landau model. While this framework is applicable to any interacting condensed-matter mixture of a charged and a neutral component, we focus on nuclear matter in neutron star cores, where proton and neutron condensates are coupled via non-dissipative entrainment. We employ the Skyrme interaction to determine the neutron star's equilibrium composition, and hence obtain realistic coefficients for our Ginzburg-Landau model at each depth within the star's core. We then use the Ginzburg-Landau model to determine the ground state in the presence of a magnetic field. In this way, we obtain superconducting phase diagrams for six representative Skyrme models, revealing the microphysical magnetic flux distribution throughout the neutron star core. The phase diagrams are rather complex and the locations of most of the phase transitions can only be determined through numerical calculations. Nonetheless, we find that for all equations of state considered in this work, much of the outer core exhibits type-1.5 superconductivity, rather than type-II superconductivity as is generally assumed. For local magnetic field strengths $\lesssim 10^{14} \, {\rm G}$, the magnetic flux is distributed inhomogeneously, with bundles of magnetic fluxtubes separated by flux-free Meissner regions. We provide an approximate criterion to determine the transition between this type-1.5 phase and the type-I region in the inner core.
△ Less
Submitted 26 February, 2021; v1 submitted 5 November, 2020;
originally announced November 2020.
-
The nuclear symmetry energy from neutron skins and pure neutron matter in a Bayesian framework
Authors:
William G. Newton,
Gabriel Crocombe
Abstract:
We present an inference of the nuclear symmetry energy magnitude $J$, the slope $L$ and the curvature $K_{\rm sym}$ by combining neutron skin data on Ca, Pb and Sn isotopes and our best theoretical information about pure neutron matter (PNM). A Bayesian framework is used to consistently incorporate prior knowledge of the PNM equation of state from chiral effective field theory calculations. Neutro…
▽ More
We present an inference of the nuclear symmetry energy magnitude $J$, the slope $L$ and the curvature $K_{\rm sym}$ by combining neutron skin data on Ca, Pb and Sn isotopes and our best theoretical information about pure neutron matter (PNM). A Bayesian framework is used to consistently incorporate prior knowledge of the PNM equation of state from chiral effective field theory calculations. Neutron skins are modeled in a Hartree-Fock approach using an extended Skyrme energy-density functional which allows for independent variation of $J$, $L$ and $K_{\rm sym}$ without affecting the symmetric nuclear matter equation of state. We discuss the choice of neutron skin data sets, and combining errors in quadrature we obtain 95\% credible values of $J=31.3\substack{+4.2 \\ -5.9}$ MeV, $L=40\substack{+34 \\ -26}$ MeV and $K_τ = L - 6K_{\rm sym}= -444\substack{+100 \\ -84}$ MeV using uninformative priors in $J$, $L$ and $K_{\rm sym}$, and $J=31.9\substack{+1.3 \\ -1.3}$ MeV, $L=37\substack{+9 \\ -8}$ MeV and $K_τ = -480\substack{+25 \\ -26}$ MeV using PNM priors. The correlations between symmetry energy parameters induced by neutron skin data is discussed and compared with the droplet model. Neutron skin data alone is shown to place limits on the symmetry energy parameters as stringent as those obtained from chiral effective field theory alone, and when combined the 95\% credible intervals are reduced by a factor of 4-5. Ahead of new measurements of lead and calcium neutron skins from parity-violating electron scattering experiments at Jefferson Lab and Mainz Superconducting Accelerator, we make predictions based on existing data on neutron skins of tin for the neutron skins of calcium and lead of 0.166$\pm$0.008 fm and $0.169 \pm 0.014$ fm respectively, using uninformative priors, and 0.167$\pm$0.008 fm and $0.172 \pm 0.015$ fm respectively, using PNM priors.
△ Less
Submitted 23 April, 2021; v1 submitted 31 July, 2020;
originally announced August 2020.
-
The state of matter in simulations of core-collapse supernovae -- Reflections and recent developments
Authors:
Tobias Fischer,
Niels-Uwe Bastian,
David Blaschke,
Mateusz Cierniak,
Matthias Hempel,
Thomas Klähn,
Gabriel Martínez-Pinedo,
William G. Newton,
Gerd Röpke,
Stefan Typel
Abstract:
In this review article we discuss selected developments regarding the role of the equation of state (EOS) in simulations of core-collapse supernovae. There are no first-principle calculations of the state of matter under supernova conditions since a wide range of conditions is covered, in terms of density, temperature and isospin asymmetry. Instead, model EOS are commonly employed in supernova stu…
▽ More
In this review article we discuss selected developments regarding the role of the equation of state (EOS) in simulations of core-collapse supernovae. There are no first-principle calculations of the state of matter under supernova conditions since a wide range of conditions is covered, in terms of density, temperature and isospin asymmetry. Instead, model EOS are commonly employed in supernova studies. These can be divided into regimes with intrinsically different degrees of freedom: heavy nuclei at low temperatures, inhomogeneous nuclear matter where light and heavy nuclei coexist together with unbound nucleons, and the transition to homogeneous matter at high densities and temperatures. In this article we discuss each of these phases with particular view on their role in supernova simulations.
△ Less
Submitted 8 February, 2018; v1 submitted 20 November, 2017;
originally announced November 2017.
-
Nuclear pasta and supernova neutrinos at late times
Authors:
C. J. Horowitz,
D. K. Berry,
M. E. Caplan,
T. Fischer,
Zidu Lin,
W. G. Newton,
E. O'Connor,
L. F. Roberts
Abstract:
Nuclear pasta, with nucleons arranged into tubes, sheets, or other complex shapes, is expected in core collapse supernovae (SNe) at just below nuclear density. We calculate the additional opacity from neutrino-pasta coherent scattering using molecular dynamics simulations. We approximately include this opacity in simulations of SNe. We find that pasta slows neutrino diffusion and greatly increases…
▽ More
Nuclear pasta, with nucleons arranged into tubes, sheets, or other complex shapes, is expected in core collapse supernovae (SNe) at just below nuclear density. We calculate the additional opacity from neutrino-pasta coherent scattering using molecular dynamics simulations. We approximately include this opacity in simulations of SNe. We find that pasta slows neutrino diffusion and greatly increases the neutrino signal at late times of 10 or more seconds after stellar core collapse. This signal, for a galactic SN, should be clearly visible in large detectors such as Super-Kamiokande.
△ Less
Submitted 30 November, 2016;
originally announced November 2016.
-
Testing the formation scenarios of binary neutron star systems with measurements of the neutron star moment of inertia
Authors:
William G. Newton,
Andrew W. Steiner,
Kent Yagi
Abstract:
Two low mass neutron stars, J0737-3039B and the companion to J1756-2251, show strong evidence of being formed from the collapse of an ONeMg core in an electron capture supernova (ECSN) or in an ultra-stripped iron core collapse supernova (FeCCSN). Using three different systematically generated sets of equations of state we explore the relationship between the moment of inertia of J0737-3039A and t…
▽ More
Two low mass neutron stars, J0737-3039B and the companion to J1756-2251, show strong evidence of being formed from the collapse of an ONeMg core in an electron capture supernova (ECSN) or in an ultra-stripped iron core collapse supernova (FeCCSN). Using three different systematically generated sets of equations of state we explore the relationship between the moment of inertia of J0737-3039A and the binding energy of the two low mass neutron stars. We find this relationship, a less strict variant of the recently discovered I-Love-Q relations, is nevertheless more robust than a previously explored correlation between the binding energy and the slope of the nuclear symmetry energy L. We find that, if either J0737-3039B or the J1756-2251 companion were formed in an ECSN, no more than 0.06 solar masses could have been lost from the progenitor core, more than four times the mass loss predicted by current supernova modeling. A measurement of the moment of inertia of J0737-3039A to within 10% accuracy from pulsar timing, possible within a decade, can discriminate between formation scenarios such as ECSN or ultra-stripped FeCCSN and, given current constraints on the predicted core mass loss, potentially rule them out. Using the I-Love-Q relations we find that an Advanced LIGO can potentially measure the neutron star tidal polarizability to equivalent accuracy in a neutron star-neutron star merger at a distance of 200 Mpc, thus obtaining similar constraints on the formation scenarios. Such information on the occurrence of ECSNe is important for population synthesis calculations, especially for estimating the rate of binary neutron star mergers and resulting electromagnetic and gravitational wave signals. Further progress needs to be made modeling the core collapse process that leads to low-mass neutron stars, particularly in making robust predictions for the mass loss from the progenitor core.
△ Less
Submitted 28 November, 2016;
originally announced November 2016.
-
Shedding Light on the EOS-Gravity Degeneracy and Constraining the Nuclear Symmetry Energy from the Gravitational Binding Energy of Neutron Stars
Authors:
Xiao-Tao He,
F. J. Fattoyev,
Bao-An Li,
W. G. Newton
Abstract:
A thorough understanding of properties of neutron stars requires both a reliable knowledge of the equation of state (EOS) of super-dense nuclear matter and the strong-field gravity theories simultaneously. To provide information that may help break this EOS-gravity degeneracy, we investigate effects of nuclear symmetry energy on the gravitational binding energy of neutron stars within GR and the s…
▽ More
A thorough understanding of properties of neutron stars requires both a reliable knowledge of the equation of state (EOS) of super-dense nuclear matter and the strong-field gravity theories simultaneously. To provide information that may help break this EOS-gravity degeneracy, we investigate effects of nuclear symmetry energy on the gravitational binding energy of neutron stars within GR and the scalar-tensor subset of alternative gravity models. We focus on effects of the slope $L$ of nuclear symmetry energy at saturation density and the high-density behavior of nuclear symmetry energy. We find that the variation of either the density slope $L$ or the high-density behavior of nuclear symmetry energy leads to large changes in the binding energy of neutron stars. The difference in predictions using the GR and the scalar-tensor theory appears only for massive neutron stars, and even then is significantly smaller than the difference resulting from variations in the symmetry energy.
△ Less
Submitted 14 October, 2015;
originally announced October 2015.
-
Constraints on the symmetry energy from observational probes of the neutron star crust
Authors:
William G. Newton,
Joshua Hooker,
Michael Gearheart,
Kyleah Murphy,
De-Hua Wen,
Farrukh Fattoyev,
Bao-An Li
Abstract:
A number of observed phenomena associated with individual neutron star systems or neutron star populations find explanations in models in which the neutron star crust plays an important role. We review recent work examining the sensitivity to the slope of the symmetry energy $L$ of such models, and constraints extracted on $L$ from confronting them with observations. We focus on six sets of observ…
▽ More
A number of observed phenomena associated with individual neutron star systems or neutron star populations find explanations in models in which the neutron star crust plays an important role. We review recent work examining the sensitivity to the slope of the symmetry energy $L$ of such models, and constraints extracted on $L$ from confronting them with observations. We focus on six sets of observations and proposed explanations: (i) The cooling rate of the neutron star in Cassiopeia A, confronting cooling models which include enhanced cooling in the nuclear pasta regions of the inner crust, (ii) the upper limit of the observed periods of young X-ray pulsars, confronting models of magnetic field decay in the crust caused by the high resistivity of the nuclear pasta layer, (iii) glitches from the Vela pulsar, confronting the paradigm that they arise due to a sudden re-coupling of the crustal neutron superfluid to the crustal lattice after a period during which they were decoupled due to vortex pinning, (iv) The frequencies of quasi-periodic oscillations in the X-ray tail of light curves from giant flares from soft gamma-ray repeaters, confronting models of torsional crust oscillations, (v) the upper limit on the frequency to which millisecond pulsars can be spun-up due to accretion from a binary companion, confronting models of the r-mode instability arising above a threshold frequency determined in part by the viscous dissipation timescale at the crust-core boundary, and (vi) the observations of precursor electromagnetic flares a few seconds before short gamma-ray bursts, confronting a model of crust shattering caused by resonant excitation of a crustal oscillation mode by the tidal gravitational field of a companion neutron star just before merger.
△ Less
Submitted 6 June, 2015;
originally announced June 2015.
-
Observational constraints on neutron star crust-core coupling during glitches
Authors:
William G. Newton,
Sabrina Berger,
Brynmor Haskell
Abstract:
We demonstrate that observations of glitches in the Vela pulsar can be used to investigate the strength of the crust-core coupling in a neutron star, and suggest that recovery from the glitch is dominated by torque exerted by the re-coupling of superfluid components of the core that were decoupled from the crust during the glitch. Assuming that the recoupling is mediated by mutual friction between…
▽ More
We demonstrate that observations of glitches in the Vela pulsar can be used to investigate the strength of the crust-core coupling in a neutron star, and suggest that recovery from the glitch is dominated by torque exerted by the re-coupling of superfluid components of the core that were decoupled from the crust during the glitch. Assuming that the recoupling is mediated by mutual friction between the superfluid neutrons and the charged components of the core, we use the observed magnitudes and timescales of the shortest timescale components of the recoveries from two recent glitches in the Vela pulsar to infer the fraction of the core that is coupled to the crust during the glitch, and hence spun up by the glitch event. Within the framework of a two-fluid hydrodynamic model of glitches, we analyze whether crustal neutrons alone are sufficient to drive the glitch activity observed in the Vela pulsar. We use two sets of neutron star equations of state (EOSs), both of which span crust and core consistently and cover a range of the slope of the symmetry energy at saturation density $30 < L <120$ MeV. One set produces maximum masses $\approx$2.0$M_{\odot}$, the second $\approx$2.6$M_{\odot}$. We also include the effects of entrainment of crustal neutrons by the superfluid lattice. We find that for medium to stiff EOSs, observations imply $>70\%$ of the moment of inertia of the core is coupled to the crust during the glitch, though for softer EOSs $L\approx 30$MeV as little as $5\%$ could be coupled. No EOS is able to reproduce the observed glitch activity with crust neutrons alone, but extending the region where superfluid vortices are strongly pinned into the core by densities as little as 0.016fm$^{-3}$ above the crust-core transition density restores agreement with the observed glitch activity.
△ Less
Submitted 3 June, 2015;
originally announced June 2015.
-
Critical Density and Impact of $Δ(1232)$ Resonance Formation in Neutron Stars
Authors:
Bao-Jun Cai,
Farrukh J. Fattoyev,
Bao-An Li,
William G. Newton
Abstract:
The critical densities and impact of forming \D resonances in neutron stars are investigated within an extended nonlinear relativistic mean-field (RMF) model. The critical densities for the formation of four different charge states of \D are found to depend differently on the separate kinetic and potential parts of nuclear symmetry energy, the first example of a microphysical property of neutron s…
▽ More
The critical densities and impact of forming \D resonances in neutron stars are investigated within an extended nonlinear relativistic mean-field (RMF) model. The critical densities for the formation of four different charge states of \D are found to depend differently on the separate kinetic and potential parts of nuclear symmetry energy, the first example of a microphysical property of neutron stars to do so. Moreover, they are sensitive to the in-medium Delta mass $m_Δ$ and the completely unknown $Δ$-$ρ$ coupling strength $g_{ρΔ}$. In the universal baryon-meson coupling scheme where the respective $Δ$-meson and nucleon-meson coupling constants are assumed to be the same, the critical density for the first $Δ^-(1232)$ to appear is found to be \rc=$(2.08\pm0.02)ρ_0$ using RMF model parameters consistent with current constraints on all seven macroscopic parameters usually used to characterize the equation of state (EoS) of isospin-asymmetric nuclear matter (ANM) at saturation density $ρ_0$. Moreover, the composition and the mass-radius relation of neutron stars are found to depend significantly on the values of the $g_{ρΔ}$ and $m_Δ$.
△ Less
Submitted 22 June, 2015; v1 submitted 7 January, 2015;
originally announced January 2015.
-
Impact of the equation-of-state -- gravity degeneracy on constraining the nuclear symmetry energy from astrophysical observables
Authors:
Xiao-Tao He,
F. J. Fattoyev,
Bao-An Li,
W. G. Newton
Abstract:
There is a degeneracy between the equation of state (EOS) of superdense neutron-rich nuclear matter and the strong-field gravity in understanding properties of neutron stars. While the EOS is still poorly known, there are also longstanding ambiguities in choosing General Relativity or alternative gravity theories in the not-so-well tested strong-field regime. Besides possible appearance of hyperon…
▽ More
There is a degeneracy between the equation of state (EOS) of superdense neutron-rich nuclear matter and the strong-field gravity in understanding properties of neutron stars. While the EOS is still poorly known, there are also longstanding ambiguities in choosing General Relativity or alternative gravity theories in the not-so-well tested strong-field regime. Besides possible appearance of hyperons and new phases, the most uncertain part of the nucleonic EOS is currently the density dependence of nuclear symmetry energy. To provide information that may help break the EOS-gravity degeneracy, we investigate effects of symmetry energy within its uncertain range determined by terrestrial nuclear laboratory experiments on the gravitational binding energy and spacetime curvature of neutron stars within GR and the scalar-tensor (ST) theory of gravity. In particular, we focus on effects of the following parameters characterizing the EOS of neutron-rich nucleonic matter: (1) the incompressibility $K_0$ of symmetric nuclear matter, (2) the slope $L$ of symmetry energy at saturation density and (3) the high-density behavior of symmetry energy. We find that the variation of either the density slope $L$ or the high-density behavior of symmetry energy leads to large changes in both the binding energy and curvature of neutron stars while effects of varying the more constrained $K_0$ are negligibly small. The difference in predictions using the GR and the ST theory appears only for massive neutron stars, and is significantly smaller than the differences resulting from variations in the symmetry energy. We conclude that within the ST theory of gravity, the EOS-gravity degeneracy has been broken by the recent relativistic pulsar measurements, and that measurements of neutron star properties sensitive to the compactness constrain mainly the density dependence of the symmetry energy.
△ Less
Submitted 15 December, 2014; v1 submitted 4 August, 2014;
originally announced August 2014.
-
Quantifying Correlations Between Isovector Observables and the Density Dependence of Nuclear Symmetry Energy away from Saturation Density
Authors:
F. J. Fattoyev,
W. G. Newton,
Bao-An Li
Abstract:
According to the Hugenholtz-Van Hove theorem, the nuclear symmetry energy $S(ρ)$ and its slope $L(ρ)$ at arbitrary densities can be decomposed in terms of the density and momentum dependence of the single-nucleon potentials in isospin-asymmetric nuclear matter which are potentially accessible to experiment. We quantify the correlations between several well-known isovector observables and $L(ρ)$ to…
▽ More
According to the Hugenholtz-Van Hove theorem, the nuclear symmetry energy $S(ρ)$ and its slope $L(ρ)$ at arbitrary densities can be decomposed in terms of the density and momentum dependence of the single-nucleon potentials in isospin-asymmetric nuclear matter which are potentially accessible to experiment. We quantify the correlations between several well-known isovector observables and $L(ρ)$ to locate the density range in which each isovector observable is most sensitive to the density dependence of the $S(ρ)$. We then study the correlation coefficients between those isovector observables and all the components of the $L(ρ)$. The neutron skin thickness of $^{208}$Pb is found to be strongly correlated with the $L(ρ)$ at a subsaturation density of $ρ= 0.59 ρ_0$ through the density dependence of the first-order symmetry potential. Neutron star radii are found to be strongly correlated with the $L(ρ)$ over a wide range of supra-saturation densities mainly through both the density and momentum dependence of the first-order symmetry potential. Finally, we find that although the crust-core transition pressure has a complex correlation with the $L(ρ)$, it is strongly correlated with the momentum derivative of the first-order symmetry potential, and the density dependence of the second-order symmetry potential.
△ Less
Submitted 4 May, 2014;
originally announced May 2014.
-
Using Neutron Star Observations to Determine Crust Thicknesses, Moments of Inertia, and Tidal Deformabilities
Authors:
Andrew W. Steiner,
Stefano Gandolfi,
Farrukh J. Fattoyev,
William G. Newton
Abstract:
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…
▽ More
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.
△ Less
Submitted 24 January, 2015; v1 submitted 28 March, 2014;
originally announced March 2014.
-
Probing the high-density behavior of symmetry energy with gravitational waves
Authors:
F. J. Fattoyev,
W. G. Newton,
Bao-An Li
Abstract:
Gravitational wave (GW) astronomy opens up an entirely new window on the Universe to probe the equations of state (EOS) of neutron-rich matter. With the advent of next generation GW detectors, measuring the gravitational radiation from coalescing binary neutron star systems, mountains on rotating neutron stars, and stellar oscillation modes may become possible in the near future. Using a set of mo…
▽ More
Gravitational wave (GW) astronomy opens up an entirely new window on the Universe to probe the equations of state (EOS) of neutron-rich matter. With the advent of next generation GW detectors, measuring the gravitational radiation from coalescing binary neutron star systems, mountains on rotating neutron stars, and stellar oscillation modes may become possible in the near future. Using a set of model EOSs satisfying the latest constraints from terrestrial nuclear experiments, state of the art nuclear many-body calculations of the pure neutron matter EOS, and astrophysical observations consistently, we study various GW signatures of the high-density behavior of the nuclear symmetry energy, which is considered among the most uncertain properties of dense neutron-rich nucleonic matter. In particular, we find the tidal polarizability of neutron stars, potentially measurable in binary systems just prior to merger, is more sensitive to the high density component of the nuclear symmetry energy than the symmetry energy at nuclear saturation density. We also find that the upper limit on the GW strain amplitude from elliptically deformed stars is very sensitive to the density dependence of the symmetry energy. This suggests that future developments in modeling of the neutron star crust, and direct gravitational wave signals from accreting binaries will provide a wealth of information on the EOS of neutron-rich matter. We also review the sensitivity of the $r$-mode instability window to the density dependence of the symmetry energy. Whereas models with larger values of the density slope of the symmetry energy at saturation seem to be disfavored by the current observational data, within a simple $r$-mode model, we point out that a subsequent softer behavior of the symmetry energy at high densities (hinted at by recent observational interpretations) could rule them in.
△ Less
Submitted 19 September, 2013;
originally announced September 2013.
-
The cooling of the Cassiopeia A neutron star as a probe of the nuclear symmetry energy and nuclear pasta
Authors:
William G. Newton,
Kyleah Murphy,
Joshua Hooker,
Bao-An Li
Abstract:
X-ray observations of the neutron star in the Cas A supernova remnant over the past decade suggest the star is undergoing a rapid drop in surface temperature of $\approx$ $2-5.5\%$. One explanation suggests the rapid cooling is triggered by the onset of neutron superfluidity in the core of the star, causing enhanced neutrino emission from neutron Cooper pair breaking and formation (PBF). Using con…
▽ More
X-ray observations of the neutron star in the Cas A supernova remnant over the past decade suggest the star is undergoing a rapid drop in surface temperature of $\approx$ $2-5.5\%$. One explanation suggests the rapid cooling is triggered by the onset of neutron superfluidity in the core of the star, causing enhanced neutrino emission from neutron Cooper pair breaking and formation (PBF). Using consistent neutron star crust and core equations of state (EOSs) and compositions, we explore the sensitivity of this interpretation to the density dependence of the symmetry energy $L$ of the EOS used, and to the presence of enhanced neutrino cooling in the bubble phases of crustal "nuclear pasta". Modeling cooling over a conservative range of neutron star masses and envelope compositions, we find $L\lesssim70$ MeV, competitive with terrestrial experimental constraints and other astrophysical observations. For masses near the most likely mass of $M\gtrsim 1.65 M_{\odot}$, the constraint becomes more restrictive $35\lesssim L\lesssim 55$ MeV. The inclusion of the bubble cooling processes decreases the cooling rate of the star during the PBF phase, matching the observed rate only when $L\lesssim45$ MeV, taking all masses into consideration, corresponding to neutron star radii $\lesssim 11$km.
△ Less
Submitted 22 November, 2013; v1 submitted 9 August, 2013;
originally announced August 2013.
-
Efficacy of crustal superfluid neutrons in pulsar glitch models
Authors:
J. Hooker,
W. G. Newton,
Bao-An Li
Abstract:
In order to assess the ability of purely crust-driven glitch models to match the observed glitch activity in the Vela pulsar, we conduct a systematic analysis of the dependence of the fractional moment of inertia of the inner crustal neutrons on the stiffness of the nuclear symmetry energy at saturation density $L$. We take into account both crustal entrainment and the fact that only a fraction…
▽ More
In order to assess the ability of purely crust-driven glitch models to match the observed glitch activity in the Vela pulsar, we conduct a systematic analysis of the dependence of the fractional moment of inertia of the inner crustal neutrons on the stiffness of the nuclear symmetry energy at saturation density $L$. We take into account both crustal entrainment and the fact that only a fraction $Y_{\rm g}$ of the core neutrons may couple to the crust on the glitch-rise timescale. We use a set of consistently-generated crust and core compositions and equations-of-state which are fit to results of low-density pure neutron matter calculations. When entrainment is included at the level suggested by recent microscopic calculations and the core is fully coupled to the crust, the model is only able to account for the Vela glitch activity for a 1.4$M_{\odot}$ star if the equation of state is particularly stiff $L>100$ MeV. However, an uncertainty of about 10\% in the crust-core transition density and pressure allows for the Vela glitch activity to be marginally accounted for in the range $L\approx30-60$MeV consistent with a range of experimental results. Alternatively, only a small amount of core neutrons need be involved. If less than 50\% of the core neutrons are coupled to the crust during the glitch, we can also account for the Vela glitch activity using crustal neutrons alone for EOSs consistent with the inferred range of $L$. We also explore the possibility of Vela being a high-mass neutron star, and of crustal entrainment being reduced or enhanced relative to its currently predicted values.
△ Less
Submitted 6 June, 2015; v1 submitted 31 July, 2013;
originally announced August 2013.
-
Constraints on the symmetry energy from neutron star observations
Authors:
W. G. Newton,
M. Gearheart,
De-Hua Wen,
Bao-An Li
Abstract:
The modeling of many neutron star observables incorporates the microphysics of both the stellar crust and core, which is tied intimately to the properties of the nuclear matter equation of state (EoS). We explore the predictions of such models over the range of experimentally constrained nuclear matter parameters, focusing on the slope of the symmetry energy at nuclear saturation density $L$. We u…
▽ More
The modeling of many neutron star observables incorporates the microphysics of both the stellar crust and core, which is tied intimately to the properties of the nuclear matter equation of state (EoS). We explore the predictions of such models over the range of experimentally constrained nuclear matter parameters, focusing on the slope of the symmetry energy at nuclear saturation density $L$. We use a consistent model of the composition and EoS of neutron star crust and core matter to model the binding energy of pulsar B of the double pulsar system J0737-3039, the frequencies of torsional oscillations of the neutron star crust and the instability region for r-modes in the neutron star core damped by electron-electron viscosity at the crust-core interface. By confronting these models with observations, we illustrate the potential of astrophysical observables to offer constraints on poorly known nuclear matter parameters complementary to terrestrial experiments, and demonstrate that our models consistently predict $L<70$ MeV.
△ Less
Submitted 18 December, 2012;
originally announced December 2012.
-
Probing Nuclear Symmetry Energy and its Imprints on Properties of Nuclei, Nuclear Reactions, Neutron Stars and Gravitational Waves
Authors:
Bao-An Li,
Lie-Wen Chen,
Farrukh J. Fattoyev,
William G. Newton,
Chang Xu
Abstract:
Significant progress has been made in recent years in constraining nuclear symmetry energy at and below the saturation density of nuclear matter using data from both terrestrial nuclear experiments and astrophysical observations. However, many interesting questions remain to be studied especially at supra-saturation densities. In this lecture note, after a brief summary of the currently available…
▽ More
Significant progress has been made in recent years in constraining nuclear symmetry energy at and below the saturation density of nuclear matter using data from both terrestrial nuclear experiments and astrophysical observations. However, many interesting questions remain to be studied especially at supra-saturation densities. In this lecture note, after a brief summary of the currently available constraints on nuclear symmetry energy near the saturation density we first discuss the relationship between the symmetry energy and the isopin and momentum dependence of the single-nucleon potential in isospin-asymmetric nuclear medium. We then discuss several open issues regarding effects of the tensor force induced neutron-proton short-range correlation (SRC) on nuclear symmetry energy. Finally, as an example of the impacts of nuclear symmetry energy on properties of neutron stars and gravitational waves, we illustrate effects of the high-density symmetry energy on the tidal polarizability of neutron stars in coalescing binaries.
△ Less
Submitted 5 December, 2012;
originally announced December 2012.
-
Constraining the High-Density Behavior of Nuclear Symmetry Energy with the Tidal Polarizability of Neutron Stars
Authors:
F. J. Fattoyev,
J. Carvajal,
W. G. Newton,
Bao-An Li
Abstract:
Using a set of model equations of state satisfying the latest constraints from both terrestrial nuclear experiments and astrophysical observations as well as state-of-the-art nuclear many-body calculations of the pure neutron matter equation of state, the tidal polarizability of canonical neutron stars in coalescing binaries is found to be a very sensitive probe of the high-density behavior of nuc…
▽ More
Using a set of model equations of state satisfying the latest constraints from both terrestrial nuclear experiments and astrophysical observations as well as state-of-the-art nuclear many-body calculations of the pure neutron matter equation of state, the tidal polarizability of canonical neutron stars in coalescing binaries is found to be a very sensitive probe of the high-density behavior of nuclear symmetry energy which is among the most uncertain properties of dense neutron-rich nucleonic matter. Moreover, it changes less than $\pm 10%$ by varying various properties of symmetric nuclear matter and symmetry energy around the saturation density within their respective ranges of remaining uncertainty.
△ Less
Submitted 11 October, 2012;
originally announced October 2012.
-
Pure Neutron Matter Constraints and Nuclear Symmetry Energy
Authors:
F. J. Fattoyev,
W. G. Newton,
Jun Xu,
Bao-An Li
Abstract:
In this review, we will discuss the results of our recent work to study the general optimization of the pure isovector parameters of the popular relativistic mean-field (RMF) and Skyrme-Hartree-Fock (SHF) nuclear energy-density functionals (EDFs), using constraints on the pure neutron matter (PNM) equation of state (EoS) from recent {\sl ab initio} calculations. By using RMF and SHF parameterizati…
▽ More
In this review, we will discuss the results of our recent work to study the general optimization of the pure isovector parameters of the popular relativistic mean-field (RMF) and Skyrme-Hartree-Fock (SHF) nuclear energy-density functionals (EDFs), using constraints on the pure neutron matter (PNM) equation of state (EoS) from recent {\sl ab initio} calculations. By using RMF and SHF parameterizations that give equivalent predictions for ground-state properties of doubly magic nuclei and properties of symmetric nuclear matter (SNM) and PNM, we found that such optimization leads to broadly consistent symmetry energy $J$ and its slope parameter $L$ at saturation density within a tight range of $σ(J) < 2$ MeV and $σ(L) < 6$ MeV. We demonstrate that a clear model dependence shows up (a) in the curvature parameter of the symmetry energy $K_{\rm sym}$, (b) the symmetry energy at supra-saturation densities, and (c) the radius of neutron stars.
△ Less
Submitted 12 September, 2012;
originally announced September 2012.
-
Generic Constraints on the Relativistic Mean-Field and Skyrme-Hartree-Fock Models from the Pure Neutron Matter Equation of State
Authors:
F. J. Fattoyev,
W. G. Newton,
Jun Xu,
Bao-An Li
Abstract:
We study the nuclear symmetry energy S(rho) and related quantities of nuclear physics and nuclear astrophysics predicted generically by relativistic mean-field (RMF) and Skyrme-Hartree-Fock (SHF) models. We establish a simple prescription for preparing equivalent RMF and SHF parametrizations starting from a minimal set of empirical constraints on symmetric nuclear matter, nuclear binding energy an…
▽ More
We study the nuclear symmetry energy S(rho) and related quantities of nuclear physics and nuclear astrophysics predicted generically by relativistic mean-field (RMF) and Skyrme-Hartree-Fock (SHF) models. We establish a simple prescription for preparing equivalent RMF and SHF parametrizations starting from a minimal set of empirical constraints on symmetric nuclear matter, nuclear binding energy and charge radii, enforcing equivalence of their Lorenz effective masses, and then using the pure neutron matter (PNM) equation of state (EoS) obtained from ab-initio calculations to optimize the pure isovector parameters in the RMF and SHF models. We find the resulting RMF and SHF parametrizations give broadly consistent predictions of the symmetry energy J and its slope parameter L at saturation density within a tight range of <~2 MeV and <~6 MeV respectively, but that clear model dependence shows up in the predictions of higher-order symmetry energy parameters, leading to important differences in (a) the slope of the correlation between J and L from the confidence ellipse, (b) the isospin-dependent part of the incompressibility of nuclear matter K_tau, (c) the symmetry energy at supra-saturation densities, and (d) the predicted neutron star radii. The model dependence can lead to about 1-2 km difference in predictions of the neutron star radius given identical predicted values of J, L and symmetric nuclear matter (SNM) saturation properties. Allowing the full freedom in the effective masses in both models leads to constraints of 30<~J<~31.5 MeV, 35<~L<~60 MeV, -330<~K_tau<~-216 MeV for the RMF model as a whole and 30<~J<~33 MeV, 28<~L<~65 MeV, -420<~K_tau<~-325 MeV for the SHF model as a whole. Notably, given PNM constraints, these results place RMF and SHF models as a whole at odds with some constraints on K_tau inferred from giant monopole resonance and neutron skin experimental results.
△ Less
Submitted 23 August, 2012; v1 submitted 3 May, 2012;
originally announced May 2012.
-
The nuclear symmetry energy, the inner crust, and global neutron star modeling
Authors:
William G. Newton,
Michael Gearheart,
Josh Hooker,
Bao-An Li
Abstract:
The structure and composition of the inner crust of neutron stars, as well as global stellar properties such as radius and moment of inertia, have been shown to correlate with parameters characterizing the symmetry energy of nuclear matter such as its magnitude J and density dependence L at saturation density. It is thus mutually beneficial to nuclear physicists and astrophysicists to examine the…
▽ More
The structure and composition of the inner crust of neutron stars, as well as global stellar properties such as radius and moment of inertia, have been shown to correlate with parameters characterizing the symmetry energy of nuclear matter such as its magnitude J and density dependence L at saturation density. It is thus mutually beneficial to nuclear physicists and astrophysicists to examine the combined effects of such correlations on potential neutron star observables in the light of recent experimental and theoretical constraints on J, L, and relationships between them. We review some basic correlations between these nuclear and astrophysical observables, and illustrate the impact of recent progress in constraining the J-L parameter space on the composition of the inner crust, crust-core transition density and pressure, and extent of the hypothesized pasta region. We use a simple compressible liquid drop model in conjunction with a simple model of nuclear matter which allows for independent, smooth, variation of the J and L. We extend the model into the core using the same nuclear matter model to explore the effects on global crust and core properties, and on potential observables such as crust oscillation frequencies and mechanically supported crust deformation. Throughout we illustrate the importance of the relationship between J and L implicit in a particular model of nuclear matter to the predictions of neutron star properties.
△ Less
Submitted 8 December, 2011;
originally announced December 2011.
-
Sensitivity of the neutron star r-mode instability window to the density dependence of the nuclear symmetry energy
Authors:
De-Hua Wen,
W. G. Newton,
Bao-An Li
Abstract:
Using a simple model of a neutron star with a perfectly rigid crust constructed and a set of crust and core equations of state that span the range of nuclear experimental uncertainty in the density dependence of the symmetry energy from 25 MeV (soft EOS) to 115 MeV (stiff EOS), we calculate the instability window for the onset of the Chandrasekhar-Friedmann-Schutz (CFS) instability in r-mode oscil…
▽ More
Using a simple model of a neutron star with a perfectly rigid crust constructed and a set of crust and core equations of state that span the range of nuclear experimental uncertainty in the density dependence of the symmetry energy from 25 MeV (soft EOS) to 115 MeV (stiff EOS), we calculate the instability window for the onset of the Chandrasekhar-Friedmann-Schutz (CFS) instability in r-mode oscillations for canonical neutron stars (1.4 M_{\odot}) and massive neutron stars (2.0 M_{\odot}). In these models the crust-core transition density, and thus crustal thickness, is calculated consistently with the core equation of state (EOS). For the canonical neutron star, the lower bound of the r-mode instability window is reduced in frequency by \approx150 Hz from the softest to the stiffest symmetry energy used, independent of mass and temperature. The instability window also drops by \approx 100 Hz independent of EOS when the mass is raised from 1.4 M_{\odot} to 2.0 M_{\odot}. Where temperature estimates are available, the observed neutron stars in low mass X-ray binaries (LMXBs) have frequencies below the instability window for the 1.4 M_{\odot} models, while some LMXBs fall within the instability window for 2.0 M_{\odot} stars if the symmetry energy is relatively stiff, indicating that a softer symmetry energy is more consistent with observations within this model. Thus we conclude that smaller values of L help stabilize neutron stars against runaway r-mode oscillations. The critical temperature, below which no star can reach the instability window without exceeding its Kepler frequency, varies by nearly an order of magnitude from soft to stiff symmetry energies. When the crust thickness and core EOS are treated consistently, a thicker crust corresponds to a lower critical temperature, the opposite result to previous studies in which the transition density was independent of the core EOS.
△ Less
Submitted 1 January, 2012; v1 submitted 27 October, 2011;
originally announced October 2011.
-
A survey of the parameter space of the compressible liquid drop model as applied to the neutron star inner crust
Authors:
W. G. Newton,
M. Gearheart,
Bao-An Li
Abstract:
We present a systematic survey the range of predictions of the neutron star inner crust composition, crust-core transition densities and pressures, and density range of the nuclear `pasta' phases at the bottom of the crust provided by the compressible liquid drop model in the light of current experimental and theoretical constraints on model parameters. Using a Skyrme-like model for nuclear matter…
▽ More
We present a systematic survey the range of predictions of the neutron star inner crust composition, crust-core transition densities and pressures, and density range of the nuclear `pasta' phases at the bottom of the crust provided by the compressible liquid drop model in the light of current experimental and theoretical constraints on model parameters. Using a Skyrme-like model for nuclear matter, we construct baseline sequences of crust models by consistently varying the density dependence of the bulk symmetry energy at nuclear saturation density, $L$, under two conditions: (i) that the magnitude of the symmetry energy at saturation density $J$ is held constant, and (ii) $J$ correlates with $L$ under the constraint that the pure neutron matter (PNM) EoS satisfies the results of ab-initio calculations at low densities. Such baseline crust models facilitate consistent exploration of the $L$ dependence of crustal properties. The remaining surface energy and symmetric nuclear matter parameters are systematically varied around the baseline, and different functional forms of the PNM EoS at sub-saturation densities implemented, to estimate theoretical `error bars' for the baseline predictions. Inner crust composition and transition densities are shown to be most sensitive to the surface energy at very low proton fractions and to the behavior of the sub-saturation PNM EoS. Recent calculations of the energies of neutron drops suggest that the low-proton-fraction surface energy might be higher than predicted in Skyrme-like models, which our study suggests may result in a greatly reduced volume of pasta in the crust than conventionally predicted.
△ Less
Submitted 18 December, 2012; v1 submitted 18 October, 2011;
originally announced October 2011.
-
Upper limits on the observational effects of nuclear pasta in neutron stars
Authors:
M. Gearheart,
W. G. Newton,
J. Hooker,
Bao-An Li
Abstract:
The effects of the existence of exotic nuclear shapes at the bottom of the neutron star inner crust - nuclear `pasta' - on observational phenomena are estimated by comparing the limiting cases that those phases have a vanishing shear modulus and that they have the shear modulus of a crystalline solid . We estimate the effect on torsional crustal vibrations and on the maximum quadrupole ellipticity…
▽ More
The effects of the existence of exotic nuclear shapes at the bottom of the neutron star inner crust - nuclear `pasta' - on observational phenomena are estimated by comparing the limiting cases that those phases have a vanishing shear modulus and that they have the shear modulus of a crystalline solid . We estimate the effect on torsional crustal vibrations and on the maximum quadrupole ellipticity sustainable by the crust. The crust composition and transition densities are calculated consistently with the global properties, using a liquid drop model with a bulk nuclear equation of state (EoS) which allows a systematic variation of the nuclear symmetry energy. The symmetry energy J and its density dependence L at nuclear saturation density are the dominant nuclear inputs which determine the thickness of the crust, the range of densities at which pasta might appear, as well as global properties such as the radius and moment of inertia. We show the importance of calculating the global neutron star properties on the same footing as the crust EoS, and demonstrate that in the range of experimentally acceptable values of L, the pasta phase can alter the crust frequencies by up to a factor of three, exceeding the effects of superfluidity on the crust modes, and decrease the maximum quadrupole ellipticity sustainable by the crust by up to an order of magnitude. The signature of the pasta phases and the density dependence of the symmetry energy on the potential observables highlights the possibility of constraining the EoS of dense, neutron-rich matter and the properties of the pasta phases using astrophysical observations.
△ Less
Submitted 11 August, 2011; v1 submitted 23 June, 2011;
originally announced June 2011.
-
Imprints of Nuclear Symmetry Energy on Properties of Neutron Stars
Authors:
Bao-An Li,
Lie-Wen Chen,
Michael Gearheart,
Joshua Hooker,
Che Ming Ko,
Plamen G. Krastev,
Wei-Kang Lin,
William G. Newton,
De-Hua Wen,
Chang Xu,
Jun Xu
Abstract:
Significant progress has been made in recent years in constraining the density dependence of nuclear symmetry energy using terrestrial nuclear laboratory data. Around and below the nuclear matter saturation density, the experimental constraints start to merge in a relatively narrow region. At supra-saturation densities, there are, however, still large uncertainties. After summarizing the latest ex…
▽ More
Significant progress has been made in recent years in constraining the density dependence of nuclear symmetry energy using terrestrial nuclear laboratory data. Around and below the nuclear matter saturation density, the experimental constraints start to merge in a relatively narrow region. At supra-saturation densities, there are, however, still large uncertainties. After summarizing the latest experimental constraints on the density dependence of nuclear symmetry energy, we highlight a few recent studies examining imprints of nuclear symmetry energy on the binding energy, energy release during hadron-quark phase transitions as well as the $w$-mode frequency and damping time of gravitational wave emission of neutron stars.
△ Less
Submitted 23 March, 2011;
originally announced March 2011.
-
Constraining the gravitational binding energy of PSR J0737-3039B using terrestrial nuclear data
Authors:
W. G. Newton,
Bao-an Li
Abstract:
We show that the gravitational binding energy of a neutron star of a given mass is correlated with the slope of the nuclear symmetry energy at 1-2 times nuclear saturation density for equations of state without significant softening (i.e., those that predict maximum masses $M_{\rm max} > 1.44M_{\odot}$ in line with the largest accurately measured neutron star mass). Applying recent laboratory co…
▽ More
We show that the gravitational binding energy of a neutron star of a given mass is correlated with the slope of the nuclear symmetry energy at 1-2 times nuclear saturation density for equations of state without significant softening (i.e., those that predict maximum masses $M_{\rm max} > 1.44M_{\odot}$ in line with the largest accurately measured neutron star mass). Applying recent laboratory constraints on the slope of the symmetry energy to this correlation we extract a constraint on the baryon mass of the lower mass member of the double pulsar binary system, PSR J0737-3039B. We compare with independent constraints derived from modeling the progenitor star of J0737-3039B up to and through its collapse under the assumption that it formed in an electron capture supernova. The two sets of constraints are consistent only if $L \lesssim$ 70 MeV.
△ Less
Submitted 10 January, 2010; v1 submitted 12 August, 2009;
originally announced August 2009.
-
Modeling Nuclear Pasta and the Transition to Uniform Nuclear Matter with the 3D-Skyrme-Hartree-Fock Method
Authors:
W. G. Newton
Abstract:
The first results of a new three-dimensional, finite temperature Skyrme-Hartree-Fock+BCS study of the properties of inhomogeneous nuclear matter at densities and temperatures leading to the transition to uniform nuclear matter are presented. A constraint is placed on the two independent components of the quadrupole moment in order to self-consistently explore the shape phase space of nuclear con…
▽ More
The first results of a new three-dimensional, finite temperature Skyrme-Hartree-Fock+BCS study of the properties of inhomogeneous nuclear matter at densities and temperatures leading to the transition to uniform nuclear matter are presented. A constraint is placed on the two independent components of the quadrupole moment in order to self-consistently explore the shape phase space of nuclear configurations. The scheme employed naturally allows effects such as (i) neutron drip, which results in an external neutron gas, (ii) the variety of exotic nuclear shapes expected for extremely neutron heavy nuclei, and (iii) the subsequent dissolution of these nuclei into nuclear matter. In this way, the equation of state can be calculated across phase transitions from lower densities (where one dimensional Hartree-Fock suffices) through to uniform nuclear matter without recourse to interpolation techniques between density regimes described by different physical models.
△ Less
Submitted 8 March, 2009;
originally announced March 2009.
-
A New Study of the Transition to Uniform Nuclear Matter in Neutron Stars and Supernovae
Authors:
W. G. Newton
Abstract:
A comprehensive microscopic study of the properties of bulk matter at densities just below nuclear saturation $ρ_s = 2.5 \sim 10^{14}$ g cm$^{-3}$, zero and finite temperature and high neutron fraction, is outlined, and preliminary results presented. Such matter is expected to exist in the inner crust of neutron stars and during the core collapse of massive stars with $M \gtrsim 8M_{\odot}
A comprehensive microscopic study of the properties of bulk matter at densities just below nuclear saturation $ρ_s = 2.5 \sim 10^{14}$ g cm$^{-3}$, zero and finite temperature and high neutron fraction, is outlined, and preliminary results presented. Such matter is expected to exist in the inner crust of neutron stars and during the core collapse of massive stars with $M \gtrsim 8M_{\odot}
△ Less
Submitted 23 August, 2007;
originally announced August 2007.
-
From Microscales to Macroscales in 3D: Selfconsistent Equation of State for Supernova and Neutron Star Models
Authors:
W. G. Newton,
J. R. Stone,
A. Mezzacappa
Abstract:
First results from a fully self-consistent, temperature-dependent equation of state that spans the whole density range of neutron stars and supernova cores are presented. The equation of state (EoS) is calculated using a mean-field Hartree-Fock method in three dimensions (3D). The nuclear interaction is represented by the phenomenological Skyrme model in this work, but the EoS can be obtained in…
▽ More
First results from a fully self-consistent, temperature-dependent equation of state that spans the whole density range of neutron stars and supernova cores are presented. The equation of state (EoS) is calculated using a mean-field Hartree-Fock method in three dimensions (3D). The nuclear interaction is represented by the phenomenological Skyrme model in this work, but the EoS can be obtained in our framework for any suitable form of the nucleon-nucleon effective interaction. The scheme we employ naturally allows effects such as (i) neutron drip, which results in an external neutron gas, (ii) the variety of exotic nuclear shapes expected for extremely neutron heavy nuclei, and (iii) the subsequent dissolution of these nuclei into nuclear matter. In this way, the equation of state is calculated across phase transitions without recourse to interpolation techniques between density regimes described by different physical models. EoS tables are calculated in the wide range of densities, temperature and proton/neutron ratios on the ORNL NCCS XT3, using up to 2000 processors simultaneously.
△ Less
Submitted 23 August, 2007;
originally announced August 2007.
-
The Double Pulsar J0737--3039: Testing the Neutron Star Equation of State
Authors:
Ph. Podsiadlowski,
J. D. M. Dewi,
P. Lesaffre,
J. C. Miller,
W. G. Newton,
J. R. Stone
Abstract:
The double pulsar J0737--3039 has become an important astrophysical laboratory for testing fundamental physics. Here we demonstrate that the low measured mass of Pulsar B can be used to constrain the equation of state of neutron star matter {\em under the assumption} that it formed in an electron-capture supernova. We show that the observed orbital parameters as well as the likely evolutionary h…
▽ More
The double pulsar J0737--3039 has become an important astrophysical laboratory for testing fundamental physics. Here we demonstrate that the low measured mass of Pulsar B can be used to constrain the equation of state of neutron star matter {\em under the assumption} that it formed in an electron-capture supernova. We show that the observed orbital parameters as well as the likely evolutionary history of the system support such a hypothesis and discuss future refinements that will improve the constraints this test may provide.
△ Less
Submitted 23 June, 2005;
originally announced June 2005.