Researcher profile

Boaz Katz

Boaz Katz contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
3topics
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

5 published item(s)

preprint2021arXiv

Polarization signatures of the head-on collision model for Type Ia supernovae: How much asymmetry is too much?

In a previous paper, we showed that the asymmetric ejecta produced by (zero impact parameter) head-on collisions of carbon-oxygen white dwarfs allow these progenitor models for Type Ia supernovae (SNe Ia) to cover the observed two-dimensional (2D) distribution of Si II line depths (Branch plot). In this paper, we study the polarization signature associated with the 2D asymmetric ejecta of the collision model and a double-detonation model using similar TARDIS radiative transfer simulations along different lines of sight with a spherical photosphere, combined with a new 3D Monte Carlo polarization code. We show that the polarization $Q$ can be parametrized as a product $Q=Q_{\max}Q_{\rm{x}}$ of a radial structure component $Q_{\max}$ which is insensitive to the model specifics and is shown to be universally around $Q_{\max}\sim 5\%$, and a cancellation component $Q_{\rm{x}}$ which depends on the asymmetry details. The continuum polarization is found to be low for both the collision and double-detonation models with $Q\sim 0.5\%$. However, the irregular Si distribution in the 2D head-on collision model results in Si II line polarization reaching $Q\sim 3\%$ ($Q_{\rm{x}} \lesssim 50\%$) in tension with observations (mostly $\lesssim 1.2\%$). In contrast, we show that the double-detonation model also covers the Branch plot, and yet results in low line polarization $Q\lesssim 0.7\%$ ($Q_{\rm{x}} \sim 10\%$) consistent with previous results and most SNe Ia. These results strengthen the case for asymmetric explosions as progenitors of SNe Ia, emphasizing an additional requirement for large polarization cancellations to account for the low observed line polarizations.

preprint2020arXiv

A Simple Random-Walk Model Explains the Disruption Process of Hierarchical, Eccentric 3-Body Systems

We study the disruption process of hierarchical 3-body systems with bodies of comparable mass. Such systems have long survival times that vary by orders of magnitude depending on the initial conditions. By comparing with 3-body numerical integrations, we show that the evolution and disruption of such systems can be statistically described as a simple random-walk process in the outer-orbit's energy, where the energy-exchange per pericenter passage (step-size) is calculated from the initial conditions. In our derivation of the step-size, we use previous analytic results for parabolic encounters, and average over the (Kozai-Lidov) oscillations in orbital parameters, which are faster then the energy diffusion timescale. While similar random-walk models were studied before, this work differs in two manners: (a) this is the first time that the Kozai-Lidov averaged step-size is derived from first principles and demonstrated to reproduce the statistical evolution of numerical ensembles without fitting parameters, and (b) it provides a characteristic life-time, instead of answering the binary question (stable/unstable), set by case-specific criteria.

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

An asymmetric explosion mechanism may explain the diversity of Si II line widths in Type Ia supernovae

Near maximum brightness, the spectra of Type Ia supernovae (SNe Ia) present typical absorption features of Silicon II observed at roughly 6100A and 5750A. The 2-D distribution of the pseudo-equivalent widths (pEWs) of these features is a useful tool for classifying SNe Ia spectra (Branch plot). Comparing the observed distribution of SNe on the Branch plot to results of simulated explosion models, we find that 1-D models fail to cover most of the distribution. In contrast, we find that TARDIS radiative transfer simulations of the WD head-on collision models along different lines of sight almost fully cover the distribution. We use several simplified approaches to explain this result. We perform order-of-magnitude analysis and model the opacity of the Si lines using LTE and NLTE approximations. Introducing a simple toy model of spectral feature formation, we show that the pEW is a good tracer for the extent of the absorption region in the ejecta. Using radiative transfer simulations of synthetic SNe ejecta, we reproduce the observed Branch plot distribution by varying the luminosity of the SN and the Si density profile of the ejecta. We deduce that the success of the collision model in covering the Branch plot is a result of its asymmetry, which allows for a significant range of Si density profiles along different viewing angles, uncorrelated with a range of $^{56}$Ni yields that cover the observed range of SNe Ia luminosity. We use our results to explain the shape and boundaries of the Branch plot distribution.

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.