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 $\gamma$-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.

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. The expansion velocities of up to 0.2 c are marginally consistent with the observational estimates and we discuss various effects neglected here that could further improve the current agreement. Overall, our results show that the magnetically driven baryon wind is a viable explanation for the blue kilonova.

The recent multimessenger observation of the short gamma-ray burst (SGRB) GRB 170817A together with the gravitational wave (GW) event GW170817 provides evidence for the long-standing hypothesis associating SGRBs with binary neutron star (BNS) mergers. The nature of the remnant object powering the SGRB, which could have been either an accreting black hole (BH) or a long-lived magnetized neutron star (NS), is, however, still uncertain. General relativistic magnetohydrodynamic (GRMHD) simulations of the merger process represent a powerful tool to unravel the jet launching mechanism, but so far most simulations focused the attention on a BH as the central engine, while the long-lived NS scenario remains poorly investigated. Here, we explore the latter by performing a GRMHD BNS merger simulation extending up to ~100 ms after merger, much longer than any previous simulation of this kind. This allows us to (i) study the emerging structure and amplification of the magnetic field and observe a clear saturation at magnetic energy $E_\mathrm{mag} \sim 10^{51}$ erg, (ii) follow the magnetically supported expansion of the outer layers of the remnant NS and its evolution into an ellipsoidal shape without any surrounding torus, and (iii) monitor density, magnetization, and velocity along the axis, observing no signs of jet formation. We also argue that the conditions at the end of the simulation disfavor later jet formation on subsecond timescales if no BH is formed. Furthermore, we examine the rotation profile of the remnant, the conversion of rotational energy associated with differential rotation, the overall energy budget of the system, and the evolution of the GW frequency spectrum. Finally, we perform an additional simulation where we induce the collapse to a BH ~70 ms after merger, in order to gain insights on the prospects for massive accretion tori in case of a late collapse. We find that...