Researcher profile

Doron Kushnir

Doron Kushnir contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
1topics
3close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

6 published item(s)

preprint2020arXiv

All known Type Ia supernovae models fail to reproduce the observed $t_0-M_\text{Ni56}$ correlation

Type Ia supernovae (SNe Ia) are likely the thermonuclear explosions of carbon-oxygen white-dwarf stars, but their progenitor systems remain elusive. A few theoretical scenarios for the progenitor systems have been suggested, which have been shown to agree with some observational properties of SNe Ia. However, several computational challenges prohibit a robust comparison to the observations. We focus on the observed $t_0-M_\text{Ni56}$ relation, where $t_0$ (the $γ$-rays' escape time from the ejecta) is positively correlated with $M_\text{Ni56}$ (the synthesized $^{56}$Ni mass). Comparing to the $t_0-M_\text{Ni56}$ relation bypasses the need for radiation transfer calculations, as the value of $t_0$ can be directly inferred from the ejecta. We show that all known SNe Ia models fail to reproduce the observed $t_0-M_\text{Ni56}$ correlation.

preprint2020arXiv

An accurate and efficient numerical calculation of detonation waves in multidimensional supernova simulations using a burning limiter and adaptive quasi-statistical equilibrium

Resolving the small length-scale of thermonuclear detonation waves (TNDWs) in supernovae is currently not possible in multidimensional full-star simulations. Additionally, multidimensional simulations usually use small, oversimplistic reaction networks and adopt an ad hoc transition criterion to nuclear statistical equilibrium (NSE). The errors due to the applied approximations are not well understood. We present here a new accurate and efficient numerical scheme that accelerates the calculations by orders of magnitudes and allows the structure of TNDWs to be resolved. The numerical scheme has two important ingredients: (1) a burning limiter that broadens the width of the TNDW while accurately preserving its internal structure, and (2) an adaptive separation of isotopes into groups that are in nuclear statistical quasi-equilibrium, which resolves the time-consuming burning calculation of reactions that are nearly balanced out. Burning is calculated in situ employing the required large networks without the use of post-processing or pre-describing the conditions behind the TNDW. In particular, the approach to and deviation from NSE are calculated self-consistently. The scheme can be easily implemented in multidimensional codes. We test our scheme against accurate solutions of the structure of TNDWs and against homogeneous expansion from NSE. We show that with resolutions that are typical for multidimensional full-star simulations, we reproduce the accurate thermodynamic trajectory (density, temperature, etc.) to an accuracy that is better than a percent for the resolved scales (where the burning limiter is not applied), while keeping the error for unresolved scales (broadened by the burning limiter) within a few percent.

preprint2020arXiv

Constraints on the density distribution of type Ia supernovae ejecta inferred from late-time light-curve flattening

The finite time, $τ_{\rm dep}$, over which positrons from $β^{+}$ decays of $^{56}$Co deposit energy in type Ia supernovae ejecta lead, in case the positrons are trapped, to a slower decay of the bolometric luminosity compared to an exponential decline. Significant light-curve flattening is obtained when the ejecta density drops below the value for which $τ_{\rm dep}$ equals the $^{56}$Co life-time. We provide a simple method to accurately describe this "delayed deposition" effect, which is straightforward to use for analysis of observed light curves. We find that the ejecta heating is dominated by delayed deposition typically from 600 to 1200~day, and only later by longer lived isotopes $^{57}$Co and $^{55}$Fe decay (assuming solar abundance). For the relatively narrow $^{56}$Ni velocity distributions of commonly studied explosion models, the modification of the light curve depends mainly on the $^{56}$Ni mass-weighted average density, $\langle ρ\rangle t^{3}$. Accurate late-time bolometric light curves, which may be obtained with JWST far-infrared (far-IR) measurements, will thus enable to discriminate between explosion models by determining $\langle ρ\rangle t^3$ (and the $^{57}$Co and $^{55}$Fe abundances). The flattening of light curves inferred from recent observations, which is uncertain due to the lack of far-IR data, is readily explained by delayed deposition in models with $\langle ρ\rangle t^{3} \approx 0.2\,M_{\odot}\,(10^{4}\, \textrm{km}\,\textrm{s}^{-1})^{-3}$, and does not imply supersolar $^{57}$Co and $^{55}$Fe abundances.

preprint2020arXiv

The $γ$-ray deposition histories of core-collapse supernovae

The $γ$-ray deposition history in an expanding supernova (SN) ejecta has been mostly used to constrain models for Type Ia SN. Here we expand this methodology to core-collapse SNe, including stripped envelope (SE; Type Ib/Ic/IIb) and Type IIP SNe. We construct bolometric light curves using photometry from the literature and we use the Katz integral to extract the $γ$-ray deposition history. We recover the tight range of $γ$-ray escape times, $t_0\approx30-45\,\textrm{d}$, for Type Ia SNe, and we find a new tight range $t_0\approx80-140\,\textrm{d}$, for SE SNe. Type IIP SNe are clearly separated from other SNe types with $t_0\gtrsim400\,\textrm{d}$, and there is a possible negative correlation between $t_0$ and the synthesized $^{56}$Ni mass. We find that the typical masses of the synthesized $^{56}$Ni in SE SNe are larger than those in Type IIP SNe, in agreement with the results of Kushnir. This disfavours progenitors with the same initial mass range for these explosions. We recover the observed values of $ET$, the time-weighted integrated luminosity from cooling emission, for Type IIP, and we find hints of non-zero $ET$ values in some SE SNe. We apply a simple $ γ$-ray radiation transfer code to calculate the $γ$-ray deposition histories of models from the literature, and we show that the observed histories are a powerful tool for constraining models.

preprint2020arXiv

Towards an accurate description of an accretion induced collapse and the associated ejected mass

We revisit the accretion-induced collapse (AIC) process, in which a white dwarf collapses into a neutron star. We are motivated by the persistent radio source associated with the fast radio burst FRB 121102, which was explained by Waxman as a weak stellar explosion with a small ($\sim 10^{-5}M_{\odot}$) mildly relativistic mass ejection that may be consistent with AIC. Additionally, the interaction of the relatively low ejected mass with a pre-collapse wind might be related to fast optical transients. The AIC is simulated with a one-dimensional, Lagrangian, Newtonian hydrodynamic code. We put an emphasis on accurately treating the equation of state and the nuclear burning, which is required for any study that attempts to accurately simulate AIC. We leave subjects such as neutrino physics and general relativity corrections for future work. Using an existing initial profile and our own initial profiles, we find that the ejected mass is $\sim 10^{-2}$ to $10^{-1}M_{\odot}$ over a wide range of parameters, and we construct a simple model to explain our results.

preprint2019arXiv

Comments on "Numerical Stability of Detonations in White Dwarf Simulations"

Katz & Zingale (2019, KZ19) recently studied a one-dimensional test problem, intended to mimic the process of detonation ignition in head-on collisions of two carbon--oxygen (CO) white dwarfs. They do not obtain ignition of a detonation in pure CO compositions unless the temperature is artificially increased or 5% He is included. In both of these cases they obtain converged ignition only for spatial resolutions better than 0.1 km, which are beyond the capability of multidimensional simulations. This is in a contradiction with the claims of Kushnir et al. (2013, K13), that a convergence to $\sim10\%$ is achieved for a resolution of a few km. Using Eulerian and Lagrangian codes we show that a converged and resolved ignition is obtained for pure CO in this test problem without the need for He or increasing the temperature. The two codes agree to within 1% and convergence is obtained at resolutions of several km. We calculate the case that includes He and obtain a similar slow convergence, but find that it is due to a boundary numerical artifact that can (and should) be avoided. Correcting the boundary conditions allows convergence with resolution of $\sim10\,\textrm{km}$ in an agreement with the claims of K13. It is likely that the slow convergence obtained by KZ19 in this case is because of a similar boundary numerical artifact, but we are unable to verify this. KZ19 further recommended to avoid the use of the burning limiter introduced by K13. We show that their recommendation is not justified.