We present panchromatic observations and modeling of the Calcium-rich supernova 2019ehk in the star-forming galaxy M100 (d$\approx$16.2 Mpc) starting 10 hours after explosion and continuing for ~300 days. SN 2019ehk shows a double-peaked optical light curve peaking at $t = 3$ and $15$ days. The first peak is coincident with luminous, rapidly decaying $\textit{Swift}$-XRT discovered X-ray emission ($L_x\approx10^{41}~\rm{erg~s^{-1}}$ at 3 days; $L_x \propto t^{-3}$), and a Shane/Kast spectral detection of narrow H$\alpha$ and He II emission lines ($v \approx 500$ km/s) originating from pre-existent circumstellar material. We attribute this phenomenology to radiation from shock interaction with extended, dense material surrounding the progenitor star at $r<10^{15}$ cm and the resulting cooling emission. We calculate a total CSM mass of $\sim$ $7\times10^{-3}$ $\rm{M_{\odot}}$ with particle density $n\approx10^{9}\,\rm{cm^{-3}}$. Radio observations indicate a significantly lower density $n < 10^{4}\,\rm{cm^{-3}}$ at larger radii. The photometric and spectroscopic properties during the second light curve peak are consistent with those of Ca-rich transients (rise-time of $t_r =13.4\pm0.210$ days and a peak B-band magnitude of $M_B =-15.1\pm0.200$ mag). We find that SN 2019ehk synthesized $(3.1\pm0.11)\times10^{-2} ~ \rm{M_{\odot}}$ of ${}^{56}\textrm{Ni}$ and ejected $M_{\rm ej} = (0.72\pm 0.040)~\rm{M_{\odot}}$ total with a kinetic energy $E_{\rm k}=(1.8\pm0.10)\times10^{50}~\rm{erg}$. Finally, deep $\textit{HST}$ pre-explosion imaging at the SN site constrains the parameter space of viable stellar progenitors to massive stars in the lowest mass bin (~10 $\rm{M_{\odot}}$) in binaries that lost most of their He envelope or white dwarfs. The explosion and environment properties of SN 2019ehk further restrict the potential WD progenitor systems to low-mass hybrid HeCO WD + CO WD binaries.

Mergers of neutron stars (NS) and white dwarfs (WD) may give rise to observable explosive transient events. We use 3D hydrodynamical (SPH) simulations, as well as 2D hydrodynamical-thermonuclear simulations (using the FLASH AMR code) to model the disruption of CO-WDs by NSs, which produce faint transient events. We post-process the simulations using a large nuclear network and make use of the SuperNu radiation-transfer code to predict the observational signatures and detailed properties of these transients. We calculate the light-curves (LC) and spectra for five models of NS - CO-WD mergers. The small yields of Ni56 (few x 0.001Msun) results in faint, rapidly-evolving reddened transients (RRTs) with B (R) - peak magnitudes of ~ -12 (-13) to ~ -13 (-15), much shorter and fainter than both regular and faint/peculiar type-Ia SNe. We show that the spectra of RRTs share some similarities with rapidly - evolving transients such as SN2010x, though RRTs are significantly fainter, especially in the I/R bands, and show far stronger Si lines. We estimate that the upcoming Large Synoptic Survey Telescope could detect RRTs at a rate of ~ 10 - 70 yr^-1, through observations in the R/I bands.

Nuclear burning plays a key role in a wide range of astrophysical stellar transients, including thermonuclear, pair instability, and core collapse supernovae, as well as kilonovae and collapsars. Turbulence is now understood to also play a key role in these astrophysical transients. Here we demonstrate that turbulent nuclear burning may lead to large enhancements above the uniform background burning rate, since turbulent dissipation gives rise to temperature fluctuations, and in general the nuclear burning rates are highly sensitive to temperature. We derive results for the turbulent enhancement of the nuclear burning rate under the influence of strong turbulence in the distributed burning regime in homogeneous isotropic turbulence, using probability distribution function (PDF) methods. We demonstrate that the turbulent enhancement obeys a universal scaling law in the limit of weak turbulence. We further demonstrate that, for a wide range of key nuclear reactions, such as C$^{12}$(O$^{16}$, $\alpha$)Mg$^{24}$ and triple-$\alpha$, even relatively modest temperature fluctuations, of the order ten percent, can lead to enhancements of 1 - 3 orders of magnitude in the turbulent nuclear burning rate. We verify the predicted turbulent enhancement directly against numerical simulations, and find very good agreement. We also present an estimation for the onset of turbulent detonation initiation, and discuss implications of our results for the modeling of stellar transients.

Type Ia supernovae (SNe) are thought to originate from the thermonuclear explosions of carbon-oxygen (CO) white dwarfs (WDs). The proposed progenitors of standard type Ia SNe have been studied for decades and can be, generally, divided into explosions of CO WDs accreting material from stellar non-degenerate companions (single-degenerate; SD models), and those arising from the explosive interaction of two CO WDs (double-degenerate; DD models). However, current models for the progenitors of such SNe fail to reproduce the diverse properties of the observed explosions, nor do they explain the inferred rates and the characteristics of the observed populations of type Ia SNe and their expected progenitors. Here we show that the little-studied mergers of CO-WDs with hybrid Helium-CO (He-CO) WDs can provide for a significant fraction of the normal type Ia SNe. Here we use detailed thermonuclear-hydrodynamical and radiative-transfer models to show that a wide range of mergers of CO WDs with hybrid He-CO WDs can give rise to normal type Ia SNe. We find that such He-enriched mergers give rise to explosions for which the synthetic light-curves and spectra resemble those of observed type Ia SNe, and in particular, they can produce a wide range of peak-luminosities, MB(MR)~ 18.4 to 19.2 (~ 18.5 to 19:45), consistent with those observed for normal type Ia SNe. Moreover, our population synthesis models show that, together with the contribution from mergers of massive double CO-WDs (producing the more luminous SNe), they can potentially reproduce the full range of type Ia SNe, their rate and delay-time distribution.

Neutron-star (NS) - white-dwarf (WD) mergers may give rise to observable explosive transients, but have been little explored. We use 2D coupled hydrodynamical-thermonuclear FLASH-code simulations to study the evolution of WD debris-disks formed following WD-disruptions by NSs. We use a 19-elements nuclear-network and a detailed equation-of-state to follow the evolution, complemented by a post-process analysis using a larger 125-isotopes nuclear-network. We consider a wide range of initial conditions and study the dependence of the results on the NS/WD masses ($1.4-2{\rm M_{\odot}}$;$\,{\rm 0.375-0.7\,M_{\odot}}$, respectively), WD-composition (CO/He/hybrid-He-CO) and the accretion-disk structure. We find that viscous inflow in the disk gives rise to continuous wind-outflow of mostly C/O material mixed with nuclear-burning products arising from a weak detonation occurring in the inner-region of the disk. We find that such transients are energetically weak ($10^{48}-10^{49}$ergs) compared with thermonuclear-supernovae (SNe), and are dominated by the (gravitational) accretion-energy. Although thermonuclear-detonations occur robustly in all of our simulations (besides the He-WD) they produce only little energy $(1-10\%$ of the kinetic energy) and $^{56}{\rm Ni}$ ejecta (few$\times10^{-4}-10^{-3}{\rm M_{\odot}})$, with overall low ejecta masses of $\sim0.01-0.1{\rm M_{\odot}}$. Such explosions may produce rapidly-evolving transients, much shorter and fainter than regular type-Ia SNe. The composition and demographics of such SNe appear to be inconsistent with those of Ca-rich type Ib SNe. Though they might be related to the various classes of rapidly evolving SNe observed in recent years, they are likely to be fainter than the typical ones, and may therefore give rise a different class of potentially observable transients.

The mergers of neutron stars (NSs) and white dwarfs (WDs) could give rise to explosive transients, potentially observable with current and future transient surveys. However, the expected properties and distribution of such events is not well understood. Here we characterize the rates of such events, their delay time distribution, their progenitors and the distribution of their properties. We use binary populations synthesis models and consider a wide range of initial conditions and physical processes. In particular we consider different common-envelope evolution models and different NS natal kick distributions. We provide detailed predictions arising from each of the models considered. We find that the majority of NS-WD mergers are born in systems in which mass-transfer played an important role, and the WD formed before the NS. For the majority of the mergers the WDs have a carbon-oxygen composition (60-80%) and most of the rest are with oxygen-neon WDs. The rates of NS-WD mergers are in the range of 3-15% of the type Ia supernovae (SNe) rate. Their delay time distribution is very similar to that of type Ia SNe, but slightly biased towards earlier times. They typically explode in young 0.1-1Gyr environments, but have a tail distribution extending to long, Gyrs-timescales. Models including significant kicks give rise to relatively wide offset distribution extending to hundreds of kpcs. The demographic and physical properties of NS-WD mergers suggest they are likely to be peculiar type Ic-like SNe, mostly exploding in late type galaxies. Their overall properties could be related to a class of rapidly evolving SNe recently observed, while they are less likely to be related to the class of Ca-rich SNe.

White dwarfs (WDs) are the stellar core remnants of low mass stars. They are typically divided into three main composition groups: Oxygen Neon (ONe), Carbon Oxygen (CO) and Helium (He) WDs. The evolution of binary systems can significantly change the evolution of the binary stellar components. In particular, striping the envelope of an evolved star can give rise to a core remnant, which can later evolve into a WD with significantly different composition. Here we focus on the formation and evolution of hybrid HeCO WDs. We follow the formation and stellar evolution of such WDs for a range of initial conditions and provide their detailed structure, mass-radius relation and luminosity temperature evolution. We find that both low mass WDs (< 0.45M , typically thought to be He WDs) and intermediate-mass WDs (0.45 < MWD < 0.7, typically thought to be CO WDs) could in fact be hybrid HeCO WDs, with 5-25 (75 -95)% of their mass in He (CO). We use population synthesis calculations to infer the birth rate and properties of such WDs. We find that hybrid HeCO WD comprise the majority of young (< 2Gyr) WDs in binaries, but are more rare among older WDs in binaries. The high frequency and large He content of such WDs could have an important role in WD WD mergers, and may give rise to sub Chandrasekhar thermonuclear supernova explosions.

A fundamental aspect of the three-body problem is its stability. Most stability studies have focused on the co-planar three-body problem, deriving analytic criteria for the dynamical stability of such pro/retrograde systems. Numerical studies of inclined systems phenomenologically mapped their stability regions, but neither complement it by theoretical framework, nor provided satisfactory fit for their dependence on mutual inclinations. Here we present a novel approach to study the stability of hierarchical three-body systems at arbitrary inclinations, which accounts not only for the instantaneous stability of such systems, but also for the secular stability and evolution through Lidov-Kozai cycles and evection. We generalize the Hill-stability criteria to arbitrarily inclined triple systems, explain the existence of quasi-stable regimes and characterize the inclination dependence of their stability. We complement the analytic treatment with an extensive numerical study, to test our analytic results. We find excellent correspondence up to high inclinations $(\sim120^{\circ}$), beyond which the agreement is marginal. At such high inclinations the stability radius is larger, the ratio between the outer and inner periods becomes comparable, and our secular averaging approach is no longer strictly valid. We therefore combine our analytic results with polynomial fits to the numerical results to obtain a generalized stability formula for triple systems at arbitrary inclinations. Besides providing a generalized secular-based physical explanation for the stability of non co-planar systems, our results have direct implications for any triple systems, and in particular binary planets and moon/satellite systems; we briefly discuss the latter as a test case for our models.