sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Kushnir, Doron

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.