QCD phase transition drives supernova explosion of a very massive star

The nature of core-collapse supernova (SN) explosions is yet incompletely understood. The present article revisits the scenario in which the release of latent heat due to a first-order phase transition, from normal nuclear matter to the quark-gluon plasma, liberates the necessary energy to explain observed SN explosions. Here, the role of the metallicity of the stellar progenitor is investigated, comparing a solar metallicity and a low-metallicity case, both having a zero-age main sequence (ZAMS) mass of 75 $M_\odot$. It is found that low-metallicity models belong exclusively to the failed SN branch, featuring the formation of black holes without explosions. It excludes this class of massive star explosions as possible site for the nucleosynthesis of heavy elements at extremely low metallicity, usually associated with the early universe.


Introduction
All stars more massive than about 8 times the mass of the sun (M ⊙ ) end their life as a core-collapse SN. These events belong to the most powerful outbursts in the modern universe, where most of the energy is carried in form of neutrinos. SN explosions are associated with the revival of the stalled bounce shock. The latter forms when the collapsing stellar of a massive star reaches central densities slightly in excess of normal nuclear matter density (ρ sat ≃ 2.5 × 10 14 g cm −3 , equivalent 0.15 fm −3 in nuclear units), when the short-range repulsive nuclear force balances gravity such that the collapsing core bounces back. The extreme conditions reached at the SN interior, supersaturation densities and temperatures on the order of 10-50 MeV as well as large neutron excess given by the electron abundance of Y e = 0.01 − 0.3, make the core-collapse SN ideal highenergy nuclear and particle physics laboratories [1].
Several SN explosions scenarios, i.e. the liberation of energy from the nascent compact central hot and dense object-the proto-neutron star (PNS)-to the bounce shock have been proposed; the magneto-rotational mechanism [2] and the presently considered standard neutrinoheating mechanism [3]. The latter, however, requires multi-dimensional simulations where the neutrino heating efficiency increases in the presence of convection and hydrodynamics instabilities, to obtain massive star explosions for iron-core progenitors more massive than about 9 M ⊙ [4,5]. The exceptions are the slightly less massive progenitors in the zero-age main sequence mass range of 8-9 M ⊙ [6,7,8], for which neutrino-driven SN explosions can be obtained in spherical symmetric simulations.
The present article focuses on an additional SN explosion mechanism. It is related to the energy, more precisely latent heat release, due to a first-order phase transition from normal nuclear matter to the quarkgluon plasma at high density [9,10,11]. First principle calculations of Quantum Chromodynamics (QCD), the theory of strong interactions with quarks and gluons as the degrees of freedom, are only possible at vanishing baryon densities, where the QCD equations are solved numerically being mapped onto a space-time lattice. These computations predict a smooth cross-over transition at a pseudo-critical temperature in the range of 150-160 MeV [12,13,14]. At finite and high baryon den-sity, phenomenological quark matter models have commonly been employed in astrophysical studies of compact stellar objects.
The classical but simplistic quark-bag model EOS [15] has been the foundation in the core-collapse SN study of Ref. [11]. It is based on the equation of state (EOS) of non-interacting quarks, shifted by a constant bag pressure. One of the most striking problems for this class of quark-matter hybrid EOS in applications to astrophysical studies, featuring hot and dense matter, is the difficulty to obtain maximum neutron star masses compatible with the current observational constraint of about 2 M ⊙ [16,17,18,19]. This caveat has been overcome recently, with the introduction of repulsive quark matter interactions [20,21,22], providing sufficient stiffness with increasing density. The extension to finite temperatures and arbitrary isospin asymmetry is seemingly trivial for the class of quark-bag model EOS [23]. It represents a challenge for microscopic quark-matter EOS, such as the Nambu-Jona-Lasinio models [24,25,26] as well as Schwinger-Dyson models [27,28,29,30,31,32]. These phenomenological approaches are applicable to the non-perturbative regime of QCD, however, they are built on a certain truncation of QCD. Most importantly, they are lacking confinement. On the other hand, perturbative QCD is valid only near the limit of asymptotic freedom, where quarks are no more strongly coupled [33,34,35].
The recently developed phenomenological string-flip quark-matter EOS is the foundation of the present work. It features not only a first-order phase transition from normal nuclear matter, it is also built on a mechanism for (de)confinement [22]. It has been extended to arbitrary isospin asymmetry and finite temperatures [36,37]. With model parameters selected such that the conditions for the hadron-quark phase transition to occur at the interior of massive compact stars, i.e. densities for the onset of quark-matter on the order of 2-5 ×ρ sat , one is led to the conclusions that the associated SN explosion mechanism is likely to operate for the high-mass end of core-collapse SN progenitors, associated with the ZAMS mass of > 30 M ⊙ [36], whereas less massive progenitor stars would be subject to other SN explosion scenarios. Note that also the failed SN branch, in which a black hole forms instead of an explosion, is considered within the class of hadron-quark phase transition SN. It depends on the compactness, more precisely, on the mass enclosed inside the PNS at the onset of the hadron-quark phase transition. If the latter exceeds the maximum mass of the hadron-quark hybrid EOS, then an explosion is likely to fail.
The string-flip class of hadron-quark hybrid EOS has also been the foundation for studies of binary neu-tron star mergers [38,39,40]. One of the key goals has been the identification of a first-order phase transition through possibly observable signatures, in the neutrino signal from SN and gravitational waves [41], the latter also in the context of binary neutron star mergers [74].
One important aspect of the SN explosion scenario driven by the QCE phase transition is the role of metallicity 1 , which has not yet been addressed. It strongly affects the evolution of massive stars through mass loss at the star's surface occurring during the different nuclear burning stages, and hence the core structure at the onset of core collapse. This article reports the SN explosion of a very massive ZAMS progenitor star of 75 M ⊙ of solar metallicity, driven by the hadron-quark phase transition, in comparison with those launched from a 75 M ⊙ of extremely low metallicity.
The manuscript is organised as follows: In Sect. 2 the hadron-quark matter hybrid EOS is revisited. Section 3 reviews briefly the SN model employed in this study together with a discussion of the stellar progenitors. The SN simulations are discussed in Sect. 4, and the manuscript closes with the summary and conclusions in Sect. 5.

Dense matter equation of state with phase transition
For the present study the reference hadronic DD2F EOS is implemented. It is a nuclear relativistic meanfield (RMF) model with density-dependent meson-nucleon couplings [42,43], which has been modified from the original DD2 RMF parametrization in order to be consistent with the flow constraint [44]. As a consequence, DD2F is slightly softer than DD2 at densities in excess of saturation density. At low densities, below the saturation density, and at temperatures below about 15-20 MeV, nuclear clusters are taken into account based on the nuclear statistical equilibrium model of Ref. [45], with several 1000 nuclear species with tabulated and partly calculated nuclear masses. It implements the dissolving of all clusters through a geometric excluded volume approach, which can be adjusted to calculations of in-medium nuclear properties in order to overcome the limitations of the statistical model [46].
The DD2F hadronic EOS has been extended to take into account a phase transition to the quark-gluon plasma. Therefore, the string-flip EOS for quark matter is employed [22]. It is based on the relativistic density functional (RDF) approach to quark matter. The present work uses the parametrization denoted as RDF 1.2 from Table 1 DD2F-RDF 1.2 hybrid EOS properties, onset density for the phase transition (ρ onset ), the density for reaching the pure quark-matter phase (ρ final ), as well as inset neutron star mass (M onset ) and maximum mass (M max ) together with the central density of the maximum mass configuration (ρ central | M max ), for T = 0 and a constant entropy per particle of s = 3 k B . the comprehensive catalog of Ref. [37] 2 . For the phasetransition construction between the DD2F hadronic EOS and the string-flip EOS, Maxwell's condition of phase equilibrium has been applied, i.e kinetic equilibrium for the pressures with respect to the baryon chemical potential, for fixed temperature, T , and charge chemical potential, µ Q . The latter is defined as the difference between proton and neutron chemical potentials, µ p and µ n , µ Q = µ p − µ n in the hadronic phase, and as the difference between up and down chemical potentials, µ u and µ d , µ Q = µ u − µ d in the quark matter phase.
Maxwell's construction results in a jump condition for all thermodynamic quantities, from the pure hadronic phase to the pure quark matter phase. However, for practical purposes in applications of astrophysical simulations, data are needed in the coexistence region in between. Therefore, a quark-volume fraction has been employed, with a linear dependence on the baryon density. Further details can be found in Ref. [37]. The class of string-flip EOS is build on a confining mechanism through divergent particle masses towards vanishing density. It still remains to be shown in how far it matches the (dynamical) restoration of chiral symmetry. In addition, in order to obtain sufficient stiffness at high densities, repulsive interactions are taken into account. The original string-flip EOS includes besides linear terms higher-order repulsive interactions, which are omitted here in the RDF 1.2 version, such that the problem of causality-when the speed of sound approaches the speed of light-does not occur in the density range encountered at the interior of compact stellar objects. The string-flip model has been extended to include an isospin mean field, i.e. the vector-isovector interaction channel, which are comparable to the ρ-meson interactions in nuclear RMF-type models. The coupling can be adjusted to a smooth behavior of the symmetry energy across the phase transition boundary (a detailed discussion of this aspect is given in Sect. II C3. of Ref. [37]).
Selected properties of the resulting DD2F-RMF 1.2 hybrid EOS are listed in Table 1, for zero temperature and finite entropy per particle of s = 3 k B . Particularly interesting is the reduction of the onset density for the hadron-quark phase transition, ρ onset , with increasing entropy, viz. temperature, while the final density for reaching the pure quark matter phase, ρ final , remains nearly constant for the temperatures of relevance here on the order of T = 0 − 50 MeV. These aspects will become relevant in the SN simulations which will be discussed below.
The DD2F-RDF 1.2 hadron-quark hybrid EOS in β-equilibrium is illustrated in Fig. 1a for T = 0 and in Fig. 1b for finite entropy per particle of s = 3 k B , in comparison to the hadronic DD2F reference EOS. Note that the pressure, P , contains all contributions, hadrons and/or quarks as well as electrons, positrons and photons. The last two are added following Ref. [48]. The vertical lines mark ρ onset (solid lines) and ρ final (dash-dotted lines), where it is interesting to note that while for T = 0 the pressure is nearly constant in the hadron-quark mixed region, for finite entropy case there is a finite slope. The latter aspect is attributed to the fact that there are different temperatures for the phase transition, T (ρ onset ), and for reaching the pure quark-matter phase, T (ρ final ). In order to preserve the entropy per particle, the temperature decreases throughout the mixed phase, for entropies less than about s ≃ 6 − 7 k B , while the temperature increases for higher entropies. This feature is already illustrated in Fig. 1 in Ref. [36] (see also the solid green lines in Fig. 7 which will be discussed below).
It is also interesting to note that the aforementioned Maxwell construction for the phase transition results in gradients for various quantities, e.g., the proton fraction, in general the charge fraction across the coexistence region, in order to conserve the charge chemical potential (see the middle panels in Fig. 1). This behavior is imposed by the linear dependence on density in the coexistence region. For a real Maxwell construction, both pressure and charge fraction would jump from ρ onset to ρ final . Note that in quark matter we construct proton and neutron abundances, Y p and Y n , from the abundances of up-and down-quark, Y u and Y d , as follows: showing the pressure P , the proton abundance Y p and the temperature T (only for finite entropy), with respect to the restmass density ρ, for T = 0 in graph (a) and constant entropy per particle of s = 3 k B in graph (b), as well as the corresponding mass-radius relations, comparing the hadronic DD2F reference model (red dached lines) and SF2 (solid gray lines). The thin vertical solid and dash-dotted lines mark the location for the onset of the phase transition and when reaching the pure quark-matter phase, respectively, whereas the horizontal line in the massradius relation graphs marks the corresponding onset mass, together with the present maximum neutron star constraint of 2.08 ± 0.07 M ⊙ from the measurement of the pulsar PSR J0740+6620 [16,19].
The corresponding mass-radius relations are shown the right panels of Fig. 1. For the DD2F-RDF 1.2 parameters, at T = 0, the maximum mass of the hybrid star exceeds the maximum mass of the DD2F reference hadronic EOS; the exact values are given in Table 1. The onset mass for the hybrid branch, i.e. the third family of compact stellar objects, is at M onset = 1.37 M ⊙ . These hydrostatic properties change when considering hot configurations, in particular at entropy of s = 3 k B , which represents a canonical value at the PNS interior.
Note that it has been found that the universal relations [49,50,51] hold only for constant but low values of the entropy [52,53], since the entropy is a conserved quantity. Since the latter is also the case for SN, it is suitable to consider the situation of constant entropy, rather than constant temperature, in order to understand what happens in simulations of core-collapse SN regarding thermal effects. The first observation is that the coexistence region extends due to the different behaviors of hadronic and quark matter EOS with increas- Fig. 2 Comparison of the progenitors with a ZAMS masses of 75 M ⊙ but different metallicities, from Ref. [63], showing the density ρ and entropy per particle s (left panels) and the abundance of helium denoted as He (dotted lines), carbon and oxygen denoted as C+O (thin dash-dotted lines) as well as silicon and sulphur denoted as Si+S (dashed lines) and iron-group nuclei collectively denoted as 'Fe' (thick solid lines), with respect to the mass coordinate. In addition, for the low-metallicity progenitor, 56 Ni is shown as thin solid line.
(a) Solar metallicity Fig. 1b), besides the significantly larger gradients between ρ onset (vertical solid lines) and ρ final (vertical dash-dotted lines). Secondly, the radii of the hydrostatic configurations become significantly larger, on the order of 20-40 km. Secondly, while the maximum mass of the DD2F reference hadronic EOS increases, the maximum mass of the DD2F-RDF 1.2 hybrid configuration reduces. Moreover, an extended region of gravitational instability develops after the onset of quark matter at the transition from the stable hadronic to the stable quark matter branch [54]. This twin phenomenon has long been studied in the context of cold neutron stars, i.e. the existence of two compact stars with the same gravitational mass but different radii and hence belonging to different families of compact stars (cf. Ref. [20] and references therein). It is interesting to note that the twin phenomenon is only present in case of finite, and in particular high temperatures, illustrated by the temperature jump around T ≃ 50 MeV at the phase transition (middle panel in Fig. 1b).

Supernova model and stellar progenitors
The present numerical study of core-collapse SN is based on AGILE-BOLTZTRAN (see Ref. [55] and references therein), which is a general relativistic neutrino-radiation hydrodynamics model in spherical symmetry featuring accu- rate 6-species Boltzmann neutrino transport [56]. The set of non-muonic weak interactions considered here can be found in Table 1 of Ref. [57], with the references given accordingly, and the muonic weak processes are introduced and discussed in Refs. [58,56]. AGILE-BOLTZTRAN has an adaptive baryon mass mesh refinement [59,60]. In the present simulations 207 radial grid points are used, and the neutrino phase space is discretized following Refs. [61,62], with 6 neutrino lateral angle bins, cos θ ∈ {−1, +1}, and 36 neutrino energy bins, from 0.5-300 MeV.
For the present investigation of the phase transition in core-collapse SN as explosion mechanism, one of the yet-incompletely addressed questions is the dependence of the stellar progenitor model. In particular, what is the maximum progenitor mass that the hadron-quark phase transition is able to drive the explosion, and how does it depend on metallicity, denoted as Z? Therefore, the most massive stellar progenitor is selected from the comprehensive stellar evolution series of Ref. [63], with a ZAMS mass of 75 M ⊙ , while a stellar model will ZAMS mass of 50 M ⊙ has been explored previously [64,36].
Radial profiles, as a function of the enclosed mass, for selected quantities are shown in Fig. 2a for the solar metallicity model. The current value of the solar surface metallicity is Z ⊙ = 0.0187 ± 0.0009, [65,66]. This features a rather compact stellar core at the onset of stellar collapse, with an iron-core mass of about 1.5 M ⊙ and with central densities in excess of 5 × 10 9 g cm −3 , as well as a narrow silicon-sulfur shell above the iron core (dashed line in the right panel of Fig. 2a. Typical for solar-metallicity stars is a large mass-loss rate, which results in a sharp transition from the stellar core, i.e. carbon-oxygen shell (dash-dotted line in the right panel of Fig. 2a) and the transition to the heliumrich envelope, with a sharp density drop over many orders of magnitude. The resulting total enclosed mass is about 6.36 M ⊙ at the onset of stellar core collapse [63], i.e. about 68.64 M ⊙ has been lost during the previous nuclear burning stages, mostly during hydrogen and helium burning [63]. The later collision of the SN ejecta with the material ejected pre-collapse may power the brightest SN explosions known as super-luminous SN [36].
The situation is very different for stars with lower metallicity. Figure 2b shows the same radial profiles for the same 75 M ⊙ ZAMS mass but at a metallicity of 10 −4 of the solar value, from the same stellar evolution series [63]. This stellar model has experienced nearly zero mass loss. It it related with the photon opacity, which scales with the number of protons and is hence substantially higher for ions. The lack of efficient mass loss leaves a nearly as massive progenitor star, with a total mass at the onset of stellar collapse of 75 M ⊙ [63], featuring a central density of about 10 9 g cm −3 and a higher central entropy per baryon of about s ≃ 2 k B , in comparison to the solar metalicty progenitor with s ≃ 1.5 k B (see Fig. 2). Note also that the nascent iron core contains nearly 3 M ⊙ (the abundance of 56 Ni is shown as thin solid line in the right panel of Fig. 2b).

Supernova simulations
In this section, results will be discussed from simulations of core-collapse SN, launched from the stellar progenitors introduced in Sect. 3, with a ZAMS mass of 75 M ⊙ of solar metallicity and 10 −4 of the solar metallicity [63]. Note that the SN simulations of the solar metallicity progenitor contains the entire remaining mass of 6.36 M ⊙ . Both SN simulations employ the DD2F-RDF-1.2 hadron-quark hybrid EOS, reviewed in Sect. 2.
Both SN simulations proceed qualitatively similarly through the stellar core collapse, featuring slightly higher central temperatures for the low-metallicity model as well as lower central density, than the solar-metallicity one, as illustrated in graphs (a) and (b) in Fig. 3. However, of particular interest is the post-bounce evolution, during which differences grow substantially between the two models. These, in turn, can be related to the mass accretion rates, shown in Fig 4, sampled at the PNS surface. The mass accretion rate,Ṁ = 4πr 2 ρv, with rest-mass density ρ and velocity v at the PNS surface, is significantly larger for the low-metallicity case as a direct consequence of the higher density of the far more massive iron core and of the extended silicon-sulphur layer above the iron core (see the left panels in Fig 2). The mass enclosed inside the PNS grows very rapidly shortly after core bounce, reaching shortly 2 M ⊙ , as illustrated in Fig. 3c. It leads to a rapid rise of the central density and of the central temperature, and in particular to maximum temperatures at the PNS interior, denoted as T max , in excess of 30-50 MeV. The very massive PNS gives rise to high neutrino luminosities, on the order of L ν ≃ 1 − 2 × 10 53 erg s −1 (see Fig. 5) during the evolution after the ν e -deleptonization burst has been launched, about 50 ms post-bounce. This is nearly one order of magnitude higher than the one obtained in normal SN simulations after the high-density siliconsulfur layer has fallen onto the SN shock (c.f. Refs. [67,68,5,69,70,71,72] and references therein). Exactly this phenomenon has not been the case for the low metallicty model here, for the entire time of the post-bounce evolution the neutrino luminosities are powered by ex-traordinary high mass accretion rates, being approximated as accretion luminosity for ν andν e as follows, The values for the PNS mass, M PNS , is larger than for the solar metallicity case by about one third, and the mass accretion rate is larger my more than a factor of two (see Figs. 3 and 4). The high neutrino luminosities are the reflection of enhanced neutrino heating rates, which, in turn, result in larger neutrinosphere radii and also a systematically larger shock radius, as illustrated in Fig. 3d. Such phenomenon has been reported based on the SN simulation of a zero-metallicity stellar progenitor with ZAMS mass of 70 M ⊙ [73], which, however, eventually lead to a failed SN explosion and the subsequent black hole formation at about 253 ms post bounce for this low-metallicity model. It is associated with the sudden rise of the maximum temperature and central density (see graphs (a) and (b) in Fig. 3) and the abrupt cease of the mass accretion rate (see Fig. 4).
Here, for the low-metallicity 75 M ⊙ progenitor, the high mass accretion rate during the early post-bounce evolution results in a critical behavior when the conditions are reached for the hadron-quark phase transition. These are illustrated in Fig. 3 (b) through thin horizontal lines, showing the critical density for the onset of the phase transition, ρ onset = 2 × ρ sat , for the corresponding temperature, shown in Fig. 3 (a), and a representative proton fraction of Y p = 0.3. Once the phase transition takes place, the central quark-matter volume fraction starts to rise rapidly for this model and the PNS becomes gravitationally unstable due to the reduced pressure gradient in the co-existence region, between the two stable hadronic and quark matter phases, as was discussed in Sect. 2. The PNS collapse proceeds supersonically. It halts when a sufficient mass at the PNS interior is converted into the quark matter phase, where the EOS stiffens again. As a consequence, a second shock wave forms which propagates to increasingly larger radii where it turns into an accretion front. However, before the accretion front can accelerate along the decreasing density gradient towards the PNS surface, the PNS collapses behind this accretion front. This is related to the fact that the enclosed PNS mass of M ≃ 2.6 M ⊙ (see graph (c) of Fig. 3) exceeds the maximum mass of the DD2F-RDF-1.2 hybrid EOS (see Table 1). This phenomenon has already been reported in Ref. [36] for the most compact stellar progenitor with ZAMS mass of 25 M ⊙ explored in that study. The SN post-bounce evolution proceeds differently for the progenitor of solar metallicity, featuring a substantially lower mass accretion during the entire postbounce evolution (solid line in Fig. 4). It results in an extended post-bounce mass accretion phase which lasts for nearly 1.0 second, until the central density reaches a value of about twice saturation density, in comparison to the low metallicity case when the same central density is reached already at about 253 ms post bounce (see graph (b) in Fig. 3). The low mass accretion rate for the solar metallicity case, withṀ ≤ 0.1 M ⊙ s −1 already after about 500 ms post bounce, leads to a long timescale for the PNS mass to grow. In fact the PNS mass never reaches the maximum mass of the DD2F-RDF-1.2 hybrid EOS, it stays around M PNS ≃ 1.8 M ⊙ (see graph (c) in Fig. 3). Moreover, the SN shock contracts below R shock < 100−80 km once the silicon-sulfur layer has been fallen onto the shock, and the neutrino luminosities decrease to values on the order of L νx ≃ 1 × 10 51 erg s −1 (ν x denote collectively all heavy lepton flavors, here ν x = ν µ ) and L (νe)νe ≃ 2 × 10 51 erg s −1 (see Fig. 6).
Subsequently, the conditions for the hadron-quark phase transition are reached only after about 1.795 s post bounce for the solar-metallicity case, corresponding to an onset density for the phase transition of ρ onset = 5.6 × 10 14 g cm −3 , indicated by the thin horizontal line in graph (b) of Fig. 3, for the central temperature of T central ≃ 30 MeV, shown in graph (a). Note that also for this model, maximum temperatures reach as high as T max ≃ 40−50 MeV. However, unlike in the case for the low-metallicity model discussed above, the second shock wave, which is associated with the PNS collapse and transition into the hybrid branch, expands to continuously larger radii along the decreasing density at the PNS surface while the PNS remains in a metastable stable. The enclosed mass is well below the maximum mass of the hybrid EOS. During the PNS collapse, the central density rises slightly above ρ central = 5 × ρ sat , and settles back to a value of about ρ central = 4.2×ρ sat once the second shock wave accelerates. The latter quickly reaches the location of the bounce shock, at a radius of about 50 km (see graph (d) in Fig. 3), taking over the bounce shock and expanding continuously to increasingly larger radii. It reaches a radius of 1000 km about 10 ms after its formation (thin red dash-dotted line in graph (d) of Fig. 3). It defines the onset of the SN explosion triggered by the PNS collapse due to the hadron-quark phase transition, indicated by the sudden rise of the diagnostic energy, E diagnostic , shown in graph (e) of Fig. 3, which is the baryon mass integration, da, of the total specific energy, E specific , integrated from the stellar surface, M, towards the center. The specific energy is given as follows, E specific (t, a) = Γ e + 2 Γ + 1 with metric function, Γ = 1 + m/r − u, specific internal energy, e, velocity, u, and the gravitational mass, m. The integration is performed until the mass cut, a cut . The latter defines the mass coordinate inside of which material is not ejected, which is determined at the end of the simulation time of about 12 s post bounce. Note that it includes the neutrino-driven wind phase. Note further that the diagnostic energy contains the gravitational binding energy of the stellar envelop, since the entire progenitor mass is included in the SN simulations.
The diagnostic energy evolves during the SN simulations (see graph (e) in Fig. 3). It is often used as indicator for a successful SN explosion when E diagnostic > 0, referring to the SN explosion energy which is the kinetic energy of the ejecta which should match the asymptotic value obtained for E diagnostic (see the inlay of graph (e) in Fig. 3).
The entire SN evolution for the solar-metallicity model is illustrated in the temperature-density phase diagram, Fig. 7, with the color coding of the electron fraction. The evolution of the central values is given as thick black line. It is interesting to note that the central evolution across the phase co-existence, marked by the blue thick solid (onset density) and thin dash-dotted lines (reaching the pure quark matter phase), proceeds along curves of constant entropy per particle of about s = 1.5 − 3 k B , which indicates the adiabatic collapse of the PNS. As the entropy is conserved, i.e. neutrinos cannot leave the hot and dense PNS interior, the temperature decreases during the phase transition slightly. This phenomenon has already been realized in Ref. [36]. Consequently, when considering the stability of hybrid matter PNS, one should take into account the maximum mass of compact stars at finite entropy per particle (see Table 1). It is interesting to note that the SN evolution in the phase diagram is similar as those of binary neutron star mergers, which consider the hadron quark phase transition [38,74,40], however, reaching somewhat higher temperatures.
The rapid SN shock expansion of the second shock across the neutrinospheres releases an outburst of neutrinos of all flavors, however, dominated byν e and heavy lepton flavor neutrinos (see Fig. 6). This burst lasts for only a few tens of a second. The associated luminosities and average neutrino energies rise to values of several times 10 53 erg s −1 and 40 MeV, respectively. It has been pointed out that this non-standard neutrino signal-not present in neutrino-driven SN simu-lation models-is a smoking-gun signature for a sufficiently strong first-order phase transition. It is observable at the present generation of neutrino detectors, e.g., Super-Kamiokande [36].

Summary and conclusions
The present article reports about results of two corecollapse SN simulations of very massive progenitor stars with the same ZAMS mass of 75 M ⊙ but different metallicities. This important aspect was not explored in previous studies. The metallicity has important consequences for the mass loss during the stellar evolution; high(low) metallicity results in high(low) mass loss rates and consequently in high(low) compactness of the nascent iron core at the onset of stellar core collapse. This, in turn, results in substantial differences during the SN post- Fig. 7 Evolution of the SN simulation of the 75 M ⊙ progenitor, of solar metallicity, in the T −ρ phase diagram where the color-coding is due to the electron fraction, Y e . Indications for the phase transition are set as onset of the quark matter phase, ρ onset , and reaching the pure quark matter phase, ρ final . The evolution of the central temperature and density are also shown (thick solid black line), together with curves of constant entropy per particle for selected values of s = 1 − 7 k B at fixed proton fraction of Y p = 0.3 (solid green lines). bounce evolution, featuring low(high) mass accretion rates. This study considers an EOS with a first-order phase transition from normal nuclear matter to the quark-gluon plasma at high baryon density. In the subsequent SN simulations, there are two competing timescales: (i) the growth of the central density and temperature to reach the critical conditions for the onset of the phase transition, and (ii) the growth of the mass enclosed inside the PNS. Both are determined by the mass accretion rate of the still collapsing stellar core onto the bounce shock and by the compression behavior of the PNS, which, in turn is determined by the high-density behavior of the nuclear EOS. Here, the previous finding are confirmed that a high mass accretion rate may result in a massive PNS at the onset of the hadronquark phase transition, exceeding the maximum mass of the EOS. Only when the timescale for the PNS compression, i.e. the timescale associated with the rising central density and temperature, exceeding the growth timescale of the PNS mass, the nascent PNS will be in a meta-stable state after the phase transition. The formation of a second shock wave determines the onset of the SN explosion. In the present work, this is the case for the progenitor of solar metallicity only. Contrary, for the low-metallicity model, the enclosed mass exceeds the maximum hybrid EOS mass and hence a black hole forms. The SN explosion of the solar-metallicity model confirms the unique observable signature in the neutrino signal.
These findings have important consequences for the dependence on metallicity of this class of massive star explosion, driven by a sufficiently strong hadron-quark phase transition. While the first generation of stars had zero metal content, the metal enrichment is considered to rise continuously during the evolving universe. Such considerations lead to the age-metallicity relation [75]. Here, a metallicity of 10 −4 of the solar value is considered as representative example of extremely metal poor stars [76,77]. Observations of metal poor stars with surface abundances enriched with heavy r-process elements, such as strontium and europium, have been observed frequently [78,79,80,81,82]. It is yet incompletely understood which astrophysical sites can account for these observations, especially at metalicities of Z = 10 −4 × Z ⊙ and below. Specifically, whether binary neutron star mergers can account therefore is presently unclear [83,84,85]. On the other hand, even though SN explosions could have enriched the galaxy at low metallicity, canonical, neutrino-driven core-collapse SN explosions cannot account for the production of heavy neutron-capture elements heavier than molybdenum, with atomic number of Z = 42 [86]. The reason for this is the absence of sufficient neutron rich ejecta. It points to rare and, perhaps speculative sites for the r process associated with massive star explosions, such as magneto-rotational driven jet-like explosions [87,88,89,90]. In addition, SN driven by the hadron-quark phase transition have been suggested as possible but rare rprocess site [91], with an abundance pattern featuring an overproduction of nuclei in the vicinity of the second r-process peak in the region with nuclear mass numbers 140 < A < 180. It was suggested that this site may operate at low metallicity and hence account for the observed metal-poor star data. However, the findings of the present study points to the absence of this class of massive star explosions driven by the hadron-qaurk phase transition at a metallicity of Z ≤ 10 −4 × Z ⊙ .
Zero-and low-metalicty stellar evolution calculations predict qualitatively the same evolutionary pattern for massive stars with ZAMS masses > 30 M ⊙ [63,76,64,77], namely an extended high-density silicon-sulfur layer above the iron core due to negligible mass loss. Consequently, all these massive stars would belong to the failed SN branch with the failed SN shock revival and hence the formation of black holes. It excludes SN explosions of massive stars driven by the hadronquark phase transition as r-process site operating at low metallicity. At which value of the increasing metallicity during the galactic evolution they might have started to contribute remains to be explored. Data Availability Statement This manuscript has associated data in a data repository. The data associated with the equation of state and the supernova simulations are available from the corresponding author on reasonable request.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0.