Showing 1–50 of 183 results for author: Waxman, E

Shock cooling emission from explosions of massive stars: III. Blue Super Giants

Authors: J. Morag, N. Sapir, E. Waxman

Abstract: Light emission in the first hours and days following core-collapse supernovae is dominated by the escape of photons from the expanding shock-heated envelope. In preceding papers, we provided a simple analytic description of the time-dependent luminosity, $L$, and color temperature, $T_{\rm col}$, as well as of the small ($\simeq10\%$) deviations of the spectrum from blackbody at low frequencies,… ▽ More Light emission in the first hours and days following core-collapse supernovae is dominated by the escape of photons from the expanding shock-heated envelope. In preceding papers, we provided a simple analytic description of the time-dependent luminosity, $L$, and color temperature, $T_{\rm col}$, as well as of the small ($\simeq10\%$) deviations of the spectrum from blackbody at low frequencies, $hν< 3T_{\rm col}$, and of `line dampening' at $hν> 3T_{\rm col}$, for explosions of red supergiants (RSGs) with convective polytropic envelopes (without significant circum-stellar medium). Here, we extend our work to provide similar analytic formulae for explosions of blue supergiants with radiative polytropic envelopes. The analytic formulae are calibrated against a large set of spherically symmetric multi-group (frequency-dependent) calculations for a wide range of progenitor parameters (mass, radius, core/envelope mass ratios) and explosion energies. In these calculations we use the opacity tables we constructed (and made publicly available), that include the contributions of bound-bound and bound-free transitions. They reproduce the numeric $L$ and $T_{\rm col}$ to within 10\% and 5\% accuracy, and the spectral energy distribution to within $\sim20-40\%$. The accuracy is similar to that achieved for RSG explosions. △ Less

Submitted 17 July, 2024; originally announced July 2024.

Comments: 12 pages, 10 figures

Preliminary results of the Single Event Effect testing for the ULTRASAT sensors

Authors: Vlad Dumitru Berlea, Arooj Asif, Merlin F. Barschke, David Berge, Juan Maria Haces Crespo, Gianluca Giavitto, Shashank Kumar, Andrea Porelli, Nicola de Simone, Jason Watson, Steven Worm, Francesco Zappon, Adi Birman, Shay Alfassi, Amos Feningstein, Eli Waxman, Udi Netzer, Tuvia Liran, Ofer Lapid, Viktor M. Algranatti, Yossi Schvartzvald

Abstract: ULTRASAT (ULtra-violet TRansient Astronomy SATellite) is a wide-angle space telescope that will perform a deep time-resolved all-sky survey in the near-ultraviolet (NUV) spectrum. The science objectives are the detection of counterparts to short-lived transient astronomical events such as gravitational wave sources and supernovae. The mission is led by the Weizmann Institute of Science and is plan… ▽ More ULTRASAT (ULtra-violet TRansient Astronomy SATellite) is a wide-angle space telescope that will perform a deep time-resolved all-sky survey in the near-ultraviolet (NUV) spectrum. The science objectives are the detection of counterparts to short-lived transient astronomical events such as gravitational wave sources and supernovae. The mission is led by the Weizmann Institute of Science and is planned for launch in 2026 in collaboration with the Israeli Space Agency and NASA. DESY will provide the UV camera, composed by the detector assembly located in the telescope focal plane and the remote electronics unit. The camera is composed out of four back-metallized CMOS Image Sensors (CIS) manufactured in the 4T, dual gain Tower process. As part of the radiation qualification of the camera, Single Event Effect (SEE) testing has been performed by irradiating the sensor with heavy ions at the RADEF, Jyvaskyla facility. Preliminary results of both Single Event Upset (SEU) and Single Event Latch-up (SEL) occurrence rate in the sensor are presented. Additionally, an in-orbit SEE rate simulation has been performed in order to gain preliminary knowledge about the expected effect of SEE on the mission. △ Less

Submitted 4 July, 2024; originally announced July 2024.

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

The thermalization of $γ$-rays in radioactive expanding ejecta: A simple model and its application for Kilonovae and Ia SNe

Authors: Or Guttman, Ben Shenhar, Arnab Sarkar, Eli Waxman

Abstract: A semi-analytic approximation is derived for the time-dependent fraction $f_γ(t)$ of the energy deposited by radioactive decay $γ$-rays in a homologously expanding plasma of general structure. An analytic approximation is given for spherically symmetric plasma distributions. Applied to Kilonovae (KNe) associated with neutron stars mergers and Type Ia supernovae, our semi-analytic and analytic appr… ▽ More A semi-analytic approximation is derived for the time-dependent fraction $f_γ(t)$ of the energy deposited by radioactive decay $γ$-rays in a homologously expanding plasma of general structure. An analytic approximation is given for spherically symmetric plasma distributions. Applied to Kilonovae (KNe) associated with neutron stars mergers and Type Ia supernovae, our semi-analytic and analytic approximations reproduce, with a few percent and 10% accuracy, respectively, the energy deposition rates, $\dot{Q}_\text{dep}$, obtained in numeric Monte Carlo calculations. The time $t_γ$ beyond which $γ$-ray deposition is inefficient is determined by an effective frequency-independent $γ$-ray opacity $κ_{γ,\text{eff}}$, $t_γ= \sqrt{κ_{γ,\text{eff}}\langleΣ\rangle t^2}$, where $\langleΣ\rangle\propto t^{-2}$ is the average plasma column density. For $β$-decay dominated energy release, $κ_{γ,\text{eff}}$ is typically close to the effective Compton scattering opacity, $κ_{γ,\text{eff}} \approx 0.025~{\rm {cm}^{2}\,g^{-1}}$ with a weak dependence on composition. For KNe, $κ_{γ,\text{eff}}$ depends mainly on the initial electron fraction $Y_e$, $κ_{γ,\text{eff}} \approx 0.03(0.05)~{\rm {cm}^{2}\,g^{-1}}$ for $Y_e \gtrsim (\lesssim) 0.25$ (in contrast with earlier work that found $κ_{γ,\text{eff}}$ larger by 1-2 orders of magnitude for low $Y_e$), and is insensitive to the (large) nuclear physics uncertainties. Determining $t_γ$ from observations will therefore measure the ejecta $\langleΣ\rangle t^2$, providing a stringent test of models. For $\langleΣ\rangle t^2=2\times10^{11}~{\rm g\,{cm}^{-2}\,s^2}$, a typical value expected for KNe, $t_γ\approx1$ d. △ Less

Submitted 13 March, 2024; originally announced March 2024.

Comments: 18 pages, 15 figures, 1 table. submitted to MNRAS

An Analytic Description of Electron Thermalization in Kilonovae Ejecta

Authors: Ben Shenhar, Or Guttman, Eli Waxman

Abstract: A simple analytic description is provided of the rate of energy deposition by $β$-decay electrons in the homologously expanding radioactive plasma ejected in neutron star mergers, valid for a wide range of ejecta parameters -- initial entropy, electron fraction $\{s_0,Y_e\}$ and density $ρt^3$. The formulae are derived using detailed numerical calculations following the time-dependent composition… ▽ More A simple analytic description is provided of the rate of energy deposition by $β$-decay electrons in the homologously expanding radioactive plasma ejected in neutron star mergers, valid for a wide range of ejecta parameters -- initial entropy, electron fraction $\{s_0,Y_e\}$ and density $ρt^3$. The formulae are derived using detailed numerical calculations following the time-dependent composition and $β$-decay emission spectra (including the effect of delayed deposition). The deposition efficiency depends mainly on $ρt^3$ and only weakly on $\{s_0,Y_e\}$. The time $t_e$ at which the ratio between the rates of electron energy deposition and energy production drops to $1-e^{-1}$, is given by $t_e=t_{0e}\Big(\frac{ρt^3}{0.5(ρt^3)_0}\Big)^a$, where $(ρt^3)_0=\frac{0.05M_{\odot}}{4π(0.2c)^3}$, $t_{0e}(s_0,Y_e)\approx17$ days and $0.4\le a(s_0,Y_e)\le0.5$. The fractional uncertainty in $t_e$ due to nuclear physics uncertainties is $\approx10\%$. The result $a\le0.5$ reflects the fact that the characteristic $β$-decay electron energies do not decrease with time (largely due to "inverted decay chains" in which a slowly-decaying isotope decays to a rapidly-decaying isotope with higher end-point energy). We provide an analytic approximation for the time-dependent electron energy deposition rate, reproducing the numerical results to better than $50\%$ (typically $<30\%$, well within the energy production rate uncertainty due to nuclear physics uncertainties) over a 3-4 orders-of-magnitude deposition rate decrease with time. Our results may be easily incorporated in calculations of kilonovae light curves (with general density and composition structures), eliminating the need to numerically follow the time-dependent electron spectra. Identifying $t_e$, e.g. in the bolometric light curve, will constrain the (properly averaged) ejecta $ρt^3$. △ Less

Submitted 13 March, 2024; originally announced March 2024.

Comments: 10 pages, 10 figures, 1 table. submitted to MNRAS

Non-thermal emission from mildly relativistic dynamical ejecta of neutron star mergers: spectrum and sky image

Authors: Gilad Sadeh, Noya Linder, Eli Waxman

Abstract: Binary neutron star mergers are expected to produce fast dynamical ejecta, with mildly relativistic velocities extending to $β=v/c>0.6$. In a preceding paper, we derived an analytic description of the time-dependent radio to X-ray synchrotron flux produced by collisionless shocks driven by such fast ejecta into the interstellar medium, for spherical ejecta with broken power-law mass (or energy) di… ▽ More Binary neutron star mergers are expected to produce fast dynamical ejecta, with mildly relativistic velocities extending to $β=v/c>0.6$. In a preceding paper, we derived an analytic description of the time-dependent radio to X-ray synchrotron flux produced by collisionless shocks driven by such fast ejecta into the interstellar medium, for spherical ejecta with broken power-law mass (or energy) distributions, $M(>γβ)\propto(γβ)^{-s}$ with $s=s_\text{KN}$ at $γβ<γ_0β_0$ and $s=s_\text{ft}$ at $γβ>γ_0β_0$ (where $γ$ is the Lorentz factor). Here, we extend our analysis and provide analytic expressions for the self-absorption frequency, the cooling frequency, and the observed angular size of the emitting region (which appears as a ring in the sky). For parameter values characteristic of merger calculation results -- a "shallow" mass distribution, $1<s_\text{KN}<3$, for the bulk of the ejecta (at $γβ\approx0.2$), and a steep, $s_\text{ft}>5$, "fast tail" mass distribution -- the analytic results reproduce well (to tens of percent accuracy) the results of detailed numeric calculations, a significant improvement over earlier order-of-magnitude estimates (based on extrapolations of results valid for $γβ\ll1$). △ Less

Submitted 24 May, 2024; v1 submitted 11 March, 2024; originally announced March 2024.

Comments: 8 pages, 7 figures, accepted to MNRAS

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

The Early Ultraviolet Light-Curves of Type II Supernovae and the Radii of Their Progenitor Stars

Authors: Ido Irani, Jonathan Morag, Avishay Gal-Yam, Eli Waxman, Steve Schulze, Jesper Sollerman, K-Ryan Hinds, Daniel A. Perley, Ping Chen, Nora L. Strotjohann, Ofer Yaron, Erez A. Zimmerman, Rachel Bruch, Eran O. Ofek, Maayane T. Soumagnac, Yi Yang, Steven L. Groom, Frank J. Masci, Reed Riddle, Eric C. Bellm, David Hale

Abstract: We present a sample of 34 normal SNe II detected with the Zwicky Transient Facility, with multi-band UV light-curves starting at $t \leq 4$ days after explosion, as well as X-ray detections and upper limits. We characterize the early UV-optical colors and provide prescriptions for empirical host-extinction corrections. We show that the $t > 2\,$days UV-optical colors and the blackbody evolution of… ▽ More We present a sample of 34 normal SNe II detected with the Zwicky Transient Facility, with multi-band UV light-curves starting at $t \leq 4$ days after explosion, as well as X-ray detections and upper limits. We characterize the early UV-optical colors and provide prescriptions for empirical host-extinction corrections. We show that the $t > 2\,$days UV-optical colors and the blackbody evolution of the sample are consistent with the predictions of spherical phase shock-cooling (SC), independently of the presence of `flash ionization" features. We present a framework for fitting SC models which can reproduce the parameters of a set of multi-group simulations without a significant bias up to 20% in radius and velocity. Observations of about half of the SNe II in the sample are well-fit by models with breakout radii $<10^{14}\,$cm. The other half are typically more luminous, with observations from day 1 onward that are better fit by a model with a large $>10^{14}\,$cm breakout radius. However, these fits predict an early rise during the first day that is too slow. We suggest these large-breakout events are explosions of stars with an inflated envelope or a confined CSM with a steep density profile, at which breakout occurs. Using the X-ray data, we derive constraints on the extended ($\sim10^{15}$ cm) CSM density independent of spectral modeling, and find most SNe II progenitors lose $<10^{-4} M_{\odot}\, \rm yr^{-1}$ a few years before explosion. This provides independent evidence the CSM around many SNe II progenitors is confined. We show that the overall observed breakout radius distribution is skewed to higher radii due to a luminosity bias. We argue that the $66^{+11}_{-22}\%$ of red supergiants (RSG) explode as SNe II with breakout radii consistent with the observed distribution of field RSG, with a tail extending to large radii, likely due to the presence of CSM. △ Less

Submitted 14 April, 2024; v1 submitted 25 October, 2023; originally announced October 2023.

Comments: Submitted to ApJ. Comments are welcome at ido.irani@weizmann.ac.il or idoirani@gmail.com. Accepted version

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

Shock cooling emission from explosions of red super-giants: II. An analytic model of deviations from blackbody emission

Authors: Jonathan Morag, Ido Irani, Nir Sapir, Eli Waxman

Abstract: Light emission in the first hours and days following core-collapse supernovae (SNe) is dominated by the escape of photons from the expanding shock heated envelope. In a preceding paper, Paper I, we provided a simple analytic description of the time dependent luminosity, $L$, and color temperature, $T_{\rm col}$, valid up to H recombination ($T\approx0.7$ eV), for explosions of red supergiants with… ▽ More Light emission in the first hours and days following core-collapse supernovae (SNe) is dominated by the escape of photons from the expanding shock heated envelope. In a preceding paper, Paper I, we provided a simple analytic description of the time dependent luminosity, $L$, and color temperature, $T_{\rm col}$, valid up to H recombination ($T\approx0.7$ eV), for explosions of red supergiants with convective polytropic envelopes without significant circum-stellar medium (CSM). The analytic description was calibrated against "gray" (frequency-independent) photon diffusion numeric calculations. Here we present the results of a large set of 1D multi-group (frequency-dependent) calculations, for a wide range of progenitor parameters (mass, radius, core/envelope mass ratios, metalicity) and explosion energies, using opacity tables that we constructed (and made publicly available), including the contributions of bound-bound and bound-free transitions. We provide an analytic description of the small, $\simeq10\%$ deviations of the spectrum from blackbody at low frequencies, $hν< 3T_{\rm col}$, and an improved (over Paper I) description of `line dampening' for $hν> 3T_{\rm col}$. We show that the effects of deviations from initial polytropic density distribution are small, and so are the effects of `expansion opacity' and deviations from LTE ionization and excitation (within our model assumptions). A recent study of a large set of type II SN observations finds that our model accounts well for the early multi-band data of more than 50\% of observed SNe (the others are likely affected by thick CSM), enabling the inference of progenitor properties, explosion velocity, and relative extinction. △ Less

Submitted 2 February, 2024; v1 submitted 10 July, 2023; originally announced July 2023.

Comments: 20 pages, 14 figures, 2 tables

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 scientific payload of the Ultraviolet Transient Astronomy Satellite (ULTRASAT)

Authors: Sagi Ben-Ami, Yossi Shvartzvald, Eli Waxman, Udi Netzer, Yoram Yaniv, Viktor M. Algranatti, Avishay Gal-Yam, Ofer Lapid, Eran Ofek, Jeremy Topaz, Iair Arcavi, Arooj Asif, Shlomi Azaria, Eran Bahalul, Merlin F. Barschke, Benjamin Bastian-Querner, David Berge, Vlad D. Berlea, Rolf Buhler, Louise Dittmar, Anatoly Gelman, Gianluca Giavitto, Or Guttman, Juan M. Haces Crespo, Daniel Heilbrunn , et al. (23 additional authors not shown)

Abstract: The Ultraviolet Transient Astronomy Satellite (ULTRASAT) is a space-borne near UV telescope with an unprecedented large field of view (200 sq. deg.). The mission, led by the Weizmann Institute of Science and the Israel Space Agency in collaboration with DESY (Helmholtz association, Germany) and NASA (USA), is fully funded and expected to be launched to a geostationary transfer orbit in Q2/3 of 202… ▽ More The Ultraviolet Transient Astronomy Satellite (ULTRASAT) is a space-borne near UV telescope with an unprecedented large field of view (200 sq. deg.). The mission, led by the Weizmann Institute of Science and the Israel Space Agency in collaboration with DESY (Helmholtz association, Germany) and NASA (USA), is fully funded and expected to be launched to a geostationary transfer orbit in Q2/3 of 2025. With a grasp 300 times larger than GALEX, the most sensitive UV satellite to date, ULTRASAT will revolutionize our understanding of the hot transient universe, as well as of flaring galactic sources. We describe the mission payload, the optical design and the choice of materials allowing us to achieve a point spread function of ~10arcsec across the FoV, and the detector assembly. We detail the mitigation techniques implemented to suppress out-of-band flux and reduce stray light, detector properties including measured quantum efficiency of scout (prototype) detectors, and expected performance (limiting magnitude) for various objects. △ Less

Submitted 11 March, 2023; v1 submitted 30 July, 2022; originally announced August 2022.

Comments: Presented in the SPIE Astronomical Telescopes + Instrumentation 2022

Shock cooling emission from explosions of red super-giants: I. A numerically calibrated analytic model

Authors: Jonathan Morag, Nir Sapir, Eli Waxman

Abstract: Supernova light curves are dominated at early time, hours to days, by the escape of photons from the expanding shock heated envelope. We provide a simple analytic description of the time dependent luminosity, $L$, and color temperature, $T_{\rm col}$, for explosions of red supergiants (with convective polytropic envelopes), valid up to H recombination ($T\approx0.7$ eV). The analytic description i… ▽ More Supernova light curves are dominated at early time, hours to days, by the escape of photons from the expanding shock heated envelope. We provide a simple analytic description of the time dependent luminosity, $L$, and color temperature, $T_{\rm col}$, for explosions of red supergiants (with convective polytropic envelopes), valid up to H recombination ($T\approx0.7$ eV). The analytic description is based on an interpolation between earlier analytic expressions valid at different (initial planar and later spherical) stages of the expansion, calibrated against the results of numerical hydrodynamic diffusion calculations for a wide range of progenitor parameters (mass, radius, core/envelope mass and radius ratios, metalicity), and explosion energies. The numerically derived $L$ and $T_{\rm col}$ are described by the analytic expressions with 10\% and 5\% accuracy respectively. $T_{\rm col}$ is inferred from the hydrodynamic profiles using (time and space dependent) effective "gray" (frequency independent) opacity, based on opacity tables that we have constructed for this purpose (and will be made publicly available) including the contributions of bound-bound and bound-free transitions. In an accompanying paper (Paper II) we show, using a large set of multi-group photon diffusion calculations, that the spectral energy distribution is well described by a Planck spectrum with $T=T_{\rm col}$, except at UV frequencies (beyond the spectral peak at $3T_{\rm col}$), where the flux is significantly suppressed due to the presence of strong line absorption. We defer the full discussion of the multi-group results to paper II, but provide here for completeness an analytic description also of the UV suppression. Our analytic results are a useful tool for inferring progenitor properties, explosion velocity, and also relative extinction based on early multi-band shock cooling observations of supernovae. △ Less

Submitted 26 March, 2023; v1 submitted 13 July, 2022; originally announced July 2022.

Comments: 12 pages, 11 figures. Accepted version

Non-thermal emission from mildly relativistic dynamical ejecta of neutron star mergers

Authors: Gilad Sadeh, Or Guttman, Eli Waxman

Abstract: Binary neutron star mergers are expected to produce fast dynamical ejecta, with mildly relativistic velocities extending to $β=v/c>0.6$. We consider the radio to X-ray synchrotron emission produced by collisionless shocks driven by such fast ejecta into the interstellar medium. Analytic expressions are given for spherical ejecta with broken power-law mass (or energy) distributions,… ▽ More Binary neutron star mergers are expected to produce fast dynamical ejecta, with mildly relativistic velocities extending to $β=v/c>0.6$. We consider the radio to X-ray synchrotron emission produced by collisionless shocks driven by such fast ejecta into the interstellar medium. Analytic expressions are given for spherical ejecta with broken power-law mass (or energy) distributions, $M(>γβ)\propto(γβ)^{-s}$ with $s=s_{\rm KN}$ at $γβ<γ_0β_0$ and $s=s_{\rm ft}$ at $γβ>γ_0β_0$ (where $γ$ is the Lorentz factor). For parameter values characteristic of merger calculation results -- a "shallow" mass distribution, $1<s_{\rm KN}<3$, for the bulk of the ejecta (at $γβ\approx 0.2$), and a steep, $s_{\rm ft}>5$, "fast tail" mass distribution -- our model provides an accurate (to 10's of percent) description of the evolution of the flux, including at the phase of deceleration to sub-relativistic expansion. This is a significant improvement over earlier results, based on extrapolations of results valid for $γβ\gg1$ or $\ll1$ to $γβ\approx1$, which overestimate the flux by an order of magnitude for typical parameter values. It will enable a more reliable inference of ejecta parameters from future measurements of the non-thermal emission. For the merger event GW170817, the existence of a "fast tail" is expected to produce detectable radio and X-ray fluxes over a time scale of $\sim10^4$days. △ Less

Submitted 28 November, 2022; v1 submitted 12 July, 2022; originally announced July 2022.

Comments: 11 pages, 12 figures. Minor revisions following referee's report, accepted to MNRAS

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

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

Sensor characterization for the ULTRASAT space telescope

Authors: Benjamin Bastian-Querner, Nirmal Kaipachery, Daniel Küsters, Julian Schliwinski, Shay Alfassi, Arooj Asif, Merlin F. Barschke, Sagi Ben-Ami, David Berge, Adi Birman, Rolf Bühler, Nicola De Simone, Amos Fenigstein, Avishay Gal-Yam, Gianluca Giavitto, Juan M. Haces Crespo, Dmitri Ivanov, Omer Katz, Marek Kowalski, Shrinivasrao R. Kulkarni, Ofer Lapid, Tuvia Liran, Ehud Netzer, Eran O. Ofek, Sebastian Philipp , et al. (9 additional authors not shown)

Abstract: The Ultraviolet Transient Astronomical Satellite is a scientific space mission carrying an astronomical telescope. The mission is led by the Weizmann Institute of Science in Israel and the Israel Space Agency, while the camera in the focal plane is designed and built by Deutsches Elektronen Synchrotron in Germany. Two key science goals of the mission are the detection of counterparts to gravitatio… ▽ More The Ultraviolet Transient Astronomical Satellite is a scientific space mission carrying an astronomical telescope. The mission is led by the Weizmann Institute of Science in Israel and the Israel Space Agency, while the camera in the focal plane is designed and built by Deutsches Elektronen Synchrotron in Germany. Two key science goals of the mission are the detection of counterparts to gravitational wave sources and supernovae. The launch to geostationary orbit is planned for 2024. The telescope with a field-of-view of $\approx200$deg$^2$, is optimized to work in the near-ultraviolet band between $220$ and $280$nm. The focal plane array is composed of four $22.4$-megapixel, backside-illuminated CMOS sensors with a total active area of 90x90mm$^2$. Prior to sensor production, smaller test sensors have been tested to support critical design decisions for the final flight sensor. These test sensors share the design of epitaxial layer and anti-reflective coatings (ARC) with the flight sensors. Here, we present a characterization of these test sensors. Dark current and read noise are characterized as a function of the device temperature. A temperature-independent noise level is attributed to on-die infrared emission and the read-out electronics` self-heating. We utilize a high-precision photometric calibration setup to obtain the test sensors` quantum efficiency (QE) relative to PTB/NIST-calibrated transfer standards ($220$-$1100$nm), the quantum yield for $λ< 300$nm, the non-linearity of the system, and the conversion gain. The uncertainties are discussed in the context of the newest results on the setup`s performance parameters. From three ARC options, Tstd, T1 and T2, the latter optimizes out-of-band rejection and peaks in the mid of the ULTRASAT operational waveband (max. QE $\approx80\%$ at $245\mathrm{nm}$). We recommend ARC option T2 for the final ULTRASAT UV sensor. △ Less

Submitted 5 August, 2021; originally announced August 2021.

Design of the ULTRASAT UV camera

Authors: Arooj Asif, Merlin Barschke, Benjamin Bastian-Querner, David Berge, Rolf Bühler, Nicola De Simone, Gianluca Giavitto, Juan M. Haces Crespo, Nirmal Kaipachery, Marek Kowalski, Shrinivasrao R. Kulkarni, Daniel Küsters, Sebastian Philipp, Heike Prokoph, Julian Schliwinski, Mikhail Vasilev, Jason J. Watson, Steven Worm, Francesco Zappon, Shay Alfassi, Sagi Ben-Ami, Adi Birman, Kasey Boggs, Greg Bredthauer, Amos Fenigstein , et al. (12 additional authors not shown)

Abstract: The Ultraviolet Transient Astronomical Satellite (ULTRASAT) is a scientific UV space telescope that will operate in geostationary orbit. The mission, targeted to launch in 2024, is led by the Weizmann Institute of Science (WIS) in Israel and the Israel Space Agency (ISA). Deutsches Elektronen Synchrotron (DESY) in Germany is tasked with the development of the UV-sensitive camera at the heart of th… ▽ More The Ultraviolet Transient Astronomical Satellite (ULTRASAT) is a scientific UV space telescope that will operate in geostationary orbit. The mission, targeted to launch in 2024, is led by the Weizmann Institute of Science (WIS) in Israel and the Israel Space Agency (ISA). Deutsches Elektronen Synchrotron (DESY) in Germany is tasked with the development of the UV-sensitive camera at the heart of the telescope. The camera's total sensitive area of ~90mm x 90mm is built up by four back-side illuminated CMOS sensors, which image a field of view of ~200 deg2. Each sensor has 22.4 megapixels. The Schmidt design of the telescope locates the detector inside the optical path, limiting the overall size of the assembly. As a result, the readout electronics is located in a remote unit outside the telescope. The short focal length of the telescope requires an accurate positioning of the sensors within +-50 mu along the optical axis, with a flatness of +-10 mu. While the telescope will be at around 295K during operations, the sensors are required to be cooled to 200K for dark current reduction. At the same time, the ability to heat the sensors to 343K is required for decontamination. In this paper, we present the preliminary design of the UV sensitive ULTRASAT camera. △ Less

Submitted 3 August, 2021; originally announced August 2021.

Comments: "UV, X-Ray, and Gamma-Ray Space Instrumentation for Astronomy XXII", Proc. SPIE 11821, Optics and Photonics, 11821-21 (August, 2021). Contact authors Rolf Buehler (rolf.buehler@desy.de) and Merlin Barschke (merlin.barschke@desy.de)

The GALEX-PTF experiment: II. supernova progenitor radius and energetics via shock-cooling modeling

Authors: Noam Ganot, Eran O. Ofek, Avishay Gal-Yam, Maayane T. Soumagnac, Jonathan Morag, Eli Waxman, Shrinivas R. Kulkarni, Mansi M. Kasliwal, James Neill

Abstract: The radius and surface composition of an exploding massive star, as well as the explosion energy per unit mass, can be measured using early ultraviolet (UV) observations of core-collapse supernovae (CC SNe). We present the results from a simultaneous \GALEX and Palomar Transient Factory (PTF) search for early UV emission from SNe. We analyze five CC SNe for which we obtained $NUV$ measurements bef… ▽ More The radius and surface composition of an exploding massive star, as well as the explosion energy per unit mass, can be measured using early ultraviolet (UV) observations of core-collapse supernovae (CC SNe). We present the results from a simultaneous \GALEX and Palomar Transient Factory (PTF) search for early UV emission from SNe. We analyze five CC SNe for which we obtained $NUV$ measurements before the first ground-based $R$-band detection. We introduce SOPRANOS, a new maximum likelihood fitting tool for models with variable temporal validity windows, and use it to fit the \citet{SapirWaxman2017} shock cooling model to the data. We report four Type II SNe with progenitor radii in the range of $R_*\approx600-1100R_\odot$ and a shock velocity parameter in the range of $v_{s*}\approx 2700-6000 \,\rm km\,s^{-1}$ ($E/M\approx2-8\times10^{50}\,\rm erg/M_\odot$) and one type IIb SN with $R_*\approx210R_\odot$ and $v_{s*}\approx11000 \rm\, km\,s^{-1}$ ($E/M\approx1.8\times10^{51}\,\rm erg/M_\odot$). Our pilot GALEX/PTF project thus suggests that a dedicated, systematic SN survey in the $NUV$ band, such as the wide-field UV explorer \textit{ULTRASAT} mission, is a compelling method to study the properties of SN progenitors and SN energetics. △ Less

Submitted 24 November, 2020; originally announced November 2020.

Implications of CHIME repeating fast radio bursts

Authors: Wenbin Lu, Anthony L. Piro, Eli Waxman

Abstract: CHIME has now detected 18 repeating fast radio bursts (FRBs). We explore what can be learned about the energy distribution and activity level of the repeaters by fitting realistic FRB population models to the data. For a power-law energy distribution dN/dE ~ E^{-γ} for the repeating bursts, there is a critical index γ_c that controls whether the dispersion measure (DM, a proxy for source distance)… ▽ More CHIME has now detected 18 repeating fast radio bursts (FRBs). We explore what can be learned about the energy distribution and activity level of the repeaters by fitting realistic FRB population models to the data. For a power-law energy distribution dN/dE ~ E^{-γ} for the repeating bursts, there is a critical index γ_c that controls whether the dispersion measure (DM, a proxy for source distance) distribution of repeaters is bottom or top-heavy. We find γ_c = 7/4 for Poisson wait-time distribution of repeaters in Euclidean space and further demonstrate how it is affected by temporal clustering of repetitions and cosmological effects. It is especially interesting that two of the CHIME repeaters (FRB 181017 and 190417) have large DM ~ 1000 pc/cm^3. These can be understood if: (i) the energy distribution is shallow γ= 1.7^{+0.3}_{-0.1} (68% confidence) or (ii) a small fraction of sources are extremely active. In the second scenario, these two high-DM sources should be repeating more than 100 times more frequently than FRB 121102, and the energy index is constrained to be γ= 1.9^{+0.3}_{-0.2} (68% confidence). In either case, this γis consistent with the energy dependence of the non-repeating ASKAP sample, which suggests that they are drawn from the same population. Finally, our model predicts how the CHIME repeating fraction should decrease with redshift, and this can be compared with observations to infer the distribution of activity level in the whole population. △ Less

Submitted 7 August, 2020; v1 submitted 27 March, 2020; originally announced March 2020.

Comments: 11 pages, 5 figures, MNRAS accepted after moderate revision, conclusions unchanged

High-energy particles and radiation in star-forming regions

Authors: A. M. Bykov, A. Marcowith, E. Amato, M. E. Kalyashova, J. M. D. Kruijssen, E. Waxman

Abstract: Non-thermal particles and high-energy radiation can play a role in the dynamical processes in star-forming regions and provide an important piece of the multiwavelength observational picture of their structure and components. Powerful stellar winds and supernovae in compact clusters of massive stars and OB associations are known to be favourable sites of high-energy particle acceleration and sourc… ▽ More Non-thermal particles and high-energy radiation can play a role in the dynamical processes in star-forming regions and provide an important piece of the multiwavelength observational picture of their structure and components. Powerful stellar winds and supernovae in compact clusters of massive stars and OB associations are known to be favourable sites of high-energy particle acceleration and sources of non-thermal radiation and neutrinos. Namely, young massive stellar clusters are likely sources of the PeV (petaelectronvolt) regime cosmic rays (CRs). They can also be responsible for the cosmic ray composition, e.g., 22Ne/20Ne anomalous isotopic ratio in CRs. Efficient particle acceleration can be accompanied by super-adiabatic amplification of the fluctuating magnetic fields in the systems converting a part of kinetic power of the winds and supernovae into the magnetic energy through the CR-driven instabilities. The escape and CR propagation in the vicinity of the sources are affected by the non-linear CR feedback. These effects are expected to be important in starburst galaxies, which produce high-energy neutrinos and gamma-rays. We give a brief review of the theoretical models and observational data on high-energy particle acceleration and their radiation in star-forming regions with young stellar population. △ Less

Submitted 1 October, 2020; v1 submitted 25 March, 2020; originally announced March 2020.

Comments: 41 pages, Space Science Reviews v.216, id 42, topical collection Star Formation

Journal ref: Space Science Reviews, Volume 216, Issue 3, article id.42, 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

SN 2018fif: The Explosion of a Large Red Supergiant Discovered in Its Infancy by the Zwicky Transient Facility

Authors: Maayane T. Soumagnac, Noam Ganot, Ido Irani, Avishay Gal-yam, Eran O. Ofek, Eli Waxman, Jonathan Morag, Ofer Yaron, Steve Schulze, Yi Yang, Adam Rubin, S. Bradley Cenko, Jesper Sollerman, Daniel A. Perley, Christoffer Fremling, Peter Nugent, James D. Neill, Emir Karamehmetoglu, Eric C. Bellm, Rachel J. Bruch, Rick Burruss, Virginia Cunningham, Richard Dekany, V. Zach Golkhou, Mansi M. Kasliwal , et al. (10 additional authors not shown)

Abstract: High cadence transient surveys are able to capture supernovae closer to their first light than before. Applying analytical models to such early emission, we can constrain the progenitor stars properties. In this paper, we present observations of SN2018fif (ZTF18abokyfk). The supernova was discovered close to first light and monitored by the Zwicky Transient Facility (ZTF) and the Neil Gehrels Swif… ▽ More High cadence transient surveys are able to capture supernovae closer to their first light than before. Applying analytical models to such early emission, we can constrain the progenitor stars properties. In this paper, we present observations of SN2018fif (ZTF18abokyfk). The supernova was discovered close to first light and monitored by the Zwicky Transient Facility (ZTF) and the Neil Gehrels Swift Observatory. Early spectroscopic observations suggest that the progenitor of SN2018fif was surrounded by relatively small amounts of circumstellar material (CSM) compared to all previous cases. This particularity, coupled with the high cadence multiple-band coverage, makes it a good candidate to investigate using shock-cooling models. We employ the SOPRANOS code, an implementation of the model by Sapir & Waxman and its extension to early times by Morag, Sapir & Waxman. Compared with previous implementations, SOPRANOS has the advantage of including a careful account of the limited temporal validity domain of the shock-cooling model as well as allowing usage of the entirety of the early UV data. We find that the progenitor of SN2018fif was a large red supergiant, with a radius of R=744.0_{-128.0}^{+183.0} solar radii and an ejected mass of Mej=9.3_{-5.8}^{+0.4} solar masses. Our model also gives information on the explosion epoch, the progenitor inner structure, the shock velocity and the extinction. The distribution of radii is double-peaked, with lower radii corresponding to lower values of the extinction, earlier recombination times and better match to the early UV data. If these correlations persist in future objects, denser spectroscopic monitoring constraining the time of recombination, as well as accurate UV observations (e.g. with ULTRASAT), will help break the radius-extinction degeneracy and independently determine both. △ Less

Submitted 18 October, 2020; v1 submitted 25 July, 2019; originally announced July 2019.

Comments: Accepted by ApJ

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

Maser and other instabilities in a weakly magnetized relativistic plasma: Theory and the astrophysical relevance of the maser

Authors: Andrei Gruzinov, Eli Waxman

Abstract: A sufficient condition for maser instability in a weakly magnetized relativistic plasma with an isotropic particle distribution function is given. The maser growth rates and polarizations are computed starting from the exact dielectric permittivity tensor of a magnetized plasma. For very weak magnetic fields, our results confirm the approximate validity of the 'standard maser theory', which is bas… ▽ More A sufficient condition for maser instability in a weakly magnetized relativistic plasma with an isotropic particle distribution function is given. The maser growth rates and polarizations are computed starting from the exact dielectric permittivity tensor of a magnetized plasma. For very weak magnetic fields, our results confirm the approximate validity of the 'standard maser theory', which is based on the Einstein coefficients method, with one significant exception. For inclined propagation and realistic (small but finite) field, the growth rates of the two (nearly circular) polarizations differ significantly, while the standard theory predicts two (nearly circular) polarizations with similar growth rates. We show that this deviation is due to circularly polarized synchrotron emission, which is neglected in the standard theory. The maser is shown to grow slower than Langmuir waves. Nevertheless, significant generation of EM waves is seen in (highly simplified) direct numerical simulations. We study the nonlinear saturation of the maser instability and find that it offers a mechanism for the conversion of a significant fraction of the plasma energy into radio waves. We briefly discuss the conditions under which the maser instability may operate in astrophysical sources, and provide rough estimates that may be used as a guidance when studying particular astrophysical sources/phenomena. △ Less

Submitted 7 January, 2019; originally announced January 2019.

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

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

Cosmic-ray Antimatter

Authors: Kfir Blum, Ryosuke Sato, Eli Waxman

Abstract: In recent years, space-born experiments have delivered new measurements of high energy cosmic-ray (CR) $\bar p$ and $e^+$. In addition, unprecedented sensitivity to CR composite anti-nuclei anti-d and anti-He is expected to be achieved in the near future. We report on the theoretical interpretation of these measurements. While CR antimatter is a promising discovery tool for new physics or exotic a… ▽ More In recent years, space-born experiments have delivered new measurements of high energy cosmic-ray (CR) $\bar p$ and $e^+$. In addition, unprecedented sensitivity to CR composite anti-nuclei anti-d and anti-He is expected to be achieved in the near future. We report on the theoretical interpretation of these measurements. While CR antimatter is a promising discovery tool for new physics or exotic astrophysical phenomena, an irreducible background arises from secondary production by primary CR collisions with interstellar matter. Understanding this irreducible background or constraining it from first principles is an interesting challenge. We review the attempt to obtain such understanding and apply it to CR $\bar p,\, e^+,$ anti-d and anti-He. Based on state of the art Galactic cosmic ray measurements, dominated currently by the AMS-02 experiment, we show that: (i) CR $\bar p$ most likely come from CR-gas collisions; (ii) $e^+$ data is consistent with, and suggestive of the same secondary astrophysical production mechanism responsible for $\bar p$ and dominated by proton-proton collisions. In addition, based on recent accelerator analyses we show that the flux of secondary high energy anti-He may be observable with a few years exposure of AMS-02. We highlight key open questions, as well as the role played by recent and upcoming space and accelerator data in clarifying the origins of CR antimatter. △ Less

Submitted 19 September, 2017; originally announced September 2017.

Comments: 42 pgs, 18 figs. Invited paper for RMP

On the origin of fast radio bursts (FRBs)

Authors: Eli Waxman

Abstract: We derive stringent constraints on the persistent source associated with FRB 121102: Size $10^{17}$ cm $<R<10^{18}$ cm, age $<300$ yr, characteristic electron energy $\varepsilon_e\sim0.3$ GeV, total energy $\sim10^{49}$ erg. The hot radiating plasma is confined by a cold plasma of mass $M_c<0.03 (R/10^{17.5}{\rm cm})^4 M_\odot$. The source is nearly resolved, and may be resolved by 10 GHz observa… ▽ More We derive stringent constraints on the persistent source associated with FRB 121102: Size $10^{17}$ cm $<R<10^{18}$ cm, age $<300$ yr, characteristic electron energy $\varepsilon_e\sim0.3$ GeV, total energy $\sim10^{49}$ erg. The hot radiating plasma is confined by a cold plasma of mass $M_c<0.03 (R/10^{17.5}{\rm cm})^4 M_\odot$. The source is nearly resolved, and may be resolved by 10 GHz observations. The fact that $\varepsilon_e\sim m_p c^2$ suggests that the hot plasma was created by the ejection of a mildly relativistic, $M\sim10^{-5} M_\odot$ shell, which propagated into an extended ambient medium or collided with a pre-ejected shell of mass $M_c$. The inferred plasma properties are inconsistent with typical "magnetar wind nebulae" model predictions. We suggest a physical mechanism for the generation of FRBs (independent of the persistent source model): Ejection from an underlying compact object, $R_s\sim10^{6}$ cm, of highly relativistic shells, with energy $E_s=10^{41}$ erg and Lorentz factor $γ_s$~$10^3$, into a surrounding e-p plasma with density $n\sim0.1/cm^3$ (consistent with that inferred for the plasma producing the persistent emission associated with FRB 121102). Such shell ejections with energy typical for FRB events lead to plasma conditions appropriate for strong synchrotron maser emission at the GHz range, $ν_{ coh.}\sim0.5(E/10^{41}erg)^{1/4}$ GHz. In this model, a significant fraction of the deposited energy is converted to an FRB with duration $R_s/c$, accompanied by ~10 MeV photons carrying less energy than the FRB. The inferred energy and mass associated with the source are low compared to those of typical supernova ejecta. This may suggest some type of a "weak stellar explosion", where a neutron star is formed with relatively low mass and energy ejection. However, the current upper limit on R does not allow one to rule out $M_c\sim1M_\odot$. △ Less

Submitted 29 May, 2017; v1 submitted 20 March, 2017; originally announced March 2017.

Comments: Accepted to ApJ. 7 pages, 1 figure. Some explanations expanded. Discussion of beaming added. Some references updated

Journal ref: 2017 ApJ 842:34

UV/Optical emission from the expanding envelopes of type II supernovae

Authors: Nir Sapir, Eli Waxman

Abstract: The early part of a supernova (SN) light-curve is dominated by radiation escaping from the expanding shock-heated progenitor envelope. For polytropic Hydrogen envelopes, the properties of the emitted radiation are described by simple analytic expressions and are nearly independent of the polytropic index, $n$. This analytic description holds at early time, $t<$~few days, during which radiation esc… ▽ More The early part of a supernova (SN) light-curve is dominated by radiation escaping from the expanding shock-heated progenitor envelope. For polytropic Hydrogen envelopes, the properties of the emitted radiation are described by simple analytic expressions and are nearly independent of the polytropic index, $n$. This analytic description holds at early time, $t<$~few days, during which radiation escapes from shells initially lying near the stellar surface. We use numerical solutions to address two issues. First, we show that the analytic description holds at early time also for non-polytropic density profiles. Second, we extend the solutions to later times, when the emission emerges from deep within the envelope and depends on the progenitor's density profile. Examining the late time behavior of polytropic envelopes with a wide range of core to envelope mass and radius ratios, $0.1\le M_{\rm c}/M_{\rm env}\le10$ and $10^{-3}\le R_{\rm c}/R\le10^{-1}$, we find that the effective temperature is well described by the analytic solution also at late time, while the luminosity $L$ is suppressed by a factor, which may be approximated to better than $20[30]\%$ accuracy up to $t=t_{\rm tr}/a$ by $A\exp[-(at/t_{\rm tr})^α]$ with $t_{\rm tr} = 15 (M_{\rm env}/M_\odot)^{3/4}(E/10^{51}{\rm erg})^{-1/4}~\rm d$, $A=0.9[0.8]$, $a=1.7[4.6]$ and $α=0.8[0.7]$ for $n=3/2[3]$. This description holds as long as the opacity is approximately that of a fully ionized gas, i.e. for $T>0.7~\rm eV$, $t<14(R/10^{13.5}{\rm cm})^{0.55}~\rm d$. The suppression of $L$ at $t_{\rm tr}/a$ obtained for standard polytropic envelopes may account for the first optical peak of double-peaked SN light curves, with first peak at a few days for $M_{\rm env}<1M_\odot$. △ Less

Submitted 5 March, 2017; v1 submitted 13 July, 2016; originally announced July 2016.

Comments: 11 pages, 18 figures. Accepted to ApJ. Expanded comparison with previous work. Expanded discussion

Constraining High-Energy Cosmic Neutrino Sources: Implications and Prospects

Authors: Kohta Murase, Eli Waxman

Abstract: We consider limits on the local ($z=0$) density ($n_0$) of extragalactic neutrino sources set by the nondetection of steady high-energy neutrino sources producing $\gtrsim50$ TeV muon multiplets in the present IceCube data, taking into account the redshift evolution, luminosity function and neutrino spectrum of the sources. We show that the lower limit depends moderately on source spectra and stro… ▽ More We consider limits on the local ($z=0$) density ($n_0$) of extragalactic neutrino sources set by the nondetection of steady high-energy neutrino sources producing $\gtrsim50$ TeV muon multiplets in the present IceCube data, taking into account the redshift evolution, luminosity function and neutrino spectrum of the sources. We show that the lower limit depends moderately on source spectra and strongly on redshift evolution. We find $n_0\gtrsim{10}^{-8}-{10}^{-7}~{\rm Mpc}^{-3}$ for standard candle sources evolving rapidly, $n_s\propto{(1+z)}^3$, and $n_0\gtrsim{10}^{-6}-{10}^{-5}~{\rm Mpc}^{-3}$ for nonevolving sources. The corresponding upper limits on their neutrino luminosity are $L_{ν_μ}^{\rm eff}\lesssim10^{42}-10^{43}~{\rm erg}~{\rm s}^{-1}$ and $L_{ν_μ}^{\rm eff}\lesssim10^{41}-10^{42}~{\rm erg}~{\rm s}^{-1}$, respectively. Applying these results to a wide range of classes of potential sources, we show that powerful blazar jets associated with active galactic nuclei are unlikely to be the dominant sources. For almost all other steady candidate source classes (including starbursts, radio galaxies, and galaxy clusters and groups), an order of magnitude increase in the detector sensitivity at $\sim0.1-1$ PeV will enable a detection (as point sources) of the few brightest objects. Such an increase, which may be provided by next-generation detectors like IceCube-Gen2 and an upgraded KM3NET, can improve the limit on $n_0$ by more than two orders of magnitude. Future gamma-ray observations (by Fermi, HAWC and CTA) will play a key role in confirming the association of the neutrinos with their sources. △ Less

Submitted 25 October, 2016; v1 submitted 6 July, 2016; originally announced July 2016.

Comments: 14 pages, 7 figures, 1 table, accepted for publication in PRD, minor changes (references and discussions added)

Journal ref: Phys.Rev.D 94 (2016) 103006

Shock breakout theory

Authors: Eli Waxman, Boaz Katz

Abstract: The earliest supernova (SN) emission is produced when the optical depth of the plasma lying ahead of the shock, which ejects the envelope, drops below c/v, where v is the shock velocity. This "breakout" may occur when the shock reaches the edge of the star, producing a bright X-ray/UV flash on time scales of seconds to a fraction of an hour, followed by UV/optical "cooling" emission from the expan… ▽ More The earliest supernova (SN) emission is produced when the optical depth of the plasma lying ahead of the shock, which ejects the envelope, drops below c/v, where v is the shock velocity. This "breakout" may occur when the shock reaches the edge of the star, producing a bright X-ray/UV flash on time scales of seconds to a fraction of an hour, followed by UV/optical "cooling" emission from the expanding cooling envelope on a day time-scale. If the optical depth of circumstellar material (CSM) ejected from the progenitor star prior to the explosion is larger than c/v, the breakout will take place at larger radii, within the CSM, extending its duration to days time scale. The properties of the early, breakout and cooling, emission carry unique signatures of the structure of the progenitor star (e.g. its radius and surface composition) and of its mass-loss history. The recent progress of wide-field transient surveys enable SN detections on a day time scale, and are being used to set unique constraints on the progenitors of SNe of all types. This chapter includes a pedagogical description of SN breakout theory, and a concise overview of what we have learned from observations so far, and of advances in observational capabilities that are required in order to make further significant progress. △ Less

Submitted 24 April, 2017; v1 submitted 5 July, 2016; originally announced July 2016.

Comments: review chapter, to appear in Handbook of Supernovae, edited by Athem W. Alsabti and Paul Murdin (Springer)

The origin of IceCube's neutrinos: Cosmic ray accelerators embedded in star forming calorimeters

Authors: E. Waxman

Abstract: The IceCube collaboration reports a detection of extra-terrestrial neutrinos. The isotropy and flavor content of the signal, and the coincidence, within current uncertainties, of the 50 TeV to 2 PeV flux and the spectrum with the Waxman-Bahcall bound, suggest a cosmological origin of the neutrinos, related to the sources of ultra-high energy, $>10^{10}$ GeV, cosmic-rays (UHECR). The most natural e… ▽ More The IceCube collaboration reports a detection of extra-terrestrial neutrinos. The isotropy and flavor content of the signal, and the coincidence, within current uncertainties, of the 50 TeV to 2 PeV flux and the spectrum with the Waxman-Bahcall bound, suggest a cosmological origin of the neutrinos, related to the sources of ultra-high energy, $>10^{10}$ GeV, cosmic-rays (UHECR). The most natural explanation of the UHECR and neutrino signals is that both are produced by the same population of cosmological sources, producing CRs (likely protons) at a similar rate, $E^2d\dot{n}/dE\propto E^{0}$, over the [$1$ PeV,$10^{11}$ GeV] energy range, and residing in "calorimetric" environments, like galaxies with high star formation rate, in which $E/Z<100$ PeV CRs lose much of their energy to pion production. A tenfold increase in the effective mass of the detector at $\gtrsim100$ TeV is required in order to significantly improve the accuracy of current measurements, to enable the detection of a few bright nearby starburst "calorimeters", and to open the possibility of identifying the CR sources embedded within the calorimeters, by associating neutrinos with photons accompanying transient events responsible for their generation. Source identification and a large neutrino sample may enable one to use astrophysical neutrinos to constrain new physics models. △ Less

Submitted 3 November, 2015; originally announced November 2015.

Comments: To appear in "Neutrino Astronomy- Current status, future prospects", Eds. T. Gaisser & A. Karle (World Scientific)

The energy budget of GRBs based on updated prompt \& afterglow observations

Authors: Nahliel Wygoda, Dafne Guetta, Marc-Adrien Mandich, Eli Waxman

Abstract: We compare the isotropic equivalent 15-2000 keV gamma-ray energy, E_gamma, emitted by a sample of 91 swift Gamma-Ray Bursts (GRBs) with known redshifts, with the isotropic equivalent fireball energy, E_fb, as estimated within the fireball model framework from X-ray afterglow observations of these bursts. The uncertainty in E_gamma, which spans the range of ~10^51 erg to ~10^53.5 erg, is approximat… ▽ More We compare the isotropic equivalent 15-2000 keV gamma-ray energy, E_gamma, emitted by a sample of 91 swift Gamma-Ray Bursts (GRBs) with known redshifts, with the isotropic equivalent fireball energy, E_fb, as estimated within the fireball model framework from X-ray afterglow observations of these bursts. The uncertainty in E_gamma, which spans the range of ~10^51 erg to ~10^53.5 erg, is approximately 25% on average, due mainly to the extrapolation from the BAT detector band to the 15-2000 keV band. The uncertainty in E_fb is approximately a factor of 2, due mainly to the X-ray measurements' scatter. We find E_gamma and E_fb to be tightly correlated. The average(std) of η^11hr_gamma is approximately log_10(E_gamma/(3ε _eE^11hr_fb)) are -0.34(0.60), and the upper limit on the intrinsic spread of η_gamma is approximately 0.5 (ε_e is the fraction of shocked plasma energy carried by electrons and E^x hr_fb is inferred from the X-ray flux at x hours). We also find that E_fb inferred from X-ray observations at 3 and 11 hours are similar, with an average(std) of log_10(E^3hr_fb/E^11hr_fb) of 0.04(0.28). The small variance of η_gamma implies that burst-to-burst variations in ε_e and in the efficiency of fireball energy conversion to gamma-rays are small, and suggests that both are of order unity. The small variance of η_gamma and the similarity of E^3hr_fb and E^11hr_fb further imply that ε_e does not vary significantly with shock Lorentz factor, and that for most bursts the modification of fireball energy during the afterglow phase, by processes such as radiative losses or extended duration energy injection, are not significant. Finally, our results imply that if fireballs are indeed jets, then the jet opening angle satisfies θ>0.1 for most cases. [abridged] △ Less

Submitted 23 November, 2015; v1 submitted 4 April, 2015; originally announced April 2015.

Comments: An expanded discussion of the absence of energy injection and 'flat X-ray phases' in the majority of bursts

The detection rate of early UV emission from supernovae: A dedicated GALEX/PTF survey and calibrated theoretical estimates

Authors: Noam Ganot, Avishay Gal-Yam, Eran O. Ofek, Ilan Sagiv, Eli Waxman, Ofer Lapid, Shrinivas R. Kulkarni, Sagi Ben-Ami, Mansi M. Kasliwal, Doron Chelouche, Stephen Rafter, Ehud Behar, Ari Laor, Dovi Poznanski, Udi Nakar, Dan Maoz, Benny Trakhtenbrot, James D. Neill, Thomas A. Barlow, Christofer D. Martin, Suvi Gezari, Iair Arcavi, Joshua s. Bloom, Peter E. Nugent, Mark Sullivan

Abstract: The radius and surface composition of an exploding massive star,as well as the explosion energy per unit mass, can be measured using early UV observations of core collapse supernovae (SNe). We present the first results from a simultaneous GALEX/PTF search for early UV emission from SNe. Six Type II SNe and one Type II superluminous SN (SLSN-II) are clearly detected in the GALEX NUV data. We compar… ▽ More The radius and surface composition of an exploding massive star,as well as the explosion energy per unit mass, can be measured using early UV observations of core collapse supernovae (SNe). We present the first results from a simultaneous GALEX/PTF search for early UV emission from SNe. Six Type II SNe and one Type II superluminous SN (SLSN-II) are clearly detected in the GALEX NUV data. We compare our detection rate with theoretical estimates based on early, shock-cooling UV light curves calculated from models that fit existing Swift and GALEX observations well, combined with volumetric SN rates. We find that our observations are in good agreement with calculated rates assuming that red supergiants (RSGs) explode with fiducial radii of 500 solar, explosion energies of 10^51 erg, and ejecta masses of 10 solar masses. Exploding blue supergiants and Wolf-Rayet stars are poorly constrained. We describe how such observations can be used to derive the progenitor radius, surface composition and explosion energy per unit mass of such SN events, and we demonstrate why UV observations are critical for such measurements. We use the fiducial RSG parameters to estimate the detection rate of SNe during the shock-cooling phase (<1d after explosion) for several ground-based surveys (PTF, ZTF, and LSST). We show that the proposed wide-field UV explorer ULTRASAT mission, is expected to find >100 SNe per year (~0.5 SN per deg^2), independent of host galaxy extinction, down to an NUV detection limit of 21.5 mag AB. Our pilot GALEX/PTF project thus convincingly demonstrates that a dedicated, systematic SN survey at the NUV band is a compelling method to study how massive stars end their life. △ Less

Submitted 25 February, 2015; v1 submitted 12 December, 2014; originally announced December 2014.

Comments: See additional information including animations on http://www.weizmann.ac.il/astrophysics/ultrasat

Journal ref: 2016 ApJ, 820, 57

IceCube's Neutrinos: The beginning of extra-Galactic neutrino astrophysics?

Authors: E. Waxman

Abstract: The flux, spectrum and angular distribution of the excess neutrino signal detected by IceCube between 50TeV and 2PeV are inconsistent with those expected for Galactic sources. The coincidence of the excess, $E_ν^2Φ_ν=3.6\pm1.2\times10^{-8}(GeV/ cm^2 sr s)$, with the Waxman-Bahcall (WB) bound, $E_ν^2Φ_{WB}=3.4\times10^{-8}(GeV/cm^2 sr s)$, is probably a clue to the origin of IceCube's neutrinos. Th… ▽ More The flux, spectrum and angular distribution of the excess neutrino signal detected by IceCube between 50TeV and 2PeV are inconsistent with those expected for Galactic sources. The coincidence of the excess, $E_ν^2Φ_ν=3.6\pm1.2\times10^{-8}(GeV/ cm^2 sr s)$, with the Waxman-Bahcall (WB) bound, $E_ν^2Φ_{WB}=3.4\times10^{-8}(GeV/cm^2 sr s)$, is probably a clue to the origin of IceCube's neutrinos. The most natural explanation of this coincidence is that both the neutrino excess and the ultra-high energy, $>10^{19}$ eV, cosmic-ray (UHECR) flux are produced by the same population of cosmologically distributed sources, producing CRs, likely protons, at a similar rate, $E^2 dQ/dE=0.5\times10^{44}(erg/Mpc^3yr)$ (at z=0), across a wide range of energies, from $10^{15}$ eV to $>10^{20}$ eV, and residing in environments (such as starburst galaxies) in which CRs of rigidity $E/Z< 10^{17}$ eV lose much of their energy to pion production. Identification of the neutrino sources will allow one to identify the UHECR accelerators, to resolve open questions related to the accelerator models, and to study neutrino properties (related e.g. to flavor oscillations and coupling to gravity) with an accuracy many orders of magnitude better than is currently possible. The most promising method for identifying the sources is by association of a neutrino with an electromagnetic signal accompanying a transient event responsible for its generation. The neutrino flux that is produced within the sources, and that may thus be directly associated with transient events, may be significantly lower than the total observed neutrino flux, which may be dominated by neutrino production at the environment in which the sources reside. △ Less

Submitted 2 December, 2013; originally announced December 2013.

Comments: 8 pages, 2 figures, to appear in the Proc. of the 9th Rencontres du Vietnam : Windows on the universe

The energy production rate density of cosmic rays in the local universe is $\sim10^{44-45}\rm erg~Mpc^{-3}~yr^{-1}$ at all particle energies

Authors: Boaz Katz, Eli Waxman, Todd Thompson, Abraham Loeb

Abstract: The energy output (per logarithmic interval of particle energies) of Cosmic Rays (CRs) with energies $10{\rm GeV}\lesssim\varepsilon_p\lesssim100{\rm GeV}$ is $\sim 10^{47}\rm erg$ per solar mass of star$-$formation, based on the CR production rate in the Milky Way and in starburst galaxies, implying a generation rate of $\varepsilon_p^2Q\sim 10^{45}\rm erg~Mpc^{-3}~yr^{-1}$ in the local universe.… ▽ More The energy output (per logarithmic interval of particle energies) of Cosmic Rays (CRs) with energies $10{\rm GeV}\lesssim\varepsilon_p\lesssim100{\rm GeV}$ is $\sim 10^{47}\rm erg$ per solar mass of star$-$formation, based on the CR production rate in the Milky Way and in starburst galaxies, implying a generation rate of $\varepsilon_p^2Q\sim 10^{45}\rm erg~Mpc^{-3}~yr^{-1}$ in the local universe. It is only $\sim 10$ times larger than the output, $\varepsilon_p^2 Q=0.5\pm0.2\times 10^{44}\rm erg~Mpc^{-3}~yr^{-1}$, of Ultra High Energy CRs (UHECRs) at energies $10^{10.5}{\rm GeV}<\varepsilon_p<10^{12}\rm GeV$ (obtained assuming they are mostly protons), which in turn is comparable to the lower limit of $\varepsilon_p^2 Q\ge 0.5\times 10^{44}\rm erg~Mpc^{-3}~yr^{-1}$ of high energy CRs with $10^6{\rm GeV}\lesssim\varepsilon_p\lesssim 10^{8}\rm GeV$ implied by the saturation of the Waxman-Bahcall bound by the neutrino excess recently discovered by IceCube. These similarities are consistent with a flat production spectrum, $\varepsilon_p^2 Q\sim \text{const}$ for CRs at all observed energies. If a flat production spectrum is generated by our galaxy, the observed CR flux in the range $10^{6.5}-10^{9.5}\rm GeV$, above the "knee", is suppressed compared to lower energies due to propagation effects rather than acceleration upper limits. As suggested by Parizot and Aublin, the most exciting possibility is that cosmic rays at all energies are emitted from a single type of (unknown) sources, which can not be supernova remnants. △ Less

Submitted 1 November, 2013; originally announced November 2013.

Comments: 5 pages, 1 figure

Constraints on the source of ultra-high energy cosmic rays using anisotropy vs chemical composition

Authors: Ruo-Yu Liu, Andrew Taylor, Martin Lemoine, Xiang-Yu Wang, Eli Waxman

Abstract: The joint analysis of anisotropy signals and chemical composition of ultra-high energy cosmic rays offers strong potential for shedding light on the sources of these particles. Following up on an earlier idea, this paper studies the anisotropies produced by protons of energy >E/Z, assuming that anisotropies at energy >E have been produced by nuclei of charge Z, which share the same magnetic rigidi… ▽ More The joint analysis of anisotropy signals and chemical composition of ultra-high energy cosmic rays offers strong potential for shedding light on the sources of these particles. Following up on an earlier idea, this paper studies the anisotropies produced by protons of energy >E/Z, assuming that anisotropies at energy >E have been produced by nuclei of charge Z, which share the same magnetic rigidity. We calculate the number of secondary protons produced through photodisintegration of the primary heavy nuclei. Making the extreme assumption that the source does not inject any proton, we find that the source(s) responsible for anisotropies such as reported by the Pierre Auger Observatory should lie closer than ~20-30, 80-100 and 180-200 Mpc if the anisotropy signal is mainly composed of oxygen, silicon and iron nuclei respectively. A violation of this constraint would otherwise result in the secondary protons forming a more significant anisotropy signal at lower energies. Even if the source were located closer than this distance, it would require an extraordinary metallicity >120, 1600, 1100 times solar metallicity in the acceleration zone of the source, for oxygen, silicon and iron respectively, to ensure that the concomitantly injected protons not to produce a more significant low energy anisotropy. This offers interesting prospects for constraining the nature and the source of ultra-high energy cosmic rays with the increase in statistics expected from next generation detectors. △ Less

Submitted 18 October, 2013; v1 submitted 26 August, 2013; originally announced August 2013.

Comments: 16 pages, 3 figures; published in ApJ; comments on Auger's paper modified

Journal ref: The Astrophysical Journal, Volume 776, Issue 2, article id. 88, 10 pp. (2013)

AMS02 results support the secondary origin of cosmic ray positrons

Authors: Kfir Blum, Boaz Katz, Eli Waxman

Abstract: We show that the recent AMS02 positron fraction measurement is consistent with a secondary origin for positrons, and does not require additional primary sources such as pulsars or dark matter. The measured positron fraction at high energy saturates the previously predicted upper bound for secondary production (Katz et al 2009), obtained by neglecting radiative losses. This coincidence, which will… ▽ More We show that the recent AMS02 positron fraction measurement is consistent with a secondary origin for positrons, and does not require additional primary sources such as pulsars or dark matter. The measured positron fraction at high energy saturates the previously predicted upper bound for secondary production (Katz et al 2009), obtained by neglecting radiative losses. This coincidence, which will be further tested by upcoming AMS02 data at higher energy, is a compelling indication for a secondary source. Within the secondary model the AMS02 data imply a cosmic ray propagation time in the Galaxy of < Myr and an average traversed interstellar matter density of order 1/cc, comparable to the density of the Milky Way gaseous disk, at a rigidity of 300 GV. △ Less

Submitted 9 November, 2013; v1 submitted 6 May, 2013; originally announced May 2013.

Comments: 6 pages, 7 figures. v2: include updated results from AMS02 (B/C, e+ flux), match journal version

Journal ref: Phys. Rev. Lett. 111, 211101 (2013)

Non-relativistic radiation mediated shock breakouts: III. Spectral properties of SN shock breakout

Authors: Nir Sapir, Boaz Katz, Eli Waxman

Abstract: The spectrum of radiation emitted following shock breakout from a star's surface with a power-law density profile $ρ\propto x^n$ is investigated. Assuming planar geometry, local Compton equilibrium and bremsstrahlung emission as the dominant photon production mechanism, numerical solutions are obtained for the photon number density and temperature profiles as a function of time, for hydrogen-heliu… ▽ More The spectrum of radiation emitted following shock breakout from a star's surface with a power-law density profile $ρ\propto x^n$ is investigated. Assuming planar geometry, local Compton equilibrium and bremsstrahlung emission as the dominant photon production mechanism, numerical solutions are obtained for the photon number density and temperature profiles as a function of time, for hydrogen-helium envelopes. The temperature solutions are determined by the breakout shock velocity $v_0$ and the pre-shock breakout density $ρ_0$, and depend weakly on the value of n. Fitting formulas for the peak surface temperature at breakout as a function of $v_0$ and $ρ_0$ are provided, with $T_{peak}\approx 9.44\exp{[12.63(v_0/c)^{1/2}]}$ eV, and the time dependence of the surface temperature is tabulated. The time integrated emitted spectrum is a robust prediction of the model, determined by $\mathcal{T}_{\rm peak}$ and $v_0$ alone and insensitive to details of light travel time or slight deviations from spherical symmetry. Adopting commonly assumed progenitor parameters, breakout luminosities of ~10^45 erg/s and ~10^44 erg/s in the 0.3-10 keV band are expected for BSG and RSG/He-WR progenitors respectively (T_{peak} is well below the band for RSGs, unless their radius is ~10^13 cm). >30 detections of SN1987A-like (BSG) breakouts are expected over the lifetime of ROSAT and XMM-Newton. An absence of such detections would imply that either the typical parameters assumed for BSG progenitors are grossly incorrect or that their envelopes are not hydrostatic. The observed spectrum and duration of XRF 080109/SN2008D are in tension with a non-relativistic breakout from a stellar surface interpretation. △ Less

Submitted 20 August, 2013; v1 submitted 23 April, 2013; originally announced April 2013.

Comments: 14 pages, 8 figures, 3 tables

Journal ref: ApJ 774 (2013) 79

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

Science with a wide-field UV transient explorer

Authors: I. Sagiv, A. Gal-Yam, E. O. Ofek, E. Waxman, O. Aharonson, E. Nakar, D. Maoz, B. Trakhtenbrot, S. R. Kulkarni, E. S. Phinney, J. Topaz, C. Beichman, J. Murthy, S. P. Worden

Abstract: The time-variable electromagnetic sky has been well-explored at a wide range of wavelengths. Numerous high-energy space missions take advantage of the dark Gamma-ray and X-ray sky and utilize very wide field detectors to provide almost continuous monitoring of the entire celestial sphere. In visible light, new wide-field ground-based surveys cover wide patches of sky with ever decreasing cadence,… ▽ More The time-variable electromagnetic sky has been well-explored at a wide range of wavelengths. Numerous high-energy space missions take advantage of the dark Gamma-ray and X-ray sky and utilize very wide field detectors to provide almost continuous monitoring of the entire celestial sphere. In visible light, new wide-field ground-based surveys cover wide patches of sky with ever decreasing cadence, progressing from monthly-weekly time scale surveys to sub-night sampling. In the radio, new powerful instrumentation offers unprecedented sensitivity over wide fields of view, with pathfinder experiments for even more ambitious programs underway. In contrast, the ultra-violet (UV) variable sky is relatively poorly explored, even though it offers exciting scientific prospects. Here, we review the potential scientific impact of a wide-field UV survey on the study of explosive and other transient events, as well as known classes of variable objects, such as active galactic nuclei and variable stars. We quantify our predictions using a fiducial set of observational parameters which are similar to those envisaged for the proposed ULTRASAT mission. We show that such a mission would be able to revolutionize our knowledge about massive star explosions by measuring the early UV emission from hundreds of events, revealing key physical parameters of the exploding progenitor stars. Such a mission would also detect the UV emission from many tens of tidal-disruption events of stars by super massive black holes at galactic nuclei and enable a measurement of the rate of such events. The overlap of such a wide-field UV mission with existing and planned gravitational-wave and high-energy neutrino telescopes makes it especially timely. △ Less

Submitted 25 March, 2013; originally announced March 2013.

Preliminary evidence for a virial shock around the Coma galaxy cluster

Authors: Uri Keshet, Doron Kushnir, Abraham Loeb, Eli Waxman

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

Submitted 21 September, 2017; v1 submitted 4 October, 2012; originally announced October 2012.

Comments: Replaced with published version

Journal ref: The Astrophysical Journal, Volume 845, Issue 1, article id. 24, 14 pp. (2017)

A Search for Correlation of Ultra-High Energy Cosmic Rays with IRAS-PSCz and 2MASS-6dF Galaxies

Authors: Foteini Oikonomou, Amy Connolly, Filipe B. Abdalla, Ofer Lahav, Shaun A. Thomas, David Waters, Eli Waxman

Abstract: We study the arrival directions of 69 ultra-high energy cosmic rays (UHECRs) observed at the Pierre Auger Observatory (PAO) with energies exceeding 55 EeV. We investigate whether the UHECRs exhibit the anisotropy signal expected if the primary particles are protons that originate in galaxies in the local universe, or in sources correlated with these galaxies. We cross-correlate the UHECR arrival d… ▽ More We study the arrival directions of 69 ultra-high energy cosmic rays (UHECRs) observed at the Pierre Auger Observatory (PAO) with energies exceeding 55 EeV. We investigate whether the UHECRs exhibit the anisotropy signal expected if the primary particles are protons that originate in galaxies in the local universe, or in sources correlated with these galaxies. We cross-correlate the UHECR arrival directions with the positions of IRAS-PSCz and 2MASS-6dF galaxies taking into account particle energy losses during propagation. This is the first time that the 6dF survey is used in a search for the sources of UHECRs and the first time that the PSCz survey is used with the full 69 PAO events. The observed cross-correlation signal is larger for the PAO UHECRs than for 94% (98%) of realisations from an isotropic distribution when cross-correlated with the PSCz (6dF). On the other hand the observed cross-correlation signal is lower than that expected from 85% of realisations, had the UHECRs originated in galaxies in either survey. The observed cross-correlation signal does exceed that expected by 50% of the realisations if the UHECRs are randomly deflected by intervening magnetic fields by 5 degrees or more. We propose a new method of analysing the expected anisotropy signal, by dividing the predicted UHECR source distribution into equal predicted flux radial shells, which can help localise and constrain the properties of UHECR sources. We find that the 69 PAO events are consistent with isotropy in the nearest of three shells we define, whereas there is weak evidence for correlation with the predicted source distribution in the two more distant shells in which the galaxy distribution is less anisotropic. △ Less

Submitted 18 May, 2013; v1 submitted 17 July, 2012; originally announced July 2012.

Comments: 23 pages, version published in JCAP

Journal ref: JCAP05(2013)015

Multimessenger astronomy with gravitational waves and high-energy neutrinos

Authors: S. Ando, B. Baret, B. Bouhou, E. Chassande-Mottin, A. Kouchner, L. Moscoso, V. Van Elewyck, I. Bartos, S. Márka, Z. Márka, A. Corsi, I. Di Palma, M. A. Papa, A. Dietz, C. Donzaud, D. Eichler, C. Finley, D. Guetta, F. Halzen, G. Jones, P. J. Sutton, S. Kandhasamy, V. Mandic, E. Thrane, K. Kotake , et al. (4 additional authors not shown)

Abstract: Many of the astrophysical sources and violent phenomena observed in our Universe are potential emitters of gravitational waves (GW) and high-energy neutrinos (HEN). Both GWs and HENs may escape very dense media and travel unaffected over cosmological distances, carrying information from the innermost regions of the astrophysical engines. Such messengers could also reveal new, hidden sources that h… ▽ More Many of the astrophysical sources and violent phenomena observed in our Universe are potential emitters of gravitational waves (GW) and high-energy neutrinos (HEN). Both GWs and HENs may escape very dense media and travel unaffected over cosmological distances, carrying information from the innermost regions of the astrophysical engines. Such messengers could also reveal new, hidden sources that have not been observed by conventional photon-based astronomy. Coincident observation of GWs and HENs may thus play a critical role in multimessenger astronomy. This is particularly true at the present time owing to the advent of a new generation of dedicated detectors: IceCube, ANTARES, VIRGO and LIGO. Given the complexity of the instruments, a successful joint analysis of this data set will be possible only if the expertise and knowledge of the data is shared between the two communities. This review aims at providing an overview of both theoretical and experimental state-of-the-art and perspectives for such a GW+HEN multimessenger astronomy. △ Less

Submitted 25 July, 2013; v1 submitted 23 March, 2012; originally announced March 2012.

Comments: Revised version, to appear in Reviews of Modern Physics (Colloquium Series). Section II on Potential GW+HEN emitters thoroughly updated. Section IV expanded with recent experimental results from joint GW+HEN searches with ANTARES, IceCube, LIGO and VIRGO, and with an extended discussion of the reach and implications of future searches

Journal ref: Reviews of Modern Physics 85 (2013) 1401-1420

Early emission from type Ia supernovae

Authors: Itay Rabinak, Eli Livne, Eli Waxman

Abstract: A unique feature of deflagration-to-detonation (DDT) white dwarf explosion models of SNe of type Ia is the presence of a strong shock wave propagating through the outer envelope. We consider the early emission expected in such models, which is produced by the expanding shock-heated outer part of the ejecta and precedes the emission driven by radioactive decay. We expand on earlier analyses by cons… ▽ More A unique feature of deflagration-to-detonation (DDT) white dwarf explosion models of SNe of type Ia is the presence of a strong shock wave propagating through the outer envelope. We consider the early emission expected in such models, which is produced by the expanding shock-heated outer part of the ejecta and precedes the emission driven by radioactive decay. We expand on earlier analyses by considering the modification of the pre-detonation density profile by the weak-shocks generated during the deflagration phase, the time evolution of the opacity, and the deviation of the post-shock equation of state from that obtained for radiation pressure domination. A simple analytic model is presented and shown to provide an acceptable approximation to the results of 1D numerical DDT simulations. Our analysis predicts a thousand second long UV/optical flash with a luminosity of ~1 to 3*1e39 erg/s. Lower luminosity corresponds to faster (turbulent) deflagration velocity. The predicted luminosity of the UV flash is an order of magnitude lower than that of earlier estimates, and is expected to be strongly suppressed at times longer than an hour due to the deviation from pure radiation domination. △ Less

Submitted 29 August, 2011; originally announced August 2011.

Comments: 10 pages, 4 figures

Imploding ignition waves: I. one dimensional analysis

Authors: Doron Kushnir, Eli Livne, Eli Waxman

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

Submitted 2 June, 2012; v1 submitted 23 August, 2011; originally announced August 2011.

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

PTF11eon/SN2011dh: Discovery of a Type IIb Supernova From a Compact Progenitor in the Nearby Galaxy M51

Authors: Iair Arcavi, Avishay Gal-Yam, Ofer Yaron, Assaf Sternberg, Itay Rabinak, Eli Waxman, Mansi M. Kasliwal, Robert M. Quimby, Eran O. Ofek, Assaf Horesh, Shrinivas R. Kulkarni, Alexei V. Filippenko, Jeffrey M. Silverman, S. Bradley Cenko, Weidong Li, Joshua S. Bloom, Mark Sullivan, Derek B. Fox, Peter E. Nugent, Dovi Poznanski, Evgeny Gorbikov, Amedee Riou, Stephane Lamotte-Bailey, Thomas Griga, Judith G. Cohen , et al. (18 additional authors not shown)

Abstract: On May 31, 2011 UT a supernova (SN) exploded in the nearby galaxy M51 (the Whirlpool Galaxy). We discovered this event using small telescopes equipped with CCD cameras, as well as by the Palomar Transient Factory (PTF) survey, and rapidly confirmed it to be a Type II supernova. Our early light curve and spectroscopy indicates that PTF11eon resulted from the explosion of a relatively compact progen… ▽ More On May 31, 2011 UT a supernova (SN) exploded in the nearby galaxy M51 (the Whirlpool Galaxy). We discovered this event using small telescopes equipped with CCD cameras, as well as by the Palomar Transient Factory (PTF) survey, and rapidly confirmed it to be a Type II supernova. Our early light curve and spectroscopy indicates that PTF11eon resulted from the explosion of a relatively compact progenitor star as evidenced by the rapid shock-breakout cooling seen in the light curve, the relatively low temperature in early-time spectra and the prompt appearance of low-ionization spectral features. The spectra of PTF11eon are dominated by H lines out to day 10 after explosion, but initial signs of He appear to be present. Assuming that He lines continue to develop in the near future, this SN is likely a member of the cIIb (compact IIb; Chevalier and Soderberg 2010) class, with progenitor radius larger than that of SN 2008ax and smaller than the eIIb (extended IIb) SN 1993J progenitor. Our data imply that the object identified in pre-explosion Hubble Space Telescope images at the SN location is possibly a companion to the progenitor or a blended source, and not the progenitor star itself, as its radius (~10^13 cm) would be highly inconsistent with constraints from our post-explosion photometric and spectroscopic data. △ Less

Submitted 17 June, 2011; originally announced June 2011.