We present a systematic study of magnetised neutron star head on collisions. We investigate the resulting magnetic field geometries as the two neutron stars merge. Furthermore, we analyze the luminosity produced in these collisions and monitor the evolution of the magnetic fields from the time of merger until the subsequent production of a black hole. At the time of black hole formation the luminosity peaks and rings-down following the decay of the electromagnetic fields. A comparison is presented for three different cases, one where the initial magnetic field in both neutron stars is aligned, one where they are anti-aligned and also one case where they initially have unequal magnetic field strength. We identify regions and set limits so that pair creation and magnetic reconnection would occur in this scenario, and further discuss limits and differences in the radiated energy. This study should be regarded as a toy model of the case where the remnant, of a binary neutron star merger, undergoes a prompt collapse to a black hole with a negligible surrounding disk. We note that the generated electromagnetic pulses resembles the fast radio bursts phenomenology. We consider implications on the high mass mergers leading to a fast prompt collapse to a black hole and the expected flux to be observed at a distance similar to the binary neutron star gravitational wave detection GW190425.

The observations of GW170817/GRB170817A have confirmed that the coalescence of a neutron-star binary is the progenitor of a short gamma-ray burst. In the standard picture of a short gamma-ray burst, a collimated highly relativistic outflow is launched after merger and it successfully breaks out from the surrounding ejected matter. Using initial conditions inspired from numerical-relativity binary neutron-star merger simulations, we have performed general-relativistic hydrodynamic (HD) and magnetohydrodynamic (MHD) simulations in which the jet is launched and propagates self-consistently. The complete set of simulations suggests that: (i) MHD jets have an intrinsic energy and velocity polar structure with a ``hollow core'' subtending an angle $\theta_{\rm core}\approx4^{\circ}-5^{\circ}$ and an opening angle of $\theta_{\rm jet}\gtrsim10^{\circ}$; (ii) MHD jets eject significant amounts of matter and two orders of magnitude more than HD jets; (iii) the energy stratification in MHD jets naturally yields the power-law energy scaling $E(>\Gamma\beta)\propto(\Gamma\beta)^{-4.5}$; (iv) MHD jets provide fits to the afterglow data from GRB170817A that are comparatively better than those of the HD jets and without free parameters; (v) finally, both of the best-fit HD/MHD models suggest an observation angle $\theta_{\rm obs} \simeq 21^{\circ}$ for GRB170817A.

One of the main dissipation processes acting on all scales in relativistic jets is thought to be governed by magnetic reconnection. Such dissipation processes have been studied in idealized environments, such as reconnection layers, which evolve in merging islands and lead to the production of plasmoids, ultimately resulting in efficient particle acceleration. In accretion flows onto black holes, reconnection layers can be developed and destroyed rapidly during the turbulent evolution of the flow. We present a series of two-dimensional general-relativistic magnetohydrodynamic simulations of tori accreting onto rotating black holes focusing our attention on the formation and evolution of current sheets. Initially, the tori are endowed with a poloidal magnetic field having a multi-loop structure along the radial direction and with an alternating polarity. During reconnection processes, plasmoids and plasmoid chains are developed leading to a flaring activity and hence to a variable electromagnetic luminosity. We describe the methods developed to track automatically the plasmoids that are generated and ejected during the simulation, contrasting the behaviour of multi-loop initial data with that encountered in typical simulations of accreting black holes having initial dipolar field composed of one loop only. Finally, we discuss the implications that our results have on the variability to be expected in accreting supermassive black holes.

Binary neutron-star mergers have long been associated with short-duration gamma-ray bursts (GRBs). This connection was confirmed with the first coincident detection of gravitational waves together with electromagnetic radiation from GW170817. The basic paradigm for short-duration GRBs includes an ultra-relativistic jet, but the low-luminosity prompt emission together with follow-up radio and X-ray observations have hinted that this picture may be different in the case of GW170817. In particular, it has been proposed that large amounts of the magnetic energy that is amplified after the merger, can be released when the remnant collapses to a black hole, giving rise to a quasi-spherical explosion impacting on the merger ejecta. Through numerical simulations we investigate this scenario for a range of viewing angles, injected energies and matter densities at the time of the collapse. Depending on the magnitude of the energy injection and the remnant density, we find two types of outflows: one with a narrow relativistic core and one with a wide-angle, but mildly relativistic outflow. Furthermore, very wide outflows are possible, but require energy releases in excess of 10^52 erg.

Recent advances in black hole astrophysics, particularly the first visual evidence of a supermassive black hole at the center of the galaxy M87 by the Event Horizon Telescope (EHT), and the detection of an orbiting "hot spot" nearby the event horizon of Sgr A* in the Galactic center by the Gravity Collaboration, require the development of novel numerical methods to understand the underlying plasma microphysics. Non-thermal emission related to such hot spots is conjectured to originate from plasmoids that form due to magnetic reconnection in thin current layers in the innermost accretion zone. Resistivity plays a crucial role in current sheet formation, magnetic reconnection, and plasmoid growth in black hole accretion disks and jets. We included resistivity in the three-dimensional general-relativistic magnetohydrodynamics (GRMHD) code BHAC and present the implementation of an Implicit-Explicit scheme to treat the stiff resistive source terms of the GRMHD equations. The algorithm is tested in combination with adaptive mesh refinement to resolve the resistive scales and a constrained transport method to keep the magnetic field solenoidal. Several novel methods for primitive variable recovery, a key part in relativistic magnetohydrodynamics codes, are presented and compared for accuracy, robustness, and efficiency. We propose a new inversion strategy that allows for resistive-GRMHD simulations of low gas-to-magnetic pressure ratio and highly magnetized regimes as applicable for black hole accretion disks, jets, and neutron star magnetospheres. We apply the new scheme to study the effect of resistivity on accreting black holes, accounting for dissipative effects as reconnection.

The main hard pulse of prompt gamma-ray emission in GRB$\,$170817A had a duration of $\sim0.5\,{\rm s}$ and its onset was delayed with respect to the gravitational-wave chirp signal by $t_{\rm del} \approx 1.74\,{\rm s}$. Detailed follow-up of the subsequent broadband kilonova emission revealed a two-component ejecta -- a lanthanide-poor ejecta with mass $M_{\rm ej,blue}\approx0.025\,M_\odot$ that powered the early but rapidly fading blue emission and a lanthanide-rich ejecta with mass $M_{\rm ej,red}\approx 0.04\,M_\odot$ that powered the longer lasting redder emission. Both the prompt gamma-ray onset delay and the existence of the blue ejecta with modest electron fraction, $0.2\lesssim Y_e\lesssim0.3$, can be explained if the collapse to a black hole was delayed by the formation of a hypermassive neutron star (HMNS). Here, we determine the survival time of the merger remnant by combining two different constraints, namely, the time needed to produce the requisite blue-ejecta mass and that necessary for the relativistic jet to bore its way out of the expanding ejecta. In this way, we determine that the remnant of GW170817 must have collapsed to a black hole after $t_{\rm coll}=0.98_{-0.26}^{+0.31}\,{\rm s}$. We also discuss how future detections and the delays between the gravitational and electromagnetic emissions can be used to constrain the properties of the merged object.

High-resolution Very-Long-Baseline Interferometry observations of NGC 1052 show a two sided jet with several regions of enhanced emission and a clear emission gap between the two jets.This gap shrinks with increasing frequency and vanishes around $\nu\sim43$ GHz. The observed structures are due to both the macroscopic fluid dynamics interacting with the surrounding ambient medium including an obscuring torus and the radiation microphysics. In this paper we investigate the possible physical conditions in relativistic jets of NGC 1052 by directly modelling the observed emission and spectra via state-of-the-art special-relativistic hydrodynamic (SRHD) simulations and radiative transfer calculations. To investigate the physical conditions in the relativistic jet we coupled our radiative transfer code to evolutionary algorithms and performed simultaneous modelling of the observed jet structure and the broadband radio spectrum. During the calculation of the radiation we consider both thermal and non-thermal emission. In order to compare our model to VLBI observations we take into account the sparse sampling of the u-v plane, the array properties and the imaging algorithm. We present for the first time an end-to-end pipeline for fitting numerical simulations to VLBI observations of relativistic jets taking into account the macrophysics including fluid dynamics and ambient medium configurations together with thermal/non-thermal emission and the properties of the observing array. The detailed analysis of our simulations shows that the structure and properties of the observed relativistic jets in NGC 1052 can be reconstructed by a slightly over-pressured jet ($d_k\sim1.5$) embedded in a decreasing pressure ambient medium

We examine the range of rest-mass densities, temperatures and magnetic fields involved in simulations of binary neutron-star mergers and identify the conditions under which the ideal-magnetohydrodynamics approximation breaks down and hence the magnetic-field decay should be accounted for. We use recent calculations of the conductivities of warm correlated plasma in envelopes of compact stars and find that the magnetic-field decay timescales are much larger than the characteristic timescales of the merger process for lengthscales down to a meter. Because these are smaller than the currently available resolution in numerical simulations, the ideal-magnetohydrodynamics approximation is effectively valid for all realistic simulations. At the same time, we find that the Hall effect can be important at low densities and low temperatures, where it can induce a non-dissipative rearrangement of the magnetic field. Finally, we mark the region in temperature and density where the hydrodynamic description breaks down.

It has been suggested that a non-repeating fast radio burst (FRB) represents the final signal of a magnetized neutron star collapsing to a black hole. In this model, a supramassive neutron star supported by rapid rotation, will collapse to a black hole several thousand to million years after its birth as a result of spin down. The collapse violently snaps the magnetic-field lines anchored on the stellar surface, thus producing an electromagnetic pulse that will propagate outwards and accelerate electrons producing a massive radio burst, i.e. a "blitzar". We present a systematic study of the gravitational collapse of rotating and magnetised neutron stars with special attention to far-field evolution at late times after the collapse. By considering a series of neutron stars with rotation ranging from zero to millisecond periods and different magnetic-field strengths, we show that the blitzar emission is very robust and always characterised by a series sub-millisecond pulses decaying exponentially in amplitude. The luminosity and energy released when the magnetosphere is destroyed are well reproduced by a simple expression in terms of the stellar magnetic field and radius. Finally, we assess the occurrence of pair production during a blitzar scenario, concluding that for typical magnetic-field strengths of $10^{12}\,{\rm G}$ and spin frequencies of a few Hz, pair production is suppressed. Overall, the very good match between the results of the simulations and the luminosities normally observed for FRBs lends credibility to the blitzar model as a simple and yet plausible explanation for the phenomenology of non-repeating FRBs.

In the dawn of the multi-messenger era including gravitational waves, which was marked by the first ever coincident detection of gravitational waves and electromagnetic radiation it is important to lay back and think about established knowledge. Numerical simulations of binary neutron star mergers and simulations of short GRB jets have to combine efforts in order to understand such complicated and phenomenologicaly rich explosions. We review the status of numerical relativity simulations with respect to any jet or magnetized outflow produced after merger. We compare what is known from such simulations, with what is used and obtained from short GRB jet simulations propagating through the BNS ejecta. We point out facts that are established and can be considered known, and things that need to be further revised and/or clarified.

It is known that a binary neutron star merger produces a hypermassive neutron star. The lifetime of this compact remnant depends on the total mass and the equation of state. The collapse of this compact remnant to a black- hole-torus system is expected to give rise to a powerful jet and a short gamma-ray burst. Nevertheless, if the collapse is delayed half a second or so, the surrounding matter would be already accreted and/or expelled and hence no torus will be formed. However, the collapse itself will give rise to a quasi-isotropic magnetized fireball. This magnetic bomb will dissipate much of its energy due to magnetic re-connection and will produce the prompt emission, when the fireball will become transparent to gamma-rays. The energy range of such an explosion depends on the initial magnetic field strength of the two neutron stars and the amplification of the magnetic energy during merger. We discuss the production of a quasi-isotropic magnetized fireball and its subsequent interaction with the ejected matter during merger, as the outcome of the coalescence of a binary neutron star system. We further discuss a possible origin for the duration of the burst and the radial stratification of the flow following the quasi-normal modes of the black hole.

One of the most luminous explosions detected, gamma-ray bursts, especially the so-called long-duration bursts, most probably consist of an intrinsic core-collapse to a black hole inside a super massive star. We point out that this collapse alone will give a generic gravitational wave burst. It has been shown that the strength of this burst depends on the dimensionless spin parameter of the collapsing object. Under descent assumptions the gamma-ray burst's central engine powers the explosion electromagnetically due to the rotation of the newly formed black hole. We argue that the peak luminosity and the isotropic energy of the gamma-ray burst can be associated with the spin of the black hole, due to this mechanism. Since, both gravitational and electromagnetic emission depend on the spin, they can be correlated and thus give a straight estimate for the gravitational wave burst, when we have in hand a gamma-ray burst with known distance. We discuss detectability limits for present and future detectors.

We implement the KORAL code to perform two sets of very long general relativistic radiation magnetohydrodynamic simulations of an axisymmetric optically thin magnetized flow around a non-rotating black hole: one with a new term in the electromagnetic field tensor due to the radiation pressure felt by the plasma electrons on the comoving frame of the electron-proton plasma, and one without. The source of the radiation is the accretion flow itself. Without the new term, the system evolves to a standard accretion flow due to the development of the magneto-rotational instability (MRI). With the new term, however, the system eventually evolves to a magnetically arrested state (MAD) in which a large scale jet-like magnetic field threads the black hole horizon. Our results confirm the secular action of the Cosmic Battery in accretion flows around astrophysical black holes.

We present the first systematic study of the gravitational collapse of rotating and magnetised neutron stars to charged and rotating (Kerr-Newman) black holes. In particular, we consider the collapse of magnetised and rotating neutron stars assuming that no pair-creation takes place and that the charge density in the magnetosphere is so low that the stellar exterior can be described as an electrovacuum. Under these assumptions, which are rather reasonable for a pulsar that has crossed the 'death line', we show that when the star is rotating, it acquires a net initial electrical charge, which is then trapped inside the apparent horizon of the newly formed back hole. We analyse a number of different quantities to validate that the black hole produced is indeed a Kerr-Newman one and show that, in the absence of rotation or magnetic field, the end result of the collapse is a Schwarzschild or Kerr black hole, respectively.

Gamma-ray bursts (GRBs) are believed to be powered by the electromagnetic extraction of spin energy from a black hole endowed with a magnetic field supported by electric currents in a surrounding disk (Blandford & Znajek 1977). A generic feature of this mechanism is that, under certain fairly general assumptions, the energy loss rate decays exponentially. In this work, we are looking precisely for such exponential decay in the lightcurves of long duration GRBs observed with the XRT instrument on the Swift satellite. We found out that almost 30 % of XRT lightcurves show such behavior before they reach the afterglow plateau. According to Blandford & Znajek, the duration of the burst depends on the magnetic flux accumulated on the event horizon. This allows us to estimate the surface magnetic field of a possible progenitor. Our estimations are consistent with magnetic fields observed in Wolf-Rayet stars.

Gamma-ray bursts (GRBs) are violent explosions, coming from cosmological distances. They are detected in gamma-rays (also X-rays, UV, optical, radio) almost every day, and have typical durations of a few seconds to a few minutes. Some GRBs have been reported with extraordinary duration of 10^4 sec, the so-called Ultra Long GRBs. It has been debated whether these form a new distinct class of events or whether they are similar to long GRBs. According to Blandford & Znajek (1977), the spin energy of a rotating black hole can be extracted electromagnetically, should the hole be endowed with a magnetic field supported by electric currents in a surrounding disk. We argue that this can be the case for the central engines of GRBs and we show that the duration of the burst depends on the magnetic flux accumulated on the event horizon of the black hole. We thus estimate the surface magnetic field of a possible progenitor star, and we conclude that an Ultra Long GRB may originate from a progenitor star with a relatively low magnetic field.

The aberrated radiation pressure at the inner edge of the accretion disk around an astrophysical black hole imparts a relative azimuthal velocity on the electrons with respect to the ions which gives rise to a ring electric current that generates large scale poloidal magnetic field loops. This is the Cosmic Battery established by Contopoulos and Kazanas in 1998. In the present work we perform realistic numerical simulations of this important astrophysical mechanism in advection-dominated accretion flows-ADAF. We confirm the original prediction that the inner parts of the loops are continuously advected toward the central black hole and contribute to the growth of the large scale magnetic field, whereas the outer parts of the loops are continuously diffusing outward through the turbulent accretion flow. This process of inward advection of the axial field and outward diffusion of the return field proceeds all the way to equipartition, thus generating astrophysically significant magnetic fields on astrophysically relevant timescales. We confirm that there exists a critical value of the magnetic Prandtl number between unity and 10 in the outer disk above which the Cosmic Battery mechanism is suppressed.

According to Blandford & Znajek (1977), energy can be extracted electromagnetically from a rotating black hole, should the latter be endowed with a magnetic field supported by electric currents in a surrounding disk. We show that exact models of black hole magnetospheres produce Poynting flux that decreases almost exponentially with time. We went through the Swift BAT-XRT lightcurves and identified a subclass of GRBs that exhibits a clear exponential decay over more than three orders of magnitude in flux (EDOHS GRBs). We estimate the energy given-off in the X-rays and discuss a possible correlation between the peak brightness of the X-ray prompt emission and its decay time. We investigate a possible application of this result in high redshift Cosmology.

We investigate the structure of the steady-state force-free magnetosphere around a Kerr black hole in various astrophysical settings. The solution Psi(r,theta) depends on the distributions of the magnetic field line angular velocity omega(Psi) and the poloidal electric current I(Psi). These are obtained self-consistently as eigenfunctions that allow the solution to smoothly cross the two singular surfaces of the problem, the Inner Light Surface (ILS) inside the ergosphere, and the Outer Light Surface (OLS), which is the generalization of the pulsar light cylinder. Magnetic field configurations that cross both singular surfaces (e.g. monopole, paraboloidal) are uniquely determined. Configurations that cross only one light surface e.g. the artificial case of a rotating black hole embedded in a vertical magnetic field) are degenerate. We show that, similarly to pulsars, black hole magnetospheres naturally develop an electric current sheet that potentially plays a very important role in the dissipation of black hole rotational energy and in the emission of high-energy radiation.

We explore the analogy between a rotating magnetized black hole and an axisymmetric pulsar and derive its electromagnetic spindown after its formation in the core collapse of a supermassive star. The spindown shows two characteristic phases, an early Blandford-Znajek phase that lasts a few hundred seconds, and a late pulsar-like afterglow phase that lasts much longer. During the first phase, the spindown luminosity decreases almost exponentially, whereas during the afterglow phase it decreases as t^{-a} with 1<a<1,5. We associate our findings with long duration gamma-ray bursts (GRB) and compare with observations.