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.

We present SphericalNR, a new framework for the publicly available Einstein Toolkit that numerically solves the Einstein field equations coupled to the equations of general relativistic magnetohydrodynamics (GRMHD) in a 3+1 split of spacetime in spherical coordinates without symmetry assumptions. The spacetime evolution is performed using reference-metric versions of either the Baumgarte-Shapiro-Shibata-Nakamura equations or the fully covariant and conformal Z4 system with constraint damping. We have developed a reference-metric version of the Valencia formulation of GRMHD with a vector potential method, guaranteeing the absence of magnetic monopoles during the evolution. In our framework, every dynamical field (both spacetime and matter) is evolved using its components in an orthonormal basis with respect to the spherical reference-metric. Furthermore, all geometric information about the spherical coordinate system is encoded in source terms appearing in the evolution equations. This allows for the straightforward extension of Cartesian high-resolution shock-capturing finite volume codes to use spherical coordinates with our framework. To this end, we have adapted GRHydro, a Cartesian finite volume GRMHD code already available in the Einstein Toolkit, to use spherical coordinates. We present the full evolution equations of the framework, as well as details of its implementation in the Einstein Toolkit. We validate SphericalNR by demonstrating it passes a variety of challenging code tests in static and dynamical spacetimes.

We compute the binding energy and angular momentum of a test-particle at the last stable circular orbit (LSO) on the equatorial plane around a general relativistic, rotating neutron star (NS). We present simple, analytic, but accurate formulas for these quantities that fit the numerical results and which can be used in several astrophysical applications. We demonstrate the accuracy of these formulas for three different equations of state (EOS) based on nuclear relativistic mean-field theory models and argue that they should remain still valid for any NS EOS that satisfy current astrophysical constraints. We compare and contrast our numerical results with the corresponding ones for the Kerr metric characterized by the same mass and angular momentum.

The induced gravitational collapse (IGC) paradigm explains a class of energetic, $E_{\rm iso}\gtrsim 10^{52}$~erg, long-duration gamma-ray bursts (GRBs) associated with Ic supernovae, recently named binary-driven hypernovae (BdHNe). The progenitor is a tight binary system formed of a carbon-oxygen (CO) core and a neutron star companion. The supernova ejecta of the exploding CO core triggers a hypercritical accretion process onto the neutron star, which reaches in a few seconds the critical mass, and gravitationally collapses to a black hole emitting a GRB. In our previous simulations of this process we adopted a spherically symmetric approximation to compute the features of the hypercritical accretion process. We here present the first estimates of the angular momentum transported by the supernova ejecta, $L_{\rm acc}$, and perform numerical simulations of the angular momentum transfer to the neutron star during the hyperaccretion process in full general relativity. We show that the neutron star: i) reaches in a few seconds either mass-shedding limit or the secular axisymmetric instability depending on its initial mass; ii) reaches a maximum dimensionless angular momentum value, $[c J/(G M^2)]_{\rm max}\approx 0.7$; iii) can support less angular momentum than the one transported by supernova ejecta, $L_{\rm acc} > J_{\rm NS,max}$, hence there is an angular momentum excess which necessarily leads to jetted emission.

We construct equilibrium configurations of uniformly rotating neutron stars for selected relativistic mean-field nuclear matter equations of state (EOS). We compute in particular the gravitational mass ($M$), equatorial ($R_{\rm eq}$) and polar ($R_{\rm pol}$) radii, eccentricity, angular momentum ($J$), moment of inertia ($I$) and quadrupole moment ($M_2$) of neutron stars stable against mass-shedding and secular axisymmetric instability. By constructing the constant frequency sequence $f=716$ Hz of the fastest observed pulsar, PSR J1748-2446ad, and constraining it to be within the stability region, we obtain a lower mass bound for the pulsar, $M_{\rm min}=[1.2$-$1.4] M_\odot$, for the EOS employed. Moreover we give a fitting formula relating the baryonic mass ($M_b$) and gravitational mass of non-rotating neutron stars, $M_b/M_\odot=M/M_\odot+(13/200)(M/M_\odot)^2$ [or $M/M_\odot=M_b/M_\odot-(1/20)(M_b/M_\odot)^2$], which is independent on the EOS. We also obtain a fitting formula, although not EOS independent, relating the gravitational mass and the angular momentum of neutron stars along the secular axisymmetric instability line for each EOS. We compute the maximum value of the dimensionless angular momentum, $a/M\equiv c J/(G M^2)$ (or "Kerr parameter"), $(a/M)_{\rm max}\approx 0.7$, found to be also independent on the EOS. We compare and contrast then the quadrupole moment of rotating neutron stars with the one predicted by the Kerr exterior solution for the same values of mass and angular momentum. Finally we show that, although the mass quadrupole moment of realistic neutron stars never reaches the Kerr value, the latter is closely approached from above at the maximum mass value, as physically expected from the no-hair theorem. In particular the stiffer the EOS is, the closer the Kerr solution is approached.