BZPEERReading, trust and collaboration for researchDiscover
Workspaces
Network
Opportunities
Account
Researcher profile
Luciano Rezzolla contributes to research discovery and scholarly infrastructure.
Trust snapshot
Actions
Identity and collaboration
Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.
Log in to claimDirect collaboration
Claim this author entity first to unlock direct invitations.
Research graph
Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
In this short note, we briefly comment on the analytical bounds that must be imposed on the parameter space of the Rezzolla-Zhidenko (RZ) metric-parametrization approach introduced in Ref. [1]. We hope this will clarify some of the confusion recently emerged on this issue [2].
Imaging a supermassive black hole and extracting physical information requires good knowledge of both the gravitational and the astrophysical conditions near the black hole. When the geometrical properties of the black hole are well understood, extracting information on the emission properties is possible. Similarly, when the emission properties are well understood, extracting information on the black-hole geometry is possible. At present however, uncertainties are present both in the geometry and in the emission, and this inevitably leads to degeneracies in the interpretation of the observations. We explore here the impact of varying geometry and emission coefficient when modelling the imaging of a spherically-accreting black hole. Adopting the Rezzolla-Zhidenko parametric metric to model arbitrary static black-holes, we first demonstrate how shadow-size measurements leave degeneracies in the multidimensional space of metric-deviation parameters, even in the limit of infinite-precision measurements. Then, at finite precision, we show that these degenerate regions can be constrained when multiple pieces of information, such as the shadow-size and the peak image intensity contrast, are combined. Such degeneracies can potentially be eliminated with measurements at increased angular-resolution and flux-sensitivity. While our approach is restricted to spherical symmetry and hence idealised, we expect our results to hold also when more complex geometries and emission processes are considered.
Determining the phase structure of Quantum Chromodynamics (QCD) and its Equation of State (EOS) at densities and temperatures realized inside neutron stars and their mergers is a long-standing open problem. The holographic V-QCD framework provides a model for the EOS of dense and hot QCD, which describes the deconfinement phase transition between a dense baryonic and a quark matter phase. We use this model in fully general relativistic hydrodynamic (GRHD) simulations to study the formation of quark matter and the emitted gravitational wave signal of binary systems that are similar to the first ever observed neutron star merger event GW170817.
The gravitational-wave events GW170817 and GW190425 have led to a number of important insights on the equation of state of dense matter and the properties of neutron stars, such as their radii and the maximum mass. Some of these conclusions have been drawn on the basis of numerical-relativity simulations of binary neutron-star mergers with vanishing initial spins. While this may be a reasonable assumption in equal-mass systems, it may be violated in the presence of large mass asymmetries accompanied by the presence of high spins. To quantify the impact of high spins on multi-messenger gravitational-wave events, we have carried out a series of high-mass binary neutron-star mergers with a highly spinning primary star and large mass asymmetries that have been modelled self-consistently using two temperature-dependent equations of state. We show that, when compared with equal-mass, irrotational binaries, these systems can lead to significant differences in the remnant lifetime, in the dynamical ejecta, in the remnant disc masses, in the secular ejecta, and on the bulk kilonova properties. These differences could be exploited to remove the degeneracy between low- and high-spin priors in the detection of gravitational waves from binary neutron-star mergers.
The maximum mass of a nonrotating neutron star, $M_{\rm TOV}$, plays a very important role in deciphering the structure and composition of neutron stars and in revealing the equation of state (EOS) of nuclear matter. Although with a large-error bar, the recent mass estimate for the black-widow binary pulsar PSR J0952-0607, i.e. $M=2.35\pm0.17~M_\odot$, provides the strongest lower bound on $M_{\rm TOV}$ and suggests that neutron stars with very large masses can in principle be observed. Adopting an agnostic modelling of the EOS, we study the impact that large masses have on the neutron-star properties. In particular, we show that assuming $M_{\rm TOV}\gtrsim 2.35\,M_\odot$ constrains tightly the behaviour of the pressure as a function of the energy density and moves the lower bounds for the stellar radii to values that are significantly larger than those constrained by the NICER measurements, rendering the latter ineffective in constraining the EOS. We also provide updated analytic expressions for the lower bound on the binary tidal deformability in terms of the chirp mass and show how larger bounds on $M_{\rm TOV}$ lead to tighter constraints for this quantity. In addition, we point out a novel quasi-universal relation for the pressure profile inside neutron stars that is only weakly dependent from the EOS and the maximum-mass constraint. Finally, we study how the sound speed and the conformal anomaly are distributed inside neutron stars and show how these quantities depend on the imposed maximum-mass constraints.
Magnetic reconnection is thought to be one of the main energy-dissipation mechanisms fueling energy to the plasma in the vicinity of a black hole. Indeed, plasmoids formed through magnetic reconnection may play a key role in $γ$-ray, X-ray and near-infrared flares from the black hole at the center of our galaxy, SgrA*. We report the results of three-dimensional general-relativistic ideal and resistive magnetohydrodynamics simulations modelling magnetic reconnection in accretion flows around astrophysical black holes. As an important difference with similar works, our accretion discs have an initial dipolar magnetic-field configuration with loops of alternating polarity. We show that current sheets are formed and destroyed rapidly in the turbulent environment of black-hole accretion. Plasmoids are formed from current sheets close to the event horizon, in a region of $\sim2-15$ gravitational radii. We further quantify the magnetic dissipation and the process of energy transfer to the plasmoids, reporting the reconnection rate, the relative current density with respect to the local magnetic field, and the size of the plasmoids. We find that plasmoids gain energy through reconnection and heat up to relativistic temperatures, with the largest ones being sufficiently energetic to leave the black hole near the polar regions. During their evolution, plasmoids are stretched and elongated, becoming disrupted when the shear is sufficiently large, although some plasmoids survive as well-distinguished structures at distances of $\sim30-40$ gravitational radii from the black hole. Finally, we find that in some cases the plasmoids acquire a super-Keplerian azimuthal velocity, as suggested by recent observations of flares from Sgr~A*.
We report measurements of the gravitationally lensed secondary image -- the first in an infinite series of so-called "photon rings" -- around the supermassive black hole M87* via simultaneous modeling and imaging of the 2017 Event Horizon Telescope (EHT) observations. The inferred ring size remains constant across the seven days of the 2017 EHT observing campaign and is consistent with theoretical expectations, providing clear evidence that such measurements probe spacetime and a striking confirmation of the models underlying the first set of EHT results. The residual diffuse emission evolves on timescales comparable to one week. We are able to detect with high significance a southwestern extension consistent with that expected from the base of a jet that is rapidly rotating in the clockwise direction. This result adds further support to the identification of the jet in M87* with a black hole spin-driven outflow, launched via the Blandford-Znajek process. We present three revised estimates for the mass of M87* based on identifying the modeled thin ring component with the bright ringlike features seen in simulated images, one of which is only weakly sensitive to the astrophysics of the emission region. All three estimates agree with each other and previously reported values. Our strongest mass constraint combines information from both the ring and the diffuse emission region, which together imply a mass-to-distance ratio of $4.20^{+0.12}_{-0.06}~μ{\rm as}$ and a corresponding black hole mass of $(7.13\pm0.39)\times10^9M_\odot$, where the error on the latter is now dominated by the systematic uncertainty arising from the uncertain distance to M87*.
Ultra-high angular resolution in astronomy has always been an important vehicle for making fundamental discoveries. Recent results in direct imaging of the vicinity of the supermassive black hole in the nucleus of the radio galaxy M87 by the millimeter VLBI system Event Horizon Telescope and various pioneering results of the Space VLBI mission RadioAstron provided new momentum in high angular resolution astrophysics. In both mentioned cases, the angular resolution reached the values of about 10-20 microrcseconds. Further developments toward at least an order of magnitude "sharper" values are dictated by the needs of astrophysical studies and can only be achieved by placing millimeter and submillimeter wavelength interferometric systems in space. A concept of such the system, called Terahertz Exploration and Zooming-in for Astrophysics (THEZA), has been proposed in the framework of the ESA Call for White Papers for the Voayage 2050 long term plan in 2019. In the current paper we discuss several approaches for addressing technological challenges of the THEZA concept. In particular, we consider a novel configuration of a space-borne millimeter/sub-millimeter antenna which might resolve several bottlenecks in creating large precise mechanical structures. The paper also presents an overview of prospective space-qualified technologies of low-noise analogue front-end instrumentation for millimeter/sub-millimeter telescopes, data handling and processing. The paper briefly discusses approaches to the interferometric baseline state vector determination and synchronisation and heterodyning system. In combination with the original ESA Voyage 2050 White Paper, the current work sharpens the case for the next generation microarcsceond-level imaging instruments and provides starting points for further in-depth technology trade-off studies.
We perform three-dimensional (3D) general-relativistic magnetohydrodynamic simulations to model the jet break-out from the ejecta expected to be produced in a binary neutron-star merger. The structure of the relativistic outflow from the 3D simulation confirms our previous results from 2D simulations, namely, that a relativistic magnetized outflow breaking out from the merger ejecta exhibits a hollow core of $θ_{\rm core}\approx4^{\circ}$, an opening angle of $θ_{\rm jet}\gtrsim10^{\circ}$, and is accompanied by a wind of ejected matter that will contribute to the kilonova emission. We also compute the non-thermal afterglow emission of the relativistic outflow and fit it to the panchromatic afterglow from GRB170817A, together with the superluminal motion reported from VLBI observations. In this way, we deduce an observer angle of $θ_{\rm obs}= 35.7^{\circ \,\,+1.8}_{\phantom{\circ \,\,}-2.2}$. We further compute the afterglow emission from the ejected matter and constrain the parameter space for a scenario in which the matter responsible for the thermal kilonova emission will also lead to a non-thermal emission yet to be observed.
The detection of the binary events GW170817 and GW190814 has provided invaluable constraints on the maximum mass of nonrotating configurations of neutron stars, $M_{_{\rm TOV}}$. However, the large differences in the neutron-star masses measured in GW170817 and GW190814 has also lead to a significant tension between the predictions for such maximum masses, with GW170817 suggesting that $M_{_{\rm TOV}} \lesssim 2.3\,M_{\odot}$, and GW190814 requiring $M_{_{\rm TOV}} \gtrsim 2.5\,M_{\odot}$ if the secondary was a (non- or slowly rotating) neutron star at merger. Using a genetic algorithm, we sample the multidimensional space of parameters spanned by gravitational-wave and astronomical observations associated with GW170817. Consistent with previous estimates, we find that all of the physical quantities are in agreement with the observations if the maximum mass is in the range $M_{_{\rm TOV}} = 2.210^{+0.116}_{-0.123} \,M_{\odot}$ within a $2\textrm{-}σ$ confidence level. By contrast, maximum masses with $M_{_{\rm TOV}} \gtrsim 2.5\,M_{\odot}$, not only require efficiencies in the gravitational-wave emission that are well above the numerical-relativity estimates, but they also lead to a significant under-production of the ejected mass. Hence, the tension can be released by assuming that the secondary in GW190814 was a black hole at merger, although it could have been a rotating neutron star before.
The recent 230 GHz observations of the Event Horizon Telescope (EHT) are able to image the innermost structure of the M87 and show a ring-like structure which is in agreement with thermal synchrotron emission generated in a torus surrounding a supermassive black hole. However, at lower frequencies M87 is characterised by a large-scale and edge-brightened jet with clear signatures of non-thermal emission. In order to bridge the gap between these scales and to provide a theoretical interpretation of these observations we perform general relativistic magnetohydrodynamic simulations of accretion on to black holes and jet launching. M87 has been the target for multiple observations across the entire electromagnetic spectrum. Among these VLBI observations provide unique details on the collimation profile of the jet down to several gravitational radii. In this work we aim to model the observed broad-band spectrum of M87 from the radio to the NIR regime and at the same time fit the jet structure as observed with Global mm-VLBI at 86 GHz. We use general relativistic magnetohydrodynamics and simulate the accretion of the magnetised plasma onto Kerr-black holes in 3D. The radiative signatures of these simulations are computed taking different electron distribution functions into account and a detailed parameter survey is performed in order to match the observations. The results of our simulations show that magnetically arrested disks around fast spinning black holes ($a_\star\geq0.5$) together with a mixture of thermal and non-thermal particle distributions are able to model simultaneously the broad-band spectrum and the innermost jet structure of M87
We investigate the impact of using high-order numerical methods to study the merger of magnetised neutron stars with finite-temperature microphysics and neutrino cooling in full general relativity. By implementing a fourth-order accurate conservative finite-difference scheme we model the inspiral together with the early post-merger and highlight the differences to traditional second-order approaches at the various stages of the simulation. We find that even for finite-temperature equations of state, convergence orders higher than second order can be achieved in the inspiral and post-merger for the gravitational-wave phase. We further demonstrate that the second-order scheme overestimates the amount of proton-rich shock-heated ejecta, which can have an impact on the modelling of the dynamical part of the kilonova emission. Finally, we show that already at low resolution the growth rate of the magnetic energy is consistently resolved by using a fourth-order scheme.
With the goal of providing more accurate and realistic estimates of the secular behavior of the mass accretion and drag rates in the "common-envelope" scenario encountered when a black hole or a neutron star moves in the stellar envelope of a red supergiant star, we have carried out the first general relativistic simulations of the accretion flow onto a nonrotating black hole moving supersonically in a medium with regular but different density gradients. The simulations reveal that the supersonic motion always rapidly reaches a stationary state and it produces a shock cone in the downstream part of the flow. In the absence of density gradients we recover the phenomenology already observed in the well-known Bondi--Hoyle--Lyttleton accretion problem, with super-Eddington mass accretion rate and a shock cone whose axis is stably aligned with the direction of motion. However, as the density gradient is made stronger, the accretion rate also increases and the shock cone is progressively and stably dragged toward the direction of motion. With sufficiently large gradients, the shock-cone axis can become orthogonal to the direction, or even move in the upstream region of the flow in the case of the largest density gradient. Together with the phenomenological aspects of the accretion flow, we have also quantified the rates of accretion of mass and momentum onto the black hole. Simple analytic expressions have been found for the rates of accretion of mass, momentum, drag force, and bremsstrahlung luminosity, all of which have been employed in the astrophysical modelling of the secular evolution of a binary system experiencing a common-envelope evolution...
We report on the first detection of pulsed radio emission from a radio pulsar with the ALMA telescope. The detection was made in the Band-3 frequency range (85-101 GHz) using ALMA in the phased-array mode developed for VLBI observations. A software pipeline has been implemented to enable a regular pulsar observing mode in the future. We describe the pipeline and demonstrate the capability of ALMA to perform pulsar timing and searching. We also measure the flux density and polarization properties of the Vela pulsar (PSR J0835$-$4510) at mm-wavelengths, providing the first polarimetric study of any ordinary pulsar at frequencies above 32 GHz. Finally, we discuss the lessons learned from the Vela observations for future pulsar studies with ALMA, particularly for searches near the supermassive black hole in the Galactic Center, and the potential of using pulsars for polarization calibration of ALMA.
ExaHyPE ("An Exascale Hyperbolic PDE Engine") is a software engine for solving systems of first-order hyperbolic partial differential equations (PDEs). Hyperbolic PDEs are typically derived from the conservation laws of physics and are useful in a wide range of application areas. Applications powered by ExaHyPE can be run on a student's laptop, but are also able to exploit thousands of processor cores on state-of-the-art supercomputers. The engine is able to dynamically increase the accuracy of the simulation using adaptive mesh refinement where required. Due to the robustness and shock capturing abilities of ExaHyPE's numerical methods, users of the engine can simulate linear and non-linear hyperbolic PDEs with very high accuracy. Users can tailor the engine to their particular PDE by specifying evolved quantities, fluxes, and source terms. A complete simulation code for a new hyperbolic PDE can often be realised within a few hours - a task that, traditionally, can take weeks, months, often years for researchers starting from scratch. In this paper, we showcase ExaHyPE's workflow and capabilities through real-world scenarios from our two main application areas: seismology and astrophysics.
The theory of gauge-invariant non-spherical metric perturbations of Schwarzschild black hole spacetimes is now well established. Yet, as different notations and conventions have been used throughout the years, the literature on the subject is often confusing and sometimes confused. The purpose of this paper is to review and collect the relevant expressions related to the Regge-Wheeler and Zerilli equations for the odd and even-parity perturbations of a Schwarzschild spacetime. Special attention is paid to the form they assume in the presence of matter-sources and, for the two most popular conventions in the literature, to the asymptotic expressions and gravitational-wave amplitudes. Besides pointing out some inconsistencies in the literature, the expressions collected here could serve as a quick reference for the calculation of the perturbations of Schwarzschild black hole spacetimes driven by generic sources and for those approaches in which gravitational waves are extracted from numerically generated spacetimes.
The capability of the Event Horizon Telescope (EHT) to image the nearest supermassive black hole candidates at horizon-scale resolutions offers a novel means to study gravity in its strongest regimes and to test different models for these objects. Here, we study the observational appearance at 230 GHz of a surfaceless black hole mimicker, namely a non-rotating boson star, in a scenario consistent with the properties of the accretion flow onto Sgr A*. To this end, we perform general relativistic magnetohydrodynamic simulations followed by general relativistic radiative transfer calculations in the boson star space-time. Synthetic reconstructed images considering realistic astronomical observing conditions show that, despite qualitative similarities, the differences in the appearance of a black hole -- either rotating or not -- and a boson star of the type considered here are large enough to be detectable. These differences arise from dynamical effects directly related to the absence of an event horizon, in particular, the accumulation of matter in the form of a small torus or a spheroidal cloud in the interior of the boson star, and the absence of an evacuated high-magnetization funnel in the polar regions. The mechanism behind these effects is general enough to apply to other horizonless and surfaceless black hole mimickers, strengthening confidence in the ability of the EHT to identify such objects via radio observations.
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 $θ_{\rm core}\approx4^{\circ}-5^{\circ}$ and an opening angle of $θ_{\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(>Γβ)\propto(Γβ)^{-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 $θ_{\rm obs} \simeq 21^{\circ}$ for GRB170817A.
Very Long Baseline Interferometry (VLBI) at sub-millimeter waves has the potential to image the shadow of the black hole in the Galactic Center, Sagittarius A* (Sgr A*), and thereby test basic predictions of the theory of general relativity. We investigate the imaging prospects of a new Space VLBI mission concept. The setup consists of two satellites in polar or equatorial circular Medium-Earth Orbits with slightly different radii, resulting in a dense spiral-shaped uv-coverage with long baselines, allowing for extremely high-resolution and high-fidelity imaging of radio sources. We simulate observations of a general relativistic magnetohydrodynamics model of Sgr A* for this configuration with noise calculated from model system parameters. After gridding the $uv$-plane and averaging visibilities accumulated over multiple months of integration, images of Sgr A* with a resolution of up to 4 $μ$as could be reconstructed, allowing for stronger tests of general relativity and accretion models than with ground-based VLBI.
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.
With the first detection of gravitational waves from a binary system of neutron stars, GW170817, a new window was opened to study the properties of matter at and above nuclear-saturation density. Reaching densities a few times that of nuclear matter and temperatures up to $100\,\rm{MeV}$, such mergers also represent potential sites for a phase transition (PT) from confined hadronic matter to deconfined quark matter. While the lack of a postmerger signal in GW170817 has prevented us from assessing experimentally this scenario, two theoretical studies have explored the postmerger gravitational-wave signatures of PTs in mergers of binary systems of neutron stars. We here extend and complete the picture by presenting a novel signature of the occurrence of a PT. More specifically, using fully general-relativistic hydrodynamic simulations and employing a suitably constructed equation of state that includes a PT, we present the occurrence of a "delayed PT", i.e. a PT that develops only some time after the merger and produces a metastable object with a quark-matter core, i.e. a hypermassive hybrid star. Because in this scenario, the postmerger signal exhibits two distinct fundamental gravitational-wave frequencies -- before and after the PT -- the associated signature promises to be the strongest and cleanest among those considered so far, and one of the best signatures of the production of quark matter in the present Universe.
The first binary neutron-star merger event, GW170817, and its bright electromagnetic counterpart have provided a remarkable amount of information. By contrast, the second event, GW190425, with $M_{\rm tot}=3.4^{+0.3}_{-0.1}\,M_{\odot}$ and the lack of an electromagnetic counterpart, has hardly improved our understanding of neutron-star physics. While GW190425 is compatible with a scenario in which the merger has lead to a prompt collapse to a black hole and little ejected matter to power a counterpart, determining the mass ratio and the effective spin $\tildeχ$ of the binary remains difficult. This is because gravitational waveforms cannot yet well constrain the component spins of the binary. However, since the mass of GW190425 is significantly larger than the maximum mass for nonrotating neutron stars, $M_{_{\rm TOV}}$, the mass ratio $q$ cannot be too small, as the heavier star would not be gravitationally stable. Making use of universal relations and a large number of equations of state, we provide limits in the $(\tildeχ,q)$ plane for GW190425, namely: $q_{\rm min} \geq 0.38$ and $ \tildeχ_{\rm max}\leq 0.20$, assuming $M_\mathrm{tot} \simeq 3.4\, M_\odot$. Finally, we show how future observations of high-mass binaries can provide a lower bound on $M_{_{\rm TOV}}$.
We provide a systematic description of the steps necessary -- and of the potential pitfalls to be encountered -- when implementing a two-moment scheme within an Implicit-Explicit (IMEX) scheme to include radiative-transfer contributions in numerical simulations of general-relativistic (magneto-)hydrodynamics. We make use of the M1 closure, which provides an exact solution for the optically thin and thick limit, and an interpolation between these limits. Special attention is paid to the efficient solution of the emerging set of implicit conservation equations. In particular, we present an efficient method for solving these equations via the inversion of a $4\times 4$-matrix within an IMEX scheme. While this method relies on a few approximations, it offers a very good compromise between accuracy and computational efficiency. After a large number of tests in special relativity, we couple our new radiation code, \texttt{FRAC}, with the general-relativistic magnetohydrodynamics code \texttt{BHAC} to investigate the radiative Michel solution, namely, the problem of spherical accretion onto a black hole in the presence of a radiative field. By performing the most extensive exploration of the parameter space for this problem, we find that the accretion's efficiency can be expressed in terms of physical quantities such as temperature, $T$, luminosity, $L$, and black-hole mass, $M$, via the expression $\varepsilon=(L/L_{\rm Edd})/(\dot{M}/\dot{M}_{\rm Edd})= 7.41\times 10^{-7}\left(T/10^6\,\mathrm{K}\right)^{0.22} \left(L/L_\odot\right)^{0.48} \left(M/M_\odot\right)^{0.48}$, where $L_{\mathrm{Edd}}$ and $\dot{M}_{\mathrm{Edd}}$ are the Eddington luminosity and accretion rate, respectively. Finally, we also consider the accretion problem away from spherical symmetry, finding that the solution is stable under perturbations in the radiation field.
We study in detail the nuclear aspects of a neutron-star merger in which deconfinement to quark matter takes place. For this purpose, we make use of the Chiral Mean Field (CMF) model, an effective relativistic model that includes self-consistent chiral symmetry restoration and deconfinement to quark matter and, for this reason, predicts the existence of different degrees of freedom depending on the local density/chemical potential and temperature. We then use the out-of-chemical-equilibrium finite-temperature CMF equation of state in full general-relativistic simulations to analyze which regions of different QCD phase diagrams are probed and which conditions, such as strangeness and entropy, are generated when a strong first-order phase transition appears. We also investigate the amount of electrons present in different stages of the merger and discuss how far from chemical equilibrium they can be and, finally, draw some comparisons with matter created in supernova explosions and heavy-ion collisions.