Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “NEUTRON STARS”

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.

At least 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The Binary–Host Connection: Astrophysics of Gravitational-Wave Binaries from Host Galaxy Properties

Gravitational waves from the merger of binary neutron stars (BNSs) are accompanied by electromagnetic counterparts, making it possible to identify the associated host galaxy. In this work, we explore how properties of the hosts relate to the astrophysical processes leading to the mergers. It is thought that the BNS merger rate within a galaxy at a given epoch depends primarily on the galaxy's star formation history, as well as the underlying merger time-delay distribution of the binary systems. The stellar history of a galaxy, meanwhile, depends on the cosmological evolution of the galaxy through time, and is tied to the growth of structure in the universe. We study the hosts of BNS mergers in the context of structure formation by populating the UniverseMachine simulations with gravitational wave (GW) events, based on a simple time-delay model. We find that different time-delay distributions predict different properties of the associated host galaxies, including the distributions of stellar mass, star formation rate, halo mass, and local and large-scale clustering of hosts. Moreover, BNSs merging today with short delay times occur preferentially in hosts with high star formation rates, while those with long delay times live in dense regions within massive halos that have low star formation. Here, we show that with ${ \mathcal O }(10)$ events from current GW detector networks, it is possible to make preliminary distinctions between formation channels which trace stellar mass, halo mass, or star formation rate. We also find that strategies to follow-up GW events with electromagnetic telescopes can be significantly optimized using the clustering properties of their hosts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Nucleosynthesis of an 11.8 M ⊙ Supernova with 3D Simulation of the Inner Ejecta: Overall Yields and Implications for Short-lived Radionuclides in the Early Solar System

Based on a 3D supernova simulation of an 11.8 M ⊙ progenitor model with an initial solar composition, we study the nucleosynthesis using tracers covering the innermost 0.1 M ⊙ of the ejecta. These ejecta are mostly proton-rich and contribute significant amounts of 45 Sc and 64 Zn. The production of heavier isotopes is sensitive to the electron fraction and hence the neutrino emission from the proto-neutron star. The yields of these isotopes are rather uncertain due to the approximate neutrino transport used in the simulation. In order to obtain the total yields for the whole supernova, we combine the results from the tracers with those for the outer layers from a suitable 1D model. Using the yields of short-lived radionuclides (SLRs), we explore the possibility that an 11.8 M ⊙ supernova might have triggered the formation of the solar system and provided some of the SLRs measured in meteorites. In particular, we discuss two new scenarios that can account for at least the data on 41 Ca, 53 Mn, and 60 Fe without exceeding those on the other SLRs.

Astronomy & Astrophysics↗
DOE OSTI2021

The Multipolar Magnetic Field of the Millisecond Pulsar PSR J0030+0451

Modeling of the NICER X-ray waveform of the pulsar PSR J0030+0451, aimed at constraining the neutron star mass and radius, has inferred surface hot spots (the magnetic polar caps) that imply significantly nondipolar magnetic fields. To this end, we investigate magnetic field configurations that comprise offset dipole-plus-quadrupole components using a static vacuum field and force-free global magnetosphere models. Taking into account the compactness and observer angle values provided by Miller et al. and Riley et al., we compute geodesics from the observer plane to the polar caps to compute the resulting X-ray light curve. We explore, through Markov Chain Monte Carlo techniques, the detailed magnetic field configurations that can reproduce the observed X-ray light curve and have discovered degeneracies, i.e., diverse field configurations, which can provide sufficient descriptions of the NICER X-ray waveforms. Having obtained the force-free field structures, we then compute the corresponding synchronous γ-ray light curves following Kalapotharakos et al.; these we compare to those obtained by Fermi-LAT, to provide models consistent with both the X-ray and the γ-ray data, thereby restricting further the multipole field parameters. An essential aspect of this approach is the proper computation of the relative phase between the synchronous X- and γ-ray light curves. We conclude with a discussion of the broader implications of our study.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Axisymmetric Radiative Transfer Models of Kilonovae

The detailed observations of GW170817 proved for the first time directly that neutron star mergers are a major production site of heavy elements. The observations could be fit by a number of simulations that qualitatively agree, but can quantitatively differ (e.g., in total r-process mass) by an order of magnitude. We categorize kilonova ejecta into several typical morphologies motivated by numerical simulations, and apply a radiative transfer Monte Carlo code to study how the geometric distribution of the ejecta shapes the emitted radiation. We find major impacts on both spectra and light curves. The peak bolometric luminosity can vary by two orders of magnitude and the timing of its peak by a factor of five. These findings provide the crucial implication that the ejecta masses inferred from observations around the peak brightness are uncertain by at least an order of magnitude. Mixed two-component models with lanthanide-rich ejecta are particularly sensitive to geometric distribution. A subset of mixed models shows very strong viewing angle dependence due to lanthanide "curtaining," which persists even if the relative mass of lanthanide-rich component is small. The angular dependence is weak in the rest of our models, but different geometric combinations of the two components lead to a highly diverse set of light curves. We identify geometry-dependent P Cygni features in late spectra that directly map out strong lines in the simulated opacity of neodymium, which can help to constrain the ejecta geometry and to directly probe the r-process abundances.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Stellar Mass Black Hole Formation and Multimessenger Signals from Three-dimensional Rotating Core-collapse Supernova Simulations

We present self-consistent 3D core-collapse supernova simulations of a 40 ${M}_{\odot }$ progenitor model using the isotropic diffusion source approximation for neutrino transport and an effective general relativistic potential up to `0.9 s postbounce. We consider three different rotational speeds with initial angular velocities of ${{\rm{\Omega }}}_{0}=0$, 0.5, and 1 rad s -1 and investigate the impact of rotation on shock dynamics, black hole (BH) formation, and gravitational wave (GW) signals. The rapidly rotating model undergoes an early explosion at ~250 ms postbounce and shows signs of the low $T/| W| $ instability. We do not find BH formation in this model within ~460 ms postbounce. In contrast, we find BH formation at 776 ms postbounce and 936 ms postbounce for the nonrotating and slowly rotating models, respectively. The slowly rotating model explodes at ~650 ms postbounce, and the subsequent fallback accretion onto the proto–neutron star (PNS) results in BH formation. In addition, the standing accretion shock instability induces rotation of the PNS in the model that started with a nonrotating progenitor. Assuming conservation of specific angular momentum during BH formation, this corresponds to a BH spin parameter of a = J/M = 0.046. However, if no explosion sets in, all the angular momentum will eventually be accreted by the BH, resulting in a nonspinning BH. The successful explosion of the slowly rotating model drastically slows down the accretion onto the PNS, allowing continued cooling and contraction that results in an extremely high GW frequency (f ~ 3000 Hz) at BH formation, while the nonrotating model generates GW signals similar to our corresponding 2D simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

r-Process Radioisotopes from Near-Earth Supernovae and Kilonovae

The astrophysical sites where r-process elements are synthesized remain mysterious: it is clear that neutron star mergers (kilonovae (KNe)) contribute, and some classes of core-collapse supernovae (SNe) are also possible sources of at least the lighter r-process species. The discovery of 60Fe on the Earth and Moon implies that one or more astrophysical explosions have occurred near the Earth within the last few million years, probably SNe. Intriguingly, 244 Pu has now been detected, mostly overlapping with 60 Fe pulses. However, the 244 Pu flux may extend to before 12 Myr ago, pointing to a different origin. Here, motivated by these observations and difficulties for r-process nucleosynthesis in SN models, we propose that ejecta from a KN enriched the giant molecular cloud that gave rise to the Local Bubble, where the Sun resides. Accelerator mass spectrometry (AMS) measurements of 244 Pu and searches for other live isotopes could probe the origins of the r-process and the history of the solar neighborhood, including triggers for mass extinctions, e.g., that at the end of the Devonian epoch, motivating the calculations of the abundances of live r-process radioisotopes produced in SNe and KNe that we present here. Given the presence of 244 Pu, other r-process species such as 93 Zr, 107 Pd, 129 I, 135 Cs, 182 Hf, 236 U, 237 Np, and 247 Cm should be present. Their abundances and well-resolved time histories could distinguish between the SN and KN scenarios, and we discuss prospects for their detection in deep-ocean deposits and the lunar regolith. We show that AMS 129 I measurements in Fe–Mn crusts already constrain a possible nearby KN scenario.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Radio and X-Ray Observations of the Luminous Fast Blue Optical Transient AT 2020xnd

We present deep X-ray and radio observations of the fast blue optical transient (FBOT) AT 2020xnd/ZTF 20acigmel at z = 0.2433 from 13 days to 269 days after explosion. AT 2020xnd belongs to the category of optically luminous FBOTs with similarities to the archetypal event AT 2018cow. AT 2020xnd shows luminous radio emission reaching L ν ≈ 8 × 10 29 erg s -1 Hz -1 at 20 GHz and 75 days post-explosion, accompanied by luminous and rapidly fading soft X-ray emission peaking at L X ≈ 6 × 10 42 erg s -1 . Interpreting the radio emission in the context of synchrotron radiation from the explosion's shock interaction with the environment, we find that AT 2020xnd launched a high-velocity outflow (v ~ 0.1c–0.2c) propagating into a dense circumstellar medium (effective $\dot{M}\approx {10}^{-3}\,{M}_{\odot }$ yr -1 for an assumed wind velocity of v w = 1000 km s -1 ). Similar to AT 2018cow, the detected X-ray emission is in excess compared to the extrapolated synchrotron spectrum and constitutes a different emission component, possibly powered by accretion onto a newly formed black hole or neutron star. These properties make AT 2020xnd a high-redshift analog to AT 2018cow, and establish AT 2020xnd as the fourth member of the class of optically luminous FBOTs with luminous multiwavelength counterparts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

IceCube Search for Neutrinos Coincident with Gravitational Wave Events from LIGO/Virgo Run O3

Using data from the IceCube Neutrino Observatory, we searched for high-energy neutrino emission from the gravitational-wave events detected by the advanced LIGO and Virgo detectors during their third observing run. We did a low-latency follow-up on the public candidate events released during the detectors' third observing run and an archival search on the 80 confident events reported in the GWTC-2.1 and GWTC-3 catalogs. An extended search was also conducted for neutrino emission on longer timescales from neutron star containing mergers. Follow-up searches on the candidate optical counterpart of GW190521 were also conducted. We used two methods; an unbinned maximum likelihood analysis and a Bayesian analysis using astrophysical priors, both of which were previously used to search for high-energy neutrino emission from gravitational-wave events. No significant neutrino emission was observed by any analysis, and upper limits were placed on the time-integrated neutrino flux as well as the total isotropic equivalent energy emitted in high-energy neutrinos.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Search for Ultra-high-energy Photons from Gravitational Wave Sources with the Pierre Auger Observatory

A search for time-directional coincidences of ultra-high-energy (UHE) photons above 10 EeV with gravitational wave (GW) events from the LIGO/Virgo runs O1 to O3 is conducted with the Pierre Auger Observatory. Due to the distinctive properties of photon interactions and to the background expected from hadronic showers, a subset of the most interesting GW events is selected based on their localization quality and distance. Time periods of 1000 s around and 1 day after the GW events are analyzed. No coincidences are observed. Upper limits on the UHE photon fluence from a GW event are derived that are typically at ~7 MeV cm -2 (time period 1000 s) and ~35 MeV cm -2 (time period 1 day). Due to the proximity of the binary neutron star merger GW170817, the energy of the source transferred into UHE photons above 40 EeV is constrained to be less than 20% of its total GW energy. These are the first limits on UHE photons from GW sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Comparing Early Evolution of Flames in X-Ray Bursts in Two and Three Dimensions

We explore the early evolution of flame ignition and spreading on the surface of a neutron star in three dimensions, in the context of X-ray bursts. We look at the nucleosynthesis and morphology of the burning front and compare to two-dimensional axisymmetric simulations to gauge how important a full three-dimensional treatment of the flame is for the early dynamics. Finally, we discuss the progress toward full-star resolved flame simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Multimessenger Diagnostics of the Engine behind Core-collapse Supernovae

Abstract Core-collapse supernova explosions play a wide role in astrophysics by producing compact remnants (neutron stars or black holes) and the synthesis and injection of many heavy elements into their host galaxy. Because they are produced in some of the most extreme conditions in the universe, they can also probe physics in extreme conditions (matter at nuclear densities and extreme temperatures and magnetic fields). To quantify the impact of supernovae on both fundamental physics and our understanding of the universe, we must leverage a broad set of observables of this engine. In this paper, we study a subset of these probes using a suite of one-dimensional, parameterized mixing models: ejecta remnants from supernovae, ultraviolet, optical and infrared light curves, and transient gamma-ray emission. We review the other diagnostics and show how the different probes tie together to provide a more clear picture of the supernova engine. Join us in improving and evolving this document through active community engagement. Instructions are provided at this link: https://github.com/clfryer/MM-SNe .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Designing an Optimal Kilonova Search Using DECam for Gravitational-wave Events

We address the problem of optimally identifying all kilonovae detected via gravitational-wave emission in the upcoming LIGO/Virgo/KAGRA observing run, O4, which is expected to be sensitive to a factor of ~7 more binary neutron star (BNS) alerts than previously. Electromagnetic follow-up of all but the brightest of these new events will require >1 m telescopes, for which limited time is available. We present an optimized observing strategy for the DECam during O4. We base our study on simulations of gravitational-wave events expected for O4 and wide-prior kilonova simulations. We derive the detectabilities of events for realistic observing conditions. We optimize our strategy for confirming a kilonova while minimizing telescope time. For a wide range of kilonova parameters, corresponding to a fainter kilonova compared to GW170817/AT 2017gfo, we find that, with this optimal strategy, the discovery probability for electromagnetic counterparts with the DECam is ~80% at the nominal BNS gravitational-wave detection limit for O4 (190 Mpc), which corresponds to an ~30% improvement compared to the strategy adopted during the previous observing run. For more distant events (~330 Mpc), we reach an ~60% probability of detection, a factor of ~2 increase. For a brighter kilonova model dominated by the blue component that reproduces the observations of GW170817/AT 2017gfo, we find that we can reach ~90% probability of detection out to 330 Mpc, representing an increase of ~20%, while also reducing the total telescope time required to follow up events by ~20%.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Could a Kilonova Kill: A Threat Assessment

Abstract Binary neutron star mergers produce high-energy emissions from several physically different sources, including a gamma-ray burst (GRB) and its afterglow, a kilonova (KN), and, at late times, a remnant many parsecs in size. Ionizing radiation from these sources can be dangerous for life on Earth-like planets when located too close. Work to date has explored the substantial danger posed by the GRB to on-axis observers; here we focus instead on the potential threats posed to nearby off-axis observers. Our analysis is based largely on observations of the GW170817/GRB 170817A multi-messenger event, as well as theoretical predictions. For baseline KN parameters, we find that the X-ray emission from the afterglow may be lethal out to ∼1 pc and the off-axis gamma-ray emission may threaten a range out to ∼4 pc, whereas the greatest threat comes years after the explosion, from the cosmic rays accelerated by the KN blast, which can be lethal out to distances up to ∼11 pc. The distances quoted here are typical, but the values have significant uncertainties and depend on the viewing angle, ejected mass, and explosion energy in ways we quantify. Assessing the overall threat to Earth-like planets, KNe have a similar kill distance to supernovae, but are far less common. However, our results rely on the scant available KN data, and multi-messenger observations will clarify the danger posed by such events.

Astronomy & Astrophysics↗
DOE OSTI2024

Emergent Nucleosynthesis from a 1.2 s Long Simulation of a Black Hole Accretion Disk

Abstract We simulate a black hole accretion disk system with full-transport general relativistic neutrino radiation magnetohydrodynamics for 1.2 s. This system is likely to form after the merger of two compact objects and is thought to be a robust site of r -process nucleosynthesis. We consider the case of a black hole accretion disk arising from the merger of two neutron stars. Our simulation time coincides with the nucleosynthesis timescale of the r -process (∼1 s). Because these simulations are time-consuming, it is common practice to run for a “short” duration of approximately 0.1–0.3 s. We analyze the nucleosynthetic outflow from this system and compare the results of stopping at 0.12 and 1.2 s. We find that the addition of mass ejected in the longer simulation as well as more favorable thermodynamic conditions from emergent viscous ejecta greatly impacts the nucleosynthetic outcome. We quantify the error in nucleosynthetic outcomes between short and long cuts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A Parametric Study of the SASI Comparing General Relativistic and Nonrelativistic Treatments

We present numerical results from a parameter study of the standing accretion shock instability (SASI), investigating the impact of general relativity (GR) on the dynamics. Using GR hydrodynamics with GR gravity, and nonrelativistic (NR) hydrodynamics with Newtonian gravity, in an idealized model setting, we vary the initial radius of the shock, and by varying its mass and radius in concert, the proto-neutron star compactness. We investigate four compactnesses expected in a post-bounce core-collapse supernova (CCSN). We find that GR leads to a longer SASI oscillation period, with ratios between the GR and NR cases as large as 1.29 for the highest-compactness suite. We also find that GR leads to a slower SASI growth rate, with ratios between the GR and NR cases as low as 0.47 for the highest-compactness suite. We discuss implications of our results for CCSN simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Safety First: Stability and Dissipation of Line-tied Force-free Flux Tubes in Magnetized Coronae

Abstract Magnetized plasma columns and extended magnetic structures with both footpoints anchored to a surface layer are an important building block of astrophysical dissipation models. Current loops shining in X-rays during the growth of plasma instabilities are observed in the corona of the Sun and are expected to exist in highly magnetized neutron star magnetospheres and accretion disk coronae. For varying twist and system sizes, we investigate the stability of line-tied force-free flux tubes and the dissipation of twist energy during instabilities using linear analysis and time-dependent force-free electrodynamics simulations. Kink modes ( m = 1) and efficient magnetic energy dissipation develop for plasma safety factors q ≲ 1, where q is the inverse of the number of magnetic field line windings per column length. Higher-order fluting modes ( m > 1) can distort equilibrium flux tubes for q > 1 but induce significantly less dissipation. In our analysis, the characteristic pitch μ ˜ 0 of flux-tube field lines determines the growth rate ( ∝ μ ˜ 0 3 ) and minimum wavelength of the kink instability ( ∝ μ ˜ 0 − 1 ). We use these scalings to determine a minimum flux tube length for the growth of the kink instability for any given μ ˜ 0 . By drawing analogies to idealized magnetar magnetospheres with varying regimes of boundary shearing rates, we discuss the expected impact of the pitch-dependent growth rates for magnetospheric dissipation in magnetar conditions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

On a Spectral Method for β -particle Bound Excitation Collisions in Kilonovae

The interaction of β-particles with the weakly ionized plasma background is an important mechanism for powering the kilonova (KN) transient signal from neutron star mergers. For this purpose, we present an implementation of the approximate fast-particle collision kernel, described by Inokuti following the seminal formulation of Bethe, in a spectral solver of the Vlasov–Maxwell–Boltzmann equation. In particular, we expand the fast-particle plane-wave atomic excitation kernel into coefficients of the Hermite basis, and derive the relevant discrete spectral system. In this fast-particle limit, the approach permits the direct use of atomic data, including optical oscillator strengths, normally applied to photon–matter interaction. The resulting spectral matrix is implemented in the MASS-APP spectral solver framework, in a way that avoids full matrix storage per spatial zone. We numerically verify aspects of the matrix construction, and present a proof-of-principle 3D simulation of a 2D axisymmetric KN ejecta snapshot. Our preliminary numerical results indicate that a reasonable choice of Hermite basis parameters for β-particles in the KN is a bulk velocity parameter u = 0, a thermal velocity parameter α = 0.5c, and a 9 × 9 × 9 mode velocity basis set (Hermite orders of 0–8 in each dimension). For interior-ejecta sample zones, we estimate that the ratio of thermalization from large-angle (≳2fdg5) bound excitation scattering to total thermalization is ~0.002–0.003.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Force-free Wave Interaction in Magnetar Magnetospheres: Computational Modeling in Axisymmetry

Abstract Crustal quakes of highly magnetized neutron stars can disrupt their magnetospheres, triggering energetic phenomena like X-ray and fast radio bursts. Understanding plasma wave dynamics in these extreme environments is vital for predicting energy transport across scales to the radiation length. This study models relativistic plasma wave interaction in magnetar magnetospheres with force-free electrodynamics simulations. For propagation along curved magnetic field lines, we observe the continuous conversion of Alfvén waves to fast magnetosonic (FMS) waves. The conversion efficiency can be up to three times higher when counter-propagating Alfvén waves interact in the equatorial region. Alfvén waves generate FMS waves of twice their frequency during their first crossing of the magnetosphere. After the initial transient burst of FMS waves, Alfvén waves convert to FMS waves periodically, generating variations on timescales of the magnetospheric Alfvén wave crossing time. This decaying FMS wave tail carries a significant portion (half) of the total energy emitted. Plastic damping of “bouncing” Alfvén waves by the magnetar crust has minimal impact on the FMS efficiency. We discuss the implications of the identified wave phenomena for magnetar observations. Outgoing FMS waves can develop electric zones, potential sources of coherent radiation. Long wavelength FMS waves could generate FRBs through reconnection beyond the light cylinder.

Astronomy & Astrophysics↗