Source author record

Eli Waxman

Eli Waxman appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

49works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

49 published item(s)

preprint2020arXiv

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

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.

preprint2020arXiv

Implications of CHIME repeating fast radio bursts

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.

preprint2016arXiv

Constraining High-Energy Cosmic Neutrino Sources: Implications and Prospects

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.

preprint2015arXiv

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

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.

preprint2015arXiv

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

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]

preprint2013arXiv

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

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.

preprint2013arXiv

AMS02 results support the secondary origin of cosmic ray positrons

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.

preprint2013arXiv

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

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.

preprint2013arXiv

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

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.

preprint2013arXiv

The early UV/Optical emission from core-collapse supernovae

We derive a simple approximate model describing the early, hours to days, UV/optical supernova emission, which is produced by the expansion of the outer <~0.01 solar mass part of the shock-heated envelope, and precedes the optical emission driven by radioactive decay. Our model includes an approximate description of the time dependence of the opacity (due mainly to recombination), and of the deviation of the emitted spectrum from a black body spectrum. We show that the characteristics of the early UV/O emission constrain the radius of the progenitor star, its envelope composition, and the ratio of the ejecta energy to its mass, E/M. For He envelopes, neglecting the effect of recombination may lead to an over estimate of progenitor radius by more than an order of magnitude. We also show that the relative extinction at different wavelengths may be inferred from the light-curves at these wave-lengths, removing the uncertainty in the estimate of progenitor radius due to reddening (but not the uncertainty in E/M due to uncertainty in absolute extinction). The early UV/O observations of the type Ib SN2008D and of the type IIp SNLS-04D2dc are consistent with our model predictions. For SN2008D we find progenitor radius to be approx. 10^11 cm, and an indication that the He envelope contains a significant C/O fraction.

preprint2013arXiv

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

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.

preprint2013arXiv

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

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.

preprint2012arXiv

Imploding ignition waves: I. one dimensional analysis

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.

preprint2011arXiv

Early emission from type Ia supernovae

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.

preprint2011arXiv

Non-relativistic radiation mediated shock breakouts: I. Exact bolometric planar breakout solutions

The problem of a non-steady planar radiation mediated shock (RMS) breaking out from a surface with a power-law density profile, ρ\propto x^n, is numerically solved in the approximation of diffusion with constant opacity. For an appropriate choice of time, length and energy scales, determined by the breakout opacity, velocity and density, the solution is universal, i.e. depends only on the density power law index n. The resulting luminosity depends weakly on the value of n. An approximate analytic solution, based on the self-similar hydrodynamic solutions and on the steady RMS solutions, is constructed and shown to agree with the numerical solutions as long as the shock is far from the surface, τ>> c/v_{sh}. Approximate analytic expressions, calibrated based on the exact solutions, are provided, that describe the escaping luminosity as a function of time. These results can be used to calculate the bolometric properties of the bursts of radiation produced during supernova (SN) shock breakouts. For completeness, we also use the exact breakout solutions to provide an analytic approximation for the maximum surface temperature for fast (v_{sh}>~0.1) non-thermal breakouts, and show that it is few times smaller than inferred based on steady-state RMS solutions.

preprint2011arXiv

Non-relativistic radiation mediated shock breakouts: II. Bolometric properties of SN shock breakout

Exact bolometric light curves of supernova shock breakouts are derived based on the universal, non relativistic, planar breakout solutions (Sapir et al. 2011), assuming spherical symmetry, constant Thomson scattering opacity, κ, and angular intensity corresponding to the steady state planar limit. These approximations are accurate for progenitors with a scale height much smaller than the radius. The light curves are insensitive to the density profile and are determined by the progenitor radius R, and the breakout velocity and density, v_0 and ρ_0 respectively, and κ. The total breakout energy, E_BO, and the maximal ejecta velocity, v_max, are shown to be E_BO=8.0πR^2κ^-1cv_0 and v_max=2.0v_0 respectively, to an accuracy of about 10%. The calculated light curves are valid up to the time of transition to spherical expansion, t_sph\approx R/4v_0. Approximate analytic expressions for the light curves are provided for breakouts in which the shock crossing time at breakout, t_0=c/κρ_0v_0^2, is << R/c (valid for R<10^14 cm). Modifications of the flux angular intensity distribution and differences in shock arrival times to the surface, Δt_asym, due to moderately asymmetric explosions, affect the early light curve but do not affect v_max and E_BO. For 4v_0<<c, valid for large (RSG) progenitors, L\propto t^{-4/3} at max(Δt_asym,R/c)< t<t_sph and R may be accurately estimated from R\approx 2*10^13 (L/10^43 erg s^-1)^{2/5}(t/1 hr)^{8/15}.

preprint2011arXiv

On The GeV & TeV Detections of the Starburst Galaxies M82 & NGC 253

The GeV and TeV emission from M82 and NGC 253 observed by Fermi, HESS, and VERITAS constrains the physics of cosmic rays (CRs) in these dense starbursts. We argue that the gamma rays are predominantly hadronic in origin, as expected by previous studies. The measured fluxes imply that pionic losses are efficient for CR protons in both galaxies: we show that a fraction F_cal ~ 0.2 - 0.4 of the energy injected in high energy primary CR protons is lost to inelastic proton-proton collisions (pion production) before escape, producing gamma rays, neutrinos, and secondary electrons and positrons. We discuss the factor ~2 uncertainties in this estimate, including supernova rate and leptonic contributions to the GeV-TeV emission. We argue that gamma-ray data on ULIRGs like Arp 220 can test whether M82 and NGC 253 are truly calorimetric, and we present upper limits on Arp 220 from the Fermi data. We show that the observed ratio of the GeV to GHz fluxes of the starbursts suggests that non-synchrotron cooling processes are important for cooling the CR electron/positron population. We briefly reconsider previous predictions in light of the gamma-ray detections, including the starburst contribution to the gamma-ray background and CR energy densities. Finally, as a guide for future studies, we list the brightest star-forming galaxies on the sky and present updated predictions for their gamma-ray and neutrino fluxes.

preprint2011arXiv

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

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.

preprint2011arXiv

X-rays, gamma-rays and neutrinos from collisoinless shocks in supernova wind breakouts

We show that a collisionless shock necessarily forms during the shock breakout of a supernova (SN) surrounded by an optically thick wind. An intense non-thermal flash of <~ MeV gamma rays, hard X-rays and multi-TeV neutrinos is produced simultaneously with and following the soft X-ray breakout emission, carrying similar or larger energy than the soft emission. The non-thermal flash is detectable by current X-ray telescopes and may be detectable out to 10's of Mpc by km-scale neutrino telescopes.

preprint2010arXiv

Asymptotic self-similar solutions with a characteristic time-scale

For a wide variety of initial and boundary conditions, adiabatic one dimensional flows of an ideal gas approach self-similar behavior when the characteristic length scale over which the flow takes place, $R$, diverges or tends to zero. It is commonly assumed that self-similarity is approached since in the $R\to\infty(0)$ limit the flow becomes independent of any characteristic length or time scales. In this case the flow fields $f(r,t)$ must be of the form $f(r,t)=t^{α_f}F(r/R)$ with $R\propto(\pm t)^α$. We show that requiring the asymptotic flow to be independent only of characteristic length scales imply a more general form of self-similar solutions, $f(r,t)=R^{δ_f}F(r/R)$ with $\dot{R}\propto R^δ$, which includes the exponential ($δ=1$) solutions, $R\propto e^{t/τ}$. We demonstrate that the latter, less restrictive, requirement is the physically relevant one by showing that the asymptotic behavior of accelerating blast-waves, driven by the release of energy at the center of a cold gas sphere of initial density $ρ\propto r^{-ω}$, changes its character at large $ω$: The flow is described by $0\leδ<1$, $R\propto t^{1/(1-δ)}$, solutions for $ω<ω_c$, by $δ>1$ solutions with $R\propto (-t)^{1/(δ-1)}$ diverging at finite time ($t=0$) for $ω>ω_c$, and by exponential solutions for $ω=ω_c$ ($ω_c$ depends on the adiabatic index of the gas, $ω_c\sim8$ for $4/3<γ<5/3$). The properties of the new solutions obtained here for $ω\geω_c$ are analyzed, and self-similar solutions describing the $t>0$ behavior for $ω>ω_c$ are also derived.

preprint2010arXiv

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

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

preprint2010arXiv

Fast radiation mediated shocks and supernova shock breakouts

We present a simple analytic model for the structure of non-relativistic and relativistic radiation mediated shocks. At shock velocities β_s\equiv v_s/c\gtrsim 0.1, the shock transition region is far from thermal equilibrium, since the transition crossing time is too short for the production of a black-body photon density (by Bremsstrahlung emission). In this region, electrons and photons (and positrons) are in Compton (pair) equilibrium at temperatures T_s significantly exceeding the far downstream temperature, T_s\gg T_d\approx 2(\varepsilon n_u \hbar^3c^3)^{1/4}. T_s\gtrsim 10 keV is reached at shock velocities β_s\approx 0.2. At higher velocities, β_s\gtrsim0.6, the plasma is dominated in the transition region by e^\pm pairs and 60 keV\lesssim T_s \lesssim 200 keV. We argue that the spectrum emitted during the breaking out of supernova shocks from the stellar envelopes (or the surrounding winds) of Blue Super Giants and Wolf-Rayet stars, which reach β_s>0.1 for reasonable stellar parameters, may include a hard component with photon energies reaching tens or even hundreds of keV. This may account for the X-ray outburst associated with SN2008D, and possibly for other SN-associated outbursts with spectra not extending beyond few 100 keV (e.g. XRF060218/SN2006aj).

preprint2010arXiv

FERMI constraints on the high energy, ~1 GeV, emission of long GRBs

We investigate the constraints imposed on the luminosity function (LF) of long duration Gamma Ray Bursts (LGRBs) by the flux distribution of bursts detected by the GBM at ~1 MeV, and the implications of the non detection of the vast majority, ~95%, of the LGRBs at higher energy, ~1 GeV, by the LAT detector. We find a LF that is consistent with those determined by BATSE and Swift. The non detections by LAT set upper limits on the ratio R of the prompt fluence at ~1 GeV to that at ~1 MeV. The upper limits are more stringent for brighter bursts, with R<{0.1,0.3,1} for {5,30,60}% of the bursts. This implies that for most bursts the prompt ~1 GeV emission may be comparable to the ~1 MeV emission, but can not dominate it. The value of R is not universal, with a spread of (at least) an order of magnitude around R~10^(-1). For several bright bursts with reliable determination of the photon spectral index at ~1 MeV, the LAT non detection implies an upper limit to the ~100 MeV flux which is <0.1 of the flux obtained by extrapolating the ~1 MeV flux to high energy. For the widely accepted models, in which the ~1 MeV power-law photon spectrum reflects the power-law energy distribution of fast cooling electrons, this suggests that either the electron energy distribution does not follow a power-law over a wide energy range, or that the high energy photons are absorbed. Requiring an order unity pair production optical depth at ~100 MeV sets an upper limit for the Lorentz factor, Gamma<=10^(2.5).

preprint2010arXiv

Hard X-ray emission from accretion shocks around galaxy clusters

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

preprint2010arXiv

Long wavelength unstable modes in the far upstream of relativistic collisionless shocks

The growth rate of long wavelength kinetic instabilities arising due to the interaction of a collimated beam of relativistic particles and a cold unmagnetized plasma are calculated in the ultra relativistic limit. For sufficiently culminated beams, all long wave-length modes are shown to be Weibel-unstable, and a simple analytic expression for their growth rate is derived. For large transverse velocity spreads, these modes become stable. An analytic condition for stability is given. These analytic results, which generalize earlier ones given in the literature, are shown to be in agreement with numerical solutions of the dispersion equation and with the results of novel PIC simulations in which the electro-magnetic fields are restricted to a given k-mode. The results may describe the interaction of energetic cosmic rays, propagating into the far upstream of a relativistic collisionless shock, with a cold unmagnetized upstream. The long wavelength modes considered may be efficient in deflecting particles and could be important for diffusive shock acceleration. It is shown that while these modes grow in relativistic shocks propagating into electron-positron pair plasmas, they are damped in relativistic shocks propagating into electron-proton plasmas with moderate Lorenz factors Γ_{sh}\lesssim 100. If these modes dominate the deflection of energetic cosmic rays in electron-positron shocks, it is argued that particle acceleration is suppressed at shock frame energies that are larger than the downstream thermal energy by a factor greater than the shock Lorentz factor.

preprint2010arXiv

Magnetic fields and cosmic rays in clusters of galaxies

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

preprint2010arXiv

Nonthermal emission from clusters of galaxies

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

preprint2010arXiv

Relativistic Radiation Mediated Shocks

The structure of relativistic radiation mediated shocks (RRMS) propagating into a cold electron-proton plasma is calculated and analyzed. A qualitative discussion of the physics of relativistic and non relativistic shocks, including order of magnitude estimates for the relevant temperature and length scales, is presented. Detailed numerical solutions are derived for shock Lorentz factors $Γ_u$ in the range $6\leΓ_u\le30$, using a novel iteration technique solving the hydrodynamics and radiation transport equations (the protons, electrons and positrons are argued to be coupled by collective plasma processes and are treated as a fluid). The shock transition (deceleration) region, where the Lorentz factor $ Γ$ drops from $ Γ_u $ to $ \sim 1 $, is characterized by high plasma temperatures $ T\sim Γm_ec^2 $ and highly anisotropic radiation, with characteristic shock-frame energy of upstream and downstream going photons of a few~$\times\, m_ec^2$ and $\sim Γ^2 m_ec^2$, respectively.Photon scattering is dominated by e$^\pm$ pairs, with pair to proton density ratio reaching $\approx10^2Γ_u$. The width of the deceleration region, in terms of Thomson optical depths for upstream going photons, is large, $Δτ\simΓ_u^2$ ($Δτ\sim1$ neglecting the contribution of pairs) due to Klein Nishina suppression of the scattering cross section. A high energy photon component, narrowly beamed in the downstream direction, with a nearly flat power-law like spectrum, $νI_ν\proptoν^0$, and an energy cutoff at $ \sim Γ_u^2 m_ec^2 $ carries a fair fraction of the energy flux at the end of the deceleration region. An approximate analytic model of RRMS, reproducing the main features of the numerical results, is provided.

preprint2010arXiv

The local energy production rates of GRB photons and of UHECRs

In a recent analysis it was found that the local (z=0) rate at which gamma-ray bursts (GRBs) produce energy in 1 MeV photons, Q_GRB(z=0), is 300 times lower than the local energy production rate in ultra-high energy cosmic-rays. This may appear to be in contradiction with earlier results, according to which Q_GRB(z=0) is similar to the local energy production rate in >10^{19} eV cosmic-rays, Q_{10EeV}(z=0). This short (1 page) note identifies the origin of the apparent discrepancy and shows that Q_GRB(z=0) \sim Q_{10EeV}(z=0) holds.

preprint2010arXiv

The spectrum of Cosmic Rays escaping from relativistic shocks

We derive expressions for the time integrated spectrum of Cosmic Rays (CRs) that are accelerated in a decelerating relativistic shock wave and escape ahead of the shock. It is assumed that at any given time the CRs have a power law form, carry a constant fraction of the energy E_tot of the shocked plasma, and escape continuously at the maximal energy attainable. The spectrum of escaping particles is highly sensitive to the instantaneous spectral index due to the fact that the minimal energy, E_min ~ Γ^2 m_pc^2 where Γis the shock Lorentz factor, changes with time. In particular, the escaping spectrum may be considerably harder than the canonical N(E)\propto E^-2 spectrum. For a shock expanding into a plasma of density n, a spectral break is expected at the maximal energy attainable at the transition to non relativistic velocities, E ~ 10^19 (ε_B/0.1)(n/1 cm^-3)^(1/6)(E_tot/10^51 erg)^(1/3) eV where ε_B is the fraction of the energy flux carried by the magnetic field. If ultra-high energy CRs are generated in decelerating relativistic blast waves arising from the explosion of stellar mass objects, their generation spectrum may therefore be different than the canonical N(E)\propto E^-2.

preprint2010arXiv

What can we really learn from positron flux 'anomalies'?

We present a critical analysis of the observational constraints on, and of the theoretical modeling of, aspects of cosmic ray (CR) generation and propagation in the Galaxy, which are relevant for the interpretation of recent positron and anti-proton measurements. We give simple, analytic, model independent expressions for the secondary pbar flux, and an upper limit for the secondary e+ flux, obtained by neglecting e+ radiative losses, e+/(e+ + e-)<0.2\pm0.1 up to ~300 GeV. These expressions are completely determined by the rigidity dependent grammage, which is measured from stable CR secondaries up to ~150 GeV/nuc, and by nuclear cross sections measured in the laboratory. pbar and e+ measurements, available up to ~100 GeV, are consistent with these estimates, implying that there is no need for new, non-secondary, pbar or e+ sources. The radiative loss suppression factor f_{s,e+} of the e+ flux depends on the e+ propagation in the Galaxy, which is not understood theoretically. A rough, model independent estimate of f_{s,e+} 1/3 can be obtained at a single energy, E\sim20 GeV, from unstable secondary decay and is found to be consistent with e+ measurements, including the positron fraction measured by PAMELA. We show that specific detailed models, that agree with compositional CR data, agree with our simple expressions for the e+ and pbar flux, and that the claims that the positron fraction measured by PAMELA requires new primary e+ sources are based on assumptions, that are not supported by observations. If PAMELA results are correct, they suggest that f_{s,e+} is slightly increasing with energy, which provides an interesting constraint on CR propagation models. We argue that measurements of the e+ to pbar ratio are more useful for challenging secondary production models than the positron fraction.

preprint2009arXiv

The energy production rate & the generation spectrum of UHECRs

We derive simple analytic expressions for the flux and spectrum of ultra-high energy cosmic-rays (UHECRs) predicted in models where the CRs are protons produced by extra-Galactic sources. For a power-law scaling of the CR production rate with redshift and energy, d\dot{n} /dE\propto E^-α(1+z)^m, our results are accurate at high energy, E>10^18.7 eV, to better than 15%, providing a simple and straightforward method for inferring d\dot{n}/dE from the observed flux at E. We show that current measurements of the UHECR spectrum, including the latest Auger data, imply E^2d\dot{n}/dE(z=0)=(0.45\pm0.15)(α-1) 10^44 erg Mpc^-3 yr^-1 at E<10^19.5 eV with αroughly confined to 2\lesseqα<2.7. The uncertainty is dominated by the systematic and statistic errors in the experimental determination of individual CR event energy, (ΔE/E)_{sys} (ΔE/E)_{stat} ~20%. At lower energy, d\dot{n}/dE is uncertain due to the unknown Galactic contribution. Simple models in which α\simeq 2 and the transition from Galactic to extra-Galactic sources takes place at the "ankle", E ~10^19 eV, are consistent with the data. Models in which the transition occurs at lower energies require a high degree of fine tuning and a steep spectrum, α\simeq 2.7, which is disfavored by the data. We point out that in the absence of accurate composition measurements, the (all particle) energy spectrum alone cannot be used to infer the detailed spectral shapes of the Galactic and extra-Galactic contributions.

preprint2007arXiv

EeV neutrinos associated with UHECR sources

Electromagnetic energy losses of charged pions and muons suppress the expected high energy, >1E18 eV, neutrino emission from sources of ultrahigh energy, >1E19 eV, cosmic-rays. We show here that >1E19 eV photons produced in such sources by neutral pion decay may escape the sources, thanks to the Klein-Nishina suppression of the pair production cross section, and produce muon pairs in interactions with the cosmic microwave background. The flux of muon decay neutrinos, which are expected to be associated in time and direction with the electromagnetic emission from the sources, may reach a few percent of the Waxman-Bahcall bound. Their detection may allow one to directly identify the sources of >1E19 eV cosmic-rays, and will provide the most stringent constraints on quantum-gravity-induced Lorentz violation.

preprint2007arXiv

Properties of the Radio-Emitting Gas Around SgrA*

We show that the radial profiles of the temperature and density of the electrons as well as the magnetic field strength around the massive black hole at the Galactic center, SgrA*, may be constrained directly from existing radio data without any need to make prior assumptions about the dynamics of the emitting gas. The observed spectrum and wavelength-dependent angular size of SgrA* indicate that the synchrotron emission originates from an optically-thick plasma of quasi-thermal electrons. We find that the electron temperature rises above the virial temperature within tens of Schwarzschild radii from the black hole, suggesting that the emitting plasma may be outflowing. Constraints on the electron density profile are derived from polarization measurements. Our best-fit results differ from expectations based on existing theoretical models. However, these models cannot be ruled out as of yet due to uncertainties in the source size measurements. Our constraints could tighten considerably with future improvements in the size determination and simultaneous polarization measurements at multiple wavelengths.

preprint2006arXiv

Self-Similar Collisionless Shocks

Observations of gamma-ray burst afterglows suggest that the correlation length of magnetic field fluctuations downstream of relativistic non-magnetized collisionless shocks grows with distance from the shock to scales much larger than the plasma skin depth. We argue that this indicates that the plasma properties are described by a self-similar solution, and derive constraints on the scaling properties of the solution. For example, we find that the scaling of the characteristic magnetic field amplitude with distance from the shock is B \propto D^{s_B} with -1<s_B<=0, that the spectrum of accelerated particles is dn/dE \propto E^{-2/(s_B+1)}, and that the scaling of the magnetic correlation function is <B_i(x)B_j(x+dx)> \propto x^{2s_B} (for x>>D). We show that the plasma may be approximated as a combination of two self-similar components: a kinetic component of energetic particles and an MHD-like component representing "thermal" particles. We argue that the latter may be considered as infinitely conducting, in which case s_B=0 and the scalings are completely determined (e.g. dn/dE \propto E^{-2} and B \propto D^0). Similar claims apply to non- relativistic shocks such as in supernova remnants, if the upstream magnetic field can be neglected. Self-similarity has important implications for any model of particle acceleration and/or field generation. For example, we show that the diffusion function in the angle μof momentum p in diffusive shock acceleration models must satisfy D_{μμ}(p,D) = D^{-1}D'_{μμ}(p/D), and that a previously suggested model for the generation of large scale magnetic fields through a hierarchical merger of current-filaments should be generalized. A numerical experiment testing our analysis is outlined (Abridged).

preprint2006arXiv

The redshift distribution of short gamma-ray bursts from dynamically formed neutron star binaries

Short-hard gamma-ray bursts (SHBs) may arise from gravitational wave (GW) driven mergers of double neutron star (DNS) systems. DNSs may be "primordial" or can form dynamically by binary exchange interactions in globular clusters during core-collapse. For primordial binaries, the time delay between formation and merger is expected to be short, tau~0.1 Gyr, implying that the redshift distribution of merger events should follow that of star-formation. We point out here that for dynamically formed DNSs, the time delay between star-formation and merger is dominated by the cluster core-collapse time, rather than by the GW inspiral time, yielding delays comparable to the Hubble time. We derive the redshift distribution of merger events of dynamically formed DNSs, and find it to differ significantly from that typically expected for primordial binaries. The observed redshift distribution of SHBs favors dynamical formation, although a primordial origin cannot be ruled out due to possible detection biases. Future red-shift observations of SHBs may allow to determine whether they are dominated by primordial or dynamically formed DNSs.

preprint2006arXiv

The upstream magnetic field of GRB shocks

Gamma-ray burst (GRB) afterglow emission is believed to be produced by synchrotron emission of electrons accelerated to high energy by a relativistic collisionless shock propagating into a weakly magnetized plasma. Afterglow observations have been used to constrain the post-shock magnetic field and structure, as well as the accelerated electron energy distribution. Here we show that X-ray afterglow observations on day time scale constrain the pre-shock magnetic field to satisfy B>0.2[n/(1/cc)]^{5/8} mG, where n is the pre-shock density. This suggests that either the shock propagates into a highly magnetized fast, v~10^3 km/s, wind, or that the pre-shock magnetic field is strongly amplified, most likely by the streaming of high energy shock accelerated particles. More stringent constraints may be obtained by afterglow observations at high photon energy at late, >1 d, times.

preprint2005arXiv

A Complete Survey of the Transient Radio Sky and Implications for Gamma-Ray Bursts, Supernovae, and other Relativistic Explosions

We had previously reported on a survey for radio transients, used to set an upper limit on the number of orphan gamma-ray burst (GRB) radio afterglows, and thus a lower limit on the typical GRB beaming factor. Here we report radio and optical follow-up observations of these possible transients, achieving the first full characterization of the transient radio sky. We find that only two source are likely to be real radio transients, an optically obscured radio supernova (SN) in the nearby galaxy NGC 4216, and a source not associated with a bright host galaxy, which is too radio luminous to be a GRB afterglow. We speculate that this may be a flare from a peculiar active galactic nucleus, or a burst from an unusual Galactic compact object. We place an upper limit of 65 radio transients above 6 mJy over the entire sky at the 95% confidence level. The implications are as follows. First, we derive a limit on the typical beaming of GRBs; we find f_b^{-1} >~ 60, ~5 times higher than our earlier results. Second, we impose an upper limit on the rate of events that eject >~ 10^{51} erg in unconfined relativistic ejecta, whether or not accompanied by detectable emission in wavebands other than the radio. Our estimated rate, <=1000/y/Gpc, is about two orders of magnitude smaller than the rate of core-collapse SNe (and type Ib/c events in particular), indicating that only a minority of such events eject significant amounts of relativistic material, which are required by fireball models of long-soft GRBs. Finally, we show that future wider and/or deeper radio variability surveys are expected to detect numerous orphan radio GRB afterglows. Our survey also illustrates the great potential of sensitive surveys with new instruments to revolutionize the study of nearby SNe (abridged).

preprint2004arXiv

Probing the Magnetic Field Structure in Gamma-Ray Bursts through Dispersive Plasma Effects on the Afterglow Polarization

(Abr) The origin and structure of magnetic fields in Gamma-Ray Burst (GRB) fireball plasmas are two of the most important open questions in all GRB models. We show that the structure and strength of the magnetic field may be constrained by radio and IR observations of the early afterglow, where plasma effects on the polarization of propagating radiation are significant. We calculate these propagation effects for cold and relativistic plasmas, and find that in the presence of a uniform equipartition field the degree of linear polarization is suppressed, and circular polarization prevails at low frequencies, nu < 1-3 GHz, (2x10^11 Hz < nu < few x 10^14 Hz) in the forward (reverse) shock. At higher frequencies linear polarization dominates. At the frequency of the transition between circular and linear polarization, the net level of polarization is minimal, ~10-20%. These features are nearly independent of the circumburst density. The transition frequency is smaller by a factor of ~10 when the uniform field is much weaker than equipartition. The dependence of these results on viewing geometry, outflow collimation and magnetic field orientation is discussed. When the configuration of the field is entangled over length scales much smaller than the extent of the emitting plasma, the aforementioned effects should not be observed and a linear polarization at the few % level is expected. Polarimetric observations during the early afterglow, and particularly of the reverse shock emission, may therefore place strong constraints on the structure and strength of the magnetic field within the fireball plasma.

preprint2004arXiv

Prompt GRB spectra: detailed calculations and the effect of pair production

We present detailed calculations of the prompt spectrum of gamma-ray bursts (GRBs) predicted within the fireball model framework, where emission is due to internal shocks in an expanding relativistic wind. Our time dependent numerical model describes cyclo-synchrotron emission and absorption, inverse and direct Compton scattering, and e^\pm pair production and annihilation (including the evolution of high energy electro-magnetic cascades). It allows, in particular, a self-consistent calculation of the energy distribution of e^\pm pairs produced by photon annihilation, and hence a calculation of the spectra resulting when the scattering optical depth due to pairs, τ_\pm, is high. We show that emission peaks at ~1MeV for moderate to large τ_\pm, reaching τ_\pm ~ 10^2. In this regime of large compactness we find that (i) A large fraction of shock energy can escape as radiation even for large τ_\pm; (ii) The spectrum depends only weakly on the magnetic field energy fraction; (iii) The spectrum is hard, epsilon^2 dN/depsilon \propto epsilon^αwith 0.5<α<1, between the self absorption (epsilon_{ssa}= 10^{0.5\pm0.5} keV) and peak (epsilon_{peak}= 10^{0.5\pm0.5} MeV) photon energy, (iv) and shows a sharp cutoff at ~10 MeV; (v) Thermal Comptonization leads to emission peaking at epsilon_peak>=30 MeV, and can not therefore account for observed GRB spectra. For small compactness, spectra extend to >10 GeV with flux detectable by GLAST, and the spectrum at low energy depends on the magnetic field energy fraction. Comparison of the flux at ~1 GeV and ~ 100 keV may therefore allow to determine the magnetic field strength. For both small and large compactness, the spectra depend only weakly on the spectral index of the energy distribution of accelerated electrons.

preprint2003arXiv

Has the GZK suppression been discovered?

The energy spectra of ultra high energy cosmic rays reported by the AGASA, Fly's Eye, Haverah Park, HiRes, and Yakutsk experiments are all shown to be in agreement with each other for energies below 10^{20} eV (after small adjustments, within the known uncertainties, of the absolute energy scales). The data from HiRes, Fly's Eye, and Yakutsk are consistent with the expected flux suppression above 5\times 10^{19} eV due to interactions of cosmic rays with the cosmic microwave background, the Greisen-Zatsepin-Kuzmin (GZK) "supression," and are inconsistent with a smooth extrapolation of the observed cosmic ray energy spectrum to energies > 5\times 10^{19}$ eV. AGASA data show an excess of events above 10^{20} eV, compared to the predicted GZK suppression and to the flux measured by the other experiments.

preprint2000arXiv

Gamma-Ray Background from Structure Formation in the Intergalactic Medium

The universe is filled with a diffuse and isotropic extragalactic background of gamma-ray radiation, containing roughly equal energy flux per decade in photon energy between 3 MeV-100 GeV. The origin of this background is one of the unsolved puzzles in cosmology. Less than a quarter of the gamma-ray flux can be attributed to unresolved discrete sources, but the remainder appears to constitute a truly diffuse background whose origin has hitherto been mysterious. Here we show that the shock waves induced by gravity during the formation of large-scale structure in the intergalactic medium, produce a population of highly-relativistic electrons with a maximum Lorentz factor above 10^7. These electrons scatter a small fraction of the microwave background photons in the present-day universe up to gamma-ray energies, thereby providing the gamma-ray background. The predicted diffuse flux agrees with the observed background over more than four decades in photon energy, and implies a mean cosmological density of baryons which is consistent with Big-Bang nucleosynthesis.

preprint2000arXiv

High Energy Astrophysical Neutrinos: the Upper Bound is Robust

We elucidate the physical basis for the upper bound on high energy neutrino fluxes implied by the observed cosmic ray flux. We stress that the bound is valid for neutrinos produced either by p,gamma reactions or by p-p(n) reactions in sources which are optically thin for high energy protons to photo-meson and nucleon-meson interactions. We show that the upper bound is robust and conservative. The Waxman-Bahcall bound overestimates the most likely neutrino flux by a factor ~ 5/tau, for small optical depths tau. The upper limit cannot be plausibly evaded by invoking magnetic fields, optically thick AGNs, or large hidden fluxes of extragalactic protons. We describe the implications of the bound for future experiments including the AMANDA, ANTARES, Auger, ICECUBE, NESTOR, and OWL/AIRWATCH detectors.

preprint1998arXiv

Can high energy neutrino annihilation on relic neutrinos generate the observed highest energy cosmic-rays?

Annihilation of high energy, $\sim 10^{21}$eV, neutrinos on big bang relic neutrinos of $\sim 1$eV mass, clustered in the Galactic halo or in a nearby galaxy cluster halo, has been suggested to generate, through hadronic Z decay, high energy nucleons and photons which may account for the detected flux of >10^{20}eV cosmic-rays. We show that the flux of high energy nucleons and photons produced by this process is dominated by annihilation on the uniform, non-clustered, neutrino background, and that the energy generation rate of 10^{21}eV neutrinos required to account for the detected flux of >10^{20}eV particles is >10^{48} erg/Mpc^3 yr. This energy generation rate, comparable to the total luminosity of the universe, is 4 orders of magnitude larger than the rate of production of high energy nucleons required to account for the flux of >10^{19}eV cosmic-rays. Thus, in order for neutrino annihilation to contribute significantly to the detected flux of >10^{20}eV cosmic-rays, the existence of a new class of high-energy neutrino sources, likely unrelated to the sources of >10^{19}eV cosmic-rays, must be invoked.

preprint1998arXiv

Gamma-Ray Burst Afterglow: Polarization and Analytic Light Curves

GRB afterglow polarization is discussed. We find an observable, up to 10%, polarization, if the magnetic field coherence length grows at about the speed of light after the field is generated at the shock front. Detection of a polarized afterglow would show that collisionless ultrarelativistic shocks can generate strong large scale magnetic fields and confirm the synchrotron afterglow model. Non-detection, at a 1% level, would imply that either the synchrotron emission model is incorrect, or that strong magnetic fields, after they are generated in the shock, somehow manage to stay un-dissipated at ``microscopic'', skin depth, scales. Analytic lightcurves of synchrotron emission from an ultrarelativistic self-similar blast wave are obtained for an arbitrary electron distribution function, taking into account the effects of synchrotron cooling. The peak synchrotron flux and the flux at frequencies much smaller than the peak frequency are insensitive to the details of the electron distribution function; hence their observational determination would provide strong constraints on blast wave parameters.

preprint1998arXiv

High Energy Neutrinos from Astrophysical Sources: An Upper Bound

We show that cosmic-ray observations set a model-independent upper bound to the flux of high-energy, > 10^14 eV, neutrinos produced by photo-meson (or p-p) interactions in sources of size not much larger than the proton photo-meson (or pp) mean-free-path. The bound applies, in particular, to neutrino production by either AGN jets or GRBs. This upper limit is two orders of magnitude below the flux predicted in some popular AGN jet models, but is consistent with our predictions from GRB models. We discuss the implications of these results for future km^2 high-energy neutrino detectors.

preprint1997arXiv

Gamma-ray burst after-glow: Confirming the cosmological fireball model

The recent detection of delayed X-ray and optical emission, ``afterglow,'' associated with gamma-ray bursts (GRBs) supports models, where the bursts are produced by relativistic expanding blastwaves, ``fireballs,'' at cosmological distances. The detection of absorption lines in the optical afterglow of the GRB of 8 May 1997 confirms that the sources lie at cosmological distance. We show here that the new features detected in GRB970508 afterglow, radio emission one week following the burst and a 2 day increase in optical flux, are consistent with the blastwave model. The fireball optical depth at radio frequencies is much smaller than previously estimated, which accounts for the observed radio emission. The initial suppression of optical flux is consistent with that predicted due to electron cooling. The combined radio and optical data imply that the fireball energy is $\sim10^{52}$erg, and that the density of the medium into which the blastwave expands is $\sim 1cm^{-3}$, a value typical for gas within galaxies. We predict the time dependence of the radio flux and the absorption frequency, which constitute tests of the fireball model as described in this paper.

preprint1997arXiv

Implications of the radio afterglow from the gamma-ray burst of May 8, 1997

Radio observations of the afterglow of the gamma-ray burst GRB970508 provide unique new constraints on afterglow models. The quenching of diffractive scintillation at 4 week delay provides the first direct estimate of source size and expansion rate. It implies an apparent size $R\sim10^{17}$cm and expansion at a speed comparable to that of light at 4 weeks delay, in agreement with the fireball model prediction $R=10^{17}(t/week)^{5/8}$cm. The radio flux and its dependence on time and frequency at 1--5 week delay are in agreement with the model and imply a fireball energy (assuming spherical symmetry) of order 10^{52}erg, consistent with the value inferred from observations at shorter delay. The observed radio behavior deviates from model predictions at delays larger than 5 weeks. This is expected, since at this delay the fireball is in transition from highly-relativistic to sub-relativistic expansion, with Lorentz factor γ<2. Deviation may be due to a change in the physical processes associated with the shock wave as it becomes sub-relativistic (e.g. a decrease in the fraction of energy carried by magnetic field), or to the fireball being a cone of opening angle $1/γ\sim0.5$. We predict the future behavior of the radio flux assuming that the latter interpretation is valid. These predictions may be tested by radio observations in the frequency range 0.1--10GHz on time scale of months.