Dying stars - Neutron stars, black holes, etcetera
Theory of stellar evolution applied to dying stars, neutron stars, dwarf stars, and black holes
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Theory of stellar evolution applied to dying stars, neutron stars, dwarf stars, and black holes
Since August, 2008, the Fermi Gamma-ray Space Telescope has been scanning the sky, producing a full-sky image every three hours. These cosmic gamma-rays come from extreme astrophysical phenomena, many related to exploding stars (supernovae) or what these explosions leave behind: supernova remnants, neutron stars, and black holes. This talk uses sample Fermi results, plus simple demonstrations, to illustrate the exotic properties of these endpoints of stellar evolution.
Neutron stars were discovered almost 40 years ago, and yet many of their most fundamental properties remain mysteries. There have been many attempts to measure the mass and radius of a neutron star and thereby constrain the equation of state of the dense nuclear matter at their cores. These have been complicated by unknown parameters such as the source distance and burning fractions. A clean, straightforward way to access the neutron star parameters is with high-resolution spectroscopy. I will present the results of searches for gravitationally red-shifted absorption lines from the neutron star atmosphere using XMM-Newton and Chandra.
Neutron star structure and evolution from presupernovae stars including gravitational collapse and neutrino emission processes
Neutron star rotational flattening after formation in supernova implosion
In May 2023, the LIGO-Virgo-KAGRA (LVK) Collaboration reported the likely black hole-neutron star (BHNS) merger GW230529_181500. The signal was observed with high significance in only one detector, limiting constraints on the black hole spin and motivating our study of disruption in this mass regime. That event is expected to be the merger of a 2.5–4.5 M ⊙ primary with a secondary compact object of mass between 1.2–2.0 M ⊙ . This makes it the first BHNS merger with a significant potential for the production of electromagnetic (EM) counterparts, and provides further evidence for compact objects existing within the suspected lower mass gap. To produce post-merger EM transients, the component of the black hole spin aligned with the orbital angular momentum must be sufficiently high, allowing the neutron star to be tidally disrupted. The disrupting BHNS binary may then eject a few percent of a solar mass of matter, leading to an observable kilonova driven by radioactive decays in ejecta, and/or a compact-binary gamma-ray burst (cbGRB) resulting from the formation of an accretion disk and relativistic jet. Determining which mergers lead to disruption of the neutron star is necessary to predict the prevalence of EM signals from BHNS mergers, yet most BHNS simulations so far have been performed far from the minimum spin required for tidal disruption. Here, we use the Spectral Einstein Code to explore the behavior of BHNS mergers in a mass range consistent with GW230529_181500 close to that critical spin, and compare our results against the mass remnant model currently used by the LVK Collaboration to predict the probability of tidal disruption. Furthermore, our numerical results reveal the emergence of non-zero accretion disks even below the predicted NS disruption limit, of low mass but capable of powering cbGRBs. Our results also demonstrate that the remnant mass model underpredicts the disk mass for the DD2 equation of state, while they are within expected modeling errors for SFHo. The disruption limit itself, however, is not found to significantly differ from the predictions of the analytical model, unless remnant masses M rem ≲ 0.001 M ⊙ prove interesting observationally. In all of our simulations, any kilonova signal would be dim and most likely dominated by post-merger disk outflows.
Neutron stars (NSs) probe the high-density regime of the nuclear equation of state (EOS). However, inferring the EOS from observations of NSs is a computationally challenging task. Here, in this work, we efficiently solve this inverse problem by leveraging differential programming in two ways. First, we enable full Bayesian inference in under one hour of wall time on a GPU by using gradient-based samplers, without requiring pretrained machine learning emulators. Moreover, we demonstrate efficient scaling to high-dimensional parameter spaces. Second, we introduce a novel gradient-based optimization scheme that recovers the EOS of a given NS mass-radius curve. We demonstrate how our framework can reveal consistencies or tensions between nuclear physics and astrophysics. First, we show how the breakdown density of a metamodel description of the EOS can be determined from NS observations. Second, we demonstrate how degeneracies in EOS modeling using nuclear empirical parameters can influence the inverse problem during gradient-based optimization. Looking ahead, our approach opens up new theoretical studies of the relation between NS properties and the EOS, while effectively tackling the data analysis challenges brought by future detectors.
Neutron star vibration period and thermal, rotational, vibrational and magnetic energy storage mechanisms in supernova after envelope ejection
Neutron stars, discussing equations of state and stellar interiors, models, atmosphere, cooling, vibration, rotation and magnetic fields
Neutron star superfluid turbulent state, applying dynamic model to conditions in pulsars
Neutron star models, discussing crystalline crust, particle superfluid properties, rotational vortices and relation to Crab Nebula pulsar
We present 3D general-relativistic neutrino-radiation hydrodynamics simulations of two asymmetric binary neutron star mergers producing long-lived neutron stars remnants and spanning a fraction of their cooling time scale. The mergers are characterized by significant tidal disruption with neutron rich material forming a massive disc around the remnant. The latter develops one-armed dynamics that is imprinted in the emitted kilo-Hertz gravitational waves. Angular momentum transport to the disc is initially driven by spiral-density waves and enhanced by turbulent viscosity and neutrino heating on longer timescales. The mass outflows are composed by neutron-rich dynamical ejecta of mass ∼10 −3 –10 −2 M ⊙ followed by a persistent spiral-wave/neutrino-driven wind of ≳ 10 −2 M ⊙ with material spanning a wide range of electron fractions, ∼0.1–0.55. Dynamical ejecta (winds) have fast velocity tails up to ∼0.8 (∼0.4) c. The outflows are further evolved to days timescale using 2D ray-by-ray radiation-hydrodynamics simulations that include an online nuclear network. We find complete r-process yields and identify the production of 56 Ni and the subsequent decay chain to 56 Co and 56 Fe. Synthetic kilonova light curves predict an extended (near-) infrared peak a few days postmerger originating from r-process in the neutron-rich/high-opacity ejecta and UV/optical peaks at a few hours (ten minutes) postmerger originating from weak r-process (free-neutron decay) in the faster ejecta components. Additionally, the fast tail of tidal origin generates kilonova afterglows potentially detectable in radio and X band on a few to ten years time scale. Quantitative effects originating from the tidal disruption merger dynamics are reflected in the multimessenger emissions.
A neutron star is not a place most would want to visit. This dense remnant of a collapsed star has a magnetic field billions of times stronger than Earth's, enough to shuffle your body's molecules long before you even land. The featureless surface is no fun either. Crushing gravity ensures that the star is a near perfect sphere, compressing all matter so that a sand-grain-sized scoop of neutron star material would weigh as much as a battleship on Earth. At least black holes offer the promise of funky singularity, time warps, and the Odyssean temptation to venture beyond a point of no return. What s a journey to a neutron star good for, one might ask? Well, for starters, it offers the possibility of confirming a theorized state of matter called quark-gluon plasma, which likely existed for a moment after the Big Bang and now might only exist in the superdense interiors of neutron stars. Beneath the neutron star crust, a kilometer-thick plate of crystalline matter, lies the great unknown. The popular theory is that the neutron star interior is made up of a neutron superfluid - a fluid without friction. With the help of two NASA satellites - the Rossi X-Ray Timing Explorer and the Chandra X-Ray Observatory - scientists are journeying to the center of a neutron star. Matter might be so compressed there that it breaks down into quarks, the building blocks of protons and neutrons, and gluons, the carrier of the strong nuclear force. To dig inside a neutron star, no simple drill bit will do. Scientists gain insight into the interior through events called glitches, a sudden change in the neutron star s precise spin rate. 'Glitches are one of the few ways we have to study the neutron star interior,' says Frank Marshall of NASA s Goddard Space Flight Center, who has used the Rossi Explorer to follow the escapades of the glitchiest of all neutron stars, dubbed the Big Glitcher and known scientifically as PSR J0537-6910.
Both the core collapse of rotating massive stars, and the coalescence of neutron star (NS) binaries result in the formation of a hot, differentially rotating NS remnant. The timescales over which differential rotation is removed by internal angular-momentum transport processes (viscosity) have key implications for the remnant’s long-term stability and the NS equation of state (EOS). Guided by a nonrotating model of a cooling proto-NS, we estimate the dominant sources of viscosity using an externally imposed angular-velocity profile Ω(r). Although the magneto-rotational instability provides the dominant source of effective viscosity at large radii, convection and/or the Tayler–Spruit dynamo dominate in the core of merger remnants where dΩ/dr ≥ 0. Furthermore, the viscous timescale in the remnant core is sufficiently short that solid-body rotation will be enforced faster than matter is accreted from rotationally supported outer layers. Guided by these results, we develop a toy model for how the merger remnant core grows in mass and angular momentum due to accretion. We find that merger remnants with sufficiently massive and slowly rotating initial cores may collapse to black holes via envelope accretion, even when the total remnant mass is less than the usually considered threshold ≈1.2 M TOV for forming a stable solid-body rotating NS remnant (where M TOV is the maximum nonrotating NS mass supported by the EOS). This qualitatively new picture of the post-merger remnant evolution and stability criterion has important implications for the expected electromagnetic counterparts from binary NS mergers and for multimessenger constraints on the NS EOS
The fate of the binary neutron star (NS) merger remnants hinges sensitively upon the NS equation of state and the threshold mass, M ls , that separates a long-lived from a short-lived NS remnant. The nature of the electromagnetic counterparts is also influenced by the remnant type, particularly in determining whether a gamma-ray burst from a compact binary merger (cbGRB) is of short or long duration. We propose a novel approach to probe Mls by linking it to the estimated observed ratio of long to short cbGRBs. We find that current observations broadly favor a relatively high value for this transition, M ls ≃ 1.3M TOV , for which M TOV ≲ 2.6M ⊙ , consistent with numerical simulations, as also shown here. Furthermore, our results disfavor nuclear physics scenarios that would lead to catastrophic pressure loss at a few times nuclear density and temperatures of tens of MeV, leading to a rapid gravitational collapse of binaries with total mass M ≲ 1.3M TOV . Future individual gravitational wave events with on-axis cbGRBs can further bound Mls.
Neutron stars explore matter at the highest densities in the universe such that their inner cores might undergo a phase transition from hadronic to exotic phases, e.g., quark matter. Such a transition could be associated with nontrivial structures in the density behavior of the speed of sound, such as jumps and sharp peaks. Here, in this work, we employ a physics-agnostic approach to model the density dependence of the speed of sound in neutron stars and study to what extent the existence of nontrivial structures can be inferred from existing astrophysical observations of neutron stars. For this, we exhaustively study different equations of state, including as well those with explicit first-order phase transitions. We obtain a large number of different equation of states which reproduce the same astrophysical observations and obey the same physical constraints, such as mechanical stability and causality. Among them, some have nontrivial structures in the sound speed whereas others do not. We conclude that astrophysical information to date do not require the existence of a phase transition to quark matter in the density range explored in the core of neutron stars.
Triaxial neutron stars can be sources of continuous gravitational radiation detectable by ground-based interferometers. The amplitude of the emitted gravitational wave can be greatly affected by the state of the hydrodynamical fluid flow inside the neutron star. In this work, we examine the most triaxial models along two sequences of constant rest mass, confirming their dynamical stability. We also study the response of a triaxial figure of quasiequilibrium under a variety of perturbations that lead to different fluid flows. Starting from the general relativistic compressible analog of the Newtonian Jacobi ellipsoid, we perform simulations of Dedekind-type flows. We find that in some cases the triaxial neutron star resembles a Riemann-S-type ellipsoid with minor rotation and gravitational wave emission as it evolves towards axisymmetry. The present results highlight the importance of understanding the fluid flow in the interior of a neutron star in terms of its gravitational wave content.