Neutron Stars as Laboratories for Nuclear and Particle Physics
Neutron Stars as Laboratories for Nuclear and Particle Physics
Engineering topics
Publications and source records attributed to Reddy, Sanjay.
Neutron Stars as Laboratories for Nuclear and Particle Physics
Nuclear Physics at m_π=82 MeV and Axion Condensation
The major goals of the TEAMS project are (i) to improve the modeling of many of the suggested astrophysical sites of the r‐process, as well as some potential p‐process sites, (ii) to use our models to support and guide the experimental efforts at FRIB and similar facilities, and (iii) to use our tools to cement our understanding of the roles that massive stars play in galactic chemical evolution. As part of the TEAMS project, Mezzacappa (UTK, Reddy (UW), Surman (ND), and Steiner (UTK) are (i) generating equations of state and neutrino opacities for simulations of core‐collapse supernovae and neutron star mergers which are more faithful to current understanding of the nucleon‐nucleon interaction, (ii) integrating those new equations of state and neutrino opacities into new simulation codes being developed, (iii) performing detailed nuclear reaction network post‐processing to determine abundances from simulations and (iv) quantifying the impact of nuclear physics uncertainties on r‐process abundances.
The accumulation of certain types of dark matter particles in neutron star cores due to accretion over long timescales can lead to the formation of a mini black hole. In this scenario, the neutron star is destabilized and implodes to form a black hole without significantly increasing its mass. When this process occurs in neutron stars in coalescing binaries, one or both stars might be converted to a black hole before they merge. Thus, in the mass range of ~1–2M ⊙ , the Universe might contain three distinct populations of compact binaries: one containing only neutron stars, the second population of only black holes, and a third, mixed population consisting of a neutron star and a black hole. However, it is unlikely to have a mixed population as the various timescales allow for both neutron stars to remain or collapse within a short timescale. In this paper, we explore the capability of future gravitational-wave detector networks, including upgrades of Advanced LIGO and Virgo, and new facilities such as the Cosmic Explorer and Einstein Telescope (XG network), to discriminate between different populations by measuring the effective tidal deformability of the binary, which is zero for binary black holes but nonzero for binary neutron stars. Furthermore, we show that observing the relative abundances of the different populations can be used to infer the timescale for neutron stars to implode into black holes, and in turn, provide constraints on the particle nature of dark matter. Furthermore, the XG network will infer the implosion timescale to within an accuracy of 0.01 Gyr at 90% credible interval and determine the dark matter mass and interaction cross section to within a factor of 2 GeV and 10 cm –2 , respectively.
With the improved accuracy of neutron star observational data, it is necessary to derive new equation of state where the crust and the core are consistently calculated within a unified approach. For this purpose we describe non-uniform matter in the crust of neutron stars employing a compressible liquid-drop model, where the bulk and the neutron fluid terms are given from the same model as the one describing uniform matter present in the core. We then generate a set of fifteen unified equations of state for cold catalyzed neutron stars built on realistic modelings of the nuclear interaction, which belongs to two main groups: the first one derives from the phenomenological Skyrme interaction and the second one from χEFT Hamiltonians. The confrontation of these model predictions allows us to investigate the model dependence for the crust properties, and in particular the effect of neutron matter at low density. Here, the new set of unified equations of state is available at the CompOSE repository.
Here, the NASA telescope NICER has recently measured x-ray emissions from the heaviest of the precisely known two-solar mass neutron stars, PSR J0740 + 6620. Analysis of the data suggests that PSR J0740 + 6620 has a radius in the range of R 2.0 ≈ (11.4 – 16.1) km at the 68% credibility level. In this article, we study the implications of this analysis for the sound speed in the high-density inner cores by using recent chiral effective field theory (χEFT) calculations of the equation of state at next-to-next-to-next-to-leading order to describe outer regions of the star at modest density. We find that the lower bound on the maximum speed of sound in the inner core, min {c$^{2}_{s,max}$}, increases rapidly with the radius of massive neutron stars. If χ EFT remains an efficient expansion for nuclear interactions up to about twice the nuclear saturation density, R 2.0 ≥ 13 km requires min {c$^{2}_{s,max}$} ≥ 0.562 and 0.442 at the 68% and 95% credibility level, respectively.
Here we perform general relativistic one-dimensional supernova (SN) simulations to identify observable signatures of enhanced axion emission from the pion-induced reaction π - + p →n + a inside a newly born protoneutron star (PNS). We focus on the early evolution after the onset of the supernova explosion to predict the temporal and spectral features of the neutrino and axion emission during the first 10 s. Pions are included as explicit new degrees of freedom in hot and dense matter. Their thermal population and their role in axion production are both determined consistently to include effects due to their interactions with nucleons. For a wide range of ambient conditions encountered inside a PNS, we find that the pion-induced axion production dominates over nucleon-nucleon bremsstrahlung processes. By consistently including the role of pions on the dense matter equation of state and on the energy loss, our simulations predict robust discernible features of neutrino and axion emission from a galactic supernova that can be observed in terrestrial detectors. For axion couplings that are compatible with current bounds, we find a significant suppression with time of the neutrino luminosity during the first 10 s. This suggests that current bounds derived from the neutrino signal from SN 1987A can be improved and that future galactic supernovae may provide significantly more stringent constraints.
We combine the equation of state of dense matter up to twice nuclear saturation density 𝑛 sat obtained using chiral effective field theory (𝜒EFT) and recent observations of neutron stars to gain insights about the high-density matter encountered in their cores. A key element in our study is the recent Bayesian analysis of correlated EFT truncation errors based on order-by-order calculations up to next-to-next-to-next-to-leading order in the 𝜒EFT expansion. We refine the bounds on the maximum mass imposed by causality at high densities and provide stringent limits on the maximum and minimum radii of ∼1.4M ⊙ and ∼2.0M ⊙ stars. Including 𝜒EFT predictions from 𝑛sat to 2𝑛 sat reduces the permitted ranges of the radius of a ~1.4M ⊙ star, 𝑅 1.4 , by ∼3.5km. If observations indicate 𝑅 1.4 < 11.2km, then our study implies that either the squared speed of sound 𝑐$^{2}_{𝑠}$ > 1/2 for densities above 2𝑛 sat or that 𝜒EFT breaks down below 2𝑛 sat . Here, we also comment on the nature of the secondary compact object in GW190814 with mass ≃ 2.6M ⊙ and discuss the implications of massive neutron stars > 2.1M ⊙ (2.6M ⊙ ) in future radio and gravitational-wave searches. Some form of strongly interacting matter with 𝑐$^{2}_{𝑠}$ > 0.35(0.55) must be realized in the cores of such massive neutron stars. In the absence of phase transitions below 2𝑛 sat , the small tidal deformability inferred from GW170817 lends support for the relatively small pressure predicted by 𝜒EFT for the baryon density 𝑛 B in the range 1–2𝑛 sat . Together they imply that the rapid stiffening required to support a high maximum mass should occur only when 𝑛 B ≳ 1.5–1.8 𝑛 sat .
Recent observations of neutron stars with gravitational waves and X-ray timing provide unprecedented access to the equation of state (EoS) of cold dense matter at densities difficult to realize in terrestrial experiments. At the same time, predictions for the EoS equipped with reliable uncertainty estimates from chiral effective field theory ( x EFT) allow us to bound our theoretical ignorance. In this work, we analyze astrophysical data using a nonparametric representation of the neutron-star EoS conditioned on x EFT to directly constrain the underlying physical properties of the compact objects without introducing modeling systematics. We discuss how the data alone constrain the EoS at high densities when we condition on x EFT at low densities. Here, we also demonstrate how to exploit astrophysical data to directly test the predictions of x EFT for the EoS up to twice nuclear saturation density, in order to estimate the density at which these predictions might break down. We nd that the existence of massive pulsars, gravitational waves from GW170817, and NICER observations of PSR J0030+0451 favor x EFT predictions for the EoS up to nuclear saturation density over a more agnostic analysis by as much as a factor of 7 for the quantum Monte Carlo (QMC) calculations used in this work. While x EFT predictions using QMC are fully consistent with gravitational-wave data up to twice nuclear saturation density, NICER observations suggest that the EoS stiffens relative to these predictions at or slightly above nuclear saturation density. Additionally, for these QMC calculations, we marginalize over the uncertainty in the density at which x EFT begins to break down, constraining the radius of a 1.4M $\bigodot$ neutron star to R 1.4 = $11.40^{+1.38}_{–1.04}$ ($12.54^{+0.71}_{–0.63}$) km and the pressure at twice nuclear saturation density to p(2n sat ) = $14.2^{+18.1}_{–8.4}$ ($28.7^{+15.3}_{–15.0}$) MeV=fm 3 with massive pulsar and gravitational-wave (and NICER) data.
Here we study the role of pions in hot dense matter encountered in astrophysics. We find that strong interactions enhance the number density of negatively charged pions and that this enhancement can be calculated reliably for a relevant range of density and temperature by using the virial expansion. We assess the influence of pions and muons on the equation of state (EOS) and weak-interaction rates in hot dense matter. We find that thermal pions increase the proton fraction and soften the EOS. We also find that charged current weak reactions involving pions and muons $\bar{ν}_ μ + μ^– → π$ and $ν_μ + π^– → μ^–$ make an important contribution to the opacity of muon neutrinos. This could influence the dynamics of core-collapse supernovae and neutron-star mergers. Finally, we note that pion-nucleon reactions can alter the evolution of the proton fraction when weak interactions are not in equilibrium.
The properties of neutron stars are determined by the nature of the matter that they contain. These properties can be constrained by measurements of the star’s size. We obtain stringent constraints on neutron-star radii by combining multimessenger observations of the binary neutron-star merger GW170817 with nuclear theory that best accounts for density-dependent uncertainties in the equation of state. We construct equations of state constrained by chiral effective field theory and marginalize over these using the gravitational-wave observations. Combining this with the electromagnetic observations of the merger remnant that imply the presence of a short-lived hypermassive neutron star, we find that the radius of a 1.4 M⊙ neutron star is $R_{1.4M_⊙} = 11.0^{+0.9}_{-0.6} km$ (90% credible interval). Finally, using this constraint, we show that neutron stars are unlikely to be disrupted in neutron star–black hole mergers; subsequently, such events will not produce observable electromagnetic emission.
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 AM 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.