Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “neutron star mergers”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Impact of Newly Measured 𝛽-Delayed Neutron Emitters around 78 Ni on Light Element Nucleosynthesis in the Neutrino Wind Following a Neutron Star Merger

Neutron emission probabilities and half-lives of 37 𝛽-delayed neutron emitters from 75 Ni to 92 Br were measured at the RIKEN Nishina Center in Japan, including 11 one-neutron and 13 two-neutron emission probabilities and six half-lives for the first time that supersede theoretical estimates. These nuclei lie in the path of the weak 𝑟 process occurring in neutrino-driven winds from the accretion disk formed after the merger of two neutron stars synthesizing elements in the 𝐴∼80 abundance peak. The presence of such elements dominates the accompanying kilonova emission over the first few days and have been identified in the AT2017gfo event, associated to the gravitational wave detection GW170817. Abundance calculations based on over 17,000 simulated trajectories describing the evolution of matter properties in the merger outflows show that the new data lead to an increase of 50%–70% in the abundance of Y, Zr, Nb, and Mo. This enhancement is large compared to the scatter of relative abundances observed in old very metal poor stars and thus is significant in the comparison with other possible astrophysical processes contributing to the light-element production. These results underline the importance of including experimental decay data for very neutron-rich 𝛽 -delayed neutron emitters into 𝑟 -process models.

59 ≤ A ≤ 89↗
DOE OSTI2025

Impact of Magnetic-field-driven Anisotropies on the Equation of State Probed in Neutron Star Mergers

Binary neutron star mergers can produce extreme magnetic fields, some of which can lead to strong magnetar-like remnants. While strong magnetic fields have been shown to affect the dynamics of outflows and angular momentum transport in the remnant, they can also crucially alter the properties of nuclear matter probed in the merger. In this work, we provide a first assessment of the latter, determining the strength of the pressure anisotropy caused by Landau-level quantization and the anomalous magnetic moment. To this end, we perform the first numerical relativity simulation with a magnetic polarization tensor and a magnetic-field-dependent equation of state using a new algorithm we present here, which also incorporates a mean-field dynamo model to control the magnetic field strength present in the merger remnant. Our results show that—in the most optimistic case—corrections to the anisotropy can be in excess of 10% and are potentially largest in the outer layers of the remnant. This work paves the way for a systematic investigation of these effects.

General relativity↗
DOE OSTI2025

Connecting GRBs from binary neutron star mergers to nuclear properties of neutron stars

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.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Nuclear Astrophysics through simulations of neutron star mergers using Monte‐Carlo neutrino radiation transport (DE-SC0020435 Final Technical Report)

Neutron star mergers are an important source of information for nuclear physics. Except for black holes, neutron stars are the densest macroscopic objects known to exist in the Universe. They provide us with a remarkable laboratory to study the poorly understood high-density states of matter, complementing the study of heavy nuclei performed by nuclear physics experiments on Earth. Collisions of neutron stars are of particular interest. They are among the most energetic events observable in the Universe, powering a broad range of signals across the electromagnetic spectrum as well as gravitational wave signals. These signals contain important information about the properties of extremely dense matter. Colliding neutron stars additionally eject large amounts of neutron-rich material into the surrounding interstellar medium -- material that then undergoes rapid neutron-capture (r-process) nucleosynthesis, the mechanism thought to be responsible for the production of about half of the heavy nuclei. In order to understand current and future observations of neutron star mergers, we need reliable models for the signals that they power. An important component in the construction of such models are numerical simulations of colliding neutron stars. These simulations are costly, running for multiple months on supercomputers, and require the inclusion of complex physics (general relativity, magnetohydrodynamics, neutrino physics, nuclear reactions). In this document, we describe the development of new methods for the treatment of neutrinos in merger simulations, as well as the first simulations capable of evolving the equations of neutrino radiation transport directly and their impact on our understanding of neutron star mergers and more broadly in nuclear astrophysics. We also discuss additional work performed in the study of nucleosynthesis and neutrino physics in neutron star mergers as part of Early Career Award DE-SC0020435.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Impact of in situ nuclear networks and atomic opacities on neutron star merger ejecta dynamics, nucleosynthesis, and kilonovae

Context. Binary neutron star merger (BNSM) ejecta are key sites of rapid neutron capture (r-process) nucleosynthesis and they produce kilonovae powered by the radioactive decay of freshly synthesized nuclei. Modeling their evolution requires multi-physics simulations involving hydrodynamics, nuclear reactions, and radiative processes. The impact of nuclear burning and atomic opacity is poorly understood and often treated with simplified prescriptions. Aims. We systematically investigate different treatments of nuclear heating, particle thermalization, and atomic opacities in radiation-hydrodynamics simulations of BNSM ejecta and kilonova light curves. Methods. Ejecta profiles from long-term numerical-relativity simulations of asymmetric neutron star binaries with a massive neutron star remnant were evolved to ∼30 days using a 2D ray-by-ray approach. We compared simplified heating-rate and thermalization prescriptions with in situ Nuclear reaction Network (NN) calculations that track nuclear energy deposition and include a composition-dependent thermalization scheme. We also contrasted various gray opacity models with a frequency-dependent treatment based on atomic calculations. Results. Coupling NN and hydrodynamics significantly affects nucleosynthesis and kilonova emission. Assuming homologous expansion alters abundance evolution and produces a narrower, less populated second r-process peak and a third peak shifted to higher mass numbers. The back-reaction of nuclear heating affects the temperature evolution enough to delay and redden the early (t∼ hours) kilonova peaks. A constant thermalization efficiency underestimates and reddens the early emission while overestimating the late-time luminosity compared to the composition-dependent treatment. Analytical opacity prescriptions yield a more extended, colder photosphere, resulting in dimmer, redder kilonovae at early times (t≲ hour), while the delayed recession of the photosphere prolongs the red emission at t ≳ 5 days. Conclusions. Coupling hydrodynamics to an in situ NN is crucial for reliable nucleosynthesis and kilonova predictions. Resolving the first several hundred milliseconds of the hydrodynamics is essential for robust nucleosynthesis calculations. Composition-dependent thermalization and frequency-dependent, atomic-physics-based opacities are needed to accurately capture the temperature evolution of the ejecta and the brightness and color evolution of the kilonova. Calibrated analytic nuclear-power fits with simplified thermalization and opacity prescriptions can still reproduce the density and temperature evolution of the ejecta.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Signature of hadron-quark crossover in binary-neutron-star mergers

Here, we study observational signatures of the hadron-quark crossover in binary-neutron-star mergers by numerical-relativity simulations with various mass configurations. We employ two equations of state (EoSs) for matter consistent with inference from the observational data. In the crossover scenario the EoS is softened in a density realized in binary-neutron-star mergers and is smoothly continued to quark matter. In the phase transition scenario without crossover, the EoS remains stiff and a first-order phase transition takes place in a density out of reach of mergers. A GW170817-like system forms a remnant massive neutron star in both scenarios, and it collapses into a black hole only in the crossover scenario due to the softening while gravitational-wave emission is strong. This difference is clearly reflected in the sudden shutdown of gravitational waves. For a given EoS, the lifetime of the merger remnant is determined primarily by the total mass of the system. Identifying these features in a variety of future events with the next generation of ground-based gravitational-wave detectors will enable us to clarify details of hadron-quark transition. The mass of the accretion disk surrounding the remnant black hole is affected not only by the lifetime of the remnant but also by the mass ratio of the system. Electromagnetic emission associated with the disk outflow will also be useful for detailed investigation of the hadron-quark transition.

QCD phase transitions↗
DOE OSTI2025

Assessing the difficulty of capturing the distribution function of neutrinos in neutron star merger simulations

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.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Distinguishing prompt-collapse binary neutron star mergers from binary black Holes: Tidal effects and remnant properties

We study the properties of remnants formed in prompt-collapse binary neutron star mergers. We consider nonspinning neutron star binaries over a range of total masses and mass ratios across a set of 22 equations of state, totaling 107 numerical relativity simulations. We report the final mass and spin of the systems (including the accretion disk and ejecta) to be constrained in a narrow range—0.98 ≲ 𝑀 𝑓 /𝑀 ≲ 0.99 for the mass and 0.85 ≲ 𝑎 𝑓 ≲ 0.95 for the dimensionless spin—regardless of the binary configuration and matter effects. This sets them apart from binary black hole merger remnants. We assess the detectability of the postmerger signal in a future 40 km Cosmic Explorer observatory and find that the signal-to-noise ratio in the postmerger of an optimally located and oriented binary at a distance of 100 Mpc can range from <1 to 8, depending on the binary configuration and equation of state, with a majority of them greater than 4 in the set of simulations that we consider. We also consider the distinguishability between prompt-collapse binary neutron star and binary black hole mergers with the same masses and spins. We find that Cosmic Explorer will be able to distinguish such systems primarily via the measurement of tidal effects in the late inspiral. Neutron star binaries with reduced tidal deformability $\tilde{Λ}$ as small as ∼ 3.5 can be identified up to a distance of 100 Mpc, while neutron star binaries with $\tilde{Λ}$ ∼ 22 can be identified to distances greater than 250 Mpc. This is larger than the distance up to which the postmerger will be visible. Finally, we discuss the possible implications of our findings for the equation of state of neutron stars from the gravitational wave event GW230529.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Investigating the impact of higher-order phase transitions in binary neutron-star mergers

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.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Inference of Multichannel r -process Element Enrichment in the Milky Way Using Binary Neutron Star Merger Observations

Observations of GW170817 strongly suggest that binary neutron star (BNS) mergers produce rapid neutron-capture nucleosynthesis ( r -process) elements. However, it remains an open question whether these mergers can account for all the r -process element enrichment in the Milky Way’s history. Here, we constrain the contributions of the BNS channel using astrophysical neutron star observations. The rate and mass distributions are constrained by LIGO/Virgo/Kagra through the latest catalog GWTC-3, the neutron star equation of state by gravitational-wave, radio, and X-ray observations, and the delay time distribution by short gamma-ray burst (GRB) host galaxy associations. We present a Bayesian framework to consistently combine these observations with abundance information to quantify the contribution and uncertainties of single and multiple astrophysical enrichment sources, and obtain a distribution of per-event BNS r -process element yields consistent with geophysical and astrophysical abundance constraints. We then adopt a Galactic chemical evolution model assuming an instantaneous and fixed amount of Fe enrichment from core-collapse supernovae, and show that BNS-only enrichment scenarios remain inconsistent with the observed r-process abundance trend of disk stars in the Galaxy even with the uncertainties in BNS merger observations. Using stellar abundance observations instead of the short GRB constraints, we can infer a shorter BNS delay time distribution with power-law index α ≤ −2.0 and minimum delay time ${t}_{{\rm{\min }}}\leqslant 40$ Myr at 90% confidence, consistent with detailed Galactic chemical evolution models. Such delay times are in tension with those predicted by standard BNS formation models. Alternatively, we confirm that a two-channel scenario, in which the second channel tracks the star formation history without significant delay, can account for both Galactic stellar and short GRB observations. We estimate that 45%–90% of the r -process abundance in the Milky Way today would have been produced by this star formation-tracking channel, rather than BNS mergers with significant delay times.

gravitational wave astronomy↗
DOE OSTI2025

Long-lived neutron-star remnants from asymmetric binary neutron star mergers: element formation, kilonova signals and gravitational waves

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.

abundances↗
DOE OSTI2025

Angle-dependent in Situ Fast Flavor Transformations in Post-neutron-star-merger Disks

Abstract The remnant black hole–accretion disk system resulting from binary neutron star mergers has proven to be a promising site for synthesizing the heaviest elements via rapid neutron capture ( r -process). A critical factor in determining the full r -process pattern in these environments is the neutron richness of the ejecta, which is strongly influenced by neutrino interactions. One key ingredient shaping these interactions is fast neutrino flavor conversions (FFCs), which arise due to angular crossings in neutrino distributions and occur on nanosecond timescales. We present the first three-dimensional in situ angle-dependent modeling of FFCs in postmerger disks, implemented within general relativistic magnetohydrodynamics with Monte Carlo neutrino transport. Our results reveal that, by suppressing electron neutrinos, FFCs more efficiently cool the disk and weaken the early thermally driven wind. Less releptonization due to electron neutrino absorption makes this cooler wind more neutron rich, producing a more robust r -process at higher latitudes of the outflow. This study underscores the necessity of incorporating FFCs in realistic simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Kilonova Emissions from Neutron Star Merger Remnants: Implications for the Nuclear Equation of State

Multimessenger observations of binary neutron star mergers can provide valuable information on the nuclear equation of state (EOS). Here, we investigate the extent to which electromagnetic observations of the associated kilonovae allow us to place constraints on the EOS. For this, we use state-of-the-art three-dimensional general-relativistic magnetohydrodynamics simulations and detailed nucleosynthesis modeling to connect properties of observed light curves to properties of the accretion disk, and hence, the EOS. Using our general approach, we use multimessenger observations of GW170817/AT2017gfo to study the impact of various sources of uncertainty on inferences of the EOS. We constrain the radius of a 1.4M ⊙ neutron star to lie within 10.30 ≤ R 1.4 ≤ 13.0 km and the maximum mass to be M TOV ≤ 3.06M ⊙ .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Impact of neutrino flavor conversions on neutron star merger dynamics, ejecta, nucleosynthesis, and multimessenger signals

We present numerical relativity simulations of binary neutron star mergers incorporating neutrino flavor transformations triggered by fast flavor instability, quantum many-body effects, or potential beyond standard model physics. In both long-lived and short-lived remnant scenarios, neutrino flavor conversions modify species-dependent neutrino luminosities and mean energies, and drive the matter towards more neutron rich conditions. They produce up to more neutron rich ejecta and significantly boost the r-process yields, especially in low-density, near-equatorial outflows. We identify regions unstable to fast flavor instabilities and find that these instabilities persist despite flavor conversions. We further test the sensitivity to the equilibration timescale of the flavor conversions, finding that slower flavor conversions can interact with thermodynamic equilibration, and increase the neutron richness of the ejecta. Flavor conversions may also contribute to stronger gravitational wave and neutrino emissions, pointing to a correlation between neutrino transport and merger dynamics. Furthermore, these results highlight the potential impact of flavor conversions while motivating future work to improve on theoretical understanding of flavor instabilities in global simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Binary neutron star mergers as potential sources for ultrahigh-energy cosmic rays and high-energy neutrinos

Recent studies suggest that the most energetic cosmic rays, exceeding 100 EeV, may primarily consist of 𝑟-process nuclei. This highlights binary neutron star mergers and collapsars as promising sources of ultrahigh-energy cosmic rays (UHECRs). Building on these insights, we examine the conditions that facilitate the efficient production of UHE 𝑟-process nuclei during the prompt radiation (PR), extended emission (EE), and plateau emission phases of short gamma-ray bursts (sGRBs) following neutron star mergers. Our study reveals that jets associated with the PR phase, characterized by typical bulk Lorentz factors (≳ 400–500), dissipation radii, and magnetic field strengths, can accelerate 𝑟-process nuclei to energies ≳ 100 EeV, while preserving them during propagation within the source. Additionally, we investigate the production of HE neutrinos from photomeson and hadronic interactions, as well as from the 𝛽 decay of accelerated 𝑟-process nuclei. We find that the HE neutrino fluxes from sGRBs, mainly produced via photomeson interactions, are significantly limited to preserve the accelerated heavy nuclei, leading to lower fluxes than the predictions without allowing for contributions to UHECRs. Our results suggest that sGRBs may potentially contribute to UHECRs during the PR phase and to HE neutrinos during the EE phase—a scenario that can be tested by future neutrino observatories.

Cosmic ray sources↗
DOE OSTI2025

Open-source library for performance-portable neutrino reaction rates: Application to neutron star mergers

A realistic and detailed description of neutrinos in binary neutron star (BNS) mergers is essential to build reliable models of such systems. To this end, we present bns_nurates, a novel open-source numerical library designed for the efficient on-the-fly computation of neutrino interactions, with particular focus on regimes relevant to BNS mergers. bns_nurates targets a higher level of accuracy and realism in the implementation of commonly employed reactions by accounting for relevant microphysics effects on the interactions, such as weak magnetism and mean field effects. It also includes the contributions of inelastic neutrino scattering off electrons and positrons and (inverse) nucleon decays. Finally, it offers a way to reconstruct the neutrino distribution function in the framework of moment-based transport schemes. As a first application, we compute both energy-dependent and energy-integrated neutrino emissivities and opacities for conditions extracted from a BNS merger simulation with m1 transport scheme. We find some qualitative differences in the results when considering the impact of the additional relevant reactions and of microphysics effects. For example, neutrino-electron/positron scattering reactions are important for the energy exchange of heavy-type neutrinos as they do not undergo semileptonic charged-current processes, when μ± are not accounted for. Moreover, weak magnetism and mean field effects can significantly modify the contribution of β processes for electron-type (anti)neutrinos, increasing at the same time the importance of (inverse) neutron decays. Here, the improved treatment for the reaction rates also modifies the conditions at which neutrinos decouple from matter in the system, potentially affecting their emission spectra.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Effect of magnetic fields on Urca rates in neutron star mergers

Isospin-equilibrating weak processes, called “Urca” processes, are of fundamental importance in astrophysical environments like (proto-)neutron stars, neutron star mergers, and supernovae. In these environments, matter can reach high temperatures of tens of MeVs and be subject to large magnetic fields. We thus investigate Urca rates at different temperatures and field strengths by performing the full temperature and magnetic-fielddependent rate integrals for different equations of state. We find that the magnetic fields play an important role at temperatures of a few MeV, especially close to or below the direct Urca threshold, which is softened by the magnetic field. At higher temperatures, the effect of the magnetic fields can be overshadowed by the thermal effects. Finally, we observe that the magnetic field influences the neutron decay rates more strongly than the electron-capture rates, leading to a shift in the flavor equilibrium.

electroweak interactions in nuclear physics↗
DOE OSTI2025

Cocoon shock breakout emission from binary neutron star mergers

Shock breakout emission is among the first observable signals in a wide variety of astrophysical phenomena, including neutron star (NS) mergers, and it can be the dominant component in low-luminosity short γ-ray bursts, as exemplified by GRB 170817A. In this work, we investigate the cocoon shock breakout emission in NS mergers and how its signal depends on the outermost layers of the ejecta profile, which we derive from general relativistic radiation hydrodynamic simulations. We study the formation of the cocoon as a consequence of a relativistic jet propagating through the ejecta. To explore the influence of the outermost layers of the ejecta on the breakout emission, we explore cases where the ejecta has a sharp cutoff or an extended smooth tail. We find that the shock breakout emission is strongly influenced by the shape of the ejecta outer layers, with extended tails yielding results consistent with the observed properties of GRB 170817A, whereas sharp cutoffs overestimate the radiated energy. Using a Bayesian analysis, we estimate the best fit parameters of the central engine, considering both accreting black hole and magnetized neutron star scenarios. Our findings indicate a slight preference for the scenarios where the engine is a black hole. Furthermore, our work probes the nature of neutron star mergers and highlights the importance of the shape of the ejecta profile in modeling early electromagnetic counterparts to these mergers.

79 ASTRONOMY AND ASTROPHYSICS↗