Showing 1–50 of 60 results for author: Kushnir, D

A Cepheid systematics-free test of $H_0$ to $\lesssim2.5\%$ accuracy using SH0ES photometry

Authors: Doron Kushnir, Amir Sharon

Abstract: The recent SH0ES determination of the Hubble constant, $H_0=73.04\pm1.04$ km/s/Mpc, deviates significantly by $\approx5σ$ from the \textit{Planck} value, stimulating discussions on cosmological model extensions. To minimize statistical uncertainty and mitigate sensitivity to systematic errors in any single anchor distance determination, SH0ES combines Cepheids from various observations, including… ▽ More The recent SH0ES determination of the Hubble constant, $H_0=73.04\pm1.04$ km/s/Mpc, deviates significantly by $\approx5σ$ from the \textit{Planck} value, stimulating discussions on cosmological model extensions. To minimize statistical uncertainty and mitigate sensitivity to systematic errors in any single anchor distance determination, SH0ES combines Cepheids from various observations, including those from Type Ia supernova (SNe Ia) host galaxies, NGC 4258, and closer galaxies (MW, LMC, SMC, and M31), although this mixed sample may introduce unknown or subtle systematic errors due to comparing distant and closer Cepheids. To address this, we propose a subset excluding Cepheids from the closer galaxies, retaining only the NGC 4258 water megamasers as a single anchor, circumventing potential systematic errors associated with observational methods and reduction techniques. Focusing solely on these Cepheids yields competitive statistical errors, approximately $2.5\%$, sufficient to identify a $\approx3σ$ tension with the \textit{Planck} $H_0$ value. Our approach offers an opportunity to utilize optical photometry with systematic uncertainty smaller than the statistical uncertainty, potentially achieving higher precision than NIR photometry, given the lower optical background. However, currently the optical photometry sample's fidelity does not match that of NIR photometry. The significant Hubble tension obtained is unrelated to Cepheids and we discuss other options. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Comments: 13 pages, 8 figures

Coma cluster $γ$-ray and radio emission is consistent with a secondary electron origin for the radio halo

Authors: Doron Kushnir, Uri Keshet, Eli Waxman

Abstract: Observations of diffuse, non-thermal radio emission spanning several megaparsecs have been documented in over 100 galaxy clusters. This emission, classified as giant radio halos (GHs), mini halos, and radio relics based mainly on their location and morphology, is interpreted as synchrotron radiation and implies the presence of relativistic electrons and magnetic fields in the intra-cluster medium… ▽ More Observations of diffuse, non-thermal radio emission spanning several megaparsecs have been documented in over 100 galaxy clusters. This emission, classified as giant radio halos (GHs), mini halos, and radio relics based mainly on their location and morphology, is interpreted as synchrotron radiation and implies the presence of relativistic electrons and magnetic fields in the intra-cluster medium (ICM). GHs were initially thought to be generated by secondary electrons resulting from inelastic $p+p\rightarrow X+π^{\pm}$ collisions. However, recent literature has leaned towards primary-electron turbulent (re)acceleration models, partly due to claimed upper limits on the $γ$-ray emission from $π^0$ decay. We demonstrate that the observed GH and $γ$-ray flux in the Coma cluster are consistent with a secondary origin for the GH across a broad range of magnetic field values. Although the constraints on magnetic field configuration are not stringent, they align well with previous estimates for Coma. Within this magnetic field range, the energy density of cosmic-ray protons (CRp) constitutes a few percent to tens of percent of the ICM energy density, as predicted and observed for a sample of radio-emitting galaxy clusters. Notably, we detect a rise in the ratio of CRp to ICM energy densities towards the outer regions of the cluster. This phenomenon was anticipated to arise from either adiabatic compression of CRp accelerated by accretion shocks or, more likely, from strong CRp diffusion. △ Less

Submitted 19 April, 2024; originally announced April 2024.

Comments: 7 pages, 3 figures

Asteroid collisions: expected visibility and rate

Authors: E. O. Ofek, D. Polishook, D. Kushnir, G. Nir, S. Ben-Ami, Y. Shvartzvald, N. L. Strotjohann, E. Segre, A. Blumenzweig, M. Engel, D. Bodewits, J. W. Noonan

Abstract: Asteroid collisions are one of the main processes responsible for the evolution of bodies in the main belt. Using observations of the Dimorphos impact by the DART spacecraft, we estimate how asteroid collisions in the main belt may look in the first hours after the impact. If the DART event is representative of asteroid collisions with a ~1m size impactor, then the light curves of these collisions… ▽ More Asteroid collisions are one of the main processes responsible for the evolution of bodies in the main belt. Using observations of the Dimorphos impact by the DART spacecraft, we estimate how asteroid collisions in the main belt may look in the first hours after the impact. If the DART event is representative of asteroid collisions with a ~1m size impactor, then the light curves of these collisions will rise on time scales of about >100s and will remain bright for about one hour. Next, the light curve will decay on a few hours time scale to an intermediate luminosity level in which it will remain for several weeks, before slowly returning to its baseline magnitude. This estimate suffers from several uncertainties due to, e.g., the diversity of asteroid composition, their material strength, and spread in collision velocities. We estimate that the rate of collisions in the main belt with energy similar or larger than the DART impact is of the order of 7000 per year (+/-1dex). The large range is due to the uncertainty in the abundance of ~1-m size asteroids. We estimate the magnitude distribution of such events in the main belt, and we show that ~6% of these events may peak at magnitudes brighter than 21. The detection of these events requires a survey with <1hr cadence and may contribute to our understanding of the asteroids' size distribution, collisional physics, and dust production. With an adequate survey strategy, new survey telescopes may regularly detect asteroid collisions. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 7 pages, AJ in press

UV to near-IR observations of the DART-Dimorphos collision

Authors: E. O. Ofek, D. Kushnir, D. Polishook, E. Waxman, A. Tohuvavohu, S. Ben-Ami, B. Katz, O. Gnat, N. L. Strotjohann, E. Segre, A. Blumenzweig, Y. Sofer-Rimalt, O. Yaron, A. Gal-Yam, Y. Shvartzvald, M. Engel, S. B. Cenko, O. Hershko

Abstract: The impact of the Double Asteroid Redirection Test (DART) spacecraft with Dimorphos allows us to study asteroid collision physics, including momentum transfer, the ejecta properties, and the visibility of such events in the Solar System. We report observations of the DART impact in the ultraviolet (UV), visible light, and near-infrared (IR) wavelengths. The observations support the existence of at… ▽ More The impact of the Double Asteroid Redirection Test (DART) spacecraft with Dimorphos allows us to study asteroid collision physics, including momentum transfer, the ejecta properties, and the visibility of such events in the Solar System. We report observations of the DART impact in the ultraviolet (UV), visible light, and near-infrared (IR) wavelengths. The observations support the existence of at least two separate components of the ejecta: a fast and a slow component. The fast-ejecta component is composed of a gaseous phase, moving at about 1.6 km/s with a mass of <10^4 kg. The fast ejecta is detected in the UV and visible light, but not in the near-IR $z$-band observations. Fitting a simplified optical thickness model to these observations allows us to constrain some of the properties of the fast ejecta, including its scattering efficiency and the opacity of the gas. The slow ejecta component is moving at typical velocities of up to about 10 m/s. It is composed of micrometer-size particles, that have a scattering efficiency, at the direction of the observer, of the order of 10^-3 and a total mass of about 10^6 kg. The larger particles in the slow ejecta, whose size is bound to be in the range between ~1 mm to ~1 m, likely have a scattering efficiency larger than that of the pre-impact Didymos system. △ Less

Submitted 20 November, 2023; originally announced November 2023.

Comments: Submitted to MNRAS

Photometric prioritization of neutron star merger candidates

Authors: E. O. Ofek, N L. Strotjohann, I. Arcavi, A. Gal-Yam, D. Kushnir, E. Waxman, M. M. Kasliwal, A. Drake, M. Graham, J. Purdum, B. Rusholme, Y. Sharma, R. Smith, A. Wold, B. F. Healy

Abstract: Rapid identification of the optical counterparts of Neutron Star (NS) merger events discovered by gravitational wave detectors may require observing a large error region and sifting through a large number of transients to identify the object of interest. Given the expense of spectroscopic observations, a question arises: How can we utilize photometric observations for candidate prioritization, and… ▽ More Rapid identification of the optical counterparts of Neutron Star (NS) merger events discovered by gravitational wave detectors may require observing a large error region and sifting through a large number of transients to identify the object of interest. Given the expense of spectroscopic observations, a question arises: How can we utilize photometric observations for candidate prioritization, and what kinds of photometric observations are needed to achieve this goal? NS merger kilonova exhibits low ejecta mass (~5x10^-2 solar mass) and a rapidly evolving photospheric radius (with a velocity ~0.2c). As a consequence, these sources display rapid optical-flux evolution. Indeed, selection based on fast flux variations is commonly used for young supernovae and NS mergers. In this study, we leverage the best currently available flux-limited transient survey - the Zwicky Transient Facility Bright Transient Survey - to extend and quantify this approach. We focus on selecting transients detected in a 3-day cadence survey and observed at a one-day cadence. We explore their distribution in the phase space defined by g-r, g-dot, and r-dot. Our analysis demonstrates that for a significant portion of the time during the first week, the kilonova AT 2017gfo stands out in this phase space. It is important to note that this investigation is subject to various biases and challenges; nevertheless, it suggests that certain photometric observations can be leveraged to identify transients with the highest probability of being fast-evolving events. We also find that a large fraction (~0.75) of the transient candidates with |g-dot|>0.7 mag/day, are cataclysmic variables or active galactic nuclei with radio counterparts. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: 9 pages, 4 figures, MNRAS in press

Resolving the explosion of supernova 2023ixf in Messier 101 within its complex circumstellar environment

Authors: E. A. Zimmerman, I. Irani, P. Chen, A. Gal-Yam, S. Schulze, D. A. Perley, J. Sollerman, A. V. Filippenko, T. Shenar, O. Yaron, S. Shahaf, R. J. Bruch, E. O. Ofek, A. De Cia, T. G. Brink, Y. Yang, S. S. Vasylyev, S. Ben Ami, M. Aubert, A. Badash, J. S. Bloom, P. J. Brown, K. De, G. Dimitriadis, C. Fransson , et al. (32 additional authors not shown)

Abstract: Observing a supernova explosion shortly after it occurs can reveal important information about the physics of stellar explosions and the nature of the progenitor stars of supernovae (SNe). When a star with a well-defined edge explodes in vacuum, the first photons to escape from its surface appear as a brief shock-breakout flare. The duration of this flare can extend to at most a few hours even for… ▽ More Observing a supernova explosion shortly after it occurs can reveal important information about the physics of stellar explosions and the nature of the progenitor stars of supernovae (SNe). When a star with a well-defined edge explodes in vacuum, the first photons to escape from its surface appear as a brief shock-breakout flare. The duration of this flare can extend to at most a few hours even for nonspherical breakouts from supergiant stars, after which the explosion ejecta should expand and cool. Alternatively, for stars exploding within a distribution of sufficiently dense optically thick circumstellar material, the first photons escape from the material beyond the stellar edge, and the duration of the initial flare can extend to several days, during which the escaping emission indicates photospheric heating. The difficulty in detecting SN explosions promptly after the event has so far limited data regarding supergiant stellar explosions mostly to serendipitous observations that, owing to the lack of ultraviolet (UV) data, were unable to determine whether the early emission is heating or cooling, and hence the nature of the early explosion event. Here, we report observations of SN 2023ixf in the nearby galaxy M101, covering the early days of the event. Using UV spectroscopy from the Hubble Space Telescope (HST) as well as a comprehensive set of additional multiwavelength observations, we trace the photometric and spectroscopic evolution of the event and are able to temporally resolve the emergence and evolution of the SN emission. △ Less

Submitted 27 March, 2024; v1 submitted 16 October, 2023; originally announced October 2023.

Journal ref: Nature 627, 759 (2024)

A 12.4 day periodicity in a close binary system after a supernova

Authors: Ping Chen, Avishay Gal-Yam, Jesper Sollerman, Steve Schulze, Richard S. Post, Chang Liu, Eran O. Ofek, Kaustav K. Das, Christoffer Fremling, Assaf Horesh, Boaz Katz, Doron Kushnir, Mansi M. Kasliwal, Shri R. Kulkarni, Dezi Liu, Xiangkun Liu, Adam A. Miller, Kovi Rose, Eli Waxman, Sheng Yang, Yuhan Yao, Barak Zackay, Eric C. Bellm, Richard Dekany, Andrew J. Drake , et al. (15 additional authors not shown)

Abstract: Neutron stars and stellar-mass black holes are the remnants of massive star explosions. Most massive stars reside in close binary systems, and the interplay between the companion star and the newly formed compact object has been theoretically explored, but signatures for binarity or evidence for the formation of a compact object during a supernova explosion are still lacking. Here we report a stri… ▽ More Neutron stars and stellar-mass black holes are the remnants of massive star explosions. Most massive stars reside in close binary systems, and the interplay between the companion star and the newly formed compact object has been theoretically explored, but signatures for binarity or evidence for the formation of a compact object during a supernova explosion are still lacking. Here we report a stripped-envelope supernova, SN 2022jli, which shows 12.4-day periodic undulations during the declining light curve. Narrow H$α$ emission is detected in late-time spectra with concordant periodic velocity shifts, likely arising from hydrogen gas stripped from a companion and accreted onto the compact remnant. A new Fermi/LAT $γ$-ray source is temporally and positionally consistent with SN 2022jli. The observed properties of SN 2022jli, including periodic undulations in the optical light curve, coherent H$α$ emission shifting, and evidence for association with a $γ$-ray source, point to the explosion of a massive star in a binary system leaving behind a bound compact remnant. Mass accretion from the companion star onto the compact object powers the light curve of the supernova and generates the $γ$-ray emission. △ Less

Submitted 11 October, 2023; originally announced October 2023.

Comments: Accepted for publication in Nature

Reassessing the Constraints from SH0ES Extragalactic Cepheid Amplitudes on Systematic Blending Bias

Authors: Amir Sharon, Doron Kushnir, Wenlong Yuan, Lucas Macri, Adam Riess

Abstract: The SH0ES collaboration Hubble constant determination is in a ${\sim}5σ$ difference with the $Planck$ value, known as the Hubble tension. The accuracy of the Hubble constant measured with extragalactic Cepheids depends on robust stellar-crowding background estimation. Riess et al. (R20) compared the light curves amplitudes of extragalactic and MW Cepheids to constrain an unaccounted systematic ble… ▽ More The SH0ES collaboration Hubble constant determination is in a ${\sim}5σ$ difference with the $Planck$ value, known as the Hubble tension. The accuracy of the Hubble constant measured with extragalactic Cepheids depends on robust stellar-crowding background estimation. Riess et al. (R20) compared the light curves amplitudes of extragalactic and MW Cepheids to constrain an unaccounted systematic blending bias, $γ=-0.029\pm0.037\,\rm{mag}$, which cannot explain the required, $γ=0.24\pm0.05\,\rm{mag}$, to resolve the Hubble tension. Further checks by Riess et al. demonstrate that a possible blending is not likely related to the size of the crowding correction. We repeat the R20 analysis, with the following main differences: (1) we limit the extragalactic and MW Cepheids comparison to periods $P\lesssim50\,\rm{d}$, since the number of MW Cepheids with longer periods is minimal; (2) we use publicly available data to recalibrate amplitude ratios of MW Cepheids in standard passbands; (3) we remeasure the amplitudes of Cepheids in NGC 5584 and NGC 4258 in two HST filters ($F555W$ and $F350LP$) to improve the empirical constraint on their amplitude ratio $A^{555}/A^{350}$. We show that the filter transformations introduce an ${\approx}0.04\,\rm{mag}$ uncertainty in determining $γ$, not included by R20. While our final estimate, $γ=0.013\pm0.057\,\rm{mag}$, is consistent with the value derived by R20 and is consistent with no bias, the error is somewhat larger, and the best fitting value is shifted by ${\approx}0.04\,\rm{mag}$ and closer to zero. Future observations, especially with JWST, would allow better calibration of $γ$. △ Less

Submitted 19 February, 2024; v1 submitted 23 May, 2023; originally announced May 2023.

Comments: 19 pages, 10 figures. Accepted for publication in MNRAS

ULTRASAT: A wide-field time-domain UV space telescope

Authors: Y. Shvartzvald, E. Waxman, A. Gal-Yam, E. O. Ofek, S. Ben-Ami, D. Berge, M. Kowalski, R. Bühler, S. Worm, J. E. Rhoads, I. Arcavi, D. Maoz, D. Polishook, N. Stone, B. Trakhtenbrot, M. Ackermann, O. Aharonson, O. Birnholtz, D. Chelouche, D. Guetta, N. Hallakoun, A. Horesh, D. Kushnir, T. Mazeh, J. Nordin , et al. (19 additional authors not shown)

Abstract: The Ultraviolet Transient Astronomy Satellite (ULTRASAT) is scheduled to be launched to geostationary orbit in 2026. It will carry a telescope with an unprecedentedly large field of view (204 deg$^2$) and NUV (230-290nm) sensitivity (22.5 mag, 5$σ$, at 900s). ULTRASAT will conduct the first wide-field survey of transient and variable NUV sources and will revolutionize our ability to study the hot… ▽ More The Ultraviolet Transient Astronomy Satellite (ULTRASAT) is scheduled to be launched to geostationary orbit in 2026. It will carry a telescope with an unprecedentedly large field of view (204 deg$^2$) and NUV (230-290nm) sensitivity (22.5 mag, 5$σ$, at 900s). ULTRASAT will conduct the first wide-field survey of transient and variable NUV sources and will revolutionize our ability to study the hot transient universe: It will explore a new parameter space in energy and time-scale (months long light-curves with minutes cadence), with an extra-Galactic volume accessible for the discovery of transient sources that is $>$300 times larger than that of GALEX and comparable to that of LSST. ULTRASAT data will be transmitted to the ground in real-time, and transient alerts will be distributed to the community in $<$15 min, enabling a vigorous ground-based follow-up of ULTRASAT sources. ULTRASAT will also provide an all-sky NUV image to $>$23.5 AB mag, over 10 times deeper than the GALEX map. Two key science goals of ULTRASAT are the study of mergers of binaries involving neutron stars, and supernovae: With a large fraction ($>$50%) of the sky instantaneously accessible, fast (minutes) slewing capability and a field-of-view that covers the error ellipses expected from GW detectors beyond 2025, ULTRASAT will rapidly detect the electromagnetic emission following BNS/NS-BH mergers identified by GW detectors, and will provide continuous NUV light-curves of the events; ULTRASAT will provide early (hour) detection and continuous high (minutes) cadence NUV light curves for hundreds of core-collapse supernovae, including for rarer supernova progenitor types. △ Less

Submitted 27 April, 2023; originally announced April 2023.

Comments: 40 pages, 16 figures, 3 tables. Submitted to the AAS journals

The Large Array Survey Telescope -- System Overview and Performances

Authors: E. O. Ofek, S. Ben-Ami, D. Polishook, E. Segre, A. Blumenzweig, N. L. Strotjohann, O. Yaron, Y. M. Shani, S. Nachshon, Y. Shvartzvald, O. Hershko, M. Engel, M. Segre, N. Segev, E. Zimmerman, G. Nir, Y. Judkovsky, A. Gal-Yam, B. Zackay, E. Waxman, D. Kushnir, P. Chen, R. Azaria, I. Manulis, O. Diner , et al. (16 additional authors not shown)

Abstract: The Large Array Survey Telescope (LAST) is a wide-field visible-light telescope array designed to explore the variable and transient sky with a high cadence. LAST will be composed of 48, 28-cm f/2.2 telescopes (32 already installed) equipped with full-frame backside-illuminated cooled CMOS detectors. Each telescope provides a field of view (FoV) of 7.4 deg^2 with 1.25 arcsec/pix, while the system… ▽ More The Large Array Survey Telescope (LAST) is a wide-field visible-light telescope array designed to explore the variable and transient sky with a high cadence. LAST will be composed of 48, 28-cm f/2.2 telescopes (32 already installed) equipped with full-frame backside-illuminated cooled CMOS detectors. Each telescope provides a field of view (FoV) of 7.4 deg^2 with 1.25 arcsec/pix, while the system FoV is 355 deg^2 in 2.9 Gpix. The total collecting area of LAST, with 48 telescopes, is equivalent to a 1.9-m telescope. The cost-effectiveness of the system (i.e., probed volume of space per unit time per unit cost) is about an order of magnitude higher than most existing and under-construction sky surveys. The telescopes are mounted on 12 separate mounts, each carrying four telescopes. This provides significant flexibility in operating the system. The first LAST system is under construction in the Israeli Negev Desert, with 32 telescopes already deployed. We present the system overview and performances based on the system commissioning data. The Bp 5-sigma limiting magnitude of a single 28-cm telescope is about 19.6 (21.0), in 20 s (20x20 s). Astrometric two-axes precision (rms) at the bright-end is about 60 (30)\,mas in 20\,s (20x20 s), while absolute photometric calibration, relative to GAIA, provides ~10 millimag accuracy. Relative photometric precision, in a single 20 s (320 s) image, at the bright-end measured over a time scale of about 60 min is about 3 (1) millimag. We discuss the system science goals, data pipelines, and the observatory control system in companion publications. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: Submitted to PASP, 15pp

The $γ$-ray deposition histories of Calcium-rich supernovae

Authors: Amir Sharon, Doron Kushnir

Abstract: Calcium-rich supernovae (Ca-rich SNe) are faint, rapidly evolving transients whose progenitor system is yet to be determined. We derive the $γ$-ray deposition histories of five Ca-rich SNe from the literature in order to place constraints on possible progenitor systems. We find that the $ γ$-ray escape time, $ t_0 $, of the Ca-rich SNe sample is $\approx35$-$65 \,\rm{d}$, within the unoccupied reg… ▽ More Calcium-rich supernovae (Ca-rich SNe) are faint, rapidly evolving transients whose progenitor system is yet to be determined. We derive the $γ$-ray deposition histories of five Ca-rich SNe from the literature in order to place constraints on possible progenitor systems. We find that the $ γ$-ray escape time, $ t_0 $, of the Ca-rich SNe sample is $\approx35$-$65 \,\rm{d}$, within the unoccupied region between Type Ia SNe and stripped envelope supernovae (SESNe). The $ t_0$-$M_\mathrm{Ni56} $ distribution of these SNe, where $M_\mathrm{Ni56}$ is the synthesised $^{56}$Ni mass in the explosion, creates a continuum between the Type Ia and SESNe $ t_0$-$M_\mathrm{Ni56} $ distribution, hinting at a possible connection between all the events. By comparing our results to models from the literature, we were able to determine that helium shell detonation models and core-collapse models of ultra-stripped stars are unlikely to explain Ca-rich SNe, since the gamma-ray escape time in these models is smaller than the observed values. Models that agree with the observed $ t_0$-$M_\mathrm{Ni56} $ distribution are explosions of low mass, $M\approx0.75$-$0.8\,M_\odot $, white dwarfs and core-collapse models of stripped stars with an ejecta mass of $M\approx1$-$3\,M_{\odot}$. △ Less

Submitted 29 December, 2022; originally announced December 2022.

Comments: 12 pages, submitted to MNRAS

Data Collection and Analysis of French Dialects

Authors: Omar Shaur Choudhry, Paul Omara Odida, Joshua Reiner, Keiron Appleyard, Danielle Kushnir, William Toon

Abstract: This paper discusses creating and analysing a new dataset for data mining and text analytics research, contributing to a joint Leeds University research project for the Corpus of National Dialects. This report investigates machine learning classifiers to classify samples of French dialect text across various French-speaking countries. Following the steps of the CRISP-DM methodology, this report ex… ▽ More This paper discusses creating and analysing a new dataset for data mining and text analytics research, contributing to a joint Leeds University research project for the Corpus of National Dialects. This report investigates machine learning classifiers to classify samples of French dialect text across various French-speaking countries. Following the steps of the CRISP-DM methodology, this report explores the data collection process, data quality issues and data conversion for text analysis. Finally, after applying suitable data mining techniques, the evaluation methods, best overall features and classifiers and conclusions are discussed. △ Less

Submitted 1 August, 2022; originally announced August 2022.

Comments: 4 pages plus 1 page for references. 4 figures including 1 image

Strong NIR emission following the long duration GRB 211211A: Dust heating as an alternative to a kilonova

Authors: Eli Waxman, Eran O. Ofek, Doron Kushnir

Abstract: The prolonged near infrared (NIR) emission observed following the long duration GRB 211211A is inconsistent with afterglow emission from the shock driven into the circum-stellar medium (CSM), and with emission from a possible underlying supernova. It has therefore been suggested that the observed NIR flux is the signature of a kilonova -- a radioactive ejecta that is similar to the outcome of the… ▽ More The prolonged near infrared (NIR) emission observed following the long duration GRB 211211A is inconsistent with afterglow emission from the shock driven into the circum-stellar medium (CSM), and with emission from a possible underlying supernova. It has therefore been suggested that the observed NIR flux is the signature of a kilonova -- a radioactive ejecta that is similar to the outcome of the binary neutron star merger GW170817. We propose here an alternative plausible explanation. We show that the NIR flux is consistent with thermal emission from dust, heated by UV and soft X-ray radiation produced by the interaction of the GRB jet plasma with the CSM. This NIR emission was predicted by Waxman & Draine for GRBs residing near or withing massive molecular clouds. The dust NIR emission scenario is consistent with a GRB at $z\lesssim1$. Inspection of the environment of GRB 211211A suggests that there are at least two host-galaxy candidates, one at $z=0.076$ and the other at $z=0.459$. The $z=0.459$ possibility is also consistent with the non-detection of a supernova signature in the light curve of the GRB afterglow, and with a typical GRB $γ$-ray energy for the fluence of GRB 211211A. △ Less

Submitted 8 February, 2023; v1 submitted 21 June, 2022; originally announced June 2022.

Comments: Expanded discussion of ionizing radiation and evidence for residual extinction

Stress overshoot, hysteresis and Bauschinger effect in sheared dense colloidal suspensions

Authors: Dmytro Kushnir, Céline Ruscher, Eckhard Bartsch, Fabrice Thalmann, Pascal Hébraud

Abstract: The mechanical non-linear response of dense Brownian suspensions of polymer gel particles is studied experimentally and by means of numerical simulations. It is shown that the response to the application of a constant shear rate depends on the previous history of the suspension. When the flow starts from a suspension at rest, it exhibits an elastic response followed by a stress overshoot and then… ▽ More The mechanical non-linear response of dense Brownian suspensions of polymer gel particles is studied experimentally and by means of numerical simulations. It is shown that the response to the application of a constant shear rate depends on the previous history of the suspension. When the flow starts from a suspension at rest, it exhibits an elastic response followed by a stress overshoot and then a plastic flow regime. Conversely, after flow reversal, the stress overshoot does not occur, and the apparent elastic modulus is reduced while numerical simulations reveal that the anisotropy of the local microstructure is delayed relative to the macroscopic stress. △ Less

Submitted 23 February, 2022; originally announced February 2022.

Comments: Submitted

Confronting double-detonation sub-Chandrasekhar models with the low-luminosity suppression of Type Ia supernovae

Authors: Arka Ghosh, Doron Kushnir

Abstract: Type Ia supernovae (SNe Ia) are likely the thermonuclear explosions of carbon-oxygen (CO) white-dwarf (WD) stars, but their progenitor systems remain elusive. Recently, Sharon & Kushnir (2022) used The Zwicky Transient Facility Bright Transient Survey to construct a synthesized $^{56}$Ni mass, $M_\text{Ni56}$, distribution of SNe Ia. They found that the rate of low-luminosity (… ▽ More Type Ia supernovae (SNe Ia) are likely the thermonuclear explosions of carbon-oxygen (CO) white-dwarf (WD) stars, but their progenitor systems remain elusive. Recently, Sharon & Kushnir (2022) used The Zwicky Transient Facility Bright Transient Survey to construct a synthesized $^{56}$Ni mass, $M_\text{Ni56}$, distribution of SNe Ia. They found that the rate of low-luminosity ($M_\text{Ni56}\approx0.15\,M_{\odot}$) SNe Ia is lower by a factor of $\sim10$ than the more common $M_\text{Ni56}\approx0.7\,M_{\odot}$ events. We here show that in order for the double-detonation model (DDM, in which a propagating thermonuclear detonation wave, TNDW, within a thin helium shell surrounding a sub-Chandrasekhar mass CO core triggers a TNDW within the core) to explain this low-luminosity suppression, the probability of a low-mass ($\approx0.85\,M_{\odot}$) WD explosion should be $\sim100$-fold lower than that of a high-mass ($\approx1.05\,M_{\odot}$) WD. One possible explanation is that the ignition of low-mass CO cores is somehow suppressed. We use accurate one-dimensional numerical simulations to show that if a TNDW is able to propagate within the helium shell, then the ignition within the CO core is guaranteed (resolved here for the first time in a full-star simulation), even for $0.7\,M_{\odot}$ WDs, providing no natural explanation for the low-luminosity suppression. DDM could explain the low-luminosity suppression if the mass distribution of primary WDs in close binaries is dramatically different from the field distribution; if the Helium shell ignition probability is suppressed for low-mass WDs; or if multidimensional perturbations significantly change our results. △ Less

Submitted 22 February, 2023; v1 submitted 25 October, 2021; originally announced October 2021.

Comments: 7 pages, 3 figures. Revised following referee report

AT2018lqh and the nature of the emerging population of day-scale duration optical transients

Authors: E. O. Ofek, S. M. Adams, E. Waxman, A. Sharon, D. Kushnir, A. Horesh, A. Ho, M. M. Kasliwal, O. Yaron, A. Gal-Yam, S. R. Kulkarni, E. Bellm, F. Masci, D. Shupe, R. Dekany, M. Graham, R. Riddle, D. Duev, I. Andreoni, A. Mahabal, A. Drake

Abstract: We report on the discovery of AT2018lqh (ZTF18abfzgpl) -- a rapidly-evolving extra-galactic transient in a star-forming host at 242 Mpc. The transient g-band light curve's duration above half-maximum light is about 2.1 days, where 0.4/1.7 days are spent on the rise/decay, respectively. The estimated bolometric light curve of this object peaked at about 7x10^42 erg/s -- roughly seven times brighter… ▽ More We report on the discovery of AT2018lqh (ZTF18abfzgpl) -- a rapidly-evolving extra-galactic transient in a star-forming host at 242 Mpc. The transient g-band light curve's duration above half-maximum light is about 2.1 days, where 0.4/1.7 days are spent on the rise/decay, respectively. The estimated bolometric light curve of this object peaked at about 7x10^42 erg/s -- roughly seven times brighter than AT2017gfo. We show that this event can be explained by an explosion with a fast (v~0.08 c) low-mass (~0.07 Msun) ejecta, composed mostly of radioactive elements. For example, ejecta dominated by Ni-56 with a time scale of t_0=1.6 days for the ejecta to become optically thin for gamma-rays fits the data well. Such a scenario requires burning at densities that are typically found in the envelopes of neutron stars or the cores of white dwarfs. A combination of circumstellar material (CSM) interaction power at early times and shock cooling at late times is consistent with the photometric observations, but the observed spectrum of the event may pose some challenges for this scenario. The observations are not consistent with a shock breakout from a stellar envelope, while a model involving a low-mass ejecta ramming into low-mass CSM cannot explain both the early- and late-time observations. △ Less

Submitted 22 September, 2021; originally announced September 2021.

Comments: 10 pages, ApJ in press

The ZTF-BTS Type Ia supernovae luminosity function is consistent with a single progenitor channel for the explosions

Authors: Amir Sharon, Doron Kushnir

Abstract: We construct the Type Ia supernovae (SNe Ia) luminosity function (LF) using the Zwicky Transient Facility Bright Transient Survey (BTS) catalogue. While this magnitude-limited survey has an unprecedented number of objects, it suffers from large distance uncertainties and lacks an estimation of host extinction. We bypass these issues by calculating the intrinsic luminosities from the shape paramete… ▽ More We construct the Type Ia supernovae (SNe Ia) luminosity function (LF) using the Zwicky Transient Facility Bright Transient Survey (BTS) catalogue. While this magnitude-limited survey has an unprecedented number of objects, it suffers from large distance uncertainties and lacks an estimation of host extinction. We bypass these issues by calculating the intrinsic luminosities from the shape parameters of the light curve's $ g $ and $ r $ bands, with the luminosities calibrated from the well observed SNe Ia sample of the Carnegie Supernova Project, allowing us to construct, for the first time, the intrinsic LF of SNe Ia. We then use a novel tight relation between the color stretch and the synthesized $^{56}$Ni mass, $M_\mathrm{Ni56}$, to determine the $M_\mathrm{Ni56}$ distribution of SNe Ia. We find that the LFs are unimodal, with their peaks in line with previous results, but have a much lower rate of dim events and luminous events. We show that the features on top of the unimodal LF-derived distributions are all compatible with statistical noise, consistent with a single progenitor channel for the explosions. We further derive, for the first time, the SNe Ia distribution of host galaxy extinction, and find a mean selective extinction of $E(B-V)\approx0.1$ and a non-negligible fraction with large, $ >1\,\text{mag} $, extinction in the optical bands. The high extinction is typical for luminous SNe, supporting their young population origin. △ Less

Submitted 13 December, 2021; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: 24 pages, 20 figures. Accepted for publication in MNRAS

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

Authors: Amir Sharon, Doron Kushnir

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

Submitted 24 August, 2020; originally announced August 2020.

Comments: 4 pages, 1 figure. arXiv admin note: substantial text overlap with arXiv:2008.08592

Sub-Chandrasekhar-mass detonations are in tension with the observed $t_0-M_\text{Ni56}$ relation of type Ia supernovae

Authors: Doron Kushnir, Nahliel Wygoda, Amir Sharon

Abstract: Type Ia supernovae (SNe Ia) are likely the thermonuclear explosions of carbon-oxygen (CO) white-dwarf (WD) stars, but their progenitor systems remain elusive. Recent studies have suggested that a propagating detonation within a thin helium shell surrounding a sub-Chandrasekhar mass CO core can subsequently trigger a detonation within the core (the double-detonation model, DDM). The outcome of this… ▽ More Type Ia supernovae (SNe Ia) are likely the thermonuclear explosions of carbon-oxygen (CO) white-dwarf (WD) stars, but their progenitor systems remain elusive. Recent studies have suggested that a propagating detonation within a thin helium shell surrounding a sub-Chandrasekhar mass CO core can subsequently trigger a detonation within the core (the double-detonation model, DDM). The outcome of this explosion is similar to a central ignition of a sub-Chandrasekhar mass CO WD (SCD). While SCD is consistent with some observational properties of SNe Ia, several computational challenges prohibit a robust comparison to the observations. We focus on the observed $t_0-M_\text{Ni56}$ relation, where $t_0$ (the $γ$-rays' escape time from the ejecta) is positively correlated with $M_\text{Ni56}$ (the synthesized $^{56}$Ni mass). We apply our recently developed numerical scheme to calculate SCD and show that the calculated $t_0-M_\text{Ni56}$ relation, which does not require radiation transfer calculations, converges to an accuracy of a few percent. We find a clear tension between our calculations and the observed $t_0-M_\text{Ni56}$ relation. SCD predicts an anti-correlation between $t_0$ and $M_\text{Ni56}$, with $t_0\approx30\,\textrm{day}$ for luminous ($M_\text{Ni56}\gtrsim0.5\,M_{\odot}$) SNe Ia, while the observed $t_0$ is in the range of $35-45\,\textrm{day}$. We show that this tension is larger than the uncertainty of the results, and that it exists in all previous studies of the problem. Our results hint that more complicated models are required, but we argue that DDM is unlikely to resolve the tension with the observations. △ Less

Submitted 26 October, 2020; v1 submitted 19 August, 2020; originally announced August 2020.

Comments: 24 pages, 20 figures. Main results in the first 3 pages with 1 figure. Minor revisions following referee's report

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

Authors: Amir Sharon, Doron Kushnir

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

Submitted 23 July, 2020; v1 submitted 15 April, 2020; originally announced April 2020.

Comments: 29 pages, 12 figures, accepted to MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 496, Issue 4, August 2020, Pages 4517-4545

Diffusion-based Deep Active Learning

Authors: Dan Kushnir, Luca Venturi

Abstract: The remarkable performance of deep neural networks depends on the availability of massive labeled data. To alleviate the load of data annotation, active deep learning aims to select a minimal set of training points to be labelled which yields maximal model accuracy. Most existing approaches implement either an `exploration'-type selection criterion, which aims at exploring the joint distribution o… ▽ More The remarkable performance of deep neural networks depends on the availability of massive labeled data. To alleviate the load of data annotation, active deep learning aims to select a minimal set of training points to be labelled which yields maximal model accuracy. Most existing approaches implement either an `exploration'-type selection criterion, which aims at exploring the joint distribution of data and labels, or a `refinement'-type criterion which aims at localizing the detected decision boundaries. We propose a versatile and efficient criterion that automatically switches from exploration to refinement when the distribution has been sufficiently mapped. Our criterion relies on a process of diffusing the existing label information over a graph constructed from the hidden representation of the data set as provided by the neural network. This graph representation captures the intrinsic geometry of the approximated labeling function. The diffusion-based criterion is shown to be advantageous as it outperforms existing criteria for deep active learning. △ Less

Submitted 23 March, 2020; originally announced March 2020.

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

Authors: Doron Kushnir, Eli Waxman

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

Submitted 22 March, 2020; v1 submitted 27 January, 2020; originally announced January 2020.

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

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

Authors: Doron Kushnir, Boaz Katz

Abstract: 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 unders… ▽ More 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. △ Less

Submitted 20 March, 2020; v1 submitted 12 December, 2019; originally announced December 2019.

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

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

Authors: Doron Kushnir, Boaz Katz

Abstract: 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:… ▽ More 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/(βt_p)^2\int_0 ^{βt_{p}} t'Q(t')dt'$, where $Q(t)$ is the heating rate of the ejecta, and $β$ 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 $β$ 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. △ Less

Submitted 16 October, 2019; originally announced October 2019.

Comments: 3 pages, 1 figure

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

Authors: Doron Kushnir, Boaz Katz

Abstract: 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 spati… ▽ More 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. △ Less

Submitted 12 December, 2019; v1 submitted 22 April, 2019; originally announced April 2019.

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

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

Authors: Amir Sharon, Doron Kushnir

Abstract: 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 interactio… ▽ More 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. △ Less

Submitted 28 April, 2020; v1 submitted 17 April, 2019; originally announced April 2019.

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

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

Authors: Juna A. Kollmeier, Ping Chen, Subo Dong, Nidia Morrell, M. M. Phillips, Doron Kushnir, J. L. Prieto, Anthony L. Piro, Joshua D. Simon

Abstract: 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 prominen… ▽ More 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$α$ emission line of $L_{{\rm H}α}=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$α$ emission ($L_{{\rm H}α}\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$α$ 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$α$ emission could be as high as $\sim 10\%$ level among sub-luminous SNe Ia. △ Less

Submitted 6 February, 2019; originally announced February 2019.

Comments: Submitted to MNRAS

Late time kilonova light curves and implications to GW 170817

Authors: Eli Waxman, Eran O. Ofek, Doron Kushnir

Abstract: 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 lum… ▽ More 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. △ Less

Submitted 26 February, 2019; v1 submitted 4 February, 2019; originally announced February 2019.

Comments: Improved treatment of synchrotron subtraction

The structure of detonation waves in supernovae revisited

Authors: Doron Kushnir

Abstract: 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 plasm… ▽ More 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. △ Less

Submitted 28 November, 2018; v1 submitted 4 June, 2018; originally announced June 2018.

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

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

Authors: Doron Kushnir, Eli Waxman, Andrey I. Chugunov

Abstract: 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… ▽ More 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 $+ν_{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^{+}+ν_{e}$ is enhanced. This effect may be measured in the near future. △ Less

Submitted 2 April, 2019; v1 submitted 22 May, 2018; originally announced May 2018.

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

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

Authors: Subo Dong, Boaz Katz, Juna A. Kollmeier, Doron Kushnir, N. Elias-Rosa, Subhash Bose, Nidia Morrell, J. L. Prieto, Ping Chen, C. S. Kochanek, G. M. Brandt, T. W. -S. Holoien, Avishay Gal-Yam, Antonia Morales-Garoffolo, Stuart Parker, M. M. Phillips, Anthony L. Piro, B. J. Shappee, Joshua D. Simon, K. Z. Stanek

Abstract: 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 compare… ▽ More 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. △ Less

Submitted 31 May, 2018; v1 submitted 30 April, 2018; originally announced May 2018.

Comments: Minor Changes. Accepted by MNRAS Letter

Journal ref: Published in Mon.Not.Roy.Astron.Soc. 479 (2018) no.1, L70-L75

Active Community Detection with Maximal Expected Model Change

Authors: Dan Kushnir, Benjamin Mirabelli

Abstract: We present a novel active learning algorithm for community detection on networks. Our proposed algorithm uses a Maximal Expected Model Change (MEMC) criterion for querying network nodes label assignments. MEMC detects nodes that maximally change the community assignment likelihood model following a query. Our method is inspired by detection in the benchmark Stochastic Block Model (SBM), where we p… ▽ More We present a novel active learning algorithm for community detection on networks. Our proposed algorithm uses a Maximal Expected Model Change (MEMC) criterion for querying network nodes label assignments. MEMC detects nodes that maximally change the community assignment likelihood model following a query. Our method is inspired by detection in the benchmark Stochastic Block Model (SBM), where we provide sample complexity analysis and empirical study with SBM and real network data for binary as well as for the multi-class settings. The analysis also covers the most challenging case of sparse degree and below-detection-threshold SBMs, where we observe a super-linear error reduction. MEMC is shown to be superior to the random selection baseline and other state-of-the-art active learners. △ Less

Submitted 20 March, 2020; v1 submitted 10 January, 2018; originally announced January 2018.

Linear Time Clustering for High Dimensional Mixtures of Gaussian Clouds

Authors: Dan Kushnir, Shirin Jalali, Iraj Saniee

Abstract: Clustering mixtures of Gaussian distributions is a fundamental and challenging problem that is ubiquitous in various high-dimensional data processing tasks. While state-of-the-art work on learning Gaussian mixture models has focused primarily on improving separation bounds and their generalization to arbitrary classes of mixture models, less emphasis has been paid to practical computational effici… ▽ More Clustering mixtures of Gaussian distributions is a fundamental and challenging problem that is ubiquitous in various high-dimensional data processing tasks. While state-of-the-art work on learning Gaussian mixture models has focused primarily on improving separation bounds and their generalization to arbitrary classes of mixture models, less emphasis has been paid to practical computational efficiency of the proposed solutions. In this paper, we propose a novel and highly efficient clustering algorithm for $n$ points drawn from a mixture of two arbitrary Gaussian distributions in $\mathbb{R}^p$. The algorithm involves performing random 1-dimensional projections until a direction is found that yields a user-specified clustering error $e$. For a 1-dimensional separation parameter $γ$ satisfying $γ=Q^{-1}(e)$, the expected number of such projections is shown to be bounded by $o(\ln p)$, when $γ$ satisfies $γ\leq c\sqrt{\ln{\ln{p}}}$, with $c$ as the separability parameter of the two Gaussians in $\mathbb{R}^p$. Consequently, the expected overall running time of the algorithm is linear in $n$ and quasi-linear in $p$ at $o(\ln{p})O(np)$, and the sample complexity is independent of $p$. This result stands in contrast to prior works which provide polynomial, with at-best quadratic, running time in $p$ and $n$. We show that our bound on the expected number of 1-dimensional projections extends to the case of three or more Gaussian components, and we present a generalization of our results to mixture distributions beyond the Gaussian model. △ Less

Submitted 1 March, 2018; v1 submitted 19 December, 2017; originally announced December 2017.

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

Authors: Eli Waxman, Eran Ofek, Doron Kushnir, Avishay Gal-Yam

Abstract: 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$ componen… ▽ More 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, {$κ<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 $κ_ν\approx0.1$\,cm$^2$\,g$^{-1}$ at 1 to $2μ$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. △ Less

Submitted 4 September, 2018; v1 submitted 27 November, 2017; originally announced November 2017.

Comments: Accepted to MNRAS. Expanded explanations, improved figures

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

Authors: Matias Zaldarriaga, Doron Kushnir, Juna A. Kollmeier

Abstract: 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 sp… ▽ More 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. △ Less

Submitted 2 February, 2017; originally announced February 2017.

Comments: 5 pages, 4 figures

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

Authors: Doron Kushnir, Matias Zaldarriaga, Juna A. Kollmeier, Roni Waldman

Abstract: 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 explor… ▽ More 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. △ Less

Submitted 27 July, 2016; v1 submitted 12 May, 2016; originally announced May 2016.

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

Dynamical Tides Reexpressed

Authors: Doron Kushnir, Matias Zaldarriaga, Juna A. Kollmeier, Roni Waldman

Abstract: 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 ve… ▽ More 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, $β_{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 $β_{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. △ Less

Submitted 13 September, 2017; v1 submitted 12 May, 2016; originally announced May 2016.

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

Economic and Technological Complexity: A Model Study of Indicators of Knowledge-based Innovation Systems

Authors: Inga Ivanova, Oivind Strand, Duncan Kushnir, Loet Leydesdorff

Abstract: The Economic Complexity Index (ECI; Hidalgo & Hausmann, 2009) measures the complexity of national economies in terms of product groups. Analogously to ECI, a Patent Complexity Index (PatCI) can be developed on the basis of a matrix of nations versus patent classes. Using linear algebra, the three dimensions: countries, product groups, and patent classes can be combined into a measure of "Triple He… ▽ More The Economic Complexity Index (ECI; Hidalgo & Hausmann, 2009) measures the complexity of national economies in terms of product groups. Analogously to ECI, a Patent Complexity Index (PatCI) can be developed on the basis of a matrix of nations versus patent classes. Using linear algebra, the three dimensions: countries, product groups, and patent classes can be combined into a measure of "Triple Helix" complexity (THCI) including the trilateral interaction terms between knowledge production, wealth generation, and (national) control. THCI can be expected to capture the extent of systems integration between the global dynamics of markets (ECI) and technologies (PatCI) in each national system of innovation. We measure ECI, PatCI, and THCI during the period 2000-2014 for the 34 OECD member states, the BRICS countries, and a group of emerging and affiliated economies (Argentina, Hong Kong, Indonesia, Malaysia, Romania, and Singapore). The three complexity indicators are correlated between themselves; but the correlations with GDP per capita are virtually absent. Of the world's major economies, Japan scores highest on all three indicators, while China has been increasingly successful in combining economic and technological complexity. We could not reproduce the correlation between ECI and average income that has been central to the argument about the fruitfulness of the economic complexity approach. △ Less

Submitted 7 December, 2016; v1 submitted 7 February, 2016; originally announced February 2016.

Journal ref: Technological Forecasting and Social Change 120 (July 2017) 77-89

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

Authors: Kfir Blum, Doron Kushnir

Abstract: 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… ▽ More 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. △ Less

Submitted 13 January, 2016; originally announced January 2016.

Comments: 11 pages, 6 figures

The Globalization of Academic Entrepreneurship? The Recent Growth (2009-2014) in University Patenting Decomposed

Authors: Loet Leydesdorff, Henry Etzkowitz, Duncan Kushnir

Abstract: The contribution of academia to US patents has become increasingly global. Following a pause, with a relatively flat rate, from 1998 to 2008, the long-term trend of university patenting rising as a share of all patenting has resumed, driven by the internationalization of academic entrepreneurship and the persistence of US university technology transfer. We disaggregate this recent growth in univer… ▽ More The contribution of academia to US patents has become increasingly global. Following a pause, with a relatively flat rate, from 1998 to 2008, the long-term trend of university patenting rising as a share of all patenting has resumed, driven by the internationalization of academic entrepreneurship and the persistence of US university technology transfer. We disaggregate this recent growth in university patenting at the US Patent and Trademark Organization (USPTO) in terms of nations and patent classes. Foreign patenting in the US has almost doubled during the period 2009-2014, mainly due to patenting by universities in Taiwan, Korea, China, and Japan. These nations compete with the US in terms of patent portfolios, whereas most European countries--with the exception of the UK--have more specific portfolios, mainly in the bio-medical fields. In the case of China, Tsinghua University holds 63% of the university patents in USPTO, followed by King Fahd University with 55.2% of the national portfolio. △ Less

Submitted 14 December, 2015; originally announced December 2015.

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

Authors: Cole Holcomb, Doron Kushnir

Abstract: 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 employe… ▽ More 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. △ Less

Submitted 26 October, 2015; originally announced October 2015.

Comments: 5 Pages, 2 Figure, 3 Tables

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

Authors: Doron Kushnir

Abstract: 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… ▽ More 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. △ Less

Submitted 8 June, 2015; originally announced June 2015.

Comments: 4 pages, 3 figures, 2 tables

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

Authors: Doron Kushnir

Abstract: 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 thermon… ▽ More 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. △ Less

Submitted 10 February, 2015; originally announced February 2015.

Comments: 4 pages, 5 figures, 3 tables

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

Authors: Doron Kushnir, Boaz Katz

Abstract: 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) calculation… ▽ More 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. △ Less

Submitted 14 September, 2015; v1 submitted 2 December, 2014; originally announced December 2014.

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

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

Authors: Boaz Katz, Subo Dong, Doron Kushnir

Abstract: 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 fin… ▽ More 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. △ Less

Submitted 27 February, 2014; originally announced February 2014.

Comments: 3 pages, 2 figures

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

Authors: Subo Dong, Boaz Katz, Doron Kushnir, Jose L. Prieto

Abstract: 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 velo… ▽ More 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 (Δ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. △ Less

Submitted 1 September, 2015; v1 submitted 14 January, 2014; originally announced January 2014.

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

Journal ref: MNRAS (November 21, 2015) 454 (1): L61-L65

Early Hydrodynamic Evolution of a Stellar Collision

Authors: Doron Kushnir, Boaz Katz

Abstract: 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)… ▽ More 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). △ Less

Submitted 5 November, 2013; originally announced November 2013.

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

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

Authors: Gilad Rave, Doron Kushnir, Eli Waxman

Abstract: 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 w… ▽ More 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. △ Less

Submitted 18 September, 2013; v1 submitted 15 April, 2013; originally announced April 2013.

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

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

Authors: Doron Kushnir, Boaz Katz, Subo Dong, Eli Livne, Rodrigo Fernández

Abstract: 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 whic… ▽ More 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. △ Less

Submitted 19 November, 2013; v1 submitted 5 March, 2013; originally announced March 2013.

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

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

Authors: Boaz Katz, Doron Kushnir, Subo Dong

Abstract: 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 appr… ▽ More 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. △ Less

Submitted 28 January, 2013; originally announced January 2013.

Comments: 1 paragraph