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.

Models of pair-instability supernovae (PISNe) predict a gap in black hole (BH) masses between $\sim 45M_\odot-120M_\odot$, which is referred to as the upper BH mass-gap. With the advent of gravitational-wave astrophysics it has become possible to test this prediction, and there is an important associated effort to understand what theoretical uncertainties modify the boundaries of this gap. In this work we study the impact of rotation on the hydrodynamics of PISNe, which leave no compact remnant, as well as the evolution of pulsational-PISNe (PPISNe), which undergo thermonuclear eruptions before forming a compact object. We perform simulations of non-rotating and rapidly-rotating stripped helium stars in a metal poor environment $(Z_\odot/50)$ in order to resolve the lower edge of the upper mass-gap. We find that the outcome of our simulations is dependent on the efficiency of angular momentum transport, with models that include efficient coupling through the Spruit-Tayler dynamo shifting the lower edge of the mass-gap upwards by $\sim 4\%$, while simulations that do not include this effect shift it upwards by $\sim 15\%$. From this, we expect the lower edge of the upper mass-gap to be dependent on BH spin, which can be tested as the number of observed BH mergers increases. Moreover, we show that stars undergoing PPISNe have extended envelopes ($R\sim 10-1000~R_\odot$) at iron-core collapse, making them promising progenitors for ultra-long 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 present optical and near-infrared photometry and spectroscopy of the Type IIn supernova (SN) 2014ab, obtained by the Carnegie Supernova Project II (CSP-II) and initiated immediately after its optical discovery. We also present mid-infrared photometry obtained by the Wide-field Infrared Survey Explorer (WISE) satellite extending from 56 days prior to the optical discovery to over 1600 days. The light curve of SN 2014ab evolves slowly, while the spectra exhibit strong emission features produced from the interaction between rapidly expanding ejecta and dense circumstellar matter. The light curve and spectral properties are very similar to those of SN 2010jl. The estimated mass-loss rate of the progenitor of SN 2014ab is of the order of 0.1 Msun/yr under the assumption of spherically symmetric circumstellar matter and steady mass loss. Although the mid-infrared luminosity increases due to emission from dust, which is characterized by a blackbody temperature close to the dust evaporation temperature (~ 2000 K), no clear signatures of in situ dust formation within the cold dense shell located behind the forward shock are observed in SN 2014ab in early phases. Mid-infrared emission of SN 2014ab may originate from pre-existing dust located within dense circumstellar matter that is heated by the SN shock or shock-driven radiation. Finally, for the benefit of the community, we also present in an Appendix five near-infrared spectra of SN 2010jl obtained between 450 to 1300 days post discovery.

A number of astrophyical transients originating from stellar explosions are powered by the collision of the ejected material with the circumstellar medium, which efficiently produces thermal radiation via shock dissipation. We investigate how such interaction-powered transients are distributed in the peak bolometric luminosity vs the rise time phase space. Taking the advantage of less time-consuming one-dimensional simulations with spherical symmetry, we calculated more than 500 models with different circumstellar mass and radius, ejecta mass and energy, and chemical compositions. The peak bolometric luminosity, the total radiated energy, and the rise time of the interaction-powered emission are measured for each simulated light curve. We consider how these characteristic quantities are determined as a function of the model parameters and discuss possible implications for the observed populations of (potential) interaction-powered transients, such as type IIn supernovae and fast blue optical transients.

We present a sample of supernovae Type IIn (SNe IIn) from the untargeted, magnitude-limited surveys of the Palomar Transient Factory (PTF) and its successor, the intermediate PTF (iPTF). The SNe IIn found and followed by the PTF/iPTF were used to select a sample of 42 events with useful constraints on the rise times as well as with available post-peak photometry. The sample SNe were discovered in 2009-2016 and have at least one low-resolution classification spectrum, as well as photometry from the P48 and P60 telescopes at Palomar Observatory. We study the light-curve properties of these SNe IIn using spline fits (for the peak and the declining portion) and template matching (for the rising portion). We find that the typical rise times are divided into fast and slow risers at $20\pm6$ d and $50\pm11$ d, respectively. The decline rates are possibly divided into two clusters, but this division has weak statistical significance. We find no significant correlation between the peak luminosity of SNe IIn and their rise times, but the more luminous SNe IIn are generally found to be more long-lasting. Slowly rising SNe IIn are generally found to decline slowly. The SNe in our sample were hosted by galaxies of absolute magnitude $-22 \lesssim M_g \lesssim -13$ mag. The K-corrections at light-curve peak of the SNe IIn in our sample are found to be within 0.2 mag for the observer's frame $r$-band, for SNe at redshifts $z < 0.25$. By applying K-corrections and also including ostensibly "superluminous" SNe IIn, we find that the peak magnitudes are $M_{\rm peak}^{r} = -19.18\pm1.32$ mag. We conclude that the occurrence of conspicuous light-curve bumps in SNe IIn, such as in iPTF13z, are limited to $1.4^{+14.6}_{-1.0} \%$ of the SNe IIn. We also investigate a possible sub-type of SNe IIn with a fast rise to a $\gtrsim 50$ d plateau followed by a slow, linear decline.

We present optical and near-infrared observations of SN~Ib~2019ehk. We show that it evolved to a Ca-rich transient according to its spectral properties and evolution in late phases. It, however, shows a few distinguishable properties from the canonical Ca-rich transients: a short-duration first peak in the light curve, high peak luminosity, and association with a star-forming environment. Indeed, some of these features are shared with iPTF14gqr and iPTF16hgs, which are candidates for a special class of core-collapse SNe (CCSNe): the so-called ultra-stripped envelope SNe, i.e., a relatively low-mass He (or C+O) star explosion in a binary as a precursor of double neutron star binaries. The estimated ejecta mass ($0.43 M_\odot$) and explosion energy ($1.7 \times 10^{50} $~erg) are consistent with this scenario. The analysis of the first peak suggests existence of dense circumstellar material in the vicinity of the progenitor, implying a CCSN origin. Based on these analyses, we suggest SN 2019ehk is another candidate for an ultra-stripped envelope SN. These ultra-stripped envelope SN candidates seem to form a subpopulation among Ca-rich transients, associated with young population. We propose that the key to distinguishing this population is the early first peak in their light curves.

We present multi-wavelength observations of two gap transients followed by the Carnegie Supernova Project-II and supplemented with data obtained by a number of different programs. Here in the first of two papers, we focus on the intermediate luminosity red transient (ILRT) designated SNhunt120, while in a companion paper we examine the luminous red novae AT 2014ej. Our data set for SNhunt120 consists of an early optical discovery, estimated to be within 3 days after outburst, the subsequent optical and near-infrared broadband followup extending over a $\sim$2 month period, two visual- and two near-infrared wavelength spectra, and Spitzer Space Telescope observations extending from early ($+$28 d) to late ($+$1155 d) phases. SNhunt120 resembles other ILRTs such as NGC 300-2008-OT and SN 2008S, and like these other ILRTs, SNhunt120 exhibits prevalent mid-infrared emission at both early and late phases. From the comparison of SNhunt120 and other ILRTs to electron-capture supernova simulations, we find that the current models underestimate the explosion kinetic energy and thereby produce synthetic light curves that over-estimate the luminosity. Finally, examination of pre-outburst Hubble Space Telescope images yields no progenitor detection.

We present optical and near-infrared broadband photometry and optical spectra of AT 2014ej from the the Carnegie Supernova Project-II. These observations are complemented with data from the CHilean Automatic Supernova sEarch, the Public ESO Spectroscopic Survey of Transient Objects, and from the Backyard Observatory Supernova Search. Observational signatures of AT 2014ej reveal that it is similar to other members of the gap-transient subclass known as luminous red novae (LRNe), including the ubiquitous double hump light curve and spectral properties akin to the LRN SN 2017jfs. A medium-dispersion, visual-wavelength spectrum of AT 2014ej taken the Magellan Clay telescope, exhibits a P Cygni H$\alpha$ feature characterized by a blue velocity at zero intensity of $\approx 110$ km s$^{-1}$ and a P Cygni minimum velocity of $\approx70$ km s$^{-1}$, and which we attribute to emission from a circumstellar wind. Inspection of pre-outbust Hubble Space Telescope images yields no conclusive progenitor detection. In comparison with a sample of LRNe from the literature, AT 2014ej lies among the brighter end of the luminosity distribution. Comparison of the ultra-violet, optical, infrared (UVOIR) light curves of well-observed LRNe to common-envelope evolution models from the literature, indicates the models under predict the luminosity of the comparison sample at all phases and also produce inconsistent time-scales of the secondary peak. Future efforts to model LRNe should expand upon the current parameter space explored and therefore may consider more massive systems and a wider range of dynamical timescales.

We present optical and near-IR data of three Type II supernovae (SNe II), SN 2008bm, SN 2009aj, and SN 2009au. These SNe display the following common characteristics: signs of early interaction of the ejecta with circumstellar material (CSM), blue $B\!-\!V$ colours, weakness of metal lines, low expansion velocities, and $V$-band absolute magnitudes 2-3 mag brighter than those expected for normal SNe II based on their expansion velocities. Two more SNe reported in the literature (SN 1983K and LSQ13fn) share properties similar to our sample. Analysing this set of five SNe II, which are luminous for their low expansion velocities (LLEV), we find that their properties can be reproduced assuming ejecta-CSM interaction that lasts between 4-11 weeks post explosion. The contribution of this interaction to the radiation field seems to be the dominant component determining the observed weakness of metal lines in the spectra rather than the progenitor metallicity. Based on hydrodynamic simulations, we find that the interaction of the ejecta with a CSM of ~3.6 M$_\odot$ can reproduce the light curves and expansion velocities of SN 2009aj. Using data collected by the Chilean Automatic Supernova Search, we estimate an upper limit for the LLEV SNe II fraction to be 2-4 per cent of all normal SNe II. With the current data-set, it is not clear whether the LLEV events are a separated class of SNe II with a different progenitor system, or whether they are the extreme of a continuum mediated by CSM interaction with the rest of the normal SN II population.

We present our spectropolarimetric observations of SN 2017egm, a Type I superluminous supernova (SLSN-I) in a nearby galaxy NGC 3191, with Subaru telescope at 185.0 days after the g-band maximum light. This is the first spectropolarimetric observation for SLSNe at late phases. We find that the degree of the polarization in the late phase significantly changes from that measured at the earlier phase. The spectrum at the late phase shows a strong Ca emission line and therefore we reliably estimate the interstellar polarization component assuming that the emission line is intrinsically unpolarized. By subtracting the estimated interstellar polarization, we find that the intrinsic polarization at the early phase is only ~ 0.2 %, which indicates an almost spherical photosphere, with an axial ratio ~ 1.05. The intrinsic polarization at the late phase increases to ~ 0.8 %, which corresponds to the photosphere with an axial ratio 1.2. A nearly constant position angle of the polarization suggests the inner ejecta are almost axisymmetric. By these observations, we conclude that the inner ejecta are more aspherical than the outer ejecta. It may suggest the presence of a central energy source producing aspherical inner ejecta.

We present ultra-violet to mid-infrared observations of the long-lasting Type IIn supernova (SN) 2013L obtained by the Carnegie Supernova Project II (CSP-II). The spectra of SN 2013L are dominated by H emission features characterized by three components attributed to different regions. A unique feature of this Type IIn SN is that the blue shifted line profile is dominated by the macroscopic velocity of the expanding shock wave of the SN. We are therefore able to trace the evolution of the shock velocity in the dense and partially opaque circumstellar medium (CSM), from $\sim 4800~km~s^{-1}$ at +48 d, decreasing as $t^{-0.23}$ to $\sim 2700~km~s^{-1}$ after a year. We perform spectral modeling of both the broad- and intermediate-velocity components of the H$\alpha$ line profile. The high-velocity component is consistent with emission from a radially thin, spherical shell located behind the expanding shock with emission wings broadened by electron scattering. We propose that the intermediate component originates from pre-ionized gas from the unshocked dense CSM with the same velocity as the narrow component, $\sim 100~km~s^{-1}$, but also broadened by electron scattering. The spectral energy distributions (SEDs) of SN 2013L after +132 d are well reproduced by a two-component black-body (BB) model. The circumstellar-interaction model of the bolometric light curve reveals a mass-loss rate history with large values ($1.7\times 10^{-2} - 0.15~M_\odot~yr^{-1}$) over the $\sim $25 - 40 years before explosion. The drop in the light curve at $\sim 350$ days and presence of electron scattering wings at late epochs indicate an anisotropic CSM. The mass-loss rate values and the unshocked CSM velocity are consistent with the characteristics of a massive star, such as a luminous blue variable (LBV) undergoing strong eruptions, similar to $\eta$ Carina.

With a booming number of Type Ia supernovae (SNe Ia) discovered within a few days of their explosions, a fraction of SNe Ia that show luminosity excess in the early phase (early-excess SNe Ia) have been confirmed. In this article, we report early-phase observations of seven photometrically normal SNe Ia (six early detections and one deep non-detection limit) at the COSMOS field through a half-year transient survey as a part of the Hyper Suprime-Cam Subaru Strategic Program (HSC SSP). In particular, a blue light-curve excess was discovered for HSC17bmhk, a normal SN Ia with rise time longer than 18.8 days, during the first four days after the discovery. The blue early excess in optical wavelength can be explained not only by interactions with a non-degenerate companion or surrounding dense circumstellar matter but also radiation powered by radioactive decays of $^{56}$Ni at the surface of the SN ejecta. Given the growing evidence of the early-excess discoveries in normal SNe Ia that have longer rise times than the average and a similarity in the nature of the blue excess to a luminous SN Ia subclass, we infer that early excess discovered in HSC17bmhk and other normal SNe Ia are most likely attributed to radioactive $^{56}$Ni decay at the surface of the SN ejecta. In order to successfully identify normal SNe Ia with early excess similar to that of HSC17bmhk, early UV photometries or high-cadence blue-band surveys are necessary.

Rapidly evolving transients form a new class of transients which show shorter timescales of the light curves than those of typical core-collapse and thermonuclear supernovae. We performed a systematic search for rapidly evolving transients using the deep data taken with the Hyper Suprime-Cam Subaru Strategic Program Transient Survey. By measuring the timescales of the light curves of 1824 transients, we identified 5 rapidly evolving transients. Our samples are found in a wide range of redshifts (0.3 $\le$ z $\le$ 1.5) and peak absolute magnitudes ($-$17 $\ge$ $M_i$ $\ge$ $-$20). The properties of the light curves are similar to those of the previously discovered rapidly evolving transients. They show a relatively blue spectral energy distribution, with the best-fit blackbody of 8,000 - 18,000 K. We show that some of the transients require power sources other than the radioactive decays of $^{56}$Ni because of their high peak luminosities and short timescales. The host galaxies of all the samples are star-forming galaxies, suggesting a massive star origin for the rapidly evolving transients. The event rate is roughly estimated to be $\sim$4,000 events yr$^{-1}$ Gpc$^{-3}$, which is about 1 $\%$ of core-collapse supernovae.

During the first three observing runs of the Advanced gravitational-wave detector network, the LIGO/Virgo collaboration detected several black hole binary (BHBH) mergers. As the population of detected BHBH mergers grows, it will become possible to constrain different channels for their formation. Here we consider the chemically homogeneous evolution (CHE) channel in close binaries, by performing population synthesis simulations that combine realistic binary models with detailed cosmological calculations of the chemical and star-formation history of the Universe. This allows us to constrain population properties, as well as cosmological and aLIGO detection rates of BHBH mergers formed through this pathway. We predict a BHBH merger rate at redshift zero of $5.8 \hspace{1mm} \textrm{Gpc}^{-3} \textrm{yr}^{-1}$ through the CHE channel, to be compared with aLIGO's measured rate of ${53.2}_{-28.2}^{+55.8} \hspace{1mm} \text{Gpc}^{-3}\text{yr}^{-1}$, and find that eventual merger systems have BH masses in the range $17 - 43 \hspace{1mm} \textrm{M}_{\odot}$ below the pair-instability supernova (PISN) gap, and $>124 \hspace{1mm} \textrm{M}_{\odot}$ above the PISN gap. We further investigate the effects of momentum kicks during black hole formation, calculate cosmological and magnitude limited PISN rates and investigate the effects of high-redshift deviations in the star formation rate. We find that momentum kicks tend to increase delay times of BHBH systems, and our magnitude limited PISN rate estimates indicate that current deep surveys should be able to detect such events. Lastly, we find that our cosmological merger rate estimates change by at most $\sim 8\%$ for mild deviations of the star formation rate in the early Universe, and by up to $\sim 40\%$ for extreme deviations.

Searches for optical transients are usually performed with a cadence of days to weeks, optimised for supernova discovery. The optical fast transient sky is still largely unexplored, with only a few surveys to date having placed meaningful constraints on the detection of extragalactic transients evolving at sub-hour timescales. Here, we present the results of deep searches for dim, minute-timescale extragalactic fast transients using the Dark Energy Camera, a core facility of our all-wavelength and all-messenger Deeper, Wider, Faster programme. We used continuous 20s exposures to systematically probe timescales down to 1.17 minutes at magnitude limits $g > 23$ (AB), detecting hundreds of transient and variable sources. Nine candidates passed our strict criteria on duration and non-stellarity, all of which could be classified as flare stars based on deep multi-band imaging. Searches for fast radio burst and gamma-ray counterparts during simultaneous multi-facility observations yielded no counterparts to the optical transients. Also, no long-term variability was detected with pre-imaging and follow-up observations using the SkyMapper optical telescope. We place upper limits for minute-timescale fast optical transient rates for a range of depths and timescales. Finally, we demonstrate that optical $g$-band light curve behaviour alone cannot discriminate between confirmed extragalactic fast transients such as prompt GRB flashes and Galactic stellar flares.

We perform a series of two-dimensional radiation-hydrodynamic simulations of the collision between supernova ejecta and circumstellar media (CSM). The hydrodynamic interaction of a fast flow and the surrounding media efficiently dissipates the kinetic energy of the fast flow and considered as a dominant energy source for a specific class of core-collapse supernovae. Despite some observational evidence for aspherical ejecta and/or CSM structure, multi-dimensional effects in the ejecta-CSM interaction are relatively unexplored. Our numerical simulations equipped with an adaptive mesh refinement technique successfully reproduce hydrodynamic instabilities developing around the ejecta-CSM interface. We also investigate effects of disk-like CSM on the dynamical evolution of supernova ejecta and bolometric light curves. We find that emission powered by ejecta-disk interaction exhibits significant viewing angle dependence. For a line of sight close to the symmetry axis, the observer directly sees the supernova ejecta, leading to a short brightening timescale. For an observer seeing the emission through the CSM disk, thermal photons diffuse throughout the CSM and thus the light curve is severely smeared out.

The origin of iPTF14hls, which had Type IIP supernova-like spectra but kept bright for almost two years with little spectral evolution, is still unclear. We here propose that iPTF14hls was not a sudden outburst like supernovae but rather a long-term outflow similar to stellar winds. The properties of iPTF14hls, which are at odds with a supernova scenario, become natural when interpreted as a stellar wind with variable mass-loss rate. Based on the wind hypothesis, we estimate the mass-loss rates of iPTF14hls in the bright phase. We find that the instantaneous mass-loss rate of iPTF14hls during the 2-year bright phase was more than a few Msun/yr ("hyper-wind") and it reached as much as 10 Msun/yr. The total mass lost over two years was about 10 Msun. Interestingly, we find that the light curve of iPTF14hls has a very similar shape to that of eta Carinae during the Great Eruption, which also experienced a similar but less extreme brightening accompanied by extraordinary mass loss, shedding more than 10 Msun in 10 years. The progenitor of iPTF14hls is less than 150 Msun if it still exists, which is similar to eta Carinae. The two phenomena may be related to a continuum-driven extreme wind from very massive stars.

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.

We investigate a possibility that a recently reported radio transient in M81, VTC J095517.5+690813, is caused by an accretion-induced collapse of a white dwarf. It became bright in radio but no associated optical transient was discovered. The accretion-induced collapse is predicted to be radio bright but optically faint, satisfying the observed property. We compare predicted radio emission from the accretion-induced collapse with that of VTC J095517.5+690813 and show that it can be reasonably explained by the accretion-induced collapse. Although it is difficult to firmly conclude that VTC J095517.5+690813 is an accretion-induced collapse, our study shows that radio-bright transients without an optical counterpart could still be related to stellar deaths.

We investigate predicted circumstellar properties of Type Ia supernova progenitor systems with non-degenerate helium star donors. It has been suggested that systems consisting of a carbon+oxygen white dwarf and a helium star can lead to Type Ia supernova explosions. Binary evolution calculations for the helium star donor channel predict that such a progenitor system is in either a stable helium-shell burning phase or a weak helium-shell flash phase at the time of the Type Ia supernova explosion. By taking the binary evolution models from our previous study, we show that a large fraction of the progenitor systems with a helium star donor have low enough density to explain the current non-detection of radio emission from Type Ia supernovae. Most of the progenitor systems in the weak helium-shell flash phase at the time of the Type Ia supernova explosions, which may dominate the prompt (short delay time) Type Ia supernova population, have both low circumstellar density and a faint helium star donor to account for the non-detection of radio emission and a pre-explosion companion star in SN 2011fe and SN 2014J. We also find some progenitor systems that are consistent with the properties of the companion star candidate identified at the explosion location of Type Iax SN 2012Z.

The optical and ultra-violet broadband photometric and spectroscopic observations of the Type II supernova (SN) 2016gfy are presented. The $V$-band light curve (LC) shows a distinct plateau phase with a slope, $s_2$ $\sim$ 0.12 mag / 100 d and a duration of 90 $\pm$ 5 d. Detailed analysis of SN 2016gfy provided a mean $\rm^{56}Ni$ mass of 0.033 $\pm$ 0.003 $\rm M_{\odot}$, a progenitor radius of $\sim$ 350-700 $\rm R_{\odot}$, a progenitor mass of $\sim$ 12-15 $\rm M_{\odot}$ and an explosion energy of 0.9-1.4$\rm \times 10^{51}\ erg\ s^{-1}$. The P-Cygni profile of H$\rm \alpha$ in the early phase spectra ($\sim$ 11-21 d) shows a boxy emission. Assuming that this profile arises from the interaction of the SN ejecta with the pre-existing circumstellar material (CSM), it is inferred that the progenitor underwent a recent episode (30-80 years prior to the explosion) of enhanced mass loss. Numerical modeling suggests that the early LC peak is reproduced better with an existing CSM of 0.15 $\rm M_{\odot}$ spread out to $\sim$ 70 AU. A late-plateau bump is seen in the $VRI$ LCs during $\sim$ 50-95 d. This bump is explained as a result of the CSM interaction and/or partial mixing of radioactive $\rm ^{56}Ni$ in the SN ejecta. Using strong-line diagnostics, a sub-solar oxygen abundance is estimated for the supernova H II region (12 + log(O/H) = 8.50 $\pm$ 0.11), indicating an average metallicity for the host of a Type II SN. A star formation rate of $\sim$ 8.5 $\rm M_{\odot}\ yr^{-1}$ is estimated for NGC 2276 using the archival $GALEX$ FUV data.

We report our observations of HSC16aayt (SN 2016jiu), which was discovered by the Subaru/Hyper Suprime-Cam (HSC) transient survey conducted as part of Subaru Strategic Program (SSP). It shows very slow photometric evolution and its rise time is more than 100 days. The optical magnitude change in 400 days remains within 0.6 mag. Spectra of HSC16aayt show a strong narrow emission line and we classify it as a Type IIn supernova. The redshift of HSC16aayt is 0.6814 +/- 0.0002 from the spectra. Its host galaxy center is at 5 kpc from the supernova location and HSC16aayt might be another example of isolated Type IIn supernovae, although the possible existence of underlying star forming activity of the host galaxy at the supernova location is not excluded.

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 an overview of a deep transient survey of the COSMOS field with the Subaru Hyper Suprime-Cam (HSC). The survey was performed for the 1.77 deg$^2$ ultra-deep layer and 5.78 deg$^2$ deep layer in the Subaru Strategic Program over 6- and 4-month periods from 2016 to 2017, respectively. The ultra-deep layer shows a median depth per epoch of 26.4, 26.3, 26.0, 25.6, and 24.6 mag in $g$, $r$, $i$, $z$, and $y$ bands, respectively; the deep layer is $\sim0.6$ mag shallower. In total, 1,824 supernova candidates were identified. Based on light curve fitting and derived light curve shape parameter, we classified 433 objects as Type Ia supernovae (SNe); among these candidates, 129 objects have spectroscopic or COSMOS2015 photometric redshifts and 58 objects are located at $z > 1$. Our unique dataset doubles the number of Type Ia SNe at $z > 1$ and enables various time-domain analyses of Type II SNe, high redshift superluminous SNe, variable stars, and active galactic nuclei.

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.

As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations.

Superluminous supernovae have been proposed to arise from Population III progenitors that explode as pair-instability supernovae. Pop III stars are the first generation of stars in the Universe, and are thought to form as late as $z \sim 6$. Future near-infrared imaging facilities such as ULTIMATE-Subaru can potentially detect and identify these PISNe with a dedicated survey. Gravitational lensing by intervening structure in the Universe can aid in the detection of these rare objects by magnifying the high-$z$ source population into detectability. We perform a mock survey with ULTIMATE-Subaru, taking into account lensing by line-of-sight structure to evaluate its impact on the predicted detection rate. We compare a LOS mass reconstruction using observational data from the Hyper Suprime Cam survey to results from cosmological simulations to test their consistency in calculating the magnification distribution in the Universe to high-$z$, but find that the data-based method is still limited by an inability to accurately characterize structure beyond $z \sim1.2$. We also evaluate a survey strategy of targeting massive galaxy clusters to take advantage of their large areas of high magnification. We find that targeting clusters can result in a gain of a factor of $\sim$two in the predicted number of detected PISNe at $z > 5$, and even higher gains with increasing redshift, given our assumed survey parameters. For the highest-redshift sources at $z \sim 7-9$, blank field surveys will not detect any sources, and lensing magnification by massive clusters will be necessary to observe this population.

ULTIMATE-Subaru (Ultra-wide Laser Tomographic Imager and MOS with AO for Transcendent Exploration on Subaru) and WFIRST (Wide Field Infra-Red Survey Telescope) are the next generation near-infrared instruments that have a large field-of-view. They allow us to conduct deep and wide transient surveys in near-infrared. Such a near-infrared transient survey enables us to find very distant supernovae that are redshifted to the near-infrared wavelengths. We have performed the mock transient surveys with ULTIMATE-Subaru and WFIRST to investigate their ability to discover Population III pair-instability supernovae. We found that a 5-year 1 deg2 K-band transient survey with the point-source limiting magnitude of 26.5 mag with ULTIMATE-Subaru may find about 2 Population III pair-instability supernovae beyond the redshift of 6. A 5-year 10 deg2 survey with WFIRST reaching 26.5 mag in the F184 band may find about 7 Population III pair-instability supernovae beyond the redshift of 6. We also find that the expected numbers of the Population III pair-instability supernova detections increase about a factor of 2 if the near-infrared transient surveys are performed towards clusters of galaxies. Other supernovae such as Population II pair-instability supernovae would also be detected in the same survey. This study demonstrates that the future wide-field near-infrared instruments allow us to investigate the explosions of the first generation supernovae by performing the deep and wide near-infrared transient surveys.

We present Keck spectroscopic observations of three probable high redshift superluminous supernovae (SLSNe) from the Subaru HIgh-Z sUpernova CAmpaign (SHIZUCA), confirming redshifts of 1.851, 1.965 and 2.399. The host galaxies were selected for transient monitoring from multi-band photometric redshifts. The supernovae are detected during their rise, and the classically scheduled spectra are collected near maximum light. The restframe far-ultraviolet (FUV;$\sim$1000A-2500A) spectra include a significant host galaxy flux contribution and we compare our host galaxy subtracted spectra to UV-luminous SNe from the literature. While the signal-to-noise ratios of the spectra presented here are sufficient for redshift confirmation, supernova spectroscopic type confirmation remains inconclusive. The success of the first SHIZUCA Keck spectroscopic follow-up program demonstrates that campaigns such as SHIZUCA are capable of identifying high redshift SLSNe with sufficient accuracy, speed and depth for rapid, well-cadenced and informative follow-up.

SN 2007bi and SN 1999as are among the first superluminous supernovae discovered. SN 2007bi was suggested to be powered by the radioactive decay of a large amount (5 - 10 Msun) of 56Ni. SN 1999as has a similar spectrum to SN 2007bi. One suggested way to synthesize such a large amount of 56Ni is through energetic core-collapse supernovae from very massive progenitors. Although the synthetic light curves of extremely energetic core-collapse supernovae have been shown to be consistent with SN 2007bi, no synthetic spectra have been reported. Here, we present synthetic spectra of extremely energetic core-collapse supernovae during the photospheric phases. We find that the ejecta density structure above 13,000 - 16,000 km/s needs to be cut in order to explain the co-existing broad and narrow line absorptions in SN 2007bi and SN 1999as. The density cut is likely caused by the interaction between the supernova ejecta and a dense circumstellar medium. Our results indicate that about 3 Msun of hydrogen-free dense circumstellar media might exist near the progenitors of SN 2007bi and SN 1999as. These massive circumstellar media would significantly affect the light-curve and spectral properties of the supernovae. The precursors that are sometimes observed in superluminous supernovae might be related to the collision of the ejecta with such dense circumstellar media. We also confirm results of previous studies that synthetic spectra from pair-instability supernova models do not match the early spectra of SN 2007bi and SN 1999as.

We report our first discoveries of high-redshift supernovae from the Subaru HIgh-Z sUpernova CAmpaign (SHIZUCA), the transient survey using Subaru/Hyper Suprime-Cam. We report the discovery of three supernovae at the spectroscopically-confirmed redshifts of 2.399 (HSC16adga), 1.965 (HSC17auzg), and 1.851 (HSC17dbpf), and two supernova candidates with the host-galaxy photometric redshifts of 3.2 (HSC16apuo) and 4.2 (HSC17dsid), respectively. In this paper, we present their photometric properties and the spectroscopic properties of the confirmed high-redshift supernovae are presented in the accompanying paper (Curtin et al. 2019). The supernovae with the confirmed redshifts of z ~ 2 have the rest ultraviolet peak magnitudes close to -21 mag and they are likely superluminous supernovae. The discovery of three supernovae at z ~ 2 roughly corresponds to the approximate event rate of ~ 900 +/- 520 Gpc-3 yr-1 with Poisson error, which is consistent with the total superluminous supernova rate estimated by extrapolating the local rate based on the cosmic star-formation history. Adding unconfirmed superluminous supernova candidates would increase the event rate. Our superluminous supernova candidates at the redshifts of around 3 and 4 indicate the approximate superluminous supernova rates of ~ 400 +/- 400 Gpc-3 yr-1 (z ~ 3) and ~ 500 +/- 500 Gpc-3 yr-1 (z ~ 4) with Poisson errors. Our initial results demonstrate the outstanding capability of Hyper Suprime-Cam to discover high-redshift supernovae.

A group of super-luminous supernovae (SL-SNe) characterised by broad light curves have been suggested to be Pair Instability SNe (PISNe). Nebular spectra computed using PISN models have failed to reproduce the broad emission lines observed in these SNe, casting doubts on their true nature. Here, models of both PISNe and the explosion following the collapse of the core of a very massive star (100 Msun) are used to compute nebular spectra, which are compared to the spectrum of the prototypical PISN candidate, SN 2007bi. PISN models are confirmed to produce synthetic spectra showing narrow emission lines, resulting from the confinement of 56Ni to the lowest velocities (<~ 2000 km/s) and in clear disagreement with the spectrum of SN 2007bi. Spectra more closely resembling SN 2007bi are obtained if the PISN models are fully mixed in abundance. Massive core-collapse models produce enough 56Ni to power the light curve of PISN candidates, but their spectra are also not adequate. The nebular spectrum of SN 2007bi can be successfully reproduced if the inner region is artificially filled with oxygen-rich, low-velocity ejecta. This most likely requires a grossly aspherical explosion. A major difference between PISN and massive collapse models is that the former emit much more strongly in the NIR. It is concluded that: a) current PISN candidates, in particular SN 2007bi, are more likely the result of the collapse and explosion of massive stars below the PI limit; b) significant asymmetry is required to reproduce the late-time spectrum of SN 2007bi.

We report the discovery of an SN1988Z-like type IIn supernova KISS15s found in a low-mass star-forming galaxy at redshift z=0.038 during the course of the Kiso Supernova Survey (KISS). KISS15s shows long-duration optical continuum and emission line light curves, indicating that KISS15s is powered by a continuous interaction between the expanding ejecta and dense circumstellar medium (CSM). The H$\alpha$ emission line profile can be decomposed into four Gaussians of narrow, intermediate, blue-shifted intermediate, and broad velocity width components, with a full width at half maximum of $\lesssim 100$, $\sim 2,000$, and $\sim 14,000$ km s${}^{-1}$ for the narrow, intermediate, and broad components, respectively. The presence of the blue-shifted intermediate component, of which the line-of-sight velocity relative to the systemic velocity is about $-5,000$ km s${}^{-1}$, suggests that the ejecta-CSM interaction region has an inhomogeneous morphology and anisotropic expansion velocity. We found that KISS15s shows increasing infrared continuum emission, which can be interpreted as hot dust thermal emission of $T \sim 1,200$ K from newly formed dust in a cool, dense shell in the ejecta-CSM interaction region. The progenitor mass-loss rate, inferred from bolometric luminosity, is $\dot{M} \sim 0.4 M_{\odot} \text{yr}^{-1} (v_{w}/40 \text{km}~\text{s}^{-1})$, where $v_{w}$ is the progenitor's stellar wind velocity. This implies that the progenitor of KISS15s was a red supergiant star or a luminous blue variable that had experienced a large mass-loss in the centuries before the explosion.

We present observations of supernova (SN) 2017ens, discovered by the ATLAS survey and identified as a hot blue object through the GREAT program. The redshift z=0.1086 implies a peak brightness of M_g=-21.1 mag, placing the object within the regime of superluminous supernovae. We observe a dramatic spectral evolution, from initially being blue and featureless, to later developing features similar to those of the broadlined Type Ic SN 1998bw, and finally showing ~2000 km s^-1 wide H-alpha and H-beta emission. Relatively narrow Balmer emission (reminiscent of a SN IIn) is present at all times. We also detect coronal lines, indicative of a dense circumstellar medium. We constrain the progenitor wind velocity to ~50-60 km s^-1 based on P-Cygni profiles, which is far slower than those present in Wolf-Rayet stars. This may suggest that the progenitor passed through a luminous blue variable phase, or that the wind is instead from a binary companion red supergiant star. At late times we see the ~2000 km s^-1 wide H-alpha emission persisting at high luminosity (~3x10^40 erg s^-1) for at least 100 day, perhaps indicative of additional mass loss at high velocities that could have been ejected by a pulsational pair instability.

The progenitors of many core-collapse supernovae (CCSNe) are expected to be in binary systems. By performing a series of three-dimensional hydrodynamical simulations, we investigate how CCSN explosions affect their binary companion. We find that the amount of removed stellar mass, the resulting impact velocity, and the chemical contamination of the companion that results from the impact of the SN ejecta, strongly increases with decreasing binary separation and increasing explosion energy. Also, it is foud that the impact effects of CCSN ejecta on the structure of main-sequence (MS) companions, and thus their long term post-explosion evolution, is in general not be dramatic.

Compact neutron star binary systems are produced from binary massive stars through stellar evolution involving up to two supernova explosions. The final stages in the formation of these systems have not been directly observed. We report the discovery of iPTF 14gqr (SN 2014ft), a Type Ic supernova with a fast evolving light curve indicating an extremely low ejecta mass ($\approx 0.2$ solar masses) and low kinetic energy ($\approx 2 \times 10^{50}$ ergs). Early photometry and spectroscopy reveal evidence of shock cooling of an extended He-rich envelope, likely ejected in an intense pre-explosion mass loss episode of the progenitor. Taken together, we interpret iPTF 14gqr as evidence for ultra-stripped supernovae that form neutron stars in compact binary systems.

The energy liberated by fallback accretion has been suggested as a possible engine to power hydrogen-poor superluminous supernovae. We systematically investigate this model using the Bayesian light-curve fitting code MOSFiT (Modular Open Source Fitter for Transients), fitting the light curves of 37 hydrogen-poor superluminous supernovae assuming a fallback accretion central engine. We find that this model can yield good fits to their light curves, with a fit quality that rivals the popular magnetar engine models. Examining our derived parameters for the fallback model, we find the total energy requirements from the accretion disk are estimated to be 0.002 - 0.7 Msun c^2. If we adopt a typical conversion efficiency ~ 1e-3, the required mass to accrete is thus 2 - 700 Msun. Many superluminous supernovae, therefore, require an unrealistic accretion mass, and so only a fraction of these events could be powered by fallback accretion unless the true efficiency is much greater than our fiducial value. The superluminous supernovae that require the smallest amounts of fallback mass still remain to be the fallback accretion powered supernova candidates, but they are difficult to be distinguished solely by their light curve properties.

Context: It has been suggested that some supernovae (SNe) may be powered by a magnetar formed at the moment of the explosion. While this scenario has mostly been applied to hydrogen-free events, it may be possible also for hydrogen-rich objects. Aims: We explore the effect of including a magnetar on the light curves of supernovae with H-rich progenitors. Methods: We have applied a version of our one-dimensional LTE radiation hydrodynamics code that takes into account the relativistic motion of the ejecta caused by the extra energy provided by the magnetar. For a fixed red supergiant (RSG) progenitor, we have obtained a set of light curves that corresponds to different values of the magnetar initial rotation energy and the spin-down timescale. The model is applied to SN~2004em and OGLE-2014-SN-073, two peculiar Type II SNe with long-rising SN1987A-like light curves, although with much larger luminosities. Results: The presence of a plateau phase in either normal or superluminous supernovae is one possible outcome, even if a magnetar is continuously injecting energy into the ejecta.In other cases, the light curve shows a peak but not a plateau. Also, there are intermediate events with a first peak followed by a slow decline and a late break of the declining slope. Our models show that bright and long rising morphologies are possible even assuming RSG structures. Conclusions: A large number of supernova discoveries per year reveal unexpected new types of explosions. According to our results, SLSNe II-P are to be expected, as well as a variety of light curve morphologies that can all be possible signs of a newly born magnetar.

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.

We investigate a progenitor mass distribution of core-collapse supernova remnants (CCSNRs) in our Galaxy and the Large and Small Magellanic Clouds, for the first time. We count the number of CCSNRs in three mass ranges divided by the zero-age main-sequence mass, $M_{\rm ZAMS}$; A: $M_{\rm ZAMS} < 15\ {\rm M}_\odot$, B: $15\ {\rm M}_\odot < M_{\rm ZAMS} < 22.5\ {\rm M}_\odot$, C: $M_{\rm ZAMS} > 22.5\ {\rm M}_\odot$. Simple compilation of progenitor masses in the literature yields a progenitor mass distribution of $f_{\rm A}: f_{\rm B}: f_{\rm C} =0.24:0.28:0.48$, where $f$ is the number fraction of the progenitors. The distribution is inconsistent with any standard initial mass functions. We notice, however, that previous mass estimates are subject to large systematic uncertainties because most of the relative abundances (X/Si) are not really good probe for the progenitor masses. Instead, we propose to rely only on the Fe/Si ratio which is sensitive to the CO core mass ($M_{\rm COcore}$) and $M_{\rm ZAMS}$. Comparing Fe/Si ratios in SNRs in the literature with the newest theoretical model, we estimate 33 $M_{\rm COcore}$ and $M_{\rm ZAMS}$, leading to a revised progenitor mass distribution of $f_{\rm A}: f_{\rm B}: f_{\rm C} = 0.47: 0.32 : 0.21$. This is consistent with the standard Salpeter initial mass function. However, the relation between $M_{\rm COcore}$ and $M_{\rm ZAMS}$ could be affected by binary evolution, which is not taken into account in this study and should be considered in the future work to derive a better progenitor mass distribution estimate.

We model the late evolution and mass loss history of rapidly rotating Wolf-Rayet stars in the mass range $5\,\rm{M}_{\odot}\dots 100\,\rm{M}_{\odot}$. We find that quasi-chemically homogeneously evolving single stars computed with enhanced mixing retain very little or no helium and are compatible with Type\,Ic supernovae. The more efficient removal of core angular momentum and the expected smaller compact object mass in our lower mass models lead to core spins in the range suggested for magnetar driven superluminous supernovae. Our more massive models retain larger specific core angular momenta, expected for long-duration gamma-ray bursts in the collapsar scenario. Due to the absence of a significant He envelope, the rapidly increasing neutrino emission after core helium exhaustion leads to an accelerated contraction of the whole star, inducing a strong spin-up, and centrifugally driven mass loss at rates of up to $10^{-2}\,\rm{M}_{\odot}~\rm{yr^{-1}}$ in the last years to decades before core collapse. Since the angular momentum transport in our lower mass models enhances the envelope spin-up, they show the largest relative amounts of centrifugally enforced mass loss, i.e., up to 25\% of the expected ejecta mass. Our most massive models evolve into the pulsational pair-instability regime. We would thus expect signatures of interaction with a C/O-rich circumstellar medium for Type~Ic superluminous supernovae with ejecta masses below $\sim 10\,\rm{M}_{\odot}$ and for the most massive engine-driven explosions with ejecta masses above $\sim 30\,\rm{M}_{\odot}$. Signs of such interaction should be observable at early epochs of the supernova explosion, and may be related to bumps observed in the light curves of superluminous supernovae, or to the massive circumstellar CO-shell proposed for Type~Ic superluminous supernova Gaia16apd.

We present optical and near-infrared observations of a low-luminosity Type IIP supernova (SN) 2016bkv from the initial rising phase to the plateau phase. Our observations show that the end of the plateau is extended to $\gtrsim 140$ days since the explosion, indicating that this SN takes one of the longest time to finish the plateau phase. among Type IIP SNe (SNe IIP), including low-luminosity (LL) SNe IIP. The line velocities of various ions at the middle of the plateau phase are as low as 1,000--1,500 km s$^{-1}$, which is the lowest even among LL SNe IIP. These measurements imply that the ejecta mass in SN 2016bkv is larger than that of the well-studied LL IIP SN 2003Z. In the early phase, SN 2016bkv shows a strong bump in the light curve. In addition, the optical spectra in this bump phase exhibit a blue continuum accompanied with a narrow H$\alpha$ emission line. These features indicate an interaction between the SN ejecta and the circumstellar matter (CSM) as in SNe IIn. Assuming the ejecta-CSM interaction scenario, the mass loss rate is estimated to be $\sim 1.7 \times 10^{-2} M_{\odot}$ yr$^{-1}$ within a few years before the SN explosion. This is comparable to or even larger than the largest mass loss rate observed for the Galactic red supergiants ($\sim 10^{-3} M_{\odot}$ yr$^{-1}$ for VY CMa). We suggest that the progenitor star of SN 2016bkv experienced a violent mass loss just before the SN explosion.

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

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.

Superluminous supernovae are a new class of supernovae that were recognized about a decade ago. Both observational and theoretical progress has been significant in the last decade. In this review, we first briefly summarize the observational properties of superluminous supernovae. We then introduce the three major suggested luminosity sources to explain the huge luminosities of superluminous supernovae, i.e., the nuclear decay of 56Ni, the interaction between supernova ejecta and dense circumstellar media, and the spin down of magnetars. We compare these models and discuss their strengths and weaknesses.

It is difficult to establish the properties of massive stars that explode as supernovae. The electromagnetic emission during the first minutes to hours after the emergence of the shock from the stellar surface conveys important information about the final evolution and structure of the exploding star. However, the unpredictable nature of supernova events hinders the detection of this brief initial phase. Here we report the serendipitous discovery of a newly born, normal type IIb supernova (SN 2016gkg), which reveals a rapid brightening at optical wavelengths of about 40 magnitudes per day. The very frequent sampling of the observations allowed us to study in detail the outermost structure of the progenitor of the supernova and the physics of the emergence of the shock. We develop hydrodynamical models of the explosion that naturally account for the complete evolution of the supernova over distinct phases regulated by different physical processes. This result suggests that it is appropriate to decouple the treatment of the shock propagation from the unknown mechanism that triggers the explosion.

This contribution presents the results of workshop 8 (supernovae) held during the IAU Symposium 339 "Southern Horizons in Time-Domain Astronomy" at Stellenbosch University, South Africa. Workshop 8 covered a cornucopia of topics with each one having a short presentation by a pre-determined participant followed by a round table discussion. During the first of two sessions, G. Hosseinzadeh and H. Kuncarayakti presented results of their recent works on interacting supernovae. This includes both the intriguing Type Ibn supernova subclass, as well as SN 2017dio, which appears to be the first Type Ic supernova to exhibit signatures of hydrogen-rich circumstellar interaction at all phases. During the second session, M. Sullivan provided an excellent summary related to the future of transient science in the era of Big Data, and the participants discussed strategies to determine which targets and fields should be selected for spectroscopic followup. Workshop 8 concluded with a rather heated discussion regarding the need for the IAU supernova working group to consider to modify the current criteria for a confirmed supernova to receive an official IAU designation.

We propose that massive Thorne-Zytkow objects can explode. A Thorne-Zytkow object is a theoretically predicted star that has a neutron core. When nuclear reactions supporting a massive Thorne-Zytkow object terminate, a strong accretion occurs towards the central neutron core. The accretion rate is large enough to sustain a super-Eddington accretion towards the neutron core. The neutron core may collapse to a black hole after a while. A strong large-scale outflow or a jet can be launched from the super-Eddington accretion disk and the collapsing Thorne-Zytkow object can be turned into an explosion. The ejecta have about 10 Msun but the explosion energy depends on when the accretion is suppressed. We presume that the explosion energy could be as low as ~ 1e47 erg and such a low-energy explosion could be observed like a failed supernova. The maximum possible explosion energy is ~ 1e52 erg and such a high-energy explosion could be observed as an energetic Type II supernova or a superluminous supernova. Explosions of Thorne-Zytkow objects may provide a new path to spread lithium and other heavy elements produced through the irp process such as molybdenum in the Universe.

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.

The progenitor systems of the class of "Ca-rich transients" is a key open issue in time domain astrophysics. These intriguing objects exhibit unusually strong calcium line emissions months after explosion, fall within an intermediate luminosity range, are often found at large projected distances from their host galaxies, and may play a vital role in enriching galaxies and the intergalactic medium. Here we present multi-wavelength observations of iPTF15eqv in NGC 3430, which exhibits a unique combination of properties that bridge those observed in Ca-rich transients and Type Ib/c supernovae. iPTF15eqv has among the highest [Ca II]/[O I] emission line ratios observed to date, yet is more luminous and decays more slowly than other Ca-rich transients. Optical and near-infrared photometry and spectroscopy reveal signatures consistent with the supernova explosion of a < 10 solar mass star that was stripped of its H-rich envelope via binary interaction. Distinct chemical abundances and ejecta kinematics suggest that the core collapse occurred through electron capture processes. Deep limits on possible radio emission made with the Jansky Very Large Array imply a clean environment ($n <$ 0.1 cm$^{-3}$) within a radius of $\sim 10^{17}$ cm. Chandra X-ray Observatory observations rule out alternative scenarios involving tidal disruption of a white dwarf by a black hole, for masses > 100 solar masses). Our results challenge the notion that spectroscopically classified Ca-rich transients only originate from white dwarf progenitor systems, complicate the view that they are all associated with large ejection velocities, and indicate that their chemical abundances may vary widely between events.

We present our study of OGLE-2014-SN-073, one of the brightest Type II SN ever discovered, with an unusually broad lightcurve combined with high ejecta velocities. From our hydrodynamical modelling we infer a remarkable ejecta mass of $60^{+42}_{-16}$~M$_\odot$, and a relatively high explosion energy of $12.4^{+13.0}_{-5.9} \times10^{51}$~erg. We show that this object belongs, with a very small number of other hydrogen-rich SNe, to an energy regime that is not explained by standard core-collapse (CC) neutrino-driven explosions. We compare the quantities inferred by the hydrodynamical modelling with the expectations of various exploding scenarios, trying to explain the high energy and luminosity released. We find some qualitative similarities with pair-instabilities SNe, although a prompt injection of energy by a magnetar seems also a viable alternative to explain such extreme event.

A Type IIn supernova (SN) is dominated by the interaction of SN ejecta with the circumstellar medium (CSM). Some SNe IIn (e.g., SN 2006jd) have episodes of re-brightening ("bumps") in their light curves. We present iPTF13z, a SN IIn discovered by the intermediate Palomar Transient Factory (iPTF) and characterised by several bumps in its light curve. We analyse this peculiar behaviour trying to infer the properties of the CSM and of the SN explosion, as well as the nature of its progenitor star. We obtained multi-band optical photometry for over 1000 days after discovery with the P48 and P60 telescopes at Palomar Observatory. We obtained low-resolution optical spectra in the same period. We did an archival search for progenitor outbursts. We analyse our photometry and spectra, and compare iPTF13z to other SNe IIn. A simple analytical model is used to estimate properties of the CSM. iPTF13z was a SN IIn showing a light curve with five bumps during its decline phase. The bumps had amplitudes between 0.4 and 0.9 mag and durations between 20 and 120 days. The most prominent bumps appeared in all our different optical bands. The spectra showed typical SN IIn characteristics, with emission lines of H$\alpha$ (with broad component FWHM ~$10^{3}-10^{4} ~{\rm ~km ~s^{-1}}$ and narrow component FWHM ~$10^2 \rm ~km ~s^{-1}$) and He I, but also with Fe II, Ca II, Na I D and H$\beta$ P-Cygni profiles (with velocities of ~$10^{3}$ ${\rm ~km ~s^{-1}}$). A pre-explosion outburst was identified lasting $\gtrsim 50$ days, with $M_r \approx -15$ mag around 210 days before discovery. Large, variable progenitor mass-loss rates (~> 0.01 $M_{\odot} \rm ~yr^{-1}$) and CSM densities (~> 10$^{-16}$ g cm$^{-3}$) are derived. We suggest that the light curve bumps of iPTF13z arose from SN ejecta interacting with denser regions in the CSM, possibly produced by the eruptions of a luminous blue variable star.

Although Type Ia supernova cosmology has now reached a mature state, it is important to develop as many independent methods as possible to understand the true nature of dark energy. Recent studies have shown that Type II supernovae (SNe II) offer such a path and could be used as alternative distance indicators. However, the majority of these studies were unable to extend the Hubble diagram above redshift $z=0.3$ because of observational limitations. Here, we show that we are now ready to move beyond low redshifts and attempt high-redshift ($z \gtrsim 0.3$) SN~II cosmology as a result of new-generation deep surveys such as the Subaru/Hyper Suprime-Cam (HSC) survey. Applying the "standard candle method" to SN$\sim$2016jhj ($z=0.3398 \pm 0.0002$; discovered by HSC) together with a low-redshift sample, we are able to construct the highest-redshift SN II Hubble diagram to date with an observed dispersion of 0.27 mag (i.e., 12-13% in distance). This work demonstrates the bright future of SN~II cosmology in the coming era of large, wide-field surveys like that of the Large Synoptic Survey Telescope.

The Crab nebula originated from a core-collapse supernova (SN) explosion observed in 1054 A.D. When viewed as a supernova remnant (SNR), it has an anomalously low observed ejecta mass and kinetic energy for an Fe-core collapse SN. Intensive searches were made for a massive shell that solves this discrepancy, but none has been detected. An alternative idea is that the SN1054 is an electron-capture (EC) explosion with a lower explosion energy by an order of magnitude than Fe-core collapse SNe. In the X-rays, imaging searches were performed for the plasma emission from the shell in the Crab outskirts to set a stringent upper limit to the X-ray emitting mass. However, the extreme brightness of the source hampers access to its vicinity. We thus employed spectroscopic technique using the X-ray micro-calorimeter onboard the Hitomi satellite. By exploiting its superb energy resolution, we set an upper limit for emission or absorption features from yet undetected thermal plasma in the 2-12 keV range. We also re-evaluated the existing Chandra and XMM-Newton data. By assembling these results, a new upper limit was obtained for the X-ray plasma mass of <~ 1Mo for a wide range of assumed shell radius, size, and plasma temperature both in and out of the collisional equilibrium. To compare with the observation, we further performed hydrodynamic simulations of the Crab SNR for two SN models (Fe-core versus EC) under two SN environments (uniform ISM versus progenitor wind). We found that the observed mass limit can be compatible with both SN models if the SN environment has a low density of <~ 0.03 cm-3 (Fe core) or <~ 0.1 cm-3 (EC) for the uniform density, or a progenitor wind density somewhat less than that provided by a mass loss rate of 10-5 Mo yr-1 at 20 km s-1 for the wind environment.

Some ancient, dim, metal-poor stars may have formed in the ashes of the first supernovae (SNe). If their chemical abundances can be reconciled with the elemental yields of specific Population III (Pop III) explosions, they could reveal the properties of primordial stars. But multidimensional simulations of such explosions are required to predict their yields because dynamical instabilities can dredge material up from deep in the ejecta that would otherwise be predicted to fall back on to the central remnant and be lost in one-dimensional (1D) models. We have performed two-dimensional (2D) numerical simulations of two low-energy Pop III SNe, a 12.4 Msun explosion and a 60 Msun explosion, and find that they produce elemental yields that are a good fit to those measured in the most iron-poor star discovered to date, SMSS J031300.36-670839.3 (J031300). Fallback on to the compact remnant in these weak explosions accounts for the lack of measurable iron in J031300 and its low iron-group abundances in general. Our 2D explosions produce higher abundances of heavy elements (atomic number Z > 20) than their 1D counterparts due to dredge-up by fluid instabilities. Since almost no Ni is ejected by these weak SNe, their low luminosities will prevent their detection in the near-infrared with the James Webb Space Telescope and future 30-m telescopes on the ground. The only evidence that they ever occurred will be in the fossil abundance record.

Nascent neutron stars with millisecond periods and magnetic fields in excess of $10^{16}$ Gauss can drive highly energetic and asymmetric explosions known as magnetar-powered supernovae. These exotic explosions are one theoretical interpretation for supernovae Ic-BL which are sometimes associated with long gamma-ray bursts. Twisted magnetic field lines extract the rotational energy of the neutron star and release it as a disk wind or a jet with energies greater than 10$^{52}$ erg over $\sim 20$ sec. What fractions of the energy of the central engine go into the wind and the jet remain unclear. We have performed two-dimensional hydrodynamical simulations of magnetar-powered supernovae (SNe) driven by disk winds and jets with the CASTRO code to investigate the effect of the central engine on nucleosynthetic yields, mixing, and light curves. We find that these explosions synthesize less than 0.05 Msun of Ni and that this mass is not very sensitive to central engine type. The morphology of the explosion can provide a powerful diagnostic of the properties of the central engine. In the absence of a circumstellar medium these events are not very luminous, with peak bolometric magnitudes $M_b \sim -16.5 $ due to low Ni production.

We propose that superluminous transients that appear at central regions of active galactic nuclei (AGNs) such as CSS100217:102913+404220 (CSS100217) and PS16dtm, which reach near or super-Eddington luminosities of the central black holes, are powered by the interaction between accretion disk winds and clouds in broad-line regions (BLRs) surrounding them. If the disk luminosity temporary increases by, e.g., limit-cycle oscillations, leading to a powerful radiatively driven wind, strong shock waves propagate in the BLR. Because the dense clouds in the AGN BLRs typically have similar densities to those found in Type IIn supernovae, strong radiative shocks emerge and efficiently convert the ejecta kinetic energy to radiation. As a result, transients similar to Type IIn supernovae can be observed at AGN central regions. Since a typical black-hole disk wind velocity is ~ 0.1c where c is the speed of light, the ejecta kinetic energy is expected to be ~ 1e52 erg when ~ 1 Msun is ejected. This kinetic energy is transformed to radiation energy in a timescale for the wind to sweep up a similar mass to itself in the BLR, which is a few hundred days. Therefore, both luminosities (~ 1e44 erg/s) and timescales (~ 100 days) of the superluminous transients from AGN central regions match to those expected in our interaction model. If CSS100217 and PS16dtm are related to the AGN activities triggered by limit-cycle oscillations, they become bright again in coming years or decades.

Gravitational waves from the binary black hole (BH) merger GW150914 may enlighten our understanding of ultra-luminous X-ray sources (ULXs), as BHs>30Msun can reach luminosities>4x10^39 erg s^-1 without exceeding their Eddington limit. It is then important to study variations of evolutionary channels for merging BHs, which might instead form accreting BHs and become ULXs. It was recently shown that massive binaries with mass ratios close to unity and tight orbits can undergo efficient rotational mixing and evolve chemically homogeneously, resulting in a compact BH binary. We study similar systems by computing ~120000 detailed binary models with the MESA code covering a wide range of initial parameters. For initial mass ratios M2/M1~0.1-0.4, primaries >40Msun can evolve chemically homogeneously, remaining compact and forming a BH without undergoing Roche-lobe overflow. The secondary then expands and transfers mass to the BH, initiating a ULX phase. We predict that ~1 out of 10^4 massive stars evolves this way, and that in the local universe 0.13 ULXs per Msun yr^-1 of star-formation rate are observable, with a strong preference for low-metallicities. At metallicities log Z>-3, BH masses in ULXs are limited to 60Msun due to the occurrence of pair-instability supernovae which leave no remnant, resulting in an X-ray luminosity cut-off. At lower metallicities, very massive stars can avoid exploding as pair-instability supernovae and instead form BHs with masses above 130Msun, producing a gap in the ULX luminosity distribution. After the ULX phase, neutron-star-BH binaries that merge in less than a Hubble time are produced with a formation rate <0.2 Gpc^-3 yr^-1. We expect that upcoming X-ray observatories will test these predictions, which together with additional gravitational wave detections will provide strict constraints on the origin of the most massive BHs that can be produced by stars.

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.

We present and analyse an extensive dataset of the superluminous supernova (SLSN) LSQ14mo (z = 0.256), consisting of a multi-colour lightcurve from -30 d to +70 d in the rest-frame and a series of 6 spectra from PESSTO covering -7 d to +50 d. This is among the densest spectroscopic coverage, and best-constrained rising lightcurve, for a fast-declining hydrogen-poor SLSN. The bolometric lightcurve can be reproduced with a millisecond magnetar model with ~ 4 M_sol ejecta mass, and the temperature and velocity evolution is also suggestive of a magnetar as the power source. Spectral modelling indicates that the SN ejected ~ 6 M_sol of CO-rich material with a kinetic energy of ~ 7 x 10^51 erg, and suggests a partially thermalised additional source of luminosity between -2 d and +22 d. This may be due to interaction with a shell of material originating from pre-explosion mass loss. We further present a detailed analysis of the host galaxy system of LSQ14mo. PESSTO and GROND imaging show three spatially resolved bright regions, and we used the VLT and FORS2 to obtain a deep (five-hour exposure) spectra of the SN position and the three star-forming regions, which are at a similar redshift. The FORS spectrum at +300 days shows no trace of SN emission lines and we place limits on the strength of [O I] from comparisons with other Ic SNe. The deep spectra provides a unique chance to investigate spatial variations in the host star-formation activity and metallicity. The specific star-formation rate is similar in all three components, as is the presence of a young stellar population. However, the position of LSQ14mo exhibits a lower metallicity, with 12 + log(O/H) = 8.2 in both the R23 and N2 scales (corresponding to ~ 0.3 Z_sol). We propose that the three bright regions in the host system are interacting, which thus triggers star-formation and forms young stellar populations.

Many Type Ic superluminous supernovae have light-curve decline rates after their luminosity peak which are close to the nuclear decay rate of 56Co, consistent with the interpretation that they are powered by 56Ni and possibly pair-instability supernovae. However, their rise times are typically shorter than those expected from pair-instability supernovae, and Type Ic superluminous supernovae are often suggested to be powered by magnetar spin-down. If magnetar spin-down is actually a major mechanism to power Type Ic superluminous supernovae, it should be able to produce decline rates similar to the 56Co decay rate rather easily. In this study, we investigate the conditions for magnetars under which their spin-down energy input can behave like the 56Ni nuclear decay energy input. We find that an initial magnetic field strength within a certain range is sufficient to keep the magnetar energy deposition within a factor of a few of the 56Co decay energy for several hundreds of days. Magnetar spin-down needs to be by almost pure dipole radiation with the braking index close to 3 to mimic 56Ni in a wide parameter range. Not only late-phase 56Co-decay-like light curves, but also rise time and peak luminosity of most 56Ni-powered light curves can be reproduced by magnetars. Bolometric light curves for more than 700 days are required to distinguish the two energy sources solely by them. We expect that more slowly-declining superluminous supernovae with short rise times should be found if they are mainly powered by magnetar spin-down.

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.

We investigate observational properties of accretion-induced collapse (AIC) of white dwarfs in radio frequencies. If AIC is triggered by accretion from a companion star, a dense circumstellar medium can be formed around the progenitor system. Then, the ejecta from AIC collide to the dense circumstellar medium, making a strong shock. The strong shock can produce synchrotron emission which can be observed in radio frequencies. Even if AIC occurs as a result of white dwarf mergers, we argue that AIC may cause fast radio bursts if a certain condition is satisfied. If AIC forms neutron stars which are so massive that rotation is required to support themselves (i.e., supramassive neutron stars), the supramassive neutron stars may immediately lose their rotational energy by the r-mode instability and collapse to black holes. If the collapsing supramassive neutron stars are strongly magnetized, they may emit fast radio bursts as previously suggested. The AIC radio transients from the single-degenerate systems may be detected in the future radio transient surveys like Very Large Array Sky Survey or the Square Kilometer Array transient survey. Because AIC is suggested to be a gravitational wave source, gravitational waves from AIC may be accompanied by radio-bright transients which can be used to confirm the AIC origin of observed gravitational waves.

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).

We investigate the expected rates and bolometric light-curve properties of stripped-envelope electron-capture supernovae (ECSNe) using stellar models from the Binary Population and Spectral Synthesis (BPASS) code. We find that 0.8 per cent (Z = 0.020) and 1.2 per cent (Z = 0.004) of core-collapse supernovae are stripped-envelope ECSNe. Their typical ejecta masses are estimated to be about 0.3 Msun (Z = 0.020) and 0.6 Msun (Z = 0.004). Assuming ECSN explosion properties from numerical explosion simulations, an explosion energy of 1.5e50 erg and a 56Ni mass of 2.5e-3 Msun, we find that stripped-envelope ECSNe have a typical rise time of around 7 days (Z = 0.020) or 13 days (Z = 0.004) and peak luminosity of around 1e41 erg/s (-13.8 mag, Z = 0.020) or 7e40 erg/s (-13.4 mag, Z = 0.004). Their typical ejecta velocities are around 7000 km/s (Z = 0.020) or 5000 km/s (Z = 0.004). Thus, stripped-envelope ECSNe are observed as rapidly-evolving faint transients with relatively small velocities. SN 2008ha-like supernovae, which are the faintest kind of SN 2002cx-like (a.k.a. Type Iax) supernovae, may be related to stripped-envelope ECSNe.

Superluminous supernovae (SLSNe) have been suggested to be powered by strongly magnetized, rapidly rotating neutron stars which are often called magnetars. In this process, rotational energy of the magnetar is radiated via magnetic dipole radiation and heats the supernova ejecta. However, if magnetars are highly distorted in their geometric shape, rotational energy is mainly lost as gravitational wave radiation and thus such magnetars cannot power SLSNe. By simply comparing electromagnetic and gravitational wave emission timescales, we constrain upper limits to the ellipticity of magnetars by assuming that they power the observed SLSNe. We find that their ellipticity typically needs to be less than about a few 1e-3. This indicates that the toroidal magnetic field strengths in these magnetars are typically less than a few 1e16 G so that their distortions remain small. Because light-curve modelling of SLSNe shows that their dipole magnetic field strengths are of the order of 1e14 G, the ratio of poloidal to toroidal magnetic field strengths is found to be larger than ~ 0.01 in magnetars powering SLSNe.

We investigate circumstellar and explosion properties of Type Ibn supernovae (SNe) by analyzing their bolometric light curves. Bolometric light curves of Type Ibn SNe generally have a large contrast between peak luminosity and late-phase luminosity, which is much larger than those of 56Ni-powered SNe. Thus, most of them are likely powered by the interaction between SN ejecta and dense circumstellar media. In addition, Type Ibn SNe decline much faster than Type IIn SNe, and this indicates that the interaction in Type Ibn SNe ceases earlier than in Type IIn SNe. Thus, we argue that Type Ibn SN progenitors experience high mass-loss rates in a short period just before explosion, while Type IIn SN progenitors have high mass-loss rates sustained for a long time. Furthermore, we show that rise time and peak luminosity of Type Ibn and Type IIn SNe are similar and thus, they have similar explosion properties and circumstellar density. The similar circumstellar density in the two kinds of SNe may indicate that mass-loss rates of Type Ibn SN progenitors are generally higher than those of Type IIn as the wind velocities inferred from narrow spectral components are generally higher in Type Ibn SNe. We also show that 56Ni mass and explosion energy of Type Ibn SNe may be smaller than those of other stripped-envelope SNe, probably because they tend to suffer large fallback or some of them may not even be terminal stellar explosions.

Recent advances in gravitational-wave astronomy make the direct detection of gravitational waves from the merger of two stellar-mass compact objects a realistic prospect. Evolutionary scenarios towards mergers of double compact objects generally invoke common-envelope evolution which is poorly understood, leading to large uncertainties in merger rates. We explore the alternative scenario of massive overcontact binary (MOB) evolution, which involves two very massive stars in a very tight binary which remain fully mixed due to their tidally induced high spin. We use the public stellar-evolution code MESA to systematically study this channel by means of detailed simulations. We find that, at low metallicity, MOBs produce double-black-hole (BH+BH) systems that will merge within a Hubble time with mass ratios close to one, in two mass ranges, ~25...60msun and >~ 130msun, with pair instability supernovae (PISNe) being produced in-between. Our models are also able to reproduce counterparts of various stages in the MOB scenario in the local Universe, providing direct support for it. We map the initial parameter space that produces BH+BH mergers, determine the expected chirp mass distribution, merger times, Kerr parameters and predict event rates. We typically find that for Z~<Z_sun/10, there is one BH+BH merger for ~1000 core-collapse supernovae. The advanced LIGO (aLIGO) detection rate is more uncertain and depends on the metallicity evolution. Deriving upper and lower limits from a local and a global approximation for the metallicity distribution of massive stars, we estimate aLIGO detection rates (at design limit) of ~19-550 yr^(-1) for BH+BH mergers below the PISN gap and of ~2.1-370 yr^(-1) above the PISN gap. Even with conservative assumptions, we find that aLIGO should soon detect BH+BH mergers from the MOB scenario and that these could be the dominant source for aLIGO detections.

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.

We propose that nebular H-alpha emission as detected in the Type Ic superluminous supernova iPTF13ehe stems from matter which is stripped from a companion star when the supernova ejecta collide with it. The temporal evolution, the line broadening, and the overall blueshift of the emission are consistent with this interpretation. We scale the nebular H-alpha luminosity predicted for Type Ia supernovae in single-degenerate systems to derive the stripped mass required to explain the H-alpha luminosity of iPTF13ehe. We find a stripped mass of 0.1 - 0.9 solar masses, assuming that the supernova luminosity is powered by radioactivity or magnetar spin down. Because a central heating source is required to excite the H-alpha emission, an interaction-powered model is not favored for iPTF13ehe if the H-alpha emission is from stripped matter. We derive a companion mass of more than 20 solar masses and a binary separation of less than about 20 companion radii based on the stripping efficiency during the collision, indicating that the supernova progenitor and the companion formed a massive close binary system. If Type Ic superluminous supernovae generally occur in massive close binary systems, the early brightening observed previously in several Type Ic superluminous supernovae may also be due to the collision with a close companion. Observations of nebular hydrogen emission in future Type Ic superluminous supernovae will enable us to test this interpretation.

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.

The progenitors of many CCSNe are expected to be in binary systems. After the SN explosion, the companion may suffer from mass stripping and be shock heated as a result of the impact of the SN ejecta. If the binary system is disrupted, the companion is ejected as a runaway and hypervelocity star. By performing a series of 3D hydrodynamical simulations of the collision of SN ejecta with the companion star, we investigate how CCSN explosions affect their companions. We use the BEC code to construct the detailed companion structure at the time of SN explosion. The impact of the SN blast wave on the companion is followed by means of 3D SPH simulations using the Stellar GADGET code. For main-sequence (MS) companions, we find that the amount of removed mass, impact velocity, and chemical contamination of the companion that results from the impact of the SN ejecta, strongly increases with decreasing binary separation and increasing explosion energy. Their relationship can be approximately fitted by power laws, which is consistent with the results obtained from impact simulations of SNe~Ia. However, we find that the impact velocity is sensitive to the momentum profile of the outer SN ejecta and, in fact, may decrease with increasing ejecta mass, depending on the modeling of the ejecta. Because most companions to Ib/c CCSNe are in their MS phase at the moment of the explosion, combined with the strongly decaying impact effects with increasing binary separation, we argue that the majority of these SNe lead to inefficient mass stripping and shock heating of the companion star following the impact of the ejecta. Our simulations show that the impact effects of Ib/c SN ejecta on the structure of MS companions, and thus their long-term post-explosion evolution, is in general not dramatic. We find that at most 10% of their mass is lost, and their resulting impact velocities are less than 100 km/s.

The progenitors of SNe Ia are not yet fully understood. The two leading progenitor scenarios are the single-degenerate (SD) scenario and the double-degenerate scenario. In the SD scenario, the collision of the SN Ia ejecta with its companion star is expected to produce detectable ultraviolet (UV) emission in the first few days after the SN explosion within certain viewing angles. A strong UV flash has recently been detected in an SN 2002es-like peculiar SN Ia iPTF14atg by Cao et al., which is interpreted as evidence of an early-time UV signature due to SN ejecta interacting with its companion star, supporting the SD scenario. In this paper, we present the expected luminosity distributions of early-time UV emission arising from SN Ia ejecta-companion interaction by performing binary population synthesis calculations for different progenitor systems in the SD scenario. Our theoretical predictions will be helpful for future early-time observations of SNe Ia to constrain their possible progenitors. Assuming the observed strong UV pulse of iPTF14atg was indeed produced by the SN ejecta-companion interaction, our population synthesis model suggests that the progenitor system of iPTF14atg is most likely a red-giant donor binary system, and it is unlikely to have been a main-sequence or helium-star donor system.

Massive, luminous stars reaching the Eddington limit in their interiors develop very dilute, extended envelopes. This effect is called envelope inflation. If the progenitors of Type Ib/c supernovae, which are believed to be Wolf-Rayet (WR) stars, have inflated envelopes then the shock breakout signals diffuse in them and can extend their rise times significantly. We show that our inflated, hydrogen-free, WR stellar models with a radius of ~Rsun can have shock breakout signals longer than ~60 s. The puzzlingly long shock breakout signal observed in the Type Ib SN 2008D can be explained by an inflated progenitor envelope, and more such events might argue in favour of existence of inflated envelopes in general.

We report the late-time evolution of Type IIb Supernova (SN IIb) 2013df. SN 2013df showed a dramatic change in its spectral features at ~1 year after the explosion. Early on it showed typical characteristics shared by SNe IIb/Ib/Ic dominated by metal emission lines, while later on it was dominated by broad and flat-topped Halpha and He I emissions. The late-time spectra are strikingly similar to SN IIb 1993J, which is the only previous example clearly showing the same transition. This late-time evolution is fully explained by a change in the energy input from the $^{56}$Co decay to the interaction between the SN ejecta and dense circumstellar matter (CSM). The mass loss rate is derived to be ~(5.4 +- 3.2) x 10^{-5} Msun/yr (for the wind velocity of ~20 km/s), similar to SN 1993J but larger than SN IIb 2011dh by an order of magnitude. The striking similarity between SNe 2013df and 1993J in the (candidate) progenitors and the CSM environments, and the contrast in these natures to SN 2011dh, infer that there is a link between the natures of the progenitor and the mass loss: SNe IIb with a more extended progenitor have experienced a much stronger mass loss in the final centuries toward the explosion. It might indicate that SNe IIb from a more extended progenitor are the explosions during a strong binary interaction phase, while those from a less extended progenitor have a delay between the strong binary interaction and the explosion.

Many core-collapse supernova progenitors are presumed to be in binary systems. If a star explodes in a binary system, the early supernova light curve can be brightened by the collision of the supernova ejecta with the companion star. The early brightening can be observed when the observer is in the direction of the hole created by the collision. Based on a population synthesis model, we estimate the fractions of core-collapse supernovae in which the light-curve brightening by the collision can be observed. We find that 0.19% of core-collapse supernova light curves can be observed with the collisional brightening. Type Ibc supernova light curves are more likely to be brightened by the collision (0.53%) because of the high fraction of the progenitors being in binary systems and their proximity to the companion stars. Type II and IIb supernova light curves are less affected (~1e-3% and ~1e-2%, respectively). Although the early, slow light-curve declines of some Type IIb and Ibc supernovae are argued to be caused by the collision with the companion star (e.g. SN 2008D), the small expected fraction, as well as the unrealistically small separation required, disfavour the argument. The future transient survey by the Large Synoptic Survey Telescope is expected to detect ~10 Type Ibc supernovae with the early collisional brightening per year, and they will be able to provide information on supernova progenitors in binary systems.

We construct spectra of supernovae interacting strongly with a circumstellar medium (CSM) by adding SN templates, a black-body continuum and an emission-line spectrum. In a Monte Carlo simulation we generate more than 800 spectra, distribute them to 10 different classifiers and study how the different simulation parameters affect the appearance of the spectra. SNe IIn showing some structure over the continuum were characterized as SNe IInS. We demonstrate that the flux ratio of the underlying SN to the continuum fv is the most important parameter determining the spectral classification. Thermonuclear SNe get progressively classified as Ia-CSM, IInS and IIn as fv decreases. The transition between Ia-CSM and IInS occurs at fv~0.2-0.3. It shown that SNe Ia-CSM are found at the magnitude range -19.5> M >-21.6, in good agreement with observations, and that the faintest SN IIn that can hide a SN Ia has M = -20.1. The sample of SNe Ia-CSM shows an association with 91T-like SNe Ia. Our experiment does not support that this association is due to a luminosity bias (91T-like being brighter than normal Ia). We therefore conclude that this association has real physical origins and we propose that 91T-like explosions result from single degenerate progenitors. Despite the similarities between SNe Ibc and SNe Ia, the number of misclassifications between these types was small and mostly at low S/N. Combined with the luminosity function needed to reproduce the observed SN Ia-CSM luminosities, it is unlikely that SNe Ibc constitute an important contaminant within this sample. We show how Type II spectra transition to IIn and how the Ha profiles vary with fv. SNe IIn fainter than M = -17.2 are unable to mask SNe IIP brighter than M = -15. A simulation including radiative transfer shows that our simplified model is a good first order approximation. The spectra obtained are in good agreement with real data.

We constrain the explosion and circumstellar properties at the 2012b event of SN 2009ip based on its late-phase bolometric light curve recently reported. The explosion energy and ejected mass at the 2012b event are estimated as 0.01 Msun and 2e49 erg, respectively. The circumstellar medium is assumed to have two components: an inner shell and an outer wind. The inner shell which is likely created at the 2012a event has 0.2 Msun. The outer wind is created by the wind mass loss before the 2012a mass ejection, and the progenitor is estimated to have had the mass-loss rate about 0.1 Msun/yr with the wind velocity 550 km/s before the 2012a event. The estimated explosion energy and ejected mass indicate that the 2012b event is not caused by a regular supernova.

Stars close to the Eddington luminosity can have large low-density inflated envelopes. We show that the rise times of shock breakout signals from supernovae can be extended significantly if supernova progenitors have an inflated stellar envelope. If the shock breakout occurs in such inflated envelopes, the shock breakout signals diffuse in them, and their rise time can be significantly extended. Then, the rise times of the shock breakout signals are dominated by the diffusion time in the inflated envelope rather than the light-crossing time of the progenitors. We show that our inflated Wolf-Rayet star models whose radii are of the order of the solar radius can have shock breakout signals which are longer than ~100 sec. The existence of inflated envelopes in Wolf-Rayet supernova progenitors may be related to the mysterious long shock breakout signal observed in Type Ib SN 2008D. Extended shock breakout signals may provide evidence for the existence of inflated stellar envelopes and can be used to constrain the physical properties of these enigmatic structures.

SNe Iax are proposed as one new sub-class of SNe Ia. SNe Iax have been estimated to account for ~5-30% of the total SN Ia rate, and most SNe Iax have been discovered in late-type galaxies. In addition, observations constrain the progenitor systems of some SN Iax progenitors have ages of <80Myr. Recently, the weak deflagration explosions of Ch-mass CO WDs seem to provide a viable physical scenario for SNe Iax. Here, comparing theoretical predictions from BPS calculations with SN Iax observations, we put constraints on the SD Ch-mass model as a possible Iax progenitor. A long delay times of >3Gyr and low SN rates of ~3e-5/yr are found in our RG donor channel, indicating that this channel is unlikely to produce SNe Iax. We predict that the Galactic rate from the MS (He) donor channel is ~1.5e-3/yr (3e-4/yr), which is consistent with the observed Iax rate. The short delay times in the He channel support the young host environments of SNe Iax. However, the relative long delay times in the MS donor channel are less favourable for the observational constraints on the Iax progenitor ages. Finally, we set an upper limit on the pre-SN mass-loss rate at < 10e-4 M_sun/yr (for V_wind=100 km/s). The delay times in the SD Ch-mass model do not account for that most SNe Iax are located in late-type galaxies. However, at least one Iax event (SN 2008ge) is hosted by a S0 galaxy with no signs of star formation. Moreover, current X-ray observations cannot rule out the SD Ch-mass model. Taking all these into account and considering the uncertainty of the observed rate for SNe Iax, we suggest that some SNe Iax may be produced from weak deflagrations of Ch-mass CO WDs in SD progenitors, especially in the He-star channel. This is consistent with recent analysis of HST observations for the SN Iax SN2012Z. However, this SD deflagration model is still unlikely to be the most common progenitor scenario for SNe Iax.

Recent stellar evolution models show consistently that very massive metal-free stars evolve into red supergiants shortly before they explode. We argue that the envelopes of these stars, which will form pair-instability supernovae, become pulsationally unstable and that this will lead to extreme mass-loss rates despite the tiny metal content of the envelopes. We investigate the pulsational properties of such models and derive pulsationally induced mass-loss rates, which take the damping effects of the mass loss on the pulsations selfconsistently into account. We find that the pulsations may induce mass-loss rates of ~ 1e-4 - 1e-2 Msun/yr shortly before the explosions, which may create a dense circumstellar medium. Our results show that very massive stars with dense circumstellar media may stem from a wider initial mass range than pulsational-pair instability supernovae. The extreme mass loss will cease when so much of the hydrogen-rich envelope is lost that the star becomes more compact and stops pulsating. The helium core of these stars therefore remains unaffected, and their fate as pair-instability supernovae remains unaltered. The existence of dense circumstellar media around metal-free pair-instability supernovae can make them brighter and bluer, and they may be easier to detect at high redshifts than previously expected. We argue that the mass-loss enhancement in pair-instability supernova progenitors can naturally explain some observational properties of superluminous supernovae: the energetic explosions of stars within hydrogen-rich dense circumstellar media with little 56Ni production and the lack of a hydrogen-rich envelope in pair-instability supernova candidates with large 56Ni production.

Betelgeuse, a nearby red supergiant, is a runaway star with a powerful stellar wind that drives a bow shock into its surroundings. This picture has been challenged by the discovery of a dense and almost static shell that is three times closer to the star than the bow shock and has been decelerated by some external force. The two physically distinct structures cannot both be formed by the hydrodynamic interaction of the wind with the interstellar medium. Here we report that a model in which Betelgeuse's wind is photoionized by radiation from external sources can explain the static shell without requiring a new understanding of the bow shock. Pressure from the photoionized wind generates a standing shock in the neutral part of the wind and forms an almost static, photoionization-confined shell. Other red supergiants should have significantly more massive shells than Betelgeuse, because the photoionization-confined shell traps up to 35 per cent of all mass lost during the red supergiant phase, confining this gas close to the star until it explodes. After the supernova explosion, massive shells dramatically affect the supernova lightcurve, providing a natural explanation for the many supernovae that have signatures of circumstellar interaction.

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.

We investigate the diversity in the wind density, supernova ejecta energy, and ejecta mass in Type IIn supernovae based on their rise times and peak luminosities. We show that the wind density and supernova ejecta properties can be estimated independently if both the rise time and peak luminosity are observed. The peak luminosity is mostly determined by the supernova properties and the rise time can be used to estimate the wind density. We find that the ejecta energy of Type IIn supernovae needs to vary by factors of 0.2-5 from the average if their ejecta mass is similar. The diversity in the observed rise times indicates that their wind density varies by factors of 0.2-2 from the average. We show that Type IIn superluminous supernovae should have not only large wind density but also large ejecta energy and/or small ejecta mass to explain their large luminosities and the rise times at the same time. We also note that the shock breakout does not necessarily occur in the wind even if it is optically thick, except for the case of superluminous supernovae, and we analyze the observational data both with and without assuming that the shock breakout occurs in the dense wind of Type IIn supernovae.

Objects of known brightness, like Type Ia supernovae (SNIa), can be used to measure distances. If a massive object warps spacetime to form multiple images of a background SNIa, a direct test of cosmic expansion is also possible. However, these lensing events must first be distinguished from other rare phenomena. Recently, a supernova was found to shine much brighter than normal for its distance, which resulted in a debate: was it a new type of superluminous supernova or a normal SNIa magnified by a hidden gravitational lens? Here we report that a spectrum obtained after the supernova faded away shows the presence of a foreground galaxy--the first found to strongly magnify a SNIa. We discuss how more lensed SNIa may be found than previously predicted.

We present a novel mechanism to enhance the mass-loss rates of massive stars shortly before their explosion. The neutrino luminosities of the stellar core of massive stars become larger as they get closer to the time of the core collapse. As emitted neutrinos escape freely from the core, the core mass is significantly reduced when the neutrino luminosity is large. If a star is near the Eddington luminosity when the neutrino luminosity is large, the star can exceed the Eddington luminosity because of the core neutrino mass loss. We suggest that the stellar surface mass-loss rates due to the core neutrino emission can be higher than 1e-4 Msun/yr from ~ 1 year before the core collapse. The mass-loss rates can exceed 1e-2 Msun/yr in ~ 10 days before the core collapse. This mass-loss mechanism may be able to explain the enhanced mass loss observed in some supernova progenitors shortly before their explosion. Even if the star is not close enough to the Eddington luminosity to enhance mass loss, the star can still expand because of the reduced gravitational force. This mechanism can be activated in Wolf-Rayet stars and it can make hydrogen-poor, as well as hydrogen-rich, dense circumstellar media observed in some supernovae.

Recently, Moriya et al. (2013) developed an analytic bolometric light curve model for supernovae interacting with dense circumstellar media. Because of the dense circumstellar medium, the shocked region is assumed to be radiative and make a thin dense shell. The model is based on the conservation of momentum in the shocked dense shell. However, the analytic model was mentioned to neglect the 'snow-plow' phase of the shocked dense shell by Ofek et al. (2014). The 'snow-plow' or momentum-conserving phase refers to the period in which the momentum injection from the supernova ejecta is almost terminated and the radiative shocked dense shell keeps moving due only to the momentum previously provided by the supernova ejecta. In this Note, I clarify that the analytic model of Moriya et al. (2013) does take the 'snow-plow' phase into account and the criticism of Ofek et al. (2014) is incorrect. In addition, Ofek et al. (2014) related the sudden luminosity break observed in the light curve of Type IIn SN 2010jl to the transition to the 'snow-plow' phase. However, I argue that the sudden transition to the 'snow-plow' phase is not consistent with the luminosity break observed in SN 2010jl. The luminosity break is likely to be related to other phenomena like the dense shell exiting the dense part of the circumstellar medium.

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.

We discuss four questions dealing with massive star evolution. The first one is about the origin of slowly rotating, non-evolved, nitrogen rich stars. We propose that these stars may originate from initially fast rotating stars whose surface has been braked down. The second question is about the evolutionary status of alpha-Cygni variables. According to their pulsation properties, these stars should be post red supergiant stars. However, some stars at least present surface abundances indicating that they should be pre red supergiant stars. How to reconcile these two contradictory requirements? The third one concerns the various supernova types which are the end point of the evolution of stars with initial masses between 18 and 30 Msun, i.e. the most massive stars which go through a red supergiant phase during their lifetime. Do they produce types IIP, IIL, IIn, IIb or Ib supernovae or do they end without producing any SN event? Finally, we shall discuss reasons why so few progenitors of type Ibc supernovae have yet been detected?

Observational identification of the first stars is one of the great challenges in the modern astronomy. Although a single first star is too faint to be detected, supernova explosions of the first stars can be bright enough. An important question is whether such supernovae can be detected in the limited observational area with realistic observational resources. We perform detailed simulations to study the detectability of superluminous supernovae (SLSNe) at high redshifts, using the observationally-calibrated star formation rate density and supernova occurrence rate. We show that a 100 deg^2 survey with the limiting magnitude of 26 mag in near-infrared wavelengths will be able to discover about 10 SLSNe at z > 10. If the survey is extended to 200 deg^2 with 27 mag depth, about 10 SLSNe can be discovered at z > 15. We emphasize that the observations at > 3 um are crucial to detect and select SLSNe at z > 10. Our simulations are also applied to the planned survey with Euclid, WFIRST, and WISH. These surveys will be able to detect about 1000, 400, and 3000 SLSNe up to z ~ 5, 7, and 12, respectively. We conclude that detection of SLSNe at z > 10 is in fact achievable in the near future.

Supernovae (SNe) have been proposed to be the main production sites of dust grains in the Universe. Our knowledge on their importance to dust production is, however, limited by observationally poor constraints on the nature and amount of dust particles produced by individual SNe. In this paper, we present a spectrum covering optical through near-Infrared (NIR) light of the luminous Type IIn supernova (SN IIn) 2010jl around one and half years after the explosion. This unique data set reveals multiple signatures of newly formed dust particles. The NIR portion of the spectrum provides a rare example where thermal emission from newly formed hot dust grains is clearly detected. We determine the main population of the dust species to be carbon grains at a temperature of ~1,350 - 1,450K at this epoch. The mass of the dust grains is derived to be ~(7.5 - 8.5) x 10^{-4} Msun. Hydrogen emission lines show wavelength-dependent absorption, which provides a good estimate on the typical size of the newly formed dust grains (~0.1 micron, and most likely <~0.01 micron). We attribute the dust grains to have been formed in a dense cooling shell as a result of a strong SN-circumstellar media (CSM) interaction. The dust grains occupy ~10% of the emitting volume, suggesting an inhomogeneous, clumpy structure. The average CSM density is required to be >~3 x 10^{7} cm^{-3}, corresponding to a mass loss rate of >~0.02 Msun yr^{-1} (for a mass loss wind velocity of ~100 km s^{-1}). This strongly supports a scenario that SN 2010jl and probably other luminous SNe IIn are powered by strong interactions within very dense CSM, perhaps created by Luminous Blue Variable (LBV)-like eruptions within the last century before the explosion.

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.

<Context> Ideally, one would like to know which type of core-collapse SNe is produced by different progenitors and the channels of stellar evolution leading to these progenitors. These links have to be very well known to use the observed frequency of different types of SN events for probing the star formation rate and massive star evolution in different types of galaxies. <Aims> We investigate the link between LBV as SN progenitors and the appearance of episodic radio light curve modulations of the SN event. <Methods> We use the 20Msun and 25Msun models with rotation at solar metallicity, part of an extended grid of stellar models computed by the Geneva team. At their pre-SN stage, these two models have recently been shown to have spectra similar to those of LBV stars and possibly explode as Type IIb SNe. Based on the wind properties before the explosion, we derive the density structure of their circumstellar medium. This structure is used as input for computing the SN radio light curve. <Results> We find that the 20Msun model shows radio light curves with episodic luminosity modulations, similar to those observed in some Type IIb SNe. This occurs because the evolution of the 20Msun model terminates in a region of the HR diagram where radiative stellar winds present strong density variations, caused by the bistability limit. The 25Msun model, ending its evolution in a zone of the HR diagram where no change of the mass-loss rates is expected, presents no such modulations in its radio SN light curve. <Conclusions> Our results reinforce the link between SN progenitors and LBV stars. We also confirm the existence of a physical mechanism for a single star to have episodic radio light curve modulations. In the case of the 25Msun progenitors, we do not obtain modulations in the radio light curve, but our models may miss some outbursting behavior in the late stages of massive stars.

Recently, Chornock and co-workers announced the Pan-STARRS discovery of a transient source reaching an apparent peak luminosity of ~4x10^44 erg s^-1. We show that the spectra of this transient source are well fit by normal Type Ia supernova (SNIa) templates. The multi-band colors and light-curve shapes are also consistent with normal SNeIa at the spectroscopically determined redshift of z=1.3883; however, the observed flux is a constant factor of ~30 times too bright in each band over time as compared to the templates. At minimum, this shows that the peak luminosities inferred from the light-curve widths of some SNeIa will deviate significantly from the established, empirical relation used by cosmologists. We argue on physical grounds that the observed fluxes do not reflect an intrinsically luminous SNIa, but rather PS1-10afx is a normal SNIa whose flux has been magnified by an external source. The only known astrophysical source capable of such magnification is a gravitational lens. Given the lack of obvious lens candidates, such as galaxy clusters, in the vicinity, we further argue that the lens is a supermassive black hole or a comparatively low-mass dark matter halo. In this case, the lens continues to magnify the underlying host galaxy light. If confirmed, this discovery could impact a broad range of topics including cosmology, gamma-ray bursts, and dark matter halos.

We present a study of SN 2009js in NGC 918. Multi-band Kanata optical photometry covering the first ~120 days show the source to be a Type IIP SN. Reddening is dominated by that due to our Galaxy. One-year-post-explosion photometry with the NTT, and a Subaru optical spectrum 16 days post-discovery, both imply a good match with the well-studied subluminous SN 2005cs. The plateau phase luminosity of SN 2009js and its plateau duration are more similar to the intermediate luminosity IIP SN 2008in. Thus, SN 2009js shares characteristics with both subluminous and intermediate luminosity SNe. Its radioactive tail luminosity lies between SN 2005cs and SN 2008in, whereas its quasi-bolometric luminosity decline from peak to plateau (quantified by a newly-defined parameter Delta[logL] measuring adiabatic cooling following shock breakout) is much smaller than both the others. We estimate the ejected mass of 56Ni to be low (~0.007 Msun). The SN explosion energy appears to have been small, similar to SN 2005cs. SN 2009js is the first subluminous SN IIP to be studied in the mid-infrared. It was serendipitously caught by Spitzer at very early times. In addition, it was detected by WISE 105 days later with a significant 4.6 micron flux excess above the photosphere. The infrared excess luminosity relative to the photosphere is clearly smaller than that of SN 2004dj extensively studied in the mid-infrared. The excess may be tentatively assigned to heated dust with mass ~3e-5 Msun, or to CO fundamental emission as a precursor to dust formation.

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.

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.

The origin of super-luminous supernovae (SLSNe), especially the source of their huge luminosities, has not been clarified yet. While a strong interaction between SN ejecta and dense circumstellar media (CSM) is a leading scenario, alternative models have been proposed. In this Letter, we suggest new diagnostics to discriminate the strong SN-CSM interaction scenario from the others: a decline in the luminosity ('dip') before the main peak of the light curve. This dip is an unavoidable consequence of having a dense CSM within which the shock breakout occurs. If a dense CSM shell is located far at large radii from the progenitor inside, it takes time for the SN ejecta to reach it and the early light curve can be powered by the SN ejecta before the collision. Once the SN ejecta collides with the dense CSM, the electron density and thus the Thomson scattering opacity suddenly increase. Photons become unable to go out of the shock even if there is a source of emission inside, which results in the dip in the light curve. This dip is a solid prediction from the strong interaction scenario irrespective of a power source for the early emission. Eventually the forward shock breaks out from within the dense CSM, and the luminosity increases by the continuous strong SN-CSM interaction, resulting in an SLSN. The possible dip observed in the hydrogen-poor SLSN, 2006oz, could be the first example of this signature and give support to the SN-CSM interaction as the power source of SLSN 2006oz.

Usual supernova remnants have either ionizing plasma or plasma in collisional ionization equilibrium, i.e., the ionization temperature is lower than or equal to the electron temperature. However, the existence of recombining supernova remnants, i.e., supernova remnants with the ionization temperature higher than the electron temperature, is recently confirmed. One suggested way to have recombining plasma in a supernova remnant is to have a dense circumstellar medium at the time of the supernova explosion. If the circumstellar medium is dense enough, collisional ionization equilibrium can be established in the early stage of the evolution of the supernova remnant and subsequent adiabatic cooling which occurs after the shock wave gets out of the dense circumstellar medium makes the electron temperature lower than the ionization temperature. We study the circumstellar medium around several supernova progenitors and show which supernova progenitors can have a circumstellar medium which is dense enough to establish collisional ionization equilibrium soon after the explosion. We find that the circumstellar medium around red supergiants (especially massive ones) and the circumstellar medium which is dense enough to make Type IIn supernovae can establish collisional ionization equilibrium soon after the explosion and can evolve to recombining supernova remnants. Wolf-Rayet stars and white dwarfs have possibility to be recombining supernova remnants but the fraction is expected to be very small. As the occurrence rate of the explosions of red supergiants is much higher than that of Type IIn supernovae, the major progenitors of recombining supernova remnants are likely to be red supergiants.

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.

Observations of high-redshift supernovae (SNe) open a novel opportunity to study the massive star population in the early Universe. We study the detectability of superluminous SNe with upcoming optical and near-infrared (NIR) surveys. Our calculations are based on the cosmic star formation history, the SN occurence rate, the characteristic colour and the light curve of the SNe that are all calibrated by available observations. We show that 15-150 SNe up to z ~ 4 will be discovered by the proposed Subaru/Hyper Suprime-Cam deep survey: 30 deg^2 survey with 24.5 AB mag depth in z-band for 3 months. With its ultra-deep layer (3.5 deg^2 with 25.6 AB mag depth in z-band for 4 months), the highest redshift can be extended to z ~ 5. We further explore the detectability by upcoming NIR survey utilizing future satellites such as Euclid, WFIRST, and WISH. The wide-field NIR surveys are very efficient to detect high-redshift SNe. With a hypothetical deep NIR survey for 100 deg^2 with 26 AB mag depth at 1-4 um, at least ~ 50 SNe will be discovered at z>3 in half a year. The number of the detected SNe can place a strong constraint on the stellar initial mass function or its slope especially at the high-mass end. Superluminous SNe at high redshifts can be distinguished from other types of SNe by the long time-scale of their light curves in the observer's frame, the optical colours redder than other core-collapse SNe and the NIR colours redder than any other types of SNe.

We show that the diversity in the density slope of the dense wind due to non-steady mass loss can be one way to explain the spectral diversity of Type II luminous supernovae (LSNe). The interaction of SN ejecta and wind surrounding it is considered to be a power source to illuminate LSNe because many LSNe show the wind signature in their spectra (Type IIn LSNe). However, there also exist LSNe without the spectral features caused by the wind (Type IIL LSNe). We show that, even if LSNe are illuminated by the interaction, it is possible that they do not show the narrow spectra from the wind if we take into account of non-steady mass loss of their progenitors. When the shock breakout takes place in the dense wind with the density structure \rho\propto r^{-w}, the ratio of the diffusion timescale in the optically thick region of the wind (td) and the shock propagation timescale of the entire wind after the shock breakout (ts) strongly depends on w. For the case w<~1, both timescales are comparable (td/ts \simeq 1) and td/ts gets smaller as w gets larger. For the case td/ts\simeq 1, the shock goes through the entire wind just after the light curve (LC) peak and narrow spectral lines from the wind cannot be observed after the LC peak (Type IIL LSNe). If td/ts is much smaller, the shock wave continues to propagate in the wind after the LC peak and unshocked wind remains (Type IIn LSNe). This difference can be obtained only through a careful treatment of the shock breakout condition in a dense wind. The lack of narrow Lorentzian line profiles in Type IIL LSNe before the LC peak can also be explained by the difference in the density slope. Furthermore, we apply our model to Type IIn LSN 2006gy and Type IIL LSN 2008es and find that our model is consistent with the observations.

We perform hydrodynamical calculations of core-collapse supernovae (SNe) with low explosion energies. These SNe do not have enough energy to eject the whole progenitor and most of the progenitor falls back to the central remnant. We show that such fallback SNe can have a variety of light curves (LCs) but their photospheric velocities can only have some limited values with lower limits. We also perform calculations of nucleosynthesis and LCs of several fallback SN model, and find that a fallback SN from the progenitor with a main-sequence mass of 13 Msun can account for the properties of the peculiar Type Ia supernova SN 2008ha. The kinetic energy and ejecta mass of the model are 1.2*10^{48} erg and 0.074 Msun, respectively, and the ejected 56Ni mass is 0.003 Msun. Thus, SN 2008ha can be a core-collapse SN with a large amount of fallback. We also suggest that SN 2008ha could have been accompanied with long gamma-ray bursts and long gamma-ray bursts without associated SNe may be accompanied with very faint SNe with significant amount of fallback which are similar to SN 2008ha.

We present the results of near-infrared (NIR) multi-epoch observations of the optical transient in the nearby galaxy NGC300 (NGC300-OT) at 398 and 582 days after the discovery with the Infrared Camera (IRC) onboard AKARI. NIR spectra (2--5 um) of NGC300-OT were obtained for the first time. They show no prominent emission nor absorption features, but are dominated by continuum thermal emission from the dust around NGC300-OT. NIR images were taken in the 2.4, 3.2, and 4.1 um bands. The spectral energy distributions (SED) of NGC300-OT indicate the dust temperature of 810 (+-14) K at 398 days and 670 (+-12) K at 582 days. We attribute the observed NIR emission to the thermal emission from dust grains formed in the ejecta of NGC300-OT. The multi-epoch observations enable us to estimate the dust optical depth as larger than about 12 at 398 days and larger than about 6 at 582 days at 2.4 um, by assuming an isothermal dust cloud. The observed NIR emission must be optically thick, unless the amount of dust grains increases with time. Little extinction at visible wavelengths reported in earlier observations suggests that the dust cloud around NGC300-OT should be distributed inhomogeneously so as to not screen the radiation from the ejecta gas and the central star. The present results suggest the dust grains are not formed in spherically symmetric geometry, but rather in a torus, a bipolar outflow, or clumpy cloudlets.

We present a core-collapse supernova model for the extremely luminous Type Ic supernova 2007bi. By performing numerical calculations of hydrodynamics, nucleosynthesis, and radiation transport, we find that SN 2007bi is consistent with the core-collapse supernova explosion of a 43 Msun carbon and oxygen core obtained from the evolution of a progenitor star with a main sequence mass of 100 Msun and metallicity of Z = Zsun/200, from which its hydrogen and helium envelopes are artificially stripped. The ejecta mass and the ejecta kinetic energy of the models are 40 Msun and 3.6*10^{52} erg. The ejected 56Ni mass is as large as 6.1 Msun, which results from the explosive nucleosynthesis with large explosion energy. We also confirm that SN 2007bi is consistent with a pair-instability supernova model as has recently been claimed. We show that the earlier light curve data can discriminate between the models for such luminous supernovae.

A new mechanism for enhancing the helium abundance in the blue main sequence stars of omega Centauri and NGC 2808 is investigated. We suggest that helium enhancement was caused by the inhomogeneous big bang nucleosynthesis. Regions with extremely high baryon-to-photon ratios are assumed to be caused by the baryogenesis. Its mass scale is also assumed to be 10^6 Msun. An example of the mechanisms to realize these two things was already proposed as the Affleck-Dine baryogenesis. As the baryon-to-photon ratio becomes larger, the primordial helium abundance is enhanced. We calculated the big bang nucleosynthesis and found that there exists a parameter region yielding enough helium to account for the split of the main sequence in the aforementioned globular clusters while keeping the abundances of other elements compatible with observations. Our mechanism predicts that heavy elements with the mass number of around 100 is enhanced in the blue main sequence stars. We estimate the timescales of diffusion of the enhanced helium and mass accretion in several stages after the nucleosynthesis to investigate whether these processes diminish the enhancement of helium. We found that the diffusion does not influence the helium content. A cloud with a sufficiently large baryon-to-photon ratio to account for the multiple main sequence collapsed immediately after the recombination. Subsequently, the cloud accreted the ambient matter with the normal helium content. If the star formation occurred both in the collapsed core and the accreted envelope, then the resultant star cluster has a double main sequence.