Researcher profile

Riccardo Ciolfi

Riccardo Ciolfi contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2026arXiv

Evolution and afterglow emission of gamma-ray burst jets from binary neutron star mergers

Relativistic jets launched in binary neutron star (BNS) mergers are widely accepted as the engines powering most of the population of short gamma-ray bursts (GRBs). Understanding their structure and dynamics-particularly during and after breakout from the merger ejecta-is crucial for interpreting GRB afterglows, especially for off-axis observers. Traditional models often assume simple angular or radial jet profiles, potentially missing key features emerging for jets piercing through realistic environments. This work aims to investigate the formation and evolution of the jet structure as it propagates through a non-homogeneous, anisotropic BNS merger environment. We focus on how the interaction with the ambient medium shapes the jet's angular and velocity distributions and assess the impact of this realistic structure on the resulting afterglow light curves. We perform a series of 3D relativistic magnetohydrodynamic simulations of jets launched in post-merger environments, exploring different injection conditions. Simulations are evolved to late times, approaching the ballistic regime, where further dynamical evolution becomes negligible. From the resulting outflows, we extract energy and velocity profiles and compute multi-wavelength afterglow light curves using a semi-analytic model that includes radial stratification and the full 3D jet geometry. More energetic or earlier-launched jets drill more efficiently through the ejecta, but all develop asymmetries that leave clear imprints in the off-axis afterglow light curves. All models exhibit a complex multi-shock breakout structure responsible for an early, dimmer peak in the afterglow. Despite structural differences, all simulated jets are consistent with the observational data of the multi-messenger BNS merger event GW170817.

preprint2022arXiv

HARM3D+NUC: A new method for simulating the post-merger phase of binary neutron star mergers with GRMHD, tabulated EOS and neutrino leakage

The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $γ$-ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merger disk. Simulating the disk and predicting the composition of the outflows requires general relativistic magnetohydrodynamical (GRMHD) simulations that include a realistic, finite-temperature equation of state (EOS) and self-consistently calculating the impact of neutrinos. In this work, we detail the implementation of a finite-temperature EOS and the treatment of neutrinos in the GRMHD code HARM3D+NUC, based on HARM3D. We include formal tests of both the finite-temperature EOS and the neutrino leakage scheme. We further test the code by showing that, given conditions similar to those of published remnant disks following neutron star mergers, it reproduces both recombination of free nucleons to a neutron-rich composition and excitation of a thermal wind.

preprint2022arXiv

Implementing a new recovery scheme for primitive variables in the general relativistic magnetohydrodynamic code Spritz

General relativistic magnetohydrodynamic (GRMHD) simulations represent a fundamental tool to probe various underlying mechanisms at play during binary neutron star (BNS) and neutron star (NS) - black hole (BH) mergers. Contemporary flux-conservative GRMHD codes numerically evolve a set of conservative equations based on `conserved' variables which then need to be converted back into the fundamental (`primitive') variables. The corresponding conservative-to-primitive variable recovery procedure, based on root-finding algorithms, constitutes one of the core elements of such GRMHD codes. Recently, a new robust, accurate and efficient recovery scheme called RePrimAnd was introduced, which has demonstrated the ability to always converge to a unique solution. The scheme provides fine-grained error policies to handle invalid states caused by evolution errors, and also provides analytical bounds for the error of all primitive variables. In this work, we describe the technical aspects of implementing the RePrimAnd scheme into the GRMHD code Spritz. To check our implementation as well as to assess the various features of the scheme, we perform a number of GRMHD tests in three dimensions. Our tests, which include critical cases such as a NS collapse to a BH as well as the early evolution (~50 ms) of a Fishbone-Moncrief BH-accrection disk system, show that RePrimAnd is able to support magnetized, low density environments with magnetic-to-fluid pressure ratios as high as 10^4, in situations where the previously used recovery scheme fails.

preprint2021arXiv

Robust Recovery of Primitive Variables in Relativistic Ideal Magnetohydrodynamics

Modern simulation codes for general relativistic ideal magnetohydrodynamics are all facing a long standing technical problem given by the need to recover fundamental variables from those variables that are evolved in time. In the relativistic case, this requires the numerical solution of a system of nonlinear equations. Although several approaches are available, none has proven completely reliable. A recent study comparing different methods showed that all can fail, in particular for the important case of strong magnetization and moderate Lorentz factors. Here, we propose a new robust, efficient, and accurate solution scheme, along with a proof for the existence and uniqueness of a solution, and analytic bounds for the accuracy. Further, the scheme allows us to reliably detect evolution errors leading to unphysical states and automatically applies corrections for typical harmless cases. A reference implementation of the method is made publicly available as a software library. The aim of this library is to improve the reliability of binary neutron star merger simulations, in particular in the investigation of jet formation and magnetically driven winds.

preprint2020arXiv

Binary neutron star mergers after GW170817

The first combined detection of gravitational waves and electromagnetic signals from a binary neutron star (BNS) merger in August 2017 (event named GW170817) represents a major landmark for the ongoing investigation on these extraordinary systems. In this short review, we introduce BNS mergers as events of the utmost importance for astrophysics and fundamental physics and discuss the main discoveries enabled by this first multimessenger observation, which include compelling evidence that such mergers produce a copious amount of heavy r-process elements and can power short gamma-ray bursts. We further discuss key open questions left behind on this event and BNS mergers in general, focussing the attention on the current status and limitations of theoretical models and numerical simulations.

preprint2020arXiv

Collimated outflows from long-lived binary neutron star merger remnants

The connection between short gamma-ray bursts (SGRBs) and binary neutron star (BNS) mergers was recently confirmed by the association of GRB 170817A with the merger event GW170817. However, no conclusive indications were obtained on whether the merger remnant that powered the SGRB jet was an accreting black hole (BH) or a long-lived massive neutron star (NS). Here, we explore the latter case via BNS merger simulations covering up to 250 ms after merger. We report, for the first time in a full merger simulation, the formation of a magnetically-driven collimated outflow along the spin axis of the NS remnant. For the system at hand, the properties of such an outflow are found largely incompatible with a SGRB jet. With due consideration of the limitations and caveats of our present investigation, our results favour a BH origin for GRB 170817A and SGRBs in general. Even though this conclusion needs to be confirmed by exploring a larger variety of physical conditions, we briefly discuss possible consequences of all SGRB jets being powered by accreting BHs.

preprint2020arXiv

Intrinsic properties of the engine and jet that powered the short gamma-ray burst associated with GW170817

GRB 170817A was a subluminous short gamma-ray burst detected about 1.74 s after the gravitational wave signal GW170817 from a binary neutron star (BNS) merger. It is now understood as an off-axis event powered by the cocoon of a relativistic jet pointing 15 to 30 degrees away from the direction of observation. The cocoon was energized by the interaction of the incipient jet with the non-relativistic baryon wind from the merger remnant, resulting in a structured outflow with a narrow core and broad wings. In this paper, we couple the observational constraints on the structured outflow with a model for the jet-wind interaction to constrain the intrinsic properties with which the jet was launched by the central engine, including its time delay from the merger event. Using wind prescriptions inspired by magnetized BNS merger simulations, we find that the jet was launched within about 0.4 s from the merger, implying that the 1.74 s observed delay was dominated by the fireball propagation up to the photospheric radius. We also constrain, for the first time for any gamma-ray burst, the jet opening angle at injection and set a lower limit to its asymptotic Lorentz factor. These findings suggest an initially Poynting-flux dominated jet, launched via electromagnetic processes. If the jet was powered by an accreting black hole, they also provide a significant constraint on the survival time of the metastable neutron star remnant.

preprint2020arXiv

Magnetically driven baryon winds from binary neutron star merger remnants and the blue kilonova of August 2017

The observation of a radioactively-powered kilonova associated with the first binary neutron star (BNS) merger detected in gravitational waves proved that these events are ideal sites for the production of heavy r-process elements. However, the physical origin of the ejected material responsible for the early ("blue") and late ("red") components of this kilonova is still debated. Here, we investigate the possibility that the early/blue kilonova originated from the magnetically driven baryon wind launched after merger by the metastable neutron star remnant. Exploiting a magnetized BNS merger simulation with over 250 ms of post-merger evolution, we can follow for the first time the full mass ejection process up to its final decline. We find that the baryon wind carries $\simeq\!0.010\!-\!0.028\,M_\odot$ of unbound material, proving that the high mass estimated for the blue kilonova can be achieved. We also find expansion velocities of up to $\sim\!0.2\,c$, consistent with the lower end of the observational estimates, and we discuss possible effects neglected here that could further increase the final ejecta velocity. Overall, our results show that the magnetically driven baryon wind represents a viable channel to explain the blue kilonova.

preprint2020arXiv

SN 2017gci: a nearby Type I Superluminous Supernova with a bumpy tail

We present and discuss the optical spectro-photometric observations of the nearby (z=0.087) Type I superluminous supernova (SLSN I) SN 2017gci, whose peak K-corrected absolute magnitude reaches Mg=-21.5 mag. Its photometric and spectroscopic evolution includes features of both slow and of fast evolving SLSN I, thus favoring a continuum distribution between the two SLSN-I subclasses. In particular, similarly to other SLSNe I, the multi-band light curves of SN 2017gci show two rebrightenings at about 103 and 142 days after the maximum light. Interestingly, this broadly agrees with a broad emission feature emerging around 6520 A after 51 days from the maximum light, which is followed by a sharp knee in the light curve. If we interpret this feature as Halpha, this could support the fact that the bumps are the signature of late interactions of the ejecta with a (hydrogen rich) circumstellar material. Then we fitted magnetar and CSM-interaction powered synthetic light curves onto the bolometric one of SN 2017gci. In the magnetar case, the fit suggests a polar magnetic field Bp = 6 x 1e14 G, an initial period of the magnetar Pinitial=2.8 ms, an ejecta mass Mejecta=9 Msun and an ejecta opacity k = 0.08 cm g^{-1} . A CSM interaction scenario would imply a CSM mass of 5 Msun and an ejecta mass of 12 Msun. Finally, the nebular spectrum of phase 187 days was modeled, deriving a mass of 10 Msun for the ejecta. Our models suggest that either a magnetar or CSM interaction might be the power sources for SN 2017gci and that its progenitor was a massive (40 Msun) star.

preprint2020arXiv

Spritz: a new fully general-relativistic magnetohydrodynamic code

The new era of multimessenger astrophysics requires the capability of studying different aspects of the evolution of compact objects. In particular, the merger of neutron star binaries is a strong source of gravitational waves and electromagnetic radiation, from radio to $γ$-rays, as demonstrated by the detection of GW170817 and its electromagnetic counterparts. In order to understand the physical mechanisms involved in such systems, it is necessary to employ fully general relativistic magnetohydrodynamic (GRMHD) simulations able to include the effects of a composition and temperature dependent equation of state describing neutron star matter as well as neutrino emission and reabsorption. Here, we present our new code named $\texttt{Spritz}$ that solves the GRMHD equations in 3D Cartesian coordinates and on a dynamical spacetime. The code can support tabulated equations of state, taking into account finite temperature effects and allowing for the inclusion of neutrino radiation. In this first paper, we present the general features of the code and a series of tests performed in special and general relativity to assess the robustness of the basic GRMHD algorithms implemented. Among these tests, we also present the first comparison between a non-staggered and a staggered formulation of the vector potential evolution, which is used to guarantee the divergence-less character of the magnetic field. With respect to other publicly available GRMHD codes, $\texttt{Spritz}$ combines the robust approach of a staggered formulation of the vector potential together with the use of an equation of state driver ($\texttt{EOS_Omni}$) that can allow the code to use finite temperature equations of state. A next version of the code will fully test the $\texttt{EOS_Omni}$ driver by coupling it with a neutrino leakage scheme.

preprint2020arXiv

The key role of magnetic fields in binary neutron star mergers

The first multimessenger observation of a binary neutron star (BNS) merger in August 2017 demonstrated the huge scientific potential of these extraordinary events. This breakthrough led to a number of discoveries and provided the best evidence that BNS mergers can launch short gamma-ray burst (SGRB) jets and are responsible for a copious production of heavy r-process elements. On the other hand, the details of the merger and post-merger dynamics remain only poorly constrained, leaving behind important open questions. Numerical relativity simulations are a powerful tool to unveil the physical processes at work in a BNS merger and as such they offer the best chance to improve our ability to interpret the corresponding gravitational wave (GW) and electromagnetic emission. Here, we review the current theoretical investigation on BNS mergers based on general relativistic magnetohydrodynamics simulations, paying special attention to the magnetic field as a crucial ingredient. First, we discuss the evolution, amplification, and emerging structure of magnetic fields in BNS mergers. Then, we consider their impact on various critical aspects: (i) jet formation and the connection with SGRBs, (ii) matter ejection, r-process nucleosynthesis, and radiocatively-powered kilonova transients, and (iii) post-merger GW emission.