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.
Density and electron fraction of neutron star stuff leftover after a black hole eats the neutron star. Events like this one are believed to power gamma ray bursts, some of the most energetic events in the universe, and be the source of heavy elements like gold and platinum.
The synthesis of hyper-heavy elements is investigated under conditions simulating neutron star environment. The constrained molecular dynamics approach is used to simulate low energy collisions of extremely n-rich nuclei. A new type of the fusion barrier due to a “neutron wind” is observed when the effect of neutron star environment (screening of Coulomb interaction) is introduced implicitly. When introducing also a background of surrounding nuclei, the nuclear fusion becomes possible down to temperatures of 10 8 K and synthesis of extremely heavy and n-rich nuclei appears feasible. A possible existence of hyper-heavy nuclei in a neutron star environment could provide a mechanism of extra coherent neutrino scattering or an additional mechanism, resulting in x-ray burst or a gravitational wave signal and, thus, becoming another crucial process adding new information to the suggested models on neutron star evolution.
Due to their unique set of multimessenger signals, neutron star mergers have emerged as a novel environment for studies of new physics beyond the Standard Model (SM). As a case study, we consider the simplest extension of the SM scalar sector involving a light CP-even scalar singlet S mixing with the SM Higgs boson. These S particles can be produced abundantly in neutron star mergers via the nucleon bremsstrahlung process. Here, we show that the S particles may either be trapped in or stream freely out of the merger remnant, depending on the S mass, its mixing with the SM Higgs boson, and the temperature and baryon density in the merger. In the free-streaming region, the scalar S will provide an extra channel to cool down the merger remnant, with cooling timescales as small as Script O(ms). On the other hand, in the trapped region, the Bose gas of S particles could contribute a larger thermal conductivity than the trapped neutrinos in some parts of the parameter space, thus leading to faster thermal equilibration than expected. Therefore, future observations of the early postmerger phase of a neutron star merger could effectively probe a unique range of the S parameter space, largely complementary to the existing and future laboratory and supernova limits. In view of these results, we hope the merger simulation community will be motivated to implement the effects of light CP-even scalars into their simulations in both the free-streaming and trapped regimes.
In recent years our understanding of the dense matter equation of state (EOS) of neutron stars has significantly improved by analyzing multimessenger data from radio/X-ray pulsars, gravitational wave events, and from nuclear physics constraints. Here we study the additional impact on the EOS from the jointly estimated mass and radius of PSR J0740+6620, presented in Riley et al. by analyzing a combined data set from X-ray telescopes NICER and XMM-Newton. We employ two different high-density EOS parameterizations: a piecewise-polytropic (PP) model and a model based on the speed of sound in a neutron star (CS). Here, at nuclear densities these are connected to microscopic calculations of neutron matter based on chiral effective field theory (EFT) interactions. In addition to the new NICER data for this heavy neutron star, we separately study constraints from the radio timing mass measurement of PSR J0740+6620, the gravitational wave events of binary neutron stars GW190425 and GW170817, and for the latter the associated kilonova AT2017gfo. By combining all these, and the NICER mass-radius estimate of PSR J0030+0451, we find the radius of a 1.4 M ⊙ neutron star to be constrained to the 95% credible ranges 12.33$^{+0.76}_{–0.81}$km (PP model) and 12.18$^{+0.56}_{–0.79}$km (CS model). In addition, we explore different chiral EFT calculations and show that the new NICER results provide tight constraints for the pressure of neutron star matter at around twice saturation density, which shows the power of these observations to constrain dense matter interactions at intermediate densities.
We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run (O3). We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization area of 4480 deg 2 , a median distance of 267 Mpc, and false-alarm rates ranging from 1.5 to 10 -25 yr -1 . The ZTF coverage in the g and r bands had a median enclosed probability of 39%, median depth of 20.8 mag, and median time lag between merger and the start of observations of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UltraViolet/Optical/InfraRed (UVOIR) photometric points, 64 OIR spectra, and three radio images using 17 different telescopes. We find no promising kilonovae (radioactivity-powered counterparts), and we show how to convert the upper limits to constrain the underlying kilonova luminosity function. Initially, we assume that all GW triggers are bona fide astrophysical events regardless of false-alarm rate and that kilonovae accompanying BNS and NSBH mergers are drawn from a common population; later, we relax these assumptions. Assuming that all kilonovae are at least as luminous as the discovery magnitude of GW170817 (-16.1 mag), we calculate that our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than-16.6 mag (the extrapolated peak magnitude of GW170817) and fade at a rate of 1 mag day -1 (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations based on the online classifications, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6 mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, no more than <57% (<89%) of putative kilonovae could be brighter than -16.6 mag assuming flat evolution (fading by 1 mag day -1 ) at the 90% confidence level. If we further take into account the online terrestrial probability for each GW trigger, we find that no more than <68% of putative kilonovae could be brighter than -16.6 mag. Comparing to model grids, we find that some kilonovae must have M ej < 0.03 M ⊙ , X lan > 10 -4 , or φ > 30° to be consistent with our limits. We look forward to searches in the fourth GW observing run; even 17 neutron star mergers with only 50% coverage to a depth of-16 mag would constrain the maximum fraction of bright kilonovae to <25%.
Binary neutron star mergers play an important role in nuclear astrophysics: Their gravitational wave and electromagnetic signals carry information about the equation of state of cold matter above nuclear saturation density, and they may be one of the main sources of r-process elements in the Universe. Neutrino-matter interactions during and after merger impact the properties of these electromagnetic signals and the relative abundances of the produced r-process elements. Existing merger simulations are, however, limited in their ability to realistically model neutrino transport and neutrino-matter interactions. Here, we perform a comparison of the impact of the use of state-of-the art two-moment or Monte Carlo transport schemes on the outcome of merger simulations, for a single binary neutron star system with a short-lived neutron star remnant [(5–10) ms]. We also investigate the use of different reaction rates in the simulations. While the best transport schemes generally agree well on the qualitative impact of neutrinos on the system, differences in the behavior of the high-density regions can significantly impact the collapse time and the properties of the hot tidal arms in this metastable merger remnant. The chosen interaction rates and transport algorithm as well as recent improvements by Radice et al. to the two-moment algorithms can all contribute to changes at the (10–30)% level in the global properties of the merger remnant and outflows. Furthermore, the limitations of previous moment schemes fixed by Radice et al. also appear sufficient to explain the large difference that we observed in the production of heavy-lepton neutrinos in a previous comparison of Monte Carlo and moment schemes in the context of a low-mass binary neutron star system.
The binary-driven hypernova (BdHN) model proposes long gamma-ray bursts (GRBs) originate in binaries composed of a carbon–oxygen (CO) star and a neutron star (NS) companion. The CO core collapse generates a newborn NS and a supernova that triggers the GRB by accreting onto the NSs, rapidly transferring mass and angular momentum to them. This article aims to determine the conditions under which a black hole (BH) forms from NS collapse induced by the accretion and the impact on the GRB's observational properties and taxonomy. We perform three-dimensional, smoothed particle hydrodynamics simulations of BdHNe using up-to-date NS nuclear equations of state, with and without hyperons, and calculate the structure evolution in full general relativity. We assess the binary parameters leading either NS in the binary to the critical mass for gravitational collapse into a BH and its occurrence time, t col . We include a nonzero angular momentum of the NSs and find that t col ranges from a few tens of seconds to hours for decreasing NS initial angular momentum values. BdHNe I are the most compact (about 5 minute orbital period), promptly form a BH, and release ≳10 52 erg of energy. They form NS–BH binaries with tens of kiloyears merger timescales by gravitational-wave emission. BdHNe II and III do not form BHs, and release ~10 50 –10 52 erg and ≲10 50 erg of energy, respectively. They form NS–NS binaries with a range of merger timescales larger than for NS–BH binaries. In some compact BdHNe II, either NS can become supramassive, i.e., above the critical mass of a nonrotating NS. Magnetic braking by a 10 13 G field can delay BH formation, leading to BH–BH or NS–BH with tens of kiloyears merger timescales.
By inverting the observational data of several neutron star observables in the three-dimensional parameter space of the constant speed of sound (CSS) model while fixing all hadronic equation-of-state parameters at their currently known most probable values, we constrain the three parameters of the CSS model and their correlations. Using two lower radius limits of $R$ $2.01$ = $11.41$ km and $R$ $2.01$ = $12.2$ km for PSR J0740+6620 obtained from two independent analyses using different approaches by the Neutron Star Interior Composition Explorer (NICER) Collaboration, the speed of sound squared $c$ $^2_{QM}$ in quark matter is found to have a lower limit of 0.35 or 0.43 in units of $c$ $2$ , above its conformal limit of $c$ $^2_{QM}$ < $1/3$. Furthermore, an approximately linear correlation between the first-order hadron-quark transition density $ρ$ $t$ and its strength $Δε$ is found. Moreover, the presence of twin stars is deemed improbable by the present work.
As neutron stars merge they can approach very high nuclear density. Here, we summarized recent results for the evolution and gravitational wave emission from binary-neutron star mergers using a a variety of nuclear equations of state with and without a crossover transition to quark matter. We discuss how the late time gravitational wave emission from binary neutron star mergers may possibly reveal the existence of a crossover transition to quark matter.
Predicting the properties of the matter ejected during and after a neutron star merger is crucial to our ability to use electromagnetic observations of these mergers to constrain the masses of the neutron stars, the equation of state of dense matter, and the role of neutron star mergers in the enrichment of the Universe in heavy elements. Furthermore, our ability to reliably provide such predictions is however limited by a broad range of factors, including the finite resolution of numerical simulations, their treatment of magnetic fields, neutrinos, and neutrino-matter interactions, and the approximate modeling of the equation of state of dense matter. In this manuscript, we study specifically the role that a small residual eccentricity and different implementations of the same equation of state have on the matter ejected during the merger of a 1.3M ⊙ –1.4M ⊙ binary neutron star system. We find that a residual eccentricity e ~ 0.01, as measured ~ 4–6 orbits before merger, causes O(25%–30%) changes in the amount of ejected mass, mainly due to changes in the amount of matter ejected as a result of core bounces during merger. We note that O(1%) residual eccentricities have regularly been used in binary neutron star merger simulations as proxy for circular binaries, potentially creating an additional source of error in predictions for the mass of the dynamical ejecta.
In this paper we investigate quark deconfinement in neutrons stars and their mergers, focusing on the effects of higher orders for the phase transition between hadronic and quark matter. The different descriptions we use to describe matter microscopically contain varying particle degrees of freedom, including nucleons, hyperons, Delta baryons, and light and strange quarks. We use tabulated equations of state from the CompOSE database in which the quark deconfinement phase transition is described as being first order, and then smooth it out by introducing a percolation, replacing the single first-order phase transition with two transitions of second or third order. We then perform binary neutron-star merger simulations using these new equations of state, focusing on groups of binaries with the same single-star mass, radius, and tidal deformability, but different equations of state. Here, we go on to discuss differences in their evolution, and the ramifications for interpreting future gravitational wave observations and the potential to learn about dense matter.
The interpretation of observations of cooling neutron star crusts in quasipersistent x-ray transients is affected by predictions of the strength of neutrino cooling via crust Urca processes. The strength of crust Urca neutrino cooling depends sensitively on the electron-capture and β -decay ground-state-to-ground-state transition strengths of neutron-rich rare isotopes. Nuclei with a mass number of A = 61 are predicted to be among the most abundant in accreted crusts, and the last remaining experimentally undetermined ground-state-to-ground-state transition strength was the β decay of 61 V. This Letter reports the first experimental determination of this transition strength, a ground-state branching of $8.1^{+4.0}_{- 3.1}%$, corresponding to a log ft value of $5.5^{+0.2}_{-0.2}$. This result was achieved through the measurement of the β -delayed γ rays using the total absorption spectrometer SuN and the measurement of the β -delayed neutron branch using the neutron long counter system NERO at the National Superconducting Cyclotron Laboratory at Michigan State University. This method helps to mitigate the impact of the pandemonium effect in extremely neutron-rich nuclei on experimental results. The result implies that A = 61 nuclei do not provide the strongest cooling in accreted neutron star crusts as expected by some predictions, but that their cooling is still larger compared to most other mass numbers. Finally, only nuclei with mass numbers 31, 33, and 55 are predicted to be cooling more strongly. However, the theoretical predictions for the transition strengths of these nuclei are not consistently accurate enough to draw conclusions on crust cooling. With the experimental approach developed in this work, all relevant transitions are within reach to be studied in the future.
The theory governing the strong nuclear force—quantum chromodynamics—predicts that at sufficiently high energy densities, hadronic nuclear matter undergoes a deconfinement transition to a new phase of quarks and gluons. Although this has been observed in ultrarelativistic heavy-ion collisions, it is currently an open question whether quark matter exists inside neutron stars. By combining astrophysical observations and theoretical ab initio calculations in a model-independent way, we find that the inferred properties of matter in the cores of neutron stars with mass corresponding to 1.4 solar masses (M ⊙ ) are compatible with nuclear model calculations. However, the matter in the interior of maximally massive stable neutron stars exhibits characteristics of the deconfined phase, which we interpret as evidence for the presence of quark-matter cores. For the heaviest reliably observed neutron stars with mass M ≈ 2M ⊙ , the presence of quark matter is found to be linked to the behaviour of the speed of sound cs in strongly interacting matter. If the conformal bound ${c}_{\rm{s}}^{2}\le 1/3$ is not strongly violated, massive neutron stars are predicted to have sizable quark-matter cores. This finding has important implications for the phenomenology of neutron stars and affects the dynamics of neutron star mergers with at least one sufficiently massive participant.
The collision of two neutron stars is a rich source of information about nuclear physics. In particular, the kilonova signal following a merger can help us elucidate the role of neutron stars in nucleosynthesis, and informs us about the properties of matter above nuclear saturation. Approximate modeling of neutrinos remains an important limitation to our ability to make predictions for these observables. Part of the problem is the fermionic nature of neutrinos. By the exclusion principle, the expected value 𝑓 𝜈 for the number of neutrinos in a quantum state is at most 1. Any process producing neutrinos is suppressed by a blocking factor (1 −𝑓 𝜈 ). Recent simulations focused on neutrino physics mostly use a gray two-moment scheme to evolve neutrinos. This evolves integrals of 𝑓 𝜈 over momentum space, preventing direct calculations of blocking factors. Monte Carlo methods may be an attractive alternative, providing access to the full distribution of neutrinos. Their current implementation is, however, inadequate to estimate 𝑓 𝜈 : in our most recent simulations, a single Monte Carlo packet causes, in the worst cases, estimates of 𝑓 𝜈 to jump from 𝑓 𝜈 =0 to 𝑓 𝜈 ∼10 5 . While this is concerning, this brazen violation of the fermionic nature of neutrinos has been largely inconsequential, as the interactions used in simulations avoid direct calculations of 𝑓 𝜈 . We are, however, reaching a level of modeling at which this problem can no longer be ignored. Here, we discuss the relatively simple origin of this issue. We then show that very rough estimates of 𝑓 𝜈 can in theory be obtained in merger simulations, but that they will require a combination of unintuitive weighting schemes for Monte Carlo packets and smoothing of the neutrino distribution at coarser resolution than what the merger simulation uses.
We introduce here mass excesses (ME) of neutron-rich isotopes of Ar through Fe, obtained via time of flight B ρ mass spectrometry at the National Superconducting Cyclotron Laboratory. Our new results have significantly reduced systematic uncertainties relative to a prior analysis, enabling the first determination of ME for Ti 58 , 59 , V 62 , Cr 65 , Mn 67 , 68 , and Fe 69 , 70 . Our results show the N = 34 subshell weaken at Sc and vanish at Ti, along with the absence of an N = 40 subshell at Mn. This further leads to a cooler accreted neutron star crust, highlighting the connection between the structure of nuclei and neutron stars.
We study out-of-thermodynamic-equilibrium effects in neutron-star mergers with 3D general-relativistic neutrino-radiation large-eddy simulations. During mergers, the cores of the neutron stars remain cold (T ~ a few MeV) and out of thermodynamic equilibrium with trapped neutrinos originating from the hot collisional interface between the stars. However, within ~2 to 3 ms matter and neutrinos reach equilibrium everywhere in the remnant massive neutron star. Furthermore, our results show that dissipative effects, such as bulk viscosity, if present, are only active for a short window of time after the merger.
We revisit the problem of inertial r-modes in stratified stars, drawing on a more precise description of the composition stratification in a mature neutron star. The results highlight issues with the traditional approach to the problem, leading us to rethink the computational strategy for the r-modes of nonbarotropic neutron stars. We outline two strategies for dealing with the problem. For moderate to slowly rotating neutron stars the only viable alternative may be to approach the problem numerically from the outset, while a meaningful slow-rotation calculation can be carried out for the fastest known spinning stars (which may be close to being driven unstable by the emission of gravitational waves). We demonstrate that the latter approach leads to a problem close, but not identical, to that for barotropic inertial modes. We also suggest that these reformulations of the problem likely resolve the long-standing problem of singular behavior associated with a corotation point in rotating relativistic neutron stars. This issue needs to be resolved in order to guide future gravitational-wave searches.
Compact binary mergers involving neutron stars can eject a fraction of their mass to space. Being extremely neutron rich, this material undergoes rapid neutron capture nucleosynthesis, and the resulting radioactivity powers fast, short-lived electromagnetic transients known as kilonova or macronova. Such transients are exciting probes of the most extreme physical conditions and their observation signals the enrichment of the Universe with heavy elements. Here the current understanding of the mass ejection mechanisms, the properties of the ejecta, and the resulting radioactive transients are reviewed. The first well-observed event in the aftermath of GW170817 delivered a wealth of insights, but much of today's picture of such events is still based on a patchwork of theoretical studies. Apart from summarizing the current understanding, questions where no consensus has been reached yet are also pointed out, and possible directions for the future research are sketched. In an appendix, a publicly available heating rate library based on the WinNet nuclear reaction network is described, and a simple fit formula to alleviate the implementation in hydrodynamic simulations is provided.