We present deep (9 hours) Gemini-N/GNIRS near-infrared spectroscopic observations of ULAS J1342+0928, a luminous quasar at z=7.54. Various broad emission lines were detected, as well as the underlying continuum and iron forests over the rest-frame wavelength 970-2930A. There is a clear trend that higher-ionization emission lines show larger blueshifts with CIV1549 exhibiting 5510^{+240}_{-110} km s-1 blueshift with respect to the systematic redshift from the far-infrared [CII] 158um emission line. Those high ionization lines have wide profiles with FWHM more than 10000 km s-1. A modest blueshift of 340^{+110}_{-80} km s-1 is also seen in MgII, the lowest ionization line identified in the spectrum. The updated MgII-based black hole mass of M_BH=9.1_{-1.3}^{+1.4} x 10^8 M_sun and the Eddington ratio of L_bol/L_Edd=1.1_{-0.2}^{+0.2} confirm that ULAS J1342+0928 is powered by a massive and actively accreting black hole. There is no significant difference in the emission line ratios such as SiIV/CIV and AlIII/CIV when compared to lower-redshift quasars in a similar luminosity range, which suggests early metal pollution of the broad-line-region clouds. This trend also holds for the FeII/MgII line ratio, known as a cosmic clock that traces the iron enrichment in the early universe. Different iron templates and continuum fitting ranges were used to explore how the FeII/MgII measurement changes as a function of spectral modeling. Quasars at even higher redshift or at fainter luminosity range (L_bol<10^46 erg s-1) are needed to probe the sites of early metal enrichment and a corresponding change in the FeII/MgII ratio.

Some metal-poor stars have abundance patterns which are midway between the slow (s) and rapid (r) neutron capture processes. We show that the helium shell of a fast rotating massive star experiencing a jet-like explosion undergoes two efficient neutron capture processes: one during stellar evolution and one during the explosion. It eventually provides a material whose chemical composition is midway between the s- and r-process. A low metallicity 40~$M_{\odot}$ model with an initial rotational velocity of $\sim 700$~km~s$^{-1}$ was computed from birth to pre-supernova with a nuclear network following the slow neutron capture process. A 2D hydrodynamic relativistic code was used to model a $E = 10^{52}$~erg relativistic jet-like explosion hitting the stellar mantle. The jet-induced nucleosynthesis was calculated in post-processing with a network of 1812 nuclei. During the star's life, heavy elements from $30 \lesssim Z \lesssim 82$ are produced thanks to an efficient s-process, which is boosted by rotation. At the end of evolution, the helium shell is largely enriched in trans-iron elements and in (unburnt) $^{22}$Ne, whose abundance is $\sim 20$ times higher than in a non-rotating model. During the explosion, the jet heats the helium shell up to $\sim 1.5$ GK. It efficiently activates ($\alpha,n$) reactions, such as $^{22}$Ne($\alpha,n$), and leads to a strong n-process with neutron densities of $\sim 10^{19} - 10^{20}$~cm$^{-3}$ during $0.1$~second. This has the effect of shifting the s-process pattern towards heavier elements (e.g. Eu). The resulting chemical pattern is consistent with the abundances of the carbon-enhanced metal-poor r/s star CS29528-028, provided the ejecta of the jet model is not homogeneously mixed. This is a new astrophysical site which can explain at least some of the metal-poor stars showing abundance patterns midway between the s- and r-process.

We present the results of CO(1-0) observations of the host galaxy of a Type I superluminous supernova (SLSN-I), SN2017egm, one of the closest SLSNe-I at z = 0.03063, by using the Atacama Large Millimeter/submillimeter Array. The molecular gas mass of the host galaxy is $M_{\rm gas} = (4.8 \pm 0.3) \times 10^9$ $M_{\odot}$, placing it on the sequence of normal star-forming galaxies in an $M_{\rm gas}$-star-formation rate (SFR) plane. The molecular hydrogen column density at the location of SN2017egm is higher than that of the Type II SN PTF10bgl, which is also located in the same host galaxy, and those of other Type II and Ia SNe located in different galaxies, suggesting that SLSNe-I have a preference for a dense molecular gas environment. On the other hand, the column density at the location of SN2017egm is comparable to those of Type Ibc SNe. The surface densities of molecular gas and the SFR at the location of SN2017egm are consistent with those of spatially resolved local star-forming galaxies and follow the Schmidt-Kennicutt relation. These facts suggest that SLSNe-I can occur in environments with the same star-formation mechanism as in normal star-forming galaxies.

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.

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.

The nature of the first massive stars may be inferred by investigating the origin of the extremely metal-poor (EMP) stars, likely formed from the ejecta of one or a few previous massive stars. We investigate the rotational properties of early massive stars by comparing the abundance patterns of EMP stars with rotating massive stellar models. Low metallicity 20 $M_{\odot}$ stellar models with initial rotation rates between 0 and $70~\%$ of the critical velocity are computed. Explosions with strong fallback are assumed. The ejected material is considered to fit individually the abundance patterns of 272 EMP stars with $-4<$ [Fe/H] $<-3$. With increasing initial rotation, the [C/H], [N/H], [O/H], [Na/H], [Mg/H] and [Al/H] ratios in the massive star ejecta are gradually increased. Among the 272 EMP stars considered, $\sim 40-50~\%$ are consistent with our models. About $60 - 70~\%$ of the CEMP star sample is reproduced against $\sim 20 - 30~\%$ for the C-normal EMP star sample. The CEMP stars are preferentially reproduced with a material coming from mid to fast rotating massive stars. The velocity distribution derived from the best massive star models increases from no rotation to fast rotation. The maximum is reached for massive stars having initial equatorial velocities of $\sim 550 - 640$ km~s$^{-1}$. Although subject to significant uncertainties, these results suggest that the rotational mixing operating in between the H-burning shell and the He-burning core of early massive stars played an important role in the early chemical enrichment of the Universe. The comparison of the velocity distribution derived from the best massive star models with velocity distributions of nearby OB stars suggests a greater amount of massive fast rotators in the early Universe. This may have important consequences for reionization or integrated light from high redshift galaxies.

Using a prototype of the Tomo-e Gozen wide-field CMOS mosaic camera, we acquire wide-field optical images at a cadence of 2 Hz and search them for transient sources of duration 1.5 to 11.5 seconds. Over the course of eight nights, our survey encompasses the equivalent of roughly two days on one square degree, to a fluence equivalent to a limiting magnitude about $V = 15.6$ in a 1-second exposure. After examining by eye the candidates identified by a software pipeline, we find no sources which meet all our criteria. We compute upper limits to the rate of optical transients consistent with our survey, and compare those to the rates expected and observed for representative sources of ephemeral optical light.

We report observations of a stellar occultation by the classical Kuiper belt object (50000) Quaoar occurred on 28 June 2019. A single-chord high-cadence (2 Hz) photometry dataset was obtained with the Tomo-e Gozen CMOS camera mounted on the 1.05 m Schmidt telescope at Kiso Observatory. The obtained ingress and egress data do not show any indication of atmospheric refraction and allow to set new $1\sigma$ and $3\sigma$ upper limits of 6 and 16 nbar, respectively, for the surface pressure of a pure methane atmosphere. These upper limits are lower than the saturation vapor pressure of methane at Quaoar's expected mean surface temperature ($T \sim 44$ K) and imply the absence of a $\sim$10 nbar-level global atmosphere formed by methane ice on Quaoar's surface.

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

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

We 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 investigate the molecular gas in, and star-formation properties of, the host galaxy (CGCG 137-068) of a mysterious transient, AT2018cow, at kpc and larger scales, using archival band-3 data from the Atacama Large Millimeter/submillimeter Array (ALMA). AT2018cow is the nearest Fast-Evolving Luminous Transient (FELT), and this is the very first study unveiling molecular-gas properties of FELTs. The achieved rms and beam size are 0.21 mJy beam$^{-1}$ at a velocity resolution of $40$ km s$^{-1}$ and $3".66\times2".71$ ($1.1~{\rm kpc} \times 0.8~{\rm kpc}$), respectively. CO($J$=1-0) emission is successfully detected. The total molecular gas mass inferred from the CO data is $(1.85\pm0.04)\times10^8$ M$_\odot$ with the Milky Way CO-to-H$_2$ conversion factor. The H$_2$ column density at the AT2018cow site is estimated to be $8.6\times10^{20}$ cm$^{-2}$. The ALMA data reveal that (1) CGCG 137-068 is a normal star-forming (SF) dwarf galaxy in terms of its molecular gas and star-formation properties and (2) AT2018cow is located between a CO peak and a blue star cluster. These properties suggest on-going star formation and favor the explosion of a massive star as the progenitor of AT2018cow. We also find that CGCG 137-068 has a solar or super-solar metallicity. If the metallicity of the other FELT hosts is not higher than average, then some property of SF dwarf galaxies other than metallicity may be related to FELTs.

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

We present 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 observational evidence that an aspherical supernova explosion could have occurred in the First stars in the early universe. Our results are based on the First determination of a Zn abundance in a Hubble Space Telescope/Cosmic Origins Spectrograph high-resolution UV spectrum of a hyper metal-poor (HMP) star, HE1327-2326, with [Fe/H](NLTE) = -5.2. We determine [Zn/Fe] = 0.80$\pm$0.25 from a UV Zn I line at 2138 detected at $3.4\sigma$. Yields of a 25M$_{\odot}$ aspherical supernova model with artificially modified densities exploding with E = 5x10$^{51}$ ergs best match the entire abundance pattern of HE1327-2326. Such high-entropy hypernova explosions are expected to produce bipolar outfows which could facilitate the external enrichment of small neighboring galaxies. This has already been predicted by theoretical studies of the earliest star forming minihalos. Such a scenario would have significant implications for the chemical enrichment across the early universe as HMP Carbon Enhanced Metal-Poor (CEMP) stars such as HE1327-2326 might have formed in such externally enriched environments.

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

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

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.

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

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.

Wide-field optical surveys have begun to uncover large samples of fast (t_rise < 5d), luminous (M_peak < -18), blue transients. While commonly attributed to the breakout of a supernova shock into a dense wind, the great distances to the transients of this class found so far have hampered detailed investigation of their properties. We present photometry and spectroscopy from a comprehensive worldwide campaign to observe AT2018cow (ATLAS18qqn), the first fast-luminous optical transient to be found in real time at low redshift. Our first spectra (<2 days after discovery) are entirely featureless. A very broad absorption feature suggestive of near-relativistic velocities develops between 3-8 days, then disappears. Broad emission features of H and He develop after >10 days. The spectrum remains extremely hot throughout its evolution, and the photospheric radius contracts with time (receding below R<10^14 cm after 1 month). This behaviour does not match that of any known supernova, although a relativistic jet within a fallback supernova could explain some of the observed features. Alternatively, the transient could originate from the disruption of a star by an intermediate-mass black hole, although this would require long-lasting emission of highly super-Eddington thermal radiation. In either case, AT2018cow suggests that the population of fast luminous transients represents a new class of astrophysical event. Intensive follow-up of this event in its late phases, and of any future events found at comparable distance, will be essential to better constrain their origins.

We present optical and near-infrared observations of the nearby Type Iax supernova (SN) 2014dt from 14 to 410 days after the maximum light. The velocities of the iron absorption lines in the early phase indicated that SN 2014dt showed slower expansion than the well-observed Type Iax SNe 2002cx, 2005hk and 2012Z. In the late phase, the evolution of the light curve and that of the spectra were considerably slower. The spectral energy distribution kept roughly the same shape after ~100 days, and the bolometric light curve flattened during the same period. These observations suggest the existence of an optically thick component that almost fully trapped the {\gamma}-ray energy from 56 Co decay. These findings are consistent with the predictions of the weak deflagration model, leaving a bound white dwarf remnant after the explosion.

Imaging observations of faint meteors were carried out on April 11 and 14, 2016 with a wide-field CMOS mosaic camera, Tomo-e PM, mounted on the 105-cm Schmidt telescope at Kiso Observatory, the University of Tokyo. Tomo-e PM, which is a prototype model of Tomo-e Gozen, can monitor a sky of ${\sim}1.98\,\mathrm{deg^2}$ at 2\,Hz. The numbers of detected meteors are 1514 and 706 on April 11 and 14, respectively. The detected meteors are attributed to sporadic meteors. Their absolute magnitudes range from $+4$ to $+10\,\mathrm{mag}$ in the $V$-band, corresponding to about $8.3{\times}10^{-2}$ to $3.3{\times}10^{-4}\,\mathrm{g}$ in mass. The present magnitude distributions we obtained are well explained by a single power-law luminosity function with a slope parameter $r = 3.1{\pm}0.4$ and a meteor rate $\log_{10}N_0 = -5.5{\pm}0.5$. The results demonstrate a high performance of telescopic observations with a wide-field video camera to constrain the luminosity function of faint meteors. The performance of Tomo-e Gozen is about two times higher than that of Tomo-e PM. A survey with Tomo-e Gozen will provide a more robust measurement of the luminosity function.

We have searched for optical variability in the host galaxy of the radio variable source which is possibly associated with fast radio burst (FRB) 150418. We compare images of the galaxy taken 1 day after the burst using Subaru/Suprime-Cam with images taken $\sim$ 1 year after the burst using Gemini-South/GMOS. No optical variability is found between the two epochs with a limiting absolute magnitude $\gtrsim -19$ (AB). This limit applies to optical variability of the putative active galactic nucleus in the galaxy on a timescale of $\sim$ 1 year, and also to the luminosity of an optical counterpart of FRB~150418 one day after the burst should it have occurred in this galaxy.

The origin of fast radio bursts (FRBs), bright millisecond radio transients, is still somewhat of a mystery. Several theoretical models expect that the FRB accompanies an optical afterglow (e.g., Totani 2013; Kashiyama et al. 2013). In order to investigate the origin of FRBs, we perform $gri$-band follow-up observations of FRB~151230 (estimated $z \lesssim 0.8$) with Subaru/Hyper Suprime-Cam at $8$, $11$, and $14$~days after discovery. The follow-up observation reaches a $50\%$ completeness magnitude of $26.5$~mag for point sources, which is the deepest optical follow-up of FRBs to date. We find $13$ counterpart candidates with variabilities during the observation. We investigate their properties with multicolor and multi-wavelength observations and archival catalogs. Two candidates are excluded by the non-detection of FRB~151230 in the other radio feed horns that operated simultaneously to the detection, as well as the inconsistency between the photometric redshift and that derived from the dispersion measure of FRB~151230. Eight further candidates are consistent with optical variability seen in AGNs. Two more candidates are well fitted with transient templates (Type IIn supernovae), and the final candidate is poorly fitted with all of our transient templates and is located off-center of an extended source. It can only be reproduced with rapid transients with a faint peak and rapid decline and the probability of chance coincidence is $\sim3.6\%$. We also find that none of our candidates are consistent with Type Ia supernovae, which rules out the association of Type Ia supernovae to FRB~151230 at $z\leq0.6$ and limits the dispersion measure of the host galaxy to $\lesssim300$~pc~cm$^{-3}$ in a Type Ia supernova scenario.

Extremely metal-poor (EMP) stars in the Milky Way (MW) allow us to infer the properties of their progenitors by comparing their chemical composition to the metal yields of the first supernovae. This method is most powerful when applied to mono-enriched stars, i.e. stars that formed from gas that was enriched by only one previous supernova. We present a novel diagnostic to identify this subclass of EMP stars. We model the first generations of star formation semi-analytically, based on dark matter halo merger trees that yield MW-like halos at the present day. Radiative and chemical feedback are included self-consistently and we trace all elements up to zinc. Mono-enriched stars account for only $\sim 1\%$ of second generation stars in our fiducial model and we provide an analytical formula for this probability. We also present a novel analytical diagnostic to identify mono-enriched stars, based on the metal yields of the first supernovae. This new diagnostic allows us to derive our main results independently from the specific assumptions made regarding Pop III star formation, and we apply it to a set of observed EMP stars to demonstrate its strengths and limitations. Our results may provide selection criteria for current and future surveys and therefore contribute to a deeper understanding of EMP stars and their progenitors.

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.

A pulsar wind is a relativistic outflow dominated by Poynting energy at its base. Based on the standard ideal magnetohydrodynamic (MHD) model of pulsar wind nebulae (PWNe) with the ordered magnetic field, the observed slow expansion $v_{\rm PWN} \ll c$ requires the wind to be dominated by kinetic energy at the upstream of its termination shock, which conflicts with the pulsar wind theory ($\sigma$-problem). In this paper, we extend the standard model of PWNe by phenomenologically taking into account conversion of the ordered to turbulent magnetic field and dissipation of the turbulent magnetic field. Disordering of the magnetic structure is inferred from the recent three-dimensional relativistic ideal MHD simulations, while magnetic dissipation is a non-ideal MHD effect requiring a finite resistivity. We apply this model to the Crab Nebula and find that the conversion effect is important for the flow deceleration, while the dissipation effect is not. Even for Poynting-dominated pulsar wind, we obtain the Crab Nebula's $v_{\rm PWN}$ by adopting a finite conversion time-scale of $\sim 0.3$ yr. Magnetic dissipation primarily affects the synchrotron radiation properties. Any values of the pulsar wind magnetization $\sigma_{\rm w}$ are allowed within the present model of the PWN dynamics alone, and even a small termination shock radius of $\ll 0.1$ pc reproduces the observed dynamical features of the Crab Nebula. In order to establish a high-$\sigma_{\rm w}$ model of PWNe, it is important to extend the present model by taking into account the broadband spectrum and its spacial profiles.

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.

We compare elemental abundance patterns of $\sim 200$ extremely metal-poor (EMP; [Fe/H]$<-3$) stars with supernova yields of metal-free stars in order to obtain insights into the characteristic masses of the first (Population III or Pop III) stars in the Universe. Supernova yields are prepared with nucleosynthesis calculations of metal-free stars with various initial masses ($M=$13, 15, 25, 40 and 100 $M_{\odot}$) and explosion energies ($E_{51}=E/10^{51}$[erg]$=0.5-60$) to include low-energy, normal-energy, and high-energy explosions. We adopt the mixing-fallback model to take into account possible asymmetry in the supernova explosions and the yields that best-fit the observed abundance patterns of the EMP stars are searched by varying the model parameters. We find that the abundance patterns of the EMP stars are predominantly best-fitted with the supernova yields with initial masses $M<40 M_{\odot}$, and that more than than half of the stars are best fitted with the $M=25 M_\odot$ hypernova ($E_{51}=10$) models. The results also indicate that the majority of the primordial supernovae have ejected $10^{-2}-10^{-1} M_\odot$ of $^{56}$Ni leaving behind a compact remnant, either a neutron star or a black hole, with mass in a range of $\sim 1.5-5 M_{\odot}$. The results suggest that the masses of the first stars responsible for the first metal-enrichment are predominantly $< 40 M_{\odot}$. This implies that the higher mass first stars were either less abundant or directly collapsing into a blackhole without ejecting heavy elements or that a supernova explosion of a higher-mass first star inhibits the formation of the next generation of low-mass stars at [Fe/H]$<-3$.

We present the results of 3 GHz radio continuum observations of the 8 host galaxies of super-luminous supernovae (SLSNe) at $0.1 < z < 0.3$ by using the Karl G. Jansky Very Large Array. Four host galaxies are detected significantly, and two of them are found to have high star-formation rates (SFRs $>$ 20 $M_{\odot}$ yr$^{-1}$) derived from radio emission, making them the most intensely star-forming host galaxies among SLSN host galaxies. We compare radio SFRs and optical SFRs, and find that three host galaxies have an excess in radio SFRs by a factor of $>$2, suggesting the existence of dust-obscured star formation, which cannot be traced by optical studies. Two of the three host galaxies, which are located in the galaxy main sequence based on optical SFRs, are found to be above the main sequence based on their radio SFRs. This suggests a higher fraction of starburst galaxies in SLSN hosts than estimated in previous studies. We calculate extinction from the ratio between radio SFRs and dust-uncorrected optical SFRs and find that the hosts are on the trend of increasing extinction with metallicity, which is consistent with the relation in local star-forming galaxies. We also place a constraint on a pulsar-driven SN model, which predicts quasi-steady synchrotron radio emission.

The birth properties of neutron stars yield important information on the still debated physical processes that trigger the explosion and on intrinsic neutron-star physics. These properties include the high space velocities of young neutron stars with average values of several 100 km/s, whose underlying "kick" mechanism is not finally clarified. There are two competing possibilities that could accelerate NSs during their birth: anisotropic ejection of either stellar debris or neutrinos. We here present new evidence from X-ray measurements that chemical elements between silicon and calcium in six young gaseous supernova remnants are preferentially expelled opposite to the direction of neutron star motion. There is no correlation between the kick velocities and magnetic field strengths of these neutron stars. Our results support a hydrodynamic origin of neutron-star kicks connected to asymmetric explosive mass ejection, and they conflict with neutron-star acceleration scenarios that invoke anisotropic neutrino emission caused by particle and nuclear physics in combination with very strong neutron-star magnetic fields.

We perform a $z$-band survey for an optical counterpart of a binary neutron star coalescence GW170817 with Subaru/Hyper Suprime-Cam. Our untargeted transient search covers $23.6$ deg$^2$ corresponding to the $56.6\%$ credible region of GW170817 and reaches the $50\%$ completeness magnitude of $20.6$ mag on average. As a result, we find 60 candidates of extragalactic transients, including J-GEM17btc (a.k.a. SSS17a/DLT17ck). While J-GEM17btc is associated with NGC 4993 that is firmly located inside the 3D skymap of GW170817, the other 59 candidates do not have distance information in the GLADE v2 catalog or NASA/IPAC Extragalactic Database (NED). Among 59 candidates, 58 are located at the center of extended objects in the Pan-STARRS1 catalog, while one candidate has an offset. We present location, $z$-band apparent magnitude, and time variability of the candidates and evaluate the probabilities that they are located inside of the 3D skymap of GW170817. The probability for J-GEM17btc is $64\%$ being much higher than those for the other 59 candidates ($9.3\times10^{-3}-2.1\times10^{-1}\%$). Furthermore, the possibility, that at least one of the other 59 candidates is located within the 3D skymap, is only $3.2\%$. Therefore, we conclude that J-GEM17btc is the most-likely and distinguished candidate as the optical counterpart of GW170817.

We report the discovery of four Fast Radio Bursts (FRBs) in the ongoing SUrvey for Pulsars and Extragalactic Radio Bursts (SUPERB) at the Parkes Radio Telescope: FRBs 150610, 151206, 151230 and 160102. Our real-time discoveries have enabled us to conduct extensive, rapid multi-messenger follow-up at 12 major facilities sensitive to radio, optical, X-ray, gamma-ray photons and neutrinos on time scales ranging from an hour to a few months post-burst. No counterparts to the FRBs were found and we provide upper limits on afterglow luminosities. None of the FRBs were seen to repeat. Formal fits to all FRBs show hints of scattering while their intrinsic widths are unresolved in time. FRB 151206 is at low Galactic latitude, FRB 151230 shows a sharp spectral cutoff, and FRB 160102 has the highest dispersion measure (DM = $2596.1\pm0.3$ pc cm$^{-3}$) detected to date. Three of the FRBs have high dispersion measures (DM >$1500$ pc cm$^{-3}$), favouring a scenario where the DM is dominated by contributions from the Intergalactic Medium. The slope of the Parkes FRB source counts distribution with fluences $>2$ Jyms is $\alpha=-2.2^{+0.6}_{-1.2}$ and still consistent with a Euclidean distribution ($\alpha=-3/2$). We also find that the all-sky rate is $1.7^{+1.5}_{-0.9}\times10^3$FRBs/($4\pi$ sr)/day above $\sim2$ Jyms and there is currently no strong evidence for a latitude-dependent FRB sky-rate.

The first detected gravitational wave from a neutron star merger was GW170817. In this study, we present J-GEM follow-up observations of SSS17a, an electromagnetic counterpart of GW170817. SSS17a shows a 2.5-mag decline in the $z$-band from 1.7 days to 7.7 days after the merger. Such a rapid decline is not comparable with supernovae light curves at any epoch. The color of SSS17a also evolves rapidly and becomes redder for later epochs; the $z-H$ color changed by approximately 2.5 mag in the period of 0.7 days to 7.7 days. The rapid evolution of both the optical brightness and the color are consistent with the expected properties of a kilonova that is powered by the radioactive decay of newly synthesized $r$-process nuclei. Kilonova models with Lanthanide elements can reproduce the aforementioned observed properties well, which suggests that $r$-process nucleosynthesis beyond the second peak takes place in SSS17a. However, the absolute magnitude of SSS17a is brighter than the expected brightness of the kilonova models with the ejecta mass of 0.01 $\Msun$, which suggests a more intense mass ejection ($\sim 0.03 \Msun$) or possibly an additional energy source.

Recent detection of gravitational waves from a neutron star (NS) merger event GW170817 and identification of an electromagnetic counterpart provide a unique opportunity to study the physical processes in NS mergers. To derive properties of ejected material from the NS merger, we perform radiative transfer simulations of kilonova, optical and near-infrared emissions powered by radioactive decays of r-process nuclei synthesized in the merger. We find that the observed near-infrared emission lasting for > 10 days is explained by 0.03 Msun of ejecta containing lanthanide elements. However, the blue optical component observed at the initial phases requires an ejecta component with a relatively high electron fraction (Ye). We show that both optical and near-infrared emissions are simultaneously reproduced by the ejecta with a medium Ye of ~ 0.25. We suggest that a dominant component powering the emission is post-merger ejecta, which exhibits that mass ejection after the first dynamical ejection is quite efficient. Our results indicate that NS mergers synthesize a wide range of r-process elements and strengthen the hypothesis that NS mergers are the origin of r-process elements in the Universe.

Extremely metal-poor (EMP) stars are the living fossils with records of chemical enrichment history at the early epoch of galaxy formation. By the recent large observation campaigns, statistical samples of EMP stars have been obtained. This motivates us to reconsider their classification and formation conditions. From the observed lower-limits of carbon and iron abundances of Acr(C) ~ 6 and [Fe/H]cr ~ -5 for C-enhanced EMP (CE-EMP) and C-normal EMP (CN-EMP) stars, we confirm that gas cooling by dust thermal emission is indispensable for the fragmentation of their parent clouds to form such low-mass, i.e., long-lived stars, and that the dominant grain species are carbon and silicate, respectively. We constrain the grain radius r_i^{cool} of a species i and condensation efficiency f_{ij} of a key element j as r_C^{cool}/f_{C, C} = 10 um and r_Sil^{cool}/f_{Sil, Mg} = 0.1 um to reproduce Acr(C) and [Fe/H]cr, which give a universal condition 10^{[C/H]-2.30} + 10^[Fe/H] > 10^{-5.07} for the formation of every EMP star. Instead of the conventional boundary [C/Fe] = 0.7 between CE- and CN-EMP stars, this condition suggests a physically meaningful boundary [C/Fe]_{b} = 2.30 above and below which carbon and silicate grains are dominant coolants, respectively.

Type Ia supernovae (SNe Ia) arise from the thermonuclear explosion of carbon-oxygen white dwarfs. Though the uniformity of their light curves makes them powerful cosmological distance indicators, long-standing issues remain regarding their progenitors and explosion mechanisms. Recent detection of the early ultraviolet pulse of a peculiar subluminous SN Ia has been claimed as new evidence for the companion-ejecta interaction through the single-degenerate channel. Here, we report the discovery of a prominent but red optical flash at $\sim$ 0.5 days after the explosion of a SN Ia which shows hybrid features of different SN Ia sub-classes: a light curve typical of normal-brightness SNe Ia, but with strong titanium absorptions, commonly seen in the spectra of subluminous ones. We argue that the early flash of such a hybrid SN Ia is different from predictions of previously suggested scenarios such as the companion-ejecta interaction. Instead it can be naturally explained by a SN explosion triggered by a detonation of a thin helium shell either on a near-Chandrasekhar-mass white dwarf ($\gtrsim$ 1.3 M$_{\odot}$) with low-yield $^{56}$Ni or on a sub-Chandrasekhar-mass white dwarf ($\sim$ 1.0 M$_{\odot}$) merging with a less massive white dwarf. This finding provides compelling evidence that one branch of the previously proposed explosion models, the helium-ignition scenario, does exist in nature, and such a scenario may account for explosions of white dwarfs in a wider mass range in contrast to what was previously supposed.

We present the results of the detailed analysis of an optical imaging survey conducted using the Subaru / Hyper Suprime-Cam (HSC), which aims to identify an optical counterpart to the gravitational wave event GW151226. In half a night, the $i$- and $z$-band imaging survey by HSC covers 63.5deg$^2$ of the error region, which contains about 7\% of the LIGO localization probability, and the same field is observed in three different epochs. The detectable magnitude of the candidates in a differenced image is evaluated as $i \sim 23.2$ mag for the requirement of at least two 5$\sigma$ detections, and 1744 candidates are discovered. Assuming a kilonova as an optical counterpart, we compared the optical properties of the candidates with model predictions. A red and rapidly declining light curve condition enables the discrimination of a kilonova from other transients, and a small number of candidates satisfy this condition. The presence of stellar-like counterparts in the reference frame suggests that the surviving candidates are likely to be flare stars. The fact that most of those candidates are in galactic plane, $|b|<5^{\circ}$, supports this interpretation. We also checked whether the candidates are associated with the nearby GLADE galaxies, which reduces the number of contaminants even with a looser color cut. When a better probability map (with localization accuracy of $\sim50{\rm deg}^2$) is available, kilonova searches of up to approximately $200$ Mpc will become feasible by conducting immediate follow-up observations with an interval of 3--6 days.

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.

We present monitoring campaign observations at optical and near-infrared (NIR) wavelengths for a radio-loud active galactic nucleus (AGN) at z=0.840, SDSS~J110006.07+442144.3 (hereafter, J1100+4421), which was identified during a flare phase in late February, 2014. The campaigns consist of three intensive observing runs from the discovery to March, 2015, mostly within the scheme of the OISTER collaboration. Optical-NIR light curves and simultaneous spectral energy distributions (SEDs) are obtained. Our measurements show the strongest brightening in March, 2015. We found that the optical-NIR SEDs of J1100+4421 show an almost steady shape despite the large and rapid intranight variability. This constant SED shape is confirmed to extend to $\sim5~\mu$m in the observed frame using the archival WISE data. Given the lack of absorption lines and the steep power-law spectrum of $\alpha_{\nu}\sim-1.4$, where $f_{\nu}\propto\nu^{\alpha_{\nu}}$, synchrotron radiation by a relativistic jet with no or small contributions from the host galaxy and the accretion disk seems most plausible as an optical-NIR emission mechanism. The steep optical-NIR spectral shape and the large amplitude of variability are consistent with this object being a low $\nu_{\rm{peak}}$ jet-dominated AGN. In addition, sub-arcsec resolution optical imaging data taken with Subaru Hyper Suprime-Cam does not show a clear extended component and the spatial scales are significantly smaller than the large extensions detected at radio wavelengths. The optical spectrum of a possible faint companion galaxy does not show any emission lines at the same redshift and hence a merging hypothesis for this AGN-related activity is not supported by our observations.

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.

Low-mass population III (PopIII) stars of $\lesssim 0.8 M_{\odot}$ could survive up until the present. Non-detection of low-mass PopIII stars in our Galaxy has already put a stringent constraint on the initial mass function (IMF) of PopIII stars, suggesting that PopIII stars have a top-heavy IMF. On the other hand, some claims that the lack of such stars stems from metal enrichment of their surface by accretion of heavy elements from interstellar medium (ISM). We investigate effects of the stellar wind on the metal accretion onto low-mass PopIII stars because accretion of the local ISM onto the Sun is prevented by the solar wind even for neutrals. The stellar wind and radiation of low-mass PopIII stars are modeled based on knowledge of nearby low-mass stellar systems including our Sun. We find that low-mass PopIII stars traveling across the Galaxy forms the stellar magnetosphere in most of their life. Once the magnetosphere is formed, most of neutral interstellar particles are photoionized before reaching to the stellar surface and are blown away by the wind. Especially, the accretion abundance of iron will be reduced by a factor of $< 10^{-12}$ compared with Bondi-Hoyle-Lyttleton accretion. The metal accretion can enhance iron abundance [Fe/H] only up to $\sim -14$. This demonstrates that low-mass PopIII stars remain pristine and will be found as metal free stars and that further searches for them are valuable to constrain the IMF of PopIII stars.

$^{56}$Ni is an important indicator of the supernova explosions, which characterizes light curves. Nevertheless, rather than $^{56}$Ni, the explosion energy has often been paid attention from the explosion mechanism community, since it is easier to estimate from numerical data than the amount of $^{56}$Ni. The final explosion energy, however, is difficult to estimate by detailed numerical simulations because current simulations cannot reach typical timescale of saturation of explosion energy. Instead, the amount of $^{56}$Ni converges within a short timescale so that it would be a better probe of the explosion mechanism. We investigated the amount of $^{56}$Ni synthesized by explosive nucleosynthesis in supernova ejecta by means of numerical simulations and an analytic model. For numerical simulations, we employ Lagrangian hydrodynamics code in which neutrino heating and cooling terms are taken into account by light-bulb approximation. Initial conditions are taken from Woosley & Hegel (2007), which have 12, 15, 20, and 25 $M_\odot$ in zero age main sequence. We additionally develop an analytic model, which gives a reasonable estimate of the amount of $^{56}$Ni. We found that, in order to produce enough amount of $^{56}$Ni, $\mathcal{O}(1)$ Bethe s$^{-1}$ of growth rate of the explosion energy is needed, which is much larger than that found in recent exploding simulations, typically $\mathcal{O}(0.1)$ Bethe s$^{-1}$.

We report on the detection of a remarkable new fast high-energy transient found in the Chandra Deep Field-South, robustly associated with a faint ($m_{\rm R}=27.5$ mag, $z_{\rm ph}$$\sim$2.2) host in the CANDELS survey. The X-ray event is comprised of 115$^{+12}_{-11}$ net 0.3-7.0 keV counts, with a light curve characterised by a $\approx$100 s rise time, a peak 0.3-10 keV flux of $\approx$5$\times$10$^{-12}$ erg s$^{-1}$ cm$^{-2}$, and a power-law decay time slope of $-1.53\pm0.27$. The average spectral slope is $\Gamma=1.43^{+0.23}_{-0.13}$, with no clear spectral variations. The \hbox{X-ray} and multi-wavelength properties effectively rule out the vast majority of previously observed high-energy transients. A few theoretical possibilities remain: an "orphan" X-ray afterglow from an off-axis short-duration Gamma-ray Burst (GRB) with weak optical emission; a low-luminosity GRB at high redshift with no prompt emission below $\sim$20 keV rest-frame; or a highly beamed Tidal Disruption Event (TDE) involving an intermediate-mass black hole and a white dwarf with little variability. However, none of the above scenarios can completely explain all observed properties. Although large uncertainties exist, the implied rate of such events is comparable to those of orphan and low-luminosity GRBs as well as rare TDEs, implying the discovery of an untapped regime for a known transient class, or a new type of variable phenomena whose nature remains to be determined.

We present the early-phase spectra and the light curves of the broad-lined supernova (SN) 2016coi from $t=7$ to $67$ days after the estimated explosion date. This SN was initially reported as a broad-lined Type SN Ic (SN Ic-BL). However, we found that spectra up to $t=12$ days exhibited the He~{\sc i} $\lambda$5876, $\lambda$6678, and $\lambda$7065 absorption lines. We show that the smoothed and blueshifted spectra of normal SNe Ib are remarkably similar to the observed spectrum of SN 2016coi. The line velocities of SN 2016coi were similar to those of SNe Ic-BL and significantly faster than those of SNe Ib. Analyses of the line velocity and light curve suggest that the kinetic energy and the total ejecta mass of SN 2016coi are similar to those of SNe Ic-BL. Together with broad-lined SNe 2009bb and 2012ap for which the detection of He~{\sc i} were also reported, these SNe could be transitional objects between SNe Ic-BL and SNe Ib, and be classified as broad-lined Type `Ib' SNe (SNe `Ib'-BL). Our work demonstrates the diversity of the outermost layer in broad-lined SNe, which should be related to the variety of the evolutionary paths.

We study the number and the distribution of low mass Pop III stars in the Milky Way. In our numerical model, hierarchical formation of dark matter minihalos and Milky Way sized halos are followed by a high resolution cosmological simulation. We model the Pop III formation in H2 cooling minihalos without metal under UV radiation of the Lyman-Werner bands. Assuming a Kroupa IMF from 0.15 to 1.0 Msun for low mass Pop III stars, as a working hypothesis, we try to constrain the theoretical models in reverse by current and future observations. We find that the survivors tend to concentrate on the center of halo and subhalos. We also evaluate the observability of Pop III survivors in the Milky Way and dwarf galaxies, and constraints on the number of Pop III survivors per minihalo. The higher latitude fields require lower sample sizes because of the high number density of stars in the galactic disk, the required sample sizes are comparable in the high and middle latitude fields by photometrically selecting low metallicity stars with optimized narrow band filters, and the required number of dwarf galaxies to find one Pop III survivor is less than ten at <100 kpc for the tip of redgiant stars. Provided that available observations have not detected any survivors, the formation models of low mass Pop III stars with more than ten stars per minihalo are already excluded. Furthermore, we discuss the way to constrain the IMF of Pop III star at a high mass range of > 10 Msun.

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 report the results of optical--infrared follow-up observations of the gravitational wave (GW) event GW151226 detected by the Advanced LIGO in the framework of J-GEM (Japanese collaboration for Gravitational wave ElectroMagnetic follow-up). We performed wide-field optical imaging surveys with Kiso Wide Field Camera (KWFC), Hyper Suprime-Cam (HSC), and MOA-cam3. The KWFC survey started at 2.26 days after the GW event and covered 778 deg$^2$ centered at the high Galactic region of the skymap of GW151226. We started the HSC follow-up observations from 12 days after the event and covered an area of 63.5 deg$^2$ of the highest probability region of the northern sky with the limiting magnitudes of 24.6 and 23.8 for i band and z band, respectively. MOA-cam3 covered 145 deg$^2$ of the skymap with MOA-red filter 2.5 months after the GW alert. Total area covered by the wide-field surveys was 986.5 deg$^2$. The integrated detection probability of all the observed area was $\sim$29%. We also performed galaxy-targeted observations with six optical and near-infrared telescopes from 1.61 days after the event. Total of 238 nearby (<100 Mpc) galaxies were observed with the typical I band limiting magnitude of $\sim$19.5. We detected 13 supernova candidates with the KWFC survey, and 60 extragalactic transients with the HSC survey. Two third of the HSC transients were likely supernovae and the remaining one third were possible active galactic nuclei. With our observational campaign, we found no transients that are likely to be associated with GW151226.

We present an application of machine-learning (ML) techniques to source selection in the optical transient survey data with Hyper Suprime-Cam (HSC) on the Subaru telescope. Our goal is to select real transient events accurately and in a timely manner out of a large number of false candidates, obtained with the standard difference-imaging method. We have developed the transient selector which is based on majority voting of three ML machines of AUC Boosting, Random Forest, and Deep Neural Network. We applied it to our observing runs of Subaru-HSC in 2015 May and August, and proved it to be efficient in selecting optical transients. The false positive rate was 1.0% at the true positive rate of 90% in the magnitude range of 22.0--25.0 mag for the former data. For the latter run, we successfully detected and reported ten candidates of supernovae within the same day as the observation. From these runs, we learned the following lessons: (1) the training using artificial objects is effective in filtering out false candidates, especially for faint objects, and (2) combination of ML by majority voting is advantageous.

A fading radio source, coincident in time and position with the fast radio burst FRB150418, has been associated with the galaxy WISE J071634.59-190039.2. Subsequent observations of this galaxy have revealed that it contains a persistent, but variable, radio source. We present e-MERLIN, VLBA, and ATCA radio observations and Subaru optical observations of WISE J071634.59-190039.2 and find that the persistent radio source is unresolved and must be compact (<0.01 kpc), and that its location is consistent with the optical centre of the galaxy. We conclude that it is likely that WISE J071634.59-190039.2 contains a weak radio AGN.

We present our optical follow-up observations to search for an electromagnetic counterpart of the first gravitational wave source GW150914 in the framework of the Japanese collaboration for Gravitational wave ElectroMagnetic follow-up (J-GEM), which is an observing group utilizing optical and radio telescopes in Japan, as well as those in New Zealand, China, South Africa, Chile, and Hawaii. We carried out a wide-field imaging survey with Kiso Wide Field Camera (KWFC) on the 1.05-m Kiso Schmidt telescope in Japan and a galaxy-targeted survey with Tripole5 on the B&C 61-cm telescope in New Zealand. Approximately 24 deg2 regions in total were surveyed in i-band with KWFC and 18 nearby galaxies were observed with Tripole5 in g-, r-, and i-bands 4-12 days after the gravitational wave detection. Median 5-sigma depths are i~18.9 mag for the KWFC data and g~18.9 mag, r~18.7 mag, and i~18.3 mag for the Tripole5 data. Probability for a counterpart to be in the observed area is 1.2% in the initial skymap and 0.1% in the final skymap. We do not find any transient source associated to an external galaxy with spatial offset from its center, which is consistent with the local supernova rate. We summarize future prospects and ongoing efforts to pin down electromagnetic counterparts of binary black hole mergers as well as neutron star mergers.

We present extensively dense observations of the super-Chandrasekhar supernova (SC SN) candidate SN 2012dn from $-11$ to $+140$ days after the date of its $B$-band maximum in the optical and near-infrared (NIR) wavelengths conducted through the OISTER ToO program. The NIR light curves and color evolutions up to 35 days after the $B$-band maximum provided an excellent match with those of another SC SN 2009dc, providing a further support to the nature of SN 2012dn as a SC SN. We found that SN 2012dn exhibited strong excesses in the NIR wavelengths from $30$ days after the $B$-band maximum. The $H$ and $K_{s}$-band light curves exhibited much later maximum dates at $40$ and $70$ days after the $B$-band maximum, respectively, compared with those of normal SNe Ia. The $H$ and $K_{s}$-band light curves subtracted by those of SN 2009dc displayed plateaued evolutions, indicating a NIR echo from the surrounding dust. The distance to the inner boundary of the dust shell is limited to be $4.8 - 6.4\times10^{-2}$ pc. No emission lines were found in its early phase spectrum, suggesting that the ejecta-CSM interaction could not occur. On the other hand, we found no signature that strongly supports the scenario of dust formation. The mass loss rate of the pre-explosion system is estimated to be $10^{-6}-10^{-5}$ M$_{\odot}$ yr$^{-1}$, assuming that the wind velocity of the system is $10-100$ km~s$^{-1}$, which suggests that the progenitor of SN 2012dn could be a recurrent nova system. We conclude that the progenitor of this SC SN could be explained by the single-degenerate scenario.

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

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

In recent years, millisecond duration radio signals originating from distant galaxies appear to have been discovered in the so-called Fast Radio Bursts. These signals are dispersed according to a precise physical law and this dispersion is a key observable quantity which, in tandem with a redshift measurement, can be used for fundamental physical investigations. While every fast radio burst has a dispersion measurement, none before now have had a redshift measurement, due to the difficulty in pinpointing their celestial coordinates. Here we present the discovery of a fast radio burst and the identification of a fading radio transient lasting $\sim 6$ days after the event, which we use to identify the host galaxy; we measure the galaxy's redshift to be $z=0.492\pm0.008$. The dispersion measure and redshift, in combination, provide a direct measurement of the cosmic density of ionised baryons in the intergalactic medium of $\Omega_{\mathrm{IGM}}=4.9 \pm 1.3\%$, in agreement with the expectation from WMAP, and including all of the so-called "missing baryons". The $\sim6$-day transient is largely consistent with a short gamma-ray burst radio afterglow, and its existence and timescale do not support progenitor models such as giant pulses from pulsars, and supernovae. This contrasts with the interpretation of another recently discovered fast radio burst, suggesting there are at least two classes of bursts.

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 report the first detection of thermal X-ray line emission from the supernova remnant (SNR) RX J1713.7-3946, the prototype of the small class of synchrotron dominated SNRs. A softness-ratio map generated using XMM-Newton data shows that faint interior regions are softer than bright shell regions. Using Suzaku and deep XMM-Newton observations, we have extracted X-ray spectra from the softest area, finding clear line features at 1 keV and 1.35 keV. These lines can be best explained as Ne Ly-alpha and Mg He-alpha from a thermal emission component. Since the abundance ratios of metals to Fe are much higher than solar values in the thermal component, we attribute the thermal emission to reverse-shocked SN ejecta. The measured Mg/Ne, Si/Ne, and Fe/Ne ratios of 2.0-2.6, 1.5-2.0, and <0.05 solar suggest that the progenitor star of RX J1713.7-3946 was a relatively low-mass star (<~20 M_sun), consistent with a previous inference based on the effect of stellar winds of the progenitor star on the surrounding medium. Since the mean blastwave speed of ~6000 km/s (the radius of 9.6 pc divided by the age of 1600 yr) is relatively fast compared with other core-collapse SNRs, we propose that RX J1713.7-3946 is a result of a Type Ib/c supernova whose progenitor was a member of an interacting binary. While our analysis provides strong evidence for X-ray line emission, our interpretation of its nature as thermal emission from SN ejecta requires further confirmation especially through future precision spectroscopic measurements using ASTRO-H.

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

The TMT Detailed Science Case describes the transformational science that the Thirty Meter Telescope will enable. Planned to begin science operations in 2024, TMT will open up opportunities for revolutionary discoveries in essentially every field of astronomy, astrophysics and cosmology, seeing much fainter objects much more clearly than existing telescopes. Per this capability, TMT's science agenda fills all of space and time, from nearby comets and asteroids, to exoplanets, to the most distant galaxies, and all the way back to the very first sources of light in the Universe. More than 150 astronomers from within the TMT partnership and beyond offered input in compiling the new 2015 Detailed Science Case. The contributing astronomers represent the entire TMT partnership, including the California Institute of Technology (Caltech), the Indian Institute of Astrophysics (IIA), the National Astronomical Observatories of the Chinese Academy of Sciences (NAOC), the National Astronomical Observatory of Japan (NAOJ), the University of California, the Association of Canadian Universities for Research in Astronomy (ACURA) and US associate partner, the Association of Universities for Research in Astronomy (AURA).

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

We report observations of the Type Iax supernova (SN Iax) 2012Z at optical and near-infrared wavelengths from immediately after the explosion until $\sim$ $260$ days after the maximum luminosity using the Optical and Infrared Synergetic Telescopes for Education and Research (OISTER) Target-of-Opportunity (ToO) program and the Subaru telescope. We found that the near-infrared (NIR) light curve evolutions and color evolutions are similar to those of SNe Iax 2005hk and 2008ha. The NIR absolute magnitudes ($M_{J}\sim-18.1$ mag and $M_{H}\sim-18.3$ mag) and the rate of decline of the light curve ($\Delta$ $m_{15}$($B$)$=1.6 \pm 0.1$ mag) are very similar to those of SN 2005hk ($M_{J}\sim-17.7$ mag, $M_{H}\sim$$-18.0$ mag, and $\Delta$ $m_{15}$($B$)$\sim1.6$ mag), yet differ significantly from SNe 2008ha and 2010ae ($M_{J}\sim-14 - -15$ mag and $\Delta$ $m_{15}$($B$)$\sim2.4-2.7$ mag). The estimated rise time is $12.0 \pm 3.0$ days, which is significantly shorter than that of SN 2005hk or any other Ia SNe. The rapid rise indicates that the $^{56}$Ni distribution may extend into the outer layer or that the effective opacity may be lower than that in normal SNe Ia. The late-phase spectrum exhibits broader emission lines than those of SN 2005hk by a factor of 6--8. Such high velocities of the emission lines indicate that the density profile of the inner ejecta extends more than that of SN 2005hk. We argue that the most favored explosion scenario is a `failed deflagration' model, although the pulsational delayed detonations is not excluded.

A rapidly rotating neutron star with strong magnetic fields, called magnetar, is a possible candidate for the central engine of long gamma-ray bursts and hypernovae (HNe). We solve the evolution of a shock wave driven by the wind from magnetar and evaluate the temperature evolution, by which we estimate the amount of $^{56}$Ni that produces a bright emission of HNe. We obtain a constraint on the magnetar parameters, namely the poloidal magnetic field strength ($B_p$) and initial angular velocity ($\Omega_i$), for synthesizing enough $^{56}$Ni mass to explain HNe ($M_{^{56}\mathrm{Ni}}\gtrsim 0.2M_\odot$), i.e. $(B_p/10^{16}~\mathrm{G})^{1/2}(\Omega_i/10^4~\mathrm{rad~s}^{-1})\gtrsim 0.7$.

We investigate the fraction of metal nuclei in the relativistic jets of gamma-ray bursts associated with core-collapse supernovae. We simulate the fallback in jet-induced explosions with two-dimensional relativistic hydrodynamics calculations and the jet acceleration with steady, radial, relativistic magnetohydrodynamics calculations, and derive detail nuclear composition of the jet by postprocessing calculation. We found that if the temperature at the jet launch site is above $4.7\times 10^9$K, quasi-statistical equilibrium (QSE) is established and heavy nuclei are dissociated to light particles such as $^4$He during the acceleration of the jets. The criterion for the survival of metal nuclei is written in terms of the isotropic jet luminosity as $L_{\rm j}^{\rm iso} \lesssim 3.9\times 10^{50}(R_{\rm i}/10^7{\rm cm})^2 (1+\sigma_{\rm i})~{\rm erg~s^{-1}}$, where $R_{\rm i}$ and $\sigma_{\rm i}$ are the initial radius of the jets and the initial magnetization parameter, respectively. If the jet is initially dominated by radiation field (i.e., $\sigma_{\rm i} \ll 1$) and the isotropic luminosity is relatively high ($L_{\rm j}^{\rm iso}\gtrsim 4\times 10^{52}~{\rm erg s^{-1}}$), the metal nuclei cannot survive in the jet. On the other hand, if the jet is mainly accelerated by magnetic field (i.e., $\sigma_{\rm i} \gg1$), metal nuclei initially contained in the jet can survive without serious dissociation even for the case of high luminosity jet. If the jet contains metal nuclei, the dominant nuclei are $^{28}$Si, $^{16}$O, and $^{32}$S and the mean mass number can be $\langle$A$\rangle\sim25$.

The Kiso Supernova Survey (KISS) is a high-cadence optical wide-field supernova (SN) survey. The primary goal of the survey is to catch the very early light of a SN, during the shock breakout phase. Detection of SN shock breakouts combined with multi-band photometry obtained with other facilities would provide detailed physical information on the progenitor stars of SNe. The survey is performed using a 2.2x2.2 deg field-of-view instrument on the 1.05-m Kiso Schmidt telescope, the Kiso Wide Field Camera (KWFC). We take a three-minute exposure in g-band once every hour in our survey, reaching magnitude g~20-21. About 100 nights of telescope time per year have been spent on the survey since April 2012. The number of the shock breakout detections is estimated to be of order of 1 during our 3-year project. This paper summarizes the KISS project including the KWFC observing setup, the survey strategy, the data reduction system, and CBET-reported SNe discovered so far by KISS.

We present our discovery of dramatic variability in SDSS J1100+4421 by the high-cadence transient survey Kiso Supernova Survey (KISS). The source brightened in the optical by at least a factor of three within about half a day. Spectroscopic observations suggest that this object is likely a narrow-line Seyfert 1 galaxy (NLS1) at z=0.840, however with unusually strong narrow emission lines. The estimated black hole mass of ~ 10^7 Msun implies bolometric nuclear luminosity close to the Eddington limit. SDSS J1100+4421 is also extremely radio-loud, with a radio loudness parameter of R ~ 4 x 10^2 - 3 x 10^3, which implies the presence of relativistic jets. Rapid and large-amplitude optical variability of the target, reminiscent of that found in a few radio- and gamma-ray loud NLS1s, is therefore produced most likely in a blazar-like core. The 1.4 GHz radio image of the source shows an extended structure with a linear size of about 100 kpc. If SDSS J1100+4421 is a genuine NLS1, as suggested here, this radio structure would then be the largest ever discovered in this type of active galaxies.

The most iron-deficient stars in the Milky Way provide important observational constraints on the nature of astrophysical objects that have enriched the primordial gas with heavy elements from which these stars were formed. Among them, the recently discovered iron-deficient star SMSS J031300.36-670839.3 shows a remarkable chemical composition with non-detection of iron ([Fe/H]$<-7.1$) and large enhancement of carbon and magnesium relative to calcium. We investigate the supernova yields of metal-free (Population III) stars to interpret the observed abundance pattern for this star. We report that the high [C/Ca] and [C/Mg] ratios and upper limits determined for other elemental abundances are well reproduced with the yields of core-collapse supernovae (that have normal kinetic energies of explosion $E$ of $E_{51}=E/10^{51}$erg$=1$) or hypernovae ($E_{51}\geq 10$) of the Population III 25$M_{\odot}$ or 40$M_{\odot}$ stars. The best-fit model assumes that the explosion of the Population III progenitor undergoes extensive matter mixing and fallback, leaving behind a black hole remnant. In these models, Ca is produced by static/explosive O burning and incomplete Si burning in the Population III supernova/hypernova, in contrast to the suggestion that Ca is originated from the hot CNO cycle during the presupernova evolution. Among the five most iron-poor stars with [Fe/H]$<-4.5$, four carbon-enhanced stars are consistent with the faint supernova models with the ejected mass of $^{56}$Ni less than 10$^{-3}M_{\odot}$.

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 perform cosmological hydrodynamics simulations with non-equilibrium primordial chemistry to obtain 59 minihalos that host first stars. The obtained minihalos are used as initial conditions of local three dimensional radiation hydrodynamics simulations to investigate the formation of the first stars. We find two-thirds of the minihalos host multiple stars, while the rest of them have single stars. The mass of the stars found in our simulations are in the range of 1 Msun \la M \la 300 Msun, peaking at several x 10 Msun. Most of the very massive stars of \ga 140 Msun are born as single stars, although not all of the single stars are very massive. We also find a few stars of \la 1 Msun that are kicked by the gravitational three body interactions to the position distant from the center of mass. The frequency that a star forming minihalo contains a binary system is \sim 50%. We also investigate the abundance pattern of the stellar remnants by summing up the contributions from the first stars in the simulations. Consequently, the pattern is compatible with that of the low metallicity Damped Lyman-alpha systems or the Extremely Metal Poor (EMP) stars, if the mass spectrum obtained in our experiment is shifted to the low mass side by 0.2 dex. If we consider the case that an EMP star is born in the remnant of the individual minihalo without mixing with others, the chemical signature of the pair instability supernova is more prominent, because most of them are born as single stars.

We investigate the random walk process in relativistic flow. In the relativistic flow, photon propagation is concentrated in the directions of the flow velocity due to relativistic beaming effect. We show that, in the pure scattering case, the number of scatterings is proportional to the size parameter $\xi\equiv L/l_0$ if the flow velocity $\beta\equiv v/c$ satisfies $\beta/\Gamma\gg \xi^{-1}$, while it is proportional to $\xi^2$ if $\beta/\Gamma\ll \xi^{-1}$ where $L$ and $l_0$ are the size of the system in the observer frame and the mean free path in the comoving frame, respectively. We also examine the photon propagation in the scattering and absorptive medium. We find that, if the optical depth for absorption $\tau_{\rm a}$ is considerably smaller than the optical depth for scattering $\tau_{\rm s}$ ($\tau_{\rm a}/\tau_{\rm s} \ll 1$) and the flow velocity satisfies $\beta\gg \sqrt{2\tau_{\rm a}/\tau_{\rm s}}$, the effective optical depth is approximated by $\tau_*\simeq\tau_{\rm a}(1+\beta)/\beta$. Furthermore, we perform Monte Carlo simulations of radiative transfer and compare the results with the analytic expression for the number of scattering. The analytic expression is consistent with the results of the numerical simulations. The expression derived in this Letter can be used to estimate the photon production site in relativistic phenomena, e.g., gamma-ray burst and active galactic nuclei.

We show that the low ratios of $\alpha$ elements (Mg, Si, and Ca) to Fe recently found for a small fraction of extremely metal-poor stars can be naturally explained with the nucleosynthesis yields of core-collapse supernovae, i.e., $13-25M_\odot$ supernovae, or hypernovae. For the case without carbon enhancement, the ejected iron mass is normal, consistent with observed light curves and spectra of nearby supernovae. On the other hand, the carbon enhancement requires much smaller iron production, and the low [$\alpha$/Fe] of carbon enhanced metal-poor stars can also be reproduced with $13-25M_\odot$ faint supernovae or faint hypernovae. Iron-peak element abundances, in particular Zn abundances, are important to put further constraints on the enrichment sources from galactic archaeology surveys.

After the Big Bang nucleosynthesis, the first heavy element enrichment in the Universe was made by a supernova (SN) explosion of a population (Pop) III star (Pop III SN). The abundance ratios of elements produced from Pop III SNe are recorded in abundance patterns of extremely metal-poor (EMP) stars. The observations of the increasing number of EMP stars have made it possible to statistically constrain the explosion properties of Pop III SNe. We present Pop III SN models whose nucleosynthesis yields well-reproduce individually the abundance patterns of 48 such metal-poor stars as [Fe/H] $\mathrel{\rlap{\lower 4pt \hbox{$\sim$}}\raise 1pt \hbox {$<$}}-3.5$. We then derive relations between the abundance ratios of EMP stars and certain explosion properties of Pop III SNe: the higher [(C+N)/Fe] and [(C+N)/Mg] ratios correspond to the smaller ejected Fe mass and the larger compact remnant mass, respectively. Using these relations, the distributions of the abundance ratios of EMP stars are converted to those of the explosion properties of Pop III SNe. Such distributions are compared with those of the explosion properties of present day SNe: The distribution of the ejected Fe mass of Pop III SNe has the same peak as that of the resent day SNe but shows an extended tail down to $\sim10^{-2}-10^{-5}M_\odot$, and the distribution of the mass of the compact remnant of Pop III SNe is as wide as that of the present day stellar-mass black holes. Our results demonstrate the importance of large samples of EMP stars obtained by ongoing and future EMP star surveys and subsequent high-dispersion spectroscopic observations in clarifying the nature of Pop III SNe in the early Universe.

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

We present detailed chemical abundances for the bright carbon-enhanced metal-poor (CEMP) star BD+44 493, previously reported on by Ito et al. Our measurements confirm that BD+44 493 is an extremely metal-poor ([Fe/H]=-3.8) subgiant star with excesses of carbon and oxygen. No significant excesses are found for nitrogen and neutron-capture elements (the latter of which place it in the CEMP-no class of stars). Other elements that we measure exhibit abundance patterns that are typical for non-CEMP extremely metal-poor stars. No evidence for variations of radial velocity have been found for this star. These results strongly suggest that the carbon enhancement in BD+44 493 is unlikely to have been produced by a companion asymptotic giant-branch star and transferred to the presently observed star, nor by pollution of its natal molecular cloud by rapidly-rotating, massive, mega metal-poor ([Fe/H] < -6.0) stars. A more likely possibility is that this star formed from gas polluted by the elements produced in a "faint" supernova, which underwent mixing and fallback, and only ejected small amounts of elements of metals beyond the lighter elements. The Li abundance of BD+44 493 (A(Li)=log(Li/H)+12=1.0) is lower than the Spite plateau value, as found in other metal-poor subgiants. The upper limit on Be abundance (A(Be)=log(Be/H)+12<-1.8) is as low as those found for stars with similarly extremely-low metallicity, indicating that the progenitors of carbon- (and oxygen-) enhanced stars are not significant sources of Be, or that Be is depleted in metal-poor subgiants with effective temperatures of ~5400K.

We propose an alternative explanation for the post-breakout emission of SN 2008D associated with the X-ray transient 080109. Observations of this object show a very small contrast of 0.35 dex between the light-curve minimum occurring soon after the breakout, and the main luminosity peak that is due to radioactive heating of the ejecta. Hydrodynamical models show that the cooling of a shocked Wolf-Rayet star leads to a much greater difference (> 0.9 dex). Our proposed scenario is that of a jet produced during the explosion which deposits 56Ni-rich material in the outer layers of the ejecta. The presence of high-velocity radioactive material allows us to reproduce the complete luminosity evolution of the object. Without outer 56Ni it could be possible to reproduce the early emission purely from cooling of the shocked envelope by assuming a larger progenitor than a Wolf-Rayet star, but that would require an initial density structure significantly different from what is predicted by stellar evolution models. Analytic models of the cooling phase have been proposed reproduce the early emission of SN 2008D with an extended progenitor. However, we found that the models are valid only until 1.5 days after the explosion where only two data of SN 2008D are available. We also discuss the possibility of the interaction of the ejecta with a binary companion, based on published analytic expressions. However, the binary separation required to fit the early emission should be < 3 Rsun which is too small for a system containing two massive stars.

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

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

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

Shock breakout is the brightest radiative phenomenon in a supernova (SN) but is difficult to be observed owing to the short duration and X-ray/ultraviolet (UV)-peaked spectra. After the first observation from the rising phase reported in 2008, its observability at high redshift is attracting enormous attention. We perform multigroup radiation hydrodynamics calculations of explosions for evolutionary presupernova models with various main-sequence masses $M_{\rm MS}$, metallicities $Z$, and explosion energies $E$. We present multicolor light curves of shock breakout in Type II plateau SNe, being the most frequent core-collapse SNe, and predict apparent multicolor light curves of shock breakout at various redshifts $z$. We derive the observable SN rate and reachable redshift as functions of filter $x$ and limiting magnitude $m_{x,{\rm lim}}$ by taking into account an initial mass function, cosmic star formation history, intergalactic absorption, and host galaxy extinction. We propose a realistic survey strategy optimized for shock breakout. For example, the $g'$-band observable SN rate for $m_{g',{\rm lim}}=27.5$ mag is 3.3 SNe degree$^{-2}$ day$^{-1}$ and a half of them locates at $z\geq1.2$. It is clear that the shock breakout is a beneficial clue to probe high-$z$ core-collapse SNe. We also establish ways to identify shock breakout and constrain SN properties from the observations of shock breakout, brightness, time scale, and color. We emphasize that the multicolor observations in blue optical bands with $\sim$ hour intervals, preferably over $\geq2$ continuous nights, are essential to efficiently detect, identify, and interpret shock breakout.

We show that the recently observed elemental abundance pattern of the carbon-rich metal-poor Damped Lyman $\alpha$ (DLA) system is in excellent agreement with the nucleosynthesis yields of faint core-collapse supernovae of primordial stars. The observed abundance pattern is not consistent with the nucleosynthesis yields of pair-instability supernovae. The DLA abundance pattern is very similar to that of carbon-rich extremely metal-poor (EMP) stars, and the contributions from low-mass stars and/or binary effects should be very small in DLAs. This suggests that chemical enrichment by the first stars in the first galaxies is driven by core-collapse supernovae from $\sim 20-50 M_\odot$ stars, and also supports the supernova scenario as the enrichment source of EMP stars in the Milky Way Galaxy.

We observed the Cygnus Loop with XMM-Newton (9 pointings) and Suzaku (32 pointings) between 2002 and 2008. The total effective exposure time is 670.2 ks. By using all of the available data, we intended to improve a signal-to-noise ratio of the spectrum. Accordingly, the accumulated spectra obtained by the XIS and the EPIC show some line features around 3 keV that are attributed to the S He$\beta$ and Ar He$\alpha$ lines, respectively. Since the Cygnus Loop is an evolved ($\sim$10,000 yr) supernova remnant whose temperature is relatively low ($<$1 keV) compared with other young remnants, its spectrum is generally faint above 3.0 keV, no emission lines, such as the Ar-K line have ever been detected. The detection of the Ar-K line is the first time and we found that its abundance is significantly higher than that of the solar value; 9.0$^{+4.0}_{-3.8}$ and 8.4$^{+2.5}_{-2.7}$ (in units of solar), estimated from the XIS and the EPIC spectra, respectively. We conclude that the Ar-K line originated from the ejecta of the Cygnus Loop. Follow-up X-ray observations to tightly constrain the abundances of Ar-rich ejecta will be useful to accurately estimate the progenitor's mass.

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.

We investigate the formation of dust grains in the ejecta of a SN IIb and their evolution in the shocked gas in the SNR by considering the uniform and power-law density structures for the CSM. Based on these calculations, we also simulate the time evolution of thermal emission from the shock-heated dust in the SNR and compare the results with the observations of Cas A SNR. We find that the total mass of dust formed in the SN IIb is as large as 0.167 M_sun but the average radius of dust is smaller than 0.01 mum and is significantly different from those in SNe II-P with the massive H-envelope. In the explosion with the small-mass H-envelope, the expanding He core undergoes little deceleration, so that the gas density in the He core is too low for large-sized grains to form. In addition, the low-mass H-envelope of the SN IIb leads to the early arrival of the reverse shock at the dust-forming region. If the CSM is more or less spherical, therefore, the newly formed grains would be completely destroyed in the relatively dense shocked gas for the CSM density of n_H > 0.1 cm^{-3}. However, the actual CSM is likely to be non-spherical, so that a part of grains could be ejected into the ISM without being shocked. We demonstrate that the time evolution of the SED by thermal dust emission is sensitive to the ambient gas density and structure that affects the passage of the reverse shock into the ejecta. Thus, the SED evolution well reflects the evolution of dust through erosion by sputtering and stochastic heating. For Cas A, we consider the CSM produced by the steady mass loss of ~8x10^{-5} M_sun/yr during the supergiant phase. Then we find the infrared SED of Cas A is reasonably reproduced by thermal emission from the newly formed dust of 0.08 M_sun, which consists of 0.008 M_sun shocked warm dust and 0.072 M_sun unshocked cold dust.

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

We present optical spectroscopic and photometric observations of supernova (SN) 2008D, associated with the luminous X-ray transient 080109, at >300 days after the explosion (nebular phases). We also give flux measurements of emission lines from the H II region at the site of the SN, and estimates of the local metallicity. The brightness of the SN at nebular phases is consistent with the prediction of the explosion models with an ejected 56Ni mass of 0.07 Msun, which explains the light curve at early phases. The [O I] line in the nebular spectrum shows a double-peaked profile while the [Ca II] line does not. The double-peaked [O I] profile strongly indicates that SN 2008D is an aspherical explosion. The profile can be explained by a torus-like distribution of oxygen viewed from near the plane of the torus. We suggest that SN 2008D is a side-viewed, bipolar explosion with a viewing angle of > 50^{\circ} from the polar direction.

We present an optical spectropolarimetric observation of the unique Type Ib supernova (SN) 2005bf at 8 days after the second maximum. The data, combined with the polarization spectrum taken at 6 days before the second maximum (Maund et al. 2007a), enable us to closely examine the intrinsic properties of the SN. The polarization percentage is smaller at the later epoch over a wide wavelength range, while the position angle is similar at the two epochs. We find that an assumption of complete depolarization of strong lines at the emission peak is not necessarily correct. The intrinsic polarization of the SN is larger, and thus, the ejecta of SN 2005bf would be more asymmetric than previously expected. The axis ratio of the photosphere projected on the sky deviates from unity by at least 20 %. If the position angle of interstellar polarization is aligned with the spiral arm of the host galaxy, the deviation is larger than 25 %. The line polarization at the He I, Ca II and Fe II lines is also smaller at the later epoch. The changes in the position angle across these lines, which were observed at the earlier epoch, are still marginally present at the later epoch. The similar polarimetric behavior suggests that the distributions of these ions are correlated. Properties of polarization, as well as the light curve and the spectra both at photospheric and nebular phases, can be explained by an aspherical, possibly unipolar explosion of a WN star in which the blob of 56Ni penetrates C+O core and stops within the He layer.

The very peculiar abundance patterns observed in extremely metal-poor (EMP) stars can not be explained by ordinary supernova nucleosynthesis but can be well-reproduced by nucleosynthesis in hyper-energetic and hyper-aspherical explosions, i.e., Hypernovae (HNe). Previously, such HNe have been observed only as Type Ic supernovae. Here, we examine the properties of recent Type Ib supernovae (SNe Ib). In particular, SN Ib 2008D associated with the luminous X-ray transient 080109 is found to be a more energetic explosion than normal core-collapse supernovae. We estimate that the progenitor's main sequence mass is 20--25 M_sun and a kinetic energy of explosion is ~ 6 x 10^{51} erg. These properties are intermediate between those of normal SNe and hypernovae associated with gamma-ray bursts. Such energetic SNe Ib can make important contribution to the chemical enrichment in the early Universe.

We review the final stages of stellar evolution, supernova properties, and chemical yields as a function of the progenitor's mass M. (1) 8 - 10 Ms stars are super-AGB stars when the O+Ne+Mg core collapses due to electron capture. These AGB-supernovae may constitute an SN 2008S-like sub-class of Type IIn supernovae. These stars produce little alpha-elements and Fe-peak elements, but are important sources of Zn and light p-nuclei. (2) 10 - 90 Ms stars undergo Fe-core collapse. Nucleosynthesis in aspherical explosions is important, as it can well reproduce the abundance patterns observed in extremely metal-poor stars. (3) 90 - 140 Ms stars undergo pulsational nuclear instabilities at various nuclear burning stages, including O and Si-burning. (4) 140 - 300 Ms stars become pair-instability supernovae, if the mass loss is small enough. (5) Stars more massive than 300 Ms undergo core-collapse to form intermediate mass black holes.

We present recent advances in theoretical studies of the formation and evolution of dust in primordial supernovae (SNe) that are considered to be the main sources of dust in the early universe. Being combined with the results of calculations of dust formation in the ejecta of Population III SNe, the investigations of the evolution of newly formed dust within supernova remnants (SNRs) show that smaller grains are predominantly destroyed by sputtering in the shocked gas, while larger grains are injected into the ambient medium. The mass of dust grains surviving the destruction in SNRs reaches up to 0.1--15 $M_\odot$, which is high enough to account for the content of dust observed for the host galaxies of quasars at $z > 5$. In addition, the transport of dust formed in the ejecta causes the formation of low-mass stars in the dense shell of primordial SNRs and affects the elemental composition of those stars. We also show that the flat extinction curve is expected in the high-redshift universe where SNe are the possible sources of dust.

Core--collapsed supernovae (CCSNe) have been considered to be one of sources of dust in the universe. What kind and how much mass of dust are formed in the ejecta and are injected into the interstellar medium (ISM) depend on the type of CCSNe, through the difference in the thickness (mass) of outer envelope. In this review, after summarizing the existing results of observations on dust formation in CCSNe, we investigate formation of dust in the ejecta and its evolution in the supernova remnants (SNRs) of Type II--P and Type IIb SNe. Then, the time evolution of thermal emission from dust in the SNR of Type IIb SN is demonstrated and compared with the observation of Cas A. We find that the total dust mass formed in the ejecta does not so much depend on the type; $\sim 0.3-0.7 M_{\odot}$ in Type II--P SNe and $\sim 0.13 M_{\odot}$ in Type IIb SN. However the size of dust sensitively depends on the type, being affected by the difference in the gas density in the ejecta: the dust mass is dominated by grains with radii larger than 0.03 $\mu$m in Type II-P, and less than 0.006 $\mu$m in Type IIb, which decides the fate of dust in the SNR. The surviving dust mass is $\sim 0.04-0.2 M_{\odot}$ in the SNRs of Type II--P SNe for the ambient hydrogen density of $n_{\rm H}=10.0-1.0$ cm$^{-3}$, while almost all dust grains are destroyed in the SNR of Type IIb. The spectral energy distribution (SED) of thermal emission from dust in SNR well reflects the evolution of dust grains in SNR through erosion by sputtering and stochastic heating. The observed SED of Cas A SNR is reasonably reproduced by the model of dust formation and evolution for Type IIb SN.

We present a theoretical model for supernova (SN) 2008D associated with the luminous X-ray transient 080109. The bolometric light curve and optical spectra of the SN are modelled based on the progenitor models and the explosion models obtained from hydrodynamic/nucleosynthetic calculations. We find that SN 2008D is a more energetic explosion than normal core-collapse supernovae, with an ejecta mass of Mej = 5.3 +- 1.0 Msun and a kinetic energy of E = 6.0 +- 2.5 x 10^{51} erg. The progenitor star of the SN has a 6-8 Msun He core with essentially no H envelope (< 5 x 10^{-4} Msun) prior to the explosion. The main-sequence mass of the progenitor is estimated to be Mms =20-25 Msun, with additional systematic uncertainties due to convection, mass loss, rotation, and binary effects. These properties are intermediate between those of normal SNe and hypernovae associated with gamma-ray bursts. The mass of the central remnant is estimated as 1.6 - 1.8 Msun, which is near the boundary between neutron star and black hole formation.

We performed optical spectroscopy and photometry of SN 2006gy at late time, ~400 days after the explosion, with the Subaru/FOCAS in a good seeing condition. We found that the SN faded by ~3 mag from ~200 to ~400 days after the explosion (i.e., by ~5 mag from peak to ~400 days) in R band. The overall light curve is marginally consistent with the 56Ni heating model, although the flattening around 200 days suggests the optical flux declined more steeply between ~200 and ~400 days. The late time spectrum was quite peculiar among all types of SNe. It showed many intermediate width (~2000 km/s FWHM) emission lines, e.g., [Fe II], [Ca II], and Ca II. The absence of the broad [O I] 6300, 6364 line and weakness of [Fe II] and [Ca II] lines compared with Ca II IR triplet would be explained by a moderately high electron density in the line emitting region. This high density assumption seems to be consistent with the large amount of ejecta and low expansion velocity of SN 2006gy. The H-alpha line luminosity was as small as ~1x10^39 erg/s, being comparable with those of normal Type II SNe at similar epochs. Our observation indicates that the strong CSM interaction had almost finished by ~400 days. If the late time optical flux is purely powered by radioactive decay, at least M_Ni ~ 3 M_sun should be produced at the SN explosion. In the late phase spectrum, there were several unusual emission lines at 7400--8800 AA and some of them might be due to Ti or Ni synthesized at the explosion. (abridged)

Theoretically expected natures of a supernova driven by a wind/jet are discussed. Approximate analytical formulations are derived to clarify basic physical processes involved in the wind/jet-driven explosions, and it is shown that the explosion properties are characterized by the energy injection rate (Edot_iso) and the mass injection rate (Mdot_iso). To explain observations of SN 1998bw associated with Gamma-Ray Burst (GRB) 980425, the following conditions are required: Edot_iso Mdot_iso > ~ 10^{51} erg M_sun s^{-2} and Edot_iso > ~ 2 x 10^{52} erg s^{-1} (if the wind Lorentz factor Gamma_w ~ 1) or Edot_iso > ~ 7 x 10^{52} erg s^{-1} (if Gamma_w >> 1). In SN 1998bw, 56Ni (~ 0.4M_sun) is probably produced in the shocked stellar mantle, not in the wind. The expected natures of SNe, e.g., ejected 56Ni masses and ejecta masses, vary depending on Edot_iso and Mdot_iso. The sequence of the SN properties from high Edot_iso and Mdot_iso to low Edot_iso and Mdot_iso is the following: SN 1998bw-like -- intermediate case -- low mass ejecta (< ~ 1M_sun$) where 56Ni is from the wind -- whole collapse. This diversity may explain the diversity of supernovae associated with GRBs. Our result can be used to constrain natures of the wind/jet, which is linked to the central engine of GRBs, by studying properties of the associated supernovae.

We present optical UBVRI photometric and spectroscopic data of the type Ibn supernova SN 2006jc, until the onset of the dust forming phase. The optical spectrum shows a blue continuum and is dominated by the presence of moderately narrow (velocity ~2500 km/s) He I emission lines superimposed over a relatively weak supernova spectrum. The helium lines are produced in a pre-existing He rich circumstellar shell. The observed helium line fluxes indicate the circumstellar shell is dense, with a density of ~10^9 - 10^{10} cm^{-3}. The helium mass in this shell is estimated to be <~0.07 Msun. The optical light curves show a clear signature of dust formation, indicated by a sharp decrease in the magnitudes around day 50, accompanied by a reddening of the colours. The evolution of the optical light curves during the early phase and that of the uvoir bolometric light curve at all phases is reasonably similar to normal Ib/c supernovae.

There have been attempts to fit the abundance patterns of extremely metal-poor stars with supernova nucleosynthesis models for the lighter elements than Zn. On the other hand, observations have revealed that the presence of EMP stars with peculiarly high ratio of "weak r-process elements" Sr, Y and Zr. Although several possible processes were suggested for the origin of these elements, the complete solution for reproducing those ratios is not found yet. In order to reproduce the abundance patterns of such stars, we investigate a model with neutron rich matter ejection from the inner region of the conventional mass-cut. We find that explosive nucleosynthesis in a high energy supernova (or "hypernova") can reproduce the high abundances of Sr, Y and Zr but that the enhancements of Sr, Y and Zr are not achieved by nucleosynthesis in a normal supernova. Our results imply that, if these elements are ejected from a normal supernova, nucleosynthesis in higher entropy flow than that of the supernova shock is required.

We investigate hydrodynamical and nucleosynthetic properties of the jet-induced explosion of a population III $40M_\odot$ star and compare the abundance patterns of the yields with those of the metal-poor stars. We conclude that (1) the ejection of Fe-peak products and the fallback of unprocessed materials can account for the abundance patterns of the extremely metal-poor (EMP) stars and that (2) the jet-induced explosion with different energy deposition rates can explain the diversity of the abundance patterns of the metal-poor stars. Furthermore, the abundance distribution after the explosion and the angular dependence of the yield are shown for the models with high and low energy deposition rates $\dot{E}_{\rm dep}=120\times10^{51} {\rm ergs s^{-1}}$ and $1.5\times10^{51} {\rm ergs s^{-1}}$. We also find that the peculiar abundance pattern of a Si-deficient metal-poor star HE 1424--0241 can be reproduced by the angle-delimited yield for $\theta=30^\circ-35^\circ$ of the model with $\dot{E}_{\rm dep}=120\times10^{51} {\rm ergs s^{-1}}$.

We investigate the evolution of dust formed in Population III supernovae (SNe) by considering its transport and processing by sputtering within the SN remnants (SNRs). We find that the fates of dust grains within SNRs heavily depend on their initial radii $a_{\rm ini}$. For Type II SNRs expanding into the ambient medium with density of $n_{\rm H,0} = 1$ cm$^{-3}$, grains of $a_{\rm ini} < 0.05$ $\mu$m are detained in the shocked hot gas and are completely destroyed, while grains of $a_{\rm ini} > 0.2$ $\mu$m are injected into the surrounding medium without being destroyed significantly. Grains with $a_{\rm ini}$ = 0.05-0.2 $\mu$m are finally trapped in the dense shell behind the forward shock. We show that the grains piled up in the dense shell enrich the gas up to 10$^{-6}-10^{-4}$ $Z_\odot$, high enough to form low-mass stars with 0.1-1 $M_\odot$. In addition, [Fe/H] in the dense shell ranges from -6 to -4.5, which is in good agreement with the ultra-metal-poor stars with [Fe/H] < -4. We suggest that newly formed dust in a Population III SN can have great impacts on the stellar mass and elemental composition of Population II.5 stars formed in the shell of the SNR.

The first metal enrichment in the universe was made by supernova (SN) explosions of population (Pop) III stars. The trace remains in abundance patterns of extremely metal-poor (EMP) stars. We investigate the properties of nucleosynthesis in Pop III SNe by means of comparing their yields with the abundance patterns of the EMP stars. We focus on (1) jet-induced SNe with various energy deposition rates [$\dot{E}_{\rm dep}=(0.3-1500)\times10^{51}{\rm ergs s^{-1}}$], and (2) SNe of stars with various main-sequence masses ($M_{\rm ms}=13-50M_\odot$) and explosion energies [$E=(1-40)\times10^{51}$ergs]. The varieties of Pop III SNe can explain varieties of the EMP stars: (1) higher [C/Fe] for lower [Fe/H] and (2) trends of abundance ratios [X/Fe] against [Fe/H].

Jet-induced supernovae (SNe) have been suggested to occur in gamma-ray bursts (GRBs) and highly-energetic SNe (hypernovae). I investigate hydrodynamical and nucleosynthetic properties of the jet-induced explosion of a population III $40_\odot$ star with a two-dimensional special relativistic hydrodynamical code. The abundance distribution after the explosion and the angular dependence of the yield are obtained for the models with high and low energy deposition rates $\dot{E}_{\rm dep}=120\times10^{51} {\rm ergs s^{-1}}$ and $1.5\times10^{51} {\rm ergs s^{-1}}$. The ejection of Fe-peak products and the fallback of unprocessed materials in the jet-induced SNe account for the abundance patterns of the extremely metal-poor (EMP) stars. It is also found that the peculiar abundance pattern of a Si-deficient metal-poor star HE 1424--0241 is reproduced by the angle-delimited yield for $\theta=30^\circ-35^\circ$ of the model with $\dot{E}_{\rm dep}=120\times10^{51} {\rm ergs s^{-1}}$. Furthermore, I compare the yield of the jet-induced explosion with that of the spherical explosion and confirm the ejection and fallback in the jet-induced explosion is almost equivalent to the "mixing-fallback" in spherical explosions. In contrast to the spherical models, however, the high-entropy environment is realized in the jet-induced explosion and thus [(Sc, Ti, V, Cr, Co, Zn)/Fe] are enhanced. The enhancements of [Sc/Fe] and [Ti/Fe] improve agreements with the abundance patterns of the EMP stars.

Several short-lived radionuclides (SLRs) were present in the early solar system, some of which should have formed just prior to or soon after the solar system formation. Stellar nucleosynthesis has been proposed as the mechanism for production of SLRs in the solar system, but no appropriate stellar source has been found to explain the abundances of all solar system SLRs. In this study, we propose a faint supernova with mixing and fallback as a stellar source of SLRs with mean lives of <5 Myr (26Al, 41Ca, 53Mn, and 60Fe) in the solar system. In such a supernova, the inner region of the exploding star experiences mixing, a small fraction of mixed materials is ejected, and the rest undergoes fallback onto the core. The modeled SLR abundances agree well with their solar system abundances if mixing-fallback occurs within the C/O-burning layer. In some cases, the initial solar system abundances of the SLRs can be reproduced within a factor of 2. The dilution factor of supernova ejecta to the solar system materials is ~10E-4 and the time interval between the supernova explosion and the formation of oldest solid materials in the solar system is ~1 Myr. If the dilution occurred due to spherically symmetric expansion, a faint supernova should have occurred nearby the solar system forming region in a star cluster.

We present a theoretical model for Type Ib supernova (SN) 2006jc. We calculate the evolution of the progenitor star, hydrodynamics and nucleosynthesis of the SN explosion, and the SN bolometric light curve (LC). The synthetic bolometric LC is compared with the observed bolometric LC constructed by integrating the UV, optical, near-infrared (NIR), and mid-infrared (MIR) fluxes. The progenitor is assumed to be as massive as $40M_\odot$ on the zero-age main-sequence. The star undergoes extensive mass loss to reduce its mass down to as small as $6.9M_\odot$, thus becoming a WCO Wolf-Rayet star. The WCO star model has a thick carbon-rich layer, in which amorphous carbon grains can be formed. This could explain the NIR brightening and the dust feature seen in the MIR spectrum. We suggest that the progenitor of SN 2006jc is a WCO Wolf-Rayet star having undergone strong mass loss and such massive stars are the important sites of dust formation. We derive the parameters of the explosion model in order to reproduce the bolometric LC of SN 2006jc by the radioactive decays: the ejecta mass $4.9M_\odot$, hypernova-like explosion energy $10^{52}$ ergs, and ejected $^{56}$Ni mass $0.22M_\odot$. We also calculate the circumstellar interaction and find that a CSM with a flat density structure is required to reproduce the X-ray LC of SN 2006jc. This suggests a drastic change of the mass-loss rate and/or the wind velocity that is consistent with the past luminous blue variable (LBV)-like event.

The only supernovae (SNe) to have shown early gamma-ray or X-ray emission thus far are overenergetic, broad-lined Type Ic SNe (Hypernovae - HNe). Recently, SN 2008D shows several novel features: (i) weak XRF, (ii) an early, narrow optical peak, (iii) disappearance of the broad lines typical of SNIc HNe, (iv) development of He lines as in SNeIb. Detailed analysis shows that SN 2008D was not a normal SN: its explosion energy (KE ~ 6*10^{51} erg) and ejected mass (~7 Msun) are intermediate between normal SNeIbc and HNe. We derive that SN 2008D was originally a ~30Msun star. When it collapsed a black hole formed and a weak, mildly relativistic jet was produced, which caused the XRF. SN 2008D is probably among the weakest explosions that produce relativistic jets. Inner engine activity appears to be present whenever massive stars collapse to black holes.

SN 2006jc is a peculiar supernova (SN), in which the formation of dust has been confirmed at an early epoch of ~50 days after the explosion. We investigate the possibility of such an earlier formation of dust grains in the expanding ejecta of SN 2006jc, applying the Type Ib SN model that is developed to reproduce the observed light curve. We find that the rapid decrease of the gas temperature in SN 2006jc enables the condensation of C grains in the C-rich layer at 40-60 days after the explosion, which is followed by the condensation of silicate and oxide grains until ~200 days. The average radius of each grain species is confined to be less than 0.01 micron due to the low gas density at the condensation time. The calculated total dust mass reaches ~1.5 Msun, of which C dust shares 0.7 Msun. On the other hand, based on the calculated dust temperature, we show that the dust species and mass evaluated to reproduce the spectral energy distribution observed by AKARI and MAGNUM at day 200 are different from those obtained by the dust formation calculations; the dust species contributing to the observed flux are hot C and FeS grains with masses of $5.6 \times 10^{-4}$ Msun and $2.0 \times 10^{-3}$ Msun, respectively, though we cannot defy the presence of a large amount of cold dust such as silicate and oxide grains up to 0.5 Msun. One of the physical processes responsible for the difference between calculated and evaluated masses of C and FeS grains could be considered to be the destruction of small-sized clusters by energetic photons and electrons prevailing within the ejecta at the earlier epoch.

The connection between the long GRBs and Type Ic Supernovae (SNe) has revealed the interesting diversity: (i) GRB-SNe, (ii) Non-GRB Hypernovae (HNe), (iii) X-Ray Flash (XRF)-SNe, and (iv) Non-SN GRBs (or dark HNe). We show that nucleosynthetic properties found in the above diversity are connected to the variation of the abundance patterns of extremely-metal-poor (EMP) stars, such as the excess of C, Co, Zn relative to Fe. We explain such a connection in a unified manner as nucleosynthesis of hyper-aspherical (jet-induced) explosions Pop III core-collapse SNe. We show that (1) the explosions with large energy deposition rate, $\dot{E}_{\rm dep}$, are observed as GRB-HNe and their yields can explain the abundances of normal EMP stars, and (2) the explosions with small $\dot{E}_{\rm dep}$ are observed as GRBs without bright SNe and can be responsible for the formation of the C-rich EMP (CEMP) and the hyper metal-poor (HMP) stars. We thus propose that GRB-HNe and the Non-SN GRBs (dark HNe) belong to a continuous series of BH-forming stellar deaths with the relativistic jets of different $\dot{E}_{\rm dep}$.

$UBVRI$ photometry and medium resolution optical spectroscopy of peculiar Type Ia supernova SN 2005hk are presented and analysed, covering the pre-maximum phase to around 400 days after explosion. The supernova is found to be underluminous compared to "normal" Type Ia supernovae. The photometric and spectroscopic evolution of SN 2005hk is remarkably similar to the peculiar Type Ia event SN 2002cx. The expansion velocity of the supernova ejecta is found to be lower than normal Type Ia events. The spectra obtained $\gsim 200$ days since explosion do not show the presence of forbidden [\ion{Fe}{ii}], [\ion{Fe}{iii}] and [\ion{Co}{iii}] lines, but are dominated by narrow, permitted \ion{Fe}{ii}, NIR \ion{Ca}{ii} and \ion{Na}{i} lines with P-Cygni profiles. Thermonuclear explosion model with Chandrasekhar mass ejecta and a kinetic energy smaller ($\KE = 0.3 \times 10^{51} {\rm ergs}$) than that of canonical Type Ia supernovae is found to well explain the observed bolometric light curve. The mass of \Nifs synthesized in this explosion is $0.18 \Msun$. The early spectra are successfully modeled with this less energetic model with some modifications of the abundance distribution. The late spectrum is explained as a combination of a photospheric component and a nebular component.

Following the recent abundance measurements of Mg, Al, Ca, Fe, and Ni in the black hole X-ray binary XTE J1118+480 using medium-resolution Keck II/ESI spectra of the secondary star (Gonz\'alez Hern\'andez et al. 2006), we perform a detailed abundance analysis including the abundances of Si and Ti. These element abundances, higher than solar, indicate that the black hole in this system formed in a supernova event, whose nucleosynthetic products could pollute the atmosphere of the secondary star, providing clues on the possible formation region of the system, either Galactic halo, thick disk, or thin disk. We explore a grid of explosion models with different He core masses, metallicities, and geometries. Metal-poor models associated with a formation scenario in the Galactic halo provide unacceptable fits to the observed abundances, allowing us to reject a halo origin for this X-ray binary. The thick-disk scenario produces better fits, although they require substantial fallback and very efficient mixing processes between the inner layers of the explosion and the ejecta, making quite unlikely an origin in the thick disk. The best agreement between the model predictions and the observed abundances is obtained for metal-rich progenitor models. In particular, non-spherically symmetric models are able to explain, without strong assumptions of extensive fallback and mixing, the observed abundances. Moreover, asymmetric mass ejection in a supernova explosion could account for the required impulse necessary to launch the system from its formation region in the Galactic thin disk to its current halo orbit.

We present our latest results on near- to mid- infrared observation of SN2006jc at 200 days after the discovery using the Infrared Camera (IRC) on board $AKARI$. The near-infrared (2--5$\mu$m) spectrum of SN2006jc is obtained for the first time and is found to be well interpreted in terms of the thermal emission from amorphous carbon of 800$\pm 10$K with the mass of $6.9\pm 0.5 \times 10^{-5}M_{\odot}$ that was formed in the supernova ejecta. This dust mass newly formed in the ejecta of SN 2006jc is in a range similar to those obtained for other several dust forming core collapse supernovae based on recent observations (i.e., $10^{-3}$--$10^{-5}$$M_{\odot}$). Mid-infrared photometric data with {\it{AKARI}}/IRC MIR-S/S7, S9W, and S11 bands have shown excess emission over the thermal emission by hot amorphous carbon of 800K. This mid-infrared excess emission is likely to be accounted for by the emission from warm amorphous carbon dust of 320$\pm 10$K with the mass of 2.7$^{+0.7}_{-0.5} \times 10^{-3}M_{\odot}$ rather than by the band emission of astronomical silicate and/or silica grains. This warm amorphous carbon dust is expected to have been formed in the mass loss wind associated with the Wolf-Rayet stellar activity before the SN explosion. Our result suggests that a significant amount of dust is condensed in the mass loss wind prior to the SN explosion. A possible contribution of emission bands by precursory SiO molecules in 7.5--9.5$\mu$m is also suggested.

We review the nucleosynthesis yields of core-collapse supernovae (SNe) for various stellar masses, explosion energies, and metallicities. Comparison with the abundance patterns of metal-poor stars provides excellent opportunities to test the explosion models and their nucleosynthesis. We show that the abundance patterns of extremely metal-poor (EMP) stars, e.g., the excess of C, Co, Zn relative to Fe, are in better agreement with the yields of hyper-energetic explosions (Hypernovae, HNe) rather than normal supernovae. We note that the variation of the abundance patterns of EMP stars are related to the diversity of the Supernova-GRB connection. We summarize the diverse properties of (1) GRB-SNe, (2) Non-GRB HNe/SNe, (3) XRF-SN, and (4) Non-SN GRB. In particular, the Non-SN GRBs (dark hypernovae) have been predicted in order to explain the origin of C-rich EMP stars. We show that these variations and the connection can be modeled in a unified manner with the explosions induced by relativistic jets. Finally, we examine whether the most luminous supernova 2006gy can be consistently explained with the pair-instability supernova model.

The connection between long Gamma Ray Bursts (GRBs) and Supernovae (SNe), have been established through the well observed cases of GRB980425/SN 1998bw, GRB030329/SN 2003dh and GRB031203/SN 2003lw. These events can be explained as the prompt collapse to a black hole (BH) of the core of a massive star (M ~ 40 Msun) that had lost its outer hydrogen and helium envelopes. All these SNe exhibited strong oxygen lines, and their energies were much larger than those of typical SNe, thus these SNe are called Hypernovae (HNe). The case of SN 2006aj/GRB060218 appears different: the GRB was weak and soft (an X-Ray Flash, XRF); the SN is dimmer and has very weak oxygen lines. The explosion energy of SN 2006aj was smaller, as was the ejected mass. In our model, the progenitor star had a smaller mass than other GRB/SNe (M ~ 20 Msun), suggesting that a neutron star (NS) rather than a black hole was formed. If the nascent neutron star was strongly magnetized (a so-called magnetar) and rapidly spinning, it may launch a weak GRB or an XRF. The final fate of 20-30 Msun stars show interesting variety, as seen in the very peculiar Type Ib/c SN 2005bf. This mass range corresponds to the NS to BH transition. We also compare the nucleosynthesis feature of HNe with the metal-poor stars and suggest the Hypernova-First Star connection.

We investigate the evolution of dust that formed at Population III supernova (SN) explosions and its processing through the collisions with the reverse shocks resulting from the interaction of the SN ejecta with the ambient medium. In particular, we investigate the transport of the shocked dust within the SN remnant (SNR), and its effect on the chemical composition, the size distribution, and the total mass of dust surviving in SNR. We find that the evolution of the reverse shock, and hence its effect on the processing of the dust depends on the thickness of the envelope retained by the progenitor star. Furthermore, the transport and survival of the dust grains depend on their initial radius, a_{ini}, and composition: For Type II SNRs expanding into the interstellar medium (ISM) with a density of n_{H,0}=1 cm^{-3}, small grains with a_{ini} < ~ 0.05 micron are completely destroyed by sputtering in the postshock flow, while grains with a_{ini}= 0.05--0.2 micron are trapped into the dense shell behind the forward shock. Very large grains of a_{ini} > ~ 0.2 micron are ejected into the ISM without decreasing their sizes significantly. We find that the total mass fraction of dust that is destroyed by the reverse shock ranges from 0.2 to 1.0, depending on the energy of the explosion and the density of the ambient ISM. The results of our calculations have significant impact on the abundance pattern of subsequent generation of stars that form in the dense shell of primordial SNRs.

Late phase nebular spectra and photometry of Type Ib Supernova (SN) 2005bf taken by the Subaru telescope at ~ 270 and ~ 310 days since the explosion are presented. Emission lines ([OI]6300, 6363, [CaII]7291, 7324, [FeII]7155) show the blueshift of ~ 1,500 - 2,000 km s-1. The [OI] doublet shows a doubly-peaked profile. The line luminosities can be interpreted as coming from a blob or jet containing only ~ 0.1 - 0.4 Msun, in which ~ 0.02 - 0.06 Msun is 56Ni synthesized at the explosion. To explain the blueshift, the blob should either be of unipolar moving at the center-of-mass velocity v ~ 2,000 - 5,000 km s-1, or suffer from self-absorption within the ejecta as seen in SN 1990I. In both interpretations, the low-mass blob component dominates the optical output both at the first peak (~ 20 days) and at the late phase (~ 300 days). The low luminosity at the late phase (the absolute R magnitude M_R ~ -10.2 mag at ~ 270 days) sets the upper limit for the mass of 56Ni < ~ 0.08 Msun, which is in contradiction to the value necessary to explain the second, main peak luminosity (M_R ~ -18.3 mag at ~ 40 days). Encountered by this difficulty in the 56Ni heating model, we suggest an alternative scenario in which the heating source is a newly born, strongly magnetized neutron star (a magnetar) with the surface magnetic field Bmag ~ 10^{14-15} gauss and the initial spin period P0 ~ 10 ms. Then, SN 2005bf could be a link between normal SNe Ib/c and an X-Ray Flash associated SN 2006aj, connected in terms of Bmag and/or P0.

Optical spectroscopy and photometry of SN 2006aj have been performed with the Subaru telescope at t > 200 days after GRB060218, the X-ray Flash with which it was associated. Strong nebular emission-lines with an expansion velocity of v ~ 7,300 km/s were detected. The peaked but relatively broad [OI]6300,6363 suggests the existence of ~ 2 Msun of materials in which ~1.3 Msun is oxygen. The core might be produced by a mildly asymmetric explosion. The spectra are unique among SNe Ic in (1) the absence of [CaII]7291,7324 emission, and (2) a strong emission feature at ~ 7400A, which requires ~ 0.05 Msun of newly-synthesized 58Ni. Such a large amount of stable neutron-rich Ni strongly indicates the formation of a neutron star. The progenitor and the explosion energy are constrained to 18 Msun < Mms < 22 Msun and E ~ (1 - 3) 10^{51} erg, respectively.

The connection between the long Gamma Ray Bursts (GRBs) and Type Ic Supernovae (SNe) has revealed interesting diversity. We review the following types of the GRB-SN connection. (1) GRB-SNe: The three SNe all explode with energies much larger than those of typical SNe, thus being called Hypernovae (HNe). They are massive enough for forming black holes. (2) Non-GRB HNe/SNe: Some HNe are not associated with GRBs. (3) XRF-SN: SN 2006aj associated with X-Ray Flash 060218 is dimmer than GRB-SNe and has very weak oxygen lines. Its progenitor mass is estimated to be small enough to form a neutron star rather than a black hole. (4) Non-SN GRB: Two nearby long GRBs were not associated SNe. Such ``dark HNe'' have been predicted in this talk (i.e., just before the discoveries) in order to explain the origin of C-rich (hyper) metal-poor stars. This would be an important confirmation of the Hypernova-First Star connection. We will show our attempt to explain the diversity in a unified manner with the jet-induced explosion model.

Long-duration gamma-ray bursts (GRBs) are thought to be connected to luminous and energetic supernovae (SNe), called hypernovae (HNe), resulting from the black-hole (BH) forming collapse of massive stars. For recent nearby GRBs~060505 and 060614, however, the expected SNe have not been detected. The upper limits to the SN brightness are about 100 times fainter than GRB-associated HNe (GRB-HNe), corresponding to the upper limits to the ejected $^{56}$Ni masses of $M({\rm ^{56}Ni})\sim 10^{-3}M_\odot$. SNe with a small amount of $^{56}$Ni ejection are observed as faint Type II SNe. HNe and faint SNe are thought to be responsible for the formation of extremely metal-poor (EMP) stars. In this Letter, a relativistic jet-induced BH forming explosion of the 40 $M_\odot$ star is investigated and hydrodynamic and nucleosynthetic models are presented. These models can explain both GRB-HNe and GRBs without bright SNe in a unified manner. Their connection to EMP stars is also discussed. We suggest that GRBs without bright SNe are likely to synthesize $\Mni\sim 10^{-4}$ to $10^{-3}M_\odot$ or $\sim 10^{-6}M_\odot$.

We perform hydrodynamical and nucleosynthesis calculations of core-collapse supernovae (SNe) and hypernovae (HNe) of Population (Pop) III stars. We provide new yields for the main-sequence mass of $M_{\rm MS}=13-50$ $M_{\odot}$ and the explosion energy of $E=1-40\times10^{51}$ ergs to apply for chemical evolution studies. Our HN yields based on the mixing-fallback model of explosions reproduce the observed abundance patterns of extremely metal-poor (EMP) stars ($-4< {\rm [Fe/H]}< -3$), while those of very metal-poor (VMP) stars ($-3<{\rm [Fe/H]} < -2$) are reproduced by the normal SN yields integrated over the Salpeter initial mass function. Moreover, the observed trends of abundance ratios [X/Fe] against [Fe/H] with small dispersions for the EMP stars can be reproduced as a sequence resulting from the various combination of $M_{\rm MS}$ and $E$. This is because we adopt the empirical relation that a larger amount of Fe is ejected by more massive HNe. Our results imply that the observed trends with small dispersions do not necessarily mean the rapid homogeneous mixing in the early galactic halo at [Fe/H] $< -3$, but can be reproduced by the ``inhomogeneous'' chemical evolution model. In addition, we examine how the modifications of the distributions of the electron mole fraction $Y_{\rm e}$ and the density in the presupernova models improve the agreement with observations. In this connection, we discuss possible contributions of nucleosynthesis in the neutrino-driven wind and the accretion disk.

A detailed elemental abundance analysis has been carried out for the very metal-poor ([Fe/H]=-2.7) star BS16934-002, which was identified in our previous work as a star exhibiting large overabundances of Mg and Sc. A comparison of the abundance pattern of this star with that of the well-studied metal-poor star HD122563 indicates excesses of O, Na, Mg, Al, and Sc in BS16934-002. Of particular interest, no excess of C or N is found in this object, in contrast to CS22949-037 and CS29498-043, two previously known carbon-rich, extremely metal-poor stars with excesses of the alpha elements. No established nucleosynthesis model exists that explains the observed abundance pattern of BS16934-002. A supernova model, including mixing and fallback, assuming severe mass loss before explosion, is discussed as a candidate progenitor of BS16934-002.

We present near-infrared spectra of late phase (>200 d) Type Ia supernovae (SNe Ia) taken at the Subaru telescope. The [Fe II] line of SN 2003hv shows a clear flat-topped feature, while that of SN 2005W show less prominent flatness. In addition, a large shift in their line center, varying from -3000 to 1000 (km/s) with respect to the host galaxies, is seen. Such a shift suggests the occurrence of an off-center, non-spherical explosion in the central region, and provides important, new constraints on the explosion models of SNe Ia.

We calculate the evolution of heavy element abundances from C to Zn in the solar neighborhood adopting our new nucleosynthesis yields. Our yields are calculated for wide ranges of metallicity (Z=0-Z_\odot) and the explosion energy (normal supernovae and hypernovae), based on the light curve and spectra fitting of individual supernovae. The elemental abundance ratios are in good agreement with observations. Among the alpha-elements, O, Mg, Si, S, and Ca show a plateau at [Fe/H] < -1, while Ti is underabundant overall. The observed abundance of Zn ([Zn/Fe] ~ 0) can be explained only by the high energy explosion models, which requires a large contribution of hypernovae. The observed decrease in the odd-Z elements (Na, Al, and Cu) toward low [Fe/H] is reproduced by the metallicity effect on nucleosynthesis. The iron-peak elements (Cr, Mn, Co, and Ni) are consistent with the observed mean values at -2.5 < [Fe/H] < -1$, and the observed trend at the lower metallicity can be explained by the energy effect. We also show the abundance ratios and the metallicity distribution functions of the Galactic bulge, halo, and thick disk. Our results suggest that the formation timescale of the thick disk is ~ 1-3 Gyr.

Supernovae connected with long-duration gamma-ray bursts (GRBs) are hyper-energetic explosion resulting from the collapse of very massive stars (about 40Mo, where Mo is the mass of the Sun) stripped of their outher hydrogen and helium envelopes. A very massive progenitor, collapsing to a black hole, was thought to be a requirement for the launch of a GRB. Here we report the results of modelling the spectra and light curve of SN 2006aj, which demonstrate that the supernova had a much smaller explosion energy and ejected much less mass than the other GRB-sueprnovae, suggesting that it was produced by a star whose initial mass was only about 20Mo. A star of this mass is expected to form a neutron star rather than a black hole when its core collapses. The smaller explosion energy of SN 2006aj is matched by the weakness and softness of GRB 060218 (an X-ray flash), and the weakness of the radio flux of the supernova. Our results indicate that the supernova-GRB connection extends to a much broader range of stellar masses than previously thought, possibly involving different physical mechanisms: a `collapsar' for the more massive stars collapsing to a black hole, and magnetic activity of the nascent neutron star for the less massive stars.

We report on detailed spectroscopic studies performed for the secondary star in the black hole binary (micro-quasar) V4641 Sgr in order to examine its surface chemical composition and to see if its surface shows any signature of pollution by ejecta from a supernova explosion. High-resolution spectra of V4641 Sgr observed in the quiescent state in the blue-visual region are compared with those of the two bright well-studied B9 stars (14 Cyg and $\nu$ Cap) observed with the same instrument. The effective temperature of V4641 Sgr (10500 $\pm$ 200 K) is estimated from the strengths of He~{\sc i} lines, while its rotational velocity, $\it v$ sin $\it i$ (95 $\pm$ 10 km s${}^{-1}$), is estimated from the profile of the Mg~{\sc ii} line at 4481 \AA. We obtain abundances of 10 elements and find definite over-abundances of N (by 0.8 dex or more) and Na (by 0.8 dex) in V4641 Sgr. From line-by-line comparisons of eight other elements (C, O, Mg, Al, Si, Ti, Cr, and Fe) between V4641 Sgr and the two reference stars, we conclude that there is no apparent difference in the abundances of these elements between V4641 Sgr and the two normal late B-type stars, which have been reported to have solar abundances. An evolutionary model of a massive close binary system has been constructed to explain the abundances observed in V4641 Sgr. The model suggests that the progenitor of the black hole forming supernova was as massive as ~ 35 Msun on the main-sequence and, after becoming a ~ 10 Msun He star, underwent "dark" explosion which ejected only N and Na-rich outer layer of the He star without radioactive $^{56}$Ni.

We present new nucleosynthesis yields as functions of the stellar mass, metallicity, and explosion energy (corresponding to normal supernovae and Hypernovae). We apply the results to the chemical evolution of the solar neighborhood. Our new yields are based on the new developments in the observational/theoretical studies of supernovae (SNe) and extremely metal-poor (EMP) stars in the halo, which have provided excellent opportunities to test the explosion models and their nucleosynthesis. We use the light curve and spectra fitting of individual SN to estimate the mass of the progenitor, explosion energy, and produced 56Ni mass. Comparison with the abundance patterns of EMP stars has made it possible to determine the model parameters of core-collapse SNe, such as mixing-fallback parameters. More specifically, we take into account the two distinct new classes of massive SNe: 1) very energetic Hypernovae, whose kinetic energy (KE) is more than 10 times the KE of normal core-collapse SNe, and 2) very faint and low energy SNe (Faint SNe). These two new classes of SNe are likely to be ``black-hole-forming'' SNe with rotating or non-rotating black holes. Nucleosynthesis in Hypernovae is characterized by larger abundance ratios (Zn,Co,V,Ti)/Fe and smaller (Mn,Cr)/Fe than normal SNe, which can explain the observed trends of these ratios in EMP stars. Nucleosynthesis in Faint SNe is characterized by a large amount of fall-back, which explains the abundance pattern of the most Fe-poor stars. These comparisons suggest that black-hole-forming SNe made important contributions to the early Galactic (and cosmic) chemical evolution.

The Gamma-Ray Burst 031203 at a redshift z=0.1055 revealed a highly reddened Type Ic Supernova, SN 2003lw, in its afterglow light. This is the third well established case of a link between a long-duration GRB and a type Ic SN. The SN light curve is obtained subtracting the galaxy contribution and is modelled together with two spectra at near-maximum epochs. A red VLT grism 150I spectrum of the SN near peak is used to extend the spectral coverage, and in particular to constrain the uncertain reddening, the most likely value for which is E_{G+H}(B-V) about 1.07 +/- 0.05. Accounting for reddening, SN 2003lw is about 0.3 mag brighter than the prototypical GRB-SN 1998bw. Light curve models yield a 56Ni mass of about 0.55 solar mass. The optimal explosion model is somewhat more massive (ejecta mass about 13 solar mass) and energetic (kinetic energy about 6 times 10^52 erg) than the model for SN 1998bw, implying a massive progenitor (40 - 50 solar mass). The mass at high velocity is not very large (1.4 solar mass above 30000 km/s, but only 0.1 solar mass above 60000 km/s), but is sufficient to cause the observed broad lines. The similarity of SNe 2003lw and 1998bw and the weakness of their related GRBs, GRB031203 and GRB980425, suggest that both GRBs may be normal events viewed slightly off-axis or a weaker but possibly more frequent type of GRB.

The abundance patterns of extremely metal-poor (EMP) stars provide us with important information on nucleosynthesis in supernovae (SNe) formed in a Pop III or EMP environment, and thus on the nature of the first stars in the Universe. We review nucleosynthesis yields of various types of those SNe, focusing on core-collapse (black-hole-forming) SNe with various progenitor masses, explosion energies (including Hypernovae), and asphericity. We discuss the implications of the observed trends in the abundance ratios among iron-peak elements, and the large C/Fe ratio observed in certain EMP stars with particular attention to recently discovered hyper metal-poor (HMP) stars. We show that the abundance pattern of the HMP stars with [Fe/H] < -5 and other EMP stars are in good accord with those of black-hole-forming supernovae, but not pair-instability supernovae. This suggests that black-hole-forming supernovae made important contributions to the early Galactic (and cosmic) chemical evolution. Finally we discuss the nature of First (Pop III) Stars.

Observations and modeling for the light curve (LC) and spectra of supernova (SN) 2005bf are reported. This SN showed unique features: the LC had two maxima, and declined rapidly after the second maximum, while the spectra showed strengthening He lines whose velocity increased with time. The double-peaked LC can be reproduced by a double-peaked $^{56}$Ni distribution, with most $^{56}$Ni at low velocity and a small amount at high velocity. The rapid post-maximum decline requires a large fraction of the $\gamma$-rays to escape from the $^{56}$Ni-dominated region, possibly because of low-density ``holes''. The presence of Balmer lines in the spectrum suggests that the He layer of the progenitor was substantially intact. Increasing $\gamma$-ray deposition in the He layer due to enhanced $\gamma$-ray escape from the $^{56}$Ni-dominated region may explain both the delayed strengthening and the increasing velocity of the He lines. The SN has massive ejecta ($\sim6-7\Msun$), normal kinetic energy ($\sim 1.0-1.5\times 10^{51}$ ergs), high peak bolometric luminosity ($\sim 5\times 10^{42}$ erg s$^{-1}$) for an epoch as late as $\sim$ 40 days, and a large $^{56}$Ni mass ($\sim0.32\Msun$). These properties, and the presence of a small amount of H suggest that the progenitor was initially massive (M$\sim 25-30 \Msun$) and had lost most of its H envelope, and was possibly a WN star. The double-peaked $^{56}$Ni distribution suggests that the explosion may have formed jets that did not reach the He layer. The properties of SN 2005bf resemble those of the explosion of Cassiopeia A.

We investigate a "hypernova" model for gamma-ray bursts (GRBs), i.e., massive C+O star model with relativistic jets. In this model, non-thermal precursors can be produced by the "first" relativistic shell ejected from the star. Main GRBs are produced behind the "first"-shell by the collisions of several relativistic shells. They become visible to distant observers after the colliding region becomes optically thin. We examine six selected conditions using relativistic hydrodynamical simulations and simple analyses. Interestingly, our simulations show that sub-relativistic $(v \sim 0.8c)$ jets from the central engine is sufficient to produce highly-relativistic $(\Gamma > 100)$ shells. We find that the relativistic shells from such a star can reproduce observed GRBs with certain conditions. Two conditions are especially important. One is the sufficiently long duration of the central engine $ \gsim 100$ sec. The other is the existence of a dense-shell somewhere behind the "first"-shell. Under these conditions, both the existence and non-existence of precursors, and long delay between precursors and main GRBs can be explained.

We present BVRI photometry and optical spectroscopy of SN 2005bf near light maximum. The maximum phase is broad and occurred around 2005 May 7, about forty days after the shock breakout. SN 2005bf has a peak bolometric magnitude M_{bol}=-18.0\pm 0.2: while this is not particularly bright, it occurred at an epoch significantly later than other SNe Ibc, indicating that the SN possibly ejected ~0.31 M_{sun} of 56Ni, which is more than the typical amount. The spectra of SN 2005bf around maximum are very similar to those of the Type Ib SNe 1999ex and 1984L about 25-35 days after explosion, displaying prominent He I, Fe II, Ca II H & K and the near-IR triplet P Cygni lines. Except for the strongest lines, He I absorptions are blueshifted by <~6500 km/s, and Fe II by \~7500-8000 km/s. No other SNe Ib have been reported to have their Fe II absorptions blueshifted more than their He I absorptions. Relatively weak H-alpha and very weak H-beta may also exist, blueshifted by ~15,000 km/s. We suggest that SN 2005bf was the explosion of a massive He star, possibly with a trace of a hydrogen envelope.

Recent studies of core-collapse supernovae have revealed the existence of two distinct classes of massive supernovae (SNe): 1) very energetic SNe (Hypernovae), whose kinetic energy (KE) exceeds $10^{52}$ erg, about 10 times the KE of normal core-collapse SNe, and 2) very faint and low energy SNe (E < 0.5 $\times$ $10^{51}$ erg; Faint supernovae). These two classes of supernovae are likely to be "black-hole-forming" supernovae with rotating or non-rotating black holes. We compare their nucleosynthesis yields with the abundances of extremely metal-poor (EMP) stars to identify the Pop III (or first) supernovae. We show that the EMP stars, especially the C-rich type, are likely to be enriched by black-hole-forming supernovae.

The recent discovery of a hyper metal-poor (HMP) star, whose metallicity Fe/H is smaller than 1/100,000 of the solar ratio, together with one earlier HMP star, has raised a challenging question if these HMP stars are the actual first generation, low mass stars in the Universe. We argue that these HMP stars are the second generation stars being formed from gases which were chemically enriched by the first generation supernovae. The key to this solution is the very unusual abundance patterns of these HMP stars with important similarities and differences. We can reproduce these abundance features with the core-collapse ``faint'' supernova models which undergo extensive matter mixing and fallback during the explosion.

Type Ic supernovae, the explosions following the core collapse of massive stars that have previously lost their hydrogen and helium envelopes, are particularly interesting because of the link with long-duration gamma-ray bursts. Although indications exist that these explosions are aspherical, direct evidence has still been missing. Late-time observations of SN 2003jd, a luminous Type Ic supernova, provide such evidence. Recent Subaru and Keck spectra reveal double-peaked profiles in the nebular lines of neutral oxygen and magnesium. These profiles are different from those of known Type Ic supernovae, with or without a gamma-ray burst, and they can be understood if SN 2003jd was an aspherical, axisymmetric explosion viewed from near the equatorial plane. If SN 2003jd was associated with a gamma-ray burst, we missed the burst as it was pointing away from us.

The net optical light curves and spectra of the supernova (SN) 2003dh are obtained from the published spectra of GRB 030329, covering about 6 days before SN maximum to about 60 days after. The bulk of the U-band flux is subtracted from the observed spectra using early-time afterglow templates, because strong line blanketing greatly depresses the UV and U-band SN flux in a metal-rich, fast-moving SN atmosphere. The blue-end spectra of the gamma-ray burst (GRB)connected hypernova SN 1998bw is used to determine the amount of subtraction. The subtraction of a host galaxy template affects the late-time results. The derived SN 2003dh light curves are narrower than those of SN 1998bw, rising as fast before maximum, reaching a possibly fainter maximum, and then declining ~ 1.2-1.4 times faster. We then build UVOIR bolometric SN light curve. Allowing for uncertainties, it can be reproduced with a spherical ejecta model of Mej ~ 7+/-3 Msun, KE ~ (3.5+/-1.5)E52 ergs, with KE/Mej ~ 5 following previous spectrum modelling, and M(Ni56) ~ (0.4 +0.15/-0.1) Msun. This suggests a progenitor main-sequence mass of about 25-40 Msun, lower than SN 1998bw but significantly higher than normal Type Ic SNe and the GRB-unrelated hypernova SN 2002ap.

Late-time near-infrared and optical spectra are presented for the normal-bright SN2003du. At about 300 days after the explosion, the emission profiles of well isolated [FeII] lines (in particular that of the strong 1.644mu feature) trace out the global kinematic distribution of radioactive material in the expanding. In SN2003du, the 1.644 mu [FeII] line shows a flat-topped, profile, indicative of a thick but hollow-centered expanding shell, rather than a strongly-peaked profile that would be expected from a ``center-filled'' distribution.Based on detailed models for exploding Chandrasekhar mass white dwarfs, we show that the feature is consistent with spherical explosion models.Our model predicts central region of non-radioactive electron-capture elements up to 2500--3000km/s as a consequence of burning under high density, and an extended region of 56Ni up to 9,000--10,000km/s. Furthermore our analysis indicates that the 1.644mu [FeII] profile is not consistent with strong mixing between the regions of electron- capture isotopes and the 56Ni layers as is predicted by detailed 3D models for nuclear deflagration fronts. We discuss the possibility that the flat-topped profile could be produced as a result of an infrared catastrophe and conclude that such an explanation is unlikely. We put our results in context to other SNeIa and briefly discuss the implications of our result for the use of SNe Ia as cosmological standard candles.

XRF 030723 is the first X-ray flash (XRF) to show in its optical light curve (LC) a bump that has been interpreted as the signature of a supernova (SN). After subtracting the afterglow component from the observed optical LC of the XRF counterpart, the properties of the putative SN are constrained by means of synthetic LCs of core-collapse SNe. For the redshift range $z \sim 0.3$ -- 1, all possible models require a rather small mass of synthesized \Nifs, i.e. $M$(\Nifs) $\sim$ 0.01 -- 0.3 $\Msun$. The models used to describe the energetic SNe Ic associated with gamma-ray bursts (SNe 1998bw and 2003dh) are too massive for the observed LC. If the relation between ejected \Nifs mass and total ejecta mass established from models of various Type Ic SNe also holds for the putative SN in XRF 030723, the ejecta mass is constrained to be $\sim$ 1 -- 3 $\Msun$ and the kinetic energy $\lsim 1\times 10^{52}$ erg. This corresponds to a progenitor with $15\Msun \lesssim M_{\rm MS}\lesssim 25\Msun$. The SN therefore appears to have properties intermediate between a normal SN Ic like SN 1994I and a more energetic object like SN 2002ap.

The spectra of SN 2003dh, identified in the afterglow of GRB030329, are modeled using radiation transport codes. It is shown that SN 2003dh had a high explosion kinetic energy ($\sim 4 \times 10^{52}$ erg in spherical symmetry), making it one of the most powerful hypernovae observed so far, and supporting the case for association between hypernovae and Gamma Ray Bursts. However, the light curve derived from fitting the spectra suggests that SN 2003dh was not as bright as SN 1998bw, ejecting only $\sim 0.35\Msun$ of \Nifs. The spectra of SN 2003dh resemble those of SN 1998bw around maximum, but later they look more like those of the less energetic hypernova SN 1997ef. The spectra and the inferred light curve can be modeled adopting a density distribution similar to that used for SN 1998bw at $ v > 25,000$\kms but more like that of SN 1997ef at lower velocities. The mass of the ejecta is $\sim 8\Msun$, somewhat less than in the other two hypernovae. The progenitor must have been a massive star ($M \sim 35-40\Msun$), as for other hypernovae. The need to combine different one-dimensional explosion models strongly indicates that SN 2003dh was an asymmetric explosion.

Spectroscopic and spectropolarimetric observations of SN 2003dh/GRB 030329 obtained in 2003 May using the Subaru 8.2 m telescope are presented. The properties of the SN are investigated through a comparison with spectra of the Type Ic hypernovae SNe 1997ef and 1998bw. (Hypernovae being a tentatively defined class of SNe with very broad absorption features: these features suggest a large velocity of the ejected material and possibly a large explosion kinetic energy.) Comparison with spectra of other hypernovae shows that the spectrum of SN 2003dh obtained on 2003 May 8 and 9, i.e., 34-35 rest-frame days after the GRB (for z=0.1685), are similar to those of SN 1997ef obtained ~34-42 days after the fiducial time of explosion of that SN. The match with SN 1998bw spectra is not as good (at rest 7300-8000 A, but again spectra obtained ~33-43 days after GRB 980425 are preferred. This indicates that the SN may have intermediate properties between SNe 1997ef and 1998bw. Based on the analogy with the other hypernovae, the time of explosion of SN 2003dh is then constrained to be between -8 and +2 days of the GRB. The Si and O P-Cygni lines of SN 2003dh seem comparable to those of SN 1997ef, which suggests that the ejected mass in SN 2003dh may match that in SN 1997ef. Polarization was marginally detected at optical wavelengths. This is consistent with measurements of the late afterglow, implying that it mostly originated in the interstellar medium of the host galaxy.

Stars more massive than $\sim$ 20 - 25 \ms form a black hole at the end of their evolution. Stars with non-rotating black holes are likely to collapse "quietly" ejecting a small amount of heavy elements (Faint supernovae). In contrast, stars with rotating black holes are likely to give rise to very energetic supernovae (Hypernovae). We present distinct nucleosynthesis features of these two types of "black-hole-forming" supernovae. Nucleosynthesis in Hypernovae is characterized by larger abundance ratios (Zn,Co,V,Ti)/Fe and smaller (Mn,Cr)/Fe than normal supernovae, which can explain the observed trend of these ratios in extremely metal-poor stars. Nucleosynthesis in Faint supernovae is characterized by a large amount of fall-back. We show that the abundance pattern of the recently discovered most Fe-poor star, HE0107-5240, and other extremely metal-poor carbon-rich stars are in good accord with those of black-hole-forming supernovae, but not pair-instability supernovae. This suggests that black-hole-forming supernovae made important contributions to the early Galactic (and cosmic) chemical evolution. Finally we discuss the nature of First (Pop III) Stars.