We report the discovery of a nearby massive white dwarf with He-H atmosphere. The white dwarf is located at a distance of $74.5\pm0.9$ pc. Its radius, mass, effective temperature, H/He ratio and age are $R=2500\pm100$ km, $M=1.33\pm0.01$ $\rm M_{\odot}$, $T_{\mathrm{eff}}=31200\pm 1200$ K, ${\mathrm{H/He}}\sim0.1$ and $330\pm40$ Myr, respectively. The observed spectrum is redshifted by $V_{\mathrm{r}}=+240\pm15$ km s$^{-1},$ which is mostly attributed to the gravitational redshift. The white dwarf shows a regular stable photometric variability with amplitude $\Delta g\approx0.06^{\mathrm{m}}$ and period $P=353.456$ s suggesting rapid rotation. This massive, hot and rapidly rotating white dwarf is likely to originate from the merging of close binary white dwarf system that avoided explosion in a thermonuclear type Ia supernova at the Carboniferous Period of the Earth history.

We show that the luminous supernovae associated with ultra-long gamma-ray bursts can be caused by the slow adiabatic cooling of hydrogen-free progenitors extended by pulsational pair-instability. In the accompanying paper (Marchant & Moriya 2020), we have shown that some rapidly-rotating hydrogen-free gamma-ray burst progenitors that experience pulsational pair-instability can keep an extended structure caused by pulsational pair-instability until the core collapse. Such progenitors have large radii exceeding 10 Rsun and they sometimes reach beyond 1000 Rsun at the time of the core collapse. They are, therefore, promising progenitors of ultra-long gamma-ray bursts. We here perform the light-curve modeling of the explosions of one extended hydrogen-free progenitor with a radius of 1962 Rsun. Thanks to the large progenitor radius, the ejecta experience slow adiabatic cooling and they become rapidly evolving (<~ 10 days) luminous (>~ 1e43 erg/s) supernovae even without the energy input from the 56Ni nuclear decay. They also have fast photospheric velocity above 10,000 km/s and hot photospheric temperature above 10,000 K at around the peak luminosity. We find that SN 2011kl, which is associated with the ultra-long gamma-ray burst GRB 111209A, can be explained by the adiabatic cooling phase of the extended progenitor. The ultra-long gamma-ray burst progenitors proposed in Marchant & Moriya 2020 can explain both the ultra-long gamma-ray burst duration and the accompanying supernova properties. When the gamma-ray burst jet is off-axis or choked, the luminous supernovae could be observed as fast blue optical transients without accompanying gamma-ray bursts.

Mixing of 56Ni, whose nuclear decay energy is a major luminosity source in stripped-envelope supernovae, is known to affect the observational properties of stripped-envelope supernovae such as light-curve and color evolution. Here we systematically investigate the effect of 56Ni mixing on the photospheric velocity evolution in stripped-envelope supernovae. We show that 56Ni mixing significantly affects the early photospheric velocity evolution. The photospheric velocity, which is often used to constrain the ejecta mass and explosion energy, significantly varies by just changing the degree of 56Ni mixing. In addition, the models with a small degree of 56Ni mixing show a flattening in the early photospheric velocity evolution, while the fully mixed models show a monotonic decrease. The velocity flattening appears in both helium and carbon+oxygen progenitor explosions with a variety of ejecta mass, explosion energy, and 56Ni mass. Some stripped-envelope supernovae with early photospheric velocity information do show such a flattening. We find that Type Ib SN 2007Y, which has early photospheric velocity information, has a signature of a moderate degree of 56Ni mixing in the photospheric velocity evolution and about a half of the ejecta is mixed in it. The immediate spectroscopic follow-up observations of stripped-envelope supernovae shortly after the explosion providing the early photospheric evolution give an important clue to constrain 56Ni mixing in the ejecta.

We explore the unique supernova (SN) "Refsdal" - the first discovered gravitationally lensed SN with multiple images. A large magnification provided by the galactic-scale lens, augmented by the cluster lens, gave us the opportunity to perform a detailed modelling of a distant SN at $z \simeq 1.5$. % We present results of radiation hydrodynamics modelling of \snref. According to our calculations, the \snref progenitor is likely to be a more massive and energetic version of \sna, i.e. a blue supergiant star with a radius of ${\rm R} = 50 R_{\odot}$ before an explosion and a total mass of ${\rm M}= 26 M_{\odot}$. The total energy release is ${\rm E}=5\times10^{51}$ erg. Reconstruction of SN light curves allowed us to obtain time delays and magnifications for the images S2-S4 relative to S1 with higher accuracy than previous estimates of Rodney et al. We estimate the Hubble constant $H_0 = 68.6^{+13.6}_{-9.7}$ km s$^{-1}$ Mpc$^{-1}$ via re-scaling the time delays predicted by different lens models to match the values obtained in this work. With more photometric data on the fifth image SX, we will be able to further refine the time delay and magnification estimates for SX and obtain competitive constraints on $H_0$.

The properties of low-mass neutron stars with rigid rotation are considered. The possible evolution paths of such stars in a close binary system with mass transfer are calculated. The properties of the gamma-ray burst GRB 170817A interpreted in terms of the stripping model, a natural ingredient of which is the explosion of a low-mass neutron star - a binary component, are briefly discussed.

We present nearly self-consistent stellar-halo models of Fornax dwarf spheroidal galaxy associated with the Milky Way galaxy. Such satellites are dominated by dark matter and have almost no gas in the system. Therefore, they are excellent objects for N-body modelling that takes into account visible and dark matter halo components. We perform seven dynamical models in total based on two N-body source code sets. Three of them are made by self-consistent realizations containing a spheroidal body modeled by truncated King or Plummer models for stellar components and the lowered Evans dark halo model. These components are given by the distribution functions which are the functions of energy $E$ and angular momentum $L_z$ only. The others are made by the AGAMA code based on the calculation of orbits. To guess the parameters which make the N-body models close to the visible object we use hydrodynamic stellar-dark model of the Fornax galaxy taking into account the velocity anisotropy parameter. The AGAMA models show the best agreement of the resulting velocity dispersion profiles with the observed data.

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.

We calculate the evolution of massive stars, which undergo pulsational pair-instability (PPI) when the O-rich core is formed. The evolution from the main-sequence through the onset of PPI is calculated for stars with the initial masses of $80 - 140$ $M_{\odot}$ and metallicities of $Z = 10^{-3} - 1.0$ $Z_\odot$. Because of mass loss, $Z \leq 0.5$ $Z_\odot$ is necessary for stars to form He cores massive enough (i.e., mass $>40 ~M_\odot$) to undergo PPI. The hydrodynamical phase of evolution from PPI through the beginning of Fe core collapse is calculated for the He cores with masses of $40 - 62 ~M_\odot$ and $Z = 0$. During PPI, electron-positron pair production causes a rapid contraction of the O-rich core which triggers explosive O-burning and a pulsation of the core. We study the mass dependence of the pulsation dynamics, thermodynamics, and nucleosynthesis. The pulsations are stronger for more massive He cores and result in such a large amount of mass ejection such as $3 - 13$ $M_\odot$ for $40 - 62 ~M_\odot$ He cores. These He cores eventually undergo Fe-core collapse. The $64 ~M_\odot$ He core undergoes complete disruption and becomes a pair-instability supernova. The H-free circumstellar matter ejected around these He cores is massive enough for to explain the observed light curve of Type I (H-free) superluminous supernovae with circumstellar interaction. We also note that the mass ejection sets the maximum mass of black holes (BHs) to be $\sim 50$ $M_{\odot}$, which is consistent with the masses of BHs recently detected by VIRGO and aLIGO.

We perform a high-cadence transient survey with Subaru Hyper Suprime-Cam (HSC), which we call the Subaru HSC survey Optimized for Optical Transients (SHOOT). We conduct HSC imaging observations with time intervals of about one hour on two successive nights, and spectroscopic and photometric follow-up observations. A rapidly declining blue transient SHOOT14di at $z=0.4229$ is found in observations on two successive nights with an image subtraction technique. The rate of brightness change is $+1.28^{+0.40}_{-0.27}~{\rm mag~day^{-1}}$ ($+1.83^{+0.57}_{-0.39}~{\rm mag~day^{-1}}$) in the observer (rest) frame and the rest-frame color between $3400$ and $4400~\unicode[.8,0]{x212B}$ is $M_{\rm 3400\unicode[.8,0]{x212B}}-M_{\rm 4400\unicode[.8,0]{x212B}}=-0.4$. The nature of the object is investigated by comparing its peak luminosity, decline rate, and color with those of transients and variables previously observed, and those of theoretical models. None of the transients or variables share the same properties as SHOOT14di. Comparisons with theoretical models demonstrate that, while the emission from the cooling envelope of a Type IIb supernova shows a slower decline rate than SHOOT14di, and the explosion of a red supergiant star with a dense circumstellar wind shows a redder color than SHOOT14di, the shock breakout at the stellar surface of the explosion of a $25M_{\odot}$ red supergiant star with a small explosion energy of $\leq0.4\times10^{51}$ erg reproduces the multicolor light curve of SHOOT14di. This discovery shows that a high-cadence, multicolor optical transient survey at intervals of about one hour, and continuous and immediate follow-up observations, is important for studies of normal core-collapse supernovae at high redshifts.

High cadence ultraviolet, optical and near-infrared photometric and low-resolution spectroscopic observations of the peculiar Type II supernova (SN) 2018hna are presented. The early phase multiband light curves exhibit the adiabatic cooling envelope emission following the shock breakout up to ~14 days from the explosion. SN~2018hna has a rise time of $\sim$\,88 days in the V-band, similar to SN 1987A. A $\rm^{56}Ni$ mass of ~0.087$\pm$0.004 $\rm M_{\odot}$ is inferred for SN 2018hna from its bolometric light curve. Hydrodynamical modelling of the cooling phase suggests a progenitor with a radius ~50 $\rm R_{\odot}$, a mass of ~14-20 $\rm M_{\odot}$ and explosion energy of ~1.7-2.9$\rm \times$ $\rm 10^{51}\ erg$. The smaller inferred radius of the progenitor than a standard red supergiant is indicative of a blue supergiant progenitor of SN 2018hna. A sub-solar metallicity (~0.3 $\rm Z_{\odot}$) is inferred for the host galaxy UGC 07534, concurrent with the low-metallicity environments of 1987A-like events.

We investigate a simplified kinetic system of the hydrogen atom (two levels plus continuum) under conditions of a type IIP supernova at the plateau phase that realistically describes the basic properties of the complete system. We have found the Lyapunov function for the reduced system using which we have analytically obtained the ionization freeze-out effect on long time scales. Since the system completely recombines in the equilibrium approximation on long time scales, which does not occur in reality, this result confirms the necessity of allowance for the time-dependent effect in the kinetics during the photospheric phase in a supernova explosion.

Recent supernova searches revealed a number of fast-evolving luminous transients. We perform radiation hydrodynamic simulations of light curves of several models of supernova from super--asymptotic giant branch (super-AGB, SAGB) stars with low mass envelopes ($M_{\rm env}$ $\sim$ $0.05-1.25$ $M_{\odot}$). The differences in the light curves found among the models are used to link the observed events to the most appropriate models. In particular, we propose that KSN 2015K is an electron-capture supernova. We assume "optically" thick CSM around SAGB and the circumstellar interaction powers the peak luminosity of the light curve with a short rise time. The faint tail might be influenced by the spin-down luminosity of a newborn Crab-like pulsar. Our fits indicate an ejecta mass of $0.02-0.05$ $M_{\odot}$, a circumstellar medium (CSM) mass of $0.10-0.12$ $M_{\odot}$, a radius of the CSM photosphere $\sim$ 10$^{14}$ cm, a kinetic energy of $\sim$ 3 $\times$ 10$^{50}$ erg, a photospheric velocity $v_{\rm ph} \gtrsim 10,000$ km s$^{-1}$ and a pulsar total spin energy ($2.5-4$) $\times$ 10$^{49}$ erg.

We present synthetic light curves of fallback-powered supernovae based on a neutrino-driven explosion of a 40 Msun zero-metallicity star with significant fallback accretion onto a black hole that was previously simulated by Chan et al. (2018) until shock breakout. Here, we investigate the light curve properties of the explosion after shock breakout for various fallback models. Without extra power from fallback accretion, the light curve is that of a Type IIP supernova with a plateau magnitude of around -14 mag and a plateau duration of 40 days. With extra power for the light curve from fallback accretion, however, we find that the transient can be significantly more luminous. The light-curve shape can be SN 1987A-like or Type IIP-like, depending on the efficiency of the fallback engine. If the accretion disk forms soon after the collapse and more than 1% of the accretion energy is released as the central engine, fallback accretion powered supernovae become as luminous as superluminous supernovae. We suggest that Type II superluminous supernovae with broad hydrogen features could be related to such hydrogen-rich supernovae powered by fallback accretion. In the future, such superluminous supernovae powered by fallback accretion might be found among the supernovae from the first stars in addition to pair-instability supernovae and pulstational pair-instability supernovae.

We present optical and near-infrared observations of the rapidly evolving supernova (SN) 2017czd that shows hydrogen features. The optical light curves exhibit a short plateau phase ($\sim 13$ days in the $R$-band) followed by a rapid decline by $4.5$ mag in $\sim 20 \mathrm{days}$ after the plateau. The decline rate is larger than those of any standard SNe, and close to those of rapidly evolving transients. The peak absolute magnitude is $-16.8$ mag in the $V$-band, which is within the observed range for SNe IIP and rapidly evolving transients. The spectra of SN 2017czd clearly show the hydrogen features and resemble those of SNe IIP at first. The H$\alpha$ line, however, does not evolve much with time and it becomes similar to those in SNe IIb at decline phase. We calculate the synthetic light curves using a SN IIb progenitor which has 16 M$_{\odot}$ at the zero-age main sequence and evolves in a binary system. The model with a low explosion energy ($5\times 10^{50}$ erg) and a low ${}^{56}$Ni mass ($0.003 \mathrm{M}_{\odot}$) can reproduce the short plateau phase as well as the sudden drop of the light curve as observed in SN 2017czd. We conclude that SN 2017czd might be the first identified weak explosion from a SN IIb progenitor. We suggest that some rapidly evolving transients can be explained by such a weak explosion of the progenitors with little hydrogen-rich envelope.

We draw attention to recent high explosive (HE) experiments which provide compression of macroscopic amount of matter to high, even record, values of pressure in comparison with other HE experiments. The observed bounce after the compression corresponds to processes in core-collapse supernova explosions after neutrino trapping. Conditions provided in the experiments resemble those in core-collapse supernovae, permitting their use for laboratory astrophysics. A unique feature of the experiments is compression at low entropy. The values of specific entropy are close to those obtained in numerical simulations during the process of collapse in supernova explosions, and much lower than those obtained at laser ignition facilities, another type of high-compression experiment. Both in supernovae and HE experiments the bounce occurs at low entropy, so the HE experiments provide a new platform to realize some supernova collapse effects in laboratory, especially to study hydrodynamics of collapsing flows and the bounce. Due to the good resolution of diagnostics in the compression of macroscopic amounts of material with essential effects of nonideal plasma in EOS, and observed development of 3D instabilities, these experiments may serve as a useful benchmark for astrophysical hydrodynamic codes.

We perform hydrodynamical simulations of the interaction between supernova (SN) ejecta and circumstellar medium (CSM) for SN 1993J and SN 2011dh, and calculate the radio and X-ray emissions expected from the shocked gas at late epochs ($t$). Considering the ejecta structure from multi-group radiation hydrodynamics simulation, we find that the observed rapid drop in radio and X-ray light curves of SN 1993J at $t>$3000 days can be due to a change in the mass-loss rate ($\dot M$) around $\sim$6500 years prior to the explosion of the SN. The exact epoch scales inversely with the assumed wind velocity of $v_{\rm w}=10~ km~s^{-1}$. The progenitor of this SN very likely belonged to a binary system, where, during its evolution, the primary had transferred material to the secondary. It is argued in the paper that the change in $\dot M$ can happen because of a change in the mass accretion efficiency ($\eta$) of the companion star. It is possible that before $\sim6500~(v_{\rm w}/10~km~s^{-1})^{-1}$years prior to the explosion, $\eta$ was high, thus the CSM was tenuous, which causes the late time downturn in fluxes. In the case of SN 2011dh, the late time evolution is found to be consistent with a wind medium with $\dot M/v_{\rm w}=4\times10^{-6 }~M_{\odot}~ yr^{-1}/10 ~{km ~s^{-1}}$. It is difficult from our analysis to predict whether the progenitor of this SN had a binary companion, however, if future observations show similar decrease in radio and X-ray fluxes, then that would give strong support to a scenario where both SNe had undergone similar kind of binary evolution before explosion.

We explore the appearance of light clusters at high densities of collapsing stellar cores. Special attention is paid to the unstable isotope H4, which was not included in previous studies. The importance of light clusters in the calculation of rates for neutrino matter interaction is discussed. The main conclusion is that thermodynamic quantities are only weakly sensitive to the chemical composition. The change in pressure and hence the direct change in collapse dynamics will be minor. But the change in neutrino heating and neutronization processes can be significant.

We present optical and ultraviolet observations of nearby type Ic supernova SN 2017ein as well as detailed analysis of its progenitor properties from both the early-time observations and the prediscovery Hubble Space Telescope (HST) images. The optical light curves started from within one day to $\sim$275 days after explosion, and optical spectra range from $\sim$2 days to $\sim$90 days after explosion. Compared to other normal SNe Ic like SN 2007gr and SN 2013ge, \mbox{SN 2017ein} seems to have more prominent C{\footnotesize II} absorption and higher expansion velocities in early phases, suggestive of relatively lower ejecta mass. The earliest photometry obtained for \mbox{SN 2017ein} show indications of shock cooling. The best-fit obtained by including a shock cooling component gives an estimate of the envelope mass as $\sim$0.02 M$_{\odot}$ and stellar radius as 8$\pm$4 R$_{\odot}$. Examining the pre-explosion images taken with the HST WFPC2, we find that the SN position coincides with a luminous and blue point-like source, with an extinction-corrected absolute magnitude of M$_V$$\sim$$-$8.2 mag and M$_I$$\sim$$-$7.7 mag.Comparisons of the observations to the theoretical models indicate that the counterpart source was either a single WR star or a binary with whose members had high initial masses, or a young compact star cluster. To further distinguish between different scenarios requires revisiting the site of the progenitor with HST after the SN fades away.

MASTER OT J033744.97+723159.0/SN 2017gpn was discovered in the error-box of LIGO/Virgo G299232 gravitational wave alert. The spectrum of SN 2017gpn is consistent with a Type IIb supernova. In this work we present the photometry of 20 epochs of observations performed with CCD-photometer on the Zeiss-1000 telescope. The light curves in B and Rc filters were obtained. The multicolour light curves were also modelled numerically using the one-dimensional radiation hydro-dynamical code STELLA.

Blue-supergiant stars develop into core-collapse supernovae --- one of the most energetic outbursts in the universe --- when all nuclear burning fuel is exhausted in the stellar core. Previous attempts failed to explain observed explosions of such stars which have a zero-age main sequence mass of 50~M$_\odot$ or more. Here we exploit the largely uncertain state of matter at high density, and connect the modeling of such stellar explosions with a first-order phase transition from nuclear matter to the quark-gluon plasma. The resulting energetic supernova explosions can account for a large variety of lightcurves, from peculiar type II to super-luminous events. The remnants are neutron stars with quark matter core, known as hybrid stars, of about 2~M$_\odot$ at birth. A galactic event of this kind could be observable due to the release of a second neutrino burst. Its observation would confirm such a first-order phase transition at densities relevant for astrophysics.

We investigate the effect of mixing of radioactive nickel ($^{56}$Ni) on the early-time color evolution of Type Ib and Ic supernovae (SNe Ib/Ic) using multi-group radiation hydrodynamics simulations. We consider both helium-rich and helium-poor progenitors. Mixing of $^{56}$Ni is parameterized using a Gaussian distribution function. We find that the early-time color evolution with a weak $^{56}$Ni mixing is characterized by three different phases: initial rapid reddening, blueward evolution due to the delayed effect of $^{56}$Ni heating, and redward evolution thereafter until the transition to the nebular phase. With a strong $^{56}$Ni mixing, the second phase disappears. We compare our models with the early-time color evolution of several SNe Ib/Ic (SN1999ex, SN 2008D, SN 2009jf, iPTF13bvn, SN 1994I, SN 2007gr, SN 2013ge, and 2017ein) and find signatures of relatively weak and strong $^{56}$Ni mixing for SNe Ib and SNe Ic, respectively. This suggests that SNe Ib progenitors are distinct from SN Ic progenitors in terms of helium content and that $^{56}$Ni mixing is generally stronger in the carbon-oxygen core and weaker in the helium-rich envelope. We conclude that the early-time color evolution is a powerful probe of $^{56}$Ni mixing in SNe Ib/Ic.

Type II supernovae (SNe) originate from the explosion of hydrogen-rich supergiant massive stars. Their first electromagnetic signature is the shock breakout, a short-lived phenomenon which can last from hours to days depending on the density at shock emergence. We present 26 rising optical light curves of SN II candidates discovered shortly after explosion by the High cadence Transient Survey (HiTS) and derive physical parameters based on hydrodynamical models using a Bayesian approach. We observe a steep rise of a few days in 24 out of 26 SN II candidates, indicating the systematic detection of shock breakouts in a dense circumstellar matter consistent with a mass loss rate $\dot{M} > 10^{-4} M_\odot yr^{-1}$ or a dense atmosphere. This implies that the characteristic hour timescale signature of stellar envelope SBOs may be rare in nature and could be delayed into longer-lived circumstellar material shock breakouts in most Type II SNe.

The discovery of GW signal from merging neutron stars by LIGO on 17th August 2017 was followed by a short GRB170817A discovered by FERMI and INTEGRAL 1.7 seconds after the loss of the GW signal when it just reached its maximum. Here we present a reproduction of the first paper (published by us in 1984) predicting a short GRB after GW signal of merging neutron stars. Our paper followed the scenario by Clark and Eardley (1977) who predicted a catastrophic disruption of a neutron star in a binary 1.7 seconds after the peak of GW signal. Our next paper in 1990 predicted all the main properties of the short GRB with quite a reasonable accuracy. Typos in English translation are corrected and a few comments are added in the current publication as numbered footnotes (the only footnote from the original paper is marked by an asterisk).

We present an approach to finding distances to young supernova remnants. Our method is based on hydrodynamical simulations of Tycho's SNR using the SUPREMNA code. For the explosion models, we use the classical W7 deflagration model and the delayed detonation model that was previously shown to provide good fits to the X-ray emission of Tycho's SNR. Combining our hydrodynamical simulations with the Chandra high-resolution images of the remnant obtained in 2015, we derive a distance estimate to Tycho's SNR of $2.8\pm 0.4$ kpc.

We introduce the first synthetic light-curve model set of Type IIP supernovae exploded within circumstellar media in which the acceleration of the red supergiant winds is taken into account. Because wind acceleration makes the wind velocities near the progenitors low, the density of the immediate vicinity of the red supergiant supernova progenitors can be higher than that extrapolated by using a constant terminal wind velocity. Therefore, even if the mass-loss rate of the progenitor is relatively low, it can have a dense circumstellar medium at the immediate stellar vicinity and the early light curves of Type IIP supernovae are significantly affected by it. We adopt a simple beta velocity law to formulate the wind acceleration. We provide bolometric and multicolor light curves of Type IIP supernovae exploding within such accelerated winds from the combinations of three progenitors, 12 - 16 Msun; five beta, 1-5; seven mass-loss rates, 1e-5 - 1e-2 Msun/yr; and four explosion energies, 0.5e51 - 2e51 erg. All the light curve models are available at https://goo.gl/o5phYb. When the circumstellar density is sufficiently high, our models do not show a classical shock breakout as a consequence of the interaction with the dense and optically-thick circumstellar media. Instead, they show a delayed 'wind breakout', substantially affecting early light curves of Type IIP supernovae. We find that the mass-loss rates of the progenitors need to be 1e-3 - 1e-2 Msun/yr to explain typical rise times of 5 - 10 days in Type IIP supernovae assuming a dense circumstellar radius of 1e15 cm.

We update the capabilities of the software instrument Modules for Experiments in Stellar Astrophysics (MESA) and enhance its ease of use and availability. Our new approach to locating convective boundaries is consistent with the physics of convection, and yields reliable values of the convective core mass during both hydrogen and helium burning phases. Stars with $M<8\,{\rm M_\odot}$ become white dwarfs and cool to the point where the electrons are degenerate and the ions are strongly coupled, a realm now available to study with MESA due to improved treatments of element diffusion, latent heat release, and blending of equations of state. Studies of the final fates of massive stars are extended in MESA by our addition of an approximate Riemann solver that captures shocks and conserves energy to high accuracy during dynamic epochs. We also introduce a 1D capability for modeling the effects of Rayleigh-Taylor instabilities that, in combination with the coupling to a public version of the STELLA radiation transfer instrument, creates new avenues for exploring Type II supernovae properties. These capabilities are exhibited with exploratory models of pair-instability supernova, pulsational pair-instability supernova, and the formation of stellar mass black holes. The applicability of MESA is now widened by the capability of importing multi-dimensional hydrodynamic models into MESA. We close by introducing software modules for handling floating point exceptions and stellar model optimization, and four new software tools -- MESAWeb, MESA-Docker, pyMESA, and mesastar.org -- to enhance MESA's education and research impact.

We study the evolution and final outcome of long-lived (${\approx}10^5$ years) remnants from the merger of a He white dwarf (WD) with a more massive C/O or O/Ne WD. Using Modules for Experiments in Stellar Astrophysics ($\texttt{MESA}$), we show that these remnants have a red giant configuration supported by steady helium burning, adding mass to the WD core until it reaches $M_{\rm core}\approx 1.12-1.20 M_\odot$. At that point, the base of the surface convection zone extends into the burning layer, mixing the helium burning products (primarily carbon and magnesium) throughout the convective envelope. Further evolution depletes the convective envelope of helium, and dramatically slows the mass increase of the underlying WD core. The WD core mass growth re-initiates after helium depletion, as then an uncoupled carbon burning shell is ignited and proceeds to burn the fuel from the remaining metal-rich extended envelope. For large enough initial total merger masses, O/Ne WD cores would experience electron-capture triggered collapse to neutron stars (NSs) after growing to near Chandrasekhar mass ($M_{\rm Ch}$). Massive C/O WD cores could suffer the same fate after a carbon-burning flame converts them to O/Ne. The NS formation would release ${\approx}10^{50}$ ergs into the remaining extended low mass envelope. Using the STELLA radiative transfer code, we predict the resulting optical light curves from these exploded envelopes. Reaching absolute magnitudes of $M_V\approx -17$, these transients are bright for about one week, and have many features of the class of luminous, rapidly evolving transients studied by Drout and collaborators.

We investigate the possibility that the energetic Type II supernova OGLE-2014-SN-073 is powered by a fallback accretion following the failed explosion of a massive star. Taking massive hydrogen-rich supernova progenitor models, we estimate the fallback accretion rate and calculate the light curve evolution of supernovae powered by the fallback accretion. We find that such fallback accretion powered models can reproduce the overall observational properties of OGLE-2014-SN-073. It may imply that some failed supernovae could be observed as energetic supernovae like OGLE-2014-SN-073 instead of faint supernovae as previously proposed.

HD 49798/RX J0648.0--4418 is a peculiar X-ray binary with a hot subdwarf (sdO) mass donor. The nature of the accreting compact object is not known, but its spin period $P=13.2$~s and $\dot P =-2.15 \times 10^{-15}$s~s$^{-1}$, prove that it can be only either a white dwarf or a neutron star. The spin-up has been very stable for more than 20 years. We demonstrate that the continuous stable spin-up of the compact companion of HD 49798 can be best explained by contraction of a young white dwarf with an age $\sim 2$~Myrs. This allows us to interpret all the basic parameters of the system in the framework of an accreting white dwarf. We present examples of binary evolution which result in such systems. If correct, this is the first direct evidence for a white dwarf contraction on early evolutionary stages.

Upcoming high-cadence transient survey programmes will produce a wealth of observational data for Type Ia supernovae. These data sets will contain numerous events detected very early in their evolution, shortly after explosion. Here, we present synthetic light curves, calculated with the radiation hydrodynamical approach Stella for a number of different explosion models, specifically focusing on these first few days after explosion. We show that overall the early light curve evolution is similar for most of the investigated models. Characteristic imprints are induced by radioactive material located close to the surface. However, these are very similar to the signatures expected from ejecta-CSM or ejecta-companion interaction. Apart from the pure deflagration explosion models, none of our synthetic light curves exhibit the commonly assumed power-law rise. We demonstrate that this can lead to substantial errors in the determination of the time of explosion. In summary, we illustrate with our calculations that even with very early data an identification of specific explosion scenarios is challenging, if only photometric observations are available.

Observations of Gaia16apd revealed extremely luminous ultraviolet emission among superluminous supernovae (SLSNe). Using radiation hydrodynamics simulations we perform a comparison of UV light curves, color temperatures and photospheric velocities between the most popular SLSN models: pair-instability supernova, magnetar and interaction with circumstellar medium. We find that the interaction model is the most promising to explain the extreme UV luminosity of Gaia16apd. The differences in late-time UV emission and in color evolution found between the models can be used to link an observed SLSN event to the most appropriate model. Observations at UV wavelengths can be used to clarify the nature of SLSNe and more attention should be paid to them in future follow-up observations.

The importance of allowance for the time-dependent effect in the kinetics at the photospheric phase during a supernova explosion has been confirmed by several independent research groups. The time-dependent effect provides a higher degree of hydrogen ionization in comparison with the steady-state solutions and strengthens the H$\alpha$ line in the resulting simulated spectrum, with the intensity of the effect increasing with time. However, some researchers argue that the time-dependent ionization effect is unimportant. Its allowance leads to an insignificant strengthening of H$\alpha$ in their modeling only in the first days after explosion. We have demonstrated the importance of the time-dependent effect with the models of SN 1999em as an example using the new original LEVELS software package. The role of a number of factors that can weaken the time-dependent effect has been checked. We have confirmed that the intensity of the effect is affected by the abundance of metal admixtures in the envelope, while the addition of extra levels to themodel hydrogen atom weakens the time-dependent effect to a lesser degree and never removes it completely.

We report the discovery and multicolor (VRIW) photometry of a rare explosive star MASTER OT J004207.99+405501.1 - a luminous red nova - in the Andromeda galaxy M31N2015-01a. We use our original light curve acquired with identical MASTER Global Robotic Net telescopes in one photometric system: VRI during first 30 days and W (unfiltered) during 70 days. Also we added publishied multicolor photometry data to estimate the mass and energy of the ejected shell, and discuss the likely formation scenarios of outbursts of this type. We propose the interpretation of the explosion, that is consistent with the evolutionary scenario where star merger is a natural stage of the evolution of close-mass stars and may serve as an extra channel for the formation of nova outbursts.

Type IIP supernova progenitors are often surrounded by dense circumstellar media that may result from mass-loss enhancement of the progenitors shortly before their explosions. Previous light-curve studies suggest that the mass-loss rates are enhanced up to ~ 0.1 Msun/yr, assuming a constant wind velocity. However, density of circumstellar media at the immediate stellar vicinity can be much higher than previously inferred for a given mass-loss rate if wind acceleration is taken into account. We show that the wind acceleration has a huge impact when we estimate mass-loss rates from early light curves of Type IIP supernovae by taking SN 2013fs as an example. We perform numerical calculations of the interaction between supernova ejecta and circumstellar media with a constant mass-loss rate but with a beta-law wind velocity profile. We find that the mass-loss rate of the progenitor of SN 2013fs shortly before the explosion, which was inferred to be ~ 0.1 Msun/yr with a constant wind velocity of 10 km/s by a previous light-curve modeling, can be as low as ~ 1e-3 Msun/yr with the same terminal wind velocity of 10 km/s but with a wind velocity profile with beta ~ 5. In both cases, the mass of the circumstellar medium is similar (~ 0.5 Msun). Therefore, the beginning of the progenitor's mass-loss enhancement in our interpretation is ~ 100 years before the explosion, not several years. Our result indicates that the immediate dense environment of Type II supernova progenitors may be significantly influenced by wind acceleration.

Being a superluminous supernova (SLSN), PTF12dam can be explained by a $^{56}$Ni-powered model, a magnetar-powered model or an interaction model. We propose that PTF12dam is a pulsational pair instability supernova, where the outer envelope of a progenitor is ejected during the pulsations. Thus, it is powered by double energy source: radioactive decay of $^{56}$Ni and a radiative shock in a dense circumstellar medium. To describe multicolor light curves and spectra we use radiation hydrodynamics calculations of STELLA code. We found that light curves are well described in the model with 40M$_{\odot}$ ejecta and 20-40M$_{\odot}$ circumstellar medium. The ejected $^{56}$Ni mass is about 6M$_{\odot}$ which results from explosive nucleosynthesis with large explosion energy (2-3)$\cdot$10$^{52}$ ergs. In comparison with alternative scenarios of pair-instability supernova and magnetar-powered supernova, in interaction model all the observed main photometric characteristics are well reproduced: multicolor light curves, color temperatures, and photospheric velocities.

We present optical observations of SN 2013dx, related to the Fermi burst GRB 130702A occurred at a redshift z = 0.145. It is the second-best sampled GRB-SN after SN~1998bw: the observational light curves contain more than 280 data points in uBgrRiz filters until 88 day after the burst, and the data were collected from our observational collaboration (Maidanak Observatory, Abastumani Observatory, Crimean Astrophysical Observatory, Mondy Observatory, National Observatory of Turkey, Observatorio del Roque de los Muchachos) and from the literature. We model numerically the multicolour light curves using the one-dimensional radiation hydrodynamical code STELLA, previously widely implemented for the modelling of typical non-GRB SNe. The best-fitted model has the following parameters: pre-supernova star mass M = 25 M_Sun, mass of a compact remnant M_CR = 6 M_Sun, total energy of the outburst E_oburst = 3.5 x 10^(52) erg, pre-supernova star radius R = 100 R_Sun, M_56Ni = 0.2 M_Sun which is totally mixed through the ejecta; M_O = 16.6 M_Sun, M_Si = 1.2 M_Sun, and M_Fe = 1.2 M_Sun, and the radiative efficiency of the SN is 0.1 per cent.

We investigate light-curve and spectral properties of ultra-stripped core-collapse supernovae. Ultra-stripped supernovae are the explosions of heavily stripped massive stars which lost their envelopes via binary interactions with a compact companion star. They eject only ~ 0.1 Msun and may be the main way to form double neutron-star systems which eventually merge emitting strong gravitational waves. We follow the evolution of an ultra-stripped supernova progenitor until iron core collapse and perform explosive nucleosynthesis calculations. We then synthesize light curves and spectra of ultra-stripped supernovae using the nucleosynthesis results and present their expected properties. Ultra-stripped supernovae synthesize ~ 0.01 Msun of radioactive 56Ni, and their typical peak luminosity is around 1e42 erg/s or -16 mag. Their typical rise time is 5 - 10 days. Comparing synthesized and observed spectra, we find that SN 2005ek, some of the so-called calcium-rich gap transients, and SN 2010X may be related to ultra-stripped supernovae. If these supernovae are actually ultra-stripped supernovae, their event rate is expected to be about 1 per cent of core-collapse supernovae. Comparing the double neutron-star merger rate obtained by future gravitational-wave observations and the ultra-stripped supernova rate obtained by optical transient surveys identified with our synthesized light-curve and spectral models, we will be able to judge whether ultra-stripped supernovae are actually a major contributor to the binary neutron star population and provide constraints on binary stellar evolution.

The black hole binary properties inferred from the LIGO gravitational wave signal GW150914 posed several serious problems. The high masses and low effective spin of black hole binary can be explained if they are primordial (PBH) rather than the products of the stellar binary evolution. Such PBH properties are postulated ad hoc but not derived from fundamental theory. We show that the necessary features of PBHs naturally follow from the slightly modified Affleck-Dine (AD) mechanism of baryogenesis. The log-normal distribution of PBHs, predicted within the AD paradigm, is adjusted to provide an abundant population of low-spin stellar mass black holes. The same distribution gives a sufficient number of quickly growing seeds of supermassive black holes observed at high redshifts and may comprise an appreciable fraction of Dark Matter which does not contradict any existing observational limits. Testable predictions of this scenario are discussed.

A brief pedagogical review for students on basic physics of Superluminous Supernovae. Almost no references to literature are given, but important formulas are derived in some detail.

We investigate the consequences of fairly normal Type Ia supernovae being embedded in compact and dense envelopes of carbon and oxygen rich circumstellar material by means of detailed radiation hydrodynamic simulations. Our main focus rests on exploring the effects of the interaction between ejecta and circumstellar material on the ejecta evolution and the broad-band light curve. In our calculations, we find that a strong reverse shock efficiently decelerates and compresses the ejecta material. This leads to a significant broadening of the optical light curve, a longer rise to maximum and a slower decline in the tail phase. During the interaction, substantial radiative energy is generated, which mostly emerges in the extreme ultraviolet and X-ray regime. Only if reprocessing due to radiation-matter interactions is very efficient, a significant boost in the optical light curve is observed. We discuss these findings in particular in the context of the super-luminous event SN 2009dc. As our calculations are able to reproduce a number of its peculiar properties, we conclude that the flavour of the interaction scenario investigated in this work constitutes a promising candidate to explain such 'Super-Chandrasekhar' supernovae.

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.

Rapidly rotating, strongly magnetized neutron stars (magnetars) can release their enormous rotational energy via magnetic spin-down, providing a power source for bright transients such as superluminous supernovae. On the other hand, particularly massive (so-called supramassive) neutron stars require a minimum rotation rate to support their mass against gravitational collapse, below which the neutron star collapses to a black hole. We model the light curves of supernovae powered by magnetars which transform into black holes. Although the peak luminosities can reach high values in the range of superluminous supernovae, their post maximum light curves can decline very rapidly because of the sudden loss of the central energy input. Early black hole transformation also enhances the shock breakout signal from the magnetar-driven bubble relative to the main supernova peak. Our synthetic light curves of supernovae powered by magnetars transforming to black holes are consistent with those of some rapidly evolving bright transients recently reported by Arcavi et al. (2016).

A number of Type I (hydrogenless) superluminous supernova (SLSN) events have been discovered recently. However, their nature remains debatable. One of the most promising ideas is the shock-interaction mechanism, but only simplified semi-analytical models have been applied so far. We simulate light curves for several Type I SLSN (SLSN-I) models enshrouded by dense, non-hydrogen circumstellar envelopes, using a multi-group radiation hydrodynamics code that predicts not only bolometric, but also multicolor light curves. We demonstrate that the bulk of SLSNe-I including those with relatively narrow light curves like SN 2010gx or broad ones like PTF09cnd can be explained by the interaction of the SN ejecta with he CS envelope, though the range of parameters for these models is rather wide. Moderate explosion energy ($\sim (2 - 4)\cdot 10^{51}$ ergs) is sufficient to explain both narrow and broad SLSN-I light curves, but ejected mass and envelope mass differ for those two cases. Only 5 to 10 $M_\odot$ of non-hydrogen material is needed to reproduce the light curve of SN 2010gx, while the best model for PTF09cnd is very massive: it contains almost $ 50 M_\odot $ in the CS envelope and only $ 5 M_\odot $ in the ejecta. The CS envelope for each case extends from 10 $R_\odot$ to $\sim 10^5R_\odot$ ($7\cdot 10^{15} $ cm), which is about an order of magnitude larger than typical photospheric radii of standard SNe near the maximum light. We briefly discuss possible ways to form such unusual envelopes.

We describe our modelling of the radiatively cooling shocks and their thin shells with various numerical tools in different physical and calculational setups. We inspect structure of the dense shell, its formation and evolution, pointing out physical and numerical factors that sustain its shape and also may lead to instabilities. We have found that under certain physical conditions, the circular shaped shells show a strong bending instability and successive fragmentation on Cartesian grids soon after their formation, while remain almost unperturbed when simulated on polar meshes. We explain this by physical Rayleigh--Taylor like instabilities triggered by corrugation of the dense shell surfaces by numerical noise. Conditions for these instabilities follow from both the shell structure itself and from episodes of transient acceleration during re-establishing of dynamical pressure balance after sudden radiative cooling onset. They are also easily excited by physical perturbations of the ambient medium. The widely mentioned Non-linear Thin Shell Instability, in contrast, in tests with physical perturbations is shown to have only limited chances to develop in real radiative shocks, as it seems to require a special spatial arrangement of fluctuations to be excited efficiently. The described phenomena also set new requirements on further simulations of the radiatively cooling shocks in order to be physically correct and free of numerical artefacts.

We present a new method to effectively select objects which may be low-mass active black holes (BHs) at galaxy centers using high-cadence optical imaging data, and our first spectroscopic identification of an active 2.7x10^6 Msun BH at z=0.164. This active BH was originally selected due to its rapid optical variability, from a few hours to a day, based on Subaru Hyper Suprime-Cam~(HSC) g-band imaging data taken with 1-hour cadence. Broad and narrow H-alpha and many other emission lines are detected in our optical spectra taken with Subaru FOCAS, and the BH mass is measured via the broad H-alpha emission line width (1,880 km s^{-1}) and luminosity (4.2x10^{40} erg s^{-1}) after careful correction for the atmospheric absorption around 7,580-7,720A. We measure the Eddington ratio to be as low as 0.05, considerably smaller than those in a previous SDSS sample with similar BH mass and redshift, which indicates one of the strong potentials of our Subaru survey. The g-r color and morphology of the extended component indicate that the host galaxy is a star-forming galaxy. We also show effectiveness of our variability selection for low-mass active BHs.

The discovery of the most luminous supernova ASASSN-15lh triggered a shock-wave in the supernova community. The three possible mechanisms proposed for the majority of other superluminous supernovae do not produce a realistic physical model for this particular supernova. In the present study we show the limiting luminosity available from a nickel-powered pair-instability supernova. We computed a few exotic nickel-powered explosions with a total mass of nickel up to 1500 solar masses. We used the hydrostatic configurations prepared with the GENEVA and MESA codes, and the STELLA radiative-transfer code for following the explosion of these models. We show that 1500 solar masses of radioactive nickel is needed to power a luminosity of 2x10^45 erg/s. The resulting light curve is very broad and incompatible with the shorter ASASSN-15lh time-scale. This rules out a nickel-powered origin of ASASSN-15lh. In addition, we derive a simple peak luminosity - nickel mass relation from our data, which may serve to estimate of nickel mass from observed peak luminosities.

For an accurate treatment of the shock wave propagation in high-energy astrophysical phenomena, such as supernova shock breakouts, gamma-ray bursts and accretion disks, knowledge of radiative transfer plays a crucial role. In this paper we consider one-dimensional (1D) special relativistic radiation hydrodynamics by solving the Boltzmann equation for radiative transfer. The structure of a radiative shock is calculated for a number of shock tube problems, including strong shock waves, and relativistic- and radiation-dominated cases. Calculations are performed using an iterative technique that consistently solves the equations of relativistic hydrodynamics and relativistic comoving radiative transfer. A comparison of radiative transfer solutions with the Eddington approximation and the M1 closure is made. A qualitative analysis of moment equations for radiation is performed and the conditions for the existence of jump discontinuity for non-relativistic cases are investigated numerically.

The properties of the first generation of stars and their supernova (SN) explosions remains unknown due to the lack of their actual observations. Recently many transient surveys are conducted and the feasibility of the detection of supernovae (SNe) of Pop III stars is growing. In this paper we study the multicolor light curves for a number of metal-free core-collapse SN models (25-100 M$_{\odot}$) to provide the indicators for finding and identification of first generation SNe. We use mixing-fallback supernova explosion models which explain the observed abundance patterns of metal poor stars. Numerical calculations of the multicolor light curves are performed using multigroup radiation hydrodynamic code STELLA. The calculated light curves of metal-free SNe are compared with non-zero metallicity models and several observed SNe. We have found that the shock breakout characteristics, the evolution of the photosphere's velocity, the luminosity, the duration and color evolution of the plateau - all the SN phases are helpful to estimate the parameters of SN progenitor: the mass, the radius, the explosion energy and the metallicity. We conclude that the multicolor light curves can be potentially used to identify first generation SNe in the current (Subaru/HSC) and future transient surveys (LSST, JWST). They are also suitable for identification of the low-metallicity SNe in the nearby Universe (PTF, Pan-STARRS, Gaia).

We present rapidly rising transients discovered by a high-cadence transient survey with Subaru telescope and Hyper Suprime-Cam. We discovered five transients at z=0.384-0.821 showing the rising rate faster than 1 mag per 1 day in the restframe near-ultraviolet wavelengths. The fast rising rate and brightness are the most similar to SN 2010aq and PS1-13arp, for which the ultraviolet emission within a few days after the shock breakout was detected. The lower limit of the event rate of rapidly rising transients is ~9 % of core-collapse supernova rates, assuming a duration of rapid rise to be 1 day. We show that the light curves of the three faint objects agree with the cooling envelope emission from the explosion of red supergiants. The other two luminous objects are, however, brighter and faster than the cooling envelope emission. We interpret these two objects to be the shock breakout from dense wind with the mass loss rate of ~10^{-3} Msun yr^{-1}, as also proposed for PS1-13arp. This mass loss rate is higher than that typically observed for red supergiants. The event rate of these luminous objects is >~1 % of core-collapse supernova rate, and thus, our study implies that more than ~1 % of massive stars can experience an intensive mass loss at a few years before the explosion.

It has been suggested that Type II supernovae with rapidly fading light curves (a.k.a. Type IIL supernovae) are explosions of progenitors with low-mass hydrogen-rich envelopes which are of the order of 1 Msun. We investigate light-curve properties of supernovae from such progenitors. We confirm that such progenitors lead to rapidly fading Type II supernovae. We find that the luminosity of supernovae from such progenitors with the canonical explosion energy of 1e51 erg and 56Ni mass of 0.05 Msun can increase temporarily shortly before all the hydrogen in the envelope recombines. As a result, a bump appears in their light curves. The bump appears because the heating from the nuclear decay of 56Ni can keep the bottom of hydrogen-rich layers in the ejecta ionized, and thus the photosphere can stay there for a while. We find that the light-curve bump becomes less significant when we make explosion energy larger (>~ 2e51 erg), 56Ni mass smaller (<~ 0.01 Msun), 56Ni mixed in the ejecta, or the progenitor radius larger. Helium mixing in hydrogen-rich layers makes the light-curve decline rates large but does not help reducing the light-curve bump. Because the light-curve bump we found in our light-curve models has not been observed in rapidly fading Type II supernovae, they may be characterized by not only low-mass hydrogen-rich envelopes but also higher explosion energy, larger degrees of 56Ni mixing, and/or larger progenitor radii than slowly fading Type II supernovae, so that the light-curve bump does not become significant.

An increasing number of so-called superluminous supernovae (SLSNe) are discovered. It is believed that at least some of them with slowly fading light curves originate in stellar explosions induced by the pair instability mechanism. Recent stellar evolution models naturally predict pair instability supernovae (PISNe) from very massive stars at wide range of metallicities (up to Z=0.006, Yusof et al. 2013). In the scope of this study we analyse whether PISN models can match the observational properties of SLSNe with various light curve shapes. Specifically, we explore the influence of different degrees of macroscopic chemical mixing in PISN explosive products on the resulting observational properties. We artificially apply mixing to the 250 Msun PISN evolutionary model from Kozyreva et al. (2014) and explore its supernova evolution with the one-dimensional radiation hydrodynamics code STELLA. The greatest success in matching SLSN observations is achieved in the case of an extreme macroscopic mixing, where all radioactive material is ejected into the hydrogen-helium outer layer. Such an extreme macroscopic redistribution of chemicals produces events with faster light curves with high photospheric temperatures and high photospheric velocities. These properties fit a wider range of SLSNe than non-mixed PISN model. Our mixed models match the light curves, colour temperature and photospheric velocity evolution of two well-observed SLSNe PTF12dam and LSQ12dlf. However, these models' extreme chemical redistribution may be hard to realise in massive PISNe. Therefore, alternative models such as the magnetar mechanism or wind-interaction may still to be favourable to interpret rapidly rising SLSNe.

We develop a time-dependent multi-group multidimensional relativistic radiative transfer code, which is required to numerically investigate radiation from relativistic fluids involved in, e.g., gamma-ray bursts and active galactic nuclei. The code is based on the spherical harmonic discrete ordinate method (SHDOM) that evaluates a source function including anisotropic scattering in spherical harmonics and implicitly solves the static radiative transfer equation with a ray tracing in discrete ordinates. We implement treatments of time dependence, multi-frequency bins, Lorentz transformation, and elastic Thomson and inelastic Compton scattering to the publicly available SHDOM code. Our code adopts a mixed frame approach; the source function is evaluated in the comoving frame whereas the radiative transfer equation is solved in the laboratory frame. This implementation is validated with various test problems and comparisons with results of a relativistic Monte Carlo code. These validations confirm that the code correctly calculates intensity and its evolution in the computational domain. The code enables us to obtain an Eddington tensor that relates first and third moments of intensity (energy density and radiation pressure) and is frequently used as a closure relation in radiation hydrodynamics calculations.

A model of a supernova explosion inside a dense extended hydrogenless envelope is proposed to explain the properties of the light curve for one of the superluminous supernovae PTF12dam. It is argued in the literature that the flux of this supernova rises too fast to be explained by the explosion model due to the instability associated with the electron-positron pair production (pair-instability supernova, PISNe), but it is well described by the models with energy input by a magnetar. We show that the PTF12dam-type supernovae can be explained without a magnetar in a model with a radiative shock in a dense circumstellar envelope that does not require an excessively large explosion energy.

We investigate the early-time light-curves of a large sample of 223 type II supernovae (SNe) from the Sloan Digital Sky Survey and the Supernova Legacy Survey. Having a cadence of a few days and sufficient non-detections prior to explosion, we constrain rise-times, i.e. the durations from estimated first to maximum light, as a function of effective wavelength. At restframe g-band (4722A), we find a distribution of fast rise-times with median of (7.5+/-0.3) days. Comparing these durations with analytical shock models of Rabinak and Waxman (2013); Nakar and Sari (2010) and hydrodynamical models of Tominaga et al. (2009), which are mostly sensitive to progenitor radius at these epochs, we find a median characteristic radius of less than 400 solar radii. The inferred radii are on average much smaller than the radii obtained for observed red supergiants (RSG). Investigating the post-maximum slopes as a function of effective wavelength in the light of theoretical models, we find that massive hydrogen envelopes are still needed to explain the plateaus of SNe II. We therefore argue that the SN II rise-times we observe are either a) the shock cooling resulting from the core collapse of RSG with small and dense envelopes, or b) the delayed and prolonged shock breakout of the collapse of a RSG with an extended atmosphere or embedded within pre-SN circumstellar material.

Recent progress in numerical simulations of thermonuclear supernova explosions brings up a unique opportunity in studying the progenitors of Type Ia supernovae. Coupling state-of-the-art explosion models with detailed hydrodynamical simulations of the supernova remnant evolution and the most up-to-date atomic data for X-ray emission calculations makes it possible to create realistic synthetic X-ray spectra for the supernova remnant phase. Comparing such spectra with high quality observations of supernova remnants could allow to constrain the explosion mechanism and the progenitor of the supernova. The present study focuses in particular on the oxygen emission line properties in young supernova remnants, since different explosion scenarios predict a different amount and distribution of this element. Analysis of the soft X-ray spectra from supernova remnants in the Large Magellanic Cloud and confrontation with remnant models for different explosion scenarios suggests that SNR 0509-67.5 could originate from a delayed detonation explosion and SNR 0519-69.0 from an oxygen-rich merger.

The dense shell method for the determination of distances to type IIn supernovae has been briefly reviewed. Applying our method to SN 2006gy, SN 2009ip, and SN 2010jl supernovae, we have obtained distances in excellent agreement with the previously known distances to the parent galaxies. The dense shell method is based on the radiation hydrodynamic model of a supernova. The method of the blackbody model, as well as the correctness of its application for simple estimates of distances from observation data, has been justified.

We consider consequences of hypothetical existence of baryo-dense stars created in the very early universe within an extension of Affleck-Dine scenario of baryogenesis. New constraints on the possible number of compact antimatter objects are derived. The contemporary observational data do not exclude significant amount of antimatter in the Galaxy (and in other galaxies) in the form of the baryo-dense low-massive stars.

The most massive stars on the asymptotic giant branch (AGB), so called super-AGB stars, are thought to produce supernovae (SNe) triggered by electron captures in their degenerate O+Ne+Mg cores. Super-AGB stars are expected to have slow winds with high mass-loss rates, so their wind density is high. The explosions of super-AGB stars are therefore presumed to occur in this dense wind. We provide the first synthetic light curves (LCs) for such events by exploding realistic electron-capture supernova (ecSN) progenitors within their super-AGB winds. We find that the early LC, i.e. before the recombination wave reaches the bottom of the H-rich envelope of SN ejecta (the plateau phase), is not affected by the dense wind. However, after the plateau phase, the luminosity remains higher when the super-AGB wind is taken into account. We compare our results to the historical LC of SN 1054, the progenitor of the Crab Nebula, and show that the explosion of an ecSN within an ordinary super-AGB wind can explain the LC features. We conclude that SN 1054 could have been a Type IIn SN without any extra extreme mass loss which was previously suggested to be necessary to account for its early high luminosity. We also show that our LCs match Type IIn SNe with an early plateau phase (`Type IIn-P') and suggest that they are ecSNe within super-AGB winds. Although some ecSNe can be bright in the optical spectral range due to the large progenitor radius, their X-ray luminosity from the interaction does not necessarily get as bright as other Type IIn SNe whose optical luminosities are also powered by the interaction. Thus, we suggest that optically-bright X-ray-faint Type IIn SNe can emerge from ecSNe. Optically-faint Type IIn SNe, such as SN 2008S, can also originate from ecSNe if their H-rich envelope masses are small. Some of them can be observed as `Type IIn-b' SNe due to the small H-rich envelope mass.

So called superluminous supernovae have been recently discovered in the local Universe. It appears possible that some of them originate from stellar explosions induced by the pair instability mechanism. Recent stellar evolution models also predict pair instability supernovae (PISNe) from very massive stars at fairly high metallicities (i.e. Z~0.004). We provide supernova (SN) models and synthetic light curves for two progenitor models, a 150 Msun red-supergiant and a 250 Msun yellow-supergiant at a metallicity of Z=0.001, for which the evolution from the main sequence to collapse, and the initiation of the PISN itself, has been previously computed in a realistic and self-consistent way. We are using the radiation hydrodynamics code STELLA to describe the SN evolution of both models over a time frame of about 500 days. We describe the shock-breakout phases of both SNe which are characterized by a higher luminosity, a longer duration and a lower effective temperature than those of ordinary SNeIIP. We derive the bolometric as well as the U, B, V, R and I light curves of our PISN models, which show a long-lasting plateau phase with maxima at Mbol=-19.3 mag and -21.3 mag for our lower and higher mass model, respectively. We also describe the photospheric composition and velocity as function of time. We conclude that the light curve of the explosion of our initially 150 Msun star resembles those of relatively bright SNeIIP, whereas its photospheric velocity at early times is smaller. Its 56Ni mass of 0.04 Msun also falls well into the range found in ordinary core collapse SNe. The light curve and photospheric velocity of our 250 Msun models has a striking resemblance with that of the superluminous SN2007bi, strengthening its interpretation as PISN. We conclude that PISNe may occur more frequently in the local universe than previously assumed.

We present results of a systematic study of the mass-loss properties of Type IIn supernova progenitors within decades before their explosion. We apply an analytic light curve model to 11 Type IIn supernova bolometric light curves to derive the circumstellar medium properties. We reconstruct the mass-loss histories based on the estimated circumstellar medium properties. The estimated mass-loss rates are mostly higher than 1e-3 Msun/yr and they are consistent with those obtained by other methods. The mass-loss rates are often found to be constantly high within decades before their explosion. This indicates that there exists some mechanism to sustain the high mass-loss rates of Type IIn supernova progenitors for at least decades before their explosion. Thus, the shorter eruptive mass loss events observed in some Type IIn supernova progenitors are not always responsible for creating their dense circumstellar media. In addition, we find that Type IIn supernova progenitors may tend to increase their mass-loss rates as they approach to the time of their explosion. We also show a detailed comparison between our analytic prediction and numerical results.

Recent discoveries of weak and fast optical transients raise the question of their origin. We investigate the minimum ejecta mass associated with core-collapse supernovae (SNe) of Type Ic. We show that mass transfer from a helium star to a compact companion can produce an ultra-stripped core which undergoes iron core collapse and leads to an extremely fast and faint SN Ic. In this Letter, a detailed example is presented in which the pre-SN stellar mass is barely above the Chandrasekhar limit, resulting in the ejection of only ~0.05-0.20 M_sun of material and the formation of a low-mass neutron star. We compute synthetic light curves of this case and demonstrate that SN 2005ek could be explained by our model. We estimate that the fraction of such ultra-stripped to all SNe could be as high as 0.001-0.01. Finally, we argue that the second explosion in some double neutron star systems (for example, the double pulsar PSR J0737-3039B) was likely associated with an ultra-stripped SN Ic.

A mechanism of creation of stellar-like objects in the very early universe, from the QCD phase transition till BBN and somewhat later, is studied. It is argued that in the considered process primordial black holes with masses above a few solar masses up to super-heavy ones could be created. This may explain an early quasar creation with evolved chemistry in surrounding medium and the low mass cutoff of the observed black holes. It is also shown that dense primordial stars can be created at the considered epoch. Such stars could later become very early supernovae and in particular high redshift gamma-bursters. In a version of the model some of the created objects can consist of antimatter.

The structure of the deflagration burning front in type Ia supernovae is considered. The parameters of the flame are obtained: its normal velocity and thickness. The results are in good agreement with the previous works of different authors. The problem of pulsation instability of the flame, subject to plane perturbations, is studied. First, with the artificial system with switched-off hydrodynamics the possibility of secondary reactions to stabilize the front is shown. Second, with account of hydrodynamics, realistic EOS and thermal conduction we can obtain pulsations when Zeldovich number was artificially increased. The critical Zeldovich numbers are presented. These results show the stability of the flame in type Ia supernovae against pulsations because its effective Zeldovich number is small.

We present an analytic model for bolometric light curves which are powered by the interaction between supernova ejecta and a dense circumstellar medium. This model is aimed at modeling Type IIn supernovae to determine the properties of their supernova ejecta and circumstellar medium. Our model is not restricted to the case of steady mass loss and can be applied broadly. We only consider the case in which the optical depth of the unshocked circumstellar medium is not high enough to affect the light curves. We derive the luminosity evolution based on an analytic solution for the evolution of a dense shell created by the interaction. We compare our model bolometric light curves to observed bolometric light curves of three Type IIn supernovae (2005ip, 2006jd, 2010jl) and show that our model can constrain their supernova ejecta and circumstellar medium properties. Our analytic model is supported by numerical light curves from the same initial conditions.

An electron-capture supernova (ECSN) is a core-collapse supernova (CCSN) explosion of a super-asymptotic giant branch (SAGB) star with a main-sequence mass $M_{\rm ms}\sim7-9.5M_\odot$. The explosion takes place in accordance with core bounce and subsequent neutrino heating and is a unique example successfully produced by first-principle simulations. This allows us to derive a first self-consistent multicolor light curves of a CCSN. Adopting the explosion properties derived by the first-principle simulation, i.e., the low explosion energy of $1.5\times10^{50}$ erg and the small $^{56}$Ni mass of $2.5\times10^{-3}M_\odot$, we perform a multigroup radiation hydrodynamics calculation of ECSNe and present multicolor light curves of ECSNe of SAGB stars with various envelope mass and hydrogen abundance. We demonstrate that a shock breakout has peak luminosity of $L\sim2\times10^{44}$ erg/s and can evaporate circumstellar dust up to $R\sim10^{17}$ cm for a case of carbon dust, that plateau luminosity and plateau duration of ECSNe are $L\sim10^{42}$ erg/s and $t\sim60-100$ days, respectively, and that a plateau is followed by a tail with a luminosity drop by $\sim4$ mag. The ECSN shows a bright and short plateau that is as bright as typical Type II plateau supernovae, and a faint tail that might be influenced by spin-down luminosity of a newborn pulsar. Furthermore, the theoretical models are compared with ECSN candidates: SN 1054 and SN 2008S. We find that SN 1054 shares the characteristics of the ECSNe. For SN 2008S, we find that its faint plateau requires an ECSN model with a significantly low explosion energy of $E\sim10^{48}$ erg.

We study thermal emission from circumstellar structures heated by gamma-ray burst (GRB) radiation and ejecta and calculate its contribution to GRB optical and X-ray afterglows using the modified radiation hydro-code small STELLA. It is shown that thermal emission originating in heated dense shells around the GRB progenitor star can reproduce X-ray plateaus (like observed in GRB 050904, 070110) as well as deviations from a power law fading observed in optical afterglows of some GRBs (e.g. 020124, 030328, 030429X, 050904). Thermal radiation pressure in the heated circumburst shell dominates the gas pressure, producing rapid expansion of matter similar to supenova-like explosions close to opacity or radiation flux density jumps in the circumburst medium. This phenomenon can be responsible for so-called supernova bumps in optical afterglows of several GRBs. Such a `quasi-supernova' suggests interpretation of the GRB-SN connection which does not directly involve the explosion of the GRB progenitor star.

We demonstrate the applicability of our new method (the Dense Shell Method or DSM) for the determination of astronomical distances by calculating the distance to SN2009ip. The distance to this supernova has been accurately determined in the standard approach via the cosmic distance ladder and has been found to be 20.4 Mpc. Our direct method, assuming the most reasonable parameter values, gives a very close result, namely 20.1+/-0.8 (68% CL) Mpc to SN2009ip.

We numerically investigate light curves (LCs) of shocked circumstellar shells which are suggested to reproduce the observed LC of superluminous SN 2006gy analytically. In the previous analytical model, the effects of the recombination and the bolometric correction on LCs are not taken into account. To see the effects, we perform numerical radiation hydrodynamic calculations of shocked shells by using STELLA, which can numerically treat multigroup radiation transfer with realistic opacities. We show that the effects of the recombination and the bolometric correction are significant and the analytical model should be compare to the bolometric LC instead of a single band LC. We find that shocked circumstellar shells have a rapid LC decline initially because of the adiabatic expansion rather than the luminosity increase and the shocked shells fail to explain the LC properties of SN 2006gy. However, our synthetic LCs are qualitatively similar to those of superluminous SN 2003ma and SN 1988Z and they may be related to shocked circumstellar shells.

Some features of the physics of radiation-dominated shock waves are discussed with emphasis on the peculiarities which are important for correct numerical modeling of shock breakouts in supernova. With account of those peculiarities, a number of models for different supernova types is constructed based on multigroup radiation transfer coupled to hydrodynamics. We describe the implementation of a new algorithm RADA, designed for modeling photon transfer at extremely-relativistic motions of matter, into our older code STELLA. The results of numerical simulations of light curves, and continuum spectra are presented. The influence of effects of photon scattering on electrons, of thermalization depth and of special relativity in transfer equation is considered. Some cases are demonstrated, when the appearance of hard X-ray emission is possible at the shock breakout. The necessary refinements in numerical algorithms for radiative transfer and hydrodynamics are pointed out. Prospects for using the results of numerical simulation to analyze and interpret available and future data from space observatories are discussed.

We show model light curves of superluminous supernova 2006gy on the assumption that the supernova is powered by the collision of supernova ejecta and its dense circumstellar medium. The initial conditions are constructed based on the shock breakout condition, assuming that the circumstellar medium is dense enough to cause the shock breakout within it. We perform a set of numerical light curve calculations by using a one-dimensional multigroup radiation hydrodynamics code STELLA. We succeeded in reproducing the overall features of the early light curve of SN 2006gy with the circumstellar medium whose mass is about 15 Msun (the average mass-loss rate ~ 0.1 Msun/yr). Thus, the progenitor of SN 2006gy is likely a very massive star. The density profile of the circumstellar medium is not well constrained by the light curve modeling only, but our modeling disfavors the circumstellar medium formed by steady mass loss. The ejecta mass is estimated to be comparable to or less than 15 Msun and the explosion energy is expected to be more than 4e51 erg. No 56Ni is required to explain the early light curve. We find that the multidimensional effect, e.g., the Rayleigh-Taylor instability, which is expected to take place in the cool dense shell between the supernova ejecta and the dense circumstellar medium, is important in understanding supernovae powered by the shock interaction. We also show the evolution of the optical and near-infrared model light curves of high-redshift superluminous supernovae. They can be potentially used to identify SN 2006gy-like superluminous supernovae in the future optical and near-infrared transient surveys.

We introduce a novel approach, a Dense Shell Method (DSM), for measuring distances for cosmology. It is based on original Baade idea to relate absolute difference of photospheric radii with photospheric velocity. We demonstrate that this idea works: the new method does not rely on the Cosmic Distance Ladder and gives satisfactory results for the most luminous Type IIn Supernovae. This allows one to make them good primary distance indicators for cosmology. Fixing correction factors for illustration, we obtain with this method the median distance of 68^{+19}_{-15} (68%CL) Mpc to SN 2006gy and median Hubble parameter 79^{+23}_{-17} (68%CL) km/s/Mpc.

Several highly luminous Type Ia supernovae (SNe Ia) have been discovered. Their high luminosities are difficult to explain with the thermonuclear explosions of the Chandrasekhar-mass white dwarfs (WDs). In the present study, we estimate the progenitor mass of SN 2009dc, one of the extremely luminous SNe Ia, using the hydrodynamical models as follows. Explosion models of super-Chandrasekhar-mass (super-Ch-mass) WDs are constructed, and multi-color light curves (LCs) are calculated. The comparison between our calculations and the observations of SN 2009dc suggests that the exploding WD has a super-Ch mass of 2.2-2.4 solar masses, producing 1.2-1.4 solar masses of Ni-56, if the extinction by its host galaxy is negligible. If the extinction is significant, the exploding WD is as massive as \sim2.8 solar masses, and \sim1.8 solar masses of Ni-56 is necessary to account for the observations. Whether the host-galaxy extinction is significant or not, the progenitor WD must have a thick carbon-oxygen layer in the outermost zone (20-30% of the WD mass), which explains the observed low expansion velocity of the ejecta and the presence of carbon. Our estimate on the mass of the progenitor WD, especially for the extinction-corrected case, is challenging to the current scenarios of SNe Ia. Implications on the progenitor scenarios are also discussed.

CCD UBVRI photometry is presented for type IIb SN 2011dh for about 300 days. The main photometric parameters are derived and the comparison with SNe of similar types is reported. The light curves are similar to those for SN IIb 2008ax, but the initial flash is stronger and very short, and there are humps on the light curves in U and B at the onset of linear decline. Preliminary modeling is carried out, and the results are compared to the quasi-bolometric light curve and to the light curves in UBVRI bands.

We calculate multicolor light curves (LCs) of supernovae (SNe) from red supergiants (RSGs) exploded within dense circumstellar medium (CSM). Multicolor LCs are calculated by using a multi-group radiation hydrodynamics code STELLA. If CSM is dense enough, the shock breakout signal is delayed and smeared by CSM and kinetic energy of SN ejecta is efficiently converted to thermal energy which is eventually released as radiation. We find that explosions of RSGs are affected by CSM in early epochs when mass-loss rate just before the explosions is higher than 10^{-4} Msun/yr. Their characteristic features are that the LC has a luminous round peak followed by a flat LC, that multicolor LCs are simultaneously bright in ultraviolet and optical at the peak, and that photospheric velocity is very low at these epochs. We calculate LCs for various CSM conditions and explosion properties, i.e., mass-loss rates, radii of CSM, density slopes of CSM, explosion energies of SN ejecta, and SN progenitors inside, to see their influence on LCs. We compare our model LCs to those of ultraviolet-bright Type IIP SN 2009kf and show that the mass-loss rate of the progenitor of SN 2009kf just before the explosion is likely to be higher than 10^{-4} Msun/yr. Combined with the fact that SN 2009kf is likely to be an energetic explosion and has large 56Ni production, which implies that the progenitor of SN 2009kf is a massive RSG, our results indicate that there could be some mechanism to induce extensive mass loss in massive RSGs just before their explosions.

We present numerical models for supernova remnant evolution, using a new version of the hydrodynamical code SUPREMNA. We added cosmic ray diffusion equation to the code scheme, employing two-fluid approximation. We investigate the dynamics of the simulated supernova remnants with different values of cosmic ray acceleration efficiency and diffusion coefficient. We compare the numerical models with observational data of Tycho's and SN1006 supernova remnants. We find models which reproduce the observed locations of the blast wave, contact discontinuity, and reverse shock for the both remnants, thus allowing us to estimate the contribution of cosmic ray particles into total pressure and cosmic-ray energy losses in these supernova remnants. We derive that the energy losses due to cosmic rays escape in Tycho's supernova remnant are 10-20% of the kinetic energy flux and 20-50% in SN1006.

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 $z\geq1.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.

The propagation of thermonuclear flame in presupernovae Ia is considered. Front parameters are obtained, some speculations on front stability are presented.

Motivated by a recent discovery of Supernova 2010gx and numerical results of Fryer et al.(2010), we simulate light curves for several type I supernova models, enshrouded by dense circumstellar shells, or "super-wind", rich in carbon and oxygen and having no hydrogen. We demonstrate that the most luminous events like SN2010gx can be explained by those models at moderate explosion energies (2-3) foe if the total mass of SN ejecta and a shell is (3-5) Msun and the radius of the shell is ~10^{16} cm.

Electrons and electron neutrinos in the inner core of the core-collapse supernova are highly degenerate and therefore numerous during a few seconds of explosion. In contrast, leptons of other flavors are non-degenerate and therefore relatively scarce. This is due to lepton flavor conservation. If this conservation law is broken by some non-standard interactions, electron neutrinos are converted to muon and tau-neutrinos, and electrons -- to muons. This affects the supernova dynamics and the supernova neutrino signal in several ways. In particular, the total neutrino luminosity in the first second of the collapse is increased due to the larger free path of the non-electron neutrinos. This effect may have important consequences as the increase of the neutrino luminosity is known to facilitate the explosion. We consider an extension of the Standard Model by scalar bileptons which mediate lepton flavor violation. It is shown that in case of TeV-mass bileptons the electron fermi gas is equilibrated with non-electron species inside the inner supernova core at a time-scale of order of (1-100) ms. In particular, a scalar triplet which generates neutrino masses through the see-saw type II mechanism is considered. Non-observation of rare decays and data on neutrino mixing and neutrino masses restrict possible lepton flavor violation effects in this case. However a region in the parameter space of the model exists which fits all the experimental constraints and provides lepton flavor violation sufficient for observable effects in supernova.

Extending the population synthesis method to isolated young cooling white dwarfs we are able to confront our model assumptions with observations made in ROSAT All-Sky Survey (Fleming et al., 1996). This allows us to check model parameters such as evolution of spectra and separation of heavy elements in DA WD envelopes. It seems like X-ray spectrum temperature of these objects is given by the formula T_{X-ray} = min(T_eff, T_max). We have obtained DA WD's birth rate and upper limit of the X-ray spectrum temperature: DA birth rate $= 0.61\times 10^{-12}$ in cubic parsec per year and T_max = 41000 K. These values are in good correspondence with values obtained by other authors (Liebert et al., 2004; Wolff et al., 1996). From this fact we also conclude that our population synthesis method is applicable to the population of close-by isolated cooling white dwarfs as well as to the population of the isolated cooling neutron stars.

We present methodology and results of numerical modeling of the interaction of GRB prompt emission with the circumburst medium using a modified version of the multi-group radiation hydrocode STELLA. The modification includes the nonstationary photoionization, the photoionization heating and the Compton heating along with the hydrodynamics and radiation transfer. The lightcurves and spectra of the outcoming gamma-ray, X-ray and optical emission are calculated for a set of models (shells) of the circumburst environment, which differ in dimensions, density, density profile, composition, temperature. In some cases total bolometric and optical luminosities can reach 10^47 and 10^43 erg/s respectively. These effects can be responsible for irregularities which are seen on lightcurves of some GRB's X-ray and optical afterglows.

Electrons and electron neutrinos in the inner core of the core-collapse supernova are highly degenerate and therefore numerous during a few seconds of explosion. In contrast, leptons of other flavors are non-degenerate and therefore relatively scarce. This is due to lepton flavor conservation. If this conservation law is broken by some non-standard interactions, electron neutrinos are converted to muon and tau-neutrinos, and electrons - to muons. This affects the supernova dynamics and the supernova neutrino signal. We consider lepton flavor violating interactions mediated by scalar bileptons, i.e. heavy scalars with lepton number 2. It is shown that in case of TeV-mass bileptons the electron fermi gas is equilibrated with non-electron species inside the inner supernova core at a time-scale of order of (1-100) ms. In particular, a scalar triplet which generates neutrino masses through the see-saw type II mechanism is considered. It is found that supernova core is sensitive to yet unprobed values of masses and couplings of the triplet.

We point out that confusion sometimes arises when using a chemical potential in plasma with Coulomb interactions. The results of our consideration are applied to the discussion of nuclear reactions screening. Finally, we present a transparent derivation of the thermal conductivity coefficient of a degenerate electron gas.

The theory of radiative transfer provides the link between the physical conditions in an astrophysical object and the observable radiation which it emits. Thus accurately modelling radiative transfer is often a necessary part of testing theoretical models by comparison with observations. We describe a new radiative transfer code which employs Monte Carlo methods for the numerical simulation of radiation transport in expanding media. We discuss the application of this code to the calculation of synthetic spectra and light curves for a Type Ia supernova explosion model and describe the sensitivity of the results to certain approximations made in the simulations.

CCD UBVRI photometry is presented for type IIb SN 2008ax for about 320 days. The photometric behavior is typical for core-collapse SNe with low amount of hydrogen. The main photometric parameters are derived and the comparison with SNe of similar types is reported. Preliminary modeling is carried out, and the results are compared to the observed light curves. The main parameters of the hydrodynamical models are close to those used for SN IIb 1993J.

Shock breakout is the brightest radiative phenomenon in a Type II supernova (SN). Although it was predicted to be bright, the direct observation is difficult due to the short duration and X-ray/ultraviolet-peaked spectra. First entire observations of the shock breakouts of Type II Plateau SNe (SNe IIP) were reported in 2008 by ultraviolet and optical observations by the {\it GALEX} satellite and supernova legacy survey (SNLS), named SNLS-04D2dc and SNLS-06D1jd. We present multicolor light curves of a SN IIP, including the shock breakout and plateau, calculated with a multigroup radiation hydrodynamical code {\sc STELLA} and an evolutionary progenitor model. The synthetic multicolor light curves reproduce well the observations of SNLS-04D2dc. This is the first study to reproduce the ultraviolet light curve of the shock breakout and the optical light curve of the plateau consistently. We conclude that SNLS-04D2dc is the explosion with a canonical explosion energy $1.2\times10^{51}$ ergs and that its progenitor is a star with a zero-age main-sequence mass $20M_\odot$ and a presupernova radius $800R_\odot$. The model demonstrates that the peak apparent $B$-band magnitude of the shock breakout would be $m_{\rm B}\sim26.4$ mag if a SN being identical to SNLS-04D2dc occurs at a redshift $z=1$, which can be reached by 8m-class telescopes. The result evidences that the shock breakout has a great potential to detect SNe IIP at $z\gsim1$.

Neutrino with magnetic moment can experience a chirality flip while scattering off charged particles. This effect may lead to important consequences for the dynamics and the neutrino signal of the core-collapse supernova. It is known that if neutrino is a Dirac fermion, then nu_L->nu_R transition induced by the chirality flip leads to the emission of sterile right-handed neutrinos. The typical energies of these neutrinos are rather high, E ~ (100-200)MeV. Neutrino spin precession in the magnetic field either inside the collapsing star or in the interstellar space may lead to the backward transition, nu_R->nu_L. Both possibilities are known to be interesting. In the former case high-energy neutrinos can deliver additional energy to the supernova envelope, which can help the supernova to explode. In the latter case high-energy neutrinos may be detected simultaneously with the "normal" supernova neutrino signal, which would be a smoking gun for the Dirac neutrino magnetic moment. We report the results of the calculation of the supernova right-handed neutrino luminosity up to 250 ms after bounce. They allow to refine the estimates of the energy emitted in right-handed neutrinos. Also the sensitivity of water Cherenkov detectors to the Dirac neutrino magnetic moment is estimated. For mu_Dirac=10^{-13}mu_B Super-Kamiokande is expected to detect at least few high-energy events from a galactic supernova explosion. Also we briefly discuss the case of Majorana neutrino magnetic moment. It is pointed out that spin flips may quickly equilibrate electron neutrinos with non-electron antineutrinos if mu_Majorana~10^{-12}mu_B. This may lead to various consequences for supernova physics.

A few remarks on Dark Matter (DM) models are presented. An example is Mirror Matter which is the oldest but still viable DM candidate, perhaps not in the purest form. It can serve as a test-bench for other analogous DM models, since the properties of macroscopic objects are quite firmly fixed for Mirror Matter. A pedagogical derivation of virial theorem is given and it is pointed out that concepts of virial velocity or virial temperature are misleading for some cases. It is shown that the limits on self-interaction cross-sections derived from observations of colliding clusters of galaxies are not real limits for individual particles if they form macroscopic bodies. The effect of the heating of interstellar medium by Mirror Matter compact stars is very weak but may be observable. The effect of neutron star heating by accretion of M-baryons may be negligible. Problems of MACHOs as Mirror Matter stars are touched upon.

The equation of state and composition of matter are calculated for conditions typical for pre-collapse and early collapse stages in core collapse supernovae. The composition is evaluated under the assumption of nuclear statistical equilibrium, when the matter is considered as an `almost' ideal gas with corrections due to thermal excitations of nuclei, to free nucleon degeneracy, and to Coulomb and surface energy corrections. The account of these corrections allows us to obtain the composition for densities a bit below the nuclear matter density. Through comparisons with the equation of state (EOS) developed by Shen et al. we examine the approximation of one representative nucleus used in most of recent supernova EOS's. We find that widely distributed compositions in the nuclear chart are different, depending on the mass formula, while the thermodynamical quantities are quite close to those in the Shen's EOS.

The extremely luminous supernova SN2006gy is explained in the same way as other SNIIn events: light is produced by a radiative shock propagating in a dense circumstellar envelope formed by a previous weak explosion. The problems in the theory and observations of multiple-explosion SNe IIn are briefly reviewed.

We report on X-ray observations of the supernova remnant 0509-67.5 in the Large Magellanic Cloud with XMM-Newton X-ray observatory. We use the imaging spectroscopy (EPIC) and Reflective Grating Spectrometer (RGS) data to investigate properties of the remnant and its environment. The X-ray spectra were analyzed with SPEX software package. In addition to this we performed a numerical hydrodynamic simulation of the remnant. The EPIC data show prominent Fe K line emission, but the deduced overall amount of iron in the shocked ejecta is low. The data also show that the remnant has an asymmetric ejecta structure: the bright southwest region of the remnant shows an overabundance of metals. The analysis of the RGS spectrum shows that the remnant has a high lines velocity broadening of ~5000 km/s. We found a hydrodynamical model for the remnant with basic hydrodynamical and spectral parameters similar to the observed ones. The data analysis show that the reverse shock just recently reached iron layers of the ejecta. The brightness enhancement in the southwest region could be a sign of an asymmetric explosion or it could be the result of a density enhancement of the interstellar medium. We constructed numerical models which are in good agreement with the observations, with circumstellar density of 3e-25 g/cm^3, age of ~400 years, velocities of ~5000 km/s and an electron to ion temperature ratio of 0.01.

The extremely luminous supernova SN 2006gy challenges the traditional view that the collapse of a stellar core is the only mechanism by which a massive star makes a supernova, because it seems too luminous by more than a factor of ten. Here we report that the brightest supernovae in the modern Universe arise from collisions between shells of matter ejected by massive stars that undergo an interior instability arising from the production of electron-positron pairs. This "pair instability" leads to explosive burning that is insufficient to unbind the star, but ejects many solar masses of the envelope. After the first explosion, the remaining core contracts and searches for a stable burning state. When the next explosion occurs, several solar masses of material are again ejected, which collide with the earlier ejecta. This collision can radiate 1E50 erg of light, about a factor of ten more than an ordinary supernova. Our model is in good agreement with the observed light curve for SN 2006gy and also shows that some massive stars can produce more than one supernova-like outburst.

A simulation of the thermonuclear explosion of a Chandrasekhar-mass C+O white dwarf, the most popular scenario of a type Ia supernova (SN Ia), is presented. The underlying modeling is pursued in a self-consistent way, treating the combustion wave as a turbulent deflagration using well tested methods developed for laboratory combustion and based on the concept of `large eddy simulations' (LES). Such consistency requires to capture the onset of the turbulent cascade on resolved scales. This is achieved by computing the dynamical evolution on a 1024$^3$ moving grid, which resulted in the best-resolved three-dimensional SN Ia simulation carried out thus far, reaching the limits of what can be done on present supercomputers. Consequently, the model has no free parameters other than the initial conditions at the onset of the explosion, and therefore it has considerable predictive power. Our main objective is to determine to which extent such a simulation can account for the observations of normal SNe Ia. Guided by previous simulations with less resolution and a less sophisticated flame model, initial conditions were chosen that yield a reasonably strong explosion and a sufficient amount of radioactive nickel for a bright display. We show that observables are indeed matched to a reasonable degree. In particular, good agreement is found with the light curves of normal SNe Ia. Moreover, the model reproduces the general features of the abundance stratification as inferred from the analysis of spectra. This indicates that it captures the main features of the explosion mechanism of SNe Ia. However, we also show that even a seemingly best-choice pure deflagration model has shortcomings that indicate the need for a different mode of nuclear burning at late times, perhaps the transition to a detonation at low density.

We present extensive u'g'r'i'BVRIYJHKs photometry and optical spectroscopy of SN 2005hk. These data reveal that SN 2005hk was nearly identical in its observed properties to SN 2002cx, which has been called ``the most peculiar known type Ia supernova.'' Both supernovae exhibited high ionization SN 1991T-like pre-maximum spectra, yet low peak luminosities like SN 1991bg. The spectra reveal that SN 2005hk, like SN 2002cx, exhibited expansion velocities that were roughly half those of typical type Ia supernovae. The R and I light curves of both supernovae were also peculiar in not displaying the secondary maximum observed for normal type Ia supernovae. Our YJH photometry of SN 2005hk reveals the same peculiarity in the near-infrared. By combining our optical and near-infrared photometry of SN 2005hk with published ultraviolet light curves obtained with the Swift satellite, we are able to construct a bolometric light curve from ~10 days before to ~60 days after B maximum. The shape and unusually low peak luminosity of this light curve, plus the low expansion velocities and absence of a secondary maximum at red and near-infrared wavelengths, are all in reasonable agreement with model calculations of a 3D deflagration which produces ~0.25 M_sun of 56Ni.

Type Ia supernovae (SNe Ia) are one of the major tools to determine the cosmological parameters. Utilizing them as distance indicators, it is possible to geometrically survey the universe. To this end, the intrinsic scatter in the luminosities of these events needs to be calibrated using empirical relations between observables. A theoretical explanation for these relations is still lacking and can only be provided by a sound understanding of the mechanism of SNe Ia. Recently there has been significant progress in modeling SNe Ia. We report on numerical simulations of the explosion process, compare the results with observations of nearby SNe Ia, and discuss current uncertainties of the models. The presented simulations shed some light on the origin of the diversity of SNe Ia. Such simulations will pave the way towards an understanding of SN Ia diversities and correlations of their properties and ultimately provide a tool to validate the cosmological implications of SN Ia distance measurements.

The diversity of Type Ia supernova (SN Ia) photometry is explored using a grid of 130 one-dimensional models. It is shown that the observable properties of SNe Ia resulting from Chandrasekhar-mass explosions are chiefly determined by their final composition and some measure of ``mixing'' in the explosion. A grid of final compositions is explored including essentially all combinations of 56Ni, stable ``iron'', and intermediate mass elements that result in an unbound white dwarf. Light curves (and in some cases spectra) are calculated for each model using two different approaches to the radiation transport problem. Within the resulting templates are models that provide good photometric matches to essentially the entire range of observed SNe Ia. On the whole, the grid of models spans a wide range in B-band peak magnitudes and decline rates, and does not obey a Phillips relation. In particular, models with the same mass of 56Ni show large variations in their light curve decline rates. We identify the physical parameters responsible for this dispersion, and consider physically motivated ``cuts'' of the models that agree better with the Phillips relation. For example, models that produce a constant total mass of burned material of 1.1 +/- Msun do give a crude Phillips relation, albeit with much scatter. The scatter is further reduced if one restricts that set to models that make 0.1 to 0.3 Msun of stable iron and nickel isotopes, and then mix the ejecta strongly between the center and 0.8 Msun. We conclude that the supernovae that occur most frequently in nature are highly constrained by the Phillips relation and that a large part of the currently observed scatter in the relation is likely a consequence of the intrinsic diversity of these objects.

We present synthetic bolometric and broad-band UBVRI light curves of SNe Ia, for four selected 3-D deflagration models of thermonuclear supernovae. The light curves are computed with the 1-D hydro code STELLA, which models (multi-group time-dependent) non-equilibrium radiative transfer inside SN ejecta. Angle-averaged results from 3-D hydrodynamical explosion simulations with the composition determined in a nucleosynthetic postprocessing step served as the input to the radiative transfer model. The predicted model UBV light curves do agree reasonably well with the observed ones for SNe Ia in the range of low to normal luminosities, although the underlying hydrodynamical explosion models produced only a modest amount of radioactive Ni56 and relatively low kinetic energy in the explosion. The evolution of predicted B and V fluxes in the model with a Ni56 mass of 0.42 M_sun follows the observed decline rate after the maximum very well, although the behavior of fluxes in other filters somewhat deviates from observations, and the bolometric decline rate is a bit slow. Using our models, we check the validity of Arnett's rule and the accuracy of the procedure for extracting the Ni56 mass from the observed light curves. We find that the comparison between theoretical light curves and observations provides a useful tool to validate SN Ia models. The steps necessary to improve the agreement between theory and observations are set out.

We present u'g'r'i'BV photometry and optical spectroscopy of the Type Ib/Ic SN 2005bf covering the first ~100 days following discovery. The u'g'BV light curves displayed a double-peaked morphology, which is among Type Ib/Ic supernovae. The bolometric light curve indicates that SN 2005bf was a remarkably luminous event. Spectroscopically, SN 2005bf underwent a unique transformation from a Type Ic-like event at early times to a typical Type Ib supernova at later phases. The initial maximum in u'g'BV was accompanied by the presence in the spectrum of high velocity absorption lines of Fe II, Ca II, and H I. The photospheric velocity derived from spectra at early epochs was unusually low compared with ordinary Type Ib supernovae. We describe one-dimensional computer simulations which attempt to account for these remarkable properties. The most favored model is that of a very energetic (2 x 10^{51} erg), asymmetric explosion of a massive (8.3 M_sun) Wolf-Rayet WN star that had lost most of its hydrogen envelope. We speculate that an unobserved relativistic jet was launched producing a two-component explosion consisting of 1) a polar explosion containing a small fraction of the total mass and moving at high velocity, and 2) the explosion of the rest of the star. At first, only the polar explosion is observed, producing the initial maximum and the high velocity absorption-line spectrum resembling a Type Ic event. At late times, this fast-moving component becomes optically-thin, revealing the slower-moving explosion of the rest of the star and transforming the observed spectrum to that of a typical Type Ib supernova. If this scenario is correct, then SN 2005bf is the best example to date of a transition object between normal Type Ib/Ic supernovae and gamma ray bursts.

We investigated the kinematics of the pulsar wind nebula (PWN) associated with PSR B1951+32 in the old supernova remnant CTB 80 using the Fabry-Perot interferometer of the 6m Special Astrophysical Observatory telescope. In addition to the previously known expansion of the system of bright filaments with a velocity of 100-200km/s, we detected weak high-velocity features in the H-alpha line at least up to velocities of 400-450km/s. We analyzed the morphology of the PWN in the H-alpha, [SII], and [OIII] lines using HST data and discuss its nature. The shape of the central filamentary shell, which is determined by the emission in the [OIII] line and in the radio continuum, is shown to be consistent with the bow-shock model for a significant (about 60 degrees) inclination of the pulsar's velocity vector to the plane of the sky. In this case, the space velocity of the pulsar is twice higher than its tangential velocity, i.e., it reaches ~500 km/s, and PSR B1951+32 is the first pulsar whose line-of-sight velocity (of about 400 km/s) has been estimated from the PWN observations. The shell-like H-alpha-structures outside the bow shock front in the east and the west may be associated with both the pulsar's jets and the pulsar-wind breakthrough due to the layered structure of the extended CTB 80 shell.

We examine and analyze the physical processes that should be taken into account when modeling young type-Ia SNRs, with ages of several hundred years. It is shown, that energy losses in the metal-rich ejecta can be essential for remnants already at this stage of evolution. The influence of electron thermal conduction and the rate of the energy exchange between electrons and ions on the temperature distribution and the X-radiation from such remnants is studied. The data for Tycho SNR from the XMM-Newton X-ray telescope have been employed for the comparison of calculations with observations.

We compare parameters of well-observed type II SN1999em derived by M.Hamuy and D.Nadyozhin based on Litvinova-Nadyozhin (1985) analytic fits with those found from the simulations with our radiative hydro code Stella. The difference of SN parameters is quite large for the long distance scale. The same code applied to models of SN1993J allows us to estimate systematic errors of extracting foreground extinction toward SN1993J suggested by Clocchiatti et al. (1995). A new implicit two-temperature hydro code code Supremna is introduced which self-consistently takes into account the kinetics of ionization, electron thermal conduction, and radiative losses for predicting X-ray spectra of young supernova remnants such as Tycho and Kepler.

Extensive light and colour curves for the Type Ia supernova SN 2002er are presented as part of the European Supernova Collaboration. We have collected UBVRI photometry from ten different telescopes covering the phases from 7 days before until 619 days after maximum light. Corrections for the different instrumental systems and the non-thermal spectrum of the supernova (S-corrections) have been applied. With the densely sampled light curves we can make detailed comparisons to other well-observed objects. SN 2002er most closely resembles SN 1996X after maximum, but clearly shows a different colour evolution before peak light and a stronger shoulder in V and R bands compared to other well-observed SNe Ia. In particular, the rise time appears to be longer than what is expected from rise-time vs.decline-rate relation. We use several methods to determine the reddening towards SN 2002er based on the colour evolution at near peak and at late phases. The uvoir (bolometric) light curve shows great similarity with SN 1996X, but also indications of a higher luminosity, longer rise time and a more pronounced shoulder 25 days past maximum. The interpretation of the light curves was done with two independent light curve codes. Both find that given the luminosity of SN 2002er the 56Ni mass exceeds 0.6 Msun with prefered values near 0.7 Msun. Uncertainties in the exact distance to SN 2002er are the most serious limitation of this measurement. The light curve modelling also indicates a high level of mixing of the nickel in the explosion of SN 2002er.

We present and analyse spectra of the Type IIn supernova 1994W obtained between 18 and 203 days after explosion. During the luminous phase (first 100 d) the line profiles are composed of three major components: (i) narrow P-Cygni lines with the absorption minima at -700 km/s; (ii) broad emission lines with BVZI ~4000 km/s; and (iii) broad, smooth wings, most apparent in H-alpha. These components are identified with an expanding circumstellar (CS) envelope, shocked cool gas in the forward post-shock region, and multiple Thomson scattering in the CS envelope, respectively. The absence of broad P-Cygni lines from the supernova is the result of the formation of an optically thick, cool, dense shell at the interface of the ejecta and the CS envelope. We model the supernova deceleration and Thomson scattering wings to recover the density, radial extent and Thomson optical depth of the CS envelope during the first month. We reproduce the light curve with a hydrodynamical model and find it to be powered by a combination of internal energy leakage after the explosion of an extended pre-supernova (~10^15 cm) and luminosity from circumstellar interaction. We recover the pre-explosion kinematics of the CS envelope: it is close to homologous expansion with outer velocity ~1100 km/s and a kinematic age of ~1.5 yr. The CS envelope's high mass and kinetic energy, combined with its small age, strongly suggest that the CS envelope was explosively ejected about 1.5 yr before the supernova explosion.

We performed accurate numerical calculations of angle-, time-, and frequency-dependent radiative transfer for the relativistic motion of matter in gamma-ray burst (GRB) models. Our technique for solving the transfer equation, which is based on the method of characteristics, can be applied to the motion of matter with a Lorentz factor up to 1000. The effect of synchrotron self-absorption is taken into account. We computed the spectra and light curves from electrons with a power-law energy distribution in an expanding relativistic shock and compare them with available analytic estimates. The behavior of the optical afterglows from GRB 990510 and GRB 000301c is discussed qualitatively.

Time-dependent thermal effects should accompany standard non-thermal afterglows of GRB when gamma-rays pass through inhomogeneous surroundings of the GRB site. Thermal relaxation of an optically thin plasma is calculated using time-dependent collisional ionization of the plasma ion species. X-ray emission lines are similar to those found in the fading X-ray afterglow of GRB 011211. Thermal relaxation of clouds or shells around the GRB site could also contribute to the varying late optical GRB afterglows, such as in GRB021004 and GRB030329.

We present and analyse spectra of the Type IIn supernova 1994W obtained between 18 and 202 days after explosion. During the first 100 days the line profiles are composed of three major components: (i) narrow P Cygni lines with absorption minima at -700 km/s; (ii) broad emission lines with BVZI ~4000 km/s; (iii) broad wings most apparent in H-alpha. These components are identified with the expanding circumstellar (CS) envelope (Sollerman, Cumming & Lundqvist 1998), shocked cool gas in the forward postshock region, and multiple Thomson scattering in the CS envelope, respectively. The absence of broad P Cygni lines from the supernova (SN) is the result of the formation of an optically thick, cool, dense shell at the interface of the ejecta and the CS envelope. Models of the SN deceleration and Thomson scattering wings are used to recover the Thomson optical depth of the CS envelope and its radial extent, 4E15 cm. The light curve, which we reproduce by a hydrodynamical model, is powered by a combination of internal energy leakage after the explosion of an extended presupernova (~1E15 cm) and luminosity from circumstellar interaction. We recover the pre-explosion kinematics of the CS envelope and find it to be close to homologous expansion with outmost velocity ~1100 km/s and a kinematic age of \~1.5 yr. The high mass and kinetic energy of the CS envelope strongly suggest that the CS envelope was explosively ejected only a few years before explosion.

The successful theoretical supernova explosion models should be able to explain any features of the emission from supernovae at any evolutionary stage. We check several models from two different points of view. With the multi-frequency radiation hydro code STELLA we calculate gamma-ray, bolometric and broad-band UBVI light curves. Then we use the same models to calculate the emission from young supernova remnants. Here we present new plots for gamma-ray luminosity from several SN Ia models and recomputations of bolometric and UBVRI light curves of model 13C for SN 1993J.

Time-dependent thermal X-ray spectra are calculated from physically plausible conditions around GRB. It is shown that account for time-dependent ionization processes strongly affects the observed spectra of hot rarefied plasma. These calculations may provide an alternative explanation to the observed X-ray lines of early GRBs afterglows (such as GRB 011211). Our technique will allow one to obtain independent constraints on the GRB collimation angle and on the clumpiness of circumstellar matter.

Supernovae of type Ia (SNe Ia) are very important for cosmography. To exclude systematic effects in linking the observed light of distant SNe Ia to the parameters of cosmological models, one has to understand the nature of supernova outbursts and to build accurate algorithms for predicting their emission. We review the recent progress of modeling the propagation of nuclear flame subject to numerous hydrodynamic instabilities inherent to the flame front. The Rayleigh-Taylor (RT) instability is the main process governing the corrugation of the front on the largest scales, while on the smallest scales the front propagation is controlled by the Landau-Darrieus instability. Based on several hydrodynamic explosion models, we predict the broad-band UBVI and bolometric light curves of SNe Ia, using our 1D-hydro code which models multi-group time-dependent non-equilibrium radiative transfer inside SN ejecta. We employ our new corrected treatment for line opacity in the expanding medium, which is important especially in UV and IR bands. The results are compared with the observed light curves. Especially interesting is a recent 3D-deflagration model computed at MPA, Garching, by M. Reinecke et al.

We simulate shock propagation in a wide range of core-collapsing presupernovae: from compact WR stars exploding as SNe Ib/c through very extended envelopes of the narrow-line SNe IIn. We find that the same physical phenomenon of radiating shocks can produce outbursts of X-ray radiation (with photon energy $3kT \sim 1$ keV) lasting only a second in SNe Ib/c, as well as a very high flux of visual light, lasting for months, in SNe IIn.

We present theoretical UBVI- and bolometric light curves of SNe Ia for several explosion models, computed with our multi-group radiation hydro code. We employ our new corrected treatment for line opacity in the expanding medium. The results are compared with observed light curves. Our goal is to find the most viable thermonuclear SN model that gives good fits not only to a typical SN Ia light curves, but also to X-ray observations of young SN Ia remnants. It appears that classical 1D SN Ia models, such as deflagration model W7 and delayed detonation one DD4, fit the light curves not so good as a new 3D deflagration model by Reinecke et al. (which is averaged over angles for our LC modelling). This model seems good also in reproducing X-ray observations of Tycho SNR. We believe that the main feature of this model which allows us to get correct radiation during the first month, as well as after a few hundred years, when an SNR forms, is strong mixing that pushes material enriched in iron and nickel to the outermost layers of SN ejecta.

Numerical modeling of shock breakout in compact presupernovae of type Ib/c is done by two independent methods. Peak color temperatures reach 4 to 8 million K for explosion energies 1 to 9 foe with a typical local time-scale of 0.03 seconds. In the obverver's frame the burst is smeared over few seconds due to light travel time correction R/c.

Thermonuclear supernovae are valuable for cosmology but their physics is not yet fully understood. Modeling the development and propagation of nuclear flame is complicated by numerous instabilities. The predictions of supernova light curves still involve some simplifying assumptions, but one can use the comparison of the computed fluxes with observations to constrain the explosion mechanism. In spite of great progress in recent years, a number of issues remains unsolved both in flame physics and light curve modeling.

We present results of hydrodynamical simulations of young supernova remnants. To model the ejecta, we use several models (discussed in literature) of type Ia supernova explosions with different abundances. Our hydro models are one-dimensional and spherically symmetrical, but they take into account ionization kinetics with all important processes. We include detailed calculations for the X-ray emission, allowing for time-dependent ionization and recombination. In particular, we compare the computed X-ray spectra with recent XMM-Newton observations of the Tycho SN remnant. Our goal is to find the most viable thermonuclear SN model that gives good fits to both these X-ray observations and typical SN Ia light curves.

A treatment of line opacity in expanding medium is most crucial for the light curve (LC) modeling of Type Ia supernovae (SNe Ia). Spectral lines are the main source of opacity inside SN Ia ejecta from ultraviolet through infrared range. Here we focus on the mean opacity for the energy equation. We solve the Boltzmann equation for photons in the comoving frame for a spherically-symmetrical flow. For rectangle line profiles we find an analytical expression for frequency averaged intensity and absorptive opacity. The results differ from previously known heuristic solutions. The LCs in the I-band are in better agreement with observations.

Supernova 1987A remains the most well-studied supernova to date. Observations produced excellent broad-band photometric and spectroscopic coverage over a wide wavelength range at all epochs. We model the observed spectra from Day 1 to Day 81 using a hydrodynamical model. We show that good agreement can be obtained at times up to about 60 days, if we allow for extended nickel mixing. Later than about 60 days the observed Balmer lines become stronger than our models can reproduce. We show that this is likely due to a more complicated distribution of gamma-rays than we allow for in our spherically symmetric calculations. We present synthetic light curves in UBVRIJHK and a synthetic bolometric light curve. Using this broad baseline of detailed spectroscopic models we find a distance modulus mu = 18.5 +/- 0.2 using the SEAM method of determining distances to supernovae. We find that the explosion time agrees with that of the neutrino burst and is constrained at 68% confidence to within +/- 0.9 days. We argue that the weak Balmer lines of our detailed model calculations casts doubt on the accuracy of the purely photometric EPM method. We also suggest that Type IIP supernovae will be most useful as distance indicators at early times due to a variety of effects.

We define the critical coordinate velocity v_c. A particle moving radially in Schwarzschild background with this velocity, v_c= c/sqrt 3, is neither accelerated, nor decelerated if gravitational field is weak, r_g << r, where r_g is the gravitational radius, while r is the current one. We find that the numerical coincidence of v_c with velocity of sound in ultrarelativistic plasma, u_s, is accidental, since two velocities are different if the number of spatial dimensions is not equal to 3

Modelling of high resolution Balmer line profiles in the early-time spectra of SN 1998S shows that the inferred fast (roughly 400 km/s) circumstellar (CS) gas on days 23 and 42 post-explosion is confined to a narrow, negative velocity gradient shell just above the photosphere. This gas may be identified with a slow (v < 40 km/s) progenitor wind accelerated at the ejecta-wind interface. In this scenario, the photosphere coincides with a cool dense shell formed in the reverse shock. Acceleration by radiation from the supernova or by a shock-accelerated relativistic particle precursor are both possible explanations for the observed fast CS gas. An alternative, equally plausible scenario is that the fast CS gas is accelerated within shocked clouds engulfed by the outer shock, as it propagates through the intercloud wind.

A model of a carbon-oxygen (C--O) presupernova core with an initial mass 1.33 M_\odot, an initial carbon mass fraction 0.27, and with an average mass growth-rate 5 x 10^{-7} M_\odot/yr due to accretion in a binary system was evolved from initial central density 10^9 g/cm^3, and temperature 2.05 x 10^8 K through convective core formation and its subsequent expansion to the carbon runaway at the center. The only thermonuclear reaction contained in the equations of evolution and runaway was the carbon burning reaction 12C + 12C with an energy release corresponding to the full transition of carbon and oxygen (with the same rate as carbon) into 56Ni. As a parameter we take \alpha_c - a ratio of a mixing length to the size of the convective zone. In spite of the crude assumptions, we obtained a pattern of the runaway acceptable for the supernova theory with the strong dependence of its duration on \alpha_c. In the variants with large enough values of \alpha_c=4.0 x 10^{-3} and 3.0 x 10^{-3} the fuel combustion occurred from the very beginning as a prompt detonation. In the range of 2.0 x 10^{-3} >= \alpha_c >= 3.0 x 10^{-4} the burning started as a deflagration with excitation of stellar pulsations with growing amplitude. Eventually, the detonation set in, which was activated near the surface layers of the presupernova (with m about 1.33 M_\odot) and penetrated into the star down to the deflagration front. Excitation of model pulsations and formation of a detonation front are described in detail for the variant with \alpha_c=1.0 x 10^{-3}.

Supernova 1987A remains the most well-observed and well-studied supernova to date. Observations produced excellent broad-band photometric and spectroscopic coverage over a wide wavelength range at all epochs. Here, we focus on the very early spectroscopic observations. Only recently have numerical models been of sufficient detail to accurately explain the observed spectra. In SN 1987A, good agreement has been found between observed and synthetic spectra for day one, but by day four, the predicted Balmer lines become much weaker than the observed lines. We present the results of work based on a radiation-hydrodynamic model by Blinnikov and collaborators. Synthetic non-LTE spectra generated from this model by the general radiation transfer code PHOENIX strongly support the theory that significant mixing of nickel into the outer envelope is required to maintain strong Balmer lines. Preliminary results suggest a lower limit to the average nickel mass of 1.0 \times 10^{-5} solar masses is required above 5000 \kmps by day four. PHOENIX models thus have the potential to be a sensitive probe for nickel mixing in the outer layers of a supernova.

The data on the location of gamma-ray bursts (GRBs) relative to their host galaxies are used to derive the distribution of surface density of GRBs along the galaxy radius. It is shown that the gradient of GRB surface density changes abruptly near the half-light radius. In the central parts of galaxies the distribution of GRBs resembles closely the luminosity distribution, while in the outer parts the galactic surface brightness falls much steeper than the GRBs density. The radial distribution of type Ib/c supernovae is investigated on the basis of enlarged statistics. It is shown that SNe Ib/c do not differ significantly from other types of supernovae and their distribution is more similar to the one for recent star formation sites than that of GRBs. In spite of the poor statistics of GRBs, the difference in the distributions of active star formation regions and GRBs appears to be significant. We get the Kolmogorov-Smirnov probability P_ks of only 4% that GRBs and star-forming sites belong to the same distribution. The correlation of GRBs with the distribution of dark matter in the outer parts of galaxies is not excluded.

Recently, there appeared many papers devoted to the modeling of X-ray properties of Cygnus X-1 and other black hole accretion disk candidates: e.g. J.Poutanen, J.Krolik, F.Ryde, MNRAS 292 (1997) L21; E.Agol, J.Krolik, ApJ 507 (1998) 304; A.Beloborodov, in ``High Energy Processes in Accreting Black Holes'', ASP Conf. Series, 161, (1999), p.295. The goal of this electronic publication is to draw attention to our old papers where many ideas of recent discussions were anticipated (hot coronae, Comptonization, photon damping of waves, particle acceleration and matter ejection from accretion disks with large scale poloidal magnetic fields, etc.)

Ultraviolet light curves are calculated for several thermonuclear supernova models using a multifrequency radiation hydrodynamic code. It is found that Chandrasekhar-mass models produce very similar light curves both for detonation and deflagration. Sub-Chandrasekhar-mass models essentially differ from ``normal'' Chandrasekhar ones regarding behaviour of their UV fluxes. Differences in absolute brightness and in shape of light curves of thermonuclear supernovae could be detectable up to 300 Mpc with modern UV space telescopes.

The optical/UV light curves of SN 1987A are analyzed with the multi-energy group radiation hydrodynamics code STELLA. The calculated monochromatic and bolometric light curves are compared with observations shortly after shock breakout, during the early plateau, through the broad second maximum, and during the earliest phase of the radioactive tail. We have concentrated on a progenitor model calculated by Nomoto & Hashimoto and Saio, Nomoto, & Kato, which assumes that 14 solar masses of the stellar mass is ejected. Using this model, we have updated constraints on the explosion energy and the extent of mixing in the ejecta. In particular, we determine the most likely range of E/M (explosion energy over ejecta mass) and R_0 (radius of the progenitor). In general, our best models have energies in the range E = (1.1 +/- 0.3) x 10^{51} ergs, and the agreement is better than in earlier, flux-limited diffusion calculations for the same explosion energy. Our modeled B and V fluxes compare well with observations, while the flux in U undershoots after about 10 days by a factor of a few, presumably due to NLTE and line transfer effects. We also compare our results with IUE observations, and a very good quantitative agreement is found for the first days, and for one IUE band (2500-3000 A) as long as for 3 months. We point out that the V flux estimated by McNaught & Zoltowski should probably be revised to a lower value.

The properties of the cosmic Gamma-ray Bursts (GRBs) are briefly summarized. A detailed bibliography is given with titles of the papers. Two fundamental theoretical problems are pointed out: the problem of the energy source, and the problem of compactness. I demonstrate some inconsistencies in the estimates of the fireball optical thickness that are widely used in the discussion of the latter problem. The possible connection of GRBs with the Dark Matter candidates is mentioned. I argue that GRBs can be produced by collapses or mergers of stars made of one probable Dark Matter candidate, namely the mirror particles. I speculate on the impact that the parameters of the neutrino oscillations might have on the observed properties of GRBs if the latter are the products of mirror star deaths.

The calculations of the light curves of thermonuclear supernovae are carried out by a method of multi-group radiation hydrodynamics. The effects of spectral lines and expansion opacity are taken into account. The predictions for UBVI fluxes are given. The values of rise time for B and V bands found in our calculations are in good agreement with the observed values. We explain why our results for the rise time have more solid physical justification than those obtained by other authors. It is shown that small variations in the chemical composition of the ejecta, produced in the explosions with different regimes of nuclear burning, can influence drastically the light curve decline in the B band and, to a lesser extent, in the V band. We argue that recent results on positive cosmological constant Lambda, found from the high redshift supernova observations, could be wrong in the case of possible variations of the preferred mode of nuclear burning in the earlier Universe.

We argue that a nonthermally looking spectrum of a gamma-ray burst (GRB) can be formed as a superposition of a set of thermal blackbody spectra. This superposition may be done by time integration which is present even in `time resolved' GRB spectroscopy. A nonthermal spectrum can be obtained also by the space integration which should take place unless all the emission comes from a plane front moving exactly towards the observer. We propose a model of the gamma-ray burst spectrum formation based on this idea. This model allows the GRB radiation to be optically thick and to have higher values of baryon load. Thus the latter is limited by the energy considerations only, and not by the condition of a small optical depth.

I argue that cosmic Gamma-ray Bursts (GRB) may be produced by collapses or mergers of stars made of `mirror' matter. The mirror neutrinos (which are sterile for our matter) produced at these events can oscillate into ordinary neutrinos. The annihilations or decays of the latter create an electron-positron plasma and subsequent relativistic fireball with a very low baryon load needed for GRBs. The concept of mirror matter is able to explain several key problems of modern astrophysics: neutrino anomalies, the missing mass, MACHO microlensing events and GRBs. Thus this concept becomes very appealing and should be considered quite seriously and attentively.

A hydro code coupled with radiation transfer was applied to produce monochromatic light curves of two models of type II supernovae (SN II) simulating SN II-P and SN IIb (SN 1993J-like). We then used these template light curves to evaluate the possibility of detecting SNe II at different redshifts. With a 5 hour exposure at VLT/FORS the SN II-P model may be detected at z=1. However, since our model of SN II-P is underluminous at early phase (t < 10 days) by roughly 1.5 mag a detection at z=2 is quite plausible. SN IIb can be detected as far as at z=4. For 100% detection efficiency up to z=2 one expects to find roughly 1 SN II yr^{-1} arcmin^{-2}.

The huge neutrino pulse that occurs during the collapse of a massive stellar core, is expected to contribute to the origination of a number of isotopes both of light chemical elements and heavy ones. It is shown that, in general, the heating of stellar matter due to the neutrino scattering off electrons and the heat released from the neutrino-helium breakup followed by the thermonuclear reactions should be taken into account. On the base of kinetic network, using all the important reactions up to Z=8, the main features and the time-dependent character of the neutrino- driven neutron flux are investigated. The time-dependent densities of free neutrons produced in helium breakup, Y_n(t), were used to calculate the r-process nucleosynthesis with another full kinetic network for 3200 nuclides. It was found that in the case of metal-deficient stars, Z < 0.01 Z(solar), the resulting density of free neutrons seems to be high enough to drive the r-process efficiently under favorable conditions. But it is impossible to obtain a sufficient amount of heavy nuclei in neutrino-induced r-process in a helium shell at radii R > R_cr \approx 10^9 cm. We speculate that to make the neutrino-induced r-process work efficiently in the shell, one has to invoke nonstandard presupernova models in which helium hopefully is closer to the collapsed core owing, for instance, to a large scale mixing or/and rotation and magnetic fields. Apart from this exotic possibility, the neutrino-induced nucleosynthesis in the helium shell is certainly not strong enough to explain the observed solar r-process abundances.

High quality spectra to be produced by SUBARU 8 meter telescope will be extremely valuable for numerous astronomical projects. Here I concentrate on two classical types of weak broad line objects -- distant supernovae and white dwarfs -- with emphasis of their importance for fundamental physics, not only for astronomy. With regard to recent results of MACHO experiment and HST weak star counts I discuss one exotic possibility: invisible stars made of the so called ``mirror'' matter.

Various types of expansions in series of Chebyshev-Hermite polynomials currently used in astrophysics for weakly non-normal distributions are compared, namely the Gram-Charlier, Gauss-Hermite and Edgeworth expansions. It is shown that the Gram-Charlier series is most suspect because of its poor convergence properties. The Gauss-Hermite expansion is better but it has no intrinsic measure of accuracy. The best results are achieved with the asymptotic Edgeworth expansion. We draw attention to the form of this expansion found by Petrov for arbitrary order of the asymptotic parameter and present a simple algorithm realizing Petrov's prescription for the Edgeworth expansion. The results are illustrated by examples similar to the problems arising when fitting spectral line profiles of galaxies, supernovae, or other stars, and for the case of approximating the probability distribution of peculiar velocities in the cosmic string model of structure formation.

The light curve of Supernova (SN) 1993J is calculated using two approaches to radiation transport as exemplified by the two computer codes, STELLA and EDDINGTON. Particular attention is paid to shock breakout and the photometry in the U, B, and V bands during the first 120 days. The hydrodynamical model, the explosion of a 13 Msun star which had lost most of its hydrogenic envelope to a companion, is the same in each calculation. The comparison elucidates differences between the approaches and also serves to validate the results of both. STELLA includes implicit hydrodynamics and is able to model supernova evolution at early times, before the expansion is homologous. STELLA also employs multi-group photonics and is able to follow the radiation as it decouples from the matter. EDDINGTON uses a different algorithm for integrating the transport equation, assumes homologous expansion, and uses a finer frequency resolution. Good agreement is achieved between the two codes only when compatible physical assumptions are made about the opacity. A new result for SN 1993J is a prediction of the continuum spectrum near the shock breakout (calculated by STELLA) which is superior to the results of other standard single energy group hydrocodes such as VISPHOT or TITAN. Based on the results of our independent codes, we discuss the uncertainties involved in the current time dependent models of supernova light curves.

An energy deposition of \sim 10^{50} ergs into the exterior 10^{-3} M_\sun layers of a red giant is calculated to produce an optical phenomenon similar to afterglows of gamma-ray bursts (GRB) recently observed. This model can be realized if a GRB is generated by some mechanism in a close binary system. In contrast to a ``hypernova'' scenario for GRB recently proposed by Paczy\'nski (1997), this model does not require a huge kinetic energy of the expanding shell to explain optical afterglows of GRB.