Shock breakout (SBO) signal is the first signature of the supernova explosion apart from gravitational waves and neutrinos. Observational properties of SBO, such as bolometric luminosity and colour temperature, connect to the supernova progenitor and explosion parameters. Detecting SBO or SBO-cooling will constrain the progenitor and explosion models of collapsing stars. In the light of recently launched eROSITA telescope, the rate for detection of SBO is a few events during a year. In the current study, we examine the analytic formulae derived by Shussman et al. (2016). We use four red supergiant models from their study, while running explosions with the radiation hydrodynamics code STELLA. We conclude that there is a good agreement between analytic and numerical approaches for bolometric luminosity and colour temperature during SBO. The analytic formulae for the SBO signal based on the global supernova parameters can be used instead of running time-consuming numerical simulations. We define spectral range where analytic formulae for the SBO spectra are valid. We provide improved analytical expression for the SBO spectral energy distribution. We confirm dependence of colour temperature on radius derived by analytical studies and suggest to use early time observations to confine the progenitor radius. Additionally we show the prediction for the SBO signal from red supergiants as seen by eROSITA instrument.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

In the present study, we systematically explore the effect of the radioactive 56Ni and its mixing properties in the ejecta on the plateau of Type IIP supernovae (SNe). We evaluate the importance of 56Ni in shaping light curves of SNe IIP by simulating light curves for two red supergiant models using different amounts of 56Ni and with different types of mixing: uniform distribution of 56Ni out to different fractions of the envelope and "boxcar" distribution of 56Ni. We find, similarly to previous studies, that 56Ni extends duration of the plateau. We find a formula to estimate the extension based on the observed bolometric light curves and show that for most SNe IIP 56Ni extends the plateau by about 20%. Another effect of 56Ni consists in reduction of the plateau decline rate, i.e. 56Ni presented in the ejecta flattens the plateau. Our simulations suggest that for typical SNe IIP it can reduce the decline rate by about 1 mag/100 day. We find that for the contribution of 56Ni seen in most SNe our simulated bolometric light curves resemble observed ones for various types of 56Ni mixing. We thereby cannot determine the level of 56Ni mixing in these SNe based on the light curve alone. However, for SN2009ib we find that only a model where 56Ni is mixed significantly throughout most of the hydrogen envelope is consistent with the observed light curve. Our light curves are available via link https://wwwmpa.mpa-garching.mpg.de/ccsnarchive/data/Kozyreva2018/.

Massive stars with a core-halo structure are interesting objects for stellar physics and hydrodynamics. Using simulations for stellar evolution, radiation hydrodynamics, and radiative transfer, we study the explosion of stars with an extended and tenuous envelope (i.e., stars in which 95% of the mass is contained within 10% of the surface radius or less). We consider both H-rich supergiant and He-giant progenitors resulting from close-binary evolution and dying with a final mass of 2.8-5Msun. An extended envelope causes the supernova (SN) shock to brake and a reverse shock to form, sweeping core material into a dense shell. The shock deposited energy, which suffers little degradation from expansion, is trapped in ejecta layers of moderate optical depth, thereby enhancing the SN luminosity at early times. With the delayed 56Ni heating, we find that the resulting optical and near-IR light curves all exhibit a double-peak morphology. We show how an extended progenitor can explain the blue and featureless optical spectra of some Type IIb and Ib SNe. The dense shell formed by the reverse shock leads to line profiles with a smaller and near-constant width. This ejecta property can explain the statistically narrower profiles of Type IIb compared to Type Ib SNe, as well as the peculiar Halpha profile seen in SN1993J. At early times, our He-giant star explosion model shows a high luminosity, a blue color, and featureless spectra reminiscent of the Type Ib SN2008D, suggesting a low-mass progenitor.

The progenitors of low-luminosity Type II-Plateau supernovae (SNe II-P) are believed to be red supergiant (RSG) stars, but there is much disparity in the literature concerning their mass at core collapse and therefore on the main sequence. Here, we model the SN radiation arising from the low-energy explosion of RSG stars of 12, 25, and 27 M$_{\odot}$ on the main sequence and formed through single star evolution. Despite the narrow range in ejecta kinetic energy (2.5$-$4.2$\times$10$^{50}$ erg) in our model set, the SN observables from our three models are significantly distinct, reflecting the differences in progenitor structure (e.g., surface radius, H-rich envelope mass, He-core mass). Our higher mass RSG stars give rise to Type II SNe that tend to have bluer colors at early times, a shorter photospheric phase, and a faster declining $V$-band light curve (LC) more typical of Type II-linear SNe, in conflict with the LC plateau observed for low-luminosity SNe II. The complete fallback of the CO core in the low-energy explosions of our high mass RSG stars prevents the ejection of any ${}^{56}$Ni (nor any core O or Si), in contrast to low-luminosity SNe II-P, which eject at least 0.001 M$_{\odot}$ of ${}^{56}$Ni. In contrast to observations, type II SN models from higher mass RSGs tend to show an H$\alpha$ absorption that remains broad at late times (due to a larger velocity at the base of the H-rich envelope). In agreement with the analyses of pre-explosion photometry, we conclude that low-luminosity SNe II-P likely arise from low-mass rather than high-mass RSG stars.

Zahn (1975) first put forward and calculated in detail the torque experienced by stars in a close binary systems due to dynamical tides. His widely used formula for stars with radiative envelopes and convective cores is expressed in terms of the stellar radius, even though the torque is actually being applied to the convective core at the core radius. This results in a large prefactor, which is very sensitive to the global properties of the star, that multiplies the torque. This large factor is compensated by a very small multiplicative factor, $E_{2}$. Although this is mathematically accurate, depending on the application this can lead to significant errors. The problem is even more severe, since the calculation of $E_{2}$ itself is non-trivial, and different authors have obtained inconsistent values of $E_{2}$. Moreover, many codes (e.g. BSE, StarTrack, MESA) interpolate (and sometimes extrapolate) a fit of $E_{2}$ values to the stellar mass, often in regimes where this is not sound practice. We express the torque in an alternate form, cast in terms of parameters at the envelope-core boundary and a dimensionless coefficient, $\beta_{2}$. Previous attempts to express the torque in such a form are either missing an important factor, which depends on the density profile of the star, or are not easy to implement. We show that $\beta_{2}$ is almost independent of the properties of the star and its value is approximately unity. Our formula for the torque is simple to implement and avoids the difficulties associated with the classic expression.

Using 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations, we study the ejecta properties required to match the early and late-time photometric and spectroscopic properties of supernovae (SNe) associated with long-duration gamma-ray bursts (LGRBs). To match the short rise time, narrow light curve peak, and extremely broad spectral lines of SN1998bw requires a model with <3Msun ejecta but a high explosion energy of a few 1e52erg and 0.5Msun of Ni56. However the relatively high luminosity, the presence of narrow spectral lines of intermediate mass elements, and the low ionization at the nebular stage are matched with a more standard C-rich Wolf-Rayet (WR) star explosion, with an ejecta of >10Msun, an explosion energy >1e51erg, and only 0.1Msun of Ni56. As the two models are mutually exclusive, the breaking of spherical symmetry is essential to match the early/late photometric/spectroscopic properties of SN1998bw. This conclusion confirms the notion that the ejecta of SN1998bw is aspherical on large scales. More generally, with asphericity, the energetics and Ni56 mass of LGRB/SNe are reduced and their ejecta mass is increased, favoring a massive fast-rotating Wolf-Rayet star progenitor. Contrary to persisting claims in favor of the proto-magnetar model for LGRB/SNe, such progenitor/ejecta properties are compatible with collapsar formation. Ejecta properties of LGRB/SNe inferred from 1D radiative-transfer modeling are fundamentally flawed.

We present the light curves of the hydrogen-poor superluminous supernovae (SLSNe-I) PTF12dam and iPTF13dcc, discovered by the (intermediate) Palomar Transient Factory. Both show excess emission at early times and a slowly declining light curve at late times. The early bump in PTF12dam is very similar in duration (~10 days) and brightness relative to the main peak (2-3 mag fainter) compared to those observed in other SLSNe-I. In contrast, the long-duration (>30 days) early excess emission in iPTF13dcc, whose brightness competes with that of the main peak, appears to be of a different nature. We construct bolometric light curves for both targets, and fit a variety of light-curve models to both the early bump and main peak in an attempt to understand the nature of these explosions. Even though the slope of the late-time light-curve decline in both SLSNe is suggestively close to that expected from the radioactive decay of $^{56}$Ni and $^{56}$Co, the amount of nickel required to power the full light curves is too large considering the estimated ejecta mass. The magnetar model including an increasing escape fraction provides a reasonable description of the PTF12dam observations. However, neither the basic nor the double-peaked magnetar model is capable of reproducing the iPTF13dcc light curve. A model combining a shock breakout in an extended envelope with late-time magnetar energy injection provides a reasonable fit to the iPTF13dcc observations. Finally, we find that the light curves of both PTF12dam and iPTF13dcc can be adequately fit with the circumstellar medium (CSM) interaction model.

Supernova (SN) 2008bk is a well observed low-luminosity Type II event visually associated with a low-mass red-supergiant progenitor. To model SN 2008bk, we evolve a 12Msun star from the main sequence until core collapse, when it has a total mass of 9.88Msun, a He-core mass of 3.22Msun, and a radius of 502Rsun. We then artificially trigger an explosion that produces 8.29Msun of ejecta with a total energy of 2.5x10^50erg and ~0.009Msun of 56Ni. We model the subsequent evolution of the ejecta with non-Local-Thermodynamic-Equilibrium time-dependent radiative transfer. Although somewhat too luminous and energetic, this model reproduces satisfactorily the multi-band light curves and multi-epoch spectra of SN 2008bk, confirming the suitability of a low-mass massive star progenitor. As in other low-luminosity SNe II, the structured Halpha profile at the end of the plateau phase is probably caused by BaII 6496.9A rather than asphericity. We discuss the sensitivity of our results to changes in progenitor radius and mass, as well as chemical mixing. A 15% increase in progenitor radius causes a 15% increase in luminosity and a 0.2mag V-band brightening of the plateau but leaves its length unaffected. An increase in ejecta mass by 10% lengthens the plateau by ~10d. Chemical mixing introduces slight changes to the bolometric light curve, limited to the end of the plateau, but has a large impact on colours and spectra at nebular times.

Observations of type II supernova early light, from breakout until recombination, can be used to constrain the explosion energy and progenitor properties. Currently available for this purpose are purely analytic models, which are accurate only to within an order of magnitude, and detailed numerical simulations, which are more accurate but are applied to any event separately. In this paper we derive an analytic model that is calibrated by numerical simulations. This model is much more accurate than previous analytic models, yet it is as simple to use. To derive the model we analyze simulated light curves from numerical explosion of $124$ red supergiant progenitors, calculated using the stellar evolution code MESA. We find that although the structure of the progenitors we consider varies, the resulting light curves can be described rather well based only on the explosion energy, ejecta mass and progenitor radius. Our calibrated analytic model, which is based on these three parameters, reproduces the bolometric luminosity within $25\%-35\%$ accuracy and the observed temperature within $15\%$ accuracy (compared to previous analytic models which are indeed found to be accurate only to within an order of magnitude). We also consider deviations of the early time spectrum from blackbody, and find that the Rayleigh-Jeans regime is slightly shallower (roughly $L_\nu \propto \nu^{1.4}$). This modified spectrum affects the optical/near-UV light curve mostly during the first day when the typical observed temperature is $\gg 10^4 ~^\circ$K. We use our results to study the optical and near-UV early light curves from first light until recombination and briefly discuss what can be learned from current and future observations. Light curves generated using our calibrated model can be downloaded at http://www.astro.tau.ac.il/~tomersh/.

With an increasing number of superluminous supernovae (SLSNe) discovered the question of their origin remains open and causes heated debates in the supernova community. Currently, there are three proposed mechanisms for SLSNe: (1) pair-instability supernovae (PISN), (2) magnetar-driven supernovae, and (3) models in which the supernova ejecta interacts with a circumstellar material ejected before the explosion. Based on current observations of SLSNe, the PISN origin has been disfavoured for a number of reasons. Many PISN models provide overly broad light curves and too reddened spectra, because of massive ejecta and a high amount of nickel. In the current study we re-examine PISN properties using progenitor models computed with the GENEC code. We calculate supernova explosions with FLASH and light curve evolution with the radiation hydrodynamics code STELLA. We find that high-mass models (200 and 250 solar masses) at relatively high metallicity (Z=0.001) do not retain hydrogen in the outer layers and produce relatively fast evolving PISNe Type I and might be suitable to explain some SLSNe. We also investigate uncertainties in light curve modelling due to codes, opacities, the nickel-bubble effect and progenitor structure and composition.

We explore the implications of the observed low spin of GW150914 within the context of stellar astrophysics and progenitor models. We conclude that many of the recently proposed scenarios are in marked tension with this observation. We derive a simple model for the observed spin in the case that the progenitor system was a field binary composed of a black hole (BH) and a Wolf-Rayet star and explore the implications of the observed spin for this model. The spin observation allows us to place a lower limit for the delay time between the formation of the BH+BH binary and the actual merger, $t_{\textrm{merge}}$. We use typical values for these systems to derive $t_{\textrm{merge}}\gtrsim10^{8}\,\textrm{yr}$, which proves to be an important diagnostic for different progenitor models. We anticipate the next series of events, and the associated spin parameters, will ultimately yield critical constraints on formation scenarios and on stellar parameters describing the late-stage evolution of massive stars.

It was recently shown that the bolometric light curves of type II supernovae (SNe) allow an accurate and robust measurement of the product of the radiation energy in the ejecta, $E_r$, and the time since the explosion, $t$, at early phases ($t\lesssim 10d$) of the homologous expansion. This observable, denoted here $ET \equiv E_rt$ is constant during that time and depends only on the progenitor structure and explosion energy. We use a 1D hydrodynamic code to find $ET$ of simulated explosions of 145 red supergiant progenitors obtained using the stellar evolution code MESA, and relate this observable to the properties of the progenitor and the explosion energy. We show that $ET$ probes only the properties of the envelope (velocity, mass and initial structure), similarly to other observables that rely on the photospheric phase emission. Nevertheless, for explosions where the envelope dominates the ejected mass, $M_{env}/M_{ej} \gtrsim 0.6$, $ET$ is directly related to the explosion energy $E_{exp}$ and ejected mass $M_{ej}$ through the relation $ET \approx 0.15 E_{exp}^{1/2} R_* M_{ej}^{1/2}$, where $R_*$ is the progenitor radius, to an accuracy better than $30\%$. We also provide relations between $ET$ and the envelope properties that are accurate (to within 20\%) for all the progenitors in our sample, including those that lost most of their envelope. We show that when the envelope velocity can be reasonably measured by line shifts in observed spectra, the envelope is directly constrained from the bolometric light curve (independent of $E_{exp}$). We use that to compare observations of 11 SNe with measured $ET$ and envelope velocity to our sample of numerical progenitors. This comparison suggests that many SNe progenitors have radii that are $\lesssim 500~R_\odot$. In the framework of our simulations this indicates, most likely, a rather high value of the mixing length parameter.

We present 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations for a large grid of supernovae (SNe) IIb/Ib/Ic that result from the terminal explosion of the mass donor in a close-binary system. Our sample covers ejecta masses $M_{\rm e}$ of 1.7$-$5.2M$_\odot$, kinetic energies $E_{\rm kin}$ of 0.6$-$5.0$\times$10$^{51}$erg, and $^{56}$Ni masses of 0.05$-$0.30M$_\odot$. We find a strong correlation between the $^{56}$Ni mass and the photometric properties at maximum, and between the rise time to bolometric maximum and the post-maximum decline rate. We confirm the small scatter in ($V-R$) at 10d past $R$-band maximum. The quantity $V_{\rm m} \equiv \sqrt{2E_{\rm kin}/M_{\rm e}}$ is comparable to the Doppler velocity measured from HeI 5875\AA\ at maximum in SNe IIb/Ib, although some scatter arises from the uncertain level of chemical mixing. The OI7772\AA\ line may be used for SNe Ic, but the correspondence deteriorates with higher ejecta mass/energy. We identify a temporal reversal of the Doppler velocity at maximum absorption in the $\sim$1.05$\mu$m feature in all models. The reversal is due to HeI alone and could serve as a test for the presence of helium in SNe Ic. Because of variations in composition and ionisation, the ejecta opacity shows substantial variations with both velocity and time. This is in part the origin of the offset between our model light curves and the predictions from the Arnett model.

Using radiation-hydrodynamics and radiative-transfer simulations, we explore the origin of the spectral diversity of interacting supernovae (SNe) of type IIn. We revisit SN1994W and investigate the dynamical configurations that can give rise to spectra with narrow lines at all times. We find that a standard ~10Msun 10^51erg SN ejecta ramming into a 0.4Msun dense CSM is inadequate for SN1994W, as it leads to the appearance of broad lines at late times. This structure, however, generates spectra that exhibit the key morphological changes seen in SN1998S. For SN1994W, we consider a completely different configuration, which involves the interaction at a large radius of a low mass inner shell with a high mass outer shell. Such a structure may arise in an 8-12Msun star from a nuclear flash (e.g., of Ne) followed within a few years by core collapse. Our simulations show that the large mass of the outer shell leads to the complete braking of the inner shell material, the formation of a slow dense shell, and the powering of a luminous SN IIn, even for a low inner shell energy. Early on, our model line profiles are typical of SNe IIn, exhibiting narrow cores and broad electron-scattering wings. As observed in SN1994W, they also remain narrow at late times. Our SN1994W model invokes two low energy ejections, both atypical of observed massive stars, and illustrates the diversity of configurations leading to SNe IIn. These results also highlight the importance of spectra to constrain the dynamical properties and understand the origin of SNe IIn.

We present 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations for supernovae (SNe) of type IIb, Ib, and Ic that result from the terminal explosion of the mass donor in a close-binary system. Here, we select three ejecta with a total kinetic energy of ~1.2e51erg, but characterised by different ejecta masses (2-5Msun), composition, and chemical mixing. The type IIb/Ib models correspond to the progenitors that have retained their He-rich shell at the time of explosion. The type Ic model arises from a progenitor that has lost its helium shell, but retains 0.32Msun of helium in a CO-rich core of 5.11Msun. We discuss their photometric and spectroscopic properties during the first 2-3 months after explosion, and connect these to their progenitor and ejecta properties including chemical stratification. For these three models, Arnett's rule overestimates the 56Ni mass by ~50% while the procedure of Katz et al., based on an energy argument, yields a more reliable estimate. The presence of strong CI lines around 9000A prior to maximum is an indicator that the pre-SN star was under-abundant in helium. As noted by others, the 1.08micron feature is a complex blend of CI, MgII, and HeI lines, which makes the identification of He uncertain in SNe Ibc unless other HeI lines can be identified. Our models show little scatter in (V-R) colour 10d after R-band maximum. We also address a number of radiative transfer properties of SNe Ibc, including the notion of a photosphere, the inference of a representative ejecta expansion rate, spectrum formation, blackbody fits and "correction factors".

We explore the properties of Type II-Plateau (II-P) supernovae (SNe) together with their red-supergiant (RSG) star progenitors. Using MESA STAR, we modulate the parameters (e.g., mixing length, overshoot, rotation, metallicity) that control the evolution of a 15Msun main-sequence star to produce a variety of physical pre-SN models and SN II-P ejecta. We extend previous modeling of SN II-P radiation to include photospheric and nebular phases, as well as multi-band light curves and spectra. Our treatment does not assume local thermodynamic equilibrium, is time dependent, treats explicitly the effects of line blanketing, and incorporates non-thermal processes. We find that the color properties of SNe II-P require large model atoms for FeI and FeII, much larger than adopted in Dessart & Hillier (2011). The color properties also imply RSG progenitors of limited extent (~500Rsun) --- larger progenitor stars produce a SN II-P radiation that remains too blue for too long. This finding calls for a reduction of RSG radii, perhaps through a strengthening of convective energy transport in RSG envelopes. Increased overshoot and rotation reduce the ratio of ejecta to helium-core mass, similarly to an increase in main-sequence mass, and thus complicate the inference ofprogenitor masses. In contrast to the great sensitivity on progenitor radius, SN II-P color evolution appears insensitive to variations in explosion energy. Finally, we document the numerous SN II-P signatures that vary with progenitor metallicity, revealing their potential for metallicity determinations in the nearby and distant Universe.

We present non-LTE time-dependent radiative-transfer simulations of pair-instability supernovae (PISNe) stemming from red-supergiant (RSG), blue-supergiant (BSG) and Wolf-Rayet (WR) star rotation-free progenitors born in the mass range 160-230Msun, at 10^-4 Zsun. Although subject to uncertainties in convection and stellar mass-loss rates, our initial conditions come from physically-consistent models that treat evolution from the main-sequence, the onset of the pair-production instability, and the explosion phase. With our set of input models characterized by large 56Ni and ejecta masses, and large kinetic energies, we recover qualitatively the Type II-Plateau, II-peculiar, and Ib/c light-curve morphologies, although they have larger peak bolometric luminosities (~10^9 to 10^10 Lsun) and a longer duration (~200d). We discuss the spectral properties for each model during the photospheric and nebular phases, including Balmer lines in II-P and II-pec at early times, the dominance of lines from intermediate-mass-elements (IMEs) near the bolometric maximum, and the strengthening of metal line blanketing thereafter. Having similar He-core properties, all models exhibit similar post-peak spectra that are strongly blanketed by FeII and FeI lines, characterized by red colors, and that arise from photospheres/ejecta with a temperature of <4000K. Combined with the modest line widths after bolometric peak, these properties contrast with those of known super-luminous SNe suggesting that PISNe are yet to be discovered. Being reddish, PISNe will be difficult to observe at high redshift except when they stem from RSG explosions, in which case they could be used as metallicity probes and distance indicators.

Much uncertainty surrounds the origin of super-luminous supernovae (SNe). Motivated by the discovery of the Type Ic SN2007bi, we study its proposed association with a pair-instability SN (PISN). We compute stellar-evolution models for primordial ~200Msun stars, simulating the implosion/explosion due to the pair-production instability, and use them as inputs for detailed non-LTE time-dependent radiative-transfer simulations that include non-local energy deposition and non-thermal processes. We retrieve the basic morphology of PISN light curves from red-supergiant, blue-supergiant, and Wolf-Rayet (WR) star progenitors. Although we confirm that a progenitor 100Msun helium core (PISN model He100) fits well the SN2007bi light curve, the low ratios of its kinetic energy and 56Ni mass to the ejecta mass, similar to standard core-collapse SNe, conspire to produce cool photospheres, red spectra subject to strong line blanketing, and narrow line profiles, all conflicting with SN2007bi observations. He-core models of increasing 56Ni-to-ejecta mass ratio have bluer spectra, but still too red to match SN2007bi, even for model He125 -- the effect of 56Ni heating is offset by the associated increase in blanketing. In contrast, the delayed injection of energy by a magnetar represents a more attractive alternative to reproduce the blue, weakly-blanketed, and broad-lined spectra of super-luminous SNe. The extra heat source is free of blanketing and is not explicitly tied to the ejecta. Experimenting with a ~9Msun WR-star progenitor, initially exploded to yield a ~1.6B SN Ib/c ejecta but later influenced by tunable magnetar-like radiation, we produce a diversity of blue spectral morphologies reminiscent of SN2007bi, the peculiar Type Ib SN2005bf, and super-luminous SN2005ap-like events.

Dark matter particles may be captured by a star and then thermalized in the star's core. At the end of its life a massive star collapses suddenly and a compact object is formed. The dark matter particles redistribute accordingly. In the inelastic dark matter model, an extended dense dark matter mini-halo surrounding the neutron star may be formed. Such mini-halos may be common in the Galaxy. The electron/positron flux resulting in the annihilation of dark matter particles, however, is unable to give rise to observable signal unless a nascent mini-halo is within a distance \sim a few 0.1 pc from the Earth.

We present non-LTE time-dependent radiative-transfer simulations of supernova (SN) IIb/Ib/Ic spectra and light curves, based on ~1B-energy piston-driven ejecta, with and without 56Ni, produced from single and binary Wolf-Rayet (W-R) stars evolved at solar and sub-solar metallicities. Our bolometric light curves show a 10-day long post-breakout plateau with a luminosity of 1-5x10^7Lsun. In our 56Ni-rich models, with ~3Msun ejecta masses, this plateau precedes a 20-30-day long re-brightening phase initiated by the outward-diffusing heat wave powered by radioactive decay at depth. In low ejecta-mass models with moderate mixing, Gamma-ray leakage starts as early as ~50d after explosion and causes the nebular luminosity to steeply decline by ~0.02mag/d. Such signatures, which are observed in standard SNe IIb/Ib/Ic, are consistent with low-mass progenitors derived from a binary-star population. We propose that the majority of stars with an initial mass ~<20Msun yield SNe II-P if 'effectively" single, SNe IIb/Ib/Ic if part of a close binary system, and SN-less black holes if more massive. Our ejecta, with outer hydrogen mass fractions as low as ~>0.01 and a total hydrogen mass of ~>0.001Msun, yield the characteristic SN IIb spectral morphology at early times. However, by ~15d after the explosion, only Halpha may remain as a weak absorption feature. Our binary models, characterised by helium surface mass fractions of ~>0.85, systematically show HeI lines during the post-breakout plateau, irrespective of the 56Ni abundance. Synthetic spectra show a strong sensitivity to metallicity, which offers the possibility to constrain it directly from SN spectroscopic modelling.

Recently several type Ib supernovae (SNe; with the prototypical SN 2005E) have been shown to have atypical properties. These SNe are faint (absolute peak magnitude of ~ -15) and fast SNe that show unique composition. They are inferred to have low ejecta mass (a few tenths of a solar mass) and to be highly enriched in calcium, but poor in silicon elements and nickel. These SNe were therefore suggested to belong to a new class of calcium-rich faint SNe explosions. Their properties were proposed to be the result of helium detonations that may occur on helium accreting white dwarfs. In this paper we theoretically study the scenario of helium detonations, and focus on the results of detonations in accreted helium layers on low mass carbon-oxygen (CO) cores. We present new results from one dimensional simulations of such explosions, including their light curves and spectra. We find that when the density of the helium layer is low enough the helium detonation produces large amounts of intermediate elements, such as calcium and titanium, together with a large amount of unburnt helium. Our results suggest that the properties of calcium-rich faint SNe could indeed be consistent with the helium-detonation scenario on small CO cores. Above a certain density (larger CO cores) the detonation leaves mainly 56Ni and unburnt helium, and the predicted spectrum will unlikely fit the unique features of this class of SNe. Finally, none of our studied models reproduces the bright, fast evolving light curves of another type of peculiar SNe suggested to originate in helium detonations (SNe 1885A, 1939B and 2002bj).

We present radiation-hydrodynamics simulations of core-collapse supernova (SN) explosions, artificially generated by driving a piston at the base of the envelope of a rotating or non-rotating red-supergiant progenitor star. We search for trends in ejecta kinematics in the resulting Type II-Plateau (II-P) SN, exploring dependencies with explosion energy and pre-SN stellar-evolution model. We recover the trivial result that larger explosion energies yield larger ejecta velocities in a given progenitor. However, we emphasise that for a given explosion energy, the increasing helium-core mass with main-sequence mass of such Type II-P SN progenitors leads to ejection of core-embedded oxygen-rich material at larger velocities. We find that the photospheric velocity at 15d after shock breakout is a good and simple indicator of the explosion energy in our selected set of pre-SN models. This measurement, combined with the width of the nebular-phase OI6303-6363A line, can be used to place an upper-limit on the progenitor main-sequence mass. Using the results from our simulations, we find that the current, but remarkably scant, late-time spectra of Type II-P SNe support progenitor main-sequence masses inferior to ~20Msun and thus, corroborate the inferences based on the direct, but difficult, progenitor identification in pre-explosion images. The narrow width of OI6303-6363A in Type II-P SNe with nebular spectra does not support high-mass progenitors in the range 25-30Msun. Combined with quantitative spectroscopic modelling, such diagnostics offer a means to constrain the main-sequence mass of the progenitor, the mass fraction of the core ejected, and thus, the mass of the compact remnant formed.

Observations of transient phenomena in the Universe reveal a spectrum of mass-ejection properties associated with massive stars, covering from Type II/Ib/Ic core-collapse supernovae (SNe) to giant eruptions of Luminous Blue Variables (LBV) and optical transients. Here, we hypothesize that a fraction of these phenomena may have an explosive origin, the distinguishing ingredient being the ratio of the prompt energy release E_dep to the envelope binding energy E_binding. Using one-dimensional one-group radiation hydrodynamics and a set of 10-25Msun, massive-star models, we explore the dynamical response of a stellar envelope subject to a strong, sudden, and deeply-rooted energy release. Following energy deposition, a shock systematically forms, crosses the progenitor envelope on a day timescale, and breaks-out with a signal of hour-to-days duration and a 10^5-10^11 Lsun luminosity. For E_dep > E_binding, full envelope ejection results with a SN-like bolometric luminosity and kinetic energy, modulations being commensurate to the energy deposited and echoing the diversity of Type II-Plateau SNe. For E_dep ~ E_binding, partial envelope ejection results with a small expansion speed, and a more modest but year-long luminosity plateau, reminiscent of LBV eruptions or so-called SN impostors. For E_dep < E_binding, we obtain a "puffed-up" star, secularly relaxing back to thermal equilibrium. In parallel with gravitational collapse and Type II SNe, we argue that the thermonuclear combustion of merely a few 0.01Msun of C/O could power a wide range of explosions/eruptions in loosely-bound massive stars, as those in the 8-12Msun range, or in more massive ones owing to their proximity to the Eddington limit and/or critical rotation.

Stars with initial masses 10 M_{solar} < M_{initial} < 100 M_{solar} fuse progressively heavier elements in their centres, up to inert iron. The core then gravitationally collapses to a neutron star or a black hole, leading to an explosion -- an iron-core-collapse supernova (SN). In contrast, extremely massive stars (M_{initial} > 140 M_{solar}), if such exist, have oxygen cores which exceed M_{core} = 50 M_{solar}. There, high temperatures are reached at relatively low densities. Conversion of energetic, pressure-supporting photons into electron-positron pairs occurs prior to oxygen ignition, and leads to a violent contraction that triggers a catastrophic nuclear explosion. Tremendous energies (>~ 10^{52} erg) are released, completely unbinding the star in a pair-instability SN (PISN), with no compact remnant. Transitional objects with 100 M_{solar} < M_{initial} < 140 M_{solar}, which end up as iron-core-collapse supernovae following violent mass ejections, perhaps due to short instances of the pair instability, may have been identified. However, genuine PISNe, perhaps common in the early Universe, have not been observed to date. Here, we present our discovery of SN 2007bi, a luminous, slowly evolving supernova located within a dwarf galaxy (~1% the size of the Milky Way). We measure the exploding core mass to be likely ~100 M_{solar}, in which case theory unambiguously predicts a PISN outcome. We show that >3 M_{solar} of radioactive 56Ni were synthesized, and that our observations are well fit by PISN models. A PISN explosion in the local Universe indicates that nearby dwarf galaxies probably host extremely massive stars, above the apparent Galactic limit, perhaps resulting from star formation processes similar to those that created the first stars in the Universe.

The discovery of the extremely luminous supernova SN 2006gy, possibly interpreted as a pair instability supernova, renewed the interest in very massive stars. We explore the evolution of these objects, which end their life as pair instability supernovae or as core collapse supernovae with relatively massive iron cores, up to about $3 M_\odot$.

The discovery of the extremely luminous supernova SN 2006gy, possibly interpreted as a pair instability supernova, renewed the interest in very massive stars. We explore the evolution of these objects, which end their life as pair instability supernovae or as core collapse supernovae with relatively massive iron cores, up to about $3 M_\odot$.

We follow the evolution of helium stars of initial mass $(2.2 - 2.5) M_\odot$, and show that they undergo off-center carbon burning, which leaves behind ${\mathbf \sim 0.01 M_\odot}$ of unburnt carbon in the inner part of the core. When the carbon-oxygen core grows to Chandrasekhar mass, the amount of left-over carbon is sufficient to ignite thermonuclear runaway. At the moment of explosion, the star will possess an envelope of several $0.1 M_{\odot}$, consisting of He, C, and possibly some H, perhaps producing a kind of peculiar SN. Based on the results of Waldman and Barkat (2007) for accreting white dwarfs, we expect to get thermonuclear runaway at a broad range of $\rho_c \approx (1 - 6) \times 10^9 \mathrm{g cm^{-3}}$, depending on the amount of residual carbon. We verified the feasibility of this scenario by showing that in a close binary system with initial masses $(8.5 + 7.7) M_{\odot}$ and initial period of 150 day the primary produces a helium remnant of $2.3 M_{\odot}$ that evolves further like the model we considered.

Exploring stellar models which ignite carbon off-center (in the mass range of about 1.05 - 1.25 Msun, depending on the carbon mass fraction) we find that they may present an interesting SN I progenitor scenario, since whereas in the standard scenario runaway always takes place at the same density of about 2 X 10^9 gr/cm^3, in our case, due to the small amount of carbon ignited, we get a whole range of densities from 1 X 10^9 up to 6 X 10^9 gr/cm^3. These results could contribute in resolving the emerging recognition that at least some diversity among SNe I exists, since runaway at various central densities is expected to yield various outcomes in terms of the velocities and composition of the ejecta, which should be modeled and compared to observations.

We explore the hypothesis, that helium stars in a certain mass range can evolve to a carbon core explosion similar to what is widely accepted as an explanation for the SN I phenomenon. This should happen when their carbon-oxygen core grows thanks to the helium shell burning above the core. We found that in the mass range of about 1.7-2.2 Msun, indeed this can happen. The main new insight we believe we gained is the crucial importance of an "early" off-center ignition of carbon, which at a later stage prevents the carbon which forms below the helium burning shell and ignites, from burning the carbon all the way to the center. When helium is almost depleted in the convective envelope by the helium burning shell at its bottom, the now super-Chandrasekhar mass carbon-oxygen core contracts, and the residual degenerate carbon at the center is ignited, resulting in a runaway similar to the classical SN I scenario. Since the structure and behavior of the carbon-oxygen core of the helium stars of our interest is very similar to that of a mass accreting carbon-oxygen star, we also thoroughly examined the behavior of carbon-oxygen stars. We discovered that the models which ignite carbon off-center (in the mass range of about 1.05-1.18 Msun, depending on the carbon mass fraction) present an interesting SN I progenitor scenario of their own, since whereas in the standard scenario runaway always takes place at the same density of about 2E9 gr/cm3, in our case, due to the small amount of carbon ignited, we get a whole range of densities from 1E9 up to 6E9 gr/cm3.