No Arabic abstract
We present theoretical expectations for non-thermal emission due to the bulk Comptonization at the ultra-relativistic shock breakout. We calculate the transfer of photons emitted from the shocked matter with a Monte Carlo code fully taking into account special relativity. As a hydrodynamical model, we use a self-similar solution of Nakayama & Shigeyama (2005). Our calculations reveal that the spectral shape exhibits a double peak or a single peak depending on the shock temperature at the shock breakout. If it is significantly smaller than the rest energy of an electron, the spectrum has a double peak. We also display a few example of light curves, and estimate the total radiation energy. In comparison with observations of gamma-ray bursts, a part of the higher energy component in the spectra and the total energy can be reproduced by some parameter sets. Meanwhile, the lower energy counterpart in the Band function is not reproduced by our results and the duration time seems too short to represent the entire event of a gamma-ray burst. Therefore the subsequent phase will constitute the lower energy part in the spectrum.
The brief transient emitted as a shock wave erupts through the surface of a presupernova star carries information about the stellar radius and explosion energy. Here the CASTRO code, which treats radiation transport using multigroup flux-limited diffusion, is used to simulate the light curves and spectra of shock breakout in very low-energy supernovae (VLE SNe), explosions in giant stars with final kinetic energy much less than 10$^{51}$ erg. VLE SNe light curves, computed here with the KEPLER code, are distinctively faint, red, and long-lived, making them challenging to find with transient surveys. The accompanying shock breakouts are brighter, though briefer, and potentially easier to detect. Previous analytic work provides general guidance, but numerical simulations are challenging due to the range of conditions and lack of equilibration between color and effective temperatures. We consider previous analytic work and extend discussions of color temperature and opacity to the lower energy range explored by these events. Since this is the first application of the CASTRO code to shock breakout, test simulations of normal energy shock breakout of SN1987A are carried out and compared with the literature. A set of breakout light curves and spectra are then calculated for VLE SNe with final kinetic energies in the range $10^{47} - 10^{50}$ ergs for red supergiants with main sequence masses 15 Msun and 25 Msun. The importance of uncertainties in stellar atmosphere model, opacity, and ambient medium is discussed, as are observational prospects with current and forthcoming missions.
Shock breakout is the brightest radiative phenomenon in a supernova (SN) but is difficult to be observed owing to the short duration and X-ray/ultraviolet (UV)-peaked spectra. After the first observation from the rising phase reported in 2008, its observability at high redshift is attracting enormous attention. We perform multigroup radiation hydrodynamics calculations of explosions for evolutionary presupernova models with various main-sequence masses $M_{rm MS}$, metallicities $Z$, and explosion energies $E$. We present multicolor light curves of shock breakout in Type II plateau SNe, being the most frequent core-collapse SNe, and predict apparent multicolor light curves of shock breakout at various redshifts $z$. We derive the observable SN rate and reachable redshift as functions of filter $x$ and limiting magnitude $m_{x,{rm lim}}$ by taking into account an initial mass function, cosmic star formation history, intergalactic absorption, and host galaxy extinction. We propose a realistic survey strategy optimized for shock breakout. For example, the $g$-band observable SN rate for $m_{g,{rm lim}}=27.5$ mag is 3.3 SNe degree$^{-2}$ day$^{-1}$ and a half of them locates at $zgeq1.2$. It is clear that the shock breakout is a beneficial clue to probe high-$z$ core-collapse SNe. We also establish ways to identify shock breakout and constrain SN properties from the observations of shock breakout, brightness, time scale, and color. We emphasize that the multicolor observations in blue optical bands with $sim$ hour intervals, preferably over $geq2$ continuous nights, are essential to efficiently detect, identify, and interpret shock breakout.
We present 3D simulations of core-collapse supernovae from blast-wave initiation by the neutrino-driven mechanism to shock breakout from the stellar surface, considering two 15 Msun red supergiants (RSG) and two blue supergiants (BSG) of 15 Msun and 20 Msun. We demonstrate that the metal-rich ejecta in homologous expansion still carry fingerprints of asymmetries at the beginning of the explosion, but the final metal distribution is massively affected by the detailed progenitor structure. The most extended and fastest metal fingers and clumps are correlated with the biggest and fastest-rising plumes of neutrino-heated matter, because these plumes most effectively seed the growth of Rayleigh-Taylor (RT) instabilities at the C+O/He and He/H composition-shell interfaces after the passage of the SN shock. The extent of radial mixing, global asymmetry of the metal-rich ejecta, RT-induced fragmentation of initial plumes to smaller-scale fingers, and maximal Ni and minimal H velocities do not only depend on the initial asphericity and explosion energy (which determine the shock and initial Ni velocities) but also on the density profiles and widths of C+O core and He shell and on the density gradient at the He/H transition, which lead to unsteady shock propagation and the formation of reverse shocks. Both RSG explosions retain a great global metal asymmetry with pronounced clumpiness and substructure, deep penetration of Ni fingers into the H-envelope (with maximum velocities of 4000-5000 km/s for an explosion energy around 1.5 bethe) and efficient inward H-mixing. While the 15 Msun BSG shares these properties (maximum Ni speeds up to ~3500 km/s), the 20 Msun BSG develops a much more roundish geometry without pronounced metal fingers (maximum Ni velocities only ~2200 km/s) because of reverse-shock deceleration and insufficient time for strong RT growth and fragmentation at the He/H interface.
The maximum energy to which cosmic rays can be accelerated at weakly-magnetised ultra-relativistic shocks is investigated. We demonstrate that for such shocks, in which the scattering of energetic particles is mediated exclusively by ion skin-depth scale structures, as might be expected for a Weibel-mediated shock, there is an intrinsic limit on the maximum energy to which particles can be accelerated. This maximum energy is determined from the requirement that particles must be isotropised in the downstream plasma frame before the mean field transports them far downstream, and falls considerably short of what is required to produce ultra-high-energy cosmic rays. To circumvent this limit, a highly disorganised field is required on larger scales. The growth of cosmic-ray induced instabilities on wavelengths much longer than the ion-plasma skin depth, both upstream and downstream of the shock, is considered. While these instabilities may play an important role in magnetic field amplification at relativistic shocks, on scales comparable to the gyroradius of the most energetic particles, the calculated growth-rates have insufficient time to modify the scattering. Since strong modification is a necessary condition for particles in the downstream region to re-cross the shock, in the absence of an alternative scattering mechanism, these results imply that acceleration to higher energies is ruled out. If weakly-magnetised ultra-relativistic shocks are disfavoured as high energy particle accelerators in general, the search for potential sources of ultra-high-energy cosmic rays can be narrowed.
Ultra-high energy (UHE) photons play an important role as an independent probe of the photo-pion production mechanism by UHE cosmic rays. Their observation, or non-observation, may constrain astrophysical scenarios for the origin of UHECRs and help to understand the nature of the flux suppression observed by several experiments at energies above $10^{19.5}$ eV. Whereas the interaction length of UHE photons above $10^{17}$ eV ranges from a few hundred kpc up to tenths of Mpc, photons can interact with the extragalactic background radiation initiating the development of electromagnetic cascades which affect the fluxes of photons observed at Earth. The interpretation of the current experimental results rely on the simulations of the UHE photon propagation. In this paper, we present the novel Monte Carlo code EleCa to simulate the $Ele$ctromagnetic $Ca$scading initiated by high-energy photons and electrons. We provide an estimation of the surviving probability for photons inducing electromagnetic cascades as a function of their distance from the observer and we calculate the distances within which we expect to observe UHE photons with energy between $10^{17}$ and $10^{19}$ eV. Furthermore, the flux of GZK photons at Earth is investigated in several astrophysical scenarios where we vary both injection spectrum and composition at the source and the intensity of the intervening extragalactic magnetic field. Although the photon propagation depends on several astrophysical factors, our numerical predictions combined with future experimental observations (or non-observations) of UHE photons -- in the energy range between $10^{17.5}$ eV and $10^{20}$ eV -- can help to constrain these scenarios.