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})^\alpha]$ 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 $\alpha=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$.

With the advent of new wide-field, high-cadence optical transient surveys, our understanding of the diversity of core-collapse supernovae has grown tremendously in the last decade. However, the pre-supernova evolution of massive stars, that sets the physical backdrop to these violent events, is theoretically not well understood and difficult to probe observationally. Here we report the discovery of the supernova iPTF13dqy = SN 2013fs, a mere ~3 hr after explosion. Our rapid follow-up observations, which include multiwavelength photometry and extremely early (beginning at ~6 hr post-explosion) spectra, map the distribution of material in the immediate environment (<~ 10^15 cm) of the exploding star and establish that it was surrounded by circumstellar material (CSM) that was ejected during the final ~1 yr prior to explosion at a high rate, around 10^-3 solar masses per year. The complete disappearance of flash-ionised emission lines within the first several days requires that the dense CSM be confined to within <~ 10^15 cm, consistent with radio non-detections at 70--100 days. The observations indicate that iPTF13dqy was a regular Type II SN; thus, the finding that the probable red supergiant (RSG) progenitor of this common explosion ejected material at a highly elevated rate just prior to its demise suggests that pre-supernova instabilities may be common among exploding massive stars.

We present here an efficient numerical scheme for solving the non-relativistic 1D radiation-hydrodynamics equations including inelastic Compton scattering, which is not included in most codes and is crucial for solving problems such as shock breakout. The devised code is applied to the problems of a steady-state planar radiation mediated shock (RMS) and RMS breakout from a stellar envelope. The results are in agreement with those of a previous work on shock breakout \citep{Sapir13}, in which Compton equilibrium between matter and radiation was assumed and the "effective photon" approximation was used to describe the radiation spectrum. In particular, we show that the luminosity and its temporal dependence, the peak temperature at breakout, and the universal shape of the spectral fluence derived in this earlier work are all accurate. Although there is a discrepancy between the spectral calculations and the effective photon approximation due to the inaccuracy of the effective photon approximation estimate of the effective photon production rate, which grows with lower densities and higher velocities, the difference in peak temperature reaches only $30\%$ for the most discrepant cases of fast shocks in blue supergiants. The incompatibility of the stellar envelope shock breakout model results with observed properties of X-ray flashes and the discrepancy between the predicted and observed rates of X-ray flashes \citep{Sapir13} remain unexplained.

The spectrum of radiation emitted following shock breakout from a star's surface with a power-law density profile $\rho \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 $\rho_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 $\rho_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.

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.

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, \kappa, 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 \rho_0 respectively, and \kappa. The total breakout energy, E_BO, and the maximal ejecta velocity, v_max, are shown to be E_BO=8.0\pi R^2\kappa^-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/\kappa\rho_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, \Delta 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(\Delta 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}.

The problem of a non-steady planar radiation mediated shock (RMS) breaking out from a surface with a power-law density profile, \rho\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, \tau>> 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.