Hadron-quark combustion in dense matter is a central topic in the study of phases in compact stars and their high-energy astrophysics. We critically review the literature on hadron-quark combustion, dividing them into a "first wave" that treats the problem as a steady-state burning with or without constraints of mechanical equilibrium, and a "second wave" which uses numerical techniques to resolve the burning front and solves the underlying Partial Differential Equations for the chemistry of the burning front under less restrictive conditions. We detail the inaccuracies that the second wave amends over the first wave, and highlight crucial differences between various approaches in the second wave. We also include results from time-dependent simulations of the reaction zone that include a hadronic EOS, neutrinos, and self-consistent thermodynamics without using parameterized shortcuts.

The hadron-quark combustion front is a system that couples various processes, such as chemical reactions, hydrodynamics, diffusion, and neutrino transport. Previous numerical work has shown that this system is very nonlinear, and can be very sensitive to some of these processes. In these proceedings, we contextualize the hadron-quark combustion as a nonlinear system, subject to dramatic feedback triggered by leptonic weak decays and neutrino transport.

We simulate and study the microphysics of combustion (flame burning) of two flavored quark matter (u, d) to three flavoured quark matter (u,d,s) in a trapped neutrino regime applicable to conditions prevailing in a hot proto-neutron star. The reaction-diffusion-advection equations for (u,d) to (u,d,s) combustion are coupled with neutrino transport, which is modelled through a flux-limited diffusion scheme. The flame speed is proportional to initial lepton fraction because of the release of electron chemical potential as heat, and reaches a steady-state burning speed of (0.001-0.008)c. We find that the burning speed is ultimately driven by the neutrino pressure gradient, given that the pressure gradient induced by quarks is opposed by the pressure gradients induced by electrons. This suggests, somewhat counter-intuitively, that the pressure gradients that drive the interface are controlled primarily by leptonic weak decays rather than by the quark Equation of State (EOS). In other words, the effects of the leptonic weak interaction, including the corresponding weak decay rates and the EOS of electrons and neutrinos, are at least as important as the uncertainties related to the EOS of high density matter. We find that for baryon number densities nB <= 0.35 fm-3, strong pressure gradients induced by leptonic weak decays drastically slow down the burning speed, which is thereafter controlled by the much slower burning process driven by backflowing downstream matter. We discuss the implications of our findings to proto-neutron stars.

FRBs are puzzling, millisecond, energetic radio transients with no discernible source; observations show no counterparts in other frequency bands. The birth of a quark star from a parent neutron star experiencing a quark nova - previously thought undetectable when born in isolation - provides a natural explanation for the emission characteristics of FRBs. The generation of unstable r-process elements in the quark nova ejecta provides millisecond exponential injection of electrons into the surrounding strong magnetic field at the parent neutron star's light cylinder via $\beta$-decay. This radio synchrotron emission has a total duration of hundreds of milliseconds and matches the observed spectrum while reducing the inferred dispersion measure by approximately 200 cm$^{-3}$ pc. The model allows indirect measurement of neutron star magnetic fields and periods in addition to providing astronomical measurements of $\beta$-decay chains of unstable neutron rich nuclei. Using this model, we can calculate expected FRB average energies ($\sim$ 10$^{41}$ ergs), spectra shapes and provide a theoretical framework for determining distances.

Burn-UD is a hydrodynamic combustion code used to model the phase transition of hadronic to quark matter with particular application to the interior of neutron stars. Burn-UD models the flame micro-physics for different equations of state (EoS) on both sides of the interface, i.e. for both the ash (up-down-strange quark phase) and the fuel (up-down quark phase). It also allows the user to explore strange quark seeding produced by different processes including DM annihilation inside neutron stars. The simulations provide a physical window to diagnose whether the combustion process will simmer quietly and slowly, lead to a transition from deflagration to detonation or a (quark) core-collapse explosion. Such an energetic phase transition (a Quark-Nova) would have consequences in high-energy astrophysics and could aid in our understanding of many still enigmatic astrophysical transients. Furthermore, having a precise understanding of the phase transition dynamics for different EoSs could aid further in constraining the nature of the non-perturbative regimes of QCD in general. We hope that Burn-UD will evolve into a platform/software to be used and shared by the QCD community exploring the phases of Quark Matter and astrophysicists working on Compact Stars.

Current r-process models under-produce A<130 nuclei. P-process models either underproduce A<100 and 150<A<165 p-nuclei, or have other limitations. We argue that the puzzles in these two mechanisms can be solved by heavy-ion spallation. Spallation happens in the explosive phase transition from a neutron star to a quark star - the quark-nova. The quark-nova triggers the r-process which creates abundant A>130 nuclei, and spallation fragments these isotopes into A<130 nuclei and all 35 p-nuclei. Our model is universal in relation to a star's age, metallicity, and chemistry.

Current cc-SN models predict overproduction of 44Ti compared to observations. We present a model for an alternative channel where a cc-SN explosion is followed by a neutron star detonation (Quark Nova or QN), resulting in a spallation reaction of SN ejecta that produces 44Ti. We can achieve a 44Ti production of ~ 10^-4 Msun with our model under the right time delay between the QN and the SN. Our model also produces unique signals not found in standard, cc-SN nucleosynthesis models. Some of these unique signals include a significantly large production of 7Be and 22Na. We discuss some of these signals by analyzing the late time light curve and gamma spectroscopy of our model.

Titanium-rich subluminous supernovae are rare and challenge current SN nucleosynthesis models. We present a model in which ejecta from a standard Supernova is impacted by a second explosion of the neutron star (a Quark-nova), resulting in spallation reactions that lead to 56Ni destruction and 44Ti creation under the right conditions. Basic calculations of the spallation products shows that a delay between the two explosions of ~ 5 days reproduces the observed abundance of 44Ti in Cas A and explains its low luminosity as a result of the destruction of 56Ni. Our results could have important implications for lightcurves of subluminous as well as superluminous supernovae.