Normalized to: Kushnir, D.

**33** article(s) in total. **33** co-authors, from **1** to **12** common article(s). Median position in authors list is **1,0**.

[1] oai:arXiv.org:1904.08427 [pdf] - 2099184

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

Comments: Accepted for publication in ApJ, 15 pages, 9 figures

Submitted: **2019-04-17**, last modified: **2020-04-28**

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.

[2] oai:arXiv.org:2004.07244 [pdf] - 2078769

The $\gamma$-ray deposition histories of core-collapse supernovae

Comments: 28 pages, 12 figures

Submitted: **2020-04-15**

The $\gamma$-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 curve using
photometry from the literature and we use the Katz integral to extract the
$\gamma$-ray deposition history. We recover the tight range of $\gamma$-ray
escape times, $t_0\approx30-45\,\textrm{day}$, for Type Ia SNe, and we find a
new tight range $t_0\approx80-140\,\textrm{day}$, for SE SNe. Type IIP SNe are
clearly separated from other SNe types with $t_0\gtrsim400\,\textrm{day}$, 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 (2015). This disfavors 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 $ \gamma$-ray radiation
transfer code to calculate the $\gamma$-ray deposition histories of models from
the literature, and we show that the observed histories are a powerful tool for
constraining models.

[3] oai:arXiv.org:2001.10005 [pdf] - 2068199

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

Comments: 9 pages, 5 figures. Minor revisions following referee's report

Submitted: **2020-01-27**, last modified: **2020-03-22**

The finite time, $\tau_{\rm dep}$, over which positrons from $\beta^{+}$
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 $\tau_{\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 \rho \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 \rho \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 \rho\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.

[4] oai:arXiv.org:1912.06151 [pdf] - 2067535

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

Comments: 21 pages, 15 figures. The main results are summarized in the first 5
pages (with 3 figures). Revised following referee's review

Submitted: **2019-12-12**, last modified: **2020-03-20**

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.

[5] oai:arXiv.org:1904.10001 [pdf] - 2025478

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

Comments: 7 pages, 6 figures. Modified following referee report

Submitted: **2019-04-22**, last modified: **2019-12-12**

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.

[6] oai:arXiv.org:1910.07537 [pdf] - 1981815

Analytical calculation of the numerical results of Khatami and Kasen for
transient peak time and luminosity

Comments: 3 pages, 1 figure

Submitted: **2019-10-16**

The diffusion approximation is often used to study supernovae light-curves
around peak light, where it is applicable. By analytic arguments and numerical
studies of toy models, Khatami & Kasen (2019) recently argued for a new
approximate relation between peak bolometric Luminosity, $L_p$, and the time of
peak since explosion, $t_p$, for transients involving homologous expansion:
$L_p=2/(\beta t_p)^2\int_0 ^{\beta t_{p}} t'Q(t')dt'$, where $Q(t)$ is the
heating rate of the ejecta, and $\beta$ is an order unity parameter that is
calibrated from numerical calculations. Khatami & Kasen (2019) demonstrated its
validity using Monte-Carlo radiation transfer simulations of ejecta with
homogenous density and (for most cases considered) constant opacity.
Interestingly, constant values of $\beta$ accurately reproduce the numerical
calculations for different heating distributions and over a wide range of
energy release times. Here we show that the diffusion and the adiabatic loss of
energy in homologous expansion is equivalent to a static diffusion equation and
provide an analytic solution for the case of uniform density and opacity
(extending the results of Pinto & Eastman 2000). Our accurate analytical
solutions reproduce and extend the results of Khatami & Kasen (2019) for this
case, allowing clarification for the universality of their peak time-luminosity
relation as well as new limitations to its use.

[7] oai:arXiv.org:1805.08788 [pdf] - 1863811

Screening of fusion reactions from the principle of detailed balance and
application to the $pep$ reaction

Comments: 4 pages, no figures. Minor revisions following referee report

Submitted: **2018-05-22**, last modified: **2019-04-02**

Dewitt et al. (1973) derived a useful relation between the plasma screening
factor for a reaction of two fusing ions and their chemical potentials, based
on the plasma pair distribution functions. We show that their result can be
derived in a simpler, more straightforward way, by applying the principle of
detailed balance, which also enables us to generalize the relation to reactions
involving $N$ fusing ions. In order to demonstrate the usefulness of applying
the principle of detailed balance, we calculate the screening factor for the
$pep$ reaction, $p+e+p\rightarrow ^{2}$D $+\nu_{e}$. For the plasma conditions
near the centre of the Sun, the reaction is suppressed by roughly the same
amount by which the reaction $p+p\rightarrow ^{2}$D $+$ $e^{+}+\nu_{e}$ is
enhanced. This effect may be measured in the near future.

[8] oai:arXiv.org:1902.01197 [pdf] - 1905676

Late time kilonova light curves and implications to GW 170817

Comments: Improved treatment of synchrotron subtraction

Submitted: **2019-02-04**, last modified: **2019-02-26**

We discuss the late time (tens of days) emission from the radioactive ejecta
of mergers involving neutron stars, when the ionization energy loss time of
beta-decay electrons and positrons exceeds the expansion time. We show that if
the e$^\pm$ are confined to the plasma (by magnetic fields), then the time
dependence of the plasma heating rate, $\dot{\varepsilon}_d$, and hence of the
bolometric luminosity $L=\dot{\varepsilon}_d$, are given by $d\log L/d\log
t\simeq-2.8$, nearly independent of the composition and of the instantaneous
radioactive energy release rate, $\dot{\varepsilon}$. This universality of the
late time behavior is due to the weak dependence of the ionization loss rate on
composition and on e$^\pm$ energy. The late time IR and optical measurements of
GW 170817 are consistent with this expected behavior provided that the
ionization loss time exceeds the expansion time at $t>t_\varepsilon\approx
7$~d, as predicted based on the early (few day) electromagnetic emission.

[9] oai:arXiv.org:1902.02251 [pdf] - 1863976

H-$\alpha$ emission in the nebular spectrum of the Type Ia supernova
ASASSN-18tb

Comments: Submitted to MNRAS

Submitted: **2019-02-06**

As part of the 100IAS survey, a program aimed to obtain nebular-phase spectra
for a volume-limited and homogeneous sample of Type Ia supernovae (SNe Ia), we
observed ASASSN-18tb (SN 2018fhw) at 139 days past maximum light. ASASSN-18tb
was a fast-declining, sub-luminous event that fits well within the observed
photometric and spectroscopic distributions of the SN Ia population. We detect
a prominent H$\alpha$ emission line of $L_{{\rm
H}\alpha}=2.2\pm0.2\times10^{38}$ ergs s$^{-1}$ with FWHM $\approx1100$ km
s$^{-1}$ in the nebular-phase spectrum of this SN Ia. High luminosity H$\alpha$
emission ($L_{{\rm H}\alpha}\gtrsim 10^{40}$ ergs~s$^{-1}$) has previously been
discovered in a rare class of SNe Ia-like objects showing CSM interactions (SNe
Ia-CSM). They predominantly belong to over-luminous ($M_{\rm max}<-19$ mag in
optical) 1991T-like SNe Ia and are exclusively found in star-forming galaxies.
By contrast, ASASSN-18tb is a sub-luminous SN Ia ($M_{B, {\rm max}}\sim -17.7$
mag) found in an early-type galaxy dominated by old stellar populations. We
discuss possible origins for the observed hydrogen. Out of 75 SNe Ia for which
we have so far obtained nebular spectra in 100IAS, no other SN shows a $\sim
1000 \,{\rm km s^{-1}}$ H$\alpha$ emission line with comparable line luminosity
as ASASSN-18tb, emphasizing the rarity of such emission in the nebular phase.
Based on preliminary results from our survey, the rate for ASASSN-18tb-like
nebular H$\alpha$ emission could be as high as $\sim 10\%$ level among
sub-luminous SNe Ia.

[10] oai:arXiv.org:1806.01286 [pdf] - 1795718

The structure of detonation waves in supernovae revisited

Comments: 33 pages, 27 figures, 14 tables. Revised following a referee report

Submitted: **2018-06-04**, last modified: **2018-11-28**

The structure of a thermonuclear detonation wave can be solved accurately
and, thus, may serve as a test bed for studying different approximations that
are included in multidimensional hydrodynamical simulations of supernova. We
present the structure of thermonuclear detonations for the equal mass fraction
of $^{12}$C and $^{16}$O (CO) and for pure $^{4}$He (He) over a wide range of
upstream plasma conditions. The lists of isotopes we constructed allow us to
determine the detonation speeds, as well as the final states for these
detonations, with an uncertainty of the percent level (obtained here for the
first time). We provide our results with a numerical accuracy of $\sim0.1\%$,
which provides an efficient benchmark for future studies. We further show that
CO detonations are pathological for all upstream density values, which differs
from previous studies, which concluded that for low upstream densities CO
detonations are of the Chapman-Jouget (CJ) type. We provide an approximate
condition, independent of reaction rates, that allows to estimate whether
arbitrary upstream values will support a detonation wave of the CJ type. Using
this argument, we are able to show that CO detonations are pathological and to
verify that He detonations are of the CJ type, as was previously claimed for
He. Our analysis of the reactions that control the approach to nuclear
statistical equilibrium, which determines the length-scale of this stage,
reveals that at high densities, the reactions
$^{11}$B$+p\leftrightarrow3^{4}$He plays a significant role, which was
previously unknown.

[11] oai:arXiv.org:1711.09638 [pdf] - 1747841

Constraints on the ejecta of the GW170817 neutron-star merger from its
electromagnetic emission

Comments: Accepted to MNRAS. Expanded explanations, improved figures

Submitted: **2017-11-27**, last modified: **2018-09-04**

We present a simple analytic model, that captures the key features of the
emission of radiation from material ejected by the merger of neutron stars
(NS), and construct the multi-band and bolometric luminosity light curves of
the transient associated with GW170817, AT\,2017gfo, using all available data.
The UV to IR emission is shown to be consistent with a single
$\approx0.05$\,M$_\odot$ component ejecta, with a power-law velocity
distribution between $\approx 0.1\,c$ and $>0.3\,c$, a low opacity,
{$\kappa<1$\,cm$^2$\,g$^{-1}$}, and a radioactive energy release rate
consistent with an initial $Y_{\rm e}<0.4$. The late time spectra require an
opacity of $\kappa_\nu\approx0.1$\,cm$^2$\,g$^{-1}$ at 1 to $2\mu$m. If this
opacity is provided entirely by Lanthanides, their implied mass fraction is
$X_{\rm Ln}\approx10^{-3}$, approximately 30 times below the value required to
account for the solar abundance. The inferred value of $X_{\rm Ln}$ is
uncertain due to uncertainties in the estimates of IR opacities of heavy
elements, which also do not allow the exclusion of a significant contribution
to the opacity by other elements (the existence of a slower ejecta rich in
Lanthanides, that does not contribute significantly to the luminosity, can also
not be ruled out). The existence of a relatively massive, $\approx
0.05$\,M$_\odot$, ejecta with high velocity and low opacity is in tension with
the results of numerical simulations of NS mergers.

[12] oai:arXiv.org:1805.00010 [pdf] - 1960272

A Significantly off-center Ni56 Distribution for the Low-Luminosity Type
Ia Supernova SN 2016brx from the 100IAS survey

Comments: Minor Changes. Accepted by MNRAS Letter

Submitted: **2018-04-30**, last modified: **2018-05-31**

We present nebular-phase spectra of the Type Ia supernova (SN Ia) 2016brx, a
member of the 1991bg-like subclass that lies at the faint end of the SN Ia
luminosity function. Nebular spectra are available for only three other
1991bg-like SNe, and their Co line centers are all within <~ 500 km/s of each
other. In contrast, the nebular Co line center of SN 2016brx is blue-shifted by
>1500 km/s compared to them and by ~1200 km/s compared to the rest frame. This
is a significant shift relative to the narrow nebular line velocity dispersion
of <~ 2000 km/s of these SNe. The large range of nebular line shifts implies
that the Ni56 in the ejecta of SN 1991bg-like events is off-center by ~1000
km/s rather than universally centrally confined as previously suggested. With
the addition of SN 2016brx, the Co nebular line shapes of 1991bg-like objects
appear to connect with the brighter SNe Ia that show double-peak profiles,
hinting at a continuous distribution of line profiles among SNe Ia. One class
of models to produce both off-center and bi-modal Ni56 distributions is
collisions of white dwarfs with unequal and equal masses.

[13] oai:arXiv.org:1210.1574 [pdf] - 1579430

Preliminary evidence for a virial shock around the Coma galaxy cluster

Comments: Replaced with published version

Submitted: **2012-10-04**, last modified: **2017-09-21**

Galaxy clusters, the largest gravitationally bound objects in the Universe,
are thought to grow by accreting mass from their surroundings through
large-scale virial shocks. Due to electron acceleration in such a shock, it
should appear as a $\gamma$-ray, hard X-ray, and radio ring, elongated towards
the large-scale filaments feeding the cluster, coincident with a cutoff in the
thermal Sunyaev-Zel'dovich (SZ) signal. However, no such signature was found
until now, and the very existence of cluster virial shocks has remained a
theory. We find preliminary evidence for a large, $\sim 5$ Mpc minor axis
$\gamma$-ray ring around the Coma cluster, elongated towards the large scale
filament connecting Coma and Abell 1367, detected at the nominal $2.7\sigma$
confidence level ($5.1\sigma$ using control signal simulations). The
$\gamma$-ray ring correlates both with a synchrotron signal and with the SZ
cutoff, but not with Galactic tracers. The $\gamma$-ray and radio signatures
agree with analytic and numerical predictions, if the shock deposits $\sim 1\%$
of the thermal energy in relativistic electrons over a Hubble time, and $\sim
1\%$ in magnetic fields. The implied inverse-Compton and synchrotron cumulative
emission from similar shocks can significantly contribute to the diffuse
extragalactic $\gamma$-ray and low frequency radio backgrounds. Our results, if
confirmed, reveal the prolate structure of the hot gas in Coma, the feeding
pattern of the cluster, and properties of the surrounding large scale voids and
filaments. The anticipated detection of such shocks around other clusters would
provide a powerful new cosmological probe.

[14] oai:arXiv.org:1605.03810 [pdf] - 1579932

Dynamical Tides Reexpressed

Comments: 7 pages, 2 figures. The results of Claret & Cunha in Figure 1 were
corrected, and are now consistent with the other stellar models in the
literature. Other results and conclusions are unchanged

Submitted: **2016-05-12**, last modified: **2017-09-13**

Zahn (1975) first put forward and calculated in detail the torque experienced
by stars in a close binary systems due to dynamical tides. His widely used
formula for stars with radiative envelopes and convective cores is expressed in
terms of the stellar radius, even though the torque is actually being applied
to the convective core at the core radius. This results in a large prefactor,
which is very sensitive to the global properties of the star, that multiplies
the torque. This large factor is compensated by a very small multiplicative
factor, $E_{2}$. Although this is mathematically accurate, depending on the
application this can lead to significant errors. The problem is even more
severe, since the calculation of $E_{2}$ itself is non-trivial, and different
authors have obtained inconsistent values of $E_{2}$. Moreover, many codes
(e.g. BSE, StarTrack, MESA) interpolate (and sometimes extrapolate) a fit of
$E_{2}$ values to the stellar mass, often in regimes where this is not sound
practice. We express the torque in an alternate form, cast in terms of
parameters at the envelope-core boundary and a dimensionless coefficient,
$\beta_{2}$. Previous attempts to express the torque in such a form are either
missing an important factor, which depends on the density profile of the star,
or are not easy to implement. We show that $\beta_{2}$ is almost independent of
the properties of the star and its value is approximately unity. Our formula
for the torque is simple to implement and avoids the difficulties associated
with the classic expression.

[15] oai:arXiv.org:1702.00885 [pdf] - 1598007

Research Note: The Expected Spins of Gravitational Wave Sources With
Isolated Field Binary Progenitors

Comments: 5 pages, 4 figures

Submitted: **2017-02-02**

We explore the consequences of dynamical evolution of field binaries composed
of a primary black hole (BH) and a Wolf-Rayet (WR) star in the context of
gravitational wave (GW) source progenitors. We argue, from general
considerations, that the spin of the WR-descendent BH will be maximal in a
significant number of cases due to dynamical effects. In other cases, the spin
should reflect the natal spin of the primary BH which are currently
theoretically unconstrained. We argue that the three currently published LIGO
systems (GW150914, GW151226, LVT151012) suggest that this spin is small. The
resultant effective spin distribution of gravitational wave sources should thus
be bi-model if this classic GW progenitor channel is indeed dominant. While
this is consistent with the LIGO detections thus far, it is in contrast to the
three best-measured high-mass x-ray binary (HMXB) systems. A comparison of the
spin distribution of HMXBs and GW sources should ultimately reveal whether or
not these systems arise from similar astrophysical channels.

[16] oai:arXiv.org:1605.03839 [pdf] - 1457294

GW150914: Spin based constraints on the merger time of the progenitor
system

Comments: 7 pages, 4 figures. Minor changes after referee report

Submitted: **2016-05-12**, last modified: **2016-07-27**

We explore the implications of the observed low spin of GW150914 within the
context of stellar astrophysics and progenitor models. We conclude that many of
the recently proposed scenarios are in marked tension with this observation. We
derive a simple model for the observed spin in the case that the progenitor
system was a field binary composed of a black hole (BH) and a Wolf-Rayet star
and explore the implications of the observed spin for this model. The spin
observation allows us to place a lower limit for the delay time between the
formation of the BH+BH binary and the actual merger, $t_{\textrm{merge}}$. We
use typical values for these systems to derive
$t_{\textrm{merge}}\gtrsim10^{8}\,\textrm{yr}$, which proves to be an important
diagnostic for different progenitor models. We anticipate the next series of
events, and the associated spin parameters, will ultimately yield critical
constraints on formation scenarios and on stellar parameters describing the
late-stage evolution of massive stars.

[17] oai:arXiv.org:1601.03422 [pdf] - 1467558

Neutrino Signal of Collapse-Induced Thermonuclear Supernovae: The Case
for Prompt Black Hole Formation in SN1987A

Comments: 11 pages, 6 figures

Submitted: **2016-01-13**

Collapse-induced thermonuclear explosion (CITE) may explain core-collapse
supernovae (CCSNe). We present a preliminary analysis of the neutrino signal
predicted by CITE and compare it to the neutrino burst of SN1987A. For strong
CCSNe, as SN1987A, CITE predicts a proto-neutron star (PNS) accretion phase,
accompanied by the corresponding neutrino luminosity, that can last a few
seconds and that is cut-off abruptly by black hole (BH) formation. The neutrino
luminosity can later be revived by accretion disc emission after a dead time of
few to a few ten seconds. In contrast, the neutrino mechanism for CCSNe
predicts a shorter PNS accretion phase, followed by a slowly declining PNS
cooling luminosity. We repeat statistical analyses used in the literature to
interpret the neutrino mechanism, and apply them to CITE. The first 1-2 sec of
the neutrino burst are equally compatible with CITE and with the neutrino
mechanism. However, the data hints to a luminosity drop at t=2-3 sec, in some
tension with the neutrino mechanism while being naturally attributed to BH
formation in CITE. The occurrence of neutrino events at 5 sec in SN1987A
suggests that the accretion disc formed by that time. We perform 2D numerical
simulations, showing that CITE may be able to accommodate this disc formation
time while reproducing the ejected $^{56}$Ni mass and ejecta kinetic energy
within factors 2-3 of observations. We estimate the disc neutrino luminosity
and show that it can roughly match the data. This suggests that direct BH
formation is compatible with the neutrino burst of SN1987A. With current
neutrino detectors, the neutrino burst of the next Galactic CCSN may give us
front-row seats to the formation of an event horizon in real time. Access to
phenomena near the event horizon motivates the construction of a few Megaton
neutrino detector that should observe extragalactic CCSNe on a yearly basis.

[18] oai:arXiv.org:1510.07649 [pdf] - 1457120

Can helium envelopes change the outcome of direct white dwarf
collisions?

Comments: 5 Pages, 2 Figure, 3 Tables

Submitted: **2015-10-26**

Collisions of white dwarfs (WDs) have recently been invoked as a possible
mechanism for type Ia supernovae (SNIa). A pivotal feature for the viability of
WD collisions as SNIa progenitors is that a significant fraction of the mass is
highly compressed to the densities required for efficient $^{56}$Ni production
before the ignition of the detonation wave. Previous studies have predominantly
employed model WDs composed entirely of carbon-oxygen (CO), whereas WDs are
expected to have a non-negligible helium envelope. Given that helium is more
susceptible to explosive burning than CO under the conditions characteristic of
WD collision, a legitimate concern is whether or not early time He detonation
ignition can translate to early time CO detonation, thereby drastically
reducing $^{56}$Ni synthesis. We investigate the role of He in determining the
fate of WD collisions by performing a series of two-dimensional hydrodynamics
calculations. We find that a necessary condition for non-trivial reduction of
the CO ignition time is that the He detonation birthed in the contact region
successfully propagates into the unshocked shell. We determine the minimal He
shell mass as a function of the total WD mass that upholds this condition.
Although we utilize a simplified reaction network similar to those used in
previous studies, our findings are in good agreement with detailed
investigations concerning the impact of network size on He shell detonations.
This allows us to extend our results to the case with more realistic burning
physics. Based on the comparison of these findings against evolutionary
calculations of WD compositions, we conclude that most, if not all, WD
collisions will not be drastically impacted by their intrinsic He components.

[19] oai:arXiv.org:1412.1096 [pdf] - 1288308

Failure of a neutrino-driven explosion after core-collapse may lead to a
thermonuclear supernova

Comments: 13 pages, 8 figures. ApJ accepted. Some minor revisions

Submitted: **2014-12-02**, last modified: **2015-09-14**

We demonstrate that $\sim10\,\textrm{s}$ after the core-collapse of a massive
star, a thermonuclear explosion of the outer shells is possible for some
(tuned) initial density and composition profiles, assuming that the neutrinos
failed to explode the star. The explosion may lead to a successful supernova,
as first suggested by Burbidge et al. We perform a series of one-dimensional
(1D) calculations of collapsing massive stars with simplified initial density
profiles (similar to the results of stellar evolution calculations) and various
compositions (not similar to 1D stellar evolution calculations). We assume that
the neutrinos escaped with a negligible effect on the outer layers, which
inevitably collapse. As the shells collapse, they compress and heat up
adiabatically, enhancing the rate of thermonuclear burning. In some cases,
where significant shells of mixed helium and oxygen are present with
pre-collapsed burning times of $\lesssim100\,\textrm{s}$ ($\approx10$ times the
free-fall time), a thermonuclear detonation wave is ignited, which unbinds the
outer layers of the star, leading to a supernova. The energy released is small,
$\lesssim10^{50}\,\textrm{erg}$, and negligible amounts of synthesized material
(including $^{56}$Ni) are ejected, implying that these 1D simulations are
unlikely to represent typical core-collapse supernovae. However, they do serve
as a proof of concept that the core-collapse-induced thermonuclear explosions
are possible, and more realistic two-dimensional and three-dimensional
simulations are within current computational capabilities.

[20] oai:arXiv.org:1401.3347 [pdf] - 1286901

Type Ia Supernovae with Bi-Modal Explosions Are Common -- Possible
Smoking Gun for Direct Collisions of White-Dwarfs

Comments: To be published by MNRAS Letters. Minor changes of the main text, an
Appendix added

Submitted: **2014-01-14**, last modified: **2015-09-01**

We discover clear doubly-peaked line profiles in 3 out of ~20 type Ia
supernovae (SNe Ia) with high-quality nebular-phase spectra. The profiles are
consistently present in three well-separated Co/Fe emission features. The two
peaks are respectively blue-shifted and red-shifted relative to the host
galaxies and are separated by ~5000 km/s. The doubly-peaked profiles directly
reflect a bi-modal velocity distribution of the radioactive Ni56 in the ejecta
that powers the emission of these SNe. Due to their random orientations, only a
fraction of SNe with intrinsically bi-modal velocity distributions will appear
as doubly-peaked spectra. Therefore SNe with intrinsic bi-modality are likely
common, especially among the SNe in the low-luminosity part on the Philips
relation (\Delta m15(B) >~ 1.3; ~40% of all SNe Ia). Such bi-modality is
naturally expected from direct collisions of white dwarfs (WDs) due to the
detonation of both WDs and is demonstrated in a 3D 0.64 M_Sun-0.64 M_Sun WD
collision simulation. In the future, with a large sample of nebular spectra and
a comprehensive set of numerical simulations, the collision model can be
unambiguously tested as the primary channel for type Ia SNe, and the
distribution of nebular line profiles will either be a smoking gun or rule it
out.

[21] oai:arXiv.org:1506.02655 [pdf] - 1129182

The progenitors of core-collapse supernovae suggest thermonuclear origin
for the explosions

Comments: 4 pages, 3 figures, 2 tables

Submitted: **2015-06-08**

Core-collapse supernovae (CCSNe) are the explosions of massive stars
following the collapse of the stars' iron cores. Poznanski (2013) has recently
suggested an observational correlation between the ejecta velocities and the
inferred masses of the red supergiant progenitors of type II-P explosions,
which implies that the kinetic energy of the ejecta ($E_{\textrm{kin}}$)
increases with the mass of the progenitor. I point out that the same conclusion
can be reached from the model-free observed correlation between the ejected
$^{56}$Ni masses ($M_{\textrm{Ni}}$) and the luminosities of the progenitors
for type II supernovae, which was reported by Fraser et al. (2011). This
correlation is in an agreement with the predictions of the collapse-induced
thermonuclear explosions (CITE) for CCSNe and in a possible contradiction with
the predictions of the neutrino mechanism. I show that a correlation between
$M_{\textrm{Ni}}$ and $E_{\textrm{kin}}$ holds for all types of CCSNe
(including type Ibc). This correlation suggests a common mechanism for all
CCSNe, which is predicted for CITE, but is not produced by current simulations
of the neutrino mechanism. Furthermore, the typical values of
$E_{\textrm{kin}}$ and $M_{\textrm{Ni}}$ for type Ibc explosions are larger by
an order of a magnitude than the typical values for II-P explosions, a fact
which disfavors progenitors with the same initial mass range for these
explosions. Instead, the progenitors of type Ibc explosions could be massive
Wolf-Rayet stars, which are predicted to yield strong explosions with low
ejecta masses (as observed) according to CITE. In this case, there is no
deficit of high mass progenitors for CCSNe, which was suggested under the
assumption of a similar mass range for the progenitors of types II-P and Ibc
supernovae.

[22] oai:arXiv.org:1502.03111 [pdf] - 934394

Thermonuclear explosion of rotating massive stars could explain
core-collapse supernovae

Comments: 4 pages, 5 figures, 3 tables

Submitted: **2015-02-10**

It is widely thought that core-collapse supernovae (CCSNe), the explosions of
massive stars following the collapse of the stars' iron cores, is obtained due
to energy deposition by neutrinos. So far, this scenario was not demonstrated
from first principles. Kushnir and Katz (2014) have recently shown, by using
one-dimensional simulations, that if the neutrinos failed to explode the star,
a thermonuclear explosion of the outer shells is possible for some (tuned)
initial profiles. However, the energy released was small and negligible amounts
of ejected $^{56}$Ni were obtained, implying that these one-dimensional
collapse induced thermonuclear explosions (CITE) are unlikely to represent
typical CCSNe. Here I provide evidence supporting a scenario in which the
majority of CCSNe are the result of CITE. I use two-dimensional simulations to
show that collapse of stars that include slowly (few percent of breakup)
rotating $\sim0.1-10\,M_{\odot}$ shells of mixed helium-oxygen, leads to an
ignition of a thermonuclear detonation wave that unbinds the stars' outer
layers. Simulations of massive stars with different properties show that CITE
is a robust process, and results in explosions with kinetic energies in the
range of $10^{49}-10^{52}\,\textrm{erg}$, and $^{56}$Ni yields of up to
$\sim\,M_{\odot}$, which are correlated, in agreement with observations for the
majority of CCSNe. Stronger explosions are predicted from higher mass
progenitors that leave more massive remnants, in contrast to the neutrino
mechanism. Neutron stars are produced in weak ($\lt10^{51}\,\textrm{erg}$)
explosions, while strong ($\gt10^{51}\,\textrm{erg}$) explosions leave black
hole remnants.

[23] oai:arXiv.org:1402.7083 [pdf] - 790216

Luminosity function suggests up to 100 white dwarfs within 20 pc may be
hiding in multiple systems

Comments: 3 pages, 2 figures

Submitted: **2014-02-27**

We examine the luminosity function of white dwarfs (WDs) in the local
``complete'' WD sample ($d<20$ pc) of Holberg et. al. 2008. We find that the
fraction of bright and young WDs is anomalously high among the WDs detected in
multiple systems with main sequence (MS) companions compared to that of the
single WDs and theoretical expectations. This indicates a significant
observation bias against finding relatively faint WDs in multiple systems. At
the bright end ($M_V<11.5$), the amount of WDs with MS companions is
approximately equal to that of single white dwarfs, indicating that $\gtrsim
50\%$ of WDs have MS companions, consistent with the high multiplicity fraction
of early-type MS stars. If true, a significant fraction of WDs in multiple
systems within 20 pc may have not been detected yet, and the number density of
WDs in the solar neighborhood and elsewhere may be up to twice as much as
presently believed.

[24] oai:arXiv.org:1303.1180 [pdf] - 1165023

Head-on collisions of white dwarfs in triple systems could explain type
Ia supernova

Comments: 10 pages, 5 figures, published in ApJL. Revised manuscript has
calculations with higher resolution and an expanded discussion of previous
work

Submitted: **2013-03-05**, last modified: **2013-11-19**

Type Ia supernovae (SNe Ia), thermonuclear explosions of carbon-oxygen white
dwarfs (CO-WDs), are currently the best cosmological "standard candles", but
the triggering mechanism of the explosion is unknown. It was recently shown
that the rate of head-on collisions of typical field CO-WDs in triple systems
may be comparable to the SNe Ia rate. Here we provide evidence supporting a
scenario in which the majority of SNe Ia are the result of such head-on
collisions of CO-WDs. In this case, the nuclear detonation is due to a well
understood shock ignition, devoid of commonly introduced free parameters such
as the deflagration velocity or transition to detonation criteria. By using
two-dimensional hydrodynamical simulations with a fully resolved ignition
process, we show that zero-impact-parameter collisions of typical CO-WDs with
masses $0.5-1\,M_{\odot}$ result in explosions that synthesize $^{56}$Ni masses
in the range of $\sim0.1-1\,M_{\odot}$, spanning the wide distribution of
yields observed for the majority of SNe Ia. All collision models yield the same
late-time ($>60$ days since explosion) bolometric light curve when normalized
by $^{56}$Ni masses (to better than $30\%$), in agreement with observations.
The calculated widths of the $^{56}$Ni-mass-weighted-line-of-sight velocity
distributions are correlated with the calculated $^{56}$Ni yield, agreeing with
the observed correlation. The strong correlation, shown here for the first
time, between $^{56}$Ni yield and total mass of the colliding CO-WDs
(insensitive to their mass ratio), is suggestive as the source for the
continuous distribution of observed SN Ia features, possibly including the
Philips relation.

[25] oai:arXiv.org:1311.1209 [pdf] - 1180507

Early Hydrodynamic Evolution of a Stellar Collision

Comments: 7 pages, 5 figures (without appendices). Submitted to ApJ

Submitted: **2013-11-05**

The early phase of the hydrodynamic evolution following collision of two
stars is analyzed. Two strong shocks propagate at a constant velocity (which is
a small fraction of the velocity of the approaching stars) from the contact
surface toward the center of each star. The shocked region near the contact
surface has a planar symmetry and a uniform pressure. The density vanishes at
the (Lagrangian) surface of contact and the speed of sound diverges there. The
temperature, however, reaches a finite value, since as the density vanishes,
the finite pressure is radiation dominated. For Carbon-Oxygen white dwarfs
collisions this temperature is too low for any appreciable nuclear burning at
early times. The divergence of the speed of sound limits numerical studies of
stellar collisions, as it makes convergence tests exceedingly expensive unless
dedicated schemes are used. We provide a new one-dimensional Lagrangian
numerical scheme to achieve this. Self-similar planar solutions are derived for
zero-impact parameter collisions between two identical stars, under some
simplifying assumptions. These solutions provide rough approximations that
capture the main features of the flow and allow a general study as well as a
detailed numerical verification test problem. The self-similar solution in the
upstream frame is the planar version of previous piston problems that were
studied in cylindrical and spherical symmetries. We found it timely to present
a global picture of self similar piston problems. In particular, we derive new
results regarding the non trivial transition to accelerating shocks at
sufficiently declining densities (not relevant for collisions).

[26] oai:arXiv.org:1304.4234 [pdf] - 721118

What can we really learn about Magnetic Fields in Galaxy Clusters from
Faraday Rotation observations?

Comments: This version includes an error in the analysis of sources, for which
it was concluded that the RM observations may be accounted for by a
homogeneous field (pointed to us by Simon White and Larry Rudnick). A revised
version will soon be posted. The properties of fields inferred by our
analysis for sources for which inhomogeneous fields are required are not
expected to be modified by the correction

Submitted: **2013-04-15**, last modified: **2013-09-18**

We construct a simple and robust approach for deriving constraints on
magnetic fields in galaxy clusters from rotation measure (RM) maps. Relaxing
the commonly used assumptions of a correlation between the magnetic field
strength and the plasma density and of a power-law (in wave number) magnetic
field power spectrum, and using an efficient numerical analysis method, we test
the consistency of a wide range of magnetic field models with RM maps of 11
extended sources in 5 clusters, for which the data were made available to us.
We show that the data reveal no indication for a radial dependence of the
average magnetic field strength, and in particular no indication for a
correlation between the gas density and the field strength. The RM maps of a
considerable fraction of the sources either require or are consistent with the
presence of a spatially uniform magnetic field of a relatively small strength,
0.02-0.3 muG, which contributes significantly to the RM. The RM maps of all but
one source do not require a power-law magnetic field power spectrum, and most
are consistent with a power spectrum dominated by a single wave length. The
uncertainties in the magnetic field strengths (and spatial correlation lengths)
derived from RM maps exceed an order of magnitude (and often more). These
uncertainties imply, in particular, that there is no indication in current RM
data for a systematic difference between the magnetic field strengths in
radio-halo clusters and in radio-quiet clusters. With the improvement expected
in the near future of the quality and quantity of RM data, our analysis method
will enable one to derive more accurate constraints on magnetic fields in
galaxy clusters.

[27] oai:arXiv.org:1301.6766 [pdf] - 620014

An exact integral relation between the Ni56 mass and the bolometric
light curve of a type Ia supernova

Comments: 1 paragraph

Submitted: **2013-01-28**

An exact relation between the Ni56 mass and the bolometric light curve of a
type Ia supernova can be derived as follows, using the following excellent
approximations: 1. the emission is powered solely by Ni56-> Co56 ->Fe56; 2.
each mass element propagates at a non-relativistic velocity which is constant
in time (free coasting); and 3. the internal energy is dominated by radiation.
Under these approximations, the energy E(t) carried by radiation in the ejecta
satisfies: dE/dt=-E(t)/t-L(t)+Q(t), where Q(t) is the deposition of energy by
the decay which is precisely known and L(t) is the bolometric luminosity. By
multiplying this equation by time and integrating over time we find:
E(t)*t=\int_0^t Q(t')t'dt' -\int_0^t L(t')t'dt'. At late time, t>> t_peak, the
energy inside the ejecta decreases rapidly due to its escape, and thus we have
\int_0^t Q(t')t'dt'=\int_0^t L(t')t'dt'. This relation is correct regardless of
the opacities, density distribution or Ni56 deposition distribution in the
ejecta and is very different from "Arnett's rule", L_peak ~ Q(t_peak). By
comparing \int_0^t Q(t')t'dt' with \int_0^t L(t')t'dt' at t~40 day after the
explosion, the mass of Ni56 can be found directly from UV, optical and infrared
observations with modest corrections due to the unobserved gamma-rays and due
to the small residual energy in the ejecta, E(t)*t>0.

[28] oai:arXiv.org:1108.4690 [pdf] - 1083585

Imploding ignition waves: I. one dimensional analysis

Comments: 12 pages, 10 figures. Somewhat modified, published in ApJ

Submitted: **2011-08-23**, last modified: **2012-06-02**

We show that converging spherical and cylindrical shock waves may ignite a
detonation wave in a combustible medium, provided the radius at which the
shocks become strong exceeds a critical radius, R_c. An approximate analytic
expression for R_c is derived for an ideal gas equation of state and a simple
(power-law-Arrhenius) reaction law, and shown to reproduce the results of
numerical solutions. For typical acetylene-air experiments we find R_c~0.1 mm
(spherical) and R_c~1 mm (cylindrical). We suggest that the deflagration to
detonation transition (DDT) observed in these systems may be due to converging
shocks produced by the turbulent deflagration flow, which reaches sub (but
near) sonic velocities on scales >>R_c. Our suggested mechanism differs from
that proposed by Zel'dovich et al., in which a fine-tuned spatial gradient in
the chemical induction time is required to be maintained within the turbulent
deflagration flow. Our analysis may be readily extended to more complicated
equations of state and reaction laws. An order of magnitude estimate of R_c
within a white dwarf at the pre-detonation conditions believed to lead to Type
Ia supernova explosions is 0.1 km, suggesting that our proposed mechanism may
be relevant for DDT initiation in these systems. The relevance of our proposed
ignition mechanism to DDT initiation may be tested by both experiments and
numerical simulations.

[29] oai:arXiv.org:1002.3869 [pdf] - 1025301

Closing the gap in the solutions of the strong explosion problem: An
expansion of the family of second-type self-similar solutions

Comments: 10 pages, 11 figures, somewhat revised, version accepted to ApJ

Submitted: **2010-02-20**, last modified: **2010-08-27**

Shock waves driven by the release of energy at the center of a cold ideal gas
sphere of initial density rho\propto r^{-omega} approach a self-similar (SLS)
behavior, with velocity \dot{R}\propto R^delta, as R->\infty. For omega>3 the
solutions are of the second-type, i.e., delta is determined by the requirement
that the flow should include a sonic point. No solution satisfying this
requirement exists, however, in the 3\leq omega\leq omega_{g}(gamma) ``gap''
(\omega_{g}=3.26 for adiabatic index gamma=5/3). We argue that second-type
solutions should not be required in general to include a sonic point. Rather,
it is sufficient to require the existence of a characteristic line r_c(t), such
that the energy in the region r_c(t)<r<R approaches a constant as R->\infty,
and an asymptotic solution given by the SLS solution at r_c(t)<r<R and
deviating from it at r<r_c may be constructed. The two requirements coincide
for omega>omega_g and the latter identifies delta=0 solutions as the asymptotic
solutions for 3\leq omega\leq omega_{g} (as suggested by Gruzinov03). In these
solutions, r_c is a C_0 characteristic. It is difficult to check, using
numerical solutions of the hydrodynamic equations, whether the flow indeed
approaches a delta=0 SLS behavior as R->\infty, due to the slow convergence to
SLS for omega~3. We show that in this case the flow may be described by a
modified SLS solution, d\ln\dot{R}/d\ln R=delta with slowly varying delta(R),
eta\equiv d delta/d\ln R<<1, and spatial profiles given by a sum of the SLS
solution corresponding to the instantaneous value of delta and a SLS correction
linear in eta. The modified SLS solutions provide an excellent approximation to
numerical solutions obtained for omega~3 at large R, with delta->0 (and
eta\neq0) for 3\leq omega\leq omega_{g}. (abridged)

[30] oai:arXiv.org:0903.2275 [pdf] - 122903

Magnetic fields and cosmic rays in clusters of galaxies

Comments: 7 pages, 2 figures, somewhat revised, published in JCAP

Submitted: **2009-03-12**, last modified: **2010-02-22**

We argue that the observed correlation between the radio luminosity and the
X-ray luminosity in radio emitting galaxy clusters implies that the radio
emission is due to secondary electrons that are produced by p-p interactions
and lose their energy by emitting synchrotron radiation in a strong magnetic
field, B>(8\pi a T_{CMB}^4)^{1/2}\simeq 3\muG. We construct a simple model that
naturally explains the correlation, and show that the observations provide
stringent constraints on cluster magnetic fields and cosmic rays (CRs): Within
the cores of clusters, the ratio beta_{core} between the CR energy (per
logarithmic particle energy interval) and the thermal energy is beta_{core}\sim
2*10^{-4}; The source of these CRs is most likely the cluster accretion shock,
which is inferred to deposit in CRs ~ 0.1 of the thermal energy it generates;
The diffusion time of 100 GeV CRs over scales \gtrsim100 kpc is not short
compared to the Hubble time; Cluster magnetic fields are enhanced by mergers to
\gtrsim 1 % of equipartition, and decay (to <1 muG) on 1 Gyr time scales. The
inferred value of beta_{core} implies that high energy gamma-ray emission from
secondaries at cluster cores will be difficult to detect with existing and
planned instruments.

[31] oai:arXiv.org:0903.2271 [pdf] - 122902

Nonthermal emission from clusters of galaxies

Comments: 22 pages, 16 figures, somewhat revised, published in JCAP

Submitted: **2009-03-12**, last modified: **2010-02-22**

We show that the spectral and radial distribution of the nonthermal emission
of massive, M>10^{14.5}M_sun, galaxy clusters (GCs) may be approximately
described by simple analytic expressions, which depend on the GC thermal X-ray
properties and on two model parameter, beta_{core} and eta_e. beta_{core} is
the ratio of CR energy density (within a logarithmic CR energy interval) and
the thermal energy density at the GC core, and eta_{e(p)} is the fraction of
the thermal energy generated in strong collisionless shocks, which is deposited
in CR electrons (protons). Using a simple analytic model for the evolution of
ICM CRs, which are produced by accretion shocks (primary CRs), we find that
beta_{core} ~ eta_{p}/200, nearly independent of GC mass and with a scatter
Delta ln(beta_{core}) ~ 1 between GCs of given mass. We show that the HXR and
gamma-ray luminosities produced by IC scattering of CMB photons by primary
electrons exceed the luminosities produced by secondary particles (generated in
hadronic interactions within the GC) by factors ~500(eta_e/eta_p)(T/10
keV)^{-1/2} and ~150(eta_e/eta_p)(T/10 keV)^{-1/2} respectively, where T is the
GC temperature. Secondary particle emission may dominate at the radio and VHE
(> 1 TeV) gamma-ray bands. Our model predicts, in contrast with some earlier
work, that the HXR and gamma-ray emission from GCs are extended, since the
emission is dominated at these energies by primary electrons. Our predictions
are consistent with the observed nonthermal emission of the Coma cluster for
eta_peta_e ~ 0.1. The implications of our predictions to future HXR
observations (e.g. by NuStar, Simbol-X) and to (space/ground based) gamma-ray
observations (e.g. by Fermi, HESS, MAGIC, VERITAS) are discussed. Finally, we
show that our model's results agree with results of detailed numerical
calculations.

[32] oai:arXiv.org:0905.1950 [pdf] - 122909

Hard X-ray emission from accretion shocks around galaxy clusters

Comments: 7 pages, 3 figures, somewhat revised, published in JCAP

Submitted: **2009-05-12**, last modified: **2010-02-22**

We show that the hard X-ray (HXR) emission observed from several galaxy
clusters is naturally explained by a simple model, in which the nonthermal
emission is produced by inverse Compton scattering of cosmic microwave
background photons by electrons accelerated in cluster accretion shocks: The
dependence of HXR surface brightness on cluster temperature is consistent with
that predicted by the model, and the observed HXR luminosity is consistent with
the fraction of shock thermal energy deposited in relativistic electrons being
\lesssim 0.1. Alternative models, where the HXR emission is predicted to be
correlated with the cluster thermal emission, are disfavored by the data. The
implications of our predictions to future HXR observations (e.g. by NuStar,
Simbol-X) and to (space/ground based) gamma-ray observations (e.g. by Fermi,
HESS, MAGIC, VERITAS) are discussed.

[33] oai:arXiv.org:astro-ph/0502125 [pdf] - 70916

The stability of decelerating shocks revisited

Comments: 10 pages, 9 figures. Accepted to ApJ; Expanded discussion of boundary
conditions

Submitted: **2005-02-06**, last modified: **2005-08-21**

We present a new method for analyzing the global stability of the Sedov-von
Neumann-Taylor self-similar solutions, describing the asymptotic behavior of
spherical decelerating shock waves, expanding into ideal gas with density
\propto r^{-\omega}. Our method allows to overcome the difficulties associated
with the non-physical divergences of the solutions at the origin. We show that
while the growth rates of global modes derived by previous analyses are
accurate in the large wave number (small wavelength) limit, they do not
correctly describe the small wave number behavior for small values of the
adiabatic index \gamma. Our method furthermore allows to analyze the stability
properties of the flow at early times, when the flow deviates significantly
from the asymptotic self-similar behavior. We find that at this stage the
perturbation growth rates are larger than those obtained for unstable
asymptotic solutions at similar [\gamma,\omega]. Our results reduce the
discrepancy that exists between theoretical predictions and experimental
results.